fast spectral rainflow fatigue damage assessment under wideband multipeak loading
DESCRIPTION
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. - PowerPoint PPT PresentationTRANSCRIPT
www.ifremer.fr
lfre
mer
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
www.ifremer.fr
lfre
mer
ContextFine description of sea states
climate as partition into waves systems
Complex fatigue damage assessment a large set of operational sea states (combinations of all the wave systems components) has to be considered
1
www.ifremer.fr
lfre
mer
metocean database
discretization+
statistics
scatter diagram (jPDF of HS, Tp, )
dimension = 3 102 to 103 fatigue
loadings
unimodal loadings
metocean database
discretization+
statistics
scatter diagrams (jPDF of HS
(i), Tp(i), (i) ,
i=1,2,3)
dimension = 9 106 to 109 fatigue
loadings.
partition
wave systems data (swells, wind sea)
multimodal loadings
2
www.ifremer.fr
lfre
mer
ObjectiveUse 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
3
www.ifremer.fr
lfre
mer
Outline
1. ICA FROMULA
2. APPLICATION: FPSO ON WEST AFRICA AREA
3. CONCLUSION & PERSSPECTIVES
4
www.ifremer.fr
lfre
mer
BASIC IDEA UNDER ICA FORMULA
sum signal
LF signalMaxA ; minA
MaxB ; minB
Partition of the set of turning points
LF narrow band
Miner damage with S-N curve: N = K Sm
HF & LF clearly separated
Assumptions
From the mathematical formulation of the rainflow counting [Rychlik, 1987], each subset is stable by rainflow counting
DS = DA + DB
1. ICA FORMULA
5
www.ifremer.fr
lfre
mer
EXPRESSION OF DAMAGES DA & DB
1. ICA FORMULA
HFm
aNA DkkD S
LH
S
AN N
N
N
Nk ;1
factor due to reduction in number of cycles
scaling effect due to addition of low-frequency signal
S
HFBa
BLFAASS
HF
Aa
M
MMk
MNMNMN
MM
Mk
11 :obtain we
from
mean value range :;
mBmLF
BM
K
ND
E2
(NB)
(NB = NLF)
1
22
m
K
ND
m
lfLF
LF
LF
mB
mB D
M
mD
E
12
1
1
hf
E[MBm] & E[MB] are estimated with a slightly conservative
approximation of the distribution of MB in terms of some spectral moments
6
www.ifremer.fr
lfre
mer
1. ICA FORMULA
Approximation of MB distribution
Ff ,The local maxima of the composite signal follow a Rice distribution [1945],
NS/NLH = 1/ : mean number of peaks between two successive zero-upcrossings. The peaks are assumed independent (conservative assumption)
sum signal
LF signalMaxA ; minA
MaxB ; minB
MB is the maximum among the local maxima between two successive zero upcrossings of the low-frequency signal
111
11 )()(
1)(,)()(
xFxfxyxFxY
7
www.ifremer.fr
lfre
mer
1. ICA FORMULA
Conservatism level
The conservatism level (Clevel) is higher for larger and lower For m = 3, Clevel > 130% for ; and ; For m = 5, Clevel > 130% for
8
0.1 mh,0/mt,0 0.9 ; 1.2 Tp,h/Tp;l 20 with Tp,h = 5s
+ m = 3+ m = 5
130%
www.ifremer.fr
lfre
mer
Briefly,
ICA is 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).
12
3
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).
1. ICA FORMULA
9
www.ifremer.fr
lfre
mer
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
2. APPLICATION
10
www.ifremer.fr
lfre
mer
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)
2. APPLICATION
Partitions of the sea states of the metocean databasecombinations number frequency [%]
MS only 1212 15,07MS + SS 2740 34,08MS + WS 1536 19,1
MS + SS + WS 2549 31,7no MS 3 0,05total 8040 100
11
www.ifremer.fr
lfre
mer
12
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
www.ifremer.fr
lfre
mer
Comparison of statistical properties of HS under assumption (H1)
metocean databasere-built metocean climate
13
2. APPLICATION
stats. prop.
mean HS 1,59 0,99 1,4 1,27 1,25 1,2 1,22 1,43
std HS 0,48 0,49 0,4 0,5 0,32 0,46 0,32 0,47
HS,1/3 2,14 1,25 1,85 1,77 1,6 1,46 1,58 1,92
HS,1/10 2,63 1,75 2,23 2,15 1,9 1,9 1,91 2,28
MS only MS + SS MS + WS MS + SS + WS
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).
www.ifremer.fr
lfre
mer
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
2. APPLICATION
stats. prop.
mean HS 1,59 1,64 1,4 1,32 1,25 1,27 1,22 1,15
std HS 0,48 0,5 0,4 0,44 0,32 0,36 0,32 0,35
HS,1/3 2,14 2,25 1,85 1,77 1,6 1,46 1,58 1,48
HS,1/10 2,63 2,75 2,23 2,15 1,9 1,9 1,91 1,79
MS only MS + SS MS + WS MS + SS + WS
Comparison of statistical properties of HS under assumption (H2)
metocean databasere-built metocean climate
www.ifremer.fr
lfre
mer
ICA-BASED DAMAGE ASSESSMENT PROCEDURE
1. Partition & Discretization of the wave spectra of the metocean database
2. Determination of the wave systems combinations probabilities of occurrence
3. Calculation of the responses spectra to the wave systems components
4. Responses partition and identification of the « response systems »
5. Calculation of the damages and other characteristics (spectral moments) of the response systems
6. Combination within each possible sea state (i.e. each wave systems combination) using ICA formula
7. Summation with the probabilities of occurrence
15
2. APPLICATION
www.ifremer.fr
lfre
mer
valeur Ratio ./Dt valeur Ratio ./Dt valeur Ratio ./Dt
MS only 15,07 0,2863 0,2684 93,75% 0,2557 89,31% 0,1492 52,11%MS + SS 34,08 0,2542 0,2263 89,02% 0,2314 91,04% 0,1618 63,66%MS + WS 19,01 0,0604 0,0548 90,73% 0,0567 93,89% 0,0438 72,55%
MS + SS + WS 31,7 0,0748 0,0613 81,95% 0,0692 92,57% 0,0547 73,17%total 100 0,6757 0,6108 90,40% 0,6358 94,09% 0,4772 70,62%
Dt* DHS,1/3 DHS,1/10sea states Prob. Occ. [%] Dt
Fatigue damage from metocean database (Reference)
* Partition with « triangle » spectral shape model
16
2. APPLICATION
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.
www.ifremer.fr
lfre
mer
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
2. APPLICATION
valeur Ratio ./Dt* valeur Ratio ./Dt*
MS only 15,07 0,2684 0,0808 30,11% 0,3615 134,70%MS + SS 34,08 0,2263 0,2483 109,73% 0,217 95,91%MS + WS 19,01 0,0548 0,1005 183,42% 0,0721 131,52%
MS + SS + WS 31,7 0,0613 0,233 380,06% 0,0638 104,15%total 100 0,6108 0,6626 108,49% 0,7145 116,98%
sea states Prob. Occ. [%] Dt*DICA - (H1) DICA - (H2)
Fatigue damage from ICA-based assessment method
* Partition with « triangle » spectral shape model
www.ifremer.fr
lfre
mer
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
www.ifremer.fr
lfre
mer
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
www.ifremer.fr
lfre
mer
Thank you for your attention !!!