Global composite element iteration for analysis of seepage free surface

ＩｎｔｒｏｄｕｃｔｉｏｎＳｅｅｐａｇｅｆｒｅｅｓｕｒｆａｃｅｉｓｏｆｔｅｎｅｎｃｏｕｎｔｅｒｅｄｉｎｔｈｅｃｏｕｒｓｅｏｆｓｅｅｐａｇｅｐｒｏｂｌｅｍｔｈｒｏｕｇｈｍｅｄｉａｏｆｓｌｏｐｅ,ｅａｒｔｈｄａｍａｎｄｔｕｎｎｅｌｅｎｇｉｎｅｅｒｉｎｇ．Ｔｈｅｐｈｅｎｏｍｅｎａｏｆｓｅｅｐａｇｅｆｒｅｅｓｕｒｆａｃｅｉｓｏｎｅｏｆｔｈｅｉｍｐｏｒｔａｎｔａｎｄｄｉｆｆｉｃｕｌｔｔｏｐｉｃｓｏｆｒｏｃｋｍａｓｓｈｙｄｒｏｍｅｃｈａｎｉｃｓ［１～３］,ａｎｄｑｕｉｔｅａｎｉｍｐｏｒｔａｎｔｒｅｓｅａｒｃｈ…     

隧道洞口段非均质及各向异性土质裂缝仰坡稳定性上限解

针对忽略岩土体非均质和各向异性将导致边坡稳定性评价产生误差的问题,应用极限分析上限理论及抗剪强度系数折减法,推导土体强度非均质和各向异性影响下隧道洞口含裂缝仰坡稳定性解析式,探究土体强度非均质和各向异性对仰坡稳定性系数、坡顶裂缝位置、隧道拱顶失稳范围及仰坡安全系数的影响。结果表明,裂缝深度及坡角越大,仰坡稳定性系数越小;非均质系数越大和各向异性系数越小,维持仰坡稳定的临界坡高越大;非均质系数及各向异性系数越大,裂缝距坡顶边缘越远,隧道拱顶失稳范围越大;非均质系数增大有利于仰坡稳定,而各向异性系数越大仰坡越易失稳。     

对蓬安县三元水库大坝外坡裂缝产生的原因及整治工程的探讨

本文通过对蓬安县三元水库大坝外坡裂缝产生的原因分析，进而运用岩、土力学的理论与方法反算推论出了裂缝的产生不是深层坝体土滑动引起的，而是局部浅层坝体土在水的诱发作用下作蠕滑-拉裂变形所产生的。在明确了大坝裂缝产生的原因后，提出了一系列相应整治工程措施，以达到彻底根治裂缝和稳固大坝的目的。     

MODELLING OF WAVE PROPAGATION IN THE NEARSHORE REGION USING THE MILD SLOPE EQUATION WITH GMRES-BASED ITERATIVE SOLVERS

The mild slope equation in its linear and non-linear forms is used for the modelling of nearshore wave propagation. The finite difference method is used to descretize the governing elliptic equations and the resulting system of equations is solved using GMRES-based iterative method. The original GMRES solution technique of Saad and Schultz is not directly applicable to the present case owing to the complex coefficient matrix. The simpler GMRES algorithm of Walker and Zhou is used as the core solver, making the upper Hessenberg factorization unneccessary when solving the least squares problem. Several preconditioning-based acceleration strategies are tested and the results show that the GMRES-based iteration scheme performs very well and leads to monotonic convergence for all the test-cases considered.     

考虑剪切效应的旋转FGM楔形梁刚柔耦合动力学建模与仿真

本文对一类中心刚体-柔性梁系统在大范围转动下的刚柔耦合动力学问题进行了研究. 柔性梁为功能梯度材料(functionally graded materials, FGM)楔形变截面梁,材料体积分数在梁轴向呈幂律分布变化. 以弧长坐标来描述柔性FGM梁的几何位移关系,分别使用倾角和拉伸应变变量描述柔性梁的横向弯曲和纵向拉伸变形,并计及剪切效应. 采用假设模态法离散变形场,运用第二类拉格朗日方程进行方程推导,得到系统考虑剪切效应的刚柔耦合动力学模型. 基于全新的刚柔耦合动力学建模理论,研究不同轴向材料梯度分布的FGM楔形梁,通过数值仿真计算,分析讨论不同的转速、梯度分布规律以及变截面参数对系统动力学特性的影响. 结果表明,剪切效应对大高跨比的FGM楔形梁的变形影响较为明显,不容忽略;材料梯度分布规律和截面参数的选取均会对旋转FGM楔形梁的动力学响应和频率产生较大影响. 本文提出的考虑剪切效应的倾角刚柔耦合动力学模型是对以往非剪切模型的进一步完善,可应用于工程中的 Timoshenko梁结构的动力学问题求解.      

EAST液态锂限制器分配盒参数对锂流动性影响

为考察溢流孔结构几何参数及壁面电导率对液态锂流动的影响,通过数值模拟方法对不同尺寸、形状及壁面电导率的溢流孔中液态锂的流动进行了数值模拟。结果表明：溢流孔的截面宽度越小,其出口的速度分布越均匀,速度峰值也越小,但溢流孔进出口之间的压力降会显著增加;且随着壁面电导率增大,MHD压降也随之增大。对于圆形及椭圆形截面的溢流孔,其出口的速度分布远比矩形溢流孔均匀。另外,关于限制器热平衡的计算有助于保证限制器的安全稳定运行。结果对分配盒几何参数的设计、电磁泵及冷却系统的设计具有重要的指导意义。     

凹边坡参数对边坡稳定性的影响

抽象出凹边坡的几何模型,采用量纲分析的方法,给出凹边坡几何参数和物理力学参数对边坡稳定系数的影响通式．在FLAC定义的边坡稳定系数意义下,经过数万次计算,采用二次回归的方法,确定了通式中的影响指数,形成经验公式．结果表明,受Mohr-Coulomb和拉应力两个强度准则控制的材料,形成的凹边坡的稳定性受无量纲粘聚力、边坡角、内摩擦角的控制．凹形边坡有利于边坡的稳定,边坡稳定性与岩土的抗拉强度关系不大．这一稳定系数公式对于形状规则均匀性好且符合Mohr-Coulomb破坏准则的凹边坡可以给出精度高的估计．对于其它边坡的稳定性评价还需要做有针对性的工作．     

Flux Limiter Lattice Boltzmann Scheme Approach to Compressible Flows with Flexible Specific-Heat Ratio and Prandtl Number

We further develop the lattice Boltzmann (LB) model [Physica A 382 (2007) 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar--Gross-Krook (BGK) collision term in the LB equation, which makes the model suitable for simulating flows with different Prandtl numbers. Secondly, the flux limiter finite difference (FLFD) scheme is employed to calculate the convection term of the LB equation, which makes the unphysical oscillations at
discontinuities be effectively suppressed and the numerical dissipations be significantly diminished. The proposed model is validated by recovering results of some well-known benchmarks, including (i) The thermal Couette flow; (ii) One- and two-dimensional Riemann problems. Good agreements are obtained
between LB results and the exact ones or previously reported solutions. The flexibility, together with the high accuracy of the new model, endows the proposed model considerable potential for tracking some long-standing problems and for investigating nonlinear nonequilibrium complex systems.     

Extension of the DSMC method to high pressure flows

A collision-limiter method, designated as equilibrium direct simulation Monte Carlo (eDSMC), is proposed to extend the DSMC technique to high pressure flows. The method is similar to collision-limiter schemes considered in the past with the important distinction that for inviscid flows, equilibrium is enforced in the entire flow by providing a sufficient number of collisions, based on pre-simulation testing. To test the method with standard DSMC and Navier–Stokes (NS) methods, axi-symmetric nozzle and embedded-channel flows are simulated and compared with experimental temperature data and pre-existing calculations, respectively. The method is shown to agree with third-order Eulerian nozzle flows and first-order channel flows. Chapman–Enskog theory is utilized to predict the range of initial conditions where eDSMC is potentially useful for modeling flows that contain viscous boundary layer regions. Comparison with supersonic nozzle data suggests that the eDSMC method is not adequate for capturing the large variation in flow length scales occurring in supersonic expansions into a vacuum. However, when eDSMC is used in combination with the baseline-DSMC method a near-exact solution is obtained with a considerable computational savings compared to the exact DSMC solution. Viscous flow channel calculations are found to agree well with an exact Navier–Stokes (NS) calculation for a small Knudsen number case as predicted by Chapman–Enskog theory.     

Impactless biped walking on a slope

Walking without impacts has been considered in dynamics as a motion/force control problem. In order to avoid impacts, an approach for both the specified motion of the biped and its ground reaction forces was presented yielding a combined motion and force control problem. As an application, a walker on a horizontal plane has been considered. In this paper, it is shown how the control of the ground reaction forces and the energy consumption depend on the gradient of a slope. The biped dynamics and the constraints within the biped system and on the ground are discussed. A motion control synthesis is developed using the inverse dynamics principle proven to be most efficient for human walking research, too. The impactless walking with controlled legs is illustrated by a seven-link biped. The “flying” biped has nine degrees of freedom, with six control inputs. During locomotion, the standing leg has three scleronomic constraints, and the trunk has three rheonomic constraints. However, there are three rheonomic constraints for the prescribed leg motion or three scleronomic constraints for reaction forces of the trailing leg, respectively. The nominal control action for impactless walking can be precomputed and stored. The model proposed allows the investigation of several problems: uphill and downhill walking, optimization of step length, stiction of the feet on the slope and many more. All these findings are also of interest in biomechanics.