ex5.10(1).inp

来自「ansys有限元分析程序」· INP 代码 · 共 99 行

! 5.10  圆筒在一对横向集中力作用下的变形
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除数据, 不加载初始化文件
! （1）设置工程选项、分析类型、实常数、材料参数
/FILNAME, EX5.10(1)	! 定义工程文件名称
/PREP7	! 加载前处理模块
SMRT, OFF	! 关闭智能化网格功能
/TITLE, EX5.10(1), PINCHED CYLINDER UNDER CONCENTRATED FORCE.
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, SHELL93	! 定义第1类单元为8节点结构壳体单元SHELL93
R, , 2.3876	! 定义板厚为2.3876mm
MP, EX, , 72.4E3	! 定义第1类材料的弹性模量EX＝72.4 GPa
MP, NUXY, , 0.3125	! 定义第1类材料的泊松比NUXY＝0.3125
! （2）在柱坐标系下定义关键点、柱面并对其进行单元剖分
CSYS, 1	! 切换坐标系到柱坐标系
K, 1, 125.8	! 在坐标(125.8mm, 0, 0)处定义第1个关键点
K, 2, 125.8, , 131.5	! 在坐标(125.8mm, 0, 131.5mm)处定义第2个关键点
KGEN, 2, 1, 2, 1, , 90	! 按照1号节点和2号节点模式, Y坐标增量90度, 生成3号和4号节点
A, 1, 2, 4, 3	! 过节点1, 2, 4, 3定义面
ESIZE, , 8	! 定义单元剖分时剖分段数为8
/View, 1, 1, 1, 1	! 将窗口1的视点切换到等轴测方向
AMESH, 1	! 对1号面执行面单元剖分操作
! (3)在直角坐标系下定义约束
CSYS, 0	! 切换坐标系到整体直角坐标系
NSEL, S, LOC, X, 0	! 选择X=0位置的节点
DSYM, SYMM, X, 0	! 定义关于X轴(YZ平面)的位移对称条件
NSEL, S, LOC, Y, 0	! 选择Y=0位置的节点
DSYM, SYMM, Y, 0	! 定义关于Y轴(XZ平面)的位移对称条件
NSEL, S, LOC, Z, 0	! 选择Z=0位置的节点
DSYM, SYMM, Z, 0	! 定义关于Z轴(XY平面)的位移对称条件
NSEL, ALL	! 选择所有结点
FK, 3, FY, -111.2	! 在3号关键点施加沿着-Y方向的集中力111.2N
FINISH	! 退出前处理模块PREP7
! (4)在求解模块求解
/SOLU	! 加载求解模块
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLUTION
! (5)在后处理模块显示节点位移计算结果
/POST1	! 进入一般后处理模块
NSEL, S, LOC, Y, 125.8	! 选择Y=125.8mm位置的节点
NSEL, R, LOC, Z, 0	! 继续选择Z=0位置的节点
NSEL, R, LOC, X, 0	! 继续选择X=0位置的节点
PRNSOL, U, COMP	! 打印位移结果和位移矢量和
TOP_NODE = NODE(125.8, 90, 0)	! 选择距离(125.8mm, 90度, 0)最近的节点, 存入变量TOP_NODE
*GET, DISP, NODE, TOP_NODE, U, Y	! 提取节点TOP_NODE的位移UY, 存入变量DISP
! （6）将计算结果保存到文件
SAVE, TABLE_1	! 保存数据到数据库TABLE_1
FINISH	! 退出一般后处理模块
! (7)改变单元重新计算
/CLEAR, NOSTART	! 清除所有数据, 不加载初始化文件
/FILNAME, EX5.10(2)	! 定义工程文件名称
/PREP7	! 加载前处理模块
/TITLE, EX5.10(2), PINCHED CYLINDER UNDER CONCENTRATED FORCE(SHELL150).
SMRT, OFF	! 关闭智能化网格功能
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, SHELL150	! 定义第1类单元为8节点结构壳体单元SHELL150
R, , 2.3876	! 定义壳体单元厚度为2.3876mm
MP, EX, , 72.4E3	! 定义第1类材料的弹性模量EX=72.4 GPa
MP, NUXY, , 0.3125	! 定义第1类材料的泊松比NUXY=0.3125
! （8）在柱坐标系下定义关键点, 柱面并对其进行单元剖分
CSYS, 1	! 切换坐标系到柱坐标系
K, 1, 125.8	! 在坐标(125.8mm, 0, 0)处定义第1个关键点
K, 2, 125.8, , 131.5	! 在坐标(125.8mm, 0, 131.5mm)处定义第2个关键点
KGEN, 2, 1, 2, 1, , 90	! 按照1号节点和2号节点模式, Y坐标增量90度, 生成3号和4号节点
A, 1, 2, 4, 3	! 过节点1, 2, 4, 3定义面
ESIZE, , 8	! 定义单元剖分时剖分段数为8
/View, 1, 1, 1, 1	! 将窗口1的视点切换到等轴测方向
AMESH, 1	! 对1号面执行面单元剖分
! （9）在直角坐标系下定义约束
CSYS, 0	! 切换坐标系到整体直角坐标系
NSEL, S, LOC, X, 0	! 选择X=0位置的节点
DSYM, SYMM, X, 0	! 定义关于X轴(YZ平面)的位移对称条件
NSEL, S, LOC, Y, 0	! 选择Y=0位置的节点
DSYM, SYMM, Y, 0	! 定义关于Y轴(XZ平面)的位移对称条件
NSEL, S, LOC, Z, 0	! 选择Z=0位置的节点
DSYM, SYMM, Z, 0	! 定义关于Z轴(XY平面)的位移对称条件
NSEL, ALL	! 选择所有结点
FK, 3, FY, -111.2	! 在3号关键点施加沿着-Y方向的集中力111.2N
FINISH	! 退出前处理模块PREP7
! (10) 在求解模块求解
/SOLU	! 加载求解模块
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLUTION
! (11) 在后处理模块显示节点位移计算结果
/POST1	! 进入一般后处理模块
SET, 1	! 读入第1荷载步的计算结果
NSEL, S, LOC, Y, 125.8	! 选择Y=125.8mm位置的节点
NSEL, R, LOC, Z, 0	! 继续选择Z=0位置的节点
NSEL, R, LOC, X, 0	! 继续选择X=0位置的节点
PRNSOL, U, COMP	! 打印位移结果和位移矢量和
TOP_NODE = NODE(125.8, 90, 0)	! 选择距离(125.8mm, 90度, 0)最近的节点, 存入变量TOP_NODE
*GET, DISP, NODE, TOP_NODE, U, Y	! 提取节点TOP_NODE的位移UY, 存入变量DISP
*STATUS, PARM	! 显示保存在变量中的结果
Resume, Table_1	! 从数据库文件中恢复结果
*STATUS, PARM	! 显示保存在变量中的结果
FINISH	! 退出后处理模块