An inverse model and mathematical solution for inferring viscoelastic properties and dynamic deformations of heterogeneous structures

Abstract We proposed an adjoint-based inverse method and mathematical solutions using Lagrange multiplier theorem for inferring dynamic viscoelastic properties of heterogeneous structures with displacement observations to improve computation reliability and efficiency. Existing methods such as the one using finite-difference approximated gradient may not be efficient and accurate enough for considering complex situations of the coupled effects of dynamic loading, material deformation memory, and structural heterogeneity. The proposed method derives both gradient and Hessian mathematically to satisfy the first-order necessary and second-order sufficient optimal conditions, which has great advantages compared to other approaches. We also proposed robust numerical algorithms for accurate and fast computations, including a regularization to control reasonable parameter ranges, a Newton’s method with Hessian function to determine search direction for fast convenience, and a modified Armijo rule to find a stable step length. We developed a Galerkin time-domain finite-element method for numerical solutions of resulting partial differential equations, and validated the model for a layered structure under dynamic loading tests. Results indicate that the proposed method has greatly improved computation accuracy and speed as compared to existing approaches. It may be applied to a broad range of heterogeneous materials and structures at different length and time scales.

[1]  Louay N. Mohammad,et al.  Modeling Asphalt Pavement Rutting under Accelerated Testing , 2008 .

[2]  Jianlin Xia,et al.  New Efficient and Robust HSS Cholesky Factorization of SPD Matrices , 2012, SIAM J. Matrix Anal. Appl..

[3]  Jintai Chung,et al.  A family of single-step Houbolt time integration algorithms for structural dynamics , 1994 .

[4]  Robert Y. Liang,et al.  Efficient computational algorithms for forward and backward analysis of a dynamic pavement system , 1998 .

[5]  George K. Chang,et al.  Experimental and numerical study of asphalt material geospatial heterogeneity with intelligent compaction technology on roads , 2014 .

[6]  R. Al-Khoury,et al.  Spectral element technique for efficient parameter identification of layered media. Part III: viscoelastic aspects , 2002 .

[7]  Jim E. Jones,et al.  Newton–Krylov-multigrid solvers for large-scale, highly heterogeneous, variably saturated flow problems , 2001 .

[8]  Jose Herskovits,et al.  Characterisation by Inverse Techniques of Elastic, Viscoelastic and Piezoelectric Properties of Anisotropic Sandwich Adaptive Structures , 2010 .

[9]  Panayot S. Vassilevski,et al.  Approved for Public Release; Further Dissemination Unlimited on Mesh-independent Convergence of an Inexact Newton-multigrid Algorithm on Mesh–independent Convergence of an Inexact Newton–multigrid Algorithm , 2022 .

[10]  A. Bhasin,et al.  Constitutive modeling of the nonlinearly viscoelastic response of asphalt binders; incorporating three-dimensional effects , 2013 .

[11]  N. Petra,et al.  Model Variational Inverse Problems Governed by Partial Differential Equations , 2011 .

[12]  S. W. Park,et al.  Methods of interconversion between linear viscoelastic material functions. Part I-a numerical method based on Prony series , 1999 .

[13]  José Herskovits,et al.  Interior point algorithms for nonlinear constrained least squares problems , 2004 .

[14]  Jorge A Prozzi,et al.  Static versus viscoelastic wave propagation approach for simulating loading effects on flexible pavement structure , 2014 .

[15]  Robert L. Lytton,et al.  COMPARISON OF COMPUTER PREDICTIONS AND FIELD DATA FOR DYNAMIC ANALYSIS OF FALLING WEIGHT DEFLECTOMETER DATA , 1991 .

[16]  Omar Ghattas,et al.  Site characterization using full waveform inversion , 2013 .

[17]  O. Ghattas,et al.  An inexact Gauss-Newton method for inversion of basal sliding and rheology parameters in a nonlinear Stokes ice sheet model , 2012, Journal of Glaciology.

[18]  M. Mitchell,et al.  Measurement and evaluation of asphalt concrete thermal expansion and contraction , 2008 .

[19]  M. Solaimanian,et al.  Modeling temperature distribution and thermal property of asphalt concrete for laboratory testing applications , 2010 .

[20]  Jorge A Prozzi,et al.  A time-domain finite element method for dynamic viscoelastic solution of layered-half-space responses under loading pulses , 2015 .

[21]  Mansour Solaimanian,et al.  Modelling linear viscoelastic properties of asphalt concrete by the Huet–Sayegh model , 2009 .

[22]  Georg Stadler,et al.  Parallel scalable adjoint-based adaptive solution of variable-viscosity Stokes flow problems , 2009 .

[23]  M Fink,et al.  Measurement of viscoelastic properties of homogeneous soft solid using transient elastography: An inverse problem approach , 2004 .

[24]  J. Tronto,et al.  Conducting Polymers / Layered Double Hydroxides Intercalated Nanocomposites , 2013 .

[25]  K. Schindhelm,et al.  Elastic and viscoelastic properties of porcine subdermal fat using MRI and inverse FEA , 2010, Biomechanics and modeling in mechanobiology.

[26]  Eyal Levenberg,et al.  Inverse analysis of viscoelastic pavement properties using data from embedded instrumentation , 2013 .

[27]  Souhir Zghal,et al.  Model reduction methods for viscoelastic sandwich structures in frequency and time domains , 2015 .

[28]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[29]  Ruogang Zhao,et al.  Comparison of analytical and inverse finite element approaches to estimate cell viscoelastic properties by micropipette aspiration. , 2009, Journal of biomechanics.

[30]  Qinwu Xu,et al.  A finite-element and Newton-Raphson method for inverse computing multilayer moduli , 2014 .

[31]  Arash Fathi Full-waveform inversion in three-dimensional PML-truncated elastic media : theory, computations, and field experiments , 2015 .

[32]  A. Fathi,et al.  Full-waveform inversion in three-dimensional PML-truncated elastic media , 2015, 1504.08340.

[33]  Fulin Lei,et al.  Inverse analysis of constitutive models: biological soft tissues. , 2007, Journal of biomechanics.

[35]  R. Christensen,et al.  Theory of Viscoelasticity , 1971 .

[36]  James F. Greenleaf,et al.  Inverse estimation of viscoelastic material properties for solids immersed in fluids using vibroacoustic techniques , 2007 .

[37]  Santiago Badia,et al.  A Highly Scalable Parallel Implementation of Balancing Domain Decomposition by Constraints , 2014, SIAM J. Sci. Comput..

[38]  Peter Deuflhard,et al.  Scientific Computing with Ordinary Differential Equations , 2002 .

[39]  Victor M. Calo,et al.  Mathematical modeling of coupled drug and drug-encapsulated nanoparticle transport in patient-specific coronary artery walls , 2012 .

[40]  Eyal Levenberg,et al.  Viscoelastic Genetic Algorithm for Inverse Analysis of Asphalt Layer Properties from Falling Weight Deflections , 2013 .

[41]  José Herskovits,et al.  Estimation of piezoelastic and viscoelastic properties in laminated structures , 2009 .

[42]  T. Lu,et al.  Biothermomechanics of skin tissues , 2008 .

[43]  Loukas F. Kallivokas,et al.  Time‐domain hybrid formulations for wave simulations in three‐dimensional PML‐truncated heterogeneous media , 2015 .

[44]  Robert L. Lytton,et al.  SYSTEM IDENTIFICATION METHOD FOR BACKCALCULATING PAVEMENT LAYER PROPERTIES , 1993 .