CN105571752A - 一种光弹性实验中超声波应力定量测量方法 - Google Patents

一种光弹性实验中超声波应力定量测量方法 Download PDF

Info

Publication number
CN105571752A
CN105571752A CN201610010134.0A CN201610010134A CN105571752A CN 105571752 A CN105571752 A CN 105571752A CN 201610010134 A CN201610010134 A CN 201610010134A CN 105571752 A CN105571752 A CN 105571752A
Authority
CN
China
Prior art keywords
photoelastic
gray
vibration
image
sample
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
CN201610010134.0A
Other languages
English (en)
Other versions
CN105571752B (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.)
Institute of Acoustics CAS
Original Assignee
Institute of Acoustics CAS
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 Institute of Acoustics CAS filed Critical Institute of Acoustics CAS
Priority to CN201610010134.0A priority Critical patent/CN105571752B/zh
Publication of CN105571752A publication Critical patent/CN105571752A/zh
Application granted granted Critical
Publication of CN105571752B publication Critical patent/CN105571752B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00—Measuring force or stress, in general
    • G01L1/25—Measuring force or stress, in general using wave or particle radiation, e.g. X-rays, microwaves, neutrons
    • G01L1/255—Measuring force or stress, in general using wave or particle radiation, e.g. X-rays, microwaves, neutrons using acoustic waves, or acoustic emission

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Health & Medical Sciences (AREA)
  • Toxicology (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明涉及一种光弹性实验中超声波应力定量测量方法，该方法包括以下步骤：利用动态光弹性成像系统获取超声纵波靠近样品表面处的光弹图像；利用激光测振仪测量样品的振动速度；以及利用计算机读取光弹图像中的灰度值；调节动态光弹性成像系统的超声波强度，重复上述步骤以获取样品多个振动速度和多个光弹图像的灰度值；根据样品的多个振动速度的平方与多个光弹图像的灰度值绘制灰度值-振速平方曲线；读取光弹图像中任一点图像灰度值Bx，通过灰度值-振速平方曲线，确定相应的振动速度Ax，可计算得出光弹图像中该点的超声波应力值。本发明提供的超声波应力定量测量方法，可快速精确定量测量光弹性实验中超声波的应力。

Description

一种光弹性实验中超声波应力定量测量方法
技术领域
本发明涉及一种光弹性实验中应力测量方法，特别是涉及一种光弹性实验中超声波应力定量测量方法。
背景技术
光弹性实验是以暂时双折射现象为物理基础来研究弹性力学问题的一种实验应力分析方法。近年来，光弹性实验在汽车制造、航天工业以及土木工程等领域得到了广泛应用。光弹性实验能够直接观察超声波的传播和散射，所以也应用到了固体超声波散射问题的研究中。
目前光弹性实验中应力测量的方法主要有：拉伸试件补偿法、石英楔块补偿法、检偏镜补偿法和光度计补偿法。
拉伸试件补偿法，该方法的拉力试件与被测点的主应力方向平行，并使拉力试件与被测点重叠投影在成像屏幕上，给拉力试件加力，直至测点在屏幕变黑为止，再读取拉力值，从而获得条纹值。但是，相对于光弹实验中施加的其他类型的外力，超声波的应力较小，在加力使测点变黑的过程中易受自身系统的衍射光以及环境光污染干扰，难以精确定量。此外，该方法操作繁琐，需要多次成像以确定测点是否变黑。
石英楔块法是采用两块直角三角形石英楔块调节光程差，直至测点在屏幕变黑为止，再读取拉力值，从而获得条纹值。与拉伸试件补偿法类似，也面临难以精确定量和多次成像的问题。
检偏镜补偿法操作较复杂，不够简便。而对于光度计补偿法，则面临难以准确调整偏正片的问题。
发明内容
本发明的目的是针对现有技术中光弹性实验中超声波应力测量法方法在测量速度、精确度以及操作便利性方面存在的缺陷，提供一种光弹性实验中超声波应力测量方法，以快速精确测量光弹性实验中的超声波应力。
为实现上述目的，本发明提供了一种光弹性实验中超声波应力定量测量方法，该方法应用于由动态光弹性成像系统和激光测振仪构成的应力测量系统中，该方法包括以下步骤：
利用动态光弹法成像系统获取超声纵波靠近样品表面处的光弹图像；利用激光测振仪测量样品的振动速度；以及利用计算机读取光弹图像中的灰度值。通过调节动态光弹性成像系统的超声波强度，以获取样品表面处的多个光弹图像；并利用激光测振仪测量获取样品的多个振动速度；以及利用计算机读取多个光弹图像的灰度值。根据样品的多个振动速度的平方与多个光弹图像的灰度值绘制灰度值-振速平方曲线。
读取光弹图像中任一点图像灰度值Bx，通过灰度值-振速平方曲线，确定相应的振动速度Ax。根据该振动速度Ax，计算该任一点的超声波应力值，具体为：通过公式σ1x＝(λ+2μ)Ax/cl(其中λ、μ为样品材料的拉梅系数，cl为材料的纵波声速)，计算得出该任一点沿超声波传播方向的超声波主应力值；以及通过公式σ2x＝λAx/cl，计算得出该任一点垂直于超声波传播方向的超声波主应力值。
优选地，上述样品表面处的光弹图像为距离样品表面4mm处的光弹图像。
优选地，上述振动速度为样品底部最大振动速度。
优选地，上述灰度值为测振仪测点所对应的光弹图像的最大灰度值。
优选地，上述样品为K9玻璃制品。
本发明提供的光弹性实验中超声波应力定量测量方法绘制了灰度值-振速平方对应曲线，根据光弹图像上任一点的灰度值，通过灰度值-振速平方对应曲线得出该点的振动速度，然后，根据通过公式σ1x＝(λ+2μ)Ax/cl和σ2x＝λAx/cl计算出该点的超声波应力值，具有快速、精确和操作便利的优点。
附图说明
图1为本发明实施例提供的动态光弹性成像系统示意图；
图2为本发明实施例提供的激光测振仪原理图；
图3本发明实施例提供的灰度值-振速平方对应曲线。
具体实施方式
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
绘制灰度值-振速平方对应曲线。
1、选用K9玻璃制成的样品。
2、利用图1所示的动态光弹性成像系统获取超声纵波距离样品表面4mm处的光弹图像。
3、利用图2所示的激光测振仪测量样品底部最大振动速度A1。
4、利用计算机读取测振仪测点所对应的光弹图像的最大灰度值B1。
5、调节动态光弹性成像系统超声波强度，重复步骤2-4，可随机设定重复次数，本实施例重复次数为12次，获得样品底部最大振动速度A2、A3…A12和振仪测点所对应的光弹图像的最大灰度值B2、B3…B12，并计算振速A的平方，所得数据如表1所示。
表1
振速A的平方 0.4624 0.8464 1.1664 1.5376 2.1904 2.56 2.6896 3.5344 6.1504 6.9696 9.9856 16
灰度值B 14.0079 15.5635 16.3333 25.7024 24.119 35.9524 27.6627 35.873 73.5476 67.246 112.0881 118.6333
6、振动速度A1-A12的平方与灰度值B1-B12线性拟合，得到图3所示的灰度值-振速平方对应曲线。
利用灰度值-振速平方曲线计算得出光弹图像中任一点的超声波应力值。
1、利用计算机读取光弹图像中任一点图像灰度值，在本实施例中，该任一点的灰度值为30.
2、通过图3所示的灰度值-振速平方曲线获得灰度值30相应的振动速度平方2.25cm2/s2，则振速为0.015m/s。
3、通过公式σ1x＝(λ+2μ)A/cl＝2.241×105Pa(其中λ、μ为样品材料的拉梅系数，cl为材料的纵波声速，K9玻璃样品的λ＝2.2825×1010、μ＝3.308×1010、cl＝5957m/s)，即可得到该点沿超声波传播方向的主应力值；通过公式σ2x＝λA/cl＝5.747×104Pa，即可得到该点垂直于超声波传播方向的超声波主应力值。
本实施例提供的光弹性实验中超声波应力定量测量方法首先绘制了灰度值-振速平方对应曲线，其次，读取光弹图像上任一点的灰度值，通过灰度值-振速平方对应曲线得出该点的振动速度，再次，根据通过公式σ1x＝(λ+2μ)Ax/cl和σ2x＝λAx/cl计算出该点的超声波应力值，具有快速、精确和操作便利的优点。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.一种光弹性实验中超声波应力定量测量方法，应用于由动态光弹性成像系统和激光测振仪构成的应力测量系统中，其特征在于，所述方法包括以下步骤：
利用动态光弹性成像系统获取超声纵波靠近样品表面处的光弹图像；利用激光测振仪测量样品的振动速度；以及利用计算机读取所述光弹图像中的灰度值；
通过调节所述动态光弹性成像系统的超声波强度，以获取样品表面处的多个光弹图像；并利用激光测振仪测量获取所述样品的多个振动速度；以及利用计算机读取所述多个光弹图像的灰度值；
根据所述样品的多个振动速度的平方与所述多个光弹图像的灰度值绘制灰度值-振速平方曲线；
读取光弹图像中任一点图像灰度值Bx，通过所述灰度值-振速平方曲线，确定相应的振动速度Ax；
根据所述振动速度Ax，计算光弹图像中所述任一点的超声波应力值。
2.根据权利要求1所述的方法，其特征在于，所述根据所述振动速度Ax，计算光弹图像中所述任一点的超声波应力值，具体为：
通过公式σ1x＝(λ+2μ)Ax/cl(其中λ、μ为样品材料的拉梅系数，cl为材料的纵波声速)，计算得出所述任一点沿超声波传播方向的超声波主应力值；
以及通过公式σ2x＝λAx/cl，计算得出所述任一点垂直于超声波传播方向的超声波主应力值。
3.根据权利要求1所述的方法，其特征在于，所述样品表面处的光弹图像为距离样品表面4mm处的光弹图像。
4.根据权利要求1所述的方法，其特征在于，所述振动速度为所述样品底部最大振动速度。
5.根据权利要求1所述的方法，其特征在于，所述灰度值为所述测振仪测点所对应的光弹图像的最大灰度值。
6.根据权利要求1所述的方法，其特征在于，所述样品为K9玻璃制品。
CN201610010134.0A 2016-01-08 2016-01-08 一种光弹性实验中超声波应力定量测量方法 Expired - Fee Related CN105571752B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610010134.0A CN105571752B (zh) 2016-01-08 2016-01-08 一种光弹性实验中超声波应力定量测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610010134.0A CN105571752B (zh) 2016-01-08 2016-01-08 一种光弹性实验中超声波应力定量测量方法

Publications (2)

Publication Number Publication Date
CN105571752A true CN105571752A (zh) 2016-05-11
CN105571752B CN105571752B (zh) 2018-12-21

Family

ID=55882141

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610010134.0A Expired - Fee Related CN105571752B (zh) 2016-01-08 2016-01-08 一种光弹性实验中超声波应力定量测量方法

Country Status (1)

Country Link
CN (1) CN105571752B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107884061A (zh) * 2017-11-30 2018-04-06 中国科学院声学研究所 一种动态光弹超声成像方法和系统
CN109579971A (zh) * 2018-10-17 2019-04-05 中国科学院声学研究所 一种利用激光多普勒效应进行超声定量测量的系统及方法
CN109959477A (zh) * 2019-03-29 2019-07-02 华南理工大学 一种gis盆式绝缘子环氧试块内应力超声纵波检测方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4033182A (en) * 1976-06-28 1977-07-05 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Method for measuring biaxial stress in a body subjected to stress inducing loads
CN1834623A (zh) * 2005-03-17 2006-09-20 中国科学院半导体研究所 半导体材料残余应力的测试装置及方法
CN103140754A (zh) * 2010-07-13 2013-06-05 普乐福尼克·迪特·布什股份公司 用于利用热成像预测旋转机器的零件上的故障的方法和系统
CZ305286B6 (cs) * 2014-07-17 2015-07-15 České Vysoké Učení Technické V Praze - Fakulta Jaderná A Fyzikálně Inženýrská Soustava pro měření zbytkových napětí v polykrystalických materiálech metodou rentgenové difrakce

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4033182A (en) * 1976-06-28 1977-07-05 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Method for measuring biaxial stress in a body subjected to stress inducing loads
CN1834623A (zh) * 2005-03-17 2006-09-20 中国科学院半导体研究所 半导体材料残余应力的测试装置及方法
CN103140754A (zh) * 2010-07-13 2013-06-05 普乐福尼克·迪特·布什股份公司 用于利用热成像预测旋转机器的零件上的故障的方法和系统
CZ305286B6 (cs) * 2014-07-17 2015-07-15 České Vysoké Učení Technické V Praze - Fakulta Jaderná A Fyzikálně Inženýrská Soustava pro měření zbytkových napětí v polykrystalických materiálech metodou rentgenové difrakce

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
毕谦等: "动态光弹性法在穿爆效应及边坡稳定中的初步应用", 《云南冶金》 *
金士杰等: "光弹法测量超声换能器声场", 《应用声学》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107884061A (zh) * 2017-11-30 2018-04-06 中国科学院声学研究所 一种动态光弹超声成像方法和系统
CN107884061B (zh) * 2017-11-30 2020-04-14 中国科学院声学研究所 一种动态光弹超声成像方法和系统
CN109579971A (zh) * 2018-10-17 2019-04-05 中国科学院声学研究所 一种利用激光多普勒效应进行超声定量测量的系统及方法
CN109579971B (zh) * 2018-10-17 2020-08-04 中国科学院声学研究所 一种利用激光多普勒效应进行超声定量测量的系统及方法
CN109959477A (zh) * 2019-03-29 2019-07-02 华南理工大学 一种gis盆式绝缘子环氧试块内应力超声纵波检测方法及系统
CN109959477B (zh) * 2019-03-29 2021-02-19 华南理工大学 一种gis盆式绝缘子环氧试块内应力超声纵波检测方法及系统

Also Published As

Publication number Publication date
CN105571752B (zh) 2018-12-21

Similar Documents

Publication Publication Date Title
Lopez-Crespo et al. Study of a crack at a fastener hole by digital image correlation
JP2011053157A (ja) ひずみ計測方法、ひずみ計測装置およびプログラム
Huang et al. Development of a novel laser-based measuring system for the thread profile of ballscrew
CN105571752A (zh) 一种光弹性实验中超声波应力定量测量方法
Xavier et al. Measuring displacement fields by cross-correlation and a differential technique: experimental validation
Yue et al. Applications of digital image correlation (DIC) and the strain gage method for measuring dynamic mode I fracture parameters of the white marble specimen
Kumpenza et al. Measuring Poisson’s ratio: mechanical characterization of spruce wood by means of non-contact optical gauging techniques
Yu et al. Modeling of systematic errors in stereo-digital image correlation due to camera self-heating
Mohamed et al. Photogrammetric evaluation of shear modulus of glulam timber using torsion test method and dual stereo vision system
Molˈkov Application of the method of digital image correlation to the construction of stress–strain diagrams
Philippon et al. A device enhancement for the dry sliding friction coefficient measurement between steel 1080 and vascomax with respect to surface roughness changes
He Experimental investigation on large deformation of silicon rubber
Mousavi et al. Non-equilibrium split Hopkinson pressure bar procedure for non-parametric identification of complex modulus
Gharavi et al. End effect on determining shear modulus of timber beams in torsion tests
Nežerka et al. Use of open source DIC tools for analysis of multiple cracking in fiber-reinforced concrete
Zhang et al. Mixed mode fracture study of PMMA using digital gradient sensing method
Othman et al. Testing aluminum alloy from quasi-static to dynamic strain-rates with a modified split Hopkinson bar method
Li et al. Measurements of velocity distributions in the deformation zone in cold rolling by a scanning LDV
Ivanyts’ kyi et al. The choice of spacing in measuring displacements for the evaluation of strains by the method of optical-digital image correlation
RU2535645C1 (ru) Способ определения характеристик изгибной жесткости протяженных объектов с помощью кривизномера
Richter‐Trummer et al. Methodology for in situ stress intensity factor determination on cracked structures by digital image correlation
Ju et al. Determining 2D notch SIFs by the image-correlation method
Holzer et al. Development of the bulge test equipment for measuring mechanical properties of thin films
Liu et al. Parallax correction for a digital array high-speed imaging system for dynamic deformation measurements
KR101033031B1 (ko) 변형률 측정 장치

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20181221

Termination date: 20200108