15 March 2025, Volume 44 Issue 5

  • Select all
    |
    VIBRATION AND MECHANICS SCIENCE
  • GAO Jianshe1, WU Jiabao1, WU Luji2, YANG Linjie2, RAO Xiaobo1
    Journal of Vibration and Shock. 2025, 44(5): 1-11.
    Abstract ( 116 ) Download PDF ( 100 )   Knowledge map   Save
    Here, in order to deeply investigate phenomenon of grazing collision in non-orthogonal face gear transmission system, a dynamic model of non-orthogonal face gear transmission system considering multiple strong nonlinear factors was constructed by dividing phase space into three parts with tooth face collision surface and tooth back collision surface as interfaces.Position Poincaré maps at two interfaces were established, and conditions for grazing bifurcation were calculated.Then, periodic variation diagram and phase diagram of meshing force were combined to find grazing bifurcation point, and solve Floquet multiplier near bifurcation point.By analyzing effects of grazing bifurcation on system dynamic behavior, grazing collision characteristics were revealed.Finally, the software Adams was used for simulation verification.The study showed that grazing collision can cause changes in impact state of gear system, but usually can’t affect motion state of gear system; when grazing bifurcation and period-doubling bifurcation occur simultaneously, gear system can have codimension-two grazing bifurcation; this moment, grazing collision can cause changes in both system’s motion state and impact state; the results can provide a theoretical basis for improving the reliability and service life of face gear transmission systems in engineering.
  • YUAN Pingping1, DING Yaxin2, YANG Chen1, REN Weixin3, FANG Wang2
    Journal of Vibration and Shock. 2025, 44(5): 12-17.
    Abstract ( 81 ) Download PDF ( 38 )   Knowledge map   Save
    To improve the accuracy of estimating instantaneous frequency with chirplet transform (CT), a new method for identifying structural instantaneous frequency was proposed based on CT combined with the pied kingfisher optimizer (PKO) and the radial basis function moving least squares (RBFMLS) algorithms.This method could use a positive definite tightly supported radial basis function as the weight function of moving least squares approximation to estimate the energy ridge of CT, and optimize RBFMLS node support radius and CT window function width using PKO.The effectiveness of the proposed method was verified through a case of numerical analytical signal and a time-varying cable test.The study results showed that the proposed method can effectively improve the energy aggregation of signal analysis and the identification accuracy of instantaneous frequency.
  • ZHANG Runqiang1, NAN Lingbo1, CHEN Diyi1, HAN Wei2, HUANG Weining3
    Journal of Vibration and Shock. 2025, 44(5): 18-25.
    Abstract ( 94 ) Download PDF ( 35 )   Knowledge map   Save
    Centrifugal pumps are important equipment in liquid transportation, they are widely used in fields of industry, agriculture and engineering.Understanding complex flow inside pumps is of great significance for pump design and optimization.Here, dynamic mode decompositionDMD)and reconstruction methods were introduced, and complex flow characteristics inside centrifugal pumps were studied.Firstly, the relation between modal energy and frequency was studied aiming at decomposition results of impeller.It was shown that low-frequency modes have higher energy, and modal energy decreases with increase in frequency.Subsequently, distribution of various modes in impeller flow channel was analyzed, and it was shown that main flow structures affecting complex flow inside impeller are distributed at inlet and outlet of impeller; increase in number of modes causes a more fragmented distribution of dynamic modes, so higher-order modes can capture more detailed flow structures.Temporal changes of various modes were studied, it was shown that effects of higher-order modes on flow field gradually decrease with increase in order number.Modal reconstruction methods were studied, it was shown that modal reconstruction can better construct the original flow field.Comparisons of reduced order reconstruction and full order reconstruction with the original flow field showed that the reduced order reconstruction method can also effectively reflect flow’s main characteristics, but some flow information may be lost, this part of information mainly is the amplitude distribution in flow field, it is stored in higher-order modes and related to complex flow at impeller outlet; DMD method can effectively extract main flow characteristics inside centrifugal pump, and provide a study basis for design optimization of vibration and noise reduction in centrifugal pumps.
  • REN Zhaopeng1, HUANG Zhe1, WEI Yi1, WANG Shenghai1, 2, SUN Yuqing1, 2, CHEN Haiquan1, 2
    Journal of Vibration and Shock. 2025, 44(5): 26-37.
    Abstract ( 118 ) Download PDF ( 22 )   Knowledge map   Save
    Considering the problem that the offshore-cranes cannot finish the hoisting assignment under the complex marine environment. Firstly, the dynamic model of the offshore-cranes under the stochastic ocean waves excitation was established and analyzed. Then the parallel cable-driven payload limiting-position mechanism (PCD-PLM) was designed and the adaptive anti-swing tension control method (PCD-AAT) was proposed to realize the stable control of payload. The payload remains stable state, even if the offshore-crane is affected by external excitation. The swing problem of payload is solved when the offshore-crane working under the heavy ocean conditions. Ensure that the hoisting assignment can be finished safely and efficiently. Finally, the offshore-crane anti-swing experiment platform was built to verify the effectiveness and practicability of the PCD-PLM and PCD-AAT. The simulation and experimental results show that the swing-angle of the payload reduced 99% and 80%, the stability control effects exceed 99% and 80%. The PCD-AAT show excellent dynamic tracking performance and robustness.
  • XIN Bingchen1, JIANG Bin1, CHENG Yuanqing2, BI Gang2, FAN Lili1, ZHAO Peiyi1
    Journal of Vibration and Shock. 2025, 44(5): 38-49.
    Abstract ( 55 ) Download PDF ( 20 )   Knowledge map   Save
    The instantaneous cutting energy consumption of the cutter teeth is an important indicator for evaluating the dynamic cutting performance of efficient milling cutters. The existed milling cutter cutting energy consumption solution model lacks of the distribution information of the instantaneous cutting energy consumption on the cutting edge of the cutter teeth, and cannot accurately identify the complexity and variability of the instantaneous cutting state at different positions of the cutting edge. Aiming at the problem that the distribution area of the instantaneous cutting energy consumption of the cutting edge cannot be accurately identified due to the uncertainty of the cutting boundary, a solution method for the instantaneous boundary of the main and secondary cutting edges is proposed to obtain the equations of the main and secondary cutting edges under vibration and wear. A solution model for the instantaneous cutting energy consumption distribution and its fractal dimension is established to obtain the distribution of the instantaneous cutting energy consumption on the cutter teeth. A method for identifying the dynamic distribution characteristics of cutting energy consumption is proposed. The extreme value, coefficient of variation and first-order difference of the fractal dimension are used to reveal the dynamic characteristics of the instantaneous cutting energy consumption distribution of different cutter teeth, and the experimental results of milling surface morphology and tooth flank wear depth are used to verify the proposed method.
  • LI Xianghe, ZHANG Zhimeng, JI Chunning, ZHAO Hanghao, ZHANG Yan
    Journal of Vibration and Shock. 2025, 44(5): 50-59.
    Abstract ( 90 ) Download PDF ( 39 )   Knowledge map   Save
    Two-dimensional numerical simulations of transverse flow-induced vibrations of cylinder arrays in uniform flow field are carried out using the iterative immersed boundary method. The cylinders are arranged orthogonally with the spacing ratios of L/D = 3.2, 8, the reduced velocity Ur = 2-11, Re = 100, and the mass ratio of the cylinders m* = 2. It is found that the cylinders in the front row and the rear row of the arrays exhibit the response characteristics of vortex-induced vibration and wake-induced galloping, respectively, for different spacing ratios. The amplitudes of the rear row of cylinders are generally larger than that of the front row of cylinders, and the maximum amplitudes of the cylinder arrays with small spacing ratios are higher than that of the case with large spacing ratios. For different spacing ratios, the vortex-shedding patterns of the neighboring cylinders in the front row (in the cross-flow direction) are stable in-phase-synchronized pattern and anti-phase-synchronized pattern, or the switching of the two patterns. In the in-line direction, the front-row neighboring cylinders at small spacing ratios show extended-body and reattachment modes, respectively, whereas they were all in co-shedding mode at large spacing ratios. There is a strong connection between the dominant frequency f*, amplitude and lift-displacement phase difference of each cylinder. With increasing Ur, the cylinder lift-displacement phase difference jumps from 0° to 180°, which corresponds to the excitation and suppression effects of the fluid on the cylinder vibration, respectively, where the rear row of cylinders need to transition to a larger reduced velocity interval to complete the phase transition. There is a small vibration frequency difference between neighboring side-by-side cylinders, and the phase difference between their displacements is unstable. With the increase of Ur, the range of the stabilized phase difference in the group of cylinders with large spacing ratio generally narrows upstream, while it expands downstream with small spacing ratio.
  • LUO Weirong, HAN Guangdong, ZHAO Shilong, HUO Ziteng, LI Jian, WANG Shenghai
    Journal of Vibration and Shock. 2025, 44(5): 60-66.
    Abstract ( 66 ) Download PDF ( 26 )   Knowledge map   Save
    Here, aiming at suspension cable being simplified to a straight line in current study on large-span flexible cable-driven parallel mechanism to cause its workspace being incorrect, considering flexible cable suspension effect, the workspace of the mechanism was studied.Firstly, a catenary suspension model and a parabolic suspension model were established.Comparison between test and simulation showed that under the same external conditions, the suspension length obtained with the parabolic suspension model is closer to test result than that obtained with the catenary suspension model.Based on the parabolic suspension model, static equations of a facade flexible cable-driven parallel mechanism were established, workspace and total power at each point in space considering suspension sag were analyzed.The simulation results showed that compared to the facade flexible cable-driven parallel mechanism under the linear suspension model, the workspace range with tension conditions based on the parabolic suspension model decreases by 12.6%, and the low-energy consumption area during lateral motion of mobile platform is significantly smaller than that during vertical motion.
  • DING Wangcai1, 2, MA Shuai1, JIN Zhongyuan3, LI Deyang1, WU Shaopei1, LI Guofang1
    Journal of Vibration and Shock. 2025, 44(5): 67-79.
    Abstract ( 77 ) Download PDF ( 40 )   Knowledge map   Save
    Taking the typical rail corrugation disease as the research object, a three-dimensional wheel-rail transient rolling contact finite element model considering the under-rail structure and wheel-rail elastic-plastic contact was constructed. The three-dimensional wear distribution calculation method based on the finite element model was used to study the wheel-rail contact characteristics and rail wear behavior under uniform speed and traction. The difference, as well as the influence of rail corrugation with different disease parameters on wheel-rail contact characteristics and wear under traction. The results show that the longitudinal force of wheel-rail under traction state is more significantly affected by the inherent stiffness irregularity. The influence of corrugation on relative slip is related to the longitudinal force, and it needs to be superimposed on the influence of inherent stiffness irregularity. The superposition effect is related to the wavelength and valley depth of corrugation. In the valley depth limit of corrugation, when the corrugation wavelength is short or the valley depth is large, the longitudinal force will have a vibration mode similar to the corrugation geometry, resulting in the relative slip distribution related to the corrugation geometry. The change of contact pressure further affects the wear depth, so that the wear near the crest increases, the wear near the trough decreases, and the corrugation further develops. Due to the small phase difference between the wear increment phase and the geometric phase of the corrugation, the rail corrugation tends to slow down, and the corrugation slows down faster under larger wavelength and valley depth. The research conclusions are expected to provide a theoretical reference for the rail maintenance strategy in high-speed wheel-rail corrugation section.
  • HU Bo1, 2, TAN Qingjie1, DONG Jianxiong1.2, XIAO Zeliang1, 2, YIN Lairong1, 2
    Journal of Vibration and Shock. 2025, 44(5): 80-87.
    Abstract ( 93 ) Download PDF ( 17 )   Knowledge map   Save
    Time-varying mesh stiffness is a crucial parameter in gear dynamics design and analysis. The lubricant filling the meshing tooth surface will change the contact characteristics of the gear pair, thereby affecting the time-varying mesh stiffness. The center distance of the small module gear pair is very small, and the angle of engagement, contact ratio and mesh stiffness are very sensitive to the center distance error. In light of the influence of the grease lubrication and center distance error on the contact characteristics of small module gears, the formula for calculating the oil film bearing stiffness was deduced in this study based on the elastohydrodynamic lubrication theory. Combining the potential energy method, a time-varying mesh stiffness calculation model for grease-lubricated gear pair was constructed to investigate the effects of the speed, viscosity and center distance error on the time-varying mesh stiffness. The research findings reveal that the secondary pressure peak and film thickness contraction phenomenon are evident in grease-lubricated small module gears, with the position of the peak shifting towards the exit with increasing load. The lubricating oil film filling the meshing tooth surfaces increases the contact area between meshing tooth surfaces, thereby enhancing the time-varying mesh stiffness of small module gear pair. On the contrary, the increase of center distance error will weaken the time-varying meshing stiffness of the gear pair. Oil viscosity enhances the strengthening effect of the oil film on the mesh stiffness of gear pairs with grease lubrication, but the increase of the gear speed weakens the mesh stiffness.
  • ZHAO Guoshou1, CAO Linlin2, WU Dazhuan2
    Journal of Vibration and Shock. 2025, 44(5): 88-96.
    Abstract ( 54 ) Download PDF ( 19 )   Knowledge map   Save
    Cavitation-induced loading oscillations are mainly responsible for the pump vibration, the excitation source identification of which is critical for the control of pump vibration. The water-tunnel experiments and numerical simulations are conducted to study the pump cavitation under the wake flow perturbation of inlet guide vanes. Meanwhile, an excitation source identification method based on second-order statistics is proposed, which connects the blade loading distribution and the overall thrust considering physical causality. Combined with the mode decomposition, the loading distribution of pump cavitation with a dominant energy contribution can be effectively revealed. The results showed that the temporal-spatial characteristics of inflow incidence angle largely regulate the development and distribution of blade cavitation. The correlation modes with amplitude can represent the energy contribution. The normalized correlation modes can represent the similarity between variables. For the cavitation condition, it is found that the whole part from the blade middle to the hub contributes much to the thrust oscillations through excitation source identification based on the second-order statistics, which is dominated by low frequency. However, only the cavitation-induced loading oscillations can be extracted by first-order statistics.
  • TANG Dalun1, LIU Yangzhao1, DAI Kaoshan1, SHI Yuanfeng1, PU Qiong2, Del ChungMing Yang2, David Yáez3
    Journal of Vibration and Shock. 2025, 44(5): 97-105.
    Abstract ( 107 ) Download PDF ( 39 )   Knowledge map   Save
    In order to explore the wind vibration characteristics of bladeless wind turbines, a bladeless wind turbine with a height of 3 m and a rated power of 100 W was taken as the research object. Firstly, the finite element model of the wind turbine is established, the natural mode information is calculated, and the wind-induced vibration is predicted. Subsequently, a full-scale model of the turbine was tested in a wind tunnel. Finally, the experimental results are compared with the theoretical prediction results, and the wind-induced vibration law is comprehensively analyzed. The results show that the cross-wind displacement at the top of the fan increases first and then decreases with the increase of wind speed, showing a typical vortex phenomenon, so there is no need to consider safety measures such as braking system in high wind speed environment. At the same time, when the wind speed is stable, the displacement time history of the top of the fan is close to the standard sinusoidal curve, showing stable wind-induced vibration. In addition, the aerodynamic characteristics of the wind turbine in all directions are also about the same, with good adaptability to all wind directions, and there is no need to rely on an additional yaw system to cope with wind direction changes.
  • PENG Ziyan, ZHU Hongtao, YANG Fufeng
    Journal of Vibration and Shock. 2025, 44(5): 106-115.
    Abstract ( 99 ) Download PDF ( 28 )   Knowledge map   Save
    An internal suspension sensor was designed for the rotary magnetorheological damper in armored vehicle suspension systems. This sensor integrates the measurement of the damper arm rotation angle, single-axis acceleration, and temperature into a single-board system. The three sensor signals are digitized through AD sampling, modulation, transmission, and demodulation, enabling the digital transmission of multiple sensor signals. Through thermal and magnetic field simulations of the installation structure, the insulation, heat dissipation, and magnetic circuit structures of the sensor were determined. Simulation results indicate that the sensor can still operate stably within the working range of the magnetorheological damper and under external magnetic field interference, demonstrating good environmental adaptability, as well as the advantages of being compact and highly resistant to interference.
  • QIU Xintao, LIU Ning, ZHANG Yabin
    Journal of Vibration and Shock. 2025, 44(5): 116-121.
    Abstract ( 60 ) Download PDF ( 17 )   Knowledge map   Save
    This paper presents the design of an automatic loading system for mortars equipped with dual coordinating arms and establishes a dynamic model of the coordinating arm automatic loading system. An experimental scheme for a single coordinating arm automatic loading system is designed to verify the accuracy of the simulation results. A comparative study of the dynamic simulations between the dual coordinating arm and single coordinating arm automatic loading systems is conducted, exploring the dynamic characteristics of the dual coordinating arm automatic loading system over different time intervals. The results indicate that the dual coordinating arm significantly mitigates vibrations compared to the single coordinating arm; as the system stiffness increases, the ammunition supply process becomes smoother. The acceleration of the dual coordinating arm exhibits convergence, with longer system lengths resulting in extended loading times. The optimal structural dimensions for the dual coordinating arm are identified, leading to a 34% increase in loading speed compared to the single coordinating arm. This study investigates the dynamic characteristics of the coordinating arm automatic loading system and the impact of structural variations on loading speed, laying the foundation for further structural optimization and precision control.
  • SHI Shengbo, YUAN Xinyi, SU Yan
    Journal of Vibration and Shock. 2025, 44(5): 122-130.
    Abstract ( 83 ) Download PDF ( 23 )   Knowledge map   Save
    Under the action of horizontal harmonic excitation, the amplitude response curve of the internal fluid free surface motion of the shallow-water tuned liquid damper (TLD) will show a sudden drop at the jump frequency point. In this paper, based on the third-order nonlinear amplitude time domain equation, the approximate calculation formula of jump frequency considering the influence of viscous damping is obtained by ignoring the high-order term through dimensional analysis. In order to verify the validity of the formula, a series of experiments were carried out through a six-degree-of-freedom motion simulator and a shallow-water TLD tank. The calculation results of the proposed formula were compared with the calculation results of the existing formula and the experimental measurement results. It is found that the new formula is more consistent with the experimental measurement results as a whole, which is better than the existing calculation formula. At the same time, the new formula comprehensively considers the influence of water depth, characteristic length, excitation amplitude and viscosity, and has better applicability in the field with greater damping influence.
  • CIVIL ENGINEERING
  • YAO Xiaojun1, SUN Shoupeng1, WANG Qiang1, YANG Xiaomei2
    Journal of Vibration and Shock. 2025, 44(5): 131-139.
    Abstract ( 104 ) Download PDF ( 44 )   Knowledge map   Save
    To locate the instant and location of abrupt structural damage accurately, this paper proposes a damage identification method based on variational mode decomposition (VMD) and the autoregressive integrated moving average model (ARIMA). Initially, the method accurately selects the initial central frequencies of high-energy modes and the number of modes to be decomposed using the autoregressive model power spectrum. Then, the variational mode decomposition method is applied to decompose the vibration non-stationary signal into several stationary signal components, followed by fitting each order signal component using the ARIMA model to obtain model residuals. The specific time of damage is determined by analyzing the residuals obtained from fitting the signal components with ARIMA. Finally, the structural modal shapes are obtained using principal component analysis (PCA), and a damage index that integrates frequency and modal shape is constructed. Combined with a damage threshold, the location of the damage is pinpointed. Furthermore, the method is analyzed using numerical examples of a ten-degree-of-freedom structure under seismic excitation and the monitoring acceleration data of an actual steel bridge. The results demonstrate that the proposed method can accurately locate the occurrence time and location of structural damage under both stationary and non-stationary excitations.
  • ZHU Jin1, CHEN Yuanwen1, WU Mengxue2, ZHANG Yu3, HENG Junlin4, LI Yongle1
    Journal of Vibration and Shock. 2025, 44(5): 140-148.
    Abstract ( 77 ) Download PDF ( 27 )   Knowledge map   Save
    To investigate the fatigue issues of welds in steel bridges under severe corrosion environments and heavy loading, a three-stage probabilistic corrosion fatigue damage evolution model was developed based on previous experimental results. The model is divided into three key stages: corrosion-driven growth, competition between corrosion and fatigue, and fatigue-driven growth. The model also incorporates the effects of multiple failure paths, considering the initiation, coalescence, and propagation of multi-source fatigue cracks in the welds. Using a coastal bridge as the engineering background, a reliability analysis was conducted through multi-scale finite element simulations combined with the probabilistic corrosion fatigue damage evolution model, integrating the measured vehicle load spectrum and environmental characteristics of the bridge. The results indicate that, compared to the weld toe failure mode, the corrosion fatigue life of the weld root failure mode is reduced by 11.3%, and the dual failure mode life is reduced by 16.3%, with the weld root failure mode being dominant. Under a C5 environment, the corrosion fatigue life is reduced by 26.9% compared to the C4 environment, indicating that the higher the environmental corrosion grade, the shorter the corrosion fatigue life. When the traffic volume coefficients are 1.2% and 2.4%, the corrosion fatigue life decreases by 29.9% and 38.8%, respectively, compared to constant traffic conditions, indicating that the higher the traffic volume, the shorter the corrosion fatigue life. Under the combined effects of severe environmental corrosion and heavy traffic loads, the corrosion fatigue life is reduced by 47.8%, demonstrating significant degradation of the fatigue resistance of welds in steel bridges under the coupled effects of severe corrosion and heavy loads.
  • ZHOU Wenjun1, ZHENG Shixiong2, ZHU Donghong1, SHEN Guohui2
    Journal of Vibration and Shock. 2025, 44(5): 149-154.
    Abstract ( 51 ) Download PDF ( 23 )   Knowledge map   Save
    Wind tunnel tests were conducted on a complex urban landform to obtain the mean wind velocity profile of 24 wind directions and fit the roughness power index to determine the ground roughness categories for each wind direction angle. Meanwhile, the influence of the upwind sector angle on the test results was investigated. The results indicate that upwind areas with mixed urban and mountainous landforms are more likely to develop into Category C landforms. The goodness of fit of the roughness power index values for 24 wind directions is greater than 0.95, indicating that the exponential law can effectively express the mean wind velocity profile of complex urban landforms. The upwind sector angle has a significant impact on the classification of urban ground roughness. The average wind speed profile and fitted power index are very close when the sector angle is between 70°and 180°, while the index decreases when the angle is between 30°and 60°. It can be determined that the minimum angle of the upstream sector is 70°, which is close to the 90° specified in ASCE and EN specifications. The upstream sector value of 180°, which is equal to a semicircle, regulated in the Load Code for the Design of Building Structures, seems to be too large.
  • WANG Li1, PAN Qiren1, GU Haowei1, HU Qi1, 2, YU Lusong1, LI Ziqi1, 2
    Journal of Vibration and Shock. 2025, 44(5): 155-168.
    Abstract ( 120 ) Download PDF ( 35 )   Knowledge map   Save
    To study the seismic performance of steel tube confined concrete columns under different ambient temperatures, nine quasi-static tests of steel tube confined concrete columns under different temperature conditions were conducted. The test phenomena and failure modes of each specimen under different temperatures were analyzed, revealing the mechanism of the influence of ambient temperature on the seismic performance of steel tube confined concrete columns. In view of the lack of consideration of environmental effects in the current calculation method for the theoretical skeleton curve model of steel tube confined concrete, based on experimental and finite element calculation results, a temperature correction formula for the characteristic points of the theoretical skeleton curve of steel tube confined concrete considering the influence of environmental temperature is proposed. The results show that the material variation of steel tube confined concrete caused by the change of ambient temperature and the Temperature-induced additional constraint effect have a significant effect on the horizontal bearing capacity and ductility of the specimens. Under the condition of high temperature (60 °C), the bearing capacity and ductility of the specimens decrease by 9.4 % and 12.2 % respectively due to the decrease of the strength of the core concrete and the debonding of the steel-concrete interface. Under low-temperature (-40℃) condition, due to the increase in the strength of the core concrete, and at the same time, the outer wall steel pipe strengthens the restraining effect on the core concrete, so that the specimen bearing capacity increases by 18.7%, but its ductility decreases by a maximum of 47.2%, which has a negative impact on seismic resistance. The theoretical skeleton curve of steel tube confined concrete considering the influence of ambient temperature proposed in this paper has good reliability and accuracy, which can provide the necessary basis for the design of steel tube confined concrete in high cold and large temperature difference areas.
  • TAN Chao1, QUAN Yong1, ZHANG Zhengwei2, YAO Bo3, FU Guoqiang1
    Journal of Vibration and Shock. 2025, 44(5): 169-175.
    Abstract ( 123 ) Download PDF ( 36 )   Knowledge map   Save
    Compared with aeroelastic model wind tunnel test, forced-vibration model test offers a high degree of controllability in terms of parameters and is a significant means to study the aeroelastic effects of super high-rise buildings. However, the manufacturing difficulties of forced-vibration devices limit the widespread application of this method. In this study, Large Eddy Simulation (LES) is used to simulate the forced-vibration model wind tunnel test of super high-rise buildings, in order to investigate their aeroelastic effects. The results show that the proposed method in this study can effectively simulate the aerodynamic forces and flow field of forced-vibration model for super high-rise buildings. When the incoming wind speed is within the range of vortex-induced resonance (VIR) speeds, the increase in cross-wind amplitude of super high-rise buildings will significantly aggravate the vortex shedding on the side of the structure, resulting in a obvious increase in the side wind pressure and the root mean square values of the lift coefficient. However, when the reduced wind speed is far above the VIR speed, the cross-wind vibration of the structure has a little impact on its aerodynamic forces. The changes of amplitude and reduced wind speed have significant effects on the cross-wind aerodynamic stiffness and aerodynamic damping of buildings.
  • YANG Han1, 2, ZHENG Shixiong1, 2, MA Cunming1, 2
    Journal of Vibration and Shock. 2025, 44(5): 176-183.
    Abstract ( 134 ) Download PDF ( 34 )   Knowledge map   Save
    Here, to study effects of attachments on vortex induced vibration (VIV) of a triple-box girder bridge, taking a certain highway and rail way dual-use separated triple-box girder bridge as background, wind tunnel tests were conducted in XNJD-1 reflux series wind tunnel.Multiple sets of operating conditions were set up in tests to study effects of rail maintenance, single beam attachments and combinations of attachments on VIVs of triple-box girder.After analyzing the test results, it was shown that effects of rail maintenance on VIV characteristics of triple-box girder are limited; highway wind barriers can deteriorate both vertical and torsional VIVs, while railway wind barriers can suppress VIVs; crash barriers can suppress vertical VIVs but can also induce torsional VIVs; compared with crash barriers, railway wind barriers perform better in suppressing VIVs, while highway wind barriers have a larger deteriorating effect on vertical VIVs compared to suppressing VIV effects of railway wind barriers and crash barriers; highway wind barriers are main excitation source of torsional VIV, and working conditions with highway wind barriers produce severe torsional VIVs.Finally, by using computational fluid dynamics  numerical simulation, flow field characteristics of triple-box girder were analyzed, the results showed that attachments can promote generation and development of vortices in main beam gap, and intensify VIVs of triple-box.
  • EARTHQUAKE SCIENCE AND STRUCTURE SEISMIC RESILIENCE
  • WANG Changsheng1, LIN Jianhao1, YANG Yan1, XU Jiayun2
    Journal of Vibration and Shock. 2025, 44(5): 184-190.
    Abstract ( 71 ) Download PDF ( 28 )   Knowledge map   Save
    Aiming at the problem that the coupled modal damping matrix cannot be accurately decoupled in the real number domain in the motion equation of soil-structure interaction system. Firstly, the forced decoupling method is used to decouple the non-classical damping matrix, and then the error caused by the forced decoupling method is analyzed. The real mode approximate decoupling method is proposed to solve the structural system. At the same time, the system response is expressed by the linear combination of displacement and velocity of a series of standard oscillators combined with Laplace transform. Through the example analysis, it can be seen that the seismic response of the structure obtained by the real mode approximate decoupling method is in good agreement with the results obtained by the accurate complex mode method, and its accuracy is higher than that of the forced decoupling method. Especially in the analysis of the dynamic response of the superstructure of the soil-structure interaction system, its advantages are more prominent. The proposed real mode approximate decoupling method has high accuracy, avoids the operation in the complex domain, and is easy to understand in engineering significance. It can be applied to other structural systems with non-classical damping characteristics.
  • ZHENG Xingqun, ZHAO Boming
    Journal of Vibration and Shock. 2025, 44(5): 191-196.
    Abstract ( 74 ) Download PDF ( 21 )   Knowledge map   Save
    All strong ground motion records in southwest China from 2009 to 2022 were collected, and small ground motion records of strong earthquake data below magnitude 5 were selected to form a data set. The geometric spreading term and the anelastic attenuation term were determined separately. Firstly, the least square method was used to piecewise fit the geometric segment points value, and the calculated results are 74 km for R1 and 125 km for R2. Then the four regional crustal medium parameters were obtained by genetic algorithm inversion: the stress drop is 77.06 MPa, the quality factor Q0 and η are 223.13 and 0.5697 respectively, and the near-surface high frequency attenuation parameter κ0 is 0.0686 s. Based on the inversion of crustal parameters, the acceleration Fourier spectrum was calculated, and the time history of ground motion under a certain magnitude and epicenter was obtained with random phase spectrum. The attenuation relationship of peak acceleration of ground motion in southwest China was established by extracting PGA. The results show that the attenuation relationship established is in good agreement with the actual strong earthquake records, especially in the near field area, which will provide important reference for the engineering construction in southwest China.
  • LIU Xiaojuan1, 2, HUANG Zhexuan2, LIU Yang1, 2, GUO Zixiong1, 2, JIANG Huanjun3, CHEN Yujie2
    Journal of Vibration and Shock. 2025, 44(5): 197-206.
    Abstract ( 66 ) Download PDF ( 19 )   Knowledge map   Save
    Due to the influence of external environment, the steel corrosion ratio in the reinforced concrete (RC) column is often inconsistent, and the influence mechanism of tensile and compressive longitudinal reinforcement corrosion on the compressive-flexural performance of RC columns is also different. Quasi-static loading tests were carried out on 1 uncorroded RC column and 5 corroded RC columns to investigate the effect of longitudinal reinforcement corrosion on the seismic performance of RC columns. The effect of corrosion ratio of steel reinforcement, corrosion location and axial compression ratio on the seismic behavior of RC columns were analyzed. The results show that the damage distribution and mechanical properties of RC columns with corroded longitudinal reinforcements in single side are asymmetrical. For the RC column with axial compressive ratio of 0.1, the fracture of corroded tensile longitudinal reinforcement is the main reason for the mechanical degradation of the column, while for RC columns with large axial compressive ratio, the compressive damage of corrosive concrete and the buckling of corroded compressive longitudinal reinforcement are the main causes of mechanical performance deterioration of the column. Compared with the uncorroded RC column, the positive loading capacity of the columns corroded with corrosion ratio of longitudinal reinforcement in single side not exceeding 15% decreased by less than 10%, the reverse loading capacity decreased by less than 5%, and the cumulative energy consumption at ultimate limit state decreased by about 6.48% to 15.21%. In comparing with the load-carrying capacity and energy dissipation capacity, the deformation capacity of RC columns is effected by the reinforcement corrosion the most significantly. With the increase of steel corrosion ratio, the proportion of bending deformation of RC columns decreases, and the proportion of shear deformation and slip deformation increases. For the column with tensile longitudinal reinforcement premature fracture due to corrosion, the positive limit deformation decreased by about 30% to 37%, while the reverse limit deformation is reduced by about 10%~17%. The seismic evaluation of corroded RC columns should properly consider the non-uniformity corrosion of tensile and compressive longitudinal reinforcement in the column. 
  • QI Liangjie1, 2, CAO Yi1, YUAN Zhen1, XUE Jianyang1, 2, LIU Mengda1
    Journal of Vibration and Shock. 2025, 44(5): 207-217.
    Abstract ( 76 ) Download PDF ( 17 )   Knowledge map   Save
    To improve the rapid post-earthquake reparability of steel structure antique buildings, this paper proposes a shape memory alloy (SMA) bar-friction damper series device suitable for the characteristics of antique buildings. The cyclic loading tests were conducted on beam-column joints in steel structure antique buildings with and without the proposed device. Properties such as damage characteristics, load-displacement curves, stiffness, hysteretic energy dissipation and recovering capacity of these two types of joints were analyzed. The results showed that only the SMA bars are stretched when the loading displacement angle does not exceed 1.8%, the friction damper and the SMA bars work simultaneously with the increasing of the loading displacement. Installing the SMA bar-friction damper series device around beam-column joints can effectively improve the seismic performance of the steel antique building, the ultimate strength and the cumulative energy dissipation is increased by 151% and 86%, the residual displacement is decreased by 27%, the damage and deformation of the specimen mainly occurs at the steel angles. The finite element models of beam-column joints with SMA bar-friction damper series device in steel antique buildings was established and its parameters are analyzed. The results showed that the bearing capacity of the proposed joint is increasing with the increase of sliding force in friction dampers, and it decreases with the increase of column axial compression ratio, while the effect of the SMA bar preload is negligible.
  • CHEN Yuliang1, 2, 3, ZHANG Shaosong1, WANG Shuangyi1, YE Peihuan1, 2, JI Yunpeng1
    Journal of Vibration and Shock. 2025, 44(5): 218-229.
    Abstract ( 51 ) Download PDF ( 20 )   Knowledge map   Save
    In order to explore the calculation method of torsional bearing capacity and hysteretic behavior of steel reinforced concrete columns with welded studs (SSRC) columns under the combined action of compression, bending, shear and torsion, the quasi-static loading tests of 12 SSRC columns and 1 control column were designed and completed with the parameters of steel ratio, longitudinal reinforcement ratio, volume stirrup ratio, stud spacing and stud position. The failure modes of the specimens were observed, and the seismic performance indexes such as hysteresis curve, skeleton curve, bearing capacity, energy dissipation and deformation performance of each specimen were obtained. The influence of different parameters on the above indexes was analyzed, and the calculation formula of torsional bearing capacity of this kind of members was proposed. The results show that the variation parameters have little effect on the failure mode of the specimens, and all the specimens show obvious bending and torsion failure. The torque-torsion angle hysteresis curve of the specimen has obvious pinching phenomenon, which is shown as an inverse 'S' shape, but the bending moment-displacement hysteresis curve is relatively full, which is shown as a spindle shape. The steel ratio and stirrup ratio have little effect on the torsional and flexural capacity of SSRC columns. The longitudinal reinforcement ratio and stud spacing have a great influence on the torsional bearing capacity of SSRC columns, and the maximum increase is 12.29% and 18.92% respectively, but the influence on the flexural bearing capacity of SSRC columns is not significant. Welded studs can significantly enhance the ductility and energy dissipation capacity of steel reinforced concrete columns. Considering the seismic indexes such as bearing capacity, energy dissipation and ductility, it is recommended to use X-type welding for studs. Finally, a practical formula for calculating the torsional bearing capacity of SSRC columns is proposed. The calculated results are in good agreement with the test results.
  • CHEN Shitong1, 2, 3, FAN Xin1, 2, GAN Pingping4, GUO Chaoqun2, XIANG Min4
    Journal of Vibration and Shock. 2025, 44(5): 230-242.
    Abstract ( 58 ) Download PDF ( 33 )   Knowledge map   Save
    Compared with conventional bridges, H-type double-column super-high pier continuous rigid frame bridge has more complicated stress, significant influence of high-order vibration modes and great difference in longitudinal and transverse stiffness. The incident angle of ground motion is one of the key factors affecting its seismic response. In order to explore the influence of incident angle on its fragility, a simulation model is established to analyze the seismic response of a three-span super-high pier continuous rigid frame bridge of a high-speed railway in western mountainous areas. Based on the longitudinal and transverse bridge directions, the bi-directional Fragility analysis of components is carried out, and the contribution rate of super-high piers and bearings to the bridge system is considered. Based on the first-order reliability theory, the Fragility model of bridge composite system is constructed, and the fragility surface of bridge system is established, and the influence law of earthquake incidence angle on the seismic fragility of bridge system is discussed. The results show that the damage probability of the middle part of the ultra-high pier under longitudinal earthquake excitation is obviously higher than that of the middle and low piers, and the peak curvature of the transverse bridge is out of sync with the capacity demand ratio, so it is not reasonable to evaluate the seismic performance of the ultra-high pier only by judging the weak part of the peak curvature. The most unfavorable incident angles at the bottom of pier are 0~15 and 165~180, and the most unfavorable ones at the cross beam of pier body are 60~120, and the most 
    unfavorable ones at the top of pier vary under different earthquake intensities. The seismic fragility of bearing is sensitive to the change of earthquake intensity and incident angle, and it presents different damage trends under different working conditions, and it is more prone to damage than pier, so the mechanical properties of transverse stop should be paid attention to in design. The seismic fragility of bridge system shows a strong correlation with the incidence angle of ground motion, and the most unfavorable incidence angles are 75~105 and 255~285.
  • SHOCK AND EXPLOSION
  • PENG Birong1, XUE Xiaochun1, CAO Yongjie2, HUANG Lei1, YU Yonggang1, WANG Hongjin2
    Journal of Vibration and Shock. 2025, 44(5): 243-252.
    Abstract ( 97 ) Download PDF ( 39 )   Knowledge map   Save
     In order to solve the problem of excessive system recoil caused by the large kinetic energy of the muzzle of a 40mm submerged bullet artillery gun, a composite recoil reduction technology that combines the muzzle brake and expanding wave technology is proposed. Based on this technology, the internal ballistic equation of the submerged bullet expanding wave artillery was established, and the change rule of the chamber pressure, velocity and body tube force with time during the internal ballistic period was analyzed. On this basis, based on the three-dimensional non-constant Navier-Stocks equation, combined with the dynamic mesh technology, the numerical simulation of the breech flow field during the after-effect period of the submerged projectile gun is carried out, which focuses on revealing the influence of the gunpowder and gas flow on the force characteristics of the retractor, and comparing and analyzing the recoil and recoil impulse change characteristics of the two submerged projectile guns containing the composite recoil reduction technology and the two types of submerged projectile guns without this technology in the entire internal ballistic path and the after-effect period. The results of the study show that the recoil and recoil impulse of the two types of submerged projectiles with the composite recoil reduction technology and without the technology are analyzed in comparison. The results of the study show that the composite recoil reduction technology can effectively reduce the recoil force of the submerged projectile gun and reduce the recoil impulse from 1720.36 N-s to 306.04 N-s, with a recoil reduction efficiency of 82.2%, while keeping the muzzle velocity of the projectile unchanged.
  • KONG Linghao1, 2, YAN Peng1, 2, LIU Xiao1, 2, ZHOU Chao1, 2, ZHANG Xiangyu1, 2, LU Wenbo1, 2
    Journal of Vibration and Shock. 2025, 44(5): 253-262.
    Abstract ( 118 ) Download PDF ( 36 )   Knowledge map   Save
    Identification of the first arrivals of different waveform components in blasting seismic waves rapidly and correctly is crucial for the analysis of energy and frequency characteristics, wave velocity calculation, precise and control of blasting vibration effects and inversion of rock dynamics parameters. In this paper, by analysing the propagation law and energy and frequency characteristics of blasting seismic waves, improving the existing P-wave first arrivals identification method - STA/LTA algorithm and S-wave first arrivals identification method - eigenfunction method, a new method of automatic identification of first arrivals moments of blasting seismic waves is proposed and verified by carrying out blasting tests. The results show that: compared with the manual identification of P-wave arrival results, the maximum error of P-wave first-to-moment identification in this paper is 0.3ms, and the average error is 0.1ms; the peak value of S-wave first-to-moment identification eigenfunction is more obvious, which can improve the automatic identification efficiency. The recognition method has lower error, can effectively identify the blasting seismic wave P-wave and S-wave first-arrive-moment.
  • ZHANG Xin1, 2, LIU Zegong1, 2, CHANG Shuai2, CHEN Xiangsheng2, XUE Yonglin2 ZHU Jialiang2, SONG Xin2
    Journal of Vibration and Shock. 2025, 44(5): 263-277.
    Abstract ( 103 ) Download PDF ( 30 )   Knowledge map   Save
    In order to investigate the performance of Riedel-Hiermaier-Thoma(RHT) and Holomquist-Johnson-Cook(HJC) constitutive models in the numerical simulation of blasting in deep coal and rock bodies, the parameters of RHT and HJC constitutive models were calibrated based on the deep coal in Huaibei mining area; then the validity of the calibrated parameters was verified through the indoor tests and numerical simulations of single-hole blasting, and the composite failure criterion of maximal tensile stress and maximal shear stress was added to improve the HJC model. Finally, the numerical simulation of single-hole blasting under different ground stresses is carried out based on the two constitutive models to explore the applicability of the two constitutive models in simulating the range of damage zones and rupture patterns. The numerical simulation results of single-hole blasting show that both the RHT model and the improved HJC model are able to describe the formation of the crushed zone and crack evolution under no ground stress loading condition; the numerical simulation results of the RHT and HJC models under different ground stress loading conditions show that the ground stress has a restraining effect on the crack expansion, and the quantitative indexes of the expansion range of the crack zone show that the RHT model is closer to the theoretical value than the HJC model. The RHT model is closer to the theoretical calculation value than the HJC model. According to the results of the study, the RHT model and the improved HJC constitutive model are both suitable for the numerical simulation of the blasting problem in the absence of ground stress, and the RHT constitutive model should be the first choice for the blasting in the high ground stress condition.
  • ZHU Feixiang, GAO Fei, LIU Chenkang, DENG Shuxin
    Journal of Vibration and Shock. 2025, 44(5): 278-288.
    Abstract ( 130 ) Download PDF ( 31 )   Knowledge map   Save
    In order to explore the law of explosive damage of shell charges in concrete, the drill and blast method is used to carry out explosion tests and numerical simulation of shallow and deep shell charges under different charge depths and shell thickens. The macro morphology changes of shallow explosive projectile crater and the plastic damage evolution law of deep buried closed explosive concrete are analyzed. The results show that under shallow-buried explosive condition with shell charge, the energy acting on concrete destruction increases with the increase of charge depth, and the macro size of crater increases. According to the macroscopic morphology of the crater, the crater forming area can be divided into funnel crater area, expanding hole area and warhead action area. The influence of shell thickness on concrete pitting size is basically the same under different charging depths. When the ratio of shell thickness to charge diameter δ/de is with in 0.075~0.125, the damage effect increases with the increase of shell thickness, and the damage effect is the best and tends to be stable at 0.125~0.15, and the damage effect decreases when the shell thickness continues to increase. On the one hand, the existence of shell can produce agglomeration effect of detonation products and enhance the damage effect, on the other hand, the thickness of shell will consume the detonation energy and weaken the damage effect, so there is a certain shell thickness to make the damage effect best.
  • CHEN Gongqing1, WU Hao1, OU Yuan2, WANG Fan2, L Jinxian1
    Journal of Vibration and Shock. 2025, 44(5): 289-301.
    Abstract ( 87 ) Download PDF ( 31 )   Knowledge map   Save
    Damage analysis of ground building under blast loading has significant practical guidance for developing combat strike strategies and engineering protection design. The LS-DYNA finite element analysis software is used to reproduce the existing near-field explosion test of the reinforced concrete (RC) masonry-infilled frame structure, which fully verifies the applicability of the adopted refined numerical simulation method. Combined with the hybrid element modeling method of building structures, a simulation analysis is carried out on the dynamic response of the three-story masonry-infilled RC frame under the explosion of typical warheads (100kg and 200kg TNT equivalent). The propagation of blast waves inside the structure and structural damage characteristics are investigated. Based on the equivalent single-degree-of-freedom (SDOF)method, the damage degrees of RC beams, columns, slabs and infill walls under blast loads are predicted, and a damage assessment method for building structures under internal explosions is established. Its applicability is verified by comparing with the refined numerical simulation results. It indicates that, under 100kg and 200kg TNT explosion conditions, the overall functionality and structural damage degrees of the building in the refined numerical simulation are both moderate and slight. The corresponding damage degrees obtained by the equivalent SDOF simplified assessment are consistent with the refined numerical simulation results. In addition, it can be seen from the damage degree of members that compared with load-bearing components, i.e. slabs, beams and columns, masonry infill walls are more prone to failure, resulting in a larger damage range of rooms in the horizontal direction within the floor.
  • FAULT DIAGNOSIS ANALYSIS
  • XU Huanwei, ZHAO Zewei, XIAO Xuyuan, WANG Zhonglai
    Journal of Vibration and Shock. 2025, 44(5): 302-313.
    Abstract ( 101 ) Download PDF ( 55 )   Knowledge map   Save
    In actual industrial production, different operating conditions lead to variations in data distribution, posing a challenge for bearing fault diagnosis under different working conditions. To address this issue, a fault diagnosis method based on multi-adversarial and balanced distribution adaptation was proposed. Firstly, an improved residual network was used to directly extract domain-invariant features from the original vibration signals, enhancing feature extraction efficiency while preserving rich fault feature information. Secondly, a domain adaptation method combining correlation alignment and multi-adversarial domain adaptation was proposed, which can simultaneously align marginal distribution and conditional distribution of source domain and target domain to minimize data distribution differences between domains.Thirdly, the balanced distribution adaptation method was improved with designing a balance factor to allocate weights to the marginal distribution and conditional distribution in the adaptation process, so as to enhance cross-domain fault diagnosis effect. Finally, the effectiveness of the proposed method was validated using publicly available bearing fault datasets. Experimental results show that compared to popular domain adaptation methods, the proposed method achieves higher fault diagnosis accuracy, showing practical application value in bearing fault diagnosis tasks under different working conditions.
  • YUE Xia, LI Zhibin, ZHANG Chunliang, WANG Yadong, WANG Yuhua, LONG Shangbin
    Journal of Vibration and Shock. 2025, 44(5): 314-322.
    Abstract ( 81 ) Download PDF ( 25 )   Knowledge map   Save
    The joint robot is used in all kinds of production links, real-time monitoring of the load is the premise to ensure the safe operation of the robot. However, in some special scenarios, the load cannot be measured directly, and the dynamic method is usually used to solve the problem indirectly. Due to its obvious nonlinear characteristics and the uncertainty of model parameters, the accuracy and efficiency of load identification are not high. Therefore, this paper proposes an improved model based on Fourier neural network to realize load identification, so as to improve the prediction accuracy and timeliness of system load parameters. The proposed method uses the convolution and frequency domain truncation mechanism of Fourier neural network to obtain the characteristic signal quickly, and fuses the output result of feedforward neural network to obtain the identification result. Compared with the dynamic model solving method, the proposed method has higher precision and faster calculation speed. It only needs to learn several interphase sample sets within the prediction range to identify any result within the prediction range, and has good generalization ability. At the same time, the network sensitive parameters are analyzed, and the performance is compared with the mature neural network algorithm. In this method, two kinds of neural network models cooperate with each other, which can effectively identify different feature sets in high-dimensional data, so as to realize parameter identification, and provide reference for parameter identification of complex nonlinear systems.
  • LIU Shiya1, LIANG Wensheng1, HE Jun1, CHEN Zhiwen2, DAI Lei3
    Journal of Vibration and Shock. 2025, 44(5): 323-330.
    Abstract ( 71 ) Download PDF ( 28 )   Knowledge map   Save
    Most of existing domain adaptation-based scenarios used to estimate inter-domain distribution discrepancy can ignore the impact aroused by conditional distribution discrepancy between classes in two domains on the measurement of inter-domain distribution discrepancy, which can make the accuracy of knowledge transfer drop dramatically. A fault diagnosis method based on one-dimensional convolutional auto-encoder and conditional distribution domain adaptation is proposed for rolling bearings across rotating machinery. An inter-domain dissimilarity criterion (T-LCORAL) based on threshold local subdomain correlation alignment is constructed, which is embedded into the output layer of the autoencoder to guide the model for feature extraction effectively. Then, pseudo labels are employed to divide the source and target domains into multiple subdomains for fine-grained alignment within the domain, Moreover, the pseudo labels are screened by reliability thresholds to improve the discriminability of the boundary distribution between classes in the domain. Finally, three bearing datasets were used for cross mechanical bearing fault transfer diagnosis experiments, and the results showed that the proposed method can achieved better performance.