A contact crack problem in an infinite plate

The problem of an infinite plate with crack of length 2a loaded by the remote tensile stress P and a pair of concentrated forces Q is discussed. The value of the force Q for the initial contact of crack face is investigated and the contact length elevated, while the Q force increases. The problem is solved assuming that the stress intensity factor vanishes at the end point of the contact portion. By the Fredholm integral equation for the multiple cracks, the reduction of stress intensity factor due to Q is found. (C) 1999 Elsevier Science Ltd. All rights reserved.

Torsional impact response of a penny-shaped interface crack in bonded materials with a graded material interlayer

In this paper, the dynamic response of a penny-shaped interface crack in bonded dissimilar homogeneous half-spaces is studied. It is assumed that the two materials are bonded together with such a inhomogeneous interlayer that makes the elastic modulus in the direction perpendicular to the crack surface is continuous throughout the space. The crack surfaces art assumed to be subjected to torsional impact loading. Laplace and Hankel integral transforms are applied combining with a dislocation density,function to reduce the mixed boundary value problem into a singular integral equation with a generalized Cauchy kernel in Laplace domain. By solving the singular integral equation numerically, and using a numerical Laplace inversion technique, the dynamic stress intensity factors art obtained. The influences of material properties and interlayer thickness on the dynamic stress intensity factor are investigated.

A rigorous analytical method for doubly periodic cylindrical inclusions under longitudinal shear and its application

An infinite elastic solid containing a doubly periodic parallelogrammic array of cylindrical inclusions under longitudinal shear is studied. A rigorous and effective analytical method for exact solution is developed by using Eshelby's equivalent inclusion concept integrated with the new results from the doubly quasi-periodic Riemann boundary value problems. Numerical results show the dependence of the stress concentrations in such heterogeneous materials on the periodic microstructure parameters. The overall longitudinal shear modulus of composites with periodic distributed fibers is also studied. Several problems of practical importance, such as those of doubly periodic holes or rigid inclusions, singly periodic inclusions and single inclusion, are solved or resolved as special cases. The present method can provide benchmark results for other numerical and approximate methods. (C) 2003 Elsevier Ltd. All rights reserved.

Critical serrated flow features during nanoindentation in La-based bulk metallic glasses

La57.6Al17.5(Cu,Ni)(24.9) and La64Al14(Cu,Ni)(22) bulk metallic glasses (BMGs) were prepared by copper-mould casting method. Plastic deformation behavior of the two BMGs at various loading rates was studied by nanoindentation. The results showed that the La57.6Al17.5(Cu,Ni)(24.9) BMG with a glass transition temperature of 423 K exhibited prominent serrated flow at low loading rates, whereas less pronounced serrated flow at high rates during nanoindentation. In contrast, the La64Al14(Cu,Ni)(22) BMG with a glass transition temperature of 401 K exhibited prominent serrated flow at high loading rates. The different rate dependency of serrated flow in the two La-based BMGs is related to the different glass transition temperature, and consequently the degree of viscous flow during indentation at room temperature. A smoother flow occurs in the alloy with relatively lower glass transition temperature, due to the relaxation of stress concentration.

Three-dimensional crack problem analysis using boundary element method with finite-part integrals

A set of hypersingular integral equations of a three-dimensional finite elastic solid with an embedded planar crack subjected to arbitrary loads is derived. Then a new numerical method for these equations is proposed by using the boundary element method combined with the finite-part integral method. According to the analytical theory of the hypersingular integral equations of planar crack problems, the square root models of the displacement discontinuities in elements near the crack front are applied, and thus the stress intensity factors can be directly calculated from these. Finally, the stress intensity factor solutions to several typical planar crack problems in a finite body are evaluated.

On crack path stability in a layered material

Crack paths in an elastic layer on top of a substrate are considered. Crack growth is initiated from an edge crack in the layer. The plane of the initially straight crack forms an angle to the free surface. The load consists of a pair of forces applied at the crack mouth and parallel to the interface. Crack paths are calculated using a boundary element method. Crack growth is assumed to proceed along a path for which the mode II stress intensity factor vanishes. The inclination and the length of the initial crack are varied. The effect of two different substrates on the crack path evolution is demonstrated. A crack path initially leading perpendicularly to the interface is shown to be directionally unstable for a rigid substrate. Irrespective of its initial angle, the crack does not reach the interface, but reaches the free surface if the layer is infinitely long. At finite layer length the crack reaches the upper free surface if the initial crack inclination to the surface is small enough. For an inextendable flexible substrate, on the other hand, the crack reaches the interface if its initial inclination is large enough. For the flexible substrate an unstable path parallel with the sides of an infinitely long layer is identified. The results are compared with experimental results and discussed in view of characterisation of directionally unstable crack paths. The energy release rate for an inclined edge crack is determined analytically.

Molecular dynamics simulation of interaction of a dislocation array from a crack tip with grain boundaries

The interaction of a dislocation array emitted from a crack tip under mode II loading with asymmetric tilt grain boundaries (GBs) is analysed by the molecular dynamics method. The GBs can generally be described by planar and linear matching zones and unmatching zones. All GBs are observed to emit dislocations. The GBs migrated easily due to their planar and linear matching structure and asymmetrical type. The diffusion induced by stress concentration is found to promote the GB migration. The transmissions of dislocations are either along the matched plane or along another plane depending on tilt angle theta. Alternate processes of stress concentration and stress relaxation take place ahead of the pileup. The stress concentration can be released either by transmission of dislocations, by atom diffusion along GBs, or by migration of GBs by formation of twinning bands. The simulated results also unequivocally demonstrate two processes, i.e. asymmetrical GBs evolving into symmetrical ones and unmatching zones evolving into matching ones during the loading process.

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 effect of thermal-activation on dislocation processes at an atomistic crack-tip

The effects of thermal activation on the dislocation emission from an atomistic crack tip are discussed, Molecular dynamics simulations at different constant temperatures are carried out to investigate the thermal effects. The simulated results show that the processes of the partial dislocation generation and emission are temperature dependent. As the temperature increases, the incipient duration of the partial dislocation nucleation becomes longer, the critical stress intensity factor for partial dislocation emission is reduced and, at the same loading level, more dislocations are emitted. The dislocation velocity moving away from the crack tip and the separations of partial dislocations are apparently not temperature dependent. The simulated results also show that, as the temperature increases, the stress distribution along the crack increases slightly. Therefore stress softening at the crack tip induced by thermal activation does not exist in the present simulation. A simple model is proposed to evaluate the relation of the critical stress intensity factor versus temperature. The obtained relation is in good agreement with our molecular dynamics results.

3-dimensional transient wave response in a cracked elastic solid

The three-dimensional transient wave response problem is presented for an infinite elastic medium weakened by a plane crack of infinite length and finite width. Tractions are applied suddenly to the crack, which simulates the case of impact loading. The integral transforms are utilized to reduce the problem to a standard Fredholm integral equation in the Laplace transform variable and sequentially invert the Laplace transforms of the stress components by numerical inversion method. The dynamic mode I stress intensity factors at the crack tip are obtained and some numerical results are presented in graphical form.

Orthotropy rescaling and implications for fracture in composites

For most practically important plane elasticity problems of orthotropic materials, stresses depend on elastic constants through two nondimensional combinations. A spatial rescaling has been found to reduce the orthotropic problems to equivalent problems in materials with cubic symmetry. The latter, under favorable conditions, may be approximated by isotropic materials. Consequently, solutions for orthotropic materials can be constructed approximately from isotropic material solutions or rigorously from cubic ones. The concept is developed to gain insight into the interplay between anisotropy and finite geometry. The inherent simplicity of the solutions allows a variety of technical problems to be addressed efficiently. Included are stress concentration related cracking, effective contraction of orthotropic material specimens, crack deflection onto easy fracture planes, and surface flaw induced delamination.

Experimental relaxation behavior of chopped-strand mat glass-fiber reinforced polyester beams in 3 point bending

A dimensionless relation of the form for collating fatigue crack starting growth data is proposed in which Δkth 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.

Molecular dynamics studies of brittle fracture of SiO2

By applying for molecular dynamics (MD) simulation and Griffith fracture criterion, the brittle behavior of crack extension of mode I type is investigated. The critical stress intensity factor (SIF)K-Ic(MD) of crack extension is calculated, and the evolution of atoms near crack tip is observed. It is found that K-Ic(MD) is in good agreement with the Griffith ftacture criterion K-Ic(Griffith).

涂层-梯度层-基底结构中的热残余应力和应力强度因子研究

热喷涂构件中的热残余应力分布对其使用寿命的影响很大，在涂层和基底之间加入功能梯度材料可以降低构件中的热残余应力。本文用数值方法研究了热喷涂涂层-基底结构中孔隙率的随机分布对结构中热残余应力的影响，以及结构中热残余应力和界面裂纹的应力强度因子与功能梯度层层数之间的关系，并与实验结果相验证。论文由四部分组成。首先，介绍了热喷涂技术的发展历史和目前的主要应用、热喷涂产品中存在的问题以及研究和解决这些问题的方法。其次，研究了涂层中随机分布的孔隙对构件中热残余应力的影响，第三，分析了功能梯度层分层的层数对构件中热残余应力的影响，最后，分析了功能梯度层分层的层数对构件界面裂纹应力强度因子的影响。主要结论为：在期望和方差相同的情况下，涂层中的孔隙率呈高斯分布还是呈均匀分布，对涂层-功能梯度层-基底中的热残余应力的分布影响不大；在基底、涂层和功能梯度层的相对厚度不变的情况下，功能梯度层的层数愈多（即材料性质的过渡愈平缓），涂层-（功能梯度层）-基底结构中的热残余应力愈小；在基底、涂层和功能梯度层的相对厚度不变的情况下，功能梯度层分层的层数愈多，涂层-（功能梯度层）-基底结构的界面裂纹的应力强度因子愈小。

线性压电材料的有限元分析

本文利用线性压电学理论，编制了线性压电材料四结点等参有限元程序，进行了校核，并对PZT-5A材料压电智能元件和压电材料标准断裂试进行了计算，计算包括：①采用实际工程应用的压电智能元件尺寸，计算了元件的压电响应；并针对元件内部电极尖端区域容易引起破坏的现象，计算了该区域的奇异应力、应变场及电场。②计算了加力和加电两种情况下压电材料标准断裂试件应力强度因子影响系数F_I和电位移强度因子影响系数F_D。裂纹面边界条件采用D-P条件，试件包括紧凑拉伸标准试件和三点弯曲标准试件。③应用Lagrange乘子法将Parton裂纹边界条件加于有限元程序中，计算了上述两种标准断裂试件的F_I和F_D，计算结果与采用D-P裂纹边界条件的计算结果有很大差异。