CN109033725B - 一种定床河工模型试验大面积床面切应力的估算方法 - Google Patents
一种定床河工模型试验大面积床面切应力的估算方法 Download PDFInfo
- Publication number
- CN109033725B CN109033725B CN201811079731.4A CN201811079731A CN109033725B CN 109033725 B CN109033725 B CN 109033725B CN 201811079731 A CN201811079731 A CN 201811079731A CN 109033725 B CN109033725 B CN 109033725B
- Authority
- CN
- China
- Prior art keywords
- bed surface
- shear stress
- bed
- small ball
- 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.)
- Active
Links
Images
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
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 (4)
1.一种定床河工模型试验大面积床面切应力的估算方法,其特征在于,包括以下步骤:
S1、在定床河工模型中放置多个小球;
S2、通过摄像机连续拍摄多幅该定床河工模型上小球的位置图像,并对每幅图上小球位置的坐标进行拟合,得到每个小球的运动轨迹;
S3、根据每个小球的运动轨迹计算每个小球的加速度;
S4、根据每个小球的加速度计算不同部位的床面切应力;
所述定床河工模型包括床面和设置在床面两侧的河岸,水流由床面的上游横截面流入,从床面的下游横截面流出,所述床面上方设有摄像机。
2.根据权利要求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可得到该小球的运动轨迹。
3.根据权利要求1所述的定床河工模型试验大面积床面切应力的估算方法,其特征在于,所述步骤S3中加速度的计算公式为:
a=(v2-v1)/dt (3)
上式中,v1为小球在当前时刻与前一时刻的运动速度,v2为小球在当前时刻与后一时刻的运动速度,dt为两幅图像的间隔时间,其中v1和v2的计算公式分别为:
v1=ds1/dt1 (4)
v2=ds2/dt2 (5)
上式中,ds1为小球在当前时刻与前一时刻的运动距离,dt1为当前时刻与前一时刻的时间间隔,ds2为小球在当前时刻与后一时刻的运动距离,dt2当前时刻与后一时刻的时间间隔。
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 CN109033725A (zh) | 2018-12-18 |
CN109033725B true 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) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111400974B (zh) * | 2020-04-27 | 2020-12-08 | 中国水利水电科学研究院 | 一种估算矩形河渠壁面和床面切应力的方法 |
CN112033641B (zh) * | 2020-08-26 | 2021-05-18 | 大连理工大学 | 一种海岸地貌演变物理模型试验方法 |
Citations (8)
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 | 魏润杰 | 一种高速粒子图像测速装置 |
CN102359862A (zh) * | 2011-08-12 | 2012-02-22 | 河海大学 | 粉沙质和淤泥质海岸泥沙运动数值模拟方法 |
WO2013026934A3 (en) * | 2011-08-24 | 2013-09-19 | 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 |
CN106503421A (zh) * | 2016-09-28 | 2017-03-15 | 西南交通大学 | 一种边坡渐进破坏全过程计算方法 |
CN107145678A (zh) * | 2017-05-22 | 2017-09-08 | 中国水利水电科学研究院 | 一种河道平面二维水流模型糙率的率定方法 |
CN107436977A (zh) * | 2017-07-24 | 2017-12-05 | 珠江水利委员会珠江水利科学研究院 | 复杂河网分流的数值模拟方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101776449A (zh) * | 2010-01-11 | 2010-07-14 | 中国科学院力学研究所 | 利用超声波测量砂土床面冲刷形貌的方法及装置 |
US10108760B2 (en) * | 2014-09-05 | 2018-10-23 | Chevron U.S.A. Inc. | Sediment transport simulation with parameterized templates for depth profiling |
-
2018
- 2018-09-14 CN CN201811079731.4A patent/CN109033725B/zh active Active
Patent Citations (8)
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 | 魏润杰 | 一种高速粒子图像测速装置 |
CN102359862A (zh) * | 2011-08-12 | 2012-02-22 | 河海大学 | 粉沙质和淤泥质海岸泥沙运动数值模拟方法 |
WO2013026934A3 (en) * | 2011-08-24 | 2013-09-19 | 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 |
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 |
---|
COMPARISON OF DIRECT AND INDIRECT BOUNDARY SHEAR STRESS MEASUREMENTS ALONG VEGETATED STREAMBANKS;L.C.HOPKINSON 等;《RIVER RESEARCH AND APPLICATIONS》;20160126;第1755-1764页 * |
Directmeasurementofbottomshearstressunderhigh-velocity flow conditions;Jae HyeonPark 等;《Flow Measurementand Instrumentation》;20151130;第121-127页 * |
二维波流共线作用下全直接法床面切应力测量仪设计研究;黄海龙 等;《水道港口》;20161031;第37卷(第5期);第563-568页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109033725A (zh) | 2018-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105913454B (zh) | 一种视频图像中运动目标的像素坐标轨迹预测方法 | |
CN109033725B (zh) | 一种定床河工模型试验大面积床面切应力的估算方法 | |
CN104575003B (zh) | 一种基于道路监控视频的车辆速度检测方法 | |
CN103698553B (zh) | 基于红外摄像的新型表面流场测速系统和测速方法 | |
Evangelista et al. | Simulation of dam-break waves on movable beds using a multi-stage centered scheme | |
CN106204638A (zh) | 一种基于尺度自适应和遮挡处理的航拍目标跟踪方法 | |
Ramesh et al. | Near-bed particle motion over transitionally-rough bed | |
CN110263075B (zh) | Ladcp与usbl组合观测及数据资料处理方法 | |
CN110132247B (zh) | 一种多空间尺度的水体流场观测互补方法 | |
CN102288778B (zh) | 基于网格对象跟踪的实时泥石流表面速度测量方法 | |
CN105632170A (zh) | 一种基于Mean Shift跟踪算法的交通流检测方法 | |
CN104766335A (zh) | 一种岩土材料变形数字图像相关分析优化方法 | |
CN105354863A (zh) | 基于特征滤波和快速运动检测模板预测的自适应尺度图像序列目标跟踪方法 | |
CN108625337B (zh) | 一种确定潮流界以下沙质河床河段整治水位的方法 | |
CN102879301A (zh) | 液态金属流动性瞬时流速、瞬时流长检测方法 | |
CN103778646B (zh) | 高速运动标记点图像快速搜索方法 | |
CN103745118B (zh) | 一种基于磁偶极子等效源法的地磁异常数据网格化方法 | |
Admiraal et al. | Case study: Particle velocimetry in a model of lake Ogallala | |
CN115588027B (zh) | 一种落石冲击力非接触式重构方法及系统 | |
EP1178436A3 (en) | Image measurement method, image measurement apparatus and image measurement program storage medium | |
Yang et al. | The 3-D spread of saltation sand over a flat bed surface in aeolian sand transport | |
CN114972409A (zh) | 一种运动粒子的二维及三维示踪轨迹续连方法 | |
Kumar et al. | Experimental study of flow field and movement of sediment over a ramp | |
Ge et al. | Analysis of overestimation in historical ice flow velocity maps in Western Pacific Ocean Sector, Antarctica | |
CN103243675A (zh) | 一种适用于水库库湾的水流垂向紊动测量装置和方法 |
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 |