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

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

水利工程论文-碾压混凝土坝渗流场分析的缝面渗流平面单元模拟法.doc水利工程论文-碾压混凝土坝渗流场分析的缝面渗流平面单元模拟法.doc -- 2 元

宽屏显示 收藏 分享

资源预览需要最新版本的Flash Player支持。
您尚未安装或版本过低,建议您

水利工程论文碾压混凝土坝渗流场分析的缝面渗流平面单元模拟法摘要针对碾压混凝土坝成层施工的结构特点及其渗流特性,笔者在非均质各向异性等效连续体模型的基础上,提出可精细而又方便地模拟大坝中层面、缝面和裂缝等渗流行为的缝面渗流平面单元模型。在没有增加计算网格结点的前提下,不但能高效率地刻画出层面和缝面所带来的坝体强渗透各向异性的特点,而且又能极容易地对那些局部零星分布的有强集中渗透能力的冷缝和裂缝等缝隙进行逐个细致地模拟,坝体的渗流实际情况在数值模型中得到了甚为客观的展现,提高了大坝渗流场的计算分析精度,基本上完整地解决了碾压混凝土坝的渗流建模问题。最后,较详细地进行了一个高102m的大坝的算例分析。关键词碾压混凝土坝渗流层面缝面裂缝缝面渗流平面单元在碾压混凝土成层结构中,混凝土本体为弱透水介质,而层面和缝面为相对强透水结构面,通常可对碾压混凝土坝体用水力等效连续体模型和非连续体裂隙模型进行渗流特性建模[1~4],其中对冷缝面或裂缝采用空间薄层单元进行模拟。由于正常的施工层面和缝面的间距相对于大坝的特征尺寸很小,相对简单的已有成熟的理论基础和丰富运用经验的等效连续体模型已得到了广泛的推广应用,但这种模型中存在的很大的弱点是对上述那些单个地出现的冷缝面和裂缝有时不能进行很好地模拟,薄层单元的引入在一定程度上不但扩大了解题的规模,而且还会影响解题的精度。笔者针对这个问题和碾压混凝土坝中存在的局部集中渗漏现象较普遍的事实,提出无厚度缝面平面单元,用来专门模拟坝体中的冷缝面和裂缝等集中渗流通道,并融入等效连续体模型中,从而达到既整体又突出局部地刻画和模拟碾压混凝土坝的渗流行为和高精度地进行渗流场的计算分析。可在没有增加网格结点的前提下,能高效率地反映出碾压混凝土坝中正常施工层面和缝面的渗流影响,而且又能对那些具有强透水性零星分布的冷缝面或结构性裂缝的渗流行为进行单独细致的模拟。1基本理论1.1渗流基本理论非均质各向异性多孔隙连续体介质中的稳定饱和渗流连续微分控制方程为1式中xi为坐标,i1,2,3kij为达西渗透系数张量,刻画坝体的渗透各向异性程度hx3p/γ为总水头,x3为位置水头,p/γ为压力水头。边界条件为2式中h1为已知水头函数ni为渗流边界面外法线方向余弦i1,2,3Γ1、Γ2、Γ3和Γ4分别为第一类和第二类渗流边界,以及渗流自由面和渗流逸出面qn为法向流量,流出为正。1.2立方定律若将有集中渗流能力的缝面和裂缝视为具有一定水力隙宽的缝隙结构面,并认为结构面中的水流流态为层流,则可以用层流缝隙流立方定律[5]来描述结构面的渗流行为,即3qνdfkfdfI4式中ν为结构面中的平均流速df为结构面的水力等效隙宽I为结构面中的水力梯度q为结构面中的单宽流量kf为水力等效渗透系数。1.3缝面平面单元模型因碾压混凝土本体的透水能力很小而缝面冷缝和裂缝的透水能力与其水力隙宽的立立成正比,缝面无论在其切向或法向上的透水能力均为相对很大,尤其是在缝面法向方向上其渗透系数为一相对大值,通常缝面需单独划分成具有一定厚度的薄层单元。但是缝面的隙宽一般很小,只有几十μm,甚至更小,单元在缝面法向方向的尺寸相对于缝面切向方向的尺度甚小,这种薄层单元在理论上精度差,有时甚至会严重降低整个渗流场的求解精度。另外,薄层单元需单独划分出来,增加了解题的规模。这里提出缝面缝隙渗流的无厚度缝面平面单元,可用来专门解决具有集中渗流能力的零散性缝面的渗流模拟问题。因缝面的水力隙宽很小,法向的透水能力又极大以及混凝土本体透水能力极小,对于工程渗流问题而言,完全可以认为在缝面内缝面法向的水头损失为零,缝面中的水流呈准二维渗流状态,因此有5式中为与缝面相关的局部坐标为缝面平面单元的二维渗透系数张量,反映缝面的渗透各向异性和缝面的透水能力,其中可借助上述立方定律缝面的水力隙宽df也被考虑在中。局部坐标与整体坐标xi之间的关系为垂直于缝面6式中、[Te]和{xi}分别为局部坐标向量,缝面平面单元的坐标转换矩阵和整体坐标向量。任何一个缝面平面单元e的结点局部坐标得到确定后,就可进行这个单元的传导矩阵元素的计算7式中sf为缝面单元域Nr和Ns为缝面单元的插值函数m为缝面单元的结点数。这种缝隙渗流缝面平面单元对坝体内渗流场特性的影响是通过下式的水头连续条件来实现的。8即缝面sf任一点处碾压混凝土壁面上的水头与无厚度缝面平面单元同一位置处的水头是相同的。再按通常的有限单元法要求,据计算域中任一结点处的流量平衡条件[7]式12就可以组装成常规形式的求解整个渗流场的有限单元法支配方程。9式中n为总结点数QRCC和Qf分别为由碾压混凝土坝三维等效连续体单元和二维缝面平面单元对结点i所作贡献的等效结点流量式12须对环绕i结点的所有常规单元和缝面单元求和。用上述缝面平面单元模拟冷缝面或裂缝渗流行为的优点有1二维缝面平面单元没有厚度,网格剖分时无需对它进行专门剖分,只需额外形成一个简单的所有这类二维单元的单元信息表即可,不增加额外结点和自由度,几乎不增加网格剖分和求解时的计算工作量2避开了在缝面法向隙宽方向上的微分及积分运算,不因对缝面或裂缝渗流行为的处理而给整个渗流场求解带来大的误差3能极方便地正确模拟碾压混凝土坝中局部零星散落的裂缝、层面和缝面的集中渗流行为4此方法的引入,碾压混凝土坝的渗流数值建模问题已得到了基本完整的解决5能考虑缝面的渗流各向异性。2碾压混凝土坝渗流场分析的有限单元法求解2.1渗流场求解的有限单元法就固定网格求解有自由面的渗流问题的有限单元法的结点虚流量法可详见文献[8,9]。此外,在实际工程中某些渗流控制措施、边界条件和渗透系数张量的确定也往往是求解精度的主要影响因素,需分别专门对待。因混凝土大坝渗流场的水头分布主要受控于排水措施的作用,解题时必须对计算域中的所有排水措施比如排水孔、排水洞等进行精细的模拟。2.2排水孔渗流开关器排水孔排水幕的渗流行为可严密地用排水子结构来处理,但不同类型的排水孔在渗流边界条件上的不同表现导致程序中对它们处理方式的不同。因工程渗流问题的复杂性,事先不知道坝基每一个溢流型排水孔的顶面高程一定是低于或高于邻域自由面的高程,就象是确定逸流型排水孔内渗流逸出线位置那样,也需在迭代求解过程中据中间解逐步给予确认,为此需引进排水孔开关器的概念[10]。在每一个溢流型排水孔顶面虚构一个数学开关器,当开关器打开时排水孔工作,孔壁面为一个等水头面当开关器关闭时排水孔完全不工作,即此时孔顶面高程高于邻域自由面位置,排水孔被自由面穿过。这样解题就可以确保每一个溢流型和逸流型排水孔的真实工作状态都得到正确的甄别和模拟。3计算分析3.1概况某碾压混凝土重力坝,坝高102m,底宽70m,顶宽12m,上游水位100m,下游水位20m,在坝体中距上游面23m处设置孔径15cm,孔距4m的垂直排水幕,坝基设置2道深分别为30m、15m的主排水幕,中间2道深为8m的抽排水幕,其孔口高程为基础廊道底版高程,孔距、孔径都与坝体排水孔相同。假定坝体内自下而上每隔20m存在一个冷缝面,共4个,坝上、下游面分别设置03m厚的变态混凝土,上游面变态混凝土后有厚2m的二级配混凝土,建基面上设有2m厚的常态混凝土垫层。计算中当考虑防渗体水平开裂时,裂缝都贯穿于坝上游面的变态混凝土和其后的二级配碾压混凝土坝上游面垂直劈头缝位置假定发生在两相邻排水孔的中间。图1为大坝的剖面图,利用水力对称性计算域在坝轴线向厚2m。计算域材料的渗透系数见表1所示,其中k∥和k⊥分别为坝体层面切向和法向的主渗透系数。表1计算域材料的渗透系数单位cm/s渗透各向性变态混凝土二级配混凝土三级配混凝土常态混凝土坝基帷幕k1109k∥1107k⊥=1109k∥1106k⊥=1109k5109k1104Δ60m以上k1105Δ60m以下k11053.2计算工况这里主要研究坝体变态混凝土防渗体的水平向或竖直向的开裂及4个典型集中渗水缝面不同水力隙宽时坝体渗流场的变化情况。具体的计算工况见表2所示,裂缝均贯穿二级配混凝土。表2计算工况工况号结果图号工况说明12345图2图3图4图5图6上游面防渗体无裂缝和坝体无缝面上游面防渗体有隙宽0.2mm的水平裂缝,坝体有隙宽0.01mm缝面但无排水幕上游面防渗体有隙宽0.2mm的水平裂缝,坝体有隙宽1mm缝面且有排水幕上游面防渗体有隙宽0.2mm的垂直劈头裂缝,坝体有隙宽1mm缝面但无排水幕上游面防渗体有隙宽0.2mm的垂直劈头裂缝,坝体有隙宽1mm缝面且有排水幕3.3计算结果分析图2渗流场水头线分布的主要特点有1坝上游面变态混凝土防渗体起到一定的防渗作用,但仅占坝上、下游面总水头的10~202尽管二级配碾压混凝土的层面切向渗透系数比三级配碾压混凝土的小一个数量级,但防渗效果并不明显3坝体渗透各向异性比使得等水头线向下游折弯,因三级配碾压混凝土各向异性比达3个数量级,在三级配碾压混凝土区等水头线就更呈现得水平向状态,水平向的渗透水力阻力很小4下游面变态混凝土的阻渗作用非常明显,若没有它等水头线会呈现几乎完全水平的状态,渗透水流会直接从坝下游面逸出,逸出线位置会很高5坝下游面有变态混凝土后,在实际工程中坝内排水设施务必要布置得充分一些6因变态混凝土很不透水,坝下游面有蒸发作用,一般坝下游面不会留有渗流逸出水迹现象,改善了大坝运行期的观瞻性。图1碾压混凝土坝剖面及渗流场计算域图2坝体无裂缝和冷缝面时等水头线分布单位m图3与图2的结果对比,尽管缝面的水力隙宽只0.01mm,亦即赤眼是完全看不到缝面的存在,但是由于碾压混凝土本体的渗透性极小,这样的缝面已明显地形成了坝中的主要集中渗流
编号:201312122224246401    大小:52.37KB    格式:DOC    上传时间:2013-12-12
  【编辑】
2
关 键 词:
专业文献 建筑水利 精品文档 水利工程
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 人人文库网仅提供交流平台,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

当前资源信息

4.0
 
(2人评价)
浏览:6次
奋斗不息上传于2013-12-12

官方联系方式

客服手机:13961746681   
2:不支持迅雷下载,请使用浏览器下载   
3:不支持QQ浏览器下载,请用其他浏览器   
4:下载后的文档和图纸-无水印   
5:文档经过压缩,下载后原文更清晰   

相关资源

相关资源

相关搜索

专业文献   建筑水利   精品文档   水利工程  
关于我们 - 网站声明 - 网站地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5