Vibrational Energy Flow Model for a High Damping Beam with Constant Axial Force

)e energy flow analysis (EFA) method is developed to predict the energy density of a high damping beam with constant axial force in the high-frequency range. )e energy density and intensity of the beam are associated with high structural damping loss factor and axial force and introduced to derive the energy transmission equation. For high damping situation, the energy loss equation is derived by considering the relationship between potential energy and total energy.)en, the energy density governing equation is obtained. Finally, the feasibility of the EFA approach is validated by comparing the EFA results with the modal solutions for various frequencies and structural damping loss factors. )e effects of structural damping loss factor and axial force on the energy density distribution are also discussed in detail.


Introduction
With the development of high-speed aircrafts and transportation vehicles, high-frequency vibration of components is of great concern for both academic and industrial communities in recent years. As one of the typical structural components, beams with high damping treatment for vibrational reduction are extensively used in mechanical and aerospace engineering. Moreover, these beam-type structures sometimes experience axial forces arising from initial stress and temperature variation, which can significantly influence the dynamic characteristics of structures. Consequently, high-frequency dynamic response prediction of a high damping beam with axial force is of great significance for a beam-type structure design.
Both the statistical energy analysis (SEA) [1] and the energy flow analysis (EFA) [2] are popular methods to predict structural time and space averaged energy densities in high-frequency range. SEA only can provide a single averaged energy value for each subsystem of a built-up structure. As an alternative method to SEA, EFA can provide the locally smoothed energy value at every interested point in a built-up structure. As a consequence, EFA is more advantageous over SEA in solving high-frequency dynamic responses due to its ability to obtain detailed local response information.
Wohlever et al. [3] derived the energy density governing equation of hysterically damped rods and beams from traditional displacement solutions. Chen et al. [4] calculated the transient response of beams under high-frequency shock load by energy finite element method and virtual mode synthesis method. Liu et al. [5] developed the vibrational energy flow model for functionally graded beams for highfrequency response prediction. Kong et al. [6] theoretically tested the validity region and criterion of EFA from the formation of reverberant plane wave field. Nokhbatolfoghahai et al. [7] and Navazi et al. [8] verified the vibrational energy flow model of a uniform beam and plate from experimental data. Zhu et al. [9] developed a hybrid method for midfrequency vibrational problems, which combined FEA for the stiff member and EFA for the reverberant field of the flexible member. Sadeghmanesh et al. [10] presented a method for defining a criterion to select a proper structural theory based on the order of shear deformation and rotary inertia for low-to high-frequency vibration analyses. Furthermore, the influences of the fluid loading on the dynamic behaviors of structures were considered such that EFA was applied to model the interaction of the fluid with structures [11][12][13]. ere usually exist the discontinuities, which refer to material and structural configuration variation for the coupled structures in engineering. e energy density is not continuous on coupling boundaries where wave reflection and transmission arise from encountering the discontinuities. In the conventional EFA coupling approach, FEM is used to solve the energy density governing equation and the wave-based method is utilized to obtain the power transferring coefficients [14][15][16]. Park and Hong [17] and Zhihui et al. [18] developed a hybrid power flow analysis method using coupling loss factor of SEA to solve the coupling problems in EFA. Kwon et al. [19] and Zheng et al. [20] developed an energy flow model for high-frequency vibroacoustic analysis of complex structures using EFA method.
For structures with initial stress, Zhang et al. [21] and Di et al. [22] investigated the effects of the axial force due to thermal expansion on the high-frequency vibration of beams and plates. Recently Zhang et al. [23] have derived the energy density governing equation by considering the energy density and intensity associated with axial force, which is neglected in [21,22]. Unfortunately the developed EFA can only be used to analyze the high-frequency vibration of lightly damping beams. In this paper, we extend the EFA model developed by Zhang et al. [23] to high damping beams with axial force. is paper is organized as follows. Firstly, the energy density governing equation is derived from the energy transmission equation and the energy loss equation for a high damping beam with axial force. Particularly, the relationship between energy density and intensity of structures and the dissipated energy equation are derived from the wavenumber obtained without approximation of the real and imaginary term. en the accuracy of the presented formulation is verified by comparing the developed EFA results with modal analysis results of a pinnedpinned beam with high damping loss factor and axial force. e effects of damping loss factor and axial force on the energy density distribution are also discussed. Finally, the present EFA method is applied to beams with different boundary conditions and conclusions on the new EFA model are given.

Energy Transmission Equation.
For a beam with constant axial force, driven by a harmonic point force as depicted in Figure 1, the flexible motion equation can derived as [24,25] where w is the transverse displacement of the vibrating beam, D c � D(1 + jη) is the complex bending stiffness with D representing the real part of D c and η the structural damping loss factor, N c � N(1 + jη) is the complex axial load with N representing the real constant part of N c , ρ is the mass density, S is the cross-section area, F 0 e jωt δ(x − x 0 ) is the transverse harmonic load applied at point x 0 with F 0 denoting its amplitude, and ω is its circular frequency, δ(x − x 0 ) is the Delta function, L is the length of the beam. Substituting the general solution w(x, t) � Ce j(kx+ωt) into equation (1) gives the dispersion relation: en the complex wavenumber, k, can be solved from equation (2) and be presented as where k c1 � k 1 + jk 2 is the complex wavenumber, k 1 and k 2 are the corresponding real and imaginary parts, respectively. In terms of equation (2), the group velocity c g can be expressed as equation (4) by implicit derivative with respect to the real part k 1 of the complex wavenumber k c1 : In a lightly damped system, from equations (3) and (4), the higher-order terms with respect to the structural damping can be neglected and the flexural wavenumber and group velocity can be approximated as [21,23] As a result, the general solution of equation (1) is expressed as w(x, t) � Ae −jk c1 x + Be jk c1 x + Ge − k c2 x + He k c2 x e jωt , (7) where A, B, G, and H are complex coefficients that can be calculated from the boundary conditions. e right first two terms of equation (7) are far-field solutions representing the propagating wave, while the right last two ones of equation (7) are near-field solutions representing the evanescent wave. For high-frequency vibration, the flexible wavelength is quite short. e near-field solutions are neglected beyond one wavelength away from the boundaries and driving points where they are dissipated rapidly.
us, only the far-field solutions are taken into account for EFA and the expression of displacement is rewritten as e energy density in a vibrating beam is the sum of the kinetic and potential energy densities which are associated with bending strain and constant axial force, respectively. e time averaged energy density can be represented as [3,26] 〈e〉 � 1 4 Re ρS zw zt where the bracket operator 〈·〉 denotes the time average over one period and the superscript * denotes the complex conjugate operator. e total intensity is associated with the shear force, the bending moment, and the constant axial force. e time averaged intensity of the vibrating beam is presented as [3,26] 〈I〉 � Substituting equation (8) into equations (9) and (10), and space averaging over a wavelength, the time and space averaged energy density and intensity can be obtained as where 〈·〉 denotes the time and space averaged operator to eliminate the spatially harmonic terms, α 1 and α 2 are the space averaging values of e 2k 2 x and e − 2k 2 x . From equations (11) and (12), the relationship between the time and space averaged energy density and intensity, socalled energy transmission equation, is derived as

Energy Loss Equation.
In the case of hysteric damping, the dissipated power of per volume smoothed by time and space averaging is proportional to the potential energy density smoothed by time and space averaging, for an elastic system harmonically vibrating with frequency ω where e p is the potential energy density. In a lightly damped system, the potential energy density is equal to the kinetic energy, and hence the dissipated power is expressed as In a vibrating beam with constant axial force, the time averaged potential energy density can be expressed as [27] Substituting equation (8) into equation (15) and smoothed by local space averaging, equation (15) can be rewritten as Combining equations (11), (14), and (16), the relationship between time and locally space averaged energy density and the dissipated power of per volume can be represented as

Energy Density Governing Equation.
For the elastic medium, the power balance equation at the steady state is expressed as Mathematical Problems in Engineering where π in refers to the input power at the driving point Substituting equations (13) and (17) into equation (18), the energy density governing equation for the high damping beam with axial force can be derived as Unlike the lightly damping system in [23], the energy transition equation is derived using the exact wavenumber without neglecting higher-order damping terms. e energy transition equation represents the relationship between time and space averaged energy density and intensity without group velocity because the damping effect is not neglected. On account of high structural damping, the potential energy density is not the same as the kinetic energy density of the structure.
us, the energy loss equation is derived according to the relationship between potential energy and total energy.

Solution of Energy Density Governing Equation
As shown in Figure 1, the beam is divided into two regions at x 0 . From equation (19), the general solution of the flexural energy density governing equation is expressed as where the subscript i denotes the i th region of the beam. e constants C i1 , C i2 can be determined by the energy boundary conditions. From equations (13) and (20), the energy intensity of the vibrating beam is represented as Because there are no energy outflow and inflow at both ends of the pinned-pinned beam, the following energy intensity boundary conditions can be determined as At the interface between the regions ① and ②, the energy density is continuous and the energy intensity is subject to conservation of energy. e resulting boundary condition can be written as e input power of the vibrating beam is obtained by the impedance method and expressed as [23,28] where Z is the impedance at the driving point of the vibrating beam with axial force. Substituting equation (20) into equations (24)∼(26), the matrix equation can be obtained as According to equation (29), the flexural energy density can be analytically calculated for the vibrating beam with high structural damping loss factor and constant axial force.

Verification and Discussion
To validate the proposed EFA formulation and investigate the effects of structural damping loss factor and axial force, various energy flow analyses are performed for a pinnedpinned beam at both ends shown in Figure 1. e results from analytical solutions of EFA governing equation are compared with exact modal solutions (Appendix A and Appendix B), which are defined as time averaged energy density obtained from the displacement solutions. For each modal analysis, the first 5000 modes are extracted to ensure the accuracy of the solution. In addition, the present vibrational energy flow method is applied to beams with different boundary conditions.

Model Verification.
e beam is made of aluminum alloy with density ρ � 2700kg/m 3 and elastic modulus E � 70 GPa, and its dimensions are L � 3mS � 2 × 10 −4 m 2 , e structural damping loss factor η ranges from 0.05 to 0.8 and the axial force N from 0 to 20 kN, respectively. e pinned-pinned beam, as shown in Figure 1 is excited by a transverse harmonic force located at x 0 � L/2. Figure 2 shows the energy density distribution of pinned-pinned beams for the case with the axial force N � 20 kN, the amplitude of transverse harmonic force F 0 � 1 N, and the structural damping loss factor η � 0.5 when the analysis frequencies are 2, 4, 6, and 8 kHz. Except for that the developed EFA results differ near both ends of the beam, the developed EFA results is analogous to the modal analysis results because the evanescent wave is neglected in the energy flow models. Moreover, it can be observed that the difference of the EFA solution from the modal solution is smaller at higher frequencies, and hence the accuracy of the EFA model is improved at higher frequencies. e energy density distributions with different amplitudes of transverse harmonic force (F 0 � 0.1, 1, 10 kN), as shown in Figure 3, are calculated when the excitation frequency is f � 6 kHz. It is also found that the EFA solutions are smooth representing the overall trend of modal solutions and in good agreement with modal solutions except the vicinity of boundaries. In addition, the global variation of energy density increases with larger amplitudes of harmonic point force. e exact energy density distribution obtained by modal analysis method oscillates spatially without space averaging. At high-frequency ranges, the variance of classical solutions is uniform resulting from short wavelength. erefore, the developed EFA results, obtained with space averaging, are analogous to the classical solutions as the analysis frequency increases. Figures 4 and 5 show the energy density distributions of pinned-pinned beams with N � 20 kN, F 0 � 1 N, and f � 6 kHz when the structural damping loss factor is η � 0.05 and 0.5, respectively. If the structural damping loss factor is small, the wavenumber and group velocity given in equations (3) and (4) in the developed energy flow model are well approximated by those in equations (5) and (6) in the traditional energy flow model. Hence, in a low damping beam the developed EFA solution is close to the traditional EFA solution. However, in a high damping beam, the developed EFA result represents the trend of exact solution from modal analysis, and the developed EFA result concurs more with the exact solution than the traditional EFA result because higher-order damping terms are no longer neglected. Figure 6 shows the energy density distributions of pinned-pinned beams with different axial forces 0kN ≤ N ≤ 20kN, amplitude of transverse harmonic force F 0 � 1 N, and structural damping loss factor η � 0.5 when the analysis frequency is 6 kHz. e energy density varies more dramatically along x position with increasing axial force. e total energy can be calculated from integrating the energy density over the length of the beam. Figure 7 shows that the total energy is higher at larger axial force. e effects of constant axial force can be explained by the input power and group velocity. Figures 8 and 9 show the input power and group velocity, respectively. It is observed that the input power increases and the group velocity decreases with increasing axial force. Because the total energy is proportional to the input power, the total energy is higher for larger axial force arousing larger input power. e lower group velocity makes energy propagate more slowly from the driving point to the boundary. It means that more energy dissipates near the driving point and less energy propagates to the boundary. As a result, the energy density changes more dramatically with larger axial force along the x position.    Figure 10 shows the energy density distributions of pinned-pinned beams with different structural damping loss factor 0 ≤ η ≤ 0.8, amplitude of transverse harmonic force F 0 � 1 N, and axial force N � 20 kN when the analysis frequency is 6 kHz. e developed EFA results represent the global variation of the modal solutions at each structural damping loss factor. Owing to the exclusion of the evanescent wave, the developed energy flow solutions more closely approach the modal solutions as the structural damping loss factor increases. e energy density changes more suddenly along x position with increasing structural damping loss factor. Figure 11 shows that the total energy is lower at larger structural damping loss factor. Figures 12 and 13 show the input power and group velocity at various structural damping loss factors, respectively. It is found that the input power decreases and the group velocity increases with an increase of structural damping loss factor. Since the total energy is proportional to the input power, the total energy is lower for larger structural damping loss factor leading to less input power. In fact, the larger group velocity contributes to flatter change of the energy density along the x position. But the energy density distribution changes more dramatically with larger structural damping loss factor along the x position. It attributes to that the energy density is mainly affected by the input power while the effect of group velocity is minor.

Application to Beams with Different Boundary
Conditions. e analysis model consists of a beam fixed at both ends as shown in Figure 14 and a beam fixed at one end and free at the other as shown in Figure 15, respectively. Both beams share the common dimensions of L � 3mS � 2 × 10 −4 m 2 I b � 2 × 10 −10 m 4 , material properties of ρ � 2700kg/m 3 and E � 70GPa, and the structural damping loss factor η � 0.5. Since the energy outflow at the fixed point is zero and the energy flow satisfies the equilibrium condition at the driving point, the energy intensity boundary conditions for the fixed-fixed beam is identical to those for the pinned-pinned beam. However, the impedance of fixed-fixed beam at the driving point is obtained by substituting equations (B.5) and (B.6) in Appendix B into equation (A.8) in Appendix A. For fixedfree beam, the subdomain ② is absent so that the coefficients C 21 , C 22 in equation (29) vanish, the following matrix equation is rewritten as Similarly the impedance of fixed-free beam at driving point is determined by substituting the equation (B.8) and (B.9) in Appendix B into equation (A.8) in Appendix A. From equation (30), the energy density for the fixed-fixed beam can be evaluated. Figure 14 shows a fixed-fixed beam with constant axial force. Figures 16 and 17 show the energy density distribution of fixed-fixed beams and fixed-free beams with F 0 � 1 N, N � 20 kN, and η � 0.5 when the analysis frequencies are 2 kHz, 4 kHz, 6 kHz, and 8 kHz, respectively. e developed EFA results coincide well with the overall trend of the exact solutions from modal analysis at each frequency. Similar to previous cases, the developed energy flow model gives more accurate result at high frequencies despite the exclusion of the evanescent wave. As a consequence, it is feasible to predict the vibrational energy density of beams with various boundary conditions in the high-frequency range by the proposed energy flow model.

Conclusion
An energy flow model for a beam with axial force that includes high damping effect is established to evaluate the energy density in the high-frequency range. e wavenumber is expressed as the function of the damping loss factor and the axial force. Particularly, the energy transmission equation and the energy dissipation equation are derived from the wavenumber without approximation of its real and imaginary term. en, the energy density governing    To verify the developed energy flow model, numerical analyses are performed for the simply supported beam with various damping loss factor and axial force. e EFA solutions are greatly consistent with the modal solutions in all cases. Furthermore, it is found that both the high damping factor and large axial force can significantly alter the level of energy density as well as the distribution. Due to the resulting deceasing input power, the increasing damping loss factor makes the energy density decay more dramatically along axial position. However, that is the contrary case for the axial force. e proposed method is expected to be useful for the prediction of the high-frequency vibration of high damping structures with initial stress.
represents the so-called "modal force" and are the so-called "modal mass" and "modal stiffness". e effect of structural damping is modeled by assuming a hysteretic damping model which introduces an additional complex term (1 + jη) to modal stiffness. e natural frequency of the ith mode, ω i , can be determined from the corresponding wavenumber in equation (2), which can be written as e forced vibration response of the beam can be expressed as the linear combination of the normal modes e impedance at the driving point of the finite beam with axial force can be obtained from its definition: According to the general solution in equation (7) and boundary conditions with one fixed end and the other free end, the frequency equation is derived as cos k i L cosh k i L + 1 � 0.
(B.8) e vibration mode of the beam is written as ϕ i (x) � cos k i x − cosh k i x + cos k i L + cosh k i L sin k i L + sinh k i L · sin k i x − sinh k i x .

(B.9)
Data Availability e data used to support the finding of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there no conflicts of interest.