Search PubMed⌕ Search

Biomedical subjects

V C Mow

Publications and source records attributed to V C Mow.

At least 91 records · Page 5Linked to original sources

Correlating magnetic resonance imaging with the biochemical content of the normal human intervertebral disc.

Magnetic resonance imaging was used to determine the T2 relaxation times of prepared proteoglycan solutions and of normal human intervertebral disc tissue from the annulus fibrosus (AF) and nucleus pulposus (NP). The collagen, proteoglycan, and water contents of the disc tissue samples were determined by biochemical assays after they were scanned. Correlations among 1/T2, collagen, proteoglycan, and water contents of the tissue samples and among 1/T2, water, and proteoglycan contents of the proteoglycan solutions were calculated. A moderate negative correlation between 1/T2 and water content was noted for the tissue samples, and a very high negative correlation was found between 1/T2 and water content for the proteoglycan solutions. The very high positive correlation between 1/T2 and proteoglycan content of the proteoglycan solutions is probably due to this negative correlation between 1/T2 and water content. There was no significant correlation between 1/T2 and proteoglycan content of the tissues. The moderate positive correlation between 1/T2 and collagen content is probably due to the high negative correlation between collagen content and water content. No significant correlation was found between the collagen and proteoglycan contents of the tissues. Thus it appears that the data confirm previous reports in the literature that the collagen of the disc tissue functions to control its water content.

Adult↗

An analysis of the squeeze-film lubrication mechanism for articular cartilage.

An asymptotic analysis of a lubrication problem is presented for a model of articular cartilage and synovial fluid under the squeeze-film condition. This model is based upon the following constitutive assumptions: (1) articular cartilage is a linear porous-permeable biphasic material filled with a linearly viscous fluid (i.e. Newtonian fluid); (2) synovial fluid is also a linearly viscous fluid. The geometry of the problem is defined by assuming that (1) cartilage is a uniform layer of thickness H; (2) synovial fluid is a very thin layer compared to H; (3) the radius R of the load-supporting area (or the effective radius of curvature of joint surface, Ri) is large compared to H. Squeeze-film action is generated in the lubricant by a step loading function applied onto the two bearing surfaces. The model assumptions and the material properties yield two small parameters in the mathematical formulation. Based on these two small parameters, two coupled nonlinear partial differential equations were derived from an asymptotic analysis of the problem: one for the lubricant (analogous to the Reynolds equation) and one for the cartilage. For known properties of normal cartilage, our calculations show: (1) the cartilage layer deforms to enlarge the load-supporting area; (2) cartilage deformation acts to reduce the lateral fluid speed in the lubricant, thus prolonging the squeeze-film time which ranges from 1 to 10 s; (3) lubricant fluid in the gap is forced from the central high-pressure region into cartilage, and expelled from the tissue at the low-pressure periphery of the load-bearing region; and (4) tensile hoop stress exists at the cartilage surface despite the compressive squeeze-film loading condition. This hoop stress results directly from the radial flow of the interstitial fluid in the cartilage layer.

Biomechanical Phenomena↗

A transversely isotropic biphasic finite element model of the meniscus.

A finite element model of the meniscus is presented, based on an axisymmetric geometric approximation of the menisci and a biphasic description of the tissue as a mixture of solid and fluid components. The highly fibrous nature of the meniscal tissue is accounted for by using a fiber-reinforced, transversely isotropic description of the solid phase. This model is used to study the response of a meniscus resting on a perfectly lubricated tibial surface and subjected to distributed loads applied to the femoral surface, and to examine the effects of changes in loading conditions at the femoral and tibial interfaces. Quantities of interest include the stress, pressure and strain distributions at discrete times early in the meniscal response, and the flow of the fluid phase relative to the solid phase. Of particular interest are regions of large tensile strains which could lead to meniscal failures such as the bucket-handle tear. We show that all components of strain are positive in regions of the outer third of the meniscus, and that the maximum tensile strain perpendicular to the circumferentially arranged fibers (largest principal strain in the axisymmetric cross section) is positive throughout most of the cross section. Changing the partition of the load on the femoral surface and the permeability at the tibial surface changes the time-dependent response, but has little effect on the strain distributions at times of the order of 5 s considered in this study. The inclusion of a transversely isotropic, fibrous representation of the solid phases is shown to be essential to proper meniscal simulation. The results demonstrate the importance of the biphasic representation since the fluid phase is shown to carry a significant part of the applied load.

Biomechanical Phenomena↗

Curvature characteristics and congruence of the thumb carpometacarpal joint: differences between female and male joints.

Three-dimensional geometric models of the articular surfaces of the thumb carpometacarpal (CMC) joint were constructed using precise data obtained from stereophotogrammetry (SPG). It was demonstrated that by using a least-squares surface-fitting technique, the SPG data on the surface can accurately be described by a single parametric biquintic spline function. From this mathematical description, curvature maps of the surfaces were calculated for 13 CMC joints (eight females, average 64 yr old, five males, average 70 yr old). The surface geometry of each joint was analyzed, comparisons were made between trapezial and metacarpal surfaces of the joint and differences determined between males and females. With regard to joint surface areas, the female trapezium is significantly smaller than that of the metacarpal. The shape of the female trapezial surface is also fundamentally different from that of males. No gender-related difference exists regarding the shape of the metacarpal surface. Congruence of the two opposing articular surfaces was defined by their relative principal curvatures. From these definitions, congruence in the radioulnar and dorsovolar anatomic directions, as well as the global congruence of the joint, were calculated. Most CMC joints were found to be more congruent along the radioulnar direction than the dorsovolar direction and, globally, female joints were found to be less congruent than male joints. The concept of joint congruence has played a central role in a number of hypotheses relating to the etiology of CMC joint osteoarthritis (OA), although conflicting hypotheses do exist. The precise quantitative findings of this study may lead to an improved understanding of CMC joint OA, and perhaps explain its prevalence in the female population over 55.

Aged↗

Cartilage and diarthrodial joints as paradigms for hierarchical materials and structures.

The anatomic forms of diarthrodial joints are important structural features which provide and limit the motions required for the joint. Typically, the length scale of topographic variation of anatomic forms ranges from 0.5 to 15 cm. Articular cartilage is the thin layer of hydrated soft tissue (0.5-5.0 mm thick) covering the articulating bony ends in diarthrodial joints. This tissue has a set of unique mechanical and physicochemical properties which are responsible for its load-carrying capabilities and near-frictionless qualities. The mechanical properties of articular cartilage are determined at the tissue-scale level and these properties depend on the composition of the tissue, mainly collagen and proteoglycan, and their molecular and ultrastructural organization (ultra-scale: 10(-8)-10(-6) m). Because proteoglycans possess a high density of fixed negative charges, articular cartilage exhibits a significant Donnan osmotic pressure effect. This physicochemically derived osmotic pressure is an important component of the total swelling pressure; the other component of the total swelling pressure stems from the charge-to-charge repulsive force exerted by the closely spaced (1-1.5 nm) negative charge groups along the proteoglycan molecules. Thus these interactions take place at a nano-scale level: 10(-10)-10(-9) m. Finally, cartilage biochemistry and organization are maintained by the chondrocytes which exist at a micro-scale level (10(-7)-10(-6) m). Significant mechanoelectrochemical transduction occurs within the extracellular matrix at the micro-scale level which affects and modulates cellular anabolic and catabolic activities. At present, the exact details of these transduction mechanisms are unknown. In this review, we present a summary of the hierarchical features for articular cartilage and diarthrodial joints and tables of known material properties for cartilage. Also we summarize how the multi-scale interactions in articular cartilage provide for its unique material properties and tribological characteristics.

Animals↗

A finite element analysis of the indentation stress-relaxation response of linear biphasic articular cartilage.

The indentation problem of a thin layer of hydrated soft tissue such as cartilage or meniscus by a circular plane-ended indenter is investigated. The tissue is represented by a biphasic continuum model consisting of a solid phase (collagen and proteoglycan) and a fluid phase (interstitial water). A finite element formulation of the linear biphasic continuum equations is used to solve an axisymmetric approximation of the indentation problem. We consider stress-relaxation problems for which analytic solution is intractable; where the indenter is impermeable (solid) and/or when the interface between the indenter and tissue is perfectly adhesive. Thicknesses corresponding to a thin and thick specimen are considered to examine the effects of tissue thickness. The different flow, pressure, stress and strain fields which are predicted within the tissue, over time periods typically used in the mechanical testing of soft tissues, will be presented. Results are compared with the case of a porous free-draining indenter with a perfectly lubricated tissue-indenter interface, for which an analytic solution is available, to show the effects of friction at the tissue-indenter interface, and the effects of an impermeable indenter. While these effects are present for both thin and thick tissues, they are shown to be more significant for the thin tissue. We also examine the effects of the stiffness of the subchondral bone on the response of the soft tissue and demonstrate that the subchondral bone substrate can be modeled as a rigid, impermeable boundary. The effects of a curved tissue-subchondral bone interface, and the early time response are also studied. For physiologically reasonable levels of curvature, we will show that the curved tissue-subchrondal bone interface has negligible influence on the tissue response away from the interface. In addition, the short-time stress-relaxation responses of the tissue (e.g., at times less than 1s) demonstrate the essential role of the fluid phase in supporting the load applied to the tissue, and by extrapolation to shorter times characteristics of normal joint motion, suggest the essential role of a biphasic model in representing soft tissue behavior in joint response.

Biomechanical Phenomena↗

Articular geometry of the glenohumeral joint.

Little quantitative data exists defining the true shape of the humeral head and glenoid articular surfaces. This study uses a precise stereophotogrammetry (SPG) technique and provides highly accurate quantitative results for determining the three-dimensional geometry of glenohumeral joint articular surfaces, including their "sphericity", surface areas, cartilage thickness, and the difference in these quantities between the genders. Results indicate that glenohumeral joint surfaces may be approximated by a section of a sphere with small deviations from sphericity of less than 1% of the radius. Furthermore, results indicate that mating humeral head and glenoid articular surfaces are quite congruent and have radii within 2 mm in 88% of cases, and within 3 mm in all cases measured. The lack of anatomic stability of this joint is therefore not attributable to the relative shallowness or lack of congruence of the joint but rather to the small surface area of the glenoid which does not enclose the humeral head. Cartilage thickness results may partially explain perceptions of glenohumeral incongruity obtained from roentgenographic measurements where the glenoid appears flatter than the corresponding humeral head. Only when the actual articular cartilage surfaces are analyzed is it determined that the actual articulating surfaces do conform.

Aged↗

Interspecies comparisons of in situ intrinsic mechanical properties of distal femoral cartilage.

We measured the in situ biomechanical properties of knee joint cartilage from five species (bovine, canine, human, monkey, and rabbit) to examine the biomechanical relevance of animal models of human knee joint injuries and osteoarthritis. In situ biphasic creep indentation experiments were performed to simultaneously determine all three intrinsic material coefficients (aggregate modulus, Poisson's ratio, and permeability) of the cartilage as represented by the linear KLM biphasic model. In addition, we also assessed the effects of load bearing on these intrinsic properties at "high" and "low" weight-bearing regions on the distal femur. Our results indicate that significant differences exist in some of these material properties among species and sites. The aggregate modulus of the anterior patellar groove within each species is the lowest among all sites tested, and the permeability of the patellar groove cartilage is the highest and does not vary among species. Similarly, the Poison's ratio in the patellar groove is the lowest in all species, except in the rabbit. These results lead to the conclusion that patellar groove cartilage can undergo greater and faster compression. Thus, under high compressive loads, the cartilage of the patellar groove surface can more rapidly compress to create a congruent patellofemoral joint articulation. For any given location, no differences were found in the aggregate modulus among all the species, and no correlation was found between aggregate modulus and thickness at the test site. Thus, in the process of selecting a suitable experimental animal model of human articular cartilage, it is essential to consider the significant interspecies differences of the mechanical properties.

Analysis of Variance↗

The density and strength of proteoglycan-proteoglycan interaction sites in concentrated solutions.

Rheological flow properties of link-stable and link-free proteoglycan (PG) aggregates in concentrated solutions were measured using a cone-on-plate viscometer. A second-order constitutive model, based upon the statistical-network theories of Lodge, [Rheol. Acta 7, 379-392 (1968)] and De Kee and Carreau [J. Non-Newtonian Fluid Mech. 6, 127-143 (1979)], was developed to describe the measured steady and transient flow responses exhibited by the PG solutions. Our measurements confirmed previous experimental findings that the complex shear modulus of PG solutions depends on the frequency of the imposed small-amplitude oscillatory shear, and the apparent viscosity and primary normal-stress difference depend nonlinearly on the shear rate under steady-shear flow conditions [Mow et al., J. Biomechanics 17, 325-338 (1984b); Hardingham et al., J. orthop. Res. 5, 36-46 (1987)]. In the present study, we found that PG solutions exhibit pronounced stress overshoot responses and large hysteresis loop effects. These transient responses were shown to be sensitive to acceleration strain (i.e. the second rate of strain) as well as PG structure (i.e. link-protein stabilization). The model parameters were determined by curvefitting of the second-order constitutive model and experimental data from steady, oscillatory and transient shear flow measurements. Using this network model, we calculated the density of the idealized interaction sites existing in the PG network, and the average strength of these interaction sites. The results indicate that link-protein stabilization of PG aggregates does not change the density of interaction sites formed in the PG network, rather, it increases the average strength of these interaction sites.

Acceleration↗

Quantitation of articular surface topography and cartilage thickness in knee joints using stereophotogrammetry.

An analytical stereophotogrammetry (SPG) technique has been developed based upon some of the pioneering work of Selvik [Ph.D. thesis, University of Lund, Sweden (1974)] and Huiskes and coworkers [J. Biomechanics 18, 559-570 (1985)], and represents a fundamental step in the construction of biomechanical models of diarthrodial joints. Using this technique, the precise three-dimensional topography of the cartilage surfaces of various diarthrodial joints has been obtained. The system presented in this paper delivers an accuracy of 90 microns in the least favorable conditions with 95% coverage using the same calibration method as Huiskes et al. (1985). In addition, a method has been developed, using SPG, to quantitatively map the cartilage thickness over the entire articular surface of a joint with a precision of 134 microns (95% coverage). In the present study, our SPG system has been used to quantify the topography, including surface area, of the articular surfaces of the patella, distal femur, tibial plateau, and menisci of the human knee. Furthermore, examples of cartilage thickness maps and corresponding thickness data including coefficient of variation, minimum, maximum, and mean cartilage thickness are also provided for the cartilage surfaces of the knee. These maps illustrate significant variations over the joint surfaces which are important in the determination of the stresses and strains within the cartilage during diarthrodial joint function. In addition, these cartilage surface topographies and thickness data are essential for the development of anatomically accurate finite element models of diarthrodial joints.(ABSTRACT TRUNCATED AT 250 WORDS)

Adult↗

Experimental studies on repair of large osteochondral defects at a high weight bearing area of the knee joint: a tissue engineering study.

There is a vast clinical need for the development of an animal model to study the fundamentals of healing of injured or diseased diarthrodial joints (knee, hip, shoulder, wrist, etc). Current prosthetic replacements do not offer acceptable treatment for injuries and diseases of these joints in young active individuals. New clinical treatment modalities, based on sound biologic principles, are sought for the development of repair or healing tissues engineered to have similar biomechanical properties as normal articular cartilage. In this paper we present a brief review of this need, and propose a grafting procedure which may lead to a successful animal model for studies of long term repair of major osteochondral defects. This grafting procedure uses an autologous periosteum-bone graft or an autologous-synthetic bone replacement graft. We have applied these grafts for in vivo repair of large surgically created defects in the high weight bearing area of the distal femoral condyle of mature New Zealand white rabbits. Further, an interdisciplinary study, including histochemistry, biochemistry (composition and metabolic activities), and biomechanics (biphasic properties), was performed to assess the feasibility of our animal model to generate viable repair tissues. We found our grafting procedure produced, 8 weeks postoperatively, tissues which were very similar to those found in normal articular cartilage. However, our histological studies indicate incomplete bonding between the repair tissue and the adjacent cartilage, and lack of an appropriate superficial zone at the articular surface. These deficiencies may cause long term failure of the repair tissue. Further studies must be undertaken to enhance development of a strong bond and a collagen-rich surface zone. This may require the use of growth factors (e.g., transforming growth factors beta) capable of simulating extra collagen production, or the use of serum derived tissue glue for bonding. At present, we are pursuing these studies.

Animals↗

A triphasic theory for the swelling and deformation behaviors of articular cartilage.

Swelling of articular cartilage depends on its fixed charge density and distribution, the stiffness of its collagen-proteoglycan matrix, and the ion concentrations in the interstitium. A theory for a tertiary mixture has been developed, including the two fluid-solid phases (biphasic), and an ion phase, representing cation and anion of a single salt, to describe the deformation and stress fields for cartilage under chemical and/or mechanical loads. This triphasic theory combines the physico-chemical theory for ionic and polyionic (proteoglycan) solutions with the biphasic theory for cartilage. The present model assumes the fixed charge groups to remain unchanged, and that the counter-ions are the cations of a single-salt of the bathing solution. The momentum equation for the neutral salt and for the intersitial water are expressed in terms of their chemical potentials whose gradients are the driving forces for their movements. These chemical potentials depend on fluid pressure p, salt concentration c, solid matrix dilatation e and fixed charge density cF. For a uni-uni valent salt such as NaCl, they are given by mu i = mu io + (RT/Mi)ln[gamma 2 +/- c(c + cF)] and mu w = mu wo + [p-RT phi (2c + cF) + Bwe]/pwT, where R, T, Mi, gamma +/-, phi, pwT and Bw are universal gas constant, absolute temperature, molecular weight, mean activity coefficient of salt, osmotic coefficient, true density of water, and a coupling material coefficient, respectively. For infinitesimal strains and material isotropy, the stress-strain relationship for the total mixture stress is sigma = - pI-TcI + lambda s(trE)I + 2 musE, where E is the strain tensor and (lambda s, mu s) are the Lamé constants of the elastic solid matrix. The chemical-expansion stress (-Tc) derives from the charge-to-charge repulsive forces within the solid matrix. This theory can be applied to both equilibrium and non-equilibrium problems. For equilibrium free swelling problems, the theory yields the well known Donnan equilibrium ion distribution and osmotic pressure equations, along with an analytical expression for the "pre-stress" in the solid matrix. For the confined-compression swelling problem, it predicts that the applied compressive stress is shared by three load support mechanisms: 1) the Donnan osmotic pressure; 2) the chemical-expansion stress; and 3) the solid matrix elastic stress. Numerical calculations have been made, based on a set of equilibrium free-swelling and confined-compression data, to assess the relative contribution of each mechanism to load support. Our results show that all three mechanisms are important in determining the overall compressive stiffness of cartilage.

Biomechanical Phenomena↗

[Relationship of acromial architecture and diseases of the rotator cuff].

Variations in the architecture of the coraco-acromial arch can lead to a clinically symptomatic rotator cuff lesion. Differences in the shape and slope of the acromion, anterior acromial spurs and inferior protruding acromio-clavicular osteophytes decrease the volume of the subacromial space, leading to impingement. Recent anatomic, radiographic, biomechanic and stereophotogrammetric studies have confirmed these architectural variations and their effect on the contents of the subacromial space. Abnormal contact between the acromion and these soft tissues can result in pathological lesions. Surgical procedures should be aimed at increasing the space beneath the coraco-acromial arch to reduce wear on the rotator cuff.

Acromioclavicular Joint↗

The relationship of acromial architecture to rotator cuff disease.

Variations in the architecture of the coracoacromial arch can lead to a clinically symptomatic rotator cuff lesion. Differences in the development and morphology of the acromion, and the presence of anterior acromial spurs and inferior acromioclavicular osteophytes decrease the volume of the subacromial space, leading to impingement. Recent anatomic, radiographic, biomechanical, and SPG studies have confirmed these architectural variations and their effects on the contents of the subacromial space. Abnormal contact between the acromion and these soft tissues can lead to pathologic lesions. Surgical procedures should be directed at increasing the space beneath the coracoacromial arch to reduce wear on the rotator cuff.

Acromioclavicular Joint↗

Effects of proteoglycan extraction on the tensile behavior of articular cartilage.

We undertook an interdisciplinary biomechanical and biochemical study to explore the extent and manner in which the total pool of proteoglycans influences the kinetic and static behavior of bovine articular cartilage in tension. Two biomechanical tests were used: (a) the viscoelastic creep test and (b) a slow constant-rate uniaxial tension test; and two enzymatic proteoglycan extraction procedures were used: (a) chondroitinase ABC treatment and (b) a sequential enzymatic treatment with chondroitinase ABC, trypsin, and Streptomyces hyaluronidase. We found that the viscoelastic creep response of all cartilage specimens may be divided into two distinct phases: an initial phase (less than 15 s), characterized by a rapid increase in strain following load application, and a late phase (15 s less than or equal to t less than 25,000 s), characterized by a more gradual increase in strain. A major finding of this study is that the kinetics of the creep response is greatly influenced by the glycosaminoglycan content of the tissue. For untreated and control specimens, the initial response comprises about 50% of the total strain, while for chondroitinase ABC and sequentially extracted specimens, the initial response comprises up to 83% of the total strain. Furthermore, most untreated and control specimens did not reach equilibrium within the 25,000 s test period, while enzymatically digested specimens often reached equilibrium in less than 100 s. Thus, we conclude that through their physical restraints on collagen, the bulk of proteoglycan present in the tissue acts to retard fibrillar reorganization and alignment under tensile loading, thereby effectively preventing sudden extension of the collagen network. In contrast, the results of our slow constant-rate uniaxial tension experiment show that essentially complete extraction of proteoglycan glycosaminoglycans does not affect the intrinsic tensile stiffness and strength of cartilage specimens or the collagen network in a significant manner. Hence, an important function of the bulk proteoglycans (i.e., the large aggregating type) in cartilage is to retard the rate of stretch and alignment when a tensile load is suddenly applied. This mechanism may be useful in protecting the cartilage collagen network during physiological situations, where sudden impact forces are imposed on a joint.

Animals↗

The nonlinear characteristics of soft gels and hydrated connective tissues in ultrafiltration.

A one-dimensional ultrafiltration problem of fluid flow through a soft permeable tissue or gel under high pressure and compressive strain is solved. A finite deformation biphasic theory is used to model the behavior of the soft porous permeable solid matrix. This theory includes a Helmholtz free energy function which depends on the three principal invariants (I, II, III) of the right Cauchy-Green tensor and which satisfies the Baker-Ericksen inequalities on the principal stresses and strains. The dependence of the porosity phi f and the solidity phi s on deformation is deduced and a generalization of the exponential strain-dependent functional form for the permeability, k = k0 exp (M epsilon), of Lai and Mow (Biorheology 103, 111-123, 1980) is proposed. In this one-dimensional problem, we show that the dependence of the permeability on phi f, phi s, and III is equivalent to its dependence on hydration as proposed by Fatt and Goldstick (J. Colloid Sci. 20, 962-988, 1965). The exact solution of the ultrafiltration problem is derived and asymptotic and numerical methods are used to evaluate it. For high pressures and finite strains, the solution provides some surprising effects. The theory predicts that a material starting with a homogeneous porosity will have a strongly non-homogeneous porosity throughout the column during ultrafiltration. The resulting change in pore size through the filtration column may be very important in understanding its filtration characteristics. It is also found that there is a long delay time, up to 10 to 15 min, before the filtration velocity reaches an equilibrium. In filtration experiments where the rate of mass transport across the tissue or column of gel is important, sufficient time must be allowed for the steady state to be reached.

Animals↗

A finite deformation theory for cartilage and other soft hydrated connective tissues--I. Equilibrium results.

The determination of valid stress-strain relations for articular cartilage under finite deformation conditions is a prerequisite for constructing models for synovial joint lubrication. Under physiological conditions of high strain rates and/or high stresses in the joint, large strains occur in cartilage. A finite deformation theory valid for describing cartilage, as well as other soft hydrated connective tissues under large loads, has been developed. This theory is based on the choice of a specific Helmholtz energy function which satisfies the generalized Coleman-Noll (GCN0) condition and the Baker-Ericksen (B-E) inequalities established in finite elasticity theory. In addition, the finite deformation biphasic theory includes the effects of strain-dependent porosity and permeability. These nonlinear effects are essential for properly describing the biomechanical behavior of articular cartilage, even when strain rates are low and strains are infinitesimal. The finite deformation theory describes the large strain behavior of cartilage observed in one-dimensional confined compression experiments at equilibrium, and it reduces to the linear biphasic theory under infinitesimal strain and slow strain rate conditions. Using this theory, we have determined the material coefficients of both human and bovine articular cartilages under large strain conditions at equilibrium. The theory compares very well with experimental results.

Animals↗

Effects of friction on the unconfined compressive response of articular cartilage: a finite element analysis.

A finite element analysis is used to study a previously unresolved issue of the effects of platen-specimen friction on the response of the unconfined compression test; effects of platen permeability are also determined. The finite element formulation is based on the linear KLM biphasic model for articular cartilage and other hydrated soft tissues. A Galerkin weighted residual method is applied to both the solid phase and the fluid phase, and the continuity equation for the intrinsically incompressible binary mixture is introduced via a penalty method. The solid phase displacements and fluid phase velocities are interpolated for each element in terms of unknown nodal values, producing a system of first order differential equations which are solved using a standard numerical finite difference technique. An axisymmetric element of quadrilateral cross-section is developed and applied to the mechanical test problem of a cylindrical specimen of soft tissue in unconfined compression. These studies show that interfacial friction plays a major role in the unconfined compression response of articular cartilage specimens with small thickness to diameter ratios.

Biomechanical Phenomena↗