EshelbyFiniteKinematics

2D-serial

Two-dimensional

Serial

Validated using check script

./bin/mechanics-2d-g++ tests/EshelbyFiniteKinematics/input
Input file (../../tests/EshelbyFiniteKinematics/input)

#@  [2D-serial]
#@  exe = mechanics
#@  dim    = 2
#@  nprocs = 1

alamo.program = mechanics
alamo.program.mechanics.model = finite.neohookeanpre

plot_file               = tests/EshelbyFiniteKinematics/output

# this is not a time integration, so do
# exactly one timestep and then quit
timestep             = 0.1
stop_time                   = 0.1

# amr parameters
amr.plot_int             = 1
amr.max_level                    = 5
amr.n_cell               = 32 32 32
amr.blocking_factor      = 8
amr.regrid_int                 = 1
amr.grid_eff             = 1.0
amr.cell.all                = 1

# geometry
geometry.prob_lo          = -8 -8 -8 
geometry.prob_hi          = 8 8 8
geometry.is_periodic     = 0 0 0

# ellipse configuration
ic.type        = ellipse
ic.ellipse.a   = 1.0 0.75 0.5  # ellipse radii
ic.ellipse.x0  = 0 0 0 # location of ellipse center
ic.ellipse.eps = 0.1 # diffuse boundary

# elastic moduli
nmodels = 2
model1.E = 210 
model1.nu = 0.3
model1.eps0  = 0.05 0 0 0 0.05 0 0 0 0.05 # eigenstrain
model2.E = 210 
model2.nu = 0.3
model2.eps0  = 0 0 0 0 0 0 0 0 0 # eigenstrain

solver.verbose = 3
solver.nriters = 1000
solver.nrtolerance = 1E-6
solver.max_iter = 20

bc.type = constant