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


*NSET, NSET=ORIG
 23,123
*NSET, NSET=CONS
 25,125
***********************************************
** regular generalized plane strain elements **
***********************************************
*ELEMENT, TYPE=CPEG3T
 1, 1,2,12
 2, 1,12,11
*ELGEN,ELSET=EL1
 1, 4,1,2,4,10,8
 2, 4,1,2,4,10,8
**********************************************
** hybrid generalized plane strain elements **
**********************************************
*ELEMENT, TYPE=CPEG3HT
 101, 101,102,112
 102, 101,112,111
*ELGEN,ELSET=EL2
 101, 4,1,2,4,10,8
 102, 4,1,2,4,10,8
***************************************
** material and property definitions **
***************************************
*SOLID SECTION, ELSET=EL1, MATERIAL=MAT, REFNODE=998
*SOLID SECTION, ELSET=EL2, 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