声场在固体中传播问题

ansys平安

在模拟声波在圆柱体中传播时遇到这样的错误


Radius of the FLUID130 Acoustic Absorbing element 40325 should be positive;instead found a radius of 0
Please verify the element real constant values.

但是在设置R时，设置为r,3,30,0,0时，又会出现  12162 of the FLUID130 Acoustic Absorbing element 40325 should lie on a SPHERE of radius 30.
instead lies at a radius of 30.016662.

请大家给分析分析。
命令流如下：

FINI
/FILENAME,JK31
/TITLE, GAS OF COAL SEEPAGE
/PLOPTS,DATE,0
/PREP7
ET,1,CPT216
KEYOPT,1,2,1
ET,2,FLUID30
ET,3,FLUID130
R,3,30,0,0
K=2.2e-6
MP,EX,1,185E3
MP,PRXY,1,0.32
MP,DENS,1,1.32e-6
TB,DP,1
TBDATA,1,0.3E3
TBDATA,2,32
TB,PM,1,,,PERM
TBDATA,1,k
TB,PM,1,,,BIOT
TBDATA,1,0.8
TB,PM,1,,,SP
TBDATA,1,118E3,,0.8
MP,DENS,2,1030
MP,SONC,2,1460
PCIRC,25,0,0,360
ADELE,ALL,,,0
K,5,12.5,12.5,0
K,6,12.5,-12.5,0
K,7,-12.5,12.5,0
K,8,-12.5,-12.5,0
K,9,-30,-30,0
K,10,30,30,0
K,11,30,-30,0
K,12,-30,30,0
K,13,0,0,0
K,14,17.5,0,0
K,15,-17.5,0,0
K,16,0,17.5,0
K,17,0,-17.5,0
L,5,10
L,5,16
L,16,7
L,7,15
L,8,15
L,7,12
L,8,17
L,17,6
L,6,14
L,14,5
L,8,9
L,6,11
L,3,15
L,15,13
L,17,4
L,14,1
L,16,2
L,16,13
L,14,13
L,17,13
LSBL,3,15
L,18,8
LSBL,4,16
L,6,9
LSBL,1,5
LSBL,2,10
L,5,11
L,7,10
AL,1,10,7,21
AL,10,8,17,5
AL,17,25,9,3
AL,3,26,19,11
AL,19,27,4,12
AL,4,15,20,13
AL,20,14,16,2
AL,2,6,21,28
AL,7,8,18,22
AL,18,9,11,24
AL,12,13,23,24
AL,23,14,6,22
AADD,9,12
AADD,10,11
AADD,9,13
VEXT,all,,,0,0,100
VGLUE,ALL
cylind,25,30,0,100,0,180
cylind,25,30,0,100,180,360
VGLUE,10,11
allsel,all
vsel,s,volu,,10
vsel,a,volu,,12
vsel,inve
vatt,1,1,1,0
aslv,s,1
lsla,s,1
lesize,all,,,3
vsweep,all
allsel,all
vsel,s,volu,,10
vsel,a,volu,,12
vatt,2,1,2,0
lesize,71,,,50
lesize,69,,,50
lesize,92,,,50
lesize,91,,,50
lesize,72,,,4
lesize,70,,,4
lesize,68,,,50
lesize,66,,,50
lesize,90,,,50
lesize,89,,,50
lesize,67,,,4
lesize,65,,,4
esize,75,,,100
lesize,73,,,100
lesize,74,,,100
lesize,76,,,100
vsweep,all
allsel,all

!define pres
ASEL,S,AREA,,45
ASEL,A,AREA,,57
NSLA,S,1
TYPE,3
REAL,3
MAT,2
ESURF
ALLSEL,ALL
ASEL,S,AREA,,46
ASEL,A,AREA,,58
NSLA,S,1
ESEL,S,TYPE,,2
SF,ALL,FSI,1
ALLSEL,ALL
ASEL,S,AREA,,42
ASEL,A,AREA,,15
ASEL,A,AREA,,19
ASEL,A,AREA,,23
ASEL,A,AREA,,27
ASEL,A,AREA,,31
ASEL,A,AREA,,35
ASEL,A,AREA,,39
ASEL,A,AREA,,9
NSLA,S,1
D,ALL,PRES,1.5E3
ALLSEL,ALL
ASEL,S,AREA,,10
ASEL,A,AREA,,1
ASEL,A,AREA,,2
ASEL,A,AREA,,3
ASEL,A,AREA,,4
ASEL,A,AREA,,5
ASEL,A,AREA,,6
ASEL,A,AREA,,7
ASEL,A,AREA,,8
NSLA,S,1
D,ALL,PRES,0.1E3
D,ALL,UZ,0
SF,ALL,PRES,%F%
ALLSEL,ALL
ASEL,S,AREA,,40
ASEL,A,AREA,,11
ASEL,A,AREA,,16
ASEL,A,AREA,,22
ASEL,A,AREA,,26
ASEL,A,AREA,,28
ASEL,A,AREA,,32
ASEL,A,AREA,,36
NSLA,S,1
SF,ALL,PRES,3E3
ALLSEL,ALL

/SOLU
ANTYPE,HARMIC
NROPT,UNSYMM
OUTRES,ALL,ALL
TIME,100
NSUBS,10,10,10
SOLVE

Alfred_Leo

我初步判断：这个分析的问题在于你将FLUID130用在一个柱面上了。在ANSYS帮助文档Mechanical APDL->Fluids Analysis Guide-> II. Acoustics 这部分中有这么一段：
When using FLUID129 and FLUID130, the boundary of the underlying finite element mesh must be circular (2-D and axisymmetric) or spherical (3-D), and the radius of the circular or spherical boundary of the finite domain must be specified as real constant RAD. (See Figure 1.2.) The center of the circle or sphere is also specified using real constants: R,3,RAD,X0,Y0!REAL set 3 for FLUID129
R,3,RAD,X0,Y0,Z0!REAL set 3 for FLUID130

tzhdacy

FLUID129 and FLUID130单元目前是用在2d的圆周和3d球面边界。如果是圆柱形边界免设定无反射，就别用这两个单元。