# Simplified computational method for mooring forces of mooring systems

## Z. W. Wu1, J. K. Liu2, Z. Q. Liu3, Z. R. Lu4

1, 2, 3, 4Department of Mechanics, Sun Yat-sen University, Guangzhou 510006, P. R. China

4Corresponding author

Vibroengineering PROCEDIA, Vol. 5, 2015, p. 498-502.
Accepted 12 August 2015; published 18 September 2015

Copyright © 2015 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Views 166
Abstract.

The most common marine mooring systems of floating structures are taken as the numerical models. The effects of wave parameters and structural parameters on the wave forces are studied. A simplified formula of wave forces is established and the valuing ranges of key parameters are determined. Results obtained by the simplified computational method show great agreement with those obtained by numerical simulations and physical modelling experiments. The mooring forces can be obtained easily and efficiently by the proposed method.

Keywords: mooring forces, wave forces, numerical analysis, simplified method.

#### 1. Introduction

It is hard to determine the magnitude of mooring forces efficiently due to the lack of a simplified formula of wave forces. Frequently, the mooring forces are estimated by the empirical methods, which may be over-estimated or under-estimated. There are many numerical methods for determining the magnitude of mooring forces [1, 2]. However, they are generally expensive and time consuming for the use of professional software. Therefore, it is necessary and significant to develop a new type of simplified methods for calculating mooring forces of mooring systems.

The compositions of mooring forces are mainly current force, wind force and wave force. The Ref.  proposed simplified formulas of current force, wind force except wave force. Many researchers [4, 5] investigated extensively the effects of wave parameters and structural parameters on wave forces.

#### 2. Theory and models

Fig. 1 shows some numerical models of most common marine mooring floating structures. The equation of motion for moored floating body systems can be written as :

(1)
$\left(M+m\right)\stackrel{¨}{x}\left(t\right)+{\int }_{-\mathrm{\infty }}^{i}K\left(t-\tau \right)\stackrel{˙}{x}\left(t\right)d\tau +Cx\left(t\right)={F}_{WA}\left(t\right)+{F}_{M}\left(t\right)+{F}_{WI}+{F}_{CU},$

where $M$ and $m$ are the mass and added mass matrices, respectively, $K\left(t-\tau \right)$ is the delay function matrix, $C$ is the hydrostatic resilience matrix, $\stackrel{¨}{x}\left(t\right)$, $\stackrel{˙}{x}\left(t\right)$ and $x\left(t\right)$ are the acceleration, velocity and displacement matrices, respectively, and ${F}_{WA}\left(t\right)$, ${F}_{M}\left(t\right)$, ${F}_{WI}$ and ${F}_{CU}$ are the wave force, cable tension, wind force and current force, respectively.

The wind force and current force can be given by Ref. , and the active wave forces can be written as :

(2)

where $\rho$ is the density of water, ${F}_{j}$ is the active force in the $j$-th direction, ${n}_{j}$ is the generalized surface normal for the $j$th direction, $S$ is the wet surface of the body, ${F}_{fk}$ is the Froude-Krylov force, $\mathrm{\Phi }$ is the wave potential , and ${F}_{df}$ is the diffraction force.

In time domain, the wave forces can be calculated directly from the integral of the wave pressure acting on all the elements. The instantaneous wave force on each element is:

(3)

where $p$ is the dynamic wave pressure, $\stackrel{-}{N}$ is the outward normal vector of the body surface, and $S\left(t\right)$ is the wet surface of the body at an instant in time. In the case of the second order finite depth, $p$ can be obtained by the Pinkster approximation method.

Fig. 1. Numerical models a) Rectangular box b) Cylinder c) Immersed tunnel element d) Ship

#### 3. Simplified formula for calculating wave forces

The wave forces acting on the floating body can be decomposed into the lateral wave force ${F}_{x}$ and the longitudinal wave force ${F}_{y}$. All the factors affecting wave forces are taken into consideration to establish a simplified formula of wave forces. The considered factors of wave height, wave direction, wave period, structural length, and draft are expressed as influence factor ${k}_{i}$$\left(i=$ 1, 2, 3…$\right)$, while all the other unconsidered factors are expressed as synthesis factor $x$ and $y$, respectively.

#### 3.1.1. Wave height

In Fig. 2 (a), $y$-axis is the dimensionless lateral wave force ${F}_{x}$ divided by $\text{0.5}\rho g{a}^{2}L$, $a$ is the wave amplitude, $L$ is the structural length. It can be seen that the dimensionless lateral wave force is increased with the increase of wave height, which is basically in direct proportion to the wave height. Therefore the influence factor ${k}_{1}$ can be given as:

(4)
${k}_{1}=H,$

where $H$ is the incident wave height.

Fig. 2. The effect of wave direction and wave height on lateral wave force a) Wave height b) Wave direction

#### 3.1.2. Wave direction

In Fig. 2(b), one can see that the value of lateral wave force varies in sinusoidal curve with the increase of wave direction. Multiplying $\mathrm{s}\mathrm{i}\mathrm{n}\theta \left(1+\sqrt{\left|\mathrm{s}\mathrm{i}\mathrm{n}\theta \right|}\right)$ by each maximum value of numerical computing results to obtain fitted curves, it can be seen that the numerical curves coincide exactly with the fitted curves. Therefore the influence factor ${k}_{2}$ for symmetrical structure can be given as:

(5)

where $\theta$ is the incident wave direction.

#### 3.1.3. Wave period

According to the theoretical study and experimental verification [6, 7] it can easily cause parametric rolling when the rolling period is nearly double the wave period. In Fig. 3, it can be seen that with the increase of $T/{T}_{x}$, the dimensionless lateral wave force is increased, which reaches a maximum value at the point of 0.5, after that the its value attenuate. Making use of superposition principle to solve the differential equation by introducing frequency response function, the rolling amplification factor can be obtained. Thus, the influence factor ${k}_{3}$ can be given as:

(6)

where $T$ is the wave period, ${T}_{x}$ is the rolling period, ${\mu }_{x}$ is the rolling non-dimensional damping coefficient.

#### 3.1.4. Structural length and draft

It is known that the lateral wave force is directly proportional to structural length and draft , thus the influence factor ${k}_{4}$ can be given as:

(7)
${k}_{4}=LD,$

where $L$ is the length of structure, $D$ is the draft of structure.

By using Eqs. (4)-(7), and all the influence factors are multiplied by $\rho g$ to maintain dimensional consistency, ${F}_{x}$ can be obtained as:

(8)

where $x$ is the dimensionless lateral synthesis coefficient.

#### 3.2. Longitudinal wave force

Similarly, ${F}_{y}$ can be obtained as:

(9)

where $y$ is the dimensionless longitudinal synthesis coefficient, $B$ is the structural width, ${T}_{y}$ is the pitching period, and ${\mu }_{y}$ is the dimensionless pitching damping coefficient.

#### 3.3. Determination of key parameters of simplified formula

It is noted that  ${\mu }_{x}=\text{0.035}\text{-}\text{0.1}$ and ${\mu }_{y}=\text{0.3-0.5}$. A smaller value is chosen to make sure that the result is more safety, and then the fitting analysis for curves is carried to determine the value of $x$.

Fig. 3. Lateral colligation factor of simplified formula Fig. 4. Validation of the simplified method In Fig. 3, $y$-axis is the dimensionless lateral wave force ${F}_{x}$ which is divided by $H\mathrm{s}\mathrm{i}\mathrm{n}\theta \left(1+\sqrt{\left|\mathrm{s}\mathrm{i}\mathrm{n}\theta \right|}\right)LD\rho g$. ${F}_{x}$ is increased accordingly, which reaches a maximum value at the point of 0.5, then the value of force attenuates. It is clear that the variation of the fitted value curves agree well with the numerical value curves. When the value of $x$ varies from 0.1 to 0.7, the numerical value curves can be wrapped by fitted value curves basically. So the value of $x$ varies from 0.1 to 0.7. Similarly, the value of $y$ varies from 0.1 to 0.7. In most cases, the values of $x$ and $y$ are recommended as 0.5. In the conservative case, they are recommended as 0.7.

#### 4. Simplified method for calculating mooring forces

Then the standard value of mooring forces of mooring systems can be obtained as:

(10)
$N=\frac{K}{n}\left(\frac{\sum {F}_{x}}{\mathrm{s}\mathrm{i}\mathrm{n}\alpha \mathrm{c}\mathrm{o}\mathrm{s}\beta }+\frac{\sum {F}_{y}}{\mathrm{c}\mathrm{o}\mathrm{s}\alpha \mathrm{c}\mathrm{o}\mathrm{s}\beta }\right),$

where $N$ is the standard value of mooring forces, $K$ is the non-uniform coefficient of distribution of force, $\alpha$ is the angle between horizontal projection of mooring cable and $x$-axis, $\beta$ is the angle between mooring cable and horizontal plane, $\sum {F}_{x}$ and $\sum {F}_{y}$ are the total lateral force and the total longitudinal force induced by wind, current and wave, respectively, and $n$ is the number of loaded mooring cables.

#### 5. Validation of the simplified method

In Fig. 4, it can be shown that the results obtained by the proposed simplified method agree well with those obtained by numerical simulations and physical modelling experiments. In practice, the key parameters of formula can be chosen according to the specific circumstances of engineering application. Generally speaking, the conservative parameters are recommended.

#### 6. Conclusions

Parameters study is made for the effects of wave parameters and structural parameters on wave forces. A simplified method for calculation of mooring forces is proposed. The results from the simplified method agree well with those from numerical simulations and physical modelling experiments. The proposed simplified method can be used to obtain the mooring forces easily and efficiently, but it still needs to be further improved.

#### Acknowledgements

The study is supported by the National Natural Science Foundation of China (11172333, 11272361). Such financial aids are gratefully acknowledged.

1. Koo W. C., Kim M. H. Current effects on nonlinear wave-body interactions by a 2D fully nonlinear numerical wave tank. Journal of Waterway, Port, Coastal, and Ocean Engineering, Vol. 133, Issue 2, 2007, p. 136-146. [Search CrossRef]
2. Li B. A 3-D model based on Navier-Stokes equations for regular and irregular water wave propagation. Ocean Engineering, Vol. 35, 2008, p. 1842-1853. [Search CrossRef]
3. Prediction of Wind and Current Loads on VLCCs, 2nd Edition. OCIMF, Wither by and Co Ltd., London, 1994. [Search CrossRef]
4. Ming-Chung Fang, Gung-Rong Chen On the nonlinear hydrodynamic forces for a ship advancing in waves. Ocean Engineering, Vol. 33, 2006, p. 2119-2134. [Search CrossRef]
5. Dong G. H., Sun L., Zong Z., Zhao Y. P. Numerical analysis of the forces exerted on offshore structures by ship waves. Ocean Engineering, Vol. 36, 2009, p. 468-476. [Search CrossRef]
6. Buchner B., Dijk A., Wilde J. Numerical multiple-body simulation of side-by-side mooring to an FPSO. Proceedings of 11th ISOPE, 2001, p. 343-353. [Search CrossRef]
7. Sheng Zhen-bang, Liu Ying-Zhong Principles of Ship. Shanghai Jiao Tong University Press, Shanghai, 2009, (in Chinese). [Search CrossRef]