CN106501182A - Method of nondestructive elasticity measurement by photoacoustic eigen spectrum analysis - Google Patents

Method of nondestructive elasticity measurement by photoacoustic eigen spectrum analysis Download PDF

Info

Publication number
CN106501182A
CN106501182A CN201610844345.4A CN201610844345A CN106501182A CN 106501182 A CN106501182 A CN 106501182A CN 201610844345 A CN201610844345 A CN 201610844345A CN 106501182 A CN106501182 A CN 106501182A
Authority
CN
China
Prior art keywords
elastomer
eigenfrequency
step
method
equation
Prior art date
Application number
CN201610844345.4A
Other languages
Chinese (zh)
Other versions
CN106501182B (en
Inventor
陶超
郜晓翔
刘晓峻
王学鼎
Original Assignee
南京大学
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 南京大学 filed Critical 南京大学
Priority to CN201610844345.4A priority Critical patent/CN106501182B/en
Priority claimed from CN201610844345.4A external-priority patent/CN106501182B/en
Publication of CN106501182A publication Critical patent/CN106501182A/en
Application granted granted Critical
Publication of CN106501182B publication Critical patent/CN106501182B/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using infra-red, visible or ultra-violet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/1702Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using infra-red, visible or ultra-violet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/1702Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
    • G01N2021/1706Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids in solids

Abstract

The invention discloses a method of nondestructive elasticity measurement by photoacoustic eigen spectrum analysis. The method provided by the invention utilizes laser pulse to irradiate an elastomer to make it radiate sound wave under the action of photoacoustic effect, and through analysis of the photoacoubtic signal radiated by the elastomer, a photoacoustic eigen spectrum of the elastomer can be obtained, the eigen frequency of the elastomer can be extracted from the photoacoustic eigen spectrum, on the basis of the eigen frequency, inversion algorithm is also utilized to evaluate the elastic properties of the elastomer. The method of nondestructive elasticity measurement by photoacoustic eigen spectrum analysis provided by the invention has no need for mutual contact of an acoustic transducer and the elastomer during measurement, and also has no need for any dissection treatment on the elastomer. The method can realize non-contact and nondestructive evaluation of elastic parameters, and has high safety and usability.

Description

一种利用光声本征谱分析法无损测量弹性的方法 A method of spectral analysis of the intrinsic elasticity of the nondestructive measurement using the photoacoustic

技术领域 FIELD

[0001] 本发明涉及一种弹性性质的无损测量方法,具体说是一种利用光声效应以及最优化算法无损测量弹性的方法。 [0001] The present invention relates to a method for nondestructive measurement of elastic properties, in particular to a photoacoustic effect, and a method using a non-destructive measurement algorithm optimized elasticity.

背景技术 Background technique

[0002] 非侵入性地、非接触式地测量弹性不仅仅是材料科学的基础研究内容,在生物医学领域也有重要的意义。 [0002] non-invasive, non-contact measurement of elastic materials science basic research is not only content, but also in the field of biomedical importance. 在许多疾病中,包括动脉硬化、肝硬化、肿瘤等,组织弹性是一个重要的生物医学指标。 In many diseases, including atherosclerosis, cirrhosis, cancer and other tissue elasticity is an important biomedical indicators. 超声共振谱分析法可用于测量物体弹性,但一般情况下需要样品与测量换能器接触。 Ultrasonic resonance spectral analysis can be used to measure the elastic body, but generally the need to sample the measurement transducer contacts. 原子力显微镜是一种能对局部弹性精确测量的技术,但也仅局限于物体的表面或者亚表面检测。 It is an atomic force microscope capable of accurate measurement of the local elastic technology, but is limited to the surface or sub-surface detection object.

[0003] 基于光声效应的技术在材料检测方面已经引起广泛的关注。 [0003] Based on the photoacoustic effect technology has attracted attention in terms of the detection material. 尤其是在过去一二十年内,光声成像由于高对比度、深层组织高分辨率、高安全性这些优点,已经在生物医学领域被广泛应用。 Especially in the past two decades, photoacoustic imaging due to the high contrast, deep tissue high-resolution, high-security these advantages, has been widely used in biomedical applications. 利用光声技术,我们可以成功地提取许多生物信息,包括血氧饱和度、血流速度、温度、粘弹性等等。 Using photoacoustic technique, we can successfully extracted a number of biological information, including oxygen saturation, blood flow rate, temperature, viscoelasticity and the like.

发明内容 SUMMARY

[0004] 发明目的:为了克服现有技术中存在的不足,本发明提供一种利用光声本征谱分析法无损测量弹性的方法,该方法利用光声效应检测弹性体的本征频率,再通过最优化算法反演计算出弹性参数,提供了一种非侵入性的、非接触式的测量弹性的方法。 [0004] Object of the invention: In order to overcome the disadvantages present in the prior art, the present invention provides a method using photoacoustic spectroscopy intrinsic elasticity nondestructive measurement method, the method utilizes a photoacoustic effect eigenfrequency detected elastomer, then by optimizing the inversion algorithm to calculate the elastic parameters, there is provided a non-invasive, non-contact method of measuring the elastic type.

[0005] 技术方案:为实现上述目的,本发明采用的技术方案为: [0005] Technical Solution: To attain the above object, the technical solution adopted by the invention is:

[0006] -种利用光声本征谱分析法无损测量弹性的方法,包括以下步骤: [0006] - methods using photoacoustic spectroscopy intrinsic elasticity nondestructive measurement method, comprising the steps of:

[0007] 步骤1,使用脉冲激光器对待测弹性体照射脉冲激光。 [0007] Step 1, measured using a pulsed laser to treat elastomeric pulsed laser. 待测弹性体在光声效应激发后辐射声波。 Elastomeric acoustic radiation measured after excitation photoacoustic effect.

[0008] 步骤2,利用超声换能器接收待测弹性体辐射出的声波信号。 [0008] Step 2, using the ultrasonic transducer receives the acoustic test signal radiated elastomer.

[0009] 步骤3,对步骤2中超声换能器采集的声波信号进行处理,计算出声波信号的时频图。 [0009] Step 3, step 2 ultrasonic transducer acoustic signal acquisition processing, calculating the spectrogram of the acoustic signal. 在时频图的后段得到待测弹性体的本征频率谱线,提取本征频率谱线得到待测弹性体的本征频率,该本征频率为测得的本征频率。 When the frequency of the segment of FIG eigenfrequency spectrum obtained test elastomer obtained extract eigenfrequency spectrum measured eigenfrequency elastomer, the intrinsic frequency of the measured eigenfrequency.

[0010] 步骤4,根据弹性体密度以及质点位移矢量建立弹性体满足的波动方程及其边界条件。 [0010] Step 4, to establish the elastic wave equation and the boundary conditions of the elastic member satisfy the bulk density and the particle displacement vector. 根据边界条件得到弹性体振动的特征方程,求解特征方程可得到弹性体的本征频率, 该本征频率为特征方程计算出的本征频率。 Boundary conditions vibration characteristic equation of the elastic body according to solving the characteristic equation is obtained elastomer eigenfrequency, the eigenfrequency of the characteristic equation of the eigenfrequency.

[0011] 步骤5,根据步骤3中测得的本征频率和步骤4中特征方程计算出的本征频率建立误差函数F(E,〇): [0011] Step 5 is calculated based on Equation 4 in step 3 wherein measured eigenfrequency and establishing steps eigenfrequency error function F (E, square):

Figure CN106501182AD00041

[0012] (1) [0012] (1)

[0013] 兵ψ,b73物K候松tQ止频竿的个数,gi为测得的本征频率,fi为特征方程计算出的本征频率,Wi为加权系数。 [0013] Bing ψ, b73 K was phenology tQ loose stopper pole number frequency, gi is the measured intrinsic frequency, fi is the characteristic equation calculated eigenfrequencies, Wi is a weighting factor. 当误差函数F (;Ε,σ)取最小值时,利用Levenberg- Marquardt算法,反演得到的E和σ为评估值,求得反演结果。 When the error function F.; Time (Ε, σ) takes a minimum value, using Levenberg- Marquardt algorithm, inversion and [sigma] is the evaluation value E, determined inversion results.

[0014] 所述步骤4中特征方程计算出的本征频率的获取方法如下: [0014] The method of obtaining the calculated eigenfrequency wherein said step of Equation 4 as follows:

[0015] 步骤41,弹性体的波动方程满足如下公式: [0015] Step 41, the wave equation elastomer satisfies the following formula:

[0016 [0016

Figure CN106501182AD00051

[0017] 其中,Ps为弹性体密度,E为杨氏模量,〇为泊松比,U= (x,y,z)为质点位移矢量, ▽ = (3/θχ, (5/¾) , t表示时间,X、y、z表示质点空间坐标。 [0017] wherein, Ps is the density of the elastomer, E is Young's modulus, Poisson's ratio is square, U = (x, y, z) is the particle displacement vector, ▽ = (3 / θχ, (5 / ¾) , t represents time, X, y, z spatial coordinates represents a particle.

[0018] 同时,处于无粘性流体中的球形弹性体,在边界上满足3个边界条件:i)球体内的法向应力与流体内的压强相等。 [0018] Meanwhile, in the spherical elastomeric inviscid fluid, three boundary conditions on the boundary: i) the normal stress and pressure of the fluid inside the sphere equal. ii)球体内的法向位移与流体内的法向位移相等。 ii) the normal sphere fluid displacement equal to the displacement method. iii)球体内切变应力的切向成分为零。 iii) cutting a sphere to zero shear stress component.

[0019] 步骤42,根据步骤41中的边界条件,可以得到弹性体振动的特征方程,所述弹性体振动的特征方程为|D|=0,其中,D= {du},1 = 1,2或3,J=1,2或3, [0019] Step 42, based on the boundary conditions in step 41, it is possible to obtain the characteristic equation of the elastic vibration member, wherein the elastic body vibration equation is | D | = 0, where, D = {du}, 1 = 1, 2 or 3, J = 1,2 or 3,

Figure CN106501182AD00052

[0024] 其中,Pf为周围介质密度,PsS球形弹性体的密度,a为球形弹性体半径,1^和九分别表示η阶汉克尔函数和η阶球贝塞尔函数,η表示汉克尔函数和球贝塞尔函数的阶数;kt = 2JTf/ct,kf = 2JTf/Cf,ki = 2JTf/Cl为波数,f为振动频率;Cf,Ct,Cl分别为介质的声速、弹性体的横波速度和纵波速度。 [0024] where, Pf is the density of the surrounding medium, the spherical elastomeric PsS density, a is the radius of the spherical elastomer, and nine, respectively represent 1 ^ [eta] [eta] order Hankel function and order spherical Bessel functions, [eta] represents Hank Er order spherical Bessel functions and the functions; kt = 2JTf / ct, kf = 2JTf / Cf, ki = 2JTf / Cl is the wave number, f is the frequency of vibration; Cf, Ct, Cl sound velocity medium, respectively, the elastic body the P-wave velocity and shear wave velocity.

[0025] 步骤43,根据弹性体密度求解步骤42中的特征方程可得到弹性体的本征频率,该本征频率即为特征方程计算出的本征频率A。 [0025] Step 43, the step of solving the elastic body 42 can be obtained density eigenfrequency elastomeric characteristic equation, which is the characteristic eigenfrequency of the equation eigenfrequency A.

[0026] 优选的:所述脉冲激光器的脉宽8_12ns。 [0026] Preferably: the pulse width of the pulsed laser 8_12ns.

[0027] 优选的:所述步骤2中超声换能器采集到的声波信号通过小信号放大器放大以及数字采集卡采样之后存储在计算机中。 [0027] Preferred: Step 2 of the acoustic signal in the ultrasonic transducer is collected through a small signal amplifier amplifying the digital acquisition card samples and thereafter stored in the computer.

[0028] 优选的:所述加权系数取fl/f。 [0028] The preferred: the weighting coefficient takes fl / f.

[0029] 有益效果:本发明与超声共振谱分析法、原子力显微镜相比,具有以下优点: [0029] The beneficial effects: Compared with the ultrasonic resonance spectral analysis, atomic force microscope, has the following advantages:

[0030] (1)本发明利用光声效应激发待测物体振动,然后利用声学换能器接收物体辐射出的声波信号进行分析。 [0030] (1) The present invention utilizes a photoacoustic effect excited vibration object to be measured, and then using the acoustic transducer receiving acoustic signals radiated from the object is analyzed. 整个过程不需要物体与换能器接触,所以此方法是非接触式的。 The whole process does not require contact with the object with the transducer, so this method is contactless.

[0031] (2)本发明不需要换能器与待测物体相互接触,所以可以方便地实现非侵入性的深层组织的检测。 [0031] (2) detecting the present invention does not need the transducer and object to be measured in contact with each other, it is possible to easily achieve non-invasive deep tissue.

附图说明 BRIEF DESCRIPTION

[0032]图1本发明光声本征谱分析法无损测量弹性的方法的系统示意图。 [0032] The system diagram of FIG method for measuring a light elastic acoustic present invention Nondestructive intrinsic spectral analysis.

[0033]图2金属球样品的本振频率提取示意图,图2a为黄铜球受激光激发后辐射出的声波信号图,图2b为整段声信号的时频图,图2c为头波与尾波的频谱图。 [0033] FIG. 2 the local oscillator frequency sample extraction schematic metal balls, Figure 2a is a view of a brass ball acoustic signal radiated by the laser excitation, Figure 2b is the whole FIG frequency acoustic signals, and Figure 2c for the first wave spectrum of coda.

[0034]图3弹性评估值与实际值比较。 [0034] FIG. 3 elasticity evaluation value compared with the actual value.

具体实施方式 Detailed ways

[0035] 下面结合附图和具体实施例,进一步阐明本发明,应理解这些实例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。 [0035] The accompanying drawings and the following specific examples further illustrate the invention, it should be understood that these examples are illustrative only of the present invention is not intended to limit the scope of the present invention, after reading this disclosure, those skilled in the art of the present invention. various equivalent modifications are forms fall within the present application as defined in the appended claims scope.

[0036] -种利用光声本征谱分析法无损测量弹性的方法,以测量一直径为1.0mm的黄铜球的弹性参数为例,具体包括以下步骤: [0036] - Method photoacoustic spectroscopy eigen analysis method using the nondestructive measurement of elastic kind, to measure a 1.0mm diameter brass ball elastic parameters as an example, includes the following steps:

[0037] 步骤1,如图1所示,一个直径为1 · Omm的黄铜球埋入圆柱体状琼脂中,深度为2厘米,琼脂半径约为2厘米。 [0037] Step 1 shown in Figure 1, having a diameter of 1 · Omm brass balls embedded in agar cylindrical shape, a depth of 2 cm and a radius of about 2 cm agar. 利用Nd: YAG脉冲激光器对待测弹性体(黄铜球)照射脉冲激光,激光脉宽约为l〇ns,并由激光器给数字采集卡一个触发;脉冲激光器的脉宽在8-12ns即可。 Using a pulsed Nd: YAG laser treatment test elastomer (brass ball) irradiating the pulsed laser, laser pulse duration of about l〇ns by the laser to trigger a data acquisition card; laser pulse width can be in 8-12ns. 待测弹性体在光声效应激发后辐射声波。 Elastomeric acoustic radiation measured after excitation photoacoustic effect.

[0038] 步骤2,黄铜球辐射声波之后,由利用超声换能器接收待测弹性体辐射出的声波信号。 [0038] Step 2, after the brass ball radiating sound waves, acoustic signals received radiation measured by the elastic body using an ultrasonic transducer. 超声换能器中心频率为4.39MHz,-6dB带宽为4.4MHz。 Ultrasonic transducer center frequency of 4.39MHz, -6dB bandwidth of 4.4MHz. 声波信号经过信号放大器(小信号放大器)放大,并经采样率为60MHz的采集卡采样,最终保存到计算机中。 Acoustic signal through a signal amplifier (small signal amplifier) ​​amplified and sampled by the sampling rate is 60MHz acquisition card, eventually saved to the PC.

[0039]如图2所示,图2a表示黄铜球受激光激发后辐射出的声波信号。 [0039] As shown, FIG. 2a represents the brass ball receiving acoustic signals radiated laser excitation. 信号的前部幅度较强,这段信号称为头波。 The front portion of the amplitude of the signal is strong, this is called the first signal wave. 头波是金属球直接受激光激励辐射出的光声信号,头波的放大图如图所示。 The direct wave is a metal ball head photoacoustic signal radiated by laser excitation, as shown in an enlarged view of the head wave in FIG. 头波之后,是一段时间相对较长、幅度不断衰减的声信号,这段信号称为尾波。 After the first wave, is the relatively long period of time, the amplitude of the continuous acoustic signal is attenuated, this signal is called coda. 尾波是在激光消失之后,金属球由于自身惯性保持振动辐射出的声波。 Coda are disappeared after the laser, due to its inertia metal ball holding radiated acoustic vibration.

[0040] 整段声信号的时频图如图2b所示。 [0040] When the whole audio frequency signal as shown in FIG. 2b. 时频图前段宽带的信号对应于头波,而后段较多的窄带共振谱线则对应于尾波。 FIG preceding wideband time-frequency signals corresponding to the first wave, then more segments corresponding to the narrow-band resonance line coda. 尾波信号在频域上由多个本征频率组成,表明尾波信号不是随机的,而是周期的。 Coda signal in the frequency domain by a plurality of eigenfrequencies composition showed coda signal is not random, but the cycle. 尾波信号的共振谱线即对应金属球的本征频率。 Coda resonance line signal corresponding to the eigenfrequency, i.e. metal balls.

[0041] 如图2c为头波与尾波的频谱。 [0041] Figure 2c for the first and coda-wave spectrum. 头波为宽带谱,而尾波由多个共振峰组成,即对应时频图中的共振谱线。 A broad-spectrum of the first wave, and a plurality of formant coda composition, the resonance lines in FIG i.e. the corresponding frequency.

[0042] 步骤3,通过傅里叶变换对步骤2中超声换能器采集的声波信号进行处理,计算出声波信号的时频图。 When [0042] Step 3, the processing of acoustic signals by a Fourier transform step 2 ultrasonic transducer acquired, the calculated frequency of the acoustic signal in FIG. 在时频图的后段(尾部)得到待测弹性体的本征频率谱线,提取本征频率谱线得到待测弹性体的本征频率,该本征频率为测得的本征频率。 When the rear section of FIG frequency (tails) obtained eigenfrequency spectrum measured elastomer extraction eigenfrequency measured spectrum obtained eigenfrequency elastomer, the intrinsic frequency of the measured eigenfrequency.

[0043] 步骤4,根据弹性体密度以及质点位移矢量建立弹性体满足的波动方程及其边界条件。 [0043] Step 4, to establish the elastic wave equation and the boundary conditions of the elastic member satisfy the bulk density and the particle displacement vector. 根据边界条件得到弹性体振动的特征方程,求解特征方程可得到弹性体的本征频率, 该本征频率为特征方程计算出的本征频率。 Boundary conditions vibration characteristic equation of the elastic body according to solving the characteristic equation is obtained elastomer eigenfrequency, the eigenfrequency of the characteristic equation of the eigenfrequency.

[0044] 步骤41,弹性体的波动方程满足如下公式: [0044] Step 41, the wave equation elastomer satisfies the following formula:

[0045] [0045]

Figure CN106501182AD00061

[0046] 其中,PsS弹性体密度,E为杨氏模量,σ为泊松比,U= (x,y,z)为质点位移矢量, ▽ = (3/次r, 表示时间,x、y、z表示质点空间坐标。 [0046] wherein, PsS elastomer density, E is Young's modulus, σ is Poisson's ratio, U = (x, y, z) is the particle displacement vector, ▽ = (3 / time r, represents time, x, y, z spatial coordinates represents a particle.

[0047] 同时,处于无粘性流体中的球形弹性体,在边界上满足3个边界条件:i)球体内的法向应力与流体内的压强相等。 [0047] Meanwhile, in the spherical elastomeric inviscid fluid, three boundary conditions on the boundary: i) the normal stress and pressure of the fluid inside the sphere equal. ii)球体内的法向位移与流体内的法向位移相等。 ii) the normal sphere fluid displacement equal to the displacement method. iii)球体内切变应力的切向成分为零。 iii) cutting a sphere to zero shear stress component.

[0048] 步骤42,根据步骤41中的边界条件,可以得到弹性体振动的特征方程,所述弹性体振动的特征方程为|D|=0,其中,D= {du},1 = 1,2或3,J=1,2或3。 [0048] Step 42, based on the boundary conditions in step 41, it is possible to obtain the characteristic equation of the elastic vibration member, wherein the elastic body vibration equation is | D | = 0, where, D = {du}, 1 = 1, 2 or 3, J = 1,2 or 3.

Figure CN106501182AD00071

[0049] [0049]

[0050] [0050]

[0051] [0051]

[0052] [0052]

[0053] 其中,Pf为周围介质密度,Ps为球形弹性体的密度,a为球形弹性体半径,1^和九分别表示η阶汉克尔函数和η阶球贝塞尔函数,η表示汉克尔函数和球贝塞尔函数的阶数,也就是说汉克尔函数和球贝塞尔函数是同阶的;kt = 2iif/ct,kf = 2iif/cf,ki = 2iif/ci为波数,f为振动频率;Cf,Ct,C1分别为介质的声速、弹性体的横波速度和纵波速度。 [0053] where, Pf is the density of the surrounding medium density, Ps spherical elastomer, a is the radius of the spherical elastomer, and nine, respectively represent 1 ^ [eta] [eta] order Hankel function and order spherical Bessel functions, [eta] represents Han Kerr order spherical Bessel functions and function, i.e. spherical Bessel functions and Hankel function is of the same order; kt = 2iif / ct, kf = 2iif / cf, ki = 2iif / ci is the wave number , f is the frequency of vibration; Cf, Ct, C1 sound velocity medium, respectively, the longitudinal wave velocity and shear wave velocity elastomer.

[0054] 步骤43,根据弹性体密度求解步骤42中的特征方程可得到弹性体的本征频率,该本征频率即为特征方程计算出的本征频率A。 [0054] Step 43, the step of solving the elastic body 42 can be obtained density eigenfrequency elastomeric characteristic equation, which is the characteristic eigenfrequency of the equation eigenfrequency A.

[0055] 步骤5,根据步骤3中测得的本征频率和步骤4中特征方程计算出的本征频率建立误差函数F(E,〇): [0055] Step 5 is calculated based on Equation 4 in step 3 wherein measured eigenfrequency and establishing steps eigenfrequency error function F (E, square):

[0056] [0056]

Figure CN106501182AD00072

(2) (2)

[0057] 其中,E为杨氏模量,〇为泊松比,N为本征频率的个数,gi为测得的本征频率,fi为特征方程计算出的本征频率,Wi为加权系数,加权系数一般取E=IZa2s当误差函数F (E,〇)取最小值时,利用Levenberg-Marquardt算法,反演得到的E和σ为评估值,求得反演结果。 [0057] wherein, E is Young's modulus, Poisson's ratio for the square, the number N of intrinsic frequency, gi eigenfrequency is measured, calculated as Fi eigenfrequency characteristic equation, Wi is a weighting factor, the weighting factor generally taken when E = IZa2s error function F. (E, square) takes a minimum value, using the Levenberg-Marquardt algorithm, inversion E and σ is an evaluation value, inversion results obtained.

[0058] 将图2中求得的本征频率gi代入公式(2),利用Levenberg-Marquardt算法,求得误差函数最小时的杨氏模量以及泊松比。 [0058] FIG. 2 is obtained in the eigenfrequency gi into the formula (2), using the Levenberg-Marquardt algorithm, to obtain an error ratio of the Young's modulus and Poisson's function is minimum. 杨氏模量以及泊松比的评估值和实际值,如图3所示。 Young's modulus and Poisson's evaluation value and the actual value of the ratio, as shown in FIG. 结果表明,光声本征谱分析法无损测量弹性的方法具有较高的准确性。 The results show that the method nondestructive measurement of the intrinsic elasticity of the photoacoustic spectrum method with high accuracy.

[0059] 本发明利用激光脉冲照射弹性体使其在光声效应作用下辐射声波,通过分析弹性体辐射出的光声信号,可以得到弹性体的光声本征谱,从光声本征谱中提取出弹性体的本征频率,基于这些本征频率,再利用反演算法可以评估弹性体的弹性性质。 [0059] The present invention is irradiated with a laser pulse radiating sound waves in the elastic body so that the action of the photoacoustic effect, by analyzing light radiated elastomeric acoustic signals can be obtained photoacoustic intrinsic spectrum elastomer intrinsic spectrum from the photoacoustic extracted eigenfrequency elastomer, based on these intrinsic frequency reuse inversion algorithm can evaluate the elastic properties of the elastomer. 本发明提出的利用光声本征谱分析法无损测量弹性的方法,无须声学换能器与弹性体在测量时相互接触, 更不需要对弹性体做任何解剖处理,此方法可实现对弹性参数的非接触、无损评估,具有较高的安全性和易用性。 The method of using photoacoustic spectroscopy intrinsic elasticity nondestructive measurement method proposed by the present invention, without the acoustic transducer and the elastomer in contact with each other in the measurement, but do not need to make any anatomical elastomer processing, this method can be realized on the elastic parameters non-contact, non-destructive evaluation, with high security and ease of use.

[0060] 以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。 [0060] The above are only preferred embodiments of the present invention, it should be noted: to those of ordinary skill in the art, in the present invention without departing from the principles of the premise, can make various improvements and modifications, such modifications and modifications should also be regarded as the protection scope of the present invention.

Claims (5)

1. 一种利用光声本征谱分析法无损测量弹性的方法,其特征在于,包括以下步骤: 步骤1,使用脉冲激光器对待测弹性体照射脉冲激光;待测弹性体在光声效应激发后辐射声波; 步骤2,利用超声换能器接收待测弹性体辐射出的声波信号; 步骤3,对步骤2中超声换能器采集的声波信号进行处理,计算出声波信号的时频图;在时频图的后段得到待测弹性体的本征频率谱线,提取本征频率谱线得到待测弹性体的本征频率,该本征频率为测得的本征频率; 步骤4,根据弹性体密度以及质点位移矢量建立弹性体满足的波动方程及其边界条件; 根据边界条件得到弹性体振动的特征方程,求解特征方程可得到弹性体的本征频率,该本征频率为特征方程计算出的本征频率; 步骤5,根据步骤3中测得的本征频率和步骤4中特征方程计算出的本征频率建立误差函数F(E,〇): 1. A method for photoacoustic spectroscopy eigen analysis method using the nondestructive measurement of elasticity, characterized by comprising the following steps: Step 1, a pulsed laser irradiation of a pulse laser measuring treated elastomer; measured at excitation elastomer photoacoustic effect after radiating sound waves; step 2, using the ultrasonic transducer receives the acoustic test signal radiated elastomer; step 3, step 2 ultrasonic transducer acoustic signal acquisition process is performed to calculate the acoustic signal in the frequency map; in FIG time-frequency segment obtained after the test elastomer eigenfrequency spectrum, obtained extract eigenfrequency spectrum measured eigenfrequency elastomer, the intrinsic frequency of the measured intrinsic frequency; step 4. the elastomer particle displacement vector building density, and the elastic wave equation and satisfies the boundary condition; the characteristic equation boundary conditions elastomeric vibration, solving the characteristic equation is obtained elastomer eigenfrequency, the eigenfrequency characterized equation the eigenfrequency; step 5, calculated from equation 4 in step 3 wherein measured eigenfrequency and establishing steps eigenfrequency error function F (E, square):
Figure CN106501182AC00021
其中,E为杨氏模量,σ为泊松比,N为本征频率的个数,gi为测得的本征频率,fi为特征方程计算出的本征频率,Wi为加权系数;当误差函数F (Ε,σ)取最小值时,利用Levenberg-Marquardt算法,反演得到的E和σ为评估值,求得反演结果。 Wherein, E is Young's modulus, Poisson's ratio of [sigma], the number N of intrinsic frequency, gi is the measured intrinsic frequency, Fi is the characteristic equation calculated eigenfrequency, Wi is a weighting coefficient; if error function F (Ε, σ) takes a minimum value when using the Levenberg-Marquardt algorithm, inversion and [sigma] is the evaluation value E, determined inversion results.
2. 根据权利要求1所述的利用光声本征谱分析法无损测量弹性的方法,其特征在于:所述步骤4中特征方程计算出的本征频率的获取方法如下: 步骤41,弹性体的波动方程满足如下公式: 2. The method of claim 1 using photoacoustic spectroscopy analysis of the intrinsic elasticity of claim nondestructive measurement, wherein: said step of acquiring method characterized Equation 4 eigenfrequency calculated as follows: Step 41, the elastic member the wave equation satisfy the following equation:
Figure CN106501182AC00022
其中,Ps为弹性体密度,Ε为杨氏模量,σ为泊松比,U= (x,y,z)为质点位移矢量,x、y、z表示质点空间坐标,V = (θ/δχ, δ/狀3/3z) , t表示时间; 同时,处于无粘性流体中的球形弹性体,在边界上满足3个边界条件:i)球体内的法向应力与流体内的压强相等;ii)球体内的法向位移与流体内的法向位移相等;iii)球体内切变应力的切向成分为零; 步骤42,根据步骤41中的边界条件,可以得到弹性体振动的特征方程,所述弹性体振动的特征方程为|D|=0,其中,D= {如},1 = 1,2或3,J=1,2或3; Wherein, Ps is the density of the elastomer, Ε Young's modulus, σ particles representing spatial coordinates, V = (θ is the particle displacement vector, x, y, z is Poisson's ratio, U = (x, y, z) / δχ, δ / shape 3 / 3z), t represents time; while, in the spherical elastomeric inviscid fluid, to meet the three boundary condition on the boundary: i) within the sphere method is equal to the pressure stress within the fluid; ii) the normal sphere of displacement fluid displaced is equal to the method; iii) cutting sphere to zero shear stress component; step 42, the boundary conditions in step 41, it is possible to obtain elastomeric vibration characteristic equation the elastomeric vibration characteristic equation of | D | = 0, where, D = {as}, = 1, 2 or 3, J = 1, 2 or 3;
Figure CN106501182AC00023
其中,Pf为周围介质密度,PS为球形弹性体的密度,a为球形弹性体半径,1^和九分别表示η阶汉克尔函数和η阶球贝塞尔函数,η表示汉克尔函数和球贝塞尔函数的阶数;kt = 2Jif/ct, kf = 2Jif/Cf,ki = 2Jif/ci为波数,f为振动频率;Cf,ct,ci分别为介质的声速、弹性体的横波速度和纵波速度; 步骤43,根据弹性体密度求解步骤42中的特征方程可得到弹性体的本征频率,该本征频率即为特征方程计算出的本征频率h。 Wherein, Pf is the density of the surrounding medium, PS is the density of the spherical elastomer, a is the radius of the spherical elastomer, and nine, respectively represent 1 ^ [eta] [eta] order Hankel function and order spherical Bessel functions, [eta] represents a Hankel function and spherical Bessel functions of order; kt = 2Jif / ct, kf = 2Jif / Cf, ki = 2Jif / ci is the wave number, f is the frequency of vibration; Cf, ct, ci are transverse acoustic velocity medium, the elastomer velocity and compressional wave velocity; step 43, the step of solving the density of the elastic body 42 in the characteristic equation is obtained elastomer eigenfrequency, the eigenfrequency of the equation is the characteristic eigenfrequency h.
3. 根据权利要求1所述的利用光声本征谱分析法无损测量弹性的方法,其特征在于:所述脉冲激光器的脉宽8-12ns。 The use of photoacoustic spectroscopy eigen analysis method according to a non-destructive method for measuring the elastic claim, wherein: said pulse width 8-12ns laser.
4. 根据权利要求1所述的利用光声本征谱分析法无损测量弹性的方法,其特征在于:所述步骤2中超声换能器采集到的声波信号通过小信号放大器放大以及数字采集卡采样之后存储在计算机中。 4. The method of claim 1 using photoacoustic spectroscopy analysis of the intrinsic elasticity of claim nondestructive measurement, wherein: said step 2 ultrasonic transducer acoustic signal acquired by amplifying a small signal amplifier and a digital acquisition card after the samples stored in the computer.
5. 根据权利要求1所述的利用光声本征谱分析法无损测量弹性的方法,其特征在于:所述加权系数取wi=l/gi2。 5. A non-destructive method for measuring the elastic claim eigen using photoacoustic spectrum method of claim 1, wherein: the weighting coefficient takes wi = l / gi2.
CN201610844345.4A 2016-09-22 A method of utilizing the intrinsic spectrum analysis method nondestructive measurement elasticity of optoacoustic CN106501182B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610844345.4A CN106501182B (en) 2016-09-22 A method of utilizing the intrinsic spectrum analysis method nondestructive measurement elasticity of optoacoustic

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610844345.4A CN106501182B (en) 2016-09-22 A method of utilizing the intrinsic spectrum analysis method nondestructive measurement elasticity of optoacoustic

Publications (2)

Publication Number Publication Date
CN106501182A true CN106501182A (en) 2017-03-15
CN106501182B CN106501182B (en) 2019-06-04

Family

ID=

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6041020A (en) * 1997-04-21 2000-03-21 University Of Delaware Gas-coupled laser acoustic detection
EP1998161A1 (en) * 2007-05-29 2008-12-03 The Technical University of Denmark (DTU) Acoustic resonator cell for spectroscopic analysis of a compound in a fluid
CN102226752A (en) * 2011-04-08 2011-10-26 常熟舒茨电子科技发展有限公司 Method for detecting concentration of sulfur dioxide by using photoacoustic spectrometry
CN103324786A (en) * 2013-06-03 2013-09-25 北京工业大学 Two-dimensional solid-phase phonon crystal XY module cord gap optimization method based on genetic algorithm
US20140073907A1 (en) * 2012-09-12 2014-03-13 Convergent Life Sciences, Inc. System and method for image guided medical procedures
CN103884774A (en) * 2014-02-27 2014-06-25 浙江大学 Device and method for measuring viscoelasticity of material by using ultrasound water printing pressure
US8896838B2 (en) * 2010-03-05 2014-11-25 The General Hospital Corporation Systems, methods and computer-accessible medium which provide microscopic images of at least one anatomical structure at a particular resolution
CN104490365A (en) * 2015-01-12 2015-04-08 南京大学 Method for rebuilding broad-band photo-acoustic image by using narrow-band sensor
JP2015161505A (en) * 2014-02-26 2015-09-07 株式会社トプコン Calibration method of interferometer and interferometer using calibration method
CN105203461A (en) * 2015-09-17 2015-12-30 天津大学 Method for detecting piezoelectric material surface layer Young modulus through laser acoustic surface waves
CN105572048A (en) * 2015-12-16 2016-05-11 华南师范大学 Opto-acoustic quantitative elasticity and viscosity double-parameter imaging method and device
CN105572049A (en) * 2015-12-16 2016-05-11 华南师范大学 Opto-acoustic quantitative elasticity imaging method and device
US20160157828A1 (en) * 2014-06-05 2016-06-09 Chikayoshi Sumi Beamforming method, measurement and imaging instruments, and communication instruments
CN105769128A (en) * 2016-05-06 2016-07-20 华南师范大学 Integrated opto-acoustic, ultrasonic and opto-acoustic elastic endoscopic imaging device and method thereof

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6041020A (en) * 1997-04-21 2000-03-21 University Of Delaware Gas-coupled laser acoustic detection
EP1998161A1 (en) * 2007-05-29 2008-12-03 The Technical University of Denmark (DTU) Acoustic resonator cell for spectroscopic analysis of a compound in a fluid
US8896838B2 (en) * 2010-03-05 2014-11-25 The General Hospital Corporation Systems, methods and computer-accessible medium which provide microscopic images of at least one anatomical structure at a particular resolution
CN102226752A (en) * 2011-04-08 2011-10-26 常熟舒茨电子科技发展有限公司 Method for detecting concentration of sulfur dioxide by using photoacoustic spectrometry
US20140073907A1 (en) * 2012-09-12 2014-03-13 Convergent Life Sciences, Inc. System and method for image guided medical procedures
CN103324786A (en) * 2013-06-03 2013-09-25 北京工业大学 Two-dimensional solid-phase phonon crystal XY module cord gap optimization method based on genetic algorithm
JP2015161505A (en) * 2014-02-26 2015-09-07 株式会社トプコン Calibration method of interferometer and interferometer using calibration method
CN103884774A (en) * 2014-02-27 2014-06-25 浙江大学 Device and method for measuring viscoelasticity of material by using ultrasound water printing pressure
US20160157828A1 (en) * 2014-06-05 2016-06-09 Chikayoshi Sumi Beamforming method, measurement and imaging instruments, and communication instruments
CN104490365A (en) * 2015-01-12 2015-04-08 南京大学 Method for rebuilding broad-band photo-acoustic image by using narrow-band sensor
CN105203461A (en) * 2015-09-17 2015-12-30 天津大学 Method for detecting piezoelectric material surface layer Young modulus through laser acoustic surface waves
CN105572048A (en) * 2015-12-16 2016-05-11 华南师范大学 Opto-acoustic quantitative elasticity and viscosity double-parameter imaging method and device
CN105572049A (en) * 2015-12-16 2016-05-11 华南师范大学 Opto-acoustic quantitative elasticity imaging method and device
CN105769128A (en) * 2016-05-06 2016-07-20 华南师范大学 Integrated opto-acoustic, ultrasonic and opto-acoustic elastic endoscopic imaging device and method thereof

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
IRVING ITZKAN.ETAL: "The thermoelastic basis of short pulsed laser ablation of biological tissue", 《PROC NATL.ACAD.SCI》 *
SHAOHUA WANG.ETAL: "Quantitative detection of stochastic microstructure in turbid media by photoacoustic spectral matching", 《APPLIED PHYSICS LETTERS》 *
YUE ZHAO.ETAL: "Simultaneous optical absorption and viscoelasticity imaging based on photoacoustic lock-in measurement", 《OPTICS LETTERS》 *
王敬时等: "利用激光超声技术研究表面微裂纹缺陷材料的低通滤波效应", 《物理学报》 *
程曦等: "利用激光超声研究功能梯度材料中声表面波的传播特性", 《声学学报》 *

Similar Documents

Publication Publication Date Title
Lemistre et al. Structural health monitoring system based on diffracted Lamb wave analysis by multiresolution processing
Abbate et al. Signal detection and noise suppression using a wavelet transform signal processor: application to ultrasonic flaw detection
Gennisson et al. Assessment of elastic parameters of human skin using dynamic elastography
CN1774210B (en) Elastography device and method for determining and imaging of mechanical and elastic parameters of an examination object
Greenleaf et al. Detection of tissue harmonic motion induced by ultrasonic radiation force using pulse-echo ultrasound and Kalman filter
Rose et al. Determination of flaw characteristics from ultrasonic scattering data
US20100191109A1 (en) Biological information processing apparatus and biological information processing method
CA2647283C (en) A method and a device for imaging a visco-elastic medium
Gore et al. Ultrasonic backscattering from human tissue: a realistic model
Legendre et al. Wavelet-transform-based method of analysis for Lamb-wave ultrasonic NDE signals
Yu et al. Experimental ultrasound characterization of red blood cell aggregation using the structure factor size estimator
Lazaro et al. Influence of thresholding procedures in ultrasonic grain noise reduction using wavelets
US20110239766A1 (en) Photoacoustic imaging apparatus and photoacoustic imaging method
CN101595396A (en) High-resolution ultrasound displacement measurement apparatus and method
Hachiya et al. Determination of sound speed in biological tissues based on frequency analysis of pulse response
CA2309916C (en) Ultrasonic system for grading meat
CN102066921A (en) Method for the non-destructive testing of a test object by way of ultrasound and corresponding device
Ta et al. Analysis of frequency dependence of ultrasonic backscatter coefficient in cancellous bone
Socco et al. Feasibility of ultrasonic tomography for nondestructive testing of decay on living trees
Labyed et al. A theoretical comparison of attenuation measurement techniques from backscattered ultrasound echoes
Wang et al. Pear dynamic characteristics and firmness detection
CN102283679B (en) Ultrasonic imaging system for elasticity measurement and method for measuring elasticity of biological tissue
Bigelow et al. In vivo ultrasonic attenuation slope estimates for detecting cervical ripening in rats: Preliminary results
JP5850633B2 (en) Object information acquiring apparatus
US20020035872A1 (en) Apparatus and method for evaluating the physical properties of a sample using ultrasonics

Legal Events

Date Code Title Description
C06 Publication
SE01