*HEADING
: BENDING OF CYLINDER DUE TO ASYMMETRIC PORE PRESSURE FIELD  [C3D20P]
** 6-IN LONG CYLINDER OF 2-IN & 6-IN INNER AND OUTER RADII (2X6 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
*NCOPY,NEWSET=P234,OLDSET=P001,SHIFT,CHANGENUMBER=100,MULTIPLE=24

0.,0.,0., 0.,0.,1., 7.5
*NSET,NSET=NALL
P001,P234
*NSET,NSET=YSYMM,GEN
1,9
11,19
21,29
2401,2409
2411,2419
2421,2429
*ELEMENT,TYPE=C3D20P,ELSET=ALL
1,1,21,25,5,201,221,225,205, 11,23,15,3,211,223,215,
203, 101,121,125,105
*ELGEN,ELSET=ALL
1,2,4,1,12,200,2
*TRANSFORM,NSET=NALL,TYPE=C
0.,0.,0., 0.,0.,1.
*ORIENTATION,NAME=CYLORI,SYS=CYL
0.,0.,0., 0.,0.,1.
3,0.
*SOLID SECTION,ELSET=ALL,MATERIAL=MATPROP,ORIENTATION=CYLORI
*MATERIAL,NAME=MATPROP
*ELASTIC
 1.E8,0.3
*PERMEABILITY,SPECIFIC=1.0
 1.E-5, 1.
*INITIAL CONDITIONS,TYPE=RATIO
NALL,1.
*NSET,NSET=BOTTOM,GEN
1,2401,100
3,2403,200
5,2405,100
7,2407,200
9,2409,100
*NSET,NSET=L1-000
1,21
*NSET,NSET=L1-015
201,221       
*NSET,NSET=L1-030
401,421
*NSET,NSET=L1-045
601,621
*NSET,NSET=L1-060
801,821
*NSET,NSET=L1-075
1001,1021
*NSET,NSET=L1-090
1201,1221
*NSET,NSET=L1-105
1401,1421
*NSET,NSET=L1-120
1601,1621
*NSET,NSET=L1-135
1801,1821
*NSET,NSET=L1-150
2001,2021
*NSET,NSET=L1-165
2201,2221
*NSET,NSET=L1-180
2401,2421
*NSET,NSET=L3-000
9,29,10
*NSET,NSET=L3-015
209,229
*NSET,NSET=L3-030
409,429
*NSET,NSET=L3-045
609,629
*NSET,NSET=L3-060
809,829
*NSET,NSET=L3-075
1009,1029
*NSET,NSET=L3-090
1209,1229
*NSET,NSET=L3-105
1409,1429
*NSET,NSET=L3-120
1609,1629
*NSET,NSET=L3-135
1809,1829
*NSET,NSET=L3-150
2009,2029
*NSET,NSET=L3-165
2209,2229
*NSET,NSET=L3-180
2409,2429
*NSET,NSET=L1
L1-000,L1-180
*NSET,NSET=L3
L3-000,L3-180
*NSET,NSET=NFILE,GEN
1,9
11,19
21,29
1201,1209
1211,1219
1221,1229
2401,2409
2411,2419
2421,2429
*ELSET,ELSET=EFILE
1,2,23,24
*BOUNDARY
1,1
BOTTOM,3
YSYMM,2
*RESTART,WRITE,FREQ=999
*STEP
*SOILS
 1., 1.
*BOUNDARY
L1-000,8,,1.0E6 
L1-015,8,,0.965926E6 
L1-030,8,,0.866E6 
L1-045,8,,0.707107E6 
L1-060,8,,0.5E6 
L1-075,8,,0.258819E6 
L1-090,8,,0.0
L1-105,8,,-0.258819E6 
L1-120,8,,-0.5E6 
L1-135,8,,-0.707107E6 
L1-150,8,,-0.866E6 
L1-165,8,,-0.965926E6 
L1-180,8,,-1.0E6 
L3-000,8,,3.0E6 
L3-015,8,,2.897778E6 
L3-030,8,,2.598E6 
L3-045,8,,2.121321E6 
L3-060,8,,1.5E6 
L3-075,8,,0.776457E6 
L3-090,8,,0.0
L3-105,8,,-0.776457E6 
L3-120,8,,-1.5E6 
L3-135,8,,-2.121321E6  
L3-150,8,,-2.598E6 
L3-165,8,,-2.897778E6 
L3-180,8,,-3.0E6 
*EL PRINT,ELSET=EFILE
S, 
*EL PRINT,ELSET=EFILE,POSITION=AVERAGED AT NODES
S, 
E, 
*NODE PRINT,NSET=NFILE,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
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=EFILE
S,
*OUTPUT,FIELD
*NODE OUTPUT,NSET=NFILE
U,POR,RVT,RF
*END STEP