CN105571752B - 一种光弹性实验中超声波应力定量测量方法 - Google Patents
一种光弹性实验中超声波应力定量测量方法 Download PDFInfo
- Publication number
- CN105571752B CN105571752B CN201610010134.0A CN201610010134A CN105571752B CN 105571752 B CN105571752 B CN 105571752B CN 201610010134 A CN201610010134 A CN 201610010134A CN 105571752 B CN105571752 B CN 105571752B
- Authority
- CN
- China
- Prior art keywords
- photoelastic
- vibration
- image
- gray value
- 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.)
- Expired - Fee Related
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
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 (5)
1.一种光弹性实验中超声波应力定量测量方法，应用于由动态光弹性成像系统和激光测振仪构成的应力测量系统中，其特征在于，所述方法包括以下步骤：
利用动态光弹性成像系统获取超声纵波靠近样品表面处的光弹图像；利用激光测振仪测量样品的振动速度；以及利用计算机读取所述光弹图像中的灰度值；
通过调节所述动态光弹性成像系统的超声波强度，以获取样品表面处的多个光弹图像；并利用激光测振仪测量获取所述样品的多个振动速度；以及利用计算机读取所述多个光弹图像的灰度值；
根据所述样品的多个振动速度的平方与所述多个光弹图像的灰度值绘制灰度值-振速平方曲线；
读取光弹图像中任一点图像灰度值Bx，通过所述灰度值-振速平方曲线，确定相应的振动速度Ax；
根据所述振动速度Ax，计算光弹图像中所述任一点的超声波应力值；
所述根据所述振动速度Ax，计算光弹图像中所述任一点的超声波应力值，具体为：
通过公式σ1x＝(λ+2μ)Ax/cl，计算得出所述任一点沿超声波传播方向的超声波主应力值；
以及通过公式σ2x＝λAx/cl，计算得出所述任一点垂直于超声波传播方向的超声波主应力值；
其中λ、μ为样品材料的拉梅系数，cl为材料的纵波声速。
2.根据权利要求1所述的方法，其特征在于，所述样品表面处的光弹图像为距离样品表面4mm处的光弹图像。
3.根据权利要求1所述的方法，其特征在于，所述振动速度为所述样品底部最大振动速度。
4.根据权利要求1所述的方法，其特征在于，所述灰度值为所述测振仪测点所对应的光弹图像的最大灰度值。
5.根据权利要求1所述的方法，其特征在于，所述样品为K9玻璃制品。
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 CN105571752A (zh) | 2016-05-11 |
| CN105571752B true 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) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107884061B (zh) * | 2017-11-30 | 2020-04-14 | 中国科学院声学研究所 | 一种动态光弹超声成像方法和系统 |
| CN109579971B (zh) * | 2018-10-17 | 2020-08-04 | 中国科学院声学研究所 | 一种利用激光多普勒效应进行超声定量测量的系统及方法 |
| CN109959477B (zh) * | 2019-03-29 | 2021-02-19 | 华南理工大学 | 一种gis盆式绝缘子环氧试块内应力超声纵波检测方法及系统 |
Family Cites Families (4)
| 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 |
| CN100468044C (zh) * | 2005-03-17 | 2009-03-11 | 中国科学院半导体研究所 | 半导体材料残余应力的测试装置及方法 |
| DE102010027072A1 (de) * | 2010-07-13 | 2012-01-19 | Prüftechnik Dieter Busch AG | Verfahren und System zur Vorhersage von Fehlern an Komponenten rotierender Maschinen mit Thermographie |
| CZ2014500A3 (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 |
-
2016
- 2016-01-08 CN CN201610010134.0A patent/CN105571752B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN105571752A (zh) | 2016-05-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Haddadi et al. | Use of rigid-body motion for the investigation and estimation of the measurement errors related to digital image correlation technique | |
| Lopez-Crespo et al. | Study of a crack at a fastener hole by digital image correlation | |
| JP4568888B2 (ja) | 応力解析方法及び応力解析装置 | |
| CN105571752B (zh) | 一种光弹性实验中超声波应力定量测量方法 | |
| Xavier et al. | Direct evaluation of cohesive law in mode I of Pinus pinaster by digital image correlation | |
| CN103076347B (zh) | 基于原位x射线断层照相的脆性材料力学损伤的测量方法 | |
| PI et al. | Measure strain distribution using digital image correlation (DIC) for tensile tests | |
| Ayatollahi et al. | Digital image correlation method for calculating coefficients of Williams expansion in compact tension specimen | |
| Kumpenza et al. | Measuring Poisson’s ratio: mechanical characterization of spruce wood by means of non-contact optical gauging techniques | |
| Brabec et al. | Utilization of digital image correlation in determining of both longitudinal shear moduli of wood at single torsion test | |
| KR101447833B1 (ko) | 응력 확대 계수 측정 방법 | |
| CN111307347A (zh) | 一种试件表面主应力测试装置及方法 | |
| Gharavi et al. | End effect on determining shear modulus of timber beams in torsion tests | |
| Molˈkov | Application of the method of digital image correlation to the construction of stress–strain diagrams | |
| Yin et al. | High-temperature fracture mechanics parameter measurement and yielding zone analysis of superalloy GH4169 based on single-lens 3D digital image correlation | |
| Dahl et al. | Planar strain measurements on wood specimens | |
| CN109632487A (zh) | 一种室内岩石在载荷作用下内部应变场的定性测试方法 | |
| Nežerka et al. | Use of open source DIC tools for analysis of multiple cracking in fiber-reinforced concrete | |
| Gonzales et al. | Evaluation of the effect of water content on the stress optical coefficient in PMMA | |
| JP2008139273A (ja) | ひずみの測定方法及び測定システム | |
| CN211401505U (zh) | 一种试件表面主应力测试装置 | |
| CN106370727B (zh) | 基于应力波波场检测波速、材料状态及损伤程度的方法 | |
| CN101929847A (zh) | 基于碳纳米管为传感媒介的应变分量无损检测技术 | |
| Holzer et al. | Development of the bulge test equipment for measuring mechanical properties of thin films | |
| CN117571506B (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 | ||
| 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 |