! 5.6  懸臂長板的大撓度彎曲
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除數據，不加載初始化文件
! (1)設置工程參數、分析類型、單元類型、實常數和材料參數
/FILNAM, EX5.6(1)	! 定義工程文件名稱
/PREP7	! 加載前處理模塊
/TITLE, EX5.6(1), LARGE DEFLECTION OF A CANTILEVERED PLATE.
/NOPR	! 關閉提示信息
SMRT, OFF	! 關閉智能化功能
ANTYPE, STATIC	! 定義分析類型為靜力分析
NLGEOM, ON	! 打開大變形選項
ET, 1, SHELL43	! 定義第1類單元為4節點塑性大應變殼體單元SHELL43
R, 1, 1	! 定義板厚1mm
MP, EX, 1, 1800	! 定義第1類材料的彈性模量EX=1800 N/mm^2
MP, NUXY, , 0	! 定義第1類材料的泊松比NUXY=0
! (2)定義關鍵點、線和面, 并剖分為面單元
K, 1	! 在坐標原點定義第1個關鍵點
K, 2, 12	! 在坐標(12mm, 0)定義第2個關鍵點
K, 3, 12, 1	! 在坐標(12mm, 1mm)定義第3個關鍵點
K, 4, , 1	! 在坐標(0, 1mm)定義第4個關鍵點
L, 1, 2	! 過關鍵點1, 2定義直線
L, 3, 4	! 過關鍵點3, 4定義直線
LESIZE, ALL, , , 2	! 定義所有線剖分單元時剖分為2段
ESIZE, , 1	! 沒有指定剖分段數的線剖分時剖分為1段
A, 1, 2, 3, 4	! 過關鍵點1, 2, 3, 4定義面
/VIEW, 1, 0.61, -0.64, 0.47	! 設置觀察模型的視點
/VUP, 1, Z	! 設置Z軸向上
AMESH, 1	! 對1號面生成面單元網格
! (3)施加位移約束和荷載力
NSEL, S, LOC, X	! 選擇X=0位置的節點
D, ALL, ALL	! 固定所有節點的所有自由度
NSEL, S, LOC, X, 12	! 選擇X=12mm位置的節點
CP, 1, ROTY, ALL	! 定義自由端節點的繞Y軸的轉動自由度為耦合自由度
TORQ=15.708	! 定義集中力偶荷載15.708 N?mm
F, 2, MY, TORQ	! 在2號節點的繞Y軸方向施加力偶15.708 N?mm
NSEL, ALL	! 選擇所有結點
FINISH	! 退出前處理模塊PREP7
! (4)進入求解選模塊, 設置求解選項后, 開始求解
/SOLU	! 進入求解模塊
AUTOTS, ON	! 使用自動荷載步長
NSUBST, 20	! 在每一個荷載步最多使用10個荷載子步
LNSRCH, ON	! 使用直線搜索方法
OUTPR, BASIC, LAST	! 輸出最后一個荷載子步的基本數據結果
OUTRES, ALL, ALL	! 在求解輸出文件中保存每一個荷載子步的結果
SOLVE	! 開始求解
FINISH	! 退出求解模塊SOLU
! (5)在后處理模塊, 顯示應力和位移計算結果
/POST1	! 進入一般后處理模塊
RSYS, SOLU	! 選擇剛生成的坐標系
SET, 1	! 讀入荷載步1的計算結果
SHELL, TOP	! 選擇應力輸出為殼體上表面
PRNSOL, S, COMP	! 列表顯示節點應力結果
PRNSOL, DOF	! 列表顯示節點位移結果
! (6)提取結果后存入變量
*GET, UX_N4, NODE, 4, U, X	! 提取4號節點X方向的位移UX, 存入變量UX_N4
*GET, UZ_N4, NODE, 4, U, Z	! 提取4號節點Z方向的位移UZ, 存入變量UZ_N4
*GET, ROTY_N4, NODE, 4, ROT, Y	! 提取4號節點繞Y軸的轉角位移ROTY, 存入變量ROTY_N4
*GET, STRSS_N1, NODE, 1, S, X	! 提取1號節點沿著X軸方向的應力SX, 存入變量STRSS_N1
! (7)保存計算結果
SAVE, TABLE_1	! 保存數據到數據庫TABLE_1
FINISH	! 退出后處理模塊
! (8)重用改用有限應變層合殼單元SHELL181計算
/CLEAR, NOSTART	! 清除所有數據, 不加載初始化文件
/FILNAM, EX5.6(2)	! 定義工程文件名稱
/PREP7	! 加載前處理模塊
/TITLE, EX5.6(2), LARGE DEFLECTION OF A CANTILEVERED PLATE(SHELL181).
SMRT, OFF	! 關閉智能化功能
/NOPR	! 關閉提示信息
ANTYPE, STATIC	! 定義分析類型為靜力分析
NLGEOM, ON	! 打開大變形選項
ET, 1, SHELL181	! 定義第1類單元為有限應變層合殼單元SHELL181
R, 1, 1	! 定義單元厚度為1mm
RMORE, , , , , 10.0	! 定義薄膜單元沙漏參數10.0
MP, EX, 1, 1800	! 定義第1類材料的彈性模量EX=1800 N/mm^2
MP, NUXY, , 0	! 定義第1類材料的泊松比NUXY=0
! (9)定義關鍵點、線和面, 并將面剖分為面單元
K, 1	! 在坐標原點定義第1個關鍵點
K, 2, 12	! 在坐標(12mm, 0)定義第2個關鍵點
K, 3, 12, 1	! 在坐標(12mm, 1mm)定義第3個關鍵點
K, 4, , 1	! 在坐標(0, 1mm)定義第4個關鍵點
L, 1, 2	! 過關鍵點1, 2定義直線
L, 3, 4	! 過關鍵點3, 4定義直線
LESIZE, ALL, , , 2	! 定義所有線剖分單元時截為2段
ESIZE, , 1	! 沒有指定剖分段數的線剖分時截為1段
A, 1, 2, 3, 4	! 過關鍵點1, 2, 3, 4定義面
AMESH, 1	! 對1號面生成面單元網格
! (10) 施加位移約束和荷載力
NSEL, S, LOC, X	! 選擇X=0位置的節點
D, ALL, ALL	! 固定所有節點的所有自由度
NSEL, S, LOC, X, 12	! 選擇X=12mm位置的節點
CP, 1, ROTY, ALL	! 定義自由端節點的繞Y軸的轉動自由度為耦合自由度
TORQ=7.854	! 定義一半荷載7.854 N?mm
F, 2, MY, TORQ	! 在2號節點的繞Y軸方向施加力偶7.854 N?mm
NSEL, ALL	! 選擇所有結點
FINISH	! 退出前處理模塊PREP7
! (11) 進入求解模塊, 設置求解選項后, 開始求解
/SOLU	! 進入求解模塊
AUTOTS, ON	! 使用自動荷載步長
NSUBST, 10	! 在每一個荷載步最多使用10個荷載子步
LNSRCH, ON	! 使用直線搜索方法
OUTPR, BASIC, LAST	! 輸出最后一個荷載子步的基本數據結果
OUTRES, ALL, ALL	! 在求解輸出文件中保存每一個荷載子步的結果
SOLVE	! 開始求解
FINISH	! 退出求解模塊SOLUTION
! (12)在求解模塊重啟動分析
/SOLU	! 加載求解模塊
ANTYPE, , REST	! 重啟動分析
F, 2, MY, TORQ*2	! 給2號節點施加全部荷載2*7.854 N?mm
SOLVE	! 開始求解
FINISH	! 退出求解模塊SOLU
! (13)在后處理模塊, 顯示計算結果
/POST1	! 進入一般后處理模塊
RSYS, SOLU	! 選擇剛生成的坐標系
SET, 2 	! 讀入荷載步2(從重啟動分析)的計算結果
SHELL, TOP	! 選擇應力輸出的殼體上表面
PRNSOL, S, COMP	! 列表顯示節點應力結果
PRNSOL, DOF	! 列表顯示節點位移結果
*GET, UX_N4, NODE, 4, U, X	! 提取4號節點X方向的位移UX, 存入變量UX_N4
*GET, UZ_N4, NODE, 4, U, Z	! 提取4號節點Z方向的位移UZ, 存入變量UZ_N4
*GET, ROTY_N4, NODE, 4, ROT, Y	! 提取4號節點繞Y軸的轉角位移ROTY, 存入變量ROTY_N4
*GET, STRSS_N1, NODE, 1, S, X	! 提取1號節點沿著X軸方向的應力SX, 存入變量STRSS_N1
! (14) 顯示計算結果
*STATUS, PARM	! 顯示保存在變量中的結果
Resume, TABLE_1	! 從數據庫TABLE_1恢復數據
*STATUS, PARM	! 顯示保存在變量中的結果
Finish