Next Article in Journal
Bone Bruises and Concomitant Meniscus and Cartilage Damage in Anterior Cruciate Ligament Injuries: A Systematic Review and Meta-Analysis
Previous Article in Journal
Modulation of Canine Adipose-Derived Mesenchymal Stem/Medicinal Signalling Cells with Ascorbic Acid: Effect on Proliferation and Chondrogenic Differentiation on Standard Plastic and Silk Fibroin Surfaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Analysis of Stress–Stretch Curves in Canine Lumbar Vertebrae Using Modified Logistic Functions

by
Ernest Kostenko
1,*,
Rimantas Stonkus
2,
Jakov Šengaut
3,
Nikolaj Višniakov
4 and
Algirdas Maknickas
4,5,*
1
Department of Veterinary, Faculty of Agrotechnologies, Vilniaus Kolegija/Higher Education Institution, 08105 Vilnius, Lithuania
2
Department of Mechatronics, Robotics and Digital Manufacturing, Vilnius Gediminas Technical University, 10105 Vilnius, Lithuania
3
Jakov’s Veterinary Centre, 03147 Vilnius, Lithuania
4
Institute of Mechanical Science, Vilnius Gediminas Technical University, 10105 Vilnius, Lithuania
5
Department of Biomechanical Engineering, Vilnius Gediminas Technical University, 10105 Vilnius, Lithuania
*
Authors to whom correspondence should be addressed.
Bioengineering 2024, 11(5), 516; https://doi.org/10.3390/bioengineering11050516
Submission received: 16 April 2024 / Revised: 3 May 2024 / Accepted: 15 May 2024 / Published: 20 May 2024
(This article belongs to the Section Biomechanics and Sports Medicine)

Abstract

:
Background: The mechanical characteristics of bone are crucial for comprehending its functionality and response to different load conditions, which are essential for advancing medical treatments, implants, and prosthetics. By employing mathematical modeling to analyze the mechanical properties of bone, we can assess stress and deformation under both normal and abnormal conditions. This analysis offers valuable perspectives on potential fracture risks, the effects of diseases, and the effectiveness of various treatments. Therefore, researchers are attempting to find an adequate mathematical description of the mechanical properties of bone. Methods: Experimental stress–stretch external loading curves were obtained through investigations of canine vertebrae. The obtained experimental curves were fitted using the SciPy Python library with a slightly modified logistic function (logistic function plus additional const). Results: The resulting coefficient of determination  R 2  (R squared) for most curves was near 0.999, indicating that an appropriate fitting function was selected for the description of the experimental stress–stretch curves. Conclusions: The stress–stretch behavior of canine vertebrae can be described using a logistic function modified by adding additional parameters for the most accurate fitting results.

Graphical Abstract

1. Introduction

The mechanical properties of bone are fundamental to understanding its function and behavior under various loading conditions, as well as to the development of medical treatments, implants, and prosthetics. Mathematical modeling of the mechanical properties of bone allows for the analysis of stress and stretch under physiological and pathological conditions, providing insights into fracture risks, disease impacts, and treatment outcomes. Previous research has explored various mathematical models that describe the stress–stretch behavior of bone, each catering to different aspects of bone mechanics and applicable to a range of scenarios from everyday activities to traumatic impacts. Depending on the error tolerance of the model considered, it can be useful to differentiate between linear and non-linear mathematical models of mechanical stress–strain (or stretch for large deformation) of bone, including isotropic or orthotropic elastic, elastic–plastic, viscoelastic, and poroelastic models [1,2,3,4,5,6,7,8].
For non-linear modeling of biological tissues, particularly when considering large deformations that exceed the limits of traditional linear elastic models, several hyperelastic models can be applied effectively. These models are capable of capturing the complex, non-linear stress–stretch behavior of bone and other biological tissues under various loading conditions. The choice of model depends on the specific characteristics of the bone tissue studied and the nature of the deformation. The neo-Hookean model [9] is often used as a first approximation in biological tissue modeling. Its simplicity makes it appropriate for initial studies, although it may not capture the full complexity of the mechanical responses of bone. The Mooney–Rivlin model extends the neo-Hookean framework by adding an additional term to better accommodate the non-linear characteristics of materials under compression and tension. This model is more flexible and can more accurately represent the mechanical behavior of soft tissues and rubber-like materials. Its ability to fit experimental data for various deformation modes makes it a valuable tool for more detailed studies of bone and other biological tissues [10]. The Ogden model is particularly useful for capturing the complex, non-linear elastic behavior of anisotropic materials, such as bone, where the response to loading can vary significantly in different directions [10]. The Arruda–Boyce model [11] offers a micro-mechanical approach to hyperelasticity based on the statistical mechanics of polymer chains. Although this model was originally developed for polymers, its principles can be applied to modeling the non-linear elastic behavior of biological tissues, including bone, especially in contexts where the microstructure of the material plays a crucial role in its mechanical response. The Yeoh model [12] starts with a simple form that accurately captures the initial stiffness and adds higher-order terms as required to account for material responses under larger stretches. This model is useful for capturing the non-linear behavior of bone, particularly when focusing on compressive loading conditions.
The weakness of the above-mentioned non-linear bone tissue models is that they overlap in the lower, middle, or upper part of the stretch range but not over the full stretch interval from one to the value at which ultimate stress occurs. This has impelled researchers to identify the best-fitting curve for an accurate description of vertebra stress–stretch curves in external load experiments using canine vertebrae. The first such attempt involved the application of seventh-order polynomials for stress–time and second-order polynomials for stretch–time relations in a parametric form of the stress–stretch curve [13]. During the investigation of soft biological tissues, such as a canine aorta, in opposite tension tests, we have found that the stress–stretch relations of the aorta can be approximated by a logistic function [14]. The aim of the present research is to propose a new non-linear model based on a logistic function with four parameters appropriate for the approximation of stress–stretch relations in canine vertebra load tests.

2. Materials and Methods

2.1. Research Objects

In this study, we utilized lumbar vertebrae samples from three spayed mongrel canines, aged eight, nine, and ten years, with weights of 28, 26, and 20 kg, respectively. The first canine had been diagnosed with a mammary tumor, and the second with a pulmonary thromboembolism. Following these diagnoses, the owners of the dogs consented to euthanasia, and all necessary consent forms were completed accordingly. The third canine exhibited a mass located on the right side of the groin region. The owners of this dog expressed a lack of interest in determining the structure and cause of this tumor and consequently made the decision to euthanize the animal.

2.2. Vertebrae Preparation and Processing

The animals underwent an autopsy, wherein the lumbar vertebrae from L1 to L7 were extracted, as shown in Figure 1. The dissected segment was then stored at a temperature of  20 ° C  in a freezer. For biomechanical analysis, this segment was slowly returned to room temperature ( 22 ° C ). After thawing, the segment underwent processing, during which the surrounding muscular and adipose connective tissues were carefully removed. Each vertebra was individually isolated, and the spinal cord was excised.

2.3. Mechanical Experiment

Two types of equipment were used in our experiments. For load until fracture of vertebra experiments, we employed a 2055 P-5 universal tensile testing machine (Tochpribor, Ivanovo, Russia) equipped with a compression testing attachment. The canine lumbar vertebrae were subjected to vertical compressive forces through a mechanical loading system controlled by “LabVIEW version 16” (National Instruments, Austin, TX, USA). The vertebrae were aligned such that their longitudinal axes were perpendicular to the direction of the applied load, facilitating the accurate and direct force application necessary for biomechanical evaluation. Next, the compression tests utilized PXI system hardware, including the chassis NI PXIe-1073 and controller PXIe-4330 (both from Austin, TX, USA). An S-type tension/compression load cell with a maximum capacity of 1 kN was utilized. We varied the testing load velocity for different tests over the range  1 10 m m m i n . Bone samples were frozen for storage and preservation, then thawed before testing; accordingly, our experiments used thawed wet specimens.
For the second group of tests, lumbar vertebrae were examined using a Mecmesin MultiTest 2.5-i micro-compression machine (Mecmesin Limited, Slinfold, UK) (Figure 2). The Mecmesin AFG25 cell, which controls deformation, achieved a measurement precision of within  ± 0.01  mm. It applied compressive forces ranging from 2 N to 2500 N with a precision of  ± 0.1 %  and maintained a crosshead speed with the same level of accuracy. During this process, each vertebra underwent compression at a force of 950 N and at a rate of 1 mm per minute. The compression process was repeated ten times for each vertebra. Post-experimentation, the axes of each vertebra were measured, and cyclic loading was applied when the specimen was in place.

2.4. Error Estimation

The coefficient of determination, often denoted as  R 2 , is a statistical measure used in the context of statistical models. Its main purpose is to predict future outcomes or test hypotheses based on related information. The coefficient of determination,  R 2 , is defined as follows:
R 2 = 1 S S res S S tot
where  S S res  is the sum of squares of residuals, also known as the residual sum of squares:
S S res = i = 1 n ( y i y ^ i ) 2
and  S S tot  is the total sum of squares (a measure of the total variance in the observed data):
S S tot = i = 1 n ( y i y ¯ ) 2
where  y i  is the observed outcome,  y ^ i  is the predicted value for the i-th observation, and  y ¯  is the mean value of the observed data.
In practice, the coefficient of determination is used as an indicator of the goodness of fit of a model. A higher  R 2  value indicates a better fit and suggests that the model can better explain the variation in the output using different independent variables. However, it should not be the sole criterion for model selection, especially because it can increase with the number of predictors without improving the predictive power of the model.
The confidence intervals for the lengths and loading surfaces of the measured vertebrae were calculated using the formula provided in reference [15], as follows:
δ x = ± t s / n
where t represents the Student’s t-value, s is the standard deviation, and n denotes the degrees of freedom plus one. For two degrees of freedom, the Student’s t-value is set at 4.303, corresponding to a 95% confidence interval [15]. Lastly, the error propagation for functions involving two variables is applied to the external pressure-induced stress, denoted by  σ , as follows:
σ = F / S ,
where F is the measured load force and S is the vertebra loading surface.  λ  is the stretch, which is calculated as
λ = ( L l ) / L ,
where L is vertebra length between loading surfaces and can be expressed as follows [16]:
δ σ / | σ | = ( δ F / F ) 2 + ( δ S / S ) 2
δ λ = δ l / L 2

2.5. Approximation of Stress–Stretch Curves

Parametric equations are often employed to describe the coordinates of points forming a geometric entity, such as a curve. In this context, the equations are jointly referred to as a parametric representation [17,18]. The experimental data-gathered stress  σ k ( t )  on the increasing and decreasing external load and displacement,  λ k ( t ) , of the vertebra can be effectively represented using the following parametric formulas:
σ k ( t ) = C 0 k + C 1 k 1 + exp ( C 2 k · ( t C 3 k ) ) ,
λ k ( t ) = D 2 k · t 2 + D 1 k · t + D 0 k ,
where t represents the independent variable of time, and  C i k  and  D i k  are constants specific to the material for each loading/unloading cycle denoted by k. The proposed fitting curve can be characterized by a tangent at point  t 0 = C 3 k , as follows:
tan ( α ) = L k ( t ) t = C 3 k = C 1 k · C 2 k 4
where angle  α  is tilt angle relative to the x-axis.
We investigated the mechanical properties of vertebrae in two types of experiments. The first type of experiment was initiated by subjecting the vertebrae to cyclic loads using a Mecmesin MultiTest 2.5-i micro-compression machine (Mecmesin Limited, Slinfold, UK), which resulted in the formation of hysteresis loops, which were observed after ten compression cycles, as shown in Figure 3a. In the second type of experiment, we observed vertebrae under a load value, eventually cracking the vertebrae using a 2055 P-5 universal tensile testing machine, as shown in Figure 3b. We then analyzed these loops to identify their minimum and maximum values. Subsequently, we extracted the peak values for stress and stretch over time. Utilizing these peaks, we were able to characterize each distinct compression–relaxation cycle, necessitating individual approximations for each (see Figure 4a,b). Subsequently, we approximated the time indexes of the peak stress and stretch values using a modified logistic function with four parameters and a polynomial with three parameters for stress and stretch, respectively. To achieve this, we employed numerical optimization using the SciPy v 1.10.1 optimize curve_fit function, as explained in [19].
In living organisms, the hysteresis loop varies because of the compensatory and decompensatory mechanisms at play. During the first cycle (as shown in Figure 3a), the vertebrae adjusted to the applied load. Clear hysteresis loops consistent with our proposed approximations were observed in the subsequent second to tenth cycles, as is described in the next section.

3. Results and Discussion

All curve-fitting coefficients  C i k  and  D i k  for vertebra loading and unloading in cyclic load experiments and in vertebra load until fracture experiments are collected into the coefficients data file. This file, with 491 rows of coefficient data, can be found in the Supplementary Material [20]. All experimental and fitted stress–stretch curve images (Figures S1–S128) can be downloaded from the same open data cloud server. Some of these images, in particular Figure 5, Figure 6 and Figure 7, represent comparisons of stress–stretch relations for experimental and fitted curves. Figure 5 presents cycle 1 and vertebra L3’s loading and unloading stress–stretch relations; Figure 6 presents cycle 1 and vertebra L4’s loading and unloading stress–stretch relations; Figure 7 presents cycle 1 and vertebra L5’s loading and unloading stress–stretch relations. Examples of the fitted curves’ coefficients and their  R 2 ’s (R squared) for the above-mentioned figures can also be found in Table 1 and Table 2. The resulting coefficient of determination for most curves was greater than 0.999, indicating that an appropriate fitting function was selected for the description of the experimental stress–stretch curves. The number of experimental cases in Table 1 and Table 2 is different because the fitting in the first two experiments was optimized for the stretch dependence time parameter t, whereas the third experiment included only one loading cycle, for which we used stretch as the independent parameter in the logistic curve-fitting algorithm. In general, this approach is applicable when the force cell moves at a constant velocity. For all experimentally observed cases, this is true with one exception: vertebra L5 in experiment 1.
In Table 1, the columns are as follows: the first column denotes the vertebra, which was used as the experiment name; the second column denotes the experiment number; the next four columns show the values of the fitted coefficients  C i ; the  R 2  column shows the ‘R squared’ values for each fitted curve; and the final column indicates the experiment type, i.e., ’load’ or ’unload’. Coefficients  C i  can be interpreted as follows:  C 0  is the shifting value of the logistic function in the direction of the stress axis;  C 1  is a theoretical limit of the ultimate strength;  C 2  is the exponential increasing/decreasing scale ratio; and  C 3  denotes the middle of the logistic curve (i.e., the maximum  α  value).
The  D 1  coefficient values for vertebrae L3, L4, and L5 for the first loading–unloading cycle in two different experiments are presented in Table 2. All coefficient values can be found in Tables S1–S3. As shown in the results tables (Table 2 and Table S3), the stretch in most cases was linear and a function of time, with the exception of vertebra L5 in experiment 1. In all other cases, the proposed logistic function can be expressed as a logistic function of stretch with transformed linear coefficients  C 2 , C 3 .  as functions of  C i  and  D i .
Finally, the fitting function behavior of the system can be considered. For a sufficiently small t, logistic functions can be approximated as the exponential function  σ k ( t ) C 0 k + C 1 k e C 2 k ( t C 3 k ) . This behavior exhibits good agreement with the findings of other studies concerning the stress–stretch exponential behavior of elongating stretches [21,22]. For near-maximum values of t, the logistic function flattens and reaches its constant threshold values  C 0 k + C 1 k ; in mechanical terms, this indicates that the ultimate strength of the biological material has been achieved. On the other hand, the stretch parametric curves in most measured cases were linear. Therefore, in most cases, after the rearrangement of the stretch–time relation with t a function of stretch and the insertion of this new function into the logistic function, we obtain a stress–stretch logistic function. This was demonstrated in the third experimental stress–stretch fitting curve, in which good R squared values were obtained.

4. Conclusions

Investigations of three different canine vertebrae, L1, …, L7, under an external load have allowed experimental stress–stretch curves to be obtained. The experimental stress–stretch dependence of canine vertebrae was modeled as a logistic function. To achieve better fitting results, the logistic function was modified using additional constant parameters. The resulting coefficient of determination,  R 2 , for almost all curves was near to or greater than 0.999, indicating that an appropriate fitting function was selected for the description of experimental stress–stretch curves. The proposed mathematical model accurately describes the investigated experimental stress–stretch dependence over the whole range of the experiment, from the beginning of loading until the ultimate strength values were reached or the vertebra fractured.

Supplementary Materials

The following supporting information can be downloaded at: https://doi.org/10.18279/MIDAS.MathCurves.243873, accessed on 3 May 2024. Figures S1–S128: Comparisons of all stress–stretch relations for experimental and fitted curves; Tables S1–S3: Tables of all  C i  and  D i  coefficient values.

Author Contributions

Conceptualization, E.K. and A.M.; methodology, A.M.; software, E.K.; validation, E.K., N.V. and R.S.; formal analysis, E.K.; investigation, A.M.; resources, J.Š.; data curation, A.M., R.S. and N.V.; writing—original draft preparation, E.K. and A.M.; writing—review and editing, N.V., R.S. and J.Š.; visualization, E.K. and A.M.; supervision, A.M.; project administration, A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

All bioethical and animal welfare standards adhered to Regulation (EC) No. 1069/2009 of the European Parliament and of the Council of 21 October 2009 and Regulation (EC) No. 1774/2002 (animal byproducts regulation), 17 articles. This study did not require ethical approval.

Informed Consent Statement

Not applicable.

Data Availability Statement

The following experimental data files can be downloaded https://doi.org/10.18279/MIDAS.MathCurves.243873, accessed on 3 May 2024.

Conflicts of Interest

Jakov Šengaut is employed by Jakov’s Veterinary Centre. The authors declare no conflicts of interest.

References

  1. Wiczenbach, T.; Pachocki, L.; Daszkiewicz, K.; Łuczkiewicz, P.; Witkowski, W. Development and validation of lumbar spine finite element model. PeerJ. 2023, 11, e15805. [Google Scholar] [CrossRef] [PubMed]
  2. Imai, K.; Ohnishi, I.; Bessho, M.; Nakamura, K. Nonlinear finite element model predicts vertebral bone strength and fracture site. Spine 2006, 31, 1789–1794. [Google Scholar] [CrossRef] [PubMed]
  3. Matsuura, Y.; Giambini, H.; Ogawa, Y.; Fang, Z.; Thoreson, A.; Jaszemski, Y.; Lu, L.; An, K. Specimen-specific nonlinear finite element modeling to predict vertebrae fracture loads after vertebroplasty. Spine 2014, 39, E1291–E1296. [Google Scholar] [CrossRef] [PubMed]
  4. Chen, C.; Cheng, C.; Liu, C.; Lo, W. Stress analysis of the disc adjacent to interbody fusion in lumbar spine. Med. Eng. Phys. 2001, 23, 485–493. [Google Scholar] [CrossRef] [PubMed]
  5. Schmidt, H.; Heuer, F.; Drumm, J.; Klezl, Z.; Claes, L.; Wilke, H. Application of a calibration method provides more realistic results for a finite element model of a lumbar spinal segment. Clin. Biomech. 2007, 22, 377–384. [Google Scholar] [CrossRef] [PubMed]
  6. Zhong, Z.; Chen, S.; Hung, C. Load- and displacement-controlled finite element analyses on fusion and non-fusion spinal implants. Proc. Inst. Mech. Eng. Part H 2009, 223, 143–157. [Google Scholar] [CrossRef] [PubMed]
  7. Hight, T.K.; Brandeau, J.F. Mathematical modeling of the stress strain-strain rate behavior of bone using the Ramberg-Osgood equation. J. Biomech. 1983, 16, 445–450. [Google Scholar] [CrossRef] [PubMed]
  8. Schmidt, H.; Bashkuev, M.; Dreischarf, M.; Rohlmann, A.; Duda, G.; Wilke, H.; Shirazi-Adl, A. Computational biomechanics of a lumbar motion segment in pure and combined shear loads. J. Biomech. 2013, 46, 2513–2521. [Google Scholar] [CrossRef] [PubMed]
  9. Ogden, R.W. Non-Linear Elastic Deformations; Courier Corporation: Chelmsford, MA, USA, 2013. [Google Scholar]
  10. Adibhatla, S.; Satyadevi, A.; Naidu, N.S. Applicability of Hyper Elastic Models for the Analysis of Femur Bone. Int. J. Recent Technol. Eng. 2019, 8, 802–805. [Google Scholar]
  11. Arrunda, E.M.; Boyce, M.C. A three-dimensional constitutive model for the large strech behavior of rubber elastic materials. J. Mech. Phys. Solids 1993, 41, 389–412. [Google Scholar] [CrossRef]
  12. Yeoh, O.H. Some forms of the strain energy function for rubber. Rubber Chem. Technol. 1993, 66, 754–771. [Google Scholar] [CrossRef]
  13. Kostenko, E.; Stonkus, R.; Šengaut, J.; Maknickas, A. Empirical case report of the mechanical properties of three spayed canine lumbar vertebrae. Open Vet J. 2022, 12, 414–415. [Google Scholar] [CrossRef] [PubMed]
  14. Petuchova, A.; Maknickas, A.; Kostenko, E.; Stonkus, R. Experimental and theoretical investigation of aortic wall tissue in tensile tests. Technol. Health Care 2023, 31, 2411–2421. [Google Scholar] [CrossRef] [PubMed]
  15. King, A.; Eckersley, R. Statistics for Biomedical Engineers and Scientists, kindle ed.; Elsevier Science: Amsterdam, The Netherlands, 2019. [Google Scholar]
  16. Ku, H.H. Notes on the use of propagation of error formulas. J. Res. Natl. Bur. Stand. 1966, 70C, 262. [Google Scholar] [CrossRef]
  17. Thomas, G.B.; Finney, R.L. Calculus and Analytic Geometry, 5th ed.; Addison-Wesley: Boston, MA, USA, 1979. [Google Scholar]
  18. Lapshin, R. An improved parametric model for hysteresis loop approximation. Rev. Sci. Instrum. 2020, 91, 065106. [Google Scholar] [CrossRef] [PubMed]
  19. Virtanen, P.; Gommers, R.; Oliphant, T.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python. Nat. Methods 2020, 17, 261–272. [Google Scholar] [CrossRef] [PubMed]
  20. Kostenko, E.; Stonkus, R.; Šengaut, J.; Višniakov, N.; Maknickas, A. Mathematical Modeling of Stress-Strain Curves in Canine Lumbar Vertebrae, Supplemental Data; The National Open Access Research Data Archive (MIDAS): Vilnius, Lithuania, 2024. [Google Scholar] [CrossRef]
  21. Fung, Y. Elasticity of soft tissues in simple elongation. Am. J. Physiol. 1967, 213, 1532–1544. [Google Scholar] [CrossRef] [PubMed]
  22. Fung, Y.-C. Biomechanics: Mechanical Properties of Living Tissues; Springer Science & Business Media: Berlin, Germany, 2013. [Google Scholar]
Figure 1. Canine vertebrae after vertebrae separation.
Figure 1. Canine vertebrae after vertebrae separation.
Bioengineering 11 00516 g001
Figure 2. Mecmesin MultiTest 2.5-i micro-compression machine.
Figure 2. Mecmesin MultiTest 2.5-i micro-compression machine.
Bioengineering 11 00516 g002
Figure 3. (a) Minima (green x) and maxima (orange x) of experimental stress values of a cyclic load; (b) minima and maxima of experimental stress values of ultimate load.
Figure 3. (a) Minima (green x) and maxima (orange x) of experimental stress values of a cyclic load; (b) minima and maxima of experimental stress values of ultimate load.
Bioengineering 11 00516 g003
Figure 4. (a) Stress approximations by logistic function for each cycle; (b) stretch approximations by second-order polynomials for each cycle.
Figure 4. (a) Stress approximations by logistic function for each cycle; (b) stretch approximations by second-order polynomials for each cycle.
Bioengineering 11 00516 g004
Figure 5. Comparison between the stress–stretch relationships for experimental and fitted curves: (a) cycle 1, vertebra L3, experiment 1, (b) cycle 1, vertebra L3, experiment 2, (c) cycle 1, vertebra L3, experiment 3.
Figure 5. Comparison between the stress–stretch relationships for experimental and fitted curves: (a) cycle 1, vertebra L3, experiment 1, (b) cycle 1, vertebra L3, experiment 2, (c) cycle 1, vertebra L3, experiment 3.
Bioengineering 11 00516 g005
Figure 6. Comparison between the stress–stretch relationships for experimental and fitted curves: (a) cycle 1, vertebra L4, experiment 1, (b) cycle 1, vertebra L4, experiment 2, (c) cycle 1, vertebra L4, experiment 3.
Figure 6. Comparison between the stress–stretch relationships for experimental and fitted curves: (a) cycle 1, vertebra L4, experiment 1, (b) cycle 1, vertebra L4, experiment 2, (c) cycle 1, vertebra L4, experiment 3.
Bioengineering 11 00516 g006
Figure 7. Comparison between the stress–stretch relationships for experimental and fitted curves: (a) cycle 1, vertebra L5, experiment 1, (b) cycle 1, vertebra L5, experiment 2, (c) cycle 1, vertebra L5, experiment 3.
Figure 7. Comparison between the stress–stretch relationships for experimental and fitted curves: (a) cycle 1, vertebra L5, experiment 1, (b) cycle 1, vertebra L5, experiment 2, (c) cycle 1, vertebra L5, experiment 3.
Bioengineering 11 00516 g007
Table 1. C i  coefficient values for vertebrae L3, L4, and L5 during the first loading and unloading cycle in three different experiments.
Table 1. C i  coefficient values for vertebrae L3, L4, and L5 during the first loading and unloading cycle in three different experiments.
VertebraExp. C 0 C 1 C 2 C 3 R 2 Stress
L310.11337.1631.46993.23920.999load
−0.0477.62143.16594.28210.9998unload
20.887510.14380.88242.99980.9934load
0.00576.94514.01483.66380.9993unload
3−2.488313.20422.1660.9320.9992load
-----unload
L410.04677.141826.67210.22540.9996load
−0.05187.640634.37040.29860.9998unload
20.01076.02523.23382.60630.9989load
−0.20397.73963.57973.43330.9995unload
3−0.779324.65117.54730.81890.9966load
-----unload
L510.03566.3102242.63060.02920.9995load
−0.05665.7441296.32360.04530.9994unload
20.21387.32331.7942.16830.9996load
−0.05847.05223.8273.10420.9994unload
3−4.377933.419.11670.78860.9996load
-----unload
Table 2. D i  coefficient values for vertebrae L3, L4, and L5 during the first loading and unloading cycle in two different experiments.
Table 2. D i  coefficient values for vertebrae L3, L4, and L5 during the first loading and unloading cycle in two different experiments.
VertebraExp. D 2 D 1 D 0 R 2 Direction
L31−0.00002−0.45441.000160.99999down
0.000550.0399990.648540.99993up
2−0.00006−0.049811.000120.99999down
0.001490.036950.688870.99941up
L410.02194−0.459411.00090.99988down
0.032570.430830.759640.99998up
20.00041−0.051051.000350.99991down
−0.002010.064620.657120.99907up
L51−21.63102−3.046091.008230.99321down
42.69855−0.831420.846510.99158up
20.00013−0.050291.00020.99999down
0.000510.046240.727080.99998up
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kostenko, E.; Stonkus, R.; Šengaut, J.; Višniakov, N.; Maknickas, A. Quantitative Analysis of Stress–Stretch Curves in Canine Lumbar Vertebrae Using Modified Logistic Functions. Bioengineering 2024, 11, 516. https://doi.org/10.3390/bioengineering11050516

AMA Style

Kostenko E, Stonkus R, Šengaut J, Višniakov N, Maknickas A. Quantitative Analysis of Stress–Stretch Curves in Canine Lumbar Vertebrae Using Modified Logistic Functions. Bioengineering. 2024; 11(5):516. https://doi.org/10.3390/bioengineering11050516

Chicago/Turabian Style

Kostenko, Ernest, Rimantas Stonkus, Jakov Šengaut, Nikolaj Višniakov, and Algirdas Maknickas. 2024. "Quantitative Analysis of Stress–Stretch Curves in Canine Lumbar Vertebrae Using Modified Logistic Functions" Bioengineering 11, no. 5: 516. https://doi.org/10.3390/bioengineering11050516

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop