JP2016162311A - Sliding mode control method and sliding mode control device - Google Patents

Abstract

PROBLEM TO BE SOLVED: To provide a sliding mode control method and a sliding mode control device which are hardly affected by uncertainty, and have high robust property.SOLUTION: A sliding mode control device comprises: a plant 11; a plant nominal model 12; and a sliding mode controller (SMC) 20. The SMC is a controller for controlling a control object indicated by the plant by sliding mode control. The sliding mode control is a control mode for converging a state amount to a target value of the control object in a state in which, the state amount indicating a state of the control object is restricted on a changeover plane on a hyperplane (sliding mode), the sliding mode control is a non-linear control method which is hardly affected by uncertainty such as disturbance, parameter variation and modeling error, and has high robust property.SELECTED DRAWING: Figure 1

Description

  The present invention relates to a sliding mode control method and a sliding mode control apparatus.

A hyperplane for sliding mode control is set by a linear function with multiple state quantities to be controlled of the controlled object as variables. The state quantity is converged on the hyperplane, and the state quantity is constrained on the hyperplane. On the other hand, in the sliding mode control method for controlling to the target state quantity by converging the state quantity to the equilibrium point on the hyperplane, when the state quantity to be controlled is x ₁ and x ₂ , these state quantities are used as variables. The hyperplane H represented by σ = 0 is designed in advance using the linear function (σ = s ₁ x ₁ + s ₂ x ₂ ) as a function representing the convergence state on the hyperplane. Then, when the state quantities x ₁ and x ₂ are σ ≠ 0, the state quantities are rapidly converged on the hyperplane H (σ = 0) by high gain control according to a so-called reaching law (mode 1). Further, by so-called equivalent control input, the state quantities x ₁ and x ₂ are constrained on the hyperplane and converged to an equilibrium point (convergence point, point x ₁ = x ₂ = 0) on the hyperplane (mode 2). Is disclosed (Patent Document 1).

JP-A-9-274504

However, the mode 1 that converges to the hyperplane H (σ = 0) from the state in which the state quantities x ₁ and x ₂ are σ ≠ 0 is likely to be affected by disturbances, modeling errors, and the like. However, the above sliding mode control method has a problem that the robustness is lowered due to the influence of uncertainty.

  The problem to be solved by the present invention is to provide a sliding mode control method and a sliding mode control device which are not easily affected by uncertainty and have high robustness.

  In the present invention, the state quantity indicating the state of the controlled object is constrained by the switching surface represented on the phase plane, and the state quantity is set at the equilibrium point on the phase plane corresponding to the target value of the controlled object. In sliding mode control for controlling a controlled object by converging, when a plane whose axes are the state quantity and the first-order differential value of the state quantity is a phase plane, the locus of the quartic curve on the phase plane The above-mentioned problem is solved by generating a switching surface represented by the following and controlling a control object based on the switching surface. Here, the quartic curve passes through the point P indicating the initial state of the state quantity and the origin indicating the equilibrium point, has a tangent at the point P, and is a function related to the first order differential value of the state quantity and the state quantity.

  The present invention establishes a sliding mode that can converge to the equilibrium point while constraining the state quantity on the trajectory of the switching surface in the process from the initial state of the state quantity to the equilibrium state indicated by the equilibrium point. As a result, it is possible to provide a sliding mode control method and a sliding mode control device that are less affected by disturbances, modeling errors, and the like, and that are less affected by uncertainty and that are highly robust.

It is a block diagram of the sliding mode control apparatus which concerns on embodiment of this invention. It is the coordinate of the phase plane which showed the Lemnis skating curve of the switching surface. It is the graph which showed the torus body. It is the conceptual diagram which showed the sliding mode control in a 1st switching surface, and the sliding mode control in a 2nd switching surface on the phase plane. It is a flowchart which shows the control flow of the sliding mode control apparatus of FIG. It is the conceptual diagram which showed the state at the time of a state quantity converging to an equilibrium point in the conventional control method. It is the coordinates of the phase plane showing the trajectory of the state quantity for each attenuation coefficient. It is the coordinates of the phase plane showing the trajectory of the state quantity for each switching surface, and the Lemniska curve, ellipse, and line. (A) It is a graph which shows the error state-quantity characteristic of sliding mode control using a remnant skate, an ellipse, and a linear switching surface. (B) It is a graph which shows the characteristic of the speed of the error state quantity of sliding mode control using a remnant skate, an ellipse, and a linear switching surface. It is a graph which shows the characteristic of the energy consumption in the sliding mode control using a remuni skate, an ellipse, and a linear switching surface. It is a block diagram of a rudder angle control system concerning this embodiment. It is a block diagram of the sliding mode control apparatus which concerns on other embodiment of this invention. It is a graph which shows the characteristic of the switching gain in the sliding mode control using a remuni skate, an ellipse, and a linear switching surface. It is a graph which shows the characteristic of the energy consumption in the sliding mode control using a remuni skate, an ellipse, and a linear switching surface. It is a graph which shows the correspondence of an attenuation coefficient, a characteristic (switching gain, energy consumption, and quick response), and a switching surface. It is a flowchart which shows the control flow of the sliding mode control apparatus which concerns on other embodiment of this invention.

  Hereinafter, embodiments of the present invention will be described with reference to the drawings.

  FIG. 1 is a block diagram of a sliding mode control apparatus according to an embodiment of the present invention. The sliding mode control device and control method of this example are applied when, for example, steering angle control or the like in automatic parking of a vehicle is performed.

  As shown in FIG. 1, the sliding mode control device includes a plant 11, a plant nominal model 12, and a sliding mode controller (SMC) 20.

  The plant 11 is an actual control object. The plant 11 outputs an output value (y) to the control input (u) of the SMC 20. The plant nominal model 12 is an ideal model of a physical model to be controlled. The plant nominal model 12 outputs an ideal output value of the controlled object with respect to the input target value (r).

The SMC 20 is a controller for controlling the control object represented by the plant 11 by sliding mode control. The sliding mode control is a control mode for converging the state quantity to the target value of the control object in a state where the state quantity indicating the state of the control object is constrained by the switching surface on the hyperplane (sliding mode). This is a non-linear control method that is not affected by uncertainties such as disturbances, parameter fluctuations, and modeling errors, and has high robustness. In the sliding control mode, a linear state feedback term (u _L ) called an equivalent control input is calculated, a nonlinear feedback term (u _NL ) called a switching control input is calculated, and the feedback terms (u _L , u _NL ) are calculated. The added value is calculated as a control input to the control object.

  As shown in FIG. 1, the SMC 20 includes a CPU, a RAM, and the like, and as a functional block for executing the sliding mode control, a subtracter 21, an error observation estimation unit 22, a control region determination unit 23, and a switching The surface generation unit 24, the equivalent control input calculation unit 25, the switching control input calculation unit 26, and an adder 27 are included.

  The subtractor 21 calculates the error (e) by subtracting the state quantity (actual output value) of the plant 11 from the ideal output value (corresponding to the target value) of the plant nominal model 12. The error (e) is an error between the target value and the actual output value, and is also referred to as an error state quantity. An actual output value of the plant 11 is fed back and input to the subtractor 21.

  The error observation estimator 22 estimates the error velocity (primary differential value) and the error acceleration (secondary differential value) from the error (e) by an arithmetic process using a Kalman filter or the like. When the high-order derivative term of the state quantity such as the speed and acceleration of the controlled object cannot be observed, the error observation estimation unit 22 can estimate the high-order derivative term through arithmetic processing, thereby realizing a highly robust control. . Further, the error observation estimation unit 22 estimates the initial state of the error (e) by estimating the speed and acceleration of the error (e) from the initial value of the error (e), and similarly, the error during the sliding mode control. The state of the controlled object being controlled is estimated from (e). The error observation estimation unit 22 calculates the error (e) state (initial state and controlled state) from the observed value when the speed and acceleration of the controlled object can be observed.

  The control region determination unit 23 determines whether or not the error (e) is within the region where the switching surface is transitioned from the state of the error (e) estimated by the error observation estimation unit 22. As will be described later, the SMC 20 changes the switching surface according to the position of the state quantity on the hyperplane in a state where the state quantity (error e) is constrained by the switching plane on the hyperplane. The switching surface includes at least two types of switching surfaces, and a switching surface represented by a locus of a quadratic curve (hereinafter also referred to as a first switching surface) and a switching surface represented by a linear function (hereinafter, (Also referred to as a second switching surface).

  A region for switching from the first switching surface to the second switching surface is represented by an error e and a first-order differential value of the error. The area is set according to the required accuracy required for controlling the control object. For example, when an error that is about three times the required accuracy is set as an error, a range defined by the state quantity corresponding to the allowable error and the speed of the state quantity is an area for switching the switching surface. Then, the control region determination unit 23 determines that the sliding mode control is to be performed on the first switching surface when the current error e and the speed are out of the region, and when the current error e and the speed are within the region. Determines to perform the sliding mode control on the second switching surface. The control region determination unit 23 outputs the determination result and the estimated value of the error observation estimation unit 22 to the switching surface generation unit 24.

  The switching surface generation unit 24 is an arithmetic unit for generating a switching surface based on the error e, and includes a first switching surface generation unit 24a and a second switching surface generation unit 24b. The first switching surface generation unit 24a calculates a first switching surface that causes the error to reach from the point indicating the initial state of the error e to the transition point on the phase plane. The second switching surface generator 24b calculates a second switching surface that converges the error from the transition point to the equilibrium point on the phase plane. A transition point is a point which makes a transition from the first switching surface to the second switching surface.

  The switching surface generation unit 24 calculates the switching surface while selecting either the first switching surface generation unit 24a or the second switching surface generation unit 24b according to the determination result by the control region determination unit 23, and calculates the equivalent. It outputs to the control input calculating part 25 and the switching control input calculating part 26.

The equivalent control input calculation unit 25 calculates an equivalent control term (u _L ) from the input error e based on the switching surface generated by the switching surface generation unit 24 and outputs it to the adder 27. The equivalent control input calculation unit 25 calculates a switching control term (u _NL ) from the input error e based on the switching surface generated by the switching surface generation unit 24, and outputs it to the adder 27.

The adder 27 calculates a control input input to the plant 11 by adding the equivalent control term (u _L ) and the switching control term (u _NL ).

  The output of the plant 11 is controlled so as to follow the output of the plant nominal model 12 by the calculation of the above functional blocks.

  Next, the sliding mode control using the first switching surface will be described. First, the trajectory that minimizes the third-order differentiation (jump) of the state quantity and the trajectory of the Remnis skate curve, which is the trajectory of the quartic curve, will be described.

_Let us consider controlling the state of a certain plant from the initial state O to the target state L at a predetermined time tf. At this time, the smooth target trajectory with the controlled variable minimized is called the jerk minimizing trajectory (Non-patent literature: lash, T., Hogan, N .: The coordination of arm movements: An experimentally F confirmed mathematical model, Journal of Neuroscience, 5, 1688-1703, 1985).

In the control for making the state quantity of the plant follow the target value (hereinafter also referred to as target follow-up control), if the state quantity is x (t), the initial position 0, the elapsed time t, the arrival time t _f , and the target value L, The jerk minimization trajectory is expressed by the following equation (1).

  The jerk minimization trajectory is a smooth target trajectory in which the control amount is minimized. However, when the SMC 20 executes the sliding mode control using the curve of the jerk minimization trajectory as a switching surface, as shown in the equation (1), higher-order arithmetic processing is required. Therefore, in the present invention, the sliding mode control is performed by the switching surface represented by the trajectory of the Remnice skate curve, using the Lemnis skate curve function as a quartic function approximating the jerk minimization trajectory.

  Hereinafter, the derivation formula of the Lemnice skate curve will be described with reference to FIG. FIG. 2 is a coordinate system of a phase plane showing a Lemnice skating curve. In the phase plane, the horizontal axis (x-axis) is the state quantity, the vertical axis (y-axis) is the first-order derivative of the state quantity, and the state quantity is in the right half plane. In the following description, it is assumed that the state quantity exists in the right half plane of the phase plane. The straight line l of the Lemnice skating curve is a tangent at the point P. The case where the state quantity exists in the left half plane of the phase plane can be explained in the same manner as described below.

The Lemnis skate curve is represented by the following formula (2), which is an equation in the polar coordinate format, and is a curve obtained by cutting the donut-shaped torus body shown in FIG. FIG. 3 is a coordinate system showing a torus body. In addition, the torus body shown in FIG. 3 is a graph when a = 1 in the following (2).

When the curve shown by Formula (2) is represented by an orthogonal coordinate system, it is represented by the following Formula (3).

When the Remnant skate curve represented by the orthogonal coordinate system is enlarged by multiplying the x axis by a and multiplying the y axis by b, the following equation (4) is obtained. The relationship between the equations (3) and (4) is the same as the relationship between the equation of the circle and the equation of the ellipse obtained by multiplying the equation of the circle by the vertical axis and the horizontal axis, respectively.

Here, when both sides of the formula (4) are differentiated by x, the following formula (5) is obtained.

The phase plane shown in FIG. 2 is a plane representing the motion of the mass point system, where position (x) is the horizontal axis and velocity ( _{x_i} ) is the vertical axis. In the following, let us consider the trajectory of the Lemnice skating curve on the plane.

The slope of the tangent at an arbitrary point on the curve is expressed by equation (6).

_{x> 0, x _i> 0} , assuming the case of _{x _ii} <0 is considered the first quadrant. In addition, since the Lemnice skating curve is a line symmetrical figure with respect to the x-axis and the y-axis, the second quadrant, the third quadrant, and the fourth quadrant can be proved in the same manner.

As shown in FIG. 2, the Remnice skate curve passing through the point P (x, _{x_i} ) in the first quadrant is expressed by the following equation (7).

The slope of the tangent l of the curve at point P is expressed by equation (8) from equations (5) and (6). In the following formula, the point P is set to (x ₀ , x _{0 — i} ) as an initial state.

Equation (9) is obtained by transforming Equation (7).

By substituting equation (9) into equation (8), equation (10) is obtained.

Here, when b = as and equation (10) is transformed, equation (11) is obtained.

Further, by setting s ² = T and applying the solution formula of the quadratic equation, Expression (12) is obtained.

Further, the Lemnice skate curve has a tangent slope at the origin of ± 1. Therefore, the tangential slope (s) passing through the origin in the Remnant skate curve multiplied by a and b on the x-axis and y-axis, respectively, is expressed as s = b / a. The inclination near the origin can also be expressed as s = − (2ζ / ω) using the attenuation coefficient (ζ) and the angular frequency (ω) of the plant 11. Then, the following formulas (13) and (14) are obtained from the relationship between the formula (12) and the tangential slope.

The relationship as described above, through the origin indicating a P and equilibrium point points showing the initial state of the state quantity, as the locus of the quartic curve with a tangent slope at the point _{P (x _ii / x _i)} , lemniscate curve Can be derived.

ただし、α、βは、実モデルを表す係数である。ｕは制御入力である。 Next, when the switching surface is represented by a Lemnice skate curve, sliding control using the switching surface will be described using equations. First, a control input in sliding mode control (a control input to the plant 11 in FIG. 1) is derived. Here, it is assumed that the actual plant is represented by a secondary system shown in Expression (15).

Here, α and β are coefficients representing the actual model. u is a control input.

ただし、α_ｒ、β_ｒは、ノミナルモデルを表す係数である。γは参照入力値である。 Further, the nominal model is represented by equation (16).

Here, α _r and β _r are coefficients representing a nominal model. γ is a reference input value.

The error e can be expressed by equation (17).

The switching surface σ represented by the locus of the Lemnice skate curve is defined as in Expression (18).

The first-order differential value ( _{σ_i} ) of the switching surface σ is expressed by the following equation (19).

Further, the Lyapunov function V is defined as shown in Expression (20).

If the following equation (21) is satisfied, it is guaranteed that the state quantity is restrained by the switching surface (sliding surface).

Substituting Equation (18) and Equation (19) into Equation (21) yields Equation (22).

ただし、Ｋはゲインであり、ｓｇｎは符号関数である。 The expression in the parenthesis on the left side of the left side of the equation (22) is the equation of the Lemnice skate curve, and when the error state quantity (e, _{e_i} ) exists inside the closed curve of the Lemnis skate, the expression in the parenthesis on the left side of the left side is negative. Takes the value of On the other hand, when the error state quantity (e, _{e_i} ) exists outside the closed curve of the Lemnice skate, the expression in the parentheses on the left side of the left side takes a positive value. Therefore, in order for Formula (22) to hold, the formula in the parenthesis on the right side of the left side needs to have a sign opposite to the sign of the formula in the parenthesis on the left side of the left side. Therefore, the control input can be obtained by deriving u from the equation (22) while satisfying the condition indicated by the Lyapunov function.

However, K is a gain and sgn is a sign function.

Here, the term in curly brackets of the term including the sign function of Equation (23) is the switching control term (u _NL ), and the other term is the equivalent control term (u _L ). Then, referring to FIG. 1, the equivalent control input calculation unit 25 calculates the equivalent control term (u _L ) represented by the equation (23) as a linear feedback term, and the switching control input calculation unit 16 The switching control term (u _NL ) expressed by the equation (23) is calculated as a nonlinear feedback term having a gain (K). Further, the control input to the plant 11 is calculated by adding the target value (r) to the equivalent control term (u _L ) and the switching control term (u _NL ). Thereby, a control input in the sliding mode control is derived.

  While the state quantity (point P) is constrained by the first switching surface, the SMC 20 performs the sliding mode control by the above control method until the state quantity reaches the transition point. Then, after the state quantity reaches the transition point, the SMC 20 performs the sliding mode control using the second switching surface while changing the switching surface from the first switching surface to the second switching surface.

The transition of the switching surface and the sliding mode control after switching to the second switching surface will be described with reference to FIG. FIG. 4 is a conceptual diagram showing the sliding mode control on the first switching surface and the sliding mode control on the second switching surface on the phase plane. The vertical axis is the state quantity error (e), and the vertical axis is the first-order differential value ( _{e_i} ) of the error. E indicates a switching surface by a Lemnice skate curve, a point P indicates an initial state of the state quantity error, a point Q indicates a transition point, and L indicates a linear function linear switching surface. A region S represents a region where the sliding mode control is performed using the first switching surface. The equilibrium point is the origin.

ｑは所定の係数である。遷移点Ｑは、レムニスケート曲線Ｅと領域Ｓを示す楕円との交点である。 The region S is represented by an ellipse centered on the equilibrium point. The region S is represented by the following formula (24).

q is a predetermined coefficient. The transition point Q is an intersection of the Lemnice skating curve E and the ellipse indicating the region S.

  Until the transition point Q is reached from the initial state (point P) of the error state quantity, the error state quantity moves on the Lemmis skate curve E while being constrained on the switching surface represented by the trajectory of the Lemnis skate curve. When the error state quantity reaches the equilibrium point and the error state quantity enters the region S, the SMC 20 changes the switching surface from the first switching surface to the second switching surface.

  The second switching surface is represented by a linear function (linear function) passing through a point indicating the error state quantity and the origin. The slope of the linear function indicating the second switching surface is expressed by s = b / a. That is, the inclination of the second switching surface is an inclination of a tangent line passing through the origin of the Lemnice skate curve. And it moves on the straight line L, restraining on the switching surface represented by the straight line L within the region S.

  As described above, the SMC 20 has a plurality of switching surfaces when performing the sliding mode control, and transitions the switching surfaces when the state quantity reaches the transition point Q. The Lemnis skating curve is a closed curve. Then, when the sliding mode control is performed only on the switching surface of the Lemnice skate curve, and the state quantity does not stay at the equilibrium point, the state quantity draws a circular trajectory. In addition, since the Lemnis skating curve is a closed curve, the jerk increases when converging to the equilibrium point, and vibration may be induced. According to the present invention, since a plurality of different switching surfaces can be changed from a switching surface having a Lemnice skate curve to a switching surface having another trajectory according to conditions, the state quantity remains at the equilibrium point, and system stabilization can be realized. .

  Next, the control flow of the sliding mode control by the SMC 20 will be described with reference to FIG.

In step S1, the error observation estimation unit 22 of the SMC 20 observes or estimates the initial state of the error state quantity. In step S2, the switching surface generation unit 24 generates the first switching surface by the first switching surface generation unit 24a, and sets the switching surface in the sliding mode control as the first switching surface. In step S3, the SMC 20 performs an equivalent control term (u _L ) and a switching control term (u) by the equivalent control input calculation unit 25 and the switch control input calculation unit 26 in order to perform the sliding mode control along the first switching surface. _NL ) is calculated and added to calculate the control input to the plant and output it to the plant.

  In step S4, the error observation estimation unit 22 observes and estimates the current error state quantity, and the control region determination unit 23 determines that the current error state quantity is within the region S (region for switching the switching surface). It is determined whether or not there is. If the current error state quantity is outside the region S, the process returns to step S3.

  When the current error state quantity is within the region S, in step S5, the switching surface generation unit 24 generates the second switching surface by the second switching surface generation unit 24b, and selects the switching surface in the sliding mode control. Set to the second switching surface.

  In step S6, the error observation estimation unit 22 observes and estimates the current error state quantity, and the control region determination unit 23 determines whether or not the current error state quantity is within the region S. If the current error state quantity is outside the region S, the process returns to step S3. That is, when a large disturbance or the like is applied and the error state quantity is out of the region S, the first switching surface is generated again, and the sliding mode control using the Lemnice skating curve as the switching surface is executed.

  On the other hand, when the current error state quantity is in the region S, the SMC 20 performs the sliding mode control along the second switching surface (step S7). In step S8, the SMC 20 determines whether or not the current error state quantity has converged to the equilibrium point. If the current error state quantity has not converged to the equilibrium point, the process returns to step S6, and if the current error state quantity has converged to the equilibrium point, the control flow of the sliding mode control by the SMC 20 ends.

As described above, the control apparatus according to the present invention, in the phase plane, tangent through the origin indicating a P and equilibrium value points showing the initial state of the state quantity, and the slope at the point _{P (x} _ii _/ x _{_i)} A switching surface represented on a phase plane is generated with a locus of a quartic curve having a sine curve, and sliding mode control is executed along the switching surface. Thereby, it becomes difficult to be influenced by uncertainties such as disturbances and modeling errors, and a highly robust control method and control device can be realized.

  Incidentally, the conventional control method as shown in Patent Document 1 has the following problems. Hereinafter, the problem caused by the conventional control method will be described with reference to FIG. FIG. 6 is a conceptual diagram showing a state when the state quantity converges to the equilibrium point in the conventional control method, and is represented by phase plane coordinates.

The conventional control method performs control in two modes, an arrival mode and a sliding mode. The arrival mode is a control mode in which an initial value e (0) at a position away from the switching surface is constrained by the switching surface. The initial value e (0) is a value indicating the initial state of the error state quantity. The sliding mode is a control mode in which the error state quantity is converged to the equilibrium point in a state where the error state quantity is constrained by the switching surface of the linear function ( _{e_i} = −ae) passing through the equilibrium point. However, a is an inclination and a> 0.

  First, since the arrival mode is a control mode that is easily affected by uncertainties such as disturbance, the robustness is low. In the phase plane, since the initial value is in various positions depending on the state of the controlled object, the distance between the initial value and the switching surface changes according to the state of the controlled object and is also affected by uncertainty. Therefore, the controller cannot grasp the time until the initial state is constrained by the switching surface. Therefore, the controller cannot grasp the control timing of the controlled object.

  Furthermore, the reaching mode requires an excessive amount of control so that the state quantity is forcibly constrained on the switching surface of the linear function regardless of the dynamics of the system. The switching control input in such a control system becomes a feedback term having a very high gain nonlinear characteristic. The sliding mode is affected by the high gain of the nonlinear characteristic and the delay element of the system, and therefore a high frequency vibration phenomenon called chattering occurs.

  In the present embodiment, as described above, since the sliding mode control is established in the entire process from the initial state of the state quantity to the equilibrium point, the reaching mode does not occur. Therefore, it is possible to realize a control law that is not easily affected by disturbances and modeling errors, and the robustness is enhanced. In addition, when the uncertainties called matching conditions satisfy the conditions that can be expressed by the control input channel, the influence of the uncertainties can be nullified.

  In the present embodiment, in the phase plane, the state quantity moves along the switching surface until it converges from the initial state to the equilibrium point, so the SMC 20 can also calculate the convergence time until convergence to the equilibrium point. Therefore, in this embodiment, since the analytical solution of the convergence time from the initial state to the equilibrium state can be obtained by the geometric analysis, the control performance can be easily estimated.

  Further, in the present embodiment, during sliding mode control, the influence of high gain as in the arrival mode and the influence of system delay elements are also suppressed, so that chattering can also be suppressed.

  Next, the difference in energy efficiency will be described after comparing the case where the trajectory of the switching surface is a Lemnice skate curve, an ellipse, and a linear shape.

In general, the transfer function (G ₁ (s)) of the first-order lag system is expressed by the following equation (25). The transfer function (G ₂ (s)) of the second-order lag system is a transfer function of the second-order standard form and is represented by the following formula (26).

  T is the time constant of the first-order lag system, ζ is the attenuation coefficient of the second-order lag system, and ω is the angular frequency of the second-order lag system.

  In the second-order stabilization model expressed by the equation (26), when the damping coefficient (ζ) is set to 0.4, 0.6, 1.0, and 2.0, it converges from the initial state to the equilibrium point. The trajectory up to this is a trajectory as shown in FIG.

  FIG. 7 is a coordinate system of the phase plane showing the Remnice skate curve, the ellipse, and the linearity, respectively. In the phase plane, the horizontal axis (x-axis) is the error state quantity, and the vertical axis (y-axis) is the first derivative of the error state quantity. Point P indicates the state quantity in the initial state. The graph “Lemnis skate curve” represents a trajectory until convergence from the initial state to the equilibrium point when the switching surface is a lemnis skate curve. The graph “ellipse” indicates a locus when the switching surface is an ellipse, and “linear” indicates the locus when the switching surface is a straight line. When the model to be controlled is a first-order lag system, the state quantity converges to the equilibrium point while being constrained by the switching surface represented by the “linear” trajectory.

  As shown in FIG. 7, when the damping coefficient is large (for example, when ζ is 1.0 or more), that is, when the model is a second-order lag system having an overdamping characteristic, near the origin on the phase plane. The state quantity approaches the behavior of the first-order lag system. In FIG. 7, when the attenuation coefficient is large, when looking at the coincidence with the graph indicating the “linear” trajectory, “Lemni skating” is more “linear” than “ellipse”. Close to shape. Therefore, when sliding control is performed using a switching surface represented by an elliptic curve trajectory for a second-order lag system with overdamping characteristics, sliding is performed using a switching surface represented by a trajectory of a Lemmy skate curve. Compared with the case where the control is performed, the energy consumption rate can be suppressed by the sliding control based on the Lemnice skate curve.

  In addition, when the attenuation coefficient of the second-order lag model is 1.0, sliding control based on a switching surface of “Lemni skate”, “ellipse”, and “linear” will be described. FIG. 8 is a coordinate system of the phase plane showing the Remnice skate curve, the ellipse, and the linear shape, respectively, and is the same diagram as FIG. However, since the scales of the vertical axis and the horizontal axis are different, the shape of the graph only changes with respect to the graph of FIG. 7, and the characteristics of “Lemnis skate”, “oval”, and “linear” do not change. .

  When control is performed so that the state quantity converges from the initial state to the equilibrium point using the switching surfaces represented by the “Lemni skate”, “ellipse”, and “linear” loci, the state quantity trajectory As shown in FIG. In FIG. 8, the state quantity trajectories are graphs indicated by “sliding mode control using Lemnice skating”, “sliding mode control using an ellipse”, and “sliding mode control using a linear switching surface”.

  When the trajectory representing the switching surface is “linear”, an arrival mode is required to constrain the state quantity to the switching surface, and chattering occurs near the equilibrium point. Therefore, energy consumption increases. On the other hand, when the trajectory representing the switching surface is “Lem skating”, the sliding control can be performed from the initial state to the equilibrium point, and therefore the state quantity moves along the lemnis skating curve. Similarly, when the locus representing the switching surface is “ellipse”, the state quantity moves along an elliptic curve.

  Next, in the control for converging the state quantity to the equilibrium point shown in FIG. 8, the time characteristic of the state quantity, the time characteristic of the first derivative of the state quantity, and the time characteristic of the energy consumption are changed for each switching surface. Evaluated. The evaluation results are shown in FIGS. FIG. 9A is a graph showing the time characteristic of the error state quantity, and FIG. 9B is a graph showing the time characteristic of the first-order differential value of the error state quantity. FIG. 10 is a graph showing time characteristics of energy consumption.

As shown in FIG. 9A, when the trajectory representing the switching surface is “linear”, the time (t _c ) of the error convergence point is the time (t _a ) of the error convergence point of “Lemskating” and It is longer than the time (t _b ) of the error convergence point of “ellipse”. The error convergence point indicates a point in time when the error state quantity is close to zero.

As shown in FIG. 9B, when the trajectory of the switching surface is “linear”, chattering occurs in a long time until the error state quantity converges to the error convergence point (t _c ). .

As shown in FIG. 10, energy consumption amounts until the error state quantity converges to the error convergence point are represented as E _a and E _b . E _a is the energy consumption when the trajectory of the switching surface is “Lemni skate”, and E _b is the energy consumption when the trajectory of the switching surface is “ellipse”. In addition, when the trajectory of the switching surface is “linear”, the energy consumption amount until the error state quantity converges to the error convergence point is considerably larger than the energy consumption amount of “Lemni skate” and “ellipse”. Therefore, it is not shown in FIG.

The energy consumption (E _a ) when the trajectory of the switching surface is “Lemni skate” is lower than the energy consumption (E _b ) when the trajectory of the switching surface is “linear”, and the energy consumption (E _b ) About 20% lower.

  As described above, in the present embodiment, a switching surface having a remnant skate curve trajectory in accordance with the trajectory of the leni skate curve in accordance with the trajectory minimizing trajectory that is a smooth optimal trajectory is generated, and sliding control is performed based on the switching surface. It is carried out. Thereby, energy efficiency can be improved when performing sliding control for a system having an overdamping characteristic.

  Next, a case where the sliding control device according to the present embodiment is applied to steering angle control for automatic driving will be described as an example. The steering angle control for automatic driving is performed, for example, when a vehicle is parked in a predetermined parking space in non-contact charging of an electric vehicle. In non-contact charging, a positional shift between the ground-side power feeding coil and the vehicle-side power receiving coil greatly affects the energy transmission efficiency from the ground to the vehicle. Therefore, in non-contact charging, accurate positioning of the vehicle is required.

  FIG. 11 shows a block diagram of the steering angle control system. The steering angle control system includes a target trajectory calculation unit 101, an SMC (slide mode control) control unit 102, an actuator drive unit 103, an actuator 104, a vehicle 105, a state observation unit 106, and a vehicle position estimation unit 107. It has. In FIG. 11, the vehicle 105 and other system blocks are illustrated separately, but the target trajectory calculation unit 101 and the like are provided in the vehicle.

  The target trajectory calculation unit 101 calculates a target trajectory from the current position of the vehicle 105 to the target parking position. The SMC control unit 102 follows the target track from the target track calculated by the target track calculation unit 101 and the position of the vehicle 105 (actual vehicle position) estimated by the vehicle position estimation unit 107. While generating the switching surface, the sliding mode control is performed based on the switching surface to generate a control signal.

  The actuator driving unit 103 receives the control signal generated by the SMC control unit 102 and drives the actuator 104. The actuator 104 is a device that supplies power to the steering mechanism. When the actuator 104 is driven, the steering angle of the vehicle 105 is controlled.

  The state observation unit 106 observes the state of the vehicle using the steering angle sensor and the wheel speed pulse. The vehicle position estimation unit 107 estimates the current position of the vehicle using the observation value of the state observation unit 106. The vehicle position estimated by the vehicle position estimation unit 107 is fed back to the SMC control unit 102, and the SMC control unit 102 calculates the difference between the target value corresponding to the target track and the current vehicle position as an error state quantity. Then, the sliding mode control is performed based on the error state quantity.

  As a modification of the present embodiment, when the current state quantity (error state quantity) is within the region S, the switching surface generation unit 24 may adjust the inclination of the second switching surface. The inclination of the second switching surface is the inclination of a straight line (linear function) representing the second switching surface.

When the first switching surface generation unit 24a calculates the rem skate curve of the first switching surface, the switching surface generation unit 24 calculates the transition point Q that is the intersection of the preset region S and the lemnis skate curve. The transition point Q is an intersection with the Lemnice skate curve on the boundary line of the region S. _{Let the} transition point Q be a point (e ₁ , e _{1 — i} ).

When the error state quantity converges to the equilibrium point along the first switching surface, the error state quantity is a dynamic characteristic (e (t)) of a first-order lag having a time constant that is a slope of a linear function indicating the second switching surface. And is represented by the following equation (25).

βは第２切換面の傾きの変化速度である。 When the inclination of the switching surface is changed according to time, the dynamic characteristic (e (t)) of the error state quantity can be expressed by Expression (28).

β is the change rate of the inclination of the second switching surface.

  When the current state quantity is within the region, the second switching surface generation unit 24b generates the second switching surface by calculation using Expression (28). Accordingly, the second switching surface generation unit 24b adjusts the switching surface so that the inclination of the switching surface becomes smaller as the convergence time when the error state quantity converges to the equilibrium point elapses. In the region S, as the error state quantity approaches the equilibrium point (origin), the slope of the second switching surface gradually decreases, so that the stability of control increases.

  In the above modification, the inclination of the second switching surface is adjusted according to the convergence time of the error state quantity. However, instead of the convergence time, the distance from the state quantity to the equilibrium point is used to change the second switching surface. You may adjust the inclination of. In addition, when adjusting the inclination of the second switching surface, the inclination in the initial state may be the inclination of the tangent line at the origin of the Lemnice skating curve.

  The switching surface generating unit 24, the equivalent control input calculating unit 25, the switching control input calculating unit 26, and the adder 27 correspond to the “calculator” of the present invention.

<< Second Embodiment >>
A sliding mode control apparatus and control method according to another embodiment of the present invention will be described. In the present embodiment, a switching surface represented by a remnant skate trajectory, a switching surface represented by an elliptical trajectory, and a switching surface represented by a linear function trajectory can be selected as compared to the first embodiment described above. However, the difference is that the sliding control is performed based on the selected switching surface. Other configurations are the same as those in the first embodiment described above, and the description thereof is incorporated.

  FIG. 12 is a block diagram of a sliding control apparatus according to another embodiment. The switching surface generation unit 24 includes a third switching surface generation unit 24c, a fourth switching surface generation unit 24d, and a fifth switching surface generation unit 24e. Other configurations included in the SMC 20 are the same as those in the first embodiment, and the plant 11 and the plant nominal model 12 are also the same as those in the first embodiment.

  The third switching surface generator 24c calculates a third switching surface that causes the error to reach from the point indicating the initial state of the error e to the transition point on the phase plane. The third switching surface is a switching surface represented by the trajectory of the Lemnice skate curve. The method for generating the switching surface in the third switching surface generation unit 24c is the same as in the first embodiment.

  The fourth switching surface generation unit 24d calculates a fourth switching surface that causes the error to reach from the point indicating the initial state of the error e to the transition point on the phase plane. The fourth switching surface is a switching surface represented by an elliptical locus. The elliptic curve that is the locus of the fourth switching surface is a quadratic function that approximates the jerk minimization trajectory and can be derived by the following derivation formula.

The ellipse representing the trajectory of the fourth switching surface shares either the major axis or the minor axis with the x axis on the phase plane, passes through the origin 0 and an arbitrary point P (x, _{x_i} ), Further, an ellipse having a tangent l at the point P is assumed. The slope of the tangent line l _(x _ii _/ x _{_i)} the value of the second derivative of the state quantity is a value obtained by dividing the first derivative of the state quantity. The ellipse is represented by the following formula (29).

A tangent vector at an arbitrary point on the elliptic curve can be expressed by Expression (30).

The tangent line l that touches the ellipse at point P is given by equation (31).

Expression (32) is obtained from Expression (30) and Expression (31).

Here, when the coordinates of the point P in the initial state are substituted into the equation (29) representing the ellipse and the equation is transformed, (33) is obtained from the relationship of the equation (32).

When b = as is set on the left side and the right side of equation (33) (s indicates a parameter), equation (34) is obtained by eliminating b.

Further, when Expression (34) is substituted into the middle side and the right side of Expression (33), Expression (35) is obtained.

In order to satisfy the causality of the elliptical locus in the phase plane, it is necessary to develop clockwise. Then, when having an initial state (point P) on the ellipse, in the case of x ₀ > a, the speed (x 0 — _ii ) is always negative because the direction is from the maximum speed to the minimum speed, and x ₀ < In the case of a, the acceleration (x 0 — _ii ) is always positive because the speed is from the minimum value to the maximum value.

As shown in the above formulas (29) to (35), in the coordinates of the phase plane, an arbitrary ellipse passes through the origin O and shares either the major axis or the minor axis with the x axis. If the initial state (point P) is given, an ellipse representing the fourth switching surface is uniquely given (in other words, the variables (a, b) in the equation (29) are changed to the coordinates (( e ₀ , e _{0 — i} , e _{0 — ii} ) or (x ₀ , x _{0 — i} , x _{0 — ii} )), so that an equation representing an ellipse is uniquely defined).

  The fifth switching surface generation unit 24e generates a fifth switching surface represented by a linear trajectory. The fifth switching surface is a switching surface represented by a linear function passing through the equilibrium point on the phase plane.

  The switching surface generator 24 selects a switching surface to be used for sliding mode control from the third switching surface, the fourth switching surface, and the fifth switching surface according to the dynamic characteristics of the nominal plant. A method for selecting the switching surface will be described later.

The equivalent control input calculating unit 25 and the switching control input calculating unit 26 calculate an equivalent control input (u _L ) and a switching control input (u _NL ) according to the selected switching surface. The calculation method of the equivalent control input calculation unit 25 and the switch control input calculation unit 26 when the third switching surface is selected is the same as that of the first embodiment.

When the fourth switching surface is selected, the equivalent control input calculating unit 25 and the switching control input calculating unit 26 are based on the derivation formula described below, and the equivalent control input (u _L ) and the switching control input (u _NL ) is calculated.

The actual plant, the nominal model, and the error e can be expressed by the equations (15) to (17) shown in the first embodiment, respectively, and the fourth switching surface σ expressed by the locus of the elliptic curve is expressed by the following equation: (36)

The first-order differential value ( _{σ_i} ) of the switching surface σ is expressed by the following equation (37).

And when satisfy | filling the conditions of Formula (21) shown in 1st Embodiment, it is guaranteed that a state quantity is restrained by a switching surface (sliding surface). When Expression (36) and Expression (37) are substituted into Expression (21), Expression (38) is obtained.

ただし、Ｋはゲインであり、ｓｇｎは符号関数である。 The expression in parentheses on the left side of the left side of the equation (38) is an elliptic equation, and when the error state quantity (e, _{e_i} ) exists inside the ellipse, the expression in the left parenthesis on the left side of the left side takes a negative value. . On the other hand, when the error state quantity (e, _{e_i} ) exists outside the ellipse, the expression in parentheses on the left side of the left side takes a positive value. Therefore, in order for Formula (38) to hold, the formula in the parentheses on the right side of the left side needs to have a sign opposite to the sign of the formula in the parentheses on the left side of the left side. Therefore, the control input can be obtained by deriving u of Expression (38) while satisfying Expression (38).

However, K is a gain and sgn is a sign function.

The term in the curly bracket of the first term of the equation (39) is the equivalent control term (u _L ), and the term in the curly bracket of the second term is the switching control term (u _NL ). Then, the equivalent control input calculation unit 25 calculates the equivalent control term (u _L ) represented by the equation (39) as a linear feedback term, and the switching control input calculation unit 26 has a predetermined gain (K). The switching control term (u _NL ) expressed by the equation (39) is calculated as the nonlinear feedback term.

  Next, a method for selecting a switching surface in the switching surface generation unit 24 will be described. The switching surface generation unit 24 uses the attenuation coefficient of the controlled object model as the dynamic characteristics of the nominal plant that serves as an index when selecting the switching surface. Further, the switching gain (K) and the energy consumption are dependent on the attenuation coefficient of the plant 11. Therefore, the switching surface generator 24 selects the switching surface according to the attenuation coefficient so that the switching gain for the attenuation coefficient or the energy consumption for the attenuation coefficient can be suppressed.

  13 and 14 show the characteristics of the switching gain (K) with respect to the attenuation coefficient of the plant 11 and the characteristics of energy consumption. FIG. 13 is a graph showing the characteristic of switching gain (K), and FIG. 14 is a graph showing the characteristic of energy consumption. In FIG. 13 and FIG. 14, “Lem skating” indicates characteristics when the switching surface is a Lemn skating curve and control is performed based on the switching surface. “Ellipse” indicates a characteristic when the switching surface is an elliptic curve and control is performed based on the switching surface. “Linear” indicates a characteristic when the switching surface is a linear function and control is performed based on the switching surface.

  As shown in FIG. 13, when the attenuation coefficient is 2 or less, the gain (K) when “linear” is the highest, and the gains of “Lemni skate” and “ellipse” are low. Therefore, when the attenuation coefficient is 2 or less and it is desired to suppress the gain (K), the switching surface generation unit 24 switches either the third switching surface or the fourth switching surface to the switching of the sliding control. Select as face. In addition, when the model transmission system has an underdamping characteristic (for example, when the attenuation coefficient is 1 or less), the gain (K) when “linear” is significantly higher than that of “Remni skate” and “ellipse”. It is high. Therefore, when the plant 11 has an underdamping characteristic, the switching surface generation unit 24 does not select the fifth switching surface.

  When the attenuation coefficient is greater than 2 and less than 4, the gain of “Lemni skate” is the highest. In addition, the “linear” gain is lower than “Remni skate” and “ellipse” in most of the range from 2 to 4 of the attenuation coefficient (the attenuation coefficient is about 2.5 or more). Yes. Therefore, when the attenuation coefficient is greater than 2 and less than 4, and when the highest priority is given to suppressing the gain (K), the switching surface generation unit 24 selects the fifth switching surface. In addition, since the gains of “Lemni skate” and “ellipse” are not extremely high, the switching surface generation unit 24 may select the third switching surface or the fourth switching surface.

  When the attenuation coefficient is 4 or more, the gain of “Lemni skate” and “ellipse” is higher than the gain of “linear”. Therefore, when the attenuation coefficient is greater than 2 and less than 4, and when the highest priority is given to suppressing the gain (K), the switching surface generation unit 24 selects the fifth switching surface.

  As shown in FIG. 14, when the attenuation coefficient is 2 or less, the energy consumption of “linear” is the highest, and the energy consumption of “Lemni skate” and “ellipse” is low. Therefore, when the attenuation coefficient is 2 or less and it is desired to reduce the energy consumption, the switching surface generation unit 24 replaces either the third switching surface or the fourth switching surface with the switching surface for sliding control. Select as. Further, when the dynamic characteristic of the model has an underdamping characteristic, the gain (K) when “linear” is significantly higher than that of “Lemnis skate” and “ellipse”. Therefore, when the plant 11 has an underdamping characteristic, the switching surface generation unit 24 does not select the fifth switching surface.

  When the attenuation coefficient is greater than 2 and less than 4, the energy consumption of “Lemni skate” is within the range from 2 to 4 of the attenuation coefficient (range where the attenuation coefficient is about 3.3 or more). ) Is lower than “ellipse” and “linear”. Therefore, when the attenuation coefficient is greater than 2 and less than 4, and the highest priority is to suppress the energy consumption, the switching surface generation unit 24 selects the third switching surface.

  When the attenuation coefficient is 4 or more, the energy consumption of “linear” is the lowest, and the energy consumption of “Lemni skate” and “ellipse” is high. Therefore, when the attenuation coefficient is 4 or more and it is desired to reduce the energy consumption, the switching surface generation unit 24 selects the fifth switching surface as a switching surface for sliding control. Further, since the energy consumption of “Lemni skate” and “ellipse” becomes larger when the attenuation coefficient exceeds 4, the switching surface generation unit 24 does not select the third switching surface and the fourth switching surface. .

  Moreover, the switching surface generation part 24 may select a switching surface according to an attenuation coefficient so that a quick response may be improved so that it may demonstrate below. The quick response corresponds to the time until the state quantity converges to the equilibrium point, and the shorter the time until convergence, the higher the quick response.

The state quantity speed ( _{e_i} or _{x_i} ) when the state quantity converges to the equilibrium point will be described with reference to FIG. 7 shown in the first embodiment. The speed of the state quantity corresponds to the speed at which the equilibrium point converges. The switching surface having the largest absolute value of the maximum speed of the state quantity is “ellipse” among “Lemni skate”, “ellipse”, and “linear”. For this reason, the switching surface generation unit 24 selects the fourth switching surface when priority is given to responsiveness. As for the responsiveness, even if the trajectory of the switching surface is “Lemnis skate” or “linear”, sufficient responsiveness can be expected, so the switching surface generator 24 is not limited to the fourth switching surface, The switching surface or the fifth switching surface may be selected.

  As described above, the switching surface generation unit 24 has the third switching surface, the fourth switching surface, and the second switching surface so that the switching gain is low, the energy consumption is low, or the responsiveness is high. A switching surface for sliding control is selected from the five switching surfaces.

  FIG. 15 shows the relationship between the attenuation coefficient, the feature, and the switching surface. Features are energy consumption, switching gain, and responsiveness. In FIG. 15, “◎”, “○”, and “△” are selectable switching surfaces, “◎” is most suitable for selection, and “○” is suitable for selection in the order of “△”. It shows that. “×” cannot be selected. In the description of the present embodiment, “x” cannot be selected, but a selection surface corresponding to “x” may be selected according to an actual model.

  When selecting the switching plane so that the switching gain is minimized, the switching plane generator 24 selects as follows. When the attenuation coefficient is 2 or less, the switching surface generation unit 24 selects either the third switching surface or the fourth switching surface. When the attenuation coefficient is greater than 2 and less than 4, the switching surface generator 24 selects the fifth switching surface. When the attenuation coefficient is 4 or more, the switching surface generation unit 24 selects the fifth switching surface.

  When the switching surface is selected so that the energy consumption is minimized, the switching surface generation unit 24 selects as follows. When the attenuation coefficient is 2 or less, the switching surface generation unit 24 selects either the third switching surface or the fourth switching surface. When the attenuation coefficient is greater than 2 and less than 4, the switching surface generator 24 selects the third switching surface. When the attenuation coefficient is 4 or more, the switching surface generation unit 24 selects the fifth switching surface.

  When the switching surface generation unit 24 selects the switching surface so that the quick response is maximized, the switching surface generation unit 24 selects the fourth switching surface.

  Next, the control flow of the sliding mode control by the SMC 20 will be described with reference to FIG.

  In step S21, the error observation estimation unit 22 of the SMC 20 observes or estimates the initial state of the error state quantity. In step S22, the switching surface generation unit 24 uses the third switching surface, the fourth switching surface, and the fifth switching surface by the third switching surface generation unit 24c, the fourth switching surface generation unit 24d, and the fifth switching surface generation unit 24e. Is generated.

  In step S <b> 23, the switching surface generator 24 selects a switching surface to be used for sliding mode control from the third switching surface, the fourth switching surface, and the fifth switching surface according to the attenuation coefficient of the plant 11. For example, the switching surface generation unit 24 selects the switching surface so that the energy consumption is minimized.

  In step S24, the SMC 20 performs sliding mode control along the selected switching surface. In step S25, the SMC 20 determines whether or not the current error state quantity has converged to the equilibrium point. If the current error state quantity has not converged to the equilibrium point, the process returns to step S24. If the current error state quantity has converged to the equilibrium point, the control flow of the sliding mode control by the SMC 20 ends.

  As described above, in this embodiment, a switching surface used for sliding control is selected from the third switching surface, the fourth switching surface, and the fifth switching surface. Thereby, when performing sliding control, the optimal switching surface according to the characteristic of a model can be selected.

  Further, in the present embodiment, the third switching surface, the fourth switching surface, and the fifth switching surface are controlled according to the dynamic characteristics of the nominal plant so that the energy consumed when controlling the controlled object by the sliding mode control is minimized. A switching surface is selected from the switching surfaces, and a switching surface used for sliding control is selected. Thereby, energy consumption can be reduced.

  Further, in the present embodiment, the switching used for the sliding control is selected from the third switching surface, the fourth switching surface, and the fifth switching surface according to the dynamic characteristics of the nominal plant so that the switching control gain is minimized. Select a face. Thereby, a gain can be suppressed.

  In this embodiment, the time until the state quantity converges to the equilibrium point is shortened. According to the dynamic characteristics of the nominal plant, the third switching surface, the fourth switching surface and the fifth switching surface are selected, and the switching surface used for sliding control is selected. Thereby, responsiveness can be improved.

  Note that, as a modification of the present embodiment, the switching surface generation unit 24 performs a third operation according to the dynamic characteristics of the nominal plant so that the energy consumption, the switching control gain in the sliding mode control, and the responsiveness are optimized. You may select from the switching surface, the 4th switching surface, and the 5th switching surface, and may select the switching surface used for sliding control.

  As shown in FIG. 15, for example, if the “ellipse” or “linear” is selected when the attenuation coefficient is greater than 2 and less than 4, the energy consumption becomes high. On the other hand, when the “Lemnis skate curve” is selected, the switching gain is higher than “linear”, but the increase in the switching gain is small when compared with “linear”. In addition, when the “Lemni skate curve” is selected, the responsiveness is slightly worse than the “ellipse”. When the “Lemnis skate curve” is selected, the energy consumption can be greatly reduced as compared with “ellipse” and “linear”. Therefore, in the modification, when the attenuation coefficient is greater than 2 and less than 4, “Lemni skate” is selected.

  In addition, when the attenuation coefficient is 2 or less, selecting an ellipse can obtain an optimum result with switching gain, energy consumption, and quick response. Therefore, if the attenuation coefficient is 2 or less, an ellipse is selected.

  In addition, when the attenuation coefficient is 4 or more, when “Lemnis skating” or “ellipse” is selected, the energy consumption is significantly higher than “linear”, and the switching gain is also increased. Therefore, when the attenuation coefficient is 4 or more, “linear” is selected. Thereby, energy consumption can be suppressed, sliding control can be performed with a low switching gain, and quick response can be ensured.

11 ... Plant 12 ... Plant nominal model 20 ... Sliding mode controller (SMC)
DESCRIPTION OF SYMBOLS 21 ... Subtractor 22 ... Error observation estimation part 23 ... Control area | region determination part 24 ... Switching surface production | generation part 25 ... Equivalent control input calculating part 26 ... Switching control input calculating part 27 ... Adder

Claims (12)

With the state quantity indicating the state of the controlled object being constrained by the switching surface represented on the phase plane, the state quantity is converged to the equilibrium point on the phase plane corresponding to the target value of the controlled object. In the control method by sliding mode control to control the control object,
Generating the switching surface;
Calculating an equivalent control input of the sliding mode control from the state quantity and the target value;
Calculating a switching control input of the sliding mode control from the state quantity and the target value;
Calculating a value obtained by adding the equivalent control input and the switching control input as a control input to the control object,
The phase plane is a plane with the state quantity and the first-order differential value of the state quantity as axes,
A sliding mode control method, wherein the switching surface is represented by a locus of a quartic curve having a tangent at the point P, passing through a point P indicating the initial state of the state quantity and an origin indicating the equilibrium point. .
However, the point P is expressed by the state quantity and the first-order differential value of the state quantity, and the slope of the tangent is a value obtained by dividing the second-order differential value of the state quantity by the first-order differential value of the state quantity. The quartic curve is represented by a function relating to the state quantity and a first-order differential value of the state quantity. The sliding mode control method according to claim 1,
The sliding mode control method, wherein the quartic curve is represented by a quartic function approximating a jerk minimizing trajectory. The sliding mode control method according to claim 1 or 2,
The sliding mode control method, wherein the quartic curve is a Lemnice skate curve. In the sliding mode control method according to any one of claims 1 to 3,
Determining whether the state quantity is within a predetermined region including the equilibrium point, and based on the determination result, switching to either one of the first switching surface and the second switching surface,
The switching surface includes the first switching surface and the second switching surface,
The first switching surface is represented by the quartic curve,
The second switching surface is a switching surface that passes through the equilibrium point and is represented by a locus of a function different from that of the first switching surface, and converges the state quantity to the equilibrium point within the predetermined region. The sliding mode control method characterized by this. In the sliding mode control method according to claim 4,
The second switching surface is represented by a linear function passing through the equilibrium point,
The sliding mode control method, wherein the slope of the linear function is adjusted according to a distance between the state quantity and the equilibrium point, or a convergence time when the state quantity converges to the equilibrium point. In the sliding control method according to any one of claims 1 to 5,
The sliding mode control characterized in that at least one of the first-order differential value of the state quantity and the second-order differential value of the state quantity is an estimated value of a higher-order derivative term estimated from the state quantity. Technique. In the sliding mode control method according to any one of claims 1 to 5,
Selecting a switching surface used for the sliding mode control from the third switching surface, the fourth switching surface and the fifth switching surface included in the switching surface;
The third switching surface is represented by a trajectory having the quartic curve as a remnant skate,
The fourth switching surface is represented by an elliptical trajectory passing through the point P and the equilibrium point,
The sliding mode control method according to claim 5, wherein the fifth switching surface is represented by a locus of a linear function passing through the equilibrium point. The sliding mode control method according to claim 7, wherein
The step of selecting the switching surface includes:
According to the dynamic characteristics of the nominal plant, the third switching surface, the fourth switching surface, and the fifth switching surface so that the energy consumed when controlling the controlled object by the sliding mode control is minimized. Select from
The sliding mode control method, wherein the nominal plant is a representative model of a physical model of the controlled object. The sliding mode control method according to claim 7, wherein
The step of selecting the switching surface includes:
In order to minimize the switching control gain in the sliding mode control, select from the third switching surface, the fourth switching surface and the fifth switching surface according to the dynamic characteristics of the nominal plant,
The sliding mode control method, wherein the nominal plant is a representative model of a physical model of the controlled object. The sliding mode control method according to claim 7, wherein
The step of selecting the switching surface includes:
Select from among the third switching surface, the fourth switching surface, and the fifth switching surface according to the dynamic characteristics of the nominal plant so that the time until the state quantity converges to the equilibrium point is minimized. And
The sliding mode control method, wherein the nominal plant is a representative model of a physical model of the controlled object. The sliding mode control method according to claim 7, wherein
The step of selecting the switching surface includes:
According to the dynamic characteristics of the nominal plant, the third switching surface, the fourth switching surface, and the fifth switching surface are optimized so as to optimize the energy consumption, the switching control gain in the sliding mode control, and the convergence time. Choose from
The energy consumption amount is an amount of energy consumed when controlling a controlled object by the sliding mode control,
In the switching control gain, the convergence time is a time until the state quantity converges to the equilibrium point,
The sliding mode control method, wherein the nominal plant is a representative model of a physical model of the controlled object. With the state quantity indicating the state of the controlled object being constrained by the switching surface represented on the phase plane, the state quantity is converged to the equilibrium point on the phase plane corresponding to the target value of the controlled object. In the control device for sliding mode control for controlling the control object,
A generator for generating the switching surface;
An arithmetic unit for calculating an equivalent control input of the sliding mode control from the state quantity and the target value;
A computing unit for computing a switching control input of the sliding mode control from the state quantity and the target value;
An arithmetic unit that calculates a value obtained by adding the equivalent control input and the switching control input as a control input to the control object;
The phase plane is a plane with the state quantity and the first-order differential value of the state quantity as axes,
The switching mode control device, wherein the switching surface is represented by a locus of a quartic curve passing through a point P indicating the initial state of the state quantity and an origin indicating the equilibrium point and having a tangent line at the point P. .
However, the point P is expressed by the state quantity and the first-order differential value of the state quantity, and the slope of the tangent is a value obtained by dividing the second-order differential value of the state quantity by the first-order differential value of the state quantity. The quartic curve is a function related to the state quantity and the first-order differential value of the state quantity.

Images