! Analysis a box beam section with shell elements of ANSYS
! 用ansys的殼體單元分析箱梁
! Box dimension: 10*4*4m with shell thickness of 0.04m
! By Lu Xinzheng, Depart. Civil Engineering, Tsinghua University, Beijing
! 陸新征，清華大學土木系
! Aug. 2004
! Define the Element
! 定義單元
/PREP7
!*
ET,1,SHELL93
!*
! Define the section for shells
! 定義殼單元截面
R,1,.04, , , , , ,
!*
! Define the material
! 定義材料
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,200e9
MPDATA,PRXY,1,,.3
! Setup the model
! 建模
! Define the keypoints
! 定義關鍵點
k,1,0,0,0
k,2,4,0,0
k,3,4,4,0
k,4,0,4,0
! Define the lines
! 定義線
l,1,2
l,2,3
l,3,4
l,4,1
! Define the section area
! 定義截面
a,1,2,3,4
! Extrude a volumn from area
! 從面拉伸得到體
VEXT,1, , ,0,0,10,,,,
! Delete useless volumn and areas
! 刪除不必要的體和面
VDELE, 1
ADELE,1,2,1
! Mesh
! 劃分網格
ESIZE,0.5,0,
AMESH,ALL
!*
FINISH
/SOL
! Fixed end
! 固定端部
DL,1, ,ALL,
DL,2, ,ALL,
DL,3, ,ALL,
DL,4, ,ALL,
! Give the pressure on the top surface
! 頂部施加壓力
SFA,5,2,PRES,1

SOLVE
FINISH

/POST1
! Display the deformation
! 繪制變形
PLDISP,0
! Select the mid-section of the shells
! 選定殼中部截面
shell,mid
! Define a path to plot the stress
! 定義一個路徑來繪制應力
PATH,MidPath1,4,30,20,
! Define path with nodes
! 用節點來定義路徑
PPATH,1,89
PPATH,2,34
PPATH,3,570
PPATH,4,1066
PATH,STAT
!*
AVPRIN,0,0,
!*
! Map stress Z to the path
! 將Stress Z映射到路徑
PDEF, ,S,Z,AVG
/PBC,PATH, ,0
! Plot the stress
! 繪制應力
PLPATH,SZ
FINISH