我的命令流，大家看有没有什么问题


!********参数准备********
!横梁参数
l1=14.5 $h1=0.8 $b1=0.6
!直桩参数
l2=6 $h2=0.5 $b2=0.5
rat=4 !叉桩斜率
gap1=0.7 !双直桩中心距
gap2=0.85 !叉桩顶部中心距
ee=25e9 !材料弹性模量
u=0.17 !泊松比
g=ee/2/(1+u) !剪切模量
density=2500 !材料密度
fs=1.2 !矩形截面的剪切方向系数
elth=0.3 !单元划分长度
nfreq=10 !模态分析的阶数
nof=3 !节点自由度数

!********建模********
/PREP7
!加载单元类型、材料属性和实常数
et,1,beam3
mp,ex,1,ee
mp,dens,1,density
mp,prxy,1,u
mp,gxy,1,g
r,1,b1*h1,b1*h1**3/12,h1,fs !横梁断面
r,2,b2*h2,b2*h2**3/12,h2,fs !桩断面
!建立几何模型
k,1,0 $k,2,l1 !横梁端点
k,3,2-gap1/2 $k,4,2+gap1/2 !双直桩顶点
k,5,12.5-gap2/2 $k,6,12.5+gap2/2 !叉桩顶点
k,7,2-gap1/2,-l2 $k,8,2+gap1/2,-l2 !双直桩底部
k,9,12.5-gap2/2-l2/rat,-l2 $k,10,12.5+gap2/2+l2/rat,-l2 !叉桩底部
k,11,7.25 $k,12,7.25,-l2 !单直桩上下端点（单桩位于排架中点）

l,1,2
l,3,7 $l,4,8
l,5,9 $l,6,10
l,11,12

lglue,all !连接所有相连的线为一个整体

!将几何模型转换为有限元模型
nummrg,all
numcmp,all

lsum=lsinqr(0,12) !获取定义的线数目

lesize,all,elth

lsel,s,line,,6,lsum !对横梁划分单元
latt,1,1,1
lmesh,all

lsel,s,line,,1,5 !对桩划分单元
latt,1,2,1
lmesh,all

allsel

nummrg,all
numcmp,all

finish
!********求解模态********
/solu
antype,modal,new
modopt,subsp,nfreq
mxpand,nfreq,,,1
nsel,s,loc,y,-l2
d,all,all,0
allsel
solve
finish

急用，拜托了！