據統計有超過八成的成年人發生過背痛，其原因大部分來自腰部的病變，在醫學臨床上，有一種常見的疾病為脊椎解離症。脊椎解離症約有90%發生在第五腰椎。當第五節腰椎椎弓產生斷裂一開始為脊椎解離，其使患者移動性及活動性降低，進而發展為脊椎滑脫，產生脊椎體的不穩定。在我們人體脊柱體發生脊椎滑脫後，會使得L5~S1椎間盤及椎間關節承受到較大的力，加上身體的重量和脊椎的活動，長期下來會造成椎間盤含水量降低，使椎間盤高度下降，產生椎間退化的情形。
是故配合醫學實務上之需求，本研究的目的就在於針對在不同腰薦角下、不同滑脫等級(採用Meyerding 分類方式)，探討在不同運動情形下的椎間盤最大應力、椎間轉角及應變能之分佈，以期在脊椎病變治療上針對不同病人有更正確的治療方向，以利醫學研究。
而在有限元素模型的建立方面則依照本研究團隊先前所使用的方法，由電腦斷層掃描(CT)取得DICOM檔，接著由醫學影像軟體3DDoctor處理堆疊成STL檔，再交由MSC.Patran有限前處理軟體建立實體模型，最後由ABAQUS有限元素軟體求解。

According to statistics, backache occurr in more then 80% adults and most of the lesions come from the waist. In clinic, a common spinal disease is spondylolysis. About 90% of spondylolysis occurs in the fifth lumbar vertebra. The fifth lumbar vertebra (L5) pars starting fracture is called spondylolysis. That causes patients with reduced mobility and activity and to develop spondylolisthesis, resulting in vertebral instability. Spondylolisthesis will allow L5 ~ S1 disc and the intervertebral joints (facet joint) to bear a larger force under the body weight and the activities of the spine. Over the time, it will cause the water content of the disc to decrease, so that the disc height decreases, resulting in the case of intervertebral disk degeneration.
Therefore, coordinating with the demand of medical practice, the purpose of this study is targeted at different lumbosacral angle, different Meyerding levels, the effects of different exercise conditions the maximum stress under the disc, intervertebral angle and the strain energy distribution, with a view to the treatment of spinal disease patients for different treatment of the right direction more to facilitate medical research.
In the finite element model of side research team in accordance with the method used previously by computed tomography (CT) to achieve DICOM file, followed by a stack of medical image processing software 3DDoctor into STL files are then referred to MSC.Patran pre-processing software to create solid model, and finally solved by the ABAQUS finite element program.

參考文獻
[1] 謝牧鄰，『以CT斷層掃描影像為基礎之脊椎有限元素分析』，國立成功大學土木工程研究所，碩士論文，中華民國九十四年六月。
[2] 林柏君，『電腦輔助腰椎之有限元素分析』，國立成功大學土木工程研究所，碩士論文，中華民國九十五年六月。
[3] Chen, C.-H., “A Finite Element study of the Biomechanical Behavior of the Nonlinear Ligamentous Thoracic and Lumbar Spine”, National Cheng Kung University, Department of Civil Engineering, Dissertation of Master Degree, June 2007.
[4] 林冠瑋,『電腦輔助脊椎之有限元素分析』，國立成功大學土木工程研究所，碩士論文，中華民國九十七年六月。
[5] Nordin, M., Frankel, V. H., ”Basic Biomechanics of the Musculoskeletal System. Lippincott Williams and Wilkins”, Baltimore, 2001
[6] Huang, K.-Y., Lee, C.-H. and Chen, P.-Q. “Basic Research and Future Treatment of Intervertebral Disc Disorder”, Formosan Journal of Medicine, Vol. 7, No. 5, 2003
[7] 劉哲榮,『後方腰椎間融合手術後應力重新分配之有限元素分析』，國立成功大學土木工程研究所，碩士論文，中華民國九十八年六月。
[8] 王修卓、彭完全，『腰椎內固定術後對鄰近節段的影響』，中國醫師進修雜誌，2007年，第30卷，第17期。
[9] Lee, C.-K., “Accelerated degeneration of the segment adjacent to a lumbar fusion”, Spine, Vol. 13, pp. 375–7, 1988
[10] Lin, R.-M., Huang, K.-Y., Lee, Y.-L. and Li, J.-D., “Factors affecting disability and physical function in degenerative lumbar spondylolisthesis of L4–5: evaluation with axially loaded MRI”, European Spine Journal, Published online: 14 June 2009.
[11] Gray Henry, “Grey's Anatomy”, New York, 1917
[12] Tittel, K. “Beschreibende und funktionelle Anatomie des Menschen”,
8. Ed. Jena, VEB Gustav-Fischer Verlag, 1990.
[13] Mow, V. and Hayes, W. C., “Basic orthopaedic biomechanics”, New York, Raven press Ltd, 1991.
[14] White III, A. A. and Panjabi, M. M., “Clinical biomechanics of the spine”, second edition, J.B. Lippincott Company, Philadelphia, 1990
[15] Andrea Strayer, ” Lumbar Spine Surgery-A Guide to Preoperative and Postoperative Patient Care”, AANN Reference Series for Clinical Practice, 4700 W. Lake Avenue, Glenview, IL 60025-1485
[16] Cook, R. D., Malkus, D. S., Plesha, M. E. and Witt, R. J., “Concepts and applications of finite element analysis”, fourth edition, John Wiley & Sons, New York, 2002
[17] Lee, C., Kim, Y.-E., Lee, C.-S., Hong, Y.-M., Jung, J. and Goel, V. K., “Impact response of the intevertebral disc in a finite element model”, Spine, Vol. 25, pp. 2431-2439, 2000.

[18] Cheung, J.T.-M., Zhang, M. and Chow, D.H.-K.,”Biomechanical responses of the intervertebral joints to static and vibrational loading: a finite element study”, Clinical Biomechanics, Vol.18, pp790-799, 2003.
[19] Baroud, G., Nemes, J., Heini, P., and Steffen, T.,”Load shift of the intervertebral disc after a vertebroplasty: a finite-element study”, European Spine Journal, Vol. 12, No. 4, pp. 421-426, 2003.
[20] Polikeit, A., Ferguson S. J., Nolte L. P. and Orr T. E.,“Factors influencing stresses in the lumbar spine after the insertion of intervertebral cages: finite element analysis”, European Spine Journal, Vol. 12, No. 4, pp. 413-420, 2003.
[21] Denozière, G.,”Numerical modeling of a ligamentous lumbar motion segment”, M.S. thesis, Department of Mechanical Engineering, Georgia Institute of Technology, Georgia, U.S.A., 2004.
[22] Sairyo, K., Goel, V. K., Masuda, A., Vishnubhotla, S, Faizan, A, Biyani, A., Ebraheim, N, Yonekura, D, Murakami, R. I. and Terai., T.
,”Three-dimensional finite element analysis of the pediatric lumbar spine”, European Spine Journal, Vol. 15, pp.923-929, 2006.
[23] Rohlmann, A., Burra, N.K., Zander, T. and Bergmann, G., “Comparison of the effects of bilateral posterior dynamic and rigid fixation devices on the loads in the lumbar spine”, European Spine Journal, 2007.
[24] Wang, J.-L., Parnianpour, M., Shirazi-Adl, A. and Engin, A. E.,
“Viscoelastic Finite-Element Analysis of a Lumbar Motion Segment in Combined Compression and Sagittal Flexion: Effect of Loading Rate”, Spine, Vol. 25, No. 3, pp. 310-318, 2000.

[25] Pitzen, T., Geisler, F. H., Matthis, D.,Müller-Storz, H., Pedersen, K. and Steudel, W.-I., “The influence of cancellous bone density on load sharing in human lumbar spine: a comparison between an intact and a surgically altered motion segment”, European Spine Journal, Vol. 10, No. 1, pp. 23-29, 2001.
[26] Goto, K., Tajima, N., Chosa, E., Totoribe, K., Kuroki, H. and Arizumi, Y., “Mechanical analysis of the lumbar vertebrae in a three-dimensional finite element method model in which intradiscal pressure in the nucleus pulposus was used to establish the model”, Journal of Orthopaedic Science, Vol. 7, No. 2, pp.243-246., 2002.
[27] Natarajan, R. N., and Andersson, G.B.J.,” Modeling the annular incision in a herniated lumbar intervertebral disc to study its effect on disc stability”, Computers Structures, Vol. 64, No. 5-6, pp.1291-1297, 1997.
[28] Ng, H.-W., Teo, E.-C., “Nonlinear Finite-Element Analysis of the Lower Cervical Spine(C4–C6) Under Axial Loading”, Journal of Spinal Disorders, Vol. 14, No. 3, pp. 201–210, 2001.
[29] Gwanseob, Shin,”Viscoelastic responses of the lumbar spine during prolonged stooping”, Ph.D. dissertation, NCSU, USA, 2005.
[30] Glema, A., Lodygowski, T., Kakol, W., Wierszycki, M. and Ogurkowska, M. B., “Modeling of intervertebral discs in the numerical analysis of spinal segment”, ECCOMAS, pp. 24-28, 2004.
[31] Adams, M. A., Hutton, W. C., “The effect of posture on the role of the apophysial joints in resisting intervertebral compressive forces”, The Journal of Bone and Joint Surgery, Vol. 62-B, No. 3, pp. 368-362, 1980.

[32] Sharma, M., Langrana, N. A. and Rodriguez, J., “Role of ligaments and facets in lumbar spinal stability”, Spine, Vol. 20, No. 8, pp. 887-900, 1995.
[33] Sharma, M., Langrana, N. A. and Rodriguez, J., “Modeling of facet articulation as a nonlinear moving contact problem: sensitivity study on lumbar facet response”, Journal of Biomechanical Engineering, Vol. 120, pp. 118-125, 1998.
[34]Chazal, J., Tanguy, A., Bourges, M., Gaurel, G., Escande, G, Guillot, M. and Vanneuville, G., ” Biomechanical properties of spinal ligaments and a histological study of the supraspinal ligament in traction”, Journal of Biomech. , Vol.18(3), pp.167-176, 1985.
[35] Shirazi-Adl, S. A., Ahmed, A. M. and Shrivastava, S. C., “A finite element study of a lumbar motion segment subjected to pure sagittal plane moments”, Journal of Biomech. , Vol.19, pp.331-350, 1986.
[36] Naira H.C.K., “Computational analysis of the time-dependent biomechanical behavior of the lumbar spine”, Ph.D. thesis, Ohio State University, USA, 2004.
[37] Goel, V. K., Park, H. and Kong, W.,”Investigation of vibration characteristics of ligamentous lumbar spine using the finite element approach”, Journal of Biomechanical Engineering, Vol.116, pp.377-383, 1994.
[38] Manuel, J. G. R., Leidy Y. S. G., “Comparison of hyperelastic material models in the analysis of fabrics”, International Journal of Clothing Science and Technology, Vol. 18 No. 5, pp. 314-325, 2006.
[39] ABAQUS CAE, Material Property Evaluation.

[40] Holte, D. C., O'Brien, J. P., Renton, P., “Anterior lumbar fusion using a hybrid interbody graft.A preliminary radiographic report”, European Spine Journal, Vol. 3,No. 1, pp. 32-38, 1994.
[41]Spine Fusion Surgery Video, Spine-health.com.
[42] Peter, F. U. Jr., “Bone Graft Substitutes”, Spine-health.com.
[43] Sengupta,D. K., Mehdian, SMH , Mulholland, R. C.,Webb, J. K. and Ohnmeiss,D. D., “Biomechanical evaluation of immediate stability with rectangular versus cylindrical interbody cages in stabilization of the lumbar spine”, BMC Musculoskeletal Disorders, 2002.
[44] Peter F. U. Jr.,” Lumbar Spine Surgery”, Spine-health.com.
[45] Bray, R. S., “Interplate spine fusion device, subsidence control without stress shielding”, Orthopaedic Product News, September/October 2006, pp. 22-25.
[46] MSC.Patran, Documentation.
[47] Yamamoto, I., Panjabi, M. M., Crisco, T. and Oxland, T., “Three-Dimensional Movements of the Whole Lumbar Spine and Lumbosacral Joint”, Spine, Vol. 14, No.11, pp.1256-1260.
[48] Chen, C.-S. , Cheng, C.-K., Liu, C.-L. and Lo, W.-H., “Stress analysis of the disc adjacent to inter-body fusion in lumbar”, Medical Engineering & Physics, Vol. 23, pp. 483–491, 2001.
[49] Jack Zigler, “Total Disc Replacement-ProDisc”, Spine-health.com