! 4.7  一對集中力作用下的圓環
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!．用二維幾何體素建立圓環的求解程序。
FINISH        ! 退出以前模塊
/CLEAR,NOSTART       ! 清除內存中的所有數據,不讀入初始化文件
! （1）定義文件名、標題、在前處理模塊中通過循環在柱坐標系下定義節點
/FilName,EX4.7(2)       ! 定義工程文件名稱
/TITLE,EX4.7(2), Round Ring under concentrated force by geometry modeling.
/PREP7      ! 進入前處理模塊PREP7
CYL4, , ,50, ,100 ! 定義以坐標原點（0，0）為圓心，內半徑50mm，外半徑100mm的圓環
! (2)定義單元類型、實常數和材料參數,利用循環定義單元
ET,1,PLANE82,,,3      ! 定義第1類單元為帶厚度的四節點平面應力單元PLANE82
R,1,1      ! 定義單元的第1類實常數:單元厚度為1mm
MP,EX,1,2e3      ! 定義第1類材料的彈性模量EX＝2 GPa=2e3 N/mm2
MP,PRXY,0.3      ! 定義第1類材料的泊松比PRXY＝0.3
esize,5        ! 設置單元剖分時的默認單元尺寸為5mm
amesh,1      ! 對圓環面進行單元剖分
FINISH      ! 退出前處理模塊
! (3)施加位移邊界條件后，求解
/SOLU       ! 進入求解模塊
dk,3,all      ! 限制最左側的3號關鍵點的所有位移
dk,1,uy       ! 限制最右側點1號關鍵點的Y方向位移
dk,1,ux,-1   ! 給最右側點1號關鍵點施加-X方向的大小為1的強制位移UX=1
SOLVE       ! 開始靜力分析求解
FINISH      ! 退出求解模塊
! (4)在后處理模塊,顯示變形和應力結果
/POST1      ! 進入后處理模塊
PLDISP,2      ! 帶原模型圖,顯示變形圖
PLNSOL,S,X      ! 顯示水平方向的應力SX
PLNSOL,S,Y      ! 顯示豎直方向應力SY
PLNSOL,S,XY      ! 顯示XY平面內的剪應力SXY
PLNSOL,S,1      ! 顯示第1主應力S1
PLNSOL,S,3      ! 顯示第3主應力S3
PLNSOL,S,EQV     ! 顯示Mises應力SEQV
FINISH       ! 退出后處理模塊