*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=100, NEWSET=ALL2


*NCOPY, SHIFT, OLDSET=ALL1, CHANGE=200, NEWSET=ALL2


*NCOPY, SHIFT, OLDSET=ALL1, CHANGE=300, NEWSET=ALL2


*NSET, NSET=ORIG
 23,123,223,323
*NSET, NSET=CONS
 25,125,225,325
*********************************************************
** 6 noded - regular generalized plane strain elements **
*********************************************************
*ELEMENT, TYPE=CPEG6
 1, 1,3,23,2,13,12
 2, 1,23,21,12,22,11
*ELGEN,ELSET=EL1
 1, 2,2,2,2,20,4
 2, 2,2,2,2,20,4
********************************************************
** 6 noded - hybrid generalized plane strain elements **
********************************************************
*ELEMENT, TYPE=CPEG6H
 101, 101,103,123, 102,113,112
 102, 101,123,121, 112,122,111
*ELGEN,ELSET=EL2
 101, 2,2,2,2,20,4
 102, 2,2,2,2,20,4
*******************************************************************
** 6 noded - regular modified, generalized plane strain elements **
*******************************************************************
*ELEMENT, TYPE=CPEG6M
 201, 201,203,223,202,213,212
 202, 201,223,221,212,222,211
*ELGEN,ELSET=EL3
 201, 2,2,2,2,20,4
 202, 2,2,2,2,20,4
******************************************************************
** 6 noded - hybrid modified, generalized plane strain elements **
******************************************************************
*ELEMENT, TYPE=CPEG6MH
 301, 301,303,323,302,313,312
 302, 301,323,321,312,322,311
*ELGEN,ELSET=EL4
 301, 2,2,2,2,20,4
 302, 2,2,2,2,20,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
*MATERIAL, NAME=MAT
*ELASTIC
 30.0E6, 0.3
************************
** history definition **
************************
*STEP, PERTURBATION
*STATIC
*NODE PRINT, F=0
*EL PRINT,F=0
*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,0.0
 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