ex4.11.inp

ansys有限元分析程序

! 4.11  纯弯曲悬臂曲梁的二维静力分析
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH	! 退出以前模块
/CLEAR, NOSTART	! 清除数据，不加载初始化文件
! (1)设置工程选项，定义分析类型、单元类型和材料参数
/FILNAME, EX4.11	! 定义工程文件名称
/PREP7	! 加载前处理模块
SMRT, OFF	! 关闭智能化功能
/TITLE, EX4.11, BENDING OF A CURVED BEAM.
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, PLANE2, , , 3	! 定义第1类单元为带厚度的平面应力单元PLANE2
ET, 2, BEAM3	! 定义第2类单元为2D弹性梁单元BEAM3
R, 1, 25.4	! 定义第1类实常数, 单元厚度25.4mm
R, 2, 1.0, 1, 1.0	! 定义第2类实常数, 梁的截面积是1.0, 惯性矩是1, 高度是1.0
MP, EX, 1, 68950	! 定义第1类材料的弹性模量EX = 68.95 GPa
MP, NUXY, 1, 0.25	! 定义第1类材料的泊松比NUXY＝0.25
! (2)在柱坐标系下定义曲梁上的4个关键点, 自由端处的线, 通过线单元剖分生成梁单元
CSYS, 1	! 切换到柱坐标系
K, 1, 89	! 在坐标(89mm，0)处定义第1个关键点
K, 2, 89, 90	! 在坐标(89mm，90度)处定义第2个关键点
KGEN, 2, 1, 2, 1, 25.4	! 循环2次, 按照1号到2号节点, 增量为1的节点范围, 每个节点的半径坐标增加25.4mm, 生成3号节点和4号节点.
L, 2, 4	! 过2号节点和4号节点定义线, 以定义刚性梁单元
LESIZE, 1, , , 4	! 指定1号线的单元剖分段数为4
TYPE, 2	! 切换单元类型为第2类单元2D弹性梁单元BEAM3
REAL, 2	! 切换实常数类型为第2类梁的实常数
LMESH, 1	! 将1号线剖分为梁单元
! (3)过关键点定义曲梁的内侧和外侧圆弧线，面单元剖分曲梁面区域后得到有限元模型
L, 1, 2	! 过关键点1, 2定义曲梁所在的内弧线
LESIZE, 2, , , 25	! 指定2号圆弧线的单元剖分段数为25
L, 3, 4	! 过3号节点和4号节点定义曲梁所在的外弧线
LESIZE, 3, , , 25	! 指定3号线的单元剖分段数为25
ESIZE, , 4	! 指定单元剖分段数为4
A, 1, 2, 4, 3	! 过节点1, 2, 4和3定义面
AATT, 1, 1, 1	! 将第1类材料参数, 第1类实常数和第1类单元等属性赋值给前面定义的面
/DEVICE, VECTOR, 1	! 打开矢量图显示开关, 便于单元网格的显示
AMESH, 1	! 对1号面执行单元剖分
! (4)通过关键点定义位移约束和荷载
NROTAT, ALL	! 将所有节点的节点坐标系旋转到柱坐标系方向
DK, 1, ALL, , , 1	! 约束第1个关键点的所有自由度, 并指定为约束展开点
DK, 3, UY, , , 1	! 约束第3个关键点的UY环向自由度, 并指定为约束展开点
FK, 2, FY, 452	! 在第2号关键点的环向逆时针方向施加集中力452N
FK, 4, FY, -452	! 在第4号关键点的环向顺时针方向施加集中力452N
FINISH	! 退出前处理模块PREP7
! (5)进入求解模块求解
/SOLU	! 加载求解模块
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLU
! (6)在后处理模块中显示指定结点位置的应力计算结果
/POST1	! 进入一般后处理模块
RSYS, 1	! 将结果坐标系切换到柱坐标系
NSEL, R, LOC, X, 89	! 选择曲梁内侧R＝89mm位置的节点
NSEL, U, LOC, Y, 90	! 取消选择曲梁自由端Y＝90度的节点
PRNSOL, S, COMP	! 列表显示曲梁内侧结点应力计算结果
NSORT, S, Y	! 将节点计算结果按曲梁内侧环向应力SY排序
*GET, SI, SORT, , MAX	! 提取排序后的最大环向应力, 存入变量SI
NSEL, ALL	! 选择所有节点
NSEL, R, LOC, X, 114	! 选择曲梁外侧R＝114mm的节点
NSEL, U, LOC, Y, 90	! 取消选择曲梁自由端Y＝90度的节点
PRNSOL, S, COMP	! 列表显示曲梁外侧结点应力计算结果
NSORT, S, Y	! 将节点计算结果按曲梁外侧环向应力SY排序
*GET, SO, SORT, , MIN	! 提取排序后的最小环向应力, 存入变量SO
/DEVICE, VECTOR, 0	! 关闭矢量图显示开关, 便于节点应力色度图的显示
PLNSOL, S, Y	! 图形显示环向应力分布
FINISH	! 退出后处理模块