On Asymptotics of Solutions of Nonlinear Differential Equations of the Second Order

We study the problem of asymptotics of unbounded solutions of differential equations of the form y″ = α₀ p(t)ϕ(y), where α₀ ∈ {−1, 1}, p: [a, ω[→]0, +∞[, −∞ < a < ω ≤ +∞, is a continuous function, and ϕ: [y ₀, +∞[→]0, +∞[ is a twice continuously differentiable function close to a power function in a certain sense.__________Translated from Neliniini Kolyvannya, Vol. 8, No. 1, pp. 18–28, January–March, 2005.     

A field model interpretation of crack initiation and growth behavior in ferroelectric ceramics: change of poling direction and boundary condition

Strain energy density expressions are obtained from a field model that can qualitatively exhibit how the electrical and mechanical disturbances would affect the crack growth behavior in ferroelectric ceramics. Simplification is achieved by considering only three material constants to account for elastic, piezoelectric and dielectric effects. Cross interaction of electric field (or displacement) with mechanical stress (or strain) is identified with the piezoelectric effect; it occurs only when the pole is aligned normal to the crack. Switching of the pole axis by 90° and 180° is examined for possible connection with domain switching. Opposing crack growth behavior can be obtained when the specification of mechanical stress σ^∞ and electric field E^∞ or (σ^∞,E^∞) is replaced by strain ε^∞ and electric displacement D^∞ or (ε^∞,D^∞). Mixed conditions (σ^∞,D^∞) and (ε^∞,E^∞) are also considered. In general, crack growth is found to be larger when compared to that without the application of electric disturbances. This includes both the electric field and displacement. For the eight possible boundary conditions, crack growth retardation is identified only with (E_y^∞,σ_y^∞) for negative E_y^∞ and (D_y^∞,ε_y^∞) for positive D_y^∞ while the mechanical conditions σ_y^∞ or ε_y^∞ are not changed. Suitable combinations of the elastic, piezoelectric and dielectric material constants could also be made to suppress crack growth.     

Crack size and speed interaction characteristics at micro-, meso- and macro-scale

The motivation to examine physical events at even smaller size scale arises from the development of use-specific materials where information transfer from one micro- or macro-element to another could be pre-assigned. There is the growing belief that the cumulated macroscopic experiences could be related to those at the lower size scales. Otherwise, there serves little purpose to examine material behavior at the different scale levels. Size scale, however, is intimately associated with time, not to mention temperature. As the size and time scales are shifted, different physical events may be identified. Dislocations with the movements of atoms, shear and rotation of clusters of molecules with inhomogeneity of polycrystals; and yielding/fracture with bulk properties of continuum specimens. Piecemeal results at the different scale levels are vulnerable to the possibility that they may be incompatible. The attention should therefore be focused on a single formulation that has the characteristics of multiscaling in size and time. The fact that the task may be overwhelmingly difficult cannot be used as an excuse for ignoring the fundamental aspects of the problem.Local nonlinearity is smeared into a small zone ahead of the crack. A “restrain stress” is introduced to also account for cracking at the meso-scale.The major emphasis is placed on developing a model that could exhibit the evolution characteristics of change in cracking behavior due to size and speed. Material inhomogeneity is assumed to favor self-similar crack growth although this may not always be the case. For relatively high restrain stress, the possible nucleation of micro-, meso- and macro-crack can be distinguished near the crack tip region. This distinction quickly disappears after a small distance after which scaling is no longer possible. This character prevails for Mode I and II cracking at different speeds. Special efforts are made to confine discussions within the framework of assumed conditions. To be kept in mind are the words of Isaac Newton in the Fourth Regula Philosophandi:

“Men are often led into error by the love of simplicity which disposes us to reduce things to few principles, and to conceive a greater simplicity in nature than there really is … We may learn something of the way in which nature operates from fact and observation; but if we conclude that it operates in such a manner, only because to our understanding that operates to be the best and simplest manner, we shall always go wrong.”––Isaac Newton

Effect of material nonhomogeneities on the HRR dominance

This work studies the asymptotic stress and displacement fields near the tip of a stationary crack in an elastic–plastic nonhomogeneous material with the emphasis on the effect of material nonhomogeneities on the dominance of the crack tip field. While the HRR singular field still prevails near the crack tip if the material properties are continuous and piecewise continuously differentiable, a simple asymptotic analysis shows that the size of the HRR dominance zone decreases with increasing magnitude of material property gradients. The HRR field dominates at points that satisfy |α⁻¹ ∂α/∂x_δ|1/r, |α⁻¹ ∂²α/(∂x_δ ∂x_γ)|1/r², |n⁻¹ ∂n/∂x_δ|1/[r|ln(r/A)|] and |n⁻¹ ∂²n/(∂x_δ ∂x_γ)|1/[r²|ln(r/A)|], in addition to other general requirements for asymptotic solutions, where α is a material property in the Ramberg–Osgood model, n is the strain hardening exponent, r is the distance from the crack tip, x_δ are Cartesian coordinates, and A is a length parameter. For linear hardening materials, the crack tip field dominates at points that satisfy |E_tan⁻¹ ∂E_tan/∂x_δ|1/r, |E_tan⁻¹ ∂²E_tan/(∂x_δ ∂x_γ)|1/r², |E⁻¹ ∂E/∂x_δ|1/r, and |E⁻¹ ∂²E/(∂x_δ ∂x_γ)|1/r², where E_tan is the tangent modulus and E is Young’s modulus.     

Influence of singularity dominated zone for propagating cracks in finite size specimens

The singularity dominated zones for straight as well as curved cracks propagating in finite size specimens were determined experimentally by using the optical method of dynamic photoelasticity using the near-field stress equations. Experimental data was carefully analyzed using improved numerical schemes to get the complete stress field around the propagating crack. This stress field was critically examined to evaluate the size of singularity dominated zones for cracks propagating in straight as well as curved paths. For this purpose, the exact solution was compared with the singular solution using stress components σ_x, σ_y, τ_xy and the maximum shear stress τ_max as a criterion respectively. For straight cracks where the stress field is symmetric about the crack path, the singularity dominated zones can be determined by using any one of the stresses. However, for a curved crack, the zones were unsymmetric. This study shows that σ_y, the crack opening stress, yields the best result for characterizing the singularity dominated zone around a running crack tip.     

Influence of stress triaxiality and temperature on void growth and ductile/brittle transition behavior of 40 Cr steel

Examined experimentally are the influence of stress triaxiality and temperature on the growth of microvoids and the ductile/brittle transition (DBT) macrobehavior of 40 Cr steel subjected to two different heat treatments. This is accomplished by testing more than 300 smooth and notched specimens over a temperature range of 20°C to −196°C. Changes in the microstructure morphology are examined by scanning electron microscopy (SEM) and identified with fracture data on a surface constructed from the uniaxial strain ε_c at fracture, the stress triaxiality R_σ and the temperature T. While stress triaxiality has a significant influence on the DBT temperature T_c, it does not affect the ratio of the average radius of voids R_o to that of inclusions R_i. The ratio R_o/R_i is found to increase with temperature and remains constant in specimens with different notch radii regardless of the temperature. Empirical relations between T_c and R_σ⁻ and R_o/R_i and T are proposed to better understand how macrofracture parameters are influenced by microstructure entities.     

Effect of grain size on slow fatigue crack propagation and plastic deformation near crack tip

Fatigue crack growth and its threshold are investigated at a stress ratio of 0.5 for the three-point bend specimen made of Austenitic stainless steel. The effect of grain size on the crack tip plastic deformation is investigated. The results show that the threshold value Δk_th increases linearly with the square root of grain size d and the growth rate is slower for materials with larger grain size. The plastic zone size and ratio for different grain sizes are different at the threshold. The maximum stress intensity factor is k_max and σ_ys is the yield strength. At the same time, the characteristics of the plastic deformation development is discontinuous and anti-symmetric as the growth rate is increased from 2·10^—8 to 10⁻⁷ mm/cycle.A dimensionless relation of the form for collating fatigue crack starting growth data is proposed in which Δk_th represents the stress intensity factor range at the threshold. Based on experimental results, this relation attains the value of 0.6 for a fatigue crack to start growth in the Austenitic stainless steel investigated in this work. Metallurgical examinations were also carried out to show a transgranular shear mode of cyclic cleavage and plastic shear.     

On perfectly plastic flow in porous material

Experimental observations suggest that for perfectly-plastic materials containing pores, the (small) strain at which significant macroscopic yielding occurs is relatively insensitive to porosity, for volume fractions below approximately 15–20% (although the yield stress drops significantly with increasing porosity). Another observation is that, at these porosity levels, the stress–strain curve remains approximately linear almost up to the yield point. Based on these observations, Sevostianov and Kachanov constructed yield surfaces that explicitly reflect the shapes of the pores and their orientation. The underlying microscale mechanism is that local plastic “pockets” near pores blunt the stress concentrations; as a result, they remain limited in size and well contained in the elastic field until they connect and almost the entire matrix plasticizes within a narrow interval of stresses that can be idealized as the yield point. The present paper provides direct insight into the micromechanics of poroplasticity through direct microscale numerical simulation. Besides confirming the basic microscale mechanism, these simulations reveal that the reduction of the macroscopic poroplastic yield stress is approximated quite closely by 1−v₂ times the dense nonporous yield stress, where v₂ is the volume fraction of the pores.     

Effects of deformed volume, volume fraction and particle size on the deformation behaviour of W/Cu composites

Tungsten/copper (W/Cu) particle reinforced composites were used to investigate the scaling effects on the deformation and fracture behaviour. The effects of the volume fraction and the particle size of the reinforcement (tungsten particles) were studied. W/Cu-80/20, 70/30 and 60/40 wt.% each with tungsten particle size of 10 μm and 30 μm were tested under compression and shear loading. Cylindrical compression specimens with different volumes (D_S = H) were investigated with strain rates between 0.001 s⁻¹ and about 5750 s⁻¹ at temperatures from 20 °C to 800 °C. Axis-symmetric hat-shaped shear specimens with different shear zone widths were examined at different strain rates as well. A clear dependence of the flow stress on the deformed volume and the particle size was found under compression and shear loading. Metallographic investigation was carried out to show a relation between the deformation of the tungsten particles and the global deformation of the specimens. The size of the deformed zone under either compression or shear loading has shown a clear size effect on the fracture of the hat-shaped specimens.The quasi-static flow curves were described with the material law from Swift. The parameters of the material law were presented as a function of the temperature and the specimen size. The mechanical behaviour of the composite materials were numerically computed for an idealized axis-symmetric hat-shaped specimen to verify the determined material law.     

Plane-strain crack tip field and a dual-parameter fracture-criterion for non-linear fracture mechanics

Hancock and Cowling measured the critical crack tip opening displacements, δ_f, at fracture initiation in HY-80 steel specimens of six different configurations. δ_f varied from 90 μm in a deeply double-edge-cracked tensile panel to 900 μm in a single-edge-cracked tensile panel.McMeeking and Parks, and Shih and German have shown by their finite element calculations that the characteristics of the plane strain crack tip fields in both large scale yielding and general yielding are strongly dependent on specimen geometry and load level.In this study, the plane strain crack tip fields in the specimens tested by Hancock and Cowling were calculated using the finite element method. The crack tip triaxial tensile stress field is strongly affected by specimen geometric constraint, and the state of the triaxial tensile stress in a crack tip region is monitored by the ratio between the local tensile stress and the effective stress, i.e., ( ), at a distance x=2δ from the crack tip. The values of ( ) vary from 3.1 for the double-edge-cracked tensile panel to 1.7 for the single-edge-cracked tensile panel. The δ_f measured by Hancock and Cowling correlates very well with the ratio ( ). δ_f is a measure of the fracture ductility of the material ahead of the crack tip, and the ductility decreases with an increase in the triaxial tensile stress, i.e., the ratio ( ).     

Damage model for viscoplastic-softening behavior: cement paste and concrete

A viscoplastic-softening model is developed; it invokes damage accumulation depending on the viscous strain and stress rates. For deformation beyond the peak on the uniaxial stress-strain curve, the softening behavior is modelled by applying the accounting for loss in stiffness due to localized material damage by cracking. Predicted are the hardening/softening behavior of cement paste. The results for applied strain rates of 3 × 10⁻³, 3 × 10⁻² and 3 × 10⁻¹ s⁻¹ agreed well with the test data. Similar success was obtained for the creep of two types of concrete under compression.     

Fatigue and fretting fatigue of ion-nitrided 34CrNiMo6 steel

This work is concerned with the surface treatment (ion nitriding) of fretting fatigue and fatigue resistance of 34CrNiMo6. Tests are made on a servo-hydraulic machine under tension for both treated and non-treated specimens. The test parameters involve the applied displacements δ±80–±170 μm; fretting pressure σ_n=1000–1400 MPa; fatigue stress amplitude σ_a=380–680 MPa and stress ratio R=−1. The ion nitriding process improves both fatigue and fretting fatigue lives. Subsurface crack initiation from internal discontinuities was found for ion-nitrided specimens.     

Loading rate spectra for fracture initiation in metals

In the first part of the paper a review of experimental results is presented on the effects of loading rate of the fracture initiation toughness of structural alloys. Recent progress in new experimental techniques enables for measuring the fracture initiation properties over ten orders of magnitude in , where loading rate for plane strain fracture toughness. The range covered in experiments is

Full-size image (1K)

.Experimental results for aluminum and titanium alloys, as well as for some steels are discussed. Those results are shown as loading rate spectra where the fracture toughness K_Ic is plotted versus log at constant temperature. The patterns of the loading rate spectra demonstrate a negative role of the strain rate sensitivity in the process of fracture initiation, at least within a certain ranges of the loading rate .In the second part an attempt is presented to model the observed phenomena. Consequently, a new model has been proposed which combines an semi-empirical correlation between the critical fracture stress σ_F, yield stress σ_y from one side and a simple thermally activated potential of plastic flow from the other. It has been demonstrated that this simple approach can predict properly a fundamental changes in the loading rate spectra for steels.     

Stretch and twist of rotating prismatic beam

A turbine blade is modelled as a uniform isotropic prismatic beam of general cross-section and “setting angle” rotating about one end, and is analysed according to the linear theory of elasticity. A semi-inverse solution is presented which reduces the three-dimensional problem to one of two dimensions, and explicit stress and strain components given for the mathematically amenable elliptic cross-section. As expected, the planar stresses σ_x,σ_y, and τ_xy arising from the two-dimensional problem are found to be small. For the general section, the theory predicts small curvature of the blade centre line, and a twisting moment which tends to reduce the “angle of set”.     

Mechanics of adhesive contact on a power-law graded elastic half-space

We consider adhesive contact between a rigid sphere of radius R and a graded elastic half-space with Young's modulus varying with depth according to a power law E=E₀(z/c₀)^k (0<k<1) while Poisson's ratio ν remaining a constant. Closed-form analytical solutions are established for the critical force, the critical radius of contact area and the critical interfacial stress at pull-off. We highlight that the pull-off force has a simple solution of P_cr=−(k+3)πRΔγ/2 where Δγ is the work of adhesion and make further discussions with respect to three interesting limits: the classical JKR solution when k=0, the Gibson solid when k→1 and ν=0.5, and the strength limit in which the interfacial stress reaches the theoretical strength of adhesion at pull-off.     

Spreading of inkjet droplet of non-Newtonian fluid on solid surface with controlled contact angle at low Weber and Reynolds numbers

The dynamics of inkjet droplet of non-Newtonian fluid on glass substrates was investigated experimentally and compared with that of Newtonian fluid. The non-Newtonian fluids used here were 100 ppm solutions of polyethylene oxide (300k, 600k and 900k) dissolved in the 1:1 mixture of water and glycerin. Weber number (We) was 2–35 and Ohnesorge number was fixed at 0.057 ± 0.003. The wettability of solid substrate was also varied. The diameter of inkjet droplets in the present study was about 50 μm and was much smaller than the size of the previous studies on drop impact. Due to the development of a thin and long thread at the rear of the main drop the jetting window of polymer solution was much narrower than that of Newtonian fluid, and hence the experimental range of Weber number was restricted. The impact scenarios of non-Newtonian inkjet droplets were found to be qualitatively different from those of Newtonian droplets during the receding phase while they were almost the same as the Newtonian fluid case during the kinematic phase. The spreading diameter at the equilibrium was well correlated with the modified Weber number (We′ = We/(1 − cos θ_eq)) as in the case of Newtonian fluid, where θ_eq is the equilibrium contact angle. The similarity or disparity between the Newtonian and non-Newtonian cases was discussed considering the conformation of polymer chains during each stage of drop deformation.     

Phenomenological study of parabolic and spherical indentation of elastic-ideally plastic material

A phenomenological study of parabolic and spherical indentation of elastic ideally plastic materials was carried out by using precise results of finite elements calculations. The study shows that no “pseudo-Hertzian” regime occurs during spherical indentation. As soon as the yield stress of the indented material is exceeded, a deviation from the, purely elastic Hertzian contact behaviour is found. Two elastic–plastic regimes and two plastic regimes are observed for materials of very large Young modulus to Yield stress ratio, E/σ_y. The first elastic–plastic regime corresponds to a strong evolution of the indented plastic zone. The first plastic regime corresponds to the commonly called “fully plastic regime”, in which the average indentation pressure is constant and equal to about three times the yield stress of the indented material. In this regime, the contact depth to penetration depth ratio tends toward a constant value, i.e. h_c/h = 1.47. h_c/h is only constant for very low values of yield strain (σ_y/E lower than 5 × 10^?6) when aE¹/Rσ_y is higher than 10,000. The second plastic regime corresponds to a decrease in the average indentation pressure and to a steeper increase in the pile-up. For materials with very large E/σ_y ratio, the second plastic regime appears when the value of the non-dimensional contact radius a/R is lower than 0.01. In the case of spherical and parabolic indentation, results show that the first plastic regime exists only for elastic-ideally plastic materials having an E/σ_y ratio higher than approximately 2.000.     

A centrally cracked thin circular disk, Part I: 3D Elastic–plastic finite element analysis

A full field solution, based on small deformation, three-dimensional elastic–plastic finite element analysis of the centrally cracked thin disk under mode I loading has been performed. The solution for the stresses under small-scale yielding and lo!cally fully plastic state has been compared with the HRR plane stress solution. At the outside of the 3D zone, within a distance of rσ_o/J=18, HRR dominance is maintained in the presence of a significant amount of compressive stress along the crack flanks. Ahead of this region, the HRR field overestimate the stresses. These results demonstrate a completely reversed state of stress in the near crack front compared to that in the plane strain case. The combined effect of geometry and finite thickness of the specimen on elastic–plastic crack tip stress field has been explored. To the best of our knowledge, such an attempt in the published literature has not been made yet. For the qualitative assessment of the results some of the field parameters have been compared to the available experimental results of K, gives a fair estimate of the crack opening stress near the crack front at a distance of order 10⁻² in. On the basis of this analysis, the Linear Elastic Fracture Mechanics approach has been adopted in analyzing the fatigue crack extension experiments performed in the disk (Part II).     

Near tip field for pseudo Mode II crack growth: Power law hardening material

The asymptotic behavior of stress and strain near the tip of a Mode II crack growing in power law hardening material is analyzed by assuming that the crack grows straight ahead even though tests show otherwise. The results show that the stress and strain possess the singularities of (ln r)^2/(n−1) and (ln r)^2n/(n−1) respectively. The distance from the crack tip is r, and n is the hardening exponent, i.e. σⁿ. The amplitudes of the stress and strain near the crack tip are determined by the asymptotic analysis.