*HEADING
 :  GENERALIZED PLANE STRAIN ELEMENTS
            RELATIVE MOTION OF BOUNDING PLANES
*NODE
 1,
 5,1.
 41,0., 1.
 45,1., 1.
*NODE,NSET=EXT
 998,.5,.5
*NGEN, NSET=BOT
 1,5
*NGEN,NSET=TOP
 41,45
*NFILL, NSET=ALL1
BOT,TOP,4,10
*NCOPY, SHIFT, OLDSET=ALL1, CHANGE=500, NEWSET=ALL2


*NFILL
 ALL1, ALL2, 5, 100
*NSET, NSET=ORIG, GEN
 23,523,100
*NSET, NSET=CONS, GEN
 25,525,100
***********************************************
** regular generalized plane strain elements **
***********************************************
*ELEMENT, TYPE=CPEG4
 1, 1,2,12,11
*ELGEN,ELSET=EL1
 1, 4,1,1,4,10,4
*********************************************************
** incompatible mode generalized plane strain elements **
*********************************************************
*ELEMENT, TYPE=CPEG4I
 101, 101,102,112,111
*ELGEN,ELSET=EL2
 101, 4,1,1,4,10,4
************************************************************
** reduced intergration generalized plane strain elements **
************************************************************
*ELEMENT, TYPE=CPEG4R
 201, 201,202,212,211
*ELGEN,ELSET=EL3
 201, 4,1,1,4,10,4
**********************************************
** hybrid generalized plane strain elements **
**********************************************
*ELEMENT, TYPE=CPEG4H
 301, 301,302,312,311
*ELGEN,ELSET=EL4
 301, 4,1,1,4,10,4
*****************************************************************
** hybrid, incompatible mode generalized plane strain elements **
*****************************************************************
*ELEMENT, TYPE=CPEG4IH
 401, 401,402,412,411
*ELGEN,ELSET=EL5
 401, 4,1,1,4,10,4
********************************************************************
** hybrid, reduced intergration generalized plane strain elements **
********************************************************************
*ELEMENT, TYPE=CPEG4RH
 501, 501,502,512,511
*ELGEN,ELSET=EL6
 501, 4,1,1,4,10,4
***************************************
** material and property definitions **
***************************************
*SOLID SECTION, ELSET=EL1, MATERIAL=MAT, REFNODE=998
*SOLID SECTION, ELSET=EL2, MATERIAL=MAT, REFNODE=998
*SOLID SECTION, ELSET=EL3, MATERIAL=MAT, REFNODE=998
*SOLID SECTION, ELSET=EL4, MATERIAL=MAT, REFNODE=998
*SOLID SECTION, ELSET=EL5, MATERIAL=MAT, REFNODE=998
*SOLID SECTION, ELSET=EL6, MATERIAL=MAT, REFNODE=998
*MATERIAL, NAME=MAT
*ELASTIC
 30.0E6, 0.3
************************
** history definition **
************************
*STEP, PERTURBATION
 AXIAL TENSION
*STATIC
*NODE PRINT, F=0
*EL PRINT, F=0, POSITION=NODES
*OUTPUT,FIELD
*NODE OUTPUT
 U,RF
*ELEMENT OUTPUT
 S,E
*LOAD CASE,NAME="AXIAL TENSION"
*BOUNDARY,OP=NEW
 998,3,3,0.01
 998,4,5
 ORIG,1,2
 CONS,2
*END LOAD CASE
*LOAD CASE,NAME="RELATIVE ROTATION"
*BOUNDARY,OP=NEW
 998,3,4
 998,5,5,0.01
 ORIG,1,2
 CONS,2
*END LOAD CASE
*END STEP