减振装置环形阀片局部变载的变形研究

提出了减振器环形节流阀片受局部变载荷作用时挠曲变形量的一种研究方法.根据节流阀片物理模型以及薄板力学的相关理论,推导阀片受局部变载荷作用时挠曲变形量的解析式.与有限元计算结果对比验证了公式以及推导过程的正确性,研究了不同载荷、载荷作用半径以及节流阀片相关结构参数和材料特性对其任意半径处变形量的影响规律.得出的解析式和结论为减振装置阻尼阀的精确设计提供了依据.     

深槽表面气体静压润滑低Reynolds数层流节流效应实验研究

气体介质在润滑间隙流动过程中,沿气体流动方向开设的表面矩形深槽结构内部产生旋涡阻碍气体流动,形成节流效应.通过矩形深沟槽表面静压润滑实验,开展深槽表面和光滑表面的流量特性对比测试,研究了低Reynolds数层流状态下深槽表面间隙的气体节流效应.结果表明:深槽表面产生显著的节流效应,并且随Reynolds数增加而增强;在低Reynolds数层流状态下,深槽表面可以产生节流效应,但是节流效应强度不稳定,当处于完全湍流状态时,节流效应维持定值;间隙尺寸、槽深等结构参数对节流效应影响明显.     

?????????????????????

针对国家标准GB/T18215.1-2000中规范的"差压式流量计示值修正公式",从以下方面进行了讨论: (1)干气体流量检测示值修正公式; (2)湿气体干部分流量检测示值修正公式; (3)气体组分改变时流量检测示值修正公式.国标中规范的修正公式是建立在: "当气体的温度、压力变化时,其密度也随着发生改变,相同的差压值(设计状态与工作状态)所对应的流量值是不同的".本文认为:当气体工作温度、压力等参数偏离设计条件时,同一质量流量在节流装置前后产生的设计状态与工作状态差压是不同的.这一基础要领的分岐,影响着所有修正公式.     

带有环形阀片的阻尼阀可靠性研究

针对跑车试验时油气弹簧中节流阀片受冲击失效的现象,提出了研究阻尼阀水击压强的必要性,建立了带有环形节流阀片的阻尼阀物理模型. 运用水击理论及解偏微分方程的特征线法,创造性地对给定阻尼阀结构进行了水击建模计算. 编程分析节流阀片失效的原因,系统研究了常通孔结构参数变化对阻尼阀所受瞬态冲击及其阻力值的影响规律,提出需从屈服和疲劳强度两个方面入手对阀门可靠性进行校核,并给出计算实例,完善了现有阀系设计的算法. 所得结论可以作为可靠性研究的参考,同时避免了阀片断裂的再次发生.      

半潜式海上浮式风机气动阻尼特性研究

由于海上漂浮式风机具有较大的支撑平台运动,气动阻尼效应对海上漂浮式风机的运动响应带来了重要的影响, 日渐受到相关国内外学者的关注. 为了研究海上浮式风机的气动阻尼特性,本文推导了海上浮式风机气动阻尼力的数学模型,并借助模型实验和数值计算的方法,研究了半潜式海上浮式风机的气动阻尼特性及其作用规律. 结果表明,浮式风机的风轮旋转时的气动阻尼比风轮非旋转状态时更加明显;在作业工况下,气动阻尼对半潜式浮式风机平台的纵荡、纵摇、机舱的运动有明显的抑制作用,且主要体现为对半潜式浮式风机的平台运动固有频率响应的抑制作用,对波频范围的平台运动作用甚微. 其变化规律与风速大小、波浪载荷等有关,在风机的额定工况之前,气动阻尼通常与风速呈正相关关系,但是增长率有逐渐减小的趋势;在控制系统作用下,当入流风速接近或超过风机额定风速时,容易出现气动负阻尼现象,反而进一步强化浮式风机的运动响应,此时通过降低变桨距控制器的比例系数,即降低变桨距控制器的灵敏度,有助于增加海上浮式风机的气动阻尼效果,并且在一定程度上减缓负的气动阻尼的发生,改善海上浮式风机的运动响应.      

Mindlin 矩形微板的热弹性阻尼解析解

首次给出了四边简支的 Mindlin 矩形微板热弹性阻尼的解析解. 基于考虑一阶剪切变形的 Mindlin 板理论和单向耦合热传导理论建立了微板热弹性耦合自由振动控制微分方程. 忽略温度梯度在面内的变化,在上下表面绝热边界条件下求得了用变形几何量表示的温度场的解析解. 进一步将包含热弯曲内力的结构振动方程转化为只包含挠度振幅的四阶偏微分方程. 利用特征值问题之间在数学上的相似性,在四边简支条件下给出了用无阻尼 Kirchhoff 微板的固有频率表示的 Mindlin 矩形微板的复频率解析解,从而利用复频率法求得了反映热弹性阻尼水平的逆品质因子. 最后,通过数值结果定量地分析了剪切变形、材料以及几何参数对热弹性阻尼的影响 规律. 结果表明,Mindlin 板理论预测的热弹性阻尼小于 Kirchhoff 板理论预测的热弹性阻尼. 两种理论预测的热弹性阻尼之间的差值在临界厚度附近十分显著. 另外,随着微板的边/厚比增大,Mindlin 微板的热弹性阻尼最大值单调增大,而 Kirchhoff 微板的热弹性阻尼最大值却保持不变.      

一类阻尼控制半主动隔振系统的解析研究

研究了一类基于相对速度反馈的含立方刚度的单自由度非线性半主动隔振系统.通过平均法得到了系统分别在基于加速度-相对速度反馈的加速度驱动阻尼控制策略、速度-相对速度反馈的天棚阻尼控制策略和位移-相对速度反馈的地棚阻尼控制策略下主共振响应的近似解析解,并利用数值解验证了近似解析解的准确性.通过 Lyapunov 理论对不同控制策略下系统的稳定性进行了分析,讨论了系统参数对控制效果的影响.分析结果表明,对 3 种基于相对速度反馈的控制策略进行解析研究时,切换条件中的控制参数具有相同的表达式;在抑制共振响应振幅方面,基于速度-相对速度反馈的天棚阻尼控制策略在低频时的减振效果最好,而基于加速度-相对速度反馈的加速度驱动阻尼控制策略在高频时的减振效果最优;在抑制瞬态响应振幅方面,基于速度-相对速度反馈的天棚阻尼控制策略的减振效果最好.此类解析研究方法可应用到其他半主动开关控制策略中,为半主动隔振系统的控制策略研究提供了有效的方法和手段.      

ANLYTICAL SOLUTION OF THERMOELASTIC DAMPING IN RECTANGULAR MINDLIN MICRO PLATES 1)

首次给出了四边简支的 Mindlin 矩形微板热弹性阻尼的解析解. 基于考虑一阶剪切变形的 Mindlin 板理论和单向耦合热传导理论建立了微板热弹性耦合自由振动控制微分方程. 忽略温度梯度在面内的变化,在上下表面绝热边界条件下求得了用变形几何量表示的温度场的解析解. 进一步将包含热弯曲内力的结构振动方程转化为只包含挠度振幅的四阶偏微分方程. 利用特征值问题之间在数学上的相似性,在四边简支条件下给出了用无阻尼 Kirchhoff 微板的固有频率表示的 Mindlin 矩形微板的复频率解析解,从而利用复频率法求得了反映热弹性阻尼水平的逆品质因子. 最后,通过数值结果定量地分析了剪切变形、材料以及几何参数对热弹性阻尼的影响 规律. 结果表明,Mindlin 板理论预测的热弹性阻尼小于 Kirchhoff 板理论预测的热弹性阻尼. 两种理论预测的热弹性阻尼之间的差值在临界厚度附近十分显著. 另外,随着微板的边/厚比增大,Mindlin 微板的热弹性阻尼最大值单调增大,而 Kirchhoff 微板的热弹性阻尼最大值却保持不变.     

DYNAMICAL ANALYSIS ON A KIND OF SEMI-ACTIVE VIBRATION ISOLATION SYSTEMS WITH DAMPING CONTROL 1)

研究了一类基于相对速度反馈的含立方刚度的单自由度非线性半主动隔振系统.通过平均法得到了系统分别在基于加速度-相对速度反馈的加速度驱动阻尼控制策略、速度-相对速度反馈的天棚阻尼控制策略和位移-相对速度反馈的地棚阻尼控制策略下主共振响应的近似解析解,并利用数值解验证了近似解析解的准确性.通过 Lyapunov 理论对不同控制策略下系统的稳定性进行了分析,讨论了系统参数对控制效果的影响.分析结果表明,对 3 种基于相对速度反馈的控制策略进行解析研究时,切换条件中的控制参数具有相同的表达式;在抑制共振响应振幅方面,基于速度-相对速度反馈的天棚阻尼控制策略在低频时的减振效果最好,而基于加速度-相对速度反馈的加速度驱动阻尼控制策略在高频时的减振效果最优;在抑制瞬态响应振幅方面,基于速度-相对速度反馈的天棚阻尼控制策略的减振效果最好.此类解析研究方法可应用到其他半主动开关控制策略中,为半主动隔振系统的控制策略研究提供了有效的方法和手段.     

消能锥阀结构改进的三维数值模拟研究

对常规消能锥阀体型进行了改进,提出了一种新型的嵌套式消能锥阀.该新阀将原来通过阀体的一股环状射流变为两股互不平行的环状射流,通过这两股高速水流在阀口的碰撞提高消能率、改善流态.采用可压缩的两相流模型,辅以Realizable k–ε湍流模型来模拟锥阀附近水流的水力特性,对于自由水面的处理采用了VOF法,并进行压强场和流速场及消能情况的对比分析.结果显示:在上游进口端控制断面的水流流量相同的情况下使用嵌套式消能锥阀,阀前进水管道的压强可由原来的25kPa减小到15kPa,减小了40%;阀体所承受的压力由35kPa减小到10kPa,减小了57%;下游管道中水流落点之后的沿程断面平均流速减小了约0.1m/s.在上游进口端具有相同的总水头的情况下,嵌套式消能锥阀对流速水头及压强水头的消能能力都要比使用常规消能锥阀好.     

Simulation and Modelling of a Skewed Turbulent Channel Flow

A time-dependent three-dimensionally skewed flow is investigated using direct numerical simulations of the incompressible Navier-Stokes equations. The effect on the instantaneous and mean turbulent field is investigated. Instantaneous flowfields reveal that the skewing has the effect of initially reducing the strength and height of quasi-streamwise vortices of both signs of rotation with respect to the skewing. A mechanism for this process is put forward. The mean flowfields show drops in turbulence quantities such as turbulence kinetic energy. In addition to this, two-equation turbulence modelling of the flow is carried out. This highlights a deficiency, in that the standard turbulence models are unable to capture the drop in turbulence intensity due to the skewing. A modification based on the exact dissipation equation is found to significantly improve the model behaviour for this flow. This revised version was published online in July 2006 with corrections to the Cover Date.     

THE HISTORICAL CONTRIBUTION OF BERNOULLI'S EQUATION TO THE ESTABLISHMENT OF FLUID MECHANICS THEORY

著名的理想流体定常流动的能量方程即伯努利方程,自建立以来在流体力学领域中贡献卓著。本文依据伯努利方程的建立内涵,阐述了其在流体静力学、定常孔口出流、皮托管测速、文丘里管流量和翼型绕流等具体流动中的成功应用。同时,进一步说明了由伯努利方程建立提出的局部跟随流体质点的建模思想,被欧拉概括为描述流体运动的流场法,是建立欧拉方程组和N-S方程组的基本依据,也为后来湍流理论、边界层理论、气动噪声等理论的建立奠定了基础。     

Turbulence energy balance in reacting turbulent flows

The turbulence accompanying combustion and the propagation of detonation waves in gases has been studied theoretically and experimentally in many papers [1–8]. The attention of researchers has been concentrated on essential questions like how the turbulent flow field interacts with the kinetics of the chemical reaction and to what extent the process of chemical change is intensified, and how the turbulence itself is deformed by the heat released and the accompanying expansion of the gases. The various mechanisms proposed for these phenomena are based on various hypotheses concerning the structure of the combusion zone and the determinative stage of the interaction of the turbulence with the chemical-reaction kinetics. The mechanism of turbulence generation by combustion proposed in a number of papers [3–6] is based on the observation in turbulent flow of a weakly curved flickering laminar flame. This gives rise to a nonuniform flow field of the gas, part of the energy of which goes over into the energy of turbulent fluctuations. Other authors [7, 8] considered the turbulence field to interact with the chemical-reaction kinetics via a volume mechanism and suggested a criterion of turbulence intensification based on certain physical considerations, e.g., the condition for the intensification of thermogaskinetic oscillations proposed by Rayleigh [9]. In the present paper the problem is analyzed by introducing Kolmogorov's general equation for the turbulence energy balance in reacting turbulent flows [10]. In accordance with, this equation the turbulence energy can vary due to energy exchange between the turbulent motion and the mean gas flow as a result of the work on turbulent mass transport in the acceleration field of the mean flow, and due to the effect of pressure fluctuations on the rate of thermal expansion from the chemical reaction. Each of these effects is considered and analyzed.     

Measurements of temperature,density, pressure,and their fluctuations in supersonic turbulence using laser-induced fluorescence

A laser-induced fluorescence (LIF) method has been developed that provides simultaneous measurements of temperature, density, and their fluctuations owing to turbulence in unheated compressible flows. Pressure and its fluctuations are also deduced using the equation of state. Fluorescence is induced in nitric oxide that has been seeded into a nitrogen flow in concentrations of 100 ppm. Measurements are obtained from each laser pulse, with a spatial resolution of 1 mm and a temporal resolution of 125 ns. The method was applied to a supersonic, turbulent, boundary-layer flow with a free-stream Mach number of 2. For stream conditions in the range from 150–300 K and 0.3–1 atm, temperature is measured with an uncertainty of approximately 1% rms, while density and pressure uncertainties are approximately 2% rms.     

Structure of turbulent pipe flow

The problem of turbulent flow in a straight circular pipe is solved. We consider a system consisting of the equation of motion, the equation for the turbulence energy, the expression relating the turbulence coefficient with the turbulence scale, and the integral formula for determining the turbulence scale. A numerical solution is presented for this closed system of equations for turbulent flow. The results of calculations are compared with experimental data.     

不可压N-S方程高效算法及二维槽道湍流分析

构造了基于非等距网格的迎风紧致格式,并将其与三阶精度的Adams半隐方法相结合,构造了求解不可压N－S方程高效算法。该算法利用基于交错网格的离散形式的压力Poisson方程求解压力项,解决了边界处的残余散度问题;同时还利用快速Fourier变换将方程的隐式部分解耦,离散后的代数方程组利用追赶法求解,大大减少了计算量。通过对二维槽道流动的数值模拟,证实了所构造的数值方法具有精度高,稳定性好,能抑制混淆误差等优点,同时具有很高的计算效率,是进行壁湍流直接数值模拟的有效方法。在数值模拟的基础上对二维槽道流动进行了分析,得到了Reynolds数从6000到15000的二维流动饱和态解（所谓“二维槽道湍流”）;定性及定量结果均与他人的数值计算结果吻合十分理想。对流场进行了分析,指出了“二维湍流”与三维湍流统计特性的区别。     

一类平衡大气边界层边界条件研究

基于标准k-ε湍流模型,首先利用湍流粘度方程和剪切应力在整个边界层内恒定的假设,推导出一类耗散率表达式,并根据常用的湍动能入口剖面方程以及平均风速剖面方程,计算获得相应的耗散率方程;然后在输运方程中添加自定义源项,通过已经确定的平均速度方程、湍动能方程、耗散率方程计算得到相应输运方程的自定义源项表达式,并进行空风洞数值模拟,从而得到了一类满足平衡大气边界层的来流边界条件．通过将这种边界条件与由湍流平衡条件得到的边界条件进行比较,表明本方法获得的边界条件更适用．并且,本方法无需考虑修正壁面函数和修正湍流模型常数,因而计算更为简单,可为平衡大气边界层的研究提供一种新的思路．     

An implementation of the Spalart–Allmaras DES model in an implicit unstructured hybrid finite volume/element solver for incompressible turbulent flow

In this paper, we describe an implicit hybrid finite volume (FV)/element (FE) incompressible Navier–Stokes solver for turbulent flows based on the Spalart–Allmaras detached eddy simulation (SA‐DES). The hybrid FV/FE solver is based on the segregated pressure correction or projection method. The intermediate velocity field is first obtained by solving the original momentum equations with the matrix‐free implicit cell‐centered FV method. The pressure Poisson equation is solved by the node‐based Galerkin FE method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centers and the auxiliary variable at vertices, making the current solver a staggered‐mesh scheme. The SA‐DES turbulence equation is solved after the velocity and the pressure fields have been updated at the end of each time step. The same matrix‐free FV method as the one used for momentum equations is used to solve the turbulence equation. The turbulence equation provides the eddy viscosity, which is added to the molecular viscosity when solving the momentum equation. In our implementation, we focus on the accuracy, efficiency and robustness of the SA‐DES model in a hybrid flow solver. This paper will address important implementation issues for high‐Reynolds number flows where highly stretched elements are typically used. In addition, some aspects of implementing the SA‐DES model will be described to ensure the robustness of the turbulence model. Several numerical examples including a turbulent flow past a flat plate and a high‐Reynolds number flow around a high angle‐of‐attack NACA0015 airfoil will be presented to demonstrate the accuracy and efficiency of our current implementation. Copyright © 2008 John Wiley & Sons, Ltd.