ex4.8.inp

ansys有限元分析程序

! 4.8  用实体单元分析变截面杆的拉伸
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH	! 退出以前模块
/CLEAR, NOSTART	! 清除数据，不加载初始化文件
! (1)设置工程选项, 定义分析类型、单元类型和材料参数
/FILNAME, EX4.8	! 定义工程文件名称
/PREP7	! 加载前处理模块
SMRT, OFF	! 关闭智能化功能
/TITLE, EX4.8, ELONGATION OF A SOLID BAR WITH VARIED SECTION.
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, SOLID45	! 定义第1类单元为三维结构固体单元SOLID45
MP, EX, 1, 70.71E3	! 定义第1类材料的弹性模量EX=70.71 GPa
MP, NUXY, 1, 0.3	! 定义第1类材料的泊松比NUXY=0.3
! (2)定义固定端4个点和自由端4个点，过这8个点定义体，对体剖分单元
/VIEW, 1, 1, 1, 1	! 改变视点为(1, 1, 1), 沿着等轴测方向
K, 1, 25.4, , 25.4	! 在坐标（25.4mm，0，25.4mm）处定义第1个关键点
K, 2, -25.4, , 25.4	! 在坐标（-25.4mm，0，25.4mm）处定义第2个关键点
K, 3, -25.4, , -25.4	! 在坐标（-25.4mm，0，-25.4mm）处定义第3个关键点
K, 4, 25.4, , -25.4	! 在坐标（25.4mm，0，-25.4mm）处定义第4个关键点
K, 5, 12.7, 254, 12.7	! 在坐标（12.7mm，254mm，12.7mm）处定义第5个关键点
K, 6, -12.7, 254, 12.7	! 在坐标（-12.7mm，254mm，12.7mm）处定义第6个关键点
K, 7, -12.7, 254, -12.7	! 在坐标（-12.7mm，254mm，-12.7mm）处定义第7个关键点
K, 8, 12.7, 254, -12.7	! 在坐标（12.7mm，254mm，-12.7mm）处定义第8个关键点
V, 1, 2, 3, 4, 5, 6, 7, 8	! 过前面定义的8个关键点定义体
LSEL, S, LINE, , 5, 11, 2	! 从5号线到11号线，增量2选择5，7，9和11号线
LESIZE, ALL, , , 7	! 指定前面选中的线在剖分单元是剖分为7段
LSEL, ALL 	! 选择所有的线
ESIZE, , 1 	! 单元剖分段数为1
/OUT, SCRATCH	! 将输出结果重定向到SCRATCH文件，以避免屏幕输出
VMESH, 1	! 将1号体剖分为体单元
/OUT	! 结束数据重定向，关闭输出文件
! (3)施加位移约束和荷载
OUTPR, BASIC, ALL	! 设置输出所有荷载步的基本信息
NSEL, S, LOC, Y, 0	! 选择Y=0位置的所有节点，即固定端位置的节点
D, ALL, ALL 	! 位移约束选择集中的所有节点的所有自由度
NSEL, ALL	! 选择所有节点
NSEL, S, LOC, Y, 254	! 选择Y=254mm位置的所有节点，即自由端节点
SF, , PRES, -68.949	! 给节点指定表面荷载44483N/25.4mm^2=68.949MPa
NSEL, ALL	! 选择所有节点
FINISH	! 退出前处理模块PREP7
! (4)在求解模块，开始静力分析的求解
/SOLU	! 加载求解模块
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLUTION
! (5)在后处理模块中，显示Y方向应力云图
/POST1	! 进入一般后处理模块
ETABLE, SIGY, S, Y	! 以Y方向的正应力为内容，定义单元表SIGY
/VIEW, 1, 1	! 改变视点为(1, 0, 0), 沿着X轴的负方向
/VUP, 1, -Y	! 设置-Y轴向上
NSLE, S    	! 选择和单元相连接的节点
PLNSOL, S, Y  	! 显示轴向应力
ESEL, S, ELEM, , 4	! 选择中间的第4号单元
PRETAB, SIGY 	! 显示单元表中的单元应力计算结果
PRNSOL, S, COMP	! 显示节点应力计算结果
ESEL, ALL  	! 选择所有单元
NSEL, S, LOC, Y, 254	! 选择位于Y＝254mm的节点（自由端）
PRNSOL, DOF  	! 显示节点位移计算结果
NSEL, ALL	! 选择所有节点
! (6)选择指定位置的节点和单元，并提取它们的位移和应力
/NOPR	! 关闭输入反馈提示信息
MID_NODE = NODE(0, 128, 0)	! 选择距离（0, 128mm, 0）最近的节点，命名为MID_NODE
MID_ELM = ENEARN(MID_NODE)	! 选择距离MID_NODE节点最近的单元，命名为MID_ELM
BOT_NODE = NODE (0, 254, 0)	! 选择距离（0, 254mm, 0）最近的节点，命名为BOT_NODE
*GET, DEF, NODE, BOT_NODE, U, Y	! 提取BOT_NODE节点的位移UY，存入变量DEF
*GET, STRSS, ELEM, MID_ELM, ETAB, SIGY	! 在单元表SIGY中提取MID_ELM单元的应力，存入变量STRSS
*STATUS, PARM	! 显示所有内存变量的值
FINISH	! 退出后处理模块