Dear Members

I am trying to perform a static analysis of 3D cylinder subjected to an internal pressure using Ansys APDL. in the next, the code i did but it seems it didn't work properly

could anyone help

/title, static cylinder 

fini

/clear

! MATERIAL BOT

*set, E_m, 207.78e9 

*set, nu_m, 0.3177

*set,dens_m,8166 ! [Kg/m3]

cylin_len=1 ! [m] length

cylin_radout=cylin_len/3 ! [m] outer radius

cylin_radin=cylin_radout/2 ! [m] inner radius

R1 = cylin_radin ! [m] inner radius

R2 = cylin_radout ! [m] outer radius

diff=(R2-R1)/10

Z1 = 0

Z2 = cylin_len

THETA1 = 0

THETA2 = 180

LENGTH = Z2-Z1 ! LENGTH OF THE FLUID ELEMENT (M)

DIAMETER = 2*R1 ! DIAMETER (M)

/PREP7

ET,1,SOLID185

mp, ex, 1, E_m

mp, prxy, 1, nu_m

mp, dens, 1, den_m

K,1,0,0,Z1

K,2,0,0,Z2

L,1,2

LESIZE,1,,,3

CYLIND,R1,R2,Z1,Z2,THETA1,THETA2

CYLIND,R1,R2,Z1,Z2,THETA1+180,THETA2+180

VGLUE,1,2

MSHAP,TRIOPT,3D

LSEL,S,LINE,,2,4,2

LSEL,A,LINE,,7,9,2

LESIZE,ALL,,,10

LSEL,S,LINE,,10,13,1

LESIZE,ALL,,,18

LSEL,S,LINE,,6,8,2

LSEL,A,LINE,,28,29,1

LSEL,A,LINE,,3,5,2

LSEL,A,LINE,,26,27,1

LESIZE,ALL,,,24

LDELE, 1

TYPE,1

MAT,1

REAL,1

VMESH,ALL

ALLSEL,ALL

/SOLU

antype,static

p=2e9

ASEL,S,LOC,x,R1

SFA,all,1,PRES,P

NSEL,ALL

OUTPR,,ALL

SOLVE

save ! LOAD STEP 1 - INTERNAL PRESSURE

FINISH