CN112172190A - 一种无损检测指导复合材料铺放工艺在线调整系统及方法 - Google Patents

一种无损检测指导复合材料铺放工艺在线调整系统及方法 Download PDF

Info

Publication number
CN112172190A
CN112172190A CN202010936787.8A CN202010936787A CN112172190A CN 112172190 A CN112172190 A CN 112172190A CN 202010936787 A CN202010936787 A CN 202010936787A CN 112172190 A CN112172190 A CN 112172190A
Authority
CN
China
Prior art keywords
laying
signal
defect
laser
nondestructive testing
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
CN202010936787.8A
Other languages
English (en)
Other versions
CN112172190B (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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN202010936787.8A priority Critical patent/CN112172190B/zh
Publication of CN112172190A publication Critical patent/CN112172190A/zh
Application granted granted Critical
Publication of CN112172190B publication Critical patent/CN112172190B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B29WORKING OF PLASTICS; WORKING OF SUBSTANCES IN A PLASTIC STATE IN GENERAL
    • B29CSHAPING OR JOINING OF PLASTICS; SHAPING OF MATERIAL IN A PLASTIC STATE, NOT OTHERWISE PROVIDED FOR; AFTER-TREATMENT OF THE SHAPED PRODUCTS, e.g. REPAIRING
    • B29C70/00Shaping composites, i.e. plastics material comprising reinforcements, fillers or preformed parts, e.g. inserts
    • B29C70/04Shaping composites, i.e. plastics material comprising reinforcements, fillers or preformed parts, e.g. inserts comprising reinforcements only, e.g. self-reinforcing plastics
    • B29C70/28Shaping operations therefor
    • B29C70/30Shaping by lay-up, i.e. applying fibres, tape or broadsheet on a mould, former or core; Shaping by spray-up, i.e. spraying of fibres on a mould, former or core
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B29WORKING OF PLASTICS; WORKING OF SUBSTANCES IN A PLASTIC STATE IN GENERAL
    • B29CSHAPING OR JOINING OF PLASTICS; SHAPING OF MATERIAL IN A PLASTIC STATE, NOT OTHERWISE PROVIDED FOR; AFTER-TREATMENT OF THE SHAPED PRODUCTS, e.g. REPAIRING
    • B29C70/00Shaping composites, i.e. plastics material comprising reinforcements, fillers or preformed parts, e.g. inserts
    • B29C70/04Shaping composites, i.e. plastics material comprising reinforcements, fillers or preformed parts, e.g. inserts comprising reinforcements only, e.g. self-reinforcing plastics
    • B29C70/28Shaping operations therefor
    • B29C70/54Component parts, details or accessories; Auxiliary operations, e.g. feeding or storage of prepregs or SMC after impregnation or during ageing
    • 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 sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • G01N21/88Investigating the presence of flaws or contamination
    • G01N21/8851Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges
    • 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 sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • G01N21/88Investigating the presence of flaws or contamination
    • G01N21/8851Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges
    • G01N2021/8887Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges based on image processing techniques
    • 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 sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • G01N21/88Investigating the presence of flaws or contamination
    • G01N21/8851Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges
    • G01N2021/8896Circuits specially adapted for system specific signal conditioning

Landscapes

  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Composite Materials (AREA)
  • Biochemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Signal Processing (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)

Abstract

本发明涉及一种无损检测指导复合材料铺放工艺在线调整系统及方法,包括无损检测系统,所述的无损检测系统包括对工控机、示波器和无损检测装置;所述的无损检测装置包括设置在铺放样品前方的YAG激光器、在YAG激光器的前方先后放有滤光透镜和偏振分束器、在偏振分束器前方先后放有反光镜和聚光透镜;样品后方设有双波混合干涉仪和光电探测器。根据无损检测系统检测结果的图像处理,获得最新铺放层与下一层之间的层内缺陷,进而对铺放工艺进行调整。本发明具有如下优点:1)实现了复合材料铺放过程中最新铺放层与下一之间的层内缺陷在线检测;2)实现实时根据复合材料铺放层内缺陷对铺放工艺的在线调整。

Description

一种无损检测指导复合材料铺放工艺在线调整系统及方法
技术领域
本发明涉及一种复合材料铺放工艺参数调整系统及方法,尤其是涉及一种无损检测指导复合材料铺放工艺在线调整系统及方法。
背景技术
碳纤维复合材料具有比强度大、弹性模量高、耐疲劳和耐化学腐蚀等优良性能,被广泛地应用在航空航天和汽车船舶等领域。复合材料自动铺放具有自动化和高效率的特点,是目前较为常用的复合材料成形技术,复合材料铺放工艺中常需要控制的工艺参数包括铺放温度、铺放压力和铺放速度。复合材料采用铺放成形时受到铺放工艺参数组合的影响会出现褶皱、重叠、气泡和间隙等缺陷,影响构件力学性能。
激光超声是一种非接触、高精度、无损伤的新型检测技术,检测可以不用耦合介质,且检测速度快,非常适合用于对复合材料构件的缺陷进行无损检测,但目前激光超声无损检测主要用来检测铺放完成后的构件的整体缺陷。由于铺放是逐层进行的,后续的铺层会受到先前铺放层的影响,铺放环境如环境湿度和温度等变化也会对铺放质量产生影响,并且会存在干扰使得某项工艺参数无法达到设定值,使得当前的工艺参数组合不再是最佳铺放工艺参数组合,进而需要通过其它工艺参数的在线调整,重新获得最佳工艺参数组合,避免出现铺放缺陷。目前控制工艺参数达到设定目标的控制方法,无法自适应工艺参数在铺放过程中在线调整的需求。
发明内容
本发明主要是解决现有技术所存在的技术问题;提供了一种用于复合材料在线铺放过程中实时获取铺放质量的无损检测系统及方法,且不影响铺放进程,可实现对复合材料自动铺放过程中的层内缺陷进行快速检测和成像,获得铺放工艺参数对铺放缺陷的影响,并可指导铺放工艺参数的在线调整。
本发明的上述技术问题主要是通过下述技术方案得以解决的:
本发明的技术方案如下:
一种无损检测指导复合材料铺放工艺在线调整系统,其特征在于,包括无损检测系统,无损检测系统包括无损检测装置以及与无损检测装置连接的工控机,控制无损检测装置纵向移动的纵向移动滑台。
在上述的一种无损检测指导复合材料铺放工艺在线调整系统,所述的无损检测装置包括YAG激光器、位于YAG激光器前方的偏振分束器、位于YAG激光器和偏振分束器之间的滤光透镜、位于YAG激光器前面的反光镜、位于铺放样件与聚光透镜之间的反光镜、位于铺放样件前方的双波混合干涉仪和对双波混合干涉仪进行控制的控制器,工控机对偏振分束器、光电探测器和控制器进行控制。
在上述的一种无损检测指导复合材料铺放工艺在线调整系统,激光超声无损对刚铺放完成的复合材料构件的检测光路走向依次是:YAG激光器发射的激光通过滤光透镜后通过偏振分束器分成两束激光,其中一束激光先通过反光镜反射后射到聚光透镜,然后射到刚刚完成的铺放样件的表面,铺放样件的表面反射的激光束即反射光,带有缺陷声波的反射光与偏振分束器分出的另一束激光即参考光,一起射到双波混合干涉仪中进行干涉;光电探测器将反射光中由于超声振动引起的相位解调出来,并将该信号传送到工控机。
一种无损检测指导复合材料铺放工艺在线调整方法,其特征在于,包括:
步骤A.纵向移动滑台底部与常规的复合材料自动铺放装置连接,从而实现无损检测装置随着复合材料铺放过程一起运动,实现无损检测装置对铺放完成区域层内的缺陷在线检测。
步骤B.在复合材料自动铺放时,工控机控制YAG激光器沿着垂直于复合材料自动铺放方向发射激光,激光先后通过滤光透镜和偏振分束器之后变成两束激光,其中一束先后通过反光镜照射到铺放样件表面,其反射光与另一束激光在双波混合干涉仪发生干涉。
步骤C.控制器对光电探测器进行控制,光电探测器将反射光中带有的超声振动信号解调出来,并将解调的信号传送到工控机。
步骤D.采用信号匹配追踪法对信号噪声进行初步去除噪声。
步骤D1、选用由调制后的高斯窗函数ψr(t)构成Gabor时频原子构成过完备原子库,ψr(t)的表达式为:
Figure BDA0002672230370000021
其中,t表示时间参数,s、u、v、w分别表示尺度参数(尺度因子)、原子的位移参数(位移因子)、原子的频率参数(频率)和原子的相位参数(相位),这些参数构成信号的时频参数r=(s,u,v,w),并将其离散化:
r={(αj,pαjΔu,kα-jΔv,iΔw)}
其中,时频参数r的离散化中的参数α,j,p,k分别为:α=2,0<j≤log2N-1,0≤p≤2-j+1N,0≤k≤2j+1,Δu、Δv和Δw取值分别为:Δu=1/2,Δv=π,Δw=π/6,N为信号长度,0≤i≤12;
步骤D2、从过程完备原子库中选出与带分解信号x最为匹配的原子ψr0,Γ为参数r的集合,其中ψr0满足以下条件:
Figure BDA0002672230370000031
将信号分解为在最佳原子
Figure BDA0002672230370000034
上的分量和残余两部分,用R1x表示最佳原子对原信号进行最佳匹配后的残余,即为:
Figure BDA0002672230370000035
步骤D3、对最佳匹配后的残余按上式的形式进行重复的k次迭代分解:
Figure BDA0002672230370000036
其中
Figure BDA0002672230370000037
是每次迭代后的最佳原子与上一次信号残余的相关系数,且
Figure BDA0002672230370000038
满足
Figure BDA0002672230370000039
经n步分解后,信号被分解为:
Figure BDA0002672230370000032
信号长度有限时,经过n次分解后,残差信号的能量递减最后收敛到零,此时信号中有用的成分已经被全部提取,余下的部分为噪声,此时完成对信号的初步去噪,并且使得含有缺陷振动的波形特征较好的保存下来,最终信号被分解为:
Figure BDA0002672230370000033
步骤E、采用分谱处理法对匹配信号法对去除噪声后保留下来的缺陷信号进行分谱处理,进一步去除噪声,并对信号进行特征提取;
步骤E1、利用快速傅里叶变换获得原始信号的x(t)的频率—能量谱:
Figure BDA0002672230370000041
其中X(ω)为x(t)的傅里叶变换结果,w代表频率,t代表时间,e-iωt为复变函数。
将能量谱最大值的一半所对应的起始频率和截止频率分别用于指定信号的主要频率带,然后利用一组带通滤波器把所制定的主要频率带切分为一系列窄的频率带Yl(ω),l=1~M,Yl(ω)为第l个频率带,设定有总数为M的一组高斯带通滤波器,它们具有不同的中心频率fl,l=1~M,相邻的两个带通滤波器的中心频率间隔为Δf,相同的带宽2×Δf;
Δf=1/T;M=BT+1
其中T是信号的采样时间,即信号总的时间长度;B为主要频率带的带宽,为使得相邻的带通滤波器之间有足够的混叠,取B=4/T;
步骤E2、将所切分的M个窄频率带Yl(ω),l=1~M,各自进行傅里叶逆变换,重建各个时间域信号yl(t):
Figure BDA0002672230370000042
步骤E3、结合极性阈值法与最小值法恢复信号,设恢复后信号为f(t):
Figure BDA0002672230370000043
步骤F、在完成对缺陷信号处理后,采用合成孔径聚焦成像方法对缺陷成像;
步骤F1、对刚刚铺放完成区域内部点C进行重建
Figure BDA0002672230370000044
Om,m=1,2,3,…,P为激发光的位置,mesh为三维网格;d0和dm分别为检测光即为反射光和激发光与缺陷之间的距离;xe和xd分别为初始激发点和固定监测点与缺陷之间的横向距离,用v表示纵波的速度,探测光在
Figure BDA0002672230370000045
时刻探测到的信号S(Om,t)将会出现一个有缺陷引起的反射波峰,将各S(Om,t)累加起来,m=1,2,3…P,将重建图形中C点所在位置的值增强,任意一点的传播距离D为:
Figure BDA0002672230370000051
其中
Figure BDA0002672230370000052
为激发光和缺陷之间的距离,
Figure BDA0002672230370000053
为检测光和缺陷之间的距离,用q(q=0,1,2,…,P-1)表示激发点扫描过的步数,Δd表示每一步之间的间隔,xt=nΔd,xe为初始激发点与缺陷之间的横向距离,xd为检测点与缺陷之间的横向距离,z为缺陷的深度;
步骤F2、当激光在刚铺放的表面检测时,被缺陷反射的体纵波将在横轴为扫查步数q,q=0,1,2,…,P-1,纵轴为时间t的B扫图上呈现双曲线的一支,将不同双曲线的各个值进行相加,所得结果即可反应C点的性质,获得垂直于铺放方向的刚刚铺放完成的缺陷的横向截面图,按照时间序列将横向截面图叠加获得空间图像,根据复合材料铺放过程中厚度逐渐累加的特性,在厚度方向上提取一定厚度的空间图像,即为最新铺层与下一层之间的缺陷图像;
步骤G、在大量试验的基础上获得铺放工艺参数包括铺放温度、压力和速度与铺放缺陷之间的关系,获得较好铺放质量时的工艺参数组合,选取一组铺放质量较好的工艺参数组合作为铺放工艺参数初始值,采用上述的激光无损检测方式获得在线铺放过程中最新铺层与下一层之间的缺陷信息,当出现缺陷超过铺放质量可接受的缺陷范围时,当前的工艺参数组合已经不再适合铺放;
根据铺放过程中压力传感器、温度传感器和速度传感器等获得的工艺参数值实时大小和铺放工艺参数组合与铺放质量之间的关系,对铺放工艺参数进行在线调整,重新获得适合当前铺放条件的工艺参数组合,减小铺放缺陷。
因此,本发明具有如下优点:1.实现了对复合材料自动铺放过程中层间铺放缺陷的在线检测;2.实现根据层间缺陷对铺放工艺参数在线调整。
附图说明
图1是本发明所提供的采用激光超声无损检测方式的复合材料铺放缺陷在线检测系统示意图;
图2为本发明所提供的用激光超声无损在线检测复合材料层内缺陷指导工艺参数实时调整流程图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。图中,无损检测系统1、工控机2、YAG激光器3、滤光透镜4、偏振分束器5、反光镜6、聚光透镜7、双波混合干涉仪8、控制器9、光电探测器10、无损检测装置11、纵向移动滑台12、铺放样件13。
实施例:
参照图1,一种无损检测指导复合材料铺放工艺在线调整方法,其包括无损检测系统1,所述的无损检测系统1包括无损检测装置11、工控机2和控制无损检测装置11纵向移动的纵向移动滑台12。
无损检测装置11包括YAG激光器3、位于YAG激光器3前方的偏振分束器5、位于YAG激光器3和偏振分束器5之间的滤光透镜4、位于YAG激光器前面的反光镜6、位于铺放样件13与聚光透镜7之间的反光镜6、位于铺放样件13前方的双波混合干涉仪8。
YAG激光器3发射的激光通过滤光透镜4后通过偏振分束器5分成两束激光,其中一束激光先后通过反光镜6和聚光透镜7后射到自动铺放装置11刚刚完成的铺放样件13的表面,其激光束的反射光(即信号光)与偏振分束器5分出的另一束激光(即参考光),一起射到双波混合干涉仪8中,干涉形成动态光栅;所述的光电探测器10将接收到的超声振动信号解调出来。
参照图2,一种用激光超声无损在线检测指导复合材料铺放工艺调整方法,包括以下步骤:
A.纵向移动滑台12对无损检测装置11与复合材料样件的相对位置调整到合适的距离,纵向移动滑台12的底部与常规的自动铺放装置连接在一起,实现无损检测装置11随着复合材料铺放过程沿着复合材料铺放方向移动,实现对刚刚铺放完成区域层内的缺陷实时检测。
B.在复合材料自动铺放时,工控机2控制YAG激光器3发射沿着垂直于复合材料自动铺放方向的激光束,激光先后通过滤光透镜4和偏振分束器5之后变成两束激光,其中一束先后通过反光镜6后照射到铺放样件表面,其反射光携带缺陷信息与另一束激光在双波混合干涉仪8发生干涉。
C.控制器9对光电探测器10进行控制,同时光电探测器10将反射光中带有的超声振动信号解调出来,并将解调出来的信号传送到工控机2。
D.采用信号匹配追踪法对信号噪声进行初步去除噪声。
D1.选用由调制后的高斯窗函数ψr(t)构成Gabor时频原子构成过完备原子库,ψr(t)的表达式为:
Figure BDA0002672230370000071
其中,s、u、v、w分别表示尺度参数(尺度因子)、原子的位移参数(位移因子)、原子的频率参数(频率)和相位参数(相位),这些参数构成信号的时频参数r=(s,u,v,w),并将其离散化:
r={(αj,pαjΔu,kα-jΔv,iΔw)}
其中,α=2,Δu=1/2,Δv=π,Δw=π/6,0<j≤log2N-1,0≤p≤2-j+1N,0≤k≤2j +1,0≤i≤12,N为信号的长度。
D2.从过程完备原子库中选出与带分解信号x最为匹配的原子ψr0,Γ为参数r的集合,其中ψr0满足以下条件:
Figure BDA0002672230370000072
将信号分解为在最佳原子
Figure BDA0002672230370000073
上的分量和残余两部分,用R1x表示最佳原子对原信号进行最佳匹配后的残余,即为:
Figure BDA0002672230370000074
D3.对最佳匹配后的残余按上式的形式进行重复的k次迭代分解:
Figure BDA0002672230370000075
其中
Figure BDA0002672230370000076
是每次迭代后的最佳原子与上一次信号残余的相关系数,且
Figure BDA0002672230370000077
满足
Figure BDA0002672230370000078
经n步分解后,信号被分解为:
Figure BDA0002672230370000079
信号长度有限时,经过n次分解后,残差信号的能量递减最后收敛到零,此时信号中有用的成分已经被全部提取,余下的部分为噪声,此时完成对信号的初步去噪,并且使得含有缺陷的波形特征较好的保存下来,最终信号被分解为:
Figure BDA0002672230370000081
E.采用分谱处理法对匹配信号法对去除噪声后保留下来的缺陷信号进行分谱处理,进一步去除噪声,并对信号进行特征提取。
E1.利用快速傅里叶变换获得原始信号的x(t)的频率—能量谱:
Figure BDA0002672230370000082
其中X(ω)为x(t)的傅里叶变换结果,w代表频率,t代表时间,e-iωt为复变函数。
将能量谱最大值的一半所对应的起始频率和截止频率分别用于指定信号的主要频率带,然后利用一组带通滤波器把所制定的主要频率带切分为一系列窄的频率带(Yl(ω),l=1~M),Yl(ω)为第l个频率带,设定有总数为M的一组高斯带通滤波器,它们具有不同的中心频率(fl,l=1~M),相邻的两个带通滤波器的中心频率间隔为Δf,相同的带宽2×Δf。其中Δf=1/T;M=BT+1,T是信号的采样时间,即信号总的时间长度;B为主要频率带的带宽,为使得相邻的带通滤波器之间有足够的混叠,取B=4/T。
E2.将所切分的M各窄频率带(Yl(ω),l=1~M)各自进行傅里叶逆变换,重建各个时间域信号yl(t):
Figure BDA0002672230370000083
E3.结合极性阈值法与最小值法恢复缺陷信号,设恢复后缺陷信号为f(t):
Figure BDA0002672230370000084
F.在完成对缺陷信号处理后,采用合成孔径聚焦成像方法对缺陷成像。
F1.对刚刚铺放完成区域内部点C进行重建
Figure BDA0002672230370000085
Om(m=1,2,3,…,P)为激发光的位置;d0和dm分别为检测光即为反射光和激发光与缺陷之间的距离;xe和xd分别为初始激发点和固定监测点与缺陷之间的横向距离。用v表示纵波的速度,探测光在
Figure BDA0002672230370000091
时刻探测到的信号S(Om,t)将会出现一个有缺陷引起的反射波峰,将各S(Om,t)(m=1,2,3…P)累加起来,将重建图形中C点所在位置的值增强。任意一点的传播距离D为:
Figure BDA0002672230370000092
其中
Figure BDA0002672230370000093
为激发光和缺陷之间的距离,
Figure BDA0002672230370000094
为检测光和缺陷之间的距离,用q(q=0,1,2,…,P-1)表示激发点扫描过的步数,Δd表示每一步之间的间隔,xt=nΔd,xe为初始激发点与缺陷之间的横向距离,xd为检测点与缺陷之间的横向距离,z为缺陷的深度;
F2.当激光在刚铺放的表面检测时,被缺陷反射的体纵波将在横轴为扫查步数q(q=0,1,2,…,P-1,)纵轴为时间t的B扫图上呈现双曲线的一支。将不同双曲线的各个值进行相加,所得结果即可反应C点的性质,获得垂直于铺放方向的刚刚铺放完成的缺陷的横向截面图,按照时间序列将横向截面图叠加获得空间图像。根据复合材料铺放过程中厚度逐渐累加的特性,在厚度方向上提取一定厚度的空间图像,即为最新铺层与下一层之间的缺陷图像,在工控机2中的显示器中进行显示。
G.在大量试验的基础上获得铺放工艺参数包括铺放温度、压力和速度与铺放缺陷之间的关系,获得较好铺放质量时的工艺参数组合。选取一组铺放质量较好的工艺参数组合作为铺放工艺参数初始值。采用上述的激光超声无损检测方式获得在线铺放过程中最新铺层与下一层之间的缺陷信息,当出现缺陷超过铺放质量可接受的缺陷范围时,当前的工艺参数组合已经不再适合铺放。
根据铺放过程中压力传感器、温度传感器和速度传感器等获得的工艺参数值实时大小和铺放工艺参数组合与铺放质量之间的关系,对铺放工艺参数进行在线调整,重新获得适合当前铺放条件的工艺参数组合,减小铺放缺陷。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

Claims (4)

1.一种无损检测指导复合材料铺放工艺在线调整系统,其特征在于,包括无损检测系统(1),无损检测系统(1)包括无损检测装置(13)以及与无损检测装置(11)连接的工控机(2),控制无损检测装置(13)纵向移动的纵向移动滑台(12)。
2.根据权利要求1所述的系统,其特征在于,无损检测装置(11)包括YAG激光器(3)、位于YAG激光器(3)前方的偏振分束器(5)、位于YAG激光器(3)和偏振分束器(5)之间的滤光透镜(4)、位于YAG激光器前面的反光镜(6)、位于铺放样件(13)与聚光透镜(7)之间的反光镜(6)、位于铺放样件(13)前方的双波混合干涉仪(8)和对双波混合干涉仪(8)进行控制的控制器(9),工控机(2)对偏振分束器(5)、光电探测器(10)和控制器(9)进行控制。
3.根据权利要求2的系统,其特征在于,激光无损对刚铺放完成的复合材料构件的检测光路走向依次是:YAG激光器(3)发射的激光通过滤光透镜(4)后通过偏振分束器(5)分成两束激光,其中一束激光先通过反光镜(6)反射后射到聚光透镜(7),然后射到刚刚完成的铺放样件(13)的表面,铺放样件(13)的表面反射的激光束即反射光,带有缺陷声波的反射光与偏振分束器(5)分出的另一束激光即参考光,一起射到双波混合干涉仪(8)中进行干涉;光电探测器(10)将反射光中由于超声振动引起的相位解调出来,并将该信号传送到工控机(2)。
4.一种无损检测指导复合材料铺放工艺在线调整方法,其特征在于,包括:
步骤A、纵向移动滑台底部与常规的复合材料自动铺放装置连接,从而实现无损检测装置(11)随着复合材料铺放过程运动,实现无损检测装置对铺放完成区域层内的缺陷实时检测;
步骤B、在复合材料自动铺放时,工控机(2)控制YAG激光器(3)沿着垂直于复合材料自动铺放方向发射激光,激光先后通过滤光透镜(4)和偏振分束器(5)之后变成两束激光,其中一束先后通过反光镜(6)照射到铺放样件表面,其反射光与另一束激光在双波混合干涉仪(8)发生干涉;
步骤C.控制器(9)对光电探测器(10)进行控制,光电探测器(10)将反射光中带有的超声振动信号解调出来,并将解调的信号传送到工控机(2);
步骤D、采用信号匹配追踪法对信号噪声进行初步去除噪声;
步骤D1、选用由调制后的高斯窗函数ψr(t)构成Gabor时频原子构成过完备原子库,ψr(t)的表达式为:
Figure FDA0002672230360000021
其中,t表示时间参数,s、u、v、w分别表示尺度参数、位移参数、频率参数和相位参数,这些参数构成信号的时频参数r=(s,u,v,w),并将其离散化:
r={(αj,pαjΔu,kα-jΔv,iΔw)}
其中,时频参数r的离散化中的参数α,j,p,k分别为:α=2,0<j≤log2N-1,0≤p≤2-j+ 1N,0≤k≤2j+1,Δu、Δv和Δw取值分别为:Δu=1/2,Δv=π,Δw=π/6,N为信号长度,0≤i≤12;
步骤D2、从过程完备原子库中选出与带分解信号x最为匹配的原子ψr0,Γ为参数r的集合,其中ψr0满足以下条件:
Figure FDA0002672230360000022
将信号分解为在最佳原子
Figure FDA0002672230360000023
上的分量和残余两部分,用R1x表示最佳原子对原信号进行最佳匹配后的残余,即为:
Figure FDA0002672230360000024
步骤D3、对最佳匹配后的残余按上式的形式进行重复的k次迭代分解:
Figure FDA0002672230360000025
其中
Figure FDA0002672230360000026
是每次迭代后的最佳原子与上一次信号残余的相关系数,且
Figure FDA0002672230360000027
满足
Figure FDA0002672230360000028
经n步分解后,信号被分解为:
Figure FDA0002672230360000029
信号长度有限时,经过n次分解后,残差信号的能量递减最后收敛到零,此时信号中有用的成分已经被全部提取,余下的部分为噪声,此时完成对信号的初步去噪,并且使得含有缺陷振动的波形特征较好的保存下来,最终信号被分解为:
Figure FDA00026722303600000210
步骤E、采用分谱处理法对匹配信号法对去除噪声后保留下来的缺陷信号进行分谱处理,进一步去除噪声,并对信号进行特征提取;
步骤E1、利用快速傅里叶变换获得原始信号的x(t)的频率—能量谱:
Figure FDA0002672230360000031
其中X(ω)为x(t)的傅里叶变换结果,w代表频率,t代表时间,e-iωt为复变函数;
将能量谱最大值的一半所对应的起始频率和截止频率分别用于指定信号的主要频率带,然后利用一组带通滤波器把所制定的主要频率带切分为一系列窄的频率带Yl(ω),l=1~M,Yl(ω)为第l个频率带,设定有总数为M的一组高斯带通滤波器,它们具有不同的中心频率fl,l=1~M,相邻的两个带通滤波器的中心频率间隔为Δf,相同的带宽2×Δf;
Δf=1/T;M=BT+1
其中T是信号的采样时间,即信号总的时间长度;B为主要频率带的带宽,为使得相邻的带通滤波器之间有足够的混叠,取B=4/T;
步骤E2、将所切分的M个窄频率带Yl(ω),l=1~M,各自进行傅里叶逆变换,重建各个时间域信号yl(t):
Figure FDA0002672230360000032
步骤E3、结合极性阈值法与最小值法恢复信号,设恢复后信号为f(t):
Figure FDA0002672230360000033
步骤F、在完成对缺陷信号处理后,采用合成孔径聚焦成像方法对缺陷成像;
步骤F1、对刚刚铺放完成区域内部点C进行重建
Figure FDA0002672230360000034
Om,m=1,2,3,…,P为激发光的位置,mesh为三维网格;d0和dm分别为检测光即为反射光和激发光与缺陷之间的距离;xe和xd分别为初始激发点和固定监测点与缺陷之间的横向距离,用v表示纵波的速度,探测光在
Figure FDA0002672230360000035
时刻探测到的信号S(Om,t)将会出现一个有缺陷引起的反射波峰,将各S(Om,t)累加起来,m=1,2,3…P,将重建图形中C点所在位置的值增强,任意一点的传播距离D为:
Figure FDA0002672230360000041
其中
Figure FDA0002672230360000042
为激发光和缺陷之间的距离,
Figure FDA0002672230360000043
为检测光和缺陷之间的距离,用q(q=0,1,2,…,P-1)表示激发点扫描过的步数,Δd表示每一步之间的间隔,xt=nΔd,xe为初始激发点与缺陷之间的横向距离,xd为检测点与缺陷之间的横向距离,z为缺陷的深度;
步骤F2、当激光在刚铺放的表面检测时,被缺陷反射的体纵波将在横轴为扫查步数q,q=0,1,2,…,P-1,纵轴为时间t的B扫图上呈现双曲线的一支,将不同双曲线的各个值进行相加,所得结果即可反应C点的性质,获得垂直于铺放方向的刚刚铺放完成的缺陷的横向截面图,按照时间序列将横向截面图叠加获得空间图像,根据复合材料铺放过程中厚度逐渐累加的特性,在厚度方向上提取一定厚度的空间图像,即为最新铺层与下一层之间的缺陷图像;
步骤G、在大量试验的基础上获得铺放工艺参数包括铺放温度、压力和速度与铺放缺陷之间的关系,获得较好铺放质量时的工艺参数组合,选取一组铺放质量较好的工艺参数组合作为铺放工艺参数初始值,采用上述的激光无损检测方式获得在线铺放过程中最新铺层与下一层之间的缺陷信息,当出现缺陷超过铺放质量可接受的缺陷范围时,当前的工艺参数组合已经不再适合铺放;
根据铺放过程中压力传感器、温度传感器和速度传感器等获得的工艺参数值实时大小和铺放工艺参数组合与铺放质量之间的关系,对铺放工艺参数进行在线调整,重新获得适合当前铺放条件的工艺参数组合,减小铺放缺陷。
CN202010936787.8A 2020-09-08 2020-09-08 一种无损检测指导复合材料铺放工艺在线调整系统及方法 Active CN112172190B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010936787.8A CN112172190B (zh) 2020-09-08 2020-09-08 一种无损检测指导复合材料铺放工艺在线调整系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010936787.8A CN112172190B (zh) 2020-09-08 2020-09-08 一种无损检测指导复合材料铺放工艺在线调整系统及方法

Publications (2)

Publication Number Publication Date
CN112172190A true CN112172190A (zh) 2021-01-05
CN112172190B CN112172190B (zh) 2022-02-01

Family

ID=73925403

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010936787.8A Active CN112172190B (zh) 2020-09-08 2020-09-08 一种无损检测指导复合材料铺放工艺在线调整系统及方法

Country Status (1)

Country Link
CN (1) CN112172190B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113406003A (zh) * 2021-04-24 2021-09-17 南京理工大学 基于环形光束激光超声合成孔径聚焦成像装置及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102445416A (zh) * 2011-10-11 2012-05-09 东华大学 一种复合材料实时在线无损检测装置
CN106950280A (zh) * 2017-03-16 2017-07-14 东北大学 基于激光无损扫描的纤维增强复合材料参数辨识方法及装置
CN109374627A (zh) * 2018-08-30 2019-02-22 中国工程物理研究院机械制造工艺研究所 一种材料内部缺陷的激光超声透射时延检测方法
CN111257236A (zh) * 2020-01-07 2020-06-09 温州大学 双脉冲激光超声波检测装置及其检测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102445416A (zh) * 2011-10-11 2012-05-09 东华大学 一种复合材料实时在线无损检测装置
CN106950280A (zh) * 2017-03-16 2017-07-14 东北大学 基于激光无损扫描的纤维增强复合材料参数辨识方法及装置
CN109374627A (zh) * 2018-08-30 2019-02-22 中国工程物理研究院机械制造工艺研究所 一种材料内部缺陷的激光超声透射时延检测方法
CN111257236A (zh) * 2020-01-07 2020-06-09 温州大学 双脉冲激光超声波检测装置及其检测方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113406003A (zh) * 2021-04-24 2021-09-17 南京理工大学 基于环形光束激光超声合成孔径聚焦成像装置及方法

Also Published As

Publication number Publication date
CN112172190B (zh) 2022-02-01

Similar Documents

Publication Publication Date Title
US6078047A (en) Method and apparatus for terahertz tomographic imaging
CN112098526B (zh) 基于激光超声技术的增材制件近表面缺陷特征提取方法
Hou et al. Automatic multi-mode Lamb wave arrival time extraction for improved tomographic reconstruction
JP4854288B2 (ja) 超音波検査のための方法及びシステム
Selim et al. Defect reconstruction by non-destructive testing with laser induced ultrasonic detection
Baskaran et al. TOFD imaging: Ultrasonic TOFD flaw sizing and imaging in thin plates using embedded signal identification technique (ESIT)
US20130088724A1 (en) Method and apparatus for the inspection of sandwich structures using laser-induced resonant frequencies
CN112172190B (zh) 一种无损检测指导复合材料铺放工艺在线调整系统及方法
EP3940338A1 (en) Characterizing internal structures via ultrasound
Malo et al. Wave mode discrimination of coded ultrasonic guided waves using two-dimensional compressed pulse analysis
Safaeinili et al. Air-coupled ultrasonic estimation of viscoelastic stiffnesses in plates
US6176135B1 (en) System and method for laser-ultrasonic frequency control using optimal wavelength tuning
Antonio Jr et al. Ultrasonic imaging of concrete by synthetic aperture focusing technique based on hilbert-huang transform of time domain data
Barzegar et al. Pulse-echo and pitch-catch: Ultrasonic C-scan of adhesively bonded single lap-joints
Selim et al. Material defect reconstruction by non-destructive testing with laser induced ultrasonics
Yu et al. Multi-damage detection with embedded ultrasonic structural radar algorithm using piezoelectric wafer active sensors through advanced signal processing
Lian et al. Advances in Post-processing Technology for Laser Ultrasound Detection Signals: A Review
CN112149508B (zh) 基于剪切波特性融合用于焊缝缺陷回波信号特征提取稀疏表示的方法
JP2002168838A (ja) 内部欠陥検出方法、およびそれを用いて内部欠陥を検出する超音波顕微鏡
CN113953271B (zh) 一种激光清洗和无损检测系统及方法
Selim et al. Fully Non-Contact Hybrid NDT Inspection for 3D Defect Reconstruction Using an Improved SAFT Algorithm
Chiou et al. Modeling and processing of terahertz imaging in space shuttle external tank foam inspection
CN114487115B (zh) 一种基于Canny算子与超声平面波成像相结合的高分辨缺陷无损检测方法
CN116203128A (zh) 基于激光超声saft缺陷检测的优化方法及系统
Peng et al. Investigation of Signal Separation for Ultrasonic Nondestructive Evaluation

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