盾构壁后注浆对地表沉降影响模拟研究

结合宁波轨道交通某区间隧道施工, 研究了盾构壁后注浆对地表沉降的影响, 考虑壁后注浆压力消散和浆液固结硬化过程, 采用FLAC3D有限元软件对开挖过程模拟计算, 分析了不同注浆量和不同注浆压力下的地表沉降变化规律以及浆液硬化作用对地表沉降和管片位移的影响规律, 并将地表沉降和拱顶位移的计算值和实测值进行对比分析. 结果表明: 本工程最佳注浆量4.0~5.0m3, 最佳注浆压力约0.3MPa; 硬化作用引起的地表沉降占总沉降约5%, 引起的管片位移占总位移约15%. 针对宁波土质, 合理确定注浆量和注浆压力能够有效控制地表沉降, 合理配制浆液, 实时监控注浆量、注浆压力和浆液强度, 能够有效控制管片的上浮和变形.     

A New Higher Order Shear Deformation Model for Static Behavior of Functionally Graded Plates

In this paper, a new displacement based high-order shear deformation theory is introduced for the static response of functionally graded plate. Unlike any other theory, the number of unknown functions involved is only four, as against five in case of other shear deformation theories. The theory presented is variationally consistent, has strong similarity with classical plate theory in many aspects, does not require shear correction factor, and gives rise to transverse shear stress variation such that the transverse shear stresses vary parabolically across the thickness satisfying shear stress free surface conditions. The mechanical properties of the plate are assumed to vary continuously in the thickness direction by a simple power-law distribution in terms of the volume fractions of the constituents. Numerical illustrations concerned flexural behavior of FG plates with Metal-Ceramic composition. Parametric studies are performed for varying ceramic volume fraction, volume fraction profiles, aspect ratios and length to thickness ratios. The validity of the present theory is investigated by comparing some of the present results with those of the classical, the first-order and the other higher-order theories. It can be concluded that the proposed theory is accurate and simple in solving the static behavior of functionally graded plates.     

Flow Characteristics of Flapping Motion of a Plane Water Jet Impinging onto Free Surface

A submerged turbulent plane jet in shallow water impinging vertically onto the free surface will produce a large-scale flapping motion when the jet exit velocity is larger than a critical one. The flapping phenomenon is verified in this paper through a large eddy simulation where the free surface is modeled by volume of fluid approach. The quantitative results for flapping jet are found to be in good agreement with available experimental data in terms of mean velocity, flapping-induced velocity and turbulence intensity. Results show that the flapping motion is a new flow pattern with characteristic flapping frequency for submerged turbulent plane jets, the mean centerline velocity decay is considerably faster than that of the stable impinging jet without flapping motion, and the flapping-induced velocities are as important as the turbulent fluctuations.     

MULTI-OBJECTIVE OPTIMIZATION AND AERODYNAMIC PERFORMANCE ANALYSIS OF THE UPPER SURFACE FOR HYPERSONIC VEHICLES

为分析小攻角巡航条件下吸气式高超声速飞行器上壁面的变化对其气动性能和容积的影响, 以参数化后的飞行器上壁面对称面型线为设计变量, 在飞行马赫数6.5, 飞行高度27 km, 飞行攻角为4°的条件下, 采用计算流体力学为性能分析工具, Pareto多目标遗传算法为优化设计方法, 开展了二维条件下的升阻比/容积双目标优化设计. 在此基础上, 选择典型的二维优化结果, 重构生成对应的三维构型并进行数值分析, 获得了飞行器气动性能和容积间的相互关系. 结果表明在巡航条件下, 尽管二维/三维条件下飞行器的气动参数数值有较大差别, 但在这2种条件下, 飞行器的升阻比和容积间的关系均近似呈线性反比例关系. 同时, 对于三维构型而言, 在给定容积不变的条件下, 通过改变上壁面对称面型线的形状仅能使升阻比获得较小的增量(约0.36%). 相比之下, 当给定升阻比基本不变的条件下, 飞行器容积可调空间相对较大, 约为1.93%. 此外, 计算结果还表明, 在飞行器的容积基本不变情况下, 通过调节上壁面对称面型线, 可使飞行器的俯仰力矩获得5%左右的调节空间, 且其升阻比基本不变.     

Assessment of existing and introduction of a new and robust efficient definition of the representative volume element

Accurate numerical homogenization necessitates the thorough determination of the Representative Volume Element (RVE). There exists several seminal works on the notion of the RVE in homogenization, its definitions and methods of determination for efficient computation of composite effective properties. The objective of the current work is to assess the ability of numerical RVE determination methods to deliver accurate effective properties of composite materials. This paper demonstrates that common and well-established RVE determination methods, based on studying the convergence rate of the effective properties with respect to the volume element size, are invalid for the case of composites reinforced by randomly oriented fibers and yield erroneous estimates of their effective properties. Following the failure of traditional RVE determination methods, we proposed a new RVE determination criterion that is not based on the average property stability, but its statistical variations. Our new proposed criterion has been shown to be more accurate than other criteria in computing the effective properties of composites for aspect ratios up to 60. Moreover, the proposed criterion does not necessitate a convergence study over the volume element size, hence reducing considerably the RVE determination cost. Finally, our work questions the validity of many published works dealing with composites including heterogeneities of high aspect ratios.     

About the formulation,verification and validation of the hypersonic flow solver Eilmer

We describe the formulation of the gas dynamics and high‐temperature thermochemical modules of the Eilmer code, an open‐source Navier–Stokes solver for transient compressible flow in two and three dimensions. The core gas dynamics formulation is based on finite‐volume cells, and the thermochemical effects are handled with specialised updating schemes that are coupled into the overall time‐stepping scheme. Verification of the code is explored via a number of case studies that use analytic and semi‐analytic solutions as comparison. These include both smooth and shocked flows and are used to demonstrate the order of spatial accuracy of the code. Cases include manufactured solutions for rather abstract inviscid and viscous flow, an idealised detonation wave supported by a curved body, and the transient flow of an idealised but high‐performance shock tube. Validation of the inviscid gas dynamics and thermochemical models is then explored using data from a selection of experimental studies. These studies include ballistic range experiments with chemically‐inert noble gases and high‐temperature chemically‐reacting air. These comparisons show that the code performs well and they provide a lesson in considering a range of experimental data rather than relying upon isolated data points for validation. These verification and validation cases are described in full detail and will be useful for other code developers of high‐temperature compressible flow solvers. Copyright © 2013 John Wiley & Sons, Ltd.     

Local pressure preconditioning method for steady incompressible flows

The convergence and accuracy characteristics of the preconditioned incompressible Euler and Navier–Stokes equations are studied. An object-oriented C++ numerical code has been developed for solving the inviscid and viscous, steady, incompressible flows problems. The code is based on the cell-centred finite volume method. In this scheme, two-dimensional incompressible Euler and Navier–Stokes equations are modified by a robust artificial compressibility (AC) and a local preconditioning matrix of pressure-sensor type. The preconditioned equations are solved with the Jameson's numerical approach, i.e. artificial dissipation and artificial viscosity terms under the form of a fourth- and second-order derivative, respectively. An explicit four-stage Runge–Kutta integration algorithm is applied to obtain the steady-state condition. The computed results include the steady-state solution of flow past the NACA-hydrofoils and a circular cylinder in free stream, for which the numerical results are compared with numerical works of other researchers. Good agreement is observed. The effects of AC parameter, artificial viscosity and dissipation factor, and local preconditioning coefficient on convergence rate and solution accuracy are tested by computing flow over the NACA0012 hydrofoil. In addition, some important design criteria of a preconditioner, such as stiffness reduction, hyperbolicity, symmetrisability, accuracy preservation for M → 0, and M-property have been examined analytically.     

缝洞型多孔介质中多相流的有限体积法数值模拟

本文研究的碳酸盐岩油藏储集体属于缝洞型多孔介质.这类缝洞型多孔介质由裂缝、溶蚀孔洞和低孔隙度低渗透率的基岩组成.裂缝是空隙流体流动的主要通道;溶蚀孔洞大小从几厘米到数米不等,渗透率和孔隙度都很高,是流体主要的储集空间.由于缝洞型多孔介质空隙空间的复杂性和强非均质性,数值计算中基本控制方程的空间离散应采用非结构化网格的计算模型.本文采用有限体积法模拟缝洞型多孔介质中多相流体的流动,并给出了相应的单元中心格式有限体积法的计算公式.裂缝介质和溶洞介质中单元间多相流体的流动考虑为高速非达西流,其质量通量采用Forchheimer定律计算.非线性方程的离散选取全隐式格式,并采用Newton-Raphson迭代进行求解.通过两个二维模型注水驱油的数值模拟,验证了本文方法的有效性.     

Transient flow around an impulsively started cylinder using a dynamic mesh method

The finite volume method (FVM) with a dynamic mesh method (DMM) to deal with the moving boundary was applied to the simulation of two-dimensional incompressible viscous flow past a circular cylinder that is impulsively started into rotation and translation. The non-dimensional rotating to translating speed ratio α is varied from 0.28 to 2.07, with the Reynolds number being 200 for the range of α. The computation covers a period, during which the cylinder translates seven times its diameter. The current scheme handles the impulsively moving boundary directly by DMM, which is implemented using both mesh deforming and local remeshing. The instantaneous asymmetrical flow configurations for various α are presented and compared with the experimental visualizations. Quantitatively, the velocity distributions with drag and lift coefficients are also compared with the experimental and numerical results. Results show that the flow is strongly influenced by the rotation. Comparisons are found to be satisfactory.     

三维空间二阶精度流体体积界面重构方法

给出三维空间网格模板含81个单元的最小二乘流体体积界面重构方法,并和Youngs方法及网格模板含125个单元的最小二乘流体体积界面重构方法进行比较.静态和动态的测试例子均表明:该方法能精确重构任意方向的平面界面,对C²光滑曲面它能达到二阶收敛精度.和网格模板含125个单元的最小二乘流体体积界面重构方法相比,在达到同样网格精度的条件下,减少了计算量,节省了计算时间,提高了计算效率.