Research Papers

Numerical Approximation of Elasticity Tensor Associated With Green-Naghdi Rate

[+] Author and Article Information
Haofei Liu

Department of Mechanics,
Tianjin University,
92 Weijin Road,
Tianjin 300072, China

Wei Sun

Tissue Mechanics Laboratory,
The Wallace H. Coulter Department of Biomedical Engineering,
Georgia Institute of Technology,
Technology Enterprise Park,
Room 206, 387 Technology Circle,
Atlanta, GA 30313-2412
e-mail: wei.sun@bme.gatech.edu

Manuscript received July 26, 2016; final manuscript received May 3, 2017; published online June 16, 2017. Assoc. Editor: Hai-Chao Han.

J Biomech Eng 139(8), 081007 (Jun 16, 2017) (8 pages) Paper No: BIO-16-1317; doi: 10.1115/1.4036829 History: Received July 26, 2016; Revised May 03, 2017

Objective stress rates are often used in commercial finite element (FE) programs. However, deriving a consistent tangent modulus tensor (also known as elasticity tensor or material Jacobian) associated with the objective stress rates is challenging when complex material models are utilized. In this paper, an approximation method for the tangent modulus tensor associated with the Green-Naghdi rate of the Kirchhoff stress is employed to simplify the evaluation process. The effectiveness of the approach is demonstrated through the implementation of two user-defined fiber-reinforced hyperelastic material models. Comparisons between the approximation method and the closed-form analytical method demonstrate that the former can simplify the material Jacobian evaluation with satisfactory accuracy while retaining its computational efficiency. Moreover, since the approximation method is independent of material models, it can facilitate the implementation of complex material models in FE analysis using shell/membrane elements in abaqus.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.


Gasser, T. C. , Ogden, R. W. , and Holzapfel, G. A. , 2006, “ Hyperelastic Modelling of Arterial Layers With Distributed Collagen Fibre Orientations,” J. R. Soc. Interface, 3(6), pp. 15–35. [CrossRef] [PubMed]
Fan, R. , and Sacks, M. S. , 2014, “ Simulation of Planar Soft Tissues Using a Structural Constitutive Model: Finite Element Implementation and Validation,” J. Biomech., 47(9), pp. 2043–2054. [CrossRef] [PubMed]
Holzapfel, G. A. , and Gasser, T. C. , 2001, “ A Viscoelastic Model for Fiber-Reinforced Composites at Finite Strains: Continuum Basis, Computational Aspects and Applications,” Comput. Methods Appl. Mech. Eng., 190(34), pp. 4379–4403. [CrossRef]
Genet, M. , Rausch, M. K. , Lee, L. C. , Choy, S. , Zhao, X. , Kassab, G. S. , Kozerke, S. , Guccione, J. M. , and Kuhl, E. , 2015, “ Heterogeneous Growth-Induced Prestrain in the Heart,” J. Biomech., 48(10), pp. 2080–2089. [CrossRef] [PubMed]
Lubarda, V. A. , and Hoger, A. , 2002, “ On the Mechanics of Solids With a Growing Mass,” Int. J. Solids Struct., 39(18), pp. 4627–4664. [CrossRef]
Nordsletten, D. , McCormick, M. , Kilner, P. J. , Hunter, P. , Kay, D. , and Smith, N. P. , 2011, “ Fluid–Solid Coupling for the Investigation of Diastolic and Systolic Human Left Ventricular Function,” Int. J. Numer. Methods Biomed. Eng., 27(7), pp. 1017–1039. [CrossRef]
Martin, C. , and Sun, W. , 2013, “ Modeling of Long-Term Fatigue Damage of Soft Tissue With Stress Softening and Permanent Set Effects,” Biomech. Model. Mechanobiol., 12(4), pp. 645–655. [CrossRef] [PubMed]
Haj-Ali, R. , Dasi, L. P. , Kim, H. S. , Choi, J. , Leo, H. , and Yoganathan, A. P. , 2008, “ Structural Simulations of Prosthetic Tri-Leaflet Aortic Heart Valves,” J. Biomech., 41(7), pp. 1510–1519. [CrossRef] [PubMed]
Kim, H. , Lu, J. , Sacks, M. S. , and Chandran, K. B. , 2008, “ Dynamic Simulation of Bioprosthetic Heart Valves Using a Stress Resultant Shell Model,” Ann. Biomed. Eng., 36(2), pp. 262–275. [CrossRef] [PubMed]
Shi, Y. , Yao, J. , Young, J. M. , Fee, J. A. , Perucchio, R. , and Taber, L. A. , 2014, “ Bending and Twisting the Embryonic Heart: A Computational Model for C-Looping Based on Realistic Geometry,” Front. Physiol., 5, p. 297. [CrossRef] [PubMed]
Sun, W. , Li, K. , and Sirois, E. , 2010, “ Simulated Elliptical Bioprosthetic Valve Deformation: Implications for Asymmetric Transcatheter Valve Deployment,” J. Biomech., 43(16), pp. 3085–3090. [CrossRef] [PubMed]
Venkatasubramaniam, A. , Fagan, M. , Mehta, T. , Mylankal, K. , Ray, B. , Kuhan, G. , Chetter, I. C. , and McCollum, P. T. , 2004, “ A Comparative Study of Aortic Wall Stress Using Finite Element Analysis for Ruptured and Non-Ruptured Abdominal Aortic Aneurysms,” Eur. J. Vasc. Endovasc. Surg., 28(2), pp. 168–176. [PubMed]
ABAQUS, 2011, “ ABAQUS/Standard: User's Manual,” Dassault Systèmes Simulia, Johnston, RI.
Simo, J. , and Hughes, T. J. R. , 1998, Computational Inelasticity, Springer, New York.
Prot, V. , Skallerud, B. , and Holzapfel, G. , 2007, “ Transversely Isotropic Membrane Shells With Application to Mitral Valve Mechanics. Constitutive Modelling and Finite Element Implementation,” Int. J. Numer. Methods Eng., 71(8), pp. 987–1008. [CrossRef]
ABAQUS, 2011, “ ABAQUS Theory Manual, ABAQUS 611 Documentation,” Dassault Systèmes Simulia, Johnston, RI.
Miehe, C. , 1996, “ Numerical Computation of Algorithmic (Consistent) Tangent Moduli in Large-Strain Computational Inelasticity,” Comput. Methods Appl. Mech. Eng., 134(3–4), pp. 223–240. [CrossRef]
Sun, W. , Chaikof, E. L. , and Levenston, M. E. , 2008, “ Numerical Approximation of Tangent Moduli for Finite Element Implementations of Nonlinear Hyperelastic Material Models,” ASME J. Biomech. Eng., 130(6), p. 061003. [CrossRef]
Liu, H. , and Sun, W. , 2015, “ Computational Efficiency of Numerical Approximations of Tangent Moduli for Finite Element Implementation of a Fiber-Reinforced Hyperelastic Material Model,” Comput. Methods Biomech. Biomed. Eng., 19(11), pp. 1–10.
Tanaka, M. , and Fujikawa, M. , 2011, “ Numerical Approximation of Consistent Tangent Moduli Using Complex-Step Derivative and Its Application to Finite Deformation Problems,” Trans. Jpn. Soc. Mech. Eng., Ser. A, 77(773), pp. 27–38. [CrossRef]
Holzapfel, G. A. , 2000, Nonlinear Solid Mechanics: A Continuum Approach for Engineering, Wiley, Chichester, UK.
Davis, F. M. , Luo, Y. , Avril, S. , Duprey, A. , and Lu, J. , 2015, “ Pointwise Characterization of the Elastic Properties of Planar Soft Tissues: Application to Ascending Thoracic Aneurysms,” Biomech. Model. Mechanobiol., 14(5), pp. 967–978. [CrossRef] [PubMed]
Schiavone, A. , and Zhao, L. G. , 2015, “ A Study of Balloon Type, System Constraint and Artery Constitutive Model Used in Finite Element Simulation of Stent Deployment,” Mech. Adv. Mater. Modern Processes, 1(1), pp. 1–15.
Smoljkić, M. , Vander Sloten, J. , Segers, P. , and Famaey, N. , 2015, “ Non-Invasive, Energy-Based Assessment of Patient-Specific Material Properties of Arterial Tissue,” Biomech. Model. Mechanobiol., 14(5), pp. 1045–1056. [CrossRef] [PubMed]
Sun, W. , Chaikof, E. L. , and Levenston, M. E. , 2008, “ Development and Finite Element Implementation of a Nearly Incompressible Structural Constitutive Model for Artery Substitute Design,” ASME Paper No. SBC2008-193164.
Sun, W. , and Sacks, M. S. , 2005, “ Finite Element Implementation of a Generalized Fung-Elastic Constitutive Model for Planar Soft Tissues,” Biomech. Model. Mechanobiol., 4(2), pp. 190–199. [CrossRef] [PubMed]
Sacks, M. , 1999, “ A Method for Planar Biaxial Mechanical Testing That Includes In-Plane Shear,” ASME J. Biomech. Eng., 121(5), pp. 551–555. [CrossRef]
Sacks, M. S. , Smith, D. B. , and Hiester, E. D. , 1997, “ A Small Angle Light Scattering Device for Planar Connective Tissue Microstructural Analysis,” Ann. Biomed. Eng., 25(4), pp. 678–689. [CrossRef] [PubMed]
Tanaka, M. , Fujikawa, M. , Balzani, D. , and Schröder, J. , 2014, “ Robust Numerical Calculation of Tangent Moduli at Finite Strains Based on Complex-Step Derivative Approximation and Its Application to Localization Analysis,” Comput. Methods Appl. Mech. Eng., 269, pp. 454–470. [CrossRef]
Zulliger, M. A. , Fridez, P. , Hayashi, K. , and Stergiopulos, N. , 2004, “ A Strain-Energy Function for Arteries Accounting for Wall Composition and Structure,”J. Biomech., 37(7), pp. 989–1000. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

FE meshes used in the biaxial test simulation: the arrows denote the loading directions with their round roots showing the sites of loading, the box enclosed region in the center represents the area delineated by the four markers in the biaxial test (a), Von Mises stress contour in the deformed configuration (b), and curves of Cauchy stress versus Green strain resulted from experiment (dashed line) and numerical simulation (solid line) (c)

Grahic Jump Location
Fig. 2

Relative error of the Jacobian approximation at different logarithmic perturbations (a), the CPU time (b), and the iteration number (c) at different logarithmic perturbations using themodified HGO model in the biaxial testing simulation. Dashed lines (—) represent the level achieved using UHYPER subroutine.

Grahic Jump Location
Fig. 3

FE meshes used in the vessel inflation simulation, before (a) and after (b) deformation; loading pressure versus the outermost radius of the deformed vessel, observed in experiment (dotted line) and numerical simulations. Material models with various mean fiber directions while the other parameters are kept the same as baseline values (c); material models with various standard deviations while the other parameters are kept the same as baseline values (d); and the iteration number at different logarithmic perturbations using the distributed fiber model in the rat carotid artery simulation (e).



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In