Dimensional parameters of the TLSRG.

## Abstract

This chapter addresses the dynamic simulation and control of linear switched reluctance generators for direct drive conversion systems. The electromechanical energy conversion principles of linear switched reluctance machines are briefly explained. A detailed mathematical model is developed for linear switched reluctance generators. The different types of control strategies adopted for switched reluctance generators are referred. The hysteresis controller is applied to control the conversion system with constant damping load. The proposed control strategy also includes a DC/DC isolated converter to control the system DC bus voltage by adjusting the energy flow between the conversion system and the resistive load. The mathematical model is applied to simulate the behavior of a tubular linear switched reluctance generator as power take-off system in an ocean wave point absorber device. To accomplish this task, the dynamic equations of the point absorber are presented and integrated with the linear switched reluctance generator dynamic model. In the simulation process, the system is driven by a regular ocean wave and operates with constant damping load. The system performance is evaluated for different load values, and the simulation results are presented for the optimal damping load case scenario.

### Keywords

- tubular linear generator
- switched reluctance machine
- direct drive conversion
- ocean wave energy

## 1. Introduction

The switched reluctance machine (SRM) is an electromechanical structure with salient poles defined by a fixed and a movable part as the main constitutive elements. This machine is characterized by the absence of permanent magnets and by the electric phase coils only being present in just one part (usually the fixed part). The part where the phase coils are located is designed as primary and the other part, which works as passive element, as secondary. The SRMs can be rotating or linear. Since these two configurations are homologous, they share the same operation principles with the exception of the electromagnetic force direction [1]. In SRMs, a force is developed when the magnetic structure tends to minimize the magnetic energy by displacing the movable part to achieve a configuration with minimum reluctance. This type of machine is characterized by a nonlinear electromagnetic behavior and can work as actuator or generator when the electric phase is energized during the period of increasing or decreasing inductance, respectively [2]. The SRM may have a flat structure [3] or a tubular one [4]. The tubular topology may overcome some disadvantages of the flat machine. The former has no edge effects due to the closed geometry. Also, the radial attraction forces are canceled due to its axisymmetric configuration [5].

This type of machine is suitable for low-speed applications and is a robust and low-cost solution since no permanent magnets are required for its operation. Also, the energy losses from Joule effect occurs in just one part and can, therefore, be more easily dissipated [6]. These characteristics make the SRM a strong candidate for generator in direct drive conversion systems. It has already been adopted as linear generator for ocean wave point absorbers [3]. Due to their unique advantages, the application of linear switched reluctance generators (LSRG) in direct drive conversion systems is of significant importance. Thus, design tools are necessary to develop and evaluate this type of machines, and suitable control techniques must be applied to ensure its reliable operation.

The aim of this chapter is to present the mathematical model to assist the design and assessment of LSRG. Also, control techniques are proposed for low-velocity operation. The mathematical formulation is applied to simulate the performance of a tubular linear switched generator (TLSRG) operating as a power take-off system in an ocean wave point absorber converter.

## 2. Electromagnetic conversion principles

In SRM, the magnetic circuit of each electric phase may be characterized with different values of magnetic reluctance for distinct relative positions of the movable part, in respect to the static part. When a magnetic field is established in this circuit, an electromagnetic force is developed to displace the movable part to the position with minimum magnetic reluctance in order to minimize the energy in the system. The referred position, designed as alignment position

The typical linear inductance profile for SRMs is displayed in Figure 1. The inductance is given as function of the relative position between the movable and fixed parts. As consequence of the structural pattern inherent to these machines, when the movable part is moved in a single direction, the inductance is characterized by a periodic and symmetrical profile resulting from successive movements, which are approaching (the rate of change of the inductance

By energizing the electric phases in the approximation region or separation region, the SRM may operate as an actuator or a generator, respectively. In the first case, magnetic flux established by the electric phase will tend to minimize its energy developing an electromechanical force to achieve structural alignment. As a result, a linear force will act on the machine’s movable part. In the presence of an external load, which can overcome the referred electromagnetic force, the generator is displaced from the alignment position, increasing the reluctance of the active phase, which will reduce the respective magnetic flux. As a consequence, a back-electromotive force is developed, seeking to increase the electric current to restore the magnetic flux. During this procedure, mechanical energy is extracted from the movable part of the generator and converted to electrical energy [7].

The electric phase equation is given as:

where for phase

with

These conversion processes may only occur with appropriate power electronic converter to control the energy flow. With proper switch commands, the conversion periods may be established, and electric phase current intensity may be regulated by the converter [8].

In [9, 10], several electronic converters are identified that are suitable for electric generation with SRMs. The asymmetric H-bridge converter is the typical choice. This converter topology can be found, for one phase, in Figure 2. It is a practical and simple solution, but it is characterized with a variable output voltage due to the self-excitation process. However, this drawback may be minimized with an external voltage source [11]. This converter has the less apparent power required to operate and is classified as an economical solution [10].

With the converter illustrated in Figure 2, the conversion cycle may be defined by the excitation period, the generation period, and the free-wheeling period, an intermediate stage between the first two. The circuit configurations for the different periods are illustrated in Figure 3. Assuming that the capacitor is already at its nominal operating voltage, the excitation period is initiated when the switches

When the current intensity reaches a certain value, the switches are opened and the generation period begins. During this period, the back-electromotive force increases the electric current as a consequence of the magnetic flux reduction. Thus, the electric current is maintained through the diodes

where

The electric energy returned by the phase to the converter is:

with

## 3. Dynamic model

To assess the dynamic characteristics of switched reluctance generator, a mathematical model must be formulated to describe the system behavior in transient state. The mathematical model can be solved with numerical methods and computational calculus in order to estimate the system response for the given operating conditions. The mathematical model of SRMs is obtained from Eq. (1) that describes the transient phenomena involved in the electromagnetic conversion of each electric phase. The solution for this equation may be obtained by applying time integration to the linkage magnetic flux [13] or to the electric current derivative with time [14]. For both methods, it is required to relate the machine electromagnetic characteristics with the phase electric current and the movable part position. Several approaches have been proposed in the literature to include the machine nonlinear nature in the mathematical model. In [15], look-up tables are used to model the electromagnetic characteristics of the SRM. These are obtained from a static analysis and are used as data base to build the look-up tables. Analytic functions are also used to represent the machine magnetic characteristics within the mathematical model. The use of analytical expressions simplifies the computation process related to the differentiation and integration of the electromagnetic entities [16]. This approach can be achieved by fitting an appropriate function to discrete data [17]. Fourier series expansion may also provide analytical expressions to represent nonlinear electromagnetic characteristics in the SRM mathematical model [18]. A simpler method, but less precise, comprises the use of piecewise functions to establish a linear [19] or nonlinear [20] relation with the independent variables. Artificial intelligence-based methods have already been adopted in the dynamical analysis of SRMs. These methods comprise the use of neural networks and fuzzy logic with real data to develop a mathematical representation of SRM electromagnetic characteristics [21].

All the referred approaches require the representative data of the electromagnetic characteristics, usually expressed as a function of the electric phase current intensity and movable part relative position. With look-up tables, these data are used directly. Methods based on analytic expressions or artificial intelligence models apply these data to develop appropriate mathematical expressions.

The required discrete data may be obtained through experimental measurements. This process provides realistic curves, but a physical model of the machine is needed for it. In electric machine design, where several structural possibilities must be assessed, it becomes unpractical to get the machine electromagnetic characteristics with experimental evaluation, as numerical or analytical methods are needed to perform this task. As an alternative to experimental tests, the finite element method (FEM) is used for electromagnetic characterization of SRMs, providing results with great precision, but needing for a sophisticated mathematical implementation and large amount of computational resources to process the solution [22]. However, with the existence of FEM-based commercial software, it is possible to evaluate electromagnetic systems without a deep knowledge of electromagnetism. Thus, FEM is getting more adopted in the design of electric machines.

The mathematical model of switched reluctance generator results from the analysis of the associated power converter. Since all electric phases are defined by identic magnetic and electric circuits, they are represented by the same dynamic equations. Its analysis is generalized following the notation specified in Figure 2. The voltage across each electric phase

where for each phase

with

Applying the chain rule to Eq. (8), one has:

where

The electromotive force, which results from the change of magnetic flux in phase

During the conversion period, electric energy is exchanged between the capacitor of the converter and the electric phases of the generator. In the excitation period, the capacitor supplies energy to the phase and in the generation period receives energy from it. The voltage across capacitor

The total net value of the electric energy flowing between the electric phases and the capacitor

The bus voltage

In order to fully define the model of the switched reluctance generator system, the voltage across each electric phase must be known. Its value depends on the different circuit configurations for the electric energy flow in the phase. According to the possible combinations of switch states, there are three different configurations, as illustrated in Figure 3, each of them corresponding to a distinct period of electromagnetic conversion.

The voltage across each electric phase, according to the switch states, is:

with

The linear force exerted by the generator

where

All the formulation presented so far was defined for a phase electric position falling within the two opposing nonalignment positions. As the movable part of the generator is displaced, the electric position describes a periodic profile which is identical in all electric phases of the SRM. However, these electric positions have distinct values for the same absolute position because the phases are shifted. The relation between the electric position of one phase

where

The graphical representation of Eq. (18) can be found in Figure 4.

The mathematical model of the LSRG system is defined by a nonlinear problem of initial value and is schematized in Figure 5. To obtain the respective solution, it is necessary to perform the time integration of the differential equations that govern the system electromechanical behavior. Also, a controller must be included to provide the appropriate switching commands to the electronic power converter.

The bus voltage of the LSRG converter should be kept as constant as possible to allow for a proper self-excitation operation. Also, since high voltage levels are required to energize each phase in due time, it may be desirable to reduce the output voltage applied to the electric load. To achieve these requirements, an H-bridge isolated DC/DC converter is admitted in the conversion system to control the energy flow between the conversion system and the electric load, as schematized in Figure 6. With the ideal model of the H-bridge isolated DC/DC converter, it is possible to establish a relation between the mean values of the input and output voltages (

The duty cycle

where, respectively,

## 4. Control

The SRM operation relies on the switching electric positions. Fixed values for these positions may cause system instability, especially when operating at variable velocity. As a consequence, the converter bus voltage may change considerably, depending on the system electric load [24]. Therefore, a closed loop control is required. As stated in [25], when operating as generator, the control of SRMs must be applied to preserve the converter output voltage, conditioned by the current flow in the electric phases and electric load.

The amount of energy extracted for conversion relies on the electromechanical force exerted by the generator which is also affected by the electric phase current. Thus, with a proper switching strategy, it is possible to control the electric current intensity to attain a desirable voltage level as well as to improve the system conversion efficiency. With the additional DC/DC converter, the bus voltage may also be kept near a nominal value by modifying the energy flow between the LSRG converter and the load. This method is proposed in [24] to achieve maximum energy conversion by controlling the voltage level according to the velocity of the generator.

At low operation velocities, the back-electromotive force is inferior to the bus voltage and the current is forced to decrease gradually. In this situation, the phase must be submitted to successive commutations to adjust the applied voltage and achieve the desirable current intensity [7]. This process can be accomplished with hysteresis band control [26]. This control is characterized with a variable switching frequency because it is conditioned by the rate of change of the electric phase current [27].

The control of the switched reluctance generator is the selection of appropriate parameter values that are responsible for its behavior. For velocities superior to the base velocity, it is only necessary to account the electric positions that define the conversion cycle. For lower velocities, it is also required to define a reference value for the phase current. With the mathematical model, it is possible to establish an optimal relation between the values of the control parameters and the physical entities that need to be controlled and include it in the control process through look-tables or fitted analytical expressions. Thus, the appropriate control parameters can be defined as function of the operation variables [28].

The proportional integral (PI) control may be applied to control the switched reluctance generator. It has been used in real-time optimal control of rotating generators, where the reference current and commutation angles are computed from the error of the output voltage [29]. PI control was adopted in [30] to adjust the phase current intensity and minimize the converter output voltage ripple in linear generators.

### 4.1 Hysteresis control

With proper switching commands, the current is maintained within the hysteresis band

The conversion cycle is characterized by successive excitation and generation periods that occur while the phase electric position is between

### 4.2 Proportional integral (PI) control

The H-bridge isolated DC/DC converter is used to maintain the system DC bus voltage level close to a reference value. In order to accomplish it, a PI control is used to properly compute the duty cycle value for the applied switching commands. The control variable

with

where

with

## 5. Practical case study scenario

The described formulation will be applied to evaluate the performance of a tubular linear SRM working as generator in an ocean wave point absorber device schematized in Figure 9.

The system comprises a floating body that drives the generator by action of the incoming ocean waves. The linear generator is a three-phase machine with coils located in the inner part. The outer part is rigidly coupled to a floating body that only allows a vertical motion. Mechanical springs are used to link the generator outer part to the reference system, which is fixed to the ocean bottom. The tubular LSRG is illustrated in Figure 10.

### 5.1 Point absorber mathematical model

The mathematical model of an ocean wave direct drive converter is also needed to fully define the conversion system and to compute the mechanical entities that are used as input in the generator dynamic model.

The dynamic behavior of the point absorber device can be found in [32] and is described by the expression:

where

The excitation force is given by:

with

The radiation force is calculated as follows:

In Eq. (27),

The hydrostatic force is:

with

The floating body is subjected to a viscous drag force given by:

where

### 5.2 System simulation/simulation results

A sinusoidal ocean wave with a period of 7.7 s and a height of 1.3 m was considered to drive the floating body of the point absorber. The floating body is cylinder with a diameter of 3.4 m and a vertical extension of 4.9 m. The TLSRG was designed with a mean damping force of 40 kN, which is the maximum load expected for the ocean wave considered for simulation. The dimensions of the TLSRG are presented in Table 1. The electromagnetic characteristics of the machine, displayed in Figure 12, were obtained with MagNet^{®}, a commercial software that uses the finite element method for electromagnetic analysis. The control parameters

5.1 | 230 | 499.5 | 0.1 | 1.7 | 170 | 152 | 222.4 | 46.7 | 13.9 | 9.3 | 107.3 |

These values were included in the mathematical model through 2D look-up tables and were computed, as function of the TLSRG movable part velocity and phase current, by linear interpolation. The capacitor in the converter was defined with a capacitance of 0.05 F.

A value of 41,430 kg was assumed for the mass of the oscillating body ^{3}, and for the viscous drag coefficient, a value of 0.88 was defined. The excitation force profile and impulse response function were computed with NEMOH and are displayed in Figure 14. For the H-bridge isolated DC/DC converter, the ratio ^{®} and was solved by the Dormand-Prince method, with a relative tolerance of 1 × 10^{−3} and a maximum step of 8 × 10^{−4} s. The system’s mathematical model was simulated for distinct values of

Attending to the simulation results, the greatest average of electric generated power was 7.6 kW for a reference current of 35 A. For this reference, current value also had the best electric efficiency of 43%. Figure 16 shows the oscillating body absolute position and velocity obtained from the simulations with

The electric phase current and electromechanical force profiles are shown, respectively, in Figures 17 and 18. The TLSRG converter DC bus voltage level is displayed in Figure 19. With the application of the DC/DC isolated converter, the DC bus voltage was kept near its reference value with a maximum error of 1.05%.

## 6. Conclusion

In this chapter, a mathematical model to simulate the dynamic behavior of linear switched reluctance generators was presented, to apply in ocean wave direct drive converters. The model was developed according with the circuit configuration of the H-bridge asymmetric converter, adopted as power electronic converter to control the energy flow in the machine when operating as a generator. A hysteresis controller was applied to maintain the electric phase current close to a reference value during each conversion cycle, in order to allow the control of the generator damping force. Also, a DC/DC isolated converter was admitted to adjust the energy flow between the H-bridge asymmetric converter and the system electric load, in order to keep the DC bus voltage level near its nominal value. A PI controller was proposed to control the pulse width of the DC/DC conversion stage. A practical study case scenario was considered where the generator mathematical model was applied to simulate TLSRG operating as a power take-off system in an ocean wave point absorber device. The mathematical model of the TLSRG was integrated with the dynamic equations of the point absorber to evaluate the system behavior. In the simulation process, a regular ocean wave was assumed to drive the system. The system performance was evaluated for distinct values of

## Acknowledgments

The authors wish to thank Instituto de Telecomunicações for funding this work, with the investigation scholarship UID/EEA/50008/2019 and research grant BPD/N°19—29/01/2019—UID/EEA/50008/2019.