三维直流电阻率自然单元-无限元耦合数值模拟

直流电阻率数值模拟对于实测电阻率数据反演、解释地下真实电阻率结构具有重要意义.为解决复杂地电模型剖分以及处理截断边界时面临的困难,本文提出一种新的三维自然单元-无限元耦合法.该算法在自然单元区域灵活布设自然节点以填充复杂模型,并引入无限单元替代截断边界条件以促使刚度矩阵与电源位置无关.首先推导了稳定电流场基本方程;然后介绍了二、三维自然单元法及其形函数的构建过程;再引入三维多向映射无限单元并实现自然单元-无限元的有效耦合;最后通过4个数值模型,结合传统有限单元法、有限元-无限元耦合法以及COMSOL Multiphysics软件的计算结果,验证了算法的正确性及有效性.数值结果表明,本文所提算法满足精度要求,能有效减小计算区域,无需二次更新刚度矩阵,对复杂模型适应性强.     

基于数字图像测量技术的粉状煤系土微观结构分形特性分析

以武深高速广东段沿线的粉状煤系土为研究对象,利用环境扫描电镜技术(ESEM)获得了不同含水率直剪试验后粉状煤系土剪切面的微观结构SEM图像;结合MATLAB及Image Pro Plus(IPP)软件,对剪切面微观结构特征进行了分析;基于分形理论,建立了煤系土的分形模型,求出了二维空间内煤系土孔隙轮廓分维数、孔隙数量～孔径分布分维数。结果表明:煤系土微观结构多为片状颗粒集合体,接触关系主要为面-面接触和面-边接触;随着含水率的增加,剪切面粗糙度先增加后减小,力学强度参数先增大后减小,转折点在最优含水率附近(10%～15%之间);煤系土微观结构具有明显的分形特征,可用孔隙等效面积-等效周长分形模型、孔隙数量～孔径分布分形模型描述,其分维数介于1～2。