CN109033725A - 一种定床河工模型试验大面积床面切应力的估算方法 - Google Patents
一种定床河工模型试验大面积床面切应力的估算方法 Download PDFInfo
- Publication number
- CN109033725A CN109033725A CN201811079731.4A CN201811079731A CN109033725A CN 109033725 A CN109033725 A CN 109033725A CN 201811079731 A CN201811079731 A CN 201811079731A CN 109033725 A CN109033725 A CN 109033725A
- Authority
- CN
- China
- Prior art keywords
- bead
- bed
- moment
- shear stress
- fixed bed
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Abstract
本发明涉及一种定床河工模型试验大面积床面切应力的估算方法。本发明提供了一种利用底部颗粒运动轨迹估算底部床面切应力的方法,可以通过选择不同粒径的小球,减少水深对切应力测量的影响,通过本发明可以同步快速获取河工定床模型试验中大面积不同部位的床面切应力,大大提高了床面切应力估算效率。
Description
技术领域
本发明涉及水工测量技术领域,具体涉及一种定床河工模型试验大面积床面切应力的估算方法。
背景技术
水流在流动时床面底部对水流流动造成阻力,这种阻力通常可用床面切应力来表征。同时,床面切应力又是表征泥沙运动强度的重要参数之一,几乎所有的泥沙运动强度都可以用床面切应力来表示。目前床面切应力可通过局部传感器测量,或者通过底部流速分布反推,其中局部传感器法需要在床面设置传感器,较为复杂;流速分布则需要对局部流速分布进行测量,费时费力;以上两种方法通常只能测量河工模型中某一个点的床面切应力,效率较低。
发明内容
本发明所要解决的技术问题是提供一种定床河工模型试验大面积床面切应力的估算方法,解决床面切应力计算效率低的问题。
本发明解决上述技术问题的技术方案如下:一种定床河工模型试验大面积床面切应力的估算方法,包括以下步骤:
S1、在定床河工模型中放置多个小球;
S2、通过摄像机连续拍摄多幅该定床河工模型上小球的位置图像,并对每幅图上小球位置的坐标进行拟合,得到每个小球的运动轨迹;
S3、根据每个小球的运动轨迹计算每个小球的加速度;
S4、根据每个小球的加速度计算不同部位的床面切应力。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步,所述定床河工模型包括床面和设置在床面两侧的河岸,水流由床面的上游横截面流入,从床面的下游横截面流出,所述床面与实际河道形状类似,所述床面上方设有摄像机。
进一步,所述步骤S2中的坐标拟合方法为高元多项式法,具体步骤为:
S21、针对某一个小球,获取其在连续各帧图像中的位置坐标,每幅图像获得一个坐标数据(x,y),得到m幅图像的横坐标序列x1,x2,…,xm和纵坐标序列y1,y2…,ym;
S22、通过最小二乘法对坐标数据进行拟合,拟合公式为:
y=anxn+an-1xn-1+…a1x+a0 (1)
上式中,a0~an为待定参数,n为项数;
S23、根据不同的小球运动轨迹选择合适的n,拟合得到待定参数a0~an;
S24、根据横坐标序列x1~xm计算拟合后的横坐标序列x1a~x(m-2)a,计算公式为:
x(m-2)a=(x(m-2)+x(m-1)+xm)/3 (2);
S25、将拟合后的横坐标序列x1a~x(m-2)a代入公式(1)中,得到拟合后的纵坐标序列y1a~y(m-2)a,通过拟合后的横坐标序列x1a~x(m-2)a和拟合后的纵坐标序列y1a~y(m-2)a可得到该小球的运动轨迹。
进一步,所述步骤S3中加速度的计算公式为:
a=(v2-v1)/dt (3)
上式中,v1为小球在该时刻与前一时刻的运动速度,v2为小球在该时刻与后一时刻的运动速度,dt为两幅图像的间隔时间,其中v1和v2的计算公式分别为:
v1=ds1/dt1 (4)
v2=ds2/dt2 (5)
上式中,ds1为小球在该时刻与前一时刻的运动距离,dt1为该时刻与前一时刻的时间间隔,ds2为小球在该时刻与后一时刻的运动距离,dt2该时刻与后一时刻的时间间隔。
进一步,所述床面切应力的计算公式为:
上式中,p0为小球密度,d0为小球直径。
本发明的有益效果是:本发明提供了一种利用底部颗粒运动轨迹估算底部床面切应力的方法,可以通过选择不同粒径的小球,减少水深对切应力测量的影响,通过本发明可以同步快速获取河工定床模型试验中大面积不同部位的床面切应力,大大提高了床面切应力估算效率。
附图说明
图1为本发明步骤流程图;
图2为本发明中的定床河工模型示意图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
如图1所示,一种定床河工模型试验大面积床面切应力的估算方法,包括以下步骤:
S1、在定床河工模型中放置多个小球;
S2、通过摄像机拍摄多幅该定床河工模型上小球的位置图像,并对每幅图上小球位置的坐标进行拟合,得到每个小球的运动轨迹;
S3、根据每个小球的运动轨迹计算每个小球的加速度;
S4、根据每个小球的加速度计算不同部位的床面切应力。
如图2所示,定床河工模型包括床面和设置在床面两侧的河岸,水流由床面的上游横截面流入,从床面的下游横截面流出,床面与实际河道形状类似,床面上方设有摄像机。
在本发明实施例中,步骤S2中的坐标拟合方法为高元多项式法,具体步骤为:
S21、针对某一个小球,获取其在连续各帧图像中的位置坐标,每幅图像获得一个坐标数据(x,y),得到m幅图像的横坐标序列x1,x2,…,xm和纵坐标序列y1,y2…,ym;
S22、通过最小二乘法对坐标数据进行拟合,拟合公式为:
y=anxn+an-1xn-1+…a1x+a0 (1)
上式中,a0~an为待定参数,n为项数;
S23、根据不同的小球运动轨迹选择合适的n,拟合得到待定参数a0~an;
S24、根据横坐标序列x1~xm计算拟合后的横坐标序列x1a~x(m-2)a,计算公式为:
x(m-2)a=(x(m-2)+x(m-1)+xm)/3 (2);
S25、将拟合后的横坐标序列x1a~x(m-2)a代入公式(1)中,得到拟合后的纵坐标序列y1a~y(m-2)a,通过拟合后的横坐标序列x1a~x(m-2)a和拟合后的纵坐标序列y1a~y(m-2)a可得到该小球的运动轨迹。
在本发明实施例中,步骤S3中加速度的计算公式为:
a=(v2-v1)/dt (3)
上式中,v1为小球在该时刻与前一时刻的运动速度,v2为小球在该时刻与后一时刻的运动速度,dt为两幅图像的间隔时间,其中v1和v2的计算公式分别为:
v1=ds1/dt1 (4)
v2=ds2/dt2 (5)
上式中,ds1为小球在该时刻与前一时刻的运动距离,dt1为该时刻与前一时刻的时间间隔,ds2为小球在该时刻与后一时刻的运动距离,dt2该时刻与后一时刻的时间间隔。
在本发明实施例中,床面切应力的计算公式为:
上式中,p0为小球密度,d0为小球直径。
在本发明实施例中,在上游横截面上放入多个已知粒径的白色玻璃小球,白色玻璃小球的大小通过定床河工模型的水深选择(玻璃小球的直径d小于水深h的1/5),摄像机所拍摄的图像应确保每个玻璃小球在单个方向上有2-3个像素点,摄像机的拍摄频率为可根据小球大小、水流运动速度选择,小球运动速度越快则所需摄像机拍摄帧频越高,通常可选用25fps,即每秒钟至少拍摄25幅画面。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种定床河工模型试验大面积床面切应力的估算方法,其特征在于,包括以下步骤:
S1、在定床河工模型中放置多个小球;
S2、通过摄像机连续拍摄多幅该定床河工模型上小球的位置图像,并对每幅图上小球位置的坐标进行拟合,得到每个小球的运动轨迹;
S3、根据每个小球的运动轨迹计算每个小球的加速度;
S4、根据每个小球的加速度计算不同部位的床面切应力。
2.根据权利要求1所述的定床河工模型试验大面积床面切应力的估算方法,其特征在于,所述定床河工模型包括床面和设置在床面两侧的河岸,水流由床面的上游横截面流入,从床面的下游横截面流出,所述床面与实际河道形状类似,所述床面上方设有摄像机。
3.根据权利要求1所述的定床河工模型试验大面积床面切应力的估算方法,其特征在于,所述步骤S2中的坐标拟合方法为高元多项式法,具体步骤为:
S21、针对某一个小球,获取其在连续各帧图像中的位置坐标,每幅图像获得一个坐标数据(x,y),得到m幅图像的横坐标序列x1,x2,…,xm和纵坐标序列y1,y2…,ym;
S22、通过最小二乘法对坐标数据进行拟合,拟合公式为:
y=anxn+an-1xn-1+…a1x+a0 (1)
上式中,a0~an为待定参数,n为项数;
S23、根据不同的小球运动轨迹选择合适的n,拟合得到待定参数a0~an;
S24、根据横坐标序列x1~xm计算拟合后的横坐标序列x1a~x(m-2)a,计算公式为:
x(m-2)a=(x(m-2)+x(m-1)+xm)/3 (2);
S25、将拟合后的横坐标序列x1a~x(m-2)a代入公式(1)中,得到拟合后的纵坐标序列y1a~y(m-2)a,通过拟合后的横坐标序列x1a~x(m-2)a和拟合后的纵坐标序列y1a~y(m-2)a可得到该小球的运动轨迹。
4.根据权利要求1所述的定床河工模型试验大面积床面切应力的估算方法,其特征在于,所述步骤S3中加速度的计算公式为:
a=(v2-v1)/dt (3)
上式中,v1为小球在该时刻与前一时刻的运动速度,v2为小球在该时刻与后一时刻的运动速度,dt为两幅图像的间隔时间,其中v1和v2的计算公式分别为:
v1=ds1/dt1 (4)
v2=ds2/dt2 (5)
上式中,ds1为小球在该时刻与前一时刻的运动距离,dt1为该时刻与前一时刻的时间间隔,ds2为小球在该时刻与后一时刻的运动距离,dt2该时刻与后一时刻的时间间隔。
5.根据权利要求1所述的定床河工模型试验大面积床面切应力的估算方法,其特征在于,所述床面切应力的计算公式为:
上式中,p0为小球密度,d0为小球直径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811079731.4A CN109033725B (zh) | 2018-09-14 | 2018-09-14 | 一种定床河工模型试验大面积床面切应力的估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811079731.4A CN109033725B (zh) | 2018-09-14 | 2018-09-14 | 一种定床河工模型试验大面积床面切应力的估算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109033725A true CN109033725A (zh) | 2018-12-18 |
CN109033725B CN109033725B (zh) | 2020-05-05 |
Family
ID=64622092
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811079731.4A Active CN109033725B (zh) | 2018-09-14 | 2018-09-14 | 一种定床河工模型试验大面积床面切应力的估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109033725B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111400974A (zh) * | 2020-04-27 | 2020-07-10 | 中国水利水电科学研究院 | 一种估算矩形河渠壁面和床面切应力的方法 |
CN112033641A (zh) * | 2020-08-26 | 2020-12-04 | 大连理工大学 | 一种海岸地貌演变物理模型试验方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999046448A1 (en) * | 1998-03-12 | 1999-09-16 | The Texas A & M University System | Apparatus and methods for prediction of scour related information in soils |
CN201037868Y (zh) * | 2007-03-13 | 2008-03-19 | 魏润杰 | 一种高速粒子图像测速装置 |
CN101776449A (zh) * | 2010-01-11 | 2010-07-14 | 中国科学院力学研究所 | 利用超声波测量砂土床面冲刷形貌的方法及装置 |
CN102359862A (zh) * | 2011-08-12 | 2012-02-22 | 河海大学 | 粉沙质和淤泥质海岸泥沙运动数值模拟方法 |
WO2013026934A2 (en) * | 2011-08-24 | 2013-02-28 | Institute Of Technology Sligo | An apparatus and method for determining the shear strength of sediments on the floor of a body of water |
WO2014096424A1 (de) * | 2012-12-21 | 2014-06-26 | Universität Innsbruck | Geschiebemessung mittels positions- oder formveränderlicher störkörper |
US20160070829A1 (en) * | 2014-09-05 | 2016-03-10 | Chevron U.S.A. Inc. | Sediment transport simulation with parameterized templates for depth profiling |
CN106503421A (zh) * | 2016-09-28 | 2017-03-15 | 西南交通大学 | 一种边坡渐进破坏全过程计算方法 |
CN107145678A (zh) * | 2017-05-22 | 2017-09-08 | 中国水利水电科学研究院 | 一种河道平面二维水流模型糙率的率定方法 |
CN107436977A (zh) * | 2017-07-24 | 2017-12-05 | 珠江水利委员会珠江水利科学研究院 | 复杂河网分流的数值模拟方法 |
-
2018
- 2018-09-14 CN CN201811079731.4A patent/CN109033725B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999046448A1 (en) * | 1998-03-12 | 1999-09-16 | The Texas A & M University System | Apparatus and methods for prediction of scour related information in soils |
CN201037868Y (zh) * | 2007-03-13 | 2008-03-19 | 魏润杰 | 一种高速粒子图像测速装置 |
CN101776449A (zh) * | 2010-01-11 | 2010-07-14 | 中国科学院力学研究所 | 利用超声波测量砂土床面冲刷形貌的方法及装置 |
CN102359862A (zh) * | 2011-08-12 | 2012-02-22 | 河海大学 | 粉沙质和淤泥质海岸泥沙运动数值模拟方法 |
WO2013026934A2 (en) * | 2011-08-24 | 2013-02-28 | Institute Of Technology Sligo | An apparatus and method for determining the shear strength of sediments on the floor of a body of water |
WO2014096424A1 (de) * | 2012-12-21 | 2014-06-26 | Universität Innsbruck | Geschiebemessung mittels positions- oder formveränderlicher störkörper |
US20160070829A1 (en) * | 2014-09-05 | 2016-03-10 | Chevron U.S.A. Inc. | Sediment transport simulation with parameterized templates for depth profiling |
CN106503421A (zh) * | 2016-09-28 | 2017-03-15 | 西南交通大学 | 一种边坡渐进破坏全过程计算方法 |
CN107145678A (zh) * | 2017-05-22 | 2017-09-08 | 中国水利水电科学研究院 | 一种河道平面二维水流模型糙率的率定方法 |
CN107436977A (zh) * | 2017-07-24 | 2017-12-05 | 珠江水利委员会珠江水利科学研究院 | 复杂河网分流的数值模拟方法 |
Non-Patent Citations (3)
Title |
---|
JAE HYEONPARK 等: "Directmeasurementofbottomshearstressunderhigh-velocity flow conditions", 《FLOW MEASUREMENTAND INSTRUMENTATION》 * |
L.C.HOPKINSON 等: "COMPARISON OF DIRECT AND INDIRECT BOUNDARY SHEAR STRESS MEASUREMENTS ALONG VEGETATED STREAMBANKS", 《RIVER RESEARCH AND APPLICATIONS》 * |
黄海龙 等: "二维波流共线作用下全直接法床面切应力测量仪设计研究", 《水道港口》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111400974A (zh) * | 2020-04-27 | 2020-07-10 | 中国水利水电科学研究院 | 一种估算矩形河渠壁面和床面切应力的方法 |
CN112033641A (zh) * | 2020-08-26 | 2020-12-04 | 大连理工大学 | 一种海岸地貌演变物理模型试验方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109033725B (zh) | 2020-05-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103308714B (zh) | 一种水流中彩色粒子图像测速系统及其测量方法 | |
CN105913454B (zh) | 一种视频图像中运动目标的像素坐标轨迹预测方法 | |
CN109033725A (zh) | 一种定床河工模型试验大面积床面切应力的估算方法 | |
US20050018882A1 (en) | Controlled surface wave image velocimetry | |
CN104410855A (zh) | 一种监控视频的抖动检测方法 | |
CN102944507B (zh) | 一种轻质异形颗粒曳力系数的测量装置及测量方法 | |
CN102331511B (zh) | Piv图像高频采集方法 | |
CN102288778B (zh) | 基于网格对象跟踪的实时泥石流表面速度测量方法 | |
CN103096117B (zh) | 视频噪声检测方法及装置 | |
Li et al. | Large-scale particle tracking velocimetry with multi-channel CCD cameras | |
CN103208125A (zh) | 视频帧图像中颜色和运动全局对比度的视觉显著性算法 | |
CN102620903A (zh) | 一种水流中流迹线的动态显示系统及其显示方法 | |
CN104568376A (zh) | 通过图像分析卵砾石瞬时输沙强度的方法以及系统 | |
Li et al. | Blinkflow: A dataset to push the limits of event-based optical flow estimation | |
Bung et al. | Image processing techniques for velocity estimation in highly aerated flows: Bubble Image Velocimetry vs. Optical Flow | |
CN105741277A (zh) | 基于ViBe算法与SLIC超像素的背景差分法 | |
Admiraal et al. | Case study: Particle velocimetry in a model of lake Ogallala | |
Yang et al. | The 3-D spread of saltation sand over a flat bed surface in aeolian sand transport | |
CN114119670A (zh) | 一种基于智能手机获取河流视频的流速测量方法 | |
CN103020982B (zh) | 一种基于跟踪反馈的非参数运动目标检测方法 | |
CN104992149A (zh) | 拌线检测方法、系统、网络摄像机及监控系统 | |
CN116129362B (zh) | 一种基于坐标横向断面图的河流漂浮污染物的监测方法 | |
Cauteruccio et al. | Wind Tunnel Validation of a Particle Tracking Model to Evaluate the Wind-Induced Bias of Precipitation Measurements. | |
CN116129363B (zh) | 一种河流漂浮污染物的监测方法及存储介质 | |
Miranda-Tello et al. | Characterization of the level fluctuations in a physical model of the steel continuous casting mold through image processing |
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 |