*HEADING
: BENDING OF CYLINDER DUE TO ASYMMETRIC PORE PRESSURE FIELD  [CAXA8RP4]
** 6-IN LONG CYLINDER OF 2-IN & 6-IN INNER AND OUTER RADII (2X1 MESH)
*NODE
1,2.,0.
9,6.,0.
21,2.,6.
29,6.,6.
*NGEN,NSET=BOT1
1,9
*NGEN,NSET=TOP1
21,29
*NFILL,NSET=P001
BOT1,TOP1,2,10
*ELEMENT,TYPE=CAXA8RP4,OFFSET=1200,ELSET=ALL
1,1,5,25,21,3,15,23,11
*ELGEN,ELSET=ALL
1,2,4,1
*SOLID SECTION,ELSET=ALL,MATERIAL=MATPROP
*MATERIAL,NAME=MATPROP
*ELASTIC
 1.E8,0.3
*PERMEABILITY,SPECIFIC=1.0
 1.E-5, 1.
*INITIAL CONDITIONS,TYPE=RATIO
NALL,1.
*NSET,NSET=NALL,GEN
1,29
1201,1229
2401,2429
3601,3629
4801,4829
*NSET,NSET=BOTTOM,GEN
1,9
1201,1209
2401,2409
3601,3609
4801,4809
*NSET,NSET=L1-000
1,21
*NSET,NSET=L1-045
1201,1221
*NSET,NSET=L1-090
2401,2421
*NSET,NSET=L1-135
3601,3621
*NSET,NSET=L1-180
4801,4821
*NSET,NSET=L3-000
9,29
*NSET,NSET=L3-045
1209,1229
*NSET,NSET=L3-090
2409,2429
*NSET,NSET=L3-135
3609,3629
*NSET,NSET=L3-180
4809,4829
*NSET,NSET=L1
L1-000,L1-045,L1-090,L1-135,L1-180
*NSET,NSET=L3
L3-000,L3-045,L3-090,L3-135,L3-180
*NSET,NSET=NFILE
NALL, 
*ELSET,ELSET=EFILE
1,2
*BOUNDARY
1,1
BOTTOM,2,2
*RESTART,WRITE,FREQ=999
*STEP
*SOILS
 1., 1.
*BOUNDARY
L1-000,8,,1.E6
L1-045,8,,0.707107E6
L1-090,8,,0.0
L1-135,8,,-0.707107E6
L1-180,8,,-1.E6
L3-000,8,,3.E6
L3-045,8,,2.121321E6
L1-090,8,,0.0
L3-135,8,,-2.121321E6
L3-180,8,,-3.E6
*EL PRINT,ELSET=EFILE
S, 
*EL PRINT,ELSET=EFILE,POSITION=AVERAGED AT NODES
S, 
E, 
*NODE PRINT,NSET=NALL,TOTAL=NO
U,POR,RVT,RF
** *NODE PRINT,NSET=L1,TOTAL=NO
** COORD,U,POR,RVT,RF
** *NODE PRINT,NSET=L5,TOTAL=NO
** COORD,U,POR,RVT,RF
** *NODE PRINT,NSET=L9,TOTAL=NO
** COORD,U,POR,RVT,RF
*EL FILE,ELSET=EFILE,POSITION=AVERAGED AT NODES
E, 
*NODE FILE,NSET=NFILE
U,POR,RVT,RF
*OUTPUT,FIELD
*NODE OUTPUT,NSET=NFILE
U,POR,RVT,RF
*END STEP