1 / 21

FAST SPECTRAL RAINFLOW FATIGUE DAMAGE ASSESSMENT UNDER WIDEBAND MULTIPEAK LOADING

FAST SPECTRAL RAINFLOW FATIGUE DAMAGE ASSESSMENT UNDER WIDEBAND MULTIPEAK LOADING. Michel OLAGNON & Zakoua GU É D É IFREMER, Brest (France). Fatigue Workshop - 24 th February 2010 Broadband spectral fatigue: from gaussian to non-gaussian, from research to industry.

argus
Download Presentation

FAST SPECTRAL RAINFLOW FATIGUE DAMAGE ASSESSMENT UNDER WIDEBAND MULTIPEAK LOADING

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. FAST SPECTRAL RAINFLOW FATIGUE DAMAGE ASSESSMENT UNDER WIDEBAND MULTIPEAK LOADING Michel OLAGNON & Zakoua GUÉDÉ IFREMER, Brest (France) Fatigue Workshop - 24th February 2010 Broadband spectral fatigue: from gaussian to non-gaussian, from research to industry

  2. Fine description of sea states climate as partition into waves systems Context Complex fatigue damage assessment a large set of operational sea states (combinations of all the wave systems components) has to be considered 1

  3. unimodal loadings multimodal loadings metocean database metocean database partition wave systems data (swells, wind sea) discretization + statistics discretization + statistics scatter diagram (jPDF of HS, Tp, q) scatter diagrams (jPDF of HS(i), Tp(i), q(i) , i=1,2,3) dimension = 3 102 to 103 fatigue loadings dimension = 9 106 to 109 fatigue loadings. 2

  4. Use the Iterative Components Addition (ICA) formulas set up earlier [Olagnon & Guédé, 2008] to simplify the damage computation. ICA formulas allow to compute the damage of a multimodal loading spectra in terms of the damage of its individual components taken separately, keeping a low level of conservatism. D( S1 + S2 + S3 … ) = ICA( D(S1) , D(S2) , D(S3) …) Full damage computation [Stress – Rainflow – Miner] needed only for the wave systems components and ICA-based damage computation for the large set of all their combinations by ICA formula. drastic reduction of computation time Objective 3

  5. Outline ICA FROMULA APPLICATION: FPSO ON WEST AFRICA AREA CONCLUSION & PERSSPECTIVES 4

  6. 1. ICA FORMULA MaxB; minB MaxA; minA LF signal sum signal Assumptions Partition of the set of turning points • HF & LF clearly separated • LF narrow band • Miner damage with S-N curve: N = K Sm From the mathematical formulation of the rainflow counting [Rychlik, 1987], each subset is stable by rainflow counting DS = DA + DB BASIC IDEA UNDER ICA FORMULA 5

  7. 1. ICA FORMULA scaling effect due to addition of low-frequency signal factor due to reduction in number of cycles (NB) (NB = NLF) EXPRESSION OF DAMAGES DA & DB E[MBm] & E[MB] are estimated with a slightly conservative approximation of the distribution of MB in terms of some spectral moments 6

  8. 1. ICA FORMULA MB is the maximum among the local maxima between two successive zeroupcrossings of the low-frequency signal MaxB; minB MaxA; minA LF signal sum signal The local maxima of the composite signal follow a Rice distribution [1945], Approximation of MB distribution NS/NLH = 1/b : mean number of peaks between two successive zero-upcrossings. The peaks are assumed independent (conservative assumption) 7

  9. 1. ICA FORMULA Conservatism level 0.1  a=mh,0/mt,0  0.9 ; 1.2  x =Tp,h/Tp;l  20 with Tp,h = 5s x +m = 3 +m = 5 130% The conservatism level (Clevel) is higher for larger xand lower a For m = 3, Clevel > 130% for a = 0.1 ; x > 5 and a = 0.2 ; x > 7 For m = 5, Clevel > 130% for x > 7. 8

  10. 1. ICA FORMULA 3 2 1 Briefly, • ICAis simple and depends only on the main spectral properties of the signal components, namely their variance, irregularity factor and zero-crossing period(m0,m2,m4). • ICA can be used recursively when the signal has more than two modes (no narrow-band assumption on high-frequency component). • ICA becomes inconvenient with two-slopes S-N curve.Nevertheless, a formula which allows to obtain the two-slopes S-N damage in terms of the simple S-N has been set up (Olagnon, 2009 – in submission). 9

  11. 2. APPLICATION • Industrial application • Wave loadings: from metocean database of measured sea states in the West Africa area (almost 3 years) • Mechanical modeling: Wave bending moment of an FPSO Hull girder  linear hydrodynamic response (with RAO) • Fatigue design requirement: double-slopes S-N curve from Bureau Veritas requirements, 100 years design lifetime Work • Assessment of the occurrence probabilities of all the operational sea states given under some assumptions • Assessment of the total damage • from the metocean database with time-domain simulations • from the operator’s specifications with a ICA-based method 10

  12. 2. APPLICATION • Partitions of the sea states of the metocean database Metocean database partition & assumptions on metocean climate • Assumptions made for the re-construction of the metocean climate • (H1) : independent wave systems • (H2) : wave systems conditioned on the type of combination • Including a requirement to meet the proportions observed in the database and exclusion criteria on the ratio of the peak periods and the discrepancy between the directions (e.g. close wave systems components, unrealistic combinations) 11

  13. 2. APPLICATION Joint occurrence probabilities 1. Directional scatter diagrams of the Wave systems • Suit the discretization to the criteria: • log-scale (multiplicative) classes for the periods • half criterion threshold as class width for directions and the log(period). 2. Joint occurrence probabilities of all the possible combinations Given under independence assumption and truncated to the set of possible combinations. product of the wave systems frequencies from the scatter diagrams normalized by the sum of the frequencies of the possible combinations As a result, observed MS fall into 169 classes, SS into 148, WS into 165. (H1) :  2 millions combinations allowed to occur (H2) :  800 thousand combinations allowed to occur 12

  14. 2. APPLICATION Comparison of statistical properties of HS under assumption (H1) metocean database re-built metocean climate Significant discrepancies between statistics derived from the hypothesis (H1) and that of the database. The fact that the HS are higher for MS only is not reflected by (H1). 13

  15. 2. APPLICATION Comparison of statistical properties of HS under assumption (H2) metocean database re-built metocean climate Smaller discrepancies between statistics derived from the assumption (H2) and that of the database. However, under (H2) the HS are greater for combinations without SS and lower for combinations with SS. Both assumptions do not reflect the trends observed in the metocean database. 14

  16. 2. APPLICATION ICA-BASED DAMAGE ASSESSMENT PROCEDURE • Partition & Discretization of the wave spectra of the metocean database • Determination of the wave systems combinations probabilities of occurrence • Calculation of the responses spectra to the wave systems components • Responses partition and identification of the « response systems » • Calculation of the damages and other characteristics (spectral moments) of the response systems • Combination within each possible sea state (i.e. each wave systems combination) using ICA formula • Summation with the probabilities of occurrence 15

  17. 2. APPLICATION Fatigue damage from metocean database (Reference) * Partition with « triangle » spectral shape model • The sea states with higher HS are responsible for a greater proportion of the damage. DHS,1/3 represents 95% of the damage and DHS,1/10 represents 70%. But this result depends on the structural response under consideration. • The triangle spectral shape considered here (which does not take into account the tail of the spectra) yields a non-gaussian process and gives lower damages. A « wallops » spectral shape model (with a high frequency tail) would provide larger results. 16

  18. 2. APPLICATION Fatigue damage from ICA-based assessment method * Partition with « triangle » spectral shape model • The damage computed with the ICA-based procedure suit the effects of the assumptions made for the metocean climate re-construction. • Concerning the computation efficiency. The step n°7 of the ICA-based procedure last some seconds with a FORTRAN program on a fast computer ( 15s for 2 millions combinations). 17

  19. CONCLUSION • ICA formula is simply implemented in a general procedure for fatigue damage assessment under multimodal sea states with good performance; • This formula provided reasonably conservative results for the actual application. • Problems highlighted: • The use of two-slopes S-N curve. Problem solved using a formula to obtain the damage with double slope S-N curve from those with single slope S-N curve. • the metocean climate constructed under the independence assumption is not satisfactory. (H1) shows major discrepancies with the database. (H2) is better but not totally satisfactory. 18

  20. PERSPECTIVES • In the short term, other applications to be considered : • a very low-frequency response in addition to the wave systems responses • structural responses which are significantly direction-dependent • non-linear responses which can be represented by a linear response corrected with a factor depending on some parameters (e.g. HS,Tp) • In a longer term, extend the approximation of the distribution of MB to non-gaussian responses. • Improve the re-construction of the metocean climate, using instead on the overall statistics of environmental parameters, their evolution over time (events statistics). 19

  21. Thank you for your attention !!!

More Related