*HEADING
: BENDING OF CYLINDER DUE TO ASYMMETRIC PORE PRESSURE FIELD  [CAXA8RP3]
** 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=CAXA8RP3,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
*NSET,NSET=BOTTOM,GEN
1,9
1201,1209
2401,2409
3601,3609
*NSET,NSET=L1-000
1,21
*NSET,NSET=L1-060
1201,1221
*NSET,NSET=L1-120
2401,2421
*NSET,NSET=L1-180
3601,3621
*NSET,NSET=L3-000
9,29
*NSET,NSET=L3-060
1209,1229
*NSET,NSET=L3-120
2409,2429
*NSET,NSET=L3-180
3609,3629
*NSET,NSET=L1
L1-000,L1-060,L1-120,L1-180
*NSET,NSET=L3
L3-000,L3-060,L3-120,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-060,8,,0.5E6
L1-120,8,,-0.5E6
L1-180,8,,-1.E6
L3-000,8,,3.E6
L3-060,8,,1.5E6
L3-120,8,,-1.5E6
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
*END STEP