The Effect of Perforation on the Dynamics of a Flexible Panel

Introduction of holes into plate-like structures is commonly found as one of the practical noise control measures to reduce sound radiation. However, perforation also reduces the panel stiffness and hence increases its vibration. The discussion on this effect is lacking and hence this paper discusses the dynamics of a perforated panel from the results obtained from Finite Element (FE) model. Different hole geometries and arrangement are simulated to investigate their effect on the plate mobility. In general, it is found that increasing the perforation ratio increases the plate mobility. For a fixed perforation ratio, the mobility increases at high frequency (above 1 kHz) for a smaller hole density in the plate. The plate with holes concentrated at the middle shows the largest increase of vibration around the plate centre compared to those uniformly distributed or away from themiddle and concentrated at the plate edges. This is because as the hole separation becomes smaller, the reduction of the global stiffness around the mid area of the plate becomes greater. This also corresponds to the finding here that the mobility is greater at the vicinity of the hole. Different conditions of the plate edges are found to give consistent trend of the effect of perforation.


Introduction
The vibration of engineering structures, particularly those consisting of thin plate-like members, can be a significant source of noise in many situations.It is common in noise control technique to reduce the sound radiation of such structures directly by constructing them from perforates.This technique is known to be capable of reducing considerable noise radiation and has found many practical applications, including safety guard enclosures over flywheels and producing collection hoppers.
However, the effect of perforation on the plate dynamics is rarely discussed, at least in terms of change in the vibration level due to perforation.The recent models [1,2] to calculate the sound radiation from a perforated panel also ignore this effect.The investigation of the effect of perforation on dynamic properties of plates began in the early 1960s in order to determine an accurate stress analysis of perforated panels used to support the tubes in a heat exchanger [3,4].To obtain equivalence with the actual properties corresponding to a solid plate, the effective material properties were defined, namely, effective Young's modulus, effective Poisson's ratio, and effective density.
Soler and Hill [5] later proposed an analytical formula to calculate the bending stiffness of a perforated panel as a function of hole geometry.The formula, however, is only valid for small perforation ratio.Forskitt et al. [6] used the Finite Element Method (FEM) to obtain the dynamic properties, namely, Young's modulus and Poisson's ratio of a perforation with circular holes.The density is calculated based on the fraction of a solid plate.These properties are then used to calculate the natural frequencies of the perforated panel.
According to Burgemeister and Hansen [7], the equation proposed in [6] does not provide correct resonant frequencies when the effective material properties are used.The FEM was then implemented to model the modal response of range of plates with varying perforation geometries, where, from here, the resonant frequencies of the perforated plates were obtained.The resonant frequencies were then compared with that from the FEM model for the solid plate.The results of these comparisons were then fitted to a cubic expression as a general equation for other perforation geometries which is independent of the mode order of the panel.Patil et al. [8] also used a curve fitting technique to develop an effective resonant frequency formula from simulation data for five perforated plates using FEM.The formula was a function of the so-called mass remnant ratio (MR), that is, the ratio of surface area of a perforated plate to that of the solid plate for the same dimension.The formula was then expressed as a seventh order polynomial function of MR.However, the formula was applicable only for the first three modes of vibration.
Recently, Mali and Singru [9] proposed a technique using concentrated negative mass in the Galerkin method to model the natural frequency of a perforated plate.It is shown that the analytical results are in good agreement with those from a FEM model, but only for small size of holes.
As previous works were concerned with only the dynamic properties of the perforation, this paper investigates the change of vibration level due to the introduction of holes into a solid panel.As no established analytical model is available to calculate the mobility of the perforation, the numerical model is therefore implemented using Finite Element Method (FEM).The numerical calculation of mobility of the perforation is then compared against that of the solid plate in one-third octave band frequency.

Theoretical Model.
Mobility is defined as the ratio between the vibration velocity V and the excitation force  given as a function of frequency, () = V()/().For a finite rectangular plate with dimensions  × , the mobility at an arbitrary point (, ) subjected to a point force  at (  ,   ) is given by [10] where Φ  is the th mass-normalized mode shape of the structure,   is the th natural frequency, and  is the damping loss factor.A case often considered is that of a rectangular plate with simply supported boundary conditions as this system provides a simple analytical solution.For a simply supported rectangular plate, the mode shape can be represented by sine functions.This and the natural frequency for mode (, ) are, respectively, where  is the total mass of the plate,  is the mass per unit area, and  =

Validation with Theory.
The model of the perforated panel is developed here using PATRAN and NASTRAN software programs.In order to ensure the FE model used is valid in terms of the element size and setting of the boundary condition, the result is validated against the theoretical mobility as in (1).
A FE model of a solid panel having dimensions of 0.3 × 0.2 m and thickness 1 mm is developed with simply supported condition along all the edges.The panels are made of aluminium having Young's modulus of 70 GPa, density of 2700 kg m −3 , and Poisson's ratio of 0.3.The 2D shell elements are used to model plate structures where a shell structure carries loads in all directions and therefore undergoes bending and twisting, as well as in-plane deformation [11].
The mesh element type used in the FE model is the basic two-dimensional triangular element (Tria shape and Paver).The element size is set to be much smaller than half of the structural wavelength corresponding to the maximum frequency of analysis to give accurate results at high frequencies.The global edge length of each element is set to be 0.003 m which is sufficient to provide accurate result up to 5 kHz.The plate was set to be an ideal simply supported boundary condition, where the displacement is restrained at the plate edges, but the rotation is set to be free.The location of the excitation force and the response are the same for all simulated plates for consistency of plate mobility analysis.
An excitation force of 1 N is applied slightly away from the centroid of the panel at coordinates (0.18, 0.12) where the (0,0) is at the bottom left corner of the plate as illustrated in Figure 1.This is to generate optimum numbers of vibration modes particularly at low frequencies.The damping loss factor is set to be 0.03 which is a typically small value for most engineering structures.The vibration level of the panel is analyzed in the frequency range of 10 Hz to 5 kHz.From the FE calculation, by extracting the vibration velocity from the response point and as a unit force was used, this result gives the mobility.The result is postprocessed using MATLAB.For the theory, the calculation was also performed using MATLAB for 400 modes ( = 20,  = 20) up to 5 kHz with the same location of excitation force and response as well as the value of the damping factor.The comparison of the results is shown in Figure 2  In order to ensure the right meshing is used in this study, the effect of the size of the mesh is also investigated by simulating the mobility of a perforated panel with different global edge length (GEL) of the mesh.Figure 3 shows two examples of the FE model with different meshing size.Figure 4 plots the mobility for four different GELs.It can be observed that the results for perforated plates with GEL = 3 mm and 5 mm have the same results which collapse together for the corresponding frequency range.For larger GEL, that is, GEL = 10 mm and GEL = 20 mm, the mobility shifts to higher frequency.Thus, for the rest of the simulation, GEL = 5 mm is used in this study.

FE Model for Perforated Panels.
The simulated samples of perforated panels are determined by varying the perforation ratio with fixed diameter or number of holes and also by varying the diameter and number of holes with fixed perforation ratio.For comprehensive study, this also includes the arrangement of holes in the plate, which can be a triangular pattern (usually referred to as diagonal pattern) or a square array as well as the hole distributions in the plate, that is, those concentrated on the middle or along the plate edges.The choice of these parameters is summarized in Tables 1, 2, 3, 4, and 5.For convenience of referring of

Advances in Acoustics and Vibration
Frequency (Hz)  the plate, each perforated plate is named according to its parameter, where the first part of the numbers indicates the number of holes, the second part is the diameter, and the third part is the perforation ratio.The plate parameters are chosen to represent any possibilities of perforation configuration in practice.Figure 5 shows examples of the Finite Element (FE) model of perforated panels from PATRAN software with the mesh condition as described in Section 2.2.1.

Validation with Experiment.
A modal test is carried out to validate the FE result.The experimental setup for measuring the mobility is shown in Figure 6.To realise the free-free boundary conditions, the plate was hung vertically using light-rope.At the same location where the Kistler instrumented hammer stroked, a DYTRAN "teardrop" accelerometer was attached to the plate to measure Frequency (Hz) the plate acceleration.The signal was recorded and processed using Data Physics signal analyzer.Figure 7 shows the results from the FE models and the experiment where satisfactory agreement is achieved.Small discrepancies occur which might be due to the nonuniformity of the panel in the experiment which is different from the ideal plate defined in the FE model.

Simulation Results and Discussion
The numerical results from the FE model are presented in one-third octave bands to observe the "average" level of mobility in a frequency band rather than to deal with details of discrepancies in the narrow band.Also for convenience of analysis, it is of interest to observe the amount of mobility level which increases due to perforation.This can be expressed with "the effect of perforation" in dB unit by where   is the mobility of the perforated plate and   is for the solid plate.The results are presented for both mobility Y and effect of perforation X positioned side-by-side.

Effect of Hole Diameter and Number of Holes at Fixed
Perforation Ratio.In Figure 8, it is shown that, by comparing the solid and the perforated plates, introduction of perforation increases the mobility from low to high frequency.

Advances in Acoustics and Vibration
Frequency (Hz)  The level of increment becomes greater as the perforation ratio increases.However, among the perforated plates with fixed perforation ratio, varying the hole diameter as well as the number of holes does not have significant change on the mobility except at high frequencies (above 1 kHz).This can also be seen in the effect of perforation  results.The variation of mobility for different hole diameters (with fixed perforation ratio) can be clearly observed as the perforation ratio is increased.The largest hole diameter gives the highest mobility at high frequencies.This indicates that, at high frequencies, the mobility of the perforations depends on the hole density in the panel where large hole density (many holes, but smaller) is preferred to reduce vibration.Larger hole density has also been shown to produce the greatest noise reduction effect [1,2].
It should be noted here that the high peaks shown in the effect of perforation at low frequencies below 1 kHz, for example, at the fundamental frequency of the plate, that is, 250 Hz, are the effects due to the natural frequency shift because of perforation [7], although the level at the fundamental frequency as seen in Figure 8(e) is also increased by roughly 10 dB compared to the level at the fundamental frequency of the solid plate.ratio, by maintaining the hole diameter and by increasing the number of holes, is plotted in Figure 9 where, as already discussed, the mobility increases as the perforation ratio increases.Figure 9(b) shows that, for 20% perforation ratio, the mobility can increase to almost 10 dB.

Effect of Hole Arrays at Fixed Perforation Ratio.
For different perforation arrays, that is, square and triangular, at fixed perforation ratio in Figure 10, both perforation arrays also show less significant effect on the plate mobility with each other.The mobility again increases as the perforation ratio is increased (cf.Figures 10(b) and 10(d)).This again confirms that it is the number of holes per unit area (hole density) which affects the plate mobility, not the hole arrays.

Effect of Holes Arrangement at Fixed Perforation Ratio.
For different arrangement of holes, that is, uniform across the plate surface, at the edges of the plate and at the center of plate, the result gives different mobility level for fixed perforation.Figure 11 shows that the panel with centralized arrangement of holes shows the highest increment of mobility level, followed by the uniformly distributed holes across the to investigate the effect of the boundary conditions on the mobility of the perforated panels.For the former, in the FE model, the displacement (vertical, z-direction; see again Figure 1) and the rotation at the plate edges are allowed and, for the latter, both the displacement and rotation are constrained.These represent two extreme conditions: the plate with completely free edges and that with completely restrained edges.
The results are shown in Figures 12, 13, 14, and 15 where all the results have the same trend as those from the plate with simply supported boundary condition; that is, the mobility increases when the perforation ratio increases, when the plate has smaller hole density, and when the holes are concentrated at the center of the plate.

Effect of Response Location.
The previous results are calculated for the response point (and excitation location) at a fixed location and therefore the effect of where the holes become closer to the response location is implicit in the results.Here, the effect of the response location at the vicinity of hole is discussed.In Figure 16, it can be seen that the change of mobility is obvious for small perforation ratio where the

Conclusion
The mobility of perforated plates with varying holes geometries and perforation patterns has been investigated by using Finite Element Method.It is found that the mobility increases as the perforation ratio increases.For a fixed perforation ratio, however, the small hole density (small number, but large diameter holes) gives the most effect to increase vibration at high frequency.The effect can also be seen for holes concentrated in the center of the panel, where vibration can be significantly increased compared to those distributed uniformly or arranged around the plate edges at a fixed perforation ratio.The perforation arrays (either triangular or square) were shown to have insignificant effect on the plate mobility for the same perforation ratio.The same phenomenon is also found for different boundary conditions, that is, free-free and clamped-clamped edges, where the results show a similar trend of the effect of perforation with that from the simply supported edges.The forcing and X Y Max 3.52 + 007 @Nd Min 1.77 + 006 @Nd

Figure 1 :
Figure 1: Schematic diagram of the panel and the excitation point.
(a)  with good agreement particularly below 1 kHz.Above 1 kHz, small discrepancies

Figure 2 :Figure 3 :
Figure 2: Comparison of mobility of a solid panel from FE model and from analytical calculation: (a) narrow band and (b) one-third octave band.

Figure 4 :
Figure 4: Mobility of a solid plate with different global edge length of the mesh.

Figure 7 :
Figure 7: Comparison of the mobility for solid and perforated panels from FE model and experiment: (a) solid and (b)-(c) perforated panels at different perforation ratio.(d) Transfer mobility for fixed perforation ratio of 20%.

Figure 8 :
Figure 8: Mobility ((a), (c), and (e)) and effect of perforation ((b), (d), (f)) of the perforated panels with different hole diameters and number of holes at fixed perforation ratio.

Figure 9 :
Figure 9: Mobility ((a), (c)) and effect of perforation ((b), (d)) of the perforated panels with fixed hole diameter and fixed number of holes.
Fixed Hole Diameter or Fixed Number of Holes.The effect of increasing the perforation

Figure 10 :
Figure 10: Mobility ((a), (c)) and effect of perforation ((b), (d)) of the perforated panels with different hole arrays and at fixed perforation ratio.

Figure 12 :
Figure 12: Mobility ((a), (c)) and effect of perforation ((b), (d)) of the perforated panels for free-free edges for (a)-(b) fixed perforation ratio and (c)-(d) different arrays of holes at fixed perforation ratio.

Figure 17 :
Figure 17: Stress contour of the perforated panels with different hole densities: (a) small hole density, (b) large hole density, and (c) hole concentrated at the center.
3/12(1 − V 2 ) is the plate bending stiffness for  is Young's modulus,  is the plate thickness, and  is Poisson's ratio.The next section discusses the FE model used to calculate the mobility of the perforated plate.

Table 1 :
Parameters of the perforated panel with fixed hole diameter and different number of holes and perforation ratio.

Table 2 :
Parameters of the perforated panel with different hole diameter and number of holes and fixed perforation ratios.

Table 3 :
Parameters of the perforated panel with fixed number of holes and different hole diameter and perforation ratio.

Table 4 :
Parameters of the perforated panel with different hole distributions.

Table 5 :
Parameters of the perforated panel with different hole arrays.