横观各向同性饱和两相介质动力反应计算分析

针对横观各向同性饱和两相介质的弹性波动方程组,应用基于显式逐步积分计算格式的时域显式有限元方法对其在输入地震波作用下的动力反应进行计算和分析,并将计算结果与完全各向同性饱和两相介质的计算结果进行对比研究。结果显示横观各向同性饱和两相介质与完全各向同性饱和两相介质的动力反应具有较为显著的差异。同时表明时域显式有限元方法是进行横观各向同性饱和两相介质动力反应计算分析的一种有效方法。     

强度折减有限元在堤防稳定分析中的应用研究

采用强度折减有限元法,研究非均值堤防塑性区的开展特征和失稳破坏过程,并根据堤防应力场分布和临界滑动面的形成来分析堤岸整体失稳破坏的机理.结果表明:把强度折减有限元应用到非均质土层的堤防边坡稳定分析,在理论上、数值模拟实现上都是可行的;有限元静力平衡和位移计算不收敛作为堤防边坡整体失稳的标志,同时考察滑移面上某些特征点的...     

渭河盆地基底伸展与地裂缝成因关系探讨

渭河盆地是我国乃至世界上地裂缝最为发育,灾害最为严重的地区。本文以渭河盆地为原型,采用有限元数值分析方法,分析了基底伸展作用下盆地浅表层岩土介质和多级破裂系统的应力和变形响应特征,从而揭示基底伸展变形与该区域地裂缝之间的成因关系。结果表明,盆地基底的伸展作用可以引起近地表岩土介质的拉张变形,当这种拉张作用与盆地断裂的上盘正断倾滑式伸展拉张叠加时,必然形成或加剧地表土层的张剪性破裂,从而为地裂缝的发育奠定了构造基础。     

不排水不可压缩条件下两相介质的两重网格算法

对于不排水、不可压缩饱和软土地基的固结问题的有限元分析,可以用Biot固结方程来考虑土体颗粒与孔隙水间的相互作用。由于受Babuska-Brezzi稳定条件的限制,用常规的等插值u-p混合有限元法求解将导致孔隙压力出现紊乱的结果。提出了基于位移和压力线性等插值函数的两重网格,但位移独立变量总数大于独立压力变量总数的计算方法,可以满足Babuska-Brezzi稳定条件,使得位移场和压力场单元插值阶数保持一致。通过几个简单算例验证了提出方法的正确性。     

顶板诱导崩落模式选择时变数值分析

顶板诱导崩落是一种新的空区处理技术,其实施效果与诱导崩落施工路径、方法和顺序有关,同时矿岩体作为一种弹塑性体,具有非线性时变力学特征。运用时变力学的基本理论,建立了顶板诱导崩落的时变力学有限元基本方程。针对大厂铜坑矿92号矿体试验采场的地质特征,构建了顶板诱导崩落的时变力学有限元模型,采用多步骤开挖模拟两种不同顶板诱导崩落模式,研究两种不同工序的顶板塑性区发展、东西预裂与崩顶硐室的位移及其安全系数,分析了其对采场的综合效应。结果表明,先预裂后崩顶的顶板诱导崩落模式有利于顶板崩落诱导,并综合考虑诱导崩落效果与作业安全,建议采用先预裂爆破后强制崩顶的微差爆破一次成型工艺。     

多孔介质中含均相/非均相反应传质过程数值模拟

结合化学反应方程式,并应用多孔多相介质溶混污染物输运过程的数值模型,对多孔多相介质中含均相/非均相化学反应传质过程进行了数值模拟。化学反应主要包含均相快速/慢速和非均相快速/慢速等5种化学反应过程,溶质输运行为的控制机制主要考虑对流、扩散及降解、吸附等。基于原有的隐式特征线Galerkin离散化的有限元方法,求解模型控制方程的边值初值问题,求解过程中把均相化学反应物质中按照反应物和生成物分开,非均相反应物质按照固相和液相分开,对均相反应物及非均相液相物质浓度耦合求解,而均相生成物和非均相固相物质独立求解。使方程组按照其不同类型进行分类,同时可减少未知数的个数。对于含有非线性内状态变量的右端项进行迭代求解。数值例题结果验证了所提出的数值方法的有效性、计算精度和稳定性。     

地震加速度动态分布及对高土石坝坝坡抗震稳定的影响

《水工建筑物抗震设计规范》（DL5073-1997）建议的土石坝地震加速度动态分布系数适合于150 m以下的土石坝,而前许多土石坝的设计高度已远大于150 m。与低坝相比,高坝受地基刚性的约束减弱,坝体的自振周期延长。在坝体地震反应中,高阶自振周期与地震卓越周期耦合的机率增大,高阶振型的振动易被激发放大,从而导致坝体地震加速度沿坝高的分布与低坝有所不同。据此,采用有限元法研究了高土石坝的加速度分布,提出了高于150 m土石坝的地震加速度动态分布系数图示。在此基础上,利用堆石材料的非线性强度准则,对高度超过150 m的土石坝坝坡抗震稳定性作了进一步的分析,得出随着坝体地震加速度动态分布系数的降低,坝坡临界破坏的安全系数有所提高。     

多层弹性平面问题解的精确刚度矩阵法

从平面弹性力学的基本方程出发,利用Fourier积分变换等数学手段,推导出了单层平面问题的刚度矩阵,然后按有限元法组成总体刚度矩阵。通过求解由总体刚度矩阵所构成的代数方程和Fourier积分逆变换,得到在任意静荷载作用下多层弹性平面问题的精确解。由于刚度矩阵不含有正指数项,计算时不会出现溢出现象,从而克服了传递矩阵法的缺点。推导过程中摒弃了应力函数的选择,使得问题的求解更加合理化,同时也为进一步研究此类问题如温度场、动力学等奠定了理论基础。计算实例证明了推导结果的准确性。     

桩-土间剪切应力的传递探讨

对不同情况下桩-土间剪切应力进行了初步的理论分析,并基于ADINA有限元程序,对不同情况下的桩-土间剪切应力进行了计算对比分析。结果表明,同种土体和桩体条件下,各种模式达到极限状态时桩-土间剪切应力不同,堆载作用下负摩擦力值最大,抗拔桩负摩阻力值最小,桩顶受压桩、托桩的摩阻力值介于其间;桩材泊松比在0.15～0.35范围内时,桩侧摩阻力随泊松比的增大而增大。所得结论有助于进一步认识桩-土相互作用规律以及桩-土荷载传递机制。     

面板堆石坝面板竖缝特性对面板应力变形影响分析

采用三维有限元方法,研究了面板竖缝特性对面板堆石坝中面板应力变形的影响。基于类似子结构法概念,建立了精细模拟面板及其分缝特性的三维有限元分析方法。针对马来西亚巴贡面板堆石坝,计算了多种面板竖缝宽度和填料的方 案,得到坝体和面板的应力变形。计算结果表明,面板竖缝几何尺寸和填料特性对面板的挠度、顺坡向应力、水平向压应力和周边缝张合变形影响不大,而对面板的水平向拉应力和竖缝压缩量影响很大;面板竖缝填料变软和宽度增大均会导致面板的水平向压应力减小,而竖缝压缩量明显增大;采用软木作为填充材料较之填充硬木,可以大约降低一半的压应力;面板竖缝宽度及其填料特性对同期填筑的面板、特别是邻近的面板水平向压应力影响较大。