*HEADING
: Membrane patch test [S4RS]
driven submodel: substs_s_shel_c3d8r_p_xpl.inp [C3d8R]
displacement field applied:
U_x = A (x + y/2), U_y = A (y + x/2), U_z = 0
**
*PREPRINT,ECHO=NO
*RESTART,WRITE
*NODE,NSET=ALL
10,0,0
20,.24,0.
30,.24,.12
40,0.,.12
102,.04,.02
104,.18,.03
204,.16,.08
202,.08,.08
*NSET,NSET=SIDES,GENERATE
10,40,10
*ELEMENT,TYPE=S4RS, ELSET=EALL
1,10,20,104,102
2,10,102,202,40
3,102,104,204,202
4,104,20,30,204
5,202,204,30,40
*ELSET,ELSET=EL3
3,
*SHELL SECTION,MATERIAL=A1, ELSET=EALL
0.0125,5
*MATERIAL,NAME=A1
*ELASTIC,TYPE=ISOTROPIC
1.E6,.25
*DENSITY
1.0, 
*NSET, NSET=QA_TEST
ALL,
**
*AMPLITUDE, NAME=QSTATIC, SMOOTH=0.25, TIME=TOTAL TIME
0.0,0.0,0.01,1.0,0.02,2.0,0.03,1.0
**
** Step 1: General Quasi-static Step
**
*STEP, NLGEOM=NO
*DYNAMIC, EXPLICIT
,0.01
*BOUNDARY, AMPLITUDE=QSTATIC
10,1,2
ALL,3,5
20,1,,.24E-3
20,2,,.12E-3
30,1,,.3E-3
30,2,,.24E-3
40,1,,.06E-3
40,2,,.12E-3
**
** output requests for submodeling
**
*FILE OUTPUT, NUMBER INTERVAL=20
*NODE FILE
U,
**
** output requests
**
*OUTPUT, FIELD, NUMBER INTERVAL=10
*NODE OUTPUT
U,RF
*ELEMENT OUTPUT
S,LE
*OUTPUT, HISTORY, TIME INTERVAL=0.01
*ENERGY OUTPUT, VARIABLE=PRESELECT
**
**output for QA
**
*OUTPUT, FIELD, NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TEST
U, 
*END STEP
**
** Step 2: General Quasi-static Step
**         displacement BC's of STEP 1 doubled
**
*STEP, NLGEOM=NO
*DYNAMIC, EXPLICIT
,0.01
*BOUNDARY, AMPLITUDE=QSTATIC
10,1,2
ALL,3,5
20,1,,.24E-3
20,2,,.12E-3
30,1,,.3E-3
30,2,,.24E-3
40,1,,.06E-3
40,2,,.12E-3
*END STEP
**
** Step 3: General Quasi-static Step
**         Displacement BC's same as STEP 1
**
*STEP, NLGEOM=NO
*DYNAMIC, EXPLICIT
,0.01
*BOUNDARY, AMPLITUDE=QSTATIC
10,1,2
ALL,3,5
20,1,,.24E-3
20,2,,.12E-3
30,1,,.3E-3
30,2,,.24E-3
40,1,,.06E-3
40,2,,.12E-3
*END STEP