含裂隙岩体单轴压缩裂纹扩展机制离散元分析

将由室内试验总结得到的岩石微观胶结模型嵌入离散元软件,对Lac du Bonnet花岗岩石进行预制单裂隙单轴压缩试验DEM数值模拟,分析了压缩过程中裂隙试样中应力的分布,并与理论计算结果进行对比分析,同时对各种断裂判据中裂纹起裂角的预测值进行了适用性的对比分析。结果表明,离散元模拟试样破坏形态与试验结果相近;离散元分析得到的应力分布与理论解在定性上相似;当预制角度较小时,侧向应力都处于拉压状态;由于裂隙左右两端压应变的集中造成了裂隙上下面拉应变的产生,造成了裂隙周围特殊的应力分布;当裂隙角度较大时,应力集中现象已不明显,因而,理论值与试验值有偏差;在断裂判据中最大周应力准则和最大能量释放率准则得到的裂纹扩展角与室内试验与DEM结果中的数值较为吻合。     

不同开采方法下深海能源土离散元模拟

天然气水合物以胶结及孔隙填充等形式存在于深海能源土中,开采时因其分解会劣化地层力学特性进而引发海底事故,使得人们对能源土开采过程进行中力学特性的变化愈发重视。在前期室内试验的基础上,将一个温度-水压-力学二维微观胶结模型引入离散元商业软件PFC2D中,通过对排气、排水性较好的土体进行升温及降压法开采进行数值模拟,并将模拟结果与相同条件下的室内试验结果对比,验证了该胶结模型的适用性。进一步分析了颗粒接触分布与颗粒平均纯转动率（averaged pure rotation rate, APR）在水合物分解时的变化情况。升温分解过程中随温度升高,颗粒总接触分布各向异性程度增大;胶结接触逐渐减少并始终保持主方向为水平方向,无胶结接触增多并始终保持主方向为竖直方向;APR值逐渐增大且正负值分布逐渐趋于集中。降压分解过程中随反（水）压降低,颗粒总接触由各向同性分布逐渐发展为主方向为竖直方向的各向异性,APR值较小且分布均匀;恢复反压后,试样进一步破坏,颗粒总接触各向异性更加明显,APR值增大且正负值呈集中分布。     

天山乌鲁木齐河源1号冰川致冷效应的小气候观测

高山冰川以其下垫面的致冷效应形成独特的冰川小气候.为研究冰川小气候特征,2007年7月在天山乌鲁木齐河源1号冰川表面及末端冰碛上架设5台自动气象站,并进行了为期一个月的基本气象要素的观测.以观测数据为基础,描述和分析了与冰川致冷效应有关的冰川区温度与湿度变化特征、冰面逆温、温跃现象、冰川风现象,并就冰川致冷效应对局地对...     

基于Morphing的三维地质界面生成

利用Morphing技术,参照四周地质界线的形态特征并以之为约束,插值生成一系列形态渐变的过渡曲线,采用轮廓线算法构建相邻曲线间的不规则三角网,并最终生成较为光滑的三维地质界面。将实际复合问题转化为4种基本问题进行求解,并详述了四边问题的实现细节。实验表明,该方法能够生成光滑三维地质界面,并较好地保留了已知地质界线的形态特征。     

砂土中单桩竖向抗压承载机制的离散元分析

将单桩视为平面问题,采用二维离散单元法分析了砂土中单桩竖向抗压承载机制。首先利用离散元模拟地基的形成、挖孔灌注桩的成桩过程以及地基土的双轴压缩试验;然后通过离散元模拟单桩竖向抗压静载试验,分析其承载机制。结果表明：随着桩的沉降增大,桩端阻力一直增大,桩侧摩阻力先增大后达到稳定值;单位桩长侧摩阻力从桩顶到桩端呈非线性增长趋势;当荷载达到桩的极限承载力时,桩端发生刺入破坏,导致桩侧与桩端附近土体的转动场梯度、应力场梯度增大,主应力发生大角度旋转。     

单粒组密砂剪切带的直剪试验离散元数值分析

为分析砂土的剪切力学特性,采用离散元商业软件PFC2D对单粒组密砂在直剪试验中出现的剪切带进行了数值分析。改进传统的分条带观察方式,采用矩形格分割法观察试样内部不同区域的变形特性。对离散元商业软件进行了二次开发,实现了大小主应力及其应力主方向角的可视化,分析了试样内部的应力偏转情况。同时,以纯转动率和颗粒速度等微观变量为中心,观察了试样内部颗粒的运动状态,解析了直剪试验中砂土剪切带形成的微观机制。研究表明,直剪试验中应变局部化区域集中在剪切面附近的一个条带内。通过转动场和速度场的分析可知,剪切带宽度约为10～15倍的砂土平均粒径,带内颗粒转动明显,带边缘出现大的速度和转动变化梯度。对组构和接触力分布的分析可知,剪切过程中粒间接触点和接触力的主轴方向发生了相一致的偏转,偏转后的主轴方向为60°左右。     

孔隙填充型能源土的宏微观力学特性真三轴试验离散元分析

天然气水合物沉积物因其作为绿色新型能源且具有广阔的开发前景而备受全球瞩目。水合物在水合物沉积物（俗称能源土试样）中有不同的赋存型式,如孔隙填充型和胶结型等。针对孔隙填充型水合物的赋存形态,生成特定饱和度的能源土试样。开展了同一 平面上不同中主应力系数（b = 0、0.25、0.50、0.75、1.00）的真三轴排水试验的三维离散元模拟,将微观参数接触组构及剪切滑移率的演化规律与材料的宏观力学响应相结合,分析了中主应力对孔隙填充型能源土的宏微观力学响应的影响。结果表明：强接触处组构张量的大、中、小主值随b的变化规律与大、中、小主应力及大、中、小主应变随b的变化规律表现出良好的相关性。三向应力状态下的破坏强度接近Lade-Duncan破坏准则。能源土试样剪切滑移率随着中主应力系数的增大而增大;当处于三轴拉伸状态时,试样的剪切滑移率最大。     

基于TIN的体布尔算法及其地质应用

提出一种基于TIN的体布尔操作算法,即采用快速碰撞检测方法、改进的三角形求交算法和三角形包含测试算法,能解决复杂地质体的交、并、差等布尔操作。该方法不仅可以有效处理具有侵入、交错等复杂空间关系的地质体建模问题,还适用于基于地质模型的工程分析应用,如开挖模拟等,并在北京市某岩土工程信息系统建设项目中得到了检验。     

结构性砂土海床管道上浮破坏的离散元模拟

借助离散单元法（DEM）建立海底管道上浮模型,模拟分析了不同胶结强度下胶结试样的力学特性;通过微观力学响应（局部变形、胶结点破坏分布和位移场）对海底管道上浮破坏机制进行分析;对比分析了海床砂土胶结强度对管道上浮土体破坏的影响（上浮抗力和地表隆起量）。结果表明：伴随土体变形和胶结破坏,上覆土体的渐进破坏主要集中在上覆土体区域内;土表隆起量随着管道上浮位移的增长而增长;胶结强度较高时,管道周围土体无回填行为,上覆土体有拉裂带形成,仅管道附近土体产生胶结破坏,峰值上浮抗力随海床砂土胶结强度的增加有显著提高     

钢板桩挡墙主动土压力分布的形状效应

采用有限元数值试验,研究帽型钢板桩的截面形状对墙后主动土压力分布规律的影响。首先,对室内缩尺模型试验进行数值模拟,对比实测数据验证数值模型的有效性;然后,建立钢板桩挡墙数值模型,模拟该挡墙在不同位移模式下土压力变化和分布的规律,与平板挡墙计算结果进行对比,分析钢板桩截面形状对土压力分布的影响,进一步探讨形状效应的可能影响因素及机制。分析结果表明,钢板桩的截面形状影响墙后主动土压力的分布形式,影响程度与位移模式有关;墙体平动和绕墙底转动情况下,钢板桩挡墙凸出部分的主动土压力值大于凹处,但墙体绕墙顶转动情况下差异不明显;主动土压力的形状效应由土拱效应引起,截面高宽比对其影响显著,墙土摩擦角影响有限。