内损耗因子对铝合金壁板高频疲劳寿命影响的研究

基于统计能量分析方法计算高超声速飞行器铝合金壁板高频声振响应,分析了在不同内损耗因子情况下壁板均方根应力和危险点响应集中因子。进而采用频域疲劳寿命方法,结合零阶矩应力谱法计算了壁板的高频疲劳寿命,研究了在不同内损耗因子情况下铝合金壁板高频疲劳寿命的变化规律。研究结果表明:随着内损耗因子的增大,壁板的均方响应逐渐减小,但响应集中因子逐渐增大,最终导致壁板疲劳寿命出现先变长后缩短的变化趋势;随着内损耗因子的增大,高频段应力响应对疲劳寿命的贡献量明显增大。     

宽带噪声作用下黏弹性板的矩Lyapunov指数

主要研究了在超音速流中受宽带噪声作用的黏弹性板随机振动系统的矩Lyapunov指数.首先, 采用vonKarman板弯理论, 活塞理论以及Galerkin近似法建立了两个自由度耦合的系统运动的随机微分方程. 其次, 应用随机平均法将四维系统降为二维系统. 接着, 对系统采用极坐标变换,通过Girsanov定理和Feynmann-Kac公式得到后向微分算子. 通过对特征函数进行正交Fourier余弦级数展开得到系统矩Lyapunov指数的近似解析式. 并通过MonteCarle仿真得到系统矩Lyapunov指数的数值解验证了近似解析式的可信性. 最后研究了系统参数、气动力参数以及随机噪声谱密度对黏弹性板稳定性的影响.      

基于精细积分技术的非线性动力学方程的同伦摄动法

将精细积分技术（PIM）和同伦摄动方法（HPM）相结合,给出了一种求解非线性动力学方程的新的渐近数值方法。采用精细积分法求解非线性问题时,需要将非线性项对时间参数按Taylor级数展开,在展开项少时,计算精度对时间步长敏感;随着展开项的增加,计算格式会变得越来越复杂。采用同伦摄动法,则具有相对筒单的计算格式,但计算精度较差,应用范围也限于低维非线性微分方程。将这两种方法相结合得到的新的渐近数值方法则同时具备了两者的优点,既使同伦摄动方法的应用范围推广到高维非线性动力学方程的求解,又使精细积分方法在求解非线性问题时具有较简单的计算格式。数值算例表明,该方法具有较高的数值精度和计算效率。     

顺桨风力机塔叶干涉翼型非稳气动力时频特性

由于风力机叶片与塔筒流场相互干涉,实际气动力与理想情况存在较大差异,这种干涉作用造成的气动力差异给叶片与塔筒结构可靠性带来不可忽视的影响.以翼型DU91-W2-250为研究对象,采用瞬态数值分析与本征正交分解方法,考虑叶片和塔筒流场相互干涉作用,分析顺桨工况翼型非稳气动力时频特性及其影响规律,量化不同雷诺数下塔叶相对位置及几何参数对气动力均值、波动幅度和频率的影响程度,通过流场模态能量分布形态分析,揭示流场干涉对气动力的影响机制.结果表明,翼型气动中心至塔筒几何中心的垂直距离、水平距离以及塔筒直径相对于翼型弦长的无量纲参数y*,x*和D*对气动力均有不同程度影响,其中y*对升阻力系数均值影响最大,对频率无明显影响,y*绝对值越大,C_l均值越接近单翼型C_l值,y*绝对值越小升阻力系数波动幅度越大,y*从-12增大到12,升力系数均值最小值为-0.48,最大值为1.16;x*减小和D*增大,反向阻力均值增大,波动幅度增大,波动频率略有下降,当x*小于临界值5时,带塔翼型阻力均值反向;在计算范围内,带塔翼型升力系数均值相对于单翼型升力系数最大偏差为...     

浮动冲击平台横向冲击谱受外形结构的影响

基于ABAQUS软件中的声固耦合法,采用设备、平台一体化分析方法,展开对中型浮动冲击平台的研究,探讨平台外形结构形式对平台冲击环境的影响,并提出了提高平台横向冲击谱谱值的设计方案。首先对外形结构进行初步设计,分析其对平台冲击谱的影响,并找到决定影响程度的关键因素。然后针对结构进行优化,使其更大幅度的提高平台冲击谱谱值。计算表明:在平台外部下方加装挡板结构不会明显影响平台垂向冲击谱,但可以提高平台横向冲击谱;舷侧底端向下延伸加装竖直挡板由于受到冲击波绕射和阻力的影响,增加横向谱值的效果不是很明显;平台底部流线型挡板可以有效增加平台对爆炸载荷的接收效果,同时可以尽量减小阻力影响,从而明显提高平台横向冲击谱谱值。     

同轴圆柱间滑移旋转流动问题的同伦解析解

从理论上研究了具有延伸柱面的同轴圆柱间滑移流动问题.通过引入适当的相似变换将控制方程组转化为一类非线性边值问题, 利用同伦分析方法首次获得问题的近似解析解, 分析讨论滑移参数、雷诺数和内外筒半径比对流动的影响.结果表明: 滑移边界参数对剪切应力有较大的影响, 滑移边界小的流体对壁面和边界层内流场施加了更大的剪切力; 增大雷诺数Re, 能增大内筒的壁面剪切力和扭矩; 内外筒半径比对流动结构也有较大的影响, 内筒半径固定, 增大外筒半径能够减少内筒的纵向剪切力.      

填隙二阶流体下圆球平行于壁平移的粘性阻力

离散元法是分析散体力学行为的数值方法。存在填隙流体时，颗粒之间或颗粒与壁之间产生的法向挤压力和切向阻力、阻力矩，是湿颗粒离散元法的理论基础。二阶流体是以微小偏离牛顿流体本构而考虑时间影响的一种流体。它具有常粘度，并且第一和第二法向应力差正比于剪切率的平方。根据Reynolds润滑理论，采用小参数法，导出了存在填隙二阶流体时，圆球沿平行于平壁缓慢移动时流体的速度场和压力方程，进而求出切向阻力和阻力矩的解析解。有趣的是在推导时所得的速度场和压力方程形式比牛顿流体要复杂得多，但最终结果表明圆球沿平行于平壁移动时因填隙二阶流体引起的切向阻力和阻力矩与牛顿流体时的结果相同。     

确定可靠性分析Weibull分布参数的概率加权矩法

用概率加权矩理论来确定描述可靠性分析中疲劳寿命Weibull分布的参数。推导了Weibull分布参数与样本概率加权矩的关系，给出了计算样本概率加权矩的MATLAB语言程序片段。根据疲劳寿命的样本值，计算各阶样本概率加权矩，然后估计疲劳寿命分布的参数。蒙特卡罗模拟的结果和应用实例表明，概率加权矩法计算结果精度较高，具有理论和实用价值。     

巷道内传质对通风阻力的影响

通风阻力的计算与测定是保障煤矿安全生产必须进行的日常工作。通常认为巷道内造成风流能量损失的途径主要是沿程摩擦损失、局部阻力损失、热阻力。传热导致流体物性参数会发生大变化,产生热阻力。巷道内不断有水蒸气和瓦斯等气体掺混进入主流空气中,形成复杂的对流传质现象,同样也会影响物性参数,从而产生附加通风阻力。这是一种被人们忽视的能量损失。本文对这种能量损失形式进行了探讨,提出了传质附加阻力概念,并运用CFD模拟技术研究了不同传热、不同散湿强度下通风阻力的变化规律。结果表明:当巷道风流速度低于0.5m/s时,传质使巷道阻力增加15%以上;当风速增大到2m/s以后,传质对巷道阻力的影响减弱,使巷道阻力约增大5%左右;当传质与传热同时存在时,传质会削弱热阻力,其综合效果不是简单的代数相加,传热与传质对阻力的影响是同一数量级。文中给出了传质附加阻力拟合计算式,可用于实际巷道通风设计和日常矿井通风能力核定。     

线性随机结构在随机激励下动力响应分析

利用虚拟激励法对随机结构正交展开理论进行扩展,并在Ritz向量子空间中对扩阶系统方程进行动力聚缩,提出了一类可以快速高效地进行线性随机结构复合随机振动分析的计算方法．算例分析表明,该法可以方便地分析随机结构在平稳或非平稳随机激励下的复合随机振动问题,且分析结果与 Monte Carlo模拟分析结果符合良好;与均值参数确定性结构传统随机振动分析计算结果相比,随机结构在相同随机激励下响应自谱密度曲线具有峰值降低、谱宽增大的特点．     

Large eddy simulation and spectrum analysis of flow around two square cylinders arranged side by side

Large eddy simulation cooperated with the second order full extension ETG ( Euler-Taylor-Galerkin ) finite element method was applied to simulate the flow around two square cylinders arranged side by side at a spacing ratio of 1.5. The second order full extension ETG finite element method was developed by Wang and He. By means of Taylor expansion of terms containing time derivative, time derivative is replaced by space derivative. The function of it is equal to introducing an artificial viscosity term. The streamlines of the flow at different moments were obtained. The time history of drag coefficient, lift coefficient and the streamwise velocity on the symmetrical points were presented. Furthermore, the symmetrical problem of the frequency spectrum of flow around two square cylinders arranged side by side were studied by using the spectral analysis technology. The data obtained at the initialstage are excluded in order to avoid the influence of initial condition on the results. The power spectrums of drag coefficient, lift coefficient, the streamwise velocity on the symmetrical points were analyzed respectively. The results show that although the time domain process of dynamic parameters is nonsymmetrical, the frequency domain process of them is symmetrical under the symmetrical boundary conditions.     

High‐fidelity aerodynamic shape optimization of modern transport wing using efficient hierarchical parameterization

Aerodynamic shape optimization technology is presented, using an efficient domain element parameterization approach. This provides a method that allows geometries to be parameterized at various levels, ranging from gross three‐dimensional planform alterations to detailed local surface changes. Design parameters control the domain element point locations and, through efficient global interpolation functions, deform both the surface geometry and corresponding computational fluid dynamics volume mesh, in a fast, high quality, and robust fashion. This results in total independence from the mesh type (structured or unstructured), and optimization independence from the flow‐solver is achieved by obtaining gradient information for an advanced gradient‐based optimizer by finite‐differences. Hence, the optimization tool can be used in conjunction with any flow‐solver and/or mesh generator. Results have been presented recently for two‐dimensional aerofoil cases, and shown impressive results; drag reductions of up to 45% were demonstrated using only 22 active design parameters. This paper presents the extension of these methods to three dimensions, with results for highly constrained optimization of a modern aircraft wing in transonic cruise. The optimization uses combined global and local parameters, giving 388 design variables, and produces a shock‐free geometry with an 18% reduction in drag, with the added advantage of significantly reduced root moments. Copyright © 2009 John Wiley & Sons, Ltd.     

Sensitivity and randomness in homogenization of periodic fiber-reinforced composites via the response function method

The main issue this paper addresses is the derivation and implementation of a general homogenization method, including the simultaneous determination of sensitivity gradients and probabilistic moments of the effective elasticity tensor. This is possible with an application of the perturbation method based on Taylor expansion and with the effective modules method. The computational procedure is implemented using plane strain analysis carried out with the finite element method (program MCCEFF) and the symbolic computations system MAPLE. The sensitivity gradients and probabilistic moments are commonly determined on the basis of partial derivatives for the homogenized elasticity tensor, calculated using the response function method with respect to some composite parameters. They are subjected separately to a normalization procedure (in deterministic analysis) and the relevant algebraic combinations (for the stochastic case). This enriched homogenization procedure is tested on a periodic fiber-reinforced two component composite, where the material parameters are taken as design variables and then, the input random quantities. The results of computational analysis are compared against the results of the central finite difference approach in the case of sensitivity gradients determination as well as the direct Monte-Carlo simulation approach. This numerical methodology may be further applied not only in the context of the homogenization method, but also to extend various discrete computational techniques, such as Boundary/Finite element and finite difference together with various meshless methods.     

薄板统计分析的随机边界元法

本文建立了分析含随机材料参数的薄板弯曲问题的随机边界元法。基于Ｔａｙｌｏｒ级数展开技术，分别得到了位移的均值和一阶偏差的边界积分方程，发现材料参数的随机性可作为一个等效的随机荷载处理，从而得到边界位移或边界力的均值和协方差，并进一步求出内点位移和力矩的均值和协方差，最后用本文方法计算了两个算例，并对结果进行了必要的分析。     

Effect of boundary representation on viscous,separated flows in a discontinuous-Galerkin Navier–Stokes solver

The effect of curved-boundary representation on the physics of the separated flow over a NACA 65(1)-412 airfoil is thoroughly investigated. A method is presented to approximate curved boundaries with a high-order discontinuous-Galerkin spectral element method for the solution of the Navier–Stokes equations. Multiblock quadrilateral element meshes are constructed with the grid generation software GridPro. The boundary of a NACA 65(1)-412 airfoil, defined by a cubic natural spline, is piecewise-approximated by isoparametric polynomial interpolants that represent the edges of boundary-fitted elements. Direct numerical simulation of the airfoil is performed on a coarse mesh and fine mesh with polynomial orders ranging from four to twelve. The accuracy of the curve fitting is investigated by comparing the flows computed on curved-sided meshes with those given by straight-sided meshes. Straight-sided meshes yield irregular wakes, whereas curved-sided meshes produce a regular Karman street wake. Straight-sided meshes also produce lower lift and higher viscous drag as compared with curved-sided meshes. When the mesh is refined by reducing the sizes of the elements, the lift decrease and viscous drag increase are less pronounced. The differences in the aerodynamic performance between the straight-sided meshes and the curved-sided meshes are concluded to be the result of artificial surface roughness introduced by the piecewise-linear boundary approximation provided by the straight-sided meshes.     

槽道湍流展向振荡电磁力减阻的机理研究

对槽道湍流的展向振荡电磁力控制进行了实验和数值研究. 实验通过PIV系统和浮动床阻力测试系统记录近壁区的条带变化和壁面阻力变化. 计算时, 利用谱方法直接模拟电磁力控制下的近壁流场. 实验和计算结果定性一致, 皆表明展向振荡电磁力可以减少壁面阻力, 并使条带倾斜. 计算结果还进一步揭示了电磁力减阻的机理. 电磁力诱导产生的流向涡与壁湍流的相互作用, 在近壁处形成负的脉动展向涡, 该涡将导致流向涡的倾斜和振荡, 从而抑制湍流, 减少壁面阻力.      

Mechanism of drag reduction by spanwise oscillating Lorentz force in turbulent channel flow

Numerical simulations and experimental research are both carried out to investigate the controlled effect of spanwise oscillating Lorentz force on a turbulent channel flow. The variations of the streaks and the skin friction drag are obtained through the PIV system and the drag measurement system, respectively. The flow field in the near-wall region is shown through direct numerical simulations utilizing spectral method. The experimental results are consistent with the numerical simulation results qualitatively, and both the results indicate that the streaks are tilted into the spanwise direction and the drag reduction utilizing spanwise oscillating Lorentz forces can be realized. The numerical simulation results reveal more detail of the drag reduction mechanism which can be explained, since the spanwise vorticity generated from the interaction between the induced Stokes layer and intrinsic turbulent flow in the near-wall region can make the longitudinal vortices tilt and oscillate, and leads to turbulence suppression and drag reduction.     

Drag force on a sphere in creeping motion throug a carreau model fluid

In the report of Chhabra and Uhlherr [1], a theoretical formulation for the drag force on a sphere in a fluid governed by the Carreau viscosity equation was presented. The formulation gave a good representation of the drag force when compared with experimental measurements, although the experiments were carried out on fluids exhibiting both shear-thinning and elastic behaviour. In the present work the ‘boundary element method’ has been used to obtain a fully numerical solution for comparison with the theoretical predictions. The two sets of results agree closely. We also present a simple empirical relation which gives a good first approximation to the drag force (～ 10%).     

用初参数法解C型波纹管在子午面内整体弯曲

波纹管是一种呈波纹状的管形薄壁结构,作为弹性敏感元件、密封元件和位移补偿器,在现代工业、国防和民用设备中的应用日益广泛.波纹管的计算,特别是在非轴对称受力和变形方面亟待完善.本文讨论了C型波纹管在子午面内整体弯曲的应力分布和端面角位移.从柔性圆环壳在子午面内整体弯曲的复变量方程出发,化边值问题为初值问题,用S.Gill法,得到了该波纹管的数值解,适用于0＜α＜1的范围.通过3个例子,将本文的数值解和已有的细环壳一般解(解析解)作了比较.结果表明,两者基本一致,并且,当环壳较细(α较小)时,两种计算结果更为接近.     

Viscoelastic flow past confined objects using a micro–macro approach

This paper is concerned with the numerical prediction of viscoelastic flow past a cylinder in a channel and a sphere in a cylinder using molecular-based models. The basis of the numerical method employed is a micro–macro model in which the polymer dynamics is described by the evolution of an ensemble of Brownian configuration fields. The spectral element method is used to discretize the equations in space. Comparisons are made between the macroscopic simulations based on the Oldroyd B constitutive model and microscopic simulations based on Hookean dumbbells, and excellent agreement is found. The micro–macro approach can be used to simulate models, such as the finitely extensible nonlinear elastic (FENE) dumbbell model, which do not possess a closed-form constitutive equation. Numerical simulations are performed for the FENE model. The influence of the model parameters on the flow is described and, in particular, the dependence of the drag as a function of the Weissenberg number.