Location: Cardiovascular Circulation @ dc0f84d01fd9 / HeartModel / main_cv.cellml
- Author:
- soroush <ssaf006@aucklanduni.ac.nz>
- Date:
- 2021-06-13 20:55:55+12:00
- Desc:
- Added a simple heart model.
- Permanent Source URI:
- https://models.cellml.org/workspace/4ac/rawfile/dc0f84d01fd93162db1ad739f7a0d8e9e3f90a33/HeartModel/main_cv.cellml
<?xml version='1.0' encoding='UTF-8'?>
<model name="CardiovascularSystem" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:xlink="http://www.w3.org/1999/xlink">
<import xlink:href="Units.cellml">
<units name="m3" units_ref="m3"/>
<units name="m3_per_s" units_ref="m3_per_s"/>
<units name="m3_per_s2" units_ref="m3_per_s2"/>
<units name="m6_per_J" units_ref="m6_per_J"/>
<units name="per_m" units_ref="per_m"/>
<units name="J_per_m3" units_ref="J_per_m3"/>
<units name="Js_per_m6" units_ref="Js_per_m6"/>
<units name="Js2_per_m6" units_ref="Js2_per_m6"/>
<units name="Js2_per_m5" units_ref="Js2_per_m5"/>
<units name="J_per_s" units_ref="J_per_s"/>
<units name="J_per_m6" units_ref="J_per_m6"/>
<units name="Js_per_m3" units_ref="Js_per_m3"/>
<units name="UnitValve" units_ref="UnitValve"/>
</import>
<import xlink:href="Parameters_cv.cellml">
<component component_ref="Parameters_Systemic" name="Parameters_Systemic"/>
<component component_ref="Parameters_Pulmonary" name="Parameters_Pulmonary"/>
<component component_ref="Parameters_Heart" name="Parameters_Heart"/>
</import>
<component name="environment">
<variable name="time" public_interface="out" units="second"/>
</component>
<!-------------- Systemic Circulation ------------//-->
<component name="Systemic">
<variable name="t" private_interface="out" public_interface="in" units="second"/>
<!-- Constitutive Parameters-->
<variable name="C_aa" private_interface="in" units="m6_per_J"/>
<variable name="C_vc" private_interface="in" units="m6_per_J"/>
<variable name="R_aa" private_interface="in" units="Js_per_m6"/>
<variable name="R_vc" private_interface="in" units="Js_per_m6"/>
<variable name="I_aa" private_interface="in" units="Js2_per_m6"/>
<variable name="I_vc" private_interface="in" units="Js2_per_m6"/>
<!-- State Variables-->
<variable name="v_aov" private_interface="out" public_interface="in" units="m3_per_s"/>
<variable name="u_ra" public_interface="in" units="J_per_m3"/>
<variable initial_value="0.0" name="u_aa" private_interface="out" units="J_per_m3"/>
<variable initial_value="0.0" name="u_vc" private_interface="out" units="J_per_m3"/>
<variable initial_value="0.0" name="v_aa" private_interface="out" units="m3_per_s"/>
<variable initial_value="0.0" name="v_vc" private_interface="out" units="m3_per_s"/>
<!-- Constitutive Relations-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>u_aa</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>v_aov</ci>
<ci>v_aa</ci>
</apply>
<ci>C_aa</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>u_vc</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>v_aa</ci>
<ci>v_vc</ci>
</apply>
<ci>C_vc</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>v_aa</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>u_aa</ci>
<ci>u_vc</ci>
</apply>
<apply>
<times/>
<ci>v_aa</ci>
<ci>R_aa</ci>
</apply>
</apply>
<ci>I_aa</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>v_vc</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>u_vc</ci>
<ci>u_ra</ci>
</apply>
<apply>
<times/>
<ci>v_vc</ci>
<ci>R_vc</ci>
</apply>
</apply>
<ci>I_vc</ci>
</apply>
</apply>
</math>
</component>
<!-------------- Pulmonary Circulation ------------//-->
<component name="Pulmonary">
<variable name="t" public_interface="in" units="second"/>
<!-- Constitutive Parameters-->
<variable name="C_pa" private_interface="in" units="m6_per_J"/>
<variable name="C_pv" private_interface="in" units="m6_per_J"/>
<variable name="R_pa" private_interface="in" units="Js_per_m6"/>
<variable name="R_pv" private_interface="in" units="Js_per_m6"/>
<variable name="I_pa" private_interface="in" units="Js2_per_m6"/>
<variable name="I_pv" private_interface="in" units="Js2_per_m6"/>
<!-- State Variables-->
<variable initial_value="4000.0" name="u_pa" public_interface="out" units="J_per_m3"/>
<variable initial_value="0.0" name="u_pv" public_interface="out" units="J_per_m3"/>
<variable name="u_la" public_interface="in" units="J_per_m3"/>
<variable initial_value="0.0" name="v_pa" public_interface="out" units="m3_per_s"/>
<variable initial_value="0.0" name="v_pv" public_interface="out" units="m3_per_s"/>
<variable name="v_puv" public_interface="in" units="m3_per_s"/>
<!-- Constitutive Relations-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>u_pa</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>v_puv</ci>
<ci>v_pa</ci>
</apply>
<ci>C_pa</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>u_pv</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>v_pa</ci>
<ci>v_pv</ci>
</apply>
<ci>C_pv</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>v_pa</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>u_pa</ci>
<ci>u_pv</ci>
</apply>
<apply>
<times/>
<ci>v_pa</ci>
<ci>R_pa</ci>
</apply>
</apply>
<ci>I_pa</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>v_pv</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>u_pv</ci>
<ci>u_la</ci>
</apply>
<apply>
<times/>
<ci>v_pv</ci>
<ci>R_pv</ci>
</apply>
</apply>
<ci>I_pv</ci>
</apply>
</apply>
</math>
</component>
<!-------------- Heart ------------//-->
<component name="Heart">
<variable name="t" public_interface="in" units="second"/>
<!-- Constitutive Parameters-->
<variable name="mt" units="second"/>
<variable name="mta" units="second"/>
<variable name="T" private_interface="in" units="second"/>
<variable name="CQ_trv" private_interface="in" units="UnitValve"/>
<variable name="CQ_puv" private_interface="in" units="UnitValve"/>
<variable name="CQ_miv" private_interface="in" units="UnitValve"/>
<variable name="CQ_aov" private_interface="in" units="UnitValve"/>
<variable name="q_ra_0" private_interface="in" units="m3"/>
<variable name="q_rv_0" private_interface="in" units="m3"/>
<variable name="q_la_0" private_interface="in" units="m3"/>
<variable name="q_lv_0" private_interface="in" units="m3"/>
<variable name="E_lv_max" private_interface="in" units="J_per_m6"/>
<variable name="E_lv_min" private_interface="in" units="J_per_m6"/>
<variable name="E_la_max" private_interface="in" units="J_per_m6"/>
<variable name="E_la_min" private_interface="in" units="J_per_m6"/>
<variable name="E_rv_max" private_interface="in" units="J_per_m6"/>
<variable name="E_rv_min" private_interface="in" units="J_per_m6"/>
<variable name="E_ra_max" private_interface="in" units="J_per_m6"/>
<variable name="E_ra_min" private_interface="in" units="J_per_m6"/>
<variable name="tau_1_lv" private_interface="in" units="dimensionless"/>
<variable name="tau_2_lv" private_interface="in" units="dimensionless"/>
<variable name="tau_1_la" private_interface="in" units="dimensionless"/>
<variable name="tau_2_la" private_interface="in" units="dimensionless"/>
<variable name="tau_1_rv" private_interface="in" units="dimensionless"/>
<variable name="tau_2_rv" private_interface="in" units="dimensionless"/>
<variable name="tau_1_ra" private_interface="in" units="dimensionless"/>
<variable name="tau_2_ra" private_interface="in" units="dimensionless"/>
<variable name="m_1_lv" private_interface="in" units="dimensionless"/>
<variable name="m_2_lv" private_interface="in" units="dimensionless"/>
<variable name="m_1_la" private_interface="in" units="dimensionless"/>
<variable name="m_2_la" private_interface="in" units="dimensionless"/>
<variable name="m_1_rv" private_interface="in" units="dimensionless"/>
<variable name="m_2_rv" private_interface="in" units="dimensionless"/>
<variable name="m_1_ra" private_interface="in" units="dimensionless"/>
<variable name="m_2_ra" private_interface="in" units="dimensionless"/>
<variable name="g_1_lv" units="dimensionless"/>
<variable name="g_2_lv" units="dimensionless"/>
<variable name="g_1_la" units="dimensionless"/>
<variable name="g_2_la" units="dimensionless"/>
<variable name="g_1_rv" units="dimensionless"/>
<variable name="g_2_rv" units="dimensionless"/>
<variable name="g_1_ra" units="dimensionless"/>
<variable name="g_2_ra" units="dimensionless"/>
<variable name="t_onset_a" private_interface="in" units="dimensionless"/>
<variable name="E_lv" units="J_per_m6"/>
<variable name="E_la" units="J_per_m6"/>
<variable name="E_rv" units="J_per_m6"/>
<variable name="E_ra" units="J_per_m6"/>
<variable name="H_la" units="dimensionless"/>
<variable name="H_lv" units="dimensionless"/>
<variable name="H_ra" units="dimensionless"/>
<variable name="H_rv" units="dimensionless"/>
<!-- State Variables-->
<variable name="u_ra" public_interface="out" units="J_per_m3"/>
<variable name="u_rv" units="J_per_m3"/>
<variable name="u_la" public_interface="out" units="J_per_m3"/>
<variable name="u_lv" units="J_per_m3"/>
<variable name="u_aa" public_interface="in" units="J_per_m3"/>
<variable name="u_pa" public_interface="in" units="J_per_m3"/>
<variable name="v_trv" units="m3_per_s"/>
<variable name="v_puv" public_interface="out" units="m3_per_s"/>
<variable name="v_miv" units="m3_per_s"/>
<variable name="v_aov" public_interface="out" units="m3_per_s"/>
<variable name="v_vc" public_interface="in" units="m3_per_s"/>
<variable name="v_pv" public_interface="in" units="m3_per_s"/>
<variable initial_value="3000e-6" name="q_lv" units="m3"/>
<variable initial_value="3000e-6" name="q_rv" units="m3"/>
<variable initial_value="0e-6" name="q_la" units="m3"/>
<variable initial_value="0e-6" name="q_ra" units="m3"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>H_lv</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>tau_1_lv</ci>
<ci>tau_2_lv</ci>
</apply>
<apply>
<plus/>
<ci>tau_1_lv</ci>
<ci>tau_1_lv</ci>
</apply>
</apply>
<ci>m_1_lv</ci>
</apply>
<apply>
<plus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>tau_1_lv</ci>
<ci>tau_2_lv</ci>
</apply>
<apply>
<plus/>
<ci>tau_1_lv</ci>
<ci>tau_1_lv</ci>
</apply>
</apply>
<ci>m_1_lv</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>tau_1_lv</ci>
<ci>tau_2_lv</ci>
</apply>
<apply>
<plus/>
<ci>tau_2_lv</ci>
<ci>tau_2_lv</ci>
</apply>
</apply>
<ci>m_2_lv</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>H_la</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>tau_1_la</ci>
<ci>tau_2_la</ci>
</apply>
<apply>
<plus/>
<ci>tau_1_la</ci>
<ci>tau_1_la</ci>
</apply>
</apply>
<ci>m_1_la</ci>
</apply>
<apply>
<plus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>tau_1_la</ci>
<ci>tau_2_la</ci>
</apply>
<apply>
<plus/>
<ci>tau_1_la</ci>
<ci>tau_1_la</ci>
</apply>
</apply>
<ci>m_1_la</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>tau_1_la</ci>
<ci>tau_2_la</ci>
</apply>
<apply>
<plus/>
<ci>tau_2_la</ci>
<ci>tau_2_la</ci>
</apply>
</apply>
<ci>m_2_la</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>H_rv</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>tau_1_rv</ci>
<ci>tau_2_rv</ci>
</apply>
<apply>
<plus/>
<ci>tau_1_rv</ci>
<ci>tau_1_rv</ci>
</apply>
</apply>
<ci>m_1_rv</ci>
</apply>
<apply>
<plus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>tau_1_rv</ci>
<ci>tau_2_rv</ci>
</apply>
<apply>
<plus/>
<ci>tau_1_rv</ci>
<ci>tau_1_rv</ci>
</apply>
</apply>
<ci>m_1_rv</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>tau_1_rv</ci>
<ci>tau_2_rv</ci>
</apply>
<apply>
<plus/>
<ci>tau_2_rv</ci>
<ci>tau_2_rv</ci>
</apply>
</apply>
<ci>m_2_rv</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>H_ra</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>tau_1_ra</ci>
<ci>tau_2_ra</ci>
</apply>
<apply>
<plus/>
<ci>tau_1_ra</ci>
<ci>tau_1_ra</ci>
</apply>
</apply>
<ci>m_1_ra</ci>
</apply>
<apply>
<plus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>tau_1_ra</ci>
<ci>tau_2_ra</ci>
</apply>
<apply>
<plus/>
<ci>tau_1_ra</ci>
<ci>tau_1_ra</ci>
</apply>
</apply>
<ci>m_1_ra</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>tau_1_ra</ci>
<ci>tau_2_ra</ci>
</apply>
<apply>
<plus/>
<ci>tau_2_ra</ci>
<ci>tau_2_ra</ci>
</apply>
</apply>
<ci>m_2_ra</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_lv</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>E_lv_max</ci>
<ci>E_lv_min</ci>
</apply>
<ci>g_1_lv</ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci>g_1_lv</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<ci>g_2_lv</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>H_lv</ci>
</apply>
</apply>
<ci>E_lv_min</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>g_1_lv</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>mt</ci>
<apply>
<times/>
<ci>tau_1_lv</ci>
<ci>T</ci>
</apply>
</apply>
<ci>m_1_lv</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>g_2_lv</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>mt</ci>
<apply>
<times/>
<ci>tau_2_lv</ci>
<ci>T</ci>
</apply>
</apply>
<ci>m_2_lv</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>E_la</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>E_la_max</ci>
<ci>E_la_min</ci>
</apply>
<ci>g_1_la</ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci>g_1_la</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<ci>g_2_la</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>H_la</ci>
</apply>
</apply>
<ci>E_la_min</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>g_1_la</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>mta</ci>
<apply>
<times/>
<ci>tau_1_la</ci>
<ci>T</ci>
</apply>
</apply>
<ci>m_1_la</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>g_2_la</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>mta</ci>
<apply>
<times/>
<ci>tau_2_la</ci>
<ci>T</ci>
</apply>
</apply>
<ci>m_2_la</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>E_rv</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>E_rv_max</ci>
<ci>E_rv_min</ci>
</apply>
<ci>g_1_rv</ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci>g_1_rv</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<ci>g_2_rv</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>H_rv</ci>
</apply>
</apply>
<ci>E_rv_min</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>g_1_rv</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>mt</ci>
<apply>
<times/>
<ci>tau_1_rv</ci>
<ci>T</ci>
</apply>
</apply>
<ci>m_1_rv</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>g_2_rv</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>mt</ci>
<apply>
<times/>
<ci>tau_2_rv</ci>
<ci>T</ci>
</apply>
</apply>
<ci>m_2_rv</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>E_ra</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>E_ra_max</ci>
<ci>E_ra_min</ci>
</apply>
<ci>g_1_ra</ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci>g_1_ra</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<ci>g_2_ra</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>H_ra</ci>
</apply>
</apply>
<ci>E_ra_min</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>g_1_ra</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>mta</ci>
<apply>
<times/>
<ci>tau_1_ra</ci>
<ci>T</ci>
</apply>
</apply>
<ci>m_1_ra</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>g_2_ra</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>mta</ci>
<apply>
<times/>
<ci>tau_2_ra</ci>
<ci>T</ci>
</apply>
</apply>
<ci>m_2_ra</ci>
</apply>
</apply>
<!-- Elastance Functions-->
<apply>
<eq/>
<ci>mt</ci>
<apply>
<minus/>
<ci>t</ci>
<apply>
<times/>
<ci>T</ci>
<apply>
<floor/>
<apply>
<divide/>
<ci>t</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mta</ci>
<apply>
<minus/>
<apply>
<minus/>
<ci>t</ci>
<apply>
<times/>
<ci>t_onset_a</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<times/>
<ci>T</ci>
<apply>
<floor/>
<apply>
<divide/>
<apply>
<minus/>
<ci>t</ci>
<apply>
<times/>
<ci>t_onset_a</ci>
<ci>T</ci>
</apply>
</apply>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<!-- Constitutive Relations-->
<apply>
<eq/>
<ci>v_trv</ci>
<piecewise>
<piece>
<apply>
<times/>
<ci>CQ_trv</ci>
<apply>
<root/>
<apply>
<minus/>
<ci>u_ra</ci>
<ci>u_rv</ci>
</apply>
</apply>
</apply>
<apply>
<geq/>
<ci>u_ra</ci>
<ci>u_rv</ci>
</apply>
</piece>
<piece>
<cn cellml:units="m3_per_s">0.0</cn>
<apply>
<lt/>
<ci>u_ra</ci>
<ci>u_rv</ci>
</apply>
</piece>
</piecewise>
</apply>
<apply>
<eq/>
<ci>v_puv</ci>
<piecewise>
<piece>
<apply>
<times/>
<ci>CQ_puv</ci>
<apply>
<root/>
<apply>
<minus/>
<ci>u_rv</ci>
<ci>u_pa</ci>
</apply>
</apply>
</apply>
<apply>
<geq/>
<ci>u_rv</ci>
<ci>u_pa</ci>
</apply>
</piece>
<piece>
<cn cellml:units="m3_per_s">0.0</cn>
<apply>
<lt/>
<ci>u_rv</ci>
<ci>u_pa</ci>
</apply>
</piece>
</piecewise>
</apply>
<apply>
<eq/>
<ci>v_miv</ci>
<piecewise>
<piece>
<apply>
<times/>
<ci>CQ_miv</ci>
<apply>
<root/>
<apply>
<minus/>
<ci>u_la</ci>
<ci>u_lv</ci>
</apply>
</apply>
</apply>
<apply>
<geq/>
<ci>u_la</ci>
<ci>u_lv</ci>
</apply>
</piece>
<piece>
<cn cellml:units="m3_per_s">0.0</cn>
<apply>
<lt/>
<ci>u_la</ci>
<ci>u_lv</ci>
</apply>
</piece>
</piecewise>
</apply>
<apply>
<eq/>
<ci>v_aov</ci>
<piecewise>
<piece>
<apply>
<times/>
<ci>CQ_aov</ci>
<apply>
<root/>
<apply>
<minus/>
<ci>u_lv</ci>
<ci>u_aa</ci>
</apply>
</apply>
</apply>
<apply>
<geq/>
<ci>u_lv</ci>
<ci>u_aa</ci>
</apply>
</piece>
<piece>
<cn cellml:units="m3_per_s">0.0</cn>
<apply>
<lt/>
<ci>u_lv</ci>
<ci>u_aa</ci>
</apply>
</piece>
</piecewise>
</apply>
<apply>
<eq/>
<ci>u_ra</ci>
<apply>
<times/>
<ci>E_ra</ci>
<apply>
<minus/>
<ci>q_ra</ci>
<ci>q_ra_0</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_rv</ci>
<apply>
<times/>
<ci>E_rv</ci>
<apply>
<minus/>
<ci>q_rv</ci>
<ci>q_rv_0</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_la</ci>
<apply>
<times/>
<ci>E_la</ci>
<apply>
<minus/>
<ci>q_la</ci>
<ci>q_la_0</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_lv</ci>
<apply>
<times/>
<ci>E_lv</ci>
<apply>
<minus/>
<ci>q_lv</ci>
<ci>q_lv_0</ci>
</apply>
</apply>
</apply>
<!-- Conservation Laws-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_ra</ci>
</apply>
<apply>
<minus/>
<ci>v_vc</ci>
<ci>v_trv</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_rv</ci>
</apply>
<apply>
<minus/>
<ci>v_trv</ci>
<ci>v_puv</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_la</ci>
</apply>
<apply>
<minus/>
<ci>v_pv</ci>
<ci>v_miv</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_lv</ci>
</apply>
<apply>
<minus/>
<ci>v_miv</ci>
<ci>v_aov</ci>
</apply>
</apply>
</math>
</component>
<!-------------- Mapping ------------//-->
<connection>
<map_components component_1="Heart" component_2="Systemic"/>
<map_variables variable_1="v_vc" variable_2="v_vc"/>
<map_variables variable_1="u_ra" variable_2="u_ra"/>
<map_variables variable_1="u_aa" variable_2="u_aa"/>
<map_variables variable_1="v_aov" variable_2="v_aov"/>
</connection>
<connection>
<map_components component_1="Heart" component_2="Pulmonary"/>
<map_variables variable_1="u_pa" variable_2="u_pa"/>
<map_variables variable_1="v_pv" variable_2="v_pv"/>
<map_variables variable_1="u_la" variable_2="u_la"/>
<map_variables variable_1="v_puv" variable_2="v_puv"/>
</connection>
<connection>
<map_components component_1="Parameters_Systemic" component_2="Systemic"/>
<map_variables variable_1="C_aa" variable_2="C_aa"/>
<map_variables variable_1="C_vc" variable_2="C_vc"/>
<map_variables variable_1="R_aa" variable_2="R_aa"/>
<map_variables variable_1="R_vc" variable_2="R_vc"/>
<map_variables variable_1="I_aa" variable_2="I_aa"/>
<map_variables variable_1="I_vc" variable_2="I_vc"/>
</connection>
<connection>
<map_components component_1="Parameters_Pulmonary" component_2="Pulmonary"/>
<map_variables variable_1="C_pa" variable_2="C_pa"/>
<map_variables variable_1="C_pv" variable_2="C_pv"/>
<map_variables variable_1="R_pa" variable_2="R_pa"/>
<map_variables variable_1="R_pv" variable_2="R_pv"/>
<map_variables variable_1="I_pa" variable_2="I_pa"/>
<map_variables variable_1="I_pv" variable_2="I_pv"/>
</connection>
<connection>
<map_components component_1="Parameters_Heart" component_2="Heart"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="CQ_trv" variable_2="CQ_trv"/>
<map_variables variable_1="CQ_puv" variable_2="CQ_puv"/>
<map_variables variable_1="CQ_miv" variable_2="CQ_miv"/>
<map_variables variable_1="CQ_aov" variable_2="CQ_aov"/>
<map_variables variable_1="E_lv_max" variable_2="E_lv_max"/>
<map_variables variable_1="E_lv_min" variable_2="E_lv_min"/>
<map_variables variable_1="E_la_max" variable_2="E_la_max"/>
<map_variables variable_1="E_la_min" variable_2="E_la_min"/>
<map_variables variable_1="E_rv_max" variable_2="E_rv_max"/>
<map_variables variable_1="E_rv_min" variable_2="E_rv_min"/>
<map_variables variable_1="E_ra_max" variable_2="E_ra_max"/>
<map_variables variable_1="E_ra_min" variable_2="E_ra_min"/>
<map_variables variable_1="tau_1_lv" variable_2="tau_1_lv"/>
<map_variables variable_1="tau_2_lv" variable_2="tau_2_lv"/>
<map_variables variable_1="tau_1_la" variable_2="tau_1_la"/>
<map_variables variable_1="tau_2_la" variable_2="tau_2_la"/>
<map_variables variable_1="tau_1_rv" variable_2="tau_1_rv"/>
<map_variables variable_1="tau_2_rv" variable_2="tau_2_rv"/>
<map_variables variable_1="tau_1_ra" variable_2="tau_1_ra"/>
<map_variables variable_1="tau_2_ra" variable_2="tau_2_ra"/>
<map_variables variable_1="m_1_lv" variable_2="m_1_lv"/>
<map_variables variable_1="m_2_lv" variable_2="m_2_lv"/>
<map_variables variable_1="m_1_la" variable_2="m_1_la"/>
<map_variables variable_1="m_2_la" variable_2="m_2_la"/>
<map_variables variable_1="m_1_rv" variable_2="m_1_rv"/>
<map_variables variable_1="m_2_rv" variable_2="m_2_rv"/>
<map_variables variable_1="m_1_ra" variable_2="m_1_ra"/>
<map_variables variable_1="m_2_ra" variable_2="m_2_ra"/>
<map_variables variable_1="t_onset_a" variable_2="t_onset_a"/>
<map_variables variable_1="q_ra_0" variable_2="q_ra_0"/>
<map_variables variable_1="q_rv_0" variable_2="q_rv_0"/>
<map_variables variable_1="q_la_0" variable_2="q_la_0"/>
<map_variables variable_1="q_lv_0" variable_2="q_lv_0"/>
</connection>
<connection>
<map_components component_1="environment" component_2="Heart"/>
<map_variables variable_1="time" variable_2="t"/>
</connection>
<connection>
<map_components component_1="environment" component_2="Systemic"/>
<map_variables variable_1="time" variable_2="t"/>
</connection>
<connection>
<map_components component_1="environment" component_2="Pulmonary"/>
<map_variables variable_1="time" variable_2="t"/>
</connection>
<!--//////////////////////////////////////////////////////////////////////////////////////////////////////
/////////////////////////////////////////// Time Mapping /////////////////////////////////////////////
//////////////////////////////////////////////////////////////////////////////////////////////////////-->
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="Heart">
<component_ref component="Parameters_Heart"/>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="Pulmonary">
<component_ref component="Parameters_Pulmonary"/>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="Systemic">
<component_ref component="Parameters_Systemic"/>
</component_ref>
</group>
</model>