PubMed Health⌕ Search

Biomedical subjects

J M Huyghe

Publications and source records attributed to J M Huyghe.

At least 19 recordsLinked to original sources

Depth-dependent compressive equilibrium properties of articular cartilage explained by its composition.

For this study, we hypothesized that the depth-dependent compressive equilibrium properties of articular cartilage are the inherent consequence of its depth-dependent composition, and not the result of depth-dependent material properties. To test this hypothesis, our recently developed fibril-reinforced poroviscoelastic swelling model was expanded to include the influence of intra- and extra-fibrillar water content, and the influence of the solid fraction on the compressive properties of the tissue. With this model, the depth-dependent compressive equilibrium properties of articular cartilage were determined, and compared with experimental data from the literature. The typical depth-dependent behavior of articular cartilage was predicted by this model. The effective aggregate modulus was highly strain-dependent. It decreased with increasing strain for low strains, and increases with increasing strain for high strains. This effect was more pronounced with increasing distance from the articular surface. The main insight from this study is that the depth-dependent material behavior of articular cartilage can be obtained from its depth-dependent composition only. This eliminates the need for the assumption that the material properties of the different constituents themselves vary with depth. Such insights are important for understanding cartilage mechanical behavior, cartilage damage mechanisms and tissue engineering studies.

Cartilage, Articular↗

A composition-based cartilage model for the assessment of compositional changes during cartilage damage and adaptation.

OBJECTIVE: The composition of articular cartilage changes with progression of osteoarthritis. Since compositional changes are associated with changes in the mechanical properties of the tissue, they are relevant for understanding how mechanical loading induces progression. The objective of this study is to present a computational model of articular cartilage which enables to study the interaction between composition and mechanics. METHODS: Our previously developed fibril-reinforced poroviscoelastic swelling model for articular cartilage was combined with our tissue composition-based model. In the combined model both the depth- and strain-dependencies of the permeability are governed by tissue composition. All local mechanical properties in the combined model are directly related to the local composition of the tissue, i.e., to the local amounts of proteoglycans and collagens and to tissue anisotropy. RESULTS: Solely based on the composition of the cartilage, we were able to predict the equilibrium and transient response of articular cartilage during confined compression, unconfined compression, indentation and two different 1D-swelling tests, simultaneously. CONCLUSION: Since both the static and the time-dependent mechanical properties have now become fully dependent on tissue composition, the model allows assessing the mechanical consequences of compositional changes seen during osteoarthritis without further assumptions. This is a major step forward in quantitative evaluations of osteoarthritis progression.

Adaptation, Physiological↗

A comparison between mechano-electrochemical and biphasic swelling theories for soft hydrated tissues.

Biological tissues like intervertebral discs and articular cartilage primarily consist of interstitial fluid, collagen fibrils and negatively charged proteoglycans. Due to the fixed charges of the proteoglycans, 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). This excess of ion particles leads to an osmotic pressure difference, which causes swelling of the tissue. In the last decade several mechano-electrochemical models, which include this mechanism, have been developed. As these models are complex and computationally expensive, it is only possible to analyze geometrically relatively small problems. Furthermore, there is still no commercial finite element tool that includes such a mechano-electrochemical theory. Lanir (Biorheology, 24, pp. 173-187, 1987) hypothesized that electrolyte flux in articular cartilage can be neglected in mechanical studies. Lanir's hypothesis implies that the swelling behavior of cartilage is only determined by deformation of the solid and by fluid flow. Hence, the response could be described by adding a deformation-dependent pressure term to the standard biphasic equations. Based on this theory we developed a biphasic swelling model. The goal of the study was to test Lanir's hypothesis for a range of material properties. We compared the deformation behavior predicted by the biphasic swelling model and a full mechano-electrochemical model for confined compression and 1D swelling. It was shown that, depending on the material properties, the biphasic swelling model behaves largely the same as the mechano-electrochemical model, with regard to stresses and strains in the tissue following either mechanical or chemical perturbations. Hence, the biphasic swelling model could be an alternative for the more complex mechano-electrochemical model, in those cases where the ion flux itself is not the subject of the study. We propose thumbrules to estimate the correlation between the two models for specific problems.

Animals↗

An ionised/non-ionised dual porosity model of intervertebral disc tissue.

The volume of the intrafibrillar water space--i.e. the water contained inside the collagen fibres--is a key parameter that is relevant to concepts of connective tissue structure and function. Confined compression and swelling experiments on annulus fibrosus samples are interpreted in terms of a dual porosity model that distinguishes between a non-ionised intrafibrillar porosity and an ionised extrafibrillar porosity. Both porosities intercommunicate and are saturated with a monovalent ionic solution, i.c. NaCl. The extrafibrillar fixed charge density of the samples is assessed using radiotracer techniques and the collagen content is evaluated by measurement of hydroxyproline concentration. The interpretation of the experimental data yields values for the intrafibrillar water content, the average activity coefficient of the ions, the Donnan osmotic coefficient, the fraction of intrafibrillar water, the stress-free deformation state, and an effective stress-strain relationship as a function of the radial position in the disc. A linear fit between the second Piola-Kirchhoff effective stress and Green-Lagrange strain yielded an effective stiffness: H(e)=1.087 +/- 0.657 MPa. The average fraction of intrafibrillar water was 1.16 g/g collagen. The results were sensitive to changes in the activity and osmotic coefficients and the fraction of intrafibrillar water. The fixed charge density increased with distance from the outer edge of the annulus, whereas the hydroxyproline decreased.

Animals↗

Remodelling of continuously distributed collagen fibres in soft connective tissues.

Extracellular matrix remodelling plays an essential role in tissue engineering of load-bearing structures. The goal of this study is to model changes in collagen fibre content and orientation in soft connective tissues due to mechanical stimuli. A theory is presented describing the mechanical condition within the tissue and accounting for the effects of collagen fibre alignment and changes in fibre content. A fibre orientation tensor is defined to represent the continuous distribution of collagen fibre directions. A constitutive model is introduced to relate the fibre configuration to the macroscopic stress within the material. The constitutive model is extended with a structural parameter, the fibre volume fraction, to account for the amount of fibres present within the material. It is hypothesised that collagen fibre reorientation is induced by macroscopic deformations and the amount of collagen fibres is assumed to increase with the mean fibre stretch. The capabilities of the model are demonstrated by considering remodelling within a biaxially stretched cube. The model is then applied to analyse remodelling within a closed stented aortic heart valve. The computed preferred fibre orientation runs from commissure to commissure and resembles the fibre directions in the native aortic valve.

Adaptation, Physiological↗

Numerical simulation of deformations and electrical potentials in a cartilage substitute.

Cartilage exhibits a swelling and shrinking behaviour that influences the function of the cells inside the tissue. This behaviour is caused by mechanical, chemical and electrical loads. It is described by the electrochemomechanical mixture theory, in which the tissue is represented by four components: a charged porous solid, a fluid, cations and anions. By distinguishing between the cations and anions, electrical phenomena can be modelled. This mixture theory is verified by fitting the deformations and the electrical potentials in a uniaxial confined swelling and compression experiment to a mixed finite element simulation. The fitted stiffness, permeability, diffusion coefficients, and osmotic coefficients are in the same range as reported in literature.

Cartilage↗

Measuring principles of frictional coefficients in cartilaginous tissues and its substitutes.

The frictional properties of cartilaginous tissues, such as the hydraulic permeability, the electro-osmotic permeability, the diffusion coefficients of various ions and solutes, and the electrical conductance, are vital data to characterise the extracellular environment in which chondrocytes reside. This paper analyses one-dimensional measurement principles of these coefficients. Particular attention is given to the deformation dependence of them and the highly deformable nature of the tissues. A suggested strategy is the combination of a diffusion experiment using radiotracer methods, an electro-osmotic flow experiment and an electro-osmotic pressure experiment at low electric current.

Animals↗

Spatial interaction between tissue pressure and skeletal muscle perfusion during contraction.

The vascular waterfall theory attributes decreased muscle perfusion during contraction to increased intramuscular pressure (P(IM)) and concomitant increase in venous resistance. Although P(IM) is distributed during contractions, this theory does not account for heterogeneity. This study hypothesises that pressure heterogeneity could affect the interaction between P(IM) rise and perfusion. Regional tissue perfusion during submaximum (100kPa) tetanic contraction is studied, using a finite element model of perfused contracting skeletal muscle. Capillary flow in muscles with one proximal artery and vein (SIM(1)) and with an additional distal artery and vein (SIM(2)) is compared. Blood flow and pressures at rest and P(IM) during contraction ( approximately 25kPa maximally) are similar between simulations, but capillary flow and venous pressure differ. In SIM(2), venous pressure and capillary flow correspond to P(IM) distribution, whereas capillary flow in SIM(1) is less than 10% of flow in SIM(2), in the muscle half without draining vein. This difference is caused by a high central P(IM), followed by central venous pressure rise, in agreement with the waterfall theory. The high central pressure (SIM(1)), obstructs outflow from the distal veins. Distal venous pressure rises until central blood pressure is reached, although local P(IM) is low. Adding a distal vein (SIM(2)) restores the perfusion. It is concluded that regional effects contribute to the interaction between P(IM) and perfusion during contraction. Unlike stated by the vascular waterfall theory, venous pressure may locally exceed P(IM). Although this can be explained by the principles of this theory, the theory does not include this phenomenon as such.

Animals↗

Requirements for an artificial intervertebral disc.

Intervertebral disc degeneration is an important social and economic problem. Presently available artificial intervertebral discs (AIDs) are insufficient and the main surgical intervention is still spinal fusion. The objective of the present study is to present a list of requirements for the development of an AID which could replace the human lumbar intervertebral disc and restore its function. The list addresses geometry, stiffness, range of motion, strength, facet joint function, center of rotation, fixation, failsafety and implantation technique. Date are obtained from the literature, quantified where possible and checked for consistency. Existing AIDs are evaluated according to the presented list of requirements. Endplate size is a weak point in existing AIDs. These should be large and fit vertebral bodies to prevent migration. Disc height and wedge angle should be restored, unless this would overstretch ligaments. Finally, stiffness and range of motion in all directions should equal those of the healthy disc, except for the axial rotation to relieve the facet joints.

Biomechanical Phenomena↗

Mechanical blood-tissue interaction in contracting muscles: a model study.

A finite element (FE) model of blood perfused biological tissue has been developed. Blood perfusion is described by fluid flow through a series of 5 intercommunicating vascular compartments that are embedded in the tissue. Each compartment is characterized by a blood flow permeability tensor, blood volume fraction and vessel compliance. Local non-linear relationships between intra-extra vascular pressure difference and blood volume fraction, and between blood volume fraction and the permeability tensor, are included in the FE model. To test the implementation of these non-linear relations, FE results of blood perfusion in a piece of tissue that is subject to increased intramuscular pressure, are compared to results that are calculated with a lumped parameter (LP) model of blood perfusion. FE simulation of blood flow through a contracting rat calf muscle is performed. The FE model used in this simulation contains a transversely isotropic, non-linearly elastic description of deforming muscle tissue, in which local contraction stress is prescribed as a function of time. FE results of muscle tension, total arterial inflow and total venous outflow of the muscle during contraction, correspond to experimental results of an isometrically and tetanically contracting rat calf muscle.

Animals↗

Partition and diffusion of sodium and chloride ions in soft charged foam: the effect of external salt concentration and mechanical deformation.

The partition and diffusion characteristics of an acrylic acid/acrylamide hydrogel, copolymerized in the pores of a polyurethane foam with sodium and chloride ions, were studied by radiochemical methodologies. The hydrogel foam swells by 51%, 80%% and 260% relative to its raw state under bath salt concentrations of 2.0, 1.0, and 0.15 M, respectively. The corresponding partition coefficients are 1.13, 1.29, and 1.99 for sodium (Na+) and 0.89, 0. 85, and 0.65 for chloride (Cl-). The diffusion coefficients are independent of bath concentration and increase linearly with hydration towards their values in water. Deformation affects partition and diffusion solely by dilatation, which determines the swelling and hydration. Comparison of the hydrogel foam with cartilage and intervertebral disc shows considerable similarities and suggests that the same mechanisms control their function.

Acrylamide↗

Nonhomogeneous permeability of canine anulus fibrosus.

STUDY DESIGN: This report examines the permeability coefficient and aggregate modulus of slices of anulus cut from canine lumbar intervertebral discs. OBJECTIVES: To examine the influence of radial position on the properties of these materials, including outer samples with intact anulus edge. SUMMARY OF BACKGROUND DATA: The outer edge of anulus fibrosus shows radial bulge during axial compression of motion segments. The radial bulge increases monotonically when the axial compression is sustained for several hours, until a plateau is reached. Triphasic modeling of axial compression shows that this time course of radial bulge can not be obtained using a uniform permeability coefficient according to values in the literature. METHODS: Confined consolidation experiments (controlled load) were designed to measure the time course of uniaxial deformation of samples of anulus that were 4 mm in diameter and 1 mm tall. The rotation symmetry axis of the samples was defined in the radial direction of the disc. The radial permeability coefficient and the aggregate modulus were determined using the consolidation data and the linear biphasic theory. RESULTS: The permeability coefficient was lower at the periphery than in deeper layers of the anulus. Outer samples with outer surfaces that were 0.0-0.5 mm from the anulus edge had an average permeability coefficient of (1.02 +/- 0.57) x 10(-16) m4/Ns (n = 24). Inner samples that were 2.0-2.5 mm from the anulus edge had an average permeability coefficient of (2.81 +/- 0.98) x 10(-16) m4/Ns (n = 13). The aggregate modulus HA of outer samples was significantly higher (HA = 1.56 +/- 0.34 MPa) than that of inner samples (HA = 1.31 +/- 0.47 MPa). CONCLUSIONS: The fact that the outer anulus is less permeable than the inner anulus may explain why radial bulge of anulus fibrosus increases monotonically in time to an equilibrium value during sustained axial compression of a motion segment.

Animals↗

Finite-element simulation of blood perfusion in muscle tissue during compression and sustained contraction.

Mechanical interaction between tissue stress and blood perfusion in skeletal muscles plays an important role in blood flow impediment during sustained contraction. The exact mechanism of this interaction is not clear, and experimental investigation of this mechanism is difficult. We developed a finite-element model of the mechanical behavior of blood-perfused muscle tissue, which accounts for mechanical blood-tissue interaction in maximally vasodilated vasculature. Verification of the model was performed by comparing finite-element results of blood pressure and flow with experimental measurements in a muscle that is subject to well-controlled mechanical loading conditions. In addition, we performed simulations of blood perfusion during tetanic, isometric contraction and maximal vasodilation in a simplified, two-dimensional finite-element model of a rat calf muscle. A vascular waterfall in the venous compartment was identified as the main cause for blood flow impediment both in the experiment and in the finite-element simulations. The validated finite-element model offers possibilities for detailed analysis of blood perfusion in three-dimensional muscle models under complicated loading conditions.

Animals↗

Strain distribution on rat medial gastrocnemius (MG) during passive stretch.

Deformation of the surface of passive medial gastrocnemius muscle (MG) was measured in vivo while performing a hysteresis test. The gastrocnemius muscle of male rats were dissected free and the distal tendon was cut. The lateral head was separated from the medial head. The muscle origins were left intact. 60-70 fluorescent, polystyrene spheres (diameter 0.7 mm) were attached to the surface of the MG. During the experiment, two-dimensional video recordings of the movements of the MG were made. The coordinates of the marker centroids were obtained by computer processing of digitized images and marker displacements as a function of time were calculated. Green-Lagrange strains in two principal directions were calculated (epsilon 1, epsilon 2) for three specimens. epsilon 1 had approximately the same direction as the muscle fibers. The longitudinal strain of the fibers (20-30%) was larger than the strain of the aponeurosis (1-5%); p < 0.001. No significant difference was found between the values of the transverse strains of muscle fibers and aponeurosis; the value of epsilon 2 was -6 to -9% for both tissue structures.

Analysis of Variance↗

A 3-D finite element model of blood perfused rat gastrocnemius medialis muscle.

A finite element description of blood perfusion has been developed, and is applied to skeletal muscles. Three-dimensional distributions of blood pressures and flows in deforming muscles are calculated. The muscle tissue is considered as a fluid-saturated porous solid. The blood is modeled as a series of five intercommunicating compartmental fluids, representing arterial, arteriolar, capillary, venular and venous blood, that reside in the pores (blood vessels) of the muscle tissue. The blood vessels are modeled as distensible tubes, embedded in the muscle tissue. A 3-D finite element mesh has been mapped on a reconstructed geometry of a gastrocnemius medialis muscle of the rat. Blood perfused linear elastic muscle material behaviour has been assigned to this mesh. A simulation of blood perfusion, resulting from a constant arterio-venous pressure difference, through the reconstructed muscle has been performed. Calculated blood pressure and flow distributions were within physiological range.

Animals↗

Regional wall mechanics in the ischemic left ventricle: numerical modeling and dog experiments.

The mechanics of the ischemic left ventricle during a complete cardiac cycle were simulated using a finite-element model accounting for the thick-walled ventricular geometry, the fibrous nature of the myocardial tissue, and the dependency of active muscle fiber stress on time, strain, and strain rate. Ischemia was modeled by disabling the generation of active stress in a region comprising approximately 12% of total wall volume. In the model simulations, the approximately 12% reduction in the amount of normally contracting tissue resulted in an approximately 25% reduction in stroke work compared with the normal situation. The more-than-proportional loss of stroke work may partly be attributed to storage of elastic energy in the bulging ischemic region. Furthermore the mechanical performance in the nonischemic border zone deteriorated because of reduced systolic fiber stress (if fibers were in series with those in the ischemic region) or reduced fiber shortening (if fibers were parallel). The deformation pattern of the ventricle was asymmetric with respect to the ischemic region because of the anisotropy of the myocardial tissue. Epicardial fiber shortening in and around the ischemic region, as predicted from the model simulations, was in qualitative agreement with shortening, as measured in four dogs in which ischemia was induced by occlusion of the distal part of the left anterior interventricular coronary artery.

Animals↗