- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Computational and Mathematical Methods in Medicine
Volume 2013 (2013), Article ID 326150, 23 pages
A Review of the Combination of Experimental Measurements and Fibril-Reinforced Modeling for Investigation of Articular Cartilage and Chondrocyte Response to Loading
1Department of Clinical Neurophysiology, Kuopio University Hospital, FI-70211 Kuopio, Finland
2Department of Biomedical Engineering, Eindhoven University of Technology, 5600 MB, Eindhoven, The Netherlands
3Department of Applied Physics, University of Eastern Finland, FI-70211 Kuopio, Finland
4Division of Solid Mechanics, Department of Orthopaedics, Lund University, 221 00 Lund, Sweden
5Human Performance Laboratory, Faculty of Kinesiology, University of Calgary, Alberta, BC, Calgary, T2N 1N4, Canada
Received 20 September 2012; Revised 11 January 2013; Accepted 23 February 2013
Academic Editor: Leping Li
Copyright © 2013 Petro Julkunen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The function of articular cartilage depends on its structure and composition, sensitively impaired in disease (e.g. osteoarthritis, OA). Responses of chondrocytes to tissue loading are modulated by the structure. Altered cell responses as an effect of OA may regulate cartilage mechanotransduction and cell biosynthesis. To be able to evaluate cell responses and factors affecting the onset and progression of OA, local tissue and cell stresses and strains in cartilage need to be characterized. This is extremely challenging with the presently available experimental techniques and therefore computational modeling is required. Modern models of articular cartilage are inhomogeneous and anisotropic, and they include many aspects of the real tissue structure and composition. In this paper, we provide an overview of the computational applications that have been developed for modeling the mechanics of articular cartilage at the tissue and cellular level. We concentrate on the use of fibril-reinforced models of cartilage. Furthermore, we introduce practical considerations for modeling applications, including also experimental tests that can be combined with the modeling approach. At the end, we discuss the prospects for patient-specific models when aiming to use finite element modeling analysis and evaluation of articular cartilage function, cellular responses, failure points, OA progression, and rehabilitation.
As a fluid-saturated composite tissue, articular cartilage provides smooth sliding at joint surfaces. Articular cartilage also absorbs forces and reduces stresses experienced by bones. The extracellular matrix (ECM) of cartilage can be divided in two phases: solid and fluid. The solid phase is primarily composed of cartilage matrix proteins: proteoglycans (PGs) and collagen (mainly type II). Cartilage also contains chondrocytes, that is, articular cartilage cells, which are surrounded by a thin layer called the pericellular matrix (PCM). The characteristic fibrous structure of the ECM is a result of depth-wise remodeling of the collagen fibril network during maturation [1–4] and represents a three-layer laminar architecture in adults (Figure 1) . The collagen network is known to be arranged into parallel planes, revealing split-line patterns in cartilage surfaces . The well-organized collagen fibril bundles link to each other with smaller cross-link chains.
The mechanical response of articular cartilage is determined by the tissue composition and structure [6–10]. Due to the fixed negative charges of the PGs, the total ion concentration inside the tissue is higher than in the surrounding synovial fluid (cation concentration is higher and the anion concentration is lower), which leads to a difference in osmotic pressure that results in tissue swelling . The swelling pressure modulates the shape and volume of the tissue by controlling the fluid content of the tissue . The collagen network of the ECM stabilizes the matrix by balancing the expansion caused by the swelling pressure [13, 14].
The collagen network is believed to mainly affect the response of cartilage under dynamic or instantaneous loading [6, 7, 19–22], while its direct impact on the equilibrium response is minor [6, 7, 19, 21]. Oloyede and Broom utilized a consolidation theory and experiments to understand the role of interstitial fluid flow for the mechanical response and deformation behavior of loaded cartilage tissue [23–25]. During loading the interstitial fluid flows within the cartilage solid matrix. The flow velocity is determined by the induced fluid pressure and matrix permeability. Since the fluid flow is restricted by the cartilage matrix during mechanical loading, the internal pressure of the tissue increases. Therefore, the interstitial fluid also directly affects the dynamic and instantaneous mechanical responses of cartilage.
Chondrocytes produce the ECM matrix components. The activity of chondrocytes is affected by genetic and environmental factors, electrokinetic forces, mechanical forces, and fluid pressurization. Chondrocytes change shape and volume in response to mechanical loading of cartilage [26–28]. The ECM and PCM possess organized collagenous structures, and the collagen fibril orientation is affected by local tissue strains , which further modulate cell shape and volume . Consequently, these factors modulate the PCM and ECM composition [29–31]. Through a sequence of mechanobiological events, the mechanical signals from the ECM are converted into intracellular responses, which serve to maintain or alter cartilage ECM. The PCM around the cells acts as a transducer for the biomechanical and biophysical signals perceived by the cells and thus alters the cell activity .
Osteoarthritis (OA) is a joint disease, characterized by the progressive degeneration of articular cartilage. PG depletion occurs in the superficial cartilage during the early stage of OA [32–35]. Simultaneous alterations in the collagen network also take place [32, 35–37]. These alterations include a dramatic change in orientation of the superficial collagen fibrils [32, 34, 38, 39] and a reduction in the collagen content . Thereby, these changes increase cartilage hydration via intake of water. Subsequently, tissue permeability increases. Together, these changes impair the mechanical function of cartilage by decreasing its mechanical stiffness , which may further accelerate the progression of OA. As a result of the OA changes in tissue composition and structure, the extracellular osmolarity increases. This results in increased cell volume, which may in turn affect tissue biosynthesis and cartilage mechanotransduction.
Changes in the cartilage structure and composition can be revealed by using in vivo imaging techniques, such as magnetic resonance imaging (MRI) or contrast-enhanced computed tomography (CECT) or in vitro techniques such as quantitative microscopy [2, 34, 40–46]. However, these imaging techniques alone are unable to quantify mechanical characteristics of cartilage that are critical to the demanding mechanical function in the joint. Instead, computational modeling in conjunction with imaging techniques is needed to predict the mechanical performance of cartilage and other tissues during joint loading [17, 32, 47–53].
Chondrocyte responses to osmotic or mechanical loading have been experimentally characterized by using confocal or dual-photon microscopy [54–56]. In earlier studies, both mechanical responses, for example, volume and morphology changes, and biological responses, such as Ca2+ signaling, have been characterized [54–57]. The measurements have traditionally been conducted on isolated cells or chondrocytes have been measured through the cut surface of tissue explants. Recently, cell responses to loading have also been conducted in situ using fully intact tissue samples, providing more physiological environment for cells [55, 56]. Like at the tissue level, the stresses, strains, and fluid flow in and around the cells cannot be measured experimentally. It is also experimentally challenging, if not impossible, to specifically investigate the roles of different cartilage constituents (e.g., collagen, PGs, fluid, and ions) on the cell responses. Therefore, computational models have been employed [58, 59].
The modern fibril reinforced computational models of articular cartilage may be implemented to include inhomogeneous tissue composition and complex structure [8, 32, 53, 59–68]. These models can simulate nonlinear behavior of the tissues, caused by tissue anisotropy and inhomogeneity. Cartilage models have been employed to evaluate the mechanical behavior of cartilage at the cell, tissue, and joint level, to evaluate static and dynamic tissue behavior, to explore the effects of mechanical and biochemical loading, and to predict tissue remodeling, growth, and adaptation over time. As a next step, the 3D models may be developed and validated to include patient-specific tissue characteristics. Ultimately, they may then become clinical tools for predicting the progression of OA and thus for identifying or optimizing patient-specific treatment strategies. To achieve this goal, a smooth transition is needed from joint-level models to cell-level imaging and modeling. This could help to characterize the effects of joint loads on tissues, and cells, to investigate possible failure sites in tissues, and to predict cell-based tissue adaptation to loading in normal and diseased joints. In this paper, we concentrate on cell- and tissue-level models and point out their future possibilities. Theories of present fibril reinforced models are reviewed, and major challenges in validation and application of these models are critically discussed.
2. Review of Fibril Reinforced Computational Models of Articular Cartilage
This paper focuses on the application of fibril reinforced mixture models. The application of fibril reinforcement in the models is justified by the heavily influential collagen fibril network for the cartilage mechanical behavior [69, 70]. Oloyede and Broom considered that there are six aspects influencing the deformation behavior of cartilage which a theoretical or physical model should incorporate : (1) the rate of fluid outflow from the tissue matrix, (2) the diffusion process or movement of interstitial fluid between adjacent regions of the matrix, (3) deformation dependent permeability, (4) loading-rate-dependent viscous drag, (5) physico-chemistry and osmotic swelling, and (6) the nonlinearity of the stress-strain characteristic of the solid skeleton or structural framework.
It is known that articular cartilage exhibits viscoelastic behavior [72–77]. The viscoelastic behavior is often divided to flow-dependent and flow-independent behavior [74, 78]. Prior to the development of fibril reinforced models, the biphasic mixture theory (fluid and solid phase) was introduced for cartilage by Mow et al. . Use of the biphasic (and poroelastic) material models can at least partly account for the flow-dependent viscoelastic behavior, which arises from fluid flow within a porous material with low permeability [75, 76, 79–84]. The biphasic mixture theory has been demonstrated to agree well with direct experimental measures of fluid pressurization [23, 24, 83, 85]. While the theory was considered successful in many aspects, it failed to predict certain mechanical behavior observed under mechanical tests [86, 87]. The models were improved by transversely isotropic elastic  and conewise linear elastic  implementation of the solid matrix. Use of biphasic transversely isotropic models produced better curve fits as compared to biphasic linear elastic models [88, 89]. The transversely isotropic models were able to capture one aspect of the compression-tension nonlinearity of a cartilage having different mechanical properties in the direction of compression and perpendicular to that direction. However, those models could not account for the experimentally observed compression-tension nonlinearity, caused by the collagen fibrils, measured between compressive and tensile tests [21, 90, 91]. The conewise linear elastic model overcame this issue via a bimodular approach setting the material properties in the model different during compressive and tensile strain . Soulhat et al.  implemented a composite fibril reinforced biphasic model in which the collagen fibrils were represented as springs over the solid matrix, resisting deformation only in tension. These are evolutionary steps taken towards present models used to simulate the mechanical behavior of articular cartilage. Further development of modern fibril reinforced models from that introduced by Soulhat et al. is reviewed later in this section.
Since it has been difficult to separate the flow-independent from the flow-dependent viscoelasticity, it has been questioned whether cartilage possesses flow-independent viscoelastic properties at all [80, 81, 93]. Some experimental evidence supports the existence of flow-independent viscoelastic properties in cartilage [94–99], which can partly be attributed to the collagen fibrils and partly to a lesser extent to the nonfibrillar matrix or their interaction [97, 100–105]. For modeling cartilage as a biphasic viscoelastic material, Mak  combined the biphasic theory with a quasilinear viscoelastic theory for the solid matrix by Fung . This allowed for the separation of flow-dependent and flow-independent viscoelasticity. The more recent fibril reinforced models developed for articular cartilage have often implemented the fluid-independent viscoelastic behavior of the nonfibrillar solid, viscoelastic fibrils, or both [8, 12, 18, 68, 81, 108].
Recent mechanical simulation studies including articular cartilage have often implemented a fibril reinforced material model. However, for more complex geometries, the implementation of fibril reinforced articular cartilage is often omitted. Mostly, these approaches are used when implementation of fibril reinforcement is considered unnecessary (justified simplification) for the aim of the study. The inclusion of fibril reinforcement may sometimes be unnecessary due to the role of articular cartilage for the simulated result. Excessive computational cost may provide reasoning for simplification of the cartilage material model. For instance, it is still customary to model articular cartilage as isotropic elastic material during complete joint simulations [109–113]. However, recent models of a joint have included fibril reinforcement when simulating joint mechanics with focus on cartilage [53, 114–117]. Adouni et al. found that use of a fibril reinforced model in a total knee joint model exhibited larger contact pressures at knee joint contact surfaces when comparing to a simplified isotropic material model . Often, the use of simplified material implementation is justified due to the short-term type of simulation during which viscoelastic behavior does not play a great role and articular cartilage appears as an incompressible and elastic solid . Simplifying the material behavior of articular cartilage as an incompressible, single-phase material is appropriate when rapid loading with short duration is investigated. The equivalence of the short-time biphasic and incompressible elastic responses has been demonstrated . In the following, we present fibril reinforced mixture models more in detail.
2.1. Tissue-Level Fibril Reinforced Models
The fibril reinforced modeling of articular cartilage at the tissue level started in the late nineties. Then for the first time, a model included specifically both the collagen network and the poroelastic matrix (PGs and fluid) . Since then, development of computational models of cartilage has advanced quickly to include complex material characteristics models with more realistic behavior (Table 1). The main features of the fibril reinforced models of cartilage tissue are presented below.
The first fibril reinforced models were composed of the fibrillar part, representing the collagen fibrils, and the nonfibrillar part, mimicking the porous material (mainly PGs) filled with fluid [7, 119]. Total stress () in this material can be expressed as: where is the nonfibrillar matrix stress, is the fibril network stress, is the fluid pressure, and is the unit tensor.
2.1.1. The Nonfibrillar Matrix
The nonfibrillar part was modeled as linear elastic material in the original models obeying Hooke’s law, while at larger strains a hyperelastic, Neo-Hookean material description is preferable [17, 18]. The constitutive equation for the Neo-Hookean material is where and are temperature-dependent material parameters, is Cauchy-Green deformation tensor, and is the elastic volume ratio. The material parameters can be expressed as where is the initial shear modulus and is Poisson’s ratio of the nonfibrillar matrix. The shear modulus can be defined via Young’s modulus of the nonfibrillar matrix () and :
Fluid flow in the nonfibrillar matrix has been generally modeled with Darcy’s law: where is the flow rate, is the permeability, and is the pressure gradient across the region. The permeability can be defined to be dependent on the porosity and void ratio, that is, ratio of fluid to solid fraction, for example, according to the following equation given by van der Voet : where is the initial permeability, and are initial and current void ratios, and is a positive constant [18, 120]. To be consistent with experimental findings, the fluid fraction that is used to compute the void ratios has also been implemented in the models in a depth-dependent manner over the cartilage thickness [18, 61, 64, 121–123]. In the end, the required material parameters for the nonfibrillar matrix of cartilage are Young’s modulus (), Poisson’s ratio (), and permeability (). To include the effect of fibril reinforcement in the permeability, Federico and Herzog implemented the anisotropy of permeability by considering the orientation of the collagen fibrils in any given location in the tissue in the undeformed state . Few recent studies have also considered intrinsic anisotropy of permeability due to the collagen fibrils [67, 125, 126]. Imaging studies measuring diffusion of water within the cartilage matrix have observed anisotropy of diffusion, mainly affected by the general collagen fibril orientation [127–129]. In addition, the variation in excess pore pressure in the different directions parallel and perpendicular to the surface of cartilage provides experimental support for the implementation . If we were to include the anisotropy of permeability in (6), the results would be where and are the initial and current permeability in the direction of collagen fibrils, respectively, and are the initial and current permeability across the collagen fibrils, and and are material parameters. Ratio of and in any given point in the model would then depend on the orientation of the collagen fibrils.
2.1.2. Collagen Fibril Network
We divide the fibril reinforcement models into three types: (1) transversely isotropic [69, 88, 134] or conewise linear elastic , (2) fibril reinforcement with membrane elements or springs [7, 119, 131], and (3) vector-based fibril reinforcement [8, 12, 18, 68, 108, 135]. These types of models have been applied to simulate articular cartilage mechanics during the past decade. We will mainly focus our paper on the vector-based fibril reinforcement.
In the early fibril reinforced models, the fibrillar matrix mimicked the mechanical effect of the collagen fibrils of cartilage . The properties of the fibril network were controlled by Young’s modulus of the fibril network (). The elastic properties of the fibril network were further characterized with a nonlinear relation as such relation has been reported in the literature : where is the initial fibril network modulus, is the strain-dependent fibril network modulus, and is the fibril strain [7, 119]. Fibrils were presented with spring elements [7, 119]. Another approach was to use nonlinear membrane elements . As noticed from (8) and (9), the characteristic feature of the fibrils is that they resist only tensile forces. This is a major difference between transversely isotropic and spring-based models, the former having the same stiffness in compression and tension. Collagen fibrils have been demonstrated to possess flow-independent viscoelasticity [100–102]. Therefore, the collagen fibril stresses () were later modeled as viscoelastic: where is the viscoelastic damping coefficient, and and are the stress and strain rates, respectively [17, 18]. Unlike in type 1 and 2 models, the fibrils in the vector-based fibril reinforced models (type 3) can be constructed into any structural arrangement with fibril vectors (Figure 2) [18, 47, 53, 63, 68, 108]. This allows for the implementation of realistic fibril orientation to be independent of the mesh in use. The fibril vectors () can be used for determining logarithmic fibril strain used in (10) : where is the deformation tensor.
The fibrillar matrix in the model can be divided into primary and secondary fibrils . The primary fibrils may be assumed to represent the arcade-like collagen architecture, as detected with the polarized light microscopy [69, 137]. This network creates a depth-dependent tensile modulus for the tissue. As is commonly thought the fibrils are oriented parallel to the cartilage surface in the superficial zone, curve in the middle zone, and turn to perpendicular orientation with the cartilage surface in the deep zone . The secondary fibrils mimic the less organized part of the collagen network which is observed in scanning electron microscopy . When looking closely, the collagen fibrillar network appears random but exhibits highly directional configuration overall [37, 139]. The stresses for the primary and secondary fibrils ( and , resp.) can be formulated as follows: where represents the fibril volume fraction and is the density ratio of primary and secondary fibrils. The stress of the fibril network () is finally determined as the sum of the stresses in each individual fibril (): where is the total number of fibrils used to implement the fibrillar structure in an integration point.
2.1.3. Extensions to Fibrillar and Nonfibrillar Matrices
Tissue swelling is an essential part of cartilage mechanical behavior. Swelling properties of the fixed charge density in PGs have been presented in previous swelling cartilage models [6, 8, 12, 16, 59, 60, 64, 132, 140, 141]. Swelling of cartilage is influenced by osmotic swelling due to an excess amount of ion particles inside the tissue and chemical expansion due to repulsion of closely spaced negatively charged PGs [141, 142]. Lai et al.  introduced a triphasic mixture model, which includes three phases: an incompressible solid, an incompressible fluid, and a monovalent ionic phase. In the triphasic model, swelling behavior is explained by three mechanisms: Donnan osmosis, excluded volume effect, and chemical expansion . Simon et al.  introduced an alternative model by including the swelling effects in fluid equilibrium. The fundamental difference between these two approaches is that the model by Lai et al. causes fluid flow within the tissue, while the model by Simon et al. considers swelling to be an integral part of the fluid in equilibrium . Olsen et al. formulated a swelling model based on Biot’s principle of effective stress by including the swelling directly in the effective stress [25, 141, 146]. Huyghe and Janssen extended the triphasic theory to a quadriphasic mixture theory which includes electric flux and potential gradients at finite deformations [147, 148]. A biphasic swelling model introduced by Wilson et al.  simplifies the triphasic and quadriphasic models by neglecting the electrolyte flux [143, 147] and assumes that the ion concentration remains always in equilibrium.
After including the osmotic swelling in (1), the total stress of the material becomes where is the osmotic pressure gradient and is the chemical potential of fluid [12, 59, 147, 149]. The osmotic pressure gradient is caused by the difference in ion concentration of the cartilage and that of the surrounding fluid [12, 147, 149]. It is also referred to as the Donnan swelling pressure gradient. However, the mere osmotic swelling theory seems to underestimate the amount of swelling . Two possible explanations for that are (1) the distinction of intra- and extrafibrillar water  and (2) chemical expansion . The inclusion of chemical expansion further extends (15) to where is the chemical expansion stress. This model has been applied earlier specifically for cartilage since its swelling properties due to the fixed charge density have a significant role for the deformation behavior of the tissue, especially under static loading. For the implementation of swelling properties, values of the fixed charge density can be obtained from experimental measurements [152, 153]. The fixed charge density () affects both osmotic swelling and chemical expansion [142, 143]. In the biphasic swelling model, chemical expansion is only a function of . However, the inclusion of the chemical expansion term as proposed by Lai et al.  may not be appropriate, as it conflicts with the second law of thermodynamics. Huyghe et al. demonstrated that the chemical expansion stress part could produce free energy during closed loading cycles if the theory were true . Nevertheless, osmotic swelling alone may not fully simulate the swelling behavior in cartilage . However, when splitting the fluid phase into intra- and extrafibrillar portions [8, 150], implementation of chemical potential may not be needed.
As discussed earlier, the swelling of cartilage is resisted by the collagen network, inducing prestresses in the collagen fibrils. There are also other representations to the mechanical behavior of collagen fibrils than those presented above. Specifically, the nonlinear stress-strain tensile behavior of the collagen fibrils has been presented as follows: where is the first Piola-Kirchhoff fibril stress, is the total fibril strain, is the strain of the spring , and , , , and are constants [8, 64]. Again, the stress is computed only when the fibrils are in tension (i.e., and ; otherwise ).
2.1.4. Optional Solutions for Fibril Reinforced Material Simulation
The above-listed models present only the development of one branch of fibril reinforced models, applied in several studies. However, several other models used for understanding cartilage mechanics and etiology of OA have been proposed and applied as well. The principles of a few of these optional solutions are outlined below.
Li et al. implemented spring elements to simulate fibril reinforcement . This implementation combined two types of elements (springs and poroelastic matrix) overlaid on top of each other. Such implementation allowed for easy description of compression-tension nonlinearity in the cartilage, with collagen fibrils affecting cartilage total stress only in tension. In a quite similar manner, Shirazi and Shirazi-Adl simulated fibrillar matrix using directional membrane elements overlaid over poroelastic matrix . The models provide similar results .
Olsen and Oloyede implemented an overlay model , which was designed to study articular cartilage fibril reinforcement in a swollen construct. They later formulated and applied the model for simulating transient responses . In the model, cartilage solid was considered as a multilayer matrix, the layers of which all undergo the same deformation and contribute to the total stress. Using the overlay method, several different components could be applied in a single model, although the interactions between the overlays could not be implemented.
Ateshian et al.  implemented the solid matrix of cartilage with a continuous fibril angle distribution including osmotic swelling behavior in the nonfibrillar matrix. They only simulated the equilibrium responses subsiding interstitial fluid flow and flow-independent viscoelasticity. The application of continuous fibril angle distribution was able to correctly simulate the compression-tension nonlinearities in tissue stress and deformation. Ateshian et al. found that most of the conducted predictions could not be obtained from models using only three orthogonal fibril bundles. However, they considered it possible that some of their simulated features may also be captured by discrete (vector-based) fibril models with more than three orthogonal fibril bundles.
Holzapfel and Gasser described a viscoelastic fibril reinforced structural model  applicable also for simulating mechanics of articular cartilage with three-dimensional stress and deformation. This model falls into the vector-based fibril reinforced models with two fiber directions embedded in an isotropic ground matrix describing the composite fibril reinforcement. Gasser et al. further developed the model to incorporate a structural parameter characterizing the dispersion of collagen fibrils , which was later applied for articular cartilage .
Pierce et al.  proposed a structurally motivated 3D fibril reinforced mixture model. Their model was able to simultaneously address the dependence of fibril reinforcement and permeability from the fibrillar structure but omitted flow-independent viscoelasticity. The model included a dispersion parameter representing the degree of anisotropy within the fibrillar matrix, which could potentially be adjusted to account for collagen fibrillation associated with tissue degeneration. The model was able to reproduce mechanical behavior of articular cartilage under a variety of loading scenarios and demonstrated good agreement with experimental results.
Quinn and Morel implemented an intuitive model for simulation of collagen network mechanics and interactions between the fibrillar and nonfibrillar matrices . In their model, the fibril density and orientation can be utilized. Direction of the collagen fibrils is determined much like in the vector-based models presented above. This model separates the hydrated and dry portions of collagen and PGs, while also including “free” water and considering the structure of the fibrillar matrix. Although the solution presented by Quinn and Morel was restricted to equilibrium state, it could be extended to more general applications.
2.1.5. Controversies and Limitations Related to Tissue-Level Fibril Reinforced Models
Not all aspects of the biological cartilage tissue have been considered when attempting to simulate its mechanical behavior, and more specifically that of its fibrillar matrix. Limitations and challenges for implementation are caused by the different scales defining the mechanical appearance in the structure of the collagen fibrils and the coupling between chemistry, biology, and mechanical behavior. Aspects not often considered at tissue level include physical phenomena such as buckling of the fibrils [156–158], fibril-fibril interactions , fibril-PG interactions , or length and width of the collagen fibrils. The present fibril reinforced models also neglect some basic physiological phenomena such as fluid-pressurization-induced contact lubrication effects .
Buckling behavior of the collagen fibrils is rarely considered in fibril reinforced models of articular cartilage. However, it is likely that fibril buckling occurs in a confined space . Schinagl et al. suggested that a microstructurally motivated model considering buckling of the collagen fibrils at a critical load might account for the experimentally observed softening of the material with increasing load, following an initially stiff response . Bursac et al. modeled cartilage through a network of cables under pretension, representing collagen submitted to osmotic pressure. Under small deformations, the fibrils are in tension, which is reflected as high network stiffness. Instead, at larger deformations, the forces within the fibrils (represented as cables) become compressive, eventually leading to buckling with the result of reduced stiffness . Hence, buckling of the fibrils in cartilage will likely contribute to the compressive stiffness of the tissue at large deformations.
Fibril-fibril interactions have not been considered in the tissue-level fibril reinforced models. Nevertheless, the characteristics of the fibrillar matrix are likely highly dependent on the micromechanical interactions of the collagen fibrils. For instance, cross-linking has been demonstrated to affect the properties of the collagen fibril network [136, 162]. Buehler demonstrated that high cross-link density is associated with high yield strength and stress and brittle failure, while low cross-link density is associated with more ductile behavior with strain softening following peak stress . Chen and Broom suggested that breakdown in the interconnectivity of neighboring fibrils could lead to strongly aligned structure [163, 164]. Hence, any degenerative change in the fibrillar matrix could lead to rearrangement of fibrils along with varying degrees of fibrillar alignment. The loss of interconnectivity in the fibrillar network leaving the PGs intact could lead to reduction in PG entrapment accompanied by increased swelling tendency and decreased matrix stiffness . Broom and Silyn-Roberts reported that matrix cohesion is a consequence of fibril-fibril linkage independent of PGs, and biomechanical models involving shear stress transfer between the fibrillar and nonfibrillar matrix may therefore be inappropriate .
Fibril network and nonfibrillar solid matrices have obvious mechanical interaction considering that the osmotic swelling due to fixed charges within the nonfibrillar matrix is balanced by the tensile properties of the fibrillar network . Quinn and Morel highlighted the interactions between the collagen network and the PG gel . They explained the collagen fibril prestress in free-swelling cartilage emerging naturally and without introduction of additional artificial model parameters. Prestressing the collagen fibrils can be simulated during free swelling.
According to the implementation of the nonfibrillar matrix as presented above (Section 2.1.1), the nonfibrillar matrix can resist load on its own without the fibrillar matrix confining it. However, this may not be correct [14, 139, 167]. Through the interaction between the PG gel and the fibrillar matrix, the nonfibrillar matrix is known to contribute to the compressive tissue stiffness [7, 70, 90, 168] and to greatly resist fluid flow [169, 170]. In the aforementioned mixture-type models, the fibrillar matrix could be theoretically removed from the mixture leaving only the fluid-saturated nonfibrillar matrix, which would still be able to resist load as a porous elastic or hyperelastic matrix. This controversy in the interaction of the fibrillar and nonfibrillar matrix is not often considered in the implementation of the fibril reinforced models. Instead, the nonfibrillar matrix is actually considered to include the effects of fibrillar matrix confinement of the PG gel. This is due to implementation of infinitesimally thin fibrils (and fibril bundles), using for instance, springs. Hence, the fibril reinforced models do not fully separate the collagen and PG contributions to cartilage mechanics, as the drained nonfibrillar matrix possesses a nonphysical Poisson’s ratio giving the appearance that collagen fibrils support the nonfibrillar matrix even at zero tensile stress [14, 155]. Therefore the fibril reinforced models are still limited in their potential for estimation of cartilage mechanical properties from independently determined microstructural data and molecular physics .
Microstructural data such as fibril length and diameter are rarely considered in the implementation of tissue-level fibril reinforced models. However, it must be acknowledged that fibril length and diameter have impact for their mechanical performance affecting buckling and reorganization of the fibrils within cartilage tissue during external loading. Fibril diameter may relate directly to the mechanical properties of the tissue [171–173]. The length of the fibril also hypothetically affects the arrangement and reorientation of the fibrillar matrix during loading.
Past experimental measurements and theoretical predictions of the fluid pressurization of cartilage have demonstrated that the load support provided by interstitial fluid approaches the applied load immediately after loading, while recovering to zero over time during static loading [23, 24, 83, 85]. This excess fluid pressurization causes an issue with cartilage lubrication modifying the contact friction . This aspect has been poorly incorporated into the presently applied (fibril reinforced) models, although a theory applied by Ateshian et al.  has been shown to agree well with experimental measurements .
2.2. Cell-Level Models and Cell-Tissue Interactions
Biphasic multiscale models were first introduced in early 2000 . The goal of these models was to evaluate cell behavior in articular cartilage, cell-tissue interactions, and the effect of altered ECM and PCM properties for the deformation behavior of cells. The models have been applied to simulate cell-matrix interactions under steady state  and transient loading [59, 177, 178]. In the multiscale models, those parameters of the macroscopic tissue model that associate with the global loading (displacement, fluid pressure) serve as input parameters for the local submodel, including the chondrocyte and its local mechanical environment (Figure 3).
Recently, these multiscale models have incorporated the fibril reinforced constitutive laws of cartilage tissue presented in Section 2.1 [59, 132]. In addition to the ECM, the PCM around cells has been modeled. The PCM has also been modeled as fibril reinforced biphasic tissue and the fibrils have been implemented according to microscopic studies for tissue structure [132, 133]. This is because the loading-induced changes in the fibrillar patterns surrounding the chondrocytes likely transmit information from the matrix to cells causing metabolic changes in the cells . Hence, incorporation of fibril reinforced models to study mechanical behavior of chondrocytes is appropriate. Chondrocytes have been typically considered using the same constitutive equations as the ECM and PCM, except that the fibrils have been neglected. See the typical material properties of the ECM, PCM, and cell typically implemented in multiscale models of cartilage (Table 2).
Fibril reinforced multiscale models of cartilage have emphasized the importance of cartilage structure, composition, and mechanical properties on cell responses [16, 59, 132]. Specifically, the arcade-like orientation of the collagen network in the ECM has been suggested to modulate the cell shape differently in different zones . The fixed charge density of the ECM has been shown to increase cell aspect ratio especially in the superficial zone, while the fluid has been suggested to alter the transient behavior of cells . The collagen fibrils in the PCM have been found to modulate the mechanical signals experienced by the cells and protect cells, while the fixed charge density in the PCM may significantly alter cell morphology and deformation behavior [58, 133].
2.3. Adaptation Models
Mechanobiology investigates the biological responses of the cell to physical loading. Hence, it involves understanding how the tissue adapts the geometry, composition, and structure to the mechanical loading it is exposed to. With respect to articular cartilage, it is believed that the specific anisotropic organization of the collagen fibers is a result of the mechanical stimulation in the joint during development and growth . However, the rules by which the spatial organization of the collagen fibers is related to mechanical and chemical signals are not yet elucidated. With that in mind, algorithms that relate collagen remodeling and adaptation to stress and strain have been developed [179–181]. In the collagen remodeling algorithm, the collagen fibrils align along the preferred fibril directions that are situated between the positive principal strain directions [179, 180]. The simulation is conducted over time. After each step the collagen fibril orientation is updated and followed by an iteration where the new tissue strains are computed. This process is repeated until homeostasis is achieved. By combining the collagen remodeling algorithm with a fibril reinforced swelling model, Wilson et al. were able to predict the development of the typical Benninghoff-type collagen fiber orientation . Nagel and Kelly applied the remodeling model to observe effects in collagen architecture of cartilage in the vicinity of transplanted tissue .
Mechanobiological models have frequently been used to investigate changes in tissue composition and tissue volumetric growth in bone remodeling, fracture healing, and growth plate development [182–185]. Thus far, these models have not been extended to articular cartilage although van Donkelaar and Wilson simulated the effects of PG and collagen synthesis on chondrocyte hypertrophy . Such models may become useful in future studies addressing the mechanical influence on the progression of OA. In particular, combining mechanobiological adaptive models that evaluate tissue biosynthesis based on responses experienced by cells with imaging of the tissue geometry and composition may enable the development of patient-specific models to predict articular cartilage degradation and OA progression with time.
3. Practical Considerations in Modeling Applications
At present fibril reinforced models of articular cartilage are inhomogeneous and anisotropic, requiring several components to describe the overall mechanical behavior. These complex models are required for advanced simulations in order to mimic cartilage behavior at the same time considering the computational limitations. With a great number of tissue components and material parameters, the models become prone to uncertainties caused by the assumptions, for example, in the interactions between the model components. These uncertainties need to be addressed and the model implementation and function should be verified. With a carefully designed and verified computational model, experimental tests can be simulated quickly, and cases which are not possible in an experimental testing setup can be investigated, for example, through parametric analysis. Although a model may be plausible in theory, a thorough validation is required before the model can be applied. Subsequently, computational fibril reinforced models can be applied for retrieving intrinsic properties of the tissue, which describe its mechanical behavior by combining experimental tests and model simulations [7, 17, 47, 64, 153, 187–190]. In the following, some general considerations, also necessary when applying fibril reinforced materials of articular cartilage, are reviewed and discussed.
3.1. Parametric Analysis and Design of Experiments (DOE)
Sensitivity of different model components for the mechanical response of cartilage can be evaluated using parametric analysis. In general, the evaluation can be important either for simulations and experimental designs or for optimizing model behavior. Fibril reinforced models can have several material parameters describing stress of the collagen fibril matrix in addition to stress of the remaining solid matrix and fluid flow. To successfully predict the experimental data, the model should have a limited number of unknown parameter values in order to optimize the inverse problem uniquely. With a low number of parameters in optimization, the optimization routine becomes more reliable and efficient. In addition, the optimized variables are required to have distinguishable effects on the simulated outcome. To limit the number of optimized parameters, it is important to know which model parameters most sensitively affect the outcome of the simulations and which have a minor effect on the outcome. This way, the critical parameters can be optimized, while the parameters with lower impact can be assumed and fixed during optimization.
One common method of conducting parametric analysis is to vary one parameter at a time while keeping other parameters constant. This parametric analysis requires a basis, which is the set of reference parameter values which are then adjusted separately. The effect of these adjustments is then compared with the basis simulation. Several previous finite element analysis studies have used simple parametric analysis to either test the sensitivity of model parameters on the simulation outcome or to predict how the changes of model parameters would translate into observed changes in real-life physiology or the experiment outcome [20, 132, 191–193]. However, with such a simple method the interactions between the parameters cannot be captured or effectively predicted. Instead, a full factorial design could be used to test the effects of multiple parameters concurrently. In such case, the results of the parametric analysis are not affected by the basis simulation. Full factorial designs are suitable for studying systems with a low number of parameters and levels. They also provide more information on the interactions between the parameters. However, when the number of parameters and levels becomes higher, a full factorial design becomes unpractical due to the high number of required simulation runs.
One method for performing such multidimensional parametric analysis is the DOE approach using factorial analysis, which can also evaluate certain interactions between different model parameters . A full factorial design is the most reliable and comes with the greatest computational cost. In a full factorial design, all unknown parameters are evaluated as a function of each other. In the case of the more advanced models, it is clear that all of the models have a large amount of model parameters, some of which have values that have been assumed. For such cases, the full factorial design may become too extensive, and a more efficient way to determine the most significant parameters for the simulation outcome is required. One such approach is the fractional factorial design, which can significantly reduce the number of factorial test runs [195, 196]. Fractional factorial design is therefore suitable when determining the model parameters that affect the simulation outcome the most in models with a large number of independent parameters. This may be suitable for some of the more complex fibril reinforced models of cartilage . Different DOE methods should generally be used to reveal mechanisms that control the simulation outcome in advanced computational models [6, 183, 194, 197].
Although the applications of DOE for finite element analysis are numerous, it is still a quite new and seldom used technique in orthopaedic biomechanics [6, 183, 198–200]. The most popular design and method is the one by Taguchi et al. , which generally emphasizes the optimization of model performance by selection of the best values of the controllable parameters. We expect DOEs to gain more interest in the future also in finite element analysis of articular cartilage. The fibril reinforced models are becoming more complex and are being applied in more complex situations. Further, the need for optimized model design in material and geometry aspects becomes more important.
3.2. Considerations for Optimization
Curve fitting through optimization is an essential part of the application of fibril reinforced or any other biomechanical model, since material parameters describing model behavior can be derived to produce a simulation which corresponds with the experimental tests. These optimized material parameters then provide insight into the intrinsic tissue mechanical properties. Since there are various different constitutive models that can be used to simulate a similar outcome and agreement with certain experimental tests, it must be kept in mind that the optimized material parameters only describe the mechanical behavior of the tissue with the applied material model. Thus, the optimized material parameter values do not necessarily represent intrinsic tissue properties which can be generalized. Hence, the comparison of those material parameter values retrieved using fibril reinforced models should be conducted only between studies conducted using the same model. However, similarities in the overall principles between the different constitutive models may allow for comparison of the phenomena and trends observed between different studies and models.
As mentioned earlier, the modern models of articular cartilage include the description of many material parameters that cannot be measured directly. Due to the complexity of the fibril reinforced models, the fitting of the model to experimental data is often difficult through selection of model parameter values using the trial and error method. Therefore, optimization of parameter values against, for example, experimental data is conducted. During optimization, the model parameter values are often adjusted in a manner that will make the simulation result correspond with the experimental data as well as possible based on a defined objective function. In such a case, an objective function is minimized in case it represents an agreement error or maximized if the objective function represents the agreement between the experimental and simulated data [8, 12, 17–19, 47, 64, 65, 68, 98, 188, 189, 202–205]. Choice of the objective function is a major factor in determining the success of a fit. Most types of commercial simulation software, such as ABAQUS, COMSOL Multiphysics, provide a feasible platform to conduct modeling, and optimization can be done in conjunction with external software like MATLAB [189, 206] or Isight [207–209].
Prior to multidimensional optimization, it is suitable to run parametric analysis in order to choose only the sensitive model parameters for optimization. This will help to reach a unique solution. The optimized variables are required to have distinguishable effects on the simulated outcome for the solution to be unique. Generally, the lower the number of variables to be optimized, the more reliable and efficient the optimization routine. However, this may come with the cost of reduced agreement between the experimental and simulated data. The convergence of the optimization becomes quicker if the initial values are close to optimal . Therefore, initial values should not be selected in random, but within reasonable range, for instance based on literature source, or screening with parametric analysis [188, 210]. Randomly selected or unrealistic initial values could potentially cause unnecessary problems in the convergence of the model.
Considering that the experimental data includes some error, for example, due to uncertainties related to experimental testing, an error in the optimization should be accepted within that error margin. Therefore, a satisfactory optimization could result in various sets of optimized values depending on the initial values of the model parameters prior to optimization. However, changing the initial values and conducting the optimization again should result in a quite similar set of optimized parameter values, provided that the parameters chosen for optimization were considered to affect the simulation outcome sensitively. Therefore, it may be beneficial to use several initial value sets to verify uniqueness of fit and to determine the range at which the parameter values lie within satisfactory objective function values [68, 189, 202, 211].
Local minima are problems for the optimization (Figure 4), since they may result in fulfillment of some cutoff criteria for the optimization, providing one with a faulty set of model parameter values. A common practice in multidimensional optimization is to repeat the optimization routine in case of unsatisfactory objective function outcome (e.g., in a local minimum) by setting the optimum parameter value set as initial values for the following optimization (Figure 4). One must keep in mind that the chosen objective function has an effect on the optimized set of parameter values provided that complete agreement cannot be achieved . Hence, the most appropriate objective function should be selected. Generally, when comparing experimental and simulated data, an error function comparing the two for absolute agreement is considered suitable. Instead, if a trend between the experimental and simulation data is optimized, it may be appropriate to use correlation between the two as an objective function. One could apply weight functions for each data point in use when assessing the success of a fit. The weight functions could be based on standard deviation of the experimental data at a specific data point to consider the uncertainty in accuracy of the experimental data.
3.3. Convergence Tests
The rather complex geometries required to model the joint behavior are prone to ineffective computing, provided that the mesh is not optimized. To optimize the mesh density and biases after determining the simulation geometry, a convergence test should be performed to determine the coarsest acceptable mesh, using the element type that can produce an outcome with acceptable error, compared with the ideal; that is, the converged outcome. The localized anisotropies within the material model, like in the fibril reinforced models may affect the convergence, and optimal mesh design (Figure 5). In the test, mesh density is increased and the simulated outcome is recorded. As the simulation sequence reaches critical mesh density, that is, denser mesh will not change the simulation outcome, the mesh with an acceptable error (e.g., <5%) can be used to simulate the overall outcome. An example of such test is presented in Figure 5. It is also noteworthy that the CPU cost increases with denser finite element mesh, which is an extremely important time-consumption issue when the sample-specific properties are solved through optimization.
When a model is intended to simulate localized phenomena in a more complex geometry (e.g., complete knee joint) or in multiscale, it is necessary to use a mesh dense enough to have the model converge also for the localized simulations. Therefore, a similar convergence should be run to confirm the accuracy of the simulations specifically at regions of interest.
3.4. Model Validation
All models should be validated for feasibly corresponding with reality. The overall goal for validation is to make the model applicable to address a problem of interest with sufficient accuracy and confidence. In model calibration/optimization, the model should be compared with a real set of experimental data, and acceptable limits of discrepancies between the simulated and experimental data should be obtained [8, 12, 18, 64, 65, 78, 186, 212]. Considering fibril reinforced material models of cartilage tissue, validation of the model may come with some unique challenges, such as how to validate the fibril reinforcement in a mixture model and the interactions between the fibrillar and nonfibrillar matrices. Validation against experimental data can be made with three different approaches: direct, indirect, and trend measurements .
Direct experimental measurements and comparison with simulation results serve as the most reliable validation of the model. Further confidence can be gained through the direct parametric analysis of experimental tests and then comparing the results with those predicted by the simulations. A good example of a direct validation of a model theory to predict internal fluid pressurization was conducted by Soltz and Ateshian [83, 85] and Olsen et al. . However, direct validation is often difficult (e.g., pressure distribution in articular cartilage in a knee joint), since the simulated experiments cannot be replicated in controlled measurements with live subjects. For this reason, it may be required that indirect measurements are used instead. In such case, an experimentally measurable property can be compared with the simulation results, even though the measurable property is not needed. The measurable property must be selected so that it has a direct link with the model behavior being validated. However, indirect measurements validate the potential of a model to perform as required. Therefore, the importance of indirect validation should not be overestimated .
Sometimes the model is used for describing a certain kind of physiological behavior and it may not require direct or indirect validation for absolute agreement. In such case, trend measurements may be used to explain whether the model exhibits similar trends, for example, during parametric analysis. Trend measurements are especially important when validating the interactions between the model components (e.g., how the composition changes affect the measured mechanical response in articular cartilage). Although validation provides confidence in a model, it does not confirm that it agrees with reality. Each new model should be validated with the direct or indirect evidence of its function.
3.5. Experimental Tests
In the following subsection, we will discuss what factors should be taken into account when designing experimental mechanical tests at tissue and cell level in order to determine the material properties through optimization. We will also show how the models can be used to explain the experimentally observed phenomena.
3.5.1. Tissue-Level Experiments
With regard to the fibril reinforced model of cartilage, unique material parameters can be obtained provided that the role of each model parameter for the measured output parameter has been taken into account in the experimental design. A simple way of characterizing three parameters of the fibril reinforced poroelastic model (fibril network modulus, nonfibrillar matrix modulus, and permeability) is to fit the model to a stress-relaxation experiment under unconfined or indentation geometry (Figure 6). The fibril network modulus controls primarily the peak forces, the nonfibrillar matrix modulus modulates the equilibrium forces, and the permeability affects the rate of relaxation. Thus, all three parameters have their own distinct effect on the measured stress-relaxation response.
When adding strain-dependent parameters of the collagen fibrils and permeability in the model ((6) and (8)), the number of model parameters increases and the nonlinearity has to be taken into account for example by increasing the numbers of steps in stress-relaxation experiments (Figure 6). Other option would be the inclusion of different experiments, for example, in other loading geometries [8, 12, 98]. This way the nonlinearity caused by the added parameters can be captured by the model. As an example, permeability and nonfibrillar matrix properties could be fitted to confined compression data, after which the fibril properties could be fitted on unconfined compression and/or indentation experiments. If further swelling parameters are added in the models (16), swelling tests could be performed simply by recording cartilage swelling and shrinking in media with different osmolarities [214, 215].
3.5.2. Cell-Level Experiments
The cell-level biomechanical tests alone cannot be used for the optimization of model parameters. Instead, for model validation, experimental cell deformation behavior is compared with the corresponding computational analysis . The model parameters of the ECM for such model validation should come from tissue-level experiments. Ideally, in order to have sample-specific material parameters also for the PCM and cell, microscopic-scale experimental tests should be employed. The mechanical properties of the PCM have been characterized with microaspiration technique . The same technique can be used for measuring the mechanical properties of single cells . Atomic force microscopy with nanoindentation can be used to characterize the small-scale biomechanical properties of the ECM and PCM, and it has also been used for the characterization of cell properties .
In order to test cell deformation within cartilage tissue, osmotic and mechanical loading of the tissue have been utilized simultaneously, concurrently recording the cells using confocal or dual-photon laser scanning microscopy [26, 55, 220]. In the osmotic loading experiments, diffusion of fluid in or out of the tissue is carried out by altering the osmotic environment of the immersion media. The hypotonic medium increases the osmotic pressure difference in and out of the tissue, creating tissue swelling. This causes swelling of cells as well. Cells have, however, the ability to recover their original state after the initial swelling . This rapid recovery can be specifically seen using isolated cells. Some recent studies suggest that this recovery may be partly prevented or delayed by the ECM and/or PCM in the intact tissue , while chondrocytes in tissue explants and collagenase-degraded samples recovered back the original volume (Figure 7). Specific reasons for this phenomenon are, however, not known and computational modeling could provide an answer to such a problem.
Mechanical loading of articular cartilage, combined with simultaneous imaging of cells using confocal microscope, has been applied in unconfined compression and indentation geometries [26, 28]. Indentation combined with cell imaging can be used to characterize cell volume and morphology during compression. In this test, cartilage is intact, providing the physiological condition for cells. With this technique cell volumes were found to increase as a result of the mechanical loading of osteoarthritic rabbit cartilage, while the cells in normal joint cartilage were reduced in volume (Figure 7) [26, 27]. This technique, even using dual-photon microscope, can only reach some hundreds of microns into the tissue from the cartilage surface, limiting scanning to the superficial tissue layer. A technique in which chondrocytes are measured through the cut surface of tissue explants simultaneously with the compression of the cartilage surface can provide a way to characterize cell responses in the superficial, middle, and deep zones . A limitation of this technique may be that cutting of the samples can damage the integrity of the samples and loosening of the collagen fibril tension around cells. This may eventually lead to different responses of cells.
4. Visions for the Future
4.1. New Challenges
With the present cartilage models, supported by their further development, the new challenges for the future can be expected to advance from the generalization of theories and disease etiology more towards patient- and sample-specific simulations. At the same time, onset and progression of OA may be more specifically considered at a single-patient level as well as at a general level. For future applications, it is important to be able to predict how cartilage will react to changes in its mechanical environment. Important questions will arise, including: will cartilage be damaged and when, how will the damage progress, and can we predict the outflow of PGs? For future model development, it will be important to incorporate the production of the cartilage matrix  and incorporate tissue differentiation into cartilage models.
4.2. Combining In Vivo Imaging with Modeling
In order to be able to evaluate the cell responses of a patient with a model, joint and tissue forces need to be obtained. In order to obtain the most realistic evaluation of those forces, two key points need to be implemented in the models in a patient-specific manner: geometry and material properties. The former can be done using for instance MRI (Figure 8). The latter must be taken from the literature because the mechanical properties cannot be measured from patients, even though some studies have characterized the relationships between the MRI parameters and biomechanical properties of cartilage [47, 221]. However, cartilage structure, directly related to the mechanical properties of the tissue, can be estimated from in vivo imaging techniques. T2 relaxation time mapping serves as a technique to quantify the arrangement of the collagenous network in cartilage [40, 222], and the effect of T2-based collagen orientation on stresses and strains of cartilage in a knee joint was modeled recently . Diffusion tensor imaging may also be used to evaluate the collagenous structure of cartilage in undeformed and deformed state [63, 127, 129]. On the other hand, Gd-DTPA2−-enhanced T1 imaging (dGEMRIC) assesses the spatial changes in tissue PG content . Contrast-enhanced cartilage tomography can be used to characterize the fixed charge density distribution of cartilage .
4.3. Toward Patient-Specific Estimation of Disease Progression and Treatment
The models presented in this paper and future adaptations of new fibril reinforced models could lead towards patient-specific estimation of the loading effects on cartilage degeneration in OA. By incorporating adaptive algorithms and phenomena behind tissue formation and degradation in the models, and combining these models with in vivo imaging and human movement analysis (Figure 8), it may be possible to develop computational models of entire joints. These models could be used, for example, to estimate the development of OA in patients with joint injuries (e.g., cartilage damage, meniscus tear, and ligament injury).
Some computer simulation work has already been conducted with fibril reinforced material models of articular cartilage to understand and perhaps identify risk factors of OA development and progression [53, 115–117, 225]. In the future, modeling could be used for example for the evaluation of the effect of different clinical operations on the onset and progression of OA, in the first stage through stresses and strains of cartilage and in the later stages through adaptation models [179, 226, 227]. Modeling could also be used to evaluate the effect of different conservative strategies (for instance change in exercise) on the disease progression. Thus, the models could be used in the clinical decision making. Quantitative evaluation and prediction of the joint condition for the future would certainly aid in choosing the best treatment for the disease. When combined with potential repair materials, the computational models may have potential for optimizing and preparing growth protocols for those materials mimicking cartilage . Furthermore, computational models should aid in the design of loading protocols for joint rehabilitation with degenerated cartilage, for example, when transplanting repair tissues that lack the mechanical requirements of the host tissue . Recently Khoshgoftar et al.  investigated the potential of using a numerical model for stimulating formation of physiological collagen architecture in tissue-engineered cartilage. Such applications combined with the most modern models of tissue adaptation will certainly help to utilize mechanical loading regimes to benefit the optimizing of tissue-engineered cartilage .
The research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP/2007-2013)/ERC Grant Agreement no. 281180, and from the Academy of Finland (Grant 218038). Research has been supported by the strategic funding of the University of Eastern Finland. Also, support was obtained from the AI-HS Team grant on Osteoarthritis, CIHR, and the Killam Foundation. The CSC-IT Center for Science is acknowledged for technical support.
- P. Julkunen, J. Iivarinen, P. A. Brama, J. Arokoski, J. S. Jurvelin, and H. J. Helminen, “Maturation of collagen fibril network structure in tibial and femoral cartilage of rabbits,” Osteoarthritis and Cartilage, vol. 18, no. 3, pp. 406–415, 2010.
- W. Gründer, “MRI assessment of cartilage ultrastructure,” NMR in Biomedicine, vol. 19, no. 7, pp. 855–876, 2006.
- J. Rieppo, M. M. Hyttinen, E. Halmesmaki et al., “Changes in spatial collagen content and collagen network architecture in porcine articular cartilage during growth and maturation,” Osteoarthritis and Cartilage, vol. 17, no. 4, pp. 448–455, 2009.
- M. C. van Turnhout, H. Schipper, B. Engel, W. Buist, S. Kranenbarg, and J. L. van Leeuwen, “Postnatal development of collagen structure in ovine articular cartilage,” BMC Developmental Biology, vol. 10, article 62, 2010.
- A. Benninghoff, “Form und Bau der Gelenkknorpel in ihren Beziehungen zur Funktion,” Zeitschrift für Zellforschung und Mikroskopische Anatomie, vol. 2, no. 5, pp. 783–862, 1925.
- P. Julkunen, J. S. Jurvelin, and H. Isaksson, “Contribution of tissue composition and structure to mechanical response of articular cartilage under different loading geometries and strain rates,” Biomechanics and Modeling in Mechanobiology, vol. 9, no. 2, pp. 237–245, 2010.
- R. K. Korhonen, M. S. Laasanen, J. Töyräs, R. Lappalainen, H. J. Helminen, and J. S. Jurvelin, “Fibril reinforced poroelastic model predicts specifically mechanical behavior of normal, proteoglycan depleted and collagen degraded articular cartilage,” Journal of Biomechanics, vol. 36, no. 9, pp. 1373–1379, 2003.
- W. Wilson, J. M. Huyghe, and C. C. van Donkelaar, “A composition-based cartilage model for the assessment of compositional changes during cartilage damage and adaptation,” Osteoarthritis and Cartilage, vol. 14, no. 6, pp. 554–560, 2006.
- N. D. Broom, “Further insights into the structural principles governing the function of articular cartilage,” Journal of Anatomy, vol. 139, no. 2, pp. 275–294, 1984.
- N. D. Broom and D. B. Myers, “A study of the structural response of wet hyaline cartilage to various loading situations,” Connective Tissue Research, vol. 7, no. 4, pp. 227–237, 1980.
- J. P. G. Urban, A. Maroudas, M. T. Bayliss, and J. Dillon, “Swelling pressures of proteoglycans at the concentrations found in cartilaginous tissues,” Biorheology, vol. 16, no. 6, pp. 447–464, 1979.
- W. Wilson, C. C. van Donkelaar, B. van Rietbergen, and R. Huiskes, “A fibril-reinforced poroviscoelastic swelling model for articular cartilage,” Journal of Biomechanics, vol. 38, no. 6, pp. 1195–1204, 2005.
- V. C. Mow, W. Zhu, and A. Ratcliffe, “Structure and function of articular cartilage and meniscus,” in Basic Orthopaedic Biomechanics, V. C. Mow and W. C. Hayes, Eds., pp. 143–198, Raven Press, New York, NY, USA, 1991.
- A. Maroudas, “Balance between swelling pressure and collagen tension in normal and degenerate cartilage,” Nature, vol. 260, no. 5554, pp. 808–809, 1976.
- P. Julkunen, E. P. Halmesmäki, J. Iivarinen et al., “Effects of growth and exercise on composition, structural maturation and appearance of osteoarthritis in articular cartilage of hamsters,” Journal of Anatomy, vol. 217, no. 3, pp. 262–274, 2010.
- R. K. Korhonen, P. Julkunen, J. S. Jurvelin, and S. Saarakkala, “Structural and compositional changes in peri- and extracellular matrix of osteoarthritic cartilage modulate chondrocyte morphology,” Cellular and Molecular Bioengineering, vol. 4, no. 3, pp. 484–494, 2011.
- P. Julkunen, P. Kiviranta, W. Wilson, J. S. Jurvelin, and R. K. Korhonen, “Characterization of articular cartilage by combining microscopic analysis with a fibril-reinforced finite-element model,” Journal of Biomechanics, vol. 40, no. 8, pp. 1862–1870, 2007.
- W. Wilson, C. C. van Donkelaar, B. van Rietbergen, K. Ito, and R. Huiskes, “Stresses in the local collagen network of articular cartilage: a poroviscoelastic fibril-reinforced finite element study,” Journal of Biomechanics, vol. 37, no. 3, pp. 357–366, 2004.
- P. Julkunen, T. Harjula, J. Marjanen, H. J. Helminen, and J. S. Jurvelin, “Comparison of single-phase isotropic elastic and fibril-reinforced poroelastic models for indentation of rabbit articular cartilage,” Journal of Biomechanics, vol. 42, no. 5, pp. 652–656, 2009.
- P. Julkunen, R. K. Korhonen, W. Herzog, and J. S. Jurvelin, “Uncertainties in indentation testing of articular cartilage: a fibril-reinforced poroviscoelastic study,” Medical Engineering and Physics, vol. 30, no. 4, pp. 506–515, 2008.
- M. S. Laasanen, J. Töyräs, R. K. Korhonen et al., “Biomechanical properties of knee articular cartilage,” Biorheology, vol. 40, no. 1–3, pp. 133–140, 2002.
- D. L. Bader, G. E. Kempson, J. Egan, W. Gilbey, and A. J. Barrett, “The effects of selective matrix degradation on the short-term compressive properties of adult human articular cartilage,” Biochimica et Biophysica Acta, vol. 1116, no. 2, pp. 147–154, 1992.
- A. Oloyede and N. D. Broom, “Is classical consolidation theory applicable to articular cartilage deformation?” Clinical Biomechanics, vol. 6, no. 4, pp. 206–212, 1991.
- A. Oloyede and N. Broom, “Stress-sharing between the fluid and solid components of articular cartilage under varying rates of compression,” Connective Tissue Research, vol. 30, no. 2, pp. 127–141, 1993.
- A. Oloyede and N. D. Broom, “Complex nature of stress inside loaded articular cartilage,” Clinical Biomechanics, vol. 9, no. 3, pp. 149–156, 1994.
- S. K. Han, R. Seerattan, and W. Herzog, “Mechanical loading of in situ chondrocytes in lapine retropatellar cartilage after anterior cruciate ligament transection,” Journal of the Royal Society Interface, vol. 7, no. 47, pp. 895–903, 2010.
- S. M. Turunen, S. K. Han, W. Herzog, and R. K. Korhonen, “Cell deformation behavior in mechanically loaded rabbit articular cartilage 4 weeks after anterior cruciate ligament transection,” Osteoarthritis and Cartilage, vol. 21, no. 3, pp. 505–513, 2012.
- F. Guilak, A. Ratcliffe, and V. C. Mow, “Chondrocyte deformation and local tissue strain in articular cartilage: a confocal microscopy study,” Journal of Orthopaedic Research, vol. 13, no. 3, pp. 410–421, 1995.
- M. D. Buschmann, Y. A. Gluzband, A. J. Grodzinsky, J. H. Kimura, and E. B. Hunziker, “Chondrocytes in agarose culture synthesize a mechanically functional extracellular matrix,” Journal of Orthopaedic Research, vol. 10, no. 6, pp. 745–758, 1992.
- F. Guilak, L. G. Alexopoulos, M. L. Upton et al., “The pericellular matrix as a transducer of biomechanical and biochemical signals in articular cartilage,” Annals of the New York Academy of Sciences, vol. 1068, no. 1, pp. 498–512, 2006.
- R. A. Stockwell, “Cartilage failure in osteoarthritis: relevance of normal structure and function. A review,” Clinical Anatomy, vol. 4, no. 3, pp. 161–191, 1991.
- S. Saarakkala, P. Julkunen, P. Kiviranta, J. Mäkitalo, J. S. Jurvelin, and R. K. Korhonen, “Depth-wise progression of osteoarthritis in human articular cartilage: investigation of composition, structure and biomechanics,” Osteoarthritis and Cartilage, vol. 18, no. 1, pp. 73–81, 2010.
- J. P. A. Arokoski, J. S. Jurvelin, U. Väätäinen, and H. J. Helminen, “Normal and pathological adaptations of articular cartilage to joint loading,” Scandinavian Journal of Medicine and Science in Sports, vol. 10, no. 4, pp. 186–198, 2000.
- X. Bi, X. Yang, M. P. G. Bostrom, and N. P. Camacho, “Fourier transform infrared imaging spectroscopy investigations in the pathogenesis and repair of cartilage,” Biochimica et Biophysica Acta, vol. 1758, no. 7, pp. 934–941, 2006.
- J. A. Buckwalter and H. J. Mankin, “Articular cartilage. Part II: degeneration and osteoarthrosis, repair, regeneration, and transplantation,” Journal of Bone and Joint Surgery. American, vol. 79, no. 4, pp. 612–632, 1997.
- N. D. Broom, T. Ngo, and E. Tham, “Traversing the intact/fibrillated joint surface: a biomechanical interpretation,” Journal of Anatomy, vol. 206, no. 1, pp. 55–67, 2005.
- N. D. Broom, “The collagenous architecture of articular cartilage—a synthesis of ultrastructure and mechanical function,” Journal of Rheumatology, vol. 13, no. 1, pp. 142–152, 1986.
- X. Bi, G. Li, S. B. Doty, and N. P. Camacho, “A novel method for determination of collagen orientation in cartilage by Fourier transform infrared imaging spectroscopy (FT-IRIS),” Osteoarthritis and Cartilage, vol. 13, no. 12, pp. 1050–1058, 2005.
- H. E. Panula, M. M. Hyttinen, J. P. A. Arokoski et al., “Articular cartilage superficial zone collagen birefringence reduced and cartilage thickness increased before surface fibrillation in experimental osteoarthritis,” Annals of the Rheumatic Diseases, vol. 57, no. 4, pp. 237–245, 1998.
- M. T. Nieminen, J. Rieppo, J. Töyräs et al., “T2 relaxation reveals spatial collagen architecture in articular cartilage: a comparative quantitative MRI and polarized light microscopic study,” Magnetic Resonance in Medicine, vol. 46, no. 3, pp. 487–493, 2001.
- X. Bi, X. Yang, M. P. G. Bostrom et al., “Fourier transform infrared imaging and MR microscopy studies detect compositional and structural changes in cartilage in a rabbit model of osteoarthritis,” Analytical and Bioanalytical Chemistry, vol. 387, no. 5, pp. 1601–1612, 2007.
- H. A. Alhadtaq, Y. Xia, J. B. Moody, and J. R. Matyas, “Detecting structural changes in early experimental osteoarthritis of tibial cartilage by microscopic magnetic resonance imaging and polarised light microscopy,” Annals of the Rheumatic Diseases, vol. 63, no. 6, pp. 709–717, 2004.
- M. L. Gray, D. Burstein, and Y. Xia, “Biochemical (and functional) imaging of articular cartilage,” Seminars in Musculoskeletal Radiology, vol. 5, no. 4, pp. 329–343, 2001.
- J. D. Rubenstein, J. K. Kim, I. Morava-Protzner, P. L. Stanchev, and R. M. Henkelman, “Effects of collagen orientation on MR imaging characteristics of bovine articular cartilage,” Radiology, vol. 188, no. 1, pp. 219–226, 1993.
- A. Bashir, M. L. Gray, and D. Burstein, “Gd-DTPA2- as a measure of cartilage degradation,” Magnetic Resonance in Medicine, vol. 36, no. 5, pp. 665–673, 1996.
- A. Bashir, M. L. Gray, J. Hartke, and D. Burstein, “Nondestructive imaging of human cartilage glycosaminoglycan concentration by MRI,” Magnetic Resonance in Medicine, vol. 41, no. 5, pp. 857–865, 1999.
- P. Julkunen, R. K. Korhonen, M. J. Nissi, and J. S. Jurvelin, “Mechanical characterization of articular cartilage by combining magnetic resonance imaging and finite-element analysis—a potential functional imaging technique,” Physics in Medicine and Biology, vol. 53, no. 9, pp. 2425–2438, 2008.
- P. Kiviranta, E. Lammentausta, J. Töyräs, I. Kiviranta, and J. S. Jurvelin, “Indentation diagnostics of cartilage degeneration,” Osteoarthritis and Cartilage, vol. 16, no. 7, pp. 796–804, 2008.
- R. L. Sah, A. S. Yang, A. C. Chen et al., “Physical properties of rabbit articular cartilage after transection of the anterior cruciate ligament,” Journal of Orthopaedic Research, vol. 15, no. 2, pp. 197–203, 1997.
- J. S. Wayne, K. A. Kraft, K. J. Shields, C. Yin, J. R. Owen, and D. G. Disler, “MR imaging of normal and matrix-depleted cartilage: correlation with biomechanical function and biochemical composition,” Radiology, vol. 228, no. 2, pp. 493–499, 2003.
- J. R. Parsons and J. Black, “Mechanical behavior of articular cartilage quantitative changes with enzymatic alteration of the proteoglycan fraction,” Bulletin of the Hospital for Joint Diseases Orthopaedic Institute, vol. 47, no. 1, pp. 13–30, 1987.
- T. Franz, E. M. Hasler, R. Hagg, C. Weiler, R. P. Jakob, and P. Mainil-Varlet, “In situ compressive stiffness, biochemical composition, and structural integrity of articular cartilage of the human knee joint,” Osteoarthritis and Cartilage, vol. 9, no. 6, pp. 582–592, 2001.
- M. E. Mononen, M. T. Mikkola, P. Julkunen et al., “Effect of superficial collagen patterns and fibrillation of femoral articular cartilage on knee joint mechanics-a 3D finite element analysis,” Journal of Biomechanics, vol. 45, no. 3, pp. 579–587, 2012.
- P. G. Bush and A. C. Hall, “Regulatory volume decrease (RVD) by isolated and in situ bovine articular chondrocytes,” Journal of Cellular Physiology, vol. 187, no. 3, pp. 304–314, 2001.
- S. M. Turunen, M. J. Lammi, S. Saarakkala, A. Koistinen, and R. K. Korhonen, “Hypotonic challenge modulates cell volumes differently in the superficial zone of intact articular cartilage and cartilage explant,” Biomechanics and Modeling in Mechanobiology, vol. 11, no. 5, pp. 665–675, 2012.
- S. K. Han, P. Colarusso, and W. Herzog, “Confocal microscopy indentation system for studying in situ chondrocyte mechanics,” Medical Engineering and Physics, vol. 31, no. 8, pp. 1038–1042, 2009.
- S. K. Han, W. Wouters, A. Clark, and W. Herzog, “Mechanically induced calcium signaling in chondrocytes in situ,” Journal of Orthopaedic Research, vol. 30, no. 3, pp. 475–481, 2012.
- M. Likhitpanichkul, X. E. Guo, and V. C. Mow, “The effect of matrix tension-compression nonlinearity and fixed negative charges on chondrocyte responses in cartilage,” MCB Molecular and Cellular Biomechanics, vol. 2, no. 4, pp. 191–204, 2005.
- R. K. Korhonen, P. Julkunen, W. Wilson, and W. Herzog, “Importance of collagen orientation and depth-dependent fixed charge densities of cartilage on mechanical behavior of chondrocytes,” Journal of Biomechanical Engineering, vol. 130, no. 2, Article ID 021003, 2008.
- M. C. van Turnhout, S. Kranenbarg, and J. L. van Leeuwen, “Contribution of postnatal collagen reorientation to depth-dependent mechanical properties of articular cartilage,” Biomechanics and Modeling in Mechanobiology, vol. 10, no. 2, pp. 269–279, 2011.
- W. Wilson, J. M. Huyghe, and C. C. van Donkelaar, “Depth-dependent compressive equilibrium properties of articular cartilage explained by its composition,” Biomechanics and Modeling in Mechanobiology, vol. 6, no. 1-2, pp. 43–53, 2007.
- L. P. Räsänen, M. E. Mononen, M. T. Nieminen, E. Lammentausta, J. S. Jurvelin, and R. K. Korhonen, “Implementation of subject-specific collagen architecture of cartilage into a 2D computational model of a knee joint—data from the Osteoarthritis Initiative (OAI),” Journal of Orthopaedic Research, vol. 31, no. 1, pp. 10–22, 2012.
- D. M. Pierce, W. Trobin, J. G. Raya et al., “DT-MRI based computation of collagen fiber deformation in human articular cartilage: a feasibility study,” Annals of Biomedical Engineering, vol. 38, no. 7, pp. 2447–2463, 2010.
- P. Julkunen, W. Wilson, J. S. Jurvelin et al., “Stress-relaxation of human patellar articular cartilage in unconfined compression: prediction of mechanical response by tissue composition and structure,” Journal of Biomechanics, vol. 41, no. 9, pp. 1978–1986, 2008.
- J. J. García and D. H. Cortés, “A biphasic viscohyperelastic fibril-reinforced model for articular cartilage: formulation and comparison with experimental data,” Journal of Biomechanics, vol. 40, no. 8, pp. 1737–1744, 2007.
- G. A. Ateshian, V. Rajan, N. O. Chahine, C. E. Canal, and C. T. Hung, “Modeling the matrix of articular cartilage using a continuous fiber angular distribution predicts many observed phenomena,” Journal of Biomechanical Engineering, vol. 131, no. 6, 2009.
- D. M. Pierce, T. Ricken, and G. A. Holzapfel, “A hyperelastic biphasic fibre-reinforced model of articular cartilage considering distributed collagen fibre orientations: continuum basis, computational aspects and applications,” Computer Methods in Biomechanics and Biomedical Engineering, 2012.
- A. Seifzadeh, J. Wang, D. C. Oguamanam, and M. Papini, “A nonlinear biphasic fiber-reinforced porohyperviscoelastic model of articular cartilage incorporating fiber reorientation and dispersion,” Journal of Biomechanical Engineering, vol. 133, no. 8, Article ID 081004, 2011.
- R. K. Korhonen, M. Wong, J. Arokoski et al., “Importance of the superficial tissue layer for the indentation stiffness of articular cartilage,” Medical Engineering and Physics, vol. 24, no. 2, pp. 99–108, 2002.
- J. Rieppo, J. Töyräs, M. T. Nieminen et al., “Structure-function relationships in enzymatically modified articular cartilage,” Cells Tissues Organs, vol. 175, no. 3, pp. 121–132, 2003.
- A. Oloyede and N. D. Broom, “A physical model for the time-dependent deformation of articular cartilage,” Connective Tissue Research, vol. 29, no. 4, pp. 251–261, 1993.
- J. R. Parsons and J. Black, “The viscoelastic shear behavior of normal rabbit articular cartilage,” Journal of Biomechanics, vol. 10, no. 1, pp. 21–29, 1977.
- J. R. Parsons and J. Black, “Mechanical behavior of articular cartilage: quantitative changes with alteration of ionic environment,” Journal of Biomechanics, vol. 12, no. 10, pp. 765–773, 1979.
- W. C. Hayes and A. J. Bodine, “Flow-independent viscoelastic properties of articular cartilage matrix,” Journal of Biomechanics, vol. 11, no. 8-9, pp. 407–419, 1978.
- P. A. Torzilli and V. C. Mow, “On the fundamental fluid transport mechanisms through normal and pathological articular cartilage during function. II. The analysis, solution and conclusions,” Journal of Biomechanics, vol. 9, no. 9, pp. 587–606, 1976.
- P. A. Torzilli and V. C. Mow, “On the fundamental fluid transport mechanisms through normal and pathological articular cartilage during function. I the formulation,” Journal of Biomechanics, vol. 9, no. 8, pp. 541–552, 1976.
- V. C. Mow and J. M. Mansour, “The nonlinear interaction between cartilage deformation and interstitial fluid flow,” Journal of Biomechanics, vol. 10, no. 1, pp. 31–39, 1977.
- V. C. Mow, S. C. Kuei, W. M. Lai, and C. G. Armstrong, “Biphasic creep and stress relaxation of articular cartilage in compression: theory and experiments,” Journal of Biomechanical Engineering, vol. 102, no. 1, pp. 73–84, 1980.
- J. M. Mansour and V. C. Mow, “The permeability of articular cartilage under compressive strain and at high pressures,” Journal of Bone and Joint Surgery. American, vol. 58, no. 4, pp. 509–516, 1976.
- C. W. McCutchen, “Cartilage is poroelastic, not viscoelastic (including an exact theorem about strain energy and viscous loss, and an order of magnitude relation for equilibration time,” Journal of Biomechanics, vol. 15, no. 4, pp. 325–327, 1982.
- L. Li and W. Herzog, “The role of viscoelasticity of collagen fibers in articular cartilage: theory and numerical formulation,” Biorheology, vol. 41, no. 3-4, pp. 181–194, 2004.
- T. Macirowski, S. Tepic, and R. W. Mann, “Cartilage stresses in the human hip joint,” Journal of Biomechanical Engineering, vol. 116, no. 1, pp. 10–18, 1994.
- M. A. Soltz and G. A. Ateshian, “Experimental verification and theoretical prediction of cartilage interstitial fluid pressurization at an impermeable contact interface in confined compression,” Journal of Biomechanics, vol. 31, no. 10, pp. 927–934, 1998.
- L. A. Setton, W. Zhu, and V. C. Mow, “The biphasic poroviscoelastic behavior of articular cartilage: role of the surface zone in governing the compressive behavior,” Journal of Biomechanics, vol. 26, no. 4-5, pp. 581–592, 1993.
- M. A. Soltz and G. A. Ateshian, “A conewise linear elasticity mixture model for the analysis of tension-compression nonlinearity in articular cartilage,” Journal of Biomechanical Engineering, vol. 122, no. 6, pp. 576–586, 2000.
- T. D. Brown and R. J. Singerman, “Experimental determination of the linear biphasic constitutive coefficients of human fetal proximal femoral chondroepiphysis,” Journal of Biomechanics, vol. 19, no. 8, pp. 597–605, 1986.
- C. G. Armstrong, W. M. Lai, and V. C. Mow, “An analysis of the unconfined compression of articular cartilage,” Journal of Biomechanical Engineering, vol. 106, no. 2, pp. 165–173, 1984.
- B. Cohen, W. M. Lai, and V. C. Mow, “A transversely isotropic biphasic model for unconfined compression of growth plate and chondroepiphysis,” Journal of Biomechanical Engineering, vol. 120, no. 4, pp. 491–496, 1998.
- P. M. Bursać, T. W. Obitz, S. R. Eisenberg, and D. Stamenović, “Confined and unconfined stress relaxation of cartilage: appropriateness of a transversely isotropic analysis,” Journal of Biomechanics, vol. 32, no. 10, pp. 1125–1130, 1999.
- V. C. Mow, D. C. Fithian, and M. A. Kelly, “Fundamentals of articular cartilage and meniscus biomechanics,” in Articular Cartilage and Knee Joint Function: Basic Science and Arthroscopy, J. W. Ewing, Ed., pp. 1–18, Raven Press, New York, NY, USA, 1990.
- F. Boschetti and G. M. Peretti, “Tensile and compressive properties of healthy and osteoarthritic human articular cartilage,” Biorheology, vol. 45, no. 3-4, pp. 337–344, 2008.
- J. Soulhat, M. D. Buschmann, and A. Shirazi-Adl, “A fibril-network-reinforced biphasic model of cartilage in unconfined compression,” Journal of Biomechanical Engineering, vol. 121, no. 3, pp. 340–347, 1999.
- D. K. Miller, “Technical note: modelling soft tissue using biphasic theory—a word of caution,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 1, no. 3, pp. 261–263, 1998.
- C. Y. Huang, V. C. Mow, and G. A. Ateshian, “The role of flow-independent viscoelasticity in the biphasic tensile and compressive responses of articular cartilage,” Journal of Biomechanical Engineering, vol. 123, no. 5, pp. 410–417, 2001.
- C. Y. Huang, M. A. Soltz, M. Kopacz, V. C. Mow, and G. A. Ateshian, “Experimental verification of the roles of intrinsic matrix viscoelasticity and tension-compression nonlinearity in the biphasic response of cartilage,” Journal of Biomechanical Engineering, vol. 125, no. 1, pp. 84–93, 2003.
- R. K. June, S. Ly, and D. P. Fyhrie, “Cartilage stress-relaxation proceeds slower at higher compressive strains,” Archives of Biochemistry and Biophysics, vol. 483, no. 1, pp. 75–80, 2009.
- W. Zhu, V. C. Mow, T. J. Koob, and D. R. Eyre, “Viscoelastic shear properties of articular cartilage and the effects of glycosidase treatments,” Journal of Orthopaedic Research, vol. 11, no. 6, pp. 771–781, 1993.
- M. R. DiSilvestro and J. K. F. Suh, “A cross-validation of the biphasic poroviscoelastic model of articular cartilage in unconfined compression, indentation, and confined compression,” Journal of Biomechanics, vol. 34, no. 4, pp. 519–525, 2001.
- M. R. DiSilvestro, Q. Zhu, M. Wong, J. S. Jurvelin, and J. K. F. Suh, “Biphasic poroviscoelastic simulation of the unconfined compression of articular cartilage: I—simultaneous prediction of reaction force and lateral displacement,” Journal of Biomechanical Engineering, vol. 123, no. 2, pp. 191–197, 2001.
- D. F. Betsch and E. Baer, “Structure and mechanical properties of rat tail tendon,” Biorheology, vol. 17, no. 1-2, pp. 83–94, 1980.
- J. L. Wang, M. Parnianpour, A. Shirazi-Adl, and A. E. Engin, “Failure criterion of collagen fiber: viscoelastic behavior simulated by using load control data,” Theoretical and Applied Fracture Mechanics, vol. 27, no. 1, pp. 1–12, 1997.
- R. Sanjeevi, N. Somanathan, and D. Ramaswamy, “Viscoelastic model for collagen fibres,” Journal of Biomechanics, vol. 15, no. 3, pp. 181–183, 1982.
- L. Han, E. H. Frank, J. J. Greene et al., “Time-dependent nanomechanics of cartilage,” Biophysical Journal, vol. 100, no. 7, pp. 1846–1854, 2011.
- M. B. Schmidt, V. C. Mow, L. E. Chun, and D. R. Eyre, “Effects of proteoglycan extraction on the tensile behavior of articular cartilage,” Journal of Orthopaedic Research, vol. 8, no. 3, pp. 353–363, 1990.
- W. Zhu, J. C. Iatridis, V. Hlibczuk, A. Ratcliffe, and V. C. Mow, “Determination of collagen-proteoglycan interactions in vitro,” Journal of Biomechanics, vol. 29, no. 6, pp. 773–783, 1996.
- A. F. Mak, “The apparent viscoelastic behavior of articular cartilage—the contributions from the intrinsic matrix viscoelasticity and interstitial fluid flows,” Journal of Biomechanical Engineering, vol. 108, no. 2, pp. 123–130, 1986.
- Y. C. Fung, Biomechanics: Mechanical Properties of Living Tissues, Springer, New York, NY, USA, 1981.
- D. M. Pierce, W. Trobin, S. Trattnig, H. Bischof, and G. A. Holzapfel, “A phenomenological approach toward patient-specific computational modeling of articular cartilage including collagen fiber tracking,” Journal of Biomechanical Engineering, vol. 131, no. 9, Article ID 091006, 2009.
- W. Mesfar and A. Shirazi-Adl, “Biomechanics of the knee joint in flexion under various quadriceps forces,” Knee, vol. 12, no. 6, pp. 424–434, 2005.
- C. K. Fitzpatrick and P. J. Rullkoetter, “Influence of patellofemoral articular geometry and material on mechanics of the unresurfaced patella,” Journal of Biomechanics, vol. 45, no. 11, pp. 1909–1915, 2012.
- A. Homyk, A. Orsi, S. Wibby, N. Yang, H. Nayeb-Hashemi, and P. K. Canavan, “Failure locus of the anterior cruciate ligament: 3D finite element analysis,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 15, no. 8, pp. 865–874, 2012.
- M. Adouni and A. Shirazi-Adl, “Knee joint biomechanics in closed-kinetic-chain exercises,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 12, no. 6, pp. 661–670, 2009.
- E. Peña, B. Calvo, M. A. Martínez, and M. Doblaré, “Computer simulation of damage on distal femoral articular cartilage after meniscectomies,” Computers in Biology and Medicine, vol. 38, no. 1, pp. 69–81, 2008.
- M. Adouni, A. Shirazi-Adl, and R. Shirazi, “Computational biodynamics of human knee joint in gait: from muscle forces to cartilage stresses,” Journal of Biomechanics, vol. 45, no. 12, pp. 2149–2156, 2012.
- R. Shirazi and A. Shirazi-Adl, “Computational biomechanics of articular cartilage of human knee joint: effect of osteochondral defects,” Journal of Biomechanics, vol. 42, no. 15, pp. 2458–2465, 2009.
- R. Shirazi and A. Shirazi-Adl, “Analysis of partial meniscectomy and ACL reconstruction in knee joint biomechanics under a combined loading,” Clinical Biomechanics, vol. 24, no. 9, pp. 755–761, 2009.
- M. Kazemi, L. P. Li, M. D. Buschmann, and P. Savard, “Partial meniscectomy changes fluid pressurization in articular cartilage in human knees,” Journal of Biomechanical Engineering, vol. 134, no. 2, Article ID 021001, 2012.
- G. A. Ateshian, B. J. Ellis, and J. A. Weiss, “Equivalence between short-time biphasic and incompressible elastic material responses,” Journal of Biomechanical Engineering, vol. 129, no. 3, pp. 405–412, 2007.
- L. P. Li, J. Soulhat, M. D. Buschmann, and A. Shirazi-Adl, “Nonlinear analysis of cartilage in unconfined ramp compression using a fibril reinforced poroelastic model,” Clinical Biomechanics, vol. 14, no. 9, pp. 673–682, 1999.
- A. van der Voet, “A comparison of finite element codes for the solution of biphasic poroelastic problems,” Proceedings of the Institution of Mechanical Engineers. Part H, vol. 211, no. 2, pp. 209–211, 1997.
- R. Shirazi and A. Shirazi-Adl, “Deep vertical collagen fibrils play a significant role in mechanics of articular cartilage,” Journal of Orthopaedic Research, vol. 26, no. 5, pp. 608–615, 2008.
- L. P. Li, M. D. Buschmann, and A. Shirazi-Adl, “A fibril reinforced nonhomogeneous poroelastic model for articular cartilage: inhomogeneous response in unconfined compression,” Journal of Biomechanics, vol. 33, no. 12, pp. 1533–1541, 2000.
- L. P. Li, M. D. Buschmann, and A. Shirazi-Adl, “Strain-rate dependent stiffness of articular cartilage in unconfined compression,” Journal of Biomechanical Engineering, vol. 125, no. 2, pp. 161–168, 2003.
- S. Federico and W. Herzog, “On the anisotropy and inhomogeneity of permeability in articular cartilage,” Biomechanics and Modeling in Mechanobiology, vol. 7, no. 5, pp. 367–378, 2008.
- M. B. Lilledahl, D. M. Pierce, T. Ricken, G. A. Holzapfel, and C. de Lange Davies, “Structural analysis of articular cartilage using multiphoton microscopy: input for biomechanical modeling,” IEEE Transactions on Medical Imaging, vol. 30, no. 9, pp. 1635–1648, 2011.
- K. B. Gu and L. P. Li, “A human knee joint model considering fluid pressure and fiber orientation in cartilages and menisci,” Medical Engineering and Physics, vol. 33, no. 4, pp. 497–503, 2011.
- S. K. de Visser, J. C. Bowden, E. Wentrup-Byrne et al., “Anisotropy of collagen fibre alignment in bovine cartilage: comparison of polarised light microscopy and spatially resolved diffusion-tensor measurements,” Osteoarthritis and Cartilage, vol. 16, no. 6, pp. 689–697, 2008.
- L. Filidoro, O. Dietrich, J. Weber et al., “High-resolution diffusion tensor imaging of human patellar cartilage: feasibility and preliminary findings,” Magnetic Resonance in Medicine, vol. 53, no. 5, pp. 993–998, 2005.
- R. Meder, S. K. de Visser, J. C. Bowden, T. Bostrom, and J. M. Pope, “Diffusion tensor imaging of articular cartilage as a measure of tissue microstructure,” Osteoarthritis and Cartilage, vol. 14, no. 9, pp. 875–881, 2006.
- A. Oloyede and N. D. Broom, “The generalized consolidation of articular cartilage: an investigation of its near-physiological response to static load,” Connective Tissue Research, vol. 31, no. 1, pp. 75–86, 1994.
- R. Shirazi and A. Shirazi-Adl, “Analysis of articular cartilage as a composite using nonlinear membrane elements for collagen fibrils,” Medical Engineering and Physics, vol. 27, no. 10, pp. 827–835, 2005.
- P. Julkunen, W. Wilson, J. S. Jurvelin, and R. K. Korhonen, “Composition of the pericellular matrix modulates the deformation behaviour of chondrocytes in articular cartilage under static loading,” Medical and Biological Engineering and Computing, vol. 47, no. 12, pp. 1281–1290, 2009.
- R. K. Korhonen and W. Herzog, “Depth-dependent analysis of the role of collagen fibrils, fixed charges and fluid in the pericellular matrix of articular cartilage on chondrocyte mechanics,” Journal of Biomechanics, vol. 41, no. 2, pp. 480–485, 2008.
- P. S. Donzelli, R. L. Spilker, G. A. Ateshian, and V. C. Mow, “Contact analysis of biphasic transversely isotropic cartilage layers and correlations with tissue failure,” Journal of Biomechanics, vol. 32, no. 10, pp. 1037–1047, 1999.
- G. A. Holzapfel and T. C. Gasser, “A viscoelastic model for fiber-reinforced composites at finite stains: continuum basis, computational aspects and applications,” Computer Methods in Applied Mechanics and Engineering, vol. 190, no. 34, pp. 4379–4403, 2001.
- G. D. Pins, E. K. Huang, D. L. Christiansen, and F. H. Silver, “Effects of static axial strain on the tensile properties and failure mechanisms of self-assembled collagen fibers,” Journal of Applied Polymer Science, vol. 63, no. 11, pp. 1429–1440, 1997.
- J. P. A. Arokoski, M. M. Hyttinen, T. Lapveteläinen et al., “Decreased birefringence of the superficial zone collagen network in the canine knee (stifle) articular cartilage after long distance running training, detected by quantitative polarised light microscopy,” Annals of the Rheumatic Diseases, vol. 55, no. 4, pp. 253–264, 1996.
- M. J. Kääb, R. G. Richards, K. Ito, I. Ap Gwynn, and H. P. Nötzli, “Deformation of chondrocytes in articular cartilage under compressive load: a morphological study,” Cells Tissues Organs, vol. 175, no. 3, pp. 133–139, 2003.
- N. D. Broom and D. L. Marra, “Ultrastructural evidence for fibril-to-fibril associations in articular cartilage and their functional implication,” Journal of Anatomy, vol. 146, pp. 185–200, 1986.
- S. Olsen and A. Oloyede, “A finite element analysis methodology for representing the articular cartilage functional structure,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 5, no. 6, pp. 377–386, 2002.
- S. Olsen, A. Oloyede, and C. Adam, “A finite element formulation and program to study transient swelling and load-carriage in healthy and degenerate articular cartilage,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 7, no. 2, pp. 111–120, 2004.
- W. Wilson, C. C. van Donkelaar, R. van Rietbergen, and R. Huiskes, “The role of computational models in the search for the mechanical behavior and damage mechanisms of articular cartilage,” Medical Engineering and Physics, vol. 27, no. 10, pp. 810–826, 2005.
- W. M. Lai, J. S. Hou, and V. C. Mow, “A triphasic theory for the swelling and deformation behaviors of articular cartilage,” Journal of Biomechanical Engineering, vol. 113, no. 3, pp. 245–258, 1991.
- J. M. Huyghe, W. Wilson, and K. Malakpoor, “On the thermodynamical admissibility of the triphasic theory of charged hydrated tissues,” Journal of Biomechanical Engineering, vol. 131, no. 4, Article ID 044504, 2009.
- B. R. Simon, J. P. Liable, D. Pflaster, Y. Yuan, and M. H. Krag, “A poroelastic finite element formulation including transport and swelling in soft tissue structures,” Journal of Biomechanical Engineering, vol. 118, no. 1, pp. 1–9, 1996.
- M. A. Biot, “General theory of three-dimensional consolidation,” Journal of Applied Physics, vol. 12, no. 2, pp. 155–164, 1941.
- J. M. Huyghe and J. D. Janssen, “Quadriphasic mechanics of swelling incompressible porous media,” International Journal of Engineering Science, vol. 35, no. 8, pp. 793–802, 1997.
- A. J. H. Frijns, J. M. Huyghe, and J. D. Janssen, “A validation of the quadriphasic mixture theory for intervertebral disc tissue,” International Journal of Engineering Science, vol. 35, no. 15, pp. 1419–1429, 1997.
- W. Wilson, C. C. van Donkelaar, and J. M. Huyghe, “A comparison between mechano-electrochemical and biphasic swelling theories for soft hydrated tissues,” Journal of Biomechanical Engineering, vol. 127, no. 1, pp. 158–165, 2005.
- A. Maroudas and C. Bannon, “Measurement of swelling pressure in cartilage and comparison with osmotic pressure of constituent proteoglycans,” Biorheology, vol. 18, no. 3-6, pp. 619–632, 1981.
- S. R. Eisenberg and A. J. Grodzinsky, “The kinetics of chemically induced nonequilibrium swelling of articular cartilage and corneal stroma,” Journal of Biomechanical Engineering, vol. 109, no. 1, pp. 79–89, 1987.
- A. Maroudas, “Physicochemical properties of cartilage in the light of ion exchange theory,” Biophysical Journal, vol. 8, no. 5, pp. 575–595, 1968.
- A. C. Chen, W. C. Bae, R. M. Schinagl, and R. L. Sah, “Depth- and strain-dependent mechanical and electromechanical properties of full-thickness bovine articular cartilage in confined compression,” Journal of Biomechanics, vol. 34, no. 1, pp. 1–12, 2001.
- T. C. Gasser, R. W. Ogden, and G. A. Holzapfel, “Hyperelastic modelling of arterial layers with distributed collagen fibre orientations,” Journal of the Royal Society Interface, vol. 3, no. 6, pp. 15–35, 2006.
- T. M. Quinn and V. Morel, “Microstructural modeling of collagen network mechanics and interactions with the proteoglycan gel in articular cartilage,” Biomechanics and Modeling in Mechanobiology, vol. 6, no. 1-2, pp. 73–82, 2007.
- M. R. Buckley, J. P. Gleghorn, L. J. Bonassar, and I. Cohen, “Mapping the depth dependence of shear properties in articular cartilage,” Journal of Biomechanics, vol. 41, no. 11, pp. 2430–2437, 2008.
- R. M. Schinagl, D. Gurskis, A. C. Chen, and R. L. Sah, “Depth-dependent confined compression modulus of full-thickness bovine articular cartilage,” Journal of Orthopaedic Research, vol. 15, no. 4, pp. 499–506, 1997.
- H. Nötzli and J. Clark, “Deformation of loaded articular cartilage prepared for scanning electron microscopy with rapid freezing and freeze-substitution fixation,” Journal of Orthopaedic Research, vol. 15, no. 1, pp. 76–86, 1997.
- N. D. Broom and H. Silyn-Roberts, “The three-dimensional 'knit' of collagen fibrils in articular cartilage,” Connective Tissue Research, vol. 23, no. 4, pp. 261–277, 1989.
- G. A. Ateshian, “The role of interstitial fluid pressurization in articular cartilage lubrication,” Journal of Biomechanics, vol. 42, no. 9, pp. 1163–1176, 2009.
- P. Bursac, C. V. McGrath, S. R. Eisenberg, and D. Stamenovic, “A microstructural model of elastostatic properties of articular cartilage in confined compression,” Journal of Biomechanical Engineering, vol. 122, no. 4, pp. 347–353, 2000.
- M. J. Buehler, “Nanomechanics of collagen fibrils under varying cross-link densities: atomistic and continuum studies,” Journal of the Mechanical Behavior of Biomedical Materials, vol. 1, no. 1, pp. 59–67, 2008.
- M. H. Chen and N. Broom, “On the ultrastructure of softened cartilage: a possible model for structural transformation,” Journal of Anatomy, vol. 192, no. 3, pp. 329–341, 1998.
- M. H. Chen and N. D. Broom, “Concerning the ultrastructural origin of large-scale swelling in articular cartilage,” Journal of Anatomy, vol. 194, no. 3, pp. 445–461, 1999.
- N. Broom, M. H. Chen, and A. Hardy, “A degeneration-based hypothesis for interpreting fibrillar changes in the osteoarthritic cartilage matrix,” Journal of Anatomy, vol. 199, no. 6, pp. 683–698, 2001.
- N. D. Broom and H. Silyn-Roberts, “Collagen-collagen versus collagen-proteoglycan interactions in the determination of cartilage strength,” Arthritis and Rheumatism, vol. 33, no. 10, pp. 1512–1517, 1990.
- V. C. Mow, A. F. Mak, and W. M. Lai, “Viscoelastic properties of proteoglycan subunits and aggregates in varying solution concentrations,” Journal of Biomechanics, vol. 17, no. 5, pp. 325–338, 1984.
- J. Seog, D. Dean, A. H. K. Plaas, S. Wong-Palms, A. J. Grodzinsky, and C. Ortiz, “Direct measurement of glycosaminoglycan intermolecular interactions via high-resolution force spectroscopy,” Macromolecules, vol. 35, no. 14, pp. 5601–5615, 2002.
- A. Maroudas, “Biophysical chemistry of cartilaginous tissues with special reference to solute and fluid transport,” Biorheology, vol. 12, no. 3-4, pp. 233–248, 1975.
- T. M. Quinn, P. Dierickx, and A. J. Grodzinsky, “Glycosaminoglycan network geometry may contribute to anisotropic hydraulic permeability in cartilage under compression,” Journal of Biomechanics, vol. 34, no. 11, pp. 1483–1490, 2001.
- D. A. D. Parry, G. R. G. Barnes, and A. S. Craig, “A comparison of the size distribution of collagen fibrils in connective tissues as a function of age and a possible relation between fibril size distribution and mechanical properties,” Proceedings of the Royal Society of London B, vol. 203, no. 1152, pp. 305–321, 1978.
- C. J. Doillon, M. G. Dunn, and F. H. Silver, “Relationship between mechanical properties and collagen structure of closed and open wounds,” Journal of Biomechanical Engineering, vol. 110, no. 4, pp. 352–356, 1988.
- D. L. Christiansen, E. K. Huang, and F. H. Silver, “Assembly of type I collagen: fusion of fibril subunits and the influence of fibril diameter on mechanical properties,” Matrix Biology, vol. 19, no. 5, pp. 409–420, 2000.
- R. Krishnan, M. Kopacz, and G. A. Ateshian, “Experimental verification of the role of interstitial fluid pressurization in cartilage lubrication,” Journal of Orthopaedic Research, vol. 22, no. 3, pp. 565–570, 2004.
- G. A. Ateshian, H. Wang, and W. M. Lai, “The role of interstitial fluid pressurization and surface porosities on the boundary friction of articular cartilage,” Journal of Tribology, vol. 120, no. 2, pp. 241–248, 1998.
- R. Krishnan, E. N. Mariner, and G. A. Ateshian, “Effect of dynamic loading on the frictional response of bovine articular cartilage,” Journal of Biomechanics, vol. 38, no. 8, pp. 1665–1673, 2005.
- F. Guilak and V. C. Mow, “The mechanical environment of the chondrocyte: a biphasic finite element model of cell-matrix interactions in articular cartilage,” Journal of Biomechanics, vol. 33, no. 12, pp. 1663–1673, 2000.
- S. K. Han, S. Federico, A. Grillo, G. Giaquinta, and W. Herzog, “The mechanical behaviour of chondrocytes predicted with a micro-structural model of articular cartilage,” Biomechanics and Modeling in Mechanobiology, vol. 6, no. 3, pp. 139–150, 2007.
- W. Wilson, N. J. B. Driessen, C. C. van Donkelaar, and K. Ito, “Prediction of collagen orientation in articular cartilage by a collagen remodeling algorithm,” Osteoarthritis and Cartilage, vol. 14, no. 11, pp. 1196–1202, 2006.
- N. J. B. Driessen, W. Wilson, C. V. C. Bouten, and F. P. T. Baaijens, “A computational model for collagen fibre remodelling in the arterial wall,” Journal of Theoretical Biology, vol. 226, no. 1, pp. 53–64, 2004.
- T. Nagel and D. J. Kelly, “The composition of engineered cartilage at the time of implantation determines the likelihood of regenerating tissue with a normal collagen architecture,” Tissue Engineering Part A, 2012.
- H. Isaksson, I. Gröngröft, W. Wilson et al., “Remodeling of fracture callus in mice is consistent with mechanical loading and bone remodeling theory,” Journal of Orthopaedic Research, vol. 27, no. 5, pp. 664–672, 2009.
- H. Isaksson, C. C. van Donkelaar, R. Huiskes, J. Yao, and K. Ito, “Determining the most important cellular characteristics for fracture healing using design of experiments methods,” Journal of Theoretical Biology, vol. 255, no. 1, pp. 26–39, 2008.
- C. C. van Donkelaar and R. Huiskes, “The PTHrP-Ihh feedback loop in the embryonic growth plate allows PTHrP to control hypertrophy and Ihh to regulate proliferation,” Biomechanics and Modeling in Mechanobiology, vol. 6, no. 1-2, pp. 55–62, 2007.
- H. Isaksson, “Recent advances in mechanobiological modeling of bone regeneration,” Mechanics Research Communications, vol. 42, pp. 22–31, 2012.
- C. C. van Donkelaar and W. Wilson, “Mechanics of chondrocyte hypertrophy,” Biomechanics and Modeling in Mechanobiology, vol. 11, no. 5, pp. 655–664, 2012.
- P. Julkunen, T. Harjula, J. Iivarinen et al., “Biomechanical, biochemical and structural correlations in immature and mature rabbit articular cartilage,” Osteoarthritis and Cartilage, vol. 17, no. 12, pp. 1628–1638, 2009.
- M. A. Sweigart, C. F. Zhu, D. M. Burt et al., “Intraspecies and interspecies comparison of the compressive properties of the medial meniscus,” Annals of Biomedical Engineering, vol. 32, no. 11, pp. 1569–1579, 2004.
- F. Lei and A. Z. Szeri, “Inverse analysis of constitutive models: biological soft tissues,” Journal of Biomechanics, vol. 40, no. 4, pp. 936–940, 2007.
- L. G. Alexopoulos, G. M. Williams, M. L. Upton, L. A. Setton, and F. Guilak, “Osteoarthritic changes in the biphasic mechanical properties of the chondrocyte pericellular matrix in articular cartilage,” Journal of Biomechanics, vol. 38, no. 3, pp. 509–517, 2005.
- H. Yao and W. Y. Gu, “Physical signals and solute transport in human intervertebral disc during compressive stress relaxation: 3D finite element analysis,” Biorheology, vol. 43, no. 3-4, pp. 323–335, 2006.
- R. K. Korhonen, P. Julkunen, J. Rieppo, R. Lappalainen, Y. T. Konttinen, and J. S. Jurvelin, “Collagen network of articular cartilage modulates fluid flow and mechanical stresses in chondrocyte,” Biomechanics and Modeling in Mechanobiology, vol. 5, no. 2-3, pp. 150–159, 2006.
- P. Kiviranta, J. Rieppo, R. K. Korhonen, P. Julkunen, J. Töyräs, and J. S. Jurvelin, “Collagen network primarily controls poisson's ratio of bovine articular cartilage in compression,” Journal of Orthopaedic Research, vol. 24, no. 4, pp. 690–699, 2006.
- J. Yao, P. D. Funkenbusch, J. Snibbe, M. Maloney, and A. L. Lerner, “Sensitivities of medial meniscal motion and deformation to material properties of articular cartilage, meniscus and meniscal attachments using design of experiments methods,” Journal of Biomechanical Engineering, vol. 128, no. 3, pp. 399–408, 2006.
- F. H. Dar, J. R. Meakin, and R. M. Aspden, “Statistical methods in finite element analysis,” Journal of Biomechanics, vol. 35, no. 9, pp. 1155–1161, 2002.
- M. S. Phadge, Quality Engineering Using Robust Design, Prentice Hall, Englewood Cliffs, NJ, USA, 1989.
- H. Isaksson, C. C. van Donkelaar, and K. Ito, “Sensitivity of tissue differentiation and bone healing predictions to tissue properties,” Journal of Biomechanics, vol. 42, no. 5, pp. 555–564, 2009.
- J. Yao, A. D. Salo, J. Lee, and A. L. Lerner, “Sensitivity of tibio-menisco-femoral joint contact behavior to variations in knee kinematics,” Journal of Biomechanics, vol. 41, no. 2, pp. 390–398, 2008.
- F. H. Dar and R. M. Aspden, “A finite element model of an idealized diarthrodial joint to investigate the effects of variation in the mechanical properties of the tissues,” Proceedings of the Institution of Mechanical Engineers, Part H, vol. 217, no. 5, pp. 341–348, 2003.
- A. Chagnon, C. E. Aubin, and I. Villemure, “Biomechanical influence of disk properties on the load transfer of healthy and degenerated disks using a poroelastic finite element model,” Journal of Biomechanical Engineering, vol. 132, no. 11, p. 111006, 2010.
- G. Taguchi, S. Chowdhury, and S. Taguchi, Robust Engineering, McGraw-Hill, New York, NY, USA, 2000.
- E. Kim, F. Guilak, and M. A. Haider, “An axisymmetric boundary element model for determination of articular cartilage pericellular matrix properties in situ via inverse analysis of chondron deformation,” Journal of Biomechanical Engineering, vol. 132, no. 3, 2010.
- J. H. Plochocki, C. V. Ward, and D. E. Smith, “Evaluation of the chondral modeling theory using fe-simulation and numeric shape optimization,” Journal of Anatomy, vol. 214, no. 5, pp. 768–777, 2009.
- L. Cao, I. Youn, F. Guilak, and L. A. Setton, “Compressive properties of mouse articular cartilage determined in a novel micro-indentation test method and biphasic finite element model,” Journal of Biomechanical Engineering, vol. 128, no. 5, pp. 766–771, 2006.
- V. C. Mow, M. C. Gibbs, W. M. Lai, W. B. Zhu, and K. A. Athanasiou, “Biphasic indentation of articular cartilage—II. A numerical algorithm and an experimental study,” Journal of Biomechanics, vol. 22, no. 8-9, pp. 853–861, 1989.
- K. A. M. Kulmala, R. K. Korhonen, P. Julkunen et al., “Diffusion coefficients of articular cartilage for different CT and MRI contrast agents,” Medical Engineering and Physics, vol. 32, no. 8, pp. 878–882, 2010.
- M. A. Baldwin, C. W. Clary, C. K. Fitzpatrick, J. S. Deacy, L. P. Maletsky, and P. J. Rullkoetter, “Dynamic finite element knee simulation for evaluation of knee replacement mechanics,” Journal of Biomechanics, vol. 45, no. 3, pp. 474–483, 2012.
- C. K. Fitzpatrick, M. A. Baldwin, C. W. Clary, L. P. Maletsky, and P. J. Rullkoetter, “Evaluating knee replacement mechanics during ADL with PID-controlled dynamic finite element analysis,” Computer Methods in Biomechanics and Biomedical Engineering, 2012.
- P. J. Laz and M. Browne, “A review of probabilistic analysis in orthopaedic biomechanics,” Proceedings of the Institution of Mechanical Engineers, Part H, vol. 224, no. 8, pp. 927–943, 2010.
- T. M. Guess, G. Thiagarajan, M. Kia, and M. Mishra, “A subject specific multibody model of the knee with menisci,” Medical Engineering and Physics, vol. 32, no. 5, pp. 505–515, 2010.
- A. Seifzadeh, D. C. D. Oguamanam, N. Trutiak, M. Hurtig, and M. Papini, “Determination of nonlinear fibre-reinforced biphasic poroviscoelastic constitutive parameters of articular cartilage using stress relaxation indentation testing and an optimizing finite element analysis,” Computer Methods and Programs in Biomedicine, vol. 107, no. 2, pp. 315–326, 2012.
- E. J. Miller, R. F. Riemer, T. L. Haut Donahue, and K. R. Kaufman, “Experimental validation of a tibiofemoral model for analyzing joint force distribution,” Journal of Biomechanics, vol. 42, no. 9, pp. 1355–1359, 2009.
- B. M. Nigg, “General comments about modelling,” in Biomechanics of the Musculo-Skeletal System, B. M. Nigg and W. Herzog, Eds., pp. 513–522, John Wiley & Sons, Toronto, Canada, 2007.
- S. R. Eisenberg and A. J. Grodzinsky, “Swelling of articular cartilage and other connective tissues: electromechanochemical forces,” Journal of Orthopaedic Research, vol. 3, no. 2, pp. 148–159, 1985.
- R. K. Korhonen and J. S. Jurvelin, “Compressive and tensile properties of articular cartilage in axial loading are modulated differently by osmotic environment,” Medical Engineering and Physics, vol. 32, no. 2, pp. 155–160, 2010.
- R. K. Korhonen, S. K. Han, and W. Herzog, “Osmotic loading of articular cartilage modulates cell deformations along primary collagen fibril directions,” Journal of Biomechanics, vol. 43, no. 4, pp. 783–787, 2010.
- L. G. Alexopoulos, M. A. Haider, T. P. Vail, and F. Guilak, “Alterations in the mechanical properties of the human chondrocyte pericellular matrix with osteoarthritis,” Journal of Biomechanical Engineering, vol. 125, no. 3, pp. 323–333, 2003.
- A. J. Bidhendi and R. K. Korhonen, “A finite element study of micropipette aspiration of single cells: effect of compressibility,” Computational and Mathematical Methods in Medicine, vol. 2012, Article ID 192618, 9 pages, 2012.
- L. Ng, H. H. Hung, A. Sprunt, S. Chubinskaya, C. Ortiz, and A. Grodzinsky, “Nanomechanical properties of individual chondrocytes and their developing growth factor-stimulated pericellular matrix,” Journal of Biomechanics, vol. 40, no. 5, pp. 1011–1023, 2007.
- S. M. Turunen, M. J. Lammi, S. Saarakkala et al., “The effect of collagen degradation on chondrocyte volume and morphology in bovine articular cartilage following a hypotonic challenge,” Biomechanics and Modeling in Mechanobiology, 2012.
- M. T. Nieminen, J. Töyräs, M. S. Laasanen, J. Silvennoinen, H. J. Helminen, and J. S. Jurvelin, “Prediction of biomechanical properties of articular cartilage with quantitative magnetic resonance imaging,” Journal of Biomechanics, vol. 37, no. 3, pp. 321–328, 2004.
- Y. Xia, “Relaxation anisotropy in cartilage by NMR microscopy (μMRI) at 14-μm resolution,” Magnetic Resonance in Medicine, vol. 39, no. 6, pp. 941–949, 1998.
- M. J. Nissi, J. Töyräs, M. S. Laasanen et al., “Proteoglycan and collagen sensitive MRI evaluation of normal and degenerated articular cartilage,” Journal of Orthopaedic Research, vol. 22, no. 3, pp. 557–564, 2004.
- T. S. Silvast, J. S. Jurvelin, A. S. Aula, M. J. Lammi, and J. Töyräs, “Contrast agent-enhanced computed tomography of articular cartilage: association with tissue composition and properties,” Acta Radiologica, vol. 50, no. 1, pp. 78–85, 2009.
- J. T. Mäkelä, M. R. Huttu, and R. K. Korhonen, “Structure-function relationships in osteoarthritic human hip joint articular cartilage,” Osteoarthritis and Cartilage, vol. 20, no. 11, pp. 1268–1277, 2012.
- M. Kazemi, L. P. Li, P. Savard, and M. D. Buschmann, “Creep behavior of the intact and meniscectomy knee joints,” Journal of the Mechanical Behavior of Biomedical Materials, vol. 4, no. 7, pp. 1351–1358, 2011.
- S. P. Vadher, H. Nayeb-Hashemi, P. K. Canavan, and G. M. Warner, “Finite element modeling following partial meniscectomy: effect of various size of resection,” in Proceedings of the 28th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '06), pp. 2098–2101, September 2006.
- L. Kock, C. C. van Donkelaar, and K. Ito, “Tissue engineering of functional articular cartilage: the current status,” Cell and Tissue Research, vol. 347, no. 3, pp. 613–627, 2012.
- A. I. Vasara, M. T. Nieminen, J. S. Jurvelin, L. Peterson, A. Lindahl, and I. Kiviranta, “Indentation stiffness of repair tissue after autologous chondrocyte transplantation,” Clinical Orthopaedics and Related Research, no. 433, pp. 233–242, 2005.
- M. Khoshgoftar, C. C. van Donkelaar, and K. Ito, “Mechanical stimulation to stimulate formation of a physiological collagen architecture in tissue-engineered cartilage: a numerical study,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 14, no. 2, pp. 135–144, 2011.