The Influence of Fiber Orientation on the Equilibrium Properties of Neutral and Charged Biphasic Tissues

Constitutive models facilitate investigation into load bearing mechanisms of biological tissues and may aid attempts to engineer tissue replacements. In soft tissue models, a commonly made assumption is that collagen fibers can only bear tensile loads. Previous computational studies have demonstrated that radially aligned fibers stiffen a material in unconfined compression most by limiting lateral expansion while vertically aligned fibers buckle under the compressive loads. In this short communication, we show that in conjunction with swelling, these intuitive statements can be violated at small strains. Under such conditions, a tissue with fibers aligned parallel to the direction of load initially provides the greatest resistance to compression. The results are further put into the context of a Benninghoff architecture for articular cartilage. The predictions of this computational study demonstrate the effects of varying fiber orientations and an initial tare strain on the apparent material parameters obtained from unconfined compression tests of charged tissues.

A Penetration-Based Finite Element Method for Hyperelastic 3D Biphasic Tissues in Contact. Part II: Finite Element Simulations

The penetration method allows for the efficient finite element simulation of contact between soft hydrated biphasic tissues in diarthrodial joints. Efficiency of the method is achieved by separating the intrinsically nonlinear contact problem into a pair of linked biphasic finite element analyses, in which an approximate, spatially and temporally varying contact traction is applied to each of the contacting tissues. In Part I of this study, we extended the penetration method to contact involving nonlinear biphasic tissue layers, and demonstrated how to derive the approximate contact traction boundary conditions. The traction derivation involves time and space dependent natural boundary conditions, and requires special numerical treatment. This paper (Part II) describes how we obtain an efficient nonlinear finite element procedure to solve for the biphasic response of the individual contacting layers. In particular, alternate linearization of the nonlinear weak form, as well as both velocity-pressure, v‐p, and displacement-pressure, u‐p, mixed formulations are considered. We conclude that the u‐p approach, with linearization of both the material law and the deformation gradients, performs best for the problem at hand. The nonlinear biphasic contact solution will be demonstrated for the motion of the glenohumeral joint of the human shoulder joint.

A Penetration-Based Finite Element Method for Hyperelastic 3D Biphasic Tissues in Contact: Part 1-Derivation of Contact Boundary Conditions

In this study, we extend the penetration method, previously introduced to simulate contact of linear hydrated tissues in an efficient manner with the finite element method, to problems of nonlinear biphasic tissues in contact. This paper presents the derivation of contact boundary conditions for a biphasic tissue with hyperelastic solid phase using experimental kinematics data. Validation of the method for calculating these boundary conditions is demonstrated using a canonical biphasic contact problem. The method is then demonstrated on a shoulder joint model with contacting humerus and glenoid tissues. In both the canonical and shoulder examples, the resulting boundary conditions are found to satisfy the kinetic continuity requirements of biphasic contact. These boundary conditions represent input to a three-dimensional nonlinear biphasic finite element analysis; details of that finite element analysis will be presented in a manuscript to follow.