In the field of computational biomechanics, investigators have primarily used commercial software that is neither geared toward biological applications nor sufficiently flexible to follow the latest developments in the field. This lack of a tailored software environment has hampered research progress, as well as dissemination of models and results. To address these issues, we developed the FEBio software suite (http://mrl.sci.utah.edu/software/febio), a nonlinear implicit finite element (FE) framework, designed specifically for analysis in computational solid biomechanics. This paper provides an overview of the theoretical basis of FEBio and its main features. FEBio offers modeling scenarios, constitutive models, and boundary conditions, which are relevant to numerous applications in biomechanics. The open-source FEBio software is written in C++, with particular attention to scalar and parallel performance on modern computer architectures. Software verification is a large part of the development and maintenance of FEBio, and to demonstrate the general approach, the description and results of several problems from the FEBio Verification Suite are presented and compared to analytical solutions or results from other established and verified FE codes. An additional simulation is described that illustrates the application of FEBio to a research problem in biomechanics. Together with the pre- and postprocessing software PREVIEW and POSTVIEW, FEBio provides a tailored solution for research and development in computational biomechanics.

Due to its avascular nature, articular cartilage exhibits a very limited capacity to regenerate and to repair. Although much of the tissue-engineered cartilage in existence has been successful in mimicking the morphological and biochemical appearance of hyaline cartilage, it is generally mechanically inferior to the natural tissue. In this study, we tested the hypothesis that the application of dynamic deformational loading at physiological strain levels enhances chondrocyte matrix elaboration in cell-seeded agarose scaffolds to produce a more functional engineered tissue construct than in free swelling controls. A custom-designed bioreactor was used to load cell-seeded agarose disks dynamically in unconfined compression with a peak-to-peak compressive strain amplitude of 10 percent, at a frequency of 1 Hz, 3 x (1 hour on, 1 hour off)/day, 5 days/week for 4 weeks. Results demonstrated that dynamically loaded disks yielded a sixfold increase in the equilibrium aggregate modulus over free swelling controls after 28 days of loading (100 +/- 16 kPa versus 15 +/- 8 kPa, p < 0.0001). This represented a 21-fold increase over the equilibrium modulus of day 0 (4.8 +/- 2.3 kPa). Sulfated glycosaminoglycan content and hydroxyproline content was also found to be greater in dynamically loaded disks compared to free swelling controls at day 21 (p < 0.0001 and p = 0.002, respectively).

Over the last two decades, considerable progress has been reported in the field of cartilage mechanics that impacts our understanding of the role of interstitial fluid pressurization on cartilage lubrication. Theoretical and experimental studies have demonstrated that the interstitial fluid of cartilage pressurizes considerably under loading, potentially supporting most of the applied load under various transient or steady-state conditions. The fraction of the total load supported by fluid pressurization has been called the fluid load support. Experimental studies have demonstrated that the friction coefficient of cartilage correlates negatively with this variable, achieving remarkably low values when the fluid load support is greatest. A theoretical framework that embodies this relationship has been validated against experiments, predicting and explaining various outcomes, and demonstrating that a low friction coefficient can be maintained for prolonged loading durations under normal physiological function. This paper reviews salient aspects of this topic, as well as its implications for improving our understanding of boundary lubrication by molecular species in synovial fluid and the cartilage superficial zone. Effects of cartilage degeneration on its frictional response are also reviewed.

The continuum mechanical treatment of biological growth and remodeling has attracted considerable attention over the past fifteen years. Many aspects of these problems are now wellunderstood, yet there remain areas in need of significant development from the standpoint of experiments, theory, and computation. In this perspective paper we review the state of the field and highlight open questions, challenges, and avenues for further development.

No abstract

A biphasic mixture model is developed which can account for the observed tension-compression nonlinearity of cartilage by employing the continuum-based Conewise Linear Elasticity (CLE) model of Curnier et al. (J Elasticity 37:1-38, 1995) to describe the solid phase of the mixture. In this first investigation, the orthotropic octantwise linear elasticity model was reduced to the more specialized case of cubic symmetry, to reduce the number of elastic constants from twelve to four. Confined and unconfined compression stress-relaxation, and torsional shear testing were performed on each of nine bovine humeral head articular cartilage cylindrical plugs from 6 month old calves. Using the CLE model with cubic symmetry, the aggregate modulus in compression and axial permeability were obtained from confined compression (H −A =0.64±0.22 MPa, k z = 3.62 ± .97 × 10 −16 m 4 /N.s, r 2 =0.95 ±0.03), the tensile modulus, compressive Poisson ratio and radial permeability were obtained from unconfined compression (E +Y = 12.75 ± 1.56 MPa, ν − =0.03±0.01, k r =6.06±2.10×10 −16 m 4 /N.s, r 2 =0.99±0.00), and the shear modulus was obtained from torsional shear (µ=0.17±0.06 MPa). The model was also employed to successfully predict the interstitial fluid pressure at the center of the cartilage plug in unconfined compression (r 2 =0.98±0.01). The results of this study demonstrate that the integration of the CLE model with the biphasic mixture theory can provide a model of cartilage which can successfully curvefit three distinct testing configurations while producing material parameters consistent with previous reports in the literature.

An optimal strategy for the functional tissue engineering of articular cartilage, particularly to accelerate construct development, may incorporate sequential application of different growth factors and applied deformational loading.

Taken together, these results point to the utility of dynamic deformational loading in the mechanical preconditioning of engineered articular cartilage constructs and the necessity for increasing feed media volume and serum supplementation with increasing cell seeding densities.