J. Ocean Eng. Technol. Search

CLOSE


J. Ocean Eng. Technol. > Volume 34(6); 2020 > Article
Park and Park: Energetics of In-plane Motions in Coupled Plate Structures

Abstract

Energy flow analysis (EFA) has been used to predict the frequency-averaged vibrational responses of built-up structures at high frequencies. In this study, the frequency-averaged exact energetics of the in-plane motions of the plate were derived for the first time by solving coupled partial differential equations. To verify the EFA for the in-plane waves of the plate, numerical analyses were performed on various coupled plate structures. The prediction results of the EFA for coupled plate structures were shown to be accurate approximations of the frequency-averaged exact energetics, which were obtained from classical displacement solutions. The accuracy of the results predicted via the EFA increased with an increase in the modal density, regardless of various structural parameters. Therefore, EFA is an effective technique for predicting the frequency-averaged vibrational responses of built-up structures in the high frequency range.

Nomenclature

<e>: Energy density (J/m2)

cg: Group velocity (m/s)

ω: Angular frequency (rad/s)

k: Wave number (rad/m)

q: Intensity (J/m/s)

1. Introduction

The demand for ecofriendly and lightweight systems has been gradually increasing in the transport machinery industries, e.g., shipbuilding, automobiles, aviation, and railroads, with the growing interest in environmental issues caused by global warming and the development of industrial technology. These industrial changes have increased the proportion of vibration noise in the high frequency range, exacerbating the vibration-noise issue of transportation machinery systems. Therefore, developing technology for predicting the vibration noise in the corresponding frequency range is becoming important for ensuring competitiveness in the transport machinery industries. Typically, deterministic approaches based on displacement methods, such as the conventional finite element method and boundary element method, have been utilized for vibration noise prediction method in the low frequency range. However, because these methods lack efficiency and reliability, statistical approaches such as statistical energy analysis (SEA) and energy flow analysis (EFA) are generally applied in the high frequency range. Numerical methods such as the finite element technique can be applied in EFA because it has an energy governing equation in the form of a differential equation, in contrast to SEA. Thus, EFA has recently emerged as a new alternative to high-frequency range vibration noise analysis, because the modeling efficiency is high and a detailed design review of the design parameters can be performed. In the early days after EFA was proposed by Belov et al. (1977), an energy flow model for out-of-plane motion was developed, which was mainly based on various simple structural element theories, e.g., those of the membrane, Euler beam, and Kirchhoff plate. Recently, however, Park and Hong (2006a), Park and Hong (2006b), and Park and Hong (2008) developed an energy flow model of the Timoshenko beam and Mindlin plate that reflects the shear-deformation effect and rotatory-inertia effect to increase the accuracy of the previously developed EFA model in the high frequency range, thereby expanding the analysis area. When beams or plates are coupled at arbitrary angles in actual structures, such as ships and offshore structures, the out-of-plane motion and in-plane motion are coupled, and the modal density of the in-plane motion increases with the frequency, increasing the importance of the model for in-plane motion in EFA. Therefore, Park et al. (2001) developed an energy flow model for in-plane motion for the EFA analysis of a coupled plate structure in a high frequency range. However, a study on the effectiveness of the energy flow model according to the frequency is needed, because no systematic comparison with the exact solution for the in-plane motion of the plate was performed in the study. In the present study, an exact solution of the in-plane motion for the plate coupled at an arbitrary angle was newly derived, and the effectiveness of the energy flow model for the in-plane motion of the plate was examined through comparison with the analysis results of the energy flow model.

2. Energy Flow Model for In-plane Motion of Plate

The equation of motion for the in-plane wave of a damped plate without a load can be expressed as follows (Dym and Shames, 2013):
(1)
2ux2+(1ν)22uy2+(1+ν)22vxy=(1ν2)ρEc2ut2
(2)
2vy2+(1ν)22vx2+(1+ν)22uxy=(1ν2)ρEc2vt2
where u and v represent the in-plane displacements of the plate in the x and y directions, respectively, Ec = E(1 + ) represents the complex elastic modulus, η represents the structural damping loss factor, ν represents the Poisson’s ratio, and ρ represents the material density of the plate. The energy-flow model for the in-plane motion of the foregoing equation of motion can be derived by separating it into an in-plane longitudinal wave and an in-plane shear wave, as follows (Park et al., 2001):
(3)
cg,l2ηω(2x2+2x2)el+ηωel=Πin,l
(4)
cg,s2ηω(2x2+2x2)es+ηωes=Πin,s
where cg,l and cg,s represent the group velocities of the in-plane longitudinal wave and in-plane shear wave of a thin plate, respectively; 〈el〉 and 〈es〉 represent the energy densities of the in-plane longitudinal wave and in-plane shear wave averaged in space and time respectively; and Πin,l and Πin,s represent the vibrational input power for the in-plane longitudinal wave and in-plane shear wave, respectively.

3. Energy Flow Model for Out-of-plane Motion of Plate

The equation of motion for the out-of-plane motion of a damped Kirchhoff plate without a load can be expressed as follows (Dym and Shames, 2013):
(5)
Dc(2wx4+24wx2y2+4wy4)+ρh2wt2=0
where Dc = D(1+) represents the complex flexural rigidity. The energy flow model for the flexural wave of the Kirchhoff plate derived from the foregoing equation of motion can be expressed as follows (Bouthier and Bernhard, 1995):
(6)
cg,f2ηω(2x2+2x2)ef+ηωef=Πin,f
where cg,f represents the group velocity of the out-of-plane flexural wave of a thin plate, 〈ef〉 represents the energy density of the out-of-plane flexural wave averaged in space and time and Πin,f represents the vibrational input power of the out-of-plane flexural wave.

4. Analysis of Vibration Energy of Plate Coupled at Arbitrary Angle

The plate model coupled at an arbitrary angle was considered for an analytical study of the energy flow model for the in-plane wave of the plate.

4.1 Exact Solution of Coupled Plate Structures

In the exact solution of a coupled plate for which all the outer boundaries are simply supported, the displacement can be calculated in the form of a Levy series that satisfies the y-direction boundary condition equations, i.e., Eqs. (7)(10). The displacement in the x- and y-directions in the ith region can be expressed using Eqs. (11)(13).
(7)
u1(0,y)=u3(Lx,y)=0
(8)
v1(0,y)x=v3(Lx,y)x=0
(9)
w1(0,y)=w3(Lx,y)=0
(10)
2w1(0,y)x2=2w3(Lx,y)x2=0
(11)
ui(xi,y,t)=[m=1Um,i(xi)sinkmy]ejωt
(12)
vi(xi,y,t)=[m=1Vm,i(xi)coskmy]ejωt
(13)
wi(xi,y,t)=[m=1Wm,i(xi)sinkmy]ejωt
Here, km = /Ly and Um,i(xi), Vm,i(xi), and Wm,i(xi) represent the wave solutions in the xi direction that satisfy the equation of motion. The following equation can be obtained by substituting the out-of-plane displacement expressed by Eq. (13) into Eq.(5), i.e., the equation of motion (Park et al., 2001).
(14)
d4Wm,i(xi)dxi42km2d2Wm,i(xi)dxi2+d4Wm,i(xi)dy4+(km4kf4)Wm,i(xi)=0
Here, kf,i=(ρihiω2)/Dc,i4) represents the flexural wavenumber of the plate in the ith plate section, and ρi and hi represent the density and thickness of the plate, respectively. Because Eq. (14) is a fourth-order homogeneous ordinary differential equation, the wave solution of the displacement w in the xi direction can be obtained as follows:
(15)
Wm,i(xi,t)=[Af1+exp(αm,ixi)+Af1exp(αm,ixi)+[Bf1+exp(βm,ixi)+Bf1exp(βm,ixi)]×ejωt
where αm,i2=km2+kf,i2 and βm,i2=km2+kf,i2. Although such an exact solution of the flexural wave, which is an out-of-plane displacement, was used to validate the energy flow model for the flexural wave of a plate, the exact solution of the in-plane displacement and the energy flow solution were not directly compared in previous studies (Park et al., 2001). Park and Hong (2008) derived an exact out-of-plane motion solution from the equation of motion implemented as the simultaneous differential equations for the Mindlin plate. The matrix equation can be obtained as follows by substituting Eqs. (11)(12) into the simultaneous differential Eqs. (1)(2) to obtain the wave solution in the xi direction of the in-plane motion of a plate in a similar manner (Park, 2018).
(16)
xi{Mm,i}={Nm,i}{Mm,i}
Here, the matrices {Mm,i} and [Nm,i] are given by Eqs. (17) and (18), respectively.
(17)
{Mm,i}=U{Um,iUm,ixiVm,iVm,ixi}T
(18)
[Nm,i(ω)]=[0100(1-νi)2km+(1-νi2)Ec,iρiω200-(1+νi)2km00010(1+νi)2km(1+νi)2km-(1+νi2)2ρiω20]
The general solution of the above differential-equation system can be obtained as follows:
(19)
{Mm,i}=n=14Cm,i,nHm,i,neλm,i,nx
where Cm,i,n is the complex coefficient of the nth-order mode, and H→m,i,n and λm,i,n represent the nth eigenvector and eigenvalue of matrices [Nm,i(ω)], respectively. The number of unknowns for each plate boundary of the plate structure in Fig. 1 is 4 flexural waves and 4 in-plane waves, requiring a total of 24 boundary conditions. The boundary conditions at x1 = 0 are expressed as follows.
(20)
u1(0,y)=0
(21)
v1(0,y)x1=0
(22)
w1(0,y)=0
(23)
2w1(0,y)x12=0
The boundary conditions at the excitation-point locations of x1 = x0 and x2 =− Lx1 +x0 are expressed as follows.
(24)
u1(x0,y)=u2(x0Lx1,y)
(25)
v1(x0,y)=v2(x0Lx1,y)
(26)
w1(x0,y)=w2(x0Lx1,y)
(27)
w1(0,y)x1=w2(x0Lx1,y)x1
(28)
Nxx,1(x0,y)=Nxx,2(x0Lx1,y)
(29)
Nxy,1(x0,y)=Nxy,2(x0Lx1,y)
(30)
Vxz,1(x0,y)+(yy0)=Vxz,2(x0Lx1,y)
(31)
Mxx,1(x0,y)=Mxx,2(x0Lx1,y)
Here,
KSOE-2020-063ueq1.jpg
The boundary conditions at x2 = x3 = 0, where the two plates are coupled at an arbitrary angle, are expressed as follows.
(32)
u2(0,y)=w3(0,y)sinθ+u3(0,y)cosθ
(33)
v2(0,y)=v3(0,y)
(34)
w2(0,y)=w3(0,y)cosθ+u3(0,y)sinθ
(35)
Nxx,2(0,y)=Vxz,3(0,y)sinθ+Nxx,3(0,y)cosθ
(36)
Nxy,2(0,y)=Nxy,3(0,y)
(37)
Vxx,2(0,y)=Vxz,3(0,y)cosθ+Nxx,3(0,y)sinθ
(38)
w2(0,y)x2=w3(0,y)x3
(39)
M2(0,y)=M3(0,y)
Lastly, the four boundary conditions at x3 = Lx2 are expressed as follows (similar to the boundary conditions at x1 = 0).
(40)
u3(Lx2,y)=0
(41)
v3(Lx2,y)x3=0
(42)
w3(Lx2,y)=0
(43)
2w3(Lx2,y)x32=0
If the displacement solution is obtained using the foregoing boundary conditions, the energy density and vibration intensity of the flexural wave and the in-plane wave of the plate can be obtained using the following equations:
(44)
efi=Di4{|2wixi2|2+|2wiy2|2+2νi2wixi2(2wiy2)*+2(1νi)2wixiy(2wixiy)*}+ρihi4|wit|2
(45)
qfxi=12Re{Qxz,i(wit)2+Mx,i(2wixit)*+Mxy,i(2wiyit)*}
(46)
qfyi=12Re{Qyz,i(wit)2+My,i(2wiyit)*+Myx,i(2wixit)*}
Here,
(47)
Mx,i=-Dc,i{2wixi2+νi2wiyi2},My,i=-Dc,i{2wiyi2+νi2wixi2},Mxy,i=Myx,i=-Dc,i(1-νi)2wixiy,Qxz,i=-Dc,i{3wixi3+νi3wixiy2},Qyz,i=-Dc,i{3wiyi3+νi3wixi2y}eIi=14Re[Gc,ihi|uiy+vixi|2+Kc,i(uixi+νiviy)(uixi)*+Kc,i(viy+νiuixi)(viy)*+ρihi(|uit|2+|vit|2)]
(48)
qIxi=12Re[Gc,ihi(uiy+vixi)(vit)*+Kc,i(uixi+νiviy)(uit)*]
(49)
qIyi=12Re[Gc,ihi(uiy+vixi)(vit)*+Kc,i(uiy+νivixi)(uit)*]
where Gc,i = Ec,i/{2(1+νi)} and Kc,i=Ec,ihi/(1νi2)

4.2 Energy Flow Solution of Coupled-plate Structures

The energy flow solution of the coupled plate can be obtained as in previous studies (Park et al., 2001). The Levy series solution that satisfies the y-direction boundary condition of the energy governing equation in Eqs. (3), (4), and (6) was expressed, and the energy flow solution was calculated using the boundary condition of the wave solution in the xi direction.
(55)
eni(xi,y)=m=0Em,n,i(xi)cos(kmy)(n=f,l,s)
Here, the wave solution of the energy density solution for each n-type wave is given by Eq. (45).
(56)
Em,n,i(xi)=Am,n,i+eαm,n,ixi+Am,n,ieαm,n,ixi
Here, Am,n,i+ represents the mth-order mode coefficient of the wave solution traveling in the n-type ± direction in the i-plate region, and Am,n,i+.

4.3 Energy Density Comparison of Coupled-plate Structures

It was assumed that the size of the plate structure in Fig. 1 was Lx1 = Lx2 = Ly =1 m and that the coupled angle of the plate was 30°. The material of the plate was aluminum (E = 7.1 ×1010 N/m2 and ρ = 2700 kg/m3), and the structural damping loss factor of the plate was set as η = 0.01 to express the spatial damping to a certain degree. All the plate thicknesses were h = 0.001 m. The vertical excitation force was applied to the center of the first plate position (Lx1, Ly) = (0.5 m, 0.5 m), and the size was 10 N. The number of modes for obtaining an approximate solution was set as 200 to improve the accuracy of the analysis. The analysis results for an exact solution at an excitation frequency of 10 kHz are shown in Figs. 2,5. The energy density of the flexural wave was the highest at the excitation point because the plate was excited with a vertical load, and it tended to decrease as the distance from the excitation point increased. Additionally, the plate had the highest energy density at the boundary where the plates were coupled at an arbitrary angle, because the in-plane wave was generated at the position where the two plates were coupled at an arbitrary angle (in contrast to the flexural wave), and the energy density decreased owing to the increased structural damping toward the inside of each plate. The propagation tendency of the vibration energy can be observed in the excitation part of the vibrational intensity distribution, as shown in Figs. 4 and 5. The modal density of the plate vibration mode with respect to the frequency is presented in Table 1. Although the theoretical modal density of a flexural wave is constant with respect to the frequency, the in-plane wave tended to increase in proportion to the frequency.
The effectiveness of the EFA, which is a statistical approach, increased in the high frequency range with a high modal density. The energy density distributions of in-plane waves for exact solutions at excitation frequencies of 30 and 100 kHz are shown in Figs. 6 and 7, respectively. The modal density increased with the frequency, increasing the mode superposition.
The comparison results for the energy density of the exact solution and the energy flow solution at different frequencies, with Ly = = 0.5 m, are shown in Figs. 8,10.
Compared with the flexural wave, the modal density was relatively low for the in-plane wave at the same frequency, as indicated by the foregoing analysis results. Therefore, although the exact solution and energy flow solution for the flexural wave in the previous studies had similar energy density distributions owing to mode superposition even in a relatively low frequency range, the spatial distributions of the energy density of the exact solution and the energy flow solution became similar at a significantly high frequency. However, the in-plane energy flow solution successfully showed the average distribution of the exact solution even in a relatively low frequency range, indicating that the EFA is extremely effective for predicting the overall vibration energy level in the middle–high frequency broadband. The size of the plate was doubled (Lx1 = Lx2 = Ly = 2 m) to examine the effects on the plate size, plate coupled angle, and structural damping loss factor, as shown in Figs. 11,13, which present comparisons of the energy density between the exact solution and the EFA solution of the model with a plate coupled angle of 45° and a structural damping loss factor η = 0.001 for the plate material. Overall, the results were similar to those for the first example, even when the plate size, coupled angle, and structural damping loss factor were different.

5. Conclusion

In this study, the EFA, which is attracting attention as an effective vibration noise prediction method in the middle–high frequency ranges, was used to examine in-plane motions in coupled-plate structures. In previous studies, are view of the energetic characteristics depending on the frequency was performed by comparing the energy flow solution and the exact solution for the out-of-plane motion of the plate. However, a study on the reliability of the energy flow solution for in-plane motion was required, because no comparative review of in-plane motion had been conducted. In the present study, an EFA was performed to predict the average energy density level of the exact solution successfully for both out-of-plane and in-plane motions even in relatively low frequency ranges with a high out-of-plane modal density and low in-plane modal density. The EFA method is considered to be extremely effective for vibration analysis in the middle–high frequency ranges for engineers who are interested in the spatial distribution of the overall vibration energy level inside the structure rather than the spatial mode characteristics. In the future, the in-plane wave in the coupled model of the Mindlin plate, which has improved reliability in the high frequency range, must be verified.

Acknowledgements

This research was supported by Changwon National University 2019–2020

Fig. 1
Simply supported plates coupled at an arbitrary angle
KSOE-2020-063f1.jpg
Fig. 2
Flexural energy density level (dB) distribution of exact solutions with f = 10 kHz and η = 0.01
KSOE-2020-063f2.jpg
Fig. 3
In-plane energy density level (dB) distribution of exact solutions with f = 10 kHz and η = 0.01
KSOE-2020-063f3.jpg
Fig. 4
Flexural intensity distribution of exact solutions with f = 10 kHz and η = 0.01
KSOE-2020-063f4.jpg
Fig. 5
In-plane intensity distribution of exact solutions with f = 10 kHz and η = 0.01
KSOE-2020-063f5.jpg
Fig. 6
In-plane energy density level (dB) distribution of exact solutions with f = 30 kHz and η = 0.01
KSOE-2020-063f6.jpg
Fig. 7
In-plane energy density level (dB) distribution of exact solutions with f = 100 kHz and η = 0.01
KSOE-2020-063f7.jpg
Fig. 8
Comparison between the exact solution and the energy flow solution with f = 10 kHz and η = 0.01: ✸, exact solution; □, energy flow solution
KSOE-2020-063f8.jpg
Fig. 9
Comparison between the exact and energy flow solutions with f = 30 kHz and η = 0.01: ✸, exact solution; □, energy flow solution
KSOE-2020-063f9.jpg
Fig. 10
Comparison between the exact and energy flow solutions with f = 100 kHz and η = 0.01: ✸, exact solution; □, energy flow solution
KSOE-2020-063f10.jpg
Fig. 11
Comparison between the exact and energy flow solutions with f = 10 kHz and η = 0.001: ✸, exact solution; □, energy flow solution
KSOE-2020-063f11.jpg
Fig. 12
Comparison between the exact and energy flow solutions with f = 30 kHz and η = 0.001: ✸, exact solution; □, energy flow solution
KSOE-2020-063f12.jpg
Fig. 13
Comparison between exact and energy flow solutions when f = 100 kHz, η = 0.001: ✸: exact solution; □ : energy flow solution
KSOE-2020-063f13.jpg
Table 1
Theoretical modal density of plate according to frequency change
Frequency (Hz) Modal density (modes/Hz)

Flexural Longitudinal Shear
10,000 0.3188 0.0064 0.0021
30,000 0.3188 0.0192 0.0064
100,000 0.3188 0.0636 0.0213

References

Belov, VD., Rybak, SA., & Tartakovskii, BD. (1977). Propagation of Vibrational Energy in Absorbing Structures. Journal of Soviet Physics-Acoustics, 23, 115-119.

Bouthier, OM., & Bernhard, RJ. (1995). Simple Models of the Energetics of Transversely Vibrating Plates. Journal of Sound and Vibration, 182(1), 149-164. https://doi.org/10.1006/jsvi.1995.0187
crossref
Dym, CL., & Shames, IH. (2013). Solid Mechanics: A Variational Approach, Augmented Edition. New Yorkm USA. Springer-Verlag: https://doi.org/10.1007/978-1-4614-6034-3

Park, CH. (2018). A Study on High-Frequency Vibrational Energetics of Kirchhoff Plate Structures Coupled at Arbitary Angles (Master’s Thesis). Changwon National University; Changwon, Korea.

Park, D.-H., Hong, S.-Y., Kil, HG., & Jeon, J.-J. (2001). Power Flow Models and Analysis of In-plane Waves in Finite Coupled Thin Plates. Journal of Sound and Vibration, 244(4), 651-668. https://doi.org/10.1006/jsvi.2000.3517
crossref
Park, Y.-H., & Hong, S.-Y. (2006a). Vibrational Energy Flow Analysis of Corrected Flexural Waves in Timoshenko Beam – Part I: Theory of an Energetic Model. Shock and Vibration, 13, 137-165. https://doi.org/10.1155/2006/308715
crossref pdf
Park, Y.-H., & Hong, S.-Y. (2006b). Vibrational Energy Flow Analysis of Corrected Flexural Waves in Timoshenko Beam – Part II: Application to Coupled Timoshenko Beams. Shock and Vibration, 13, 167-196. https://doi.org/10.1155/2006/562762
crossref pdf
Park, Y.-H., & Hong, S.-Y. (2008). Vibrational Power Flow Models for Transversely Vibrating Finite Mindlin Plate. Journal of Sound and Vibration, 317(3–5), 800-840. https://doi.org/10.1016/j.jsv.2008.03.049
crossref
TOOLS
Share :
Facebook Twitter Linked In Google+ Line it
METRICS Graph View
  • 0 Crossref
  •   Scopus
  • 2,775 View
  • 56 Download


ABOUT
BROWSE ARTICLES
ARTICLE CATEGORY

Browse all articles >

PUBLICATION ETHICS
FOR CONTRIBUTORS
Editorial Office
President Office BD Rm. 1302, 13 Jungang-daero 180beon-gil, Dong-gu, Busan 48821, Republic of Korea
Tel: +82-51-759-0656    Fax: +82-51-759-0656    E-mail: ksoehj@ksoe.or.kr                

Copyright © 2024 by The Korean Society of Ocean Engineers.

Developed in M2PI

Close layer
prev next