! 5.7  用殼體單元分析受均布荷載作用的固支圓板大撓度彎曲
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除數據，不加載初始化文件
! （1）設置工程參數、分析類型、單元類型、實常數和材料參數
/FILNAME, EX5.7	! 定義工程文件名稱
/PREP7	! 加載前處理模塊
/TITLE, EX5.7, LARGE DEFLECTION BENDING OF A CIRCULAR PLATE.
ANTYPE, STATIC	! 定義分析類型為靜力分析
NLGEOM, ON	! 打開大變形選項
SSTIF, ON	! 打開應力強化選項
ET, 1, SHELL51, , , , , , , 1	! 定義第1類單元為軸對稱殼體單元SHELL51，抑制單元輸出
R, 1, 0.0025	! 定義殼體厚度為0.0025m
MP, EX, 1, 2E11	! 定義第1類材料的彈性模量EX=200GPa
MP, NUXY, 1, 0.3	! 定義第1類材料的泊松比NUXY=0.3
! (2)定義節點和單元
N, 1	! 在坐標原點(0, 0)定義第1號結點
N, 11, 0.25	! 在坐標(0.25m, 0)定義第11號結點
FILL	! 在1號節點和11號節點之間填充生成其余節點
E, 1, 2	! 過節點1和節點2定義單元
EGEN, 10, 1, 1	! 按照前面單元模式, 循環10次生成其余單元, 每次節點號增加1
! (3)設置求解選項, 位移約束和荷載
OUTPR, , LAST	! 只輸出最后一個荷載步的基本結果
NSEL, S, LOC, X, 0	! 選擇X=0位置的節點
DSYM, SYMM, X	! 施加關于X軸(YZ平面)的對稱位移條件
NSEL, ALL	! 選擇所有結點
D, 11, ALL	! 完全固定11號節點
SFE, 1, 1, PRES, , 6585.2	! 給1號單元施加單元表面荷載6585.2 N/m2
*REPEAT, 10, 1	! 重復循環10次, 給其余單元施加表面荷載
FINISH	! 退出前處理模塊PREP7
! (4)進入求解模塊后, 開始求解
/SOLU	! 加載求解模塊
SOLVE	! 開始非線性迭代求解
FINISH	! 退出求解模塊SOLUTION
! (5)在后處理模塊中, 顯示計算結果
/POST1	! 進入一般后處理模塊
*GET, UY, NODE, 1, U, Y	! 提取1號節點的Y方向位移UY, 存入變量UY
*STATUS, PARM	! 顯示內存變量的值