2 set duplicate block elements on
6 create Displacement on surface 1 dof 3 fix 0
7 create Displacement on vertex 1 dof 3 dof 2 fix 0
8 create Displacement on vertex 3 dof 3 dof 1 fix 0
9 create Displacement on vertex 4 dof 1 dof 2 dof 3 fix 0
11 create pressure on surface 2 magnitude 2
12 #create force on surface 2 force value 1 direction y
15 block 1 name 'MAT_ELASTIC'
16 block 1 attribute count 2
17 block 1 attribute index 1 1e1 #young modulus
18 block 1 attribute index 2 0.25 #poisson ratio
21 block 2 name "BODY_FORCES"
22 block 2 attribute count 4
24 block 2 attribute index 1 0.01 #material density
26 block 2 attribute index 2 0. #constant acceleration in x-direction
27 block 2 attribute index 3 0. #constant acceleration in y-direction
28 block 2 attribute index 4 0. #constant acceleration in z-direction
30 volume all scheme Tetmesh
31 volume all size auto factor 9
34 # mpirun -np 4 ../../nonlinear_dynamics \
36 # -ts_dt 0.01 -ts_max_time 4 -ts_max_snes_failures -1 \
37 # -my_output_prt -1 -my_max_post_proc_ref_level 0 \
38 # -my_disp_order 2 -my_time_data_file rod_history.in -default_material NEOHOOKEAN 2>&1 | tee log