| 查看: 666 | 回复: 1 | ||
[求助]
模拟脉动球体声辐射时用到的surf153单元中的实常数ADMSUA(单位面积附加质量)怎么填
|
|
我做的一个例题，是使用fluid29单元和surf153单元模拟水中脉动球体声辐射。其中surf153单元的实常数ADMSUA，即单位面积附加质量例题中直接填了749989.484/(2*acos(-1)*160)，但是我不知道这个数字是怎么来的，我猜测是辐射阻抗里的同振质量Xr/w，但是算出来不是例题中的数据。请教一下这是怎么回事。下面贴出前处理的代码。加粗的那行就是实常数设置。 /clear,start /title,radiating sphere,f=160Hz,ka=6.702064328E-02,SURF153(fluid29) /filname,radiating sphere by using fluid29 /prep7 et,1,fluid29 keyopt,1,2,0 keyopt,1,3,1 et,2,fluid29 keyopt,2,2,1 keyopt,2,3,1 et,3,surf153 keyopt,3,3,1 keyopt,3,4,1 r,1,1e-6 r,2,1e-6 r,3, rmodif,3,6,749989.484/(2*acos(-1)*160) mp,sonc,1,1500 mp,dens,1,1000 mp,sonc,2,1500 mp,dens,2,1000 mp,visc,3,753971.665 pcirc,0.1,1.5,89,90 esize,(1500/3e4)/15 mshape,0 mshkey,1 aatt,2,2,2 amesh,all /pnum,type,1 /num,1 /auto eplot csys,1 nrotat,all type,1 real,1 mat,1 nsel,s,loc,x,0.1-1e-8,0.1+1e-8 nsel,a,loc,x,1.5-1e-8,1.5+1e-8 esln sf,all,fsi nsle emodif,all d,all,uy nsel,u,loc,x,0.1-1e-8,0.1+1e-8 nsel,u,loc,x,1.5-1e-8,1.5+1e-8 d,all,ux nsel,s,loc,x,0.1 d,all,ux,0,-0.005/(2*acos(-1)*160) type,3 real,3 mat,3 esel,all nsel,s,loc,x,1.5-1e-8,1.5+1e-8 esurf,all allsel,all fini |
回复此楼
2楼2014-08-28 18:13:31
