Research Papers

A Coupled Fiber-Matrix Model Demonstrates Highly Inhomogeneous Microstructural Interactions in Soft Tissues Under Tensile Load

[+] Author and Article Information
Lijuan Zhang

Scientific Computation Research Center,
Rensselaer Polytechnic Institute,
Low Center for Industrial Innovation, CII-4011,
110 8th Street,
Troy, NY 12180

Spencer P. Lake, Victor H. Barocas

Department of Biomedical Engineering,
University of Minnesota,
7-105 Nils Hasselmo Hall,
312 Church Street SE,
Minneapolis, MN 55455

Victor K. Lai

Department of Chemical Engineering
and Materials Science,
University of Minnesota,
421 Washington Ave SE,
Minneapolis, MN 55455

Catalin R. Picu

Scientific Computation Research Center,
Rensselaer Polytechnic Institute,
Low Center for Industrial Innovation,
CII-4011, 110 8th Street,
Troy, NY 12180;
Department of Mechanical,
Aerospace and Nuclear Engineering,
Rensselaer Polytechnic Institute,
Jonsson Engineering Center,
Rm. 2049, 110 8th Street,
Troy, NY 12180

Mark S. Shephard

Scientific Computation Research Center,
Rensselaer Polytechnic Institute,
Low Center for Industrial Innovation,
CII-4011, 110 8th Street,
Troy, NY 12180
e-mail: shephard@rpi.edu

1Corresponding author.

Contributed by the Bioengineering Division of ASME for publication in the Journal of Biomechanical Engineering. Manuscript received July 26, 2012; final manuscript received November 16, 2012; accepted manuscript posted December 8, 2012; published online December 27, 2012. Assoc. Editor: James C. Iatridis.

J Biomech Eng 135(1), 011008 (Dec 27, 2012) (9 pages) Paper No: BIO-12-1323; doi: 10.1115/1.4023136 History: Received July 26, 2012; Revised November 16, 2012; Accepted December 08, 2012

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.

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


Fung, Y. C., 1967, “Elasticity of Soft Tissues in Simple Elongation,” Am. J. Physiol., 213(6), pp. 1532–1544. [PubMed]
Humphrey, J. D., 1995, “Mechanics of the Arterial Wall: Review and Directions,” Crit. Rev. Biomed. Eng., 23(1–2), pp. 1–162. [PubMed]
Horgan, C. O. and Saccomandi, G., 2003, “A Description of Arterial Wall Mechanics Using Limiting Chain Extensibility Constitutive Models,” Biomech. Model. Mechanobiol., 1(4), pp. 251–266. [CrossRef] [PubMed]
Lanir, Y., 1979, “A Structural Theory for the Homogeneous Biaxial Stress-Strain Relationships in Flat Collagenous Tissues,” J. Biomech., 12(6), pp. 423–436. [CrossRef] [PubMed]
Decraemer, W. F., Maes, M. A., and Vanhuyse, V. J., 1980, “An Elastic Stress-Strain Relation for Soft Biological Tissues Based on a Structural Model,” J. Biomech., 13(6), pp. 463–468. [CrossRef] [PubMed]
Kwan, M. K. and Woo, S. L., 1989, “A Structural Model to Describe the Nonlinear Stress-Strain Behavior for Parallel-Fibered Collagenous Tissues,” ASME J. Biomech. Eng., 111(4), pp. 361–363. [CrossRef]
Cortes, D. H., Lake, S. P., Kadlowec, J. A., Soslowsky, L. J., and Elliott, D. M., 2010, “Characterizing the Mechanical Contribution of Fiber Angular Distribution in Connective Tissue: Comparison of Two Modeling Approaches,” Biomech. Model. Mechanobiol., 9(5), pp. 651–658. [CrossRef] [PubMed]
Wagner, H. P. and Humphrey, J. D., 2011, “Differential Passive and Active Biaxial Mechanical Behaviors of Muscular and Elastic Arteries: Basilar Versus Common Carotid,” ASME J. Biomech. Eng., 133(5), p. 051009. [CrossRef]
Hollander, Y., Durban, D., Lu, X., Kassab, G. S., and Lanir, Y., 2011, “Experimentally Validated Microstructural 3D Constitutive Model of Coronary Arterial Media,” ASME J. Biomech. Eng., 133(3), p. 031007. [CrossRef]
Kao, P. H., Lammers, S. R., Tian, L., Hunter, K., Stenmark, K. R., Shandas, R., and Qi, H. J., 2011, “A Microstructurally Driven Model for Pulmonary Artery Tissue,” ASME J. Biomech. Eng., 133(5), p. 051002. [CrossRef]
Soares, A. L., Stekelenburg, M., and Baaijens, F. P., 2011, “Remodeling of the Collagen Fiber Architecture Due to Compaction in Small Vessels Under Tissue Engineered Conditions,” ASME J. Biomech. Eng., 133(7), p. 071002. [CrossRef]
Szczesny, S. E., Peloquin, J. M., Cortes, D. H., Kadlowec, J. A., Soslowsky, L. J., and Elliott, D. M., 2012, “Biaxial Tensile Testing and Constitutive Modeling of Human Supraspinatus Tendon,” ASME J. Biomech. Eng., 134(2), p. 021004. [CrossRef]
Holzapfel, G. A., Gasser, T. C., and Ogden, R. W., 2000, “A New Constitutive Framework for Arterial Wall Mechanics and a Comparative Study of Material Models,” J. Elast., 61, pp. 1–48. [CrossRef]
Driessen, N. J., Bouten, C. V., and Baaijens, F. P., 2005, “A Structural Constitutive Model for Collagenous Cardiovascular Tissues Incorporating the Angular Fiber Distribution,” ASME J. Biomech. Eng., 127(3), pp. 494–503. [CrossRef]
Tang, H., Buehler, M. J., and Moran, B., 2009, “A Constitutive Model of Soft Tissue: From Nanoscale Collagen to Tissue Continuum,” Ann. Biomed. Eng., 37(6), pp. 1117–1130. [CrossRef] [PubMed]
Nagel, T. and Kelly, D. J., 2012, “Remodelling of Collagen Fibre Transition Stretch and Angular Distribution in Soft Biological Tissues and Cell-Seeded Hydrogels,” Biomech. Model. Mechanobiol., 11(3–4), pp. 325–339. [CrossRef] [PubMed]
Quapp, K. M. and Weiss, J. A., 1998, “Material Characterization of Human Medial Collateral Ligament,” ASME J. Biomech. Eng., 120(6), pp. 757–763. [CrossRef]
Guerin, H. L. and Elliott, D. M., 2007, “Quantifying the Contributions of Structure to Annulus Fibrosus Mechanical Function Using a Nonlinear, Anisotropic, Hyperelastic Model,” J. Orthop. Res., 25(4), pp. 508–516. [CrossRef] [PubMed]
Abraham, A. C., Moyer, J. T., Villegas, D. F., Odegard, G. M., and Haut Donahue, T. L., 2011, “Hyperelastic Properties of Human Meniscal Attachments,” J. Biomech., 44(3), pp. 413–418. [CrossRef] [PubMed]
Cortes, D. H. and Elliott, D. M., 2012, “Extra-Fibrillar Matrix Mechanics of Annulus Fibrosus in Tension and Compression,” Biomech. Model. Mechanobiol., 11(6), pp. 781–790. [CrossRef] [PubMed]
Wagner, D. R. and Lotz, J. C., 2004, “Theoretical Model and Experimental Results for the Nonlinear Elastic Behavior of Human Annulus Fibrosus,” J. Orthop. Res., 22(4), pp. 901–909. [CrossRef] [PubMed]
Peng, X. Q., Guo, Z. Y., and Moran, B., 2006, “An Anisotropic Hyperelastic Constitutive Model With Fiber-Matrix Shear Interaction for the Human Annulus Fibrosus,” ASME J. Appl. Mech., 73, pp. 815–824. [CrossRef]
O'Connell, G. D., Guerin, H. L., and Elliott, D. M., 2009, “Theoretical and Uniaxial Experimental Evaluation of Human Annulus Fibrosus Degeneration,” ASME J. Biomech. Eng., 131(11), p. 111007. [CrossRef]
Lake, S. P., Hadi, M. F., Lai, V. K., and Barocas, V. H., 2012, “Mechanics of a Fiber Network Within a Non-Fibrillar Matrix: Model and Comparison With Collagen-Agarose Co-gels,” Ann. Biomed. Eng., 40(10), pp. 2111–2121. [CrossRef] [PubMed]
Lake, S. P. and Barocas, V. H., 2011, “Mechanical and Structural Contribution of Non-Fibrillar Matrix in Uniaxial Tension: A Collagen-Agarose Co-Gel Model,” Ann. Biomed. Eng., 39(7), pp. 1891–1903. [CrossRef] [PubMed]
Chandran, P. L. and Barocas, V. H., 2006, “Affine Versus Non-Affine Fibril Kinematics in Collagen Networks: Theoretical Studies of Network Behavior,” ASME J. Biomech. Eng., 128(2), pp. 259–270. [CrossRef]
Stylianopoulos, T. and Barocas, V. H., 2007, “Multiscale, Structure-Based Modeling for the Elastic Mechanical Behavior of Arterial Walls,” ASME J. Biomech. Eng., 129(4), pp. 611–618. [CrossRef]
Sander, E. A., Stylianopoulos, T., Tranquillo, R. T., and Barocas, V. H., 2009, “Image-Based Multiscale Modeling Predicts Tissue-Level and Network-Level Fiber Reorganization in Stretched Cell-Compacted Collagen Gels,” Proc. Natl. Acad. Sci. U.S.A., 106(42), pp. 17675–17680. [CrossRef] [PubMed]
Kang, J., Steward, R. L., Kim, Y., Schwartz, R. S., LeDuc, P. R., and Puskar, K. M., 2011, “Response of an Actin Filament Network Model Under Cyclic Stretching Through a Coarse Grained Monte Carlo Approach,” J. Theor. Biol., 274(1), pp. 109–119. [CrossRef] [PubMed]
Nachtrab, S., Kapfer, S. C., Arns, C. H., Madadi, M., Mecke, K., and Schroder-Turk, G. E., 2011, “Morphology and Linear-Elastic Moduli of Random Network Solids,” Adv. Mater., 23(22–23), pp. 2633–2637. [CrossRef] [PubMed]
Stein, A. M., Vader, D. A., Weitz, D. A., and Sander, L. M., 2011, “The Micromechanics of Three-Dimensional Collagen-I Gels,” Complexity, 16(4), pp. 22–28. [CrossRef]
Billiar, K. L. and Sacks, M. S., 2000, “Biaxial Mechanical Properties of the Native and Glutaraldehyde-Treated Aortic Valve Cusp—Part II: A Structural Constitutive Model,” ASME J. Biomech. Eng., 122(4), pp. 327–335. [CrossRef]
Holzapfel, G. A., 2000, Nonlinear Solid Mechanics: A Contiuum Approach for Engineering, John Wiley and Sons, New York.
Nagel, T. and Kelly, D. J., 2010, “The Influence of Fiber Orientation on the Equilibrium Properties of Neutral and Charged Biphasic Tissues,” ASME J. Biomech. Eng., 132(11), p. 114506. [CrossRef]
Benkherourou, M., Rochas, C., Tracqui, P., Tranqui, L., and Gumery, P. Y., 1999, “Standardization of a Method for Characterizing Low-Concentration Biogels: Elastic Properties of Low-Concentration Agarose Gels,” ASME J. Biomech. Eng., 121(2), pp. 184–187. [CrossRef]
Lake, S. P., Hald, E. S., and Barocas, V. H., 2011, “Collagen-Agarose Co-Gels as a Model for Collagen-Matrix Interaction in Soft Tissues Subjected to Indentation,” J. Biomed. Mater. Res. Part A, 99(4), pp. 507–515. [CrossRef]
Chandran, P. L. and Barocas, V. H., 2007, “Deterministic Material-Based Averaging Theory Model of Collagen Gel Micromechanics,” ASME J. Biomech. Eng., 129(2), pp. 137–147. [CrossRef]
Stylianopoulos, T. and Barocas, V. H., 2007, “Volume-Averaging Theory for the Study of the Mechanics of Collagen Networks,” Comput. Methods Biomech. Biomed. Eng., 196(31–32), pp. 2981–2990. [CrossRef]
Weiler, K. J.1988, “The Radial-Edge Structure: A Topological Representation for Non-Manifold Geometric Boundary Representations,” Geometric Modeling for CAD Applications, M. J.Wozny, H. W.McLaughlin, and J. L.Encarnacao, eds., North Holland, Amsterdam, pp. 3–36.
Simmetrix web page, http://www.simmetrix.com
Shephard, M. S., 2000, “Meshing Environment for Geometry-Based Analysis,” Int. J. Numer. Methods Eng., 47(1–3), pp. 169–190. [CrossRef]
Franck, C., Maskarinec, S. A., Tirrell, D. A., and Ravichandran, G., 2011, “Three-Dimensional Traction Force Microscopy: A New Tool for Quantifying Cell-Matrix Interactions,” PLoS ONE, 6(3), p. e17833. [CrossRef] [PubMed]
Hadi, M. F., Sander, E. A., and Barocas, V. H., 2012, “Multiscale Model Predicts Tissue-Level Failure From Collagen Fiber-Level Damage,” ASME J. Biomech. Eng., 134(9), p. 091005. [CrossRef]
Hadi, M. F., Sander, E. A., Ruberti, J. W., and Barocas, V. H., 2012, “Simulated Remodeling of Loaded Collagen Networks via Strain-Dependent Enzymatic Degradation and Constant-Rate Fiber Growth,” Mech. Mater., 44, pp. 72–82. [CrossRef] [PubMed]
Pence, T. J., Monroe, R. J., and Wright, N. T., 2008, “On the Computation of Stress in Affine Versus Nonaffine Fibril Kinematics Within Planar Collagen Network Models,” ASME J. Biomech. Eng., 130(4), p. 041009. [CrossRef]
Sopakayang, R., De Vita, R., Kwansa, A., and Freeman, J. W., 2012, “Elastic and Viscoelastic Properties of a Type I Collagen Fiber,” J. Theor. Biol., 293, pp. 197–205. [CrossRef] [PubMed]
Silver, F. H., Ebrahimi, A., and Snowhill, P. B., 2002, “Viscoelastic Properties of Self-Assembled Type I Collagen Fibers: Molecular Basis of Elastic and Viscous Behaviors,” Connect. Tissue Res., 43(4), pp. 569–580. [CrossRef] [PubMed]
Broedersz, C. P., Sheinman, M., and MacKintosh, F. C., 2012, “Filament-Length-Controlled Elasticity in 3D Fiber Networks,” Phys. Rev. Lett., 108(6), p. 078102. [CrossRef] [PubMed]
Head, D. A., Levine, A. J., and MacKintosh, F. C., 2003, “Distinct Regimes of Elastic Response and Deformation Modes of Cross-Linked Cytoskeletal and Semiflexible Polymer Networks,” Phys. Rev. E, 68(6), p. 061907. [CrossRef]
Chandran, P. L. and Barocas, V. H., 2004, “Microstructural Mechanics of Collagen Gels in Confined Compression: Poroelasticity, Viscoelasticity, and Collapse,” ASME J. Biomech. Eng., 126(2), pp. 152–166. [CrossRef]
Knapp, D., Barocas, V. H., Moon, A., Yoo, K., Petzold, L., and Tranquillo, R. T., 1997, “Rheology of Reconstituted Type I Collagen Gel in Confined Compression,” J. Rheol., 41, p. 971. [CrossRef]
Hsu, S., Jamieson, A. M., and Blackwell, J., 1994, “Viscoelastic Studies of Extracellular Matrix Interactions in a Model Native Collagen Gel System,” Biorheology, 31(1), pp. 21–36. [PubMed]
Meghezi, S., Couet, F., Chevallier, P., and Mantovani, D., 2012, “Effects of a Pseudophysiological Environment on the Elastic and Viscoelastic Properties of Collagen Gels,” Int. J. Biomater., 2012, p. 319290. [PubMed]
Krishnan, L., Weiss, J. A., Wessman, M. D., and Hoying, J. B., 2004, “Design and Application of a Test System for Viscoelastic Characterization of Collagen Gels,” Tissue Eng., 10(1–2), pp. 241–252. [CrossRef] [PubMed]
Pryse, K. M., Nekouzadeh, A., Genin, G. M., Elson, E. L., and Zahalak, G. I., 2003, “Incremental Mechanics of Collagen Gels: New Experiments and a New Viscoelastic Model,” Ann. Biomed. Eng., 31(10), pp. 1287–1296. [CrossRef] [PubMed]
Wagenseil, J. E., Wakatsuki, T., Okamoto, R. J., Zahalak, G. I., and Elson, E. L., 2003, “One-Dimensional Viscoelastic Behavior of Fibroblast Populated Collagen Matrices,” ASME J. Biomech. Eng., 125(5), pp. 719–725. [CrossRef]


Grahic Jump Location
Fig. 2

Development of the coupled fiber-matrix model: (a) work-flow demonstrating construction protocol for the coupled model; (b) illustration of the conforming multidimensional mesh and the mesh classification to the geometric model; (c) schematic showing the interior of the multidimensional mesh showing fibers (black lines) and meshed matrix (yellow elements with blue borders)

Grahic Jump Location
Fig. 1

Computational representations of collagen fiber networks were constructed by randomly placing a set of seed points in a representative volume element, constructing a Voronoi tessellation about these points, defining fibers as the edges of the Voronoi elements, removing seed points, and placing pin-joint nodes at each edge-edge and edge-boundary intersection

Grahic Jump Location
Fig. 4

The magnitude of stress contributions differed depending on whether simulations were evaluated at constant total strain (solid lines) or constant total stress (dashed lines); (a) while matrix stress values increased in both cases, (b) fiber stresses show opposite trends for the two cases, where decreasing values for the constant-stress case demonstrates stress-shielding (by the matrix) at high shear moduli

Grahic Jump Location
Fig. 5

Average (a) matrix stress, (b) fiber stress, and (c) fraction of total stress at 10% strain and with νNFM = 0.1 show good agreement between the parallel and coupled models; stress values at a larger Poisson's ratio (i.e., νNFM = 0.45) at G = 110 Pa show a small shift of stress from the matrix to fibers (mean ± 95% CI; n = 5)

Grahic Jump Location
Fig. 6

Average (a) fiber stretch, (b) RVE Poisson's ratio, and (c) Ω11 (representing fiber orientation in the loading direction) at 10% strain and with νNFM = 0.1 show decreased values for the coupled model compared to the parallel model, but similar qualitative changes as a function of increasing G; for the case where νNFM = 0.1 and G = 110 Pa, the fiber stretch and Poisson's ratio increased and decreased, respectively, with no change in Ω11 (mean ± 95% CI; n = 5)

Grahic Jump Location
Fig. 7

Interior normal and shear stress fields at 10% strain on the mid-section slice for a representative network (G = 720 Pa; νNFM = 0.1) demonstrates a highly inhomogeneous distribution for all six independent stress components; slices were cut normal to the loading (1-) direction in the 2-3 plane (represented by the dashed lines in the RVE schematic) and black dots indicate locations where fibers intersect the cutting plane; examples of two regions of high stress concentrations are indicated by arrows in the σ23 plot

Grahic Jump Location
Fig. 3

(a) Fiber stress, (b) matrix stress, (c) total stress, and (d) RVE Poisson's ratio versus engineering strain for the coupled fiber-matrix model at varying values of the NFM shear modulus; the stress values increased with increasing G, particularly for the matrix, while Poisson's ratio decreased (data are representative results from a single network)

Grahic Jump Location
Fig. 8

(a) Maximum, and (b) minimum principal stress (normalized by the shear modulus of the matrix material) distributions at 10% strain over all mesh elements of the matrix material; the mean value increased and the standard deviation decreased as G increased, with curves moving more towards the affine model (the values from five networks are lumped and plotted)

Grahic Jump Location
Fig. 9

Distribution of fiber stretch values at 10% strain and with νNFM = 0.1; values concentrated near ∼1 when G was low (i.e., G = 10 Pa and 110 Pa), but fibers were stretched to a greater extent at higher G values, similar to what was shown in the affine model (the values from five networks are lumped and plotted)



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