基于全试验设计方法的含裂纹铝合金板复合材料修补参数优化

建立含穿透裂纹铝合金板复合材料单面胶接修补板条的三维有限元模型,基于位移外推法对裂纹尖端的应力强度因子(SIF)进行求解。使用全试验设计的方法对不同修补参数下修补板条的单向拉伸试验进行仿真模拟,利用二次方程描述并研究了补片长度、补片厚度及胶层弹性模量共同作用时对SIF的影响,确定了以SIF为评价指标时对修补效果影响最大的修补参数,优化了修补设计,并应用优化修补参数进行单向静拉伸试验。结果表明,当三类修补参数共同作用时,补片长度对修补效果影响最大;使用优化修补参数单面修补试验件的破坏强度比未修补板的提高了12.1%,恢复到完好板的90.5%。     

Extended Isogeometric Analysis (XIGA) of Fatigue Life in Attachment Lug

In the present paper, extended isogeometric analysis (XIGA) is used to compute the stress intensity factor (SIF) in straight lugs of Aluminum 7075-T6. XIGA uses the non-uniform rational B-spline (NURBS) functions and enrichment functions through the partition of unity. The Heaviside function is enriched to capture jump at the crack faces, while the analytical asymptotic solutions are incorporate with NURBS to perform the crack tip singularity. The XIGA-based SIF of edge-cracked plate and straight attachment lug are compared with analytical solution and extended finite element method capability available in ABAQUS. Also, crack growth and fatigue life of single crack in attachment lug are estimated and then compared with the available experimental data for two different load ratios equal to 0.1 and 0.5. The SIF calculated from XIGA are in reasonable agreement with the available data.     

Multiple penny-shaped cracks interaction in a finite body and their effect on stress intensity factor

In this paper we investigate the stress intensity factors (SIFs) of multiple penny-shaped cracks in an elastic solid cylinder under mode I (axial tension) loading. The cracks are located symmetrically and in parallel to one another in the isotropic cylinder. The fractal-like finite element method (FFEM) is employed to study the interaction of multiple cracks and to demonstrate the efficiency of the FFEM for multiple crack problems. The results show that the SIF values of the inner cracks, which are denoted as crack number 1,2,3,…,(n+1)/2 of a stack of n parallel cracks, are lower than the SIF values of a single crack by between 16% and 48%. Also, the outermost crack, that is the crack closest to the boundaries of a multiple cracked body, has the highest SIF values and is, therefore, likely to fail first.     

新的估算表面裂纹应力强度因子经验公式

该文给出了新的估算拉伸和纯弯曲载荷下表面裂纹应力强度因子的经验公式。根据疲劳裂纹扩展的数值模拟结果确定强度因子分布函数;利用按已知应力强度因子分布函数求裂纹形状及相应应力强度因子的方法计算给定尺寸的表面裂纹的应力强度因子;通过对数值结果的曲线回归得到估算表面裂纹应力强度因子经验公式。利用该公式对有限厚度和宽度平板内表面裂纹的应力强度因子进行了估算,并与已知的半椭圆形表面裂纹的应力强度因子解进行了比较。该文结果为估算表面裂纹应力强度因子提供了一种新的途径。     

A new method of crack-tip opening displacement determined based on maximum crack opening displacement

In this paper a new method is presented to determine the crack-tip opening displacement (CTOD) for the center cracked plate with uniaxial uniform tension load. The maximum crack opening displacement (MCOD) is adopted to estimate CTOD. Based on the series of calculation results by elastic–plastic finite element simulation, an explicit function expression for the CTOD versus MCOD is determined, which enables to consider the influence effects of crack geometries, plate sizes, applied loads, plane state and material properties. Hence, the presented method of CTOD determined by MCOD is suitable to any center crack finite plate of any material under uniaxial tension.     

Effect of Specimen Geometry on Stress Intensity Factors of Inclined Crack by Finite Element Method

A rectangular plate with inclined cracks of different crack lengths at different crack inclination angles under biaxial loading condition are being analyzed in mixed mode condition using finite element method (FEM) for the determination of stress intensity factors (SIFs). With increases of the width of the plate and the crack length ratio, SIF increases up to 45° of crack inclination angle and then decreases; maximum value is obtained at 45° of crack inclination angle. With the increasing value of size factor, the value of SIF starts decreasing. The accuracy of the results of the proposed method is validated by comparing with the previously obtained results by theoretical and experimental methods. The FEM results give significant result for the two-dimensional mixed mode loading conditions.     

Determination of Fracture Parameters for Multiple Cracks of Laminated Composite Finite Plate

A predictive method for estimation of stress state at zone of crack tip and assessment of remaining component lifetime depend on the stress intensity factor (SIF). This paper discusses the numerical approach for prediction of first ply failure load (F_L), progressive failure load, SIF and critical SIF for multiple cracks configurations of laminated composite finite plate using finite element method (FEM). The Hashin and Chang failure criterion are incorporated in ABAQUS using subroutine approach user defined field variables (USDFLD) for prediction of progressive fracture response of laminated composite finite plate, which is not directly available in the software. A tensile experiment on laminated composite finite plate with stress concentration is performed to validate the numerically predicted subroutine results, shows excellent agreement. The typical results are presented to examine effect of changing the crack tip distance (S), crack offset distance (H), and stacking fiber angle (θ) on F_L, and SIF .     

Modeling arbitrary crack propagation in coupled shell/solid structures with X‐FEM

In this paper, the extended finite element method (X‐FEM) formulation for the modeling of arbitrary crack propagation in coupled shell/solid structures is developed based on the large deformation continuum‐based (CB) shell theory. The main features of the new method are as follows: (1) different kinematic equations are derived for different fibers in CB shell elements, including the fibers enriched by shifted jump function or crack tip functions and the fibers cut into two segments by the crack surface or connecting with solid elements. So the crack tip can locate inside the element, and the crack surface is not necessarily perpendicular to the middle surface. (2) The enhanced CB shell element is developed to realize the seamless transition of crack propagation between shell and solid structures. (3) A revised interaction integral is used to calculate the stress intensity factor (SIF) for shells, which avoids that the auxiliary fields for cracks in Mindlin–Reissner plates cannot satisfy exactly the equilibrium equations. Several numerical examples, including the calculation of SIF for the cracked plate under uniform bending and crack propagation between solid and shell structures are presented to demonstrate the performance of the developed method. Copyright © 2015 John Wiley & Sons, Ltd.     

Approximate weight function methods using a new condition

A method is presented to derive the approximate weight functions, by using a new condition on the crack mouth, for edge cracks in the semi-infinite or finite plate under Mode I loadings. This method requires only one information, say, the reference stress intensity factor (SIF), K_r. Numerical examples show that the present method is efficient in evaluation of the SIFs for edge cracks subjected to the polynomial loadings on the crack faces.     

Double noding technique for mixed mode crack propagation studies

A simple dynamic finite element algorithm for analysing a propagating mixed mode crack tip is presented. A double noding technique, which can be easily incorporated into existing dynamic finite element codes, is used together with a corrected ? integral to extract modes I and II dynamic stress intensity factors of a propagating crack. The utility of the procedure is demonstrated by analysing test problems involving a mode I central crack propagating in a plate subjected to uniaxial tension, a stationary slanted central crack in a plate subjected to uniaxial inpact loading and an extending slanted-edge crack in plate subjected to uniaxial tension.     

Stress intensity factors for large aspect ratio surface and corner cracks at a semi-circular notch in a tension specimen

Stress intensity factor solutions for semi-elliptic surface and quarter-elliptic corner cracks emanating from a semi-circular notch in a tension specimen are presented. A threedimensional finite-element analysis in conjunction with the equivalent domain integral was used to calculate stress intensity factors (SIF). SIF solutions for surface or corner crack (crack length to depth ratio of 2) at a notch are presented for a wide range of crack sizes and notch radii. Results showed that the SIF are larger for larger crack lengths and for larger notch radii. The SIF are nearly constant all along the crack front for deep surface cracks and for all corner cracks analysed.     

Effects of CFRP bond locations on the Mode I stress intensity factor of centre‐cracked tensile steel plates

The fatigue life of cracked steel members can be greatly extended by externally attached carbon fibre reinforced plastics (CFRP), which reduces the stress intensity factors (SIFs) at the crack tip. Access to cracks is sometimes limited and the CFRP has to be attached away from the cracks. There is a lack of knowledge on SIFs for such strengthening scheme. This paper presents the effects of CFRP bond locations on the Mode I SIF of centre‐cracked tensile (CCT) steel plate. The Mode I SIF at the crack tip is calculated using the finite element (FE) models. A correction factor is introduced as a function of CFRP bond location and crack length. The FE results are compared and agree well with experimental tests conducted by the authors. By combining with another two factors (one considering CFRP mechanical properties and the other considering CFRP bond width) derived previously by the authors, SIF formulae are proposed for CFRP reinforced CCT steel plates.     

Evaluation of C* integral for interacting cracks in plates under tension

A closed form solution for C* integral of two interacting cracks in plates under tension is developed on the basis of reference stress method. Comprehensive finite element (FE) creep analyses are carried out to provide the benchmark of the interaction evaluation of multiple cracks. Results indicate that more pronounced interaction is observed between the C* of double cracks and that of a single crack compared to that denoted by stress intensity factor (SIF). Overall good agreement is achieved between the proposed method for C* of multiple crack interaction and the FE results which provides confidence in practical application.     

On the through-thickness crack with a curve front in center-cracked tension specimens

Three-dimensional finite element analyses are performed on through-thickness cracks with slightly wavy front in center-cracked plates. Considering there is an inherent relationship between the crack shape and the corresponding stress intensity factor (SIF) distribution of a crack, the curved configuration of the crack is determined using a heuristically derived iterative procedure if the SIF distribution function is known. Several simple SIF distribution functions, for instance the constant SIF distribution along the crack front, are assumed to determine the crack shape. Under the assumption that the rate of fatigue crack growth depends on the SIF range or the effective SIF range, possible effects of plate thickness, crack length and crack closure level gradient on the behaviour of crack tunneling are investigated. The stability of the curved shape of a through-thickness crack in fatigue is also discussed, i.e. whether a crack can maintain its shape satisfying the conditions of constant SIF distribution or other distribution along the crack front during fatigue growth. This study will be useful for a better understanding of the behaviour of crack tunneling and help to evaluate the validity of the two-dimensional linear elastic fracture mechanics in cracked plates.     

Effect of biaxial loads on elastic-plastic <Emphasis Type="Italic">J</Emphasis> and crack tip constraint for cracked plates: finite element study

Based on detailed two-dimensional (2-D) and three-dimensional (3-D) finite element (FE) analyses, this paper attempts to quantify in-plane and out-of-plane constraint effects on elastic-plastic J and crack tip stresses for a plate with a through-thickness crack and semi-elliptical surface crack under positive biaxial loading. For the plate with a through-thickness crack, plate thickness and relative crack length are systematically varied, whereas for the plate with a semi-elliptical surface crack, the relative crack depth and aspect ratio of the semi-elliptical crack are systematically varied. It is found that the reference stress based approach for uniaxial loading can be applied to estimate J under biaxial loading, provided that the limit load specific to biaxial loading is used, implying that quantification of the biaxiality effect on the limit load is important. Investigation on the effect of biaxiality on the limit load suggests that for relatively thin plates with small cracks, in particular with semi-elliptical surface cracks, the effect of biaxiality on the limit load can be neglected for positive biaxial loading, and thus elastic-plastic J for a biaxially loaded plate could be estimated, assuming that such plate is subject to uniaxial load. Regarding the effect of biaxiality on crack tip stress triaxiality, it is found that such effect is more pronounced for a thicker plate. For plates with semi-elliptical surface cracks, the crack aspect ratio is found to be more important than the relative crack depth, and the effect of biaxiality on crack tip stress triaxiality is found to be more pronounced near the surface points along the crack front.     

Sustained Mode I and Mode II fatigue crack growth in flat plates

This study was conducted to contribute to the understanding of fatigue crack growth under mixed mode loading. This was accomplished by developing and analyzing a flat plate specimen capable of maintaining crack growth on a plane oblique to the direction of the applied load. Several specimens were built and exposed to controlled fatigue loading in the laboratory. These specimens were then modeled using finite elements to determine the stress intensity factors (SIF). For the “Mode I/Mode II” specimens developed, the crack was forced to grow in a direction other than perpendicular to the load. The resulting crack front did not remain straight and flat, but stabilized into a curved or warped shape. Based on finite element analyses of these curved specimen cracks, it is concluded that the SIR were predominantly Mode I, with the Mode II and III SIR being negligible.     

Near-threshold fatigue propagation of physically through-thickness short and long cracks in a low alloy steel

In this paper, the near-threshold fatigue behavior of physically through-thickness short cracks and of long cracks in a low alloy steel is investigated by experiments in ambient air. Physically through-thickness short fatigue cracks are created by gradually removing the plastic wake of long cracks in compact tension specimens. The crack closure is systematically measured using the compliance variation technique with numerical data acquisition and filtering for accurate detection of the stress intensity factor (SIF) at the crack opening. Based on the experimental results, the nominal threshold SIF range is shown to be dependent on the crack length and the characteristic of the crack wake which is strongly dependent on the loading history. The effective threshold SIF range and the relation between the crack propagation rate and the effective SIF range after the crack closure correction are shown to be independent on crack length and loading history. The shielding effect of the crack closure is shown to be related to the wake length and load history. The effective threshold SIF range and the relationship between the crack growth rate and the effective SIF range appear to be unique for this material in ambient air. These properties can be considered as specific fatigue properties of the couple material/ambient air environment.     

Stress intensity factor for center-cracked plate with crack surface interference

This paper presents a simple and physically acceptable analysis of stress intensity factor (SIF) for the center-cracked infinite and finite-width plates. The analysis includes the effect of crack surface interference (i.e., the upper and lower crack surfaces are not allowed to overlap) that influences both the SIF at the tension-side crack tip and the crack opening displacement (COD) profile. For an infinite plate, exact solutions are obtained by superimposing the classical (overlapping) solutions. For a finite-width plate, where the SIF solutions cannot be found in closed form, the solutions are carried out numerically. The overlapping SIF solutions from the weight function method are used. An example is given for the case of a finite-width plate under bending. It was found that the overlapping solutions underestimate the stress intensity factor at the tension-side crack tip up to 15%. The analysis results are also compared with the finite element solutions for verification purpose.     

Stress intensity factor for a crack emanating from the root of a semi-circular edge notch in a tension plate

The stress intensity factor (SIF) for a crack radially emanating from the root of a semi-circular edge cut-out is determined theoretically. Firstly, the advantage of using finite Mellin type transform is demonstrated mainly to resolve the boundary conditions for a plate with an edge notch. Secondly, this method has been used to calculate the SIF of a crack in a semi-infinite tension plate, which involves the solution of a Weiner-Hopf integral equation. Finally these analytical results are verified through transmission photoelasticity experiments on model specimens. It is concluded that the SIF increases less rapidly compared to Bowie's problem [1] and there is no appreciable effect of the cut-out on the stress field close to the tip region when the crack size is about seven times the edge-notch radius (Fig. 1).     

Crack face pressure loading of semielliptical cracks located along the bore of a hole

The finite element-alternating method is used to obtain Mode I stress intensity factors for semielliptical surface cracks centered along the bore of a hole in a large plate of finite thickness. The variation in stress intensity factor around the crack perimeter is provided for two loading configurations: remote uniaxial tension and a crack face pressure described by a third order polynomial. It is suggested that the tabulated solutions for the crack face pressure loading represent “general” results which can be used with linear superposition to compute stress intensity factors for many other practical loading configurations. Two example problems describe application of the superposition procedure with the crack face pressure results.