Xiang Liu

Flow-induced vibration of high-speed pantographs becomes increasingly important as train speed approaches and exceeds 400 km/h. This paper develops an efficient computational framework that couples component-resolved unsteady aerodynamic loads from improved delayed detached eddy simulation (IDDES) with a dynamic-stiffness-method (DSM) flexible model of a high-speed pantograph. Two operating orientations, namely, knuckle-downstream and knuckle-upstream, are compared at 400 km/h, and the more unfavorable knuckle-upstream orientation is further investigated over 400 km/h to 600 km/h. The DSM model contains 49 beam elements and 42 nodes and shows good agreement with a refined three-dimensional solid-element finite element model in the low-order frequency range. For a 2 s transient analysis, the proposed model predicts the panhead displacement response with peak errors below 5% relative to the finite element model while reducing the computational time from 53 min 22 s to 35 s on the same platform. The results show that vertical vibration dominates the structural response, with the panhead peak vertical displacement reaching about 20 mm in the studied 400 km/h open-line case. Frequency-domain inspection of the panhead aerodynamic lift and vertical displacement shows that broadband aerodynamic excitation mainly activates the low-order structural modes, with a low-frequency aerodynamic component around 3 Hz to 4 Hz and additional energy mainly over the 20 Hz to 30 Hz range. The knuckle-upstream orientation increases the standard deviation of the equivalent contact-force response by 46% compared with the knuckle-downstream orientation at 400 km/h. For the knuckle-upstream orientation, increasing speed from 400 km/h to 600 km/h raises the standard deviation by 189%. The proposed framework provides an efficient tool for rapid comparative evaluation of pantograph flow-induced vibration under multiple operating conditions.

An efficient modal and critical speed analysis method is proposed for railway catenary systems based on analytical dynamic stiffness (DS) method empowered by substructure technique and the extended Wittrick-Williams algorithm. First, a systematic and versatile DS modeling framework combined with substructure technique is proposed for refined modeling of different types of catenaries. Consequently, the broadband vibrations are described exactly by very few degrees of freedom. Then, the Wittrick–Williams algorithm is extended for the condensed DS model of a multi-span catenary to compute exact natural frequencies. Also, the DS model for a single-span is also established following Bloch's theorem. Exact dispersion curves are computed, and the critical speed of the entire catenary is determined as the intersecting points of the group and phase velocity curves. Therefore, the contributions of all catenary components and constraints on the critical speed are considered, offering a more thorough and accurate assessment than existing methods computing that of an isolate contact wire. The proposed method is extremely efficient, taking only 4.2% and 2.4% of computational time by COMSOL for modal and critical speed analyses, respectively. Finally, the influence of key parameters such as pretension, material and geometry parameters is systematically analysed.

Closed-form dynamic stiffness (DS) formulations are developed for exact dynamic analysis of complicated damped beam structures with internal hinges and generalized boundary conditions (BCs). The damping effects are considered by introducing both the strain rate dependent viscous and the velocity-dependent viscous damping factors. A range of theories are adopted for the beam element, that is, the classical and Rayleigh-Love theories for the axial vibration and the Euler–Bernoulli and Timoshenko theories for the bending vibration. Based on the governing differential equations considering the damping effect, the elemental DS matrices are derived by substituting the exact shape functions into the generalized displacement and force BCs. Then, a hinge joint for a beam element is simulated by implementing an end release. Next, assembling the elemental DS matrices and introducing generalized BCs leads to the final DS matrix used for dynamic analysis of beam structures. The complex transcendental frequency equation of the beam structures considered herein is solved by a global complex roots and poles finding algorithm with a self-adaptive mesh generator. The proposed method is verified against the finite element method. The results obtained by using different theories are also compared. This research provides an exact tool for dynamic analysis of complicated damped beam structures.