运用实验结果及断裂理论分析研究1996年丽江7级地震孕震过程

通过追踪1996年丽江7级地震的孕震过程,发现其震前30余年前兆性中强地震的发生顺序和迁移特征与前人实验结果相吻合。从断裂理论分析,中强震前兆性地震序列及主震发生的先后顺序由孕震区内断裂构造及断裂端部应力核区内介质的E·γ综合决定。这种运用实验结果和相关理论来分析和研究某些天然地震孕震过程的方式,有可能为地震研究和地震预报尝试一条新途径。     

Adaptive integration of constitutive rate equations

In finite element calculations the constitutive model plays a key role. The evaluation of the stress response of the constitutive relation for a given strain increment, which is a time integration in the case of models of the rate type, is a typical sub task in such calculations. Adaptive behaviour of the time integration is essential to assure numerical stability and to control the accuracy of the solution. An adaptive second order semi-implicit method is developed in this paper. Its numerical behaviour is compared with an adaptive second order explicit scheme. The two proposed methods control the local error and guarantee numerical stability of the time integration. We include several numerical geotechnical element tests using hypoplasticity with intergranular strain. The element tests simulate the behaviour of a finite element method based on the displacement formulation.     

基于离散元的冻结砂土三轴力学特性研究

采用离散元软件中的平行黏结模型模拟冻土中冰-土颗粒间相互作用。通过调整模型细观参数,模拟了特定围压不同试验温度和应变速率条件下冻结砂土的三轴试验规律。进一步对比分析计算和试验结果表明,模型的细观参数中的刚度系数和黏结强度对材料的宏观力学行为影响较大,设定的黏结强度参数与温度成反比而和应变速率成正比。在宏观上,峰值强度随温度降低而增大,随应变速率的增大而增大。进一步分析离散元中颗粒的运动产生的剪切带倾角表明：剪切带倾角与温度成反比,而应变速率对其的影响较小;根据已有数据结果,剪切带与内摩擦角间的统计关系符合经典土力学中的摩尔-库伦解。     

考虑填土强度的加筋土挡墙动位移计算

加筋土挡墙在地震荷载作用下的位移大小对结构的抗震性能影响显著。为了计算地震荷载作用下加筋土挡墙的位移,Newmark滑块法通常被用于设计中。由于传统的Newmark滑块法在计算中忽略了填土强度的变化,因而采用单一峰值或残余强度的计算将可能导致计算得到的位移小于或大于实际位移值。在二楔块破坏模式的假定条件下,根据楔块的力学平衡条件,建立了加筋土挡墙滑动安全系数计算式,同时通过引入位移阈值考虑了填土的应变软化特点。通过将计算结果与模型试验结果对比,得到以下结论：相较于单楔块法,二楔块法更能真实地反映出墙体的实际破坏模式,且计算得到的屈服加速度系数更接近试验值;相较于采用单一峰值或残余强度计算的位移,考虑填土应变软化的计算解更接近于模型测试值。     

饱和软土大变形非线性自重固结模型

基于分段线性方法,建立了饱和软土一维自重固结模型（简称SWC模型）。该模型能考虑自重固结过程中土体的大变形效应和材料参数的非线性变化。将该模型的计算结果与相关解析解、现场试验及室内试验结果进行了对比验证,证明了SWC模型能准确计算出大变形和非线性条件下饱和软土的自重固结过程,包括沉降量、平均固结度、孔隙比分布和超孔隙水压力分布等参数随时间的变化过程。随后,以现场试验为基础,采用SWC模型对饱和软土自重固结的4个主要影响因素（即土体初始高度、边界排水条件、初始孔隙比和土粒相对密度）进行了参数分析。结果表明,上述4个参数对软土自重固结过程均具有重要影响：土体初始高度越高,则自重固结沉降量和最终平均应变值越大;边界排水条件对土体自重固结的速度有重要影响,但不影响自重固结的最终沉降量;初始孔隙比越大,则自重固结沉降量越大,其完成自重固结所需时间越短;土粒相对密度越大,则土体的最终沉降量越大,完成自重固结所需时间越短。     

一种适用于孔隙体积应变的有效应力方程

土力学奠基石Terzaghi有效应力原理被广泛应用于油藏孔隙和渗透率应力敏感研究中,然而其对于岩石孔隙体积应变的适用性存在争议。对颗粒不可压缩和颗粒可压缩的多孔介质分别进行了受力分析,推导了总体积、颗粒骨架、孔隙体积的有效应力表达式,与Biot、Skepmton有效应力方程对比,建立了适用于孔隙体积应变的新型有效应力方程,并进行了试验论证和应用举例。研究表明:在颗粒不可压缩多孔介质中,有效应力为超出平衡孔隙流压之外的颗粒间宏观等效应力;在颗粒可压缩变形多孔介质中,有效应力为其相同应变下的等效应力,有3种有效应力分别适用于总体积应变、颗粒体积应变、孔隙体积应变;新提出的孔隙体积有效应力方程与孔隙度、岩石总体积压缩系数、颗粒压缩系数、总应力和流压相关,4个理论计算式计算结果在3种多孔介质试验测试结果中的偏差均在5%以内;孔隙体积有效应力系数解决了如何定量增总应力来等效模拟储层降流压生产过程这一关键问题,3个压缩系数关系式理论计算准确方便。     

有明黏土各向异性固结特性的试验研究

由于大多数黏土矿物是片状的,在沉积和随后的固结过程中黏土颗粒和组构单元趋向水平定向排列,使天然沉积黏土呈现各向异性的固结特性。通过垂直切和水平切试样的恒应变速率固结试验研究了日本有明黏土的各向异性固结特性。试验结果表明,水平切试样的先期固结压力pch与垂直切试样的先期固结压力pcv的比值在0.5～1.0之间,主要是因为天然沉积黏土各向异性的屈服轨迹以及垂直切和水平切试样所经历的应力路径不同;在正常固结阶段,水平切试样的固结系数ch与垂直切试样的固结系数cv之比约为1.43,水平切试样的渗透系数kh与垂直切试样的渗透系数kv之比约为1.40,水平切试样的体积压缩系数mh与垂直切试样的体积压缩系数mv基本相等。因此,有明黏土固结系数的各向异性主要是因为渗透系数的各向异性,而渗透系数的各向异性本质上是因为其微观结构的各向异性。     

水−岩作用下砂岩力学特性及微细观结构损伤演化

在库水位反复升降作用下,库岸边坡消落带部分岩体处于浸泡―风干循环状态,其力学特性的劣化将直接影响库岸边坡的变形稳定。以三峡库区库岸边坡消落带砂岩为研究对象,设计并进行了考虑浸泡―风干循环和水压力升降变化的水?岩作用试验。结合SEM电镜扫描,分析了水?岩作用对砂岩微细观结构的影响,采用PFC2D离散元软件对水?岩作用下砂岩的劣化过程进行了模拟分析。研究结果发现：（1）水?岩作用下砂岩的力学性能劣化效应明显,弹性模量和单轴抗压强度总体呈现先陡后缓的劣化趋势,其中前6次的浸泡―风干循环作用导致的劣化效应尤为明显。（2）水?岩作用下砂岩的微细观结构变化特征明显,由较致密的孔隙式胶结结构逐渐向不规则的蜂窝状结构发展,宏观上就导致了砂岩变形和强度特性的劣化。（3）由PFC2D离散元模拟分析发现：水?岩作用导致岩石内部颗粒强度和颗粒间黏结强度逐渐劣化,非均匀性增强,荷载作用下的应变局部化特征逐渐明显,整体承载能力逐渐降低,破坏模式由典型的张拉破坏逐渐转变为剪切破坏。     

3D-DIC系统在岩石力学试验中的应用

介绍了全场3D-DIC系统,在试样周围均匀布置三组3D采集元,利用散斑图像的采集,尽可能在较大范围内捕捉试样表面的变形并进行分析。为了验证3D-DIC系统测量结果的可靠性,以江持安山岩为试验材料,进行了单轴压缩条件下的恒定应变速率和交替应变速率试验。通过3D-DIC系统的测量结果与应变片测量结果的对比发现,两者吻合度较好,表明3D-DIC系统能够满足岩石力学试验非接触式、可视化变形测量的需求。另外以田下凝灰岩为试验材料,通过图像分析得到了加载过程中岩石表面的全场变形信息,发现岩石变形的不均匀性和应变局部化现象,表明利用3D-DIC系统对岩石破裂过程变形测量比传统变形测量方法更具有优越性,为岩石力学试验研究提供了一种新的、便捷的测试方法。     

考虑超固结效应的不同温度路径下饱和粉质黏土的热固结

通过室内试验,研究了超固结比OCR对不同温度路径下中空圆柱饱和粉质黏土热固结特性的影响。研究表明,同一升温过程、不同超固结比试样所产生的最大孔压基本相同,但在相同的降温过程中,随着OCR的增大,产生的负孔压的绝对值有增大的趋势;升温或降温过程产生的体应变随着OCR的增大而减小,呈指数形式变化,并且温度幅度越大,这种变化幅度也越大;经过升温-降温过程后的最终体应变也随着OCR的增大而减小;升温或降温到同一温度值,其过程设置的温度等级越多,固结体应变也越大,并且随着OCR的增大,不同等级之间的体应变差值减小,而OCR对不同等级升温-降温后的最终体应变间的差值影响较小。