- Author:
- Catherine Lloyd <c.lloyd@auckland.ac.nz>
- Date:
- 2010-08-06 02:01:44+12:00
- Desc:
- Updated documentation and curation comment.
- Permanent Source URI:
- https://models.cellml.org/workspace/oomens_maenhout_vanoijen_drost_baaijens_2003/rawfile/25a6ffcede519bae774b567da710dc6b64d7542f/oomens_maenhout_vanoijen_drost_baaijens_2003.cellml
<?xml version='1.0' encoding='utf-8'?>
<!-- FILE : oomens_model_2003.xml
CREATED : 8th September 2004
LAST MODIFIED : 8th September 2004
AUTHOR : Catherine Lloyd
The Bioengineering Institute
The University of Auckland
MODEL STATUS : This model conforms to the CellML 1.0 Specification released on
10th August 2001, and the 16/01/2002 CellML Metadata 1.0 Specification.
DESCRIPTION : This file contains a CellML description of Oomens et al.'s 2003 finite element model of contracting skeletal muscle.
CHANGES:
--><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="oomens_2003" name="oomens_2003">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Finite Element Modelling of Contracting Skeletal 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 cannot be run in either OpenCell or COR due to the presence of the MathML operator "sum" which is not yet supported by the CellML API.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
ABSTRACT: To describe the mechanical behaviour of biological tissues and transport processes in biological tissues, conservation laws such as conservation of mass, momentum and energy play a central role. Mathematically these are cast into the form of partial differential equations. Because of nonlinear material behaviour, inhomogeneous properties and usually a complex geometry, it is impossible to find closed-form analytical solutions for these sets of equations. The objective of the finite element method is to find approximate solutions for these problems. The concepts of the finite element method are explained on a finite element continuum model of skeletal muscle. In this case, the momentum equations have to be solved with an extra constraint, because the material behaves as nearly incompressible. The material behaviour consists of a highly nonlinear passive part and an active part. The latter is described with a two-state Huxley model. This means that an extra nonlinear partial differential equation has to be solved. The problems and solutions involved with this procedure are explained. The model is used to describe the mechanical behaviour of a tibialis anterior of a rat. The results have been compared with experimentally determined strains at the surface of the muscle. Qualitatively there is good agreement between measured and calculated strains, but the measured strains were higher.
</para>
<para>
The original paper reference is cited below:
</para>
<para>
Finite element modelling of contracting skeletal muscle, C.W.J. Oomens, M. Maenhout, C.H. van Oijen, M.R. Drost, and F.P. Baaijens, 2003, <emphasis>Phil. Trans. R. Soc. Lond. B </emphasis>, 358, 1453-1460. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/14561336">PubMed ID: 14561336</ulink>
</para>
<para>
Muscle Activation and Contraction: Constitutive Relations Based Directly on Cross-Bridge Kinetics, George I. Zahalak and Shi-Ping Ma, 1990, <emphasis>Journal of Biomedical Engineering</emphasis>, 122, 52-62. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&list_uids=2308304&dopt=Abstract">PubMed ID: 2308304</ulink>
</para>
</sect1>
</article>
</documentation>
<units name="micromolar">
<unit units="mole" prefix="micro"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="millisecond">
<unit units="second" prefix="milli"/>
</units>
<units name="flux">
<unit units="micromolar"/>
<unit units="second" exponent="-1"/>
</units>
<units name="first_order_rate_constant">
<unit units="second" exponent="-1"/>
</units>
<units name="second_order_rate_constant">
<unit units="micromolar" exponent="-1"/>
<unit units="second" exponent="-1"/>
</units>
<units name="micromolar2_per_second">
<unit units="micromolar" exponent="2"/>
<unit units="second" exponent="-1"/>
</units>
<component name="environment">
<variable units="second" public_interface="out" name="time"/>
</component>
<component cmeta:id="r" name="r">
<variable units="dimensionless" public_interface="out" name="r"/>
<variable units="dimensionless" public_interface="in" name="c"/>
<variable units="dimensionless" public_interface="in" name="mu"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>r</ci>
<apply>
<divide/>
<apply>
<power/>
<ci> c </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci> c </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<times/>
<ci> mu </ci>
<ci> c </ci>
</apply>
<apply>
<power/>
<ci> mu </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="Ct" name="Ct">
<variable units="micromolar" public_interface="out" name="Ct"/>
<variable units="dimensionless" name="Ct_initial"/>
<variable units="dimensionless" name="b"/>
<variable units="dimensionless" name="km" initial_value="0.006"/>
<variable units="millisecond" name="tau0" initial_value="9.2"/>
<variable units="dimensionless" name="Q0"/>
<variable units="micromolar" public_interface="in" name="c_"/>
<variable units="dimensionless" public_interface="in" name="rho"/>
<variable units="dimensionless" public_interface="in" name="c"/>
<variable units="dimensionless" public_interface="in" name="r"/>
<variable units="dimensionless" public_interface="in" name="mu"/>
<variable units="dimensionless" public_interface="in" name="chi"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ct</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> rho </ci>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> c </ci>
<ci> c_ </ci>
</apply>
</apply>
<ci> chi </ci>
</apply>
<apply>
<times/>
<apply>
<power/>
<ci> tau0 </ci>
<cn cellml:units="dimensionless"> -1.0 </cn>
</apply>
<apply>
<divide/>
<ci> c </ci>
<apply>
<plus/>
<ci> c </ci>
<ci> km </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ct_initial</ci>
<apply>
<plus/>
<ci> c </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> b </ci>
<ci> Q0 </ci>
</apply>
<apply>
<times/>
<ci> r </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<apply>
<divide/>
<ci> mu </ci>
<ci> c </ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<times/>
<ci> b </ci>
<ci> Q0 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="Ca_in" name="Ca_in">
<variable units="flux" public_interface="out" name="Ca_in"/>
<variable units="dimensionless" public_interface="in" name="rho"/>
<variable units="dimensionless" public_interface="in" name="chi"/>
<variable units="micromolar" public_interface="in" name="c_"/>
<variable units="dimensionless" public_interface="in" name="c"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_in</ci>
</apply>
<apply>
<times/>
<ci> rho </ci>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<apply>
<divide/>
<ci> c </ci>
<ci> c_ </ci>
</apply>
</apply>
<ci> chi </ci>
</apply>
</apply>
</math>
</component>
<component cmeta:id="chi" name="chi">
<variable units="dimensionless" public_interface="out" name="chi"/>
<variable units="dimensionless" name="chi_"/>
<variable units="dimensionless" name="i"/>
<variable units="second" name="time_i"/>
<variable units="second" name="tau1" initial_value="0.005"/>
<variable units="second" name="tau2" initial_value="0.001"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>chi</ci>
<apply>
<sum/>
<domainofapplication>
<ci type="set"> i </ci>
</domainofapplication>
<apply>
<times/>
<ci> chi_ </ci>
<apply>
<minus/>
<ci> time </ci>
<ci> time_i </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>chi_</ci>
<apply>
<times/>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> time </ci>
</apply>
<ci> tau1 </ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> time </ci>
</apply>
<ci> tau2 </ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci> tau1 </ci>
<ci> tau2 </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="model_parameters">
<variable units="dimensionless" public_interface="out" name="mu" initial_value="0.2"/>
<variable units="dimensionless" public_interface="out" name="c"/>
<variable units="micromolar" public_interface="out" name="c_"/>
<variable units="dimensionless" public_interface="out" name="rho" initial_value="2.9"/>
</component>
<connection>
<map_components component_2="environment" component_1="Ct"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="chi"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="Ca_in"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="r"/>
<map_variables variable_2="c" variable_1="c"/>
<map_variables variable_2="mu" variable_1="mu"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="Ct"/>
<map_variables variable_2="c" variable_1="c"/>
<map_variables variable_2="mu" variable_1="mu"/>
<map_variables variable_2="rho" variable_1="rho"/>
<map_variables variable_2="c_" variable_1="c_"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="Ca_in"/>
<map_variables variable_2="c" variable_1="c"/>
<map_variables variable_2="c_" variable_1="c_"/>
<map_variables variable_2="rho" variable_1="rho"/>
</connection>
<connection>
<map_components component_2="r" component_1="Ct"/>
<map_variables variable_2="r" variable_1="r"/>
</connection>
<connection>
<map_components component_2="chi" component_1="Ct"/>
<map_variables variable_2="chi" variable_1="chi"/>
</connection>
<connection>
<map_components component_2="chi" component_1="Ca_in"/>
<map_variables variable_2="chi" variable_1="chi"/>
</connection>
<rdf:RDF>
<rdf:Bag rdf:about="rdf:#0cecd8d8-b025-4068-acbf-dbf847408519">
<rdf:li>myofilament mechanics</rdf:li>
<rdf:li>skeletal muscle</rdf:li>
<rdf:li>myosin</rdf:li>
<rdf:li>calcium</rdf:li>
<rdf:li>finite element</rdf:li>
<rdf:li>muscle contraction</rdf:li>
<rdf:li>mechanical constitutive laws</rdf:li>
</rdf:Bag>
<rdf:Seq rdf:about="rdf:#7bc88b18-413c-422b-8f9c-8a23af933ca2">
<rdf:li rdf:resource="rdf:#15bc880e-47a5-4214-8e44-04906428d188"/>
<rdf:li rdf:resource="rdf:#0f173c3c-bf52-4040-9ce2-f21cb4d4aceb"/>
<rdf:li rdf:resource="rdf:#ee14d508-e6b4-4fd7-90a0-4dd5bc35651d"/>
<rdf:li rdf:resource="rdf:#7558b790-0721-4471-92ca-80e2b52a19a6"/>
<rdf:li rdf:resource="rdf:#a45336c9-54fe-48a3-8b4d-8c2d86edb512"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#ee14d508-e6b4-4fd7-90a0-4dd5bc35651d">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#0a65125e-ac03-4cae-a044-2d2eadbc8e86"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#132f745a-d1f8-4bce-b572-2090ca84ef9e">
<dcterms:W3CDTF>2003-09-29</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:publisher>
The University of Auckland, Auckland Bioengineering Institute
</dc:publisher>
<dcterms:created rdf:resource="rdf:#af16453c-8231-4d86-9ef5-dae86f85142c"/>
<dc:creator rdf:resource="rdf:#fefcb4d1-1ed2-4363-85ab-f57edc684ab5"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#36c1aa74-a6e4-4b28-b441-0fa2970ef365">
<vCard:Given>C</vCard:Given>
<vCard:Family>Oomens</vCard:Family>
<vCard:Other>J</vCard:Other>
<vCard:Other>W</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#bda8451b-6dcc-4866-958a-ff57528c2f1b">
<dc:creator rdf:resource="rdf:#7bc88b18-413c-422b-8f9c-8a23af933ca2"/>
<dc:title>
Finite element modelling of contracting skeletal muscle
</dc:title>
<bqs:volume>358</bqs:volume>
<bqs:first_page>1453</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#7f942b46-56bf-4ebb-ab0f-ab66b05262c5"/>
<dcterms:issued rdf:resource="rdf:#132f745a-d1f8-4bce-b572-2090ca84ef9e"/>
<bqs:last_page>1460</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#af16453c-8231-4d86-9ef5-dae86f85142c">
<dcterms:W3CDTF>2004-09-08</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#df9a66e5-aa84-4005-803f-964f98a2ed70">
<dc:subject rdf:resource="rdf:#608a9111-9eb7-4a16-9696-a77ac986ac83"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d1425af8-eba5-455f-9d42-e075e3164209">
<vCard:Given>F</vCard:Given>
<vCard:Family>Baaijens</vCard:Family>
<vCard:Other>P</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="#oomens_2003">
<dc:title>
Oomens et al.'s 2003 finite element model of contracting skeletal muscle.
</dc:title>
<cmeta:comment rdf:resource="rdf:#c960ec88-1716-4878-b666-2a611863212f"/>
<bqs:reference rdf:resource="rdf:#8b90fd97-869f-4168-a032-147107170c06"/>
<bqs:reference rdf:resource="rdf:#df9a66e5-aa84-4005-803f-964f98a2ed70"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#801c103a-1a28-4ec6-b16f-bbfc6d43946f">
<vCard:Given>M</vCard:Given>
<vCard:Family>Maenhout</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0f173c3c-bf52-4040-9ce2-f21cb4d4aceb">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#801c103a-1a28-4ec6-b16f-bbfc6d43946f"/>
</rdf:Description>
<rdf:Description rdf:about="#Ct">
<dcterms:alternative>total calcium in the myofibrillar space</dcterms:alternative>
<dc:title>Ct</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7f942b46-56bf-4ebb-ab0f-ab66b05262c5">
<dc:title>Phil. Trans. R. Soc. Lond. B</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a45336c9-54fe-48a3-8b4d-8c2d86edb512">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#d1425af8-eba5-455f-9d42-e075e3164209"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0a65125e-ac03-4cae-a044-2d2eadbc8e86">
<vCard:Given>C</vCard:Given>
<vCard:Family>van Oijen</vCard:Family>
<vCard:Other>H</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="#r">
<dcterms:alternative>fraction of activated actin</dcterms:alternative>
<dc:title>r</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e8c42e4f-eaf3-45bf-b763-a4e833dac39c">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#19fe7ae6-e427-4a12-967e-53f0effbd3a1">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a4a51e86-c2f4-4014-8e59-1c700a3f3d37">
<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
<rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#c960ec88-1716-4878-b666-2a611863212f">
<dc:creator rdf:resource="rdf:#e8c42e4f-eaf3-45bf-b763-a4e833dac39c"/>
<rdf:value>
This is the CellML description of Oomens et al.'s 2003 finite element model of contracting skeletal muscle.
</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="#chi">
<dcterms:alternative>Sum of the normalised impulse functions</dcterms:alternative>
<dc:title>chi</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#15bc880e-47a5-4214-8e44-04906428d188">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#36c1aa74-a6e4-4b28-b441-0fa2970ef365"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#608a9111-9eb7-4a16-9696-a77ac986ac83">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value rdf:resource="rdf:#0cecd8d8-b025-4068-acbf-dbf847408519"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7558b790-0721-4471-92ca-80e2b52a19a6">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#6eb797b8-1d0a-401c-87f1-a394ef210116"/>
</rdf:Description>
<rdf:Description rdf:about="#Ca_in">
<dcterms:alternative>rate of injection of calcium from the SR</dcterms:alternative>
<dc:title>Ca_in</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#6eb797b8-1d0a-401c-87f1-a394ef210116">
<vCard:Given>M</vCard:Given>
<vCard:Family>Drost</vCard:Family>
<vCard:Other>R</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#85734511-0e70-4745-b5c6-ff51aabefdc7">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fefcb4d1-1ed2-4363-85ab-f57edc684ab5">
<vCard:ORG rdf:resource="rdf:#85734511-0e70-4745-b5c6-ff51aabefdc7"/>
<vCard:EMAIL rdf:resource="rdf:#a4a51e86-c2f4-4014-8e59-1c700a3f3d37"/>
<vCard:N rdf:resource="rdf:#19fe7ae6-e427-4a12-967e-53f0effbd3a1"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#8b90fd97-869f-4168-a032-147107170c06">
<bqs:Pubmed_id>14561336</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#bda8451b-6dcc-4866-958a-ff57528c2f1b"/>
</rdf:Description>
</rdf:RDF>
</model>