ex6.3.inp

ansys有限元分析程序

! 6.3  悬臂压杆的过曲屈分析
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除数据, 不加载初始化文件
! (1)设置工程选项、单元类型、实常数和材料参数
/FILNAME, EX6.3	! 定义工程文件名称
/PREP7	! 加载前处理模块
/TITLE, EX6.3, LARGE DEFLECTIONS OF A BUCKLED BAR (THE ELASTICA).
ANTYPE, STATIC	! 定义分析类型为静力分析
NLGEOM, ON  	! 激活大变形选项
ET, 1, BEAM3, , , , , , , 1	! 定义第1类单元为二维弹性梁单元BEAM3, 关闭输出选项
R, 1, 161.29, 2167.86, 12.7	! 定义梁的截面积为161.29mm^2, 惯性矩为2167.86mm^4, 截面高度为12.7mm
MP, EX, 1, 207000	! 定义第1类材料的弹性模量EX = 207000 N/mm^2
MP, PRXY, , 0.3	! 定义第1类材料的泊松比PRXY = 0.3
! （2）定义节点和单元
N, 1	! 在坐标原点(0, 0, 0)定义第1号结点
N, 11, , 2540	! 在坐标(0, 2540mm, 0)定义第11号结点
FILL	! 在第1号节点和第11号节点之间填充生成其余节点
E, 1, 2	! 过节点1和节点2定义单元
EGEN, 10, 1, 1	! 按照前面单元的模式循环10次, 生成其余的9个单元
FINISH	! 退出前处理模块PREP7
! （3）设置求解选项, 位移约束和荷载, 开始求解
/SOLU	! 加载求解模块
SOLCONTROL, 0	! 指定不使用（用参数0来表示）最优化的非线性求解器
NEQIT, 250	! 指定最大的非线性叠代次数为150次
OUTPR, BASIC, LAST	! 只输出最后一个荷载步的基本内容
D, 1, ALL	! 约束1号节点的所有自由度
FCR=-171.5	! 设置压杆失稳的临界荷载值171.5N
PI=3.14159265359	! 定义参数Pi
F, 11, FY, FCR*1.015	! 给11号节点施加荷载
F, 11, FX, 1 	! 施加非常小的横向荷载1N
SOLVE	! 开始求解荷载步1
! (4)增加荷载, 继续求解
FDEL, 11, FX	! 移去前面施加的横向荷载
F, 11, FY, FCR*1.063	! 给11号节点施加更大的轴向荷载
SOLVE	! 开始求解荷载步2
F, 11, FY, FCR*1.152	! 给11号节点施加更大的轴向荷载
SOLVE	! 开始求解荷载步3
F, 11, FY, FCR*1.293	! 给11号节点施加更大的轴向荷载
SOLVE	! 开始求解荷载步4
F, 11, FY, FCR*1.518	! 给11号节点施加更大的轴向荷载
SOLVE	! 开始求解荷载步5
F, 11, FY, FCR*1.884	! 给11号节点施加更大的轴向荷载
SOLVE	! 开始求解荷载步6
FINISH	! 退出求解模块SOLUTION
! (5)在一般后处理模块中, 绘制后屈曲构型图
/POST1	! 进入一般后处理模块
/USER	! 定义用户自定义的焦点和观察距离
/FOCUS, , 1270, 1270	! 所有图形对象居中
/DIST, , 1397  	! 设置观察距离
/DSCALE, , 1	! 按照模型的实际比例
SET, 1, 0	! 读入第1个荷载步的计算结果
PLDISP, 1	! 带原来模型图, 显示变形后的模型
/NOERASE	! 不清除屏幕, 重叠绘制下一个图形
SET, 2, 0	! 读入第2个荷载步的计算结果
PLDISP	! 显示变形后的模型
SET, 3, 0	! 读入第3个荷载步的计算结果
PLDISP	! 显示变形后的模型
*GET, UX3, NODE, 11, U, X	! 提取11号节点的X方向的位移UX, 存入变量UX3
*GET, UY3, NODE, 11, U, Y	! 提取11号节点的Y方向的位移UY, 存入变量UY3
*GET, ROT3A, NODE, 11, ROT, Z	! 提取11号节点绕Z轴的转角ROTZ, 存入变量ROT3A
SET, 4, 0	! 读入第4个荷载步的计算结果
PLDISP	! 显示变形后的模型
*GET, UX4, NODE, 11, U, X	! 提取11号节点的X方向的位移UX, 存入变量UX4
*GET, UY4, NODE, 11, U, Y	! 提取11号节点的Y方向的位移UY, 存入变量UY4
*GET, ROT4A, NODE, 11, ROT, Z	! 提取11号节点绕Z轴的转角ROTZ, 存入变量ROT4A
SET, 5, 0	! 读入第5个荷载步的计算结果
PLDISP	! 显示变形后的模型
*GET, UX5, NODE, 11, U, X	! 提取11号节点的X方向的位移UX, 存入变量UX5
*GET, UY5, NODE, 11, U, Y	! 提取11号节点的Y方向的位移UY, 存入变量UY5
*GET, ROT5A, NODE, 11, ROT, Z	! 提取11号节点绕Z轴的转角ROTZ, 存入变量ROT5A
SET, 6, 0	! 读入第6个荷载步的计算结果
PLDISP	! 显示变形后的模型
*GET, UX6, NODE, 11, U, X	! 提取11号节点的X方向的位移UX, 存入变量UX6
*GET, UY6, NODE, 11, U, Y	! 提取11号节点的Y方向的位移UY, 存入变量UY6
*GET, ROT6A, NODE, 11, ROT, Z	! 提取11号节点绕Z轴的转角ROTZ, 存入变量ROT6A
*STATUS, PARM	! 显示内存变量的值
FINISH