求自振频率怎么求啊。（附命令流）

zhshle1010

我的模型中间部分是两块材料不同的长方体，上部是一根梁结构，通过竖向弹簧连接。下部也是竖向弹簧固定弹簧最下端的节点。我求了好多次实在没求出来，请各位看下（对于梁和长方体的两端边界条件暂时固定吧）命令流如下：请问错在哪里呢？？？？？？？？？？？？？？？

/prep7
et,1,beam4
et,2,solid45
et,3,combin14,,2
!
mp,ex,1,2.06e11                                                                        
mp,prxy,1,0.3                                                                           
mp,dens,1,7800
r,1,7745e-6,3217e-8,                                                                     
!
r,2,25e6                                                                             
!**
mp,ex,2,3.4e10                                                                          
mp,prxy,2,0.2                                                                          
mp,dens,2,2.8e3                                                                        
!**
mp,ex,3,8e9                                                                             
mp,prxy,3,0.2                                                                           
mp,dens,3,2.6e3                                                                       
!
r,4,80e6                                                                 


!
wpoffs,,,-1.4
blc4,,,0.625*40,0.24,2.8
wpoffs,,-0.3,-0.3
blc4,,,0.625*40,0.3,3.4
wpcsys
vglue,all
wpoffs,,,1.435/2
vsbw,all
wpoffs,,,(2.8-1.435)/2
vsbw,all
wpoffs,,,-(2.8+1.435)/2
vsbw,all
wpoffs,,,-(2.8-1.435)/2
vsbw,all
wpcsys
vsel,s,loc,y,0,0.26
vatt,2,,2
vsel,s,loc,y,-0.3,0
vatt,3,,2
allsel
lesize,31,,,40*2
esize,0.3
vmesn,all
nsel,s,loc,y,0.24
nsel,r,loc,z,1.435/2
ngen,2,4000,all,,,,0.2
allsel
type,3
real,2
*do,i,1,81
e,1134+i,1134+i+4000
*enddo

nsel,s,loc,y,0.24
nsel,r,loc,z,-1.435/2
ngen,2,6000,all,,,,0.2
allsel
type,3
real,2
*do,i,1,81
e,2025+i,2025+i+6000
*enddo

type,1
real,1
mat,1
*do,i,1,80
e,1134+i+4000,1135+i+4000
*enddo
*do,i,1,80
e,2025+i+6000,2026+i+6000
*enddo

nsel,s,loc,y,-0.3
ngen,2,9000,all,,,,-0.2
type,3
real,4
*do,i,2,82
e,i,i+9000
*enddo
*do,i,244,324
e,i,i+9000
*enddo
e,326,326+9000
*do,i,407,567
e,i,i+9000
*enddo
*do,i,649,729
e,i,i+9000
*enddo
*do,i,813,974
e,i,i+9000
*enddo
*do,i,1462,1544
e,i,i+9000
*enddo
*do,i,1785,1944
e,i,i+9000
*enddo
*do,i,2917,3240
e,i,i+9000
*enddo
allsel
!约束
nsel,s,loc,y,-0.5
d,all,all
allsel


nsel,s,loc,x,0
d,all,all
nsel,s,loc,x,40*0.625
d,all,all
allsel




/solu
antype,modal
modopt,lanb,10
mxpand,10
solve
望各位不吝指教，真的纠缠好久了

zhshle1010

咋没人回答列

Chelsea

现在用的是笔记本，没办法运行你的命令流
有什么错误信息或者你觉得有什么地方不对吗

zhshle1010

你好，谢谢回答，计算的结果都为0，没有结果，已查明是梁单元实常数没输完整导致的，但计算出来的结果应该还是存在问题，计算结果好像不正确，振型图也挺怪的