Continues Model for Vertical Vibration of Tension Leg ...

2 downloads 0 Views 1MB Size Report
Key-Words: - TLP, Axial vibration, Continues model, Wave load. 1 Introduction. Tension leg platform (TLP) is a well-known structure for oil exploitation in deep ...
Proceedings of the 9th WSEAS International Conference on Applied Mathematics, Istanbul, Turkey, May 27-29, 2006 (pp518-523)

Continues Model for Vertical Vibration of Tension Leg Platform M. R. TABESHPOUR, A. A. GOLAFSHANI, M. S. SEIF Department of Civil Engineering Sharif University of Technology Azadi Ave., Tehran, P.O. Box: 11365-9313 IRAN

Abstract: -In this paper the dynamic response of the leg of a tension leg platform (tether) subjected to the load simulated as ocean wave at the top of the leg is presented. The structural model is very simple but several complicated factors such as foundation effect, buoyancy and simulated ocean wave load are considered. Two continuous models are proposed to present the structural system and the mentioned effects. The problem is solved by means of non-harmonic Fourier expansion in terms of eigenfunctions obtained from a non-regular Sturm-Liouville system. Key-Words: - TLP, Axial vibration, Continues model, Wave load

1 Introduction Tension leg platform (TLP) is a well-known structure for oil exploitation in deep water. Many studies have been carried out to understand the structural behavior of TLP and to determine the effect of several parameters on dynamic response and average life time of the structure [1-4]. The most important point in the design of TLP is the pretension of the legs. The pretension causes that the platform behaves like a stiff structure with respect to the vertical degrees of freedom (heave, pitch and roll), whereas with respect to the horizontal degrees of freedom (surge, sway and yaw) it behaves as a floating structure. Among the various degrees of freedom, vertical motion (heave) is very important because of the direct effect on the stress fluctuation that leads to fatigue and fracture of tethers. Therefore the conceptual studies to understand the dynamic vertical response of TLP can be useful for designers. Rossit et al. (1996) presented an analytical solution for the dynamic response of the leg of TLP subjected to an axial suddenly applied load at one end [5]. The applied load was constant and the effect of the buoyancy was not considered. The aim of this paper is the solution of the mentioned

problem using two models. The structural models are very simple but several complicated factors such as foundation, buoyancy and simulated ocean wave loading are considered. At the first model the foundation assumed rigid but at the second model it is assumed that the foundation is embedded in the ocean bottom, which acts as a Winkler-type foundation. The buoyancy is modeled as a spring at the top of the leg. A concentrated force is applied at the top of the leg as simulated ocean wave load. The problem is solved by means of non-harmonic Fourier expansion in terms of eigenfunctions obtained from a non-regular Sturm-Liouville system [6, 7]. Tabeshpour et al., (2004a) have investigated the effect of added mass fluctuation on the heave response of tension leg platform for a discrete model by using perturbation method [8]. A continues model for vertical motion of TLP considering the effect of continues foundation has been reported by Tabeshpour et al., (2004b), [9]. Also Tabeshpour et al., (2004c) have presented a closed form formulation for the effect of added mass fluctuation on the heave response of tension leg platform considering continues model [10]. General configuration of TLP is shown in Fig. 1.

Proceedings of the 9th WSEAS International Conference on Applied Mathematics, Istanbul, Turkey, May 27-29, 2006 (pp518-523)

Fig. 2 - Dynamic structural model foundation material respectively, lt and l f are the length of tether and foundation respectively and δ denotes the Dirac delta function. The applied vertical load subjected to the mass m, is the N

generated wave load,

2 Analytical Solution of the Model The structural model of the system is shown in Fig. 2. The behavior of the system is described by the following differential equation

[ [

]

[ [

]

∑ F sin(Ω t + φ ) j

j

j

j =1

Fig. 1 - Configuration of TLP

⎛ u( y) − u( y − l f ) E f A f + ⎞ ∂ 2 v ⎜ ⎟ + ⎜ u ( y − l f ) − u ( y − l ) E t At ⎟ ∂y 2 ⎝ ⎠ Fh (t ) δ( y − l ) =

Fh (t ) =

]

obtained from the wave spectrum. The system is linear, therefore the solution of the equation (1) is carried out considering a single term input, Fh (t ) = F0 sin(Ωt ) , and then the overall response of the system is evaluated by summation of all responses. The initial conditions are v( y,0) = 0 ,

(1)

⎛ u( y) − u( y − l f ) ρ f A f + ⎞ ⎟ ∂ 2v ⎜ ⎜ u( y − l f ) − u( y − l ) ρt At + ⎟ 2 ⎟ ∂t ⎜ ⎜ M f δ( y − l f ) + M δ ( y − l ) ⎟ ⎠ ⎝

]

where u is step function, v is the axial deformation, E is the Young modulus of the tether material, At and A f are the cross sectional areas of the tether and foundation respectively, ρt and ρ f are the

∂v ( y ,0 ) = 0 ∂t

(2)

The mass distribution functions are defined as

[

]

m( y ) = u ( y ) − u ( y − l f ) ρ f A f +

[ u( y − l

f

]

) − u ( y − l ) ρ t A t + M f δ( y − l f ) +

(3)

M δ( y − l )

In the case of free vibration, Eq. (1) becomes

[ [

]

⎛ u( y) − u( y − l f ) E f A f + ⎞ ∂ 2 v ∂ 2v ⎟ ⎜ m ( y ) = ⎜ u ( y − l f ) − u ( y − l ) E t At ⎟ ∂y 2 ∂t 2 ⎠ ⎝

]

(4)

density of tether and Equation

[

(4)

can

]

be

solved

assuming

m( y ) = u ( y ) − u ( y − l f ) ρ f A f

[

]

+ u( y − l f ) − u( y − l ) ρ t At

subjected to the boundary conditions v(0, t ) = 0

(5)

Proceedings of the 9th WSEAS International Conference on Applied Mathematics, Istanbul, Turkey, May 27-29, 2006 (pp518-523)

∂v ∂ 2v (l f , t ) − M f (l f , t ) = ∂y ∂t 2 ∂v E f Af (l f , t ) ∂y

where α nt is the separation constant and and ct α nt = ω nt is the angular frequency.

− E t At

− k b v(l , t ) − M

∂ 2v ∂t

2

(l , t ) = Et At

(6)

It is clear that displacement and force are continues at y1 = l f or y 2 = 0 . From the continuity in

∂v (l , t ) ∂y

displacement one has

(7)

Yn1 (l f ) = Yn 2 (l t )

Equation (4) can be solved for two parts of the bar. In 0 ≤ y = y1 ≤ l f , one has E f Af

∂ 2v ∂y 2

= m( y )

∂ 2v ∂t 2

or

∂ 2v ∂y 2

=

1 ∂ 2v

or ⎡ ⎤ ⎛k ⎞ 1 M − α nt l t ⎟⎟ sin α nt l t ⎥ A⎢cos α nt l t + ⎜⎜ b ⎝ k t α nt l t m t ⎠ ⎣⎢ ⎦⎥ − B sin α nf l f = 0

(8)

c 2f ∂t 2

where m( y ) = ρ f A f , c 2f = E f ρ f . Using separation of variables, the eigenfunctions are determined as

Also the continuity in force gives

Yn1 = B sin α nf y1

E t At

where α nf

(9)

is the separation constant and and

∂v1 (l t , t ) = ∂y1

c f α nf = ω nf is the angular frequency.

− E f Af

In l f ≤ y = y 2 ≤ l , one has

or

Et At

∂ 2v ∂y 2

= m( y )

∂ 2v

or

∂t 2

where m( y ) = ρ t At ,

ct2

∂ 2v ∂y 2

=

1 ∂ 2v

ct2 ∂t 2

= Et ρ t .

Similarly using separation of eigenfunctions are determined as

variables,

⎛Mf 2 2 + B⎜ α l sin α nf l f − α nf l f cos α nf l f ⎜ m f nf f ⎝

the

(11)

⎞ ⎛k M 1 cos α nt lt + ⎜⎜ b − α nt lt ⎟⎟ sin α nt l t k l m α t ⎠ ⎝ t nt t ⎡ ⎤ ⎛k ⎞ kt M 1 α nt lt ⎢sin α nt lt − ⎜⎜ b − α nt lt ⎟⎟ cos α nt l t ⎥ kf ⎢⎣ ⎥⎦ ⎝ k t α nt lt mt ⎠

Now the frequency equation is resulted as

(14)

⎞ ⎟=0 ⎟ ⎠

The coefficients A and B can be determined by solving Eqs. (12) and (14). There are nonzero solutions if

Yn 2 = ⎡ ⎤ ⎛k ⎞ 1 M − α nt l t ⎟⎟ sin α nt y 2 ⎥ A⎢cos α nt y 2 + ⎜⎜ b ⎝ k t α nt l t m t ⎠ ⎣⎢ ⎦⎥

(13)

∂v 2 ∂ 2 v1 (l f , t ) − M f (l t , t ) ∂y 2 ∂t 2

⎡sin α nt l t − ⎤ kt ⎢ ⎥ α nt l t ⎢⎛ k b 1 A ⎞ M ⎜ ⎟ − α nt l t ⎟ cos α nt l t ⎥ kf ⎢⎜ k α l ⎥ ⎠ ⎣⎝ t nt t mt ⎦

(10)

(12)

− sin α nf l f Mf mf

=0 2 2 α nf lf

sin α nf l f − α nf l f cos α nf l f

(15)

Proceedings of the 9th WSEAS International Conference on Applied Mathematics, Istanbul, Turkey, May 27-29, 2006 (pp518-523)

⎞ ⎡⎛ k b 1 ⎤⎛ M f ⎞ M − α nt lt ⎟⎟ tan α nt l t +1⎥⎜ − α nf l f tan α nf l f + 1⎟ − ⎢⎜⎜ ⎜ ⎟ ⎠ ⎣⎢⎝ k t α nt l t mt ⎦⎥⎝ m f ⎠ ⎡ ⎛k ⎞⎤ k t α nt l t 1 M − α nt l t ⎟⎟⎥ = 0 tan α nf l f ⎢ tan α nt lt − ⎜⎜ b k f α nf l f ⎝ k t α nt lt mt ⎠⎦⎥ ⎣⎢

where:

ρ t At lt = mt : total mass of the tether;



ρ f A f l f = m f : total mass of the foundation;



Et At lt = k t : the axial stiffness of the tether, and E f A f l f = k f : axial stiffness of the foundation. The response of the tether subjected to axial load, can be expressed in terms of normal modes of the system ∞

v( y, t ) =

∑ Y ( y)T (t ) n

f

]

]

) − u ( y − l ) Y 2( y )

[

(18)

]

M ( y) = u ( y) − u( y − l f ) ρ f A f +

[ u( y − l

f

]

) − u ( y − l ) ρ t A t + M f δ( y − l f )

(19)

+ M δ( y − l )





0

l 0

0

lf 0

lt 0

M 2 ( y 2 )Yn 2 ( y 2 )Yr 2 ( y 2 ) = 0 (n ≠ r )

M ( y )Yn ( y )Yr (t ) = 0 (n ≠ r )

M ( y )Yn ( y )Yr (t ) = H r (n = r )

(20a) (20b)

M 2 ( y 2 )Yn 2 ( y 2 )Yr 2 ( y 2 ) = H r (n = r )

H r = H r1 + H r 2

(24a)

⎛Mf + Yr21 (l f )⎜⎜ 2 ⎝ 2

mf

⎞ ⎟ ⎟ ⎠

(24b) 2

⎞ mt mt ⎛k M 1 + + H r 2 = ⎜⎜ b − α nt l t ⎟⎟ k l m α 2 t ⎠ 2 ⎝ t nt t ⎛ kb ⎛Mf ⎞ M⎞ ⎟ + Yr22 (0)⎜ Yr22 (l t )⎜⎜ + ⎟⎟ 2 2 ⎟ ⎜ 2 ⎠ ⎝ 2 ⎠ ⎝ 2α rt c t

[

(24c)

]

Yn ( y )dy = u ( y ) − u ( y − l f ) Y1 ( y )dy +

[ u( y − l

f

]

) − u ( y − l ) Y 2( y )dy

and integrating between 0 and l , one obtains

[ [

]

⎛ u( y) − u( y − l f ) E f A f + ⎞ ⎜ ⎟× ⎜ u ( y − l f ) − u ( y − l ) E t At ⎟ ⎝ ⎠

]

∫ Y ( ∑ Y ′′T )dy + 0

r

n

(21)

Fh (t )

∫ Y δ( y)dy = ∫ r

0

(25)

n

l

l 0

M ( y )Yr

( ∑ Y T&&) dy n n

or

and M 2 ( y 2 ) = ρ t A t + M δ( y 2 ) + M f δ( y 2 − l t )

(23b)

l

Defining M 1 ( y1 ) = ρ f A f + M f δ( y1 − l f )

(23a)

M 1 ( y1 )Yn1 ( y1 )Yr1 ( y1 ) +

Multiplying Eq.(1) by

Because of the orthogonality of the normal modes, it can be shown that l

lt

M 1 ( y1 )Yn1 ( y1 )Yr1 ( y1 ) +

H r1 =

Yn ( y ) = u ( y ) − u ( y − l f ) Y1 ( y ) +

[ u( y − l



0

(17)

n

n =1

[



lf

(16)

(22)

∫ Y ( ∑ Y ′′ T )dy + E A ∫ Y ( ∑ Y ′′ T )dy + lf

E f Af

0

r1

n1

n

1

l

Equation (20) can be rewritten as

t

t

lf

r2

n2

n

2

(26)

& Fh (t ) Yr 2 (0) = H r T& r

since Yn1 satisfies (8) and Yn 2 satisfies (10), one has

Proceedings of the 9th WSEAS International Conference on Applied Mathematics, Istanbul, Turkey, May 27-29, 2006 (pp518-523)

M 1 ( y1 ) 2 2 c f α nf Yn1 E f Af

Yn′′1 = −

where

(27)

k n2 = c 2f α 2nf M (y ) Yn′′2 = − 2 2 ct2 α 2nt Yn 2 E t At

(28)

H n1 H + ct2 α 2nt n 2 Hn Hn

(31)

Initial conditions results in

Substituting (27) and (28) in (26) and applying (23) results in & + ⎛⎜ c 2 α 2 H n1 + c 2 α 2 H n 2 T& n t nt ⎜ f nf H Hn n ⎝ Yn 2 (0) F0 sin(Ωt ) Hn

⎞ ⎟Tn = ⎟ ⎠

An = 0 ; Bn = −

Ω Yn 2 (0) − Ω kn H n

F0 k n2

2

(32)

and (29)

Tn (t ) =

F0 Ω Yn 2 (0) × 2 k − Ω kn H n 2 n

⎛ Ω ⎞ ⎜ − sin k n t + sin(Ωt ) ⎟ ⎜ k ⎟ ⎝ n ⎠

Solving the above differential equation, one has

(33)

Tn (t ) = An cos k n t + B n sin k n t +

or

(30)

Yn 2 (0) sin(Ωt ) Hn k −Ω F0

2 n

2

⎞ Ω Yn 2 (0) ⎛ Ω ⎜− sin k n t + sin(Ωt ) ⎟⎟ ⎜ k H k −Ω n n ⎝ n ⎠

F0 k n2

Tn (t ) =

(

)

⎛ kb ⎞ ⎛ kb 1 ⎞ mt mt M M ⎞⎟ 2 ⎜ ⎟+⎜ ⎟ ⎟ ⎜ k α l − m α nt lt ⎟ 2 + 2 + Yr 2 (0)⎜ 2α 2 c 2 + 2 ⎟ t ⎠ ⎠ ⎝ t nt t ⎝ rt t ⎠

Fj

Ω j Yn 2 (0) ⎛ Ω j ⎞ ⎜− sin k n t + sin(Ω j t + φ j ) ⎟⎟ k n H n ⎜⎝ k n ⎠

⎛Mf + Yr21 (l f ) + Yr22 (lt ) ⎜⎜ 2 ⎝ 2

mf

2

2

(34)

and

k n2

Tnj (t ) =

mf 2

+

(

Yr21 (l f

− Ω 2j

) + Yr22 (lt )

)

⎛Mf ⎜ ⎜ 2 ⎝

2

⎛ kb ⎞ ⎛ kb 1 ⎞ mt mt M ⎞⎟ M 2 ⎜ ⎟+⎜ ⎟ ⎟ ⎜ k α l − m α nt lt ⎟ 2 + 2 + Yr 2 (0)⎜ 2α 2 c 2 + 2 ⎟ t ⎠ ⎠ ⎝ t nt t ⎝ rt t ⎠

(35)

and

[

3 Conclusion

]

Yn ( y ) = u ( y ) − u ( y − l f ) sin α nf y1 +

[ u( y − l

f

]

) − u( y − l ) ×

(36)

⎡ ⎤ ⎛ kb 1 ⎞ M − α nt l t ⎟⎟ sin α nt y 2 ⎥ ⎢cos α nt y 2 + ⎜⎜ ⎝ k t α nt l t m t ⎠ ⎣⎢ ⎦⎥

Now the dynamic response of the tether is ∞

v( y, t ) =

N

∑∑ Y ( y)T n

nj (t )

(37)

The analytical solution of the tether response of TLP was presented for a simple continuous model. The applied load is a simulation of ocean wave. Some complicated factors such as foundation effect and buoyancy were considered. The presented solutions give a conceptual view of the heave response of TLP under sea wave loads. The formulation presented herein can be used in analytical study on fatigue life of tethers.

n =1 j =1

The dynamic stress of the tether becomes σ( y, t ) = EA

∂v ( y, t ) ∂t

(38)

References

Proceedings of the 9th WSEAS International Conference on Applied Mathematics, Istanbul, Turkey, May 27-29, 2006 (pp518-523)

[1] Ahmad, S., 1996, Stochastic TLP response under long crested random sea, Journal of Computers and Structures, 61(6), pp. 975-993. [2] Jain, A.K., 1997, Nonlinear coupled response of offshore tension leg platforms to regular wave forces, Ocean Engineering, 24 (7), pp. 577-592. [3] Faltinsen, O.M., Van Hooff, R.W., Fylling, I.J., and Teigen, P.S., 1982, Theoretical and experimental investigations tension leg platform behavior, Proceedings of BOSS 1, pp. 411-423. [4] Lee, C.P., 1994. Dragged surge motion of a tension leg structure. International Journal of Ocean Engineering 21 (3), 311-328. [5] Rossit, C.A., Laura, P.A.A., and Bambill, D.V., 1996, Dynamic response of the leg of a TLP subjected to an axial suddenly applied load at one end, Ocean Engineering, 23(3), pp. 219-224. [6] Benedek, A., Guichal, E., and Panzone, R., 1974, On certain non-harmonic Fourier expansions as eigenfunctions of non-regular Sturm-Liouville systems, Instituto de

Matematica, Publication No. 4, Universidad Nacional del Sur, Bahia Blanca, Argentina. [7] Laura, P.A.A., Reyes, J.A., and Rossi, R.E., 1974, Analysis of a cable-like system suddenly stopped at one end, J. Sound Vibr. 37, pp. 195204. [8] Tabeshpour, MR, Golafshani, AA, and Seif, MS (2004a). The Effect of Added Mass Fluctuation on Heave Response of Compliant Offshore Structures, Proc 8th International Conference on Mechanical Engineering, ISME, Iran, (In Persian). [9] Tabeshpour, MR, Golafshani, AA, and Seif, MS (2004b). Simple Models for Heave Response of TLP under Harmonic Vertical Load, Proc 8th International Conference on Mechanical Engineering, ISME, Iran. [10] Tabeshpour, M.R., Seif, M.S. and Golafshani, A.A., (2004c), Vertical response of TLP with the effect of added mass fluctuation, 16th Symposium on Theory and Practice of Shipbuilding, Crovatia.