Hermitian Mindlin Plate Wavelet Finite Element Method for Load Identification

A new Hermitian Mindlin plate wavelet element is proposed. The two-dimensional Hermitian cubic spline interpolation wavelet is substituted into finite element functions to construct frequency response function (FRF). It uses a system’s FRF and response spectrums to calculate load spectrums and then derives loads in the time domain via the inverse fast Fourier transform. By simulating different excitation cases, Hermitian cubic spline wavelets on the interval (HCSWI) finite elements are used to reverse load identification in the Mindlin plate. The singular value decomposition (SVD) method is adopted to solve the ill-posed inverse problem. Compared with ANSYS results, HCSWIMindlin plate element can accurately identify the applied load. Numerical results show that the algorithmofHCSWIMindlin plate element is effective.The accuracy ofHCSWI can be verified by comparing the FRF of HCSWI and ANSYS elements with the experiment data. The experiment proves that the load identification of HCSWI Mindlin plate is effective and precise by using the FRF and response spectrums to calculate the loads.


Introduction
For many practical engineering problems, such as structural strength analysis, health monitoring, fault diagnosis, and vibration isolation, the knowledge of the dynamic loads acting on the structure is always required and important [1].The use of plate and shell in the industry becomes increasingly important [2].As shear deformation Mindlin plate has a very wide range of application [3], load identification for Mindlin plate is very meaningful.However, it is difficult to directly measure the loads in a dynamic system, especially the interaction force in complex structures [4].An alternative is to obtain these loads by the inverse force identification, which is an important topic in the identification of Mindlin plate under operating conditions.Therefore, load identification methods are necessary when the direct measurement is unfeasible [5].The main idea of load identification is to use the measured structural response signals to determine the excitation load.However, it is an ill-posed inverse problem [6].Different methods have been developed to make use of the response measurements of a structure to estimate the input excitations [4].Frequency-domain method is one of the most commonly used load identification methods.Structure FRF can be constructed by numerical simulation.It uses a system's FRF and response spectrums to calculate load spectrums and then derives loads in the time domain via the inverse fast Fourier transform [7,8].Through the experiment, some scholars studied load identification using the method of FRF [9].As for an ill-posed inverse problem, regularization methods such as the SVD, Chebyshev-Based Method, cubic B-spline-based method [10], and Tiknohov regularization are often used to transform it into a wellposed problem for great identification accuracy.However, the frequency-domain requires a long stationary time to calculate the inversion of the FRF matrix at each frequency [11][12][13].
It should be noted that the load position or distribution are usually assumed to be known in the load identification procedure [6].Liu analyzed the truncated SVD to overcome ill-posed FRF matrix inversion [14].
Many vibroacoustic studies use finite element method (FEM) model to predict the dynamic behavior of structures.The aim is to reverse the load identification of the model by calculating sources from the measurements of their effects.Current commercial software and computing resources allow engineers to create numerical models that can be applied for more complex structures [15].Zhao [16] presented the sliding loads identification based on the FEM and the strain measurement point selection procedure without information on their initial positions and magnitudes.Bahra [17] adopted the frame frequencies and mode shapes in a FEM model and solved the axial loads until the difference between measured and model frequencies is minimized, thus inferring the member axial loads.Lage [18] proposed a twostep loads identification method: the first step identifies the number of forces and their locations based on the response transmissibility concept; the second step reconstructs the load vector.Berry [19] presented the theoretical developments based on the virtual fields method and processed numerically simulated data to validate the identification algorithm.
However, for many complicated problems, traditional FEM has some disadvantages, such as low efficiency, insufficient accuracy, and slow convergence to correct solutions.Load identification, an ill-posed inverse problem, needs multiresolution analysis to improve the accuracy and reduce the computational work.Due to the property of multiresolution analysis, wavelet provides a natural mechanism for decomposing the solution into a set of coefficients.Wavelet functions can be viewed as interpolating functions, which are similar to those used in signal and image processing.Basu [20] maintained that the finite difference and Ritztype methods have been largely replaced by the FEM, the boundary element method, and the mesh-less method, and in the near future it might be the turn of the waveletbased numerical method.Because wavelet finite element has the advantages of reasonable interpolating functions configuration and higher precision.One can use fewer elements to solve the load identification.HCSWI elements embody a prominent advantage of improving precision by adding the appropriate wavelet functions [21,22].By using the modified one-dimensional HCSWI elements, a multiscale wavelet-based numerical method can solve load identification in rod and Timoshenko beam.Compared with traditional FEM and B-spline wavelet on interval (BSWI) [23][24][25][26][27][28][29] finite element, the method of HCSWI rod and beam elements has the advantage of higher precision.The tensor product of the modified Hermitian wavelets expanded at each coordinate is used to construct two-dimensional Mindlin plate Hermitian wavelet interpolation function.The load identification based on FRF is restricted as the result is inaccurate when the frequency content is close to certain resonance frequency of the structure [30].Hermitian wavelet finite element with high precision can alleviate this problem effectively.
In the present study, an effective new wavelet numerical method based on two-dimensional Mindlin plate Hermitian wavelets interpolation function is proposed to construct the FRF and recover the loads.The paper is organized as follows.Section 2 addresses Hermitian Mindlin plate wavelet element.Section 3 describes load identification of HCSWI.Next, Section 4 details the numerical implementation and analysis.Section 5 presents the experimental verification.Concluding remarks are included in Section 6.

Hermitian Mindlin Plate Wavelet Element
A new Hermitian Mindlin plate wavelet element is constructed by the Hermitian cubic spline wavelets on interval interpolation functions.It can be substituted into finite element functions to solve the Hermitian wavelet stiffness and mass matrix for the purpose of obtaining FRF matrix.
Through the inverse matrix of FRF matrix, one can solve the load identification.The method of SVD is adopted to solve the ill-posed inverse problem.

Hermitian Wavelet Interpolation Functions.
Hermitian scaling functions  1, are shown in Figure 1 and the equations are where Shock and Vibration

3
In order to solve the boundary problem of Hermitian interpolation wavelet, the modified Hermitian scaling functions with all the good performance of Hermitian wavelet [22] are proposed.
The modified scaling functions   1, and the corresponding wavelets  , in the multiresolution approximation space H 1 0 (0, 1) are constructed through one-dimensional multiresolution analysis of HCSWI on the interval [0, 1].A separable wavelet semiorthonormal basis at scale  of  2 (R) is constructed with tensor products of the multiresolution approximation space H 1 0 (0, 1). 1 is the initial scaling space, and   ( = 1, 2, . ..) indicates the wavelet space at different levels.Therefore, the subspace of tensor product is spanned by   = Φ −1 ⊗ Φ −1 .The scaling functions are  1 =  1 ⊗  1 ; the appropriate wavelet functions are As Figure 4 shows, the tensor product of HCSWI scaling functions and wavelet functions are generated by the onedimensional HCSWI at scales  = 1, 2.

Mindlin Plate Theory. Hermitian wavelet interpolation function can be substituted into the potential function for
Mindlin plate to solve the stiffness and mass matrix for load identification.According to Mindlin plate theory, the elemental generalized function of potential energy for Mindlin plate bending problem is [31]        where ] in which   = ℎ 3 /12(1 −  2 ) and   = 5ℎ/12(1 + ), where  is Young's modulus, ℎ is the plate thickness, and  is Poisson's ratio.Substituting , , G D , G C into (8), we can obtain When we use the two-dimensional tensor product of HCSWI scaling functions for scale  to construct element on the elemental solving domain Ω  as shown in Figure 5, each node has three degrees of freedom (DOFs), that is,   ,   ,   ( = 1, 2, . . ., ( + 1) 2 ); the displacement function in the standard elemental solving domain Ω  can be replaced by where physical DOFs of the whole elemental nodes are For (10), first mapping the elemental solving domain Ω  to the standard elemental solving domain Ω  and then substituting ( 11) into (10) and according to variation principle, let Π  = 0; we can obtain the elemental FEM solution where Due to the independent interpolation of   ,   , and , we can change the location of matrix element of the elemental FEM in solving (13) to get  [

Shock and Vibration
where is the consistent mass matrix as follows: where  is the density and  is the area of the Mindlin plate.
In order to satisfy the continuity and compatibility at the interface between neighboring elements and the boundary condition, the stiffness matrix for dynamic analysis should be transformed into physical space.The corresponding elemental DOFs can be solved from wavelet coefficients into unknown field functions.The one-dimensional transformation matrices  1 and  2 can construct a twodimensional transformation matrix   =  1 ⊗  2 .The other procedures can follow those of the traditional FEM [32].
The first-order derivative of HCSWI has the property of continuity, which is suitable for the analysis of plate and shell structure.Therefore, Hermitian wavelet with high precision as a finite element interpolation function can be used in the Mindlin plate.

Integral Solution about Mindlin Plate Stiffness Matrix.
Hermitian interpolation function wavelet can be substituted into the potential function for Mindlin plate to solve the stiffness matrix; the integral solution equation is as follows: The key to realize the multiscale solution of mechanical structures is to compute the integral function.
The integral solution equation is as follows: For the integral Γ 0,1 = ∫ 1 0 Φ(Φ T /), the multiscale constant matrix is where the submatrix of Γ 0,1 can be calculated by Equation ( 23) reveals that only the submatrices C , , ( =  1 , . . .,  −1 ;  =   ) will be calculated while other submatrices remain the same when scale  is lifted to  + 1, thus increasing the calculating efficiency.
The integral solution equation is as follows: For the integral Γ 1,0 = ∫ 1 0 (Φ T /)Φ, the multiscale constant matrix is where the submatrix of Γ 1,0 can be calculated by Equation ( 26) reveals that only the submatrices B , , ( =  1 , . . .,  −1 ;  =   ) can be calculated while other submatrices remain the same when scale  is lifted to +1, thus increasing the calculating efficiency.
Based on the orthogonal characteristic of the wavelet bases with respect to the given inner product, the corresponding multiscale equation will be decoupled across total or partial levels and the wavelet bases suit the nesting approximation.

Load Identification
One can make use of Mindlin plate stiffness matrix and mass matrix and Laplace transform technique to obtain the FRF and then calculate the load.
Equation of motion can be easily discretized and described as follows: where C is the global damping matrix, F is a global excitation load vector, which is associated with time, and ẍ , ẋ , and  stand for the acceleration, velocity, and displacement vectors, respectively.Usually, the damping matrix C can be expressed as a Rayleigh damping if necessary.By using the Laplace transform technique, one can obtain where where H  () is the FRF matrix of the structure Through the modal shape and natural frequency of the structure, the FRF can be obtained.Due to the orthogonality in modal shape, we can obtain where diag is the diagonal matrix, Φ is the modal shape,   is the damping ratio, and   is the natural frequency.According to (35), we can obtain  Substituting ( 36) into (34), we can obtain Assuming  =  ( = √ −1), then the FRF matrix can be obtained: Every element H   () represents the  point excitation and the  point response of the FRF [7].The expression is And load identification equation can be written as In this way, the load identification problem of Mindlin plate can be solved.The main idea of load identification is to use the FRF and the measured structural response signals to determine the excitation load.However, the inverse matrix of the FRF is ill-posed.However, select singular values are the key factors in solving the inverse problem.In this paper, the Lcurve of regularization methods is used to choose the singular    values [33].Therefore, the method of SVD is adopted to solve the problem.SVD is a useful tool that will help one to prove the nature of the deconvolution problem.The SVD of H  is defined as follows:

Shock and Vibration
where U is the matrix of left singular vectors, V is the matrix of right singular vectors, and diagonal elements of matrix Σ are the singular values of H  .The solution of (40) can be rewritten in a term of SVD: The singular values will help one to determine the nature of the problems.The FRF and the measured structural response signals can solve the excitation load.
The equation of absolute error of the identified load is calculated as follows: The root mean square (RMS) error value of the identified load can be given by

Numerical Examples
It is well known that the vibration structures are complex in engineering.In this section, we involve different excitations cases to simulate the load identification in the Mindlin plate.The Mindlin plate dimension is as follows: length 1 m, width 1 m, and thickness 0.01 m.The material is assumed as aluminum, with Young's modulus 70 GPa, Poisson ratio 0.3, and density 2730 kg/m 3 .We fixed the right of the Mindlin plate.A computer is used with an Intel CPU of 3.1 GHz and a memory of 4 G bytes, 64-bit Windows 7 operating system, and Matlab.The responses of different nodes are substituted into (40) to reconstruct the load.Load identification is implemented under an undamped condition.The excitations are used as exact solutions.
Commercial software ANSYS 14.5 is used for modeling, and the shell 43 elements are selected.The mesh of the Mindlin plate is shown in Figure 6.The geometric modeling and material parameters of the Mindlin plate are input into the software ANSYS.Harwell-Boeing method is utilized to extract the stiffness and mass matrix.Three DOFs at each node are retained, that is,   ,   ,   ( = 1, 2, . . ., ( + 1) 2 ).Substituting the stiffness and mass matrix into (40), we can obtain the load identification results of ANSYS.The load identification results of grids 17 × 17 and 35 × 35 are compared with that of HCSWI element.

Load Identification with Single Excitation.
Figure 7 shows that the location of the excitation has been applied at node 1.The sampling frequency studied is 1000 Hz and the time duration of measurement is 1 s.The signal-to-noise ratio (SNR) of the excitation signal is 130 dB.The applied load is used to obtain the displacement of nodes.The displacements of nodes 2, 9, 18, 19, 20, and 145 can be substituted into a HCSWI element and 17 × 17, 35 × 35 ANSYS elements for load reconstruction.Comparison between the precision of HCSWI element and those of ANSYS element verifies the superiority of the former.The excitation  = sin(6 × ) is assumed to act in the Mindlin plate as Figure 7 shows.Figure 8 shows the load identification results at nodes 2, 9, 18, 19, 20, and 145.The high precision is obtained at nodes 2, 19, and 20.When the response nodes such as 9, 145, and 18 are far away from node 1 where the load is applied, the error is great.It can be seen that when the response node is close to node 1, the load identification has a high precision.Obviously, HCSWI element has the highest precision.Figure 11 shows the FRF between the different response nodes and applied node.In this paper we construct the high precision FRF by using HCSWI element.The FRF has a high precision and can calculate the load in engineering practice.

Load Identification with Two Excitations.
Figure 12 shows the location of the excitation signals applied at nodes 1 and 9. Compared with ANSYS results, HCSWI element can accurately identify the applied load, which verifies its superiority.The excitations  1 = sin(6 × ) + sin(12 × ) and  2 = sin(8 × ) are assumed to act in the Mindlin plate as shown in Figure 12.
Figure 13 shows the load identification at nodes 2, 10,     that HCSWI element has higher accuracy than ANSYS elements.
Figure 16 shows the location of the excitation signals applied at nodes 1 and 18. Figure 17 shows the load identification at nodes 2, 5, 10, and 17.Because node 2 is close to the location of applied load  1 and away from that of  2 , the response at node 2 can accurately identify  1 but cannot identify  2 .For the same reason, the response node 17 can accurately identify  2 but cannot identify  1 .Figures 18 and  19 show that the accuracy of HCSWI element is higher than ANSYS elements.

Experimental Verification
The Mindlin plate structure for load identification is set up as shown in Figure 20.The plate has been used by Qiao et al. for impact and harmonic force identification [34].The Mindlin plate dimension is as follows: length 0.5 m, width 0.5 m, and height 0.005 m.The material is number 45 steel with Young's modulus 206 GPa, Poisson ratio 0.3, and density 7917 kg/m 3 .The left of the Mindlin plate is fixed.In the impact experiment, the excitation which is marked as  in Figure 21(a  response, the forces can be identified by HCSWI and ANSYS elements. A modal impact hammer (PCB 086C01 with sensitivity 12.37 mV/N) is used to impose impact forces.In the impact experiment, the Mindlin plate is impacted many times.Accelerometers sensor and force channels are synchronously sampled by LMSSCADASIII data acquisition system with a sampling frequency of 4096 Hz and spectral line of 1025.At each excitation, the FRF is measured.After the inverse fast Fourier transform, the data length of the discrete impulse response function is 2050.In this paper, we use SVD to solve the ill-posed problem since it is an efficient algorithm.
The structure is assumed to exhibit Rayleigh damping as follows: where  1 and  2 are the damping coefficients.During the impact experiment, impact excitations are continuously recorded.Each case is impacted using a modal impact hammer to apply and measure the load, while the acceleration response is measured using an accelerometer sensor described previously.Figure 21 shows the original impact loads acting on the location of node 9 and the acceleration response at node 8.It can be seen from Figure 21 that impact loads have a typical impulsive profile and the  response of each impact case displays a rapid attenuation tendency.
One can correct the elastic moduli; those are deduced from the characteristic value of the Mindlin plate [35][36][37][38] The experiment measures the acceleration frequency response functions.The following equation is the acceleration frequency response functions solved by FEM.The accuracy of HCSWI can be verified by comparing the FRF of HCSWI and ANSYS elements with the experiment data: Figure 22 shows the acceleration FRF at nodes 8, 9, 8, and 26 after the model correction.As shown in Figure 22, the result of HCSWI is the closest to the experimental acceleration FRF.With the increase of ANSYS grid, the frequency responses of ANSYS are close to the experimental results.
Figure 23 shows the load identification results at nodes 9 and 26.When the excitation node such as 26 is far away from node 8 where the accelerometer sensor is mounted, the error is great.It can be seen that when the excitation node is close to node 8, the load identification has a high precision.Obviously, HCSWI element has the highest precision.
Figures 24 and 25 show the absolute error and the RMS error value by using HCSWI and ANSYS 17 × 17 and 35 × 35.It can be seen that HCSWI element error line is lower than result of ANSYS.And ANSYS 35 × 35 error is lower than that of 17 × 17.The accuracy of ANSYS is gradually improved with the increasing number of grids.The result verifies the superiority of HCSWI element.

Conclusion
The present study constructs a Hermitian Mindlin plate wavelet element.It is substituted into finite element functions to solve the FRF and then derive the loads.Numerical calculations are carried out with different excitation cases in the Mindlin plate.The SVD is adopted to solve ill-posed inverse problem of the load identification.The results of ANSYS shell 43 elements with grids 17 × 17 and 35 × 35 are compared with that of HCSWI Mindlin plate element.It shows that the accuracy of ANSYS gradually improves with the increasing number of grids.Compared with ANSYS, HCSWI Mindlin plate element can accurately identify the applied load.The load identification has desirable results when the response node is close to the node where the load is applied.The accuracy of HCSWI can be verified by comparing the FRF of HCSWI and ANSYS elements with the experiment data.The FRFs and response spectrums can calculate the loads.The experiment proves that the load identification of HCSWI Mindlin plate is effective and precise.But the frequency-domain requires a long stationary time to calculate the inversion of the FRF matrix at each frequency.The load identification based on FRF is restricted as the result is inaccurate when the frequency content is close to certain resonance frequency of the structure.Hermitian wavelet finite element with high precision can alleviate this problem but cannot fundamentally solve the problem.

Figure 3 :
Figure 3: Wavelets in wavelet space : (a) wavelet in wavelet space  1 and (b) wavelet in wavelet space  2 .

Figure 7 :
Figure 7: The location of  in the Mindlin plate.

Figure 10 :
Figure 10: The root mean square error value with the excitation  applied at node 1.

Figure 14 :
Figure 14: The absolute error with excitation  1 and excitation  2 applied at nodes 1 and 9.

Figure 15 :
Figure 15: The root mean square error value with excitation  1 and excitation  2 applied at nodes 1 and 9.

Figure 16 :
Figure 16: The location of  1 and  2 in the Mindlin plate.

Figure 18 :
Figure 18: The absolute error with excitation  1 and excitation  2 applied at nodes 1 and 18.

Figures 9 and 10
Figure8shows the load identification results at nodes 2, 9, 18, 19, 20, and 145.The high precision is obtained at nodes 2, 19, and 20.When the response nodes such as 9, 145, and 18 are far away from node 1 where the load is applied, the error is great.It can be seen that when the response node is close to node 1, the load identification has a high precision.Obviously, HCSWI element has the highest precision.Figures 9 and 10 show the absolute error and the RMS error value by using HCSWI and ANSYS 17 × 17 and 35 × 35.It can be seen that HCSWI element error line is almost horizontal.And ANSYS 35 × 35 error is lower than that of 17 × 17.The accuracy of ANSYS is gradually improved with the increasing number of grids.The result verifies the superiority of HCSWI element.Figure11shows the FRF between the different response nodes and applied node.In this paper we construct the high

Figure 19 :
Figure 19: The root mean square error value with excitation  1 and excitation  2 applied at nodes 1 and 18.

Figure 20 :
Figure 20: The Mindlin plate structure applied for load identification.

Figure 21 :
Figure 21: The measurement data from force transducer and accelerometer: (a) impact loads and (b) acceleration responses.

Figure 22 :
Figure 22: Acceleration FRF with the excitation at nodes 9 and 26 and the response at node 8.

Figure 24 :
Figure 24: The absolute error with excitation  applied at nodes 9 and 26.
Identified load  at node 26

Figure 25 :
Figure 25: The root mean square error value with excitation  applied at nodes 9 and 26.

2 M
. The corresponding elastic modulus value    ( = 1, 2, . ..) is solved for each frequency, which makes the experiments frequency values equal to calculating frequency values of FEM.The equation of the corrected elastic modulus can be given by 2 Identified load f 1 at node 17 (d2) Identified load f 2 at node 17 Figure 17: Load identification at different nodes with excitation  1 and excitation  2 applied at nodes 1 and 18.
18, and 19.Because node 2 or 19 is close to the location of applied load  1 and away from that of  2 , the response node 2 or 19 can accurately identify  1 but cannot identify  2 .For the same reason, the response at node 10 can accurately identify  2 but cannot identify  1 .Figures 14 and 15 show (d) Identified load  1 at node 19 ) is continuously impacted.The Mindlin plate is divided into 5 × 5 elements (6 × 6 nodes).The excitation locations are at nodes 9 and 26.The accelerometer sensor is mounted at node 8.For the purposes of load identification validations, the force signals are recorded by the force transducers inserted in the impact hammer.The response signals are recorded by accelerometer sensor at the selected node 8. Therefore, the forces acting on the location at nodes 9 and 26 can be constructed from the response of the selected node 8.With the acceleration Identified load  at node 26 Figure 23: Load identification at node 8 with excitation  applied at nodes 9 and 26.