一种计算溢流环换能器阵声辐射特性的方法
【专利摘要】本发明公开了一种计算溢流环换能器阵声辐射特性的方法，包括如下步骤：(1)建立溢流环换能器阵的有限元模型；(2)划分模型有限元网格；(3)进行有限元计算求解；(4)建立溢流环换能器阵的边界元模型；(5)结合有限元计算结果利用边界元方法计算声辐射特性。该计算溢流环换能器阵声辐射特性的方法，利用有限元方法结合边界元方法对溢流环换能器阵的声辐射特性进行计算，在溢流环换能器阵的设计阶段就对其性能进行预测，从而指导我们进行换能器阵的优化设计，得到性能优良的发射换能器阵，不仅计算精确，而且节省了大量的人力、物力和财力。
【专利说明】一种计算溢流环换能器阵声辐射特性的方法
【技术领域】
[0001]本发明涉及换能器基阵【技术领域】，尤其是涉及一种计算溢流环换能器阵声辐射特性的方法。
【背景技术】
[0002]溢流式圆环换能器，简称溢流环换能器，是一种水平无方向性的换能器，利用其液腔振动和径向振动的耦合，使得换能器具有低的谐振频率、宽的频带和高的效率，并且体积小、重量轻、功率大，还有优异的深水性能和良好的指向性。其结构简单紧凑，工作稳定可靠，在水声领域得到了广泛的应用。为了提高水下声源的发射声源级、发射指向性等性能，可以把多个溢流环换能器组阵发射。从而，我们需要深入研究溢流环换能器阵的发射性能，对其声福射特性进行计算。McMahon给出了溢流环换能器液腔谐振频率和径向谐振频率的近似计算公式，但是McMahon没有给出溢流环换能器其它参数的计算方法，例如换能器的发射电压响应、阻抗特性以及发射指向性等。为了得到溢流环换能器及基阵的发射电压响应、阻抗特性以及发射指向性等特性，通常采用实验测量的方法。利用阻抗分析仪测量其阻抗特性，利用测量换能器及基阵在水下工作时发射声压的大小来得到其发射电压响应及发射指向性等特性。
[0003]采用McMahon给出的理论公式计算的溢流环换能器液腔谐振频率和径向谐振频率，只是一个近似计算公式，不是很准确。并且理论公式只能计算谐振频率，对发射电压响应、阻抗特性以及发射指向性等其它特性无法计算，更不能计算溢流环换能器阵的发射特性。而采用实验测量的方法得到溢流环换能器及基阵的发射电压响应、阻抗特性以及发射指向性等特性。这种实验测量的方法需要耗费大量的人力、物力、财力，同时测量受到水声环境、人为操作误差和实验系统误差的影响。并且，实验测量方法不能在基阵设计阶段对其性能进行预测，只能制作出来后再测量其性能好坏，性能不好还需要重新设计、制作，浪费大量的人力和财力。
【发明内容】
[0004]本发明所要解决的技术问题是:提供一种计算溢流环换能器阵声辐射特性的方法，利用有限元方法结合边界元方法对溢流环换能器阵的声辐射特性进行计算，在溢流环换能器阵的设计阶段就对其性能进行预测，从而指导我们进行换能器阵的优化设计，得到性能优良的发射换能器阵，不仅计算精确，而且节省了大量的人力物力和财力。
[0005]为解决上述技术问题，本发明的技术方案是:一种计算溢流环换能器阵声辐射特性的方法，包括如下步骤:
[0006](I)建立溢流环换能器阵的有限元模型:基于溢流环换能器及基阵的结构为上下对称和轴对称，只需建立换能器阵上半部分的轴对称有限元模型，首先给出溢流环换能器的内半径尺寸、外半径尺寸、高度尺寸以及换能器阵的阵元间距，使用ANSYS软件建立溢流环换能器阵的有限元模型；[0007](2)划分模型有限元网格:在ANSYS软件中，指定所要划分的溢流环换能器阵及其所在流体域有限元网格的单元类型和材料参数，设定网格的大小，然后用自由划分网格的方法对所建立的换能器阵的有限元模型进行网格划分，从而得到换能器阵的有限元模型网格；
[0008](3)进行有限元计算求解:在ANSYS软件中利用模态分析、谐波响应分析方法计算出溢流环换能器及基阵的振动模态、谐振频率、导纳特性、发射电压响应、近场声特性以及振动位移分布特性；
[0009](4)建立溢流环换能器阵的边界元模型:在ANSYS软件中，提取溢流环换能器阵表面的单元网格，导入到SYSNOISE软件中，从而得到换能器阵的边界元模型网格；
[0010](5)结合有限元计算结果利用边界元方法计算声辐射特性:把在ANSYS软件中用有限元方法对溢流环换能器阵进行谐波响应分析得到的换能器阵的表面振动位移分布导入到SYSNOISE软件中，利用边界元方法计算溢流环换能器阵的辐射声场远场特性和辐射指向性特性。
[0011]优选的，步骤⑴中，在ANYSYS软件中利用画图操作画出溢流环换能器阵的结构模型图。
[0012]优选的，步骤(2)中，需要设定溢流环换能器阵压电材料的刚度矩阵、压电矩阵及介电矩阵，在ANSYS软件中，用来解决压电问题的压电关系的方程为:
[0013]{T} = [cE] {S}-[e] {E},
[0014]{D} = [e]T{S} + [es] {E},
[0015]其中，{T}为应力向量；
[0016]{S}为应变向量；
[0017]{D}为电位移向量；
[0018]{E}为电场向量；
[0019][cE]为压电材料在常数电场下的刚度矩阵；
[0020][e]为压电材料的压电矩阵；
[0021][ ε s]为压电材料在常数应变下的介电矩阵；
[0022]在ANSYS软件中利用画图操作对溢流环换能器阵的结构模型进行网格划分。
[0023]优选的,步骤(3)中,把溢流环换能器及基阵当作一个振动体结构,结构周围流体域为水，对于流固耦合的声学问题，结构上的有限元方程为:
[0024]Mii+ Cii+ Ku-RP = F
[0025]其中M、C和K分别为振动结构的质量矩阵、阻尼矩阵和刚度矩阵，Ι!、?和U分别为结构有限元节点的位移向量、速度向量和加速度向量，R为流固耦合界面上的耦合矩阵，P为流固耦合界面上流体域节点上的声压向量，F为结构节点上的力向量；
[0026]对于声传播的流体域，由声学波动方程考虑流固耦合界面上的相互作用及流体域边界上的衰减，可得到流体域的有限元方程为:
[0027]MfP H- CfP + KfP + /70rT? = O
[0028]其中Mf、Cf和Kf分别为流体的质量矩阵、阻尼矩阵和刚度矩阵，Pj和？分别为流体域有限元节点的声压向量、声压向量对时间的一阶导数和二阶导数，P ^为流体的密度，R为流固耦合界面上的耦合矩阵，?为流固耦合界面上结构节点的加速度向量；[0029]给定振动结构的参数(包括几何参数和材料参数等)及流体域的参数(包括流体的大小、声速、密度等)，有限元网格生成后，M、C、K、Mf、Cf、Kf和R等矩阵就可以完全唯一确定了，再给出结构节点上的力向量F(对于压电有限元给出所施加的电压)，联合两式就可以求解出结构节点的位移向量和流体域节点的声压向量。当结构节点上的力向量F为零(压电有限元所施加的电压为零)时联合两式求解可对结构体进行模态分析，即分析结构不施加外力时的谐振频率和振动形状；
[0030]在溢流环换能器的内外壁上施加电压，对换能器施加IV的电压。换能器的导纳为电流与电压之比，即为:
[0031]
【权利要求】
1.一种计算溢流环换能器阵声辐射特性的方法，其特征在于，包括如下步骤: (1)建立溢流环换能器阵的有限元模型:基于溢流环换能器及基阵的结构为上下对称和轴对称，只需建立换能器阵上半部分的轴对称有限元模型，首先给出溢流环换能器的内半径尺寸、外半径尺寸、高度尺寸以及换能器阵的阵元间距，使用ANSYS软件建立溢流环换能器阵的有限元模型； (2)划分模型有限元网格:在ANSYS软件中，指定所要划分的溢流环换能器阵及其所在流体域有限元网格的单元类型和材料参数，设定网格的大小，然后用自由划分网格的方法对所建立的换能器阵的有限元模型进行网格划分，从而得到换能器阵的有限元模型网格； (3)进行有限元计算求解:在ANSYS软件中利用模态分析、谐波响应分析方法计算出溢流环换能器及基阵的振动模态、谐振频率、导纳特性、发射电压响应、近场声特性以及振动位移分布特性； (4)建立溢流环换能器阵的边界元模型:在ANSYS软件中，提取溢流环换能器阵表面的单元网格，导入到SYSNOISE软件中，从而得到换能器阵的边界元模型网格； (5)结合有限元计算结果利用边界元方法计算声辐射特性:把在ANSYS软件中用有限元方法对溢流环换能器阵进行谐波响应分析得到的换能器阵的表面振动位移分布导入到SYSNOISE软件中，利用边界元方法计算溢流环换能器阵的辐射声场远场特性和辐射指向性特性。
2.如权利要求所述的计算溢流环换能器阵声辐射特性的方法，其特征在于，步骤(1)中，在ANYSYS软件中利用画图操作画出溢流环换能器阵的结构模型图。
3.如权利要求所述的计算溢流环换能器阵声辐射特性的方法，其特征在于，步骤(2)中，需要设定溢流环换能器阵压电材料的刚度矩阵、压电矩阵及介电矩阵，在ANSYS软件中，用来解决压电问题的压电关系的方程为:
{T} = [cE] {s}-[e] {E}，
{0} = [e]T{S} + [es] {E}， 其中，{T}为应力向量； {S}为应变向量； {D}为电位移向量； {E}为电场向量； [cE]为压电材料在常数电场下的刚度矩阵； [e]为压电材料的压电矩阵； [ε S]为压电材料在常数应变下的介电矩阵； 在ANSYS软件中利用画图操作对溢流环换能器阵的结构模型进行网格划分。
4.如权利要求所述的计算溢流环换能器阵声辐射特性的方法，其特征在于，步骤(3)中，把溢流环换能器及基阵当作一个振动体结构，结构周围流体域为水，对于流固耦合的声学问题，结构上的有限元方程为:
Mu + Cu + Ku-RP = F 其中M、C和K分别为振动结构的质量矩阵、阻尼矩阵和刚度矩阵，U、?和U分别为结构有限元节点的位移向量、速度向量和加速度向量，R为流固耦合界面上的耦合矩阵，P为流固耦合界面上流体域节点上的声压向量，F为结构节点上的力向量；对于声传播的流体域，由声学波动方程考虑流固耦合界面上的相互作用及流体域边界上的衰减，可得到流体域的有限元方程为:
MfP + CfP + KfP + ^0RtU = O 其中Mf、Cf和Kf分别为流体的质量矩阵、阻尼矩阵和刚度矩阵，Pj和护分别为流体域有限元节点的声压向量、声压向量对时间的一阶导数和二阶导数，P ^为流体的密度，R为流固耦合界面上的耦合矩阵，?为流固耦合界面上结构节点的加速度向量； 给定振动结构的参数(包括几何参数和材料参数等)及流体域的参数(包括流体的大小、声速、密度等)，有限元网格生成后^、(:、1(為、(^、1^和1?等矩阵就可以完全唯一确定了，再给出结构节点上的力向量F(对于压电有限元给出所施加的电压)，联合两式就可以求解出结构节点的位移向量和流体域节点的声压向量。当结构节点上的力向量F为零(压电有限元所施加的电压为零)时联合两式求解可对结构体进行模态分析，即分析结构不施加外力时的谐振频率和振动形状； 在溢流环换能器的内外壁上施加电压，对换能器施加IV的电压。换能器的导纳为电流与电压之比，即为:
Y = L = Q + jB 其中，Y表不换能器的导纳，V表不在换能器换能材料上所施加的电压，I表不所产生的电流，G为电导，B为电纳； 电流通过换能材料电极上所集聚的电荷计算得到:
I = j2 Ji fQ 其中f为谐波响应分析的频率，Q为换能材料电极上所集聚的总电荷； 换能器的发射电压响应为:在换能器的声轴方向上，距离其有效声中心Cltl米距离上产生的球面波自由场声压Pf与加到换能器输入端的激励电压V的比值:
Sv = Pf.d0/V 当用分贝(dB)表示时，则称为发射电压响应级，即:
S1vZ = 20 Ig-^~ = 20 lg(-^-) + 120
Wv/re/^ 其中基准值(Sv)rtf取I μ Pa.m/V, d0取lm。 在ANSYS软件中提取出不同频率下溢流环换能器及基阵在施加一定电压下，与溢流环径向方向一致的声轴方向上的远场辐射声压值，即可计算出换能器及基阵的发射电压响应。
5.如权利要求所述的计算溢流环换能器阵声辐射特性的方法，其特征在于，步骤(4)中，在ANYSYS软件中利用画图操作得到溢流环换能器阵表面的的结构模型图，并进行网格划分，把结果导入到SYSNOISE软件中得到边界元模型。
6.如权利要求所述的计算溢流环换能器阵声辐射特性的方法，其特征在于，步骤(5)中，把Helmholtz边界积分方程的积分形式离散化后转化为矩阵形式。对于外部积分方程和表面积分方程，对积分曲面Stl进行离散，划分成有N个节点的边界元，这些节点上的声压构成列向量记为{P,}，法向振速也构成列向量记为IvJ，在声场中选择M个点,这些点上的声压构成列向量记为{pQ}。满足方程:{pQ} = [A] {pJ + [B] {vj
[H] {pj = [G] {vn} 其中[A]，[B], [H]，[G]为系数矩阵，都可以通过积分求得； 还可以得到:
{pj = [Hr1M Iv1J
{pQ} = ([A] [Hr1M+ [B]) KJ
令[Cl = [A] [Hr1EG]+ [B]代入上式可得
{pQ} = [C] IvJ 溢流环换能器及基阵表面上以及辐射声场中任意一点的声压都可以表示成表面法向振速向量的线性组合，其组合系数由矩阵[H]、[G]、[C]来决定，而这几个矩阵与振速无关，只与振动系统本身有关，包括振动面的几何形状、边界面阻抗特性、声传播媒质的物理特性(声速、密度)、频率以及场点位置。
【文档编号】G06F17/50GK104008237SQ201410212447
【公开日】2014年8月27日 申请日期:2014年5月16日 优先权日:2013年8月31日
【发明者】何正耀, 马远良, 杨坤德, 雷波, 韩一娜 申请人:西北工业大学