水平地震时斜面水坝受到的非线性动水压力的全过程的解析解

本文给出在非线性Ｒａｙｌｅｉｇｈ阻尼作用下，由水平地震激发的水库内动水压力变化全过程的解析解：在激发初期其幅值较小，可用小参数法使控制方程线性化，从而求得各阶渐近控制方程的解析解；当动水压力的幅值增大到小参数展开不适用时，用ｖａｎ ｄｅｒｐｏｌ法求得非线性控制方程的渐近解，从而显示了动水压力的幅值从有序变化到发生突然跳跃的多值性的非线性本质。     

动压推力轴承的边界元分析

本文利用Green第二公式,将Reynolds方程转化为沿边界的积分方程,并将非线性项作为自由项的一部分处理,采用常单元离散边界Γ,用迭代技术求出油膜压力分布,与有限差分法和有限元法比较,边界元法的结果更接近解析解.     

流动稳定性弱非线性理论中的问题及其改进

流动稳定性弱非线性理论存在三个主要问题:(1)其摄动参数收敛半径太小且无具体估值;(2)其解有特定结构,一般不能满足很多实际问题的初始条件;(3)当其线性部分不是中性时,Landau系数的确定法不唯一,不知何者应是最好的方法。本文以平面Poiseuille流为例,从理论上解决了问题(1)和(2),并提出了改进办法。又通过将理论和数值模拟结果相比较的方法,澄清了问题(3)。改进了的理论,将可解决许多实际问题。     

Nonlinear rheological behavior of a concentrated spherical silica suspension

Linear and nonlinear viscoelastic properties were examined for a 50 wt% suspension of spherical silica particles (with radius of 40 nm) in a viscous medium, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture. The effective volume fraction of the particles evaluated from zero-shear viscosities of the suspension and medium was 0.53. At a quiescent state the particles had a liquid-like, isotropic spatial distribution in the medium. Dynamic moduli G^* obtained for small oscillatory strain (in the linear viscoelastic regime) exhibited a relaxation process that reflected the equilibrium Brownian motion of those particles. In the stress relaxation experiments, the linear relaxation modulus G(t) was obtained for small step strain (0.2) while the nonlinear relaxation modulus G(t, ) characterizing strong stress damping behavior was obtained for large (>0.2). G(t, ) obeyed the time-strain separability at long time scales, and the damping function h() (–G(t, )/G(t)) was determined. Steady flow measurements revealed shear-thinning of the steady state viscosity () for small shear rates (< ^–1; = linear viscoelastic relaxation time) and shear-thickening for larger (>^–1). Corresponding changes were observed also for the viscosity growth and decay functions on start up and cessation of flow, + (t, ) and _– (t, ). In the shear-thinning regime, the and dependence of ₊(t,) and _–(t,) as well as the dependence of () were well described by a BKZ-type constitutive equation using the G(t) and h() data. On the other hand, this equation completely failed in describing the behavior in the shear-thickening regime. These applicabilities of the BKZ equation were utilized to discuss the shearthinning and shear-thickening mechanisms in relation to shear effects on the structure (spatial distribution) and motion of the suspended particles.Dedicated to the memory of Prof. Dale S. Parson     

Anomalous rolling of spheres down an inclined plane

A sphere in air will roll down a plane that is tilted away from the vertical. The only couple acting about the point of contact between the sphere and the plane is due to the component of the weight of the sphere along the plane, provided that air friction is negligible. If on the other hand the sphere is immersed in a liquid, hydrodynamic forces will enter into the couples that turn the sphere, and the rotation of the sphere can be anomalous, i.e., as if rolling up the plane while it falls. In this paper we shall show that anomalous rolling is a characteristic phenomenon that can be observed in every viscoelastic liquid tested so far. Anomalous rolling is normal for hydrodynamically levitated spheres, both in Newtonian and viscoelastic liquids. Normal and anomalous rolling are different names for dry and hydrodynamic rolling. Spheres dropped at a vertical wall in Newtonian liquids are forced into anomalous rotation and are pushed away from the wall while in viscoelastic liquids, they are forced into anomalous rotation, but are pushed toward the wall. If the wall is inclined and the fluid is Newtonian, the spheres will rotate normally for dry rolling, but the same spheres rotate anomalously in viscoelastic liquids when the angle of inclination from the vertical is less than some critical value. The hydrodynamic mechanisms underway in the settling of circular particles in a Newtonian fluid at a vertical wall are revealed by an exact numerical simulation based on a finite-element solution of the Navier-Stokes equations and Newton's equations of motion for a rigid body.     

Natural oscillations of a jet membrane under hydrostatic pressure

The natural oscillations of a submerged elastic jet membrane modeling an axisymmetric hydrodynamic direct-flow radiator are analyzed. The fundamental frequency is calculated as a function of the geometric and hydrodynamic characteristics of the jet. The principle possibility of adjusting the frequency by changing the elasticity of the fluid, which is a function of hydrostatic pressure, is demonstrated. Numerical and experimental data are compared.__________Translated from Prikladnaya Mekhanika, Vol. 40, No. 12, pp. 92–98, December 2004.     

Brownian dynamics simulations on CPU and GPU with BD_BOX

There has been growing interest in simulating biological processes under in vivo conditions due to recent advances in experimental techniques dedicated to study single particle behavior in crowded environments. We have developed a software package, BD_BOX, for multiscale Brownian dynamics simulations. BD_BOX can simulate either single molecules or multicomponent systems of diverse, interacting molecular species using flexible, coarse-grained bead models. BD_BOX is written in C and employs modern computer architectures and technologies; these include MPI for distributed-memory architectures, OpenMP for shared-memory platforms, NVIDIA CUDA framework for GPGPU, and SSE vectorization for CPU.     

Bifurcation of limit cycles in fluid film bearings

Rotors supported by journal bearings may become unstable due to self-excited vibrations when a critical rotor speed is exceeded. Linearised analysis is usually used to determine the stability boundaries. Non-linear bifurcation theory or numerical integration is required to predict stable or unstable periodic oscillations close to the critical speed. In this paper, a dynamic model of a short journal bearing is used to analyse the bifurcation of the steady state equilibrium point of the journal centre. Numerical continuation is applied to determine stable or unstable limit cycles bifurcating from the equilibrium point at the critical speed. Under certain working conditions, limit cycles themselves are shown to disappear beyond a certain rotor speed and to exhibit a fold bifurcation giving birth to unstable limit cycles surrounding the stable supercritical limit cycles. Numerical integration of the system of equations is used to support the results obtained by numerical continuation. Numerical simulation permitted a partial validation of the analytical investigation.     

爆炸焊接界面波的数值模拟

借助动力学分析软件ANSYS/LS-DYNA 11.0,运用光滑粒子流体动力学方法(SPH方法),建立以Johnson-Cook材料模型和Grüneisen状态方程为基础的热塑性流体力学模型,对同种钢板爆炸焊接界面波进行了数值模拟.结果表明,运用SPH方法可以得到清晰的爆炸焊接界面波形貌.与实验结果的比较表明:模拟误...