! 4.7 一对集中力作用下的圆环
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH ! 退出以前模块
/CLEAR, NOSTART ! 清除内存中的所有数据, 不读入初始化文件
! （1）定义文件名、标题、在前处理模块中通过循环在柱坐标系下定义节点
/FILNAME, EX4.7(1) ! 定义工程文件名称
/TITLE, EX4.7(1), Round Ring under concentrated force by direct modeling.
/PREP7 ! 进入前处理模块PREP7
CSYS, 1 ! 切换当前坐标系为柱坐标系
*DO, i, 1, 36, 1 ! 环向循环36次生成节点
*DO, j, 1, 11, 1 ! 半径方向循环11次生成节点
N, i+(j-1)*36, 50+5*(j-1), (i-1)*10 ! 定义所有节点的半径位置和角度位置
*ENDDO ! 结束循环
*ENDDO ! 结束循环
! (2)定义单元类型、实常数和材料参数, 利用循环定义单元
ET, 1, PLANE82,,, 3 ! 定义第1类单元为带厚度的四节点平面应力单元PLANE82
R, 1, 1 ! 定义单元的第1类实常数:单元厚度为1mm
MP, EX, 1, 2e3 ! 定义第1类材料的弹性模量EX＝2 GPa=2e3 N/mm^2
MP, PRXY, 0.3 ! 定义第1类材料的泊松比PRXY＝0.3
*DO, j, 1, 10, 1 ! 半径方向循环10次生成单元
*DO, i, 1, 35, 1 ! 环向循环35次生成单元
E, i+(j-1)*36, i+(j-1)*36+1, i+(j-0)*36+1, i+(j-0)*36 ! 定义单元的4个节点
*ENDDO ! 结束内层循环
E, 36+(j-1)*36, 1+(j-1)*36, 1+(j-0)*36, 36+(j-0)*36 ! 定义最后一列单元
*ENDDO ! 结束外层循环
FINISH ! 退出前处理模块
! (3)施加位移边界条件后，求解
/SOLU ! 进入求解模块
D, 379, ALL, 0 ! 限制最左侧的379号节点的所有位移
D, 361, uy, 0 ! 限制最右侧361号节点的Y方向位移
D, 361, ux, -1 ! 给最右侧361号节点施加-X方向的大小为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 ! 退出后处理模块
