Biphasic Finite Element Modeling of Hydrated Soft Tissue Contact Using an Augmented Lagrangian Method

Guo, Haijie; Spilker, Robert L.

doi:10.1115/1.4005378

Cited by 18 publications

(28 citation statements)

References 20 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…An excellent agreement between the predictions of ITERATIVE, FEBio and the literature (Guo and Spilker 2011) was found along the contact boundary of the glenoid cartilage, regarding both the pore pressure distribution (Figure 11(a)) and the total stress (Figure 11(b)). A higher sensitivity to the mesh size was found in ITERATIVE compared to FEBio.…”

Section: Unconfined Compression Test -Glenohumeral Jointmentioning

confidence: 58%

“…Pore pressure (a) and total stress (b) distributions along the glenoid contact surface predicted for the unconfined compression test of the glenohumeral joint with the three algorithms STANDARD, FEBio and ITERATIVE. Results of another algorithm reported in the literature (Guo and Spilker 2011) are also shown.…”

Section: Discussionmentioning

confidence: 89%

“…The glenohumeral joint contact of the human shoulder (Figure 3(c)) has been simulated with data from a previous work (Guo and Spilker 2011). Since in the reference paper a substantial equivalence between axisymmetric and threedimensional finite element solutions was found, the model …”

Section: Unconfined Compression Test -Glenohumeral Jointmentioning

confidence: 99%

See 2 more Smart Citations

Comparison of various contact algorithms for poroelastic tissues

Galbusera¹,

Bashkuev²,

Wilke³

et al. 2012

Computer Methods in Biomechanics and Biomedical Engineering

View full text Add to dashboard Cite

Capabilities of the commercial finite element package ABAQUS in simulating frictionless contact between two saturated porous structures were evaluated and compared with those of an open source code, FEBio. In ABAQUS, both the default contact implementation and another algorithm based on an iterative approach requiring script programming were considered. Test simulations included a patch test of two cylindrical slabs in a gapless contact and confined compression conditions; a confined compression test of a porous cylindrical slab with a spherical porous indenter; and finally two unconfined compression tests of soft tissues mimicking diarthrodial joints. The patch test showed almost identical results for all algorithms. On the contrary, the confined and unconfined compression tests demonstrated large differences related to distinct physical and boundary conditions considered in each of the three contact algorithms investigated in this study. In general, contact with non-uniform gaps between fluid-filled porous structures could be effectively simulated with either ABAQUS or FEBio. The user should be aware of the parameter definitions, assumptions and limitations in each case, and take into consideration the physics and boundary conditions of the problem of interest when searching for the most appropriate model.

show abstract

Section: Unconfined Compression Test -Glenohumeral Jointmentioning

confidence: 58%

Section: Discussionmentioning

confidence: 89%

See 1 more Smart Citation

Comparison of various contact algorithms for poroelastic tissues

Galbusera¹,

Bashkuev²,

Wilke³

et al. 2012

Computer Methods in Biomechanics and Biomedical Engineering

View full text Add to dashboard Cite

show abstract

“…Solid mechanics in the Structural Mechanics Module and Darcy’s Law in the Earth Science Module were used (Guo et al, 2013; Guo et al, 2012; Guo and Spilker, 2011; Guo and Spilker, 2014). The user defined strain energy density function in COMSOL was used to input the strain energy density function (Eqn.…”

Section: Methodsmentioning

confidence: 99%

A biphasic finite element study on the role of the articular cartilage superficial zone in confined compression

Guo

Maher

Torzilli

2015

Journal of Biomechanics

Self Cite

View full text Add to dashboard Cite

The aim of this study was to investigate the role of the superficial zone on the mechanical behavior of articular cartilage. Confined compression of articular cartilage was modeled using a biphasic finite element analysis to calculate the one-dimensional deformation of the extracellular matrix (ECM) and movement of the interstitial fluid through the ECM and articular surface. The articular cartilage was modeled as an inhomogeneous, nonlinear hyperelastic biphasic material with depth and strain-dependent material properties. Two loading conditions were simulated, one where the superficial zone was loaded with a porous platen (normal test) and the other where the deep zone was loaded with the porous platen (upside down test). Compressing the intact articular cartilage with 0.2 MPa stress reduced the surface permeability by 88%. Removing the superficial zone increased the rate of change for all mechanical parameters and decreased the fluid support ratio of the tissue, resulting in increased tissue deformation. Apparent permeability linearly increased after superficial removal in the normal test, yet it did not change in the upside down test. Orientation of the specimen affected the time-dependent biomechanical behavior of the articular cartilage, but not equilibrium behavior. The two tests with different specimen orientations resulted in very different apparent permeabilities, suggesting that in an experimental study which quantifies material properties of an inhomogeneous material, the specimen orientation should be stated along with the permeability result. The current study provides new insights into the role of the superficial zone on mechanical behavior of the articular cartilage.

show abstract

“…All other geometric values and material properties are considered as variable inputs with values chosen to span the ranges reported in the literature for normal knees and listed in Table 1. All bFEMs were solved using COMSOL Multiphysics (COSMOL Inc., Burlington, MA) on a Mac workstation with two 3.06 GHz 6-Core Intel Xeon processors and 24 GB 1333 MHz DDR3 Memory [19][20][21][22].…”

Section: The Finite Element Simulatormentioning

confidence: 99%

Using a Statistically Calibrated Biphasic Finite Element Model of the Human Knee Joint to Identify Robust Designs for a Meniscal Substitute

Leatherman

Guo

Gilbert

et al. 2014

Journal of Biomechanical Engineering

Self Cite

View full text Add to dashboard Cite

This paper describes a methodology for selecting a set of biomechanical engineering design variables to optimize the performance of an engineered meniscal substitute when implanted in a population of subjects whose characteristics can be specified stochastically. For the meniscal design problem where engineering variables include aspects of meniscal geometry and meniscal material properties, this method shows that meniscal designs having simultaneously large radial modulus and large circumferential modulus provide both low mean peak contact stress and small variability in peak contact stress when used in the specified subject population. The method also shows that the mean peak contact stress is relatively insensitive to meniscal permeability, so the permeability used in the manufacture of a meniscal substitute can be selected on the basis of manufacturing ease or cost. This is a multiple objective problem with the mean peak contact stress over the population of subjects and its variability both desired to be small. The problem is solved by using a predictor of the mean peak contact stress across the tibial plateau that was developed from experimentally measured peak contact stresses from two modalities. The first experimental modality provided computed peak contact stresses using a finite element computational simulator of the dynamic tibial contact stress during axial dynamic loading. A small number of meniscal designs with specified subject environmental inputs were selected to make computational runs and to provide training data for the predictor developed below. The second experimental modality consisted of measured peak contact stress from a set of cadaver knees. The cadaver measurements were used to bias-correct and calibrate the simulator output. Because the finite element simulator is expensive to evaluate, a rapidly computable (calibrated) Kriging predictor was used to explore extensively the contact stresses for a wide range of meniscal engineering inputs and subject variables. The predicted values were used to determine the Pareto optimal set of engineering inputs to minimize peak contact stresses in the targeted population of subjects.

show abstract

Biphasic Finite Element Modeling of Hydrated Soft Tissue Contact Using an Augmented Lagrangian Method

Cited by 18 publications

References 20 publications

Comparison of various contact algorithms for poroelastic tissues

Comparison of various contact algorithms for poroelastic tissues

A biphasic finite element study on the role of the articular cartilage superficial zone in confined compression

Using a Statistically Calibrated Biphasic Finite Element Model of the Human Knee Joint to Identify Robust Designs for a Meniscal Substitute

Contact Info

Product

Resources

About