应用于含硬质基底粘弹性薄板的粘弹性测量方法及粘弹性测量系统技术方案

本发明专利技术涉及一种应用于含硬质基底粘弹性薄板的粘弹性测量方法及粘弹性测量系统，该方法包括：利用不同频率的外部激励在待测材料中激励出导波；利用预设成像方法记录导波的传播信息，确定近距范围及获得平均测量距离；根据近距范围内导波的传播信息，计算出近距区域内波速与频率的关系，以给出频散曲线；根据频散曲线和测量范围确定临界频率，并根据表面波模型拟合计算出待测材料的粘弹性。本发明专利技术应用于含硬质基底粘弹性薄板的粘弹性测量方法，利用近距区域内的粘弹性表面波模型，获得薄板的粘弹性。相比于兰姆波法，降低了对测量频率的要求；又因为表面波模型的简单，降低了对测量数据后处理的要求。

【技术实现步骤摘要】
应用于含硬质基底粘弹性薄板的粘弹性测量方法及粘弹性测量系统
本专利技术涉及一种应用于含硬质基底粘弹性薄板的粘弹性测量方法及粘弹性测量系统。
技术介绍
组织病变通常会导致生物组织的力学特性发生改变。因此，通过弹性成像方法在体测量生物组织的力学特性，对某些重大疾病的诊断监控以及药物疗效评价都具有重要意义。对于关节软骨结构，临床上希望有无创定量的力学特性诊断技术。总结现有的软骨弹性测量方法，从原理上看大致有如下两种：1)基于应力应变法的弹性测量原理是测量软骨上的应力σ和应变ε，然后根据应力应变关系E＝σ/ε计算弹性模量E。该方法的难点在于，难以无创的测量软骨应力。在具体实现上，要么在关节软骨上开微小的窗口，从而有创(或微创)的测量软骨受力，如基于印压原理的软骨弹性测量；要么只测量应变，实现软骨弹性的定性测量，如基于核磁的关节软骨应变弹性成像。因此，基于应力应变法无法定量无创测量关节软骨的弹性。2)基于波速法的弹性测量原理是测量组织中的波速，然后根据波速与组织弹性的对应关系(即波传播规律)，计算组织弹性。目前基于波速法的超声弹性仪器，已经可以定量无创测量肝脏、乳腺、血管等生物组织的弹性。肝脏和乳腺组织的尺寸较大，可以近似为三维结构，波传播规律满足剪切波模型；血管壁较薄，周边组织的弹性较低，血管壁可以近似为薄壁结构，波传播规律满足单层板中的兰姆波模型。但关节软骨是有别于肝脏和血管的特殊结构。软骨的厚度较薄，与软骨相连的骨弹性模量远高于软骨。因此从力学上看，关节软骨是一种双层粘弹性板结构，且具有硬质基底。目前已可以基于双层板中的兰姆波模型，测量离体牛关节软骨上的软骨粘弹性。但是双层结构中粘弹性波传播规律非常复杂，波速与软骨厚度相关，测量和反演计算软骨粘弹性的难度均较大。
技术实现思路
本专利技术的目的在于提供一种应用于含硬质基底粘弹性薄板的粘弹性测量方法，其降低了对测量频率的要求，又降低了对测量数据后处理的要求为达到上述目的，本专利技术提供如下技术方案：一种应用于含硬质基底粘弹性薄板的粘弹性测量方法，所述方法包括：S1、利用不同频率的外部激励在待测材料中激励出导波；S2、利用预设成像方法记录导波的传播信息，计算近距范围及获得平均测量距离；S3、根据近距范围内导波的传播信息，计算出近距区域内波速与频率的关系，以给出频散曲线；S4、根据频散曲线及平均测量距离确定临界频率；S5、利用频散曲线，并根据表面波模型和临界频率，拟合计算出待测材料的粘弹性。进一步地，所述待测材料包括关节软骨或具有硬质基底的粘弹性薄板。进一步地，所述外部激励为机械激励或超声辐射激励。进一步地，所述预设成像包括超声成像、核磁成像和光学成像中的一种。进一步地，所述表面波模型为：其中，其中,ρ为密度，ω为角频率，υ为泊松比。进一步地，在所述步骤S4中，所利用的临界频率fcritical可表示为：进一步地，所述近距区域为靠近外部激励作用在待测材料上的驱动点的区域。进一步地，所述步骤S2还包括：测量薄板厚度h，所述近距的取值范围技术方法如下：xcritical≈1.35*h，其中，h为薄板厚度。本专利技术还提供了一种粘弹性测量系统，所述粘弹性测量系统采用上述应用于含硬质基底粘弹性薄板的粘弹性测量方法。本专利技术的有益效果在于：本专利技术的应用于含硬质基底粘弹性薄板的粘弹性测量方法利用近距区域内的临界频率，根据表面波模型，以获得粘弹性薄板的粘弹性。相比于兰姆波法，降低了对测量频率的要求；又因为表面波模型的简单，降低了对测量数据后处理的要求。上述说明仅是本专利技术技术方案的概述，为了能够更清楚了解本专利技术的技术手段，并可依照说明书的内容予以实施，以下以本专利技术的较佳实施例并配合附图详细说明如后。附图说明图1为应用于含硬质基底粘弹性薄板的粘弹性测量方法的流程图；图2为关节软骨的模型示意图；图3为关节软骨在外部激励后测试到的无量纲速度云图；图4a为近距的无量纲速度云图；图4b为根据速度云图(图4a)计算的二维FFT结果；图5a为波速与频率的关系；图5b为根据5(a)计算得到的频率与波长之间的关系。具体实施方式下面结合附图和实施例，对本专利技术的具体实施方式作进一步详细描述。以下实施例用于说明本专利技术，但不用来限制本专利技术的范围。请参见图1，本专利技术的应用于含硬质基底粘弹性薄板的粘弹性测量方法包括：S1、利用不同频率的外部激励在待测材料中激励出导波，该待测材料包括关节软骨或具有硬质基底的粘弹性薄板，所述外部激励为机械激励或超声辐射激励。S2、利用预设成像方法记录导波的传播信息，并测量薄板厚度h，计算近距范围及获得平均测量距离，其中，预设成像包括超声成像、核磁成像和光学成像中的一种；S3、根据近距范围内导波的传播信息，计算出近距区域内波速与频率的关系，以给出频散曲线；S4、根据频散曲线及平均测量距离确定临界频率；S5、利用频散曲线，并根据表面波模型和临界频率，拟合计算出待测材料的粘弹性。在步骤S5中，所述表面波模型为：其中ρ为密度，ω为角频率，υ为泊松比。此表面波模型基于Voigt材料模型。该表面波模型也可基于其它粘弹性材料模型建立，如基于Maxwell材料模型、Zener材料模型等。也可预先不假设材料模型，而根据有限窗口2DFFT方法计算出材料的粘弹性规律。所述“根据频散曲线及平均测量距离确定临界频率”中所利用的临界频率fcritical可表示为：在所述步骤S2中，所述近距区域为靠近外部激励作用在待测材料上的驱动点的区域，具体的，所述近距的取值范围技术方法如下：xcritical≈1.35*h，其中，h为薄板厚度。上述应用于含硬质基底粘弹性薄板的粘弹性测量方法可应用于电脑或计算载体中以形成粘弹性测量系统。上述应用于含硬质基底粘弹性薄板的粘弹性测量方法利用近距区域内的频散曲线，根据表面波模型获得粘弹性薄板的粘弹性。相比于兰姆波法，降低了对测量频率的要求；又因为表面波模型的简单，降低了对测量数据后处理的要求。为了更清楚的说明本专利技术的方案的优势，下面通过一详细实施方案进行比对说明。采用商用有限元软件Abaqus对关节软骨中的波传播过程进行数值模拟，需要说明的是：从力学视角看，关节软骨可简化为双层板(如图2)，其中软骨可视为粘弹性薄板1，骨可视为弹性硬质基底2，软骨和骨均采用轴对称单元(cax4r)。其中，软骨的厚度为2.5mm，密度为1100kg/m3，泊松比为0.25，粘弹性本构关系为Voigt模型(μ＝μ1+iωμ2)，剪切模量μ1为1MPa，粘性系数μ2为30Pa.s。骨的厚度为150mm，密度为743kg/m3，泊松比为0.25，本构关系为线弹性模型，杨氏模量为100MPa。在关节软骨上激发不同频率的正弦波(6000-10000Hz本文档来自技高网...

【技术保护点】
1.一种应用于含硬质基底粘弹性薄板的粘弹性测量方法，其特征在于，所述方法包括：/nS1、利用不同频率的外部激励在待测材料中激励出导波；/nS2、利用预设成像方法记录导波的传播信息，计算近距范围及获得平均测量距离；/nS3、根据近距范围内导波的传播信息，计算出近距区域内波速与频率的关系，以给出频散曲线；/nS4、根据频散曲线及平均测量距离确定临界频率；/nS5、利用频散曲线，并根据表面波模型和临界频率，拟合计算出待测材料的粘弹性。/n

【技术特征摘要】
1.一种应用于含硬质基底粘弹性薄板的粘弹性测量方法，其特征在于，所述方法包括：
S1、利用不同频率的外部激励在待测材料中激励出导波；
S2、利用预设成像方法记录导波的传播信息，计算近距范围及获得平均测量距离；
S3、根据近距范围内导波的传播信息，计算出近距区域内波速与频率的关系，以给出频散曲线；
S4、根据频散曲线及平均测量距离确定临界频率；
S5、利用频散曲线，并根据表面波模型和临界频率，拟合计算出待测材料的粘弹性。


2.如权利要求1所述的应用于含硬质基底粘弹性薄板的粘弹性测量方法，其特征在于，所述待测材料包括关节软骨或具有硬质基底的粘弹性薄板。


3.如权利要求1所述的应用于含硬质基底粘弹性薄板的粘弹性测量方法，其特征在于，所述外部激励为机械激励或超声辐射激励。


4.如权利要求1所述的应用于含硬质基底粘弹性薄板的粘弹性测量方法，其特征在于，所述预设成像包括超声成像、核磁成像和光学成像中的一种。

【专利技术属性】
技术研发人员：徐浩，罗宗平，
申请(专利权)人：苏州大学，
类型：发明
国别省市：江苏;32