! 4.10  在端部集中力下的變截面懸臂梁
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!
FINISH
/CLEAR, NOSTART
! (1)設置工程選項，定義分析類型、單元類型和材料參數
/FILNAME, EX4.10(1)
/PREP7
/TITLE, EX4.10(1), LATERALLY LOADED TAPERED SUPPORT STRUCTURE (QUAD. ELEMENTS)
! 使用 PLANE42 單元
ANTYPE, STATIC
ET, 1, PLANE42, , , 3, , , 2 以提取中間節點的結果
R, 1, 50.8
MP, EX, 1, 207E3
MP, NUXY, 1, 0.0
! (2)定義節點位置和單元連接
N, 1, 635
N, 7, 1905
FILL
N, 8, 635, -76.2
N, 14, 1905, -228.6
FILL
E, 2, 1, 8, 9
EGEN, 6, 1, 1
! (3)選擇節點，施加位移邊界條件和荷載
NSEL, S, LOC, X, 1905
D, ALL, ALL
NSEL, ALL
F, 1, FY, -17793
FINISH
! (4)在求解模塊，開始求解
/SOLU
OUTPR, , 1
SOLVE
FINISH
! (5)提取指定位置節點的應力，繪制變形圖
/POST1
END_NODE = NODE (1905, 0, 0)
*GET, STS_E_42, NODE, END_NODE, S, X
PLDISP, 2
MID_NODE = NODE (1270, 0, 0)
*GET, STS_M_42, NODE, MID_NODE, S, EQV
PLNSOL, S, X
PLNSOL, S, 1
PLNSOL, S, 3
FINISH
! (6)改用更高階的平面應力單元PLANE82重新計算
/PREP7
/TITLE, EX4.10(2), LATERALLY LOADED TAPERED SUPPORT STRUCTURE (PLANE82).
ET, 1, PLANE82, , , 3, , , 2 間節點的結果
EMID
NSEL, R, LOC, X, 1905
NSEL, R, LOC, Y, -114.3
D, ALL, ALL
NSEL, ALL
FINISH
! (7)在求解模塊求解
/SOLU
SOLVE
FINISH
! (8)提取指定位置的節點的應力計算結果
/POST1
END_NODE = NODE (1905, 0, 0)
*GET, STS_E_82, NODE, END_NODE, S, X
MID_NODE = NODE (1270, 0, 0)
*GET, STS_M_82, NODE, MID_NODE, S, EQV
PLNSOL, S, X
PLNSOL, S, 1
PLNSOL, S, 3
FINISH