A Computationally Efficient Optimization Kernel for Material Parameter Estimation Procedures

[+] Author and Article Information
H. Schmid1

Bioengineering Institute, University of Auckland, Private Bag 92019, Auckland, 1001 New Zealandh.schmid@auckland.ac.nz

M. P. Nash, A. A. Young, O. Röhrle, P. J. Hunter

Bioengineering Institute, University of Auckland, Private Bag 92019, Auckland, 1001 New Zealand


Corresponding author.

J Biomech Eng 129(2), 279-283 (Sep 21, 2006) (5 pages) doi:10.1115/1.2540860 History: Received September 19, 2005; Revised September 21, 2006

Estimating material parameters is an important part in the study of soft tissue mechanics. Computational time can easily run to days, especially when all available experimental data are taken into account. The material parameter estimation procedure is examplified on a set of homogeneous simple shear experiments to estimate the orthotropic constitutive parameters of myocardium. The modification consists of changing the traditional least-squares approach to a weighted least-squares. This objective function resembles a L2-norm type integral which is approximated using Gaussian quadrature. This reduces the computational time of the material parameter estimation by two orders of magnitude.

Copyright © 2007 by American Society of Mechanical Engineers
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Figure 1

A sketch of all six possible shear modes for myocardium. The first letter indicates the normal vector of the face that is shifted, and the second letter indicates the direction in which it is shifted. The condensed lines indicate the fibers and the other lines indicate the sheets. With courtesy of American Journal of Physiology, Dokos (see Ref. 8).

Grahic Jump Location
Figure 2

Convergence analysis of the objective function ∣ΔΩrelGLS∣ and the difference vector ∣ΔϑrelGLSk∣ plotted against the number of Gauss points G

Grahic Jump Location
Figure 3

The figure depicts the experimental (dotted) and fitted force-displacement curves (solid) of the Costa law for all six modes for one experiment. The overall error ΩrelGLS is 1.5%. Note the different scales on each graph. The abscissa shows the displacement in mm, whereas the ordinate shows the top face force in mN.



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