ABAQUS程序在某混合坝不同坝型接触面渗流分析中的应用.docxVIP

第３２卷第２期水 电 能 源 科 学Ｖｏｌ．３２Ｎｏ．２２ ０ １ ４年２月Ｗａｔｅｒ Ｒｅｓｏｕｒｃｅｓ ａｎｄ ＰｏｗｅｒＦｅｂ．２ ０ １ ４文章编号：（ ）１０００－７７０９２０１４０２－００８６－０４ＡＢＡＱＵＳ程序在某混合坝不同坝型 接触面渗流分析中的应用冯龙龙１，２，３，包腾飞１，２，３，李晓晨１，２，３（１．河海大学 水文水资源与水利工程科学国家重点实验室，江苏 南京 ２１００９８；２．河海大学 水资源高效利用 与工程安全国家工程研究中心，江苏 南京２１００９８；３．河海大学 水利水电学院，江苏 南京２１００９８）摘要：针对混合坝坝体边界较复杂致使边界条件施加困难，从而增加渗流计算难度的问题，利用大型通用有，，限元软件ＡＢＡＱＵＳ提供的用户子程序 对ＡＢＡＱＵＳ进行二次开发 对某混合坝不同坝型接触面进行渗流计，。，算时 采用Ｆｏｒｔｒａｎ子程序对其施加孔压边界条件并进行渗流场分析结果表明 该混合坝的防渗安全满足要求，进而证明提出利用二次开发施加边界条件可行。关键词：渗流分析；ＡＢＡＱＵＳ程序；二次开发；三维混合坝；防渗安全中图分类号：ＴＶ１３９．１４文献标志码：Ａ１ 引言２ 渗流计算方程渗流在水利工程中普遍存在，渗流问题直接达西定律表达式为：关系着水利工程的投资与安全。许多水利工程的（）ｖｉ ＝－ｋｉｈ，ｉ１失事均与渗流有关，如年鲍德温山坝由于与式中， （，，）为三个渗透主轴方向的渗流１９６４ｖｉ ｉ＝１２３基础接触面发生渗流破坏而失事；年堤堂坝速度； （，，）为三个渗透主轴方向的渗透１９７６ｋｉ ｉ＝１２３由于右岸一个窄断层发生渗透破坏而垮坝［１］。正系数； （，，）为渗流场中各点的测压管水头。ｈｉｉ＝１２３确模拟和分析渗流问题对水工建筑物的设计和安稳定渗流场连续性方程为：全运行尤为重要，因此水利工程中尤为重视渗流（）ｖｉ，ｉ ＝０２问题。，、、将式（）代入式（）中，即可获得稳定渗流的目前 渗流分析常用ＡＮＳＹＳＡＤＩＮＡ１２、，、、微分方程［２］的一般表达式：ＭＡＲＣＡＢＡＱＵＳ等软件 但ＡＮＳＹＳＡＤＩＮＡ３，ＭＡＲＣ无专门的渗流计算模块 需通过温度场模拟渗流场。相比而言，ＡＢＡＱＵＳ软件则有渗流（ ）（）∑ｋｉｈ，ｉ ，ｉ＝０３ｉ＝１计算模块，提供了用于渗流计算的孔压单元，且能（ （））将渗流基本微分方程 式 ２与定解条件相较好地适应复杂边界和复杂渗流介质的空间渗流结合即可确定渗流数学模型。渗流计算过程中控，。但 ＡＢＡＱＵＳ在渗制方程常用的定解条件有：①第一类边界条件（水问题 使渗流计算变得简便流计算时还有一些不足，如在一定情况下边界条头边界条件或Ｄｉｒｉｃｈｌｅｔ条件）：件、初始条件的施加较繁琐。因此，本文结合（，，） （，，，）（）ｈｘｙｚ ＝ｈ０ ｘｙｚｔ（ｘ，，ｚ）４ｙ ∈Γ１， （，，）；为已知水，式中 ｈ０ ｘｙｚ 为已知水头函数Γ１ＡＢＡＱＵＳ提供的二次开发平台 通过子程序， ，，。②第二类边界条件，头边界 ｘｙｚ位于Γ１ 上ＤＩＳＰ实现孔压边界条件的施加 进而对某混合坝不同坝型接触面进行渗流分析，结果表明可满足（）：流量边界条件或 Ｎｅｕｍａｎｎ条件防渗安全需要。（，，，）（）－ｋｈｉ，ｉｎｉ ＝ｑ０ ｘｙｚｔ（ｘ，，ｚ）５ｙ ∈Γ２收稿日期：２０１３－０６－０３，修回日期：２０１３－０７－３１基金项目：国家自然科学基金项目（）；新世纪优秀人才支持计划项目（）；高等学校博士学ＣＥＴ－１１－０６２８（）；（）点专项科研项目 ２０１２００９４１１０００５中央高校基本科研业务费项目 ２０１２Ｂ０７２１４作者简介：冯龙龙（），男，硕士研究生，研究方向为大坝安全监测， ：１９９０－Ｅ－ｍａｉｌｆｅｎｇｌｏｎｇｌｏｎｇ５＠１６３．ｃｏｍ通讯作者：包腾飞（），男，教授、博导，研究方向为大坝安全监控理论及分析， ：１９７４－Ｅ－ｍａｉｌｂａｏｔｆ＠ｈｈｕ．ｅｄｕ．ｃｎ第卷第期冯龙龙等：程序在某混合坝不同坝型接触面渗流分析中的应用··３２２ＡＢＡＱＵＳ８７式中， （，，，）为已知流量函数；为具有高，心墙堆石坝部分最大坝高ｑ０ｘｙｚｔΓ２１５９．００ ｍ７１．００给定流量的边界，，， 位于上；为的外。坝顶总长，其中混凝土坝部分长ｘｙｚΓ２ｎｉΓ２ｍ１ １５８．００ｍ法线方向。，，８３８．０３ｍ心墙堆石坝部分长３１９．９７ｍ正常蓄３ 稳定渗流分析在 ＡＢＡＱＵＳ中的水位１ １３４．００ｍ。混凝土坝与心墙堆石坝之间采用插入式（包实现裹式）软接头布置形式。心墙堆石坝的防渗体为土质心墙，采用自上而下逐渐加厚的形式。心墙　　多孔介质渗流分析的关键问题是要确定自由顶宽６ｍ，两侧坡比１∶０．２。为保护心墙不发生面（浸润面）和逸出点的位置。目前，求解自由面渗透变形，在心墙上下游侧均布置了两层反滤层，的