岩土类软化材料的柱形孔扩张统一解问题

用应力一次跌落应变软化模型来模拟岩土软化材料特性，采用双剪统一强度理论的屈眼函数形式，推导了柱形孔扩张问题应力场及位移场的统一计算公式，分析了不同模型、不同软化程度对柱形孔扩张时应力、位移及塑性区开展规律的影响。     

BDS空间信号监测精度分析

为了验证BDS空间信号精度,借鉴Galileo空间信号监测精度SISMA的原理,采用亚太区域仿真的测站数据评估BDS空间信号性能。分析结果表明,利用亚太局域完好性监测网络对经过亚太上空的BDS在轨健康卫星的SISMA监测结果,满足生命安全服务SOL(Safety-Of-Life)要求,可信度为100%(除C12是97.92%以外)。     

动渗耦合作用下软黏土渗透系数演化机制分析

针对天津地区滨海软黏土,开展了一系列动渗耦合三轴试验,系统分析了长期循环作用下软黏土渗透系数的变化规律。结合扫描电镜（SEM）和压汞试验（MIP）对不同变形阶段软黏土的微观结构特征进行分析,阐述了渗透系数的演化规律及机制。研究表明：在动渗耦合作用下,软黏土渗透系数呈现三阶段变化特征;渗透系数随渗流力呈现先增大后减小的变化规律,随动应力和频率呈现先减小后增大的特征。微观试验表明,渗透系数的演化是微结构孔隙大小、形态、分布特征渐进调整的宏观体现：加载初期,颗粒间的架空超大孔隙（孔隙直径D>2.5 μm）大幅压缩,孔隙形态变化不大,渗透系数呈线性衰减;加载中期,大颗粒破碎生成的小颗粒密布填充于大颗粒间,生成大量“密小而曲折”的小孔隙（0.05 μm<D<0.1 μm）,渗透系数呈现缓慢下降;加载末期,结构致密稳定,渗透系数基本不变。对于非稳定变形,颗粒团（聚）体间重新生成“条带”状孔隙,渗透系数反弹增大。该研究成果可为流固耦合分析中渗透系数的确定提供理论依据。     

胶结颗粒接触力学特性测试装置研制

为验证天然结构性砂土离散元模拟中接触模型及其参数的合理性,设计了一套用于理想胶结颗粒成型及实现不同加载条件下接触力学特性测试装置。通过胶结颗粒成型装置在两大小相同的铝棒间形成具有特定几何尺寸的胶结物,随后,采用一系列辅助加载装置实现简单及复杂加载条件下胶结颗粒接触力学特性的测试。试验结果表明：该装置可用于胶结颗粒在不同加载条件下接触力学特性的测试,实测胶结颗粒接触力学响应与天然砂土离散元中接触模型基本相符,且其抗剪和抗扭强度均随着法向压力的增大而增大,在三维应力空间中胶结颗粒强度包线呈椭圆抛物面状。     

考虑简化胶结模型的深海能源土宏观力学性质离散元数值模拟分析

根据蒋明镜等所提出天然结构性砂土微观胶结模型[1–2]及微观胶结试验结果[3–5],将该模型引入离散元商业软件PFC2D,进行能源土双轴试验离散元数值模拟分析,并同Masui等[6]能源土三轴试验结果进行对比分析,结果表明,蒋明镜等所提出胶结模型能够较好的模拟水合物的微观胶结力学行为,水合物胶结的存在对能源土强度具有一定的贡献。     

基于颗粒抗转动模型的刚性挡墙被动土压力临界状态离散元分析

简要介绍了颗粒抗转动模型,并将其引入离散元程序中,通过建立挡墙地基模型和合理选取模型参数,分别考虑了地基填土不同密实度和挡墙不同位移模式（被动T模式、RB模式、RT模式）情况下,刚性挡墙被动土压力随挡墙位移增长发展到达临界状态时,土压力系数 随位移发展的变化规律及墙后填土剪切带的形成规律,并与其他学者的研究成果进行对比分析。研究结果表明,土压力系数 随着挡墙位移增长的变化规律与填土的孔隙比（或相对密实度）和挡墙的位移模式紧密相关。随着孔隙比的减小或相对密实度的增大,土压力系数 会逐渐由位移硬化特性过渡为位移软化特性。尽管中密试样在双轴压缩试验中呈现出应变软化特性,而中密样的土压力系数 随着挡墙平动位移的增长可能呈现出位移软化特性,也可能呈现位移硬化特性。随着刚性挡墙向墙后土体推移,试样中的剪应变随之增大,并会在墙后形成应变局部化,即剪切带的出现。与室内试验剪应变云图相似,离散元较好地模拟了土压力临界状态时剪切带分布规律。同时,墙后土体表面不再是光滑的平面,而是逐渐隆起的凹凸面;随着挡墙位移增长,土体表面隆起量越来越大,直至土体破坏。     

人工岛沉降的离心机模型试验与数值分析

采用离心机模型试验与数值方法,对某一人工岛深厚软黏土地基的沉降特性进行了研究。在模拟现场施工条件下,通过离心机模型试验,研究了60 m厚软黏土地基的沉降特性。基于离心机模型试验,建立了二维有限元模型,应用弹黏塑性本构模型对软黏土地基的沉降规律进行了分析。通过模型试验与数值计算结果的对比分析,结果表明,模型试验与数值计算得到的工后沉降曲线基本吻合,初期沉降速率较大,后期沉降速率相对稳定。在浅部土层,计算得到的超孔压消散快于模型试验结果,而深部土层两者超孔压消散规律基本一致。二者规律的一致性验证了所用数值方法的可靠性。     

结构性黄土一维湿陷特性的离散元数值模拟

针对结构性湿陷性黄土大孔隙和胶结特性,应用离散元生成了不同含水率结构性黄土试样,研究试样的一维湿陷特性。首先,根据已有的结构性黄土试验资料和胶结颗粒材料离散元数值试验成果,建立胶结强度和初始饱和度之间的关系。其次,采用蒋明镜等提出的分层欠压法[1]和胶结模型[2]制得不同含水率结构性黄土离散元试样,然后进行不同含水率双线法和同一含水率4个压力下单线法湿陷试验的离散元数值模拟。数值模拟结果表明,提出的离散元分析方法能模拟天然结构性湿陷性黄土的主要力学性质,随着含水率的减少,结构屈服应力和最大湿陷压力增加,湿陷系数随着压力先增加后减小,湿陷起始压力为饱和试样的结构屈服应力,单线法湿陷后压缩曲线与饱和试样的压缩曲线接近。此外,模拟结果还表明,不同含水率结构性黄土离散元试样的最大湿陷系数与天然结构性湿陷性黄土相差较远,但在最大湿陷系数与孔隙比的比值上相接近;结构屈服对应着胶结的逐步破坏,湿陷伴随着大量的胶结破坏。提出了基于胶结点数目的损伤变量,研究了其在加载和湿陷过程中的变化规律。研究成果为认识黄土复杂力学特性和建立其本构理论提供了基础。     

一种静力触探模型箱试验装置的研制及其试验研究

设计了一套静力触探模型箱试验装置,该装置可以在模型箱水平面内任意位置进行贯入试验,贯入速率在一定范围可控,贯入倾角可调节。将该装置用于TJ-1模拟月壤静力触探力学特性的研究,试验结果表明：地基表层模拟月壤发生刺入剪切破坏;贯入阻力随贯入深度的增加而增加;探杆周围同一深度处土体的水平土压力表现出轴对称特征,其值在探头靠近土压力盒过程中逐渐增大至最大值,后随探头远离土压力盒又逐渐减小至0,且不同深度土体水平土压力峰值随贯入深度增加而增大;贯入点正下方土体的垂直土压力随贯入深度的增加而增加,主要影响范围约为距探头200 mm处。基于以上试验研究,表明该模型箱的设计是合理的,且TJ-1模拟月壤力学特性的研究成果能够为将来涉及到TJ-1模拟月壤的航天试验提供必要的技术参考。     

应力路径试验前后不同黄土的孔隙分形特征

为了解不同黄土孔隙形状复杂程度的差异和应力路径对孔隙形状复杂程度的影响,对两种黄土应力路径前后孔隙分形特征进行了研究。首先比较了3种分维模型所得孔隙分形维数的可靠性,然后选用热力学关系模型,由进汞、退汞试验得到地裂缝区黄土、充填黄土初始样和三轴应力路径试验后的孔隙分布,据此分析了两种土体进汞孔隙和退汞孔隙分形维数的差异和受载后的分形维数变化,根据退汞过程仅管形孔内的汞流出及孔隙由管形孔和球形孔组成的假定,得到了应力路径试验前后孔隙形状的改变。结果表明,热力学关系模型得到的孔隙分形维数合理、可靠。初始状态,原状和重塑充填黄土的孔隙形状比原状和重塑黄土复杂;常规三轴试验后试样孔隙比减少,孔隙分形维数增加,球形孔向管形孔转变;相对于常规三轴压缩试验,减围压三轴压缩试验后试样的孔隙分形维数较小,管形孔占总孔隙的体积比例较少;总体上管形孔的分形维数比球形孔大,且基本不受应力路径的影响,其占总孔隙的体积比例随着试样宏观孔隙比的减少而增加。