软土蠕变固结特性研究

通过大量对比试验系统研究了广东典型软土的蠕变固结特性和计算模型。试验结果表明：软土的固结和次固结在变形中所占比例和发生受多种因素影响,包括初始荷载、加荷比、土体固结度、含水率、孔隙比;次固结系数与土体的应力历史、当前固结状态及加荷比有关,在预压荷载下会明显减小;次固结系数与压缩指数之比一般为常数,其值在0.025?0.10之间;剪切蠕变特性与土体所承受的固结压力有关,土样蠕变规律呈应变硬化特征;三轴蠕变特性与排水条件关系密切,固结作用可以削弱蠕变效应;同等加荷条件下,偏应力-轴向应变关系在排水剪切中呈线性关系。而在不排水剪切中呈非线性,并有明显的屈服特征;不排水剪切条件下持续加荷的加荷方式可以适当提高不排水强度。根据上述试验成果分别从考虑非线性本构关系和考虑排水条件分析了蠕变固结模型的建立路径。     

自适应自然单元法研究——自适应细化

对h型自适应自然单元法的自适应细化方案进行了初步研究。在ZZ误差分析的基础上实现了节点的自动加密,使得随着自适应细化的进行,求解误差减小,而且误差分布趋于均匀。在数值分析中,主要有两种误差来源--插值误差和积分误差。随着节点的加密,Delaunay三角形的尺寸随之减小,三角形内的应力场趋于线性分布,那么插值误差和积分误差也都会随之减小。因而,h型自适应分析可以同时减小上述两种误差而达到不断提高求解精度的目的。由于自然单元法求解依赖于求解域内离散节点的Voronoi结构,建议的细化方案中新节点的引入只需局部调整Delaunay结构,算法的实现极为容易,程序实现简单、高效。研究表明,建议的自适应方案是可行的,自然单元法特别适合进行h型自适应分析。     

一种真空预压区地下水位计算方法

根据抽真空时加固区内真空渗流场的形成机制和变化特点,把真空渗流场的形成过程分为非稳定渗流和稳定渗流两个阶段。利用渗流量计算公式,计算加固区内渗流量随时间变化过程曲线。结合两个阶段加固区地下渗流场的变化特性,分别利用Theis和Dupuit计算模型计算不同时刻的地下水位,并用实测进行验证。     

液化土层上建筑物倾斜计算要点初步分析

地震中饱和砂土地基液化会引起结构物倾斜并导致其功能丧失,但目前缺乏相应的数值模拟方法。通过振动台简单模型试验,寻找输入波-基底竖向动应力-基底孔压-结构震陷之间的关系,提出发展液化土层上建筑物倾斜数值模拟方法所需要考虑的要点和应满足的条件。结果表明：（1）分析方法中孔压增长模型应适于不规则波计算,能准确地计算出孔压增长过程,准确地计算出峰值一样但不同波形下孔压增长的差别;（2）孔压增长模型应能反映土的各向异性特性对孔压增长的影响,能合理地计算出拉、压不同应力作用下孔压增长的差异;（3）孔压增长模型应能反映非均等固结条件对孔压增长过程的影响,能合理地计算出结构底部不同固结比土体中孔压增长过程;（4）分析方法中土体变形的计算应能跟踪液化过程中的变形发展,且具备大变形计算能力。     

Backward Calculation Based on the Advection and Diffusion of Oil Spills on the Sea Surface

Advection and diffusion equation is an effective means of investigating the motional state of oil spills. Common study is based on the known initial and boundary conditions to calculate the equation onwards. Unlike previous works, in this paper we focus on backward calculating advection and diffusion equation, namely, starting computation from the moment when the oil spill is found back to the moment when the incident occurs. Through such a simulation, we hope to trace the track of oil spills …     

射弹在岩石介质中的侵彻实用计算方法

根据岩石（混凝土）类材料在射弹侵彻时所产生的变形和破坏现象,考虑靶体材料的可压缩性,建立靶体材料变形的数学模型,得出射弹在靶体中侵彻时所受的阻力,再根据牛顿运动定律建立射弹的运动方程,用增量方法求解方程,得到弹体侵彻过程中的参数及侵彻深度。与经验公式相比较,说明该方法具有实用性和可靠性。     

三维地形大地电磁场的边界元模拟方法

提出了一种用边界元法计算大地电磁场三维地形影响的数值模拟方法.首先用矢量积分理论和电磁场边界条件, 将上半空间(空气)和下半空间(地下介质)两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程, 其中一个计算磁场, 称磁场方程; 另一个计算电场, 称电场方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中, 假设单元中电磁场为水平均匀大地空间电磁场与地形影响的迭加, 并假设地形影响为常项, 这样既保证了计算精度又使得计算方法简便.通过分解和计算, 每一个矢量面积分方程分解为对应3个坐标方向的3个常量线性方程, 这些线性方程组成了对角占优的线性方程组, 可用SSOR方法求解.文中给出了2个三维地形上大地电磁视电阻率曲线的计算结果.      

外壁恒温条件下冻结管壁热流密度变化规律数值计算研究

应用相似理论将冻结温度场方程和参数无量纲化,然后建立无量纲的数值计算模型,研究了在外表面恒温情况下冻结管壁的无量纲热流密度与各无量纲影响参数间的关系.结果表明:无量纲热流密度随未冻土与冻土的导热系数之比线性增加,与冻土与未冻土的容积比热之比成线性关系,随无量纲结冰潜热线性增长,与无量纲盐水温度成线性关系;可表示为无量纲时间的负三分之一次方的线性函数.最后对167706次数值计算结果进行了回归分析,得到具有较高精度的回归计算公式.     

岩石爆破的相对损伤与损伤累积计算

为了便于工程应用,利用损伤变量定义,确定既真实反映力学作用过程,又易于获取和进行客观度量损伤变量的一种方法。通过提出的相对损伤概念,用相对损伤变量表示含损伤材料在载荷作用下损伤扩展程度和不同部位的相对损伤程度,使衡量受力作用程度的损伤变量确定更易实现。同时,提出了一种基于裂纹扩展的损伤累积计算方法,并用实例计算说明了其合理可行。     

三维地层的三棱柱剖分与土方计算

三维地层剖分算法是实现地层三维可视化与土方计算的重要基础。采用三棱柱体建模方法构造三维地层,分析任意平面切割三棱柱的不同情形,得出了针对三棱柱地层的剖分算法;提出了以切割面为分界面,将单个三棱柱剖分成多个小三棱柱的三角投影分割法--将三棱柱顶面、底面及切割截面按一定规则进行合理三角化划分,使划分后三角形在水平面上的投影相重合;将该3层三角形组合成以切割面为分界面的2层小三棱柱,从而完成剖分。土方计算中首先给出开挖体的边界控制方程,根据该控制方程所限定的内边界得到开挖区域,结合上述剖分算法,对该区域内所有小三棱柱体积求和,得到开挖土方量。将以上算法应用到广东省某重大水利工程中,实现了基坑模拟开挖及其土方计算。