循环荷载下砂土液化特性颗粒流数值模拟

利用PFC2D常体积循环双轴试验条件,对砂土在不排水循环荷载作用下的液化特性进行了颗粒流数值模拟,数值模拟按等应力幅加荷方式进行。颗粒流数值模拟的优点在于得到试样液化宏观力学表现的同时,通过不同循环加荷时刻试样内细观组构参量（包括配位数、接触法向分布、粒间法向接触力、粒间切向接触力）的演化规律,分析砂土液化过程中细观组构变化与宏观力学响应之间的内在联系,从而可进一步探讨砂土液化的细观力学机制。数值模拟研究结果表明,砂土液化现象在宏观力学表现上反映为超静孔隙水压力的累积上升和平均有效主应力的不断减小,在细观组构上对应于配位数的累积损失和粒间接触力的不断减小。砂土液化细观机制分析表明,试样配位数的减少与循环加荷过程中组构各向异性滞后于应力各向异性有关。     

颗粒流强度折减法和重力增加法的边坡安全系数研究

将强度折减法和重力增加法的思路引入离散单元法,对土坡安全系数的评价作尝试性研究。运用颗粒流软件,采用强度折减法和重力增加法对边坡进行稳定分析。结果表明颗粒流计算得到的结果与有限元和条分法计算结果比较接近,为计算边坡稳定开辟了一条新的途径。采用颗粒流求解边坡的安全系数不需要条分,也即不需要对条间力做假定,同时不需要假定滑移面的位置和形状,颗粒根据所受到的接触力调整其位置,最终从抗剪强度最弱面发生剪切破坏。最后对强度折减法的不等比例折减提出了不同观点,为今后的研究提供思路。     

冻结砂土力学性质的离散元模拟

基于离散单元法颗粒流理论,土体颗粒单元间采用接触黏结模型中来考虑冻土中冰的胶结作用,建立了冻结砂土的颗粒流模型。通过改变计算模型中颗粒单元的参数,模拟了在不同冻结温度以及不同围压下冻结砂土的宏观力学性质,并与冻结砂土的室内试验结果进行了比较,结果表明：颗粒流方法可以较好地模拟冻结砂土的应力-应变关系以及剪切带的发展变化过程,颗粒流细观参数对温度具有显著的依赖性。研究结果对离散单元法在特殊土中的应用具有一定的理论和应用价值。     

船舶在平整冰区行进过程的离散元分析

船舶在冰区航行中,冰荷载是船舶结构设计和安全航行的重要影响因素。采用具有黏结破碎特征的离散单元法对船舶在平整冰区的航行过程进行了数值分析。在船舶破冰过程的离散元模拟中,平整冰由球形颗粒黏结而成,同时考虑了海水浮力和拖曳力对海冰单元的作用,颗粒间的黏结作用可在海冰与船体的相互作用过程中发生破碎;船体由三角形单元构造,分析了在海水浮力、推进器推力及海冰共同作用下船体的六自由度运动特性。通过船体与海冰相互作用的离散元分析,确定了船舶在破冰过程中海冰作用在船体上的冰荷载,探讨了船舶推进器推力和海冰厚度对冰荷载的影响,并与Lindqvist船体冰阻力经验公式进行了对比。     

基于DDEM的自然崩落采矿法崩落规律的数值模拟

利用DDEM（可变形离散单元法）原理和方法,研究了自然崩落法矿体崩落规律。针对具体地质条件和赋存情况及矿岩物理力学性质,采用计算机数值模拟分析了自然崩落法的崩落规律,得出了DDEM在用于模拟矿块自然崩落法中崩落规律和分析其崩落机理时是较适合的结论,不仅可以从力学角度对崩落的力学机理进行研究,而且能直观地用图形给出矿体在崩落过程中的具体形态变化,有效地指导放矿,现场试验验证效果良好。     

Microseismic Event Detection Kalman Filter: Derivation of the Noise Covariance Matrix and Automated First Break Determination for Accurate Source Location Estimation

Since 1972, Weir-Jones Engineering Consultants (WJEC) has been involved in the development and installation of microseismic monitoring systems for the mining, heavy construction and oil/gas industries. To be of practical value in an industrial environment, microseismic monitoring systems must produce information which is both reliable and timely. The most critical parameters obtained from a microseismic monitoring system are the real-time location and magnitude of the seismic events. Location and magnitude are derived using source location algorithms that typically utilize forward modeling and iterative optimal estimation techniques to determine the location of the global minimum of a predefined cost function in a three-dimensional solution space. Generally, this cost function is defined as the RMS difference between measured seismic time series information and synthetic measurements generated by assuming a velocity structure for the area under investigation (forward modeling). The seismic data typically used in the source location algorithm includes P- and S-wave arrival times, and raypath angles of incidence obtained from P-wave hodogram analysis and P-wave first break identification. In order to obtain accurate and timely source location estimates it is of paramount importance that the extraction of accurate P-wave and S-wave information from the recorded time series be automated—in this way consistent data can be made available with minimal delay. WJEC has invested considerable resources in the development of real-time digital filters to optimize extraction, and this paper outlines some of the enhancements made to existing Kalman Filter designs to facilitate the automation of P-wave first break identification.     

一种通用球面等积离散网格生成算法

球面离散网格模型是管理海量数据的有效途径,如何构建球面网格是这类模型的基础问题.文中利用Snyder等积多面体投影建立了平面与球面的对应关系,设计了球面等积三角形和四边形离散网格的生成算法.结合实验分析了两种网格的最大/最小边长比和单元周长均方差两项指标,结果表明全部网格的指标均收敛,适合用作多分辨率地理空间信息的管理.     

Simple and Ordinary Multigaussian Kriging for Estimating Recoverable Reserves

Multigaussian kriging is used in geostatistical applications to assess the recoverable reserves in ore deposits, or the probability for a contaminant to exceed a critical threshold. However, in general, the estimates have to be calculated by a numerical integration (Monte Carlo approach). In this paper, we propose analytical expressions to compute the multigaussian kriging estimator and its estimation variance, thanks to polynomial expansions. Three extensions are then considered, which are essential for mining and environmental applications: accounting for an unknown and locally varying mean (local stationarity), accounting for a block-support correction, and estimating spatial averages. All these extensions can be combined; they generalize several known techniques like ordinary lognormal kriging and uniform conditioning by a Gaussian value. An application of the concepts to a porphyry copper deposit shows that the proposed “ordinary multigaussian kriging” approach leads to more realistic estimates of the recoverable reserves than the conventional methods (disjunctive and simple multigaussian krigings), in particular in the nonmineralized undersampled areas.