Dynamic Analysis of Rectangular Plate Stiffened by Any Number of Beams with Different Lengths and Orientations

The present work is concerned with dynamic characteristics of beam-stiﬀened rectangular plate by an improved Fourier series method (IFSM), including mobility characteristics, structural intensity, and transient response. The artiﬁcial coupling spring technology is introduced to establish the clamped or elastic connections at the interface between the plate and beams. According to IFSM, the displacement ﬁeld of the plate and the stiﬀening beams are expressed as a combination of the Fourier cosine series and its auxiliary functions. Then, the Rayleigh–Ritz method is applied to solve the unknown Fourier coeﬃcients, which determines the dynamic characteristics of the coupled structure. The Newmark method is adopted to obtain the transient response of the coupled structure, where the Rayleigh damping is taken into consideration. The rapid convergence of the current method is shown, and good agreement between the predicted results and FEM results is also revealed. On this basis, the eﬀects of the factors related to the stiﬀening beam (including the length, orientations, and arrangement spacing of beams) and elastic parameters, as well as damping coeﬃcients on the dynamic characteristics of the stiﬀened plate are investigated.


Introduction
e stiffened plate component can be regarded as a coupling structure that is composed of a plate and several beams. In practical engineering, the connection between the plate and beam does not only involve the classical coupling, but the elastic coupling is also frequently encountered. In addition, due to the complexity of the actual working conditions, the stiffened plates will be often subjected to the elastic boundary restrictions. A good understanding of the dynamics of the stiffened plate with general boundary restrains can provide direct benefits to the structural design of complex systems. However, there are only a few literatures on the dynamic analysis of the rectangular plate with some beams of arbitrary lengths and orientations. e current work is to present the dynamic analysis of the beam-stiffened rectangular plate with general boundary restraints.
In the past decades, many scholars have made a lot of efforts on free vibration characteristics of the stiffened plate. Mukherjee and Mukhopadhyay [1] presented the free vibration of the isoparametric stiffened plate element, which could accommodate irregular boundaries. Holopainen [2] proposed a new finite element model to consider free vibration of eccentrically stiffened plate. In this research, a plate with any number of arbitrarily oriented stiffeners was presented. Based on the Rayleigh-Ritz method, Liew et al. [3] and Lee and Ng [4] both investigated the free vibration of the stiffened plate with arbitrarily oriented stiffeners. e difference was that Liew studied the stiffened plate with various shapes, while Lee only examined the stiffened rectangular plate. A hierarchical finite element in conjunction with local trigonometric interpolation functions was employed by Barrette et al. [5] to investigate the vibration performances of the stiffened plate. Employing the differential quadrature technique, free vibration of plate with eccentric stiffeners was performed by Zeng and Bert [6]. And their theoretical solutions were in good agreement with the experimental results. Using the finite element method, Srivastava et al. [7] explored vibration of the stiffened plate subjected to partial edge loading. However, only the simply supported and clamped cases were considered in their research.
On the basis of the first-order shear deformable theory, Peng et al. [8] used mesh-free Galerkin method to carry out the free vibration and stability analyses of the stiffened plate. On the framework of the assumed-modes method, low-frequency free vibration analysis of the thin rectangular plate with a few light stiffeners was made by Dozio and Ricciardi [9]. Based on the improved Fourier series method, Xu et al. [10] examined the free vibration of the beam-stiffened rectangular plate, where beams of arbitrary lengths and placement angles were taken into consideration. ree different modelling approaches were reported by Bhar et al. [11] to carry out the component-wise vibration analysis of the stiffened plate, and comparative studies between different models were also given. e static, free vibration and buckling analysis of the stiffened plate were carried out by Nguyen-oi et al. [12] using CS-FEM-DSG3. In Ref. [13], the CS-FEM-DSG3 method was also extended to analyze the static and free vibration of stiffened folded plates, where three-node triangular element was adopted. Four hierarchical models were developed by Bhaskar and Pydah [14] to obtain the natural frequencies of the isotropic and orthotropic stiffened plate with simply supported boundary condition. On the basis of the elastic approach, Pydah and Bhaskar [15] established an analytical model to investigate the statics and dynamics of the stiffened plate with simply supported case. Based on the assumed mode method, Cho et al. [16] investigated natural frequencies and modes of various opening stiffened panels under different boundary conditions, where the effect of the parameters of the rib on the frequency was checked. en, he and other co-authors [17,18] carried out the free vibration analysis of the stiffened panel with lumped mass and stiffness attachments, including rectangular plates without opening and with various opening shapes. Utilizing the NURBS based on isogeometric analysis approach, Qin et al. [19] examined the frequency characteristics of curved stiffened plate under in-plane loading. Considering a beam-stiffened plate, Yin et al. [20] examined its vibration transmission by applying the dynamic stiffness method, where the transmission mode was presented. Meanwhile, Damnjanović et al. [21] adopted the same approach to perform the free vibration analysis of the stiffened composite plate based on the highorder shear deformation theory, and the results were in good agreement with the available data and FEM results.
rough the abovementioned review of literatures, the predecessors have conducted in-depth research on the free vibration of the stiffened plate. Nevertheless, many scholars are not satisfied with research in the field of free vibration. Harik and Salamoun [22] applied an analytical strip method to examine bending deflection of the stiffened plate, where the type of load and the number of ribs were used to determine the number of strips. Transient response of the simply supported panel reinforced by stiffeners was predicted by Louca et al. [23] using the finite element technique, where the effect of boundary condition on transient response was considered. Sheikh and Mukhopadhyay [24] applied the spline finite strip method to analyze the linear and nonlinear transient vibration of the stiffened plate with arbitrary shapes, in which stiffeners had arbitrary orientation and eccentricity. Based on finite element analysis, Srivastava et al. [25] presented the dynamic instability of the stiffened plate, where both in-plane partial and concentrated loadings at the edge were examined. Xu et al. [26]. developed the structural intensity technique to analyze energy flow for the stiffened plate applied in the field of marine structures. Later, the same authors [27] investigated the structural intensity of the stiffened plate by means of the finite element method. In this research, they draw a conclusion that the energy flow of the whole structure was closely related to the natural frequency. Employing the harmonic compound strip method, Borković et al. [28] obtained the transient response of the beamstiffened plate, where Bernoulli-Euler beam and Kirchhoff-Love thin plate were both applied. Combining the Mindlin plate theory and Timoshenko beam theory, Cho et al. [29] extended the assumed mode method to predict forced vibration of the stiffened plate and rectangular plate with various boundary restrains. Tian and Jie [30] carried out the research on the dynamic response of the finite ribbed plate under point force moment excitations, and he concluded that the input mobilities were affected by the distance between the excitation point and the beam. en, he [31] performed the vibration response analysis and experiment of the stiffened plates subjected to the clamped support at edge. In 2018, he [32] introduced a double cosine integral transform technique to study the free vibration and forced response of the stiffened plate with free boundary restrains. Meanwhile, the energy flow of stiffened panels was reported by Cho et al. [33], where the structural intensity technique was applied. In order to efficiently solve the dynamic responses of the stiffened plate in wide frequency domain, Jia et al. [34] proposed a Composite B-spline Wavelet Elements Method (CBWEM) and verified the reliability of the method by comparison with experimental results. Applying the Mindlin plate and Timoshenko beam theory, Liu et al. [35] analyzed the dynamic power transmission characteristics of a finite stiffened plate with various classical cases, where the effects number and geometric dimensions on the stop band in the low-frequency domain were shown.
rough the review of the abovementioned literature, it can be found that the existing researches are limited to the dynamic of beam-stiffened plate with classical boundary conditions (like clamped, simply supported, and free). Only the free vibration analysis related to elastic boundary conditions and nonfixed connections between a plate and beams is carried out in Xu's et al. [10] research. However, the investigations on the mobility characteristic, structural intensity, and transient response analysis of the stiffened plate subjected to the elastic boundary are still blank. In an attempt to fill this gap, the current work presents the dynamic analysis of the rectangular beam-reinforced plate with general boundary restraints, including classical and the elastic restraints as well as elastic coupling of beams and plate. IFSM is extended to construct displacement field of the beam-stiffened plate, which has proven to be a very effective method from existing research [10,[36][37][38]. e artificial spring technology is introduced to obtain general boundary conditions and various coupling connections between a plate and beams. Based on Xu's et al. [10] research, the energy expressions of the coupling system is obtained, which can be solved by the Rayleigh-Ritz method. Transient 2 Shock and Vibration response of the stiffened plate with classical or elastic boundary restrains is also calculated by the Newmark method. e consistency of the current data and FEM results is revealed. Also some novel numerical results with respect to the structural intensity, structural mobility, and transient response of the stiffened plate are also reported.

Theoretical Formulations
2.1. Descriptions of the Model. e geometry of a typical beamstiffened plate and the coordinate system of the coupling structure are depicted in Figure 1. For the sake of illustration, only a thin plate reinforced by a beam with any lengths and angles is shown in Figure 1(a), where a, b, and h denote the length, width, and thickness of the plate, respectively. In addition, the primary coordinate system (o-x-y-z) is placed on the middle surface of the plate, where u, v, and w represent the midsurface displacements of the plate in the primary coordinate system, respectively. e width and height of the cross section for a stiffening beam are denoted by b b and h b , respectively, which is in the local coordinate system (o'-x'-y'-z'). Four types of boundary springs are used on each edge of the plate to obtain arbitrary boundary conditions, in which three types (k w , k u , and k v ) of linear springs and one type of rotational spring (K w ) are included. For a thin plate, the clamped boundary condition will be obtained when the spring stiffness value of each side of the plate is infinity. Conversely, the free boundary is easily achieved when the stiffness value is 0. In current work, the stiffening beam and plate are described in separate systems. As shown in Figure 1(a), six sets of coupling springs are used to achieve the connection between the beam and plate, including three linear springs (k pb1 , k pb2 , and k pb3 ) and three types of rotational springs (K pb1 , K pb2 , and K pb3 ). By setting stiffness values of coupling springs, different coupling relationships between the beams and plate can be achieved, including various elastic coupling and clamped coupling. It should be noted that the coupling between two intersecting beams is ignored in this paper. From Figure 1(b), the orientation between the local coordinate system (o'-x'-y'-z') attached to each beam and the primary coordinate system (o-x-y-z) is denoted by the symbol of φ when the length of the beam is expressed by L. And the coordinates of the starting end of the beam in the primary coordinate system are represented by (x 0 , y 0 ).

Energy Functionals.
In this study, both the in-plane displacement components (u, v) and the out-of-plane displacement component (w) of the plate are considered. erefore, the total strain potential energy (V p ) of the plate is composed of the out-of-plane strain potential energy (V out p ) and the in-plane strain potential energy (V in p ). Specifically, the expression of the total potential energy is as follows: in which Similarly, the kinetic energy (T p ) expression of the plate can be written as In equations (2)∼(4), D P � E p h 3 /(12(1 − μ p )) is the bending stiffness of the plate and E p , G p , μ p , ρ p , and ω represent Young's modulus, extensional rigidity, Poisson's ratio, mass density, and circular frequency, respectively [10].
Since the boundary conditions of the plate are realized by the artificial virtual spring technology, the elastic potential energy (V sp ) stored in the boundary springs here should be taken into the potential energy of plate. V sp can be expressed as

Shock and Vibration
In this paper, four degrees of freedom for a single beam are considered, which are two bending displacement components (wy b and wz b) in the y′ and z′ directions and the axial displacement component (u b ) and torsional displace- erefore, the strain potential energy (V bi ) of ith beam can be written as e kinetic energy of ith beam can be expressed as In equations (6)∼(7), Dz bi and Dy bi are, respectively, the bending rigidity in the x'-z' plane and x'-y' plane when torsional is denoted by the symbol of J bi . In addition, Young's modulus, shear modulus, mass density, and cross-sectional area of ith beam are taken place of E bi , G bi , ρ bi , and A bi , respectively.
Since the sti ening beam and plate is connected by a coupling spring, the coupling potential (Vpb Ci) of the coupling spring needs to be considered in the coupling system. According to references [10,39], its speci c expression can be written as It should be noted that six degrees of freedom between the beam and the plate are considered here. As shown in equation (8), three sets of linear springs (k pb1 , k pb2 , and k pb3 ) are used for linear displacement constraints and three sets of rotational springs (K pb1 , K pb2 , and K pb3 ) are used for angular displacement constraints. Besides, w z bi and w y bi are shown as follows

Displacement Function.
In the present research, the admissible displacement functions of a plate and a beam are established by IFSM, which has been successfully applied to solve the vibration of beams [40,41], plates [42,43], and shells [44]. e improved Fourier series consists of main functions and auxiliary functions. For the auxiliary functions, theoretically, they may have in nite forms. Considering programming and computational e ciency, a set of sine functions is used as auxiliary functions. Hence, the displacement functions of the plate can be written as follows:  Shock and Vibration in which cos λ 0 x sin λ n y , . . . , cos λ m x sin λ n y , . . . , cos λ M x sin λ n y , . . . , where λ m � mπ/a λ n � nπ/b, and M, N represent the number of truncations. Fourier coefficients of displacement functions for the plate are expressed by Au mn, a ml , b ln , Bv mn, c ml , d ln , Cw mn, e mp , and f pn (l � 1, 2; p � 1, 2, 3, 4), respectively. Similarly, the displacement field functions of the ith beam can be obtained by IFSM as follows: Shock and Vibration 5 where In the above equations (14a)∼(15d), λ m1 � m 1 π/L i and M1 indicates the truncated number. D mi , E mi , F mi , and G mi represent the unknown coefficient of ith beam's displacement functions.

Steady-State Vibration Solution.
For the calculation of the steady-state response for the beam-stiffened plate under external excitation, it is necessary to consider the work done by the external force. In terms of the coupling structure shown in Figure 1(a), the Lagrange function (L) can be expressed as follows: in which the number of stiffening beams is denoted by N, i � 0 means unreinforced bare plate, and W exc represents the work done by the external excitation force acting on the plate: In the above equations (14a) and (14b), the external load distribution function on the plate is denoted by the symbols of f u , f v , and f w . In the current research, f w is taken into consideration, which expresses the lateral load associated with out-of-plane displacement (w) of the plate. In particular, for external point excitation, f w is represented by the magnitude (F) of the point excitation force and the 2D delta function (δ): Substituting equations (1)∼(8), (10), (13), and (17) into equation (16), linear equations in the matrix form are available by employing the Rayleigh-Ritz method: where the vector H For the external excitation force with any frequency, the unknown series expansion coefficient related to displacement functions of the coupled structure can be directly calculated from equation (19). at is, Substituting equation (22) into equation (10), the displacement response of plate reinforced by beams under some specific excitation will be obtained. And vibration speed can also be obtained by v � jωw. Particularly, the complex Young's modulus (E) is introduced in the current numerical solution process. E � E p (1 + jη) is composed of Young's modulus (E p ) and structural loss factor (η), where is the imaginary unit. 6 Shock and Vibration e structure mobility can quantitatively describe the law of power flow in the structure. erefore, understanding the admittance characteristics of the stiffened plate is of great significance for structural vibration control. e structure mobility can be calculated by the following equation: in which Y nk represents the transfer mobility from point n to point k when V k is utilized to express the velocity at point k.
And the magnitude of the excitation force on point n is presented by F n . e mobilities of the drive point n can be achieved by setting n � k. e vibration energy of the structure can be described by the structural intensity vector, which helps the designers understand the energy distribution of the structure. e definition of structural intensity can be found in references [26,37,45]. In this paper, the structural intensity vector I (x, y) of any point in the stiffened plate can be obtained by vector superposition of the component I x (x, y) along the xaxis and component I y (x, y) along the y-axis.
eir relationship is seen in the following equation: In equation (24), the structural intensity at any point includes both bending vibration components and in-plane vibration components, as shown in equations (25a) and (25b): In the above equations (25a) and (25b), the structural intensity components associated with structural bending vibrations in both directions are in which the symbol of ( * ) represents the conjugate of a complex number. And the transverse shear, bending moment, and torque in the two directions are, respectively, e structural intensity associated with the in-plane vibration components in equations (25a) and (25b) is presented as follows: Substituting equations (10), (22), and (26a)∼(28b) into (25a) and (25b), the structural intensity at any point of the stiffened plate can be obtained.

Transient Vibration Solution.
e mass matrix M and stiffness matrix K of the coupled structure have been presented in equations (20)∼(21) in Section 2.4. Hence, the dynamic equation of the structural system in the time domain can be described as [46] M € x where F p is the external excitation vector, which is the function of time, and x denotes the displacement at time t.
Generally, damping constants have frequency-varying characteristics. erefore, it is difficult to accurately define the damping matrix in numerical simulation. For the sake of simplicity, only Rayleigh damping of the coupling structure is considered in this study, namely, C � α 0 M + β 0 K. And α 0 Shock and Vibration 7 and β 0 are the Rayleigh damping coefficients independent of frequency.
In order to solve the transient response of the system, the Newmark method is adopted, which is an extensive implicit algorithm. Its integration method is as follows: In equations (30)∼(31), when the velocity function and the displacement function are expanded by Taylor series, the expansion is retained to the second derivative and equations (32a)∼(33) can be calculated as follows: Combining equations (29), (32a), and (32b), the following matrix equation (33) can be obtained: in which e accuracy and stability of the Newmark method depends on β and c. In this study, β � 0.25 and c � 0.5 are used, which is known as the average acceleration Newmark method. In addition, the force in the above equation (29) can be either point force or pressure.

Results and Discussions
In this section, a series of numerical results for dynamic behaviors of the stiffened plate with general boundary restraints are carried out based on the theoretical model established in Section 2. Further analysis will be given in the following sections. Particularly, unless otherwise stated, the beam and the plate are assumed to have the same material properties, i.e., E p � E bi � 2.07 GPa, ρ p � ρ bi � 7800 kg/m 3 , μ p � μ bi � 0.3, and η � 0.01. It should also be noted that for the convenience of research, two reinforced plates of different geometric parameters are included in the current research, which are rectangular plate and square plate, respectively. For a rectangular plate, its geometric parameters are b � 1 m, a/b � 2, and h � 0.01 m, where the measured points can be expressed as follows according to the o-x-y coordinate in Figure 1 In the following analysis, the geometric and material parameters of the two plates are not repeatedly described. It is worth emphasizing that the response at the above-measured points is the lateral vibration response of the stiffened plate.
In addition, for the sake of simplicity, the symbols of C, S, F, and E are utilized to represent the clamped, simply supported, free, and elastic boundary conditions, respectively, all of which can be achieved by changing the spring stiffness. e spring stiffness values for different boundary conditions used in this paper are presented as follows: At x � constant or y � constant C: k u � 10 12 N/m, k v � 10 12 N/m, k w � 10 12 N/m, and K w � 10 12 N·m/rad. S: k u � 10 12 N/m, k v � 10 12 N/m, k w � 10 12 N/m, and K w � 0 N·m/rad. F: k u � 0 N/m, k v � 0 N/m, k w � 0 N/m, and K w � 0 N·m/rad. E 2 : k u � 10 12 N/m, k v � 10 12 N/m, k w � 10 12 N/m, and K w � 10 2 N·m/rad. E 5 : k u � 10 12 N/m, k v � 10 12 N/m, k w � 10 12 N/m, and K w � 10 5 N·m/rad. E 7 : k u � 10 12 N/m, k v � 10 12 N/m, k w � 10 12 N/m, and K w � 10 7 N·m/rad. E 55 : k u � 10 5 N/m, k v � 10 5 N/m, k w � 10 5 N/m, and K w � 10 5 N·m/rad. E 77 : k u � 10 7 N/m, k v � 10 7 N/m, k w � 10 7 N/m, and K w � 10 7 N·m/rad.
Another special note is the type and magnitude of the external load used in this paper. e excitation forces applied in the current research are point excitation and uniform pressure in the normal direction of the stiffened plate, which are denoted by the symbol of q w and p w , respectively. And the location of the concentrated point force is at point Q in Figure 2(a), when locally acting area of uniform pressure is shown by Ω in Figure 2(b). What needs special mention is q w � 1 and p w � 1 in the following calculation.

Convergence and Correctness Analysis.
In this section, the stiffened rectangular plate with CCCC case is first selected to check the convergence of the current method. e cross-section parameters of the beam are b b /h � 1 and h b / h � 1. Figure 3 shows convergence of velocity response for the rectangular plate with one central x-wise beam and one central y-wise beam under q w � 1, when the similar curves for the same coupling structure under p w � 1 are presented in Next, a square plate with one beam located at its diagonal line is utilized to implement a comparative study between the current method and FEM result ( Ansys15.0) . It is clarified that the configurations of the reinforcing component are e comparison examples are depicted in Figures 5∼6 which present the comparison of velocity response for the stiffened square plate under q w � 1 and under p w � 1, respectively. It is easy to be observed from the two figures that the current calculation method and the finite element result are in good agreement, which proves that the construction of the steady-state response dynamic analysis model by IFSM is correct. Simultaneously, compared with the finite element method, element mesh is not required in the current calculation model, which greatly saves the memory resources and improves the solving efficiency.

Structural Mobility Analysis.
In this part, the effect of the stiffening beam with different lengths on mobility characteristics of the reinforced plate subjected to SSSS will be investigated at first, where q w � 1 is adopted. As shown in Figure 7, four different lengths of beams are selected for the analysis of the mobility characteristics of the coupled structure, including L 1 � 0 m, 0. Shock and Vibration 9 plate will increase, which may be because the beam changes the local stiffness and mass distribution of the plate; Secondly, when the reinforcing beam passes through the transfer point D or the drive point C corresponding to L 1 � 0.5 m or L 1 � 1 m, the peak of their mobility curves will shift to the high frequency and magnitude of some peaks will appear to decrease in the range of 200-800 Hz, especially L 1 � 1 m, which shows that the rib increases the ratio of stiffness to mass of the plate in high frequency. e above phenomena also reveals that in structural vibration control, reinforcement helps to reduce the response peak of the plate in the midhigh-frequency range and changes the structural resonance frequency, but has a little effect on the response in low-frequency domain. Figure 8 studies the mobility characteristics of the diagonally stiffened square plate subjected to local uniform pressure with various elastic boundary conditions (i.e., E 2 SE 2 S, E 2 E 5 E 2 E 5 , E 5 E 5 E 5 E 5 , and CE 5 CE 5 ). e points C and D are selected as the observed points, and a local uniform pressure p w acts on the coupling system containing two reinforcing beams. For convenience of explanation, (x 0 , y 0 , φ) in Figure 1 and h b2 /h � 1. It can be seen from Figure 8 that as the boundary restrains vary from E 2 SE 2 S to CE 5 CE 5 , the resonant peak of the structural mobility response moves to the high-frequency direction and the peak value of the response curves will show a downward trend, especially under CE 5 CE 5 . erefore, it can be known that the frequency position corresponding to the vibration formant can be adjusted by modifying the boundary constraint spring stiffness to achieve structural vibration control. Another interesting finding is that the number of formants of the mobility responses increases when the stiffened plate is subjected to E 2 E 5 E 2 E 5 and CE 5 CE 5 boundary conditions. In order to investigate the effect of coupling spring stiffnesses on the structural mobility characteristics, a  stiffened rectangular plate with E 2 E 2 E 2 E 2 case is presented in Figure 9. e analytical model used in Figure 9 is a rectangular plate with one central x-wise beam and three y-wise evenly distributed beams. e geometric parameters of the reinforced beam used are L 1 � 2 m in the x direction, And point A and point B are selected as the drive point and transfer point, respectively. In addition, it should be noted that the coupling springs (including k pb1 , k pb2 , k pb3 , K pb1 , K pb2 , and K pb3 ) shown in Section 2 are replaced by the symbol of k b . e curves in Figure 9 show that with the increase of the coupled spring stiffness, the resonant peaks of the structural velocity mobility responses move faster and faster in the high-frequency direction. And the resonant peaks of the structural mobility move faster in the range of 150 Hz to 500 Hz than the range of 0 to 100 Hz, which indicates that the structural mobility responses in the high frequency are more sensitive to the coupled springs. In addition, we can also see that an increase in coupling spring stiffness will result in an increase in the peak values of the mobility response in the range of 0 to 150 Hz, which indicates that the reinforcement of the beam and the plate will have a negative effect on the low-frequency vibration reduction.
In the analytical model constructed in this paper, structural damping is introduced by means of complex Young's modulus. Damping effects on the mobility characteristics of a rectangular plate stiffened by three y-wise evenly distributed beams are illustrated in Figure 10. Except for removing a beam in the x direction, the remaining geometric parameters of the current model with E 5 E 5 E 5 E 5 boundary restrains and measured points follow the same values of the cases studied in Figure 9. It can be seen from Figure 10 that due to the existence of structural damping, the peak values of the vibration mobility responses are reduced. It is found after further observation that the magnitude of resonance peaks for mobility responses shows a decreasing trend as the structural damping increases, and the peak reduction in the range of 150∼500 Hz is more obvious than that in the range of 0∼150 Hz. Besides, the structural damping of the stiffened plate can only effectively reduce the mobility response amplitude of the resonance peak, while the effect of structural damping on the mobility response in the nonresonant region is not obvious.

Structural Intensity Analysis.
In the previous section, the mobility characteristics of the stiffened plates are clearly described. Undoubtedly, these numerical results provide reference data for further research. However, for the specific structural design, it is extremely necessary to understand the distribution and transmission of vibration energy in the structural system. To this end, the following analysis of structural intensity will be performed to show the flow strength and direction of the vibrational energy in the structure.
Firstly, structural intensity of a rectangular plate attached by beams with different configurations under twopoint excitation is carried out, where four symmetrical boundary conditions are employed, namely, CSCS, SFSF, E 77 E 77 E 77 E 77 , and E 55 E 55 E 55 E 55 . e geometric parameters of the beams (see Figure 11) are L 1 � L 2 � L 3 � 1, h b1 /h � 2, h b2 � h b3 � h b1 , and b b1 � b b2 � b b3 � h. In this example, regardless of the arrangement spacing of the beams, the orientation φ of all beams and their coordinates y 0 in the y direction (see Figure 1(b)) are the same, φ � 90 ∘ and y 0 � 0. erefore, the coordinates x 0 in the x direction will be used Mobility (dB) 200 400 600 800 0  Figure 11 indicate that as the stiffness of the boundary spring weakens, the presence of the reinforcement will intensify the effect of the vibrational energy gathering toward the edges of the plate and the reflection and backflow of vibrational energy occur at the edges of the plate. Moreover, regardless of the boundary conditions, the increased spacing of the ribs will have a positive increasing effect on the structural intensity of the plate between the reinforcements, especially for elastic boundary cases (E 5 E 5 E 5 E 5 and E 7 E 7 E 7 E 7 ). e above phenomenon may be because reinforcements hinder energy flow on the plate, resulting in energy accumulation in some areas. By further observing, it is found that in the case of symmetrical boundary cases and symmetrical point forces, there is always at least a section where the vibrational energy cancels each other out, which is called the power-insulation section. Furthermore, the number of such power-insulation section can be increased by setting an appropriate reinforcement spacing. Also, it can be intuitively found that the vibration energy is transmitted from the excitation to the periphery of the plate, but the force source is not always output source of the vibration energy (see Figure 11(d)). Next, structural intensity of an SSSS square plate stiffened by one beam with various orientations is shown in Figure 12 Figure 8. Moreover, the excitation frequency is 122.44 Hz (see Figure 12(A∼D)) and 195.90 Hz (see Figure 12(E∼H)), respectively, which represent the 2 nd frequency and 4 th frequency for the square plate. From Figure 12, it is no difficult to find that as the azimuth angle of the reinforcement varies from 0 ∘ to 45 ∘ , the energy flow near the excitation source S first decreases and then increases and the structural intensity reaches a minimum when the reinforcement passes through the excitation source (that is, φ � 20 ∘ ). Besides, there are significant differences in the influence of the azimuth variation of the reinforcement on the structural density distribution of the plate at different excitation frequencies, which reflects that the effect of the reinforcement on the structural density of the plate depends dramatically on the excitation frequency.

Transient Response Analysis.
is section is concerned with the transient vibration analysis of the reinforced plate under point force and local uniform pressure. In addition, it can be seen through further observation that the present results are highly consistent with FEM results, which verifies the correctness of the transient response model of the stiffened plates established in this paper. In order to further verify the correctness of the damping treatment method in this paper, a comparison between the current method and FEM results is implemented in Figure 15, where only the Rayleigh damping is added and the other parameters are consistent with these from Figure 13. And the selected Rayleigh damping coefficients are set as α 0 � 0.999 and β 0 � 0.001. As can be seen from Figure 15, a good agreement between the current results and those obtained by FEM is achieved, which indicates that the presented damping treatment of stiffened plates is reasonable.
Next, transient responses of an SSSS square plate stiffened by a beam with various orientations are presented in Figure 16. For the present analysis, the time increment is taken to be 0.33 ms and normal point load q w located at point C is applied. e geometric parameters and arrangement of the beam are the same as those in Figure 12, but the arrangement of the beam with (0, 0, 20 ∘ ) is replaced by (0, 0, 15 ∘ ) here. e curves in Figure 16 reveal that regardless of the measured point A or B, when the orientation of the reinforcing beam changes from 0 ∘ to 45 ∘ , the oscillation periods of the time response curve are shortened and * * * * * * * * the amplitude of the displacement responses shows a decreasing trend.
is is undoubtedly expected because the closer the stiffening beam is to the measured points, the greater the bending stiffness at the measured points is.
As mentioned in Section 2.5, the Rayleigh damping is selected to implement the transient vibration analysis of the beam-stiffened plate and it can be obtained by setting the values of α 0 and β 0 . erefore, the influence of the Rayleigh damping coefficients on transient performances of the stiffened plate will receive attention in the following example. As shown in Figure 17  parameters of the used beams are set as In this investigation, the same form and incentive position of normal load as Figure 17 is utilized, where the observed point is point C. From Figure 17, it can be intuitively seen that due to the existence of damping, the vibration energy of the coupling structure will be continuously consumed and the transient vibration of the coupled structure exhibits a completely different behavior than the undamped transient vibration. Secondly, the change of the mass damping coefficient α 0 has a little effect on the transient vibration of the stiffened plate (see Figure 17(a)). However, a small change in the stiffness damping coefficient β 0 will cause a significant change in the transient response (see Figure 17(b)). Specifically, when the parameter β 0 is smaller than a certain value denoted by the symbol β 1 0 , the transient amplitude of the coupled structure will gradually decrease and finally stabilize. When β 0 exceeds the value β 1 0 , the vibration response of the reinforced structure will gradually increase and finally stabilize at a certain value. Also, as the value of β 0 increases, the time for the transient response to reach the steady state will be extended.  In the last transient study, transient responses of a rectangular plate with three x-wise and one y-wise evenly distributed beams are examined in Figure 18, in which normal local uniform pressure (p w �1) is considered. Except for the difference of the arrangement of the beams, the parameters of the beams employed in Figure 3 are used here and the measured point is point B. What needs to be explained here is that when the same type of springs at edges of the stiffened plate is studied, the stiffness values of the remaining types of springs are set to infinity replaced by 10 12 ; for instance, while the elastic parameters k w varies from 0 to 10 12 , the stiffness values of k u , k v, and K w are designated as 10 12 . Figure 18 shows that as the stiffness value of transverse springs k w varies from 10 5 to 10 11 , the amplitude and phase of the lateral displacement response will change significantly. Specifically, as the stiffness value increases, the oscillation period of the response curve will be shortened and the vibration amplitude will decrease. When the stiffness value of spring k w is greater than 10 11 , the transient displacement response remains unchanged and is close to the response under CCCC case. For a rotating spring K w , the lateral displacement response will only change significantly if its stiffness value varies from 10 3 to 10 8 , but its effect is not as obvious as that of the transverse spring k w . In addition, it can be seen from Figures 18(a) and 18(b) that the lateral displacement response is not affected by the variation of the in-plane linear spring k u and k v , which is consistent with the current classical thin plate theory. In addition, the results in Figure 18 indicate that the definition of the boundary conditions at the beginning of this section is reasonable.

Conclusions
An analysis model is established for investigating the linear dynamic analysis of rectangular plates reinforced by beams of any lengths and orientations. Employing IFSM, displacement functions of the plate and reinforcements are obtained, respectively. Based on Rayleigh-Ritz method, dynamic response of the coupling structure is obtained, in which various boundary restrains and beam-plate coupling relationships are achieved using the artificial spring technology. Several numerical analyses verify the effectiveness of the current method, and a series of novel numerical results are also given. Besides, it should also be noted that even though only the linear dynamic performances of the beam-stiffened plate are reported here, the presented methods can also be extended to perform other dynamics analysis related to the stiffened plate, such as the vibration control of the stiffened plate in electrothermal-magnetic field and vibration analysis of reinforced coupling plates.
In addition, there are some conclusions that are stated as follows: (1) e mobility characteristics of the beam-stiffened plate are not only related to the length of the beams, but also to coupling spring stiffness and structural damping, especially in the high-frequency range, which provides a design basis for structural vibration control. (2) e arrangement spacing and orientation of stiffening beams have a significant effect on the structural intensity distribution of the plate. Moreover, the structural intensity distribution of the stiffened plates is sensitive to the frequency of the excitation force and boundary conditions, which give a reference for discriminating the position of the load and the transmission path of the vibration energy. beams or increasing the spring stiffness (k w and K w ) will have a negative effect on the amplitude and oscillation period of the transverse transient response of the measured point. However, the in-plane springs (k u and k v ) have no effect on the lateral transient vibration of the stiffened plate. (4) When the Rayleigh damping is considered in the transient analysis of the stiffened plate, the transient response is heavily affected by the stiffness damping coefficient rather than the mass damping coefficient.

Data Availability
All the underlying data related to this article are available upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.  Shock and Vibration