- Open Access
On the protective capacity of a safety vest for the thoracic injury caused by falling down
© The Author(s) 2019
- Received: 12 November 2018
- Accepted: 14 March 2019
- Published: 2 April 2019
Aged people all over the world are prone to fall down accidentally and be injured with fracture, such as the rib fracture. To protect the elderly, the safety vest has been developed to protect them from being injured when falling down. To effectively protect the elderly, more analysis on the protective capacity of a safety vest under different situation are needed.
Herein, a finite element model based on the computed tomography CT scanning data of a Chinese old female was built, and then used to simulate the process of falling down at different velocities. Analysis and comparison were done on the maximum shear stress, kinetic energy curves and internal energy curves with and without safety vest. The maximum shear stress indicated that the Abbreviated Injury Scale (AIS) 2+ injury risks of rib were 8%, 100% and 100% at the velocities of 1.5 m/s, 2.0 m/s and 2.5 m/s, respectively. The corresponding risks were lowered to 0%, 0% and 60% by the vest, respectively. Furthermore, the vest could absorb the internal energy resulted by the deformation of the thoracic osseous tissue by about 20%, thus decreasing the shear stress and the injury risk.
It is concluded that the safety vest decreases the injury risk when the elderly fall down, thus protects them from being injured.
- Finite element method
- Safety vest
Population aging has been a common problem all over the world. In 2008, the portion of the population older than 65 in the world was about 13%, and by 2056, it is predicted to be 24–25% . By the end of the 21 century, in every three people in the western world there would be one older than 65 . In China, the portion of the elderly has increased sustainably from 7.0% in 2000 to 10.8% in 2016 . Because of the functionality degeneration of the body, the elderly are prone to fall down, which is defined as “unexpected events in which the participant unintentionally comes to rest on the ground, floor, or lower level” . Campbell  found that about one-third of the elderly population older than 65 fall each year, and one half of those older than 85. During all the circumstances where the elderly fall down, falling down indoor occurs the most commonly . Based on the fact that women are easier to fall down in the indoor situation, the women population are more common to become the objective when researching the topic of the falling down of the elderly [5, 7–9, 11]. Besides, occurrences of falling down of women are more than those of men for all age groups , and a research on the fracture reported that the incidence of fracture is greater among women .
Since falling down in the elderly usually results in bad consequence, they are afraid and restricted to do exercise, lowering their life of quality. More seriously, they might suffer from injury or death by falling down accidently [12–14]. Hayes  and Elkington  reported that injury caused by falling down is the most common and most serious medical care problem for the elderly. According to Rubenstein , about 10% of the elderly falling down results in bone fracture, among which the rib fracture is a very common situation.
There have been several researches working on developing devices to protect the elderly from being injured when falling down [18–21]. Almost all of these devices are protective airbags with a system to detect the falling action through the sudden variation of the accelerate . When the falling action is detected, the airbag can be quickly filled with compressed air and inflated within about 200–300 ms before the elderly fall down onto the ground. However, the current using algorithms are difficult to recognize the “falling” action precisely, and moreover, it is not convenient for the airbag to be reused. Especially, the research on protection of the thorax of the elderly when falling down has never been reported.
In the equestrian sports the riders can fall occasionally, and for the riders older than 50, the most frequent injury was rib fractures [22, 23]. To provide the riders a good protection, there were several kinds of safety vest designed for them, which has been proved to have good protective capacity in some cases . We considered that the safety vest can potentially have the capacity of protecting the rib of the elderly from being injured when falling down.
The purpose of this study was to assess the risk of the rib fracture for the Chinese elderly when falling down, and the protective capacity of the safety vest by a finite element method to better describe the structure of the Chinese population, we constructed a new finite element model of the thoracic osseous tissue based on the computed tomography (CT) scanning data of a Chinese old female, instead of using the available finite element model like the total human model for safety (THUMS). Through our finite element model, the falling process was simulated to assess the risk of the rib fracture. Lastly, through comparing the computed results with and without safety vest, its protective capacity was then analyzed.
Creation of the finite element model
The CT scanning data was based on a 65-year-old Chinese female. The slice thickness was 0.620 mm and a total of 545 slices were scanned using a Philips/CT 256 scanner as shown in Fig. 1a. No lesion signals were found in the osseous tissue. The geometry containing the ribs, the rib cartilage, the sternum, the thoracic vertebrae (from L1 to L10) was obtained as shown in Fig. 1b. Some simplifications were done for the cartilage. According to the true anatomy structure of human, the 1st–7th ribs were connected with the sternum directly through the cartilage, while for the 8th–10th ribs they were not connected directly with the cartilage. For simplicity, the cartilage connecting the 8th–10th ribs were discarded. Besides, it is difficult to recognize the CT value of the cartilage precisely, therefore we drew the cartilage according to the anatomy structure. The intervertebral disks were built during the mesh process directly.
Many attempts have been made to use hexahedral elements to mesh the thoracic osseous tissue [24–27]. However, although the hexahedral elements do have better computing accuracy and efficiency, it is difficult and quite time consuming to get a finite element model with good quality only using the hexahedral elements to mesh a complicated and irregular geometry. Furthermore, many of the complex details are not accurately expressed, especially for the human skeleton. Therefore, we preferred the tetrahedral elements, even though there might be a concern that the tetrahedral elements may produce a model with lower computing accuracy and worse computing results. Actually, this concern can now be compensated by employing the rapidly advancing computing capacity. In 1992, Cifuentes found that tetrahedral elements can produce results that are equivalent to those produced with hexahedral elements in both accuracy and CPU time .
Each part of the geometry was meshed with tetrahedral elements as the cortical bone. Some researchers used the validated finite element models of the thorax with the element size of 4–5 mm [29–32], using the hexahedral elements. Here we decreased the element size to about 2 mm. For the ribs, sternum and thoracic vertebrae, triangle shell elements of the outer surface were extracted according to the obtained tetrahedral elements as the cancellous bone. For the cartilage we did not distinguish the cortical and the cancellous part, and there was only one component. The mesh of the disks was constructed according to the mesh of the vertebrae, and the elements were organized into two components: the nucleus pulposus and the fibrous ring.
Summary of the elements
Material parameters used in the model
Elastic modulus (MPa)
Load and the boundary conditions
To evaluate the protective capacity of the vest, the processes of falling down with and without vest at different velocities were simulated and compared, with the anterior part of the thorax impacting onto the floor. According to the literature , the velocity of falling can be varied among the range of 1.2–4.79 m/s due to the height and weight. Here three different velocities, 1.5 m/s, 2.0 m/s and 2.5 m/s, were chosen to be the falling velocities. In this study the model was assumed to have only translational velocity along the falling direction, and all the other motion components, such as gravitational acceleration, angular velocity and angular acceleration, were not considered. Therefore, the velocities remained constant prior to the impact with the floor, and the distance between the thorax and floor would not affect the results. The floor was assumed to be a rigid wall when the process of falling down within the indoor circumstance was simulated. The nodes in the boundary of vertebra and disc, sternum and rib were equivalent, so any boundary conditions were not needed. The tied-node-to-surface contact type was used to simulate the joints between ribs and their relating vertebrae. A contact was automatically built when we set a rigid wall and all we need to do was to select all the nodes as the slave nodes. To consider the effect of the whole body during the falling process, the upper nodes of the 1st vertebra and the lower nodes of the 10th vertebra were constrained so that they had only the displacement degree of freedom in the falling direction.
Six different simulations were accomplished in total, three without vest, and three with vest. During the process of falling down, the thoracic osseous tissue impacted the rigid wall at different velocities, and then rebound.
After that, the KE began to decrease. Then, as the rebound of the system began, the KE values increased slowly. The minimum KE values were nearly zero, for instances without vest, the values were 0.194 J, 0.652 J and 0.761 J as the initial velocity increased. For instances with vest, they were 0.305 J, 0.697 J and 1.146 J. After that, as the system rebounded, the KE began to increase It is predicted that the rebound velocity would reach an approximate constant value lower than the initial value.
The change of the IE values was just the opposite. At the beginning, the IE values were zero, and they began to increase when the contact between the system and the rigid wall occurred. At the maximum stress point, they reach the maximum value and then decreased. The maximum IE values for instances without vest were 7.981 J, 17.564 J and 26.032 J. As the vest was on, they were lowered to 6.271 J, 14.244 J and 21.788 J, respectively.
During the biomechanical simulation of the human body, it is not realistic to describe material model and its parameters accurately. Therefore in simulations where fractures were possibly to occur. A better method was to assess the fracture and injury risk according to the stress distributions and deformations, rather than considering the fitted failure criterion. The maximum shear stress is an important factor when assessing the rib injury under impact [22, 31]. In this study the AIS 2+ injury risks were obtained according to the maximum shear stress, and then compared to evaluate the protective of the safety vest. Through the results, we could conclude that the maximum shear stress in the thoracic osseous tissue was lowered by wearing a safety vest when falling down. There was a low risk of rib fracture even without vest when the falling velocity was lower than 1.5 m/s. As the velocity increased, the risk increased rapidly and wearing the vest can significantly reduce the injury risk to a very low level to protect the target from being injured. However, the risk of injury remained at high level (60%) even with the safety vest when the velocity is higher than 2.5 m/s.
The high stress zones were mainly at the 1st, 4th–6th ribs where it is much easier to occur rib fractures. Li  found that fracture occurred more commonly in the 4th and 6th ribs through experiments and simulations, and the locations corresponded to the maximum curvature points. The lung and the heart would be damaged and lead to other serious injury when the 4th–6th ribs fractured. Furthermore, they compared and studied the effect of the rib cortical thickness on the fracture locations . In reality, the cortical thickness of rib varied along the rib, and different human body had different cortical thickness [41, 42]. They found that although different thickness indicated different stress and the risk possibility are different. However, the fracture locations were the same even the rib were with different thickness. Therefore, we used a constant cortical thickness in the models. Future work will be needed to find an optimized thickness.
The influence of the vest on the value of KE was very limited, but that on the IE values was relatively significant. The IE indicated the elastic deformation resulted from the impact. Thanks to the vest, the IE value increased more slowly, and the maximum IE value at the three velocities were lowered by 22.42%, 18.90% and 16.30%, respectively. The decreased part was absorbed by the vest. Because the IE value indicated the deformation extent and thus the stress level, lower IE value means lower injury risk. We noted that the curve reaches a constant value rather than 0. It is because the motions of the anterior and posterior part of the system were not coherent, and the velocities of the different part of the model will not reach 0 at the same time. Therefore, the minimum value of the KE curve was not 0. In our models it was the elastic deformation and would disappear after an enough long time. But in reality there will be plastic deformation even if no fracture occurred, so this kind of energy difference will still exist and never disappear.
From the curves without vest we could also see that at a relative high velocity (2.0 m/s and 2.5 m/s), there was a sudden “jump” in a short time before the maximum or minimum value reached. This might be the vibration caused by the contact between the bone and the rigid wall. This situation was not found at the low velocity of 1.5 m/s. Furthermore, the vest could be a buffer between the bone and the rigid wall and thus eliminate this vibration when it was on. In computing dynamics, it is possible to be unstable when a rigid body contact with a relative “soft” body at a high velocity. The vest was proven to be able to avoid this problem, but we still do not know whether this is valid at a higher velocity. Therefore, we think that avoiding the contact with rigid body could improve our results. Of course, the results we obtained now were still good enough under the reasonable simplifications.
In conclusion, the vest does have the capacity of decreasing the injury risk when the elderly fall down. From the stress distribution we can find that the stress is localized. To improve the protective capacity of the vest, we can increase the thickness of the locations relating to the high stress, and decrease the thickness at the low stress location to make the vest more flexible.
It should be noted that the results we got in this work could be used to predict the protective capacity of the safety vest even some simplifications have been used in our model. Further work will be done to optimize our model.
There were several limitations that would affect the final results in this study. Firstly, to simulate the falling down process of a person, there were many factors should be considered, such as the body weight, accelerate and the angular accelerate caused by the gravity. In this study we simplified the models to have only the translation degree of freedom and a constant initial velocity, the other factors mentioned were not considered. In our future work we will consider these factors in our models progressively. Secondly, our models were constructed according to one single Chinese old female. Therefore, to be more accurate, more comparisons based on different objectives containing different properties are needed. Thirdly, during the falling process the organs and other soft tissues within the torso can also serve as the buffer, and therefore decrease the stress on the rib. We planned to add these soft tissue to improve the models. Lastly, the tied contact between the rib and the vertebra was just the approximate of the truth, leading to the virtual high stress at this location. Using better method to describe the joint will improve the result.
In this study we constructed the finite element model based on the CT scanning data of a Chinese old female, and then used this model to simulate the process of falling down at different velocities. Analysis and comparisons were done between the results with and without a safety vest. The maximum shear stress indicated that the AIS2+ injury risks of rib were 8%, 100% and 100% at the velocities of 1.5 m/s, 2.0 m/s and 2.5 m/s. Those risks were lowered to 0%, 0% and 60%, respectively by the protection from a safety vest. Furthermore, the vest could absorb the internal energy resulted from the deformation of the thoracic osseous tissue by about 20%, thus reducing the shear stress and the injury risk. It is believed that the safety vest could decrease the injury risk when the elderly fall down, and thus protect them from being injured. The result provided a way of studying how to protect the Chinese elderly during the process of falling down.
JL and DC initiated the research. JL and HL performed the modelling and calculations. XT supervised the project. All authors discussed and prepared the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
All data analyzed during this study are included in this published article.
Consent for publication
Ethics approval and consent to participate
This work was supported by the National Natural Science Foundation of China (Grant No. 81471743).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Day L, Finch CF, Hill KD, Haines TP, Clemson L, Thomas M, et al. A protocol for evidence-based targeting and evaluation of statewide strategies for preventing falls among community-dwelling older people in Victoria, Australia. Injury Prev. 2011;17:e3.View ArticleGoogle Scholar
- Lamb SE, Jorstad-Stein EC, Hauer K, et al. Development of a common outcome data set for fall injury prevention trials: the prevention of falls network Europe consensus. J Am Geriatr Soc. 2005;53(9):1618–22.View ArticleGoogle Scholar
- National Bureau of Statistics of China. China statistical yearbook. Beijing: China Statistics Press; 2017.Google Scholar
- Skelton DA, Kennedy J, Rutherford OM. Explosive power and asymmetry in leg muscle function in frequent fallers and no fallers aged over 65. Age Aging. 2002;31(2):119–25.View ArticleGoogle Scholar
- Campbell AJ, Robertson MC, Gardner MM, Norton RN, Buchner DM. Falls prevention over 2 years: a randomized controlled trial in women 80 years and older. Age Ageing. 1999;28:513–8.View ArticleGoogle Scholar
- Santanu M, Amit R, Subrata P. Effects of trochanteric soft tissue thickness and hip impact velocity on hip fracture in sideway fall through 3D finite element simulations. J Biomech. 2008;41:2834–42.View ArticleGoogle Scholar
- Patil R, Uusi-Rasi K, Tokola K, Karinkanta S, Kannus P, Sievänen H. Effects of a multimodal exercise program on physical function, falls, and injuries in older women: a 2-year community-based, randomized controlled trial. J Am Geriatr Soc. 2015;63(7):1306–13.View ArticleGoogle Scholar
- Tiedemann A. A balance retraining program reduces injurious falls and improves physical function in older women at risk of falling. J Physiother. 2016;62(1):48.View ArticleGoogle Scholar
- Sengul Aycicek G, Arik G, Canbaz B, Kara O, Sumer F, Ulger Z. Effect of multimodal exercise program on physical function, falls, and injuries in older women. J Am Geriatr Soc. 2016;64(2):458.View ArticleGoogle Scholar
- World Health Organization. WHO global report on falls prevention in older age. 2007. https://www.who.int/ageing/publications/Falls_prevention7March.pdf.
- Fransen M. Excess mortality or institutionalization after hip fracture: men are greater risk than women. J Am Geriatr Soc. 2002;50(4):685–90.View ArticleGoogle Scholar
- Delbaere K, Close J, Brodaty H, et al. Determinants of disparities between perceived and physiological risk of falling among elderly people: Cohort study. BMJ. 2010;341:c4165.View ArticleGoogle Scholar
- Austin N, Devine A, Dick I, et al. Fear of falling in older women: a longitudinal study of incidence, persistence, and predictors. J Am Geriatr Soc. 2007;55:1598–603.View ArticleGoogle Scholar
- Department of Health. Falls and fractures: effective interventions in health and social care. London: Department of Health; 2009.Google Scholar
- Hayes WC, Myers ER, Robinovitch SN, Van Den Kroonenberg A, Courtney AC, McMahon TA. Etiology and prevention of age-related hip fracture. Bone. 1996;18(1 Suppl):77–86.View ArticleGoogle Scholar
- Elkington J. The big falls issue. New South Wales Public Health Bull. 2002;13:2.View ArticleGoogle Scholar
- Rubenstein LZ, Josephson KR. The epidemiology of falls and syncope. Clin Geriatr Med. 2002;18:141–58.View ArticleGoogle Scholar
- Shi GY, Chan CS, Li WJ. Mobile human airbag system for fall protection using MEMS sensors and embedded SVM classifier. IEEE Sensor J. 2009;5(9):495–503.View ArticleGoogle Scholar
- Tamura T, Yoshimura T, Sekine M. A wearable airbag to prevent fall injuries. IEEE Trans Inf Technol Biomed. 2009;13(6):910–4.View ArticleGoogle Scholar
- Fukuda Y, Feng MQ, Zheng J. Rapid detection of IED light emission patterns for activating a head protective system. IEEE Sens J. 2016;16(2):498–505.View ArticleGoogle Scholar
- Majumder S, Roychowdhury A, Pal S. Effects of trochanteric soft tissue thickness and hip impact velocity on hip fracture in sideways fall through 3D finite element simulations. J Biomech. 2008;41(13):2834–42.View ArticleGoogle Scholar
- Brolin K, Jacob W. Explicit finite element methods for equestrian applications. Proc Eng. 2016;147:275–80.View ArticleGoogle Scholar
- Bilaniuk JW, Adams JM, DiFazio LT, Siegel BK, Allegra JR, Lujan JJ, Durling-Grover R, Pawar J, Rolandelli RH, Nemeth ZH. Equestrian trauma: injury patterns vary among age groups. Am Surg. 2014;80:396–402.Google Scholar
- Shin J, Untaroiu C, Lessley D. Thoracic response to shoulder belt loading: investigation of chest stiffness and longitudinal strain pattern of ribs. Sae Tech Papers. 2009;57:2775–81.Google Scholar
- Lauren WZ, Cong C, Jingwen H, Matthew P, Jonathan R. The influence of pre-existing rib fractures on global human body models consortium thorax response in frontal and oblique impact. J Biomech. 2018;69:54–63.View ArticleGoogle Scholar
- Jack CR, James VO, Emily EW. Modeling the effect of non penetrating ballistic impact as a means of detecting behind-armor blunt trauma. J Trauma. 2005;58:1241–51.View ArticleGoogle Scholar
- Yuqing N, Weixin S, James HS. Finite element models of rib as an inhomogeneous beam structure under high-speed impacts. Med Eng Phys. 2007;29:788–98.View ArticleGoogle Scholar
- Cifuentes AO, Kalbag AA. Performance study of tetrahedral and hexahedral elements in 3-D finite element structural analysis. Finite Elements Anal Des. 1992;12:313–8.View ArticleGoogle Scholar
- Belytschko T, Lin JI, Tsay C-S. Explicit algorithms for the nonlinear dynamics of shells. Comput Methods Appl Mech Eng. 1984;42(2):225–51.View ArticleGoogle Scholar
- Belytschko T, Bindeman LP. Assumed strain stabilization of the eight node hexahedral element. Comput Methods Appl Mech Eng. 1993;105(2):225–60.View ArticleGoogle Scholar
- Mendoza-Vazquez M, Davidsson J, Brolin K. Construction and evaluation of thoracic injury risk curves for a finite element human body model in frontal car crashes. Accid Anal Prev. 2015;85:73–82.View ArticleGoogle Scholar
- Mendoza-Vazquez M, Brolin K, Davidsson J, Wismans J. Human rib response to different restraint systems in frontal impacts: a study using a human body model. Int J Crashworthiness. 2013;18(5):516–29.View ArticleGoogle Scholar
- Kroell CK, Schneider DC. Impact tolerance and response of the human thorax. Stapp Car Crash Journal. 1971;15:84–134.Google Scholar
- Yang KH, Hu J, White NA. Development of numerical models for injury biomechanics research: a review of 50 years of publications in the stapp car crash conference. Stapp Car Crash J. 2006;50:429–90.Google Scholar
- Ruan J, El-Jawahri R, Chai L. Prediction and analysis of human thoracic impact responses and injuries in cadaver impacts using a full human body finite element model. Stapp Car Crash J. 2003;47:299–321.Google Scholar
- Kimpara H, Lee JB, Yang KH. Development of a three-dimensional finite element chest model for the 5th percentile female. Stapp Car Crash J. 2005;49:251–69.Google Scholar
- Ls-dyna examples of foam material modelling. Dynaexamples.com/sph.foam.Google Scholar
- Material property data for nylon 6. unreinforced. MatWeb.com.Google Scholar
- Robinovitch SN, McMahon TA, Hayes WC. Force attenuation in trochanteric soft tissues during impact from a fall. J Orthop Res. 1995;13:956–62.View ArticleGoogle Scholar
- Li ZP, Matthew WK, Damien S, Richard WK. Influence of mesh density, cortical thickness and material properties on human rib fracture prediction. Med Eng Phys. 2010;32:998–1008.View ArticleGoogle Scholar
- Haug E, Choi HY, Robin S, Beaugonim M. Handbook of numerical analysis. Computational models for the human body-human models for crash and impact simulation, vol. 12. New York: Elsevier; 2004. p. 297–361.Google Scholar
- Song E, Trosseille X, Baudrit P. Evaluation of thoracic deflection as injury criterion for side impact using a finite elements thorax model. Stapp Car Crash J. 2009;53:155–91.Google Scholar