! 6.6  矩形截面懸壁梁的側向穩定性分析
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!
FINISH
/CLEAR
/FILNAME, EX6.6
/TITLE, EX6.6, Lateral Torsional Buckling Analysis of a Cantilever Beam.
! （1）設置圖形選項并定義幾何模型
/PREP7
/GRA, POWER	! 打開PowerGraphics選項
K, 1, 0, 0, 0,	! 定義關鍵點1(0, 0, 0)
K, 2, 1000, 0, 0,	! 定義關鍵點2(1000mm, 0, 0)
K, 3, 500, 100, 0,	! 定義關鍵點3(500mm, 100mm, 0)
LSTR, 1, 2	! 過關鍵點1和2定義直線
! （2）定義單元類型和橫截面信息
ET, 1, BEAM189	! 定義第1類單元為空間有限應變梁Beam189
SECTYPE, 1, BEAM, RECT,	! 定義第1類梁為矩形截面梁
SECDATA, 10, 50 	! 指定梁的寬度10mm，高度50mm
! （3）定義材料特性并定位結點
MP, EX, 1, 210E3	! 定義第1類材料的彈性模量EX=210 GPa
MP, NUXY, 1, 0.0	! 定義第1類材料的泊松比NUXY=0.0
LSEL, S, , , 1, 1, 1	! 選擇1號線
LATT, 1, , 1, 0, 3, , 1	! 將1號材料，默認單元類型，3號關鍵點作為主慣性軸和默認截面類型等屬性賦給1號線。
! （4）對線剖分網格并確認梁的定位
LESIZE, ALL, , , 10	! 指定所有線剖分單元時分為10份
SECNum, 1	! 指向第1類截面類型
LMESH, ALL	! 將所有線剖分為線單元
/VIEW, , 1, 1, 1	! 改變視點為等軸測位置
/ESHAPE, 1	! 用實體模式顯示單元
EPLOT	! 顯示單元
! （5）定義邊界條件和荷載
DK, 1, ALL	! 約束1號關鍵點的所有自由度
FK, 2, FY, 1	! 在2號關鍵點施加y方向的單位集中力1N
FINISH	! 退出前處理模塊
! （6）打開預應力選項，開始靜力求解
/SOLU  	! 進入求解模塊
PSTRES, ON	! 打開預應力開關
EQSLV, SPARSE	! 選擇稀疏矩陣求解器
SOLVE  	! 開始靜力分析求解
FINISH 	! 退出求解模塊
! （7）重新進行線性屈曲求解
/SOLU	! 進入求解模塊
ANTYPE, BUCKLE	! 指定求解類型為屈曲分析
BUCOPT, LANB, 3	! 指定為分塊Lanczos方法求解前3階屈曲構型
MXPAND, 3, , , YES	! 展開前3階模態，并計算單元結果和反力
SOLVE	! 開始求解
FINISH	! 退出求解模塊
! （8）在后處理模塊顯示分析結果
/POST1	! 進入一般后處理模塊
SET, LIST	! 列表顯示結果匯總信息
SET, 1, 1	! 讀入第一個荷載步的結果
PLDisp, 2	! 顯示變形
FINISH	! 退出后處理模塊