首页 / 专利库 / 单位和数量 / 弹性模量 / 一种超声波测量材料内部各处弹性模量及泊松比的方法

一种声波测量材料内部各处弹性模量及泊松比的方法

阅读:1023发布:2020-06-07

专利汇可以提供一种声波测量材料内部各处弹性模量及泊松比的方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种 超 声波 测量材料内部各处 弹性模量 及泊松比的方法,包括以下步骤:S1.利用 超声波 无损检测 对待测材料进行扫查,测量获得超声波在材料中的单程传播时间,并给出材料中超声波的单程传播时间与纵波 波速 VL及横波波速VT的关系;S2.将横波波速VT用 剪切模量 进行表示,将纵波波速VL用剪切模量和泊松比进行表示;S3.确定材料内部各 位置 的剪切模量;S4.确定材料内部各位置的泊松比;S5.计算材料内部各位置的拉伸模量。本发明基于超声无损检测,能够对材料内部空间任意一点位置的弹性模量和泊松比进行测量与表征,能够对材料内部各位置性能的不均匀性进行评估。,下面是一种声波测量材料内部各处弹性模量及泊松比的方法专利的具体信息内容。

1.一种声波测量材料内部各处弹性模量及泊松比的方法,其特征在于:包括以下步骤:
S1.利用超声波无损检测对待测材料进行扫查,获得超声波在材料中的单程传播时间,并给出材料中超声波的单程传播时间与纵波波速VL及横波波速VT的关系:
其中,ΔtT/2为超声波横波的单程传播时间;ΔtL/2为超声波纵波的单程传播时间,d为超声波声束方向的材料厚度;
S2.将横波波速VT用剪切弹性模量表示,将纵波波速VL用剪切弹性模量和泊松比进行表示:
其中,G为剪切弹性模量,ρ为材料密度,ν为泊松比;
S3.确定材料内部各位置的剪切弹性模量;
所述步骤S3包括以下子步骤:
S301.给定超声单程传播路径长度为d,将d离散为N个长度为Δs的一维网格,考虑材料弹性模量的空间不均匀性,对于剪切弹性模量G,令Gk,k∈1,…,N为每个网格位置材料的弹性模量,则超声波横波在d路径上传播的总时间ΔtT为:
表示为积分形式得到:
G(x)是x位置处的剪切弹性模量;
S302.扩展到二维平面,用Γ表示声波半程传播路径,对半程传播路径进行积分得:
其中,Γ为超声在空间中的半程传播路径,二维平面上位置(x,y)处的弹性模量表示为G(x,y),位置(x,y)处的声速为 对半程传播路径Γ积分得半程传播时间ΔtT/2;
S303.对于一个二维扫查平面,给定一个扫查位置k,该位置发出超声的传播路径上点(x,y)利用(r,θ)参数表示:
r=xcosθ+ysinθ
其中r是从原点到超声的传播路径所在直线的垂直距离,θ是该垂直距离所在直线与平轴的度,应用上式对超声的传播路径上所有点的线积分为:
其中δ(·)是狄拉克-δ函数;
对 进行逆变换得到对应的剪
切弹性模量G(x,y);
S304.在待测材料表面沿表面移动探头,超声波声束发射和接收的每一个位置即为一个采样位置,从第一个采样位置开始至最后一个采样位置结束,得到材料内部不同位置的弹性模量的(r,θ)映射二维图像,按照步骤S303得到材料内部不同位置对应的 逆变换后得到材料体内部任意位置处的剪切弹性模量G(x,y);
S4.确定材料内部各位置的泊松比;
所述步骤S4包括:
由步骤S2得到超声纵波波速与剪切弹性模量G和泊松比ν的关系为:
考虑材料内部各处弹性模量的不均匀性,材料内部每一个位置测量得到的(r,θ)对应声束的纵波声速用二维积分表示为:
将材料内部各处的剪切弹性模量G(x,y)在带入
后进行积分
逆变换,得材料内部任一位置处的泊松比v(x,y);
S5.计算材料内部各位置的拉伸模量。
2.根据权利要求1所述的一种超声波测量材料内部各处弹性模量及泊松比的方法,其特征在于:所述步骤S1包括以下子步骤:
S101.给定超声波探头发射位置及声束方向的材料厚度d;对待测材料进行超声波扫描,采集超声波波包的发射时间ts和材料外壁边界反射回的波包接收时间te;
S102.计算出波包在总传播距离下的传播时间:
Δt=te-ts;
则波包的单程传播时间为Δt/2;
S103.波包单程传播时间与超声波波速的关系为:
S104.在超声波探头产生横波和纵波的条件下,分别按照步骤S101~S103获取波包单程传播时间,并给出材料中超声波的单程传播时间与纵波波速VL及横波波速VT的关系:
设超声波探头产生横波时,超声波波包的发射为时间ts1,材料外壁边界反射回的波包接收时间te1,则超声波横波的单程传播时间为 超声波横波单程传播时间与材料中的横波波速关系为:
设超声波探头产生纵波时,超声波波包的发射为时间ts2,材料外壁边界反射回的波包接收时间te2,则超声波纵波的单程传播时间为 超声波纵波单程传播时间与材料中的纵波波速关系为:
3.根据权利要求1所述的一种超声波测量材料内部各处弹性模量及泊松比的方法,其特征在于:所述步骤S2包括以下子步骤:
S201.根据超声波在理想均匀弹性介质中传播的横波波速和泊松比的关系,将横波声速VT用剪切弹性模量G和材料密度ρ表示为:
S202.根据超声波在理想均匀弹性介质中传播的纵波波速和材料弹性模量与泊松比的关系,将纵波声速VL用剪切弹性模量G和泊松比ν及材料密度ρ表示为
S203.根据材料弹性模量之间的关联关系E=2G(1+v),已知剪切弹性模量G和泊松比ν,就能够得到拉伸模量E。
4.根据权利要求1所述的一种超声波测量材料内部各处弹性模量及泊松比的方法,其特征在于:所述步骤S5包括:
根据步骤S3~S4中计算得到的材料内部各位置对应的剪切弹性模量G(x,y)和泊松比v(x,y),利用材料参数间关系E=2G(1+v),得到材料内部任意位置的弹性模量E(x,y):
E(x,y)=2G(x,y)·[1+v(x,y)]。

说明书全文

一种声波测量材料内部各处弹性模量及泊松比的方法

技术领域

[0001] 本发明涉及超声波无损检测,特别是涉及一种超声波测量材料内部各处弹性模量及泊松比的方法。

背景技术

[0002] 弹性模量可视为衡量材料产生弹性变形难易程度的指标,其值越大,使材料发生一定弹性变形的应也越大,即材料刚度越大,亦即在一定应力作用下,发生弹性变形越小。弹性材料的一种最重要、最具特征的力学性质。
[0003] 但目前,弹性模量的测量较为复杂,很难针对到材料内部空间任意一点位置的弹性模量进行测量,为材料性能的评估带来了诸多不便。

发明内容

[0004] 本发明的目的在于克服现有技术的不足,提供一种超声波测量材料内部各处弹性模量及泊松比的方法,能够针对材料内部空间任意一点位置的弹性模量进行测量。
[0005] 本发明的目的是通过以下技术方案来实现的:一种超声波测量材料内部各处弹性模量及泊松比的方法,包括以下步骤:
[0006] S1.利用超声波无损检测对待测材料进行扫查,测量获得超声波在材料中的单程传播时间,并给出材料中超声波的单程传播时间与纵波波速VL及横波波速VT的关系;
[0007] S2.将横波波速VT分别用剪切模量表示,将纵波波速VL用剪切模量和泊松比进行表示;
[0008] S3.确定材料内部各位置的剪切模量;
[0009] S4.确定材料内部各位置的泊松比;
[0010] S5.计算材料内部各位置的拉伸模量。
[0011] 其中,所述步骤S1包括以下子步骤:
[0012] S101.给定超声波探头发射位置及声束方向的材料厚度d;对待测材料进行超声波扫描,采集超声波波包的发射时间ts和材料外壁边界反射回的波包接收时间te;
[0013] S102.计算出波包在总传播距离下的传播时间:
[0014] Δt=te-ts;
[0015] 则波包的单程传播时间为Δt/2;
[0016] S103.波包单程传播时间与超声波波速的关系为:
[0017]
[0018] S104.在超声波探头产生横波和纵波的条件下,分别按照步骤S101~S103进行波包单程传播时间采集,并给出材料中超声波的单程传播时间与纵波波速VL及横波波速VT的关系:
[0019] 设超声波探头产生横波时,超声波波包的发射为时间ts1,材料外壁边界反射回的波包接收时间te1,则超声波横波的单程传播时间为 超声波横波单程传播时间与材料中的横波波速关系为:
[0020]
[0021] 设超声波探头产生纵波时,超声波波包的发射为时间ts2,材料外壁边界反射回的波包接收时间te2,则超声波纵波的单程传播时间为 超声波纵波单程传播时间与材料中的纵波波速关系为:
[0022]
[0023] 其中,所述步骤S2包括以下子步骤:
[0024] S201.根据超声波在理想均匀弹性介质中传播的横波波速和泊松比的关系,将横波声速VT用剪切弹性模量G和材料密度ρ表示为:
[0025]
[0026] S202.根据超声波在理想均匀弹性介质中传播的纵波波速和材料弹性模量与泊松比的关系,将纵波声速VL用剪切弹性模量G和泊松比ν及材料密度ρ表示为[0027]
[0028] S203.根据材料弹性模量之间的关联关系E=2G(1+v),已知剪切模量G和泊松比ν,就能够得到拉伸模量E。
[0029] 其中,所述步骤S3包括以下子步骤:
[0030] S301.给定超声单程传播路径长度为d,将d离散为N个长度为Δs的一维网格,考虑材料弹性模量的空间不均匀性,对于剪切弹性模量G,令Gk,k∈1,…,N为每个网格位置材料的弹性模量,则超声波横波在d路径上传播的总时间ΔtT为:
[0031]
[0032] 表示为积分形式得到:
[0033]
[0034] G(x)是x位置处的剪切弹性模量;
[0035] S302.扩展到二维平面,用Γ表示声波半程传播路径,对传播路径半程进行积分得:
[0036]
[0037] 其中,Γ为超声在空间中的半程传播路径,二维平面上位置(x,y)处的弹性模量表示为G(x, y),位置(x,y)处的声速为 对传播路径Γ积分得半程传播时间ΔtT/2;
[0038] S303.对于一个二维扫查平面,给定一个扫查位置k,该位置发出超声的传播路径上点(x, y)利用(r,θ)参数表示:
[0039] r=xcosθ+ysinθ
[0040] 其中r是从原点到声速所在直线的垂直距离,θ是该垂直距离所在直线与平轴的度,应用上式对声束传播路径上所有点的线积分为:
[0041]
[0042] 其中δ(·)是狄拉克-δ函数;对上式进行逆变换得到对应的剪切弹性模量G(x,y);
[0043] S304.在待测材料表面沿表面移动探头,超声波声束发射和接收的每一个位置即为一个采样点,从第一个采样位置开始至最后一个采样位置结束,得到材料内部不同位置的弹性模量的(r,θ)映射二维图像,按照步骤S303得到材料内部不同位置对应的 表达式,结合各个采样点下的横波边到边传播时间(由超声波探头按照步骤S1测得,即),逆变换后得到材料体内部任意位置处的剪切弹性模量G(x,y);
[0044] 其中,所述步骤S4包括:
[0045] 由步骤S2得到超声纵波波速与剪切模量G和泊松比v的关系为:
[0046]
[0047] 考虑材料内部各处弹性模量的不均匀性,材料内部每一个位置测量得到的(r,θ)对应声束的纵波声速用二维积分表示为:
[0048]
[0049] 将材料内部各处的剪切模量G(x,y)在带入上式进行积分逆变换,得材料内部任一位置处的泊松比v(x,y), 也就是各个采样点下的纵波边到边传播时间(由超声波探头按照步骤S1测得)。
[0050] 其中,所述步骤S5包括:
[0051] 根据步骤S3~S4中计算得到的材料内部各位置对应的剪切模量G(x,y)和泊松比v(x, y),利用材料参数间关系E=2G(1+v),得到材料内部任意位置的弹性模量E(x,y):
[0052] E(x,y)=2G(x,y)·[1+v(x,y)]。
[0053] 本发明的有益效果是:本发明基于超声声速无损检测和材料泊松比,能够对针对材料内部空间任意一点位置的弹性模量进行测量,有助于对材料性能的综合评估。附图说明
[0054] 图1为本发明的方法流程图
[0055] 图2为超声波无损检测的示意图;
[0056] 图3为波包在总传播距离下的传播时间计算示意图;
[0057] 图4为将超声单程传播路径长度为d离散为N个长度为Δs的一维网格示意图;
[0058] 图5为二维平面上声波半程传播路径的示意图;
[0059] 图6为二维扫查平面的示意图;
[0060] 图7为半径为50cm的圆柱体横波播时间的试验结果示意图;
[0061] 图8为试验数据进行积分逆变换得到的直径50cm棒体一个扫查截面的材料内部各位置处弹性模量值。

具体实施方式

[0062] 下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
[0063] 如图1所示,一种超声波测量材料内部各处弹性模量及泊松比的方法,包括以下步骤:
[0064] S1.利用超声波无损检测对待测材料进行扫查,测量获得超声波在材料中的单程传播时间,并给出材料中超声波的单程传播时间与纵波波速VL及横波波速VT的关系;
[0065] S2.将横波波速VT用剪切模量表示,将纵波波速VL用剪切模量和泊松比进行表示;
[0066] S3.确定材料内部各位置的剪切模量;
[0067] S4.确定材料内部各位置的泊松比;
[0068] S5.计算材料内部各位置的拉伸模量。
[0069] 其中,所述步骤S1包括以下子步骤:
[0070] S101.如图2所示,给定超声波探头发射位置及声束方向的材料厚度d;对待测材料进行超声波扫描,采集超声波波包的发射时间ts和材料外壁边界反射回的波包接收时间te;
[0071] S102.计算出波包在总传播距离下的传播时间:
[0072] Δt=te-ts;
[0073] 如图3所示,为Δt的计算示意图。则波包的单程传播时间为Δt/2;
[0074] S103.波包单程传播时间与超声波波速的关系为:
[0075]
[0076] S104.在超声波探头产生横波和纵波的条件下,分别按照步骤S101~S103进行波包单程传播时间采集,并给出材料中超声波的单程传播时间与纵波波速VL及横波波速VT的关系:
[0077] 设超声波探头产生横波时,超声波波包的发射为时间ts1,材料外壁边界反射回的波包接收时间te1,则超声波横波的单程传播时间为 超声波横波单程传播时间与材料中的横波波速关系为:
[0078]
[0079] 设超声波探头产生纵波时,超声波波包的发射为时间ts2,材料外壁边界反射回的波包接收时间te2,则超声波纵波的单程传播时间为 超声波纵波单程传播时间与材料中的纵波波速关系为:
[0080]
[0081] 其中,所述步骤S2包括以下子步骤:
[0082] S201.根据超声波在理想均匀弹性介质中传播的横波波速和泊松比的关系,将横波声速VT用剪切弹性模量G和材料密度ρ表示为:
[0083]
[0084] S202.根据超声波在理想均匀弹性介质中传播的纵波波速和材料弹性模量与泊松比的关系,将纵波声速VL用剪切弹性模量G和泊松比v及材料密度ρ表示为[0085]
[0086] S203.根据材料弹性模量之间的关联关系E=2G(1+v),已知剪切模量G和泊松比v,就能够得到拉伸模量E。
[0087] 其中,所述步骤S3包括以下子步骤:
[0088] S301.给定超声单程传播路径长度为d,将d离散为N个长度为Δs的一维网格,考虑材料弹性模量的空间不均匀性,对于剪切弹性模量G,令Gk,k∈1,...,N为每个网格位置材料的弹性模量,则超声波横波在d路径上传播的总时间ΔtT为:
[0089]
[0090] 表示为积分形式得到:
[0091]
[0092] G(x)是x位置处的剪切弹性模量;
[0093] S302.扩展到二维平面,用Γ表示声波半程传播路径,对传播路径半程进行积分得:
[0094]
[0095] 其中,Γ为超声在空间中的半程传播路径,二维平面上位置(x,y)处的弹性模量表示为G(x, y),位置(x,y)处的声速为 对传播路径Γ积分得半程传播时间ΔtT/2;
[0096] S303.对于一个二维扫查平面,给定一个扫查位置k,该位置发出超声的传播路径上点(x, y)利用(r,θ)参数表示:
[0097] r=xcosθ+ysinθ
[0098] 其中r是从原点到声速所在直线的垂直距离,θ是该垂直距离所在直线与水平轴的角度,应用上式对声束传播路径上所有点的线积分为:
[0099]
[0100] 其中δ(·)是狄拉克-δ函数;对上式进行逆变换得到对应的剪切弹性模量G(x,y);
[0101] S304.在待测材料表面沿表面移动探头,超声波声束发射和接收的每一个位置即为一个采样点,从第一个采样位置开始至最后一个采样位置结束,得到材料内部不同位置的弹性模量的(r,θ)映射二维图像,按照步骤S303得到材料内部不同位置对应的 逆变换后得到材料体内部任意位置处的剪切弹性模量G(x,y);
[0102] 如图7所示,为半径为50cm的圆柱体横波播时间的试验结果,对一个圆柱体一个轴向位置所在的平面扫查结果。横坐标为超声波束正交方向的角度θ,纵坐标为超声波束与原点的正交距离ρ。图像中每一个点的颜色即代表该位置对应(r,θ)的横波边到边传播时间(也就是横波在总传播距离下的传播时间,在超声波探头探测过程中,按照步骤S1测得)的相对大小;需要说明的是,图7的试验数据结果,对应的 的图像,利用试验数据,对 应用逆变换得到材料体内部任意位置处的弹性模量G(x,y)。在具体实施过程中,逆变换的本质是逆积分变换,可以采用逆Radon变换,或其他方法实现,例如基于傅立叶切片理论(Fourier slice theorem,也叫Projection slice theorem、central slice theorem)、滤波反投影(Filtered  backprojection)、希尔伯特变换(Hillbert transform)、反投影后滤波(Filtering after backprojection)等方法。利用图7中显示的试验数据进行逆变换后得到G(x,y)在网格各单元中的值,可以用图像表示,如图8所示,不同灰度颜色表示了材料弹性模量G(x,y)市目对于平均值的百分比大小。
[0103] 其中,所述步骤S4包括:
[0104] 由步骤S2得到超声纵波波速与剪切模量G和泊松比v的关系为:
[0105]
[0106] 考虑材料内部各处弹性模量的不均匀性,材料内部每一个位置测量得到的(r,θ)对应声束的纵波声速用二维积分表示为:
[0107]
[0108] 将材料内部各处的剪切模量G(x,y)在带入上式进行积分逆变换,得材料内部任一位置处的泊松比v(x,y), 也就是各个采样点下的纵波边到边传播时间(也就是纵波在总传播距离下的传播时间,在超声波探头探测过程中,按照步骤S1测得)。
[0109] 其中,所述步骤S5包括:
[0110] 根据步骤S3~S4中计算得到的材料内部各位置对应的剪切模量G(x,y)和泊松比v(x, y),利用材料参数间关系E=2G(1+v),得到材料内部任意位置的弹性模量E(x,y):
[0111] E(x,y)=2G(x,y)·[1+v(x,y)]。
[0112] 以上所述是本发明的优选实施方式,应当理解本发明并非局限于本文所披露的形式,不应该看作是对其他实施例的排除,而可用于其他组合、修改和环境,并能够在本文所述构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈