大型渡槽动力设计方法研究

根据大型渡槽结构的特点,在缺乏适用的抗震分析模型前提下,提出支撑-渡槽-水伪三维流固耦合动力分析方法。模型基于结构动力分析理论和流固耦合分析理论,整体结构分析采用Newmark方法进行动力时程分析,分析中考虑渡槽结构整体受力性能,包括下部支撑结构、地基等多种因素的影响作用;渡槽-水耦合体子结构分析采用任意拉格朗日-欧拉(Arbitrary Lagrangian Eulerian,ALE)有限元方法进行动力分析,考虑渡槽-水的耦合相互作用,模拟渡槽中水体的大幅晃动作用,相应得出渡槽中动力压力值及其分布。以某排架支撑矩形渡槽为例,采用本文提出的伪三维模型进行地震时程分析,计算结果与三维模型计算结果吻合良好。本文提出的伪三维动力分析方法计算效率高,满足渡槽结构抗震设计要求,便于工程应用,是大型渡槽结构动力分析研究和抗震设计的实用静、动力分析方法。     

An efficient formulation based on the Lagrangian method for contact–impact analysis of flexible multi-body system

In this paper,an efficien formulation based on the Lagrangian method is presented to investigate the contact–impact problems of f exible multi-body systems.Generally,the penalty method and the Hertz contact law are the most commonly used methods in engineering applications.However,these methods are highly dependent on various non-physical parameters,which have great effects on the simulation results.Moreover,a tremendous number of degrees of freedom in the contact–impact problems will influenc thenumericalefficien ysignificantl.Withtheconsideration of these two problems,a formulation combining the component mode synthesis method and the Lagrangian method is presented to investigate the contact–impact problems in fl xible multi-body system numerically.Meanwhile,the finit element meshing laws of the contact bodies will be studied preliminarily.A numerical example with experimental verificatio will certify the reliability of the presented formulationincontact–impactanalysis.Furthermore,aseries of numerical investigations explain how great the influenc of the finit element meshing has on the simulation results.Finally the limitations of the element size in different regions are summarized to satisfy both the accuracy and efficien y.     

Characterization of particle-laden, confined swirling flows by phase-doppler anemometry and numerical calculation

The particle dispersion characteristics in a confined swirling flow with a swirl number of approx. 0.5 were studied in detail by performing measurements using phase-Doppler anemometry (PDA) and numerical predictions. A mixture of gas and particles was injected without swirl into the test section, while the swirling airstream was provided through a co-flowing annular inlet. Two cases with different primary jet exit velocities were considered. For these flow conditions, a closed central recirculation bubble was established just downstream of the inlet.

The PDA measurements allowed the correlation between particle size and velocity to be obtained and also the spatial change in the particle size distribution throughout the flow field. For these results, the behaviour of different size classes in the entire particle size spectrum, ranging from about 15 to 80 μm, could be studied, and the response of the particles to the mean flow and the gas turbulence could be characterized. Due to the response characteristics of particles with different diameters to the mean flow and the flow turbulence, a considerable separation of the particles was observed which resulted in a streamwise increase in the particle mean number diameter in the core region of the central recirculation bubble. For the lower particle inlet velocity (i.e. low primary jet exit velocity), this effect is more pronounced, since here the particles have more time to respond to the flow reversal and the swirl velocity component. This also gave a higher mass of recirculating particle material.

The numerical predictions of the gas flow were performed by solving the time-averaged Navier-Stokes equations in connection with the well known kε turbulence model. Although this turbulence model is based on the assumption of isotropic turbulence, the agreement of the calculated mean velocity profiles compared to the measured gas velocities is very good. The gas-phase turbulent kinetic energy, however, is considerably underpredicted in the initial mixing region. The particle dispersion characteristics were calculated by using the Lagrangian approach, where the influence of the particulate phase on the gas flow could be neglected, since only very low mass loadings were considered. The calculated results for the particle mean velocity and the mass flux are also in good agreement with the experiments. Furthermore, the change in the particle mean diameter throughout the flow field was predicted approximately, which shows that the applied simple stochastic dispersion model also gives good results for such very complex flows. The variation of the gas and particle velocity in the primary inlet had a considerable impact on the particle dispersion behaviour in the swirling flow and the particle residence time in the central recirculation bubble, which could be determined from the numerical calculations. For the lower particle inlet velocity, the maximum particle size-dependence residence time within the recirculation region was considerably shifted towards larger particles.     

Analysis of unsteady gas–liquid flows in a rectangular tank: Comparison of Euler–Eulerian and Euler–Lagrangian simulations

We study the dynamics of gas–liquid flows experimentally and computationally in a rectangular bubble column where the gas source is introduced at the corner. The flow in this reactor is complex and inherently unsteady in nature. The two-dimensional liquid phase velocity field is calculated by an Eulerian approach solving the unsteady Reynolds Averaged Navier Stokes equations. The conservation equations are closed using a two parameter turbulence model. The two-way coupling was accounted for by adding source terms in the conservation equations of the continuous phase to take into account the interaction with the dispersed phase. Bubble tracking is achieved through a Lagrangian approach. Here the equations of motion are solved taking into account the drag, pressure, buoyancy and gravity forces. The time-averaged flows along with the variables which characterize turbulence are analyzed for a wide range of gas flow-rates using Euler–Lagrangian simulations. These simulation predictions are validated with Euler–Eulerian simulations where the gas-phase distribution is captured as a void fraction and PIV experiments. The motion of bubbles induces turbulence in the flow. The applicability of two parameter models for turbulence like the standard k–ε model on time-averaged flow properties is addressed. From the results of the time averaged velocity field, turbulence intensity, turbulent viscosity and gas hold-up profiles, it is concluded that the Euler–Lagrangian model is applicable at lower gas flow-rates. The Euler–Eulerian approach was found to be valid at lower as well as higher gas flow-rates.     

Kronecker product and a new matrix form of Lagrangian equations with multipliers for constrained multibody systems

The automatic derivation of motion equations is an important problem of multibody system dynamics. Firstly, an overview of the matrix calculus related to Kronecker product of two matrices is presented. A new matrix form of Lagrangian equations with multipliers for constrained multibody systems is then developed to demonstrate the usefulness of Kronecker product of two matrices in the study of dynamics of multibody systems. Finally, the equations of motion of mechanisms are derived using the proposed matrix form of Lagrangian equations as application examples.     

Free‐surface flow over a semi‐circular obstruction

The fully non‐linear free‐surface flow over a semi‐circular bottom obstruction was studied numerically in two dimensions using a mixed Eulerian–Lagrangian formulation. The problem was solved in the time domain that allows the prediction of a number of transient phenomena, such as the generation of upstream advancing solitary waves, as well as the simulation of wave breaking. A parametric study was performed for a range of values of the depth‐based Froude number up to 2.5 and non‐dimensional obstacle heights, α up to 0.9. When wave breaking does not occur, three distinct flow regimes were identified: subcritical, transcritical and supercritical. When breaking occurs it may be of any type: spilling, plunging or surging. In addition, for values of the Froude number close to 1, the upstream solitary waves break. A systematic study was undertaken to define the boundaries of each type of breaking and non‐breaking pattern and to determine the drag and lift coefficients, free‐surface profile characteristics and transient behavior. Copyright © 1999 John Wiley & Sons, Ltd.     

The Kinematics and the Full Minimal Dynamic Model of a 6-DOF Parallel Robot Manipulator

In this paper we present a particular architecture of parallel robots which has six-degrees-of-freedom (6-DOF) with only three limbs. The particular properties of the geometric and kinematic models with respect to that of a classical parallel robot are presented. We show that inverse problems have an analytical solution. However, to solve the direct problems, an efficient numerical procedure which needs to inverse only a 3 × 3 passive Jacobian matrix is proposed. In a second step, dynamic equations are derived using the Lagrangian formalism where the joint variables are passive and active joint coordinates. Based on the geometrical properties of the robot, the equations of motion are derived in terms of only nine coordinates related by three kinematic constraints instead of 18 joint coordinates. The computational cost of the dynamic model obtained is reduced by using a minimum set of base inertial parameters.     

Investigation on time scales in a low Reynolds number jet flow using particle-tracking velocimetry

The near flow field of an axisymmetric water jet at Reynolds numbers between 2000 and 5000 is investigated using Particle-Tracking Velocimetry. Measurements are taken in the longitudinal section (along the mean flow) and in cross-sections (orthogonal to the mean flow). From the former, correlation coefficients of the two in-plane velocity components in a Lagrangian framework are obtained: thus Lagrangian integral scales can be computed. Those of the streamwise velocity (axial) component increase on moving away from the centreline, whereas the opposite happens for the vertical velocity (radial) component: integral time scales of the two components are almost equal at the interface between jet and ambient fluids. On the other hand, integral scales are almost constant or increase slightly with the axial direction. In cross-sections, fluid ejection and injection from the jet centreline are observed to be connected to counter-rotating vortices (mushroom): their number and size change with Reynolds number in agreement with results from other authors. The maximum ejection velocity (orthogonal to the mean jet flow), at 3 nozzle diameters downstream of the outlet, is found to be one half of the mean outlet velocity.     

High-order Lagrangian cell-centered conservative scheme on unstructured meshes

A high-order Lagrangian cell-centered conservative gas dynamics scheme is presented on unstructured meshes. A high-order piecewise pressure of the cell is intro- duced. With the high-order piecewise pressure of the cell, the high-order spatial discretiza- tion fluxes are constructed. The time discretization of the spatial fluxes is performed by means of the Taylor expansions of the spatial discretization fluxes. The vertex velocities are evaluated in a consistent manner due to an original solver located at the nodes by means of momentum conservation. Many numerical tests are presented to demonstrate the robustness and the accuracy of the scheme.