会员注册 | 登录 | 微信快捷登录 支付宝快捷登录 QQ登录 微博登录 | 帮助中心 人人文库renrendoc.com美如初恋!
站内搜索 百度文库

热门搜索: 直缝焊接机 矿井提升机 循环球式转向器图纸 机器人手爪发展史 管道机器人dwg 动平衡试验台设计

   首页 人人文库网 > 资源分类 > DOC文档下载

水利工程论文-扰动淤泥与沉积淤泥的流变特性研究.doc

  • 资源星级:
  • 资源大小:26.51KB   全文页数:8页
  • 资源格式: DOC        下载权限:注册会员/VIP会员
您还没有登陆,请先登录。登陆后即可下载此文档。
  合作网站登录: 微信快捷登录 支付宝快捷登录   QQ登录   微博登录
友情提示
2:本站资源不支持迅雷下载,请使用浏览器直接下载(不支持QQ浏览器)
3:本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

水利工程论文-扰动淤泥与沉积淤泥的流变特性研究.doc

水利工程论文扰动淤泥与沉积淤泥的流变特性研究摘要为了分析不同沉积状态的淤泥流变特性,本文将淤泥作为线性粘弹性体,提出了根据泥水系统控制方程解答和波浪在泥床上传播时的实测变化规律求取淤泥流变参数的反分析方法,并利用遗传算法分别求解搅匀扰动与自然沉积淤泥的流变参数.结果表明,淤泥平均密度相同时,自然沉积淤泥的弹性模量明显大于人工搅匀淤泥的弹性模量,两者的粘性系数随淤泥密度的变化规律也有明显区别,这主要是由于沉积淤泥与搅匀淤泥的内部结构不同所造成的.关键词自然沉积淤泥扰动淤泥流变特性波浪传播在河口海岸粘性泥沙运动研究中,淤泥的流变特性和许多物理过程密切相关,如波浪在淤泥质海床上的衰减特性和波浪作用下的泥床质量输移、淤泥液化等1~3,是影响波浪与底泥床相互作用的主要因素.淤泥流变特性的变化和水动力作用下床面稳定性、底部泥沙的悬扬规律等也有着直接和间接联系.实验和计算结果表明,淤泥处于不同的沉积状态时,由于流变参数的变化,波浪在泥床上传播时波高衰减率可以相差1~2个量级4.因此,了解淤泥的流变特性是深入研究淤泥质海岸泥沙运动规律的关键问题之一.关于淤泥的流变关系,目前已建立了相当多的模型4,然而由于模型参数特别是现场条件下的模型参数难以确定,使得各种模型在应用于描述现场条件下波浪和淤泥质海床的相互作用规律时受到限制.由于现有的测量仪器,如旋转同心圆筒流变仪等实际上只能测量扰动淤泥的流变特性,对于自然沉积的底部泥床,目前还没有很好的方法来测定其流变参数,对其流变特性的变化规律有待深入研究.基于以上原因,本文假定淤泥作为线性粘弹性体,在文献5的基础上,采用根据实验结果反求模型参数的方法,比较了实验室内人工搅匀和自然沉积淤泥的流变特性,讨论了淤泥流变特性变化对波浪衰减规律的影响.1理论分析1.1关于淤泥特性的假设近年来一系列的研究表明,淤泥在波浪作用下表现出复杂的非线性粘弹性体特征,其流变参数是应变或应变率历史的函数.但为了简化问题,通过反分析方法确定淤泥的流变参数,并比较人工搅匀淤泥与沉积淤泥流变特性,这里仍然假定淤泥作为线性粘弹性体.在应变较小的情况下,上述假设精确地描述了淤泥的运动规律在应变较大的情况下,这种假设相当于对淤泥的本构方程进行等价线性化.根据线性粘弹性体假设,在频率为σ的循环振荡荷载作用下,淤泥的本构关系可表示为6对于天然或实验室内形成的自然沉积泥床,沿泥床表面向下的泥密度分布及其流变特性都是变化的,为了简化问题而方便估计沉积淤泥的流变特性,我们假定泥床具有均匀密度和流变参数,水波与搅匀或沉积泥床的相互作用就都可以用密度均匀的两层介质线性系统模型来描述.1.2泥水系统运动控制方程根据Maa和Mehta7,波浪在泥床上传播时,上层水体和下层淤泥运动的连续性方程和线性化运动方程可表示为假定水面和泥水交界面位移表达式为ηjx,tajexpikxσt5式中kkriki为待求的复波数,实部与波长L的关系为kr2π/L,虚部ki表示波浪衰减率,aj为第j层位移振幅.方程2~4中uj,wj,pj的解可表示为在波浪周期T、波高HH2a1、水深d1、泥厚d2、泥密度ρ2以及淤泥流变参数已知的条件下,上述表达式中包括k,a2和Aj,Bj,Cj,Djj1,2总共10个未知量,这些未知量可通过自由表面的运动学和动力学边界条件,泥水交界面的速度和应力连续条件以及底部非滑移边界条件来确定7,8.将式5~8代入上述边界条件并经过化简,可得到变量为A1,B1,C1,D1,A2,B2,C2,D2,a2的线性齐次方程组C09式中C为99矩阵,矩阵中各元素为复波数k的函数.式9有非零解的条件为C的行列式必须为零,即detC010式10是关于复波数k的隐式非线性方程.由于水的密度与粘性系数可作为常数,式10相当于复波数k可以表示为kfμ,G,d1,d2,T,ρ211在波浪泥床系统的基本参数和淤泥流变参数已知时,式11可以通过Muller方法求解.1.3淤泥流变参数的反求法在淤泥流变参数确定的情况下,根据方程11可以求出波浪在泥床上传播时的复波数,然后可进一步得到速度和压力分布,这是分析波浪和泥床相互作用问题的一般过程.同样,如果已知波浪在泥床上传播的复波数而淤泥流变参数未知,由方程11也可以反求出淤泥的表观复粘度,从而得到淤泥的弹性系数和粘性系数.按照把淤泥作为线性粘弹性体的假设,对于同一种淤泥,泥密度保持不变时,由不同波浪条件的实验结果反求出的淤泥粘弹性系数值应该相同.但实际上,由于淤泥的非线性性质,其粘弹性特征随受力状态而变化8,9,再加上实验测量误差的限制,根据不同实验条件反求出的淤泥粘弹性系数是不会相同的.为了求得假定淤泥为线性粘弹性体时代表性的粘弹性参数值,我们根据每组淤泥的多组实验结果构造优化问题来反求参数.设对于某一密度的淤泥共进行了M组规则波浪在淤泥床上传播的实验.对于第mm1,,M组实验,已知波浪周期T、波高H、水深d1、水密度ρ1、水的粘性系数μ1、泥厚d2、泥密度ρ2,根据实验测量结果可求得复波数kem.假定粘弹性参数已知,则根据式11可以求得第m组实验条件下的复波数kcm.由于不存在统一的μe满足所有实验条件下m1,,M的式11,我们构造一个极值问题来确定最具代表性的淤泥流变参数值.令式中kcm为μe的函数,则φ最小时的μe值即可认为淤泥作为线性粘弹性体时的代表值.文献5采用有约束的外部惩罚函数方法来反求μe值,存在对初值的依赖性高,最后收敛解难以确定等问题.这里我们采用遗传算法10来求解式12,最终得到的μe值与初值无关.1.4遗传算法简介遗传算法是通过模拟生物进化来求解问题空间解的一种方法,在处理极值优化问题时,特别是目标函数具有局部极值或导数不连续时,非常有效.遗传算法的主要思路如下在求解问题时,遗传算法首先从一定数量的随机解染色体即种群出发,开始进行进化计算,种群中第一代的每一个染色体经过一定的遗传操作,如变异、杂交和选择等,形成第二代种群,上述过程重复进行直至得到预定代数的种群,在最后的种群中,可能有一个染色体占有绝对优势,它就代表了所求问题的优化解.近年来遗传算法在各个学科得到了广泛应用10,我们采用了Houck等建立的MATLAB遗传算法工具箱求解了式12.2实验测量在天津大学海岸工程实验室内的波浪水槽内进行了波浪在搅匀和沉积泥床上的波浪传播实验.实验中的搅匀泥床指的是根据实验所需的淤泥密度,向初始密度比较大的淤泥中配置不同的自来水,进行充分搅拌形成的从泥床表面到泥床底部密度均匀的淤泥床沉积泥床是把淤泥和水进行充分搅拌,形成密度在1120kg/m2左右的均匀泥水混合液,经过自然沉积不同天数而形成的泥床,本次试验中沉积泥床最少沉积天数为4天,最大沉积天数达53天.实验用泥取自天津新港,关于实验用泥的处理、实验布置及测量方法详见文献8.根据水槽内波高传感器的实验结果可以得到波浪在泥床上传播时的波长和波高指数衰减率,由此便可以求得实测复波数ke,从而反求淤泥粘弹性参数.表1列出了搅匀与沉积泥床上的实验条件.表1搅匀与沉积泥床实验条件3搅匀和沉积淤泥流变特性的比较无论搅匀泥床还是自然沉积泥床,当淤泥密度较大时,波浪在床面上的衰减都很小,通过反分析方法难以合理确定淤泥粘弹性参数,所以对于搅匀和自然沉积淤泥,我们分别采用泥密度在1401kg/m3和1336kg/m3以下的数据来反求淤泥粘弹性参数.图1和图2分别显示了把淤泥作为线性粘弹性体时,根据波浪在泥床上传播的实验结果反求出的搅匀泥床和自然沉积泥床的弹性

注意事项

本文(水利工程论文-扰动淤泥与沉积淤泥的流变特性研究.doc)为本站会员(奋斗不息)主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网([email protected]),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。

copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5