| Sign In to gain access to subscriptions and/or personal tools. |
© 2003 International and American Associations for Dental Research THE MECHANICAL PROPERTIES OF HUMAN DENTIN: A CRITICAL REVIEW AND RE-EVALUATION OF THE DENTAL LITERATUREDivision of Biomaterials and Bioengineering, Department of Preventive and Restorative Dental Sciences, Mail Stop 0758, University of California, San Francisco, San Francisco, CA 94143-0758; Correspondence: *corresponding author, kinney3{at}llnl.gov
The past 50 years of research on the mechanical properties of human dentin are reviewed. Since the body of work in this field is highly inconsistent, it was often necessary to re-analyze prior studies, when possible, and to re-assess them within the framework of composite mechanics and dentin structure. A critical re-evaluation of the literature indicates that the magnitudes of the elastic constants of dentin must be revised considerably upward. The Youngs and shear moduli lie between 20-25 GPa and 7-10 GPa, respectively. Viscoelastic behavior (time-dependent stress relaxation) measurably reduces these values at strain rates of physiological relevance; the reduced modulus (infinite relaxation time) is about 12 GPa. Furthermore, it appears as if the elastic properties are anisotropic (not the same in all directions); sonic methods detect hexagonal anisotropy, although its magnitude appears to be small. Strength data are re-interpreted within the framework of the Weibull distribution function. The large coefficients of variation cited in all strength studies can then be understood in terms of a distribution of flaws within the dentin specimens. The apparent size-effect in the tensile and shear strength data has its origins in this flaw distribution, and can be quantified by the Weibull analysis. Finally, the relatively few fracture mechanics and fatigue studies are discussed. Dentin has a fatigue limit. For stresses smaller than the normal stresses of mastication, 30 MPa, a flaw-free dentin specimen apparently will not fail. However, a more conservative approach based on fatigue crack growth rates indicates that if there is a pre-existing flaw of sufficient size ( 0.3-1.0 mm), it can grow to catastrophic proportion with cyclic loading at stresses below 30 MPa.
Key Words: Dentin calcified tissues mechanical properties fatigue fracture toughness
Dentin is the most abundant mineralized tissue in the human tooth. Therefore, knowledge of its mechanical properties is essential for predicting the effects of microstructural alterations due to caries, sclerosis, and aging on tooth strength. In clinical dentistry, knowledge of dentin properties is important for understanding the effects of the wide variety of restorative dental procedures that range from the design of preparations to choice of bonding methods. In spite of this importance, over half a century of research has failed to provide consistent values of dentins mechanical properties. The Youngs modulus is unknown to within a factor of three; the shear modulus is uncertain by a factor of four; and the other elastic constants have not been measured. There is a six-fold uncertainty in the ultimate strength, and there have been few studies of fracture toughness or fatigue. This is a review of the last 50+ years of research into the mechanical properties of dentin. It is a critical review, in the sense that the prior literature on the mechanical properties of dentin is often re-examined and re-interpreted within the context of dentin microstructure. For example, the widely varying tensile strength data are re-examined in terms of the Weibull distribution function; the wide variations in previously reported tensile strengths can then be explained by a specimen size effect whose origin lies in the existence of a population of flaws in the dentin. The bulk of this review focuses on the elastic properties of dentin. This emphasis is a necessary consequence of both the paucity of published data on other mechanical properties, and the fact that an understanding of elastic behavior is essential to the proper interpretation of physical measurements of failure. It was often necessary to re-evaluate the data on the elastic properties to reconcile the contradictions in the literature and reach consensus. Because the physical data were frequently not available in the more recent publications, indirect methods were used to test the meaningfulness of the results. This involved consideration of the entire elastic stiffness matrix to check the constraints forced upon all of the elastic constants from the measurement of a single property. In this way, it was possible to check for self-consistency with all of the elastic constants as well as the known microstructure. After a brief discussion of the composition and microstructure of dentin, the review first considers its elastic and viscoelastic properties. Then, hardness, strength, fracture toughness, and fatigue are discussed in order. This organization loosely follows that of the last comprehensive review of dentin properties (Waters, 1980), with the exception that a more thorough treatment of fatigue and fracture toughness has been added. The review concludes by proposing reliable ranges for the magnitudes of mechanical properties. It is hoped that these recommendations, and the evidence on which they are based, spark additional discourse and research on the properties of dentin.
The microstructure of dentin suggests a hierarchical approach to the understanding of its mechanical properties. At the smallest length scales are the constituent materials: a carbonated nanocrystalline apatite mineral phase (approximately 50% by volume) and a felt-work of type I collagen fibrils (see Fig. 1 5 nm, is invariant with location.
At a higher level of organization is an intermediate, or composite, length scale. At this length scale, dentin can be modeled as a continuous fiber-reinforced composite, with the intertubular dentin forming the matrix and the tubule lumens with their associated cuffs of peritubular dentin forming the cylindrical fiber reinforcement (see Fig. 2
At the greatest length scale are the effective, or continuum, properties of dentin. The effective properties of dentin describe the response of the tooth to applied loads, and allow for predictions of tooth strength and fracture properties. Youngs modulus, tensile and compressive strength, and fracture toughness are examples of these effective properties, and reflect the complex interactions of the constituent materials and the microstructure. At this largest length scale, we anticipate that the effective properties will depend on tubule density, orientation, and the local average density of the mineral phase. Ultimately, measurement of the effective properties, particularly the properties of altered forms of dentin (such as carious or transparent), will benefit from an understanding of the material behavior at all of these length scales.
The elastic properties of dentin are of paramount importance in all discussions of tooth strength. The elastic constants, usually defined in terms of a stiffness (Cij) or compliance (Sij) matrix, include the Youngs modulus, shear modulus, and Poissons ratio. Depending on the symmetry of the microstructure, the elastic constants have different degrees of independence. For example, an isotropic structure has only two independent elastic constants, while an orthotropic structure has nine independent elastic constants. Therefore, any study of the elastic properties of dentin must consider its symmetry.
THE YOUNGS MODULUS
Constituent materials models An early attempt to provide reasonable bounds for the continuum Youngs modulus of dentin was proposed by Katz (1971). Dentin was considered to be a two-phase mixture of apatite mineral and collagen, with the volume fraction and Youngs modulus of the mineral phase given by VA and EA, respectively, and the corresponding values for the collagen given by VC and EC. The effective Youngs modulus of the mixture (Ee) was rigorously bounded from above and below by:
The elastic properties of the constituent materials in their bulk form have been measured with reasonable accuracy: EA
Micromechanics models Because of the large difference in modulus between the mineral and organic phases, we do not anticipate that improvement can be made by further refinements of bounds. Therefore, it is worthwhile considering a micromechanics approach at the composite hierarchy of dentin microstructure. At this length scale, the dentinal tubules with their surrounding peritubular cuffs are like fibers, providing reinforcement to the intertubular matrix. Assuming that the moduli of the matrix (EM) and the "fibers" (the lumens plus peritubular cuffs, EF) are known, rough estimates of the bounds can be made by appropriate substitutions into Eqs. 1 and 2. However, direct substitution of EM and EF ignores the possibly complex interactions between the "fibers" and the matrix. The interactions between fiber and matrix have been accounted for in self-consistent generalized micromechanics approaches to aligned fiber composites (Christensen, 1990). Analytical solutions have been derived for the case of aligned cylindrical fibers in an isotropic matrix. The fibers introduce a transverse anisotropy and require five equations for the independent elastic constants to be defined. These equations have been reformulated for dentin in terms of the Youngs modulus of intertubular (Ei) and peritubular (Ept) dentin (Kinney et al., 1999). The expression for the effective longitudinal (in the direction of the tubules) Youngs modulus, Eel, in terms of the tubule concentration, c, was derived:
In the above equation, Only when the ratio of the peritubular and intertubular dentin moduli is greater than 3 is there a measurable effect of the tubules on stiffness for physiologically relevant tubule densities (5-15%). Thus, a theoretical framework now exists for estimating the effects of tubule orientation on the elastic properties: For physiologically relevant tubule densities, the anisotropy caused by the tubules is insignificant (Kinney et al., 1999). This does not mean that dentin is elastically isotropic; rather, if dentin proves to be anisotropic, it is because the intertubular dentin matrix is anisotropic and not because of the tubules. Based on observations that the collagen fibrils in the intertubular matrix are aligned perpendicular to the tubule axis, we anticipate that dentin might have a transverse isotropic (hexagonal) symmetry (stiffer in the plane of the collagen fibrils). The micromechanics model offers several advantages over bounding methods. Key is the ability of the model to make accurate estimates of the effective elastic constants based on a few simple measurements of the tubule concentration. The model has its drawbacks, however. Accurate values of the elastic constants of peritubular and intertubular dentin must be available, and this requires careful mechanical testing, which we will now review in greater detail.
Tensile and compressive measurements of Youngs modulus The earliest "modern" measurement of Youngs modulus in tension was published in 1962 (Bowen and Rodriguez, 1962). In that study, specimens of dentin, nominally 2 mm thick, were cut freehand with a rotary diamond disk under constant water irrigation. The specimens were narrowed in the middle, with a "radiused" shoulder (fillet) maintained so that stress concentration would be prevented. The knife-edges of an optical strain gauge (length, 6 mm) were attached directly to the specimens. The specimens were hydrated at least 1 hr in distilled water before being tested. The stress-strain curves were linear to failure, and the mean modulus of elasticity was 19.3 GPa (2.8 x 106 psi) with a coefficient of variation of 28%. Five years later, Lehman (1967) measured the tensile properties of dentin with hollowed, cylindrical specimens from the root. As in the earlier study by Bowen and Rodriguez, the tensile stress-strain curves were linear to failure. However, the Youngs modulus was 11.0 GPa, almost half that of the earlier study, and the coefficient of variation (53%) was almost double. There were three important differences between the two studies. First, the latter study did not use strain gauges affixed to the specimen, thereby increasing the likelihood of grip effects like tow-in. Second, the specimens were hollowed along the root canal by means of a dental bur, increasing the probability of undetected flaws in the interior of the specimen. Third, the gauge sections were not "radiused", leaving a stress concentration that increased the chance of failure at grip ends, thereby invalidating the test. Though it is difficult to reconstruct an experiment decades after the fact, we are fortunate that Lehman reported all of his data. A close inspection shows that the modulus and ultimate tensile strength were strongly correlated (p < 0.001), with the higher values of elastic modulus associated with the specimens that had the highest tensile strength. Furthermore, 75% of the specimens failed in tension at or below 40 MPa, a low value suggesting that flaws introduced during specimen preparation may have affected the linear stress-strain behavior. When we restrict our focus to those specimens that failed above 40 MPa, we obtain a tensile modulus of 16.9 GPa with a coefficient of variation of 26%. This is more in line with the earlier work of Bowen and Rodriguez (1962). Around the same time that the tensile measurements were being performed, other researchers were establishing the compressive properties of dentin. The earliest of these measurements was reported by Peyton et al. (1952). This study examined the compressive behavior of 1.8 mm in cross-section by 4.5-mm-long dentin specimens. Strain gauges were affixed to steel rods that were then used to apply load to the specimens. A Youngs modulus of 11.6 GPa was reported. Concerned that the low value of the Youngs modulus obtained from their earlier study might have been caused by a combination of non-parallel alignment of the load platens (tow-in) and possible stress relaxation effects, the same group repeated their earlier study with greater attention to experimental variables (Craig and Peyton, 1958). Placing the strain gauges directly on the specimen eliminated tow-in; stress relaxation effects were measured by careful cycling of the compressive load well below the proportional limit. These corrections raised the compressive Youngs modulus to 18.5 GPa, in excellent agreement with the tensile measurements of the time.
Measurements of the modulus are sensitive to specimen preparation, experimental design, and stress relaxation; not accounting for these experimental variables leads to underestimation of the elastic constants. Therefore, the comments of Waters (1980) in his review article were surprising: "Craig and Peyton (1958) obtained values for the proportional limit and compressive strength in reasonable agreement with other workers, but their value of the modulus is, for some unaccountable reason, considerably higher." This comment reflects a bias toward favoring a low value for the Youngs modulus of dentin—a bias perhaps enforced by the large number of low values reported in subsequent years (see Fig. 3 Testing of dentin in compression and tension was performed only infrequently in the years since Craig and Peyton. Two of these studies are particularly noteworthy. In the work reported by Stanford et al. (1960), dentin specimens were tested in compression. Even with corrections for platen deformation, the compressive modulus was 13.8 GPa, significantly less than reported in the earlier studies. Though it could be argued that the data were not adjusted for the affects of non-parallel alignment (tow-in), it is unlikely that this could account for all of the discrepancy. Furthermore, in tensile testing performed by Sano et al. (1994), similarly low values of Youngs modulus were reported (13-15 GPa). Both of the above studies were performed with great attention to detail, so instrumentation artifacts or flaws in specimen preparation cannot readily explain the lower values of elastic modulus. However, one detail worth mentioning is that, in both studies, the specimens were stored for an undisclosed amount of time in water. In the work reported by Sano et al. (1994), the tensile specimens were stored in 0.9% NaCl water at 4°C for about 24 hrs prior to being tested (Pashley, 2001, personal communication). Earlier studies have showed that short-term storage in saline solution degrades bond strengths (Goodis et al., 1993), and that storage in water can reduce bend strength in bone (Gustafson et al., 1996). Given the small size of dentin specimens, it is possible that even short-term storage in water or saline might reduce the elastic modulus by dissolution of the mineral phase. The effects of water storage will be considered in greater detail later. In more recent work, Palamara et al. (2000) recorded the response of a grid pattern coated on the surfaces of dentin specimens during compressive loading. The results are of particular interest because, while the Youngs moduli measured in orientations both parallel and orthogonal to the tubule axes (the principal structural directions) were identical, the modulus at 45° off-axis to the tubules was determined to be lower. The observation that the off-axis modulus was smaller than that measured in either principal axis requires that the intertubular dentin matrix be anisotropic. Since this is a significant finding, we must consider its veracity in greater detail. Our analysis begins by a consideration of the simplest deviation from isotropic symmetry, cubic, which requires three independent elastic constants. The reciprocal Youngs modulus in a cubic system along the direction of the unit vector li can be expressed in terms of the compliance matrix Sij (Nye, 1972):
The variation of the Youngs modulus with orientation depends only on the terms in l, and is zero for the principal orthogonal directions and a maximum of 1/3 in the <111> directions. For the Youngs modulus to be less in the off-axis directions, (S11-S12-S44/2) must be less than zero. This requires the shear modulus, G, to be less than would be expected in an isotropic system:
With the values provided by Palamara et al. (2000), [E11 = 10.4 GPa and E(45°) = 7.7 GPa], we can calculate a rigorous upper bound for the shear modulus. For a cubic symmetry model, G < 2.9 GPa. This is less than would be required from isotropic symmetry for v = 0.25 (4.1 GPa). One can argue that, from an analysis of the dentin microstructure, an orthotropic symmetry would be a more reasonable alternative to cubic. However, we find that the situation is only slightly altered from the analysis for cubic symmetry. For orthotropic symmetry, we begin with the general case of a lamina loaded in plane at an arbitrary angle q with respect to a principal axis, l1, in this case parallel to the tubule axes. The Youngs modulus as a function of q can be expressed as (Jones, 1975)
In Eq. 6
Here, G12 is the shear modulus and v12 is Poissons ratio. For the data presented in Palamara et al. (2000), a = 1. The observation that Eq was less than both E1 and E2 places a rigorous restriction on the value of the shear modulus:
The above restriction on the shear modulus is identical to the one derived for cubic symmetry. With the values reported in Palamara et al. (2000) for the off-axis Youngs modulus, Eq. 6 Though the largest source of error in measurements with grids lies in the accurate determination of displacement (strain), the absence of any independent measure of displacement (e.g., no strain gauges or cross-head displacements were reported) prevents us from assessing this as a possible source of error. However, the compression tests were conducted between 50 and 150 MPa, which is within the range of the proportional limit in compression (80-140 MPa) that has been established in earlier work (Craig and Peyton, 1958; Stanford et al., 1960). It is therefore possible that Palamara et al. (2000) were operating beyond the yield stress when recording their grid displacements. This factor would explain the lower values of the Youngs modulus, and could easily explain the off-axis anisotropy. Also, considerable relaxation (creep) could have occurred during the time intervals required to obtain the images of the grids. However, in the absence of additional data, the possibility that dentin is elastically anisotropic cannot be ignored. The question of elastic symmetry will be taken up again in a later section.
Indentation measurements of Youngs modulus
The indentation modulus E*, sometimes referred to as the reduced modulus, is determined from the corrected contact stiffness, Sc, and the contact areas for each indentation, A, that are obtained from a tip shape calibration procedure (Doerner and Nix, 1986):
From Eqs. 9 and 10, it is possible to obtain the Youngs modulus of the probed specimen, Es:
In Eq. 11 Among the first to apply nanoindentation to the study of dentin and dental materials was van Meerbeek et al. (1993), who measured a Youngs modulus of dentin of 19.3 GPa. This was followed by Kinney et al. (1996), who used nanoindentation to measure the Youngs modulus of intertubular (18-21 GPa) and peritubular (29.8 GPa). Since then, nanoindentation has become a common technique for the determination of local mechanical properties of structural features in biological hard tissues (Rho et al., 1999). A significant drawback to the widespread use of nanoindentation in mineralized tissues is the inability to indent in water. To overcome this drawback, Balooch et al. (1998) applied the atomic force microscope (AFM) with a specially designed attachment called the Triboscope to perform indentations on fully hydrated specimens. In addition to allowing indentations to be made in water, the device also functioned as an AFM, allowing for precise positioning of the indenter and subsequent imaging of the indentation with microscopic spatial resolution. With this technique, Kinney et al. (1999) measured a pronounced decline in E with submersion in water. These investigators attributed this decrease to a combination of softening of the collagen phase and partial surface demineralization. Without knowing the effects of surface demineralization, they assumed that a plasticizing of the collagen fibrils and other noncollagenous proteins in the tissue caused the majority of the decrease in modulus with time in water.
Although nanoindentation probes only a thin surface layer (< 1 µm), the mechanical properties obtained are assumed to be representative of the bulk material. Chemical changes in the surface layer of mineralized tissues resulting from storage solutions are, thus, important considerations for accurate determination of nanomechanical properties. This aspect of nanoindentation has been recently explored by Habelitz et al. (2002a), who studied changes in nanomechanical properties of dentin and enamel during storage in de-ionized water, calcium-chloride-buffered saline solution, and Hanks Balanced Salt Solution (HBSS). The investigators were able to show that storing teeth in de-ionized water or CaCl2 solution resulted in a large decrease in elastic modulus and hardness. After one day of storage, a decrease in the Youngs modulus and hardness of up to 30% was observed. By one week of solution storage, mechanical properties dropped to below half of their initial values. In contrast, storing the specimens in HBSS did not significantly alter the mechanical properties for a time interval of two weeks. The behavior of the Youngs modulus of specimens stored in water and HBSS is reproduced in Fig. 6
The decrease in modulus and hardness of dentin when stored in de-ionized water or CaCl2 solution was attributed to partial demineralization of the near-surface layer. Because de-ionized water lacks calcium and phosphate ions, the chemical potential for dissolution of the mineral phase of dentin and enamel is high, and was assumed to be the major reason for the demineralization and softening of the tissues. Addition of calcium chloride to a saline solution also did not prevent demineralization. Since the pH of the CaCl2 solution was slightly acidic (pH = 5.9), it was more likely to dissolve the calcium phosphate minerals. HBSS, on the other hand, was slightly basic (pH = 8.0) and was unable to dissolve the mineral phases by acidic attack. It is highly concentrated in Ca2+, Mg2+, Na+, PO43-, and Cl- ions and has a composition comparable with that of the dental mineral phases. Therefore, the chemical potential of HBSS to dissolve the calcium phosphate phases in teeth is low, and surface demineralization was retarded. By eliminating near-surface demineralization through proper specimen storage procedures, one can now obtain consistent values of the Youngs modulus of hydrated peritubular and intertubular dentin. This is an important development, since it will allow for the study of the differences between wet and dry tissues, and provide elastic moduli more representative of in vivo conditions. The nanoindentation method is limited in that it cannot be used to determine any of the other elastic constants, and it does not provide a direct measure of the continuum Youngs modulus, since micromechanics arguments must be used to combine the separate measures of the inter- and peritubular dentin. Other techniques will now be described that have the potential to determine all of the second-order elastic constants of dentin completely.
Sonic measurements of the elastic constants
Using the technique of ultrasonic interferometry, Gilmore et al. (1969) established the sound speeds in bovine dentin. With these data, and assuming isotropic elasticity, the investigators used Eqs. 12 to derive the second-order elastic constants. To facilitate comparisons with other measurements, we have reproduced the entire range in values of the elastic moduli from this work (Table 1
The data from Gilmore et al. did not consider the possibility of anisotropy in the elastic constants. Lees and Rollins (1972) explored this possibility in a later study by measuring the longitudinal and shear wave velocities along orthogonal directions in the plane of the tubules. To facilitate these measurements, the investigators used the technique of critical angle reflection, where the angle of incidence of the sonic wave is scanned while the intensity of the reflected wave is monitored simultaneously. At certain critical angles, determined by the shear and longitudinal sound speeds in the specimen, the refraction vanishes, and there is a local maximum in the reflected signal. With Snells law, it is possible to calculate the longitudinal and shear wave speeds from the values of these critical angles.
The main difficulty with this method lies in determining the principal symmetry directions in the tooth from only two orientations. Lees and Rollins assumed hexagonal (transverse isotropic) symmetry based on the microstructure of dentin. However, since the elastic properties were measured in only a single plane, the solution to the stiffness tensor for hexagonal symmetry was underdetermined. Therefore, the investigators assumed that the lower-sound speed values reported earlier by Gilmore et al. were obtained from measurements 45° with respect to the tubule axes. With this assumption, the complete elastic stiffness matrix was estimated (Table 2
Lees and Rollins (1972) reported only the magnitudes of the elements in their derived stiffness tensor (Cij). By matrix inversion, we can obtain the elements of the strain tensor, from which it is a simple task to calculate the Youngs modulus. According to the hexagonal symmetry model and the stiffness tensor derived by Lees and Rollins, the Youngs modulus in the tubule direction was 36 GPa, and 29 GPa in the perpendicular axis.
The Youngs modulus, which can be derived from the stiffness matrix of Lees and Rollins (1972), is graphed in Fig. 7
In spite of the practicality of critical angle reflection methods, we must be skeptical of the conclusions. The chief criticism, one that was acknowledged by the investigators, was that the assumption that prior measurement made in different laboratories and using different methods corresponded to the local minimums in the hexagonal symmetry model. This assumption forces an anisotropy that may not be representative of the actual symmetry of the specimen. In fact, from our knowledge of bone, the modulus should be greatest in the direction of the collagen fibrils (Pidaparti et al., 1996). In dentin, this would mean that the elastic modulus should be greatest in the direction perpendicular to the tubules, the complete opposite of the critical angle results. Therefore, we must consider the possibility that Lees and Rollins (1972) made an incorrect assumption regarding the orientation of the specimens in the earlier work of Gilmore et al. Clearly, what is needed is an ultrasonic method suited to small specimens, and which allows for the determination of all of the elastic constants from a single orientation of the specimen and without assumptions as to the symmetry. Indeed, these requirements are met by a technique known as resonant ultrasound spectroscopy. Resonant ultrasound spectroscopy (RUS) is a method for measuring single-crystal elastic constants with great precision (Migliori et al., 1993). RUS makes use of Hookes law and Newtons second law to predict the resonant modes of vibration of a specimen of known shape (Migliori et al., 1993; Maynard, 1996). The entire stiffness tensor, Cij, can be determined by comparison of the frequency spectrum produced by the resulting eigenvalue problem with the measured resonant frequencies (vibrational eigenmodes) of the specimen. It is important to note a significant difference between RUS and other sonic methods. In RUS, sound speed is not measured. Rather, the resonant frequencies of mechanical vibration are determined. This means that the stiffness tensor can be determined for small specimens from one measurement at a single orientation. Since the development of RUS, the applications of this technique have been extended to include geological structures and complex particulate and fiber-reinforced composites (Jung et al., 1999). Recently, this technique has been applied to small specimens of human coronal dentin, with interesting results (Kinney et al., 2002). Cubes of dentin, approximately 2 mm on an edge, were mounted on opposing corners between two transducers in the RUS system. The resonant frequencies between 0.5 and 1.4 MHz were measured. The values of the resonant frequencies were calculated from an initial approximation of the stiffness tensor. The experimental and predicted frequency distributions were then compared, and the residuals, F, were calculated. The elastic constants were adjusted from their initial values to minimize these residuals. Isotropic, cubic, and hexagonal material symmetry groups were modeled.
The elastic constants of hydrated and dry dentin, as measured by RUS, are listed in Table 1
The measured anisotropy of the elastic constants of dentin can now be reconciled with its structural anisotropy. As in bone, the symmetry of the elastic constants is determined by the orientation of the collagen fibrils. However, the magnitude of the anisotropy (
Viscoelastic behavior of dentin
In a viscoelastic solid, the stress is out of phase with the strain. The amount of this phase shift is defined by a phase angle,
In a perfectly elastic solid, there would be no energy lost to creep deformation (E'' = 0), and the stress and strain would be in phase ( The classic analysis of viscoelasticity is based on the Maxwell element, a spring and dashpot arranged in series. The time-dependent relaxation modulus—for example, E(t)—can be written for a single element, i, as:
where Ei represents the stiffness of the spring element i, and
For an infinite number of Maxwell elements, the relaxation modulus can be written as a function of the continuous relaxation spectrum, H(
Similarly:
A few studies have explored stress relaxation in dentin (Korostoff et al., 1975; Tengrove et al., 1995). Still, little is known about its time-dependent behavior. In what appears to have been the most detailed study of viscoelasticity in dentin, Korostoff et al. (1975) measured stress relaxation in specimens prepared from the roots of human incisors, cuspids, and premolars. The cylindrical specimens were machined to dimensions approximately 6.4 mm tall, and had a 3.8-mm outer and a 1.1-mm inner diameter. The specimens were loaded in compression at 37°C to a constant strain,
In Eq. 18
Korostoff et al. (1975) observed an exponential decline in the relaxation modulus with time from the shortest time measured (t It is well-known that the indenter stylus continues to displace inward during the constant load portion of an indentation measurement of dentin. This continued displacement is evidence that creep relaxation is occurring; the standard procedure for analyzing the indenter load/displacement data does not account for this relaxation. Recently, Feng and Ngan (2002) derived a correction formula for creep during indentation measurements. The correction is applied to the contact stiffness, Sm, measured at the beginning of unload:
In Eq. 19
This information has never been reported for indentation measurements on calcified tissues. However, approximate values for the correction factor can be obtained by estimation from the graphed representations of load and displacement with time. A graphed representation of both the stylus displacement and the load as a function of time is shown in Fig. 8
It now appears as if creep relaxation can explain the small differences between AFM indentation and RUS measurements of the modulus. However, the creep mechanisms, and their possible dependence on the magnitude of the applied stress, remain unknown. An important question, given the absence of a well-defined time-dependent equation of state, is, what values of the elastic constants should be used to describe elastic deformation that occurs under physiological loading conditions? For example, what magnitudes of the elastic constants should be used in the frequency range of physiological interest (v = 1 Hz 1/t)? This question can be answered only partially at this time.
A complete analysis of the time-dependent equation of state for dentin requires knowledge of the relaxation spectrum, H(
With this approximation, Korostoff et al. (1975) evaluated the relaxation spectrum for root dentin for relaxation times from 30 sec to There were admitted shortcomings with the study. Most critical, no information was provided about stress relaxation at time scales shorter than 30 sec. Physiological time scales of greatest relevance lie in the range of 0.1 to 10 sec. There is a serious concern with extrapolating the above results to shorter time scales. In bone, deviation from simple exponential relaxation has been observed; stress relaxation at short time scales has followed a Kohlrausch-Williams-Watts (KWW) functional form (Sasaki et al., 1993). Much more research, particularly at short time scales, is warranted.
No one has measured the frequency dependence of the loss tangent (Eq. 13
SUMMARY OF ELASTIC PROPERTIES
Sonic methods provide the most direct and precise determination of the elastic constants. They are least affected by stress relaxation. Because of this, we propose that the data from the RUS (Kinney et al., 2002) and sound speed measurements (Gilmore et al., 1969) be used to define the most probable range for all of the elastic constants of dentin in the absence of viscoelastic effects (Table 1 The intrinsic elastic constants provide only part of the story. The viscoelastic behavior of dentin must be explored in greater detail if a constitutive equation of state is to be developed for dentin. In the frequency range of physiological interest (0.1-10 Hz), the phase angle between stress and strain is most likely small, so that dentin can be modeled as a perfectly elastic solid. However, it is likely that the effective elastic moduli are reduced by about 10-20% from the values obtained with sonic measurements. This is consistent with the differences between the RUS measurements and those obtained by indentation and uni-axial testing. However, this observation is still speculative: The mechanisms controlling viscoelastic behavior, and their possible dependence on the stress amplitude, remain unknown in dentin. As our knowledge of the elastic properties of normal dentin improves, other questions are raised that have significant implications for dentin pathologies. Is there a functional relationship between the mineral density and the elastic moduli? Does age-related transparency alter the elastic properties of root dentin? Does reparative dentin that forms at the dentin/pulp interface have the same elastic properties as the secondary dentin it adjoins? Answers to these questions will be important in the continued development of minimally invasive approaches in restorative dentistry.
Hardness tests measure the resistance of the dentin to deformation caused by penetration of an indenting stylus. Since its introduction, the Knoop (Knoop et al., 1939) indenter has proved to be the workhorse in studies of dentin. The long aspect ratio (7.11 times longer in one dimension) allows for accurate measurements of area, even with shallow indentations. Knoop, therefore, is extremely useful for the small, thin specimens typical of studies of mineralized tissues. Hardness is defined in units of pressure, or force per unit area of indentation. Unlike Vickers or Brinell methods, which use the contact area of the indenter stylus, the Knoop method uses the projected area, Ap, in the calculation of hardness (KNH):
where P is the applied load (in kg) and l is the length (in mm) across the long axis of the remnant impression. In SI units, 1 kg/mm2 Many investigators have determined that dentin hardness depends on mineral concentration. Featherstone et al. (1983) developed an analytic expression relating Knoop hardness to the volume percent of mineral (Vm):
Though the validity of Eq. 22
In addition to its association with the mineral concentration, hardness has been correlated with location in the tooth. The works of Ogawa (1983) and Wang and Weiner (1998b) show that the mantle dentin immediately subjacent to the enamel (KNH
Clearly, Eq. 22 It is tempting to try to relate dentin hardness to other physical properties like yield stress, tensile strength, or Youngs modulus. In ductile metals, for example, the yield stress and tensile strengths are often observed to scale with hardness. However enticing these scalings, we must remember that they were derived from plasticity theory for materials that display significant yielding. In contrast, mineralized tissues are more brittle, showing little if any yielding prior to failure.
Though we should not anticipate similar relationships with dentin hardness, there is a possibility that its hardness might scale with the Youngs modulus. Eq. 22
Failure data have largely been limited to measurements of shear strength or ultimate compressive or tensile strength. For dentin specimens failed in tension, reported magnitudes of the ultimate strength range from 52 MPa (Bowen and Rodriguez, 1962) to 105 MPa (Sano et al., 1994). It is likely that the two-fold differences in the tensile strength were related to flaws in the specimens, since both a size effect and an improvement in strength with surface sanding of the specimens have been noted (Sano et al., 1994). Moreover, measurements of compressive strength, which are less likely to be affected by flaws, were more consistent. Values of compressive strength range from 275 to 300 MPa (Craig and Peyton, 1958). Shear strength, measured either by punching or lap-shear, also has been highly variable. Using a shear punch apparatus, Cooper and Smith (1968) obtained values for shear strength that ranged from 64 to 132 MPa, not too dissimilar from later measurements by Roydhouse (1970) (69-147 MPa). More recent measurements by single-plane lap shear produced shear strengths of 36 MPa; this low value may have been due to the specimens having been from dentin closer to the pulp, or to problems with the experimental design, or with bending of the specimen (Gwinnett, 1994). Our group (Watanabe et al., 1996) demonstrated that some of the disparity in reported values of the shear strength could be attributed to tubule orientation and location within the tooth. The lowest values of lap shear strength were obtained in less mature dentin nearer the pulp, with the tubules oriented parallel to the shear plane (54 MPa). The highest values were obtained in mature cuspal dentin with the shear plane perpendicular to the tubules (92 MPa). Others also have observed a similar orientation dependence with the tensile strength (Lertchirakarn et al., 2001).
The large standard deviations common to all measurements of dentin strength suggest that strength is controlled by the flaw distribution in a specimen. A flaw is similar to the weakest link in a chain: Variations in the flaw size lead to variations in the failure strength. Specimens with large flaws will fail at lower stresses than will specimens containing smaller flaws. For a random distribution of flaws, the tensile strength should obey a Weibull probability distribution function (pdf), where the probability, Pf, that a specimen fails at a stress level
In Eq. 23
The Weibull modulus, m, determines the size-sensitivity of the failure strength; smaller values of m lead to greater sensitivity to specimen size.
Most studies of tensile strength have recorded only the means and standard deviations of the experimental measurements. An exception was the early work of Lehman (1967), who tabulated the test results of all 100 of his tensile strength measurements of root dentin. These results are recast in a Weibull distribution in Fig. 9
The previous discussion of dentin strength suggests that a fracture mechanics approach would be more appropriate than a strength-of-materials approach for the study of tooth failure. Yet, surprisingly few studies have taken this approach. Important exceptions to an otherwise-absence of quantitative fracture data are the works by Rasmussen et al. (Rasmussen et al., 1976; Rasmussen, 1984; Rasmussen and Patchin, 1984), El Mowafy and Watts (1986), and Imbeni et al. (2002). Rasmussen et al. determined the work-of-fracture of normal dentin in directions parallel and perpendicular to the tubule axis. This work was important in that it detected a small directional dependence in the work-of-fracture; it was energetically more favorable to fracture dentin perpendicular to the axis of the tubules than parallel to them. However, the relatively large works-of-fracture implied a toughening mechanism and the possibility of significant yielding ahead of the crack. Perhaps because of this, the authors stopped short of equating the work-of-fracture with any intrinsic material property such as fracture toughness.
The work by El Mowafy and Watts (1986) was of considerable importance. It was the first attempt to measure the intrinsic fracture toughness, Kc, of dentin with an ASTM standard specimen geometry: the compact tension specimen. For normal coronal dentin, these investigators reported a value of 3.08 MPa Three valuable findings came from the work of El Mowafy and Watts (1986): (1) Compact tension specimens could be fabricated from coronal dentin; (2) fracture toughness was high for a brittle material, indicating that either the collagen fibrils or the tubule lumens provide a toughening mechanism; and (3) fracture toughness was independent of temperature through a range from 0 to 60°C. However, since the fracture plane was parallel to the tubules (i.e., in the orientation associated with the highest work of fracture), it is possible that the fracture toughness could be considerably lower in the orthogonal plane. Furthermore, El Mowafy and Watts (1986) used a blunt notch instead of a fatigue precrack, so that notch toughening might have elevated their results. To explore these possibilities further, Imbeni et al. (2002) considered the effects of notch root radius on the measurement of fracture toughness in a direction perpendicular to the tubules.
Imbeni et al. (200) prepared 10-mm-long beam-shaped specimens from human molars. Tests were conducted in three-point bend, with a span between the loading points of approximately 5 mm. Two protocols were followed. For the first protocol, sharp notches were prepared in the specimens with the root radius ranging from 30 to 50 µm. For the second protocol, a fatigue precrack was grown out of the notch by fatigue cycling at 2 Hz with a maximum stress intensity of about 1 MPa Was the lower value of the fracture toughness measured in Imbenis work entirely the consequence of their having used a fatigue precrack, or was it also affected by the orientation of the fracture plane? This is important, since the critical flaw size, ac, is a strong function of the fracture toughness Kc:
If the difference in fracture toughness was an artifact of the notch radius, then for the same far-field applied stress,
Teeth are subject to cyclic loads during mastication. The frequency of the loading is nominally 1 Hz, and stress amplitudes of 20 MPa can be estimated at the cervical margins (Anderson, 1956). Therefore, it is important to know the response of the tooth to cyclic loading: its fatigue behavior. Characterizing fatigue behavior involves determining the total life to failure in terms of a cyclic stress range. This method is often referred to as the "S/N" (stress/life) approach, and requires the measurement of the number of cycles required to induce failure in a flaw-free material at a given alternating stress. Some materials exhibit a fatigue limit: a stress below which failure does not occur. Many materials, however, do not have a fatigue limit. For those cases, it is convenient to define a material lifetime in terms of the endurance strength, which is the alternating stress below which the material will not fail before a predefined number of cycles. Does dentin exhibit classic S/N behavior, and, if so, are the cyclic stresses of mastication below its endurance strength? A partial answer to this question was provided in a study of tensile fatigue behavior in bovine dentin (Tonami and Takahashi, 1997). In that study, tensile specimens with 1.5 x 1 x 1-mm-gauge sections were subjected to a sinusoidal cyclic unipolar load at 5 Hz. The load ratio, R, which is defined as the minimum load divided by the maximum load, was allowed to vary between 0.15 and 0.25. A staircase method was used to determine the stress at which the specimen failed at 105 cycles. For example, if a specimen did not fail at 105 cycles at a stress of 45 MPa, another specimen was tested at a stress level 5 MPa higher. With this procedure, Tonami and Takahashi (1997) established the endurance strength, for failure at 105 cycles, to lie between 45 and 50 MPa. This would suggest that fatigue is not a factor in tooth failure under conditions of normal masticatory loading. The work of Tonami and Takahashi was important in establishing the existence of fatigue failure in dentin. However, there were several shortcomings that limit its usefulness. First, the endurance strength was established for failure at 105 cycles. This was short, especially considering that 106 cycles are expected in a year from mastication. Second, the load frequency used in the study (5 Hz) was five times greater than that observed in physiological loading (1 Hz); the investigators assumed that there would be no frequency dependence on the fatigue life. This was a questionable assumption, however, especially considering the pronounced frequency dependence of fatigue in bone (Zioupos et al. 2001). Finally, the work was performed on bovine dentin. The recent work of Nalla et al. (2002) was designed to address many of these issues. In this study, specimens of human dentin were subjected to fatigue loading at a constant load ratio, R = 0.10, and at two different cyclic frequencies: 2 and 20 Hz. The authors observed the existence of a classic S/N fatigue response to repetitive loading, and identified an apparent fatigue limit at 106-107 cycles, which was estimated to lie between 25 and 45 MPa at frequencies of 2 and 20 Hz, respectively. Failure appeared to be caused by the initiation and growth of a single, dominant crack. In the absence of intrinsic flaws, it appears that fatigue failure will not occur in human teeth. However, the Weibull nature of tensile failure indicates that there is an inherent population of pre-existing flaws in normal dentin. Thus, the existence of a fatigue limit in a nominally flaw-free specimen of dentin cannot be generalized to the whole tooth. The classic S/N approach is not conservative enough; it does not address the question of whether a pre-existing flaw can grow to a critical size during normal loading conditions. This requires the measurement of fatigue growth crack rates, da/dN, for a fracture-mechanics-based approach to failure of human dentin.
Nalla et al. (2002) determined the fatigue crack growth rate by measuring stiffness loss in the specimen during cyclic loading. They observed a Paris power law relationship of the crack growth rate with the stress intensity,
The Paris law exponent, m = 8.76, was smaller than that of carbonated apatite bone substitute (m = 17; Morgan et al. 1997), but larger than that frequently observed in ductile metals (m
When a crack reaches a critical size, ac, determined by equating the fracture toughness, Kc, with the stress intensity at the crack tip, the structure will fail. Nalla et al. (2002) used lifetime modeling to estimate the number of cycles, Nf, required for an incipient crack of length ao to grow to the critical size under a far-field applied stress
In Eq. 27
Eq. 27
Two major assumptions have gone into the development of Eq. 27
Many of our conceptions of basic biomechanical properties of dentin have changed since the last comprehensive review (Waters, 1980). The magnitudes of the elastic constants must be revised upward, and viscoelastic effects must be taken into account. Also, the concept of strength as an important engineering quantity must be replaced by a fracture mechanics approach to dentin failure, since pre-existing flaws can cause teeth to fail at stresses far less than their theoretical strength. In particular, a lifetime modeling, or damage-tolerant, approach appears promising for the development of a clinical predictor of tooth failure. These paradigm shifts have been facilitated by advances in measurement science combined with a better understanding of dentin microstructure. It now appears that the elastic constants of dentin have hexagonal symmetry, with the stiffest direction oriented in the plane perpendicular to the tubules. This finding is consistent with a micromechanics model that suggests that the peritubular dentin has negligible influence on material symmetry. Furthermore, it is also consistent with the observations in other tissues that the stiffest orientation is in the direction of the mineralized collagen fibrils. We can conclude, then, that the elastic properties of dentin are explainable in terms of the microstructure of the intertubular dentin matrix, and that any correlation of the elastic properties with the tubule direction is a necessary consequence of the orthogonal relationship between the tubules and collagen fibrils. Therefore, research must focus on the mineralized collagen scaffold, and explore the coupling between the mineral and collagen phases and its effects on biomechanical properties. The Weibull behavior of tensile strength data suggests that failure initiates at flaws. These flaws may be intrinsic, perhaps regions of altered mineralization, or extrinsic, caused by cavity preparation, wear, or damage. Clearly, identifying the flaws responsible for failure is important, as is the understanding of how sub-critical flaws can grow and coalesce with cyclic loading. Though the cyclic forces of mastication are insufficient to cause failure in perfect dentin, it now appears that pre-existing flaws in normal dentin can grow to catastrophic size with application of small cyclic loads. Furthermore, understanding how shifts in the chemical environment affect the crack growth rate will be essential for developing lifetime models. At present, many of these data are lacking. Finally, little is known about the biomechanical properties of altered forms of dentin. Studies seem to suggest that there are at least two forms of transparent dentin, a form associated with caries and a form associated with age-related changes in the root. Each type of transparency may exhibit unique biomechanical signatures. In addition, there are caries lesions and genetic disorders, such as dentinogenesis imperfecta, which appears to lack the phosphoproteins that promote binding between the mineral and the collagen. It is absolutely essential that the properties of these altered forms of dentin be obtained before an understanding of how dentin pathologies affect tooth strength can be developed. Current conservative dental practice is increasingly oriented toward efforts to arrest and remineralize caries lesions with either solution chemistry or more fundamental tissue engineering approaches. An accurate understanding of the properties of the dental hard tissues is a requirement for assessment of the effectiveness of these approaches. Without knowledge of mineral-collagen coupling and its biomechanical consequences, we may not truly be restoring the affected tissues. This shift in emphasis from traditional filling procedures to tissue restoration requires a firm understanding of the structure/properties relationships in dentin at many length scales. It is hoped that this review has clarified what is known, and has raised awareness of what remains to be learned about the mechanical properties of dentin.
This work was supported in part by the National Institutes of Health, National Institute of Dental and Craniofacial Research (grants PO1 DE09859, DE 11526, and DE13029). The authors thank many investigators for sharing their time, comments, and data prior to publication. In particular, we thank Prof. David Pashley, Prof. Lawrence Katz, Prof. Robert O. Ritchie, Dr. James S. Stolken, Dr. Stefan Habelitz, and R.K. Nalla for valuable discussions.
Critical Reviews in Oral Biology & Medicine, Vol. 14, No. 1,
13-29 (2003) This article has been cited by other articles:
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
30 MPa, a flaw-free dentin specimen apparently will not fail. However, a more conservative approach based on fatigue crack growth rates indicates that if there is a pre-existing flaw of sufficient size (







is a "clamping" factor that accounts for mismatch in Poissons ratio between the peritubular and intertubular dentins, and Al and Apt are the area fractions of the intertubular and peritubular dentin, respectively. The area ratio in Eq. 3











, by (



, between the applied stress and the resulting strain. A detailed analysis defines the loss tangent of the phase angle in terms of a storage modulus, E', and a loss modulus, E'': 

i is the relaxation time of the element. A parallel arrangement of a group of Maxwell elements can be treated as a simple sum: 



c, of 0.6%. The relaxation modulus, Er(t), was determined from the experimentally observed drop in applied load, L(t), and the change in specimen height,
hc: 





is given by the expression: 


m with a 10% standard deviation. Based upon an unpublished value of the yield stress, the authors assumed that the test resulted in a valid measure of Kc. They could therefore equate their fracture toughness measures with the energy release rate, G, and obtained good agreement with the earlier data of Rasmussen (







