Inhomogeneity effects on the crack driving force in elastic and elastic-plastic materials

This article evaluates the effect of material inhomogeneities on the crack-tip driving force in general inhomogeneous bodies and reports results for bimaterial composites. The theoretical model, based on Eshelby material forces, makes no assumptions about the distribution of the inhomogeneities or the constitutive properties of the materials. Inhomogeneities are modeled by making the stored energy have an explicit dependence on the reference coordinates. Then the material inhomogeneity effect on the crack-tip driving force is quantified by the term C_inh, which is the integral of the gradient of the stored energy in the direction of crack growth. The model is demonstrated by two model problems: (i) bimaterial elastic composite using asymptotic solutions and (ii) graded elastic and elastic-plastic compact tension specimen using numerical methods for stress analysis.     

Mechanics of rubber shear springs

FEA calculations have been carried out for a model rubber shear spring, consisting of a block of a highly elastic material, bonded between two rigid parallel plates and sheared by displacing one of the plates parallel to the other in its own plane. The block was prevented from deforming in the perpendicular direction, and thus was deformed in plane strain. Stress distributions along the bond-line and the center-line are reported and compared with those expected from the theory of large elastic deformations. Unexpected tensile stresses were found to develop in the interior of the sheared block. They are attributed to the absence on the end surfaces of the stresses needed to maintain a simple shear, causing a pronounced change in the reference pressure—a consequence that is usually overlooked. Because the internal stresses are governed by the boundary conditions, they were strongly affected by the shape of the end surfaces. In addition, they were reduced markedly by assigning values to Poisson's ratio slightly lower than 0.5, thus allowing some volume expansion of the rubber. Strain energy release rates were also evaluated for growth of a crack along the bond-line, starting at the edges, and compared with those reported previously by Lindley and Teo [Energy for crack growth at the bonds of rubber springs, Plast. Rubber Mat. Appl. 4 (1979) 29-37], Muhr et al. [A fracture mechanics study of natural rubber-to-metal bond failure, J. Adhes. Sci. Technol. 10 (1996) 593-616], Gregory and Muhr [Stiffness and fracture analysis of bonded rubber blocks in simple shear, in: D. Boast, V.A. Coveny (Eds.), Finite Element Analysis of Elastomers, Professional Engineering Publications, Bury St. Edmunds, UK, 1999, pp. 265-274] and Gough and Muhr [Initiation of failure of rubber close to bondlines, in: Proceedings of the International Rubber Conference, Maastricht, Netherlands, June 2005, IOM Communications Ltd., London, 2005, pp. 165-174]. They confirm that a long crack at the compression edge will grow faster than one at the tension edge, but the results for short cracks were inconclusive.     

Limitation on crack speed for a strip-craze model applied to PMMA

A strip-craze model is proposed to study crack propagation in polymers. A nonlinear differential equation is derived to govern the dynamic process of crack propagation. The viscous feature of the material in the craze zone is taken into account by means of an experimentally determined relationship between the craze stress and crack speed. By fitting experimental data of PMMA into the model, some parameters including the strip-craze length are deduced. A non-singular stress is introduced to control the crack propagation with a strip craze at its tip. Variations of the crack length and the crack speed with time are computed and their dependence on the non-singular stress is investigated. For PMMA, three stages of crack propagation are identified in terms of initial non-singular stress σ_ns0. When σ_ns0<60 MPa, the crack speed mm/s and the crack is basically stationary; when 60 <σ_ns0<95 MPa, then mm/s the crack is in slow propagation; when σ_ns0>95 MPa, then mm/s and the crack is in rapid propagation. The proposed model is applicable only in slow crack propagation.     

Measurement of energy release rate and energy flux of rapidly bifurcating crack in Homalite 100 and Araldite B by high-speed holographic microscopy

High-speed holographic microscopy is applied to take three successive photographs of fast propagating cracks in Homalite 100 or in Araldite B at the moment of bifurcation. Crack speed at bifurcation is about 540 m/s on Homalite 100, and about 450 m/s on Araldite B. From the photographs, crack speeds immediately before and after bifurcation are obtained, and it is found that discontinuous change of crack speed does not exist at the moment of bifurcation in the case of Homalite 100, but exists in the case of Araldite B. From the photographs, crack opening displacement (COD) is also measured along the cracks as a function of distance r from the crack tips. The measurement results show that the CODs are proportional to √r before bifurcation. After bifurcation, the CODs of mother cracks are proportional to √r, though the CODs of branch cracks are not always proportional to √r. The energy release rate is obtained from the measured CODs, and it is found that energy release rate is continuous at bifurcation point in both cases of Homalite 100 and Araldite B. Energy flux that shows the energy flow toward a crack tip is also obtained.     

Experiments on crack identification in cantilever beams

In this paper we present an experimental investigation of the identification of crack location and size. By providing the first three natural frequencies through vibration measurements, curves of crack equivalent stiffness versus crack location are plotted, and the intersection of the three curves predicts the crack location and size. In the experiments, the cracked specimens were made using a wire-cut electrical discharge machine, and the cantilever beams were excited next to the free end by means of an impulse force hammer. In order to obtain the accurate natural frequencies from the transient signal measured, the method of zoom fast Fourier transform is adopted to improve frequency resolution. From experimental results, it is observed that the identification errors of crack location and size are less than 2% and 4%, respectively. The effectiveness of crack identification through vibration measurements is verified.     

Ultrasonic sizing of embedded vertical cracks in ferritic steel welds

Embedded vertical cracks are very often encountered in structural components necessitating reliable procedure for determining their size. Specialized ultrasonic techniques are required for estimation of through thickness height of these cracks [Mater. Eval. (1970) 28; J. NDI 26 (1977) 320; NND (1975) 146]. The conventional pulse-echo technique has serious limitations for quantitative estimation of vertical cracks because it is subject to variations due to reflectivity, coupling factors etc. The present study aims to improve the accuracy of embedded crack height measurement based on transit time rather than the reflected amplitude [NDT Int. (1982) 315; J. JSNDT 27 (1978) 118]. A unique method has been developed for generating known sizes of fatigue cracks ranging from 2.5 to 18 mm in height and 15 mm in length. These cracks were generated in rectangular steel specimens and subsequently embedded at predetermined locations using welding. These welded specimens containing the embedded vertical cracks were examined using ultrasonic techniques based on transit time like (a) flaw tip echo method and (b) mode conversion method. Using both the techniques excellent agreement was obtained between the estimated crack height and actual crack height with an accuracy of ±1 mm.     

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.     

嵌岩桩桩端承载力探讨

以往嵌岩桩的试验资料表明,当上覆岩土层较薄、嵌岩桩的入岩深度较小时,嵌岩桩的极限端承力对桩的承载能力起主导作用。为求出桩端极限承载力的表达式,本文将嵌岩段以上岩土层自重视为作用在基岩半平面上的均布荷载,把桩端合力分解为相对简单的荷载作用于基岩平面上,利用滑移线理论并通过叠加的方法分析出桩端基岩破坏时的最危险点。然后根据Mindlin课题求出该点的应力,运用Griffith准则的Murrell推广导出桩端极限承载力的公式,并指出了此公式的适用范围。最后结合工程实际进行了对比分析,表明本文方法具有一定的实用性。     

Molecular dynamics simulation of crack-tip processes in copper

The crack tip processes in copper under mode II loading have been simulated by a molecular dynamics method. The nucleation, emission, dislocation free zone (DFZ) and pile-up of the dislocations are analyzed by using a suitable atom lattice configuration and Finnis & Sinclair potential. The simulated results show that the dislocation emitted always exhibits a dissociated fashion. The stress intensity factor for dislocation nucleation, DFZ and dissociated width of partial dislocations are strongly dependent on the loading rate. The stress distributions are in agreement with the elasticity solution before the dislocation emission, but are not in agreement after the emission. The dislocation can move at subsonic wave speed (less than the shear wave speed) or at transonic speed (greater than the shear wave speed but less than the longitudinal wave speed), but at the longitudinal wave speed the atom lattice breaks down. The project supported by the National Natural Science Foundation of China     

Post-buckling of cantilever columns having variable cross-section under a combined load

Post-buckling of a cantilever column is examined under a combined load consisting of a tip-concentrated load and a distributed axial load, through dynamic formulation. The formulation of the problem is based on the moment–curvature relationship. The two-point boundary value problem described by the governing equations is dependent on the frequency parameter and the two load parameters. The buckling loads are those loads at which the eigencurve, namely, the load versus frequency curve of the column meets the load axis. A simple and reliable iterative procedure to convert the two-point boundary value problem into an initial value problem is followed and solved the non-linear differential equations utilizing a fourth-order Runge–Kutta integration scheme. To demonstrate the potentiality of the adopted numerical scheme, linear vibration frequencies of truncated, tapered cantilever wedges and cones are determined and compared with the published analytical and test results. Buckling and post-buckling loads of a simply supported stepped column are obtained and compared with the published test results. The loads and deflections of non-uniform cantilever columns are obtained for various slopes at the tip. The interaction of load parameters for a free–free truncated conical column has also been examined. The numerical results indicate that the path represented by the two load parameters turns out to be nearly a straight line.