An Optimization Method for Maximizing the Low Frequency Sound Insulation of Plate Structures

A combined approach based on finite element method, boundary element method, and genetic algorithm (FEM-BEM-GA) is proposed for optimizing the low frequency sound (LFS) insulation performance of plate structures. This approach can identify the optimal structural parameters (especially concerning the effects of arbitrary boundary conditions) so as to maximize the structural overall LFS insulation. The basic ideas of this approach are as follows: (1) the sound transmission loss (TL) analysis of a plate with arbitrary boundary conditions is conducted by the coupled FEM-BEM method; (2) the single-number rating method (such as low frequency sound transmission class) is used to assess the plate’s overall LFS insulation; and (3) the genetic algorithm (GA) is employed for searching the optimal solutions of the multiple-parameter optimization problem. The proposed approach is subsequently illustrated by numerical studies. The results show the effectiveness of consideration of the effects of boundary condition in the plate’s LFS insulation optimization and demonstrate the feasibility and effectiveness of this approach as a structure design tool.


Introduction
There are numerous types of noise around buildings that seriously influence people's living quality [1].Among these noises, low frequency noise is a major component of many occupational and community noises which is emitted by numerous sources in the society [2][3][4].The lack of attenuation of low frequency noise by walls, windows, and other structures and its pervasive ambient levels make low frequency noise a factor of critical importance to people's living quality and health.An widely accepted frequency range for low frequency sound (LFS) is from 20Hz to 250Hz which covers 12 one-third octave bands [4].The noise within this frequency range can cause not only audible annoyance but also nonaudible annoyance (such as vibration/rattle annoyance) [5]; and it even has negative impacts on people's health and well-being [6].
The study of plate structures can serve as a first step in understanding and manipulating the dynamic and acoustic behavior of more complicated constructions [7].Examples include the walls and windows of buildings, factory machinery casings, parts of vehicle shells, and the hulls and bulkheads of ships.These structures however often result in poor sound insulation, especially at low frequency [8].The structural parameters, including plate material, size, thickness, and boundary condition, are the key factors that determine the plate's sound insulation performance.Particularly, the effects of arbitrary boundary conditions have recently received great attention.Both analysis and experimental studies [9,10] have demonstrated that boundary condition significantly affects the plate's sound insulation, especially in the low frequency domain.Although the plate's sound insulation (especially the LFS insulation) is highly dependent on the plate's actual boundary condition, according to the authors' best knowledge, no optimization technique is yet available for maximizing the sound insulation of plate structures by considering the effects of arbitrary boundary conditions.This kind of methods however has many potential advantages because it only requires the modification of boundary supports and does not need to make any change to the main body of the structure.
The purpose of this study is to develop an optimization method for maximizing the plate's LFS insulation.This method is supposed to be capable of considering arbitrary plate parameters simultaneously as variables (especially concerning the effects of arbitrary boundary conditions).To this end, a finite element method, boundary element method, and genetic algorithm (FEM-BEM-GA) combined method is proposed, in which the coupled FEM-BEM model can accurately predict the low frequency sound transmission loss of the plate with arbitrary boundary conditions while the genetic algorithm (GA) can efficiently handle the multipleparameter optimization problems.The single-number rating of low frequency sound transmission class is used as the objective function in the optimization.

Sound Transmission Loss (TL) Analysis Model
Consider a rectangular plate of length L x , width L y , and thickness h, with arbitrary elastic boundary supports along the four edges, as shown in Figure 1.The plate is supposed to be mounted on an infinite rigid baffle, which separates the fluid medium (air medium) into source section Π + and receiver section Π − .Arbitrary incident sound waves varying harmonically are incident to the plate in the source section, and the vibration of the plate induced by the incidence waves then radiates sound waves into the receiver section.The vibroacoustic behavior of the plate can be determined by the coupled finite element and boundary element method (FEM-BEM), in which the plate vibration response is given as [11] and the sound pressures on the plate surfaces are given as {} is the global nodal displacement vector of the plate.{ + } and { − } are the radiated sound pressures on the front (in section Π + ) and back (in section Π − ) plate surfaces, respectively.{ 0 } is the vector representing the pressure of the incidence sound wave.{}, {}, and [] are the mass, damping, and stiffness matrices of the plate structure.{} is a square matrix formed by the "collocation" procedure.
{R} is a global transformation matrix converting the nodal displacement vector to the transverse deflection vector.{T} is a global transformation matrix converting the fluid pressure to point forces that act on the nodes of the plate.The damping matrix {} is supposed to be proportional to the stiffness matrix {} and is written as {} = {}, where  = 2/ 0 ,  is the damping factor, and  0 is the fundamental natural frequency of the plate structure.Particularly, in order to consider arbitrary elastic boundary supports, the stiffness matrix {} of the whole plate structure is decomposed into plate and boundary supports [11], which is expressed as {} = {  } + {  }. {  } and {  } are the stiffness matrices for the plate and boundary supports, respectively.The arbitrary elastic boundary supports are modeled as a combination of translational and rotational springs, with   and   being the translation stiffness and rotational stiffness, respectively.The element stiffness matrix of {  } can then be expressed as where {  } is the shape function vector for the plate element and  → n  is the normal unit vector of the element boundary contour Γ  .It is noted that the elastic parameters (  and   ) along the contour can arbitrarily be varied to reproduce simply supported (  = ∞ and   = 0), clamped (  = ∞ and   = ∞), free (  = 0 and   = 0), and guided edges (  = 0 and   = ∞), or, more important, any intermediate situation (i.e., general boundary condition).Moreover, these parameters can spatially vary along each edge to represent arbitrary nonuniform elastic restraint.The sound transmission loss (TL) of the plate can then be defined as where  + and  − are the sound power at the front plate surface (in the source section Π + ) and the radiated sound power at the back plate surface (in the receiver section Π − ).They can be expressed as where Δ is the area of the plate element, { + } = −{ − } = {}, and the superscript asterisk denotes the complex conjugate.
Combining ( 1)-( 6), the sound transmission loss of the plate with arbitrary elastic boundary conditions can be determined.
The accuracy of this TL model has been validated by the experiment [11].The TL measurements were conducted on an aluminum (Al) plate structure.The Al plate is 1mm thick and is installed in a steel frame with actual radiated area of 24cm by 24cm.The plate is with uniform elastic boundary supports along the four edges and the actual boundary parameters are   = 3201 and   = 13.28,estimated by the boundary condition identification (BCI) method [11].  and   are the dimensionless forms of   and   [10,11].Figure 2 shows a comparison of the measured and calculated TL results of the plate structure [11].The predicted results agree well with the experimental data, which demonstrates the effectiveness of the TL prediction model in studying the low frequency sound (LFS) insulation of plate structures with elastic boundary supports.

Single-Number Rating Indicators for Low Frequency Sound Insulation
TL is a powerful descriptor that gives comprehensive information about the sound insulation performance and allows the engineers to efficiently analyze this performance at different frequency ranges.In practice, a simpler characterization, namely, by one single number, is also demanded, especially for the purpose of easy ranking of structures and verification of the compliance with legal requirements [12].
Tokita and Nakanura Thresholds are validated to be good indicators of the potential for low frequency noise annoyance [5].These thresholds present the sound pressure levels at which subjective response to direct exposure to low frequency sound changes as functions of frequency from one characterization to another.Based on the Tokita and Nakanura Thresholds, a single-number rating method called low frequency sound transmission class (LFSTC) is recently developed [13], which is useful for assessing the structural sound insulation in the low frequency range (20Hz∼250Hz).The rating indicator, LFSTC O&N , is defined by this singlenumber rating method, which is developed based on the "oppression and noisy" threshold and is suitable for the case when high level LFS is present.The LFSTC O&N is employed in this paper to assess the plate's overall LFS insulation.The classification procedure of LFSTCO&N is easy to understand and implement, which is quite similar to that given in the international standard (ASTM E413-10 [14]).The reference contour used for calculating LFSTC O&N is given in Table 1 (in 1/3 octave bands).More details about this method can be found in [13].
The well-known C-weighted metric is another good indicator of the potential for low frequency noise annoyance [5].The C-weighted metric is considered to be applicable to most industrial sources of steady low frequency noise [15] and is used by World Health Organization (WHO) as an important indicator for providing guidance on the health effects of low frequency noise [16].The rating indicator, LFSTC C , is also defined in this paper to assess the plate's overall LFS insulation, which is based on the reference contour of Cweighted metric (also given in Table 1).The classification procedure of LFSTC C is the same as that of LFSTC O&N .

Optimization Methodology
The objective is to identify the optimal plate parameters (especially the optimal boundary parameters) so as to maximize the LFS insulation of the plate.This is a multipleparameter optimization problem, which requires the optimization method to be capable of considering multiple plate parameters simultaneously as variables.Even if only the boundary condition is set as the design variable, there can be multiple boundary parameters for a practical plate and each boundary parameter can take on an infinite number of different values.This makes it impossible to apply an exhaustive search to the task of structural optimization.Therefore, global search algorithms that can search intelligently for the optimal solution within the search space are needed.In this study, a genetic algorithm (GA) is employed and combined with the coupled FEM-BEM method for the development of structural optimization strategy.
4.1.Genetic Algorithm.Genetic algorithms (GAs) are heuristic search techniques based on principles present in natural evolution.Compared to other stochastic methods, GAs have a distinct advantage that it is extremely easy to parallelize the algorithm since the calculations of each iteration are independent of one another [17].The general procedure of a GA consists of four bioinspired operators, including initialization, crossover, selection, and mutation.This procedure is inspired by the processes of natural selection and genetics and creates "environments" where the best quality individuals tend to pass on their data to the next generation.The quality of the population can therefore increase over time.In the search process, the fitness function (objective function) is a key concept which is used for the assessment of the fitness of individuals (i.e., the quality of individuals).

FEM-BEM-GA Combined Optimization
Approach.The flowchart of the optimization strategy is given in Figure 3, which can be described as follows: (1) Input the known variables of the plate.
(2) Input the constraints of the design variables.
(3) Run the FEM-BEM-GA combined method to search the optimal values of the design variables.In this combined method, the low frequency sound transmission loss of the plate with arbitrary elastic boundary conditions is calculated by the coupled FEM-BEM method.The GA is applied to search the optimal structural parameters so as to maximize the LFS insulation of the plate.The singlenumber rating indicator (such as LFSTC O&N and LFSTC C ) is employed to assess the plate's overall LFS insulation and is used as the fitness function in GA.

Illustrative Examples
The proposed method is applied to the LFS insulation optimization of two plate cases.The element number used in the coupled FEM-BEM model is 100 (10×10).The subscript n of   and   represents the edge number (see Figure 1).In the following illustrative examples, the single-number ratings, LFSTC O&N and LFSTC C , are used as the fitness functions.The optimization algorithm is run 5 times independently in each case to ensure finding the best results; and in each run the initial population in GA is set to be 100 and the optimization algorithm terminates if no improvement in the best solution is observed for 10 consecutive generations.

Case 1 (Boundary Condition Optimization
).The known parameters, design constraints, optimization target, and final optimization results are shown in Table 2.The plate's boundary condition along Edges 1 and 4 is set to be clamped; and the boundary parameters of the other two edges (Edges 2 and 3) are taken as design variables.For comparison, the single-number rating indicators, LFSTC O&N and LFSTC C , of the same plate with two classic boundary conditions, SSSS (simply supported on all edges) and CCCC (clamped on all edges), are also calculated (see Table 2).It is proved that when the boundary parameter's value (  or   ) is larger than 1 × 10 10 , increasing the value of this boundary parameter does not affect the plate's TL results anymore [9,13].For simplicity, unless stated otherwise, the infinite large value of boundary parameter is represented by 1 × 10 10 .For example, the boundary condition of clamped supported can be represented as   = 1 × 10 10 and   = 1 × 10 10 .
From Table 2, for LFSTC O&N , it can be seen that the optimal boundary parameters are found which can maximize LFSTC O&N (LFSTC O&N =45 dB).As a result of this optimal boundary condition, increases of 10 dB and 8 dB in LFSTC O&N are obtained in comparison with those of SSSS and CCCC boundary conditions, respectively.For LFSTC C , it also can be seen that the optimal boundary parameters are found which can maximize LFSTC C (LFSTC C =42 dB).As a result of this optimal boundary condition, increases of 9 dB and 8 dB in LFSTC C are obtained in comparison with those of SSSS and CCCC boundary conditions, respectively.

Case 2 (Boundary Condition and Thickness Optimization).
The known parameters, design constraints, optimization target, and the final optimization results are shown in Table 3.The plate's boundary parameters along the four edges are taken as design variables, as well as the plate's thickness.It is assumed that the plate is with uniform boundary supports along the four edges and that the thickness is an integer multiple of 1 mm.Such situations are common in practical applications.For comparison, the LFSTC O&N and LFSTC C results of a reference plate (the same plate but with 15 mm thickness, i.e., the upper limit of thickness constraint) with SSSS and CCCC boundary conditions are also calculated (see Table 3).
For LFSTC O&N , it can be seen that the optimal boundary parameters and thickness are found which can maximize LFSTC O&N (LFSTC O&N =31 dB).As a result of these optimal structural parameters, increases of 6 dB and 3 dB in LFSTC O&N are obtained in comparison with those of the reference plate with SSSS and CCCC boundary conditions, respectively.For LFSTC C , it also can be seen that the optimal boundary parameters and thickness are found which can maximize LFSTC C (LFSTC C =32 dB).As a result of this optimal boundary condition, increases of 7 dB and 4 dB in LFSTC C are obtained in comparison with those of the reference plate with SSSS and CCCC boundary conditions, respectively.

Computation Time of Optimization.
A normal personal desktop computer manufactured in year 2015, with 8 GB of memory and Intel (R) Core (TM) i7-4790 processor (8 cores, 3.6 GHz), was used for running the simulation cases.During the whole calculation process, the average usage of

Results
Optimal Boundary condition Note: for comparison, the LFSTC O&N results of the same plate with two classic boundary conditions, SSSS (simply supported on all edges) and CCCC (clamped on all edges), are 35 dB and 37 dB, respectively. (d)

Results
Optimal Boundary condition Note: for comparison, the LFSTC C results of the same plate with two classic boundary conditions, SSSS (simply supported on all edges) and CCCC (clamped on all edges), are 33 dB and 34 dB, respectively. 6

Shock and Vibration
Meet termination criteria?(terminates if no improvement in the best solution is observed for 10 consecutive generations) CPU and memory was less than 45% and 35%, respectively.The average time of calculating TL (see ( 1)-( 6)) is about 0.01s per frequency.The average time of calculating the fitness function for each combination of design variables is about 3.21s, including the time of calculating TL for the whole frequency range of interest (10Hz-300Hz) and the time for calculating LFSTC O&N (or LFSTC C ).

Mutation
Even if only one boundary parameter is set as the design variable, it can take on an infinite number of different values; therefore it is impossible to apply an exhaustive search for the optimization.Even if only integer numbers are considered for the design variables, the number of possible combinations is still very large.Taking Case 2 as an example, this number can be estimated as (10 10 − 10 3 ) × 10 10 × (15 − 5) ≈ 1 × 10 21 .It also takes almost infinite time (about 3.21×10 21 s) for performing an exhaustive search and therefore cannot be solved without using intelligent search algorithms.
Table 4 shows the computation time of the two cases (in Section 5) by using the proposed optimization method.It can be seen that the time of running the optimization algorithm (i.e., the whole procedure of Figure 3) is about 20∼ 26 minutes.The proposed method indeed can significantly reduce the searching time and make the optimization much more efficient and applicable.

Performance of Optimization. The examples given in
Section 5 demonstrate the boundary condition has significant effects on the plate's LFS insulation and is worthy of consideration in structural optimization.The optimal results obtained based on both LFSTC O&N and LFSTC C indicators show the following: (1) For a given plate, the LFS insulation performance can be effectively improved even if only modifying the plate's boundary parameters (such as LFSTC O&N and LFSTC C of Case 1).This is useful in practical applications especially when the design constraints do not allow modifying the body and appearance of the given structure (such as the plate's material, size, shape, mass, and surface property).(2) The plate with a proper boundary condition could have even better LFS insulation performance than that of the thicker plate.In such case, it is possible to design a plate with both better LFS insulation and smaller thickness (smaller mass) within the design constraints (such as LFSTC O&N of Case 2).
It can also be noted from the examples that different fitness functions (i.e., different rating indicators, LFSTC O&N , and LFSTC C ) lead to different optimal results.The optimal values of plate parameters are highly dependent on the selected LFS insulation indicators.Users can choose freely the LFS insulation rating indicator (LFSTC O&N or LFSTC C ) according to their actual needs.

Applications of Optimization.
With the development of elastic material and structural assembly technology, it is proved that elastic plate boundary conditions can be realized in practical applications.For example, the plate structure with elastic mounting can be achieved by placing elastic materials along the border of the plate [18].Lower-density elastic materials, such as Getzner Sylomer and Sylodyn, are of a type suitable for use in an elastic boundary setup [18,19].The

Target
Maximize LFSTC C  parameters of these materials, such as their stiffness, can be obtained from their suppliers [19] or checked and corrected by the boundary condition identification (BCI) method [11].Also, it should be admitted that, at the current stage, it is not easy to exactly set the plate boundary condition to be arbitrary optimal values in practical applications.The experimental technique for controlling the plate's boundary parameters to exactly match the optimal results is still being studied (it is also our focus).Nevertheless, the optimization idea and technique proposed in this study can serve as a fundamental step in understanding and optimizing the LFS insulation of plate structures.

Conclusion
A finite element method, boundary element method, and genetic algorithm (FEM-BEM-GA) combined optimization approach is developed for maximizing the low frequency sound (LFS) insulation of plate structures.In this approach, the low frequency sound transmission loss of the plate with arbitrary elastic boundary conditions is calculated by the coupled FEM-BEM method, while GA is employed for searching the optimal structural parameters so as to maximize the plate's LFS insulation.Meanwhile, the single-number rating of low frequency sound transmission class is applied to assess the plate's overall LFS insulation.
The proposed approach is illustrated by numerical studies.The results demonstrate the effectiveness of considering boundary parameters in the plate's LFS insulation optimization, as well as the effectiveness of this approach as a structural design tool.According to the authors' best knowledge, no optimization method is yet available in the literature for maximizing the LFS insulation of a plate structure by considering the effects of arbitrary boundary conditions; Shock and Vibration therefore the proposed idea as well as the optimization model can be valuable for both academic and design practitioners.

Figure 1 :
Figure 1: Schematic illustration of a baffled rectangular plate structure: (a) a plate mounted on an infinite rigid baffle; (b) elastic boundary supports along the edges.

Figure 2 :
Figure 2: Comparison of the predictions and measured data.

Table 1 :
Reference sound insulation contours for calculation of single-number ratings.

Table 2 :
Known parameters and optimization of Case 1.

Table 3 :
Known parameters and optimization of Case 2.  1 =  2 =  3 =  4 = 11707.61, 1 =  2 =  3 =  4 = 0.37 : for comparison, the LFSTCO&N results of the reference plate (the same plate but with 15 mm thickness) with SSSS and CCCC boundary conditions are 25 dB and 28 dB, respectively. Note 91 Note: for comparison, the LFSTC C results of the reference plate (the same plate but with 15 mm thickness) with SSSS and CCCC boundary conditions are 25 dB and 28 dB, respectively.

Table 4 :
Computation time of optimization.