A finite volume approach with local adaptation scheme for the simulation of free surface flow in porous media

We present a method for solving steady‐state flow with a free surface in porous media. This method is based on a finite volume approach and is halfway between a fixed and an adaptive mesh method, taking advantage of both approaches: computational efficiency and localization accuracy. Most of the mesh remains fixed during the iterative process, while the cells in contact with the free surface (free surface cells) are being reshaped. Based on this idea, we developed two methods. In the first one, only the volumes of the free surface cells are adapted. In the second one, the computational nodes of the free surface cells are relocated exactly at the free surface. Both adaptations are designed for a better application of the free surface boundary conditions. Implementation details are given on a regular finite volume mesh for the case of homogeneous and heterogeneous rectangular dams in 2D and 3D. Accuracy and convergence properties of the proposed approach are demonstrated by comparison with an analytical solution and with existing references. Copyright © 2011 John Wiley & Sons, Ltd.     

云南丽江台OSG-066重力潮汐观测结果及应用

基于云南丽江台超导重力仪OSG-066超过5 a的连续重力潮汐观测资料，获取高精度潮汐参数并构建地方重力潮汐模型，为丽江地区精确重力潮汐改正提供参考。研究9个全球海潮模型在丽江台的重力影响改正有效性，用时域内回归分析法和频域内小波分析法分别研究大气重力导纳值的变化特征，并利用观测潮汐数据估算自由核章动本征参数。结果表明，9个海潮模型在丽江台的重力影响改正效果基本一致，NAO.99b模型的改正精度较高；大气重力导纳值在时域和频域内均显示出较为明显的季节性变化特征，且与降雨量变化存在负相关关系。自由核章动本征周期为431.5(411.2, 454.0)恒星日，品质因子为 -2 796。     

淋溶时间对红黏土物理力学特性的影响

以湖北省鄂州市庙岭镇的红黏土为研究对象,以化学选择溶解法结合淋溶作用为处理手段,将原状土样中的游离氧化铁按梯度去除,获得连二亚硫酸钠-柠檬酸钠-碳酸氢钠(DCB)溶液淋溶时间与除铁率的关系,并分析了不同淋溶时间对红黏土物理力学特性及微观孔隙结构的影响。试验结果表明,DCB溶液淋溶时间与除铁率相关性较强,除铁率与时间呈先快速增长后逐渐稳定的关系。由于游离氧化铁在土颗粒间主要起胶结作用,故淋溶时间对红黏土的物理力学特性有较大影响。随着淋溶时间的增长,黏粒与胶粒含量逐渐增大,耐热性小幅提升,无侧限抗压强度显著减小,且无侧限抗压强度变化整体趋势为前期急剧降低至后期趋于平稳。核磁共振、差热分析和电镜扫描的结果表明,随着淋溶时间的增大,庙岭红黏土内部孔隙增大,自由水减少,结合水增多,微观形态上团粒结构遭到破坏,胶结物质明显减少,结构形式从紧凑的粒团堆叠结构逐渐转化为聚集体-松散颗粒状结构。     

利用体应变观测资料检测地球自由振荡

选用通河地震台体应变观测数据,采用功率谱密度估计方法,获得2015年5月30日日本小笠原群岛地区8.0级大地震和2017年9月8日墨西哥沿岸近海8.2级大地震激发的_0S_4—_0S_(35)基频球型自由振荡频率,并与地球初步参考模型(PREM)的理论自由振荡频率进行对比,结果发现:前者共被检测到11个基频信号,有7个信号与理论值存在偏差,1个存在频谱分裂现象;后者共被检测到17个基频信号,有9个信号与理论值存在偏差,1个存在频谱分裂现象。     

激光跟踪仪用于大型预埋板形位误差的检测

根据激光跟踪测量系统的特点,结合工程实例,用跟踪仪自由设站测得预埋板所要检测的数据;用平面拟合、直线拟合、坐标转换等数据处理方法,得出偏差数据,很好地对预埋板的形位误差进行了检测,以便设计、施工的需要。     

A new parabolic variational inequality formulation of Signorini's condition for non‐steady seepage problems with complex seepage control systems

By extending Darcy's law to the dry domain above the free surface and specifying the boundary condition on the potential seepage surfaces as Signorini's type, a partial differential equation (PDE) defined in the entire domain of interest is formulated for non‐steady seepage flow problems with free surfaces. A new parabolic variational inequality (PVI) formulation equivalent to the PDE formulation is then proposed, in which the flux part of the complementary condition of Signorini's type in the PDE formulation is transformed into natural boundary condition. Consequently, the singularity at the seepage points is eliminated and the difficulty in selecting the trial functions is significantly reduced. By introducing an adaptive penalized Heaviside function in the finite element analysis, the numerical stability of the discrete PVI formulation is well guaranteed. The proposed approach is validated by the existing laboratory tests with sudden rise and dropdown of water heads, and then applied to capture the non‐steady seepage flow behaviors in a homogeneous rectangular dam with five drainage tunnels during a linear dropdown of upstream water head. The non‐steady seepage flow in the surrounding rocks of the underground powerhouse in the Shuibuya Hydropower Project is further modeled, in which a complex seepage control system is involved. Comparisons with the in situ monitoring data show that the calculation results well illustrate the non‐steady seepage flow process during impounding and the operation of the reservoir as well as the seepage control effects of the drainage hole arrays and drainage tunnels. Copyright © 2010 John Wiley & Sons, Ltd.     

B‐bar based algorithm applied to meshfree numerical schemes to solve unconfined seepage problems through porous media

The object of this work is to establish a meshfree framework for solving coupled, steady and transient problems for unconfined seepage through porous media. The Biot's equations are formulated in displacements (or u − w) assuming an elastic solid skeleton. The free surface location and its evolution in time are obtained by interpolation of pore water pressures throughout the domain. Shape functions based on the principle of local maximum entropy are chosen for the meshfree approximation schemes. In order to avoid the locking involved in the fluid phase of the porous media, a B‐bar based algorithm is devised to compute the average volumetric strain in a patch composed of various integration points. The efficiency of such an implementation for one phase problems is shown through the Benchmark problem, Cook's membrane loaded by a distributive shear load. The proposed methodology is firstly applied to various classical examples in unconfined steady seepage problems through earth dams, then to the dynamic consolidation of a soil column. The results obtained for both problems are quite satisfactory and demonstrate the feasibility of the proposed method in solving coupled problems in porous media. Copyright © 2015 John Wiley & Sons, Ltd.     

计算机辅助浮雕设计与加工系统

介绍一个具有一定特色的ＣＡＤ／ＣＡＭ系统──ＤＭＲ系统。它采用模块化程序设计思想，各模块相互独立，用数据文件传递信息。在创成型面时有机地结合ＰＢ样条曲线理论和Ｂｅｚｉｅｒ曲面理论，简化了其中的一些计算；在计算刀具轨迹时又采用简单快速的Ｂｅｚｉｅｒ曲面离散法，计算速度快、精度高，能在微机上运行。     

用主成分法求解秩亏自由网的研究

结合多数平差和主成分估计理论,导出误差方程中含多重共线性关系时求解未知数的公式,并证明主成分估计的解是最小二乘最小范效解。由此,将主成分估计推广到秩亏自由网平差中,导出求解未知多数及其协因数的公式,并据主成分估计演绎出用附有条件的参数平差法求解秩亏自由网的条件方程（基准在达式）。     

Asymptotic eigenfrequencies of the Earth's normal modes

We derive asymptotic formulae for the toroidal and spheroidal eigenfrequencies of a SNREI earth model with two discontinuities, by considering the constructive interference of propagating SH and P-SV body waves. For a model with a smooth solid inner core, fluid outer core and mantle, there are four SH and 10 P-SV ray parameters regimes, each of which must be examined separately. The asymptotic eigenfrequency equations in each of these regimes depend only on the intercept times of the propagating wave types and the reflection and transmission coefficients of the waves at the free surface and the two discontinuities. If the classical geometrical plane-wave reflection and transmission coefficients are used, the final eigenfrequency equations are all real. In general, the asymptotic eigenfrequencies agree extremely well with the exact numerical eigenfrequencies; to illustrate this, we present comparisons for a crustless version of earth model 1066A.