本研究提出一种耦合实验与有限元模拟的方法,用于获取软生物材料(脑、肝、肌腱、脂肪)的单轴动态力学响应 等等。通过生物材料有限元分析的迭代优化,将分离式霍普金森压杆试验中因试样鼓胀而获得的多轴实验结果转化为单轴真应力-应变行为。
需要JoVE订阅才能观看此内容。 请登录或开始免费试用
本研究提出一种耦合实验与有限元模拟的方法,用于获取软生物材料(脑、肝、肌腱、脂肪)的单轴动态力学响应 等等。通过生物材料有限元分析的迭代优化,将分离式霍普金森压杆试验中因试样鼓胀而获得的多轴实验结果转化为单轴真应力-应变行为。
本研究提供了一种结合实验与有限元(FE)模拟的方法,用于研究软生物材料的力学行为例如 脑、肝、肌腱、脂肪 等。) 在受到高应变率作用时。本研究采用分离式霍普金森压杆(Split-Hopkinson Pressure Bar, SHPB)产生 100–1,500 sec⁻¹ 的应变率。-1. 采用由粘弹性材料(聚碳酸酯)制成的撞击杆的分离式霍普金森压杆(SHPB)装置。在死亡后不久获取生物材料样本,并制备用于SHPB测试。将试样置于入射杆和透射杆之间,启动SHPB的气动组件,驱动撞击杆撞击入射杆。由此产生的冲击生成一个压缩应力波(即 入射波)穿过入射杆。当压缩应力波到达入射杆末端时,一部分继续向前穿过样品和透射杆(即 透射波),而另一部分则在入射杆中以拉伸波的形式反向传播(即 反射波)。这些波通过安装在入射杆和透射杆上的应变片进行测量。试样的真实应力-应变行为是基于波传播理论和动态力平衡方程计算得到的。由于试样在加载过程中发生鼓胀,实验测得的应力-应变响应本质上具有三维特性,因此采用静水应力(第一不变量)来构建应力-应变关系。为了提取组织的单轴(一维)力学响应,结合实验结果与包含组织内部状态变量(ISV)材料模型的有限元分析(FEA),进行了迭代耦合优化。该优化过程中,用于模拟实验装置的有限元仿真中的ISV材料模型被反复校准(即 优化后的实验数据与有限元分析(FEA)的应变片数值及应力的第一不变量吻合良好。
动机
软生物材料(如脑、肝、肌腱、脂肪)的耦合分离式霍普金森压杆(SHPB)实验与有限元建模研究的核心目标 等等。目的是提取其单轴力学行为,以便进一步应用于人体在损伤性机械载荷下的有限元(FE)模拟。该人体有限元(FE)模型包含详细的人体网格以及适用于多种人体器官的历史相关多尺度黏弹性-黏塑性内变量(ISV)材料模型。该人体模型可作为构建更优损伤防护标准、设计创新防护装备以及实现以乘员为中心的车辆设计的框架。
在人类创伤中,已广泛观察到两种高速损伤模式:爆炸冲击和钝性撞击。爆炸性武器造成的损伤是创伤性损伤(TI)的主要来源,也是战场上的首要致死原因1。当这些炸药引爆时,会产生向外传播的冲击波,引发剧烈且突然的加速度和形变,对暴露个体构成严重威胁。尽管人体任何解剖部位都可能受到冲击波伤害,但最受关注的部位主要有两个:(1)下肢,因其靠近地面;(2)头部,因为头部损伤可能影响正常脑功能并危及生命2,3。根据所受损伤类型的不同,这些损伤可分为原发性、继发性或三级损伤。由于爆炸强度由其重量或尺寸、 standoff 距离、正压脉冲持续时间以及传播介质等因素决定,因此准确分类这些损伤具有挑战性3-6。国会报告指出,从2000年至2010年3月,在伊拉克和阿富汗战争中,军事人员因爆炸性武器和车辆碰撞已遭受近17.9万例创伤性损伤2。由于现代战争的性质和作战地点特点,头部损伤已成为军方和民众共同关注的重点问题3。
除战斗场景外,创伤性损伤(TI)还由多种原因引起,包括机动车创伤、竞技运动、摩托车和家庭事故以及运动损伤。例如,尽管安全设备和防护措施已有所改进,但机械性创伤性脑损伤(TBI)在美国仍是导致死亡和终身病残的主要原因之一。美国疾病控制与预防中心(CDC)报告称,每年约发生140万例TBI事件,其中近5万例导致死亡。仅美式足球一项每年就造成超过30万例TBI病例7。此类损伤的幸存者面临与感觉、认知和沟通相关的长期神经系统并发症风险。目前,约有530万美国人正承受着这些慢性功能障碍和残疾。2000年至2010年间,美国与TBI相关的直接和间接医疗费用总计达600亿美元8。然而,这些数据并未计入非医疗支出和损失,也未包括TBI患者家庭成员和亲友所承担的负担。除单纯的经济分析外,TBI所致残疾还会显著降低患者生活质量,并对家庭和社会造成沉重负担。
进一步理解创伤性脑损伤(TI)的形成机制、特征描述及预防措施的需求显而易见。针对导致TI的潜在机制开展生物力学研究,有助于深入认识并提供减少高危人群暴露风险或改进安全防护措施的机会。此外,对TI形成机制的更深入认识,可能改进诊断方法与诊断标准,为治疗TI的医疗专业人员提供更有效的手段,以改善临床预后并挽救生命。
要制定有效的人体防护措施,需要更深入地了解损伤机制以及损伤发展过程中的生物力学特性。从历史角度来看,旨在预测损伤的模拟研究一直受到计算能力限制以及所采用的解剖和材料模型精度不足的制约。全身模拟通常关注各身体部位所承受的整体载荷,但对各个器官、肌肉、骨骼中的局部应力、应变和损伤情况的研究仍显不足。 等等。 尚未观察到此类现象。例如,肩部力矩模型利用手臂的尺寸、载荷以及施加的角度来查找表格中的数值,以判断特定场景是否存在危险。这类计算对于快速估算很有帮助,但无法捕捉从手部到肩部局部区域的实际状况,尤其是当损伤和伤害本质上具有局部性时。其次,有限元(FE)模拟已被用于获取局部响应。这些研究的局限性并不在于有限元分析(FEA)本身,而在于定义人体各部位在爆炸载荷下行为的材料模型。以往采用的材料模型多借鉴自较简单的材料,未能充分描述生物组织所表现出的多种复杂力学行为。因此,采用内状态变量(ISV)材料模型的高保真计算模型,对人體器官进行建模,是研究创伤性损伤(TIs)物理机制与生物力学特性、设计创新防护装备以及建立更优损伤评估标准的最真实可行的方法。
分离式霍普金森压杆(SHPB)与内变量状态(ISV)材料模型背景介绍
由于涉及伦理问题 在 体内 人体器官测试及大规模人体遗骸测试相关的物流问题,当前研究工作涉及机械实验 在 体外 使用从动物替代品中提取的器官制备的标本(例如,猪作为最常用的替代模型)。聚合物SHPB一直是首选方法 体外 在高应变率下测试软生物材料。将SHPB测试获得的相关变形行为,以及从组织微结构特征中提取的与组织损伤相关的信息,整合到我们的ISV材料模型中,用于器官力学特性的描述9-10。然后将这些材料模型导入虚拟人体模型中,以开展针对多种损伤的有限元分析。该过程使我们能够朝着准确预测特定器官在不同力学载荷条件下损伤的物理机制与本质这一目标迈进。例如 爆炸冲击、车祸及钝性撞击)提供精确的生物力学响应,而无需进一步的物理实验。为了准确描述用于人体有限元(FE)仿真中生物材料的现象学力学特性,特别是其在高应变率下的依赖性,研究人员对这些生物材料开展了分离式霍普金森压杆(SHPB)实验,以获取与人体创伤性损伤(TIs)相关应变率条件下的动态力学响应。本文展示了位于密西西比州立大学(MSU)先进车辆系统中心(CAVS)的SHPB实验装置的总体概况。 图1.
先前的研究表明,SHPB测试存在三个主要缺陷12-18第一个也是最主要的问题是材料的惯性效应,该效应在生物材料试样高应变率力学响应中表现为初始的尖峰。为解决这一问题,先前的研究建议将试样几何形状从圆柱形改为立方体形或环形。由于试样几何形状影响了应力波传播、波的相互作用以及力学响应,这些研究所得的力学行为彼此之间存在差异。此类对试样几何形状的修改导致了对生物材料力学响应(多轴和非均匀应力状态)的错误表征。第二个主要缺陷是在测试过程中无法维持动态力平衡。研究人员通过减小试样厚度与直径的比值和/或在测试前冷冻组织来克服这一问题。虽然减小厚度与直径的比值解决了动态力平衡问题,但冷冻组织因组织中水分结晶而改变了材料性能,进一步复杂化了测试流程。一些研究完全放弃了使用SHPB,以避免上述缺陷,转而采用激波管在多种动物模型(大鼠、猪)中获取压力-时间响应 等等。然而,这些动物模型无法提供有限元模拟中所用材料模型所需的单轴一维应力-应变行为。第三个缺陷是,由于材料较软以及试样含水量较高,导致试样在SHPB实验中发生鼓胀(barreling),从而无法获得一维的应力-应变结果。
因此,分离式霍普金森压杆(SHPB)是一种获取高应变率数据的有效测试装置。然而,对于软材料而言,SHPB会引发鼓胀现象,从而主要由静水压力产生三维应力状态,而实际所需的是单维应力-应变数据。本文展示了如何仍可利用SHPB获取用于材料模型标定的单轴真实应力-应变曲线;但获得该单轴真实应力-应变曲线的过程较为复杂。该过程既包含多轴实验数据,也涉及有限元(FE)模拟结果,并需要对材料模型常数进行迭代重标定。在MATLAB中实现的一维ISV材料模型(亦称材料点模拟器)需要一维实验数据进行标定。因此,ISV材料模型通过系统化的标定流程进行了优化。在此过程中,SHPB实验数据结合波理论公式化方法与动态力平衡分析(MSU高应变率软件)进行处理。为考虑聚合物SHPB的粘弹性色散效应,Zhao et al.(2007)报道的粘弹性色散方程被集成至MSU高应变率软件中,以确保测试过程中的动态力平衡。随后,在“实验-有限元建模”耦合方法的框架下,对一维材料点模拟器进行调整,直至两种方法的结果达到良好一致性。这些数据被用于通过比较MATLAB材料响应模拟器(一维)的力学响应与SHPB有限元模型(一维)试样中心线应力来调整ISV模型的材料常数。此处,有限元模型中试样的应力分量沿波加载方向。然后,通过迭代执行有限元模拟并调整ISV常数,对有限元模型试样的三维行为进行标定,使加载方向上的体积平均应力与实验测得的真实应力-应变响应高度吻合。因此,开展了一套在实验数据、有限元结果与一维ISV材料模型之间进行迭代优化的流程。表1总结了ISV材料模型(MSU TP Ver. 1.1)11的变量。
该方法学中最关键的环节是获取生物材料的一维力学响应及其针对ISV材料模型的材料参数,从而规避了SHPB测试中应力状态不均匀的问题。该方法还能分离出由惯性效应引起的生物材料初始非线性响应,从而获得仅反映材料本征特性的力学响应。耦合方法还表明,试样几何形状的改变会完全改变边界值问题(BVP)以及试样在加载方向上的真实应力-应变响应。因此,上述方法可与任何材料模型(无论是唯象模型还是基于微结构的模型)结合使用,用于标定并在模拟中再现人体器官在损伤性力学载荷下的高应变率行为。
访问受限。请登录或开始试用以查看此内容。
注意:伦理声明:本研究工作符合本机构的研究政策,并严格遵守相应的生物安全和监管合规办公室(ORC)指南。
1. 生物材料样本采集
2. 生物材料样品制备
3. 分离式霍普金森压杆测试
4. SHPB 数据后处理
5. SHPB 有限元建模
访问受限。请登录或开始试用以查看此内容。
耦合方法的有效性在图3中得以体现。在此图中,脑组织的SHPB实验应力-应变响应处于较低的应力状态(峰值应力为0.32 MPa),相较于一维材料点模拟器的应力状态(峰值为0.74 MPa)更低,后者近似于有限元样本中心线(单元)的平均值。这归因于软性生物材料所表现出的变形特性。由于应变率较高,而生物材料的波速和强度较低,材料内部的惯性效应和应力波传播导致了非均匀变形。该现象在样本边缘最为显著,在中心区域最弱。由于试样为圆柱形,其中心区域无法像边缘那样发生非均匀的径向膨胀。因此,在经过一定时间后,样本中心线的变形可被近似视为单轴变形。
由于在经历一段初始的“上升时间”后,可观察到样品中心线呈现单轴变形,因此可利用有限元分析(FEA)提取实验装置无法捕捉的中心线数据。此处的“上升时间”是指在分离式霍普金森压杆(SHPB)实验初始阶段达到应力状态平衡所需的时间。为此,需将虚拟的FEA应变片结果与实验应变片结果进行比较,并调整材料常数,直至两者达到良好一致性。表2给出了通过耦合SHPB实验与有限元模拟方法获得的脑组织代表性材...
访问受限。请登录或开始试用以查看此内容。
报道的方法将SHPB实验与SHPB的有限元(FE)建模相结合,提供了一种新颖且独特的方法,用于评估生物材料在高应变率下的单轴真应力-应变响应。为了获得固有于原生组织的力学性能,在SHPB测试前必须将生物材料试样保持在5.56–7.22 °C之间。如果试样冷却至5.56 °C以下,组织中的水分会开始结晶成冰,从而改变组织的力学性能。尽管其他研究人员15-18出于防止力学性能退化的保存目的而将样品冷冻,但由此SHPB测试获得的结果显示出显著不同的力学响应。此外,Van Ee和Myers的一份报告23表明,死后5小时内测试的软生物材料可获得最佳实验结果。另外,选择磷酸盐缓冲液(PBS)溶液来储存生物材料试样和圆柱形样品,是因为其渗透压和离子浓度与生物体液相似9。
基于 Gray 和 Blumenthal 的研究工作24 在《ASM 软材料高应变率测试手册》中,针对不同类型的生物材料(猪脑、肝、肌腱或脂肪)进行测试时,确定的最佳试样长径比(即试样厚度与直径之比)为 0.5 或更低。Gray 和 Blum...
访问受限。请登录或开始试用以查看此内容。
作者声明,本出版物所涉及的所有材料均不存在利益冲突。
作者谨此感谢密西西比州立大学先进车辆系统中心(CAVS)及农业与生物工程系对本研究工作的支持。本研究工作由美国陆军坦克车辆与装备司令部(TACOM Life Cycle Command)资助,合同号为 W56HZV-08-C-0236,通过密西西比州立大学的分包合同执行,属于“基于仿真的可靠性与安全性”(Simulation Based Reliability and Safety, SimBRS)研究项目的一部分。此外,本研究还得到美国国家核安全管理局(NNSA,隶属于能源部)的资助,资助编号为 [DE-FC26-06NT42755]。最后,作者谨向 David Adams 先生、Michael McCollum 先生以及 Erin Colebeck 女士在本研究中所付出的努力表示衷心感谢。
访问受限。请登录或开始试用以查看此内容。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 高压316不锈钢螺纹管接头,1/2英寸外螺纹 × 1/4英寸内螺纹管径,六角减径衬套 | McMaster-Carr | 2 | |
| 316型不锈钢螺纹管接头,3/4英寸外螺纹 × 1/4英寸内螺纹,六角减径衬套,150 psi | McMaster-Carr | 2 | |
| 易维护型316SS球阀,带316不锈钢端部,1/2" NPT内螺纹 | McMaster-Carr | 2 | |
| 易维护型316SS球阀,带316不锈钢端部,3/4" NPT内螺纹 | McMaster-Carr | 2 | |
| ASME标准不锈钢弹出式安全阀,1/4 NPT外螺纹,300 psi | McMaster-Carr | 2 | |
| 精密超高压316SS管接头,1/2 × 1/2管径,1-7/8"长度,六角短接头 | McMaster-Carr | 8 | |
| 316型不锈钢螺纹管接头,1/2管径,三通,150 psi | McMaster-Carr | 2 | |
| 带安全外壳的测试压力表,聚酯外壳,标准型,干式,600 psi | McMaster-Carr | 2 | |
| 数字压力表,塑料外壳,2-1/2"表盘,1/4底部接口,300 psi | McMaster-Carr | 2 | |
| 316型不锈钢37度扩口管接头,用于1/4"管外径 × 1/8" NPT外螺纹管 | McMaster-Carr | 12 | |
| 303不锈钢37度JIC旋转接头,用于3/16"内径 | McMaster-Carr | 12 | |
| 高压化学软管,3/16"内径,0.312"外径,3,000 psi | McMaster-Carr | 6 | |
| 高纯度单级气体减压器,氮气,0–125 psi,CGA #580 | McMaster-Carr | 2 | |
| 氮气、氩气和氧气用软管,黄铜内螺纹接头,PTFE软管,长度3英尺,1/4"内径,3,600 psi | McMaster-Carr | 2 | |
| [header] | |||
| 超高压316不锈钢螺纹管接头,1/4 × 1/4管径,六角短接头 | McMaster-Carr | 4 | |
| 超高压316不锈钢螺纹管接头,3/4 × 3/4管径,六角短接头 | McMaster-Carr | 2 | |
| 超高压316不锈钢螺纹管接头,1/4英寸外螺纹 × 1/8英寸内螺纹管径,六角衬套 | McMaster-Carr | 2 | |
| 标准黄铜压缩管接头转换头,用于1/4"管外径 × 1/4" NPTF外螺纹管 | McMaster-Carr | 4 | |
| Kobalt 1/4英寸微型调节器带压力表 | Lowes | 2 | |
| 1/4" × 25英尺聚乙烯管 | Lowes | 2 | |
| 直径1-1/2"聚碳酸酯(PC)棒材 | McMaster-Carr | 2 | |
| LTV-35 四通阀,Mead流体动力 | Motion Industries | 2 | |
| 气动双作用执行器 | Valtronic | 2 | |
| 不锈钢球阀 1/2" | Valtronic | 2 | |
| Buckeye压力容器 | Buckeye | 2 | |
| SR-4通用FAE-25-35SX应变片 | Micro-Measurement Vishay Precision Group | 2 | |
| M-M信号调理放大器2310A | Micro-Measurement Vishay Precision Group | 1 | |
| 激光ROLS-W光学传感器 | Monarch Instruments | 1 |
访问受限。请登录或开始试用以查看此内容。