直流电阻率法2.5维正演的外推瀑布式多重网格法

引入外推瀑布式多重网格法(EXCMG)求解2.5维直流电阻率有限元计算形成的大型稀疏线性方程组,结合基于地址矩阵的压缩存贮方式以及最优化离散波数,使得2.5维电阻率正演程序的计算速度大大提高而内存需求大大减小.研究结果表明:EXCMG法的收敛速度与网格尺寸无关,计算速度明显优于不完全Cholesky共轭梯度(ICCG)...     

考虑非线性徐变的混凝土温度裂缝扩展过程模拟

混凝土温度应力超过相应龄期的抗拉强度一半时,混凝土的徐变变形不再与应力呈线性关系。采用一个混凝土非线性徐变模型来进行混凝土的温度应力计算,并给出了相应的算法。同时,吸取无网格方法在裂缝扩展模拟方面的优势,采用无网格方法对混凝土的温度裂缝过程进行了仿真计算。数值算例计算结果表明,所编制的无网格法程序具有较好的收敛性和可靠性,能够有效地模拟混凝土块在温度应力作用下的开裂发展过程,考虑混凝土非线性徐变效应后,在开裂时刻混凝土急剧增加的非线性徐变变形缓和了缝端的温度拉应力,使得相同龄期混凝土的裂缝深度和扩展速率要比只考虑混凝土线性徐变要小一些。     

砂土液化变形的有限元-无网格耦合方法

通过构造Biot固结理论u-p方程的无网格伽辽金-有限元耦合方法,对砂土液化变形问题进行了数值模拟。对于饱和砂土,采用Oka等提出的弹塑性本构模型,同时采用更新的Lagrange计算格式推导了控制方程。耦合方法能够发挥有限元和无网格各自的优点,既避免了由于单元变形扭曲而引起的计算中断,也可节约计算时间,算例验证了该方法在地震液化问题中的有效性。     

基于无网格法的含缺陷岩石变形局部化数值模拟研究

利用自主开发的基于无网格法的岩石力学计算分析程序,采用Mohr-Coulomb破坏准则,考虑材料塑性屈服后的软化特性,研究单轴压缩条件下含一个缺陷的岩石试样的变形特性,探讨含缺陷材料的岩石变形局部化的演化规律。计算结果表明,岩石试件的破坏是一个渐进破坏的过程;加载过程中试件中的缺陷弱粒子首先出现粒子破坏并产生声发射,随加载进行试件产生微裂隙并扩张形成一个变形局部化带,局部化带扩张、发展,最终形成一个明显的剪切破坏带,缺陷粒子对试件的破坏起到了控制作用。岩石的声发射记录了岩石材料的塑性破坏过程,岩石试件的剪切破坏带与塑性区的发展具有非常类似的规律,但塑性区面积要大于剪切破坏带,从所处位置来看剪切破坏带位于塑性区之内,通过研究声发射特征信息发现变形局部化发生在应力-加载步（应变）达到峰值强度前。     

基于三维非线性有限元的边坡稳定分析方法

刚体极限平衡法不能反映岩体中实际的应力分布,而基于有限元的强度折减系数法在判断收敛性方面存在一些问题。为了解决这些问题,采用多重网格法,分别建立用于有限元计算的结构网格和用于计算滑面稳定安全系数的滑面网格,可以方便地获得任意滑面或滑块的稳定安全系数,从而将非线性有限元和极限平衡分析结合起来。为了提高计算规模和计算精度,采用有限元并行计算程序TFINE.Pfem进行计算,分析了网格密度对计算结果精度的影响,并应用于锦屏高边坡的稳定分析中。与刚体极限平衡法结果的对比分析表明,由于考虑了计算过程中的非线性应力调整,该方法的计算结果比刚体极限平衡法偏大,而且更符合实际情况。     

模拟摩擦接触问题的新型无网格数值方法

针对摩擦接触问题,给出一种新型的无网格数值方法。该方法基于单位分解思想,在标准无网格Galerkin法的位移模式中嵌入不连续函数和裂尖奇异函数,分别反映接触面的不连续性和接触面端点的奇异性;结合接触摩擦定律,继而构造出求解摩擦接触问题的无网格线性互补模型。在该方法中,接触面方程使用接触点对离散,全局离散系统方程则转化为标准的线性互补问题,可以很方便地使用Lemke算法求解。算例分析证明,本文方法的正确性和有效性。     

试验模型位移场中的几种光学测量方法比较

详细地分析了几种光学测试方法在模型试验中的使用情况。自动网格法在模型表面布置网格点容易,由于开挖或裂缝等原因去掉网格点,改进后的程序仍能自动编码,必要时可通过人机对话手工编码。散斑互相关法,在模型表面布置散斑较容易,改进后的程序在开挖情况下也能自动处理,但在散斑场遭受严重破坏情况下无法比较相关性。其他光学测试技术由于受试验条件的限制,而不适合使用。     

多重网格法在MM5自适应模式中的应用

将多重网格法应用于MM5中来加速自适应网格的生成,并用这种方法来模拟2002年6月17-18日发生在江南和华南的暴雨过程。试验表明：在收敛精度要求高的情况下自适应模式模拟结果的精度高于均匀网格模式,且多重网格法有助于提高运算效率。     

海岸线分维数计算方法及其比较研究

海岸线是分形领域最传统的研究课题 ,但是在目前的海岸线研究中却忽略了一个重要且根本性的问题 ,那就是海岸线分维数不同计算方法之间是否存在着一定的差异 ,根据不同计算方法计算所得的分维数之间是否具有可比性。在概述海岸线分维数计算方法的基础上 ,分析了它们之间的异同 ,并进行了实证研究。指出 :使用量规法、网格法计算所得海岸线分维数结果之间存在着客观的差异 ,文中利用量规法计算所得的海岸线分维数结果大于使用网格法计算所得的海岸线分维数结果。     

海岸线分维数计算方法及其比较研究

海岸线是分形领域最传统的研究课题，但是在目前的海岸线研究中却忽略了一个重要且根本性的问题，那就是海岸线分维数不同计算方法之间是否存在着一定的差异，根据不同计算方法计算所得的分维数之间是否具有可比性，在概述海岸线分维数计算方法的基础上，分析了它们之间的异同，并进行了实证研究，指出：使用量规法、网格法计算所得海岸线分维数结果之间存在着客观的差异，文中利用量规法计算所得的海岸线分维数结果大于使用网格法计算所得的海岸线分维数结果。