UniaxialTensionPeriodic

2d-periodicy-amr0

	Two-dimensional
	Serial
	Validated using check script
	./bin/mechanics-2d-g++ tests/UniaxialTensionPeriodic/input amr.max_level="0"

2d-periodicy-amr1

	Two-dimensional
	Serial
	Validated using check script
	./bin/mechanics-2d-g++ tests/UniaxialTensionPeriodic/input amr.max_level="1"

2d-periodicy-amr2

	Two-dimensional
	Serial
	Validated using check script
	./bin/mechanics-2d-g++ tests/UniaxialTensionPeriodic/input amr.max_level="2"

Input file (../../tests/UniaxialTensionPeriodic/input)

#@  [2d-periodicy-amr0]
#@  exe=mechanics
#@  dim=2
#@  args=amr.max_level=0
#@
#@  [2d-periodicy-amr1]
#@  exe=mechanics
#@  dim=2
#@  args=amr.max_level=1
#@
#@  [2d-periodicy-amr2]
#@  exe=mechanics
#@  dim=2
#@  args=amr.max_level=2
#@


alamo.program                 = mechanics
alamo.program.mechanics.model = linear.isotropic
plot_file		      = tests/UniaxialTensionPeriodic/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              = 0 0 0
explicitmesh.hi1              = 47 47 47

explicitmesh.lo2              = 32 32 32
explicitmesh.hi2              = 81 81 81

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

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

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

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

### --- Faces
bc.constant.type.xlo = disp disp disp
bc.constant.type.xhi = disp disp disp
bc.constant.val.xhi = 0.0 0.1 0.0
bc.constant.val.xhiylo = 0.0 0.1 0.0
bc.constant.val.xhiyhi = 0.0 0.1 0.0
bc.constant.type.ylo=periodic periodic periodic
bc.constant.type.yhi=periodic periodic periodic
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.0 0.1 0.0
bc.constant.val.xhiyhi = 0.0 0.1 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


print_model = 1
print_residual = 1