*heading
*restart,write
*NODE,system=c
 1,10,0
 101,20,0
 111,20,90
 11,10,90
**
 10000,0,0
 10001,0,0
*ngen,line=c,nset=ID,system=rc
 1,11,5,,0,0
*ngen,line=c,nset=OD,system=rc
 101,111,5,,0,0
*nfill
 ID,OD,4,25
*node,nset=xtra
 99998,0.,0.
*element,type=cpeg4t
 1,1,26,31,6
*elgen,elset=eall
 1,4,25,1,2,5,4
*SOLID SECTION,ELSET=eall,MATERIAL=MAT,refnode=xtra
 1., 
*MATERIAL,NAME=MAT
*ELASTIC
 10.,0.3
*density
 1.e2, 
*expansion
 1.e-3, 
*CONDUCTIVITY
 1.e-2, 
*BOUNDARY
 ID,1,2
 ID,11,11
*NSET,NSET=down,generate
 1,101,25
*NSET,NSET=up,generate
 11,111,25
*NSET,NSET=down1
 101, 
*NSET,NSET=up1
 111, 
*MPC
 CYCLSYM,down,up
**CYCLSYM,1,11
**CYCLSYM,26,36
**CYCLSYM,51,61
**CYCLSYM,76,86
**CYCLSYM,down1,up1
**
**STEP,PERT
*STEP,inc=10000
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
 0.01,1
*cload
 106,1,1.e2
*BOUNDARY
 101,11,11,100
 111,11,11,100
 106,11,11,200
 99998,3,5
*NODE PRINT
 U,
 NT,
*NODE FILE
 U,NT
*OUTPUT,FIELD
*NODE OUTPUT
 U,NT
*END STEP