【压电声场求助】仿照网友帖子做的，出现错误，请教版主及各位老师(内附代码)

shenda1304

我仿照网友icecreamhan的帖子（http://forum.simwe.com/viewthread.php?tid=855600&;highlight=%C9%F9%D1%B9）中给出的例子做的，但是出现了下面这个错误，
Radius of the FLUID130 Acoustic Absorbing element 61731 should be positive; instead found a radius of 0. Please verify the element Real Constant values.

我设置的FLUID130的半径是正的啊，为什么会出现这种错误？

下面是我的代码，请大虾帮我看看，谢谢！
finish
/clear
/prep7
!PZT-5H材料特性参数
mp,DENS,1,7500                 
mp,perx,1,1704                  
mp,pery,1,1704
mp,perz,1,1433

tb,ANEL,1                              
tbdata,1,1.2720E11,8.0212E10,8.4670E10  
tbdata,7,1.2720E11,8.4670E10            
tbdata,12,1.1744E11                     
tbdata,16,2.3496E10                    
tbdata,19,2.2989E10                    
tbdata,21,2.2989E10                  

tb,PIEZ,1                             
tbdata,3,-6.6228                  
tbdata,6,-6.6228               
tbdata,9,23.2403
tbdata,14,17.0345               
tbdata,16,17.0345            

!定义Si层材料参数  
mp,dens,2,2330                        
MP,EX,2,120e9
MP,nuxy,2,0.42

!材料3为空气介质
MP,DENS,3,1.21
MP,SONC,3,344

!单元类型
et,1,solid5,3
et,2,solid45

ET,3,FLUID30

ET,4,FLUID130
r,4,0.002,0,0  !这里是正数的啊
            
ET,5,FLUID30,,1,0

!设置尺寸参数
*set,pzt_x_y,2e-3
*set,pzt_z,30e-6
*set,si_x_y,2e-3
*set,si_z,-10e-6

!!1/4个振膜
!!上半部分
block,0,pzt_x_y/2,0,pzt_x_y/2,0,pzt_z
block,-pzt_x_y/2,0,0,pzt_x_y/2,0,pzt_z

block,0,si_x_y/2,0,si_x_y/2,si_z,0
block,-si_x_y/2,0,0,si_x_y/2,si_z,0

!!下半部分
block,0,pzt_x_y/2,0,-pzt_x_y/2,0,pzt_z
block,-pzt_x_y/2,0,0,-pzt_x_y/2,0,pzt_z

block,0,si_x_y/2,0,-si_x_y/2,si_z,0
block,-si_x_y/2,0,0,-si_x_y/2,si_z,0

numcmp,all

vglue,all
SPHERE,,0.002,0,90
VOVLAP,all
!vglue,all
numcmp,all
alls

!对编号为1_6_7_8的体划分网格PZT
lesize,5,,,10,1
lesize,6,,,10,1
lesize,7,,,10,1
lesize,8,,,10,1

lesize,45,,,10,1
lesize,19,,,10,1
lesize,44,,,10,1

lesize,51,,,10,1
lesize,53,,,10,1
lesize,40,,,10,1

lesize,29,,,10,1
lesize,41,,,10,1
!!
lesize,38,,,10,1
lesize,52,,,10,1
lesize,50,,,10,1
lesize,43,,,10,1

lesize,1,,,10,1
lesize,13,,,10,1
lesize,42,,,10,1

lesize,2,,,10,1
lesize,3,,,10,1
lesize,4,,,10,1

lesize,28,,,10,1
lesize,39,,,10,1
!!
lesize,11,,,2,1
lesize,12,,,2,1
lesize,22,,,2,1
lesize,31,,,2,1

lesize,30,,,2,1
lesize,32,,,2,1
lesize,20,,,2,1
lesize,9,,,2,1

VSEL,s,,,1
VSEL,a,,,6
VSEL,a,,,7
VSEL,a,,,8
vatt,1,,1
mshkey,1
VMESH,all
alls

!对编号为2_3_4_9的体划分网格si
lesize,49,,,10,1
lesize,27,,,10,1
lesize,23,,,10,1
lesize,48,,,10,1

lesize,24,,,10,1
lesize,25,,,10,1
lesize,26,,,10,1

lesize,55,,,10,1
lesize,53,,,10,1
lesize,54,,,10,1

lesize,60,,,10,1
lesize,58,,,10,1
!!
lesize,36,,,1,1
lesize,34,,,1,1
lesize,46,,,1,1
lesize,47,,,1,1
lesize,59,,,1,1
lesize,56,,,1,1
lesize,57,,,1,1
lesize,35,,,1,1

VSEL,s,,,2
VSEL,a,,,3
VSEL,a,,,4
VSEL,a,,,9
vatt,2,,2
mshkey,1
VMESH,all
alls

!对编号为5的体划分网格(!!!!!此步骤划分网格需要手动进行!!!!)
VSEL,s,,,5
VATT,3,,3  ! FLUID30和空气 指定的材料属性为空气介质!!!(与结构耦合)!!!
SMRT,1
VMESH,5
alls

!修改不与平板（结构）接触单元的材料类型类型为空气介质!!!（不与结构耦合）!!!
esel,s,type,,1
esel,a,type,,2
nsle,s !p404
esln,s,0  
nsle,s
esel,inve !p402
nsle,s
emodif,all,type,5 !p214
alls

!在!!边界!!生成吸收单元。（使用单元fluid130）
Asel,s,,,7
NSLA,S,1
type,4   !FLUID130
real,4   
mat,3    !空气
esurf    !p208            
alls

!建立流－固接触面                        
asel,s,,,2
asel,a,,,4
asel,a,,,6

asel,a,,,11
asel,a,,,15
asel,a,,,16

nsla,s,1                                

esel,s,type,,3                   !选择附属于这些节点上的流体单元。
sf,all,fsi,1 !p275               !指定所选择的节点作为流体-结构接触面上的节点。
alls                             !注意—在对分析进行求解前确保重选所有的节点。

!施加约束                    
asel,s,,,4
asel,a,,,15
asel,a,,,24
asel,a,,,27
asel,a,,,6
asel,a,,,16
asel,a,,,20
asel,a,,,36

asel,a,,,14
asel,a,,,34
asel,a,,,30
asel,a,,,39
asel,a,,,29
asel,a,,,37
asel,a,,,9
asel,a,,,25

nsla,s,1
d,all,ux,0
d,all,uy,0
d,all,uz,0
alls

!施加法向对称约束!!!!!!!没有信心?????????
asel,s,,,12
nsla,s,1
dsym,symm,y,
alls

asel,s,,,13
nsla,s,1
dsym,symm,x,
alls

/solu
antype,harmic
hropt,full

!施加电压
asel,s,,,2
asel,a,,,22
asel,a,,,18
asel,a,,,32
nsla,s,1
!!为什么不能耦合节点电压?????????
!cp,1,volt,all
!*get,n_top,node,0,num,min
d,all,volt,10
alls

asel,s,,,1
asel,a,,,17
asel,a,,,21
asel,a,,,31
nsla,s,1
!cp,1,volt,all
!*get,n_bot,node,0,num,min
d,all,volt,0
alls

!d,n_bot,volt,0
!d,n_top,volt,10
!alls

!alls
nsubst,5
kbc,1
HARF,10000,100000            
SOLVE

!对与结构耦合的空气单元和未与结构耦合的空气单元进行后处理，得到声压分布
/post1
SET,FIRST
esel,s,type,,3
esel,a,type,,5
PLESOL,NMISC,4,0,1.0  !p321

shenda1304

另外就是
1.我在1/4球面的两个对称面上施加法向对称约束，不知道这样做对不对？
2.我在施加电压载荷时发现，不能先耦合再加电压，而只能选择全部节点然后在全部节点上施加，不知这是为什么？
3.我最后是想把指向性圆图整出来，现在才只做到这一步

shenda1304

很奇怪,我用1/2的薄板和1/2的空气建模可以计算,没弹出之前的那个错误,但是为什么用整个薄板和1/4的空气就不行呢?其他的参数我都没改啊,搞不懂。

tzhdacy

我现在机子上没有ansys加上你上传出错信息图片太小。我只是猜测，你1/4模型是不是流体是1/4，固体是半个或者1个？如果是，注意你的流固耦合单元设定

shenda1304

谢谢版主老师，
楼顶的模型：全部固体薄膜，1/4流体，计算出现如下错误
*** ERROR ***     SUPPRESSED MESSAGE    CP =      18.359   TIME= 13:33:18 Radius of the FLUID130 Acoustic Absorbing element 11198 should be        positive; instead found a radius of 0.                                    Please verify the element Real Constant values.                        
*** ERROR ***     SUPPRESSED MESSAGE    CP =      18.359   TIME= 13:33:18 Real constant set 4 undefined but referenced by element 11198.
3楼的模型：1/2固体薄膜，1/2流体，没有错误可以计算。

问题是，
3楼的模型和1楼的模型相比只是尺寸上做了修改，其他的参数都没变，3楼的可以计算，为什么1楼的会出现上述错误呢？我的130单元的实常数半径是0.002m啊，怎么还出现了上述错误中的“ instead found a radius of 0”。

tzhdacy

看错误提示：应该还是流固耦合设置时候出现些错误。
你是不是1/2变成1/4模型的时候 流固耦合界面没有修改，仍然是以前的1/2模型选取范围？

shenda1304

谢谢版主老师，
流固耦合界面应该是已经修改了的，请您看下面这幅图，这是1/4模型时流固耦合界面设置后的显示，

另外请您看一下这幅图

这是论坛上一个网友做的，您看他是用什么方法划分网格的？谢谢

tzhdacy

从这几个角度看的话，好像没有什么问题。你计算一下，看看结果。如果有时间可以把自己做的经验总结下，让后人学习。

shenda1304

谢谢版主老师，流固耦合界面没问题的话那是什么问题呢？
全部固体薄膜和1/4流体，计算时出现的错误
*** ERROR ***     SUPPRESSED MESSAGE    CP =      18.359   TIME= 13:33:18 Radius of the FLUID130 Acoustic Absorbing element 11198 should be        positive; instead found a radius of 0.                                    Please verify the element Real Constant values.                        
*** ERROR ***     SUPPRESSED MESSAGE    CP =      18.359   TIME= 13:33:18 Real constant set 4 undefined but referenced by element 11198.

tzhdacy

我感觉你的流固耦合的单元还是没有选择好，下面一个网友划分的网格属于映射划分，整体的感觉比较清晰。你用esel,s,,,,11198,找到这个单元,仔细看看他的位置

tzhdacy

最外层的边界条件怎么没有看见？就是sf,all,impd这个命令。最近工作比较繁忙，不能在电脑上运行你命令流。你自己先找找错误。

tzhdacy

你的吸收边界设定有问题，还有是不是没有的设定mat,,,mu 吸收系数。如果在最外面mu=1和fluid130互相配合使用模拟吸收边界。其他空气可以设定mu=0 表示无吸收

shenda1304

好的，谢谢版主老师，我再试试

shenda1304

版主老师我重新划分了网格按照您的建议做几遍还是不行，还是出现这个错误：
Radius of the FLUID130 Acoustic Absorbing element 61731 should be positive; instead found a radius of 0. Please verify the element Real Constant values.
我查看了报错的所有单元，显示出来发现只是是球面的一部分单元大约占整个球面的2/3左右，我按照您的建议把fluid30的mu也加上了，fluid130貌似没有mu这个属性哦。
请您再帮我看看到底是什么问题，另外如果不用fluid130改用mu=1的fluid30您看行吗，也就是球面的节点属性换成mu=1的fluid30。谢谢您啦！

shenda1304

另外我把出现错误的这些单元都选了出来，用emodif,all更改他们的实常数半径还是不行，这到底是什么原因啊？

shenda1304

刚刚用mu=1的fluid30（用来代替fluid130）作为球面全吸收边界做了一个声压级的图，

感觉比较混乱，也不知道这样做对不对？

hy1966

可以试试：
1、再增加个实常数参考声压，对应于fluid30，编号从1开始，现在实常数编号就是2，不要用4，应该可以解决半径为0的问题；
2、划分密一些、较规则的六面体网格，把大空间化成多个小块，然后设置单元的尺寸控制，提高计算精度，平滑显示；
3、如果用1/4模型，固体流体均可用1/4，这个模型用1/16应该也可以吧，减少计算量。

tzhdacy

对， 换成fluid30是可以的。当年我比较过结果和fluid130，差别不大。上面hy1996说的节点坐标转换诗歌需要注意问题。最近在忙一个项目，比较忙，有机会我抽空帮你运行下命令流

tzhdacy

声压得beam pattern,我只看见了main lobe,他的周围是不是还应该有minor lobe?

活耗子

是不是网格划分不好引起的错误？