ex6.2.inp

ansys有限元分析程序

! 6.2  利用实体单元计算压杆稳定性
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除数据, 不加载初始化文件
! (1)设置工程选项、 单元类型、实常数和材料参数
/FILNAME, EX6.2	! 定义工程文件名称
/PREP7	! 加载前处理模块
/TITLE, EX6.2, BUCKLING OF A BAR WITH HINGED BY PLANE42.
ET, 1, PLANE42, , , 3	! 定义第1类单元为二维固体单元PLANE42
R, 1, 12.7 	! 定义方板厚度为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号节点之间填充生成其余节点
! NGEN, ITIME, INC, NODE1, NODE2, NINC, DX, DY, DZ, SPACE
NGEN, 2, 11, 1, 11, , 12.7	! 循环2次, 节点号增量11, 按1号节点到11号节点的范围, X坐标增加12.7mm。生成12号到22号节点。
E, 12, 13, 2, 1	! 过节点12, 13, 2和1定义单元
EGEN, 10, 1, 1	! 循环10次, 每次单元的节点号增加1
FINISH	! 退出前处理模块PREP7
! （3）定义求解选项, 位移约束和荷载力, 并开始静力分析的求解
/SOLU	! 加载求解模块
PSTRES, ON  	! 考虑预应力效应
D, 1, ALL, , , 12, 11	! 位移约束所有自由度
F, 11, FY, -0.5, , 22, 11	! 将单位荷载平均分配到自由端的节点上
OUTPR, , 1	! 设置只输出第1荷载步的基本计算结果
SOLVE	! 开始静力分析的求解
FINISH	! 退出求解模块SOLUTION
! (4)再次进入求解模块, 进行屈曲分析
/SOLU	! 加载求解模块
ANTYPE, BUCKLE	! 定义分析类型为屈曲分析
BUCOPT, SUBSP, 1	! 使用子空间求解方法, 展开1阶失稳模态
MXPAND, 1  	! 展开第1阶失稳模态
SOLVE	! 开始特征值问题的求解
! (5)提取计算结果
*GET, FCR1, MODE, 1, FREQ	! 提取第1阶频率计算结果, 并保存到变量FCR1
FINISH
! (6)改用其他特征值求解方法
/SOLU	! 进入求解模块
BUCOPT, LANB , 1	! 使用分块LANCZOS求解方法, 展开1阶失稳模态
SOLVE	! 开始模态求解
! (7)显示计算结果
*GET, FCR2, MODE, 1, FREQ	! 提取第1阶频率计算结果, 并保存到变量FCR2
*STATUS, PARM	! 显示内存变量的值
FINISH