Finite Element Solutions for Plane Strain Mode I Crack with Strain Grading Effect

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom omega(i) are introduced in addition to the conventional three translational degrees of freedom mu(i). omega(i) is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l(1). Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale. (C) 2002 Elsevier Science Ltd. All rights reserved.

Torsional impact of transversely isotropic solid with functionally graded shear moduli and a penny-shaped crack

The torsional impact response of a penny-shaped crack in an unbounded transversely isotropic solid is considered. The shear moduli are assumed to be functionally graded such that the mathematics is tractable. Laplace transform and Hankel transform are used to reduce the problem to solving a Fredholm integral equation. The crack tip stress fields are obtained. Investigated are the influence of material nonhomogeneity and orthotropy on the dynamic stress intensity factor. The peak value of the dynamic stress intensity factor can be suppressed by increasing the shear moduli's gradient and/or increasing the shear modulus in a direction perpendicular to the crack surface.

Mode I and Mode II Crack Tip Asymptotic Fields with Strain Gradient Effects

The strain gradient effect becomes significant when the size of fracture process zone around a crack tip is comparable to the intrinsic material length l, typically of the order of microns. Using the new strain gradient deformation theory given by Chen and Wang, the asymptotic fields near a crack tip in an elastic-plastic material with strain gradient effects are investigated. It is established that the dominant strain field is irrotational. For mode I plane stress crack tip asymptotic field, the stress asymptotic field and the couple stress asymptotic field can not exist simultaneously. In the stress dominated asymptotic field, the angular distributions of stresses are consistent with the classical plane stress HRR field; In the couple stress dominated asymptotic field, the angular distributions of couple stresses are consistent with that obtained by Huang et al. For mode II plane stress and plane strain crack tip asymptotic fields, only the stress-dominated asymptotic fields exist. The couple stress asymptotic field is less singular than the stress asymptotic fields. The stress asymptotic fields are the same as mode II plane stress and plane strain HRR fields, respectively. The increase in stresses is not observed in strain gradient plasticity for mode I and mode II, because the present theory is based only on the rotational gradient of deformation and the crack tip asymptotic fields are irrotational and dominated by the stretching gradient.

The effect of straining mode on void growth

Large strain finite element method is employed to investigate the effect of straining mode on void growth. Axisymmetric cell model embedded with spherical void is controlled by constant triaxiality: loading,while plane-stress model containing a circular void is loaded by constant ratio of straining. Elastic-plastic material is used for the matrix in both cases. It is concluded that, besides the known effect of triaxiality, the straining mode which intensifies the plastic concentration around the void is also a void growth stimulator. Experimental results are cited to justify the computation results.

Moment Tensor Analysis of the Acoustic Emission Source in the Rock Damage Process

To further investigate the mechanism of acoustic emission (AE) in the rock fracture experiment, moment tensor analysis was carried out. The AE sources characterized by crack sizes, orientations and fracture modes, are represented by a time-dependent momen

Mode I crack tip field with strain gradient effects

A plane strain mode I crack tip field with strain gradient effects is investigated. A new strain gradient theory is used. An elastic-power law hardening strain gradient material is considered and two hardening laws, i.e. a separation law and an integration Law are used respectively. As for the material with the separation law hardening, the angular distributions of stresses are consistent with the HRR field, which differs from the stress results([19]); the angular distributions of couple stresses are the same as the couple stress results([19]). For the material with the integration law hardening, the stress field and the couple stress field can not exist simultaneously, which is the same as the conclusion([19]), but for the stress dominated field, the angular distributions of stresses are consistent with the HRR field; for the couple stress dominated field, the angular distributions of couple stresses are consistent with those in Ref. [19]. However, the increase in stresses is not observed in strain gradient plasticity because the present theory is based on the rotation gradient of the deformation only, while the crack tip field of mode I is dominated by the tension gradient, which will be shown in another paper.

Investigation Of Mode Ii Crack Growth Following A Very High Speed Impact

A recoverable plate impact testing technology has been developed for studying fracture mechanisms of mode II crack. With this technology, a single duration stress pulse with submicrosecond duration and high loading rates, up to 10(8) MPam(1/2)s(-1), can be produced. Dynamic failure tests of Hard-C 60# steel were carried out under asymmetrical impacting conditions with short stress-pulse loading. Experimental results show that the nucleation and growth of several microcracks ahead of the crack tip, and the interactions between them, induce unsteady crack growth. Failure mode transitions during crack growth, both from mode I crack to mode II and from brittle to ductile fracture, were observed. Based on experimental observations, a discontinuous crack growth model was established. Analysis of the crack growth mechanisms using our model shows that the shear crack extension is unsteady when the extending speed is between the Rayleigh wave speed c(R) and the shear wave speed c(S). However, when the crack advancing speed is beyond c(S), the crack grows at a steady intersonic speed approaching root 2c(S). It also shows that the transient mechanisms, such as nucleation, growth, interaction and coalescence among microcracks, make the main crack speed jump from subsonic to intersonic and the steady growth of all the subcracks causes the main crack to grow at a stable intersonic speed.

Concentrated-Mass Cantilever Enhances Multiple Harmonics In Tapping-Mode Atomic Force Microscopy

The natural frequencies of a cantilever probe can be tuned with an attached concentrated mass to coincide with the higher harmonics generated in a tapping-mode atomic force microscopy by the nonlinear tip-sample interaction force. We provide a comprehensive map to guide the choice of the mass and the position of the attached particle in order to significantly enhance the higher harmonic signals containing information on the material properties. The first three eigenmodes can be simultaneously excited with only one carefully positioned particle of specific mass to enhance multiple harmonics. Accessing the interaction force qualitatively based on the high-sensitive harmonic signals combines the real-time material characterization with the imaging capability. (C) 2008 American Institute of Physics.

Non-Slipping Jkr Model For Transversely Isotropic Materials

A generalized plane strain JKR model is established for non-slipping adhesive contact between an elastic transversely isotropic cylinder and a dissimilar elastic transversely isotropic half plane, in which a pulling force acts on the cylinder with the pulling direction at an angle inclined to the contact interface. Full-coupled solutions are obtained through the Griffith energy balance between elastic and surface energies. The analysis shows that, for a special case, i.e., the direction of pulling normal to the contact interface, the full-coupled solution can be approximated by a non-oscillatory one, in which the critical pull-off force, pull-off contact half-width and adhesion strength can be expressed explicitly. For the other cases, i.e., the direction of pulling inclined to the contact interface, tangential tractions have significant effects on the pull-off process, it should be described by an exact full-coupled solution. The elastic anisotropy leads to an orientation-dependent pull-off force and adhesion strength. This study could not only supply an exact solution to the generalized JKR model of transversely isotropic materials, but also suggest a reversible adhesion sensor designed by transversely isotropic materials, such as PZT or fiber-reinforced materials with parallel fibers. (c) 2007 Elsevier Ltd. All rights reserved.

Dynamic mode I perturbation solution for a moving crack unsteadily

Rice et al. (Jounal of Mechanics and Physics of Solids 42, 813-843) analyze the propagation of a planar crack with a nominally straight front in a model elastic solid with a single displacement component. Using the form of Willis er al. (Journal of the Mechanics and Physics of Solids 43, 319-341), of dynamic mode I weight functions for a moving crack, we address that problem solved by Rice ei al. in the 3D context of elastodynamic theory. Oscillatory crack tip motion results from constructive-destructive interference of stress intensity waves. Those waves, including system of the dilatational, shear and Rayleigh waves, interact on each other and with moving edge of crack, can lead to continuing fluctuations of the crack front and propagation velocity. (C) 1997 Elsevier Science Ltd.

Stress intensity factor histories of a steadily propagating crack under 3-D combined mode traction

Elastodynamic stress intensity factor histories of an unbounded solid containing a semi-infinite plane crack that propagates at a constant velocity under 3-D time-independent combined mode loading are considered. The fundamental solution, which is the response of point loading, is obtained. Then, stress intensity factor histories of a general loading system are written out in terms of superposition integrals. The methods used here are the Laplace transform methods in conjunction with the Wiener-Hopf technique.

On the mechanism of void growth and the effect of straining mode in ductile materials

Finite element analysis is employed to investigate void growth embedded in elastic-plastic matrix material. Axisymmetric and plane stress conditions are considered. The simulation of void growth in a unit cell model is carried out over a wide range of triaxial tensile stressing or large plastic straining for various strain hardening materials to study the mechanism of void growth in ductile materials. Triaxial tension and large plastic strain encircling around the void are found to be of most importance for driving void growth. The straining mode of incremental loading which favors the necessary strain concentration around void for its growth can be characterized by the vanishing condition of a parameter called "the third invariant of generalized strain rate". Under this condition, it accentuates the internal strain concentration and the strain energy stored/dissipated within the material layer surrounding the void. Experimental results are cited to justify the effect of this loading parameter. (C) 2000 Elsevier Science Ltd. All rights reserved.

Dynamical symmetry of screened Coulomb potential and isotropic harmonic oscillator

It is shown that for the screened Coulomb potential and isotropic harmonic oscillator, there exists an infinite number of closed orbits for suitable angular momentum values. At the aphelion (perihelion) points of classical orbits, an extended Runge-Lenz vector for the screened Coulomb potential and an extended quadrupole tensor for the screened isotropic harmonic oscillator are still conserved. For the screened two-dimensional (2D) Coulomb potential and isotropic harmonic oscillator, the dynamical symmetries SO3 and SU(2) are still preserved at the aphelion (perihelion) points of classical orbits, respectively. For the screened 3D Coulomb potential, the dynamical symmetry SO4 is also preserved at the aphelion (perihelion) points of classical orbits. But for the screened 3D isotropic harmonic oscillator, the dynamical symmetry SU(2) is only preserved at the aphelion (perihelion) points of classical orbits in the eigencoordinate system. For the screened Coulomb potential and isotropic harmonic oscillator, only the energy (but not angular momentum) raising and lowering operators can be constructed from a factorization of the radial Schrodinger equation.

Effect of dentin tubules on the mechanical properties of dentin - part III: Numerical, analysis

Imitating a real tooth and the periodontal supporting tissues, we have established a 2D finite element model and carried out a numerical analysis based on the inhomogeneous and anisotropic (IA) stress-strain relation and strength model of dentin proposed in the preceding Parts I and II, and the conventional homogeneous and isotropic (III) model, respectively. Quite a few cases of loadings for a non-defected and a defected tooth are considered. The numerical results show that the stress level predicted by the IA model is remarkably higher than that by the III model, revealing that the effect of the dentin tubules should be taken into a serious consideration from the viewpoint of biomechanics.

Long-range automaton models of earthquakes: Power-law accelerations, correlation evolution, and mode-switching

We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r(-p) dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.