Research on modeling joining and joining modification method of hybrid FE-SEA model

Since the connection method at the boundary of the hybrid FE-SEA (finite element-statistical energy analysis, FE-SEA) model affects the overall calculation accuracy of the model, the missing or inaccurate connection relation will lead to a large error in the model calculation. In order to effectively solve the calculation error caused by the connection problem of hybrid FE-SEA model, this paper starts from the modeling methods and connection relations of “hybrid-point” connection, “hybrid-line” connection and “hybrid-surface” connection. In order to solve the modeling method and correction problems of the “hybrid-point” connection in the hybrid FE-SEA model, the “out-of-plane” wave motion equation at the “hybrid-point” connection was established by using the superposition principle of plane waveforms in polar coordinate system. The “hybrid-point” connection and wave-number relation in Cartesian coordinate system are studied. The radiation radius correction method of “hybrid-point” connection is proposed. An example is given to verify the effectiveness of the method. In order to solve the modeling and correction problem of “hybrid-line” in hybrid model, a “hybrid-line” connected triangular waveform function model was established by using the method of linear difference. The direct field dynamic stiffness matrix of “hybrid-line” connection in node coordinate system is studied. According to the shape function of the “hybrid-line” connection wavenumber space, a method to correct the shape function of the “hybrid-line” connection is proposed, and the validity of the method is verified.


Wave theory
According to the hybrid FE-SEA theory, the response of the stochastic subsystem is composed of the direct field and the reverberation field.The energy reciprocity principle exists between the direct field and reverberation field, and the average response of the subsystem can be solved by wave theory and modal method.The direct field is defined as the physical field in which the energy input to the subsystem does not pass through any boundary reflection.The reverberation field or diffusion field is defined as the physical field in which the energy input to the subsystem is reflected many times.The hybrid FE-SEA model carries out the system energy transfer and exchange between the "direct-reverberation field" through different connection methods.A simple composite system (consisting of two thin plates embedded in a beam) is used to describe the direct field-reverberation reciprocity relationship of the hybrid model system.As shown in Fig. 1, the most direct way to analyze the dynamic response of a system is to generate a finite element model.The degree of freedom consists of a set of displacements distributed throughout the system.Thin plate structures have relatively short vibration wavelengths, the idea to solve this problem can be realized by considering a single thin plate.As shown in Fig. 2, the displacement degree of freedom on the boundary is q.For harmonic motion of frequency ω , the relationship between displacement degree of freedom on the boundary and a set of external forces f acting on the boundary can be expressed as follows: where D is the dynamic stiffness matrix, q is the displacement degree of freedom on the boundary, and f is the external force at the boundary.The dynamic characteristics of the whole system can be obtained from the dynamic characteristics of the composite beam frame of two thin plate structures in Fig. 1.The harmonic motion at the boundary is realized by wave theory to generate elastic waves.These elastic waves propagate throughout the system and bounce back each time a boundary is encountered.
The direct field dynamic stiffness matrix D dir is defined as the matrix D without reflection.The difference between the boundary force D dir q generated by this matrix and the actual boundary force Dq is called the "rever- beration" force f rev = D dir q − Dq .Then, Eq. (1) can be expressed as: Assuming that the reverberation wave field is a "diffusion field", then: (1)  where E[] represents the average value of the random structure, E is the vibration energy of the plate, and n is the modal density of the plate.Equation ( 4) is called "diffusion field reciprocity relation".However, there are cases where components may be coupled to other components across domains, rather than just at the boundary, where a modal approach is needed for synthesis.

Modal method
Based on modal theory, a modal parameter method is proposed.This method is used to describe the degree of freedom q of the thin plate in Figs. 1 and 2 as a subset of all the degrees of freedom.The direct field dynamic stiffness matrix D dir is replaced by the statistical mean of the dynamic stiffness matrix E[D] .The reverberation force vector f rev is replaced by a randomly varying force component f ran .Therefore, f ran = E[D]q − Dq , Eq. ( 2) can be expressed as 21 : Suppose, the degree of freedom q is at the component boundary, then the relation E[D] = D dir is obtained.In addition, the receiver matrix H dir is used to represent the inverse matrix of D dir .According to these defini- tions, if the natural frequency and mode shape of the component conform to Gaussian orthogonality ensemble (GOE), then: where m = ωηn represents the modal overlap of components, and E 1 represents the exponential integral.

The equations of motion of point connection in hybrid model
Point connection is one of the important connection methods of hybrid model.Point connection generally exists between thin plate structure and beam structure.Energy transfer between hybrid models is carried out through nodes.Its connection relationship is shown in Fig. 3. Figure 3a shows the point connection relationship between plate and beam structure, and Angle α represents the radiation Angle of point connection.Figure 3b represents the point connection mass and moment of inertia of the two connection structures, and M represents the connection mass at the point connection.Xg, Yg and Zg respectively represent the coordinate system at the point connection, and Ixx, Iyy and Izz respectively represent the moment of inertia in three directions at the point connection.Figure 3c shows the point-connection offset distance of the two-connection structure.dx, dy and dz respectively represent the offset distance of the point connection structure on the three coordinate axes X, Y and Z. Figure 3d shows the stiffness effect of the point connection structure, and B1 and B2 represent two independent subsystems respectively.Kx and Kz represent the connection stiffness on the X axis and Z axis respectively.The local coordinate system of point connection is established, as shown in Fig. 4. It is assumed that the point connection mode of infinite thin plates is a massless rigid circular plane.The infinite thin plate contains three kinds of propagation waveform and three kinds of dissipation waveform: plane displacement wave ω , plane bending wave θ x and plane bending wave θ y .And the tensile and shear waves associated with the in-plane displacements u and υ .And the shear wave associated with the Angle θ z .Under the condition of polar coordinates (r, θ, z) with point connection as the origin, the out-of-plane displacement of an infinite thin plate can be expressed as the superposition of six waveforms related to ω , θ x and θ y , then it can be deduced 6 : where u z is the plane external displacement of the thin plate in the polar coordinate system, and H (2) n is the Hankel function.a 1c are the amplitudes of the dissipative wavelength respectively.The wave number in the bending direction of the sheet is k b .The in-plane motion of a thin plate can be described as a tensile wave in-plane and a shear wave in-plane, and the relationship is as follows 20 : where a  1 (k e r) k e r ] cos(θ) 1s H (2) 1 (k e r) cos(θ) − a (e) 1c H (2) 1 (k e r) sin(θ)] + a (s) 0 H (2)   www.nature.com/scientificreports/associated with u and υ , respectively.k e and k s are respectively the tensile and shear wave numbers of the thin plate.

The relationship between point connection and wave number in hybrid model
According to the bending theory of circular thin plate, the expressions of bending moment M x , M y and external force F z in Cartesian coordinate system are obtained as follows 21 : where in the Cartesian coordinate system, the expressions of external forces F x ,F y and bending moment M z are as follows: Substituting the above Eqs.(10), (11) and (12) into Eqs.(13) to (18), the displacement equation of thin plate structure expressed by waveform coefficient can be obtained as follows.Where S represents the transverse force.

Line connection of hybrid models
A line connection is a form of connection used to describe a thin plate component in a hybrid model.If the coupling between the structural subsystem and the acoustic cavity subsystem is not considered in the system, the connection mode of the system is mainly line connection and point connection.In general, the line connection mode describing the system is shown in Figs. 5 and 6. Figure 5 shows the local coordinate system ( x j , y j , z j ) of the line connection of the j-block plate of the hybrid model.The connection is located at ( y j = 0 ).Assuming ( j ) is a thin plate, the equation of motion at the line junction can be described by translational motion in three directions ( (u j , υ j , ω j ) T ) and rotation ( θ j ) about the X-axis.According to the isotropic theory of thin plate, the governing equation of transverse free motion of thin plate can be described as 22 : (13) where, , D o is the bending stiffness of the thin plate.m is the mass per unit area.According to the Snell theorem (as shown in Fig. 6), the transmitted wave and the reflected wave have the same spatial wave term.The waveform of the structure at the junction line has the same trajectory velocity to satisfy the phase compatibility condition at the junction.Assuming that both the transmitted and reflected waves are bending waves, the wave term is exp(−ikx + µ B y + iωt) .Where µ B is the undefined wavenumber in the y direction, and k is the wavenumber in the x direction.If the Angle between the incident wave and the normal is set as φ ′ and the wave number is k ′ B , then k can be expressed as 21 : Substitute exp(−ikx + µ B y + iωt) into Eq.( 20) to get µ B .The following relationship is satisfied: where k B is the bending wave number.If the two values of µ B are µ B1 and µ B2 respectively, then the transverse displacement of the thin plate can be expressed as: The Angle θ x in the x direction is expressed as θ x = ∂ω ∂y .Substitute it into Eq.( 23), and the following relationship can be obtained: Considering the boundary, y = 0, and Eqs. ( 23) and ( 24), the following relation can be obtained: By omitting exp(−ikx + iωt) , a 1 and a 2 can be obtained from Eq. (25) : According to the classical thin shell plate theory, the X-axis moment and the transverse force can be expressed as:  26) and ( 27), the following relationship can be obtained: where , the dynamic stiffness matrix of transverse vibration of thin plate line connection can be written as:

Surface connection of hybrid models
In the hybrid model, the boundary between the plate and the volume of the cavity is related to the out-of-plane displacement of the whole plate.This coupling connection mode is called "surface connection" or "area connection", and the diagram of surface connection relationship is shown in Fig. 7. Figure 7a shows the diagram of surface connection between the thin plate and the two sound cavities, where L represents the distance between the two sound cavities.TL represents the acoustic transmission loss between two sound cavities.R represents the radiation efficiency between the thin plate structure and the two sound cavities.NR stands for noise emission between two cavities.Figure 7b shows a schematic diagram of the surface connection between the panels.The arrows in the figure represent the sound wave transmission at the interface between plates, P represents the thin plate structure, and A represents the sound wave of the thin plate structure.When the subsystems in the hybrid model have plane connections, the boundary degree of freedom Q can be defined as the out-of-plane response covering the point grid on the thin plate surface.In order to calculate the system response of Eq. ( 4), it is necessary to identify the dynamic stiffness matrix D dir .In this case, the jk term of the receive matrix can be expressed as: where r jk is the distance from grid node J to K, and G is Green's function of an infinite plate.The dynamic stiffness matrix D dir is obtained by the inverse matrix of H dir .Under the condition that the thin plate is directly excited, it can be deduced from Eq. ( 5) : where the first term on the right side of the equation is the response of the infinite thin plate, and the second term is the additional response caused by the system.For weakly damped systems, the second term dominates.By combining Eq. (35) with Eq. ( 6), the following relation can be obtained 21 : In Eq. ( 6), assuming that the second term on the right side of the equation can be ignored, the following matrix identity applies to any matrix A: Therefore, Eq. ( 36) can be transformed into: Equation (38) represents the correlation between the response at node j and the response at node k.It's proportional to the imaginary part of Green's function between two points.Based on the above equation, the system response at the hybrid surface joint can be determined.The degree of freedom in the hybrid model is composed of a free grid q i connected by each region.It is expressed as a deterministic set of degrees of freedom q m that forms associated components with deterministic finite element (FE) models.A particular collection q i has multiple SEA subsystems.The thin plate and the sound cavity share a common grid, so that the general structure of the equation has the form of a complete set of determined degrees of freedom 20-22 : where D m is the dynamic stiffness matrix of the (FE) model, including the "direct field" dynamic stiffness matrix of any subsystem connected to the model by non-regional nodes.Ni is the number of subsystems connected to the grid of the i node, and D i dir,j is the dynamic stiffness matrix of the jth subsystem connected to the ith grid node in the direct field.C i is the deterministic coupling matrix (assuming no power is consumed in the coupling, the coupling is the Hermitian matrix).F i represents the external force acting directly on q i degree of freedom and the reverberation force generated by the additional subsystem.One of the solving difficulties of the current method is that q i may contain multiple degrees of freedom, so it can be further reduced by directly applying the hybrid method in Eq. (39) to obtain: Further, the matrices involved are reduced by rearrangement: Therefore, Eq. (40) can be written as: Equation ( 44) is a simplified form of deterministic system equation.Where, the degree of freedom of q 1 does not reflect, and the energy flow generated by the reverberation field of subsystem j is considered.According to Eq. ( 44), the deterministic subsystem response F 1 is provided by the following equation: According to the principle of mutuality of the hybrid model, the following can be obtained 21 : www.nature.com/scientificreports/ The above results can be rewritten as: where, the reduced direct field dynamic stiffness matrix is expressed as: Assuming that the dynamic stiffness matrix is symmetric, it can be concluded that: Therefore, the input power of subsystem k caused by the response of Eq. ( 48) is: According to Eq. ( 40) and the condition that F 1 = 0 , the power equation is: According to the above method, the standard hybrid equation can be realized only by using deterministic DOF without solving the grid DOF of the reference region.The direct field dynamic stiffness matrix associated with subsystem j is determined by Eq. (42).This method can greatly reduce the scale of the hybrid surface connection matrix.

The existing hybrid model connection correction method
Based on the wave theory and the assumption of infinite thin plate structure, the boundary of the direct field dynamic stiffness matrix of point-line connection in the hybrid model needs to be evaluated by parameter correction.Generally, the modified radiation Angle α is used to describe the influence of point-line connection on boundary to boundary 6 , and the radiation Angle correction diagram is shown in Fig. 8. Figure 8a shows the radiation Angle diagram of the beam structure and the thin plate structure.The parameter α represents the radiation Angle at the point junction, and the connection relation attributes can be corrected by modifying the radiation Angle at the point junction.Figure 8b shows the schematic diagram of radiation Angle of line connection between thin plates.The radiation Angle between the two plates is α.X, Y and Z respectively represent the connection coordinate system at the line connection, and dz and dy represent the offset distance between the two thin plate structures and the line connection.The radiative impedance matrix of point-line connection in the hybrid model after parameter modification can be expressed as 20 : where β α is the correction factor of the point-line connection mode of the hybrid model.The subscript α repre- sents the correction factor defined by the radiation Angle.Modifying the radiation Angle can modify the radiation impedance matrix of the point-line connection, so as to modify the precision of the point-line connection model.However, the radiation Angle correction method also has some limitations: 1)The radiation Angle parameter is defined as a constant term over the entire frequency range.The correction Angle with respect to frequency cannot be defined.It has certain limitations.
2)It is a difficult problem to determine the parameter range of radiation Angle when modeling the hybrid model with connection forms such as point and line.In the existing literature, the parameter selection range of radiation Angle has not been elaborated in detail.

Modification method of hybrid point connection
Considering the displacement of the circular surface connected by points at r = r 0 , the relation of displacement outside the plane in the polar coordinate system and the Cartesian coordinate system is expressed as follows: www.nature.com/scientificreports/where f is the frequency of the system, p is the density of the thin plate, h is the thickness of the thin plate, and k w is the bending stiffness.Further, the wavelength b in the thin plate structure system can be deduced, which can be expressed as: The point connection position is taken as the center of the circle, and the wavelength b of the thin plate structure is taken as the radiation correction radius to establish a circular region.The area where the thin plate structure region overlaps with the circular region is defined as the correction region of the point connection under the condition of frequency f.The correction factor of the radiation radius can be expressed as: where ⌢ β is the ratio of the dynamic stiffness matrix solved by the finite element method, which is used to correct the error term.Therefore, the direct field dynamic stiffness matrix of the modified hybrid connection based on the radiation radius is expressed as: The correction factor of the radiation radius of the hybrid connection is related to the location of the hybrid connection, the analysis frequency and the inherent physical properties of the thin plate structure.

Modification method of hybrid line connection
The dynamic equation was established based on the wavenumber space coordinate system.The direct field dynamic stiffness matrix of the coordinate system of nodes is obtained by Fourier transform.Let the displacement ω(x) = n a n u n (x) , where n = 1,2… N, u n (x) is a set of shape function values.Then the force generated by the node motion ω at the connection position x can be expressed in the convolution form 21,22 : where d(x − x ′ ) represents Green's function of the line connection model.Its physical meaning is the force gener- ated at x by applying a unit pulse displacement at position x ′ .Substitute the shape function ω(x) = n a n u n (x) into Eq.( 73), and the following relation can be obtained: The load obtained by Eq. ( 74) is projected onto the generalized coordinate a m .Then the generalized force under the projection can be expressed as: Substituting Eq. (74) into Eq.( 73), the following relation can be obtained: Therefore, the dynamic stiffness at the line connection is defined as: According to the above equation, it can be deduced: u n (x′) and d(x) are expressed in inverse Fourier transform form: where the value range of R is {−∞ ≤ K ≤ + − ∞} .Substitute Eqs. ( 79) and (80) into the expression of D mn and the following relationship can be obtained 22 : The shape function of the linearly connected triangular wave is further determined.The triangular waveform function obtained by linear interpolation is shown in Fig. 10.The lateral displacement at the connection can be expressed as: where u n (x) is the nth shape function, and a n is the participating factor of the shape function.The node dis- placement can be solved by linear interpolation.The mth shape function is defined as follows: the displacement of the (m-1)th node and the (m + 1)th node is zero, and the intermediate displacement is calculated by linear interpolation.Then, the triangular waveform function as shown in Fig. 11 is constructed, and its function expression is as follows: where l is the distance between nodes of the line connection model.The hybrid line connection model is modified by modifying the triangle wave function of the line connection in the hybrid model.Line connection correction of the hybrid model was carried out according to the empirical parameters of line connection of the shape function in Eq. ( 83) 20 .The direct field dynamic stiffness matrix of the model is modified by modifying the (75) where Sinc(x) = sin(x)/x , for the shape function of any node, the expression of wavenumber space can be obtained by Fourier transform: The shape function above is used to modify the parameters to establish the motion equation of the thin plate, and the displacement can be expressed as: According to the above shape function, the dynamic stiffness matrix of the direct field connected by the line in the node coordinate system can be obtained:

Algorithm verification analysis Validation analysis of hybrid point connection correction method
In order to verify the validity of the correction factor model of point connection, the structure shown in Fig. 12 was analyzed.The model is composed of a beam structure and two thin plate structures.The material properties of the hybrid model are shown in Table 1.The thin plates are connected to the beam structure through joints.The material of the beam structure is steel, the cross section is hollow square, the height of the beam is 25.4 mm, the thickness of the beam is 3.2 mm, and the combined structure is composed of 12 beams.The structure size of the thin plate is 600 mm × 800 mm, the thickness of the thin plate is 1 mm, the material is aluminum, and the analysis frequency is 1-1300 Hz.Loads are applied laterally to one of the plate structures.Since the point connection is located inside the thin plate structure, the radiation Angle (360°) of the point connection under ideal conditions is converted to a correction factor of 1.The connection method of the model is modeled by the hybrid point connection model.Then the radiation radius correction factor is used to correct it.In the process of hybrid model modeling, the in-plane motion of the thin plate structure and the beam structure are used to (FEA) determine the structural modeling.The out-plane motion of the thin plate is modeled by the statistical energy (SEA) due to the relatively small wavelength.Figure 13 shows the finite element model of the structure.Considering the uncertainty of the medium frequency response of the hybrid model, the finite element energy flow method and Monte Carlo test method were used as the reference solution of the numerical simulation.The method is used to simulate the perturbation of structural parameters by applying 20 random distributed masses on each thin plate with a total mass of 15%.The radiation radius method proposed in this paper is compared with the radiation Angle method of the existing technology method.The radiation Angle is selected by the setting method in Literature 22 .The radiation angles were 360°, 180° and 90°, respectively.According to the reference 20 and Eq.(69), the equivalent correction factors were β = 1, β = 0.5 and β = 0.25, respectively.Figure 14 shows the numerical simulation results of vibration velocity of plate 1 by different calculation methods.As can be seen from Fig. 14a, the error between the existing method (β = 1) and the reference solution is large in the frequency range of 800-1200 Hz, and the maximum error occurs in the frequency range of 800-1000 Hz.The proposed method is used to modify the point connection modeling, and the fitting degree between the modified curve and the reference solution is good in the whole frequency band.Within the range of 0-600 Hz, the method in this paper is basically consistent with the curve of reference solution.Within the frequency range of 800-1200 Hz, the error range between the proposed method and the reference solution is also within the reasonable engineering error range.The peak correspondence of the curve is good in the whole frequency range, which indicates that the correction factor method based on the radiation radius is more suitable for the correction of the intermediate frequency response of the hybrid model than the radiation Angle method of the prior art.Compared with the existing literature 20 , the correction accuracy of the proposed method for the whole frequency band is higher, and it is suitable for the point connection model correction of hybrid FE-SEA model.Compared with the reference 20 , the hybrid model point connection correction method proposed in this   www.nature.com/scientificreports/ of the correction factor, and makes a mathematical correction to the error amount of the correction factor.Figure 14b,c show the comparison between the existing method (β = 0.5, β = 0.25) and the proposed method.It can be seen from the correction results that although the radiation Angle correction by the existing method can reduce the response error of the whole system, the correction accuracy is still far from the method presented in this paper.In order to further verify the differences between the proposed method and the existing methods, Fig. 15 shows the numerical calculation results of vibration energy of plate 1 by different calculation methods.From Fig. 15a-c, it can be seen that the error between the existing method (β = 1, β = 0.5, β = 0.25) and the reference solution is large in the intermediate frequency range.The modified method of radiation radius in this paper is in good agreement with the Monte Carlo reference solution.In the low frequency range of 0-600 Hz, the proposed method has a high degree of curve fit with the reference solution.The peak-peak correspondence of the curve is good in the whole low frequency range.In the intermediate frequency range of 600-1200 Hz, the method presented in this paper is consistent with the peak trend of the reference solution curve.The peak-peak value is within a reasonable error range, and the simulation accuracy is higher than that of the existing methods.In this paper, the radius-based mixing point connection model modification method can be used as a mixing model modification method, which provides a solution for accurate mixing model modeling and modification.

Validation analysis of hybrid line connection correction method
The line connection of the hybrid model can be modified by modifying the shape function of the line connection.
In order to verify the effectiveness of the modified hybrid line connection model, the structure shown in Fig. 16 was verified.The thin plate structure of the hybrid line connection model is composed of two plate elements of the same size and connecting plate.Its material parameters are shown in Table 2.The sheet material is aluminum sheet, with a length of 600 mm, a width of 700 mm and a thickness of 1 mm.The connecting plate is made of steel structure, with the length of 600 mm, width of 70 mm and thickness of 4 mm.The analysis frequency range is 0-1000 Hz.Apply unit sweep excitation on plate 1. 20 nodes were randomly selected from each thin plate as Montecarlo response reference points 20 to solve the vibration velocity and coupling loss factor of the system.Considering the dynamic characteristics of the whole system, the displacement out-plane with high modal density is modeled by the statistical energy method.The finite element model is used for the joint plates with small modal density and large in-plane displacement and stiffness.The two hybrid models are coupled through two hybrid lines, and the hybrid model is shown in Fig. 17.
The reference solution is a hybrid finite element energy flow method combined with Monte Carlo test.The finite element method is a widely used dynamic analysis method at present.At the same time, Monte Carlo simulation is used to introduce parameter perturbation to describe the influence of uncertainty on intermediate frequency response.Since the finite element method does not use a lot of assumptions, the results of the energy flow analysis-Monte Carlo simulation method can be used as the basis for improving the statistical energy method or the prediction method of the medium frequency mechanical environment.In the existing methods, the setting method in the literature 22 is used to correct the radiation Angle, and the radiation Angle is 360°, 180° and 90° respectively.According to the literature 20 and Eq.(69), the equivalent correction factors are β = 1, β = 0.5 and β = 0.25, respectively.Figure 18 shows the vibration velocity curves of plate 1 calculated by different methods.As can be seen from Fig. 18a, the curve fitting degree between the proposed method and the Montecarlo reference solution is relatively high in the frequency range of 0-1000 Hz.The peak-peak correspondence of the curve is good in the frequency range of 0-500 Hz, and the maximum error of the whole frequency band is less than 10 -1 .Compared with literature 20,21 , the method in this paper has a higher degree of curve fitting in the middle frequency band.Figure 18b shows the comparison between the existing method (β = 1) and the reference solution.It can be seen from the figure that the existing method (β = 1) has a large error with the reference solution in frequency bands of 100 Hz, 370 Hz, 700 Hz, etc., and the peak-peak correspondence of the curve is poor.Similarly, the existing methods (β = 0.5, β = 0.25) in Fig. 18c,d have large errors in multiple peak frequencies with the reference solution, and the curve fitting degree is not high.
In order to further verify the modified method of hybrid line connection, the coupling loss factor is verified for the hybrid model of line connection.Figure 19 shows the calculation results of coupling loss factors of plate 1 with different calculation methods.It can be seen from Fig. 19a that the proposed method has a high degree of fit with the Montecarlo reference solution.In the frequency range of 0-500 Hz, the peak value correspondence of the coupling loss factor curves obtained by the proposed method and the reference solution method is good, and the maximum error is less than 10 -1 , which shows a high precision of linear connection calculation.When the frequency range is 500-1000 Hz, the other frequencies have good calculation accuracy except for the weak fluctuation of individual peaks.Compared with reference 20,21 , the modeling and correction accuracy of hybrid line connection using the method in this paper can effectively improve the calculation accuracy of Intermediate Frequency and the corresponding frequency of peak-peak value of curve.Based on reference 20,21 , this paper synthesizes the existing theoretical methods of radiation Angle correction and shape function correction.The conversion method of radiation Angle is optimized, the modified method of mixed line connection based on form function is optimized, and the node form function and the node dynamic matrix based on form function are optimized, which are further optimization of the existing theoretical methods.Figure 19b is a calculation comparison between the existing method (β = 1) and the reference solution.It can be seen from the figure that the existing method (β = 1) has poor correspondence with the trend of the reference solution curve, and the peak error in the whole calculated frequency band is large.Similarly, Fig. 19c,d show a comparison between the existing method (β = 0.5, β = 0.25) and the reference solution.It can be seen from the figure that the correction error of the line connection model of the radiation Angle in the existing method is large, and it is not consistent with the curve of the reference solution, and the peak-peak frequency does not have good correspondence in the whole calculated frequency band.The above numerical examples show that the existing methods have large  with the existing methods, it has a good practical value and can provide a solution for the modification of hybrid wire connection model.

Conclusion
In order to solve the problem of the modeling method and correction of the "hybrid-point" connection in the hybrid FE-SEA model, the "out-plane" tension wave and shear wave motion equations of the "hybrid-point" connection are established by using the superposition principle of plane waves in polar coordinates.The relation between "hybrid-point" connection and wave number in Cartesian coordinate system is studied.According to the wave propagation characteristics of thin plate structures, a "hybrid-point" connection radiation radius correction method is proposed.Compared with the prior art, the modification of the hybrid FE-SEA model with the radiation radius correction method is in good agreement with the Monte Carlo reference solution results.An example study shows that the error between the existing method (β = 1) and the reference solution is large   in the frequency range of 800 Hz to 1200 Hz.The proposed method is used to modify the "hybrid-point" connection modeling, and the fitting degree between the modified curve and the reference solution is good in the whole frequency band.In the range of 0-600 Hz, the method in this paper is basically consistent with the curve of reference solution.In the frequency range of 800-1200 Hz, the error between the proposed method and the reference solution is within a reasonable range.Compared with reference 20,21 , the correction accuracy of the proposed method for the spectrum in the whole frequency band is higher.In order to solve the modeling method and correction problems of the "hybrid-line" connection in the hybrid FE-SEA model, the triangular waveform function model of the "hybrid-line" connection was established by using the linear difference method.According to the "hybrid-line" connection shape function in wavenumeral space, a modification method of "hybrid-line" form function is proposed.The numerical example shows that the modified method based on the shape function has a good fit with the Monte Carlo reference solution.In the range of 0-500 Hz, the method presented in this paper has a good correspondence with the peak value of the coupling loss factor curve obtained from the reference solution, and the maximum errors are all less than 10 -1 .In the frequency range of 500-1000 Hz, except for the weak fluctuation of individual peaks, the other frequencies have good calculation accuracy.Compared with reference [20][21][22] , the accuracy of the "hybrid-line" connection modeling and correction is effectively improved.
1c are the amplitudes of the propagation waveform associated with (ω, θ x , θ y ) respectively.a(bn) 0 , a (bn) 1s and a (bn) 1s are the amplitudes of the plane tensile propagation waveform related to u and υ respectively.a (s) 0 is the amplitude of the shear plane wave in-plane.a (s) 1c and a (s) 1s are the amplitudes of the plane shear waveforms

Figure 3 .
Figure 3. Schematic diagram of point connection relationship.

Figure 5 .
Figure 5. Local coordinate system of line connection mode.

Figure 10 .
Figure 10.Graph of displacement function of line connection.

Figure 11 .
Figure 11.Shape function of the line connection.
https://doi.org/10.1038/s41598-023-43250-3www.nature.com/scientificreports/errors in the correction of the hybrid line model, and the selection range of radiation Angle correction by the existing techniques has certain limitations.The method in this paper is modified based on the shape function of the hybrid line connection model, so as to modify the direct field dynamic stiffness matrix at the joint.Compared

Figure 15 .
Figure 15.Vibration energy of plate 1.(a) Comparison of the improved method with the existing method (β = 1).(b) Comparison of the improved method with the existing method (β = 0.5).(c) Comparison of the improved method with the existing method (β = 0.25).

Figure 18 .
Figure 18.The vibration velocity calculated by different methods.

Table 1 .
Hybrid model material parameters.

Table 2 .
Hybrid model material parameters.