Author: |
@@ -18758,27 +18707,26 @@
{120,20}}), iconTransformation(extent={{80,-20},{120,20}})));
equation
ArtyFractGz = artyFractGz[Status_Posture];
- VeinFractGz = veinFractGz[Status_Posture]; //in legs veins are valves
-
- TorsoArtyGradient = TorsoHeight * GravityAcceleration * BloodDensity * ArtyFractGz; //Hydrostatic pressure: P=ro.g.h
- //cm * 0.01 m/cm * 9.81 m/s² * 1060 kg/m3 * 1/133.322 mmHg/Pa
- //= height[cm] * 0.77996129671022036873134216408395 [mmHg]
- //Blood density = 1060 kg/m3: Cutnell, John & Johnson, Kenneth. Physics, Fourth Edition. Wiley, 1998: 308.
-
- //bloodDensity = 1060 kg/m3
- //gravityAcceleration = 9.81 m/s2
- //Pa2mmHg = 1/133.322 mmHg/Pa
- //cm2m = 0.01 m/cm
- //averageHydrostaticHeight = TorsoCM * postureCoef[Status_Posture]
- //hydrostaticPressure = averageHydrostaticHeight * cm2m * bloodDensity * gravityAcceleration * Pa2mmHg
- // = 0.77996 * averageHydrostaticHeight
-
- TorsoVeinGradient = TorsoHeight * GravityAcceleration * BloodDensity * VeinFractGz;
+ VeinFractGz = veinFractGz[Status_Posture];
+//in legs veins are valves
+ TorsoArtyGradient = TorsoHeight * GravityAcceleration * BloodDensity * ArtyFractGz;
+//Hydrostatic pressure: P=ro.g.h
+//cm * 0.01 m/cm * 9.81 m/s² * 1060 kg/m3 * 1/133.322 mmHg/Pa
+//= height[cm] * 0.77996129671022036873134216408395 [mmHg]
+//Blood density = 1060 kg/m3: Cutnell, John & Johnson, Kenneth. Physics, Fourth Edition. Wiley, 1998: 308.
+//bloodDensity = 1060 kg/m3
+//gravityAcceleration = 9.81 m/s2
+//Pa2mmHg = 1/133.322 mmHg/Pa
+//cm2m = 0.01 m/cm
+//averageHydrostaticHeight = TorsoCM * postureCoef[Status_Posture]
+//hydrostaticPressure = averageHydrostaticHeight * cm2m * bloodDensity * gravityAcceleration * Pa2mmHg
+// = 0.77996 * averageHydrostaticHeight
+ TorsoVeinGradient = TorsoHeight * GravityAcceleration * BloodDensity * VeinFractGz;
Arty = max( (fromPressure + TorsoArtyGradient), 0);
- Vein = max( toPressure + ( TorsoVeinGradient * Pump_Effect), 0); //only lower torso has LegMusclePump_Effect
- Capy = max( ( Alpha * Arty) + ( 1.0 - Alpha) * Vein, 0);
-
+ Vein = max( toPressure + ( TorsoVeinGradient * Pump_Effect), 0);
+//only lower torso has LegMusclePump_Effect
+ Capy = max(Alpha * Arty + (1.0 - Alpha) * Vein, 0);
annotation (
Documentation(info="
@@ -18800,25 +18748,25 @@
parameter Real TiltTable_Degrees( final quantity="Angle", final displayUnit="Deg")= 0;
TorsoHydrostatics UpperTorsoHydrostatics(
- artyFractGz={0,1,1,sin(Modelica.SIunits.Conversions.from_deg(
+ artyFractGz={0,1,1,sin(Modelica.Units.Conversions.from_deg(
TiltTable_Degrees)),1.76947},
- veinFractGz={0,1,1,sin(Modelica.SIunits.Conversions.from_deg(
+ veinFractGz={0,1,1,sin(Modelica.Units.Conversions.from_deg(
TiltTable_Degrees)),1.76947},
TorsoHeight=-0.1)
annotation (Placement(transformation(extent={{-12,34},{16,62}})));
TorsoHydrostatics LowerTorsoHydrostatics(
- artyFractGz={0,0.7,1,sin(Modelica.SIunits.Conversions.from_deg(
+ artyFractGz={0,0.7,1,sin(Modelica.Units.Conversions.from_deg(
TiltTable_Degrees)),0.0190301},
- veinFractGz={0.2,0.7,1,0.2 + sin(Modelica.SIunits.Conversions.from_deg(
+ veinFractGz={0.2,0.7,1,0.2 + sin(Modelica.Units.Conversions.from_deg(
TiltTable_Degrees)),0.0190301},
TorsoHeight=0.5)
annotation (Placement(transformation(extent={{-12,-42},{16,-14}})));
TorsoHydrostatics MiddleTorsoHydrostatics(
- artyFractGz={0,1,1,sin(Modelica.SIunits.Conversions.from_deg(
+ artyFractGz={0,1,1,sin(Modelica.Units.Conversions.from_deg(
TiltTable_Degrees)),-0.00024891},
- veinFractGz={0,1,1,sin(Modelica.SIunits.Conversions.from_deg(
+ veinFractGz={0,1,1,sin(Modelica.Units.Conversions.from_deg(
TiltTable_Degrees)),-0.00024891},
TorsoHeight=0.04)
annotation (Placement(transformation(extent={{-12,-4},{16,24}})));
@@ -19078,8 +19026,8 @@
Physiolibrary.Types.AmountOfSubstance ECFVActiveOsmoles;
Physiolibrary.Types.AmountOfSubstance ICFVActiveOsmoles;
Physiolibrary.Types.AmountOfSubstance ActiveOsmoles(start=11.6831341496947);
- // Real Osm_conc_CellWall;
- //Real Osm_conc_Osmoreceptors;
+ // Real Osm_conc_CellWall;
+ //Real Osm_conc_Osmoreceptors;
Physiolibrary.Types.RealIO.OsmolarityOutput
OsmBody_Osm_conc_CellWalls annotation (
Placement(transformation(extent={{60,-20},{100,20}}), iconTransformation(
@@ -19110,11 +19058,12 @@
ICFVActiveOsmoles = ICFVActiveElectrolytes + OsmCell_NonElectrolytes;
ActiveOsmoles = ECFVActiveOsmoles + ICFVActiveOsmoles;
OsmBody_Osm_conc_CellWalls = ActiveOsmoles / BodyH2O_Vol;
- Osmoreceptors = ActiveElectrolytes / BodyH2O_Vol; //, simplified=0.25331);
- ICFV = ( ICFVActiveOsmoles / ActiveOsmoles) * BodyH2O_Vol;
- ECFV = BodyH2O_Vol - ICFV; // == ( ECFVActiveOsmoles / ActiveOsmoles) * BodyH2O_Vol
-
- ECFVOsmolarity = ECFVActiveOsmoles/ECFV; // == OsmBody_Osm_conc_CellWalls
+ Osmoreceptors = ActiveElectrolytes / BodyH2O_Vol;
+//, simplified=0.25331);
+ ICFV = ICFVActiveOsmoles / ActiveOsmoles * BodyH2O_Vol;
+ ECFV = BodyH2O_Vol - ICFV;
+// == ( ECFVActiveOsmoles / ActiveOsmoles) * BodyH2O_Vol
+ ECFVOsmolarity = ECFVActiveOsmoles / ECFV; // == OsmBody_Osm_conc_CellWalls
annotation (
Documentation(info="
QHP 2008 / OsmBody
@@ -19336,8 +19285,7 @@
Placement(transformation(extent={{-120,-20},{-80,20}}),
iconTransformation(extent={{-110,-10},{-90,10}})));
- //parameter Real H2OTarget(final quantity="VolumeFlowRate",final displayUnit="ml/min") = 0;
-
+//parameter Real H2OTarget(final quantity="VolumeFlowRate",final displayUnit="ml/min") = 0;
parameter Real H2OMassEffect[ :,3] = {{ 0.0, 0.0, 0.0}, { 50.0, 1.0, 0.0}}
"volume to outtake effect";
@@ -19433,17 +19381,16 @@
Placement(transformation(extent={{-120,-60},{-80,-20}}),
iconTransformation(extent={{-120,-60},{-80,-20}})));
*/
- //parameter Physiolibrary.Types.Osmolarity Fiber_mass = 0.043;
- //parameter Real Na_EqToAllConnectedOsm(final displayUnit="mOsm/mEq") = 2;
- //parameter Real K_EqToAllConnectedOsm(final displayUnit="mOsm/mEq") = 2;
- //parameter Real initialVolume(final quantity="Volume", final displayUnit="ml") = initialValue;
+ //parameter Physiolibrary.Types.Osmolarity Fiber_mass = 0.043;
+ //parameter Real Na_EqToAllConnectedOsm(final displayUnit="mOsm/mEq") = 2;
+ //parameter Real K_EqToAllConnectedOsm(final displayUnit="mOsm/mEq") = 2;
+ //parameter Real initialVolume(final quantity="Volume", final displayUnit="ml") = initialValue;
parameter Physiolibrary.Types.Volume mass_start= 0.000949201;
Physiolibrary.Types.Volume mass(start= mass_start)
"water volume in gastro intestinal lumen";
Physiolibrary.Types.Osmolarity OsmNa;
Physiolibrary.Types.Osmolarity OsmK;
- //Physiolibrary.Types.Osmolarity Fiber;
-
+ //Physiolibrary.Types.Osmolarity Fiber;
Physiolibrary.Types.RealIO.VolumeOutput
Vol
annotation (Placement(transformation(
@@ -19453,18 +19400,16 @@
transformation(extent={{120,60},{80,100}}), iconTransformation(extent={{
-100,-50},{-80,-30}})));
initial equation
- //mass = 949.201;
- equation
- // der(mass) = (intake.q + absorbtion.q + outtake.q)/Library.SecPerMin;
-
- // env.pressure = Vol; //used in diarrea or vomitus calculation
- // OsmNa=Na_EqToAllConnectedOsm*GILumenSodium_Mass/mass;
- // OsmK = K_EqToAllConnectedOsm*GILumenPotassium_Mass/mass;
- OsmNa=2*GILumenSodium/mass;
- OsmK =2*GILumenPotassium/mass;
-
- // Fiber =Fiber_mass; ///mass;
- absorbtion.o = (Fiber + OsmNa + OsmK);
+//mass = 949.201;
+ equation
+// der(mass) = (intake.q + absorbtion.q + outtake.q)/Library.SecPerMin;
+// env.pressure = Vol; //used in diarrea or vomitus calculation
+// OsmNa=Na_EqToAllConnectedOsm*GILumenSodium_Mass/mass;
+// OsmK = K_EqToAllConnectedOsm*GILumenPotassium_Mass/mass;
+ OsmNa = 2 * GILumenSodium / mass;
+ OsmK =2*GILumenPotassium/mass;
+// Fiber =Fiber_mass; ///mass;
+ absorbtion.o = Fiber + OsmNa + OsmK;
Vol = mass;
@@ -19548,12 +19493,11 @@
annotation (Placement(transformation(extent={{80,44},{100,
64}}), iconTransformation(extent={{90,-10},{110,10}})));
parameter Physiolibrary.Types.Osmolarity Fiber = 0.043;
- //parameter Real Na_EqToAllConnectedOsm(final displayUnit="mOsm/mEq") = 2;
- //parameter Real K_EqToAllConnectedOsm(final displayUnit="mOsm/mEq") = 2;
- //parameter Real initialVolume(final quantity="Volume", final displayUnit="ml") = 1000;
-
- // Na_EqToAllConnectedOsm=Na_EqToAllConnectedOsm,
- // K_EqToAllConnectedOsm=K_EqToAllConnectedOsm,
+ //parameter Real Na_EqToAllConnectedOsm(final displayUnit="mOsm/mEq") = 2;
+ //parameter Real K_EqToAllConnectedOsm(final displayUnit="mOsm/mEq") = 2;
+ //parameter Real initialVolume(final quantity="Volume", final displayUnit="ml") = 1000;
+ // Na_EqToAllConnectedOsm=Na_EqToAllConnectedOsm,
+ // K_EqToAllConnectedOsm=K_EqToAllConnectedOsm,
Outtake vomitus
annotation (Placement(transformation(extent={{-7,-7},{7,7}},
rotation=180,
@@ -20096,8 +20040,7 @@
import Physiomodel;
extends Physiolibrary.Icons.Torso;
- // parameter Real interstitiumProteins = 2.6;
-
+// parameter Real interstitiumProteins = 2.6;
parameter Physiolibrary.Types.Volume InterstitialWater_start
"2270,5670,3400 ml";
parameter Physiolibrary.Types.Volume IntracellularWater_start
@@ -20131,8 +20074,8 @@
extent={{-6,-6},{6,6}},
origin={48,-82})));
- // Physiolibrary.Types.Volume volume;
- // Physiolibrary.Types.VolumeFlowRate change;
+// Physiolibrary.Types.Volume volume;
+ // Physiolibrary.Types.VolumeFlowRate change;
Physiolibrary.Osmotic.Sources.SolventInflux metabolicH2O(useSolutionFlowInput=true)
annotation (Placement(transformation(extent={{68,-98},{88,-78}})));
Physiolibrary.Osmotic.Sensors.FlowMeasure flowMeasure1 annotation (
@@ -20230,9 +20173,8 @@
points={{54.6,-82},{78,-82},{78,-81}},
color={0,0,127},
smooth=Smooth.None));
- // volume = extravascularH2O.WaterVolume;
- // change = extravascularH2O.q_out.q;
-
+// volume = extravascularH2O.WaterVolume;
+// change = extravascularH2O.q_out.q;
connect(vascularH2O, capyMembrane.q_in) annotation (Line(
points={{-82,-32},{-82,-14},{-60,-14}},
color={127,127,0},
@@ -20624,7 +20566,7 @@
end CerebrospinalFluid;
model Bladder
- // extends Library.PressureFlow.VolumeCompartement;
+ // extends Library.PressureFlow.VolumeCompartement;
extends Physiolibrary.Icons.Bladder;
Physiolibrary.Hydraulic.Components.ElasticVessel volumeCompartement(
stateName="BladderVolume.Mass", volume_start=0.0002)
@@ -20640,12 +20582,11 @@
transformation(extent={{-100,69},{-80,89}}), iconTransformation(
extent={{-100,-100},{-80,-80}})));
- // Real volume;
- // Real change;
+// Real volume;
+ // Real change;
equation
- // volume = volumeCompartement.Volume;
- // change = volumeCompartement.con.q;
-
+// volume = volumeCompartement.Volume;
+// change = volumeCompartement.con.q;
connect(busConnector.BladderVoidFlow, bladderVoidFlow.solutionFlow)
annotation (Line(
points={{-90,79},{-90,0},{-65.5,0},{-65.5,-10.25}},
@@ -20751,8 +20692,7 @@
equation
q_out.q + CD_H2O_Reab.q = 0;
q_out.o = CD_H2O_Reab.o;
-
- // der(Inactive) = ((1/InactivateKinv) * CD_H2O_Reab.q - ReactivateK * Inactive) / Library.SecPerMin;
+// der(Inactive) = ((1/InactivateKinv) * CD_H2O_Reab.q - ReactivateK * Inactive) / Library.SecPerMin;
Active = 2 - Inactive;
state = Inactive;
@@ -20944,9 +20884,7 @@
points={{48,-8},{48,2}},
color={0,0,127},
smooth=Smooth.None));
-
- //!!!
-
+//!!!
connect(busConnector.DT_Na_Outflow, ADHEffect.yBase) annotation (Line(
points={{-92,86},{24,86},{24,62}},
color={0,0,127},
@@ -21268,7 +21206,7 @@
end Peritoneum_const;
model Bladder_steady
- // extends Library.PressureFlow.VolumeCompartement;
+ // extends Library.PressureFlow.VolumeCompartement;
extends Physiolibrary.Icons.Bladder;
Physiolibrary.Hydraulic.Interfaces.HydraulicPort_a con annotation (
Placement(transformation(extent={{-84,18},{-44,58}}),
@@ -21343,7 +21281,7 @@
end Bladder_steady;
model Bladder_steady2
- // extends Library.PressureFlow.VolumeCompartement;
+ // extends Library.PressureFlow.VolumeCompartement;
extends Physiolibrary.Icons.Bladder;
Physiolibrary.Osmotic.Interfaces.OsmoticPort_a con annotation (
Placement(transformation(extent={{-84,18},{-44,58}}),
@@ -22129,9 +22067,9 @@
string="%second",
index=1,
extent={{5,2},{5,2}}));
- //ifv=bone.FractIFV +brain.FractIFV +fat.FractIFV +GITract.FractIFV +kidney.FractIFV +leftHeart.FractIFV +liver.FractIFV +otherTissue.FractIFV +respiratoryMuscle.FractIFV +rightHeart.FractIFV +skin.FractIFV +skeletalMuscle.FractIFV;
- //organ=bone.FractOrganH2O +brain.FractOrganH2O +fat.FractOrganH2O +GITract.FractOrganH2O +kidney.FractOrganH2O +leftHeart.FractOrganH2O +liver.FractOrganH2O +otherTissue.FractOrganH2O +respiratoryMuscle.FractOrganH2O +rightHeart.FractOrganH2O +skin.FractOrganH2O +skeletalMuscle.FractOrganH2O;
- /*
+//ifv=bone.FractIFV +brain.FractIFV +fat.FractIFV +GITract.FractIFV +kidney.FractIFV +leftHeart.FractIFV +liver.FractIFV +otherTissue.FractIFV +respiratoryMuscle.FractIFV +rightHeart.FractIFV +skin.FractIFV +skeletalMuscle.FractIFV;
+//organ=bone.FractOrganH2O +brain.FractOrganH2O +fat.FractOrganH2O +GITract.FractOrganH2O +kidney.FractOrganH2O +leftHeart.FractOrganH2O +liver.FractOrganH2O +otherTissue.FractOrganH2O +respiratoryMuscle.FractOrganH2O +rightHeart.FractOrganH2O +skin.FractOrganH2O +skeletalMuscle.FractOrganH2O;
+/*
assert(bone.FractIFV +
brain.FractIFV +
fat.FractIFV +
@@ -23572,8 +23510,7 @@
model Water
import QHP = Physiomodel;
extends Physiolibrary.Icons.Water;
- // extends Physiomodel.Library.Utilities.DynamicState(stateName="BodyH2O.Vol", initType=Library.Utilities.Init.NoInit, STEADY=false);
-
+ // extends Physiomodel.Library.Utilities.DynamicState(stateName="BodyH2O.Vol", initType=Library.Utilities.Init.NoInit, STEADY=false);
Physiolibrary.Osmotic.Sources.SolventOutflux Hemorrhage(
useSolutionFlowInput=true)
annotation (Placement(transformation(extent={{16,47},{32,63}})));
@@ -23599,7 +23536,7 @@
SweatFract=0.33,
SkinFract=0.33)
annotation (Placement(transformation(extent={{21,-31},{41,-11}})));
- // interstitiumProteins=3.3,
+ // interstitiumProteins=3.3,
QHP.Water.WaterCompartments.GI_Absorption GILumen
annotation (Placement(transformation(extent={{-56,54},{-36,74}})));
QHP.Water.WaterCompartments.Kidney.Kidney Kidney
@@ -23609,15 +23546,14 @@
QHP.Water.WaterProperties waterProperties
annotation (Placement(transformation(extent={{-80,-89},{-60,-69}})));
- // Real bodyH2O;
- /*initial equation
- if STEADY then
- waterProperties.BodyH2O_Change.y = 0;
- end if;
-*/
-
- //initial equation
- // bodyH2O = waterProperties.BodyH2O.y;
+// Real bodyH2O;
+ /*initial equation
+ if STEADY then
+ waterProperties.BodyH2O_Change.y = 0;
+ end if;
+ */
+ //initial equation
+ // bodyH2O = waterProperties.BodyH2O.y;
Physiolibrary.Osmotic.Components.OsmoticCell plasma(
useImpermeableSolutesInput=true,
NumberOfMembraneTypes=2,
@@ -23652,7 +23588,7 @@
SweatFract=0.33,
SkinFract=0.33)
annotation (Placement(transformation(extent={{21,-90},{41,-70}})));
- // interstitiumProteins=2.2,
+ // interstitiumProteins=2.2,
QHP.Water.WaterCompartments.Peritoneum_const peritoneum_const(
initialVolume=1e-09)
annotation (Placement(transformation(extent={{73,-16},{93,4}})));
@@ -23669,9 +23605,8 @@
QHP.Water.WaterCompartments.LungEdema_const lungEdema_const
annotation (Placement(transformation(extent={{72,14},{92,34}})));
equation
- // changePerMin = waterProperties.BodyH2O_Change.y;
- // stateValue = bodyH2O;
-
+// changePerMin = waterProperties.BodyH2O_Change.y;
+// stateValue = bodyH2O;
connect(IVDrip.solutionFlow, busConnector.IVDrip_H2ORate) annotation (Line(
points={{24,34.6},{61,34.6},{61,90},{-89,90}},
color={0,0,127},
@@ -24798,7 +24733,7 @@
extends Physiolibrary.Icons.ProteinsLib;
model Synthesis
- // parameter Physiolibrary.Types.MassFlowRate SynthesisBasic "10 mg/min";
+ // parameter Physiolibrary.Types.MassFlowRate SynthesisBasic "10 mg/min";
parameter Physiolibrary.Types.MolarFlowRate SynthesisBasic = 2.75753e-09
"10 mg/min";
parameter Real[:,3] data = {{ 20.0, 3.0, 0.0}, { 28.0, 1.0, -0.2}, { 40.0, 0.0, 0.0}}
@@ -24813,20 +24748,19 @@
-10,-110; 10,-90], Placement(transformation(extent={{90,-10},{110,10}})));
Physiolibrary.Types.Pressure COP;
- // Physiolibrary.Types.AmountOfSubstance synthetizedAmount(start=0);
- // Physiolibrary.Types.Mass synthetizedMass(start=0);
- //protected
- // constant Physiolibrary.Types.Time sec=1;
- // constant Physiolibrary.Types.Volume ghostPlasmaVol=3.02e-3
- // "Strange dependence derived from original HumMod";
+ // Physiolibrary.Types.AmountOfSubstance synthetizedAmount(start=0);
+ // Physiolibrary.Types.Mass synthetizedMass(start=0);
+ //protected
+ // constant Physiolibrary.Types.Time sec=1;
+ // constant Physiolibrary.Types.Volume ghostPlasmaVol=3.02e-3
+ // "Strange dependence derived from original HumMod";
equation
COP = q_out.conc * Modelica.Constants.R * 310.15;
c.u=COP;
q_out.q = -SynthesisBasic * c.val;
-
- //TODO: state
- //der(synthetizedAmount) = -q_out.q;
- // ProteinsMass2AmountOfSubstance(synthetizedMass,ghostPlasmaVol) = synthetizedAmount;
+//TODO: state
+//der(synthetizedAmount) = -q_out.q;
+// ProteinsMass2AmountOfSubstance(synthetizedMass,ghostPlasmaVol) = synthetizedAmount;
annotation (
defaultComponentName="synthesis",
Icon(coordinateSystem(preserveAspectRatio=true, extent={{-100,-100},{
@@ -24844,9 +24778,9 @@
end Synthesis;
model Degradation
- // parameter Physiolibrary.Types.MassFlowRate DegradationBasic "10 mg/min";
- // parameter Real[:,3] data = {{ 0.00, 0.0, 0.0}, { 0.07, 1.0, 40.0}, { 0.09, 6.0, 0.0}}
- // "ProteinEffect";
+ // parameter Physiolibrary.Types.MassFlowRate DegradationBasic "10 mg/min";
+ // parameter Real[:,3] data = {{ 0.00, 0.0, 0.0}, { 0.07, 1.0, 40.0}, { 0.09, 6.0, 0.0}}
+ // "ProteinEffect";
parameter Physiolibrary.Types.MolarFlowRate DegradationBasic = 2.75753e-09
"10 mg/min";
parameter Real[:,3] data = {{ 0.00, 0.0, 0.0}, { 1.45, 1.0, 1.59}, { 1.97, 6.0, 0.0}}
@@ -24860,14 +24794,14 @@
transformation(extent={{-100,0},{-60,40}}), iconTransformation(extent=
{{-110,-10},{-90,10}})));
- // Physiolibrary.Types.AmountOfSubstance degradedAmount(start=0);
- // Physiolibrary.Types.Mass degradedMass(start=0);
- //protected
- // constant Physiolibrary.Types.Time sec=1;
- // constant Physiolibrary.Types.Volume ghostPlasmaVol=3.02e-3
- // "Strange dependence derived from original HumMod";
+// Physiolibrary.Types.AmountOfSubstance degradedAmount(start=0);
+ // Physiolibrary.Types.Mass degradedMass(start=0);
+ //protected
+ // constant Physiolibrary.Types.Time sec=1;
+ // constant Physiolibrary.Types.Volume ghostPlasmaVol=3.02e-3
+ // "Strange dependence derived from original HumMod";
equation
- // ProteinsMassConcentration2Concentration(c.u*1000) = q_in.conc;
+// ProteinsMassConcentration2Concentration(c.u*1000) = q_in.conc;
c.u = q_in.conc;
q_in.q = DegradationBasic * c.val;
// q_in.q =ProteinsMass2AmountOfSubstance(DegradationBasic*c.val*sec,ghostPlasmaVol)/sec;
@@ -24935,10 +24869,9 @@
albumin = totalProteins * (0.63/1.45);
nonAlbuminMolarConc = totalProteins - albumin;
nonAlbumin = nonAlbuminMM * nonAlbuminMolarConc;
-
- // inversion of totalProteins=(320*101325/760)/(310.15*8.314) *0.001*tProtMassConc + (1160*101325/760)/(310.15*8.314)* (0.001*tProtMassConc)^2;
- // tProtMassConc = 0.0000170159 * (-8.106e6 + 63.6632 * ((1.6212e10 + 1.4208e10 * totalProteins)^0.5));
- // linear aproximation at point totalProteins = 1.45 mmol/l :
+// inversion of totalProteins=(320*101325/760)/(310.15*8.314) *0.001*tProtMassConc + (1160*101325/760)/(310.15*8.314)* (0.001*tProtMassConc)^2;
+// tProtMassConc = 0.0000170159 * (-8.106e6 + 63.6632 * ((1.6212e10 + 1.4208e10 * totalProteins)^0.5));
+// linear aproximation at point totalProteins = 1.45 mmol/l :
nonAlbuminMM = 34.16-10*(totalProteins-1.45);
end ProteinDivision;
@@ -25616,9 +25549,9 @@
package Metabolism "Body Nutrients, Metabolism and Heat Systems"
extends Physiolibrary.Icons.NutrientsMetabolismLib;
model LiverMetabolism
- // import QHP = Physiomodel;
+ // import QHP = Physiomodel;
extends Physiolibrary.Icons.Liver;
- //outer parameter Real EFC_Vol(final displayUnit="ml");
+ //outer parameter Real EFC_Vol(final displayUnit="ml");
Physiolibrary.Chemical.Components.Substance Glycogen(solute_start=120e-6/
Substances.Glucose.mw)
"amount of bound glucose molecules in glycogen (initial mass 120 mg of bound glucose molecules)"
@@ -26310,7 +26243,7 @@
model Glucose
extends Physiolibrary.Icons.Glucose;
- //outer parameter Real EFC_Vol(final displayUnit="ml");
+ //outer parameter Real EFC_Vol(final displayUnit="ml");
Physiolibrary.Chemical.Components.Substance GlucosePool(
stateName="GlucosePool.Mass",
useNormalizedVolume=false,
@@ -26955,7 +26888,7 @@
model Lipids
extends Physiolibrary.Icons.Lipids;
- //outer parameter Real EFC_Vol(final displayUnit="ml");
+ //outer parameter Real EFC_Vol(final displayUnit="ml");
Physiolibrary.Chemical.Components.Substance LipidDepositsFA(stateName=
"LipidDeposits.Mass", solute_start=12.0324/Substances.FattyAcids.mw)
"amount of fatty acids stored as lipid deposits (initial 12.0324 kg)"
@@ -27307,7 +27240,7 @@
model Lactate
extends Physiolibrary.Icons.Lactate;
- //outer parameter Real EFC_Vol(final displayUnit="ml");
+ //outer parameter Real EFC_Vol(final displayUnit="ml");
Physiolibrary.Chemical.Interfaces.ChemicalPort_b Lactate(q(final
displayUnit="mEq/min"), conc(
final displayUnit="mEq/ml",
@@ -27396,7 +27329,7 @@
model AminoAcids
extends Physiolibrary.Icons.AminoAcids;
- //outer parameter Real EFC_Vol(final displayUnit="ml");
+ //outer parameter Real EFC_Vol(final displayUnit="ml");
Physiolibrary.Chemical.Interfaces.ChemicalPort_b AminoAcids(conc(final
displayUnit="mg/ml"), q(final displayUnit="mg/min"))
"extracellular storage" annotation (Placement(transformation(extent={{80,34},
@@ -27523,15 +27456,15 @@
parameter Boolean canBurnFattyAcids = true;
parameter Physiolibrary.Types.MolarEnergy glucoseEnergy = (1/0.2439e-6)*4.1864*Substances.Glucose.mw
- "Metabolism-Tools.CarboAerobic(mg/cal)=0.2439"; // 1/(kg/cal)*(J/cal)*(kg/mol) = J/mol
+ "Metabolism-Tools.CarboAerobic(mg/cal)=0.2439"; // 1/(kg/cal)*(J/cal)*(kg/mol) = J/mol
parameter Physiolibrary.Types.MolarEnergy fattyAcidsEnergy = (1/0.1075e-6)*4.1864*Substances.FattyAcids.mw
"Metabolism-Tools.Fat(mg/cal)=0.1075";
parameter Physiolibrary.Types.MolarEnergy lactateEnergy = (1/0.2538e-6)*4.1864*Substances.Lactate.mw
"Metabolism-Tools.Lac-(mg/cal)=0.2438";
parameter Physiolibrary.Types.MolarEnergy ketoAcidsEnergy = (1/0.1075e-6)*4.1864*Substances.KetoAcids.mw
"Metabolism-Tools.Fat(mg/cal)=0.1075";
- parameter Physiolibrary.Types.MolarEnergy oxygenEnergy = (1/0.2093)*4.1864*22710.95322615; // 22710.95322615 ml/mol of oxygen at 0 degC
-
+ parameter Physiolibrary.Types.MolarEnergy oxygenEnergy = (1/0.2093)*4.1864*22710.95322615;
+// 22710.95322615 ml/mol of oxygen at 0 degC
parameter Physiolibrary.Types.MolarEnergy glucoseAnaerobicEnergy = glucoseEnergy - lactateEnergy;
parameter Physiolibrary.Types.MolarEnergy glycogenAnaerobicEnergy = glucoseAnaerobicEnergy;
/*
@@ -27576,8 +27509,8 @@
parameter Real[:,3] LacFractionData
"fraction of oxygen to be use for lactate metabolism (depends on lactate concentration)";
- // parameter Real[:,3] PO2OnAerobicFractionData
- // "Aerobic Fraction of O2 tissue use depents on pO2";
+// parameter Real[:,3] PO2OnAerobicFractionData
+ // "Aerobic Fraction of O2 tissue use depents on pO2";
parameter Real[:,3] PO2OnAerobicFractionData = {{2,0,0},{10,1,0}}
"Aerobic Fraction of O2 tissue use depents on pO2";
@@ -27590,12 +27523,12 @@
origin={-10,52})));
Physiolibrary.Chemical.Interfaces.ChemicalPort_a lactate(q(nominal=1e-6))
annotation (Placement(transformation(extent={{46,47},{66,67}}),
- iconTransformation(extent={{-76,90},{-56,110}}))); //(q(final displayUnit="mEq/min"), conc(final displayUnit="mEq/ml")) "in mEq/ml"
-
+ iconTransformation(extent={{-76,90},{-56,110}})));
+//(q(final displayUnit="mEq/min"), conc(final displayUnit="mEq/ml")) "in mEq/ml"
Physiolibrary.Chemical.Interfaces.ChemicalPort_a glucose(q(nominal=1e-6))
annotation (Placement(transformation(extent={{45,-4},{65,16}}),
- iconTransformation(extent={{-110,70},{-90,90}}))); //(q(final displayUnit="mg/min"), conc(final displayUnit="mg/ml")) "in mg/ml"
-
+ iconTransformation(extent={{-110,70},{-90,90}})));
+//(q(final displayUnit="mg/min"), conc(final displayUnit="mg/ml")) "in mg/ml"
Physiolibrary.Types.RealIO.MolarFlowRateOutput
Tissue_CO2FromMetabolism
annotation (
@@ -27632,15 +27565,15 @@
iconTransformation(
extent={{-11.4284,-11.4285},{28.5715,28.5715}},
origin={-104.572,-48.5715})));
- // Real O2Use;
+ // Real O2Use;
Physiolibrary.Types.MolarFlowRate O2UseByGlu,O2UseByLac,O2UseByFA, O2UseByKA;
Physiolibrary.Types.Fraction FAfraction, KAfraction, Ratio;
Physiolibrary.Types.Power AnaerobicCals;
Physiolibrary.Chemical.Interfaces.ChemicalPort_a fattyAcids
annotation (Placement(transformation(extent={{45,-60},{65,-40}}),
- iconTransformation(extent={{-110,-10},{-90,10}}))); //(q(final displayUnit="mg/min"), conc(final displayUnit="mg/ml")) "in mg/ml"
-
+ iconTransformation(extent={{-110,-10},{-90,10}})));
+//(q(final displayUnit="mg/min"), conc(final displayUnit="mg/ml")) "in mg/ml"
Physiolibrary.Chemical.Interfaces.ChemicalPort_a glycogen "glycogen flow"
annotation (Placement(transformation(extent={{4,-50},{26,-28}}),
iconTransformation(extent={{-110,30},{-90,50}})));
@@ -27662,72 +27595,61 @@
annotation (Placement(transformation(extent={{67,-40},{87,-20}}),
iconTransformation(extent={{-38,90},{-18,110}}))); //(q(final displayUnit="mg/min"), conc(final displayUnit="mg/ml")) "in mg/ml"
- Modelica.SIunits.Power gluW=glucose.q*glucoseEnergy, gluO2W=O2UseByGlu*oxygenEnergy, lacW=lactate.q*lactateEnergy, lacO2W=LacFraction.val * CalsUse, anaerobicGluW=glucoseAnaerobicEnergy *(glucose.q - (O2UseByGlu*oxygenEnergy)/glucoseEnergy), anaerobicGlyW=glycogenAnaerobicEnergy*glycogen.q;
- Modelica.SIunits.MolarFlowRate anaerobicGluGly=((glucose.q - (O2UseByGlu*oxygenEnergy)/glucoseEnergy) + glycogen.q), anaerobicProducedLac=(lactate.q - LacFraction.val * CalsUse/lactateEnergy);
+ Modelica.Units.SI.Power gluW=glucose.q*glucoseEnergy, gluO2W=O2UseByGlu*oxygenEnergy, lacW=lactate.q*lactateEnergy, lacO2W=LacFraction.val * CalsUse, anaerobicGluW=glucoseAnaerobicEnergy *(glucose.q - (O2UseByGlu*oxygenEnergy)/glucoseEnergy), anaerobicGlyW=glycogenAnaerobicEnergy*glycogen.q;
+ Modelica.Units.SI.MolarFlowRate anaerobicGluGly=((glucose.q - (O2UseByGlu*oxygenEnergy)/glucoseEnergy) + glycogen.q), anaerobicProducedLac=(lactate.q - LacFraction.val * CalsUse/lactateEnergy);
equation
AerobicFraction.u = pO2;
- LacFraction.u = lactate.conc; // * lacDensity * 100.0; // [mg/dl]
- O2Use*oxygenEnergy = AerobicFraction.val * CalsUse; //O2Use = AerobicFraction.val * eTOo2_coef * CalsUse;
- O2UseByLac = LacFraction.val * O2Use;
-
- Tissue_MetabolicH2ORate = (O2Use/2)*Substances.Water.mw / WaterDensity;
-
- if canBurnFattyAcids then
- ketoAcids.q=0;
- KAfraction=0;
- O2UseByKA=0;
-
- O2UseByFA = FAfraction * O2Use;
- O2UseByGlu + O2UseByLac + O2UseByFA = O2Use;
- FAfraction = (1-LacFraction.val)*(Ratio/(Ratio+0.026*(Substances.Glucose.mw/Substances.FattyAcids.mw)));
- Ratio = (fattyAcids.conc) / (glucose.conc); // Ratio = fattyAcids.conc / glucose.conc;
-
- //Tissue_CO2FromMetabolism = 0.0446 * (O2UseByGlu*RQglu + O2UseByLac*RQlac + O2UseByFA*RQfat);
- Tissue_CO2FromMetabolism = (O2UseByGlu*RQglu + O2UseByLac*RQlac + O2UseByFA*RQfat);
-
- // Tissue_MetabolicH2ORate = 0.000176 * CalsUse;
- // Tissue_MetabolicH2ORate = (0.176e-9/4.1864) * CalsUse; // 0.000176 ml/cal = 0.176e-9/4.1864 m3/J
-
- fattyAcids.q*fattyAcidsEnergy = FAfraction * CalsUse;
- else
- fattyAcids.q=0;
- FAfraction=0;
- O2UseByFA=0;
-
- O2UseByKA = KAfraction * O2Use;
- O2UseByGlu + O2UseByLac + O2UseByKA = O2Use;
- KAfraction = (1-LacFraction.val)*(Ratio/(Ratio+0.222*(Substances.Glucose.mw/Substances.KetoAcids.mw)));
- Ratio = (ketoAcids.conc) / (glucose.conc);
-
- Tissue_CO2FromMetabolism = (O2UseByGlu*RQglu + O2UseByLac*RQlac + O2UseByKA*RQketo);
- //Tissue_MetabolicH2ORate = (0.176e-9/4186.4) * CalsUse; // 0.000176 ml/cal = 0.176e-9/4186.4 m3/J
-
- ketoAcids.q*ketoAcidsEnergy = KAfraction * CalsUse;
- end if;
-
- GlycogenAvailabilityCurve.u = glycogen.conc;
-
- if AerobicFraction.val<1-Modelica.Constants.eps then
- // with anaerobic metabolism
- AnaerobicCals = (1-AerobicFraction.val)*CalsUse; //anaerobic callories
-
- glycogen.q = GlycogenAvailabilityCurve.val * AnaerobicCals;
-
- // anaerobic_glu2energy *(glucose.q - (O2UseByGlu/eTOo2_coef) * eTOglu_coef) + gly2energy*glycogen.q = AnaerobicCals;
- // (glucose.q - (O2UseByGlu/eTOo2_coef) * eTOglu_coef) + gly2lac*glycogen.q = (lactate.q*lacDensity - LacFraction.val * CalsUse * eTOlac_coef);
-
- glucoseAnaerobicEnergy *(glucose.q - (O2UseByGlu*oxygenEnergy)/glucoseEnergy) + glycogenAnaerobicEnergy*glycogen.q = AnaerobicCals;
-
- ((glucose.q - (O2UseByGlu*oxygenEnergy)/glucoseEnergy) + glycogen.q) = -2*(lactate.q - LacFraction.val * CalsUse/lactateEnergy);
-
- else
- // only aerobic metabolism
- AnaerobicCals = 0.0;
- glycogen.q = 0.0;
- lactate.q*lactateEnergy = LacFraction.val * CalsUse;
- glucose.q*glucoseEnergy = O2UseByGlu*oxygenEnergy;
+ LacFraction.u = lactate.conc;
+// * lacDensity * 100.0; // [mg/dl]
+ O2Use * oxygenEnergy = AerobicFraction.val * CalsUse;
+//O2Use = AerobicFraction.val * eTOo2_coef * CalsUse;
+ O2UseByLac = LacFraction.val * O2Use;
+
+ Tissue_MetabolicH2ORate = (O2Use/2)*Substances.Water.mw / WaterDensity;
+ if canBurnFattyAcids then
+ ketoAcids.q = 0;
+ KAfraction = 0;
+ O2UseByKA = 0;
+ O2UseByFA = FAfraction * O2Use;
+ O2UseByGlu + O2UseByLac + O2UseByFA = O2Use;
+ FAfraction = (1 - LacFraction.val) * (Ratio / (Ratio + 0.026 * (Substances.Glucose.mw / Substances.FattyAcids.mw)));
+ Ratio = fattyAcids.conc / glucose.conc;
+// Ratio = fattyAcids.conc / glucose.conc;
+//Tissue_CO2FromMetabolism = 0.0446 * (O2UseByGlu*RQglu + O2UseByLac*RQlac + O2UseByFA*RQfat);
+ Tissue_CO2FromMetabolism = O2UseByGlu * RQglu + O2UseByLac * RQlac + O2UseByFA * RQfat;
+// Tissue_MetabolicH2ORate = 0.000176 * CalsUse;
+// Tissue_MetabolicH2ORate = (0.176e-9/4.1864) * CalsUse; // 0.000176 ml/cal = 0.176e-9/4.1864 m3/J
+ fattyAcids.q * fattyAcidsEnergy = FAfraction * CalsUse;
+ else
+ fattyAcids.q = 0;
+ FAfraction = 0;
+ O2UseByFA = 0;
+ O2UseByKA = KAfraction * O2Use;
+ O2UseByGlu + O2UseByLac + O2UseByKA = O2Use;
+ KAfraction = (1 - LacFraction.val) * (Ratio / (Ratio + 0.222 * (Substances.Glucose.mw / Substances.KetoAcids.mw)));
+ Ratio = ketoAcids.conc / glucose.conc;
+ Tissue_CO2FromMetabolism = O2UseByGlu * RQglu + O2UseByLac * RQlac + O2UseByKA * RQketo;
+//Tissue_MetabolicH2ORate = (0.176e-9/4186.4) * CalsUse; // 0.000176 ml/cal = 0.176e-9/4186.4 m3/J
+ ketoAcids.q * ketoAcidsEnergy = KAfraction * CalsUse;
end if;
+ GlycogenAvailabilityCurve.u = glycogen.conc;
+ if AerobicFraction.val < 1 - Modelica.Constants.eps then
+// with anaerobic metabolism
+ AnaerobicCals = (1 - AerobicFraction.val) * CalsUse;
+//anaerobic callories
+ glycogen.q = GlycogenAvailabilityCurve.val * AnaerobicCals;
+// anaerobic_glu2energy *(glucose.q - (O2UseByGlu/eTOo2_coef) * eTOglu_coef) + gly2energy*glycogen.q = AnaerobicCals;
+// (glucose.q - (O2UseByGlu/eTOo2_coef) * eTOglu_coef) + gly2lac*glycogen.q = (lactate.q*lacDensity - LacFraction.val * CalsUse * eTOlac_coef);
+ glucoseAnaerobicEnergy * (glucose.q - O2UseByGlu * oxygenEnergy / glucoseEnergy) + glycogenAnaerobicEnergy * glycogen.q = AnaerobicCals;
+ glucose.q - O2UseByGlu * oxygenEnergy / glucoseEnergy + glycogen.q = -2 * (lactate.q - LacFraction.val * CalsUse / lactateEnergy);
+ else
+// only aerobic metabolism
+ AnaerobicCals = 0.0;
+ glycogen.q = 0.0;
+ lactate.q * lactateEnergy = LacFraction.val * CalsUse;
+ glucose.q * glucoseEnergy = O2UseByGlu * oxygenEnergy;
+ end if;
annotation ( Documentation(revisions="
@@ -28068,19 +27990,18 @@
partial model TissueMetabolismBase
extends Physiolibrary.Icons.MetabolismPart;
- parameter Physiolibrary.Types.AmountOfSubstance initialTissueO2(displayUnit="ml_STP")=1.7612646902879852858386051702718e-4; //pO2*solubility*LiquidVolume
- // parameter Physiolibrary.Types.GasSolubility O2solubility(final displayUnit="ml/mmHg"); //=0.00003;
- // parameter Real O2fromBloodtoTissueConductance(final displayUnit="ml/(min.mmHg)")=1000000;
-
+ parameter Physiolibrary.Types.AmountOfSubstance initialTissueO2(displayUnit="ml_STP")=1.7612646902879852858386051702718e-4; //pO2*solubility*LiquidVolume
+ // parameter Physiolibrary.Types.GasSolubility O2solubility(final displayUnit="ml/mmHg"); //=0.00003;
+ // parameter Real O2fromBloodtoTissueConductance(final displayUnit="ml/(min.mmHg)")=1000000;
parameter Physiolibrary.Types.AmountOfSubstance initialLactateMass(displayUnit="mEq") = 0;
parameter Real NormalCalsUsed
"basal energy per kilogram of tissue [cal/min]";
parameter Physiolibrary.Types.PowerPerMass BasalCalsUsed
- "basal energy per kilogram of tissue"; //= NormalCalsUsed*(4.1864/60)
-
- //BasalCalsUsedPerKG(final displayUnit="kcal/min")
+ "basal energy per kilogram of tissue";
+//= NormalCalsUsed*(4.1864/60)
+ //BasalCalsUsedPerKG(final displayUnit="kcal/min")
parameter Physiolibrary.Types.Mass TissueSize "mass of tissue";
parameter Physiolibrary.Types.VolumeFlowRate DC(displayUnit="(mEq/min)/(mEq/ml)");
@@ -28242,8 +28163,7 @@
Physiolibrary.Types.Constants.MassConst size(k=TissueSize)
annotation (Placement(transformation(extent={{98,90},{90,98}})));
initial equation
- // BasalCalsUsedPerKG = NormalCalsUsed*(4186.4/60)/mass;
-
+// BasalCalsUsedPerKG = NormalCalsUsed*(4186.4/60)/mass;
equation
connect(Thyroid.y,HeatMetabolism_Skin. yBase) annotation (Line(
points={{70,80},{70,82}},
@@ -28460,13 +28380,13 @@
HeatMetabolism_Skin(data={{10.0+273.15,0.0,0},{37.1+273.15,1.0,0.12},{40.0+273.15,1.5,0},
{46.0+273.15,0.0,0}}),
PartCaloriesUse(nin=4));
- // Lactate(stateName="SkeletalMuscle-Lactate.Mass"),
+ // Lactate(stateName="SkeletalMuscle-Lactate.Mass"),
extends Physiolibrary.Icons.SkeletalMuscle;
Physiolibrary.Chemical.Interfaces.ChemicalPort_a fattyAcids
annotation (
Placement(transformation(extent={{-103,-78},{-83,-58}}),
- iconTransformation(extent={{90,-70},{110,-50}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
+ iconTransformation(extent={{90,-70},{110,-50}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
Physiomodel.Metabolism.CellularMetabolism
cellularMetabolism(
LacFractionData={{10,0,0},{100,0.3,0}},
@@ -28722,7 +28642,7 @@
extends Physiomodel.Metabolism.TissueMetabolism.TissueMetabolismBase(
TotalCaloriesUse(nin=4), HeatMetabolism_Skin(data={{10.0+273.15,0.0,0},{
37.1+273.15,1.0,0.12},{40.0+273.15,1.5,0},{46.0+273.15,0.0,0}}));
- // ,Lactate(stateName="RespiratoryMuscle-Lactate.Mass")
+ // ,Lactate(stateName="RespiratoryMuscle-Lactate.Mass")
extends Physiolibrary.Icons.RespiratoryMuscle;
Physiolibrary.Types.Constants.HeatFlowRateConst Shivering_Cals(k=0)
@@ -28766,7 +28686,7 @@
Physiolibrary.Chemical.Interfaces.ChemicalPort_a fattyAcids
annotation (
Placement(transformation(extent={{-101,-78},{-81,-58}}),
- iconTransformation(extent={{90,-70},{110,-50}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
+ iconTransformation(extent={{90,-70},{110,-50}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
Physiomodel.Metabolism.CellularMetabolism
cellularMetabolism(
LacFractionData={{10,0,0},{100,0.3,0}},
@@ -28971,7 +28891,7 @@
model Brain
extends Physiomodel.Metabolism.TissueMetabolism.TissueMetabolismBase;
- // (Lactate( stateName="Brain-Lactate.Mass"));
+ // (Lactate( stateName="Brain-Lactate.Mass"));
extends Physiolibrary.Icons.Brain;
Physiomodel.Metabolism.CellularMetabolism
cellularMetabolism(
@@ -29088,7 +29008,8 @@
model Kidney
extends Physiomodel.Metabolism.TissueMetabolism.TissueMetabolism(
- TotalCaloriesUse(nin=4));//, Lactate(stateName="Kidney-Lactate.Mass"));
+ TotalCaloriesUse(nin=4));
+ //, Lactate(stateName="Kidney-Lactate.Mass"));
extends Physiolibrary.Icons.Kidney;
Physiolibrary.Types.RealIO.MolarFlowRateInput PT_Na_Reab( displayUnit="mmol/min")
@@ -29157,13 +29078,14 @@
end HeartMuscle;
model GITract
- extends TissueMetabolismBase;//(Lactate(stateName="GITract-Lactate.Mass"));
+ extends TissueMetabolismBase;
+ //(Lactate(stateName="GITract-Lactate.Mass"));
extends Physiolibrary.Icons.GITract;
Physiolibrary.Chemical.Interfaces.ChemicalPort_a fattyAcids
annotation (
Placement(transformation(extent={{-103,-80},{-83,-60}}),
- iconTransformation(extent={{90,-70},{110,-50}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
+ iconTransformation(extent={{90,-70},{110,-50}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
Physiomodel.Metabolism.CellularMetabolism
cellularMetabolism(LacFractionData={{10,0,0},{100,0.3,0}},
PO2OnAerobicFractionData=PO2OnAerobicFractionData)
@@ -29313,24 +29235,25 @@
end GITract;
model Liver
- extends TissueMetabolismBase;//(Lactate(stateName="Liver-Lactate.Mass"));
+ extends TissueMetabolismBase;
+ //(Lactate(stateName="Liver-Lactate.Mass"));
extends Physiolibrary.Icons.Liver;
Physiolibrary.Chemical.Interfaces.ChemicalPort_a fattyAcids
annotation (
Placement(transformation(extent={{-103,-80},{-83,-60}}),
- iconTransformation(extent={{90,-70},{110,-50}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
+ iconTransformation(extent={{90,-70},{110,-50}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
Physiomodel.Metabolism.CellularMetabolism
cellularMetabolism(LacFractionData={{10,0,0},{100,0.3,0}},
PO2OnAerobicFractionData=PO2OnAerobicFractionData,
Ratio(start=0.087))
annotation (Placement(transformation(extent={{62,-76},{92,-46}})));
- // glucose(conc(start=1.43)),
+ // glucose(conc(start=1.43)),
NutrientDelivery_Fat fattyAcidsDelivery
annotation (Placement(transformation(extent={{-74,-80},{-54,-60}})));
- //(neededFlow(q(start=125.166135276543)))
+ //(neededFlow(q(start=125.166135276543)))
NutrientDelivery_2 glucoseDelivery
annotation (Placement(transformation(extent={{-72,-44},{-52,-24}})));
- //(delivery(start=991.319))
+ //(delivery(start=991.319))
Physiolibrary.Chemical.Interfaces.ChemicalPort_a GILumenCarbohydrates
annotation (Placement(transformation(extent={{-96,-58},{-76,-38}}),
iconTransformation(extent={{-30,90},{-10,110}})));
@@ -29540,22 +29463,25 @@
model Skin
extends TissueMetabolism(HeatMetabolism_Skin(data={{10.0+273.15,0.0,0},{29.5+273.15,
1.0,0.12},{40.0+273.15,1.5,0},{46.0+273.15,0.0,0}}));
- //(Lactate(stateName="Skin-Lactate.Mass"));
+ //(Lactate(stateName="Skin-Lactate.Mass"));
extends Physiolibrary.Icons.Skin;
end Skin;
model Bone
- extends TissueMetabolism;//(Lactate(stateName="Bone-Lactate.Mass"));
+ extends TissueMetabolism;
+ //(Lactate(stateName="Bone-Lactate.Mass"));
extends Physiolibrary.Icons.Bone;
end Bone;
model Fat
- extends TissueMetabolism;//(Lactate(stateName="Fat-Lactate.Mass"));
+ extends TissueMetabolism;
+ //(Lactate(stateName="Fat-Lactate.Mass"));
extends Physiolibrary.Icons.Fat;
end Fat;
model OtherTissue
- extends TissueMetabolism;//(Lactate(stateName="OtherTissue-Lactate.Mass"));
+ extends TissueMetabolism;
+ //(Lactate(stateName="OtherTissue-Lactate.Mass"));
extends Physiolibrary.Icons.OtherTissue;
end OtherTissue;
@@ -29569,7 +29495,7 @@
TissueSize=8.558,
DC(displayUnit="ml/min") = 3e-06)
annotation (Placement(transformation(extent={{66,88},{40,114}})));
- //initialLactateMass=3.00477
+ //initialLactateMass=3.00477
Physiomodel.Metabolism.TissueMetabolism.Brain brain(
PO2OnAerobicFractionData={{2.0,0.0,0},{20.0,1.0,0}},
NormalCalsUsed=193.709,
@@ -29577,7 +29503,7 @@
BasalCalsUsed=11.492766,
TissueSize=1.1572)
annotation (Placement(transformation(extent={{-60,86},{-32,114}})));
- //initialLactateMass=0.850703,
+ //initialLactateMass=0.850703,
Physiomodel.Metabolism.TissueMetabolism.Fat fat(
NormalCalsUsed=29.9658,
initialTissueO2(displayUnit="ml_STP") = 0.00096297146941496,
@@ -29585,8 +29511,7 @@
TissueSize=16.6,
DC(displayUnit="ml/min") = 4.5e-06)
annotation (Placement(transformation(extent={{66,-54},{40,-28}})));
- // initialLactateMass=2.73208,
-
+ // initialLactateMass=2.73208,
Physiomodel.Metabolism.TissueMetabolism.GITract gITract(
NormalCalsUsed=92.7945,
initialTissueO2(displayUnit="ml_STP") = 8.0401733111647e-05,
@@ -29594,7 +29519,7 @@
TissueSize=1.2716,
DC(displayUnit="ml/min") = 3e-06)
annotation (Placement(transformation(extent={{66,-22},{40,4}})));
- // initialLactateMass=0.934802)
+ // initialLactateMass=0.934802)
Physiomodel.Metabolism.TissueMetabolism.Kidney kidney(
PO2OnAerobicFractionData={{2.0,0.0,0},{20.0,1.0,0}},
NormalCalsUsed=36.8312,
@@ -29603,7 +29528,7 @@
BasalCalsUsed=9.90876,
TissueSize=0.2552)
annotation (Placement(transformation(extent={{66,10},{40,36}})));
- // initialLactateMass=0.187607)
+ // initialLactateMass=0.187607)
Physiomodel.Metabolism.TissueMetabolism.Liver liver(
PO2OnAerobicFractionData={{2.0,0.0,0},{10.0,1.0,0}},
NormalCalsUsed=147.338,
@@ -29612,7 +29537,7 @@
TissueSize=1.5422,
DC(displayUnit="ml/min") = 1.6666666666667e-06)
annotation (Placement(transformation(extent={{66,44},{40,70}})));
- // initialLactateMass=1.13373)
+ // initialLactateMass=1.13373)
Physiomodel.Metabolism.TissueMetabolism.OtherTissue otherTissue(
NormalCalsUsed=39.1075,
initialTissueO2(displayUnit="ml_STP") = 0.00022235966714886,
@@ -29620,7 +29545,7 @@
BasalCalsUsed=0.739668,
TissueSize=3.63)
annotation (Placement(transformation(extent={{66,-88},{40,-62}})));
- // initialLactateMass=2.66855)
+ // initialLactateMass=2.66855)
Physiomodel.Metabolism.TissueMetabolism.RespiratoryMuscle respiratoryMuscle(
NormalCalsUsed=6.27882,
DC(displayUnit="ml/min") = 4.5e-06,
@@ -29628,7 +29553,7 @@
BasalCalsUsed=0.13956,
TissueSize=3.08)
annotation (Placement(transformation(extent={{-60,52},{-32,80}})));
- // initialLactateMass=2.66976)
+ // initialLactateMass=2.66976)
Physiomodel.Metabolism.TissueMetabolism.RightHeartMuscle rightHeart(
NormalCalsUsed=2.34779,
DC(displayUnit="ml/min") = 5e-08,
@@ -29638,8 +29563,8 @@
BasalCalsUsed=4.1868,
TissueSize=0.0385)
annotation (Placement(transformation(extent={{-60,14},{-32,42}})));
- // ,Lactate(stateName="RightHeart-Lactate.Mass")
- // initialLactateMass=2.8e-002,
+ // ,Lactate(stateName="RightHeart-Lactate.Mass")
+ // initialLactateMass=2.8e-002,
Physiomodel.Metabolism.TissueMetabolism.LeftHeartMuscle leftHeart(
DC=18,
NormalCalsUsed=15.7067,
@@ -29649,8 +29574,8 @@
MotionCals(displayUnit="cal/min") = 1.67472,
HeatCals(displayUnit="cal/min") = 6.07086)
annotation (Placement(transformation(extent={{-58,-20},{-30,8}})));
- // ,Lactate(stateName="LeftHeart-Lactate.Mass")
- // initialLactateMass=0.169817,
+ // ,Lactate(stateName="LeftHeart-Lactate.Mass")
+ // initialLactateMass=0.169817,
Physiomodel.Metabolism.TissueMetabolism.SkeletalMuscleMetabolism skeletalMuscle(
PO2OnAerobicFractionData={{0.0,0.0,0},{15.0,0.2,0.04},{20.0,1.0,0}},
NormalCalsUsed=145.228,
@@ -29659,7 +29584,7 @@
BasalCalsUsed=0.362856,
TissueSize=27.4)
annotation (Placement(transformation(extent={{-58,-58},{-30,-30}})));
- // initialLactateMass=23.7503)
+ // initialLactateMass=23.7503)
Physiomodel.Metabolism.TissueMetabolism.Skin skin(
PO2OnAerobicFractionData={{2.0,0.0,0},{20.0,1.0,0}},
HeatMetabolism_Skin(data={{10.0+273.15,0.0,0},{29.5+273.15,1.0,0.12},{40.0+273.15,1.5,0},{
@@ -29670,30 +29595,31 @@
BasalCalsUsed=0.955986,
TissueSize=2.244)
annotation (Placement(transformation(extent={{-60,-96},{-32,-68}})));
- // initialLactateMass=1.64965)
+ // initialLactateMass=1.64965)
Physiolibrary.Chemical.Interfaces.ChemicalPort_a glucose
annotation (
Placement(transformation(extent={{-12,2},{8,22}}), iconTransformation(
- extent={{90,-20},{110,0}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
+ extent={{90,-20},{110,0}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
Physiolibrary.Chemical.Interfaces.ChemicalPort_a fattyAcids
annotation (
Placement(transformation(extent={{-10,-16},{10,4}}),
- iconTransformation(extent={{90,70},{110,90}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
+ iconTransformation(extent={{90,70},{110,90}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
Physiolibrary.Chemical.Interfaces.ChemicalPort_a ketoAcids
annotation (
Placement(transformation(extent={{-10,76},{10,96}}),
- iconTransformation(extent={{90,20},{110,40}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
+ iconTransformation(extent={{90,20},{110,40}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
Physiolibrary.Types.BusConnector busConnector annotation (Placement(
transformation(extent={{-92,106},{-72,126}}), iconTransformation(
extent={{-92,80},{-72,100}})));
Physiolibrary.Chemical.Interfaces.ChemicalPort_a GILumenCarbohydrates
annotation (Placement(transformation(extent={{70,72},{90,92}}),
iconTransformation(extent={{-70,-110},{-50,-90}})));
- //( conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
+ //( conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
Physiolibrary.Chemical.Interfaces.ChemicalPort_a lactate
annotation (
Placement(transformation(extent={{-10,56},{10,76}}),
- iconTransformation(extent={{90,-70},{110,-50}})));//(conc(final displayUnit="mEq/ml"), q(final displayUnit="mEq/min"))
+ iconTransformation(extent={{90,-70},{110,-50}})));
+ //(conc(final displayUnit="mEq/ml"), q(final displayUnit="mEq/min"))
Modelica.Blocks.Math.Feedback Heat annotation (Placement(transformation(
extent={{-5,-5},{5,5}},
rotation=180,
@@ -31598,17 +31524,17 @@
origin={-80,60})));
Physiolibrary.Types.MolarFlowRate consumption;
- // Real fromGILumen_conc;
- Physiolibrary.Types.MolarFlowRate delivery; //(start=300);
-
- equation
- // delivered.conc = neededFlow.conc + concChange;
- if ( -neededFlow.q > delivery) then
- consumption = delivery;
- FuelFractUseDelay = consumption/(-neededFlow.q);
- else
- consumption + neededFlow.q = 0;
- FuelFractUseDelay = 1;
+ // Real fromGILumen_conc;
+ Physiolibrary.Types.MolarFlowRate delivery;
+//(start=300);
+ equation
+// delivered.conc = neededFlow.conc + concChange;
+ if (-neededFlow.q) > delivery then
+ consumption = delivery;
+ FuelFractUseDelay = consumption / (-neededFlow.q);
+ else
+ consumption + neededFlow.q = 0;
+ FuelFractUseDelay = 1;
end if;
delivery = q_in.q + fromGILumen.q + GITUsed;
@@ -31616,12 +31542,10 @@
neededFlow.conc * solventFlow = - q_out.q;
- q_in.q = q_in.conc * solventFlow;
-
- //calculation of glucose concentration in gartro intestinal tract: fromGILumen.conc
+ q_in.q = q_in.conc * solventFlow;
+//calculation of glucose concentration in gartro intestinal tract: fromGILumen.conc
q_in.conc = fromGILumen.conc;
- // fromGILumen.q + q_in.conc * (solventFlow-HepaticArty) + GITUsed = fromGILumen.conc * (solventFlow-HepaticArty);
-
+// fromGILumen.q + q_in.conc * (solventFlow-HepaticArty) + GITUsed = fromGILumen.conc * (solventFlow-HepaticArty);
annotation (
Icon(coordinateSystem(preserveAspectRatio=true, extent={{-100,-100},{
100,100}}), graphics={
@@ -31654,7 +31578,7 @@
BasalCalsUsed=0.6126684,
TissueSize=8.558)
annotation (Placement(transformation(extent={{66,88},{40,114}})));
- //initialLactateMass=3.00477
+ //initialLactateMass=3.00477
Physiomodel.Metabolism.TissueMetabolism.Brain brain(
PO2OnAerobicFractionData={{2.0,0.0,0},{20.0,1.0,0}},
NormalCalsUsed=193.709,
@@ -31662,7 +31586,7 @@
BasalCalsUsed=12.4431696,
TissueSize=1.1572)
annotation (Placement(transformation(extent={{-60,86},{-32,114}})));
- //initialLactateMass=0.850703,
+ //initialLactateMass=0.850703,
Physiomodel.Metabolism.TissueMetabolism.Fat fat(
NormalCalsUsed=29.9658,
initialTissueO2(displayUnit="ml_STP") = 0.00096297146941496,
@@ -31670,8 +31594,7 @@
BasalCalsUsed=0.136071,
TissueSize=16.6)
annotation (Placement(transformation(extent={{66,-54},{40,-28}})));
- // initialLactateMass=2.73208,
-
+ // initialLactateMass=2.73208,
Physiomodel.Metabolism.TissueMetabolism.GITract gITract(
NormalCalsUsed=92.7945,
initialTissueO2(displayUnit="ml_STP") = 8.0401733111647e-05,
@@ -31679,7 +31602,7 @@
BasalCalsUsed=5.4191148,
TissueSize=1.2716)
annotation (Placement(transformation(extent={{66,-22},{40,4}})));
- // initialLactateMass=0.934802)
+ // initialLactateMass=0.934802)
Physiomodel.Metabolism.TissueMetabolism.Kidney kidney(
PO2OnAerobicFractionData={{2.0,0.0,0},{20.0,1.0,0}},
NormalCalsUsed=36.8312,
@@ -31688,7 +31611,7 @@
BasalCalsUsed=11.40903,
TissueSize=0.2552)
annotation (Placement(transformation(extent={{66,10},{40,36}})));
- // initialLactateMass=0.187607)
+ // initialLactateMass=0.187607)
Physiomodel.Metabolism.TissueMetabolism.Liver liver(
PO2OnAerobicFractionData={{2.0,0.0,0},{10.0,1.0,0}},
NormalCalsUsed=147.338,
@@ -31697,7 +31620,7 @@
BasalCalsUsed=7.1357028,
TissueSize=1.5422)
annotation (Placement(transformation(extent={{66,44},{40,70}})));
- // initialLactateMass=1.13373)
+ // initialLactateMass=1.13373)
Physiomodel.Metabolism.TissueMetabolism.OtherTissue otherTissue(
NormalCalsUsed=39.1075,
initialTissueO2(displayUnit="ml_STP") = 0.00022235966714886,
@@ -31705,7 +31628,7 @@
BasalCalsUsed=0.7933986,
TissueSize=3.63)
annotation (Placement(transformation(extent={{66,-88},{40,-62}})));
- // initialLactateMass=2.66855)
+ // initialLactateMass=2.66855)
Physiomodel.Metabolism.TissueMetabolism.RespiratoryMuscle respiratoryMuscle(
NormalCalsUsed=6.27882,
DC(displayUnit="ml/min") = 4.5e-06,
@@ -31713,7 +31636,7 @@
BasalCalsUsed=0.1584006,
TissueSize=3.08)
annotation (Placement(transformation(extent={{-60,52},{-32,80}})));
- // initialLactateMass=2.66976)
+ // initialLactateMass=2.66976)
Physiomodel.Metabolism.TissueMetabolism.RightHeartMuscle rightHeart(
NormalCalsUsed=2.34779,
DC(displayUnit="ml/min") = 5e-08,
@@ -31723,8 +31646,8 @@
BasalCalsUsed=5.5398342,
TissueSize=0.0385)
annotation (Placement(transformation(extent={{-60,14},{-32,42}})));
- // ,Lactate(stateName="RightHeart-Lactate.Mass")
- // initialLactateMass=2.8e-002,
+ // ,Lactate(stateName="RightHeart-Lactate.Mass")
+ // initialLactateMass=2.8e-002,
Physiomodel.Metabolism.TissueMetabolism.LeftHeartMuscle leftHeart(
DC=18,
NormalCalsUsed=15.7067,
@@ -31734,8 +31657,8 @@
BasalCalsUsed=6.0959808,
TissueSize=0.231)
annotation (Placement(transformation(extent={{-58,-20},{-30,8}})));
- // ,Lactate(stateName="LeftHeart-Lactate.Mass")
- // initialLactateMass=0.169817,
+ // ,Lactate(stateName="LeftHeart-Lactate.Mass")
+ // initialLactateMass=0.169817,
Physiomodel.Metabolism.TissueMetabolism.SkeletalMuscleMetabolism skeletalMuscle(
PO2OnAerobicFractionData={{0.0,0.0,0},{15.0,0.2,0.04},{20.0,1.0,0}},
NormalCalsUsed=145.228,
@@ -31744,7 +31667,7 @@
BasalCalsUsed=0.411702,
TissueSize=27.4)
annotation (Placement(transformation(extent={{-58,-58},{-30,-30}})));
- // initialLactateMass=23.7503)
+ // initialLactateMass=23.7503)
Physiomodel.Metabolism.TissueMetabolism.Skin skin(
PO2OnAerobicFractionData={{2.0,0.0,0},{20.0,1.0,0}},
HeatMetabolism_Skin(data={{10.0+273.15,0.0,0},{29.5+273.15,1.0,0.12},{40.0+273.15,1.5,0},{
@@ -31755,30 +31678,31 @@
BasalCalsUsed=0.9622662,
TissueSize=2.244)
annotation (Placement(transformation(extent={{-60,-96},{-32,-68}})));
- // initialLactateMass=1.64965)
+ // initialLactateMass=1.64965)
Physiolibrary.Chemical.Interfaces.ChemicalPort_a glucose
annotation (
Placement(transformation(extent={{-12,2},{8,22}}), iconTransformation(
- extent={{90,-20},{110,0}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
+ extent={{90,-20},{110,0}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
Physiolibrary.Chemical.Interfaces.ChemicalPort_a fattyAcids
annotation (
Placement(transformation(extent={{-10,-16},{10,4}}),
- iconTransformation(extent={{90,70},{110,90}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
+ iconTransformation(extent={{90,70},{110,90}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
Physiolibrary.Chemical.Interfaces.ChemicalPort_a ketoAcids
annotation (
Placement(transformation(extent={{-10,76},{10,96}}),
- iconTransformation(extent={{90,20},{110,40}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
+ iconTransformation(extent={{90,20},{110,40}}))); //(conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
Physiolibrary.Types.BusConnector busConnector annotation (Placement(
transformation(extent={{-92,106},{-72,126}}), iconTransformation(
extent={{-92,80},{-72,100}})));
Physiolibrary.Chemical.Interfaces.ChemicalPort_a GILumenCarbohydrates
annotation (Placement(transformation(extent={{70,72},{90,92}}),
iconTransformation(extent={{-70,-110},{-50,-90}})));
- //( conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
+ //( conc(final displayUnit="mg/ml"), q(final displayUnit="mg/min"))
Physiolibrary.Chemical.Interfaces.ChemicalPort_a lactate
annotation (
Placement(transformation(extent={{-10,56},{10,76}}),
- iconTransformation(extent={{90,-70},{110,-50}})));//(conc(final displayUnit="mEq/ml"), q(final displayUnit="mEq/min"))
+ iconTransformation(extent={{90,-70},{110,-50}})));
+ //(conc(final displayUnit="mEq/ml"), q(final displayUnit="mEq/min"))
Modelica.Blocks.Math.Feedback Heat annotation (Placement(transformation(
extent={{-5,-5},{5,5}},
rotation=180,
@@ -33352,7 +33276,7 @@
Physiomodel.Metabolism.TissueMetabolism.Metabolism
tissuesMetabolism
annotation (Placement(transformation(extent={{90,-6},{50,34}})));
- //(glucose(conc(start=1.05)))
+ //(glucose(conc(start=1.05)))
Physiolibrary.Types.BusConnector busConnector annotation (Placement(
transformation(extent={{60,60},{100,100}}), iconTransformation(extent=
{{60,60},{100,100}})));
@@ -33605,9 +33529,9 @@
T.MolarFlowRate brain_O2Use(varName="Brain-Flow.O2Use", storeUnit="ml_STP/min")
"Brain oxygen consumption. "
annotation (Placement(transformation(extent={{-36,70},{-30,76}})));
- //T.Fraction BrainFuel_FractUseDelay(varName="Brain-Fuel.FractUseDelay")
- // "Brain deficiency of nutrients needed by cells metabolism. "
- //annotation (Placement(transformation(extent={{-36,80},{-30,86}})));
+ //T.Fraction BrainFuel_FractUseDelay(varName="Brain-Fuel.FractUseDelay")
+ // "Brain deficiency of nutrients needed by cells metabolism. "
+ //annotation (Placement(transformation(extent={{-36,80},{-30,86}})));
T.Fraction brain_Fuel_FractUseDelay(varName="Brain-Fuel.FractUseDelay")
"Brain deficiency of nutrients needed by cells metabolism. "
annotation (Placement(transformation(extent={{-36,90},{-30,96}})));
@@ -33659,9 +33583,9 @@
T.Concentration Glucose(varName="GlucosePool.[Glucose(mG/dL)]", storeUnit="mg/dl", unitConversions = Substances.Glucose.unitConversions)
"Extracellular glucose concentration."
annotation (Placement(transformation(extent={{-36,260},{-30,266}})));
- //Variable KA_Change(varName="KAPool.Change")
- // "Change of keto-acids extracellular mass."
- //annotation (Placement(transformation(extent={{-34,270},{-28,276}})));
+ //Variable KA_Change(varName="KAPool.Change")
+ // "Change of keto-acids extracellular mass."
+ //annotation (Placement(transformation(extent={{-34,270},{-28,276}})));
T.AmountOfSubstance KAPool_Osmoles(varName="KAPool.Osmoles", storeUnit="mosm")
"Keto-acids extracellular osmoles."
annotation (Placement(transformation(extent={{-36,280},{-30,286}})));
@@ -33876,7 +33800,7 @@
"Triiodothyronine and thyroxine effect to metabolism."
annotation (Placement(transformation(extent={{-82,4},{-76,10}})));
- ////to debug:
+////to debug:
T.Concentration FAPool(varName="FAPool.[FA(mMol/L)]")
"Fatty-acids extracellular concentration."
annotation (Placement(transformation(extent={{-100,-42},{-94,-36}})));
@@ -34937,9 +34861,9 @@
end IO_Bus;
model LiverMetabolism2
- // import QHP = Physiomodel;
+ // import QHP = Physiomodel;
extends Physiolibrary.Icons.Liver;
- //outer parameter Real EFC_Vol(final displayUnit="ml");
+ //outer parameter Real EFC_Vol(final displayUnit="ml");
Physiolibrary.Chemical.Components.Substance Glycogen(solute_start=120e-6/
Substances.Glucose.mw)
"amount of bound glucose molecules in glycogen (initial mass 120 mg of bound glucose molecules)"
@@ -35706,14 +35630,14 @@
extent=[-10,-110; 10,-90], Placement(transformation(extent={{90,20},{
110,40}}), iconTransformation(extent={{90,20},{110,40}})));
- parameter Modelica.SIunits.StoichiometricNumber a=1
+ parameter Modelica.Units.SI.StoichiometricNumber a=1
"Stoichiometric number of solute A";
- parameter Modelica.SIunits.StoichiometricNumber b=1
+ parameter Modelica.Units.SI.StoichiometricNumber b=1
"Stoichiometric number of solute B";
equation
a*A.q + b*B.q = 0;
- a*A.conc = b*B.conc; // sound strange, dissociation constand should be here
+ a*A.conc = b*B.conc; // sound strange, dissociation constand should be here
annotation (
Icon(coordinateSystem(preserveAspectRatio=true, extent={{-100,-100},{
100,100}}), graphics={
@@ -35762,10 +35686,9 @@
rotation=270,
origin={60,40})));
- // parameter Modelica.SIunits.StoichiometricNumber a_in=1 "Stoichiometric number of solute A";
- // parameter Modelica.SIunits.StoichiometricNumber a_out=1 "Stoichiometric number of solute B";
- // parameter Modelica.SIunits.StoichiometricNumber a_out2=1 "Stoichiometric number of solute C";
-
+// parameter Modelica.SIunits.StoichiometricNumber a_in=1 "Stoichiometric number of solute A";
+ // parameter Modelica.SIunits.StoichiometricNumber a_out=1 "Stoichiometric number of solute B";
+ // parameter Modelica.SIunits.StoichiometricNumber a_out2=1 "Stoichiometric number of solute C";
equation
q_out.q + coef*q_in.q = 0;
q_out2.q + coef2*q_in.q = 0;
@@ -35867,8 +35790,7 @@
model Insulin
extends Physiolibrary.Icons.Pancreas;
- // outer parameter Real ECF_Vol(final displayUnit="ml");
-
+ // outer parameter Real ECF_Vol(final displayUnit="ml");
Insulin_Physiolibrary.Chemical.Components.Substance InsulinPool(
stateName="InsulinPool.Mass",
useNormalizedVolume=false,
@@ -36099,8 +36021,7 @@
model Glucagon
extends Physiolibrary.Icons.Pancreas;
- // outer parameter Real ECF_Vol(final displayUnit="ml");
-
+ // outer parameter Real ECF_Vol(final displayUnit="ml");
Modelica.Blocks.Math.Gain gain(k=0.05/60)
annotation (Placement(transformation(extent={{-44,-32},{-38,-26}})));
Physiolibrary.Chemical.Sources.UnlimitedSolutePumpOut clearance(
@@ -36264,8 +36185,7 @@
model ADH
extends Physiolibrary.Icons.Hypophysis;
- // outer parameter Real ECF_Vol(final displayUnit="ml");
-
+ // outer parameter Real ECF_Vol(final displayUnit="ml");
constant Physiolibrary.Types.Volume NormalMedullaVolume = 31e-6
"Normal kidney medulla volume (31ml)";
constant Physiolibrary.Types.Volume NormalECFVolume = 14.3610097207699e-3
@@ -36548,7 +36468,7 @@
annotation (Placement(
transformation(
extent={{-10,-10},{10,10}},
- origin={-44,-2}))); //(toAnotherUnitCoef=1000, unitsString="pmol/l")
+ origin={-44,-2}))); //(toAnotherUnitCoef=1000, unitsString="pmol/l")
Physiolibrary.Chemical.Sources.UnlimitedSolutePumpOut Degradation(
useSoluteFlowInput=true)
annotation (Placement(transformation(extent={{0,-68},{20,-48}})));
@@ -36668,8 +36588,7 @@
model ANP
extends Physiolibrary.Icons.Heart;
- // outer parameter Real ECF_Vol(final displayUnit="ml");
-
+ // outer parameter Real ECF_Vol(final displayUnit="ml");
Physiolibrary.Chemical.Components.Substance ANPPool(
stateName="ANPPool.Mass",
useNormalizedVolume=false,
@@ -36892,16 +36811,14 @@
model EPO
extends Physiolibrary.Icons.Kidney;
- // outer parameter Real ECF_Vol(final displayUnit="ml");
-
+ // outer parameter Real ECF_Vol(final displayUnit="ml");
EPO_Physiolibrary.Chemical.Components.Substance EPOPool(
stateName="EPOPool.Mass",
useNormalizedVolume=false,
solute_start=(20e3*14361.0097207699*0.4)*1e-6)
"default = 20 iu/l = 20e3 iu/m3 "
annotation (Placement(transformation(extent={{-50,-54},{-30,-34}})));
- //NominalSolute=Substances.Erythropoietin.molpIU
-
+ //NominalSolute=Substances.Erythropoietin.molpIU
EPO_Physiolibrary.Chemical.Sources.UnlimitedSolutePump secretion(
useSoluteFlowInput=true)
annotation (Placement(transformation(extent={{76,-54},{56,-34}})));
@@ -36917,7 +36834,7 @@
annotation (Placement(
transformation(
extent={{-10,-10},{10,10}},
- origin={-38,-32}))); //(toAnotherUnitCoef=1000, unitsString="U/l")
+ origin={-38,-32}))); //(toAnotherUnitCoef=1000, unitsString="U/l")
EPO_Physiolibrary.Chemical.Sources.UnlimitedSolutePumpOut Clearance(
useSoluteFlowInput=true)
annotation (Placement(transformation(extent={{-30,-98},{-10,-78}})));
@@ -37055,8 +36972,7 @@
model Thyroxine
extends Physiolibrary.Icons.ThyroidGland;
- // outer parameter Real ECF_Vol(final displayUnit="ml");
-
+ // outer parameter Real ECF_Vol(final displayUnit="ml");
Physiolibrary.Chemical.Components.Substance ThyroidPool(
stateName="ThyroidPool.Mass",
useNormalizedVolume=false,
@@ -37172,8 +37088,7 @@
model Thyrotropin
extends Physiolibrary.Icons.Hypophysis;
- // outer parameter Real ECF_Vol(final displayUnit="ml");
-
+ // outer parameter Real ECF_Vol(final displayUnit="ml");
TSH_Physiolibrary.Types.Constants.ConcentrationConst hormoneConstant(k=4)
"4 uU/ml"
annotation (Placement(transformation(extent={{-8,30},{-20,42}})));
@@ -37256,8 +37171,7 @@
model Renin
extends Physiolibrary.Icons.Kidney;
- // outer parameter Real ECF_Vol(final displayUnit="ml");
-
+ // outer parameter Real ECF_Vol(final displayUnit="ml");
Renin_Physiolibrary.Chemical.Components.Substance ReninPool(
stateName="ReninPool.Mass",
useNormalizedVolume=false,
@@ -37279,7 +37193,7 @@
Renin_Physiolibrary.Chemical.Sensors.ConcentrationMeasure concentrationMeasure
annotation (Placement(transformation(
extent={{-10,-10},{10,10}},
- origin={-62,-32}))); //(unitsString="GU/ml")
+ origin={-62,-32}))); //(unitsString="GU/ml")
Renin_Physiolibrary.Chemical.Sources.UnlimitedSolutePumpOut Clearance(
useSoluteFlowInput=true)
annotation (Placement(transformation(extent={{-30,-98},{-10,-78}})));
@@ -37649,12 +37563,11 @@
end Renin;
model Angiotensine2
- // extends Library.Interfaces.Tissues.Liver;
- // extends Library.Interfaces.Tissues.Kidney;
+ // extends Library.Interfaces.Tissues.Liver;
+ // extends Library.Interfaces.Tissues.Kidney;
extends Physiolibrary.Icons.Lungs;
- // outer parameter Real ECF_Vol(final displayUnit="ml");
-
+// outer parameter Real ECF_Vol(final displayUnit="ml");
Physiolibrary.Blocks.Factors.Normalization A2 annotation (Placement(transformation(extent={{-2,-2},{18,18}})));
Physiolibrary.Types.Constants.ConcentrationConst A2CONC(k=0.3333e-9/
Substances.AngiotensinII.mw) "0.3333 pg/ml = 0.3333*1e-15/1e-6 kg/m3"
@@ -37724,8 +37637,7 @@
model Leptin
extends Physiolibrary.Icons.Fat;
- // outer parameter Real ECF_Vol(final displayUnit="ml");
-
+ // outer parameter Real ECF_Vol(final displayUnit="ml");
Modelica.Blocks.Math.Gain gain( y(unit="mol/s"), k(unit=
"1/s") = 0.01/60)
annotation (Placement(transformation(extent={{-44,-32},{-38,-26}})));
@@ -37908,7 +37820,7 @@
"initial pool of PT gland";
initial equation
- // PT_pool = initialPT;
+// PT_pool = initialPT;
der( PT_pool) = 0;
equation
der( PT_pool) = (1-PT_pool)*0.01*(0.85*T64m+0.15) - PT_pool*0.01*(0.85*T64p+0.15);
@@ -38095,11 +38007,11 @@
iconTransformation(extent={{60,-20},{100,20}})));
initial equation
- //One_alpha_OH = 126;
+//One_alpha_OH = 126;
der( One_alpha_OH) = 0;
equation
der( One_alpha_OH) = 6.3*H(PTH, 1.90371, 1.549,0.1112)*( if PO4<=16.8/14 then 1 else H_minus(PO4,1.525,1.302,8.252,-1.78249)) -0.05 *One_alpha_OH;
- //H7,9:alpha was calculated,
+//H7,9:alpha was calculated,
end One_alpha_OH;
One_alpha_OH one_alpha_OH
annotation (Placement(transformation(extent={{-26,16},{-6,36}})));
@@ -38203,8 +38115,7 @@
model Epinephrine
extends Physiolibrary.Icons.AdrenalGland;
- // outer parameter Real ECF_Vol(final displayUnit="ml");
-
+ // outer parameter Real ECF_Vol(final displayUnit="ml");
Physiolibrary.Chemical.Components.Substance EpiPool(
stateName="EpiPool.Mass",
useNormalizedVolume=false,
@@ -38306,8 +38217,7 @@
model Norepinephrine
extends Physiolibrary.Icons.AdrenalGland;
- // outer parameter Real ECF_Vol(final displayUnit="ml");
-
+ // outer parameter Real ECF_Vol(final displayUnit="ml");
Physiolibrary.Chemical.Components.Substance NEPool(
stateName="NEPool.Mass",
useNormalizedVolume=false,
@@ -39157,19 +39067,19 @@
model OvarianCycle
Modelica.StateGraph.InitialStepWithSignal FOLLICULAR(
- nIn=3)
+ nOut = 1, nIn = 3)
annotation (Placement(transformation(extent={{-94,18},{-74,38}})));
Modelica.StateGraph.TransitionWithSignal Follicle
annotation (Placement(transformation(extent={{-68,38},{-48,18}})));
Modelica.StateGraph.StepWithSignal OVULATORY(
- nOut=2)
+ nIn = 1, nOut = 2)
annotation (Placement(transformation(extent={{-42,18},{-22,38}})));
Modelica.StateGraph.TransitionWithSignal Ovulation
annotation (Placement(transformation(extent={{-10,38},{10,18}})));
Modelica.StateGraph.StepWithSignal LUTEAL(
- nOut=2)
+ nIn = 1, nOut = 2)
annotation (Placement(transformation(extent={{16,18},{36,38}})));
- Modelica.StateGraph.StepWithSignal IDLE
+ Modelica.StateGraph.StepWithSignal IDLE(nIn = 1, nOut = 1)
annotation (Placement(transformation(extent={{74,18},{94,38}})));
Modelica.StateGraph.TransitionWithSignal CorpusLuteum
annotation (Placement(transformation(extent={{48,38},{68,18}})));
@@ -39576,8 +39486,7 @@
end Progesterone;
model LH
- //Library.Interfaces.Tissues.Hypothalamus;
-
+ //Library.Interfaces.Tissues.Hypothalamus;
Physiolibrary.Chemical.Components.Substance LH_Circulating_Mass(
q_out(conc(final displayUnit="IU/ml"), q(final displayUnit="IU/min")),
stateName="LH-Circulating.Mass",
@@ -40611,10 +40520,9 @@
T.Fraction BetaPool_Effect(varName="BetaPool.Effect")
"Humoral and neural effect on beta receptor stimulation."
annotation (Placement(transformation(extent={{-40,38},{-34,44}})));
- //Variable NephronADH_conc(varName="NephronADH.[ADH]")
- // "Nephron vasopresine concentration. [ng/l]"
- //annotation (Placement(transformation(extent={{-36,170},{-30,176}})));
-
+ //Variable NephronADH_conc(varName="NephronADH.[ADH]")
+ // "Nephron vasopresine concentration. [ng/l]"
+ //annotation (Placement(transformation(extent={{-36,170},{-30,176}})));
T.Concentration Epinephrine(varName="EpiPool.[Epi(nMol/L)]", storeUnit=
"nmol/l") "Epinephrine extracelular concentration."
annotation (Placement(transformation(extent={{-40,-28},{-34,-22}})));
@@ -40889,14 +40797,14 @@
redeclare model extends Variables
extends Physiolibrary.Icons.Hormones;
- //Variable NephronADH_conc(varName="NephronADH.[ADH]")
- // "Nephron vasopresine concentration. [ng/l]"
- //annotation (Placement(transformation(extent={{-36,170},{-30,176}})));
-
+//Variable NephronADH_conc(varName="NephronADH.[ADH]")
+ // "Nephron vasopresine concentration. [ng/l]"
+ //annotation (Placement(transformation(extent={{-36,170},{-30,176}})));
T.Concentration EPOPool(
varName="EPOPool.[EPO]", storeUnit="iu/l")
"Erythropoetin concentration."
- annotation (Placement(transformation(extent={{-40,48},{-34,54}})));
+ annotation (Placement(transformation(extent={{-40,48},{-34,54}})));
+
equation
connect(EPOPool.y, busConnector.Erythropoietin) annotation (Line(
@@ -41116,23 +41024,20 @@
import Modelica.Math;
Physiolibrary.Types.Concentration betaX "buffer value of blood";
- //(unit = "mEq/l")
+ //(unit = "mEq/l")
Physiolibrary.Types.Concentration betaP "buffer value of plasma";
- //(unit = "mEq/l")
-
+ //(unit = "mEq/l")
Physiolibrary.Types.Concentration _cTH
"total concentration of tiratable hydrogen ions";
- //( displayUnit="mEq/l")
+ //( displayUnit="mEq/l")
Physiolibrary.Types.Concentration _BE "base excess";
- //( displayUnit="mEq/l")
-
+ //( displayUnit="mEq/l")
Physiolibrary.Types.Concentration _BEox
"base excess in fully oxygenated blood";
- //( displayUnit="mEq/l")
+ //( displayUnit="mEq/l")
Physiolibrary.Types.Concentration _cTHox
"total concentration of tiratable hydrogen ions in fully oxygenated blood";
- //( displayUnit="mEq/l")
-
+ //( displayUnit="mEq/l")
Physiolibrary.Types.RealIO.pHOutput pH( start=7.4) "plasma pH" annotation (Placement(transformation(
extent={{30,-98},{70,-58}}), iconTransformation(
extent={{-20,-20},{20,20}},
@@ -41190,26 +41095,20 @@
rotation=180,
origin={90,-80})));
- equation
- //Van Slyke equation: - T=37
- _cTH = -(1-ctHb/cHb)*((cHCO3-24.5)+betaX*(pH-7.4));
- // _cTH = -(1-(1-0.57)*(0.44))*((cHCO3-24.5)+betaX*(pH-7.4));
-
- betaX = 2.3*ctHb + betaP;
- betaP = 8*ctAlb + 0.075*ctGlb + 0.309*ctPO4;
-
- // _betaX = 0.11 * plasmaProteinConc(displayUnit="g/l"); //mmol/l
-
- //Haldane effect:
- _cTH = homotopy(_cTHox - 0.3 * (1-sO2),_cTHox);
-
- //base excess
- _BEox = - _cTHox;
- _BE = - _cTH;
-
- //erythrocytes:
- pH_ery = homotopy(7.19 + 0.77*(pH-7.4) + 0.035*(1-sO2),7.19 + 0.77*(pH-7.4));
-
+ equation
+//Van Slyke equation: - T=37
+ _cTH = -(1 - ctHb / cHb) * (cHCO3 - 24.5 + betaX * (pH - 7.4));
+// _cTH = -(1-(1-0.57)*(0.44))*((cHCO3-24.5)+betaX*(pH-7.4));
+ betaX = 2.3 * ctHb + betaP;
+ betaP = 8*ctAlb + 0.075*ctGlb + 0.309*ctPO4;
+// _betaX = 0.11 * plasmaProteinConc(displayUnit="g/l"); //mmol/l
+//Haldane effect:
+ _cTH = homotopy(_cTHox - 0.3 * (1 - sO2), _cTHox);
+//base excess
+ _BEox = -_cTHox;
+ _BE = - _cTH;
+//erythrocytes:
+ pH_ery = homotopy(7.19 + 0.77 * (pH - 7.4) + 0.035 * (1 - sO2), 7.19 + 0.77 * (pH - 7.4));
annotation (
preferredView = text,
Documentation(info="
@@ -41268,9 +41167,7 @@
origin={-80,100})));
equation
_BEox = BEox / Modelica.Constants.F;
-
- // cTHox.conc*1000 = _cTHox;
-
+// cTHox.conc*1000 = _cTHox;
annotation ( Documentation(revisions="
@@ -41297,8 +41194,7 @@
end BloodPhMeassure;
model NormalSID
- // extends BloodPhBase;
-
+ // extends BloodPhBase;
Physiolibrary.Types.RealIO.ConcentrationInput ctAlb(
final displayUnit="mmol/l")
"concentration of total plasma albumins(dofault=0.65mmol/l)" annotation (Placement(transformation(extent={{-120,72},
@@ -41524,10 +41420,8 @@
points={{-100,62},{-69,62},{-69,64},{-38,64}},
color={0,0,127},
smooth=Smooth.None));
-
- // NBB = (1-Hct)( [HCO3]P(40mmHg,7.4) + tAlbP*-zAlbP(7.4) + tGlbP*-zGlb(7.4) + tPO4P*-zPO4(7.4) + tORGP*-zORG(7.4) + .. ) + Hct ( [HCO3]E(40mmHg,NpHE) + tHb*-zHb(sO2,sCO2,NpHE) + tPO4E*-zPO4(NpHE) ..)
- // NBB = 0;//_cTHox;
-
+// NBB = (1-Hct)( [HCO3]P(40mmHg,7.4) + tAlbP*-zAlbP(7.4) + tGlbP*-zGlb(7.4) + tPO4P*-zPO4(7.4) + tORGP*-zORG(7.4) + .. ) + Hct ( [HCO3]E(40mmHg,NpHE) + tHb*-zHb(sO2,sCO2,NpHE) + tPO4E*-zPO4(NpHE) ..)
+// NBB = 0;//_cTHox;
connect(n_sO2.y, haemoglobin.sO2) annotation (Line(
points={{-72.5,-60},{-38,-60}},
color={0,0,127},
@@ -41588,16 +41482,13 @@
parameter Real pKaHPO4=11.78 "HPO4 acid dissociation constant";
parameter Real pKaKA=4.3
"average Keto-acids acid dissociation constant";
- //constant Real ml2l = 0.001;
-
+ //constant Real ml2l = 0.001;
Real zPO4 "average charge of one phosphate in urine";
- //(displayUnit="Eq/Mol")
-
+ //(displayUnit="Eq/Mol")
Real zKA "average charge of one keto acid in purine";
- //(displayUnit="Eq/Mol")
-
+ //(displayUnit="Eq/Mol")
Real zCO2_qCO2 "charge outflow in (bi)carbonates to urine";
- //(displayUnit="Eq/min")
+ //(displayUnit="Eq/min")
Physiolibrary.Types.RealIO.PressureInput
pCO2(displayUnit="mmHg")
"partial CO2 pressure" annotation (Placement(transformation(
@@ -41660,8 +41551,7 @@
zCO2_qCO2 = -CO2_solubility*pCO2*(10^(PHU-pKaCO2))*(qH2O);
- qX + 1*qNa + 1*qK + 2*qMg + 2*qCa + 1*qNH4 + (-1)*qCl + (-2)*qSO4 + zPO4*qPO4 + zKA*qKA + zCO2_qCO2 = 0; // electroneutrality of urine outflow
-
+ qX + 1*qNa + 1*qK + 2*qMg + 2*qCa + 1*qNH4 + (-1)*qCl + (-2)*qSO4 + zPO4*qPO4 + zKA*qKA + zCO2_qCO2 = 0; // electroneutrality of urine outflow
annotation (
Icon(coordinateSystem(preserveAspectRatio=true, extent={{-100,-120},{100,100}}),
graphics={
@@ -42141,7 +42031,7 @@
"tissue interstitial bicarbonate concentration"
annotation (Placement(transformation(extent={{-106,-4},
{-90,12}}), iconTransformation(extent={{-10,-10},{10,10}},
- origin={-90,-40}))); //(final displayUnit="mEq/ml")
+ origin={-90,-40}))); //(final displayUnit="mEq/ml")
Physiolibrary.Types.RealIO.pHOutput pH_interstitial
"tissue interstitial pH" annotation (Placement(
transformation(extent={{24,-16},{44,4}}), iconTransformation(
@@ -43459,8 +43349,7 @@
{34,106}}), iconTransformation(extent={{-10,-10},{10,10}},
origin={-90,40}))); //(final displayUnit="mEq/l")
equation
-
- //base excess
+//base excess
ctHox_lessLactate/Modelica.Constants.F + cLactate = _cTHox;
// cTH_interstitial = -(1-((ctHb_ery/3)/43))*((cHCO3_interstitial-24.5)+betaX*(pH_interstitial-7.42));
@@ -45214,7 +45103,8 @@
Placement(transformation(extent={{-98,66},{-84,78}}),
iconTransformation(extent={{90,-10},{110,10}})));
- Physiolibrary.Types.MolarFlowRate cTHfromTissues;//(displayUnit="mEq/min");
+ Physiolibrary.Types.MolarFlowRate cTHfromTissues;
+ //(displayUnit="mEq/min");
equation
cTHfromTissues = -( q_out.q + q_in.q);
@@ -45735,8 +45625,7 @@
{40,2}}), iconTransformation(extent={{-10,-10},{10,10}},
origin={-90,-72})));
equation
-
- //base excess
+//base excess
ctHox = _cTHox;
// cTH_interstitial = -(1-((ctHb_ery/3)/43))*((cHCO3_interstitial-24.5)+betaX*(pH_interstitial-7.42));
@@ -45916,7 +45805,7 @@
origin={24,24})));
Physiolibrary.Types.RealIO.MolarFlowRateInput LactateFromMetabolism annotation (Placement(transformation(extent={{-82,-16},
{-70,-4}}), iconTransformation(extent={{-10,-10},{10,10}}, rotation=180, origin={90,-60})));
- //(displayUnit="mEq/min")
+ //(displayUnit="mEq/min")
Physiolibrary.Chemical.Interfaces.ChemicalPort_b q_out annotation (
Placement(transformation(extent={{74,14},{94,34}}),
iconTransformation(extent={{90,-10},{110,10}})));
@@ -45928,7 +45817,7 @@
transformation(extent={{-4,-4},{4,4}}), iconTransformation(
extent={{-20,-20},{20,20}},
rotation=270,
- origin={20,-120}))); //(displayUnit="mEq/l")
+ origin={20,-120}))); //(displayUnit="mEq/l")
Physiolibrary.Chemical.Interfaces.ChemicalPort_a q_in annotation (
Placement(transformation(extent={{-16,16},{0,32}}),
iconTransformation(extent={{-110,-10},{-90,10}})));
@@ -46202,16 +46091,15 @@
Physiolibrary.Types.Concentration cdCO2;
Physiolibrary.Types.VolumeDensityOfCharge y2;
- equation
-
- //Henderson-Hasselbalch equation:
- pK = 6.1 + (-0.0026)*(T-310.15);
- aCO2 = 0.00023 * 10^(-0.0092*(T-310.15)); //solubility depends on temperature
- cdCO2 = aCO2*pCO2;
- y = -cdCO2*10^(pH-pK)* Modelica.Constants.F; //bicarbonate
-
- y2 = -cdCO2*10^(pH-6.35)* Modelica.Constants.F; //bicarbonate
-
+ equation
+//Henderson-Hasselbalch equation:
+ pK = 6.1 + (-0.0026) * (T - 310.15);
+ aCO2 = 0.00023 * 10^(-0.0092*(T-310.15));
+//solubility depends on temperature
+ cdCO2 = aCO2 * pCO2;
+ y = -cdCO2*10^(pH-pK)* Modelica.Constants.F;
+//bicarbonate
+ y2 = -cdCO2 * 10 ^ (pH - 6.35) * Modelica.Constants.F; //bicarbonate
annotation (Icon(graphics={Text(
extent={{-79,-27},{79,27}},
lineColor={0,0,255},
@@ -46229,7 +46117,7 @@
Physiolibrary.Types.MassConcentration tAlb_mass;
parameter Physiolibrary.Types.MolarMass Alb_MolarMass(displayUnit="g/mol") = 66.463;
equation
- tAlb_mass=tAlb*Alb_MolarMass; //kg/m3 = g/l
+ tAlb_mass=tAlb*Alb_MolarMass; //kg/m3 = g/l
y = - tAlb_mass * ( 0.123 * pH - 0.631) * Modelica.Constants.F;
//where [ Albx- ] is in mEq / L and [ Albumin ] is in g / dL. Therefore, at pH 7.40, the charge contributed by 4.40 g / dL of albumin is approximately -12.3 mEq / L.
@@ -46254,7 +46142,7 @@
ctGlb
annotation (Placement(transformation(extent={{-120,-40},{-80,0}})));
equation
- // y=-ctGlb*(0.075/0.81)*(pH-3.6); //sklon -0.075/0.81 so siggaarda, pri 7.4 chcem -10 aby to pokrilo nastavene SID u Colemana
+// y=-ctGlb*(0.075/0.81)*(pH-3.6); //sklon -0.075/0.81 so siggaarda, pri 7.4 chcem -10 aby to pokrilo nastavene SID u Colemana
y=-ctGlb*((0.075/0.77)*(pH-7.4) + 2.5/28)* Modelica.Constants.F; //chcem naboj -2.5 mEq/l pre normalnu koncentraciu 28 mg/l a pH=7.4
// y=-ctGlb*((0.075/0.77)*(pH-7.4) + 15.6/28)* Modelica.Constants.F; //chcem naboj -15.6 mEq/l pre normalnu koncentraciu 28 mg/l a pH=7.4
annotation (Icon(graphics={Text(
@@ -46342,18 +46230,21 @@
T
annotation (Placement(transformation(extent={{-120,-120},{-80,-80}})));
equation
- aCO2 = 0.23 * 10^(-0.0092*(T-37)); //solubility depends on temperature
- cdCO2 = aCO2*pCO2*0.001; //from mmol to mol
-
- carbaminohaemoglobin = tHb_E*cdCO2*( ( H*H*((1-sO2)/(KzO*KcO)+sO2/(KzR*KcR)) + H*((1-sO2)/KcO + sO2/KcR) + cdCO2) / ((H*H/(KzR*KcR) + H/KcR + cdCO2)*(H*H/(KzO*KcO) + H/KcO + cdCO2))); //, where cdCO2=free disolved CO2 concentration in erythrocytes in mol/l
- sCO2 = carbaminohaemoglobin/ tHb_E;
+ aCO2 = 0.23 * 10^(-0.0092*(T-37));
+//solubility depends on temperature
+ cdCO2 = aCO2 * pCO2 * 0.001;
+//from mmol to mol
+ carbaminohaemoglobin = tHb_E * cdCO2 * ((H * H * ((1 - sO2) / (KzO * KcO) + sO2 / (KzR * KcR)) + H * ((1 - sO2) / KcO + sO2 / KcR) + cdCO2) / ((H * H / (KzR * KcR) + H / KcR + cdCO2) * (H * H / (KzO * KcO) + H / KcO + cdCO2)));
+//, where cdCO2=free disolved CO2 concentration in erythrocytes in mol/l
+ sCO2 = carbaminohaemoglobin / tHb_E;
y = -tHb_E*(sO2*(1-sCO2)*zOxyHb + (1-sO2)*(1-sCO2)*zDeoxyHb + sO2*sCO2*zOxyCarbaminoHb + (1-sO2)*sCO2*zDeoxyCarbaminoHb)* Modelica.Constants.F;
zOxyHb = betaOxyHb * (pH-pIo);
zDeoxyHb = betaDeoxyHb * (pH-pIr);
- zOxyCarbaminoHb = zOxyHb + (1+2*10^(pKzO-pH))/(1+10^(pKzO-pH)); //= zOxyHb + -[O2HbNH2]/([O2HbNH2]+[O2HbNH3+]) + -2[O2HbNH3+]/([O2HbNH2]+[O2HbNH3+]), where pKZO=8.08
- zDeoxyCarbaminoHb = zDeoxyHb + (1+2*10^(pKzR-pH))/(1+10^(pKzR-pH)); //betaOxyHb(pH) * (pH-pIO) + -[HbNH2]/([HbNH2]+[HbNH3+]) + -2[HbNH3+]/([HbNH2]+[HbNH3+]), where pKZR=7.14
+ zOxyCarbaminoHb = zOxyHb + (1+2*10^(pKzO-pH))/(1+10^(pKzO-pH));
+//= zOxyHb + -[O2HbNH2]/([O2HbNH2]+[O2HbNH3+]) + -2[O2HbNH3+]/([O2HbNH2]+[O2HbNH3+]), where pKZO=8.08
+ zDeoxyCarbaminoHb = zDeoxyHb + (1 + 2 * 10 ^ (pKzR - pH)) / (1 + 10 ^ (pKzR - pH)); //betaOxyHb(pH) * (pH-pIO) + -[HbNH2]/([HbNH2]+[HbNH3+]) + -2[HbNH3+]/([HbNH2]+[HbNH3+]), where pKZR=7.14
annotation (Icon(graphics={Text(
extent={{-79,-27},{79,27}},
lineColor={0,0,255},
@@ -46399,9 +46290,8 @@
transformation(extent={{-44,72},{-4,112}}), iconTransformation(
extent={{-20,-18},{20,22}})));
- //initial equation
- // ArtysPh.pH=7.4;
-
+//initial equation
+ // ArtysPh.pH=7.4;
BloodPhMeassure lungsPH
annotation (Placement(transformation(extent={{28,28},{0,56}})));
Tissues.Bloods tissues_blood_acidity
@@ -56140,9 +56030,9 @@
rotation=270,
origin={-20,50})));
Physiolibrary.Types.Fraction KAdjustment;
- //Real Cations( final quantity="Concentration", final unit = "mEq/l");
+ //Real Cations( final quantity="Concentration", final unit = "mEq/l");
Physiolibrary.Types.VolumeDensityOfCharge Anions;
- // Real ProteinAnions( final quantity="Concentration", final unit = "mEq/l");
+ // Real ProteinAnions( final quantity="Concentration", final unit = "mEq/l");
/* Library.Interfaces.RealInput_ HCO3(final quantity="Concentration",
final displayUnit="mEq/l") annotation (Placement(
transformation(extent={{0,30},{40,70}}), iconTransformation(
@@ -56151,11 +56041,11 @@
origin={20,50})));
*/
equation
- q_in.q + q_out.q = 0;
- //Cations = q_in.conc*1000+otherCations;
- Anions = Cations;
- // ProteinAnions = Anions - otherStrongAnions - q_in.conc*1000 - HCO3;
- KAdjustment = (Cations-(Anions-ProteinAnions))/(Cations+(Anions-ProteinAnions));
+ q_in.q + q_out.q = 0;
+//Cations = q_in.conc*1000+otherCations;
+ Anions = Cations;
+// ProteinAnions = Anions - otherStrongAnions - q_in.conc*1000 - HCO3;
+ KAdjustment = (Cations - (Anions - ProteinAnions)) / (Cations + Anions - ProteinAnions);
q_out.conc = (1+KAdjustment)*q_in.conc;
annotation (
@@ -56829,9 +56719,9 @@
rotation=270,
origin={-20,50})));
Physiolibrary.Types.Fraction KAdjustment;
- //Real Cations( final quantity="Concentration", final unit = "mEq/l");
+ //Real Cations( final quantity="Concentration", final unit = "mEq/l");
Physiolibrary.Types.VolumeDensityOfCharge Anions;
- // Real ProteinAnions( final quantity="Concentration", final unit = "mEq/l");
+ // Real ProteinAnions( final quantity="Concentration", final unit = "mEq/l");
/* Library.Interfaces.RealInput_ HCO3(final quantity="Concentration",
final displayUnit="mEq/l") annotation (Placement(
transformation(extent={{0,30},{40,70}}), iconTransformation(
@@ -56840,11 +56730,11 @@
origin={20,50})));
*/
equation
- q_in.q + q_out.q = 0;
- //Cations = q_in.conc*1000+otherCations;
- Anions = Cations;
- // ProteinAnions = Anions - otherStrongAnions - q_in.conc*1000 - HCO3;
- KAdjustment = (Cations-(Anions-ProteinAnions))/(Cations+(Anions-ProteinAnions));
+ q_in.q + q_out.q = 0;
+//Cations = q_in.conc*1000+otherCations;
+ Anions = Cations;
+// ProteinAnions = Anions - otherStrongAnions - q_in.conc*1000 - HCO3;
+ KAdjustment = (Cations - (Anions - ProteinAnions)) / (Cations + Anions - ProteinAnions);
q_out.conc = (1+KAdjustment)*q_in.conc;
annotation (
@@ -58397,10 +58287,13 @@
parameter Real GILumenVolume(displayUnit="l") = 1.3;
initial equation
- A[1] = 1.3; //Ca gut
- A[2] = 0.5; //Ca transporters
+ A[1] = 1.3;
+//Ca gut
+ A[2] = 0.5;
+//Ca transporters
A[3] = 0.839;
- A[4] = 32.9; //Ca pool
+ A[4] = 32.9;
+//Ca pool
A[5] = 16.8;
A[6] = 1260;
A[7] = 53.9;
@@ -58408,38 +58301,36 @@
A[9] = 126;
A[10] = 0.5;
A[11] = 1;
- equation
- //gut
- v_intake = D[1] * H_plus(A[2],0.9,0.75,4,0.25); //alpha was calculated == 0.35
- der( A[1]) = v_intake - v1_4;
+ equation
+//gut
+ v_intake = D[1] * H_plus(A[2], 0.9, 0.75, 4, 0.25);
+//alpha was calculated == 0.35
+ der(A[1]) = v_intake - v1_4;
v1_4 = v1_4a + v1_4p;
v1_4a = (H(A[1]/GILumenVolume,0.9,1,1) * (A[2]/0.5) / (A[1] + 0.75)) * A[1];
- v1_4p = (0.0495) * A[1];
-
- //transporters
- H62p = H_plus(A[6]/ECFV2,0.037,90,1,0.003);
+ v1_4p = (0.0495) * A[1];
+//transporters
+ H62p = H_plus(A[6] / ECFV2, 0.037, 90, 1, 0.003);
H62m = H_minus(A[6]/ECFV2,0.037,90,1,0.003);
- der( A[2]) = H62p*(1-A[2])- H62m*A[2];
-
- //Ca pool
- der(A[4]) = 0; //v1_4 - v4_u;
-
- //to urine
- v4_u = (2-H(A[6]/ECFV2,2,90,1)) * (0.3*GFR*A[4] - H(A[4]/ECFV2,(-0.38+59.22)/0.317936,1.573,1)*H(A[7]/ECFV2,1.0615,53.9/ECFV2,1));
- //!!it's used concentration of ecf Ca, not the mass such as in the article!!, H_4-u:alpha was calculated, H7,4-u:delta was assumed as normal ecf PTH concentration
- //a = (2-H(A[6]/ECFV2,2,90,1)) * (0.3*GFR*A[4]);
- //b = (2-H(A[6]/ECFV2,2,90,1)) * (H(A[4]/ECFV2,1,1.573,1)*H(A[7]/ECFV2,1.0615,53.9/ECFV2,1));
- //v4_u = v1_4;
-
- //calcitriol
- der( A[6]) = A[9] -0.1 *A[6];
- der( A[9]) = 6.3*H(A[7]/ECFV2, 1.90371, 1.549,0.1112)*( if A[5]<=16.8 then 1 else H_minus(A[5]/ECFV2,1.525,1.302,8.252,-1.78249)) -0.05 *A[9]; //H7,9:alpha was calculated,
- //a = 6.3*H(A[7]/ECFV2, 1.90371 ,1.549,0.1112);
- //c = if A[5]<16.8 then 1 else H_minus(A[5]/ECFV2,1.525,1.302,8.252,1.181);
- //b = 0.05 *A[9];
-
- //PTH
- der( A[7]) = H_minus(A[4]/ECFV2,96.25, 3.037, 11.74,6249)*(A[10]/0.5)*A[11] - 7.143*A[7];
+ der( A[2]) = H62p*(1-A[2])- H62m*A[2];
+//Ca pool
+ der(A[4]) = 0;
+//v1_4 - v4_u;
+//to urine
+ v4_u = (2 - H(A[6] / ECFV2, 2, 90, 1)) * (0.3 * GFR * A[4] - H(A[4] / ECFV2, ((-0.38) + 59.22) / 0.317936, 1.573, 1) * H(A[7] / ECFV2, 1.0615, 53.9 / ECFV2, 1));
+//!!it's used concentration of ecf Ca, not the mass such as in the article!!, H_4-u:alpha was calculated, H7,4-u:delta was assumed as normal ecf PTH concentration
+//a = (2-H(A[6]/ECFV2,2,90,1)) * (0.3*GFR*A[4]);
+//b = (2-H(A[6]/ECFV2,2,90,1)) * (H(A[4]/ECFV2,1,1.573,1)*H(A[7]/ECFV2,1.0615,53.9/ECFV2,1));
+//v4_u = v1_4;
+//calcitriol
+ der(A[6]) = A[9] - 0.1 * A[6];
+ der( A[9]) = 6.3*H(A[7]/ECFV2, 1.90371, 1.549,0.1112)*( if A[5]<=16.8 then 1 else H_minus(A[5]/ECFV2,1.525,1.302,8.252,-1.78249)) -0.05 *A[9];
+//H7,9:alpha was calculated,
+//a = 6.3*H(A[7]/ECFV2, 1.90371 ,1.549,0.1112);
+//c = if A[5]<16.8 then 1 else H_minus(A[5]/ECFV2,1.525,1.302,8.252,1.181);
+//b = 0.05 *A[9];
+//PTH
+ der(A[7]) = H_minus(A[4] / ECFV2, 96.25, 3.037, 11.74, 6249) * (A[10] / 0.5) * A[11] - 7.143 * A[7];
a= H_minus(A[4]/ECFV2,96.25, 3.037, 11.74,6249)*(A[10]/0.5)*A[11];
b=7.143*A[7];
c=b/((A[10]/0.5)*A[11]);
@@ -58448,10 +58339,9 @@
T64m = 1-(exp(0.03*(A[6]/14 - 90*(32.9/A[4]))) - exp(-0.03*(A[6]/14 - 90*(32.9/A[4])))) / (exp(0.03*(A[6]/14 - 90*(32.9/A[4]))) + exp(-0.03*(A[6]/14 - 90*(32.9/A[4]))));
T64p = 1+(exp(0.03*(A[6]/14 - 90*(32.9/A[4]))) - exp(-0.03*(A[6]/14 - 90*(32.9/A[4])))) / (exp(0.03*(A[6]/14 - 90*(32.9/A[4]))) + exp(-0.03*(A[6]/14 - 90*(32.9/A[4]))));
- der( A[11]) =0.0001604 *H_minus(A[6]/ECFV2,4.103,69,12.5,0.9) -0.0001604*A[11];
-
- //phosphate
- der(A[3])=0;
+ der( A[11]) =0.0001604 *H_minus(A[6]/ECFV2,4.103,69,12.5,0.9) -0.0001604*A[11];
+//phosphate
+ der(A[3]) = 0;
der(A[5])=0;
der(A[8])=0;
end Calcium;
@@ -61154,8 +61044,7 @@
end GangliaKidney;
model BetaReceptorsActivityFactor
- //extends Library.BaseFactor2;
-
+ //extends Library.BaseFactor2;
Modelica.Blocks.Interfaces.RealInput GangliaGeneral_NA(final quantity=
"Frequency", final displayUnit="Hz") annotation (Placement(transformation(extent={{-100,40},
{-80,60}}), iconTransformation(extent={{-108,-70},{-88,-50}})));
@@ -61261,8 +61150,7 @@
end BetaReceptorsActivityFactor;
model AplhaReceptorsActivityFactor
- //extends Library.BaseFactor2;
-
+ //extends Library.BaseFactor2;
Physiolibrary.Types.RealIO.FractionInput GangliaGeneral_NA(displayUnit="1") annotation (Placement(transformation(extent={{-100,40},
{-80,60}}), iconTransformation(extent={{-108,-70},{-88,-50}})));
Physiolibrary.Types.RealIO.FractionInput AlphaPool_Effect
@@ -62512,14 +62400,14 @@
model BladderHeat
extends Physiolibrary.Icons.Bladder;
- parameter Physiolibrary.Types.Density urineDensity = 1000; // kg/m3
+ parameter Physiolibrary.Types.Density urineDensity = 1000; // kg/m3
Physiolibrary.Thermal.Components.HeatAccumulation Bladder(
stateName="BladderTemperature.Mass",
useMassInput=true,
relativeHeat_start=0)
annotation (Placement(transformation(extent={{6,-80},{26,-60}})));
- //(259707.204) - 310.15*weight*(4186.8))
- // initialTemperature(displayUnit="degC"),
+ //(259707.204) - 310.15*weight*(4186.8))
+ // initialTemperature(displayUnit="degC"),
Physiolibrary.Types.BusConnector busConnector annotation (Placement(
transformation(extent={{-72,-44},{-52,-24}}), iconTransformation(
extent={{-104,86},{-84,106}})));
@@ -62579,9 +62467,8 @@
useMassInput=true,
relativeHeat_start=0)
annotation (Placement(transformation(extent={{-50,-82},{-30,-62}})));
- //(1293302.52) - 310.15*weight*(4186.8))
-
- // initialTemperature(displayUnit="degC") = 310,
+ //(1293302.52) - 310.15*weight*(4186.8))
+ // initialTemperature(displayUnit="degC") = 310,
Physiolibrary.Thermal.Components.Conductor resistorWithCondParam(
Conductance(displayUnit="kcal/(min.degC)") = 0.6978)
annotation (Placement(transformation(
@@ -62713,9 +62600,8 @@
useMassInput=true,
relativeHeat_start=0)
annotation (Placement(transformation(extent={{-8,-8},{12,12}})));
- //(29579490.792) - 310.15*weight*(3475.044))
-
- // initialTemperature=310.166,
+ //(29579490.792) - 310.15*weight*(3475.044))
+ // initialTemperature=310.166,
Physiolibrary.Thermal.Components.IdealRadiator muscleFlux(SpecificHeat=
3851.856, useMassFlowInput=true)
annotation (Placement(transformation(extent={{40,-8},{20,12}})));
@@ -62807,9 +62693,8 @@
useMassInput=true,
relativeHeat_start=-75362.4)
annotation (Placement(transformation(extent={{30,-72},{50,-52}})));
- //(2363163.8976) - 310.15*weight*(3475.044))
-
- //initialTemperature=300.731,
+ //(2363163.8976) - 310.15*weight*(3475.044))
+ //initialTemperature=300.731,
Physiolibrary.Thermal.Components.IdealRadiator skinFlux(
useMassFlowInput=true, SpecificHeat(displayUnit="kcal/(kg.K)") = 3851.856)
annotation (Placement(transformation(extent={{38,-40},{58,-20}})));
@@ -63070,7 +62955,7 @@
model Heat2
extends Physiolibrary.Icons.Heat;
- // initialTemperature(displayUnit="degC") = 310.15,
+ // initialTemperature(displayUnit="degC") = 310.15,
Physiolibrary.Types.BusConnector busConnector annotation (Placement(
transformation(extent={{-104,86},{-84,106}}), iconTransformation(
extent={{-104,86},{-84,106}})));
@@ -64141,8 +64026,7 @@
useMassInput=true,
relativeHeat_start=0)
annotation (Placement(transformation(extent={{-6,-54},{14,-34}})));
- //(45005587.92) - 310.15*weight*(3475.044))
-
+ //(45005587.92) - 310.15*weight*(3475.044))
Physiolibrary.Types.RealIO.TemperatureOutput
core_T
annotation (Placement(
@@ -64611,7 +64495,7 @@
model Heat3
extends Physiolibrary.Icons.Heat;
- // initialTemperature(displayUnit="degC") = 310.15,
+ // initialTemperature(displayUnit="degC") = 310.15,
Physiolibrary.Types.BusConnector busConnector annotation (Placement(
transformation(extent={{-104,86},{-84,106}}), iconTransformation(
extent={{-104,86},{-84,106}})));
@@ -69001,8 +68885,8 @@
Boolean a(start=false);
Boolean b(start=true);
Boolean c(start=false);
- //initial equation
- // c=false;
+ //initial equation
+ // c=false;
equation
a=FunctionEffecteffectAtRevitaling;
@@ -69381,37 +69265,37 @@
end PatientStatus;
block Normal
- extends Modelica.StateGraph.InitialStepWithSignal;
+ extends Modelica.StateGraph.InitialStepWithSignal(nIn = 1, nOut = 1);
extends Physiolibrary.Icons.StatusNormal;
end Normal;
block Confused
- extends Modelica.StateGraph.StepWithSignal;
+ extends Modelica.StateGraph.StepWithSignal(nIn = 1, nOut = 1);
extends Physiolibrary.Icons.StatusConfused;
end Confused;
block Impaired
- extends Modelica.StateGraph.StepWithSignal;
+ extends Modelica.StateGraph.StepWithSignal(nIn = 1, nOut = 1);
extends Physiolibrary.Icons.StatusImpaired;
end Impaired;
block Comatose
- extends Modelica.StateGraph.StepWithSignal;
+ extends Modelica.StateGraph.StepWithSignal(nIn = 1, nOut = 1);
extends Physiolibrary.Icons.StatusComatose;
end Comatose;
block NotBreathing
- extends Modelica.StateGraph.StepWithSignal;
+ extends Modelica.StateGraph.StepWithSignal(nIn = 1, nOut = 1);
extends Physiolibrary.Icons.StatusNotBreathing;
end NotBreathing;
block MayBeDead
- extends Modelica.StateGraph.StepWithSignal;
+ extends Modelica.StateGraph.StepWithSignal(nIn = 1, nOut = 1);
extends Physiolibrary.Icons.StatusMayBeDead;
end MayBeDead;
block IsReallyDead
- extends Modelica.StateGraph.StepWithSignal;
+ extends Modelica.StateGraph.StepWithSignal(nIn = 1, nOut = 1);
extends Physiolibrary.Icons.StatusIsReallyDead;
end IsReallyDead;
@@ -69663,6 +69547,7 @@
end RightHeart;
end Tissues;
+
package IO_Bus
extends Physiolibrary.Types.IO_Bus;
redeclare model extends Variables
@@ -69691,9 +69576,9 @@
T.Fraction KidneyFunctionEffect(varName="Kidney-Function.Effect")
"Kidney tissue function effect. []"
annotation (Placement(transformation(extent={{-36,80},{-30,86}})));
- //T.Fraction KidneyFunction_Effect(varName="Kidney-Function.Effect")
- // "Kidney tissue function effect. []"
- //annotation (Placement(transformation(extent={{-36,90},{-30,96}})));
+ //T.Fraction KidneyFunction_Effect(varName="Kidney-Function.Effect")
+ // "Kidney tissue function effect. []"
+ //annotation (Placement(transformation(extent={{-36,90},{-30,96}})));
T.Fraction Kidney_StructureEffect(varName="Kidney-Structure.Effect")
"Kidney tissue structural damage effect. []"
annotation (Placement(transformation(extent={{-36,100},{-30,106}})));
@@ -70195,7 +70080,6 @@
smooth=Smooth.None));
end Status_test_SI;
end IO_Bus;
-
annotation (Documentation(revisions="
Author: |
@@ -70220,7 +70104,7 @@
end Status;
- annotation (uses(Modelica(version="3.2.1"), Physiolibrary(version="2.3.1")), Documentation(revisions="
+ annotation (uses(Modelica(version = "4.0.0")), Documentation(revisions="