Solid

linear-isotropic

Three-dimensional

Serial

Not validated

./bin/mechanics-3d-g++ tests/Solid/input

linear-cubic

Three-dimensional

Serial

Validated using check script

./bin/mechanics-3d-g++ tests/Solid/input  alamo.program.mechanics.model="linear.cubic" model1.C11="168.3" model1.C12="1.221" model1.C44="0.757" ignore="model1.E model1.nu"

affine-cubic

Three-dimensional

Serial

Validated using check script

./bin/mechanics-3d-g++ tests/Solid/input  alamo.program.mechanics.model="affine.cubic" model1.C11="168.3" model1.C12="1.221" model1.C44="0.757" model1.F0="0.001 0 0 0 0.001 0 0 0 0.001" ignore="model1.E model1.nu"

affine-hexagonal

Three-dimensional

Serial

Validated using check script

./bin/mechanics-3d-g++ tests/Solid/input  alamo.program.mechanics.model="affine.hexagonal" model1.C11="0.597" model1.C12="0.262" model1.C13="0.217" model1.C33="0.617" model1.C44="0.164" model1.F0="0.001 0 0 0 0.001 0 0 0 0.001" ignore="model1.E model1.nu"

j2

Three-dimensional

Serial

Validated using check script

10.42s (beaker) 7.02s (statler) 11.09s (github)

./bin/mechanics-3d-g++ tests/Solid/input  timestep="0.001" alamo.program.mechanics.model="affine.j2" model1.E="210" model1.nu="0.3" model1.sigma0="0.2" bc.tension_test.disp="(0,0.25,0.75,1.0:0,0.002,-0.002,0)"

affine-isotropic

Three-dimensional

Serial

Not validated

./bin/mechanics-3d-g++ tests/Solid/input  alamo.program.mechanics.model="affine.isotropic" model1.F0="0.001 0 0 0 0.001 0 0 0 0.001" model1.E="210" model1.nu="0.3"

neo-hookean

Three-dimensional

Serial

Validated using check script

10.73s (beaker) 8.17s (statler) 11.35s (github)

./bin/mechanics-3d-g++ tests/Solid/input  timestep="0.01" alamo.program.mechanics.model="finite.neohookean" model1.mu="3.0" model1.kappa="6.5" bc.tension_test.disp="(0,1:0,1)" solver.nriters="10" ignore="model1.E model1.nu"

neo-hookean-2d

Two-dimensional

Serial

Validated using check script

./bin/mechanics-2d-g++ tests/Solid/input  timestep="0.01" alamo.program.mechanics.model="finite.neohookean" model1.mu="3.0" model1.kappa="6.5" bc.tension_test.disp="(0,1:0,1)" solver.nriters="10" ignore="model1.E model1.nu"

pseudolinear-cubic

Three-dimensional

Serial

Validated using check script

./bin/mechanics-3d-g++ tests/Solid/input  alamo.program.mechanics.model="finite.pseudolinear.cubic" model1.C11="168.3" model1.C12="1.221" model1.C44="0.757" ignore="model1.E model1.nu"
Input file (../../tests/Solid/input)
#@
#@  [linear-isotropic]
#@  exe=mechanics
#@  dim = 3
#@  check = false
#@
#@  [linear-cubic]
#@  exe=mechanics
#@  dim = 3
#@  check-file = reference/linear-cubic.dat
#@  args = alamo.program.mechanics.model=linear.cubic
#@  args = model1.C11 = 168.3
#@  args = model1.C12 = 1.221
#@  args = model1.C44 = 0.757
#@  ignore = model1.E model1.nu
#@  coverage = true
#@  
#@  [affine-cubic]
#@  exe=mechanics
#@  dim = 3
#@  check-file = reference/affine-cubic.dat
#@  args = alamo.program.mechanics.model=affine.cubic
#@  args = model1.C11 = 168.3
#@  args = model1.C12 = 1.221
#@  args = model1.C44 = 0.757
#@  args = model1.F0=0.001 0 0 0 0.001 0 0 0 0.001
#@  ignore = model1.E model1.nu
#@  coverage = true
#@
#@
#@  [affine-hexagonal]
#@  exe=mechanics
#@  dim = 3
#@  check-file = reference/affine-hexagonal.dat
#@  args = alamo.program.mechanics.model=affine.hexagonal
#@  args = model1.C11 = 0.597
#@  args = model1.C12 = 0.262
#@  args = model1.C13 = 0.217
#@  args = model1.C33 = 0.617
#@  args = model1.C44 = 0.164
#@  args = model1.F0=0.001 0 0 0 0.001 0 0 0 0.001
#@  ignore = model1.E model1.nu
#@  coverage = true
#@
#@  [j2]
#@  exe=mechanics
#@  dim = 3
#@  check-file = reference/j2.dat
#@  args = timestep = 0.001
#@  args = alamo.program.mechanics.model=affine.j2
#@  args = model1.E=210
#@  args = model1.nu=0.3
#@  args = model1.sigma0=0.2
#@  args = bc.tension_test.disp = (0,0.25,0.75,1.0:0,0.002,-0.002,0)
#@  benchmark-beaker = 10.42
#@  benchmark-statler = 7.02
#@  benchmark-github = 11.09
#@
#@  [affine-isotropic]
#@  exe=mechanics
#@  dim = 3
#@  check = false
#@  args = alamo.program.mechanics.model=affine.isotropic
#@  args = model1.F0=0.001 0 0 0 0.001 0 0 0 0.001
#@  args = model1.E=210 
#@  args = model1.nu=0.3
#@  coverage=true
#@
#@  [neo-hookean]
#@  exe=mechanics
#@  dim=3
#@  check-file = reference/neo-hookean.dat
#@  args = timestep=0.01
#@  args = alamo.program.mechanics.model=finite.neohookean
#@  args = model1.mu=3.0
#@  args = model1.kappa=6.5
#@  args = bc.tension_test.disp=(0,1:0,1)
#@  args = solver.nriters=10
#@  ignore = model1.E model1.nu
#@  benchmark-beaker = 10.73
#@  benchmark-statler = 8.17
#@  benchmark-github = 11.35
#@ 
#@  [neo-hookean-2d]
#@  exe=mechanics
#@  dim=2
#@  check-file = reference/neo-hookean-2d.dat
#@  args = timestep=0.01
#@  args = alamo.program.mechanics.model=finite.neohookean
#@  args = model1.mu=3.0
#@  args = model1.kappa=6.5
#@  args = bc.tension_test.disp=(0,1:0,1)
#@  args = solver.nriters=10
#@  ignore = model1.E model1.nu
#@  coverage = true
#@
#@  [pseudolinear-cubic]
#@  exe=mechanics
#@  dim = 3
#@  check-file = reference/pseudolinear-cubic.dat
#@  args = alamo.program.mechanics.model=finite.pseudolinear.cubic
#@  args = model1.C11=168.3
#@  args = model1.C12=1.221
#@  args = model1.C44=0.757
#@  ignore = model1.E model1.nu
#@  coverage=true
#@

alamo.program = mechanics


plot_file		    = tests/Solid/output

type=static

timestep		    = 0.01
stop_time		    = 1.0

# amr parameters
amr.plot_dt		    = 0.1
amr.max_level		    = 0
amr.n_cell		    = 4 4 4
amr.blocking_factor         = 2

amr.thermo.int = 1
amr.thermo.plot_int = 1

# geometry
geometry.prob_lo	    = 0 0 0
geometry.prob_hi	    = 1 1 1

nmodels = 1
model1.E=210 
model1.nu=0.3

solver.verbose = 3
solver.nriters = 1
solver.max_iter = 30 #30

bc.type = tension_test
bc.tension_test.type = uniaxial_stress
bc.tension_test.disp=(0,1:0,0.01)