ex4.13(1).inp

ansys有限元分析程序

! 4.13  均匀拉力作用下含圆孔板的孔边应力集中
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH	! 退出以前模块
/CLEAR, NOSTART	! 清除数据, 不加载初始化文件
! （1）设置工程选项、分析类型、单元类型和材料参数
/FILNAME, EX4.13(1)	! 定义工程文件名称
/PREP7	! 加载前处理模块
SMRT, OFF	! 关闭智能化功能
/TITLE, EX4.13(1), STRESS CONCENTRATION AT A HOLE IN A PLATE with PLANE2.
/NOPR	! 关闭输入信息反馈提示
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, PLANE2	! 定义第1类单元为二维6节点三角形固体单元PLANE2
MP, EX, 1, 207E3	! 定义第1类材料的弹性模量EX ＝ 207 GPa
MP, NUXY, 1, 0.3	! 定义第1类材料的泊松比NUXY ＝ 0.3
! （2）定义关键点
K, 1, 152.4	! 在坐标(152.4mm, 0)定义第1个关键点
K, 2, 152.4, 152.4	! 在坐标(152.4mm, 152.4mm)定义第2个关键点
K, 3, , 152.4	! 在坐标(0, 152.4mm)定义第3个关键点
K, 4, , 12.7	! 在坐标(0, 12.7mm)定义第4个关键点
K, 5, 12.7	! 在坐标(12.7mm, 0)定义第5个关键点
K, 6	! 在坐标(0, 0)定义第6个关键点
! （3）定义线和线单元剖分的尺寸
L, 1, 2	! 过关键点1, 2定义直线
L, 2, 3	! 过关键点2, 3定义直线
L, 3, 4	! 过关键点3, 4定义直线
LESIZE, 3, , , 4, .25	! 给3号线定义剖分段数为4和间距比例0.25
LARC, 4, 5, 6, 12.7	! 过4号点和5号点, 以6号点为圆心, 以12.7mm为半径定义圆弧
LESIZE, 4, , , 6	! 给4号线定义剖分段数为6
L, 5, 1	! 过关键点5, 1定义直线
LESIZE, 5, , , 4, 4	! 给5号线定义剖分段数为4和间距比例为4
! （4）定义面并将其剖分为面单元, 图形显示单元布置
AL, 1, 2, 3, 4, 5	! 用线1, 2, 3, 4和5定义面
ESIZE, , 4 	! 单元剖分尺寸为分为4段
AMESH, ALL	! 将所有的面剖分为面单元
/AUTO, 1	! 设置模型显示的最佳比例
/PLOPTS, INFO, 0	! 关闭图形窗口中图例信息显示栏
/PLOPTS, WINS, 0	! 关闭图形窗口的自适应调整
/WINDOW, , LTOP	! 在图形窗口的左上角定义窗口1
EPLOT	! 图形显示单元布置
! （5）选择相应的线, 施加位移和荷载边界条件
LSEL, S, LINE, , 3, 5, 2	! 选择3号线和5号线
DL, ALL, , SYMM	! 对所有的选择集中的线施加对称位移约束
LSEL, S, LINE, , 1	! 选择1号线
NSLL, , 1	! 选择1号线上的节点
SF, ALL, PRES, -6.895	! 施加板端的拉力荷载6.895 MPa
LSEL, ALL	! 选择所有的线
NSEL, ALL	! 选择所有节点
CSYS, 1	! 将坐标系切换到柱坐标系
FINISH	! 退出前处理模块PREP7
! （6）进入求解模块求解
/SOLU	! 加载求解模块
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLUTION
SAVE	! 保存数据到数据库EX4.13(1).DB
! （7）进入后处理, 显示计算结果
/POST1	! 进入一般后处理模块
SET, 1, 1	! 读入第1个荷载步的计算结果
NSORT, S, X, , , 3	! 将计算结果按SX排序, 只保留前3项最高的结果
PRNSOL, S, COMP	! 列表显示应力分量计算结果
/WINDOW, 1, OFF	! 关闭1号窗口的显示
/NOERASE	! 不擦除图形窗口
/DSCALE, 2, 1	! 用缩放位移约束命令, 缩放位移约束实部2倍, 虚部1倍
/WINDOW, 2, RTOP	! 在图形窗口的右上角定义窗口2
PLNSOL, S, X	! 图形显示水平方向的节点应力SX分布图
PLNSOL, S, Y	! 图形显示竖直方向的节点应力SY分布图
*GET, CRSESTR, NODE, 18, S, X	! 提取18号节点的应力SX, 存入变量CRSESTR
*STATUS, PARM	! 显示内存变量中的值
! （8）将计算结果保存到文件TABLE_1
SAVE, TABLE_1	! 保存数据到数据库TABLE_1
! （9）改变单元类型为PLANE42, 用子模型技术重新计算
FINISH	! 退出后处理模块
/CLEAR, NOSTART	! 清除数据, 不加载初始化文件! 
/FILNAME, 4.13(2)	! 定义子结构分析文件名
/PREP7	! 加载前处理模块
SMRT, OFF	! 关闭智能化网格功能
/NOPR	! 关闭输入信息反馈提示
/TITLE, EX4.13(2), STRESS CONCENTRATION AT A HOLE IN A PLATE With PLANE42 SUB-MODEL.
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, PLANE42	! 定义第1类单元为二维等参固体单元PLANE42
MP, EX, 1, 207E3	! 定义第1类材料的弹性模量EX＝207 GPa
MP, NUXY, 1, 0.3	! 定义第1类材料的泊松比NUXY＝0.3
! (10) 在柱坐标系下定义孔边子模型上的关键点
CSYS, 1	! 切换到柱坐标系下
K, 10, 12.7, 45	! 为子模型在坐标（12.7mm, 45度）定义第10个关键点
K, 11, 12.7, 90	! 在坐标（12.7mm, 90度）定义第11个关键点
K, 12, 38.1, 45	! 在坐标（38.1mm, 45度）定义第12个关键点
K, 13, 38.1, 90	! 在坐标（38.1mm, 90度）定义第13个关键点
A, 10, 12, 13, 11	! 关键点10, 12, 13和11定义面
ESIZE, , 8	! 单元剖分为8格
MSHK, 1 	! 设定映射网格选项
MSHA, 0, 2D	! 使用四边形网格
AMESH, 1	! 按照上面选项, 对1号面积进行网格剖分
! (11) 设置窗口选项后, 绘制单元布置
/WINDOW, 2, OFF	! 关闭2号窗口的显示
/ERASE 	! 擦除屏幕
/PLOPTS, INFO, 0	! 关闭图形窗口中图例信息显示栏
/PLOPTS, WINS, 0	! 关闭图形窗口的自适应调整
/WINDOW, , LTOP	! 在图形窗口的左上角定义窗口1
/USER   	! 使用用户自定义缩放比例
/DIST, 1, 18	! 用户自定义1号窗口的观察距离18mm
/FOCUS, 1, 14, 25	! 用户自定义1号窗口的焦点(14mm, 25mm)
EPLOT	! 图形显示单元布置
LSEL, S, LINE, , 1, 2	! 选择1号线和2号线
NSLL, , 1 	! 选择1号线上的节点
NWRITE 	! 将节点数据写入到文件 4.13(2).NODB
LSEL, ALL	! 选择所有的线
NSEL, ALL	! 选择所有节点
FINISH	! 退出前处理模块PREP7
SAVE	! 保存子模型数据到数据库4.13(2).DB
! (12) 在后处理模块恢复模型数据, 读入计算结果
/POST1	! 进入一般后处理模块
RESUME, EX4.13(1).db	! 从数据库EX4.13(1).DB中恢复平板模型数据
FILE, EX4.13(1), rst	! 定义结果文件
CBDOF, , , , , , , 0, , 0	! 激活剖切边界插值
! (13) 重新读模型, 施加不同的边界条件
FINISH	! 退出一般后处理模块
/PREP7	! 加载前处理模块
SMRT, OFF	! 关闭智能化功能
RESUME, 4.13(2).DB	! 从子模型文件4.13(2).DB中恢复子模型数据
/NOPR	! 关闭提示反馈
/INPUT, , cbdo, , :cb1	! 从文件4.13(2).CBDO中读入插值边界信息
/GOPR	! 打开提示反馈
LSEL, S, LINE, , 3	! 选择3号线, 以应用其余的边界条件
DL, ALL, , SYMM	! 对选择集中的线应用对称边界条件
FINISH	! 退出前处理模块PREP7
! (14) 进入求解模块求解
/SOLU	! 加载求解模块
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLU
! (15) 在后处理模块显示应力计算结果, 
/POST1	! 进入一般后处理模块
SET, 1, 1	! 读入第1个荷载步的计算结果
NSORT, S, X, , , 3	! 将计算结果按SX排序, 只保留前3项最高的结果
PRNSOL, S, COMP	! 列表显示应力分量显示
/WINDOW, 1, OFF	! 关闭1号窗口的显示
/AUTO, 3	! 指定3号窗口为以最佳比例显示图形
/WINDOW, 3, BOT	! 在图形窗口的下半部分定义窗口3
/NOERASE	! 不擦除图形窗口
/PLOPTS, MINM, 1	! 显示最大和最小图例
/USER, 3	! 指定用户自定义的视点距离和焦点位置
/DIST, 3, 5.08	! 指定用户自定义的视点距离为5.08mm
/FOCUS, 3, 5.08, 12.7	! 指定用户自定义的焦点位置(5.08mm, 12.7mm)
/CONTOUR, 3, , AUTO	! 给窗口3应用自动梯度比例
PLNSOL, SX	! 列表显示应力分量SX计算结果
*GET, SUBSTR, NODE, 18, S, X	! 提取18号节点的X方向应力SX, 存入变量SUBSTR
! (16) 将计算结果保存到文件TABLE_2
SAVE, TABLE_2	! 保存数据到默认数据库TABLE_2
! (17) 改变单元类型为PLANE146, 重新计算
FINISH	! 退出后处理模块
/CLEAR, NOSTART	! 清除所有数据, 不加载初始化文件
/FILNAME, 4.13(3)	! 定义子结构分析文件名
/PREP7	! 加载前处理模块
SMRT, OFF	! 关闭智能化网格功能
/TITLE, EX4.13(3), STRESS CONCENTRATION AT A HOLE IN A PLATE With PLANE146 Element.
/NOPR	! 关闭输入信息反馈提示
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, PLANE146	! 定义第1类单元为二维结构固体自适应单元PLANE146
MP, EX, 1, 207E3	! 定义第1类材料的弹性模量EX ＝ 207 GPa
MP, NUXY, 1, 0.3	! 定义第1类材料的泊松比NUXY ＝ 0.3
! (18)定义关键点
K, 1, 152.4	! 在坐标(152.4mm, 0)定义第1个关键点
K, 2, 152.4, 152.4	! 在坐标(152.4mm, 152.4mm)定义第2个关键点
K, 3, , 152.4	! 在坐标(0, 152.4mm)定义第3个关键点
K, 4, , 12.7	! 在坐标(0, 12.7mm)定义第4个关键点
K, 5, 12.7	! 在坐标(12.7mm, 0)定义第5个关键点
K, 6	! 在坐标(0, 0)定义第6个关键点
! (19)定义线和单元剖分的段数
L, 1, 2	! 过关键点1, 2定义直线
L, 2, 3	! 过关键点2, 3定义直线
L, 3, 4	! 过关键点3, 4定义直线
LESIZE, 3, , , 4, .25	! 给3号线定义剖分段数为4和间距比例0.25
LARC, 4, 5, 6, 12.7	! 以6号点为圆心, 过5号点和6号点, 以12.7mm为半径定义圆弧
LESIZE, 4, , , 6	! 给4号线定义剖分段数为6
L, 5, 1	! 过关键点5, 1定义直线
LESIZE, 5, , , 4, 4	! 给3号线定义剖分段数为4和间距比例为4
! (20) 定义面并将其剖分为面单元, 图形显示单元布置
AL, 1, 2, 3, 4, 5	! 用线1, 2, 3, 4和5定义面
ESIZE, , 4  	! 单元剖分尺寸为分为4段
AMESH, ALL 	! 将所有的面剖分为面单元
/AUTO, 1	! 设置模型显示的最佳比例
/PLOPTS, INFO, 0	! 关闭图形窗口中图例信息显示栏
/PLOPTS, WINS, 0	! 关闭图形窗口的自适应调整
/WINDOW, , LTOP	! 在图形窗口的左上角定义窗口1
EPLOT	! 图形显示单元布置
! (21) 选择相应的线, 施加位移和荷载边界条件
LSEL, S, LINE, , 3, 5, 2	! 选择3号线和5号线
DL, ALL, , SYMM	! 对所有的选择集中的线施加对称位移约束
LSEL, S, LINE, , 1	! 选择1号线
NSLL, , 1	! 选择1号线上的节点
SF, ALL, PRES, -6.895	! 施加板端的拉力荷载6.895 MPa
LSEL, ALL	! 选择所有的线
NSEL, ALL	! 选择所有节点
CSYS, 1	! 将坐标系切换到柱坐标系
FINISH	! 退出前处理模块PREP7
! (22) 进入求解模块求解
/SOLU	! 加载求解模块
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLU
! (23) 进入后处理, 显示计算结果
/POST1	! 进入一般后处理模块
SET, 1, 1	! 读入第1个荷载步的计算结果
/WINDOW, 1, OFF	! 关闭1号窗口的显示
/NOERASE	! 不擦除图形窗口
/DSCALE, 2, 1	! 用缩放位移约束命令, 缩放位移约束实部2倍, 虚部1倍
/WINDOW, 2, RTOP	! 在图形窗口的右上角定义窗口2
PLNSOL, S, X	! 图形应力SX的分布图
*GET, SUBSTR, NODE, 18, S, X	! 提取18号节点的应力SX, 存入变量SUBSTR
! (24) 将结果数据保存到数据库TABLE_3
SAVE, TABLE_3	! 保存数据到默认数据库TABLE_3
FINISH	! 退出后处理模块
! (25) 从各个数据库中恢复数据, 依次输出显示
RESUME, TABLE_1
*STATUS, PARM	! 列表显示变量参数的值
RESUME, TABLE_2
*STATUS, PARM	! 列表显示变量参数的值
RESUME, TABLE_3
*STATUS, PARM	! 列表显示变量参数的值
FINISH	! 退出后处理模块
! (26) 删除中间数据文件
/DEL, TABLE_1
/DEL, TABLE_2
/DEL, TABLE_3