圆柱壳全频段声振特性研究

1、圆柱壳全频段声振特性研究指导老师：张亚辉老师杨文健2015.12.26该论文的主要内容有： 基于解析法软件的圆柱壳中低频段声振响应分析 基于VA One的圆柱壳全频段声振响应分析 对未来的工作展望1.基于解析法软件的圆柱壳中低频段声振响应分析u 解析法软件理论基础 解析法求解单、双层圆柱壳振动与声辐射响应的主要原理 解析法求解单、双层圆柱壳固有频率的主要原理u 基于Matlab 平台的软件制作方法简介u 解析法软件正确性的验证解析法求解单、双层圆柱壳振动与声辐射响应的主要原理解析法求解单、双层圆柱壳固有频率的主要原理将圆柱壳运动方程中的力激励项去掉，求剩下项的行列式等于0时的频率值。由于该方程

2、组是超越方程所以只能搜索求得方程的解。具体方法具体方法:逐个搜索每个频率，如果搜索到的频率的行列式值与下一个频率的行列式值反号，则表示搜索到的频率的行列式值近似为零，取该频率为固有频率。基于Matlab平台的软件制作方法简介在Matlab平台上开发软件界面通常分为三个部分：软件静态界面制作软件静态界面制作、内核程序编写内核程序编写和程序封装程序封装。p 软件静态界面制作 静态界面的制作就是在空白的GUI模板上加入各种空间的过程，一般来说添加控件有2种方法： 1. 基于GUI的方法：简单说就是手动在图形窗口的所需的位置上添加所需大小的控件，并且双击控件出现属性菜单，手动设置控件的属性值。 2.

3、基于命令行的方式：通过直接在M文件中编写添加控件的程序代码 来完成静态界面的制作。命令格式为： H为uicontrol(hfig,属性名，属性值，)H为创建对象的句柄值，简单理解就是控件的代号Hfig为某个图形窗口的句柄值p 软件内核程序编写静态界面制作完成后，就要编写动态程序，通俗的说，就是编写控件对应的响应函数，也就是在原来无界面的程序的基础上进行修改，添加一些控制语句来配合界面控件的操作，从而实现界面的功能。每一个控件都对应有一个响应函数，当用户在界面上对某个控件进行操作时，matlab就会调用其响应函数。p 将程序封装成exe文件软件的外观和内核程序代码制作好以后，就要将这些内在外在都

4、封装到一个exe文件中，方便提交用户使用。封装的步骤如下:1.首先需要配置自己的Matlab Compiler, Matlab Compiler的作用是将程序编译成为机器可以直接执行的程序。 配置Compiler的方法是在Matlab命令窗口输入:mbuild一setup按提示选择matlab自带编译器LCCo2.用编译器编译所有文件 在matlab环境中编译文件，命令是: mcc -m main funl fun2funn其中main为主函数，funl到funn是子函数。执行完毕后，在MATLAB的CurrentDirectory目录下，会生成封装后的exe文件，至此，软件就制作成功了。解析

5、法软件正确性的验证选取单层环肋圆柱壳为计算模型：单双层圆柱壳声振试验实验模型三种方法结果对比验证将数值法数值法、实验法实验法的结果与解析法软件解析法软件计算的结果进行对比：l 由均方振速曲线可以看出，解析法和数值法计算结果趋势基本相同，峰值位置、大小接近。在150Hz-500Hz，两种计算方法吻合的较好，峰值位置出现较小的偏移;在500Hz-1000Hz，解析法计算结果曲线整体向左偏移，数值略大。l 由辐射声功率曲线可以看出，解析法和数值法吻合的较好。在150Hz-500Hz 两条曲线基本重合;在500Hz-1000Hz，两条曲线峰值位置出现偏移。l 由加速度曲线可以看出，在150Hz-500

6、Hz三条曲线吻合的较好;在500Hz-1000Hz，试验结果曲线相对于解析法和数值法出现偏移。误差来源l 解析法误差主要来源： 采用模态叠加法进行求解，模态截断及忽略模态间的耦合都会使得计算结果出现误差。l 边界元方法误差主要来源： 网格划分、积分奇异性的处理。固有频率计算结果对比本程序计算结果与数值法在模态阶数较低时固有频率吻合情况较好，模态阶数升高后，数值法结果略有偏高。本章小结p 介绍了解析法的基本原理p 在MATLAB平台上开发软件的过程p 将该软件的计算结果与传统数值计算方法和试验数据进行对比，验证了该软件的正确性 2. 基于VA One 的圆柱壳全频段声振响应分析u VA One

7、理论基础u VA One 的核心功能和简要的操作流程u 单双层圆柱壳声振响应全频段算例计算u 结构参数对单双层圆柱壳全频段声振响应的影响u VA One 计算结果与解析法的对比u 本章小结 VA One 理论基础VA One软件低频段中频段高频段有限元和边界元方法Hybrid方法统计能量分析方法有限元方程为：可求得结构各节点的位移 统计能量分析方法dwwdNwn)()(估算各子系统的动力响应划分子系统模态相似原则：阵型要有着相同的动力学特性（相同的阻尼、振动模态、耦合损耗因子）一根梁有纵向振动子系统（即纵向振动相似模态群）、横向振动子系统（横向振动相似模态群）和扭转振动子系统（扭转振动相似模态

8、群）确定各子系统与各子系统间的统计能量分析参数计算各子系统振动能量模态密度dKKdNKn)()(模态数CK2波数振动偏微分方程振型函数频率方程模态数N(K)模态密度分离变量法边界条件波数空间单位频率内的模态数一维杆纵向自由振动微分方程：分离变量法：杆件为两端固支二维平板横向振动方程为：同样的，我们可以通过分离变量法得到振型函数满足的方程：半径为K的1/4圆面积内模态数目为：其中：输入功率VFP统计能量分析主要应用于受随机激励的系统，但这里应强调的一点是，统计能量分析最主要的特征之一是只建立被激励系统的统计模型，不专门对随机激励建立统计模型。考虑线性弹簧共振子受到一个稳态随机激励的外力：声振耦合

9、问题中：声快振动声慢振动由于有限板的二维性及边界条件产生驻波的影响内损耗因子模态内损耗因子测量方法 模态内损耗因子测量方法（稳态法） 带均（频带平均）内损耗因子测量方法 现场测量方法耦合损耗因子两个共振子间的耦合损耗因子写成对称的形式由外部供给的总输入功率，都被二个共振子的阻尼元件耗散掉了共振子间的传递功率虚拟激励法解得利用时域中的均方值与频域中的功率潜密度关系式:通过虚拟激励法可得：A只与共振子参只与共振子参数及耦合参数有数及耦合参数有关关两个线性共振子间功率流正比于两个共振子解耦能量之差二振型群间的耦合损耗因子假定：Hybrid 方法整个系统的总体运动方程系统的平均互谱响应通过有限元、边界

10、元的方法对其进行求解可得extffS通过混响场的功率流平衡方程VA One 的核心功能 低频结构振动和噪声分析的综合性功能 高频结构振动及噪声分析功能 中频结构振动及噪声混合求解功能模型生成参数设定加载、建立连接、计算结果查看l 各个子系统的平均能量输入l 各子系统之间能量传递系数l 结构均方振速l 振动加速度l 辐射声功率l 辐射效率l 辐射声压l 模态能量操作流程单层壳算例空气中计算结果 不论是均方速度还是升功率曲线，在低频段都呈现多峰值形态，而中高频段没有明显峰值，随着频率增加，曲线缓慢下降。 均方速度曲线和声功率曲线在全频段内变化趋势都非常一致，表明了单层壳在空气中的振动与声辐射响应在

11、全频段内都呈现单调关系。 辐射效率曲线符合一般规律，从低频开始逐渐升高，到达临界频率后逐渐趋近于1，成为有限辐射体。空气中计算结果 从图中可以看出，单层壳在水中的振动响应和空气中的振动响应低频段峰值位置不同，水中的振动响应要小一些，但整个曲线的趋势还是一致的。 从（b）中可以看出，虽然空气中的振动均方速度要比水中的振动均方速度大，但高频段水中的辐射声辐射响应比辐射声功率反而大。 从（c）中可得，从辐射效率的角度来揭示2中提到的现象。低频段，重流体负载的存在使得结构的辐射效率下降，而水中的均方振速比空气中的小，因此该频段水中的辐射声功率比空气中的小；在达到临界频率以上，流体负载影响很小，空气中和

12、水中辐射效率趋近一致，而高频段水中均方速度比空气中相差不多，且水中特性阻抗比空气的特性阻抗大得多，所以高频段水中的辐射声功率就会比空气中的大。双层壳算例空气中计算结果水中计算结果 从( a)和 (b)可以看到，双层壳在水中的振动响应和空气中的振动响应低频段峰值位置不同，水中的振动响应要小一些，但整个曲线的趋势还是一致的，这个特性跟单层壳体是一致的。 从 (d),(e)上可以看出，相比于空气中不同的是，在全频段范围内，水中的辐射声功率是随频率增加缓慢升高的。这是由于临界频率与c2c成正比，水中的临界频率比空气中大很多，如图3-11 (f)所示，水中还没达到临界频率，辐射效率一直在升高。结构参数对单双层圆柱壳全频段声振响应的影响 单层板算例 改变肋距的影响 改变壳体厚度的影响 双层板算例 改变实肋板肋距的影响 改变内壳厚度的影响 改变外壳厚度的影响 改变肋距的影响单层板 改变壳体厚度的影响 改变实肋板肋距的影响双层板 改变内壳厚度的影响 改变外壳厚度的影响VA One 计算结果与解析法的对比本章小结p 使用VA One 软件对计算模型在空气中和水中的振动声辐射响应进行计算，并得出了一些结论p 分析了