<!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>Method of the Multi-UAV Formation Flight Control</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Chepizhenko Valeriy</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Skyrda Ihor</string-name>
          <email>skyrda2@gmail.com</email>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>National Aviation University</institution>
          ,
          <addr-line>Kyiv</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Ukrainian State Air Traffic Services Enterprise (UkSATSE)</institution>
          ,
          <addr-line>Boryspil</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
      </contrib-group>
      <abstract>
        <p>The main issue that arises when considering multi-UAV formation flight in a group is that of collision probability. In this case, without human control involved, it is artificial intelligence that is responsible for flight performance in the airspace in such a way that collision is avoided. Taking into account a rapid increase in civil and military applications of UAVs, a collision avoidance algorithm is proposed based on artificial potential field method. This method makes it possible to detect a potential conflict between multiple vehicles and other static or moving obstacles found in airspace, to provide collision resolution by changing UAVs flight parameters through maintaining minimum separation distance, including cases when manned vehicles are found in the same airspace. There can be distinguished a wide range of obstacles: static non-moving objects or vehicles having different sizes or flight parameters (multi-rotor, fixed wing and single rotor UAVs), or a few UAVs of one type but with different types of hardware configuration, at the same time considering the possibility of flight performance in the same airspace with manned aircraft. Group formation keeps shape on the flight path, taking into account some ground speed restrictions and turn bank angle values according to UAV's flight performance characteristics. The proposed method is used for multi-UAV control without any leader and provide multiple conflicts resolution, where each UAV is characterized by its protection zone.</p>
      </abstract>
      <kwd-group>
        <kwd>autonomous unmanned aerial vehicle</kwd>
        <kwd>artificial potential field</kwd>
        <kwd>synergetic</kwd>
        <kwd>formation flight control</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>Remotely Piloted Aircraft System (RPAS) or Unmanned Aircraft System (UAS),
colloquially known as ‘drones’, are aerial vehicles that fly without an on-board pilot, as
well as the systems that support them to do so. RPAS refers to a system, extending
beyond the Remotely Piloted Aircraft (RPA) or Unmanned Aerial Vehicle (UAV) to
include ground stations (where control units and remote pilots are based) and
communications infrastructures. Within the broad definition of UAS lies a diverse range of
systems and UAVs. Some differences between these UAS are immediately apparent,
such as the size or weight of UAVs. Other differences are more subtle, such as the
medium of communication between the vehicle and the ground station. These systems
have varying degrees of automation and autonomy, but usually include human remote
pilots who control the vehicle from meters, kilometers or continents away. Perhaps the
most established and visible applications of UAS are for military purposes, including
combat and surveillance operations, but many applications have been identified for
domestic uses such as environmental monitoring, security, emergency response,
surveillance and recreation. In addition to the significant functional and economic benefits of
these civil UAS, the technologies required for civil UAS operations are ready for
market and the principal barriers to development in the sector are regulatory. In response
to demand, the European Commission (EC) has published strategies to allow the
gradual integration of UAS into normal airspace.</p>
      <p>The main technical peculiarity of UAS is defined by the extent of autonomy and
automation delegated from the pilot to the system. Automation levels range from those
that are fully piloted from a remote location to those that are fully automated. There are
also several points in-between, with some maneuvers triggered automatically through
autonomous monitoring of conditions. Depending upon system priorities, autonomous
maneuvers may have priority over, or be overridden by, the commands of a remote
pilot. The International Civil Aviation Organization (ICAO) and current EC plans will
only permit autonomous maneuvers to override pilot command in extraordinary
circumstances such as communication failure or imminent collision risk. The UAS
technologies beyond this definition, featuring greater autonomy, are also quite well
developed and, while integration is not currently planned, it could plausibly follow a
successful period of development in the UAS sector [1].
2</p>
    </sec>
    <sec id="sec-2">
      <title>Analysis of Researches and Publications</title>
      <p>The results of analysis show that most of known methods for multi-UAV control have
a number of significant limitations that are connected with multiple conflicts resolution
and group formation. Particularly, the main disadvantage that potential conflicts can be
solved pairwise, when this issue needs to be done in a global way. For example, the
system called Traffic Alert and Collision Avoidance System (TCAS) that is already
installed aboard uses range measurements and range-rate estimates to determine if a
conflict exists [2]. Methods developed for group control in robotics do not include such
feature, so UAVs must deal with constant movement and limited turning ability, which
makes collision avoidance much more complicated [3].</p>
      <p>The classical approach is called geometric, in which aircraft trajectory predictions
are based on linear projections of the current vehicle states [4-5]. The major
disadvantage is that prediction errors are negligible only for short time periods and require
high rate of surveillance information update. The class of stochastic approaches is
related to the problem of probabilistic conflict detection in the presence of various
uncertainties during the flight. The aircraft dynamics are described by using stochastic
differential equations, and the future aircraft’s trajectory is determined by solving the
stochastic trajectory optimization problem, it could be applied for conflict definition at
rather big distances [6], so stochastic approach can hardly be applied in order to control
a group of UAVs flying close to each other.</p>
      <p>Linear programming is a mathematical method [7] where optimal control problem
lies in finding trajectories that minimize objective function. The drawback of such
approach is the flyability of the optimal trajectories as far as safety and performance
aspects of a given flight route are concerned.</p>
      <p>The common disadvantage of all these methods is that they do not meet the main
requirements with respect to autonomous UAVs: the absence of any communication
links with the appropriate ground stations, with on-board computational and power
sources being limited.</p>
      <p>The summarized disadvantages of the analyzed methods make it impossible to
simultaneously use a combination of such parameters as heading, speed and altitude change
maneuvers to resolve multiple potential conflicts. Therefore, it is necessary to develop
some new methods for multiple autonomous UAVs control in a group in a
three-dimensional space. The method developed in this article is the evolution of potential field
method proposed in article [8]. A potential fields approach is based on assigning
magnetic or electrical charges of the same sign to UAVs, while the opposite charges are
assigned to destinations, with the principle being based on the laws of physics according
to which the like particles will repel each other, while the destinations having opposite
charges will attract them. The main feature of such approach is that UAVs do not
necessarily need to know the positions of all other aircraft, so artificial force generated by
each UAV allows them to avoid each other spontaneously, at the same time keeping a
group form [9]. According to [10], this approach is scalable and can be applied to a big
number of UAVs, even in case of multiple conflicts.
3</p>
    </sec>
    <sec id="sec-3">
      <title>Problem Statement</title>
      <p>To solve this problem, a potential field approach is used. This method uses the property
of the real world charged particles to generate a force field (electric or magnetic), which
causes attraction and repulsion forces when these particles interact. The matter itself is
a typical example of the self-organization principle in nature. UAVs are considered as
dynamic objects with the same sign, with the point of destination having the opposite
sign, it is analogous to the free movement of the aircraft autonomous motion where
they constantly have potential conflicts, and it is required to avoid collisions with other
dynamic objects or static/dynamic obstacles. In this case, the term ‘potential conflict’
is a situation, when the minimum separation standard between dynamic objects is
violated. The protection zone of dynamic objects is generally defined as follows: the
minimum allowed horizontal separation and the vertical separation requirement depending
on the sizes of dynamic objects. The dynamic objects collision is the process of
interaction between the dynamic objects or obstacles at a distance in which the dynamic
objects change their direction of motion and the speed module.</p>
      <p>The dynamic objects interact similarly to the particles of substances that are found
in other aggregate states of matter (solid, liquid). The forces act simultaneously. For
different dynamic objects, the general character of the force of gravity from distance is
qualitatively the same: the force of attraction between dynamic objects dominates at
large distances, while the force of repulsion acts at short distances. Fig.1 shows the
qualitative dependence of interaction of forces between two dynamic objects found at
distance r between two dynamic objects is presented, where F+ and F- – are the
dependence of the attraction and repulsion forces respectively, and F++F- - is a resultant force.
At a critical distance r = rcr the resultant force is equal to zero, i.e. the forces of
attraction and repulsion are counterbalanced. This distance rcr corresponds to the equilibrium
distance between the dynamic objects.</p>
      <p>This article considers a group system consisting of n autonomous UAVs, with a
point-mass model used to describe UAV formation movement. The related variables
are defined with respect to the inertial coordinate system and are shown in Fig. 2.</p>
      <p>The point-mass UAV model captures most of the dynamical effects encountered in
civil aviation aircraft. The point-mass equations of motion are formulated with respect
to a coordinate system shown in Fig. 2. The point-mass model assumes that the UAV
thrust is directed along the velocity vector, and that the UAV always performs
coordinated maneuvers. It further assumes a flat, non-rotating earth. These assumptions are
reasonable for UAVs operating within different ranges, therefore, this method can be
used in conflict resolution between different types of UAVs, with the fidelity provided
by the point-mass model being adequate for formulating these problems.</p>
      <p>Point-mass models applicable for spherical earth approximations can also be
developed. The fuel expenditure is negligible, i.e. the center of mass is time-invariant [11].
Under these assumptions, the motion equations of the i-th UAV can be described as
follows:
 ̇ = 
 ̇ = 
 cos  cos  ;
 cos  sin  ;
ℎ̇ =   sin  ;
 ̇ =   cos   −   cos   ;
   
  sin  ;
    cos 
 ̇ =   −   −  sin  ;
 ̇ =</p>
      <p>where: i=1, 2, …, n is the index of multiple UAVs under consideration. xi, yi, hi denote
the components of UAV gravity center position. For i-th UAV, xi is down range; yi is
cross range; hi is altitude; Vi is ground speed; γi is flight path angle; χi is heading angle;
Ti is engine thrust; Di is drag; mi is mass; g is acceleration due to gravity; φi is bank
angle; Li is vehicle lift. Bank angle φi and engine thrust Ti are control variables for an
aircraft. Bank angle is commanded via combining rudder and aileron trims, thrust is
commanded by engine throttle. The g-load ni=Li/gm is controlled by elevator, though
it refers only to UAV construction characteristics having higher limits due to the
absence of crew on board an aircraft in comparison to traditional application. Throughout
the multi-UAV control process, these control variables will be constrained to remain
within their respective limits. The most common constraints considered are upper and
lower bounds on ground speed (Vi), altitude (hi), g-load (ni), thrust (Ti), bank angle (φi)
and climb or descent rates.</p>
      <p>
        Heading angle χi and flight path angle γi are computed as:
tan   =
tan   =
 ̇
 ̇
ℎ̇

 
(
        <xref ref-type="bibr" rid="ref1">1</xref>
        )
(
        <xref ref-type="bibr" rid="ref2">2</xref>
        )
(
        <xref ref-type="bibr" rid="ref3">3</xref>
        )
      </p>
      <p>In air traffic, conflict resolution is determined by separation constraints, forming the
so-called conflict envelopes or ‘protection zones’ so that UAVs flight trajectories do
not overlap during the flight. The conflict between two UAVs or an UAV with the
above-mentioned obstacles implies that their altitude should differ in value hpr given in
UAV flight performance characteristics, or they should not get closer in the horizontal
plane than indicated by value rpr. The protection zone can be visualized for each UAV
as shown in Fig. 3.</p>
    </sec>
    <sec id="sec-4">
      <title>Method of the Multi-UAV Formation Flight Control</title>
      <p>In order to apply this approach it is required to transfer the real world properties of
UAVs and their position coordinates to the virtual world with its synergetic properties,
with the potential conflicts that may occur on the flight path being taken into account.
[12-13].</p>
      <sec id="sec-4-1">
        <title>This process includes the following steps:</title>
        <p>- structural and parametric synthesis of the virtual world;
- structure formation and parameters of virtual measuring systems that provide
conflict free trajectories calculation.</p>
        <p>UAVs are transferred from real to virtual world as dynamic objects, with mass,
attraction and repulsion potentials values being assigned to them. So, the equilibrium
state can be represented as:</p>
        <p>+(  ,   ,  ,    ) =  −(  ,   ,  ,    )
where mi, mj – masses of i-th and j-th dynamic bodies, G – gravitational constant,</p>
      </sec>
      <sec id="sec-4-2">
        <title>Attraction and repulsion forces can be calculated as:</title>
        <p>Projections of attraction and repulsion forces between i-th and j-th bodies on axes X
and Y are calculated by the formulas:
  + =
  − =</p>
        <p>;


   кр;

α ϵ {2, 3, …};
β ϵ {3, 4, …};
+
Fijx = Fij</p>
        <p>+
Fijy = Fij
+ |xi−xj|
+ |yi−yj|
rij
rij</p>
        <p>−
Fijx = Fij</p>
        <p>
          −
Fijy = Fij
− |xi−xj|
− |yi−yj|
rij
rij
rij = √(xi − xj)2 + (yi − yj)2
(
          <xref ref-type="bibr" rid="ref4">4</xref>
          )
(
          <xref ref-type="bibr" rid="ref5">5</xref>
          )
(
          <xref ref-type="bibr" rid="ref6">6</xref>
          )
(
          <xref ref-type="bibr" rid="ref7">7</xref>
          )
(
          <xref ref-type="bibr" rid="ref8">8</xref>
          )
        </p>
        <p>2 2 = 1 ∑ ≠ (Fi+jx − Fi−jx +   )</p>
        <p>1 ∑ ≠ (Fi+jy − Fi−jy +   )</p>
        <p>∈  ,  ∈  .</p>
        <p>
          In equations (
          <xref ref-type="bibr" rid="ref5">5</xref>
          ) and (
          <xref ref-type="bibr" rid="ref6">6</xref>
          ), the aggregate state of the environment of the virtual world
(solid, liquid, gas) is chosen by the ratio α/β, which characterizes the degree of
selforganization of the dynamic objects. Analogy of the aggregate state of a virtual
environment can serve as an aggregate state of matter - gaseous, liquid, crystalline, etc.
        </p>
        <p>The resultant vector at each point of dynamic object location consists of the sum of
attraction and repulsion forces</p>
        <p>+ +   −, but can perform a group formation, so to
produce dynamic objects movement there should be present one more force which takes
into account thrust force Pijx, Pijy direction with projection on axes X and Y (Fig. 4):</p>
        <p>The main condition for dynamic object motion should be satisfied in the following
way:</p>
        <p>
          + +   − &lt;  (  ). The group consists of n dynamic objects and each of them can
be described by the system of equations:
(
          <xref ref-type="bibr" rid="ref10">10</xref>
          )
(11)
(12)
(13)
(14)
Fig. 4. The scheme of forces with four dynamic objects in the original position
The main advantage of the virtual world that was formed is that when the dynamic
objects approach the critical distance rpr, the resultant force acting on them is zero, i.e.
the forces of attraction and repulsion balance each other. Thus, rpr allows to set the size
of the dynamic objects protection zone.
        </p>
        <p>Fi+j = Fij
−
(15)</p>
        <p>The absence of intersections of such zones, taking into account the uncertainty of
the forecasted position of the dynamic objects, allows maintaining a guaranteed level
of traffic safety in the multi-UAV formation flight control (Fig. 5).</p>
        <p>If a static obstacle occurs on a multi-UAV path, the group interacts with it through
applying attraction   + and repulsion   − forces (Fig. 6). This type of maneuver can be
conducted provided   − is neglected, because the obstacle is static:
  + &lt; Fi+j + Fi−j +   (  )
(16)</p>
        <p>The values of heading angle χi and ground speed Vi may change depending on
dynamic objects location relative to the obstacle and destination point.
In order to find out if the potential field approach can be applied in the solution of the
problem of multi-UAV formation flight control, Matlab simulators were used. All in
all, 2 cases were simulated with a different number of dynamic objects, with UAV being
referred to as a dynamic object. In Experiment 1 (see Fig. 7), 8 dynamic objects were
considered with the point-mass of 1 kg and protection radius 3 m, with only one 6
mradius obstacle to overcome. In Experiment 2 (see Fig. 8), 12 dynamic objects were
considered whose point-mass was 1 kg and protection radius was 3 m, with three
obstacles in the way whose radii varied from 3.5 to 4.5 m. The path was divided into 3
main stages of flight: 1) group formation; 2) obstacle avoidance; 3) straight line flight
in a group to the destination. Figures represent dynamic objects movement trajectory
(a), distance between moving dynamic objects, with dotted line showing protection
zone with radius 3 m (b), heading angle χi (c) and change in ground speed Vi (d).
tan   =  ̇̇ or tan   = FFiijjyx
  = √ ̇ +  ̇
(17)
(18)</p>
        <p>
          The dynamic objects are in their original positions with the starting speed being
equal to zero. At the first stage of modelling, due to the action of attraction (
          <xref ref-type="bibr" rid="ref5">5</xref>
          ) and
repulsion (
          <xref ref-type="bibr" rid="ref6">6</xref>
          ) forces the process of group formation begins, which depends on the
distance between them (
          <xref ref-type="bibr" rid="ref9">9</xref>
          ). Heading angle χi has the same direction as vector Fij, which is
projected on axes X (
          <xref ref-type="bibr" rid="ref10">10</xref>
          ), Y (11) and is formed by their sum, including thrust force (12).
At the same time, the shape of group formation is regulated be the equilibrium state (
          <xref ref-type="bibr" rid="ref4">4</xref>
          ),
(15).
        </p>
        <p>Fig. 7. Experiment 1: a) trajectory of movement; b) distance between objects; c) ground
speed; d) heading angles</p>
      </sec>
    </sec>
    <sec id="sec-5">
      <title>Conclusions</title>
      <p>1. UAVs are widely used in different areas of human activity, and multi-UAV
performance has many advantages compared with the performance of an individual UAV.
Research institutions and groups are currently developing an algorithm for a group of
UAV autonomous control since manual control is not available.</p>
      <p>2. For multi-UAV formation control, the artificial potential field approach is used,
where UAVs are denoted as interacting dynamic objects influenced by attraction and
repulsion forces. The movement of each dynamic object is described by a system of
equations, with the direction of movement coinciding with thrust force angle projected
on each of axes.</p>
      <p>3. To check the potential field approach applicability, two simulations were
performed for 8 and 12 dynamic objects. The tasks were to form a group, avoid obstacles,
and continue movement in the given direction with no change in the shape of the group.
The results show that in this form the approach can be applied to a group formation and
multi-UAV flight control. All dynamic objects moved within the allowable range
determined by heading angle χi and ground speed Vi keeping within protection zones.
11. Xue, R., Cai, G.: Formation Flight Control of Multi-UAV System with Communication</p>
      <p>
        Constraints. Journal of Aerospace Technology and Management 8(
        <xref ref-type="bibr" rid="ref2">2</xref>
        ), 203-210 (2016).
12. Pavlova, S.V., Pavlov, V.V., Chepizhenko, V.I.: Virtual Einstein force fields in synergy
of navigation environment of difficult ergatic systems. Proceedings of the National Aviation
University 3, 15–27 (2012).
13. Chepizhenko, V.I.: Synthesis of artificial gravitational fields virtual meters for the
polyconflicts resolution in the aeronavigation environment. Proceedings of the National Aviation
University 2, 60–69 (2012).
      </p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Boucher</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          :
          <article-title>Domesticating the Drone: The Demilitarisation of Unmanned Aircraft for Civil Markets</article-title>
          .
          <source>Science and Engineering Ethics</source>
          <volume>21</volume>
          ,
          <fpage>1393</fpage>
          -
          <lpage>1412</lpage>
          (
          <year>2015</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Kuchar</surname>
            ,
            <given-names>J. K.</given-names>
          </string-name>
          , &amp;
          <string-name>
            <surname>Yang</surname>
            ,
            <given-names>L. C.</given-names>
          </string-name>
          :
          <article-title>A Review of Conflict Detection and Resolution Modeling Methods</article-title>
          .
          <source>IEEE Transactions on Intelligent Transportation Systems</source>
          <volume>1</volume>
          (
          <issue>4</issue>
          ),
          <fpage>179</fpage>
          -
          <lpage>189</lpage>
          (
          <year>2000</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Murray</surname>
            ,
            <given-names>R. M.</given-names>
          </string-name>
          :
          <article-title>Recent Research in Cooperative Control of Multivehicle Systems</article-title>
          .
          <source>Journal of Dynamic Systems, Measurement, and Control</source>
          <volume>129</volume>
          (
          <issue>5</issue>
          ),
          <fpage>571</fpage>
          -
          <lpage>598</lpage>
          (
          <year>2007</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Geser</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Muñoz</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          :
          <article-title>A Geometric Approach to Strategic Conflict Detection and Resolution</article-title>
          .
          <source>In: Proceedings of the 21st Digital Avionics Systems Conference</source>
          , vol.
          <volume>1</volume>
          , pp.
          <fpage>6B1</fpage>
          -
          <lpage>1</lpage>
          - 6B1-11, Irvine, CA, USA (
          <year>2002</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Park</surname>
          </string-name>
          , J.-W.,
          <string-name>
            <surname>Oh</surname>
            ,
            <given-names>H.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Tahk</surname>
          </string-name>
          , M.-J.:
          <article-title>UAV collision avoidance based on geometric approach</article-title>
          .
          <source>In: Proceedings of the 2008 SICE Annual Conference</source>
          , pp.
          <fpage>2122</fpage>
          -
          <lpage>2126</lpage>
          , Tokyo, (
          <year>2008</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Matsuno</surname>
            ,
            <given-names>Y.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Tsuchiya</surname>
            ,
            <given-names>T.</given-names>
          </string-name>
          :
          <article-title>Probabilistic Conflict Detection in the Presence of Uncertainty</article-title>
          . In: Electronic Navigation Research Institute (eds) Air
          <source>Traffic Management and Systems. Lecture Notes in Electrical Engineering</source>
          , vol
          <volume>290</volume>
          , pp.
          <fpage>17</fpage>
          -
          <lpage>33</lpage>
          (
          <year>2014</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Borrelli</surname>
            ,
            <given-names>F.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Subramanian</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Raghunathan</surname>
          </string-name>
          . A.,
          <string-name>
            <surname>Biegler</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          :
          <article-title>MILP and NLP techniques for centralized trajectory planning of multiple unmanned air vehicles</article-title>
          .
          <source>In: Proceedings American Control Conference</source>
          , pp.
          <fpage>5763</fpage>
          -
          <lpage>5768</lpage>
          , Minneapolis, USA (
          <year>2006</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>Chepizhenko</surname>
          </string-name>
          , V.:
          <article-title>Energy-potential method of dynamic objects polyconflicts guaranteed collision resolution</article-title>
          . Cybernetics and computer engineering 168,
          <fpage>80</fpage>
          -
          <lpage>87</lpage>
          (
          <year>2012</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <surname>Leonard</surname>
            ,
            <given-names>N.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Fiorelli</surname>
          </string-name>
          , E.:
          <article-title>Virtual leaders, artificial potentials and coordinated control of groups</article-title>
          .
          <source>In: Proceedings of the 40th IEEE Conference on Decision and Control</source>
          , vol.
          <volume>3</volume>
          , pp.
          <fpage>2968</fpage>
          -
          <lpage>2973</lpage>
          , Orlando, USA (
          <year>2001</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <surname>Nguyen</surname>
            ,
            <given-names>B.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Chuang</surname>
            ,
            <given-names>Y.-L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Tung</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Hsieh</surname>
          </string-name>
          , Ch.,
          <string-name>
            <surname>Jin</surname>
          </string-name>
          , Zh.,
          <string-name>
            <surname>Shi</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Marthaler</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Bertozzi</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Murray</surname>
            ,
            <given-names>R. M.:</given-names>
          </string-name>
          <article-title>Virtual attractive-repulsive potentials for cooperative control of second order dynamic vehicles on the Caltech MVWT</article-title>
          .
          <source>In: Proceedings of the American Control Conference</source>
          , vol.
          <volume>2</volume>
          , pp.
          <fpage>1084</fpage>
          -
          <lpage>1089</lpage>
          , Portland, USA (
          <year>2005</year>
          ).
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>