ex5.1.inp

ansys有限元分析程序

! 5.1  含椭圆孔的椭圆薄膜在外部张力作用下的静力分析
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除数据, 不加载初始化文件
! （1）设置工程选项, 分析类型、单元类型和材料参数
/FILNAME, EX5.1(1)	! 设置工程文件名称
/PREP7	! 加载前处理模块
SMRT, OFF	! 关闭智能化网格功能
/TITLE, EX5.1(1), ELLIPTIC MEMBRANE UNDER A UNIFORMLY-DISTRIBUTED LOAD(PLANE42).
/COM,   USING 2-D STRUCTURAL SOLID, PLANE42
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, PLANE42, , , 3	! 定义第1类单元为带厚度的平面应力单元PLANE42
MP, EX, 1, 210E9	! 定义第1类材料的弹性模量EX＝210GPa
MP, NUXY, 1, 0.3	! 定义第1类材料的泊松比NUXY=0.3
R, 1, 0.1  	! 设置厚度0.1m
! （2）在椭圆坐标系下, 定义关键点、线和面
! LOCAL, KCN, KCS, XC, YC, ZC, THXY, THYZ, THZX, PAR1, PAR2 
LOCAL, 11, 1, , , , , , , 0.5	! 在坐标原点定义椭圆坐标系, Y方向的短轴与X方向的长轴之比为0.5
K, 1, 2, 90	! 在坐标(2m, 90度)处定义第1个关键点
K, 2, 2, 0	! 在坐标(2m, 0)处定义关键点2
L, 1, 2	! 过关键点1, 2定义直线
LOCAL, 12, 1, , , , , , , 0.8461585	! 在坐标原点定义椭圆坐标系, Y方向的短轴与X方向的长轴之比为0.8461585
K, 3, 3.25, 90	! 在坐标(3.25m, 90度)处定义关键点3
K, 4, 3.25, 0	! 在坐标(3.25m, 0)处定义关键点4
L, 3, 4	! 过关键点3, 4定义直线
CSYS, 0	! 将当前坐标系切换到整体直角坐标系
L, 2, 4	! 过关键点2, 4定义直线
L, 1, 3	! 过关键点1, 3定义直线
AL, 2, 4, 1, 3	! 用线2, 4, 1和3定义面
! (3)在相应的线上施加位移约束
DL, 4, 1, SYMM	! 对1号面的4号线施加对称位移边界条件
DL, 3, 1, SYMM	! 对1号面的3号线施加对称位移边界条件
SFL, 2, PRES, -10E6	! 对2号线施加线压力荷载10E6 N/m
MSHK, 2   	! 如果可能, 采用映射网格剖分
MSHA, 0, 2D  	! 使用四边形网格
SAVE	! 保存模型数据到默认数据库
FINISH	! 退出前处理模块
! (4)执行自适应网格剖分和求解宏程序, 执行自定义的后处理宏程序
ADAPT, 3, 7, , , 1	! 使用ANSYS预定义的自适应网格剖分宏程序, 用NSOLN=3, STARGT=7, 和FACMX=1等参数来求解
/Post1	! 进入一般后处理模块
EPLOT	! 绘制单元
PLDISP, 1	! 显示变形图
PLNSOL, S, X	! 显示水平应力分布
PLNSOL, S, Y	! 显示竖向应力分布
PLNSOL, S, 1	! 显示最大拉应力分布
FINISH	! 退出后处理模块
! (5)定义完成后处理功能的宏程序MAC
*CREATE, MAC	! 定义用于完成后处理的宏程序MAC
/Graphics, Full	! 关闭PowerGraphics, 以显示能量范数误差
/POST1    	! 进入一般后处理模块
PRERR      	! 打印能量范数误差(SEPC)
NSEL, S, LOC, Y, 0.0	! 选择Y=0位置的节点
NSEL, R, LOC, X, 2.0	! 继续选择X=2.0m位置的节点
*GET, MNODE, NODE, , NUM, MAX	! 提取当前选择集中最大节点号, 保存到变量MNODE
*GET, SY_D, NODE, MNODE, S, Y	! 提取节点MNODE的竖向应力SY, 存入变量SY_D
NSEL, ALL	! 选择所有节点
ESEL, ALL	! 选择所有单元
*STATUS, PARM	! 显示内存变量的值
/Graphics, Power	! 打开PowerGraphics, 以显示其他结果
*END             	! 宏程序MAC定义结束
! (6)用前定义的宏程序完成后处理, 将计算结果保存到变量
*USE, MAC   	! 执行后处理宏程序MAC
SAVE, TABLE1	! 保存数据到默认数据库TABLE1
FINISH	! 退出后处理模块
! (7)改变单元为PLANE82, 重新计算
/CLEAR, NOSTART	! 在开始第二部分求解前, 清除数据, 不加载初始化文件! 
/FILNAME, EX5.1(2)	! 设置工程文件名称
/PREP7	! 加载前处理模块
SMRT, OFF	! 关闭智能化网格功能
RESUME, EX5.1(1).db	! 重新读入模型数据
/TITLE, EX5.1(2), ELLIPTIC MEMBRANE UNDER A UNIFORMLY-DISTRIBUTED LOAD(PLANE82).
ET, 1, PLANE82, , , 3	! 定义第1类单元为带厚度的二维8节点平面应力单元PLANE82
FINISH	! 退出前处理模块PREP7
! （8）执行自适应网格宏程序, 执行自定义的后处理宏程序
ADAPT, 2, 5, , , 1	! 使用ANSYS预定义的自适应网格剖分宏程序, 用NSOLN=2, STARGT=5, 和FACMX=1等参数来求解
/Post1	! 进入一般后处理模块
/WINDOW, 5, OFF	! 关闭5号窗口的显示
/WINDOW, 1, FULL	! 将1号窗口的显示设置为整个图形窗口
PLNSOL, S, EQV	! 显示等效MISES应力分布
PLDISP, 1	! 显示变形图
PLNSOL, S, X	! 显示水平应力分布
PLNSOL, S, Y	! 显示竖向应力分布
PLNSOL, S, 1	! 显示最大拉应力分布
FINISH	! 退出后处理模块
! （9）用前定义的宏程序完成后处理, 将计算结果保存到变量
*USE, MAC 	! 执行后处理宏程序MAC
FINISH	! 退出后处理模块