! 5.4  圓柱形薄殼在均勻內(nèi)壓作用下的靜力分析
! 本程序來(lái)源于邢靜忠編著的《ANSYS應(yīng)用實(shí)例與分析》，科學(xué)出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除數(shù)據(jù)，不加載初始化文件
/FILNAME, EX5.4	! 工程工件名稱(chēng)
/TITLE, EX5.4, CYLINDRICAL MEMBRANE UNDER PRESSURE.
! （1）定義工程選項(xiàng)、分析類(lèi)型、單元類(lèi)型、實(shí)常數(shù)和材料參數(shù)
/PREP7	! 加載前處理模塊
ANTYPE, STATIC	! 定義分析類(lèi)型為靜力分析
ET, 1, SHELL41	! 定義第1類(lèi)單元薄殼單元SHELL41
R, 1, 25.4	! 定義厚度25.4 mm
MP, EX, 1, 207E3	! 定義第1類(lèi)材料的彈性模量EX=207 GPa
MP, NUXY, 1, 0.3	! 定義第1類(lèi)材料的泊松比NUXY=0.3
! （2）在柱坐標(biāo)系下, 定義節(jié)點(diǎn)和單元
CSYS, 1	! 切換為柱坐標(biāo)系
/View, 1, 1, 1, 1	! 將窗口1切換到等軸側(cè)視點(diǎn)方向
N, 1, 1524 	! 在(1524mm, 0, 0)位置定義節(jié)點(diǎn)1
N, 2, 1524, , 254	! 在(1524mm, 0, 254mm)位置定義節(jié)點(diǎn)
NGEN, 2, 2, 1, 2, 1, , 10	! 按照1, 2號(hào)節(jié)點(diǎn), 循環(huán)2次, 節(jié)點(diǎn)號(hào)加2, Y坐標(biāo)加10度, 生成3, 4節(jié)點(diǎn)
NROTAT, ALL	! 旋轉(zhuǎn)節(jié)點(diǎn)坐標(biāo)系到柱坐標(biāo)系方向
E, 1, 2, 4, 3	! 過(guò)節(jié)點(diǎn)1, 2, 4, 3定義單元
! （3）定義耦合自由度、位移約束和荷載力
CP, 1, UX, 1, 2, 3, 4	! 將1, 2, 3, 4號(hào)節(jié)點(diǎn)的UX徑向自由度定義為第1個(gè)耦合自由度
CP, 2, UZ, 2, 4	! 將2, 4號(hào)節(jié)點(diǎn)的UZ軸向自由度定義為第2個(gè)耦合自由度
D, 1, UZ, , , 3, 2	! 位移約束1號(hào)節(jié)點(diǎn)的UZ自由度, 按增量2循環(huán)到3號(hào)節(jié)點(diǎn)
D, ALL, UY  	! 位移約束所有節(jié)點(diǎn)的UY自由度
SFE, 1, 4, PRES, , -103.43*25.4	! 定義1號(hào)單元的4號(hào)面上的軸向拉應(yīng)力103.43 MPa
SFE, 1, 1, PRES, , -3.447575	! 定義內(nèi)壓3.447575 MPa
FINISH	! 退出前處理模塊PREP7
! (4)在求解模塊, 設(shè)置求解選項(xiàng)后求解
/SOLU	! 加載求解模塊
OUTPR, , 1	! 輸出第1荷載步的基本結(jié)果
SOLVE	! 開(kāi)始求解
FINISH	! 退出求解模塊SOLUTION
! (5)在后處理模塊中提取應(yīng)力計(jì)算結(jié)果
/POST1	! 進(jìn)入一般后處理模塊
*GET, STRS_HOP, NODE, 1, S, 2	! 提取1號(hào)節(jié)點(diǎn)的第2主應(yīng)力, 存入變量STRS_HOP
*GET, STRS_AX, NODE, 1, S, 1	! 提取1號(hào)節(jié)點(diǎn)的第1主應(yīng)力, 存入變量STRS_AX
FINISH