Name: 一个耦合试验，有限元建模方法以评估生物材料软的高应变率力学响应 (视频) | JoVE Journal
Uploaded: 2015-05-18T13:00:00.000+00:00
Duration: 11 min 28 s
Description: 了解耦合试验与有限元建模方法，用于评估脑组织和肝脏等软生物材料在高应变率下的动态力学响应。该方法结合分离式霍普金森压杆测试与迭代优化仿真，提高材料模型精度。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

playerTrialDesktop

作者谨此感谢密西西比州立大学先进车辆系统中心（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 女士在本研究中所付出的努力表示衷心感谢。

注意：伦理声明：本研究工作符合本机构的研究政策，并严格遵守相应的生物安全和监管合规办公室（ORC）指南。
1. 生物材料样本采集
<ol>
	<li>根据实验室和/或机构的标准生物安全规程穿戴个人防护装备。在处理猪组织和进行测试时，需穿着包趾鞋、长裤，佩戴实验服、外科手套、防护口罩和安全护目镜。</li>
	<li>从当地屠宰场获取健康猪只死后1至2小时内取出的猪组织（头部、腹部或后腿）。</li>
	<li>将猪组织存放于生物危害安全袋中，然后置于冰容器中（约5.56–7.22 °C）。 
	注意：使用温度计确认猪组织样本温度不低于7.22 °C。</li>
	<li>将猪组织运送至最近的实验室（密西西比州立大学兽医学院）进行解剖。</li>
	<li>在兽医学院兽医的监督下，手术取出猪器官（脑、肝、肌肉、脂肪或肌腱），并将其放入盛有磷酸盐缓冲液（PBS）的容器中临时保存（pH 7.4）。</li>
	<li>将装有PBS的容器存放在冰盒中（约5.56–7.22 °C），并立即运送至测试场所进行样品制备和SHPB测试。</li>
</ol>
2. 生物材料样品制备
<ol>
	<li>从 PBS 容器中取出猪器官，并将其放置在无菌表面上。</li>
	<li>注意：确定主要纤维取向以及每个测试样本的取样位置。使用内径为 30 mm 的圆柱形冲模从猪器官上切取测试样本。</li>
	<li>如果测试样本卡在圆柱形冲模内部，通过解剖工具的另一端注入 PBS，使测试样本完整滑出。将取出的测试样本放置在无菌表面的另一区域。</li>
	<li>使用手术刀将样本修整至规定的厚度和长径比。 
	注意：对于猪组织样本的 SHPB 测试，样本厚度为 10–15 mm，长径比（厚度/直径）为 0.33–0.50（图 2）。</li>
	<li>使用游标卡尺在三个不同位置测量样本的厚度和直径。</li>
	<li>将所有测试样本保存在新鲜 PBS 中，直至 SHPB 设备准备就绪进行测试。 
	注意：确保样本在屠宰后 4 小时内完成测试。</li>
	<li>因切口误差或横截面不规则而无法形成圆柱形的样本应予弃用，并将废弃样本放入生物危害安全袋中。重复步骤 2.2–2.6 以获取额外的测试样本。</li>
</ol>
3. 分离式霍普金森压杆测试
<ol>
	<li>将冲杆、入射杆和透射杆放置在SHPB测试用的金属支架中。 
	注意：确保各杆件可自由移动，且其接触面相互对齐。为透射杆设置止动装置以确保安全。</li>
	<li>将粘贴在入射杆和透射杆上的应变片连接至信号放大器。打开信号调理放大器和数据采集模块计算机。</li>
	<li>启动高速数据采集软件。</li>
	<li>检查信号的实时采集情况，确认其处于正常范围内，并通过点击归零图标消除噪声信号。</li>
	<li>输入触发水平和数据采集速率（2 MHz）。</li>
	<li>设置软件在达到触发水平后开始记录数据。</li>
	<li>将冲杆置于压力腔附近。将压力腔充至所需压力。 
	注意：典型压力范围为5–25 psi。</li>
	<li>按下归零按钮对激光测速仪归零，并通过在冲杆后方激光传感器位置设置反射条，使测速仪读取冲杆速度。</li>
	<li>放置样品约束腔，确保其不妨碍入射杆和反射杆的运动。使入射杆与透射杆接触。</li>
	<li>为校准目的，在冲杆的压力腔上打开触发开关，进行一次无样品的测试。</li>
	<li>数据采集完成后，保存并分析SHPB应变片数据（下节将讨论），以确认测试流程正常运行。</li>
	<li>将圆柱形样品置于入射杆和透射杆之间，然后关闭样品约束腔。 
	注意：确保未对样品进行任何预处理。</li>
	<li>在样品置于入射杆和透射杆之间的情况下，执行步骤3.4–3.7的操作。 
	注意：确保样品中心线与杆件中心线一致。在继续操作前，还需检查样品未受压，且保持与取出时相同的几何形态。</li>
	<li>测试完成后，使用一次性卫生擦拭巾清除入射杆、透射杆及样品约束腔中的样品残留物。将所有残留物和擦拭巾丢弃至生物危害安全袋中。</li>
	<li>使用70%乙醇清洁液和卫生擦拭巾对杆件及样品约束腔进行消毒。</li>
</ol>
4. SHPB 数据后处理
<ol>
	<li>打开“MSU High Rate Software19”以分析霍普金森杆波形。</li>
	<li>启动软件后，先进入设置窗口，在模式选项卡中选择“拉伸/压缩”模式以进行单轴测试；同时在应变片选项卡中选择“2个应变片”，然后点击“继续”。</li>
	<li>在主窗口中，选择“打开文件1”选项卡，导航至入射杆上应变片记录的入射波数据；再选择“打开文件2”选项卡以导入透射杆上的应变片记录。</li>
	<li>在主窗口中选择“参数”选项卡，输入试验装置的物理参数，包括：杆的尺寸、电压-应变转换系数、应变片位置以及粘弹性弥散常数，然后点击“继续”。</li>
	<li>接着在主窗口中选择“选择数据”选项卡，使用光标条将数据集截取为仅包含入射波、反射波和透射波的有效部分，然后点击“继续”。</li>
	<li>然后在主窗口中选择“选择波形”选项卡，使用光标条在“入射波形图”中限定入射波，在“反射波形图”中限定反射波，在“透射波形图”中限定透射波，完成后点击“继续”。</li>
	<li>之后在主窗口中选择“校正”选项卡，使软件对粘弹性弥散20-21进行校正。</li>
	<li>现在选择主窗口中的“平移”选项卡，在波形图中，通过“波形选择”选项卡依次单独选中入射波、反射波和透射波，并使用光标将其拖动至相同的时间起始位置；在“数据图”中查看所有波形对齐情况；完成后点击“继续”。</li>
	<li>在结果文件中，点击“另存为”以保存载荷、位移、位置和速度曲线。</li>
	<li>使用 Microsoft Excel（或其他电子表格软件）中的常规方法，结合霍普金森杆试验前测量的试样尺寸，计算真实应力和真实应变。</li>
</ol>
5. SHPB 有限元建模
<ol>
	<li>使用商用有限元（FE）软件，建立SHPB装置的有限元模型。 
	注意：使用相同的几何形状和材料属性。</li>
	<li>为冲杆的有限元模型赋予初始速度，以启动有限元仿真。 
	注意：冲杆的速度应与特定应变率下SHPB实验中的冲杆速度相对应9。</li>
	<li>建立一个在入射杆和透射杆之间未放置试样的SHPB装置有限元模型，并运行有限元仿真。 
	注意：模拟的冲杆速度应与“无试样”条件下实验中的冲杆速度一致。为聚合物杆指定表1中给出的材料属性。</li>
	<li>验证实验与有限元仿真中应变片测量结果（应变随时间变化）是否具有良好一致性。</li>
	<li>将生物材料试样纳入SHPB装置的有限元模型中，并为该生物材料试样分配ISV材料模型的三维实现（以vumat文件格式22）11。</li>
	<li>采用三种不同的网格尺寸进行网格细化研究，然后分析结果以确定解是否收敛。 
	注意：网格尺寸对应于构成有限元模型的六面体和/或四面体单元的总数。在进一步仿真后，选择最低网格尺寸且解已收敛的有限元模型9。</li>
	<li>执行两步有限元模型校准。第一步，将实验数据导入ISV材料模型的一维实现中。</li>
	<li>通过调整ISV材料模型的参数（见表1），使实验的真实应力-应变曲线与模型的真实应力-应变曲线相匹配。 
	注意：由于实验SHPB数据本质上是三维的，而材料模型为一维，因此需要进一步迭代。</li>
	<li>将ISV材料常数赋给SHPB装置有限元模型中的生物材料试样。</li>
	<li>运行有限元仿真，使冲杆速度和试样变形应变率与相同应变率下的SHPB测试相对应。</li>
	<li>比较实验与有限元仿真中的应变片测量结果，确保两者具有良好一致性（应变随时间变化）。 
	注意：如果有限元仿真与实验的应变片数值具有良好一致性，则进入模型校准过程的第二步；否则重复任务5.7–5.11。</li>
	<li>在有限元模型校准的第二步中，运行有限元仿真，并使用SHPB实验后处理软件MSU High Rate软件19–21处理应变片数据。 
	注意：如果模拟的真实应力-应变响应与实验的真实应力-应变响应相符，则两步有限元模型校准已完成；否则重复任务5.7–5.12。</li>
	<li>对有限元模型试样中心线单元的加载方向（Σ33）应力进行体积平均。 
	注意：如果该应力与一维ISV材料模型结果的应力-应变曲线具有良好一致性，则通过任务5.7–5.12获得的结果已完全校准；否则重复任务5.7–5.13。通过ISV材料模型一维实现所获取的真实应力-应变响应，代表了在SHPB装置中测试的生物材料的单轴真实应力-应变响应。</li>
</ol>

<ol><li>Champion, H. R., Holcomb, J. B., Young, L. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Injuries+from+explosions%3A+physics%2C+biophysics%2C+pathology%2C+and+required+research+focus%5BTitle%5D)+AND+%22J+Trauma%22%5BJournal%5D)">Injuries from explosions: physics, biophysics, pathology, and required research focus.</a> J Trauma. 66 (5), 1468-1477 (2009).</li>
<li>Aubry, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Summary+and+agreement+statement+of+the+First+International+Conference+on+Concussion+in+Sport%2C+Vienna+2001.+Recommendations+for+the+improvement+of+safety+and+health+of+athletes+who+may+suffer+concussive+injuries%5BTitle%5D)+AND+%22Br+J+Sports+Med%22%5BJournal%5D)">Summary and agreement statement of the First International Conference on Concussion in Sport, Vienna 2001. Recommendations for the improvement of safety and health of athletes who may suffer concussive injuries.</a> Br J Sports Med. 36 (1), 6-10 (2002).</li>
<li>Born, C. T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Blast+trauma%3A+the+fourth+weapon+of+mass+destruction%5BTitle%5D)+AND+%22Scand+J+Surg%22%5BJournal%5D)">Blast trauma: the fourth weapon of mass destruction.</a> Scand J Surg. 94 (4), 279-285 (2005).</li>
<li>Cullis, I. G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Blast+waves+and+how+they+interact+with+structures%5BTitle%5D)+AND+%22J+R+Army+Med+Corps%22%5BJournal%5D)">Blast waves and how they interact with structures.</a> J R Army Med Corps. 147, 16-26 (2001).</li>
<li>Ngo, T., Mendis, P., Gupta, A., Ramsay, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Blast+Loading+and+Blast+Effects+on+Structures%E2%80%93An+Overview%5BTitle%5D)+AND+%22Electronic+Journal+of+Structural+Engineering%22%5BJournal%5D)">Blast Loading and Blast Effects on Structures–An Overview.</a> Electronic Journal of Structural Engineering. 7, 76-91 (2007).</li>
<li>Usmani, Z. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Intelligent+Agents+in+Extreme+Conditions+%E2%80%93+Modeling+and+Simulation+of+Suicide+Bombing+for+Risk+Assessment%5BTitle%5D)+AND+%22Web+Intelligence+and+Intelligent+Agents%22%5BJournal%5D)">Intelligent Agents in Extreme Conditions – Modeling and Simulation of Suicide Bombing for Risk Assessment.</a> Web Intelligence and Intelligent Agents. , (2010).</li>
<li>Guskiewicz, K. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Cumulative+effects+associated+with+recurrent+concussion+in+collegiate+football+players+the+NCAA+Concussion+Study%5BTitle%5D)+AND+%22JAMA%22%5BJournal%5D)">Cumulative effects associated with recurrent concussion in collegiate football players the NCAA Concussion Study.</a> JAMA. 290 (19), 2549-2555 (2003).</li>
<li>Finkelstein, E., Corso, P., Miller, T. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aFinkelstein%2C+E+%22The+Incidence+and+Economic+Burden+of+Injuries+in+the+United+States%22">The Incidence and Economic Burden of Injuries in the United States</a>. , Oxford University Press. New York (NY). (2006).</li>
<li>Prabhu, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Coupled+experiment%2Ffinite+element+analysis+on+the+mechanical+response+of+porcine+brain+under+high+strain+rates%5BTitle%5D)+AND+%22JMech+Behav+Biomed+Mater%22%5BJournal%5D)">Coupled experiment/finite element analysis on the mechanical response of porcine brain under high strain rates.</a> JMech Behav Biomed Mater. 4 (7), 1067-1080 (2011).</li>
<li>Horstemeyer, M. F. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aHorstemeyer%2C+MF+%22Integrated+Computational+Materials+Engineering+%28ICME%29%3A+Using+Multiscale+Modeling+to+Invigorate+Engineering+Design+with+Science%22">Integrated Computational Materials Engineering (ICME): Using Multiscale Modeling to Invigorate Engineering Design with Science</a>. , Wiley Press. (2012).</li>
<li>Bouvard, J. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+general+inelastic+internal+state+variable+model+for+amorphous+glassy+polymers%5BTitle%5D)+AND+%22Acta+Mechanica%22%5BJournal%5D)">A general inelastic internal state variable model for amorphous glassy polymers.</a> Acta Mechanica. 213, 1-2 (2010).</li>
<li>Kenner, V. H., Goldsmith, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Impact+on+a+simple+physical+model+of+the+head%5BTitle%5D)+AND+%22J+Biomech%22%5BJournal%5D)">Impact on a simple physical model of the head.</a> J Biomech. 6 (1), 1-11 (1973).</li>
<li>Khalil, T. B., Viano, D. C., Smith, D. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Experimental+analysis+of+the+vibrational+characteristics+of+the+human+skull%5BTitle%5D)+AND+%22J.+Sound+Vib%22%5BJournal%5D)">Experimental analysis of the vibrational characteristics of the human skull.</a> J. Sound Vib. 63 (3), 351-376 (1979).</li>
<li>Pervin, F., Chen, W. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Dynamic+mechanical+response+of+bovine+gray+matter+and+white+matter+brain+tissues+under+compression%5BTitle%5D)+AND+%22J+Biomech%22%5BJournal%5D)">Dynamic mechanical response of bovine gray matter and white matter brain tissues under compression.</a> J Biomech. 42 (6), 731-735 (2009).</li>
<li>Prevost, T. P., Balakrishnan, A., Suresh, S., Socrate, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Biomechnics+of+brain+tissue%5BTitle%5D)+AND+%22Acta+Biomater%22%5BJournal%5D)">Biomechnics of brain tissue.</a> Acta Biomater. 7 (1), 83-95 (2011).</li>
<li>Saraf, H., Ramesh, K. T., Lennon, A. M., Merkle, A. C., Roberts, J. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Mechanical+properties+of+soft+human+tissues+under+dynamic+loading.J%5BTitle%5D)+AND+%22J+Biomech%22%5BJournal%5D)">Mechanical properties of soft human tissues under dynamic loading.J.</a> J Biomech. 40 (9), 1960-1967 (2007).</li>
<li>Van Sligtenhorst, C., Cronin, D. S., Wayne Brodland, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((High+strain+rate+compressive+properties+of+bovine+muscle+tissue+determined+using+a+split+Hopkinson+bar+apparatus%5BTitle%5D)+AND+%22J+Biomech%22%5BJournal%5D)">High strain rate compressive properties of bovine muscle tissue determined using a split Hopkinson bar apparatus.</a> J Biomech. 39 (10), 1852-1858 (2006).</li>
<li>Song, B., Chen, W., Ge, Y., Weerasooriya, Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Dynamic+and+quasi-static+compressive+response+of+porcine+muscle%5BTitle%5D)+AND+%22J+Biomech%22%5BJournal%5D)">Dynamic and quasi-static compressive response of porcine muscle.</a> J Biomech. 40 (13), 2999-3005 (2007).</li>
<li>MSU JHBT Data Processing and MSU High Rate Software Manual. , Available from: 
 <a target="_blank" href="https://icme.hpc.msstate.edu/mediawiki/index.php/File:MSU_JHBT_Data_Processing_and_MSU_High_Rate_Software_Manual.zip">https://icme.hpc.msstate.edu/mediawiki/index.php/File:MSU_JHBT_Data_Processing_and_MSU_High_Rate_Software_Manual.zip</a> (2014).</li>
<li>Zhao, H., Gary, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((On+the+use+of+SHPB+techniques+to+determine+the+dynamic+behavior+of+materials+in+the+range+of+small+strains%5BTitle%5D)+AND+%22Int+J+Solids+Struct%22%5BJournal%5D)">On the use of SHPB techniques to determine the dynamic behavior of materials in the range of small strains.</a> Int J Solids Struct. 33 (23), 3363-3375 (1996).</li>
<li>Zhao, H., Gary, G., Klepaczko, J. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((On+the+use+of+a+viscoelastic+split+hopkinson+pressure+bar%5BTitle%5D)+AND+%22Int+J+Impact+Eng%22%5BJournal%5D)">On the use of a viscoelastic split hopkinson pressure bar.</a> Int J Impact Eng. 19 (4), 319-330 (1997).</li>
<li>MSU TP Ver 1.1.. , Available from: 
 <a target="_blank" href="https://icme.hpc.msstate.edu/mediawiki/index.php/File:MSU_TP_Ver_1.1.zip">https://icme.hpc.msstate.edu/mediawiki/index.php/File:MSU_TP_Ver_1.1.zip</a> (2014).</li>
<li>Gray, G. T., Blumenthal, W. R. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aGray%2C+GT+%22ASM+Handbook%2C+Mechanical+Testing+and+Evaluation%22">ASM Handbook, Mechanical Testing and Evaluation</a>. 8, ASM International. 488-496 (2000).</li>
<li>Dharan, C. K. H., Hauser, F. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Determination+of+stress-strain+characteristics+at+very+high+strain+rates%5BTitle%5D)+AND+%22Exp.+Mech%22%5BJournal%5D)">Determination of stress-strain characteristics at very high strain rates.</a> Exp. Mech. 10 (9), 370-376 (1970).</li>
<li>Chen, J., Priddy, L. B., Prabhu, R., Marin, E. B., Horstemeyer, M. F., Williams, L. N., Liao, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Traumatic+Injury%3A+Mechanical+Response+of+Porcine+Liver+Tissue+under+High+Strain+Rate+Compression+Testing%5BTitle%5D)+AND+%22Proceedings+of+the+ASME+2009+Summer+Bioengineering+Conference+%28SBC2009%29%22%5BJournal%5D)">Traumatic Injury: Mechanical Response of Porcine Liver Tissue under High Strain Rate Compression Testing.</a> Proceedings of the ASME 2009 Summer Bioengineering Conference (SBC2009). , Resort at Squaw Creek. Lake Tahoe, CA, USA. (2009).</li>
</ol>

作者声明，本出版物所涉及的所有材料均不存在利益冲突。

报道的方法将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...

动机
软生物材料（如脑、肝、肌腱、脂肪）的耦合分离式霍普金森压杆（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）材料模型背景介绍
由于涉及伦理问题 在&#160;体内 人体器官测试及大规模人体遗骸测试相关的物流问题，当前研究工作涉及机械实验 在&#160;体外 使用从动物替代品中提取的器官制备的标本（例如，猪作为最常用的替代模型）。聚合物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）以及试样在加载方向上的真实应力-应变响应。因此，上述方法可与任何材料模型（无论是唯象模型还是基于微结构的模型）结合使用，用于标定并在模拟中再现人体器官在损伤性力学载荷下的高应变率行为。

<table><tbody><tr><td>高压316不锈钢螺纹管接头，1/2英寸外螺纹 × 1/4英寸内螺纹管径，六角减径衬套</td><td>McMaster-Carr</td><td></td><td>2</td></tr><tr><td>316型不锈钢螺纹管接头，3/4英寸外螺纹 × 1/4英寸内螺纹，六角减径衬套，150&amp;nbsp;psi</td><td>McMaster-Carr</td><td></td><td>2</td></tr><tr><td>易维护型316SS球阀，带316不锈钢端部，1/2&quot; NPT内螺纹</td><td>McMaster-Carr</td><td></td><td>2</td></tr><tr><td>易维护型316SS球阀，带316不锈钢端部，3/4&quot; NPT内螺纹</td><td>McMaster-Carr</td><td></td><td>2</td></tr><tr><td>ASME标准不锈钢弹出式安全阀，1/4 NPT外螺纹，300&amp;nbsp;psi</td><td>McMaster-Carr</td><td></td><td>2</td></tr><tr><td>精密超高压316SS管接头，1/2 × 1/2管径，1-7/8&quot;长度，六角短接头</td><td>McMaster-Carr</td><td></td><td>8</td></tr><tr><td>316型不锈钢螺纹管接头，1/2管径，三通，150&amp;nbsp;psi</td><td>McMaster-Carr</td><td></td><td>2</td></tr><tr><td>带安全外壳的测试压力表，聚酯外壳，标准型，干式，600&amp;nbsp;psi</td><td>McMaster-Carr</td><td></td><td>2</td></tr><tr><td>数字压力表，塑料外壳，2-1/2&quot;表盘，1/4底部接口，300&amp;nbsp;psi</td><td>McMaster-Carr</td><td></td><td>2</td></tr><tr><td>316型不锈钢37度扩口管接头，用于1/4&quot;管外径 × 1/8&quot; NPT外螺纹管</td><td>McMaster-Carr</td><td></td><td>12</td></tr><tr><td>303不锈钢37度JIC旋转接头，用于3/16&quot;内径</td><td>McMaster-Carr</td><td></td><td>12</td></tr><tr><td>高压化学软管，3/16&quot;内径，0.312&quot;外径，3,000 psi</td><td>McMaster-Carr</td><td></td><td>6</td></tr><tr><td>高纯度单级气体减压器，氮气，0–125 psi，CGA #580</td><td>McMaster-Carr</td><td></td><td>2</td></tr><tr><td>氮气、氩气和氧气用软管，黄铜内螺纹接头，PTFE软管，长度3英尺，1/4&quot;内径，3,600 psi</td><td>McMaster-Carr</td><td></td><td>2</td></tr><tr><td colspan="4">[header]</td></tr><tr><td>超高压316不锈钢螺纹管接头，1/4 × 1/4管径，六角短接头</td><td>McMaster-Carr</td><td></td><td>4</td></tr><tr><td>超高压316不锈钢螺纹管接头，3/4 × 3/4管径，六角短接头</td><td>McMaster-Carr</td><td></td><td>2</td></tr><tr><td>超高压316不锈钢螺纹管接头，1/4英寸外螺纹 × 1/8英寸内螺纹管径，六角衬套</td><td>McMaster-Carr</td><td></td><td>2</td></tr><tr><td>标准黄铜压缩管接头转换头，用于1/4&quot;管外径 × 1/4&quot; NPTF外螺纹管</td><td>McMaster-Carr</td><td></td><td>4</td></tr><tr><td>Kobalt 1/4英寸微型调节器带压力表</td><td>Lowes</td><td></td><td>2</td></tr><tr><td>1/4&quot; × 25英尺聚乙烯管</td><td>Lowes</td><td></td><td>2</td></tr><tr><td>直径1-1/2&quot;聚碳酸酯（PC）棒材</td><td>McMaster-Carr</td><td></td><td>2</td></tr><tr><td>LTV-35 四通阀，Mead流体动力</td><td>Motion Industries</td><td></td><td>2</td></tr><tr><td>气动双作用执行器</td><td>Valtronic</td><td></td><td>2</td></tr><tr><td>不锈钢球阀 1/2&quot;</td><td>Valtronic</td><td></td><td>2</td></tr><tr><td>Buckeye压力容器</td><td>Buckeye</td><td></td><td>2</td></tr><tr><td>SR-4通用FAE-25-35SX应变片</td><td>Micro-Measurement Vishay Precision Group</td><td></td><td>2</td></tr><tr><td>M-M信号调理放大器2310A</td><td>Micro-Measurement Vishay Precision Group</td><td></td><td>1</td></tr><tr><td>激光ROLS-W光学传感器</td><td>Monarch Instruments</td><td></td><td>1</td></tr></tbody></table>

一种用于评估软生物材料高应变率力学响应的耦合实验-有限元建模方法

A Coupled Experiment-finite Element Modeling Methodology for Assessing High Strain Rate Mechanical Response of Soft Biomaterials

本研究提供了一种结合实验与有限元（FE）模拟的方法，用于研究软生物材料的力学行为例如 脑、肝、肌腱、脂肪 等。) 在受到高应变率作用时。本研究采用分离式霍普金森压杆（Split-Hopkinson Pressure Bar, SHPB）产生 100–1,500 sec⁻¹ 的应变率。-1. 采用由粘弹性材料（聚碳酸酯）制成的撞击杆的分离式霍普金森压杆（SHPB）装置。在死亡后不久获取生物材料样本，并制备用于SHPB测试。将试样置于入射杆和透射杆之间，启动SHPB的气动组件，驱动撞击杆撞击入射杆。由此产生的冲击生成一个压缩应力波（即 入射波）穿过入射杆。当压缩应力波到达入射杆末端时，一部分继续向前穿过样品和透射杆（即 透射波），而另一部分则在入射杆中以拉伸波的形式反向传播（即 反射波）。这些波通过安装在入射杆和透射杆上的应变片进行测量。试样的真实应力-应变行为是基于波传播理论和动态力平衡方程计算得到的。由于试样在加载过程中发生鼓胀，实验测得的应力-应变响应本质上具有三维特性，因此采用静水应力（第一不变量）来构建应力-应变关系。为了提取组织的单轴（一维）力学响应，结合实验结果与包含组织内部状态变量（ISV）材料模型的有限元分析（FEA），进行了迭代耦合优化。该优化过程中，用于模拟实验装置的有限元仿真中的ISV材料模型被反复校准（即 优化后的实验数据与有限元分析（FEA）的应变片数值及应力的第一不变量吻合良好。

耦合方法的有效性在图3中得以体现。在此图中，脑组织的SHPB实验应力-应变响应处于较低的应力状态（峰值应力为0.32 MPa），相较于一维材料点模拟器的应力状态（峰值为0.74 MPa）更低，后者近似于有限元样本中心线（单元）的平均值。这归因于软性生物材料所表现出的变形特性。由于应变率较高，而生物材料的波速和强度较低，材料内部的惯性效应和应力波传播导致了非均匀变形。该现象在样本边缘最为显著，在中心区域最弱。由于试样为圆柱形，其中心区域无法像边缘那样发生非均匀的径向膨胀。因此，在经过一定时间后，样本中心线的变形可被近似视为单轴变形。
由于在经历一段初始的“上升时间”后，可观察到样品中心线呈现单轴变形，因此可利用有限元分析（FEA）提取实验装置无法捕捉的中心线数据。此处的“上升时间”是指在分离式霍普金森压杆（SHPB）实验初始阶段达到应力状态平衡所需的时间。为此，需将虚拟的FEA应变片结果与实验应变片结果进行比较，并调整材料常数，直至两者达到良好一致性。表2给出了通过耦合SHPB实验与有限元模拟方法获得的脑组织代表性材...

本研究提出一种耦合实验与有限元模拟的方法，用于获取软生物材料（脑、肝、肌腱、脂肪）的单轴动态力学响应 等等。通过生物材料有限元分析的迭代优化，将分离式霍普金森压杆试验中因试样鼓胀而获得的多轴实验结果转化为单轴真应力-应变行为。

本研究提供了一种结合实验与有限元（FE）模拟的方法，用于研究软生物材料（例如……）的力学行为。

本实验的总体目标是获取被测软生物材料在单轴高应变率条件下的真实应力-应变响应。首先，利用分离式霍普金森压杆（SHPB）装置对脑组织样品进行测试，然后根据实验获得的真实应力-应变响应，标定材料模型的一维形式。随后，对SHPB实验装置进行有限元（FE）仿真，以评估样品在平均加载方向上的应力及中心线应力。最后一步是迭代比较一维材料模型拟合结果与有限元模拟的中心线应力，以及三维有限元加载方向的应力与实验测得的真实应力-应变响应，直至上述结果之间达到良好一致性。最终，通过耦合实验与有限元方法表明，可利用分离式霍普金森压杆（SHPB）获取软生物材料的单轴真实应力-应变响应。相较于冲击管试验等现有方法，该技术的主要优势在于能够表征软生物材料在高应变率下的单轴力学响应。该技术的应用前景包括通过虚拟人体建模来诊断创伤性脑损伤。因为当前方法获得的数据能够准确捕捉由头部外伤引起局部损伤和组织破坏。我们最初产生该方法的想法是在将脑组织高应变率实验数据与甲状腺实验及建模结果进行关联研究时提出的。本次实验操作将由兽医 Dr. Jim Cooley、研究生 Mr. Soro Patnaik 和研究助理 Mr. Wilburn Wittington 演示。在兽医的监督下，外科手术提取猪的脑、肝、肌肉、脂肪或肌腱等器官。然后将其放入盛有磷酸盐缓冲液（PBS）的容器中进行临时保存。将装有PBS的容器置于冰盒中，并立即运送至检测实验室，以进行样品制备和分离式霍普金森压杆（SHPB）测试。到达检测实验室后，从PBS容器中取出猪器官，放置于无菌表面。 确定每个测试样本的主要纤维取向和位置。然后使用内径为30毫米的圆柱形染色工具从猪器官上切取测试样本。接下来，使用支架将样本修整至猪组织SHPB测试规定的厚度和长宽比。厚度为10至15毫米，而长径比为0.33至0.5。使用卡尺在三个不同位置测量厚度和直径。将撞击杆、入射杆和透射杆放置于金属延伸部分中，用于SHPV测试。确保入射杆和透射杆用手触摸时可自由移动，且两者的端面相互对齐。为安全起见，需为透射杆配备止动装置。将应变片连接至入射杆和透射杆，并将其接至信号放大器。打开信号调理放大器和数据采集模块计算机，然后初始化高速数据采集软件。检查信号的实时采集情况，确认其是否处于正常范围内，并通过点击归零图标消除噪声信号。将冲头杆安装至压力室旁边。打开氮气罐的喷嘴，然后打开压力室的进气阀，使压力室达到所需压力。按下归零按钮，将激光测速仪归零，并通过设置冲头杆上的反射条，使测速仪读取冲头杆的速度。反射条位于激光传感器后方。放置样品约束腔，确保其不妨碍入射杆和透射杆的运动。然后将入射杆与透射杆接触，以进行校准。通过打开Stryker杆上压力腔的触发开关，运行一次无样品的测试。校准完成后，将圆柱形样品置于入射杆和透射杆之间，然后关闭样品 confinement 腔室。确

a-coupled-experiment-finite-element-modeling-methodology-for

Research

科研

JoVE Journal

Bioengineering

生物工程

一个耦合试验，有限元建模方法以评估生物材料软的高应变率力学响应

讨论了一种方法，通过该方法 体内 刺激响应材料的力学行为随时间变化被监测。样品进行测试 离体 使用带有环境控制的微拉伸测试仪模拟生理环境。这项工作进一步促进了对 体内...

环境可控的微拉伸测试用于力学自适应聚合物纳米复合材料的研究 离体 表征

着地过程中，下肢骨骼会承受较大的机械载荷并发生变形。测量骨骼变形对于深入理解与冲击相关的骨应力性损伤机制至关重要。本研究采用一种新颖的方法，结合个体特异性肌肉骨骼建模与有限元分析，用于测量动态运动过程中胫骨的应变。...

用于研究动态运动中骨应变的特定受试者肌肉骨骼模型

本文介绍了一种纯单轴加载装置的构建方案。通过采用关键的设计要素，确保获得准确且可重复的测试结果。...

单轴加载机器开发中的设计要素应用

本文介绍了一种利用高速伺服液压加载框架在中等应变率下对拉伸试样进行动态表征的方法。同时，还规定了试样上应变片的安装与分析，以及数字图像相关应变测量的实验步骤。...

使用数字图像相关法进行中等应变率下的材料表征

有必要确定冠状血管中哪些动脉粥样硬化病变会进展，以便在心肌梗死发生前指导干预措施。本文概述了利用商业有限元求解器中的流固耦合技术，基于光学相干断层成像数据对动脉进行生物力学建模的方法，以帮助预测此类病变的进展。...

基于光学相干断层扫描的冠状动脉粥样硬化进展生物力学流固耦合分析

本方案介绍了一种用于软材料纳米压痕实验的完整工作流程，包括水凝胶和细胞。首先，详细阐述了获取力谱数据的实验步骤；随后，通过一款新开发的开源 Python 软件对这类数据进行分析，该软件可从 GitHub...

软物质纳米压痕的实验与数据分析流程

本研究展示了一种用于材料变形与失效的三维虚拟仿真实验，可提供可视化的实验过程。通过一系列实验，用户能够在沉浸式与交互式的学习环境中熟悉设备并掌握相关操作。...

基于扫描电子显微镜的力学虚拟仿真实验：材料变形与失效

本方案描述了利用有限元模拟研究波纹锥形管准静态压缩性能的方法。研究了壁厚梯度对压缩性能的影响。结果表明，合理的壁厚梯度设计可以改变管件的变形模式，并显著提高其能量吸收性能。...

有限元建模模拟波纹锥形管准静态压缩

基于血管内超声图像的有限元建模方法用于量化人冠状动脉的体内力学特性

在本方案中，我们演示了该技术所基于的垂直位移微柱阵列微致动器的制备方法，并介绍了如何对该基础技术进行改进，以在二维和三维培养模式下实现高通量的力学动态细胞培养。...

用于力学动态细胞培养的微加工平台

了解耦合试验与有限元建模方法，用于评估脑组织和肝脏等软生物材料在高应变率下的动态力学响应。该方法结合分离式霍普金森压杆测试与迭代优化仿真，提高材料模型精度。

一种用于评估软生物材料高应变率力学响应的耦合实验-有限元建模方法

探索更多视频