CN106338445B - 一种基于频闪照明的疲劳裂纹扩展试验图像在线采集方法 - Google Patents
一种基于频闪照明的疲劳裂纹扩展试验图像在线采集方法 Download PDFInfo
- Publication number
- CN106338445B CN106338445B CN201610724084.2A CN201610724084A CN106338445B CN 106338445 B CN106338445 B CN 106338445B CN 201610724084 A CN201610724084 A CN 201610724084A CN 106338445 B CN106338445 B CN 106338445B
- Authority
- CN
- China
- Prior art keywords
- frequency
- image
- stroboscopic
- test
- crack propagation
- 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
Links
- 238000012360 testing method Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000005286 illumination Methods 0.000 title claims abstract description 19
- 238000009825 accumulation Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 208000037656 Respiratory Sounds Diseases 0.000 claims description 4
- 230000001960 triggered effect Effects 0.000 claims description 4
- 238000003360 curve fit method Methods 0.000 claims description 3
- 238000009661 fatigue test Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 240000002853 Nelumbo nucifera Species 0.000 claims description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 2
- 238000004088 simulation Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 230000005284 excitation Effects 0.000 abstract 2
- 238000009826 distribution Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004801 process automation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000010421 standard material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/32—Investigating strength properties of solid materials by application of mechanical stress by applying repeated or pulsating forces
- G01N3/38—Investigating strength properties of solid materials by application of mechanical stress by applying repeated or pulsating forces generated by electromagnetic means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
- G01N21/88—Investigating the presence of flaws or contamination
- G01N21/8851—Scan 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
- G01N21/88—Investigating the presence of flaws or contamination
- G01N21/8851—Scan 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/8887—Scan 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0001—Type of application of the stress
- G01N2203/0005—Repeated or cyclic
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0062—Crack or flaws
- G01N2203/0066—Propagation of crack
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- Pathology (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Signal Processing (AREA)
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Electromagnetism (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种基于频闪照明的疲劳裂纹扩展试验图像在线采集方法,包括以下步骤:1)将散斑试件装夹到试验机上,通过扫频试验获得共振频率;2)在激振信号的零相位时触发频闪光源和CCD快门,设定频闪光源和CCD快门的触发频率,采集一系列散斑试件的图像;3)采用快速数字图像相关方法计算裂纹散斑图像的裂纹尖端弹性区域的应变场,提取应变场中同一特征点的应变值进行最小二乘正弦曲线快速拟合。4)计算拟合出的应变值正弦曲线与试验机激振信号之间的相位差,并将此相位差反馈给频闪控制器调节频闪光源与CCD快门的触发时间;再进行二次图像采集,采集到载荷最大点与最小点位置时的图像。本发明易于实现、成本较低、且可实现实时在线处理。
Description
技术领域
本发明涉及高频振动图像采集技术领域,特别涉及一种基于频闪照明的疲劳裂纹扩展试验图像在线采集方法。
背景技术
高频疲劳裂纹扩展试验是指标准材料试件在高频率的交变载荷作用,探究裂纹演化规律的试验。为了实现整个过程自动化和智能化,通过摄像头采集喷了散斑的试件图像,传输至计算机,通过数字图像相关方法(DIC)获得裂纹区域附近的微观变形场,研究裂纹附近变形场在单个交变载荷周期内的变化规律以及在整个裂纹扩展过程中的变化规律。高频疲劳裂纹扩展试验中的交变载荷的疲劳一般在80Hz~250Hz,要想采集到一个载荷周期内的图像,尤其是要采集到最大和最小载荷位置时刻的图像,具有较高的难度。一般可以采用使用高速摄像的方式,在一秒钟内采集上千张图像,但是高速摄像机一般价格不菲,且每秒钟要处理上千张图像,这无疑对计算机的性能提出了更高的要求,较难实现实时在线。
发明内容
为了克服已有高速摄影在高频疲劳裂纹扩展试验图像采集过程中成本较高,无法实现在线的不足,本发明提供一种易于实现、成本较低、且可实现实时在线处理的基于频闪照明的疲劳裂纹扩展试验图像在线采集方法。
本发明解决其技术问题所采用的技术方案是:
一种基于频闪照明的疲劳裂纹扩展试验图像在线采集方法,包括以下步骤:
1)疲劳裂纹扩展试验中采用电磁谐振式高频疲劳试验机,将散斑试件装夹到试验机上,通过扫频试验获得共振频率,该共振频率也是施加到试件上交变载荷的频率;
2)首先在激振信号的零相位时触发频闪光源和CCD快门,按照公式(1)设定频闪光源和CCD快门的触发频率,采集一系列散斑试件的图像;
其中,f采为采样得到的正弦曲线频率,f为待采集载荷信号的频率,N为两个采样点之间相隔的整数个波形,n为完成一个完整波形采样的采样点数;
3)采用快速数字图像相关方法计算裂纹散斑图像的裂纹尖端弹性区域的应变场,提取应变场中同一特征点的应变值进行最小二乘正弦曲线快速拟合。
4)计算拟合出的应变值正弦曲线与试验机激振信号之间的相位差,并将此相位差反馈给频闪控制器调节频闪光源与CCD快门的触发时间;接下来再进行二次图像采集,采集到载荷最大点与最小点位置时的图像。
进一步,所述步骤3)中,快速拟合正弦曲线过程使用三参数正弦波拟合法并使用和角公式将非线性问题转化为线性为题求解。
再进一步,所述步骤3)中,裂纹尖端弹性区域内特征点的正弦曲线拟合方法如下:
假设理想的正弦波信号为:
y(t)=Asin(ωt+φ)+D (2)
其中,y是裂纹尖端弹性区域某一特征点的应变值,A为振幅,φ为初相位,D为y平均值,ω=2πf为前文中扫频试验获得的共振角频率,通过正弦和角变换可变为:
y(t)=Acosφsin(ωt)+Asinφcos(ωt)+D (3)
三参数正弦波拟合过程即为输入信号的频率已知,求取A、φ、D使误差平方和函数最小:
构造矩阵:
则(4)式用矩阵表示为:
ξ=ξ(n)=(y-ψx0)T(y-ψx0) (5)
当(4)式最小时得x0的最小二乘解为:
本发明的有益效果主要表现在:可以精确采集到在高频疲劳试验过程中一个载荷周期最大值以及最小值时刻的图像,且成本较低,算法简单,易于实现实时在线。
附图说明
图1是基于频闪照明的疲劳裂纹扩展试验图像在线采集方法系统总体组成。
图2是基于频闪照明的疲劳裂纹扩展试验图像在线采集方法流程图。
图3是基于频闪照明的疲劳裂纹扩展试验图像在线采集方法中二级采样频闪照明图像采集原理,实线表示载荷曲线,黑点表示频闪图像采集点,虚线表示频闪图像采集点拟合曲线,其中,(a)表示f=100Hz,N=2,n=8,(b)表示f=100Hz,N=5,n=10,(c)表示f=80Hz,N=2,n=25。
图4是基于频闪照明的疲劳裂纹扩展试验图像在线采集方法中通过数字图像相关方法计算得到的裂纹区域应变场分布图,其中,(a)~(h)分别表示一个载荷循环周期内设定8个采样点时对应时刻的裂纹尖端区域的应变场。
图5是基于频闪照明的疲劳裂纹扩展试验图像在线采集方法中应变值正弦曲线拟合以及与激振信号间的相位差。
图6是基于频闪照明的疲劳裂纹扩展试验图像在线采集方法中激振信号、试件载荷、频闪光源、相机同步控制时序关系。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图1~图6,一种基于频闪照明的疲劳裂纹扩展试验图像在线采集方法,实现该方法的系统包括电磁谐振式疲劳试验机,CCD和镜头,光源,频闪控制器,图像采集卡,CCD控制器,同步控制器,数字载荷控制器和计算机,所述采集方法包括如下步骤:
首先将喷了散斑的试件装夹到试验机上,电磁谐振式疲劳试验机是一个振动系统,通过扫频试验得到系统的共振频率,一般为80~250Hz,向激振器发出该共振频率的激振信号,开始进行疲劳裂纹扩展试验。
接下来确定频闪采样频率,为了采集高频振动下的散斑图像,本文使用了二级采样频闪照明方法,如图3所示,具体采样频率公式如下:
其中,f采为采样得到的正弦曲线频率,f为待采集载荷信号的频率,N为两个采样点之间相隔的整数个波形,n为完成一个完整波形采样的采样点数。图3中分别显示了不同采样频率不同采样点数情况时的采样结果。本发明实施例中设定N=10,n=8,根据上述公式计算出来的采样频率触发频闪光源和CCD快门时间。为了在图像采集时将试件表面的散斑特征清晰地凸现出来,本实施例中采用LED环形频闪光源低角度照明的方式。相机选用SONY公司生产的XC-HR70工业相机,它是一款具有高分辨率、高帧速率顺序扫描的黑白CCD摄像机,并搭配Myutron公司生产的12mm定焦工业镜头LM12NC3。图像采集卡选用美国NI公司的PCI-1410图像采集卡。接下来采集试件的散斑图像,并将图像传输至计算机,采用数字图像相关方法计算裂纹尖端区域应变场,图4为一个载荷循环周期内裂纹尖端区域应变场分布情况。根据线弹性断裂力学理论,在裂纹尖端弹性区域内的应变值和施加到试件上的载荷服从同样的正弦变化规律。
为了快速拟合正弦曲线,使用了三参数正弦波拟合法并使用和角公式将非线性问题转化为线性为题求解,大大的降低了计算量从而降低了计算时间。其基本原理为假设理想的正弦波信号为:
y(t)=Asin(ωt+φ)+D (2)
其中,y是裂纹尖端弹性区域某一特征点的应变值,A为振幅,φ为初相位,D为y平均值,ω=2πf为前文中扫频试验获得的共振角频率。通过正弦和角变换可变为:
y(t)=Acosφsin(ωt)+Asinφcos(ωt)+D (3)
三参数正弦波拟合过程即为输入信号的频率已知,求取A、φ、D使误差平方和函数最小:
构造矩阵:
则(4)式可用矩阵表示为:
ξ=ξ(n)=(y-ψx0)T(y-ψx0) (5)
当(4)式最小时可得x0的最小二乘解为:
这样就得到了拟合正弦曲线的函数式,如图5显示了应变值正弦曲线拟合结果,并将其与激振信号比较求出两者间的相位差,将其该相位差反馈给频闪光源和CCD快门进行相位延时,图6是激振信号、试件载荷、频闪光源、相机同步控制时序关系。接下来再进行第二次试件散斑图像采集,便可采集到最大载荷和最小载荷处等关键载荷点处的图像。
最后说明的是,以上实施例仅仅是对于本发明专利精神作举例说明。本发明专利所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方法替代,但并不会偏离本发明专利的精神或者超越所附权利要求书所定义的范围。
Claims (3)
1.一种基于频闪照明的疲劳裂纹扩展试验图像在线采集方法,其特征在于:包括以下步骤:
1)疲劳裂纹扩展试验中采用电磁谐振式高频疲劳试验机,将散斑试件装夹到试验机上,通过扫频试验获得共振频率,该共振频率也是施加到试件上交变载荷的频率;
2)首先在激振信号的零相位时触发频闪光源和CCD快门,按照公式(1)设定频闪光源和CCD快门的触发频率,采集一系列散斑试件的图像;
其中,f采为采样得到的正弦曲线频率,f为待采集载荷信号的频率,N为两个采样点之间相隔的整数个波形,n为完成一个完整波形采样的采样点数;
3)采用快速数字图像相关方法计算裂纹散斑图像的裂纹尖端弹性区域的应变场,提取应变场中同一特征点的应变值进行最小二乘正弦曲线快速拟合;
4)计算拟合出的应变值正弦曲线与试验机激振信号之间的相位差,并将此相位差反馈给频闪控制器调节频闪光源与CCD快门的触发时间;接下来再进行二次图像采集,采集到载荷最大点与最小点位置时的图像。
2.如权利要求1所述的一种基于频闪照明的疲劳裂纹扩展试验图像在线采集方法,其特征在于:所述步骤3)中,快速拟合正弦曲线过程使用三参数正弦波拟合法并使用和角公式将非线性问题转化为线性问题求解。
3.如权利要求2所述的一种基于频闪照明的疲劳裂纹扩展试验图像在线采集方法,其特征在于:所述步骤3)中,裂纹尖端弹性区域内特征点的正弦曲线拟合方法如下:
假设理想的正弦波信号为:
y(t)=Asin(ωt+φ)+D (2)
其中,y是裂纹尖端弹性区域某一特征点的应变值,A为振幅,φ为初相位,D为y平均值,ω=2πf为所述扫频试验获得的共振角频率,通过正弦和角变换可变为:
y(t)=Acosφsin(ωt)+Asinφcos(ωt)+D (3)
三参数正弦波拟合过程即为输入信号的频率已知,求取A、φ、D使误差平方和函数最小:
构造矩阵:;
则(4)式用矩阵表示为:
当(4)式最小时得x0的最小二乘解为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610724084.2A CN106338445B (zh) | 2016-08-25 | 2016-08-25 | 一种基于频闪照明的疲劳裂纹扩展试验图像在线采集方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610724084.2A CN106338445B (zh) | 2016-08-25 | 2016-08-25 | 一种基于频闪照明的疲劳裂纹扩展试验图像在线采集方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106338445A CN106338445A (zh) | 2017-01-18 |
CN106338445B true CN106338445B (zh) | 2019-07-05 |
Family
ID=57824584
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610724084.2A Active CN106338445B (zh) | 2016-08-25 | 2016-08-25 | 一种基于频闪照明的疲劳裂纹扩展试验图像在线采集方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106338445B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108827804B (zh) * | 2018-07-12 | 2021-04-06 | 浙江工业大学 | 一种谐振式疲劳试验机动态载荷误差在线补偿方法 |
CN110987617A (zh) * | 2019-12-11 | 2020-04-10 | 上海交通大学 | 拉伸应力应变曲线在缩颈阶段时由系统变频的dic测量方法 |
CN115511881B (zh) * | 2022-11-08 | 2023-04-25 | 南京航空航天大学 | 一种数字图像相关和数字体相关中的相关性调谐法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1869640A (zh) * | 2006-05-24 | 2006-11-29 | 浙江大学 | 一种疲劳裂纹扩展测试的方法 |
CN102692188A (zh) * | 2012-05-08 | 2012-09-26 | 浙江工业大学 | 机械视觉疲劳裂纹扩展试验裂纹长度动态测量方法 |
CN105699218A (zh) * | 2016-01-26 | 2016-06-22 | 清华大学 | 一种在线测量有机玻璃疲劳裂纹扩展的方法 |
-
2016
- 2016-08-25 CN CN201610724084.2A patent/CN106338445B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1869640A (zh) * | 2006-05-24 | 2006-11-29 | 浙江大学 | 一种疲劳裂纹扩展测试的方法 |
CN102692188A (zh) * | 2012-05-08 | 2012-09-26 | 浙江工业大学 | 机械视觉疲劳裂纹扩展试验裂纹长度动态测量方法 |
CN105699218A (zh) * | 2016-01-26 | 2016-06-22 | 清华大学 | 一种在线测量有机玻璃疲劳裂纹扩展的方法 |
Non-Patent Citations (1)
Title |
---|
孔径旋转频闪散斑照相法测量物体的固有频率;陈炳泉;《光学学报》;20021130;第22卷(第11期);1358-1361 * |
Also Published As
Publication number | Publication date |
---|---|
CN106338445A (zh) | 2017-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106338445B (zh) | 一种基于频闪照明的疲劳裂纹扩展试验图像在线采集方法 | |
CN108519225B (zh) | 一种航空发动机叶片高温振动疲劳特性的检测装置及方法 | |
CN107228748A (zh) | 基于非接触式测量的卫星天线结构振动测量装置与方法 | |
CN110310237A (zh) | 去除图像摩尔纹、显示面板子像素点亮度测量、Mura缺陷修复的方法及系统 | |
EP2031374A3 (en) | Apparatus and method for obtaining information related to terahertz waves | |
WO2004010380A3 (en) | Measuring 3d deformations of an object by comparing focusing conditions for sharp capturing of said object before and after deformation | |
CN103760025A (zh) | 引伸计及其测量方法 | |
CN102359966A (zh) | 一种混凝土表面微小裂缝的定位系统 | |
CN105842062A (zh) | 一种裂纹扩展实时监测装置及方法 | |
CN101354240A (zh) | 基于显微光学的微纳米级纤维高精度测量系统 | |
CN104020086A (zh) | 一种基于数字图像处理技术的液雾特性快速测量方法 | |
Rizzuto et al. | Characterization of a digital image correlation system for dynamic strain measurements of small biological tissues | |
CN103324888B (zh) | 基于家族样本的病毒特征自动提取方法及系统 | |
CN104237031A (zh) | 基于数字图像的分离式霍普金森压杆实验的同步测量方法 | |
TWI413799B (zh) | The method of automatically looking for the focus position of the optical microscope | |
CN109271845A (zh) | 基于计算机视觉的人体动作分析与评价方法 | |
CN108093173A (zh) | 光耦合式图像采集时刻标定的系统和方法 | |
CN205301075U (zh) | 一种测定冲击试样断口纤维断面率系统 | |
CN104434139A (zh) | 一种检验科用多功能智能采血系统 | |
Tessari et al. | Advanced sensing in geotechnical centrifuge models | |
Guangmang et al. | Micro-image definition evaluation using multi-scale decomposition and gradient absolute value | |
CN101226641A (zh) | 凝视型相机运动点目标定位方法与系统 | |
WO2017113152A1 (zh) | 基于多组光源的心率测量方法及实现该方法的穿戴式设备 | |
CN109444150A (zh) | 非接触式裂纹测量方法及其装置 | |
CN104799811B (zh) | 一种用于眼底照相的照明光控制电路 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |