Microsoft Word - numero_36_art_16 R. Citarella, Frattura ed Integrità Strutturale, 36 (2016) 160-167; DOI: 10.3221/IGF-ESIS.36.16 160 Multiple crack propagation by DBEM in a riveted butt-joint: a simplified bidimensional approach R. Citarella (https://orcid.org/0000-0003-3167-019X) Dept. of Industrial Engineering, University of Salerno, Via Giovanni Paolo II, 132, Fisciano (SA), Italy rcitarella@unisa.it ABSTRACT. A Multi-Site Damage (MSD) crack growth simulation is presented, carried out by means of Dual Boundary Element Method (DBEM), in a two-dimensional analysis of a cracked butt-joint made of aluminium 2024 T3. An equivalent crack length is proposed for an approximated 2D analysis of a 3D problem where the crack front assumes a part elliptical shape due to secondary bending effects. The assumptions made to perform such simplified bidimensional analyses are validated by comparing numerical results with experimental data, the latter obtained from a fatigue tested riveted butt-joint. KEY WORDS: DBEM; MSD; LEFM. INTRODUCTION omplex engineering structures will often present microscopic flaws that can be caused by the manufacturing process or by fatigue, impact or environmental effects such as corrosion. It is therefore crucial to be able to predict how these cracks will affect the integrity of the structure, for example to determine whether a particular crack will grow at all under service loading and, in case, what the life of the component before failure. It is evident from this considerations that the field of Fracture Mechanics will be an integral part of the design process, particularly in the aircraft industry due to the high strength but low crack resistance materials used in weight critical applications. Experimental studies have shown that, if the hypothesis of Linear Elastic Fracture Mechanics (LEFM) holds true, failure occurs when the Stress Intensity Factor (SIF) reaches a critical value. Hence the predictions of crack behaviour and the integrity of a structure are dependent on the accurate calculation of SIFs. (a) (b) Figure 1: Deformed lap-joint undergoing a traction load with highlight of sites undergoing maximum bending stresses (a) and stress distribution through the thickness of the sheet (b). C R. Citarella, Frattura ed Integrità Strutturale, 36 (2016) 160-167; DOI: 10.3221/IGF-ESIS.36.16 161 In this work a Multi Site Damage (MSD) crack growth simulation is presented, carried out by means of Dual Boundary Element Method (DBEM) [1-4] as implemented in the commercial code BEASY [5], in a two-dimensional analysis of a butt-joint, made of aluminium alloy 2024 T3, undergoing a traction fatigue load. The bidimensional model represents an approximation of the real phenomena because of the secondary bending effects, illustrated in Fig. 1 with reference to a lap-joint (but the phenomenon is completely equivalent for a butt-joint) and responsible for a non-straight crack propagating front. Such phenomena would prevent a bidimensional modelling approach, but this drawback can be circumvented, in a preliminary analysis, by adopting an equivalent straight initial crack front as explained in the following. J INTEGRAL he J integral was developed by Rice and Cherepanov to characterize fracture for two dimensional configuration. The material is elastic or non-linear elastic. The J integral is defined as: 1 1 ( )ii u J Wn t d x      where W is the strain energy density, ti are components of the traction vector and ui are components of the displacement vector. The vector n is the unit normal to the contour  . The strain energy is defined as 0 ( ) ij ij ij ijW W d       i, j = 1,2 and ij is the infinitesimal strain tensor component. The J integral is integrated along a contour  surrounding the crack tip (the ends of the  contour are on opposite faces of the crack and the contour encloses the crack tip). The crack is parallel to the x1 axis. On addition or subtraction of quantities above and below the crack axis, symmetric and anti-symmetric stress and strain derivatives are found from which symmetric stress and strain products are calculated respectively. These products can be decomposed into two parts: one consisting of symmetric stress and strain derivatives (symmetric integrands), the other consisting of anti-symmetric stress and strain derivatives (anti-symmetric integrands). The symmetric integrands can be integrated over the contour area to obtain the mode I J-integral; the anti-symmetric integrands, when integrated, will produce the mode II J-integral [6]. FATIGUE CRACK GROWTH he well-known Paris law is adopted for the crack growth rate assessment: nKC dN da  (1) with the calibration parameter C and n for Al 2024 T3 given by n=2.144 and C=4.72*10-10 [7], 22 2 IIIeff KKKK  (Tanaka formula [1]) expressed in MPamm and da in mm. It has been checked that, during the propagation, the SIF’s range is included in the interval of validity of Paris law, skipping the very initial crack growth phase because of the non-modelled threshold phenomenon and the final crack growth phase, which is dominated by very high growth rates, with Kmax approaching Kc (Kc is the material fracture toughness). PROBLEM DESCRIPTION AND PRELIMINARY RESULTS ON UNDAMAGED PANEL n the following a Multi-Site damage (MSD) crack growth simulation for a butt-joint is presented, carried out by DBEM in a two-dimensional analysis. Numerical results are then compared with corresponding experimental outcomes, obtained from a fatigue tested riveted butt-joint specimen of aluminium 2024 T3 [7]. T T I R. Citarella, Frattura ed Integrità Strutturale, 36 (2016) 160-167; DOI: 10.3221/IGF-ESIS.36.16 162 The modelled specimen is representative of the riveted joint existing between the upper shell and the lower shell on the front fuselage section, manufactured by AerMacchi for the Dornier Do328 aircraft (Fig. 2a). The corresponding DBEM model is illustrated in Fig. 2b, where internal points are placed in correspondence of the strain gauge locations (Fig. 2a) in order to compare the numerical and experimental strains and assess the model accuracy. The distance from the panel edge to the first hole was chosen from the experimenter equal to 10 mm, giving a total panel width W=280 mm (Fig. 2a). The four holes close to the plate border were cold worked in order to avoid crack initiation during the fatigue test (in Fig. 2a, holes with crack tips N. 1, 28, 29, 56) or at least strongly retard crack propagation [8]. The remote applied load is t=100 MPa (Fig. 2b) and the material properties are: Young modulus E=72500 MPa and Poisson ratio .33  . 1 28 29 56 (a) Line of internal points in correspondence of strain gauge positions Hole with no modelling of inside pin Hole with explicit modelling of inside pin (b) Figure 2: Butt joint sketch (a) with highlight of strain gauge positions and crack tip numbering (N. 1:28 in the lower hole row and N. 29:56 in the upper cracked hole row); DBEM uncracked model (b) with highlight of internal spring (yellow symbols), internal points and hole constraints (green symbols). R. Citarella, Frattura ed Integrità Strutturale, 36 (2016) 160-167; DOI: 10.3221/IGF-ESIS.36.16 163 The butt-joint is designed and loaded in the experimental tests, in such a way to reproduce, as close as possible, the real in-service stress state, like, for example, a uniform distribution of longitudinal stresses yyalong a transversal section of the undamaged joint sufficiently far from the riveted area. Along such section a row of strain gauges was placed and correspondingly a row of internal points was introduced in the DBEM model (N. 6:17 in Fig. 2b) in order to monitor the level of aforementioned stresses. The same was done to cross check the longitudinal stress distribution in a longitudinal section (internal points N. 1-5 in Fig. 2b). In particular, for the undamaged panel, the DBEM model provides the longitudinal stress distribution at strain gauge locations illustrated in Figs. 3a-b: the numerical stresses are in good agreement with the experimental stress state [7]. (a) (b) Figure 3: Membrane longitudinal stresses yy on transversal section at y=54.25 mm (a); longitudinal stresses yy on longitudinal section at x=110 mm (b). It is evident that, whenever experimental strains were available from both sides of the panel they should had been averaged in order to be comparable with numerical results (strain values are slightly different on the two panel sides because of the secondary bending). This check turned out very useful in order to calibrate the model and in particular to decide how many rivets to be explicitly modelled by inserting a pin in the corresponding hole: such pin would be constrained against y translation and disconnected from the hole upper surface by means of internal spring of negligible stiffness (Fig. 2b). Alternatively, the holes not directly involved by crack propagation can be modelled by just introducing longitudinal constraints on 180 degrees of the hole boundary and skipping the explicit pin modelling (Fig. 2b). Gap elements have also been introduced, to better tackle contact conditions but the solution improvement has been judged quite negligible (less than 2%), except in case of very short cracks initiated from the holes, more sensitive to pin- hole contact conditions (but not present in this work). For this reason, and due to the computational effort of a non-linear analysis, they have not been used anymore. R. Citarella, Frattura ed Integrità Strutturale, 36 (2016) 160-167; DOI: 10.3221/IGF-ESIS.36.16 164 The J-integral technique is adopted for SIF’s evaluation, being more stable than Crack Opening Displacement method, against a variable refinement of crack mesh. On the J-integral path, 33 integration points are used (the increment of accuracy with 66 points turn out to be negligible). The mesh used for the butt-joint is based on about 327 quadratic elements: a p-convergence study has been realized showing that cubic elements provide an accuracy improvement of less than 2% and that 2 quadratic elements per 90 degrees are sufficient on the cracked hole, except for very short cracks, where 3 elements are recommendable (possibly with a scaling ratio). All the undamaged holes are modelled with 6 elements, constrained in y-direction. NUMERICAL CRACK PROPAGATION he crack growth has been simulated considering the interval ranging between 73500 and 113885 fatigue cycles. The initial crack length in the model, at left or right hole sides, is not taken equal to that (size b in Fig. 4a) visually observed on the free surfaces at 73500 cycles [7] because of the following reasons. The crack length monitored in the experimental phase is related to the external surface (“penetrated crack”), but the crack first appear on the faying surface because here undergoes the primary tensile stress and, in addition, the secondary bending stress (Fig. 1). As a matter of fact, due to the superimposed secondary bending the crack assumes a part elliptical shape (Fig. 4), in such a way that it is already abundantly extended on the faying surface when appearing on the external surface.   b (a) (b) Figure 4: Geometric sizes of rivet hole (a) and crack shape development (b) [9]. It is possible to make a forecast of the hidden part crack length, right before the appearance of the external crack front (Fig. 4b), observing that the ratios between the two ellipse semi-axis, crack depth a and crack length c and between a and t (specimen thickness), even if variable at initiation, assume always the same values for a given bending factor k=max tension/max bending when the crack is on the verge to become through the thickness (e.g. for k=1 we obtain a/c=0.575 and a/t=0.88) [9]. In our particular case, for a plate thickness t=1.2 mm, when each crack appears on the visible plate side with a length b (as reported by the experimenter), by using the aforementioned relationships, it is possible to know the size of the part elliptical crack on the faying side, obtaining c=a/0.575=0.88*t/0.575=1.84 mm (Figs. 4a-b). Keeping into account these data and the shape of the real hole (Fig. 4a) against the simplified modelled geometry (a cylinder of radius 2.4 mm), when modelling an equivalent straight initial crack we have prolonged of 1 mm the initial crack length value b measured by the experimenter. During the propagation some cracks will link-up and the criteria chosen to assess such condition is the overlapping of the crack tip plastic zones: L=rp1+rp2 where L is the residual ligament and rp=(Keq2/y2)/; y is the yield stress. When two cracks link-up a new crack initiation is postulated from the hole side opposite to the crack. T R. Citarella, Frattura ed Integrità Strutturale, 36 (2016) 160-167; DOI: 10.3221/IGF-ESIS.36.16 165 Intermediate deformed plots, during the multi side crack propagation, are presented in Figs. 5a-g, where each new crack introduced in the model comes from the experimental recordings or is postulated to appear when two cracks link-up (as previously said). a b c d e f g Constraints in the pin to model the action of the overlapped plate (not explicitly modelled) Continuity condition between hole and rivet Internal springs of negligible stiffness to model pin-hole disconnection Constraints in the hole to model the action of the overlapped plate (no pin modelling) 41 42 41 42 43 41 42 43 41 44 47 39 40 41 44 45 47 39 44 45 47 38 39 44 45 47 51 Figure 5: Deformed shape (scale factor 20) of multiple crack scenario after a number of cycles equal to: 73500 (a), with initial cracks (tip N. 41 and 42) as recorded by the experimenter; 84000 (b) with a new recorded initiation (tip N. 43); 99500 (c) with propagation of the previously indicated cracks and a new recorded initiation (tip N. 47); 105770 (d) with a link-up between tips N. 42 and 43 and corresponding new postulated initiation (tip N. 44); 108000 (e) with new recorded initiations (tip N. 39, 40 and 45); 110500 (f) with a link-up between tips N. 40 and 41 and a new recorded initiations (Tip N. 51); 112000 (g) with new recorded initiations (Tip N. 38). R. Citarella, Frattura ed Integrità Strutturale, 36 (2016) 160-167; DOI: 10.3221/IGF-ESIS.36.16 166 In Tab. 1 and Tab. 2 the crack length vs. cycles and the related equivalent SIFs are respectively shown. With the aforementioned bidimensional approach a good agreement between numerical and experimental [7] crack scenarios is obtained and it is possible to capture the essential features of the response even if the secondary bending is not modelled. TIP N 38 39 40 41 42 43 44 45 47 51 73500 1.82 1.42 78810 2.63 2.28 84030 3.36 3.04 1.16 88830 4.46 4.18 2.24 93280 5.20 4.94 2.84 97330 5.92 5.70 3.43 99500 6.54 6.37 3.94 1.19 102850 7.23 7.14 4.52 1.77 105770 7.86 7.91 5.13 1.50 2.17 108000 1.15 1.12 9.27 link-up link-up 2.28 1.19 3.20 110500 1.83 1.94 10.30 3.25 1.78 3.60 1.12 112000 1.20 2.59 link-up link-up 4.02 2.14 4.00 1.43 113885 1.81 3.57 5.02 2.61 4.36 1.69 Table 1: Crack propagation lengths (mm) versus load cycles. TIP N 38 39 40 41 42 43 44 45 47 51 73500 389 400 78810 397 403 84030 406 410 391 88830 428 435 390 93280 443 453 401 97330 456 472 415 99500 471 496 438 400 102850 474 501 476 393 105770 505 600 549 419 390 108000 502 555 611 417 468 418 110500 535 631 709 652 484 430 396 112000 542 720 729 527 500 400 113885 548 740 750 553 572 400 Table 2: Equivalent SIF’s (MPamm) versus number of cycles. R. Citarella, Frattura ed Integrità Strutturale, 36 (2016) 160-167; DOI: 10.3221/IGF-ESIS.36.16 167 CONCLUSIONS t is important to point out the need, in a bidimensional crack propagation model for a butt-joint, to model the part elliptical crack front at initiation with an “equivalent” prolonged straight crack. As a matter of fact, in such case the penetrated front experiences a lower stress, compared with the hidden front surface and consequently there is not the “catch up” behaviour that is typical of the pure tensile case, where the “penetrated” crack front reaches immediately the more prominent front surface because of the higher SIFs [9]. Even with the illustrated simplified bidimensional approach it is possible, for such kind of problems, to obtain a satisfactory agreement with experimental crack growth rates. Moreover, very short run times are needed to run the whole propagation and an easy preprocessing phase is enabled by the DBEM approach: an automatic remeshing is possible as the crack grows and the manual intervention is just necessary to initiate new cracks. A more accurate DBEM bidimensional approach (with respect to the simplified approach adopted in this work) to model the butt joint assembly is also possible when needed: each layer can be considered as an individual two-dimensional structure; individual layers can be explicitly modelled and connected with rivets; by gap elements can be used at the interface pin-hole…[10]. The very basic approach presented here, with related very short run times, becomes mandatory in case of a probabilistic approach to crack propagation simulation where hundred thousands of such simulations are to be performed (e.g. when resorting to Monte Carlo method…) [11]. REFERENCES [1] Citarella, R., Perrella, M., Multiple surface crack propagation: numerical simulations and experimental tests, Fatigue and Fracture of Engineering Material and Structures, 28 (2005) 135-148. DOI: 10.1111/j.1460-2695.2004.00842.x. [2] Citarella, R., Cricrì, G., Armentani, E., Multiple crack propagation with Dual Boundary Element Method in stiffened and reinforced full scale aeronautic panels, Key Engineering Materials, 560 (2013) 129-155. DOI: 10.4028/www.scientific.net/KEM.560.129. [3] Citarella, R., MSD Crack propagation on a repaired aeronautic panel by DBEM, Advances in Engineering Software, 42 (10) (2011) 887-901. DOI: 10.1016/j.advengsoft.2011.02.014. [4] Citarella, R., Non Linear MSD crack growth by DBEM for a riveted aeronautic reinforcement, Advances in Engineering Software, 40(4) (2009) 253–259. DOI: 10.1016/j.advengsoft.2008.04.007. [5] BEASY V10r14, Documentation, C.M. BEASY Ltd, (2011). [6] Rigby, R.H., Aliabadi, M.H., Decomposition of the mixed-mode J-integral-revisited, Int. J. Solids Structures, 35(17) (1998) 2073-2099. DOI: 10.1016/S0020-7683(97)00171-6. [7] Cattaneo, G., Cavallini, G., Galatolo, R., SMAAC (Testing of “Simple” Specimens), Document No. SMAAC-TR-3.2- 07-1.3/AEM, (1998). [8] Citarella, R., Cricrì, G., Lepore, M., Perrella, M., Assessment of Crack Growth from a Cold Worked Hole by Coupled FEM-DBEM Approach, Key Engineering Materials, 577-578 (2014) 669-672. DOI: 10.4028/www.scientific.net/KEM.577-578.669. [9] Fawaz, S.A., Fatigue Crack Growth in Riveted Joints, Doctoral Thesis, Delft University Press, The Netherlands, (1997). [10] Armentani, E., Citarella, R., DBEM and FEM analysis on non-linear multiple crack propagation in an aeronautic doubler-skin assembly, International Journal of Fatigue, 28 (2006) 598–608. DOI: 10.1016/j.ijfatigue.2005.06.050. [11] Citarella, R., Apicella, A., Advanced design concepts and maintenance by integrated risk evaluation for aerostructures, Structural Durability & Health Monitoring, 2(3) (2006) 183-196. DOI: 10.3970/sdhm.2006.002.183. I << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice