<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Mathematical Modeling of Motion of Iron Bird Target Node of Security Data Management System Sensors</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>National Aviation University</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Ukraine tachinina</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>@gmail.com</string-name>
          <email>lysenko.a.i.1952@gmail.com</email>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>National Technical University of Ukraine «Igor Sikorsky Kyiv Polytechnic Institute»</institution>
          ,
          <addr-line>Kyiv</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
      </contrib-group>
      <fpage>0000</fpage>
      <lpage>0001</lpage>
      <abstract>
        <p>The article is devoted to presentation of the results of mathematical modeling of motion of the target node, which imitates a mini-drone and is mounted on the carriage of iron bird (test bench) for testing the sensors of security data management system. It is demonstrated that mathematical modeling made it possible to predict the possibility of self-oscillations during modeling of protective maneuvers of mini-drones and to confirm effectiveness of the algorithmic method of test bench upgrading. The results of computer simulations of the bench operation show that self-oscillations and an unstable mode of operation of a hypothetical typical biaxial iron bird (test bench) may occur during iron bird modeling of both deterministic and stochastic motion of the minidrone. These self-oscillations did not occur earlier when modeling the motion of inertial objects, because the cross-links between the channels when iron bird modeling with low accelerations were not significant. The simulation results showed the existence of a scientific and technical problem that arises when trying to use the available iron birds for modeling the motion of maneuverable targets. To solve the scientific and technical problem, it is proposed to perform algorithmic adaptation of the servo drive control system.</p>
      </abstract>
      <kwd-group>
        <kwd>Mathematical Modeling of Motion</kwd>
        <kwd>Target Node</kwd>
        <kwd>Iron Bird</kwd>
        <kwd>Sensor Network</kwd>
        <kwd>Flying Robots</kwd>
        <kwd>Drones</kwd>
        <kwd>Security Data Management System</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1 Introduction</title>
      <p>
        Currently the availability of anti-drone defense is one of the most important
conditions for solving the problem of ensuring the security (sustainable functioning) of
national critical infrastructure [
        <xref ref-type="bibr" rid="ref1 ref10 ref2 ref3 ref4">1–4, 10</xref>
        ].
      </p>
      <p>The most topical task is protecting against unauthorized access into the protected
zone of individual mini-drones, its groups and a "mosquito" raid of mini-drones. One
way to solve this problem is to use a security data management system (SDMS).</p>
      <p>
        Source information about the target in such systems comes from stationary, mobile
or quasi-mobile sensor networks. This information is used in SDMS for subsequent
identification, inspection, localization (interception) and, if necessary, destruction of a
mini-drone (for example, by a narrow electromagnetic pulse) [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ].
      </p>
      <p>The critical elements of SDMS sensor networks are target sensors (TS). The
quality of information that is used to make a decision in the SDMS depends on the
efficiency of the TS. TSs can operate on one physical principle (homogeneous TSs) or
different ones (multi-sensors of a target or heterogeneous TSs). The most important
element of the TS is the target tracking unit (TTU). At the stage of validation of made
design decisions, TSs together with TTUs are sent for iron bird testing.</p>
      <p>The test benches of modeling flat (two-dimensional) or spatial (three-dimensional)
target motion are used for testing.</p>
      <p>The target’s maneuvers are modeled by the motion of bench’s target node (Target
Motion Simulators (TMS)) using tracking drives that operate along two or three
mutually perpendicular axes.</p>
      <p>
        The sensor under test is placed in a special rotary support, which is spatially
located relative to the target node, taking into account the scale of modeling. TTU of TS is
checked for stability of the target capture and quality of tracking. At the time of
writing, the corresponding specialized companies have available two- and three-axis
modeling test benches (iron birds) (see Fig. 1 [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ]).
      </p>
      <p>
        Fig. 1. Example of combination of two-axis hypothetic test benches on one support [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ]:
1– rotary support for TS; 2 – TMS of first two-axis test bench; 3 – TMS of first two-axis
second test bench.
      </p>
      <p>Existing test benches are designed for iron bird simulation of the motion of objects
that are less maneuverable compared to mini-drones. This means that the dynamics of
test benches themselves when simulating the target motion can be neglected. If you
simulate the motion of super maneuverable mini-drones using same test benches, then
significant dynamic errors arise. There are two ways to solve the problem of modeling
the motion of a target, which is a maneuverable mini-drone. The first way is the
creation of specialized test benches with improved dynamic characteristics and the second
– is the algorithmic modernization of existing test benches, which is economically
feasible. The article is devoted to the description of method for algorithmic
modernization of existing biaxial test benches.</p>
    </sec>
    <sec id="sec-2">
      <title>2 Problem statement</title>
      <p>The biaxial bench (iron bird) of simulation of the target motion is considered. The
drives of the iron bird allow moving the carriage with the mounted target assembly.
The carriage (target node) movements imitate the maneuvers of the target in relation
to SDMS target sensor, which tracks this movement. The equipment mounted on the
carriage imitates the target’s identification features (for example, unmasking acoustic
or heat radiation, reflective properties of a mini-drone when it is illuminated with
radio-frequency radiation or laser locator light waves). We assume the easiest option
for fixing SDMS target sensor - the sensor is fixed stationary. The plane of
observation for the target motion is located at a certain distance from SDMS target sensor,
taking into account the scale and angle of modeling. The position of this plane is set
by the axes of movement of the iron bird target node.</p>
      <p>Typical movements of the mini-drone are assigned to the target node, then the
target tracking device, which is part of SDMS target sensor, is experimentally checked
for functioning.</p>
      <p>The iron bird target node should adequately simulate the motion of the mini-drone
relative to SDMS target sensor. The iron bird should not introduce significant errors
in the form of trajectory of the mini-drone and the rate of movement along it. This
means that the dynamic properties of the tracking drives of target unit should not
significantly affect the typical nature of the movement that is being modeled. If the
target is inertial, then the dynamic properties of tracking drives of the target node can
be neglected.</p>
      <p>However, the modern mini-drones have small mass and main moments of inertia
and are high-powered. This allows them to perform fast and complex maneuvers.
Such maneuvers were not available for those types for which the existing iron birds
were designed.</p>
      <p>We will check the possibilities of the existing biaxial iron birds for the adequate
execution of protective maneuvers of mini-drones on the computer mathematical
model of target node motion. All calculations and construction of the computer
mathematical models are performed using MATLAB + Simulink computer mathematics
system.</p>
      <p>
        A computer mathematical model of the typical dynamics of mini-drones motion
relative to the target sensor along ОХ and ОY axes of the observation plane of target
motion is set by Transfer Fon2 and Transfer Fon6 blocks (ref. Fig. 2) [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ]. ОХ and ОY
axes are located in the plane of observation of target movement: ОY axis is the line
perpendicular to the horizon plane and passing through the target node, ОХ axis is the
line passing through the target node and parallel to the line of intersection of the
observation plane of the target motion with the horizon plane.
      </p>
      <p>On XYGraph1 and XYGraph2 plotters, we can observe the trajectories of reference
motion of the mini-drone, and result of modeling this motion using a typical
hypothetical biaxial iron bird with tracking DC drives.</p>
      <p>The reference motion of mini-drone is a test movement when the target sensor will
be checked (tested). The mathematical model of drives that operate along OX and OY
axes is assumed to be the same and is specified in the simulation scheme (Fig. 2) by
State Spase1 and State Spase2 blocks (the model parameters are shown in the
diagram).</p>
      <p>Let us simulate the movement of mini-drone in the plane of observation of the
target movement. We assume that along each axis of mobility of the hypothetical iron
bird there is a follow-up digital drive with a digital proportional integral differential
controller.</p>
      <p>
        The transmission coefficients of the proportional, integral, and differential signals
of digital controller, which was used to simulate operation of the prototype test bench,
were assumed to be the same for both channels and equal 1, 0.08, and 1.47 (the
parameters of digital proportional integral differential controller were calculated by
Ziegler-Nichols method [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ]).
      </p>
      <p>
        Three types of the mini-drone reference motion were considered: “oval”;
"Lissajous figure"; "Random walk" [
        <xref ref-type="bibr" rid="ref3 ref4">3, 4</xref>
        ]. The result of modeling these reference motions is
shown in Fig. 3 (a, c), Fig. 4 (a, c), Fig. 5 (a, c). The dimension along ОХ and ОY
axes is meters.
Fig. 3. Graphic image of the reference trajectory of oval mini-drone (a) and result of
reproducing this movement with a hypothetical typical biaxial iron bird (b): position of P1 and P2
switches (Fig. 2) «down» and «up».
      </p>
      <p>Fig. 4. Graphic image of the reference trajectory of « Lissajous figure» mini-drone (a) and
result of the reproduction of this movement by a hypothetical typical biaxial iron bird (b):
position of P1 and P2 switches (Fig. 2) is «up» and «down», and P3 and P4 – «down» and «up».</p>
      <p>The results of computer simulations show that during iron bird modeling of both
deterministic and stochastic movements of the mini-drone, self-oscillations and an
unstable mode of operation of a hypothetical typical biaxial iron bird may occur.</p>
      <p>These self-oscillations did not occur earlier when modeling the motion of inertial
objects, because the cross-links between the channels in the iron bird modeling with
low accelerations were not significant.</p>
      <p>The simulation results showed the existence of a scientific and technical problem
that arises when trying to use the available iron birds to simulate the movement for
which these test benches were not designed.</p>
      <p>The easiest way to solve a scientific and technical problem is to perform an
algorithmic adaptation of the drive control system. We show how this can be done.</p>
    </sec>
    <sec id="sec-3">
      <title>3 Problem solving</title>
      <p>To eliminate self-oscillations of the target node that arise when modeling the
protective maneuvers of mini-drones, the algorithmic adaptation of servo drives
operating along the axes of mobility of a hypothetical iron bird was performed. The
algorithmic adaptation consisted in application of the modal control algorithm to correct
dynamic properties of the servo drives. The modal controller was synthesized
according to the following procedure:</p>
      <p>
        1. A continuous mathematical model of the target unit drive of a hypothetical
typical biaxial iron bird in the form of a MIMO LTI model was identified. Identification
performed using specifications of DC motors [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ].
      </p>
      <p>
        Based on the identification results, it was assumed that the MIMO LTI (Multiple
Input Multiple Output) drive models on the OX and OY channels practically coincide
and are described by equations in the two-dimensional state space [
        <xref ref-type="bibr" rid="ref9">9</xref>
        ]
dX/dt=A·X(t)+B·U(t)·
Y(t)=C·X(t)+D·U(t)
(1)
(2)
where X(t)=[x1(t); x2(t)], x1(t)=i(t) – current in the armature circuit of a DC motor,
x2(t)=w(t) – angular velocity of the armature rotation; U(t)=[u1(t); u2(t)], u1(t) –
control signal, which is fed to the independent excitation winding, u2(t) – signal
simulating the effect of cross-coupling between ОХ and ОY channels;
A=[-25 -7.5;7.5 0], B=[5 0;0 -5], C=[1 0;0 1], D=[0 0;0 0].
      </p>
      <p>2. Parameters of the drive’s discrete mathematical model are analytically
calculated. The functions MATLAB+Simulink «ss» and «c2d» of a computer mathematics
system were used for calculations.</p>
      <p>The period of quantization in time Т0 was calculated based on the following
reasoning.</p>
      <p>It is accepted that at the input of analog to digital converters, which are part of the
digital meters of the armature circuit current and angular velocity of rotation of the
armature of DC motors of OX and OY drive channels, the low-pass filters operate.</p>
      <p>Their upper cutoff frequency coincides with the upper cutoff frequency of the
useful signal and is equal to 51.4 rad/s for a typical drive. The time sampling period is
found T0 Ј p / 51.4 = 0.0611 s by the quantization theorem. It was accepted that
T0= 0.06 s.</p>
      <p>The discrete mathematical model of the drives of ОХ and ОY channels has the
following form</p>
      <p>X(n+1)=Ad·X(n)+ Bd·U(n)
Y(n)=C·X(n)+D·U(n)
(3)
(4)
where Ad=[0.1841 –0.2256; 0.2256 0.9359], Bd=[0.1504 0.04274; 0.04274 –0.2928].</p>
      <p>
        3. The modal controller for correction of the dynamic properties of the drives of
OX and OY channels was calculated [
        <xref ref-type="bibr" rid="ref7 ref9">7, 9</xref>
        ].
      </p>
      <p>The calculation results to the values of matrix amplification factors T, K,
inv(T*Bd) for:
1) modal analyzers Gain 7, Gain 10</p>
      <p>T=[-1.1860 -0.3955;-0.3955 -1.1860] (fig.6).</p>
      <p>Each modal analyzer converts the output signals of digital sensors of the drive motor
armature current and the angular velocity of this armature into virtual coordinates in
two-dimensional space;
2) modal regulators Gain 8, Gain 11</p>
      <p>K = [ -0.2407 0; 0 0.3607] (fig.6).</p>
      <p>Each modal controller calculates a virtual control vector.</p>
      <p>3) modal synthesizers Gain 9, Gain 12
inv(T*Bd)=[ -5.7616 1.1356;-1.9215 3.4057] (fig.6).</p>
      <p>Each modal synthesizer computes a real control vector.</p>
      <p>
        After connecting the modal controller, the digital proportional integral differential
controller was reconfigured in the servo system of ОХ and ОY channels according to
Ziegler–Nichols method [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ].
      </p>
      <p>The computer mathematical model and simulation results of functioning of the
modernized hypothetical typical biaxial iron bird are presented in Fig. 6, 7 (a, b), 8 (a,
b) and 9 (a, b).</p>
      <p>Fig. 6. Computer mathematical model of motion of the target node of a hypothetical typical
biaxial iron bird, which is designed to test the sensors of security data management system.
Fig. 9. Graphic image of the reference trajectory of “random walk” mini-drone (a) and result of
reproduction this movement by a hypothetical typical biaxial iron bird (b): positions of switches
(Fig. 6) P1, P3 and P2, P4 are “up” and “down”.</p>
    </sec>
    <sec id="sec-4">
      <title>4 Conclusions</title>
      <p>1. Mathematical modeling of the target assembly motion, which was virtually
driven by a hypothetical biaxial bench for testing the sensors of security data
management system, allowed: to predict the possibility of self-oscillations in the
simulation of protective maneuvers of mini-drones; confirm the effectiveness of algorithmic
method of upgrading the test bench.</p>
      <p>To implement the algorithmic method of upgrading the test bench it is necessary to
identify: a mathematical model of movement of the carriage (target node) of iron bird,
taking into account cross-links between the motion simulation channels; mathematical
model of extreme maneuvers of mini-drones.</p>
      <p>2. The complex algorithm of servo drives operation proposed in the article consists
of two algorithms: the digital modal controller algorithm and the digital proportional
integral differential controller algorithm. The first one corrects the dynamic properties
of electric motor, and the second – provides a quasi-adaptive control with respect to
the disturbing effect of cross-links.</p>
      <p>3. The iron birds, modernized by the algorithmic method described in the article,
can be used to test the tracking nodes for individual mini-drones or clusters of
minidrones (distributed drones). These tracking nodes can be used in data management
systems with a human operator or artificial intelligence in remote control loop of
inspection and destruction means of security data management systems.</p>
      <p>4. The direction of further research should be considered construction of a
mathematical model of the motion of target node, which will be installed on the carriage of
three-axis iron bird for testing the sensors of security data management system.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Yashchenko</surname>
          </string-name>
          , V.:
          <article-title>Intelligent robotic systems to ensure life safety in emergency situations of peace and war</article-title>
          .
          <source>Mathematical machines and systems</source>
          , vol.
          <volume>1</volume>
          , pp.
          <fpage>3</fpage>
          -
          <lpage>29</lpage>
          (
          <year>2020</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          <article-title>2. Kharakteristika i klassifikatsiya chrezvychaynykh situatsiy</article-title>
          , http://www.grandars.ru/shkola/bezopasnost-zhiznedeyatelnosti/klassifikaciyachrezvychaynyh-situaciy,
          <source>last accessed</source>
          <year>2020</year>
          /2/21.
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <given-names>Combat</given-names>
            <surname>Robots</surname>
          </string-name>
          <article-title>Will Outnumber Human Soldiers in the US Army by 2025</article-title>
          , https://www.industryleadersmagazine.com/ combat-robots
          <article-title>-will-outnumber-humansoldiers-us-army-2025, last accessed</article-title>
          <year>2020</year>
          /2/21.
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4. MQ-1 Predator, https://militaryarms.ru/voennaya-texnika/aviaciya/mq-1-predator/#gallery3,
          <source>last accessed</source>
          <year>2020</year>
          /2/21.
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5. Servotechnica, http://www.servotechnica.ru/project/index.pl?id_service=22, last accessed
          <year>2020</year>
          /2/21.
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Lysenko</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Tachinina</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Chumachenko</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          :
          <article-title>Mathematical formulation of the optimization problem of movement of a group of quadrocopters</article-title>
          .
          <source>Technical mechanics</source>
          , vol.
          <volume>1</volume>
          , pp.
          <fpage>73</fpage>
          -
          <lpage>78</lpage>
          (
          <year>2016</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Goodwin</surname>
            ,
            <given-names>G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Grebe</surname>
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Salgado</surname>
            <given-names>M.:</given-names>
          </string-name>
          <article-title>Proyektirovaniye sistem upravleniya [Design Control Systems]</article-title>
          .
          <source>Laboratory of Knowledge</source>
          , Moscow (
          <year>2004</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>German-Galkin</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          :
          <source>Matlab &amp; Simulink. Proektirovanie mehatronnykh sistem na PK [Matlab &amp; Simulink. Designing Mechatronic Systems on PC]. KORONA-Vek Publ., St. Petersburg</source>
          (
          <year>2011</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <surname>Dorf</surname>
            ,
            <given-names>R.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Bishop</surname>
          </string-name>
          , R.:
          <source>Sovremennyye sistemy upravleniya [Modern Control Systems]. Laboratory of Basic Knowledge Moscow</source>
          (
          <year>2002</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <surname>Lysenko</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Tachinina</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Chumachenko</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          , and it.:
          <article-title>Modeli primeneniya informatsionno-telekommunikatsionnykh tekhnologiy na osnove bespilotnykh aviatsionnykh kompleksov v chrezvychaynykh situatsiyakh [Models of application of informationtelecommunication technologies on the basis of unmanned aerial systems in emergency situations]</article-title>
          .
          <string-name>
            <surname>Lap-Lambert</surname>
          </string-name>
          Academic-Publishing (
          <year>2018</year>
          ).
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>