*HEADING
: MEMBRANE PATCH TEST  [CPEG4RH]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*restart,write
*NODE
1000,0,0
10,0,0
20,.24,0.
30,.24,.12
40,0.,.12
50,.12,0.
60,.24,.06
70,.12,.12
80,.0,.06
102,.04,.02
104,.18,.03
204,.16,.08
202,.08,.08
*NGEN
102,104
202,204
102,202,50
104,204,50
10,102,46
40,202,81
30,204,87
20,104,42
10,40,15
*ELEMENT,TYPE=CPEG4RH,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,
*SOLID SECTION,MATERIAL=A1,ELSET=EALL,REFNODE=1000
.001, 
*MATERIAL,NAME=A1
*ELASTIC,TYPE=ISOTROPIC
1.E6,.25
*NSET,NSET=SIDES,GENERATE
10,80,10
*SURFACE, NAME=DIAG
2,S3
2,S2
1,S3
1,S2
*SURFACE, NAME=USHAPE
2,S1
3,S1
3,S2
3,S3
2,S3
*SURFACE, NAME=CLOSINT
3,S1
3,S2
3,S3
3,S4
*SURFACE, NAME=CLOSEXT
1,S3
4,S4
5,S1
2,S2
*STEP,PERTURBATION
*STATIC
*BOUNDARY
1000,3,5
10,1,2
20,1,,.24E-3
20,2,,.12E-3
30,1,,.3E-3
30,2,,.24E-3
40,1,,.06E-3
40,2,,.12E-3
*ENERGY FILE
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*ENERGY PRINT
*EL PRINT
S, 
E, 
*EL PRINT,POSITION=NODES
S, 
E, 
*EL PRINT,POSITION=AVERAGED AT NODES
S, 
E, 
*EL PRINT,POSITION=CENTROID
S, 
E, 
*EL FILE
S, 
E, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S, 
E, 
*EL FILE,POSITION=NODES
S, 
E, 
*EL FILE,POSITION=AVERAGED AT NODES
S, 
E, 
*EL FILE,POSITION=CENTROID
S, 
E, 
**
*SECTION PRINT, NAME=DIAG,   SURFACE=DIAG,   AXES=LOCAL
SOAREA,SOF
*SECTION FILE, NAME=DIAG,    SURFACE=DIAG,   AXES=LOCAL
SOAREA,SOF
*SECTION PRINT, NAME=USHAPE, SURFACE=USHAPE, AXES=LOCAL
SOAREA,SOF
*SECTION FILE, NAME=USHAPE,  SURFACE=USHAPE, AXES=LOCAL
SOAREA,SOF
*SECTION PRINT, NAME=closint,SURFACE=closint,AXES=LOCAL
SOAREA,SOF
*SECTION FILE, NAME=closint, SURFACE=closint,AXES=LOCAL
SOAREA,SOF
*SECTION PRINT, NAME=closext,SURFACE=closext,AXES=LOCAL
SOAREA,SOF
*SECTION FILE, NAME=closext, SURFACE=closext,AXES=LOCAL
SOAREA,SOF
**
*END STEP
*STEP,NLGEOM
*STATIC
*CONTROLS,PARAMETERS=FIELD
1.E-14, 1., 10000. , , 1.E-12
*DLOAD
1, P1, -10000.
2, P4, -10000.
4, P2, -10000.
5, P3, -10000.
*BOUNDARY
10,1,2
20,2
1000,3,5
*ENERGY PRINT
*ENERGY FILE
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*EL PRINT,SUMMARY=NO,ELSET=EL3
S,E
*EL FILE,ELSET=EL3
S,E
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=EL3
S,E
*NODE PRINT,SUMMARY=NO
U, 
RF, 
*NODE FILE
U,RF
*OUTPUT,FIELD
*NODE OUTPUT
U,RF
*END STEP
*STEP,PERTURBATION
*STATIC
*BOUNDARY
1000,3,5
10,1,2
20,1,,.24E-3
20,2,,.12E-3
30,1,,.3E-3
30,2,,.24E-3
40,1,,.06E-3
40,2,,.12E-3
*ENERGY PRINT
COORD
*NODE FILE
U,RF,CF,COORD
*OUTPUT,FIELD
*NODE OUTPUT
U,RF,CF,COORD
*ENERGY FILE
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*ENERGY PRINT
*EL PRINT
S, 
E, 
*EL PRINT,POSITION=NODES
S, 
E, 
*EL PRINT,POSITION=AVERAGED AT NODES
S, 
E, 
*EL PRINT,POSITION=CENTROID
S, 
E, 
*EL FILE
S, 
E, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S, 
E, 
*EL FILE,POSITION=NODES
S, 
E, 
*EL FILE,POSITION=AVERAGED AT NODES
S, 
E, 
*EL FILE,POSITION=CENTROID
S, 
E, 
*END STEP