【原创】请帮忙解决声分析中的模态分析问题

sunhaolan

就是把圆划成4X1/4圆，再划分。
你的主要问题是约束不够，自己加点约束。建模部分，我帮你做了，命令流如下：
/filnam,m  
/prep7  
/pnum,kp,1 $/pnum,line,1  
  
$/pnum,area,1  
et,1,30,,1  
et,2,29  
et,3,30  
et,4,63  
  
r,4,0.004  
mp,dens,4,1700  
mp,ex,4,2.1e10  
mp,nuxy,4,0.1  
r,1,1  
mp,sonc,1,1460  
mp,dens,1,1000  
  
!建模,利用一曲面绕y轴旋转而成，  
k,1,0,0  
k,2,0.07,0.4  
k,3,0.162,0  
k,4,0,0.4  
k,5,0.28,0  
k,6,0.2,0.4  
larc,2,3,1,0.5
larc,5,6,3,0.8
l,2,6  
l,3,5  
a,2,6,5,3  
  
type,2  
lesize,1,,,18,1
lesize,2,,,18,1
lesize,3,,,4,1
lesize,4,,,4,1
mshkey,1
mshape,0,2d
amesh,1
  
type,1  
real,1  
mat,1  
esize,0,6  
  
vrotat,1,,,,,,1,4,360 !生成体积  
  
!用壳单元划分壳体四个侧面：
  
type,4  
real,4  
mat,4  
esize,0.03,0  
amesh,5,20,5  
  
wprot,,-90   
CSYS,4   
CYL4,,,.162
wpoff,,,.4   
CYL4,,,.07
  
asel,s,area,,21,22   
lsla
ksll
L,26,24   
L,25,23   
L,21,19   
L,20,22   
lsel,s,line,,41,42   
asbl,22,all  
lsel,s,line,,43,44   
asbl,21,all  
  
!用壳单元底面进行划分,然后放样。  
type,4  
real,4  
mat,4  
lsel,s,line,,33,40
lesize,all,,,6
mshape,0,2d  
amesh,21,29
  
NUMMRG,ALL
csys,0
  
!用type1对壳底面进行放样，单元间隔
  
0.01  
  
!划一条线号为43, 101(0,0.4),102(0,0.8)  
k,101,,0.4  
k,102,,0.8
l,101,102  
  
!划一条线号为44, 103(0,0),104(0,-0.4)  
k,103,
k,104,,-.4  
l,103,104
  
type,1  
real,1  
mat,1  
esize,0.03,0  
vdrag,23,24,25,26,,,43  
vdrag,2,7,17,,,,43
vdrag,22,27,28,29,,,44  
vdrag,4,9,14,19,,,44  
  
!对面1清除平面单元29  
asel,s,,, 1  
aclear, 1  
  
nummrg,all !合并重叠点  
  
!选择壳体上的面 5,10,20,15,22-29  
!选择面上的结点，再选择结点上的单元  
!设置流固耦合面  
asel,s,,,5,20,5
asel,a,,,22,29
nsla !选取所选面上所有结点  
esln !选择结点上所有流体单元  
esel,r,type,,1
emodif,all,type,3  
  
sf,all,fsi,1  
  
!选择整个模型最外部的面及该些面的所有结点，结点的声压值设为0  
asel,s,,, 3  
asel,a,,, 8  
asel,a,,, 13  
asel,a,,, 18  
  
asel,a,,, 43  
asel,a,,, 48  
asel,a,,, 53  
asel,a,,, 56
  
asel,a,,, 74  
asel,a,,, 78  
asel,a,,, 82  
asel,a,,, 85
  
asel,a,loc,y,-.4   
asel,a,loc,y,.8  
nsla
  
d,all,pres,0.0  
allsel
fini

泡泡鱼

不愧为高手！万分感谢您！

泡泡鱼

不过还是很奇怪，1。“Coefficient ratio exceeds 1.0e8 - Check results.” 的警告依然存在，还有为什么有时候模态分析指定分析十阶，却只出现7阶呢？
2。我的机器内存不够，如果划分的过粗，是不是也有可能出现原先算不出来的警告呢

sunhaolan

那个警告不用理它。
模态分析的问题，是不是还是你的约束给的不合适？你怎么给的约束？

泡泡鱼

全部施加约束
asel,s,,,5,20,5
nsla
d,all,all
  
加约束和不加约束算下的结果肯定不会相同，计算这个圆台壳是不要求加约束的，但不加，又算不出来，加了约束会对结果影响多大？？

ynn1980

sf,all,fsi,1  
  
!选择整个模型最外部的22个面及该些面的所有结点，结点的声压值设为0  
asel,s,,, 37 $asel,a,,, 39 $asel,a,,, 34 $asel,a,,, 31 $asel,a,,, 27  
asel,a,,, 53 $asel,a,,, 55 $asel,a,,, 47 $asel,a,,, 50  
asel,a,,, 29 $asel,a,,, 38 $asel,a,,, 35 $asel,a,,, 32  
asel,a,,, 8 $asel,a,,, 13 $asel,a,,, 18 $asel,a,,, 3  
asel,a,,, 56 $asel,a,,, 48 $asel,a,,, 51 $asel,a,,, 54 $asel,a,,, 44  
nsla,,  
  
d,all,pres,0.0  
allsel,,   
fini  
  
我想问一下，为什么要加压力为零，把结构外面的流体球加上fsi不就可以了么？为什么要把所有的压力设为0?
谢谢！希望能尽快得到sunhaolan的回答。你是高手啊！