Location: Guyton, Stress, Relaxation, 2008 @ e64eb7faa3e0 / guyton_stress_relaxation_2008.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2008-11-18 04:03:27+13:00
Desc:
committing version02 of guyton_stress_relaxation_2008
Permanent Source URI:
http://models.cellml.org/workspace/guyton_stress_relaxation_2008/rawfile/e64eb7faa3e09a19b55ee6c8c38a8fa08e9b1a47/guyton_stress_relaxation_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:xlink="http://www.w3.org/1999/xlink" 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#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" cmeta:id="stress_relaxation_CellML1_0_model" name="stress_relaxation_CellML1_0_model">

<!-- ========================================	DOCUMENTATION	============================================= -->

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Guyton Model: stress_relaxation</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
        This CellML model has been validated. Due to the differences between procedural code (in this case C-code) and declarative languages (CellML), some aspects of the original model were 
		not able to be encapsulated by the CellML model (such as the damping of variables). This may affect the transient behaviour of the model, however the steady-state behaviour would remain
		the same. The equations in this file and the steady-state output from the model conform to the results from the MODSIM program.
	</para>
  </section>
  <sect1 id="sec_structure">
  <title>Model Structure</title>

<para>
Arthur Guyton (1919-2003) was an American physiologist who became famous for his 1950s experiments in which he studied the physiology of cardiac output and its relationship with the peripheral circulation.  The results of these experiments challenged the conventional wisdom that it was the heart itself that controlled cardiac output.  Instead Guyton demonstrated that it was the need of the body tissues for oxygen which was the real regulator of cardiac output.  The "Guyton Curves" describe the relationship between right atrial pressures and cardiac output, and they form a foundation for understanding the physiology of circulation. 
</para>

<para>
The Guyton model of fluid, electrolyte, and circulatory regulation is an extensive mathematical model of human circulatory physiology, capable of simulating a variety of experimental conditions, and contains a number of linked subsystems relating to circulation and its neuroendocrine control.
</para>

<para>
This is a CellML translation of the Guyton model of the regulation of the circulatory system.  The complete model consists of separate modules each of which characterise a separate physiological subsystems.  The Circulation Dynamics is the primary system, to which other modules/blocks are connected.  The other modules characterise the dynamics of the kidney, electrolytes and cell water, thirst and drinking, hormone regulation, autonomic regulation, cardiovascular system etc, and these feedback on the central circulation model.  The CellML code in these modules is based on the C code from the programme C-MODSIM created by Dr Jean-Pierre Montani.
</para>

<para>
This particular CellML model describes the effect of stress relaxation on basic venous volume (V0).  This model calculates the effect over a period of time caused by excess volume (or too little volume) in the venous tree to cause changes in the volume holding capacity of the venous tree when it is fully filled with blood but at zero pressure.  In this model, there are two separate parallel stress relaxations of the veins.  One of these has a short time constant (SRK) and the other has a long time constant (SRK2).
</para>

<informalfigure float="0" id="full_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="full_model.png"/>
  </imageobject>
</mediaobject>
<caption>A systems analysis diagram for the full Guyton model describing circulation regulation.</caption>
</informalfigure>

<informalfigure float="0" id="stress_relaxation_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="stress.png"/>
  </imageobject>
</mediaobject>
<caption>A schematic diagram of the components and processes described in the current CellML model.</caption>
</informalfigure>

<para>
There are several publications referring to the Guyton model. One of these papers is cited below:
</para>

<para>
<ulink url="http://arjournals.annualreviews.org/doi/abs/10.1146/annurev.ph.34.030172.000305">Circulation: Overall Regulation</ulink>,  A.C. Guyton, T.G. Coleman, and H.J. Granger, 1972, <ulink url="http://www.biophysj.org/">
            <emphasis>Annual Review of Physiology</emphasis>
          </ulink>, 34, 13-44.  (A <ulink url="http://arjournals.annualreviews.org/doi/pdf/10.1146/annurev.ph.34.030172.000305?cookieSet=1">PDF</ulink> version of the article are available to journal subscribers on the <emphasis>Annual Review of Physiology</emphasis> website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=4334846&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 4334846</ulink>
</para>

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

<!-- ============================================================================================================ -->

	
	


<!-- ========================================	ENVIRONMENT COMPONENT	============================================= -->
   <component name="environment">
      <variable units="second" public_interface="out" private_interface="none" name="time" cmeta:id="environment_time"/>
   </component>

<!-- ========================================	STRESS RELAXATION TOP-LEVEL COMPONENT	============================================= -->
	<component cmeta:id="stress_relaxation" name="stress_relaxation">
		

		<variable units="second" public_interface="in" private_interface="out" name="time"/>

<!-- Inputs from components in other models -->
		<variable units="dimensionless" public_interface="none" private_interface="out" name="VVE" initial_value="0.743224"/>
	</component>

<!-- INPUT CONNECTIONS -->
	<connection>
		<map_components component_2="environment" component_1="stress_relaxation"/>
		<map_variables variable_2="time" variable_1="time"/>
	</connection>

<!-- ========================================	SHORT-TERM STRESS RELAXATION	============================================= -->
	<component cmeta:id="short_term_stress_relaxation" name="short_term_stress_relaxation">
		
		
<!-- Inputs from other components -->
		<variable units="dimensionless" public_interface="in" private_interface="none" name="VVE"/>
		<variable units="second" public_interface="in" private_interface="none" name="time"/>

<!-- Outputs to other components -->
		<variable private_interface="none" name="VV7" initial_value="0.00366525" cmeta:id="short_term_stress_relaxation_VV7" units="dimensionless" public_interface="out"/>

<!-- Parameters from parameter_file -->
		<variable units="dimensionless" public_interface="in" private_interface="none" name="SR"/>
		<variable units="dimensionless" public_interface="in" private_interface="none" name="SRK"/>

		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply id="SR1_to_SR5">
				<eq/>
				<apply>
					<diff/>
					<bvar>
						<ci>time</ci>
					</bvar>
					<ci>VV7</ci>
				</apply>
				<apply>
					<divide/>
					<apply>
						<minus/>
						<apply>
							<times/>
							<apply>
								<minus/>
								<ci>VVE</ci>
								<cn cellml:units="dimensionless">0.74</cn>
							</apply>
							<ci>SR</ci>
						</apply>
						<ci>VV7</ci>
					</apply>
					<ci>SRK</ci>
				</apply>
			</apply>
		</math>
	</component>

<!-- INPUT CONNECTIONS -->
	<connection>
		<map_components component_2="stress_relaxation" component_1="short_term_stress_relaxation"/>
		<map_variables variable_2="VVE" variable_1="VVE"/>
		<map_variables variable_2="time" variable_1="time"/>
	</connection>

<!-- PARAMETER CONNECTIONS -->
	<connection>
		<map_components component_2="parameter_values" component_1="short_term_stress_relaxation"/>
		<map_variables variable_2="SR" variable_1="SR"/>
		<map_variables variable_2="SRK" variable_1="SRK"/>
	</connection>

<!-- ========================================	LONG-TERM STRESS RELAXATION	============================================= -->
	<component cmeta:id="long_term_stress_relaxation" name="long_term_stress_relaxation">
		
		
<!-- Inputs from other components -->
		<variable units="dimensionless" public_interface="in" private_interface="none" name="VVE"/>
		<variable units="second" public_interface="in" private_interface="none" name="time"/>

<!-- Outputs to other components -->
		<variable private_interface="none" name="VV6" initial_value="0.0101913" cmeta:id="long_term_stress_relaxation_VV6" units="dimensionless" public_interface="out"/>

<!-- Parameters from parameter_file -->
		<variable units="dimensionless" public_interface="in" private_interface="none" name="SR2"/>
		<variable units="dimensionless" public_interface="in" private_interface="none" name="SRK2"/>

		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply id="SR1A_to_SR5A">
				<eq/>
				<apply>
					<diff/>
					<bvar>
						<ci>time</ci>
					</bvar>
					<ci>VV6</ci>
				</apply>
				<apply>
					<divide/>
					<apply>
						<minus/>
						<apply>
							<times/>
							<apply>
								<minus/>
								<ci>VVE</ci>
								<cn cellml:units="dimensionless">0.74</cn>
							</apply>
							<ci>SR2</ci>
						</apply>
						<ci>VV6</ci>
					</apply>
					<ci>SRK2</ci>
				</apply>
			</apply>
		</math>
	</component>

<!-- INPUT CONNECTIONS -->
	<connection>
		<map_components component_2="stress_relaxation" component_1="long_term_stress_relaxation"/>
		<map_variables variable_2="VVE" variable_1="VVE"/>
		<map_variables variable_2="time" variable_1="time"/>
	</connection>

<!-- PARAMETER CONNECTIONS -->
	<connection>
		<map_components component_2="parameter_values" component_1="long_term_stress_relaxation"/>
		<map_variables variable_2="SR2" variable_1="SR2"/>
		<map_variables variable_2="SRK2" variable_1="SRK2"/>
	</connection>

<!-- ==========================================================	PARAMETER VALUES	================================================================ -->
	<component cmeta:id="parameter_values" name="parameter_values">
		<variable units="dimensionless" public_interface="out" private_interface="none" name="SR" initial_value="1"/>	<!-- sensitivity controller, short-term stress relaxation [P] -->
		<variable units="dimensionless" public_interface="out" private_interface="none" name="SR2" initial_value="1"/>	<!-- sensitivity controller, long-term stress relaxation [P] -->
		<variable units="dimensionless" public_interface="out" private_interface="none" name="SRK" initial_value="5"/>	<!-- short time constant for stress relaxation [P] -->
		<variable units="dimensionless" public_interface="out" private_interface="none" name="SRK2" initial_value="10000"/>	<!-- long time constant for stress relaxation [P] -->
	</component>
	
<!-- ============================================================		GROUPING		=============================================================== -->

	<group>
		<relationship_ref relationship="containment"/>
		<component_ref component="stress_relaxation">
			<component_ref component="short_term_stress_relaxation"/>
			<component_ref component="long_term_stress_relaxation"/>
		</component_ref>
	</group>
      
	<group>
		<relationship_ref relationship="encapsulation"/>
		<component_ref component="stress_relaxation">
			<component_ref component="parameter_values"/>
			<component_ref component="short_term_stress_relaxation"/>
			<component_ref component="long_term_stress_relaxation"/>
		</component_ref>
	</group>
	

<!-- SIMULATION METADATA -->

   	
   
<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#5c8407ba-500c-48aa-8f59-5c56c3694995">
    <rdf:li rdf:resource="rdf:#d8ba9cbc-e159-4f84-9005-45e8078a06c3"/>
    <rdf:li rdf:resource="rdf:#ad2a9a5d-5954-4ce7-8402-0a88a71266f5"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#e2d3881b-2446-4da4-8987-fa47bf6a234a">
    <dc:creator rdf:resource="rdf:#5c8407ba-500c-48aa-8f59-5c56c3694995"/>
    <dc:title/>
    <bqs:volume/>
    <bqs:first_page/>
    <bqs:Journal rdf:resource="rdf:#74d8e7e2-6e86-40d4-84f9-1c34f82fb926"/>
    <dcterms:issued rdf:resource="rdf:#715e3537-a7b4-45ad-b286-ac5da1bc458a"/>
    <bqs:last_page/>
  </rdf:Description>
  <rdf:Description rdf:about="#stress_relaxation_CellML1_0_model">
    <bqs:reference rdf:resource="rdf:#66936b4b-6377-4b0e-ab0c-d536fa1b1218"/>
    <cmeta:comment rdf:resource="rdf:#67e085bc-41db-44d6-8841-c3760fefe985"/>
    <rdf:value>
			Effect of Stress Relaxation on Basic Venous Volume (V0)
			This section calculates the effect over a period of time caused by excess volume 
			(or too little volume) in the venous tree to cause changes in the volume holding 
			capacity of the venous tree when it is fully filled with blood but at zero pressure.  
			In this model, there are two separate parallel stress relaxations of the veins.  
			One of these has a short time constant (SRK) and the other has a long time constant (SRK2).
			</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher/>
    <cmeta:comment rdf:resource="rdf:#39e3cd79-5ba8-4c2f-8282-aece51688afe"/>
    <dcterms:created rdf:resource="rdf:#a77fff93-534b-4862-a331-1e6a0287a220"/>
    <dc:creator rdf:resource="rdf:#2bbacea4-839b-4a55-94cd-a61bdb21e8c2"/>
    <cmeta:modification rdf:resource="rdf:#267d4036-bb34-4fae-bc39-bd4821da6264"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#87aac64e-8ed5-41d5-b09b-bc6e3c2d0024">
    <vCard:Given/>
    <vCard:Family>Guyton</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="#stress_relaxation">
    <rdf:value>
				Encapsulation grouping component containing all the components in the Stress Relaxation Model. The inputs and 
				outputs of the Stress Relaxation Model must be passed by this component.
				</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a77fff93-534b-4862-a331-1e6a0287a220">
    <dcterms:W3CDTF>2008-11-18T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f84fe6bf-d933-4d06-b52f-4d35c026f975">
    <vCard:Orgname>University of Auckland</vCard:Orgname>
    <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#74d8e7e2-6e86-40d4-84f9-1c34f82fb926">
    <dc:title/>
  </rdf:Description>
  <rdf:Description rdf:about="#stress_relaxation_CellML1_0_model">
    <ns7:simulation rdf:resource="rdf:#$QQv43"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d8ba9cbc-e159-4f84-9005-45e8078a06c3">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#87aac64e-8ed5-41d5-b09b-bc6e3c2d0024"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$RQv43">
    <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
    <rdf:first rdf:resource="rdf:#$SQv43"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2d425a66-033b-44f5-b3f6-7f996d9addaf">
    <vCard:N rdf:resource="rdf:#1efe3e39-192d-488b-9ff0-2d4573d9e8c7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#66936b4b-6377-4b0e-ab0c-d536fa1b1218">
    <bqs:Pubmed_id/>
    <bqs:JournalArticle rdf:resource="rdf:#e2d3881b-2446-4da4-8987-fa47bf6a234a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ad2a9a5d-5954-4ce7-8402-0a88a71266f5">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#acbe76cd-4848-45b0-b0fb-2cf2449ac9e2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$3SdBX3">
    <ns7:boundIntervals rdf:resource="rdf:#$4SdBX3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6d4b9279-d964-4a3e-9c28-a175d48e5b00">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>j.terkildsen@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#39e3cd79-5ba8-4c2f-8282-aece51688afe">
    <dc:creator rdf:resource="rdf:#b2b22de4-a14d-4aeb-a17a-68e2a8872748"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#053824a6-1920-4c95-8bbf-67f548aa4df0">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="#short_term_stress_relaxation">
    <rdf:value>
				SR1 and SR2:  
				Calculation of the ultimate degree of change in basic venous volume to 
				be caused by the short-term stress relaxation factor with input to the 
				system equal to the instantaneous excess venous volume (VVE); a multiplier 
				factor controls the degree of stress relaxation that will occur (SR).

				SR3, SR4, and SR5:  
				This is a delay circuit having a short time constant (SRK).  The output of 
				this circuit (VV7) approaches the ultimate degree of stress relaxation 
				caused by short-term stress relaxation as calculated from Block SR2.
				NB - REMOVED THE DAMPING FROM THE INTEGRAL!!!
				</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8e70348d-4fba-40f0-8204-ea015e30e073">
    <vCard:Given>Jonna</vCard:Given>
    <vCard:Family>Terkildsen</vCard:Family>
    <vCard:Other/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#715e3537-a7b4-45ad-b286-ac5da1bc458a">
    <dcterms:W3CDTF>2008-11-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2bbacea4-839b-4a55-94cd-a61bdb21e8c2">
    <vCard:ORG rdf:resource="rdf:#f84fe6bf-d933-4d06-b52f-4d35c026f975"/>
    <vCard:EMAIL rdf:resource="rdf:#6d4b9279-d964-4a3e-9c28-a175d48e5b00"/>
    <vCard:N rdf:resource="rdf:#8e70348d-4fba-40f0-8204-ea015e30e073"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$SQv43">
    <ns7:endingValue>200000</ns7:endingValue>
    <ns7:nonstandard-pointDensity>100000</ns7:nonstandard-pointDensity>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$b2Cb93">
    <ns7:boundIntervals rdf:resource="rdf:#$c2Cb93"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f0680c00-d3c7-4680-906f-4205ef68cc13">
    <dcterms:W3CDTF>2008-11-18T16:01:51+13:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$c2Cb93">
    <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
    <rdf:first rdf:resource="rdf:#$d2Cb93"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$QQv43">
    <ns7:boundIntervals rdf:resource="rdf:#$RQv43"/>
  </rdf:Description>
  <rdf:Description rdf:about="/stress_1-0.cellml#stress_relaxation_CellML1_0_model#LfY7zdZ%60EYs">
    <ns7:simulation rdf:resource="rdf:#$3SdBX3"/>
    <ns7:simulation rdf:resource="rdf:#$b2Cb93"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#67e085bc-41db-44d6-8841-c3760fefe985">
    <dc:creator rdf:resource="rdf:#053824a6-1920-4c95-8bbf-67f548aa4df0"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$4SdBX3">
    <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
    <rdf:first rdf:resource="rdf:#$5SdBX3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#acbe76cd-4848-45b0-b0fb-2cf2449ac9e2">
    <vCard:Given/>
    <vCard:Family>Stress Relaxation</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b2b22de4-a14d-4aeb-a17a-68e2a8872748">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="#long_term_stress_relaxation">
    <rdf:value>
				SR1A, SR2A, SR3A, AR4A, and SR5A:  
				Similar calculations to the above but this time with a long time constant 
				for stress relaxation (SRK2), and also having a separate variable for control 
				of the ultimate degree of the stress relaxation (SR2).  The output of this 
				long time constant stress relaxation (VV6) along with the output from the 
				short time constant stress relaxation (VV7) are subtracted from the actual 
				venous volume (VVS) in Block CD15 in the hemodynamic section of the model.
				</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#267d4036-bb34-4fae-bc39-bd4821da6264">
    <dcterms:modified rdf:resource="rdf:#f0680c00-d3c7-4680-906f-4205ef68cc13"/>
    <rdf:value>Added cmeta:id's to some variables</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#2d425a66-033b-44f5-b3f6-7f996d9addaf"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1efe3e39-192d-488b-9ff0-2d4573d9e8c7">
    <vCard:Given>Jonna</vCard:Given>
    <vCard:Family>Terkildsen</vCard:Family>
    <vCard:Other/>
  </rdf:Description>
</rdf:RDF>
</model>