Location: 12 L Platform 1 model codes @ a8a92308e217 / Airway_FTU / Airway_FTU / airway1_radius1.cellml
- Author:
- aram148 <42922407+aram148@users.noreply.github.com>
- Date:
- 2022-07-22 15:47:05+12:00
- Desc:
- Added documentation for VSM model
- Permanent Source URI:
- https://models.cellml.org/workspace/6b0/rawfile/a8a92308e217ac5626809237dd90a31240b22834/Airway_FTU/Airway_FTU/airway1_radius1.cellml
<?xml version='1.0'?>
<model name="airway2_radius" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:xlink="http://www.w3.org/1999/xlink">
<units name="mm2">
<unit exponent="2" prefix="milli" units="metre"/>
</units>
<units name="mm">
<unit prefix="milli" units="metre"/>
</units>
<units name="per_sec">
<unit exponent="-1" units="second"/>
</units>
<component name="Environment">
<variable cmeta:id="Environment_time" name="time" public_interface="out" units="second"/>
</component>
<import xlink:href="hai_1.cellml">
<component component_ref="stress" name="stress"/>
<component component_ref="AMp" name="AMp"/>
<component component_ref="AM" name="AM"/>
<component component_ref="Mp" name="Mp"/>
<component component_ref="M" name="M"/>
<component component_ref="model_parameters1" name="model_parameters1"/>
</import>
<component name="lambert_r21">
<variable name="time" public_interface="in" units="second"/>
<variable name="Ri_sq_r21" public_interface="in" units="mm2"/>
<variable name="rmax_sq_r21" public_interface="in" units="mm"/>
<variable name="N1_r21" public_interface="in" units="dimensionless"/>
<variable name="N2_r21" public_interface="in" units="dimensionless"/>
<variable name="P1_r21" public_interface="in" units="pascal"/>
<variable name="P2_r21" public_interface="in" units="pascal"/>
<variable name="rho" public_interface="in" units="dimensionless"/>
<variable initial_value="0.2960" name="r21" public_interface="out" units="mm"/>
<variable name="Ptm_r21" units="pascal"/>
<variable name="P0" public_interface="in" units="pascal"/>
<!-- var kappa: dimensionless {pub: in};-->
<variable name="rad21" public_interface="out" units="mm"/>
<variable name="rmax_r21" public_interface="in" units="mm"/>
<variable name="V_r21" units="dimensionless"/>
<variable name="gamma" public_interface="in" units="dimensionless"/>
<variable name="stress" public_interface="in" units="dimensionless"/>
<variable name="v_smc_r21" public_interface="in" units="dimensionless"/>
<variable name="f_a_r21" public_interface="in" units="pascal"/>
<variable name="Pmin_r21" public_interface="out" units="pascal"/>
<variable name="P_t_r21" public_interface="in" units="pascal"/>
<variable name="tau_r21" public_interface="in" units="dimensionless"/>
<variable name="pmidr21" public_interface="in" units="pascal"/>
<variable initial_value="0.445" name="R_ref" units="mm"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Pmin_r21</ci>
<apply>
<plus/>
<ci>P0</ci>
<ci>P_t_r21</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Ptm_r21</ci>
<apply>
<plus/>
<apply>
<minus/>
<ci>pmidr21</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>f_a_r21</ci>
<ci>R_ref</ci>
</apply>
<ci>r21</ci>
</apply>
</apply>
<ci>tau_r21</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>rad21</ci>
<piecewise>
<piece>
<apply>
<root/>
<apply>
<times/>
<ci>Ri_sq_r21</ci>
<apply>
<power/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Ptm_r21</ci>
<ci>P1_r21</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>N1_r21</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<leq/>
<ci>Ptm_r21</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<piece>
<apply>
<root/>
<apply>
<minus/>
<ci>rmax_sq_r21</ci>
<apply>
<times/>
<apply>
<minus/>
<ci>rmax_sq_r21</ci>
<ci>Ri_sq_r21</ci>
</apply>
<apply>
<power/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Ptm_r21</ci>
<ci>P2_r21</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>N2_r21</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<geq/>
<ci>Ptm_r21</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>r21</ci>
</apply>
<apply>
<times/>
<ci>rho</ci>
<apply>
<minus/>
<ci>rad21</ci>
<ci>r21</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>V_r21</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<ci>gamma</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<pi/>
<ci>rmax_r21</ci>
</apply>
</apply>
<ci>rho</ci>
<apply>
<minus/>
<ci>rad21</ci>
<ci>r21</ci>
</apply>
<ci>v_smc_r21</ci>
</apply>
</apply>
</math>
</component>
<component name="lambert_r22">
<variable name="time" public_interface="in" units="second"/>
<variable name="Ri_sq_r22" public_interface="in" units="mm2"/>
<variable name="rmax_sq_r22" public_interface="in" units="mm"/>
<variable name="N1_r22" public_interface="in" units="dimensionless"/>
<variable name="N2_r22" public_interface="in" units="dimensionless"/>
<variable name="P1_r22" public_interface="in" units="pascal"/>
<variable name="P2_r22" public_interface="in" units="pascal"/>
<variable name="rho" public_interface="in" units="dimensionless"/>
<variable initial_value="0.28" name="r22" public_interface="out" units="mm"/>
<variable name="Ptm_r22" units="pascal"/>
<variable name="P0" public_interface="in" units="pascal"/>
<!-- var kappa: dimensionless {pub: in};-->
<variable name="rad22" public_interface="out" units="mm"/>
<variable name="rmax_r22" public_interface="in" units="mm"/>
<variable name="V_r22" units="dimensionless"/>
<variable name="gamma" public_interface="in" units="dimensionless"/>
<variable name="stress" public_interface="in" units="dimensionless"/>
<variable name="v_smc_r22" public_interface="in" units="dimensionless"/>
<variable name="f_a_r22" public_interface="in" units="pascal"/>
<variable name="Pmin_r22" public_interface="out" units="pascal"/>
<variable name="P_t_r22" public_interface="in" units="pascal"/>
<variable name="tau_r22" public_interface="in" units="dimensionless"/>
<variable name="pmidr22" public_interface="in" units="pascal"/>
<variable initial_value="0.445" name="R_ref" units="mm"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Pmin_r22</ci>
<apply>
<plus/>
<ci>P0</ci>
<ci>P_t_r22</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Ptm_r22</ci>
<apply>
<plus/>
<apply>
<minus/>
<ci>pmidr22</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>f_a_r22</ci>
<ci>R_ref</ci>
</apply>
<ci>r22</ci>
</apply>
</apply>
<ci>tau_r22</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>rad22</ci>
<piecewise>
<piece>
<apply>
<root/>
<apply>
<times/>
<ci>Ri_sq_r22</ci>
<apply>
<power/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Ptm_r22</ci>
<ci>P1_r22</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>N1_r22</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<leq/>
<ci>Ptm_r22</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<piece>
<apply>
<root/>
<apply>
<minus/>
<ci>rmax_sq_r22</ci>
<apply>
<times/>
<apply>
<minus/>
<ci>rmax_sq_r22</ci>
<ci>Ri_sq_r22</ci>
</apply>
<apply>
<power/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Ptm_r22</ci>
<ci>P2_r22</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>N2_r22</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<geq/>
<ci>Ptm_r22</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>r22</ci>
</apply>
<apply>
<times/>
<ci>rho</ci>
<apply>
<minus/>
<ci>rad22</ci>
<ci>r22</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>V_r22</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<ci>gamma</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<pi/>
<ci>rmax_r22</ci>
</apply>
</apply>
<ci>rho</ci>
<apply>
<minus/>
<ci>rad22</ci>
<ci>r22</ci>
</apply>
<ci>v_smc_r22</ci>
</apply>
</apply>
</math>
</component>
<component name="airway2_vals">
<variable name="e_w_r21" units="dimensionless"/>
<variable name="e_m_r21" units="dimensionless"/>
<variable name="r_w_r21" units="mm"/>
<variable name="r_m_r21" units="mm"/>
<variable name="r_s_r21" units="mm"/>
<variable name="Rw_r21" public_interface="in" units="mm"/>
<variable name="Ri_r21" public_interface="in" units="mm"/>
<variable name="Rm_r21" public_interface="in" units="mm"/>
<variable name="fl_r21" units="dimensionless"/>
<variable name="rsmax_r21" public_interface="in" units="mm"/>
<variable name="f_a_r21" public_interface="out" units="pascal"/>
<variable name="stress" public_interface="in" units="dimensionless"/>
<variable name="D_Rm_r21" units="dimensionless"/>
<variable name="sigma_rrm_r21" units="dimensionless"/>
<variable name="sigma_rrw_r21" units="dimensionless"/>
<variable name="sigma_ttm_r21" units="dimensionless"/>
<variable name="sigma_ttw_r21" units="dimensionless"/>
<variable name="moo" public_interface="in" units="dimensionless"/>
<variable name="noo" public_interface="in" units="dimensionless"/>
<variable name="P_t_r21" public_interface="out" units="pascal"/>
<variable name="f" public_interface="in" units="per_sec"/>
<variable name="time" public_interface="in" units="second"/>
<variable name="rad21" public_interface="in" units="mm"/>
<variable name="r21" public_interface="in" units="mm"/>
<variable name="v_smc_r21" public_interface="out" units="dimensionless"/>
<variable name="rho" public_interface="in" units="dimensionless"/>
<variable name="kappa" public_interface="in" units="pascal"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>P_t_r21</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0</cn>
<apply>
<sin/>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<pi/>
<ci>f</ci>
<ci>time</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>e_w_r21</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>Rw_r21</ci>
<ci>Ri_r21</ci>
</apply>
<ci>Ri_r21</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>e_m_r21</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>Rm_r21</ci>
<ci>Rw_r21</ci>
</apply>
<ci>Ri_r21</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>r_w_r21</ci>
<apply>
<times/>
<ci>Ri_r21</ci>
<apply>
<root/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>e_w_r21</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<apply>
<divide/>
<ci>rad21</ci>
<ci>Ri_r21</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>r_m_r21</ci>
<apply>
<times/>
<ci>Ri_r21</ci>
<apply>
<root/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<power/>
<apply>
<plus/>
<ci>e_m_r21</ci>
<ci>e_w_r21</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<apply>
<divide/>
<ci>rad21</ci>
<ci>Ri_r21</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>r_s_r21</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<apply>
<plus/>
<ci>r_w_r21</ci>
<ci>r_m_r21</ci>
</apply>
</apply>
</apply>
<!-- Expt data fit for AM force at latch state-->
<!-- fl_r21 = sel
case r_s_r21 <= 2{dimensionless}*rsmax_r21:
pow(sin(pi*r_s_r21/(2{dimensionless}*rsmax_r21)), 3{dimensionless});
otherwise:
0{dimensionless};
endsel;-->
<apply>
<eq/>
<ci>fl_r21</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<!-- Total active force-->
<apply>
<eq/>
<ci>f_a_r21</ci>
<apply>
<times/>
<ci>fl_r21</ci>
<ci>kappa</ci>
<ci>stress</ci>
</apply>
</apply>
<!-- Radial stress-->
<apply>
<eq/>
<ci>D_Rm_r21</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>Rm_r21</ci>
<ci>r_m_r21</ci>
</apply>
<ci>Rm_r21</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>sigma_rrm_r21</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>moo</ci>
<apply>
<plus/>
<ci>D_Rm_r21</ci>
<apply>
<power/>
<apply>
<times/>
<ci>noo</ci>
<ci>D_Rm_r21</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<ci>P_t_r21</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>sigma_rrw_r21</ci>
<apply>
<minus/>
<ci>sigma_rrm_r21</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>f_a_r21</ci>
<apply>
<minus/>
<ci>Rm_r21</ci>
<ci>Rw_r21</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<apply>
<plus/>
<ci>r_m_r21</ci>
<ci>r_w_r21</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<!-- Hoop stress-->
<apply>
<eq/>
<ci>sigma_ttm_r21</ci>
<apply>
<plus/>
<ci>sigma_rrm_r21</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>f_a_r21</ci>
<apply>
<minus/>
<ci>Rm_r21</ci>
<ci>Rw_r21</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>r_m_r21</ci>
<ci>r_w_r21</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>sigma_ttw_r21</ci>
<apply>
<plus/>
<ci>sigma_rrw_r21</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>f_a_r21</ci>
<apply>
<minus/>
<ci>Rm_r21</ci>
<ci>Rw_r21</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>r_m_r21</ci>
<ci>r_w_r21</ci>
</apply>
</apply>
</apply>
</apply>
<!--Vel of SMC contraction-->
<apply>
<eq/>
<ci>v_smc_r21</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<minus/>
<ci>r_w_r21</ci>
<ci>r_m_r21</ci>
</apply>
</apply>
<ci>rho</ci>
<apply>
<minus/>
<ci>rad21</ci>
<ci>r21</ci>
</apply>
<ci>rad21</ci>
<apply>
<ln/>
<apply>
<divide/>
<ci>r_w_r21</ci>
<ci>r_m_r21</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="airway22_vals">
<variable name="e_w_r22" units="dimensionless"/>
<variable name="e_m_r22" units="dimensionless"/>
<variable name="r_w_r22" units="mm"/>
<variable name="r_m_r22" units="mm"/>
<variable name="r_s_r22" units="mm"/>
<variable name="Rw_r22" public_interface="in" units="mm"/>
<variable name="Ri_r22" public_interface="in" units="mm"/>
<variable name="Rm_r22" public_interface="in" units="mm"/>
<variable name="fl_r22" units="dimensionless"/>
<variable name="rsmax_r22" public_interface="in" units="mm"/>
<variable name="f_a_r22" public_interface="out" units="pascal"/>
<variable name="stress" public_interface="in" units="dimensionless"/>
<variable name="D_Rm_r22" units="dimensionless"/>
<variable name="sigma_rrm_r22" units="dimensionless"/>
<variable name="sigma_rrw_r22" units="dimensionless"/>
<variable name="sigma_ttm_r22" units="dimensionless"/>
<variable name="sigma_ttw_r22" units="dimensionless"/>
<variable name="moo" public_interface="in" units="dimensionless"/>
<variable name="noo" public_interface="in" units="dimensionless"/>
<variable name="P_t_r22" public_interface="out" units="pascal"/>
<variable name="f" public_interface="in" units="per_sec"/>
<variable name="time" public_interface="in" units="second"/>
<variable name="rad22" public_interface="in" units="mm"/>
<variable name="r22" public_interface="in" units="mm"/>
<variable name="v_smc_r22" public_interface="out" units="dimensionless"/>
<variable name="rho" public_interface="in" units="dimensionless"/>
<variable name="kappa" public_interface="in" units="pascal"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>P_t_r22</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0</cn>
<apply>
<sin/>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<pi/>
<ci>f</ci>
<ci>time</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>e_w_r22</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>Rw_r22</ci>
<ci>Ri_r22</ci>
</apply>
<ci>Ri_r22</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>e_m_r22</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>Rm_r22</ci>
<ci>Rw_r22</ci>
</apply>
<ci>Ri_r22</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>r_w_r22</ci>
<apply>
<times/>
<ci>Ri_r22</ci>
<apply>
<root/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>e_w_r22</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<apply>
<divide/>
<ci>rad22</ci>
<ci>Ri_r22</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>r_m_r22</ci>
<apply>
<times/>
<ci>Ri_r22</ci>
<apply>
<root/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<power/>
<apply>
<plus/>
<ci>e_m_r22</ci>
<ci>e_w_r22</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<apply>
<divide/>
<ci>rad22</ci>
<ci>Ri_r22</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>r_s_r22</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<apply>
<plus/>
<ci>r_w_r22</ci>
<ci>r_m_r22</ci>
</apply>
</apply>
</apply>
<!-- Expt data fit for AM force at latch state-->
<!-- fl_r22 = sel
case r_s_r22 <= 2{dimensionless}*rsmax_r22:
pow(sin(pi*r_s_r22/(2{dimensionless}*rsmax_r22)), 3{dimensionless});
otherwise:
0{dimensionless};
endsel;-->
<apply>
<eq/>
<ci>fl_r22</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<!-- Total active force-->
<apply>
<eq/>
<ci>f_a_r22</ci>
<apply>
<times/>
<ci>fl_r22</ci>
<ci>kappa</ci>
<ci>stress</ci>
</apply>
</apply>
<!-- Radial stress-->
<apply>
<eq/>
<ci>D_Rm_r22</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>Rm_r22</ci>
<ci>r_m_r22</ci>
</apply>
<ci>Rm_r22</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>sigma_rrm_r22</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>moo</ci>
<apply>
<plus/>
<ci>D_Rm_r22</ci>
<apply>
<times/>
<ci>noo</ci>
<apply>
<power/>
<ci>D_Rm_r22</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<ci>P_t_r22</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>sigma_rrw_r22</ci>
<apply>
<minus/>
<ci>sigma_rrm_r22</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>f_a_r22</ci>
<apply>
<minus/>
<ci>Rm_r22</ci>
<ci>Rw_r22</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<apply>
<plus/>
<ci>r_m_r22</ci>
<ci>r_w_r22</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<!-- Hoop stress-->
<apply>
<eq/>
<ci>sigma_ttm_r22</ci>
<apply>
<plus/>
<ci>sigma_rrm_r22</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>f_a_r22</ci>
<apply>
<minus/>
<ci>Rm_r22</ci>
<ci>Rw_r22</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>r_m_r22</ci>
<ci>r_w_r22</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>sigma_ttw_r22</ci>
<apply>
<plus/>
<ci>sigma_rrw_r22</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>f_a_r22</ci>
<apply>
<minus/>
<ci>Rm_r22</ci>
<ci>Rw_r22</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>r_m_r22</ci>
<ci>r_w_r22</ci>
</apply>
</apply>
</apply>
</apply>
<!--Vel of SMC contraction-->
<apply>
<eq/>
<ci>v_smc_r22</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<minus/>
<ci>r_w_r22</ci>
<ci>r_m_r22</ci>
</apply>
</apply>
<ci>rho</ci>
<apply>
<minus/>
<ci>rad22</ci>
<ci>r22</ci>
</apply>
<ci>rad22</ci>
<apply>
<ln/>
<apply>
<divide/>
<ci>r_w_r22</ci>
<ci>r_m_r22</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="parameters_r2">
<variable initial_value="0.003364" name="Ri_sq_r21" public_interface="out" units="mm2"/>
<variable name="Ri_sq_r22" public_interface="out" units="mm2"/>
<variable initial_value="0.0876" name="rmax_sq_r21" public_interface="out" units="mm"/>
<variable name="rmax_sq_r22" public_interface="out" units="mm"/>
<variable initial_value="1" name="N1_r21" public_interface="out" units="dimensionless"/>
<variable name="N1_r22" public_interface="out" units="dimensionless"/>
<variable initial_value="7" name="N2_r21" public_interface="out" units="dimensionless"/>
<variable name="N2_r22" public_interface="out" units="dimensionless"/>
<variable initial_value="15.728" name="P1_r21" public_interface="out" units="pascal"/>
<variable name="P1_r22" public_interface="out" units="pascal"/>
<variable initial_value="-2757.374611" name="P2_r21" public_interface="out" units="pascal"/>
<variable name="P2_r22" public_interface="out" units="pascal"/>
<variable initial_value="1" name="rho" public_interface="out" units="dimensionless"/>
<variable initial_value="1333.2239" name="P0" public_interface="out" units="pascal"/>
<!-- var kappa: dimensionless {pub:in};-->
<variable initial_value="0.296" name="rmax_r21" public_interface="out" units="mm"/>
<variable name="rmax_r22" public_interface="out" units="mm"/>
<variable initial_value="25" name="gamma" public_interface="out" units="dimensionless"/>
<variable initial_value="0.109" name="Rw_r21" public_interface="out" units="mm"/>
<variable name="Rw_r22" public_interface="out" units="mm"/>
<variable initial_value="0.058" name="Ri_r21" public_interface="out" units="mm"/>
<variable name="Ri_r22" public_interface="out" units="mm"/>
<variable initial_value="0.121" name="Rm_r21" public_interface="out" units="mm"/>
<variable name="Rm_r22" public_interface="out" units="mm"/>
<variable initial_value="0.312" name="rsmax_r21" public_interface="out" units="mm"/>
<variable name="rsmax_r22" public_interface="out" units="mm"/>
<variable initial_value="5.8" name="moo" public_interface="out" units="dimensionless"/>
<variable initial_value="1.5" name="noo" public_interface="out" units="dimensionless"/>
<variable initial_value="0.25" name="f" public_interface="out" units="per_sec"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Ri_sq_r22</ci>
<ci>Ri_sq_r21</ci>
</apply>
<apply>
<eq/>
<ci>rmax_sq_r22</ci>
<ci>rmax_sq_r21</ci>
</apply>
<apply>
<eq/>
<ci>N1_r22</ci>
<ci>N1_r21</ci>
</apply>
<apply>
<eq/>
<ci>N2_r22</ci>
<ci>N2_r21</ci>
</apply>
<apply>
<eq/>
<ci>P1_r22</ci>
<ci>P1_r21</ci>
</apply>
<apply>
<eq/>
<ci>P2_r22</ci>
<ci>P2_r21</ci>
</apply>
<apply>
<eq/>
<ci>rmax_r22</ci>
<ci>rmax_r21</ci>
</apply>
<apply>
<eq/>
<ci>Rw_r22</ci>
<ci>Rw_r21</ci>
</apply>
<apply>
<eq/>
<ci>Ri_r22</ci>
<ci>Ri_r21</ci>
</apply>
<apply>
<eq/>
<ci>Rm_r22</ci>
<ci>Rm_r21</ci>
</apply>
<apply>
<eq/>
<ci>rsmax_r22</ci>
<ci>rsmax_r21</ci>
</apply>
</math>
</component>
</model>