! 6.1 利用梁单元计算压杆稳定性
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART ! 清除数据, 不加载初始化文件
! （1）定义工程选项、单元类型、实常数和材料参数
/FILNAME, EX6.1 ! 定义工程文件名称
/PREP7 ! 加载前处理模块
/TITLE, EX6.1, BUCKLING OF A BAR WITH HINGED END BY BEAM3.
ET, 1, BEAM3 ! 定义第1类单元为二维弹性梁单元BEAM3
R, 1, 161.29, 2167.87, 12.7 ! 定义梁的截面积为161.29mm^2, 惯性矩为2167.87mm^4, 截面高度为12.7mm
MP, EX, 1, 207000 ! 定义第1类材料的弹性模量EX = 207000 N/mm2
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 ! 加载求解模块
ANTYPE, STATIC ! 定义分析类型为静力分析
PSTRES, ON ! 考虑预应力效应
D, 1, ALL ! 约束1号节点的所有自由度
F, 11, FY, -1 ! 给11号节点施加沿着－Y方向的单位力
OUTPR, , 1 ! 设置只输出第1荷载步的基本计算结果
SOLVE ! 开始静力分析求解
FINISH ! 退出求解模块
! (4)再次进入求解模块, 进行屈曲分析
/SOLU ! 加载求解模块
ANTYPE, BUCKLE ! 指定分析类型为屈曲分析
BUCOPT, LANB, 1 ! 使用分块LANCZOS特征值展开方法, 只展开第1阶失稳模态
MXPAND, 1 ! 展开第1阶失稳模态
SOLVE ! 开始屈曲分析的求解
! (5)提取计算结果
*GET, FCR, MODE, 1, FREQ ! 提取第1阶频率计算结果, 并保存到变量FCR
*STATUS, PARM ! 显示内存变量的值
FINISH
