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倍。     

基于LS-DYNA的贯通式潜孔锤反循环钻头强度的优化分析

通体花键式贯通式潜孔锤反循环钻头在实际施工中尾部较易损坏,影响反循环钻头的使用寿命。为提高反循环钻头的使用寿命,对反循环钻头尾部结构进行了优化设计,并采用非线性动力学仿真软件LS-DYNA对改进前后的反循环钻头进行了仿真分析。分析结果表明：由于在冲锤的高频冲击作用下,反循环钻头尾部易出现较大的局部应力集中,导致钻头体出现疲劳破坏;通过对钻头尾部结构设计优化后,与通体花键式反循环钻头相比,新结构钻头在无泄风槽情况下,应力集中值减小44%,在有泄风槽的情况下,应力集中值减小19%。数值模拟分析表明,优化后能大幅度提高钻头寿命。     

频谱激电法的发展与展望

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

水库异重流一维水沙耦合模型

建立了基于库区不规则断面的一维非恒定异重流数学模型,并采用明流与异重流水沙输移模型交替运算的两步模式,即用潜入条件动态判别异重流计算的上游边界位置,将潜入点上游的明流浑水段与下游异重流段计算连接起来。水流运动、泥沙输移与河床变形过程完全耦合,采用TVD(Total Variation Diminishing)形式的MUSCL-Hancock格式进行数值求解。将该模型应用于恒定流量与释放定量悬沙两种条件下的异重流水槽实验模拟,比较了有无水面梯度项对模拟精度的影响,计算结果表明该模型能较为准确地预测异重流的厚度、含沙量分布及传播过程。     

路堤荷载下碎石桩处理软基沉降计算

以沉降控制为目的的碎石桩复合地基,其沉降计算在设计、施工中具有重要地位。本文基于典型的桩土单元体模型,建立了复合地基体积压缩系数与复合地基弹性模量的关系,考虑施工中的时间因素及相应的固结度对沉降的影响,提出相应的沉降计算公式,并结合四川省遂-资高速公路软基变形监测的数据进一步验证其可行性。研究表明本算法计算值比实测沉降约大10% ～20%。填筑过程中,计算沉降历时曲线与实测曲线拟合较好,更接近实测值;在此基础上,初步探讨了影响路堤荷载下碎石桩处理软基沉降变化的因素,得出桩径、桩间距对其影响较为明显。同时结合本例,建议布桩时桩间距与桩径比值最好控制在2～4之间。     

顾及多要素影响的路网匹配改进算法

针对导航地图中现有的匹配算法在复杂路段的匹配效果一般较差的问题,本文提出一种改进的匹配算法。首先,引入分块思想,提出将路网中的每个路段按规则分为简单路段块和复杂路段块,不同路段块采取不同的匹配算法;其次,充分利用车辆行驶的连续性,综合考虑了车辆行驶方向、车速与交叉口距离等多种因素对匹配结果的影响,在选取正确的匹配路段时合理利用了历史数据信息。实验表明,该算法不但对能够改善复杂路段的匹配效果,而且还能够对匹配点沿路段方向上的误差进行实时校正。