ex4.10(1).inp

ansys有限元分析程序

! 4.10  在端部集中力下的变截面悬臂梁
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH	! 退出以前模块
/CLEAR, NOSTART	! 清除数据，不加载初始化文件
! (1)设置工程选项，定义分析类型、单元类型和材料参数
/FILNAME, EX4.10(1)	! 定义工程文件名称
/PREP7	! 加载前处理模块
/TITLE, EX4.10(1), LATERALLY LOADED TAPERED SUPPORT STRUCTURE (QUAD. ELEMENTS)
! 使用 PLANE42 单元
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, PLANE42, , , 3, , , 2	! 定义第1类单元为带厚度（参数“3”表示）的4节点平面应力单元PLANE42，打开表面结果输出（参数“2”表示），以提取中间节点的结果
R, 1, 50.8	! 定义厚度为50.8mm
MP, EX, 1, 207E3	! 定义第1类材料的弹性模量为EX＝207E3 N/mm^2
MP, NUXY, 1, 0.0	! 定义第1类材料的泊松比NUXY＝0.0
! (2)定义节点位置和单元连接
N, 1, 635	! 在(635mm, 0, 0)位置定义第1号节点
N, 7, 1905	! 在(1905mm, 0, 0)位置定义第7号节点
FILL	! 在1号节点和7号节点之间，填充生成2号到6号节点
N, 8, 635, -76.2	! 在(635mm, -76.2mm, 0)位置定义第8号节点
N, 14, 1905, -228.6	! 在(1905mm, -228.6mm, 0)位置定义第14号节点
FILL	! 在8号节点和14号节点之间，填充生成9号到13号节点
E, 2, 1, 8, 9	! 过节点2，1，8和9定义单元
EGEN, 6, 1, 1	! 按照前面的单元模式，循环6次定义单元，每次节点号增加1
! (3)选择节点，施加位移边界条件和荷载
NSEL, S, LOC, X, 1905	! 选择X＝1905mm位置的节点
D, ALL, ALL  	! 约束选择集中的所有节点的所有自由度
NSEL, ALL	! 选择所有节点
F, 1, FY, -17793	! 在1号节点的-Y方向施加17793N的集中力
FINISH	! 退出前处理模块PREP7
! (4)在求解模块，开始求解
/SOLU	! 加载求解模块
OUTPR, , 1	! 只输出第1荷载步的基本信息
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLU
! (5)提取指定位置节点的应力，绘制变形图
/POST1	! 进入一般后处理模块
END_NODE = NODE (1905, 0, 0)	! 选择距离（1905mm, 0, 0）位置最近的节点，命名为END_NODE
*GET, STS_E_42, NODE, END_NODE, S, X	! 提取节点END_NODE的X方向应力SX，存入变量TS_E_42
PLDISP, 2	! 带原来模型图轮廓，显示变形后的模型
MID_NODE = NODE (1270, 0, 0)	! 选择距离（1270mm, 0, 0）位置最近的节点，命名为MID_NODE
*GET, STS_M_42, NODE, MID_NODE, S, EQV	! 提取节点MID_NODE位置的MISES应力，存入STS_M_42
PLNSOL, S, X	! 图形显示水平方向的正应力分布
PLNSOL, S, 1	! 图形显示拉应力分布
PLNSOL, S, 3	! 图形显示压应力分布
FINISH	! 退出后处理模块
! (6)改用更高阶的平面应力单元PLANE82重新计算
/PREP7	! 加载前处理模块
/TITLE, EX4.10(2), LATERALLY LOADED TAPERED SUPPORT STRUCTURE (PLANE82).
ET, 1, PLANE82, , , 3, , , 2	! 定义第1类单元为带厚度（参数“3”表示）的8节点平面应力单元，打开表面结果输出（参数“2”表示），以提取中间节点的结果
EMID  	! 给PLANE82单元添加中间节点
NSEL, R, LOC, X, 1905	! 选择X＝1905mm位置的节点
NSEL, R, LOC, Y, -114.3	! 选择固定端处的中点
D, ALL, ALL   	! 位移约束固定端的中点
NSEL, ALL	! 选择所有节点
FINISH	! 退出前处理模块PREP7
! (7)在求解模块求解
/SOLU	! 加载求解模块
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLU
! (8)提取指定位置的节点的应力计算结果
/POST1	! 进入一般后处理模块
END_NODE = NODE (1905, 0, 0)	! 选择距离（1905mm, 0, 0）最近的节点，命名为END_NODE
*GET, STS_E_82, NODE, END_NODE, S, X	! 提取节点END_NODE的X方向应力，存入TS_E_82
MID_NODE = NODE (1270, 0, 0)	! 选择距离（1270mm, 0, 0）最近的节点，命名为MID_NODE
*GET, STS_M_82, NODE, MID_NODE, S, EQV	! 提取节点MID_NODE位置的MISES应力，存入STS_M_82
PLNSOL, S, X	! 图形显示水平方向的正应力分布
PLNSOL, S, 1	! 图形显示拉应力分布
PLNSOL, S, 3	! 图形显示压应力分布
FINISH	! 退出后处理模块