三维有限元分析腰骶椎结构的动态特.doc

www.CRTER.org武晓丹，等. 三维有限元分析腰骶椎结构的动态特性三维有限元分析腰骶椎结构的动态特性武晓丹1，张顺心1，范顺成1，李 晔2，贾少薇1，谢俊德1，韩 立3，4 (1河北工业大学机械学院，天津市 300130；2中国医学科学院北京协和医院，北京市 100005；3天津医科大学医学影像学院，天津市 300203；4美国密歇根大学安娜堡分校医学院，美国 48105)引用本文：武晓丹，张顺心，范顺成，李晔，贾少薇，谢俊德，韩立. 三维有限元分析腰骶椎结构的动态特性J.中国组织工程研究，2017，21(15):2388-2394.DOI:10.3969/j.issn.2095-4344.2017.15.017 ORCID: 0000-0002-2449-3124(武晓丹)文章快速阅读：椎弓根螺钉内固定后腰骶椎模型的有限元分析有限元分析：利用有限元分析软件ABAQUS对正常腰骶椎模型及椎弓根螺钉内固定后腰骶椎模型分别进行模态分析及谐响应分析。比较指标：(1)模态分析结果中的固有频率；(2)稳态分析得到的幅频特性曲线中位移最大幅值。建模：(1)正常腰骶椎模型；(2)椎弓根螺钉内固定后腰骶椎模型。武晓丹，女，1990年生，河北省石家庄市人，汉族，2017年河北工业大学毕业，硕士，主要从事脊柱的生物力学方面的研究。通讯作者：韩立，副教授，天津医科大学，天津市 300070 通讯作者：张顺心，教授，硕士生导师，河北工业大学，天津市 300130 中图分类号:R318文献标识码:B文章编号:2095-4344(2017)15-02388-07稿件接受：2017-03-15Wu Xiao-dan, Master, School of Mechanical Engineering, Hebei University of Technology, Tianjin 300130, ChinaCorresponding author: Han Li, Associate professor, School of Medical Imaging, Tianjin Medical University, Tianjin 300203, China; Medical College, University of Michigan, Ann Arbor, Michigan 48105, USACorresponding author: Zhang Shun-xin, Professor, Masters supervisor, School of Mechanical Engineering, Hebei University of Technology, Tianjin 300130, China结论：(1)固定后的腰骶椎模型在受到外界激励时的振幅比(2)正常腰骶椎模型减小，对外界激励的敏感性减弱；原因是螺钉内固定系统对腰骶椎部分起到了保护作用。文题释义：模态分析：模态是机械结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型。通过模态分析方法可以得到结构物在某一易受影响的频率范围内的各阶主要模态的特性，就可以预言结构在此频段内在外部或内部各种振源作用下产生的实际振动响应。因此，模态分析是结构动态设计及设备故障诊断的重要方法。谐响应分析：用于确定线性结构在承受随时间按正弦(简谐)规律变化的载荷时的稳态响应，分析过程中只计算结构的稳态受迫振动，不考虑激振开始时的瞬态振动，谐响应分析的目的在于计算出结构在几种频率下的响应值(通常是位移)对频率的曲线，从而使设计人员能预测结构的持续性动力特性，验证设计是否能克服共振、疲劳以及其他受迫振动引起的有害效果。摘要背景：有文献对人体正常脊柱腰骶椎进行了固有模态分析及谐响应分析，但尚未有研究涉及对比椎弓根螺钉内固定前后腰骶椎的固有模态分析及谐响应分析。目的：采用三维有限元方法探讨人体腰骶椎的动力学特性。方法：基于CT扫描图像建立并验证经椎弓根螺钉系统内固定前后腰骶椎L1-S1节段有限元模型，然后对固定前后腰骶椎模型分别进行有限元模态分析和谐响应分析。结果与结论：在腰骶椎模型的腰椎L1，L3，L5节段棘突处选取代表性节点，节点编号分别为A，B，C；内固定后腰骶椎模型的各个节点在Y，Z方向上位移最大峰值比正常腰骶椎模型位移最大峰值均明显减小，表明螺钉内固定系统对腰骶椎部分起到了保护作用，使得固定后的腰骶椎在受到外界激励时振幅减小，对外界激励的敏感性减弱；腰骶椎的模态分析是进一步进行动力学分析的基础，确定了腰骶椎的固有频率、振型和振幅等振动参数。通过谐响应分析研究简谐载荷对人体脊柱腰骶椎L1-S1节段的影响，对于腰骶椎的振动特性分析方面具有重要意义。关键词：骨科植入物；数字化骨科；有限元；生物力学；振动；腰骶椎；模型；内固定系统；ABAQUS；模态分析；动力学；谐响应主题词：腰椎；内固定器；生物力学；有限元分析；动力学；组织工程基金资助：天津市应用基础与前沿技术研究计划资助项目(13JCYBJC41200)3 P.O.Box 1200,Shenyang 110004 2391ISSN 2095-4344 CN 21-1581/R CODEN: ZLKHAHDynamic characteristics of the lumbosacral vertebrae based on three-dimensional finite element models Wu Xiao-dan1, Zhang Shun-xin1, Fan Shun-cheng1, Li Ye2, Jia Shao-wei1, Xie Jun-de1, Han Li 3, 4 (1School of Mechanical Engineering, Hebei University of Technology, Tianjin 300130, China; 2Peking Union Medical College Hospital, Chinese Academy of Medical Sciences, Beijing 100005, China; 3School of Medical Imaging, Tianjin Medical University, Tianjin 300203, China; 4Medical College, University of Michigan, Ann Arbor, Michigan 48105, USA)AbstractBACKGROUND: Inherent modal analysis and harmonic response analysis on the human normal lumbosacral vertebrae have been reported, but there is a lack of comparative research on their modal analysis results before and after pedicle screw fixation. OBJECTIVE: To explore the dynamic characteristics of human lumbosacral vertebrae using three-dimensional finite element method. METHODS: Finite element model of lumbosacral vertebrae (L1-S1) before and after pedicle screw fixation was developed and validated based on CT images, and the modal analysis and harmonic response analysis were then conducted. RESULTS AND CONCLUSION: (1) Representative nodes were selected at the spinous process segments of L1, L3 and L5, and numbered as A, B, and C, respectively. (2) The maximum displacement of each node in Y and Z directions of lumbosacral vertebral model after internal fixation was significantly decreased compared with those of the normal lumbosacral vertebral model, suggesting that screw fixation system plays a protective role in lumbosacral vertebrae, and reduces its amplitude under external load, thus diminishing its sensitivity to external load. (3) The lumbosacral vertebral modal analysis can provide basis for further study on dynamic analysis, and the parameters such as natural frequency, modal shape and vibration amplitude of the lumbar spine have been determined.Subject headings: Lumbar Vertebrae; Internal Fixators; Biomechanics; Finite Element Analysis; Kinetics; Tissue EngineeringFunding: the Application Basic and Advanced Technology Research Program of Tianjin City, No. 13JCYBJC41200Cite this article: Wu XD, Zhang SX, Fan SC, Li Y, Jia SW, Xie JD, Han L. Dynamic characteristics of the lumbosacral vertebrae based on three-dimensional finite element models. Zhongguo Zuzhi Gongcheng Yanjiu. 2017;21(15):2388-2394. 0 引言 Introduction腰椎是人体传导重力的重要机构，自L1顺序而下椎体逐渐增大，发病率逐渐增高，呈塔形分布，椎管横截面积逐渐变小，椎间盘荷载由小到大。其中以L5-S1为应力集中节段，该处椎间盘厚，活动范围较大，容易受到损伤。长期以来，振动一直被认为是影响脊柱健康的重要危险因素。有研究报道，曝于振动环境中的人产生背痛的概率比平常人普遍要高1.4-9.5倍。长期性振动载荷作用于腰椎将引起机械损伤、椎间盘退化和脊柱畸形1-4。中国是劳动密集型社会，体力劳动者数量庞大，脊椎类疾病的防治与治疗更应受到社会的重视。人体处于汽车、飞机等振动环境时，腰椎会随着外界环境的激振产生持续时间较长的衰减振动引起疾病。有关振动引起的危害和脊柱的累积损伤已开展了大量的实验研究、有限元数值模拟分析3-11。在有限元仿真方面12-20，用于振动分析的有限元模型有线性和非线性模型、单节段和多节段模型1, 21-23。单节段腰椎模型对全身振动引起的腰椎临床相关生物力学参数的预测具有局限性。此外，关于腰椎模态分析的研究多集中在受损腰椎的情况5, 24，尚未有研究涉及对比椎弓根螺钉内固定前后腰骶椎L1-S1节段的固有模态分析及谐响应分析。文章在以往的研究基础上，建立腰骶椎L1-S1节段的后路经椎弓根螺钉内固定前后的三维有限元模型25-31，利用有限元分析软件ABAQUS分别对固定前后腰骶椎模型进行模态分析和谐响应分析31-35，确定腰骶椎的固有频率、振型和振幅等振动参数，并通过谐响应分析研究简谐载荷对人体脊柱腰骶椎L1-S1节段的影响，将分析结果进行比较，这对于腰骶椎的振动特性分析等方面具有重要意义。 1 材料和方法 Materials and methods 1.1 设计 有限元模拟仿真试验。1.2 时间及地点 于2014年10月至2016年10月在河北工业大学实验室完成。1.3 材料 主要仪器：General Electrics 64层螺旋CT机。软件：医学图形处理软件MIMICS 16.0、逆向工程软件GEOMAGIC 12.0、三维建模软件SolidWorks和有限元分析软件ABAQUS6.12-3。1.4 方法1.4.1 腰骶椎后路螺钉内固定前后的模型建立 使用螺旋CT对1例30岁健康女性志愿者的脊柱进行连续的横断面扫描，扫描范围T12-S1，扫描层厚度0.625 mm；将DICOM格式的图像数据导入三维重建软件Mimics 16.0，根据已选取的CT图像中不同组织的灰度值，通过阈值划分来提取相应的组织，建立腰骶椎的三维实体模型。将模型进行表面光滑优化处理，使模型更加准确；将生成的腰骶椎三维有限元模型导入Mimics中的3-matic模块中建立与上下椎体紧密贴合的椎间盘模型。根据医学解剖结构，将椎体分为皮质骨、松质骨和椎体后部三部分，椎间盘分为上终板、下终板、髓核和纤维环四部分；然后对每一部分分别进行体网格划分，生成各部分的三维四面体单元。利用SolidWorks三维建模软件系统，模拟XIA脊柱内固定系统(美国STRYKER公司)，包括10枚椎弓根螺钉 (4.5 mm60 mm)、2枚连接棒(6.0 mm300 mm)；将钉棒系统以STL文件格式输入Mimics软件，在三维视窗下，组装腰椎和钉棒系统。按腰椎后路椎弓根内固定的手术确认骨钉和连接棒的位置。其中腰椎植入螺钉的进针点在上关节突外缘的切线与横突平分线的交点位置，螺钉向中心倾斜5-10，如图1所示。然后对装配好的腰骶椎部分和螺钉部分进行布尔减运算，使腰骶椎模型中出现安装螺钉所需要的孔道。最后在Mimics软件中对模型赋予材料属性后将其以inp格式输出到有限元分析软件Abaqus6.12-3中进行相关问题的计算。模型各组成元件的材料和性质：本研究采用的材料属性是大多数研究者公认的经过验证的数据2。正常脊柱腰骶椎各部位的材料属性如表1所示。图1 进钉位置图Figure 1 Nail position map椎体之间靠韧带连接，本研究用线性张力弹簧模拟出前纵韧带、后纵韧带、黄韧带、棘上韧带、棘间韧带、横突间韧带等主要韧带，定义弹簧刚度=弹性模量横截面积平均长度，各韧带结构的参数采用大多数文献公认的数据如表2所示36。 脊柱内固定系统的钉棒材质均为Ti6AI4VELI(钛合金)，参数采用大多数文献公认的数据37-39，其中椎弓根螺钉和连接棒的弹性模量均为104 800 MPa，泊松比均为0.33。模型的有效性验证：有限元模型有效性验证的常用方法有力-位移曲线的比较和力矩-旋转角度曲线的比较21， 40。但由于不同作者实验时所加载的力或力矩大小不同，导致实验所得位移值或角度值不一样，实验结果可比性差。因此，平均刚度值的可比性更强，用平均刚度值来验证有限元模型的有效性更科学。平均刚度值的计算方法是加载的力矩除以该力矩方向上的角位移。实验在L1-S1节段施加10 Nm的力矩，通过计算，本模型在前屈、后伸、侧屈、扭转作用下的平均刚度分别为1.11，2.25，2.18， 3.70 Nm/()。董凡等41用L1-L2、L3-L4各1个，L2-L3、L4-L5、L5-S1各2个，共8个FSU进行测试，加载10 Nm的力矩；Heth等42用14个L2-S1节段施加6 Nm的力矩；Yamamoto等43用10个全腰段标本，加载10 Nm的力矩。在前屈、侧屈和扭转力矩作用时，本模型的平均刚度小于董凡和Yamamoto的模型，而大于Heth的模型；在后伸力矩作用时，本模型的平均刚度小于Heth的模型。从每一种加载方式看，实验的刚度值基本处于中间位置，见表3。所得平均刚度值和尸体实测法所得结果基本符合。因此，本实验所建立的腰骶椎三维有限元模型是有效的。边界条件和载荷：考虑人体上身质量对腰骶椎的动态特性的影响，将一个简化的50 kg质量加到腰椎L1的上表面，约束骶骨底面上平移和转动的6个自由度。椎骨与椎间盘的上下终板之间，椎间盘各分块之间，骨钉和椎骨的钉道之间的接触关系均为绑定tie，上下关节突之间定义为面与面的有限滑动接触，摩擦系数为0.2。1.4.2 模态分析的理论基础 模态特性分析可以确定腰骶椎结构的自振频率和相应振型。脊柱具有固有频率，应尽量避免使外激励频率接近其固有频率，降低共振可能性，减小振动幅值。而各阶固有频率都有对应一个固有振型，准确计算出结构固有振型对于了解不同外界激励作用下的振型具有重要意义。目前最常用的模态分析方法包括：子空间迭代法(subspace iteration)和分块Lanczos方法(Block Lanczos)。使用何种模态提取方法主要取决于模型的大小(相对于计算机计算能力而言)和具体的应用场合。当模型的规模较大，且需要提取多阶振型时，Lanczos方法的速度更快。本文选择分块Lanczos法进行腰椎模态分析。模态计算是分析动力响应和其他动力学特性问题的基础。如果忽略阻尼的影响，则多自由度无阻尼振动系统的运动微分方程为：(1)M X+K X= Q 式(1)是系统的自由振动方程，又称为动力特性方程。其中：K为刚度矩阵；M为质量矩阵；Q为载荷矢量；X表示节点位移的列向量(它是时间t的函数)。在无外力作用，即Q=0，就得到系统的自由振动方程：(2)M X+K X=0 因为任何弹性体的自由振动都可以分解为一系列简谐振动的叠加，则方程(2)的简谐振动解为：X(t)=X0sint(3) 式中：X0为位移的X的振幅列向量，与时间t无关；为固有频率。将式(3)代入式(2)，由于sint不恒为0，故得：(4)K-2M X0 =0图2 正常腰骶椎的前3阶振型Figure 2 The first, second and third order vibration mode of the normal lumbosacral vertebrae图注：图A为第1阶振型，B为第2阶振型，C为第3阶振型。B A C图3 内固定后腰骶椎的前3阶振型Figure 3 The first, second and third order vibration mode of the lumbosacral vertebrae after internal fixation图注：图A为第1阶振型；B为第2阶振型；C第3阶振型。B A C CB A图4 节点A在X，Y，Z三个方向与频率对应的位移响应值Figure 4 Frequency and displacement response values of the node A in X, Y and Z directions图注：图A为X方向；B为Y方向；C 为Z方向。图5 节点B在X，Y，Z三个方向与频率对应的位移响应值Figure 5 Frequency and displacement response values of the node B in X, Y and Z directions图注：图A为X方向；B为Y方向；C 为Z方向。图6 节点C在X，Y，Z三个方向与频率对应的位移响应值Figure 6 Frequency and displacement response values of the node C in X, Y and Z directions图注：图A为X方向；B为Y方向；C 为Z方向。正常腰骶椎内固定后腰骶椎正常腰骶椎内固定后腰骶椎移位(mm)移位(mm)正常腰骶椎内固定后腰骶椎移位(mm)频率(Hz)频率(Hz)频率(Hz) CB频率(Hz) A正常腰骶椎内固定后腰骶椎正常腰骶椎内固定后腰骶椎正常腰骶椎内固定后腰骶椎移位(mm)移位(mm)移位(mm)频率(Hz)频率(Hz) AB C正常腰骶椎内固定后腰骶椎正常腰骶椎内固定后腰骶椎正常腰骶椎内固定后腰骶椎移位(mm)移位(mm)移位(mm)频率(Hz)频率(Hz)频率(Hz)表2 腰骶椎有限元模型的主要韧带材料属性Table 2 Material parameters of the main ligaments of the finite element models of lumbosacral vertebrae 表1 腰骶椎有限元模型各组成部分材料属性Table 1 Material parameters of the finite element models of lumbosacral vertebrae 主要韧带弹性模量(MPa)截面积(mm2)平均长度(mm)刚度(N/m)密度(g/cm3)前纵韧带7.822.4208.741.0后纵韧带107.0125.831.0黄韧带1714.11515.381.0横突间韧带100.6320.191.0棘间韧带1014.11310.851.0棘上韧带8.010.5222.391.0结构弹性模量(MPa)泊松比性质密度(g/cm3)皮质骨12 0000.3均质等向1.7松质骨1000.2均质等向1.1椎体后部3 5000.25均质等向1.4终板5000.25均质等向1.2纤维环60.4均质等向1.0髓核10.49均质等向1.02表3 本模型与其他模型平均刚度的比较 Nm/()Table 3 The average stiffness of this model compared with the other models 模型加载(Nm)前屈后伸侧屈扭转Yamamoto等43(L1-L5)101.753.222.445.26Heth等42(L2-S1)61.102.351.332.61董凡等41(FSU)102.353.582.868.98本实验(L1-S1)83.70在自由振动时，结构中各节点的振幅X0不全为0，故式(4)括号内的矩阵行列式之值必须为0，即： (5)K-2M =0 结构的刚度矩阵K和质量矩阵M都为n阶方阵，n为节点自由度数目。因此，式(5)是关于的n次代数方程。由此，可以求出结构的n个固有频率。对每一个固有频率，可以根据式(4)确定各节点的振幅值X0，称为结构的振型。1.5 主要观察指标 通过有限元软件对2种腰骶椎模型进行模拟仿真实验，比较由模态分析得到的频率及由谐响应分析得到的幅频特性曲线中位移幅值。2 结果 Results 2.1 固定前后腰骶椎有限元模态分析 对腰骶椎结构的模态分析可以得到腰骶椎的固有振动特性即固有频率及其振型，在此基础上可以进一步分析腰骶椎结构在外界激励下的响应。约束骶椎下表面的全部自由度，分别提取固定前后腰骶椎的前10阶模态，计算结果见表4。由表4可知，内固定后腰骶椎模型的前3阶固有频率要稍微大于正常腰骶椎模型的固有频率，内固定后腰骶椎模型的第4-10阶固有频率要明显小于正常腰骶椎模型的固有频率。人体脊柱腰骶椎的固有频率和振型直接影响到它在遭受外界载荷时的动态响应，为避免发生共振，应尽量使人体所处环境的外界振动频率远离腰骶椎固有频率。一些研究表明人体严重的振动响应均在30 Hz以下，对于腰骶椎模态分析结果前三阶模态具有实际意义，图2，3分别为正常腰骶椎、内固定后腰骶椎的前3阶模态振型。图中的振型变化大小只是一个相对值，反映在某一频率下的阶次正常腰骶椎固有频率(Hz)内固定后腰骶椎固有频率(Hz)振型描述11.331.57整体前后俯仰振型21.651.86整体侧向俯仰振型311.6412.43整体垂向拉伸振型476.9756.33整体横向扭转振型5111.45101.18L2、L3侧向弯曲振型6127.71103.32L1、L2、L3、L4弯曲振型7208.45112.81L1、L2、L3、L4扭转振型8254.45122.06L1、L2、L3、L4弯曲二阶振型9268.23151.29L1、L2、L3、L4横向扭转振型10337.92151.82L1、L2、L3、L4横向扭转二阶振型表4 内固定前后脊柱腰骶椎的前10阶固有频率Table 4 The 10 order natural frequency of the lumbosacral vertebrae before and after internal fixation 振动情况，并不反映实际振动位移变化。通过观察变形动画，调节位移放大比例，可以观察每一阶模态腰骶椎变形情况。由模态振型描述可知，低阶振型比较简单，叠加现象不明显，但对腰骶椎整体结构动态特性影响较大。随着频率的增加，振型越来越复杂。通过分析各阶振型，可知正常腰骶椎及内固定后腰骶椎的一阶振型均表现为整体的前后俯仰振动，L1后部结构棘突处振动最大；二阶振型均表现为整体向左弯曲振动，L1后部结构变形较大；三阶振型均表现垂向拉伸振动，L1后部变形最大；第4-10阶振型则表现为相对复杂的振动。由上可知，前3阶的最大振动均发生在L1椎体后部，表明L1椎体后部对外界振动较敏感，而且高阶振型比低阶振型复杂，这与实际的脊柱生物力学理论知识相符合。由模态分析得到的腰骶椎结构每一阶的固有频率和振型对避免腰骶椎与外界环境发生共振具有重要意义。2.2 稳态动力学分析 对腰骶椎结构进行稳态动力学分析是腰骶椎结构受到谐波激励下的稳态响应，通过稳态动力学分析，可以了解简谐载荷对人体腰骶椎结构的影响。 将激励力为10 N垂直施加在L1上表面，约束骶椎下表面的全部自由度使其固定，以模拟人体上身振动对腰骶椎的影响。并根据模态分析结果选取频率范围为0-5 Hz。根据腰椎解剖结构特点分别在正常腰骶椎模型和内固定后腰骶椎模型的腰椎L1，L3，L5节段棘突处选取代表性节点，节点编号分别为A，B，C，通过计算，获得选取节点X，Y，Z 三个方向的振动幅值随频率变化情况，如图4-6所示。图中横坐标为简谐载荷的频率，单位为Hz，纵坐标表示节点的振动幅值，单位为mm。图4-6分别为两种不同腰骶椎模型的腰椎L1，L3，L5棘突处节点A，B，C位移对频率的响应曲线，比较各个模型的腰椎棘突处3个节点A，B，C在各个方向上的位移最大峰值，可知节点A处位移最大峰值最大，节点B处位移最大峰值次之，节点C处位移最大峰值最小。表明腰骶椎模型的位移峰值从上往下逐渐减小。由图可知内固定后腰骶椎模型的节点A，B，C在X方向上的位移最大峰值比正常腰骶椎模型位移最大峰值都要稍微增大，但在Y，Z方向上位移最大峰值比正常腰骶椎模型位移最大峰值均明显减小。表明螺钉内固定系统对腰骶椎部分起到了保护作用，使得固定后的腰骶椎在受到外界激励时振幅减小；正常腰骶椎模型在1.0-2.0 Hz频率范围内各节点X，Y，Z方向位移均出现峰值，表明该腰骶椎结构对1.0-2.0 Hz频率范围比较敏感，这分别与模态分析得到的第一阶振型的前后俯仰模态固有频率1.33 Hz和第二阶振型的侧向俯仰模态固有频率1.65 Hz相近；内固定后腰骶椎模型在1.0-2.0 Hz频率范围内各节点X，Y，Z方向位移均出现峰值，表明该腰骶椎结构对1.0-2.0 Hz频率范围比较敏感，这分别与模态分析得到的第一阶振型的前后俯仰模态固有频率1.57 Hz和第二阶振型的侧向俯仰模态固有频率1.86 Hz相近；表明这两种腰骶椎结构在这些频率范围内容易产生共振，为避免发生共振，应尽量使人们所处环境的外界振动频率远离腰骶椎固有频率。3 讨论 Discussion脊柱是人体的重要力学机构，对于维持人体的正常生理活动起非常重要的作用。腰椎是承重和活动范围最大的，也是脊柱类疾病的主要发病部位，因为腰椎承载了人体近2/3的质量，遭受外界环境激振后会产生持续时间较长的衰减振动导致其损伤。中国是劳动密集型社会，体力劳动者数量庞大，故脊柱类疾病的防治更应受到重视。本研究基于健康人体脊柱腰骶椎的CT图像，结合医学图形处理软件Mimics 16.0进行三维重建，利用逆向工程软件Geomagic Studio12.0对模型进行光滑和优化处理， SolidWorks三维建模软件系统进行螺钉内固定系统的建模及有限元分析软件Abaqus建立了相对符合实际腰骶椎的三维有限元模型，并对其有效性进行验证，证明模型是有效的。并利用有限元分析软件Abaqus对固定前后腰骶椎模型进行模态分析和稳态动力学分析。其中腰骶椎的模态分析是进一步进行动力学分析的基础。研究结果显示，内固定后腰骶椎模型的节点A，B，C在Y，Z方向上位移最大峰值比正常腰骶椎模型位移最大峰值均明显减小，表明螺钉内固定系统对腰骶椎部分起到了保护作用，使得固定后的腰骶椎在受到外界激励时振幅减小。内固定后腰骶椎模型比正常腰骶椎模型对外界振动敏感性减弱。两种模型在1.0-2.0 Hz频率范围内各节点X，Y，Z方向位移均出现峰值，表明腰骶椎结构对1.0-2 0 Hz频率范围比较敏感，腰骶椎结构在这些频率范围内容易产生共振，为避免发生共振，应尽量使人们所处环境的外界振动频率远离腰骶椎敏感频率范围。本文建立了固定前后脊柱腰骶椎的三维有限元模型，将模型进行有限元分析，得出确定腰椎的固有频率、振型和振幅等振动参数。通过谐响应分析研究简谐载荷对人体脊柱腰骶椎L1-S1节段的影响。使用有限元进行分析，过程直观方便，模型各部分材料属性采用大多数文献经过验证合理的数值，得出的数据结果更为科学。由于建立模型和定义材料时做了一些简化处理，使得实验模型与实际人体模型还有些差距，但是利用简化的模型进行有限元分析得到的结果，代表了一定的趋势。本文建立的模型能够在一定程度上反映腰骶椎的生物力学特性，对于腰骶椎的振动特性分析方面具有重要意义。致谢：感谢河北工业大学张顺心教授及天津医科大学韩立副教授的指导。作者贡献：韩立、张顺心、范顺成、李晔、武晓丹负责课题设计，武晓丹、贾少薇、谢俊德进行资料收集，武晓丹成文，韩立、张顺心评估、审校。利益冲突：所有作者共同认可文章内容不涉及相关利益冲突。伦理问题：临床试验研究的实施符合赫尔辛基宣言和医院对人体研究的相关伦理要求。文章的撰写与编辑修改后文章遵守了随机对照临床试验研究报告指南(CONSORT指南)或观察性临床研究报告指南(STROBE指南)。人体脊柱的CT扫描数据获取经患者本人知情同意，且签署“知情同意书”文章查重：文章出版前已经过CNKI反剽窃文献检测系统进行3次查重。文章外审：文章经国内小同行外审专家双盲外审，符合本刊发稿宗旨。作者声明：文章第一作者对研究和撰写的论文中出现的不端行为承担责任。论文中涉及的原始图片、数据(包括计算机数据库)记录及样本已按照有关规定保存、分享和销毁，可接受核查。文章版权：文章出版前杂志已与全体作者授权人签署了版权相关协议。开放获取声明：这是一篇开放获取文章，文章出版前杂志已与全体作者授权人签署了版权相关协议。根据知识共享许可协议“署名-非商业性使用-相同方式共享3.0”条款，在合理引用的情况下，允许他人以非商业性目的基于原文内容编辑、调整和扩展，同时允许任何用户阅读、下载、拷贝、传递、打印、检索、超级链接该文献，并为之建立索引，用作软件的输入数据或其它任何合法用途。4 参考文献 References1 Li XF, Liu ZD, Dai LY, et al. Dynamic response of the idiopathic scoliotic spine to axial cyclic loads. Spine. 2011; 36(7): 521-528.2 Guo LX, Teo EC, Lee KK, et al. Vibration characteristics of the human spine under axial cyclic loads: effect of frequency and damping. Spine. 2005;30(6): 631-637.2393ISSN 2095-4344 CN 21-1581/R CODEN: ZLKHAH3 Kong WZ,Goel VK. Ability of the finite element models to predict response of the human spine to sinusoidal vertical vibration. Spine. 2003;28(28): 1961-1967.4 Seidel H, Blthner R, Hinz B. Application of finite-element models to predict forces acting on the lumbar spine during whole-body vibration. Clin Biomech. 2001;16Suppl 1(16 Suppl 1): S57-S63.5 Kasra M, Shirazi-Adl A, Drouin G. Dynamics of human lumbar intervertebral joints. Experimental and finite-element investigations. Spine. 1992;17(1): 93-102.6 黄凯,刘展亮,刘少喻,等. 腰椎单侧与双侧椎弓根螺钉固定的生物力学对比性研究J. 中华临床医师杂志:电子版, 2015(22): 4143-4147.7 郝剑,姚进,朴哲,等. 腰椎坚强固定后邻近节段的有限元法生物力学分析J. 中国中西医结合外科杂志,2016,22(5):475-478.8 Zhu R, Yu Y, Zeng ZL, et al. A Review of the Static Loads Applying on the Finite Element Models of the Lumbar Spine. J Med Imaging Health Informat. 2015.9 Zafarparandeh I. Application of Finite Element Method in the human Spine Biomechanics. 2016.10 Zhang Z, Yang L I, Liao Z, et al. Research Progress and Prospect of Applications of Finite Element Method in Lumbar Spine Biomechanics. J Biomed Eng. 2016.11 王端阳,奚春阳,杨辉,等. 有限元分析法在腰椎生物力学研究中的应用及新进展J. 现代生物医学进展, 2015,15(9): 1794- 1797.12 Zulkiflil A, Ariffinl AK, Rahman MM. Probabilistic finite element analysis of vertebrae of the lumbar spine under hyperextension loading. Int J Autom Mech Eng. 2011;3(1): 256-264.13 Kavitha A, Sudhir G, Ranjani T S, et al. Implant Analysis on the Lumbar-Sacrum Vertebrae Using Finite Element Method. Springer Singapore, 2017.14 Wang DY, Chun-Yang XI, Yang H, et al. New Progress and Application of Finite Element Analysis Technology in Lumbar Spine Biomechanics. Prog Mod Biomed. 2015.15 Alonso-Rasgado MT, Martinez-Lozada FM, Bailey CG, et al. The effects of L4-L5 Fusion on the Adjacent Segments in Lumbar Spine: Finite Element Analysis. Int J Stroke. 2015; 10(2):269-274. 16 Li K, Zhang J, Jiang J, et al. Lumbar spinal finite element analysis in a gravity environment. in Eighth International Conference on Digital Image Processing. 2016.17 韩硕. 有限元分析在腰椎间盘生物力学中的应用研究D. 河北医科大学, 2016.18 曹丽君.正常及退变下腰椎在不同工况下三维有限元力学分析D.河北医科大学, 2015.19 张伟,宫赫.振动载荷作用下人体腰椎多孔弹性有限元分析C/北京力学会学术年会. 2015.20 张振军,李阳,廖振华,等. 有限元法在腰椎生物力学应用中的研究进展和展望J. 生物医学工程学杂志, 2016,(6): 1196-1202.21 秦计生, 王昱, 彭雄奇, 等. 全腰椎三维有限元模型的建立及其有效性验证J. 医用生物力学, 2013, 28(3):321-325.22 黄菊英,李海云,吴浩. 腰椎间盘突出症力学特征的仿真计算方法J. 医用生物力学, 2012, 27(1): 4630-4631.23 苏晋