关于有限元法_正则化的弹性模量反求算法研究.pdfVIP

第 26 卷 　第 1 期 应 用 力 学 学 报 Vol . 26 　No . 1 2009 年 3 月 CHINESE JOURNAL OF APPL IED MECHANICS Mar . 2009 文章编号 (2009) 0 1006006 基于有限元法正则化的弹性模量反求算法研究 蔡传宝 　汤文成 ( 东南大学 　2 11189 　南京) 摘要 : 生物组织的病变往往导致其物理特性的变化 ,研究生物组织的物理特性变化对诊断生物组 织病变与否能够提供参考 。本文提出了根据已知生物组织部分边界位移信息、外部载荷 、边界条件 和有限元模型 ,利用正则化方法反求弹性模量的方法 。首先将求解弹性模量问题描述成一个力学 反问题 ,建立逆系统辨识的数学模型 ,并应用正则化原理建立生物组织弹性模量反求算法 。最后利 用仿真算例验证了算法的有效性 ,并结合图像识别技术试验获取对象部分边界位移信息 ,反求其弹 性模量 ,证明了本文反求方法的可行性 。 关键词 :弹性模量 ;反问题 ;正则化 ;有限元方法 中图分类号 : TP39 1 　　　文献标识码 : 　A 下 ,通过记录模型变形前后的影像 ,运用图像识别及 1 　引　言 摄影测量学中影像内定向的有关技术 ,实现了对模 型变形前后位移的非接触式测量 ,利用观测点的位 生物组织的病变往往导致其物理特性的变化 , 移数据 ,在已知外部载荷 、边界条件和有限元模型的 因此生物组织的物理特性参数可以有效地表征其病 基础上 ,使用计算获取生物组织弹性模量 ,最后给出 变情况 ,其量化已成为研究的热点 。国内外学者利 算例实验和仿体实验以检验方法的正确性 。 用在一定外载荷下对生物材料直接进行应变试验获 取弹性模量 ,或者基于生物材料的本构关系 ,利用光 2 　反问题的不适定性 学 、超声波学间接获取所有的试验对象的位移信息 求得弹性模量 。 当应用经典的方法去求问题的解时 ,给定的输 人体软组织弹性模量最早采用传统的力学实验 入数据不一定能保证解的存在性或唯一性 ;或者输 方法如离体或活体应力应变试验方法测取[ 14 ] 。文 入数据的微小波动会引起相应解的巨大变化 ,而且 献[56 ]等分析血管弹性成像领域的现状 ,提出所谓 这种变化已经使得用通常方法求得的对应解失去意 弹性成像实质上只是应变成像 ,在很大程度上受弹 义 ,这种现象的产生的原因是原问题的不适定性 ,不 ( ) 性模量分布 、组织形状 、位置 、外力类型等边界条件 适定性的本质上是由于信息 输入数据 、待求的解 ( ) 的限制 ,并不显示真正的血管弹性 。文献[7 ]利用脉 不足 或过定 造成的。恢复问题的适定性尤其是稳 搏波结合智能算法测算血管壁的弹性系数 。 定性的方法有添加信息、改变拓扑度量等 ,因为度量 本文结合力学反问题的研究 ,考虑在无法利用 方式在给定的应用问题中通常难以随便改变 ,所以