*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
*NGEN,NSET=BOT1
*NGEN,NSET=TOP1
*NFILL,NSET=P001
BOT1,TOP1,2,10
*NCOPY,NEWSET=P234,OLDSET=P001,SHIFT,CHANGENUMBER=100,MULTIPLE=24
*NSET,NSET=NALL
P001,P234
*NSET,NSET=YSYMM,GEN
*ELEMENT,TYPE=C3D20P,ELSET=ALL
*ELGEN,ELSET=ALL
*TRANSFORM,NSET=NALL,TYPE=C
*ORIENTATION,NAME=CYLORI,SYS=CYL
*SOLID SECTION,ELSET=ALL,MATERIAL=MATPROP,ORIENTATION=CYLORI
*MATERIAL,NAME=MATPROP
*ELASTIC
 1.E8,0.3
*PERMEABILITY
 1.E-5, 1.
*INITIAL CONDITIONS,TYPE=RATIO
NALL,1.
*NSET,NSET=BOTTOM,GEN
*NSET,NSET=L1-000
*NSET,NSET=L1-015
*NSET,NSET=L1-030
*NSET,NSET=L1-045
*NSET,NSET=L1-060
*NSET,NSET=L1-075
*NSET,NSET=L1-090
*NSET,NSET=L1-105
*NSET,NSET=L1-120
*NSET,NSET=L1-135
*NSET,NSET=L1-150
*NSET,NSET=L1-165
*NSET,NSET=L1-180
*NSET,NSET=L3-000
*NSET,NSET=L3-015
*NSET,NSET=L3-030
*NSET,NSET=L3-045
*NSET,NSET=L3-060
*NSET,NSET=L3-075
*NSET,NSET=L3-090
*NSET,NSET=L3-105
*NSET,NSET=L3-120
*NSET,NSET=L3-135
*NSET,NSET=L3-150
*NSET,NSET=L3-165
*NSET,NSET=L3-180
*NSET,NSET=L1
L1-000,L1-180
*NSET,NSET=L3
L3-000,L3-180
*NSET,NSET=NFILE,GEN
*ELSET,ELSET=EFILE
*BOUNDARY
1,1
BOTTOM,3
YSYMM,2
*RESTART,WRITE,FREQ=999
*STEP
*SOILS
*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