车桥耦合系统随机振动的时域显式解法

在桥面和轨道随机不平顺作用下,车桥耦合系统振动是一个典型的非平稳随机振动问题.笔者分别建立表征物理演变机制的车辆系统和桥梁系统的动力响应显式表达式,然后利用车桥之间的运动相容条件,建立车桥之间接触力关于桥面不平顺的显式表达式.在此基础上,即可直接利用统计矩运算法则,获得车桥接触力的统计矩演化规律,并进一步计算车辆系统和桥梁系统关键响应的演变统计矩.此外,也可以基于车桥接触力关于桥面不平顺的显式表达式,高效地进行随机模拟（即Monte Carlo模拟, MCS）,以获得车桥耦合系统关键响应的演变统计矩及其他统计信息.在上述过程中,由于实现了车桥耦合系统物理演变机制和概率演化规律的相对分离,在响应统计矩计算中,无需反复求解车桥耦合系统的运动微分方程,且可以仅针对车桥接触力及其他所关注的关键响应开展降维计算,大幅提高了车桥耦合系统随机振动的计算效率.数值算例表明,所提出的方法具有理想的计算精度和计算效率.     

A stochastic dynamic model of train-track-bridge coupled system based on probability density evolution method

An innovative stochastic dynamic model of a 3D train-track-bridge coupled system (TTBS) with refined wheel/rail interaction is established for a high-speed railway based on the random theory of probability density evolution method (PDEM). The multi-coupling effect of excitations can be simultaneously input into the new model, e.g. random track irregularity, random vehicle loads, stochastic system parameters, et al. Moreover, a new approach, named “Number theoretic method of multi-target probability functions” (NTM-mp), is developed to obtain the discrete point sets of multidimensional random parameters in hypercube space, aims to solve the point design of system uncertainty. The stochastic harmonic function (SHF) is applied to generate representative random track irregularity samples. The results of TTBS got by PDEM are verified with several typical case studies for its efficiency and reliability, which are the deterministic results in the representative publication, the Monte Carlo method (MCM) results, and the field testing results on the high-speed railway. At last, a typical case study of TTBS on a high-speed railway is presented for numerical analysis. Discussions and significant conclusions on the random dynamic responses are presented.     

A coupled model for train-track-bridge stochastic analysis with consideration of spatial variation and temporal evolution

Due to random characteristics of system parameters and excitations, the dynamic assessment and prediction for the train-track-bridge interaction systems become rather complex issues needing to be addressed, especially considering the longitudinal inhomogeneity and uncertainty of dynamic properties in physics and correspondingly their temporal evolutions. In this paper, a temporal-spatial coupled model is developed to fully deal with the deterministically/non-deterministically computational and analytical matters in the train-track-bridge interactions with a novelty, where a train-track-bridge interaction model is newly developed by effectively coupling the three-dimensional nonlinear wheel-rail contact model and the finite element theory, moreover, the Monte-Carlo method (MCM) and Karhunen–Loève expansion (KLE) are effectively united to model the random field of track-bridge systems, and a spectral evolution method accompanied by a track irregularity probabilistic model are introduced to select the most representative track irregularity sets and to characterize their random evolutions in temporal dimension. In terms of random vibration analysis, the high-efficiency and effectiveness of this developed model is validated by comparing to a robust method, i.e., MCM. Apart from validations, multi-applications of the temporal-spatial coupled model from aspects of deterministic computation, random vibration, resonant analysis and long-term dynamic prediction, etc., have been fully presented to illustrate the universality of the proposed model.     

二系悬挂条件下的车-路垂向耦合系统的动力模型

在车辆的走行过程中,上部与下部是相互作用和影响的．因此,轨道交通问题实际上就是线路下部结构和车辆系统的体系匹配问题．在一系悬挂条件下的车-路系统耦合动力分析模型基础上,考虑了包含转向架在内的车辆的实际构成和轨下基础包括路基和地基的参振特性,利用轨道维护标准模拟行走不平顺激励,通过位移相容条件,从理论上研究车辆-轨道-路基体系的动力相互作用,建立了二系悬挂条件下的车辆-轨道-路基系统的垂向动力分析模型．为高速铁路路基的动力特性分析和设计提供参考．     

A novel dynamics model for railway ballastless track with medium-thick slabs

Motivated by the requirements for elaborated slab ballastless track dynamics analysis in practical engineering application, a novel dynamic model for the railway ballastless tracks with medium-thick slabs is proposed in this work based on the Reissner–Mindlin plate theory, and it is implemented into the coupled dynamics analysis of a vehicle and the ballastless track. First, an efficient and easily programmable computational algorithm is adopted to solve the transverse deflection of the Reissner–Mindlin plate, in which the displacements and shear strains are chosen as the independent variables and subsequently constructed by spline functions, resulting in no shear-locking effect. The involved partial differential equations are transformed into ordinary ones by using the energy variation principle. Further, a mathematical model for the ballastless track dynamics analysis is established, which can consider the effects of the shear deformation and moment of inertia involved in the medium-thick track slab. Experimental verification and comparative analysis with other models demonstrate the accuracy and efficiency of the proposed model. Finally, a spatially coupled dynamics model of a vehicle and the ballastless track is developed, and it is efficiently solved by using the hybrid explicit-implicit time integration method. Compared with the widely used modelling the track slab by elastic thin plate, the reliability and advantages of the proposed vehicle-slab track coupled dynamics model are demonstrated.     

轨道结构随机场模型与车辆-轨道耦合随机动力分析

将轨道结构视为一个参数随机系统,提出并建立了轨道结构的随机场模型.利用车辆-轨道耦合动力学的基本方法,将轨道系统有限单元模型与多刚体车辆模型相结合,建立了考虑铁路线路参数空-时随机变化的车辆-轨道动力计算模型.算例表明:所提出的方法较为可靠且高效;线路参数随机性对车辆-轨道系统的动力响应有明显的影响,随线路参数离散程度的增加,可能造成行车不安全、轨道损伤加剧等一些问题.     

Dual track and segmented single track bidirectional loop guidepath layout for AGV systems

In this study, we develop integer programming models for simultaneous design of the guidepath track and location of the pickup and dropoff stations for automated guided vehicle systems. Two alternative bidirectional topologies are compared: a dual track loop, and a single track loop partitioned into non-overlapping segments, each served by a single vehicle. We show that the empty vehicle trips have a significant impact on the quality of design in both topologies. More importantly, we identify the difference between empty flow node balancing for a dual track loop with the empty flow edge balancing on a single track loop subject to segmentation. Our work presents an economic decision making model for the trade-off analysis between the costs associated with the additional shop floor space of a dual track loop and the increased vehicle fleet size of a segmented single track loop.     

Study on numerical simulation methods of the railway vehicle-track dynamic interaction

A very efficient numerical simulation method of the railway vehicle–track dynamic interaction is described. When a vehicle runs at high speed on the railway track, contact forces between a wheel and a rail vary dynamically due to the profile irregularities existing on the surface of the rail. A large variation of contact forces causes undesired deteriorations of a track and its substructures. Therefore these dynamic contact forces are of main concern of the railway engineers. However it is very difficult to measure such dynamic contact forces directly. So it is important to develop an appropriate numerical simulation model and identify structural factors having a large influence on the variation of contact forces. When a contact force is expressed by the linearized Hertzian contact spring model, the equation of motions of the system is expressed as a second–order linear time–variant differential equation which has a time–dependent stiffness coefficient. Applying a well–known Newmark direct integration method, a numerical simulation is reduced to solving iteratively a time–variant, large–scale sparse, symmetric positive–definite linear system. In this study, by defining a special vector named a contact point one, it is shown that this time–variant stiffness coefficient can be expressed simply as a product of the contact point vector and its transpose and so the Sherman–Morrison–Woodbury formula applied for updating the inverse of the coefficient matrix. As a result, the execution of numerical simulation can be carried out very efficiently. A comparison of the computational time is given. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Transversal vibration analysis of the upper span of nonlinear closed-loop track systems

In order to analyze the effect of track pre-tension on the transversal vibration of the upper track, a multibody dynamics model of a tracked vehicle needs to be established. In view of the complex structure of tracked vehicles, the computational efficiency of conventional methods cannot meet engineering design requirements. In this paper, the Riccati Transfer Matrix Method for nonlinear closed-loop multibody systems is proposed and applied to the modeling and simulation of tracked vehicle systems. It avoids the overall dynamic equations of the system and achieves efficient computation with a smaller matrix scale. By selecting accelerations and forces as elements of the state vectors instead of coordinates and forces, linearization is avoided. The track chain is cut in a revolute joint. Its effect on the closed-loop system is represented by a set of constraint equations and internal constraint forces. An improved Riccati transformation for closed-loop systems is proposed to connect the internal force and acceleration parts in the state vector with the system constraint internal forces. The track pre-tension can be adjusted by moving the position of the idler relative to the chassis. The root mean square of vertical acceleration of mid-span with respect to the chassis is selected to evaluate the transversal vibration of the upper track. Simulation of a tracked vehicle driving on different road conditions is carried out with different idler displacements. Finally, a reasonable adjustment strategy for the track tensioner is proposed corresponding to different running conditions.     

车辆-轨道垂向耦合系统求解过程的改进算法

考虑到轨道结构长度随系统响应持时的增加而增长，提出了一种改进的车辆 轨道垂向耦合系统的动力响应求解算法.该算法事先选定某一定长度的轨道结构，并获得该轨道结构的质量矩阵、阻尼矩阵和刚度矩阵；通过在求解过程中不断地对车辆子系统定位，判断是否需要对车辆子系统的位置和轨道结构的响应矩阵进行调整，以此来达到仅增加系统响应持时而不增加轨道结构长度的目的.算例表明:该改进加快算法是精确、高效的，不仅可以真实地模拟车辆在轨道上的前进运行状态，而且可以保证轨道子系统的轨道单元数量不随系统响应持时的增加而增长，这为快速求解车辆 轨道垂向耦合系统提供了一种有效的计算方法.     

Cycle-transitive comparison of independent random variables

The discrete dice model, previously introduced by the present authors, essentially amounts to the pairwise comparison of a collection of independent discrete random variables that are uniformly distributed on finite integer multisets. This pairwise comparison results in a probabilistic relation that exhibits a particular type of transitivity, called dice-transitivity. In this paper, the discrete dice model is generalized with the purpose of pairwisely comparing independent discrete or continuous random variables with arbitrary probability distributions. It is shown that the probabilistic relation generated by a collection of arbitrary independent random variables is still dice-transitive. Interestingly, this probabilistic relation can be seen as a graded alternative to the concept of stochastic dominance. Furthermore, when the marginal distributions of the random variables belong to the same parametric family of distributions, the probabilistic relation exhibits interesting types of isostochastic transitivity, such as multiplicative transitivity. Finally, the probabilistic relation generated by a collection of independent normal random variables is proven to be moderately stochastic transitive.     

Reliable control for nonlinear systems with stochastic actuators fault and random delays through a T-S fuzzy model approach

This paper considers the reliable control design for T-S fuzzy systems with probabilistic actuators faults and random time-varying delays. The faults of each actuator occurs randomly and its failure rates are governed by a set of unrelated random variables satisfying certain probabilistic distribution. In terms of the probabilistic failures of each actuator and time-varying random delays, new fault model is proposed. Based on the new fuzzy model, reliable controller is designed and sufficient conditions for the exponentially mean square stability (EMSS) of T-S fuzzy systems are derived by using Lyapunov functional method and linear matrix inequality (LMI) technique. It should be noted that the obtained criteria depend on not only the size of the delay, but also the probability distribution of it. Finally, a numerical example is given to show the effectiveness of the proposed method.     

Three-dimensional vehicle-ballasted track-subgrade interaction: Model construction and numerical analysis

In this work a three-dimensional vehicle-ballasted track-subgrade interaction model is developed, where the vehicle is modeled as a multi-rigid-body system, the track-subgrade interaction is modelled by finite element method (FEM) with the rail modelled as beam elements, the sleeper and the subgrade layers modelled as solid elements, and the vehicle and the track-subgrade system are coupled by unified matrix formulations and solved simultaneously and time-dependently. As the modelling advancement, the versatility for FEM construction of track-subgrade systems has been promoted, where the elemental scale can be arbitrarily chosen without consideration of the node-to-node matching principle as the conventional methods; besides the vehicle and the track-subgrade systems have been coupled effectively by wheel-rail contact models. The applicability of various wheel-rail contact models has been discussed, and to reveal the effectiveness of this model in solving engineering problems such as the soil elasticity unevenness and the contact break of “hanging sleepers”, numerical examples have also been presented with referencable conclusions.     

Convexity and optimization with copulæ structured probabilistic constraints

Probability constraints play a key role in optimization problems involving uncertainties. These constraints request that an inequality system depending on a random vector has to be satisfied with a high enough probability. In specific settings, copulæ can be used to model the probabilistic constraints with uncertainty on the left-hand side. In this paper, we provide eventual convexity results for the feasible set of decisions under local generalized concavity properties of the constraint mappings and involved copulæ. The results cover all Archimedean copulæ. We consider probabilistic constraints wherein the decision and random vector are separated, i.e. left/right-hand side uncertainty. In order to solve the underlying optimization problem, we propose and analyse convergence of a regularized supporting hyperplane method: a stabilized variant of generalized Benders decomposition. The algorithm is tested on a large set of instances involving several copulæ among which the Gaussian copula. A Numerical comparison with a (pure) supporting hyperplane algorithm and a general purpose solver for non-linear optimization is also presented.     

A Poisson random field model for intermittent phenomena with application to laminar-turbulent transition and material microstructure

The alternation of a physical system between two phases or states is referred to as intermittency. Examples of intermittent phenomena abound in applications and include the transition from laminar to turbulent flow over a flight vehicle and the presence of imperfections within material microstructure. It is shown that intermittent phenomena of this type can be modeled by two-state random fields with piecewise constant samples; we refer to the states of the random field as “off” and “on” or, equivalently, 0 and 1. These random fields can be calibrated to the available information, which consists of: (1) the marginal probability that the state of the system is “on”; and (2) the average number of fluctuations between states that occur within a bounded region. The proposed model is defined by a sequence of pulses of prescribed shape and unit magnitude, located at random (Poisson) points within a bounded domain. Properties of the model are discussed, and simple algorithms to generate samples of the random field are provided. Various applications are considered, including voids within material microstructure and the random vibration of a flight vehicle subjected to a transition from laminar to turbulent flow over its surface.     

Probabilistic representation of solutions of hydrodynamic equations

It is well known that the relationships between the theory of partial differential equations (PDE) and probability theory are very deep and allow one to derive new nontrivial results concerning the behavior of random processes and properties of the solutions of boundary-value problems of PDE. It is natural to make an attempt to construct a probabilistic representation of the solution of a boundary problem for the Navier-Stokes system and other hydrodynamic equations. We construct such a representation in terms of a diffusion process. Bibliography: 16 titles. To O. A. Ladyzhenskaya Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 249, 1997, pp. 77–101. Translated by Ya. I. Belopolskaya.     

Stochastically forced cardiac bidomain model

The bidomain system of degenerate reaction–diffusion equations is a well-established spatial model of electrical activity in cardiac tissue, with “reaction” linked to the cellular action potential and “diffusion” representing current flow between cells. The purpose of this paper is to introduce a “stochastically forced” version of the bidomain model that accounts for various random effects. We establish the existence of martingale (probabilistic weak) solutions to the stochastic bidomain model. The result is proved by means of an auxiliary nondegenerate system and the Faedo–Galerkin method. To prove convergence of the approximate solutions, we use the stochastic compactness method and Skorokhod–Jakubowski a.s. representations. Finally, via a pathwise uniqueness result, we conclude that the martingale solutions are pathwise (i.e., probabilistic strong) solutions.     

Adaptive multi-element polynomial chaos with discrete measure: Algorithms and application to SPDEs

We develop a multi-element probabilistic collocation method (ME-PCM) for arbitrary discrete probability measures with finite moments and apply it to solve partial differential equations with random parameters. The method is based on numerical construction of orthogonal polynomial bases in terms of a discrete probability measure. To this end, we compare the accuracy and efficiency of five different constructions. We develop an adaptive procedure for decomposition of the parametric space using the local variance criterion. We then couple the ME-PCM with sparse grids to study the Korteweg–de Vries (KdV) equation subject to random excitation, where the random parameters are associated with either a discrete or a continuous probability measure. Numerical experiments demonstrate that the proposed algorithms lead to high accuracy and efficiency for hybrid (discrete–continuous) random inputs.     

Comportamiento dinámico de viaductos cortos considerando la interacción vehículo-vía-estructura-suelo

The dynamic vehicle-track-bridge-soil interaction is studied on high speed lines. The analysis is carried out using a general and fully three dimensional multi-body-finite element-boundary element model, formulated in the time domain to predict vibrations due to the train passage over the bridge. The vehicle is modelled as a multi-body system, the track and the bridge are modelled using finite elements and the soil is considered as a homogeneous half-space by the boundary element method. Usually, moving force model and moving mass model are employed to study the dynamic response of bridges. In this work, the multi-body system allows one to consider the quasi-static and dynamic excitation mechanisms. Soil-structure interaction is taken into account on the dynamic structure behaviour on simply-supported short span bridges. The influence of soil-structure interaction is analysed in both resonant and non-resonant regimes.     

随机环境中具有随机控制函数的受控分枝过程

对随机环境中具有随机控制函数的受控分枝过程进行了更为详尽的概率描述和直观解释;证明了此过程是时齐马氏链和随机环境中的马氏链,并对其概率母函数及矩量进行了讨论。