Coupled Macroscopic and Microscopic Scale Modeling of Fibrillar Tissues and Tissue Equivalents

[+] Author and Article Information
Balaji Agoram

Department of Chemical Engineering, University of Colorado, Boulder, CO 80309-0424

Victor H. Barocas

Department of Biomedical Engineering, University of Minnesota, Minneapolis, MN 55455

J Biomech Eng 123(4), 362-369 (Mar 31, 2001) (8 pages) doi:10.1115/1.1385843 History: Received March 07, 2000; Revised March 31, 2001
Copyright © 2001 by ASME
Your Session has timed out. Please sign back in to continue.


Silver, F. H., 1987, Biological Materials: Structure, Mechanical Properties, and Modeling of Soft Tissues, New York University Press, New York.
Tuderman,  L., and Bruckner,  P., 1998, “Genetic Diseases of the Extracellular Matrix: More Than Just Connective Tissue Disorders,” J. Mol. Med., 76, pp. 226–237.
Bell,  E., Ivarsson,  B., and Merrill,  C., 1979, “Production of a Tissue-Like Structure by Contraction of Collagen Lattices by Human Fibroblasts of Different Proliferative Potential in Vivo,” Proc. Natl. Acad. Sci. U.S.A., 76, pp. 1274–1278.
Lopez Valle,  C. A., Auger,  F. A., Rompre,  P., Bouvard,  V., and Germain,  L., 1992, “Peripheral Anchorage of Dermal Equivalents,” Br. J. Dermatol., 127, pp. 365–371.
Parenteau,  N., Sabolinski,  M., Nolte,  C., Oleson,  M., Kriwet,  K., and Bilbo,  P., 1996, “Biological and Physical Factors Influencing the Successful Engraftment of a Cultured Human Skin Substitute,” Biotechnol. Bioeng., 52, pp. 3–14.
Tranquillo,  R. T., Girton,  T. S., Bromberek,  B. A., Triebes,  T. G., and Mooradian,  D. L., 1996, “Magnetically-Oriented Tissue-Equivalent Tubes: Application to a Circumferentially-Oriented Media-Equivalent,” Biomaterials, 17, p. 349.
Fung, Y. C., 1993, Biomechanics: Mechanical Properties of Living Tissues, Springer-Verlag. New York.
Lanir,  Y., 1982, “Constitutive Equations for Fibrous Connective Tissues,” J. Biomech., 18, pp. 1–12.
Shoemaker,  P. A., Schneider,  D., Lee,  M. C., and Fung,  Y., 1986, “A Constitutive Model for Two-Dimensional Soft Tissues and Its Application to Experimental Data,” J. Biomech., 19, pp. 695–702.
Lai,  W. M., Hou,  J. S., and Mow,  V. C., 1991, “A Triphasic Theory for the Swelling and Deformation Behaviors of Articular Cartilage,” ASME J. Biomech. Eng., 113, pp. 245–258.
Simon,  B. R., Liable,  J. P., Pflaster,  D., Yuan,  Y., and Krag,  M. H., 1996, “A Poroelastic Finite Element Formulation Including Transport and Swelling in Soft Tissue Structures,” ASME J. Biomech. Eng., 118, pp. 1–9.
Farquhar,  T., Dawson,  P. R., and Torzilli,  P. A., 1990, “A Microstructural Model for the Anisotropic Drained Stiffness of Articular Cartilage,” ASME J. Biomech. Eng., 112, pp. 414–425.
Soulhat,  J., Buschmann,  M. D., and Shirazi-Adl,  A., 1999, “A Fibril-Network-Reinforced Biphasic Model of Cartilage in Unconfined Compression,” ASME J. Biomech. Eng., 121, pp. 340–347.
Schwartz,  M. H., Leo,  P. H., and Lewis,  J. L., 1994, “A Microstructural Model for the Elastic Response of Articular Cartilage,” J. Biomech., 27, pp. 865–873.
Hollister,  S. J., Fyrhie,  D. P., Jepsen,  K. J., and Goldstein,  S. A., 1991, “Application of Homogenization Theory to the Study of Trabecular Bone Mechanics,” J. Biomech., 24, pp. 825–839.
Hollister,  S. J., Brennan,  J. M., and Kikuchi,  N., 1994, “A Homogenization Sampling Procedure for Calculating Trabecular Bone Effective Stiffness and Tissue Level Stress,” J. Biomech., 27, pp. 433–444.
Barocas,  V. H., and Tranquillo,  R. T., 1997, “An Anisotropic Biphasic Theory of Tissue-Equivalent Mechanics: The Interplay Among Cell Traction, Fibrillar Network Deformation, Fibril Alignment, and Cell Contact Guidance,” ASME J. Biomech. Eng., 119, pp. 137–145.
Sherratt,  J. A., and Lewis,  J., 1993, “Stress-Induced Alignment of Actin Filaments and the Mechanics of Cytogel,” Bull. Math. Biol., 55, pp. 637–654.
Dembo,  M., and Harlow,  F., 1986, “Cell Motion, Contractile Networks, and the Physics of Interpenetrating Reactive Flow,” Biophys. J., 50, pp. 109–121.
Stamenovic,  D., Fredberg,  J. J., Wang,  N., Butler,  J. P., and Ingber,  D. E., 1996, “A Microstructural Approach to Cytoskeletal Mechanics Based on Tensegrity,” J. Theor. Biol., 181, pp. 125–136.
Strang, W. G., and Fix, J. G., 1973, An Analysis of the Finite Element Method, Prentice-Hall, Englewood Cliffs, NJ.
Gill,  P. E., Jay,  L. O., Leonard,  M. W., Petzold,  L. R., and Sharma,  V., 2000, “An SQP Method for the Optimal Control of Large-Scale Dynamical Systems,” J. Comput. Appl. Math., 120, pp. 197–213.
Gill, P., Murray, W., and Saunders, M. A., 1998, “SNOPT: a FORTRAN Package for Large-Scale Nonlinear Programming.”
Draper, N. R., and Smith, H., 1966, Applied Regression Analysis, Wiley, New York.
Guido,  S., and Tranquillo,  R. T., 1993, “A Methodology for the Systematic and Quantitative Study of Cell Contact Guidance in Oriented Gels: Correlation of Fibroblast Orientation and Gel Birefringence,” J. Cell. Sci., 105, pp. 317–331.
Klebe,  R. J., Caldwell,  H., and Milam,  S., 1990, “Cells Transmit Spatial Information by Orienting Collagen-Fibers,” Matrix, 9, pp. 451–458.
Bergren, T. E., 1993, “Controlling the Material Properties of in Vitro Collagen Through Directed Orientation and Cross-Linking of Fibrils,” Ph.D. Thesis, Department of Aerospace Engineering, University of Colorado, Boulder, CO.
Hibbeler, R. C., 1994, Mechanics of Materials, MacMillan College Publishing Company.
Mow,  V. C., Kuei,  S. C., Lai,  W. M., and Armstrong,  C. G., 1980, “Biphasic Creep and Stress Relaxation of Articular Cartilage in Compression: Theory and Experiments,” ASME J. Biomech. Eng., 102, pp. 73–84.
Knapp,  D. M., Barocas,  V. H., Moon,  A. G., Yoo,  K., Petzold,  L. R., and Tranquillo,  R. T., 1997, “Rheology of Reconstituted Type I Collagen Gel in Confined Compression,” J. Rheol., 41, pp. 971–993.
Grinnell,  F., and Lamke,  C. R., 1984, “Reorganization of Hydrated Collagen Lattices by Human Skin Fibroblasts,” J. Cell. Sci., 66, pp. 51–63.


Grahic Jump Location
Flowchart of the solution strategy. Our goal is to incorporate microscopic information into the solution of a macroscopic problem. Once the macroscopic problem has been posed (Step 1), the macroscopic nodal positions are projected onto the microscopic domain within each element to generate a set of microscopic problems (Step 2). These problems are solved (independently) to determine how the microstructure must move to accommodate the macroscopic deformation and how much elastic force results from the motion (Step 3). The microscopic forces are then used to determine the macroscopic stress at every node (Step 4), allowing the evaluation of the residual equation, Eq. (4), in Step 5.
Grahic Jump Location
The macroscopic problem: Although our method could be applied generally to any stress-deformation problem, we focus here on simple extension of a rectangular (two-dimensional) tissue equivalent (a). The sample is assumed to be clamped on the left and right edges (solid lines) but free on the top and bottom edge (dashed lines). The sample is deformed (b) by moving the clamped edges, which allow no y displacement and require a specified x displacement. The free surface then deforms in response to the resulting internal stresses (deformation greatly exaggerated).
Grahic Jump Location
Representation of FEM model: A micro–macro problem with four finite elements is shown. Macroscopic (finite element) nodes are shown by solid circles. The microscopic scale networks are all the same, but this is not a requirement.
Grahic Jump Location
Schematic of fictional segments: Solid lines indicate real fiber segments. Broken lines indicate fictional segments that provide cost to bending of solid segments about a node. S represents a fiber segment, and F1 and F2 represent two fictional segments. Displacement of segment S (shown by S) causes tension in fictional segment F1 and compression in fictional segment F2.
Grahic Jump Location
Accuracy of micro–macro model: To estimate discretization error in the FEM model stresses (in arbitrary units) for the FEM model and a direct microscopic problem for same TDOF are plotted for three different TDOF. Bars to left correspond to microscopic model. An average difference of 11 percent was found between the two methods.
Grahic Jump Location
Error in micro–macro stress conversion: The percentage variation of s about the estimated value at the 95 percent confidence level is plotted against μDOF The error decreases with increasing μDOF and becomes roughly constant beyond μDOF=120.
Grahic Jump Location
Scaling of models: CPU times are plotted for (i) FEM model versus TDOF for increasing finite elements; (ii) FEM model versus TDOF for increasing μDOF; (iii) direct solution versus TDOF. Best scaling is obtained when we increase the number of finite elements each with a constant μDOF.
Grahic Jump Location
Model predictions of typical anisotropic TE behavior (Ψ=0.62): Higher forces are seen parallel to the direction of fiber orientation (∥) than perpendicular to the direction of orientation (⊥). Intermediate forces are seen for an isotropic TE (Iso).
Grahic Jump Location
Scanning electron micrographs of oriented and unoriented collagen gels: Gel (a) was 1×3 cm, and constrained compaction led to preferential alignment of the collagen fibrils in the horizontal direction. Gel (b) was 1×1 cm, leading to isotropic compaction and thus no alignment of the collagen fibrils. Scale bar in (a) is 1 μm.
Grahic Jump Location
Experimental (a) (Ψ=0.61) and simulated (b) (Ψ=0.62) behavior of isotropic and anisotropic gels. Results for an anisotropic gel tested parallel to (∥) and perpendicular (⊥) to the direction of orientation are shown, as are results for extension of an isotropic gel (Iso). At 8 percent strain, the simulation is two times higher for F, close to the same for FIso, and five times lower for F.




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