*HEADING
ELEMENTS        : B31,  PSI34
ANALYSIS        : SMALL DISPLACEMENT
MATERIAL MODEL  : NONLINEAR ISOTROPIC, TEMPERATURE DEPENDENT
*PRE PRINT,ECHO=NO,MODEL=YES,CONTACT=NO,HISTORY=YES
*NODE,NSET=NALL
 1, 0.0, 0.0, 0.0
 2, 2.0, 0.0, 0.0
 3, 2.0, 2.0, 0.0
 4, 0.0, 2.0, 0.0
*NSET,NSET=FAR-FIELD
 3, 4,
*NSET,NSET=NPSI11
 1,
*NSET,NSET=NPSI12
 2,
*NSET,NSET=FARF11
 4,
*NSET,NSET=FARF12
 3,
*ELEMENT,TYPE=B31   ,ELSET=PIPELINE1
 1, 1, 2
*ELEMENT,TYPE=PSI34 ,ELSET=PSI34
 2, 1, 2, 3, 4
*ORIENTATION,NAME=LOCALP,DEFINITION=COORDINATES,SYSTEM=RECTANGULAR
+0.8660254,+0.5000000, 0.0000000,-0.5000000,+0.8660254, 0.0000000
3,0.0
*BEAM SECTION,ELSET=PIPELINE1,SECTION=PIPE,MATERIAL=STEEL
0.3,0.01
*PIPE-SOIL INTERACTION,ELSET=PSI34,ORIENTATION=LOCALP
*PIPE-SOIL STIFFNESS,TYPE=NONLINEAR,DIRECTION=1
-1.000D+03,-0.0020000,100
 0.000D+00, 0.0000000,100
+1.000D+03,+0.0010000,100
-2.000D+03,-0.0020000,200
 0.000D+00, 0.0000000,200
+2.000D+03,+0.0010000,200
*PIPE-SOIL STIFFNESS,TYPE=NONLINEAR,DIRECTION=2
-1.000D+03,-0.0020000,100
-1.500D+03,-0.0015000,100
-1.000D+03,-0.0003333,100
 0.000D+00, 0.0000000,100
+1.000D+03,+0.0005000,100
+2.000D+03,+0.0020000,100
-2.000D+03,-0.0020000,200
-3.000D+03,-0.0015000,200
-2.000D+03,-0.0003333,200
 0.000D+00, 0.0000000,200
+2.000D+03,+0.0005000,200
+4.000D+03,+0.0020000,200
*PIPE-SOIL STIFFNESS,TYPE=NONLINEAR,DIRECTION=3
-3.000D+03,-0.0010000,100
 0.000D+00, 0.0000000,100
+4.000D+03,+0.0010000,100
-6.000D+03,-0.0010000,200
 0.000D+00, 0.0000000,200
+8.000D+03,+0.0010000,200
*MATERIAL,NAME=STEEL
*ELASTIC
2.1E11,0.3
*PLASTIC
2.0E8,0.0
*AMPLITUDE,NAME=CYCLETWO,TIME=STEP TIME
0.0,+0.00000,1.0,+0.06667,2.0,-0.06667,3.0,+0.25000
4.0,-0.25000,5.0,+0.58333,6.0,-0.58333,7.0,+1.00000
8.0,-1.00000
*AMPLITUDE,NAME=RAMP    ,TIME=STEP TIME
0.0,1.0,1.0,2.0
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,100
*STEP,INC=100
*STATIC
0.1,8.0,,0.1
*BOUNDARY
FAR-FIELD,1,3, 0.0 
NPSI11   ,1,1,+3.0E-03
NPSI11   ,2,3, 0.0 
NPSI12   ,1,1,+3.0E-03
NPSI12   ,2,3, 0.0 
NALL     ,4,6, 0.0 
*TEMPERATURE,AMPLITUDE=RAMP
NALL,100.0
*NODE PRINT,FREQ=10
 U,
RF
*EL PRINT,FREQ=10
 S11,S22,S33,E11,E22,E33
 NFORC,
 EE,
 E,
 COORD,
*OUTPUT,FIELD,FREQ=100
*ELEMENT OUTPUT
 S,E,EE,PE
*NODE OUTPUT
 U,RF
*OUTPUT,HISTORY,FREQ=100
*ELEMENT OUTPUT,ELSET=PSI34
 S,E,EE,PE
 NFORC,
*NODE OUTPUT,NSET=NALL
 U,RF
*MONITOR,NODE=1,DOF=1
*END STEP