The distributed control system synthesis of elastic structure
- Authors: Koval’ V.A1, Torgashova O.Y.1, Stepanov M.F1
-
Affiliations:
- Yuri Gagarin State Technical University of Saratov, Saratov, Russian Federation
- Issue: Vol 25, No 4 (2017)
- Pages: 30-40
- Section: Articles
- URL: https://journals.eco-vector.com/1991-8542/article/view/20305
- DOI: https://doi.org/10.14498/tech.2017.4.%25u
- ID: 20305
Cite item
Abstract
For an elastic distributed plant with parameters depending on the spatial variable, based on the spectral method of distributed systems analysis and synthesis and considering the internal resistance according to Voigt the transition is made from PDE to an infinite system of ODE in state space form. The boundary conditions are additively included to the obtained spectral representation, which enables to control from the boundaries. The control law for suppression of oscillations is synthesized and the closed system is analyzed. The obtained results can be used in the control systems synthesis for aircraft with active dynamic compensation of elastic vibrations.
Full Text
Introduction
Elastic structures are widely used in the fields of aviation and rocket engineering to increase flight speed, reduce weight, and increase aircraft length. Under the appropriate flight conditions and particularly in terms of fuel consumption, elastic oscillations of the carrying structure occur; these are commensurate with the angular oscillations of an aircraft in terms of frequency and amplitude.
Elastic oscillations affect control system sensors and consequently, the control elements. These disturbances can result in the loss of accuracy and stability of the flight control [1, 2]. This creates an issue in the development of the aircraft control law to counteract external disturbances and elastic oscillations of the aircraft’s body.
Modern spacecraft are equipped with both rigid members and elastic structures, such as antennae, solar batteries, and outboard rods with metering instruments. These devices require passive or active stabilization for normal operation of the spacecraft.
The major control operations of elastic structures and their stabilization are based on the modern theory of analytical design of the best regulators for the systems with distributed parameters [3, 4].
It should be noted that for the objects described with the help of the system of differential or integro-differential equations with partial derivatives, the optimization system of the control-determining equations comprises a non-linear system of differential and integro-differential equations with partial derivatives [3]. The resolution of this system is a rather complicated task in both the implementation of computation procedures and the search for algorithms that provide a good convergence of the obtained solutions.
This study’s objectives are first, to perform a transition from differential equations with partial derivatives to an infinite system of standard differential equations in the form of the state space. This uses a spectral method of the theory of control [5, 6] for an elastic distributed object (an aircraft fuselage or a rocket body) with parameters dependent on the spatial variable and with regard to the inner resistance according to Voigt [2]. The second objective is to synthesize the law of control for the suppression of oscillations and perform an analysis of the closed system.
1. Mathematical model of the controlled object
It is assumed in this study that the elastic oscillations of both the aircraft fuselage and the rocket body are sufficiently accurately described by the equation of the flexible beam of the variable cross section with regard to the inner resistance as per Voigt, which, according to [2], appears as follows:
(1)
where is the spatial variable,
is time,
is the beam-axis bending measured in a perpendicular direction to the non-deformed beam axis,
is the length unit weight,
is the flexural stiffness,
is the elastic modulus,
is the moment of inertia of the beam cross section relative to the section’s neutral axis perpendicular to the oscillation plane,
is the external distributed transversal load attributed to the beam-length unit, and
is the coefficient of the inner resistance, as per Voigt.
Equation (1) is the following:
(2)
Let (2) be considered as a mathematical model of the controlled object with initial conditions
(3)
and boundary conditions
(4)
Let differential equation (2), initial conditions (3), and boundary conditions (4) be reduced to a dimensionless form. The dimensionless variables are introduced as
(5)
where are some nominal values of appropriate variables.
In the new variables in (5), the differential equation of the controlled object will be the following:
(6)
The coefficients of equation (6) are determined with expressions as
(7)
The initial conditions are
(8)
The boundary conditions are
(9)
Then, based on the properties of the spectral features in [5], the expressions for the matrix of spectral representation of the controlled object is obtained.
2. Spectral representation of the task
Let it be assumed that the function describing the controlled object’s condition is a substantial one-valued limited function with square integrability in the area of the spatial variable , and boundary conditions are applied at the points and
The function with regard to the boundary conditions can be presented as follows:
(10)
where is the function coinciding with the function at the interval ,
is the value of the single jump function at the boundary ,
is the value of the single jump function at the boundary , and
The generalized variable [7] of function (10) on x will be as follows:
For the m- variable, the following expression can be written as
.
Function shall be decomposed in the Fourier series according to the system of orthonormal functions on the variation interval
(11)
Using spectral-feature properties and with regard to , the transition from the differential equation with partial variables (2) at the initial conditions (3) and the boundary conditions (4) is performed to the system of standard differential equations as follows:
(12)
where is the vector of the spectral feature of function with components
are the infinite-dimensional square operational matrices of the first multiplier of the spectral properties of functions accordingly, whose elements are calculated with the expressions
are the infinite-dimensional square operational matrices of the multipliers whose elements are determined with the expressions
is the vector of the spectral property of function with the components
is the infinite-dimensional square operational matrix of differentiation with the elements calculated according to the expression
are vectors of the spectral properties of the boundary conditions with the elements
are operational matrices describing the jumps of function at interval calculated with the expression
Expression (12) is reduced to the following:
(13)
The new variable is introduced and represented with equation (13) in the form of the system of vector–matrix equations in Cauchy’s integral formula as
(14)
The value of moment and its derivative in time are considered as control activities at the object’s right boundary. The following designations are introduced:
(15)
Expressions (15) use the designations of vectors
With regard to designations (15), the system (14) can be written in the vector–matrix form as
, (16)
where is the vector of the condition,
is the vector of control, and
is the vector of disturbance.
Matrices A, B, and M are the following:
(17)
Thus, the transition from the description of the controlled object of equation with partial derivatives (2) with the set initial and boundary conditions (3) and (4) is performed to the system of standard differential equations (16) in the form of the state space with permanent coefficients.
Expression (16) is added with the following expression:
(18)
where is the vector of measured variables – values at the points where the sensors are installed, and D is the matrix, whose lines are composed of the orthonormal functions used for splitting into the Fourier series (11) and calculated at the points of measurement.
An analysis of the spectral representation of the controlled object (14) demonstrates that the boundary conditions are included in the object’s equation, which allows control to be performed from the object’s boundaries. Both bending moment and lateral force can be applied on the beam ends as the control activities; therefore, expressions (14) consider the members proportional to the second and third derivatives of function per the spatial variable x at the object’s right boundary.
3. Calculation of the matrix of the controlled object’s spectral representation and regulator synthesis
Let the distributed controlled object (2)-(4) with coefficients that are the functions of spatial coordinates with the following baseline data be considered as:
(19)
The distribution of weight and rigidity is
(20)
The dimensionless coefficients are obtained from (7) after the selection of the following nominal values:
The relative distribution of weight and rigidity will be the following for the selected values (20):
and the numerical values of the coefficients in (7) will be the following:
The following will be used as the system of orthonormal functions:
(21)
The values of the controlled-object matrix are provided, which consider (17). The calculations of matrix D also consider the fact that measurements are conducted at point :
(22)
(23)
The synthesis of the continuous control unit is conducted for the controlled object (16) with matrices (22) and (23) based on LQ-optimization and the theory of supervising devices in accordance with the procedure described in [6]. The control unit equation is expressed as
(24)
where is the vector of the control unit conditions and constitute the numerical matrix. The calculations consider the five amplitudes of the spatial modes.
The values of the control unit matrix are provided as
Figures 1 and 2 present the results of the closed-system analysis.
Fig. 1. Value of regulated variable at point x = 0.7
The diagram of the transition process presented in Fig. 1 demonstrates that the disturbance effect is compensated with error that does not exceed 3%. The control activities applied at the elastic object’s right boundary are active for 1.5 sec and do not exceed the module’s allowable values.
Fig. 2. The control actions. а: moment u1(t) at the object’s right boundary and б: time derivative u2(t) = ∂u1(t)/∂t at the object’s right boundary
Conclusion
Based on the spectral method of the theory of control, this study performed a transition from an equation with partial derivatives to an infinite system of standard differential equations in Cauchy’s form. The partial derivatives describe an aircraft’s elastic oscillations with regard to the inner resistance as per Voigt with irregular distribution of weight and rigidity in the structure.
Using LQ-optimization and the theory of supervising devices, the control unit was synthesized, and the watch unit was constructed with correction of the recovery error.
The obtained results can be used for the construction of aircraft control systems with active dynamic compensation of elastic oscillations. The application of these findings can facilitate an improvement in aircraft dynamics, a decrease in navigation errors, and a reduction in structural loads and stresses.
About the authors
Vladimir A Koval’
Yuri Gagarin State Technical University of Saratov, Saratov, Russian Federation
Author for correspondence.
Email: journal@eco-vector.com
(Dr. Sci. (Techn.)), Professor.
77, Politechnicheskaya st., Saratov, 410054, Russian FederationOlga Yu Torgashova
Yuri Gagarin State Technical University of Saratov, Saratov, Russian Federation
Email: journal@eco-vector.com
(Dr. Sci. (Techn.)), Professor.
77, Politechnicheskaya st., Saratov, 410054, Russian FederationMichail F Stepanov
Yuri Gagarin State Technical University of Saratov, Saratov, Russian Federation
Email: journal@eco-vector.com
(Dr. Sci. (Techn.)), Professor.
77, Politechnicheskaya st., Saratov, 410054, Russian FederationReferences
- Красовский А.А. Системы автоматического управления полетом и их аналитическое конструирование. - М.: Наука, 1973. - 560 с.
- Лебедев А.А., Чернобровкин Л.С. Динамика полета беспилотных летательных аппаратов. - М.: Машиностроение, 1973. - 616 с.
- Сиразетдинов Т.К. Оптимизация систем с распределенными параметрами. - М.: Наука, 1977. - 480 с.
- Дегтярев Г.Л., Сиразетдинов Т.К. Теоретические основы оптимального управления упругими космическими аппаратами. - М.: Машиностроение, 1986. - 216 с.
- Коваль В.А. Спектральный метод анализа и синтеза распределенных систем. - Саратов: Изд-во Сарат. гос. техн. ун-та, 2010. - 148 с.
- Коваль В.А., Торгашова О.Ю. Решение задач анализа и синтеза для пространственно-двумерного распределенного объекта, представленного бесконечной системой дифференциальных уравнений // АиТ. - 2014. - № 2. - С. 54-71.
- Гельфанд И.М., Шилов Г.Е. Обобщенные функции и действия над ними. - М.: Гос. изд-во физ.-мат. литературы, 1959. - 472 с.