*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
***********************************************
** regular generalized plane strain elements **
***********************************************
*ELEMENT, TYPE=CPEG4T
 1, 1,2,12,11
*ELGEN,ELSET=EL1
 1, 4,1,1,4,10,4
***************************************************
** hybrid mode generalized plane strain elements **
***************************************************
*ELEMENT, TYPE=CPEG4HT
 101, 101,102,112,111
*ELGEN,ELSET=EL2
 101, 4,1,1,4,10,4
************************************************************
** reduced intergration generalized plane strain elements **
************************************************************
*ELEMENT, TYPE=CPEG4RT
 201, 201,202,212,211
*ELGEN,ELSET=EL3
 201, 4,1,1,4,10,4
******************************************************
** reduced hybrid generalized plane strain elements **
******************************************************
*ELEMENT, TYPE=CPEG4RHT
 301, 301,302,312,311
*ELGEN,ELSET=EL4
 301, 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
*MATERIAL, NAME=MAT
*ELASTIC
 30.0E6, 0.3
*CONDUCTIVITY
 1.0
************************
** 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