Please wait a minute...
Frontiers of Mechanical Engineering

ISSN 2095-0233

ISSN 2095-0241(Online)

CN 11-5984/TH

Postal Subscription Code 80-975

2018 Impact Factor: 0.989

Front. Mech. Eng.    2023, Vol. 18 Issue (3) : 46    https://doi.org/10.1007/s11465-023-0762-2
RESEARCH ARTICLE
Non-convex sparse optimization-based impact force identification with limited vibration measurements
Lin CHEN1,2, Yanan WANG1,2(), Baijie QIAO1,2, Junjiang LIU1,2, Wei CHENG1,2, Xuefeng CHEN1,2
1. National Key Laboratory of Aerospace Power System and Plasma Technology, Xi’an Jiaotong University, Xi’an 710049, China
2. School of Mechanical Engineering, Xi’an Jiaotong University, Xi’an 710049, China
 Download: PDF(11829 KB)   HTML
 Export: BibTeX | EndNote | Reference Manager | ProCite | RefWorks
Abstract

Impact force identification is important for structure health monitoring especially in applications involving composite structures. Different from the traditional direct measurement method, the impact force identification technique is more cost effective and feasible because it only requires a few sensors to capture the system response and infer the information about the applied forces. This technique enables the acquisition of impact locations and time histories of forces, aiding in the rapid assessment of potentially damaged areas and the extent of the damage. As a typical inverse problem, impact force reconstruction and localization is a challenging task, which has led to the development of numerous methods aimed at obtaining stable solutions. The classical 2 regularization method often struggles to generate sparse solutions. When solving the under-determined problem, 2 regularization often identifies false forces in non-loaded regions, interfering with the accurate identification of the true impact locations. The popular 1 sparse regularization, while promoting sparsity, underestimates the amplitude of impact forces, resulting in biased estimations. To alleviate such limitations, a novel non-convex sparse regularization method that uses the non-convex 12 penalty, which is the difference of the 1 and 2 norms, as a regularizer, is proposed in this paper. The principle of alternating direction method of multipliers (ADMM) is introduced to tackle the non-convex model by facilitating the decomposition of the complex original problem into easily solvable subproblems. The proposed method named 12-ADMM is applied to solve the impact force identification problem with unknown force locations, which can realize simultaneous impact localization and time history reconstruction with an under-determined, sparse sensor configuration. Simulations and experiments are performed on a composite plate to verify the identification accuracy and robustness with respect to the noise of the 12-ADMM method. Results indicate that compared with other existing regularization methods, the 12-ADMM method can simultaneously reconstruct and localize impact forces more accurately, facilitating sparser solutions, and yielding more accurate results.

Keywords impact force identification      inverse problem      sparse regularization      under-determined condition      alternating direction method of multipliers     
Corresponding Author(s): Yanan WANG   
Just Accepted Date: 04 July 2023   Issue Date: 13 October 2023
 Cite this article:   
Lin CHEN,Yanan WANG,Baijie QIAO, et al. Non-convex sparse optimization-based impact force identification with limited vibration measurements[J]. Front. Mech. Eng., 2023, 18(3): 46.
 URL:  
https://academic.hep.com.cn/fme/EN/10.1007/s11465-023-0762-2
https://academic.hep.com.cn/fme/EN/Y2023/V18/I3/46
  Fig.A1
  Fig.A2
Fig.1  3D plots of different penalties and their corresponding contour maps: (a) 3D plots of ?2 norm, ?1 norm, ?1?2 norm, and ?0 norm and (b) corresponding contour maps.
  Fig.A3
Elastic modulusShear modulusPoisson’s ratioDensityLayer-ups
E1 = 135.0 GPa, E2 = 8.8 GPa, E3 = 8.8 GPaG12 = 4.47 GPa, G23 = 3.00 GPa, G13 = 4.47 GPav12 = 0.3, v23 = 0.4, v13 = 0.3ρ = 1560 kg/m3[45°/0°/?45°/90°]
Tab.1  Material properties of the composite plate
Fig.2  Composite plate with two opposite edges clamped: (a) schematic diagram of an applied impact force on the plate and (b) distribution of strain sensors and potential impact locations. Fifteen potential impact locations are considered, and the selected measurement positions are S8, S10, S14, and S18.
Fig.3  Reconstructed impact time histories at all monitored locations at 20 dB noise level when (a) real impact acts on P3 and (b) real impact acts on P10. ADMM: alternating direction method of multipliers.
Fig.4  Impact force identification results denoted by accuracy indicators at eight randomly chosen locations at 20 dB noise level by ?1?2-alternating direction method of multipliers (ADMM), ?1-ADMM, and ?2-ADMM: (a) global relative error, (b) local relative error, (c) peak relative error, and (d) localization error.
PositionMethodGRE/%LRE/%PRE/%LE/%Computing time/s
SineTriangleGaussian
P9?1?2-ADMM10.538.643.894.583.8515.45280.07
?1-ADMM15.0311.949.337.376.6921.92292.65
?2-ADMM85.9071.8870.2973.4272.8584.4896.88
P11?1?2-ADMM11.039.595.352.801.4916.11272.45
?1-ADMM15.1413.2110.305.984.2120.78273.10
?2-ADMM94.2188.9589.2986.8387.5090.8540.14
Tab.2  Identification accuracy indicators and computing time of ?1?2-ADMM, ?1-ADMM, and ?2-ADMM for the multimorphology continuous impact force at different locations at 20 dB noise level
Fig.5  Identification results of the multimorphology continuous impact force at P9 under 20 dB noise level by ?1?2-alternating direction method of multipliers (ADMM), ?1-ADMM, and ?2-ADMM: (a) time history reconstruction and (b) localization results.
Fig.6  Identification results of the multimorphology continuous impact force at P11 under 20 dB noise level by ?1?2-alternating direction method of multipliers (ADMM), ?1-ADMM, and ?2-ADMM: (a) time–history reconstruction and (b) localization results.
Fig.7  Four accuracy indicators of identification results for the same single impact force via ?1?2-alternating direction method of multipliers (ADMM), ?1-ADMM, and ?2-ADMM at different noise levels: (a) impact at P1 and (b) impact at P9. SNR: signal-noise ratio.
Fig.8  Experimental setup of the composite plate with two opposite edges clamped. Four sensors (S1, S2, S7, and S9) are selected to monitor the 15 potential impact locations
Accuracy indicators and computing time
MethodImpact at P1Impact at P15
GRE/%LRE/%PRE/%LE/%Computing time/sGRE/%LRE/%PRE/%LE/%Computing time/s
?1?2-ADMM5.995.900.862.2693.835.284.630.836.16267.35
?1/2-ADMM6.445.441.157.22253.916.154.961.324.56376.05
?2/3-ADMM8.105.892.2312.70261.045.924.981.214.04298.46
?1-ADMM7.967.273.718.86101.117.596.063.1511.37217.71
?2-ADMM75.8961.6366.4677.9746.2486.9177.3175.7783.0971.51
Tab.3  Identification accuracy indicators and computing times of ?1?2-ADMM, ?1/2-ADMM, ?2/3-ADMM, ?1-ADMM, and ?2-ADMM for the single-impact forces applied to P1 and P15
Fig.9  Monte Carlo generalized stein unbiased risk estimate (MC-GSURE) curves of ?1?2-alternating direction method of multipliers in single-impact cases when exerted at: (a) P1 and (b) P15.
Fig.10  Identification results of the impact force applied to P1 by ?1?2-alternating direction method of multipliers (ADMM), ?1-ADMM, ?2-ADMM, ?1/2-ADMM, and ?2/3-ADMM: (a) time history reconstruction results, (b) localization results of ?1?2-ADMM, ?1-ADMM, and ?2-ADMM, (c) localization results of ?1/2-ADMM and ?2/3-ADMM, and (d) convergence curves.
Fig.11  Identification results of the impact force applied to P15 by ?1?2-alternating direction method of multipliers (ADMM), ?1-ADMM, ?2-ADMM, ?1/2-ADMM, and ?2/3-ADMM: (a) time history reconstruction results, (b) localization results of ?1?2-ADMM, ?1-ADMM, and ?2-ADMM, (c) localization results of ?1/2-ADMM and ?2/3-ADMM, and (d) convergence curves.
Accuracy indicators and computing time
MethodImpact at P5Impact at P6
GRE/%LRE/%PRE/%LE/%Computing time/sGRE/%LRE/%PRE/%LE/%Computing time/s
Peak1Peak2Peak1Peak2
?1?2-ADMM11.769.880.480.9318.17219.715.944.950.404.297.76229.85
?1/2-ADMM19.8017.746.8415.3515.60287.4311.498.944.885.2611.44282.67
?2/3-ADMM16.4614.357.8911.6721.25188.9610.178.122.815.1710.43193.90
?1-ADMM13.5711.294.564.7121.25189.588.537.224.628.838.56213.31
?2-ADMM85.3974.9677.9077.0282.1945.6982.2068.8165.5069.8581.10101.03
Tab.4  Identification accuracy indicators and computing times of ?1?2-ADMM, ?1/2-ADMM, ?2/3-ADMM, ?1-ADMM, and ?2-ADMM for the continuous-impact forces applied to P5 and P6
Fig.12  Identification results of the impact force applied to P5 by ?1?2-alternating direction method of multipliers (ADMM), ?1-ADMM, ?2-ADMM, ?1/2-ADMM, and ?2/3-ADMM: (a) time history reconstruction results, (b) localization results of ?1?2-ADMM, ?1-ADMM, and ?2-ADMM, (c) localization results of ?1/2-ADMM and ?2/3-ADMM, and (d) convergence curves.
Fig.13  Identification results of the impact force applied to P6 by ?1?2-alternating direction method of multipliers (ADMM), ?1-ADMM, ?2-ADMM, ?1/2-ADMM, and ?2/3-ADMM: (a) time history reconstruction results, (b) localization results of ?1?2-ADMM, ?1-ADMM, and ?2-ADMM, (c) localization results of ?1/2-ADMM and ?2/3-ADMM, and (d) convergence curves.
Abbreviations
ADMMAlternating direction method of multipliers
BVIDBarely visible impact damage
GREGlobal relative error
IRFImpulse response function
LELocalization error
LRELocal relative error
MC-GSUREMonte Carlo generalized stein unbiased risk estimate
PREPeak relative error
SISOSingle-input single-output
SNRSignal-noise ratio
Variables
a(t)Impulse response function
aij(t)Impulse response function between the output position i and the input location j
aij(ω)Frequency response function between the output position i and the input location j
ATransfer matrix of the multiple-input multiple-output dynamic system
AsTransfer matrix of the single-input single-output dynamic system
BAmplitude of the Gaussian-shaped impact force
cAll-ones vector
CDamping matrix
eRandom noise in measurements
EElastic modulus
fiElements in the vector f
f(t)Impact force excitation function
fForce vector of the multiple-input multiple-output dynamic system
f^Estimated vector of f
fpActual force vector at the impact position p
f^pEstimated force vector at the impact position p
f^MLMaximum likelihood estimation of the vector f
fsForce vector of the single-input single-output dynamic system
g(f)General representation function of penalty terms
GShear modulus
hAn additional vector for variable splitting
h^Estimated vector of h
iLooping variable within the summation operation
IIdentity matrix
kNumber of iterations
KStiffness matrix
mNumber of measurement responses
MMass matrix
nNumber of impact force excitations
nfTotal number of potential impact force locations
NData length of the discretized impulse response function
NmaxMaximum number of iterations
O(nN)Computational complexity of n × N
pSerial number of the location subjected to impact force
qNorm parameter defined in R+
QProjection matrix
rGaussian white noise vector
R(f)General expression for calculating the norm of vector f
s(t)System response
sResponse vector of the multiple-input multiple-output dynamic system
s~Noisy response vector
siResponse vector at a certain position i
ssResponse vector of the single-input single-output dynamic system
tTime
t0Occurrence time instant of the impact
ΔtSampling interval
TImpact duration
uSufficient statistic of the model Eq. (8)
wmaxElement with the largest absolute value in the vector·w
wIntermediate vector defined as w = f (k+1) + z(k)
xλ(u)Solution result of Eq. (8) when f = u
y(f)Proximal operator
δSmall positive parameter
δLagrange multiplier vector
εIteration termination threshold
λRegularization parameter
ρA positive penalty parameter
σStandard deviation of the vector s
σnStandard deviation of noise in the measurements
τTime delayed operator
νPossion’s ratio
ΓThreshold value defined as Γ = λ/ρ
  
1 J Wu , X B Xu , C Liu , C Deng , X Y Shao . Lamb wave-based damage detection of composite structures using deep convolutional neural network and continuous wavelet transform. Composite Structures, 2021, 276: 114590
https://doi.org/10.1016/j.compstruct.2021.114590
2 X J Gu , X Z Su , J Wang , Y J Xu , J H Zhu , W H Zhang . Improvement of impact resistance of plain-woven composite by embedding superelastic shape memory alloy wires. Frontiers of Mechanical Engineering, 2020, 15(4): 547–557
https://doi.org/10.1007/s11465-020-0595-1
3 D Han , X Jia , H J Zhang , X G Gao , X Han , L Sun , Z K Zheng , L Zhang , F Wang , Y D Song . Foreign object damage and post-impact tensile behavior of plain-woven SiC/SiC composites. Composite Structures, 2022, 295: 115767
https://doi.org/10.1016/j.compstruct.2022.115767
4 T Singh , S Sehgal . Structural health monitoring of composite materials. Archives of Computational Methods in Engineering, 2022, 29(4): 1997–2017
https://doi.org/10.1007/s11831-021-09666-8
5 X F Chen , S B Wang , B J Qiao , Q Chen . Basic research on machinery fault diagnostics: past, present, and future trends. Frontiers of Mechanical Engineering, 2018, 13(2): 264–291
https://doi.org/10.1007/s11465-018-0472-3
6 R X Liu , E Dobriban , Z C Hou , K Qian . Dynamic load identification for mechanical systems: a review. Archives of Computational Methods in Engineering, 2022, 29(2): 831–863
https://doi.org/10.1007/s11831-021-09594-7
7 L WangY R LiuH Y Xu. Review: recent developments in dynamic load identification for aerospace vehicles considering multi-source uncertainties. Transactions of Nanjing University of Aeronautics & Astronautics, 2021, 38(2): 271–287 (in Chinese)
8 E Jacquelin , A Bennani , P Hamelin . Force reconstruction: analysis and regularization of a deconvolution problem. Journal of Sound and Vibration, 2003, 265(1): 81–107
https://doi.org/10.1016/S0022-460X(02)01441-4
9 Q F Li , Q H Lu . Impact localization and identification under a constrained optimization scheme. Journal of Sound and Vibration, 2016, 366: 133–148
https://doi.org/10.1016/j.jsv.2015.12.010
10 G Yan , H Sun , O Büyüköztürk . Impact load identification for composite structures using Bayesian regularization and unscented Kalman filter. Structural Control and Health Monitoring, 2017, 24(5): e1910
https://doi.org/10.1002/stc.1910
11 M Aucejo , Smet O De , J F Deü . On a space-time regularization for force reconstruction problems. Mechanical Systems and Signal Processing, 2019, 118: 549–567
https://doi.org/10.1016/j.ymssp.2018.09.002
12 A RezayatV NassiriB De PauwJ ErtveldtS VanlanduitP Guillaume. Identification of dynamic forces using group-sparsity in frequency domain. Mechanical Systems and Signal Processing, 2016, 70–71: 756–768
13 A Candelieri , F Archetti . Sparsifying to optimize over multiple information sources: an augmented gaussian process based algorithm. Structural and Multidisciplinary Optimization, 2021, 64(1): 239–255
https://doi.org/10.1007/s00158-021-02882-7
14 M A T Figueiredo , R D Nowak , S J Wright . Gradient projection for sparse reconstruction: application to compressed sensing and other inverse problems. IEEE Journal of Selected Topics in Signal Processing, 2007, 1(4): 586–597
https://doi.org/10.1109/JSTSP.2007.910281
15 S J Kim , K Koh , M Lustig , S Boyd , D Gorinevsky . An interior-point method for large-scale ℓ1-regularized least squares. IEEE Journal of Selected Topics in Signal Processing, 2007, 1(4): 606–617
https://doi.org/10.1109/JSTSP.2007.910971
16 I Daubechies , M Defrise , C De Mol . An iterative thresholding algorithm for linear inverse problems with a sparsity constraint. Communications on Pure and Applied Mathematics, 2004, 57(11): 1413–1457
https://doi.org/10.1002/cpa.20042
17 D GinsbergC P Fritzen. New approach for impact detection by finding sparse solution. In: Proceedings of ISMA. Leuven: Katholieke Universiteit Leuven, 2014, 2043–2056
18 B J Qiao , X W Zhang , J W Gao , X F Chen . Impact-force sparse reconstruction from highly incomplete and inaccurate measurements. Journal of Sound and Vibration, 2016, 376: 72–94
https://doi.org/10.1016/j.jsv.2016.04.040
19 M Aucejo , O De Smet . A space-frequency multiplicative regularization for force reconstruction problems. Mechanical Systems and Signal Processing, 2018, 104: 1–18
https://doi.org/10.1016/j.ymssp.2017.10.027
20 C D Pan , Z P Chen . Elimination of accelerometer mass loading effects in sparse identification of impact forces. Mechanical Systems and Signal Processing, 2023, 191: 110178
https://doi.org/10.1016/j.ymssp.2023.110178
21 Z B Xu , X Y Chang , F M Xu , H Zhang . L1/2 regularization: a thresholding representation theory and a fast solver. IEEE Transactions on Neural Networks and Learning Systems, 2012, 23(7): 1013–1027
https://doi.org/10.1109/TNNLS.2012.2197412
22 D Trehan. Non-convex optimization: a review. In: Proceedings of 2020 the 4th International Conference on Intelligent Computing and Control Systems. Madurai: IEEE, 2020, 418–423
23 Z Y Zha , B H Wen , X Yuan , S Ravishankar , J T Zhou , C Zhu . Learning nonlocal sparse and low-rank models for image compressive sensing: nonlocal sparse and low-rank modeling. IEEE Signal Processing Magazine, 2023, 40(1): 32–44
https://doi.org/10.1109/MSP.2022.3217936
24 Z Y Zha , B H Wen , X Yuan , J T Zhou , C Zhu . Image restoration via reconciliation of group sparsity and low-rank models. IEEE Transactions on Image Processing, 2021, 30: 5223–5238
https://doi.org/10.1109/TIP.2021.3078329
25 R ChartrandW T Yin. Nonconvex sparse regularization and splitting algorithms. In: Glowinski R, Osher S J, Yin W T, eds. Splitting Methods in Communication, Imaging, Science, and Engineering. Cham: Springer, 2016, 237–249
26 B J Qiao , C Y Ao , Z Mao , X F Chen . Non-convex sparse regularization for impact force identification. Journal of Sound and Vibration, 2020, 477: 115311
https://doi.org/10.1016/j.jsv.2020.115311
27 M Aucejo , O De Smet . A generalized multiplicative regularization for input estimation. Mechanical Systems and Signal Processing, 2021, 157: 107637
https://doi.org/10.1016/j.ymssp.2021.107637
28 J J Liu , B J Qiao , Y N Wang , W F He , X F Chen . Non-convex sparse regularization via convex optimization for impact force identification. Mechanical Systems and Signal Processing, 2023, 191: 110191
https://doi.org/10.1016/j.ymssp.2023.110191
29 E Esser , Y F Lou , J Xin . A method for finding structured sparse solutions to nonnegative least squares problems with applications. SIAM Journal on Imaging Sciences, 2013, 6(4): 2010–2046
https://doi.org/10.1137/13090540X
30 P H Yin , Y F Lou , Q He , J Xin . Minimization of 1−2 for compressed sensing. SIAM Journal on Scientific Computing, 2015, 37(1): A536–A563
https://doi.org/10.1137/140952363
31 H Sun , D M Feng , Y Liu , M Q Feng . Statistical regularization for identification of structural parameters and external loadings using state space models. Computer-Aided Civil and Infrastructure Engineering, 2015, 30(11): 843–858
https://doi.org/10.1111/mice.12169
32 D Ginsberg , M Ruby , C P Fritzen . Load identification approach based on basis pursuit denoising algorithm. Journal of Physics: Conference Series, 2015, 628(1): 012030
https://doi.org/10.1088/1742-6596/628/1/012030
33 S Boyd , N Parikh , E Chu , B Peleato , J Eckstein . Distributed optimization and statistical learning via the alternating direction method of multipliers. Foundations and Trends in Machine Learning, 2011, 3(1): 1–122
https://doi.org/10.1561/2200000016
34 L Meirovitch. Elements of Vibration Analysis. New York: McGraw-Hill Book Co., 1986
35 L Jankowski . Off-line identification of dynamic loads. Structural and Multidisciplinary Optimization, 2009, 37(6): 609–623
https://doi.org/10.1007/s00158-008-0249-0
36 Z Y ZhaX LiuX H HuangH L ShiY Y XuQ WangL TangX G Zhang. Analyzing the group sparsity based on the rank minimization methods. In: Proceedings of 2017 IEEE International Conference on Multimedia and Expo. Hong Kong: IEEE, 2017, 883–888
37 Z Y Zha , X G Zhang , Y Wu , Q Wang , X Liu , L Tang , X Yuan . Non-convex weighted ℓp nuclear norm based ADMM framework for image restoration. Neurocomputing, 2018, 311: 209–224
https://doi.org/10.1016/j.neucom.2018.05.073
38 X W Li , H T Zhao , Z Chen , Q B Wang , J A Chen , D P Duan . Force identification based on a comprehensive approach combining Taylor formula and acceleration transmissibility. Inverse Problems in Science and Engineering, 2018, 26(11): 1612–1632
https://doi.org/10.1080/17415977.2017.1417407
39 L Wang , J K Liu , Z R Lu . Bandlimited force identification based on sinc-dictionaries and Tikhonov regularization. Journal of Sound and Vibration, 2020, 464: 114988
https://doi.org/10.1016/j.jsv.2019.114988
40 B J Qiao , J J Liu , J X Liu , Z B Yang , X F Chen . An enhanced sparse regularization method for impact force identification. Mechanical Systems and Signal Processing, 2019, 126: 341–367
https://doi.org/10.1016/j.ymssp.2019.02.039
41 Y F Lou , M Yan . Fast L1–L2 minimization via a proximal operator. Journal of Scientific Computing, 2018, 74(2): 767–785
https://doi.org/10.1007/s10915-017-0463-2
42 X T Chang , Y J Yan , Y F Wu . Study on solving the ill-posed problem of force load reconstruction. Journal of Sound and Vibration, 2019, 440: 186–201
https://doi.org/10.1016/j.jsv.2018.10.017
43 W F Cao , J Sun , Z B Xu . Fast image deconvolution using closed-form thresholding formulas of Lq (q = 1/2, 2/3) regularization. Journal of Visual Communication and Image Representation, 2013, 24(1): 31–41
https://doi.org/10.1016/j.jvcir.2012.10.006
44 Z Y Zha , X Yuan , B H Wen , J T Zhou , J C Zhang , C Zhu . A benchmark for sparse coding: when group sparsity meets rank minimization. IEEE Transactions on Image Processing, 2020, 29: 5094–5109
https://doi.org/10.1109/TIP.2020.2972109
45 Y C Eldar . Generalized SURE for exponential families: applications to regularization. IEEE Transactions on Signal Processing, 2009, 57(2): 471–481
https://doi.org/10.1109/TSP.2008.2008212
46 D L Donoho . De-noising by soft-thresholding. IEEE Transactions on Information Theory, 1995, 41(3): 613–627
https://doi.org/10.1109/18.382009
47 S Ramani , T Blu , M Unser . Monte-carlo sure: a black-box optimization of regularization parameters for general denoising algorithms. IEEE Transactions on Image Processing, 2008, 17(9): 1540–1554
https://doi.org/10.1109/TIP.2008.2001404
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed