Digital Commons @ Michigan Tech Digital Commons @ Michigan Tech Sliding Mode Control of a Nonlinear Wave Energy Converter Sliding Mode Control of a Nonlinear Wave Energy Converter Model Model

The most accurate wave energy converter models for heaving point absorbers include nonlinearities, which increase as resonance is achieved to maximize energy capture. The efficiency of wave energy converters can be enhanced by employing a control scheme that accounts for these nonlinearities. This project proposes a sliding mode control for a heaving point absorber that includes the nonlinear effects of the Froude-Krylov forces. The sliding mode controller tracks a reference velocity that matches the phase of the excitation force to ensure higher energy absorption. This control algorithm is tested in regular linear waves and is compared to a complex-conjugate control and a nonlinear variation of the complex-conjugate control. The results show that the sliding mode control successfully tracks the reference, keeps the device displacement bounded, and absorbs more energy than the other control strategies. Furthermore, the controller can accommodate disturbances and uncertainties in the dynamic model.


Abstract
The most accurate wave energy converter models for heaving point absorbers include nonlinearities, which increase as resonance is achieved to maximize energy capture. The efficiency of wave energy converters can be enhanced by employing a control scheme that accounts for these nonlinearities. This project proposes a sliding mode control for a heaving point absorber that includes the nonlinear effects of the Froude-Krylov forces.
The sliding mode controller tracks a reference velocity that matches the phase of the excitation force to ensure higher energy absorption. This control algorithm is tested in regular linear waves and is compared to a complex-conjugate control and a nonlinear variation of the complex-conjugate control. The results show that the sliding mode control successfully tracks the reference, keeps the device displacement bounded, and absorbs more energy than the other control strategies. Furthermore, the controller can accommodate disturbances and uncertainties in the dynamic model.

Introduction
Recently, the generation of electricity from ocean waves has gained special attention.
Studies of the potential global market for wave power showed that the world's wave power resource is estimated to be 2 TW [1]. Additionally, waves have a very high power density, requiring smaller devices to capture the energy carried by the incoming waves.
Wave energy can also improve energy security by complementing the output of other

Literature Review
Researchers have used well-established control methods in WECs that include Complex-Conjugate Control (CCC), Latching, and Model Predictive Control (MPC) [5]. In complex-conjugate control, optimal energy absorption is sought by tuning the power take-off (PTO) resistance and reactance to cancel the system's inherent hydrodynamic resistance and reactance [6]. This control strategy is conceptually simple and presents a low computational cost; however, it leads to excessive motions and loads of the WEC. This can be avoided by implementing alternative sub-optimal control schemes that include physical constraints on the motions and power of the device. Latching control is based on locking the device through dedicated mechanisms during a certain period of time of the wave oscillation cycle to achieve resonance of the WEC [7], [8]. Resonance is achieved when the system is being excited at its natural frequency. While latching control solves the PTO limitation problem of CCC, it is prone to failure of the mechanical clamping pieces generating extra costs and reducing its reliability. MPC strategies maximize the energy absorbed by applying the optimal control force to achieve resonance at each time step over a future time horizon [9]. However, a prediction of the wave motion is required, and its complexity of implementation in real-time demands high computational requirements.
Most of these control techniques use the linear model developed by Cummins [10] with hydrodynamic parameters obtained from boundary element methods (BEM) due to their simplicity and high computational efficiency. Although linear models are effective for most ocean engineering applications where the aim is typically to stabilize the structure, linear models are inadequate for WEC control, where the primary objective is to amplify the device's motions to maximize power absorption [11]. As a result, control strategies based on linear models are often suboptimal and lead to higher LCoE.
Nonlinear models often include the quadratic term of the pressure from the incident flow, the integration of the pressure forces over the instantaneous wetted surface of the device called Froude-Krylov (FK) forces, and the nonlinear incident flow potential [12]. Each of these nonlinearities adds an extra level of complexity and computational load. However, it is understood that nonlinear models are more accurate than linear models, particularly for buoy geometries with varying cross-sectional areas. A previous study [13] that compared different modeling options showed that the dominant source of nonlinearities is likely to be the nonlinear Froude Krylov forces for devices with a varying cross-sectional area.
While in linear modeling, the FK forces are calculated over a constant wetted area of the device, in nonlinear models, the pressure is integrated over the instantaneous wetted surface area. This is generally done by implementing a very fine mesh or a remeshing routine of the surface, which requires more computational effort. As a result, the nonlinear model is more accurate but is very computationally expensive. However, in [14], a computationally efficient analytical method was developed to compute nonlinear dynamic and static Froude-Krylov forces based on the instantaneous wetted surface area of the buoy. This analytical method presents only a 2% error compared to the more computationally expensive methods. Therefore, in this study, the nonlinear FK forces will be incorporated using a variation of the algebraic solution proposed by Giorgi et al.

Project Overview
In this research project, a Sliding Mode Control (SMC) strategy is proposed that incorporates the nonlinear Froude-Krylov force reducing the discrepancies between the mathematical model and the actual system. In addition, this SMC is designed to track the desired reference signals using low computational efforts and rejecting disturbances encountered by the WEC or parametric errors presented in the model. The controller will be investigated in a simulation environment using a standard spherical WEC and regular waves.
A sliding mode control strategy for heaving WECs was proposed in [15], where it was shown that this control strategy provides a viable solution to increase energy extraction.
The SMC was designed to track a reference velocity without violating the physical constraints of the device. The results showed that the SMC successfully tracked the reference signal even when system perturbations were present. However, this work was done using a linear approximation of the hydrodynamics of the system.

Project Objective
This project aims to improve the performance of WECs using a nonlinear model that incorporates the nonlinear FK force and SMC to drive the system to the desired velocity profile while keeping the system very robust against parameter perturbations of the model and external disturbances due to different environmental conditions.
The main contributions of this study consist of: where the SMC strategy is designed and compared to a complex-conjugate control and a nonlinear variant of complex-conjugate control. In Section 4, the simulation used to evaluate closed-loop performance is described. Section 5 shows the results, and Section 6 provides conclusions and areas for further study.

Model Description
Mathematical models for WECs are essential for the design of model-based control.
These models are typically linear because of their low computational requirements, in which they assume small motion. However, there are situations in which the movement is not small while still producing power, and the nonlinearities of the device become significant.
The incoming waves are classified as linear and nonlinear depending on the wave steepness , which is the ratio between the wave height ℎ, and the wavelength [16].
Waves are also classified as regular and irregular, where regular waves are characterized by a specific amplitude , frequency 0 , and wavelength, whereas irregular waves are composed of multiple frequency waves. In this report, linear and regular waves are used, representing the wave free surface elevation ( ) as a function of time only: A floating single degree of freedom (SDOF) buoy is shown in Figure 1, where the lower part of the body is attached to a linear generator for power take-off (PTO), which is fixed to the seabed. The PTO subsystem converts the wave energy into electrical energy and energizes the system to track the reference profile. The excitation force is shown as ( ), the buoy's displacement from its equilibrium position is ( ), and the control force implemented by the PTO system is ( ).
The general heaving WEC hydrodynamic model used in this project is described as: where, :

Diffraction Force
When the approaching wave field encounters the buoy body, the waves deform around the corners of the body. This phenomenon is referred to as the diffraction force, which is the force that the incoming waves apply on the buoy as it is held fixed in the water [6].
On the wave energy field, it is usually represented as the convolution of the product of the wave free surface elevation and the diffraction impulse response function (IRF): where is the diffraction IRF.
When the buoy diameter is much smaller than the wavelength, the disturbance field generated by the body is small enough to be ignored [17]. Therefore, the diffraction force will be included in the equations in this report but not in the simulation. This force could also be accounted as a disturbance to the system.

Radiation Force
The radiation force is the force that the body experiences when it oscillates in calm water.
It is common to divide this force into two parts: the added mass force proportional to the body's acceleration and the wave damping force proportional to the body velocity [18].
Thus, the radiation force can be represented as: where, : Radiation IRF Evaluating the convolution integral in the radiation force is very expensive to compute at every time step of the simulation, and it is not possible on a real-time application.
However, when the buoy is exposed to only harmonic waves of frequency 0 , the convolution in (4) can be represented as: where ( 0 ) is the linear damping constant.

Froude-Krylov Force
The Froude-Krylov force is the addition of the hydrostatic force with the hydrodynamic force . The hydrostatic force is the difference between the gravity force and the force caused by the hydrostatic pressure over the wetted surface of the buoy and, the hydrodynamic force is the integral of the unsteady pressure field over the wetted surface of the floating buoy caused by the incident waves.
In linear models, the integration is done over the constant mean wetted surface, while in the nonlinear approach, the instantaneous wetted surface of the body is considered. This approach typically requires significant computational effort. However, Giorgi et al.
proposed an algebraic solution that significantly reduces the computational effort for heaving axisymmetric geometries. The algebraic solution is found by defining the pressure ( ), the infinitesimal element of the surface , and the limits of integration.
The total pressure for deep water waves can be found using Airy's wave theory. where, where ( ) is the profile of revolution of the point absorber. Using the radial and tangent vectors canonical basis, and only the vertical component for heave motion restriction, the infinitesimal surface element, and the pressure becomes: Combining these equations, the Froude-Krylov forces are: The paper continues defining the limits of integration based on the free surface elevation and the draft of the buoy ℎ 0 at equilibrium as: For a spherical point absorber, the profile of revolution derived in Giorgi et al. paper is: 2 where is the radius of the buoy, and ( ) is the vertical displacement of the buoy from equilibrium.
In this project, the limits of integration and the profile of revolution have been redefined using calculus to derive the instantaneous wetted surface area of the sphere. Considering the bottom of the sphere as the lower limit of integration, the profile of revolution of the sphere at still water level shown in Figure 2 is: ( ) = √ 2 − (− + ℎ 0 + ) 2 where the SWL was used as a reference. In this case, the lower limit of integration is −ℎ 0 and the upper limit is 0. When the buoy is displaced by ( ), and the free surface elevation ( ) is not zero as shown in Figure 2. b the profile of revolution of the sphere becomes: The limits of integration are: 0 for the upper limit, and = −ℎ 0 + ( ) − ( ) for the lower limit. Combining equations (8)- (12), the integral of the Froude-Krylov force of a heaving spherical buoy can be expressed as: Solving the integral and separating the Froude-Krylov force, the resulting nonlinear hydrostatic and hydrodynamic forces of the system are: The complete nonlinear model of a single degree of freedom WEC in heave mode can then be described as: where = ( + ( ∞ )) For the subsequent control design, it will be convenient to write (16) as: where is the true model of the system and ̅ = −1 where is the true mass. In contrast, the model used for control design is approximated as: The calculation of the absorbed energy is [19]:

Control Design
This section presents a sliding mode control (SMC) design, illustrated in Figure 3, for a heaving WEC based on the nonlinear differential equation model of Equation (16).
Motivated by linear, complex conjugate control, the SMC reference trajectory is selected to match the frequency and phase of the excitation force. Since the WEC model is not linear, as assumed by complex conjugate control, it is recognized that this reference trajectory will not be optimal. Creating a real-time, energy extraction optimal reference trajectory is an open topic and well-suited for MPC applications. The reference velocity is then integrated to produce the reference displacement profile ( ). The sliding mode control is derived by first defining the sliding surface as: where is the state vector = [ ;]. The problem of tracking a given profile now becomes equivalent to maintaining the stationarity condition: Substituting (18) into (21), and solving for ̂ gives the control law:

̂= ( ( ),( ), ( )) − (̈( ) +( )) (22)
This control law will only work if the assumed form of the model is perfect and there are no disturbances. To account for these inevitable situations, Lyapunov's Direct method is used to augment the control law of (22) to ensure stability.
The Lyapunov candidate function is selected as: Furthermore, note that ̂≠ ̅ are related by ̅ =̂. The time-dependency notation will be omitted for the following equations for brevity. We now have the stability condition: As long as is sufficiently large to dominate the uncertainties on the right side of (28), the system will be stable.
The discontinuous terms in the control of (26) ( ( )) can cause chattering. This is eliminated by replacing the signum function with a hyperbolic tangent [20]. Then, the control law is interpolated by replacing ℎ( ( )) function by ( )/ as: where Φ is the boundary layer thickness.
The reference signal generator outputs the desired velocity profile of the WEC plant that is then integrated to obtain the desired displacement profile . The control law takes as inputs the formed error signals, the desired acceleration profile ̈, the FK forces and the actual buoy displacement as in (29). The control parameters and Φ were chosen during the simulation process. As aforementioned, the reference velocity profile is designed to match the frequency and phase of the excitation force. As described in [6], the optimum velocity condition requires: To achieve this condition, phase control is needed so that the oscillation velocity (t) is in phase with the excitation force and amplitude control to get a velocity amplitude of |̇( )| = | |/2 ( 0 ). However, this optimal condition leads to excessive motions of the device that would be infeasible in real-world applications. Therefore, the reference velocity profile is defined to be in phase with the excitation force with a suboptimal amplitude that varies for each sea state condition. Even though the WEC plant is nonlinear, in this study, the calculation of the reference velocity profile is based on the linear excitation force obtained from WAMIT [21].
To compare the results of the SMC, an approximate CCC [22] is used. Referring to the impedance matching principle [23], the optimal control force may be written as: where the optimal control impedance ( ) must equal the complex-conjugate of the intrinsic impedance = ( ) + ( ) composed by the intrinsic resistance ( ) and the intrinsic reactance ( ). For monochromatic incident waves, this optimal force may be written in the time-domain as: For optimal control, the controller values must equal the intrinsic values = , ( 0 ) = ( 0 ), = where is the linear approximation of the static Froude-Krylov force. By cancellation, this control force gives the optimal response velocity shown in (30) when the WEC model is assumed to be linear, and is used instead of the nonlinear hydrostatic FK force. Since in this paper the WEC model is nonlinear, the optimal control force using CCC was modified to obtain a perfect cancellation of the as: This control strategy used for comparison will be referred to as nonlinear complexconjugate control ( − ) from this point on. The configuration of this control scheme is similar to that of the SMC shown in Figure 3. However, this control law only takes as inputs ( ),( ),( ) and, ( ).

Simulation
As shown in previous studies [12], the significance of the Froude-Krylov nonlinear force for a spherical body is great due to its varying cross-sectional area. The buoy parameters used in this simulation are listed in Table 1 and were chosen to resemble the parameters of a real device such as the WAVESTAR [24] device. The density of the body was chosen to be 500 kg/m 3 to make the draft of the buoy ℎ 0 coincident with the center of gravity of the sphere. The remaining hydrodynamic parameters were found by solving the radiation problem in WAMIT.
The results shown in the following section correspond to the wave period of = 6s and 0 = 1.05rad/s. Therefore, the frequency-dependent parameters shown in Table 1 are based on this wave frequency.
The simulation was run for different monochromatic wave periods ranging from 3 to 9 seconds. Under the deepwater assumption, the wavelength is related to the wave period as: In addition, the steepness of the wave was kept constant at 0.018, which is the highest wave steepness under the linear waves assumption [16]. To maintain the wave steepness constant, the amplitude of the wave is related to the wavelength as follows: Therefore, the wavelength and wave amplitude are determined by the choice of wave period. Table 2 shows the range of wave periods and their corresponding wavelength and wave amplitude. The reference velocity profile in (30) has been proven to be the optimum velocity [6] when using a linear approximation of the WEC, enabling maximum energy absorption for a given buoy shape. However, this is not the case for this nonlinear WEC model where the is dependent on the incident wave elevation and the buoy displacement.
Additionally, it is well known that CCC magnifies the buoy's motion, driving the device to be either fully submerged or completely out of the water for certain wave periods. This undesirable behavior can be eliminated with the SMC by keeping the amplitude of the desired displacement bounded: where is a safety factor used to determine the maximum and minimum buoy height covered by water, as shown in Figure 4. Therefore, the amplitude of the desired displacement of the device was chosen to be a value within those limits. The SMC reference amplitude for each wave period is listed in Table 2.

Results
As aforementioned in Section 3, this report compares the performance of the SMC against a complex-conjugate control (CCC) and a complex-conjugate control that includes the nonlinear static FK force in the control law (NL-CCC). The reference velocity used for the SMC is a cosine wave of the same frequency and phase as the linear hydrodynamic Froude-Krylov force obtained from WAMIT. To reduce the risk of driving the buoy completely out of the water or completely submerged into the water, a safety factor of = 10% was used in the simulation, leaving 25cm of clearance at the top and bottom of the spherical body.
The case shown in the results corresponds to a case where the CCC leads to a stable steady-state solution within the physical limits of the buoy. For higher wave periods, the solution with CCC would not reach steady-state. The maximum wave period at which the WEC controlled with the complex-conjugate control reached a steady-state solution within the physical limits of the buoy was at = 6 . Therefore, that is the case that is shown in this report. Figure 5 shows the relative displacement ( ) − ( ) of the system using the different control strategies and the displacement limits to keep the buoy safe.
The relative displacement is kept within limits for the three control strategies used. The relative displacement when SMC is implemented is very well-behaved, and it takes advantage of the full range of motion available for the WEC. When CCC is implemented, there is a transient part for the first 15 seconds due to a 20 seconds ramp used to attenuate the overshoot caused by the control force. The relative displacement then reaches steady state within the allowable range of motion. Note that the relative displacement with CCC is not in phase with the SMC. If the WEC plant were linear, the CCC and SMC would be in phase. However, since the plant is nonlinear and the controller is a linear approximation, there is no perfect cancellation of the static FK force leading to a mismatch in relative displacements. Once the nonlinear static FK force is added to the CCC to form the NL-CCC, the relative displacement is in phase with the SMC. However, the amplitude of the relative displacement does not cover the full range of motion available since the minimum point at steady-state is already very close to the lower limit. The energy absorbed by the system when using the three control strategies and the control force used is shown in Figure 6. From this figure, it can be observed that the CCC performance is very poor compared to the SMC and the NL-CCC performance. This low energy absorption is caused by the very limited control force needed to maintain the WECs motion bounded. On the contrary, the energy absorbed when using SMC and NL-CCC are similar, with the SMC being slightly higher. However, the NL-CCC required a higher magnitude control force than the SMC to achieve slightly less energy absorption.
Additionally, the SMC control force is a continuous function that does not need to be clipped in order to keep the WEC safe. On the other hand, the control force of CCC and NL-CCC are both saturated, causing discontinuities in the force. Furthermore, in order to maintain the WEC bounded in the safe range of motion, the control force of CCC and NL-CCC has to be limited at a different value for each sea state. While for the SMC, this is not necessary since the motion of the WEC is determined by the reference signals.
Although, in this case, the difference between the energy absorbed by the SMC and NL-CCC is very close, that is not the case for higher wave periods. As the wave period increases, because the wave steepness is kept constant, the wave height increases. With increasing wave amplitudes, the WECs displacement is significantly different from the wave elevation. Therefore, the mean wetted surface is significantly different than the instantaneous wetted surface used in this nonlinear approach, making the nonlinear Froude-Krylov forces very significant. Hence, when using a linear control strategy such as the CCC for higher wave periods, the control force must be highly limited in order to keep the WEC safe, reducing the power absorption significantly. Similarly, when using the NL-CCC, even though the nonlinear static FK force is incorporated in the controller, the nonlinear dynamic FK becomes very relevant; therefore, the control force must be limited to keep the WEC bounded. The power absorbed by the WEC when using the three control strategies is shown in Figure 6. periods [4].
The exact model assumes 100% certainty of the three parameters mentioned above, and the perturbed model assumes 50% uncertainty of the three parameters simultaneously. In Figure 8 is shown a line with slope − in blue, and the sliding surface ( ) in red for the perturbed model. In Figure 8.a the SMC coefficient of the term tanh ( / ) in (29) was kept the same as the value used in the exact model 10 4 N.
However, as seen from the figure, the initial loop starts at (0;0) due to matching initial conditions between the system states and the reference states, then the sliding surface keeps looping around the blue line, increasing in magnitude. This indicates that the error is not bounded, and therefore, the system could go unstable. Therefore, was increased to 10 8 N to handle the uncertainty errors. The sliding surface shown in red in Figure 8.b forms a closed contour and loops around it, meaning that the error is controlled and will not go unbounded. Hence, when the model is not exact, the needs to be increased. For different wave periods, and depending on the modeling error or disturbances, the SMC coefficient might have to be adjusted to ensure the WEC is tracking the reference signals without much deviation.  SMC coefficient had to be increased. As mentioned in section 3, the higher the magnitude of this coefficient, the more robust the SMC is. However, as seen from the control force in Figure 9.b, the increase in its magnitude barely affected the control force.
This behavior is because the term ℎ ( ( )) of the control law offsets the difference in the terms with uncertainties between the exact and perturbed model.
The energy absorption shown in Figure 9.c does not seem to be influenced by the slight error in the WEC's displacement and the control force. Both the exact and perturbed models needed almost the same control force to deliver equivalent energy absorption. Therefore, the SMC has proved to work very well in systems prone to modeling errors and environmental disturbances such as wave energy converters.

Conclusion
Three different control strategies have been presented and applied to a spherical heaving point absorber wave energy converter. These control strategies have shown promising results when used in linear models of WECs. However, it has been shown that linear model approximations lead to inaccurate models when there is a significant difference between WEC displacement and wave elevation. Therefore, this research project presented a sliding mode controller applied to a nonlinear WEC model. This plant incorporates an algebraic nonlinear static and dynamic Froude-Krylov force.
The SMC was compared to a complex-conjugate control, and a variation of this complexconjugate control with a nonlinear static FK force term. The controllers aimed to obtain the higher power possible of all control strategies tested while maintaining the WEC within a safe range of motion. The simulation was done using linear and regular waves at a constant wave steepness of 0.018.
It was found that the proposed SMC successfully tracked the reference signals given at all of the wave periods tested. This was achieved by a continuous control force without constraints in magnitude, which would reduce the chance of failure in real-world applications by making the device more reliable. Furthermore, it was shown that the control force does not need to be limited at different values for each sea-state in order to keep the WEC in the safe range of motion, as long as the reference signals fall within those safe bounds.
The sliding mode controller allowed the maximum energy absorption between the control strategies tested at all of the wave periods. At lower wave periods, this was done by requiring less control force than the NL-CCC, and at higher wave periods, where the nonlinearities become very significant, the energy absorbed when using SMC proved to be substantially greater than when using CCC and NL-CCC. In addition, the sliding mode controller allows for the use of any desired reference signal, meaning that if there already exists a path that would deliver the optimal absorbed energy, it could just be used as a reference in the SMC and achieve optimal performance. Furthermore, the SMC was simulated with modeling errors of 50% in the radiation force, and 50% error in the wavenumber used to calculate the nonlinear dynamic FK force.
Despite these uncertainties, the SMC kept the wave energy converter in the tracking reference while using a continuous control force and delivering the same energy absorbed as the exact system. This robustness characteristic of sliding mode controllers is ideal for marine devices where the system parameters may vary over time due to aging, biofouling, corrosion, etc.
This paper used purely heaving point absorbers due to the computationally efficient algebraic nonlinear Froude-Krylov forces derived from Giorgi et al.'s work. However, this could be extended to more degrees of freedom in the future, where numerical approaches need to be used to solve the Froude-Krylov integrals, as explained in [14].

Future Work
In the future, more research needs to be done in order to extend the proposed sliding mode controller to work in irregular waves to test the controller in several different wave frequencies and amplitudes, and therefore, obtain a closer simulation to real-world scenarios. This work has the potential to turn into a second publication in controls of wave energy converters. Furthermore, the tracking signals used in the sliding mode controller can be optimized using different optimization solvers and obtain the reference signals that would absorb the maximum energy.
Additionally, the nonlinear effects of the Froude-Krylov forces developed in this report and the algebraic solution presented by Giorgi. et. al. can be experimentally validated in the wave tank at Michigan Tech by using a scaled-down prototype of a point absorber WEC.