ex5.6(1).inp

ansys有限元分析程序

! 5.6  悬臂长板的大挠度弯曲
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除数据，不加载初始化文件
! (1)设置工程参数、分析类型、单元类型、实常数和材料参数
/FILNAM, EX5.6(1)	! 定义工程文件名称
/PREP7	! 加载前处理模块
/TITLE, EX5.6(1), LARGE DEFLECTION OF A CANTILEVERED PLATE.
/NOPR	! 关闭提示信息
SMRT, OFF	! 关闭智能化功能
ANTYPE, STATIC	! 定义分析类型为静力分析
NLGEOM, ON	! 打开大变形选项
ET, 1, SHELL43	! 定义第1类单元为4节点塑性大应变壳体单元SHELL43
R, 1, 1	! 定义板厚1mm
MP, EX, 1, 1800	! 定义第1类材料的弹性模量EX=1800 N/mm^2
MP, NUXY, , 0	! 定义第1类材料的泊松比NUXY=0
! (2)定义关键点、线和面, 并剖分为面单元
K, 1	! 在坐标原点定义第1个关键点
K, 2, 12	! 在坐标(12mm, 0)定义第2个关键点
K, 3, 12, 1	! 在坐标(12mm, 1mm)定义第3个关键点
K, 4, , 1	! 在坐标(0, 1mm)定义第4个关键点
L, 1, 2	! 过关键点1, 2定义直线
L, 3, 4	! 过关键点3, 4定义直线
LESIZE, ALL, , , 2	! 定义所有线剖分单元时剖分为2段
ESIZE, , 1	! 没有指定剖分段数的线剖分时剖分为1段
A, 1, 2, 3, 4	! 过关键点1, 2, 3, 4定义面
/VIEW, 1, 0.61, -0.64, 0.47	! 设置观察模型的视点
/VUP, 1, Z	! 设置Z轴向上
AMESH, 1	! 对1号面生成面单元网格
! (3)施加位移约束和荷载力
NSEL, S, LOC, X	! 选择X=0位置的节点
D, ALL, ALL	! 固定所有节点的所有自由度
NSEL, S, LOC, X, 12	! 选择X=12mm位置的节点
CP, 1, ROTY, ALL	! 定义自由端节点的绕Y轴的转动自由度为耦合自由度
TORQ=15.708	! 定义集中力偶荷载15.708 N?mm
F, 2, MY, TORQ	! 在2号节点的绕Y轴方向施加力偶15.708 N?mm
NSEL, ALL	! 选择所有结点
FINISH	! 退出前处理模块PREP7
! (4)进入求解选模块, 设置求解选项后, 开始求解
/SOLU	! 进入求解模块
AUTOTS, ON	! 使用自动荷载步长
NSUBST, 20	! 在每一个荷载步最多使用10个荷载子步
LNSRCH, ON	! 使用直线搜索方法
OUTPR, BASIC, LAST	! 输出最后一个荷载子步的基本数据结果
OUTRES, ALL, ALL	! 在求解输出文件中保存每一个荷载子步的结果
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLU
! (5)在后处理模块, 显示应力和位移计算结果
/POST1	! 进入一般后处理模块
RSYS, SOLU	! 选择刚生成的坐标系
SET, 1	! 读入荷载步1的计算结果
SHELL, TOP	! 选择应力输出为壳体上表面
PRNSOL, S, COMP	! 列表显示节点应力结果
PRNSOL, DOF	! 列表显示节点位移结果
! (6)提取结果后存入变量
*GET, UX_N4, NODE, 4, U, X	! 提取4号节点X方向的位移UX, 存入变量UX_N4
*GET, UZ_N4, NODE, 4, U, Z	! 提取4号节点Z方向的位移UZ, 存入变量UZ_N4
*GET, ROTY_N4, NODE, 4, ROT, Y	! 提取4号节点绕Y轴的转角位移ROTY, 存入变量ROTY_N4
*GET, STRSS_N1, NODE, 1, S, X	! 提取1号节点沿着X轴方向的应力SX, 存入变量STRSS_N1
! (7)保存计算结果
SAVE, TABLE_1	! 保存数据到数据库TABLE_1
FINISH	! 退出后处理模块
! (8)重用改用有限应变层合壳单元SHELL181计算
/CLEAR, NOSTART	! 清除所有数据, 不加载初始化文件
/FILNAM, EX5.6(2)	! 定义工程文件名称
/PREP7	! 加载前处理模块
/TITLE, EX5.6(2), LARGE DEFLECTION OF A CANTILEVERED PLATE(SHELL181).
SMRT, OFF	! 关闭智能化功能
/NOPR	! 关闭提示信息
ANTYPE, STATIC	! 定义分析类型为静力分析
NLGEOM, ON	! 打开大变形选项
ET, 1, SHELL181	! 定义第1类单元为有限应变层合壳单元SHELL181
R, 1, 1	! 定义单元厚度为1mm
RMORE, , , , , 10.0	! 定义薄膜单元沙漏参数10.0
MP, EX, 1, 1800	! 定义第1类材料的弹性模量EX=1800 N/mm^2
MP, NUXY, , 0	! 定义第1类材料的泊松比NUXY=0
! (9)定义关键点、线和面, 并将面剖分为面单元
K, 1	! 在坐标原点定义第1个关键点
K, 2, 12	! 在坐标(12mm, 0)定义第2个关键点
K, 3, 12, 1	! 在坐标(12mm, 1mm)定义第3个关键点
K, 4, , 1	! 在坐标(0, 1mm)定义第4个关键点
L, 1, 2	! 过关键点1, 2定义直线
L, 3, 4	! 过关键点3, 4定义直线
LESIZE, ALL, , , 2	! 定义所有线剖分单元时截为2段
ESIZE, , 1	! 没有指定剖分段数的线剖分时截为1段
A, 1, 2, 3, 4	! 过关键点1, 2, 3, 4定义面
AMESH, 1	! 对1号面生成面单元网格
! (10) 施加位移约束和荷载力
NSEL, S, LOC, X	! 选择X=0位置的节点
D, ALL, ALL	! 固定所有节点的所有自由度
NSEL, S, LOC, X, 12	! 选择X=12mm位置的节点
CP, 1, ROTY, ALL	! 定义自由端节点的绕Y轴的转动自由度为耦合自由度
TORQ=7.854	! 定义一半荷载7.854 N?mm
F, 2, MY, TORQ	! 在2号节点的绕Y轴方向施加力偶7.854 N?mm
NSEL, ALL	! 选择所有结点
FINISH	! 退出前处理模块PREP7
! (11) 进入求解模块, 设置求解选项后, 开始求解
/SOLU	! 进入求解模块
AUTOTS, ON	! 使用自动荷载步长
NSUBST, 10	! 在每一个荷载步最多使用10个荷载子步
LNSRCH, ON	! 使用直线搜索方法
OUTPR, BASIC, LAST	! 输出最后一个荷载子步的基本数据结果
OUTRES, ALL, ALL	! 在求解输出文件中保存每一个荷载子步的结果
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLUTION
! (12)在求解模块重启动分析
/SOLU	! 加载求解模块
ANTYPE, , REST	! 重启动分析
F, 2, MY, TORQ*2	! 给2号节点施加全部荷载2*7.854 N?mm
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLU
! (13)在后处理模块, 显示计算结果
/POST1	! 进入一般后处理模块
RSYS, SOLU	! 选择刚生成的坐标系
SET, 2 	! 读入荷载步2(从重启动分析)的计算结果
SHELL, TOP	! 选择应力输出的壳体上表面
PRNSOL, S, COMP	! 列表显示节点应力结果
PRNSOL, DOF	! 列表显示节点位移结果
*GET, UX_N4, NODE, 4, U, X	! 提取4号节点X方向的位移UX, 存入变量UX_N4
*GET, UZ_N4, NODE, 4, U, Z	! 提取4号节点Z方向的位移UZ, 存入变量UZ_N4
*GET, ROTY_N4, NODE, 4, ROT, Y	! 提取4号节点绕Y轴的转角位移ROTY, 存入变量ROTY_N4
*GET, STRSS_N1, NODE, 1, S, X	! 提取1号节点沿着X轴方向的应力SX, 存入变量STRSS_N1
! (14) 显示计算结果
*STATUS, PARM	! 显示保存在变量中的结果
Resume, TABLE_1	! 从数据库TABLE_1恢复数据
*STATUS, PARM	! 显示保存在变量中的结果
Finish