- Author:
- Soroush Safaei <ssaf006@aucklanduni.ac.nz>
- Date:
- 2018-05-27 17:36:42+12:00
- Desc:
- finishing diagram
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/483/rawfile/7785e48d489152628716215553fc20a8eb21ec58/BondGraph/WholeBody.cellml
<?xml version='1.0'?>
<model name="main1" 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">
<!-- in:
out: INS,G_t,vglucoseuptake-->
<import xlink:href="Units.cellml">
<units name="mol_per_s" units_ref="mol_per_s"/>
<units name="per_s" units_ref="per_s"/>
<units name="s_per_mol" units_ref="s_per_mol"/>
</import>
<import xlink:href="Parameters-min.cellml">
<component component_ref="param1" name="param"/>
</import>
<component name="environment">
<variable name="time" public_interface="out" units="second"/>
</component>
<component name="state">
<variable name="t" public_interface="in" units="second"/>
<!-- Parameters-->
<!-- V_G: distribution volume of Glucose
k_1: rate parameter
k_2: rate parameter
k_p1: extrapolated EGP at zero Glucose and Insulin
k_p2: liver Glucose effectiveness
k_p3: amplitude of insulin action on liver
k_p4: amplitude of portal Insulin action on liver
k_i: delay between Insulin signal and Insulin action
k_abs: rate of intestinal absorption
k_gri: rate of grinding
f: fraction of intestinal absorption appears in plasma
BW: body weight
D: ingested Glucose mass
p_2U: rate constant of Insulin action on the peripheral Glucose utilisation
V_I: distribution volume of Insulin
gamma: transfer rate constant between portalvein & liver
beta: pancreatic responsivity to Glucose
alpha: delay between Glucose signal & Insulin secretion
K: pancreatic responsivity to Glucose rate of change
U_ii: Glucose insulin independent utilisation-->
<variable name="V_G" public_interface="in" units="m3"/>
<variable name="V_I" public_interface="in" units="m3"/>
<variable name="k_1" public_interface="in" units="per_s"/>
<variable name="k_2" public_interface="in" units="per_s"/>
<variable name="m_1" public_interface="in" units="per_s"/>
<variable name="m_2" public_interface="in" units="per_s"/>
<variable name="m_4" public_interface="in" units="per_s"/>
<variable name="m_5" public_interface="in" units="s_per_mol"/>
<variable name="m_6" public_interface="in" units="dimensionless"/>
<variable name="k_max" public_interface="in" units="per_s"/>
<variable name="k_min" public_interface="in" units="per_s"/>
<variable name="k_abs" public_interface="in" units="per_s"/>
<variable name="k_gri" public_interface="in" units="per_s"/>
<variable name="f" public_interface="in" units="dimensionless"/>
<variable name="b" public_interface="in" units="dimensionless"/>
<variable name="d" public_interface="in" units="dimensionless"/>
<variable name="BW" public_interface="in" units="kilogram"/>
<variable name="D" public_interface="in" units="kilogram"/>
<variable name="k_p1" public_interface="in" units="mol_per_s"/>
<variable name="k_p2" public_interface="in" units="per_s"/>
<variable name="k_p3" public_interface="in" units="per_s"/>
<variable name="k_p4" public_interface="in" units="per_s"/>
<variable name="k_i" public_interface="in" units="per_s"/>
<variable name="U_ii" public_interface="in" units="mol_per_s"/>
<variable name="V_m0" public_interface="in" units="mol_per_s"/>
<variable name="V_mX" public_interface="in" units="per_s"/>
<variable name="K_m0" public_interface="in" units="mole"/>
<variable name="k_p" public_interface="in" units="per_s"/>
<variable name="part" public_interface="in" units="dimensionless"/>
<variable name="K" public_interface="in" units="dimensionless"/>
<variable name="alpha" public_interface="in" units="per_s"/>
<variable name="beta" public_interface="in" units="dimensionless"/>
<variable name="gamma" public_interface="in" units="per_s"/>
<variable name="G_b" public_interface="in" units="mole"/>
<variable name="I_b" public_interface="in" units="mM"/>
<variable name="S_b" public_interface="in" units="mM_per_s"/>
<!-- State variables-->
<!-- G_p: Glucose mass in plasma & rapidly equilibrating tissues
G_t: Glucose mass in slowly equilibrating tissues
G: Glucose concentration in plasma
EGP: endogenous Glucose production
Ra: Glucose rate of appearance in plasma
E: renal excretion
U_id: Glucose insulin dependent utilisation
U_ii: Glucose insulin independent utilisation
EGP: endogenous Glucose production
I_1: delayed Insulin signal
I_d: delayed Insulin signal
I_po: Insulin amount in portal vein
Ra: Glucose rate of appearance in plasma
Q_sto1: Glucose mass in first stomach compartment
Q_sto2: Glucose mass in second stomach compartment
Q_gut: Glucose mass in intestine
Q_sto: Glucose mass in stomach
k_empt: rate constant of gastric emptying
U_id: Insulin dependent Glucose uptake
U_idm: Insulin dependent Glucose uptake by the muscle tissue
INS: Insulin in interstitial fluid
vglucoseuptake: Glucose uptake by the adipose tissue
I: Insulin concentration in plasma
I_p: Insulin mass in plasma
I_l: Insulin mass in liver
S: Insulin secretion
HE: Insulin hepatic extraction
Y: aux state for Insulin secretion
S: Insulin secretion to beta cells
S_po: Insulin secretion to portal vein-->
<!-- 180.15588 grams.-->
<variable initial_value="178" name="q_G_plasma" public_interface="out" units="mole"/>
<variable initial_value="135" name="q_G_tissue" public_interface="out" units="mole"/>
<variable initial_value="4.5" name="q_I_liver" public_interface="out" units="mole"/>
<variable initial_value="1.25" name="q_I_plasma" public_interface="out" units="mole"/>
<variable initial_value="3.6" name="q_I_portalVein" public_interface="out" units="mole"/>
<variable initial_value="0" name="q_I_interstitium" public_interface="out" units="mM"/>
<variable initial_value="25" name="q_I_1" public_interface="out" units="mM"/>
<variable initial_value="25" name="q_I_d" public_interface="out" units="mM"/>
<variable initial_value="78000" name="q_G_sto1" public_interface="out" units="mole"/>
<variable initial_value="0" name="q_G_sto2" public_interface="out" units="mole"/>
<variable initial_value="0" name="q_G_gut" public_interface="out" units="mole"/>
<variable name="v_G_ii" public_interface="in" units="mol_per_s"/>
<variable initial_value="0.0" name="v_G_adipose" public_interface="in" units="mol_per_s"/>
<variable name="v_G_muscle" public_interface="out" units="mol_per_s"/>
<variable name="v_G_RA" public_interface="out" units="mol_per_s"/>
<variable name="v_abs" public_interface="out" units="mol_per_s"/>
<variable name="v_empt" public_interface="out" units="mol_per_s"/>
<variable name="v_gri" public_interface="out" units="mol_per_s"/>
<variable name="v_G_EGP" public_interface="out" units="mol_per_s"/>
<variable initial_value="0.0" name="v_G_RE" public_interface="out" units="mol_per_s"/>
<variable name="v_I_secretion_portalVein" public_interface="out" units="mol_per_s"/>
<variable name="v_I_secretion" public_interface="out" units="mol_per_s"/>
<variable name="v_I_interstitium" public_interface="out" units="mol_per_s"/>
<variable name="v_I_plasma" public_interface="out" units="mol_per_s"/>
<variable name="v_I_periphery" public_interface="out" units="mol_per_s"/>
<variable name="v_I_hepaticExtraction" public_interface="out" units="mol_per_s"/>
<variable name="v_G_tissue" public_interface="out" units="mol_per_s"/>
<variable initial_value="0" name="Y" public_interface="out" units="mM_per_s"/>
<variable name="k_empt" public_interface="out" units="per_s"/>
<variable name="a" public_interface="out" units="per_mol"/>
<variable name="c" public_interface="out" units="per_mol"/>
<variable name="V_mmax" public_interface="out" units="mol_per_s"/>
<!-- Conservation laws-->
<!-- Glucose Kinetics-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_G_plasma</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_G_tissue</ci>
</apply>
<ci>v_G_EGP</ci>
<ci>v_G_RA</ci>
</apply>
<ci>v_G_RE</ci>
</apply>
<ci>v_G_ii</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_G_tissue</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_G_tissue</ci>
<ci>v_G_muscle</ci>
</apply>
<ci>v_G_adipose</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_G_tissue</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>k_1</ci>
<ci>q_G_plasma</ci>
</apply>
<apply>
<times/>
<ci>k_2</ci>
<ci>q_G_tissue</ci>
</apply>
</apply>
</apply>
<!-- Insulin Kinetics-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_I_plasma</ci>
</apply>
<apply>
<minus/>
<ci>v_I_plasma</ci>
<ci>v_I_periphery</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_I_liver</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_I_plasma</ci>
</apply>
<ci>v_I_secretion</ci>
</apply>
<ci>v_I_hepaticExtraction</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_I_plasma</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>m_1</ci>
<ci>q_I_liver</ci>
</apply>
<apply>
<times/>
<ci>m_2</ci>
<ci>q_I_plasma</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_I_periphery</ci>
<apply>
<times/>
<ci>m_4</ci>
<ci>q_I_plasma</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_I_hepaticExtraction</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<ci>m_5</ci>
</apply>
<ci>v_I_secretion</ci>
</apply>
<ci>m_6</ci>
</apply>
<ci>m_1</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<ci>m_5</ci>
</apply>
<ci>v_I_secretion</ci>
</apply>
<ci>m_6</ci>
</apply>
</apply>
</apply>
<ci>q_I_liver</ci>
</apply>
</apply>
<!-- Glucose Rate of Appearance-->
<apply>
<eq/>
<ci>v_G_RA</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>f</ci>
<ci>BW</ci>
</apply>
<ci>v_abs</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_abs</ci>
<apply>
<times/>
<ci>k_abs</ci>
<ci>q_G_gut</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_empt</ci>
<apply>
<times/>
<ci>k_empt</ci>
<ci>q_G_sto2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_gri</ci>
<apply>
<times/>
<ci>k_gri</ci>
<ci>q_G_sto1</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_G_sto1</ci>
</apply>
<apply>
<minus/>
<ci>v_gri</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_G_sto2</ci>
</apply>
<apply>
<minus/>
<ci>v_gri</ci>
<ci>v_empt</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_G_gut</ci>
</apply>
<apply>
<minus/>
<ci>v_empt</ci>
<ci>v_abs</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k_empt</ci>
<apply>
<plus/>
<ci>k_min</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<ci>k_max</ci>
<ci>k_min</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<tanh/>
<apply>
<times/>
<ci>a</ci>
<apply>
<minus/>
<apply>
<plus/>
<ci>q_G_sto1</ci>
<ci>q_G_sto2</ci>
</apply>
<apply>
<times/>
<ci>b</ci>
<ci>D</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<tanh/>
<apply>
<times/>
<ci>c</ci>
<apply>
<minus/>
<apply>
<plus/>
<ci>q_G_sto1</ci>
<ci>q_G_sto2</ci>
</apply>
<apply>
<times/>
<ci>d</ci>
<ci>D</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>a</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">2.5</cn>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>b</ci>
</apply>
<ci>D</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>c</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">2.5</cn>
<apply>
<times/>
<ci>d</ci>
<ci>D</ci>
</apply>
</apply>
</apply>
<!-- Endogenous Glucose Production-->
<apply>
<eq/>
<ci>v_G_EGP</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>k_p1</ci>
<apply>
<times/>
<ci>k_p2</ci>
<ci>q_G_plasma</ci>
</apply>
</apply>
<apply>
<times/>
<ci>k_p3</ci>
<ci>q_I_d</ci>
</apply>
</apply>
<apply>
<times/>
<ci>k_p4</ci>
<ci>q_I_portalVein</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_I_1</ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<ci>k_i</ci>
</apply>
<apply>
<minus/>
<ci>q_I_1</ci>
<apply>
<divide/>
<ci>q_I_plasma</ci>
<ci>V_I</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_I_d</ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<ci>k_i</ci>
</apply>
<apply>
<minus/>
<ci>q_I_d</ci>
<ci>q_I_1</ci>
</apply>
</apply>
</apply>
<!-- Glucose Uptake-->
<apply>
<eq/>
<ci>v_G_muscle</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>V_mmax</ci>
<ci>q_G_tissue</ci>
</apply>
<apply>
<plus/>
<ci>K_m0</ci>
<ci>q_G_tissue</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>V_mmax</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>part</ci>
</apply>
<apply>
<plus/>
<ci>V_m0</ci>
<apply>
<times/>
<ci>V_mX</ci>
<ci>q_I_interstitium</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_I_interstitium</ci>
</apply>
<ci>v_I_interstitium</ci>
</apply>
<apply>
<eq/>
<ci>v_I_interstitium</ci>
<apply>
<times/>
<ci>k_p</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<divide/>
<ci>q_I_plasma</ci>
<ci>V_I</ci>
</apply>
<ci>I_b</ci>
</apply>
<ci>q_I_interstitium</ci>
</apply>
</apply>
</apply>
<!-- Insulin Secretion-->
<apply>
<eq/>
<ci>v_I_secretion</ci>
<apply>
<times/>
<ci>gamma</ci>
<ci>q_I_portalVein</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_I_portalVein</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_I_secretion</ci>
</apply>
<ci>v_I_secretion_portalVein</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_I_secretion_portalVein</ci>
<apply>
<plus/>
<ci>Y</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>K</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_G_EGP</ci>
<ci>v_G_RA</ci>
</apply>
<ci>v_G_RE</ci>
</apply>
<ci>v_G_ii</ci>
</apply>
<ci>v_Re_1</ci>
</apply>
</apply>
<ci>V_G</ci>
</apply>
<ci>S_b</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>Y</ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<ci>alpha</ci>
</apply>
<apply>
<minus/>
<ci>Y</ci>
<apply>
<times/>
<ci>beta</ci>
<apply>
<minus/>
<apply>
<divide/>
<ci>q_G_plasma</ci>
<ci>V_G</ci>
</apply>
<ci>G_b</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<connection>
<map_components component_1="param" component_2="state"/>
<!-- Glucose Kinetics-->
<map_variables variable_1="V_G" variable_2="V_G"/>
<map_variables variable_1="k_1" variable_2="k_1"/>
<map_variables variable_1="k_2" variable_2="k_2"/>
<map_variables variable_1="U_ii" variable_2="v_G_ii"/>
<!-- Insulin Kinetics-->
<map_variables variable_1="V_I" variable_2="V_I"/>
<map_variables variable_1="m_1" variable_2="m_1"/>
<map_variables variable_1="m_2" variable_2="m_2"/>
<map_variables variable_1="m_4" variable_2="m_4"/>
<map_variables variable_1="m_5" variable_2="m_5"/>
<map_variables variable_1="m_6" variable_2="m_6"/>
<!-- Ra-->
<map_variables variable_1="k_max" variable_2="k_max"/>
<map_variables variable_1="k_min" variable_2="k_min"/>
<map_variables variable_1="k_abs" variable_2="k_abs"/>
<map_variables variable_1="k_gri" variable_2="k_gri"/>
<map_variables variable_1="f" variable_2="f"/>
<map_variables variable_1="b" variable_2="b"/>
<map_variables variable_1="d" variable_2="d"/>
<map_variables variable_1="BW" variable_2="BW"/>
<map_variables variable_1="D" variable_2="D"/>
<!-- EGP-->
<map_variables variable_1="k_p1" variable_2="k_p1"/>
<map_variables variable_1="k_p2" variable_2="k_p2"/>
<map_variables variable_1="k_p3" variable_2="k_p3"/>
<map_variables variable_1="k_p4" variable_2="k_p4"/>
<map_variables variable_1="k_i" variable_2="k_i"/>
<!-- Glucose Uptake-->
<map_variables variable_1="V_m0" variable_2="V_m0"/>
<map_variables variable_1="V_mX" variable_2="V_mX"/>
<map_variables variable_1="K_m0" variable_2="K_m0"/>
<map_variables variable_1="p_2U" variable_2="k_p"/>
<map_variables variable_1="part" variable_2="part"/>
<map_variables variable_1="I_b" variable_2="I_b"/>
<!-- Insulin Secretion-->
<map_variables variable_1="K" variable_2="K"/>
<map_variables variable_1="alpha" variable_2="alpha"/>
<map_variables variable_1="beta" variable_2="beta"/>
<map_variables variable_1="gamma" variable_2="gamma"/>
<map_variables variable_1="S_b" variable_2="S_b"/>
<map_variables variable_1="G_b" variable_2="G_b"/>
</connection>
<connection>
<map_components component_1="environment" component_2="state"/>
<map_variables variable_1="time" variable_2="t"/>
</connection>
</model>
