! 5.1  含橢圓孔的橢圓薄膜在外部張力作用下的靜力分析
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除數據, 不加載初始化文件
! （1）設置工程選項, 分析類型、單元類型和材料參數
/FILNAME, EX5.1(1)	! 設置工程文件名稱
/PREP7	! 加載前處理模塊
SMRT, OFF	! 關閉智能化網格功能
/TITLE, EX5.1(1), ELLIPTIC MEMBRANE UNDER A UNIFORMLY-DISTRIBUTED LOAD(PLANE42).
/COM,   USING 2-D STRUCTURAL SOLID, PLANE42
ANTYPE, STATIC	! 定義分析類型為靜力分析
ET, 1, PLANE42, , , 3	! 定義第1類單元為帶厚度的平面應力單元PLANE42
MP, EX, 1, 210E9	! 定義第1類材料的彈性模量EX＝210GPa
MP, NUXY, 1, 0.3	! 定義第1類材料的泊松比NUXY=0.3
R, 1, 0.1  	! 設置厚度0.1m
! （2）在橢圓坐標系下, 定義關鍵點、線和面
! LOCAL, KCN, KCS, XC, YC, ZC, THXY, THYZ, THZX, PAR1, PAR2 
LOCAL, 11, 1, , , , , , , 0.5	! 在坐標原點定義橢圓坐標系, Y方向的短軸與X方向的長軸之比為0.5
K, 1, 2, 90	! 在坐標(2m, 90度)處定義第1個關鍵點
K, 2, 2, 0	! 在坐標(2m, 0)處定義關鍵點2
L, 1, 2	! 過關鍵點1, 2定義直線
LOCAL, 12, 1, , , , , , , 0.8461585	! 在坐標原點定義橢圓坐標系, Y方向的短軸與X方向的長軸之比為0.8461585
K, 3, 3.25, 90	! 在坐標(3.25m, 90度)處定義關鍵點3
K, 4, 3.25, 0	! 在坐標(3.25m, 0)處定義關鍵點4
L, 3, 4	! 過關鍵點3, 4定義直線
CSYS, 0	! 將當前坐標系切換到整體直角坐標系
L, 2, 4	! 過關鍵點2, 4定義直線
L, 1, 3	! 過關鍵點1, 3定義直線
AL, 2, 4, 1, 3	! 用線2, 4, 1和3定義面
! (3)在相應的線上施加位移約束
DL, 4, 1, SYMM	! 對1號面的4號線施加對稱位移邊界條件
DL, 3, 1, SYMM	! 對1號面的3號線施加對稱位移邊界條件
SFL, 2, PRES, -10E6	! 對2號線施加線壓力荷載10E6 N/m
MSHK, 2   	! 如果可能, 采用映射網格剖分
MSHA, 0, 2D  	! 使用四邊形網格
SAVE	! 保存模型數據到默認數據庫
FINISH	! 退出前處理模塊
! (4)執(zhí)行自適應網格剖分和求解宏程序, 執(zhí)行自定義的后處理宏程序
ADAPT, 3, 7, , , 1	! 使用ANSYS預定義的自適應網格剖分宏程序, 用NSOLN=3, STARGT=7, 和FACMX=1等參數來求解
/Post1	! 進入一般后處理模塊
EPLOT	! 繪制單元
PLDISP, 1	! 顯示變形圖
PLNSOL, S, X	! 顯示水平應力分布
PLNSOL, S, Y	! 顯示豎向應力分布
PLNSOL, S, 1	! 顯示最大拉應力分布
FINISH	! 退出后處理模塊
! (5)定義完成后處理功能的宏程序MAC
*CREATE, MAC	! 定義用于完成后處理的宏程序MAC
/Graphics, Full	! 關閉PowerGraphics, 以顯示能量范數誤差
/POST1    	! 進入一般后處理模塊
PRERR      	! 打印能量范數誤差(SEPC)
NSEL, S, LOC, Y, 0.0	! 選擇Y=0位置的節(jié)點
NSEL, R, LOC, X, 2.0	! 繼續(xù)選擇X=2.0m位置的節(jié)點
*GET, MNODE, NODE, , NUM, MAX	! 提取當前選擇集中最大節(jié)點號, 保存到變量MNODE
*GET, SY_D, NODE, MNODE, S, Y	! 提取節(jié)點MNODE的豎向應力SY, 存入變量SY_D
NSEL, ALL	! 選擇所有節(jié)點
ESEL, ALL	! 選擇所有單元
*STATUS, PARM	! 顯示內存變量的值
/Graphics, Power	! 打開PowerGraphics, 以顯示其他結果
*END             	! 宏程序MAC定義結束
! (6)用前定義的宏程序完成后處理, 將計算結果保存到變量
*USE, MAC   	! 執(zhí)行后處理宏程序MAC
SAVE, TABLE1	! 保存數據到默認數據庫TABLE1
FINISH	! 退出后處理模塊
! (7)改變單元為PLANE82, 重新計算
/CLEAR, NOSTART	! 在開始第二部分求解前, 清除數據, 不加載初始化文件! 
/FILNAME, EX5.1(2)	! 設置工程文件名稱
/PREP7	! 加載前處理模塊
SMRT, OFF	! 關閉智能化網格功能
RESUME, EX5.1(1).db	! 重新讀入模型數據
/TITLE, EX5.1(2), ELLIPTIC MEMBRANE UNDER A UNIFORMLY-DISTRIBUTED LOAD(PLANE82).
ET, 1, PLANE82, , , 3	! 定義第1類單元為帶厚度的二維8節(jié)點平面應力單元PLANE82
FINISH	! 退出前處理模塊PREP7
! （8）執(zhí)行自適應網格宏程序, 執(zhí)行自定義的后處理宏程序
ADAPT, 2, 5, , , 1	! 使用ANSYS預定義的自適應網格剖分宏程序, 用NSOLN=2, STARGT=5, 和FACMX=1等參數來求解
/Post1	! 進入一般后處理模塊
/WINDOW, 5, OFF	! 關閉5號窗口的顯示
/WINDOW, 1, FULL	! 將1號窗口的顯示設置為整個圖形窗口
PLNSOL, S, EQV	! 顯示等效MISES應力分布
PLDISP, 1	! 顯示變形圖
PLNSOL, S, X	! 顯示水平應力分布
PLNSOL, S, Y	! 顯示豎向應力分布
PLNSOL, S, 1	! 顯示最大拉應力分布
FINISH	! 退出后處理模塊
! （9）用前定義的宏程序完成后處理, 將計算結果保存到變量
*USE, MAC 	! 執(zhí)行后處理宏程序MAC
FINISH	! 退出后處理模塊