Kinetics of microcrack blunting ahead of macrocrack approaching shear wave speed

The fracturing of glass and tearing of rubber both involve the separation of material but their crack growth behavior can be quite different, particularly with reference to the distance of separation of the adjacent planes of material and the speed at which they separate. Relatively speaking, the former and the latter are recognized, respectively, to be fast and slow under normal conditions. Moreover, the crack tip radius of curvature in glass can be very sharp while that in the rubber can be very blunt. These changes in the geometric features of the crack or defect, however, have not been incorporated into the modeling of running cracks because the mathematical treatment makes use of the Galilean transformation where the crack opening distance or the change in the radius of curvature of the crack does not enter into the solution. Change in crack speed is accounted for only via the modulus of elasticity and mass density. For this simple reason, many of the dynamic features of the running crack have remained unexplained although speculations are not lacking. To begin with, the process of energy dissipation due to separation is affected by the microstructure of the material that distinguishes polycrystalline from amorphous form. Energy extracted from macroscopic reaches of a solid will travel to the atomic or smaller regions at different speeds at a given instance. It is not clear how many of the succeeding size scales should be included within a given time interval for an accurate prediction of the macroscopic dynamic crack characteristics. The minimum requirement would therefore necessitate the simultaneous treatment of two scales at the same time. This means that the analysis should capture the change in the macroscopic and microscopic features of a defect as it propagates. The discussion for a dual scale model has been invoked only very recently for a stationary crack. The objective of this work is to extend this effort to a crack running at constant speed beyond that of Rayleigh wave. Developed is a dual scale moving crack model containing microscopic damage ahead of a macroscopic crack with a gradual transition. This transitory region is referred to as the mesoscopic zone where the tractions prevail on the damaged portion of the material ahead of the original crack known as the restraining stresses, the magnitude of which depends on the geometry, material and loading. This damaged or restraining zone is not assumed arbitrarily nor assumed to be intrinsically a constant in the cohesive stress approach; it is determined for each step of crack advancement. For the range of micronotch bluntness with 0 < β < 30° and 0.2 σ_∞/σ₀ 0.5, there prevails a nearly constant restraining zone size as the crack approaches the shear wave speed. Note that β is the half micronotch angle and the applied stress ratio is σ_∞/σ₀ with σ₀ being the maximum of the restraining stress. For σ_∞/σ₀ equal to or less than 0.5, the macrocrack opening displacement COD is nearly constant and starts to decrease more quickly as the crack approaches the shear wave speed. For the present dual scale model where the normalized crack speed v/c_s increases with decreasing with the one-half microcrack tip angle β. There prevails a limit of crack tip bluntness that corresponds to β 36° and v/c_s 0.15. That is a crack cannot be maintained at a constant speed if the bluntness is increased beyond this limiting value. Such a feature is manifestation of the dependency of the restraining stress on crack velocity and the applied stress or the energy pumped into the system to maintain the crack at a constant velocity. More specifically, the transitory character from macro to micro is being determined as part of the unknown solution. Using the energy density function dW/dV as the indicator, plots are made in terms of the macrodistance ahead of the original crack while the microdefect bluntness can vary depending on the tip geometry. Such a generality has not been considered previously. The macro-dW/dV behavior with distance remains as the inverse r relation yielding a perfect hyperbola for the homogeneous material. This behavior is the same as the stationary crack. The micro-dW/dV relations are expressed in terms of a single undetermined parameter. Its evaluation is beyond the scope of this investigation although the qualitative behavior is expected to be similar to that for the stationary crack. To reiterate, what has been achieved as an objective is a model that accounts for the thickness of a running crack since the surface of separation representing damage at the macroscopic and microscopic scale is different. The transitory behavior from micro to macro is described by the state of affairs in the mesoscopic zone.     

Characterization of material inhomogeneity by stationary values of strain energy density

The strain energy density theory is applied to analyze the fracture instability of a mechanical system whose behavior is governed by the interaction of geometry, load and material inhomogeneity. This is accomplished by obtaining the location of the global and local relative minima of the strain energy density function dW/dV denoted, respectively, by [(dW/dV)_min]_g and [(dW/dV)_min]_¢l. The former refers to a fixed global coordinate system for the entire solid while the latter corresponds to local coordinate systems referred to each material point. An unique length parameter “ℓ” representing the distance between [(dW/dV)_min^max]_g and [(dW/dV)_min^max]_ℓ can thus be found and serves as a measure of the degree of system instability tending toward failure by fracture.Numerical results are obtained and displayed graphically for the case of a solid containing an inclusion of dissimilar material. The changes that take place in material inhomogeneity, loading type and physical dimensions of the solid and inclusion are reflected through ℓ. The method suggests the compatibility of ℓ for each member of a multi-component structure in order to avoid premature failure of a single member.     

Piezomagnetic and piezoelectric poling effects on mode I and II crack initiation behavior of magnetoelectroelastic materials

The ferrite and ferroelectric phase of magnetoelectroelastic (MEE) material can be selected and processed to control the macroscopic behavior of electron devices using continuum mechanics models. Once macro- and/or microdefects appear, the highly intensified magnetic and electric energy localization could alter the response significantly to change the design performance. Alignment of poling directions of piezomagnetic and piezoelectric materials can add to the complexity of the MEE material behavior to which this study will be concerned with.Appropriate balance of distortional and dilatational energy density is no longer obvious when a material possesses anisotropy and/or nonhomogeneity. An excess of the former could result in unwanted geometric change while the latter may lead to unexpected fracture initiation. Such information can be evaluated quantitatively from the stationary values of the energy density function dW/dV. The maxima and minima have been known to coincide, respectively, with possible locations of permanent shape change and crack initiation regardless of material and loading type. The direction of poling with respect to a line crack and the material microstructure described by the constitutive coefficients will be specified explicitly with reference to the applied magnetic field, electric field and mechanical stress, both normal and shear. The crack initiation load and direction could be predicted by finding the direction for which the volume change is the largest. In contrast to intuition, change in poling directions can influence the cracking behavior of MEE dramatically. This will be demonstrated by the numerical results for the BaTiO₃–CoFe₂O₄ composite having different volume fractions where BaTiO₃ and CoFe₂O₄ are, respectively, the inclusion and matrix.To be emphasized is that mode I and II crack behavior will not have the same definition as that in classical fracture mechanics where load and crack extension symmetry would coincide. A striking result is found for a mode II crack. By keeping the magnetic poling fixed, a reversal of electric poling changed the crack initiation angle from θ₀=+80° to θ₀=−80° using the line extending ahead of the crack as the reference. This effect is also sensitive to the distance from the crack tip. Displayed and discussed are results for r/a=10⁻⁴ and 10⁻¹. Because the theory of magnetoelectroelasticity used in the analysis is based on the assumption of equilibrium where the influence of material microstructure is homogenized, the local space and temporal effects must be interpreted accordingly. Among them are the maximum values of (dW/dV)_max and (dW/dV)_min which refer to as possible sites of yielding and fracture. Since time and size are homogenized, it is implicitly understood that there is more time for yielding as compared to fracture being a more sudden process. This renders a higher dW/dV in contrast to that for fracture. Put it differently, a lower dW/dV with a shorter time for release could be more detrimental.     

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.     

Pseudo plane stress mode II crack growth in plastic materials

The near tip field of mode II crack that grows in thin bodies with power hardening or perfectly plastic behavior is analyzed. It is shown that for power hardening behavior, the pseudo plane stress field possesses the logarithm singularity, i.e. σ (ln r)²/(n−1), (ln r)^{2n/(n − 1)}, where r is the distance from the crack tip, n the hardening exponent is σⁿ. When n → ∞ the solution reduced to that for the perfectly plastic case.     

On corner flows of Oldroyd-B fluids

Exact series solutions for planar creeping flows of Oldroyd-B fluids in the neighbourhood of sharp corners are presented and discussed. Both reentrant and non-reentrant sectors are considered. For reentrant sectors it is shown that more than one type of series solution can exist formally, one type exhibiting Newtonian-like asymptotic behaviour at the corner, away from walls, and another type exhibiting the same kind of asymptotics as an Upper Convected Maxwell (UCM) fluid. The solutions which are Newtonian-like away from walls are shown to develop non-integrable stress singularities at the walls when the no-slip velocity boundary condition is imposed. These mathematical solutions are therefore inadmissible from the physical viewpoint under no-slip conditions. An inadmissible solution, with stress singularities which are not everywhere integrable, is identified among the solutions of UCM-type. For a 270° reentrant sector the radial behaviour of the normal stress is everywhere r^−0.613. In the viscometric region near a wall, the radial normal stress σ^rr behaves like (rε)^−0.613, where ε is the angle made with the wall. In addition σ^rθ is infinite (not integrable) at the wall even when r is non-zero. Another UCM-type solution has a normal stress behaviour away from walls which is r^−0.985 for 270° sector. Again, this solution has a non-integrable stress singularity and is therefore inadmissible. Finally, for non-reentrant sectors it is shown that the flow is always Newtonian-like away from walls.     

Fracture instability of reinforced panel with cracks emanating from hole

Initiation of failure by yielding and/or fracture depends on the magnitude of the distortion and dilatation of material elements. According to the strain energy density theory (SED), failure is assumed to initiate at the site of the local maximum of maxima [(dW/dV)^max_max]_L by yielding and the maximum of minima [(dW/dV)^max_min]_L by fracture. The fracture is assumed to start from point L where [(dW/dV)^max_min]_L appears and tends toward G where the global maximum of dW/dV minima appears, denoted by [(dW/dV)^max_min]_G. The distance l between L and G along the anticipated crack trajectory is an indication of failure instability of the system by fracture. If l is sufficiently large and [(dW/dV)^max_min]_L exceeds the threshold, fracture initiation could lead to global failure. The local and global failure instability of a composite structural component is studied by application of the strain energy density theory. The depicted configuration is that of a panel with a circular hole reinforced by two side strips made of different material. The case of two symmetric cracks emanating from the hole and normal to the applied uniaxial tensile stress is also analyzed. Results are displayed graphically to illustrate the geometry and dissimilar material properties influence the fracture instability behavior of the two examples.     

An Instrumented Indentation Method for Young’s Modulus Measurement with Accuracy Estimation

The relationships between indentation responses and Young’s modulus of an indented material were investigated by employing dimensional analysis and finite element method. Three representative tip bluntness geometries were introduced to describe the shape of a real Berkovich indenter. It was demonstrated that for each of these bluntness geometries, a set of approximate indentation relationships correlating the ratio of nominal hardness/reduced Young’s modulus H _n /E _r and the ratio of elastic work/total work W _e/W can be derived. Consequently, a method for Young’s modulus measurement combined with its accuracy estimation was established on basis of these relationships. The effectiveness of this approach was verified by performing nanoindentation tests on S45C carbon steel and 6061 aluminum alloy and microindentation tests on aluminum single crystal, GCr15 bearing steel and fused silica.     

Micro/macro-crack growth due to creep–fatigue dependency on time–temperature material behavior

The implicit character of micro-structural degradation is determined by specifying the time history of crack growth caused by creep–fatigue interaction at high temperature. A dual scale micro/macro-equivalent crack growth model is used to illustrate the underlying principle of multiscaling which can be applied equally well to nano/micro. A series of dual scale models can be connected to formulate triple or quadruple scale models. Temperature and time-dependent thermo-mechanical material properties are developed to dictate the design time history of creep–fatigue cracking that can serve as the master curve for health monitoring.In contrast to the conventional procedure of problem/solution approach by specifying the time- and temperature-dependent material properties as a priori, the desired solution is then defined for a class of anticipated loadings. A scheme for matching the loading history with the damage evolution is then obtained. The results depend on the initial crack size and the extent of creep in proportion to fatigue damage. The path dependent nature of damage is demonstrated by showing the range of the pertinent parameters that control the final destruction of the material. A possible scenario of 20 yr of life span for the 38Cr2Mo2VA ultra-high strength steel is used to develop the evolution of the micro-structural degradation. Three micro/macro-parameters μ^*, d^* and σ^* are used to exhibit the time-dependent variation of the material, geometry and load effects. They are necessary to reflect the scale transitory behavior of creep–fatigue damage. Once the algorithm is developed, the material can be tailor made to match the behavior. That is a different life span of the same material would alter the time behavior of μ^*, d^* and σ^* and hence the micro-structural degradation history. The one-to-one correspondence of the material micro-structure degradation history with that of damage by cracking is the essence of path dependency. Numerical results and graphs are obtained to demonstrate how the inherently implicit material micro-structure parameters can be evaluated from the uniaxial bulk material properties at the macroscopic scale.The combined behavior of creep and fatigue can be exhibited by specifying the parameter ξ with reference to the initial defect size a₀. Large ξ (0.90 and 0.85) gives critical crack size a_cr = 11–14 mm (at t < 20 yr) for a₀ about 1.3 mm. For small ξ (0.05 and 0.15), there results critical a_cr = 6–7 mm (at t < 20 yr) for a₀ about 0.7–0.8 mm. The initial crack is estimated to increase its length by an order of magnitude before triggering global to the instability. This also applies ξ ≈ 0.5 where creep interacts severely with fatigue. Fine tuning of a_cr and a₀ can be made to meet the condition oft = 20 yr.Trade off among load, material and geometric parameters are quantified such that the optimum conditions can be determined for the desired life qualified by the initial–final defect sizes. The scenario assumed in this work is indicative of the capability of the methodology. The initial–final defect sizes can be varied by re-designing the time–temperature material specifications. To reiterate, the uniqueness of solution requires the end result to match with the initial conditions for a given problem. This basic requirement has been accomplished by the dual scale micro/macro-crack growth model for creep and fatigue.