有限板共线多孔MSD疲劳裂纹扩展有限元模拟

介绍了二维断裂分析有限元软件FRANC2D/L在疲劳裂纹扩展模拟方面的基本步骤．利用该软件对有限板中心孔边对称裂纹的疲劳裂纹扩展进行模拟计算,对比模拟结果和试验数据,发现两者吻合良好,证明了利用该方法模拟疲劳裂纹扩展的可靠性．将FRANC2D/L应用到有限板共线多孔MSD疲劳裂纹扩展的有限元模拟上,得到了各孔边裂纹的长度和疲劳扩展寿命之间的关系曲线．模拟计算结果表明,在相同条件下,有限板中心孔边对称裂纹的裂纹扩展寿命要远远高于MSD结构中中心孔边裂纹的疲劳扩展寿命;由于MSD结构中影响各孔边裂纹的因素有所差异,各条裂纹的疲劳扩展寿命也会有所不同．另外,还给出了不含主裂纹的MSD和含主裂纹的MSD两种情况下的疲劳裂纹扩展历程,通过比较得知,含主裂纹的MSD结构更容易发生裂纹的合并和贯穿致使结构发生破坏．     

A general method for multiple crack problems in a finite plate

A novel method for the multiple crack problems in a finite plate is proposed in this paper. The basic stress functions of the solution consist of two parts. One is the Fredholm integral equation solution for the crack problem in an infinite plate, and the other is that of the weighted residual method for general plane problems. The combined stress functions are used in the analysis and the boundary conditions on the crack surfaces and the boundary are considered. After the coefficients of the functions have been determined, the stress intensity factors (SIF) at the crack tips can be calculated. Some numerical examples are given and it was observed that when the cracks are very short, the results compare very favorably with the existing results for an infinite plate. Furthermore, the influence of the boundary can be considered. This method can be used for arbitrary multiple crack problems in a finite plate.     

The effect of an elliptical inclusion on a crack

The interaction between an elliptical inclusion and a crack is analyzed by body force method. The investigated stress field is simulated by superposing the fundamental solutions for a point force applied at a point in an infinite plate containing an elliptical inclusion. Based on numerical results, effects of the inclusion shape on the crack tip stress intensity factor are discussed. It is found that for small cracks emanating from a stress-higher point on the inclusion interface the stress intensity factors are mainly determined by the stresses, occurring at the crack starting point before the crack initiation, and the inclusion root radius, besides the crack length. However, for the cracks occurring in a stress-lower region around the inclusion, it is difficult to characterize the effect of the inclusion geometry on the stress intensity factors of small cracks by the inclusion root radius alone. This revised version was published online in July 2006 with corrections to the Cover Date.     

Multiple coplanar embedded elliptical cracks in an infinite solid subject to arbitrary crack face tractions

This paper extends the earlier work by the present authors for a single embedded crack in an infinite solid and presents a solution to the problem of multiple coplanar cracks in an infinite medium. An alternating method in conjunction with an analytical solution for a single crack is used to determine the stress intensity factors for interacting multiple coplanar embedded cracks in an infinite body. The alternating method, as implemented here, leads to a highly accurate evaluation of the appropriate stress intensity factors.     

Effect of fiber damage on the overall electrical conductivity of bare carbon fiber strand

The simple method developed by Kachanov (1985) for multiple interacting cracks in homogenous medium is extended to predict complex stress intensity factor for multiple split type interface cracks. Calculations are implemented for two equal cracks and infinite row of periodic cracks at the interface between two dissimilar isotropic materials. Results for infinite row of cracks are compared against the exact analytical solution provided by Sih (1973). The approximate method leads to the results very close the exact solution for crack density up to 0.90 (relative error is less than 3.8% for real part of stress intensity factor) and material dissimilarity does not have a major influence on the error. For crack densities higher than 0.90, the influence of material dissimilarity is more evident and the error increases as material dissimilarity increases. The promising match between the approximate and exact method proves the capability of the approximate method for solving other interacting interface crack problems, such as multiple penny-shaped interface cracks, in which the solution is not obtained in the literature yet.     

含有任意多个椭圆孔复合材料层合板的应力场计算

本文应用各向异性弹性力学的复变函数理论,用多保角变换的方法,导出了含有任意多个任意位置椭圆孔的各向异性复合材料板的多复变量应力函数表达式,然后在单位圆周上进行复Fourier级数展开,用待定系数法确定应力函数的未知系数,从而计算弹性板的应力场.编制了相应的多工况运行的FORTRAN77标准化程序,进行了考题和算例分析,给出了级数的收敛状况和孔边周向应力的分布图.      

Tension of a cylindrical bar having an infinite row of circumferential cracks

In this paper, the crack problems in the case of a cylindrical bar having a circumferential crack and a cylindrical bar having an infinite row of circumferential cracks under tension are analyzed by the body force method. The stress field for a periodic array of ring forces in an infinite body is used to solve the problems. The solution is obtained by superposing the stress fields of ring forces in order to satisfy a given boundary condition. The stress intensity factors are calculated for various geometrical conditions. The obtained values of stress intensity factor of a single circumferential crack are considered to be more reliable than the results of other paper's. As the crack becomes very shallow, the stress intensity factor of a row of circumferential cracks approaches the value corresponding to that of a row of edge cracks in a semi-infinite plate under tension. As the crack becomes very deep, it approaches the values corresponding to that of a single deep circumferential crack.     

Determination of stress intensity factor solutions for cracks in finite-width functionally graded materials

A generalized method to determine the stress intensity factor equations for cracks in finite-width specimens of functionally graded materials (FGMs), based on force balance in regions ahead of the crack tip is provided. The method uses the Westergaard's stress distribution ahead of the crack in an infinite plate and is based on the requirement of isostrain deformation of layers of varying moduli ahead of the crack tip. It is shown that the modified Westergaard equation describes the normal stress distribution and the singular stress state ahead of the crack tip in a reasonably accurate manner. Based on this, closed-form analytical equations for the stress intensity factors of cracks in finite-width center cracked specimens were derived. Comparisons of the K values from the analytical equations with that obtained from FEM simulations indicate that the derived stress intensity factor equations for FGMs are reasonably accurate. For the finite-width center-cracked-tension (CCT) specimen, the errors are less than 10% for most of the crack lengths for materials with the outer layer modulus ratios varying from 0.2 to 5. The stress intensity factors were found to be sensitive to the absolute values of moduli of the layers, the modulus ratio of the outer layers as well as the nature of gradation including the increasing and the decreasing functional forms. The stress intensity factor equations are convenient for engineering estimates of stress intensity factors as well as in the experimental determinations of fracture toughness of FGMs.     

Thermal striping fatigue damage in multiple edge-cracked geometries

Thermal fatigue striping damage may be caused when incompletely mixed hot and cold fluid streams pass over the surface of a component or structure containing a defect. Stress intensity factor (SIF) fluctuations are developed in response to the surface temperature fluctuations. An existing methodology for the analysis of striping damage in geometries containing a single edge‐crack geometry is extended to such an analysis of multiple edge cracks. SIFs are calculated as functions of crack depth, when an edge‐cracked plate and semi‐infinite solid, each containing multiple cracks, are subjected to thermal striping. The effect of various restraint conditions and striping frequencies on the SIF values for a stainless steel plate is examined. The degree of conservatism is shown when an assessment of thermal fatigue striping damage is based on a single, rather than multiple, crack analysis. Accurate curve fits are developed resulting in practical weight functions for an edge‐cracked plate and semi‐infinite solid.     

A finite element alternating approach for the bending analysis of thin cracked plates

Based on the classical plate theory, the analytical solution for an infinte thin plate containing a crack subjected to arbitrary symmetric bending moments on the crack surfaces is first derived. Using this solution, an efficient and accurate finite element alternating procedure is then devised to deal with symmetric plate bending problems with single or multiple cracks. The interaction effect among cracks and the influence of the geometric boundaries on the calculation of bending stress intensity factors are also presented in detail. Several numerical examples are solved to demonstrate the validity of the approach.     

Anti-plane transient analysis of planes with multiple cracks

In this paper, the transient dynamic stress intensity factor is determined for multiple curved cracks under impact loading. The dislocation method has rarely been applied to the problems involving dynamic loading. The transient response of Volterra-type dislocation in a plane is obtained by means of the Cagniard-de Hoop method. The distributed dislocation technique is used to construct integral equations for an infinite isotropic plane weakened by cracks. These equations are of Cauchy singular type at the location of dislocation which are solved numerically to obtain the dislocation density on the faces of the cracks. The dislocation densities are employed to determine stress intensity factors for multiple smooth cracks. Numerical results are obtained to validate the formulation and illustrate its capabilities.     

Weight function,stress intensity factor and crack opening displacement solutions to periodic collinear edge hole cracks

Periodic collinear edge hole cracks and arbitrary small cracks emanating from collinear holes, which are two typical multiple site damages occurred in the aircraft structures, are studied by using the weigh function method. An explicit closed form weight function for periodic edge hole cracks in an infinite sheet is obtained and further used to calculate the stress intensity factor and crack opening displacement for various loading cases. Compared to finite element method, the present weight function is accurate and highly efficient. The interactions of the holes and cracks on the stress intensity factor and crack opening displacement are quantitatively determined by using the present weight function. An approximate weight function method is also proposed for arbitrary small cracks emanating from multiple collinear holes. This method is very useful for calculating the stress intensity factor for arbitrary small cracks.     

Mixed mode fatigue growth of curved cracks emanating from fastener holes in aircraft lap joints

The finite element alternating method is extended further for analyzing multiple arbitrarily curved cracks in an isotropic plate under plane stress loading. The required analytical solution for an arbitrarily curved crack in an infinite isotropic plate is obtained by solving the integral equations formulated by Cheung and Chen (1987a, b). With the proposed method several example problems are solved in order to check the accuracy and efficiency of the method. Curved cracks emanating from loaded fastener holes, due to mixed mode fatigue crack growth, are also analyzed. Uniform far field plane stress loading on the plate and sinusoidally distributed pin loading on the fastener hole periphery are assumed to be applied. Small cracks emanating from fastener holes are assumed as initial cracks, and the subsequent fatigue crack growth behavior is examined until long arbitrarily curved cracks are formed near the fastener holes under mixed mode loading conditions.     

Stress intensity factors for interacting cracks and complex crack configurations in linear elastic media

In this paper, an effective numerical method for analyzing interacting multiple cracks and complex crack configurations in infinite linear elastic media is presented. By extending Bueckner’s principle suited for a crack to a general system containing multiple interacting cracks, the original problem is divided into a homogeneous problem (the one without cracks) subjected to remote loads and a multiple-crack problem in an unloaded body with applied tractions on the crack surfaces. Thus, the results in terms of the stress intensity factors (SIFs) can be obtained by taking into account the latter problem, which is analyzed easily by means of the displacement discontinuity method with crack-tip elements proposed recently by the author. Test examples are included to illustrate that the method is very simple and effective for analyzing interacting multiple cracks and complex crack configurations in an infinite linear elastic media under remote uniform stresses. Specifically, analysis of perpendicular cracks under general in-plane loading is performed using the numerical approach and many numerical results are given in the form of tables.     

A relaxation technique for evaluating stress intensity factors by the finite element method

A simple two-step corrective technique is presented in this paper for evaluating stress intensity factors in crack problems. In the first step an approximate evaluation of the stress intensity factor was made by considering the cracked plate to be of infinite size. The stresses of the problem were relaxed by the stresses of the infinite body which corresponds to the approximate value of the stress intensity factor. The expected discrepancy in the value of SIF by the infinite plate approximation was corrected in the second step where the existing residual stresses are equilibrated at the cracked plate by using any of the conventional finite element techniques and the corrective value of the stress intensity factor is calculated by using an appropriate collocation formula. The method was applied to three typical plane problems of cracked plates with satisfactory results.     

Evaluation of T‐stresses in multiple crack problems of finite plate

This paper provides a solution for T‐stresses for multiple cracks in a finite plate. The results for stress intensity factors (SIFs) are also presented. The case of two cracks in a rectangular plate is taken as an example. In the problem, the crack faces are applied by some loadings, and tractions are free along edges of a rectangular plate. The whole stress field is considered as a superposition of three particular stress fields. The first and second stress fields are initiated by loadings on the first and second crack faces in an infinite plate. The third field is chosen in a polynomial form of complex potentials. After discretization, the loadings on two cracks and the undetermined coefficients in the complex potentials become the unknowns. The relevant algebraic equations are formulated. The solution of algebraic equations will lead to the results of SIFs and T‐stresses at the crack tips. Several numerical examples are presented, which were not reported previously.     

Stress intensity factors for cracks emanating from a triangular or square hole in an infinite plate by boundary elements

This paper concerns stress intensity factors of cracks emanating from a triangular or square hole in an infinite plate subjected to internal pressure calculated by means of a boundary element method, which consists of constant displacement discontinuity element presented by Crouch and Starfied and crack tip displacement discontinuity elements proposed by the author. In the boundary element implementation the left or the right crack tip displacement discontinuity element is placed locally at corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Numerical examples are included to show that the method is very efficient and accurate for calculating stress intensity factors of plane elasticity crack problems. Specifically, the numerical results of stress intensity factors of cracks emanating from a triangular or square hole in an infinite plate subjected to internal pressure are given.     

A finite element method for calculating stress intensity factors and its application to composites

A concept which allows for the development of efficient finite element techniques in the analysis of plane elastic structures containing cracks is discussed. It consists of combining a specially defined finite element in the region surrounding each crack tip with conventional CST elements describing the remaining portion of the geometry considered. For the special element a pair of displacement functions is chosen which adequately represents the singular character of the elastic solution at the crack tip. The application of this concept is illustrated through a specific numerical method developed by W. K. Wilson for the calculation of mode I stress intensity factors.

Wilson's method was coded and used to analyze an infinitely long strip under tension with a line crack perpendicular to its axis of symmetry. Circular inclusions of different material properties were assumed to be present near the tips of the crack and their effect on the mode I stress intensity factor was investigated.

It was found that more flexible inclusions increase the intensity factor while more rigid inclusions decrease it. These results are quite similar to those obtained by analytical methods in an analogous problem involving an infinite sheet, but in the case of a strip, the influence of inclusions on the intensity factor was found to be more pronounced.     

Dynamic stress intensity factor for an unbounded plate having collinear cracks

The steady-state vibration of an infinite plate with collinear cracks is considered for low frequency cyclic loading. The formulation of the mixed boundary value problem leads to a dual trigonometric series. The Schwinger's method gives an automatic perturbation scheme. The dynamic stress intensity factor is found to be higher than the corresponding static one. The inertial effect on the stress intensity factor becomes significant only when the frequency of the external load is close to that of the shear wave.     

The residual stress intensity factors for cold‐worked cracked holes: a technical note

Cold‐working of riveted holes reduces the stress intensity factor associated with cracks that may develop at the hole boundary, by creating a compressive residual stress field. The residual stress field is determined using the finite‐element method and the reduction of the stress intensity factor for different values of the interference is evaluated with the weight function method, in the case of an infinite plate made from an elastic–perfectly plastic material, and having a hole with two symmetrical cracks. Once the weight function of the structure is known, further calculation of the stress intensity factors for different loadings such as a remote uniform stress, or a point load that simulates the action of the rivet can be performed without difficulty.