一种激发极化法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.     

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

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

基于移动网格技术的热融湖动态演化过程数值模拟

基于移动网格技术建立了热融湖动态演化有限元数值模型, 研究了青藏高原多年冻土区典型热融湖动态演化过程, 分析了热融湖半径、深度的变化过程及其对湖底及周围多年冻土温度状况的影响. 结果表明:在移动边界热融湖模型中, 热融湖半径以0.7 m·a^-1的速度近线性地增大; 随着下伏高含冰量冻土的融化, 热融湖深度增加先慢后快, 最后逐渐减小趋于稳定. 热融湖深度和半径从5月末至翌年1月末增加显著, 在2-5月间基本保持稳定. 伴随着热融湖的扩展, 地表边界逐渐演变为湖底边界, 热融湖的热影响范围逐渐增大. 在固定边界热融湖模型中, 其热影响会逐渐趋于稳定, 由于初始尺寸大, 其湖底多年冻土退化速率大于移动边界模型, 而远离湖边的多年冻土退化速率要小于移动边界模型. 如果不考虑热融湖边界随时间的变化, 可能会高估热融湖对湖底多年冻土的热影响, 而低估其对附近多年冻土的热影响.     

基于非结构化网格的2.5-D直流电阻率自适应有限元数值模拟

2.5-D直流电阻率有限元数值模拟中,模型的剖分及加密主要通过手动实现.另外,采用的单元类型比较规则如矩形单元等,不易实现复杂模型的模拟.为解决上述问题,文中提出了一种自适应有限元算法.算法中采用稳健的后验误差估计来自动预测下一次网格的单元尺寸,直到设定的迭代条件满足为止.另外,采用非结构化三角形单元实现了任意复杂模型的灵活剖分.基于此,利用垂直接触面模型分析和对比了不同自适应策略的效率.通过对比发现,点源附近的单元得到了加密以消除源的奇异性.另外,对于任意一种策略,有限元结果均能最终收敛到精确解.最后,模拟了两个模型：2-D单个异常体模型和2-D地形模型.