CN107562993A - 一种提高低频隔声性能的平板结构优化设计方法 - Google Patents
一种提高低频隔声性能的平板结构优化设计方法 Download PDFInfo
- Publication number
- CN107562993A CN107562993A CN201710610891.6A CN201710610891A CN107562993A CN 107562993 A CN107562993 A CN 107562993A CN 201710610891 A CN201710610891 A CN 201710610891A CN 107562993 A CN107562993 A CN 107562993A
- Authority
- CN
- China
- Prior art keywords
- boundary
- flat plate
- plate structure
- sound insulation
- insulation performance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Building Environments (AREA)
Abstract
本发明公开一种提高低频隔声性能的平板结构优化设计方法,包括建立可适用于任意边界条件的平板结构声传递损失的计算方程,设定各项已知参数;选择一个或多个结构低频隔声性能评价指标,并建立所选指标对应的计算方法,根据所选的指标设定平板结构的优化目标,如果步骤4所选指标多于1个,则设定各个指标的对应权重,选择优化算法,针对设定的已知参数、边界参数可调范围以及优化目标进行最优解搜索,以获得最优边界参数并形成平板结构。本发明可指导平板结构边界条件的优化设计,设计过程中利用优化算法的优势,快速搜索出边界条件最优解,从而实现平板结构的低频隔声性能最大化的效果。
Description
技术领域
本发明涉及一种提高低频隔声性能的平板结构优化设计方法。
背景技术
结构隔声降噪特性一直是建筑、机械等领域所关心的问题。为了 有效评价结构的隔声性能,人们提出了结构隔声性能的“单值评价” 指标。比如,ASTM组织提出的声传递等级“Sound Transmission Class” (STC),ISO国际化标准组织提出的计权声衰减指数“Weighted Sound Reduction Index”(Rw),以及我国建筑隔声评价标准中提出的“单 值评价量”等。然而,上述单值评价指标对低频范围考虑不足,难以 有效评价结构的低频隔声性能。为此,学者们最近又针对20-250Hz 的低频区间,提出了结构低频隔声性能的单值评价方法,如“LFSTC” 单值评价指标。该指标于2015年首次发表于英国《建筑与环境》杂 志(“Low frequency sound insulation analysis and evaluation of stiffenedbuilding structures”,Building and Environment, Vol.94(2),pp 802-809)。
平板结构是工程中最常见的结构之一,广泛应用于船舶与海洋工 程、航空航天工程、土木建筑工程、车辆工程等结构中。众所周知, 平板结构的边界安装条件对结构的低频隔声性能有重要影响,如果能 通过边界条件的优化设计来实现平板结构低频隔声性能的提高,将对 工程应用有重要意义,因为这样的设计只要求改变结构的边界条件, 不需要改动平板的自身参数(如长、宽、厚、材料等等),既能保持 平板结构的原有材料、尺寸和外观不受影响,又能有效提高其隔声性 能。然而,目前还没有一种优化设计方法,能有效指导平板结构边界 条件的优化设计,从而实现低频隔声性能的最大提升。
基于这样的背景,本发明致力于寻找一种更优的平板结构低频隔 声性能优化设计方法,用于指导工程人员快速确定平板结构的最佳边 界条件,从而实现平板结构的低频隔声性能最大化。该优化设计方法 将综合运用平板结构声传递损失预测模型以及低频隔声性能评价方 法,并引入优化算法的最优解搜索理论。
优化方法和技术已经是一种以数学为基础,用于求解各种工程问 题优化解的应用技术。所谓优化算法,其实就是一种搜索过程和规则, 它基于某种思想和机制,通过一定的途径或规则来得到满足用户要求 的问题的近似最优解。就优化机制与行为而分,目前常用的优化算法 主要可分为:经典算法、构造型算法、改进型算法、基于系统动态演 化的算法和混合型算法等等。选择哪种优化算法作为工具并不是唯一 的,各种算法都存在自身的优缺点。
本发明人利用了优化算法最优解搜索的思想,提出一种使用时并 不限定具体哪种算法,应用时可以根据实际问题的特点和需要进行选 择的平板结构低频隔声性能的优化设计方法。
发明内容
本发明的目的在于提供一种提高低频隔声性能的平板结构优化 设计方法,以指导平板结构边界条件的优化设计,通过利用优化算法 的优势可快速搜索出边界条件(边界参数)最优解,从而实现平板结 构的低频隔声性能最大化的效果。
为了实现上述目的,本发明采用如下技术方案:
一种提高低频隔声性能的平板结构优化设计方法,包括以下步 骤:
步骤1、建立可适用于任意边界条件的平板结构声传递损失的计 算方程:
其中,TL表示平板结构声传递损失;
{δ}为节点位移矢量;
ω为自振角频率;
为全局转换矩阵,可将单元负载压力转换为等效节点力;
为全局转换矩阵,可将节点位移转换为横向挠度;
{P0}为入射声波;
{M}、{C}和{K}分别为平板结构的整体质量矩阵、阻尼矩阵和 整体刚度矩阵;
{A1}和{A2}均为由边界配置法获得的方形矩阵,且分别由结构内、 外表面所处流体的介质特性所决定,{A1}和{A2}的矩阵元素可表示为: Ai(l,m)=ρiω2∫SG(xl,ξm)dS,(i=1,2),其中,ρ1、ρ2分别为结构 内、外表面的流体密度,xl和ξm分别表示接收点l和源点m的位置,G 为格林函数,S为单元面积;
步骤2、设定各项已知参数:包括入射声波、平板材料、阻尼因 子、尺寸、厚度以及平板结构内、外表面的流体密度和流体中声速;
步骤3、设定边界参数的可调范围;
步骤4、选择一个或多个结构低频隔声性能评价指标,并建立所 选指标对应的计算方法;
步骤5、根据步骤4所选的指标设定平板结构的优化目标,如果 步骤4所选指标多于1个,则根据实际需要设定各个指标的对应权重;
步骤6、选择优化算法,并将步骤1的平板结构声传递损失的计 算方程、步骤4所选的结构低频隔声性能评价指标及其计算方法,以 及所选择的优化算法三者相结合,针对步骤2设定的已知参数、步骤 3设定的边界参数可调范围以及步骤5设定的优化目标进行最优解搜 索,以获得满足要求的最优边界参数;
步骤7、根据步骤6的最优边界参数形成平板结构。
优选的,所述{M}、{C}和{K}由有限单元法获得,具体如下:{M} 由平板单元等效质量矩阵{Mp}e所组成,{K}由平板单元等效刚度矩阵 {Kp}e和支撑边界单元等效刚度矩阵{Kb}e所组成,且:
{Mp}e=ρph∫∫{N}T{Hp}{N}dxdy,
{Kp}e=∫∫{Bp}T{Dp}{Bp}dxdy,
其中,ρp为平板材料密度,h为平板厚度,{N}为单元形函数,{Bp} 为应变矩阵,{Dp}为抗弯刚度矩阵,{Hp}为关于h的对角矩阵函数, 为单元边界轮廓Γb的单位法向量;ktb和krb为支撑边界的参数, 分别代表结构支撑边界的横向刚度和转动刚度。
进一步的,所述的支撑边界的参数ktb和krb可为复数、常数或位 置函数;
当参数ktb和krb为复数时,其实部表示边界支撑的弹性特性,虚 部表示边界支撑的阻尼特性,若模拟弹性边界支撑则虚部值取0。
进一步的,所述平板结构内、外表面流体类型可以是所有流体中 的任意一种,且内、外表面的流体类型可以一样也可以不同。
优选的,所述结构低频隔声性能评价指标为LFSTCA或LFSTCO&N, 指标的计算方法具体如下:
1)将步骤1的平板结构声传递损失的计算方程计算得到的1/3 倍频程结构声传递损失数据绘制成一条频谱曲线;
2)在相同坐标比例下,根据步骤4所选择的结构低频隔声性能 评价指标,绘制出相应的1/3倍频程基准曲线;
3)将基准曲线与声传递损失曲线的横坐标互相重合,纵坐标中 基准曲线0dB与频谱曲线的一个整数坐标平齐;
4)将基准曲线向上移动,每步1dB,直至不利偏差之和尽量的 大,但不利偏差之和不超过24dB且单个最大不利偏差不超过8dB为 止;
5)基准曲线上250Hz对应的绘有声传递损失曲线的坐标值上纵 坐标的整分贝数作为所选指标的计算结果。
优选的,所示选择的优化算法为遗传算法或复合形算法。
采用上述方案后,本发明具有以下优点:
1、结合平板结构声传递损失计算模型、低频隔声性能评价指标 以及优化算法思想,可有效提高平板结构的低频隔声性能,指导平板 结构的优化设计。
2、对于给定的平板结构,可以在不改变平板自身参数(如长、 宽、厚、材料、外观等)的条件下,通过边界条件的优化设计,实现 平板结构的低频隔声性能最大化。既能保持平板结构的原有材料、尺 寸和外观不受影响,又能有效提高其隔声性能。
3、优化方法对平板结构内、外表面所处环境的流体类型没有限 制,适用于一般空气中平板结构隔声性能的优化设计,也适用于液体 中平板结构隔声性能的优化设计。
4、可指导任意边界条件下平板结构的低频隔声性能优化设计, 更符合实际需要。
5、优化过程引入了低频隔声性能评价指标,能更有效衡量和提 高平板结构的对低频噪声的隔绝能力,减小低频噪声对人体的干扰和 危害。
6、结合优化算法的优势,可同时对多个低频隔声性能评价指标 进行优化,更符合实际需要。
下面结合附图和具体实施例对本发明做进一步的说明。
附图说明
图1为本发明一种提高低频隔声性能的平板结构优化设计方法 的流程图;
图2为本发明的平板结构的示意图;
图3为本发明平板结构的正视图;
图4为本发明中遗传算法流程示意图。
具体实施方式
如图1所示,本发明揭示一种提高低频隔声性能的平板结构优化 设计方法,下面通过两个实施例进行说明。
实施例一,具体优化步骤如下:
步骤1、建立可适用于任意边界条件的平板结构声传递损失的计 算方程:
其中,TL为平板结构声传递损失,{δ}为节点位移矢量,ω为自 振角频率。为全局转换矩阵,可将单元负载压力转换为等效节 点力;为全局转换矩阵,可将节点位移转换为横向挠度;{P0}为 作用于平板内表面的入射声波,{M}、{C}和{K}分别为平板结构的 整体质量矩阵、阻尼矩阵和整体刚度矩阵,可由有限单元法获得,具 体如下:{M}由平板单元等效质量矩阵{Mp}e所组成,{K}由平板单元 等效刚度矩阵{Kp}e和支撑边界单元等效刚度矩阵{Kb}e所组成:
{Mp}e=ρph∫∫{N}T{Hp}{N}dxdy,
{Kp}e=∫∫{Bp}T{Dp}{Bp}dxdy,
其中,ρp为平板材料密度,h为平板厚度,{N}为单元形函数,{Bp} 为应变矩阵,{Dp}为抗弯刚度矩阵,{Hp}为关于h的对角矩阵函数, 为单元边界轮廓Γb的单位法向量;ktb和krb为支撑边界的参数, 分别代表结构支撑边界的横向刚度和转动刚度。这些边界参数若为复 数形式,其实部表示边界支撑的弹性特性,虚部表示边界支撑的阻尼 特性,若模拟弹性边界支撑则虚部值取0;这些边界参数也可以是常 数或位置函数,通过改变这些边界参数的形式可实现对任意边界条件 的模拟,本发明涉及平板结构可参见图2和图3。
上述的{A1}和{A2}是由边界配置法获得的方形矩阵,分别由结构 内、外表面所处流体的介质特性所决定,它们的矩阵元素可表示为: Ai(l,m)=ρiω2∫SG(xl,ξm)dS,(i=1,2),其中,ρ1、ρ2分别为结构 内、外表面的流体密度,xl和ξm分别表示接收点l和源点m的位置,G 为格林函数,S为单元面积。结构内、外表面的流体类型不限,可以 是空气、水等任何流体。内、外表面的流体类型可以一样(比如两面 都为空气)也可以不同(比如一面为空气、一面为水)
上述计算方程的各个参量可由入射声波、平板结构参数和流体介 质参数推导获得,其中平板结构参数包括平板材料、尺寸、厚度、边 界参数,流体介质参数包括结构内、外表面流体密度、流体中声速。 如图2和3所示为平板结构以及其内、外表面流体介质。
步骤2、设定各项已知参数:包括入射声波,平板材料、阻尼因 子、尺寸、厚度的平板参数,以及流体密度、流体中声速的平板结构 内表面和外表面流体特性参数。本例中的入射声波为垂直入射声波 (即入射声波入射方向与平板结构法向方向一致),其他具体参数值 如下表1所示:
表1已知参数
步骤3、设定边界参数的可调范围,如下表2所示:
表2边界参数
步骤4、选择符合要求的结构低频隔声性能评价指标,并建立所选指 标对应的计算方法。所选指标的个数不限,可以选一个也可以选多个。 本例选择LFSTCO&N指标,该指标计算方法建立如下:(1)将步骤1的平 板结构声传递损失的计算方程计算得到的1/3倍频程结构声传递损失数 据绘制成一条频谱曲线(简称为“声传递损失曲线”)。(2)在相同坐 标比例下,根据步骤4所选择的结构低频隔声性能评价指标,绘制出 LFSTCO&N指标相应的1/3倍频程基准曲线,其中,LFSTCA指标对应的基 准曲线的基准值见下表3。(3)将基准曲线与声传递损失曲线的横坐标 互相重合,纵坐标中基准曲线0dB与频谱曲线的一个整数坐标平齐。(4) 将基准曲线向上移动,每步1dB,直至不利偏差之和尽量的大,但不利偏 差之和不超过24dB且单个最大不利偏差不超过8dB为止。(5)此时基 准曲线上250Hz对应的绘有声传递损失曲线的坐标值上纵坐标的整分贝 数,就是所选指标的计算结果。
表3指标的计算结果
步骤5、根据步骤4所选的指标,设定平板结构的优化目标。如果步 骤4所选指标多于1个,则根据实际需要设定各个指标的对应权重。本 例所选指标为LFSTCO&N指标,优化目标为:实现LFSTCO&N指标值最大化。 优化目标所含指标个数为1,无需设定权重。
步骤6、选择符合要求的优化算法,并将步骤1的平板结构声传递损 失的计算方程、步骤4所选的结构低频隔声性能评价指标及其计算方法, 以及优化算法三者相结合,针对步骤2设定的已知参数、步骤3设定的 边界参数可调范围以及步骤5设定的优化目标进行最优解搜索,以获得满 足要求的最优边界参数。
具体使用哪种优化算法不做限制面,可根据实际需要进行选择,下 面分别以遗传算法和复合形算法为例说明本发明设计方法的具体步骤:
A、采用遗传算法作为优化算法:
遗传算法包括选择、交叉、变异三个主要操作算子,该算法初始时 随机产生N组解,每一组解叫一个个体,这多组解的集合叫做一个种群, 然后计算每个个体的适应度,选择操作使适应度大的个体有较大的复制 概率,能加快算法的收敛速度,交叉操作通过对两父代进行基因交换而 产生更优的个体,变异操作则能给群体带来新的基因,避免陷入局部最 优。通过这三个算子的操作,优化群体一代一代地不断进化,最终收敛 于最优状态。遗传算法包括如下运算步骤:
a、初始化:选择一个群体,即选择一个串或个体的集合bi,i=1, 2,...n,这个初始的群体也就是问题假设解的集合,通常以随机方法产 生串或个体的集合bi,i=1,2,...n;问题的最优解将通过这些初始假 设解进化而求出;
b、选择:根据适者生存原则选择下一代的个体,在选择时,以适应 度为选择原则,适应度准则体现了适者生存,不适应者淘汰的自然法则, 给出目标函数f,则f(bi)称为个体bi的适应度,适应度较高的个体, 繁殖下一代的数目较多;适应度较小的个体,繁殖下一代的数目较少; 甚至被淘汰,这样,就产生了对环境适应能力较强的后代,对于问题求 解角度来讲,就是选择出和最优解较接近的中间解;
c、交叉:对于选中用于繁殖下一代的个体,随机地选择两个个体的 相同位置,按交叉概率Pc在选中的位置实行交换,这个过程反映了随机信 息交换,目的在于产生新的基因组合,也即产生新的个体,交叉时,可 实行单点交叉或多点交叉,一般而言,交叉概率Pc取值为0.25—0.75;
d、变异:根据生物遗传中基因变异的原理,以变异概率Pm对某些个 体的某些位执行变异,在变异时,对执行变异的串的对应位求反,即把1 变为0,把0变为1,变异概率Pm与生物变异极小的情况一致,所以,Pm的取值较小,一般取0.01-0.2,单靠变异不能在求解中得到好处,但是, 它能保证算法过程不会产生无法进化的单一群体,因为在所有的个体一样时,交叉是无法产生新的个体的,这时只能靠变异产生新的个体,也 就是说,变异增加了全局优化的特质;
e、全局最优收敛:当最优个体的适应度达到给定的阈值,或者最优 个体的适应度和群体适应度不再上升时,则算法的迭代过程收敛、算法 结束;否则,用经过选择、交叉、变异所得到的新一代群体取代上一代 群体,并返回到步骤b,即在选择处继续循环执行,图4表示了遗传算法 的执行过程。
本发明利用遗传算法的进行优化设计具体包括如下步骤:
第一步、初始化控制参数:设置群体规模N,交叉概率Pc和变异概率 Pm;
第二步、在变量设定范围内随机产生初始种群;
第三步、对现有群体进行如下操作:
①计算群体中每个个体的适应度f(xi),i=1,2,...,N;
②按照赌轮机制实施选择操作,适应度大的个体的被选择的概率
③随机选出两个个体xi和xj作为父代,按照概率Pc进行交叉操作,产 生两个新的个体x′i和x′j,计算四个个体的适应度,选择其中最大的两个 个体;
④对交叉后的个体以概率Pm进行变异操作,接收变异后的新解,如 果满足收敛条件,则退出进化过程,否则转至(3)。
B、采用复合形算法作为优化算法:
本发明中复合形算法的具体步骤如下:
(1)初始化:选择一个集合Xk,其中k=1,2,...,r,r>n+1,n为优化 问题的维数;
(2)适应度函数为f;求得Xg; 得Xb;
(3)若d≤ε,则输出X*=Xg,f*=f(X*)停止,否则继续步骤 (4);
(4)求f(X0)=f0,若f0<f(Xb),则进入步骤(5), 否则进入步骤(6);
(5)若gi(X0)≤0(i=1,2,...,r),则Xb的值用X0代替,并转入步骤(2), 否则进入步骤5.1;
(5.1)取η=η0,t=0;
(5.2)求若gi(X0)≤0(i=1,2,...,r),则进入步骤5.3, 否则进入步骤5.4;
(5.3)求f(X),若f(X)<f(Xb),则Xb用X代替,并进入步骤 (2),否则进入步骤(6);
(5.4)令η=η0·η,t=t+1;若t≥N(N为认为设定的一个大数), 则转入步骤(6);否则转入步骤5.2;
(6)求若则Xb用代替,并转入步 骤(2),否则进入步骤(7);
(7)令Xk=βXg+(1-β)Xk,(k=1,2,...,r),并转入步骤(2)。
本例选择了遗传优化算法,并将步骤1的平板结构声传递损失的计 算方程、步骤4所选的结构低频隔声性能评价指标以及遗传优化算法三 者相结合,针对步骤2设定的已知参数、步骤3设定的边界参数的可调 范围以及步骤5设定的优化目标,经过最优解搜索后,得到如下表4所 示的边界参数优化结果:
表4边界参数优化结果
步骤7、根据步骤6的最优边界参数形成平板结构。
实施例一设计效果:实施例一中根据上表所示的边界参数优化结果 形成平板结构。其优化效果为LFSTCO&N=28dB,该结果比不经过优化,直 接取边界条件可调范围的上限值和下限值时的平板结构的LFSTCO&N指标 值,分别增加了17dB和7dB。
实施例二:
步骤1、建立可适用于任意边界条件的平板结构声传递损失的计 算方程;
步骤2、设定各项已知参数,其中入射声波为垂直入射声波(即 入射声波入射方向与平板结构法向方向一致),其他具体参数值如下 表5所示:
表5已知参数
步骤3、设定边界参数的可调范围,如下表所示:
表6边界参数
步骤4、同时选择LFSTCA和LFSTCO&N指标,并建立这两个指标对 应的计算方法;
步骤5、根据步骤4所选的指标,设定平板结构的优化目标为: 同时实现所选LFSTCA和LFSTCO&N的指标值最大化。优化目标所含指标 个数为2,设定二者权重一致,均为1。
步骤6、选择遗传优化算法,并将步骤1的平板结构声传递损失 的计算方程、步骤4所选的结构低频隔声性能评价指标及其计算方 法,以及遗传优化算法三者相结合,针对步骤2设定的已知参数、步 骤3设定的边界参数可调范围以及步骤5设定的优化目标,经过最优 解搜索后,得到如下表所示的边界参数优化结果:
步骤7、根据步骤6的最优边界参数形成平板结构。
上述说明示出并描述了本发明的优选实施例,应当理解本发明并 非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可 用于各种其他组合、修改和环境,并能够在本文发明构想范围内,通 过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行 的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利 要求的保护范围内。
Claims (6)
1.一种提高低频隔声性能的平板结构优化设计方法,其特征在于,包括以下步骤:
步骤1、建立可适用于任意边界条件的平板结构声传递损失的计算方程:
其中,TL表示平板结构声传递损失;
{δ}为节点位移矢量;
ω为自振角频率;
为全局转换矩阵,可将单元负载压力转换为等效节点力;
为全局转换矩阵,可将节点位移转换为横向挠度;
{P0}为入射声波;
{M}、{C}和{K}分别为平板结构的整体质量矩阵、阻尼矩阵和整体刚度矩阵;
{A1}和{A2}均为由边界配置法获得的方形矩阵,且分别由结构内、外表面所处流体的介质特性所决定,{A1}和{A2}的矩阵元素可表示为:Ai(l,m)=ρiω2∫SG(xl,ξm)dS,(i=1,2),其中,ρ1、ρ2分别为结构内、外表面的流体密度,xl和ξm分别表示接收点l和源点m的位置,G为格林函数,S为单元面积;
步骤2、设定各项已知参数:包括入射声波、平板材料、阻尼因子、尺寸、厚度以及平板结构内、外表面的流体密度和流体中声速;
步骤3、设定边界参数的可调范围;
步骤4、选择一个或多个结构低频隔声性能评价指标,并建立所选指标对应的计算方法;
步骤5、根据步骤4所选的指标设定平板结构的优化目标,如果步骤4所选指标多于1个,则设定各个指标的对应权重;
步骤6、选择优化算法,并将步骤1的平板结构声传递损失的计算方程、步骤4所选的结构低频隔声性能评价指标及其计算方法,以及所选择的优化算法三者相结合,针对步骤2设定的已知参数、步骤3设定的边界参数可调范围以及步骤5设定的优化目标进行最优解搜索,以获得满足要求的最优边界参数;
步骤7、根据步骤6的最优边界参数形成平板结构。
2.如权利要求1所述的一种提高低频隔声性能的平板结构优化设计方法,其特征在于:所述{M}、{C}和{K}由有限单元法获得,具体如下:{M}由平板单元等效质量矩阵{Mp}e所组成,{K}由平板单元等效刚度矩阵{Kp}e和支撑边界单元等效刚度矩阵{Kb}e所组成,且:
{Mp}e=ρph∫∫{N}T{Hp}{N}dxdy,
{Kp}e=∫∫{Bp}T{Dp}{Bp}dxdy,
<mrow>
<msub>
<mrow>
<mo>{</mo>
<msub>
<mi>K</mi>
<mi>b</mi>
</msub>
<mo>}</mo>
</mrow>
<mi>e</mi>
</msub>
<mo>=</mo>
<mo>&Integral;</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>k</mi>
<mrow>
<mi>t</mi>
<mi>b</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>{</mo>
<msub>
<mi>N</mi>
<mi>w</mi>
</msub>
<mo>}</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>{</mo>
<msub>
<mi>N</mi>
<mi>w</mi>
</msub>
<mo>}</mo>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>k</mi>
<mrow>
<mi>r</mi>
<mi>b</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>{</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>N</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>n</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>b</mi>
</msub>
</mrow>
</mfrac>
<mo>}</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>{</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>N</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>n</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>b</mi>
</msub>
</mrow>
</mfrac>
<mo>}</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>d&Gamma;</mi>
<mi>b</mi>
</msub>
<mo>,</mo>
</mrow>
其中,ρp为平板材料密度,h为平板厚度,{N}为单元形函数,{Bp}为应变矩阵,{Dp}为抗弯刚度矩阵,{Hp}为关于h的对角矩阵函数,为单元边界轮廓Γb的单位法向量;ktb和krb为支撑边界的参数,分别代表结构支撑边界的横向刚度和转动刚度。
3.如权利要求2所述的一种提高低频隔声性能的平板结构优化设计方法,其特征在于:所述的支撑边界的参数ktb和krb可为复数、常数或位置函数;
当参数ktb和krb为复数时,其实部表示边界支撑的弹性特性,虚部表示边界支撑的阻尼特性,若模拟弹性边界支撑则虚部值取0。
4.如权利要求1所述的一种提高低频隔声性能的平板结构优化设计方法,其特征在于:所述平板结构内、外表面流体类型可以是所有流体中的任意一种,且内、外表面的流体类型可以一样也可以不同。
5.如权利要求1所述的一种提高低频隔声性能的平板结构优化设计方法,其特征在于:所述结构低频隔声性能评价指标为LFSTCA或LFSTCO&N,指标的计算方法具体如下:
1)将步骤1的平板结构声传递损失的计算方程计算得到的1/3倍频程结构声传递损失数据绘制成一条频谱曲线;
2)在相同坐标比例下,根据步骤4所选择的结构低频隔声性能评价指标,绘制出相应的1/3倍频程基准曲线;
3)将基准曲线与声传递损失曲线的横坐标互相重合,纵坐标中基准曲线0dB与频谱曲线的一个整数坐标平齐;
4)将基准曲线向上移动,每步1dB,直至不利偏差之和尽量的大,但不利偏差之和不超过24dB且单个最大不利偏差不超过8dB为止;
5)基准曲线上250Hz对应的绘有声传递损失曲线的坐标值上纵坐标的整分贝数作为所选指标的计算结果。
6.如权利要求1所述的一种提高低频隔声性能的平板结构优化设计方法,其特征在于:所示选择的优化算法为遗传算法或复合形算法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710610891.6A CN107562993B (zh) | 2017-07-25 | 2017-07-25 | 一种提高低频隔声性能的平板结构优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710610891.6A CN107562993B (zh) | 2017-07-25 | 2017-07-25 | 一种提高低频隔声性能的平板结构优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107562993A true CN107562993A (zh) | 2018-01-09 |
CN107562993B CN107562993B (zh) | 2020-10-20 |
Family
ID=60974197
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710610891.6A Active CN107562993B (zh) | 2017-07-25 | 2017-07-25 | 一种提高低频隔声性能的平板结构优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107562993B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763641A (zh) * | 2018-04-20 | 2018-11-06 | 华侨大学 | 一种提高瞬态隔振性能的平板结构优化设计方法 |
CN109117502A (zh) * | 2018-07-04 | 2019-01-01 | 华侨大学 | 一种降低平板结构瞬态声辐射的优化设计方法 |
CN110285958A (zh) * | 2019-07-25 | 2019-09-27 | 重庆长安汽车股份有限公司 | 一种汽车零部件隔声性能的测试方法 |
CN111143962A (zh) * | 2019-11-12 | 2020-05-12 | 西南交通大学 | 一种隔声测评方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103699778A (zh) * | 2013-11-29 | 2014-04-02 | 华侨大学 | 一种考虑流体负载效应的平板结构优化设计方法 |
JP2015052264A (ja) * | 2014-10-27 | 2015-03-19 | 三菱重工業株式会社 | 遮音壁の設計方法および遮音壁 |
-
2017
- 2017-07-25 CN CN201710610891.6A patent/CN107562993B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103699778A (zh) * | 2013-11-29 | 2014-04-02 | 华侨大学 | 一种考虑流体负载效应的平板结构优化设计方法 |
JP2015052264A (ja) * | 2014-10-27 | 2015-03-19 | 三菱重工業株式会社 | 遮音壁の設計方法および遮音壁 |
Non-Patent Citations (2)
Title |
---|
CHUN HAOFENG等: "Measurement and estimation on sound insulation performance of the walls and floors in Taiwan"s dwelling-housings", 《2011 INTERNATIONAL CONFERENCE ON MULTIMEDIA TECHNOLOGY》 * |
欧达毅等: "加筋板结构隔声性能研究", 《2012全国环境声学学术会议论文集》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763641A (zh) * | 2018-04-20 | 2018-11-06 | 华侨大学 | 一种提高瞬态隔振性能的平板结构优化设计方法 |
CN108763641B (zh) * | 2018-04-20 | 2022-05-06 | 华侨大学 | 一种提高瞬态隔振性能的平板结构优化设计方法 |
CN109117502A (zh) * | 2018-07-04 | 2019-01-01 | 华侨大学 | 一种降低平板结构瞬态声辐射的优化设计方法 |
CN109117502B (zh) * | 2018-07-04 | 2022-09-02 | 华侨大学 | 一种降低平板结构瞬态声辐射的优化设计方法 |
CN110285958A (zh) * | 2019-07-25 | 2019-09-27 | 重庆长安汽车股份有限公司 | 一种汽车零部件隔声性能的测试方法 |
CN111143962A (zh) * | 2019-11-12 | 2020-05-12 | 西南交通大学 | 一种隔声测评方法及装置 |
CN111143962B (zh) * | 2019-11-12 | 2022-03-22 | 西南交通大学 | 一种隔声测评方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107562993B (zh) | 2020-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107562993B (zh) | 一种提高低频隔声性能的平板结构优化设计方法 | |
Wang et al. | Optimal design of acoustical sandwich panels with a genetic algorithm | |
CN111444643B (zh) | 一种基于神经网络的复合材料螺旋桨铺层角优化方法 | |
Du et al. | H∞ control for buildings with time delay in control via linear matrix inequalities and genetic algorithms | |
CN109885061B (zh) | 一种基于改进nsga-ⅱ的动力定位多目标优化方法 | |
CN101604019A (zh) | 一种海洋环境与声场不确实性的表征和传递的快速计算方法 | |
CN102054199A (zh) | 基于bp神经网络算法对涂层老化的分析方法 | |
CN114970366B (zh) | 一种功能梯度超材料结构优化设计方法及系统 | |
CN105608295A (zh) | 焦化炉压力的多目标遗传算法与rbf神经网络优化建模方法 | |
CN106886658A (zh) | 非线性迟滞动力学模型参数识别的粒子群新算法 | |
CN115841076A (zh) | 基于bp神经网络模型的浅海分层海底地声参数反演方法 | |
CN112926025A (zh) | 一种基于交叉熵遗传算法解决武器目标分配问题的方法 | |
CN114636360A (zh) | 五模冲击隐身复合点阵环状结构及其参数优化方法 | |
CN104062678B (zh) | 一种优化气枪阵列的方法和装置 | |
Blumrich et al. | Meteorologically induced variability of sonic-boom characteristics of supersonic aircraft in cruising flight | |
CN103699778B (zh) | 一种考虑流体负载效应的平板结构优化设计方法 | |
CN114841095B (zh) | 一种不可压缩流动的扰动域推进方法 | |
CN112287590B (zh) | 一种轨道交通减振设备布置和参数优化的方法 | |
CN113094999A (zh) | 一种基于神经网络的钢管混凝土纯弯构件承载力预测方法 | |
Fang et al. | Bandgap optimization and inverse design of labyrinth metamaterials for sound insulation | |
CN109117502B (zh) | 一种降低平板结构瞬态声辐射的优化设计方法 | |
Bo et al. | Nonlinear dynamic instability and multi-objective design optimization of the GPLR-FGP plate under biaxial impacts | |
CN103235879A (zh) | 基于小生境主从式并行遗传算法的双相介质参数反演方法 | |
CN108763641B (zh) | 一种提高瞬态隔振性能的平板结构优化设计方法 | |
CN108650043A (zh) | 一种高精度水声通信信道建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |