船用燃气轮机进气滤清器惯性级内的流场计算和实验验证

本文提出了一套求解船用燃气轮机进气滤清器流道流场的数值方法,成功地计算了流道内流场的状态,给出了各种不同型号流道的气动特性,对指导滤清器的设计有较大的现实意义, 在这套方法中,我们应用上风差分来逼近二维、非定常、粘性、不可压缩流体非守恒型的N-S方程,提出了一种可计算雷诺数高达上万的粘性流的差分方程,考察了这种差分方程的稳定性,收敛性、精度和人工粘性,本文还提出了处理一些边界拐点处壁涡的计算方法,实际算例表明,使用本文提出的差分方程和壁涡处理方法给出的计算结果和实验吻合良好。     

船舶在进出船厢运动中的附加质量的计算

采用三维有限元方法，计算船在进出船厢运动中的附加质量。文中对单圆柱体在无限域和半无限域中运动的附加质量及双圆柱体同心放置时内圆柱运动的附加质量分别进行了计算，计算值与理论值吻合很好；对船在进出船厢运动中的附加质量做了大量的计算，结果表明：随着船与船厢侧壁间距的减小，船的附加质量系数增大；随着船厢中水深的减小，船的附加质量系数增大，随着船离船厢封闭端的距离的减小，附加质量系数增大。     

Numerical simulation and experimental verification of silicone oil flow over magnetic fluid under applied magnetic field

Two-layer flow of magnetic fluid and non-magnetic silicone oil was simulated numerically. The continuity equation, momentum equations, kinematic equation, and magnetic potential equation were solved in two-dimensional Cartesian coordinate. PLIC （piecewise linear integration calculation） VOF （volume of fluid） scheme was employed to track the free interface. Surface tension was treated via a continuous surface force （CSF） model that ensures robustness and accuracy. The influences of applied magnetic field, inlet velocity profile, initial surface disturbance of interface and surface tension were analyzed. The computed interface shapes at different conditions were compared with experimental observation.     

积分单元法及后台阶绕流计算

本文主要介绍解不可压纳维尔—斯托克斯方程的新方法:积分单元法,并用该方法计算了后台阶绕流问题。结果表明,积分单元法是计算后台阶绕流的一种较好的方法。     

考虑参数分布特征的隧道围岩稳定失效概率计算

基于围岩变形破坏原理,建立了围岩稳定极限状态方程的普遍形式,研究了响应面方法求解该方程的过程,得出经典响应面方法仅仅适用于围岩物理力学参数随机分布曲线偏度系数为零的情况。针对围岩物理力学参数统计分布曲线的复杂性,利用统计矩参数对随机变量分布曲线形态特征的控制作用,提出了抽样样本修正方法和具体计算公式,归纳了修正后响应面算法的运行程序。以某公路隧道为例,分别采用经典算法和修正后算法计算了围岩稳定可靠度,然后与蒙特卡洛模拟结果进行了对照,经典方法计算结果与准精确解的绝对误差为1.017%,修正方法的计算结果与准精确解的绝对误差为0.388%。     

爆炸冲击波反射流场的理论计算方法

爆炸冲击波遇到固壁,依次发生正规和非正规反射。本文中基于镜像方法,将爆炸冲击波在固壁反射等效为真实和虚拟爆炸流场的相互作用,建立了波后流场的理论计算方法。首先,假定反射波是以虚拟爆源为中心的圆弧,马赫杆是以爆心在固壁投影点为中心的圆弧。然后,根据爆炸自由场传播规律,利用基于几何近似的方法,建立流场中冲击波结构随时间演化的计算方法,确定任意时刻波后流场区域。最后,利用新发展的叠加模型LAMBR (LAMB revisied),将真实和虚拟爆炸流场进行叠加,给出波后流场中的压力、密度和速度等物理量。通过与数值模拟结果和已有数据进行对比,发现该方法得到的流场物理量分布、峰值等能够反映流场发展的主要规律,从而验证了该理论方法的合理性。而且,该理论方法所需的时间相较于数值模拟大大缩短。     

全机绕流Euler方程多重网格分区计算方法

全机三维复杂形状绕流数值求解只能采用分区求解的方法,本文采用可压缩Euler方程有限体积方法以及多重网格分区方法对流场进行分区计算。数值方法采用改进的van Leer迎风型矢通量分裂格式和MUSCL方法,基于有限体积方法和迎风型矢通量分裂方法,建立一套处理子区域内分界面的耦合条件。各个子区域之间采用显式耦合条件,区域内部采用隐式格式和局部时间步长等,以加快收敛速度。计算结果飞机表面压力分布等气动力特性与实验值进行了比较,二者基本吻合。计算结果表明采用分析“V”型多重网格方法,能提高计算效率,加快收敛速度达到接近一个量级。根据全机数值计算结果和可视化结果讨论了流场背风区域旋涡的形成过程。     

各向同性超弹性本构模型数值计算及验证

现有多种形式的橡胶本构模型试图预测橡胶力学性质,其中部分模型已写入有限元软件中用于仿真计算,还存在较多拟合性较好的模型无法在有限元材料库中直接获得。本文详述了由不变量和主伸长率描写的各向同性超弹性本构模型的数值实现方法,并结合最新的本构模型开发了UHYPER和UMAT子程序。将UHYPER用于有限元实现对多孔橡胶板的拉伸仿真,对比仿真和试验结果,验证子程序的正确性以及评估本构模型预测复杂应变场的准确性;将UMAT用于单轴、等双轴和剪切拉伸的有限元仿真,对比仿真和本构模型理论结果,验证子程序的可靠性。结果表明,有限元仿真结果与理论结果拟合较好,子程序能够契合本构模型的力学描述,所述方法可以用于超弹性材料的数值计算。     

气体泄爆压力分步计算模型及其湍流修正

通过将密闭空间内爆燃泄放过程进行微分,假设各微分时段内爆燃泄放过程均按照先燃烧、再泄放、最后压力平衡的过程独立分步进行,最终得到泄爆压力分步计算模型。同时,在尺寸为2 m×1.2 m×0.6 m的爆炸腔体一端安装击穿压力相同、泄放面积不同的泄爆构件进行泄爆实验,对分步压力计算模型进行验证。对比发现:大面积泄放条件下,2个传感器测得的压力曲线基本重合,均为单峰值曲线,此时模型计算值与实验结果吻合较好;小面积泄放条件下,腔体内压力曲线均为双峰值曲线,由于泄放截面改变加剧口部湍流扰动,使得腔体内部产生压力梯度,近泄爆口处传感器测得的第2个压力峰值要大于腔体内部传感器相应的测量值,经湍流速度修正后的压力计算模型可以较好地描述近泄爆口处的压力变化情况。     

超空泡射弹尾拍分析与计算

对超空泡射弹进行运动学和动力学分析并数学建模,求解耦合非线性微分方程组,得到水下高速超空泡射弹运动特性。数值模拟结果表明,高速超空泡射弹在航行过程中,由于弹体头部和尾部的阻力作用,水平速度随时间迅速衰减。并且射弹的角速度呈周期性往复变化,即尾拍现象。同时由于空泡尺寸的减小导致尾拍幅度逐渐变小。射弹转动惯量越小,角速度变化幅度越平稳,相同时间内尾拍次数减少。发射深度或发射速度越大,尾拍幅度衰减越快。较大的初始角速度也会使射弹角速度很快衰减。     

水下滑翔器位姿估计的改进高斯混合粒子滤波器设计与验证（英文）

近年用于水下滑翔器的低成本导航系统成为研究热点,导航器件的成本与精度之间的折中问题仍然是目前的难题。针对因使用低成本的导航元件而造成低精度位姿估计的问题,提出用于位姿估计的改进高斯混合粒子滤波(IGMPF)方法。用高斯混合模型来估计非高斯噪声,改进的粒子滤波进一步提高位姿估计精度。为了验证其效果,该方法应用于自主设计的水下滑翔器导航系统中并做了车载实验,实验结果表明所提IGMPF方法在实际应用中比传统的EKF和UKF表现更优,姿态角和位移误差比EKF和UKF减小了至少30%。     

Flow field analysis inside a gas turbine trailing edge cooling channel under static and rotating conditions: Effect of ribs

The present work is part of a wider research program which concerns the aero-thermal characterization of cooling channels for the trailing edge of gas turbine blades. The selected passage model is characterized by a trapezoidal cross-section of high aspect-ratio and coolant discharge at the blade tip and along the wedge-shaped trailing edge, where seven elongated pedestals are also installed. In this contribution, a new channel configuration provided with inclined ribs installed inside the radial development region is analyzed, extending the previous results and completing the already available data base, thus providing an overall review of the aero-thermal performance of the considered passage. The velocity field inside the channel was measured by means of 2D and Stereo-PIV techniques in multiple flow planes under static and rotating conditions. The tests were performed under engine similar conditions with respect to both Reynolds (Re = 20,000) and Rotation (Ro = 0, 0.23) numbers. Time averaged flow fields and velocity fluctuation data inside the stationary and rotating channels are analyzed and also critically compared with the data acquired without ribs. In this way the effects on the flow field induced by both rotation and ribs are clearly described. In particular, the ribs modify substantially both the flow field on the channel walls where they are installed and the 3D separation structures that surround the pedestals. If also rotation is taken into account, the relative flow field is characterized by a considerable guiding effect of the ribs coupled with a stronger flow separation on the obstacles that further enhances the heat transfer performances. This behavior was confirmed exploiting the wide thermal data base already available, obtaining a direct link between the observed flow features and the heat transfer performances.     

Mixed compatible element and mixed hybrid incompatible element variational methods in dynamics of viscous barotropic fluids

ＭＩＸＥＤＣＯＭＰＡＴＩＢＬＥＥＬＥＭＥＮＴＡＮＤＭＩＸＥＤＨＹＢＲＩＤＩＮＣＯＭＰＡＴＩＢＬＥＥＬＥＭＥＮＴＶＡＲＩＡＴＩＯＮＡＬＭＥＴＨＯＤＳＩＮＤＹＮＡＭＩＣＳＯＦＶＩＳＣＯＵＳＢＡＲＯＴＲＯＰＩＣＦＬＵＩＤＳＳｈｅｎＸｉａｏ－ｍｉｎｇ（沈孝明...     

Evaluation of an energy relaxation method for the simulation of unsteady,viscous, real gas flows

This study investigates a new energy relaxation method designed to capture the dynamics of unsteady, viscous, real gas flows governed by the compressible Navier–Stokes equations. We focus on real gas models accounting for inelastic molecular collisions and yielding temperature‐dependent heat capacities. The relaxed Navier–Stokes equations are discretized using a mixed finite volume/finite element method and a high‐order time integration scheme. The accuracy of the energy relaxation method is investigated on three test problems of increasing complexity: the advection of a periodic set of vortices, the interaction of a temperature spot with a weak shock, and finally, the interaction of a reflected shock with its trailing boundary layer in a shock tube. In all cases, the method is validated against benchmark solutions and the numerical errors resulting from both discretization and energy relaxation are assessed independently. Copyright © 2004 John Wiley & Sons, Ltd.     

埋地冲洗管线强度计算及影响因素规律分析

研究了冲洗管线在自重、内压、温差和土壤力综合作用下的强度分析方法,给出埋地冲洗管线应力计算公式。采用有限单元法对某冲洗管线重点部位进行力学特性分析,得到其真实应力分布状态,对比了理论计算结果与数值结果误差,证实轴向应力理论计算结果偏保守。对不同壁厚、不同环境温度和不同压力下的管线力学特性进行讨论,结果表明,管道的环向应力起主导作用,且埋深0m~0.75m内管线的环向应力变化明显,因此,当考虑温差及波动内压作用时,该段管线力学特性应引起足够重视;受环境温度影响,管线最大等效应力与总位移规律近似相反,即等效应力最小,总位移近似最大;考虑到等效应力随压力变化规律,不同季节宜采用不同的冲洗压力以满足安全运行要求。     

Computational method for the SPSD of dynamic systems and its application

In this paper, a new method of obtaining the stress power spectral density (SPSD) of a structure that is not through measurement but through imitative computation is presented. First, the dynamic model of a system is established. Then, considering the randomness of the load, the response spectra of the system under different random excitations and parameters are computed by numerical imitation. Last, a concept of unit load stress matrix is put forward. Based on it, the relationship between the SPSD of a structure in a system and the response spectra of the system can be obtained. In particular, the method can be used to compute the SPSD of a new type dynamic structure and to emulate dynamic fatigue life. In the process, an example of the track of a high-speed train is used to compute its SPSD, and the result is used to estimate its fatigue life and to analyze its reliability.     

New numerical schemes based on a criterion for constructing essentially stable and accurate numerical schemes for convection-dominated equations

In order to obtain stable and accurate numerical solutions for the convection-dominated steady transport equations, we propose a criterion for constructing numerical schemes for the convection term that the roots of the characteristic equation of the resulting difference equation have poles. By imposing this criterion on the difference coefficients of the convection term, we construct two numerical schemes for the convection-dominated equations. One is based on polynomial differencing and the other on locally exact differencing. The former scheme coincides with the QUICK scheme when the mesh Reynolds number (Rm) is $\mathop \[{\textstyle{{\rm 8} \over {\rm 3}}}\] $, which is the critical value for its stability, while it approaches the second-order upwind scheme as Rm goes to infinity. Hence the former scheme interpolates a stable scheme between the QUICK scheme at Rm = $\mathop \[{\textstyle{{\rm 8} \over {\rm 3}}}\] $ and the second-order upwind scheme at Rm = ∞. Numerical solutions with the present new schemes for the one-dimensional, linear, steady convection-diffusion equations showed good results.     

A weighted residual method for elastic-plastic analysis near a crack tip and the calculation of the plastic stress intensity factors

In this paper,a weighted residual method for the elastic-plastic analysis near a crack tip is systematically given by taking the model of power-law hardening under plane strain condition as a sample.The elastic-plastic solutions of the crack tip field and an approach based on the superposition of the nonlinear finite element method on the complete solution in the whole crack body field,to calculate the plastic stress intensity factors,are also developed.Therefore,a complete analysis based on the calculation both for the crack tip field and for the whole crack body field is provided.     

考虑磁头表面高度不连续性气膜润滑的数值模拟与有效算法

随着磁头滑块的飞行高度不断降低,给气体润滑方程的数值求解带来了诸如计算时间过长、甚至计算发散等方面的问题。为了获得1Tbit/in2的存储密度,磁头滑块尾部的最小飞行高度接近1.5nm。本文基于作者提出的修正气膜润滑方程的线性流率(LFR)模型,考虑磁头滑块表面高度的不连续性,建立了基于有限体积法的气膜润滑方程离散格式,并把网格自适应技术与多重网格法应用到离散方程的迭代算法中,发展了可模拟最小飞行高度为0.5nm时磁头滑块压力分布的数值模拟方法与有效算法。文中以一个具有复杂表面形状的磁头滑块为例,检验了计算方法与算法的有效性。数值结果表明:在磁头滑块最小飞行高度较低时,必须要考虑滑块表面高度的不连续性,否则就得不到收敛的数值计算结果;与FK-Boltzmann模型相比,LFR模型具有较高的计算效率,采用网格自适应技术与多重网格法能有效地提高求解气膜润滑方程的计算效率。     

输电铁塔四地脚螺栓塔脚板抗拉承载力试验与计算方法研究

针对架空输电线路铁塔中常用的四地脚螺栓塔脚板进行了抗拉承载力试验和理论计算方法研究。首先,选取8个试件进行塔脚板抗拉承载力试验研究,分析塔脚底板板厚及有无加劲板对塔脚板抗拉承载力的影响;其次,通过有限元模拟塔脚板的应力分布情况,并结合“塑性分析,弹性设计”的思想,提出了一种基于屈服线理论的四地脚螺栓塔脚板抗拉承载力计算公式;最后,与试验结果、有限元仿真结果、已有技术规定中的计算结果进行了对比。结果表明:建议公式与试验结果和仿真结果较为吻合,相关技术规定中的计算结果偏于保守,验证了本文建议公式的精确性;底板厚度和加劲板的作用对塔脚板承载力影响较大。研究结果可为输电铁塔四地脚螺栓塔脚板结构设计提供参考。