- Author:
- David Nickerson <nickerso@users.sourceforge.net>
- Date:
- 2015-11-09 11:56:14+13:00
- Desc:
- adding a version of the Guccione model from Jagir Hussan, mathematically equivalent but different parameter names and values
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/guccione_mcculloch_waldman_1991/rawfile/6066045f96649bad8274c8e044b61eae5b137e36/new-guccione.xml
<?xml version="1.0" encoding="utf-8"?>
<model name="Guccione" cmeta:id="Guccione" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" 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>Razak Jainulabdeen</vCard:Family>
<vCard:Given>Jagir Hussan</vCard:Given>
</vCard:N>
<vCard:EMAIL rdf:parseType="Resource">
<rdf:value>r.jagir@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>The Bioengineering Institute</vCard:Orgunit>
</vCard:ORG>
</dc:creator>
<dcterms:created rdf:parseType="Resource">
<dcterms:W3CDTF>2014-10-2</dcterms:W3CDTF>
</dcterms:created>
</rdf:Description>
<rdf:Description rdf:about="#PoleZero">
<cmeta:comment rdf:parseType="Resource">
<rdf:value>
This is a CellML version of the Guccione constitutive material law,
defining the relation between the six independent strain components
and the stress components. It is assumed that the strain components
will be controlled externally by the application using this CellML
model. The initial parameters are for opening angle 0.
</rdf:value>
</cmeta:comment>
<bqs:reference>
<bqs:Pubmed_id>pubmed_id</bqs:Pubmed_id>
<bqs:JournalArticle rdf:parseType="Resource">
<dc:title>
Passive material properties of intact ventricular myocardium determined from a cylindrical model
</dc:title>
<dcterms:issued rdf:parseType="Resource">
<dcterms:W3CDTF>1991</dcterms:W3CDTF>
</dcterms:issued>
<bqs:Journal rdf:parseType="Resource">
<dc:title>J Biomech Eng.</dc:title>
</bqs:Journal>
<bqs:volume>113</bqs:volume>
<bqs:first_page>1</bqs:first_page>
<bqs:last_page>55</bqs:last_page>
</bqs:JournalArticle>
</bqs:reference>
</rdf:Description>
</rdf:RDF>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<rdf:Description rdf:about="#interface">
<cmeta:comment rdf:parseType="Resource">
<rdf:value>
We'll use this component as the "interface" to the model, all
other components are hidden via encapsulation in this component.
</rdf:value>
</cmeta:comment>
</rdf:Description>
</rdf:RDF>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<rdf:Description rdf:about="#equations">
<cmeta:comment rdf:parseType="Resource">
<rdf:value>
In this simple model we only have one component, which holds the
six equations.
</rdf:value>
</cmeta:comment>
</rdf:Description>
</rdf:RDF> <units name="strain">
<unit units="dimensionless"/>
</units>
<units name="stress">
<unit units="dimensionless"/>
</units>
<component name="interface">
<variable name="E11" units="strain" private_interface="out" initial_value="0.0"/>
<variable name="E22" units="strain" private_interface="out" initial_value="0.0"/>
<variable name="E33" units="strain" private_interface="out" initial_value="0.0"/>
<variable name="E12" units="strain" private_interface="out" initial_value="0.0"/>
<variable name="E13" units="strain" private_interface="out" initial_value="0.0"/>
<variable name="E23" units="strain" private_interface="out" initial_value="0.0"/>
<variable name="C" units="strain" private_interface="out" initial_value="0.641"/>
<variable name="b1" units="strain" private_interface="out" initial_value="2.929"/>
<variable name="b2" units="strain" private_interface="out" initial_value="10.82"/>
<variable name="b3" units="strain" private_interface="out" initial_value="0.0"/>
<variable name="b4" units="strain" private_interface="out" initial_value="14.08"/>
<variable name="Tdev11" units="stress" public_interface="out" private_interface="in"/>
<variable name="Tdev22" units="stress" public_interface="out" private_interface="in"/>
<variable name="Tdev33" units="stress" public_interface="out" private_interface="in"/>
<variable name="Tdev12" units="stress" public_interface="out" private_interface="in"/>
<variable name="Tdev13" units="stress" public_interface="out" private_interface="in"/>
<variable name="Tdev23" units="stress" public_interface="out" private_interface="in"/>
</component>
<component name="equations">
<variable name="E11" units="strain" public_interface="in"/>
<variable name="E22" units="strain" public_interface="in"/>
<variable name="E33" units="strain" public_interface="in"/>
<variable name="E12" units="strain" public_interface="in"/>
<variable name="E13" units="strain" public_interface="in"/>
<variable name="E23" units="strain" public_interface="in"/>
<variable name="C" units="strain" public_interface="in"/>
<variable name="b1" units="strain" public_interface="in"/>
<variable name="b2" units="strain" public_interface="in"/>
<variable name="b3" units="strain" public_interface="in"/>
<variable name="b4" units="strain" public_interface="in"/>
<variable name="Tdev11" units="stress" public_interface="out"/>
<variable name="Tdev22" units="stress" public_interface="out"/>
<variable name="Tdev33" units="stress" public_interface="out"/>
<variable name="Tdev12" units="stress" public_interface="out"/>
<variable name="Tdev13" units="stress" public_interface="out"/>
<variable name="Tdev23" units="stress" public_interface="out"/>
<variable name="Q" units="strain"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Q</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="strain">2</cn>
<ci>b1</ci>
<apply>
<plus/>
<ci>E11</ci>
<ci>E22</ci>
<ci>E33</ci>
</apply>
</apply>
<apply>
<times/>
<ci>b2</ci>
<apply>
<power/>
<ci>E11</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<ci>b3</ci>
<apply>
<plus/>
<apply>
<power/>
<ci>E33</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>E22</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<apply>
<power/>
<ci>E23</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<cn cellml:units="strain">2</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<!-- Since symetric shear is modelled -->
<cn cellml:units="strain">2</cn>
<ci>b4</ci>
<apply>
<plus/>
<apply>
<power/>
<ci>E13</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>E12</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Tdev11</ci>
<apply>
<times/>
<cn cellml:units="strain">0.5</cn>
<ci>C</ci>
<apply>
<exp/>
<ci>Q</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="strain">2</cn>
<ci>b1</ci>
</apply>
<apply>
<times/>
<cn cellml:units="strain">2</cn>
<ci>b2</ci>
<ci>E11</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Tdev22</ci>
<apply>
<times/>
<cn cellml:units="strain">0.5</cn>
<ci>C</ci>
<apply>
<exp/>
<ci>Q</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="strain">2</cn>
<ci>b1</ci>
</apply>
<apply>
<times/>
<cn cellml:units="strain">2</cn>
<ci>b3</ci>
<ci>E22</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Tdev33</ci>
<apply>
<times/>
<cn cellml:units="strain">0.5</cn>
<ci>C</ci>
<apply>
<exp/>
<ci>Q</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="strain">2</cn>
<ci>b1</ci>
</apply>
<apply>
<times/>
<cn cellml:units="strain">2</cn>
<ci>b3</ci>
<ci>E33</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Tdev12</ci>
<apply>
<times/>
<cn cellml:units="strain">0.5</cn>
<ci>C</ci>
<apply>
<exp/>
<ci>Q</ci>
</apply>
<cn cellml:units="strain">2</cn>
<ci>b4</ci>
<ci>E12</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Tdev13</ci>
<apply>
<times/>
<cn cellml:units="strain">0.5</cn>
<ci>C</ci>
<apply>
<exp/>
<ci>Q</ci>
</apply>
<cn cellml:units="strain">2</cn>
<ci>b4</ci>
<ci>E13</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Tdev23</ci>
<apply>
<times/>
<cn cellml:units="strain">0.5</cn>
<ci>C</ci>
<apply>
<exp/>
<ci>Q</ci>
</apply>
<cn cellml:units="strain">2</cn>
<ci>b4</ci>
<ci>E23</ci>
</apply>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="interface">
<component_ref component="equations"/>
</component_ref>
</group>
<connection>
<map_components component_1="interface" component_2="equations"/>
<map_variables variable_1="E11" variable_2="E11"/>
<map_variables variable_1="E22" variable_2="E22"/>
<map_variables variable_1="E33" variable_2="E33"/>
<map_variables variable_1="E12" variable_2="E12"/>
<map_variables variable_1="E13" variable_2="E13"/>
<map_variables variable_1="E23" variable_2="E23"/>
<map_variables variable_1="Tdev11" variable_2="Tdev11"/>
<map_variables variable_1="Tdev22" variable_2="Tdev22"/>
<map_variables variable_1="Tdev33" variable_2="Tdev33"/>
<map_variables variable_1="Tdev12" variable_2="Tdev12"/>
<map_variables variable_1="Tdev13" variable_2="Tdev13"/>
<map_variables variable_1="Tdev23" variable_2="Tdev23"/>
<map_variables variable_1="C" variable_2="C"/>
<map_variables variable_1="b1" variable_2="b1"/>
<map_variables variable_1="b2" variable_2="b2"/>
<map_variables variable_1="b3" variable_2="b3"/>
<map_variables variable_1="b4" variable_2="b4"/>
</connection>
</model>