Microsoft Word - numero_38_art_25 D. Carrella-Payan et alii, Frattura ed Integrità Strutturale, 38 (2016) 184-190; DOI: 10.3221/IGF-ESIS.38.25 184 Focussed on Multiaxial Fatigue and Fracture Implementation of fatigue model for unidirectional laminate based on finite element analysis: theory and practice D. Carrella-Payan, B. Magneville, M. Hack, C. Lequesne Siemens PLM Software Delphine.carrella@siemens.com, Benoit.magneville@siemens.com, Michael.hack@siemens.com, cedric.lequesne@siemens.com T. Naito, Y. Urushiyama Honda R&D Co Ltd, Tochigi, Japan Tadashi_Naito@n.t.rd.honda.co.jp, Yuta_Urushiyama@n.t.rd.honda.co.jp W. Yamazaki, T. Yokozeki University of Tokyo, Japan (Jp) yamazaki_w@aastr.t.u-tokyo.ac.jp, yokozeki@aastr.t.u-tokyo.ac.jp W. Van Paepegem Ghent Univeristy, Belgium Wim.VanPaepegem@UGent.be ABSTRACT. The aim of this study is to deal with the simulation of intra- laminar fatigue damage in unidirectional composite under multi-axial and variable amplitude loadings. The variable amplitude and multi-axial loading is accounted for by using the damage hysteresis operator based on Brokate method [6]. The proposed damage model for fatigue is based on stiffness degradation laws from Van Paepegem combined with the ‘damage’ cycle jump approach extended to deal with unidirectional carbon fibres. The parameter identification method is here presented and parameter sensitivities are discussed. The initial static damage of the material is accounted for by using the Ladevèze damage model and the permanent shear strain accumulation based on Van Paepegem’s formulation. This approach is implemented into commercial software (Siemens PLM). The validation case is run on a bending test coupon (with arbitrary stacking sequence and load level) in order to minimise the risk of inter-laminar damages. This intra-laminar fatigue damage model combined efficient methods with a low number of tests to identify the parameters of the stiffness degradation law, this overall procedure for fatigue life prediction is demonstrated to be cost efficient at industrial level. This work concludes on the next challenges to be addressed (validation tests, multiple-loadings validation, failure criteria, inter-laminar damages…). KEYWORDS. Composite; Fatigue; Variable Amplitude; Stiffness degradation. Citation: Carrella-Payan, D., Magneville, B., Hack, M., Lequesne, C., Naito, T., Urushiyama, Y., Yamazaki, W., Yokozeki, T., Van Paepegem, W., Implementation of fatigue model for unidirectional laminate based on finite element analysis: theory and practice, Frattura ed Integrità Strutturale, 38 (2016) 184-190. Received: 10.05.2016 Accepted: 25.06.2016 Published: 01.10.2016 Copyright: © 2016 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. D. Carrella-Payan et alii, Frattura ed Integrità Strutturale, 38 (2016) 184-190; DOI: 10.3221/IGF-ESIS.38.25 185 INTRODUCTION he increase of lightweight material in transportation industries is today facing more than ever questions on fatigue life prediction of composite structures. The critical step towards accurate prediction is to reproduce the loading conditions undergone by the composite component. In automotive application, the challenge is related to the variability of those conditions: multi-axial and variable amplitude on long duration fatigue loading. This is why Siemens PLM software has developed an innovative composite fatigue CAE methodology (patent pending) keeping track of the material degradation under such conditions. Sevenois [1] reviewed and compared the state of the art for fatigue model techniques of woven and UD composite. The study concluded that out of the four modelling methodologies (fatigue life, residual strength, residual stiffness and mechanistic model) the residual stiffness models are suitable for mechanical performance using experimental data and can also be combined with residual strength approach. The presented methodology is based on residual stiffness fatigue law combined with an efficient damage operator approach to calculate the residual stiffness. This approach will be able to perform fatigue simulations for variable amplitude loads and will allow ply-stacking optimization without additional testing or material characterizations. INTRA-LAMINAR FATIGUE SOLUTION WITH DAMAGE JUMP he intralaminar fatigue model strategy is herein presented in the following order: definition of the stiffness degradation law, then calculation optimization algorithm (N-Jump and damage jump combined to the damage operator) and finally, description of parameter identification procedure used in Siemens PLM commercial software. Fatigue and Stiffness Degradation – Theory Fatigue damage laws The damage evolution law is based on the work of Van Paepegem for woven glass fibers [2]. Three intra-laminar damage variables D11, D22 and D12 are defined at ply level and linked to the stress tensor by the following behavior law, Eq.(1)  pHCH     (1) where C is the Stiffness tensor, εp is a permanent strain tensor and H is defined as D D D 11 22 12 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 H 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1                     (2) In [2] the damage variables Dij were split into a positive and a negative part ( ijd  and ijd  ), where the positive part increases when the stress is positive, and the negative one when it is negative. At the end the two parts were added, including a crack closure coefficient for combination of tension of compression. In this work, only positive stress ratios are used, which means that there is no switch between tension and compression over a cycle and simplifies the problem. Either the stress is always positive and the damage Dij is equal to ijd  ; or it is always negative and it is equal to ijd  , Eq.(3) . Besides, the formulations of Van Paepegem [2] must be adapted to unidirectional plies. First of all, to account for the high in-plane orthotropy of unidirectional plies, independent ci parameters are defined for the three components of the damage. Therefore, fifteen parameters are used (ci,jk) instead of five. For the same reason, the coupling between D11 and D12 which T T D. Carrella-Payan et alii, Frattura ed Integrità Strutturale, 38 (2016) 184-190; DOI: 10.3221/IGF-ESIS.38.25 186 was implemented for woven is removed for UD: in woven fabrics, matrix de-cohesion clearly affects the stiffness in longitudinal and transverse directions, whereas in unidirectional plies, the effect of matrix degradation on longitudinal behavior can be neglected. However, the coupling between D22 and D12 remains mandatory and has been maintained. Finally, the deletion of this coupling imposes the addition of a propagation term in the formulation of D12, so that a pure shear load in a ply remains able to lead to its collapse. With these assumptions applied to the formulations taken from [2], the evolution laws for the damage variables become:    f d d d c exp c c d exp c c dN cd d d c exp c c d exp c dN d d c D dN 211 11 1,11 11 2,11 3,11  11 11 5,11 11 4 ,11 11 3 5,11211 11 1,11 11 2,11 3,11  11 11 11 4 ,11 11 222 1,22 12 ( )   ( )   3 ( ) 1                                                        f d exp c c d exp c c D 222 22 2,22 3,22  22 22 5,22 22 4 ,22 2 2212   1                            f f cd d d c D exp c c d exp c dN D d d d c d exp c c d exp c dN d 3 5,222 222 22 1,22 22 2,22 3,22  22 22 22 4 ,2212 2 2212 2 212 12 1,12 12 12 2,12 3,12  12 12 5,12 2 12 12 ( ) 1   31 ( ) 1   2 1                                                              c d d d c d exp c c d exp c c dN d 12 4 ,12 2 212 12 1,12 12 12 2,12 3,12  12 12 5,12 12 4 ,12 2 12 12 ( ) 1   2 1                             (3) where ci,jk are the 15 fatigue material coefficients that must be identified, the fatigue failure indices Σij are the ratio between the effective stress and the ultimate strength of the material in the ij component, and ij ij ij cycle ij ij cycle ij ij ij ij ij ij ij max max where ; 0,   0     ,   0                                   (4) ij ij ij ij where 0,   0     ,   0        Accumulated permanent strain In addition to these damage evolution laws, the model takes into account the permanent strain which appears in the ply due to a cyclic shear loading. Some matrix debris formed by the shear stress is accumulated in the opening matrix cracks during tension stress [2], which leads to a non-reversible deformation of the ply. The c9 parameter drives the fatigue permanent strain accumulation following the formulation: . p t t d dd dd c max c max dN dN dN 12 12 12 9 12 9 12            . (5) D. Carrella-Payan et alii, Frattura ed Integrità Strutturale, 38 (2016) 184-190; DOI: 10.3221/IGF-ESIS.38.25 187 Initial degradation of the ply It is also important to consider the effect of the first static loading of the ply on its fatigue damage and strain behavior, as it is not included in the fatigue degradation. Therefore, the static damage and permanent strain are evaluated with a static damage model [2] by running a first non-linear analysis up to the peak load of the cyclic fatigue analysis. The resulting initial damage and permanent strain are imposed as initial state of the fatigue analysis and the first cycle can be computed with a correct stress distribution. The fatigue damage tensor Df is then added to the initial damage tensor Ds D = Ds + Df (6) The same operation is done with the permanent strain. Calculation optimization: from N-Jump to Damage Jump Block loading: N-Jump The purpose of the N-Jump algorithm is to avoid running a full FE analysis at each load cycle and to deliberately choose a few relevant load cycles only. The cycles with no significant damage growth are “jumped”. From a first FE analysis, at each Gauss point of the FE model, the theoretical number of cycles to jump NJUMP1 is estimated by extrapolating the damage, Eq. (1) and applying to Eq. (7) N NJUMP N if D D D D D if D if D 20 1 10     0     0.5    0 0.2 0.1    0.2              (7) A global cycle jump NJUMP is defined such that P% of Gauss points verify NJUMP1 < NJUMP for the three components. The value of P has been set to 5% in this study. The damage is finally extrapolated after NJUMP, using again the progressive damage formulations Eq.(3). To validate this algorithm and the value of P, three similar fatigue analyses (three points bending) have been run on the first 100 loading cycles with a 45 degrees layup; one without NJUMP, one with P=1% and one with P=5%. Fig. 1 compares the stiffness degradation observed in the three analyses and validates the accuracy of the NJUMP algorithm. Figure 1: Effect of the N-Jump algorithm on stiffness degradation. VARIABLE AMPLITUDE: DAMAGE ACCUMULATION JUMP AND DAMAGE HYSTERESIS OPERATOR n automotive industry, the synthesis of realistic fatigue loading involves complex load schedules for different roads with variable loading (Fig. 2). The aforementioned jump algorithm is given for block loading. Also, the N-Jump is calculating the damage status by extrapolation (Eq. 7). Here, the damage accumulation jump aims to accurately calculate locally the progressive damage and stiffness degradation. I D. Carrella-Payan et alii, Frattura ed Integrità Strutturale, 38 (2016) 184-190; DOI: 10.3221/IGF-ESIS.38.25 188 Figure 2: Variable Amplitude example in automotive application. For variable amplitude, traditional fatigue approaches for metallic material use SN curves, linear Miner-Palmgren damage accumulation and cycle (rainflow) based damage evaluations [3-5]. SN curves are test based on curve fitting techniques which does not take into account the loading history of the material. In 1945, Miner developed a linear damage accumulation method, based on the work of Palmgren and added the contribution of various stress amplitude loading to the damage. However, as for SN curves, the loading history of the material is not accounted for. In Rainflow counting methods the damage level depends on full closing hysteresis loops of load cycles (Fig. 3). Figure 3: Illustration of rainflow method based on stresses/strains with nested cycles. In the case of composite materials, the fatigue behavior is changing over time due to changes in the matrix damage state. When applying variable amplitude loading, the largest load cycles – that contribute to the larger amount of damage – commonly take a very long time to complete, due to the many nested cycles (Fig. 3). In this case the approach to only consider cycles when they are completed can no longer be justified. The damage hysteresis operator approach based on Brokate works [6] is able to calculate damage at defined time increment instead of ‘closed load increment’. This operator allows to up-date the damage status depending on the pre- damage and any other external factors (i.e. temperature, humidity…). These operators are therefore suited to follow the progressive damage curves and also to include damage history of the material. This approach gives good prediction when applied to temperature dependent fatigue analysis with non-linear damage accumulation [7-8-9] and for full car structures with full load histories [10]. PARAMETER IDENTIFICATION PROCEDURE Tests protocol n this study, two different tests protocols are proposed. The first one is a traditional approach with tensile tests on five layups and five load levels per layup capturing a representative span of fatigue life (based on test method [11]). This results in twenty five configurations. The second one is an application of a more innovative approach with one-I D. Carrella-Payan et alii, Frattura ed Integrità Strutturale, 38 (2016) 184-190; DOI: 10.3221/IGF-ESIS.38.25 189 sided bending tests on the same five layups, but with only one load level for each layup, so only five configurations. The idea is to assess if the load distribution on this kind of specimen can provide enough information to feed the damage model, as a three-points-bending test results in progressive load levels along the same specimen, both in tension and compression. Peak data such as load, displacement, strains from extensometers and gauges, temperature are measured from these tests at several representative cycles. The stiffness degradation of each specimen can be estimated from these evolutions. Additionally, running-in and unloading raw data are extracted to analyse the initial stiffness drop of the specimens and their final permanent strain. Parameters identification protocol The parameters identification consists in an FE-based optimization of the fifteen fatigue parameters to fit simulated stiffness degradation with experimental results. A step-by-step methodology identifying the parameters one by one from specific experimental data has been setup. As illustrated by two examples in Fig. 4, each ci parameter has a specific contribution on the numerical stiffness degradation, and therefore can be adjusted to perfectly fit with experimental results. Figure 4: Effect of c2,12 (left) and c4,12 (right) on numerical stiffness degradation Therefore, two volume finite elements models reproducing the two testing procedures have been created, and the nonlinear fatigue solver with N-Jump algorithm is used to efficiently correlate the stiffness degradation of each testing configuration by adjusting the fifteen ci,jk parameters (note that the N-Jump is herein used as tests are carried out under constant amplitude loading). An illustration of the resulting correlation between experimental and simulated stiffness evolutions is given in Fig. 5. Figure 5: Correlation between experimental and simulated stiffness evolutions ([0]20 layup in 3Pts Bending). Finally, the parameter c9 is estimated by crosschecking the raw data of the unloading of some of the specimens. D. Carrella-Payan et alii, Frattura ed Integrità Strutturale, 38 (2016) 184-190; DOI: 10.3221/IGF-ESIS.38.25 190 VALIDATION/APPLICATION he first validation has been conducted on a coupon at constant amplitude loading. The same three points bending analysis as above is run on 80 000 cycles at an imposed load level. A more complex quasi-isotropic layup [0/60/- 60]4s is now studied. The resulting stiffness degradation is compared to experimental data (Fig. 6). The good predictability illustrates the main interest of a ply-level damage law: identification is performed on specific layups, and the resulting material data remains available for any layup without additional identification. Figure 6: Predictability of stiffness reduction of a 3Pts-bending quasi-isotropic coupon. Further validation cases have been investigated (i.e. flat and V-shaped components) but the overall stiffness degradation contribution from the fatigue loading in these cases were due to the interlaminar delamination. Additional validation cases are under investigations to account for higher stiffness degradations. This brings to the next challenges to extend this methodology to a complete intralaminar and interlaminar fatigue damage solution for variable amplitude and multi-axial loadings. REFERENCES [1] Sevenois, R.D.B., Van Paepegem, W., Fatigue damage modeling techniques for textile composites: review and comparison with unidirectional composite modeling techniques, Appl. Mech. Reviews, 67 (2015). [2] Van Paepegem W., Development and finite element implementation of a damage model for fatigue of fibre- reinforced polymers, PhD Thesis, Ghent University, Belgium, (2002). [3] Miner, M.A, Cumulative damage in fatigue, Journal of Applied Mechanics, 67 (1945) A159-A164. [4] Matsuishi, M., Endo, T., Fatigue of metals subjected to varying stress, Presented at Japanese Society of Mechanical Engineers, Fukuoka, Japan, (1968). [5] Dowling, N. E., Fatigue failure predictions for complicated stress-strain histories, Journal of Materials, JMLSA, 7(1) (1972) 71-87. [6] Brokate, M., Dressler, K., Krejci, P., Rainflow counting and energy dissipation in elasto-plasticity, Eur. J. Mech. A/Solids, 15 (1996) 705-737. [7] Nagode, M., Hack, M., The damage operator approach, creep fatigue and visco-plastic modeling in thermo- mechanical fatigue, SAE International Journal of Materials & Manufacturing, 4(1) (2011) 632-637. doi:10.4271/2011- 01-0485. [8] Nagode, M., Hack, M., Fajida, M., Low cycle thermo-mechanical fatigue: Damage operator approach, Fatigue Fract Engng Mater Struct, 33(3) (2010) 149-160. [9] Šeruga, D., Hack, M., Nagode, M., Thermomechanical Fatigue Life Predictions of Exhaust System Components, MTZ worldwide, 77(3) (2016) 44-49. [10] Brune, M., et al., FEM based durability analysis of the knuckle of the 5 series BMW, Fatigue Design’98, Helsinki (1998). [11] ASTM D3479, Space Simulation, composite mat., 15.03 (2015). T << /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