首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 140 毫秒
1.
基于连续多孔介质原理与混合体原理引入流?固耦合模型,采用FISH语言编制程序,建立了初始渗流场设置函数、饱和?非饱和渗流分析模块、特殊应力修正模块以及非饱和单元抗剪强度修正模块,实现了利用FLAC3D进行饱和?非饱和渗流分析,同时基于饱和?非饱和渗流场修正非饱和区土体单元的有效应力以及抗剪强度,完成了将FLAC3D中饱和土流?固耦合计算原理扩展到非饱和土中。在二次开发的分析模块基础上对Liakopoulos砂柱排水试验进行数值模拟,验证了计算模型和计算程序的正确性。同时也研究了降雨入渗对土质边坡的渗流场、位移场以及稳定性演变过程的影响,揭示了降雨诱发边坡失稳的机制。  相似文献   

2.
三维饱和-非饱和渗流场求解及其逸出面边界条件处理   总被引:24,自引:1,他引:24       下载免费PDF全文
对三维饱和-非饱和渗流场的有限单元法求解进行了详细论述,在边界条件中明确列出了非饱和渗流逸出面边界条件的数学表达式。为了理论上能够严密地给出非饱和逸出面的边界条件和数值处理,根据渗流场固有的物理特性,定义了极限含水率和非饱和渗流逸出面蒸逸系数两个渗透基本特性参数。针对渗流场饱和与非饱和逸出面分界线(逸出线)位置事先不知道的问题,提出在这种边界非线性渗流场迭代求解过程中,应根据中间解通过边界条件类型的转化来准确地甄别这一分界线,完整地获得精细地处理饱和及非饱和渗流逸出面边界条件的思路与方法。  相似文献   

3.
提出了求解有自由面渗流问题的三维数值流形方法,通过构造任意形状流形单元的水头函数,推导了流形单元的渗透矩阵和无压渗流分析的总体控制方程,并给出了自由面的迭代求解策略和渗透体积力的计算方法。典型算例的数值分析表明,该方法采用数学网格覆盖整个材料区域,在自由面的迭代求解过程中数学网格保持不变,只考虑自由面以下渗流区的介质,只对自由面以下的流形单元形成总体渗透矩阵,具有精度高、收敛速度快、编程简单等优点,而且能够通过单纯形积分精确计算被自由面穿越单元的渗透作用力,因此,特别适用于有自由面渗流问题的模拟。  相似文献   

4.
大型地下洞室群围岩应力-损伤-渗流耦合分析   总被引:5,自引:2,他引:3  
张巍  肖明  范国邦 《岩土力学》2008,29(7):1813-1818
以渗透体力来考虑渗流场的力学效应,建立了应力-损伤-渗透系数关系方程来考虑应力和损伤对渗流场的影响,结合岩体结构的三维弹塑性损伤有限元分析,建立了大型地下洞室开挖围岩应力-损伤-渗流耦合的计算模型。该计算模型求解的难度主要体现在岩体材料弹塑性、损伤、渗流自由面边界、渗流溢出边界、应力-损伤-渗流相互影响关系等。提出分步迭代法对以上因素进行归纳后按一定顺序分别进行迭代求解,取得了良好的计算效果。将该方法应用于某水电站大型地下洞室群分析,得出了一系列有意义的结论。  相似文献   

5.
在对饱和-非饱和渗流数学模型简要分析的基础上,对影响渗流有限元计算分析结果的一些参数选取方法进行讨论,如非饱和土的水土特征曲线、渗透系数及毛细水头等。利用GeoStud io软件对降雨条件下的土质边坡内部渗流场进行模拟,分析地下水位线、坡体饱和度、压力水头随降雨时间的变化规律。模拟结果与实际情况比较吻合,可为土质边坡稳定性的分析提供参考依据。  相似文献   

6.
基于FLAC3D平台的边坡非饱和降雨入渗分析   总被引:7,自引:0,他引:7  
蒋中明  熊小虎  曾铃 《岩土力学》2014,35(3):855-861
为研究降雨入渗条件下三维边坡渗流场的变化过程,在全面研究FLAC3D软件渗流分析模块功能及算法的基础上,通过编写FISH函数完善了FLAC3D软件的非饱和渗流计算功能。基于降雨入渗机制的分析,编写了降雨入渗及停雨出渗边界的FISH函数,实现边坡三维降雨入渗过程的模拟。通过降雨条件下算例边坡的饱和非饱和渗流场变化过程的数值计算,研究了降雨入渗引起的边坡暂态饱和区变化过程,并与已有研究成果进行对比分析,从而验证了自编FISH函数实现边坡三维非饱和渗流计算结果的正确性。研究成果表明,文中提出的边坡三维非饱和降雨入渗分析是可行的,同时也为三维边坡非饱和渗流场的研究提供了一条有益的途径。  相似文献   

7.
考虑饱和-非饱和渗流的土坡极限分析   总被引:2,自引:1,他引:1  
给出了均质土坝的下游坝坡的安全系数的计算方法。均质土坝的下游坝坡有可能会通过滑裂面发生破坏。土坝中的土坡通常处于非饱和状态。非饱和土坡的安全系数计算需要考虑吸力对抗剪强度的贡献以及土坡中的非饱和渗流。给出处于饱和-非饱和渗流状态下的土坡的安全系数,有助于评价均质土坝的安全系数。下游坝坡的安全系数计算方法有:极限平衡法、上限解法和下限解法,该算法适用于非饱和土坡,且是在饱和土坡安全系数的计算方法上修正得到的。算例中非饱和土坡安全系数的计算考虑了吸力对抗剪强度的贡献。考虑非饱和渗流理论的土坡安全系数计算方法通常更加接近现场实际情况,并且对于同一坝坡,考虑非饱和渗流计算出的土坡安全系数要比饱和渗流理论计算出的安全系数大。  相似文献   

8.
探索考虑饱和-非饱和渗流场和应力场耦合的三维强度折减有限元技术,并研制成功一个考虑饱和-非饱和渗流场和应力场耦合的三维强度折减有限元程序,通过与传统极限平衡法分析结果对比研究,对抗剪强度折减有限元法分析边坡稳定问题的适用性进行了评价,得出采用三维强度折减有限元法确定考虑饱和-非饱和渗流场的边坡稳定性安全系数是可行的结论.  相似文献   

9.
在三维非恒定渗流有限元计算中,不可避免地需计算自由面上的边界积分项。建议一种高精度数值积分的方法求取自由面边界积分项。其基本思路是,基于8节点空间等参单元,根据压力为零的边界条件确定自由面满足的曲面方程,将自由面边界积分转化为ξη平面上的二重积分;然后再计算自由面与三维等参单元的交点,将交点投影到局部坐标平面ξη上,并根据点与线的相对关系,确定投影点所围成的局部坐标积分区域;再将总积分区域划分成若干个三角形子区域,并利用变步长Simpson方法计算各三角形子区域上的二重积分,从而实现了高精度的自由面边界积分。该方法避免了单元中自由面为平面的假设,可提高计算精度,特别对于单元内自由面变化剧烈情况,更为显著。将该方法应用于砂槽模型和土坝的渗流分析中,计算结果与实际情况较为接近。  相似文献   

10.
基于Ansys的渗流自由面计算方法   总被引:2,自引:0,他引:2  
陈江 《岩土工程技术》2008,22(5):236-240
对渗流自由面计算的单元传导矩阵调整法进行了改进。在局部坐标系下,将单元离散为一系列面元(三维为体元),根据面元或体元形心处的压力水头值来确定面元或体元的权系数,然后按加权平均确定单元的等效渗透系数,该方法既可以避免复合单元中自由面上下面积或体积的计算,面元或体元的离散密度也可以任意选择,简化了计算程序。通过预设的收敛参数来控制迭代过程,可避免相邻两次计算的渗透系数变化太大而难于收敛。以Ansys为开发平台,用APDL语言编写了相应的计算程序,扩展了Ansys的计算功能,为渗流自由面计算提供了新途径。算例表明:该方法的收敛性较好、计算精度较高,能满足工程需要。  相似文献   

11.
李杰  李强 《地质力学学报》2001,7(3):271-277
在水电、石油、交通等部门的地下工程中,经常会遇到地下水的渗流问题,如果处理不当,可能诱发应力恶化而影响结构稳定。本文以非饱和介质的渗流模型为基础,实现了地下洞室开挖过程中的非稳定渗流场解析。
按分期施工方案,在保持渗流场变化连续性的前提下,对各施工阶段进行了准稳定渗流分析,并在此基础上建立了地下工程非稳定渗流数值模型。其中,对于洞室系统,基于水量平衡原理,采用了子结构模型;对于防渗帷幕,采用了隔水层模型;对于排水井列,采用了杆单元模型。
根据上述模型,采用有限差分法作为数值模拟方法,编写了地下工程非稳定渗流场分析的有限元程序,实现了复杂渗流区域的渗流场的动态模拟。工程实例计算分析表明,用该模型计算的结果是合理的。   相似文献   

12.
地下工程复杂渗流场数值模拟研究与应用   总被引:3,自引:0,他引:3  
结合地下工程复杂渗流场三维有限元分析的实际 ,用渗透张量模拟非均匀各向异性 ,用改进丢单元法迭代自由水面 ,用隐性的复合材料单元模拟断层 ,用隔水层模型模拟防渗帷幕 ,用以沟代井模拟排水井列。在此理论基础上编制了渗流分析的三维有限元计算程序。且用此程序成功求解了漫湾水电站二期工程地下厂房围岩的复杂渗流场 ,并在此基础上分析渗流场分布的主要影响因素的影响。结果验证了程序的合理性和实用性  相似文献   

13.
侯晓萍  徐青  陈胜宏 《岩土力学》2015,36(8):2345-2351
依据空气单元模拟排水孔的思路,在渗流场任意可能逸出边界外布置一层渗透性较强的单元,与其他实体单元一起进行渗流计算,可以实现该边界的自动求解。与现行方法相比,不仅省去了逸出边界迭代计算的麻烦,而且避免了由于逸出点定位不准确而可能引起的渗流计算不合理或收敛困难等缺陷。算例分析表明,空气单元法的计算精度与空气单元的相对渗透系数R有关,与厚度L无关。当选取了合适的R值后,可以较好地逼近现行方法。对于因引入空气单元而造成的强弱透水介质界面处数值拟真性差的现象,还提出了加入接触面薄层单元的改进方法。方法改进后,计算精度可进一步提高。最后,以小湾水电站22#坝段渗流场计算为例,说明该法的实用性和可行性。  相似文献   

14.
降雨条件下软岩边坡渗流-软化分析方法及其灾变机制   总被引:2,自引:0,他引:2  
谢瑾荣  周翠英  程晔 《岩土力学》2014,35(1):197-203
降雨条件下软岩边坡的灾变问题是滑坡领域研究的热点与难点。采用间接耦合法进行软岩边坡降雨入渗分析,根据软岩边坡非饱和渗流-应力计算原理,建立了软岩边坡渗流效应分析方法;将软岩抗剪强度的遇水软化特征动态赋予暂态饱和区中对应的软岩体,模拟降雨渗流条件下软岩边坡的软化效应;将渗流与软化相结合,建立基于暂态饱和区的软岩边坡降雨-渗流-软化-灾变的数值分析方法,并用于实例边坡的降雨灾变机制分析。结果显示,软岩边坡潜在滑动面深度随时间逐渐变浅,安全系数-时程曲线呈反向复曲线形态,当边坡开始软化时曲线上出现拐点。边坡稳定性计算结果与实际情况基本相符,表明思路与方法具有一定正确性,可为软岩边坡治理方案设计提供理论依据。  相似文献   

15.
拼缝渗漏水是盾构隧道最常见的病害表现形式,对衬砌的服役性能有重要影响。在前人研究成果的基础上,提出基于防水密封垫失效的拼缝有效开度假定,进而推导出基于拼缝有效开度的盾构隧道渗流量及衬砌等效渗透系数的实用计算公式。以国内较为通用的常规断面盾构隧道为例,利用上述实用公式估算其渗流量,并与基于块体离散元的数值模拟结果相互验证。再以3种不同几何形态的盾构隧道为原型,完成了135种不同工况下的渗流计算,探讨了有效开度、衬砌外水压及衬砌几何形态,对盾构隧道渗流量及衬砌等效渗透系数的影响。结果表明:拼缝有效开度对渗流量及等效渗透系数的影响十分显著,未考虑拼缝间密封垫失效的传统计算方法可能高估了盾构隧道衬砌的防水性能。本文所述实用计算公式,具有表达清晰、计算简便的特点,可为盾构隧道衬砌的防水设计及后续服役性能评价提供定性和定量的参考。  相似文献   

16.
地下工程渗流排水孔数值模拟的隐式复合单元法   总被引:1,自引:0,他引:1  
倪绍虎  肖明 《岩土力学》2008,29(6):1659-1664
根据渗流基本原理,提出了排水孔在地下工程中渗流数值模拟的隐式复合单元法。将排水孔隐含于模型网格中,得到排水孔在模型中的数据信息。判断所有排水孔与单元的相对位置,通过修正排水孔所穿过单元的渗透传导矩阵来模拟排水孔的强排水效果,并通过实例对排水子结构法和隐式复合单元法进行了比较。运用编制的三维有限元计算程序对某地下水电站厂房渗流场及排水孔进行了模拟,结果表明,采用隐式复合单元法对排水孔进行模拟是可行的。  相似文献   

17.
许桂生  陈胜宏 《岩土力学》2005,26(5):745-749
针对无压渗流问题,在传导矩阵调整法的基础上,对穿越自由面的单元用复合单元来处理,每个复合单元具有两套结点水头,一套用于饱和区的水头插值;另一套则用于无水区。通过变分原理,推出了控制方程并在程序中实现。应用复合单元法可以较好地解决自由面的插值问题,消除了穿过自由面单元的水上部分对自由面附近结点水头的贡献,同时对逸出面边界条件的处理也更加方便和准确。算例表明了该方法的有效性。  相似文献   

18.
土质边坡非饱和渗流场与应力场耦合数值分析   总被引:2,自引:0,他引:2  
阐述了坡面径流-非饱和渗流分析与应力场耦合分析的重要意义和作用;介绍了在非饱和渗流分析中具有重要作用的土-水特征曲线,包括考虑固结压力条件下的土-水特征曲线,并介绍了部分试验成果以及坡面径流-非饱和渗流耦合分析的计算方法;提出坡面径流-非饱和渗流分析与应力场的耦合计算方法,并编制相应的有限元程序;通过算例对比了耦合与非耦合情况下应力场、渗流场的差异,说明耦合分析更符合实际,且所提出的方法是切实可行的。  相似文献   

19.
渗流分析中排水孔模拟的叠单元法   总被引:2,自引:0,他引:2  
排水孔是重要的渗控措施,如何准确模拟排水孔的作用是渗流分析中的关键问题之一。目前普遍采用排水子结构法,但由于子结构网格依存于整体网格,在解决复杂工程问题时网格划分难度很大。文中提出一种新的排水孔模拟方法--叠单元法。对于含排水孔的结构,将其划分为不包含排水孔的整体网格以及各排水孔附近区域的局部网格,网格之间相互独立。然后,通过排水孔局部网格的外边界面的虚拟渗流,实现整体网格与排水孔局部网格的耦合。算例分析的结果表明,当排水孔局部网格范围大小取为整体网格单元尺寸的1~3倍,外边界面渗透系数取为材料渗透系数的102~103倍时,所提方法能够准确模拟排水孔的作用。由于叠单元法的网格划分简便易行,有望广泛应用于复杂工程问题的渗流分析  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号