An extension of the J-integral to piezoelectric materials

This paper provides an analysis of the crack propagation criterion for a thin piezoelectric plate with a symmetry of order six. On the basis of Gol’denveizer’s asymptotic integration method or Destuynder’s unidirectional zoom technique, we obtain an extension of the purely mechanical J-integral to piezoelectric materials, with a dependence of the gradient of energy of the plate only on zeroth order terms of asymptotic expansions.     

淮河大堤老应段土体蠕变特性研究及工程应用

本文依据工程需要, 采取淮河大堤老应段的原状和扰动土样, 在室内对堤体及堤基土的剪切蠕变和拉伸蠕变特性进行了试验研究。粘土的蠕变剪切强度远低于瞬时剪切强度, Ｃ∞值仅是Ｃ值的４０％左右;相同含水条件下粘土的长时抗拉强度大于瞬时抗拉强度, 长时拉伸应变量远大于瞬时拉伸应变量, 约为２倍。     

Interface crack between two shear deformable elastic layers

An improved method based on the first-order shear deformable plate theory is developed to calculate the energy release rate and stress intensity factor for a crack at the interface of a bi-layer structure. By modeling the uncracked region of the structure as two separate Reissner-Mindlin plates bonded perfectly along the interface, this method is able not only to take into account the shear deformation in the cracked region, but also to capture the shear deformation in the uncracked region of the structure. A closed form solution of energy release rate and mode decomposition at the interface crack is obtained for a general loading condition, and it indicates that the energy release rate and stress intensity factor are determined by two independent loading parameters. Compared to the approach based on the classical plate theory, the proposed method provides a more accurate prediction of energy release rate as well as mode decomposition. The computational procedures introduced are relatively straightforward, and the closed form solution can be used to predict crack growth along the layered structures.     

Application of material forces to fracture of inhomogeneous materials: illustrative examples

The material forces concept has become an elegant tool in continuum mechanics for the calculation of the thermodynamic driving force of a defect. Based on this concept, we have recently shown that inhomogeneities essentially shield or anti-shield crack tips from applied far-field stresses. The goal of this paper is to illustrate this by considering the model example of a crack in a CT-type specimen that contains a bimaterial interface. The crack driving force is calculated as the sum of the far-field driving force and the crack-tip shielding or anti-shielding. Several cases of inhomogeneity in either thermal or elastic properties are considered. Rather simple hand calculations are provided in addition to numerical results to illustrate the advantages of using the material forces concept.     

Numerical plane stress elastic–perfectly plastic crack analysis under Tresca yield condition with comparison to Dugdale plastic strip model

A number of plane stress numerical analyses of the mode I elastoplastic fracture mechanics problem have been performed in the past using the Huber–Mises yield criterion. This study employs instead the Tresca yield condition using an incremental theory of plasticity for a stationary crack. A commercial finite element program is used to solve the opening mode of fracture problem (mode I) for a square plate containing a central crack under generalized plane stress loading conditions. A biaxial uniform tensile traction is applied to the edges of a thin plate composed of a linear elastic non-work hardening material under small strain assumptions. The finite element results are compared with the analytical predictions of the Dugdale plastic strip model for a crack in an infinite plate subject to a biaxial uniform load at infinity.     

Crack propagation in an orthotropic medium with coupled elastodynamic properties

In this paper the Mode-I elastodynamic problem of a crack propagating in an orthotropic medium is studied under the condition that the matrix of elastodynamic coefficients has repeated eigenvalues. It is shown that the crack is constrained in an elastodynamic state which is defined through a compulsory condition coupling its velocity with the elastic parameters of the orthotropic medium. The dynamic stress and displacement components ahead of the crack tip as well as the energy release rate are explicitly obtained.     

Mixed-mode I/II fields around a crack with a cohesive zone ahead of the crack tip

Considering a cohesive zone ahead of the crack tip, the linear elastic crack problem under Mode I, Mode II or mixed-mode conditions is formulated in an elliptic coordinate system, so that the cohesive surfaces are conveniently represented by straight line segments. It is shown that the displacement and stress fields around the crack tip and the cohesive zone, expressed in terms of elliptic coordinates, have a simple mathematical form, which does not contain a stress singularity at the crack tip due to the existence of the cohesive zone.     

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.     

From crack deflection to lattice vibrations—macro to atomistic examination of dynamic cleavage fracture

A set of cleavage experiments with strip-shaped single-crystal silicon specimens subjected to three-point bending is reported. The experiments enabled examination of the relationships between the dynamic energy release rate, the velocity, the orientation-dependent cleavage energy, and the cleavage plane of propagation.Dynamic crack propagation experiments show that when a [0 0 1] silicon single crystal is fractured under three-point bending at ‘parallel’ velocity (directly measured at the bottom surface of the specimen) of up to , it prefers to cleave along the vertical (1 1 0) plane, while when the specimen is fractured under the same conditions but at a velocity higher than , it cleaves along the inclined (1 1 1) plane. At intermediate velocities, the crack will deflect from the (1 1 0) plane to the (1 1 1) plane. Crack velocity was determined by the initial notch length. The local (calculated) velocity of deflection between the cleavage planes ranges from , for a crack propagating on the (1 1 0) plane in the direction, to about , for a crack on the (1 1 0) plane, but in the [0 0 1] direction.It is suggested that the cause of the deflection phenomenon is the anisotropic, velocity-dependent cleavage energy, resulted phonon radiation caused by anisotropic, velocity-dependent lattice vibrations. We have studied the effect of material properties and propose selection criteria to explain the deflection phenomenon: the crack will deflect to the plane of least-energy, for which G−Γ_i(V)=max, or to the plane with maximum crack tip velocity, V_i(Γ)=max.     

Multi-residue analysis of drugs of abuse in wastewater and surface water by solid-phase extraction and liquid chromatography-positive electrospray ionisation tandem mass spectrometry

A new-multi residue method was developed for the environmental monitoring of 65 stimulants, opiod and morphine derivatives, benzodiazepines, antidepressants, dissociative anaesthetics, drug precursors, human urine indicators and their metabolites in wastewater and surface water. The proposed analytical methodology offers rapid analysis for a large number of compounds, with low limits of quantification and utilises only one solid-phase extraction-ultra performance liquid chromatography-positive electrospray ionisation tandem mass spectrometry (SPE-LC-MS/MS) method, thus overcoming the drawbacks of previously published procedures. The method employed solid phase extraction with the usage of Oasis MCX sorbent and subsequent ultra performance liquid chromatography-positive electrospray ionisation tandem mass spectrometry. The usage of a 1.7 μm particle size column (1 mm×150 mm) resulted in very low flow rates (0.04 mLmin(-1)), and as a consequence gave good sensitivity, low mobile phase consumption and short retention times for all compounds (from 2.9 to 23.1 min). High SPE recoveries (>60%) were obtained for the majority of compounds. The mean correlation coefficients of the calibration curves were typically higher than 0.997 and showed good linearity in the range 0-1000 μgL(-1). The method limits of detection ranged from 0.1 ngL(-1) for compounds including cocaine, benzoylecgonine, norbenzoylecgonine and 2-oxo-3-hydroxy-LSD to 100 ngL(-1) for caffeine. Method quantification limits ranged from 0.5 to 154.2 ngL(-1). Intra- and inter-day repeatabilities were on average less than 10%. The method accuracy range was within -33.1 to 30.1%. The new multi-residue method was used to analyse drugs of abuse in wastewater and river water in the UK environment. Of the targeted 65 compounds, 46 analytes were detected at levels above the method quantification limit (MQL) in wastewater treatment plant (WWTP) influent, 43 in WWTP effluent and 36 compounds in river water.