APDL梁单元临界转速计算命令流

*DIM, L, ARRAY, 1000 ! 定义长度为1000的数组L
*DO, i, 1, 1000, 1 ! 循环1000次，每次步长为1
L(i) = (i-1)/1000 ! 计算数组L中每个元素的值
*ENDDO ! 结束循环

*DIM, E, ARRAY, 1000 ! 定义长度为1000的数组E
*DO, i, 1, 1000, 1 ! 循环1000次，每次步长为1
E(i) = L(i)/100 ! 计算数组E中每个元素的值
*ENDDO ! 结束循环

ET, 1, LINK180 ! 定义LINK180单元类型
MP, DENS, 1, 7850 ! 定义材料密度
MP, E, 1, 2.1e11 ! 定义材料弹性模量
MP, NU, 1, 0.3 ! 定义材料泊松比

L, 1, 0.5 ! 定义梁的长度为0.5m
A, 1, 0.001 ! 定义梁的截面积为0.001m^2
FINISH ! 结束几何定义

K, 1, 1, 1 ! 定义梁的节点1的约束为固定
K, 2, 2, 2 ! 定义梁的节点2的约束为固定

*DIM, W, ARRAY, 1000 ! 定义长度为1000的数组W
*DO, i, 1, 1000, 1 ! 循环1000次，每次步长为1
W(i) = i*100 ! 计算数组W中每个元素的值
*ENDDO ! 结束循环

*DIM, F, ARRAY, 1000 ! 定义长度为1000的数组F
*DO, i, 1, 1000, 1 ! 循环1000次，每次步长为1
F(i) = W(i)*L(500)/2 ! 计算数组F中每个元素的值
*ENDDO ! 结束循环

*DIM, CRIT, ARRAY, 1000 ! 定义长度为1000的数组CRIT
*DO, i, 1, 1000, 1 ! 循环1000次，每次步长为1
NSOL, 1 ! 求解节点位移和反力
F, 2, F(i) ! 在节点2施加力F(i)
CMSL, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0 ! 计算梁的传递矩阵
SFE, 1, 1, CRIT(i) ! 计算梁的临界转速
*ENDDO ! 结束循环

*CFOPEN, CRIT.TXT ! 打开输出文件
*VWRITE, CRIT(1), ',', CRIT(2), ',', CRIT(3), ',', CRIT(4), ',', CRIT(5) ! 输出前五个临界转速
*CFWRITE ! 写入文件
*CFCLOSE ! 关闭文件