Document Type : Original Article


Faculty of Engineering, Department of Mechanics, Imam Khomeini International University, Qazvin, Iran


Objective: In this paper, finite element model of the L4 vertebra subjected to combination of compression and flexion loading in isotropic and anisotropic cases is investigated. Methods: In both cases, the vertebra is considered homogeneous. Also, the body of vertebra is divided to cancellous and cortical sections in anisotropic model, but the process is assumed isotropic such as isotropic model. The maximum Von Mises stress on the fourth lumbar vertebrae is obtained. Also, the stress intensity factor is analyzed with placing a small crack on the critical region of the model from view point of fracture mechanics. Furthermore, the required force for the fracture of fourth lumbar vertebrae is obtained through increasing the applied force for assumed model. Results: The results show that the highest stress value and its position is 7.237MPa in the upper pedicle region for anisotropic property of vertebrae. At the end of this article, stress intensity factors in different aspect ratios of crack for anisotropic vertebrae under combination of flexion and compression loading are plotted.


Ahmad, Z., Arifin A.K., 2010. Probabilistic stress analysis of the Human lumbar spine extended finite element method. The 11th Asia Pacific Industrial Engineering and Management Systems Conference. Melaka, 7-10 December. 1-5.
Bonfield, W., Grynpas, M.D., Young, R.J., 1978. Crack velocity and the fracture of bone. J. Biomechanics. 1, 473-479.
Chena, S.H., Zhongb, Z.C., Chenc, C.S., 2009. Biomechanical comparison between lumbar disc arthroplasty and fusion. J. Medical Engineering & Physics. 31, 244-253.
Hernandez, C.J., Beaupre, G.S., Keller, T.S., 2001. The influence of bone volume fraction and ash fraction on bone strength and modulus. J. Bone. 29, 74-78.
Herrera, A., Panisello, J.J, Ibarz, E., 2007. Long-term study of bone remodeling after femoral stem: A comparison between dexa and finite element simulation. J. Biomechanics. 40, 15-25.
Isaac, M., Graham, A., 2011. probabilistic assessment of vertebral cortical bone fracture of intraosteonal structures. MSc Thesis, University of Wright State. USA.
Kayabasi, O., Ekici, B., 2008. Probabilistic design of a newly designed cemented hip prosthesis using finite element method. J. Materials and Design. 29, 963-971.
Kuo, C.S., Hu, H.T., Lin, R.M., 2010. Biomechanical analysis of the lumbar spine on facet joint force an intradiscal pressure-a finite element study. J. BMC Musculoskeletal Disorders. 11, 1-13.
Kurtz, S.M., Edidin, A.A., 2006. Spine Technology Handbook. Elsevier.
Nabhani, F., Wake, M., 2002. Computer modeling and stress analysis of the lumbar spine. J. Materials Processing Technology. 127, 40-47.
Norman, T.L., Vashishth, D., Burr, D. B., 1994. Fracture toughness of human bone under tension. J. Biomechanics. 28, 309-320.
Odin, G., Savoldelli, C., Bouchard, P.O., 2010. Determination of young’s modulus of mandibular bone using inverse analysis. J. Medical Engineering & Physics. 32, 630-637.
Parker, A.P., 1981. The mechanics of fracture and fatigue: an introduction. Newyork.
Weisse, B., Aiyangar, A.K., Affolter, C., 2012. Determination of the translational and rotational stiffnesses of an L4–L5 functional spinal unit using a specimen-specific finite element model. J. Mechanical behavior of biomedical materials. 13, 45- 61.
White, A., Panjabi, M., 1990. Clinical biomechanics of the spine. Second ed. J.B. Lippincott Company. 1-84.
Xia, Q.T., Tan, K.W., Lee, V.S., 2006. Investigation of Thoracolumbar T12-L1 burst fracture mechanism using finite element method. J. Medical Engineering & Physics. 28, 656-664.
Zulkifli, A., Ariffin, A.K., 2012. Monte Carlo analysis of the human vertebra based on compressive loading. J. Mechanical & Mechatronics Engineering. 12, 69-73.