!
! Example for a curved elasto-plastic spacial beam with ANSYS
! By Lu Xinzheng, Depart. Civil Engineering,
! Tsinghua University, Beijing
! 陸新征，清華大學土木系
! Aug. 2005

R1=5 ! internal radius of the beam
R2=6 ! external radius of the beam
Thick=0.5 ! Thickness of the beam
Fy=200e6 ! Yield strength of concrete
P=1e5 ! Value of pressure load

/prep7

! Define the Element
! 定義單元

ET, 1, Solid45

! Define Material 定義材料
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,200e9
MPDATA,PRXY,1,,.3
TB,BISO,1,1,2,
TBTEMP,0
TBDATA,,Fy,2e9,,,,

! Setup the model

k,1,0,0,0
k,2,0,-R1,
k,3,R1,0
k,4,0,-R2
k,5,R2,0

LARC,2,3,1,R1
LARC,4,5,1,R2
l,2,4
l,3,5

al,1,2,3,4
VEXT,1, , ,0,0,Thick,,,,

! Set the element size
esize,thick/5
vmesh,all

/solu
DA,6,all
! Define the gradient pressure
SFGRAD,PRES, ,X,0,P,

ASEL,S, , , 2
NSLA,S,1
!*
SF,all,PRES, P/2
ALLSEL,ALL

! Define the substeps
NSUBST,20,0,20
OUTRES,ERASE
OUTRES,ALL,1

solve