! 3.4  工字形截面外伸梁的平面彎曲
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!
Finish
/CLEAR, NOSTART	! 清除數據，不加載初始化文件
/FILNAME, EX3.4	! 定義文件名稱
/PREP7	! 加載前處理模塊
/TITLE, EX3.4, BEAM STRESSES AND DEFLECTIONS.	! 設置標題
! （1）定義求解類型，實常數和材料參數
ANTYPE, STATIC	! 定義分析類型為靜力分析
ET, 1, BEAM3	! 定義第1類單元為二維彈性梁單元BEAM3
KEYOPT, 1, 9, 9	! 在9個中點位置上輸出計算結果
R, 1, 35e-4, 4.5E-5, 0.20	! 定義第1類實常數、截面積35cm^2、慣性矩4500cm^4和高度20cm
MP, EX, 1, 210E9	! 定義第1類材料的彈性模量EX＝210Gpa
MP, PRXY, 1, 0.3	! 定義第1類材料的泊松比PRXY=0.3
! （2）定義節點位置
N, 1	! 在坐標原點(0，0)定義1號節點
N, 5, 4	! 在(4m，0)位置定義5號節點
FILL	! 利用節點填充FILL生成中間3個節點
! （3）定義單元連接關系
E, 1, 2	! 過1，2節點定義第1個單元
EGEN, 4, 1, 1	! 按照前面的單元模式，生成4次，每次兩端節點編號加1，生成其余3個單元
! （4）定義位移約束和荷載條件
D, 2, UX, , , , , UY	! 定義2號節點的位移邊界條件，UX=0, UY=0
D, 4, UY	! 定義4號節點的位移邊界條件，UY=0
SFBEAM, 1, 1, PRES, 5000	! 定義第1個單元的表面荷載5000N/m
SFBEAM, 4, 1, PRES, 5000	! 定義第4個單元的表面荷載5000N/m
FINISH	! 退出前處理模塊PREP7
! （5）進入求解模塊，定義求解選項并求解
/SOLU	! 加載求解模塊
OUTPR, BASIC, 1	! 只輸出第1個荷載步的基本計算結果
SOLVE	! 開始求解
FINISH	! 退出求解模塊SOLU
! （6）進入后處理模塊，列表顯示節點位移計算結果和圖形顯示變形圖
/POST1	! 進入一般后處理模塊
SET, 1, 1	! 讀入第1個荷載步的計算結果
PRNSOL, U, COMP	! 列表顯示節點線位移和節點線位移的矢量和
PRNSOL, ROT, COMP	! 列表顯示節點轉角和節點角位移的矢量和
PLDISP, 1	! 帶原來模型圖，顯示變形后的模型
! （7）提取指定位置的節點和單元，定義單元表并獲取彎曲應力
MID_NODE = NODE (2, , , )	! 選擇距離(2, 0, 0)位置最近的節點，命名為MID_NODE
*GET, DISP, NODE, MID_NODE, U, Y	! 提取名稱為MID_NODE的節點的豎向位移UY，存入變量DISP
MID_ELM = ENEARN (MID_NODE)	! 選擇距離MID_NODE節點最近的單元, 命名為MID_ELM
ETABLE, STRS, LS, 3	! 以LS, 3（-Y一側的彎曲應力SBYB）為內容定義單元表STRS
*GET, STRSS, ELEM, MID_ELM, ETAB, STRS	! 提取單元表STRS中單元MID_ELM的應力，存入STRSS
*STATUS, PARM
FINISH