MM
E L S E V I E R
Available online at www.sciencedirect.com
c ScienceDirect
Applied Ocean Research 29 (2007) 221-230Applied O c e a n
Research
www.elsevier.com/locate/apor
Extreme response prediction for nonlinear floating offshore structures by
Monte Carlo simulation
A. Naess^'*''*, O. Gaidai^ P.S. Teigen'^
^Centre for Ships and Ocean Structures, Norwegian University of Science and Technology, NO-7491 Trondiieim, Norway ^ Department of Matiiematica! Sciences, Norwegian University of Science and Technology, NO-7491 Trondheim, Norway
^ StaloilHydro Research Centre, NO-7005 Trondheim, Norway
Received 13 April 2007; received in revised form 8 October 2007; accepted 24 December 2007 Available online 8 February 2008
Abstract
T l i e paper describes a method f o r the prediction o f extreme response statistics o f floating o f f s h o r e structures subjected to random seas by M o n t e Cario simtilation. T h e partictilar case o f the horizontal surge motions o f a tension leg p l a t f o r m is considered, taking- into account both the first order, wave frequency and the second order, s l o w - d r i f t morions. The advantage o f the M o n t e Cario method is its s i m p l i c i t y and versatility, w h i c h allows us to account f o r the effect o f time-variant w a v e - d r i f t damping, as w e l l as nonlinear m o o r i n g characteristics w i t h o u t noticeable increase in the computational complexity. It is demonstrated in this paper that the c o m m o n l y assumed obstacle against using the M o n t e Carlo method f o r estimating extreme responses, i.e. excessive C P U time, can be circumvented, b r i n g i n g the computation time d o w n to quite acceptable levels. © 2008 Elsevier L t d . A l l rights reserved.
Keywords: Monte Carlo simulation; Mean crossing rate; Extreme response prediction; Tension leg platform
1. Introduction
The problem of predicting the extreme response statistics of a floating offshore structure subjected to random waves by Monte Carlo sirnulation is the challenge we address in this paper. Specifically, we shall discuss the wave forces and motion response of a tension leg platform (TLP), with particular focus on the estimation of the total extreme surge response in a random sea way. The state-of-the-art formulation of the surge response of a TLP in random seas is in terms of a second order stochastic VolteiTa series. Such models and their staristics have been studied over several decades [1-10]. In the case of a linear dynamic system there is a possibility of accurately estimating the extreme response statistics [11-13]. In these papers a numerical calculation technique based on characteristic functions and the so-called saddle point method is developed. I f the dynamic system is nonlinear, the saddle point method is not applicable. In [14] a numerical path integration method was elaborated for a dynamic system with
* Corresponding address: CeSOS/NTNU, Otto Nielsens vei 10, NO-7491 Trondheim, Norway.
E-mail address: arvidn@math.ntnu.no (A. Naess).
nonlinear damping and nonlinear' stiffness, but only the slow-drift response could be calculated.
In principle, such limitations do not exist for the Monte Carlo method. However, the standard Monte Carlo method has not been very efficient for estimating extreme response statistics. Accurate estimates of events with a very low probability of occurring, which is typical for extreme responses, require very long simulation times. We shall show in this paper that this obstacle can be effectively circumvented by adopting an appropriate extrapolation technique obtained by exploiting the properties of the mean upcrossing rate function, which, under nonrestrictive assumptions, completely determines the extreme value distribution. It will be shown by numerical examples that by combining empirical estimates of the mean upcrossing rate obtained from the standard Monte Carlo method with an appropriate linear extrapolation technique, a dramatic reduction of the CPU time needed to obtain accurate estimates of the extreme response distribution can be achieved. In fact, this reduction in the required CPU time makes the standard Monte Carlo method available in practice for the accurate estimation of extreme responses with a modern laptop computer.
222 A. Naess et al./Applied Ocean Research 29 (2007) 221-230 Over the years, several simulation based approaches ibr
the estimation of extreme response of nonlinear floating structures subjected to random waves have been proposed. For a discussion, the reader is refeiTed to [15]. The work described in the present paper, which follows a different approach from the methods described in [15], extends some previous work by the authors reported in [16].
2. Tlie dynamic model
The equarions of motion for a floating, rigid-body structure subjected to environmental forces like wind, waves and cuiTcnt would generally be written as
M Z ( 0 + H ( Z ( 0 , Z ( r ) , 0 = F(/), (1) Here M denotes a generalized 6 x 6 mass matrix, Z =
( Z i , . . . , Z 5 ) ^ = the structure's response vector, while H denotes a nonlinear vector function. F ( 0 denotes a stochastic loading process, which in general also depends on the response of the structure. This point will be discussed at some length below.
The methodology developed in this paper applies directly to a dynamic model as expressed in Eq. (1). The only requirement is that the response time histories can be simulated. Since the main purpose of this paper is to illustrate the versatility and accuracy of the proposed method, we have chosen to discuss a simplified single-degree-of-freedom (SDOF) model for the surge response of a tension leg platform (TLP) in random waves. Except for interaction effects between different motion modes, the SDOF model allows for the introduction of most of the relevant nonlinear effects that may influence the response of the TLP. Hence, in this paper the following SDOF equation of motion will be studied
MZit) + H{Z{t),Zit),t) = Fit), (2)
where Z — Z ( f ) denotes the surge response of the TLP, M is the mass of the platform, including added mass, and W is a nonlinear function to be specified.
The hydrodynamic loading process F{t) will be assumed to consist of two components: A linear, first order (wave frequency) term Fi{t), and a nonlinear, second order term
F2{t). The linear, first order wave excitation in short-crested
random seas in deep water can be expressed as a stochastic process F\ (f) as follows
r o o r n
F i ( 0 = / / ki(r; p)dpXit -r)dr (3) where d f j X ( t ) denotes the incremental directional components
adding up to give the short-crested random seas X{t), which represent the wave elevation at some reference point. Formally,
X ( 0 = f^„dpX(t), cf. [6]. X{t) is assumed to be a stationary zero-mean Gaussian process describing a short-term sea state. ki(r; P) denotes the impulse response function of the linear transfer mechanism between a long-crested wave field propagating in the direction fi and the surge excitation forces.
Similarly, the nonlinear, second order wave excitation forces in surge can be expressed as
POO p c o rJT (•JJ
Flit) = / / / / k2(T,r':p,ft')df,X(t-r) Jo Jo J-JJ J-JJ
x d / j ' X ( f - r ' ) d r d r ' (4) where kjii, r'\ p, P') is called the bidirectional quadratic
impulse response function for surge. The total surge excitation forces Fit) on the structure that are commonly accounted for is then Fit) = F\ it) + FiQ), which is a second order stochastic Volterra series model. This representation gives a description of the hydrodynamic forces accurate to second order on a large volume structure, that is, a case where inertia forces dominate and the potential theory is applicable. However, there is one aspect of the second order forces on a floating structure that has turned out to be of some significance. This is related to what is known as wave-drift damping, which can be attributed to a dependence of ^2(0 on the motion of the TLP. A discu!5sion of this is given below.
To be able to express the force process Fit) through the random directional sea way as expressed in Eqs. (3) and (4), it is required to provide the relevant hydrodynamic transfer functions for the surge mode. As already stated, for simplicity of exposition, coupling with other motion modes is neglected. However, the method we present in this paper is equally suitable for a fully coupled analysis. Transfer functions are calculated by using commercial software, e.g. WAMIT. This program was used to calculate the (directional) linear and the (bidirectional) quadratic (surge force) transfer functions tor the TLP. These are denoted by K] io); P) and K2ia), w'; p, /^i'). I'espectively. These two functions are calculated on a discrete grid in parameter space. Specifically, let 0 < C/JQ < • • • < OJL+\ be the chosen discretization of the (positive) part of the frequency axis, and —TT < Po < • • • < PM+\ < ^ the corresponding discretization
of the direction angle range i—7T, n).
The spectral representation of the short-crested random wave elevation ^ ( 0 at the chosen I'eference location is then expressed as follows
L M n
= E E V ? ^ ^ ( 1 ^ / 1 ' Pk)^cojAPk Bjk e'^>'. (5)
j=-Lk=\ '
Here SxioJj; Pk) denotes the directional one-sided spectral density of the shortcrested random seas. Awj = icoj+\
-Wj-\)ll and APj = iPj^\ - Pj-\)/2. oj-j = -ojj. {Bjk}, j — \,...,L, k = I , . . . , yW, is a set of independent,
complex Gaussian N(0,1)-variables with independent and identically distributed real and imaginary parts. These variables also satisfy the relations S . / .A = S*^,, where * signifies complex conjugation. This provides the appropriate extension to negative frequency indices. Also note lhat throughout this paper, when the summation index runs from negative to positive values, it invariably omits zero.
A. Naess el al/Applied Ocean Research 29 (2007) 221-230 223
introduce the slow-drift approximation. This amounts to putting
KiicD, -co'; p, P') = 0 for co • w' < 0. (6)
This condition clearly implies that all sum-frequency compo-nents of the second order excitation forces are neglected, re-taining only the difference frequency components, which are the only terms tiiat can give rise to the slowly-varying response.
In order to set up the proper dynamic model for the surge response Z(t), it is necessary to take into account the fact that hydrodynamic loading on a floating body depends on the motions of the body. In case of the slow-drift motions of the TLP, it has turned out that it is of some iinportance to take into account the dependence of the slow-drift force ^ 2 ( 0 on the slowly-varying surge velocity Z 2 ( f ) . Since we have a nonlinear dynamic model, we shall have to introduce a definition of the slow-drift response Z2(t). A suitable definition would seein to be the following: The slow-drift response Z 2 ( / ) is obtained from the total response Z(t) by a low-pass filter that removes all wave frequency components. In practice, this may be achieved by a running mean operator used iteratively.
To account for the dependence of the slow-drift force on the slowly-varying velocity, it is appropriate to write
F2{t,Z2{t)) rather than ^ 2 ( 0 . However, F2{t,Z2{t)) is not
directly available to us, but only F2(r) = F2{t, 0) given by Eq. (4). Since, in the context of slow-drift motions, Z 2 ( f ) is small, the following approximation is adopted,
9^2(^,0) .
F2{t, Z 2 ( 0 ) ^ F2{t, 0) -I- :\' Z2{t). (7)
0 Z 2
It is shown in [14] that for a TLP structure
-cF2(t,0) = -c F2(t) for a suitable constant c > 0, may
to some extent serve as a useful approximation to capture qualitatively the time-variant damping effect, which is the result of the expansion in Eq. (7).
The first dynamic model adopted for Z(t) is now the following,
MZ + CZ + K(Z + EZ^) = F| (f) -f- F2it, ±2)
^ Fit) - cF2{t) Z2. (8)
Here M = M + m, where m is an appropriately chosen (constant) added mass. C, K and e are suitably chosen positive constants. This equation is rewritten in the equivalent form
Z + 2cOe^Z-\-2cOeCF2{t)Z2 + wl{Z + EZ'^) = (9) M
where = K/M, ; = C/(2a)eM), c = c/{2co^lVl).
Thus, the dynamic system is nonlinear with a Duffing-type hardening stiffness nonlinearity and time-varying damping. Since for the TLP the relative damping coefficient f is usually small, the contribution from the time-varying term is nonnegligible, especially for severe seas for which the slow-drift response is significant. The third order term in the restoring force, generally refeiTcd to as the "set down" effect, is caused by the fact that the tethers will induce the TLP to act like an inverted pendulum. To lowest order, the surge restoring force
is S sin d = Sj^, where S is the total tension in the tethers, 6 is the angular displacement corresponding to a surge displacement Z, and L is the effective length of the tethers. However, due to the pendulum effect, a surge motion will also involve a vertical displacement by an amount L ( l - cosÖ). Given that the pontoons are fully subinerged, this will lead to an additional restoring force AS, so that the total surge restoring force equals ( ^ • - f / i 5 ) s i n 0 = S&me + pnD'^L{\ - c o s 0 ) s i n 6 i
The above equation gives a reasonable representation of the surge restoring force to third order in Z provided tether elongation can be ignored and the tethers can be inodelled as straight lines. This is, however, not quite coirect and additional rnodifications of the restoring term would result by considering both tether elongation, tether dynamics and weight/buoyancy effects. Note that the set-down effect will have an influence also on the hydrodynamic loading process, which becomes dependent not only on Z but on Z as well. Even i f this dependence could have been taken into account, it has been neglected in this paper as it is of minor importance.
For a TLP subjected to wave action only (i.e. in the absence of wind and current), Eq. (9) can be expanded further in order to account for other important nonlinear effects. The most important further extension of the equation of rnotion can be obtained by including viscous drag on the hull. The most basic way to accommodate a drag term is to include an expression of the type
FD=0.5PACD\Z\Z, (11) where A denotes the projected area of the TLP on a plane
normal to the wave direction, p is the density of the sea water, and Co the nondiinensional drag coefficient.
This expression has the apparent weakness that it only accounts for the velocity of the hull relative to still water, and only indirectly (through the drag coefficient) accounts for . shielding effects and diffraction. A more appropriate force term
would result i f one includes the relative velocity between the hull and the fluid and accounts for diffraction, that is,
FD = 0.5p ^i^Dd I Z - Ui\{Z - Ui), (12)
where the summation index runs over different component members of the submerged part of the TLP, and where the fluid velocities t// in the surge direction (x-direction) are computed at the location of the structural component, and in the absence of that pardcular component.
224 A. Naess et al. /Applied Ocean Research 29 (2007) 221-230 The drag coefficient can be estimated on the basis of model tests
or obtained by superposition of the drag force on individual iTiembers. In this case shielding and diffraction can only be incorporated in an approximate sense. An argument for introducing this sirnplification is that for the heavy sea states considered in this paper, it is a reasonable approximation to assume that the wavelengths are "long" relative to the dimensions of the TLP. In such a case, the phase differences in the velocity field over the TLP structure may be effectively ignored. Furthermore, the cross flow principle can be assumed to be valid for the individual structural members of the TLP (i.e. the vertical TLP columns and the horizontal pontoons). For the one-dimensional case considered in the present paper the "drag-summation" thereforej-uns over the columns and the fore and aft pontoons, and include the x-velocities in the fluid domain only. Strictly speaking, the latter choice is not fully consistent with the cross flow principle as the j-component of the drag force on the pontoons would be proportional to (Z —
f/,.) ^ ( Z - [/,.)2 + U l and to ( Z - U,) ^{Z - U,)^ +
for the columns. However, the main uncertainty here lies with the actual magnitude of the drag coefficients, and discarding the xyz coupling may therefore seem like a small sacrifice, compared to the reduction achieved in the coinputational work. Allowing for a certain amount of linear damping the proposed equation of motion then takes the form:
MZ + CZ + cF2(t)Z2 +
J2
CD, . / | Z - UjliZ - U j)+ K{Z + EZ^) = Fit), (13)
where Cpj now are dimensional drag terms, proportional to the total projected surface area of the structural members lumped together into group no. j, and where Uj = Uxj denotes the water particle velocity in the surge direction at a suitable reference point in group no. j.
3. Extreme value prediction
There is in general no simple procedure for estimating the exact statistical distribution of the extreme response values, However, for the applications we have in mind in this paper, it is possible to obtain a very good approximation, which can be expressed in terms of the so-called mean upcrossing rate. To explain how this is achieved, let N~^i^; t\, t2) denote the random number of times that the process Z ( 0 upcrosses the level f during the time interval (fi,??)- By introducing the mean rate of ^-upcrossings of Z ( f ) at time t, denoted by v'^i^;t), then under suitable regularity conditions on the response process, which can be adopted here, the following formula is obtained
v+i^; t) = lim E[N+i^;t - At/2,t + At/2)] At
f
Jo
/ z ( / ) Z ( / ) ( t . 5 ) d ^ (14)
where /z(,)z(/)(''') denotes the joint PDF of Zit) and Z ( / ) =
dZit)/dt, and E[-] denotes expectation. Eq. (14) is often
referred to as the Rice formula [18]. u+(§; 0 is assumed throughout to be finite.
Let MiT) = max{Z(0 : 0 < / < 7 ) denote the extreme value of the surge response process Z ( 0 over the long-term time interval of length T. The cumulative distribution function of M ( r ) under the assumption lhat the stream of upcrossing events constitute a Poisson process, is then given in terms of the mean upcrossing rate as follows [19].
P r o b ( M ( r ) < ?) = Prob(A' + (?; 0, T) = 0) = e x p { - E [ A ' + ( ^ ; 0 , T)]
'+i^;t)dt\.
exp i\5)
where v+(^; l) denotes the mean level upcrossing rate at time /. For practical purposes, this is rewritten as [19,20]
Prob(/V/(r) < ^ ) = exp -T / v+il^:u))fw(w)dw , (16) where /vi/(u;) denotes the long-term (ergodic) PDF of relevant sea state parameters w, and u + d ; w) equals ihe mean upcrossing rate for the short-term (stationary) condition specified by w. In our case, W = (//,,., T/,), where Hg denotes the significant wave height and Tp is the spectral peak period of a short-term sea state. This implies that the long-term PDF can be estimated from the scatter diagram of flie sea states at the specified location.
The practical long-term analysis would typically be carried out by weighting the crossing rate function over a selected group of sea states that is believed to provide a sufficient resolution over the scatter diagram. For the extreme responses, small to moderate seas do not contribute noticeably, despite their dominating rate of occurrence. Therefore, it is the large to severe sea states that are of particular importance, reducing substantially the number of sea states that need to be included. Eq. (16) can then be written as
Prob(A'/(r) < ^ ) = exp (-J^Ck v+(?, A-) T (17)
k€lC
where /C denotes the set of included sea states parameterized by the index k. u+(§, k) denotes the mean rate of upcrossings of the level ^ by the response process X^it) for sea state k, and
Ck is its relative frequency of occurrence given by the scatter
diagram. Time T may be e.g. a 20 year period, expressed in seconds T ^ 6.3 x 10^ s.
4. Empirical estimation of the mean crossing rate
A. Naess et al. /Applied Ocean Research 29 (2007) 221-230 225 the mean upcrossing rate is conveniently estimated from the
ergodic mean value. That is, it can be assumed that lim » + ( g ; 0 , T)
T (18)
where 7i+(§; 0, T) denotes a realizadon of A'+(^; 0, T), that is, 71+(^; 0, T) denotes the counted number of upcrossings during time T froin a particular simulated time history. It is assumed that initial transient effects have been removed from these time series, that is, only the stationary parts of the time series are included. In practice, m time histories of a suitable, specified length, To say, are provided by simulation. The appropriate ergodic mean value estimate of v+(§) is then
(19)
where ; i t (^; 0, TQ) denotes the counted number of upcrossings of the level ^ by tiine history no. j. This will be the approach for the estimation of the mean upcrossing rate adopted in this paper.
For a suitable nuiuber e.g. m > 20 — 30 [21], and provided that TQ is sufficienüy large, a fair approximation of the 95 % confidence interval CIo.95(^) for the value v"'"(^) can be obtained as
CIo.95(?) = v+(?) ± 1.965(^)/V^, (20) where the einpirical standard deviation s{^) is given as
in
'n + (^;0,To)
To 0+(?)
(21)
It is pointed out that the choice of in and TQ for a specific case is mainly determined by an optimization of the estimate given by Eq. (21). Hence, the initially siinulated time series may be divided into shorter blocks for the sake of optimization of this estimate. Note that the point estiinate 0+(^) is not affected by this operation. The consistency of the estimates obtained by Eq. (21) can be checked by the observadon that Var[7V+(^; 0, r ) ] = v+(^)T since N+(^; 0,T) is a Poisson random variable by assumption; here Var[-] denotes variance. This leads to the equafion s{^)^/in = v'^{^)/{ni To). The empirical counterpart of Eq. (21) i f the assumption that
N'^i^; 0, T) is a Poisson randotn variable is explicitly used, is
then s{^)'^/m ^ v+(^)/{m To). It is worth poinfing out that this last estimate is invariant with respect to blocking of the data, hence the optiinization step is avoided.
The approach to extreme value prediction presented in this paper derives froin the realization that for the dynamic models of the type considered, the short-term as well as the long-term mean ^-upcrossing rate as a function of the level ^ is in general highly regular in a specific way. In fact, it has been observed that for a wide range of dynatnical systeins the mean upcrossing rate tail behaves with good approximation like e x p ( - f l (^ - b)'^] (? > ?o) where a, b and c are suitable constants. Hence, as discussed in detail in [22], it may be assumed that
v+(?) ~ ^ ( ^ ) e x p ( - f l ( f - b)'}, ^ > $0, (22)
where the funcdon q(^) is slowly varying compared with the exponential function e x p { - f l ( ^ - b)^} forfait values of | .
Therefore, under the assumptions made, log | log(y+(^)/ '7(^))l plotted versus log(^ — b), exhibits an almost perfectly linear tail behaviour. I f the function q(^) could be replaced by a constant value, q say, which could be determined, we would imtnediately be in a position to apply a linear extrapolation strategy to find the values of v~^{^) at extreme levels of ^. In general, q(^) is not constant, but its variation in the tail region is usually sufficiently slow to allow for its replacement by a constant, cf. [22]. In practice, an initial choice of the parameter q is extracted from an average of the ratio
v'^(^)/fzi^) for large values of § [22]; here / z ( ? ) denotes the
response probability density function, which would then also be estimated empirically. Having chosen a value for q, the value of the parameter b must be based on opfimizing the linear fit, while at the same time keeping as small as possible. The last proviso comes from the desire to obtain accurate extrapolation, which is related to the accuracy in the slope of the fitted straight line. Hence the region over which the fitting is done, should be as long as the linear fit allows.
As mentioned above, and discussed in [22], the ratio
v'^(^)/fz{^) may vary with ^. The inifial choice of a value
for the paraineter q would consequently be somewhat arbitrary. It is therefore suggested to allow for optimization also with respect to the value of i f that is deemed necessary. However, in this paper q is taken as an average of the ratio f'*'(^)//z(?), empirically determined from the data tail, which seems to work well in niost cases.
The specific behaviour of v+(^) expressed by Eq. (22) is based on the underlying assumpdon that the appropriate asymptotic extretne value distribution for the stationary response data under study is the Gumbel distribution. For all relevant dynamic models for offshore structures as given by Eq. (1), with Fit) modelled as a stationary process, this would be the case. Therefore, the following asymptotic relation for the extreme value distribution is obtained [23],
P r o b ( M ( r ) < ?) ~ exp { - exp {-ap (? - br))}, T ^ oo, (23) for suitable parameters ap > 0 and
br-In practice, a commonly used approach is to assume that the observed extreiue value response data do follow a Gumbel distribution. The problein with this approach is that classical extreme value theory cannot be used to decide to what extent the asymptotic distribuUon is actually valid for a given set of extreme value data. Note that the asymptotic Gumbel distribution corresponds to an asymptotic upcrossing rate that is purely exponential, that is, with c = 1.0 in Eq. (22). Hence, by adopting a much more general class of functions, with the purely exponential functions as a sub-class, to represent the upcrossing rate, as we have done here, the ability to capture sub-asymptotic behaviour is greatly enhanced. By this, the necessity to adopt a stricUy asymptotic extreme value distribution, of questionable validity, is avoided.
226 A. Naess el al./Applied Ocean Research 29 (2007) 221-230 Table I Particulars of the TLP Column diameter D (m) 10.0 Eigenperiod surge (,s) 128.8 Relative damping f 0.05 Total mass (incl. added ma.ss) M (kg) 1 .."i x 10
specific extreme value esümation to be carried o u t with sufficient accuracy for most practical prediction purposes with inuch less numerical efforts than a direct estimation of the pertinent extreine values. Note that on the premise of a linear tail behaviour relative to a specified plotting procedure, the esümation uncertainty is related to the estimation of a slope. Hence, i f the confidence band is sufficienüy n a i T o w over a suitable range of linear behaviour, the uncertainty in the slope will tend to be correspondingly small, and consequently, so will the uncertainty of the estimated upcrossing rate at the extreme response levels. This is a very important feature of the proposed extreme value prediction method, since the assumption that linear tail behaviour can be achieved by the proposed method seems to have wide applicability [22,24]. However, in.practice, for short simulations the available data may not always allow for a very narrow confidence band, as will be illustrated in the next section. In [25] a general approach to the estimation of a confidence interval for the extrapolated estimate of the upcrossing rate at the extreme response level is proposed.
5. Numerical results
In this section empirical estimates of the mean upcrossing rate obtained b y Monte Carlo simulations combined with an opümized extrapolation procedure will be presented.
A particular model of an offshore tension leg platform (TLP) is considered and the c o i T c s p o n d i n g linear transfer funcdon (LTF) and quadratic transfer function (QTF) are computed using a second order diffraction program (WAMIT). For simplicity, unidirectional seas are used, meaning that the directional argument fi is skipped. This simplification should have no effect on the conclusions based on the comparison of accuracy. The combined first order and second order slowly-varying surge deck motion is studied applying a single-degree-of-freedom model. The TLP particulars are detailed in Table 1 and the sub-surface part of the structure is shown in Fig. 1.
The values in Table 1 are used to obtain the second order response. This means that for the second order response we have used a simplified version of Eq. (8), where mass
M, stiffness K and damping coefficient C are frequency
independent, which is a good approximation for the slow-drift motion. The time-invariant damping part ^ is considered to be 5%.
In this paper, two versions of Eq. (9) will be used. The first version is a linear, time-invariant model obtained by putting c = e = 0. The second version is the fully nonlinear model where the parameter c in (9) is chosen such that Var[2W(,cF2(0 2,2(1)] is about 10% of Var[2&)e ? Z ( 0 ] . The parameter s is estimated
fi-om the condition that 0 . 2 Z ( 0 > eZ^(t) w h e n Z{t) < 6az,
i.e. even in the extreme response region, stiffness hardening
20 20
Fig. 1. Sketch of submerged part of TLP. Units in meters.
contributes not more than 20% relative to the linear part for severe sea, which leads to e = 1.36 x 10"'*. Finally, the following approximate values were found, c = 30/(Mg), g = 9.81 m/s^ for moderate seas, and c — 90/(Mg) for severe seas. The adopted parameter values are largely arbitrary, but the choices made seem to provide a reasonable model for the chosen TLP structure.
To get an accurate representation of the response process, there is a .specific requirement that must be observed. Since the damping ratio is only 5%, the frequency resolution Aco must secure a sufficient number of frequency values over the resonance peak. This will ensure that the second order difference frequency response component captures the TLP surge dynamics with sufficient accuracy. It is commonly required that there are at least 5 discrete frequencies over the frequency range where \L(CD)\^ is equal to or higher than half of the resonance peak height niax(|L(w)|-), where
L(ü)) = (-ar + 2if OJ, W + oj~^ ' . (24)
For the surge force QTF A'2(w, co') a suitable initial frequency grid must be chosen for which the values of the force QTF are calculated. The calculation of the force QTF is generally by far the most time consuming part of the numerical analysis. Therefore the initial grid is usually rather coarse to avoid excessive computer time. For the calculations at hand the discrete frequency range was the following: co\ =
In/30.0, ...,co„ = 27r/4.0 (rad/.s), /( = 30. This necessitates
the use of an interpolation procedure to be able to provide values of the QTF on a much finer grid than the initial one to comply with the requirement of sufficient frequency resolution to capture the dynamics of slow-drift motion. In this paper cubic spline interpolation has been used. In the particular ca.se considered here, the resolution requirement led to the choice
Aw = 0.0018 rad/s and L = 760 interpolated discrete
A. Naess el al./Applied Ocean Research 29 (2007) 221-230 227
Table 2
Representative sea states along with response standard deviations, nonlinear and linear TLP
Hs (ni) Tp (s) az (ill). Ihi- az (ni), nonlin.
10.0 11 9.3 8,2
15.0 17 6.4 6.0
xlO^N/m Surge wave exciting force. Dimensional results for unit wave amplitude
2.5 2 1.5 -1 -0.5 n ] I 1 . . , , 1 0.2 0.4 0,6 0.8 1 1,2 1,4 1.6 m
Fig. 2. Wave exciting force amplitude, surge LTF,
Each random stationary sea state is specified by a JONSWAP spectrum, which is given as follows
where g = 9.81 m s"^^ c^^^ denotes the peak frequency in rad/s and where y and CT are parameters affecting the spectral shape. CT — 0.07 when to < cOp, and CT = 0.09 when co > cOp. The parameter y is chosen to be equal to 3.3. The parameter a is determined from the following empirical relationship
a = 5.06 (1 - 0.287 In y ) (26)
VP/
where Hg denotes the significant wave height and Tp = In/cOp is the spectral peak wave period. For the present analysis, two short-term JONSWAP sea states have been used — moderate and severe. Table 2 presents the sea state parameters along with the coiTcsponding response standard deviations.
Fig. 2 shows the LTF for the wave exciting force amplitude, while Fig. 3 depicts the spline interpolated QTF. Since the QTF is complex valued, only its absolute value is plotted.
For the two chosen sea states. Figs. 4-7 present the coiTcsponding response tail crossing rates obtained by Monte Carlo simulation for the linear system given by putting e = 0 and c = 0 in Eq. (9). Figs. 4 and 6 show the results obtained froin 1000 realizations, each of duration 1 h.
Fig. 3, Wave exciting force amplitude, surge QTF,
2 3 5 8,68 ija^ Fig, 4. Empirical crossing rates by Monte Carlo simulation (*) with 95% confidence bands ( ) based on 1000 h of response time histories for the case of linear dynamics (c = e = 0), Saddle point integration results (—) coincide with the optimized linear fit with b = 0.75(TZ, q = 0.0205, where
az = 9,3 m (see Table 2), Sea state with Hg = 10 m, Tp = ' 1 1 s.
The crossing rate plots are done on the transformed scale, i.e. log|log(u+(^)/5)|versus log((^ - b)/az), see Eq. (22). Extreme response prediction based on Eq. (17) will typically involve crossing rates of the orders 10~^ — 10~^, but to illustrate the achieved accuracy the extrapolated results at the crossing rate level 10""'° are highlighted. Figs. 4 and 6 also show the highly accurate results obtained by using a saddle point integration technique [13]. These results cannot be distinguished from those obtained by linear extrapolation of the mean upcrossing rate function provided by Monte Carlo simulations. Hence, linear extrapolation can be done over several orders of magnitude with high accuracy.
To illustrate the fact that good accuracy can be obtained with much shorter time simulation records. Figs. 5 and 7 show the results of using only 25 h of simulated response fime histories, which required about 10 min. of CPU time on a standard desktop PC.
228 A. Naess et al. /Applied Ocean Research 29 (2007) 221-230
Lin TLP, moderate sea H =10m, T =11s
s p Lin TLP, severe sea H 15m, T =17s p
Fig. 5. Empirical cro.s.siiig rates by Monte Carlo simulation (*) with 95% confidence bands ( ) based on 25 h of response time histoiies, and by opümized linear fit (—) with b = 0.75(TZ, q = 0.020, where az = 9.3 ni (.see Table 2), for the case of linear dynamics (c = s = 0). Sea .state with
Hs = 10 m, Ti,= \\ s.
Fig. 7. Empirical ci;ossing rates by Monte Carlo simulation (*) with 95% confidence bands ( ) ba.sed on 25 h of response time histories for the case of linear dynamics (c = 8 = 0); optimized linear fit (—) with b = 0.39(7z.
q = 0.042, where az = 6.4 m (see Table 2). Sea state with H, = 15 m.
17 s.
Lin TLP, severe sea H =15m, T =17s
s p Nonlin TLP, moderate sea H =10m, T =11s
s p
Fig. 6. Empirical crossing rates by Monte Cailo simulation (*) with 95% confidence bands ( ) based on 1000 h of response time hi.stories for the case of linear dynamics (c = s = 0). Saddle point integration results (—) coincide with the optimized linear fit with b = 0.39az, q = 0.042, where az = 6.4 m (see Table 2). Sea state with H,,
Fig. 8. Empirical crossing rates by Monte Carlo simulation (*) with 95% confidence bands ( ) based on 1000 h of response time histories for the fully nonlinear model; optimized linear fit (—) with h = 0.54az, q = 0.0248, where az = 8.2 m (see Table 2). Sea state with H, = 10 m, 1), = 11 s.
nonlinear model given by Eq. (9). In this case the Monte Carlo simulation results are the only results available for verification of the extrapolation method for the nonlinear tnodel, but the experience from the linear case indicates that the results obtained from Figs. 8 and 10 are very accurate. Again, the crossing rate plots are done on the transformed scale, see Eq. (22), and for illustration purposes the extrapolated results at the crossing rate level 10""^ are highlighted. These results can be compared with the corresponding results obtained from only 25 h of simulated response time histories shown in Figs. 9 and
11. The agreement is again good.
6. Concluding remarks
In this paper we have described a method for predicting the extreme response of floating marine structures in random waves based on Monte Carlo simulations. The prerequisite for application of the method is thus the availability of the equations of motion of the structure which makes possible the simulation of its response time liistories caused by stochasfic environmental loads.
Nonlin TLP, moderate sea H^=10m, Tp=11s
A. Naess et al./Applied Ocean Research 29 (2007) 221-230 229
Nonlin TLP, severe sea H =15m, T =17s
p
10
6.65 9
Fig. 9. Empirical ciossing rates by Monte Carlo simulation (*) with 95% confidence bands ( ) based on 25 h of response time histories for the fully nonlinear model; optimized linear fit (—) with b = 0.54az,q = 0.0234, where
cTz = 8.2 111 (see Table 2). Sea state with //., = 10 m, r „ = 11 s.
Nonlin TLP, severe sea H =15m, T =17s
s p
Fig. 10. Empirical crossing rates by Monte Carlo simulation (*) with 95% confidence bands ( ) based on 1000 h of response time histories for the fully nonlinear model; optimized linear fit (—) with b = 0.39crz, q = 0.045, where CTZ = 6.0 m (see Table 2). Sea state with « j = 15 m, 7^ = 17 s.
numerical results for the surge response of a TLP model have been presented. For two chosen sea states, complete agreement has been demonstrated between the predictions provided by the proposed method and an exact method for the case of linear dynamics. To verify the accuracy of the proposed method for the general case of nonlinear dynamics also, extensive Monte Carlo simulations were carried out for the same two sea states, and both cases lend further support to the applicability and accuracy of the method.
Acknowledgements
This work was supported by the Research Council of Norway (NFR) through the Centre for Ships and Ocean Structures (CeSOS) at the Norwegian University of Science and
Fig. 11. Empiilcal ciossing rates by Monte Carlo simulation (*) with 95% confidence bands ( ) based on 25 h of response time histories for the fully nonlinear model; optimized linear fit (—) with b = 0.39ffz. 9 = 0.045, where CTZ = 6.0 m (see Table 2). Sea state with Hs = 15 m, Tp = M s.
Technology. The authors would also like to acknowledge the support of the StatoilHydro Research Centre.
References
[ I ] Neal E. Second order hydrodynamic forces due to stochastic excitation. In: Proc. 10th ONR symposium. 1974.
[2] Vinje T. On the statistical distribution of second-order forces and motions. International Shipbuilding Progress 1983;30:58-68.
[3] Naess A. Statistical analysis of second-order response of marine structures. Journal of Ship Research 1985;29(4):270-84.
[4] Langley RS. A statistical analysis of low frequency second-order forces and motions. Applied Ocean Research 1987;9(3).
[5] Kato S, Ando S, Kinoshita T. On the statistical theory of total second order response of moored floating structures. In: Proc. 19th annual offshore technology conference. Paper No. 5579. 1987.
[6] Naess A. Statistical analysis of nonlinear, second-order forces and motions of offshore structures in short-crested random seas. Probabilistic Engineering Mechanics 1990;5(4): 192-203.
[7] Naess A, Johnsen JM. An efficient numerical method for calculating the statistical distribution of combined first-order and wave-drift response. Journal of Offshore Mechanics and Arctic Engineering 1992;114(3):
195-204.
[8] Langley RS, McWilliam S. A statistical analysis of first and second order vessel motions induced by waves and wind gusts. Applied Ocean Research 1993;15(l):I3-23.
[9] McWilliam S, Langley RS. Extreme values of first and second order wave induced vessel motions. Applied Ocean Research I993;15(3):169-8I. [10] Grime AJ, Langley RS. On the efficiency of crossing rate predicfion
methods used to determine extreme motions of moored offshoie structures. Applied Ocean Research 2003;25(3): 127-35.
[11] Naess A. Crossing rate statistics of quadratic transformafions of gaussian processes. Probabilistic Engineering Mechanics 200I;16(3):209-17. [12] Naess A, Karlsen HC. Numerical calculation of the level ciossing rate
of second order stochastic Volterra systems. Probabilistic Engineering Mechanics 2004;19(2): 155-60.
[13] Naess A , Karlsen HC, Teigen PS. Numerical methods for calculating the crossing rate o f high and extreme response levels of compliant offshore structures subjected to random waves. Applied Ocean Research 2006; 28(1): 1-8.
230 A. Naess et al./Applied Ocean Research 29 (2007) 221-230
[15] Adegeest L, Braathen A , Vada T. Evaluation of methods for estimation of extreme nonlinear ship responses based on numerical simulations and model tests. In: Proc. 22nd symposium on naval hydrodynamics (1998), vol. 1. Washington: National Acadetny Press; 1999. p. 70-84.
[16] Naess A, Gaidai O, Teigen PS. Simulation based estimation of extreme response of floating offshore structures. In: Proceedings 26th international conference on offshore mechanics and arctic engineering. New York: ASME, 2007. p. OMAE-2007-29079.
[17] Teigen PS. Nonlinear wave loads on the shafts of a GBS. In: Proceedings of ISOPE'2004. 2004 ISOPE Paper IL-01.
[18] Rice SO. Mathematical analysis of random noise. In: Wax N , editor. Selected papers on noise and stochastic proces.ses. New York: Dover Publications, Inc.; 1954. p. 133-294,
[19] Naess A. On the long-term statistics of extreines. Applied Ocean Research 1984;6(4):227-8.
[20] Naess A, Moan T. Probabilistic design of offshore structures. In: Chakrabaiti SK. editor. Handbook of offshoie engineering, vol, I , Amsterdam: Elsevier; 2005, p. 197-277 [chapter 5].
[21] Walpole RE, Myers RH, Myers SL, Ye K, Probability and statistics for engineers and scientists. Upper Saddle River (NJ): Prentice Hall; 2002. [22] Naess A. Gaidai O, Monte Carlo methods for estimating the extreme
response of dynamical systems. Journal of Engineering Mechanics, ASCE. 2008 [in press],
[23] Leadbetter MR, Lindgren G. Rootzen H. Extremes and related properties of random sequences and processes. New York: Springer-Verlag; 1983. [24] Naess A, Gaidai O. Haver S, Efficient estimation of extreme response
of drag dominated offshoi e structures by Monte Cailo simulation. Ocean Engineeiing 2007;34( 16):2188-97,
