<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing with OASIS Tables v3.0 20080202//EN" "journalpub-oasis3.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:oasis="http://docs.oasis-open.org/ns/oasis-exchange/table" xml:lang="en" dtd-version="3.0" article-type="research-article">
  <front>
    <journal-meta><journal-id journal-id-type="publisher">MS</journal-id><journal-title-group>
    <journal-title>Mechanical Sciences</journal-title>
    <abbrev-journal-title abbrev-type="publisher">MS</abbrev-journal-title><abbrev-journal-title abbrev-type="nlm-ta">Mech. Sci.</abbrev-journal-title>
  </journal-title-group><issn pub-type="epub">2191-916X</issn><publisher>
    <publisher-name>Copernicus Publications</publisher-name>
    <publisher-loc>Göttingen, Germany</publisher-loc>
  </publisher></journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.5194/ms-13-751-2022</article-id><title-group><article-title>Dynamic and sliding mode control of <?xmltex \hack{\break}?> space netted pocket system capturing and <?xmltex \hack{\break}?>
attitude maneuver non-cooperative target</article-title><alt-title>Dynamic and sliding mode control of space netted pocket system capturing</alt-title>
      </title-group><?xmltex \runningtitle{Dynamic and sliding mode control of space netted pocket system capturing}?><?xmltex \runningauthor{C. Tang et al.}?>
      <contrib-group>
        <contrib contrib-type="author" corresp="no">
          <name><surname>Tang</surname><given-names>Chao</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="yes">
          <name><surname>Huang</surname><given-names>Zhuoran</given-names></name>
          <email>huangzhuoranhzr@163.com</email>
        </contrib>
        <contrib contrib-type="author" corresp="no">
          <name><surname>Wei</surname><given-names>Cheng</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no">
          <name><surname>Zhao</surname><given-names>Yang</given-names></name>
          
        </contrib>
        <aff id="aff1"><institution>Harbin Institute of Technology, Harbin 150001, China</institution>
        </aff>
      </contrib-group>
      <author-notes><corresp id="corr1">Zhuoran Huang (huangzhuoranhzr@163.com)</corresp></author-notes><pub-date><day>30</day><month>August</month><year>2022</year></pub-date>
      
      <volume>13</volume>
      <issue>2</issue>
      <fpage>751</fpage><lpage>760</lpage>
      <history>
        <date date-type="received"><day>4</day><month>May</month><year>2022</year></date>
           <date date-type="rev-recd"><day>30</day><month>June</month><year>2022</year></date>
           <date date-type="accepted"><day>24</day><month>July</month><year>2022</year></date>
      </history>
      <permissions>
        <copyright-statement>Copyright: © 2022 Chao Tang et al.</copyright-statement>
        <copyright-year>2022</copyright-year>
      <license license-type="open-access"><license-p>This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link></license-p></license></permissions><self-uri xlink:href="https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022.html">This article is available from https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022.html</self-uri><self-uri xlink:href="https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022.pdf">The full text article is available as a PDF file from https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022.pdf</self-uri>
      <abstract><title>Abstract</title>

      <p id="d1e106">Similar to a space flying net, the capture field of the space
netted pocket system is large and it can be applied to capture space
non-cooperative targets flexibly. To maintain the stability of the space
netted pocket system, eight inflatable rods are used as the supporting
structure of the net surface. In this paper, a space netted pocket system is
designed and modeled. Based on ANCF (absolute nodal coordinate formulation),
a dynamic model of the complex space rope net system is established, and
then an accurate model of closing rope considering the variable length is
derived by introducing mass flow element. A double closed-loop sliding
control method is designed to maintain the stable attitude of the service
spacecraft. An extended observer is applied to estimate and compensate for
the disturbances due to the uncertainty of the contact and flexibility in
the system. Finally, the dynamic model and control method is verified
through the simulation of the virtual prototype. Results show that the
service spacecraft can maintain the attitude stability during target
captured process and can track the desired angle during attitude maneuver.
The flexible deformation and collision cause great disturbance to the
service spacecraft, and the extended observer can improve the control
accuracy from 10<inline-formula><mml:math id="M1" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> to 10<inline-formula><mml:math id="M2" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">4</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>.</p>
  </abstract>
    </article-meta>
  </front>
<body>
      

<sec id="Ch1.S1" sec-type="intro">
  <label>1</label><title>Introduction</title>
      <p id="d1e142">In recent years, with the increase of human space activities, threats to
space activities from space debris are
worsening (Barmin et al., 2014; Chen, 2011).
The contact capture method is a simple and feasible way to capture space debris. There are mainly two types of contact capture: rigid capture and
flexible capture (Nishida and Kawamoto,
2011; Shan et al., 2016). Like space manipulator and capture claw, the
rigid capture methods will inevitably cause collision between the capture
target and operation platform. Therefore, the rigid capture methods are more
complicated and require high control accuracy. The space net as a typical
flexible acquisition method is more lightweight and simple
(Zhang et al., 2017). The capture domain of the space net is
larger, which can reduce the control accuracy requirements and reduce the
cost. Therefore, the space net is an active space debris removal method with
broad application prospects.</p>
      <p id="d1e145">The space fly net launches a flexible rope net to cover the space debris in
a large envelope and drags space debris to the atmosphere burned or grave
orbit. In recent years, a lot of theoretical and experimental researches
has been carried out around the space fly net. In terms of large-scale
experiments, RemoveDEBRIS project has now completed the first space
rope net capture test in
orbit (Forshaw et
al., 2016, 2017) and e.Deorbit project has carried out parabolic flight
test (Biesbroek et al., 2017). At the same time, many
theories about space net dynamics modeling, space net deployment control,
space net capture performance analysis and space net structure optimization
design are deeply
researched (Bonnal
et al., 2013; Ming et al., 2017; Xu et al., 2019). However, it is difficult
for the fly net system to maintain the net shape fully deployed for a long
time. And the capture success rate will be reduced due to the space fly net
hardly maintaining the net shape for a long time. Some scholars have
proposed the motorized fly net to control the net shape
(Huang et al., 2015; Meng et
al., 2017), but multiple motorized devices around the net will greatly
increase the complexity and cost of the system. The space rope netted pocket
capture mechanism described in this paper can inherit the advantages of
lightness and simplicity. Furthermore, the inflatable rods in the system can
maintain the net shape for a long time, which can improve the success rate
of capture. The rope net surface of the netted pocket capture mechanism in
Fig. 1 is supported in an umbrella-like shape by
eight flexible rods. First, the service spacecraft approaches the captured
target until the target enters the net. Second, the closing devices at the
end of inflatable rods are tightened to complete the closing. Finally, with
attitude maneuver of the service spacecraft, the space netted pocket drags
the captured target to change orbit.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F1"><?xmltex \currentcnt{1}?><?xmltex \def\figurename{Figure}?><label>Figure 1</label><caption><p id="d1e150">Space netted pocket capture mechanism.</p></caption>
        <?xmltex \igopts{width=184.942913pt}?><graphic xlink:href="https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022-f01.png"/>

      </fig>

      <p id="d1e160">Due to the strong nonlinearity dynamics of the net and the large uncertainty
of the target, the precise control of the rope net spacecraft is more
difficult. The traditional PD (proportional differential) controller can hardly meet the accuracy
requirements of the mission. For the stability control of tethered
satellites, Huang et al. considered the effect of the flexibility of rope,
and designed an adaptive control method to achieve attitude stability
control of combinations with unknown parameters
(Huang et al., 2016); Wei et al. (2019) designed
a spacecraft attitude stability ADRC (active disturbance rejection control) controller which can estimate and
compensate uncertainties of system parameters in real
time (Wei et al., 2019).</p>
      <p id="d1e163">However, the current research generally only focuses on the attitude control
during capture or maneuver, and the process of capture and maneuver is not
continuous. Moreover, the dynamic model and simulation prototype for
controller verification is imperfect, which cannot simulate the nonlinear
properties and contact collision of the rope accurately. So, it is untrue to
verify controllers using the imperfect dynamic model for the space netted
pocket system.</p>
      <p id="d1e166">To solve these problems, this paper designs an extended observer to compensate
the uncertainty and disturbance of the model based on the principle of
active disturbance rejection control (ADRC). Combined with the robust
sliding mode control method, the attitude controller of the service
spacecraft is constructed. At the same time, a complex dynamic model of the
net pocket system is established based on the dynamic theories such as
variable flexible cable, ANCF (absolute nodal coordinate formulation) and collision theory. The closed-loop
simulation of dynamics and control can be completed through the virtual
simulation prototype.</p>
</sec>
<sec id="Ch1.S2">
  <label>2</label><title>Dynamics modeling of space netted pocket</title>
<sec id="Ch1.S2.SS1">
  <label>2.1</label><title>Netted pocket system composition and size design</title>
      <p id="d1e184">As shown in Fig. 2, the netted pocket system
consists of eight flexible support inflatable rods, which are connected by
rope net. The inflatable rods are evenly distributed in a positive octagonal
shape, with a single rod diameter of 0.1 m. The rope-retracting mechanism at
the end of each rod pulls the rope to close the net. The size and braided
shape of the single net piece between support rods are designed as shown in
Fig. 3. The upper part of the square rope net is
<inline-formula><mml:math id="M3" display="inline"><mml:mrow><mml:mn mathvariant="normal">14</mml:mn><mml:mo>×</mml:mo><mml:mn mathvariant="normal">12</mml:mn></mml:mrow></mml:math></inline-formula> (row <inline-formula><mml:math id="M4" display="inline"><mml:mo>×</mml:mo></mml:math></inline-formula> column) configuration and the lower part of the
triangular rope net is <inline-formula><mml:math id="M5" display="inline"><mml:mrow><mml:mn mathvariant="normal">10</mml:mn><mml:mo>×</mml:mo><mml:mn mathvariant="normal">12</mml:mn></mml:mrow></mml:math></inline-formula> (row <inline-formula><mml:math id="M6" display="inline"><mml:mo>×</mml:mo></mml:math></inline-formula> column) configuration.
The maximum diameter of the net pocket is 15.12 m and the maximum depth is
12 m.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F2"><?xmltex \currentcnt{2}?><?xmltex \def\figurename{Figure}?><label>Figure 2</label><caption><p id="d1e227">Space netted pocket system composition.</p></caption>
          <?xmltex \igopts{width=199.169291pt}?><graphic xlink:href="https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022-f02.png"/>

        </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F3"><?xmltex \currentcnt{3}?><?xmltex \def\figurename{Figure}?><label>Figure 3</label><caption><p id="d1e238">Space netted pocket system size.</p></caption>
          <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022-f03.png"/>

        </fig>

</sec>
<sec id="Ch1.S2.SS2">
  <label>2.2</label><title>Dynamic model of the space netted pocket</title>
<sec id="Ch1.S2.SS2.SSS1">
  <label>2.2.1</label><title>Dynamic model of the closing ropes based on variable flexible ropes</title>
      <p id="d1e262">A rope retractor is installed at the end of each inflatable rod, and the
diameter of the opened net is adjusted by controlling the length of the
closing rope. Due to the rope retractor, the system structure is complex,
and it is difficult to simulate the ropes inside and outside the retractor
at the same time. In this case, the simulation takes a long time, and the
rope inside the retractor has little impact on the rope net capture system.
Therefore, in order to simplify the dynamic model of the rope net capture
device and improve the simulation speed, this paper establishes the dynamic
model of closing ropes based on the ANCF flexible cable dynamics model.
Using <inline-formula><mml:math id="M7" display="inline"><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> as the mass flow velocity at the element boundary node,
the variable length flexible cable element model of the closing rope is
established in Fig. 4.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F4"><?xmltex \currentcnt{4}?><?xmltex \def\figurename{Figure}?><label>Figure 4</label><caption><p id="d1e278">Model of closing rope.</p></caption>
            <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022-f04.png"/>

          </fig>

      <p id="d1e287">Based on ANCF, the generalized coordinates of flexible cable element <inline-formula><mml:math id="M8" display="inline"><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold-italic">q</mml:mi></mml:mrow></mml:math></inline-formula> can
be expressed as
              <disp-formula id="Ch1.E1" content-type="numbered"><label>1</label><mml:math id="M9" display="block"><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold-italic">q</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mfenced close="]" open="["><mml:mtable class="array" columnalign="center left"><mml:mtr><mml:mtd><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msub><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msub><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:msup><mml:mfenced close="]" open="["><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mtable class="array" columnalign="center left left left"><mml:mtr><mml:mtd><mml:mrow><mml:mo>(</mml:mo><mml:mn mathvariant="normal">0</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msubsup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mi>l</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msubsup><mml:mo>(</mml:mo><mml:mn mathvariant="normal">0</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mo>(</mml:mo><mml:mi>L</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msubsup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mi>l</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msubsup><mml:mo>(</mml:mo><mml:mi>L</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
            where <inline-formula><mml:math id="M10" display="inline"><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mo>(</mml:mo><mml:mn mathvariant="normal">0</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M11" display="inline"><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mo>(</mml:mo><mml:mi>L</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> are the position vector of nodes at
both ends of the element; and <inline-formula><mml:math id="M12" display="inline"><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msubsup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mi>l</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msubsup><mml:mo>(</mml:mo><mml:mn mathvariant="normal">0</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M13" display="inline"><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msubsup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mi>l</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msubsup><mml:mo>(</mml:mo><mml:mi>L</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> are the
gradient vector nodes at both ends of the element. And the position vector
<inline-formula><mml:math id="M14" display="inline"><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mo>(</mml:mo><mml:mi>l</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> of any node on the central axis of the element can be
written as

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M15" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E2"><mml:mtd><mml:mtext>2</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mo>(</mml:mo><mml:mi>l</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mi mathvariant="bold">S</mml:mi><mml:mfenced open="(" close=")"><mml:mi>l</mml:mi></mml:mfenced><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold-italic">q</mml:mi></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E3"><mml:mtd><mml:mtext>3</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:mi mathvariant="bold">S</mml:mi><mml:mfenced open="(" close=")"><mml:mi>l</mml:mi></mml:mfenced><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mfenced open="[" close="]"><mml:mtable class="array" columnalign="center center center center"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mi>I</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:mi>I</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msub><mml:mi>I</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mn mathvariant="normal">4</mml:mn></mml:msub><mml:mi>I</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mfenced open="{" close=""><mml:mtable class="array" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mn mathvariant="normal">1</mml:mn><mml:mn mathvariant="normal">4</mml:mn></mml:mfrac></mml:mstyle><mml:msup><mml:mfenced open="(" close=")"><mml:mrow><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mfenced open="(" close=")"><mml:mrow><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo>+</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:mfenced></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:mi>L</mml:mi><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mn mathvariant="normal">8</mml:mn></mml:mfrac></mml:mstyle><mml:msup><mml:mfenced open="(" close=")"><mml:mrow><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mfenced open="(" close=")"><mml:mrow><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo>+</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mfenced></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mn mathvariant="normal">1</mml:mn><mml:mn mathvariant="normal">4</mml:mn></mml:mfrac></mml:mstyle><mml:msup><mml:mfenced open="(" close=")"><mml:mrow><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo>+</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mfenced close=")" open="("><mml:mrow><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo>+</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:mfenced></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mn mathvariant="normal">4</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:mi>L</mml:mi><mml:mfenced open="(" close=")"><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mn mathvariant="normal">8</mml:mn></mml:mfrac></mml:mstyle><mml:msup><mml:mfenced close=")" open="("><mml:mrow><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo>+</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mfenced open="(" close=")"><mml:mrow><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mfenced></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

              where <inline-formula><mml:math id="M16" display="inline"><mml:mrow><mml:mi mathvariant="bold">S</mml:mi><mml:mo>(</mml:mo><mml:mi>l</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> is the shape function; <inline-formula><mml:math id="M17" display="inline"><mml:mi mathvariant="italic">ξ</mml:mi></mml:math></inline-formula> can be written as follows:
              <disp-formula id="Ch1.E4" content-type="numbered"><label>4</label><mml:math id="M18" display="block"><mml:mrow><mml:mi mathvariant="italic">ξ</mml:mi><mml:mfenced close=")" open="("><mml:mrow><mml:mi>l</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:mfenced><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mn mathvariant="normal">2</mml:mn><mml:mi>l</mml:mi><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced><mml:mo>-</mml:mo><mml:mi>L</mml:mi><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mfenced open="(" close=")"><mml:mi>t</mml:mi></mml:mfenced></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
            The length of the rope element <inline-formula><mml:math id="M19" display="inline"><mml:mrow><mml:mi>L</mml:mi><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> and the position of the particle
<inline-formula><mml:math id="M20" display="inline"><mml:mrow><mml:mi>l</mml:mi><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> are time variables, so the variable length element shape function
<inline-formula><mml:math id="M21" display="inline"><mml:mrow><mml:mi mathvariant="bold">S</mml:mi><mml:mfenced open="(" close=")"><mml:mrow><mml:mi>l</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:mfenced></mml:mrow></mml:math></inline-formula> is a function of time. Then the velocity and
acceleration of any point are expressed as
<?xmltex \hack{\newpage}?>

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M22" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E5"><mml:mtd><mml:mtext>5</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mover accent="true"><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mfenced close=")" open="("><mml:mrow><mml:mi>l</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:mfenced><mml:mo>=</mml:mo><mml:mi mathvariant="bold">S</mml:mi><mml:mfenced open="(" close=")"><mml:mi>l</mml:mi></mml:mfenced><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mover accent="true"><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced><mml:mo>+</mml:mo><mml:mover accent="true"><mml:mi mathvariant="bold">S</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mfenced close=")" open="("><mml:mi>l</mml:mi></mml:mfenced><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E6"><mml:mtd><mml:mtext>6</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mover accent="true"><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover><mml:mfenced open="(" close=")"><mml:mrow><mml:mi>l</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:mfenced><mml:mo>=</mml:mo><mml:mi mathvariant="bold">S</mml:mi><mml:mfenced open="(" close=")"><mml:mi>l</mml:mi></mml:mfenced><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mover accent="true"><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover><mml:mfenced open="(" close=")"><mml:mi>t</mml:mi></mml:mfenced><mml:mo>+</mml:mo><mml:mn mathvariant="normal">2</mml:mn><mml:mover accent="true"><mml:mi mathvariant="bold">S</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mfenced close=")" open="("><mml:mi>l</mml:mi></mml:mfenced><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mover accent="true"><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mfenced open="(" close=")"><mml:mi>t</mml:mi></mml:mfenced><mml:mo>+</mml:mo><mml:mover accent="true"><mml:mi mathvariant="bold">S</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover><mml:mfenced open="(" close=")"><mml:mi>l</mml:mi></mml:mfenced><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E7"><mml:mtd><mml:mtext>7</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mtable class="array" columnalign="center center"><mml:mtr><mml:mtd><mml:mrow><mml:mover accent="true"><mml:mi mathvariant="bold">S</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:mi>d</mml:mi><mml:mi mathvariant="bold">S</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mover accent="true"><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mover accent="true"><mml:mi mathvariant="bold">S</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:msup><mml:mi>d</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mi mathvariant="bold">S</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:msup><mml:mi mathvariant="italic">ξ</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mstyle><mml:msup><mml:mover accent="true"><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:mi>d</mml:mi><mml:mi mathvariant="bold">S</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mover accent="true"><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

              In the case of varying element length, the Lagrange method is no longer
applicable with the assumption of the mass flow at the boundary. The
dynamical equations are derived using D'Alembert's principle, which states
that the sum of virtual work done by inertial force and acting force on any
virtual displacement is 0.
              <disp-formula id="Ch1.E8" content-type="numbered"><label>8</label><mml:math id="M23" display="block"><mml:mrow><mml:msub><mml:mo>∑</mml:mo><mml:mi>i</mml:mi></mml:msub><mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>m</mml:mi><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mo>⋅</mml:mo><mml:mi mathvariant="italic">δ</mml:mi><mml:msub><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
            <inline-formula><mml:math id="M24" display="inline"><mml:mrow><mml:mi mathvariant="italic">δ</mml:mi><mml:msub><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the particle <inline-formula><mml:math id="M25" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula> in any moment to meet the constraints of
the virtual displacement. The force on the rope element can be decomposed
into elastic force <inline-formula><mml:math id="M26" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi mathvariant="normal">E</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and external force <inline-formula><mml:math id="M27" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi mathvariant="normal">f</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, then the dynamics
equation of the cable element can be written as
              <disp-formula id="Ch1.E9" content-type="numbered"><label>9</label><mml:math id="M28" display="block"><mml:mrow><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mn mathvariant="normal">0</mml:mn><mml:mfenced open="(" close=")"><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mfenced open="(" close=")"><mml:mi>t</mml:mi></mml:mfenced></mml:mrow></mml:munderover><mml:mi mathvariant="italic">δ</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi mathvariant="normal">f</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi mathvariant="normal">E</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ρ</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mover accent="true"><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>l</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
            Considering only the axial and bending deformation of the cable, the
expression of the virtual work of each part can be written as

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M29" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E10"><mml:mtd><mml:mtext>10</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mtable class="split" rowspacing="0.2ex" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mn mathvariant="normal">0</mml:mn><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mfenced open="(" close=")"><mml:mi>t</mml:mi></mml:mfenced></mml:mrow></mml:munderover><mml:mi mathvariant="italic">δ</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi mathvariant="normal">f</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>l</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mi mathvariant="italic">δ</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msup><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mn mathvariant="normal">0</mml:mn><mml:mfenced open="(" close=")"><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced></mml:mrow></mml:munderover><mml:msup><mml:mi mathvariant="bold">S</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mfenced close=")" open="("><mml:mi>l</mml:mi></mml:mfenced><mml:mi mathvariant="bold-italic">f</mml:mi><mml:mfenced close=")" open="("><mml:mrow><mml:mi>l</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>l</mml:mi></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mi mathvariant="italic">δ</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msup><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mi>l</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mn mathvariant="normal">1</mml:mn></mml:munderover><mml:msup><mml:mi mathvariant="bold">S</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mfenced open="(" close=")"><mml:mi mathvariant="italic">ξ</mml:mi></mml:mfenced><mml:mi mathvariant="bold-italic">f</mml:mi><mml:mfenced close=")" open="("><mml:mrow><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E11"><mml:mtd><mml:mtext>11</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mtable class="split" rowspacing="0.2ex" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mn mathvariant="normal">0</mml:mn><mml:mfenced open="(" close=")"><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced></mml:mrow></mml:munderover><mml:mi mathvariant="italic">δ</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi mathvariant="normal">E</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>l</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mi mathvariant="italic">δ</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msup><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mi>l</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mn mathvariant="normal">1</mml:mn></mml:munderover><mml:mo mathsize="2.0em">(</mml:mo><mml:msup><mml:mfenced close=")" open="("><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mo>∂</mml:mo><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold-italic">q</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mfenced><mml:mi mathvariant="normal">T</mml:mi></mml:msup></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>⋅</mml:mo><mml:mi>E</mml:mi><mml:mi>A</mml:mi><mml:mfenced close=")" open="("><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:mi>c</mml:mi><mml:msub><mml:mover accent="true"><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="italic">ε</mml:mi></mml:mrow><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mn mathvariant="normal">0</mml:mn></mml:msub></mml:mrow></mml:mfenced><mml:mo>+</mml:mo><mml:msup><mml:mfenced open="(" close=")"><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="italic">κ</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold-italic">q</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mfenced><mml:mi mathvariant="normal">T</mml:mi></mml:msup></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>⋅</mml:mo><mml:mi>E</mml:mi><mml:mi>J</mml:mi><mml:mfenced close=")" open="("><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="italic">κ</mml:mi><mml:mo>+</mml:mo><mml:mi>c</mml:mi><mml:mover accent="true"><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="italic">κ</mml:mi></mml:mrow><mml:mo mathvariant="normal">˙</mml:mo></mml:mover></mml:mrow></mml:mfenced><mml:mo mathsize="2.0em">)</mml:mo><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E12"><mml:mtd><mml:mtext>12</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mn mathvariant="normal">0</mml:mn><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mfenced open="(" close=")"><mml:mi>t</mml:mi></mml:mfenced></mml:mrow></mml:munderover><mml:msup><mml:mfenced open="(" close=")"><mml:mrow><mml:mi mathvariant="italic">δ</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold-italic">r</mml:mi></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mfenced close=")" open="("><mml:mrow><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ρ</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi>A</mml:mi><mml:mover accent="true"><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>x</mml:mi><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mi mathvariant="italic">δ</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msup><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mi>l</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mn mathvariant="normal">1</mml:mn></mml:munderover><mml:mi mathvariant="italic">ρ</mml:mi><mml:mi>A</mml:mi><mml:msup><mml:mi mathvariant="bold">S</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mover accent="true"><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

              where <inline-formula><mml:math id="M30" display="inline"><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mn mathvariant="normal">0</mml:mn><mml:mi>j</mml:mi></mml:msubsup><mml:mi mathvariant="italic">κ</mml:mi></mml:mrow></mml:math></inline-formula> are axial strain and curvature, and <inline-formula><mml:math id="M31" display="inline"><mml:mi>E</mml:mi></mml:math></inline-formula> is
elasticity modulus. <inline-formula><mml:math id="M32" display="inline"><mml:mi>A</mml:mi></mml:math></inline-formula> is the cross-sectional area and <inline-formula><mml:math id="M33" display="inline"><mml:mi>J</mml:mi></mml:math></inline-formula> is the moment
of inertia of the flexible cable element. <inline-formula><mml:math id="M34" display="inline"><mml:mi>c</mml:mi></mml:math></inline-formula> is the damping coefficient and
<inline-formula><mml:math id="M35" display="inline"><mml:mrow><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mi>l</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:mi>L</mml:mi><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle></mml:mrow></mml:math></inline-formula>.</p>
      <p id="d1e1888">Combining Eqs. (10), (11) and (12), the dynamic equation of the variable
flexible cable element can be written as
              <disp-formula id="Ch1.E13" content-type="numbered"><label>13</label><mml:math id="M36" display="block"><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold">M</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mover accent="true"><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover><mml:mo>+</mml:mo><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msubsup><mml:mi mathvariant="bold">M</mml:mi><mml:mi>v</mml:mi><mml:mi>j</mml:mi></mml:msubsup><mml:mover accent="true"><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mo>+</mml:mo><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msubsup><mml:mi mathvariant="bold">M</mml:mi><mml:mi>q</mml:mi><mml:mi>j</mml:mi></mml:msubsup><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mo>+</mml:mo><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold">Q</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="bold">0</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
            where

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M37" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E14"><mml:mtd><mml:mtext>14</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mtable class="split" rowspacing="0.2ex" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold">M</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi>L</mml:mi><mml:mfenced open="(" close=")"><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mn mathvariant="normal">1</mml:mn></mml:munderover><mml:mi mathvariant="italic">ρ</mml:mi><mml:msup><mml:mi mathvariant="bold">S</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mi mathvariant="bold">S</mml:mi><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msub><mml:mi mathvariant="bold">M</mml:mi><mml:mi>v</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi>L</mml:mi><mml:mfenced open="(" close=")"><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mn mathvariant="normal">1</mml:mn></mml:munderover><mml:mi mathvariant="italic">ρ</mml:mi><mml:msup><mml:mi mathvariant="bold">S</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mover accent="true"><mml:mi mathvariant="bold">S</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msub><mml:mi mathvariant="bold">M</mml:mi><mml:mi>q</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi>l</mml:mi><mml:mfenced open="(" close=")"><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mn mathvariant="normal">1</mml:mn></mml:munderover><mml:mi mathvariant="italic">ρ</mml:mi><mml:msup><mml:mi mathvariant="bold">S</mml:mi><mml:mi mathvariant="bold">T</mml:mi></mml:msup><mml:mover accent="true"><mml:mi mathvariant="bold">S</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E15"><mml:mtd><mml:mtext>15</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold">Q</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msub><mml:mi mathvariant="bold">Q</mml:mi><mml:mi>e</mml:mi></mml:msub><mml:msup><mml:mo>+</mml:mo><mml:mi>j</mml:mi></mml:msup><mml:msubsup><mml:mi mathvariant="bold">Q</mml:mi><mml:mi>f</mml:mi><mml:mi>j</mml:mi></mml:msubsup><mml:msub><mml:mi mathvariant="bold">Q</mml:mi><mml:mi>e</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi>L</mml:mi><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mn mathvariant="normal">1</mml:mn></mml:munderover><mml:mo mathsize="2.0em">(</mml:mo><mml:msup><mml:mfenced close=")" open="("><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="bold-italic">q</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mfenced><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mi>E</mml:mi><mml:mi>A</mml:mi><mml:mfenced close=")" open="("><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:mi>c</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">ε</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mn mathvariant="normal">0</mml:mn></mml:msub></mml:mrow></mml:mfenced></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>+</mml:mo><mml:msup><mml:mfenced close=")" open="("><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="italic">κ</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:mi mathvariant="bold-italic">q</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mfenced><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mi>E</mml:mi><mml:mi>J</mml:mi><mml:mfenced close=")" open="("><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mi mathvariant="italic">κ</mml:mi><mml:mo>+</mml:mo><mml:mi>c</mml:mi><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:mover accent="true"><mml:mi mathvariant="italic">κ</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover></mml:mrow></mml:mfenced><mml:mo mathsize="2.0em">)</mml:mo><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E16"><mml:mtd><mml:mtext>16</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:msup><mml:mi/><mml:mi>j</mml:mi></mml:msup><mml:msub><mml:mi mathvariant="bold">Q</mml:mi><mml:mi>f</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi>L</mml:mi><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced></mml:mrow><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mn mathvariant="normal">1</mml:mn></mml:munderover><mml:msup><mml:mi mathvariant="bold">S</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msup><mml:mi mathvariant="bold-italic">f</mml:mi><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="italic">ξ</mml:mi><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

              Considering a long cable with <inline-formula><mml:math id="M38" display="inline"><mml:mrow><mml:mi>N</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> flexible elements, only the lengths of the
flexible elements at the ends can be varied, while the lengths of the other
elements remain constant. The equation for a variable length cable can be
obtained by combining the dynamic equations of all the elements.
              <disp-formula id="Ch1.E17" content-type="numbered"><label>17</label><mml:math id="M39" display="block"><mml:mrow><mml:mi mathvariant="bold">M</mml:mi><mml:mover accent="true"><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover><mml:mo>+</mml:mo><mml:mi mathvariant="bold">Q</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="bold">0</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
            where <inline-formula><mml:math id="M40" display="inline"><mml:mi mathvariant="bold-italic">q</mml:mi></mml:math></inline-formula> represents the generalized coordinates of the system, consisting
of the generalized coordinates of each node. <inline-formula><mml:math id="M41" display="inline"><mml:mi mathvariant="bold">M</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M42" display="inline"><mml:mi mathvariant="bold">Q</mml:mi></mml:math></inline-formula> represent the
relative generalized mass matrix and generalized forces of the cable.</p>
</sec>
<sec id="Ch1.S2.SS2.SSS2">
  <label>2.2.2</label><title>Dynamic model of space rope netted system based on ANCF</title>
      <p id="d1e2426">In Gerstmayr and Shabana (2006), the ANCF
equations for the dynamics of the space rope netted system can be expressed
as
              <disp-formula id="Ch1.E18" content-type="numbered"><label>18</label><mml:math id="M43" display="block"><mml:mrow><mml:mfenced close="" open="{"><mml:mtable class="array" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:mi mathvariant="bold">M</mml:mi><mml:mover accent="true"><mml:mi mathvariant="bold-italic">q</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover><mml:mo>+</mml:mo><mml:msubsup><mml:mi mathvariant="bold-italic">C</mml:mi><mml:mi>q</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msubsup><mml:mi mathvariant="italic">λ</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">Q</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">Q</mml:mi><mml:mi mathvariant="italic">ε</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mi mathvariant="bold-italic">C</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="bold">0</mml:mn></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
            where <inline-formula><mml:math id="M44" display="inline"><mml:mi mathvariant="bold">M</mml:mi></mml:math></inline-formula> is the mass matrix, <inline-formula><mml:math id="M45" display="inline"><mml:mi mathvariant="bold-italic">q</mml:mi></mml:math></inline-formula> is the generalized coordinates,
<inline-formula><mml:math id="M46" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">Q</mml:mi><mml:mi mathvariant="italic">ε</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> generalized force vector, <inline-formula><mml:math id="M47" display="inline"><mml:mi mathvariant="bold-italic">C</mml:mi></mml:math></inline-formula> is constraint
equation, <inline-formula><mml:math id="M48" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">Q</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msup><mml:mfenced close=")" open="("><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:mo>∂</mml:mo><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:mi mathvariant="bold-italic">q</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mfenced><mml:mi mathvariant="normal">T</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M49" display="inline"><mml:mi>U</mml:mi></mml:math></inline-formula> is
strain energy. The details in Eq. (17) can also be obtained by taking
<inline-formula><mml:math id="M50" display="inline"><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula> from the variable length flexible cable model.</p>
</sec>
<sec id="Ch1.S2.SS2.SSS3">
  <label>2.2.3</label><title>Contact collision force model</title>
      <p id="d1e2575">According to the Hertz contact model, the contact force at the collision
point consists of the normal collision force <inline-formula><mml:math id="M51" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and the tangential
friction force <inline-formula><mml:math id="M52" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi>t</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>.
              <disp-formula id="Ch1.E19" content-type="numbered"><label>19</label><mml:math id="M53" display="block"><mml:mrow><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi>t</mml:mi></mml:msub><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
            Using the nonlinear spring damping model, the normal collision force is
calculated as
              <disp-formula id="Ch1.E20" content-type="numbered"><label>20</label><mml:math id="M54" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>k</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:msub><mml:mi mathvariant="italic">δ</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mi>g</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:msub><mml:mi>v</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
            where <inline-formula><mml:math id="M55" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula> is the normal unit vector and v is the normal relative velocity.
The expressions of <inline-formula><mml:math id="M56" display="inline"><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M57" display="inline"><mml:mrow><mml:msub><mml:mi>g</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> are as follows:
              <disp-formula id="Ch1.E21" content-type="numbered"><label>21</label><mml:math id="M58" display="block"><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">4</mml:mn><mml:mn mathvariant="normal">3</mml:mn></mml:mfrac></mml:mstyle><mml:mi>E</mml:mi><mml:msqrt><mml:mrow><mml:msub><mml:mi>r</mml:mi><mml:mi mathvariant="normal">eff</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:msub><mml:mi mathvariant="italic">δ</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:msqrt><mml:msub><mml:mi>g</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2</mml:mn><mml:msqrt><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">5</mml:mn><mml:mn mathvariant="normal">6</mml:mn></mml:mfrac></mml:mstyle></mml:msqrt><mml:mo>⋅</mml:mo><mml:mi mathvariant="italic">β</mml:mi><mml:mo>⋅</mml:mo><mml:msqrt><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:msub><mml:mi>m</mml:mi><mml:mi mathvariant="normal">eff</mml:mi></mml:msub></mml:mrow></mml:msqrt><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
            where <inline-formula><mml:math id="M59" display="inline"><mml:mrow><mml:msub><mml:mi>r</mml:mi><mml:mi mathvariant="normal">eff</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the radius of curvature of the contact point,
<inline-formula><mml:math id="M60" display="inline"><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M61" display="inline"><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> are the mass of the contact pair, <inline-formula><mml:math id="M62" display="inline"><mml:mi>E</mml:mi></mml:math></inline-formula> is the elastic modulus
of the material, and other unknown quantities are as follows:
              <disp-formula id="Ch1.E22" content-type="numbered"><label>22</label><mml:math id="M63" display="block"><mml:mtable class="split" rowspacing="0.2ex" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:mi mathvariant="italic">β</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mi>ln⁡</mml:mi><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>r</mml:mi><mml:mi mathvariant="normal">eff</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mo>/</mml:mo><mml:msqrt><mml:mrow><mml:msup><mml:mfenced open="[" close="]"><mml:mrow><mml:mi>ln⁡</mml:mi><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>r</mml:mi><mml:mi mathvariant="normal">eff</mml:mi></mml:msub></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msup><mml:mi mathvariant="italic">π</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:msqrt><mml:msub><mml:mi>S</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mn mathvariant="normal">2</mml:mn><mml:mi>E</mml:mi><mml:msqrt><mml:mrow><mml:msub><mml:mi>r</mml:mi><mml:mi mathvariant="normal">eff</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:msub><mml:mi mathvariant="italic">δ</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:msqrt><mml:msub><mml:mi>m</mml:mi><mml:mi mathvariant="normal">eff</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:msub><mml:mi>m</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>m</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>m</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
            Coulomb friction is used in the tangential direction as follows:
              <disp-formula id="Ch1.E23" content-type="numbered"><label>23</label><mml:math id="M64" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi>t</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi mathvariant="italic">μ</mml:mi><mml:mo>(</mml:mo><mml:mfenced open="|" close="|"><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mi mathvariant="normal">t</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mo>)</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">F</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
            where <inline-formula><mml:math id="M65" display="inline"><mml:mi mathvariant="italic">μ</mml:mi></mml:math></inline-formula> is the coefficient of friction, which is determined by the
tangential relative velocity <inline-formula><mml:math id="M66" display="inline"><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mi mathvariant="normal">t</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>.</p>
      <p id="d1e2978">In the process of catching and towing target, the cable element collides
with the target. Based on the principle of virtual work, the contact force
can be transformed to a generalized nodal force.</p>
</sec>
</sec>
<sec id="Ch1.S2.SS3">
  <label>2.3</label><title>Attitude dynamics of the service spacecraft</title>
      <p id="d1e2990">The attitude Euler angle of the service spacecraft can be selected as
<inline-formula><mml:math id="M67" display="inline"><mml:mrow><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo>=</mml:mo><mml:mo>[</mml:mo><mml:mtable class="array" columnalign="left left left"><mml:mtr><mml:mtd><mml:mi mathvariant="italic">γ</mml:mi></mml:mtd><mml:mtd><mml:mi mathvariant="italic">ϕ</mml:mi></mml:mtd><mml:mtd><mml:mi mathvariant="italic">φ</mml:mi></mml:mtd></mml:mtr></mml:mtable><mml:msup><mml:mo>]</mml:mo><mml:mi mathvariant="normal">T</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula>, then the spacecraft attitude dynamic equation can be
written as

                <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M68" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E24"><mml:mtd><mml:mtext>24</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mfenced open="{" close=""><mml:mtable class="array" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mi mathvariant="normal">Ω</mml:mi><mml:msub><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi>q</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="bold">M</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mover accent="true"><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mi mathvariant="bold">R</mml:mi><mml:mo>(</mml:mo><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>)</mml:mo><mml:mi mathvariant="bold-italic">ω</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E25"><mml:mtd><mml:mtext>25</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mi mathvariant="bold">R</mml:mi><mml:mo>(</mml:mo><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mfenced close="]" open="["><mml:mtable class="array" columnalign="center center center"><mml:mtr><mml:mtd><mml:mn mathvariant="normal">1</mml:mn></mml:mtd><mml:mtd><mml:mrow><mml:mi>tan⁡</mml:mi><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mi>sin⁡</mml:mi><mml:mi mathvariant="italic">γ</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>tan⁡</mml:mi><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mi>cos⁡</mml:mi><mml:mi mathvariant="italic">γ</mml:mi></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mrow><mml:mi>cos⁡</mml:mi><mml:mi mathvariant="italic">γ</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mi>sin⁡</mml:mi><mml:mi mathvariant="italic">γ</mml:mi></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:mi>sin⁡</mml:mi><mml:mi mathvariant="italic">γ</mml:mi></mml:mrow><mml:mrow><mml:mi>cos⁡</mml:mi><mml:mi mathvariant="italic">ϕ</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd><mml:mtd><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:mi>cos⁡</mml:mi><mml:mi mathvariant="italic">γ</mml:mi></mml:mrow><mml:mrow><mml:mi>cos⁡</mml:mi><mml:mi mathvariant="italic">ϕ</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd></mml:mtr></mml:mtable></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

            where <inline-formula><mml:math id="M69" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold">M</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the control torque of the service spacecraft. <inline-formula><mml:math id="M70" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub><mml:mo>∈</mml:mo><mml:msup><mml:mi mathvariant="bold">R</mml:mi><mml:mrow><mml:mn mathvariant="normal">3</mml:mn><mml:mo>×</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> is the rotational inertia of the service spacecraft and
<inline-formula><mml:math id="M71" display="inline"><mml:mi mathvariant="bold-italic">ω</mml:mi></mml:math></inline-formula> is the angular velocity of the service spacecraft. <inline-formula><mml:math id="M72" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi>q</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the
system uncertainty terms such as the interference generated by the mutual
collision and the large deformation of flexible netted pocket. <inline-formula><mml:math id="M73" display="inline"><mml:mi mathvariant="bold">Ω</mml:mi></mml:math></inline-formula> is
the skew–symmetric matrix of the form as follows:
            <disp-formula id="Ch1.E26" content-type="numbered"><label>26</label><mml:math id="M74" display="block"><mml:mrow><mml:mi mathvariant="bold">Ω</mml:mi><mml:mo>=</mml:mo><mml:mfenced open="|" close="|"><mml:mtable class="array" columnalign="center center center"><mml:mtr><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mi>z</mml:mi></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mi>y</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mi>z</mml:mi></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mi>x</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mi>y</mml:mi></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mi>x</mml:mi></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd></mml:mtr></mml:mtable></mml:mfenced><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p>
</sec>
</sec>
<sec id="Ch1.S3">
  <label>3</label><title>Design of spacecraft attitude double closed-loop sliding mode controller</title>
      <p id="d1e3333">In this paper, the spacecraft attitude controller uses an extended observer
to estimate <inline-formula><mml:math id="M75" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi>q</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, and a two-loop sliding mode control method to design the
attitude control torque <inline-formula><mml:math id="M76" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold">M</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, which can ensure the service spacecraft
attitude <inline-formula><mml:math id="M77" display="inline"><mml:mi mathvariant="italic">θ</mml:mi></mml:math></inline-formula> tracking the desired attitude angle <inline-formula><mml:math id="M78" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">d</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. The
outer-loop sliding mode control law calculates the attitude angular
velocity <inline-formula><mml:math id="M79" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mi mathvariant="normal">d</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and transmits <inline-formula><mml:math id="M80" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mi mathvariant="normal">d</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> to the inner loop; the
inner-loop sliding mode control law tracks the attitude angle <inline-formula><mml:math id="M81" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mi mathvariant="normal">d</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>.
The entire two-loop control system is shown in
Fig. 5.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F5"><?xmltex \currentcnt{5}?><?xmltex \def\figurename{Figure}?><label>Figure 5</label><caption><p id="d1e3412">Spacecraft attitude double-loop sliding mode controller.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022-f05.png"/>

      </fig>

<sec id="Ch1.S3.SS1">
  <label>3.1</label><title>Design of spacecraft attitude extended observer</title>
      <p id="d1e3428">Based on Eq. (23), the system state equation can be written as
            <disp-formula id="Ch1.E27" content-type="numbered"><label>27</label><mml:math id="M82" display="block"><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:mover accent="true"><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:msubsup><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msubsup><mml:mi mathvariant="bold">Ω</mml:mi><mml:msub><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="bold">R</mml:mi><mml:mfenced close=")" open="("><mml:mi mathvariant="italic">θ</mml:mi></mml:mfenced><mml:msubsup><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msubsup><mml:msub><mml:mi mathvariant="bold">M</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:mover accent="true"><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mo stretchy="false" mathvariant="normal">̃</mml:mo></mml:mover></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mi mathvariant="bold">R</mml:mi><mml:mfenced open="(" close=")"><mml:mi mathvariant="italic">θ</mml:mi></mml:mfenced><mml:msubsup><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msubsup><mml:msub><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi>q</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:mover accent="true"><mml:mi mathvariant="bold">R</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mfenced close=")" open="("><mml:mi mathvariant="italic">θ</mml:mi></mml:mfenced><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
          Referring to Freidovich and Khalil (2008), the
spacecraft attitude extended observer is designed as
            <disp-formula id="Ch1.E28" content-type="numbered"><label>28</label><mml:math id="M83" display="block"><mml:mrow><mml:mfenced close="" open="{"><mml:mtable class="array" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">y</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">y</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:mstyle displaystyle="false"><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow><mml:mi mathvariant="italic">δ</mml:mi></mml:mfrac></mml:mstyle></mml:mstyle><mml:mfenced open="(" close=")"><mml:mrow><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">y</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:mfenced></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">y</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">y</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mn mathvariant="normal">3</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:mi mathvariant="bold">R</mml:mi><mml:mfenced close=")" open="("><mml:mi mathvariant="bold-italic">θ</mml:mi></mml:mfenced><mml:msubsup><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msubsup><mml:mi mathvariant="bold-italic">u</mml:mi></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mspace width="1em" linebreak="nobreak"/><mml:mo>-</mml:mo><mml:mi mathvariant="bold">R</mml:mi><mml:mfenced open="(" close=")"><mml:mi mathvariant="bold-italic">θ</mml:mi></mml:mfenced><mml:msubsup><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msubsup><mml:mi mathvariant="bold">Ω</mml:mi><mml:msub><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mo>+</mml:mo><mml:mstyle displaystyle="false"><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msup><mml:mi mathvariant="italic">δ</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mstyle></mml:mstyle><mml:mfenced close=")" open="("><mml:mrow><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">y</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:mfenced></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">y</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mn mathvariant="normal">3</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="false"><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msup><mml:mi mathvariant="italic">δ</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mstyle></mml:mstyle><mml:mfenced open="(" close=")"><mml:mrow><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">y</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:math></disp-formula>
          where <inline-formula><mml:math id="M84" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">y</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is the state estimation of <inline-formula><mml:math id="M85" display="inline"><mml:mi mathvariant="bold-italic">θ</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M86" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">y</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is the
state estimation of <inline-formula><mml:math id="M87" display="inline"><mml:mover accent="true"><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover></mml:math></inline-formula> and <inline-formula><mml:math id="M88" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">y</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mn mathvariant="normal">3</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is the state estimation
of <inline-formula><mml:math id="M89" display="inline"><mml:mover accent="true"><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mo stretchy="false" mathvariant="normal">̃</mml:mo></mml:mover></mml:math></inline-formula>. <inline-formula><mml:math id="M90" display="inline"><mml:mi mathvariant="italic">δ</mml:mi></mml:math></inline-formula> is taken as follows:
            <disp-formula id="Ch1.E29" content-type="numbered"><label>29</label><mml:math id="M91" display="block"><mml:mrow><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mi mathvariant="italic">δ</mml:mi></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:mfenced open="{" close=""><mml:mtable class="array" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:mn mathvariant="normal">100</mml:mn><mml:msup><mml:mi>t</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msup></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn mathvariant="normal">100</mml:mn></mml:mtd></mml:mtr></mml:mtable></mml:mfenced><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p>
</sec>
<sec id="Ch1.S3.SS2">
  <label>3.2</label><title>Design of outer loop sliding mode control</title>
      <p id="d1e3858">The integral sliding mode (Liu
and Wang, 2011) is adopted to realize the sliding mode surface design.
Taking <inline-formula><mml:math id="M92" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">e</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mi>d</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">y</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, the sliding mode surface of the
outer loop can be written as follows:
            <disp-formula id="Ch1.E30" content-type="numbered"><label>30</label><mml:math id="M93" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>w</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">e</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="bold">K</mml:mi><mml:mi>w</mml:mi></mml:msub><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mn mathvariant="normal">0</mml:mn><mml:mi>t</mml:mi></mml:munderover><mml:msub><mml:mi mathvariant="bold-italic">e</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:msub><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
          where <inline-formula><mml:math id="M94" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold">K</mml:mi><mml:mi>w</mml:mi></mml:msub><mml:mo>&gt;</mml:mo><mml:mn mathvariant="bold">0</mml:mn></mml:mrow></mml:math></inline-formula> is <inline-formula><mml:math id="M95" display="inline"><mml:mrow><mml:mn mathvariant="normal">3</mml:mn><mml:mo>×</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:math></inline-formula> diagonal matrix.</p>
      <p id="d1e3959">The derivation of Eq. (29) can be obtained as follows:
            <disp-formula id="Ch1.E31" content-type="numbered"><label>31</label><mml:math id="M96" display="block"><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi>w</mml:mi></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">e</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi mathvariant="italic">θ</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mi>w</mml:mi></mml:msub><mml:msub><mml:mi mathvariant="bold-italic">e</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi>d</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mover accent="true"><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mo>+</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mi>w</mml:mi></mml:msub><mml:msub><mml:mi mathvariant="bold-italic">e</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi>d</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi mathvariant="bold">R</mml:mi><mml:mo>(</mml:mo><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo>)</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mi>d</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:mi mathvariant="bold">R</mml:mi><mml:mo>(</mml:mo><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo>)</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mi>e</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mi>w</mml:mi></mml:msub><mml:msub><mml:mi mathvariant="bold-italic">e</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
          where <inline-formula><mml:math id="M97" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mi>e</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mi>d</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> and outer loop control law is
            <disp-formula id="Ch1.E32" content-type="numbered"><label>32</label><mml:math id="M98" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mi>d</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msup><mml:mi mathvariant="bold">R</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup><mml:mo>(</mml:mo><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo>)</mml:mo><mml:mo>(</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi>d</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mi>w</mml:mi></mml:msub><mml:msub><mml:mi mathvariant="bold-italic">e</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>+</mml:mo><mml:msup><mml:mi mathvariant="bold">R</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup><mml:mo>(</mml:mo><mml:mi mathvariant="bold-italic">θ</mml:mi><mml:mo>)</mml:mo><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mi>w</mml:mi></mml:msub><mml:msub><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>w</mml:mi></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
          where <inline-formula><mml:math id="M99" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mi>w</mml:mi></mml:msub><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula>. Taking Lyapunov function <inline-formula><mml:math id="M100" display="inline"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>w</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mn mathvariant="normal">1</mml:mn><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:msubsup><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>w</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msubsup><mml:msub><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>w</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>
(Pham et al., 2019), it can be proved that
            <disp-formula id="Ch1.E33" content-type="numbered"><label>33</label><mml:math id="M101" display="block"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>V</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi>w</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mi>w</mml:mi></mml:msub><mml:msup><mml:mfenced close="∥" open="∥"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>w</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msubsup><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>w</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msubsup><mml:mi mathvariant="bold">R</mml:mi><mml:mo>(</mml:mo><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>)</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mi>e</mml:mi></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
          when the inner loop converges quickly and <inline-formula><mml:math id="M102" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mi>e</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is sufficiently
small. There is <inline-formula><mml:math id="M103" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>V</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi>w</mml:mi></mml:msub><mml:mo>≤</mml:mo></mml:mrow></mml:math></inline-formula> and the error <inline-formula><mml:math id="M104" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">e</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is
asymptotically stable.</p>
</sec>
<sec id="Ch1.S3.SS3">
  <label>3.3</label><title>Design of inner loop sliding mode control</title>
      <p id="d1e4348">The sliding surface of the inner loop is designed as
            <disp-formula id="Ch1.E34" content-type="numbered"><label>34</label><mml:math id="M105" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mi>e</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="bold">K</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mn mathvariant="normal">0</mml:mn><mml:mi>t</mml:mi></mml:munderover><mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mi>e</mml:mi></mml:msub><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
          where <inline-formula><mml:math id="M106" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="bold">K</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula> is <inline-formula><mml:math id="M107" display="inline"><mml:mrow><mml:mn mathvariant="normal">3</mml:mn><mml:mo>×</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:math></inline-formula> diagonal matrix. The derivation of Eq. (33)
can be obtained as follows:
            <disp-formula id="Ch1.E35" content-type="numbered"><label>35</label><mml:math id="M108" display="block"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi>n</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi>e</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="bold">K</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mi>e</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi>d</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msubsup><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msubsup><mml:mo>(</mml:mo><mml:mi mathvariant="bold">Ω</mml:mi><mml:msub><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi>q</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="bold">M</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="bold">K</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mi>e</mml:mi></mml:msub><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
          The control law of the inner loop is designed as
            <disp-formula id="Ch1.E36" content-type="numbered"><label>36</label><mml:math id="M109" display="block"><mml:mrow><mml:mtable class="split" rowspacing="0.2ex" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi mathvariant="bold">M</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mi mathvariant="bold">R</mml:mi><mml:mfenced close=")" open="("><mml:mi mathvariant="italic">θ</mml:mi></mml:mfenced><mml:msubsup><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msubsup><mml:mi mathvariant="bold">Ω</mml:mi><mml:msub><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub><mml:msup><mml:mi mathvariant="bold">R</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup><mml:mfenced open="(" close=")"><mml:mi mathvariant="italic">θ</mml:mi></mml:mfenced><mml:mo>(</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi>d</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="bold">K</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:msub><mml:mi mathvariant="bold-italic">ω</mml:mi><mml:mi>e</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mi mathvariant="normal">sat</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>+</mml:mo><mml:mi mathvariant="bold-italic">λ</mml:mi><mml:msub><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
          where <inline-formula><mml:math id="M110" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>&gt;</mml:mo><mml:mo>max⁡</mml:mo><mml:mfenced close="|" open="|"><mml:mrow><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mfenced open="(" close=")"><mml:mi>i</mml:mi></mml:mfenced></mml:mrow></mml:mfenced></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M111" display="inline"><mml:mrow><mml:mi mathvariant="bold-italic">λ</mml:mi><mml:mo>=</mml:mo><mml:mfenced open="[" close="]"><mml:mtable class="array" columnalign="center center"><mml:mtr><mml:mtd><mml:mi mathvariant="italic">λ</mml:mi></mml:mtd><mml:mtd/></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mi mathvariant="italic">λ</mml:mi></mml:mtd></mml:mtr></mml:mtable></mml:mfenced><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula>. To further eliminate the “chattering”, the
saturation function <inline-formula><mml:math id="M112" display="inline"><mml:mrow><mml:mi mathvariant="normal">sat</mml:mi><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> is used as the switching function.</p>
      <p id="d1e4712">Taking Lyapunov function, <inline-formula><mml:math id="M113" display="inline"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mn mathvariant="normal">1</mml:mn><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:msubsup><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>n</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msubsup><mml:msub><mml:mi mathvariant="bold">I</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub><mml:msub><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, it can be
proved (Slotine and Sastry, 1983) that
            <disp-formula id="Ch1.E37" content-type="numbered"><label>37</label><mml:math id="M114" display="block"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>V</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi>n</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo movablelimits="false">∑</mml:mo><mml:mfenced open="|" close="|"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mo>-</mml:mo><mml:msubsup><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>n</mml:mi><mml:mi mathvariant="normal">T</mml:mi></mml:msubsup><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mo>-</mml:mo><mml:mi mathvariant="italic">λ</mml:mi><mml:msup><mml:mfenced close="∥" open="∥"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>≤</mml:mo><mml:mo>-</mml:mo><mml:mi mathvariant="italic">λ</mml:mi><mml:msup><mml:mfenced open="∥" close="∥"><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>≤</mml:mo><mml:mo>-</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
          The Eq. (36) proves that <inline-formula><mml:math id="M115" display="inline"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is exponential convergence.</p>
</sec>
</sec>
<sec id="Ch1.S4">
  <label>4</label><title>Simulation analysis of capture and towing process</title>
      <p id="d1e4856">This paper utilizes the multi-body dynamic simulation software (MBDyn)
developed by the author's laboratory to perform dynamic simulation. The
dynamic parameters of the service spacecraft and captured target are in
Table 1.</p>

<?xmltex \floatpos{t}?><table-wrap id="Ch1.T1"><?xmltex \currentcnt{1}?><label>Table 1</label><caption><p id="d1e4862">The dynamic parameters of the service spacecraft and captured target.</p></caption><oasis:table frame="topbot"><oasis:tgroup cols="3">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="right"/>
     <oasis:colspec colnum="3" colname="col3" align="right"/>
     <oasis:thead>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">Parameter name</oasis:entry>
         <oasis:entry colname="col2">Service spacecraft</oasis:entry>
         <oasis:entry colname="col3">Captured target</oasis:entry>
       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>
         <oasis:entry colname="col1">Mass (kg)</oasis:entry>
         <oasis:entry colname="col2">1210</oasis:entry>
         <oasis:entry colname="col3">200</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M116" display="inline"><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mi>x</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> (kg m<inline-formula><mml:math id="M117" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">4282</oasis:entry>
         <oasis:entry colname="col3">1210</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M118" display="inline"><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mi>y</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> (kg m<inline-formula><mml:math id="M119" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">12 736</oasis:entry>
         <oasis:entry colname="col3">1210</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M120" display="inline"><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mi>z</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> (kg m<inline-formula><mml:math id="M121" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">14 498</oasis:entry>
         <oasis:entry colname="col3">1210</oasis:entry>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup></oasis:table></table-wrap>

      <p id="d1e4999">The material parameters of the inflatable deployment rod and rope net are in
Table 2.</p>

<?xmltex \floatpos{t}?><table-wrap id="Ch1.T2"><?xmltex \currentcnt{2}?><label>Table 2</label><caption><p id="d1e5006">The material parameters of the inflatable deployment rod and rope net.</p></caption><oasis:table frame="topbot"><oasis:tgroup cols="3">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="right"/>
     <oasis:colspec colnum="3" colname="col3" align="right"/>
     <oasis:thead>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">Parameter name</oasis:entry>
         <oasis:entry colname="col2">Nets</oasis:entry>
         <oasis:entry colname="col3">Inflatable rods</oasis:entry>
       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>
         <oasis:entry colname="col1">Diameter (m)</oasis:entry>
         <oasis:entry colname="col2">0.004</oasis:entry>
         <oasis:entry colname="col3">0.01</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Density (kg m<inline-formula><mml:math id="M122" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">1430</oasis:entry>
         <oasis:entry colname="col3">164</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Poisson ratio</oasis:entry>
         <oasis:entry colname="col2">0.3</oasis:entry>
         <oasis:entry colname="col3">0.3</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Modulus of elasticity (GPa)</oasis:entry>
         <oasis:entry colname="col2">12</oasis:entry>
         <oasis:entry colname="col3">0.75</oasis:entry>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup></oasis:table></table-wrap>

      <?xmltex \floatpos{t}?><fig id="Ch1.F6" specific-use="star"><?xmltex \currentcnt{6}?><?xmltex \def\figurename{Figure}?><label>Figure 6</label><caption><p id="d1e5100">Capturing and attitude maneuver process.</p></caption>
        <?xmltex \igopts{width=398.338583pt}?><graphic xlink:href="https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022-f06.png"/>

      </fig>

<sec id="Ch1.S4.SS1">
  <label>4.1</label><title>Service spacecraft attitude maneuvering process simulation</title>
      <p id="d1e5116">The process of the service spacecraft capture target and attitude maneuver
is designed as follows: 0–20 s is the process of net retraction during which
the service spacecraft remains stationary; 20–60 s is the process of service
spacecraft attitude maneuvering 90<inline-formula><mml:math id="M123" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula> around <inline-formula><mml:math id="M124" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> axis; and 60–100 s is
the process of service spacecraft attitude maintenance. In order to ensure
the smooth and continuous motion process, the rope retrieval length in the
0–20 s and the spacecraft attitude maneuvering in 20–60 s are planned using
the 5 times interpolation, and then the angular velocity and angular
acceleration are 0 at 20 and 60 s can be determined.</p>
      <p id="d1e5135">Set the parameters of the sliding mode controller <inline-formula><mml:math id="M125" display="inline"><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>w</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M126" display="inline"><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">10</mml:mn></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M127" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.01</mml:mn></mml:mrow></mml:math></inline-formula>, the parameters in the saturation function <inline-formula><mml:math id="M128" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.5</mml:mn></mml:mrow></mml:math></inline-formula> , and
the parameters of the expansion observer <inline-formula><mml:math id="M129" display="inline"><mml:mrow><mml:mi>R</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">100</mml:mn></mml:mrow></mml:math></inline-formula>. The dynamics and control
closed-loop simulation of the attitude maneuvering process is performed, and
the calculation results are in Fig. 6.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F7" specific-use="star"><?xmltex \currentcnt{7}?><?xmltex \def\figurename{Figure}?><label>Figure 7</label><caption><p id="d1e5209">Service spacecraft attitude motion state.</p></caption>
          <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022-f07.png"/>

        </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F8" specific-use="star"><?xmltex \currentcnt{8}?><?xmltex \def\figurename{Figure}?><label>Figure 8</label><caption><p id="d1e5221">Attitude extended observer results.</p></caption>
          <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022-f08.png"/>

        </fig>

      <p id="d1e5230">The estimates of the attitude motion state and disturbance by the expansion
observer are shown in Fig. 8. Comparing
Figs. 8 and 7, it
can be seen that the estimates of the observer are consistent with the
motion process, which illustrates the validity of the observations.</p>

      <?xmltex \floatpos{p}?><fig id="Ch1.F9" specific-use="star"><?xmltex \currentcnt{9}?><?xmltex \def\figurename{Figure}?><label>Figure 9</label><caption><p id="d1e5235">Error <inline-formula><mml:math id="M130" display="inline"><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> in each direction.</p></caption>
          <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022-f09.png"/>

        </fig>

      <p id="d1e5255">As shown in Fig. 9, compared with ordinary
sliding mode control, the error of extended observer controller is reduced
by 1 order of magnitude. Therefore, the control accuracy of the
double-loop sliding mode attitude controller based on the extended observer
is significantly higher than the ordinary sliding mode control. It indicates
that the extended observer compensates the disturbance effectively, and the
service spacecraft attitude controller designed in this paper can meet the
stability requirements of the attitude maneuvering process.</p>

      <?xmltex \floatpos{p}?><fig id="Ch1.F10" specific-use="star"><?xmltex \currentcnt{10}?><?xmltex \def\figurename{Figure}?><label>Figure 10</label><caption><p id="d1e5260">Service spacecraft control torque.</p></caption>
          <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022-f10.png"/>

        </fig>

      <p id="d1e5270">As shown in Fig. 10, the net is a centrosymmetric
structure and the initial state of target is located in the center, so the
control torque <inline-formula><mml:math id="M131" display="inline"><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mi mathvariant="normal">c</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> in each direction is small during 0 to 20 s. Due to
the collision increases at 20 and 60 s, the control torque <inline-formula><mml:math id="M132" display="inline"><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">c</mml:mi><mml:mi>z</mml:mi></mml:msub></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">c</mml:mi><mml:mi>y</mml:mi></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> produces buffeting. And during 20 to 60 s, the service
spacecraft attitude maneuvers around the <inline-formula><mml:math id="M133" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> axis, so the control torque
<inline-formula><mml:math id="M134" display="inline"><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">c</mml:mi><mml:mi>x</mml:mi></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is much higher than <inline-formula><mml:math id="M135" display="inline"><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">c</mml:mi><mml:mi>y</mml:mi></mml:msub></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">c</mml:mi><mml:mi>z</mml:mi></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>. The
flexible structure of the system generates certain perturbations. So, in the
attitude holding process after 60 s, the control torque curve still fluctuates
until the system stabilizes.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F11" specific-use="star"><?xmltex \currentcnt{11}?><?xmltex \def\figurename{Figure}?><label>Figure 11</label><caption><p id="d1e5360">Service spacecraft attitude dynamic response.</p></caption>
          <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://ms.copernicus.org/articles/13/751/2022/ms-13-751-2022-f11.png"/>

        </fig>

      <p id="d1e5369">It can be seen from the collision force curve in
Fig. 11a that the collision force between the
capture target and the rope net is mainly concentrated in the <inline-formula><mml:math id="M136" display="inline"><mml:mi>y</mml:mi></mml:math></inline-formula> direction.
During 0 to 15 s, tensile force of closing rope is mainly caused by the
bending of the inflatable rod, and the value is small as shown in
Fig. 11b. During 15 to 20 s, the collision
force is large, so the tensile force of closing rope is also with large
step. As shown in Fig. 11b, the tension force
of the closing rope tends to be stable and non-zero after 20 s, indicating
that the rope net is tightly wrapped in the target and there is no relative
movement between them. It can also be seen from the collision force curve
that the close wrapping of the rope net around the target significantly
reduces the collision in the maneuvering process. Thus, the interference is
further reduced, and the control torque of spacecraft is facilitated as shown
in Fig. 10.</p>
</sec>
</sec>
<sec id="Ch1.S5" sec-type="conclusions">
  <label>5</label><title>Conclusion</title>
      <p id="d1e5389">In this paper, the dynamics and stability control of the space net pocket
capture and towing process are studied, and the following conclusions are
obtained.
<?xmltex \hack{\newpage}?>
<list list-type="order"><list-item>
      <p id="d1e5396">The dynamic model of the space netted pocket capture system is carried
out. Based on the ANCF method, the dynamic model of the space rope system
can respond to the large deformation properties of the catching mechanism.
Based on the dynamic model of the closing rope established from the ANCF
flexible cable theory, the simulation of the non-cooperative target capture
process is realized, and the dynamic response of the closing rope recovery
can be analyzed.</p></list-item><list-item>
      <p id="d1e5400">A sliding mode controller based on the extended observer is designed with
reference to the service spacecraft dynamic equations. Only the attitude of
the spacecraft is needed to complete the maneuver control, as the extended
observer can observe the velocity angular velocity of the spacecraft.
Moreover, based on spacecraft attitude extended observer, the sliding mode
controller can ensure the high accuracy and stable attitude control of the
service spacecraft.</p></list-item><list-item>
      <p id="d1e5404">Closed-loop dynamics and control simulations are performed for the
capture and towing process of the service spacecraft. The simulation results
can reflect the large deformation of the net pocket and the contact
collision between the net pocket and the target during the capture and
towing process. The closed-loop simulation results verify that the control
accuracy is improved by 1 order of magnitude under the interference
compensation of the expansion observer, which is beneficial for realizing
the stable attitude control of service spacecrafts with large deformation
structure.</p></list-item></list></p>
</sec>

      
      </body>
    <back><notes notes-type="codeavailability"><title>Code availability</title>

      <p id="d1e5411">The simulation software is jointly developed with a third party. The other party is for commercial purposes and does not want the software code to be disclosed.</p>
  </notes><notes notes-type="dataavailability"><title>Data availability</title>

      <p id="d1e5417">In this paper, the simulation conditions and relevant parameter data are listed in Tables 1 and 2. And the effective data results are shown by curves in Figs. 7–11.</p>
  </notes><notes notes-type="authorcontribution"><title>Author contributions</title>

      <p id="d1e5423">CT was mainly responsible for the calculation and data analysis of the simulation examples in this paper.
ZH was mainly responsible for the modeling of simulation examples in this paper.
CW was mainly responsible for the development of the dynamics software in this paper.
YZ was mainly responsible for the research of dynamics and control theory in this paper.</p>
  </notes><notes notes-type="competinginterests"><title>Competing interests</title>

      <p id="d1e5429">The contact author has declared that none of the authors has any competing interests.</p>
  </notes><notes notes-type="disclaimer"><title>Disclaimer</title>

      <p id="d1e5435">Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.</p>
  </notes><notes notes-type="financialsupport"><title>Financial support</title>

      <p id="d1e5441">This research has been supported by the Heilongjiang Postdoctoral Fund (grant no.  LBH-Z21141) and the National Natural Science Foundation of China (grant no. 12102316).</p>
  </notes><notes notes-type="reviewstatement"><title>Review statement</title>

      <p id="d1e5447">This paper was edited by Daniel Condurache and reviewed by two anonymous referees.</p>
  </notes><ref-list>
    <title>References</title>

      <ref id="bib1.bib1"><label>1</label><?label 1?><mixed-citation>Barmin, I. V., Dunham, D. W., Kulagin, V. P., Savinykh, V. P., and Tsvetkov,
V. Y.: Rings of debris in near-Earth space, Solar Syst. Res., 48,
593–600, <ext-link xlink:href="https://doi.org/10.1134/s0038094614070041" ext-link-type="DOI">10.1134/s0038094614070041</ext-link>, 2014.</mixed-citation></ref>
      <ref id="bib1.bib2"><label>2</label><?label 1?><mixed-citation>
Biesbroek, R., Innocenti, L., Wolahan, A., and Serrano, S. M.: E.DEORBIT –
ESA'S ACTIVE DEBRIS REMOVAL MISSION, JBIS-J. Brit. Interpla., 70, 143–151, 2017.</mixed-citation></ref>
      <ref id="bib1.bib3"><label>3</label><?label 1?><mixed-citation>Bonnal, C., Ruault, J.-M., and Desjean, M.-C.: Active debris removal: Recent
progress and current trends, Acta Astronaut., 85, 51–60,
<ext-link xlink:href="https://doi.org/10.1016/j.actaastro.2012.11.009" ext-link-type="DOI">10.1016/j.actaastro.2012.11.009</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bib4"><label>4</label><?label 1?><mixed-citation>Chen, S.: The Space Debris Problem, Asian Perspec., 35, 537–558,
<ext-link xlink:href="https://doi.org/10.1017/S1752971911000145" ext-link-type="DOI">10.1017/S1752971911000145</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bib5"><label>5</label><?label 1?><mixed-citation>Forshaw, J. L., Aglietti, G. S., Navarathinam, N., Kadhem, H., Salmon, T.,
Pisseloup, A., Joffre, E., Chabot, T., Retat, I., Axthelm, R., Barraclough,
S., Ratcliffe, A., Bernal, C., Chaumette, F., Pollini, A., and Steyn, W. H.:
RemoveDEBRIS: An in-orbit active debris removal demonstration mission, Acta
Astronaut., 127, 448–463,
<ext-link xlink:href="https://doi.org/10.1016/j.actaastro.2016.06.018" ext-link-type="DOI">10.1016/j.actaastro.2016.06.018</ext-link>, 2016.</mixed-citation></ref>
      <ref id="bib1.bib6"><label>6</label><?label 1?><mixed-citation>Forshaw, J. L., Aglietti, G. S., Salmon, T., Retat, I., Roe, M., Burgess,
C., Chabot, T., Pisseloup, A., Phipps, A., Bernal, C., Chaumette, F.,
Pollini, A., and Steyn, W. H.: Final payload test results for the
RemoveDebris active debris removal mission, Acta Astronaut., 138,
326–342, <ext-link xlink:href="https://doi.org/10.1016/j.actaastro.2017.06.003" ext-link-type="DOI">10.1016/j.actaastro.2017.06.003</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib7"><label>7</label><?label 1?><mixed-citation>Freidovich, L. B. and Khalil, H. K.: Performance Recovery of
Feedback-Linearization-Based Designs, IEEE T. Autom. Control, 53,
2324–2334, <ext-link xlink:href="https://doi.org/10.1109/TAC.2008.2006821" ext-link-type="DOI">10.1109/TAC.2008.2006821</ext-link>, 2008.</mixed-citation></ref>
      <ref id="bib1.bib8"><label>8</label><?label 1?><mixed-citation>Gerstmayr, J. and Shabana, A. A.: Analysis of thin beams and cables using
the absolute nodal co-ordinate formulation, Nonlinear Dynam., 45, 109–130,
<ext-link xlink:href="https://doi.org/10.1007/s11071-006-1856-1" ext-link-type="DOI">10.1007/s11071-006-1856-1</ext-link>, 2006.</mixed-citation></ref>
      <ref id="bib1.bib9"><label>9</label><?label 1?><mixed-citation>Huang, P., Zhang, F., Ma, J., Meng, Z., and Liu, Z.: Dynamics and
configuration control of the Maneuvering-Net Space Robot System, Adv. Space
Res., 55, 1004–1014, <ext-link xlink:href="https://doi.org/10.1016/j.asr.2014.11.009" ext-link-type="DOI">10.1016/j.asr.2014.11.009</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bib10"><label>10</label><?label 1?><mixed-citation>Huang, P., Zhang, F., Meng, Z., and Liu, Z.: Adaptive control for space
debris removal with uncertain kinematics, dynamics and states, Acta
Astronaut., 128, 416–430, <ext-link xlink:href="https://doi.org/10.1016/j.actaastro.2016.07.043" ext-link-type="DOI">10.1016/j.actaastro.2016.07.043</ext-link>,
2016.</mixed-citation></ref>
      <ref id="bib1.bib11"><label>11</label><?label 1?><mixed-citation>Liu, J. and Wang, X.: Advanced sliding mode control for mechanical systems, in: Advanced sliding mode control for mechanical systems, Higher Education Press, <ext-link xlink:href="https://doi.org/10.1007/978-3-642-20907-9" ext-link-type="DOI">10.1007/978-3-642-20907-9</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bib12"><label>12</label><?label 1?><mixed-citation>Meng, Z., Huang, P., and Guo, J.: Approach Modeling and Control of an
Autonomous Maneuverable Space Net, IEEE T. Aero. Elec. Sys., 53,
2651–2661, <ext-link xlink:href="https://doi.org/10.1109/TAES.2017.2709794" ext-link-type="DOI">10.1109/TAES.2017.2709794</ext-link>, 2017.
</mixed-citation></ref><?xmltex \hack{\newpage}?>
      <ref id="bib1.bib13"><label>13</label><?label 1?><mixed-citation>Ming, Z., Yuanwen, Z., Leping, Y., and Qingbin, Z.: The Deployment Dynamics of a New Mixed SpaceWeb, in: Conference Proceedings of 2017 3rd Ieee International Conference on Control Science and Systems Engineering (iccsse), New York, 661–666, <ext-link xlink:href="https://doi.org/10.1109/CCSSE.2017.8088015" ext-link-type="DOI">10.1109/CCSSE.2017.8088015</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib14"><label>14</label><?label 1?><mixed-citation>Nishida, S.-I. and Kawamoto, S.: Strategy for capturing of a tumbling space
debris, Acta Astronaut., 68, 113–120,
<ext-link xlink:href="https://doi.org/10.1016/j.actaastro.2010.06.045" ext-link-type="DOI">10.1016/j.actaastro.2010.06.045</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bib15"><label>15</label><?label 1?><mixed-citation>Pham, D. B., Kim, J., Lee, S.-G., and Gwak, K.-W.: Double-Loop Control with
Hierarchical Sliding Mode and Proportional Integral Loop for 2D Ridable
Ballbot, Int. J. Precis. Eng. Man., 20, 1519–1532,
<ext-link xlink:href="https://doi.org/10.1007/s12541-019-00139-4" ext-link-type="DOI">10.1007/s12541-019-00139-4</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bib16"><label>16</label><?label 1?><mixed-citation>Shan, M., Guo, J., and Gill, E.: Review and comparison of active space
debris capturing and removal methods, Prog. Aerosp. Sci., 80,
18–32, <ext-link xlink:href="https://doi.org/10.1016/j.paerosci.2015.11.001" ext-link-type="DOI">10.1016/j.paerosci.2015.11.001</ext-link>, 2016.</mixed-citation></ref>
      <ref id="bib1.bib17"><label>17</label><?label 1?><mixed-citation>Slotine, J. J. and Sastry, S. S.: Tracking control of non-linear systems using sliding surfaces with application to robot manipulators, with application to robot manipulators, International Journal of Control, 38, 465–492, <ext-link xlink:href="https://doi.org/10.1080/00207178308933088" ext-link-type="DOI">10.1080/00207178308933088</ext-link>, 1983.</mixed-citation></ref>
      <ref id="bib1.bib18"><label>18</label><?label 1?><mixed-citation>Wei, C., Liu, H., Tan, C., Liu, Y., and Zhao, Y.: Rapid attitude maneuver of the space tether net capture system using active disturbance rejection control, Mech. Sci., 10, 575–587, <ext-link xlink:href="https://doi.org/10.5194/ms-10-575-2019" ext-link-type="DOI">10.5194/ms-10-575-2019</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bib19"><label>19</label><?label 1?><mixed-citation>Xu, B., Yang, Y., Yan, Y., and Zhang, B.: Bionics design and dynamics
analysis of space webs based on spider predation, Acta Astronaut., 159,
294–307, <ext-link xlink:href="https://doi.org/10.1016/j.actaastro.2019.03.045" ext-link-type="DOI">10.1016/j.actaastro.2019.03.045</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bib20"><label>20</label><?label 1?><mixed-citation>
Zhang, F., Huang, P., Meng, Z., Zhang, Y., and Liu, Z.: Dynamics analysis
and controller design for maneuverable tethered space net robot, J.
Guid. Control Dynam., 40, 2828–2843, 2017.</mixed-citation></ref>

  </ref-list></back>
    <!--<article-title-html>Dynamic and sliding mode control of  space netted pocket system capturing and attitude maneuver non-cooperative target</article-title-html>
<abstract-html/>
<ref-html id="bib1.bib1"><label>1</label><mixed-citation>
Barmin, I. V., Dunham, D. W., Kulagin, V. P., Savinykh, V. P., and Tsvetkov,
V. Y.: Rings of debris in near-Earth space, Solar Syst. Res., 48,
593–600, <a href="https://doi.org/10.1134/s0038094614070041" target="_blank">https://doi.org/10.1134/s0038094614070041</a>, 2014.
</mixed-citation></ref-html>
<ref-html id="bib1.bib2"><label>2</label><mixed-citation>
Biesbroek, R., Innocenti, L., Wolahan, A., and Serrano, S. M.: E.DEORBIT –
ESA'S ACTIVE DEBRIS REMOVAL MISSION, JBIS-J. Brit. Interpla., 70, 143–151, 2017.
</mixed-citation></ref-html>
<ref-html id="bib1.bib3"><label>3</label><mixed-citation>
Bonnal, C., Ruault, J.-M., and Desjean, M.-C.: Active debris removal: Recent
progress and current trends, Acta Astronaut., 85, 51–60,
<a href="https://doi.org/10.1016/j.actaastro.2012.11.009" target="_blank">https://doi.org/10.1016/j.actaastro.2012.11.009</a>, 2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib4"><label>4</label><mixed-citation>
Chen, S.: The Space Debris Problem, Asian Perspec., 35, 537–558,
<a href="https://doi.org/10.1017/S1752971911000145" target="_blank">https://doi.org/10.1017/S1752971911000145</a>, 2011.
</mixed-citation></ref-html>
<ref-html id="bib1.bib5"><label>5</label><mixed-citation>
Forshaw, J. L., Aglietti, G. S., Navarathinam, N., Kadhem, H., Salmon, T.,
Pisseloup, A., Joffre, E., Chabot, T., Retat, I., Axthelm, R., Barraclough,
S., Ratcliffe, A., Bernal, C., Chaumette, F., Pollini, A., and Steyn, W. H.:
RemoveDEBRIS: An in-orbit active debris removal demonstration mission, Acta
Astronaut., 127, 448–463,
<a href="https://doi.org/10.1016/j.actaastro.2016.06.018" target="_blank">https://doi.org/10.1016/j.actaastro.2016.06.018</a>, 2016.
</mixed-citation></ref-html>
<ref-html id="bib1.bib6"><label>6</label><mixed-citation>
Forshaw, J. L., Aglietti, G. S., Salmon, T., Retat, I., Roe, M., Burgess,
C., Chabot, T., Pisseloup, A., Phipps, A., Bernal, C., Chaumette, F.,
Pollini, A., and Steyn, W. H.: Final payload test results for the
RemoveDebris active debris removal mission, Acta Astronaut., 138,
326–342, <a href="https://doi.org/10.1016/j.actaastro.2017.06.003" target="_blank">https://doi.org/10.1016/j.actaastro.2017.06.003</a>, 2017.
</mixed-citation></ref-html>
<ref-html id="bib1.bib7"><label>7</label><mixed-citation>
Freidovich, L. B. and Khalil, H. K.: Performance Recovery of
Feedback-Linearization-Based Designs, IEEE T. Autom. Control, 53,
2324–2334, <a href="https://doi.org/10.1109/TAC.2008.2006821" target="_blank">https://doi.org/10.1109/TAC.2008.2006821</a>, 2008.
</mixed-citation></ref-html>
<ref-html id="bib1.bib8"><label>8</label><mixed-citation>
Gerstmayr, J. and Shabana, A. A.: Analysis of thin beams and cables using
the absolute nodal co-ordinate formulation, Nonlinear Dynam., 45, 109–130,
<a href="https://doi.org/10.1007/s11071-006-1856-1" target="_blank">https://doi.org/10.1007/s11071-006-1856-1</a>, 2006.
</mixed-citation></ref-html>
<ref-html id="bib1.bib9"><label>9</label><mixed-citation>
Huang, P., Zhang, F., Ma, J., Meng, Z., and Liu, Z.: Dynamics and
configuration control of the Maneuvering-Net Space Robot System, Adv. Space
Res., 55, 1004–1014, <a href="https://doi.org/10.1016/j.asr.2014.11.009" target="_blank">https://doi.org/10.1016/j.asr.2014.11.009</a>, 2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib10"><label>10</label><mixed-citation>
Huang, P., Zhang, F., Meng, Z., and Liu, Z.: Adaptive control for space
debris removal with uncertain kinematics, dynamics and states, Acta
Astronaut., 128, 416–430, <a href="https://doi.org/10.1016/j.actaastro.2016.07.043" target="_blank">https://doi.org/10.1016/j.actaastro.2016.07.043</a>,
2016.
</mixed-citation></ref-html>
<ref-html id="bib1.bib11"><label>11</label><mixed-citation>
Liu, J. and Wang, X.: Advanced sliding mode control for mechanical systems, in: Advanced sliding mode control for mechanical systems, Higher Education Press, <a href="https://doi.org/10.1007/978-3-642-20907-9" target="_blank">https://doi.org/10.1007/978-3-642-20907-9</a>, 2011.
</mixed-citation></ref-html>
<ref-html id="bib1.bib12"><label>12</label><mixed-citation>
Meng, Z., Huang, P., and Guo, J.: Approach Modeling and Control of an
Autonomous Maneuverable Space Net, IEEE T. Aero. Elec. Sys., 53,
2651–2661, <a href="https://doi.org/10.1109/TAES.2017.2709794" target="_blank">https://doi.org/10.1109/TAES.2017.2709794</a>, 2017.

</mixed-citation></ref-html>
<ref-html id="bib1.bib13"><label>13</label><mixed-citation>
Ming, Z., Yuanwen, Z., Leping, Y., and Qingbin, Z.: The Deployment Dynamics of a New Mixed SpaceWeb, in: Conference Proceedings of 2017 3rd Ieee International Conference on Control Science and Systems Engineering (iccsse), New York, 661–666, <a href="https://doi.org/10.1109/CCSSE.2017.8088015" target="_blank">https://doi.org/10.1109/CCSSE.2017.8088015</a>, 2017.
</mixed-citation></ref-html>
<ref-html id="bib1.bib14"><label>14</label><mixed-citation>
Nishida, S.-I. and Kawamoto, S.: Strategy for capturing of a tumbling space
debris, Acta Astronaut., 68, 113–120,
<a href="https://doi.org/10.1016/j.actaastro.2010.06.045" target="_blank">https://doi.org/10.1016/j.actaastro.2010.06.045</a>, 2011.
</mixed-citation></ref-html>
<ref-html id="bib1.bib15"><label>15</label><mixed-citation>
Pham, D. B., Kim, J., Lee, S.-G., and Gwak, K.-W.: Double-Loop Control with
Hierarchical Sliding Mode and Proportional Integral Loop for 2D Ridable
Ballbot, Int. J. Precis. Eng. Man., 20, 1519–1532,
<a href="https://doi.org/10.1007/s12541-019-00139-4" target="_blank">https://doi.org/10.1007/s12541-019-00139-4</a>, 2019.
</mixed-citation></ref-html>
<ref-html id="bib1.bib16"><label>16</label><mixed-citation>
Shan, M., Guo, J., and Gill, E.: Review and comparison of active space
debris capturing and removal methods, Prog. Aerosp. Sci., 80,
18–32, <a href="https://doi.org/10.1016/j.paerosci.2015.11.001" target="_blank">https://doi.org/10.1016/j.paerosci.2015.11.001</a>, 2016.
</mixed-citation></ref-html>
<ref-html id="bib1.bib17"><label>17</label><mixed-citation>
Slotine, J. J. and Sastry, S. S.: Tracking control of non-linear systems using sliding surfaces with application to robot manipulators, with application to robot manipulators, International Journal of Control, 38, 465–492, <a href="https://doi.org/10.1080/00207178308933088" target="_blank">https://doi.org/10.1080/00207178308933088</a>, 1983.
</mixed-citation></ref-html>
<ref-html id="bib1.bib18"><label>18</label><mixed-citation>
Wei, C., Liu, H., Tan, C., Liu, Y., and Zhao, Y.: Rapid attitude maneuver of the space tether net capture system using active disturbance rejection control, Mech. Sci., 10, 575–587, <a href="https://doi.org/10.5194/ms-10-575-2019" target="_blank">https://doi.org/10.5194/ms-10-575-2019</a>, 2019.
</mixed-citation></ref-html>
<ref-html id="bib1.bib19"><label>19</label><mixed-citation>
Xu, B., Yang, Y., Yan, Y., and Zhang, B.: Bionics design and dynamics
analysis of space webs based on spider predation, Acta Astronaut., 159,
294–307, <a href="https://doi.org/10.1016/j.actaastro.2019.03.045" target="_blank">https://doi.org/10.1016/j.actaastro.2019.03.045</a>, 2019.
</mixed-citation></ref-html>
<ref-html id="bib1.bib20"><label>20</label><mixed-citation>
Zhang, F., Huang, P., Meng, Z., Zhang, Y., and Liu, Z.: Dynamics analysis
and controller design for maneuverable tethered space net robot, J.
Guid. Control Dynam., 40, 2828–2843, 2017.
</mixed-citation></ref-html>--></article>
