高速列车服役模拟建模与计算方法研究

高速铁路的出现,使得轮轨交通技术达到更高的层次.速度提升不仅对列车的牵引动力与动力学性能提出更高要求.而且,列车与线路、气流等运行环境的耦合作用加剧,并直接影响到了列车的运行品质和安全性.在高速列车发展初期,研究关注的是如何保证高速列车能高速、平稳和安全运行.随着运行速度的提高,系统间耦合加强,服役模拟也越来越受到重视...     

铜、银及铂原子纳米团簇熔点随尺寸非单调变化的分子动力学模拟研究

采用分子动力学方法和原子嵌入法模型势模拟了铜原子、银原子及铂原子纳米团簇的熔化过程,研究了这些金属原子纳米团簇熔点与团簇尺寸的关系,发现各自在小于一定数目原子的金属纳米团簇的熔点不再随尺寸的变小而单调下降,通过对各种团簇溶化前后结构的比较研究,分析了导致这种现象的原因。     

非线性时滞动力系统的研究进展

具有时滞的动力系统广泛存在于各工程领域．本文从动力学角度对时滞动力系统的研究进展作一综述,内容包括时滞动力系统的特点、研究方法、动力学热点问题的研究进展等．由于时滞动力系统的演化趋势不仅依赖于系统的当前状态,还依赖于系统过去某一时刻或若干时刻的状态,其运动方程要用泛国微分方程来描述,解空间是无穷维的．即使系统中的时滞非常小,在许多情况下也不能忽略不计．对于非线性时滞常微分方程,目前的研究思路基本上与常微分方程系统理论相平行．主要研究方法可分为时域法和频域法,前者包括Ｔａｙｌｏｒ级数法,中心流形法,Ｐｏｉｎｃａｒｅ映射法等,后者包括Ｎｙｑｕｉｓｔ法等．目前对这类系统的动力学研究主要集中在稳定性、Ｈｏｐｆ分岔、混沌等方面．研究表明：时滞动力系统具有非常丰富和复杂的动力学行为,如单变量的一维非线性时滞动力系统可发生混沌现象,与用常微分方程描述的系统有本质性差别．另一方面,人们可巧妙地利用时滞来控制动力系统的行为,如时滞反馈控制是控制混饨的主要方法之一．最后,本文展望了存在的一些问题以及近期值得关注的研究．     

On damping coefficient due to phase transformation

The damping coefficient of capillary waves due to the evaporation-condensation process at the interface of the two phases of a fluid is evaluated. To highlight the mechanism of the effect of heat and mass transfer across the interface between regions of liquid and vapor, potential flow of incompressible fluids are assumed. Thus other mechanisms of damping are neglected. To fascilitate the analysis, the method of multiple-scale is employed in the analysis, even though the problem is linear.     

汽轮发电机定子结构的动力分析和模型试验研究

本文用有限元方法分析了两种常见汽轮发电机定子结构模型的动力特性,建立了固有特性对集中质量、连接刚度及基础刚度的灵敏度分析方法,对这两种结构模型进行了振动实验研究,实验结果与有限元分析结果吻合较好.从避开共振的角度出发,对两种结构作了评价.本文的研究方法及实验结果对汽轮发电机结构设计、减振、隔振和降噪等方面的工作具有一定的参考价值和实际工程意义.     

A terminally anchored polymer chain in shear flow: Self-consistent velocity and segment density profiles

The behavior of a terminally anchored freely-jointed bead-rod chain, subjected to solvent shear flow, was investigated via Brownian dynamics simulations. Previous calculations have been improved by computing the segment density and fluid velocity profiles self-consistently. The segment density distributions, components of the radius of gyration, and chain attachment shear and normal stresses were found to be sensitive to low values of shear rate. Additionally, it was found that the thickness of a model polymer layer was a strong function of the shear rate, and that the functional dependence on shear rate changed dramatically as the chain length increased. For the longest chains studied, the thickness of the model polymer layer first increased as the shear rate increased, passed through a maximum, and then decreased at high shear rates, in accordance with experimental results in theta solvents. These results suggest that a dilute or semi-dilute layer model may explain hydrodynamic behavior previously thought to be due to the entanglements that occur in dense surface bound polymer layers.Nomenclature a _i acceleration of bead i - b radius of the beads - d length of the rods connecting the chain beads - d _i vector from bead i to bead i + 1 - F _i external force applied to bead i - F _i ^b external force on bead i due to Brownian motion of surrounding fluid - F _i ^h external force on bead i due to viscous drag - F _i ^s external force on bead i due to surface interactions - f Stokes drag coefficient - Boltzmann's constant - L _h effective hydrodynamic thickness - m _i mass of bead i - N number of beads on a model chain - n number of chains anchored to the surface per unit surface area - P segment density distribution P pressure - Q flow in a tube with no surface bound polymer layer - Q _a flow in a tube with a surface bound polymer layer - R _g vector representation of the radius of gyration - R tube radius - r radial coordinate in the tube geometry - S _ij pair hydrodynamic interaction tensor for beads i and j - T _i internal chain force in rod i connecting beads i and i + 1 - T _X component of the surface attachment force in the direction of the fluid flow - T _y component of the surface attachment force perpendicular to the surface - T temperature - v _i velocity of the center of mass of bead i - V _if average fluid velocity at the location of bead i - v _if ⁰ fluid velocity in the absence of a polymer chain - v _if perturbation to the fluid velocity due to hydrodynamic interactions - V _b bead volume = 4 b ³/3 - scalar fluid speed in the axial direction down the tube - x axial coordinate in the tube geometry Greek symbols _w apparent shear rate - fluid viscosity - polymer layer permeability - volume fraction of space occupied by chain beads - (_w)_a chain attachment stress perpendicular to the surface - (_w)_a chain attachment stress in the plane of the surface and in the direction of fluid flow     

多体系统动力学微分／代数方程组数值方法

多体系统动力学微分／代数混合方程组又称Ｅｕｌｅｒ－ｌａｇｒａｎｇｅ方程,是近十年来动力学和计算数学领域研究的热点之一．本文介绍这两个领域中引入的传统的数值积分方法与新的理论．      

Chaotic motions of a rigid rotor in short journal bearings

In the present paper the conditions that give rise to chaotic motions in a rigid rotor on short journal bearings are investigated and determined. A suitable symmetry was given to the rotor, to the supporting system, to the acting system of forces and to the system of initial conditions, in order to restrict the motions of the rotor to translatory whirl. For an assigned distance between the supports, the ratio between the transverse and the polar mass moments of the rotor was selected conveniently small, with the aim of avoiding conical instability. Since the theoretical analysis of a system's chaotic motions can only be carried out by means of numerical investigation, the procedure here adopted by the authors consists of numerical integration of the rotor's equations of motion, with trial and error regarding the three parameters that characterise the theoretical model of the system: m, the half non-dimensional mass of the rotor, , the modified Sommerfeld number relating to the lubricated bearings, and , the dimensionless value of rotor unbalance. In the rotor's equations of motion, the forces due to the lubricating film are written under the assumption of isothermal and laminar flow in short bearings. The number of numerical trials needed to find the system's chaotic responses has been greatly reduced by recognition of the fact that chaotic motions become possible when the value of the dimensionless static eccentricity % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbnL2yY9% 2CVzgDGmvyUnhitvMCPzgarmWu51MyVXgaruWqVvNCPvMCG4uz3bqe% fqvATv2CG4uz3bIuV1wyUbqee0evGueE0jxyaibaieYlf9irVeeu0d% Xdh9vqqj-hEeeu0xXdbba9frFf0-OqFfea0dXdd9vqaq-JfrVkFHe9% pgea0dXdar-Jb9hs0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaaca% qabeaadaabauaaaOqaaiabew7aLnaaBaaaleaacaWGZbaabeaaaaa!4046!\[\varepsilon _s \] is greater than 0.4. In these conditions, non-periodic motions can be obtained even when rotor unbalance values are not particularly high (=0.05), whereas higher values (>0.4) make the rotor motion periodic and synchronous with the driving rotation. The present investigation has also identified the route that leads an assigned rotor to chaos when its angular speed is varied with prefixed values of the dimensionless unbalance . The theoretical results obtained have then been compared with experimental data. Both the theoretical and the experimental data have pointed out that in the circumstances investigated chaotic motions deserve more attention, from a technical point of view, than is normally ascribed to behaviours of this sort. This is mainly because such behaviours are usually considered of scarce practical significance owing to the typically bounded nature of chaotic evolution. The present analysis has shown that when the rotor exhibits chaotic motions, the centres of the journals describe orbits that alternate between small and large in an unpredictable and disordered manner. In these conditions the thickness of the lubricating film can assume values that are extremely low and such as to compromise the efficiency of the bearings, whereas the rotor is affected by inertia forces that are so high as to determine severe vibrations of the supports.Nomenclature C radial clearance of bearing (m) - D diameter of bearing (m) - e dimensional eccentricity of journal (m) - e _s value of e corresponding to the static position of the journal - E dimensional static unbalance of rotor (m) - f _x, f _y =F _x/(P), F _y/(P): non-dimensional components of fluid film force - F _x, F _y dimensional components of fluid film force (N) - g acceleration of gravity (m/s²) - L axial length of bearing (m) - m % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbnL2yY9% 2CVzgDGmvyUnhitvMCPzgarmWu51MyVXgaruWqVvNCPvMCG4uz3bqe% fqvATv2CG4uz3bIuV1wyUbqee0evGueE0jxyaibaieYlf9irVeeu0d% Xdh9vqqj-hEeeu0xXdbba9frFf0-OqFfea0dXdd9vqaq-JfrVkFHe9% pgea0dXdar-Jb9hs0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaaca% qabeaadaabauaaaOqaaiabg2da9maalaaabaGaeqyYdC3aaWbaaSqa% beaacaaIYaaaaaGcbaGaeqyYdC3aa0baaSqaaGabciaa-bdaaeaaca% WFYaaaaaaakiabg2da9maalaaabaGaeqyYdC3aaWbaaSqabeaacaaI% YaaaaOGaam4qaaqaaiabeo8aZjaadEgaaaaaaa!4C14!\[ = \frac{{\omega ^2 }}{{\omega _0^2 }} = \frac{{\omega ^2 C}}{{\sigma g}}\]: half non-dimensional mass of rotor - M half mass of rotor (kg) - n angular speed of rotor (in r.p.m.=60/2) - t time     

多体系统动力学微分／代数方程组数值方法

多体系统动力学微分／代数混合方程组又称Ｅｕｌｅｒ－ｌａｇｒａｎｇｅ方程，是近十年来动力学和计算数学领域研究的热点之一．本文介绍这两个领域中引入的传统的数值积分方法与新的理论．     

具有局部非线性动力系统周期解及稳定性方法

对于具有局部非线性的多自由度动力系统，提出一种分析周期解的稳定性及其分岔的方法该方法基于模态综合技术，将线性自由度转换到模态空间中，并对其进行缩减，而非线性自由度仍保留在物理空间中在分析缩减后系统的动力特性时，基于Ｎｅｗｍａｒｋ法的预估－校正－局部迭代的求解方法，与Ｐｏｉｎｃａｒé映射法相结合，推导出一种确定周期解，并使用Ｆｌｏｑｕｅｔ乘子判定其稳定性及分岔的方法