柱塞泵多目标优化设计及CFD仿真分析

液压泵噪声是液压系统的主要噪声源,针对轴向柱塞泵的流致振动噪声,提出一种改善泵配流特性的设计方案。首先,根据柱塞泵的工作原理对柱塞腔压力特性和泵出口流量特性准确建模并求解。通过分析压力冲击和流量脉动对错配角(φ0)的响应,得φ0=4°为佳。利用一种多目标遗传算法(NSGA-II),以减小压力超调量和流量脉动率为目标,对三角槽结构进行了优化;并获得该多目标优化问题的Pareto最优解集,通过对最优解集的分析知,深度角θ1=16°且宽度角θ2=85°时较为理想。最后,为了验证模型的正确性,建立流体域计算流体动力学(CFD)模型,对比两种模型计算结果发现吻合较好,能够相互验证。利用CFD分析结果可视化的特点,从柱塞泵流场的角度,进一步分析了泵压力冲击以及流量脉动产生的原因。     

基于转捩/尺度适应模型与FW-H声学方程的气动噪声数值模拟研究

采用计算流体力学(CFD)同"声类比"相结合的方法进行噪声模拟,利用CFD数值模拟 MD30P30N多段翼型失速攻角附近流场以及其气动特性来校核近场精度,进一步通过结合可穿透数据面的FW-H声学方程进行气动噪声分析。为了准确捕捉近场流场信息,为噪声预测提供可靠的声源精度,本文基于k-ω剪应力输运(SST)湍流模型,建立了尺度适应分离流(SAS)模型,并采用γ-Re_θt转捩模型耦合k-ω SST湍流模型建立了边界层转捩数值模拟技术;充分利用尺度适应模型在边界层表现为雷诺平均(RANS)方法这一特点,将γ-Re_θt转捩模型与SAS模型结合,建立针对包含转捩、分离现象复杂流场的数值模拟技术。文中以RANS为控制方程,分别采用全湍流k-ω SST模型、γ-Re_θtSAS转捩/分离流模型对多段翼进行近场数值模拟计算,结合可穿透数据面的FW-H声学方程进行气动噪声预测,在分析其对流场及气动噪声影响的基础上,得出了几点结论。     

On near-wall hot-wire measurements

A specially constructed hot-wire probe was used to obtain very near-wall velocity measurements in both a fully developed turbulent channel flow and flat plate boundary layer flow. The near-wall hot-wire probe, having been calibrated in a specially constructed laminar flow calibration rig, was used to measure the mean streamwise velocity profile, distributions of streamwise and spanwise intensities of turbulence and turbulence kinetic energy k in the viscous sublayer and beyond; these distributions compare very favorably with available DNS results obtained for channel flow. While low Reynolds number effects were clearly evident for the channel flow, these effects are much less distinct for the boundary layer flow. By assuming the dissipating range of eddy sizes to be statistically isotropic and the validity of Taylor's hypothesis, the dissipation rate ɛ _iso in the very near-wall viscous sublayer region and beyond was determined for both the channel and boundary layer flows. It was found that if the convective velocity U _c in Taylor's hypothesis was assumed to be equal to the mean velocity Uˉ at the point of measurement, the value of (ɛ⁺ _iso)₁ thus obtained agrees well with that of (ɛ ⁺)_DNS for y ⁺ ≥ 80 for channel flow; this suggests the validity of assuming U _c=Uˉ and local isotropy for large values of y ⁺. However, if U _c was assumed to be 10.6u _τ , the value of (ɛ⁺ _iso)₂ thus obtained was found to compare reasonably well with the distribution of (ɛ⁺ _iso)_DNS for y ⁺≤ 15. Received: 31 May 1999/Accepted: 20 December 1999     

Experiments on the effect of acceleration on the drag of tapwater bubbles

Three-dimensional (3D) trajectories of spherical air bubbles passing through a converging part of a rectangular channel have been measured. Bubble diameters, d _b, were less than 1 mm and the Reynolds numbers, Re _b, for stagnant tapwater and for mean liquid velocity vˉ _L=0.25 m/s were in about same range. Received: 15 January 2001 / Accepted: 12 June 2001     

近源堆沉积土石体细观结构与抗剪强度指标关联分析

细观结构是认知土石体力学行为本质的关键科学问题。本文通过引入物理学和数学等方法，借助颗粒物质力学理论，从几何排列与接触力的空间分布来定量刻画水平固结与山前坡地堆积两种典型环境下土石体的细观结构特征，并建立其与抗剪强度指标的关联。研究表明，（1）两种环境的土石体在细观结构上存在较大差异。在几何排列上，水平固结环境下的土石体具有长程无序和短程有序的特点，坡地堆沉积环境下的土石体表现出了无序的无定形结构；在接触力与单位接触向量的空间分布上，两者较为相似，绝大多数接触力以小于均值接触力的形式存在，其概率密度曲线P（f）呈幂函数衰减；90%以上接触方位角集中在40°~160°和220°~340°范围内。（2）基于径向分布函数、接触力概率密度和单位接触向量分别定义细观结构的特征量K_a，K_s和K_o，发现三个特征量的增大对内摩擦角呈线性促进作用，对黏聚力呈非线性削弱作用。     

On the free energy and variational theorem of elastic problem in thermal shock

In some investigations on variational principle for coupled thermoelastic problems, the free energy Φ(e_ij,θ) ,where the state variables are elastic strain e_ij and temperature increment θ, is expressed as Φ(e_ij,θ)=λ/2e_kke_ij=ue_k1e_k1-γe_kkθ-c/2 p θ²/T₀(0.1) This expression is employed only under the condition of |θ|≤T₀(absolute temperature of reference) But the value of temperature increment is great, even greater than T₀ in thermal shock. And the material properties (λ ,μ ,ν ,c , etc.) will not remain constant, they vary with θ. The expression of free energy for this condition.is derived in this paper. Equation (0.1) is its special case.Euler’s equations will be nonlinear while this expression of free energy has been introduced into variational theorem. In order to linearise, the time interval of thermal shock is divided into a number of time elements Δt_k, (Δt_k=t_k-t_k-1,k=1,2…,n), which are so small that the temperature increment θ_k within it is very small, too. Thus, the material properties may be defined by temperature field T_k-1=T(x₁,x₂,x₃,t_k-1) at instant t_k-1 , and the free energy Φ_k expressed by eg. (0.1) may be employed in element Δt_k.Hence the variational theorem will be expressed partly and approximately.     

Rheological aspects of dense lignite-water suspensions; time dependence, preshear and solids loading effects

An investigation of the rheological properties of dense lignite-water suspensions is reported here in order to evaluate the parameters that affect the rheology of these suspensions. Different types of particle size distributions were achieved via wet grinding with particle sizes down to 1 μm. Polyelectrolytes and surfactants were used to improve the stability and rheology of suspensions. The achieved solid volume fraction, φ, with acceptable rheological behavior was 0.45, while the ratio of φ/φ _m was close to 0.85; φ _m is the maximum solid volume fraction. Various types of flow tests were performed on lignite-water suspensions with various preshear times and stress levels. Typical shear rates varied in the range between 10^–2 to 10² s^–1. The low shear behavior was found to be quite different from the high shear behavior for the most concentrated suspensions, with a plateau value at low shear rates indicating the appearance of a yield stress. These plateau values, however, depend on shear history, which is responsible for the development of different structures in the sample. The viscosity curves corresponding to the ascending and descending parts of the flow curve were found to be different; these flow curve parts can be described either by well established models or by modified ones. The non-Newtonian time dependent behavior of the lignite-water suspensions is attributed to the high value of the ratio φ/φ _m , the polydispersity of the particle size distribution, and the non-spherical shape of the lignite particles.     

Studying the stability of equilibrium solutions in a planar circular restricted four-body problem

The Newtonian circular restricted four-body problem is considered. We obtain nonlinear algebraic equations determining equilibrium solutions in the rotating frame and find six possible equilibrium configurations of the system. Studying the stability of equilibrium solutions, we prove that the radial equilibrium solutions are unstable, while the bisector equilibrium solutions are stable in Lyapunov’s sense if the mass parameter satisfies the conditions μ ∈ (0, μ₀, where μ₀ is a sufficiently small number, and μ ≠ μ_j, j = 1, 2, 3. We also prove that, for μ = μ₁ and μ = μ₃, the resonance conditions of the third order and the fourth order, respectively, are satisfied and, for these values of μ, the bisector equilibrium solutions are unstable and stable in Lyapunov’s sense, respectively. All symbolic and numerical calculations are done with the Mathematica computer algebra system. Published in Neliniini Kolyvannya, Vol. 10, No. 1, pp. 66–82, January–March, 2007.     

The analytical solution for Helmholtz boundary problem in non-horizontally stratified domains

There are N domains D_j(j=0,1,...,N-1) of different physical parameters in the whole space and their interfaces S_j,i+1 are non-horizontally smooth curved surfaces. The following boundary problem is called Hclinholiz boundary problem:V²H⁽¹⁾+K_jH⁽¹⁾=0 (j=0,1,…,N-1)(H⁽⁰⁾-H⁽¹⁾)=δ(S) (δ(S):generalized function)(H¹-Hⁱ⁺¹)=0 (j=0,1,…,N-2)The analytical solution of the above problem is given in this paper.     

Havelock型格林函数振荡项数值积分的稳定性研究

Havelock型格林函数的传播项被积函数是高频振荡且奇异的复变函数,文献[4]引入变量代换获得了一种兼具积分效率和精度的积分方法, 本文研究了该方法的积分稳定性,发现该方法仍存在如下的积分困难: (1) θ=γ时复函数中分母为零引起的计算溢出;(2) θ=π/2是复函数在y及z方向偏导数的无穷间断点;(3) 场点与源点横坐标相同时伪奇异点变为真奇点。针对这些积分困难,采用极限公式计算θ=γ处复函数的值避免计算溢出;在保证积分精度的前提下采用截断方法略去θ=π/2邻近区域的积分消除无穷间断处的奇异;针对(3)采用分区法处理以避开原被积函数的高频振荡,并消除奇异性。伪奇异性存在的条件是场点必须在点源传播波的传播范围内,伪奇异点最多为2个。     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).