! 4.14  兩端固定的厚壁管道在自重作用下的變形和應(yīng)力
! 本程序來源于邢靜忠編著的《ANSYS應(yīng)用實例與分析》，科學出版社，2006年
!
FINISH
/CLEAR, NOSTART
! （1）設(shè)置工程參數(shù), 分析類型、單元類型、實常數(shù)和材料參數(shù)
/FILNAME, EX4.14
/PREP7
SMRT, OFF
/TITLE, EX4.14, BENDING OF AXISYMMETRIC THICK PIPE UNDER GRAVITY LOADING.
ANTYPE, STATIC
ET, 1, PLANE25, , , , , , 2
MP, EX, 1, 206.85E9
MP, DENS, 1, 7802
MP, NUXY, 1, 0
! (2)定義軸對稱截面上的關(guān)鍵點和線
K, 1, 0.0127
K, 2, 0.0127, 2.54
KGEN, 2, 1, 2, 1, 0.0127
L, 1, 2
LESIZE, 1, , , 24
L, 2, 4
LESIZE, 2, , , 1
L, 3, 4
LESIZE, 3, , , 24
L, 1, 3
LESIZE, 4, , , 1
A, 3, 1, 2, 4
AMESH, 1
/VIEW, 1, 0, 0, 1
/VUP, 1, X
/ANGLE, 1, 0
EPLOT
! （3）定義重力加速度和位移約束條件
ACEL, 9.804, , -9.804
MODE, 1, 1
NSEL, S, LOC, Y, 0
D, ALL, ALL
NSEL, S, LOC, Y, 2.54
D, ALL, UY
NSEL, ALL
FINISH
! (4)在求解模塊, 設(shè)置求解選項后求解
/SOLU
OUTPR, BASIC, LAST
SOLVE
FINISH
! (5)在后處理模塊中顯示第1個荷載步的0度位置的位移計算結(jié)果
/POST1
SET, 1, 1, , , , 0.0
/VUP, 1, X
/WINDOW, 1, -1.1, 2, 0, 2.5
! PLDISP, 1
PRNSOL, U, COMP
*GET, DEF_X, NODE, 3, U, X
NODE1=NODE(0.0254, 0.212, 0)
*GET, STRS_Y0, NODE, NODE1, S, Y
PLNSOL, S, Y
! (6)處理第1個荷載步的90度位置的計算結(jié)果
SET, 1, 1, , , , 90.0
/WINDOW, 1, OFF
/NOERASE
/WINDOW, 2, -1, 2, -2.5, 0
/VUP, 2, X
PLNSOL, S, Y
PRNSOL, U, COMP
*GET, DEF_Z, NODE, 3, U, Z
NODE1=NODE(0.0254, 0.212, 0)
*GET, STRS_Y90, NODE, NODE1, S, Y
*STATUS, PARM
FINISH