[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN115762685B - A sound-permeable layer and its optimization method for reducing the acoustic transmission loss at the interface of dissimilar materials - Google Patents

A sound-permeable layer and its optimization method for reducing the acoustic transmission loss at the interface of dissimilar materials Download PDF

Info

Publication number
CN115762685B
CN115762685B CN202211543461.4A CN202211543461A CN115762685B CN 115762685 B CN115762685 B CN 115762685B CN 202211543461 A CN202211543461 A CN 202211543461A CN 115762685 B CN115762685 B CN 115762685B
Authority
CN
China
Prior art keywords
layer
sound
acoustic transmission
transmission loss
dielectric
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.)
Active
Application number
CN202211543461.4A
Other languages
Chinese (zh)
Other versions
CN115762685A (en
Inventor
谷军杰
赵庆坤
曲绍兴
周昊飞
尹冰轮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202211543461.4A priority Critical patent/CN115762685B/en
Publication of CN115762685A publication Critical patent/CN115762685A/en
Application granted granted Critical
Publication of CN115762685B publication Critical patent/CN115762685B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • General Details Of Gearings (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a sound transmission layer for reducing acoustic transmission loss of a dissimilar material interface and an optimization method, and belongs to the field of optimization algorithms. The invention applies the concept of digital materials to disperse the whole optimization space into small material units, thereby parameterizing the whole design space and facilitating the optimization of the whole design space. The method adopts machine learning to obtain the prediction function of the relation between the digital material configuration and the average transmission loss, so that the prediction function is used for replacing finite element analysis in a genetic algorithm, and the calculation time is greatly saved. The method can design a more novel configuration, and is particularly suitable for designing the sound transmission layer in a low-frequency region.

Description

一种降低异种材料界面声学传递损失的透声层及优化方法A sound-permeable layer and its optimization method for reducing the acoustic transmission loss at the interface of dissimilar materials

技术领域technical field

本发明属于应用声学领域,具体涉及一种降低异种材料界面声学传递损失的透声层及优化方法。The invention belongs to the field of applied acoustics, and in particular relates to a sound-permeable layer and an optimization method for reducing the acoustic transmission loss at the interface of different materials.

背景技术Background technique

声学传递损失(STL),在介质无耗散的情况下,主要是由声学传导介质的阻抗不匹配造成的。声音在阻抗不匹配的异种材料界面处发生反射,导致透射声能小于入射声能,而由反射造成的能量损失即为声学传递损失。声学传递损失是声学传感器、换能器以及吸声材料设计等应用中的关键问题。例如,传感器与被测材料之间的传递损失会造成所测声压不准;同样,声学换能器得到的能量也会明显小于输入的声能;吸声材料(如水下消声瓦)入射界面处反射造成的升学传递损失会降低材料的吸声系数。降低声学传递损失的有效方法是在阻抗不匹配的两种介质间铺设一层声学覆盖层,即透声层。透声层通常由阻抗渐变的复合材料或结构构成,例如楔形结构、阻抗梯度变化的多层结构以及梯度颗粒夹杂的复合材料。这种阻抗渐变的复合材料使得阻抗在原阻抗不匹配的材料间均匀过渡,从而降低传递损失,因此又称为阻抗匹配层。Acoustic transfer loss (STL), in the case of a non-dissipative medium, is mainly caused by the impedance mismatch of the acoustically conducting medium. The sound is reflected at the interface of dissimilar materials with mismatched impedance, resulting in the transmitted sound energy being less than the incident sound energy, and the energy loss caused by reflection is the acoustic transmission loss. Acoustic transmission loss is a critical issue in applications such as the design of acoustic sensors, transducers, and sound-absorbing materials. For example, the transmission loss between the sensor and the measured material will cause the measured sound pressure to be inaccurate; similarly, the energy obtained by the acoustic transducer will be significantly less than the input sound energy; The transmission loss caused by reflection at the interface will reduce the sound absorption coefficient of the material. An effective way to reduce the acoustic transmission loss is to lay an acoustic covering layer, that is, a sound-permeable layer, between two media with mismatched impedances. The sound-permeable layer is usually composed of composite materials or structures with gradual impedance changes, such as wedge-shaped structures, multi-layer structures with gradient impedance changes, and composite materials with gradient particle inclusions. This kind of composite material with gradual impedance makes the impedance transition evenly between the original impedance mismatched materials, thereby reducing the transmission loss, so it is also called the impedance matching layer.

虽然这些透声层已被深入研究,但透声层的最优设计方案并不十分明朗。理论分析方面,现有的理论方法,包括传递矩阵法和非均匀传输线理论仅适用于一些简单的复合材料(例如,具有离散阻抗的多层板,具有指数型阻抗分布的介质)。在数值计算方面,传统数值优化计算方法通常受限于人的经验,优化只能对特定的复合材料模型进行几何参数优化,无法探索其它的可能的复合材料模型,因此无法得到真正意义上(全域)的最优解。另外,传统数值计算方法需要耗费大量的计算资源。Although these acoustic layers have been intensively studied, the optimal design of the acoustic layers is not very clear. In terms of theoretical analysis, existing theoretical methods, including transfer matrix method and non-uniform transmission line theory, are only applicable to some simple composite materials (for example, multilayer boards with discrete impedance, media with exponential impedance distribution). In terms of numerical calculation, the traditional numerical optimization calculation method is usually limited by human experience, optimization can only optimize the geometric parameters of a specific composite material model, and cannot explore other possible composite material models, so it is impossible to obtain the real (global ) optimal solution. In addition, traditional numerical calculation methods need to consume a lot of computing resources.

发明内容Contents of the invention

现有技术中存在的问题是:1)传统降低声学传递损失的透声层的优化设计方法只能针对特定复合材料模型,如楔型结构、梯度夹杂结构,无法探索所有可能的材料复合方式,导致无法遍历整个优化空间,以致无法得到全局最优解;2)传统数值优化计算方法(例如有限元、拓扑优化)耗时多的问题。本发明的目的在于解决现有技术中的上述技术问题,并提供一种降低异种材料界面声学传递损失的透声层及优化方法。The problems existing in the prior art are: 1) The traditional optimal design method of the sound-permeable layer to reduce the acoustic transmission loss can only be aimed at specific composite material models, such as wedge-shaped structures and gradient inclusion structures, and cannot explore all possible material composite modes. As a result, the entire optimization space cannot be traversed, so that the global optimal solution cannot be obtained; 2) Traditional numerical optimization calculation methods (such as finite element, topology optimization) consume a lot of time. The purpose of the present invention is to solve the above-mentioned technical problems in the prior art, and provide a sound-permeable layer and an optimization method for reducing the acoustic transmission loss at the interface of different materials.

本发明所采用的具体技术方案如下:The concrete technical scheme that the present invention adopts is as follows:

第一方面,本发明提供了一种降低异种材料界面声学传递损失的透声层优化方法,所述透声层设置于阻抗不匹配的第一介质层和第二介质层之间,该优化方法包括:In the first aspect, the present invention provides a method for optimizing a sound-transmitting layer for reducing the acoustic transmission loss of a dissimilar material interface. include:

S1、以透声层中单个周期性单元为优化目标,将周期性单元的纵剖面离散化为二进制矩阵,以矩阵的不同二进制元素值区分代表第一介质材料和第二介质材料,其中第一介质材料与第一介质层的材料相同或具有相似的阻抗,第二介质材料与第二介质层的材料相同或具有相似的阻抗;通过随机化生成二进制矩阵,批量生成随机的透声层数码材料模型,模型中的每个数码对应于一个材料单元,各单元的具体材料由对应的二进制矩阵中的元素值决定;S1. Taking a single periodic unit in the sound-transmitting layer as the optimization target, discretize the longitudinal section of the periodic unit into a binary matrix, and use different binary element values of the matrix to distinguish between the first medium material and the second medium material, where the first The medium material is the same as the material of the first medium layer or has a similar impedance, and the second medium material is the same as or has a similar impedance to the material of the second medium layer; a binary matrix is generated by randomization, and random sound-transmitting layer digital materials are generated in batches Model, each number in the model corresponds to a material unit, and the specific material of each unit is determined by the element value in the corresponding binary matrix;

S2、针对S1中生成的每个透声层数码材料模型,利用有限元模拟得到目标频带范围内的声学传递损失STL曲线以及声学传递损失平均值ATL;S2. For each digital material model of the sound-transmitting layer generated in S1, use finite element simulation to obtain the STL curve of the acoustic transmission loss and the average value of the acoustic transmission loss ATL within the target frequency band;

S3、将每个透声层数码材料模型的二进制矩阵转换为二进制向量v,以二进制向量v为训练样本的输入,以S2中计算得到的对应的声学传递损失平均值ATL为训练样本的标签,构建训练数据集;利用训练数据集进行机器学习训练,得到用于根据二进制向量v预测声学传递损失平均值ATL的预测函数;S3. Convert the binary matrix of each sound-permeable layer digital material model into a binary vector v, take the binary vector v as the input of the training sample, and use the corresponding acoustic transmission loss average value ATL calculated in S2 as the label of the training sample, Build a training data set; use the training data set to carry out machine learning training, and obtain a prediction function for predicting the average ATL of the acoustic transfer loss according to the binary vector v;

S4、以训练数据集中的训练样本作为初始种群,以所述预测函数为适应度函数,通过遗传算法进行优化,得到声学传递损失平均值ATL最小的最优二进制向量v;S4. Using the training samples in the training data set as the initial population, and using the prediction function as the fitness function, optimize through a genetic algorithm to obtain the optimal binary vector v with the smallest average value of acoustic transfer loss ATL;

S5、将所述最优二进制向量v转化成透声层数码材料模型得到最优构型,通过有限元模拟得到目标频带范围内的声学传递损失STL曲线以及声学传递损失平均值ATL,再根据有限元模拟结果判断该最优构型是否满足预期的优化要求,若满足则完成优化,若未满足则扩大训练数据集中的训练样本后重新进行机器学习训练生成新的预测函数,并再次进行遗传算法优化,直至最优构型满足优化要求;S5. Convert the optimal binary vector v into a digital material model of the sound-permeable layer to obtain the optimal configuration, and obtain the STL curve of the acoustic transmission loss and the average value of the acoustic transmission loss ATL within the target frequency band through finite element simulation, and then according to the finite The results of the meta-simulation judge whether the optimal configuration meets the expected optimization requirements. If so, the optimization is completed. If not, the training samples in the training data set are expanded, and then the machine learning training is performed again to generate a new prediction function, and the genetic algorithm is performed again. Optimize until the optimal configuration meets the optimization requirements;

S6、针对满足优化要求的最优构型,对两种介质材料之间的边界进行平滑,形成透声层的最终构型。S6. For the optimal configuration meeting the optimization requirements, smooth the boundary between the two dielectric materials to form the final configuration of the sound-transmitting layer.

作为上述第一方面的优选,所述第二介质层为水,第二介质材料采用水凝胶。As a preference of the first aspect above, the second medium layer is water, and the second medium material is hydrogel.

作为上述第一方面的优选,所述S5中,扩大训练数据集中的训练样本的做法为:在原始随机生成的训练数据集中增加根据人为总结经验所确定的经验构型。As a preference of the first aspect above, in the above S5, the practice of enlarging the training samples in the training data set is: adding the empirical configuration determined based on human experience in the original randomly generated training data set.

作为上述第一方面的优选,所述经验构型包括随机梯度夹杂结构、楔形梯度结构和随机楔形结构,每个经验构型均是一个透声层数码材料模型,也具有将周期性单元纵剖面离散化后形成的二进制矩阵;As a preference in the first aspect above, the empirical configuration includes a random gradient inclusion structure, a wedge-shaped gradient structure and a random wedge-shaped structure, each empirical configuration is a digital material model of a sound-transmitting layer, and also has a longitudinal section of a periodic unit The binary matrix formed after discretization;

所述随机梯度夹杂结构的透声层数码材料模型中,每个数码的材料类型随机生成,但沿着第二介质层所在侧到第一介质层所在侧的方向上,数码中出现第二介质材料的概率逐行递减,数码中出现第一介质材料的概率逐行梯度递增;In the digital material model of the sound-transmitting layer of the random gradient inclusion structure, the material type of each digital is randomly generated, but along the direction from the side where the second medium layer is located to the side where the first medium layer is located, the second medium appears in the digital The probability of the material decreases line by line, and the probability of the first medium material appearing in the digital gradually increases line by line;

所述楔形梯度结构的透声层数码材料模型中,第一介质材料整体呈连片且外轮廓相对规则的楔形分布,且沿着第二介质层所在侧到第一介质层所在侧的方向上,楔形宽度逐行递增。In the digital material model of the sound-transmitting layer with a wedge-shaped gradient structure, the first medium material is in the form of a continuous wedge-shaped distribution with a relatively regular outer contour, and along the direction from the side where the second medium layer is located to the side where the first medium layer is located , the wedge width increases row by row.

所述随机楔形结构的透声层数码材料模型中,第一介质材料整体呈不连片且外轮廓不规则的楔形分布,楔形内部随机出现第二介质材料,且沿着第二介质层所在侧到第一介质层所在侧的方向上,单行数码中出现第一介质材料的数码个数逐行递增。In the digital material model of the sound-transmitting layer with a random wedge-shaped structure, the first medium material is disjointed as a whole and has an irregular wedge-shaped distribution, and the second medium material randomly appears inside the wedge, and along the side where the second medium layer is located In the direction to the side where the first medium layer is located, the number of digits appearing in the first medium material in a single row of digits increases row by row.

作为上述第一方面的优选,所述随机梯度夹杂结构中,所述梯度递增的概率增长曲线为线性、多项式、对数或指数形式。As a preference of the first aspect above, in the stochastic gradient inclusion structure, the probability growth curve of increasing gradient is in linear, polynomial, logarithmic or exponential form.

作为上述第一方面的优选,所述楔形梯度结构中,楔形分布的第一介质材料的两侧轮廓线为直线、多项式曲线、对数曲线或指数曲线。As a preference of the first aspect above, in the wedge-shaped gradient structure, the contour lines on both sides of the wedge-shaped distribution of the first dielectric material are straight lines, polynomial curves, logarithmic curves or exponential curves.

作为上述第一方面的优选,进行所述机器学习训练时采用的机器学习模型为反向传播人工神经网络、支持向量机、卷积神经网络、线性回归。As a preference of the first aspect above, the machine learning model used in the machine learning training is a backpropagation artificial neural network, support vector machine, convolutional neural network, and linear regression.

作为上述第一方面的优选,所述机器学习模型采用BP神经网络。As a preference of the first aspect above, the machine learning model adopts BP neural network.

作为上述第一方面的优选,所述初始种群为训练数据集中声学传递损失平均值低于阈值的部分训练样本。As a preference of the first aspect above, the initial population is a part of training samples whose average value of acoustic transfer loss in the training data set is lower than a threshold.

作为上述第一方面的优选,所述的目标频带范围为0~10kHz。As a preference of the first aspect above, the range of the target frequency band is 0-10 kHz.

第二方面,本发明提供了一种根据上述第一方面任一方案所述透声层优化方法优化得到的透声层,该透声层为位于阻抗不匹配的第一介质层和第二介质层之间的双层结构,其中与第二介质层接触的一侧为阻抗跳跃层,与第一介质层接触的一侧为夹杂层,其中阻抗跳跃层完全由第一介质材料组成;而夹杂层以第二介质材料为主材,主材内部具有第一介质材料分布区域,但夹杂层中的第一介质材料分布区域由第二介质材料完全围合包裹。In a second aspect, the present invention provides a sound-transmitting layer optimized according to the sound-transmitting layer optimization method described in any one of the above-mentioned first aspects, the sound-transmitting layer is the first medium layer and the second medium located at impedance mismatches A double-layer structure between layers, wherein the side in contact with the second dielectric layer is an impedance jumping layer, and the side in contact with the first dielectric layer is an intercalation layer, wherein the impedance jumping layer is completely composed of the first dielectric material; and the inclusion The main material of the layer is the second dielectric material, and there is a distribution area of the first dielectric material inside the main material, but the distribution area of the first dielectric material in the interlayer is completely surrounded by the second dielectric material.

本发明相对于现有技术而言,具有以下有益效果:Compared with the prior art, the present invention has the following beneficial effects:

1)传统的优化方法只能针对固定构型进行少量的几何参数优化,而本发明方法应用数码材料的概念,将整个优化空间离散成一个个小的材料单元,从而对整个设计空间进行参数化,便于对整个设计空间进行优化。1) The traditional optimization method can only optimize a small number of geometric parameters for a fixed configuration, but the method of the present invention uses the concept of digital materials to discretize the entire optimization space into small material units, thereby parameterizing the entire design space , which facilitates the optimization of the entire design space.

2)相较于传统的有限元计算方法结合遗传算法的优化方法去遍历所有可能的数码材料构型,本发明方法首先利用机器学习得到数码材料构型与平均传递损失关系的预测函数,以便于在遗传算法中用预测函数代替有限元分析,大大节约了计算时间。2) Compared with the traditional finite element calculation method combined with the genetic algorithm optimization method to traverse all possible digital material configurations, the method of the present invention first uses machine learning to obtain the prediction function of the relationship between the digital material configuration and the average transmission loss, so as to facilitate In the genetic algorithm, the prediction function is used instead of the finite element analysis, which greatly saves the calculation time.

3)通过本发明方法可以设计出更加新颖的构型,尤其是在低频区间,例如水与其它介质之间的透声层设计,在(0-10kHz)范围内得到了阻抗跳越的复杂构型,相较于传统阻抗匹配构型,该构型使得平均声传递损失降低31%。3) More novel configurations can be designed by the method of the present invention, especially in the low frequency range, such as the design of the sound-permeable layer between water and other media, and a complex structure of impedance jumping is obtained in the (0-10kHz) range Type, compared with the traditional impedance matching configuration, this configuration reduces the average acoustic transmission loss by 31%.

附图说明Description of drawings

图1为透声层数码材料模型的生成方法示意图;Fig. 1 is a schematic diagram of the generation method of the digital material model of the sound-permeable layer;

图2为声学传递损失的有限元分析方法以及平均声学传递损失计算方法示意图;Fig. 2 is a schematic diagram of the finite element analysis method of the acoustic transmission loss and the calculation method of the average acoustic transmission loss;

图3为人工神经网络示意图;Fig. 3 is a schematic diagram of artificial neural network;

图4为遗传算法原理图;Fig. 4 is a schematic diagram of genetic algorithm;

图5为四种不同类型的数码材料构型生成效果,前一种为完全随机的二进制构型,后三种为基于人的经验,利用可以减小STL值的特征随机生成的经验构型;Figure 5 shows the generation effects of four different types of digital material configurations. The former is a completely random binary configuration, and the latter three are empirical configurations randomly generated based on human experience and using features that can reduce the STL value;

图6为本发明实施例中的透声层优化设计方法流程图;Fig. 6 is a flow chart of a sound-permeable layer optimization design method in an embodiment of the present invention;

图7为本发明实施例中优化设计得到的透声层构型。Fig. 7 shows the configuration of the sound-permeable layer obtained by the optimized design in the embodiment of the present invention.

具体实施方式Detailed ways

下面结合附图和具体实施方式对本发明做进一步阐述和说明。本发明中各个实施方式的技术特征在没有相互冲突的前提下,均可进行相应组合。The present invention will be further elaborated and illustrated below in conjunction with the accompanying drawings and specific embodiments. The technical features of the various implementations in the present invention can be combined accordingly on the premise that there is no conflict with each other.

为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图对本发明的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似改进,因此本发明不受下面公开的具体实施例的限制。本发明各个实施例中的技术特征在没有相互冲突的前提下,均可进行相应组合。In order to make the above objects, features and advantages of the present invention more comprehensible, specific implementations of the present invention will be described in detail below in conjunction with the accompanying drawings. In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. However, the present invention can be implemented in many other ways different from those described here, and those skilled in the art can make similar improvements without departing from the connotation of the present invention, so the present invention is not limited by the specific embodiments disclosed below. The technical features in the various embodiments of the present invention can be combined accordingly on the premise that there is no mutual conflict.

在本发明的描述中,需要理解的是,当一个元件被认为是“连接”另一个元件,可以是直接连接到另一个元件或者是间接连接即存在中间元件。相反,当元件为称作“直接”与另一元件连接时,不存在中间元件。In the description of the present invention, it will be understood that when an element is referred to as being "connected" to another element, it can be directly connected to the other element or indirectly, ie, intervening elements may be present. In contrast, when an element is referred to as being "directly" connected to another element, there are no intervening elements present.

在本发明的描述中,需要理解的是,术语“第一”、“第二”仅用于区分描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。In the description of the present invention, it should be understood that the terms "first" and "second" are only used for the purpose of distinction and description, and cannot be interpreted as indicating or implying relative importance or implicitly indicating the quantity of indicated technical features . Thus, the features defined as "first" and "second" may explicitly or implicitly include at least one of these features.

本发明中提供了一种降低异种材料界面声学传递损失的透声层优化方法,该透声层设置于阻抗不匹配的第一介质层A和第二介质层B之间。由于第一介质层A和第二介质层B之间阻抗存在差异,声音在阻抗不匹配的异种材料界面处发生反射,导致透射声能小于入射声能,而由反射造成的能量损失即为声学传递损失。而本发明设置透声层的目的是使得阻抗在原阻抗不匹配的材料间存在能够降低声学传递损失的过渡,降低传递损失。The present invention provides a method for optimizing the sound transmission layer for reducing the acoustic transmission loss at the interface of dissimilar materials. The sound transmission layer is arranged between the first dielectric layer A and the second dielectric layer B whose impedances do not match. Due to the difference in impedance between the first dielectric layer A and the second dielectric layer B, the sound is reflected at the interface of dissimilar materials with mismatched impedance, resulting in the transmitted sound energy being less than the incident sound energy, and the energy loss caused by reflection is the acoustic Transmission loss. The purpose of setting the sound-permeable layer in the present invention is to make the impedance transition between the materials whose original impedances do not match, which can reduce the acoustic transmission loss, and reduce the transmission loss.

需要说明的是,由于透声层一般是平面连续的,其中的结构呈现周期性的变化,因此本发明中无需优化整个透声层,而是仅需要将其中周期性变化的单元提取出来进行优化。由于透声层的构型(亦可称为结构)变化主要体现在纵剖面上,因此本发明中可针对透声层中单个周期性单元的纵剖面进行构型优化。It should be noted that, since the sound-transmitting layer is generally continuous in a plane, and the structure therein exhibits periodic changes, it is not necessary to optimize the entire sound-transmitting layer in the present invention, but only need to extract the periodically changing units for optimization. . Since the configuration (also called structure) change of the sound-transmitting layer is mainly reflected in the longitudinal section, the present invention can optimize the configuration of the longitudinal section of a single periodic unit in the sound-transmitting layer.

在本发明的一个较佳实施例中,上述降低异种材料界面声学传递损失的透声层优化方法包括如下步骤:In a preferred embodiment of the present invention, the method for optimizing the acoustic transmission layer for reducing the acoustic transmission loss at the interface of dissimilar materials includes the following steps:

S1、以透声层中单个周期性单元为优化目标,将周期性单元的纵剖面离散化为二进制矩阵,以矩阵的不同二进制元素值区分代表第一介质材料和第二介质材料,其中第一介质材料与第一介质层的材料相同或具有相似的阻抗,第二介质材料与第二介质层的材料相同或具有相似的阻抗;通过随机化生成二进制矩阵,批量生成随机的透声层数码材料模型,模型中的每个数码对应于一个材料单元,各单元的具体材料由对应的二进制矩阵中的元素值决定。S1. Taking a single periodic unit in the sound-transmitting layer as the optimization target, discretize the longitudinal section of the periodic unit into a binary matrix, and use different binary element values of the matrix to distinguish between the first medium material and the second medium material, where the first The medium material is the same as the material of the first medium layer or has a similar impedance, and the second medium material is the same as or has a similar impedance to the material of the second medium layer; a binary matrix is generated by randomization, and random sound-transmitting layer digital materials are generated in batches Each number in the model corresponds to a material unit, and the specific material of each unit is determined by the element value in the corresponding binary matrix.

需要说明的,上述二进制矩阵是基于材料声透射问题的物理模型构建的,如图1所示,透声层位于介质材料A和介质材料B之间,该透声层的一个周期性单元纵剖面可离散化为0-1二进制矩阵。二进制矩阵中每个矩阵元素的值都是一个二进制值,即0或1,本发明中可以用1代表第一介质材料,用0代表第二介质材料,也可以用0代表第一介质材料,用1代表第二介质材料,对此不做限定。基于二进制矩阵,可以将矩阵中的元素按照其二进制元素值对应映射形成透声层中一个填充指定介质材料的单元,本发明将将这个单元称为数码,由此形成用于进行数值化模拟的透声层数码材料模型。若透声层厚度为H,二进制矩阵为n×m阶矩阵,则透声层数码材料模型中的一个数码的高度为H/m。It should be noted that the above binary matrix is constructed based on the physical model of the material sound transmission problem. As shown in Figure 1, the sound transmission layer is located between the medium material A and the medium material B, and the longitudinal section of a periodic unit of the sound transmission layer It can be discretized into a 0-1 binary matrix. The value of each matrix element in the binary matrix is a binary value, namely 0 or 1. In the present invention, 1 can be used to represent the first dielectric material, 0 can be used to represent the second dielectric material, and 0 can also be used to represent the first dielectric material. 1 represents the second dielectric material, which is not limited. Based on the binary matrix, the elements in the matrix can be correspondingly mapped according to their binary element values to form a unit filled with a specified medium material in the sound-transmitting layer. The present invention will refer to this unit as a digital, thus forming a numerical simulation. Digital material model of the acoustic layer. If the thickness of the sound-transmitting layer is H, and the binary matrix is a matrix of order n×m, then the height of a digit in the digital material model of the sound-transmitting layer is H/m.

需要说明的是,第一介质材料以与第一介质层的材料完全相同为佳,同样的第二介质材料以与第二介质层的材料为佳,但如果两者的材料阻抗差异不大,则也可以采用两种不同的材料。例如,如果第二介质层为水,由于水为液态无法用于构建层体,因此第二介质材料可以采用水凝胶替换,因为水凝胶与水的阻抗差距很小。第一介质材料也需要根据实际模拟的第一介质层材质而定,可采用橡胶、聚氨酯(PU)等。It should be noted that the first dielectric material is preferably completely the same as the material of the first dielectric layer, and the same second dielectric material is preferably the same as the material of the second dielectric layer, but if the material impedance difference between the two is not large, Two different materials can then also be used. For example, if the second dielectric layer is water, since water cannot be used to construct the layer due to its liquid state, the second dielectric material can be replaced by hydrogel, because the impedance difference between hydrogel and water is very small. The material of the first medium also needs to be determined according to the material of the first medium layer to be actually simulated, and rubber, polyurethane (PU) and the like can be used.

S2、针对S1中生成的每个透声层数码材料模型,利用有限元模拟得到目标频带范围内的声学传递损失(STL)曲线以及声学传递损失平均值(ATL)。S2. For each digital material model of the sound-transmitting layer generated in S1, use finite element simulation to obtain an acoustic transmission loss (STL) curve and an average acoustic transmission loss (ATL) within the target frequency band.

需要说明的是,通过有限元模拟来计算目标频带范围内的声学传递损失(STL)曲线属于现有技术。如图2所示,展示了声学传递损失的有限元分析方法以及平均声学传递损失计算方法流程,其中有限元模型的构建是关键步骤。基于声透射问题的物理模型构建有限元模型时,需要上下分别设置硬声场边界、低反射边界,然后两侧设置对称性边界,在四个边界为何范围内从上到下依次设置完美匹配层、第二介质材料B层、压力声学与固体耦合边界、透声层、第一介质材料A层。其中透声层即可导入前述的透声层数码材料模型来构建。在关注的频带范围内,即在目标频带范围内,可通过有限元模拟对声波的入射过程进行模拟,得到不同频道处的STL值构成的STL曲线,将目标频带范围内所有频道对应的STL值进行平均,即可得到ATL值。例如图2中,其关注的目标频带范围为0~10kHz的低频区间,因此将这个区间分成了10个频带,对10个频带的STL值进行平均得到对应的ATL值。It should be noted that the calculation of the acoustic transmission loss (STL) curve within the target frequency band through finite element simulation belongs to the prior art. As shown in Figure 2, the flow of the finite element analysis method of the acoustic transmission loss and the calculation method of the average acoustic transmission loss is shown, and the construction of the finite element model is a key step. When building a finite element model based on the physical model of the sound transmission problem, it is necessary to set up a hard sound field boundary, a low reflection boundary, and then set a symmetrical boundary on both sides, and set a perfect matching layer from top to bottom within the range of the four boundaries. The B layer of the second dielectric material, the pressure acoustics and solid coupling boundary, the sound-permeable layer, and the A layer of the first dielectric material. The acoustic layer can be constructed by importing the aforementioned digital material model of the acoustic layer. Within the frequency band of interest, that is, within the target frequency band, the incident process of sound waves can be simulated through finite element simulation to obtain the STL curve composed of the STL values at different channels, and the STL values corresponding to all channels within the target frequency band By averaging, the ATL value can be obtained. For example, in Figure 2, the target frequency range of interest is the low-frequency range of 0 to 10 kHz, so this range is divided into 10 frequency bands, and the corresponding ATL value is obtained by averaging the STL values of the 10 frequency bands.

上述每一个透声层数码材料模型,均可通过有限元模拟得到对应的ATL值,由此即可用于作为机器学习的样本。因此,批量生成透声层数码材料模型时,其数量需要保证机器学习有足够的样本进行训练。For each digital material model of the sound-transmitting layer mentioned above, the corresponding ATL value can be obtained through finite element simulation, and thus can be used as a sample for machine learning. Therefore, when generating digital material models of sound-transmitting layers in batches, the quantity needs to ensure that machine learning has enough samples for training.

S3、将每个透声层数码材料模型的二进制矩阵转换为二进制向量v,以二进制向量v为训练样本的输入,以S2中计算得到的对应的声学传递损失平均值ATL为训练样本的标签,构建训练数据集;利用训练数据集进行机器学习训练,得到用于根据二进制向量v预测声学传递损失平均值ATL的预测函数。S3. Convert the binary matrix of each sound-permeable layer digital material model into a binary vector v, take the binary vector v as the input of the training sample, and use the corresponding acoustic transmission loss average value ATL calculated in S2 as the label of the training sample, Construct a training data set; use the training data set to carry out machine learning training, and obtain a prediction function for predicting the average ATL of the acoustic transfer loss according to the binary vector v.

需要说明的是,二进制矩阵转换为二进制向量主要是为了便于输入机器学习模型。本发明具体可采用的机器学习模型不限,反向传播人工神经网络、支持向量机、卷积神经网络、线性回归等均可作为本发明中的机器学习模型。机器学习模型的训练属于现有技术,对此不再展开描述。It should be noted that the conversion of binary matrix to binary vector is mainly for the convenience of inputting into machine learning models. The machine learning model that can be used in the present invention is not limited, and the backpropagation artificial neural network, support vector machine, convolutional neural network, linear regression, etc. can all be used as the machine learning model in the present invention. The training of the machine learning model belongs to the prior art, and no further description will be given here.

如图3所示,展示了一个实施例中采用的以BP神经网络形式的机器学习模型来构建预测函数的示意,其输入层维度为的二进制向量v维度,即1*400,隐藏层有三层,维度分别是10-5-5,输出层维度为1,即ATL值。由此,该机器学习模型训练完毕后得到的模型可视为是一个根据二进制向量v预测声学传递损失平均值ATL的预测函数f,即ATL=f(v)。As shown in Figure 3, it shows a schematic diagram of constructing a prediction function using a machine learning model in the form of a BP neural network in an embodiment. The dimension of the input layer is the dimension of the binary vector v, which is 1*400, and the hidden layer has three layers. , the dimensions are 10-5-5, and the output layer dimension is 1, which is the ATL value. Therefore, the model obtained after the training of the machine learning model can be regarded as a prediction function f for predicting the average ATL of the acoustic transfer loss according to the binary vector v, that is, ATL=f(v).

S4、以训练数据集中的训练样本作为初始种群,以上述预测函数ATL=f(v)为适应度函数,通过遗传算法进行优化,得到声学传递损失平均值ATL最小的最优二进制向量v。S4. Using the training samples in the training data set as the initial population, and using the prediction function ATL=f(v) as the fitness function, optimize through the genetic algorithm to obtain the optimal binary vector v with the smallest average ATL of the acoustic transmission loss.

本发明中的遗传算法(Genetic Algorithm,简称GA)是一种自适应随机搜索启发式算法,广泛应用于复杂函数系统优化、机器学习、系统识别、故障诊断、分类系统等领域中。如图4所示,GA算法利用计算机仿真运算,将问题的求解过程转换成类似生物进化中的染色体基因的交叉、变异等过程,其中每一轮迭代需要利用适应度函数来区分群体中个体的好坏,并基于适应值对个体进行筛选,以保证适应值好的个体有机会在下一代中产生更多的子个体。因此,在使用GA求解具体问题时,适应值函数的选择对算法的收敛性以及收敛速度的影响较大。本发明中,上述预测函数ATL=f(v)记为遗传算法优化过程中的适应度函数。最终优化的目标是最小化该适应度函数,即保证声学传递损失平均值ATL最小。当优化得到上述最优二进制向量v后,即可重新将二进制向量v转化为二进制矩阵形式,进而形成最优透声层数码材料模型。The genetic algorithm (Genetic Algorithm, referred to as GA) in the present invention is an adaptive random search heuristic algorithm, which is widely used in complex function system optimization, machine learning, system identification, fault diagnosis, classification system and other fields. As shown in Figure 4, the GA algorithm uses computer simulation operations to transform the problem-solving process into a process similar to the crossover and mutation of chromosomal genes in biological evolution. Individuals are screened based on fitness values to ensure that individuals with good fitness values have the opportunity to produce more sub-individuals in the next generation. Therefore, when using GA to solve specific problems, the choice of fitness value function has a great influence on the convergence and convergence speed of the algorithm. In the present invention, the above prediction function ATL=f(v) is recorded as the fitness function in the genetic algorithm optimization process. The final optimization goal is to minimize the fitness function, that is, to ensure the minimum average value of acoustic transmission loss ATL. After the above-mentioned optimal binary vector v is obtained through optimization, the binary vector v can be converted into a binary matrix form again, and then the optimal sound-transmitting layer digital material model can be formed.

为了加快遗传算法的收敛,在选择初始种群时,可对训练数据集中的训练样本进行筛选,去除ATL值过高的部分样本,保留ATL值较小的部分样本。In order to speed up the convergence of the genetic algorithm, when selecting the initial population, the training samples in the training data set can be screened to remove some samples with high ATL values and retain some samples with small ATL values.

S5、将上述最优二进制向量v转化成透声层数码材料模型得到最优构型,通过有限元模拟得到目标频带范围内的声学传递损失STL曲线以及声学传递损失平均值ATL,再根据有限元模拟结果判断该最优构型是否满足预期的优化要求,若满足则完成优化,若未满足则扩大训练数据集中的训练样本后重新进行机器学习训练生成新的预测函数,并再次进行遗传算法优化,直至最优构型满足优化要求。S5. Convert the above-mentioned optimal binary vector v into the digital material model of the sound-permeable layer to obtain the optimal configuration, and obtain the acoustic transmission loss STL curve and the average value of the acoustic transmission loss ATL within the target frequency band through finite element simulation, and then according to the finite element The simulation results judge whether the optimal configuration meets the expected optimization requirements, and if so, the optimization is completed; if not, the training samples in the training data set are expanded, and then the machine learning training is performed again to generate a new prediction function, and the genetic algorithm is optimized again , until the optimal configuration meets the optimization requirements.

需要说明的是,本步骤中的有限元模拟与上述S2中是一样的,此处不再赘述。其目的是为了判断该最优构型是否满足预期的优化要求,此处预期的优化要求需要根据实际情况进行确定,可预先设定希望透声层最终满足的STL曲线以及需要满足的ATL值,具体取值此处无需限定。It should be noted that the finite element simulation in this step is the same as that in S2 above, and will not be repeated here. Its purpose is to judge whether the optimal configuration meets the expected optimization requirements. The expected optimization requirements here need to be determined according to the actual situation. The STL curve that the sound-transmitting layer is expected to finally meet and the ATL value that needs to be satisfied can be preset in advance. The specific value is not limited here.

但是,如果有限元模拟后发现该最优构型并不能满足预期的优化要求,则可能是初始批量生成的透声层数码材料模型样本数量不够,寻优过程陷入了局部最优解,而没有找到全局最优解。因此,这种情况下需要扩大训练数据集中的训练样本,再重新进行机器学习训练生成新的预测函数f,并再次进行上述S4和S5的遗传算法优化,然后再次判断是否达到预期。不断循环该过程,直至最优构型满足预期的优化要求为止。However, if it is found after the finite element simulation that the optimal configuration cannot meet the expected optimization requirements, it may be that the number of samples of the digital material model of the sound-transmitting layer generated in the initial batch is not enough, and the optimization process has fallen into a local optimal solution, and there is no Find the global optimal solution. Therefore, in this case, it is necessary to expand the training samples in the training data set, and then perform machine learning training again to generate a new prediction function f, and perform the genetic algorithm optimization of S4 and S5 again, and then judge whether it meets expectations again. This process is repeated continuously until the optimal configuration meets the expected optimization requirements.

需要注意的是,在扩大训练数据集中的训练样本时,可以按照S1中的做法继续采用随机生成方式生成透声层数码材料模型并通过有限元模拟获取标签,但由于这些透声层数码材料模型对应的二进制矩阵是随机生成的,可能并不一定符合降低ATL所需的构型形式。It should be noted that when expanding the training samples in the training data set, the method of random generation can continue to be used to generate the digital material model of the sound-transmitting layer and obtain the labels through finite element simulation. The corresponding binary matrix is randomly generated, and may not necessarily conform to the configuration form required to reduce ATL.

因此作为本发明实施例的一种优选方式,在扩大训练数据集中的训练样本时,采用的做法不是随机生成,而是预先根据人为总结的经验生成部分经验构型,再将这些经验构型进行有限元模拟得到ATL值后,以样本形式加入到在原始随机生成的训练数据集中。此处增加的经验构型形式不限,可根据现有技术中的相关研究或者实际试验进行选择。本实施例中,给出了三种经验构型,包括随机梯度夹杂结构、楔形梯度结构和随机楔形结构,每个经验构型均是一个透声层数码材料模型,也具有将周期性单元纵剖面离散化后形成的二进制矩阵,二进制矩阵的维度与随机生成的二进制矩阵相同。下面对这三种经验构型的具体形式进行详细描述。Therefore, as a preferred method of the embodiment of the present invention, when expanding the training samples in the training data set, the method adopted is not to randomly generate, but to generate some empirical configurations in advance based on the experience summarized by humans, and then perform these empirical configurations. After the ATL value is obtained by finite element simulation, it is added to the original randomly generated training data set in the form of a sample. The empirical configuration form added here is not limited, and can be selected according to relevant research or actual experiments in the prior art. In this example, three empirical configurations are given, including stochastic gradient inclusion structure, wedge-shaped gradient structure and random wedge-shaped structure. The binary matrix formed after discretization of the profile, the dimension of the binary matrix is the same as that of the randomly generated binary matrix. The specific forms of these three empirical configurations are described in detail below.

如图5中的第一张图所示,以完全随机生成的透声层数码材料模型为对照,由于其二进制矩阵是随机生成的,因此每个数码中的介质材料类型也是杂乱无序的。As shown in the first picture in Figure 5, compared with the completely randomly generated digital material model of the sound-transmitting layer, since its binary matrix is randomly generated, the type of medium material in each digital is also disorderly.

如图5中的第二张图所示,随机梯度夹杂结构的透声层数码材料模型中,每个数码的材料类型随机生成,即每个数码随机设定是第一介质材料还是第二介质材料,但是宏观上而言,沿着第二介质层所在侧到第一介质层所在侧的方向上,数码中出现第二介质材料的概率逐行递减,数码中出现第一介质材料的概率逐行梯度递增。由此宏观来看,沿着第二介质层所在侧到第一介质层所在侧的方向上,数码中出现第二介质材料的概率逐行递减,每一行数码中出现第一介质材料数码的数量不断增加,但其出现的位置不定。在该随机梯度夹杂结构中,由于数码中出现第一介质材料的概率需要逐行梯度递增,因此实际实现时需要有设置相应的概率增长曲线来控制每一行出现第一介质材料的概率,可选的概率增长曲线为线性、多项式、对数或指数形式,对此不做限定,可尽可能丰富地生成不同的模型以便于全局寻优。As shown in the second picture in Figure 5, in the digital material model of the sound-permeable layer with random gradient inclusion structure, the material type of each digital is randomly generated, that is, each digital is randomly set to be the first medium material or the second medium However, macroscopically speaking, along the direction from the side where the second medium layer is to the side where the first medium layer is located, the probability of the second medium material appearing in the digits decreases row by row, and the probability of the first medium material appearing in the digits gradually decreases. Row gradient increments. From a macro point of view, along the direction from the side where the second medium layer is located to the side where the first medium layer is located, the probability of the second medium material appearing in the numbers decreases line by line, and the number of the first medium material numbers appearing in each line of numbers Constantly increasing, but the location of its occurrence is uncertain. In this stochastic gradient inclusion structure, since the probability of the first dielectric material appearing in the digital needs to increase step by row, it is necessary to set a corresponding probability growth curve to control the probability of the first dielectric material appearing in each row in actual implementation. Optional The probability growth curve of is in linear, polynomial, logarithmic or exponential form, which is not limited, and different models can be generated as richly as possible to facilitate global optimization.

如图5中的第三张图所示,楔形梯度结构的透声层数码材料模型中,第一介质材料整体呈连片且外轮廓相对规则的楔形分布,且沿着第二介质层所在侧到第一介质层所在侧的方向上,楔形宽度逐行递增。这种构型相对于前两种构型,其中的介质材料分布不再杂乱无序,而是呈现相对规则的状态。在该楔形梯度结构中,楔形分布的第一介质材料的两侧轮廓线可以是根据经验指定的直线、多项式曲线、对数曲线或指数曲线,对此不做限定,可尽可能丰富地生成不同的模型以便于全局寻优。As shown in the third picture in Figure 5, in the digital material model of the sound-transmitting layer with a wedge-shaped gradient structure, the first medium material is in a continuous piece with a relatively regular wedge-shaped distribution on the outside, and along the side where the second medium layer is located In the direction to the side where the first dielectric layer is located, the width of the wedge increases row by row. Compared with the previous two configurations, the distribution of dielectric materials in this configuration is no longer chaotic, but relatively regular. In the wedge-shaped gradient structure, the contour lines on both sides of the wedge-shaped distribution of the first dielectric material can be straight lines, polynomial curves, logarithmic curves or exponential curves specified based on experience, which is not limited, and can be generated as richly as possible. model for global optimization.

如图5中的第四张图所示,随机楔形结构的透声层数码材料模型中,第一介质材料整体呈不连片且外轮廓不规则的楔形分布,楔形内部随机出现第二介质材料,且沿着第二介质层所在侧到第一介质层所在侧的方向上,单行数码中出现第一介质材料的数码个数逐行递增。这种随机楔形结构相对于上一种楔形梯度结构而言,其区别在于外轮廓不再根据经验指定,而是随机地生成,但总体依然呈楔形,但由于是随机生成的因此楔形区域内部也可能出现第二介质材料,第一介质材料的楔形分布不再完全连续。As shown in the fourth picture in Figure 5, in the digital material model of the sound-transmitting layer with a random wedge-shaped structure, the first medium material is distributed in a wedge-shaped discontinuous and irregular outer contour as a whole, and the second medium material randomly appears inside the wedge , and along the direction from the side where the second medium layer is located to the side where the first medium layer is located, the number of numbers appearing in the first medium material in a single line of numbers increases line by line. Compared with the previous wedge-shaped gradient structure, the difference of this random wedge-shaped structure is that the outer contour is no longer specified according to experience, but randomly generated, but the overall shape is still wedge-shaped, but because it is randomly generated, the inside of the wedge-shaped area is also A second dielectric material may appear where the wedge-shaped distribution of the first dielectric material is no longer completely continuous.

S6、针对满足优化要求的最优构型,对两种介质材料之间的边界进行平滑,形成透声层的最终构型。S6. For the optimal configuration meeting the optimization requirements, smooth the boundary between the two dielectric materials to form the final configuration of the sound-transmitting layer.

需要说明的是,该步骤中进行边界平滑的目的是考虑材料实际加工的技术需求,因为边界过于离散的材料边界在实际的加工过程中不具有可实现性。而对两种介质材料之间的边界进行平滑后,可以保证其边界相对平滑,有利于实现透声层的加工。It should be noted that the purpose of boundary smoothing in this step is to consider the technical requirements of actual material processing, because material boundaries with too discrete boundaries are not achievable in actual processing. After smoothing the boundary between the two dielectric materials, the boundary can be guaranteed to be relatively smooth, which is beneficial to the processing of the sound-transmitting layer.

另外,当得到透声层的最终构型后,在实际加工透声层之前,亦可对该最终构型进行一定程度的简化,去除不符合加工工艺要求的部分结构特征,或者可进一步对优化得到的最终构型中的特征进行机理研究,确认哪些结构特征属于关键特征,哪些属于非关键特征,进而实现对实际加工构型的进一步调整。In addition, after the final configuration of the sound-transmitting layer is obtained, before the actual processing of the sound-transmitting layer, the final configuration can also be simplified to a certain extent to remove some structural features that do not meet the requirements of the processing technology, or further optimize The features in the obtained final configuration are subjected to mechanism research to confirm which structural features are key features and which are non-key features, so as to further adjust the actual processing configuration.

下面将上述S1~S6所示的降低异种材料界面声学传递损失的透声层优化方法应用与一个具体的实例中,以展示其具体实现过程和技术效果。In the following, the method for optimizing the sound transmission layer for reducing the acoustic transmission loss at the interface of dissimilar materials shown in S1-S6 above is applied to a specific example to demonstrate its specific implementation process and technical effect.

实施例Example

在本实施例中,降低异种材料界面声学传递损失的透声层优化方法的过程如图6所示,具体步骤如下:In this embodiment, the process of optimizing the sound transmission layer for reducing the acoustic transmission loss at the interface of different materials is shown in Figure 6, and the specific steps are as follows:

第一步:针对阻抗不匹配的两种介质A和B以及给定的透声层厚度H,基于数码材料的概念批量生成随机的透声层数码材料模型,透声层为第一介质材料A和第二介质材料B构成的复合材料。本实施例中,透声层两侧的介质分别为水和橡胶,因此组成透声层的介质材料分别选择为丙烯酰胺水凝胶和聚氨酯(PU)。生成方法如下:首先,随机生成一系列n×m阶二进制矩阵M,矩阵中1和0分别代表介质材料A和介质材料B;然后,根据每个二进制矩阵M中1、0和A、B的对应关系,生成一个透声层的数码材料模型,其中每个数码代表一个材料单元,每个材料单元的高度为H/m。Step 1: For two media A and B with mismatched impedance and a given thickness H of the sound-transmitting layer, generate random digital material models of the sound-transmitting layer in batches based on the concept of digital materials, and the sound-transmitting layer is the first medium material A A composite material composed of a second dielectric material B. In this embodiment, the media on both sides of the sound-transmitting layer are water and rubber respectively, so the media materials constituting the sound-transmitting layer are respectively selected as acrylamide hydrogel and polyurethane (PU). The generation method is as follows: First, randomly generate a series of n×m-order binary matrices M, 1 and 0 in the matrix represent dielectric material A and dielectric material B respectively; then, according to the According to the corresponding relationship, a digital material model of the sound-permeable layer is generated, where each number represents a material unit, and the height of each material unit is H/m.

第二步:针对批量生成的各透声层数码材料构型,利用有限元软件分别模拟计算其所关心的频带范围内的声学传递损失(STL),并计算声学传递损失在频带内的平均值ATL。在本实施例中,关心的频带范围为0~10kHz的低频区间,该频带区间分成10个频带,对STL曲线上10个频带的STL值进行平均,即可得到对应的ATL值。Step 2: For the digital material configuration of each sound-transmitting layer generated in batches, use finite element software to simulate and calculate the acoustic transmission loss (STL) in the frequency band concerned, and calculate the average value of the acoustic transmission loss in the frequency band ATL. In this embodiment, the concerned frequency range is the low frequency range of 0-10 kHz, which is divided into 10 frequency bands, and the corresponding ATL value can be obtained by averaging the STL values of the 10 frequency bands on the STL curve.

第三步:将每个透声层数码材料模型的二进制矩阵转换为二进制向量v,与计算得出的ATL值一起,构建称为训练样本,所有训练样本作为训练数据集训练BP神经网络,从而得出训练后的模型,作为二进制向量v与ATL值之间的预测函数,即ATL=f(v)。The third step: convert the binary matrix of each sound-transmitting layer digital material model into a binary vector v, together with the calculated ATL value, construct a training sample, and use all training samples as a training data set to train the BP neural network, thereby The trained model is obtained as a prediction function between the binary vector v and the ATL value, ie ATL=f(v).

第四步:从训练数据集里面筛选出ATL值较低的一部分二进制向量v,将对应的训练样本作为初始种群,通过设置合适的交叉率、变异率和收敛准则,以ATL=f(v)为适应度函数进行遗传算法优化,得出最低ATL预测值以及对应的二进制向量v。在本实施例中,遗传算法直接采用MATLAB中的optimization工具包实现,求解器选择ga-Genetic Algorithm,遗传算子中交叉概率和变异概率分别设置为0.8和0.01。Step 4: Screen out a part of the binary vector v with a lower ATL value from the training data set, and use the corresponding training sample as the initial population. By setting the appropriate crossover rate, mutation rate and convergence criterion, ATL=f(v) The genetic algorithm is optimized for the fitness function, and the lowest ATL prediction value and the corresponding binary vector v are obtained. In this embodiment, the genetic algorithm is implemented directly using the optimization toolkit in MATLAB, the solver is ga-Genetic Algorithm, and the crossover probability and mutation probability in the genetic operator are set to 0.8 and 0.01, respectively.

第五步:将遗传算法得到的最低ATL预测值的二进制向量v,转化成数码材料构型,从而得到最优构型,用有限元方法计算最优构型的声学传递损失STL曲线,并计算平均值ATL。然后判断最优构型的STL曲线以及ATL值是否满足优化预期,如果是,则优化结束,否则,要对训练数据集进行扩充,增加一些根据人的经验生成的构型。相较于完全随机的构型,经验构型的ATL值范围更广,包括了更多特征,又有利于机器学习得到更加准确的预测模型。然后重新进行第二步的有限元模拟以及第三步到第五步的优化过程。Step 5: Convert the binary vector v of the lowest ATL prediction value obtained by the genetic algorithm into a digital material configuration to obtain the optimal configuration, and use the finite element method to calculate the acoustic transmission loss STL curve of the optimal configuration, and calculate Mean ATL. Then judge whether the STL curve and ATL value of the optimal configuration meet the optimization expectations. If yes, the optimization ends. Otherwise, the training data set needs to be expanded to add some configurations generated based on human experience. Compared with the completely random configuration, the empirical configuration has a wider range of ATL values, includes more features, and is conducive to machine learning to obtain a more accurate prediction model. Then repeat the finite element simulation of the second step and the optimization process of the third to fifth steps.

本实施例中采用的经验构型包括前述图5所示的随机梯度夹杂结构、楔形梯度结构和随机楔形结构,其中随机梯度夹杂结构中材料A或材料B出现的概率延厚度方向梯度变化,梯度可以是线性、多项式、对数、指数等形式的;而楔形梯度结构中可做不同梯度的楔型结构,包括直线型、抛物线型、指数型等,曲率变化的参数可以随机产生;另外随机楔形构型的材料中楔形块的外轮廓随机产生,内部在厚度方向会有概率产生空缺。The empirical configuration adopted in this embodiment includes the stochastic gradient inclusion structure, wedge-shaped gradient structure and random wedge-shape structure shown in Fig. It can be in the form of linear, polynomial, logarithmic, exponential, etc.; in the wedge-shaped gradient structure, wedge-shaped structures with different gradients can be used, including linear, parabolic, exponential, etc., and the parameters of curvature changes can be randomly generated; in addition, random wedge In the configuration material, the outer contour of the wedge block is randomly generated, and the inner part has a probability to generate vacancies in the thickness direction.

第七步:得到满足预期的最优构型后,对最优构型进行材料界面的边缘平滑,形成最终构型。Step 7: After obtaining the optimal configuration that meets expectations, smooth the edges of the material interface on the optimal configuration to form the final configuration.

如图7所示,左图展示了本实施例中上述透声层优化方法优化得到的透声层最终构型,该透声层为位于阻抗不匹配的第一介质层和第二介质层之间的双层结构,其中与第二介质层接触的一侧为阻抗跳跃层,与第一介质层接触的一侧为夹杂层,其中阻抗跳跃层完全由第一介质材料A组成;而夹杂层以第二介质材料B为主材,主材内部具有第一介质材料A分布区域,但夹杂层中的第一介质材料A分布区域由第二介质材料B完全围合包裹。As shown in Figure 7, the left figure shows the final configuration of the sound-transmitting layer optimized by the above-mentioned sound-transmitting layer optimization method in this embodiment. The sound-transmitting layer is located between the first dielectric layer and the second dielectric layer with mismatched impedances. A double-layer structure in between, wherein the side in contact with the second dielectric layer is an impedance jumping layer, and the side in contact with the first dielectric layer is an intercalation layer, wherein the impedance jumping layer is completely composed of the first dielectric material A; and the intercalation layer The second dielectric material B is used as the main material, and there is a distribution area of the first dielectric material A inside the main material, but the distribution area of the first dielectric material A in the interlayer is completely surrounded by the second dielectric material B.

由此可见,本发明的优化方法优化得到的构型,与人为经验总结得出的构型并不一致,反而在与第二介质层接触的一侧形成了一层阻抗跳跃层,这违反了人为经验中需要保证阻抗渐变的理念。相较于传统阻抗匹配构型,这种新颖的阻抗跳越的复杂构型具有较好地性能,尤其是在低频区间(0-10kHz),该构型相比于传统阻抗渐变的匹配构型,其平均声传递损失进一步降低了31%。It can be seen that the configuration optimized by the optimization method of the present invention is not consistent with the configuration obtained by artificial experience, and instead forms a layer of impedance jump layer on the side in contact with the second dielectric layer, which violates the artificial The concept of impedance gradients needs to be guaranteed empirically. Compared with the traditional impedance matching configuration, this novel complex configuration of impedance jumping has better performance, especially in the low frequency range (0-10kHz). Compared with the traditional matching configuration of impedance gradient , and its average acoustic transmission loss was further reduced by 31%.

另外,本实施例中,上述边缘平滑后的最终构型中夹杂层内的材料B分布区域近似成椭圆形。基于该最终构型,还可以对其进行结构特征分析,提取出主要结构特征,过滤次要结构特征,然后针对主要结构特征进行简化,并通过有限元验证简化后的最优构型,从而得到最终的构型。例如,参见图7中右图所示,椭圆形的材料B分布区域可以简化为圆形,阻抗跳跃层可平滑形成一层光滑层体,通过对各种参数进行参数化,包括透声层总厚度D、宽度W、阻抗跳跃层厚度D2、材料B区域的半径d、材料B区域底部与透声层底部的间距h,以及W、D2、d、h相对于D的比例wr、tr、dr、hr,即可对这些参数继续进行伸入探究,判断这些参数影响最终性能的机理,以及该构型是否存在进一步优化的可能性。In addition, in this embodiment, the distribution area of the material B in the interlayer in the final configuration after the above-mentioned edge smoothing is approximately elliptical. Based on the final configuration, structural feature analysis can be carried out to extract the main structural features, filter the secondary structural features, and then simplify the main structural features, and verify the simplified optimal configuration through finite elements, so as to obtain final configuration. For example, as shown in the right figure in Fig. 7, the elliptical material B distribution area can be simplified into a circle, and the impedance jump layer can be smoothly formed into a smooth layer body. By parameterizing various parameters, including the overall sound transmission layer Thickness D, width W, impedance jump layer thickness D 2 , radius d of material B area, distance h between the bottom of material B area and the bottom of the sound-transmitting layer, and the ratio of W, D 2 , d, h to D w r , t r , d r , hr , you can continue to explore these parameters, judge the mechanism of these parameters affecting the final performance, and whether there is a possibility of further optimization of this configuration.

以上所述的实施例只是本发明的一种较佳的方案,然其并非用以限制本发明。有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变化和变型。因此凡采取等同替换或等效变换的方式所获得的技术方案,均落在本发明的保护范围内。The above-mentioned embodiment is only a preferred solution of the present invention, but it is not intended to limit the present invention. Various changes and modifications can be made by those skilled in the relevant technical fields without departing from the spirit and scope of the present invention. Therefore, all technical solutions obtained by means of equivalent replacement or equivalent transformation fall within the protection scope of the present invention.

Claims (10)

1. The method for optimizing the sound transmission layer for reducing the acoustic transmission loss of the dissimilar material interface is characterized by comprising the following steps of:
s1, discretizing a longitudinal section of a periodic unit into a binary matrix by taking a single periodic unit in an acoustic transmission layer as an optimization target, and respectively representing a first dielectric material and a second dielectric material by different binary element values of the matrix, wherein the first dielectric material and the material of the first dielectric layer have the same impedance, and the second dielectric material and the material of the second dielectric layer have the same impedance; generating a binary matrix through randomization, generating a random sound-transmitting layer number code material model in batches, wherein each number in the model corresponds to a material unit, and the specific material of each unit is determined by the element value in the corresponding binary matrix;
s2, obtaining an acoustic transmission loss curve and an acoustic transmission loss average value in a target frequency band range by utilizing finite element simulation for each sound transmission layer number code material model generated in the S1;
s3, converting a binary matrix of each sound transmission layer number code material model into a binary vector, taking the binary vector as input of a training sample, and taking the corresponding acoustic transmission loss average value calculated in the S2 as a label of the training sample to construct a training data set; performing machine learning training by using a training data set to obtain a prediction function for predicting an acoustic transmission loss average value according to the binary vector;
s4, taking training samples in the training data set as an initial population, taking the prediction function as an fitness function, and optimizing through a genetic algorithm to obtain an optimal binary vector with the minimum acoustic transmission loss average value;
s5, converting the optimal binary vector into a sound transmission layer number material model to obtain an optimal configuration, obtaining an acoustic transmission loss curve and an acoustic transmission loss average value in a target frequency band range through finite element simulation, judging whether the optimal configuration meets expected optimization requirements according to finite element simulation results, if so, completing optimization, if not, expanding training samples in a training data set, then performing machine learning training again to generate a new prediction function, and performing genetic algorithm optimization again until the optimal configuration meets the optimization requirements;
s6, smoothing the boundary between the two dielectric materials aiming at the optimal configuration meeting the optimization requirement to form the final configuration of the sound-transmitting layer.
2. The method for optimizing an acoustic transmission layer for reducing interface acoustic transmission loss of a dissimilar material according to claim 1, wherein the second dielectric layer is water, and the second dielectric material is hydrogel.
3. The method for optimizing a sound-transmitting layer for reducing acoustic transmission loss of a dissimilar material interface according to claim 1, wherein in S5, the method for expanding training samples in a training data set is as follows: an empirical configuration determined from an artificial summary experience is added to the original randomly generated training dataset.
4. The method for optimizing a sound-transmitting layer for reducing acoustic transmission loss of a dissimilar material interface according to claim 3, wherein the empirical configurations comprise a random gradient inclusion structure, a wedge-shaped gradient structure and a random wedge-shaped structure, each empirical configuration is a sound-transmitting layer code material model and also comprises a binary matrix formed by discretizing a longitudinal section of a periodic unit;
in the sound transmission layer number code material model of the random gradient inclusion structure, the material type of each code is randomly generated, but the probability of the second medium material in the code is gradually decreased line by line along the direction from the side of the second medium layer to the side of the first medium layer, and the probability of the first medium material in the code is gradually increased line by line gradient;
in the sound-transmitting layer number material model of the wedge-shaped gradient structure, the whole first dielectric material is in continuous sheet and has regular wedge-shaped distribution of outer outline, and the wedge-shaped width gradually increases along the direction from the side of the second dielectric layer to the side of the first dielectric layer;
in the sound transmission layer number code material model of the random wedge-shaped structure, the whole first dielectric material is in wedge-shaped distribution with non-connected pieces and irregular outer contours, the second dielectric material randomly appears in the wedge-shaped interior, and the number of the numbers of the first dielectric material appearing in a single line number is gradually increased line by line along the direction from the side of the second dielectric layer to the side of the first dielectric layer.
5. The method for optimizing an acoustically transparent layer for reducing interfacial acoustic transmission loss of a dissimilar material according to claim 4, wherein said probability growth curve of gradient increment in said random gradient inclusion structure is in the form of a linear, polynomial, logarithmic or exponential.
6. The method for optimizing a sound-transmitting layer for reducing acoustic transmission loss of a dissimilar material interface according to claim 4, wherein contour lines on two sides of the first dielectric material distributed in a wedge-shaped gradient structure are straight lines, polynomial curves, logarithmic curves or exponential curves.
7. The method for optimizing a sound transmission layer for reducing acoustic transmission loss of a dissimilar material interface according to claim 1, wherein a machine learning model adopted in the machine learning training is a counter-propagating artificial neural network, a support vector machine, a convolutional neural network or a linear regression.
8. The method for optimizing an acoustic transmission layer for reducing acoustic transmission loss at a dissimilar material interface of claim 7, wherein said initial population is a portion of training samples in a training dataset having an average value of acoustic transmission loss below a threshold value.
9. The method for optimizing an acoustic transmission layer for reducing acoustic transmission loss of a dissimilar material interface according to claim 1, wherein the target frequency band range is 0-10 khz.
10. The sound-transmitting layer optimized by the sound-transmitting layer optimizing method according to any one of claims 1 to 9, wherein the sound-transmitting layer has a double-layer structure between a first dielectric layer and a second dielectric layer with unmatched impedance, wherein the side contacting the second dielectric layer is an impedance jump layer, and the side contacting the first dielectric layer is a sandwich layer, wherein the impedance jump layer is entirely made of a first dielectric material; the second dielectric material is used as a main material of the sandwich layer, a first dielectric material distribution area is arranged in the main material, and the first dielectric material distribution area in the sandwich layer is completely surrounded and wrapped by the second dielectric material.
CN202211543461.4A 2022-12-02 2022-12-02 A sound-permeable layer and its optimization method for reducing the acoustic transmission loss at the interface of dissimilar materials Active CN115762685B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211543461.4A CN115762685B (en) 2022-12-02 2022-12-02 A sound-permeable layer and its optimization method for reducing the acoustic transmission loss at the interface of dissimilar materials

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211543461.4A CN115762685B (en) 2022-12-02 2022-12-02 A sound-permeable layer and its optimization method for reducing the acoustic transmission loss at the interface of dissimilar materials

Publications (2)

Publication Number Publication Date
CN115762685A CN115762685A (en) 2023-03-07
CN115762685B true CN115762685B (en) 2023-07-28

Family

ID=85343065

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211543461.4A Active CN115762685B (en) 2022-12-02 2022-12-02 A sound-permeable layer and its optimization method for reducing the acoustic transmission loss at the interface of dissimilar materials

Country Status (1)

Country Link
CN (1) CN115762685B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103915090A (en) * 2012-12-31 2014-07-09 中国科学院声学研究所 Broadband noise reduction porous-material acoustic liner and equipment
CN105184018A (en) * 2015-10-13 2015-12-23 同济汽车设计研究院有限公司 Method for calculating mid-frequency transmission loss of subsystems and optimizing acoustic packages
CN109165764A (en) * 2018-06-26 2019-01-08 昆明理工大学 A kind of line loss calculation method of genetic algorithm optimization BP neural network
CN113361383A (en) * 2021-06-01 2021-09-07 北京工业大学 Porcelain insulator damage neural network identification method based on genetic algorithm characteristic parameter optimization

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2490150C1 (en) * 2011-12-16 2013-08-20 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тольяттинский государственный университет" Modified laminar acoustic structure of vehicle body upholstery
CN110097185B (en) * 2019-03-29 2021-03-23 北京大学 An optimization model method and application based on generative adversarial network
US20210158211A1 (en) * 2019-11-22 2021-05-27 Google Llc Linear time algorithms for privacy preserving convex optimization
CN112040382B (en) * 2020-08-10 2021-07-30 上海船舶电子设备研究所(中国船舶重工集团公司第七二六研究所) High-bandwidth underwater acoustic transducer based on acoustic impedance gradient matching layer
CN114692813A (en) * 2020-12-29 2022-07-01 清远市富盈电子有限公司 Calculation method for impedance value of PCB (printed Circuit Board)
CN114390427A (en) * 2021-12-29 2022-04-22 瑞声光电科技(常州)有限公司 A sound field optimization method, apparatus, device and readable storage medium
CN114444390A (en) * 2022-01-25 2022-05-06 东南大学 Design method of vacuum electronic device slow wave structure based on deep learning and machine learning

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103915090A (en) * 2012-12-31 2014-07-09 中国科学院声学研究所 Broadband noise reduction porous-material acoustic liner and equipment
CN105184018A (en) * 2015-10-13 2015-12-23 同济汽车设计研究院有限公司 Method for calculating mid-frequency transmission loss of subsystems and optimizing acoustic packages
CN109165764A (en) * 2018-06-26 2019-01-08 昆明理工大学 A kind of line loss calculation method of genetic algorithm optimization BP neural network
CN113361383A (en) * 2021-06-01 2021-09-07 北京工业大学 Porcelain insulator damage neural network identification method based on genetic algorithm characteristic parameter optimization

Also Published As

Publication number Publication date
CN115762685A (en) 2023-03-07

Similar Documents

Publication Publication Date Title
CN108763718B (en) Method for quickly predicting characteristic quantity of flow field when flow-around object and working condition change
CN110473592A (en) The multi-angle of view mankind for having supervision based on figure convolutional network cooperate with lethal gene prediction technique
CN111985523A (en) A 2-exponential power deep neural network quantization method based on knowledge distillation training
CN114722690B (en) Acoustic super-surface sound field rapid prediction method based on variable reliability neural network
CN114186470B (en) A design optimization system and design method of broadband transparent absorber based on multi-objective genetic algorithm
CN117610180B (en) Board shell reinforcing rib generation type design method
CN115762685B (en) A sound-permeable layer and its optimization method for reducing the acoustic transmission loss at the interface of dissimilar materials
Zeng et al. Short-term load forecasting of smart grid systems by combination of general regression neural network and least squares-support vector machine algorithm optimized by harmony search algorithm method
Gu et al. Lowering the sound transmission loss of impedance-matching structures: Data-driven optimization assisted with a priori knowledge
CN114881204B (en) Road surface flatness prediction method based on road element splitting and GA-BP neural network model
CN110779477A (en) Acoustic method for identifying shape of object in real time
CN115270325A (en) Method for reconstructing initial geometric morphology defect of cylindrical shell based on BP neural network
CN117807823B (en) Complex surface sensor layout method for digital twin modeling
CN117454537A (en) Multi-objective optimization design method for acoustic metamaterial stiffened plate structure
Zhu et al. Selection of criteria for multi-criteria decision making of reservoir flood control operation
Carlini et al. A generalized fast marching method for dislocation dynamics
CN114168750A (en) A Knowledge Graph Representation Learning Method Based on Extended Poincaré Sphere
CN108920787B (en) An Analysis Method of Fuzzy Uncertainty of Structure Based on Adaptive Collocation
CN113608264A (en) Combined prestack elastic inversion parameter and depth network target inversion technology
CN104504289A (en) Method for reducing calculating time of intelligent optimization algorithm
CN118798058B (en) Parameter inversion method and terminal equipment of mass concrete mechanism model
CN118627323B (en) Numerical simulation method of wall function based on improved DP algorithm
Yu et al. A novel version of hierarchical genetic algorithm and its application for hyperparameters optimization in CNN models for structural delamination identification
CN115146543B (en) A reliability optimization design method and system for lightweight structural products
CN115712467B (en) A method, device and storage medium for offloading edge computing PLD tasks

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