改进的HASM-AD算法及在空间变量模拟的应用分析

高精度曲面建模(HASM)可以显著提高空间曲面模拟的精度,但是计算速度低限制了该模型的进一步应用。为了提高HASM模型的计算效率,本文对HASM-AD算法作了改进,通过在计算过程中为采样点添加索引,避免了计算过程中对采样点信息的重复查找操作;同时,在遍历独立计算单元时实时计算第一类、第二类基本量及克式符号,避免了全局存储上述变量所需要的额外内存消耗。数值试验表明,由于将全局线性方程组求解问题转化为局部独立计算单元(5×5栅格)内的方程组求解,改进的HASM-AD算法显著提高了计算效率,同时降低了模型运行过程中的内存消耗。最后,本文以全国陆地降水空间分布模拟分析作为实例,验证改进的HASM-AD算法模拟精度及计算效率,模拟结果表明,改进的HASM-AD算法模拟结果精度优于其他HASM算法(以HASM-PCG为例),并且计算效率优势更为明显,实现全国10km分辨率的降水分布模拟耗时仅为4s。表明改进的HASM-AD算法提高了计算速度,并且适于大尺度的空间变量模拟应用。     

变面积差动电容式传感器的研究

介绍一款变面积式差动电容传感器.该传感器借助于两路相位相差180°的激励信号源作用于两固定极,使动摆片与定摆片之间形成电场.动摆片在进行摆动的过程中,在左右两固定摆中的面积会发生相应的变化,进而会产生电容量的改变.动摆片摆动角度的大小与动摆片在定片中的面积以及电容量的大小呈线性比例关系.通过交流电桥法可测量出该电容传感器的电容量的变化所对应的电压输出,并推算出动摆片位移角度的大小.实验验证了该设计的可行性.     

基于循环注蒸气和烟气激励的海上稠油开采技术

通过实验和油藏数值模拟,研究循环注蒸气和烟气激励的加热降黏、气体溶解降黏、气体扩大加热腔与减少热损失、蒸气与气体协同增产效应等开采机理,优化注入蒸气和烟气温度、烟气蒸气体积比、周期蒸气和气体注入量等参数,建议注入蒸气与烟气的温度高于300℃,在大气压下烟气蒸气体积比为0.4～0.6.在渤海稠油油田开展4口井循环注蒸气和烟气激励开采的现场试验.结果表明,循环注蒸气和烟气激励开采后稠油井平均日产油量达到冷采的3倍以上,最高日产油量超过100m3,周期产油量达到1.0×104 m3左右.基于循环注蒸气和烟气激励的海上开采稠油技术可以提高海上稠油开采速度和采收率,有潜力成为未来海上稠油油田的主要热采开发方式.     

基于场变量的边坡稳定分析有限元强度折减法?

有限元强度折减法在实施过程中往往要根据不同的折减系数手工修改输入文件,并不断进行试算,过程较繁琐。为了简化计算过程,提高计算效率,在有限元强度折减法的基础上,结合ABAQUS提供的场变量,提出一种适用于有限元强度折减计算的新方法--基于场变量的有限元强度折减法。通过在ABAQUS中设置岩土体黏聚力与内摩擦角为场变量的函数,场变量为增量步时间t的函数,从而通过控制增量步时间t实现黏聚力与内摩擦角的折减,通过一次计算即可求得边坡的安全系数,无需手工修改折减系数反复试算。结合有关文献对该方法用于边坡稳定分析的可行性进行了验证,并与二分法确定安全系数的计算效率进行了比较。结果表明：该方法求得的安全系数是有效可靠的,用于边坡稳定性的计算是可行的,简化了有限元强度折减法的计算过程,提高了计算效率。     

泥质砂岩蠕变特性与锚固控制效应试验研究

通过泥质砂岩单轴压缩蠕变试验,揭示岩石蠕变变形与破坏的特点,反映了不同应力路径下的岩石蠕变速率与发展过程,建立了泥质砂岩的非线性蠕变模型,并通过对试件施加侧向刚性约束和锚固约束,以反映地下工程围岩在平面应变状态和支护状态下的蠕变特性。研究结果表明,锚固约束可以显著提高岩石发生稳定蠕变的应力阀值,减缓蠕变速率。根据改进的蠕变方程进行试验数据的拟合,确定了改进的西原模型蠕变方程中的各物理量。与单轴蠕变试验结果对比表明,该模型的计算结果与试验数据吻合较好,说明改进的西原模型可较好地反映岩石加锚后的蠕变特征。研究结果对深部地下工程围岩的变形与稳定控制研究具有重要的指导意义。     

非饱和渗流问题的自适应欠松弛变量变换方法

在使用有限元方法求解非饱和土渗流问题时,土-水特征曲线和渗透率函数的强烈非线性经常会造成计算中出现迭代不收敛、计算误差大等问题。基于变量变换的思想,结合时间步长自适应技术提出了一种求解非饱和渗流问题的新方法--欠松弛RFT变换方法（ATUR1）。ATUR1方法通过变量变换,大大降低了Richards方程中未知数在空间和时间上的非线性程度,从而改善这种非线性所带来的计算收敛困难和精度差等问题。欠松弛技术的引入减少了迭代过程中的振荡现象,进一步提高了非线性迭代计算的效率。时间步长自适应技术则有效地控制整个计算过程的误差。数值算例结果说明,ATUR1可以有效地提高计算效率和精度,是一种准确有效的计算方法。     

Retrieving earthquake signature in grace gravity solutions

The GRACE satellites have been orbiting the Earth since 2002, monitoring the time variable gravity field. Some of the observed fluctuations are due to geodynamic causes, but they are often hidden in the complex signal, composed of hydrology, ocean, atmosphere, and geodynamics, the signal of geodynamic origin being usually the smallest. In addition, dealiasing residuals and noise make the separation of the signal from the different causes more difficult. We proposed a method based on the Empirical Orthogonal Function decomposition to extract the signal of physical origin, under the hypothesis that the physical signal is spatially more consistent than the noise and aliasing incomplete correction. We used synthetic geoid variations associated with earthquakes located at nearly 2000 positions at the Earth surface, based on several examples of large actual subduction events. We show that, with the present day accuracy, we can retrieve the geoid variations associated with more than 98 per cent of the earthquakes of magnitude 9 or above, around 60 per cent for magnitude 8.8, 40 per cent for magnitude 8.6 and 33 per cent for magnitude 8.3. Some events, with the right properties and location, can be detected with magnitude as low as 8. We then applied the method to the GRACE solutions, and retrieved the Hokkaido event (2003) and the Sumatra event (2004), which is in agreement with the retrieval rates mentioned here above.     

砂土地基的黏性特性及其有限元模拟

采用非线性弹塑黏性有限元法模拟了砂土地基变速率加载室内试验,再现了砂土地基荷载-沉降曲线对加载速率突变、砂土蠕变以及应力松弛等加载情况下的响应和瞬时黏性中砂土黏性随加载进行的衰减现象。在有限元分析中,以非线性三要素模型为理论框架构建了可以综合考虑砂土黏性特性的弹塑黏性本构模型,用动态松弛法进行有限元求解,最后通过有限元与室内试验的对比分析研究了砂土地基所具有的黏性特性,同时验证了非线性三要素弹塑黏性模型用于模拟砂土黏性特征的合理性和正确性。     

A comprehensive study on the parameters setting in smoothed particle hydrodynamics (SPH) method applied to hydrodynamics problems

Smoothed particle hydrodynamics (SPH) is a meshfree, Lagrangian particle method which has advantages in handling solids with extremely large deformation. Like any other numerical methods, cares must be taken to ensure its desirable accuracy and stability through considering several correction techniques in calculation. The selection of values for parameters in those correction approaches is a key step in SPH simulation, which is always difficult for new beginners to deal well with effectively. This paper examines the common inconsistency and instability problems in SPH method and studies its computational efficiency when applied to hydrodynamics problems with material strength like soil column collapse. We analyzed in detail how the correction techniques mitigate these inconsistency and instability problems. Also, the numerical testing results associate with different values for the parameters used in the correction techniques are provided for better understanding the influence of these parameters and for finding out the desirable values. It is found that (1) the SPH method is easily subjected to an inconsistency problem in the boundary area due to the boundary deficiency, and it can be treated well by adopting “virtual particles” contributing to the particle summations. (2) The numerical oscillation in SPH simulation can be mitigated effectively by artificial viscosity with the suggested parameter values. (3) The tension cracking treatment, artificial viscosity and artificial stress work well in removing the tensile instability problem in SPH method. In addition, the nearest neighboring particle searching (NNPS) algorithm, spacing ratio, smoothing length and time step influence the efficiency and accuracy of SPH method significantly. It is shown that SPH method with suggested parameters values can produce a very good result compared with the experimental result.