[go: up one dir, main page]

CN114608995B - 一种血栓弹力测量装置的检测方法 - Google Patents

一种血栓弹力测量装置的检测方法 Download PDF

Info

Publication number
CN114608995B
CN114608995B CN202210507792.6A CN202210507792A CN114608995B CN 114608995 B CN114608995 B CN 114608995B CN 202210507792 A CN202210507792 A CN 202210507792A CN 114608995 B CN114608995 B CN 114608995B
Authority
CN
China
Prior art keywords
thromboelastometry
detection method
frequency
signal
damping ratio
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210507792.6A
Other languages
English (en)
Other versions
CN114608995A (zh
Inventor
王钟周
孙海旋
钱俊
李俊坡
吴升鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzhou Institute of Biomedical Engineering and Technology of CAS
Original Assignee
Suzhou Institute of Biomedical Engineering and Technology of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Suzhou Institute of Biomedical Engineering and Technology of CAS filed Critical Suzhou Institute of Biomedical Engineering and Technology of CAS
Priority to CN202210507792.6A priority Critical patent/CN114608995B/zh
Publication of CN114608995A publication Critical patent/CN114608995A/zh
Application granted granted Critical
Publication of CN114608995B publication Critical patent/CN114608995B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
    • G01N11/10—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties by moving a body within the material
    • G01N11/16—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties by moving a body within the material by measuring damping effect upon oscillatory body
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
    • G01N2011/0006—Calibrating, controlling or cleaning viscometers

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

本发明公开了一种血栓弹力测量装置的检测方法，属于血液分析技术领域，通过对探针端施加冲击载荷，血栓弹力测量装置在冲击载荷下作自由衰减振动，采集原始信号，对原始信号预处理得到振动信号，通过对振动信号分析得到血栓弹力测量装置固有频率和阻尼比的检测值，将固有频率和阻尼比检测值与其理论值进行对比，根据变动量与预设容差值的比较判断血栓弹力测量装置状态是否需要修正，当需要修正时，通过修正系数在线修正刚度变动导致的血栓弹力测量误差，整个检测过程不需要外部设备，可以实现在线监测血栓弹力测量装置状态以及修正测量结果。

Description

一种血栓弹力测量装置的检测方法
技术领域
本发明涉及血液分析技术领域，尤其是涉及一种血栓弹力测量装置的检测方法。
背景技术
血液凝固过程中纤维蛋白、血小板和血细胞形成三维交联的网状结构，在纤溶酶的作用下，纤维蛋白溶解。在这个过程中血液黏弹性会发生改变，通过检测凝血过程中血液黏弹性的变化来定性或定量分析凝血过程，可以帮助医生了解患者的凝血功能信息，做出准确的诊断和治疗。
血栓弹力测量装置由于弹性元件出现弯曲、装配应力松弛等问题会导致测量误差，现有技术如专利CN201710481797.5血栓弹力图仪悬垂丝劲度系数的检测方法及系统，仅对检测通道弹性元件刚度进行检测，需要外部装置，无法实现仪器在线刚度检测和血栓弹力测量结果自动修正。此外，刚度仅仅反映弹性元件状态，无法体现装配应力松弛引起的阻尼变化等隐性问题，无法全面评估检测通道状态变化。
发明内容
为了克服现有技术的不足，本发明的目的之一在于提供一种血栓弹力测量装置的检测方法，能够在线检测并且全面评估血栓弹力测量装置整体的状态变化，并自动修正血栓弹力测量结果。
本发明的目的之一采用如下技术方案实现：
一种血栓弹力测量装置的检测方法，包括以下步骤：
获取血栓弹力测量装置的固有频率和阻尼比的理论值；
对探针端施加冲击载荷，血栓弹力测量装置在冲击响应下作自由衰减振动，采集 角位移
Figure 254616DEST_PATH_IMAGE001
在
Figure 912999DEST_PATH_IMAGE002
时刻的原始信号
Figure 295439DEST_PATH_IMAGE003
，其中：
Figure 505709DEST_PATH_IMAGE004
，
Figure 606389DEST_PATH_IMAGE005
为采样间隔时间，
Figure 819196DEST_PATH_IMAGE006
；
对原始信号
Figure 311269DEST_PATH_IMAGE003
进行预处理得到振动信号
Figure 630124DEST_PATH_IMAGE007
，将基线偏离所造成的波形畸变 修正，并消除混杂于信号中的高频噪声的干扰和影响，将振动信号
Figure 31149DEST_PATH_IMAGE007
作傅里叶变换获得 复数信号
Figure 765756DEST_PATH_IMAGE008
；
将复数信号
Figure 388367DEST_PATH_IMAGE008
实部为横坐标，虚部为纵坐标，绘出随频率的奈奎斯特曲线图；
冲击载荷作用下血栓弹力测量装置的运动微分方程
Figure 550227DEST_PATH_IMAGE009
，
Figure 891079DEST_PATH_IMAGE010
和
Figure 361200DEST_PATH_IMAGE011
分别为血栓弹力 测量装置和血液的刚度，
Figure 838318DEST_PATH_IMAGE012
和
Figure 639920DEST_PATH_IMAGE013
分别为血栓弹力测量装置和血液的黏性阻尼系数，
Figure 202489DEST_PATH_IMAGE014
为血 栓弹力测量装置的转动惯量；
Figure 747740DEST_PATH_IMAGE001
为血栓弹力测量装置角位移，
Figure 141681DEST_PATH_IMAGE015
为单位冲击函数，根据运 动微分方程计算血栓弹力测量装置的频响函数
Figure 114185DEST_PATH_IMAGE016
，式中：
Figure 632891DEST_PATH_IMAGE017
，
Figure 722113DEST_PATH_IMAGE018
，
Figure 704981DEST_PATH_IMAGE019
；采用该频响函数方程对复数信号
Figure 113966DEST_PATH_IMAGE008
的奈奎 斯特曲线进行最小二乘拟合，得到血栓弹力测量装置固有频率和阻尼比的检测值；
将固有频率和阻尼比的理论值与其检测值对比，根据变动量与预设容差值的比较判断血栓弹力测量装置状态是否需要修正。
进一步地，所述频响函数
Figure 447864DEST_PATH_IMAGE020
的实部为
Figure 69338DEST_PATH_IMAGE021
；虚部为
Figure 923025DEST_PATH_IMAGE022
；式中，
Figure 237331DEST_PATH_IMAGE023
。
进一步地，
Figure 58526DEST_PATH_IMAGE024
附近，
Figure 483691DEST_PATH_IMAGE025
，
Figure 661992DEST_PATH_IMAGE026
，因此
Figure 943937DEST_PATH_IMAGE027
；
Figure 455690DEST_PATH_IMAGE029
，此时，
Figure 481284DEST_PATH_IMAGE030
和
Figure 965355DEST_PATH_IMAGE031
满足
Figure 683781DEST_PATH_IMAGE032
，该方程为圆方程，圆心坐标 为
Figure 558196DEST_PATH_IMAGE033
，半径为
Figure 856322DEST_PATH_IMAGE034
。
进一步地，采用奈奎斯特圆拟合复数信号
Figure 997497DEST_PATH_IMAGE008
，通过最小二乘法求解min{
Figure 90087DEST_PATH_IMAGE035
；
Figure 763382DEST_PATH_IMAGE036
，
Figure 334041DEST_PATH_IMAGE037
；
Figure 448496DEST_PATH_IMAGE038
；
Figure 508725DEST_PATH_IMAGE039
为 采样频率，
Figure 945348DEST_PATH_IMAGE040
，
Figure 601589DEST_PATH_IMAGE005
为采样间隔时间，可以依此来确定系统固有频率和阻尼比的检 测值。
进一步地，采用折线连接复数信号
Figure 711496DEST_PATH_IMAGE008
，该曲线与虚轴交点处
Figure 145888DEST_PATH_IMAGE024
，实部为
Figure 341246DEST_PATH_IMAGE041
，虚部
Figure 316025DEST_PATH_IMAGE042
，可以依此来确定系统固有频率和阻尼比的检测 值。
进一步地，将振动信号
Figure 14859DEST_PATH_IMAGE007
作傅里叶变换获得
Figure 416891DEST_PATH_IMAGE008
具体为：
Figure 568386DEST_PATH_IMAGE043
；
Figure 555977DEST_PATH_IMAGE037
；
Figure 374897DEST_PATH_IMAGE038
；
Figure 947830DEST_PATH_IMAGE039
为采样频率，
Figure 117780DEST_PATH_IMAGE040
，
Figure 840885DEST_PATH_IMAGE005
为采样间隔时间。
进一步地，对原始信号
Figure 311050DEST_PATH_IMAGE003
进行预处理具体为：从原始信号
Figure 258146DEST_PATH_IMAGE003
中减去最小二 乘方式最佳拟合线，通过从原始信号
Figure 197283DEST_PATH_IMAGE003
中去除线性趋势，能够将分析集中在振动信号
Figure 989659DEST_PATH_IMAGE007
的波动上。
进一步地，最小二乘方式最佳拟合线方程为
Figure 597487DEST_PATH_IMAGE044
，其中，
Figure 981064DEST_PATH_IMAGE045
，
Figure 328869DEST_PATH_IMAGE046
，
Figure 252831DEST_PATH_IMAGE047
，
Figure 635271DEST_PATH_IMAGE048
，
Figure 720908DEST_PATH_IMAGE049
=0,1,⋯N-1。
进一步地，当测杯内无血液时，
Figure 556008DEST_PATH_IMAGE050
，
Figure 221345DEST_PATH_IMAGE051
，所述血栓弹力测量装置的刚 度
Figure 464151DEST_PATH_IMAGE052
。
进一步地，所述血栓弹力测量装置的检测方法还包括修正步骤，修正步骤具体为：通过测量频率变动量，计算修正系数，在线修正刚度变动导致的测量误差。
进一步地，修正系数
Figure 986268DEST_PATH_IMAGE053
，变动量
Figure 308665DEST_PATH_IMAGE054
为血栓弹力测量装置的 刚度检测值与理论值之差，
Figure 325162DEST_PATH_IMAGE055
，转动惯量
Figure 682194DEST_PATH_IMAGE014
通过三维绘图软件计算获得。
相比现有技术，本发明一种血栓弹力测量装置的检测方法不需要外部设备，通过对探针端施加冲击载荷，血栓弹力测量装置在冲击响应下作自由衰减振动，采集原始信号，对原始信号预处理得到振动信号，通过对振动信号分析得到血栓弹力测量装置检测固有频率和阻尼比，将检测固有频率和阻尼比与其理论值对比，根据变动量与预设容差值的比较判断血栓弹力测量装置状态是否需要修正，当需要修正时，通过修正系数在线修正刚度变动导致的测量误差。
附图说明
图1为本发明血栓弹力测量装置的检测方法的检测示意图；
图2为振动信号的频响函数的奈奎斯特圆。
图中：1、测杯；2、探针。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
需要说明的是，当组件被称为“固定于”另一个组件，它可以直接在另一个组件上或者也可以存在另一中间组件，通过中间组件固定。当一个组件被认为是“连接”另一个组件，它可以是直接连接到另一个组件或者可能同时存在另一中间组件。当一个组件被认为是“设置于”另一个组件，它可以是直接设置在另一个组件上或者可能同时存在另一中间组件。本文所使用的术语“垂直的”、“水平的”、“左”、“右”以及类似的表述只是为了说明的目的。
除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本发明。本文所使用的术语“及／或”包括一个或多个相关的所列项目的任意的和所有的组合。
请参阅图1，血栓弹力测量装置用于检测血液凝固过程中血液黏弹性变化。血栓弹力测量装置包括柔性轴承、探针2以及测杯1，测杯1用于盛放血液。柔性轴承的限位轴安装探针2。检测时，杯盖装载在探针2上，通过测杯1与探针2之间的旋转运动测量血液凝固过程中血液黏弹性变化。
本发明一种血栓弹力测量装置的检测方法，检测原理为：通过对探针端施加冲击载荷，测量自由衰减振动时血栓弹力测量装置角位移，根据其频域信息绘制奈奎斯特曲线，求解固有频率和阻尼比。根据固有频率可计算血栓弹力测量装置刚度。当血栓弹力测量装置出现疲劳损伤、紧固失效或杯盖与杯身不同轴时，其阻尼比和固有频率会偏离容许值。
一种血栓弹力测量装置的检测方法，包括以下步骤：
获取血栓弹力测量装置的固有频率和阻尼比的理论值；
对探针端施加冲击载荷，血栓弹力测量装置在冲击响应下作自由衰减振动，采集 角位移
Figure 109633DEST_PATH_IMAGE001
在
Figure 450485DEST_PATH_IMAGE002
时刻的原始信号
Figure 192045DEST_PATH_IMAGE007
，其中：
Figure 669162DEST_PATH_IMAGE004
，
Figure 464906DEST_PATH_IMAGE005
为采样间隔时间，
Figure 824212DEST_PATH_IMAGE006
；
对原始信号
Figure 166200DEST_PATH_IMAGE003
进行预处理得到振动信号
Figure 497825DEST_PATH_IMAGE007
，将基线偏离所造成的波形畸变 修正，并消除混杂于信号中的高频噪声的干扰和影响，将振动信号
Figure 470329DEST_PATH_IMAGE007
作傅里叶变换获得 复数信号
Figure 802084DEST_PATH_IMAGE008
；
将复数信号
Figure 151026DEST_PATH_IMAGE008
实部为横坐标，虚部为纵坐标，绘出随频率的奈奎斯特曲线图；
冲击载荷作用下血栓弹力测量装置的运动微分方程
Figure 399473DEST_PATH_IMAGE056
，
Figure 542879DEST_PATH_IMAGE010
和
Figure 617057DEST_PATH_IMAGE011
分别为血栓弹力 测量装置和血液的刚度，
Figure 504111DEST_PATH_IMAGE012
和
Figure 607065DEST_PATH_IMAGE013
分别为血栓弹力测量装置和血液的黏性阻尼系数，
Figure 468841DEST_PATH_IMAGE014
为血 栓弹力测量装置的转动惯量；
Figure 227719DEST_PATH_IMAGE001
为血栓弹力测量装置角位移，
Figure 918463DEST_PATH_IMAGE015
为单位冲击函数，根据运 动微分方程计算血栓弹力测量装置的频响函数
Figure 79186DEST_PATH_IMAGE016
，式中：
Figure 829973DEST_PATH_IMAGE017
，
Figure 872885DEST_PATH_IMAGE018
，
Figure 588160DEST_PATH_IMAGE019
；采用该频响函数方程对复数信号
Figure 665707DEST_PATH_IMAGE008
的奈奎 斯特曲线进行最小二乘拟合，得到血栓弹力测量装置固有频率和阻尼比的检测值；
将固有频率和阻尼比的理论值与其检测值对比，根据变动量与预设容差值的比较判断血栓弹力测量装置状态是否需要修正。
具体过程如下：
获取血栓弹力测量装置的固有频率和阻尼比的理论值，此时血栓弹力测量装置的固有频率和阻尼比为原始数据，即未出现弹性元件弯曲、装配应力松弛等问题，此时血液黏弹性检测结果准确。
对探针端施加冲击载荷，血栓弹力测量装置在冲击响应下作自由衰减振动，采集 原始信号
Figure 321816DEST_PATH_IMAGE003
。
对原始信号
Figure 586444DEST_PATH_IMAGE003
进行预处理，将基线偏离所造成的波形畸变加以修正，通过带通 滤波器消除混杂于信号中的高频噪声的干扰和影响。具体的，对原始信号
Figure 618991DEST_PATH_IMAGE003
进行预处理 为：消除趋势项，线性趋势表示数据的系统性增加或减少，通常是测量环节引入的系统误 差。从原始信号
Figure 285464DEST_PATH_IMAGE003
中减去最小二乘方式最佳拟合线，通过从原始信号
Figure 643633DEST_PATH_IMAGE003
中去除线性 趋势，能够将分析集中在振动信号
Figure 864399DEST_PATH_IMAGE007
的波动上。最小二乘方式最佳拟合线方程为
Figure 982528DEST_PATH_IMAGE044
（1）
其中，
Figure 243788DEST_PATH_IMAGE057
（2），
Figure 241700DEST_PATH_IMAGE046
（3）
对振动信号
Figure 684183DEST_PATH_IMAGE007
作傅里叶变换获得复数信号
Figure 120849DEST_PATH_IMAGE008
，
Figure 965177DEST_PATH_IMAGE058
（4）
Figure 930728DEST_PATH_IMAGE037
（5）
其中
Figure 657245DEST_PATH_IMAGE038
；
血栓弹力测量装置的柔性轴承在冲击响应下作自由衰减振动时，其运动微分方程为
Figure 835285DEST_PATH_IMAGE056
（6）
Figure 81590DEST_PATH_IMAGE010
和
Figure 415445DEST_PATH_IMAGE011
分别为血栓弹力测量装置和血液的刚度，
Figure 98099DEST_PATH_IMAGE012
和
Figure 610989DEST_PATH_IMAGE013
分别为血栓弹力测量装 置和血液的黏性阻尼系数，
Figure 898750DEST_PATH_IMAGE014
为测量元件的转动惯量；
Figure 471683DEST_PATH_IMAGE001
为血栓弹力测量装置角位移，
Figure 110475DEST_PATH_IMAGE059
为单位冲击函数。
其频响函数为
Figure 427056DEST_PATH_IMAGE016
（7）
实部为
Figure 366062DEST_PATH_IMAGE060
（8）
虚部为
Figure 53438DEST_PATH_IMAGE061
（9）
式中：
Figure 241843DEST_PATH_IMAGE062
（10），
Figure 299798DEST_PATH_IMAGE018
（11），空载时
Figure 93310DEST_PATH_IMAGE051
（12），
Figure 493199DEST_PATH_IMAGE063
（13）
值得注意的是，检测固有频率和阻尼比为系统固有参数，与冲击载荷幅值无关。
Figure 168899DEST_PATH_IMAGE024
附近，
Figure 30545DEST_PATH_IMAGE025
，
Figure 209722DEST_PATH_IMAGE026
，因此
Figure 764201DEST_PATH_IMAGE064
（14）
Figure 882459DEST_PATH_IMAGE065
（15）
此时，
Figure 344533DEST_PATH_IMAGE066
和
Figure 847058DEST_PATH_IMAGE031
满足
Figure 572438DEST_PATH_IMAGE067
（16）
该方程圆心坐标为
Figure 160414DEST_PATH_IMAGE033
，半径为
Figure 832704DEST_PATH_IMAGE034
。近似后的圆方程简化了方程形式，在虚轴 交点处仍有较高精度，因此可用于参数识别。
以频响函数的实部为横坐标，虚部为纵坐标，绘出频响函数矢量随频率的Nyquist（奈奎斯特）圆，如图2所示。
采用奈奎斯特圆拟合复数信号
Figure 986474DEST_PATH_IMAGE008
，通过最小二乘法求解min{
Figure 882754DEST_PATH_IMAGE035
；
Figure 963886DEST_PATH_IMAGE036
；
Figure 705446DEST_PATH_IMAGE037
；
Figure 730034DEST_PATH_IMAGE038
；
Figure 531636DEST_PATH_IMAGE039
为 采样频率，
Figure 890942DEST_PATH_IMAGE040
，
Figure 436193DEST_PATH_IMAGE005
为采样间隔时间，可以依此来确定系统固有频率和阻尼比的检 测值。
还可以采用折线连接复数信号
Figure 236659DEST_PATH_IMAGE008
，该曲线与虚轴交点处
Figure 271480DEST_PATH_IMAGE024
，实部为
Figure 55765DEST_PATH_IMAGE068
，虚部
Figure 195585DEST_PATH_IMAGE069
，可以依此来确定系统固有频率和阻尼比的检测 值。
空载测量时，测杯内无血液，
Figure 381716DEST_PATH_IMAGE050
，
Figure 56280DEST_PATH_IMAGE070
，可计算血栓弹力测量装置刚 度，
Figure 593440DEST_PATH_IMAGE071
。
将频率和阻尼比的检测与其理论值对比，根据变动量与预设容差值的比较判断血栓弹力测量装置状态是否需要修正。
检测过程中，
Figure 27964DEST_PATH_IMAGE072
为测杯1角位移，
Figure 334180DEST_PATH_IMAGE001
为血栓弹力测量装置角位移，用于表征血栓弹 力。
Figure 914066DEST_PATH_IMAGE073
（17）
当血栓弹力测量装置刚度变动量为
Figure 469681DEST_PATH_IMAGE054
时，由此导致的测量误差
Figure 25340DEST_PATH_IMAGE074
可表示为
Figure 717221DEST_PATH_IMAGE075
（18）
弹性元件的空载固有频率变动量为
Figure 468008DEST_PATH_IMAGE076
（19）
其中，
Figure 714182DEST_PATH_IMAGE051
（20）
当血栓弹力测量装置刚度变动量为
Figure 943038DEST_PATH_IMAGE054
时，血栓弹力测量装置角位移的测量结果
Figure 755005DEST_PATH_IMAGE077
可表示为
Figure 473431DEST_PATH_IMAGE078
（21）
此时，需要将测量结果
Figure 206901DEST_PATH_IMAGE077
进行修正。设修正系数
Figure 521338DEST_PATH_IMAGE079
，用于消除因刚度变动引起的测 量误差
Figure 674232DEST_PATH_IMAGE080
（22）
解得
Figure 501242DEST_PATH_IMAGE081
（23）
可通过测量固有频率变动量，在线修正刚度变动导致的测量误差。
Figure 518746DEST_PATH_IMAGE082
时，
Figure 89404DEST_PATH_IMAGE083
，此时无需修正。
相比现有技术，本发明一种血栓弹力测量装置的检测方法不需要外部设备，通过对探针端施加冲击载荷，血栓弹力测量装置在冲击响应下作自由衰减振动，采集原始信号，通过对原始信号分析得到血栓弹力测量装置检测固有频率和阻尼比，将检测固有频率和阻尼比与固有频率和阻尼比对比，根据变动量与预设容差值的比较判断血栓弹力测量装置状态是否需要修正，当需要修正时，通过修正系数在线修正刚度变动导致的测量误差。
以上实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进演变，都是依据本发明实质技术对以上实施例做的等同修饰与演变，这些都属于本发明的保护范围。

Claims (11)

1.一种血栓弹力测量装置的检测方法，其特征在于，包括以下步骤：
获取血栓弹力测量装置的固有频率和阻尼比的理论值；
对探针端施加冲击载荷，血栓弹力测量装置在冲击响应下作自由衰减振动，采集角位移
Figure DEST_PATH_IMAGE001
在
Figure 355031DEST_PATH_IMAGE002
时刻的原始信号
Figure DEST_PATH_IMAGE003
，其中：
Figure 65105DEST_PATH_IMAGE004
，
Figure DEST_PATH_IMAGE005
为采样间隔时间，
Figure 866708DEST_PATH_IMAGE006
；
对原始信号
Figure 101380DEST_PATH_IMAGE003
进行预处理得到振动信号
Figure DEST_PATH_IMAGE007
，将基线偏离所造成的波形畸变修正，并消除混杂于信号中的高频噪声的干扰和影响，将振动信号
Figure 882517DEST_PATH_IMAGE007
作傅里叶变换获得复数信号
Figure 276458DEST_PATH_IMAGE008
；
将复数信号
Figure 327590DEST_PATH_IMAGE008
实部为横坐标，虚部为纵坐标，绘出随频率变化的奈奎斯特曲线图；
冲击载荷作用下血栓弹力测量装置的运动微分方程
Figure 678587DEST_PATH_IMAGE010
，
Figure DEST_PATH_IMAGE011
和
Figure 27529DEST_PATH_IMAGE012
分别为血栓弹力测量装置和血液的刚度，
Figure DEST_PATH_IMAGE013
和
Figure 823446DEST_PATH_IMAGE014
分别为血栓弹力测量装置和血液的黏性阻尼系数，
Figure 530634DEST_PATH_IMAGE016
为血栓弹力测量装置的转动惯量；
Figure 880843DEST_PATH_IMAGE001
为血栓弹力测量装置角位移，
Figure DEST_PATH_IMAGE017
为单位冲击函数，根据运动微分方程计算血栓弹力测量装置的频响函数
Figure 502318DEST_PATH_IMAGE018
，式中：
Figure DEST_PATH_IMAGE019
，
Figure 900545DEST_PATH_IMAGE020
，
Figure DEST_PATH_IMAGE021
；采用该频响函数的方程对复数信号
Figure 949272DEST_PATH_IMAGE008
的奈奎斯特曲线进行最小二乘拟合，得到血栓弹力测量装置固有频率和阻尼比的检测值；
将固有频率和阻尼比的理论值与其检测值对比，根据变动量与预设容差值的比较判断血栓弹力测量装置状态是否需要修正。
2.根据权利要求1所述的血栓弹力测量装置的检测方法，其特征在于：所述频响函数
Figure 521199DEST_PATH_IMAGE022
的实部为
Figure DEST_PATH_IMAGE023
；虚部为
Figure 306883DEST_PATH_IMAGE024
；式中，
Figure DEST_PATH_IMAGE025
。
3.根据权利要求2所述的血栓弹力测量装置的检测方法，其特征在于：
Figure 202027DEST_PATH_IMAGE026
附近，
Figure DEST_PATH_IMAGE027
，
Figure 31443DEST_PATH_IMAGE028
，因此
Figure DEST_PATH_IMAGE029
；
Figure 504713DEST_PATH_IMAGE030
，此时，
Figure DEST_PATH_IMAGE031
和
Figure 592623DEST_PATH_IMAGE032
满足
Figure DEST_PATH_IMAGE033
，该方程为圆方程，圆心坐标为
Figure 578159DEST_PATH_IMAGE034
，半径为
Figure DEST_PATH_IMAGE035
。
4.根据权利要求3所述的血栓弹力测量装置的检测方法，其特征在于：采用奈奎斯特圆拟合复数信号
Figure 844056DEST_PATH_IMAGE008
，通过最小二乘法求解min{
Figure 374263DEST_PATH_IMAGE036
；
Figure DEST_PATH_IMAGE037
，
Figure 16597DEST_PATH_IMAGE038
；
Figure DEST_PATH_IMAGE039
；
Figure 119289DEST_PATH_IMAGE040
为采样频率，
Figure DEST_PATH_IMAGE041
，
Figure 946299DEST_PATH_IMAGE005
为采样间隔时间，可以依此来确定系统固有频率和阻尼比的检测值。
5.根据权利要求3所述的血栓弹力测量装置的检测方法，其特征在于：采用折线连接复数信号
Figure 839169DEST_PATH_IMAGE008
，该曲线与虚轴交点处
Figure 973609DEST_PATH_IMAGE026
，实部为
Figure 510901DEST_PATH_IMAGE042
，虚部
Figure DEST_PATH_IMAGE043
，可以依此来确定系统固有频率和阻尼比的检测值。
6.根据权利要求1所述的血栓弹力测量装置的检测方法，其特征在于：将振动信号
Figure 367867DEST_PATH_IMAGE007
作傅里叶变换获得
Figure 623399DEST_PATH_IMAGE008
具体为：
Figure 298881DEST_PATH_IMAGE044
；
Figure 956258DEST_PATH_IMAGE038
；
Figure 62755DEST_PATH_IMAGE039
；
Figure 54850DEST_PATH_IMAGE040
为采样频率，
Figure 780361DEST_PATH_IMAGE041
，
Figure 42977DEST_PATH_IMAGE005
为采样间隔时间。
7.根据权利要求1所述的血栓弹力测量装置的检测方法，其特征在于：对原始信号
Figure 320375DEST_PATH_IMAGE003
进行预处理具体为：从原始信号
Figure 550499DEST_PATH_IMAGE003
中减去最小二乘方式最佳拟合线，通过从原始信号
Figure 328968DEST_PATH_IMAGE003
数据中去除线性趋势得到振动信号
Figure 695358DEST_PATH_IMAGE007
，能够将分析集中在振动信号
Figure 766826DEST_PATH_IMAGE007
的波动上。
8.根据权利要求6所述的血栓弹力测量装置的检测方法，其特征在于：最小二乘方式最佳拟合线方程为
Figure DEST_PATH_IMAGE045
，其中，
Figure 280984DEST_PATH_IMAGE046
，
Figure DEST_PATH_IMAGE047
，
Figure 269669DEST_PATH_IMAGE048
，
Figure DEST_PATH_IMAGE049
，
Figure DEST_PATH_IMAGE051
。
9.根据权利要求1所述的血栓弹力测量装置的检测方法，其特征在于：当测杯内无血液时，
Figure 975719DEST_PATH_IMAGE052
，
Figure DEST_PATH_IMAGE053
，所述血栓弹力测量装置的刚度
Figure 922815DEST_PATH_IMAGE054
，转动惯量
Figure 627333DEST_PATH_IMAGE016
通过三维绘图软件计算获得。
10.根据权利要求1所述的血栓弹力测量装置的检测方法，其特征在于：所述血栓弹力测量装置的检测方法还包括修正步骤，修正步骤具体为：通过测量频率变动量，计算修正系数，在线修正刚度变动导致的测量误差。
11.根据权利要求9所述的血栓弹力测量装置的检测方法，其特征在于：修正系数
Figure DEST_PATH_IMAGE055
，变动量
Figure 560654DEST_PATH_IMAGE056
为血栓弹力测量装置的刚度检测值与理论值之差，
Figure DEST_PATH_IMAGE057
，转动惯量
Figure 823008DEST_PATH_IMAGE016
通过三维绘图软件计算获得。
CN202210507792.6A 2022-05-11 2022-05-11 一种血栓弹力测量装置的检测方法 Active CN114608995B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210507792.6A CN114608995B (zh) 2022-05-11 2022-05-11 一种血栓弹力测量装置的检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210507792.6A CN114608995B (zh) 2022-05-11 2022-05-11 一种血栓弹力测量装置的检测方法

Publications (2)

Publication Number Publication Date
CN114608995A CN114608995A (zh) 2022-06-10
CN114608995B true CN114608995B (zh) 2022-08-12

Family

ID=81870471

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210507792.6A Active CN114608995B (zh) 2022-05-11 2022-05-11 一种血栓弹力测量装置的检测方法

Country Status (1)

Country Link
CN (1) CN114608995B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1226901A (zh) * 1968-04-05 1971-03-31
US4341111A (en) * 1979-03-05 1982-07-27 Fresenius Ag Process and apparatus for determining the visco elastic characteristics of fluids
CN106442952A (zh) * 2016-08-30 2017-02-22 诺泰科生物科技（苏州）有限公司 一种血栓弹力仪及其校准方法
CN109520889A (zh) * 2012-01-16 2019-03-26 仪宝科技公司 用于测量流体物理性能的方法、装置、和系统
CN110547825A (zh) * 2019-09-03 2019-12-10 上海交通大学 一种基于多频剪切波幅值分析的超声弹性成像技术
CN113640377A (zh) * 2021-10-15 2021-11-12 中国科学院苏州生物医学工程技术研究所 一种基于振动传感器的血液黏弹性测量方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1226901A (zh) * 1968-04-05 1971-03-31
US4341111A (en) * 1979-03-05 1982-07-27 Fresenius Ag Process and apparatus for determining the visco elastic characteristics of fluids
CN109520889A (zh) * 2012-01-16 2019-03-26 仪宝科技公司 用于测量流体物理性能的方法、装置、和系统
CN106442952A (zh) * 2016-08-30 2017-02-22 诺泰科生物科技（苏州）有限公司 一种血栓弹力仪及其校准方法
CN110547825A (zh) * 2019-09-03 2019-12-10 上海交通大学 一种基于多频剪切波幅值分析的超声弹性成像技术
CN113640377A (zh) * 2021-10-15 2021-11-12 中国科学院苏州生物医学工程技术研究所 一种基于振动传感器的血液黏弹性测量方法

Also Published As

Publication number Publication date
CN114608995A (zh) 2022-06-10

Similar Documents

Publication Publication Date Title
US10317275B2 (en) Vibration monitoring systems
CN110567574A (zh) 一种旋转叶片叶端定时振动参数辨识方法与系统
CN110186549A (zh) 基于叶尖定时传感器的叶片振动识别方法
CN118670370B (zh) 一种基于半球谐振子的角速度测量系统及方法
CN115728063A (zh) 一种柔性轴承检测方法
CN113820003A (zh) 一种适用于桥梁振动监测的加速度实时重构动态位移方法
CN107014482B (zh) 振动状态的在线监测装置和方法
JP3313028B2 (ja) 張力のかかっているケーブルの曲げ剛性及び張力の測定方法
CN114608995B (zh) 一种血栓弹力测量装置的检测方法
JP2020144111A (ja) 異常検出装置、回転機械、異常検出方法、及びプログラム
Chen et al. Torsional vibration feature extraction method from lateral vibrations based on decomposed forward and backward whirl motions
CN117332205B (zh) 压电阻抗温度补偿高精度自动优化方法及装置
CN118836964A (zh) 噪音测试方法、声压信号测试系统、分析设备及存储介质
JP2009011540A (ja) 生体検査装置及び生体検査装置のプログラム
CN117949155B (zh) 一种柔性铰链刚度检测方法
CN111684213A (zh) 机器人故障诊断方法、系统及存储装置
CN108375415A (zh) 一种现场环境下输电线路微风振动监测装置校准方法
CN209589491U (zh) 一种车内压力波动监测机构
CN110779723B (zh) 一种基于霍尔信号的变速工况电机轴承精确故障诊断方法
CN114486252A (zh) 一种矢量模极大值包络的滚动轴承故障诊断方法
CN112082641A (zh) 一种振动传感器振动测量数据的修正处理方法
CN117969350B (zh) 一种通过旋转振荡测量血液黏弹性的方法
CN113567123B (zh) 旋转机械冲击类故障自动诊断方法
JPS6112690B2 (zh)
JPH0450731A (ja) 回転機故障診断システム

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant