ex4.1.inp

ansys有限元分析程序

! 4.1  带三个圆孔的平面支座分析
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH	! 退出以前模块
/CLEAR, NOSTART	! 清除内存中的所有数据, 不读入初始化文件
/FILNAME, EX4.1	! 定义工程文件名称
/TITLE, EX4.1, 2D Plane Stress Bracket.
/PREP7	! 进入前处理模块PREP7
! （1）定义外轮廓区域
BLC4, 0, 0, 80, 100	! 在（0，0）定义宽度80mm，高度100mm的主矩形区域
BLC4, -20, 20, 20, 60	! 在（-20mm，20mm）定义宽度20mm，高度60mm的小矩形区域
CYL4, 80, 50, 50	! 以（80mm，50mm）为中心，50mm为半径定义大圆形区域
CYL4, 0, 20, 20	! 以（0，20mm）为中心，20mm为半径的定义左下方小圆形区域
CYL4, 0, 80, 20	! 以（0，80mm）为中心，20mm为半径的定义左上方小圆形区域
AADD, ALL	! 将上面定义的5块面积相加得到6号面积
! （2）定义圆孔轮廓，通过剖切生成这3个孔
CYL4, 80, 50, 30	! 以（80mm，50mm）为中心，30mm为半径定义将要剖切的大圆孔区域
CYL4, 0, 20, 10	! 以（0，20mm）为中心，10mm为半径定义将要剖切的左下小圆孔区域
CYL4, 0, 80, 10	! 以（0，80mm）为中心，10mm为半径定义将要剖切的左上小圆孔区域
ASBA, 6, ALL	! 从6号面积中减去3个小孔区域
! （3）定义单元和材料的力学特性
ET, 1, PLANE82	! 定义第1类单元为平面单元PLANE82
KEYOPT, 1, 3, 3	! 定义单元选项是带厚度的平面应力单元
R, 1, 20	! 定义单元的第1类实常数: 支座的厚度20mm
MP, EX, 1, 200000	! 定义第1类材料的弹性模量EX=200GPa（指令格式：MP, Lab, MAT, C0）
MP, PRXY, 1, 0.3	! 定义第1类材料的泊松比PRXY=0.3
! （4）设置单元剖分的尺寸并进行单元剖分
AESIZE, ALL, 5	! 剖分单元尺寸设置为5mm
AMESH, ALL	! 对所有面进行单元剖分
! （5）在求解模块中定义位移约束，定义荷载并求解
/SOLU	! 进入求解模块SOLU
ANTYPE, 0	! 申明求解类型是静力分析
! 固定下面小圆孔边界上所有节点的自由度，7，8，9和10号直线就是这个孔边的4条线
DL, 7, , ALL, 0	! 固定7号直线的所有自由度
DL, 8, , ALL, 0	! 固定8号直线的所有自由度
DL, 9, , ALL, 0	! 固定9号直线的所有自由度
DL, 10, , ALL, 0	! 固定10号直线的所有自由度
! 固定上面小圆孔边界上所有节点的自由度，11，12，13和14号直线就是这个孔边的4条线
DL, 11, , ALL, 0	! 固定11号直线
DL, 12, , ALL, 0	! 固定12号直线
DL, 13, , ALL, 0	! 固定13号直线
DL, 14, , ALL, 0	! 固定14号直线
FK, 9, FY, -1000	! 在9号关键点上施加-Y方向的大小是1000的集中力，9号关键点位于大圆孔最低点
SOLVE	! 开始求解
FINISH
! （6）在后处理模块中，显示等效应力MISES应力，并用动画方式显示该应力分布随着荷载的变化
/POST1	! 进入后处理模块POST1
PLNSOL, S, EQV	! 显示所有结点上的Mises应力
ANCNTR, 10, 0.15	! 用10帧每隔0.15秒显示一帧的速度动画显示Mises应力
FINISH