! 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
