*HEADING
 :  COUPLED 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=300, NEWSET=ALL2


*NFILL
 ALL1, ALL2, 3, 100
*NSET, NSET=ORIG, GEN
 23,323,100
*NSET, NSET=CONS, GEN
 25,325,100
***********************************************
** 8-noded generalized plane strain elements **
***********************************************
*ELEMENT, TYPE=CPEG8T
 1, 1,3,23,21, 2,13,22,11
*ELGEN,ELSET=EL1
 1, 2,2,1,2,20,2
***********************************************
** 8-noded generalized plane strain elements **
***********************************************
*ELEMENT, TYPE=CPEG8T
 101, 101,103,123,121, 102,113,122,111
*ELGEN,ELSET=EL2
 101, 2,2,1,2,20,2
***********************************************
** 8-noded hybrid gen. plane strain elements **
***********************************************
*ELEMENT, TYPE=CPEG8HT
 201, 201,203,223,221, 202,213,222,211
*ELGEN,ELSET=EL3
 201, 2,2,1,2,20,2
*********************************************************
** 8-noded red. int. hybrid gen. plane strain elements **
*********************************************************
*ELEMENT, TYPE=CPEG8RHT
 301, 301,303,323,321, 302,313,322,311
*ELGEN,ELSET=EL4 
 301, 2,2,1,2,20,2
***************************************
** 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
*CONDUCTIVITY
 1.0
************************
** history definition **
************************
*STEP,PERTURBATION
*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