An open-source FEniCS-based framework for hyperelastic parameter estimation from noisy full-field data: Application to heterogeneous soft tissues

Abstract We introduce a finite-element-model-updating-based open-source framework to identify mechanical parameters of heterogeneous hyperelastic materials from in silico generated full-field data which can be downloaded here https://github.com/aflahelouneg/inverse_identification_soft_tissue. The numerical process consists in simulating an extensometer performing in vivo uniaxial tensile experiment on a soft tissue. The reaction forces and displacement fields are respectively captured by force sensor and Digital Image Correlation techniques. By means of a forward nonlinear FEM model and an inverse solver, the model parameters are estimated through a constrained optimization function with no quadratic penalty term. As a case study, our Finite Element Model Updating (FEMU) tool has been applied on a model composed of a keloid scar surrounded by healthy skin. The results show that at least 4 parameters can be accurately identified from an uniaxial test only. The originality of this work lies in two major elements. Firstly, we develop a low-cost technique able to characterize the mechanical properties of heterogeneous nonlinear hyperelastic materials. Secondly, we explore the model accuracy via a detailed study of the interplay between discretization error and the error due to measurement uncertainty. Next steps consist in identifying the real parameters and so finding the matching preferential directions of keloid scars growth.

[1]  Jung-Ryul Lee,et al.  Identification of the four orthotropic plate stiffnesses using a single open-hole tensile test , 2005 .

[2]  Francesco Ubertini,et al.  Estimation of elastic constants of thick laminated plates within a Bayesian framework , 2007 .

[3]  Zhigang Suo,et al.  Extension limit, polarization saturation, and snap-through instability of dielectric elastomers , 2011 .

[4]  J. N. Reddy,et al.  An Introduction to CONTINUUM MECHANICS with Applications , 2007 .

[5]  M. R. Molteno,et al.  Mode I–III Decomposition of the J‐integral from DIC Displacement Data , 2015 .

[6]  Josef Guttmann,et al.  In vivo characterization of mechanical tissue properties of internal organs using endoscopic microscopy and inverse finite element analysis. , 2011, Journal of biomechanics.

[7]  A. Gent A New Constitutive Relation for Rubber , 1996 .

[8]  T. Lihoreau,et al.  Multimodal investigation of a keloid scar by combining mechanical tests in vivo with diverse imaging techniques. , 2019, Journal of the mechanical behavior of biomedical materials.

[9]  Anke Stoll,et al.  Isotropic incompressible hyperelastic models for modelling the mechanical behaviour of biological tissues: a review , 2015, Biomedizinische Technik. Biomedical engineering.

[10]  Georges Limbert,et al.  Mathematical and computational modelling of skin biophysics: a review , 2017, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[11]  Gianluca Tozzi,et al.  The use of digital image correlation in the biomechanical area: a review , 2016 .

[12]  R. Ogden,et al.  A New Constitutive Framework for Arterial Wall Mechanics and a Comparative Study of Material Models , 2000 .

[13]  Kumar Vemaganti,et al.  A Bayesian approach to selecting hyperelastic constitutive models of soft tissue , 2015 .

[14]  C. Horgan,et al.  A description of arterial wall mechanics using limiting chain extensibility constitutive models , 2003, Biomechanics and modeling in mechanobiology.

[15]  Stephane Cotin,et al.  Controlling the error on target motion through real‐time mesh adaptation: Applications to deep brain stimulation , 2017, International journal for numerical methods in biomedical engineering.

[16]  A F Mak,et al.  Estimating the effective Young's modulus of soft tissues from indentation tests--nonlinear finite element analysis of effects of friction and large deformation. , 1997, Medical engineering & physics.

[17]  Herve Delingette,et al.  Real-Time Elastic Deformations of Soft Tissues for Surgery Simulation , 1999, IEEE Trans. Vis. Comput. Graph..

[18]  L. Treloar,et al.  The elasticity of a network of long-chain molecules.—III , 1943 .

[19]  Seong-Ho Kong,et al.  Handling topological changes during elastic registration , 2017, International Journal of Computer Assisted Radiology and Surgery.

[20]  C. Le Govic,et al.  Non Determined Tests as a Way to Identify Wood Elastic Parameters the Finite Element Approach , 1991 .

[21]  D. J. Montgomery,et al.  The physics of rubber elasticity , 1949 .

[22]  Erwan Verron,et al.  Comparison of Hyperelastic Models for Rubber-Like Materials , 2006 .

[23]  Stéphane Bordas,et al.  Accelerating Monte Carlo estimation with derivatives of high-level finite element models , 2017 .

[24]  P. Wolf,et al.  AN IMPROVED METHOD OF DIGITAL IMAGE CORRELATION , 1975 .

[25]  M. Giton Hyperelastic Behaviour Identification by Forward Problem Resolution : Application to a Tear Test of a Silicone-Rubber , 2006 .

[26]  E Jacquet,et al.  An open source pipeline for design of experiments for hyperelastic models of the skin with applications to keloids. , 2020, Journal of the mechanical behavior of biomedical materials.

[27]  M. Bonnet,et al.  Overview of Identification Methods of Mechanical Parameters Based on Full-field Measurements , 2008 .

[28]  Thiago Ritto,et al.  Bayesian model selection of hyperelastic models for simple and pure shear at large deformations , 2015 .

[29]  S. L. Evans,et al.  Measuring the mechanical properties of human skin in vivo using digital image correlation and finite element modelling , 2009 .

[30]  E. Pagnacco,et al.  Identification from measurements of mechanical fields by finite element model updating strategies , 2009 .

[31]  J. Kajberg,et al.  Characterisation of materials subjected to large strains by inverse modelling based on in-plane displacement fields , 2004 .

[32]  Nazim Haouchine,et al.  Impact of Soft Tissue Heterogeneity on Augmented Reality for Liver Surgery , 2015, IEEE Transactions on Visualization and Computer Graphics.

[33]  R. Meijer,et al.  Characterisation of Anisotropic and Non-linear Behaviour of Human Skin In Vivo. , 1999, Computer methods in biomechanics and biomedical engineering.

[34]  E. Somersalo,et al.  Statistical and computational inverse problems , 2004 .

[35]  Christian Duriez,et al.  Real-time simulation of contact and cutting of heterogeneous soft-tissues , 2014, Medical Image Anal..

[36]  Stéphane P. A. Bordas,et al.  Bayesian inference for the stochastic identification of elastoplastic material parameters: Introduction, misconceptions and additional insight , 2016, ArXiv.

[37]  Stéphane Roux,et al.  Estimation of Elastoplastic Parameters via Weighted FEMU and Integrated-DIC , 2015 .

[38]  Hendrik Schmidt,et al.  Discretization error when using finite element models: analysis and evaluation of an underestimated problem. , 2009, Journal of biomechanics.

[39]  Nazim Haouchine,et al.  Image-guided simulation of heterogeneous tissue deformation for augmented reality during hepatic surgery , 2013, 2013 IEEE International Symposium on Mixed and Augmented Reality (ISMAR).

[40]  Y. Gaillard,et al.  On the uniqueness of intrinsic viscoelastic properties of materials extracted from nanoindentation using FEMU , 2020, International Journal of Solids and Structures.

[41]  Mhh Marcel Meuwissen,et al.  An inverse method for the mechanical characterisation of metals , 1998 .

[42]  S. Bordas,et al.  A simple error estimator for extended finite elements , 2007 .

[43]  Zhanwei Liu,et al.  Characterization method of mechanical properties of rubber materials based on in-situ stereo finite-element-model updating , 2019, Polymer Testing.

[44]  W A Wall,et al.  Numerical identification method for the non-linear viscoelastic compressible behavior of soft tissue using uniaxial tensile tests and image registration - application to rat lung parenchyma. , 2014, Journal of the mechanical behavior of biomedical materials.

[45]  Klaus-Jürgen Bathe,et al.  On automatic mesh construction and mesh refinement in finite element analysis , 1989 .

[46]  Timon Rabczuk,et al.  Detection of material interfaces using a regularized level set method in piezoelectric structures , 2016 .

[47]  R. Ogden Large deformation isotropic elasticity – on the correlation of theory and experiment for incompressible rubberlike solids , 1972, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[48]  R. Clough,et al.  Finite element applications in the characterization of elastic solids , 1971 .

[49]  António Andrade-Campos,et al.  Comparison of inverse identification strategies for constitutive mechanical models using full-field measurements , 2018, International Journal of Mechanical Sciences.

[50]  J. Gelin,et al.  An optimization strategy for the determination of material and process parameters to avoid segregation defects during metal injection powder , 2004 .

[51]  Timon Rabczuk,et al.  An Energy Approach to the Solution of Partial Differential Equations in Computational Mechanics via Machine Learning: Concepts, Implementation and Applications , 2019, Computer Methods in Applied Mechanics and Engineering.

[52]  Jeffrey E. Bischoff,et al.  Finite element simulations of orthotropic hyperelasticity , 2002 .

[53]  Francisco J. Montáns,et al.  The mechanical behavior of skin: Structures and models for the finite element analysis , 2017 .

[54]  P. Sandoz,et al.  Ultra‐light extensometer for the assessment of the mechanical properties of the human skin in vivo , 2017, Skin research and technology : official journal of International Society for Bioengineering and the Skin (ISBS) [and] International Society for Digital Imaging of Skin (ISDIS) [and] International Society for Skin Imaging.

[55]  Martín Abadi,et al.  TensorFlow: Large-Scale Machine Learning on Heterogeneous Distributed Systems , 2016, ArXiv.

[56]  David A. Ham,et al.  Automated Derivation of the Adjoint of High-Level Transient Finite Element Programs , 2013, SIAM J. Sci. Comput..

[57]  Hugo Sol,et al.  Mixed numerical–experimental technique for orthotropic parameter identification using biaxial tensile tests on cruciform specimens , 2007 .

[58]  R. Ogawa Keloid and hypertrophic scarring may result from a mechanoreceptor or mechanosensitive nociceptor disorder. , 2008, Medical hypotheses.

[59]  Frédéric Bos,et al.  Identification of constitutive laws for wood-based panels by means of an inverse method , 2002 .

[60]  D. Arola,et al.  Applications of digital image correlation to biological tissues. , 2004, Journal of biomedical optics.

[61]  M. Gilchrist,et al.  Mechanical characterization of brain tissue in compression at dynamic strain rates. , 2012, Journal of the mechanical behavior of biomedical materials.

[62]  E M Arruda,et al.  Finite element modeling of human skin using an isotropic, nonlinear elastic constitutive model. , 2000, Journal of biomechanics.

[63]  Stéphane Bordas,et al.  Bayesian inference to identify parameters in viscoelasticity , 2018 .

[64]  Stéphane Cotin,et al.  Real-Time Error Control for Surgical Simulation , 2016, IEEE Transactions on Biomedical Engineering.

[65]  Pierre Badel,et al.  Anisotropic and hyperelastic identification of in vitro human arteries from full-field optical measurements. , 2010, Journal of biomechanics.

[66]  Dimitri P. Bertsekas,et al.  Constrained Optimization and Lagrange Multiplier Methods , 1982 .

[67]  Cornelius O. Horgan,et al.  The remarkable Gent constitutive model for hyperelastic materials , 2015 .

[68]  Jd Jan Janssen,et al.  Determination of the elasto-plastic properties of aluminium using a mixed numerical–experimental method , 1998 .

[69]  Yueqi Wang Uncertainty Quantification of Digital Image Correlation and the Impact on Material Identification , 2015 .

[70]  Edna Suarez Pozos The effect of skin tension on the formation of keloid scars , 2014 .

[71]  Nizar Harb Identification inverse de paramètres biomécaniques en hyperélasticité anisotrope , 2013 .

[72]  Franz Chouly,et al.  Quantifying discretization errors for soft-tissue simulation in computer assisted surgery: a preliminary study , 2018, Applied Mathematical Modelling.

[73]  Michel Bornert,et al.  Experimental characterization of the local strain field in a heterogeneous elastoplastic material , 1994 .

[74]  P Sandoz,et al.  Intra‐ and inter‐individual variability in the mechanical properties of the human skin from in vivo measurements on 20 volunteers , 2017, Skin research and technology : official journal of International Society for Bioengineering and the Skin (ISBS) [and] International Society for Digital Imaging of Skin (ISDIS) [and] International Society for Skin Imaging.

[75]  E Jacquet,et al.  Numerical analysis of keloid scar in the presternal area , 2012, Computer methods in biomechanics and biomedical engineering.

[76]  Jin Kim,et al.  Characterisation of failure in human aortic tissue using digital image correlation , 2011 .

[77]  R. Ogawa,et al.  The Relationship Between Keloid Growth Pattern and Stretching Tension: Visual Analysis Using the Finite Element Method , 2008, Annals of plastic surgery.

[78]  Stéphane Bordas,et al.  A Tutorial on Bayesian Inference to Identify Material Parameters in Solid Mechanics , 2019, Archives of Computational Methods in Engineering.

[79]  Stéphane Cotin,et al.  Quantifying the uncertainty in a hyperelastic soft tissue model with stochastic parameters , 2018, Applied Mathematical Modelling.

[80]  Luciano Lamberti,et al.  Mechanical characterization of hyperelastic materials with fringe projection and optimization techniques , 2006 .

[81]  O. Yeoh Some Forms of the Strain Energy Function for Rubber , 1993 .