CN111289620A - 各向异性材料的弹性常数检测方法及系统 - Google Patents

各向异性材料的弹性常数检测方法及系统 Download PDF

Info

Publication number
CN111289620A
CN111289620A CN201811493892.8A CN201811493892A CN111289620A CN 111289620 A CN111289620 A CN 111289620A CN 201811493892 A CN201811493892 A CN 201811493892A CN 111289620 A CN111289620 A CN 111289620A
Authority
CN
China
Prior art keywords
anisotropic material
ultrasonic
probe
line
focusing
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.)
Granted
Application number
CN201811493892.8A
Other languages
English (en)
Other versions
CN111289620B (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.)
AECC Commercial Aircraft Engine Co Ltd
Original Assignee
AECC Commercial Aircraft Engine Co Ltd
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 AECC Commercial Aircraft Engine Co Ltd filed Critical AECC Commercial Aircraft Engine Co Ltd
Priority to CN201811493892.8A priority Critical patent/CN111289620B/zh
Publication of CN111289620A publication Critical patent/CN111289620A/zh
Application granted granted Critical
Publication of CN111289620B publication Critical patent/CN111289620B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/14Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object using acoustic emission techniques

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (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 Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本公开涉及一种各向异性材料的弹性常数检测方法及系统,方法包括:将各向异性材料试样设置在测试台上,并在测试时驱动超声显微镜的线聚焦探头改变相对于所述各向异性材料的距离,以调整所述线聚焦探头在所述各向异性材料内部的超声聚焦位置;记录在所述各向异性材料内部的不同超声聚焦位置下的各个时间间隔下的超声反射信号;对所述超声反射信号进行傅里叶变换,并获得指定频率下所述各向异性材料的声学特征曲线的实验值;对根据所述各向异性材料的弹性张量初始值计算出的所述各向异性材料的声学特征曲线的理论值与所述声学特征曲线的实验值进行迭代拟合,以获得所述各向异性材料的实际弹性张量。本公开实施例能够提高检测效率。

Description

各向异性材料的弹性常数检测方法及系统
技术领域
本公开涉及检测技术领域,尤其涉及一种各向异性材料的弹性常数检测方法及系统。
背景技术
复合材料是由两种或两种以上不同性质的材料,通过物理或化学的方法,在宏观(微观)上组成具有新性能的材料。在复合材料各向异性的表征方面,主要是基于破坏性的力学实验测量方法,制作不同方向的拉伸、压缩、弯曲、剪切等试样,通过拉伸、压缩、弯曲、剪切等力学试验来测定复合材料的弹性参数,以实现复合材料各向异性的测量、表征。这种破坏性的力学测试方法比较直接,但是需要对试样进行破坏性的实验,而且只能获得一些主方向上的弹性参数,很难获得完整的弹性张量参数。
由于复合材料的声学特性也随着方向的变化而变化,当声波以相同的入射角在不同方向上入射到复合材料表面时,反射声波的强度会有差异。这种差异可以通过材料的刚度矩阵或弹性张量,经过弹性力学和声学的物理建模和公式推导建立复合材料声波入射方向和反射信号之间的数学关系,从而通过逆运算实现材料刚度矩阵的表征,进而实现对复合材料各向异性的表征。线聚焦超声检测技术恰好可以在一定方向上对复合材料进行检测,通过上述过程实现对复合材料各向异性的表征。
聚焦超声检测技术是利用高频线聚焦探头发射和接收超声波对物体进行无损检测的超声检测技术。这种技术可对物体内部及亚表面结构进行高精度、高灵敏度的检测,并可观察从表层到数十毫米深度范围内存在的微米到亚微米尺度结构。聚焦超声检测技术具有高灵敏度、高分辨率和图像直观等特点,目前已广泛应用于电子工业、医学、材料科学等领域。
聚焦超声检测一般采用水作耦合剂,使用水浸聚焦超声探头进行检测。聚焦超声探头是一种能实现声束聚焦的超声波换能器,根据聚焦方式分为点聚焦和线聚焦超声波换能器。采用线聚焦探头的超声检测技术称为线聚焦超声检测技术。不同于点聚焦探头,线聚焦探头能在特定的方向上发射和接收超声检测信号,进而实现对材料内部不同方向上的声学特性进行表征,在复合材料各向异性的检测与表征方面具有独特的优势。
复合材料的各向异性可以通过4阶弹性张量表示为:
C=Cijkl(i,j,k,l=1,2,3,4,5,6) (1)
根据弹性力学,其运动方程可以表示为:
Figure BDA0001896360830000021
其中,ρ表示密度,u表示质点位移。
如果做如下代换:
(11)~1;(22)~2;(33)~3;(23)=(32)~4;(13)=(31)~5;(12)=(21)~6; (3)
4阶弹性张量Cijkl可以表示为:
Figure BDA0001896360830000022
对于不同的复合材料结构,其弹性张量Cαβ是不同的,通过测定Cαβ中各个参数,就可以实现复合材料各向异性的表征。
使用超声检测方法测定复合材料各向异性,主要是通过采用2个超声探头(一个发射探头,一个接收探头)在不同方向上以不同的入射角发射和接收超声波,再通过不同方向、角度的超声波反射信号来计算不同方向的反射系数,绘制频谱曲线,再通过建立频谱曲线和弹性张量Cαβ(见式(4))的关系,通过逆运算计算出计算Cαβ,从而实现对复合材料各向异性的表征。
然而,常规的采用双探头测定和材料反射系数实现复合材料各向异性表征的方法,需要在某个方向上不同的角度上进行重复测量反射系数,绘制不同方向上的频谱曲线,效率比较低下。另一方面,由于目前针对复合材料各向异性的线聚焦探头的超声检测技术,目前研究较少,缺乏一套系统的有效的方法和数学模型建立超声检测信号和复合材料弹性参数(弹性张量)的关系式,因此计算过程非常繁琐,不利于实际的检测应用。
发明内容
有鉴于此,本公开实施例提供一种各向异性材料的弹性常数检测方法及系统,能够提高检测效率。
在本公开的一个方面,提供一种各向异性材料的弹性常数检测方法,包括:
将各向异性材料试样设置在测试台上,并在测试时驱动超声显微镜的线聚焦探头改变相对于所述各向异性材料的距离,以调整所述线聚焦探头在所述各向异性材料内部的超声聚焦位置;
记录在所述各向异性材料内部的不同超声聚焦位置下的各个时间间隔下的超声反射信号;
对所述超声反射信号进行傅里叶变换,并获得指定频率下所述各向异性材料的声学特征曲线的实验值;
对根据所述各向异性材料的弹性张量初始值计算出的所述各向异性材料的声学特征曲线的理论值与所述声学特征曲线的实验值进行迭代拟合,以获得所述各向异性材料的实际弹性张量。
在一些实施例中,将各向异性材料试样设置在测试台的操作包括:
将所述各向异性材料试样放入所述超声显微镜的工作槽内,所述工作槽采用干净且无气泡的水作耦合剂。
在一些实施例中,在测试之前还包括:
调整所述线聚焦探头,以使所述线聚焦探头的轴向与所述各向异性材料试样的反射表面垂直,并使所述线聚焦探头的聚焦线与所述反射表面平行。
在一些实施例中,驱动线聚焦探头改变相对于所述各向异性材料的距离的操作包括:
调整所述线聚焦探头,以确定所述线聚焦探头的焦距;
在所述线聚焦探头与所述各向异性材料的距离为所述焦距时,将所述线聚焦探头在垂直于所述反射表面的方向上的位置设定为零点;
沿垂直于所述反射表面的第一方向,驱动所述线聚焦探头从零点移动到预设检测距离,所述预设检测距离小于所述焦距;
驱动所述线聚焦探头以预设步长沿第二方向运动,并运动到两倍的预设检测距离,所述第二方向为所述第一方向的反方向,以调整所述线聚焦探头在各个时间间隔下在所述各向异性材料内部的超声聚焦位置。
在一些实施例中,指定频率的确定操作包括:
在所述线聚焦探头运动到两倍的预设检测距离后,获得当前时间间隔下的超声聚焦位置的超声反射信号;
对所述当前时间间隔下的超声聚焦位置的超声反射信号进行傅里叶变换,获得频率曲线,并从曲线中确定最大值点所对应的频率作为指定频率。
在一些实施例中,对所述超声反射信号进行傅里叶变换的操作包括:
对不同超声聚焦位置的各个时间间隔下的超声反射信号进行二维傅里叶变换,并获取变换结果中指定频率下的声学特征曲线的实验值。
在一些实施例中,根据所述各向异性材料的弹性张量初始值计算所述各向异性材料的声学特征曲线的理论值的操作包括:
根据所述各向异性材料的弹性张量初始值计算声反射系数;
根据计算出的所述声反射系数计算所述各向异性材料的声学特征曲线的理论值。
在本公开的一个方面,提供一种各向异性材料的弹性常数检测系统,包括:
超声显微镜,包括超声波收发装置、线聚焦探头、工作槽和机械扫描装置,所述超声波收发装置通过所述线聚焦探头发出超声波,并接收超声反射信号,所述工作槽用于盛放待检测的各向异性材料试样,所述机械扫描装置用于驱动所述线聚焦探头相对于所述工作槽运动;和
上位机,与所述超声显微镜连接,用于记录在所述各向异性材料内部的不同超声聚焦位置下的各个时间间隔下的超声反射信号,对所述超声反射信号进行傅里叶变换,并获得指定频率下所述各向异性材料的声学特征曲线的实验值,以及对根据所述各向异性材料的弹性张量初始值计算出的所述各向异性材料的声学特征曲线的理论值与所述声学特征曲线的实验值进行迭代拟合,以获得所述各向异性材料的实际弹性张量。
在一些实施例中,在所述工作槽内设有干净且无气泡的水作耦合剂。
在一些实施例中,所述机械扫描装置被配置为调整所述线聚焦探头,以使所述线聚焦探头的轴向与所述各向异性材料试样的反射表面垂直,并使所述线聚焦探头的聚焦线与所述反射表面平行,然后调整所述线聚焦探头,以确定所述线聚焦探头的焦距,在所述线聚焦探头与所述各向异性材料的距离为所述焦距时,将所述线聚焦探头在垂直于所述反射表面的方向上的位置设定为零点,沿垂直于所述反射表面的第一方向,驱动所述线聚焦探头从零点移动到预设检测距离,所述预设检测距离小于所述焦距,驱动所述线聚焦探头以预设步长沿第二方向运动,并运动到两倍的预设检测距离,所述第二方向为所述第一方向的反方向,以调整所述线聚焦探头在各个时间间隔下在所述各向异性材料内部的超声聚焦位置。
因此,根据本公开实施例,通过对各向异性材料试样的超声测试获得指定频率下各向异性材料的声学特征曲线的实验值,并将实验值与根据各向异性材料的弹性张量初始值计算出的声学特征曲线的理论值进行拟合迭代,可获得各向异型材料的实际弹性张量。这种弹性常数检测方法既不需要进行破坏性试验,也无需进行多角度重复测试,因此在检测效率方面较高。
附图说明
构成说明书的一部分的附图描述了本公开的实施例,并且连同说明书一起用于解释本公开的原理。
参照附图,根据下面的详细描述,可以更加清楚地理解本公开,其中:
图1是根据本公开各向异性材料的弹性常数检测系统的一些实施例的结构示意图;
图2是根据本公开各向异性材料的弹性常数检测方法的一些实施例中涉及的各种参数的示意图;
图3是根据本公开各向异性材料的弹性常数检测方法的一些实施例的流程示意图;
图4是根据本公开各向异性材料的弹性常数检测方法的一些实施例中在某聚焦深度下的超声反射信号图;
图5是根据本公开各向异性材料的弹性常数检测方法的另一些实施例中驱动线聚焦探头的流程示意图;
图6是根据本公开各向异性材料的弹性常数检测方法的另一些实施例中确定指定频率的流程示意图;
图7是根据本公开各向异性材料的弹性常数检测方法的另一些实施例中计算声学特征曲线的理论值的流程示意图;
图8是根据本公开各向异性材料的弹性常数检测方法的第一测试实例得到的声反射系数的曲线图;
图9是根据本公开各向异性材料的弹性常数检测方法的第一测试实例的声学特征曲线的实验值和理论值的拟合图;
图10是根据本公开各向异性材料的弹性常数检测方法的第二测试实例得到的声反射系数的曲线图;
图11是根据本公开各向异性材料的弹性常数检测方法的第二测试实例的声学特征曲线的实验值和理论值的拟合图。
应当明白,附图中所示出的各个部分的尺寸并不是按照实际的比例关系绘制的。此外,相同或类似的参考标号表示相同或类似的构件。
具体实施方式
现在将参照附图来详细描述本公开的各种示例性实施例。对示例性实施例的描述仅仅是说明性的,决不作为对本公开及其应用或使用的任何限制。本公开可以以许多不同的形式实现,不限于这里所述的实施例。提供这些实施例是为了使本公开透彻且完整,并且向本领域技术人员充分表达本公开的范围。应注意到:除非另外具体说明,否则在这些实施例中阐述的部件和步骤的相对布置、数字表达式和数值应被解释为仅仅是示例性的,而不是作为限制。
本公开中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性,而只是用来区分不同的部分。“包括”或者“包含”等类似的词语意指在该词前的要素涵盖在该词后列举的要素,并不排除也涵盖其他要素的可能。“上”、“下”、“左”、“右”等仅用于表示相对位置关系,当被描述对象的绝对位置改变后,则该相对位置关系也可能相应地改变。
在本公开中,当描述到特定器件位于第一器件和第二器件之间时,在该特定器件与第一器件或第二器件之间可以存在居间器件,也可以不存在居间器件。当描述到特定器件连接其它器件时,该特定器件可以与所述其它器件直接连接而不具有居间器件,也可以不与所述其它器件直接连接而具有居间器件。
本公开使用的所有术语(包括技术术语或者科学术语)与本公开所属领域的普通技术人员理解的含义相同,除非另外特别定义。还应当理解,在诸如通用字典中定义的术语应当被解释为具有与它们在相关技术的上下文中的含义相一致的含义,而不应用理想化或极度形式化的意义来解释,除非这里明确地这样定义。
对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论,但在适当情况下,所述技术、方法和设备应当被视为说明书的一部分。
如图1所示,是根据本公开各向异性材料的弹性常数检测系统的一些实施例的结构示意图。参考图1,在一些实施例中,各向异性材料的弹性常数检测系统包括超声显微镜和上位机30。超声显微镜可,包括超声波收发装置24、线聚焦探头22、工作槽21和机械扫描装置23。所述超声波收发装置24可通过所述线聚焦探头22发出超声波,并接收超声反射信号。所述工作槽21用于盛放待检测的各向异性材料试样40。所述机械扫描装置23用于驱动所述线聚焦探头22相对于所述工作槽21运动。在一些实施例中,超声显微镜还可以包括示波器25和数据采集装置26。
上位机30与所述超声显微镜连接,可用于记录在所述各向异性材料内部的不同超声聚焦位置下的各个时间间隔下的超声反射信号。上位机30还可以对所述超声反射信号进行傅里叶变换,并获得指定频率下所述各向异性材料的声学特征曲线的实验值,并对根据所述各向异性材料的弹性张量初始值计算出的所述各向异性材料的声学特征曲线的理论值与所述声学特征曲线的实验值进行迭代拟合,以获得所述各向异性材料的实际弹性张量。
参考图2,为了方便计算,以各向异性材料试样为参考对坐标系以及部分参数进行定义。根据线聚焦探头相对于各向异性材料试样的对应位置,可将平行于各向异性材料试样表面的平面定义为X-Y平面,其中选取正交的两个方向分别为X轴方向和Y轴方向。将垂直于X-Y平面且离开各向异性材料试样的方向定义为Z轴方向。根据超声波的入射方向可定义θ为入射角,超声波的入射方向在X-Y平面的投影Yφ与Y轴方向的夹角为φ。
为了计算根据所述各向异性材料的弹性张量初始值计算出的所述各向异性材料的声学特征曲线的理论值,需要建立声反射系数R(θ)与各向异性材料的弹性张量之间的关系表达式。
考虑一个包含N层的各向异性复合材料体系,由Stroh方程得到:
Figure BDA0001896360830000091
其中,ξ为状态向量,表示为:
Figure BDA0001896360830000092
U和T分别为位移矢量和应力矢量:
U=[ux,uy,uz]T,T=[σxzyzzz]T (8)
A(z)可以表示为:
Figure BDA0001896360830000093
式中:
Figure BDA0001896360830000094
此时根据弹性力学的传递矩阵关系,有:
B(z1,z0)=eiAh (11)
h=z1-z0,是该层的厚度;
刚度矩阵和传递矩阵的关系为:
Figure BDA0001896360830000095
通过递归刚度矩阵方法,得到此多层各向异性复合材料的总刚度矩阵K,其递归公式如下:
Figure BDA0001896360830000096
由式(13)得到,此多层各向异性材料总柔度矩阵:
Figure BDA0001896360830000097
当声波从水中入射到该多层各向异性材料表面时,其声反射系数的计算公式为:
Figure BDA0001896360830000101
其中,Λ=cosθ/(iωρfVf),ω为角频率,θ为入射角,ρf为水的密度,Vf为水中声速。
接下来,建立线聚焦探头的响应函数V(z)(即声学特征曲线)与声反射系数R(θ)的关系。对于线聚焦探头,式(5)可以表示为:
Figure BDA0001896360830000102
其中,kf为波数,kx是波数在X轴方向上的分量,β为线聚焦探头晶片弧度角,ρ是材料密度,z是线聚焦探头沿Z轴的位移量。
通过上述推导过程可获得声反射系数R(θ)与各向异性材料的弹性张量的关系式(15),以及声反射系数R(θ)与声学特征曲线V(z)的关系式(16),从而可实现各向异性材料的声学特征曲线的理论值Vcal(z)的计算。
如图3所示,是根据本公开各向异性材料的弹性常数检测方法的一些实施例的流程示意图。参考图3,在一些实施例中,各向异性材料的弹性常数检测方法包括:
步骤100、将各向异性材料试样设置在测试台上;
步骤200、在测试时驱动超声显微镜的线聚焦探头改变相对于所述各向异性材料的距离,以调整所述线聚焦探头在所述各向异性材料内部的超声聚焦位置;
步骤300、记录在所述各向异性材料内部的不同超声聚焦位置下的各个时间间隔下的超声反射信号;
步骤400、对所述超声反射信号进行傅里叶变换,并获得指定频率下所述各向异性材料的声学特征曲线的实验值;
步骤500、对根据所述各向异性材料的弹性张量初始值计算出的所述各向异性材料的声学特征曲线的理论值与所述声学特征曲线的实验值进行迭代拟合,以获得所述各向异性材料的实际弹性张量。
在本实施例中,各向异性材料试样可以为复合材料,也可以为非复合材料。在步骤100中,可将各向异性材料试样放入所述超声显微镜的工作槽内,并且在所述工作槽采用干净且无气泡的水作耦合剂。
在步骤200的测试操作之前,可先调整所述线聚焦探头,使其轴向与所述各向异性材料试样的反射表面(即图2的X-Y平面)垂直,并使所述线聚焦探头的聚焦线与所述反射表面平行。该聚焦线与Y轴方向呈角度φ。
如图5所示,是根据本公开各向异性材料的弹性常数检测方法的另一些实施例中驱动线聚焦探头的流程示意图。参考图5,在一些实施例中,步骤200中驱动超声显微镜的线聚焦探头改变相对于所述各向异性材料的距离的操作可包括:
步骤210、调整所述线聚焦探头,以确定所述线聚焦探头的焦距;
步骤220、在所述线聚焦探头与所述各向异性材料的距离为所述焦距时,将所述线聚焦探头在垂直于所述反射表面的方向上的位置设定为零点;
步骤230、沿垂直于所述反射表面的第一方向,驱动所述线聚焦探头从零点移动到预设检测距离,所述预设检测距离小于所述焦距;
步骤240、驱动所述线聚焦探头以预设步长沿第二方向运动,并运动到两倍的预设检测距离,所述第二方向为所述第一方向的反方向,以调整所述线聚焦探头在各个时间间隔下在所述各向异性材料内部的超声聚焦位置。
在步骤210中,可通过机械扫描装置在Z轴方向上调整线聚焦探头的位置,以使测得的超声反射信号最大。此时线聚焦探头与各向异性材料试样的表面的距离即为线聚焦探头的焦距D。在步骤220中确定线聚焦探头在Z轴方向的零点后,可在步骤230中通过机械扫描装置将探头向Z轴方向的反方向移动预设检测距离z0,z0的大小根据聚焦探头的焦距而定(z0<D)。
在步骤240中,可控制机械扫描装置以较小的步长量在Z轴方向上移动线聚焦探头。每移动一个步长,可稳定一个时间间隔T,待信号稳定后可自动记录该位置的超声反射信号s(z,t)。例如图4所示的超声聚焦位置在z=-2.0mm时所停留的时间间隔下的超声反射信号图,在图4中,横轴为停留时间间隔中的时间t,纵轴为超声反射信号s(z,t)的信号幅度。通过对每个超声聚焦位置的记录,可以获得一组不同位置的超声反射信号S(z,t),(-z0≤z≤z0)。
如图6所示,是根据本公开各向异性材料的弹性常数检测方法的另一些实施例中确定指定频率的流程示意图。在本实施例中,当需要确定指定频率时,可执行步骤410和步骤420。在步骤410中,可在所述线聚焦探头运动到位置z后,获得当前时间间隔下的超声聚焦位置的超声反射信号。具体来说,在线聚焦探头以步长d从-z0运动到z0时,可对超声聚焦位置为z时的超声反射信号s(z,t)进行采集。然后可在步骤420中对所述当前时间间隔下的超声聚焦位置的超声反射信号进行傅里叶变换,获得频率曲线,即F[s(z,f)],并从曲线中确定最大值点所对应的频率作为指定频率,也即线聚焦探头的中心频率。
参考图6,在一些实施例中,步骤400还可以包括步骤430,即对不同超声聚焦位置的各个时间间隔下的超声反射信号进行二维傅里叶变换。然后从变换结果中获取前述指定频率下的声学特征曲线的实验值。
在获得声学特征曲线的实验值后,可进行步骤500中声学特征曲线的理论值与实验值的迭代拟合操作。如图7所示,是根据本公开各向异性材料的弹性常数检测方法的另一些实施例中计算声学特征曲线的理论值的流程示意图。在图7中,根据所述各向异性材料的弹性张量初始值计算所述各向异性材料的声学特征曲线的理论值的操作可包括:
步骤610、根据所述各向异性材料的弹性张量初始值计算声反射系数;
步骤620、根据计算出的所述声反射系数计算所述各向异性材料的声学特征曲线的理论值。
在步骤610中,根据刚度矩阵与柔度矩阵的互逆关系,可将各向异性材料的弹性张量初始值代入前述公式(15),来计算出声反射系数R(θ)。在获得声反射系数R(θ)之后,在步骤620则将声反射系数R(θ)代入前述公式(16),计算出该各向异性材料的声学特征曲线的理论值。
上述声学特征曲线的理论值的计算可以在获得实验值之前、之后或者同时。在获得了声学特征曲线的理论值和实验值之后,则可对这两条声学特征曲线进行迭代拟合,最终获得该各向异性材料的弹性张量,以实现各向异性材料的各向异性的表征。可用的迭代拟合公式例如
Figure BDA0001896360830000131
zi为探头在第i个采样位置的坐标(i≤N),N为采样点数。
接下来将通过两个测试实例来具体说明本公开的各向异性材料的弹性常数检测方法的有效性。
测试实例1
选取一厚度为200mm的[0°/45°/90°/-45]铺层复合材料作为待测的各向异性材料试样。该铺层复合材料的弹性张量(φ=0°)如下表:
Figure BDA0001896360830000132
将该各向异性材料试样放入超声显微镜设备的工作槽内,工作槽采用干净、无气泡的水做耦合剂。调整机械扫描装置使线聚焦探头的方向,使得线聚焦探头的轴向方向与试样反射表面保持垂直,并使线聚焦探头的聚焦线与X轴方向平行(此时φ=0°)。
采用机械扫描装置在Z轴方向上调整线聚焦探头的位置,使得测得的超声反射信号最大,此时探头与标样表面的距离即为探头的焦距D。将此时线聚焦探头Z轴方向的位置设置为零,通过机械扫描装置将探头向Z轴方向的反方向移动预设检测距离z0,z的大小根据线聚焦探头的焦距而定(z0<D)。
控制机械扫查装置以预设的步长量d(例如d=0.01mm)在Z轴方向向上移动线聚焦探头,每移动一个步长量,则稳定一个时间间隔T,待信号稳定后自动记录该位置的超声反射信号s(z,t)。如此依次记录,可以获得一组不同位置的超声反射信号S(z,t),(-z0≤z≤z0)。
对超声反射信号s(z,t)进行傅里叶变换,得到F[s(z,f)]曲线,该曲线上最大值点对应的频率fc即是该探头的中心频率。对测得的S(z,t)进行二维傅里叶变换获得一组数据F[S(z,f)],并选择该频率fc下的数据,即得到该频率下该各向异性材料试样的声学特征曲线的实验值VExp(z),参见图9。
通过上表中记录的参数,将其代入公式(15)计算其声反射系数R(θ),可得到图8所示的声反射系数R(θ),再通过公式(16),可计算获得该各向异性材料试样的声学特征曲线的理论值VCal(z),参见图9。
对图9中的上述实验值VExp(z)和理论值VCal(z)进行迭代拟合,可获得该方向上复合材料的弹性张量:C11=140.5GPa;C22=15.8GPa;C33=15.8GPa;C12=7.45GPa;C13=7.45GPa;C23=8.2GPa;C44=3.8GPa;C55=6.9GPa;C66=6.9GPa。与上表中的真实值相比,最大误差仅为1.4%。
测试实例2
选取一个铺层复合材料作为待测的各向异性材料试样。该铺层复合材料的弹性张量(φ=0°)如下表:
Figure BDA0001896360830000141
将该各向异性材料试样放入超声显微镜设备的工作槽内,工作槽采用干净、无气泡的水做耦合剂。调整机械扫描装置使线聚焦探头的方向,使得线聚焦探头的轴向方向与试样反射表面保持垂直,并使线聚焦探头的聚焦线与X轴方向平行(此时φ=0°)。
采用机械扫描装置在Z轴方向上调整线聚焦探头的位置,使得测得的超声反射信号最大,此时探头与标样表面的距离即为探头的焦距D。将此时线聚焦探头Z轴方向的位置设置为零,通过机械扫描装置将探头向Z轴方向的反方向移动预设检测距离z0,z0的大小根据线聚焦探头的焦距而定(z0<D)。
控制机械扫查装置以预设的步长量d(例如d=0.01mm)在Z轴方向向上移动线聚焦探头,每移动一个步长量,则稳定一个时间间隔T,待信号稳定后自动记录该位置的超声反射信号s(z,t)。如此依次记录,可以获得一组不同位置的超声反射信号S(z,t),(-z0≤z≤z0)。
对超声反射信号s(z,t)进行傅里叶变换,得到F[s(z,f)]曲线,该曲线上最大值点对应的频率fc即是该探头的中心频率。对测得的S(z,t)进行二维傅里叶变换获得一组数据F[S(z,f)],并选择该频率fc下的数据,即得到该频率下该各向异性材料试样的声学特征曲线的实验值VExp(z),参见图11。
通过上表中记录的参数,将其代入公式(15)计算其声反射系数R(θ),可得到图10所示的声反射系数R(θ),再通过公式(16),可计算获得该各向异性材料试样的声学特征曲线的理论值VCal(z),参见图11。
对图11中的上述实验值VExp(z)和理论值VCal(z)进行迭代拟合,可获得该方向上复合材料的弹性张量:C11=12.0GPa;C22=12.0GPa;C33=130.1GPa;C12=5.5GPa;C13=5.9GPa;C23=6.9GPa;C44=6.15GPa;C55=6.21GPa;C66=3.2GPa。与上表中的真实值相比,最大误差仅为3.6%。
至此,已经详细描述了本公开的各实施例。为了避免遮蔽本公开的构思,没有描述本领域所公知的一些细节。本领域技术人员根据上面的描述,完全可以明白如何实施这里公开的技术方案。
虽然已经通过示例对本公开的一些特定实施例进行了详细说明,但是本领域的技术人员应该理解,以上示例仅是为了进行说明,而不是为了限制本公开的范围。本领域的技术人员应该理解,可在不脱离本公开的范围和精神的情况下,对以上实施例进行修改或者对部分技术特征进行等同替换。本公开的范围由所附权利要求来限定。

Claims (10)

1.一种各向异性材料的弹性常数检测方法,包括:
将各向异性材料试样设置在测试台上,并在测试时驱动超声显微镜的线聚焦探头改变相对于所述各向异性材料的距离,以调整所述线聚焦探头在所述各向异性材料内部的超声聚焦位置;
记录在所述各向异性材料内部的不同超声聚焦位置下的各个时间间隔下的超声反射信号;
对所述超声反射信号进行傅里叶变换,并获得指定频率下所述各向异性材料的声学特征曲线的实验值;
对根据所述各向异性材料的弹性张量初始值计算出的所述各向异性材料的声学特征曲线的理论值与所述声学特征曲线的实验值进行迭代拟合,以获得所述各向异性材料的实际弹性张量。
2.根据权利要求1所述的弹性常数检测方法,其中,将各向异性材料试样设置在测试台的操作包括:
将所述各向异性材料试样放入所述超声显微镜的工作槽内,所述工作槽采用干净且无气泡的水作耦合剂。
3.根据权利要求1所述的弹性常数检测方法,其中,在测试之前还包括:
调整所述线聚焦探头,以使所述线聚焦探头的轴向与所述各向异性材料试样的反射表面垂直,并使所述线聚焦探头的聚焦线与所述反射表面平行。
4.根据权利要求3所述的弹性常数检测方法,其中,驱动线聚焦探头改变相对于所述各向异性材料的距离的操作包括:
调整所述线聚焦探头,以确定所述线聚焦探头的焦距;
在所述线聚焦探头与所述各向异性材料的距离为所述焦距时,将所述线聚焦探头在垂直于所述反射表面的方向上的位置设定为零点;
沿垂直于所述反射表面的第一方向,驱动所述线聚焦探头从零点移动到预设检测距离,所述预设检测距离小于所述焦距;
驱动所述线聚焦探头以预设步长沿第二方向运动,并运动到两倍的预设检测距离,所述第二方向为所述第一方向的反方向,以调整所述线聚焦探头在各个时间间隔下在所述各向异性材料内部的超声聚焦位置。
5.根据权利要求1所述的弹性常数检测方法,其中,指定频率的确定操作包括:
在所述线聚焦探头以预设步长运动时,获得当前时间间隔下的超声聚焦位置的超声反射信号;
对所述当前时间间隔下的超声聚焦位置的超声反射信号进行傅里叶变换,获得频率曲线,并从曲线中确定最大值点所对应的频率作为指定频率。
6.根据权利要求1所述的弹性常数检测方法,其中,对所述超声反射信号进行傅里叶变换的操作包括:
对不同超声聚焦位置的各个时间间隔下的超声反射信号进行二维傅里叶变换,并获取变换结果中指定频率下的声学特征曲线的实验值。
7.根据权利要求1所述的弹性常数检测方法,其中,根据所述各向异性材料的弹性张量初始值计算所述各向异性材料的声学特征曲线的理论值的操作包括:
根据所述各向异性材料的弹性张量初始值计算声反射系数;
根据计算出的所述声反射系数计算所述各向异性材料的声学特征曲线的理论值。
8.一种各向异性材料的弹性常数检测系统,包括:
超声显微镜,包括超声波收发装置、线聚焦探头、工作槽和机械扫描装置,所述超声波收发装置通过所述线聚焦探头发出超声波,并接收超声反射信号,所述工作槽用于盛放待检测的各向异性材料试样,所述机械扫描装置用于驱动所述线聚焦探头相对于所述工作槽运动;和
上位机,与所述超声显微镜连接,用于记录在所述各向异性材料内部的不同超声聚焦位置下的各个时间间隔下的超声反射信号,对所述超声反射信号进行傅里叶变换,并获得指定频率下所述各向异性材料的声学特征曲线的实验值,以及对根据所述各向异性材料的弹性张量初始值计算出的所述各向异性材料的声学特征曲线的理论值与所述声学特征曲线的实验值进行迭代拟合,以获得所述各向异性材料的实际弹性张量。
9.根据权利要求8所述的弹性常数检测系统,其中,在所述工作槽内设有干净且无气泡的水作耦合剂。
10.根据权利要求8所述的弹性常数检测系统,其中,所述机械扫描装置被配置为调整所述线聚焦探头,以使所述线聚焦探头的轴向与所述各向异性材料试样的反射表面垂直,并使所述线聚焦探头的聚焦线与所述反射表面平行,然后调整所述线聚焦探头,以确定所述线聚焦探头的焦距,在所述线聚焦探头与所述各向异性材料的距离为所述焦距时,将所述线聚焦探头在垂直于所述反射表面的方向上的位置设定为零点,沿垂直于所述反射表面的第一方向,驱动所述线聚焦探头从零点移动到预设检测距离,所述预设检测距离小于所述焦距,驱动所述线聚焦探头以预设步长沿第二方向运动,并运动到两倍的预设检测距离,所述第二方向为所述第一方向的反方向,以调整所述线聚焦探头在各个时间间隔下在所述各向异性材料内部的超声聚焦位置。
CN201811493892.8A 2018-12-07 2018-12-07 各向异性材料的弹性常数检测方法及系统 Active CN111289620B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811493892.8A CN111289620B (zh) 2018-12-07 2018-12-07 各向异性材料的弹性常数检测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811493892.8A CN111289620B (zh) 2018-12-07 2018-12-07 各向异性材料的弹性常数检测方法及系统

Publications (2)

Publication Number Publication Date
CN111289620A true CN111289620A (zh) 2020-06-16
CN111289620B CN111289620B (zh) 2023-06-23

Family

ID=71025253

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811493892.8A Active CN111289620B (zh) 2018-12-07 2018-12-07 各向异性材料的弹性常数检测方法及系统

Country Status (1)

Country Link
CN (1) CN111289620B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113607814A (zh) * 2021-07-30 2021-11-05 广东工业大学 一种金属增材制造件弹性常数激光超声测量方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5211059A (en) * 1990-11-29 1993-05-18 Hitachi Construction Machinery Co., Ltd. Acoustic microscope system
US20070151340A1 (en) * 2005-12-30 2007-07-05 Industrial Technology Research Institute Nano-indentation ultrasonic detecting system and method thereof
CN105092595A (zh) * 2015-08-31 2015-11-25 哈尔滨工业大学(威海) 应用于钢轨探伤的光声弹性成像方法及装置
CN105486424A (zh) * 2014-09-17 2016-04-13 南京理工大学 燃烧室内壁瞬态温度场的超声波非侵入式测量方法
CN107490627A (zh) * 2016-06-12 2017-12-19 中国航发商用航空发动机有限责任公司 超声聚焦探头参数标定方法
CN107941907A (zh) * 2017-10-31 2018-04-20 武汉大学 一种基于有效超声背散射信号提取多晶材料的晶粒平均尺寸的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5211059A (en) * 1990-11-29 1993-05-18 Hitachi Construction Machinery Co., Ltd. Acoustic microscope system
US20070151340A1 (en) * 2005-12-30 2007-07-05 Industrial Technology Research Institute Nano-indentation ultrasonic detecting system and method thereof
CN105486424A (zh) * 2014-09-17 2016-04-13 南京理工大学 燃烧室内壁瞬态温度场的超声波非侵入式测量方法
CN105092595A (zh) * 2015-08-31 2015-11-25 哈尔滨工业大学(威海) 应用于钢轨探伤的光声弹性成像方法及装置
CN107490627A (zh) * 2016-06-12 2017-12-19 中国航发商用航空发动机有限责任公司 超声聚焦探头参数标定方法
CN107941907A (zh) * 2017-10-31 2018-04-20 武汉大学 一种基于有效超声背散射信号提取多晶材料的晶粒平均尺寸的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
谢刚等: "利用测井资料预测水平井压裂裂缝形态", 《测井技术》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113607814A (zh) * 2021-07-30 2021-11-05 广东工业大学 一种金属增材制造件弹性常数激光超声测量方法及系统

Also Published As

Publication number Publication date
CN111289620B (zh) 2023-06-23

Similar Documents

Publication Publication Date Title
US5549003A (en) Method and apparatus for visualization of internal stresses in solid non-transparent materials by ultrasonic techniques and ultrasonic computer tomography of stress
CA2619156C (en) Inspection system and associated method
CN107688051B (zh) 一种基于激光超声表面波的亚表面缺陷宽度的测量方法
Degertakin et al. Lamb wave excitation by Hertzian contacts with applications in NDE
Webersen et al. Guided ultrasonic waves for determining effective orthotropic material parameters of continuous-fiber reinforced thermoplastic plates
CN112362745B (zh) 一种均匀基底涂覆层弹性常数的超声测量与反演方法
US10118111B2 (en) Apparatus, system, and method for removing gas in an immersion ultrasonic process
Ya et al. Residual stresses in laser welded aluminium plate by use of ultrasonic and optical methods
Laurent et al. Laser induced zero-group velocity resonances in transversely isotropic cylinder
US5307680A (en) Method and apparatus for visualization of internal stresses in solid non-transparent materials by elastoacoustic technique
CN110133108A (zh) 一种超精密元件亚表面微纳缺陷测量系统及测量方法
CN111289620B (zh) 各向异性材料的弹性常数检测方法及系统
Ponschab et al. Simulation-based characterization of mechanical parameters and thickness of homogeneous plates using guided waves
Lobkis et al. In-plane elastic property characterization in composite plates
CN107490627B (zh) 超声聚焦探头参数标定方法
CN117169231A (zh) 一种基于声光技术的复合材料无损检测系统
Sathish et al. Quantitative imaging of Rayleigh wave velocity with a scanning acoustic microscope
Lerch et al. Characterization of spherically focused transducers using an ultrasonic measurement model approach
Marsh Finding flaws in composites
Fei et al. Single-scan elastic property estimation in plates
Barnard et al. Towards a generic manual scanner for nondestructive inspection
Georgeson Recent advances in aerospace composite NDE
Fei et al. Rapid dispersion curve mapping and material property estimation in plates
CN116256092A (zh) 一种基于纵横波融合的工件应力分布反演系统、方法及存储介质
Xu et al. Wave velocity measurement in the through-thickness direction of the anisotropic material plate with ultrasonic polar scan

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