- Author:
- Dewan Sarwar <sarwarcse@gmail.com>
- Date:
- 2019-01-22 18:17:02+13:00
- Desc:
- adding UTF-8
- Permanent Source URI:
- http://models.cellml.org/workspace/584/rawfile/771f7cf8dd59d9a21ece87e7966a93d064a6b01a/siebert_CC_2008.cellml
<?xml version="1.0" encoding="UTF-8"?>
<model xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:semsim="http://www.bhi.washington.edu/SemSim#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" name="siebert_2008" cmeta:id="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="per_mm">
<unit exponent="-1.0" units="mm" />
</units>
<units name="newton_per_mm">
<unit exponent="-1.0" units="mm" />
<unit units="newton" />
</units>
<units name="mm_per_second">
<unit exponent="-1.0" units="second" />
<unit units="mm" />
</units>
<units name="mm">
<unit prefix="milli" units="meter" />
</units>
<component name="contraction">
<variable cmeta:id="contraction.L_cc" name="L_cc" units="mm" />
<variable cmeta:id="contraction.A" name="A" units="dimensionless" />
<variable cmeta:id="contraction.k_2" name="k_2" units="per_mm" />
<variable cmeta:id="contraction.v_cc" name="v_cc" units="mm_per_second" />
<variable cmeta:id="contraction.f_v" name="f_v" units="dimensionless" />
<variable cmeta:id="contraction.k_1" name="k_1" units="newton" />
<variable cmeta:id="contraction.d_LSEC1" name="d_LSEC1" units="mm" />
<variable cmeta:id="contraction.L_mtc" name="L_mtc" units="mm" />
<variable cmeta:id="contraction.time" name="time" public_interface="in" units="second" />
<variable cmeta:id="contraction.f_sec" name="f_sec" units="newton" />
<variable cmeta:id="contraction.k" name="k" units="newton_per_mm" />
<variable cmeta:id="contraction.v_max" name="v_max" units="mm_per_second" />
<variable cmeta:id="contraction.f_c" name="f_c" units="newton" />
<variable cmeta:id="contraction.L_3" name="L_3" units="mm" />
<variable cmeta:id="contraction.L_4" name="L_4" units="mm" />
<variable cmeta:id="contraction.L_1" name="L_1" units="mm" />
<variable cmeta:id="contraction.L_2" name="L_2" units="mm" />
<variable cmeta:id="contraction.F_1" name="F_1" units="newton" />
<variable cmeta:id="contraction.k_sh" name="k_sh" units="dimensionless" />
<variable cmeta:id="contraction.L_m" name="L_m" units="mm" />
<variable cmeta:id="contraction.delta_L_sec1" name="delta_L_sec1" units="mm" />
<variable cmeta:id="contraction.F_im" name="F_im" units="newton" />
<variable cmeta:id="contraction.f_L" name="f_L" units="newton" />
<variable cmeta:id="contraction.curv" name="curv" units="dimensionless" />
<variable cmeta:id="contraction.tau" name="tau" units="second" />
<variable cmeta:id="contraction.delta_L_sec" name="delta_L_sec" units="mm" />
<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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" 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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<and />
<apply>
<gt />
<ci>delta_L_sec</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" 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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" 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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" 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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mm">0</cn>
</apply>
</apply>
</piece>
<piece>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="newton">1</cn>
<apply>
<and />
<apply>
<gt />
<ci>L_cc</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" 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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" 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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mm">0.29</cn>
<apply>
<leq />
<ci>time</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
</apply>
</piece>
<piece>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mm">0.22</cn>
<apply>
<and />
<apply>
<gt />
<ci>time</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
</apply>
<apply>
<lt />
<ci>time</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">5</cn>
</apply>
</apply>
</piece>
<otherwise>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mm">0.19</cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="environment">
<variable cmeta:id="environment.time" name="time" public_interface="out" units="second" />
</component>
<connection>
<map_components component_1="environment" component_2="contraction" />
<map_variables variable_1="time" variable_2="time" />
</connection>
<rdf:RDF xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:bqbiol="http://biomodels.net/biology-qualifiers/" xmlns:ro="http://www.obofoundry.org/ro/ro.owl#" xmlns:opb="http://bhi.washington.edu/OPB#" xmlns:j.0="http://imc.org/vCard/3.0#">
<rdf:Description rdf:about="#contraction.k_2">
<dcterms:description>Curve-shaping parameter for setting the parallel elastic component's force-velocity relationship</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.v_max">
<dcterms:description>Maximal contractile component shortening velocity</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.v_cc">
<semsim:isComputationalComponentFor>
<rdf:Description rdf:about="#property_0">
<semsim:physicalPropertyOf>
<rdf:Description rdf:about="#process_0">
<dcterms:description>Contraction of the modeled contractile element</dcterms:description>
<semsim:name>Contraction of contractile element</semsim:name>
</rdf:Description>
</semsim:physicalPropertyOf>
<semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_01658" />
</rdf:Description>
</semsim:isComputationalComponentFor>
<dcterms:description>Velocity of contractile component</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.tau">
<dcterms:description>Time constant for muscle activation</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.L_2">
<dcterms:description>Muscle length at which the ascending limb changes slope</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="http://www.cellml.org/cellml/1.0">
<dc:title>Nonlinearities make a difference: comparison of two common Hill-type models with real muscle (CC Model)</dc:title>
<dc:creator rdf:parseType="Resource">
<vCard:N vCard:Family="Harrington" vCard:Given="Paul" />
<vCard:EMAIL rdf:value="paul.harrington@auckland.ac.nz" rdf:type="http://imc.org/vCard/3.0#internet" />
<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 dcterms:W3CDTF="2010-03-26" />
</rdf:Description>
<rdf:Description rdf:about="#contraction.k">
<dcterms:description>Scaling parameter for the series elastic component's force-elongation relationship</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.k_1">
<dcterms:description>Curve-shaping parameter for setting the parallel elastic component's force-velocity relationship</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.f_v">
<dcterms:description>Velocity-based force factor in overall contractile component force</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.L_mtc">
<dcterms:description>Length of the muscle-tendon complex</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.L_1">
<dcterms:description>Curve-shaping parameter determining width of muscle force-length relationship</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.f_L">
<dcterms:description>Length-based force factor in overall contractile component force</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.L_cc">
<dcterms:description>The contractile component length minus the optimal contractile component length</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.delta_L_sec1">
<dcterms:description>Length at which the series elastic component force-elongation relation changes from exponential to linear</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#environment.time">
<semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_01023" />
<dcterms:description>Time solution domain</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.f_sec">
<dcterms:description>Force generated in series elastic component</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.d_LSEC1">
<dcterms:description>The length at which the force�elongation relation changes from exponential to linear</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.k_sh">
<dcterms:description>Shaping parameter for the series elastic component's force-elongation relationship</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.F_1">
<dcterms:description>Force at which the series elastic component force-elongation relation changes from exponential to linear</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.L_m">
<dcterms:description>The constant sum of the slack length of the series elastic component and the optimal length of the contractile component (used in the [CC] model, but not the [CC+SEC] model)</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.L_4">
<dcterms:description>Curve-shaping parameter determining width of muscle force-length relationship</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.F_im">
<dcterms:description>The maximum active isometric force of the muscle</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.curv">
<dcterms:description>Curve-shaping parameter determining relationship between f_v and v_cc</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.time">
<semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_01023" />
<dcterms:description>Time solution domain</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.f_c">
<dcterms:description>The force at which the ascending limb changes slope</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="http://www.cellml.org/cellml/1.0#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 vCard:Family="Siebert" vCard:Given="T" />
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N vCard:Family="Rode" vCard:Given="C" />
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N vCard:Family="Herzog" vCard:Given="W" />
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N vCard:Family="Till" vCard:Given="O" />
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N vCard:Family="Blickhan" vCard:Given="R" />
</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 dcterms:W3CDTF="2008-02" />
<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:Description rdf:about="#contraction.delta_L_sec">
<dcterms:description>Elongation of series elastic component</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#siebert_CC_2008" semsim:modelId="siebert_2008" semsim:modelName="siebert_2008" />
<rdf:Description rdf:about="#contraction.A">
<dcterms:description>Activation state</dcterms:description>
</rdf:Description>
<rdf:Description rdf:about="#contraction.L_3">
<dcterms:description>Curve-shaping parameter that sets the length of the force-length plateau</dcterms:description>
</rdf:Description>
</rdf:RDF>
</model>
