移动荷载下瞬态振动的三维移动坐标有限元分析

针对移动荷载引起地基振动问题的特点,在随荷载移动坐标系下建立有限元方程进行分析,能够与常规有限元方法较好地符合。建立了可用于移动坐标有限元法的黏弹性边界,从而有效地模拟波在边界的透射。以Euler-Bernoulli梁模拟轨道和路基的作用,以三维实体单元模拟黏弹性分层地基,得到的数值模拟结果能够与瑞典Ledsgård地区高速铁路实测数据较好地吻合,表明移动坐标有限元法可以有效地分析此类问题。     

中国大陆地幔对流格局和岩石圈层构造运动

六十年代以来,由于卫星重力测量和计算技术的迅速发展,有些研究者,如Runcorn(1967)和Liu(1976-1980),根据卫星重力数据计算获得地幔对流及其应力场图象来研究全球板块构造和区域构造。愈来愈多的事实表明,构造运动不仅是地壳和岩石圈层物质运动的表现,而更重要的是地幔物质运动的反映。因此,本文应用作者计算获得的岩石圈层下面地幔流运动图象与大地构造及近代构造运动的资料相对比,探索我国岩石圈层下的地幔对流格局形成及其对构造运动的影响。     

A generic approach to modelling flexible confined boundary conditions in SPH and its application

In this paper, a new approach to applying confining stress to flexible boundaries in the smoothed particle hydrodynamics (SPH) method is developed to facilitate its applications in geomechanics. Unlike the conventional SPH methods that impose confining boundary conditions by creating extra boundary particles, the proposed approach makes use of kernel truncation properties of SPH approximations that occur naturally at free-surface boundaries. Therefore, it does not require extra boundary particles and, as a consequence, can be utilised to apply confining stresses onto any boundary with arbitrary geometry without the need for tracking the curvature change during the computation. This enables more complicated problems that involve moving confining boundaries, such as confining triaxial tests, to be simulated in SPH without difficulties. To further enhance SPH applications in elasto-plastic computations of geomaterials, a robust numerical procedure to implement Mohr-Coulomb plasticity model in SPH is presented for the first time to avoid difficulties associated with corner singularities in Mohr-Coulomb model. The proposed approach was first validated against two-dimensional finite element (FE) solutions for confining biaxial compression tests to demonstrate its predictive capability at small deformation range when FE solutions are still valid. It is then further extended to three-dimensional conditions and utilised to simulate triaxial compression experiments. Simulation results predicted by SPH show good agreement with experiments, FE solutions, and other numerical results available in the literature. This suggests that the proposed approach of imposing confining stress boundaries is promising and can handle complex problems that involve moving confining boundary conditions.     

河冰数值模拟中不规则边界处理技术探讨

数值模拟是河冰问题研究的重要手段之一, 但受限于天然河道的复杂地形条件, 研究工作开展中往往面临不规则边界的处理问题.基于河冰数值模拟工作现状, 总结和阐述了贴体网格技术对其的应用.针对二维数值模拟提出了采用滑动边界法和势流理论式相结合的方法来进一步提高边界处网格的正交性; 三维网格生成则宜采用σ竖向坐标变换法; 河冰模拟中水面位置确定通常采用标高函数法; 对于非稳态河冰数值模拟, 动边界处理技术可考虑改进并应用干湿法.     

A dynamic large deformation finite element method based on mesh regeneration

In this paper, a large deformation finite element (LDFE) approach termed ‘remeshing and interpolation technique with small strain (RITSS)’ is extended from static to dynamic soil-structure interaction applications. In addition, a technique termed ‘element addition’ is developed to improve the computational efficiency of both static and dynamic LDFE analyses that involve moving boundaries. The RITSS approach is based on frequent mesh generation to avoid element distortion. In dynamic RITSS, the field variables mapped from the old to the new mesh involve not only the stresses and material properties, but also the nodal velocities and accelerations. Using the element addition technique, new soil elements are attached to the domain boundaries periodically when the soil near the boundaries becomes affected by large displacements of the structure. The procedures of this Abaqus-based dynamic LDFE analysis and element addition technique are detailed, and the robustness of the techniques is validated and assessed through three example analyses: penetration of a flat footing into a half-space and movement of rigid and deformable landslides down slopes.     

Some garnet microstructures: an illustration of the potential of orientation maps and misorientation analysis in microstructural studies

The microstructures of two contrasting garnet grains are mapped using automated electron backscatter diffraction. In both cases there is a very strong crystallographic preferred orientation, with measurements clustered round a single dominant orientation. Each garnet grain is divided into domains with similar orientations, limited by boundaries with misorientations of 2° or more. In both samples most of misorientation angles measured across orientation domain boundaries are significantly lower than those measured between random pairs of orientation domains. One sample is a deformed garnet that shows considerable distortion within the domains. Lines of orientation measurements within domains and across domain boundaries show small circle dispersions around rational crystallographic axes. The domain boundaries are likely to be subgrain boundaries formed by dislocation creep and recovery. The second sample is a porphyroblast in which the domains have no internal distortion and the orientation domain boundaries have random misorientation axes. These boundaries probably formed by coalescence of originally separate garnets. We suggest that misorientations across these boundaries were reduced by physical relative rotations driven by boundary energy. The data illustrate the potential of orientation maps and misorientation analysis in microstructural studies of any crystalline material.     

无结构网格上二维浅水流动的数值模拟

提出在无结构网格上建立有限体积高性能格式族的统一框架,通过引入跨单元界面法向数值通量的逆风分解,将一维Osher、TVD两种通量分裂格式自然地推广至二维浅水方程组。给出了各种情况下浅水方程组有限体积法边界处理的计算公式。最后,利用该格式对穹包溢流、陆地动边界和河口潮流三个问题进行计算。结果表明,这类守恒格式具有高精度、无振荡性以及处理复杂流态过渡、自动捕俘间断和模拟陆地动边界的功能。     

钱塘江口涌潮的二维数值模拟

应用二维有限体积法、Osher格式及间断拟合法,计算了钱塘江河口涌潮产生、发展到消亡的全过程。该算法能保证水量和动量守恒,且能适应水深巨大变化和动边界。根据钱塘江涌潮的特点,可将其合理概化为一线潮。在网格尺寸较大情况下,计算的一线潮仍保持高分辨率,能准确满足水力学间断条件,且计算量小,可以在微机上实现。在对某实测半日小潮进行率定后,对随后的大潮进行了验证计算,涌潮的主要特征(如涌潮高度、移速、水位和流场等)与实测资料符合良好,证实了模型的合理性和模拟能力。     

天山乌鲁木齐河源1号冰川成冰带分布特征的再研究

通过分析2006年获取的43个雪坑和2004年的28个雪坑资料,对乌鲁木齐河源1号冰川成冰带的分布特征进行了再研究.结果表明: 1号冰川东、西支消融带和渗浸-冻结带界限(平均值)分别位于海拔4 068 m和4 086 m,渗浸-冻结带和渗浸带的界限分别位于海拔4 086 m和4 122 m处;东、西支顶部的渗浸-冻结带和渗浸带界限分别位于海拔4 184 m和4 435 m处.与20世纪80年代末相比,1号冰川的成冰带谱整体上移,成冰带界限达到各研究时期的最高,且东支的变化幅度大于西支;消融带和冰川顶部的渗浸-冻结带面积增大,渗浸带和冰川中部的渗浸-冻结带面积缩小.