# Adaptive Reconstruction of a Dynamic Force Using Multiscale Wavelet Shape Functions

Authors: Wen-Yu He; Yang Wang; Songye Zhu
Journal: Shock and Vibration
 (2018)

---

## Abstract

The shape function-based method is one of the very promising time-domain methods for dynamic force reconstruction, because it can significantly reduce the number of unknowns and shorten the reconstruction time. However, it is challenging to determine the optimum time unit length that can balance the tradeoff between reconstruction accuracy and efficiency in advance. To address this challenge, this paper develops an adaptive dynamic force reconstruction method based on multiscale wavelet shape functions and time-domain deconvolution. A concentrated dynamic force is discretized into units in time domain and the local force in each unit is approximated by wavelet scale functions at an initial scale. Subsequently, the whole response matrix is formulated by assembling the responses induced by the wavelet shape function forces of all time units which are calculated by the structural finite element model (FEM). Then, the wavelet shape function-based force-response equation is established for force reconstruction. Finally, the scale of the force-response equation is lifted by refining the wavelet shape function with high-scale wavelets and dynamic responses with more point data to improve the reconstruction accuracy gradually. Numerical examples of different structural types are analyzed to verify the feasibility and effectiveness of the proposed method.


---

## Body

## 1. Introduction

Forward and inverse dynamic analyses are two typical types of structural dynamic problems, containing three basic components, namely, excitations, structures, and responses. The forward analysis refers to structural response calculation with the knowledge of excitation and structure parameters. The inverse analysis can be classified into two types: structural parameters identification using the known excitation and response and force reconstruction using the known structure parameters and response. Dynamic force reconstruction methods based on structural dynamic response have attracted great interest in the field of structural health monitoring (SHM), as it is often difficult or impractical to measure a dynamic force directly. Dynamic force reconstruction methods can generally be categorized into two groups: frequency-domain [1–4] and time-domain [5–8] methods. The frequency- and time-domain methods estimate dynamic forces by establishing the relationship between dynamic forces and structural responses based on frequency response functions [9, 10] and impulse response functions [11, 12], respectively. Compared to the frequency-domain methods, the time-domain methods have received increasing attention in the past years because of their distinct physical meaning and relatively higher accuracy [13].Most time-domain methods are based on discretization in the time domain. The sampling time interval size is one main factor that affects the reconstruction accuracy and efficiency. In general, a smaller size of sampling time interval can achieve higher reconstruction accuracy. However, a very small size would result in too many undetermined coefficients, which not only increases the computation cost but also tends to make the inverse problem ill-conditioned [13]. Therefore, it is challenging to set an optimal sampling time interval to balance the tradeoff between reconstruction accuracy and efficiency. Approaches based on various basis functions were proposed to address this challenge, in which unknown dynamic forces were approximated by basis functions, such as Gaussian basis functions [14], B-spline functions [15, 16], triangle functions [17], exponential function [18], and Daubechies wavelet [19], and subsequently they were reconstructed by identifying the coefficients in these basis functions. These approaches could significantly reduce the number of unknowns (considerably less than that of data points) and shorten the identification time.The majority of the above-mentioned basis function-based methods are whole-domain methods, which are similar to the Rize method in structural analysis. Whole-domain methods are not efficient when actual dynamic forces are complex. Therefore, Liu et al. [20] introduced the concept of “shape function” from finite element methods (FEM) into the field of force identification. A dynamic force was discretized into time units in the time domain and the local force in each time unit was approximated by shape functions. Numerical examples indicated that the shape function-based method could provide better results than traditional methods. The reconstruction accuracy of shape function-based methods depends on the approximation ability of shape functions, which is mainly determined by the length of the discrete time unit. The shape function-based methods transform the task of determining the optimal sampling time interval in traditional time-domain force identification methods into determining the optimal time unit length. In theory, shortening the time unit length can enhance the identification accuracy at the cost of an increased computation amount. However, determining the optimum time unit length in advance that can achieve the tradeoff between reconstruction accuracy and efficiency is not easy.This paper aims to develop an adaptive method for dynamic force reconstruction based on multiscale wavelet shape functions and time-domain deconvolution. The time history of a concentrated dynamic force is discretized into time units, and wavelet scale functions at an initial scale are employed to approximate the local force within each unit. Then, the responses induced by the forces expressed using wavelet shape functions are computed by structural FEM, and the whole response matrix is formulated by assembling the responses induced by the forces expressed using different wavelet shape functions in all time units. This process is similar to assembling stiffness matrices in FEM. Subsequently, the wavelet shape function-based force-response equation is established for force reconstruction. The scale of the force-response equation is lifted by refining the wavelet shape functions with high-scale wavelets and dynamic response with more point data when necessary; consequently, the reconstruction accuracy is improved gradually. Finally, numerical examples of a simply supported beam and a two-story frame are analyzed to verify the feasibility and effectiveness of the proposed method. With this adaptive method, the challenging task in determining the time unit length can be overcome.


## 2. Multiscale Wavelet Shape Functions

### 2.1. Multiresolution Wavelet Analysis

Known as a “numerical microscope,” wavelet was firstly applied in the field of signal and image processing, and later its use was extended to numerical computation in engineering [21]. Multiresolution analysis is one of the most important characteristics of wavelets [22]. A multiresolution analysis R of L2 is a sequence of closed subspaces R=Vj⊂L2∣j∈Z, such that [22](1)
Vj⊂Vj+1;(2)
⋃j∈JVj is dense in L2;(3)
for eachj∈J, Vj has a Riesz basis given by scaling functions φj,k∣k∈Kj, where j is the level of resolution, J is an integer index set associated with resolution levels, Kj is some index set associated with the scaling functions of level j, and Vj denotes the approximation spaces of level j. For each Vj, a complement of Vj exists in Vj+1, namely, Wj. Let spaces Wj be spanned by wavelets ψj,mx for every m∈Mj, Mj=Kj+1∖Kj, where M(j) is the difference set of Kj+1 and K(j). Furthermore, let l∈K(j+1) be the index at level j+1.Multiresolution analysis enables wavelets to represent a finite energy function in a dynamic multiscale manner with different precision levels. For example, functionf(x) can be approximated in V0 (low-scale space) as(1)fx≈f0x=∑kc-0,kϕ0,k,where ϕ0,k is the scaling function at Scale 0 and c-0,k represents the corresponding wavelet coefficients. (2)c-0,k=fx,φ0,kx.Approximation accuracy can be improved by adding high-scale wavelet terms in wavelet spaceW0, which is spanned by the wavelet functions ψ0,l. Therefore, the approximation in space V1 (high-scale space) is(3)fx≈f1x=∑kc-0,kϕ0,k+∑lc0,lψ0,l,where ψ0,l is the wavelet function at Scale 0 and c0,l is the corresponding wavelet coefficients in space W0.(4)c0,l=fx,ψ0,lx.By further increasing the approximation order, the wavelet representation of the function approaches the exact function whenj→∞.(5)fx≈fjx=∑kc-0,kϕ0,k+∑j∑lcj,lψj,l,fx=fj→∞x.


### 2.2. The Refinement of Lagrange Wavelet Function

Various wavelets, including spline wavelet [23], Daubechies wavelet [24], Hermite wavelet [25], Lagrange wavelet [26], and trigonometric wavelet [27], have been employed in numeral computation. In this study, the second-generation Lagrange wavelet [26] with convenient computational characteristics is adopted to fulfill the adaptive dynamic force reconstruction.When defined within the interval[0,1], the scaling functions (Scale 0) of Lagrange wavelets are(6a)ϕ0,1=1-xx∈0,10otherwise(6b)ϕ0,2=xx∈0,10otherwise.The corresponding wavelet function ψj,l(x) at scale j is(7)ψj,lx=2j+1x-2l+2x∈l-12j,l-0.52j-2j+1x+2lx∈l-0.52j,l2j0otherwise.The scale functions (V0) and wavelet functions (W0, W1, and W2) are shown in Figure 1.Figure 1
Lagrange wavelet shape functions.


## 2.1. Multiresolution Wavelet Analysis

Known as a “numerical microscope,” wavelet was firstly applied in the field of signal and image processing, and later its use was extended to numerical computation in engineering [21]. Multiresolution analysis is one of the most important characteristics of wavelets [22]. A multiresolution analysis R of L2 is a sequence of closed subspaces R=Vj⊂L2∣j∈Z, such that [22](1)
Vj⊂Vj+1;(2)
⋃j∈JVj is dense in L2;(3)
for eachj∈J, Vj has a Riesz basis given by scaling functions φj,k∣k∈Kj, where j is the level of resolution, J is an integer index set associated with resolution levels, Kj is some index set associated with the scaling functions of level j, and Vj denotes the approximation spaces of level j. For each Vj, a complement of Vj exists in Vj+1, namely, Wj. Let spaces Wj be spanned by wavelets ψj,mx for every m∈Mj, Mj=Kj+1∖Kj, where M(j) is the difference set of Kj+1 and K(j). Furthermore, let l∈K(j+1) be the index at level j+1.Multiresolution analysis enables wavelets to represent a finite energy function in a dynamic multiscale manner with different precision levels. For example, functionf(x) can be approximated in V0 (low-scale space) as(1)fx≈f0x=∑kc-0,kϕ0,k,where ϕ0,k is the scaling function at Scale 0 and c-0,k represents the corresponding wavelet coefficients. (2)c-0,k=fx,φ0,kx.Approximation accuracy can be improved by adding high-scale wavelet terms in wavelet spaceW0, which is spanned by the wavelet functions ψ0,l. Therefore, the approximation in space V1 (high-scale space) is(3)fx≈f1x=∑kc-0,kϕ0,k+∑lc0,lψ0,l,where ψ0,l is the wavelet function at Scale 0 and c0,l is the corresponding wavelet coefficients in space W0.(4)c0,l=fx,ψ0,lx.By further increasing the approximation order, the wavelet representation of the function approaches the exact function whenj→∞.(5)fx≈fjx=∑kc-0,kϕ0,k+∑j∑lcj,lψj,l,fx=fj→∞x.


## 2.2. The Refinement of Lagrange Wavelet Function

Various wavelets, including spline wavelet [23], Daubechies wavelet [24], Hermite wavelet [25], Lagrange wavelet [26], and trigonometric wavelet [27], have been employed in numeral computation. In this study, the second-generation Lagrange wavelet [26] with convenient computational characteristics is adopted to fulfill the adaptive dynamic force reconstruction.When defined within the interval[0,1], the scaling functions (Scale 0) of Lagrange wavelets are(6a)ϕ0,1=1-xx∈0,10otherwise(6b)ϕ0,2=xx∈0,10otherwise.The corresponding wavelet function ψj,l(x) at scale j is(7)ψj,lx=2j+1x-2l+2x∈l-12j,l-0.52j-2j+1x+2lx∈l-0.52j,l2j0otherwise.The scale functions (V0) and wavelet functions (W0, W1, and W2) are shown in Figure 1.Figure 1
Lagrange wavelet shape functions.


## 3. Adaptive Dynamic Force Reconstruction Method

Wavelet functions or scale functions are employed as elemental shape functions in the wavelet FEM to analyze structures in an adaptive manner according to actual needs. When the wavelet FEM is employed to analyze structure problems, initial solutions with relatively low accuracy are obtained in low-scale models, and more accurate results can be acquired by lifting the model scale [21]. With this adaptive analysis, a tradeoff between the computation accuracy and computation efficiency can be achieved. In this section, this adaptive strategy will be further extended for dynamic force reconstruction to overcome the challenge in determining time unit length.


### 3.1. Time-Domain Deconvolution Method

Time-domain deconvolution is one of the most widely used methods for structural dynamic force reconstruction [28]. By assuming zero initial conditions, the structural response yt can be expressed as a convolution integral of an external dynamic force ft and the corresponding impulse response function h(t-τ) [20]:(8)yt=∫0tht-τfτdτ.In a real dynamic test, structural dynamic responses are usually discrete. Suppose there arer sampling points in total. Equation (8) can be transformed into a discrete type:(9)yt1yt2⋮ytr=ht10…0ht2ht1…0⋮⋮⋱⋮htrhtr-1…ht1ft1ft2⋮ftrΔt,where Δt is the discrete time interval. With this equation, the time history of the dynamic force ft can be reconstructed by the measured structural response yt and the impulse response function ht-τ generally obtained from the structural FEM.The size of the sampling time intervalΔt affects both the reconstruction accuracy and the efficiency. In general, a small size of sampling time interval is required to achieve high reconstruction accuracy. However, a very small size would result in too many undetermined coefficients, which not only increase the computation cost but also tend to make the inverse problem more ill-conditioned. In view of this, Liu et al. [20] introduced the concept of “shape function”: the time history of a dynamic force was discretized into local time units and the corresponding local force was approximated by shape functions. Then, the whole response matrix was formulated by assembling the responses induced by the loads corresponding to different shape functions in all time units. Thus, the dynamic force identification model is established through the force-response equation.Suppose a dynamic force in the whole time domain is discretized intop time units, and each unit contains q data points, and the total sampling points are r=pq-1+1. Then, (9) can be further changed for shape function-based force reconstruction as follows:(10)Y=KD, where Y=yt1yt2⋯ytrT and D=d1d2⋯drT are the response vector at r sampling points and the shape function coefficients, respectively. K is the assembly of Km, where the assembling process is similar to the method of stiffness matrices assembling in FEM [20].(11)K=∑m=1rKmKm=0⋯0HNmt10⋯00⋯0HNmt2HNmt1⋯0⋮⋮⋮⋱⋮0⋯0HNmtr-m+1HNmtr-m⋯HNmt1⋮⋮⋮⋱⋮0⋯0HNmtr-1HNmtr-2⋯HNmti-10⋯0HNmtrHNmtr-1⋯HNmti00⋮HNmt1⋮HNmti-1HNmtimrowsazzzzzaa︸m-1columnsm=1,2,…,r,where HNm is the dynamic response induced by the mth shape function Nm.The coefficient matrix in (10) is generally ill-posed and little perturbation may cause tremendous variation in solutions. To obtain a stable and meaningful solution, it is often transformed into a Tikhonov regularization optimization problem [29]:(12)min⁡Y-KD22+λD22,where λ is the regularization factor. The regularization factor weighs the norm of residual and the norm of solution. It is determined by the L-curve method in this paper [30]. Generally, the relationship of the solution norm versus the residual norm is plotted for different λ, which looks like an “L” shape in the log-log scale, and the optimal regularization parameter is selected at the corner of the “L” shape [30].The dynamic force can be estimated by the following equation:(13)F=DN, where F and N=N1N2⋯NrT are the dynamic force to be identified and the shape functions, respectively.


### 3.2. Wavelet Shape Function-Based Method

The shape function-based dynamic force reconstruction methods can considerably reduce the coefficient matrix size and enhance the computation efficiency. The accuracy of the shape function-based reconstruction method depends on the approximation ability of the shape functions, which is mainly determined by the discrete time unit length. However, it is not easy to determine the optimum time unit length that can balance the tradeoff between reconstruction accuracy and efficiency in advance. This section aims to develop an adaptive method for the reconstruction of dynamic forces based on multiresolution Lagrange wavelet shape functions and time-domain deconvolution.Using the scaling and wavelet functions of Lagrange waveletsΦj=Φ0T,Ψ0T,Ψ1T,…,Ψj-1T at Scale j (described in Section 2.2) as shape functions (N), the unknown dynamic force at each unit can be approximated as(14)y=∑kc-0,kϕ0,k+∑m=0j-1∑lcm,lψm,l=Φ0c-+Ψ0c0+Ψ1c1+⋯+Ψj-1cj-1=Φjqj,where Φ0 represents the scaling functions at Scale 0, Ψj represents the wavelet functions at Scale j, and qj=c-c0c1⋯cj-1T is the undetermined vector of wavelet coefficients.Substitute (14) into (10) and (12), and then the wavelet coefficients can be obtained and the force can be reconstructed according to (13). The corresponding solution is denoted as Scale j solution. Obviously, higher scale solutions would be more accurate than lower scale ones.Based on the multiresolution Lagrange wavelet shape functions and time-domain deconvolution, an adaptive method is developed for dynamic force reconstruction. Figure2 presents the flowchart of this adaptive force reconstruction scheme. The detailed procedure is described as follows.Figure 2
Diagram of the adaptive force reconstruction scheme.Step 1.
Install a sensor on the concerned structure and measure its dynamic response (e.g., displacement, velocity, acceleration, or strain), and discretize the full dynamic response time history intop time units with an equal length.Step 2.
Use the structural FEM to calculate the response corresponding to wavelet shape functions at Scale 0.Step 3.
Establish the Scale 0 force-response equation, and obtain the corresponding solution.Step 4.
Refine the wavelet shape functions by adding high-scale wavelets, and refine dynamic responses by adding more data points. Establish the Scale 1 force-response equation and obtain the corresponding solution. Consequently, the reconstruction accuracy can be enhanced.Step 5.
Check the convergence of the solutions, and stop if the difference is smaller than a prescribed threshold. Otherwise, repeat Step4.During the scale lifting process for the force-response equation, the submatrices/vectors of the current scale can be retained, and only a few rows and columns need to be added. The new results can be quickly obtained via iteration, with the initial values equal to the results at the previous scale. With this adaptive strategy, the force reconstruction accuracy can be enhanced by lifting the scale of the force-response equation gradually when necessary, and thus determining the time unit length in advance, which is often a challenging task, is not required.


## 3.1. Time-Domain Deconvolution Method

Time-domain deconvolution is one of the most widely used methods for structural dynamic force reconstruction [28]. By assuming zero initial conditions, the structural response yt can be expressed as a convolution integral of an external dynamic force ft and the corresponding impulse response function h(t-τ) [20]:(8)yt=∫0tht-τfτdτ.In a real dynamic test, structural dynamic responses are usually discrete. Suppose there arer sampling points in total. Equation (8) can be transformed into a discrete type:(9)yt1yt2⋮ytr=ht10…0ht2ht1…0⋮⋮⋱⋮htrhtr-1…ht1ft1ft2⋮ftrΔt,where Δt is the discrete time interval. With this equation, the time history of the dynamic force ft can be reconstructed by the measured structural response yt and the impulse response function ht-τ generally obtained from the structural FEM.The size of the sampling time intervalΔt affects both the reconstruction accuracy and the efficiency. In general, a small size of sampling time interval is required to achieve high reconstruction accuracy. However, a very small size would result in too many undetermined coefficients, which not only increase the computation cost but also tend to make the inverse problem more ill-conditioned. In view of this, Liu et al. [20] introduced the concept of “shape function”: the time history of a dynamic force was discretized into local time units and the corresponding local force was approximated by shape functions. Then, the whole response matrix was formulated by assembling the responses induced by the loads corresponding to different shape functions in all time units. Thus, the dynamic force identification model is established through the force-response equation.Suppose a dynamic force in the whole time domain is discretized intop time units, and each unit contains q data points, and the total sampling points are r=pq-1+1. Then, (9) can be further changed for shape function-based force reconstruction as follows:(10)Y=KD, where Y=yt1yt2⋯ytrT and D=d1d2⋯drT are the response vector at r sampling points and the shape function coefficients, respectively. K is the assembly of Km, where the assembling process is similar to the method of stiffness matrices assembling in FEM [20].(11)K=∑m=1rKmKm=0⋯0HNmt10⋯00⋯0HNmt2HNmt1⋯0⋮⋮⋮⋱⋮0⋯0HNmtr-m+1HNmtr-m⋯HNmt1⋮⋮⋮⋱⋮0⋯0HNmtr-1HNmtr-2⋯HNmti-10⋯0HNmtrHNmtr-1⋯HNmti00⋮HNmt1⋮HNmti-1HNmtimrowsazzzzzaa︸m-1columnsm=1,2,…,r,where HNm is the dynamic response induced by the mth shape function Nm.The coefficient matrix in (10) is generally ill-posed and little perturbation may cause tremendous variation in solutions. To obtain a stable and meaningful solution, it is often transformed into a Tikhonov regularization optimization problem [29]:(12)min⁡Y-KD22+λD22,where λ is the regularization factor. The regularization factor weighs the norm of residual and the norm of solution. It is determined by the L-curve method in this paper [30]. Generally, the relationship of the solution norm versus the residual norm is plotted for different λ, which looks like an “L” shape in the log-log scale, and the optimal regularization parameter is selected at the corner of the “L” shape [30].The dynamic force can be estimated by the following equation:(13)F=DN, where F and N=N1N2⋯NrT are the dynamic force to be identified and the shape functions, respectively.


## 3.2. Wavelet Shape Function-Based Method

The shape function-based dynamic force reconstruction methods can considerably reduce the coefficient matrix size and enhance the computation efficiency. The accuracy of the shape function-based reconstruction method depends on the approximation ability of the shape functions, which is mainly determined by the discrete time unit length. However, it is not easy to determine the optimum time unit length that can balance the tradeoff between reconstruction accuracy and efficiency in advance. This section aims to develop an adaptive method for the reconstruction of dynamic forces based on multiresolution Lagrange wavelet shape functions and time-domain deconvolution.Using the scaling and wavelet functions of Lagrange waveletsΦj=Φ0T,Ψ0T,Ψ1T,…,Ψj-1T at Scale j (described in Section 2.2) as shape functions (N), the unknown dynamic force at each unit can be approximated as(14)y=∑kc-0,kϕ0,k+∑m=0j-1∑lcm,lψm,l=Φ0c-+Ψ0c0+Ψ1c1+⋯+Ψj-1cj-1=Φjqj,where Φ0 represents the scaling functions at Scale 0, Ψj represents the wavelet functions at Scale j, and qj=c-c0c1⋯cj-1T is the undetermined vector of wavelet coefficients.Substitute (14) into (10) and (12), and then the wavelet coefficients can be obtained and the force can be reconstructed according to (13). The corresponding solution is denoted as Scale j solution. Obviously, higher scale solutions would be more accurate than lower scale ones.Based on the multiresolution Lagrange wavelet shape functions and time-domain deconvolution, an adaptive method is developed for dynamic force reconstruction. Figure2 presents the flowchart of this adaptive force reconstruction scheme. The detailed procedure is described as follows.Figure 2
Diagram of the adaptive force reconstruction scheme.Step 1.
Install a sensor on the concerned structure and measure its dynamic response (e.g., displacement, velocity, acceleration, or strain), and discretize the full dynamic response time history intop time units with an equal length.Step 2.
Use the structural FEM to calculate the response corresponding to wavelet shape functions at Scale 0.Step 3.
Establish the Scale 0 force-response equation, and obtain the corresponding solution.Step 4.
Refine the wavelet shape functions by adding high-scale wavelets, and refine dynamic responses by adding more data points. Establish the Scale 1 force-response equation and obtain the corresponding solution. Consequently, the reconstruction accuracy can be enhanced.Step 5.
Check the convergence of the solutions, and stop if the difference is smaller than a prescribed threshold. Otherwise, repeat Step4.During the scale lifting process for the force-response equation, the submatrices/vectors of the current scale can be retained, and only a few rows and columns need to be added. The new results can be quickly obtained via iteration, with the initial values equal to the results at the previous scale. With this adaptive strategy, the force reconstruction accuracy can be enhanced by lifting the scale of the force-response equation gradually when necessary, and thus determining the time unit length in advance, which is often a challenging task, is not required.


## 4. Numerical Study

Numerical examples of a simply supported beam and a double-story frame are analyzed to verify the feasibility and effectiveness of the proposed adaptive dynamic force reconstruction method. The dynamic responses and wavelet shape function responses are calculated by structural FEM in which Rayleigh damping is considered. The Newmark-β method is employed for the dynamic analysis with a time increment of 0.001 s.To simulate a practical measurement condition environment, 5% noise is introduced by adding zero-mean Gaussian noise to the accurate displacement responses [31]:(15)u=ucal+5%×Tnoise×σucal,where u and ucal are the displacement responses with and without noise, respectively. Tnoise is a vector of independent random variables following a standard normal distribution, and σucal is the standard deviation of the response.A simple example is employed to illustrate the adaptive force reconstruction procedure (Figure3). The full dynamic response time history is divided into 3 time units with equal length. In Scale 0, the responses (0, t01, t02, and t03) and Scale 0 Lagrange wavelet scaling functions (ϕ0,1 and ϕ0,2) are adopted to establish the force-response equation (denoted as Scale 0 force-response equation) as (10). The size of the coefficient matrix is 4×4. Then, the Scale 0 force-response equation is changed into a Tikhonov regularization optimization problem as (13). By solving this regularization optimization problem, the Scale 0 solutions can be obtained. Subsequently, the Scale 0 equation is lifted into Scale 1 force-response equation by adding Scale 0 Lagrange wavelet functions (ψ0,1) and responses at time t1,i (i=1,2,3). The size of the coefficient matrix in Scale 1 force-response equation is extended to 7×7 accordingly. Further, the responses at time t2,i (i=1,2,…,6) and the Scale 1 Lagrange wavelet functions (ψ1,1 and ψ1,2) and the response at time t3,i (i=1,2,…,12) and the Scale 2 Lagrange wavelet functions (ψ2,1, ψ2,2, ψ2,3, and ψ2,4) are added to formulate the Scale 2 and Scale 3 force-response equations gradually. The sizes of the coefficient matrix in Scale 2 and Scale 3 force-response equation are 13×13 and 25×25, respectively. The size of the coefficient matrix and the reconstruction accuracy would be enhanced with the scale number.Figure 3
Example of the adaptive force reconstruction process.Figure4 shows a simply supported beam with a length of L=16m, Young’s modulus of 3.3×1011 N/m2, density of 2.5×103 Kg/m3, and cross-sectional area of A=1×4 m2. The scenarios, dynamic force time histories, applied locations, and measurement locations are listed in Table 1. The whole beam is divided into 16 beam elements equally, and the FEM is established accordingly.Table 1
Scenarios considered in the numerical simulations.
Structure type
Scenarios
Force
Measurement location
Time history (N)
Location
Simply supported beam
Case B1
sin
⁡
(
π
t
)
+
cos
⁡
(
2
π
t
)
+
sin
⁡
(
4
π
t
)
BF1
BP1
Case B2
sin
⁡
(
π
t
)
+
cos
⁡
(
2
π
t
)
+
cos
⁡
(
3
π
t
)
+
cos
⁡
(
4
π
t
)
BF2
BP2
Double-story frame
Case F1
sin
⁡
(
π
t
)
+
sin
⁡
(
2
π
t
)
+
cos
⁡
(
2
π
t
)
+
sin
⁡
(
4
π
t
)
FF1
FP1
Case F2
cos
⁡
(
π
t
)
+
sin
⁡
(
2
π
t
)
+
cos
⁡
(
2
π
t
)
+
sin
⁡
(
4
π
t
)
FF2
FP2Figure 4
The simulated simply supported beam.For Case B1, the beam is excited by applying the following dynamic force at the middle span (BF1) (Figure4):(16)FB1t=sin⁡πt+cos⁡2πt+sin⁡4πtKN0≤t≤10.The displacement response at Point BP1 in they-direction is used to reconstruct the dynamic force. The total response time history is discrete in 50 time units with equal length. Following the procedures presented above, the dynamic force is reconstructed adaptively. The reconstruction parameters are listed in Table 2. The solutions and relative errors with respect to the maximal real force obtained at different scales are plotted in Figures 5 and 6, respectively. It can be seen that the reconstruction accuracy increases with the scale. For example, the maximal relative errors in Scale 0 to Scale 3 are about 30%, 10%, 2.5%, and 0.5% (except several starting points), respectively.Table 2
Reconstruction parameters in numerical simulations.
Item
Scale 0
Scale 0
Scale 0
Scale 3
Time unit number
50
50
50
50
Coefficient matrix size
51
×
51
101
×
101
201
×
201
401
×
401Figure 5
The reconstructed dynamic force in Case B1.
(a)
Scale 0 and Scale 1
(b)
Scale 2 and Scale 3Figure 6
The reconstructed error in Case B1.
(a)
Scale 0 and Scale 1
(b)
Scale 2 and Scale 3Considering the fact that the Scale 3 solution has sufficient precision, only three scales are calculated in the following cases. If higher precision is required, a higher scale force-response equation could be formulated through the lifting process.For Case B2, the dynamic force (BF2) with the following expression is applied at the section which is 4  m from the left end (Figure4):(17)FB1t=sin⁡πt+cos⁡2πt+cos⁡3πt+cos⁡4πtKN0≤t≤10.The displacement response at Point BP2 is employed to perform the proposed adaptive dynamic force reconstruction method. The reconstruction parameters are listed in Table2, and the solutions and relative errors of different scales are shown in Figures 7 and 8, respectively. Similar conclusions can be drawn in Case B2 as well.Figure 7
The reconstructed dynamic force in Case B2.Figure 8
The reconstructed error in Case B2.Furthermore, a single-bay two-story frame is simulated, as shown in Figure9. Young’s modulus and the density of the material are 3.3×1011 N/m2 and 2.5×103 Kg/m3, respectively. The cross sections are A=1×1 m2 for all the beams and columns. Each beam and column is divided into 4 beam elements equally in developing the FEM. The scenarios, force time histories, applied locations, and measurement locations are listed in Table 1, and the reconstruction parameters are listed in Table 2. The solutions of different scales in Case F1 and Case F2 are shown in Figures 10 and 12, respectively. The corresponding relative errors are compared in Figures 11 and 13, respectively.Figure 9
The simulated double-story frame.Figure 10
The reconstructed dynamic force in Case F1.Figure 11
The reconstructed error in Case F1.Figure 12
The reconstructed dynamic force in Case F2.Figure 13
The reconstructed error in Case F2.The simulation results of the simply supported beam and the two-story frame indicate that the reconstruction accuracy could be enhanced by lifting the scale of the wavelet shape function-based force-response equation gradually. This adaptive strategy can avoid determining the time unit length in advance, which is often a challenge in practice. The reconstruction scale that affects both the accuracy and the efficiency can be determined adaptively according to the actual needs.


## 5. Conclusion

As usually dynamic forces can hardly be measured, reconstructing them via structural dynamic response is essential in the field of SHM. Shape function-based methods with the ability to reduce the number of unknowns and shorten the reconstruction time significantly have attracted great interest. However, determining the optimum time unit length limits often presents a challenging task. In view of this, this paper developed an adaptive method for dynamic force reconstruction based on multiscale wavelet shape functions and time-domain deconvolution. The proposed method consists of four stages: (1) discretize the dynamic force into several time units in time domain, and approximate the local force in each unit by using wavelet scale functions at Scale 0; (2) formulate the whole response matrix by assembling the responses induced by the wavelet shape function-based forces in all time units; (3) establish the wavelet shape functions based force-response equation and obtain the Scale 0 solution; (4) lift the scale of the wavelet shape functions based force-response equation and improve the reconstruction accuracy. A simply supported beam and a two-story frame subjected to different dynamic forces are simulated to verify the feasibility and effectiveness of the proposed method. The results indicate that the proposed method can significantly reduce the number of unknowns and shorten the reconstruction time, and the reconstruction accuracy can be enhanced by lifting the scale of the wavelet shape function-based force-response equation gradually. This adaptive strategy avoids the challenging task of determining the time unit length in advance. The experimental validation of the proposed method needs to be systematically investigated in the future.



---

*Source: 8213105-2018-01-31.xml*