
Investigation of the dynamic properties of implants is essential to ensure safety and compatibility with the host's natural spinal tissue. This paper presents a simplified model of a cantilever beam to investigate the effects of holes/pores on the structures. Free vibration test is one of the most effective methods to measure the dynamic response of a cantilever beam, such as natural frequency and damping ratio. In this study, the natural frequencies of cantilever beams made of polycarbonate (PC) containing various circular open holes were investigated numerically, analytically, and experimentally. The experimental data confirmed the accuracy of the natural frequencies of the cantilever beam with open holes calculated by finite element and analytical models. In addition, two finite element simulation methods, the dynamic explicit and modal dynamic methods, were applied to determine the damping ratios of cantilever beams with open holes. Finite element analysis accurately simulated the damped vibration behavior of cantilever beams with open holes when known material damping properties were applied. The damping behavior of cantilever beams with random pores was simulated, highlighting a completely different relationship between porosity, natural frequency and damping response. The latter highlights the potential of finite element methods to analyze the dynamic response of arbitrary and complex structures, towards improved implant design.
Citation: Xiaoyu Du, Yijun Zhou, Lingzhen Li, Cecilia Persson, Stephen J. Ferguson. The porous cantilever beam as a model for spinal implants: Experimental, analytical and finite element analysis of dynamic properties[J]. Mathematical Biosciences and Engineering, 2023, 20(4): 6273-6293. doi: 10.3934/mbe.2023270
[1] | Tianyong Hao, Zhengxing Huang, Likeng Liang, Buzhou Tang . Biomedical and health information processing and analysis. Mathematical Biosciences and Engineering, 2020, 17(5): 5147-5149. doi: 10.3934/mbe.2020278 |
[2] | Baohua Hu, Yong Wang, Jingsong Mu . A new fractional fuzzy dispersion entropy and its application in muscle fatigue detection. Mathematical Biosciences and Engineering, 2024, 21(1): 144-169. doi: 10.3934/mbe.2024007 |
[3] | Shaojun Zhu, Jinhui Zhao, Yating Wu, Qingshan She . Intermuscular coupling network analysis of upper limbs based on R-vine copula transfer entropy. Mathematical Biosciences and Engineering, 2022, 19(9): 9437-9456. doi: 10.3934/mbe.2022439 |
[4] | Ying-Cheng Lai, Kwangho Park . Noise-sensitive measure for stochastic resonance in biological oscillators. Mathematical Biosciences and Engineering, 2006, 3(4): 583-602. doi: 10.3934/mbe.2006.3.583 |
[5] | Bei Liu, Hongzi Bai, Wei Chen, Huaquan Chen, Zhen Zhang . Automatic detection method of epileptic seizures based on IRCMDE and PSO-SVM. Mathematical Biosciences and Engineering, 2023, 20(5): 9349-9363. doi: 10.3934/mbe.2023410 |
[6] | Erik M. Bollt, Joseph D. Skufca, Stephen J . McGregor . Control entropy: A complexity measure for nonstationary signals. Mathematical Biosciences and Engineering, 2009, 6(1): 1-25. doi: 10.3934/mbe.2009.6.1 |
[7] | Suqi Xue, Farong Gao, Xudong Wu, Qun Xu, Xuecheng Weng, Qizhong Zhang . MUNIX repeatability evaluation method based on FastICA demixing. Mathematical Biosciences and Engineering, 2023, 20(9): 16362-16382. doi: 10.3934/mbe.2023730 |
[8] | Bei Liu, Wenbin Tan, Xian Zhang, Ziqi Peng, Jing Cao . Recognition study of denatured biological tissues based on multi-scale rescaled range permutation entropy. Mathematical Biosciences and Engineering, 2022, 19(1): 102-114. doi: 10.3934/mbe.2022005 |
[9] | Dae Hyeon Kim, Jin-Oh Park, Dae-Young Lee, Young-Seok Choi . Multiscale distribution entropy analysis of short epileptic EEG signals. Mathematical Biosciences and Engineering, 2024, 21(4): 5556-5576. doi: 10.3934/mbe.2024245 |
[10] | Linchao Yang, Ying Liu, Guanglu Yang, Shi-Tong Peng . Dynamic monitoring and anomaly tracing of the quality in tobacco strip processing based on improved canonical variable analysis and transfer entropy. Mathematical Biosciences and Engineering, 2023, 20(8): 15309-15325. doi: 10.3934/mbe.2023684 |
Investigation of the dynamic properties of implants is essential to ensure safety and compatibility with the host's natural spinal tissue. This paper presents a simplified model of a cantilever beam to investigate the effects of holes/pores on the structures. Free vibration test is one of the most effective methods to measure the dynamic response of a cantilever beam, such as natural frequency and damping ratio. In this study, the natural frequencies of cantilever beams made of polycarbonate (PC) containing various circular open holes were investigated numerically, analytically, and experimentally. The experimental data confirmed the accuracy of the natural frequencies of the cantilever beam with open holes calculated by finite element and analytical models. In addition, two finite element simulation methods, the dynamic explicit and modal dynamic methods, were applied to determine the damping ratios of cantilever beams with open holes. Finite element analysis accurately simulated the damped vibration behavior of cantilever beams with open holes when known material damping properties were applied. The damping behavior of cantilever beams with random pores was simulated, highlighting a completely different relationship between porosity, natural frequency and damping response. The latter highlights the potential of finite element methods to analyze the dynamic response of arbitrary and complex structures, towards improved implant design.
Labour and delivery is routinely monitored electronically with sensors that measure and record maternal uterine pressure
(UP) and fetal heart rate (FHR), a procedure referred to as cardiotocography (CTG). The objective of this monitoring is to detect the fetus at substantial risk of hypoxic injury so that intervention can prevent its occurrence.
Clinicians' interpretation of intra-partum CTG signals relies on the temporary decreases in FHR (FHR decelerations) in response to uterine contractions. FHR decelerations are due mainly to two contraction-induced events: (1) umbilical-cord compression and (2) a decrease in oxygen delivery through an impaired utero-placental unit. There is general consensus that deceleration depth, frequency and timing with respect to contractions are indicators of both the insult and the ability of the fetus to withstand it. In particular, later timing of decelerations, referred to as "late decelerations" are generally less reassuring to clinicians than the more innocuous (and rare) "early decelerations".
Assessing synchronization of the cardiotocography (CTG) signals uterine pressure (UP) and fetal heart rate (FHR) is a challenge both because of signal non-stationarity and noise in the form of signal discontinuities and maternal heart rate interference. As labour progress the maternal and fetal states evolve: contractions get more frequent and the cumulative effect of reduced oxygen transfer means that the fetus can experience increasing levels of hypoxia, progressing from respiratory to metabolic acidosis and finally to direct myocardial depression. Measurements under these conditions show that pH levels lower and base deficit increases: the buffer stores of the fetus which allow it to combat the increased acidosis become depleted [1]. Under such conditions, the state of the fetus as indicated by the FHR is quite non-stationary during labour and delivery.
We have characterized the stimulus-response system of UP and FHR using system identification (SI)[2]. However, obtaining delay estimates from SI involved lengthy and iterative modelling with a sweep of candidate delays. It would be helpful to have a better mechanism to estimate the delay before estimating the SI models. It would also be helpful to compare SI models with other synchronization estimates. In this study we used the model-free approaches of mutual information and transfer entropy to estimate synchronization and delays using both FHR and UP signals.
While there have been numerous studies examining differences in entropy measures in fetal heart rate variability (e.g., the use of approximate entropy in [3]), few fetal studies have applied mutual entropy or transfer entropy to this domain. Authors in [4] computed the transfer entropy of maternal and fetal heart rates using single lag analysis.
They calculated TE on maternal and fetal heart rates (derived from ECG) during gestation from 31 weeks onward. They observed a decrease in TE from fetus to mother over time and an increase in TE from mother to fetus over time. Delays in the M → F direction were estimated by TE analysis to be on the order of 5s. We used mutual entropy in [5] to examine CTG synchronization and extend that analysis to more robust mutual information estimates and to transfer entropy in this paper.
We used CTG from singleton, term pregnancies having no known congenital malformations, with ≥ 3 hours of tracing just prior to delivery. The CTG records for this study consisted of 42 pathological (P), 104 metabolic acidotic (MA) and 110 normal (N) fetuses. The pathological fetuses were severe enough to have confirmed evidence of hypoxic ischemic encephalopathy (<1st percentile of fetuses). The metabolic acidotic fetuses were defined by base deficit ≥ 12 mmol/L and with no apparent neurological injury (5th percentile of fetuses). It would be advantageous to identify these MA fetuses because they can be considered "close calls", where appropriate intervention can occur prior to the onset of injury. The cesarean section rates for the N and MA cases were 19.2% and 30% respectively, while the rate for the P cases was not consistently reported in our database. Data were provided from two US hospitals in compliance with institutional regulations.
The CTG data was recorded in a clinical setting, so it was subject to specific types of noise. The loss of sensor contact can temporarily interrupt the UP or FHR signals, and interference from the (much lower) maternal heart rate can corrupt the FHR. These both appeared in the signal as a sharp drop to much lower amplitude followed by a sharp signal restoration. As described in [6], we preprocessed the data to bridge interruptions less than 15s with linear interpolation. UP and FHR with signal loss greater than 15s were not included in the processed epochs and were removed from consideration.
As in [2], we detrended the signals by a high-pass FIR filter with sufficient extent (1991 filter coefficients) to pass a long contraction or deceleration without incurring excessive filter delay (there were 1991 coefficients). We chose a filter with a cutoff frequency of 1220s=4.5×10−3 Hz as a compromise between these competing demands. Finally, we decimated the UP and FHR signals to 0.25 Hz using an anti-aliasing low-pass filter with a cutoff frequency of 0.125Hz, before downsampling by 16. This was done to reduce the computational load of the entropy estimates described below.
Longer epochs generally result in superior numerical estimates while being more subject to the effects of non-stationarity. Following [2], we endeavoured to balance these competing effects by extracted 20-min epochs with 10-min overlap between successive epochs. This epoch length typically spanned several UP contraction-FHR deceleration pairs.
To address the issue of non-stationarity we extracted 20-min epochs with 10-min overlap between successive epochs. This epoch length typically spanned several UP contraction-FHR deceleration pairs. Epochs with insufficient valid signal were discarded.
The mutual information MI(X;Y) of random variables X and Y is defined as
MI(X;Y)=H(X)−H(Y|X)=∑x∑yp(x,y)logp(x,y)p(x)p(y) | (3.1) |
where H(X) is the entropy of X and H(X|Y)=H(X,Y)−H(Y) is the conditional entropy of X given Y. Mutual information is therefore the amount of uncertainty (i.e., entropy) about Y that is resolved by observing X. Reduction of uncertainty is equivalent to information [7]. Mutual information is symmetrical: MI(X;Y)=MI(Y;X). In this work we use the natural logarithm, giving units of nats.
In a time-series context, following the notation of [8], the mutual information can be reformulated as
MI(X;Y)=MI(XdXt−τ;Yt) | (3.2) |
Here X and Y refer to the mutual information between Y at time instant t, Yt and the state XdXt−τ, encapsulates the recent past of X, where τ>0 is a time lag in the information transfer and the dimension dX is the number of the lags in the state influencing Yt. The dimension dX is data dependent and limited in size only by candidate set of time lags C.
We calculated MI for the preprocessed UP and FHR signal pair u(t) and f(t). In previous work [5], we used binned estimates of the probability densities p(u), p(f) and p(u,f) in Eq 3.1. This binned approach restricted the delay search to a single lag, i.e., dX=1, from the source signal u(t) that maximized MI(u(t−τ),f). We repeated this calculation in this work, observing the time lag and value of the maximum MI over a range of possible UP-FHR alignments for each 20 min epoch. The candidate set C was the range of lags from -20 to 90 s (the negative lags are explained in the following paragraphs).
We chose these alignments to capture the widest range of reasonable physiological delays, empirically confirmed in [2]. We note also that UP is measured by tocography, which places a belt with a strain-gauge pressure sensor on the maternal abdomen. Because this pressure sensor is an indirect measure of intrauterine pressure [9], it may have an inherent mechanical delay, and the measured UP may reflect this with delays in onset and termination of contractions [2]. For this reason, we admitted negative alignments (i.e., -20 to 0 s).
In this work, we used a technique that admits multiple source lags to contribute to the MI; that is, where dX≥1. We computed this estimate using the Information Dynamics Toolkit (IDTxl) [10]. In accordance with the possible negative delay noted previously, we shifted the UP earlier by 20 s before estimating MI. We subsequently adjusted the lags analysis to account for this shift. Additionally, given the non-Gaussian nature of the UP and FHR distributions, we used the IDTxl implementation of the more robust k-nearest-neighbour algorithm [11] for the density estimates of Eq 3.1, with k set to the default value of 4. Finally, we tested the MI estimates for significance using 200 realizations of surrogate target data (i.e., FHR), that were generated by IDTxl by random permutations of the epoch samples (p<0.05 was considered significant). That is, if the MI with the UP corresponding to the measured FHR is greater than the value of MI of at least 190 other FHR permutations, it is significant. The surrogate testing ensures that spurious interactions in the data have a low probability (i.e., 5%) of being considered significant.
IDTxl chooses an embedding set Z of dX delays in a greedy fashion, following [12]. At each step, the delay from a candidate set C generating the highest MI and satisfying statistical significance is added to set Z (and removed from C). The search is repeated using the reduced set C to add additional delays to Z that increase the estimated MI, subject to significance testing. This continues until adding more delays to the set fails to significantly increase the MI.
Schreiber introduced the concept of conditional mutual information, also known as transfer entropy TE) to address the question "What information does the past of X provide about the future of Y, that the past of Y did not already provide?" [13]. The transfer entropy from random variable X to Y is denoted TE(X→Y) and is not symmetrical, but directed: In general, TE(X→Y)≠TE(Y→X). In the time-series context, it is defined as:
TE(X→Y)=MI(Yt;XdXt−τ;YdYt−1)=H(YdYt−1,XdXt−τ)−H(Yt,YdYt−1,XdXt−τ)+H(Yt,YdYt−1)−H(YdYt−1) | (3.3) |
Like mutual information, transfer entropy measures the amount of uncertainty about Y that is resolved by observing the state XdXt−τ, but it also removes the influence of the redundancy of the state YdYt−1 on Y. This can decode influences of X on Y that are masked in the MI calculation. While this formulation of TE allows a variable transfer delay of τ in the source X (as in Eq 3.2 for MI), it fixes the target lag to 1 to fulfill what is referred to as TE self prediction optimality [14], that is, all the recent past of Y is considered in the state YdYt−1. For TE, source and target have variable state dimensions, dX and dY, respectively.
Wibral offers a formal proof that the TE is maximized when τ=δ, the true interaction delay [14]. This proof of identifiability of the true delay in the information transfer holds strictly for the case of zero observation noise. However, with numerous numerical simulations with nonlinear dynamical systems, they demonstrate that the TE precision degrades gracefully under conditions of decreasing SNR.
For each 20-min epoch of each CTG recording, we calculated the binned version of MI and the nearest-neighbour versions of MI and TE and their associated lag(s). We assessed the differences of the MI, TE and their associated delays for the three fetal classes N, MA and P at each epoch over the last three hours of labour and delivery. We found that a single lag often dominated contributions to MI and TE and so we summarized lag reporting with these single values. The Kolgomorov-Smirnov (KS) test was used to test the null hypothesis that the distributions of pairwise-selected classes were the same, using the significance threshold p<0.05.
In addition to the full-epoch estimates, we observed "local" (i.e., "instantaneous") MI and TE estimates mi(t) and te(t) provided by IdTxl that indicate the relative coupling of individual samples within the epoch. As discussed in [15], with sufficient data, overall MI and TE values are non-negative, but these local estimates can be negative and are referred to as "misinformation”, corresponding to the opposite of coupling (i.e., indicating an element of surprise between input and output).
In the results below, we refer to the binned estimates of MI as MIb and the Kraskov estimates of MI as MIk. We refer to the single delay with the maximal contribution to MIk, that is, arg maxτ MIk(τ) as τMIk_max. We refer to the overall TE with contributions from all significant lags as TE. We refer to the single delay with the maximal contribution to TE, that is arg maxτ TE(τ) as τTE_max and the TE contribution of that lag as TEmax.
While all analysis was done using lag delays with the time unit of a 0.25 Hz sample, for clarity we convert these values to seconds in the subsequent time plots.
Figure 1 shows typical histograms of UP and FHR for a 20 min epoch. These approximations indicate that the univariate and joint densities are non-Gaussian, supporting the motivation to use the Kraskov density estimator.
Figures 2 and 3 show calculations of the binned and Kraskov MIs and the Kraskov TE for typical normal and pathological fetuses, demonstrating for these two cases the longer delay of the pathological case.
The normal case has a maximal Kraskov MI delay at 0 s with secondary delays at 12, 84 and 90 s. The TE is maximal at a sole delay of -4s. Visual inspection of the delay between UP contraction and FHR deceleration onsets in this epoch indicate that a near-zero delay is plausible. On the other hand, the binned MI value has an erroneous peak at 70 s.
The pathological case has a maximal Kraskov MI delay at 64 s and a secondary delay at 36 s. Visual inspection of the delay in this epoch indicates that these two delays are plausible. The TE is maximal at the shorter delay of 36s only. The binned MI value peaks at 70 s, closer to the longer delay.
Table 1 summarizes the processing of MI and TE. MI or TE estimates were obtained when at least one significant delay was estimated. Estimates can be missing when no significant delays are found during the search. This can be the result of the absence of coupling from either signal loss and/or noise as discussed in the Preprocessing section. We obtained ME estimates for P, MA and N fetuses in 62.4%, 68.2% and 72.6% of the epochs, respectively. For TE these numbers were lower at 43.5%, 36.6% and 35.2%, respectively. The ratio of TE/MI estimates was higher for P cases: 69.7%, vs 53.6% and 48.5% for MA and N, respectively.
Fetal class | P | MA | N |
Number of recordings | 42 | 104 | 110 |
Total epochs | 756 | 1872 | 1980 |
MI: epochs estimated | 472 | 1277 | 1437 |
TE: epochs estimated | 329 | 685 | 697 |
MI: epochs estimated (%) | 62.4% | 68.2% | 72.6% |
TE: epochs estimated (%) | 43.5% | 36.6% | 35.2% |
Ratio TE/MI estimated | 69.7% | 53.6% | 48.5% |
Figure 4 shows per-class binned MI estimates MIb during the final 3 hours of delivery. The difference between N and P is statistically significant between 140 and 30 minutes before delivery. N-MA differences are significant slightly later, between 110 and 10 minutes before delivery. There is a tendency from N to MA to P fetuses to have higher MI. The delay at maximal MI (not shown) was not discriminating.
Figure 5 shows per-class Kraskov MI estimates MIK during the final 3 hours of delivery. The difference between N and P is statistically significant between 160 and 40 minutes before delivery. N-MA differences are significant in two of these epochs, at 110 minutes and 70 minutes. There is a tendency from N to MA to P fetuses to have higher MI. Furthermore, two epochs in the final 50 minutes of delivery show longer P delays τMIkmax with statistical significance, while one epoch did so at 50 minutes for N-MA comparisons. There is a tendency from N to MA to P fetuses to have longer delays.
Figure 6 shows per-class Kraskov TE estimates during the final 3 hours of delivery. The delay with the maximal contribution to TE, τTE_max shows differences between N and P that are statistically significant in 4 epochs between 100 and 20 minutes before delivery. There is a tendency from N to MA to P fetuses towards lower TE. Furthermore, one epoch in the final 20 minutes of delivery shows longer P delays with statistical significance. The overall TE value TE shows a similar pattern but with few epochs having significant differences. The delay where the contribution to TE is greatest, TEmax, shows differences in the final 50 minutes of delivery, but only 1 epoch was statistically significant (at 20 minutes before delivery). Again there is a tendency form N to MA to P fetuses to have longer delays in these final epochs, although N-MA comparisons never reached statistical significance.
The delay estimates in the typical epochs shown were plausible and in accordance with the clinical expectation that pathology tends to be associated with longer delays between the mother and the fetus. These results are borne out at the larger scale of per-class comparisons over time, with statistically significant differences occurring early for both N-P cases but crucially for clinical significance, for N-MA comparisons. Detecting MA cases early, before injury, allows time for appropriate cesarian section to occur. It was also the case that MI and TE estimates were dominated by a single delay, and this was the rationale for displaying the argmax plots of Figures 5 and 6.
The typical epochs also illustrate that even with restricting the epoch length to 20 minutes, non-stationarity is present within the epoch, with delays occurring at at least two lags in the examples shown. This clearly introduces noise into the per-class comparisons that consider a single lag, and this noise diminishes discrimination.
Furthermore, we did not explicitly compare the degree of signal artefacts in this study among groups, but the previous study [2] shared the same pathological cases and the figure of 42% artifact for pathologicals vs. 8% for normal reported there is a valid comparison for this study. This can be a source of bias in our analysis, but the fact that we excluded very long artefacts from consideration diminishes this bias.
But despite these effects of non-stationarity and noise, the estimates do succeed about two-thirds of the time for MI and less so for TE, suggesting the degradation effect shown in [14]. The robust nearest neighbour Kraskov density estimator and the surrogate tests give confidence that the results are indeed significant. The significance of the generally lower TE success rates and the higher TE/MI ratio of epochs processed for P cases is interesting and merits further investigation, especially since TE values tended to be lower for these fetuses.
The increasing MI with fetal pathology was expected from our previous study [5]. This is plausibly explained by the fact that the fetus in distress reacts strongly, with deeper decelerations, to the onslaught of maternal contractions. But the decreasing TE with pathology was an unforeseen and novel result and cannot be explained merely by the fact that decelerations increase as labour progresses. This seems to indicate that the balance between MI and TE tips in the direction of more and more signal properties being (information-theoretically) redundant between source UP and target FHR. In contrast, more complex, FHR target-state dependent processing seems to vanish, possible indicating some regulatory mechanisms may be compromised. It is consistent with a view of a fetus in distress having difficulty compensating neurologically to the effects of labour, especially persistent contractions.
We have demonstrated the utility of two new measures of perinatal fetal state, mutual entropy and transfer entropy, that may have complementary information for estimating the nature and timing of the fetal response to labour. The observed statistically significant differences indicate that these indices may be useful in future work for the development of a fetal state classifier.
Following from this study, we would like to assess and compare transfer entropy during specific baseline, acceleration and deceleration time intervals of FHR. We also hope to include maternal heart rate to obtain multivariate estimates of TE, which are currently supported by IDTxl. Further, we hope to compare fetus → maternal TE, both for its inherent interest and also to aid in significance testing. We have assumed in this study that maternal → fetal causality likely predominates, but this may be a matter of degree that is informative.
Finally, Wibral refers to "anticipative synchronization" occurring when a slave system (Y) can anticipate the dynamics of a master system X when X is subject to a long feedback loop [14]. By testing TE (X→X) it is possible to tests for "feedback" in the source signal. This may be a way to overcome the challenges of detecting phase with signals like CTG with their inherent periodicity and under conditions of small sample sizes. This could benefit TE analysis with UP as a source in particular, especially at the end of labour when contractions become more frequent.
This research was supported by PeriGen Inc. The authors would like to thank Patricia Wollstadt and Michael Wibral for providing the open-source IDTxl package and for kindly assisting in the setup of the experiments and interpretation of the results.
PeriGen Inc. develops clinical decision support tools for cardiotocography.
[1] |
Y. Matsumoto, M. J. Griffin, Dynamic response of the standing human body exposed to vertical vibration: Influence of posture and vibration magnitude, J. Sound Vib., 212 (1998), 85–107. https://doi.org/10.1006/jsvi.1997.1376 doi: 10.1006/jsvi.1997.1376
![]() |
[2] | A. Bogadi-Sare, The effect of whole-body vibration: An unrecognized medical problem, Arh. Hig. Rada Toksikol., 44 (1993), 269–279. |
[3] |
L. X. Guo, M. Zhang, Y. M. Zhang, E. C. Teo, Vibration modes of injured spine at resonant frequencies under vertical vibration, Spine, 34 (2009), E682–E688. https://doi.org/10.1097/BRS.0b013e3181b1fdf4 doi: 10.1097/BRS.0b013e3181b1fdf4
![]() |
[4] |
G. Marini, G. Huber, K. Puschel, S. J. Ferguson, Nonlinear dynamics of the human lumbar intervertebral disc, J. Biomech., 48 (2015), 479–488. https://doi.org/10.1016/j.jbiomech.2014.12.006 doi: 10.1016/j.jbiomech.2014.12.006
![]() |
[5] |
Y. Matsumoto, M. J. Griffin, Movement of the upper-body of seated subjects exposed to vertical whole-body vibration at the principal resonance frequency, J. Sound Vib., 215 (1998), 743–762. https://doi.org/10.1006/jsvi.1998.1595 doi: 10.1006/jsvi.1998.1595
![]() |
[6] |
Y. Matsumoto, M. J. Griffin, Modelling the dynamic mechanisms associated with the principal resonance of the seated human body, Clin. Biomech., 16 (2001), S31–S44. https://doi.org/10.1016/S0268-0033(00)00099-1 doi: 10.1016/S0268-0033(00)00099-1
![]() |
[7] |
W. Fan, L. X. Guo, The role of posterior screw fixation in single-level transforaminal lumbar interbody fusion during whole body vibration: A finite element study, World Neurosurg., 114 (2018), E1086–E1093. https://doi.org/10.1016/j.wneu.2018.03.150 doi: 10.1016/j.wneu.2018.03.150
![]() |
[8] | L. Du, X. J. Sun, T. J. Zhou, Y. C. Li, C. Chen, C. Q. Zhao, et al., The role of cage height on the flexibility and load sharing of lumbar spine after lumbar interbody fusion with unilateral and bilateral instrumentation: A biomechanical study, BMC Musculoskeletal Disord., 18 (2017), 1–8. |
[9] |
P. Schleicher, R. Gerlach, B. Schar, C. Cain, W. Achatz, R. Pflugmacher, et al., Biomechanical comparison of two different concepts for stand alone anterior lumbar interbody fusion, Eur. Spine J., 17 (2008), 1757–1765. https://doi.org/10.1007/s00586-008-0797-4 doi: 10.1007/s00586-008-0797-4
![]() |
[10] |
D. U. Erbulut, A. Kiapour, T. Oktenoglu, A. F. Ozer, V. K. Goel, A computational biomechanical investigation of posterior dynamic instrumentation: Combination of dynamic rod and hinged (dynamic) screw, J. Biomech. Eng, 136 (2014), 51007. https://doi.org/10.1115/1.4027060 doi: 10.1115/1.4027060
![]() |
[11] |
G. Marini, G. Huber, K. Puschel, S. J. Ferguson, A 1-D model of the nonlinear dynamics of the human lumbar intervertebral disc, J. Sound Vib., 387 (2017), 194–206. https://doi.org/10.1016/j.jsv.2016.09.021 doi: 10.1016/j.jsv.2016.09.021
![]() |
[12] | V. Kumar, K. K. Singh, S. Gaurav, Analysis of natural frequencies for cantilever beam with I- and T- section using Ansys, Int. Res. J. Eng. Technol., 2 (2015), 1013–1020. |
[13] |
A. Suzuki, T. Hoshiai, H. Nakata, T. Otomaru, M. Oki, H. Taniguchi, et al., Modal analysis of two different types of fixed implant-supported prostheses embedded in edentulous maxillae, J. Prosthodontic Res., 63 (2019), 327–333. https://doi.org/10.1016/j.jpor.2019.01.006 doi: 10.1016/j.jpor.2019.01.006
![]() |
[14] |
H. B. Khaniki, M. H. Ghayesh, S. Hussain, M. Amabili, Porosity, mass and geometric imperfection sensitivity in coupled vibration characteristics of CNT-strengthened beams with different boundary conditions, Eng. Comput., 38 (2022), 2313–2339. https://doi.org/10.1007/s00366-020-01208-3 doi: 10.1007/s00366-020-01208-3
![]() |
[15] |
H. B. Khaniki, M. H. Ghayesh, S. Hussain, M. Amabili, Effects of geometric nonlinearities on the coupled dynamics of CNT strengthened composite beams with porosity, mass and geometric imperfections, Eng. Comput., 38 (2022), 3463–3488. https://doi.org/10.1007/s00366-021-01474-9 doi: 10.1007/s00366-021-01474-9
![]() |
[16] |
M. Amir, M. Talha, Influence of large amplitude vibration on geometrically imperfect sandwich curved panels embedded with gradient metallic cellular core, Int. J. Appl. Mech., 12 (2020), 2050099. https://doi.org/10.1142/S1758825120500994 doi: 10.1142/S1758825120500994
![]() |
[17] |
H. S. N. Kumar, S. Kattimani, T. Nguyen-Thoi, Influence of porosity distribution on nonlinear free vibration and transient responses of porous functionally graded skew plates, Def. Technol., 17 (2021), 1918–1935. https://doi.org/10.1016/j.dt.2021.02.003 doi: 10.1016/j.dt.2021.02.003
![]() |
[18] |
H. B. Khaniki, M. H. Ghayesh, R. Chin, M. Amabili, Large amplitude vibrations of imperfect porous-hyperelastic beams via a modified strain energy, J. Sound Vib., 513 (2021), 116416. https://doi.org/10.1016/j.jsv.2021.116416 doi: 10.1016/j.jsv.2021.116416
![]() |
[19] |
D. L. Rao, Z. Q. Xu, Damping behavior of 304L stainless steel after fatigue loading, Theor. Appl. Fract. Mech., 100 (2019), 110–113. https://doi.org/10.1016/j.tafmec.2019.01.006 doi: 10.1016/j.tafmec.2019.01.006
![]() |
[20] |
A. J. Molina-Viedma, E. Lopez-Alba, L. Felipe-Sese, F. A. Diaz, Full-field modal analysis during base motion excitation using high-speed 3D digital image correlation, Meas. Sci. Technol., 28 (2017). https://doi.org/10.1088/1361-6501/aa7d87 doi: 10.1088/1361-6501/aa7d87
![]() |
[21] |
C. C. Ma, M. Y. Shao, J. C. Ma, C. C. Liu, K. M. Gao, Fluid structure interaction analysis of flexible beams vibrating in a time-varying fluid domain, Proc. Inst. Mech. Eng., Part C: J. Mech., 234 (2020), 1913–1927. https://doi.org/10.1177/0954406220902163 doi: 10.1177/0954406220902163
![]() |
[22] |
B. Yang, Theory of vibration | fundamentals, Encycl. Vib., 2001 (2001), 1290–1299. https://doi.org/10.1006/rwvb.2001.0112 doi: 10.1006/rwvb.2001.0112
![]() |
[23] | W. T. Thomson, Theory of Vibration with Applications, 4th edition, Taylor & Francis, 1993. |
[24] |
H. I. Yoon, I. S. Son, S. J. Ahn, Free vibration analysis of Euler-Bernoulli beam with double cracks, J. Mech. Sci. Technol., 21 (2007), 476–485. https://doi.org/10.1007/Bf02916309 doi: 10.1007/Bf02916309
![]() |
[25] | S. U. Sawant, S. J. Chauhan, N. N. Deshmukh, Effect of crack on natural frequency for beam type of structures, in International Conference on Functional Materials, Characterization, Solid State Physics, Power, Thermal and Combustion Energy (Fcsptc-2017), 1859 (2017), 20056. https://doi.org/10.1063/1.4990209 |
[26] |
X. Du, G. Blugan, T. Kunniger, S. S. Lee, L. Vladislavova, S. J. Ferguson, Non-linear mechanical properties and dynamic response of silicon nitride bioceramic, Ceram. Int., 47 (2021), 33525–33536. https://doi.org/10.1016/j.ceramint.2021.08.261 doi: 10.1016/j.ceramint.2021.08.261
![]() |
[27] |
Z. Kiral, B. M. Icten, B. G. Kiral, Effect of impact failure on the damping characteristics of beam-like composite structures, Compos. Part B Eng., 43 (2012), 3053–3060. https://doi.org/10.1016/j.compositesb.2012.05.005 doi: 10.1016/j.compositesb.2012.05.005
![]() |
[28] | R. A. W. Jane, K. Cullum, Lanczos Algorithms for Large Symmetric Eigenvalue Computations, Birkhaeuser, Boston, 1985. |
[29] |
C. Kyriazoglou, F. J. Guild, Quantifying the effect of homogeneous and localized damage mechanisms on the damping properties of damaged GFRP and CFRP continuous and woven composite laminates—an FEA approach, Compos. Part A Appl. Sci. Manuf., 36 (2005), 367–379. https://doi.org/10.1016/j.compositesa.2004.06.037 doi: 10.1016/j.compositesa.2004.06.037
![]() |
[30] |
Z. H. Zhu, S. A. Meguid, Nonlinear FE-based investigation of flexural damping of slacking wire cables, Int. J. Solids Struct., 44 (2007), 5122–5132. https://doi.org/10.1016/j.ijsolstr.2006.12.024 doi: 10.1016/j.ijsolstr.2006.12.024
![]() |
[31] |
T. Zhang, Y. Z. Zuo, H. W. Teng, D. Y. Ma, H. Liu, Study on Rayleigh damping in dynamic analysis for story-adding RC structure of light-weight steel, Appl. Mech. Mater., 482 (2014), 123–128. https://doi.org/10.4028/www.scientific.net/amm.482.123 doi: 10.4028/www.scientific.net/amm.482.123
![]() |
[32] |
M. Liu, D. G. Gorman, Formulation of Rayleigh damping and its extensions, Comput. Struct., 57 (1995), 277–285. https://doi.org/10.1016/0045-7949(94)00611-6 doi: 10.1016/0045-7949(94)00611-6
![]() |
[33] | V. S. Geraschenko, A. S. Grishin, N. I. Gartung, Approaches for the calculation of Rayleigh damping coefficients for a time–history analysis, in Structures Under Shock and Impact XV, (2018), 227–237. |
[34] |
L. M. Khoo, P. R. Mantena, P. Jadhav, Structural damage assessment using vibration modal analysis, Struct. Health Monit., 3 (2004), 177–194. https://doi.org/10.1177/1475921704042680 doi: 10.1177/1475921704042680
![]() |
[35] |
H. Rajoria, N. Jalili, Passive vibration damping enhancement using carbon nanotube-epoxy reinforced composites, Compos. Sci. Technol., 65 (2005), 2079–2093. https://doi.org/10.1016/j.compscitech.2005.05.015 doi: 10.1016/j.compscitech.2005.05.015
![]() |
[36] |
V. G. Geethamma, R. Asaletha, N. Kalarikkal, S. Thomas, Vibration and Sound damping in polymers, Resonance, 19 (2014), 821–833. https://doi.org/10.1007/s12045-014-0091-1 doi: 10.1007/s12045-014-0091-1
![]() |
[37] |
A. M. Puthanpurayil, O. Lavan, A. J. Carr, R. P. Dhakal, Elemental damping formulation: An alternative modelling of inherent damping in nonlinear dynamic analysis, Bull. Earthquake Eng., 14 (2016), 2405–2434. https://doi.org/10.1007/s10518-016-9904-9 doi: 10.1007/s10518-016-9904-9
![]() |
[38] |
M. Wesolowski, E. Barkanov, Improving material damping characterization of a laminated plate, J. Sound Vib., 462 (2019), 114928. https://doi.org/10.1016/j.jsv.2019.114928 doi: 10.1016/j.jsv.2019.114928
![]() |
[39] |
L. Viswanadham, R. N. Rao, C. S. Chaitanya, Effect of material addition on the vibration response of a cantilever beam, Mater. Today Proc., 18 (2019), 4537–4541. https://doi.org/10.1016/j.matpr.2019.07.425 doi: 10.1016/j.matpr.2019.07.425
![]() |
[40] |
L. Gagnon, M. Morandini, G. L. Ghiringhelli, A review of friction damping modeling and testing, Arch. Appl. Mech., 90 (2020), 107–126. https://doi.org/10.1007/s00419-019-01600-6 doi: 10.1007/s00419-019-01600-6
![]() |
[41] | I. T. C. K. M. Stamatopoulos, S. D. Panteliou, Damping associated with porosity in porous rectangular plates, in Proceedings of the Ninth International Conference on Computational Structures Technology, Civil-Comp Press, (2008), 115. |
[42] |
Q. Y. Li, G. F. Jiang, J. Dong, J. W. Hou, G. He, Damping behavior and energy absorption capability of porous magnesium, J. Alloys Compd., 680 (2016), 522–530. https://doi.org/10.1016/j.jallcom.2016.04.101 doi: 10.1016/j.jallcom.2016.04.101
![]() |
[43] |
W. G. Zheng, S. M. He, R. J. Tang, S. L. He, Damping enhancement using axially functionally graded porous structure based on acoustic black hole effect, Materials, 12 (2019), 2480. https://doi.org/10.3390/ma12152480 doi: 10.3390/ma12152480
![]() |
[44] |
I. S. Golovin, H. R. Sinning, I. K. Arhipov, S. A. Golovin, M. Bram, Damping in some cellular metalic materials due to microplasticity, Mater. Sci. Eng. A, 370 (2004), 531–536. https://doi.org/10.1016/j.msea.2003.08.089 doi: 10.1016/j.msea.2003.08.089
![]() |
[45] |
H. Matsuzaki, Y. Tokuhashi, E. Matsumoto, M. Hoshino, T. Kiuchi, S. Toriyama, Problems and solutions of pedicle screw plate fixation of lumbar spine, Spine, 15 (1990), 1159–1165. https://doi.org/10.1097/00007632-199011010-00014 doi: 10.1097/00007632-199011010-00014
![]() |
[46] |
X. L. Deng, H. X. Chen, Q. Y. Xu, F. Feng, X. Y. Chen, X. W. Lv, et al., Topology optimization design of three-dimensional multi-material and multi-body structure based on irregular cellular hybrid cellular automata method, Sci. Rep., 12 (2022), 1–20. https://doi.org/10.1038/s41598-022-09249-y doi: 10.1038/s41598-022-09249-y
![]() |
[47] |
M. Fantini, M. Curto, F. D. Crescenzio, A method to design biomimetic scaffolds for bone tissue engineering based on Voronoi lattices, Virtual Phys. Prototyping, 11 (2016), 77–90. https://doi.org/10.1080/17452759.2016.1172301 doi: 10.1080/17452759.2016.1172301
![]() |
1. | Michael Moncerrat Piña-Martínez, Adriana Cristina Pliego-Carrillo, Claudia Ivette Ledesma-Ramírez, Juan Carlos Echeverría, Hugo Mendieta-Zerón, Laura Mercedes Santiago-Fuentes, Eric Alonso Abarca-Castro, Ana Karen Talavera-Peña, José Javier Reyes-Lagos, Reduced maternal–fetal cardiac coupling in moderate preterm labor: a high-resolution symbolic dynamics analysis, 2024, 1951-6355, 10.1140/epjs/s11734-024-01426-6 |
Fetal class | P | MA | N |
Number of recordings | 42 | 104 | 110 |
Total epochs | 756 | 1872 | 1980 |
MI: epochs estimated | 472 | 1277 | 1437 |
TE: epochs estimated | 329 | 685 | 697 |
MI: epochs estimated (%) | 62.4% | 68.2% | 72.6% |
TE: epochs estimated (%) | 43.5% | 36.6% | 35.2% |
Ratio TE/MI estimated | 69.7% | 53.6% | 48.5% |