考虑河岸变形的三维水沙数值模拟研究

将粘性河岸崩塌模拟力学方法与水沙模型相结合,构建了考虑河岸变形的三维数值模型。基于非正交网格,采用局部网格可动技术处理由河岸崩塌引起的河道摆动过程。模型采用有限体积法对方程进行离散,采用与动量插值技术相结合的SIMPLEC算法进行模型求解。通过弯道发展过程的计算表明,模型能够有效地模拟出河道的复杂演变过程。     

生产堤溃口展宽过程的概化模型试验研究

为研究生产堤溃口横向展宽过程的发展规律,以黄河下游生产堤土体为材料,开展概化模型试验,分析溃口处水深、流速、流量以及溃口宽度的变化。试验结果显示,溃口横向展宽的机理主要表现为水流冲刷堤身下层土体伴随上层悬臂土体发生绕轴崩塌;随着溃口展宽,堤前水深持续下降但下降速率逐渐减小,溃口流速呈先增大后减小的趋势,流速峰值约为1.1 m/s;随着主槽流量增加,溃口宽度发展至0.7 m时,溃口的水深、泄流流量、流速、展宽速率均有所增加。基于对溃口展宽的影响因素分析,提出溃口横向展宽速率的计算关系式,计算溃口展宽过程与实测值符合较好。研究结果不仅有助于生产堤展宽机理及发展规律的研究,还可为数学模型的建立与验证提供参考资料。     

生产堤溃口展宽过程的概化模型试验研究

为研究生产堤溃口横向展宽过程的发展规律,以黄河下游生产堤土体为材料,开展概化模型试验,分析溃口处水深、流速、流量以及溃口宽度的变化。试验结果显示,溃口横向展宽的机理主要表现为水流冲刷堤身下层土体伴随上层悬臂土体发生绕轴崩塌;随着溃口展宽,堤前水深持续下降但下降速率逐渐减小,溃口流速呈先增大后减小的趋势,流速峰值约为1.1 m/s;随着主槽流量增加,溃口宽度发展至0.7 m时,溃口的水深、泄流流量、流速、展宽速率均有所增加。基于对溃口展宽的影响因素分析,提出溃口横向展宽速率的计算关系式,计算溃口展宽过程与实测值符合较好。研究结果不仅有助于生产堤展宽机理及发展规律的研究,还可为数学模型的建立与验证提供参考资料。     

无粘性土堤漫顶溃决水流运动试验研究

堤防溃决时溃口附近水流形态及溃决流量变化过程的正确描述是指导溃口堵复及时避险、预测溃堤洪水演进并进一步评估淹没损失的前提条件.通过弯道水槽中的概化模型试验,根据无粘性土堤防漫顶溃决的溃口变化和水流运动特征,将无粘性土堤防溃决过程分为漫流、冲槽、展宽及稳定4个阶段.试验系统研究了堤防溃决后河道以及溃口附近水位的变化过程,研究了河道洪水流量、洪水位和筑堤材料对溃口处水位的影响.基于大量的试验资料,进一步归纳总结出溃决流量与溃口口门宽度及溃口处水头变化之间的关系.     

下荆江二元结构河岸土体特性及崩岸机理

三峡水库蓄水后下荆江河段河床冲刷下切, 局部河段崩岸险情频繁发生。为研究下荆江二元结构河岸的土体特性及崩岸机理,结合近期该河段崩岸情况,现场查勘了6个崩岸点,并对河岸土体进行了室内土工试验。试验结果表明下荆江河岸土体的垂向组成具有典型的二元结构特征:下部非粘性土（沙土）层较厚,上部粘性土（低液限粘土）层较薄且松散。以河岸崩塌过程分析为基础,提出了二元结构河岸发生绕轴崩塌时上部土层稳定性的计算方法。结合近岸水动力条件计算及土工试验结果,定量分析了二元结构河岸的崩塌机理及其影响因素:① 下部沙土层的起动流速比近岸流速小得多,故该土层容易受水流冲刷;②上部粘性土层崩塌前抗冲强度很大,但多为低液限粘土且相对松散,崩塌后堆积在岸边容易分解并被水流带走;③ 河岸稳定安全系数在一个水文年内呈周期性变化,落水期安全系数最小,故容易引发崩岸,该计算结果与近期崩岸实际统计结果一致。     

岸滩侧蚀崩塌速率的试验量测方法

岸滩侧蚀崩塌现象普遍存在于冲积河流演变过程中,是河道治理需重点关注的问题之一,岸滩侧蚀崩塌速率的准确量测则是崩岸机理及其治理措施研究的关键基础。基于图像处理技术,建立了河岸侧蚀崩塌速率试验量测系统,通过图像追踪河岸模型上边缘示踪网格位置变迁,提出了一种河岸侧蚀崩塌速率测量方法。以试验水槽为例,对河岸侧蚀崩塌速率进行了量测,并与经验公式计算值进行了比对,两者基本吻合,可较好地反映流速、土体级配等岸滩侧蚀崩塌速率的影响因素及其变化规律。研究成果可为江河崩岸机理的深入研究提供技术支撑。     

动态实时洪水风险分析及管理系统

当前洪水风险分析按照典型设计标准洪水进行计算的模式难以满足实际防洪管理需要,为了提高洪水风险分析的实时性以及适应洪水演进的动态性,设计了动态实时洪水风险分析框架。在本框架中,先采用一维和二维动态耦合水动力学数值方法耦合溃堤模型,然后在樵桑联围防洪保护区建立洪水演进模拟模型,通过灵活处理模型计算边界条件以及动态设置溃堤功能,计算不同设计标准洪水发生时,堤防出现单一溃口或者组合溃口后保护区内洪水演进过程。按照上述框架开发了樵桑联围动态实时洪水风险图编制与管理应用系统,并利用历史洪水资料开展模型验证,验证结果表明,2008-06洪水马口站、三水站、大熬站、甘竹（一）站的实测最高水位和模型计算最高水位的绝对误差分别为-0.10、0.10、0.09、0.04 m,均满足洪水模拟精度要求。利用模型计算了西江发生200年一遇的洪水情况下,江根堤防出现溃口后的洪水流量及溃口内外洪水水位变化过程,模拟溃口宽度168 m,最大溃口洪水流量达到5 190 m³,分析了堤防溃决后3、6和24 h洪水漫延导致村落淹没情况,结果表明其满足合理性分析。     

分层岸滩侧蚀坍塌过程及其水动力响应模拟

非均质岸滩广泛分布于冲积河流中,其侧蚀过程沿垂向具有分层特点。基于全三维水沙模型及河岸侧蚀坍塌力学模式,提出动态网格跟踪技术,构建了非均质岸滩侧蚀沿垂向差异的三维动力学模拟方法。以连续弯道概化水槽为例,模拟分析了非均质河岸侧蚀过程中三维水流结构的响应特征。研究结果表明:坡脚冲刷后,底部主流向凹岸偏移,在已有弯道环流的基础上,于坡脚处出现一反向次生流;上部粘性土层坍塌后,崩塌土体堆积于坡脚处,上部主流明显左偏,下部次生流消失;随着坡脚堆积体的冲刷搬运,下部主流进一步向凹岸偏移;如此循环,主流不断向凹岸偏移,致使岸坡持续崩退、河道摆动。     

分层岸滩侧蚀坍塌过程及其水动力响应模拟

非均质岸滩广泛分布于冲积河流中,其侧蚀过程沿垂向具有分层特点。基于全三维水沙模型及河岸侧蚀坍塌力学模式,提出动态网格跟踪技术,构建了非均质岸滩侧蚀沿垂向差异的三维动力学模拟方法。以连续弯道概化水槽为例,模拟分析了非均质河岸侧蚀过程中三维水流结构的响应特征。研究结果表明:坡脚冲刷后,底部主流向凹岸偏移,在已有弯道环流的基础上,于坡脚处出现一反向次生流;上部粘性土层坍塌后,崩塌土体堆积于坡脚处,上部主流明显左偏,下部次生流消失;随着坡脚堆积体的冲刷搬运,下部主流进一步向凹岸偏移;如此循环,主流不断向凹岸偏移,致使岸坡持续崩退、河道摆动。     

三峡工程运用初期石首河弯河势演变三维数值模拟

采用考虑河岸坍塌的三维水沙数值模型,对三峡工程运用初期石首河弯的河势演变进行了模拟分析.数值模型中采用基于非正交网格的局部网格可动技术,将传统水沙模型与二元结构河岸坍塌力学模型相结合.利用三峡水库蓄水前后2002~2004年石首河段冲淤演变资料对模型进行了验证,模拟结果与实测结果吻合较好;模型不仅能模拟出河道的垂向冲淤过程,而且还可模拟出由河岸坍塌所引起的河道横向摆动过程.计算分析了三峡工程运行至2016年末河段的河势演变情况,结果可为河段的河道整治提供参考.     

小尺度断层定量预测及对注水开发的影响——以渤中34油田为例

大尺度断层往往控制了沉积盆地的形成和油气成藏,而小尺度断层则影响着注水开发效果和剩余油分布.大尺度断层可以通过二维或三维地震资料识别,而小尺度断层的识别则特别困难.本文提出了一种基于断层分形生长模式和三维地质力学模拟相结合来确定小尺度断层的数量、发育位置和方位的方法,并根据油田开发动态资料来确定小尺度断层对注水开发和剩余油分布的影响.将地震上识别的大尺度断层引入到三维数值力学模型中,模拟大尺度断层形成时期断裂带附近的应力扰动作用,然后结合破裂准则来建立最易发生破裂的方位和最大库伦剪切应力网格,以这两套网格和断层尺度的幂律分布确定的小尺度断层数量为约束条件来确定随机模型,对小尺度断层的密度、产状和发育位置进行定量预测.研究表明:利用分形理论和三维地质力学模拟可以对大尺度断层伴生小尺度断层进行有效预测;小尺度断层对注水开发效果和剩余油分布的影响取决于小尺度断层的规模(断距)以及小尺度断层方位和注采方向的关系.     

渐进坍塌型崩岸的力学机制及模拟

针对渐进坍塌型崩岸,从土力学和河流动力学两方面理论出发,建立了岸坡稳定的力学模式,结合室内概化模拟试验和数值计算,分析了岸坡稳定或破坏的力学机制,揭示了缓坡出现崩岸的原因。结果表明,岸坡坡脚未受水流冲失时,若坡内渗流出逸坡降小于渗透破坏的临界坡降,岸坡处于稳定状态,当坡脚被水流冲失后,渗流渗径缩短,水土结合处坡面出逸坡降增大,大于临界坡降时则出现渗透破坏,引起局部小幅度土体崩塌,其后部土体失去支撑而陆续产生失稳破坏,随着时间的延长,土体崩塌现象逐步向后发展,最终导致岸坡整体崩塌破坏。     

Three-dimensional discrete element method parallel computation of Cauchy stress distribution over granular materials

The paper presents Cauchy stress tensor computation over parallel grids of message passing interface (MPI) parallel three-dimensional (3D) discrete element method (DEM) simulations of granular materials, considering spherical and nonspherical particles. The stress tensor computation is studied for quasi-static and dynamic conditions, and its resulting symmetry or asymmetry is discussed within the context of classical continuum mechanics (CCM), granular materials mechanics (GMM), and micropolar continuum mechanics (MCM). The average Cauchy stress tensor computation follows Bagi's and Nicot's formulations and is verified within MPI parallel 3D DEM simulations involving dynamically adaptive compute grids. These grids allow calculation of temporal and spatial distributions of stress across granular materials under static and dynamic conditions. The vertical stress component in gravitationally deposited particle assemblies exhibits nonuniform spatial distributions under static equilibrium, and its zone of maximum value changes during the process of gravitational pluviation and collapse. These phenomena reveal a microstructural effect on stress distribution within granular materials that is attributed to their discrete particulate nature (particle size, shape, gradation, boundary conditions, etc).     

具有行蓄洪区的河道流量演算方法探讨

以马斯京根流量演算法为基础提出了具有行蓄洪区时防洪系统流量演算的一般方法,是马斯京根法在具有行蓄洪区的河道洪水演算的推广应用.对行洪区的入流处理进行了讨论并提出了计算分流比的实用方法.最后以淮河中游王家坝到鲁台子河段为例,根据工程布局、水文站和水位站的分布将该河段分成10段,对12场洪水进行模拟预报计算,精度较高     

地震动加速度幅值对土坡失稳后滑坡体大小及致灾范围的影响

地震滑坡的致灾范围是判断滑坡能否会对已有建构筑物造成损失、确定预警疏散范围的重要依据,因此对地震土坡破坏后的滑坡体大小和致灾范围进行研究具有重要的意义.本研究基于SPH动力分析方法,结合弹塑性本构模型和固体力学控制方程建立了地震土坡破坏的动力分析模型;通过设置振动边界粒子和自由场边界粒子,实现了地震动加速度的施加以及自...     

Model for prediction of sea dike breaching initiated by breaking wave impact

A computational model system is proposed for the prediction of sea dike breaching initiated from the seaward side by breaking wave impact with the focus on the application of the model system for the estimation failure probability of the defence structure. The described model system is built using a number of existing models for the calculation of grass, clay, and sand erosion. The parameters identified as those having the most significant influence on the estimation of the failure have been described stochastically. Monte Carlo simulations to account for uncertainties of the relevant input parameters and the model itself have been performed and the probabilities of the breach initiation and of the full dike breaching have been calculated. This will form the basis to assess the coastal flood risk due to dike breaching.     

探地雷达在水利设施现状及隐患探测中的应用

简述探地雷达的一些基本原理,介绍探地雷达用于水利工程设施现状及隐患探测的3个实例,包括拦洪闸底板现状探测;海滨防浪堤隐患探测以及江堤滑塌成因调查。     

Three-dimensional hydromechanical modeling of internal erosion in dike-on-foundation

Currently, numerical studies at the real scale of an entire engineering structure considering internal erosion are still rare. This paper presents a three-dimensional (3D) numerical simulation of the effects of internal erosion within a linear dike located on a foundation. A two-dimensional (2D) finite element code has been extended to 3D in order to analyze the impact of internal erosion under more realistic hydromechanical conditions. The saturated soil has been considered as a mixture of four interacting constituents: soil skeleton, erodible fines, fluidized fine particles, and fluid. The detachment and transport of the fine particles have been modeled with a mass exchange model between the solid and the fluid phases. An elastoplastic constitutive model for sand-silt mixtures has been developed to monitor the effect of the evolution of both the porosity and the fines content induced by internal erosion upon the behavior of the soil skeleton. An unsaturated flow condition has been implemented into this coupled hydromechanical model to describe more accurately the seepage within the dike and the foundation. A stabilized finite element method was used to eliminate spurious numerical oscillations in solving the convection-dominated transport of fluidized particles. This numerical tool was then applied to a specific dike-on-foundation case subjected to internal erosion induced by a leakage located at the bottom of the foundation. Different failure modes were observed and analyzed for different boundary conditions, including the significant influence of the leakage cavity size and the elevation of the water level at the upstream and downstream sides of the dike.     

Continuous velocity fields for collapse and blowout of a pressurized tunnel face in purely cohesive soil

Face stability analysis of tunnels excavated under pressurized shields is a major issue in real tunnelling projects. Most of the failure mechanisms used for the stability analysis of tunnels in purely cohesive soils were derived from rigid block failure mechanisms that were developed for frictional soils, by imposing a null friction angle. For a purely cohesive soil, this kind of mechanism is quite far from the actual velocity field. This paper aims at proposing two new continuous velocity fields for both collapse and blowout of an air‐pressurized tunnel face. These velocity fields are much more consistent with the actual failures observed in undrained clays. They are based on the normality condition, which states that any plastic deformation in a purely cohesive soil develops without any volume change. The numerical results have shown that the proposed velocity fields significantly improve the best existing bounds for collapse pressures and that their results compare reasonably well with the collapse and blowout pressures provided by a commercial finite difference software, for a much smaller computational cost. A design chart is provided for practical use in geotechnical engineering. Copyright © 2012 John Wiley & Sons, Ltd.