Determination of elastic and plastic mechanical properties of dentin based on experimental and numerical studies

The aim of this study is to investigate the change of mechanical properties of human dentin due to aging and spatial variation. Sections of coronal dentin are made from human molars in three groups: young, mid-aged, and old patients. A nanoindentation test is conducted from regions near the pulp to the dentin-enamel junction (DEJ) to evaluate the load-depth indentation response and determine Young's modulus and hardness. Based on the loading and unloading load-displacement curves in nanoindentation, a numerical model of plastic damage is used to study the plastic and the damage behaviors and the contribution to the degradation in the unloading stiffness. The experimental results show that Young's modulus of the inner dentin is significantly lower than that of outer dentin in each age group. Compared with the young dentin, the old dentin has greater hardness and Young's modulus with similar spatial variations. The magnitudes of the yield strength and the damage variable are also affected by aging and vary with spatial locations. In the same age group, the yield strength in inner dentin is lower than those in middle and outer dentin, more damage occurs with similar spatial variations, and the yield strength of young dentin is generally lower and causes more damage compared with those in both the mid-aged and old groups.     

Free convection from one thermal and solute source in a confined porous medium

This paper reports a numerical study of double diffusive natural convection in a vertical porous enclosure with localized heating and salting from one side. The physical model for the momentum conservation equation makes use of the Darcy equation, and the set of coupled equations is solved using the finite-volume methodology together with the deferred central difference scheme. An extensive series of numerical simulations is conducted in the range of −10 ⩽ N ⩽ + 10, 0 ⩽ R _t ⩽ 200, 10⁻² ⩽ Le ⩽ 200, and 0.125 ⩽ L ⩽ 0.875, where N, R _t , Le, and L are the buoyancy ratio, Darcy-modified thermal Rayleigh number, Lewis number, and the segment location. Streamlines, heatlines, masslines, isotherms, and iso-concentrations are produced for several segment locations to illustrate the flow structure transition from solutal-dominated opposing to thermal dominated and solutal-dominated aiding flows, respectively. The segment location combining with thermal Rayleigh number and Lewis number is found to influence the buoyancy ratio at which flow transition and flow reversal occurs. The computed average Nusselt and Sherwood numbers provide guidance for locating the heating and salting segment.     

Numerical calculations of supersonic underexpanded jets in a cocurrent flow with the use of parabolized Navier-Stokes equations

Results of a numerical study of three-dimensional supersonic jets propagating in a cocurrent flow are described. Averaged parabolized Navier-Stokes equations are solved numerically on the basis of a developed scheme, which allows calculations in supersonic and subsonic flow regions to be performed in a single manner. A jet flow with a cocurrent flow Mach number 0.05 ⩽ M_∞ ⩽ 7.00 is studied, and its effect on the structure of the mixing layer is demonstrated. The calculated results are compared with available experimental and numerical data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 54–63, May–June, 2008.     

Viscoelastic properties of dibenzylidene sorbitol (DBS) physical gels at high frequencies

The rheological behavior of Dibenzylidene D-Sorbitol (DBS) gels formed in ethylene glycol, glycerol, mineral oil, ethanol, and chlorobenzene was studied using oscillatory squeezing flow viscometry. The storage (G ^′) and loss (G ^″) moduli were measured as a function of gellant concentration (0.5–2 w/w) and type of solvent. As expected greater values of gel strength (G ^′) were observed for gels containing higher concentrations of DBS. In addition, both storage and loss moduli of 2% systems were mostly frequency independent over the studied range, whereas 0.5% gels did exhibit some degree of dependence. We also found that the solvent plays an important role in the properties of the gels. Among the parameters that affect the viscoelastic properties of DBS gels, the solvent polarity and its ability to form hydrogen bonding may have significant effects on the gel rheology.     

Experimental determination of the friction factor coefficient in microchannels

Experimental data used to determine the friction factor in microchannels are presented. Glass microchannels with a circular cross section, diameters of 34.5, 33.6, and 24.5 μm, and different lengths are studied in the experiments. Pure deionized degassed water is used as a working fluid. The range of the Reynolds numbers is 13 ⩽ Re ⩽ 330. A method of two channels is used to calculate the friction factor. The results obtained are in good agreement with theoretical predictions for the case of a developed laminar flow in circular channels, but the use of the method of two channels has some constraints.     

Methods for Examining the Fatigue and Fracture Behavior of Hard Tissues

An understanding of the fatigue and fracture behavior of hard tissues (e.g., bone and tissues of the human tooth) is critical to the maintenance of physical and oral health. Recent studies suggest that there are a number of mechanisms contributing to crack extension and crack arrest in these materials, and that they appear to be a function of moisture and age of the tissue. An understanding of these processes can provide new ideas that are relevant to the design of multi-functional engineering materials. As a result, we have adopted the use of microscopic Digital Image Correlation (DIC) to examine the mechanisms of crack growth resistance and near-tip displacement distribution for cracks in human dentin that are subjected to opening mode loads. We have also developed a special compact tension (CT) specimen that permits evaluation of crack extension within small portions of tissue under both quasi-static and fatigue loads. The specimen embodies a selected portion of hard tissue within a resin composite restorative and enables an examination of diseased tissue, or portion with specific physiology, that would otherwise be impossible to evaluate. In this paper we describe application of these experimental methods and present some recent results concerning fatigue crack growth and stable crack extension in dentin and across the dentin-enamel-junction (DEJ) of human teeth.

A fractal model for thermal shrinkage of oriented polymers

A relationship between the amount of thermal shrinkage and the structural parameters of polyethylene extruded in the solid state and polymerization-filled compositions based on it is studied. It is shown that thermal shrinkage is determined by the fractal dimension of the polymer structure and the degree of stretching of the molecular chain. Galkin Donetsk Physicotechnical Institute, National Academy of Sciences of Ukraine, Donetsk 340114. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 1, pp. 160–163, January–February, 1998.     

Behaviour of the Reynolds stress on rough walls

Quadrant analysis is used to study the contributions, associated to the four quadrants of the (u,v) plane, to the production of the turbulent shear stress on rough walls. The measurements are described for a fully developed turbulent flow between two rough plates with varying the parameterλ _z(span/height ratio of roughness elements). The application of this technique indicates that the (Q ₂) events (ejections) and (Q ₄) events (sweeps) cause an intense production of the Reynolds stress —ρ as compared with the (Q ₁) and (Q ₃) quadrants. The (Q ₂) contribution to the Reynolds stress depends on the geometry factorλ _z. Variation of the parameterλ _zaffects the distributions of and forH⩽3, whereH is a particular threshold. Comparison with boundary layer flow shows that the region 0.2⩽y/h⩽0.7 is characterized by a Reynolds stress production, independent of the flow nature. The third moment of the longitudinal velocity fluctuations is found to be sensitive to the surface roughness.     

Effect of transverse shears on complex nonlinear vibrations of elastic beams

Models of geometrically nonlinear Euler-Bernoulli, Timoshenko, and Sheremet’ev-Pelekh beams under alternating transverse loading were constructed using the variational principle and the hypothesis method. The obtained differential equation systems were analyzed based on nonlinear dynamics and the qualitative theory of differential equations with using the finite difference method with the approximation O(h²) and the Bubnov-Galerkin finite element method. It is shown that for a relative thickness λ ⩽ 50, accounting for the rotation and bending of the beam normal leads to a significant change in the beam vibration modes.     

Investigation of shrinkage-induced stresses in stereolithography photo-curable resins

An experimental study has been undertaken to investigate the shrinkage characteristics of acrylic-based and epoxy-based stereolithography (SL) photopolymer resin systems after they have been laser cured and post-cured under ultraviolet (UV), and thermal exposure. The induced residual stresses and strains were determined by the shadow moiré and the hole-drilling strain-gage methods. Out-of-plane displacements (warpage) of acrylic-based post-cured resin plates were recorded by means of the shadow moiré method and correlated to the shrinkage strains by theoretical analysis. The induced residual stresses in the epoxy-based cylindrical resin specimens were determined from strains of three-element strain-gage rosettes of the blind-hole drilling method. Results are presented for the shrinkage stresses and strains for both material systems as a function of the post-curing process (UV, thermal). It was found that the shrinkage strains in the acrylic-based photopolymer resin were of considerable magnitude, while thermal post-curing resulted in higher shrinkage stresses for both material systems. The values of the shrinkage stresses compare well with those of the existing literature.     

Spatial local solutions of the Navier-Stokes equations

This paper considers solutions of the Navier-Stokes equations polynomial in the coordinates, which. are called local solutions. For an incompressible fluid, all higher-order terms (sums of higher-order. monomials) of degree 2 are found and it is proved that nontrivial axisymmetric higher-order terms. of degree higher than 2 do not exist. Nonsolenoidal axisymmetric solutions are listed, which can be. treated as steady-state barotropic gas flows in a potential external-force field. All elliptic vortices. generalizing the well-known Kirchhoff solution are calculated. All solutions of degree 3 with the. higher-order term of partial form are found. Some of these solutions break down in a finite time. regardless of the value and sign of viscosity. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 2, pp. 109–119, March–April, 2009.     

Pressure Filtration of Australian Bagasse Pulp

A one-dimensional pressure filtration model that can be used to predict the behaviour of bagasse pulp has been developed and verified in this study. The dynamic filtration model uses steady state compressibility parameters determined experimentally by uniaxial loading. The compressibility parameters M and N for depithed bagasse pulp were determined to be in the ranges 3000–8000 kPa and 2.5–3.0 units, respectively. The model also incorporates experimentally determined steady state permeability data from separate experiments to predict the pulp concentration and fibre pressure throughout a pulp mat during dynamic filtration. Under steady state conditions, a variable Kozeny factor required different values for the permeability parameters when compared to a constant Kozeny factor. The specific surface area was 25–30% lower and the swelling factor was 20–25% higher when a variable Kozeny factor was used. Excellent agreement between experimental data and the dynamic filtration model was achieved when a variable Kozeny factor was used.     

Some possible alternative solutions of near-tip fields for steady crack growth in elastic perfectly-plastic media

Some possible alternative solutions of near-tip fields are studied for plane-strain Mode—I quasi-static steady crack growth in incompressible (ν=1/2) elastic perfectly-plastic media. A group of four-sector solutions and a three-sector solution in which the elastic-unloading region vanishes are given. Stress functions, plastic flow factors and plastic strains in each region are also given. Project supported by the State Education Commission under a funding program for Excellent University Young Faculties and National Natural Science Foundation of China.     

Existence and time-discretization for the finite-strain Souza–Auricchio constitutive model for shape-memory alloys

We prove the global existence of solutions for a shape-memory alloys constitutive model at finite strains. The model has been presented in Evangelista et al. (Int J Numer Methods Eng 81(6):761–785, 2010) and corresponds to a suitable finite-strain version of the celebrated Souza–Auricchio model for SMAs (Auricchio and Petrini in Int J Numer Methods Eng 55:1255–1284, 2002; Souza et al. in J Mech A Solids 17:789–806, 1998). We reformulate the model in purely variational fashion under the form of a rate-independent process. Existence of suitably weak (energetic) solutions to the model is obtained by passing to the limit within a constructive time-discretization procedure.     

The effects of age on viscoelasticity and deformation recovery of articular chondrocytes

本文旨在表征年龄对软骨细胞在自然生长过程中的黏弹性和恢复变形能力的影响. 结果表明：年龄对软骨细胞黏弹性及其恢复变形能力产生显著影响，老年组软骨细胞各 项黏弹性参数值均明显高于幼年和中年组软骨细胞(p<0.0001), 而后两组无显著差异 (p>0.05); 老年组软骨细胞蠕变达到平衡态所需时间t_E显著小于幼年和中年组 (p<0.05), 而幼年和中年组无显著差异(p>0.05). 老年组软骨细胞最大蠕变位移 L_M显著大于幼年和中年组(p<0.005), 而幼年和中年组无显著差异(p>0.05). 老年组软骨细胞恢复变形时间t_R显著大于幼年和中年组(p<0.005), 而幼年和中年组 无显著差异(p>0.05). 恢复变形前8s的分析发现，幼年组软骨细胞恢复变形率K_y 显著高于中年和老年组(p<0.005), 而中年组(K_a)和老年组(K_o)软骨细胞 的恢复变形率无显著差异(p>0.05); 此外, 实验发现老年组软骨细胞的残余变形L_R 比幼年和中年组显著增大(p<0.005), 而后两组无显著差异(p>0.05). 研究工作对 于软骨组织工程、软骨细胞与支架材料相互作用以及探讨OA发生过程中的力学生物学机制具 有理论意义.     

Behavior of Piezoelectric Beams under Axially Non-uniform Distributed Loading

The paper presents a generic solution methodology for a quasi-static homogeneous monoclinic piezoelectric beam under axially distributed electric and mechanical surface loads and body forces expressed as polynomials of degree K≥ 0 of the axis variable. (In the absence of any electrical loading, this problem is known as the Almansi–Michell problem). The stress and the electrical displacement components are presented as a set of polynomials of degree ≤K+2 of the axis variable (“solution hypothesis”) containing 4K unknown tip loading constants and 3K stress functions of two variables. The cases K=0,1 stand for uniform or linear distributed loads in the axis direction. The analysis is initiated by the Kth level and continues down to lower levels. The main result of this work generalizes the “elastic” solution given recently by O. Rand and the first author (2005). Examples of solutions for axially uniform distributed loads (K=0), and equilibrium in which the stress and the electrical displacement do not depend on the axis variable, are presented. The applications to constant body loads and a hydrostatic pressure are considered.      

Radial expansion of a granular medium in spherical and cylindrical layers

An exact solution that describes the fields of displacements and stresses in an expanding spherical layer is constructed within the framework of the theory of small strains of a granular medium with rigid particles. For finite strains, the problem reduces to a nonlinear system of ordinary differential equations, which is solved by numerical methods. Similar solutions are found in the problem for a cylindrical layer. Based on these solutions, the effect of the dilatancy of the granular medium on the stress-strain state near expanding cavities is found. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 3, pp. 190–196, May–June, 2009.     

Antiplane strain in a nonlinearly elastic incompressible body

The stress-strain state of an incompressible cylindrical elastic body with antiplane strain under the action of potential forces and surface loading constant along the body is considered in a nonlinear formulation in actual variables. The stresses are expressed via the pressure and independent strains, the pressure is expressed via the force and elastic potentials, and nonlinear boundary-value problems are posed for strains (and displacements). Various methods for solving these problems are developed. For the nonlinear equations obtained, some analytical solutions containing free parameters are given, which can be used as a basis for solving particular problems. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 6, pp. 93–101, November–December, 2006.     

Studies of quality of geometrically nonlinear elasticity theory for small strains and arbitrary displacements

Earlier it was shown in [1, 2] that the equations of classical nonlinear elasticity constructed for the case of small strains and arbitrary displacements are ill posed, because their use in specific problems may result in the appearance of “spurious” bifurcation points. A detailed analysis of these equations and the construction, in their stead, of consistent equations of geometrically nonlinear theory of elasticity can be found in [3]. Certain steps in this direction were also made in [4, 5]. In [3], it was also stated that the methods and applied program packages (APPs) based on the use of the classical relations of nonlinear elasticity require some revision and correction. In the present paper, this conclusion is justified and confirmed by numerical finite-element solutions of several three-dimensional geometrically nonlinear deformation problems and linearized problems on the stability of equilibrium of a rectilinear beam. These solutions were obtained by using two APPs developed by the authors and the well-known APP “ANSYS.” It is shown that the classical equations of the geometrically nonlinear theory of elasticity, which underly the first of the developed APP and the well-known APP “ANSYS,” often lead to overestimated buckling loads for structural members as compared with the consistent equations proposed in [1–3].     

On the measurement of elastoplastic stresses

In principle, one should be able to measure elastoplastic stresses in the same manner as one does elastic stresses; i.e., measure the strains and compute the stresses from the constitutive law. In practice, this is rarely done because of the more complicated material response and the anisotropy of the plastic behavior. Further, elastoplastic stresses should be computed incrementally in the general case. This paper presents procedures for computing stresses from elastoplastic strains measured incrementally in a test under microcomputer control. The approach is evaluated for four different materials—two obeying the assumptions of classical plasticity and two showing anisotropic behavior—by computing the stresses in a smooth specimen from measured principal strains. A useful application is presented by computing the stresses at a notch root from biaxial strains measured with laser-based interferometry. The general conclusion is that even in situations where the material is clearly anisotropic, this approach can give a reasonableestimate of the largest local principal stress. Paper was presented at the 1991 SEM Spring Conference on Experimental Mechanics held in Milwaukee, WI on June 9–13.