非连续变形分析方法的改进及应用

首先,采用位移收敛判据对非连续变形的岩石断裂分析进行改进并程序实现.其次,将非连续变形分析方法和裂隙网络渗流模型结合起来,建立了渗流应力耦合分析模型,研究了裂隙岩体变形对渗流的影响以及渗流与应力耦合作用下裂隙岩体的变形破坏特征.研究发现,主干裂隙控制着渗流场分布,考虑耦合作用时的水头值比不考虑耦合作用时大,说明考虑耦合后坡体更容易变形失稳.改进前后裂纹扩展形态基本不变,表明程序计算的收敛判据是有效的.     

非连续变形分析方法的改进及应用1）

首先,采用位移收敛判据对非连续变形的岩石断裂分析进行改进并程序实现. 其次,将非连续变形分析方法和裂隙网络渗流模型结合起来,建立了渗流应力耦合分析模型,研究了裂隙岩体变形对渗流的影响以及渗流与应力耦合作用下裂隙岩体的变形破坏特征. 研究发现,主干裂隙控制着渗流场分布,考虑耦合作用时的水头值比不考虑耦合作用时大,说明考虑耦合后坡体更容易变形失稳. 改进前后裂纹扩展形态基本不变,表明程序计算的收敛判据是有效的.     

脆性岩石破裂过程损伤与渗流耦合数值模型研究

大量的实验结果表明，脆性岩石的渗透性不是一个常量，而是应力和应力诱发损伤破裂的函数．建立了一个描述非均匀岩石渗流-应力-损伤耦合数学模型(FSD Model)，开发出岩石破裂过程渗流-应力-损伤耦合分析计算系统(F-RFPA^2D)．在该系统中，单元的力学、水力学性质根据统计分布而变化，以体现材料的随机不均质性，材料在开裂破坏过程中流体压力传递通过单元渗流-损伤耦合迭代来实现．该系统能够对岩石试件在孔隙水压力和双轴荷载作用下裂纹的萌生、扩展过程中渗透率演化规律及其渗流-应力耦合机制进行模拟分析．最后，给出两个算例：算例1模拟载荷作用下岩石应力应变-渗透率全过程．模拟结果表明，非均匀性对岩石的应力峰值强度、峰值前后其渗透性演化规律及其破裂机制影响十分明显，模拟结果和实验结果较为一致；算例2模拟孔隙水压力作用下岩石拉伸断裂过程，通过和物理实验对比验证，验证了模型计算结果的可靠性。     

岩石边坡裂隙渗流的流形元模拟

为了能够更好地考虑裂隙渗流在岩石边坡破坏中的作用,在考虑裂纹扩展的数值流形方法程序基础上,从最小势能原理出发,阐述了渗流与变形的耦合作用机理并推导了相应的耦合方程,在程序中实现了对渗流与断裂的模拟。最后利用该程序对含初始裂隙的岩石边坡在渗流作用下的破坏过程进行了模拟。模拟结果很好地再现了边坡在渗流作用下的实际破坏过程。同时通过对裂纹扩展过程中应力强度因子的计算发现,第一应力强度因子始终大于第二应力强度因子,且随着计算时间增加,二者的差值逐渐增大,裂纹扩展以Ⅰ型为主。     

平面稳态渗流的随机有限元分析程序

岩体渗流模型中存在着许多不确定的因素，将渗透参数的随机性引入渗流有限元分析计算中会使分析结果更客观、更明确。文章首先基于小参数摄动随机有限元方法，推导出相应的渗流场随机分析中响应量(水头与水力梯度)的随机响应公式，进一步得到流速以及流量的随机响应公式。由于渗流问题与平面弹性问题模型具有较大的相似性，文章在平面弹性问题的...     

岩体系统渗流场与应力场耦合的广义双重介质模型的应用研究

本文通过实验资料分析, 提出了裂隙岩体应力与渗流之间的分形几何关系式。以研究区岩体系统结构为基础, 提出了岩体系统应力场与渗流场耦合的广义双重介质模型, 并将该模型应用于某水利水电工程的渗流与控制中, 对坝区岩体进行了应力场和渗流场的耦合分析, 取得了比较满意的结果。     

岩石破裂过程THMD耦合数值模型研究

从岩石的细观非均匀性特点出发,应用损伤力学、热力学和渗流力学理论,建立了岩体热(温度)-水(渗流)-岩(应力)-损伤耦合数值模型(THMD model),把岩石(体)THM耦合问题的研究从应力状态分析深入到损伤、破坏过程分析之中。探讨了THM耦合作用下岩石材料的细观结构损伤及其诱发的材料力学性能演化机制,并运用所提出方法计算温度-渗流-应力耦合作用下井筒近场围岩的稳定性,模拟得到的岩体破坏过程、应力分布、AE特性及渗流特性变化与现场标定结果有着一致的规律性,初步证明了该数值模型的合理性和有效性。THMD模型以简单的数值模型表征了岩石(体)中热、水、岩及损伤之间复杂的作用关系,为从细观损伤演化揭示宏观岩体温度-渗流-应力耦合破坏机制提供了一种新的数值分析方法。     

可变形多孔介质中的一维非定常耦合渗流

在Ｂｉｏｔ理论的基础上,考虑到可变形多孔介质的渗透系数依赖于孔隙变形的特点,建立了耦合渗流问题的基本方程;用初始层校正法求出了一维非定常耦合渗流问题的摄动解;实例计算表明,耦合分析与非耦合分析之间的判别较大,因此耦合效应不能忽略。     

考虑膨胀力的非饱和介质热-水-应力耦合二维有限元分析

从建立应力平衡方程、水连续性方程、能量守恒方程和弹塑性矩阵入手,使用Galerk in方法,将各控制方程分别在空间域和时间域进行离散,开发出了一个可考虑膨胀力的用于分析非饱和介质中热-水-应力耦合弹塑性问题的二维有限元程序.通过对一个假定的核废料地下处置库的热-水-应力耦合问题的数值计算,比较了无、有膨胀力时的情况,在定性上验证了该程序的正确性.     

岩体中饱和渗流应力耦合模型研究进展

在大坝渗流与控制、各种水工隧洞和交通隧道以及地下各种洞室开挖、采矿工程与油藏开采中的水力劈裂、 岩坡和坝基的稳定性研究等领域中经常会遇到岩体的饱和渗流应力耦合问题.本文对这些领域中的渗流应力耦合问题分为6种力学模型: 等效连续介质模型、 裂隙网络模型、 双重介质模型、 断裂力学模型、 连续损伤力学模型和统计模型.前3种主要从经典的黏弹塑性本构着手考虑, 着重对渗流场的处理, 后3种则是从考虑岩体在耦合作用下发生损伤破裂行为方面切入, 着重于岩体结构内部发生质的改变带来的更加复杂的耦合效应.对这些模型进行了详细的介绍, 指出了这些模型的优缺点与适用范围, 给出了一些有代表性的研究成果, 并对该研究领域中未来的研究方向进行了展望.     

青藏铁路冻土路基温度场随机有限元分析

从一阶Taylor展开式入手,结合冻土温度场计算有限元法,推导了冻土路基温度场随机有限元方程及有关计算公式,并对青藏铁路北麓河试验段冻土路基温度场进行了随机有限元分析,结果表明温度标准差与温度场的分布趋势大致相同,在垂直方向上,最大值都出现在路基上部,向下逐渐减小,但两者分布范围略有不同,路基内部标准差的分布范围明显较大,在一定深度处达到均化。     

SOLUTION AND ANALYSIS OF CIRCULAR TUNNEL FOR THE STRAIN-SOFTENING ROCK MASSES CONSIDERING THE AXIAL IN SITU STRESS AND SEEPAGE FORCE

为了研究轴向应力和渗透力共同作用下软化围岩的应力与位移的变化及分布规律. 基于摩尔-库伦屈服准则及应力-应变软化模型并考虑轴向应力和渗透力的共同作用,将整个塑性区分为有限个同心圆环,以弹塑性交界面处的应力、应变为初始值,并采用微小径向应力增量逐步求出各个圆环上的应力应变及塑性区半径,据此重构了考虑渗透力和轴向力共同作用下软化围岩应力应变特性的逐步求解方法. 利用该方法,推导出软化围岩应力应变的解. 计算结果表明：在考虑轴向应力作用下,塑性区半径和隧道围岩位移都随着渗透力的增加而有所减小;当轴向应力为最小主应力时,渗透力的影响更为显著. 这说明渗透力的存在对于隧道围岩的应力应变分布以及塑性半径和围岩的位移有不可忽略的影响.     

渗流对岩体高边坡位移的影响

根据渗流场和位移场耦合机理 ,结合某大型船闸工程 ,采用有限元法进行计算分析和对比 ,结果表明船闸在开挖过程中 ,渗流作用力对开挖边坡的稳定性影响很大。并指出对高边坡进行稳定分析时 ,正确考虑渗流作用力的影响十分重要。     

黄河龙口水利枢纽库坝区渗漏三维有限元分析

应用三维等参数有限单元法,对黄河龙口水库库坝区周围一定范围的初始渗流场、万家寨蓄水后的渗流场及龙口工程蓄水后的渗流场进行了数值模拟,并计算了相应的渗流量。     

地下水地表水联合模拟的解析单元法理论及应用

本文介绍了解析单元法的基本理论及适用条件 ,并运用该方法 (GFL OW1软件 )计算了某水库蓄水过程中库首结晶岩区风化裂隙水渗流场变化 ,取得良好效果.     

THE STRESS FIELD AT MODE I CRACK TIP ON STEEL PIPE UNDER AXIAL TENSION

采用弹性断裂力学Westergaard 的方法,分析在轴向拉伸载荷作用下钢管的Ⅰ型裂纹尖端处的应力场. 基于Ⅰ型裂纹钢管应力场的特征,设定尖端处的应力艾雷函数,给出其应力场模型边界条件和应力场的解析函数,并利用裂纹尖端处的切平面研究裂纹尖端局部应力场,建立了钢管裂纹尖端应力场模型. 通过钢管与平板Ⅰ型裂纹应力场模型的对比,结果表明二者明显不同,钢管裂纹尖端处应力峰值影响范围仅与裂纹长度、拉应力相关.     

采空区流场非达西渗流一种新的迭代算法

对Bachmat非线性渗流方程求解问题，结合采空区问题的力学特点，提出基于变渗透率达西 (Darcy)渗流求解非线性渗流的迭代算法. 建立了冒落采空区非线性漏风流态的有限元数 值模型，通过对复杂边界采空区的漏风渗流的计算，得到与实际流态更接近的流动规律(风 压等值线和流函数线)和采空区漏风强度分布(速度场). 结果表明，迭代算法是振荡性收 敛的，收敛速度快；与Darcy渗流相比，采空区非线性渗流速度场趋向于平缓，在工作面附 近30\,m范围内最大风速降低到原来的0.61倍，其余大部分区域速度略有增加. 迭代方法满 足工程要求.     

边坡稳定分析的虚功率法

基于有限元计算所得应力场的改进极限平衡法, 对渗流与有效应力耦合作用强烈的边坡或地基的稳定性分析具有优势. 本文提出了边坡稳定分析的虚功率法, 即基于极限分析的上限定理, 利用机动许可的组合刚体滑动机构和有限元应力场, 用滑动机构的速度间断面上的抗滑力功率与滑动力功率的比值计算安全系数. 通过分步优化方法, 获得边坡给定滑动机制的稳定安全系数. 对2个典型折线滑动机构的边坡案例进行了分析, 比较了采用静力平衡应力场和静力许可应力场对安全系数的影响. 指出基于土体线弹性本构模型所得的有效应力场计算的稳定安全系数, 也是边坡稳定安全性的一个不错的度量. 算例中本文计算所得的边坡稳定安全系数, 与文献中推荐的答案很接近, 其滑动机制与有限元强度折减法分析所得的滑动机制基本一致, 安全系数也接近, 表明本文提出的方法是合理的边坡稳定分析新方法, 为边坡和地基的稳定性分析提供了新的选择.      

ANALYSIS OF JANUS PARTICLES IN ELECTRIC FIELD

对外加电场作用下无穷大电解质溶液中二维Janus 导电颗粒的受力特性进行研究。在量级分析的基础上,对电解质溶液电渗流动的基本方程组进行解耦,通过引入有限域确保解的存在性,进而得到圆形颗粒周围电场和流场的解析解与半解析解。基于半解析解给出颗粒所受力矩与外加电场方向角之间的变化曲线,从而得到二维Janus 圆形颗粒的平衡姿态及其稳定性特性。在此基础上,通过数值方法给出了电解质溶液中二维Janus 胶囊的受力特性,平衡姿态及其稳定性特性。     

FATIGUE PROPERTY BASED ON NOTCHED SPECIMEN TESTS

缺口件疲劳问题的研究日益引起各国学者的重视. 局部应力-应变法以其简单性在工程中得到了广泛应用, 该方法通常会得到偏于安全的结果. 引入疲劳缺口因子代替弹性应力集中因子针对缺口疲劳进行研究, 仍未能从本质上改善预测结果的准确性. 考虑到"热点应力" 附近的相对应力梯度, 提出了应力梯度法研究缺口件疲劳问题, 这一概念亦被用于应力场强度方法中, 如何准确确定损伤区域是应力场强度方法需要解决的问题. 临界距离理论可将Neuber 律、Peterson 方法及应力场强度方法进行有效统一, 同时有限元方法的发展进一步支持了该理论. 目前, 该方法在高周疲劳研究中取得了较好的效果, 但对低周疲劳寿命预测的有效性仍需进一步的验证.