基于四面体格网的3维复杂地质体重构

3维地质体的重构是3DGIS的一个重要研究方向,在地质与矿山等领域十分重要,但这一方面尚需大量细致深入的研究工作。在总结已有工作的基础上指出四面体格网是目前表达和重构地质体的最佳结构,同时提出了以简化的面约束条件重构各类复杂地质体的方法及表达具有属性渐变和属性突变特性的地质体的方法。     

河南省鲁山石墨矿浮选工艺研究

河南鲁山石墨矿属区域变质型矿床.根据矿石性质,进行了粗选条件试验,阶段磨矿阶段选别试验和阶段磨矿阶段选别预先分目试验,研究结果表明阶段磨矿是保护石墨大鳞片的有效方法,同时,结合预先分目能更好地提高石墨大片率.阶段磨矿阶段选别预先分目工艺闭路试验获得的指标为:当原矿品位为3.65%时,获得的技术指标为精矿品位C90.62...     

复杂区域三维有限元四面体网格自动生成研究

基于经典三维Delaunay四面体网格的生成算法,介绍了一种三维约束Delaunay三角化网格划分算法,通过对约束的恢复和局部退化的消除,保证了边界的一致性和网格的惟一性,实现了三维约束区域的四面体网格自动生成;通过对四面体的外接球半径与其最短边长度的比值（记为a）和四面体体积与其外接球体积的比值（记为n）的控制,避免质量不佳的单元,使生成的网格达到预期的质量要求,并提出了一种a和n的计算方法,简化了质量参数的计算过程;对几何形状复杂,特别是含有薄层的实体,可通过在实体上增加约束来提高网格的质量和缩短计算时间,使该方法能够实现岩土工程中不同地质条件下的复杂区域的三维有限元网格自动生成,并且能够保证生成的网格质量较好,为高精度的有限元计算提供了充分条件。计算实例表明,该算法能够将复杂三维实体剖分成质量较好的Delaunay四面体网格,可有效地应用于工程中有限元网格的生成。     

细网格推算流域面雨量方法应用浅析

介绍了流域细网格化的划分、在不固定测站条件下通过距离权重插补网格结点的雨 量和流域面雨量计算的方法，以及利用API函数(GetPixel)统计流域内网格结点的技术。     

基于非结构化四边形网格的VOF法

提出了一种基于非结构化四边形网格的VOF法,用以模拟二维自由表面流动问题。由于非结构化四边形网格在拟合任意复杂曲线边界上的优势,此法能够在不进行任何边界处理和不引入任何人工糙度的情况下直接、准确地求解具有任意边界形状的界面对流问题。此法采用PLIC法来重构界面,以保证界面具有二阶空间精度;采用MLER法来对流流体体积,以提高算法的体积守恒性。此外,为了提高算法的计算效率,文中针对非结构化四边形网格提出了一种确定界面截距的解析方法。大应变界面对流数值实验表明:此法具有较高的体积守恒精度和二阶的形状误差收敛率。     

河北省变网格数值预报系统设计和试验

设计了一种将原ＭＭ４模式的均匀网格变换成非均匀网格的方案,构造了一套适合河北省特点的变网格预报模式,同时建立了一套可投入业务使用的变网格数值预报系统,最后通过一次暴雨过程的模拟分析,证明变网格数值预报优于均匀网格数值预报     

An unstructured-grid, finite-volume, nonhydrostatic, parallel coastal ocean simulator

A finite-volume formulation is presented that solves the three-dimensional, nonhydrostatic Navier–Stokes equations with the Boussinesq approximation on an unstructured, staggered, z-level grid, with the goal of simulating nonhydrostatic processes in the coastal ocean with grid resolutions of tens of meters. In particular, the code has been developed to simulate the nonlinear, nonhydrostatic internal wave field in the littoral ocean. The method is based on the formulation developed by Casulli, in that the free-surface and vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly with the second-order Adams–Bashforth method, while the pressure-correction method is employed for the nonhydrostatic pressure in order to achieve overall second-order temporal accuracy. Advection of momentum is accomplished with an Eulerian discretization which conserves momentum in cells that do not contain the free surface, and scalar advection is discretized in a way that ensures consistency with continuity, thereby ensuring local and global mass conservation using a velocity field that conserves volume on a local and global basis. The nonhydrostatic pressure field is solved efficiently using a block-Jacobi preconditioner, and while stability is limited by the internal gravity wave speed and vertical advection of momentum, applications requiring relatively small time steps due to accuracy or stability constraints are run efficiently on parallel computers, since the present formulation is written entirely with the message-passing interface (MPI). The ParMETIS libraries are employed in order to achieve a load-balanced parallel partitioning that minimizes interprocessor communication, and the grid is reordered to optimize per-processor performance by limiting cache misses while accessing arrays in memory. Test cases demonstrate the ability of the code to efficiently and accurately compute the nonhydrostatic lock exchange and internal waves in idealized as well as real domains, and we evaluate the parallel efficiency of the code using up to 32 processors.     

利用网格划分软件和缓坡方程计算港内波浪场

提出了一套可方便应用于波浪场计算的方法。该方法包含了准确的边界输入输出技术,简便的地形输入技术,高质量的计算网格生成以及加速查寻、排序的数据结构。应用双曲型缓坡方程计算了港湾内的波浪场,数值计算时提出了处理建筑物表面部分反射边界条件的方法。将计算结果与边界元方法的结果和实验值进行了比较,证明该方法应用于波浪场可得到较理想的结果。对大连港区内的波浪折射绕射进行了计算,表明本方法可以应用于实际港口的波浪场计算中。     

A mesh deformation algorithm for free surface problems

Employing the simple iterative technique of adjusting the element positions using computed potentials to locate the free surface can lead to finite elements with large aspect ratios as the free surface drops towards the base of the mesh. In particular, free surface modelling of earth dams with base drains suffer from this problem. The paper suggests a number of steps which can be taken to alleviate mesh distortion problems and improve the numerical stability of the iterative finite element analysis. This leads to a mesh deformation algorithm which adjusts element widths in a simple fashion depending on the free surface height as the iterations proceed. The algorithm is specialized to the sloped earth dam problem, but may find application to other geometries. © 1997 by John Wiley & Sons, Ltd.     

加卸载作用下加筋垫层变形特征试验研究

加筋垫层由于具有诸多优点而被广泛应用于反复加卸载作用下公路、铁路路基等工程构筑物中。为研究加卸载作用下加筋垫层结构的变形特征,针对加筋与未加筋垫层2组模型,开展了单次加卸载静力试验。测试并获得了不同荷载作用下垫层竖向变形、筋材应变变化规律,对比分析了加筋与未加筋垫层的变形特征,并从能量角度对加筋垫层工作机制进行了探讨。结果表明：与未加筋情况相比,荷载作用下加筋垫层（加载板处）的沉降变形和残余变形更小。加筋后垫层表面（加载板范围外）的竖向变形及其受影响范围均增大,且在卸载过程中垫层表面竖向变形的水平分布曲线特征由“凹”型变为“凸”型,其加卸载曲线呈“∞”型。格宾网筋材应变沿横向呈现非均匀分布特点,其最大应变小于0.06%,筋材始终处于弹性变形状态。加筋垫层中筋材具有储存、释放及横向传递应变能效应,这使得加筋后垫层具有更好的承载能力与弹性性能,进而可降低循环加卸载作用下垫层的塑性变形或累计变形。