带凸缘弯曲梁应力集中区内表面裂纹应力强度因子的研究

用三维光弹性冻结应力实验技术与修正的多点超定法相结合研究了带凸缘弯曲梁应力集中区内表面裂纹的应力强度因子。分析了不同过渡圆弧的应力集中对两种表面裂纹(半圆形表面裂纹与前缘直线表面裂纹)的影响。用实验方法得到了应力强度因子放大系数的数值.结果表明,应力集中对浅裂纹的影响是更大的,是不可忽视的,但放大系数随表面裂纹的几何形状变化很小。这些对管节点的断裂力学评估提供了有价值的实验依据.     

The interaction of doubly periodic cracks

In this paper, an extremely accurate and efficient method for computing the interaction of a set of or multiple sets of general doubly periodic cracks has been presented on the basis of superposition principle, pseudo-traction method, and isolating analysis technique. A great number of typical examples are given in this paper. The stress intensity factors (SIF), the minimum strain energy density factors (SED) of crack tip and the critical stress (CRS) of crack growth are calculated with the accuracy of six significant digits for the rectangularly distributed periodic cracks and five significant digits for the general doubly periodic cracks. The relation of the interaction effect of the double periodic cracks with the periods and the ratio of crack length to crack spacing is analyzed. Also in this paper, the key technique problems for this method are discussed.     

双周期平行四边形排列裂纹反平面问题的有限元分析

提出了双周期平行四边形排列裂纹反平面问题的有限元方法,通过对单位胞元引入周期边界条件,在裂纹尖端采用奇异单元,解决了有限元分析这类问题的效率和精度问题.利用Ansys软件计算,在各种有解析解对照的情形下,应力强度因子的相对误差都在0.2%以内.与现有通常限于对称阵列的双周期裂纹的研究相比,本文发展的方法适用于一般的非对称平行四边形裂纹阵列.算例揭示了行向裂纹间的相互干涉放大应力强度因子,而叠向裂纹间的作用相互屏蔽.对于平行四边形阵列的情形,这两种相反的干涉效应使应力强度因子与裂纹错动参数间呈现非单调依赖关系.     

不同倾角预制裂纹缺陷与运动裂纹的相互作用

为了研究运动裂纹与不同倾角预制裂纹缺陷之间的作用机制,采用数字激光动态焦散线方法对含不同倾角预制裂纹缺陷的三点弯曲梁进行冲击实验。研究结果表明,在冲击载荷作用下,预制裂纹缺陷尖端均产生了次生裂纹;当运动主裂纹与预制裂纹缺陷贯通时,次生裂纹不会立即起裂,而是经过0~10 μs的能量积蓄后,次生裂纹才起裂;运动主裂纹应力强度因子峰值与次裂纹起裂时的应力强度因子值均随预制裂纹倾角的增大而增大。     

Interaction between collinear periodic cracks in an infinite piezoelectric body

The problem of collinear periodic cracks in an infinite piezoelectric body is studied. Effect of saturation strips at the crack-tips is taken into account. By means of the Stroh formalism and the conformal mapping technique, the general periodic solutions for collinear cracks are obtained. The stress intensity factors and the size of saturation strips are derived analytically, and their dependencies on the ratio of the periodicity on the half-length of the crack are analyzed in detail. Numerical results show the following two facts. (1) When h/l>4.0, the stress intensity factors become almost identical to those of a single crack in an infinite piezoelectric body. This indicates that the interaction between cracks can be ignored in establishing the criterion for the crack initiation in this case. (2) The speed of the saturation strip size of periodic cracks approaching that of a single crack depends on the electric load applied at infinity. In general, a large electric load at infinity is associated with a slow approaching speed.     

Isotropic/orthotropic two-layered strips containing cracks terminating at and going through an interface

Crack problems for isotropic/orthotropic two-layered strips have been investigated. A system of two singular integral equations can be derived by using Fourier integral transformation and boundary conditions of crack problems. After stress singularities at crack tips or other special points are determined for internal and edge cracks, and for cracks terminating at and going through the interface, the system of singular integral equations is solved numerically by Gauss-Jacobi or Gauss-Chebyshev integration formulas for stress intensity factors at the tips and other singular points of cracks. Finally, possible crack growth behavior for cracks approaching and going through the interface is discussed.     

Nonlocal theory solution of two collinear cracks in the functionally graded materials

In this paper, the interaction of two collinear cracks in functionally graded materials subjected to a uniform anti-plane shear loading is investigated by means of nonlocal theory. The traditional concepts of the nonlocal theory are extended to solve the fracture problem of functionally graded materials. To make the analysis tractable, it is assumed that the shear modulus varies exponentially with the coordinate vertical to the crack. By use of the Fourier transform, the problem can be solved with the help of a pair of triple integral equations, in which the unknown variable is the displacement on the crack surfaces. To solve the triple integral equations, the displacement on the crack surfaces is expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularity is present near the crack tips. The nonlocal elastic solutions yield a finite hoop stress at the crack tip, thus allowing us to use the maximum stress as a fracture criterion in functionally graded materials. The magnitude of the finite stress field depends on the crack length, the distance between two cracks, the parameter describing the functionally graded materials and the lattice parameter of the materials.     

基于权函数法的核主泵飞轮应力强度因子研究

采用权函数法确定了含裂纹飞轮在离心力和接触压力作用下的应力强度因子,计算了在同步转速工况下裂纹尖端应力强度因子的值,并与解析法和有限元法计算结果进行了对比。结果表明:权函数法与解析法的误差在3%以内,与有限元法的误差在1%以内,验证了权函数法计算应力强度因子的准确性高;在结构不变的情况下,权函数法可以快速求解不同载荷条件、不同长度裂纹的应力强度因子。通过控制变量法研究了不同参数对应力强度因子的影响,结果表明,飞轮裂纹尖端总应力强度因子随着裂纹尺寸、旋转转速、飞轮尺寸外径与内径比值的增大而增大。     

Dynamic behavior of two collinear interface cracks between two dissimilar functionally graded piezoelectric/piezomagnetic material strips

The dynamic interaction of two collinear interface cracks between two dissimilar functionally graded piezoelectric/piezomagnetic material strips subjected to the anti-plane shear harmonic stress waves was investigated. By using the Fourier transform, the problem can be solved with the help of a pair of triple integral equations in which the unknown variable is jump of displacement across the crack surfaces. These equations are solved using the Schmidt method. Numerical examples are provided to show the effect of the functionally graded parameter, the circular frequency of the incident waves and the thickness of the strip upon stress, electric displacement and magnetic flux intensity factors of cracks.     

THE CONSTITUTIVE RELATION OF CRACK-WEAKENED ROCK MASSES UNDER AXIAL-DIMENSIONAL UNLOADING

An accurate and efficient numerical method for solving the crack-crack interaction problem is presented. The method is mainly by means of the dislocation model, stress superposition principle and Chebyshev polynomial expansion of the pseudo-traction. This method can be applied to compute the stress intensity factors of multiple kinked cracks and multiple rows of periodic cracks as well as the overall strains of rock masses containing multiple kinked cracks under complex loads. Many complex computational examples are given. The dependence of the crack-crack interaction on the crack configuration, the geometrical and physical parameters, and loads pattern, is investigated. By comparison with numerical results under confining pressure unloading, it is shown that the crack-crack interaction under axial-dimensional unloading is weaker than those under confining pressure unloading. Numerical results for single faults and crossed faults show that the single faults are more unstable than the crossed faults. It is found from numerical results for different crack lengths and different crack spacing that the interaction among kinked cracks decreases with an increase in length of the kinked cracks and the crack spacing under axial-dimensional unloading.     

界面下主微裂纹间干涉的研究

涉及两相正交各向异性体界面干涉问题的研究,多裂纹问题被分解为只含单裂纹的子问题,利用位错理论和裂面应力自由条件,列出一组可数值求解位错密度函数的奇异积分方程,从耐 注得应力强度因子。     

Elastic interaction of interfacial spherical-cap cracks in hollow particle filled composites

This work analyzes the elastic interaction between two spherical-cap cracks present along the outer surface of a hollow particle embedded in a dissimilar medium under remote uniaxial tensile loading. A semi-analytical approach based on an enriched Galerkin method is adopted to determine stress and deformation fields as functions of particle wall thickness and cracks’ configuration. The present analysis is limited to multiple interfacial spherical-cap cracks; that is, crack propagation is restrained to the particle-matrix interface and possibility of crack kinking in the matrix is not considered. Interfacial crack growth characteristics, conditions for stable crack propagation, equal crack growth, and shielding are established through energy release rate analysis. The study is relevant to the analysis of tensile and flexural failure of syntactic foams used in marine and aerospace applications. Results specialized to glass-vinyl ester syntactic foams demonstrate that particle wall thickness can be used to control crack stability and growth characteristics as well as tailoring the magnitude of the shielding phenomenon. Predictions are compared to finite element findings for validation and to results for penny-shaped cracks to elucidate the role of crack curvature.     

The periodic array of interface cracks and their interaction

A periodic array of interface cracks is considered as a simple model for studying interface crack interaction. The novel feature of the solution is the elimination of the oscillatory singularities and of the material overlapping which is accomplished by allowing frictionless contact zones at the crack tips. The results are obtained for a combination of shear and tension-compression loads applied at infinity. It is shown that the stress intnsity factors become smaller with decreasing distance between the cracks. This indicates that the spreading of an interface crack is inhibited by the presence of adjacent cracks.     

Uniformly moving screw dislocation interacting with interface cracks in anisotropic bimaterials

This paper attempts to investigate the problem for the interaction between a uniformly subsonic moving screw dislocation and interface cracks in two dissimilar anisotropic materials. Using Riemann–Schwarz’s symmetry principle integrated with the analysis singularity of complex functions, we present the general elastic solutions of this problem and the closed form solutions for interface containing one and two cracks. The expressions of stress intensity factors at the crack tips and image force acting on moving dislocation are derived explicitly. The results show that the stress intensity factors at the crack tips decrease with increasing velocity of dislocation, and larger dislocation velocity leads to the equilibrium position of dislocation leaving from crack tips. The presented solutions contain previously known results as the special cases.     

SELF-SIMILAR CRACK EXPANSION METHOD FOR THE STRESS INTENSITY FACTOR ANALYSIS

The Self-Similar Crack Expansion (SSCE) method is proposed to evaluate stress intensi-ty factors at crack tips, whereby stress intensity factors of a crack can be determined by the crackopening displacement over the crack, not just by the local displacement around the crack tip. The crackexpansion rate is estimated by taking advantage of the crack self-similarity. Therefore, the accuracy ofthe calculation is improved. The singular integrals on crack tip elements are also analyzed and are pre-cisely evaluated in terms of a special integral analysis. Combination of these two techniques greatly in-creases the accuracy in estimating the stress distribution around the crack tip. A variety of two-dimen-sional cracks, such as subsurface cracks, edge cracks, and their interactions are calculated in terms ofthe self-similar expansion rate. Solutions are satisfied with errors less than 0.5% as compared with theanalytical solutions. Based on the calculations of the crack interactions, a theory for crack interactionsis proposed such that for a group of aligned cracks the summation of the square of SIFs at the right tipsof cracks is always equal to that at the left tips of cracks. This theory was proved by the mehtod ofSelf-Similar Crack Expansion in this paper.     

Interaction of elastic waves with a periodic array of collinear inplane cracks

A boundary integral equation method is applied to the study of the interaction of plane elastic waves with a periodic array of collinear inplane cracks. Numerical results are presented for the dynamic stress intensity factors. The effects of the wave type, wave frequency, wave incidence angle, and crack spacing on the dynamic stress intensity factors are analyzed in detail. The project supported by the Committee of Science and Technology of Shanghai and Tongji University     

Thermal shock analyses for a strip containing an infinite row of periodically distributed cracks normal to its edge

The transient thermal stress problem of a semi-infinite plate containing an infinite row of periodically distributed cracks normal to its edge is investigated in this paper. The elastic medium is assumed to be cooled suddenly on the crack-containing edge. By the superposition principle, the formulation leads to a mixed boundary value problem, with the negating tractions arisen from the thermal stresses for a crack-free semi-infinite plate. The resulting singular integral equation is solved numerically. The effects on the stress intensity factors due to the presence of periodically distributed cracks in a semi-infinite plate are illustrated. For both the edge crack and the embedded crack arrays, the stress intensity factors increase, due to the reduction of the shielding effect, as the stacking cracks are more separated. For the case of embedded crack array, one has the further conclusion that the stress intensity factors decline as the crack array shifts from the plate edge.     

Numerical analysis of fast fracture and crack arrest in a sheet of elastic perfectly-plastic material

The dynamic fields for acceleration, deceleration and arrest of a crack tip have been investigated numerically. We consider cracks which start to extend rapidly under brittle conditions. The crack-tips then enter regions of elasto-plastic constitutive behavior and they are subsequently arrested. Results have been obtained for a symmetrically expanding central crack and for an edge crack, both in thin sheets. The elasto-plastic behavior has been described by J₂ flow theory, with the von Mises yield criterion and a bilinear relation between effective stress and effective strain. Numerical results are presented for stress and strain components at a short distance ahead of the propagating and arresting crack tips.     

Transient dynamic responses of a cracked solid subjected to in-plane loadings

Transient dynamic responses of an elastic cracked solid subjected to in-plane surface loadings are investigated in this study. Two vertical cracks, a surface-breaking crack and a sub-surface crack, are considered. The frequency responses of the plane strain problem are calculated by the computational mechanics combining the finite element method with the boundary integral equation. The finite element method is used for the near-field enclosing the crack, while the boundary integral equation is applied for the far-field to satisfy the Sommerfeld radiation condition. The transient responses are then obtained using fast Fourier transform. Surface displacements, crack opening displacements, and dynamic stress intensity factors are presented to show the significant effects of the cracks. The interaction between the elastic waves and the cracks as well as the mode conversion phenomena can be observed and understood through the numerical simulations.     

Closed-form weight function for edge crack problems

An accurate crack surface opening displacement formula for edge cracks in finite bodies has been derived, and was subsequently used to develop wide-range closed-form weight functions for edge cracks. The analytical weight function is general and can be easily used to determine accurate stress intensity factors and crack surface opening displacements for arbitrarily loaded edge cracks with high efficiency. Examples have been given for problems of edge crack(s) with a number of typical load conditions for illustration.