CN114608995B - Detection method of thrombus elasticity measuring device - Google Patents

Detection method of thrombus elasticity measuring device 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
elasticity measuring
thrombus elasticity
measuring device
thrombus
signal
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
Chinese (zh)
Other versions
CN114608995A (en
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/en
Publication of CN114608995A publication Critical patent/CN114608995A/en
Application granted granted Critical
Publication of CN114608995B publication Critical patent/CN114608995B/en
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

The invention discloses a detection method of a thrombus elasticity measuring device, which belongs to the technical field of blood analysis, wherein an impact load is applied to a probe end, the thrombus elasticity measuring device freely damps vibration under the impact load, an original signal is collected, the original signal is preprocessed to obtain a vibration signal, a detection value of the natural frequency and the damping ratio of the thrombus elasticity measuring device is obtained by analyzing the vibration signal, the detection value of the natural frequency and the damping ratio is compared with a theoretical value, whether the state of the thrombus elasticity measuring device needs to be corrected or not is judged according to the comparison of a variation and a preset tolerance value, when the correction is needed, the thrombus elasticity measurement error caused by rigidity variation is corrected on line through a correction coefficient, no external equipment is needed in the whole detection process, and the online monitoring of the state of the thrombus elasticity measuring device and the correction of the measurement result can be realized.

Description

Detection method of thrombus elasticity measuring device
Technical Field
The invention relates to the technical field of blood analysis, in particular to a detection method of a thrombus elasticity measuring device.
Background
Fibrin, platelets and blood cells form a three-dimensional cross-linked reticular structure in the blood coagulation process, and the fibrin is dissolved under the action of plasmin. The blood viscoelasticity can be changed in the process, and the blood coagulation process can be qualitatively or quantitatively analyzed by detecting the change of the blood viscoelasticity in the blood coagulation process, so that doctors can know the blood coagulation function information of patients to make accurate diagnosis and treatment.
In the thrombus elasticity measuring device, measurement errors can be caused due to the problems of bending, assembly stress relaxation and the like of an elastic element, and in the prior art, for example, in a method and a system for detecting the stiffness coefficient of a suspension wire of a thrombus elastography instrument disclosed by patent CN201710481797.5, only the stiffness of the elastic element of a detection channel is detected, an external device is needed, and online stiffness detection and automatic correction of a thrombus elasticity measuring result of the instrument cannot be realized. In addition, the rigidity only reflects the state of the elastic element, so that the hidden problems such as damping change and the like caused by assembly stress relaxation cannot be reflected, and the state change of the detection channel cannot be comprehensively evaluated.
Disclosure of Invention
In order to overcome the defects of the prior art, an object of the present invention is to provide a method for detecting a thrombus elasticity measuring device, which can detect and comprehensively evaluate the overall state change of the thrombus elasticity measuring device on line, and automatically correct the thrombus elasticity measuring result.
One of the purposes of the invention is realized by adopting the following technical scheme:
a detection method of a thrombus elasticity measuring device comprises the following steps:
acquiring theoretical values of natural frequency and damping ratio of the thrombus elasticity measuring device;
applying impact load to the probe end, the thrombus elasticity measuring device performs free damping vibration under the impact response, and the angular displacement is collected
Figure 254616DEST_PATH_IMAGE001
In that
Figure 912999DEST_PATH_IMAGE002
Original signal of time of day
Figure 295439DEST_PATH_IMAGE003
Wherein:
Figure 505709DEST_PATH_IMAGE004
，
Figure 606389DEST_PATH_IMAGE005
in order to sample the interval time of the sampling,
Figure 819196DEST_PATH_IMAGE006
；
for original signal
Figure 311269DEST_PATH_IMAGE003
Preprocessing the signal to obtain a vibration signal
Figure 630124DEST_PATH_IMAGE007
Correcting the waveform distortion caused by the baseline deviation and eliminating the high frequency noise mixed in the signalAcoustic interference and influence, vibration signal
Figure 31149DEST_PATH_IMAGE007
Fourier transform to obtain complex signal
Figure 765756DEST_PATH_IMAGE008
；
A complex signal
Figure 388367DEST_PATH_IMAGE008
The real part is an abscissa, the imaginary part is an ordinate, and a Nyquist plot with frequency is drawn;
motion differential equation of thrombus elasticity measuring device under impact load effect
Figure 550227DEST_PATH_IMAGE009
，
Figure 891079DEST_PATH_IMAGE010
And
Figure 361200DEST_PATH_IMAGE011
the rigidity of the thrombus elasticity measuring device and the blood respectively,
Figure 838318DEST_PATH_IMAGE012
and
Figure 639920DEST_PATH_IMAGE013
respectively the viscosity damping coefficients of the thrombus elasticity measuring device and the blood,
Figure 202489DEST_PATH_IMAGE014
the rotational inertia of the thrombus elasticity measuring device;
Figure 747740DEST_PATH_IMAGE001
in order to measure the angular displacement of the thrombus elasticity measuring device,
Figure 141681DEST_PATH_IMAGE015
calculating blood according to differential equation of motion as unit shock functionFrequency response function of bolt elasticity measuring device
Figure 114185DEST_PATH_IMAGE016
In the formula:
Figure 632891DEST_PATH_IMAGE017
，
Figure 722113DEST_PATH_IMAGE018
，
Figure 704981DEST_PATH_IMAGE019
(ii) a Using the frequency response function equation to complex signals
Figure 113966DEST_PATH_IMAGE008
Performing least square fitting on the Nyquist curve to obtain detection values of the natural frequency and the damping ratio of the thrombus elasticity measuring device;
and comparing theoretical values of the natural frequency and the damping ratio with detection values of the natural frequency and the damping ratio, and judging whether the state of the thrombus elasticity measuring device needs to be corrected or not according to comparison of the variation and a preset tolerance value.
Further, the frequency response function
Figure 447864DEST_PATH_IMAGE020
Has a real part of
Figure 69338DEST_PATH_IMAGE021
(ii) a Imaginary part of
Figure 923025DEST_PATH_IMAGE022
(ii) a In the formula (I), the compound is shown in the specification,
Figure 237331DEST_PATH_IMAGE023
。
further, the air conditioner is provided with a fan,
Figure 58526DEST_PATH_IMAGE024
in the vicinity of the location of the mobile station,
Figure 483691DEST_PATH_IMAGE025
，
Figure 661992DEST_PATH_IMAGE026
thus, therefore, it is
Figure 943937DEST_PATH_IMAGE027
；
Figure 455690DEST_PATH_IMAGE029
At this time, the process of the present invention,
Figure 481284DEST_PATH_IMAGE030
and
Figure 965355DEST_PATH_IMAGE031
satisfy the requirement of
Figure 683781DEST_PATH_IMAGE032
The equation is a circular equation with the coordinates of the center of the circle
Figure 558196DEST_PATH_IMAGE033
Radius of
Figure 856322DEST_PATH_IMAGE034
。
Further, fitting a Nyquist circle to the complex signal
Figure 997497DEST_PATH_IMAGE008
Solving min back face by least square method
Figure 90087DEST_PATH_IMAGE035
；
Figure 763382DEST_PATH_IMAGE036
，
Figure 334041DEST_PATH_IMAGE037
；
Figure 448496DEST_PATH_IMAGE038
；
Figure 508725DEST_PATH_IMAGE039
In order to be able to sample the frequency,
Figure 945348DEST_PATH_IMAGE040
，
Figure 601589DEST_PATH_IMAGE005
the detection values of the natural frequency and the damping ratio of the system can be determined according to the sampling interval.
Further, a plurality of signals are connected by adopting a broken line
Figure 711496DEST_PATH_IMAGE008
At the intersection of the curve and the imaginary axis
Figure 145888DEST_PATH_IMAGE024
Real part of
Figure 341246DEST_PATH_IMAGE041
Imaginary part
Figure 316025DEST_PATH_IMAGE042
The detected values of the natural frequency and the damping ratio of the system can be determined accordingly.
Further, the vibration signal is transmitted
Figure 14859DEST_PATH_IMAGE007
Obtained by Fourier transform
Figure 416891DEST_PATH_IMAGE008
The method specifically comprises the following steps:
Figure 568386DEST_PATH_IMAGE043
；
Figure 555977DEST_PATH_IMAGE037
；
Figure 374897DEST_PATH_IMAGE038
；
Figure 947830DEST_PATH_IMAGE039
in order to be able to sample the frequency,
Figure 117780DEST_PATH_IMAGE040
，
Figure 840885DEST_PATH_IMAGE005
is the sampling interval time.
Further, for the original signal
Figure 311050DEST_PATH_IMAGE003
The pretreatment specifically comprises the following steps: from the original signal
Figure 258146DEST_PATH_IMAGE003
Subtracting the least squares best fit line from the original signal
Figure 197283DEST_PATH_IMAGE003
The linear trend is removed, and the analysis can be concentrated on the vibration signal
Figure 989659DEST_PATH_IMAGE007
On the fluctuation of (c).
Further, the least squares best fit line equation is
Figure 597487DEST_PATH_IMAGE044
Wherein, in the step (A),
Figure 981064DEST_PATH_IMAGE045
，
Figure 328869DEST_PATH_IMAGE046
，
Figure 252831DEST_PATH_IMAGE047
，
Figure 635271DEST_PATH_IMAGE048
，
Figure 720908DEST_PATH_IMAGE049
=0,1,⋯N-1。
further, when no blood is in the measuring cup,
Figure 556008DEST_PATH_IMAGE050
，
Figure 221345DEST_PATH_IMAGE051
the rigidity of the thrombus elasticity measuring apparatus
Figure 464151DEST_PATH_IMAGE052
。
Further, the detection method of the thrombus elasticity measurement device further comprises a correction step, wherein the correction step specifically comprises the following steps: and calculating a correction coefficient by measuring the frequency variation, and correcting the measurement error caused by the rigidity variation on line.
Further, the correction coefficient
Figure 986268DEST_PATH_IMAGE053
Amount of variation
Figure 308665DEST_PATH_IMAGE054
Is the difference between the rigidity detection value and the theoretical value of the thrombus elasticity measuring device,
Figure 325162DEST_PATH_IMAGE055
moment of inertia
Figure 682194DEST_PATH_IMAGE014
Calculated by three-dimensional drawing software.
Compared with the prior art, the detection method of the thrombus elasticity measuring device does not need external equipment, the thrombus elasticity measuring device freely attenuates and vibrates under the impact response by applying impact load to the probe end, the original signal is collected and preprocessed to obtain the vibration signal, the natural frequency and the damping ratio of the thrombus elasticity measuring device are obtained by analyzing the vibration signal, the detected natural frequency and the damping ratio are compared with the theoretical value, whether the state of the thrombus elasticity measuring device needs to be corrected or not is judged according to the comparison between the variation and the preset tolerance value, and when the state needs to be corrected, the measurement error caused by the variation of the rigidity is corrected on line through the correction coefficient.
Drawings
FIG. 1 is a schematic view showing the detection of the detection method of the thrombus elasticity measuring apparatus according to the present invention;
fig. 2 is a nyquist circle of a frequency response function of a vibration signal.
In the figure: 1. measuring the cup; 2. and (3) a probe.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It will be understood that when an element is referred to as being "secured to" another element, it can be directly on the other element or intervening elements may also be present, secured by intervening elements. When a component is referred to as being "connected" to another component, it can be directly connected to the other component or intervening components may also be present. When an element is referred to as being "disposed on" another element, it can be directly disposed on the other element or intervening elements may also be present. The terms "vertical," "horizontal," "left," "right," and the like as used herein are for illustrative purposes only.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
Referring to fig. 1, the thrombus elastic force measuring device is used for detecting the change of blood viscoelasticity during blood coagulation. Thrombus elasticity measuring device includes flexible bearing, probe 2 and surveys cup 1, surveys cup 1 and is used for holding blood. The probe 2 is arranged on the limiting shaft of the flexible bearing. During detection, the cup cover is loaded on the probe 2, and the change of the viscoelasticity of blood in the blood coagulation process is measured through the rotary motion between the measuring cup 1 and the probe 2.
The invention relates to a detection method of a thrombus elasticity measuring device, which has the following detection principle: the angular displacement of the thrombus elasticity measuring device during free damping vibration is measured by applying impact load to the probe end, a Nyquist curve is drawn according to frequency domain information of the angular displacement measuring device, and the natural frequency and the damping ratio are solved. The rigidity of the thrombus elasticity measuring device can be calculated according to the natural frequency. When the thrombus elasticity measuring device is subjected to fatigue damage, fastening failure or the cup cover and the cup body are not coaxial, the damping ratio and the natural frequency of the thrombus elasticity measuring device deviate from allowable values.
A detection method of a thrombus elasticity measuring device comprises the following steps:
acquiring the theoretical values of the natural frequency and the damping ratio of the thrombus elasticity measuring device;
applying impact load to the probe end, the thrombus elasticity measuring device performs free damping vibration under the impact response, and the angular displacement is collected
Figure 109633DEST_PATH_IMAGE001
In that
Figure 450485DEST_PATH_IMAGE002
Original signal of time of day
Figure 192045DEST_PATH_IMAGE007
Wherein:
Figure 669162DEST_PATH_IMAGE004
，
Figure 464906DEST_PATH_IMAGE005
in order to sample the interval time of the sampling,
Figure 824212DEST_PATH_IMAGE006
；
for original signal
Figure 166200DEST_PATH_IMAGE003
Preprocessing the signal to obtain a vibration signal
Figure 497825DEST_PATH_IMAGE007
Correcting the waveform distortion caused by the baseline deviation, eliminating the interference and influence of high-frequency noise mixed in the signal, and vibrating the signal
Figure 470329DEST_PATH_IMAGE007
Fourier transform to obtain complex signal
Figure 802084DEST_PATH_IMAGE008
；
A complex signal
Figure 151026DEST_PATH_IMAGE008
The real part is an abscissa, the imaginary part is an ordinate, and a Nyquist plot with frequency is drawn;
motion differential equation of thrombus elasticity measuring device under impact load effect
Figure 399473DEST_PATH_IMAGE056
，
Figure 542879DEST_PATH_IMAGE010
And
Figure 617057DEST_PATH_IMAGE011
the rigidity of the thrombus elasticity measuring device and the blood respectively,
Figure 504111DEST_PATH_IMAGE012
and
Figure 607065DEST_PATH_IMAGE013
respectively the viscosity damping coefficients of the thrombus elasticity measuring device and the blood,
Figure 468841DEST_PATH_IMAGE014
the rotational inertia of the thrombus elasticity measuring device;
Figure 227719DEST_PATH_IMAGE001
in order to measure the angular displacement of the thrombus elasticity measuring device,
Figure 918463DEST_PATH_IMAGE015
calculating the frequency response function of the thrombus elasticity measuring device according to the motion differential equation as a unit impact function
Figure 79186DEST_PATH_IMAGE016
In the formula:
Figure 829973DEST_PATH_IMAGE017
，
Figure 872885DEST_PATH_IMAGE018
，
Figure 588160DEST_PATH_IMAGE019
(ii) a Using the frequency response function equation to complex signals
Figure 665707DEST_PATH_IMAGE008
Performing least square fitting on the Nyquist curve to obtain detection values of the natural frequency and the damping ratio of the thrombus elasticity measuring device;
and comparing theoretical values of the natural frequency and the damping ratio with detection values of the natural frequency and the damping ratio, and judging whether the state of the thrombus elasticity measuring device needs to be corrected or not according to comparison of the variation and a preset tolerance value.
The specific process is as follows:
and acquiring theoretical values of the natural frequency and the damping ratio of the thrombus elasticity measuring device, wherein the natural frequency and the damping ratio of the thrombus elasticity measuring device are original data, namely, the problems of bending of an elastic element, relaxation of assembly stress and the like do not occur, and the blood viscoelasticity detection result is accurate.
Applying impact load to the probe end, the thrombus elasticity measuring device performs free damping vibration under the impact response, and original signals are collected
Figure 321816DEST_PATH_IMAGE003
。
For original signal
Figure 586444DEST_PATH_IMAGE003
Preprocessing is carried out, waveform distortion caused by baseline deviation is corrected, and interference and influence of high-frequency noise mixed in signals are eliminated through a band-pass filter. In particular, for the original signal
Figure 618991DEST_PATH_IMAGE003
The pretreatment is carried out as follows: the trend term is eliminated, and the linear trend represents the systematic increase or decrease of data, and is usually the systematic error introduced by the measuring link. From the original signal
Figure 285464DEST_PATH_IMAGE003
Subtracting the least squares best fit line from the original signal
Figure 643633DEST_PATH_IMAGE003
The linear trend is removed, and the analysis can be concentrated on the vibration signal
Figure 864399DEST_PATH_IMAGE007
On the fluctuation of (c). The least squares best fit line equation is
Figure 982528DEST_PATH_IMAGE044
（1）
Wherein the content of the first and second substances,
Figure 243788DEST_PATH_IMAGE057
（2），
Figure 241700DEST_PATH_IMAGE046
（3）
for vibration signal
Figure 684183DEST_PATH_IMAGE007
Fourier transform to obtain complex signal
Figure 120849DEST_PATH_IMAGE008
，
Figure 965177DEST_PATH_IMAGE058
（4）
Figure 930728DEST_PATH_IMAGE037
（5）
Wherein
Figure 657245DEST_PATH_IMAGE038
；
When the flexible bearing of the thrombus elasticity measuring device does free damping vibration under the impact response, the motion differential equation is
Figure 835285DEST_PATH_IMAGE056
（6）
Figure 81590DEST_PATH_IMAGE010
And
Figure 415445DEST_PATH_IMAGE011
the rigidity of the thrombus elasticity measuring device and the blood respectively,
Figure 98099DEST_PATH_IMAGE012
and
Figure 610989DEST_PATH_IMAGE013
respectively the viscosity damping coefficients of the thrombus elasticity measuring device and the blood,
Figure 898750DEST_PATH_IMAGE014
is the moment of inertia of the measuring element;
Figure 471683DEST_PATH_IMAGE001
in order to measure the angular displacement of the thrombus elasticity measuring device,
Figure 110475DEST_PATH_IMAGE059
is a unit impact function.
Having a frequency response function of
Figure 427056DEST_PATH_IMAGE016
（7）
Real part of
Figure 366062DEST_PATH_IMAGE060
（8）
Imaginary part of
Figure 53438DEST_PATH_IMAGE061
（9）
In the formula:
Figure 241843DEST_PATH_IMAGE062
（10），
Figure 299798DEST_PATH_IMAGE018
(11) at no load
Figure 93310DEST_PATH_IMAGE051
（12），
Figure 493199DEST_PATH_IMAGE063
（13）
It is noted that the natural frequency and damping ratio are detected as system-inherent parameters, independent of the impact load amplitude.
Figure 168899DEST_PATH_IMAGE024
In the vicinity of the location of the mobile station,
Figure 30545DEST_PATH_IMAGE025
，
Figure 209722DEST_PATH_IMAGE026
thus, therefore, it is
Figure 764201DEST_PATH_IMAGE064
（14）
Figure 882459DEST_PATH_IMAGE065
（15）
At this time, the process of the present invention,
Figure 344533DEST_PATH_IMAGE066
and
Figure 847058DEST_PATH_IMAGE031
satisfy the requirement of
Figure 572438DEST_PATH_IMAGE067
（16）
The center coordinates of the equation are
Figure 160414DEST_PATH_IMAGE033
Radius of
Figure 832704DEST_PATH_IMAGE034
. The approximated circular equation simplifies the equation form, and the virtual axis intersection still has higher precision, so the method can be used for parameter identification.
The real part of the frequency response function is used as an abscissa, and the imaginary part is used as an ordinate, so as to draw a Nyquist circle of the frequency response function vector with frequency, as shown in fig. 2.
Fitting complex signals using Nyquist circles
Figure 986474DEST_PATH_IMAGE008
Solving min back face by least square method
Figure 882754DEST_PATH_IMAGE035
；
Figure 963886DEST_PATH_IMAGE036
；
Figure 705446DEST_PATH_IMAGE037
；
Figure 730034DEST_PATH_IMAGE038
；
Figure 531636DEST_PATH_IMAGE039
In order to be able to sample the frequency,
Figure 890942DEST_PATH_IMAGE040
，
Figure 436193DEST_PATH_IMAGE005
the detection values of the natural frequency and the damping ratio of the system can be determined according to the sampling interval.
It can also adopt broken line to connect complex signals
Figure 236659DEST_PATH_IMAGE008
At the intersection of the curve and the imaginary axis
Figure 271480DEST_PATH_IMAGE024
Real part of
Figure 55765DEST_PATH_IMAGE068
Imaginary part
Figure 195585DEST_PATH_IMAGE069
The detected values of the natural frequency and the damping ratio of the system can be determined accordingly.
When in no-load measurement, the measuring cup has no blood,
Figure 381716DEST_PATH_IMAGE050
，
Figure 56280DEST_PATH_IMAGE070
the rigidity of the thrombus elasticity measuring device can be calculated,
Figure 593440DEST_PATH_IMAGE071
。
and comparing the detection of the frequency and the damping ratio with theoretical values of the frequency and the damping ratio, and judging whether the state of the thrombus elasticity measuring device needs to be corrected or not according to the comparison of the variation and a preset tolerance value.
In the process of the detection, the detection device is used,
Figure 27964DEST_PATH_IMAGE072
in order to measure the angular displacement of the cup 1,
Figure 334180DEST_PATH_IMAGE001
the angular displacement of the thrombus elasticity measuring device is used for representing the thrombus elasticity.
Figure 914066DEST_PATH_IMAGE073
（17）
When the rigidity of the thrombus elasticity measuring device is changed by
Figure 469681DEST_PATH_IMAGE054
Time, measurement error caused thereby
Figure 25340DEST_PATH_IMAGE074
Can be expressed as
Figure 717221DEST_PATH_IMAGE075
（18）
The no-load natural frequency variation of the elastic element is
Figure 468008DEST_PATH_IMAGE076
（19）
Wherein the content of the first and second substances,
Figure 714182DEST_PATH_IMAGE051
（20）
when the rigidity of the thrombus elasticity measuring device is changed by
Figure 943038DEST_PATH_IMAGE054
Sometimes, thrombusMeasurement result of angular displacement of elastic force measuring device
Figure 755005DEST_PATH_IMAGE077
Can be expressed as
Figure 473431DEST_PATH_IMAGE078
（21）
At this time, it is necessary to measure the results
Figure 206901DEST_PATH_IMAGE077
And (6) correcting. Setting the correction factor
Figure 521338DEST_PATH_IMAGE079
For eliminating measurement errors due to stiffness variations
Figure 674232DEST_PATH_IMAGE080
（22）
Get it solved
Figure 501242DEST_PATH_IMAGE081
（23）
The measurement error caused by the rigidity variation can be corrected on line by measuring the natural frequency variation.
Figure 518746DEST_PATH_IMAGE082
When the temperature of the water is higher than the set temperature,
Figure 89404DEST_PATH_IMAGE083
no correction is needed at this time.
Compared with the prior art, the detection method of the thrombus elasticity measuring device does not need external equipment, the thrombus elasticity measuring device does free damping vibration under the impact response by applying the impact load to the probe end, the original signal is acquired, the natural frequency and the damping ratio of the thrombus elasticity measuring device are detected by analyzing the original signal, the detected natural frequency and the damping ratio are compared with the natural frequency and the damping ratio, whether the state of the thrombus elasticity measuring device needs to be corrected or not is judged according to the comparison between the variation and the preset tolerance value, and when the correction is needed, the measurement error caused by the variation of the rigidity is corrected on line through the correction coefficient.
The above examples only show some embodiments of the present invention, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for those skilled in the art, variations and modifications can be made without departing from the spirit of the invention, and all equivalent modifications and changes can be made to the above embodiments according to the essential technology of the invention, which falls into the protection scope of the invention.

Claims (11)

1. A detection method of a thrombus elasticity measuring device is characterized by comprising the following steps:
acquiring theoretical values of natural frequency and damping ratio of the thrombus elasticity measuring device;
applying impact load to the probe end, the thrombus elasticity measuring device performs free damping vibration under the impact response, and the angular displacement is collected
Figure DEST_PATH_IMAGE001
In that
Figure 355031DEST_PATH_IMAGE002
Original signal of time of day
Figure DEST_PATH_IMAGE003
Wherein:
Figure 65105DEST_PATH_IMAGE004
，
Figure DEST_PATH_IMAGE005
in order to sample the interval time of the sampling,
Figure 866708DEST_PATH_IMAGE006
；
for original signal
Figure 101380DEST_PATH_IMAGE003
Preprocessing the signal to obtain a vibration signal
Figure DEST_PATH_IMAGE007
Correcting the waveform distortion caused by the baseline deviation, eliminating the interference and influence of high-frequency noise mixed in the signal, and vibrating the signal
Figure 882517DEST_PATH_IMAGE007
Fourier transform to obtain complex signal
Figure 276458DEST_PATH_IMAGE008
；
A complex signal
Figure 327590DEST_PATH_IMAGE008
The real part is an abscissa, the imaginary part is an ordinate, and a Nyquist curve graph which changes along with the frequency is drawn;
motion differential equation of thrombus elasticity measuring device under impact load effect
Figure 678587DEST_PATH_IMAGE010
，
Figure DEST_PATH_IMAGE011
And
Figure 27529DEST_PATH_IMAGE012
the rigidity of the thrombus elasticity measuring device and the blood respectively,
Figure DEST_PATH_IMAGE013
and
Figure 823446DEST_PATH_IMAGE014
respectively the viscosity damping coefficients of the thrombus elasticity measuring device and the blood,
Figure 530634DEST_PATH_IMAGE016
the rotational inertia of the thrombus elasticity measuring device;
Figure 880843DEST_PATH_IMAGE001
in order to measure the angular displacement of the thrombus elasticity measuring device,
Figure DEST_PATH_IMAGE017
calculating the frequency response function of the thrombus elasticity measuring device according to the motion differential equation as a unit impact function
Figure 502318DEST_PATH_IMAGE018
In the formula:
Figure DEST_PATH_IMAGE019
，
Figure 900545DEST_PATH_IMAGE020
，
Figure DEST_PATH_IMAGE021
(ii) a Using the equation of the frequency response function to the complex signal
Figure 949272DEST_PATH_IMAGE008
Performing least square fitting on the Nyquist curve to obtain detection values of the natural frequency and the damping ratio of the thrombus elasticity measuring device;
and comparing theoretical values of the natural frequency and the damping ratio with detection values of the natural frequency and the damping ratio, and judging whether the state of the thrombus elasticity measuring device needs to be corrected or not according to comparison of the variation and a preset tolerance value.
2. The method for detecting a thrombus elasticity measuring apparatus according to claim 1, wherein: said frequency response function
Figure 521199DEST_PATH_IMAGE022
Has a real part of
Figure DEST_PATH_IMAGE023
(ii) a Imaginary part of
Figure 306883DEST_PATH_IMAGE024
(ii) a In the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE025
。
3. the method for detecting a thrombus elasticity measuring apparatus according to claim 2, wherein:
Figure 202027DEST_PATH_IMAGE026
in the vicinity of the location of the mobile station,
Figure DEST_PATH_IMAGE027
，
Figure 31443DEST_PATH_IMAGE028
thus, therefore, it is
Figure DEST_PATH_IMAGE029
；
Figure 504713DEST_PATH_IMAGE030
At this time, the process of the present invention,
Figure DEST_PATH_IMAGE031
and
Figure 592623DEST_PATH_IMAGE032
satisfy the requirement of
Figure DEST_PATH_IMAGE033
The equation is a circular equation with the coordinates of the center of the circle
Figure 578159DEST_PATH_IMAGE034
Radius of
Figure DEST_PATH_IMAGE035
。
4. The method for detecting a thrombus elasticity measuring apparatus according to claim 3, wherein: fitting complex signals using Nyquist circles
Figure 844056DEST_PATH_IMAGE008
Solving min back face by least square method
Figure 374263DEST_PATH_IMAGE036
；
Figure DEST_PATH_IMAGE037
，
Figure 16597DEST_PATH_IMAGE038
；
Figure DEST_PATH_IMAGE039
；
Figure 119289DEST_PATH_IMAGE040
In order to be able to sample the frequency,
Figure DEST_PATH_IMAGE041
，
Figure 946299DEST_PATH_IMAGE005
the detection values of the natural frequency and the damping ratio of the system can be determined according to the sampling interval.
5. The method for detecting a thrombus elasticity measuring apparatus according to claim 3, wherein: connecting plural signals by fold lines
Figure 839169DEST_PATH_IMAGE008
At the intersection of the curve and the imaginary axis
Figure 973609DEST_PATH_IMAGE026
Real part of
Figure 510901DEST_PATH_IMAGE042
Imaginary part
Figure DEST_PATH_IMAGE043
The detected values of the natural frequency and the damping ratio of the system can be determined accordingly.
6. The method for detecting a thrombus elasticity measuring apparatus according to claim 1, wherein: will vibrate the signal
Figure 367867DEST_PATH_IMAGE007
Obtained by Fourier transform
Figure 623399DEST_PATH_IMAGE008
The method specifically comprises the following steps:
Figure 298881DEST_PATH_IMAGE044
；
Figure 956258DEST_PATH_IMAGE038
；
Figure 62755DEST_PATH_IMAGE039
；
Figure 54850DEST_PATH_IMAGE040
in order to be able to sample the frequency,
Figure 780361DEST_PATH_IMAGE041
，
Figure 42977DEST_PATH_IMAGE005
is the sampling interval time.
7. According to claimThe method for detecting a thrombus elasticity measuring apparatus according to claim 1, characterized in that: for original signal
Figure 320375DEST_PATH_IMAGE003
The pretreatment specifically comprises the following steps: from the original signal
Figure 550499DEST_PATH_IMAGE003
Subtracting the least squares best fit line from the original signal
Figure 328968DEST_PATH_IMAGE003
Removing linear trend in data to obtain vibration signal
Figure 695358DEST_PATH_IMAGE007
Capable of concentrating the analysis on the vibration signal
Figure 766826DEST_PATH_IMAGE007
On the fluctuation of (c).
8. The method for detecting a thrombus elasticity measuring apparatus according to claim 6, wherein: the least squares best fit line equation is
Figure DEST_PATH_IMAGE045
Wherein, in the step (A),
Figure 280984DEST_PATH_IMAGE046
，
Figure DEST_PATH_IMAGE047
，
Figure 269669DEST_PATH_IMAGE048
，
Figure DEST_PATH_IMAGE049
，
Figure DEST_PATH_IMAGE051
。
9. the method for detecting a thrombus elasticity measuring apparatus according to claim 1, wherein: when no blood is in the measuring cup,
Figure 975719DEST_PATH_IMAGE052
，
Figure DEST_PATH_IMAGE053
the rigidity of the thrombus elasticity measuring apparatus
Figure 922815DEST_PATH_IMAGE054
Moment of inertia
Figure 627333DEST_PATH_IMAGE016
Calculated by three-dimensional drawing software.
10. The method for detecting a thrombus elasticity measuring apparatus according to claim 1, wherein: the detection method of the thrombus elasticity measuring device further comprises a correction step, wherein the correction step specifically comprises the following steps: and calculating a correction coefficient by measuring the frequency variation, and correcting the measurement error caused by the rigidity variation on line.
11. The method for detecting a thrombus elasticity measuring apparatus according to claim 9, wherein: correction factor
Figure DEST_PATH_IMAGE055
Amount of variation
Figure 560654DEST_PATH_IMAGE056
Is the difference between the rigidity detection value and the theoretical value of the thrombus elasticity measuring device,
Figure DEST_PATH_IMAGE057
moment of inertia
Figure 823008DEST_PATH_IMAGE016
Calculated by three-dimensional drawing software.
CN202210507792.6A 2022-05-11 2022-05-11 Detection method of thrombus elasticity measuring device Active CN114608995B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210507792.6A CN114608995B (en) 2022-05-11 2022-05-11 Detection method of thrombus elasticity measuring device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210507792.6A CN114608995B (en) 2022-05-11 2022-05-11 Detection method of thrombus elasticity measuring device

Publications (2)

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

Family

ID=81870471

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210507792.6A Active CN114608995B (en) 2022-05-11 2022-05-11 Detection method of thrombus elasticity measuring device

Country Status (1)

Country Link
CN (1) CN114608995B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1226901A (en) * 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 (en) * 2016-08-30 2017-02-22 诺泰科生物科技（苏州）有限公司 Thrombelastometer and calibration method thereof
CN109520889A (en) * 2012-01-16 2019-03-26 仪宝科技公司 For measuring the method, apparatus and system of fluid physics performance
CN110547825A (en) * 2019-09-03 2019-12-10 上海交通大学 Ultrasonic elastography technology based on multi-frequency shear wave amplitude analysis
CN113640377A (en) * 2021-10-15 2021-11-12 中国科学院苏州生物医学工程技术研究所 Blood viscoelasticity measuring method based on vibration sensor

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1226901A (en) * 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 (en) * 2012-01-16 2019-03-26 仪宝科技公司 For measuring the method, apparatus and system of fluid physics performance
CN106442952A (en) * 2016-08-30 2017-02-22 诺泰科生物科技（苏州）有限公司 Thrombelastometer and calibration method thereof
CN110547825A (en) * 2019-09-03 2019-12-10 上海交通大学 Ultrasonic elastography technology based on multi-frequency shear wave amplitude analysis
CN113640377A (en) * 2021-10-15 2021-11-12 中国科学院苏州生物医学工程技术研究所 Blood viscoelasticity measuring method based on vibration sensor

Also Published As

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

Similar Documents

Publication Publication Date Title
CN115728063A (en) Flexible bearing detection method
EP2136189B1 (en) Method for analysing vibration in rotor blades
US10317275B2 (en) Vibration monitoring systems
CN103038620B (en) The apparatus for diagnosis of abnormality of rolling bearing
CN110567574A (en) Method and system for identifying timing vibration parameters of blade end of rotating blade
KR19990087391A (en) Multiple sensor device and method for monitoring turbomachinery
CN111183345B (en) Method and device for searching for defects that can affect a rotary machine power transmission device
CN114608995B (en) Detection method of thrombus elasticity measuring device
CN109580924A (en) The frequency domain parameter analysis method and system of thrombelastogram process of setting
CN111649886B (en) Abnormality detection device, rotating machine, abnormality detection method, and computer-readable storage medium
JP4424515B2 (en) Method for detecting solid conduction sound in rolling bearings
CN109940456A (en) Machine tool chief axis method for testing performance and system
US5710715A (en) Vibration analysis method
CN108151870B (en) Construction quality problem detection method based on frequency response function
JPH03269221A (en) Abnormal-sound diagnostic apparatus for rotary equipment
CN108037316B (en) Equipment and method for evaluating reliability of accelerometer based on performance index
CN109556895A (en) The failure analysis methods and device of rotating machinery
Zamorano et al. Analysis in the time-frequency domain of different depths of a crack located in a change of section of a shaft
JP2020144111A (en) Abnormality detection device, rotary machine, abnormality detection method and program
JPH0450731A (en) Rotary machine fault diagnostic system
RU2682561C1 (en) Method for determining technical condition of current collectors
CN110307174A (en) A kind of cleaning fan manufacture assembling quality detection device and method
KR102219422B1 (en) Detecting system for measuring operating sounds of washing machine and detecting method for noise defect there of
CN110988526B (en) Robot assembly inspection method and device and storage medium
Le et al. Dynamic balancing of rigid rotors by the influence coefficient method

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