! 4.13  均勻拉力作用下含圓孔板的孔邊應力集中
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!
FINISH	! 退出以前模塊
/CLEAR, NOSTART	! 清除數據, 不加載初始化文件
! （1）設置工程選項、分析類型、單元類型和材料參數
/FILNAME, EX4.13(1)	! 定義工程文件名稱
/PREP7	! 加載前處理模塊
SMRT, OFF	! 關閉智能化功能
/TITLE, EX4.13(1), STRESS CONCENTRATION AT A HOLE IN A PLATE with PLANE2.
/NOPR	! 關閉輸入信息反饋提示
ANTYPE, STATIC	! 定義分析類型為靜力分析
ET, 1, PLANE2	! 定義第1類單元為二維6節點三角形固體單元PLANE2
MP, EX, 1, 207E3	! 定義第1類材料的彈性模量EX ＝ 207 GPa
MP, NUXY, 1, 0.3	! 定義第1類材料的泊松比NUXY ＝ 0.3
! （2）定義關鍵點
K, 1, 152.4	! 在坐標(152.4mm, 0)定義第1個關鍵點
K, 2, 152.4, 152.4	! 在坐標(152.4mm, 152.4mm)定義第2個關鍵點
K, 3, , 152.4	! 在坐標(0, 152.4mm)定義第3個關鍵點
K, 4, , 12.7	! 在坐標(0, 12.7mm)定義第4個關鍵點
K, 5, 12.7	! 在坐標(12.7mm, 0)定義第5個關鍵點
K, 6	! 在坐標(0, 0)定義第6個關鍵點
! （3）定義線和線單元剖分的尺寸
L, 1, 2	! 過關鍵點1, 2定義直線
L, 2, 3	! 過關鍵點2, 3定義直線
L, 3, 4	! 過關鍵點3, 4定義直線
LESIZE, 3, , , 4, .25	! 給3號線定義剖分段數為4和間距比例0.25
LARC, 4, 5, 6, 12.7	! 過4號點和5號點, 以6號點為圓心, 以12.7mm為半徑定義圓弧
LESIZE, 4, , , 6	! 給4號線定義剖分段數為6
L, 5, 1	! 過關鍵點5, 1定義直線
LESIZE, 5, , , 4, 4	! 給5號線定義剖分段數為4和間距比例為4
! （4）定義面并將其剖分為面單元, 圖形顯示單元布置
AL, 1, 2, 3, 4, 5	! 用線1, 2, 3, 4和5定義面
ESIZE, , 4 	! 單元剖分尺寸為分為4段
AMESH, ALL	! 將所有的面剖分為面單元
/AUTO, 1	! 設置模型顯示的最佳比例
/PLOPTS, INFO, 0	! 關閉圖形窗口中圖例信息顯示欄
/PLOPTS, WINS, 0	! 關閉圖形窗口的自適應調整
/WINDOW, , LTOP	! 在圖形窗口的左上角定義窗口1
EPLOT	! 圖形顯示單元布置
! （5）選擇相應的線, 施加位移和荷載邊界條件
LSEL, S, LINE, , 3, 5, 2	! 選擇3號線和5號線
DL, ALL, , SYMM	! 對所有的選擇集中的線施加對稱位移約束
LSEL, S, LINE, , 1	! 選擇1號線
NSLL, , 1	! 選擇1號線上的節點
SF, ALL, PRES, -6.895	! 施加板端的拉力荷載6.895 MPa
LSEL, ALL	! 選擇所有的線
NSEL, ALL	! 選擇所有節點
CSYS, 1	! 將坐標系切換到柱坐標系
FINISH	! 退出前處理模塊PREP7
! （6）進入求解模塊求解
/SOLU	! 加載求解模塊
SOLVE	! 開始求解
FINISH	! 退出求解模塊SOLUTION
SAVE	! 保存數據到數據庫EX4.13(1).DB
! （7）進入后處理, 顯示計算結果
/POST1	! 進入一般后處理模塊
SET, 1, 1	! 讀入第1個荷載步的計算結果
NSORT, S, X, , , 3	! 將計算結果按SX排序, 只保留前3項最高的結果
PRNSOL, S, COMP	! 列表顯示應力分量計算結果
/WINDOW, 1, OFF	! 關閉1號窗口的顯示
/NOERASE	! 不擦除圖形窗口
/DSCALE, 2, 1	! 用縮放位移約束命令, 縮放位移約束實部2倍, 虛部1倍
/WINDOW, 2, RTOP	! 在圖形窗口的右上角定義窗口2
PLNSOL, S, X	! 圖形顯示水平方向的節點應力SX分布圖
PLNSOL, S, Y	! 圖形顯示豎直方向的節點應力SY分布圖
*GET, CRSESTR, NODE, 18, S, X	! 提取18號節點的應力SX, 存入變量CRSESTR
*STATUS, PARM	! 顯示內存變量中的值
! （8）將計算結果保存到文件TABLE_1
SAVE, TABLE_1	! 保存數據到數據庫TABLE_1
! （9）改變單元類型為PLANE42, 用子模型技術重新計算
FINISH	! 退出后處理模塊
/CLEAR, NOSTART	! 清除數據, 不加載初始化文件! 
/FILNAME, 4.13(2)	! 定義子結構分析文件名
/PREP7	! 加載前處理模塊
SMRT, OFF	! 關閉智能化網格功能
/NOPR	! 關閉輸入信息反饋提示
/TITLE, EX4.13(2), STRESS CONCENTRATION AT A HOLE IN A PLATE With PLANE42 SUB-MODEL.
ANTYPE, STATIC	! 定義分析類型為靜力分析
ET, 1, PLANE42	! 定義第1類單元為二維等參固體單元PLANE42
MP, EX, 1, 207E3	! 定義第1類材料的彈性模量EX＝207 GPa
MP, NUXY, 1, 0.3	! 定義第1類材料的泊松比NUXY＝0.3
! (10) 在柱坐標系下定義孔邊子模型上的關鍵點
CSYS, 1	! 切換到柱坐標系下
K, 10, 12.7, 45	! 為子模型在坐標（12.7mm, 45度）定義第10個關鍵點
K, 11, 12.7, 90	! 在坐標（12.7mm, 90度）定義第11個關鍵點
K, 12, 38.1, 45	! 在坐標（38.1mm, 45度）定義第12個關鍵點
K, 13, 38.1, 90	! 在坐標（38.1mm, 90度）定義第13個關鍵點
A, 10, 12, 13, 11	! 關鍵點10, 12, 13和11定義面
ESIZE, , 8	! 單元剖分為8格
MSHK, 1 	! 設定映射網格選項
MSHA, 0, 2D	! 使用四邊形網格
AMESH, 1	! 按照上面選項, 對1號面積進行網格剖分
! (11) 設置窗口選項后, 繪制單元布置
/WINDOW, 2, OFF	! 關閉2號窗口的顯示
/ERASE 	! 擦除屏幕
/PLOPTS, INFO, 0	! 關閉圖形窗口中圖例信息顯示欄
/PLOPTS, WINS, 0	! 關閉圖形窗口的自適應調整
/WINDOW, , LTOP	! 在圖形窗口的左上角定義窗口1
/USER   	! 使用用戶自定義縮放比例
/DIST, 1, 18	! 用戶自定義1號窗口的觀察距離18mm
/FOCUS, 1, 14, 25	! 用戶自定義1號窗口的焦點(14mm, 25mm)
EPLOT	! 圖形顯示單元布置
LSEL, S, LINE, , 1, 2	! 選擇1號線和2號線
NSLL, , 1 	! 選擇1號線上的節點
NWRITE 	! 將節點數據寫入到文件 4.13(2).NODB
LSEL, ALL	! 選擇所有的線
NSEL, ALL	! 選擇所有節點
FINISH	! 退出前處理模塊PREP7
SAVE	! 保存子模型數據到數據庫4.13(2).DB
! (12) 在后處理模塊恢復模型數據, 讀入計算結果
/POST1	! 進入一般后處理模塊
RESUME, EX4.13(1).db	! 從數據庫EX4.13(1).DB中恢復平板模型數據
FILE, EX4.13(1), rst	! 定義結果文件
CBDOF, , , , , , , 0, , 0	! 激活剖切邊界插值
! (13) 重新讀模型, 施加不同的邊界條件
FINISH	! 退出一般后處理模塊
/PREP7	! 加載前處理模塊
SMRT, OFF	! 關閉智能化功能
RESUME, 4.13(2).DB	! 從子模型文件4.13(2).DB中恢復子模型數據
/NOPR	! 關閉提示反饋
/INPUT, , cbdo, , :cb1	! 從文件4.13(2).CBDO中讀入插值邊界信息
/GOPR	! 打開提示反饋
LSEL, S, LINE, , 3	! 選擇3號線, 以應用其余的邊界條件
DL, ALL, , SYMM	! 對選擇集中的線應用對稱邊界條件
FINISH	! 退出前處理模塊PREP7
! (14) 進入求解模塊求解
/SOLU	! 加載求解模塊
SOLVE	! 開始求解
FINISH	! 退出求解模塊SOLU
! (15) 在后處理模塊顯示應力計算結果, 
/POST1	! 進入一般后處理模塊
SET, 1, 1	! 讀入第1個荷載步的計算結果
NSORT, S, X, , , 3	! 將計算結果按SX排序, 只保留前3項最高的結果
PRNSOL, S, COMP	! 列表顯示應力分量顯示
/WINDOW, 1, OFF	! 關閉1號窗口的顯示
/AUTO, 3	! 指定3號窗口為以最佳比例顯示圖形
/WINDOW, 3, BOT	! 在圖形窗口的下半部分定義窗口3
/NOERASE	! 不擦除圖形窗口
/PLOPTS, MINM, 1	! 顯示最大和最小圖例
/USER, 3	! 指定用戶自定義的視點距離和焦點位置
/DIST, 3, 5.08	! 指定用戶自定義的視點距離為5.08mm
/FOCUS, 3, 5.08, 12.7	! 指定用戶自定義的焦點位置(5.08mm, 12.7mm)
/CONTOUR, 3, , AUTO	! 給窗口3應用自動梯度比例
PLNSOL, SX	! 列表顯示應力分量SX計算結果
*GET, SUBSTR, NODE, 18, S, X	! 提取18號節點的X方向應力SX, 存入變量SUBSTR
! (16) 將計算結果保存到文件TABLE_2
SAVE, TABLE_2	! 保存數據到默認數據庫TABLE_2
! (17) 改變單元類型為PLANE146, 重新計算
FINISH	! 退出后處理模塊
/CLEAR, NOSTART	! 清除所有數據, 不加載初始化文件
/FILNAME, 4.13(3)	! 定義子結構分析文件名
/PREP7	! 加載前處理模塊
SMRT, OFF	! 關閉智能化網格功能
/TITLE, EX4.13(3), STRESS CONCENTRATION AT A HOLE IN A PLATE With PLANE146 Element.
/NOPR	! 關閉輸入信息反饋提示
ANTYPE, STATIC	! 定義分析類型為靜力分析
ET, 1, PLANE146	! 定義第1類單元為二維結構固體自適應單元PLANE146
MP, EX, 1, 207E3	! 定義第1類材料的彈性模量EX ＝ 207 GPa
MP, NUXY, 1, 0.3	! 定義第1類材料的泊松比NUXY ＝ 0.3
! (18)定義關鍵點
K, 1, 152.4	! 在坐標(152.4mm, 0)定義第1個關鍵點
K, 2, 152.4, 152.4	! 在坐標(152.4mm, 152.4mm)定義第2個關鍵點
K, 3, , 152.4	! 在坐標(0, 152.4mm)定義第3個關鍵點
K, 4, , 12.7	! 在坐標(0, 12.7mm)定義第4個關鍵點
K, 5, 12.7	! 在坐標(12.7mm, 0)定義第5個關鍵點
K, 6	! 在坐標(0, 0)定義第6個關鍵點
! (19)定義線和單元剖分的段數
L, 1, 2	! 過關鍵點1, 2定義直線
L, 2, 3	! 過關鍵點2, 3定義直線
L, 3, 4	! 過關鍵點3, 4定義直線
LESIZE, 3, , , 4, .25	! 給3號線定義剖分段數為4和間距比例0.25
LARC, 4, 5, 6, 12.7	! 以6號點為圓心, 過5號點和6號點, 以12.7mm為半徑定義圓弧
LESIZE, 4, , , 6	! 給4號線定義剖分段數為6
L, 5, 1	! 過關鍵點5, 1定義直線
LESIZE, 5, , , 4, 4	! 給3號線定義剖分段數為4和間距比例為4
! (20) 定義面并將其剖分為面單元, 圖形顯示單元布置
AL, 1, 2, 3, 4, 5	! 用線1, 2, 3, 4和5定義面
ESIZE, , 4  	! 單元剖分尺寸為分為4段
AMESH, ALL 	! 將所有的面剖分為面單元
/AUTO, 1	! 設置模型顯示的最佳比例
/PLOPTS, INFO, 0	! 關閉圖形窗口中圖例信息顯示欄
/PLOPTS, WINS, 0	! 關閉圖形窗口的自適應調整
/WINDOW, , LTOP	! 在圖形窗口的左上角定義窗口1
EPLOT	! 圖形顯示單元布置
! (21) 選擇相應的線, 施加位移和荷載邊界條件
LSEL, S, LINE, , 3, 5, 2	! 選擇3號線和5號線
DL, ALL, , SYMM	! 對所有的選擇集中的線施加對稱位移約束
LSEL, S, LINE, , 1	! 選擇1號線
NSLL, , 1	! 選擇1號線上的節點
SF, ALL, PRES, -6.895	! 施加板端的拉力荷載6.895 MPa
LSEL, ALL	! 選擇所有的線
NSEL, ALL	! 選擇所有節點
CSYS, 1	! 將坐標系切換到柱坐標系
FINISH	! 退出前處理模塊PREP7
! (22) 進入求解模塊求解
/SOLU	! 加載求解模塊
SOLVE	! 開始求解
FINISH	! 退出求解模塊SOLU
! (23) 進入后處理, 顯示計算結果
/POST1	! 進入一般后處理模塊
SET, 1, 1	! 讀入第1個荷載步的計算結果
/WINDOW, 1, OFF	! 關閉1號窗口的顯示
/NOERASE	! 不擦除圖形窗口
/DSCALE, 2, 1	! 用縮放位移約束命令, 縮放位移約束實部2倍, 虛部1倍
/WINDOW, 2, RTOP	! 在圖形窗口的右上角定義窗口2
PLNSOL, S, X	! 圖形應力SX的分布圖
*GET, SUBSTR, NODE, 18, S, X	! 提取18號節點的應力SX, 存入變量SUBSTR
! (24) 將結果數據保存到數據庫TABLE_3
SAVE, TABLE_3	! 保存數據到默認數據庫TABLE_3
FINISH	! 退出后處理模塊
! (25) 從各個數據庫中恢復數據, 依次輸出顯示
RESUME, TABLE_1
*STATUS, PARM	! 列表顯示變量參數的值
RESUME, TABLE_2
*STATUS, PARM	! 列表顯示變量參數的值
RESUME, TABLE_3
*STATUS, PARM	! 列表顯示變量參數的值
FINISH	! 退出后處理模塊
! (26) 刪除中間數據文件
/DEL, TABLE_1
/DEL, TABLE_2
/DEL, TABLE_3