! 4.2 角支座应力分析
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH ! 退出以前模块
/CLEAR, NOSTART ! 清除内存中的所有数据, 不读入初始化文件
/FILNAME, EX4.2 ! 定义工程文件名称
/TITLE, EX4.2, Static Analysis of a Corner Bracket.
/PREP7 ! 进入前处理模块PREP7
! （1）定义外轮廓区域
BLC4, 0, -25, 150, 50 ! 定义左下角（0, -25mm）, 宽度150mm, 高度50mm的矩形块（图4.7的A1）
BLC4, 100, -25, 50, -50 ! 定义左下角（100mm, -25mm）, 宽度50mm, 高度-50mm的矩形（图4.7的A2）
CYL4, 0, 0, 0, 0, 25, 360, 0 ! 定义圆心（0, 0）, 内半径0, 外半径25mm的360度的圆（图4.7的A3）
CYL4, 125, -75, 0, 0, 25, 360, 0 ! 定义圆心（125mm, -75mm）, 内半径0, 外半径25mm的360度的圆形区域（图4.7中的A4）
AADD, ALL ! 将定义的4块面A1, A2, A3和A4相加得到5号面
! （2）定义开孔位置轮廓线, 利用面积相减生成开孔
CYL4, 0, 0, 0, 0, 12.5, 360, 0 ! 定义圆心在（0, 0）, 内半径0, 外半径12.5mm的360度的圆, 定义左边开孔位置的小圆孔
CYL4, 125, -75, 0, 0, 12.5, 360, 0 ! 定义圆心在（125mm, -75mm）, 内半径0, 外半径12.5mm的360度的圆, 定义右边开孔位置的小圆孔
ASBA, 5, ALL ! 从5号大面积中减去两个小圆孔
! （3）利用线倒角生成倒角区域
LFILLT, 17, 8, 10 ! 过17号和8号直线, 定义半径为10mm的内侧位置的倒角曲线
AL, 12, 13, 14 ! 用3条线12, 13和14定义倒角后形成的小面积
AADD, ALL ! 将当前所有面相加
! （4）定义单元和材料的力学特性
ET, 1, Plane82, , , 3 ! 定义第1类单元为带厚度的平面单元PLANE82
MP, EX, 1, 210e3 ! 定义第1类材料的弹性模量EX=210GPa
MP, PRXY, 1, 0.3 ! 定义第1类材料的泊松比PRXY=0.3
R, 1, 12.5 ! 定义单元的第1类实常数:厚度 12.5mm
! （5）设置单元剖分的尺寸并进行单元剖分
Esize, 5 ! 定义单元剖分的尺寸为5mm
Amesh, ALL ! 对所有面进行单元剖分
! （6）在求解模块/SOLU中定义位移约束, 定义荷载并求解
/SOLU ! 进入求解模块SOLU
! 固定左边小圆孔边界上所有节点的自由度, 1, 4, 5和6号直线就是这个孔边的4条线
DL, 1, , ALL ! 对1号线定义位移约束
DL, 6, , ALL ! 对6号线定义位移约束
DL, 4, , ALL ! 对4号线定义位移约束
DL, 5, , ALL ! 对5号线定义位移约束
! 在右边小圆孔下面边界上施加线性变化的压力荷载
SFL, 10, PRES, 1, 10 ! 在10号线上定义线性变化的分布荷载, 集度从1N/mm变化到10N/mm
SFL, 11, PRES, 10, 1 ! 在11号线上定义线性变化的分布荷载, 集度从10N/mm变化到1N/mm
SOLVE ! 开始求解
FINISH ! 退出后模块
! （7）在一般后处理模块POST1中显示计算结果
/POST1 ! 进入后处理模块POST1
SET, 1 ! 读入第1个荷载步的计算结果
PLNSOL, S, EQV ! 用色度图显示节点上的Mises应力
ANCNTR, 6, 0.15 ! 用动画显示节点上的Mises应力
FINISH
