! 3.11  偏心受壓桿的大變形分析
! 本程序來源于邢靜忠編著的《ANSYS應(yīng)用實(shí)例與分析》，科學(xué)出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除數(shù)據(jù)，不加載初始化文件
/FILNAME, EX3.11	! 定義文件名稱
! （1）定義工程選項(xiàng)、分析類型、實(shí)常數(shù)和材料參數(shù)等
/PREP7	! 加載前處理模塊
/TITLE, EX3.11, LARGE DEFLECTION ECCENTRIC COMPRESSION OF SLENDER COLUMN.
ANTYPE, STATIC	! 定義分析類型為靜力分析
NLGEOM, ON	! 激活大變形選擇項(xiàng)
ET, 1, BEAM54	! 定義第1類單元為二維錐形非對(duì)稱梁?jiǎn)卧狟EAM54
R, 1, 2167.7, 54.106e4, 14.732, 42.671	! 定義第1類實(shí)常數(shù)：截面積2167.7mm^2，抗彎慣性矩54.106cm^4，中性軸到上下邊的距離分別為14.732mm和42.671mm
RMORE, , , , -14.732, , -14.732	! 兩端偏心距e=14.732mm（HYT2, HYB2, DX1, DY1, DX2, DY2）
MP, EX, 1, 207E3	! 定義第1類材料的彈性模量EX = 207e3 N/mm^2
MP, NUXY, 1, 0.3	! 定義第1類材料的泊松比NUXY = 0.3
! （2）定義節(jié)點(diǎn)
N, 1	! 在坐標(biāo)原點(diǎn)(0, 0, 0)定義第1號(hào)節(jié)點(diǎn)
N, 5, , 1524	! 在坐標(biāo)(0, 1524mm, 0)定義第5號(hào)節(jié)點(diǎn)
FILL	! 在第1號(hào)節(jié)點(diǎn)和第5號(hào)節(jié)點(diǎn)中間填充生成其余節(jié)點(diǎn)
! （3）定義單元
E, 1, 2	! 過節(jié)點(diǎn)1和節(jié)點(diǎn)2定義單元
EGEN, 4, 1, 1	! 用前面一個(gè)單元的模式，循環(huán)4次，生成其余的3個(gè)單元
! （4）施加位移約束和荷載
D, 1, ALL	! 完全固定1號(hào)節(jié)點(diǎn)
F, 5, FY, -17793	! 定義沿著-Y軸的大小為17793N的集中力
FINISH	! 退出前處理模塊PREP7
! （5）進(jìn)入求解模塊求解
/SOLU	! 加載求解模塊
CNVTOL, F, , 1E-4	! 設(shè)置非線性分析力的收斂誤差1E-4 N
CNVTOL, M, , 1E-1	! 設(shè)置非線性分析力偶的收斂誤差為1E-1 N.m
OUTPR, , LAST	! 設(shè)置只輸出最后一個(gè)荷載步的結(jié)果
SOLVE	! 開始大變形分析的非線性迭代求解
FINISH	! 退出求解模塊SOLUTION
! （6）在后處理模塊，提取計(jì)算結(jié)果
/POST1	! 進(jìn)入一般后處理模塊
END_NODE = NODE (0, 3048/2, 0)	! 選擇距離(0, 1524mm, 0)最近的節(jié)點(diǎn), 命名為END_NODE
*GET, DEF, NODE, END_NODE, U, X	! 提取節(jié)點(diǎn)END_NODE的位移UX, 存入變量DEF
ETABLE, TENS, NMISC, 1	! 以BEAM54的最大拉應(yīng)力(NMISC, 1)為內(nèi)容定義單元表TENS
ETABLE, COMP, NMISC, 2	! 以BEAM54的最大壓應(yīng)力(NMISC, 2)為內(nèi)容定義單元表COMP
*GET, STS_TENS, ELEM, 1, ETAB, TENS	! 提取第1號(hào)單元的最大拉應(yīng)力, 存入變量STS_TENS
*GET, STS_COMP, ELEM, 1, ETAB, COMP	! 提取第1號(hào)單元的最大壓應(yīng)力, 存入變量STS_COMP
*STATUS, PARM
FINISH