Search PubMed⌕ Search

Biomedical subjects

S C Cowin

Publications and source records attributed to S C Cowin.

At least 37 records · Page 2Linked to original sources

A model for the excitation of osteocytes by mechanical loading-induced bone fluid shear stresses.

A new experimentally testable hypothesis is advanced for the mechanosensory transduction mechanism by which communicating osteocytes sense the very small in vivo strains in the calcified matrix components of bone. We propose that the osteocytes, although not responsive to substantial fluid pressures, can be stimulated by relatively small fluid shear stresses acting on the membranes of their osteocytic processes. Biot's porous media theory is used to relate the combined axial and bending loads applied to a whole bone to the flow past the osteocytic processes in their canaliculi. In this theory, the bone pores of interest are the proteoglycan filled fluid annuli that surround the osteocytic processes in the canaliculi. We show that previously predicted fluid pore pressure relaxation times were a hundred-fold too short for the lacunar-canalicular porosity because they neglected the fluid drag associated with proteoglycan matrix on the surface membrane of the osteocyte and its cell processes. The recent theory developed in Tsay and Weinbaum [J. Fluid Mech. 226, 125-148 (1991)] for flow through cross-linked fiber filled channels is used to model the flow through this proteoglycan matrix. The predicted pore relaxation time, 1-2 s, closely corresponds to the times measured by Salzstein and Pollack [J. Biomechanics 20, 271-280 (1987)]. Furthermore, using this model, the magnitude of the predicted fluid induced shear stresses, 8-30 dyn cm-2, is shown to be similar to the fluid shear stresses measured in osteoblasts and other cells in which an intracellular Ca2+ shear stress response had been observed. This model is also used, in conjunction with anatomical data and the pore fluid pressure relaxation time data, to show that the spacing between the fibers is approximately 7 nm. The result is consistent with the notion that the canalicular pore space is filled with glycosaminoglycans that are ordered by albumin according to the model of Michel [J. Physiol. 404, 1-29 (1988)]. The new hypothesis is also shown to be consistent with the experiments of McLeod et al. [J. Biomechanics (submitted)] which suggest that high-frequency low-amplitude postural strains can maintain and even increase bone mass.

Animals↗

On the sufficiency conditions for the stability of bone remodeling equilibrium.

In this technical note a sufficiency condition is established for the stability of a strain-energy-based bone remodeling theory in the special case of a beam loaded by an axial force and a bending moment. In a previous report the same condition was shown to be a necessary condition for stability in the same situation. The remodeling scheme is one characterized by a remodeling stimulus equal to the strain energy density divided by the bulk or apparent density raised to an exponent m as well an elastic modulus proportional to bulk or apparent density raised to an exponent n. In order for a remodeling scheme to be stable for an elastic beam loaded by an axial force and a bending moment, it is established that the condition that m must be greater than n is not only necessary, but also sufficient.

Algorithms↗

Bone ingrowth: an application of the boundary element method to bone remodeling at the implant interface.

Surface bone remodeling theory and the boundary element method are employed to investigate the microstructural remodeling of bone at the bone-implant interface. Three situations are considered: remodeling-induced penetration between the screw threads of an implanted screw, penetration of bone tissue into a slot or cavity in an implant, and the interaction of individual trabeculae in the remodeling processes near an implant. For each case the bone ingrowth is determined as a function of the geometry and the applied load.

Animals↗

Chaos in the discrete-time algorithm for bone-density remodeling rate equations.

We compare the predictions of the differential equation form of a class of bone-density stress adaptation models with their associated discrete-time computational algorithms. Although our considerations apply to the class of adaptation models based on bulk or apparent bone-density remodeling, we focus attention on a particular model in this class, a model employed by Weinans et al. [Trans. Orthop. Res. Soc. 14, 310 (1989); Trans. First World Congress of Biomechanics, Vol. II, p. 75 (1990)]. We show that the discrete-time computational algorithm of that stress adaptation model has a well-known chaos mechanism for stress values of practical interest. Further, we obtain a condition on the discrete-time step that prevents the transition to chaos, and conditions that insure monotonic convergence. This chaos mechanism is only present in the discrete-time computational algorithm; we show that the corresponding differential equation form of the bone-density stress adaptation model is smooth, monotonic and nonchaotic.

Adaptation, Physiological↗

Bone stress adaptation models.

The basic concepts employed in formulating models of the process of stress adaptation in living bone tissue are reviewed. A purpose of this review is to define and separate issues in the formulation of bone remodeling theories. After discussing the rationale and objective of these models, the concepts and techniques involved in the modeling process are reviewed one by one. It is concluded that some techniques will be more successful than others in achieving the goals of computational bone remodeling. In particular, rationale is given for the preference of surface bone remodeling approaches over internal bone remodeling approaches, and for interactive multi-scale level, rather than mono-scale level, computational strategies.

Adaptation, Physiological↗

Periosteal and endosteal control of bone remodeling under torsional loading.

The shape changes that occur in the mid-diaphysis of a long bone due to adaptive remodeling induced by increasing or decreasing the axial and/or torsional loading of the bone are investigated using a simple model. In this model the mid-diaphysis of a long bone is represented as a hollow thick-walled right-circular cylinder, and different optimal strategies for bone remodeling are considered. It is shown that if such a thick-walled right-circular cylinder capable of surface remodeling is subjected to an axial compressive load and a twisting torque, then the remodeling patterns depend on whether the periosteal surface or the endosteal surface controls the limits of the remodeling process. It is shown that the effect of increasing the torque is always opposite to the effect of increasing the compressive load. Thus, similar remodeling patterns are obtained by increasing one type of loading and decreasing the other. Aside from the restriction of idealized cylindrical geometry, the only assumptions made are that the bone tissue is linearly elastic and that there exists a finite range of remodeling equilibrium stresses. Only those loading situations which maintain the bone in remodeling equilibrium are considered in this work. It follows that the results presented are independent of the specific type of rule governing the temporal evolution of the bone shape, since any such rule applies only in situations where there is active remodeling and, hence, no remodeling equilibrium.

Biomechanical Phenomena↗

An evolutionary Wolff's law for trabecular architecture.

A continuum model is proposed to describe the temporal evolution of both the density changes and the reorientation of the trabecular architecture given the applied stress state in the bone and certain material parameters of the bone. The data upon which the proposed model is to be based consist of experimentally determined remodeling rate coefficients and quantitative stereological and anisotropic elastic constant measurements of cancellous bone. The model shows that the system of differential equations governing the temporal changes in architecture is necessarily nonlinear. This nonlinearity is fundamental in that it stems from the fact that, during remodeling, the relationship between stress and strain is changing as the stress and strain variables themselves are changing. In order to preserve the remodeling property of the model, terms that are of the order strain times the changes in density and/or microstructural properties must be retained. If these terms were dropped, there would be no feedback mechanism for architectural adaptation and no adaptation of the trabecular architecture. There is, therefore, no linearized version of the model of the temporal evolution of trabecular architecture. An application of the model is illustrated by an example problem in which the temporal evolution of homogeneous trabecular architecture is predicted. A limitation of the proposed continuum model is the length scale below which it cannot be applied. The model cannot be applied in regions of cancellous bone where the trabecular bone architecture is relatively inhomogeneous or at a bone-implant interface.

Bone Density↗

Correction formulae for the misalignment of axes in the measurement of the orthotropic elastic constants.

In the experimental determination of the orthotropic elastic constants, one often encounters the situation in which the symmetry axes of the material are not coincident with specimen axes along which the material testing is accomplished. The problem of calculating the compliance coefficients in the symmetry coordinate system from measurements of the compliance coefficients made in an arbitrary, specimen fixed, coordinate system is considered here.

Biomechanical Phenomena↗

Non-interacting modes for stress, strain and energy in anisotropic hard tissue.

The six non-interacting modes for stress, strain and energy in an orthotropic elastic model of human femoral cortical bone tissue are discussed and illustrated. The stress and strain modes are illustrated using the representation of the stress and strain fields around a circular hole in a flat plate of cortical bone subjected to a uniaxial field of tension as the example. The six modes play a role in the stress analysis of orthotropic elastic materials similar to the roles played by the hydrostatic and deviatoric non-interacting stress, strain and energy modes in isotropic elasticity. The biomechanical significance of the six non-interacting modes for stress, strain and energy in hard tissue is both practical and suggestive. The modes suggest a practical scheme for the representation of stress and strain fields in hard tissue. The existence of the modes suggests physical insights, for example, possible failure mechanisms or adaptation strategies possessed by the hard tissues.

Anisotropy↗

Candidates for the mechanosensory system in bone.

Some potential mechanisms by which bone cells sense mechanical loads are described and hypotheses concerning the functioning of these mechanisms are explored. It is well known that bone tissue adapts its structure to its mechanical load environment. Recent research has illuminated the biological response of bone to mechanical loading at the cellular level, but the precise mechanosensory system that signals bone cells to deposit or resorb tissue has not been identified. The purpose of this paper is to describe the current status of this research and to suggest some possible mechanosensory systems by which bone cells might sense environmental loads.

Adaptation, Physiological↗

The fabric dependence of the orthotropic elastic constants of cancellous bone.

It has been proposed that the orthotropic elastic constants of cancellous bone depend upon a tensorial measure of anisotropy called fabric as well as the tissue's structural density. Cowin (1985, Mechanics Mater, 4, 137-147; 1986, J. biomech. Engng 108, 83-88) developed explicit relationships for the elastic constant, structural density and fabric relationship. In this study the orthotropic elastic moduli, structural density, and fabric components were measured for 11 cancellous bone specimens from five bovine femora and for 75 specimens from three human proximal tibiae and fitted to these relationships using a least squares analysis. The relationships explained between 72 and 94% of the variance in the elastic constants. The relationships between the elastic constants and squared or cubed power functions of structural density had better predictive value over the entire distribution of the data than did expressions with linear functions of structural density.

Aged↗

Errors in the orientation of the principal stress axes if bone tissue is modeled as isotropic.

The error in the prediction of the orientation of the principal axes of stress in bone tissue is determined in the case when the tissue is modeled as elastically isotropic rather than as orthotropic, the probable symmetry of bone tissue. Results are two-dimensional and assume the same underlying strain state for both the orthotropic and isotropic cases. The maximum error is 45 degrees, and the typical error is generally significant.

Biomechanical Phenomena↗

Deviatoric and hydrostatic mode interaction in hard and soft tissue.

It has been established that many hard and soft tissues have anisotropic material symmetry. It is noted here that the deviatoric and hydrostatic modes interact with each other in a general anisotropic elastic material. In the special case of isotropic, linear elastic, materials these modes are non-interactive. As a consequence of the interaction of these modes encountered in anisotropic materials, the decomposition into hydrostatic and deviatoric modes, and deviatoric mode concepts such as the von Mises effective stress are not appropriate for anisotropic materials in general. The implications of this observation for the presentation of computationally generated stress contours for hard and soft tissues are discussed. It is also pointed out that the mode coupling and mode interaction raise the question of whether anisotropic living tissues respond directly to stress or to some other physical quantity such as strain or strain energy, in view of the recent hypothesis concerning the proliferation and ossification of cartilage.

Adaptation, Physiological↗

Identification of the elastic symmetry of bone and other materials.

A simplified classification scheme for the elastic symmetries of a solid is applied to the identification of the elastic symmetry of a material by three different methods--visual, stereological and numerical algorithm. Each method is illustrated with an application to bone tissues, but the methods apply to all materials.

Algorithms↗

A resolution restriction for Wolff's law of trabecular architecture.

A rigid constructionist view of Wolff's trajectorial theory of trabecular architecture, called the "orthodox theory," requires that the trabeculae of cancellous bone intersect at right angles exactly like the stress trajectories with which they are associated by the theory. It is well established that trabeculae do not always intersect at right angles. The author discusses the fallacy in the orthodox interpretation of Wolff's law, and suggests that a "resolution length restriction" be imposed on the trajectorial theory to avoid interpretations that lead to the fallacy.

Bone Development↗

On the dependence of the elasticity and strength of cancellous bone on apparent density.

This paper presents a statistical analysis of the pooled data from a number of previous experiments concerning the dependence of the Young's moduli and strength of cancellous bone tissue upon apparent density. The results show that both the Young's moduli and the strength are proportional to the square of apparent density of the tissue and are therefore proportional to one another. It is shown that the coefficient of proportionality is different for human and bovine tissue. It is concluded that the suggestion of Wolff (Das Gesetz der Transformation der Knochen, Hirschwald, Berlin, 1892) that compact bone tissue is simply more dense cancellous bone tissue is not an accurate statement when only the mechanical properties of these two tissues are considered. It is noted that estimates for the elastic modulus of the individual trabecula of human cancellous bone vary from 1 to 20 GPa and it is suggested that this question needs further study.

Biomechanical Phenomena↗

Errors induced by off-axis measurement of the elastic properties of bone.

Misalignment between the axes of measurement and the material symmetry axes of bone causes error in anisotropic elastic property measurements. Measurements of Poisson's ratio were strongly affected by misalignment errors. The mean errors in the measured Young's moduli were 9.5 and 1.3 percent for cancellous and cortical bone, respectively, at a misalignment angle of 10 degrees. Mean errors of 1.1 and 5.0 percent in the measured shear moduli for cancellous and cortical bone, respectively, were found at a misalignment angle of 10 degrees. Although, cancellous bone tissue was assumed to have orthotropic elastic symmetry, the possibility of the greater symmetry of transverse isotropy was investigated. When the nine orthotropic elastic constants were forced to approximate the five transverse isotropic elastic constants, errors of over 60 percent were introduced. Therefore, it was concluded that cancellous bone is truly orthotropic and not transversely isotropic. A similar but less strong result for cortical bone tissue was obtained.

Animals↗