城区区域性渗流场水平面有限元简化分析方法.pdfVIP

城区区域性渗流场水平面有限元简化分析方法 ， 骆冠勇 曹洪 (广州华南理工大学土木工程系，广州510640) 摘要：城区区域性渗流场，混凝土阻水结构多，透水缝隙多，渗流场的大尺度与阻水结构的小尺寸相差好几个 数量级，用一般水平面有限元分析存在较大的困难。针对城区渗流场的特点，本文提出相应的简化方法，并分 析了各种方法的精度。对于混凝土阻水结构，结合无限相似单元和子结构法各自优点利用奇异点子结构法简化 对阻水结构的模拟；对于悬挂式阻水结构体，利用解析解和阻力系数法，引入无厚度界面接触单元考虑悬挂式 阻水块体周边冈竖向流动产生的局部水头损失，以弥补水平面渗流有限元不能考虑竖向流动产生的水头损失的 缺点；对于城区大区域渗流场中建筑物间的透水缝隙，同样利用解析解和阻力系数法，引入修正的线单元对缝 隙进行模拟。利用上述方法编制相关程序，设计了相关试验对悬挂式阻水结构模拟简化方法和建筑物间缝隙模 拟的简化方法进行验证，实测结果与计算结果相互吻合。这些方法为城市大区域渗流场的分析提供了有效途径， 在不失精度的前提下，基于应用水平面有限元进行计算，在整体上用尽量少的网格实现对城市大区域渗流场的 模拟计算。 关键词：城市地下水 区域性渗流数值分析有限元 1 前言 近20年来，广州城区地下空间发展迅速，地下混凝土结构也越来越多。广州地质条件的一个主要特点是： 作为相对不透水层的基岩埋藏较浅，一般仅为10。20m。这种情况下，地下结构直接落在基岩面上或与基岩 相距很近，堵塞了该断面上的大部分渗流通道，当这些地下结构铺成片排成行时对地下水环境的影响不容忽 视。相关的研究已经开始受到重视，但定量分析工作开展不多。 城区渗流场，范围大，几十甚至上百平方公里，是个区域性渗流场。目前，区域性渗流场计算常用的方 法有水平面有限元法、有限差分法和解析元素法…，其中水平面有限元法以其灵活边界适应性而被广泛应用。 但是城区渗流场除计算区域大之外，还具有自身的特点：一是建筑物阻水结构多；二是建筑物间的透水缝隙 多。阻水结构单体的尺寸一般为几千平方米，与流场的尺寸相比相差几个数量级，建筑物间缝隙更比建筑物 本身小得多。而且阻水结构物和缝隙出入口处的角点均为奇异点，在该处的网格必须加密才能得到较为准确 的结果。另一方面，这些阻水结构中有些是悬挂式，周边存在较强烈竖向流动，一般的水平面渗流有限元因 遵循裘布依假定忽略了竖向流动的水头损失，用它进行计算存在较大的误差，用三维有限元计算区域大而显 得十分困难。上述情况导致城区渗流场有限元计算的工作非常大，有必要进行简化计算。 针对城区区域性渗流场的特点，基于应用水平面有限元进行计算，提出相应的简化方法：对于阻水结构， 结合无限相似单元【21和子结构法【31各自优点利用奇异点子结构法简化对阻水结构的模拟，以实现用大网格对 小结构进行划分；对于悬挂式阻水结构体，利用解析解和等效阻力系数法，引入无厚度界面接触单元考虑悬 挂式阻水块体周边因竖向流动产生的局部水头损失，以弥补水平面渗流有限元不能考虑竖向流动产生的水头 损失的缺点，将三维问题简化为二维问题进行处理；对于城区大区域渗流场中建筑物间的透水缝隙，同样利用 解析解和等效阻力系数法，引入修正的线单元对缝隙进行模拟，将小块体合并成大块体进行计算。上述三种方 法都旨在不失精度的前提下，在整体上用尽量少的网格实现对城市大区域渗流场的简化计算。文献[4—6]对这 些方法做了详细的论述，本文工作是简要介绍这些方法和精度，然后再利用相关试验对后面两种方法进行验证。 2 方法介绍 2．1 城区渗流场中阻水块体有限元简化分析方法 由于城区区域性渗流场范围大、建筑物阻水块体多，奇异角点多，渗流场的大区域与阻水结构的小尺寸 第一作者简介：骆冠勇，1979年生，男，博士研究生，主要从事渗流与岩土工程研究，E．mail：luogyong@126．tom ·168· 之间相差几个数量级，有限元分析存在较大的困难，引入应用奇异点子结构法对渗流场中的阻水结构周边的 奇异点进行处理，使得大区域渗流场中的每一个建筑物阻水块体在整体网格划分时其网格尺寸取为其边长的 1／2-1／3即可(如图1)，然后对阻水结构的每一个角点作为奇异点子结构进行处理。 ^更一《槭酱}靛霉