三维有限元头颅模型参数及边界条件研究培训讲学.doc

此文档收集于网络，如有侵权请联系网站删除三维有限元头颅模型参数及边界条件研究一l94一基础研究?生物医学工程与II缶床2005年7月第9卷第4期BME&ClinMed,July2005.Vo1.9,No.4三维有限元头颅模型参数及边界条件研究何黎民,卢亦成,吴建国,刘平,丁祖泉,陈学强,王保华摘要:目的确定三维有限元模型材料阻尼参数,实现寰枕关节模拟.方法对不同材料阻尼参数下的模型进行冲击载荷加载,分析比较计算后不同材料阻尼下颅骨应力一时间曲线形态,确定模型采用的材料阻尼参数值.应用弹簧单元模拟寰枕关节,比较模拟前后,模型在冲击载荷下颅骨应力,确定弹簧单元参数.结果不同颅骨材料阻尼条件下,额部冲击区域颅骨节点的yonMises应力曲线均无第2峰值,应力峰值随阻尼增大而逐渐降低,峰值出现时间后移,枕部颅骨节点应力曲线在冲击后期应力下降趋势更显着,在0.0010.004之间颅骨应力曲线形态最佳.模拟寰枕关节,弹簧单元采用颅骨材料参数,在X,l,轴方向弹簧弹性系数10N/mm,Z轴方向的为20N/mm,颅骨应力曲线后期可下降.结论颅骨线弹性材料阻尼系数对模型颅骨应力响应有显着影响.模型脑组织应力主要受颅骨应力的影响.采用三维弹簧单元模拟寰枕关节可降低模型颅骨应力,模拟寰枕关节对模型颅骨应力的影响较颅骨材料阻尼系数的影响小.关键词:有限元模型;头颅;材料参数;边界条件中图分类号:R68文献标识码:A文章编号:10097090(2005)04019404StudyonparametricandboundaryconditionofthefmiteelementYi-chen,wuJianguo,LIUPing,DINGZu-quan2,CHENXueqiang2,WANGBanhuaNeurosurgery,ChangzhengHospital,Shanghai200003,China;2.BioscienceandMedicalEngineeringUnivers,Shanghai200072,China)Correspondingauthor:.LUYicheng,E-maihlychengCHELimin,LU(J.DepartmentofAcademy,TonalAbstract:ObjectiveTostudymaterialdampingparametersandboundaryconditioninfiniteelementmodelofthehumanhead.MethodsInthethreedimensionalhumanfiniteelementmodel,brain,skuHanddurawereassumedtobehomogeneousandisotropicmaterialwithlinearelasticbehavior.Aseriesofstructuraldampingparametersofskullandbrainwereadoptedtosimulateafrontalheadimpact.ThecalculatedresultsofyonMisesstressofskullandbrainwereobserved.Using2Dspringelements,theatlantooccipitaljointwassimulatedinthreedimensionalspace.ResultsTherewerenosecondpeakintheyonMisesstresscurveofskull,whenskulldampingparameterWasdefined.ThehigherdampingparameterWasadopted,thelowerpeakvalueofskullyonMisesstresscurvewasobserved.Itispropertodefineskullstructuraldampingfactorbetween0.001and0.004.AftertheatlantooccipitaljointWassimulated,withdifferentelasticcoefficients(10ofspringelementindifferentdirections.10N/mminX-Ydirectionand20N/mminZdirection,theyonMisesstresscurveofskullbecamemoresmoothly.ConclusionsTheskulldampingparameterinfluencetheyonMisesstressofskullgreatly.Thechangeofbrmnstresscurvesassociatewithskullstresscurvesclosely.Keywords:finiteelementmodel;humanhead;materialproperties;boundaryconditions由于头颅三维有限元模型涉及的参数较多,在模拟计算直接冲击造成的头颅结构应力变化前,需确定其边界条件和模型参数.在模型参数研究方面,根据既往颅脑生物力学研究结果,制作模型各种结构的材料均为各向同性的线弹性材料;为提高计算效率,作者单位:1.上海第二军医大学长征医院脑外科,上海2o0loo3;2_上海同济大学生命科学与医学工程研究院,上海2O0072收稿日期:20041014;修回日期:20o50l一05作者简介:何黎民(男,1972一,河南荥阳人),博士,主治医师,主要从事神经外科工作和颅脑冲击伤有限元模型研究.基金项目:上海市科技发展基金项目(0141l9o99)通信联系人:卢亦成,Email:lycheng减小计算成本,首先构建头颅模型,重点研究线弹性材料阻尼系数对应力响应的影响,以及如何实现对寰枕关节的模拟.1材料与方法1.1仪器与软件MX8000多排螺旋CT(Philips意大利),Vision_plus1.5T核磁共振(Simenz德国).Dimension650图形工作站(Dell美国).ANSYSY有限元软件.1.2构建头颅有限元模型依据正常人男性国人横轴位CT,MRI图像,半自竺堕兰堡皇I堕7月第9卷第4期BME&ClinMed,July2005,Vo1.9,No.4动识别重建结构,对皮肤,颅骨等实体结构采用灰度阈值法进行图像分割,用边界提取算法进行边界提取,再用轮廓跟踪算法和B样条曲线拟合方法,提取边界关键点并拟合边界轮廓.对其他形态复杂的结构进行子结构划分,人工绘制边界,确定关键点,用自编坐标提取程序(由上海同济大学生命科学与医学工程研究院编写)提取关键点坐标,形成有限元代码命令流文件,导人有限元软件,对边界进行空间线性拟合.采用实体单元,壳单元,划分网络.筛选不同结构相邻界面上最小距离接点组合,生成粘连节点命令流文件,在前处理中完成节点粘连.结构的头颅三维有限元模型结构包括:皮肤,颅骨,硬脑膜(含大脑镰,小脑的透露三维有限元模型),大脑半球,中脑,小脑和脑干等.模型总线数804,总面数412,总体积数24,总单元数98744,总节点数42781.1.3模型材料基本参数依据既往对头颅材料的研究结果【I】和既往模型材料参数【2-51,采用制作皮肤,颅骨,硬脑膜脑组织材料为各向同性线弹性材料,其参数有杨氏摸量(MPa),泊松比(和密度(p/kg/m3)取值见表1.表1模型线弹性材料参数1.4用二维有限单元模拟寰枕关节在模型枕骨大孔层面的枕骨上确定代表枕骨髁关节面的节点,图1显示选择的节点,左右两侧各19个,记录节点序号.选用弹簧单元的参数包括:弹性常数,阻尼系数c,弹簧单元材料与颅骨材料相同,但无质量.输出选定枕骨髁关节面节点的坐标,在y,Z,yZ三个坐标轴平面生成与这些节点对应的3组新节点,共ll4个,在对应的每两个节点间定义二维单元,单元类型定义为上述弹簧单元.生成的二维单元与颅骨枕骨大孔的关系如图2所示.即使用3组二维单元模拟寰枕关节.1.5加载冲击载荷采用Nahum61尸体实验的载荷条件.选择额部18个颅骨节点,在y,Z轴方向分别加载冲击力载荷,使合力方向等同于Nahum实验受力方向,作用力方向平行Frankfort解剖平面,指向枕部,所选18个节点围成面积约754mm.冲击载荷曲线如图3所示.冲击总时程15ms,分l1个载荷步.书写有限元分析软件载荷命令流文件.采用瞬态小变形分析模式,控制注:粗点所示范围为枕骨髁关节面图l枕骨大孔层面枕骨髁关节面节点.195-注:线条表示二维弹簧单元图2模拟寰枕关节的单元与枕骨大孔关系最大子载荷步数为5,选择稀疏矩阵算法进行求解.Z时间/ms图3Nahum实验冲击载荷时间历程2结果O2.1颅骨材料阻尼系数对应力响应的影响图4显示不同颅骨材料阻尼条件下,模拟计算获得的额部冲击区域颅骨节点的vonMises应力一时问曲线.阻尼为O时,模拟此处颅骨应力峰值与冲击峰值出现时间相似,应力峰值高,冲击后期出现第2峰.设定阻尼参数后,颅骨应力峰值随阻尼系数增大而逐渐降低,峰值出现时问后移,冲击后期无第2峰.阻尼系数为0.004时,颅骨应力在曲线末期的值最小.在相同载荷下,模拟计算对冲部位(枕部)颅骨节点yon一l96-塑垦_T程与Il缶床2005年7月第9卷第4期BblE&ClinMed,July2005,Vo1.9,No.40.1250.3750.6250.8751.1251.3751.625时间/s0.1250.3750.6250.8751.1251.3751.625时间/sc=0.0o80.1250.3750.6250.8751.1251.3751.625时间/sx10-2x10-2x10-287.26.45.64.843.22.41.60.800.1250.3750.6250.8751.1251.3751.625时间/s0.1250.3750.6250.8751.1251.3751.625时间/sc=0.010.1250.3750.6250.8751.1251.3751.625时间/s图4不同颅骨材料阻尼时额部颅骨节点yonMises应力一时间曲线Mises应力,显示材料阻尼对颅骨应力响应有类似影响.模拟寰枕关节几弹簧单元参数对颅骨应力的影响图5显示采用不同方法模拟寰枕关节时,颅骨对冲击载荷的应力一时间曲线(AD).A完全约束枕骨髁平面;B仅用Z轴方向弹簧模拟寰枕关节,弹簧单元K=10N/mm;C采用,y,Z三个轴向的弹簧单元模拟寰枕关节,y轴方向弹簧单元K=10N/mm,Z轴方向弹簧单元K=20N/mm;D在C参数基础上,增加颅骨阻尼参数c=0.00l.A的颅骨应力一时间曲线峰值最高,末期下降不佳,存在第2峰值;B的曲线峰值基本不变,曲线末期第2峰值较A的有改善;C与B结果基本相同;D的颅骨应力一时间曲线有显着改善.模拟计算结果还显示,单独设定弹簧单元的阻尼系数或弹簧单元材料阻尼应力一时间曲线影响情况相似,当K>500N/mm时,颅骨应力一时间曲x10-2x10-2X10-202505OO750l0O0l2501500时间/s注:A完全约束,固定枕骨髁关节面,B仅用Z轴弹簧单元模拟寰枕关节,K=20N/mm,C用,y,Z三轴向弹簧单元模拟寰枕关节为20N/mm,10N/mm,D在C条件下增加定义颅骨材料阻尼c=0.001图5同一颅骨节点在不同边界条件下的应力曲线963086420864200一L100一L10勰加840Bd0差L1086420864000一L1009876543210Bd目.羞LI.Bd臼.差LI.Bd0羞L10生物医学工程与l晦床2005年7月第9卷第4期BME&ClinMed,July2005,Vo1.9,No.4线与边界完全约束的应力一时间曲线趋于一致.3讨论早期材料实验证实颅骨等生物组织具有均质线弹性材料特征【lj,这种材料特性在模拟计算中较易实现,但材料阻尼参数对冲击载荷下应力响应的计算结果存在显着影响,主要表现为冲击末期结构应力存在震荡表现171,应力一时间曲线出现第2,第3峰值.此计算结果与头颅尸体冲击实验测得的颅内压曲线不符合,Shuck的材料力学研究显示脑组织具有黏弹性特征,因此为实现应用线弹性材料特征进行模拟,需要设定适当的材料阻尼参数.Shuck等人测定的脑组织阻尼系数为0.00670.0087,Churn与Huangl8模型采用的颅骨,脑阻尼系数为0.001,本实验应用线弹性材料特征定义头颅结构时也发现上述应力响应现象,因此分别对颅骨,脑组织线弹性材料阻尼系数进行计算研究.结果显示颅骨阻尼系数的设定可降低颅骨yonMises应力峰值,阻尼系数为0.004时,冲击部位下颅骨应力曲线与冲击载荷曲线的形态最相似.不同阻尼条件下,对冲部位颅骨应力曲线也有上述表现,但曲线较冲击部位在时间上延迟,考虑这种情况与冲击在颅骨上传递耗时相关.头颅在受到冲击时,通过颈限制头颅运动,所以要考虑这种边界条件.因为实验冲击载荷总作用时间为10ms,在这一短暂时间内虽可完成单突触反射,但因为头颅此时整体位移很小,颈肩部骨骼肌及其附属结构的感受器尚未完全被激活,故颈肩肌肉对头颅的约束尚不充分,所以在冲击时程中颈部对头颅的约束作用主要通过寰枕关节韧带完成.因此,采用,y,Z轴三个方向的3组二维弹簧单元与枕骨髁关节面相连,模拟整个颈部,主要是颅颈交界的韧带和颈椎的韧带对头颅的约束.正常情况下,寰枕关节在水平面上可发生的一定位移,加之寰枕关节及各颈椎关节突关节在水平面上的位移,使头颅易于在水平面上转动.水平面即模型坐标的,y轴平面.在纵向上,即模型坐标的Z轴方向,寰枕关节韧带和颈椎间的纵向韧带不发生显着位移,较强地限制头颅在纵向上的运动.所以我们对水平方向和纵向弹簧采用了不同的值,y轴方向弹簧单元的较小,Z轴方向的较大,这将更好地模拟应力在颅内的分布.计算结果显示完全约束(枕骨髁完全固定)的边界条件下颅骨应力一时间曲线峰值最高,存在第2峰值,末期曲线下降不佳.应用弹簧单元模拟寰枕关节后,颅骨应力一时间曲线峰值下降,第2峰及曲线末期改善.单用Z轴方向弹簧单元模拟寰枕关节与应用三个轴.向弹簧单元模拟关节时颅骨应力一时间曲.197-线有差异,但无统计学意义,这可能与实验冲击载荷作用于额部中线上,且冲击力方向平行于头颅矢状面有关,当载荷作用部位与冲击力方向变化后,应力三个轴向弹簧单元模拟寰枕关节将更有效.从寰枕关节解剖角度分析,三维弹簧单元模拟更接近解剖实际.单独设定弹簧单元阻尼或弹簧单元材料阻尼对颅骨应力无明显影响,表明此两个参数基本不影响颅骨的应力响应,相比较而言,颅骨材料的阻尼系数在相同条件下对颅骨应力一时间曲线影响最显着,应力曲线峰值显着下降,曲线无第2峰值,末期曲线下降良好,与冲击载荷曲线形态基本一致.弹簧单元值在一定范围内(10一l00N/mm)对颅骨应力曲线的影响较相似,当大于500N/mm后,颅骨应力一时间曲线与边界完全约束时曲线趋于一致,表明模拟弹簧硬度不宜过大,否则起不到模拟关节作用.既往部分有限元模型在冲击中假设头颈边界为自由边界【8?91,这对于时程较短的冲击而言,与不模拟寰枕关节所计算的应力相似,但我们认为,为更接近实际地模拟冲击,扩大模型模拟适用范围,应当模拟寰枕关节.参考文献:【1】RuanJS,KhalilT,KingAI.HumanheaddynamicresponsetosideimpactbyfiniteelementmodelingU】.JBiomechEng,1991,113(3):276-283.【2】KumaresanS,RadhakrishnanS.Importanceofpartitioningmem?branesofthebrainandtheinfluenceoftheneckinheadin.jurymodelingJ.MedBiolEngComput,1996,34(1):27-32.【3】MesSS,ThibaultKLInfantskullandsutttreproperties:mea?surementsandimplicationsformechanismsofpediatricbraininjuryJ.JBiomechEng,2000,122(4):364-371.【4】ZhangLY,HardyW,OmoriK,eta1.Rencentadvancesinbraininjuryresearch:anewmodelandnewexperimentaldataJ】.BioengineeringConference,2001,50:833834.【5】KleivenS.InfluenceofimpactdirectiononthehumanheadinpredictionofsubduralhematomaJ.JNeurotrauma,2003,20(4):365379.【6】NahumAM,SmithR,WardCC.Intracranialpressuredynam?icsduringheadimpact【C】.ProceedingsofStappCarCrashConference,1977.339366.【7】ChuCS,LinMS,HuangHM,eta1.FiniteelementanalysisofcerebralcontusionJ.JBiomech,1994,27(2):187-194.【8】HuangHM,LeeMC,ChiuWT,eta1.ThreedimensionalfiniteelementanalysisofsubduralhematomaJ.JTrauma,1999,47(3):538-544.【9】WillingerR,KangHS,DiawB.Three-dimensionalhumanheadfiniteelementmodelvalidationagainsttwoexperimentalim?pactsJ.AnnBiomedEng,1999,27(3):403-410.此文档仅供学习和交流