近置多裂纹相互作用的渐近分析方法

考虑到近置裂纹的强相互作用,提出了一种多裂纹相互作用的渐近分析方法. 经典Kachanov 方法将裂纹表面伪面力分解为两部分：均匀分布部分和非均匀部分,并假设裂纹的相互作用 仅由均匀部分引起,而忽略非均匀部分的影响. 该假设大大简化了分析过程,而且当裂纹间 距不是很小时,有很好的精度. 但当裂纹非常靠近或者沿主荷载方向重叠时,由于裂纹尖端 进入了其它裂纹的应力强化区或者应力屏蔽区,强相互作用使得该假设不再合理,从而带来 较大的误差. 为了提高分析近置多裂纹问题的精度,将裂纹表面伪面力分解为抛物线型分布 部分及高阶部分,考虑抛物线型分布张力对其他裂纹的影响,同样忽略高阶部分的影响. 通 过对三共线裂纹及两平行偏置裂纹两个实例的分析,验证了对于近置裂纹,新渐近方法具有 良好的精度.     

Interactions of Penny-shaped Cracks in Three-dimensional Solids

The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is developed in which the opening and the sliding displacements on each crack surface are taken as the basic unknown functions. The basic unknown functions can be expanded in series of Legendre polynomials with unknown coefficients. Based on superposition technique, a set of governing equations for the unknown coefficients are formulated from the traction free conditions on each crack surface. The boundary collocation procedure and the average method for crack-surface tractions are used for solving the governing equations. The solution can be obtained for quite closely located cracks. Numerical examples are given for several crack problems. By comparing the present results with other existing results, one can conclude that the present method provides a direct and efficient approach to deal with three-dimensional solids containing multiple cracks.The English text was polished by Keren Wang     

Relaxation procedures for an iterative MFS algorithm for the stable reconstruction of elastic fields from Cauchy data in two-dimensional isotropic linear elasticity

We investigate two numerical procedures for the Cauchy problem in linear elasticity, involving the relaxation of either the given boundary displacements (Dirichlet data) or the prescribed boundary tractions (Neumann data) on the over-specified boundary, in the alternating iterative algorithm of Kozlov et al. (1991). The two mixed direct (well-posed) problems associated with each iteration are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method, while the optimal value of the regularization parameter is chosen via the generalized cross-validation (GCV) criterion. An efficient regularizing stopping criterion which ceases the iterative procedure at the point where the accumulation of noise becomes dominant and the errors in predicting the exact solutions increase, is also presented. The MFS-based iterative algorithms with relaxation are tested for Cauchy problems for isotropic linear elastic materials in various geometries to confirm the numerical convergence, stability, accuracy and computational efficiency of the proposed method.     

Longitudinal vibrations of arbitrary non-uniform rods

Free and steady-state forced longitudinal vibrations of non-uniform rods are investigated by an iteration method, which results in a series solution. The series obtained are convergent and linearly independent. Its convergence is verified by convergence tests, its linear independence confirmed by the nonzero value of the corresponding Wronski determinant. Then, the solution obtained is an exact one reducible to a classical solution for the case of uniform rods. In order to verify the method, two examples are presented as an application of the proposed method. The results obtained are equivalent to the method in literature. In contrast to the proposed method capable of dealing with arbitrary non-uniform rods in principle, the method in literature is confined to work on special cases.     

Anti-plane shear of rectangular region with two cracks

An approach for obtaining the solution of multiple cracks in a rectangular region is presented. A pair of concentrated shear force applied to the crack surface is used in conjunction with the method of superposition for generating the results for cracks whose surfaces are free from tractions. A system of Fredholm integral equations can always be obtained to yield numerical results.     

ON A CLASS OF METHOD FOR SOLVING PROBLEMS WITH RANDOM BOUNDARY NOTCHES AND/OR CRACKS—(Ⅲ) COMPUTATIONS FOR BOUNDARY CRACKS

This paper continues the discussions to a class of method for solving problems withrandom boundary notches and for cracks in refs.[1] and [2].Using the method developed in[1].[2]with important modifications about inclusion of singularities in the formulation. wearrive at a very effective computational process for problems with random boundarycracks. Actual computations for boundary cracks with or without apptied tractions in theirsurfaces. Show that the present method is quite workable for the problems consideredwithin proper range of characteristic parameters. The results obtained here extend thecontents of “Handbook of Stress Intensity Factors” given by G. C. Sih.     

Analysis of bi-material interface cracks with complex weighting functions and non-standard quadrature

A boundary element formulation is developed to determine the complex stress intensity factors associated with cracks on the interface between dissimilar materials. This represents an extension of the methodology developed previously by the authors for determination of free-edge generalized stress intensity factors on bi-material interfaces, which employs displacements and weighted tractions as primary variables. However, in the present work, the characteristic oscillating stress singularity is addressed through the introduction of complex weighting functions for both displacements and tractions, along with corresponding non-standard numerical quadrature formulas. As a result, this boundary-only approach provides extremely accurate mesh-insensitive solutions for a range of two-dimensional interface crack problems. A number of computational examples are considered to assess the performance of the method in comparison with analytical solutions and previous work on the subject. As a final application, the method is applied to study the scaling behavior of epoxy–metal butt joints.     

Stress distribution in a two-dimensional infinite anisotropic medium with collinear cracks

The plane problem of an anisotropic material with cracks, whose surfaces are subject to surface tractions of a general kind, is studied. The medium considered if of infinite extent and the cracks are located on a single line. The Fourier transform method is employed to derive the stress and displacement components at an arbitrary point of the medium in terms of the dislocation densities and the stress discontinuities on the crack line.These formulae for stress and displacement components involve the roots of a quartic equation whose coefficients are the material constants. The cases of different roots and pairwise coincident roots are examined separately. An orthotropic medium is an important example for the case of different roots while an isotropic medium is that for the case of pairwise coincident roots. These examples are discussed in detail.As an illustration of the use of these formulae the problem of a single crack in an infinite anisotropic medium is examined in detail.Work supported in part by a grant from Council of Scientific and Industrial Research, New Delhi, India.     

平面弹性介质多孔多裂纹问题的一种数值方法

提出了平面弹性介质中多孔洞多裂纹相互作用问题的一种数值计算方 法. 通过把适于单一裂纹的Bueckner原理扩充到含有多孔洞多裂纹的一般体系，将原问题 分解为承受远处载荷不含裂纹不含孔洞的均匀问题，和在远处不承受载荷但在裂纹面上和孔 洞表面上承受面力的多孔洞多裂纹问题. 于是，以应力强度因子作为参量的问题可以通过考 虑后者(多孔洞多裂纹问题)来解决，而利用提出的杂交位移不连续法，这种多孔 洞多裂纹问题是容易数值求解的. 算例说明该数值方法对分析平面弹性介质中多孔洞多裂纹 相互作用的问题既简单又有效.     

Interacting cracks and ellipsoidal inhomogeneities by the equivalent inclusion method

Based on the Eshelby's equivalent inclusion method (EIM) and Hill's theorem on discontinuities of elastic fields across the interfaces, a theory for the determination of the stress intensity factors (SIFs) of arbitrarily oriented interacting cracks under non-uniform far-field applied stress (strain) is developed. As shown in this investigation the EIM proposed by Moschovidis and Mura can be extended for treatment of such problems, but their formulations are quite cumbersome and computationally inefficient. An alternative analytical approach is proposed that is computationally more efficient, and unlike the method of Moschovidis and Mura can easily handle complex problems of interacting inhomogeneities and cracks. It is seen that as the interaction between the inhomogeneities becomes stronger, this method yields results that are closer to the solutions reported in the literature than the solutions obtained using the extended EIM of Moschovidis and Mura, which is developed herein. Problems involving combinations of interacting elliptic and penny shape cracks and inhomogeneities are excellent candidates for demonstration of the accuracy and robustness of the present theory, for which the previous EIM produces less accurate results. Due to the limitations imposed on the existing methods, every reported treatment has been tailored for a certain category of problems, and only uniform far-field loadings have been remedied. In contrast, the present theory is more general than the previously reported theories and it encompasses interacting cracks having a variety of geometries subjected to non-uniform far-field applied stress (strain); moreover, it is applicable to modes, I, II, III, and mixed mode fracture.     

BIEM analysis of dynamically loaded anti-plane cracks in graded piezoelectric finite solids

Anti-plane cracks in finite functionally graded piezoelectric solids under time-harmonic loading are studied via a non-hypersingular traction based boundary integral equation method (BIEM). The formulation allows for a quadratic variation of the material properties in two directions. The boundary integral equation (BIE) system is treated by using the frequency dependent fundamental solution based on Radon transforms. Its numerical solution provides the displacements and tractions on the external boundary as well as the crack opening displacements from which the mechanical stress intensity factor (SIF) and the electrical displacement intensity factor (EDIF) are determined. Several examples for single and multiple straight and curved cracks demonstrate the applicability of the method and show the influence of the different system parameters.     

孔洞裂纹问题的一种数值方法

研究孔洞与裂纹的相互作用问题,通过把适于单一裂纹的Bueckner原理扩充到含有多孔洞多裂纹的一般体系,将原问题分解为承受远处载荷不含裂纹不含孔洞的均匀问题,和在远处不承受载荷但在裂纹面上和孔洞表面上承受面力的多孔洞多裂纹问题.于是,以应力强度因子作为参量的问题可以通过考虑后者来解决,而利用笔者提出的杂交位移不连续法,这种多孔洞多裂纹问题是容易数值求解的.算例说明该数值方法对分析平面弹性介质中孔洞与裂纹的相互作用既简单又有效.     

Size effect on cracked functional composite micro-plates by an XIGA-based effective approach

Failure of structures and their components is one of major problems in engineering. Studies on mechanical behavior of functionally graded (FG) microplates with defects or cracks by effective numerical methods are rarely reported in literature. In this paper, an effective numerical model is derived based on extended isogeometric analysis (XIGA) for assessment of vibration and buckling of FG microplates with cracks. Based on the modified couple stress theory, the non-classical theory of Reissner–Mindlin plate is extended to capture microstructure, and thus, the size effect. In such theory, possessing C¹-continuity is straightforward with the high-order continuity of non-uniform rational B-spline. Due to the use of enrichments in XIGA, crack geometry is independent of the computational mesh. Numerical examples are performed to illustrate the effects of microplate aspect ratio, crack length, internal material length scale parameter, material distribution, and boundary condition on the mechanical responses of cracked FG microplates. The obtained results are compared with reference solutions and that shows that the frequency and buckling loads increases with decreasing the size of FG microplates and crack length. The convergence of the present method is also studied.     

主裂纹与微裂纹相互作用问题的有效数值计算方法

提出了平面弹性介质中主裂纹与微裂纹相互作用问题的有效数值计算 方法. 通过把适于单一裂纹的Bueckner原理扩充到含有多裂纹的一般体系，将原问题分解 为承受远处载荷不含裂纹的均匀问题，和在远处不承受载荷但在裂纹面上承受面力的多裂纹 问题. 于是，以应力强度因子作为参量的问题可以通过考虑后者(多裂纹问题)来解决，而 利用提出的杂交位移不连续法，这种多裂纹问题是容易数值求解的. 列举 Cai和 Faber为评价主裂纹与微裂纹相互作用问题的近似方法而列举的算例，说明 该数值方法对分析平面弹性介质中主裂纹与微裂纹相互作用问题既简单又非常有效.     

Application of the extended Kantorovich method to the bending of clamped cylindrical panels

Applicability and performance of the extended Kantorovich method (EKM) to obtain highly accurate approximate closed form solution for bending analysis of a cylindrical panel is studied. Fully clamped panel subjected to both uniform and non-uniform loadings is considered. Based on the Love–Kirchhoff first approximation for thin shallow cylindrical panels, the governing equations of the problem in terms of three displacement components include a system of two second order and one forth order partial differential equations. The governing PDE system is converted to a double set of ODE systems by assuming separable functions for displacements together with utilization of the extended Kantorovich method. The resulted ODE systems are solved iteratively. In each iteration, exact closed form solutions are presented for both ODE systems. Rapid convergence and high accuracy of the method is shown for various examples. Both displacement and stress predictions show close agreement with other analytical and finite element analysis.     

Electro-elastostatic analysis of multiple cracks in an infinitely long piezoelectric strip: A hypersingular integral approach

The problem of an arbitrary number of arbitrarily oriented straight cracks in an infinitely long piezoelectric strip is considered here. The cracks are acted by suitably prescribed internal tractions and are assumed to be either electrically impermeable or permeable. A Green's function which satisfies the conditions on the parallel edges of the strip is derived using a Fourier transform technique and applied to formulate the electroelastic crack problem in terms of a system of hypersingular integral equations. Once the hypersingular integral equations are solved, quantities of practical interest, such as the crack tip stress and electric displacement intensity factors, can be easily computed. Some specific cases of the problem are examined.     

含基体横向损伤的黏弹性板的蠕变后屈曲分析

基于Schapery三维黏弹性损伤本构关系,引入沈为和Kachanov损伤演化方程,建立了基体横向损伤的纤维单一方向铺设黏弹性板的损伤模型;应用von Karman板理论,导出了考虑损伤效应的具初始挠度的纤维单一方向铺设黏弹性矩形板的非线性压屈平衡方程. 对未知变量在空间上采用差分法离散,时间上采用增量算法和Newton-Cotes积分法离散,控制方程被迭代求解. 算例中讨论了损伤以及有关参数对黏弹性复合材料板后屈曲行为的影响,且与已有文献的结果进行了比较. 数值结果表明：随着外载荷或者初始挠度的增大,板后屈曲趋于稳定时的挠度就愈大,损伤的影响愈明显;而随着长宽比的增大,板后屈曲趋于稳定时的挠度愈小,损伤的影响却随之增大.     

带双裂纹的椭圆孔口问题的应力分析

利用复变方法,通过构造新的保角映射,研究了带双裂纹椭圆孔口的平面弹性问 题,得到了I型与II型裂纹问题应力强度因子的解析解.在极限情形下, 不仅可以还原经典的Griffith裂纹的结果,而且可模拟出十字裂纹和带双裂纹的圆 形孔口问题.     

三维裂纹分析的主值型面力边界积分方程法

给出了一组只包含Ｃａｕｃｈｙ主值积分、不含有强奇异积分的三维静动力边界积分方程及其应用于裂纹问题的具体列式，并给出了几何轴对称问题的相应半解析边界元求解方法，将三维问题降阶为一维数值问题．文中分析了无限、半无限介质中圆裂纹、平行圆裂纹系、球面裂纹等在静载及应力波作用下的静力或瞬态动力响应问题，求得了相应的应力强度因子．     

Some new aspects of boundary conditions at cracks in piezoelectrics

Fracture mechanical investigations of piezoelectric materials as components of smart structures have become popular in the last 30?years. In the early years of research, boundary conditions at crack faces have been adopted from pure mechanical systems under the assumption that boundaries were traction free. From the electrostatic point of view, cracks have been assumed to be either free of charge or fully permeable. Later, limitedly permeable crack boundary conditions have become popular among the community, nevertheless still assuming traction-free crack faces. Recently, the theoretical framework has been extended to include electrostatically induced mechanical tractions in crack models yielding a significant crack closure effect. However, these models are still simple, neglecting, e.g., the piezoelectric field coupling. In this work, we present an extended model for crack surface tractions yielding some interesting effects. In particular, the orientation of the electrical field with respect to the poling axis becomes important. Furthermore, applying a collinear stress parallel to the crack faces influences the Mode-I stress intensity factor and a Mode-II shear loading couples to the Mode-I SIF.