Journal Search Engine
Search Advanced Search Adode Reader(link)
Download PDF Export Citaion korean bibliography PMC previewer
ISSN : 1229-3431(Print)
ISSN : 2287-3341(Online)
Journal of the Korean Society of Marine Environment and Safety Vol.26 No.5 pp.437-444
DOI : https://doi.org/10.7837/kosomes.2020.26.5.437

Active Stabilization for Surge Motion of Moored Vessel in Irregular Head Waves

Sang-Do Lee*, Truong Ngoc Cuong**, Xiao Xu**, Sam-Sang You***
*Professor, Division of Navigation & Information System, Mokpo National Maritime University, Mokpo, Republic of Korea
**Doctor candidate, Department of Logistics Engineering, Korea Maritime and Ocean University, Busan, Republic of Korea
***Professor, Division of Mechanical Engineering, Korea Maritime and Ocean University, Busan, Republic of Korea

* First author : oksangdo@mmu.ac.kr, 061-240-7257


Corresponding Author : ssyou@kmou.ac.kr, 051-410-4366

June 12, 2020 August 24, 2020 August 28, 2020

Abstract


This study was focused on the stabilization of surge motions of a moored vessel under irregular head seas. A two-point moored vessel shows strong non-linearity even in regular sea, owing to its inherent non-linear restoring force. A long-crested irregular wave is subjected to the vessel system, resulting in more complex nonlinear behavior of the displacement and velocities than in the case of regular waves. Sliding mode control (SMC) is implemented in the moored vessel to control both surge displacement and surge velocity. The SMC can provide a closed-loop system with performance and robustness against parameter uncertainties and disturbances; however, chattering is the main drawback for implementing SMC. The goal of minimizing the chattering and state convergence with accuracy is achieved using a quasi-sliding mode that approximates the discontinuous function via a continuous sigmoid function. Numerical simulations were conducted to validate the effectiveness of the proposed control algorithm.



불규칙 선수파랑 중 계류된 선박의 전후동요 제어

이 상도*, 트롱 엥곡 쿠옹**, 서 효**, 유 삼상***
*목포해양대학교 항해정보시스템학부 교수
**한국해양대학교 대학원 박사과정
***한국해양대학교 기계공학부 교수

초록


본 논문에서는 불규칙 선수파랑 조건에서 선박계류시스템의 전후동요를 억제하는 연구를 수행하였다. 두점식 계류시스템은 비 선형복원력의 특성으로 인하여 규칙파 조건에서도 강한 비선형 응답 특성을 보인다. 종방향 불규칙 장파가 외란으로 선박계류시스템에 작용하게 되면 규칙파 외력이 입사하는 경우보다 변위와 속도에서 더욱 복잡한 비선형 거동이 발생한다. 계의 종동요 변위와 속도를 동 시에 억제하기 위하여 슬라이딩모드 제어법(SMC)을 적용하였다. SMC는 매개변수의 불확실성 및 외란에 대한 강인성을 갖는 폐루프 시스 템을 제공하지만, 채터링은 이 제어법을 사용할 때 큰 단점이 된다. 본 연구에서는 준 슬라이딩모드의 시그모이드 함수를 이용하여 불규 칙 해양파의 외란 조건에서 채터링을 줄이고. 수렴의 정확성에 도달하는 목표를 달성하였다. 제어법의 유효성은 수치시뮬레이션을 통해 증명하였다.



    1. Introduction

    For the maritime safety, marine vessels in the ocean waves should be controlled even in the harsh circumstances. The dynamical responses of the vessel system are varied according to the operating conditions of exciting ocean waves. In order to obtain a comfortable working environment, the ship should be suppressed with an appropriate control method. The various vessel systems contain the inherent strong non-linearity of both two-point moored system (Lee and You, 2018a;Lee and You, 2018c) and rolling vessel system (Lee and You, 2018b;Lee et al., 2020).

    Moored vessel can mainly be implemented as a multi-point and a single-point mooring system. These mooring systems represent complex dynamical behavior including co-existing periodic (harmonic as well as sub-harmonic oscillation) and aperiodic (quasi-periodic and chaotic) motions (Gottlieb and Yim, 1992) due to nonlinear features of restoring force (Umar et al., 2010). Also, the stability of the moored vessel is governed by strong sensitivity to initial conditions (Ellermann, 2005;Banik and Datta, 2010). Even in case that the harmonic excitation is acting on the system, a two-point moored vessel represents various nonlinear motions such as chaotic motion and limit cycles due to the effect of the strong restoring term (Lee and You, 2018a). Thus, the two-point moored vessel actually resembles the duffing oscillator.

    The nonlinear duffing oscillator can be seen as a two-point moored vessel subjected to harmonic excitations. It represents the harmonic or sub-harmonic to chaotic motion (Ellermann, 2005). Duffing type motion has been restrained using several control mechanism such as PID-based controller (Loria et al., 1998), sliding mode controller (Kuo et al., 2008), modified super-twisting algorithm (Khadra, 2016), super-twisting control algorithm (Lee and You, 2018c) and time-delay state feedback (Mitra et al., 2017). Their reports show the mitigation of chattering which is actually undesirable oscillations with finite amplitude with frequency.

    The sliding mode control (SMC) is one of the effective control mechanisms for nonlinear dynamical systems ensuring performance and robustness against parameter uncertainties and exogenous disturbances. Sliding surface and active controller are designed to stay the state trajectories on the time-varying sliding surface in finite time. However, one major drawback of SMC is the chattering phenomenon, which results in the harmful failure of actuators with fatigue loads due to the high frequency of switching actions. These inherent problems of chattering in SMC can be improved by suitable measures such as a sigmoid function or the second-order sliding mode. It is difficult to eliminate totally the chattering in second-order sliding mode because it contains the signum function multiplied by some bounds of uncertainties in advance (Edwards and Shtessel, 2016). Despite shortcomings of chattering, SMC has been prominently appeared in the stage of engineering fields due to its strengths of invariability to external disturbances. Recently, higher-order SMC such as super-twisting and adaptive super-twisting algorithm (Xu et al., 2020) have been reported in controlling marine vessel system (Lee and You, 2018c;Lee et al., 2019;Lee et al., 2020). However, since there are many parameters in higher-order SMC, it requires much effort and time to apply in practical cases such as complicated ship models and irregular head waves.

    This paper investigates the active stabilization of irregular surge oscillation for a two-point moored vessel under irregular head seas. Surge motions governed by a single-degree-of-freedom (SDOF) equation. Quasi-SMC is taken into account for ensuring the control performance against uncertainties. The rest of the paper is structured as follows. Section 2 briefly introduces the system dynamics of a two-point moored vessel. In Section 3, Quasi-SMC has been applied to suppress irregular surge oscillation. Section 4 presents numerical simulations for the effectiveness of the proposed control mechanism. Finally, the conclusion will be made in Section 5.

    2. Mathematical Formulation

    2.1 Two-point Moored Vessel

    The two-point moored vessel under irregular head waves governed by SDOF can be simplified by the following nonlinear differential equation (Banik and Datta 2010;Lee and You, 2018a)

    x ¨ + γ 1 x ˙ + R ( x ) = F ¯ ( x ¨ , x ˙ , t ) / M
    (1)

    with

    R ( x ) = Θ [ x + χ m r s g n ( x ) ] { 1 1 + χ m r 2 1 1 + [ x + χ m r s g n ( x ) ] 2 } γ 1 = C / M , Θ = X s / M , χ m r = b 1 / b 2 , s g n ( x ) = { + 1 for x > 0 0 for x = 0 1 for x < 0 }

    where M = m υ + m a ; m υ denotes the original vessel mass (kg) and ma is the added mass; C represents the damping coefficient (N.s/m); R(x) is the restoring force; F represents external excitation force (N); χmr is the non-dimensional values related to the distance of each moorings; sgn refers to signum function; Xs is the stiffness (N/m) including elastic cable force defined as X s = 2 E c f b 2 2 + b 1 2 ; Ecf is the elastic cable force; bi (i = 1,2) designates distances between the center of vessel system and each mooring lines; over-dot denotes a derivative of the variable with respect to time; x, x ˙ and x ¨ represent the surge displacements (m), surge velocity (m/s) and surge acceleration (m/s2), respectively (Banik and Datta 2010;King and Yim, 2007).

    The restoring force assumes linear elastic behavior so that the non-linearity is strictly due to the geometric configuration of the system (King and Yim, 2007). Particularly, the restoring force can be approximated by the least square representation as Φ n x n , : n = 1 , 3 , N where Φn describes the exciting frequency. The governing system non-linearity N = 3 is mainly cubic term (x3). Polynomials of various orders are employed within experimental range. In practice, the restoring force R(x) can be expressed as linear and cubic polynomial form as υ1x + υ2x3 (Banik and Datta 2010;Lee and You, 2018a). The nonlinear chaos motion mainly depends on the inherent restoring force term and the initial condition in case of regular excitation. Assuming the vessel mass has just unit mass for simplicity, it can be scaled into the vessel mass or f s e a = F ¯ / M ( N / k g ) . Thus, the two-point moored vessel under irregular head seas can be written as

    x ¨ + γ 1 x ˙ + υ 1 x + υ 2 x 3 = f s e a ( t )
    (2)

    2.2 Irregular head waves

    The two-point moored vessel shows rich nonlinear dynamical behaviors due to its strong non-linearity even in a regular wave. In this section, we introduces a long-crested irregular sea to show the performance of the proposed control algorithm. A long-crested sea has the unidirectional wave crests with varying separation but remaining parallel to each other (Perez, 2005). When the storm blows for a long time, a fully developed sea is created. Such long waves form a wave spectrum. We employed the Pierson-Moskowitz spectrum (PM spectrum) on the two-point moored vessel as the real sea states. PM spectrum is a one-parameter wave spectral formulation for fully developed wind-generated seas based on the wave spectra information of North Atlantic Ocean as follows (Fossen, 2011)

    S ( ω ) = S a ω 5 exp ( S b ω 4 ) ( m 2 s )
    (3)

    with Sa = 8.1 × 10-3g2 and Sb = 0.74(g/Sc)4, where Sc means the wind speed at the height 19.4 m.

    Assuming the waves can be represented by Gaussian random process and S(ω) is narrow-banded, Sb can be rewritten as 0.0323 ( g / H s ) 2 = 3.11 / H s 2 where Hs is the significant wave height. The sea surface elevation of a long-crested irregular sea can be given as the sum of harmonic wave components i as the following equation

    ξ ( t ) = i = 1 N W i ¯ sin ( ω i t ν i x + κ i )
    (4)

    where Wi represent the wave amplitude defined as 2 S ( ω i ) Δ ω ; ωi is the circular frequency; νi is the wave number; κi is the random phase angle of wave component with uniform distribution in [0, 2π]. (Perez, 2005;Fossen, 2011).

    2.3 Robust control synthesis

    By considering the excitation fsea as a time-varying disturbance, the control system represents a forced surge system by adding the actuation input u.

    x ¨ + γ 1 x ˙ + υ 1 x + υ 2 x 3 = d + u
    (5)

    where d represents the total amounts of disturbances such as waves, wind, currents, and any impacts in longitudinal direction.

    An active controller can be designed to achieve the satisfactory anti-surging effect for marine vessel. For the state variables denoted as x = x1 and x ˙ = x2, the governing equation can be described as into state-space representation as follows

    x ˙ = T a x + f ( x , t ) + T b u + T c d
    (6)

    where state vector, system matrices, and nonlinear term are given by

    x = [ x 1 x 2 ] , T a = [ 0 1 υ 1 γ 1 ] , T b = T c = [ 0 1 ] and f ( x , t ) = [ 0 υ 2 x 1 3 ]
    (7)

    As shown in Fig. 1, the proposed controller is employed to handle the error responses of the two-point moored vessel under irregular head waves.

    By introducing the surge tracking error e = x 1 x d between the actual output x1 and desired value xd, the variables on the sliding surface are considered as

    s = λ 1 e ˙ + λ 2 e
    (8)

    where λ1 and λ2 are positive constants which are related to the control accuracy and rate of convergence for the moored vessel. From (6) and (7), the sliding mode dynamics can be derived as

    s ˙ = P 0 + λ 1 ( d + u )
    (9)

    where

    P 0 = λ 1 ( γ 1 x 2 υ 1 x 1 υ 2 x 1 3 x ¨ d ) + λ 2 ( x 2 x ˙ d ) .

    The control objective is to design a suitable controller to ensure the robust stabilization of sliding mode dynamics asymptotically or s ˙ = 0 under uncertainties. In order to achieve the goal, the active controller is designed with the equivalent part and switching part as follows

    u = P 0 λ 1 + 1 λ 1 u s w
    (10)

    where P 0 λ 1 and 1 λ 1 u s w represent the equivalent controller and switching controller, respectively. Thus, the sliding mode dynamics in equation (9) can be rewritten as

    s ˙ = λ 1 d + u s w
    (11)

    In order to solve the chattering phenomenon, Quasi-SMC technique is implemented where the discontinuous term on the controller is replaced by a continuous sigmoid function s i g n ( s ) s | s | + (Shtessel et al., 2014;Cuong et al., 2020).

    Assume that the disturbance is bounded by a positive constant ρ, satisfying | d | ρ . The switching controller is designed as

    u s w = β s | s | + λ 1 ρ
    (12)

    where β is a sufficiently large positive constant, and is a free parameter which can be chosen arbitrarily by the system designer.

    Theorem 1. The surge displacement x1 of the moored vessel (6) are regular to track the constant reference signal xd asymptotically as t→∞ for all initial x1(0) with the sliding surface (8) and the sliding mode control law (10) and (12), where β, λ are positive constants.

    Proof: To prove the closed-loop stability, the quadratic Lyapunov function candidate is considered as follow.

    V ( s ) = 1 2 s 2
    (13)

    The time-derivative of V is obtained as

    V ˙ = s s ˙
    (14)

    Substituting (11) into (14) leads to the following relation

    V ˙ = s ( λ 1 d + u s w )
    (15)

    From (12) and (15), one yields,

    V ˙ = s ( λ 1 d β s | s | + λ 1 ρ ) λ 1 | s | ( | d | ρ ) β s 2 | s | +
    (16)

    Since β, λ1 and ρ are positive, it is obvious that V ˙ is a negative definite function, i.e., V ˙ ≤0. Thus, the time convergence of tracking error lim t e ( t ) = 0 is guaranteed by Lyapunov’s stability theory. This completes the proof.

    3. Simulation tests

    Recent developments in sliding mode controls have been chosen to stabilize the surge oscillation of a two-point moored vessel subjected to the irregular head waves. Table 1 shows the main parameters for the moored vessel system (Shah et al., 2005). The dynamical model is a vessel type structure with two mooring lines in the fore and aft position. When the irregular head waves are excited about this structure in the longitudinal direction, the system exhibits the irregular complex oscillations on the displacement and the velocity. So, the proposed SMC is aiming to solve convergence problems of the two states and chattering reduction within finite time simultaneously.

    The simulation results of the controlled moored vessel will be dependent on the control parameters. In the simulation, some parameters for control synthesis are given as λ1 = 0.1, λ2 = 0.04, β = 30 and = 0.01, in which the parameters are determined through many tests. In order to find the optimum parameters, we set the minimum vale for reference point. Then, we increase the each parameter gradually. The parameters λ1 and λ2 in sliding mode variables affect the convergence rate and the chattering alleviation for surge responses. In addition, the parameters β and are also set through many trials guaranteeing the chattering mitigation. The initial conditions of the surge dynamics are selected as [ x 0 , x ˙ 0 ] = [0.5 (m), 0.1 (m/s)].

    Fig. 2 shows time history curves for a long-crested irregular sea with significant wave height (Hs) 4.82 m and wind speed 15 m/s. Complex irregular wave elevations can be observed from 0 to 100 seconds. These waves excited on the moored vessel in the longitudinal direction as disturbances. When no control action is employed, the moored vessel shows various irregular motions depending on initial conditions even in regular excitations (Lee and You, 2018a;Lee and You, 2018c).

    In this paper, irregular head waves have been introduced such as PM spectrum for a real cases. The moored vessel is still in a state of nonlinear surge oscillations with irregular excitations. An active controller has been initiated to suppress the irregular oscillations of two-point moored vessel due to the head seas at 2000 seconds.

    In order to compare the effectiveness of active control input, the uncontrolled case and controlled case are compared in Figs. 3 and 5. In addition, Figs. 3 and 4 depict the responses of surge displacement and its error, respectively. The surge motion depicts the irregular oscillations and the responses riding the waves before control action due to the effect of strong restoring and the disturbances. It can be observed that the mooring vessel stays on the near zero position with the help of an active controller. Tracking error is a measure of deviation from the benchmark, in which tracking error of surge displacement marked less than 0.5 m under control action. So, the control action achieves the design goal of providing low tracking error even in heavy weather. Simultaneously, as shown in Figs. 5 and 6, the state of surge velocity is effectively suppressed. It is noted that a peak of velocity can be happened when the control action starts. The trade-off is necessary for control performance and control gain. Also, the errors of surge velocity are decreased to near zero under control shown in Fig. 6.

    Fig. 7 represents the control activity using the proposed quasi-SMC, in which discontinuous high-frequency phenomenon has been significantly reduced under the proposed quasi-SMC. The switching part of the control law in equation (10) helps to drive the state trajectories on the sliding surface regardless of disturbances. Due to the irregular head waves, the harmonic curve can hardly be seen in the control input. Chattering is almost removed with the help of a quasi-SMC which approximates the discontinuous function by a continuous sigmoid function in equation (12). This control law can be applied to some actuators to propel or stop the ship effectively. Since the discontinuous control causes the failure of actuators with fatigue loads, chattering alleviation is important when the SMC is employed to stabilize surge oscillation of the moored vessel. The control action is highly related to the role of the sliding variable in Fig. 8. As described in SMC theory in equations (8) and (9), the sum of position error and velocity error decreases to zero in finite time using the proposed controller.

    In previous reports by Lee et al. (2020) and Shtessel et al. (2012), the approach of sigmoid function is vulnerable to the high disturbances since the s-trajectories are constrained in the vicinity of the sliding surface rather than on the surface. As see in Figs. 4 and 6, there exists some errors after control action. However, due to its simplicity compared to the control effort, the proposed quasi-SMC can be described as one of the effective mechanisms with respect to suppression of the two-point moored vessel under irregular head waves such as a long-crested irregular sea.

    4. Conclusions

    The quasi-SMC scheme has been presented for suppressing the complex surge oscillations of a two-point moored vessel under irregular head waves such as a long-crested sea. Due to the inherent characteristics of restoring forces, it is hard to expect the moored vessel to ensure the comfortable working place until the control action activates. For the purpose of safety on the moored vessel under real sea conditions, one of the effective control schemes is considered.

    The proposed nonlinear control system has exhibited the strength of SMC synthesis providing control accuracy, robustness, and simple implementation against disturbances. Simulation results revealed that the state variables of surge displacement and velocity of the two-point moored vessel eventually becomes near zero with active control action. The chattering of the SMC algorithm has been removed by the implementation of a quasi-sliding mode approximating the sign function by the continuous sigmoid function. Finally, this work has been focusing on controlling surge motions for a simple SDOF moored vessel under irregular head waves and uncertainty.

    Figure

    KOSOMES-26-5-437_F1.gif

    Block diagram of surge control for a two-point moored vessel under irregular head waves.

    KOSOMES-26-5-437_F2.gif

    Time history curves for a long-crested irregular sea: (a) PM spectrum (b) wave elevation.

    KOSOMES-26-5-437_F3.gif

    Comparison of surge suppression with/without controller (control action starts at t = 2000 sec).

    KOSOMES-26-5-437_F4.gif

    Error of surge displacement under control action.

    KOSOMES-26-5-437_F5.gif

    Comparison of surge rate suppression with/without controller (control action starts at t = 2000 sec).

    KOSOMES-26-5-437_F6.gif

    Error of surge velocity under control action.

    KOSOMES-26-5-437_F7.gif

    Control activities against irregular head waves (starting at t = 2000 sec).

    KOSOMES-26-5-437_F8.gif

    Time history curve for sliding variable before and after control action.

    Table

    Model parameters of two-point moored vessel

    Reference

    1. Banik, A. K. and T. K. Datta (2010), Stability Analysis of Two-Point Mooring System in Surge Oscillation, Journal of Computational and Nonlinear Dynamics, Vol. 5, pp. 1-8 (021005).
    2. Cuong, T. N. , X. Xu, S. D. Lee, and S. S. You (2020), Dynamic analysis and management optimization for maritime supply chains using nonlinear control theory, Journal of International Maritime Safety, Environmental Affairs, and Shipping Vol. 4, No. 2, pp. 48-55.
    3. Edwards, C. and Y. B. Shtessel (2016), Adaptive Continuous Higher Order Sliding Mode Control, Automatica, Vol. 65, pp. 183-190.
    4. Ellermann, K. (2005), Dynamics of a Moored Barge under Periodic and Randomly Disturbed Excitation, Ocean Engineering, Vol. 32, pp. 1420-1430.
    5. Fossen, T. I. (2011), Handbook of Marine Craft Hydrodynamics and Motion Control, John Wiley & Sons Ltd, pp. 192-211.
    6. Gottlieb, O. and S. C. S. Yim (1992), Nonlinear Oscillations, Bifurcations and Chaos in a Multi-Point Mooring System with a Geometric Nonlinearity, Applied Ocean Research, Vol. 14, pp. 241-257.
    7. Khadra, F. A. (2016), Super-Twisting Control of the Duffing- Holmes Chaotic System, International Journal of Modern Nonlinear Theory and Application, Vol. 5, pp. 160-170.
    8. King, P. E. and S. C. Yim (2007), Stochastic Control of Sensitive Nonlinear Motions of an Ocean Mooring System, Journal of Offshore Mechanics and Arctic Engineering, Vol. 129, pp. 29-38.
    9. Kuo, T. C. , Y. J. Huang, C. H. Chang, and C. Y. Chen (2008), Robust Suppression Sliding Mode Control for Uncertain Duffing-Holmes Chaotic Systems, Engineering Letters, Vol. 16, No. 3, pp. 316-320.
    10. Lee, S. D. , B. D. H. Phuc, X. Xu, and S. S. You (2020), Roll Suppression of Marine Vessels using Adaptive Super-twisting Sliding Mode Control Synthesis, Ocean Engineering, Vol. 195, 106724.
    11. Lee, S. D. and S. S. You (2018a), Dynamical Analysis of the Moored Vessel System under Surge Excitations, Journal of the Korean Society of Marine Environment & Safety, Vol. 24, No. 2. pp. 140-145.
    12. Lee, S. D. and S. S. You (2018b), Dynamical Rolling Analysis of a Vessel in regular beam seas, Journal of the Korean Society of Marine Environment & Safety, Vol. 24, No. 3. pp. 325-331.
    13. Lee, S. D. and S. S. You (2018c), Sliding Mode Control with Super-twisting Algorithm for Surge Oscillation of Mooring Vessel System, Journal of the Korean Society of Marine Environment & Safety, Vol. 24, No. 7. pp. 953-959.
    14. Lee, S. D. , X. Xu, H. S. Kim, and S. S. You (2019), Adaptive Sliding Mode control Synthesis of Maritime Autonomous Surface ship, Journal of the Korean Society of Marine Environment & Safety, Vol. 25, No. 3. pp. 306-312.
    15. Loria, A. , E. Panteley, and H. Nijmeijer (1998), Control of the Chaotic Duffing Equation with Uncertainty in All Parameters, IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications, Vol. 45, No. 12, pp. 1252-1255.
    16. Mitra, R. K. , A. K. Banik, and S. Chatterjee (2017), State Feedback Control of Surge Oscillations of Two-Point Mooring System, Journal of Sound and Vibration, Vol. 386, pp. 1-20.
    17. Perez, P. (2005), Ship Motion Control: Course Keeping and Roll Stabilization using Rudder and Fins, Springer-Verlag London Limited, pp. 27-44.
    18. Shah, A. A. , A. Umar, and N. A. Siddiqui (2005), A Methodology for Assessing the Reliability of Taut and Slack Mooring Systems against Instability, Ocean Engineering, Vol. 32, pp. 1216-1234.
    19. Shtessel, Y. B. , M. Taleb, and F. Plestan (2012), A Novel Adaptive-gain Super-twisting Sliding Mode Controller: Methodology and Application, Automatica, Vol. 48. pp. 759-769.
    20. Shtessel, Y. , C. Edwards, L. Fridman, and A. Levant (2014), Sliding Mode Control and Observation, Springer New York: Birkhauser, pp. 1-41.
    21. Umar, A. , T. K. Datta, and S. Ahmad (2010), Complex Dynamics of Slack Mooring System Under Wave and Wind Excitations, The Open Oceanography Journal, Vol. 4, pp. 9-31.
    22. Xu, X. , S. D. Lee, H. S. Kim, and S. S. You (2020), Management and optimization of chaotic supply chain system using adaptive sliding mode control algorithm, International Journal of Production Research,