基于离散单元法的人工骨材料微观结构研究精选.docVIP

基于离散单元法的人工骨材料微观结构研究 韩琴 1，汪焰恩 1 2，魏生民 1，魏庆华 1 （1.西北工业大学，西安，710072 （2.斯特拉斯克莱德 摘 要：人工骨支架的孔隙率是影响体液在人工骨内部循环和细胞在内部生长能力的重要参数，利用离散单元法分析人工骨支架制备过程中羟基磷灰石微球的堆积过程，根据每一时刻各颗粒间的相互作用计算接触力，再运用牛顿第二定律计算单元的运动参数，从而实现对运动情况的预测。在颗粒的碰撞过程中，运动的颗粒必然趋于稳定，最终计算出稳定状态下人工骨支架的孔隙率和配位数。利用掺杂可溶解生物微球的方法调节人工骨支架的孔隙率和配位数，对人工骨支架孔隙率进行控制，实现了按照个体差异制备人工骨支架。对比前人晶格理论计算人工骨孔隙率的方法，验证了晶格理论的可行性，同时也体现了离散单元法的准确性。 关键词：人工骨支架；孔隙率；配位数；离散单元 随着社会的进步与人们生活水平的提高，人们的平均寿命不断延长，但人体本身的老化却不可避免，世界人口老龄化问题日益突显。全球人口平均寿命已延长到66岁，全球老龄人口总数已达6.29亿。人体的老化30岁就开始了，人体的结缔组织、特别是保持人体直立姿势的骨骼组织随之开始劣化。随着年龄的增长，由于骨组织的变化，其强度逐渐降低，发生骨折的几率增大。此外，由于创伤、感染、肿瘤、事故、战争或矫形手术所致的骨缺损也需要进行骨移植。 骨移植材料大体可分为：自体骨、同种异体骨和人工替代骨，自体骨组织移植效果很好,无免疫排异反应，但是存在来源有限、需要二次手术及取骨部位常遗留慢性疼痛等并发症的缺点[1]体移植虽不需两次手术, 且具有自体骨的一些优越的组织特点, 能提供足量不同形状尺寸的皮质骨或松质骨，但移植骨在经消毒处理后会失去强度, 部分或完全损失骨诱导因子, 并有感染病毒的可能。为了克服这些局限, 1995年Crane 由于不同部位的人工骨有其特定的孔隙率，故基于精确控制孔隙率制备人工骨支架是人工骨支架的一种重要制备技术，有必要对人工骨材料微观结构进行研究和探讨。汪焰恩等人借鉴材料设计晶格理论，对晶格堆积排列的14种布拉伐格子堆积方式[5-6]，得出对于微球堆积人工骨支架而言，简单立方可作为支架单胞分析模型，并以简单立方为模型计算了羟基磷灰石微球堆积的人工骨支架的孔隙率，即配位数为6时，羟基磷灰石微球为0.3时羟基磷灰石微球人工骨支架的孔隙率为；文中拟对系统的理论和数值计算进行分析，从而得出基于晶格理论的孔隙率计算方法和基于离散单元法的孔隙率计算方法之间的内在差异与一致性。 1 基于离散单元法的孔隙率计算方法 1.1 离散单元法的单元接触模型 离散单元法[7]流方法是把分析对象看成充分多的离散单元，每个颗粒或块体作为一个单元，根据整个过程中每一时刻各颗粒间的相互作用计算接触力，再运用牛顿第二定律计算单元的运动参数，这样循环运算，实现对运动情况的预测，其计算流程框图如图1所示： 图1 颗粒流方法计算流程框图 Fig.1 Calculation flow diagram of Particle flow method 颗粒流方法的颗粒模型是将颗粒与颗粒、颗粒与边界的接触采用振动运动方程进行模拟，图2为法向和切向接触模型 （a）法向接触模型 （b）切向接触模型 图2 颗粒流方法接触模型 Fig.2 Contact model of Particle flow method 在颗粒的碰撞过程中，运动的颗粒最终必然要趋于稳定，也即动能转化为热能，因而系统中单元的运动是一不可逆的过程 ，为了耗散动能，离散元计算中引入阻尼来吸收系统的动能。阻尼系数的考虑有质量阻尼和刚度阻尼，两者的作用分别为限制颗粒的绝对运动和相对运动。 1.2 单元的运动模型 颗粒流方法只需要满足两组基本方程, 即运动方程( 牛顿第二运动定律) 和物理方程( 力-位移关系) 。采用中心差分公式进行动态松弛求解，其具体的计算过程如下： 在任意时刻t，任一单元受力作用的运动可由牛顿第二定律求得。由于单元体在运动过程中消耗功能, 需引入与运动反向与质量成比例的阻尼力, 则单元体的运动方程为： (1) (2) 式中 m、,、分别表示单元体的质量、t时刻的位移和合力;、、分别表示单元体的转动惯量、t 时刻的角位移和合力矩;a为与质量成比例的粘性阻尼系数。将式( 1) 、式( 2) 中的加速度和角加速度写成中心差分形式：