A coupled fiber-matrix model demonstrates highly inhomogeneous microstructural interactions in soft tissues under tensile load

Lijuan Zhang, Spencer P. Lake, Victor Lai, Catalin R. Picu, Victor H Barocas, Mark S. Shephard

Research output: Contribution to journalArticle

23 Citations (Scopus)

Abstract

A soft tissue's macroscopic behavior is largely determined by its microstructural components (often a collagen fiber network surrounded by a nonfibrillar matrix (NFM)). In the present study, a coupled fiber-matrix model was developed to fully quantify the internal stress field within such a tissue and to explore interactions between the collagen fiber network and nonfibrillar matrix (NFM). Voronoi tessellations (representing collagen networks) were embedded in a continuous three-dimensional NFM. Fibers were represented as one-dimensional nonlinear springs and the NFM, meshed via tetrahedra, was modeled as a compressible neo-Hookean solid. Multidimensional finite element modeling was employed in order to couple the two tissue components and uniaxial tension was applied to the composite representative volume element (RVE). In terms of the overall RVE response (average stress, fiber orientation, and Poisson's ratio), the coupled fiber-matrix model yielded results consistent with those obtained using a previously developed parallel model based upon superposition. The detailed stress field in the composite RVE demonstrated the high degree of inhomogeneity in NFM mechanics, which cannot be addressed by a parallel model. Distributions of maximum/minimum principal stresses in the NFM showed a transition from fiber-dominated to matrix-dominated behavior as the matrix shear modulus increased. The matrix-dominated behavior also included a shift in the fiber kinematics toward the affine limit. We conclude that if only gross averaged parameters are of interest, parallel-type models are suitable. If, however, one is concerned with phenomena, such as individual cell-fiber interactions or tissue failure that could be altered by local variations in the stress field, then the detailed model is necessary in spite of its higher computational cost.

Original languageEnglish (US)
Article number011008
JournalJournal of Biomechanical Engineering
Volume135
Issue number1
DOIs
StatePublished - Jan 8 2013

Fingerprint

Optical fiber coupling
Tissue
Collagen
Fibers
Stress Fibers
Response Elements
Mechanics
Biomechanical Phenomena
Cell Communication
Composite materials
Poisson ratio
Fiber reinforced materials
Costs and Cost Analysis
Residual stresses
Kinematics
Elastic moduli
Costs

Keywords

  • collagen network
  • computational biomechanics
  • fiber-matrix interaction
  • non-fibrillar matrix
  • nonaffine behavior
  • soft tissue

Cite this

A coupled fiber-matrix model demonstrates highly inhomogeneous microstructural interactions in soft tissues under tensile load. / Zhang, Lijuan; Lake, Spencer P.; Lai, Victor; Picu, Catalin R.; Barocas, Victor H; Shephard, Mark S.

In: Journal of Biomechanical Engineering, Vol. 135, No. 1, 011008, 08.01.2013.

Research output: Contribution to journalArticle

@article{ceab1938d6144731b40223ed7bd4e26b,
title = "A coupled fiber-matrix model demonstrates highly inhomogeneous microstructural interactions in soft tissues under tensile load",
abstract = "A soft tissue's macroscopic behavior is largely determined by its microstructural components (often a collagen fiber network surrounded by a nonfibrillar matrix (NFM)). In the present study, a coupled fiber-matrix model was developed to fully quantify the internal stress field within such a tissue and to explore interactions between the collagen fiber network and nonfibrillar matrix (NFM). Voronoi tessellations (representing collagen networks) were embedded in a continuous three-dimensional NFM. Fibers were represented as one-dimensional nonlinear springs and the NFM, meshed via tetrahedra, was modeled as a compressible neo-Hookean solid. Multidimensional finite element modeling was employed in order to couple the two tissue components and uniaxial tension was applied to the composite representative volume element (RVE). In terms of the overall RVE response (average stress, fiber orientation, and Poisson's ratio), the coupled fiber-matrix model yielded results consistent with those obtained using a previously developed parallel model based upon superposition. The detailed stress field in the composite RVE demonstrated the high degree of inhomogeneity in NFM mechanics, which cannot be addressed by a parallel model. Distributions of maximum/minimum principal stresses in the NFM showed a transition from fiber-dominated to matrix-dominated behavior as the matrix shear modulus increased. The matrix-dominated behavior also included a shift in the fiber kinematics toward the affine limit. We conclude that if only gross averaged parameters are of interest, parallel-type models are suitable. If, however, one is concerned with phenomena, such as individual cell-fiber interactions or tissue failure that could be altered by local variations in the stress field, then the detailed model is necessary in spite of its higher computational cost.",
keywords = "collagen network, computational biomechanics, fiber-matrix interaction, non-fibrillar matrix, nonaffine behavior, soft tissue",
author = "Lijuan Zhang and Lake, {Spencer P.} and Victor Lai and Picu, {Catalin R.} and Barocas, {Victor H} and Shephard, {Mark S.}",
year = "2013",
month = "1",
day = "8",
doi = "10.1115/1.4023136",
language = "English (US)",
volume = "135",
journal = "Journal of Biomechanical Engineering",
issn = "0148-0731",
publisher = "American Society of Mechanical Engineers(ASME)",
number = "1",

}

TY - JOUR

T1 - A coupled fiber-matrix model demonstrates highly inhomogeneous microstructural interactions in soft tissues under tensile load

AU - Zhang, Lijuan

AU - Lake, Spencer P.

AU - Lai, Victor

AU - Picu, Catalin R.

AU - Barocas, Victor H

AU - Shephard, Mark S.

PY - 2013/1/8

Y1 - 2013/1/8

N2 - A soft tissue's macroscopic behavior is largely determined by its microstructural components (often a collagen fiber network surrounded by a nonfibrillar matrix (NFM)). In the present study, a coupled fiber-matrix model was developed to fully quantify the internal stress field within such a tissue and to explore interactions between the collagen fiber network and nonfibrillar matrix (NFM). Voronoi tessellations (representing collagen networks) were embedded in a continuous three-dimensional NFM. Fibers were represented as one-dimensional nonlinear springs and the NFM, meshed via tetrahedra, was modeled as a compressible neo-Hookean solid. Multidimensional finite element modeling was employed in order to couple the two tissue components and uniaxial tension was applied to the composite representative volume element (RVE). In terms of the overall RVE response (average stress, fiber orientation, and Poisson's ratio), the coupled fiber-matrix model yielded results consistent with those obtained using a previously developed parallel model based upon superposition. The detailed stress field in the composite RVE demonstrated the high degree of inhomogeneity in NFM mechanics, which cannot be addressed by a parallel model. Distributions of maximum/minimum principal stresses in the NFM showed a transition from fiber-dominated to matrix-dominated behavior as the matrix shear modulus increased. The matrix-dominated behavior also included a shift in the fiber kinematics toward the affine limit. We conclude that if only gross averaged parameters are of interest, parallel-type models are suitable. If, however, one is concerned with phenomena, such as individual cell-fiber interactions or tissue failure that could be altered by local variations in the stress field, then the detailed model is necessary in spite of its higher computational cost.

AB - A soft tissue's macroscopic behavior is largely determined by its microstructural components (often a collagen fiber network surrounded by a nonfibrillar matrix (NFM)). In the present study, a coupled fiber-matrix model was developed to fully quantify the internal stress field within such a tissue and to explore interactions between the collagen fiber network and nonfibrillar matrix (NFM). Voronoi tessellations (representing collagen networks) were embedded in a continuous three-dimensional NFM. Fibers were represented as one-dimensional nonlinear springs and the NFM, meshed via tetrahedra, was modeled as a compressible neo-Hookean solid. Multidimensional finite element modeling was employed in order to couple the two tissue components and uniaxial tension was applied to the composite representative volume element (RVE). In terms of the overall RVE response (average stress, fiber orientation, and Poisson's ratio), the coupled fiber-matrix model yielded results consistent with those obtained using a previously developed parallel model based upon superposition. The detailed stress field in the composite RVE demonstrated the high degree of inhomogeneity in NFM mechanics, which cannot be addressed by a parallel model. Distributions of maximum/minimum principal stresses in the NFM showed a transition from fiber-dominated to matrix-dominated behavior as the matrix shear modulus increased. The matrix-dominated behavior also included a shift in the fiber kinematics toward the affine limit. We conclude that if only gross averaged parameters are of interest, parallel-type models are suitable. If, however, one is concerned with phenomena, such as individual cell-fiber interactions or tissue failure that could be altered by local variations in the stress field, then the detailed model is necessary in spite of its higher computational cost.

KW - collagen network

KW - computational biomechanics

KW - fiber-matrix interaction

KW - non-fibrillar matrix

KW - nonaffine behavior

KW - soft tissue

UR - http://www.scopus.com/inward/record.url?scp=84871900713&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84871900713&partnerID=8YFLogxK

U2 - 10.1115/1.4023136

DO - 10.1115/1.4023136

M3 - Article

VL - 135

JO - Journal of Biomechanical Engineering

JF - Journal of Biomechanical Engineering

SN - 0148-0731

IS - 1

M1 - 011008

ER -