浅基础地基荷载-变形特性的分析研究

在土体主应力方向的旋转过程中,主应力方向与塑性主应变增量方向之间存在着非共轴现象,非共轴现象对应力-应变关系具有显著的影响。通过有限元二次开发,将角点结构非共轴理论应用到有限元程序ABAQUS中;对密砂单剪试验进行数值模拟,研究了非共轴现象及其对应力-应变关系的影响;采用离心模型试验方法,对圆形浅基础作用下饱和密砂地基荷载-变形特性进行试验研究;对离心模型试验进行数值模拟,将试验结果与数值计算结果进行对比,对非共轴模型的计算结果的合理性进行验证。研究结果表明:非共轴现象对浅基础地基荷载-变形特性具有显著的影响;当选取合理的非共轴塑性模量时,非共轴模型的计算结果与离心模型试验结果比较接近。     

高海拔隧道涌突水预测分析——以四川小卡子隧道为例

小卡子隧道是连接九龙县城和石棉县城的一条便利通道,它是连接S215和S211的绿色保障,九龙和石棉之间是海拔高度4 000 m以上,九龙县城周边为高原宽谷丘陵地区,地势起伏较大。九龙至周边沿线不良地质和特殊褶皱构造发育,其中冰雪害几乎每年冬春都发生,且影响线路较长,其它地质灾害也时常暴发,致使交通事故频繁发生,成为制约九龙县发展,因而急需要对线路进行改造,拟建连接九龙县至石棉县小卡子隧道即为公路改建工程。     

爆炸地震波作用下盾构隧道动力响应分析

盾构隧道衬砌由于各种类型接头的存在而与整体式衬砌的力学特性存在较大差异。将盾构隧道衬砌结构看作由弹塑性铰链连接的刚性管片组成,考虑围岩介质的黏弹性,提出了爆炸地震波作用下盾构隧道动力分析的简化计算方法。采用该方法对南京地铁盾构段典型横断面进行了动力分析,研究了爆炸地震波入射角度、围岩介质特性及管片厚度对结构受力与变形的影响规律。分析结果表明:波入射角度对盾构隧道有很大影响,斜入射时结构的动力响应要大于垂直入射时结构的动力响应;围岩介质等级越高,围岩对隧道结构的约束越强,隧道的抗爆性能越好;管片厚度的增大会增大结构的内力,合理设置管片厚度有利于提高盾构隧道抗爆性能。     

用传输矩阵方法研究多层介质增透膜的透射特性

利用麦克斯韦方程组和边界条件,推导出由高低折射率介质组成周期性的多层增透膜系的传输矩阵,用传输性矩阵求出增透膜系的反射率和透射率表示式,采用Maple语言和数值计算方法得到4层和8层高反射膜系的反射率随入射光波长变化曲线,分析和讨论入射光的入射角对反射率和透射率的影响.得到的结果有利于对多层介质增透膜系结构的理解,为光学仪器镀膜的设计提供了理论依据.     

GLONASS卫星位置的动力学改进算法

为了减弱由于广播星历拟合引起的误差,本文提出了一种基于动力学轨道改进的GLONASS卫星位置算法。该算法仅采用GLONASS广播星历信息就能够提高GLONASS卫星位置坐标的计算精度。算例表明,该方法能够有效减弱GLONASS卫星广播星历的拟合误差,卫星位置的精度损失从0.5m降到mm级;外推2-3h,精度损失也仅在cm级水平。     

基于OpenSees的砂土本构模型对比研究

饱和砂土液化是由地震引起的一种最常见的工程地质现象,也是造成重大地震灾害的主要原因之一。由于成因的复杂性和所造成灾害的严重性,饱和砂土液化一直是土动力学和岩土地震工程研究领域的重要课题。针对饱和砂土液化问题,基于开源地震工程数值计算平台OpenSees,对材料库中的4种砂土本构模型进行数值计算。采用二维u-p单元模拟土颗粒位移和孔隙水压力,分析和对比4种模型在循环动力荷载作用下的加速度、超孔隙水压力、位移、剪应力-剪应变和平均有效应力路径方面的响应结果。研究结果表明:(1)砂土对输入加速度表现出一定的放大效应,对于不同的模型,该放大效应存在一些差异;(2) Stress Density模型在循环动力荷载作用下易产生永久变形;(3)在循环动力荷载作用下,PDMY模型和CycLiqCPSP模型的强度逐渐降低,直到完全消失;(4) Stress Density模型和Manzari Dafalias模型在循环动力荷载下表现出明显的剪胀效应。研究成果对砂土液化的数值模拟问题具有重要的理论价值,可为饱和砂土的液化模拟和砂土本构模型的选取提供参考。     

胶东金矿成矿深度的构造校正测算及成矿预测

成矿深度测算对于矿床学理论研究和深部找矿都有重要意义。经典的成矿深度“压力/比重”计算方法,缺乏考虑构造应力在成矿过程中的影响。前人按照“压力/比重”的计算方法,提出胶东蚀变岩金矿是6 000~8 000 m深的元古宙成矿,石英脉型金矿是深度在3 000 m左右的中生代成矿,并据此建立了金矿垂直五层楼的分带模式。依据该模式指导的深部预测勘查效果不好。“成矿深度构造校正测算”是近几十年逐渐成长起来的一个新方法,即先减去“构造附加压力”后再进行成矿深度测算。本文介绍“成矿深度构造校正测算”的计算方法,指出其应用条件和预测意义。开展成矿深度构造校正测算需要以下条件:(1)确定成矿模式;(2)开展野外构造变形岩相测量;(3)测量岩石矿物应变,恢复成矿构造应力场,计算构造附加静水压力;(4)测算成矿深度。根据“成矿深度构造校正测算”方法已经获得胶东多个典型金矿成矿深度的测算结果:(1)夏甸金矿成矿深度为-1 979.51~-3 014.72 m;(2)焦家金矿成矿深度为-1 632.4~-2 331.6 m;(3)大尹格庄金矿成矿深度为-2 775.4~-4 164.5 m;(4)新城金矿成矿深度为-1 781.29~-2 750.0 m;(5)玲珑金矿成矿深度为-720.55~-3 454.97 m。根据以上典型金矿成矿深度的测算结果,本文认为胶东金矿属于深-1 000~-4 500 m的浅成热液剪切带型矿床,由此推断胶东典型金矿矿体主要部分仍然赋存在深部。按照“构造校正测算”方法得到的成矿深度,结合地质、物探和化探信息,预测金矿发育“深部第二富集带”,已经得到胶东金矿勘查工作的证实。     

热力学计算模拟对初始月幔结构的约束SCIEI北大核心CSCD

月球岩浆洋结晶形成的初始月球内部结构是其后续演化过程的开端,其结晶过程受月球岩浆洋的初始深度和物质组成这两个参数的制约。由于缺少直接来自月球深部的岩石样品,目前关于月球岩浆洋演化过程的探讨主要依赖实验和计算模拟手段。岩浆洋模型中形成的月壳厚度是否与探测结果一致是月球岩浆洋演化模型合理性的重要约束。最新的GRAIL(Gravity Recovery and Interior Laboratory)探测数据推算月壳厚度为34~43km,低于阿波罗时期认为的约70km,这对已有的月球岩浆洋演化模型提出了挑战。本文采用并修正FXMOTR程序包,针对月球岩浆洋在不同的初始深度和物质组成情况下的结晶过程,进行了一系列热力学计算模拟。通过量化月球岩浆洋的初始深度和物质组成对月壳厚度的影响,结合关于月球内部微量元素分配的研究结果,对比了月球岩浆洋结晶后期的残余熔体与原始克里普组分(urKREEP)的成分。本文的模拟结果显示,一个全月幔熔融且初始成分为月球初始上月幔组成(LPUM)的岩浆洋将在其深部结晶2.5%石榴子石,形成的月壳厚度符合GRAIL的约束,并且结晶出了合适的urKREEP成分。在此模型的基础上获取了月球初始的内部成分和密度结构,并对后期月幔翻转(Overturn)的程度进行了探讨。     

分布式水文模型的GPU并行化及快速模拟技术

分布式水文模型对流域水文过程的应用深度及广度不断加深,常与数值天气及气候预报相结合,面临巨大的计算量。近年来GPU技术的进步使普通电脑能够实现高效而又廉价的并行计算。提出了资料插值、单元产流以及单元汇流采用GPU并行计算,马斯京根法河道汇流采用一种非并行的递归方法。基于笔记本电脑和NVIDIA GPU/CUDA结合C#语言,由分布式新安江模型在沂河流域的模拟应用表明,降水量空间插值及新安江产流的并行执行效率为普通CPU上C#的8~9倍。使用直接递归法实现马斯京根汇流演算比以往采用汇流次序表的执行效率提升0.5~0.9倍。     

频谱激电法的发展与展望

从频谱激电法的发展概况﹑仪器系统、数学模型、反演方法及应用等方面,介绍了频谱激电法的研究现状。介绍了频谱激电法目前常用的仪器系统:SIP-FUCHSII和V8,数据模型主要有:Cole-Cole、常相位角模型CPA、普通化的Cole-Cole模型和理论模型SNP。反演方法简要介绍了几种常用算法,反演算法也由一维、二维发展到三维。列举了频谱激电法近年来在矿产资源、水资源调查等多个领域的应用概况,展望了频谱激电法的发展方向:(1)同时考虑激电效应和电磁效应的三维电磁场正演计算技术是研究的前沿和热点;(2)频谱激电法对有机污染的探查成为未来应用研究的新领域。