岩石力学与工程学报杂志 北大期刊 CSCD期刊 统计源期刊

岩石力学与工程学报 2006年第11期杂志 文档列表

岩石力学与工程学报杂志动态

第六届《岩石力学与工程学报》编委扩大会议纪要2246-2246

摘要：2006年9月13日晚，在沈阳市金城宾馆召开了第六届《岩石力学与工程学报》（以下简称《学报》）编委扩大会议，共有60余位编委或其代表参加了本次会议。会议由《学报》主编冯夏庭研究员主持，首先由《学报》常务副主编余诗刚研究员汇报了《学报》近年来的发展态势和今后的主要工作：中国岩石力学与工程学会理事长钱七虎院士最后代表学会作了总结讲话，认为《学报》经两届编辑部班子的努力，得到了发展，影响力加大，学术权威性增强。与会代表就如何办好《学报》展开了热烈的讨论，并达成以下共识：

三维块体接触判断方法的分析与改进2247-2257

摘要：非连续变形分析方法和离散单元法都是可以计算非连续介质大位移问题的数值方法，目前在二维问题的计算中应用较多，但在工程实际问题的计算中还未能得到广泛应用，因为实际的工程问题如岩体结构都是三维问题。而在计算三维问题时存在一个困难，就是难以快速、准确地判断块体之间的接触类型。对已有的3种接触判断方法尤其是其中已用于实际计算的公共面法、侵入边法进行分析，指出各种方法在进行接触判断时存在的问题。同时在综合已有方法优点的基础上，给出建立正确接触判断方法的原则并提出切割体法。2个算例的计算结果表明，该方法的确有效地解决了已有方法存在的问题，能够正确判断三维块体之间的接触关系，从而模拟块体之间的错动滑移过程。

云驾岭煤矿无烟煤的动态本构模型研究2258-2263

摘要：大直径霍普金生压杆试验表明：云驾岭煤矿无烟煤的动态本构曲线具有显著的塑性屈服特征，初始弹性模量Eb、屈服强度σs与极限强度σb都随着应变率的增加而提高，但屈服强度最为显著。根据试验曲线的特征，在对既有岩石动态本构模型进行分析和初步的拟合后，引入一个新的损伤型黏弹性本构模型，以1根线性弹簧和2个不同松弛时间的Maxwell体并联的黏弹性模型表达无烟煤显著的塑性流动，2个Maxwell体分别表示无烟煤的低应变率响应和高应变率响应，以动态应力-应变曲线任-点与原点连线的斜率E（εi）相对于初始模量Eb的降低来度量损伤D。模型的数值拟合曲线与实测动态本构曲线具有较好的一致性。拟合参数表明，无烟煤对低应变率的响应与混凝土相同，对高应变率的敏感性远远高于混凝土。

碎石土渗透特性对滑坡稳定性的影响2264-2271

摘要：结合工程实践，通过现场工程地质调查与勘探和室内外物理力学试验，采用数理统计分析方法和不平衡推力法，对碎石土的一般物理力学特性、渗透特性以及碎石土滑坡中地下水管网状排泄系统对滑坡稳定性和失稳破坏的作用机制进行分析和研究。结果表明，碎块石粒组的含量和以粉粒、黏粒为主的细粒土粒组的含量对碎石土渗透系数影响最大且最为显著，碎石土的渗透系数分别随土中碎块砾石粒组含量和小于0．1mm粒径的细粒土粒组含量的增加而呈自然指数增大和降低；碎石土滑坡中地下水管网状排泄系统对天然边坡的稳定起关键作用；当碎石土滑坡中地下水管网状排泄系统遭受破坏和堵塞时，地下水位将明显上升，从而使潜在滑面的孔隙水压力及下滑力明显增大，导致边坡稳定性系数极大地下降，甚至使边坡失稳解体破坏。

特低渗透油藏非达西渗流模型及其应用2272-2279

摘要：根据矿场试验和室内试验的研究成果，开发低渗透和特低渗透油藏时会产生非达西渗流现象，尤其是原油渗流会存在较为明显的启动压力梯度。结合低渗透油藏注水开发的特点，考虑油、水两相渗流均存在启动压力梯度，以及岩石存在应力敏感，建立三维油、水两相非达西渗流数学模型。利用差分方法得到该模型的数值差分格式，在此基础上编制模拟软件。以大庆肇源油田源121—3区块为研究对象，首先将建立的非达西渗流的数值模拟结果与现场实测的生产动态数据进行对比验证。然后，在五点注采井网的基础上，设计一个“正方形注采单元”，在保持采油井和注水井井底流压的前提条件下，改变注采井距，通过模拟计算得到油井见水时地层流体压力、流体压力梯度、地层绝对渗透率和启动压力梯度的空间分布，最后得到有效驱动注采井距，从而为该特低渗透油藏的注采井网部署提供依据。

基于模糊理论的软岩黏弹性模型识别及参数反分析2280-2286

摘要：对岩石黏弹性模型进行总结，将复杂的黏弹性模型方程转化为统一的线性方程，基于实测数据运用回归分析法和模糊数学方法，通过对线性方程参数的反分析得到统一模型的参数，运用模糊关联度分析法对具体岩石黏弹性模型进行识别，进而根据统一模型参数可反求出黏弹性模型方程的参数。研究表明，对黏弹性本构模型采用线性化处理，用统一的线性方程表示复杂的本构方程简化模型的识别和参数的获取过程。最后利用一矿山软岩流变试验结果对其黏弹性模型进行识别和模型参数进行反分析。计算结果表明，样品Ⅰ-1，Ⅰ-2黏弹性模型服从麦钦特体，样品Ⅱ-1，Ⅱ-2黏弹性模型服从H-K体，通过计算出黏弹性模型参数，得出软岩黏弹性模型的回归方程，所得结果与实际情况比较一致。

复杂泥页岩地层地应力的确定方法研究2287-2291

摘要：掌握地层原始地应力状态是有效控制钻井过程中泥页岩地层井壁失稳的前提。在分析水压致裂法、Kaiser效应法、差应变法和多极子测井方法测量地应力的技术现状的基础上，发现利用单一方法获取复杂泥页岩地层的地应力比较困难，将多种方法组合在一起才能较好地解决问题。如果地层破裂试验层段有钻井岩芯，可通过差应变方法获得主地应力的比值；如果地层破裂试验层段没有钻井岩芯，可通过多极子测井方法得到两个水平主地应力的差值。将两个水平主地应力的相对值代入地层破裂试验的破裂压力计算模型，可分别确定出两个水平主地应力的大小。利用地层破裂试验与差应变试验组合测地应力方法和地层破裂试验与多极子测井组合测得应力的方法，分别较为成功地测试地应力。

基于可见性和后修圆的三维块体域接触算法2292-2297

摘要：在离散元等三维非连续介质数值分析方法中，快速高效的接触检索算法是控制计算效率的关键性因素之一。在直接接触检索算法的基础上引入域分解、外接球判断、可视化检测以及后修圆方法，提出一种新的三维凸多面体接触检索算法。首先阐述该算法的基本原理，分析邻近块体的识别方法以及接触类型和接触方向的确定方法，而后以承受冲击作用的双柱碰撞过程为例表明该接触检索方法计算工作量小，鲁棒性好，易于计算机实现。

锦屏一级水电站左坝肩边坡稳定性研究2298-2304

摘要：从工程区的自然稳定斜坡调查入手，利用前期勘测的横剖面结果，获得工程区的砂板岩边坡和大理岩边坡的自然稳定斜坡的坡形数据；结合中国已建和在建水电边坡的开挖情况，采用自然斜坡比拟法和工程地质类比法综合确定拱间槽边坡的地质建议开挖坡比。在坝头深裂详细调查研究的基础上，分析拱间槽边坡开挖前后的可能失稳破坏模式及失稳块体的边界条件；通过采用三维极限平衡方法对可能失稳块体的计算模型进行分析，得到块体开挖前后的稳定系数；对失稳块体的稳定性进行详细评价，并提出块体加固方案建议，对类似工程提供可以借答的经验。

溪洛渡水电站坝区初始地应力场反演分析研究2305-2312

摘要：初始地应力是岩土工程设计与施工所要考虑的重要因素之一。由实测地应力资料反演初始地应力场的方法很多，且各有优势。根据溪洛渡水电站工程坝区地应力的实测资料，采用有限单元法，结合多元线性回归方法、神经网络方法和遗传算法，分别反演求得整个坝区的初始地应力场。比较发现3种方法反演计算的结果非常接近，且均能模拟实际地应力场的分布规律。进一步的分析结果表明，由于溪洛渡水电站工程坝区在天然情况下岩体的屈服范围很小，坝区初始地应力场的非线性特征不明显，故线性回归分析方法即可满足要求。同时，研究还表明，在结果均能满足要求的前提下，回归分析方法与神经网络方法和遗传算法相比具有更方便快捷、易于掌握、唯一解的优势。

本刊2006年第7期被EI收录论文（31篇，收录率为100％）题录2312-2312

三轴应力状态下不同湿度原状黄土的结构性定量化参数2313-2319

摘要：通过对黄土结构性研究现状的分析，对与原状样干密度相同的扰动饱和土施加荷载来实现破坏结构性的目的，以三轴试验得到的原状黄土以及扰动饱和黄土的应力-应变曲线为基础，用一个定量化结构性参数综合反映土的排列和胶结特征所表现出的综合结构势（结构性强弱），探讨不同围压及含水量下结构性参数在剪切过程中的变化特性。研究结果表明：（1）定义的结构性参数能够反映围压、含水量对原状黄土剪切过程中结构性的影响，具有很好的合理性、灵敏性和稳定性：（2）结构性参数随着含水量的增大而减小，低应变时含水量对结构性的影响很明显，这种影响随着应变增大而逐渐减弱；（3）初始加荷和试样破坏时，结构性参数随着围压的增大而减小，应变达到15％时，结构性参数随着围压的增大呈先增大后减小的趋势。

土体平面应变条件下的主应力关系2320-2326

摘要：平面应变条件下土体中主应力之间的关系，其实质是特定应变路径条件下土的应力-应变关系，为土的一个重要力学特性。在平面应变条件下将土体的变形和强度特性简化为二维问题处理时，其关键就在于如何简单、合理地确定平面应变方向上的主应力大小。基于弹塑性理论及土的试验规律，在平面应变的简单加载条件下，假定主应力之间的关系为双线性函数，并利用一维固结和破坏时的主应力状态，确定双线性函数的系数，得到平面应变方向上主应力计算公式的一个较合理形式，为应力状态和材料性质的函数。建议的主应力关系体现土体在平面应变条件下的变形机制，在b-R坐标内通过3个关键点：初始点（1，1）、一维固结应力状态点（bc，R）、破坏应力状态点（bps, Rps）。通过与试验数据的比较证明所提双线性主应力函数的合理性，它可合理地描述平面应变方向上主应力的变化趋势，并且形式简单，易于工程应用。

格栅加筋土挡墙数值分析的复合材料方法2327-2334

摘要：格栅加筋土挡墙由面板、格栅和填土组成，其结构特性一直是工程界关注的焦点。基于自洽理论和应变相容观点建立筋土复合材料力学模型，用来模拟格栅及其周围一定厚度的土体，依此对一台阶式格栅加筋土挡墙的结构特性进行有限元分析，并与大尺度模型试验测试结果进行比较。研究表明，该方法计算结果和模型试验测试结果在趋势上和数值上都表现出较好的一致性。

基于2．5维有限元方法分析列车荷载产生的地基波动2335-2342

摘要：开发一种2．5维有限单元结合薄层单元的高效数值方法来研究铁路轨道和周围地基在列车运动荷载作用下的振动响应以及应力波在三维地基中的传播问题。其中轨道被简化成铺设在无限半空间上的欧拉梁，具有复杂构造和地质条件的路基近场区域则采用有限单元来模拟，并且利用动力薄层单元构建透射边界来模拟波在远场土体中的传播。通过沿轨道方向的波数变换将三维问题降维成每个节点有3个自由度的平面应变问题，对求解得到的结果进行波数扩展，获得轨道和地基在移动荷载作用下的三维动力学问题的解。考虑单一移动荷载和列车多轮重荷载两种情况下轨道的振动以及周围地基中的响应波场，并用瑞典X2000高速列车运行的现场测试数据对计算结果进行验证，分析列车荷载作用下地基的响应特性和在特定场地上运行高速列车时临界速度的确定。

岩石力学与工程学报杂志短文（研究进展与工程实录）

关于Sarma法改进算法的补充2343-2345

摘要：对作者提出的Sarma法改进算法做进一步的补充，给出保证安全系数稳定收敛的迭代程式，通过算例说明方法的有效性。

单元集成法及其在土压力计算中的应用2346-2351

摘要：基于塑性力学上限定理的土压力计算，在数值上可以通过单元集成法来实现。它采用类似于有限元网格划分的方式离散墙背土体，并设定一个机动许可的滑动机构。在此滑动模式下，可以求得每个单元贡献的外力功率和内能耗散率。所有单元的能量变化率的总和就是墙背土滑裂体的总能量变化率。然后，根据上限定理可以求得与滑动机构相对应的极限外载荷，并通过非线性数学规划法找到其最小值。采用平面和对数螺旋滑裂面的单滑块机制，对典型的二维土压力问题进行分析，其结果说明该方法的有效性。

岩石力学与工程学报杂志动态

本刊2006年第8期被EI收录论文（32篇，收录率为100％）题录2351-2351