blackno11 发表于 2009-8-18 19:25

DBEM计算TL的过程

直接边界元法计算TL直接边界元与间接边界元应用的一个区别是,计算消声器传递损失时,DBEM模型里不能有交接面,也就是所有的面都在声场的边界上,在声场内部不存在边界。IBEM模型里可有。
New Name 'c-dbem' Model 1 File c-dbem.sdb Return
Option
BEM Direct Collocation Element Interior Uncoupled Frequency
Return
Environment Section SETUP ANSYSREVISION '5.7' Return
Environment Section SETUP MEMORY '200.000000' Return
Import Mesh Format Ansys File ibem.cdb Return
Check Mesh Set Subdomain Return
!检查网格
Reverse

!调整法线方向(如果需要)

Element All

Return
Set 1 Name "inlet"
!速度进口边界条件

Element Internal 75 To 87 712

Return
Set 2 Name "outlet"

Element Internal 88 To 100
!无反射出口边界条件

Return
Material Fluid
!流体材料(默认)

Name 'air'

Sound Real 3.4000e+002
Rho Real 1.2250e+000

Return
Boundary Velocity
Real 1 Imag 0
!速度边界条件(特别注意这里是正值)

Elements Set 1

Return
Boundary Impedance
Real 416.5 Imag 0
!无反射出口

Elements Set 2

Return
Save Return
Solve

Frequency 20 To 3000 LinStep 10

Return
Show Dialog File/Read return
Read File
v-post.cmd
Point
2.049164e-3 -7.2867977e-3 8.1355395e-34 Return
!场点定义在节点上(注意这点比较关键,特别是进口场点的位置,因为计算TL时要用到进口速度,当然也可以定义的其他的点,只是速度要读取,并要注意速度方向)
Point
-4.434042e-3 -7.4579693e-3 0.21 Return
PostProcess

Points All

Frequency 20 To 3000 LinStep 10

Near 2

Far 5

Quadrature 2 2 1

Save Results Step 1

Return
Combine
Read point 2 Pressure Return
Read point 2 Pressure Return
Conjugate
Multiply
Root
Real
Write Table 20 Name 'Pout' Return
Clear
Read point 1 Pressure Return
Constant Real 416.5 Imag 0
Constant Real 1 Imag 0
Multiply
Add
Read point 1 Pressure Return
Constant Real 416.5 Imag 0
Constant Real 1 Imag 0
Multiply
Add
Conjugate
Multiply
Root
Constant Real 2 Imag 0
Divide
Real
Write Table 21 Name 'Pin' Return
Clear
Read Table 20 Return
Read Table 21 Return
Divide
Square
Log10
Constant Real 10 Imag 0
Multiply
Write Table 22 Name 'TL-HZ1' Return
Return
Curve

Table 22

Format
1

NoBackground

Marker 0

Return
页: [1]
查看完整版本: DBEM计算TL的过程