! 2.1  鉸接桿在外力作用下的變形
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除已有的數(shù)據(jù)
! （1）設置工程選項, 定義幾何參數(shù)
/FilName, EX2.1	! 定義工程文件名稱
/TITLE, EX2.1, DEFLECTION OF A HINGED SUPPORT
*AFUN, DEG 	! 將三角函數(shù)的默認單位由弧度改為角度
THETA=30  	! 定義變量Theta為桿件與水平面夾角30度
A=10000	! 定義變量A為節(jié)點3的水平X坐標10000mm
B=A/2*TAN(THETA)	! 節(jié)點2的鉛直Y坐標
! （2）進入前處理模塊, 定義節(jié)點
/PREP7	! 進入前處理模塊
N, 1	! 在坐標原點(0, 0, 0)定義第1號節(jié)點
N, 2, A/2, -B	! 在 (A/2, -B) 定義第2號節(jié)點
N, 3, A	! 在 (A, 0) 定義第3號節(jié)點
! （3）在前處理模塊中, 定義單元類型、材料參數(shù)和各個單元
ET, 1, LINK1	! 定義第1類單元為二維桿單元LINK1
R, 1, 1000	! 定義第1類實常數(shù), 桿件截面面積A=1000mm^2
MP, EX, 1, 210E3	! 定義第1類材料的彈性模量EX＝210000N/mm^2
E, 1, 2	! 過節(jié)點1, 2定義單元1
E, 2, 3	! 過節(jié)點2, 3定義單元2
! （4）在前處理模塊中, 定義位移約束
D, 1, ALL, , , 3, 2	! 約束1號節(jié)點的所有位移自由度, 按增量2循環(huán)到3來約束3號節(jié)點
F, 2, FY, -1000	! 在2號節(jié)點施加沿著-Y方向大小為1000N的集中力
OUTPR, , 1	! 輸出第1個荷載步的基本計算結果
FINISH	! 退出前處理模塊PREP7
! （5）進入求解模塊, 開始求解
/SOLU	! 進入求解模塊
SOLVE	! 開始靜力分析求解
FINISH	! 退出求解模塊SOLU
! （6）進入后處理模塊, 顯示節(jié)點位移和桿件內力
/POST1	! 進入后處理模塊
PLDISP, 2	! 顯示結構變形圖, 參數(shù)“2”表示用虛線繪制出原來結構輪廓
MID_NODE = NODE (A/2, -B, 0 )	! 尋找距離 (A/2, -B, 0 )最近的節(jié)點, 存入MID_NODE
*GET, DISP, NODE, MID_NODE, U, Y	! 提取節(jié)點MID_NODE上的位移UY, 存入DISP
LEFT_EL = ENEARN (MID_NODE)	! 尋找距離節(jié)點MID_NODE最近的單元, 存入LEFT_EL
ETABLE, STRS, LS, 1	! 用軸向應力SAXL的編號”LS, 1”定義單元表STRS
*GET, STRSS, ELEM, LEFT_EL, ETAB, STRS	! 從單元表STRS中提取LEFT_EL單元的應力, 存入STRSS
*STATUS, PARM	! 顯示所有內存變量的值
FINISH	! 退出后處理模塊