! 3.9  懸臂等強度梁的彎曲
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除數據，不加載初始化文件
/FILNAME, EX3.9	! 定義工程文件名
! （1）定義工程選項，分析類型為靜力（STATIC）分析，單元類型和實常數，材料參數等
/PREP7	! 加載前處理模塊
/TITLE, EX3.9(1), BENDING OF A TAPERED PLATE (BEAM) WITH SHELL63 ELEMENT.
ANTYPE, STATIC	! 定義分析類型為靜力分析
ET, 1, SHELL63, 2	! 定義第1類單元為考慮彎曲效應的彈性殼體單元SHELL63
R, 1, 12.7	! 定義第1類實常數，殼體厚度 12.7 mm
MP, EX, 1, 207E3	! 定義第1類材料的彈性模量為207E3 N/mm^2
MP, NUXY, 1, 0	! 定義泊松比NUXY=0
! （2）通過節點和單元建模
N, 1	! 在坐標原點(0, 0, 0)定義第1號節點
N, 8, 508, -38.1	! 在(508mm, -38.1mm, 0)位置定義第8號節點
FILL	! 利用填充生成方法在前面定義的1號節點和8號節點生成中間的6個節點
N, 11	! 在坐標原點(0, 0, 0)定義第11號節點
N, 18, 508, 38.1	! 在(508mm, 38.1mm, 0)位置定義第18號節點
FILL	! 利用填充生成方法在前面定義的11號節點和18號節點生成中間的6個節點
E, 1, 2, 12	! 過1，2，12號節點定義單元
E, 2, 3, 12	! 過2，3，12號節點定義單元
E, 13, 12, 3	! 過13，12，3號節點定義單元
E, 3, 4, 14	! 過3，4，14號節點定義單元
E, 14, 13, 3	! 過14，13，3號節點定義單元
EGEN, 3, 2, 2, 5	! 循環3次，每次節點號增量為2，按照前面定義的第2個到第5個單元的模式
! （3）定義耦合自由度，以保證變形對稱
CP, 1, UZ, 2, 12	! 定義2號節點和12號節點的UZ自由度為耦合自由度，定義為集合1
CP, 2, ROTY, 2, 12	! 定義2號節點和12號節點的ROTY自由度為耦合自由度，定義為集合2
CPSGEN, 6, 1, 1, 2	! 按照前定義耦合自由度集合1到2，循環6次生成其余的耦合自由度集合
OUTPR, ALL, ALL	! 設置輸出所有荷載步的所有計算結果
D, 8, ALL, , , 18, 10	! 約束8號節點的所有位移自由度，按增量10循環到18，約束18號節點
D, ALL, ROTX, 0	! 約束所有節點的繞X軸的轉動自由度
F, 1, FZ, -44.483	! 給1號節點施加沿著-Z方向大小為44.483N的集中力
FINISH	! 退出前處理模塊
! （4）進入求解模塊，開始求解
/SOLU	! 加載求解模塊
SOLVE	! 開始求解
FINISH	! 退出求解模塊
! （5）進入一般后處理模塊，繪制變形圖定義單元表，提取應力結果和節點位移
/POST1	! 進入一般后處理模塊
/VIEW, 1, 1, 1, 1	! 切換到等軸測視點
PLDISP, 2	! 繪制變形圖
ETABLE, STRS, S, 1	! 用單元SHELL63的最大拉應力S1為內容定義單元表STRS
ESORT, STRS	! 按最大拉應力對單元表STRS排序
*GET, SMAX1, SORT, , MAX	! 獲取排序后的最大拉應力值, 存入變量SMAX1
PRNSOL, DOF	! 列表顯示節點位移計算結果
LFT_NODE = NODE(0, 0, 0)	! 提取距離(0, 0, 0)最近的節點號，存入LFT_NODE
*GET, DEFL1, NODE, LFT_NODE, U, Z	! 獲得節點LFT_NODE的位移UZ, 存入變量DEFL1
*STATUS, PARM	! 顯示變量中第4種情況的結果
! （6）將計算結果保存到數據庫
SAVE, TABLE_1	! 保存數據到數據庫TABLE_1
FINISH	! 退出后處理模塊
! （7）改為三維梁單元BEAM44，重新計算
/CLEAR, NOSTART	! 清除數據，不加載初始化文件
/PREP7	! 加載前處理模塊
/TITLE, EX3.9(2), BENDING OF A TAPERED PLATE WITH TAPERED BEAM ELEMENT BEAM44.
! （8）定義分析類型，單元類型和實常數
ANTYPE, STATIC	! 定義分析類型為靜力分析
ET, 1, BEAM44	! 定義第1類單元為三維梁單元BEAM4
*DO, I, 1, 10	! 循環10次定義不同位置的單元不同的截面幾何參數
  R, I, 1, 1, (I-1)*76.2*12.7*12.7*12.7/120, 6.35, 1	! AREA1, IZ1, IY1, TKZB1, TKYB1, IX1
  RMORE, 1, 1, I*76.2*12.7*12.7*12.7/120, 6.35, 1	! AREA2, IZ2, IY2, TKZB2, TKYB2, IX2
*ENDDO
RMOD, 1, 3, 76.2*12.7*12.7*12.7/120/2	! 修改自由端的實常數為大于零的抗彎慣性矩
RLIST	! 列出實常數
MP, EX, 1, 207E3	! 定義第1類材料的彈性模量為207E3 N/mm^2
MP, GXY, 1, 207E3/2.6	! 定義第1類材料的剪切彈性模量為207E3/2.6 N/mm^2
! （9）定義節點位置
N, 1	! 在坐標原點(0, 0, 0)定義第1號節點
N, 11, 508	! 在坐標(508mm, 0, 0)定義第11號節點
FILL	! 在第1號節點和第11號節點中間填充生成2到10號節點
N, 12, , , 10	! 在坐標(0, 0, 10mm)定義第12號節點, 以定義梁的轉軸
NGEN, 10, 1, 12	! 在12號節點位置處，循環10次，定義從12到21號節點
! (10) 定義梁單元
E, 1, 2, 12	! 過節點1, 2和12, 定義單元. 12號節點是主軸方向
EGEN, 10, 1, 1, , , , , 1	! 每次實常數增加1, 生成各個單元
! (11) 施加位移約束條件
D, 11, ALL	! 約束11號節點的所有自由度
D, 1, UY, , , 10, , ROTX, ROTZ	! 約束從1號到10號節點的UY, ROTX和ROTZ自由度
F, 1, FZ, -44.483	! 在1號節點施加沿著-Z方向的44.483N
OUTPR, ALL, ALL	! 輸出所有荷載步的所有計算結果
FINISH	! 退出前處理模塊PREP7
! (12) 在求解模塊中求解
/SOLU	! 加載求解模塊
SOLVE	! 開始靜力分析求解
FINISH	! 退出求解模塊SOLUTION
! (13) 在后處理模塊中，繪制變形圖，定義單元表，提取應力計算結果
/POST1	! 進入一般后處理模塊
/VIEW, 1, 1, 1, 1	! 切換到等軸測視點
PLDISP, 2	! 繪制變形圖
ETABLE, STRS, NMISC, 1	! 以BEAM44的最大拉應力(NMISC, 1)為內容定義單元表
ESORT, STRS 	! 按最大拉應力S1對單元排序
*GET, SMAX2, SORT, , MAX	! 提取排序后的最大拉應力, 存入變量SMAX2
PRNSOL, DOF 	! 顯示節點位移計算結果
LFT_NODE = NODE (0, 0, 0)	! 選擇距離(0, 0, 0)最近的節點, 命名為LFT_NODE
*GET, DEFL2, NODE, LFT_NODE, U, Z	! 提取節點LFT_NODE的位移UZ, 存入變量DEFL2
SAVE, TABLE_2	! 保存數據到數據庫TABLE_2
! (14) 從數據庫中恢復數據，并顯示前后兩次的計算結果
RESUME, TABLE_1	! 從數據庫TABLE1中恢復數據
*STATUS, PARM	! 顯示變量中第1種情況的結果
RESUME, TABLE_2	! 從數據庫TABLE2中恢復數據
*STATUS, PARM	! 顯示變量中第2種情況的結果
FINISH	! 退出后處理模塊