基于局部加密非结构化网格的三维电阻率法有限元数值模拟

总结了目前常用的结构化网格及其局限性,分析了非结构化网格对复杂地质体边界的适应性和Delaunay三角化算法及其网格加密策略,提出了一种新的局部节点加密方法,实现了对复杂模型的完全非结构化四面体全自动剖分,给出了三维直流电阻率模拟中四面体网格的质量评价标准和最优指标.计算和分析表明,数值解在点电流源及附近的奇异区精度最低,网格加密策略可以有效地减少其影响,极大地提高数值解的精度.本文提出的局部节点加密策略计算量最小,对精度的改善也优于经典方法.在精度要求苛刻或模型十分复杂时,局部体积加密策略和二次单元是高精度模拟的可靠保证.     

基于局部加密非结构网格的海洋可控源电磁法三维有限元正演

本文基于非结构网格实现了海洋可控源电磁法三维有限元正演模拟.该算法采用完全非结构网格剖分,可以模拟任意起伏地形和复杂地电模型.为了避免场源的奇异性,采用一次场/二次场分解算法,一次场由基于Schelkunoff势函数的一维解析公式得到.为了提高算法的精度和效率,采用对测点附近单元和异常体区域进行体积约束加密的方法,实现了非结构网格的局部加密.一、二维模型计算和分析表明,本文采用的局部加密方法能够明显地改善算法的精度,最大相对误差基本在1%以内.对三维模型计算及对比分析,说明了该算法对三维可控源电磁正演的实用性.复杂海底地形模型的正演模拟表明,海底地形对电磁场的影响很大,在进行海洋可控源电磁资料解释时,地形的影响有必要考虑在内.     

一种激发极化法2.5维正演的自适应有限元方法（英文）

传统的基于结构化网格有限元法采用的单元比较规则如矩形等,且网格剖分和加密要靠手动实现,所以传统的基于结构化网格有限元法不能准确和灵活地模拟复杂介质。本文采用易于模拟复杂介质模型的非结构化三角形网格进行剖分,且利用对偶加权后验误差估计指导网格自动细化过程,然后在电位模拟的基础上计算雅可比偏导矩阵,并依据Seigel（1959）理论实现激发极化法2.5维自适应有限元正演模拟算法。通过对垂直接触面模型进行正演分析,接收点附近网格得到了明显加密,电位数值解平均相对误差收敛到0.4%,视极化率平均相对误差收敛到1.2%,表明经自适应网格细化后,该算法数值解最终能收敛到精确解附近。最后对两个较复杂模型进行了正演计算与分析,进一步验证了该算法的准确性和灵活性。     

土钉墙在深基坑中的设计应用

介绍了滕州贵和世纪佳苑一号楼及公寓楼地下车库基坑支护方案的设计和施工,采用了土钉、挂网锚喷等工艺方法进行了支护设计及施工,做到了经济合理和安全可行,并介绍了本工程遇到的特殊情况的处理,给类似工程有很好的借鉴作用。     

基于轮廓线三维矿体表面重建的一种改进算法

基于轮廓线的三维矿体表面建模方法是矿体建模的主流方法,但在实际应用中,传统的建模方法显现出了不足之处,笔者对该方法进行了针对性地改进。对于一组单轮廓线可自动添加矿体趋势线,并且还可以进行人工编辑修改,然后利用中间加密轮廓线的方法实现对矿体形态的控制。通过投影计算封闭轮廓线之间的最短距离自动添加分支点,利用平面的带洞限定三角剖分实现分支的自动构建,大大节省了人力资源,同时保证了分支矿体的准确性。针对初始构建的三维矿体表面模型几何质量差,引入了质量控制,实现了表面模型的重构,保证了模型质量和后续的计算。这些改进在实际的应用中取得了良好的效果。     

钻井船插桩CEL数值模拟中的若干问题分析

依据模拟钻井船在不同场地条件下贯入阻力的离心模型试验结果及具体工程实例,对利用耦合欧拉-拉格朗日（CEL）有限元方法模拟钻井船在黏土、砂土、黏土下覆砂土、砂土下覆黏土及成层土场地插桩过程时,影响贯入阻力计算结果的几个因素进行了研究。结果表明,对于不同土层场地条件,有限元模型中欧拉区域范围对贯入阻力几乎没有影响。为确保CEL数值结果的精度,有限元建模时靠近桩靴部分设置为细网格区域,以外区域设置为粗网格区域;对于不同土层场地条件,减小细网格尺寸及增大细网格范围可以减小贯入阻力的振荡情况;经比较总结,细网格尺寸建议取0.05倍桩靴直径,细网格范围建议取2倍桩靴直径。采用位移控制模拟钻井船插桩时,桩靴贯入速率对黏土场地的贯入阻力影响较小,对砂土下覆黏土场地的贯入阻力影响很大,对一般成层土场地的贯入阻力有一定影响。经比较总结,建议在研究钻井船在不同土层场地的计算贯入阻力时,有限元模型的桩靴贯入速率取0.2 m/s。     

一种大规模倾斜摄影模型三维可视化方案

提出了一种基于四叉树结构的大规模倾斜摄影模型三维可视化方案：首先采用二次误差测度（QEM）的半边折叠算法对三角网格模型进行分层几何简化;然后对分层简化的模型进行纹理映射和分块,分块结果采用Morton编码命名;最后将场景按金字塔结构进行LOD（多细节层次结构）组织。选取某一试验区域的真实倾斜模型数据进行了试验,结果表明,该方案能依据模型在计算机屏幕上所占的像素大小,动态加载倾斜模型数据,完成大规模场景的实时可视化。     

Research on Bulbous Bow Optimization Based on the Improved PSO Algorithm

In order to reduce the total resistance of a hull, an optimization framework for the bulbous bow optimization was presented. The total resistance in calm water was selected as the objective function, and the overset mesh technique was used for mesh generation. RANS method was used to calculate the total resistance of the hull. In order to improve the efficiency and smoothness of the geometric reconstruction, the arbitrary shape deformation (ASD) technique was introduced to change the shape of the bulbous bow. To improve the global search ability of the particle swarm optimization (PSO) algorithm, an improved particle swarm optimization (IPSO) algorithm was proposed to set up the optimization model. After a series of optimization analyses, the optimal hull form was found. It can be concluded that the simulation based design framework built in this paper is a promising method for bulbous bow optimization.     

Impact of water and nitrogen fracturing fluids on fracturing initiation pressure and flow pattern in anisotropic shale reservoirs

A numerical model is proposed to investigate the impact of water and nitrogen fracturing fluids on the fracturing initiation pressure and the flow pattern in anisotropic shale reservoirs. This model considers the anisotropy of shale deformation and permeability, the compressibility of fracturing fluid, and the fluid moving front. A crack initiation criterion is established with the stress intensity factor of mode I crack. Both crack initiation pressure and seepage area are verified and analyzed. These results show that shale deformation and permeability anisotropy, fracturing fluid compressibility, viscosity, and pressurization rate have significant impacts on fracturing initiation pressure and seepage area.     

有自由面渗流分析的三维数值流形方法

提出了求解有自由面渗流问题的三维数值流形方法,通过构造任意形状流形单元的水头函数,推导了流形单元的渗透矩阵和无压渗流分析的总体控制方程,并给出了自由面的迭代求解策略和渗透体积力的计算方法。典型算例的数值分析表明,该方法采用数学网格覆盖整个材料区域,在自由面的迭代求解过程中数学网格保持不变,只考虑自由面以下渗流区的介质,只对自由面以下的流形单元形成总体渗透矩阵,具有精度高、收敛速度快、编程简单等优点,而且能够通过单纯形积分精确计算被自由面穿越单元的渗透作用力,因此,特别适用于有自由面渗流问题的模拟。