Stability analysis of differential scheme for dynamic equations of mooring cable system
Ruizhi Shi^{1} , Suxia Zhang^{2} , Xijun Liu^{3} , Fujiang Cui^{4} , Anqi Zhou^{5}
^{1, 2, 3, 4, 5}Department of Mechanics, School of Mechanical Engineering, Tianjin University, Tianjin, 300354, P. R. China
^{1, 2, 3, 4, 5}Tianjin Key Laboratory of Nonlinear Dynamics and Chaos Control, Tianjin, 300354, P. R. China
^{2}Corresponding author
Vibroengineering PROCEDIA, Vol. 20, 2018, p. 267272.
https://doi.org/10.21595/vp.2018.20282
Received 4 October 2018; accepted 9 October 2018; published 19 October 2018
JVE Conferences
The mooring cable system in plane motion can be modeled as two coupled partial differential equations, which can be numerical solved by finite difference method directly. The difference scheme is analyzed, and parameters selection for timemarching of displacement and velocity are deduced. The stability condition of the scheme is analyzed through Fourier series method, and parameters range which match stable scheme is given. Then, the parameters range is verified by a numerical example.
Keywords: mooring cable system, finite difference method, difference scheme, stability condition.
1. Introduction
Cable structure is widely used in engineering because of its light weight and good flexibility, and mooring cable system is one of them. It is necessary to analyze the dynamic characteristics and tension of the cable to ensure engineering safety.
In order to accurately analyze the characteristics of cable structure, a large number of scholars have studied a lot of cable models and calculation methods. In the 1980s, Irvine [1] systematically summarized the previous studies on the static and dynamic state of different cable structures. It had ignored nonlinear factors in previous studies. Later, when considering the great influence of nonlinear factors on the system motion, scholars gradually take the nonlinear factors in the study of dynamic cable models. Perkins et al. [24] derived the threedimensional dynamic equation of the fixed elastic cable through Hamilton principle, and studied the motion along the tangential, normal and auxiliary normal direction of the cable respectively; Then based on this model, the resonance response of elastic cable in steady flow is studied; Later, a theoretical model for the vibration of floating body with cable was established. According to this model, Zhang Suxia et al. [5, 6] studied the tension characteristic and impact tension under tautslack condition. This model for characterize the spatial dynamical motion of elastic cables is still used today.
When solving the coupled system of partial differential equations, the decoupling solution can be obtained by Galerkin method; Meanwhile, it also can be solved directly by finite element or finite difference method [7], and the difference scheme used in the finite difference method must be stable. Therefore, in this paper, the stability analysis of the difference scheme used to solve the dynamic equations of the mooring cable system is carried out, which in order to obtain the condition of the absolute stability of the scheme, and parameters of the time marching function are determined.
2. Model
Dynamic equations of mooring cable system in plane can be written as [2]:
where ${U}_{1}$ and ${U}_{2}$ represent tangential and normal displacement of the cable from static equilibrium respectively; $s$ is the arc length coordinate of the unstretched cable; $\rho $ is cable’s density, and ${\rho}_{w}$ is fluid's density; $A$ is the crosssectional area; $P$ is the pretension; $E$ is the elastic modulus; $\kappa $ is the curvature of the cable from static equilibrium; $\epsilon $ is the dynamic Lagrangian strain defined as Eq. (2); $g$ is the gravitational acceleration; ${l}_{t}=\mathrm{s}\mathrm{i}\mathrm{n}\theta $ and ${l}_{n}=\mathrm{c}\mathrm{o}\mathrm{s}\theta $ are direction cosines, and $\theta $ is the angle between tangential direction of the cable equilibrium and horizontal direction based on absolute coordinates; ${F}_{1}$ and ${F}_{2}$ are tangential and normal restoring forces of the cable in motion, which can be approximated by Morrison formula as Eq. (3):
where ${C}_{at}$, ${C}_{dt}$, ${C}_{an}$ and ${C}_{dn}$ are tangential and normal drag coefficients; $d$ is the effective diameter of cable.
3. Numerical solution by finite difference method
The model is a problem of partial differential equations, which can be numerical solved by finite difference method.
3.1. Difference scheme and coefficient determination
The spatial differentiation can be obtained in reference [8]. Fourthorder formulas for first and second derivatives are expressed as Eq. (4), (5):
where $k$ is the time coordinate and $j$ is the spatial coordinate.
A timemarching algorithm for solution of the nodal displacement ${u}_{j}^{k}$ and velocity ${v}_{j}^{k}$ is expressed in reference [9] as Eq. (6), (7):
The nodal acceleration ${a}_{j}^{k+1}$ at $k+$ 1 time point can be obtained by dynamic equations of the mooring cable system, while expressions on the right hand side of the equation are at k time point. $\mathrm{\Delta}t$ is the time step; ${\alpha}_{1}$, ${\alpha}_{2}$, ${\beta}_{1}$ and ${\beta}_{2}$ are weighting parameters, and the selection condition can be proposed by using the method of undetermined coefficients on secondorder backwards differencing in time as Eq. (8):
Here is in the same spatial coordinate, so subscript can be omitted.
Substitute Eq. (6), (7) into Eq. (8) leads to Eq. (9):
Due to $\frac{{\partial}^{2}{u}^{k+1}}{\partial {t}^{2}}={a}^{k+1}$, there are three conditions for selecting ${\alpha}_{1}$, ${\alpha}_{2}$, ${\beta}_{1}$ and ${\beta}_{2}$ as Eq. (10):
3.2. Stability analysis of difference scheme
By using Fourier series method (Von Neumann method) [7], the stability condition for the above scheme is analyzed.
The magnitude of $\mathrm{\Delta}s$ is associated with the length of the mooring cable, therefore nondimensional method is carried out for Eq. (1), let:
$\widehat{E}=\frac{E}{\rho {L}^{2}}\cdot 1\mathrm{}\mathrm{s}\mathrm{e}{\mathrm{c}}^{2},\widehat{\kappa}=\kappa \cdot L,\widehat{d}=\frac{d}{L},{\widehat{F}}_{1}=\frac{{F}_{1}}{M}\cdot 1\mathrm{}\mathrm{s}\mathrm{e}{\mathrm{c}}^{2},{\widehat{F}}_{2}=\frac{{F}_{2}}{M}\cdot 1\mathrm{}\mathrm{s}\mathrm{e}{\mathrm{c}}^{2},$
${{G}_{t}=\frac{\left(\rho {\rho}_{w}\right)Ag{l}_{t}}{M}\cdot 1\mathrm{}\mathrm{s}\mathrm{e}{\mathrm{c}}^{2},G}_{n}=\frac{\left(\rho {\rho}_{w}\right)Ag{l}_{n}}{M}\cdot 1\mathrm{}\mathrm{s}\mathrm{e}{\mathrm{c}}^{2},$
where $L$ is the original length of the mooring cable; and $M$ is the cable’s mass. For Fourier series method is effective only in the linear system, so linearize the equations is necessary. Through the cable static equilibrium, a relation equation between the pretension and the gravity can be obtained; then simplified coefficients of restoring forces are substituted by $\xi $:
Finally, a nondimensional linear difference scheme can be got:
Substitute Eq. (4), (5) into Eq. (11), then the discrete Fourier Transform is used as Eq. (12):
and derive out Eq. (13) to express $\widehat{u}$ which increases with time changes:
here, $G$ is a two dimensional growth matrix as:
where:
The eigenvalue equation is $\left(\lambda a\right)\left(\lambda d\right)bc=$ 0, also can be written as ${\lambda}^{2}\left(a+d\right)\lambda +adbc=$ 0. For the difference scheme as Eq. (11) is absolutely stable, the eigenvalue must satisfy the following condition:
This inequality can be separated two sets, one is:
and the other one is:
The difference scheme as Eq. (11) will be absolutely stable, if $\mathrm{\Delta}\widehat{t}$ and $\mathrm{\Delta}\widehat{s}$ satisfy either of Eq. (19), (20).
4. Numerical example
According to reference [10], structure parameters and environmental parameters of the mooring cable are taken as Table 1. The cable motion belongs to forced vibration, and there supposes a sine displacement excitation on the upper point of the cable.
Substitute parameters into Eqs. (14)(17) and make judgments about Eqs. (19), (20).
The results show that ${\left(ad\right)}^{2}+4bc\ge $ 0 is always tenable, so $\mathrm{\Delta}\widehat{t}$ and $\mathrm{\Delta}\widehat{s}$ which match condition Eq. (19) are empty set. $\mathrm{\Delta}\widehat{t}$ and $\mathrm{\Delta}\widehat{s}$ which match condition Eq. (20) are shown in Fig. 1. Values range of $\mathrm{\Delta}\widehat{t}$ and $\mathrm{\Delta}\widehat{s}$ are from 0 to 1. $\mathrm{\Delta}\widehat{s}$ labels horizontal axis of the figure, and $\mathrm{\Delta}\widehat{t}$ marks the vertical axis. The difference scheme will be absolutely stable, if $\mathrm{\Delta}\widehat{t}$ and $\mathrm{\Delta}\widehat{s}$ are in the red area. The three black lines along the border of red area are condition Eq. (20) when the two sides are equal. If $\mathrm{\Delta}\widehat{t}$ and $\mathrm{\Delta}\widehat{s}$ are out of the red area, the result may be divergent.
Table 1. Structure parameters and environmental parameters of the mooring cable
Parameter

Data

Parameter

Data

Original length $L$ (m)

1400

Fluid’s density ${\rho}_{w}$ (kg⋅m^{3})

1000

Effective diameter $d$ (m)

0.12

Tangential drag coefficients ${C}_{at}$

1

Linear density $\rho A$ (kg⋅m^{1})

70.37

Tangential drag coefficients ${C}_{dt}$

0.05

Crosssectional area $A$ (m^{2})

0.0113

Normal drag coefficients ${C}_{an}$

1

Elastic modulus $EA$ (kN)

136064

Normal drag coefficients ${C}_{dn}$

1.2

Horizontal pretension at upper point ${P}_{0}$ (N)

8.5985×10^{5}

Vertical pretension
at upper point $V$ (N)

1.29574×10^{6}

Fig. 1. $\mathrm{\Delta}\widehat{t}$ and $\mathrm{\Delta}\widehat{s}$ which match condition Eq. (20)
Extract two sets of data, respectively in and out of the red area, for numerical solution of dynamic equations of mooring cable system. The results are shown in Fig. 2.
The time histories of tangential and normal displacement at 700 m point are shown in Fig. 2(a), (b). In these results, $\mathrm{\Delta}\widehat{t}=$ 0.001 and $\mathrm{\Delta}\widehat{s}=$ 1/280 are selected form the stable area, and we can get stable numerical solution. When $\mathrm{\Delta}\widehat{t}=$ 0.01 and $\mathrm{\Delta}\widehat{s}=$ 1/280, the results are divergent shown as Fig. 2(c), (d). That means the value selection must from the stable area, then convergent results can be got.
Furthermore, the numerical solution is more accurate if $\mathrm{\Delta}\widehat{t}$ and $\mathrm{\Delta}\widehat{s}$ are smaller. This difference scheme is absolutely stable only in linear system or some weak nonlinear system, so it’s better to solve nonlinear stability without linearization if nonlinear terms have strong effect of results.
Fig. 2. Numerical solution of dynamic equations of mooring cable system at $\mathrm{\Delta}\widehat{s}=$ 1/280: a) tangential displacement at 700 m point $\mathrm{\Delta}\widehat{t}=$ 0.001, b) normal displacement at 700 m point $\mathrm{\Delta}\widehat{t}=$ 0.001, c) tangential displacement at 700 m point $\mathrm{\Delta}\widehat{t}=$ 0.01, d) normal displacement at 700 m point $\mathrm{\Delta}\widehat{t}=$ 0.01
a)
b)
c)
d)
5. Conclusions
The mooring cable system in plane motion is modeled as two coupled partial differential equations, which can be numerical solved by finite difference method without decoupling. This numerical method needs the difference scheme is stable and parameters selection is reasonable. Through the calculation and analysis, conclusions are obtained as follow:
1) Three conditions of timemarching parameters are derived, which used for displacement and velocity in the same space point, and second order accuracy can be satisfied to describe acceleration under these conditions;
2) The stability of the dimensionless linear difference scheme is analyzed by using Fourier series method, and conditions of the stability are obtained. Through a numerical example, the results show that the difference scheme can be stable under the linear system, and accurate numerical solutions can be obtained. However, in the nonlinear system, the influence of nonlinear terms should be considered.
Acknowledgements
The authors gratefully acknowledge the support of the National Natural Science Foundation of China (51479136), the Project of Tianjin Municipal Transportation Commission (2018B33), the Project of Tianjin Natural Science Foundation (17JCYBJC18700).
References
 Irvine Max H. Cable Structures. The MIT Press, Cambridge, Massachusetts, London, England, 1981. [Search CrossRef]
 Perkins N. C., Mote C. D. Threedimensional vibration of travelling elastic cables. Journal of Sound and Vibration, Vol. 114, Issue 2, 1987, p. 325340. [Publisher]
 Kim W. J., Perkins N. C. Twodimensional vortexinduced vibration of cable suspensions. Journal of Fluid and Structures, Vol. 16, Issue 2, 2002, p. 229245. [Publisher]
 Kim W. J., Perkins N. C. Linear vibration characteristics of cablebuoy systems. Journal of Sound and Vibration, Vol. 252, Issue 3, 2002, p. 443456. [Publisher]
 Tang Yougang, Zhang Suxia, Zhang Haiyan Study on snap tension induced by tautslack in a mooring system. Journal of Vibration and Shock, Vol. 27, Issue 4, 2008, p. 7072. [Search CrossRef]
 Zhang Suxia, Tang Yougang, Hou Shujun, Liu Xijun Characteristic analysis of stress wave in mooring line in tautslack condition. Engineering Mechanics, Vol. 27, Issue 5, 2010, p. 5359. [Publisher]
 Zhang Wensheng Finite Difference Methods for Partial Differential Equations in Science Computation. Higher Education Press, Beijing, China, 2006. [Search CrossRef]
 Schiesser W. E. The Numerical Method of Lines Integration of Partial Differential Equations. Academic Press, San Diego, 1991. [Search CrossRef]
 Nejad Mohammad Behbahani Wave Propagation in Elastic Cables with and Without Fluid Interaction. Mechanic Engineering in the University of Michigan, 1997, p. 79. [Search CrossRef]
 Wang Xinggang Coupled Dynamic Analysis of Deepwater Floating Structure and the Mooring system. Shanghai Jiao Tong University Press, Shanghai, China, 2014, p. 3637. [Search CrossRef]