Innovation | Open | Published:
An empirical model for educational simulation of cervical dilation in first-stage labor
Advances in Simulationvolume 3, Article number: 9 (2018)
Several models for educational simulation of labor and delivery were published in the literature and incorporated into a commercially available training simulator (CAE Healthcare Lucina). However, the engine of this simulator does not include a model for the clinically relevant indicators: uterine contraction amplitude and frequency, and cervical dilation. In this paper, such a model is presented for the primigravida in normal labor.
The conceptual and mathematical models represent oxytocin release by the hypothalamus, oxytocin pharmacokinetics, and oxytocin effect on uterine contractions, cervical dilation, and (positive) feedback from cervical dilation to oxytocin release by the hypothalamus.
Simulation results for cervical dilation are presented, together with target data for a normal primigravida. Corresponding oxytocin concentrations and amplitude and frequency of uterine contractions are also presented.
An original empirical model for educational simulation of oxytocin concentration, uterine contractions, and cervical dilation in first-stage labor is presented. Simulation results for cervical dilation match target data for a normal patient. The model forms a basis for taking into account more independent variables and patient profiles and can thereby considerably expand the range of training scenarios that can be simulated.
Uterine contractions, fetal descent, and cervical dilation are used to assess progression of labor. These related phenomena vary greatly among parturients and may have an incidence on the condition of the fetus. Simulators allow for practice of normal and critical situations at will and at no risk to real patients. Several models for educational simulation of labor and delivery were published in the literature [1,2,3] and incorporated into a commercially available training simulator . The cited models make amplitude, frequency, duration, and resting tone of the uterine pressure waveform evolve spontaneously and under the influence of oxytocin and tocolytics. However, the relationships between these variables and cervical dilation are scripted.
A positive feedback loop involving oxytocin plays a central role in the progression of first-stage labor. The hypothalamus releases oxytocin into the bloodstream, which is transported to the uterus and stimulates uterine contractions. The contractions result in descent of the fetus. The descending fetal head causes cervical dilation. Cervical stretch receptors send a signal to the hypothalamus via afferent nerves, stimulating additional release of oxytocin, thereby closing the loop. This loop results in progressively intensifying contractions and dilation over time .
A database for a number of contraction and dilation profiles could be established and used for scripted training simulation purposes. However, realistic simulation of the interaction between oxytocin concentrations, uterine contractions, fetal descent, and dilation, for a number of patient profiles, and possibly under external influences such as administration of exogenous oxytocin and tocolytics, requires a model-based approach. Such a model would not only increase insight in the interacting processes, but also greatly expand the range of scenarios that can be simulated beyond the relatively simple phenomena that can be characterized in a database.
In this innovation paper, we lay the foundation for such a model by proposing an empirical model that simulates oxytocin concentrations and contraction and dilation profiles. The model should be able to run in real, or accelerated, time and should be easy to manipulate to reflect different patients.
Basic assumptions are made, similar to the ones in the article of Bastos et al.  and based on understanding of physiological relations, in order to keep the model as simple as possible, while approaching the target data of Zhang et al. . The subsystems and variables connecting them, as described in the introduction, are reflected in the block diagram of Fig. 1.
The model equations referring to Fig. 1 from left to right are introduced below. The rate of oxytocin release by the hypothalamus r(t) in mU/min has a fixed component, reflected in the parameter P 1 and depends linearly via a gain P 2 on the sensed dilation d(t) in centimeters:
Units of parameters can easily be derived from the units of the variables. The change in mass of oxytocin depends on the elimination rate constant P 3 , the current mass m(t) in mU, and the release rate, represented by a first-order pharmacokinetic model:
Dividing the mass by the volume of distribution P 4 results in the concentration of oxytocin in mU/mL:
The concentration dependency of contraction frequency f(t) in 1/min and contraction amplitude a(t) in mm Hg is represented by a sigmoidal curve, Eqs. (4) and (5). For example in Eq. (4), for a concentration equal to zero, the frequency is zero; for a concentration equal to P 6 , it is 50% of the maximum frequency P 5 ; and for higher concentrations, it tends toward P 5 . The slope of both pharmacodynamic relationships is governed by P 7 :
The amplitude has a non-zero baseline P 8 and a concentration-dependent part reflected in P 9 (the difference between the maximum and baseline amplitude values):
The change in dilation depends on the pressure exerted by the fetus on the cervix P 10 , which is assumed fixed in this version of the model, and on the product of frequency and amplitude. This second influence is assumed linearly proportional with gain P 11 . Thus, the resulting uterine model is
Referenced and empirically adjusted parameters are listed in Table 1.
Equations (1–6) were implemented in MATLAB (MATLAB R2017a, The MathWorks Inc.) using the Euler forward method with an integration step size of 1 min to integrate the two differential equations. Model simulation results are presented and compared to target data obtained from Zhang et al.  for cervical dilation.
Figure 2 shows simulation results from the initial conditions d(0) = 2.00 cm, and m(0) = 275 mU, using the parameters listed in Table 1. The stars in Fig. 1 correspond to dilation data for patient A in Zhang et al. . This patient was admitted at 2 cm and labor progressed to 10 cm with labor considered non-protracted.
Simulated dilation matches dilation in a normal patient in an approximation that is considered realistic enough for educational simulation. Parameters P 2 and P 11 were adjusted to achieve this match. Further research could address which parameters are patient- or population-specific and can be used to simulate different labor types. The model parameter estimation procedure could then be formalized so that a clinical instructor can easily use it to meet his or her educational objectives. Ongoing work involves expanding the set of target data and matching simulation results to them, to enhance model validity and demonstrate the possibility to match different patient profiles. Critical analysis of the pharmacokinetic model may be necessary. It would also be interesting to explore validity of the model for exogenous oxytocin administration.
An original empirical model for educational simulation of oxytocin concentrations, uterine contractions, and cervical dilation in first-stage labor is presented. Simulation results for cervical dilation match target data for a normal patient. The proposed model forms a sound, explicit basis for taking into account more independent variables and patient profiles, and thereby considerably expand the range of training scenarios that can be simulated.
Bastos LF, Lobo MF, van Meurs WL, Ayres-de-Campos D. An intrauterine pressure generator for educational simulation of labour and delivery. Med Eng Phys. 2010;32:740–5.
Bastos LF, van Meurs W, Ayres-de-Campos D. A model for educational simulation of the evolution of uterine contractions during labor. Comput Methods Prog Biomed. 2012;107:242–7.
Lobo MF, Bastos LF, van Meurs WL, Ayres-de-Campos D. A model for educational simulation of the effect of oxytocin on uterine contractions. Med Eng Phys. 2013;35:524–31.
Childbirth simulator. https://caehealthcare.com/patient-simulation/lucina. Accessed 8 June 2017.
Marieb EN, Hoehn K. Pregnancy and Human Development. In: Human anatomy and physiology. 10th ed: Pearson Education, Limited; 2015. p. 1116–7.
Zhang J, Landy H, Branch D, Burkman R. Contemporary patterns of spontaneous labor with normal neonatal outcomes. Obstet Gynecol. 2013;116:1281–7.
Fuchs A-R, Romero R, Keefe D, Parra M, Oyarzun E, Behnke E. Oxytocin secretion and human parturition: pulse frequency and duration increase during spontaneous labor in women. Am J Obstet Gynecol. 1991;165:1515–23.
Becker KL. Physiology of vasopressin, oxytocin, and thirst. In: Principles and practice of endocrinology and metabolism. 3rd ed: Lippincott Williams and Wilkins; 2001. p. 276–84.
Brandis K. Fluid Physiology. http://www.anaesthesiamcq.com/FluidBook/fl2_1.php. Accessed 1 June 2017.
Caldeyro-Barcia R, Poseiro JJ. Oxytocin and contractility of the pregnant human uterus. Ann N Y Acad Sci. 1959;75(2 The Uterus):813–30.
Antonucci MC, Pitman MC, Eid T, Steer PJ, Genevier ES. Simultaneous monitoring of head-to-cervix forces, intrauterine pressure and cervical dilatation during labour. Med Eng Phys. 1997;19:317–26.
Silvano Gefferie won the Nicolaas Westerhof first prize for a young investigator for his presentation of the work elaborated upon in this paper at the conference Applied Modeling in Acute Care, April 20–21, 2017, Enschede, The Netherlands. Part of the salary of Willem van Meurs in the 2016–2017 academic year was funded by a grant from CAE Healthcare, Montreal, Canada.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Drs. Ferreira Bastos and Van Meurs received royalties for licensing of the Childbirth simulator mentioned in the introduction by the University of Porto to CAE Healthcare, Montreal, Canada. In the 2016–2017 academic year, Willem van Meurs was a consultant to and received a grant from CAE Healthcare, Montreal, Canada. The remaining authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.