The influence of the representation of collagen fibre organisation on the cartilage contact mechanics of the hip joint

Cartilage tissue engineering using decellularized biomatrix hydrogel containing TGF-β-loaded alginate microspheres in mechanically loaded bioreactor

Physiochemical tissue inducers and mechanical stimulation are both efficient variables in cartilage tissue fabrication and regeneration. In the presence of biomolecules, decellularized extracellular matrix (ECM) may trigger and enhance stem cell proliferation and differentiation. Here, we investigated the controlled release of transforming growth factor beta (TGF-β1) as an active mediator of mesenchymal stromal cells (MSCs) in a biocompatible scaffold and mechanical stimulation for cartilage tissue engineering. ECM-derived hydrogel with TGF-β1-loaded alginate-based microspheres (MSs) was created to promote human MSC chondrogenic development. Ex vivo explants and a complicated multiaxial loading bioreactor replicated the physiological conditions. Hydrogels with/without MSs and TGF-β1 were highly cytocompatible. MSCs in ECM-derived hydrogel containing TGF-β1/MSs showed comparable chondrogenic gene expression levels as those hydrogels with TGF-β1 added in culture media or those without TGF-β1. However, constructs with TGF-β1 directly added within the hydrogel had inferior properties under unloaded conditions. The ECM-derived hydrogel group including TGF-β1/MSs under loading circumstances formed better cartilage matrix in an ex vivo osteochondral defect than control settings. This study demonstrates that controlled local delivery of TGF-β1 using MSs and mechanical loading is essential for neocartilage formation by MSCs and that further optimization is needed to prevent MSC differentiation towards hypertrophy.

Computational modelling of articular joints with biphasic cartilage: recent advances, challenges and opportunities

Biphasic models have been widely used to simulate the time-dependent biomechanical response of soft tissues. Modelling techniques of joints with biphasic weight-bearing soft tissues have been markedly improved over the last decade, enhancing our understanding of the function, degenerative mechanism and outcomes of interventions of joints. This paper reviews the recent advances, challenges and opportunities in computational models of joints with biphasic weight-bearing soft tissues. The review begins with an introduction of the function and degeneration of joints from a biomechanical aspect. Different constitutive models of articular cartilage, in particular biphasic materials, are illustrated in the context of the study of contact mechanics in joints. Approaches, advances and major findings of biphasic models of the hip and knee are presented, followed by a discussion of the challenges awaiting to be addressed, including the convergence issue, high computational cost and inadequate validation. Finally, opportunities and clinical insights in the areas of subject-specific modeling and tissue engineering are provided and discussed.

A novel knee joint model in FEBio with inhomogeneous fibril-reinforced biphasic cartilage simulating tissue mechanical responses during gait: data from the osteoarthritis initiative

We developed a novel knee joint model in FEBio to simulate walking. Knee cartilage was modeled using a fibril-reinforced biphasic (FRB) formulation with depth-wise collagen architecture and split-lines to account for cartilage structure. Under axial compression, the knee model with FRB cartilage yielded contact pressures, similar to reported experimental data. Furthermore, gait analysis with FRB cartilage simulated spatial and temporal trends in cartilage fluid pressures, stresses, and strains, comparable to those of a fibril-reinforced poroviscoelastic (FRPVE) material in Abaqus. This knee joint model in FEBio could be used for further studies of knee disorders using physiologically relevant loading.

Computational assessment on the impact of collagen fiber orientation in cartilages on healthy and arthritic knee kinetics/kinematics

BACKGROUND

The inhomogeneous distribution of collagen fiber in cartilage can substantially influence the knee kinematics. This becomes vital for understanding the mechanical response of soft tissues, and cartilage deterioration including osteoarthritis (OA). Though the conventional computational models consider geometrical heterogeneity along with fiber reinforcements in the cartilage model as material heterogeneity, the influence of fiber orientation on knee kinetics and kinematics is not fully explored. This work examines how the collagen fiber orientation in the cartilage affects the healthy (intact knee) and arthritic knee response over multiple gait activities like running and walking.

METHODS

A 3D finite element knee joint model is used to compute the articular cartilage response during the gait cycle. A fiber-reinforced porous hyper elastic (FRPHE) material is used to model the soft tissue. A split-line pattern is used to implement the fiber orientation in femoral and tibial cartilage. Four distinct intact cartilage models and three OA models are simulated to assess the impact of the orientation of collagen fibers in a depth wise direction. The cartilage models with fibers oriented in parallel, perpendicular, and inclined to the articular surface are investigated for multiple knee kinematics and kinetics.

FINDINGS

The comparison of models with fiber orientation parallel to articulating surface for walking and running gait has the highest elastic stress and fluid pressure compared with inclined and perpendicular fiber-oriented models. Also, the maximum contact pressure is observed to be higher in the case of intact models during the walking cycle than for OA models. In contrast, the maximum contact pressure is higher during running in OA models than in intact models. Additionally, parallel-oriented models produce higher maximum stresses and fluid pressure for walking and running gait than proximal-distal-oriented models. Interestingly, during the walking cycle, the maximum contact pressure with intact models is approximately three times higher than on OA models. In contrast, the OA models exhibit higher contact pressure during the running cycle.

INTERPRETATION

Overall, the study indicates that collagen orientation is crucial for tissue responsiveness. This investigation provides insights into the development of tailored implants.

Multiscale modelling for investigating the long-term time-dependent biphasic behaviour of the articular cartilage in the natural hip joint

A better understanding of the time-dependent biomechanical behaviour of the biphasic hip articular cartilage (AC) under physiological loadings is important to understand the onset of joint pathology and guide the clinical treatment. Current computational studies for the biphasic hip AC were usually limited to short-term duration or using elaborate loading. The present study aimed to develop a multiscale computational modelling to investigate the long-term biphasic behaviour of the hip AC under physiological loadings over multiple gait cycles. Two-scale computational modelling including a musculoskeletal model and a finite element model of the natural hip was created. These two models were then combined and used to investigate the biphasic behaviour of hip AC over 80 gait cycles. The results showed that the interstitial fluid pressure in the AC supported over 89% of the loading during gait. When the contact area was located at the AC centre, the contact pressure and fluid pressure increased over time from the first cycle to the 80th cycle, while when the contact area approached the edge, these pressures decreased first dramatically and then slowly over time. The peak stresses and strains in the solid matrix of the AC remained at a low level and increased over time from the first cycle to the 80th cycle. This study demonstrated that the long-term temporal variations of the biphasic behaviour of hip AC under physiological loadings are significant. The methodology has potentially important implications in the biomechanical studies of human cartilage and supporting the development of cartilage substitution.

Displacement voxelization to resolve mesh-image mismatch: Application in deriving dense white matter fiber strains

BACKGROUND AND OBJECTIVE

It is common to combine biomechanical modeling and medical images for multimodal analyses. However, mesh-image mismatch may occur that prevents direct information exchange. To eliminate mesh-image mismatch, we develop a simple but elegant displacement voxelization technique based on image voxel corner nodes to achieve voxel-wise strain. We then apply the technique to derive dense white matter fiber strains along whole-brain tractography (∼35 k fiber tracts consisting of ∼3.3 million sampling points) resulting from head impact.

METHODS

Displacements at image voxel corner nodes are first obtained from model simulation via scattered interpolation. Each voxel is then scaled linearly to form a unit hexahedral element. This allows convenient and efficient voxel-wise strain tensor calculation and displacement interpolation at arbitrary fiber sampling points via shape functions. Fiber strains from displacement interpolation are then compared with those from the commonly used strain tensor projection using either voxel- or element-wise strain tensors.

RESULTS

Based on a synthetic displacement field, fiber strains interpolated from voxelized displacement are considerably more accurate than those from strain tensor projection relative to the prescribed ground-truth (determinant of coefficient (R²) of 1.00 and root mean squared error (RMSE) of 0.01 vs. 0.87 and 0.10, respectively). For a set of real-world reconstructed head impacts (N = 53), the strain tensor projection method performs similarly poorly (R² of 0.80-0.90 and RMSE of 0.03-0.07), with overestimation strongly correlated with strain magnitude (Pearson correlation coefficient >0.9). Up to ∼15% of the fiber strains are overestimated by more than the lower bound of a conservative injury threshold of 0.09. The percentage increases to ∼37% when halving the threshold. Voxel interpolation is also significantly more efficient (15 s vs. 40 s for element strain tensor projection, without parallelization).

CONCLUSIONS

Voxelized displacement interpolation is considerably more accurate and efficient in deriving dense white matter fiber strains than strain tensor projection. The latter generally overestimates with overestimation magnitude strongly correlating with fiber strain magnitude. Displacement voxelization is an effective technique to eliminate mesh-image mismatch and generates a convenient image representation of tissue deformation. This technique can be generalized to broadly facilitate a diverse range of image-related biomechanical problems for multimodal analyses. The convenient image format may also promote and facilitate biomechanical data sharing in the future.

The effect of highly inhomogeneous biphasic properties on mechanical behaviour of articular cartilage

BACKGROUND AND OBJECTIVE

Investigating the biomechanics of cartilage could help to understand the unique load-bearing property of the cartilage and optimize the scaffold design in tissue-engineering. It is important to model the cartilage as a highly inhomogeneous fibril-reinforced biphasic material to represent its complex composition and structure. The depth-dependent and strain-dependent properties of the cartilage would also play an important role in its mechanical behaviour. However, the differences in representing the cartilage as a highly inhomogeneous model or as simplified models still remain unclear. Hence, in this study, a highly inhomogeneous fibril-reinforced biphasic cartilage model considering both the depth-dependent and strain-dependent properties was constructed; the effect of highly inhomogeneous properties on the mechanical behaviour of articular cartilage was investigated.

METHODS

A finite element model of the cartilage was developed based on a flat-ended indentation test. Compressive forces were applied to four various inhomogeneous layered models through a porous indenter (Model 1: nine layers with strain-dependent permeability; Model 2: three layers with strain-dependent permeability; Model 3: single layer with strain-dependent permeability; Model 4: nine layers with constant permeability).

RESULTS

Models 1 and 2 provided similar results with less than 3% difference in the peak effective stress, contact pressure, fluid pressure as well as fluid support ratio. However, Model 1 to Model 3 differed in stress and strain distribution patterns along depth over prolonged loads, which may provide an important insight into the highly inhomogeneous depth-dependent properties of cartilage. In addition, Model 1 with strain-dependent permeability demonstrated an enhanced capability on fluid pressurisation as compared with Model 4 which had constant permeability.

CONCLUSIONS

A highly inhomogeneous fibril-reinforced biphasic model considering both depth-dependent and strain-dependent properties was developed in this study, in order to illustrate the effect of highly inhomogeneous properties on the mechanical behaviour of the articular cartilage. The number of layers in the models with depth-dependent properties should be selected according to the research questions and clinical demands. The model with strain-dependent permeability offers an enhanced capability on fluid pressurisation. In future studies, the proposed model could be adopted in cell-models to provide more in-depth information or in tissue-engineering to optimize the depth-dependent scaffold structure.

Development and validation of a finite-element musculoskeletal model incorporating a deformable contact model of the hip joint during gait

Musculoskeletal models provide non-invasive and subject-specific biomechanical investigations of the musculoskeletal system. In a musculoskeletal model, muscle forces contribute to the deformation and kinematics of the joint which in turn would alter moment arms of muscles and ground reaction forces and thus affect the prediction of muscle forces and contact forces and contact mechanics of the joint. By far, deformable contact models of the hip have not been considered in musculoskeletal models, and the role of kinematics and deformation within the hip in muscle forces and hip contact mechanics is unknown. In this study, an FE musculoskeletal model including bones, joints and muscles of the lower extremity was developed. A deformable contact model of the hip joint was incorporated and coupled into the musculoskeletal model. Joint angles and ground reaction forces during gait were used as inputs. Optimization minimizing the sum of muscle stresses squared was performed directly to the FE musculoskeletal model in order to simultaneously solve muscle forces and contact forces and contact stresses of the hip joint within a single framework. The calculated hip contact forces corresponded well to the in vivo measurement data. The maximum hip contact stress was 6.5 MPa and occurred at weight-acceptance. The influence of kinematics and deformation in the hip on muscles forces and hip contact forces was minimal and not sensitive to variations in the thickness and properties of the joint cartilage during gait. This suggests that the uncoupled approach in which the hip contact forces and contact mechanics are simulated in separate frameworks would serve as an effective and efficient alternative for subject-specific modelling of the hip. This study provides guidance for the level of complexity needed for future hip models and can be used to evaluate biomechanical changes of the musculoskeletal system following interventions.

Simulation of surface strain in tibiofemoral cartilage during walking for the prediction of collagen fiber orientation

The collagen fibers in the superficial layer of tibiofemoral articular cartilage exhibit distinct patterns in orientation revealed by split lines. In this study, we introduce a simulation framework to predict cartilage surface loading during walking to investigate if split line orientations correspond with principal strain directions in the cartilage surface. The two-step framework uses a multibody musculoskeletal model to predict tibiofemoral kinematics which are then imposed on a deformable surface model to predict surface strains. The deformable surface model uses absolute nodal coordinate formulation (ANCF) shell elements to represent the articular surface and a system of spring-dampers and internal pressure to represent the underlying cartilage. Simulations were performed to predict surface strains due to osmotic pressure, loading induced by walking, and the combination of both loading due to pressure and walking. Time-averaged magnitude-weighted first principal strain directions agreed well with split line maps from the literature for both the osmotic pressure and combined cases. This result suggests there is indeed a connection between collagen fiber orientation and mechanical loading, and indicates the importance of accounting for the pre-strain in the cartilage surface due to osmotic pressure.

Hip chondrolabral mechanics during activities of daily living: Role of the labrum and interstitial fluid pressurization

Osteoarthritis of the hip can result from mechanical factors, which can be studied using finite element (FE) analysis. FE studies of the hip often assume there is no significant loss of fluid pressurization in the articular cartilage during simulated activities and approximate the material as incompressible and elastic. This study examined the conditions under which interstitial fluid load support remains sustained during physiological motions, as well as the role of the labrum in maintaining fluid load support and the effect of its presence on the solid phase of the surrounding cartilage. We found that dynamic motions of gait and squatting maintained consistent fluid load support between cycles, while static single-leg stance experienced slight fluid depressurization with significant reduction of solid phase stress and strain. Presence of the labrum did not significantly influence fluid load support within the articular cartilage, but prevented deformation at the cartilage edge, leading to lower stress and strain conditions in the cartilage. A morphologically accurate representation of collagen fibril orientation through the thickness of the articular cartilage was not necessary to predict fluid load support. However, comparison with simplified fibril reinforcement underscored the physiological importance. The results of this study demonstrate that an elastic incompressible material approximation is reasonable for modeling a limited number of cyclic motions of gait and squatting without significant loss of accuracy, but is not appropriate for static motions or numerous repeated motions. Additionally, effects seen from removal of the labrum motivate evaluation of labral reattachment strategies in the context of labral repair.

Ulna-humerus contact mechanics: Finite element analysis and experimental measurements using a tactile pressure sensor

Elbow articular cartilage withstands high compressive and shear forces while protecting the bone from excessive loading. Better understanding of elbow cartilage contact mechanics can provide insight into cartilage degeneration. In this study a tactile pressure sensor was used to measure the contact pressure distribution within the ulno-humeral joint of two cadaver specimens at 20° flexion angle across three different axial loads of 80 N, 110 N, and 140 N. Corresponding 3D finite element (FE) models were constructed from magnetic resonance imaging (MRI) and contact analysis was performed for each specimen with boundary and loading conditions identical to the experiment. Direct comparison between FE results and experimental measurements was conducted for the validation of the FE models and a sensitivity analysis was employed for assessing the effect of cartilage parameters on the model's outputs. The results showed a good agreement between the FE models and the experiments in terms of contact characteristics. The sensitivity analysis demonstrated that outcomes of the model, particularly peak contact pressure is more sensitive to the Poisson's ratio rather than to Young's modulus under static conditions. This result suggests that selection of Poisson's ratio is very critical for accurate prediction of contact mechanics within the ulno-humeral joint.