Location: Siebert, Rode, Herzog, Till and Blickhan, 2008 @ ea4d0fbc081b / siebert_CC_2008.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2010-04-07 11:38 +1200
Desc:
Updated model files, added a second model and re-parameterised the original model according to table 2 in the paper.
Permanent Source URI:
http://models.cellml.org/workspace/siebert_rode_herzog_till_blickhan_2008/rawfile/ea4d0fbc081bf1099f9b54fa380487004e6622eb/siebert_CC_2008.cellml

<?xml version="1.0" encoding="utf-8"?>

<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="siebert_2008" name="siebert_2008">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Nonlinearities make a difference: comparison of two common Hill-type models with real muscle</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
          This CellML model represents the [CC] model from the original published paper. Parameter values have been taken from table 2, SOL1 for the [CC] model. The CellML model runs in OpenCell to recreate the published results. The units have been checked and are consistent. The model will not run in COR due to the presence of differential algebraic equations - or "circular arguments".  
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
ABSTRACT: Compared to complex structural Huxley-type models, Hill-type models phenomenologically describe muscle contraction using only few state variables. The Hill-type models dominate in the ever expanding field of musculoskeletal simulations for simplicity and low computational cost. Reasonable parameters are required to gain insight into mechanics of movement. The two most common Hill-type muscle models used contain three components. The series elastic component is connected in series to the contractile component. A parallel elastic component is either connected in parallel to both the contractile and the series elastic component (model [CC+SEC]), or is connected in parallel only with the contractile component (model [CC]). As soon as at least one of the components exhibits substantial nonlinearities, as, e.g., the contractile component by the ability to turn on and off, the two models are mechanically different. We tested which model ([CC+SEC] or [CC]) represents the cat soleus better. Ramp experiments consisting of an isometric and an isokinetic part were performed with an in situ cat soleus preparation using supramaximal nerve stimulation. Hill-type models containing force-length and force-velocity relationship, excitation-contraction coupling and series and parallel elastic force-elongation relations were fitted to the data. To test which model might represent the muscle better, the obtained parameters were compared with experimentally determined parameters. Determined in situations with negligible passive force, the force-velocity relation and the series elastic component relation are independent of the chosen model. In contrast to model [CC+SEC], these relations predicted by model [CC] were in accordance with experimental relations. In conclusion model [CC] seemed to better represent the cat soleus contraction dynamics and should be preferred in the nonlinear regression of muscle parameters and in musculoskeletal modeling. 
</para>

<para>
The original paper is cited below:
</para>

<para>Nonlinearities make a difference: comparison of two common Hill-type models with real muscle, Tobias Siebert, Christian Rode, Walter Herzog, Olaf Till and Reinhard Blickhan, 2008, <emphasis>Biological Cybernetics</emphasis>, 98, 133-143. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/18049823">PubMed ID: 18049823</ulink> 
</para>

<informalfigure float="0" id="figure">
<mediaobject>
  <imageobject>
    <imagedata fileref="siebert_2008.png"/>
  </imageobject>
</mediaobject>
<caption>Two different variants of the Hill muscle model. Isometric force is defined by a piecewise equation.</caption>
</informalfigure>

</sect1>
</article>
</documentation>

  <units name="mm">
     <unit units="meter" prefix="milli"/>
  </units>
  <units name="mm_per_second">
     <unit units="mm"/>
     <unit units="second" exponent="-1"/>
  </units>
  <units name="per_mm">
     <unit units="mm" exponent="-1"/>
  </units>
  <units name="newton_per_mm">
     <unit units="newton"/>
     <unit units="mm" exponent="-1"/>
  </units>
  <units name="mm_newton_per_second">
     <unit units="mm"/>
	 <unit units="newton"/>
     <unit units="second" exponent="-1"/>
  </units>
  
  <component name="environment">
     <variable name="time" units="second" public_interface="out"/>
  </component>
  
  <component name="contraction">
     <variable initial_value="-24" name="L_1" units="mm"/>
     <variable initial_value="-14" name="L_2" units="mm"/>
     <variable initial_value="4" name="L_3" units="mm"/>
     <variable initial_value="14" name="L_4" units="mm"/>
     <variable initial_value="0.51" name="f_c" units="newton"/>
     <variable initial_value="-152" name="v_max" units="mm_per_second"/>
     <variable initial_value="9.5" name="curv" units="dimensionless"/>
     <variable initial_value="0.0064" name="k_1" units="newton"/>
     <variable initial_value="0.414" name="k_2" units="per_mm"/>
     <variable initial_value="8.9" name="F_1" units="newton"/>
     <variable initial_value="2.9" name="d_LSEC1" units="mm"/>
     <variable initial_value="2.5" name="k_sh" units="dimensionless"/>
     <variable initial_value="2.6" name="L_m" units="mm"/>
     <variable initial_value="20.6" name="F_im" units="newton"/>
     <variable initial_value="0.065" name="tau" units="second"/>
     <variable name="v_cc" units="mm_per_second"/>
     <variable name="f_v" units="dimensionless"/>
     <variable name="f_L" units="newton"/>
     <variable name="f_sec" units="newton"/>
     <variable name="time" public_interface="in" units="second"/>
     <variable name="delta_L_sec" units="mm"/>
     <variable initial_value="4.1" name="delta_L_sec1" units="mm"/>
     <variable initial_value="8.2" name="k" units="newton_per_mm"/>
     <variable initial_value="0.2" name="L_cc" units="mm"/>
     <variable name="L_mtc" units="mm"/>
     <variable initial_value="1" name="A" units="dimensionless"/>
     
	 <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply>
		  <eq/>
		  <ci>delta_L_sec</ci>
		  <apply>
		    <minus/>
			<apply>
			  <minus/>
			  <ci>L_mtc</ci>
			  <ci>L_cc</ci>
			</apply>
			<ci>L_m</ci>
		  </apply>
		</apply>
        
		<apply>
           <eq/>
           <ci>f_sec</ci>
           <piecewise>
              <piece>
                 <apply>
                    <times/>
                    <apply>
                       <divide/>
                       <ci>F_1</ci>
                       <apply>
                          <minus/>
                          <apply>
                             <exp/>
                             <ci>k_sh</ci>
                          </apply>
                          <cn cellml:units="dimensionless">1</cn>
                       </apply>
                    </apply>
                    <apply>
                       <minus/>
                       <apply>
                          <exp/>
                          <apply>
                             <divide/>
                             <apply>
                                <times/>
                                <ci>k_sh</ci>
                                <ci>delta_L_sec</ci>
                             </apply>
                             <ci>delta_L_sec1</ci>
                          </apply>
                       </apply>
                       <cn cellml:units="dimensionless">1</cn>
                    </apply>
                 </apply>
                 <apply>
                    <and/>
                    <apply>
                       <gt/>
                       <ci>delta_L_sec</ci>
                       <cn cellml:units="mm">0</cn>
                    </apply>
                    <apply>
                       <lt/>
                       <ci>delta_L_sec</ci>
                       <ci>delta_L_sec1</ci>
                    </apply>
                 </apply>
              </piece>
              <piece>
                 <apply>
                    <plus/>
                    <ci>F_1</ci>
                    <apply>
                       <times/>
                       <ci>k</ci>
                       <apply>
                          <minus/>
                          <ci>delta_L_sec</ci>
                          <ci>delta_L_sec1</ci>
                       </apply>
                    </apply>
                 </apply>
                 <apply>
                    <leq/>
                    <ci>delta_L_sec</ci>
                    <ci>delta_L_sec1</ci>
                 </apply>
              </piece>
           </piecewise>
        </apply>
        
		<apply>
          <eq/>
          <ci>v_cc</ci>
          <apply>
            <times/>
            <apply>
              <divide/>
              <cn cellml:units="mm_newton_per_second">1</cn>
              <ci>f_v</ci>
            </apply>
            <apply>
              <divide/>
              <ci>f_sec</ci>
              <apply>
                <times/>
                <ci>A</ci>
                <ci>f_L</ci>
                <ci>F_im</ci>
              </apply>
            </apply>
          </apply>
        </apply>

        <apply>
           <eq/>
           <apply>
              <diff/>
              <bvar>
                 <ci>time</ci>
              </bvar>
              <ci>L_cc</ci>
           </apply>
           <ci>v_cc</ci>
        </apply>
        
		<apply>
           <eq/>
           <ci>f_v</ci>
           <apply>
              <divide/>
              <apply>
                 <minus/>
                 <ci>v_max</ci>
                 <ci>v_cc</ci>
              </apply>
              <apply>
                 <plus/>
                 <ci>v_max</ci>
                 <apply>
                    <times/>
                    <ci>v_cc</ci>
                    <ci>curv</ci>
                 </apply>
              </apply>
           </apply>
        </apply>
        
		<apply>
           <eq/>
           <ci>f_L</ci>
           <piecewise>
              <piece>
                 <apply>
                    <times/>
                    <apply>
                       <divide/>
                       <ci>f_c</ci>
                       <apply>
                          <minus/>
                          <ci>L_2</ci>
                          <ci>L_1</ci>
                       </apply>
                    </apply>
                    <apply>
                       <minus/>
                       <ci>L_cc</ci>
                       <ci>L_1</ci>
                    </apply>
                 </apply>
                 <apply>
                    <and/>
                    <apply>
                       <geq/>
                       <ci>L_cc</ci>
                       <ci>L_1</ci>
                    </apply>
                    <apply>
                       <leq/>
                       <ci>L_cc</ci>
                       <ci>L_2</ci>
                    </apply>
                 </apply>
              </piece>
              <piece>
                 <apply>
                    <times/>
                    <apply>
                       <divide/>
                       <apply>
                          <minus/>
                          <cn cellml:units="newton">1</cn>
                          <ci>f_c</ci>
                       </apply>
                       <apply>
                          <minus/>
                          <ci>L_2</ci>
                       </apply>
                    </apply>
                    <apply>
                       <minus/>
                       <ci>L_cc</ci>
                       <ci>L_2</ci>
                    </apply>
                 </apply>
                 <apply>
                    <and/>
                    <apply>
                       <gt/>
                       <ci>L_cc</ci>
                       <ci>L_2</ci>
                    </apply>
                    <apply>
                       <leq/>
                       <ci>L_cc</ci>
                       <cn cellml:units="mm">0</cn>
                    </apply>
                 </apply>
              </piece>
              <piece>
                 <cn cellml:units="newton">1</cn>
                 <apply>
                    <and/>
                    <apply>
                       <gt/>
                       <ci>L_cc</ci>
                       <cn cellml:units="mm">0</cn>
                    </apply>
                    <apply>
                       <leq/>
                       <ci>L_cc</ci>
                       <ci>L_3</ci>
                    </apply>
                 </apply>
              </piece>
              <piece>
                 <apply>
                    <times/>
                    <apply>
                       <divide/>
                       <apply>
                          <minus/>
                          <cn cellml:units="newton">1</cn>
                       </apply>
                       <apply>
                          <minus/>
                          <ci>L_4</ci>
                          <ci>L_3</ci>
                       </apply>
                    </apply>
                    <apply>
                       <minus/>
                       <ci>L_cc</ci>
                       <ci>L_3</ci>
                    </apply>
                 </apply>
                 <apply>
                    <and/>
                    <apply>
                       <gt/>
                       <ci>L_cc</ci>
                       <ci>L_3</ci>
                    </apply>
                    <apply>
                       <leq/>
                       <ci>L_cc</ci>
                       <ci>L_4</ci>
                    </apply>
                 </apply>
              </piece>
           </piecewise>
        </apply>
        
		<apply>
           <eq/>
           <ci>L_mtc</ci>
           <piecewise>
              <piece>
                 <cn cellml:units="mm">0.29</cn>
                 <apply>
                    <leq/>
                    <ci>time</ci>
                    <cn cellml:units="second">1</cn>
                 </apply>
              </piece>
              <piece>
                 <cn cellml:units="mm">0.22</cn>
                 <apply>
                    <and/>
                    <apply>
					<gt/>
                       <ci>time</ci>
                       <cn cellml:units="second">1</cn>
                    </apply>
                    <apply>
                       <lt/>
                       <ci>time</ci>
                       <cn cellml:units="second">5</cn>
                    </apply>
                 </apply>
              </piece>
              <otherwise>
                 <cn cellml:units="mm">0.19</cn>
              </otherwise>
           </piecewise>
        </apply>
     </math>
  </component>
  
  <connection>
     <map_components component_1="environment" component_2="contraction"/>
     <map_variables variable_1="time" variable_2="time"/>
  </connection>
  
  <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
    
	<rdf:Description rdf:about="">
      <dc:creator rdf:parseType="Resource">
        <vCard:N rdf:parseType="Resource">
          <vCard:Family>Harrington</vCard:Family>
          <vCard:Given>Paul</vCard:Given>
		</vCard:N>
        <vCard:EMAIL rdf:parseType="Resource">
          <rdf:value>paul.harrington@auckland.ac.nz</rdf:value>
          <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
        </vCard:EMAIL>
        <vCard:ORG rdf:parseType="Resource">
          <vCard:Orgname>The University of Auckland</vCard:Orgname>
          <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
        </vCard:ORG>
      </dc:creator>
  
      <dcterms:created rdf:parseType="Resource">
        <dcterms:W3CDTF>2010-03-26</dcterms:W3CDTF>
      </dcterms:created>
    </rdf:Description>
    
    <rdf:Description rdf:about="#siebert_2008">
      <cmeta:bio_entity>skeletal muscle</cmeta:bio_entity>
	  
      <bqs:reference rdf:parseType="Resource">
        <dc:subject rdf:parseType="Resource">
          <bqs:subject_type>keyword</bqs:subject_type>
          <rdf:value>
            <rdf:Bag>
              <rdf:li>mechanical constitutive laws</rdf:li>
			  <rdf:li>skeletal muscle</rdf:li>
			  <rdf:li>muscle contraction</rdf:li>
            </rdf:Bag>
          </rdf:value>
        </dc:subject>
      </bqs:reference>

      <bqs:reference rdf:parseType="Resource">
        <bqs:Pubmed_id>18049823</bqs:Pubmed_id>
        <bqs:JournalArticle rdf:parseType="Resource">
          <dc:creator>
            <rdf:Seq>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Siebert</vCard:Family>
                    <vCard:Given>T</vCard:Given>
			      </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Rode</vCard:Family>
                    <vCard:Given>C</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Herzog</vCard:Family>
                    <vCard:Given>W</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Till</vCard:Family>
                    <vCard:Given>O</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Blickhan</vCard:Family>
                    <vCard:Given>R</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
            </rdf:Seq>
          </dc:creator>
          <dc:title>
            Nonlinearities make a difference: comparison of two common Hill-type models with real muscle
          </dc:title>
          <dcterms:issued rdf:parseType="Resource">
            <dcterms:W3CDTF>2008-02</dcterms:W3CDTF>
          </dcterms:issued>
          <bqs:Journal rdf:parseType="Resource">
            <dc:title>Biological Cybernetics</dc:title>
          </bqs:Journal>
		  <bqs:volume>98</bqs:volume>
          <bqs:first_page>133</bqs:first_page>
          <bqs:last_page>143</bqs:last_page>
        </bqs:JournalArticle>
      </bqs:reference>
    </rdf:Description>
  </rdf:RDF>
</model>