UniaxialTension

2D-serial-1level

Two-dimensional

Serial

Validated using check script

./bin/mechanics-2d-g++ tests/UniaxialTension/input  amr.max_level="0"

2D-serial-2levels

Two-dimensional

Serial

Validated using check script

./bin/mechanics-2d-g++ tests/UniaxialTension/input  amr.max_level="1"

2D-serial-3levels

Two-dimensional

Serial

Validated using check script

./bin/mechanics-2d-g++ tests/UniaxialTension/input

2D-parallel-3levels

Two-dimensional

Parallel (4 procs)

Validated using check script

mpiexec -np 4 ./bin/mechanics-2d-g++ tests/UniaxialTension/input

3D-serial-3levels

Three-dimensional

Serial

Validated using check script

8.29s (github)

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

3D-parallel-3levels

Three-dimensional

Parallel (4 procs)

Validated using check script

mpiexec -np 4 ./bin/mechanics-3d-g++ tests/UniaxialTension/input
Input file (../../tests/UniaxialTension/input)
#@  [2D-serial-1level]
#@  exe=mechanics
#@  dim=2  
#@  nprocs=1
#@  args=amr.max_level=0
#@    
#@  [2D-serial-2levels]
#@  exe=mechanics
#@  dim=2  
#@  nprocs=1
#@  args=amr.max_level=1
#@    
#@  [2D-serial-3levels]
#@  exe=mechanics
#@  dim=2  
#@  nprocs=1
#@    
#@  [2D-parallel-3levels]
#@  exe=mechanics
#@  dim=2  
#@  nprocs=4
#@    
#@  [3D-serial-3levels]
#@  exe=mechanics
#@  dim=3
#@  nprocs=1
#@  benchmark-github = 8.29
#@
#@  [3D-parallel-3levels]
#@  exe=mechanics
#@  dim=3
#@  nprocs=4
#@  

alamo.program                 = mechanics
alamo.program.mechanics.model = linear.isotropic
plot_file		      = tests/UniaxialTension/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		      = 2
amr.n_cell		      = 32 32 32
amr.blocking_factor	      = 4
amr.regrid_int		      = 1
amr.grid_eff		      = 1.0
amr.cell.all                  = 1

# use an explicit mesh (i.e. no adaptive meshing)
explicitmesh.on               = 1
explicitmesh.lo1              = 16 16 16
explicitmesh.hi1              = 47 47 47
explicitmesh.lo2              = 48 48 48
explicitmesh.hi2              = 79 79 79

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

# elastic moduli
model1.planestress = true
model1.E = 210 
model1.nu = 0.3

solver.verbose = 3
solver.nriters = 1
solver.max_iter = 50

### UNIAXIAL TENSION ###
bc.type = constant

### --- Faces
bc.constant.type.xlo = disp trac trac
bc.constant.type.xhi = disp trac trac
bc.constant.val.xhi = 0.1 0.0 0.0
bc.constant.type.ylo = trac disp trac
bc.constant.type.yhi = trac trac trac
bc.constant.type.zlo = trac trac disp
bc.constant.type.zhi = trac trac trac

### --- Edges
###     XY edges
bc.constant.type.xloylo = disp trac trac
bc.constant.type.xloyhi = disp trac trac
bc.constant.type.xhiylo = disp trac trac
bc.constant.type.xhiyhi = disp trac trac
bc.constant.val.xhiylo = 0.1 0.0 0.0
bc.constant.val.xhiyhi = 0.1 0.0 0.0
###     ZX edges
bc.constant.type.zloxlo = disp trac trac
bc.constant.type.zhixlo = disp trac trac
bc.constant.type.zloxhi = disp trac trac
bc.constant.type.zhixhi = disp trac trac
bc.constant.val.zloxhi = 0.1 0.0 0.0
bc.constant.val.zhixhi = 0.1 0.0 0.0
###     YZ edges
bc.constant.type.ylozlo = trac disp disp
bc.constant.type.ylozhi = trac disp trac
bc.constant.type.yhizlo = trac trac disp
bc.constant.type.yhizhi = trac trac trac


### --- Corners
###     on xmin face
bc.constant.type.xloylozlo = disp disp disp
bc.constant.type.xloylozhi = disp disp trac
bc.constant.type.xloyhizlo = disp trac disp
bc.constant.type.xloyhizhi = disp trac trac
###     on xmax face
bc.constant.type.xhiylozlo = disp disp disp
bc.constant.val.xhiylozlo = 0.1 0.0 0.0
bc.constant.type.xhiylozhi = disp disp trac
bc.constant.val.xhiylozhi = 0.1 0.0 0.0
bc.constant.type.xhiyhizlo = disp trac disp
bc.constant.val.xhiyhizlo = 0.1 0.0 0.0
bc.constant.type.xhiyhizhi = disp trac trac
bc.constant.val.xhiyhizhi = 0.1 0.0 0.0