<!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>Parametric synthesis of a dynamic object control system with nonlinear characteristics</article-title>
      </title-group>
      <contrib-group>
        <aff id="aff0">
          <label>0</label>
          <institution>Kharkiv National Automobile and Highway University</institution>
          ,
          <addr-line>Str. Yaroslava Mudrogo 25, Kharkiv, 61002</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
      </contrib-group>
      <fpage>0000</fpage>
      <lpage>0002</lpage>
      <abstract>
        <p>The aim of the research is to solve the problem of using the method of factor experiment for the problem of parametric synthesis of a closed dynamic object control system. The method of factor experiment allows, on the basis of modeling the behavior of the system in a random external environment, to select its parameters that satisfy the requirements of the minimum of the regression functional. A quality functional that realize the requirements for a dynamic object control system and reduces the computational resource in the parametric synthesis of the system was proposed. The behavior of a closed control system was simulated considering the random external perturbations acting on the control object. The introduction of a nonlinear link with variable amplification factor widens the stability area of a closed dynamic object control system. The parametric synthesis of the closed dynamic object control system considering nonlinear characteristics on the basis of the factor experiment method was considered.</p>
      </abstract>
      <kwd-group>
        <kwd>Car</kwd>
        <kwd>Parametric Synthesis</kwd>
        <kwd>Quality Functional</kwd>
        <kwd>Transient Process</kwd>
        <kwd>Dynamic Object</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>Recently, the simulation methods become widespread use, in particular the method of
factor experiment. This method allows, on the basis of modeling the behavior of the
system in a random external environment, to select its parameters that satisfy the
requirements of the minimum of the regression functional. Additive quality functional
of the closed system can be used as a regression function. In this case, the behavior of
a closed control system is simulated into account of the random external perturbations
acting on the control object.</p>
      <p>Formulation of the problem
The aim of the research is to solve the problem of applying the method of factor
experiment to the problem of parametric synthesis of the dynamic object control system.
Due to continuous external disturbance Med t  , the position of the dynamic object
is constantly changing. Changing the position of a dynamic object Med t  is a
random function. In this situation, random functions are also functions  c t  and
 c t  (Fig. 1). The accuracy of stabilizing a dynamic object is higher, the smaller
the area under the dynamic process curves.
3</p>
      <p>Literature review
The methods of classical automatic control theory do not work for high-order
mathematical models. If the order of the differential equations that make up the
mathematical model of the control object is higher than five, then the problem of parametric
synthesis can be solved only by the methods of modern control theory, in particular,
the method of state-space, the methods of the theory of analytic design of optimal
regulators [1-3], minimax methods, methods of Liapunov functions, methods of
simulation modeling.</p>
      <p>Minimax methods and methods of Liapunov functions, which were applied by E.E.
Alexandrov to solve the problem of parametric synthesis of the control system give
satisfactory results [4-6], because they provide for the use of information in the
control algorithm only about those components of the vector of the state of the control
object, the measurement of which is not difficult. But in the process of parametric
synthesis of a digital control system, these methods involve the transition from
differential equations of the mathematical model of the control object to equations, which
introduces a certain error in the synthesis process.</p>
      <p>In addition, the high order of the differential equations describing the perturbed
motion of the control object results in considerable computational difficulties
associated with limited computer memory when applying the minimax method and the
Liapunov function method. However, if the objects contain non-analytic
nonlinearities, then the above synthesis methods are inactive. Therefore, in recent years, the
simulation methods become widespread use, in particular the method of factor
experiment [7, 8].
4</p>
      <p>Research Methodology
Therefore, the quantitative accuracy of the stabilization of the dynamic object relative
to the base position can be estimated by the following functional
where M – is the symbol of mathematical expectation.</p>
      <p>
        The quality of stabilization processes is estimated not only by the changing of
position of the dynamic object, but also by the angular velocity of motion of the dynamic
object relative to its base position [9, 10]. The quantitative characteristic of this
movement may be functional
It was an ideal case, when closed dynamic object control system solutions, both
functionalities (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) and (
        <xref ref-type="bibr" rid="ref2">2</xref>
        ) would be minimized. But this is not possible and the minimums
of functionals (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) and (
        <xref ref-type="bibr" rid="ref2">2</xref>
        ) correspond to different values of the parameters of the
controller k and k . Therefore, it is advisable to choose the parameters k and k in
condition of a minimum of additive functionality
      </p>
      <p>I   2 I
1 с
  2 I
2 с
 M
T 
 12с2 (t )   22с2 (t) dt  ,
 0 
where 1 and 2 - are the weighting factors to be selected.</p>
      <p>
        But, many experiments [11] proved that the obtained values of the varied
parameters, which give a minimum of functionals (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) and (
        <xref ref-type="bibr" rid="ref2">2</xref>
        ) are almost not different (within
5%). Considering this fact, as well as the fact that the output signal from the control
unit is limited and does not allow the dynamic object to pass excessive speed, it is
advisable to use functional (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) or its modifications for this particular system.
      </p>
      <p>In general, a functional that evaluates the quality of the stochastic system looks like
I  M
 T
  t
 t1</p>
      <p>
        
 с t  dt  ,


(
        <xref ref-type="bibr" rid="ref1">1</xref>
        )
(
        <xref ref-type="bibr" rid="ref2">2</xref>
        )
(
        <xref ref-type="bibr" rid="ref3">3</xref>
        )
(
        <xref ref-type="bibr" rid="ref4">4</xref>
        )
(
        <xref ref-type="bibr" rid="ref5">5</xref>
        )
I  M
T 
 f e t  , x t  , y t  , t  dt  ,
 0 
where f – is the function of error, input and output signals, as well as time. Using
different combinations of system and time variables, different quality estimates can be
obtained.
      </p>
      <p>
        To reduce the contribution of a significant initial error and to account for a future
error, it is more appropriate to use the form functionality
T
and when Med t   0 functional (
        <xref ref-type="bibr" rid="ref5">5</xref>
        ) is transformed into I   t с t  dt .
t1
Choosing as variables parameters the coefficients k and k* , and as the optimization
parameters, the functional (
        <xref ref-type="bibr" rid="ref5">5</xref>
        ), we use the factor experiment theory [7, 8] to find the
values of the variable parameters of the controller that give the least functional (
        <xref ref-type="bibr" rid="ref5">5</xref>
        ) for
the dynamic object control system.
      </p>
      <p>
        The points a, b correspond to the minimum of the objective function (
        <xref ref-type="bibr" rid="ref5">5</xref>
        ) without
and taking into account WН1  AН1  respectively. Parameters for (
        <xref ref-type="bibr" rid="ref5">5</xref>
        ): t1  0.25s,
T  5 s. For the case without taking into account WН1  AН1  – I  0.005958 ,
k*  199.5 , k*  16.3 . For the case considered WН1  AН1  – I  8.186339 ,
k*  323.2 , k*  21.0 .
      </p>
      <p>In Fig. 2 the transients processes of a closed dynamic object control system with
*
the obtained values of the variable parameters of the controller k and k* for the case
without taking into account WН1  AН1  are presented. As you can see, transients
processes are smooth without significant fluctuations. The amount of overshoot was
0.5 %.
k*  16.3 : control voltage U y t  (a), dynamic object speed с t (b) and for the value
of moving a dynamic object 100mm с t  (c) and the objective function I t  (d)
In Fig. 3 the transients processes of a closed dynamic object control system with the
obtained values of the variable parameters of the controller k* and k* in the case of
the given WН1  AН1  are presented. As you can see, transients processes are smooth
without significant fluctuations. The amount of overshoot was 0.1 %.</p>
      <p>A structural diagram of a dynamic object control system, in addition to a nonlinear
link with a transfer function WН1  AН1  and characteristic, of the limited (saturation)
zone type also contains a nonlinear link with a transfer function WН2  A
Н 2
 and
characteristic with variable amplification factor.</p>
      <p>In this case, this nonlinear link is introduced into the structural diagram in order to
improve the quality of transients processes and increase the reliability of the system as
a whole by extending the stability area.
k*  21.0 : control voltage U y t  (a), dynamic object speed c t (b) and for the value
of moving a dynamic object 100mm c t  (c) and the objective function I t  (d)
From the structural diagram we find the transfer function of the closed circuit with the
nonlinear link with the transfer function WН2  AН2 </p>
      <p>W3НН  s  </p>
      <p>W3Н  s  W
Н2
  k </p>
      <p>kgs
 s</p>
      <p> s
W4НН </p>
      <p>W4НН  s
1 W4НН  s</p>
      <p>We construct the stability domain of a closed dynamic object control system in the
plane of variable parameters of the controller k and k , taking into account a
nonlinear link such as a zone of limitation (saturation) and a nonlinear link with a
variable amplification factor.</p>
      <p>Ty A3 T21 s7  Ty A3 T 2 Ty A2 T21  A3 T21  s6 
  A2 T21 Ty A2 T 2 Ty A3  A3 T 2 Ty A1 T21  s5 
 Ty T21  A2 T 2 Ty A2  A3 Ty A1 T 2  A1 T21  s4 </p>
      <p> T21  A1 T 2 Ty A1  Ty T 2  A2  s3 
 T 2  A1  Ty  k k WН1  AН1  WН2  A
 kу k2 T21  s2 
1  k WН1  AН1  WН2  A</p>
      <p>Н2  k у k2  k 2 T 2  k kgs  s 
k k WН1  AН1  WН2  A</p>
      <p>
        From the transfer function of the closed system (
        <xref ref-type="bibr" rid="ref7">7</xref>
        ) we write the characteristic
equation (
        <xref ref-type="bibr" rid="ref8">8</xref>
        ). In characteristic equation (
        <xref ref-type="bibr" rid="ref8">8</xref>
        ) we make a substitution s  j , select the real
and imaginary parts and equal them to zero. Obtained algebraic equations we solve
referring to the parameters k and k :
(
        <xref ref-type="bibr" rid="ref7">7</xref>
        )
(
        <xref ref-type="bibr" rid="ref8">8</xref>
        )
k 

k 

Ty T21  A2 T 2  Ty A2  A3  Ty A1 T 2  A1 T21   4
      </p>
      <p>k WН1  AН1  WН2  AН2  kу k2  2T21 1
Ty A3 T 2  Ty A2 T21  A3 T21   6  T 2  A1  Ty   2</p>
      <p>k WН1  AН1  WН2  AН2  kу k2  2T21 1
T21  A1 T 2  Ty A1  Ty T 2  A2   2  Ty A3 T21  6</p>
      <p>k WН1  AН1  WН2  AН2  kу k2 kgs
 A2 T21  Ty A2 T 2  Ty A3  A3 T 2  Ty A1 T21   4
k WН1  AН1  WН2  AН2  k у k2 kgs

;
</p>
      <p>
        (
        <xref ref-type="bibr" rid="ref9">9</xref>
        )

(
        <xref ref-type="bibr" rid="ref10">10</xref>
        )

k k WН1  AН1  WН2  A
k WН1  AН1  WН2  A
Н2  k у k2 T 2 1
Н2  k у k2 kgs
.
      </p>
      <p>
        In Fig. 4 and 5 are shown the limits of stability area of the closed dynamic object
control systems in the plane of variable parameters of the controller k and k taking
without and with considering WН1  AН1  , WН2  A
 respectively.
into account a nonlinear variable with amplification factor, constructed using relations
(
        <xref ref-type="bibr" rid="ref9">9</xref>
        ) and (
        <xref ref-type="bibr" rid="ref10">10</xref>
        ), where a, b, c – are the points of minimum of the objective function (
        <xref ref-type="bibr" rid="ref5">5</xref>
        )
For the case with WН1  AН1  and WН2  A
 based on the method of factor
experiment we obtain: I  8.190428 , k*  607.1, k*  49.9 . Parameters for (
        <xref ref-type="bibr" rid="ref5">5</xref>
        ):
t1  0.25 s, T  5 s.
system (with kН2 Н3
stability of the linear system).
      </p>
      <p>As you can see from Fig. 4, the introduction of a nonlinear link [12 - 14] with
variable amplification factor widens the stability area of a closed dynamic object control
 k (curve 1), the stability area coincides completely with the
The optimal point c in this case moves to the area of high coefficients k and k .
The above helps to improve the reliability and accuracy of a closed dynamic object
control system.</p>
      <p>As you can see from Fig. 5, when A b  100 the equation
WН2  AН2   kН3 , and the stability area approaches to the linear.</p>
      <p>In Fig. 6 the transients processes of a closed dynamic object control system with
*
the obtained values of the variable parameters of the controller k and k*
considering WН1  AН1  and WН2  A
 are presented. The amount of overshoot was 0.1%.
Fig. 6. Transients processes of a closed dynamic object control system when k*  607.1,
k*  49.9 : control voltage Uy t  (a), dynamic object speed c t (b) and for the value
of moving a dynamic object 100mm c t (c) and the objective function I t  (d)
For a more positive result when using a variable structures need to use a more
complicated law of the changing the structure of the system, but it can reduce the
reliability of the system as a whole, due to the additional elements and connections.
5</p>
      <p>Conclusions</p>
      <p>The aim of the research to solve the problem of using the method of factor
experiment for the problem of parametric synthesis of a closed dynamic object control
system was considered. Quality functional that realize the requirements of a dynamic
object control system and reduce the calculating resource of parametric system
synthesis was proposed. Parametric synthesis of a closed dynamic object control system
with nonlinear characteristics on the basis of the factor experiment method was made
by using three practical examples.</p>
      <p>It was proved that the method of factor experiment allows, on the basis of
modeling the behavior of the system in a random external environment, to select its
parameters that satisfy the requirements of the minimum of the regression functional. The
behavior of a closed control system was simulated considering the random external
perturbations acting on the control object. It was defined that the introduction of a
nonlinear link with variable amplification factor widens the stability area of a closed
dynamic object control system.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Letov</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          :
          <article-title>Flight dynamics and control</article-title>
          . M,
          <string-name>
            <surname>Science</surname>
          </string-name>
          (
          <year>1969</year>
          )
          <article-title>(in Russian)</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Letov</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          :
          <article-title>Mathematical theory of control processes</article-title>
          . M,
          <string-name>
            <surname>Science</surname>
          </string-name>
          (
          <year>1981</year>
          )
          <article-title>(in Russian)</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Roitenberg</surname>
            ,
            <given-names>Y.</given-names>
          </string-name>
          :
          <article-title>Automatic control</article-title>
          . M,
          <string-name>
            <surname>Science</surname>
          </string-name>
          (
          <year>1978</year>
          )
          <article-title>(in Russian)</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Aleksandrov</surname>
            ,
            <given-names>E.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Bogaenko</surname>
            ,
            <given-names>I.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kuznetsov</surname>
            ,
            <given-names>B.</given-names>
          </string-name>
          :
          <article-title>Multichannel optimal control systems</article-title>
          . K,
          <string-name>
            <surname>Tekhnika</surname>
          </string-name>
          (
          <year>1995</year>
          )
          <article-title>(in Russian)</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Aleksandrov</surname>
            ,
            <given-names>E.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kostenko</surname>
            ,
            <given-names>Y.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kuznetsov</surname>
            ,
            <given-names>B.</given-names>
          </string-name>
          :
          <article-title>Optimization of multichannel control systems</article-title>
          . Kharkiv,
          <string-name>
            <surname>Osnova</surname>
          </string-name>
          (
          <year>1996</year>
          )
          <article-title>(in Russian)</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Aleksandrov</surname>
            ,
            <given-names>E.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Borisyuk</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kuznetsov</surname>
            ,
            <given-names>B.</given-names>
          </string-name>
          :
          <article-title>Parametric optimization of multichannel automatic control systems</article-title>
          . Kharkiv,
          <string-name>
            <surname>Osnova</surname>
          </string-name>
          (
          <year>1995</year>
          )
          <article-title>(in Russian)</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Adler</surname>
            ,
            <given-names>Y.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Markova</surname>
            ,
            <given-names>E.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Granovskii</surname>
            ,
            <given-names>Y.</given-names>
          </string-name>
          :
          <article-title>Planning of the experiment in search of optimal conditions</article-title>
          . M,
          <string-name>
            <surname>Science</surname>
          </string-name>
          (
          <year>1976</year>
          )
          <article-title>(in Russian)</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>Alexandrov</surname>
            ,
            <given-names>E.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Nikonov</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Skvorchevsky</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          :
          <article-title>Structural and parametric synthesis of the system of automatic control of the car brakes</article-title>
          .
          <source>Energy and resource saving</source>
          , vol.
          <volume>6</volume>
          , pp.
          <fpage>30</fpage>
          -
          <lpage>39</lpage>
          (
          <year>2009</year>
          )
          <article-title>(in Ukrainian)</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <surname>Xu</surname>
            ,
            <given-names>G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ma</surname>
            ,
            <given-names>Z</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Lu</surname>
            ,
            <given-names>F.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Hou</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          :
          <source>Kinematic Analysis of Hydraulic Excavator Working Device Based on DH Method. In: International Conference on Applied Mechanics, Mechanical and Materials Engineering</source>
          , pp.
          <volume>8</volume>
          (
          <year>2016</year>
          ). doi:
          <volume>10</volume>
          .12783/dtmse/ammme2016/6857
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <surname>Gurko</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kyrychenko</surname>
            ,
            <given-names>I.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Yaryzhko</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          :
          <article-title>Trajectories Planning and Simulation of a Backhoe Manipulator Movement</article-title>
          . CMIS, pp.
          <fpage>771</fpage>
          -
          <lpage>785</lpage>
          (
          <year>2019</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          11.
          <string-name>
            <surname>Xu</surname>
            ,
            <given-names>J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Yoon</surname>
            ,
            <given-names>H. S.:</given-names>
          </string-name>
          <article-title>A Review on Mechanical and Hydraulic System Modeling of Excavator Manipulator System</article-title>
          .
          <source>Journal of Construction Engineering</source>
          <year>2016</year>
          ,
          <volume>9409370</volume>
          (
          <year>2016</year>
          ). doi:
          <volume>10</volume>
          .1155/
          <year>2016</year>
          /9409370.
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          12.
          <string-name>
            <surname>Uspensky</surname>
            ,
            <given-names>B.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Avramov</surname>
            ,
            <given-names>K.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Liubarskyi</surname>
            ,
            <given-names>B.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Andrieiev</surname>
            ,
            <given-names>Y.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Nikonov</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          :
          <article-title>Nonlinear torsional vibrations of electromechanical coupling of diesel engine gear system and electric generator</article-title>
          .
          <source>Journal of Sound and Vibration</source>
          , vol.
          <volume>460</volume>
          ,
          <issue>114877</issue>
          (
          <year>2019</year>
          ) doi:10.1016/j.jsv.
          <year>2019</year>
          .114877
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          13.
          <string-name>
            <surname>Uspensky</surname>
            ,
            <given-names>B.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Avramov</surname>
            ,
            <given-names>K.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Nikonov</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          :
          <article-title>Nonlinear modes of piecewise linear systems forced vibrations close to superharmonic resonances</article-title>
          .
          <source>Proceedings of the Institution of Mechanical Engineers</source>
          ,
          <string-name>
            <surname>Part</surname>
            <given-names>C</given-names>
          </string-name>
          :
          <article-title>Journal of Mechanical Engineering Science</article-title>
          , vol.
          <volume>233</volume>
          ,
          <string-name>
            <surname>Issue</surname>
          </string-name>
          23-
          <issue>24</issue>
          , pp.
          <fpage>7489</fpage>
          -
          <lpage>7497</lpage>
          (
          <year>2019</year>
          ) doi: 10.1177/0954406219869967
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          14.
          <string-name>
            <surname>Gu</surname>
            ,
            <given-names>J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ma</surname>
            ,
            <given-names>X. D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ni</surname>
            ,
            <given-names>J. F.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Sun</surname>
            ,
            <given-names>L. N.</given-names>
          </string-name>
          :
          <article-title>Linear and nonlinear control of a robotic excavator</article-title>
          .
          <source>J. Cent. South Univ. 19</source>
          , pp.
          <fpage>1823</fpage>
          −
          <lpage>1831</lpage>
          (
          <year>2012</year>
          ). doi:
          <volume>10</volume>
          .1007/s11771-012-1215-y
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>