! 4.12  端部集中力下懸臂圓環(huán)曲梁彎曲的三維分析
! 本程序來源于邢靜忠編著的《ANSYS應(yīng)用實(shí)例與分析》，科學(xué)出版社，2006年
!
FINISH	! 退出以前模塊
/CLEAR, NOSTART	! 清除數(shù)據(jù), 不加載初始化文件
! (1)定義工程選項(xiàng)、分析類型、單元類型和材料參數(shù)
/FILNAME, EX4.12(1)	! 定義工程文件名稱
/PREP7	! 加載前處理模塊
SMRT, OFF	! 關(guān)閉智能化網(wǎng)格功能
/TITLE, EX4.12(1), BENDING OF A CURVED BEAM WITH SOLID5.
ANTYPE, STATIC	! 定義分析類型為靜力分析
ET, 1, SOLID5, 2	! 定義第1類單元為六面體固體SOLID5
MP, EX, 1, 70E3	! 定義第1類材料的彈性模量EX=70 GPa
MP, NUXY, 1, 0.25	!  定義第1類材料的泊松比NUXY=0.25
! (2)在柱坐標(biāo)系定義關(guān)鍵點(diǎn)、線和體
CSYS, 1	! 切換當(dāng)坐標(biāo)系為柱坐標(biāo)系
K, 1, 104.6	! 在內(nèi)徑坐標(biāo)(104.6mm=4.12in, 0)定義第1個關(guān)鍵點(diǎn)
K, 2, 109.7	! 在外徑坐標(biāo)(109.7mm=4.32in, 0)定義第2個關(guān)鍵點(diǎn)
KGEN, 2, 1, 2, 1, , , 2.54	! 循環(huán)2次, 按照1號節(jié)點(diǎn)到2號節(jié)點(diǎn)增量為1的節(jié)點(diǎn)范圍, Z坐標(biāo)增量2.54mm, 生成節(jié)點(diǎn)3和4
KGEN, 2, 1, 4, 1, , 90	! 循環(huán)2次, 按照1號節(jié)點(diǎn)到4號節(jié)點(diǎn)增量為1的節(jié)點(diǎn)范圍, Y坐標(biāo)增量90度
L, 1, 5	! 過1號關(guān)鍵點(diǎn)和5號關(guān)鍵點(diǎn), 定義線
LESIZE, 1, , , 20	! 指定1號線的單元剖分段數(shù)為20
V, 1, 2, 4, 3, 5, 6, 8, 7	! 過關(guān)鍵點(diǎn)定義體
ESIZE, , 1	! 指定單元剖分段數(shù)為1
/VIEW, 1, 1, 1, 1	! 切換視點(diǎn)到等軸測方向
/Device, VECTOR, 1	! 打開圖形矢量顯示開關(guān), 更清晰地顯示網(wǎng)格模型
VMESH, 1	! 對1號體執(zhí)行體單元剖分
! (3)選擇節(jié)點(diǎn), 施加位移約束和荷載
/Device, VECTOR, 0	! 關(guān)閉圖形矢量顯示開關(guān)
NSEL, S, LOC, Y, 0	! 選擇位于Y=0位置的節(jié)點(diǎn)
D, ALL, ALL, 0	! 位移約束所有選擇集中節(jié)點(diǎn)的所有自由度
NSEL, ALL	! 選擇所有節(jié)點(diǎn)
FK, 5, FY, 1.1121	! 給5號關(guān)鍵點(diǎn)施加沿著Y方向的集中力, 大小為1.1121N
*REPEAT, 4, 1	! 循環(huán)4次執(zhí)行前面施加集中力的命令, 每次節(jié)點(diǎn)號碼增量1
NOORDER	! 按原始單元次序, 對單元排序
FINISH	! 退出前處理模塊PREP7
! (4)進(jìn)入求解模塊, 開始求解
/SOLU	! 加載求解模塊
SOLVE	! 開始求解
FINISH	! 退出求解模塊SOLU
! (5)定義完成后處理的宏程序
*CREATE, MAC	! 開始定義完成后處理的宏程序
/POST1	! 進(jìn)入一般后處理模塊
CSYS, 0	! 切換當(dāng)坐標(biāo)系為直角坐標(biāo)系
NSEL, S, LOC, X, 0, 0.01, , 1	! 選擇X=0到X=0.01范圍的節(jié)點(diǎn), 忽略坐標(biāo)正負(fù)號
PRNSOL, U, COMP	! 列表顯示結(jié)點(diǎn)位移值
LCDEF, 1, 1	! 生成第1種荷載工況, 按照第1荷載步
LCFACT, 1, (1/2.187)	! 給第1荷載工況定義荷載縮放因子(1/2.187)
LCASE, 1	! 將第1荷載工況計(jì)算結(jié)果讀入到數(shù)據(jù)庫
PRNSOL, U, COMP	! 列表顯示結(jié)點(diǎn)位移值
FINISH	! 退出后處理模塊
*END	! 宏程序定義結(jié)束
! (6)提取自由端位移計(jì)算結(jié)果
/POST1	! 進(jìn)入一般后處理模塊
*GET, U1, NODE, 5, U, Y	! 提取自由端的5號節(jié)點(diǎn)Y方向的位移計(jì)算結(jié)果UY
PLNSOL, S, Y	! 圖形顯示節(jié)點(diǎn)環(huán)向應(yīng)力計(jì)算結(jié)果
! （7）將計(jì)算結(jié)果保存到文件
SAVE, TABLE_1	! 保存數(shù)據(jù)到數(shù)據(jù)庫TABLE_1
*USE, MAC	! 執(zhí)行后處理宏程序MACRO
! (8)重新開始用SOLID92單元分析
FINISH	! 退出以前模塊
/CLEAR, NOSTART	! 清除數(shù)據(jù), 不加載初始化文件
/FILNAME, EX4.12(2)	! 定義工程文件名稱
/PREP7	! 加載前處理模塊
SMRT, OFF	! 關(guān)閉智能化網(wǎng)格功能
/TITLE, EX4.12(2), BENDING OF A CURVED BEAM WITH SOLID92.
ET, 1, SOLID92	! 定義第1類單元為四面體固體單元SOLID92
MP, EX, 1, 70E3	! 定義第1類材料的彈性模量EX=70 GPa
MP, NUXY, 1, 0.25	!  定義第1類材料的泊松比NUXY=0.25
! (9)通過關(guān)鍵點(diǎn)定義模型
CSYS, 1	! 切換當(dāng)前坐標(biāo)系為柱坐標(biāo)系
/VIEW, 1, 1, 1, 1	! 切換視點(diǎn)到等軸測方向
K, 1, 104.6	! 在內(nèi)徑坐標(biāo)(104.6mm=4.12in, 0)定義第1個關(guān)鍵點(diǎn)
K, 2, 109.7	! 在外徑坐標(biāo)(109.7mm=4.32in, 0)定義第2個關(guān)鍵點(diǎn)
KGEN, 2, 1, 2, 1, , , 2.54	! 循環(huán)2次, 按照1號節(jié)點(diǎn)到2號節(jié)點(diǎn)增量為1的節(jié)點(diǎn)范圍, Z坐標(biāo)增量2.54mm, 生成節(jié)點(diǎn)3和4
KGEN, 2, 1, 4, 1, , 4.5	! 循環(huán)2次, 按照1號節(jié)點(diǎn)到4號節(jié)點(diǎn)增量為1的節(jié)點(diǎn)范圍, Y坐標(biāo)增量4.5度
V, 1, 2, 4, 3, 5, 6, 8, 7	! 過關(guān)鍵點(diǎn)定義體
VGEN, 20, 1, 1, 1, , 4.5, , 4	! 循環(huán)20次, 按照1號體的樣式, X坐標(biāo)不變, Y坐標(biāo)增加4.5度, 每次關(guān)鍵點(diǎn)增量為4, 生成19個體
ESIZE, , 1	! 指定單元剖分段數(shù)為1
VMESH, ALL	! 將所有體剖分為體單元
! (10)選擇節(jié)點(diǎn), 施加位移約束和荷載
NSEL, S, LOC, Y, 0	! 選擇位于Y=0位置的節(jié)點(diǎn)
D, ALL, ALL, 0	! 位移約束所有選擇集中節(jié)點(diǎn)的所有自由度
NSEL, S, LOC, Y, 90	! 選擇位于Y=90度位置的節(jié)點(diǎn)
CP, 1, UY, ALL	! 將90度面上的所有節(jié)點(diǎn)的UY方向定義為耦合自由度
NSEL, R, LOC, X, 109.7	! 在選擇集中繼續(xù)選擇內(nèi)側(cè)X＝109.7mm的節(jié)點(diǎn)
NSEL, R, LOC, Z, 0	! 在選擇集中繼續(xù)選擇Z＝0的節(jié)點(diǎn)
F, ALL, FY, 4.4483	! 給選擇集中的所有節(jié)點(diǎn)施加沿著Y方向的大小等于4.4483N的集中力
NSEL, ALL	! 選擇所有節(jié)點(diǎn)
SAVE	! 保存數(shù)據(jù)到默認(rèn)數(shù)據(jù)庫
FINISH	! 退出前處理模塊PREP7
! (11)進(jìn)入求解模塊開始求解, 并使用前面的宏程序進(jìn)行相應(yīng)的后處理
/SOLU	! 加載求解模塊
SOLVE	! 開始求解
FINISH	! 退出求解模塊SOLU
! (12)選擇自由端節(jié)點(diǎn), 提取位移計(jì)算結(jié)果
/POST1	! 進(jìn)入一般后處理模塊
PLNSOL, S, Y	! 圖形顯示節(jié)點(diǎn)環(huán)向應(yīng)力計(jì)算結(jié)果
*USE, MAC	! 執(zhí)行后處理宏程序MAC
NDE=NODE(0, 109.7, 0)	! 選擇距離(0, 109.7mm, 0)位置最近的節(jié)點(diǎn), 命名為NDE
*GET, UY2, NODE, NDE, U, Y	! 提取節(jié)點(diǎn)NDE的Y方向位移UY, 并保存到變量UY2
*SET, U2, (UY2/(1/2.187))	! 給U2變量賦值UY2/(1/2.187)
! (13)將計(jì)算結(jié)果保存到文件
SAVE, TABLE_2	! 保存數(shù)據(jù)到數(shù)據(jù)庫TABLE_2
! (14)重新開始用SOLID98單元分析
/PREP7	! 加載前處理模塊
SMRT, OFF	! 關(guān)閉智能化功能
RESUME	! 重新載入模型數(shù)據(jù)
/TITLE, EX4.12(3), BENDING OF A CURVED BEAM WITH SOLID98.
ET, 1, SOLID98, 2	! 定義第1類單元為四面體耦合場單元SOLID98, 使用線位移自由度UX, UY和UZ
FINISH	! 退出前處理模塊PREP7
! (15)進(jìn)入求解模塊開始求解, 并使用前面的宏程序進(jìn)行相應(yīng)的后處理
/SOLU	! 加載求解模塊
SOLVE	! 開始求解
FINISH	! 退出求解模塊SOLU
! (16)選擇自由端節(jié)點(diǎn), 提取位移計(jì)算結(jié)果
/POST1	! 進(jìn)入一般后處理模塊
PLNSOL, S, Y	! 圖形顯示節(jié)點(diǎn)環(huán)向應(yīng)力計(jì)算結(jié)果
*USE, MAC 	! 執(zhí)行后處理宏程序MAC
NDE=NODE(0, 109.7, 0)	! 選擇距離(0, 109.7mm, 0)位置最近的節(jié)點(diǎn), 命名為NDE
*GET, UY3, NODE, NDE, U, Y	! 提取節(jié)點(diǎn)NDE的Y方向位移UY, 并保存到變量UY3
*SET, U3, (UY3/(1/2.187))	! 給U3變量賦值UY3/(1/2.187)
! (17)將計(jì)算結(jié)果保存到文件
SAVE, TABLE_3	! 保存數(shù)據(jù)到數(shù)據(jù)庫TABLE_3
! (18)重新開始用SOLID148單元分析
/PREP7	! 加載前處理模塊
SMRT, OFF	! 關(guān)閉智能化網(wǎng)格功能
RESUME	! 重新載入模型數(shù)據(jù)
/TITLE, EX4.12(4), BENDING OF A CURVED BEAM WITH SOLID148.
ET, 1, SOLID148, 2	! 定義第1類單元為三維四面體結(jié)構(gòu)固體自適應(yīng)單元SOLID148, 開始水平為2級
FINISH	! 退出前處理模塊PREP7
! (19)進(jìn)入求解模塊開始求解
/SOLU	! 加載求解模塊
SOLVE	! 開始求解
FINISH	! 退出求解模塊SOLU
! (20)圖形顯示環(huán)向應(yīng)力分布, 使用宏程序進(jìn)行后處理, 選擇自由端節(jié)點(diǎn), 提取位移結(jié)果
/POST1	! 進(jìn)入一般后處理模塊
SET, LAST	! 讀入最后一個荷載步的計(jì)算結(jié)果
PLNSOL, S, Y	! 圖形顯示節(jié)點(diǎn)環(huán)向應(yīng)力的結(jié)果
*USE, MAC   	! 執(zhí)行后處理宏程序MACRO
NDE=NODE(0, 109.7, 0)	! 選擇距離(0, 109.7mm, 0)位置最近的節(jié)點(diǎn), 命名為NDE
*GET, UY3, NODE, NDE, U, Y	! 提取節(jié)點(diǎn)NDE的Y方向位移UY, 并保存到變量UY3
*SET, U3, (UY3/(1/2.187))	! 給U2變量賦值UY3/(1/2.187)
*STATUS, PARM	! 列表顯示變量參數(shù)的值
Resume, Table_1
*STATUS, PARM	! 列表顯示變量參數(shù)的值
Resume, Table_2
*STATUS, PARM	! 列表顯示變量參數(shù)的值
Resume, Table_3
*STATUS, PARM	! 列表顯示變量參數(shù)的值
FINISH	! 退出后處理模塊