Estimation of Rolling Motion of Ship in Random Beam Seas by Efficient Analytical and Numerical Approaches

M. Salai Mathi Selvi L. Rajendran Marwan Abukhaled

Selvi M. Salai Mathi, Rajendran L., Abukhaled Marwan (2021). Estimation of Rolling Motion of Ship in Random Beam Seas by Efficient Analytical and Numerical Approaches. Journal of Marine Science and Application, 20(1): 55-66. https://doi.org/10.1007/s11804-020-00183-x
Citation: Selvi M. Salai Mathi, Rajendran L., Abukhaled Marwan (2021). Estimation of Rolling Motion of Ship in Random Beam Seas by Efficient Analytical and Numerical Approaches. Journal of Marine Science and Application, 20(1): 55-66. https://doi.org/10.1007/s11804-020-00183-x

Estimation of Rolling Motion of Ship in Random Beam Seas by Efficient Analytical and Numerical Approaches

https://doi.org/10.1007/s11804-020-00183-x
    Corresponding author:

    Marwan Abukhaled, mabukhaed@aus.edu

  • Abstract

    A steady-state roll motion of ships with nonlinear damping and restoring moments for all times is modeled by a second-order nonlinear differential equation. Analytical expressions for the roll angle, velocity, acceleration, and damping and restoring moments are derived using a modified approach of homotopy perturbation method (HPM). Also, the operational matrix of derivatives of ultraspherical wavelets is used to obtain a numerical solution of the governing equation. Illustrative examples are provided to examine the applicability and accuracy of the proposed methods when compared with a highly accurate numerical scheme.

     

    Article Highlights
    • The second-order nonlinear differential equation is solved using a modified approach of the homotopy perturbation method.
    • A mathematical model for rolling motion of ships in random beam seas is discussed.
    • Analytical expression of the roll angle, velocity, acceleration, and damping and restoring moments is derived.
    • Ultraspherical wavelet–based method was also employed to derive a numerical approximation of roll angle and restoring and damping moments.
  • Ships, in general, may experience three types of displacement motions (heave, sway or drift, and surge) and three angular motions (yaw, pitch, and roll) as depicted in Figure 1. Ship roll reduction or stabilization consists of two main parts: The first is the roll reduction or stabilization method, and the second is the evaluation or modeling method of the roll performance. Started with Froude (1861), roll motion and related topics have been studied since the late-mid 1800s. For example, Tanaka (1961), Kato (1965), and Ikeda et al. (2004) studied roll damping based on theory, numerical calculation, and experimentation and suggested an empirical formula including the effect of a bilge keel. A simple method for predicting the roll damping of a ship at forward speed was proposed by Himeno (1981) while Ikeda suggested a modified model (Ikeda et al. 2004).

    Figure  1  Schematic diagram of ship showing the six degrees of freedom
    Download: Full-Size Img

    The dynamic roll characteristic has been one of the most central topics in sea keeping study. However, due to the complexity induced by high nonlinearity, potential flow methods cannot solve the roll motion problems effectively. Up to date, many research groups have adopted various methodologies, including the advanced experimental process, empirical formulas, or computational fluid dynamics (CFD). One of the topics is related to the characteristics of roll motion for bare hull and with different appendages. Another branch of related research has focused on the analysis method for roll damping, which includes a higher polynomial roll damping model, finite differential model, or hyperbola model. Recently, CFD-based methods have been used for the analysis of local phenomena such as vortex shedding around bilge keels (Yeung et al. 1997).

    Several authors have studied roll damping experimentally for local flow visualization (see, for example, Aloisio and Felice 2006; Bassler et al. 2007; Oliveira and Fernandes 2012). On the other hand, mathematical modeling of roll damping has been investigated by researchers, for example, Oliveira and Fernandes (2014) used the bilinear or hyperbola fitting method, and Agarwal (2015) used a fractional differential equation model. However, modeling damping motion in terms of a second-order ordinary differential equation (quadratic model) is widely accepted. Recently, Oliveira and Fernandes (2012) discussed the nonlinear roll damping of FPSO hull. Although these methods have advantages, especially for newer hull designs or nontypical appendages in terms of shape and size, the traditional methods for a typical design have not yet been fully validated.

    Behavior analysis of a ship affected by external forces and analysis of ship stability in waves are based on more complex mathematical models as a result of strong nonlinear terms in the governing equations of these systems. Nonlinear analysis makes it possible to comprehend the theory of ship behavior for forecasting stability changes in various conditions of operation. In addition, ship turning motion is an important part of ship maneuverability and is directly related to the safety of a sailing (Lihua et al. 2018). The flow in the vicinity of 2D ship sections carrying out forced roll motions is simulated using the solution of the Navier-Stokes equations (Lavrov et al. 2017). A new coarse and fine-tuning fixed grid wavelet network is proposed for online predicting ship roll motion in regular waves (Huang et al. 2018). Kianejad et al. (2019) improved prediction of a ship roll motion in parametric roll and dead ship using CFD. Another important related topic that has been under investigation is fin stabilizer control. Demirel and Alarcin (2016) used LM-based H2 and H state-feedback controller for fin stabilizer of a fishing boat and also presented a backstepping control design procedure for nonlinear fin roll control of a trawler (Demirel et al. 2017).

    The second-order nonlinear differential equation modeling the ship roll motion does not have a closed-form solution. Despite the fact that reliable numerical and experimental methods have been implemented to find approximate solutions, the need for analytical solutions is still a necessity in order to understand the effect of parameters' variation on the governing system. Of the well-established methods that can be employed to obtain analytical expressions of roll angle, velocity, and damping and restoring moments, we mention the variational iteration method (Lu 2007; Abukhaled 2013), differential transformation method (Chen and Liu 1998), Green function based method (Khuri and Abukhaled 2017; Abukhaled 2017), finite element method (Comini et al. 1970), series method (Salomi et al. 2020), homotopy perturbation method (He 1999), and homotopy analysis method (Liao and Chwang 1998; Liao 1997; Liao 2004). Another approach that can be employed is that the differential equation may be analyzed based on the formulation of an equivalent integral equation (Jang et al. 2009; Jang 2011; Jang 2013).

    In this project, the governing equation will be solved analytically by using a modified form of the homotopy perturbation method and numerically by using a high-order ultraspherical wavelets–based method. Wavelet-based methods using operational matrix of derivatives are accurate and reliable and have been intensively used in many fields of science and engineering applications. For example, wavelet-based methods are used to find approximate solutions for nonlinear differential equations (Salai Mathi Selvi et al. 2017), for Michaelis-Menten enzymatic reaction equation Salai Mathi Selvi and Hariharan (2017), and for solving integro-differential equations (Tavassoli Kajani and Hadi Vencheh 2004). Wavelet-based methods have also been employed to solve optimal control and fractional optimal control problems (Razzaghi and Yousef 2002; Sadek et al. 2007; Abualrub et al. 2018). Recently, Salai Mathi Selvi and Hariharan (2016) applied the Chebyshev wavelet algorithm for solving steady-state concentration of packed-bed reactor model, Zheng and Wei (2016), used Legendre wavelets to estimates the error of approximation when solving integral equations, and Abualrub and Abukhaled (2015) employed wavelets to regulate cellular processes for the growth or regeneration of a tissue within an assigned terminal time.

    This paper is organized as follows: Section 2 describes the mathematical formulation of the problem. Section 3 discusses the properties of ultraspherical wavelets and a wavelet-based solution method. In Section 4, the homotopy perturbation method for obtaining an analytical solution is explained. Illustrative examples are presented in Section 5 with their discussions in Section 6. Finally, Section 7 presents our conclusions.

    Under the assumption of independent oscillation, the nonlinear rolling motion of a ship is described by the differential equation (Cardo et al. 1981)

    $$ I\ddot{\phi}+D\left(\varphi, \dot{\varphi}\right)+{M}_r\left(\varphi, t\right)={E}_W\cos \varOmega t $$ (1)

    where ϕ is the roll angle to the calm sea surface and $ \dot{\phi} $ and $ \ddot{\phi} $ are the roll velocity and acceleration, respectively. Equation (1) represents the general equation of roll motion of ships in the absolute heeling angle, where the motion is induced by regular transversal waves. I is the mass moment of inertia including added mass moment (kg · m2), $ D\left(\phi, \dot{\phi}\right) $ is the moment of the dissipative forces (N · m), Mr(φ, t) is the righting moment (N · m), and Ew represents the amplitude and Ωis the angular frequency. In this study, the nonlinear damping term is assumed to take into consideration an angular dependence of the linear term, that is,

    $$ D\left(\phi, \dot{\phi}\right)=\left({D}_{01}+{D}_{21}{\phi}^2\right)\dot{\phi}+{D}_{03}{\dot{\phi}}^3 $$ (2)

    and the righting moment is expressed in the normalized form

    $$ {M}_r^{\ast}\left(\phi \right)={\omega}_0^2\phi +{\alpha}_3{\phi}^3 $$ (3)

    Introducing the angle and time scales ϕn and Tn, the following dimensionless parameters can be defined

    $$ x=\frac{\phi }{\phi_n},\tau =\frac{t}{T_n},{\omega}_0^2=\frac{\varDelta {k}_1{T}_n^2}{I},\omega =\varOmega {T}_n,{\alpha}_i\\ =\varDelta {k}_i{T}_n^2{\phi}_n^{i-1},\mu =\frac{D_{01}{T}_n}{2I},{\delta}_1=\frac{D_{21}{T}_n{\phi}_n^2}{I},{\delta}_2\\ =\frac{D_{03}{\phi}_n^2}{I{T}_n},\varepsilon =\frac{E_w{T}_n^2}{I{\phi}_n} $$ (4)

    then, using these parameters and rearranging terms gives the following cubic damping moment of nonlinear roll motion (Cardo et al. 1981)

    $$ \ddot{x}+\left(2 \mu+\delta_{1} x^{2}\right) \dot{x}+\delta_{2} \dot{x}^{3}+\omega_{0}^{2} x+\alpha_{3} x^{3}=\varepsilon \cos (\omega t), $$ (5)

    where $ {\omega}_0^2 $ is the dimensionless angular frequency, ω is the frequency, α is strength of nonlinearity coefficient, μ is dimensionless damping coefficient, δ1 and δ2 are quadratic viscous damping coefficients, and ε is the amplitude.

    In this section, we use the function presentation via shifted ultraspherical wavelets to construct a wavelet-based method to obtain a numerical solution for the governing Eq. (5). Appendices 1 and 2 provide simple prelude to this section (see Doha et al. 2016 and the references therein for more in-depth details).

    Ultraspherical wavelets $ {\psi}_{nm}^{\left(\alpha \right)}\left(\xi \right)=\left(k,n,m,\alpha, \xi \right) $ are defined on the interval [0, 1] by Doha et al. (2016).

    $$ {\psi}_{nm}^{\left(\alpha \right)}\left(\xi \right)=\left\{\begin{array}{c}\frac{2^{k+\frac{1}{2}}}{\sqrt{B_m}}{C}_m^{\left(\alpha \right)}\left({2}^{k+1}\xi -\left(2n+1\right)\right),t\in \left[\frac{n}{2^k},\frac{n+1}{2^k}\right]\\ {}0,\kern2.5em \mathrm{otherwise}\end{array}\right., $$ (6)

    where m = 0, 1, 2, …M, n = 0, 1, 2, …2k − 1, m is the constructed polynomial order, α is a prescribed parameter, ξ the normalized time, and

    $$ {B}_m=\frac{{\left(2\alpha \right)}_m\left(m+\alpha \right)\varGamma \left(\alpha \right)}{\sqrt{\pi }m!\varGamma \left(\alpha +\frac{1}{2}\right)}. $$ (7)

    A function f(t) defined on [0, 1] may be expanded in terms of ultraspherical wavelets as

    $$ f(t)={\sum}_{n=0}^{\infty }{\sum}_{m=0}^{\infty }{c}_{nm}{\psi}_{nm}^{\left(\alpha \right)}(t), $$ (8)

    where

    $$ {c}_{nm}={\left(f(t),{\psi}_{nm}^{\left(\alpha \right)}(t)\right)}_w={\int}_0^1{\left(t-{t}^2\right)}^{\alpha -\frac{1}{2}}f(t){\psi}_{nm}^{\left(\alpha \right)}(t) dt. $$ (9)

    Using the truncated finite series, f(t) can be approximated in terms of ultraspherical wavelets as

    $$ f(t)={\sum}_{n=0}^{2^k-1}{\sum}_{m=0}^{M-1}{c}_{nm}{\psi}_{nm}^{\left(\alpha \right)}(t)={C}^T{\psi}^{\left(\alpha \right)}(t), $$ (10)

    where

    $$ C={\left[{c}_{0,0},{c}_{0,1},.\dots {c}_{0,M},{c}_{1,M},\dots {c}_{2^k-1,0},..\dots {c}_{2^k-1,M}\right]}^{\mathrm{T}}, $$ (11)
    $$ {\psi}^{\left(\alpha \right)}(t)={\left[{\psi}_{0,0}^{\left(\alpha \right)},{\psi}_{0,1}^{\left(\alpha \right)},.\dots {\psi}_{0,M}^{\left(\alpha \right)},{\psi}_{1,0}^{\left(\alpha \right)},.\dots {\psi}_{1,M}^{\left(\alpha \right)},\dots {\psi}_{2^k-1,0}^{\left(\alpha \right)},..\dots {\psi}_{2^k-1,M}^{\left(\alpha \right)}\right]}^{\mathrm{T}}. $$ (12)

    The operational matrix of derivatives is given by the following theorem.

    Theorem 1 For the ultraspherical wavelets vector (12),

    $$ \frac{d{\psi}^{\left(\alpha \right)}}{dt}(t)=\boldsymbol{D}{\psi}^{\left(\alpha \right)}(t), $$ (13)

    where D is a 2k(M + 1) × 2k(M + 1) operational matrix given by

    $$ \boldsymbol{D}=\left(\begin{array}{cccc}\boldsymbol{F}& 0& \cdots & 0\\ {}0& \boldsymbol{F}& \cdots & 0\\ {}\vdots & \vdots & \ddots & \vdots \\ {}0& 0& \cdots & 0\end{array}\right), $$ (14)

    in which F is an M + 1 square matrix whose (r, s)th entry is given by

    $$ {\boldsymbol{F}}_{r,s}=\left\{\begin{array}{cc}\frac{2^{k+1}}{\sqrt{B_{r-1}{B}_{s-1}}},& r\ge 2,r \gt s\ \mathrm{and}\ \left(r+s\right)\ \mathrm{odd}\\ {}0,& \mathrm{otherwise}\end{array}\right. $$ (15)

    Proof See Doha et al. (2016).

    The following theorem states that the ultraspherical wavelet expansion of a function f(t), with a bounded second derivative, converges uniformly tof(t).

    Theorem 2 A function $ f(t)\in {L}_{\overset{\sim }{\omega}}^2\left[0,1\right],\overset{\sim }{\omega }={\left(t-{t}^2\right)}^{\alpha -\frac{1}{2}},0 \lt \alpha \lt 1 $ can be expanded as an infinite series of ultraspherical wavelets, which converges uniformly to f(t), provided that f " (t) ≤ L.

    Moreover,

    $$ \left|{c}_{nm}\right| \lt \frac{4L{\left(1+\alpha \right)}^2\varGamma \left(m+1+\frac{3\alpha }{2}\right)\sqrt{m!\left(m+\alpha \right)}}{\varGamma \left(m+\frac{\alpha }{2}\right)\sqrt{\left(m+2\alpha \right)}{\left(m-2\right)}^4{\left(n+1\right)}^{\frac{5}{2}}},\forall n\ge 0,m \gt 2 $$ (16)

    Proof See Doha et al. (2016).

    The wavelet-based method for finding a numerical solution for a homogenous or nonhomogenous nonlinear equation describing the roll motion of a ship begins by using ultraspherical wavelets expansion along with a spectral collocation method to reduce the governing initial value problem into a system of algebraic equations with unknown coefficients.

    For finding an approximate numerical solution for Eq. (5), we begin by expanding x(t) and the function g(t) = cos ωt in terms of ultraspherical wavelets as follows

    $$ x(t)\approx {\sum}_{n=0}^{2^k-1}{\sum}_{m=0}^M{c}_{nm}{\psi}_{nm}^{\left(\alpha \right)}(t)={C}^{\mathrm{T}}\varPsi (t), $$ (17)
    $$ g(t)\approx \sum \limits_{n=0}^{2^k-1}\sum \limits_{m=0}^M{g}_{nm}{\psi}_{nm}^{\left(\alpha \right)}(t)={G}^{\mathrm{T}}\Psi (t), $$ (18)

    Using the operational matrix of derivatives gives

    $$ \dot{x}(t)\approx {C}^{\mathrm{T}}\boldsymbol{D}{\varPsi}^{\left(\alpha \right)}(t),\ddot{x}(t)\approx {C}^{\mathrm{T}}{\boldsymbol{D}}^2{\varPsi}^{\left(\alpha \right)}(t). $$ (19)

    Now substituting Eqs. (17)–(19) into Eq. (5) gives the residual written explicitly as

    $$ R(t)={C}^{\mathrm{T}}{\boldsymbol{D}}^2{\varPsi}^{\left(\alpha \right)}(t)\\ \;\;\;\;\;\;\;\;\;\;+\left(2\mu +{\delta}_1{\left({C}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(t)\right)}^2\right){C}^{\mathrm{T}}\boldsymbol{D}{\varPsi}^{\left(\alpha \right)}(t)\\ \;\;\;\;\;\;\;\;\;\;+{\delta}_2{\left({C}^{\mathrm{T}}\boldsymbol{D}{\varPsi}^{\left(\alpha \right)}(t)\right)}^3+{\omega}_0^2{C}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(t)\\ \;\;\;\;\;\;\;\;\;\;+{\alpha}_3{\left({C}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(t)\right)}^3-\varepsilon {G}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(t). $$ (20)

    A brief introduction to the basic concept of the homotopy perturbation method is delineated in Appendix 3.

    As per Eq. (4) (Appendix 3), the homotopy of Eq. (5) takes the form.

    $$ \left(1-p\right)\left[\begin{array}{c}\ddot{x}(t)+\left[2\mu +{\delta}_1{x}^2(t)\right]\dot{x}(t)+{\omega}_0^2x\left(t=0\right)\\ {}+{\alpha}_3{x}^3\left(t=0\right)-\varepsilon \cos \left(\omega t\right)\end{array}\right]\\ \;\;\;\;\;\;\;\;\;\;+p\left[\begin{array}{c}\ddot{x}(t)+\left[2\mu +{\delta}_1{x}^2(t)\right]\dot{x}(t)+{\delta}_2{\dot{x}}^3\left(t=0\right)\\ {}+{\omega}_0^2x(t)+{\alpha}_3{x}^3(t)-\varepsilon \cos \left(\omega t\right)\end{array}\right], $$ (21)

    where p ∈ [0, 1] is an embedding parameter. The approximate solution of Eq. (35), expressed in a power series, is

    $$ x={x}_0+p{\mathrm{x}}_1+{p}^2{x}_2+\dots $$ (22)

    By letting p = 1, an analytical solution, in the form of series expansion, is derived.

    In this section, we apply the analytical HPM method and the numerical ultraspherical wavelet method to cubic damping moment of nonlinear roll motion.

    Example 1 Consider the nonlinear BVP (Cardo et al. 1984)

    $$ \ddot{x}+\left(2\mu +{\delta}_1{x}^2\right)\dot{x}+{\omega}_0^2x+{\alpha}_3{x}^3=0, $$ (23)
    $$ {x}_0(0)=0.2,{\dot{x}}_0^{\prime }(0)=0, $$ (24)

    where the experimental values of the parameters are $ \mu =0.005,{\delta}_1=0.1,{\omega}_0^2=1,\mathrm{and}\ {\alpha}_3=-1.75. $

    Using the ultraspherical wavelet method, we substitute Eq. (17) and Eq. (19) into Eq. (23) to obtain

    $$ {C}^{\mathrm{T}}{\boldsymbol{D}}^2{\varPsi}^{\left(\alpha \right)}(t)+0.014\left({C}^{\mathrm{T}}\boldsymbol{D}{\varPsi}^{\left(\alpha \right)}(t)\right)+{\omega}_0^2\left({C}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(t)\right)-0.014=0. $$ (25)

    Also using Eq. (17) and Eq. (19), the initial conditions given in Eq. (38) become

    $$ {C}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(0)=0.2,{C}^{\mathrm{T}}\boldsymbol{D}{\varPsi}^{\left(\alpha \right)}(0)=0, $$ (26)

    where the operational matrix of derivatives is

    $$ \boldsymbol{D}=\left(\begin{array}{ccc}0& 0& 0\\ {}{F}_{2,1}& 0& 0\\ {}0& {F}_{3,2}& 0\end{array}\right), $$ (27)

    in which

    $$ {F}_{2,1}=\frac{\Gamma \left(1+\alpha \right)\sqrt{2\left(1+\alpha \right)}}{\sqrt{\pi}\Gamma \left(\frac{1}{2}+\alpha \right)},{F}_{3,2}=\frac{\Gamma \left(1+\alpha \right)\sqrt{2\left(1+\alpha \right)}\left(2+\alpha \right)\left(1+2\alpha \right)}{\sqrt{\pi}\Gamma \left(\frac{1}{2}+\alpha \right)}. $$ (28)

    Now, the approximate wavelet solution x(t) is

    $$ x(t)={C}^T{\varPsi}^{\left(\alpha \right)}(t), $$ (29)

    where the vector ψ(α)(t) is given by

    $$ {\psi}^{\left(\alpha \right)}(t)=\left(\begin{array}{c}1\\ {}\sqrt{2{B}_2}\left(2t-1\right)\\ {}\sqrt{2{B}_3}\frac{2\left(1+\alpha \right){\left(2t-1\right)}^2-1}{3}\end{array}\right), $$ (30)

    and the constants vector C is given by

    $$ \boldsymbol{C}=\left[\begin{array}{c}{c}_0\\ {}{c}_1\\ {}{c}_2\end{array}\right]. $$ (31)

    To find the HPM analytical solution for Eq. (23), we use the homotopy in Eq. (21) to find that the coefficient associated with p0 is the nonlinear equation

    $$ \ddot{x}+(0.014){\dot{x}}_0+{x}_0+0.0001112-0.2\cos (0.33333t)=0 $$ (32)

    Now x0 is obtained by solving Eq. (32) subject to initial conditions

    $$ {x}_0(0)=0.2,{\dot{x}}_0^{\prime }(0)=0, $$ (33)

    In a similar fashion, the nonlinear differential equation associated with p1 can be obtained and then solved for x1.

    Using only two terms of Eq. (22) with p = 1, we obtain the following approximate analytical solution of Eq. (23)

    $$ x(t)=0.014+{exp}\left(-0.007t\right)\left(0.186005{sin}\left(0.9999t+1.56381\right)\right). $$ (34)

    Example 2 Consider the cubic damping moment of nonlinear roll motion represented by the IVP (Cardo et al. 1984)

    $$ \ddot{x}+\left(2\mu +{\delta}_1{x}^2\right)\dot{x}+{\omega}_0^2x+{\alpha}_3{x}^3=\varepsilon \cos \left(\omega t\right) $$ (35)
    $$ {x}_0(0)=0.2,{\dot{x}}_0^{\prime }\ (0)=0, $$ (36)

    We will consider three sets of experimental values for the parameters.

    Case 1 Assume that

    $$ \mu =0.005,{\delta}_1=0.1,{\omega}_0^2=1,{\alpha}_3=-1.75,\varepsilon =0.2,\omega =0.33333. $$

    Using the ultraspherical wavelet method, we substitute Eqs. (17) and (19) into Eqs. (42) and (43) to obtain

    $$ {C}^{\mathrm{T}}{\boldsymbol{D}}^2{\varPsi}^{\left(\alpha \right)}(t)+0.014\left({C}^{\mathrm{T}}\boldsymbol{D}{\varPsi}^{\left(\alpha \right)}(t)\right)+{\omega}_0^2\left({C}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(t)\right)-0.014=0.2\cos \omega t, $$ (37)

    and

    $$ {\boldsymbol{C}}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(0)=0.2,{\boldsymbol{C}}^{\mathrm{T}}\boldsymbol{D}{\varPsi}^{\left(\alpha \right)}(0)=0. $$ (38)

    As detailed in the example above, the approximate numerical wavelet solution, x(t), is expressed in the form

    $$ x(t)={\boldsymbol{C}}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(t). $$ (39)

    Applying the proposed approach of the HPM to Eq. (35) gives the following analytical expression for the roll angle at any time t:

    $$ x(t)=\exp \left(-0.007t\right)\left[-0.03899\ \sin \left(0.9999t+1.55541\right)\right]\\ \;\;\;\;\;\;\;\;\;\;+\sin (0.9999t)\left[0.14999\ \sin \left(0.6666t+0.010467\right)\right]\\ \;\;\;\;\;\;\;\;\;\;+0.075\ \sin \left(1.3333t+0.0052\right)\left]\\ \;\;\;\;\;\;\;\;\;\;\;+\cos (0.9999t)\left[0.14999\ \sin \left(0.6666t+1.58126\right)\right]\\ \;\;\;\;\;\;\;\;\;\;+0.075\ \sin \left(1.3333t+1.576\right)\right]+0.014. $$ (40)

    The velocity and acceleration are, respectively, given by

    $$ \dot{x}(t)=\exp \left(-0.007t\right)\left[-0.03899\ \sin \left(0.9999t+3.13159\right)\right]-\sin (0.9999t)\\ \;\;\;\left[0.04999\ \sin \left(0.6666t-1.5812\right)\right]+0.025\ \sin \left(1.3333t+1.576\right)\left]\\ \;\;\;+\cos (0.9999t)\left[0.04999\ \sin \left(0.6666t+0.0104\right)\right]-0.025\ \sin \left(1.3333t+0.0052\right)\right], $$ (41)
    $$ \ddot{x}(t)=\exp \left(-0.007t\right)\left[-0.03899\ \sin \left(0.9999t+1.56772\right)\right]\\ \;\;+\sin (0.9999t)\left[-0.01666\ \sin \left(0.6666t-0.010204\right)\right]-0.00833\ \sin \left(1.3333t+0.0048\right)\left]\\ \;\;\;\;+\cos (0.9999t)\left[-0.01666\ \sin \left(0.6666t+1.581\right)\right]-0.00833\ \sin \left(1.3333t+1.5756\right)\right] $$ (42)

    Rearranging terms gives the cubic restoring moment formula

    $$ r(x)=x(t)-1.75x{(t)}^3, $$ (43)

    and the cubic damping moment formula

    $$ d\left(\dot{x}\right)=\left(0.01+0.1{x}^2\right)\dot{x}+0.1{\dot{x}}^3. $$ (44)

    Case 2 Assume the following experimental values of the parameters: $ \mu =0.02,{\delta}_1=0.1,{\omega}_0^2=1,{\alpha}_3=-1.75,\varepsilon =0.2,\mathrm{and}\ \omega =0.33333. $

    The numerical wavelet solution is given by Eq. (39) where Ψ(α)(t) is obtained by solving the equation

    $$ {\boldsymbol{C}}^{\mathrm{T}}{\boldsymbol{D}}^2{\varPsi}^{\left(\alpha \right)}(t)+0.004\left({\boldsymbol{C}}^{\mathrm{T}}\boldsymbol{D}{\varPsi}^{\left(\alpha \right)}(t)\right)+\left({\boldsymbol{C}}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(t)\right)-0.014=0.2\cos \omega t, $$ (45)

    subject to the initial conditions given in Eq. (38). The HPM analytical solution is

    $$ x(t)=\exp \left(-0.022t\right)\left[0.03899\ \sin \left(0.99975t+1.51707\right)\right]\\ \;\;\;\;\;\;\;\;\;\;+\sin (0.9997t)\left[0.150008\ \sin \left(0.6664t+0.03299\right)\right]\\ \;\;\;\;\;\;\;\;\;\;+0.075\ \sin \left(1.3331t+0.0165\right)\left]\\ \;\;\;\;\;\;\;\;\;\;\;+\cos (0.9997t)\left[0.07502\ \sin \left(1.3331t+1.58733\right)\right]\\ \;\;\;\;\;\;\;\;\;\;+0.150009\ \sin \left(0.6664t+1.6038\right)\right]+0.014 $$ (46)

    Case 3 Assume the following experimental values of the parameters: $ \mu =0.02,{\delta}_1=0.1,{\omega}_0^2=1,{\alpha}_3=4,\varepsilon =0.2,\mathrm{and}\ \omega =0.33333. $

    The numerical wavelet solution is given by Eq. (39) where Ψ(α)(t) is obtained by solving the equation

    $$ {C}^{\mathrm{T}}{D}^2{\varPsi}^{\left(\alpha \right)}(t)+0.044\left({C}^{\mathrm{T}}D{\varPsi}^{\left(\alpha \right)}(t)\right)+\left({C}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(t)\right)\\ \;\;\;+0.032\\ \;\;\;=0.2\cos \omega t, $$ (47)

    subject to the initial conditions given in Eq. (38).

    The analytical solution obtained by the HPM is given by

    $$ x(t)=\sin (0.9997t)\left[0.150008\ \sin \left(0.6664t+0.03299\right)\right]\\ \;\;\;\;\;\;\;\;\;\;+0.075\sin \left(1.3331t+0.0165\right)\left]\\ \;\;\;\;\;\;\;\;\;\;\;+\cos (0.9997t)\left[0.07502\ \sin \left(1.3331t+1.58733\right)\right]\\ \;\;\;\;\;\;\;\;\;\;+0.150009\ \sin \left(0.6664t+1.6038\right)\right]-0.032\\ \;\;\;\;\;\;\;\;\;\;+\exp \left(-0.022t\right)\left[0.03899\ \sin \left(0.99975t+1.51707\right)\right]. $$ (48)

    Example 3 Consider the cubic damping moment of nonlinear roll motion represented by the IVP (Liao and Chwang 1998)

    $$ \ddot{x}\left(1+4{b}^2{x}^2\right)+\eta x+4{b}^2{\dot{x}}^2x=0,\eta \gt 0, $$ (49)
    $$ {x}_0(0)=l,{\dot{x}}_0^{\prime }(0)=0, $$ (50)

    where the experimental values of the parameters are η = 0.1, b = 0.5, and l = 0.2.

    Using the ultraspherical wavelet method, we obtain the approximate solution

    $$ x(t)={C}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(t). $$ (51)

    where Ψ(α)(t) is obtained from the equation

    $$ {\boldsymbol{C}}^{\mathrm{T}}{\boldsymbol{D}}^2{\varPsi}^{\left(\alpha \right)}(t)\left(1+4{b}^2{\left({\boldsymbol{C}}^T{\varPsi}^{\left(\alpha \right)}(t)\right)}^2\right)+\\ \;\;\;\;\eta \left({\boldsymbol{C}}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(t)\right)+4{b}^2{\left({\boldsymbol{C}}^{\mathrm{T}}\boldsymbol{D}{\varPsi}^{\left(\alpha \right)}(t)\right)}^2\left({\boldsymbol{C}}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(t)\right)\\ \;\;\;\;=0, $$ (52)

    subject to the initial conditions take the from

    $$ {\boldsymbol{C}}^{\mathrm{T}}{\varPsi}^{\left(\alpha \right)}(0)=l,{\boldsymbol{C}}^{\mathrm{T}}\boldsymbol{D}{\varPsi}^{\left(\alpha \right)}(0)=0. $$ (53)

    In order to apply the homotopy perturbation method, we first construct the following homotopy equation

    $$ \ddot{x}+\eta x+p\left(4{b}^2\ddot{x}{x}^2+4{b}^2{\dot{x}}^2x\right)=0. $$ (54)

    The solution and the coefficient of the linear term are expanded in the forms in Eqs. (3) and (4)

    $$ \eta ={\omega}^2+{a}_1p+{a}_2{p}^2+..\dots, $$ (55)
    $$ x={x}_0+{x}_1p+{x}_2{p}^2+..\dots, $$ (56)

    By substituting Eqs. (55) and (56) into Eq. (54), the following linear differential equation system is obtained

    $$ {\ddot{x}}_0+{\omega}^2{x}_0=0,{x}_0(0)=l,{\dot{x}}_0(0)=0, $$ (57)
    $$ {\ddot{x}}_1+{\omega}^2{x}_1+{a}_1{x}_0+4{b}^2{\ddot{x}}_0{x}_0^2+4{b}^2{\dot{x}}_0^2{x}_0=0. $$ (58)

    Solving Eq. (57) for x0 gives

    $$ {x}_0(t)=l\cos \left(\omega t\right). $$ (59)

    Substituting Eq. (59) into Eq. (58) leads to

    $$ {\ddot{x}}_1+{\omega}^2{x}_1+{a}_1l\cos \left(\omega t\right)-4{b}^2{l}^3{\omega}^2{\cos}^3\left(\omega t\right)\\ \;\;\;\;+4{b}^2{l}^3{\omega}^2{\sin}^2\left(\omega t\right)\cos \left(\omega t\right)\\ \;\;\;\;=0. $$ (60)

    or

    $$ {\ddot{x}}_1+{\omega}^2{x}_1+\left({a}_1l-2{b}^2{l}^3{\omega}^2\right)\cos \left(\omega t\right)-2{b}^2{l}^3{\omega}^2\cos \left(3\omega t\right)\\ \;\;\;=0. $$ (61)

    No secular term in x1 requires that

    $$ {\mathrm{a}}_1l-2{b}^2{l}^3{\omega}^2=0, $$ (62)

    or

    $$ {a}_1=2{b}^2{l}^2{\omega}^2. $$ (63)

    The first-order approximate from Eq. (55) is

    $$ \eta ={\omega}^2+{a}_1={\omega}^2+2{b}^2{l}^2{\omega}^2, $$ (64)

    and hence from Eqs. (63) and (64), the following frequency is readily obtained

    $$ \omega =\sqrt{\frac{\eta }{1+2{b}^2{l}^2}}. $$ (65)

    Also, from Eq. (59), we obtain the velocity expression

    $$ {\dot{x}}_0(t)=- l\omega \sin \left(\omega t\right). $$ (66)

    The proposed methods are applied to investigate a steady-state roll motion of a ship with nonlinear damping and restoring moments. A comparison between the numerical results obtained by the ultraspherical wavelet method (UWM), the analytical results obtained by a new approach of the homotopy perturbation method (HPM), and numerical simulations obtained by the fourth-order Runge-Kutta method (numerical) is discussed.

    Roll angle decay curves in Example 1 are obtained by UWM (Eq. (29)) and the HPM (Eq. (34)). Figure 2 shows that these methods are in excellent agreement with the numerical solution obtained by RK4. Figure 3 also shows an excellent agreement between the UWM solution (Eq. (39)), the analytical solution (Eq. (40)), and RK4 for the roll angle decay of Example 2, Case 1. The curves in Figures 2 and 3 resemble the large amplitudes roll angles due to cubic damping moments and show that the motion continues indefinitely when only conservative forces act, and thus, the mechanical energy remains constant. Figure 4 shows that the velocity curves in Example 2 Case 1obtained by the modified HPM (Eq. (41)) and RK4 are identical. In Figures 5, 6, and 7, the curves of acceleration (Eq. (42)), restoring moment (Eq. (43)), and damping moment (Eq. (44)) are plotted, respectively, against time using the analytical HPM. Roll angle decay curves for Example 2, Cases 2 and 3, are shown, respectively, in Figures 8 and 9, while Figure 10 represents the roll angle decay curves for Example 3 using the proposed methods. Figure 11 depicts the displacement and velocity with weak nonlinearity and small amplitude. The envelopes of nonlinear ship roll motion, velocity, and acceleration are depicted in Figures 12 and 13. The frequency curves versus roll motion, velocity, and acceleration are given in Figures 14, 15, and 16.

    Figure  2  Comparison of UWM, HPM, and NM for roll angle decay curve in Example 1
    Download: Full-Size Img
    Figure  3  Comparison of UWM, HPM, and NM for roll angle decay curve in Example 2, Case 1
    Download: Full-Size Img
    Figure  4  Comparison between the HPM and NM for the velocity curve in Example 2, Case 1
    Download: Full-Size Img
    Figure  5  Plot of analytical expression of the acceleration curve in Example 2 Case 1
    Download: Full-Size Img
    Figure  6  Plot of restoring moment versus time, Eq. (43)
    Download: Full-Size Img
    Figure  7  Plot of damping moment versus time, Eq. (44)
    Download: Full-Size Img
    Figure  8  Comparison of UWM, HPM, and NM for roll angle decay curve in Example 2, Case 2
    Download: Full-Size Img
    Figure  9  Comparison of UWM, HPM, and NM for roll angle decay curve in Example 2, Case 3
    Download: Full-Size Img
    Figure  10  Comparison of UWM, HPM, and NM for roll angle decay curve in Example 3
    Download: Full-Size Img
    Figure  11  Displacement and velocity curves for the oscillator in Example 3 with weak nonlinearity and small-amplitude oscillations
    Download: Full-Size Img
    Figure  12  Roll angle decay curve with its envelope from Eq. (34)
    Download: Full-Size Img
    Figure  13  Roll angle decay curve with its envelope from Eq. (40)
    Download: Full-Size Img
    Figure  14  Roll motion vs. frequency for Example 3
    Download: Full-Size Img
    Figure  15  Velocity vs. frequency for Example 3
    Download: Full-Size Img
    Figure  16  Roll motion vs. frequency for Example 3
    Download: Full-Size Img

    In this paper, a mathematical model of the roll motion of ships with nonlinear damping and restoring moments and exciting moments is discussed. To predict the roll motions of ships in irregular or regular waves and to identify the damping and restoring moments in the model, an analytical solution using a modified form of the HPM and a numerical solution using ultraspherical wavelet-based method are obtained. The accuracy of the analytical and numerical solutions is confirmed by a direct comparison with the highly accurate and widely used fourth-order Runge-Kutta method. Analytical expressions for the roll angle, velocity, acceleration, and damping and restoring moments were also derived.

    The ultraspherical polynomials are special types of Jacobi polynomials that are associated with the real parameter $ \left(\alpha \gt \frac{-1}{2}\right) $. They are orthogonal polynomials on the interval (− 1, 1), with respect to the weight function $ w(x)={\left(1-{x}^2\right)}^{\alpha -\frac{1}{2}}, $ and defined by

    $$ {\int}_{-1}^1{\left(1-{x}^2\right)}^{\alpha -\frac{1}{2}}{C}_m^{\left(\alpha \right)}(x){C}_n^{\left(\alpha \right)}(x) dx=\left\{\begin{array}{cc}0,& m\ne n\\ {}{h}_n& m=n\end{array}\right. $$ (67)

    where

    $$ {h}_n=\frac{\pi {2}^{1-2\alpha}\varGamma \left(n+2\alpha \right)}{n!\left(n+\alpha \right)\varGamma {\left(\alpha \right)}^2}, $$ (68)

    are the eigenfunctions of the following singular Sturm-Liouville equation

    $$ \left(1-{x}^2\right){\ddot{\phi}}_m(x)-\left(2\alpha +1\right)x{\ddot{\phi}}_m(x)+m\left(m+2\alpha \right){\phi}_m(x)=0. $$ (69)

    The following integral formula and the theorem that follows are needed to establish the convergence of the expansion of the ultraspherical wavelets

    $$ \int {C}_n^{\left(\alpha \right)}(x)w(x) dx=-\frac{2\alpha {\left(1-{x}^2\right)}^{\alpha +\frac{1}{2}}}{n\left(n+2\alpha \right)}{C}_{n-1}^{\left(\alpha +1\right)}(x),n\ge 1. $$ (70)

    Theorem 3 The following inequality holds for ultraspherical polynomials

    $$ {\left(\sin \theta \right)}^{\alpha}\left|{C}_n^{\left(\alpha \right)}\left(\cos \theta \right)\right| \lt \frac{2^{1-\alpha}\varGamma \left(n+\frac{3\alpha }{2}\right)}{\varGamma \left(\alpha \right)\varGamma \left(n+1+\frac{\alpha }{2}\right)},0\le \theta \le \pi, 0 \lt \alpha \lt 1.\kern1em $$ (71)

    The shifted ultraspherical polynomials are defined on [0, 1] by

    $$ {\overset{\sim }{C}}_n^{\left(\alpha \right)}(x)={C}_n^{\left(\alpha \right)}\left(2x-1\right). $$ (72)

    All properties of ultraspherical polynomials remain valid for the shifted polynomials.

    The orthogonality relation for $ {\overset{\sim }{C}}_n^{\left(\alpha \right)}(x) $ with respect to the weight function $ \overset{\sim }{\omega }(x)={\left(x-{x}^2\right)}^{\alpha -\frac{1}{2}} $ is given by

    $$ \int {\left(x-{x}^2\right)}^{\alpha -\frac{1}{2}}{\overset{\sim }{C}}_m^{\left(\alpha \right)}(x){\overset{\sim }{C}}_n^{\left(\alpha \right)}(x) dx\\ \;\;\;=\left\{\begin{array}{cc}0,& m\ne n\\ {}\frac{\pi {2}^{1-4\alpha}\Gamma \left(n+2\alpha \right)}{n!\left(n+\alpha \right){\left(\Gamma \left(\alpha \right)\right)}^2},& m=n.\end{array}\right. $$ (73)

    For more properties of ultraspherical polynomials, see Rainville (1960).

    Consider the nonlinear differential equation

    $$ A(u)-f(r)=0,r\in \varOmega, $$ (74)

    with the boundary condition

    $$ B\left(u,\frac{du}{dr}\right)=0,r\in \Gamma, $$ (75)

    where A, B, f(r), and Γ are a general differential operator, a boundary operator, a known analytical function, and the boundary of the domain Ω, respectively. Expressing A(u) as the sum of linear (L) and nonlinear (N) parts, Eq. (74) becomes

    $$ L(u)+N(u)-f(r)=0.\kern1.5em $$ (76)

    The homotopy technique begins by defining

    v(r, p) : Ω × [0, 1] → R, such that

    $$ H\left(v,p\right)=\left(1-p\right)\left[L(v)-L\left({u}_0\right)\right]+p\left[A(u)-f(r)\right]=0, $$ (77)

    where p ∈ [0, 1] is an embedding parameter and u0 is an initial approximation of Eq. (74) that satisfies boundary conditions (Eq. 75). Evidently, Eq. (77) implies that

    $$ H\left(v,0\right)=L(v)-L\left({u}_0\right)=0,\kern0.5em $$ (78)
    $$ H\left(v,1\right)=A(v)-f(r)=0. $$ (79)

    As p changes from 0 to 1, v(r, p) changes from u0 to ur, a process known as a homotopy. The solution of Eq. (77) may be expressed in terms of a power series in the form:

    $$ v={v}_0+p{v}_1+{p}^2{v}_2+\cdots . $$ (80)

    An approximate solution to Eq. (77) is given by the following:

    $$ u=\underset{p\to 1}{\lim }v={v}_0+{v}_1+{v}_2+\cdots . $$ (81)
    Acknowledgements: The authors are thankful to Shri J. Ramachandran, Chancellor, Col. Dr. G. Thiruvasagam, Vice-Chancellor, Academy of Maritime Education and Training (AMET), Deemed to be University, Chennai, for their support.
  • Figure  1   Schematic diagram of ship showing the six degrees of freedom

    Download: Full-Size Img

    Figure  2   Comparison of UWM, HPM, and NM for roll angle decay curve in Example 1

    Download: Full-Size Img

    Figure  3   Comparison of UWM, HPM, and NM for roll angle decay curve in Example 2, Case 1

    Download: Full-Size Img

    Figure  4   Comparison between the HPM and NM for the velocity curve in Example 2, Case 1

    Download: Full-Size Img

    Figure  5   Plot of analytical expression of the acceleration curve in Example 2 Case 1

    Download: Full-Size Img

    Figure  6   Plot of restoring moment versus time, Eq. (43)

    Download: Full-Size Img

    Figure  7   Plot of damping moment versus time, Eq. (44)

    Download: Full-Size Img

    Figure  8   Comparison of UWM, HPM, and NM for roll angle decay curve in Example 2, Case 2

    Download: Full-Size Img

    Figure  9   Comparison of UWM, HPM, and NM for roll angle decay curve in Example 2, Case 3

    Download: Full-Size Img

    Figure  10   Comparison of UWM, HPM, and NM for roll angle decay curve in Example 3

    Download: Full-Size Img

    Figure  11   Displacement and velocity curves for the oscillator in Example 3 with weak nonlinearity and small-amplitude oscillations

    Download: Full-Size Img

    Figure  12   Roll angle decay curve with its envelope from Eq. (34)

    Download: Full-Size Img

    Figure  13   Roll angle decay curve with its envelope from Eq. (40)

    Download: Full-Size Img

    Figure  14   Roll motion vs. frequency for Example 3

    Download: Full-Size Img

    Figure  15   Velocity vs. frequency for Example 3

    Download: Full-Size Img

    Figure  16   Roll motion vs. frequency for Example 3

    Download: Full-Size Img
  • Abualrub T, Abukhaled M (2015) Wavelets approach for optimal boundary control of cellular uptake in tissue engineering. Int J Comp Math 92(7): 1402–1412. https://doi.org/10.1080/00207160.2014.941826
    Abualrub T, Abukhaled M, Jamal B (2018) Wavelets approach for the optimal control of vibrating plates by piezoelectric patches. J Vibr Cont 24(6): 1101–1108. https://doi.org/10.1177/1077546316657781
    Abukhaled M (2013) Variational iteration method for nonlinear singular two-point boundary value problems arising in human physiology. J Math: ID: 720134. https://doi.org/10.1155/2013/720134
    Abukhaled M (2017) Green's function iterative approach for solving strongly nonlinear oscillators. J Comp Nonlinear Dyn 12(5): 051021. https://doi.org/10.1115/1.4036813
    Agarwal D (2015) A study on the feasibility of using fractional differential equations for roll damping models. Master thesis, Virginia Polytechnic Institute and State University. http://hdl.handle.net/10919/52959
    Aloisio G, Felice F (2006) PIV analysis around the bilge keel of a ship model in a free roll decay. In: XIV Congresso Nazionale AI VE. LA., Rome, Italy
    Bassler CC, Carneal JB, Atsavapranee P (2007) Experimental investigation of hydrodynamic coefficients of a wavepiercing tumble home hull form. In: Proceeding of the 26th International Conference on Offshore Mechanics and Arctic Engineering, San Diego, USA
    Cardo A, Francescutto A, Nabergoj R (1981) Ultraharmonics and subharmonics in the rolling motion of a ship: steadystate solution. Int Ship building Progress 28: 234–251 https://doi.org/10.3233/ISP-1981-2832602
    Cardo A, Francescutto A, Nabergoj R (1984) Nonlinear rolling response in a regular sea. Int Shipbuild Prog 31(360): 204–206. https://doi.org/10.3233/ISP-1984-3136002
    Chen CL, Liu YC (1998) Solution of two-point boundary-value problems using the differential transformation method. J Optimization Theory Appl 99: 23–35 https://doi.org/10.1023/A:1021791909142
    Comini G, Del Guidice S, Lewis RW, Zienkiewicz OC (1970) Finite element solution of non-linear heat conduction problems with special reference to phase change. Num Meth Engg 8: 613–624
    Demirel H, Alarcin F (2016) LMI-based H2 and H1 state-feedback controller design for fin stabilizer of nonlinear roll motion of a fishing boat. Brodogradnja: Teorija i praksa brodogradnje i pomorske tehnike 67(4): 91–107. https://doi.org/10.21278/brod67407
    Demirel H, Dogrul A, Sezen S, Alarcin F (2017) Backstepping control of nonlinear roll motion for a trawler with fin stabilizer. Transactions of the Royal Institution of Naval Architects Part A: International Journal of Maritime Engineering, Part A2 159, pp. 205-212. doi: 10.3940/rina.ijme.2017.a2.420
    Doha EH, Abd-Elhameed WM, Youssri YH (2016) New ultraspherical wavelets collocation method for solving 2nth order initial and boundary value problems. J Egyp Math Soc 24: 319–327. https://doi.org/10.1016/j.joems.2015.05.002
    Froude W (1861) On the rolling of ships. Transactions of the Institution of Naval Architects 2, 1861, 180–229.
    He JH (1999) Homotopy perturbation technique. Comp Meth Appl Mech Engg 178(4): 257–262. https://doi.org/10.1016/S0045-7825(99)00018-3
    Himeno Y (1981) Prediction of ship roll damping-a state of the art. DTIC Document, Ann Arbor.
    Huang BG, Zou ZJ, Ding WW (2018) Online prediction of ship roll motion based on a coarse and fine tuning fixed grid wavelet network. Ocean Eng 160: 425–437. https://doi.org/10.1016/j.oceaneng.2018.04.065
    Ikeda Y, Ali B, Yoshida H (2004) A roll damping prediction method for a FPSO with steady drift motion. Proceeding of the 14th International Conference on Offshore and Polar Engineering Conference, Toulon, France 676–681.
    Jang TS (2011) Non-parametric simultaneous identification of both the nonlinear damping and restoring characteristics of nonlinear systems whose dampings depend on velocity alone. Mech Syst Signal Process 25(4): 1159–1173. https://doi.org/10.1016/j.ymssp.2010.11.002
    Jang TS (2013) A method for simultaneous identification of the full nonlinear damping and the phase shift and amplitude of the external harmonic excitation in a forced nonlinear oscillator. Comput Struct 120: 77–85. https://doi.org/10.1016/j.compstruc.2013.02.008
    Jang TS, Choi HS, Han SL (2009) A new method for detecting non-linear damping and restoring forces in non-linear oscillation systems from transient data. Int J Non-Linear Mech 44(7): 801–808. https://doi.org/10.1016/j.ijnonlinmec.2009.05.001
    Kato H (1965) Effects of bilge keels on the rolling of ships. J. Soc. Nav. Archit. Jpn. 117: 93–101
    Khuri SA, Abukhaled M (2017) A semi-analytical solution of amperometric enzymatic reactions based on Green's functions and fixed point iterative schemes. J Electroanalytical Chem 792: 66–71. https://doi.org/10.1016/j.jelechem.2017.03.015
    Kianejad SS, Enshaei H, Du_y J, Ansarifard N (2019) Prediction of a ship roll added mass moment of inertia using numerical simulation. Ocean Eng 173: 77–89. https://doi.org/10.1016/j.oceaneng.2018.12.049
    Lavrov A, Rodrigues JM, Gadelho JFM, GuedesSoares C (2017) Calculation of hydrodynamic coe_cients of ship sections in roll motion using Navier-Stokes equations. Ocean Eng 133: 36–46. https://doi.org/10.1016/j.oceaneng.2017.01.027
    Liao S (1997) Homotopy analysis method: a new analytical technique for nonlinear problems. Commun Nonlinear Sci Numer Simul 2(2): 95–100. https://doi.org/10.1016/S1007-5704(97)90047-2
    Liao S (2004) On the homotopy analysis method for nonlinear problems. Appl Math Comput 47(2): 499–513. https://doi.org/10.1016/S0096-3003(02)00790-7
    Liao SJ, Chwang AT (1998) Application of homotopy analysis method in nonlinar oscillations. J Appl Mech 65(4): 914–922. https://doi.org/10.1115/1.2791935
    Lihua L, Peng Z, Songtao Z, Ming J, Jia Y (2018) Simulation analysis of fin stabilizer on ship roll control during turning motion. Ocean Eng 164: 733–748. https://doi.org/10.1016/j.oceaneng.2018.07.015
    Lu J (2007) Variational iteration method for solving a nonlinear system of second-order boundary value problems. Comp Math Appl 54: 1133–1138. https://doi.org/10.1016/j.camwa.2006.12.060
    Oliveira AC, Fernandes AC (2012) An empirical nonlinear model to estimate FPSO with extended bilge keel roll linear equivalent damping in extreme seas. In: Proceeding of the 31st International Conference on Ocean, Offshore and Arctic Engineering, (OMAE). Paper OMAE2012-83360, Rio de Janeiro, Brazil.
    Oliveira AC, Fernandes AC (2014) The nonlinear roll damping of a FPSO hull. J Off shore Mech Arct Eng 136: 1–10. https://doi.org/10.1115/1.4025870
    Rainville ED (1960) Special functions. The Macmillan Co., New York, 1960. MR 0107725.
    Razzaghi M, Yousef S (2002) Legendre wavelets method for constrained optimal control problems. Math Met Appl Sci 25: 529–539. https://doi.org/10.1002/mma.299
    Sadek I, Abualrub T, Abukhaled M (2007) A computational method for solving optimal control of a system of parallel beams using Legendre wavelets. Math Comp Modell 45: 1253–1264. https://doi.org/10.1016/j.mcm.2006.10.008
    Salai Mathi Selvi M, Hariharan G (2016) Wavelet based analytical algorithm for solving steady-state concentration in immobilized isomerase of packed bed reactor model. J Memb Biol. 249(4): 559–586 https://doi.org/10.1007/s00232-016-9905-2
    Salai Mathi Selvi M, Hariharan G (2017) An improved method based on Legendre computational matrix method for time dependent Michaelis-Menten enzymatic reaction model arising in mathematical chemistry. Proc Jangjeon Mathe Soc. 20(3): 483–503. https://doi.org/10.17777/pjms2017.20.3.483
    Salai Mathi Selvi M, Hariharan G, Kannan K (2017) A reliable spectral, method to reaction-diffusion equations in entrapped-cell photobioreaction packed with gel grnules using chebyshev wavelets. J Memb Biol 250: 663–670 https://doi.org/10.1007/s00232-017-0001-z
    Tanaka N (1961) A study on the bilge keels. J Soc Nav Archit Jpn 109: 205–212
    Tavassoli Kajani M, Hadi Vencheh A (2004) Solving linear integro-differential equation with Legendre wavelets Int. J Comput Math 81(6): 719–726. https://doi.org/10.1080/00207160310001650044
    Yeung RW, Cermelli C, Liao SW (1997) Vorticity fields due to rolling bodies in a free surface - experiment and theory. Proceeding of the 21st symposium on naval hydrodynamics Trondheim Norway.
    Zheng X, Wei Z (2016) Estimates of approximation error by Legendre wavelet. Appl Math 7:694–700. https://doi.org/10.4236/am.2016.77063
    Salomi RJ, Sylvia SV, Rajendran L, Abukhaled M, (2020) Electric potential and surface oxygen ion density for planar, spherical and cylindrical metal oxide grains. Sensors and Actuators B: Chemical 321: 128576. https://doi.org/10.1016/j.snb.2020.128576
WeChat click to enlarge
Figures(16)
Publishing history
  • Received:  18 November 2019
  • Accepted:  12 May 2020

目录

    /

    Return
    Return