! 4.8  用實體單元分析變截面桿的拉伸
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!
FINISH	! 退出以前模塊
/CLEAR, NOSTART	! 清除數據，不加載初始化文件
! (1)設置工程選項, 定義分析類型、單元類型和材料參數
/FILNAME, EX4.8	! 定義工程文件名稱
/PREP7	! 加載前處理模塊
SMRT, OFF	! 關閉智能化功能
/TITLE, EX4.8, ELONGATION OF A SOLID BAR WITH VARIED SECTION.
ANTYPE, STATIC	! 定義分析類型為靜力分析
ET, 1, SOLID45	! 定義第1類單元為三維結構固體單元SOLID45
MP, EX, 1, 70.71E3	! 定義第1類材料的彈性模量EX=70.71 GPa
MP, NUXY, 1, 0.3	! 定義第1類材料的泊松比NUXY=0.3
! (2)定義固定端4個點和自由端4個點，過這8個點定義體，對體剖分單元
/VIEW, 1, 1, 1, 1	! 改變視點為(1, 1, 1), 沿著等軸測方向
K, 1, 25.4, , 25.4	! 在坐標（25.4mm，0，25.4mm）處定義第1個關鍵點
K, 2, -25.4, , 25.4	! 在坐標（-25.4mm，0，25.4mm）處定義第2個關鍵點
K, 3, -25.4, , -25.4	! 在坐標（-25.4mm，0，-25.4mm）處定義第3個關鍵點
K, 4, 25.4, , -25.4	! 在坐標（25.4mm，0，-25.4mm）處定義第4個關鍵點
K, 5, 12.7, 254, 12.7	! 在坐標（12.7mm，254mm，12.7mm）處定義第5個關鍵點
K, 6, -12.7, 254, 12.7	! 在坐標（-12.7mm，254mm，12.7mm）處定義第6個關鍵點
K, 7, -12.7, 254, -12.7	! 在坐標（-12.7mm，254mm，-12.7mm）處定義第7個關鍵點
K, 8, 12.7, 254, -12.7	! 在坐標（12.7mm，254mm，-12.7mm）處定義第8個關鍵點
V, 1, 2, 3, 4, 5, 6, 7, 8	! 過前面定義的8個關鍵點定義體
LSEL, S, LINE, , 5, 11, 2	! 從5號線到11號線，增量2選擇5，7，9和11號線
LESIZE, ALL, , , 7	! 指定前面選中的線在剖分單元是剖分為7段
LSEL, ALL 	! 選擇所有的線
ESIZE, , 1 	! 單元剖分段數為1
/OUT, SCRATCH	! 將輸出結果重定向到SCRATCH文件，以避免屏幕輸出
VMESH, 1	! 將1號體剖分為體單元
/OUT	! 結束數據重定向，關閉輸出文件
! (3)施加位移約束和荷載
OUTPR, BASIC, ALL	! 設置輸出所有荷載步的基本信息
NSEL, S, LOC, Y, 0	! 選擇Y=0位置的所有節點，即固定端位置的節點
D, ALL, ALL 	! 位移約束選擇集中的所有節點的所有自由度
NSEL, ALL	! 選擇所有節點
NSEL, S, LOC, Y, 254	! 選擇Y=254mm位置的所有節點，即自由端節點
SF, , PRES, -68.949	! 給節點指定表面荷載44483N/25.4mm^2=68.949MPa
NSEL, ALL	! 選擇所有節點
FINISH	! 退出前處理模塊PREP7
! (4)在求解模塊，開始靜力分析的求解
/SOLU	! 加載求解模塊
SOLVE	! 開始求解
FINISH	! 退出求解模塊SOLUTION
! (5)在后處理模塊中，顯示Y方向應力云圖
/POST1	! 進入一般后處理模塊
ETABLE, SIGY, S, Y	! 以Y方向的正應力為內容，定義單元表SIGY
/VIEW, 1, 1	! 改變視點為(1, 0, 0), 沿著X軸的負方向
/VUP, 1, -Y	! 設置-Y軸向上
NSLE, S    	! 選擇和單元相連接的節點
PLNSOL, S, Y  	! 顯示軸向應力
ESEL, S, ELEM, , 4	! 選擇中間的第4號單元
PRETAB, SIGY 	! 顯示單元表中的單元應力計算結果
PRNSOL, S, COMP	! 顯示節點應力計算結果
ESEL, ALL  	! 選擇所有單元
NSEL, S, LOC, Y, 254	! 選擇位于Y＝254mm的節點（自由端）
PRNSOL, DOF  	! 顯示節點位移計算結果
NSEL, ALL	! 選擇所有節點
! (6)選擇指定位置的節點和單元，并提取它們的位移和應力
/NOPR	! 關閉輸入反饋提示信息
MID_NODE = NODE(0, 128, 0)	! 選擇距離（0, 128mm, 0）最近的節點，命名為MID_NODE
MID_ELM = ENEARN(MID_NODE)	! 選擇距離MID_NODE節點最近的單元，命名為MID_ELM
BOT_NODE = NODE (0, 254, 0)	! 選擇距離（0, 254mm, 0）最近的節點，命名為BOT_NODE
*GET, DEF, NODE, BOT_NODE, U, Y	! 提取BOT_NODE節點的位移UY，存入變量DEF
*GET, STRSS, ELEM, MID_ELM, ETAB, SIGY	! 在單元表SIGY中提取MID_ELM單元的應力，存入變量STRSS
*STATUS, PARM	! 顯示所有內存變量的值
FINISH	! 退出后處理模塊