N-S方程基于投影法的特征线算子分裂有限元求解

提出了一种求解非定常不可压缩纳维-斯托克斯方程（N-S方程）的新型有限元法:基于投影法的特征线算子分裂有限元法.在每一个时间层上将N-S方程分裂成扩散项、对流项、压力修正项.对流项采用多步显式格式,且在每一个对流子时间步内采用更加精确的显式特征线-伽辽金法进行时间离散,空间离散采用标准伽辽金法.应用此算法对平面泊肃叶流、方腔流和圆柱绕流进行数值模拟,所得结果与基准解符合良好.尤其对于Re=10000的方腔流,给出了方腔中分离涡发展和运动的计算结果,并发现在该雷诺数下存在周期解,表明该算法能较好地模拟流体流动中的小尺度物理量以及流场中分离涡的运动.      

土工织物补强软基的整体稳定性计算

针对现有的应用于有加筋垫层路堤的整体稳定性计算的分条圆弧法,未能反映筋材所起的全部作用,偏于保守难以全面反映加筋效果.为了从加筋机理和计算分析两方面更合理地阐述加筋垫层的加筋效果,利用现场足尺破坏试验数据,结合边坡稳定分析的圆弧滑动条分法,分析了加筋垫层路堤的失稳破坏机制.在此基础上,指出深层滑动潜在滑弧要通过路堤的沉降中心点, 并为加筋垫层路堤的稳定性分析提出了一种合理考虑加筋垫层阻滑作用的简化实用分析方法.该方法可推广应用于工程实践.     

裂纹岩体贯通路径的有限元搜索

裂隙岩体在荷载作用下导致其内部裂纹的相互贯通,从而在裂纹尖端形成较大的塑性破坏区域,岩体工程的失稳大多是由于其内部节理、裂隙等缺陷发展导致的,应用断裂力学的方法,可以追踪岩体中节理裂隙的起裂、扩展至相互贯通使岩体局部破坏的过程,从而揭示岩体失稳的渐进破坏机制.应用有限元方法中的奇异单元模拟了岩体中不同倾角的压剪裂纹的破坏问题,根据不同倾角裂纹的力学机理分析了其受力差异,依据最短塑性区距离判据分析了相近裂尖的扩展距离,得出了不同倾角的两共线裂纹的破裂角及相近裂尖破坏后岩桥之间的有效距离.     

岩石破裂过程THMD耦合数值模型研究

从岩石的细观非均匀性特点出发,应用损伤力学、热力学和渗流力学理论,建立了岩体热(温度)-水(渗流)-岩(应力)-损伤耦合数值模型(THMD model),把岩石(体)THM耦合问题的研究从应力状态分析深入到损伤、破坏过程分析之中。探讨了THM耦合作用下岩石材料的细观结构损伤及其诱发的材料力学性能演化机制,并运用所提出方法计算温度-渗流-应力耦合作用下井筒近场围岩的稳定性,模拟得到的岩体破坏过程、应力分布、AE特性及渗流特性变化与现场标定结果有着一致的规律性,初步证明了该数值模型的合理性和有效性。THMD模型以简单的数值模型表征了岩石(体)中热、水、岩及损伤之间复杂的作用关系,为从细观损伤演化揭示宏观岩体温度-渗流-应力耦合破坏机制提供了一种新的数值分析方法。     

波浪对船剖面作用的时域计算模型

为了求港口中波浪对船的作用力,用边界元数值求解复速度势建立了波浪对船剖面作用的时域计算模型,采用满足运动学和动力学条件的半拉格朗日法跟踪流体的自由表面。推导了适用本数值模型的辐射边界条件,此边界条件和海绵层边界条件组成的出口边界条件,可以使波浪很好地被吸收或传播出去。数值计算结果和实验结果对比表明此模型能精确模拟波浪对船的作用力。     

流固耦合作用下节理岩质边坡失稳过程的RFPA模拟分析

利用基于强度折减法的RFPA-Slope对渗流与应力耦合作用下的软弱互层岩质边坡稳定性进行了数值模拟分析,数值模拟不但直观形象地给出了边坡的渗流场、应力场、破坏区分布,而且得到了边坡滑移破坏面的萌生、扩展、贯通以及坡体整体失稳的渐进破坏过程,同时求得安全系数。并与无地下水的稳定性作了比较,结果表明,地下水渗流的作用使坡体位移增大,边坡安全系数减小,明显加大了滑坡范围。对实例的分析说明,RFPA-Slope能够较为准确地预测边坡潜在破坏面的形状与位置及计算相应的稳定安全系数,本文方法对于边坡,特别是对于复杂边坡的稳定性分析具有实用性。     

潮汐循环作用下海水入侵模型的现状及展望

针对潮汐循环作用下海水入侵模型的问题,综述了目前国内外海水入侵模型的研究现状,探讨了流固耦合效应和潮汐时域运动效应对海水入侵的影响,指出了建立潮汐循环作用下海水入侵的渗流-应力-损伤-对流扩散多场耦合模型,并考虑潮汐循环加载、岩土非均匀性和岩土损伤与对海水入侵之间的关系,是今后海水入侵模型的主要研究方向。     

深部岩体圆形巷道围岩松动圈形成机理数值试验研究

采用损伤力学和统计理论的单元本构模型,探讨深部岩体巷道围岩破坏规律, 利用岩石破裂过程分析软件RFPA$^{\rm 2D}$, 对深部 岩体中的圆形巷道的变形及非线性渐进破坏特征、巷道周边关键部位的位移和应力变化进行 了分析,研究了圆形巷道围岩松动圈的形成机理. 研究表明：巷道开挖后,其周边形成应力 集中带,随着围岩压力的持续作用,巷道周边产生塑性变形区,并沿径向形成裂纹;随着裂 纹的不断扩展,巷道周边出现松动破裂区,即松动圈. 产生破裂区后,应力集中程度减弱, 应力高峰点向远处转移.     

基于细观力学的混凝土弹塑脆性损伤研究

将混凝土假定为一种由硬化水泥砂浆、粗骨料、界面粘结带所组成的三相复合材料,在满足骨料级配曲线算法的基础上,采用细观单元的弹塑脆性损伤本构关系,考虑材料的非均质特性,建立了基于细观力学的混凝土弹塑脆性损伤数值模型;分别研究了单轴受拉预置裂纹试样和单轴受压混凝土试样的细观弹塑脆性损伤破坏行为,并揭示了混凝土的宏观表征强度存在明显的尺寸效应,通过将计算结果与 Bazant 尺寸效应公式、单轴受压物理实验曲线进行对比,验证了模型的正确性。数值试验表明：该模型可以清晰地模拟混凝土细观塑性屈服和失效裂纹的萌生和扩展。骨料与水泥砂浆间的界面粘结带相对薄弱,在混凝土试件形成宏观损伤局部化带前,试件的屈服和破坏首先发生在骨料边缘处的界面位置,并沿着界面粘结带扩展、贯通;同时,导致宏观裂纹形成和发展的因素仍以细观单元的拉伸破坏为主。     

锚索预应力损失的影响因素分析及其补偿措施

为了研究锚索预应力损失的影响因素以及对其的补偿措施,对多个工程锚索的应力测试结果进行了较详细的分析和研究,得出了影响锚索预应力损失的因素:锚索材料钢绞丝的松弛、锚固对象岩土体的变形以及施工中的张拉、锁定等原因.依据力学性能优选适宜于预应力锚索的新材料,并在上述分析结果的基础上提出了补偿锚索预应力损失的措施,提高了预应力锚索体系的综合锚固效果.对预应力锚索工程具有一定的理论和指导意义.