Rectangular tensile sheet with single edge crack or edge half-circular-hole crack

By using the displacement discontinuity method with crack-tip elements (a boundary element method) proposed recently by the author, this note presents the stress intensity factors (SIFs) of a rectangular tensile plate with single edge crack. Further this note studies the SIFs of crack emanating from an edge half-circular hole. By comparing the calculated SIFs of the single edge half-circular-hole crack with those of the single edge crack, a shielding effect of the half-circular hole on the SIFs of the single edge crack is discussed. It is found that the boundary element method is simple, yet accurate for calculating the SIFs of complex crack problems in finite plate.     

多孔复合材料机械连接件弹性接触内力和应力分析

本文用有限元混合法给出了多孔复合材料机械连接件弹性接触内力和应力分析。文中详细讨论了销钉弹性、层合板铺设方式和摩擦对各孔钉载分配及孔边应力的影响,并给出了多孔连接件向单孔连接件简化的条件。      

A numerical analysis of cracks emanating from a square hole in a rectangular plate under biaxial loads

This note concerns with stress intensity factors of cracks emanating from a square hole in rectangular plate under biaxial loads by means of the boundary element method which consists of the non-singular displacement discontinuity element presented by Crouch and Starfied and the 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 boundary. The present numerical results illustrate that the present approach is very effective and accurate for calculating stress intensity factors of complicated cracks in a finite plate and can reveal the effect of the biaxial load and the cracked body geometry on stress intensity factors.     

Effect of residual stress around cold worked holes on fracture under superimposed mechanical load

Cold working is one method used to enhance the fatigue life of holes in aerospace structures. The method introduces a compressive stress field in the material around the hole and this reduces the tendency for fatigue cracks to initiate and grow under superimposed cyclic mechanical load. To include the benefit of cold working in design the stress intensity factors must be evaluated for cracks growing from the hole edge. Two-dimensional (2D) finite element analyses have been carried out to quantify the residual stresses surrounding the cold worked hole. These residual stresses have been used in a finite element calculation of the effective stress intensity factor for cracks emanating from the hole edge normal to the loading direction. The results of the 2D analysis have been compared with those derived using a weight function method. The weight function results have been shown always to underestimate the stress intensity factor. A three-dimensional (3D) FEA has been carried out using the same technique for stress intensity factor evaluation to investigate the effect of through thickness variation of residual stress. Stress intensity factors calculated with the 3D analysis are generally higher than those calculated using the 2D analysis.     

A weight function method for mixed modes hole‐edge cracks

A weight function approach is proposed to calculate the stress intensity factor and crack opening displacement for cracks emanating from a circular hole in an infinite sheet subjected to mixed modes load. The weight function for a pure mode II hole‐edge crack is given in this paper. The stress intensity factors for a mixed modes hole‐edge crack are obtained by using the present mode II weight function and existing mode I Green (weight) function for a hole‐edge crack. Without complex derivation, the weight functions for a single hole‐edge crack and a centre crack in infinite sheets are used to study 2 unequal‐length hole‐edge cracks. The stress intensity factor and crack opening displacement obtained from the present weight function method are compared well with available results from literature and finite element analysis. Compared with the alternative methods, the present weight function approach is simple, accurate, efficient, and versatile in calculating the stress intensity factor and crack opening displacement.     

Numerical analysis of the influence of crack surface roughness on the crack path

The influence of the 3D frictional crack surface interaction on the fracture mechanical parameters as well as on the crack path is numerically investigated. For the solution of the boundary value problem the 3D dual boundary element method in terms of the discontinuous formulation is utilized. This method is especially suited for contact problems because it directly deals with the discontinuities at the crack surfaces. The contact problem is solved by the application of the penalty method. Coulomb’s frictional law is utilized for the consideration of the dissipative nature of friction. For discrete steps within one load cycle the stress intensity factors are determined by an extrapolation procedure from the stress field. Based on the analysis of a load cycle, the cyclic stress intensity factors are obtained. For the simulation of crack propagation an implicit time integration scheme of a crack propagation law implemented in terms of a predictor-corrector scheme is applied. The influence of the crack surface roughness on the crack path is shown by numerical examples.     

Strain energy release rate determination of prescribed cracks in adhesively-bonded single-lap composite joints with thick bondlines

An analytical model for determining the strain energy release rate due to a prescribed crack in an adhesively-bonded, single-lap composite joint with thick bondlines and subjected to axial tension is presented. An existing analytical model for determining the adhesive stresses within the joint is used as the foundation for the strain energy release rate calculation. In the stress model, the governing equations of displacements within the adherends are formulated using the first-order laminated plate theory. In order to simulate the thick bondlines, the field equations of the adhesive are formulated using the linear elastic theory to allow non-uniform stress distributions through the thickness. Based on the adhesive stress distributions, the equivalent crack tip forces are obtained and the strain energy release rate due to the crack extension is determined by using the virtual crack closure technique (VCCT). The specimen geometry of ASTM D3165 standard test is followed in the derivation. The system of second-order differential equations is solved to provide the adherend and adhesive stresses using the symbolic computational tool, Maple 7. Finite element analyses using J-integral as well as VCCT are performed to verify the developed analytical model. Finite element analyses are conducted using the commercial finite element analysis software ABAQUS™. The strain energy release rates determined using the analytical method correlate well with the results from the finite element analyses. It can be seen that the same prescribed crack has a higher strain energy release rate for the joints with thicker bondlines. This explains the reason that joints with thick bondlines tend to have a lower load carrying capacity.     

Weight function for a single edge cracked geometry with clamped ends

A single edge cracked geometry with clamped ends is well suited for fracture toughness and fatigue crack growth testing of composites and thin materials. Analysis of fiber bridging phenomenon in the composites and determination of stress intensity factors due to non-uniform stress distributions such as residual and thermal stresses generally require the use of a weight function. This paper describes the development and verification of a weight function for the single edge cracked geometry with clamped ends. Finite element analyses were conducted to determine the stress intensity factors (K) and crack opening displacements (COD) due to different types of stress distributions. The weight function was developed using the K and COD solution for a constant stress distribution. K and COD predicted using this weight function correlated well with the finite element results for non-uniform crack surface stress distributions.     

Fatigue crack propagation simulation in an aircraft engine fan blade attachment

This paper discusses the computation of three-dimensional fatigue crack growth rates in a typical military aircraft engine fan blade attachment under centrifugal and aerodynamic loads. The three-dimensional crack growth simulations utilize FRANC3D, a state-of-the-art crack propagation software developed at Cornell University, which uses boundary elements and linear elastic fracture mechanics. With an existing three-dimensional finite element contact stress analysis with a prescribed coefficient of friction (COF) along the contact surface, the displacements and stress intensity factors are calculated on the crack leading edge to yield crack propagation trajectories and growth rates. Due to complex geometry of the fan blade attachment and loading conditions, all three-fracture modes are considered and the associated stress intensity factors (SIF) are calculated using the Crack Opening Displacement (COD) approach. Crack propagation trajectories under mixed-mode conditions are obtained using the planar and maximum tangential stress crack-extension criteria. The fatigue crack in the blade attachment is subjected to an over speed mission cycle that includes high cycle frequencies (i.e., spectrum load) and the crack growth rate is predicted utilizing the Forman–Newman–de Koning (FNK) model. Scanning Electron Microscope (SEM) images of a cracked component from an engine ASMET (Accelerated Simulated Mission Endurance Test) are used to evaluate and compare the simulation results. The calculated SIF's from the simulations indicate a strong Mode-I (K_I) and Mode-III (K_III) interaction at the edge of contact (EOC). However, on the free surface it is primarily a crack opening (K_I) condition only. The crack growth rates are determined using the planar extension criterion which correlates better with the test data than the maximum tangential stress extension criteria.     

A study of crack growth retardation due to artificially induced crack surface contact

The contact of the cracked surfaces during a part of a loading cycle generally results in a reduced crack growth rate. A critical experiment was designed to evaluate the influence of the crack surface contact on crack growth. A round compact specimen made of 1070 steel with a round hole at the wake of the fatigue crack was designed. Two mating wedges were inserted into the hole of the specimen while the external load was kept at its maximum in a loading cycle. In this way, the wedges and the hole in the specimen were in firm contact during the entire loading cycle in the subsequent loading. Experiments showed that the addition of the wedges resulted in a reduction of crack growth rate in the subsequent constant amplitude loading. However, crack growth did not arrest. With the increase in the subsequent loading cycles, crack growth rate increased. The traditional crack closure concept cannot explain the experimental phenomenon because the effective stress intensity factor range was zero after the insertion of the wedges. The detailed stress–strain responses of the material near the crack tip were analyzed by using the finite element method with the implementation of a robust cyclic plasticity theory. A multiaxial fatigue criterion was used to determine the fatigue damage based upon the detailed stresses and strains. The crack growth was simulated and the predicted results were in good agreement with the experimental observations. It was confirmed that the stresses and strains near the crack tip governed cracking behavior. Crack surface contact reduced the crack tip cyclic plasticity and the result was the observed retardation in crack growth.     

Frit salvage innovation of cathode ray tube by initiating to crack along the seal edge

This article describes a thermal stress boundary element analysis which was carried out to select a mechanical device that would innovate the frit salvage of cathode ray tubes (CRT). The panel and funnel of a CRT are fritted at the seal edge. To reuse a CRT, the panel and funnel are separated by etching and then applied a thermal shock. This procedure is known as frit salvage. Current yield of the frit salvage success rate is about 61 %, which means that 39% of the components are lost and cannot be reused. The financial burden can be reduced if the salvage rate is improved through a better understanding of the thermal shock mechanism.

During thermal shock, a crack was observed to initiate at the end of the axes and traverses towards the corner along the seal end. A CRT is lost when the crack travels along the diagonal corner. A thermal stress analysis was carried out using the boundary element method. The analysis determined that maximum stresses are located near the blend radius on the panel skirt at approximately 1–2 inches on either side of the diagonal. Also it was found that the stresses along the frit seal edge are uniform from the end of the major and minor axes towards the diagonal. The stresses at the diagonal location were reduced by approximately 29%. The direction and location of the crack obtained by the boundary element anlysis were consistent with those observed in the frit salvage procedure.

The crack's promoter along the seal edge is used to protect thermal shock on the outer face panel near the heel radius in order to prevent thermal failure on the diagonal corner. The resultant stresses' directions through the crack's promoter have two holds. The directions of the dominant high stresses at the seal edge through the crack's promoter were changed perpendicular to the seal edge. Therefore, it induces the crack to follow along the seal edge. Secondly, the stresses' directions on a diagonal corner have also been changed to be perpendicular to the seal edge. Consequently, this will prevent the crack along the seal edge from propagating to the diagonal corner, as the crack will advance normally to the maximum principal stress. The crack continues to traverse along the dominant high stress lines at the seal edge as required for a successful frit salvage process. This application is an excellent example of the advantages of using the boundary element method in an industrial setting.     

Application of the weight function method to two-dimensional elastodynamic fracture mechanics

In this paper, the weight function method is used for two-dimensional mixed-mode crack analyses of clastostatic and elastodynamic problems. By the use of the Laplace transformation method and an indirect boundary element method, the dynamic stress intensity factors for a finite sheet containing a central or an edge crack are evaluated. A Green's function method is introduced which depends on the weight function for an impulsive applied load. The Green's function can be used to determine stress intensity factors for arbitrary time dependence of the boundary conditions. The stress intensity factors obtained by the weight function method are compared where possible, with existing solutions.     

Stress analysis of single-bolt, single-lap, countersunk composite joints with variable bolt-hole clearance

Single-lap, carbon-epoxy joints with countersunk fasteners were modelled using the nonlinear finite element code Abaqus. A highly-detailed analysis of the stress distribution at the countersunk hole boundary is provided. Bolt-hole clearance, which arises due to limitations in manufacturing capabilities, is modelled extensively. Clearance levels both inside and outside typical aerospace fitting tolerances are studied and the finite element model is validated with experimental data. Plots of radial stress in each ply of the countersunk laminate show the load transfer to be severely localised, with only a few plies bearing the majority of the load. The inclusion of clearance in the model was shown to result in far higher radial stresses compared to those in the neat-fit joint model. An associated loss in joint stiffness of more than 10% was recorded for the highest clearance considered (240 μm). Finally compressive through-thickness stresses are shown to be present at the damageable region of the countersunk hole, and increase with bolt-hole clearance. These compressive stresses, which are an indicator of lateral constraint, are seen to suppress “brooming” failure in the countersunk laminate.     

Weight functions and stress intensity factors for the single crack round-ended straight lug

Investigations were performed for the round-ended straight attachment lug with a single crack emanating from the hole with the weight function method. The weight functions, covering the geometries from W/D=1.5 to W/D=4.0, were generated from the results obtained with a boundary element method using the approximate weight function technique. The results have been given both in the form of analytical weight functions and tabulated dimensionless stress intensity factors for simple normalized powers of the crack line loading. This is a simple straight forward procedure to calculate stress intensity factors once the crack line loading is approximated by a polynomial. The present method is also valid for deriving stress intensity factors and weight functions for general crack configurations.     

A wide range weight function for a single edge cracked geometry with clamped ends

A single edge cracked geometry with clamped ends is well suited for fracture toughness and fatigue crack growth testing of composites and thin materials. Stress intensity factors may be determined by the weight function method. A weight function for the single edge cracked geometry with clamped ends is developed and verified in this paper. It is based on analytical forms for the reference stress intensity factor and crack mouth opening displacement. The analytical forms are shown to be valid, by comparison with finite element results, over a wide range of crack depths and plate aspect ratios. Use of the analytical form enables the weight function to be calculated for any plate aspect ratio without the need for preliminary finite element analysis. Stress intensity factors and crack mouth opening displacements, predicted using this weight function, correlated well with finite element results for non-uniform crack surface stress distributions. This revised version was published online in August 2006 with corrections to the Cover Date.     

Stress intensity factors for cracks in thin plates

This paper presents stress intensity factor solutions for several crack configurations in plates. The loadings considered include internal pressure, and also combined bending and tension. The dual boundary element method is used to model the plate and mixed mode stress intensity factors are evaluated by a crack surface displacement extrapolation technique and the J-integral technique. Several cases including centre crack, edge crack and cracks emanating from a hole in finite width plates are presented.     

Influence of foreign object damage on fatigue crack growth of gas turbine aerofoils under complex loading conditions

Foreign object damage (FOD) has been identified as one of the main life limiting factors for aeroengine blades, with the leading edge of aerofoils particularly susceptible. In this work, a generic edge ‘aerofoil’ geometry was utilized in a study of early fatigue crack growth behaviour due to FOD under low cycle fatigue (LCF), high cycle fatigue (HCF) and combined LCF and HCF loading conditions. Residual stresses due to FOD were analyzed using the finite element method. The longitudinal residual stress component along the crack path was introduced as a nodal temperature distribution, and used in the correction of the stress intensity factor range. The crack growth was monitored using a nanodirect current potential drop (DCPD) system and crack growth rates were correlated with the corrected stress intensity factor considering the residual stresses. The results were discussed with regard to the role of residual stresses in the characterization of fatigue crack growth. Small crack growth behaviour in FODed specimens was revealed only after the residual stresses were taken into account in the calculation of the stress intensity factor, a feature common to LCF, HCF and combined LCF + HCF loading conditions.     

BOUNDARY ELEMENT ANALYSIS OF CRACK INTERSECTION, INITIATION AND GROWTH AT PIN LOADED HOLES

Abstract The dual boundary element method for the analysis of cracks in linear elastic materials has been previously generalised by the authors to allow for automatic remeshing when crack tips intersect other cracks or boundaries, and initiation and growth of small cracks at positions of high stress concentration. The new cracks are assumed to result from sudden events such as an overload or the subsequent stress redistribution when cracks intersect other cracks or holes. In this paper a crack at the edge of one hole in a row of pin-loaded holes is investigated; various values are considered for the stress at which new cracks may initiate. Two rows of aligned or staggered holes are examined also. The spacing between the holes was typical for lines of holes in overlap joints in plates. For the same load transfer between the plates, new cracks are initiated less readily and grow more slowly for a double row of aligned holes than for a single row or for staggered rows.     

螺栓夹持填充孔复合材料层压板层间应力分析

为了研究含螺栓复合材料层压板连接接头的层间应力分布规律,提出了采用填充孔形式和虚拟界面层方法求解层间应力,并建立了三维有限元模型对受面内压缩载荷的螺栓夹持填充孔层压板进行分析。结果表明:层压板层间应力集中不仅发生在孔边,还会在螺栓头边缘附近出现,且夹持力越大螺栓头边缘附近的层间应力集中越严重;合理的螺栓夹持力能改善孔边应力状态,提高孔边抵抗分层的能力,但无法改善螺栓头边缘附近的层间剪切应力集中状态。因此,在进行含螺栓夹持的层压板机械连接结构孔边层间强度设计时,也要考虑螺栓头边缘的层间剪切应力集中问题,以提高复合材料结构的安全性。     

Dual boundary element method and finite element method for mixed‐mode crack propagation simulations in a cracked hollow shaft

Three‐dimensional mixed‐mode crack propagation simulations were performed by means of the dual boundary element method code BEASY and 2 finite element method‐based crack propagation codes: ZENCRACK (ZC) and CRACKTRACER3D (CT3D). The stress intensity factors (SIFs) along the front of an initial semielliptical crack, initiated from the external surface of a shaft, were calculated for 4 different load cases: bending, press fit, shear, and torsion. The methods used for the SIF assessment along the crack front were the J‐integral for BEASY and ZC and the quarter point element stress method for CT3D. Subsequently, crack propagation simulations were performed, with the crack growth rate evaluated by using Paris' law, calibrated for the material at stake (American Society for Testing and Materials A469 steel). The kink angles were evaluated by using the minimum strain energy density and maximum tangential stress criteria for BEASY, the maximum energy release rate and maximum tangential stress for ZC, and the maximum principal asymptotic stress for CT3D. The results obtained in terms of SIFs and crack propagation life show very good agreement among the 3 codes. Also, the shape of the propagated crack, which is significantly out‐of‐plane for the shear and torsion loading, matched very well.