Using Piezoelectric Materials in a Portal Frame Structure

Buoy systems are an alternative for micropowering small devices in remote locations. Portal frames are very useful to harvest the energy of the waves into usable energy. *us, using the current models for a portal frame in the literature and the spectrum of available energy in sea waves, a nonlinear mathematical model accounting for the coupling of a nonlinear piezoelectric material is considered. *e neighbour of selected variables is analyzed and then optimized by a process utilizing the particle swarm optimization (PSO) algorithm. Furthermore, an optimal control using the linear-quadratic regulator (LQR) controller is applied to control the load resistance of the piezoelectric circuit.*e optimization process and the LQR show to be effective.*e results show a general gain due to optimization and a relatively small gain using the controller.


Introduction
Due to the large coverage area of the ocean on the Earth's surface, offshore energy sources have been very attractive, bringing the attention of scientists and engineers of the energy field. is enormous interest stands out for the power that can be used as a source for consumption in different ways [1][2][3].
is energy is determined by the mechanical waves formed by the ocean's water motion. According to [1,2], it is estimated that approximately 850 TWh of electricity is generated with the conversion of this movement. In this way, several devices work with the amplitude of the tides, such as those that lift the floats that work with the sea waves generated by the wind, devices that work with the torsion system, and the ones that work by creating air pressure and thus generating electricity [4]. e authors in [5,6] utilized an air turbine actuated by the displacement through the oscillating of the ocean waves for energy generation, where such energy is stored in the form of compressed air. e authors in [7] converted the ocean wave energy through oscillating bodies coupled to a linear electrical generator. On the same hand, the authors in [8] considered the overtopping floating structures as the Wave Dragon energy converter. erefore, with the high demand for energy consumption nowadays, the need for new power sources has been growing in the last decades, especially in remote locations that make the energy conversion on site a more attractive source, and for lowpower devices, such as underwater sensors. e use of solar panels is an alternative power source that has been used for enhancing offshore energy harvesting [9]. However, it has a discontinuous load capacity, and the panels are not recyclable and required frequent cleaning.
In an attempt to improve the offshore energy conversion, the authors in [10] presented an ocean wave energy harvesting system accounting for a direct current (DC) power generator attached at the middle top of a floating platform along with a pendulum connected to the generator's shaft. It was considered that the ocean wave motion swings the platform in the vertical direction, inducing the pendulum rotation. Consequently, electrical energy is provided due to the DC power generator. For optimizing the harvested power, the pendulum parameters were analyzed by using the PSO algorithm taking into account the mass and length variation of the pendulum. In addition, the ocean wave amplitude and frequency are also considered. e numerical results showed an efficient energy conversion with the use of the DC power supply.
Piezoelectric materials are also an alternative for enhancing offshore energy harvesting. A floating platform is generally induced to vibrate due to the ocean waves, which makes the piezoelectric material an intelligent strategy. ese materials have been studied and used as shown in [11]. ese authors have explored piezoelectric materials for energy harvesting due to the wide range of frequency utility (including high frequency); they are of low-power consumption, ease of application, and very low cost.
In [12], a piezoelectric energy harvester based on a polygon-shaped cantilever array is considered. e authors employed a multifrequency operating principle for eight cantilevers with the irregular design of the cross-sectional area. A novel V-shaped vibration energy harvester based on the conventional piezoelectric bimorph cantilevered structure is explored in [13]. Experimental results showed that the V-shaped energy harvester improved the frequency response characteristic and the output performance of the electrical energy. In [14], three different kinds of energy harvesting generators are designed: one considering the monostable piezoelectric cantilever beam structure and the other two considering bistable piezoelectric cantilever beam structures. Numerical and experimental analyses were considered for power generation and dynamic behaviors of the different structures. In [15,16], the introduction of a nonlinear piezoelectric coupling in a portal frame accounting for the lateral motion of the structure is considered. Numerical results showed that energy harvesting can be influenced by nonlinear parameters of the piezoelectric model. Furthermore, the installation of the piezoelectric sensor to the seaward position in an existing coastal structure for generating energy from sea waves is carried out in [17]. e authors used a wave breaker along with the structure. Experimental results showed correlations of generation volume, and wave conditions have been found.
Another form of application is utilizing the piezoelectric material as an actuator in flexible elements to activate or prevent deflections on it. In this case, the reverse effect is a desired effect and could be used with the genetic algorithm (GA), further improving the location and resizing the piezoelectric [18]. With that, the State-Dependent Riccati Equation (SDRE) controller is used for driving the piezoelectric material to a more efficient behavior for energy harvesting [19], which may pump vibration energy into the system to control its movement.
In [9], a combination of a fuzzy-PID controller and an LQR controller was used for determining the movement of a plate. e authors concluded that both strategies are valid and robust to control the experiment and better evaluate the performance of the structure. In addition, details about the difficulty and complexity of using piezoelectric actuators in structures were discussed.
In the offshore field, as the buoy is in constant motion, it is a logical solution to harvest this vibration energy. Due to the power output, the combination of various piezoelectric materials further increases reliability. Piezoelectric material is reliable and is not only used for collecting energy. is material is present in humidifiers, printers, and high-performance diesel engines, proving its applicability in several fields and a consolidating technology. erefore, this work proposes the analysis of the energy harvesting of an offshore system containing a portal frame structure along with a piezoelectric material for energy conversion. e mathematical modelling of the system is carried out. Particle swarm optimization (PSO) is considered to optimize the energy production of the system. In addition, the LQR control technique is used for maintaining the energy conversion at its peak production. e parametrical analysis of the system is carried out throughout numerical simulations accounting for the Lyapunov exponent study, which is used for tracking the dynamic behavior of the structure to maximize the energy conversion. e rest of the work is organized in the following way. Section 2 describes the mathematical model and the structure in question. Section 3 describes the PSO algorithm. Section 4 presents the LQR optimal control. Section 5 shows the numerical simulations of the optimization of the parameters using the PSO, the project of the control, and dynamical analysis. Finally, Section 6 presents the conclusions of this study.

The Buoy System
Buoy systems are floating devices attached to the seabed that, by themselves, can be used for generating electrical energy. Most of the buoy systems are used for large macroenergy generation. However, in this case, we are interested in micropower harvesting.
In Figures 1(a)-1(c), a simple bounce system, a portal frame structure with a tuned-mass-damper (TMD) system on top of it, and the lumped parameter scheme of the system are shown, respectively [20,21]. Note that the portal frame system (Figures 1(b) and 1(c)) is located on A 1 .

Mathematical Model.
e mathematical model of the portal frame represented in Figure 1(c) is a complex iteration between the masses m 1 and m 2 . e mode is represented as a nonlinear dynamic model modelled through Newton's second law, along with the differential equation that describes the voltage of the piezoelectric material depending on the deflection of the structure. e equations of motion of the system are given by where e nonlinear piezoelectric coupling is defined as follows: where σ 1 and σ 2 are constants. As the external excitation is provided by the ocean waves, base excitation is considered as of harmonic kind which is given by [15,16,22] where A mp is the amplitude of the ocean wave and f 1 is the main frequency of excitation. Hence, equations (2) and (3) are substituted into (1) and the new dimensionless parameters as considered as follows: Equations (1) and (3) can be represented in the following form: Transforming equation (5) into state-space notation results in e harvested power is calculated according to the following equation: where τ f is the time of the simulations.

Particle Swarm Optimization
e optimization of the energy harvesting is carried out by using particle swarm optimization (PSO). e PSO mimics the behavior of groups. us, the communication between each particle provides the survival of the crowdedness.

Shock and Vibration
Mathematically, each particle uses 3 factors, which are the local best (best for each group), the global best, and the inertia of each particle. is allows the introduction of a lag in each particle to avoid high changes in the velocity of the particles. As an iterative method, the displacement of the system can be calculated as [23][24][25] where c 0 is the inertia of the system, V n is the speed of the particle at the current state, c 1 is the self-knowledge, r 1 and r 2 are random values, X lb is the direction for the local best, X n is the current position, c 2 is the group knowledge, X gb is the direction for the group best, and i is the iteration of the system.

Optimal Control
For the controller, the linear-quadratic regular (LQR) controller is chosen because it is multiple-input multipleoutput (MIMO) control that allocates poles optimized for the given restriction impose on the states and signal control [26]. e first step is to write the controlled system into statespace notation, in matrix form, as where x is the state vector, A is the state matrix for the linear, B is the input matrix, u is the control signal vector, and A nl is the nonlinear state matrix. e control u law is given by where e � ( T and x s represents the desired states for each s state.
Being P a symmetric matrix, the algebraic Riccati equation is developed, which is denoted by e control u is optimal and transfers the nonlinear system of equation (7) from any initial state to a final state e(∞) � 0. Minimizing the cost functional, it yields where Q and R are positive definite matrices.

Numerical Simulation
e numerical simulations are carried out by integrating the system of equation (7) with the fourth-order Runge-Kutta method.
e parameters accounted for in the numerical simulations are the ones described in Table 1.
In this way, we investigate the influence of some variables for the energy harvesting from the buoy system with the portal frame as the vibrating structure. In addition, the behavior of the neighborhood of these parameters is investigated and determined.
Parametric analyses for the dynamic behavior of the system and energy harvesting are carried out by varying specific nondimensional parameters of the system, such as damping of the TMD α 2 , the amplitude of the ocean waves α 4 , the linear stiffness of the TMD β 2 , and the nonlinear stiffness of the portal frame structure β 3 . e choice of these parameters is because they have a direct and strong influence on the dynamics of the system. It is expected that the TMD acts as an energy pump device, passively boosting extra vibration to the system. e parameter α 2 directly changes the amplitude of vibration of the TMD, as well as β 2 . However, β 2 tends to shift the natural frequency of the system and also may resonate with the portal frame movement, increasing the amplitude of vibration. en, due to large deformations, nonlinear behavior can be observed as well as chaos at a low frequency depending on β 3 . Offshore platforms are subjected to high amplitudes of vibration induced by the amplitude of the ocean waves (related to α 4 ), which depends on its location in the ocean.
After the initial numerical simulations, the optimization process is carried out accounting for the optimization of the parameter vector O � [β 2 , β 3 , α 2 , α 4 ] restricted to the minimization of the function cost J at 10000/P 1000 subjected to an interval of interest of the following parameters: where P 1000 is the power converter over a 1000 s. e chosen range is based on the mass of the TMD m 2 , which affects α 2 and β 2 , the nonlinear stiffness k nl for β 3 , and the amplitude of the ocean waves (A mp ) based on data from NOAA [27]. Furthermore, numerical simulations with the optimized parameters are depicted and discussed.

Dynamics and Harvested Power
Analysis. Using the parameters in Table 1, the behavior of the buoy system is analyzed through the maximum Lyapunov exponent of the system, as illustrated in Figures 2(a) and 2(c). In addition, the average harvested power is calculated for the same region of parameters, as depicted in Figures 2(b) and 2(d). Figure 2(a) presents the maximum Lyapunov exponent related to the amplitude of the wave exerted on the base of the system versus the damping coefficient of the TMD. It is observed that the higher is the damping, the system presents periodic behavior (gray area where λ < −0.005), except when the ocean wave amplitude is small (blue area where 0 ≤ λ ≤ −0.005). However, when damping is small (α 2 ≤ 0.5), irregular motions are noted for small and higher excitations. e values within the interval 0 ≤ λ ≤ −0.005 represent an interval of transition between regular and irregular motions, which may represent from periodic, quasiperiodic, to chaos. Figure 2(b) depicts the harvested power for the same interval of the parameters of Figure 2(a). It is observed that the harvested power is dominated by the excitation energy, which increases the harvested power from almost zero to P max > 500. It is also highlighted that, although irregular motions can occur with high amplitudes of motion, for this current buoy system, it is of small amplitude (green diamond in Figures 2(a) and 2(b)). Consequently, the harvested power is very small, except for the system at low TMD damping and under high excitation. Figure 2(c)) shows the maximum Lyapunov exponent for the parameters related to the linear stiffness of the TMD β 2 versus the nonlinear stiffness of the portal frame β 3 . In this case, for the same range of parameters of Figure 2(a), there is no apparent irregular behavior. However, relating Figure 2(c) to its respective harvested power of Figure 2(d), it is noted that the harvested power increases with β 2 along with the increase of the maximum Lyapunov exponent. Although the Lyapunov exponent is less than zero, it seems the system is increasing the amplitude of displacement. e variation of the harvested power is very small within the interval of obtained power as 0 ≤ P avg ≤ −0.005.

Power Optimization.
Establishing the most important parameters that the highest amount of energy is obtained, the optimization process using the PSO algorithm is carried out for the range of the parameters depicted in equation (14). Figure 3(a) shows the maximum dimensionless harvested power obtained through the PSO algorithm for various swarm sizes. e number of swarm sizes indicates the number of particles used to find the set result. It is observed that the highest harvested power is found using only a swarm size of 24, as highlighted by the blue bar. e optimized parameters found with the PSO algorithm for the highest harvested power are β 2 � 0.1, β 3 � 0.1, α 2 � 0.01, and α 4 � 1.5. Note that this result is the best among others using higher swarm sizes.
is indicates that the highest harvested power is obtained for these found parameters and there is no other possible higher value in their  neighborhood. However, the difference among all the obtained harvested power is almost insignificant. Figure 3(b) shows the number of iterations for each swarm size. Note that the optimized result needed more iterations than the other results. However, it is expected due to the small number of particles of the PSO. Figure 4 presents the number of iterations to minimize the function described by equation (14). During the PSO optimization process, all the solutions converged to the same value. It shows that, for this case, the size of the cluster does not affect the result. For a convergence with smaller cluster size, the convergence of the parameters has a satisfactory result.
On this hand, the dynamic behavior of the system and the harvested power are analyzed and calculated using equation (8) around the neighborhood of the optimized parameters, as depicted in Figure 5. Figure 5(a) shows the map of the maximum Lyapunov exponent of the system for the variation of the dimensionless amplitude of the ocean waves α 4 versus the damping coefficient of the TMD α 2 , accounting for the optimized β 2 and β 3 , for the same intervals as in Figure 2. After the optimization, it is noted that the system becomes periodic in almost the full region of values (gray area where λ < −0.005). However, when damping is higher, the system may present quasiperiodic behavior and chaotic behavior.
is is highlighted by blue and red areas of the inset (purple diamond) of Figure 5(a), where the calculated maximum Lyapunov exponent is λ � 0 (quasiperiodic or chaos) and λ � 0.3 (chaos). e interval of the parameters of the inset is 1.8 ≤ α 2 ≤ 2.0, and 0.001 ≤ α 4 ≤ 0.005. With that, the unstable behaviors only occur for higher damping and low excitation.   Figure 5(a). It is observed that the maximum harvested power increases from around P max � 500 to P max � 2000, where the maximum power is obtained for the set of parameters represented in the blue diamond in the figure (α 2 � 0.01 and α 4 � 1.5). Note that the same situation of the maximum harvested power, being at a periodic region and with high excitation amplitude, applies for the optimized parameters. Figure 5(c) shows the maximum Lyapunov exponent for the parameters related to the linear stiffness of the TMD β 2 versus the nonlinear stiffness of the portal frame β 3 , accounting for the optimized α 2 and α 4 . In this case, for the same range of parameters of Figure 2(c), irregular motions (blue regions where 0 ≤ λ ≤ −0.005) occur for the low linear stiffness of the TMD. It seems that the nonlinear stiffness of the portal frame becomes more dominant in the response of the system.
In addition, Figure 5(d) depicts the harvested power for the same region of Figure 5(c). A notable increase of the harvested power is observed in comparison with Figure 2(d), from P max � 6.18 to P max � 333.8 for β 2 � 0.1 and β 3 � 0.1. However, in this case, the maximum harvested power is obtained under an irregular motion that may be quasiperiodic or chaotic.
It is important to highlight that the PSO algorithm found the maximum harvested power for the abovementioned parameters with a very low computation task. Although it does not account for the analysis of the behavior of the system, it shows a neighborhood to be studied where it would take a long time to be optimized using other direct parametrical analyses.

Optimal Control Results.
Considering equation (7) with the addition of a control signal (u), according to equation (10), the following matrices A and B are obtained: Note from equation (15) that the matrix A is a function of x 1 . Considering the particular case where x 1 � x 1 0 � 0, the matrix A is obtained as Next, the matrices Q and R are defined as follows: Q � Replacing the matrices of equations (16) and (17) into equation (12), the vector gain k of the control law of equation (11)  e purpose of introducing a semiactive control by the LQR control is to enhance energy harvesting. In this way, to extract the maximum amount of energy, the desired variable x 5 is set as To implement the control signal as an actuator, the introduction of a variable resistor in the electrical circuit of the piezoelectric system is considered, according to the following relationship: where R 2 is obtained from where R is the general resistance, R 1 is the resistance of the circuit, and R 2 is a load resistance, which may be a battery or circuit to charge and exploit this energy.
To further increase the energy harvesting, the control is applied to the parameters obtained through the PSO. Figure 6 shows the time histories of the TMD displacement x 1 (see Figure 6(a)), the portal frame displacement x 3 (see Figure 6(b)), and the voltage x 5 obtained through the piezoelectric material (see Figure 6(c)), accounting for the cases before the optimization, after the optimization, and with the control. For both cases with parameters obtained using the PSO and the system with control, an improvement of the voltage is obtained. It is also important to highlight that the control reduced the displacements of x 1 and x 3 and increased the variation of the   voltage x 5 , which are the desired states defined previously as x 5 � x 5 and x 1 ≈ x 2 ≈ x 3 ≈ x 4 ≈ 0, demonstrating that the control proposal is very efficient. In addition, when the system is set under the optimized parameters, it is observed the increase of the displacement of x 1 and x 3 , which may be desired in many cases to increase the energy harvesting. However, it also may compromise the integrity of the structure due to large vibrations. Figures 7(a) and 7(b) show the frequency response and the energy harvesting of the system for the three cases varying the dimensionless excitation frequency (α 5 ) around the resonance, respectively. ere is a notable difference between the harvested energy of the optimized system without and with control to the nonoptimized one, as illustrated in Figure 7(b). Note that energy harvesting increases with the increase of frequency. Due to the parametric excitation, the higher is the amplitude of excitation (α 4 ), the natural frequency shifts, and hardening nonlinear behavior is noted as the cubic nonlinearities dominate the response (see Figure 7(a) for the optimized parameters). When the control is applied, the harvested power is increased with the increase of the frequency. Although the control does not change the motion of the structure, it changes the resistance of the piezoelectric circuit, as the control is chosen to actuate it, varying the value of ρ, which is directly proportional to the harvested power.
It is also important to highlight that as the ocean waves are of low frequency due to their nature, the structure oscillates out of resonance. However, for a low frequency at f � 1 Hz (α 5 � 0.0979), the optimized values and the control still make the system harvest a large amount of energy. Figure 8 shows the variation of ρ due to the resistance R 2 variation, according to the proposed control law (equation (18)). e proposed control strategy demonstrates to be feasible for applications in a real system, and that the increase of ρ implies that the proposed control is consistent with equation (19).

Conclusions
is work proposed a model for an offshore energy harvesting application excited by ocean waves. e PSO technique was used to analyze and optimize the parameters that maximize the energy production with a piezoelectric material for different conditions. e technique showed to be very efficient to optimize the energy harvesting obtained with the parameters O best � [β 2 , β 3 , α 2 , α 4 ] � [0.1, 0.1, 0.01, 1.5]. In addition, the results were obtained below the size of the cluster. erefore, we analyzed the behavior of the system in the neighborhood of these parameters with the maximum Lyapunov exponent, establishing the existence of chaotic behavior in it. However, chaos is found for the very low region with a low Lyapunov exponent, as seen in Figures 2  and 4. In addition, chaotic behavior was found close to the optimized values (see Figure 4(c)).
A semiactive control applied by the LQR controller proved to be very effective in controlling the system and increasing the energy harvesting.

Data Availability
All data and information are available within the article. Additional data or information can be requested from the corresponding author.