HEADING
SUBMODEL FROM membsub_g_m3d4r_xpl, 
which are based on the patch test. The analytical sol:
u1=C01+C11x+C21y;u2=C02+C12x+C22y;
C01=C02=0,C11=C22=1e-3,C12=C21=0.5E-3; 
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE
1,.12,0.
9,.12,.12
801,.24,0.
809,.24,.12
*NGEN,NSET=LEFT
1,9,2
*NGEN,NSET=RIGHT
801,809,2
*NFILL,NSET=ALL
LEFT,RIGHT,4,200
*ELEMENT,TYPE=M3D4R, ELSET=EALL
1,1,201,203,3
*ELGEN,ELSET=EALL
1,4,200,1,4,2,4
*NSET, NSET=QA_TEST
ALL
*ELSET, ELSET=QA_TEST
EALL,
*MEMBRANE SECTION,MATERIAL=A1, ELSET=EALL
.001, 
*MATERIAL,NAME=A1
*ELASTIC,TYPE=ISOTROPIC
1.E6,.25
*DENSITY
100.0, 
*NSET,NSET=SIDES,GENERATE
10,80,10
*AMPLITUDE, NAME=RAMP
0,0,1., 1.0
**
**
*NSET,NSET=DRIV,GEN
1,9,2
*SUBMODEL,EXTER=0.0025
DRIV, 
*STEP
*DYNAMIC, EXPLICIT
,1.0
*BOUNDARY,TYPE=VELOCITY
801,1,,.24E-3
801,2,,.12E-3
809,1,,.3E-3
809,2,,.24E-3
*BOUNDARY,SUBMODEL,STEP=1
DRIV,1,3
*FILE OUTPUT, NUM=10
*NODE FILE
 U, 
********************************
*OUTPUT,HISTORY,TIME INTERVAL=1.0
*ENERGY OUTPUT,VAR=PRESELECT
*OUTPUT,FIELD,NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TEST
U,
*ELEMENT OUTPUT, ELSET=QA_TEST
LE,
***********************
*OUTPUT,FIELD,NUMBER INTERVAL=10
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
*OUTPUT,HISTORY,TIME INTERVAL=0.1
*NODE OUTPUT, NSET=ALL
U,
*****************************
*END STEP
**
**
*STEP
*DYNAMIC, EXPLICIT
 ,1.0
*BOUNDARY, TYPE=VELOCITY
801,1,,-.24E-3
801,2,,-.12E-3
809,1,,-.3E-3
809,2,,-.24E-3
*BOUNDARY,SUBMODEL,STEP=2
DRIV,1,3
*ENDSTEP
**
**
*STEP
*DYNAMIC, EXPLICIT
,1.0
*CLOAD,AMP=RAMP
801, 1, 0.6
809, 1, 0.6
809, 2, 1.2
*BOUNDARY,TYPE=VELOCITY
801,2,2,0.0
*BOUNDARY,SUBMODEL,STEP=3
DRIV,1,3
*END STEP
**********************