傅里叶级数法内插及外推的绘图计算程序

本文根据某些地质变量的波状变化,应用傅里叶级数拟合出地质变量的波状变化曲线,并寻找出一个最佳拟合函数,以此来进行地质变量的外推及内插.绘图计算程序采用人机对话的形式,用BASIC语言编写.绘图语言指令适用于CE-515P四色宽行(80列)绘图打印机.、绘图的比例尺,可根据实际需要,选择适当的参数(S_2和S_3的大小)来自动调节观测曲线和拟合函数曲线的比例.程序较优化.经过实际工作的多次应用及检验,应用效果较理想.     

应力波在非线性结构面介质中的传播规律

用切线刚度和法线刚度描述结构面特性,研究结构面初始刚度、频率、法向闭合量与其最大允许闭合量的比值对透射系数的影响。采用基于连续介质的块体离散元程序（CDEM）模拟结构面发生非线性变形条件下块体响应,研究结果表明,应力波在岩体中传播是一个传播和块体响应的过程,结构面的存在影响了应力波传播和响应,存在一个特征频率能够有效反映结构面刚度变化,给出了近似计算该特征频率的表达式,对岩体结构探测具有一定的指导意义。     

位移不连续模型在黏弹性节理刚度计算的应用

考虑节理质量对应力波传播的影响,运用波的位移位函数推导了谐波在厚黏弹性节理的透、反射系数计算公式,采用波形相关系数描述子波穿过黏弹性节理的波形变化,讨论具有一定厚度的黏弹性节理简化位移不连续模型的适用条件。设厚黏弹性节理模型和位移不连续模型的透射波波形相关系数为0.9时对应的节理厚度为临界厚度,岩体与节理的阻抗比对临界厚度的影响很小;临界厚度随子波中心频率增大呈负指数减小,入射角越大,临界厚度随中心频率减小得越慢。试验数据分析表明：当节理厚度为0.03 m时,采用位移不连续模型和厚黏弹性模型计算得到的节理力学参数非常接近,随节理厚度和子波中心频率增加,运用位移不连续模型的计算结果偏差越大,试验结果与理论分析是一致的。     

岩体随机结构面三维网络数值模型有效性检验

岩体随机不连续面三维网络模拟技术成熟与否的一个重要标志就是最终的网络模型是否能够较真实地再现岩体的实际情况 ,为了使三维网络模型的结果能够真正地应用于实际工程 ,就必须对模拟的结果进行必要的检验 ,检验的办法就是将模拟出的模型剖面图形与原型测量窗口的图形进行图形对比 ,或者将原型的观测数据与模型的测试数据进行数据对比 ,当这种检验达到精度时 ,模拟出的模型为有效模型 ,可以投入工程实际应用     

地基不均匀沉降对上部结构影响的三维弹性支承分析法

地基不均匀沉降是导致建筑工程事故的主要原因之一。为了研究地基不均匀沉降对上部结构内力、变形、周期等影响,提出了将地基(或基础)模拟为三维弹性支座的有限元分析方法,将上部结构和地基基础进行整体分析。以深圳市少年宫“水晶石”大厅结构为例,运用有限元软件ANSYS,将三维弹性支承分析法与常用的固定支座法和竖向弹性支承分析法进行了对比分析,验证了三维弹性支承分析法的合理性,从而得出了一些有益的结论,可供设计人员参考。     

黏弹性节理岩体的位移不连续模型研究

选用三参数标准线性固体作为岩石本构,提出了一种考虑岩体黏弹性的位移不连续模型;根据一维黏弹性波的特征线法,推导了节理处质点速度、应力和应变递推公式。首先,基于分离式霍普金森压杆（SHPB）对砂层进行试验,得到其应力-应变关系,并换算出砂层节理的法向刚度;接着,通过一维强间断黏弹性波的波速公式、高频波衰减系数以及任一频率下的衰减系数,确定出数值算法中的三参数。最后,基于自制的摆锤装置,探讨了一维应力波在节理岩体中的传播规律,试验中以两根长1 000 mm、直径为68.50 m的岩杆作为入射和透射杆,以3 mm砂层模拟节理。试验和数值结果吻合度良好,进一步验证了该方法的可靠性。     

岩体结构二维可视化模型及其应用研究

岩体结构图形可视化是岩体结构研究的一个有效手段 ,它可以为地质工程师和工程设计人员提供直观的图形资料 ,以便他们做出合理的地质判断和工程设计。我们基于随机不连续面圆盘理论提出了岩体结构二维可视化模型 ,解决了岩体结构在平面图和剖面图中的定位和绘制方法 ,并将这一方法应用于西南某大型水电站坝肩抗力体岩体结构研究中 ,取得了较好的效果 ,实践证明该方法是可行的。     

岩体随机不连续面三维网络模型的检验

岩体随机不连续面三维网络模拟技术成熟与否的一个重要标志就是最终的网络模型是否能够较真实地再现岩体的实际情况,为了使三维网络模型的结果能够真正地应用于实际工程,就必须对模拟的结果进行必要的检验,检验的办法就是将模拟出的模型剖面图形与原型测量窗口的图形进行图形对比,或者将原型的观测数据与模型的测试数据进行数据对比,当这种检验达到精度时,模拟出的模型为有效模型,可以投入工程实际应用。     

岩体随机不连续面三维网络模型的检验

岩体随机不连续面三维网络模拟技术成熟与否的一个重要标志就是最终的网络模型是否能够较真实地再现岩体的实际情况，为了使三维网络模型的结果能够真正地应用于实际工程，就必须对模拟的结果进行必要的检验，检验的办法就是将模拟出的模型剖面图形与原型测量窗口的图形进行图形对比，或者将原型的观测数据与模型的测试数据进行数据对比，当这种检验达到精度时，模拟出的模型为有效模型，可以投入工程实际应用。     

基于非连续单元法的露天矿轮斗剥离台阶片帮灾害分析

针对某露天煤矿轮斗剥离台阶片帮情况,对剥离台阶岩土体进行了物理力学试验,基于有限差分方法,从边坡内部应力应变特征及塑性区扩展情况开展了单台阶片落体失稳机理分析,基于改进的结构面软化模型,采用连续-非连续单元方法(CDEM)对片落体滑塌及堆积的全过程进行了数值模拟分析,对片落体的运动轨迹及堆积情况进行了模拟分析,结合现场地质环境特征,分析了单台阶片帮的触发因素,结论表明,改进的结构面软化模型可作为地质体的渐进破坏的判据,真实还原片落体垮塌全过程;采用连续-非连续单元方法可模拟材料由连续体到非连续体的渐进破坏过程,实现大变形滑坡的灾变全过程分析,研究成果可为类似工况的边坡的评价及治理提供参考。     

共传播纵波研究砂岩非线性弹性特性

静态和动态的非线性弹性探测实验已经清晰地揭示了岩石具有较强的非线性弹性特性,其中动态探测方法运用共振在岩石内部加载稳定的应变取得了较好的效果,但是野外探测具有一定的难度.本研究利用传播的纵波在岩石内部加载应变,发展了共同传播的纵波研究岩石非线性弹性特性的方法.该方法通过相对低频高振幅的纵波在岩石内的传播在岩石内部加载一...     

边坡非线性位移的分段最小二乘-时间序列分析

利用最小二乘法对位移时间序列进行分段拟合以提取时间序列的趋势项，从而将非稳定时间序列稳定化，对去除趋势项的稳定时间序列即残差进行常规AR（p）时序分析。结合滚动预测方法，建立边坡失稳预测的叠合模型，以长江三峡某边坡为例，对该模型进行的检验表明：新模型的预测精度较高、实时可靠。     

高边坡变形非线性时变统计模型研究

基于高边坡长期变形监测资料,深入分析了边坡变形的影响因素,建立边坡非线形时变统计监控模型,并将模型应用于锦屏一级水电站边坡工程。通过典型测点长期监测资料的分析,确定时效和降雨是影响边坡变形的2个主要因素,并给出了各自的表述形式,在Matlab平台上,采用非线性回归方法得到边坡非线性时变统计模型的参数。在此基础上,根据边坡变形拟合残差时序,采用ARMA模型方法进行拟合和预测,对边坡非线性时变模型的拟合和预测进行修正。结果表明,模型能较好地描述高边坡的变形特征,并且具有很好的拟合和预测精度。     

黏性土非线性弹性K-G模型的一种改进方法

为解决《土工试验规程》建议的K-G模型未考虑土体剪缩性问题并保持其简单性,在其基础上增加了一个剪缩模量,建议了一个3模量7参数K-G-D模型,给出了用 和 与 和 表示的完整的增量应力应变关系式;通过正常固结黏土的各向等压固结试验和等平均正应力的三轴压缩排水试验标定了模型参数;进行了黏土的常规三轴压缩排水试验,利用模型参数预测了相应路径下土样的体积与轴向变形,经与实测值对比,预测结果较好。该模型可考虑土体压硬性、剪缩性、应变强化特性,能较完善地描述欠固结和正常固结土及松砂和中密砂在等 和增 路径下的变形特性;模型并未在规程模型的基础上增加试验工作,且概念清晰、形式简单、便于实际应用。     

应力波斜入射对岩体节理端部的影响

岩体节理在受到应力波扰动时端部受力会发生不同程度的连续性变化,为具体分析这种动态变化与应力波入射角度之间的关系,引用岩石非线性法向本构关系与线性切向本构关系的组合模型以及相应的在P波斜入射节理的应力波传播方程,结合岩石断裂力学中Ⅰ、Ⅱ型裂纹尖端应力和位移场的相关理论,得出组合形式下的节理端部应力场和位移场随质点振速变化的计算公式。通过不同入射角的模拟脉冲信号作用对端部应力位移的变化分析：应力波斜入射节理导致节理端部上、下两侧应力与位移场非对称分布,随着入射波质点振速的增减变化,应力集中位置会出现变化;计算并整理模型中节理端部上、下两侧0.005 m位置的数据,节理法向刚度由入射波质点振速带来的变化直接影响到应力波的透射与反射,进而会导致端部的应力与位移出现滞后效应;节理端部横向位移值与入射角度并非单调变化,而竖向位移会随着入射角的增大呈现下降趋势。     

基于球面波势函数基本解方法的弹性波三维散射与动应力求解

针对弹性波三维散射和动应力集中问题,提出一种基于球面波势函数的基本解方法(SWP-MFS)。方法基于单层位势理论,采用膨胀波和矢量剪切波球面波势函数构造散射波场,根据边界条件建立边界积分方程并配点求解。精度检验表明,该方法具有良好的数值精度及稳定性。以无限空间中三维夹杂体及空洞对平面P、SV波的散射为例,进行方法展示,并揭示了三维夹杂体周围弹性波散射的一些重要规律。结果表明:三维夹杂体随其内部介质刚度降低,位移谱曲线震荡越加强烈。三维球形空洞在P波和SV波水平入射下应力集中规律不同,前者在顶部和底部及附近更明显,后者在纵截面两45°角线附近更明显。与以往的集中力源函数相比,新的波场构造基本解更为简洁易用,为三维弹性波动分析提供了一种新型无网格数值方法。     

RVSP弹性波数值模拟及传播特征分析

基于纵波、横波解耦的弹性波高阶有限差分方程和PML吸收边界条件,实现了RVSP观测统中弹性波数值模拟。采用完全弹性波波动方程进行数值模拟,可以得到纵波和横波的混合波场,且波场丰富,符合实际地震波的传播规律。该正演模拟方法纵波、横波自然解耦,产生全波场、纯纵波和纯横波模拟记录。通过对层状介质模型、凹陷模型,以及实际复杂介质模型的RVSP弹性波进行数值模拟,得到各模型不同分量下全波、纯纵波和纯横波的波场,并对弹性波传播特征进行了分析,为下一步的RVSP地震资料处理和解释工作奠定了理论基础。     

基于弹性波波速的降雨型滑坡预警系统

降雨是诱发滑坡最主要的因素。为减少滑坡灾害造成的人员伤亡和经济损失,滑坡早期预警系统成为了最佳选择之一。根据弹性波传播的基本原理和基于降雨型滑坡变形破坏的特点,提出利用弹性波波速来反映边坡表面土体含水率和位移的变化。开发研制了一套三轴渗流-弹性波测试三轴仪和相关系统。该装置能让水从底部渗入土体,模拟降雨入渗土体的过程,同时能测试弹性波波速。试验过程中同步测试含水率、变形和弹性波波速的变化。还进行了降雨滑坡模型试验。利用三轴弯曲元注水试验和降雨滑坡模型试验,深入分析和研究降雨引起的土体滑坡过程与弹性波波速演化规律,揭示波速与含水率和变形的耦合关系。研究结果表明,弹性波波速随着含水率的增大而缓慢减小,随着变形的增大而急剧减小,临近失稳时,波速骤然减小。根据试验结果对含水率和变形导致弹性波波速减小可能的机制进行了解释,提出弹性波在波速骤然减小时发出滑坡预警。研究结果为滑坡防灾减灾和预测预报提供新的方法和可靠的依据。     

高阶交错网格有限差分弹性波场模拟的精度分析

交错网格波场数值模拟是目前地震正演中广泛使用的方法,为对比分析不同阶数的差分格式下产生的计算效率和精度差异,重新推导了弹性波方程的4种时间4阶、空间2N阶的差分公式及系数,并计算了他们的稳定性条件。利用这4种差分格式进行弹性波场数值模拟,对比分析了波场快照、合成地震记录及CPU时间。结果表明：时间4阶、空间6+6阶精度的交错网格有限差分方法在进行地震波场数值模拟时具有较高的计算精度和计算效率。     

Scattering of elastic wave by a partially sealed cylindrical shell embedded in poroelastic medium

Scattering of an elastic wave by a cylindrical shell embedded in poroelastic medium is investigated theoretically with the assumption that the shell material is also a porous elastic medium. The porous medium is modellized via Biot's theory and the scattering by cylindrical shell is expressed by the definition of scattering matrix. The normal mode expansion technique is employed for analysing the scattering field, and the asymptotic solutions of displacements, stresses and pore pressure are derived. Two limiting cases‐scattering by a poroelastic cylinder in Biot medium and a elastic cylindrical shell in elastic medium are obtained from the general solutions. The dispersion curves of displacement amplitude at the interface of shell and medium is compared with the case of elastic shell. The scattering amplitude associated with the fast, slow and transverse waves are identified by numerical simulation. Furthermore, the influence of the poroelastic property of shell material on scattering amplitude is analysed. Copyright © 2005 John Wiley & Sons, Ltd.