CN109238553B - 一种针对水中柱状冲击波的压强分布测量方法 - Google Patents
一种针对水中柱状冲击波的压强分布测量方法 Download PDFInfo
- Publication number
- CN109238553B CN109238553B CN201811212878.6A CN201811212878A CN109238553B CN 109238553 B CN109238553 B CN 109238553B CN 201811212878 A CN201811212878 A CN 201811212878A CN 109238553 B CN109238553 B CN 109238553B
- Authority
- CN
- China
- Prior art keywords
- formula
- columnar
- shock waves
- shock wave
- light
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L11/00—Measuring steady or quasi-steady pressure of a fluid or a fluent solid material by means not provided for in group G01L7/00 or G01L9/00
- G01L11/02—Measuring steady or quasi-steady pressure of a fluid or a fluent solid material by means not provided for in group G01L7/00 or G01L9/00 by optical means
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Fluid Pressure (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
本发明公开了一种针对水中柱状冲击波的压强分布测量方法,柱状冲击波为二维冲击波,采用摄像装置获取光源发出的光线经过散斑背景、无柱状冲击波的二维水槽的第一图像集;采用摄像装置获取光源发出的光线经过散斑背景、加载柱状冲击波的二维水槽的第二图像集;基于第一公式,通过第一图像集和第二图像集的对比确定背景散斑在柱状冲击波处的偏移量;基于第二公式及第三公式,通过散斑的偏移量确定光线折射率;基于第四公式及第五公式,通过光线折射率确定柱状冲击波的密度和压强;本发明不仅可以得到小范围内冲击波传播的压强分布,而且利用较少的光学设备,可以达到观察冲击波传播的目的。
Description
技术领域
本发明涉及压强分布测量领域,特别是一种针对水中柱状冲击波的压强分布测量方法。
背景技术
水中冲击波是指当液体介质中压强、密度和温度在波阵面上发生突跃变化的压缩波。
由于水的压缩性比较低,冲击波在水中的脉宽值比较窄,常温下,冲击波在水中的传播速度大于等于1500mm/s,所以测量水中冲击波需要高频响、高分辨率的压强传感器,对于水中冲击波单点的测量,目前使用比较多的是PVDF(聚偏氟乙烯)针式传感器、膜式水听器和光导纤维水听器,在多点测量甚至是对于冲击波压强分布测量的情况,传统的测量方法需要用多个压强传感器,如果是在比较小的空间里,布置多个压强传感器是比较困难的。
发明内容
本发明要解决的问题是能够测量在小范围内冲击波传播的压强分布和观察冲击波传播的图像。
为了解决上述问题,本发明提供了一种针对水中柱状冲击波的压强分布测量方法,该方法包括如下步骤:
采用摄像装置获取光源发出的光线经过散斑背景、无柱状冲击波的二维水槽的第一图像集;
采用摄像装置获取光源发出的光线经过散斑背景、加载柱状冲击波的二维水槽的第二图像集;
基于第一公式,通过所述第一图像集和所述第二图像集的对比确定散斑在所述柱状冲击波处的偏移量;
基于第二公式及第三公式,通过所述散斑的偏移量确定所述光线折射率;
基于第四公式及第五公式,通过所述光线折射率确定所述柱状冲击波的密度和压强。
进一步地,所述柱状冲击波为二维柱状冲击波。
进一步地,所述第一公式如下:
进一步地,所述第一公式中的Δx和Δy分别是在水平和垂直方向上的散斑的偏移量,Δt为高速摄像机记录的两张连可视化图片的时间间隔,计算散斑位移时,Δt取值为1,fijk是第k图像帧中第i行和第j列的交叉处的图像亮度,(fijk)x,(fijk)y,(fijk)t分别是x,y,t图像亮度的偏导数。
进一步地,所述第二公式如下:
所述第三公式如下:
进一步地,所述第二公式和第三公式中的ZD为散斑背景到柱状冲击波圆心的距离,ZB为散斑背景到相机镜头的距离,f为镜头的焦点距离。
进一步地,所述第四公式和第五公式如下:
所述第五公式如下:
进一步地,所述第四公式和第五公式中P为柱状冲击波的压强,ρ为柱状冲击波的密度,B=2963bar,γ=7.415,P0=1.03bar。
本发明所达到的有益效果是:本发明为一种针对水中柱状冲击波的压强分布测量方法,主要是测量二维柱状冲击波的压强分布,首先通过摄像装置获取散斑背景在柱状冲击波处发生偏移过程的图像,基于第一公式确定散斑背景的偏移量;根据散斑背景的偏移量,确定光线的折射率,根据光线折射率进一步确定柱状冲击波的密度和压强;其中,图像处理方法采用的是光流法中的时空间法,该方法与PIV中粒子图像处理方法相比,可以更精确的处理水中冲击波导致的小于一个像素的散斑偏移,另一方面,本发明中的压强分布测量方法,利用较少的光学设备,还可以达到观察冲击波传播的目的。
附图说明
图1是本发明实施例提供的一种针对水中柱状冲击波的压强分布测量方法流程图;
具体实施方式
为使本发明实施例的目的、技术方案和有点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方法进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动成果前提下所获得的所有其他实施例,都属于本发明保护的范围。
本实施例主要使用的实验设备包括:光源、散斑背景和高速摄像机。
图1为本发明实施例提供的一种针对水中柱状冲击波的压强分布测量方法的流程图,该方法具体包括:
S101、采用摄像装置获取光源发出的光线经过散斑背景、无柱状冲击波的二维水槽的第一图像集;
高速摄像机获取光源发出的光线经过散斑背景、无柱状冲击波的二维水槽时的图像。
S102、采用摄像装置获取光源发出的光线经过散斑背景、加载柱状冲击波的二维水槽的第二图像集;
冲击波的产生方法有很多种,例如激光、水中***、水中放电等,在本实施例中,是采用水中放电产生的冲击波进行实验的。
首先,在二维水槽产生冲击波后,由于二维水槽中柱状冲击波处的密度变化,光线经过柱状冲击波时发生折射,从而使散斑发生偏移。
此时,高速摄像机获取光源发出的光线经过散斑背景、产生柱状冲击波的二维水槽的图像。
S103、基于第一公式,通过第一图像集和第二图像集的对比,确定散斑在柱状冲击波处的偏移量;
根据第一公式,对S101与S102中的图像进行处理和计算,可以得出散斑在柱状冲击波处的偏移量。
其中第一公式如下:
第一公式中的Δx和Δy分别是在水平和垂直方向上的散斑位移,Δt为高速摄像机获取的两张连续可视化图像的时间间隔,计算散斑的偏移量时,Δt取值为1,fijk是第k图像帧中第i行和第j列的交叉处的图像亮度,(fijk)x,(fijk)y,(fijk)t分别是x,y,t图像亮度的偏导数。
S105、基于第二公式及第三公式,通过散斑的偏移量确定光线折射率;
其中第二公式如下:
第三公式如下:
第二公式和第三公式中的ZD为散斑背景到柱状冲击波圆心的距离,ZB为散斑背景到高速摄像机的镜头的距离,f为高速摄像机镜头的焦点距离。
S106、基于第四公式及第五公式,通过光线折射率确定柱状冲击波的密度和压强;
其中第四公式如下:
第五公式如下:
第四公式和第五公式中P为柱状冲击波的压强,ρ为柱状冲击波的密度,B=2963bar,γ=7.415,P0=1.03bar。
本发明中的柱状冲击波为二维柱状冲击波。
本发明中的散斑背景不仅仅是包括有规则和无规则斑点的背景,同时也包括其他样式的背景,如格子形状,条纹形状等。
Claims (3)
1.一种针对水中柱状冲击波的压强分布测量方法,其特征在于包括如下步骤:
采用摄像装置获取光源发出的光线经过散斑背景、无柱状冲击波的二维水槽的第一图像集;
采用摄像装置获取光源发出的光线经过散斑背景、加载柱状冲击波的二维水槽的第二图像集;
基于第一公式,通过所述第一图像集和所述第二图像集的对比确定散斑在所述柱状冲击波处的偏移量;
基于第二公式及第三公式,通过所述散斑的偏移量确定所述光线的折射率;
基于第四公式及第五公式,通过所述光线折射率确定所述柱状冲击波的密度和压强;
所述第一公式如下:
所述第一公式中的Δx和Δy分别是在水平和垂直方向上的所述散斑的偏移量,Δt为两张连可视化图像的时间间隔,fijk是第k图像帧中第i行和第j列的交叉处的图像亮度,(fijk)x,(fijk)y,(fijk)t分别是x,y,t相对应的图像亮度偏导数;
所述第二公式
如下:
所述第三公式如下:
所述第二公式和第三公式中的ZD为所述散斑背景到所述柱状冲击波圆心的距离,ZB为所述散斑背景到所述摄像装置的镜头的距离,f为所述摄像装置的镜头焦点距离;
所述第四公式如下:
所述第五公式如下:
所述第四公式和第五公式中的P为柱状冲击波的压强,为柱状冲击波密度,B=2963bar,γ=7.415,P0=1.03bar。
2.根据权利要求1所述的方法,其特征在于,所述柱状冲击波为二维冲击波。
3.根据权利要求1所述的方法,其特征在于,所述散斑背景至少包括下列之一:
有规则斑点、无规则斑点、格子形状和条纹形状。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811212878.6A CN109238553B (zh) | 2018-10-18 | 2018-10-18 | 一种针对水中柱状冲击波的压强分布测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811212878.6A CN109238553B (zh) | 2018-10-18 | 2018-10-18 | 一种针对水中柱状冲击波的压强分布测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109238553A CN109238553A (zh) | 2019-01-18 |
CN109238553B true CN109238553B (zh) | 2020-01-21 |
Family
ID=65053787
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811212878.6A Active CN109238553B (zh) | 2018-10-18 | 2018-10-18 | 一种针对水中柱状冲击波的压强分布测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109238553B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN86203799U (zh) * | 1986-06-11 | 1987-10-03 | 北京光电技术研究所 | 便携式激光动态全息—散斑照相装置 |
JP2010236911A (ja) * | 2009-03-30 | 2010-10-21 | Kyoto Institute Of Technology | レーザー誘起衝撃波を利用した分子の分析方法 |
CN103123286A (zh) * | 2012-11-28 | 2013-05-29 | 江苏大学 | 一种激光冲击波时空分布特性的检测方法和装置 |
CN103424221A (zh) * | 2013-03-27 | 2013-12-04 | 常州大学 | 一种可燃气体***冲击波能量的测定方法和装置 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6756579B2 (en) * | 2002-03-27 | 2004-06-29 | The United States Of America As Represented By The Secretary Of The Navy | Optical fiber based apparatus and sensor with enhanced signal for pressure measurements |
CN102494830A (zh) * | 2011-12-31 | 2012-06-13 | 北京理工大学 | 一种***冲击波场强度的测量方法及专用压力罐 |
CN103292943A (zh) * | 2012-02-27 | 2013-09-11 | 中国人民解放军总参谋部工程兵科研三所 | 非接触式冲击波超压测量方法 |
CN106706197B (zh) * | 2016-11-16 | 2019-03-05 | 哈尔滨工程大学 | 基于改进的霍普金森杆的水下***压力测量装置 |
CN108519401B (zh) * | 2018-03-27 | 2021-03-02 | 北京航空航天大学 | 一种1200℃有氧高温环境下抑制散斑图像抖动试验装置 |
CN108444353B (zh) * | 2018-04-19 | 2023-11-28 | 宁波交通工程建设集团有限公司 | 隧道开挖***冲击波防护装置及操作方法 |
-
2018
- 2018-10-18 CN CN201811212878.6A patent/CN109238553B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN86203799U (zh) * | 1986-06-11 | 1987-10-03 | 北京光电技术研究所 | 便携式激光动态全息—散斑照相装置 |
JP2010236911A (ja) * | 2009-03-30 | 2010-10-21 | Kyoto Institute Of Technology | レーザー誘起衝撃波を利用した分子の分析方法 |
CN103123286A (zh) * | 2012-11-28 | 2013-05-29 | 江苏大学 | 一种激光冲击波时空分布特性的检测方法和装置 |
CN103424221A (zh) * | 2013-03-27 | 2013-12-04 | 常州大学 | 一种可燃气体***冲击波能量的测定方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109238553A (zh) | 2019-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104966308B (zh) | 一种计算激光光束光斑大小的方法 | |
CN104303026B (zh) | 双光束非接触式位移传感器 | |
US10281264B2 (en) | Three-dimensional measurement apparatus and control method for the same | |
Dalziel et al. | Simultaneous synthetic schlieren and PIV measurements for internal solitary waves | |
US20100165152A1 (en) | Processing Images Having Different Focus | |
US20080080852A1 (en) | Single lens auto focus system for stereo image generation and method thereof | |
US11200714B2 (en) | Virtual image distance measurement method, apparatus and device | |
DE112010005646T5 (de) | Kameraabstandsmessvorrichtung | |
CN105445492B (zh) | 一种透过散射介质的激光散斑流速检测方法和装置 | |
EP3548838B1 (en) | Wire rope measuring device and wire rope measuring method | |
CN207894591U (zh) | 一种hud前挡玻璃检测设备 | |
Martins et al. | Ray-tracing based image correction of optical distortion for PIV measurements in packed beds | |
CN102506814A (zh) | 被动式可见光单目测距装置 | |
US20130057675A1 (en) | Method and arrangement for measuring flow rate of optically non-homogeneous material | |
Tsubaki et al. | Stereoscopic measurement of a fluctuating free surface with discontinuities | |
JP2015021862A (ja) | 3次元計測装置及び3次元計測方法 | |
Vinnichenko et al. | Accuracy of background oriented schlieren for different background patterns and means of refraction index reconstruction | |
JP2008249431A (ja) | 3次元画像補正方法及びその装置 | |
Wang et al. | Quantitative pressure measurement of shock waves in water using a schlieren-based visualization technique | |
CN109238553B (zh) | 一种针对水中柱状冲击波的压强分布测量方法 | |
CN103453850B (zh) | 基于数字散斑相关技术的透明液面微形貌测量方法和*** | |
Butterworth et al. | Realtime acousto-optical QA methods for high intensity fields | |
CN105333818B (zh) | 基于单目摄像机的3d空间测量方法 | |
Kumar et al. | Directional characteristics of spatially evolving young waves under steady wind | |
Sveen et al. | A dynamic masking technique for combined measurements of PIV and synthetic schlieren applied to internal gravity waves |
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 |