• ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus
  • Indexed by Core Journals of China, Chinese S&T Journal Citation Reports
  • Chinese S&T Journal Citation Reports
  • Chinese Science Citation Database

Online First

Articles in press have been peer-reviewed and accepted, which are not yet assigned to volumes/issues, but are citable by Digital Object Identifier (DOI).
Display Method:
Review on Distributed Energy Harvesting Technology for Small Unmanned Carrier Vehicles
WANG Hao, ZHANG Zutao, CUI Shuyan, YANG Lin, ZOU Rui, KONG Lingji
, Available online  , doi: 10.3969/j.issn.0258-2724.20250475
Abstract:

To meet the requirements for long endurance and high autonomy of small unmanned carrier vehicles (unmanned aerial vehicles, unmanned ground vehicles, and unmanned surface vehicles), distributed energy harvesting technology has become a crucial approach to overcoming energy supply bottlenecks. First, the basic principles of distributed energy harvesting technologies were reviewed, mainly covering the conversion mechanisms and key material advances of ambient energy sources such as solar, wind, vibration, and wave energy. Then, the structural integration and application characteristics of photovoltaic, electromagnetic, piezoelectric, and triboelectric harvesting technologies in typical scenarios involving unmanned aerial vehicles, unmanned ground vehicles, and unmanned surface vehicles were deeply analyzed. The research status of multi-source hybrid harvesting and self-powered systems was summarized, and key challenges, such as great conversion efficiency fluctuations due to environmental conditions, low system integration, and the lack of unified evaluation standards, were highlighted. Finally, the technological evolution trends centering on multi-source fusion and intelligent scheduling were discussed, and a design paradigm based on the “intelligent energy internet” was proposed, outlining the future development direction of converting discrete ambient energy into a stable supply through cross-platform collaboration and data-driven optimization.

Energy Storage-Based Co-phase Power Supply System Considering Photovoltaic Integration and Its Cooperative Control Strategy
WANG Wei, HUANG Xiaohong, GUO Kai, QIN Ken, LI Qunzhan, LIU Bingrui, LUAN Xiaoyong
, Available online  , doi: 10.3969/j.issn.0258-2724.20250192
Abstract:

To promote energy conservation and emission reduction, as well as green and low-carbon development of electrified railways, two technical schemes for a photovoltaic-integrated energy storage-based co-phase power supply system were first proposed. Then, the operating conditions of the system were analyzed, and the negative sequence compensation principle was formulated. On this basis, according to the photovoltaic power generation output and lithium battery state of charge (SOC), a two-layer cooperative control strategy was constructed with the control objectives of achieving load peak shaving and valley filling and efficient energy utilization. The upper layer was responsible for energy management and negative sequence control, while the lower layer controlled the converter operation in real time. Finally, through case analysis, the curtailment of photovoltaic generation, charging and discharging, and negative sequence compensation were controlled under comprehensive consideration of multi-dimensional factors, such as load magnitude under traction and regeneration conditions, solar irradiance, and lithium battery charging and discharging. The results indicate that the system exhibits good adaptability and stability under different operating conditions, and the correctness and effectiveness of the proposed control strategy are verified. The photovoltaic integration scheme eliminates 50% of electrical phase separations, realizes controllable negative sequence, promotes multi-energy integration and complementarity, facilitates local and nearby consumption of photovoltaics, and mitigates severe fluctuations of traction load.

Experimental Study on Anchorage Performance of Improved Self-locking Prestressed Anchor Cable
XU Xiwu, LEI Jinshan, YANG Xiuzhu, HONG Jiamin
, Available online  , doi: 10.3969/j.issn.0258-2724.20250594
Abstract:

To address the issues of excessive deformation and prestress loss caused by stress concentration at the external anchor head of traditional prestressed anchor cables in fragmented rock slope reinforcement, an improved self-locking prestressed anchor cable was developed by innovating the clip structure of the self-locking device to enhance anchoring performance. Laboratory tensioning and locking tests were conducted to verify the mechanical locking effect of the improved anchor cable. Subsequently, a typical fragmented rock slope project in Zhangjiajie, Hunan Province, was selected for in-situ tensioning and locking tests. Distributed optical fiber sensing testing technology was employed to monitor the entire tensioning and locking process along the full length of the grout for three improved self-locking tension-type prestressed anchor cables. The axial strain distribution patterns of the grout at various stages were obtained, and the shear stress distribution was derived. Simultaneously, the strain of the steel pipe within the self-locking device was tested and analyzed. The results indicate that during the tensioning process, the inner anchorage section exhibits the typical mechanical characteristics of a tension-type anchor cable. The axial strain gradually develops from the proximal end of the inner anchorage section and continuously propagates backward as the load increases. At this stage, the self-locking device has no influence on the stress state of the inner anchorage section. Within 24 hours after unloading and locking, the strain distribution of the self-locking anchorage section remains stable without significant prestress loss, demonstrating excellent anchoring and locking efficacy. Because the test anchor cable comprises two anchorage units, the interaction between these units causes notable tensile strain, negative shear stress, and an increase in shear stress at the distal end of the MS1 unit. Affected by uncertainties such as on-site construction and measurement, the shear stress distribution is uneven, with a peak value of approximately 0.3 MPa, and the peak position of the MS2 unit is slightly posterior to that of the MS1 unit. In addition, the strain of the self-locking steel pipe starts from the end of the dedicated anchorage and gradually transitions from compressive strain to tensile strain with increasing distance; the strain value increases proportionally with the tensioning load. These findings provide a reference for further optimizing the support structures of fragmented rock slopes and enhancing the anchoring performance of self-locking prestressed anchor cables.

Evaluation of Meta-Action Maintenance Effects for Electromechanical Products Based on IVIF-DANP
RAN Yan, ZENG Xiaotong
, Available online  , doi: 10.3969/j.issn.0258-2724.20250088
Abstract:

The evaluation of maintenance effects is an important reflection of the maintenance support capability of electromechanical products and a key aspect of maintenance management. Traditional methods are limited in handling high-dimensional data and interactions among indicators, making it difficult to accurately reflect maintenance effects. To solve these problems, a new method for evaluating maintenance effects that integrates interval-valued intuitionistic fuzzy sets (IVIFS) and the DEMATEL-based analytic network process (DANP) was proposed. A multi-level maintenance effect evaluation index system was established based on the simple structure and ease of fault diagnosis of electromechanical product meta-action units. In an interval-valued intuitionistic fuzzy environment, the decision-making trial and evaluation laboratory (DEMATEL) method was applied to analyze interdependencies among indicators, and the analytic network process (ANP) was used to determine their weights. Maintenance effects were quantitatively evaluated and ranked by calculating the weighted distances between candidate maintenance plans and the ideal plan. In a case study involving the worm rotation meta-action unit of a CNC grinding machine, three maintenance methods were comprehensively evaluated and ranked. The stability and feasibility of the method were verified by varying the weights assigned to experts and comparing different methods. The results show that the evaluation model integrating IVIF-DANP can accurately differentiate the relative merits of different maintenance plans.

Regulation Control of Flexible Manipulators with Unknown Parameters
XUE Xiaoqiang, WANG Xiaoyong, LIU Wei, CHE Jian, SHEN Xiaolong, CHEN He
, Available online  , doi: 10.3969/j.issn.0258-2724.20250270
Abstract:

With the rapid development of industrial automation and intelligent manufacturing, flexible manipulators have been widely used in the aerospace and precision instrument manufacturing industries. To address the challenges of achieving precise positioning and suppressing vibrations in flexible manipulators with unknown parameters, a sliding mode control method based on a time-delay estimator (TDE) is proposed. First, the dynamic characteristics of the system are analyzed and the dynamic model is established. On this basis, a sliding mode controller is designed, and the unknown system dynamics are estimated by using a time-delay estimator. Furthermore, to enhance the system’s control performance, a term for suppressing the load deflection angle is incorporated into the control law to effectively reduce the effects of flexible vibrations. Subsequently, a rigorous stability analysis of the closed-loop system is conducted, and the convergence of the system state variables is proved. Finally, simulation experiments were conducted to verify the effectiveness of the proposed method. The results show that the proposed control method can achieve accurate positioning while effectively suppressing flexible vibrations.

Numerical Investigation of Flow-Induced Vibration Characteristics Under Typical Cavitation Conditions Around Hydrofoil
CHEN Qian, ZHAO Xiaojun, WANG Yiwen, MENG Fanxin, CHEN Yue, ZHOU Xiaonan, WU Qin
, Available online  , doi: 10.3969/j.issn.0258-2724.20250198
Abstract:

To investigate the unsteady fluid-structure interaction (FSI) characteristics of the NACA66 hydrofoil during typical cavitation stages, its hydrodynamic and flow-induced vibration characteristics were numerically calculated using an unsteady cavitating FSI model. First, an unsteady cavitating FSI numerical model was established, and the calculated results were compared with the experimental results reported by Wu et al. to verify the accuracy of the model. Second, the characteristics of typical cavitation flow pattern distributions and the evolution of the unsteady flow-field structure were analyzed using the numerical method. Finally, the hydrodynamic and flow-induced vibration characteristics during different typical cavitation stages were further investigated. The results show that, during typical cavitation stages, the maximum cavity length increases as the cavitation number decreases, and cloud cavitation exhibits pronounced periodicity. The lift and drag coefficients vary with the cavitation number. The maximum tip vibration displacement and velocity of the hydrofoil remain relatively high during incipient cavitation but remain low during supercavitation due to the effect of full cavity coverage. During cloud cavitation, the maximum vibration displacement and velocity exhibit substantial periodic fluctuations due to unsteady flow effects. The maximum tip displacement and maximum tip torsion angle of the hydrofoil are linearly related to the lift and pitching moment, respectively. In addition, the flow-induced vibration frequency of the hydrofoil is approximately equal to the cavitation shedding frequency, and the vibration frequency decreases as the cavitation number decreases.

Connection Scheme of Stray Current Collection Network for DC Traction Power Supply System
LIU Feilong, LIU Wei, YANG Shuangrui
, Available online  , doi: 10.3969/j.issn.0258-2724.20260049
Abstract:

To address DC interference in surrounding power systems and damage to urban pipeline networks caused by stray currents, a novel optimization scheme for the stray current collection network is proposed, in which a continuous cable is used to connect the trackbed reinforcement. The underlying mechanism through which this scheme influences the effectiveness of stray current protection is also systematically investigated. Based on the proposed continuous cable connection scheme, a four-layer return current model of a DC traction power supply system was developed to account for the stray current collection network. The boundary element method was then employed to couple the leakage current from the stray current collection network with the return current system model, thus improving simulation accuracy. The simulation accuracy of the proposed model was verified by scaled simulation experiments. A metro project in China was used as a case study to quantitatively analyze the effects of the conventional cross-connection method using 95 mm2 cables and the continuous cable connection method on secondary stray current and the earth potential gradient. The results show that the average simulation error is 11.76%. When a 95 mm2 cable is used directly for cross-connection and the rail-to-earth transition resistance is less than 7.72 Ω•km, the earth potential gradient along the entire line remains below 2.5 mV/m. When a 120 mm2 continuous cable is used for connection and the rail-to-earth transition resistance is less than 5.37 Ω•km, the earth potential gradient along the entire line remains below 2.5 mV/m. Using a continuous cable connection can significantly reduce stray current. Using a 120 mm2 continuous cable for connection reduces the secondary stray current by 36.59% and the earth potential gradient by 35.28% overall. Using a 200 mm2 continuous cable for connection reduces the secondary stray current by 47.30% and the earth potential gradient by 46.89% overall.

Metro Stray Current Leakage Zone Localization Method Considering Rail Resistance Uncertainty
WANG Aimin, ZHOU Wenjie, YANG Xin
, Available online  , doi: 10.3969/j.issn.0258-2724.20260086
Abstract:

Existing metro stray-current leakage location methods often treat longitudinal rail resistance as deterministic and ignore its uncertainty caused by material aging, joint resistance, and environmental factors, which reduces location accuracy and generalization performance. Based on measured longitudinal rail resistance data from a metro system, a method that accounts for the probabilistic characteristics of longitudinal rail resistance was proposed for locating stray-current leakage. Kernel density estimation was used to establish the probability distribution model of longitudinal rail resistance. Monte Carlo methods were then used to simulate uncertainty in longitudinal rail resistance, and a rail-earth resistance network model was constructed. Building on this model, a method was proposed to calculate the distribution of rail-to-earth transition resistance or conductance using a hybrid genetic algorithm-particle swarm optimization algorithm, and a stepwise strategy was developed for locating stray-current leakage. The strategy first identified the inter-station section containing the leakage position through coarse localization and estimated the distribution of longitudinal rail resistance by selecting the solutions with the best fitness values. The leakage zone within the identified section was then determined through stepwise subdivision. During the location process, the leakage position was identified by calculating the probability that rail-to-earth transition conductance along the line exceeds a threshold. Comparisons of the location results with simulated leakage zones show that the proposed method can accurately identify the inter-station section containing the leakage zone when longitudinal rail resistance is nonuniformly distributed and leakage severity and range vary. The location deviation is less than 30%, and the average relative error in estimating longitudinal rail resistance is less than 5% in the presence of measurement noise, indicating strong robustness and practical engineering applicability.

Optimal Configuration Scheme of Active Stray Current Mitigation Devices in Rail Transit
WANG Miao, MA Wenzhong, XU Hailiang, GE Pingjuan, YANG Xiaofeng, ZHENG Qionglin
, Available online  , doi: 10.3969/j.issn.0258-2724.20260062
Abstract:

Active stray current mitigation devices show considerable potential for application in rail transit systems. However, the arrangement of these mitigation devices within power supply sections cannot effectively accommodate dynamic changes in stray current distribution, making it difficult to maintain effective stray current mitigation under peak current conditions during train traction or braking. Therefore, an optimal configuration strategy for these active stray current mitigation devices within power supply sections was proposed, using a DC autotransformer (DCAT)-based traction power supply system as an example. First, the rail transit traction power supply system model was established, and the suppression ratios of stray current and rail potential were derived to reveal the influence of these active devices on both stray current and rail potential distributions. Second, based on train operating characteristics, the optimal configuration principles incorporating offset enhancement, symmetric distribution, and peak response were proposed. The optimal installation positions of these active devices were determined by defining the optimized offset ratio, and the mitigation performance for stray current and rail potential in the system was compared before and after configuration optimization. Finally, the Beijing Metro Yizhuang Line was selected as a case study to validate the proposed strategy through simulations. The results show that the optimal configuration strategy significantly improves the mitigation performance without increasing the number of active devices. The error between simulation results and theoretical calculations is controlled within 3%. After configuration optimization, under traction and braking conditions, the maximum rail potential at the traction substation drops from 19.4% and 24.5% to 5.9% and 8.1%, respectively, while the maximum stray current drops from 53.1% and 66.5% to 19.7% and 26.6%, respectively, relative to the corresponding values in the conventional system.

Detection Method for Misrecognition and Missed Captures by License Plate Recognition Devices Based on Trip Chains and Neighbor Entropy
HAN Ke, SHAN Le, PENG Zhiheng, ZHANG Lei, BAI Yun, HUANG Tao, LIU Shijia
, Available online  , doi: 10.3969/j.issn.0258-2724.20260061
Abstract:

To address the vulnerability of license plate recognition (LPR) data in intelligent transportation systems (ITS) to implicit errors, such as character misrecognition and missed vehicle captures, and the lack of methods for proactively identifying these errors, this study aims to develop a detection method that can proactively and accurately identify LPR devices exhibiting abnormal behavior, thus improving overall data quality and providing decision support for equipment operation and maintenance. First, topological neighborhood relationships between devices were defined using path data derived from vehicle trip chains connecting consecutive LPR checkpoints in the road network, and multidimensional neighbor entropy feature indicators, including total entropy, upstream and downstream neighbor entropy, and asymmetric entropy, were constructed accordingly. Second, a candidate set of abnormal devices suspected of misrecognition or missed captures was initially identified by detecting distortions in neighborhood distributions and fluctuations in vehicle counts caused by implicit errors. Then, spatiotemporal constraints were introduced for cross-verification. Misrecognition devices were identified using physical speed limits, while missed capture devices were identified by incorporating road network topology and frequency-direction consistency constraints. Finally, anomaly confidence scores were calculated for each LPR device, and the final sets of misrecognition and missed capture devices were obtained. Multi-scenario simulation experiments conducted on an actual urban road network show that, across various complex scenarios with different misrecognition and missed capture probabilities and abnormal device counts, the precision of missed capture device identification exceeds 0.8, the recall of misrecognition device identification reaches 1.0, and stable and satisfactory identification performance is maintained under mixed anomaly conditions where both errors coexist.

Dual-Parameter Compensation Method for Temperature Effects in Large-Scale Fiber Optic Balances
XU Bin, LIU Long, YU Shien, ZHANG Jianzhong
, Available online  , doi: 10.3969/j.issn.0258-2724.20250650
Abstract:

A dual-parameter compensation method based on thermo-mechanical coupling mechanisms was proposed to suppress the thermally induced drift of large-scale fiber optic balances caused by non-uniform temperature fields in long-duration wind tunnel tests and improve aerodynamic force measurement accuracy. Theoretical modeling was performed to reveal the synergistic deformation mechanism of elastic elements and floating frames under the coupling action of thermal expansion and structural constraints in non-uniform temperature fields. Constraint forces neglected by traditional models were identified as the key factor limiting compensation accuracy. Finite element simulation was adopted to analyze the differential thermal response laws of different measurement channels in uniform temperature rise and temperature difference conditions. A compensation model based on the temperature and temperature difference of symmetrically arranged elastic elements was further built, and the directly measurable dual parameters were employed to indirectly characterize internal constraint force effects. Static and dynamic temperature calibrations of fiber optic balances were conducted to obtain temperature sensitivity coefficients of each balance component, and a temperature history correction formula for elastic elements was fitted to eliminate response delay induced by thermal inertia. Thirty-second long-duration force measurement tests were conducted in a blowdown wind tunnel to validate the compensation method. Test results show that the aerodynamic force measurement errors of three measured components (axial force, normal force, and pitching moment) caused by thermal drift decrease from 18.87%, 85.37%, and 11.67% before compensation to 3.19%, 5.63%, and 1.86%, respectively. The results verify the high-accuracy compensation capability of this method under non-uniform temperature fields, providing a practically applicable temperature compensation scheme for long-duration wind tunnel force measurement tests.

Demonstration of Equivalent Grade Sections for Track Longitudinal Profile on Kilometer-Scale High-Speed Railway Bridges During Operation
LYU Tao, WANG Jinxin, TAN Shehui, WANG Weiqin, XU Weichang, CHEN Rong, WANG Ping
, Available online  , doi: 10.3969/j.issn.0258-2724.20250281
Abstract:

The concept of equivalent grade sections was proposed and its reasonable equivalent method was investigated to facilitate the representation and maintenance management of track longitudinal profiles on kilometer-scale high-speed railway bridges during operation. Additionally, the objective includes addressing following problems. There is large deviation between actual and designed longitudinal profiles caused by manufacturing deviations, construction errors, and limited ballast adjustment capacity, which necessitates modification to a continuous short-slope configuration. Temperature loads during operation induce bridge deflections larger at mid-span and smaller at both ends. The actual longitudinal profile evolves into a high-order spatial curve, posing challenges to track maintenance management. By taking the tracks on two in-service kilometer-scale high-speed railway bridges as examples, four equivalent grade sections—the single-grade section, two-grade section, multi-grade section, and 200 m continuous grade section—were first established by comparing grade section lengths and mid-span elevations, simplifying the actual high-order spatial curve profile into a conventional combination of grade sections and vertical curves. Then, key indicators such as acceleration and ride comfort in actual and equivalent longitudinal profile conditions were calculated and compared. Finally, considering the need for simplified track management, the reasonable equivalent method was determined. The results show that for the Wufengshan Yangtze River Bridge, the maximum acceleration deviation between the single-grade equivalent and measured profiles is no more than 0.2 m/s2, and the ratio of resistance difference to tractive force is no more than 15%. For the Husutong Yangtze River Bridge, the maximum acceleration deviation is no more than 0.1 m/s2, and the ratio of resistance difference to tractive force is no more than 8%. The single-grade equivalent for both bridges is basically equivalent to the measured profile and provides a more intuitive representation. Considering train traction capacity, ride comfort, and the need for simplified on-site management, the single-grade equivalent is recommended as the basis for track management during operation.

Two-Stage Optimization Dispatching Strategy for Distribution Network Considering Dispatchable Potential of Charging Station
JIANG Xiaofeng, PAN Pengyu, ZHOU Bo, JIA Shicheng, TANG Chun, XU Yunyang, WEI Wei, YANG Jianwei
, Available online  , doi: 10.3969/j.issn.0258-2724.20240420
Abstract:

To address the challenges posed to the safe and economic operation of distribution networks by the large-scale integration of electric vehicles (EVs) and wind and solar power, as well as the issue of low accuracy in traditional single timescale optimization dispatching methods, a multi-timescale coordinated distribution network optimization dispatching strategy that considered the dispatchable potential of charging stations was proposed. Firstly, an assessment method for the dispatchable capacity of charging stations based on the Minkowski sum was proposed. By aggregating EV clusters within the stations and treating them as generalized energy storage systems, the real-time dispatchable potential of each station was accurately calculated. Subsequently, a day-ahead and intra-day two-stage optimization dispatching strategy for the distribution network based on model predictive control was introduced. In the day-ahead dispatching stage, with the objective of minimizing the daily operating cost of the distribution network, the electricity purchase and sales plans, along with the output schedules of controllable distributed generators and energy storage systems, were coordinately optimized to achieve optimal allocation of global resources. In the intra-day dispatching stage, by leveraging rolling optimization and feedback correction mechanisms based on the day-ahead plan and real-time states, the day-ahead output schedules of controllable devices were dynamically corrected to cope with intra-day fluctuations in wind power, photovoltaic power, and load, ensuring the stable operation of system state variables under disturbances. Finally, the performance of the proposed strategy in fulfilling day-ahead dispatching instructions under uncertain disturbances was analyzed. Simulation results demonstrate that by exploiting the dispatchable potential of charging stations, the proposed strategy reduces the daily operating cost of the distribution network by 10.86% compared to scenarios without V2G capability. Furthermore, it achieves an additional 13.66% cost reduction compared to day-ahead optimization strategies, while limiting the deviation from the day-ahead tie-line power to merely 1.77%. This approach realizes distributed optimization dispatching for the distribution network, significantly enhancing its economic efficiency and robustness amidst uncertainties.

Theoretical Calculation Methods and Field Test Optimization of Surrounding Rock Load for Shield Tunnels in Mudstone Strata
ZHOU Ziyang, ZHANG Jingxuan, FENG Kun, LU Chunfang, HE Chuan, Ma Gaoyu
, Available online  , doi: 10.3969/j.issn.0258-2724.20240395
Abstract:

To address the challenges of unclear distribution patterns and inaccurate calculation of surrounding rock load for shield tunnels within rock stratum environments, field monitoring tests were conducted based on the Jinxiu Tunnel of the Chengdu–Zigong Railway. The temporal variation pattern of water and soil load on the external segments during the construction process was analyzed, and the distribution patterns of tunnel surrounding rock loads in a mudstone environment were summarized. The applicability of commonly used calculation theories for rock stratum load was explored, and the calculation method of the surrounding rock load was improved by introducing a lateral gradient coefficient of the surrounding rock load. The results show that during the segment assembly phase, the external load distribution is uniform, with field test results indicating that the external load on each segment is within 60 kPa. Upon stabilization, the external load generally attenuates to approximately 80% of the peak load during the construction process. In the weathered mudstone strata of the tunnel site area, the distribution pattern of the surrounding rock load more closely resembles that of loose pressure distribution, and it is more appropriate to account for the influence of groundwater load by adopting the combined water-soil calculation method. Inversion of field test results suggests that the equivalent rock column height for tunnels in the weakly weathered mudstone strata of Chengdu is 1.13 times the tunnel diameter. After optimizing the calculation of lateral rock pressure in the rock layer, the root mean square error (RMSE) compared to field monitoring data is below 36.4 kPa, indicating a significant improvement in calculation accuracy.

Dynamic Response and Deformation and Failure Modes of Reservoir Bank Deposit Slopes Based on Shaking Table Tests
REN Huajiang, GUO Tao, PAN Yuanyang, LIU Jiangfan, ZHANG Duoduo
, Available online  , doi: 10.3969/j.issn.0258-2724.20250256
Abstract:

Shaking table model tests were conducted on the Bailong deposit bank slope of the Lianghekou Hydropower Station to reveal the dynamic instability mechanism of reservoir bank deposit slopes under the combined action of impoundment and earthquakes. The dynamic response, energy evolution, and whole deformation-failure process of the bank slope in coupling conditions were analyzed. Wenchuan seismic waves with step-increased peak acceleration were adopted for excitation tests. Acceleration and pore water pressure monitoring, together with Hilbert-Huang transform time-frequency analysis, were applied to identify the internal damage evolution of the slope mass. The results show that significant local amplification occurs near the water level, besides the dynamic amplification at the slope crest in impoundment conditions. Low-frequency energy near the water level is obviously higher than that in the middle and upper slopes, and impoundment aggravates the accumulation of low-frequency energy in this zone. Pore water pressure in the weak interlayer rises sharply when seismic peak acceleration exceeds 0.5 g, which reduces the effective stress of the soil mass. Marginal spectrum analysis indicates that slope damage is initiated near the water level and slope crest, and propagates upward with the increasing earthquake strength. Under the combined action of impoundment and earthquakes, deformation and failure first take place near the water level. Continuous damage develops in the upper-slope mass, and traction-type sliding failure eventually occurs along the weak interlayer. The slope evolution includes four successive stages: erosion deformation, local damage accumulation, local deformation and failure, and overall sliding instability. The findings provide experimental support for seismic stability evaluation of similar deposit bank slopes in southwestern reservoir areas.

Principle for Clause Redundancy in Clause Sets with Equality
NING Xinran, WU Guanfeng, XU Yang
, Available online  , doi: 10.3969/j.issn.0258-2724.20240586
Abstract:

To investigate clause redundancy in clause sets with equality, a determination method for clause redundancy, namely equality implication modulo resolution (EIMR) principle, was proposed. First, the need for handling redundant clauses during the simplification process in a first-order logic theorem prover was analyzed, and the limitations of implication modulo resolution principle in handling clause sets with equality were identified. By integrating flattening and resolution methods, the EIMR principle applicable to clause sets with equality was defined. The reliability of the principle was proven, and the specific role of the principle in the determination of redundant clauses was clarified. Subsequently, on the basis of EIMR, the clause elimination methods from propositional logic and clause sets without equality were extended to clause sets with equality. Asymmetric tautologies and subsumption clauses of equality resolution were defined, and the effectiveness of these methods was demonstrated. Finally, the reliability of predicate elimination preprocessing methods in mainstream first-order theorem provers was validated using the EIMR principle. The results indicate that EIMR provides a theoretical basis for clause redundancy studies in clause sets with equality and extends the applicable scope of existing clause elimination methods.

Collaborative Braking Control of Virtual Coupling Trains Based on Soft Actor-Critic Algorithm
LIN Junting, QIU Xiaohui, LI Maolin
, Available online  , doi: 10.3969/j.issn.0258-2724.20240312
Abstract:

Traditional train operation control systems usually adopt a maximum braking control strategy during emergency braking; however, when this strategy is applied to virtual coupling, it may lead to collisions among train units. Considering the heterogeneity of virtual coupling train units, different braking rates of train units, and different safe spacings, a collaborative braking control method based on the Soft Actor-Critic (SAC) algorithm for virtual coupling trains in emergency braking scenarios is designed and analyzed. First, a train safety braking model in a relative coordinate system is constructed based on the running speed of the preceding train within the virtual coupling, the spacing between adjacent trains, and the train's own braking performance, and the safe speed limit of the train is obtained. Second, an operating environment is built according to the train safety braking model, and a braking control model is designed using a multi-train SAC algorithm. The state space and action space of each train unit are determined, and the reward function is designed with the safe speed limit generated by the safety braking model in the relative coordinate system as a constraint. Finally, taking the Beijing Capital Airport Express line and its vehicles as an example, the multi-train SAC algorithm is validated and compared with the maximum braking control strategy and the minimum braking control strategy. The results show that compared with the other two traditional control strategies, the collaborative braking model can not only stop safely before the hazardous point but also ensure spacing safety during braking, thereby improving the operational safety and efficiency of the urban rail transit system.

Icing and Deicing Jump Analysis of Transmission Tower-Line System under Heavy Ice Load
ZHANG Xin, XIE Qiang, LI Yue, YAN Cong
, Available online  , doi: 10.3969/j.issn.0258-2724.20240405
Abstract:
Objective

With the increasing demand for electricity, the transmission lines inevitably need to traverse steep mountainous terrains with significant elevation differences. Traditional high-low leg transmission towers are limited by the angle constraints of their main and diagonal elements. To overcome these limitations, new tower designs with transition sections have emerged in recent years, solving the problem of positioning the towers in terrains with large elevation differences. These new tower types have been applied to the transmission lines of 220 kV and 500 kV voltage levels. However, research on the ±800 kV ultra-high-voltage (UHV) transmission lines under heavy ice load remains a research gap. Southwestern China is rich in hydropower resources, but the natural conditions of the region, steep terrain, and severe icing, present significant challenges, with many areas classified as moderate to heavy ice regions. With the implementation of the west-to-east power transmission strategy of China, this new type of tower design with transition sections has promising application prospects for transmission corridors in this region. However, the failures of transmission lines caused by tower collapses and line breakages under ice load frequently occur in mountainous areas. Therefore, the structural and electrical safety of the new tower type, especially when applied to UHV lines under icing and ice-shedding conditions, urgently needs investigation.

Methods

A steep-span tension section of the ±800 kV UHV transmission line from Sichuan to Jiangsu was studied. The research employed a catenary model to simulate the shape of ground lines and conductor lines, and a finite element model of the transmission tower-line coupling system with transition sections was developed using ABAQUS software. The model accounts for line angles, as well as elevation differences; moreover, the structural dimensions and material property parameters were selected based on actual engineering data. The structural and electrical safety characteristics of the transmission tower-line coupling system under different icing and ice-shedding conditions were analyzed. Specifically, the stress distribution patterns of main and diagonal elements along the height of the tower during icing and ice-shedding were investigated, and the tension and safety factors of ground and conductor lines were evaluated; the displacement response and energy transfer mechanisms of ice-shedding were explored.

Results

Results indicate that under the heavy icing condition of 30 mm, the stress ratios of the main and diagonal elements of the transmission tower are less than 0.4 and 0.3, respectively. The equivalent safety factors of the conductor and ground lines exceed 2.25, and the unbalanced tension complies with the code requirements, indicating a significant safety margin for the tower-line system. Compared to a uniform ice load condition, a non-uniform ice load condition significantly increases the stress on elements of the transmission tower, posing greater risks to the stability of the main and diagonal elements in transition section towers. However, the stresses remain within safe limits, and the main and diagonal elements remain stable. The weak positions of the transmission tower-line system with transition section include the main elements in the tower legs, the transition section and the cross-arms, and the diagonal elements of the ground line supports, which require special attention to the stability under heavy ice loads. The stress ratios along the height of the tower under non-uniform ice load conditions exhibit similar patterns to those under uniform ice load conditions, albeit with higher values. During ice-shedding under a 30 mm heavy ice load condition, the main and diagonal elements of the ±800 kV UHV transmission tower-line system with transition sections remain stable. The tensions of the ground and conductor lines all meet the minimum allowable safety factor, with only a few points slightly below the design safety factor value, maintaining adequate safety margins. The maximum vertical displacements of the ground and conductor lines of the ±800 kV UHV transmission tower-line system with transition sections during ice-shedding are 7.9 m and 5.2 m, respectively, which satisfy the requirements for air insulation clearances. However, for higher voltage levels or multi-phase transmission lines, electrical safety issues such as flashover induced by ice-shedding should be carefully considered. The ice-shedding process of ground and conductor lines exhibits typical damped cosine wave characteristics, with the maximum vertical displacement of a single line occurring in the first cycle, followed by gradual attenuation. Furthermore, when multiple lines unload ice simultaneously, the maximum jump height may not occur in the first cycle and may exceed the maximum displacement value of a single line, which requires significant attention. This indicates that energy transfer occurs between different lines during ice-shedding in the transmission tower-line coupling system. The jump height of a single line is influenced not only by its own response but also by energy transfer with other ground and conductor lines in the same span.

Conclusion

In conclusion, the coupling effects between transmission towers and transmission lines, as well as the complex dynamic response mechanisms induced by ice-shedding of multiple ground and conductor lines, should be fully considered in the analysis of transmission lines under heavy ice loads. The weak position of the transmission tower-line system with transition section is the main elements in the tower legs, transition section and cross-arms, and the diagonal elements of the ground line supports. Although the structure maintains adequate safety margins, it is essential to prevent flashover hazards and element instability to ensure the long-term reliability and safety of transmission lines in heavy ice regions.

High-Efficiency Miniaturized Dipole Antenna for Single-Sided Metal Wall Environment Inside Inspection Vehicle
QIU Peiyi
, Available online  , doi: 10.3969/j.issn.0258-2724.20250536
Abstract:

To achieve high aperture efficiency and miniaturization for a dipole antenna operating in an environment with a single-sided metal wall inside an inspection vehicle, this paper analyzes and optimizes the antenna’s performance using multi-arm folding and current control methods. First, based on the effective current length theory, a four-arm folded dipole antenna fed by a vertical microstrip line is proposed. While maintaining the effective current length, the radiating arms are folded horizontally and vertically, respectively, achieving antenna miniaturization. Second, a back cavity is added beneath the antenna, and a horizontal strip is introduced to concentrate the resonant current in the radiation plane, thus improving the aperture efficiency and gain of the four-arm dipole antenna. Finally, to address the impedance mismatch and radiation pattern distortion caused by the single-sided metal wall inside the inspection vehicle, the four-arm folded dipole antenna is further redesigned. A shorting patch is introduced to suppress surface currents on the side nearest the wall, while the vertical arm nearest the wall is extended to suppress undesired current excitation, thus mitigating the effects of the metal wall on the antenna’s impedance and radiation pattern performance. Both antennas have been fabricated and measured. The results show that the antenna size is reduced from a half-wavelength to approximately 1/4 wavelength, achieving a reduction of 50%. The measured bandwidth is 410–447 MHz with good impedance matching. The antenna can effectively suppress the effects of the single-sided metal wall, with a peak gain of 5.5 dBi, cross-polarization levels below −10 dB, backlobe levels below −10 dB, and radiation efficiency above 80%, meeting the requirements of the wireless transmission system inside the inspection vehicle.

Resource Allocation Algorithm for Reconfigurable Intelligent Surface-Assisted Device-to-Device Communication Network
XIE Jianli, LI Lin, ZHANG Zepeng, LI Cuiran
, Available online  , doi: 10.3969/j.issn.0258-2724.20240278
Abstract:

To address the severe interference problem caused by non-orthogonal multiple access (NOMA) and device-to-device (D2D) technologies in enhancing the capacity of heterogeneous cellular networks, an efficient resource allocation algorithm was proposed to improve the system and rate. Firstly, a reconfigurable intelligent surface (RIS)-assisted heterogeneous cellular NOMA-D2D communication network model was constructed. Under the constraints of signal-to-noise ratio (SNR) of users, transmission power, and unit membrane of RIS phase shift, an optimization problem aiming to maximize the number of D2D users and rate was established. However, this problem was a mixed integer nonlinear problem and was difficult to solve directly. Therefore, it was decomposed into three sub-problems: D2D user-cellular user (CU) matching, D2D user power control, and RIS reflection phase shift optimization. Based on this, the D2D cluster channel allocation was completed by using the bipartite graph maximum matching algorithm, and a deep reinforcement learning algorithm based on a deep deterministic gradient strategy (DDPG) was proposed to jointly optimize the D2D transmission power and RIS phase shift matrix. Simulation results show that under the same conditions, the average value of the D2D link and rate of the proposed algorithm increases by 1.96%, 10.64%, and 14.29%, respectively, compared with that of the game algorithm, random phase shift algorithm, and RIS-assisted-free algorithm, verifying its effectiveness in interference suppression and spectral efficiency improvement.

Characterization of Dynamic Heat Transfer Characteristics of Building Envelopes Based on Energy Attenuation Coefficient
HUANG Haotian, ZHANG Xiaojing, XIE Jingchao, MA Anke
, Available online  , doi: 10.3969/j.issn.0258-2724.20250509
Abstract:

An energy attenuation coefficient (EAC) is introduced to quantify the dynamic response of building envelopes to outdoor thermal disturbances. It is defined as the ratio of the time integral of the absolute heat flux at the internal surface to that at the external surface over a given period. A higher value indicates weaker attenuation, whereas a lower value indicates stronger mitigation of outdoor thermal disturbances. Yanqi, Xinjiang, characterised by large air-temperature fluctuations and strong solar radiation, is used as a case study. Heat transfer in 25 typical envelope assemblies was simulated using a combined frequency- and time-domain method, and layer-specific EAC values were evaluated at daily, monthly, seasonal, and annual scales. The results show that the annual mean EAC values for the external insulation layer, concrete layer, and internal insulation layer are 0.85, 0.40, and 0.99, respectively, indicating that the concrete layer is the primary contributor to energy attenuation. During the heating season, the EAC of the concrete layer increases significantly to 0.78 compared with other typical periods, indicating that the building envelope is less effective at attenuating outdoor thermal disturbances. The external insulation layer mainly provides steady-state thermal resistance. At the annual scale, moderately increasing the thermal conductivity and heat capacity of the insulation layer improves the overall attenuation of outdoor thermal disturbances throughout the year.

Green Power Utilization and Zero-Carbon Feed-Through Power Supply Technology for Electrified Railways
LI Qunzhan, HUANG Xiaohong, WU Bo, XIE Shaofeng
, Available online  , doi: 10.3969/j.issn.0258-2724.20250054
Abstract:

Electrified railways rely heavily on public power grids dominated by fossil fuels, and actively carrying out the development and utilization of renewable energy in rail transit is one of the important measures to achieve China’s “carbon peaking and carbon neutrality” goals. To extensively and comprehensively use green power to construct electrified railways with zero-carbon emission, a technical scheme of zero-carbon feed-through power supply and its zero-carbon operation control strategy were proposed. First, the basic methods of new energy connection in electrified railways were analyzed, and the connection points and energy flow paths of renewable energy in the traction power supply system were clarified. On this basis, a technical framework for green power utilization based on feed-through power supply was constructed, and the natural flow of train traction and regenerative power was realized by using the feed-through traction network, thereby reducing the capacities and investment of green power and energy storage equipment. Second, a three-component power balance control strategy covering green power, traction, and energy storage was proposed. Green power was controlled according to maximum power point tracking to maximize on-site power generation; the energy storage devices in the traction substation were controlled for charging and discharging according to the active power specified in the green power purchase/sale agreement with the power grid, thereby achieving zero-carbon operation. Finally, by taking an actual railway line as an example, the effectiveness and economy of the proposed scheme were verified. The results indicate that the zero-carbon feed-through power supply system does not obtain fossil power from the power grid and causes no impact on the power grid; when the active power exchanged with the power grid is zero, the negative-sequence power and through power are also zero. The feed-through power supply traction network eliminates electrical phase separation, and the regenerative braking energy and green power are fully utilized. In the verified case, the total regenerative braking energy of the traction substation accounts for only 0.6% of the traction energy and is all effectively utilized, with a photovoltaic curtailment rate of less than 10%. Green power and energy storage adopt local independent control, which significantly reduces the complexity of communication and power flow control. Based on capacity configuration and economic analysis using actual data of heavy-haul railways, the zero-carbon power supply requirements are met when lithium iron phosphate batteries with a charge/discharge rate of 0.5 C and an energy storage configuration capacity of 145 MW•h are adopted. According to the annual average irradiance and the irradiance under extreme weather conditions, the installed photovoltaic capacities of the entire line need to be configured as 45 MWp and 83 MWp, respectively, with cost recovery periods of about 5.9 years and 6.8 years, respectively.

NGO-Based CNN-BiLSTM-AM Model for Landslide Displacement Prediction
WANG Huiqin, GUO Ruili, HE Yongqiang, LIU Bincan
, Available online  , doi: 10.3969/j.issn.0258-2724.20240550
Abstract:

A convolutional-bidirectional long short-term memory neural network-attention mechanism (CNN-BiLSTM-AM) prediction model optimized by the northern goshawk optimization (NGO) algorithm for landslide displacement was proposed to address challenges that a single prediction model fails to effectively extract complex sequence features and that manual parameter tuning tends to fall into local optima in current landslide displacement prediction research. Firstly, according to the factors affecting the landslide, the multivariate empirical mode decomposition (MEMD) algorithm was used to decompose various landslide displacement data into trend and periodic components. The trend components were predicted using the autoregressive integrated moving average (ARIMA) method. For the periodic components, influencing factors were identified through the gray correlation degree, and a CNN-BiLSTM-AM combined model was constructed for prediction. The optimal hyperparameters of this model were obtained through NGO. Then, by considering the lag of the periodic components, the Spearman correlation coefficient was used to select the optimal lagged displacement to further enhance the model’s predictive performance. Finally, the model was validated using monitoring data of the Tuojiashan Landslide in Weiyuan, Gansu Province. The results show that the root mean square error (RMSE) and mean absolute error (MAE) of the total displacement prediction of the Tuojiashan landslide are as low as 0.22 mm and 0.37 mm, respectively, showing the prediction accuracy of the correction, while the R2 reaches 0.98, which fully verifies the validity and reliability of the proposed model in landslide displacement prediction.

Text-to-Image Generation Method Based on Fine-Grained Token Competitive Attention
LIU Yi, ZHOU Jing, YUAN Ting, WANG Qian
, Available online  , doi: 10.3969/j.issn.0258-2724.20250670
Abstract:

In multi-attribute text-to-image generation, maintaining object integrity and semantic consistency in complex scenes remains challenging, primarily due to competition and interference among attributes that can lead to structural incompleteness and semantic mismatches. Existing diffusion models struggle with multiple objects and their attribute descriptions due to insufficient prompt information, semantic entanglement, and uneven attention allocation, leading to missing objects, attribute confusion, and misaligned relations. To address this, a multi-stage generation method based on fine-grained token-level competitive attention regularization is proposed, enabling precise control via prompt enhancement, semantic decoupling, and attention optimization. First, in the prompt stage, an LLM enhances the semantics of the input text and structures its representation to reduce ambiguity. Next, in the encoding stage, an attribute-object-level semantic decoupling mechanism is designed, integrating a token encoding module, positive and negative semantic vectors, and a semantic neighbor strategy to improve the stability and discriminability of embeddings while reducing semantic interference. Finally, in the generation stage, dual competitive attention regularization enhances the responses of attributes and objects in cross-modal attention while suppressing attention overlap between different objects, thus improving object independence and reducing semantic interference. Experiments on the T2I-benchmark show significant improvements in color, shape, and texture attribute binding, with bootstrapping language-image pre-training visual question answering (BLIP-VQA) scores of 0.7392, 0.6205, and 0.5013, respectively. Visualizations further confirm that the method effectively reduces attribute mixing and semantic misalignment in complex multi-object scenes while preserving object integrity and semantic consistency, demonstrating good robustness and controllability.

Optimization Model for Rolling Stock Operation Scheduling in Urban Rail Transit Under Flexible Maintenance Strategy
CHEN Zhexuan, CHEN Shaokuan, FENG Jia, MI Qi, BAI Yun
, Available online  , doi: 10.3969/j.issn.0258-2724.20250443
Abstract:

To address the conflict between high-intensity operation and frequent maintenance of urban rail transit trains under the traditional maintenance strategy, an optimization method for multi-day operation scheduling of urban rail transit rolling stock considering a flexible maintenance strategy was studied under a decentralized execution mode of maintenance tasks. Firstly, an optimization model was constructed to minimize the sum of the fixed maintenance cost, the over-maintenance cost, and the insufficient-maintenance cost of rolling stock caused by unbalanced maintenance, considering constraints such as transportation task execution, maintenance task execution, depot maintenance capacity, rolling stock maintenance demand, and standby train arrangement. Secondly, in view of the characteristics of the multi-day operation plan, such as large scale, complex constraints, and strong correlation, a rolling stock space-time connection network with transportation tasks and maintenance tasks as nodes was established, and an improved depth-first search algorithm was designed to solve the constructed model. Finally, a case study was conducted taking an urban rail transit line as an example. The results indicate that compared with the traditional maintenance strategy, the flexible maintenance strategy reduces the over-maintenance and insufficient-maintenance mileages of rolling stock by 59.42% and 49.75% on average, respectively, decreases the total maintenance cost of rolling stock by 6.60% on average under scenarios with different numbers of operating rolling stock, and saves the minimum number of rolling stock required for the operation plan by one train; under the scenario with the same number of operating rolling stock, the decentralized execution of maintenance tasks reduces the standard deviation of the operation mileage of all rolling stock by 7.08% and reduces the number of maintenance rolling stock by two trains; the insufficient-maintenance cost of rolling stock under the same operation mileage level increases with the decrease of the failure rate parameter, and thus the rolling stock tends to execute maintenance before reaching the maintenance mileage standard; the proposed method can flexibly adjust the maintenance locations and timings of rolling stock under the conditions of limited depot maintenance capacity and increased transportation tasks during holidays, which verifies the robustness of the model in extreme scenarios.

Visual Localization Method Research of Fatigue Crack in Steel Box Girder Based on Registration Assisted by Weld
GAO Tian, JI Bohai, YUAN Zhouzhiyuan, LIU Yuhang, SUN Hongbin
, Available online  , doi: 10.3969/j.issn.0258-2724.20240432
Abstract:

To address the localization problem of fatigue cracks in steel box girders with complex structures, a registration algorithm for crack and structure was proposed by leveraging the linear geometric features of welds. Consequently, the whole perception process for identifying, quantifying, and localizing the cracks was integrated and automated. Firstly, a dual-task deep neural network, CrackFormer-Hough, was established to identify the cracks and welds in two-dimensional (2D) images simultaneously. Then, the identified targets were reconstructed by stereo vision techniques to generate target point clouds and quantify their geometric features. Finally, by taking welds as the reference, the coordinate transformation relationship between the camera and structure space was obtained based on the coarse/fine point cloud registration algorithms. Based on structural coordinates, the crack segments from each scope were spliced and fused, and the complete crack morphology was plotted in the structural map. Field tests show that the identification network can overcome background interference and accurately distinguish the weld structure of specific regions by introducing the deep Hough transform to learn the scene semantics of steel box girders. By means of the obvious linear features of welds, the registration errors are controlled to about 1 mm with sparse point clouds. The fatigue crack perception performance is not sensitive to the camera poses. Under different lenses, the average localization error is within 5 mm, and the length measurement error is only about 1%.

Online Identification of Ship Motion in Different Maneuvering Conditions
MENG Yao, ZHANG Xianku, ZHANG Xiufeng, DUAN Yating
, Available online  , doi: 10.3969/j.issn.0258-2724.20240320
Abstract:

An online non-parametric model for ship motion applicable to various maneuvering conditions was developed to enhance the accuracy of ship motion identification modeling and the autonomy and safety of ships during navigation. Firstly, given the complexity of ship navigation characteristics in different maneuvering conditions and challenges of online non-parametric identification, an adaptively updated ship motion non-parametric identification method was proposed by combining the sliding time window and relevance vector machine (RVM). Secondly, the effectiveness of offline non-parametric identification models based on RVM was validated via two different training sample selection schemes, with the importance of training sample quality emphasized. Finally, based on the proposed identification method and adaptive non-parametric model updating rule, online non-parametric identification for the three-degree-of-freedom ship motion states, course, and motion trajectory was conducted, and the identification results of the proposed scheme were compared with those of the non-adaptive identification scheme. The experimental results show that the proposed scheme can adaptively update the non-parametric model according to maneuvering condition changes. The mean absolute error (MAE) and root mean square error (RMSE) of the proposed scheme’s identification results are less than 0.11 and 0.18 respectively, while the MAE and RMSE of the non-adaptive method’s identification results are below 1.43 and 2.10 respectively. This fully validates the proposed scheme’s significant advantages in generalization, demonstrating higher identification accuracy and further confirming its applicability in various maneuvering conditions.

Dual-Objective Optimization Research of Ship Cargo Flow Allocation in Multi-Mode Dam-Crossing System
LIU Shun, ZHANG Yu, TIAN Hongwei, ZHENG Qianqian, TANG Kexin
, Available online  , doi: 10.3969/j.issn.0258-2724.20240399
Abstract:

The navigation congestion problem in water control projects is often caused by the unbalanced allocation of cargo flow for dam-passing ships. To address this issue, a cargo flow distribution optimization model was established to achieve dual objectives of minimizing the total additional cost for dam-crossing ships and waiting time. Furthermore, the service charge mechanism oriented to the dam-crossing mode of the new waterway was further introduced, and its influence on the transfer law of ship cargo flow was analyzed. Then, a non-dominated sorting genetic algorithm (NSGA)-II was applied to solve the model, and the dam-crossing ship cargo flow allocation optimization scheme was proposed. Finally, by taking the Three Gorges water control project as an example, the effectiveness of the proposed model and algorithm was verified. The results demonstrate that the algorithm successfully obtains multiple sets of Pareto frontier solutions and provides diversified ship cargo flow allocation schemes. The total additional cost of dam-crossing ships is significantly negatively correlated with their waiting time. The cargo flow allocation schemes of the dam-crossing system with dual and multiple modes can achieve the objectives of effectively reducing costs and improving efficiency. The optimized scheme with a multi-mode system reduces the total additional cost for dam-crossing ships and waiting time by 67.1% and 0.5%, respectively. Compared with the initial scheduling scheme, the proposed ship cargo flow allocation scheme results in a mode shift of at least 33.4% of the ship cargo flow, which alleviates navigation congestion in water control projects. Moreover, a reasonable service charge rate for the new waterway can effectively reduce the cost and improve the efficiency of dam-crossing ships.

Improvement in Numerical Solution of Probability Density Evolution Equation and Slope Stability Analysis
LIU Lang, ZHAO Yifei, SONG Yumin
, Available online  , doi: 10.3969/j.issn.0258-2724.20250470
Abstract:

To improve the computational accuracy and efficiency of the probability density evolution method in analyzing high-dimensional stochastic responses of complex structures, a total variation diminishing (TVD) scheme using Van Leer’s flux limiter and a hypersphere-based uniform point selection method are proposed. First, for the TVD scheme in difference solving, the limitations of the Roe-Sweby flux limiter are analyzed. Van Leer’s flux limiter is selected to optimize numerical solution accuracy, and its effectiveness is verified through numerical examples. Second, based on existing mathematical findings, the probability calculation formulas for the n-dimensional standard normal distribution within hyperspheres and hypercubes are derived, by which the hypersphere-based high-dimensional point selection method is quantitatively validated. Furthermore, a hypersphere-based good lattice point method is developed to generate uniform point sets and perform probability calculations within hyperspheres. Finally, taking a high slope project under construction as a case study, a probability density evolution equation for slope stability is established, and the proposed improved method is applied to calculate the slope stability coefficient. The results show that under normal service conditions and rainfall conditions, compared with 10 000 Monte Carlo simulations, the minimum errors of the probability integral for the slope stability coefficient are 0.4% and 0.2%, and the maximum errors are 2.9% and 4.2%, respectively. When addressing high-dimensional stochastic structural problems, the proposed numerical solution method demonstrates satisfactory computational accuracy and efficiency.

Management and Monitoring Suggestions for Beam-End Integration Device in Long-Span High-Speed Railway Bridge Based on Data Analysis and Physical Modeling
HE Qing, ZHOU Siyuan, WAN Qiushi, GUO Hui, LING Liepeng, WANG Ping
, Available online  , doi: 10.3969/j.issn.0258-2724.20240272
Abstract:

Beam-end integration devices in long-span high-speed railway bridges have complex structures and are prone to generating peak structural track irregularities, which directly affect train running stability. However, targeted monitoring and management criteria for such devices are currently lacking. In this study, monitoring management values for beam-end integration devices are investigated. Monitoring data collected in recent years are analyzed, and a finite element model of the beam-end integration device is established. Based on static irregularity standards, simulation analyses are conducted to propose monitoring management values for beam-end longitudinal displacement, lateral displacement, deflection angle, and fixed steel sleeper lifting and unsupported sleeper. For an operational speed of 250 km/h, the proposed Class Ⅰ monitoring management thresholds are as follows: beam-end longitudinal compression of 160 mm, lateral displacement of 7 mm, deflection angle of 3×10−3 rad, fixed steel sleeper lifting of 3 mm, and unsupported sleeper of 2 mm. The results show that the proposed monitoring management values can effectively identify abnormal conditions of the device. Longitudinal compression at the beam end significantly affects track profile and gauge irregularities; lateral displacement primarily influences gauge irregularity; both the fixed steel sleeper unsupported section and lifting impact track profile irregularity. The proposed monitoring management values provide a reference for the operation and maintenance of long-span high-speed railway bridges.

High-Fidelity Reconstruction Method for Rolling Bearing Fault Features Based on Bessel Series Empirical Wavelet Transform and Attention-Guided Denoising Autoencoder
JIN Junjie, MA Zixuan, ZHANG Qi, LI Bo, YANG Wenhua, LI Ye, LAI Engang, SUN Feng
, Available online  , doi: 10.3969/j.issn.0258-2724.20250592
Abstract:

To address the waveform distortion of rolling bearing fault features under strong background noise interference, a method combining Bessel series empirical wavelet transform (BSEWT) and attention-guided denoising autoencoder (AGDA) was proposed. First, the time-domain signal was mapped to a data-driven energy spectrum via Bessel basis expansion. By using a multi-scale peak detection strategy, frequency band boundaries were adaptively identified to separate modal components containing fault information. Second, an AGDA model was constructed to reconstruct the background noise distribution by using an attention weighting mechanism, reversely stripping noise components from noisy modes. Finally, the feature extraction property under five rotational speed conditions and various signal-to-noise ratio (SNR) levels was evaluated by using simulated signals and the public bearing dataset of Case Western Reserve University (CWRU). The results show that BSEWT adaptively narrows the frequency bands and finely isolates complex narrow-band interference better than the traditional empirical wavelet transform (EWT) in frequency band division. In simulation experiments with five different rotational speeds, when the input SNR exceeds 10 dB, the output SNR, root mean square error, and normalized correlation coefficient of the BSEWT-AGDA method are comprehensively superior to those of variational mode decomposition (VMD), empirical mode decomposition (EMD), and EWT. For input SNRs below −5 dB, the BSEWT method achieves the optimal output SNR index. In CWRU data validation, the signal processed by the proposed method yields a kurtosis of 8.163 5 and an impulse factor of 8.074 1, surpassing comparative algorithms like VMD and EMD. BSEWT overcomes the defect of relying on fixed frequency band division in traditional methods, achieving precise complex frequency spectrum segmentation. Combined with the AGDA model, it effectively eliminates noise interference in higher SNR ranges and achieves high-fidelity reconstruction of rolling bearing fault features.

Analysis of Impact Effect of Cable Breakage in Half-Through Railway Arch Bridges with CFRP Cables
ZENG Yongping, LIU Liwei, TAO Qi, WAN Xing, ZHANG Xun, JIA Hongyu
, Available online  , doi: 10.3969/j.issn.0258-2724.20240555
Abstract:

To study the effect of cable breakage on the impact response of concrete-filled steel tube arch bridges and the difference in safety factor requirements between carbon fiber reinforced polymer (CFRP) cables and steel cables, the dynamic response of a railway bridge under accidental cable breakage was analyzed. A spatial finite element model was established by ANSYS. The force characteristic variations of the residual structure of the arch bridge under five cable breakage conditions were studied based on the equivalent unloading method. The impact sensitivity of the structure after cable breakage was evaluated by dynamic amplification factor (DDAF) and demand capacity ratio (DDCR). The effects of different cable materials, namely steel cables and carbon cables, on the dynamic response of the arch bridge were compared. The results show that the dynamic response of the main girder and the stress of the arch rib are greatly affected by the position and number of cable breakages. The redistribution ratio of the cable force is inversely proportional to the distance from the broken cable area and the cable length and directly proportional to the number of failed cables. The DDAF of the arch bridge with carbon cables is higher than that of the arch bridge with steel cables, ranging from 1.190 to 1.430. The DDCR of the remaining cable after cable breakage does not exceed 1, indicating large redundancy. Compared with bridges with steel cables, arch bridges with carbon cables require smaller safety factors under cable breakage conditions, ranging from 1.000 to 1.500.

Numerical Simulation of Geyser Process Caused by High-Pressure Entrapped Air Release in Baffle-Drop Shafts
YANG Qian, YANG Qinghua, CHEN Feng
, Available online  , doi: 10.3969/j.issn.0258-2724.20240426
Abstract:

To reduce the structural and ground safety risks caused by the high-speed air–water mixture during geyser events in baffle-drop shafts, a numerical simulation was conducted to investigate the influence of void fraction and connecting pipe diameter on shaft pressure and geyser intensity. The variation patterns of impact loads on the baffles at the shaft bottom were analyzed, proposing the installation of a throttling orifice at the shaft midpoint to control the geyser intensity. The results show that the pressure in the connecting pipe first reduces and then increases with rising void fraction, reaching a minimum within the range of 0.2–0.4. Among the three diameter ratios considered, the geyser intensity reaches its maximum when the diameter of the connecting pipe is half that of the shaft. The impact load on the baffles decreases continuously from bottom to top, and for any given baffle, the impact load near the partition wall and shaft wall is greater than that at the baffle’s edge. Installing a throttling orifice at the shaft midpoint can effectively control the geyser intensity, and the impact load on the throttling orifice is 10 times greater than that on the bottom baffle. These findings provide a reference for the safe operation of urban deep tunnel drainage systems.

Effect of Eccentric Wear of Brake Shoe on Fatigue Crack Propagation Behavior of Wheel Treads
LI Zhiqiang, SHEN Xuwen, HAN Chen, CHENG Shipeng, YANG Zhiyong
, Available online  , doi: 10.3969/j.issn.0258-2724.20250427
Abstract:

Transverse eccentric wear of the brake shoe is a co mmon and critical braking fault in freight trains, which can lead to thermal damage in wheels. The mechanism of transverse eccentric wear of the brake shoe influencing thermal fatigue crack propagation in freight train wheels was investigated to support wheel condition-based maintenance and safety early warning. First, material tests were carried out to obtain the mechanical property parameters and fatigue crack propagation parameters of ZL-C wheels at various temperatures. Second, a thermo-mechanically coupled finite element model of the wheel was developed to simulate the temperature and stress fields under emergency braking conditions with different degrees of eccentric wear of the brake shoe. Finally, based on fracture mechanics theory, a predictive model for thermal fatigue crack propagation service life was established to analyze the crack propagation path and its service life. The results show that the eccentric wear of the brake shoe will lead to a significant braking temperature rise during the emergency braking of the wheel, and a large cyclic circumferential stress will be generated. The combined effect of the two will lead to the initiation of “scale-like” cracks at the chamfer of the wheel rim. The thermal fatigue crack mainly propagates along the axial direction of the wheel and propagates slowly along the radial direction. With the increase of crack size, the depth-to-length ratio of thermal fatigue crack decreases continuously. Under the influence of the eccentric wear of the brake shoe, the propagation service life of the thermal crack is greatly reduced. Under the emergency braking condition of the eccentric wear of the brake shoe of 25 mm, the unstable fracture occurs when the crack extends to 50.0 mm. Compared with that in the condition without the eccentric wear of the brake shoe, the thermal fatigue crack is extended from 0.2 mm to 50.0 mm when the eccentric wear of the brake shoe is 25 mm. The propagation service life of the crack is reduced by 90.1%. For the wheels with thermal fatigue cracks, it is necessary to further strengthen the monitoring, avoiding the early failure of the wheels caused by the unstable propagation of the cracks.

Machine Learning Model for Vibration Prediction of Deep-Sea Hydrate Mining Risers Based on CLR-ZOA-RF
GUO Xiaoqiang, TONG Yilin, LIU Yuechao, LI Yingwei, YANG Kelun, LV Junlin, LIU Jun
, Available online  , doi: 10.3969/j.issn.0258-2724.20250468
Abstract:

To address the challenges of vibration prediction for deep-sea natural gas hydrate mining risers under complex environmental loads, a machine learning model named CLR-ZOA-RF was proposed. This model enhanced the performance of the zebra optimization algorithm (ZOA) by incorporating chaotic mapping, Levy flight, and random perturbation mutation (CLR) strategies, which are then integrated with the random forest (RF) algorithm to achieve real-time vibration response prediction of mining risers under multi-factor influences. Based on the similarity principle, a vibration simulation experimental device for the mining riser was developed to conduct multi-factor riser vibration experiments and construct an experimental dataset. The research results indicate that the proposed CLR-ZOA-RF model exhibits superior generalization performance in processing multi-dimensional non-linear data. Specifically, the determination coefficients (R2) for in-line and cross-flow vibrations are improved from 0.935 9 and 0.926 4 of the ZOA-RF model to 0.971 9 and 0.966 7, respectively. Compared with those of the ZOA-RF model, the mean absolute error and root mean square error of the CLR-ZOA-RF model decrease by 26.83% and 33.72% for in-line vibration and by 45.15% and 29.49% for cross-flow vibration. It can be seen that the model achieves high prediction accuracy and good generalization, confirming the effectiveness of the hybrid intelligent algorithm.

Influence of Coiling and Annealing Temperatures on Hole Expansion Property of CP780 Hot-Rolled Ultra-High Strength Automotive Steel Sheet
JIN Yang, WENG Haibo, ZHENG Haoqing, QIAN Xingyu, RAO Tian’en, FAN Hua, YUAN Kaitong, LI Wei
, Available online  , doi: 10.3969/j.issn.0258-2724.20250576
Abstract:

Achieving a balanced combination of high strength and high hole expansion ratio is a critical research objective for the application of high-strength steel in automotive manufacturing. A novel low-carbon Ti-Nb microalloyed hot-rolled multiphase steel with a high hole expansion of 780 MPa grade was designed. The thermomechanical controlled process (TMCP) was employed to investigate the influence of different coiling and annealing temperatures on its microstructure and mechanical property evolution, as well as to elucidate the mechanisms governing the high hole expansion property. The results indicate that increasing the annealing temperature enhances the precipitation strengthening effect of Ti/NbC particles, promotes grain refinement of the experimental steel, and leads to a progressive increase in both yield and tensile strengths. However, elongation after fracture and hole expansion ratio exhibit a gradual decline with rising annealing temperature. At a coiling temperature of 480 ℃, the hole expansion property is primarily governed by the bainite/ferrite ratio. Elevated annealing temperatures result in increased bainite content, which restricts crack deflection and bifurcation, thereby causing a monotonic reduction in hole expansion ratio. At a coiling temperature of 520 ℃ and an annealing temperature of 610 ℃, an optimal matching among the fine-grained microstructure of experimental steel, an appropriate bainite/ferrite phase ratio, and dispersed fine Ti/NbC precipitates was realized. This microstructural configuration effectively promotes crack deflection, bifurcation, and microvoid formation during hole expansion, enabling greater energy consumption prior to fracture and resulting in a hole expansion ratio exceeding 60%.

Performance Study of Novel PV/T System Integrated with Micro-Channel Heat Pipe Module
CUI Yuanlong, ZHU Kefeng
, Available online  , doi: 10.3969/j.issn.0258-2724.20250491
Abstract:

Aiming at the problem that the photoelectric and photothermal conversion performance of photovoltaic/thermal (PV/T) systems is still constrained by low winter temperatures in cold regions of northern China, this paper proposes a novel solar photovoltaic/thermal integrated system integrated with a micro-channel heat pipe array and evaluates its energy conversion performance under actual operating conditions. The system consists of a PV module, an aluminum micro-channel heat pipe array using R141b as the working fluid, and a U-shaped manifold heat exchanger. Indoor experiments determine that the optimal mass flow rate of cooling water in the U-shaped manifold heat exchanger is 0.025 kg/s, achieving a balance between energy conversion efficiency and circulating pump energy consumption. Based on multi-day outdoor experiments, the temperature distribution and energy conversion performance of the system in winter and summer are investigated, as well as its net energy conversion performance when circulating pump energy consumption is considered. The experimental results of two days with representative weather conditions show that under winter outdoor conditions, the thermal efficiency, net electrical efficiency, and net overall efficiency of the system are 27.21%, 10.41%, and 37.62%, respectively; the average daily net overall energy output is 376.53 W, including a net electrical output of 127.54 W and a thermal output of 248.98 W. Under summer operating conditions, the corresponding thermal efficiency, net electrical efficiency, and net overall efficiency can reach 47.42%, 8.55%, and 55.97%, respectively; the average daily net overall energy output is 643.71 W, including a net electrical output of 110.49 W and a thermal output of 533.22 W. The experimental results show obvious seasonal performance differences.

Vibration Damping Performance of Grid-Type Periodic Acoustic Black Holes and Reinforcement Structures
WEN Xianglong, DENG Sixian, SONG Chunsheng, JIANG Youliang
, Available online  , doi: 10.3969/j.issn.0258-2724.20250577
Abstract:

A periodic acoustic black hole is a new type of vibration-damping structure that exhibits both bandgap characteristics and broadband vibration-damping performance. The current research on periodic acoustic black holes has difficulty in forming an omnidirectional low-frequency bandgap at the two-dimensional level and low structural strength. Based on this, the analysis method of combining real and complex energy bands was adopted to achieve low-frequency broadband vibration damping, and a grid-type periodic acoustic black hole structure was proposed, verifying its favorable omnidirectional bandgap characteristics at the two-dimensional level with low-frequency bandgap characteristics and broadband vibration-damping performance. Then, two reinforcement structures of grid-type acoustic periodic black holes were proposed, and the theoretical analysis was validated by employing the vibration transmission experiment. The results show that the first bandgap onset frequency of the reinforcement structure with stiffeners increases by 170.1%, while its bandgap range is reduced to 42.1% of the original value. For the reinforcement structure with stiffening grids, the first bandgap onset frequency increases by 69.3%, with its bandgap range reducing to 42.3% of the original one. The low-frequency bandgap characteristics of the reinforcement structures are not as good as those of the original structure, but the vibration-damping performance is better than that of the original structure in the medium and high frequency bands, and the reinforcement structures have excellent broadband vibration-damping performance. Through comparison of two reinforcement structures, the one with stiffening grids has better low-frequency bandgap characteristics and broadband vibration-damping performance.

Derailment Risk Influence of Rear-End Collision of Short Electric Multiple Units with Tight-Lock Couplers
WANG Mingmeng, ZHOU Jianxin, YANG Chao, ZHANG Luze, ZHU Tao, ZHANG Weihua
, Available online  , doi: 10.3969/j.issn.0258-2724.20250572
Abstract:

In order to investigate the influence law of vehicle parameters on derailment risk of rear-end collisions of electric multiple units, a mathematical model of a speed-sensitive tight-lock coupler was innovatively proposed and integrated into the vehicle dynamics model. The train collision dynamics model was established through the coupling of adjacent vehicles. Rear-end collisions of two identical short electric multiple units were simulated on straight tracks to obtain the dynamic responses of the vehicles under different parameters. Finally, the derailment risk of rear-end collisions was evaluated by using the derailment coefficient, wheel load reduction rate, wheel–rail contact force, wheel lift, and other derailment indices. The influence law of collision speed, leading car body mass, and anti-climber steady-state force on the derailment indices was investigated. The results show that under various collision conditions, the maximum values of wheel lift and wheel load reduction rate are observed at the right wheel of the third wheelset on the leading car of the active train during rear-end collisions of short electric multiple units on straight tracks. The maximum values of derailment coefficient and wheel–rail contact force are observed at the right wheel of the second wheelset on the leading car under the same conditions. Derailment indices increase with higher collision speed under identical initial conditions. The wheel lift increases from 0.01 mm to 33.40 mm when the collision speed increases from 25 km/h to 55 km/h. The maximum vertical wheel force increases from 83.8 kN to 159.0 kN, an increase of 89.7%, under the same conditions. The peak value of the wheel load reduction rate firstly decreases and then increases with the increase in the leading car body mass. The wheel lift increases by 83.6% compared with its maximum when the anti-climber steady-state force is minimized.

Study on Magnetic Force Characteristics of Ground Coils and Side-Wall Beams in High-Speed Superconducting Maglev Transportation
HE Zhonglin, ZHAO Chunfa, YU Ruodan, LI Kai, LI Hongmei, MU Hanlin, HE Feng
, Available online  , doi: 10.3969/j.issn.0258-2724.20260126
Abstract:

Clarifying the magnetic force characteristics of null-flux coils and side-wall beams under mechanical, electrical, and magnetic multi-field coupling has significant guiding value for the design of coil modules and the construction of track structures in superconducting electrodynamic suspension systems. To this end, an electrodynamic force numerical calculation model was established based on the dynamic circuit principle, referencing the parameters of the Yamanashi line electrodynamic suspension system in Japan. The electrodynamic force responses of a single null-flux coil, coil module, and side-wall beam under the full-load operation of the maglev train were calculated, and the basic characteristics of the system magnetic force were obtained. Subsequently, the influence of motion parameters of superconducting magnets on the electrodynamic force was explored, and the variation rules of dynamic loads of magnets under different motion speeds, sinking amounts, and lateral displacements were statistically analyzed. The results indicate that the load on the track follows an increase, stabilization, and decrease trend, and it takes at least 11 consecutive null-flux coils to stabilize; the main characteristic wavelengths in the frequency-domain response of electrodynamic force are the superconducting coil pole pitch of 1.35 m and the 2.5-times wavelength of 3.375 m; the magnetic force amplitude of a single null-flux coil decreases with the increase of speed, and the average magnetic force acting on the side-wall beam gradually increases. When the magnetic force involves multiple null-flux coils, its magnitude is jointly determined by the induced current and speed, and the fundamental frequency gradually shifts from the frequency corresponding to the 1.35 m wavelength to that of the 3.375 m wavelength; at the speed of 200–600 km/h, the fundamental frequency of the side-wall beam’s magnetic force response ranges from 17 to 49 Hz. By considering that the commonly used operating speed of maglev trains is 400–600 km/h, it is recommended that the fundamental frequency of the side-wall beam structure should not exceed 33 Hz. The lateral force acting on the track is sensitive to the lateral displacement of superconducting magnets; to enhance the dynamic stability of the side-wall beam structure, significant lateral displacements of magnets should be avoided as much as possible.

Research on Shear Property of Special Shear Keys in Composite Structures of Double-Layer Plate Truss
MA Chi, REN Canqing, CHEN Tiandi, HE Jiayuan, PENG Zhuo
, Available online  , doi: 10.3969/j.issn.0258-2724.20250558
Abstract:

To study the shear mechanical properties of special shear keys in composite structures of a double-layer plate truss, the influence of concrete strength and the hole diameter of transverse perforated stiffening ribs on the shear property of special shear keys was explored through push-out tests. The refined numerical simulation analysis was verified, revealing the failure mode of specimens and the collaborative working mechanism of special shear keys. Based on numerical simulation for parameter analysis, a calculation formula for shear bearing capacity suitable for special shear keys was proposed. The results show that the composite structure of a double-layer plate truss bears the shear force by the transverse perforated stiffening ribs and studs. As the concrete strength increases, the shear bearing capacity of the shear keys gradually improves. The concrete strength increases from C30 to C50, and the shear bearing capacity of the shear key rises from 341.16 kN to 448.09 kN, an increase of 31.34%. The ductility of the shear key with a hole diameter of 25 × 20 mm is 2.06 mm, which is 89.0% higher than that without a hole. Moreover, the greater strength of the concrete leads to a greater influence on the ductility of the shear key. The stress at the root of the stud first reaches the yield strength, and then the transverse perforated stiffening rib reaches the yield strength. Parameters, such as the diameter of the stud, the ultimate tensile strength of the stud, the strength of the concrete, and the thickness of the transverse perforated stiffening rib, affect the shear bearing capacity of the special shear key. Based on the parameter analysis, a calculation formula suitable for the shear bearing capacity of the special shear key is proposed. The prediction results are in good agreement with the finite element results, providing a reference basis for engineering design.

Effect of Loading Angle on Crack Propagation in 6005A-T6 Aluminum Alloy Under Different Degradation Levels
YANG Bing, ZHANG Zhe, LIU Jinbang, XIAO Shoune, YANG Guangwu, ZHU Tao, CHEN Dongdong
, Available online  , doi: 10.3969/j.issn.0258-2724.20250429
Abstract:

To investigate the effects of loading angle on the crack propagation life and damage mechanisms of 6005A-T6 aluminum alloy under different fatigue degradation levels, a systematic fatigue crack propagation test was designed and conducted. The compact tension shear specimens were employed, and the comparative tests were performed under four degradation levels and four loading angles. The crack propagation rate was calculated and evaluated by using the seven-point incremental polynomial method. First, crack propagation curves were obtained under non-degraded and different degrees of pre-fatigue conditions. Second, the variation law of residual life at different loading angles was compared. Finally, the local strain distribution at the crack tip was analyzed by using digital image correlation (DIC) technology, and fracture morphology was observed to reveal the crack propagation mechanisms. The results show that the loading angle significantly affects the crack propagation behavior. As the loading angle increases, the crack propagation rate decreases significantly, and the residual life is greatly extended. Under the non-degraded condition, the degree of life extension can reach up to 186%. Under the 1.5 × 105 cycles of degradation, the degree of life extension drops to approximately 140%, indicating that the effect of loading angle on life extension diminishes as degradation deepens. DIC results show that increasing the loading angle reduces the local strain at the crack tip, decreases the sensitivity of the crack to the main loading direction, and exhibits a stronger propagation hysteresis effect. In terms of fracture morphology, as the loading angle increases, the fracture roughness increases, and the cleavage steps and tear ridges gradually become coarser, with the crack propagation path becoming more complex, revealing the dominant role of loading angle in the crack failure mode.

Effect of Phase Change Materials on Thermal Properties of Ceramsite Foam Concrete
LI Fuhai, YANG Bo, HUO Jiateng, LI Jin, YANG Zongchi, CHEN Zhao
, Available online  , doi: 10.3969/j.issn.0258-2724.20250187
Abstract:

To investigate the effect of phase change materials (PCMs) on the thermal properties of ceramsite foam concrete, phase change energy storage ceramsite lightweight foam concrete, namely PCESLFC was prepared by using solid–solid phase change material (SS) and phase change microencapsule (MC). The effects of different PCM dosages on the thermal conductivity, specific heat capacity, phase change temperature, and phase change enthalpy of the ceramsite foam concrete were studied. A concrete panel heat transfer test was conducted to evaluate the thermal insulation property of PCESLFC, and a miniature test house was constructed to assess its thermal regulation property under natural environments. The results show that the incorporation of MC effectively reduces the thermal conductivity and increases the specific heat capacity of the ceramsite foam concrete. The MC-100% group exhibits a reduction of 18.16% in thermal conductivity and an increase of 21.1% in specific heat capacity compared to the control group. In contrast, the SS has limited enhancement effects on the thermal properties of the ceramsite foam concrete. Differential scanning calorimetry (DSC) analysis reveals that the MC has a higher phase change enthalpy, with the MC-100% group achieving a phase change enthalpy of 9.08 J/g, which is 7.32% higher than that of the SS-100% group. Temperature damping coefficient calculations indicate that the thermal insulation property improvement of the ceramsite foam concrete in the SS group is significantly lower than that in the MC group. Temperature monitoring in the miniature test house shows that the occurrence time of peak temperature in the PCESLFC test house is delayed by 20 min, and the peak temperature reduces by 0.8 ℃ compared to that of the control group. The occurrence time of minimum temperature is delayed by 20 min, and the valley temperature increases by 0.4 ℃. MC shows excellent thermal regulation properties to ceramsite foam concrete, making it suitable for building envelope systems.

Design and Optimization of Coupling Mechanism for Wireless Energy Transmission of Hovering Unmanned Aerial Vehicles
DONG Liang, SU Yuming, FAN Zhixun, SU Zhouhao
, Available online  , doi: 10.3969/j.issn.0258-2724.20250261
Abstract:

To address the problems of limited endurance of unmanned aerial vehicles and unstable power transmission caused by positional offset during hovering charging, a T-shaped asymmetric magnetic coupling mechanism based on demagnetizing field theory was proposed. Firstly, the magnetic field concentration ability was enhanced, and the demagnetization effect was suppressed by optimizing the magnetic circuit design, thereby improving the coupling coefficient; the influences of key parameters such as the upper side length, lower side length, and middle magnetic column height of the launch mechanism on the coupling performance were analyzed. Secondly, the optimal structural dimension combination was determined using a parametric sweep method. Finally, an experimental prototype was built to conduct performance testing. The research results indicate that the mutual inductance fluctuation rate of the structure is only 2.3% within a rotational offset range of ±20° and 4.1% within an axial offset range of ±20 mm, exhibiting excellent anti-offset characteristics; the mutual inductance fluctuation rate under rotational offset is 2.4%, and the output power fluctuation rate within an axial offset range of ±20 mm is 4.2%, which are basically consistent with the simulation results; under varying load conditions, the system can achieve constant current output, and the designed coupling mechanism has good robustness and stability under dynamic working conditions. The proposed T-shaped coupling mechanism has advantages such as compact structure, strong anti-offset ability, and stable output and is suitable for unmanned aerial vehicle application scenarios requiring continuous power supply, such as relay communication and power inspection.

Vehicle Detection and Tracking Algorithm for Distant Small Targets Based on Improved YOLOv5s-CSC and Deep-Sort
HU Xin, SHAO Liangbin, ZHOU Yun-qiang, XIAO Jian, CHENG Hongliang
, Available online  , doi: 10.3969/j.issn.0258-2724.20240528
Abstract:

To address the problems of missed detection of distant small vehicle targets, as well as ID switching and jumping during tracking in traffic surveillance scenarios, a small-target detection and tracking algorithm combining improved You Only Look Once version 5 small (YOLOv5s) and improved deep simple online and realtime tracking (Deep-Sort) was proposed. For detection, based on the YOLOv5s baseline network, the convolutional block attention module (CBAM) mechanism was introduced to enhance the feature capture of small targets. The Swin Transformer Block was fused with the cross stage partial (CSP) structure in the backbone network to construct the CSPSTB module and strengthen global information attention. The content-aware reassembly of features (CARAFE) upsampling operator was adopted to reduce feature loss, forming the YOLOv5s-CSC detection network. For tracking, a deep appearance feature extraction network containing a residual squeeze-and-excitation (Res_SE) structure was designed, which was more adaptable to vehicle features during data association. Experimental results on the UA-DETRAC dataset indicate that, compared with the original YOLOv5s model, the improved detection algorithm increases the mean average precision (mAP), mean average precision at an intersection over union (IoU) threshold of 0.50 (mAP50), and mean average precision at an IoU threshold of 0.75 (mAP75) by 1.5%, 1.5%, and 1.9%, respectively, and the small-target detection accuracy improves by 4.2%. Tracking experiments on the VERI-Wild dataset show that, after combining the improved detector and tracking network, the multi-target tracking accuracy, high-order tracking accuracy, and IDF1 (identity F1) score increase by 9.6%, 3.2%, and 7.7%, respectively. The combination of YOLOv5s-CSC and the improved Deep-Sort effectively overcomes difficulties of tracking failure of distant small vehicle targets, as well as ID switching and jumping during tracking, verifying the superiority of the improved model in tracking accuracy and robustness.

From Status to Theory: Connotation of Green Design in Railway Engineering
LI Yajuan, BAO Xueying
, Available online  , doi: 10.3969/j.issn.0258-2724.20250181
Abstract:

With the national initiative for green transport development, it has become a key research focus to promote green railway engineering through green design and achieve sustainable railway development at the source. To reveal the connotation of green design in railway engineering, first, the relevant definitions of green design were systematically reviewed. A keyword cloud atlas was generated using Python to identify key elements of green design definitions. Based on practical features of green railway engineering, a definition was thus formulated for green design in railway engineering. Then, based on bibliometrics and knowledge mapping, CiteSpace was applied to visualize characteristics of keyword structures and temporal clusters in the research content of green design in railway engineering. On this basis, the fundamental connotation was summarized. Through a “railway engineering-resources & environment” bidirectional coupling, the connotation was reinterpreted from the perspective of compatible symbiosis. With a design criterion that engineering impact intensity must be less than environmental capacity, a “railway engineering-resources & environment” compatible symbiosis framework as well as its characterization model were constructed to quantify the compatibility between engineering impact intensity and environmental capacity. Therefore, green design schemes were allowed to shift from qualitative comparison to quantitative decision-making. Finally, a logical analytical framework was established for green design connotation from four aspects: green design elements, their coupling mechanisms, performance assessment, and optimization/regulation mechanisms, so as to systematically deconstruct the connotation. The results are expected to provide theoretical guidance for construction of green railway engineering.

Research on Seismic Performance of Precast Concrete Beam-Column Joints with Bolt-Plate Mechanical Connection
ZHANG Liping, WU Yi, YANG Chun, CAO Zhongming, LI Tianxian, CHEN Qingjun
, Available online  , doi: 10.3969/j.issn.0258-2724.20250104
Abstract:

To simplify the connection construction of precast concrete beam-column joints and improve construction efficiency, a precast beam-column joint with bolt-plate mechanical connection was proposed. The pseudo-static loading tests of seven specimens were carried out to analyze the effects of high-strength bolts, longitudinal reinforcement ratio of the beam, anchorage length and methods, as well as embedded steel plates on the enlarged head surface on the seismic performance of the joints. Results show that the joints exhibit excellent bearing capacity and ductility, with a reliable mechanical connection and a ductility coefficient of 3.79–5.69. The targets of plastic hinge failure at the beam end and the ductile design principle of “strong joint” can be achieved. When the anchorage length of longitudinal reinforcement of the beam is increased from 200 mm to 300 mm, the joint’s peak bearing capacity, ductility coefficient, and cumulative energy dissipation can be enhanced by 24.38%, 16.36%, and 38.08%, respectively. The peak bearing capacity and ductility coefficient of the joint with embedded steel plates on the enlarged end surface can be increased by 18.69% and 13.29%, respectively. While the reinforcement ratio of longitudinal reinforcement of the beam is increased from 0.75% to 1.18%, the failure mode of the specimen transfers from flexural-shear failure at the beam end to wedge-shaped failure at the enlarged head. The peak bearing capacity and cumulative energy dissipation can be improved by 20.65% and 20.30%, but the ductility coefficient decreases by 16.34%. Optimizing bolt edge distance and spacing, as well as enhancing the parameters including reinforcement along the beam direction on the enlarged head are recommended for subsequent research, thereby ensuring the rational formation of a plastic hinge at the beam end and ultimately achieving the ductile design principle of “strong joint”.

Research on Modular Functionality Recovery Functions for Cable-Stayed Bridges
DUAN Jiahong, LUO Tian, YIN Weitao, ZHANG Yang, LIAO Luoheng, CHEN Zhiyu, ZHAO Canhui
, Available online  , doi: 10.3969/j.issn.0258-2724.20250204
Abstract:

To quantify the functionality recovery weighting coefficient in the seismic resilience assessment of cable-stayed bridges and avoid the subjectivity relying on experience or expert questionnaires, a seismic resilience assessment method based on a modular functionality recovery function was proposed. First, a Pair Copula hierarchical model was employed to divide a cable-stayed bridge into three levels, i.e., components, subsystems, and the system, and the seismic fragility of each level was calculated. Second, based on the probability ratio of each subsystem under different damage states, probability weighting was performed to construct a quantitative model of the functionality recovery weighting coefficient. Furthermore, a modular functionality recovery function was established by comprehensively considering functionality recovery, repair time, and repair path. Then, a system-level functionality recovery function was formed by integrating the repair sequence of subsystems. Finally, a seismic resilience assessment procedure for cable-stayed bridges was established, and a seismic resilience assessment was conducted taking a cable-stayed bridge as an example. The results indicate that the functionality recovery weighting coefficient varies with the ground motion intensity and is closely related to the damage degree of each subsystem. If a constant weighting coefficient is adopted, the recovery contribution of the bearing subsystem is underestimated by approximately 50%. The seismic responses of the components of cable-stayed bridges have strong correlations; if the components are assumed to be completely independent, the calculation error of the system fragility can reach 28%–53.3%. Within the ground motion intensity range specified by codes, if the repair sequence of subsystems is neglected, the system resilience index is underestimated by approximately 9%.

Experimental Study on Axial Compressive Mechanical Performance of Concrete-Filled Short Double-Skin Composite Tubular Columns with Inner Corrugated Steel Pipe
LIU Zhonghua, SHU Ganping, YAO Zhen, BIAN Zhaowei, YUN Zuoyi, CHENG Hao, MAO Dongxu
, Available online  , doi: 10.3969/j.issn.0258-2724.20250410
Abstract:

To address the problems of premature inward buckling of the inner steel pipe and insufficient confinement to concrete in traditional concrete-filled double-skin composite tubular columns, a novel concrete-filled double-skin composite tubular column with an inner corrugated steel pipe was proposed. By taking the hollow ratio, nominal steel ratio, wall thickness of the corrugated steel pipe, and inner tube configuration as the main parameters, a systematic experimental study on the axial compressive performance of eight concrete-filled short double-skin composite tubular column specimens with inner corrugated steel pipe was conducted. Based on the experimental results, a finite element model was established to further discuss the load-bearing mechanism, confinement effect, and influence rules of key structural parameters of the components under axial loads. The results indicate that all concrete-filled double-skin composite tubular columns with inner corrugated steel pipe exhibit good ductility under axial compression. The failure mode is dominated by local outward buckling of the outer steel pipe, while no obvious instability occurs in the corrugated steel pipe. With an increase in the hollow ratio, the bearing capacity of the components decreases significantly; the nominal steel ratio has a pronounced effect on the bearing capacity, whereas the wall thickness of the corrugated steel pipe has a minor effect on the ultimate bearing capacity. Compared with traditional concrete-filled double-skin composite tubular columns, the peak bearing capacity of the novel components decreases slightly, but the initial stiffness and ductility improve significantly. The strength indices of all specimens are greater than 1.00, reaching a maximum of 1.33, which indicates that a more effective composite effect is formed between the corrugated steel pipe and the concrete. The finite element calculation results further reveal that the corrugated steel pipe provides strong lateral confinement to the sandwiched concrete mainly through the trough regions, while the helical angle has a minor effect on the axial compressive performance, and the waveform parameters have a certain effect on the bearing capacity. The research results provide a reference for the engineering application and design of concrete-filled double-skin composite tubular columns with an inner corrugated steel pipe.

Electrical Coupling Effects and Dynamic Control Strategy for Urban Rail Transit Power Supply Systems and Buried Pipelines
LI Qiaoyue, ZHENG Xin, DU Guifu, ZHOU Tongtong, ZHANG Dongliang
, Available online  , doi: 10.3969/j.issn.0258-2724.20250656
Abstract:

To investigate the interference of stray currents from the power supply system of urban rail transit lines on adjacent buried metal pipelines and their control, the electrical coupling effects and dynamic control strategies between the urban rail transit power supply system and buried pipelines were analyzed. First, a dynamic coupling parameter calculation model integrating dynamic control strategies was proposed for the urban rail transit power supply system and buried metal pipelines. Based on actual line parameters, simulations were conducted to analyze the dynamic distribution of power supply system parameters and pipeline parameters under different coupling scenarios, clarifying the influence of multi-return-path coupling interactions. Second, to address the electrical coupling interference between the urban rail transit power supply system and pipelines, a dynamic control strategy was proposed, and the structure of a dynamic control device was designed, with its control characteristics analyzed. Finally, by using parameters from an actual urban rail transit line in China as an example, simulation analysis was performed to verify the dynamic control effect on pipeline stray current and pipe-to-soil potential. The results indicate that electrical coupling between the urban rail transit power supply system and adjacent pipelines leads to an increase in both the return system parameters of the urban rail transit and the pipeline safety parameters. By taking simulation case 3 as an example, the maximum instantaneous pipeline potential reaches 16.2 V. The coupling effect is related to the insulation condition of the running rails, the insulation condition of the pipeline, and the rail potential at the coupling location. The proposed dynamic control strategy can effectively suppress the electrical coupling impact of the urban rail transit power supply system on pipelines, clamp the pipeline voltage within the standard limits, and redirect the stray current on the pipeline back to the urban rail transit power supply system.

Optimization Method for Measurement Points of Grounding Grids Based on a Composite Weighted Vertex Cover Model
PEI Feng, JIA Lulu, ZHENG Shaohua, XU Ruicheng, CHENG Hongbo
, Available online  , doi: 10.3969/j.issn.0258-2724.2026067
Abstract:

To address the problem of disconnection between topology and electrical characteristics in the optimization of measurement points of grounding grids and to enhance the accuracy and engineering applicability of fault diagnosis for grounding grids, grounding grids were transformed into an undirected graph model with node-branch relationships. The topological connectivity characteristics of nodes were measured using node degrees, and the electrical connection characteristics of nodes were measured using node sensitivities. A determination model of node composite weight integrating node degree and node sensitivity was constructed. A greedy algorithm based on dynamic tuning of composite weights was adopted to optimize the node weight model, and an optimized layout scheme of measurement points for grounding grids that can cover all branches was selected. Simulation results show that this method can select 25 measurement points from 49 nodes in a 7 × 7 grounding grid, achieving complete electrical observability of all branches of the grounding grid. Under the same number of measurement points, the average error of single-branch corrosion diagnosis is reduced by 15.6% compared with the network topology hierarchical reduction method, and the average error of multi-branch corrosion diagnosis is reduced by 13.4% compared with the network topology hierarchical reduction method. Furthermore, even with fewer measurement points, the multi-branch diagnosis error is still superior to that of the agglomerative hierarchical clustering method, achieving accurate diagnosis of single and multiple branches with different corrosion degrees (mild, moderate, and severe). The proposed method can reduce the on-site testing workload while ensuring the accuracy of corrosion diagnosis, and it has important practical significance for the precise maintenance of grounding grid corrosion.

Analysis of Characteristic Parameters of Seismic Design Spectrum in Yunnan Region Based on Pseudo Acceleration Spectrum
LI Han, BAI Yu, YE Liaoyuan, PAN Wen
, Available online  , doi: 10.3969/j.issn.0258-2724.20250355
Abstract:

To investigate the issue that the response spectrum for structures with high damping exceeds that with low damping in the long-period range in code design spectra, strong earthquake records with magnitudes greater than 4.5 in the Yunnan region since 2007 were collected based on the pseudo acceleration response spectrum theory considering the true inertial forces of structures. A classification system based on the site category classification standard, magnitude, epicentral distance, and design earthquake group was established. A differential evolution algorithm with global and neighborhood mutation was employed to calibrate and obtain the average values of the calibration parameters under each category. On this basis, the effects of various factors on the characteristic period, the plateau value of the amplification factor spectrum, and the decay index were obtained, and the recommended parameter values for the Yunnan region were ultimately proposed. The results indicate that both site category and epicentral distance have a positive correlation with the characteristic period $ {T}_{\text{g}} $. The recommended values for design earthquake groups under different site categories proposed based on the calibration results are increased by 0.05 s on average compared with the code recommended values. Epicentral distance and magnitude have a weak effect on the plateau value of the amplification factor spectrum $ {\beta }_{\max } $, and no clear regularity is shown. The calibration result of $ {\beta }_{\max } $is larger than the code value. Considering the statistical results and engineering safety redundancy comprehensively, it is suggested that the value range of $ {\beta }_{\max } $ in the Yunnan region be set to 2.575–2.675. The damping ratio is the most dominant factor affecting the decay index γ. The linear relationship between γ and the damping ratio $ \xi $ quantifies the effect of the damping ratio on the spectral decay segment, providing a basis for reasonably considering the damping effect.

Dynamic Gate Assignment Method for Improving Flight Bridge Rate
MOU Qifeng, ZHANG Yi, WU Yuankai
, Available online  , doi: 10.3969/j.issn.0258-2724.20250069
Abstract:

Flight bridge rate is a key indicator for evaluating aviation service quality and airport operational efficiency. To improve the flight bridge rate, an optimized dynamic gate assignment method based on existing operational workflows and towing conditions was proposed. By arranging the full-day gate assignment plan in stages, and comprehensively considering constraints such as parking time limits for towed flights, towing eligibility conditions, and towing gate type restrictions, a gate assignment model taking the improvement of the flight bridge rate as the core objective was constructed. The number of flights assigned to contact gates was optimized by the model based on the towing penalty cost to maximize the overall efficiency. Meanwhile, to address the problem that the traditional pre-flight gate assignment is difficult to adapt to dynamic changes in flight operations, a multi-stage dynamic assignment algorithm based on an improved genetic algorithm was designed. In this algorithm, a continuous assignment mechanism was adopted to lock the completed assignment results, and the gate assignment for newly added flights was continuously optimized. Finally, simulation experiments were carried out using the actual operational data of a busy large-scale airport in Southwest China to verify the effectiveness of the proposed model and method. The results show that compared with the traditional gate assignment model without considering towing, the model considering the towing strategy increases the flight bridge rate by 1.1%; on this basis, compared with the traditional genetic algorithm, the adopted multi-stage dynamic assignment method increases the flight bridge rate by 4.0% and decreases the flight towing frequency by 16.6%. The proposed method can effectively balance the assignment of contact gates and aircraft towing costs, further optimize the flight bridge rate on the premise of controllable towing costs, and provide a feasible optimization method for the improvement of the flight bridge rate and dynamic gate assignment at airports.

Digital Twin-Based Life Prediction of Full-Ceramic Rolling Bearings
ZHANG Ke, HONG Yang, GAO Tianhao, BAI Xu, WANG Nan, SHI Huaitao, LONG Yanze
, Available online  , doi: 10.3969/j.issn.0258-2724.20250347
Abstract:

To address two challenges in the health monitoring of full-ceramic rolling bearings, namely the difficulty in accurately characterizing the crack propagation mechanism and the limited accuracy of life prediction caused by the scarcity of monitoring data samples under extreme operating conditions, a hybrid network model for high-precision life prediction under small-sample conditions was developed. A two-degree-of-freedom dynamic model was established based on bearing operating characteristics to analyze and describe the nonlinear dynamic behavior of the bearing through vibration responses. Time-domain statistical and frequency-domain energy features were combined to conduct a joint characterization of state information. By considering the characteristic that the crack propagation process of brittle materials is difficult to directly observe, physical constraints governing stable crack propagation were integrated into the physics-informed neural network training, ensuring simulated signal evolution aligned with actual crack propagation trends. A generative adversarial neural network mapped simulated signals to real signals, producing high-fidelity simulated vibration data to expand the training sample space and alleviate small-sample limitations. A random forest regression model was applied for life prediction and identification of critical crack instability states, forming a complete analysis process of residual life assessment for ceramic bearings. Results indicate that the method accurately captures the laws of dynamic evolution from crack initiation to propagation, achieving an average life prediction accuracy of over 99.5%. The simulated signals fused with physical constraints show high consistency with real signals in the energy distribution of key feature frequency bands and peak locations of the envelope spectrum. This method improves the capability of identifying crack propagation trends and enables the recognition of the critical node corresponding to the transition from stable to unstable crack propagation.

Study on Multi-Factor Coupling Analysis of Temperature Field and Freezing Radius Prediction for Cold-Region Tunnels
YAN Qixiang, HE Wencheng, YANG Yifan, ZHAO Zechang, YANG Xiao
, Available online  , doi: 10.3969/j.issn.0258-2724.20250126
Abstract:

To clarify the evolution laws of surrounding-rock temperature field and freezing radius in cold-region tunnels under the combined action of natural wind and train piston wind, and to improve the capability of rapid frost-damage risk assessment, a joint study of numerical simulation and machine learning was conducted using the Guolashan Tunnel as an engineering background. Based on field parameters, a three-dimensional transient heat transfer model consisting of air, primary support, secondary lining, and surrounding rock was established. Non-isothermal flow was used to couple the airflow with the heat transfer process of the lining and surrounding rock, and an equivalent wind speed method was adopted to characterize the effects of train piston wind and residual wind. Seven types of factors, including natural wind speed, natural wind temperature, natural wind direction, tunnel cross-section size, initial rock temperature, train speed, and train operation frequency, were analyzed through the control variable method. The longitudinal and radial temperature distributions and the freezing radius after 90 days were extracted, and random forest, support vector machine, and XGBoost prediction models were established based on 658 sets of orthogonal samples. The results indicate that the temperature-adjustment zone at the tunnel portal expands with the increase of cold-air action time; the longitudinal temperature of the surrounding rock presents two-stage characteristics of rapid rise and stable growth, and the radial temperature gradually approaches the initial rock temperature. Natural wind temperature, natural wind direction, and initial rock temperature are the dominant factors, with sensitivity coefficients of 0.22, 0.21, and 0.21, respectively. The XGBoost model has the optimal prediction accuracy with a root mean square error of 0.118; the predicted freezing radius for the engineering case is 1.85 m, with a relative error of 1.1%. External low-temperature air conditions and the initial thermal state of the surrounding rock are the key factors for frost-damage risk control in cold-region tunnels. The established models can provide a basis for rapid estimation of the freezing radius and optimization of anti-freezing design parameters.

Research on Mechanical Properties of Glacial Tills in Purang Region of Xizang
JIA Mincai, ZHENG Yiming, HUANG Jin
, Available online  , doi: 10.3969/j.issn.0258-2724.20240406
Abstract:

To reveal the mechanical properties of glacial tills in the Purang region of Xizang, in-situ direct shear tests with normal pressures from 100 kPa to 400 kPa were carried out on surface glacial tills in a natural state, and laboratory large-scale direct shear tests with normal pressures from 100 kPa to 400 kPa and large-scale triaxial tests with confining pressures from 100 kPa to 400 kPa were carried out on glacial tills with the compaction degree of 96%. The results show that the cohesion of surface glacial tills with a compaction degree of 91.7% is 11.0 kPa, and the internal angle of friction is 41.0°. The cohesion of glacial tills with a compaction degree of 96% is between 9.4 kPa and 11.2 kPa; the internal angle of friction is between 45.3° and 46.7°; the strength parameters obtained from laboratory large-scale triaxial tests are higher than those from large-scale direct shear tests. The peak strengths of glacial tills with a compaction degree of 96% are higher than those in a natural state, but the initial moduli are lower than those in a natural state. The stress–strain curve of glacial tills exhibits softening characteristics under various confining pressures, and the peak strain shows a trend of an increase followed by a decrease with the increase of confining pressures. The modified Duncan-Chang model can well describe the relationship between deviatoric stress and axial strain of glacial tills and reflect the strain softening characteristics of glacial tills in the Purang region.

Simulation Research of Transient Cavitation Flow Characteristics of Diesel in High-Pressure Nozzles Considering Thermal Effects
LIU Jianshi, ZUO Zhaowei, LIU Weilong, ZHAO Jianhui
, Available online  , doi: 10.3969/j.issn.0258-2724.20250379
Abstract:

The fuel flow inside the high-pressure nozzle exhibits significant pressure drops and high flow velocities, leading to a high tendency towards gas-liquid two-phase flow with intense transient heat and a coupling between cavitation flow and thermal effect. Turbulence and cavitation models were corrected to incorporate thermal effects, and a non-isothermal compressible flow model for fuel within the nozzle was established. Simulation research of the transient cavitation flow of diesel in a high-pressure nozzle was conducted. Results show that 1) during the opening stage of the needle valve, cavitation develops rapidly, with its cavitation intensity influenced by the inlet pressure. At 180 MPa, the average cavitation volume fraction inside the orifice increases by approximately 24% compared to that at 60 MPa. Cavitation induces significant fluctuations in the fuel mass flow rate, thereby influencing injection stability. During the closing stage of the needle valve, the orifice is almost blocked by cavitation; the effective flow area is reduced, and the mass flow rate is significantly decreased. 2) The fuel temperature within the nozzle exhibits a non-uniform distribution. Fuel decompression, expansion, and endothermic phase transition reduce the local temperature and form a supercooled region. Temperature rise regions result from viscous friction of fuel caused by velocity gradients between the core orifice region and the near-wall region, with the maximum temperature reaching 410 K. 3) Fuel cavitation and temperature variations interact. Cavitation reduces fuel flow velocity, increasing velocity gradients between fluid layers and amplifying internal friction and kinetic energy losses of fuel, which elevates fuel temperature. The temperature rise subsequently increases saturated vapor pressure, further promoting cavitation development.

Matrix Analysis of Variable Cross-Section Curved Beam Elements Based on Energy Method
WANG Lijuan, ZHAO Lei, JIANG Ning, QIAO Junting, LIU Shizhong, ZHANG Jiangong
, Available online  , doi: 10.3969/j.issn.0258-2724.20240256
Abstract:

To obtain the analytical solution formula of the stiffness matrix of a 2-node 6-degree-of-freedom explicit variable cross-section curved beam element, the variable cross-section curved beam element was incorporated into the widely used finite element system of rod structures. Based on Castigliano’s second theorem, the stiffness matrix of the spatial variable cross-section curved beam element was derived, and the flexibility matrix of the variable cross-section curved beam element was deduced by using the principle of stationary potential energy of elastic bodies; then, the cantilever end stiffness matrix of the variable cross-section curved beam element was obtained through inverse calculation. According to the static equilibrium condition and the principle of virtual work, the element stiffness matrix of the overall variable cross-section curved beam was obtained. In addition, the finite element formulation of the variable cross-section curved beam element could be degenerated into the formulation of uniform cross-section and uniform curvature beam elements and the standard formulation of uniform cross-section straight beam elements. A static calculation program for variable cross-section curved beam bridges was developed using MATLAB, and it was compared with the ANSYS solid finite element model to verify the static analysis theory of variable cross-section curved beams. The results indicate that the developed theory analyzes the bending-torsion coupling of the variable cross-section curved beam based on the finite element theory. The maximum deflection error between the variable cross-section curved beam theory and the ANSYS solid finite element model is 3.72%, and the maximum error compared with the ANSYS beam element model is less than 0.5%; for variable cross-section curved beam, the maximum deflection error between the variable cross-section straight beam degenerated from theoretion degeneration and the ANSYS solid finite element beam model is 1.72%, and the maximum error compared with the ANSYS beam element model is less than 0.1%.

Numerical Simulation of Propagation Characteristics for Three-Dimensional Landslide Surges and Its Force on Square Pier
KANG Azhen, WU Xianbin, ZHANG Dongming, PAN Siqi, LUO Bing
, Available online  , doi: 10.3969/j.issn.0258-2724.20250172
Abstract:

To study the propagation characteristics of three-dimensional landslide surges and the mechanism of the surging force on the square pier, based on the general moving objects model (GMO) module within Flow-3D, a three-dimensional numerical model was established to simulate the fluid-solid coupling interaction among the landslide mass, the water body, and the square piers. Firstly, by taking the site area of a continuous rigid frame bridge in a certain reservoir area as the research background, the propagation process of surges triggered by the entry of a potential landslide body into the water was simulated. The propagation characteristics of landslide surges were then explored based on the time domain and frequency domain analyses. Secondly, by taking the main piers and side piers of the bridge as the research objects, the characteristics of surging forces acting on square piers of different sizes were investigated numerically. Finally, a simplified algorithm for the surging force on square piers was proposed based on the Morison equation. The results were also compared with the estimation results obtained by the harmonic superposition method. The research shows that during the propagation of landslide surges, the head wave height decreases along the way; the secondary wave absorbs the high-frequency energy of the head wave, and it increases along the way. The surging force on the square pier increases nonlinearly with the increase of structural size and exhibits strong nonlinear characteristics where the absolute value of the valley is significantly greater than the peak value. Based on numerical simulation, the Morison equation can be applied to effectively calculate the surging force on a square pier and accurately capture its temporal characteristics. However, the evaluation method combining the harmonic superposition method and the Morison equation underestimates the surging force on the square pier.

Pore Pressure Inversion Model for Tunnel Lining Concrete Spalling Driven by Fire Test Temperature
XIONG Tenggen, WANG Feng, LIU Jixin
, Available online  , doi: 10.3969/j.issn.0258-2724.20250320
Abstract:

To accurately quantify and predict the internal pore pressure of tunnel lining concrete at the critical moment of spalling under high fire temperature, fire exposure tests on tunnel lining concrete were conducted. The evolution laws of the transient temperature field and pore pressure of lining concrete specimens during the fire process were obtained. Based on the tests, a “quasi-steady-state” pore pressure inversion model for high-temperature spalling of lining concrete was constructed, using experimental temperature data as the key driver. The reliability of the model was validated by comparing it with experimental results, and the spatial distribution characteristics of pore pressure at the moment of high-temperature spalling were revealed. The results show that for typical tunnel lining concrete specimens with an initial moisture content of 3.50%, at the critical moment of spalling under the HC fire condition, the pore pressure value ranges from 0.79 MPa to 0.96 MPa, with an average value of 0.88 MPa. The pore pressure value predicted by the model ranges from 0.91 MPa to 1.13 MPa, with an average of 0.99 MPa. The peak pore pressure predicted by the model agrees well with the experimental data, with errors ranging from 2.08% to 28.41% and an average error of 15.23%, which is lower than 20%. The pore pressure along the depth direction from the heated surface at the high-temperature spalling moment exhibits a “hump-shaped” spatial distribution, which can be resolved into the thermal cracking zone, the pore pressure concentration zone, and the thermodynamic equilibrium zone.

Optimization of Wheelset Tread for High-Speed Trains Considering Curve Passing Performance of LMA/CHN60N Wheel-Track
ZHANG Lei, ZHENG Jie, ZHAO Yafei, SUN Yuanbo, LI Pengfei
, Available online  , doi: 10.3969/j.issn.0258-2724.20250186
Abstract:

To improve the passing performance of high-speed trains on curved sections, reduce wheelset wear, and extend service life, the LMA/CHN60N wheel-track combination was taken as the research object. Six dimensions of the wheel tread in close contact with the track when the train passes through the curved track section were set as design variables, and the maximum lateral displacement, maximum climb, and maximum contact stress between the wheel and track on both sides were defined as objective functions. The response surface optimization design of the wheel tread dimensions was conducted, and the transient dynamics analysis of the optimization scheme was performed. The results indicate that compared with existing wheelsets of high-speed trains, the maximum lateral displacement, maximum climb, and maximum contact stress of the optimized outer wheels decrease by 1.74%, 1.34%, and 4.49%, respectively, and the maximum lateral displacement, maximum climb, and maximum contact stress of the inner wheels decrease by 0.41%, 1.53%, and 3.92%, respectively. Under the premise of ensuring other performances of the train, the optimization scheme improves the passing performance of the vehicle on curved track sections and can reduce wheelset wear; moreover, the indicators of straight passing performance, safety, smoothness, and durability are all within allowable ranges, thereby extending its service life.

Finite Element Numerical Simulation and Experimental Verification of Frequency Characteristics of Wind-Induced Vibration of Catenary
MIN Yongzhi, ZHANG Zhihui, HU Yanwen, WANG Guo, ZHANG Tingrong, ZHAO Shanpeng
, Available online  , doi: 10.3969/j.issn.0258-2724.20250497
Abstract:

Wind-induced vibration of catenary is a typical disaster of railway power supply systems in windy areas. It is difficult for existing theoretical analyses of vibration characteristics and finite element software to accurately simulate field conditions, which makes it difficult to verify their calculation results with field measurement data, and is not conducive to improving the calculation accuracy of subsequent nonlinear models. Therefore, a small-scale catenary model was built. Firstly, a finite element numerical calculation program integrating the large-deformation geometric nonlinearity of catenary and the stretching-contracting stiffness nonlinearity of suspension string was written in MATLAB, and the vibration modes of the small-scale catenary model under lateral excitations of different frequencies were solved. Secondly, an experimental platform with adjustable lateral excitation frequency was built, and the vibration modes under horizontal airflows at various frequencies were measured by a laser displacement sensor. Finally, the vibration waveforms obtained by the two methods were comparatively studied. The calculation results indicate that in a certain frequency range, the vibration frequency of the model is the same as the lateral excitation frequency, and when the lateral excitation frequency is equal to the fundamental frequency of the model, the maximum horizontal displacement distance reaches the maximum. Under the lateral excitations of different frequencies, the trajectories of the measuring point are different: When the frequency is less than the fundamental frequency, the trajectory is mostly 8-shaped; when the frequency is greater than the fundamental frequency, it is mainly a flat oval; when the frequency is equal to the fundamental frequency, it is almost a horizontal straight line. By comparing the test results, it is revealed that the finite element numerical calculation method can accurately calculate the fundamental frequency of the model, as well as the response frequency, vibration mode, and variation trend of amplitude with frequency under excitation. However, there is a significant deviation between the calculated amplitude values and the test results. Therefore, the calculation and prediction of key amplitude values need to be corrected by combining with experiments.

Voxel-Based Multi-Feature Fusion Modeling Method for Fractured Surrounding Rock in Tunnels
ZHU Qing, JIANG Zhaoyi, WU Haoyu, DING Yulin, WANG Qiang, ZHENG Weipeng
, Available online  , doi: 10.3969/j.issn.0258-2724.20250079
Abstract:

Fractured surrounding rock zones encountered during tunnel construction can easily trigger geological hazards such as collapses, roof falls, water inrushes, and mudbursts. To comprehensively improve the accuracy and efficiency of dynamic optimization design and multi-field coupled risk analysis, a modeling method for fractured surrounding rock integrating multi-source short-range advanced geological forecasting data was proposed. Firstly, a machine learning algorithm was employed to accurately extract the feature information of fractured surrounding rock from tunnel face images and ground-penetrating radar images. Secondly, the geometric and semantic features of the two types of data were registered and fused within a unified spatiotemporal reference framework to generate a unified semantic point set. Finally, a dense voxel grid was constructed, and the semantic point sets within the neighborhood of each voxel were traversed to build a three-dimensional fractured surrounding rock voxel model using the inverse distance weighting interpolation algorithm. A typical long and deep-buried tunnel was selected as a case to conduct experimental evaluation, and the voxel model and registered deepened blast hole data were subjected to voxelization processing and attribute alignment comparative analysis. The results indicate that the average classification accuracy of the constructed voxel model is 83.31% in fractured risk zones and 82.89% in non-risk zones; the model demonstrates high consistency with the blast hole data in terms of spatiotemporal distribution at the voxel scale, effectively characterizing the spatiotemporal evolution characteristics of the fractured surrounding rock during the dynamic excavation process.

Modified p–y Curve by Cone Penetration Test for Predicting Lateral Load–Displacement Behavior of Helical Steel Piles
SHAO Kang, LIAO Wenjun, ZHAO Kewei, SHAO Guoxia, ZHANG Junwei, YANG Shuisheng
, Available online  , doi: 10.3969/j.issn.0258-2724.20250400
Abstract:

In order to investigate the lateral load–displacement behavior of helical steel piles in sand, by considering the nonlinear pile-soil interaction, three p-y curves (p denoted soil reaction per unit length, and y denoted lateral displacement of the pile shaft) were selected to predict and compare lateral bearing behaviors. First, the existing classical sand p-y curve and two CPT-based p-y curves (power function and hyperbola) were analyzed, and their shortcomings were identified. Second, a linearly varying correction factor was introduced within the helical plate influence zone to modify three p-y curves, thereby incorporating the lateral bearing action of the helical plates. Third, the finite-difference method was employed to solve the bending differential equation of the laterally loaded helical steel pile, with pile shaft displacements obtained through iterative solution. Finally, by using the field tests of lateral bearing capacity for helical steel piles conducted in the literature, the prediction results of the three modified p-y curves were subjected to error and statistical analysis against the actual field measurements. The results show that the hyperbola-modified p-y curve produces load–displacement predictions closest to the measurements, with an average predicted-to-measured ratio of 1.04 and a goodness of fit of 98.5%. Among the three methods, the hyperbola approach yields the mean absolute percentage error of 11.63%, compared to 36.32% for the power function method and 24.25% for the classical p-y curve method. Its 95% confidence interval for mean prediction error spans from −6.0% to 11.9%, and the average error of 90% confidence interval spans from −4.6% to 10.5%, which verifies the reliability of the method. These findings provide a robust reference for accurately predicting the lateral load–displacement behavior of helical steel piles.

Life Prediction and Sensitivity Analysis of Recycled Concrete Under Dry, Wet, Freeze-Thaw, and Salt Intrusion
WANG Chenxia, ZHANG Weize, CAO Fubo, HAN Huichao, SU Tian, XIAO Shengxian, QIN Shijie
, Available online  , doi: 10.3969/j.issn.0258-2724.20240403
Abstract:

To address the insufficient durability of recycled concrete structures in the Hetao irrigated saline soil region, a typical geomorphological environment in the area was simulated, and concrete specimens were placed in a dry, wet, and freeze-thaw environment to investigate the erosion regularity and mechanism of sodium sulfate. Concrete type, water-cement ratio, and sodium sulfate solution concentration were employed as experimental variables. The quality evaluation parameter ƞa and the dynamic elastic modulus evaluation parameter ƞb were selected as key evaluation parameters to establish a comprehensive damage evaluation index. Based on this, a two-parameter Weibull reliability analysis was conducted, and a sensitivity analysis of each influencing factor was performed using grey relational grade. The results show that ordinary concrete exhibits better erosion resistance than recycled concrete; after the fourth cycle, the ƞb value of ordinary concrete is four times that of recycled concrete. Furthermore, appropriately reducing the water-cement ratio can improve the erosion resistance of recycled concrete. Under the influence of multiple factors, the damage to recycled concrete is caused by the cumulative effects of sodium sulfate crystallization expansion and chemical corrosion. The sodium sulfate with a higher concentration leads to faster damage development. After three cycles, the ƞb value of the RC-0.46-0% group is still 0.051, while the ƞb value of the RC-0.46-10% group is −0.066, indicating that the specimen has already been destroyed. The Weibull distribution function can effectively describe the erosion damage process of recycled concrete. By establishing an accelerated life reliability function, the service life of the specimen can be intuitively predicted, and the prediction results are consistent with the experimental results. According to the analysis results of the grey relational grade, after four cycles, the ƞa value of the specimen with a water-cement ratio of 0.53 is 0.36 and 0.65 lower than that of the specimens with water-cement ratios of 0.46 and 0.39, respectively, while the ƞb values of the latter two groups of specimens are basically the same, Therefore, the water-cement ratio is the main factor influencing the durability of recycled concrete.

Analysis of Influence of Urban Rail Transit Metro Lines on Direct Current at Neutral Point of Transformers
TANG Zeyang, MEI Xin, RUAN Ling, DENG Wanting, CUI Yibo, LI Yiwen
, Available online  , doi: 10.3969/j.issn.0258-2724.20250658
Abstract:

To investigate the influence of stray currents from urban rail transit on power grid transformers, the characteristics and causes of the neutral point current of transformers along metro lines were analyzed based on long-term monitoring data from the Hubei Observation and Research Station of Earth Stray Currents and Energy Equipment Safety. First, the balance and difference in the neutral point current among transformers within the same power grid were statistically analyzed. Second, the current characteristic quantities in different power grids were compared based on the wavelet transform. Finally, the differences in current waveform characteristics before and after metro operations, as well as during different operation periods, were examined. The results indicate that stray currents generated by urban rail transit cause an imbalance in the amplitude of the neutral point current among different transformers in the same grid, with a maximum difference of over 10 times. Denser metro lines around a substation result in a larger amplitude of the neutral point current of its transformer. When metro lines and power grids are geographically closely interconnected, the similarity of neutral point current waveforms among different transformers in the grid is high, with a correlation coefficient of over 0.8. The direct current generated by some metro lines flows only between local stations in the grid, which significantly reduces the similarity of current waveforms within the grid. Within the same grid, the waveform similarity is high when the wavelet scales of the neutral point current waveforms of transformers are close; however, the wavelet characteristics of the neutral point current of transformers from different grids are not comparable. When there are obvious differences in the alignment or operation of metro lines, the influence of different metro lines on the neutral point current of transformers can be identified based on wavelet scales.

Study on Optimization of Return Current Nodes in Active Return Traction Power System
SUN Xiangxuan, YANG Xiaofeng, LIU Ke, ZHENG Qionglin, ZHAO Runda, WU Tianyi, CHEN Lu
, Available online  , doi: 10.3969/j.issn.0258-2724.20260081
Abstract:

The rail potential and stray current in urban rail transit are mitigated by the active return traction power system (AR-TPS) by transferring the running rail current to the return cables through the return current node (RCN). To improve the mitigation performance of the AR-TPS, an optimization layout method of RCNs was proposed. First, a two-layer distributed parameter model and an RCN current transfer model were established to reveal the influence laws of RCN number and layout position on the current transfer capability. Then, based on engineering constraints and node functions, the number of RCNs was simplified to two, and the layout positions were optimized near the midpoint of the section with the optimal mitigation performance as the objective. Finally, a simulation verification was conducted based on the typical train operation curve. The results indicate that the outer RCN working in coordination with the active resistance converter is the primary path for transferring the return current to the return cables. The mitigation performance of the AR-TPS continuously improves as the RCN approaches the midpoint of the section. When the conventional traction power system is adopted, the maximum rail potential is 34.86 V; the maximum stray current is 2.89 A, and the total leakage charge is 134.00 C. When the AR-TPS with uniformly distributed RCNs before optimization is adopted, the above indicators are reduced to 14.84 V, 0.94 A, and 55.70 C, and the suppression rates are 57.42%, 67.47%, and 58.43%, respectively, compared with the conventional traction power system. When the AR-TPS optimized by the proposed method is adopted, the above indicators are further reduced to 8.66 V, 0.66 A, and 38.67 C, and the suppression rates are increased to 75.16%, 77.16%, and 71.14%.

Localization Method for Metro Stray Current Leakage Points Based on Tellegen’s Theorem
ZHOU Qi, ZHAN Ying, WANG Aimin, HE Bin, ZHAO Liping, LIN Sheng
, Available online  , doi: 10.3969/j.issn.0258-2724.20260013
Abstract:

To locate the metro stray current leakage points from the source, an equivalent resistance network model of the metro traction power supply system was first established. According to the characteristic of invariant topology of the equivalent resistance network before and after the rail insulation damage, the relationship between the rail-to-ground transition resistance and the branch voltages and currents of the traction substations and trains was derived based on Tellegen’s theorem. Secondly, by considering the spatiotemporal correspondence between train operations and the related electrical quantities of the traction substations and trains, a multi-time-section solving equation system for the rail-to-ground transition resistance was constructed. The locating interval was progressively subdivided by using the bisection method to increase the dimensionality of the equation system, and a genetic algorithm was combined to iteratively optimize the equation system to obtain the distribution of the rail-to-ground transition resistance along the line. Finally, the calculation results of the rail-to-ground transition resistance were compared with the identification threshold of stray current leakage points to realize the localization of stray current leakage points. Simulation verification was conducted based on actual metro line parameters and measured track inspection train data, and the applicability of the proposed localization method was analyzed from three aspects: stray current leakage point location, leakage severity, and leakage section length. The results show that under the three types of simulation scenarios of different leakage point locations, different leakage severities, and different leakage section lengths, the maximum localization errors are 9.375 m, 0 m, and 6.875 m, respectively, and the relative errors are all below 0.4%. It can be seen that the localization results of the proposed localization method are basically consistent with the preset locations. In summary, the method can realize the localization of the section where the stray current leakage points are located and can provide a theoretical basis for stray current mitigation in scenarios such as rail fastener damage and localized insulation degradation.

Calculation and Analysis of Cable Metal Sheath-to-Ground Potential Under Influence of Metro Stray Current
LIN Sheng, MAO Gefan, ZHOU Qi
, Available online  , doi: 10.3969/j.issn.0258-2724.20260071
Abstract:

To analyze the influence of metro stray current on the cable metal sheath-to-ground potential, by considering its dual intrusion paths through the metal sheath and the ground, electrical quantity analysis models for the cable metal sheath under double-end grounding and cross-bonded grounding methods were respectively established. Based on the superposition theorem, the induced voltage and the voltage drop generated by stray current were superposed, and a calculation method for the sheath-to-ground potential at any position of the cable metal sheath under the two grounding methods was derived. The calculated values were verified by the measured current of a 110 kV cable under the cross-bonded grounding method, and the calculation results of the sheath-to-ground potential under the two grounding methods were verified based on the PSCAD simulation platform. On this basis, the proposed models were adopted to analyze the influence of stray current amplitude, frequency, and soil resistivity on the metal sheath-to-ground potential. The research results indicate that the trend of the calculated value of the cable metal sheath current is basically consistent with that of the measured value, with a relative error of 3.08%; the calculation error of the sheath-to-ground potential at the cross-bonded joint is 0.86%, and the calculation error of the sheath-to-ground potential at the midpoint position of the metal sheath under the double-end grounding method is 1.13%; under the double-end grounding method, the cable metal sheath-to-ground potential always remains within the safety limit, and when the stray current increases from 0 A to 300 A, the sheath-to-ground potential at the midpoint position of the metal sheath increases from 1.10 V to 30.25 V; under the cross-bonded grounding method, when the stray current exceeds 162 A, the sheath-to-ground potential at the cross-bonded joint exceeds the 50 V limit; the stray current frequency has no obvious influence on the amplitude of the cable metal sheath-to-ground potential under the two grounding methods but only makes the peak appear more frequently; the increase of soil resistivity raises the cable metal sheath-to-ground potential, and the raising amplitude under the double-end grounding method is more significant than that under the cross-bonded grounding method.

Event Causality Identification Based on Large Language Model-Constructed Graph Networks
PAN Lei, YUAN Hongxiao, ZHONG Zhun, LIAO Hongzhou, YANG Ruijia
, Available online  , doi: 10.3969/j.issn.0258-2724.20240484
Abstract:

To enhance the accuracy of document-level event causality identification, the implicit knowledge within large language models was first leveraged to filter out events related to the target event from the document, constructing a candidate event set. Next, the candidate event set was organized into a fully connected event relationship graph, which was then subjected to conditional constraints. These constraints simplified the fully connected graph into a constrained event relationship graph, reducing noise propagation from irrelevant events. Finally, a self-attention mechanism was used to compute the influence of any node on other nodes in the graph network, and the model was trained with a binary classification loss function incorporating focal loss, further mitigating false positives in causality identification. Experimental results show that in the causality identification tasks, the model achieves a precision of 77.3% for intra-sentence event pair causality identification on the Causal-TimeBank dataset. On the EventStoryLine dataset, it achieves a precision of 75.2% for intra-sentence event pair causality identification, an F1-score of 60.6% for inter-sentence event pair identification, and an F1-score of 59.6% for document-level event pair identification.

Seismic Resilience Enhancement of Railway Beam Bridges Based on Viscoelastic Dampers
ZHAO Canhui, LIANG Jinhao, LU Hao, YIN Weitao, FENG Junji, KANG Wei
, Available online  , doi: 10.3969/j.issn.0258-2724.20250318
Abstract:

As the seismic design philosophy of bridges gradually shifts from performance-based to resilience-oriented, the current research focus lies not only in ensuring structural safety under the action of earthquakes, but also in paying more attention to the preservation and recovery capacities of post-earthquake traffic function. Railway bridges serve as critical nodes in transportation networks, and their post-earthquake traffic capacity directly influences the operational efficiency of the lines and the recovery process of earthquake-stricken areas. Based on this, a shear-type viscoelastic damper was proposed to decrease track damage risk by controlling the relative transverse displacement of main beams, thereby enhancing the bridge’s post-earthquake traffic function. By taking a typical five-span railway simply-supported bridge as a case study, a comparative analysis was conducted on the traffic function vulnerability, post-earthquake function recovery capacity, and seismic resilience in the three constraint conditions of no constraint (NC), steel restrainer (SR) only, and SR combined with the shear-type viscoelastic damper (DP). The analytical results demonstrate that proposing the concept of a “component recovery weight coefficient” provides a quantitative index for determining the sequence of component repairs. Additionally, the proposed shear-type viscoelastic damper can significantly control the relative transverse displacement of main beams at beam ends, effectively reducing track alignment irregularity under the strong action of earthquakes and consequently enhancing the bridge’s post-earthquake traffic function. The addition of the viscoelastic damper to SR can notably shorten the post-earthquake repair time. In particular, under higher ground motion intensities, such as peak ground acceleration (PGA) of 0.8g, the repair time can be reduced by up to 9.1 days compared to the NC system. Furthermore, the DP configuration demonstrates significantly better post-earthquake residual traffic function and seismic resilience indices than the SR and NC configurations at various PGA levels. For instance, under PGA=0.4g, the post-earthquake residual traffic function improves by 20.0% and 56.9%, and seismic resilience increases by 13.1% and 34.4% respectively, confirming the remarkable effectiveness of this device in enhancing both post-earthquake traffic function and seismic resilience of bridges.

Construction of Model for Carbon Emissions and Analysis of Quantitative Characteristics for Whole Life Cycle of Railway Tunnel Engineering
MENG Fanqiang, GONG Xun, YE Zilin, ZHAO Liuhui, GONG Hao, WANG Dongmei, GONG Zhengjun
, Available online  , doi: 10.3969/j.issn.0258-2724.20250279
Abstract:

Given the complex construction processes and significant procedural variations in railway tunnels, a model for accounting carbon emissions covering the engineering materialization phase and the operation and maintenance phase was established to accurately assess carbon emission characteristics throughout the whole life cycle of railway tunnel engineering and support specialized analysis and emission reduction strategies. Combined with the unit rate method by procedure and the carbon factor method by material, quantitative accounting and variance analysis of carbon emissions were conducted in a case study. Furthermore, Pareto analysis was employed to identify high-emission sectors and key influencing factors, realizing the decomposition and hierarchical management of carbon emissions. The results indicate that the total carbon emissions during the whole life cycle of the case tunnel engineering are 68.107 5 × 104 tCO2e, with the engineering materialization phase generating 59.404 0 × 104 tCO2e (accounting for 87.2%) and the operation and maintenance phase generating 8.703 5 × 104 tCO2e (accounting for 12.8%). During the engineering materialization phase, material production and on-site construction processes contribute up to 95.76% of the carbon emissions, among which the carbon emission intensity of drill-and-blast section construction is 13 843 tCO2e/km, and that of shield section construction is 39.8% of the former. During the operation and maintenance phase, the maintenance of ventilation systems, waterproofing and drainage systems, and reinforced concrete precast pipes accounts for 73.54% of the carbon emissions. In the whole life cycle, Class A carbon-emitting materials include ordinary cement 42.5 grade (high-performance concrete), ordinary cement 42.5 grade, and ribbed steel bars; Class A carbon-emitting machinery includes axial flow fans, rail-type internal combustion locomotives, and dump trucks; Class A carbon-emitting processes include the lining of shield section, the support of drill-and-blast section, and the support of shield section.

Adaptive Sliding Mode Control of Shift Clutch in Hydraulic Mechanical Transmission Based on Linear Quadratic Optimization
CAO Fuyi, QIU Futao, YUAN Tianqi, YAN Xianghai
, Available online  , doi: 10.3969/j.issn.0258-2724.20250285
Abstract:

To improve the torque transmission control accuracy of the shift clutch in hydraulic mechanical transmission (HMT) and enhance its shift quality, an adaptive sliding mode control method based on linear quadratic optimization was proposed. Firstly, the shift dynamics model of HMT and the numerical model of the clutch were constructed. Secondly, the optimal torque transmission trajectory of the shift clutch was solved under the target functional based on the linear quadratic optimization model. Combined with the dynamic friction coefficient model, the trajectory variable was converted into the expected oil pressure trajectory considering the time-varying effect of the lubricating oil film and the dynamic contact characteristics of the friction interface. Finally, the adaptive sliding mode control algorithm was used to track and control the expected oil pressure trajectory. Simulation results show that compared with the direct application of the linear quadratic control scheme, the proposed algorithm effectively improves the torque transmission control accuracy of the shift clutch, shortens the shift time by 13.8%, reduces the sliding friction work by 11.2%, and decreases the maximum shock by 39.1%. Experiments further confirm that the algorithm not only effectively improves the torque transmission control accuracy of the shift clutch but also improves the shift quality of HMT. The research results can provide a reference for the formulation of shift control strategies in the engineering application of HMT devices.

Distribution and Migration Characteristics of Volatile Organic Compounds at Polluted Sites of Chemical Plants and Their Concentration Prediction Models
MA Zhiqiang, ZHAO Lin, XIAO Man
, Available online  , doi: 10.3969/j.issn.0258-2724.20240223
Abstract:

To clarify the distribution and migration characteristics of volatile organic compounds (VOCs) in organic pollution sites, 120 soil and 48 soil gas sampling holes were arranged at a site where a chemical plant was once located, and 97 pairs of soil/soil gas samples were collected from 14 holes in the heavily polluted area. The pollutant concentrations in the samples were determined through laboratory analysis. The distribution characteristics of the pollutants were analyzed using Sufer software, and the distribution and migration characteristics of trichloroethylene (TCE) were interpreted in combination with the stratigraphic lithology. The applicability of two models for predicting the content of soil gas VOCs was analyzed with predicted and measured values. The results show that as the hole depth increases, the TCE concentration in the soil in the heavily polluted area decreases with an exceedance rate of 50.4%. The maximum value is in the clayey silt layer (1610 mg/kg), exceeding the standard by 2300 times. The TCE concentration in soil gas first increases and then decreases, with the maximum value occurring in the sandy soil layer (643 mg/m3), exceeding the standard by 546 times with an exceedance rate of 97.9%. The TCE concentration in soil and soil gas outside the heavily polluted area first increases and then decreases with soil depth, and the exceedance rate is 6.6% and 54.9%, respectively. The horizontal migration of TCE is weaker than its vertical migration. The area that exceeds the standard first increases and then decreases with soil depth, and the average exceedance area is 7.3 × 103 m2. TCE in soil gas spreads rapidly from the heavily polluted area to the surrounding areas. The exceedance area covers the entire chemical plant area and pollutes the area outside the chemical plant. The average exceedance area is 2.4 × 104 m2. Based on the soil TCE concentration, both the linear model and the dual equilibrium desorption (DED) model can predict the soil gas TCE concentration. The proportion of data points with a prediction error of one order of magnitude or less for the two models is 78.6% and 71.4% for the fill and sandy soil, respectively. Regression analysis shows that the DED model predicts soil gas TCE concentrations more close to the actual values for fill (irreversible adsorption degree f = 0.1), silty clay (f = 1.0), and clayey silt (f = 0), while the linear model and DED model yield relatively close predictions of soil gas TCE concentrations for sandy soil. When a pollution assessment and secondary development of similar sites are conducted, the concentrations of soil gas pollutants should be monitored or predicted to avoid underestimating the degree and scope of site pollution.

Multi-Physical Field Coupling Calculation of Gas-Insulated Switchgear Bus Temperature Distribution Based on Radiative Cooling Technology
GAO Guoqiang, ZHOU Xianzhi, ZHANG Chuan, LIU Yijie, GUO Yujun, WU Guangning
, Available online  , doi: 10.3969/j.issn.0258-2724.20240568
Abstract:
Objective

In the outdoor environment of high solar irradiance, the long-term operation of gas-insulated switchgear (GIS) equipment will lead to uneven distribution of bus tube surface temperature, which will cause cracking at the fixed support of GIS. Radiative cooling technology has a good application prospect in the aspects of energy saving, cooling, and temperature uniformity, but there are few reports in the electrical field. Therefore, it is very important to study the radiative cooling technology to solve the influence of the uneven temperature distribution of the GIS bus barrel in an outdoor high solar irradiance environment.

Method

In this paper, the influence of gas axial flow and solar radiation in GIS was considered, and a multi-physical field coupling calculation method of GIS bus temperature distribution based on radiative cooling technology was proposed. Firstly, a multi-physical coupling model of the three-dimensional electromagnetic, thermal, and flow-radiation field of the GIS bus was established, and the electromagnetic loss of the conductor and the shell was calculated in the electromagnetic field; then, the electromagnetic loss was introduced into the subsequent three physical fields. The subsequent multi-physical field calculation results in turn affected the conductivity and changed the electromagnetic loss, so as to carry out the coupling. The electromagnetic field and thermal field were bidirectional coupling through the resistance effect of electrical conductivity; the flow field and thermal field were bidirectional coupling through the gas density and other parameters, and the radiation field and thermal field were unidirectional coupling through the temperature change caused by solar radiation, and the emissivity of the object did not change with the thermal field. Secondly, the temperature distribution of the GIS bus was compared and analyzed under the conditions with or without solar radiation and cooling with or without radiation, and a study on the all-day cooling effect of radiative cooling materials on the GIS bus when the external parameters such as solar radiation intensity, solar radiation angle, and ambient temperature change with time in a day was carried out. Then, the accuracy of the multi-physics simulation was verified by a simplified experimental model. Finally, the effects of ambient temperature, solar radiation intensity, and solar radiation angle on GIS, bus temperature, and the cooling effect of radiative cooling materials were discussed through simulation calculation and theoretical derivation.

Result

Under the influence of solar radiation, the maximum temperature difference between the upper and lower surfaces of the outer shell of outdoor GIS is 14.12 ℃, and the overall temperature of GIS increases by 10–15 ℃ compared with that without considering solar radiation. The temperature distribution of GIS after using radiative cooling materials is roughly the same as that without considering solar radiation, or in other words, the effect of solar radiation can be offset by using radiative cooling materials, reducing the maximum temperature difference between the upper and lower surfaces of the shell surface by 51.56%. In one day, the maximum cooling effect of radiative cooling materials is 14.6 ℃; the average cooling effect is 7.42 ℃; the average temperature difference between the upper and lower sides of the shell surface is reduced from 10 ℃ to 5.16 ℃, and the soaking performance is increased by 48.4%. The experimental and simulation errors of the conductor and shell are less than 10% through the simplified GIS model. The temperature of the GIS conductor and the outer shell increases approximately linearly with the increase of ambient temperature and light intensity. The slope of the curve of the conductor and shell temperature changing with the light intensity is smaller than that of the uncoated conductor and shell. In other words, the use of radiative cooling materials can reduce the average temperature rise of the conductor and shell so that the temperature distribution is more uniform. At the same time, the radiative cooling material can reduce the temperature between the conductor and the shell and the average temperature difference between the conductor and the shell. The solar radiation angle has little influence on the overall temperature distribution of GIS. The reason is that the bus bar of GIS is cylindrical, and the area exposed by the sun in each direction is roughly the same, but it has a significant influence on the temperature at the intercept point. The temperature at the intercept point is the highest at the vertical incidence and decreases at the bilateral incidence.The cooling power of radiative cooling materials mainly depends on the emissivity of 8–13 um band, and the power amplitude caused by other bands is very small. In order to improve the size of the cooling power, the emissivity of the material in the 8–13 um band can be appropriately increased, so as to improve its cooling effect. When the ambient temperature is constant, the cooling power decreases with the increase of solar radiation intensity. When the solar radiation intensity is constant, the cooling power increases with the increase of the ambient temperature. In the process of the change of the solar radiation angle of 0°–90°, a greater solar radiation angle means lower cooling power, and the size of the cooling power is proportional to the cooling difference; a greater cooling power indicates a greater temperature than the ambient temperature reduction amplitude.

Conclusion

After the application of radiative cooling material on the GIS bus barrel, the temperature of the GIS bus bar decreases, and the soaking performance of the bus barrel improves, which is of great significance for the calculation of outdoor GIS bus bar temperature and the improvement of temperature heterogeneity. However, the current research only applies the ideal condition that the emissivity of the radiative cooling material does not decrease with time. In practice, there will be some problems of contamination attached to the surface of the material or the service life of the material itself, so it is necessary to conduct more research on the application test on the physical GIS and the test of material emissivity change after long-term operation in the future.

Seismic Performance of CFSST Column-H-Shaped Steel Beam Joints with External Stiffening Rings
GAO Chunyan, WANG Jingjie, WANG Yanbin, WEN Yang, QI Jinliang, ZHOU Zhaodi
, Available online  , doi: 10.3969/j.issn.0258-2724.20250307
Abstract:

Quasi-static tests were conducted on three types of concrete-filled square steel tube (CFSST) column-H-shaped steel beam joints with external stiffening rings featuring different structures to investigate the seismic performance of prefabricated CFSST column—H-shaped steel beam joints with external stiffening rings. The bearing capacity, failure processes and modes, stress distribution in joint zones, and various seismic performance indicators were studied. The experimental results show that joint failures all occur near the beam ends, characterized by buckling and fracture failure of the beam sections outside the diaphragm. The hysteresis curves of the joints all exhibit a “spindle shape” with significant pinching effects. The loading process is essentially symmetric in the positive and negative directions, with a sliding segment of zero rotational stiffness observed in the curves. The ductility coefficient of joints ranges from 2.55 to 4.20, and the ultimate angle of rotation is between 0.032 rad and 0.062 rad, both exceeding the requirements specified in seismic codes. Compared to the integral external stiffening rings, the separate external stiffening rings exhibit increases of 60.7% in ductility and 209.0% in cumulative energy dissipation. Meanwhile, the cantilever energy-dissipating external stiffening rings show improvements of 38.1% in bearing capacity, 64.7% in ductility, and 400.3% in energy dissipation capacity over the integral connection. This indicates that the structural configuration of “separate multi-layer force transfer” and “dual-stage energy dissipation via cantilever segments and energy-dissipating cover plates” can effectively optimize the seismic performance of joints. There is sound agreement between the simulated and experimental hysteresis curves and failure modes, validating the numerical analysis reliability.

Noise Reduction Performance of 3D Printed Special-Shaped Noise Barriers for High-Speed Railways
ZHANG Xun, DONG Chaoxi, HU Guoqing, WAN Xing
, Available online  , doi: 10.3969/j.issn.0258-2724.20250179
Abstract:

To study the noise reduction performance of 3D printed special-shaped noise barrier for high-speed railways, the geometric models of the noise barrier were established by using three-dimensional modeling software, and the noise reduction model of a special-shaped noise barrier for high-speed railways under three sound source modes was established. Based on the acoustic indirect boundary element method, the optimal design parameters of a special-shaped noise barrier were obtained by comparing the noise reduction effect, sound field distribution, spectrum characters, and insertion loss characteristics of the special-shaped noise barrier under different wave distances, wave heights, and barrier heights. The effects of different barrier-line center distances and train speeds on the noise reduction effect of a special-shaped noise barrier were compared and analyzed. The results show that the special-shaped noise barrier has a significant noise reduction advantage for the middle and far field points, the field points in the sound shadow area, and the middle and high frequency noise above 630 Hz. Compared with that of the vertical noise barrier, the additional noise reduction gain of the special-shaped noise barrier can be obtained when the optimized wave distance is in the range of 3.2–4.0 m, and the additional noise reduction effect can reach 1.7 dB (A). In the case of a certain wave distance, with the increase of wave height, the noise reduction effect of the special-shaped noise barrier is enhanced, and the maximum can reach 15–17 dB (A). The additional insertion loss of the special-shaped noise barrier is 1.8 dB (A) when the wave distance is 4 m, and the wave height is 600 mm. When the special-shaped noise barrier height increases by 1 m, the sound pressure level of the four standard field points decreases by 1.6–4.8 dB (A), which mainly acts on the range of 5–10 m near the barrier area and the high frequency cut-off area in space. The additional noise reduction is attenuated from the area to both sides. In the range of 3–4 m, the barrier-line center distance has little effect on the noise reduction level of the special-shaped noise barrier, and the recommended value of the parameter is 3.4–3.6 m. As the train speed increases, the insertion loss of the special-shaped noise barrier decreases. When the train speed is increased from 300 km/h to 340 km/h, the insertion loss of the field point under different working conditions is reduced by 0.17–0.34 dB (A) and 0.55–0.63 dB (A), respectively.

Influence of Grouting Point Layout Angles on Deformation of Shield Tunnel under Load
HUANG Dawei, ZENG Xiaotao, GENG Daxin, JIANG Qiangbo, NIE Xianxiu, YOU Chenkang
, Available online  , doi: 10.3969/j.issn.0258-2724.20250336
Abstract:

Aiming at the problem of poor controllability of deformation recovery in the process of transverse elliptic deformation of shield tunnel treated by traditional stratum grouting, the experimental study on transverse elliptic deformation of shield tunnel treated by capsular bag grouting was carried out innovatively. Two different capsular bag layout schemes were adopted in the experiment. The analysis of the test results shows that the tunnel has a vertical elliptical deformation near the grouting point, and the additional earth pressure, tunnel diameter and deflection deformation near the grouting point change the most, showing a gradual decreasing trend along the two ends of the tunnel. Under the same amount of grouting, vertical grouting can better control the transverse elliptic deformation of the existing tunnel, and the side bottom grouting can realize the recovery effect of the transverse elliptic deformation of the tunnel while uplifting the tunnel and effectively reducing the settlement of the tunnel during operation. When grouting at the bottom of the tunnel side, due to the influence of the gravity of the built tunnel, the deformation of the tunnel caused by grouting is lower than that caused by vertical grouting during the expansion process of the bag. Therefore, it is suggested to increase the grouting amount appropriately in the grouting process to achieve the expected recovery effect of transverse ellipse deformation. After grouting, the capsular bag presents a cylindrical shape with wide middle and narrow sides. The capsular bag grouting method can effectively limit the slurry diffusion path, and the slurry is squeezed and diffused in the capsular bag, so that the local pressure of the tunnel increases rapidly, and the recovery effect of transverse ellipse deformation is remarkable. Vertical grouting at the side of the tunnel has a significant effect on the recovery of transverse ellipse deformation, but excessive grouting pressure and other reasons can easily induce problems such as dislocation of shield tunnel segments. It is suggested that when the transverse ellipse deformation of the shield tunnel exceeds the limit, the grouting point should be arranged at the bottom of the tunnel side, that is, about 30°~60°away from both sides of the tunnel.

Seismic Performance of Hybrid Reinforced Columns with Built-in Spiral Stirrup Core Columns Under Composite Salt Erosion
WANG Tan, YANG Fan, ZHOU Zhijie, LIU Jinling, LI Ning, WANG Kaiqi
, Available online  , doi: 10.3969/j.issn.0258-2724.20240675
Abstract:

To enhance the durability of structural members while ensuring ductility, a new glass fiber-reinforced polymer (GFRP)-steel hybrid reinforced (HRBS) column with a built-in spiral stirrup core column was proposed. Four HRBS columns were subjected to quasi-static loading tests, among which two HRBS columns underwent composite salt dry-wet cycling tests. The failure process, ultimate failure mode, hysteresis curves, skeleton curves, energy dissipation capacity, performance degradation, and residual displacement of the HRBS columns before and after composite salt dry-wet cycling were investigated. The results indicate that HRBS columns demonstrate good seismic performance both before and after composite salt dry-wet cycling, and all specimens exhibit flexural failure. Under conventional environments, the yield load and peak load of HRBS columns increase by 25.77% and 28.68%, respectively, with an increase in the core column diameter. In a composite salt environment, increasing the core column diameter improves the load-bearing capacity, displacement ductility factor, strength degradation factor, energy dissipation capacity, overall stiffness, and self-centering ability of the HRBS columns. After composite salt erosion, the yield load and peak load of specimens with a 200 mm core column diameter decrease by 17.65% and 15.77%, respectively, while the energy dissipation capacity and displacement ductility factor increase by 14.41% and 32.61%, respectively. Therefore, by designing an appropriate core column diameter, HRBS columns ensure good durability and overall seismic performance in both conventional and corrosive environments.

Experimental Study on the Mechanical Behavior of Laminated-Rubber Bearing-Concrete Interface
SHAO Changjiang, CUI Haomeng, WANG Chunyang, JIANG Pengxu, ZHUANG Weilin
, Available online  , doi: 10.3969/j.issn.0258-2724.20250217
Abstract:

In small-to-medium span highway girder bridges in China, laminated-rubber bearings are typically positioned between the steel plate at the beam bottom and the concrete padstone on the pier/abutment. While relevant specifications indicate that relative sliding can occur at the bearing-steel plate interface, under certain conditions, sliding is more prone to occur at the bearing-concrete padstone interface. To investigate the mechanical behavior of this interface, reciprocating compression-shear tests were conducted under three levels of contact pressure: low (8 MPa), medium (10 MPa), and high (12 MPa). A three-phase deformation and failure mechanism was revealed, a bulging failure criterion was established, the evolution laws of the shear modulus and friction coefficient were quantified, and a friction coefficient attenuation model was developed. The results show that the bearing hysteresis curves are characterized by a parallelogram shape with internal crescent-shaped abrupt loops. Bidirectional damage transfer occurs at the bearing-concrete interface, evidenced by bulging of the rubber cover layer on the bearing and the presence of rubber debris adhered to the concrete surface. In the sliding friction stage, the shear modulus decreases with increasing contact pressure and consistently remains below the code-specified dynamic shear modulus of 1200 kPa. The interface friction coefficient, measured to be below the code-suggested values of 0.25 and 0.30, exhibits a negative correlation with the contact pressure. It is recommended to control the bearing contact pressure below 10 MPa in practical engineering to avoid the risks of a sudden drop in the friction coefficient and bulging failure. The proposed models can provide a reference for the seismic design of small-to-medium span bridges.

Analysis of Shear Strength Parameters of Waste SoilBased on Field and Laboratory Tests
ZHANG Kunyong, XUE Tangwei, LI Junji, LI Dalong, XIE Xin
, Available online  , doi: 10.3969/j.issn.0258-2724.20240607
Abstract:

At present, most of the studies on the shear strength parameters of waste soil are carried out based on laboratory waste soil samples, and the parameters obtained through laboratory tests are often quite different from the actual parameters, resulting in a large error in the calculation of waste landfill deformation. In view of the error problem caused by the laboratory test, by relying on the closure and restoration project of the waste landfill under construction, the influencing factors and change laws of the shear strength parameters of the waste soil were analyzed by carrying out the large-scale field direct shear test and the small-size laboratory direct shear test of the waste soil and making a comparison with the literature test results. The results show that the shear capability of waste soil is higher than that of traditional soil, and the shear stress continues to rise and tends to be stable with the increase of shear displacement. The density of waste soil is linearly and positively correlated with its shear strength. When the density increases from 750 kg/m3 to 950 kg/m3, the cohesion increases from 22.5 kPa to 34.5 kPa, and the friction angle decreases from 36.7° to 30.2°. The mean cohesion value (14.91 kPa) obtained by the small-size laboratory shear box test is 6.27 kPa lower than that of the large-size field test (21.18 kPa), and the mean friction angle (35.17°) is 6.61° higher than that of the field test (28.56°).

Vibration Characteristics of Riser Induced by Gas-Liquid-Solid Three-Phase Flow in Deep-Sea Hydrate Extraction Under VIV Effect
GUO Xiaoqiang, JIANG Zhefu, YANG Kelun, XU Jie, LV Junlin, LI Xinye
, Available online  , doi: 10.3969/j.issn.0258-2724.20240534
Abstract:
Objective

Natural gas hydrate, as a low-carbon unconventional energy source, is increasingly receiving attention from the scientific community. It is considered the most promising alternative energy source in the 21st century. At present, the extraction of deep-sea hydrates is mainly carried out through riser transportation. Due to the state instability of the hydrate, it is easy to decompose during migration, forming a typical vibration phenomenon induced by gas-liquid-solid three-phase flow, which can easily lead to nonlinear vibration of the riser. This nonlinear vibration mechanism is completely different from the vibration mechanisms induced by single-phase and two-phase flows. It is particularly important to reveal the vibration mechanism of the riser in gas-liquid-solid three-phase flow.

Method

Due to external ocean loads during the process of mining riser operations, it is easy to cause vortex-induced vibrations (VIV). The riser was subjected to the action of internal gas-liquid-solid three-phase flow of hydrate. Moreover, the upper end of the riser was subjected to the heave motion of the ocean platform. These factors require the use of multiple methods in combination for the nonlinear vibration model of the riser. Therefore, the gas-liquid-solid three-phase FIV model for deep-sea hydrate extraction riser was established using the finite element method, Hamiltonian principle, and energy method. The fluid structure coupling effect between ocean flow field and mining riser was achieved through the wake oscillator model. The coupling effect between internal multiphase flow and mining riser was achieved through additional mass, collision energy loss, and flow velocity changes. Simultaneously, a dynamic decomposition model for hydrates was established to identify changes in the content of gas, solid, and liquid phases. The platform was prone to six degrees of freedom motion under the action of wind, waves, and currents. However, during the process of hydrate extraction, the platform was most prone to heave movement. With the help of previous research, a platform heave motion model was established, which can effectively obtain the displacement boundary of the upper end of the mining riser. Then, in order to achieve a numerical solution of the model, the vibration model was solved using the combined iterative method of the incrementally applied Newmark-β method and Newton-Raphson method, obtaining the vibration response of the mining riser. Due to the excessive consideration of nonlinear factors in the model, the correctness and effectiveness of its numerical solution need to be rigorously verified. However, the actual vibration data on site could not be obtained. The simulation experiment of the water tank was particularly important. Therefore, by using the principle of similarity, a nonlinear vibration simulation experiment system for a hydrate mining riser under internal and external flow excitation was developed. The root mean square of the vibration displacement and amplitude frequency curves of the mining riser were experimentally measured. At the same time, the calculation parameters of the theoretical model were set exactly the same as the experimental parameters to compare the theoretical calculation results and the experimental test results. On this basis, the root mean square of the vibration displacement and the amplitude frequency curve of the mining riser were calculated. Comparing the experimental results with the model calculation results, the correctness of the established model was verified, and the comparison accuracy could be higher than 90%. By using frequency domain and time domain analysis methods, the influence of external environmental parameters and multiphase flow parameters on the nonlinear FIV response of the riser was explored.

Result

The results show that the vibration amplitude of the riser in the cross-flow (CF) direction is higher than that in the in-line flow (IL) direction, and the variation of internal flow parameters has a more significant impact on the CF vibration of the riser. This phenomenon indicates a close relationship between the internal gas-liquid-solid three-phase flow and the vortex-induced effect. The axial vibration of the riser mainly consists of two parts: One part is dominated by gravity and platform heave motion, exhibiting low-frequency and high-amplitude vibration characteristics, which may cause the strength failure of the riser. The other part is induced by internal and external flow field loads, exhibiting high-frequency and low-amplitude vibration characteristics, which may accelerate the fatigue failure of the riser. With the shear flow velocity increases, the displacement of the riser in the IL direction gradually increases, and the axial second-order main frequency amplitude increases accordingly. However, in the CF direction, the amplitude and root mean square of displacement show a decreasing trend, indicating that the increase in displacement in the IL direction has a certain inhibitory effect on the CF vibration. With the increase of internal output volume, the flow velocity of each phase increases, resulting in an increase in displacement and amplitude of the riser in the IL direction, as well as an increase in amplitude in the CF direction. This widens the frequency band, expands the displacement envelope range, and reduces the number of modes in the CF direction. As the particle size of hydrates increases, the solid-phase flow velocity decreases, and the collision energy loss increases, resulting in a decrease in the vibration amplitude of the riser in the IL direction, a narrowing of the vibration frequency band, and an increase in amplitude in the CF direction. The characteristic of more concentrated energy makes it possible to control it more effectively in the future. When the shear flow velocity reaches 1.4 m/s, the hydrate abundance reaches 80%, or the hydrate particle size reaches 7 mm, the vibration amplitude of the riser in two or three directions will significantly increase, and the vibration response of the riser will not change with the variation of other parameter values. At these specific values, the frequency of the interaction between the internal and external fluids on the riser tends to approach the natural frequency of the riser system, resulting in resonance phenomena. In the actual operation process, these parameter values should be avoided to ensure the safety and stability of the riser.

Conclusion

Once commercial exploitation of deep-sea hydrates is achieved, the nonlinear vibration model of deep-sea hydrate mining riser can effectively guide the parameter configuration in the later hydrate mining process. It can evaluate the safety and service life of deep-sea hydrate mining risers, ensuring the safe operation of commercial mining cycles of hydrate.

Study on Safety Threshold for Rail Corrugation in High-Speed Railway Lines Under Coupled Action of Wheel’s Concave Wear
QI Yayun, ZOU Rui, DAI Huanyun, CHEN Zhaowei, HE Xing, LI Liangxu
, Available online  , doi: 10.3969/j.issn.0258-2724.20250007
Abstract:

To investigate the influence of rail corrugation coupled with wheel concave wear on the dynamics characteristics of the high-speed train’s wheel-rail system and determine the safety threshold for rail corrugation in high-speed railways, based on vehicle-track coupling dynamics theory and field-measured rail corrugation data of a certain line, a vehicle-track rigid flexible coupling dynamics model of the CRH3 high-speed train was established. Rail corrugation superimposed on the Wuhan–Guangzhou track irregularity spectrum was used as the excitation input.. The influence of rail corrugation with different wavelengths and depths on wheel-rail dynamics characteristics was explored, and the coupling effect of wheel concave wear at different operation mileages with rail corrugation on vehicle subsystem’s vibration acceleration was analyzed. Then, the safety threshold for rail corrugation was proposed. Results show that rail corrugation significantly increases the wheel-rail vertical force. The coupled action of wheel concave wear and rail corrugation further amplifies the wheel-rail force response, with more severe concave wear leading to greater wheel-rail forces. Worn wheels that have traveled 200 000 km show an 10.8% increase in wheel-rail vertical force compared to new wheels. By considering wheel concave wear, for high-speed trains operating at an average speed of 300 km/h, the recommended safety thresholds for rail corrugation of 50 mm, 80 mm, 100 mm, 120 mm, and 150 mm are wave depths of 0.023 mm, 0.036 mm, 0.05 mm, 0.054 mm, and 0.069 mm, respectively. In practical applications, maintenance strategies should be adjusted according to specific operating conditions and track structures, and grinding should be carried out promptly when wave depths exceed the safety thresholds.

Influence of Strong Seepage Fractures on Stability of Water-Resistant Rock Mass at Tunnel Face
PENG Tao, QIAN Rui, YUAN Feiyun, REN Dongxing, FU Wenxi, YE Fei
, Available online  , doi: 10.3969/j.issn.0258-2724.20240625
Abstract:

Water and mud inrush is one of the common geological disasters in karst tunnels, which poses a severe threat to the safety of tunnel construction and long-term operation. Therefore, ensuring the safety of water-resistant rock mass at the face is the key to the prevention and control of water and mud inrush disasters. To ensure the stability of water-resistant rock mass at the face under the action of strong seepage, a mechanical analysis model of water-resistant rock mass in strong seepage fractures at the tunnel face was constructed. The drag force effect generated by the water flow on the wall of the penetrating fracture was fully considered, and the minimum safe thickness formula of the water-resistant rock mass in the karst tunnel under the action of strong seepage was systematically derived. Parameter sensitivity analysis and engineering examples were employed; the stability of the surrounding rock during tunnel excavation was further investigated by using the finite-difference software FLAC3D, considering penetrating fractures and water-rich karst cavities ahead of the tunnel face; the influence mechanism of key parameters on the stability of the water-resistant rock mass was discussed in depth. The results show that the increase of fracture opening significantly enhances the effect of the drag force. Under the condition of strong seepage, the fracture drag force has a significant adverse effect on the minimum safe thickness of the water-resistant rock mass, and its influence degree increases nonlinearly with the increase of fracture opening. In addition, the influence of fracture at different location distributions on the minimum safe thickness shows significant differences. It is found that there is a critical transition point (height ratio of cantilever beams Ⅰ and Ⅱ = 5.07∶4.93) at which the upper and lower beams and slabs of the fracture are destroyed synchronously, which further reveals the internal relationship between the fracture location and structural instability. The findings of this study can provide a theoretical basis and engineering guidance for designing the thickness of water-resistant rock mass in karst tunnel face, mitigating karst-related disasters (e.g., grouting sealing, fracture reinforcement, and seepage control), selecting monitoring and early-warning indicators, and conducting risk classification during construction.

Research on Long-Wave Irregularity of Metro Rail and Wheel-Rail Short-Wave Roughness Spectra
WEI Kai, WANG Sen, MA Meng, DU Siyu
, Available online  , doi: 10.3969/j.issn.0258-2724.20250026
Abstract:

To establish a wheel-rail broadband excitation model for numerical simulation and prediction of metro vibration noise, a study was conducted on the selection for long-wave irregularity of metro rail and wheel-rail short-wave roughness spectra based on measured data of long-wave irregularity of metro rail and rail surface roughness from Beijing metro lines, as well as wheel tread out-of-roundness (OOR) from Qingdao metro lines. Firstly, the Welch method was applied to the measured data of rail irregularity to estimate the power spectral density, and the expressions for the short-wave irregularity power spectrum of rail and the OOR power spectrum of wheel, as well as their fitting parameters were provided. In terms of long wavelength, typical rail spectra were selected according to the principle of similarity in individual high and low TQI values corresponding to speed classes. Regarding short wavelength, the short-wave irregularity power spectrum of the rail surface and the OOR power spectrum of the wheel were respectively divided into five grades based on quartiles, and a wheel-rail roughness spectrum was proposed by considering their coherence. The results show that the selected American sixth grade spectrum aligns well with the long-wave (1–42 m) irregularity power spectrum of metro lines with 80 km/h designed by Beijing, according to the principle of similarity in individual high and low TQI values. Both the short-wave irregularity power spectrum of the rail surface and the wheel OOR power spectrum exhibit a significantly skewed amplitude distribution under the same short wavelength. The wheel-rail roughness spectra of different levels in the short-wave range of 0.01–1.00 m are mainly dominated by the wheel OOR spectrum, and in severe cases, it even exceeds the short-wave spectrum recommended by the China Academy of Railway Sciences.

Study on Complex Mechanical Behavior of Coarse Granular Materials Based on Continuous-Discontinuous Deformation Analysis
XU Dongdong, LU Bo, WU Aiqing, ZHU Jiebing, WANG Bin
, Available online  , doi: 10.3969/j.issn.0258-2724.20240448
Abstract:

To investigate the mechanical meso-mechanism and macroscopic response of coarse granular materials in rockfill dams, and to overcome the limitations of traditional continuum-based methods in simulating force chain evolution and particle breakage, a continuous-discontinuous deformation analysis method suitable for triaxial tests on coarse granular materials was developed. This method, based on the traditional discontinuous deformation analysis (DDA) framework, introduced a hybrid displacement mode to differentiate the mechanical responses of various blocks. Critical damping was employed to accelerate computational convergence, and a continuous-discontinuous simulation technique was proposed to characterize particle breakage. Through conventional triaxial numerical simulations, the deformation evolution, force chain development, particle breakage, and shear band formation during loading were analyzed, with particular emphasis on the effects of size and end friction. The results indicate that the simulation outcomes agree well with experimental data, reflecting both the macroscopic mechanical response and meso-mechanisms of coarse granular materials. The size effect leads to a 21.3% increase in peak stress under 0.4 MPa confining pressure, whereas its influence becomes negligible at 3.0 MPa. Under 3.0 MPa confining pressure, end friction contributes to an approximately 7.4% increase in peak stress. This study provides an effective numerical tool for further understanding the mechanical properties of coarse granular materials.

Design and Test Study of New Dual-Curved Carbon Fiber-Reinforced Polymer Plate Anchorage
HOU Suwei, TANG Zhiqian, TIAN Shuheng
, Available online  , doi: 10.3969/j.issn.0258-2724.20250086
Abstract:

To address the issues of uneven transverse distribution of compressive stress in existing carbon fiber-reinforced polymer (CFRP) plate anchorage technology, which leads to tearing failure of plates under tension and the difficulty in controlling bolt preload, a dual-curved CFRP plate anchorage was developed. This device’s main construction features were outer clamping plates with transverse arc surfaces to ensure uniform transverse distribution of compressive stress sustained by the CFRP plate and inner clamping plates with longitudinal arc surfaces to fully compress the middle section of the anchorage zone and avoid shear failure at the ends. The compressive displacement of the outer clamping plates was controlled by the limit plate to anchor the CFRP plate with the required preload. The anchorage performance test for six groups of new anchorage specimens was carried out to explore their stress-bearing principle. Finite element software ANSYS was used to perform stress simulation; the key parameters affecting anchorage performance were deeply analyzed; the stress-bearing performance of CFRP plate anchorage was compared with that of flat-plate anchorage. The results indicate that the thickness of the outer clamping plates is the dominant factor influencing the non-uniform transverse distribution of compressive stress sustained by CFRP plates. When the thickness of outer clamping plates is 27 mm, the CFRP plates exhibit a higher uniform transverse distribution of compressive stress, with a compressive stress difference of only 9.3 MPa between maximum and minimum values. The curved surface design of inner and outer clamping plates effectively optimizes the compressive stress distribution state of the CFRP plate, while minimizing stress difference among carbon fibers and resulting in significantly improved synchronization of transverse deformation of CFRP plates. Although the specimens ultimately fail in a tearing state, the new anchorage achieves an exceptional anchoring efficiency coefficient of 104.16%.

Influence Analysis of Traction Rods on Locomotive Axle Load Transfer
CHEN Qinghua, WANG Jianyi, RAN Xiangrui, GONG Jingchun, WANG Kaiyun
, Available online  , doi: 10.3969/j.issn.0258-2724.20240501
Abstract:

Mitigating axle load transfer by optimizing structural parameters of locomotives is a primary strategy for enhancing adhesion utilization. To address the unclear mechanism by which the traction rod of a locomotive affects axle load transfer, a theoretical calculation model for axle load transfer given the rotation of the traction rod was established based on quasi-static equilibrium. The influence of the traction rod parameters on axle load transfer was explored based on the Sobol sensitivity analysis. Further analysis investigated the influence of traction force magnitude, stiffness of rubber sleeves, and position of the traction rod on axle load transfer. The results show that when the rotation of the traction rod is considered, the calculation results of the theoretical model are closer to those of the Simpack model, and the theoretical model exhibits significantly higher computational efficiency than the dynamics model. The initial tilt angle of the traction rod has a greater effect on axle load transfer than other traction rod parameters, and the rotation of the traction rod causes the load on each axle to vary nonlinearly with the traction force. As the radial stiffness of rubber sleeves of the traction rod exceeds 160 MN/m, the change in the locomotive’s axle load transfer tends to be gentle. As the deflection stiffness of rubber sleeves of the traction rod increases from 20 N m/rad to 500 N m/rad, the locomotive’s axle load transfer increases by 25.7%. As the height between the end of the traction rod frame and the rail surface increases from 0.05 m to 0.80 m, the locomotive’s axle load transfer increases by 84.3%. As the longitudinal distance between the end of the traction rod frame and the frame centroid increases from 0.5 m to 3.5 m, the locomotive’s axle load transfer decreases by 30.4%. The rotation angle of the traction rod approaches 0 when the initial tilt angle of the traction rod is around 11°. The axle load transfer of the second and third wheelsets approaches 0 when the initial tilt angle of the traction rod is between 13° and 14°.

Study and Evaluation of Performance of Polymer Hybrid Fiber- Reinforced Cementitious Composites
LI Fuhai, LI Zhou, ZHANG Yuxuan, HUO Jiateng, DING Yijun, CHEN Zhao
, Available online  , doi: 10.3969/j.issn.0258-2724.20250089
Abstract:

To address the issues such as the high density and corrosion susceptibility of traditional steel-polymer fiber hybrid systems, the hybrid effects of three types of polymer fibers, namely polyvinyl alcohol (PVA) fibers, polypropylene (PP) fibers, and polyoxymethylene (POM) fibers, on the workability and mechanical properties of fiber-reinforced cementitious composites (FRCC) were investigated. A total of 12 groups of specimens were designed, including single fiber and different combinations of hybrid fibers, including PVA and POM, PVA and PP, and POM and PP. Flowability, compressive, bending, uniaxial tensile, and four-point bending tests were conducted. The performance of each group was systematically analyzed and comprehensively evaluated, while the hybrid fiber mechanism was also analyzed using a high-resolution optical microscope. The results have shown that hybrid fibers can effectively compensate for the shortcomings of single fibers, significantly enhancing the overall performance of FRCC. The flowability of A0.5M1.5P0 increases by 18.62% compared to the single PVA fiber. The bending strength of A1M1P0 is 14.53 MPa, a 24.94% improvement over the single POM fiber. The compressive strength of A0M1.5P0.5 is 73.33 MPa, a 25.50% increase over the single PP fiber. A1.5M0.5P0 exhibits the best tensile performance, with a tensile strength of 3.36 MPa and strain energy density of 54.2 kJ/m3, showing increases of 54.13% and 2158.33%, respectively, compared to the single POM fiber. The bending strength of A1M1P0 is 9.03 MPa, a 51.26% improvement over the single PVA fiber, while the equivalent bending strength of 3.89 MPa shows a 77.63% increase over the single POM fiber. Based on the comprehensive evaluation of the radar chart, A1.5M0.5P0 exhibits the best overall performance. Two fibers with different physical properties produce a synergistic effect when hybridized. For example, PVA fibers limit crack propagation, while POM fibers dissipate energy through deformation during crack propagation. The synergy of PVA and POM fibers enhances the toughness and load-bearing capacity of the composite material.

Mechanical Characteristics of Large-Diameter Shield Tunnel Segment Structure under High Water Pressure
SU Ang, FENG Kun, HE Chuan, LIANG Kun, GUO Wenqi
, Available online  , doi: 10.3969/j.issn.0258-2724.20240419
Abstract:

To investigate the mechanical characteristics of large-diameter shield tunnel segment structures under ultimate loads, a multifunctional shield tunnel structure loading device was used to conduct a prototype loading failure test on the shield tunnel segment structure of the Sutong gas-insulated metal-enclosed transmission line (GIL) utility tunnel project. The mechanical response laws of deformation, internal force, joint deformation, bolt and steel bar strain of the segment structure under ultimate load were analyzed. Based on the crack morphology, the failure characteristics and mechanisms of the segment structure were revealed. The results show that the deformation of the segment structure presents a horizontal duck egg shape, the bending moment presents a butterfly shape, and the axial force presents a circular shape. The maximum single point change rate is 10.3‰, the maximum positive bending moment is 3896 kN•m, and the maximum axial force is 17099.28 kN. The maximum bolt strain is 1871 με, far less than the yield strain, and the maximum steel bar strain is 2213 με, exceeding its yield strain. The safety margin for the deformation indicator of the segment structure is between 0.92 and 1.72, while for the strength indicator, it is between 0.10 and 0.16. The strength indicator reaches the ultimate failure value before the deformation indicator. Reinforcement yielding or the main crack penetration is recommended as the judgment criterion of the ultimate limit state for the normal serviceability of the segment structure. Cracks of the segment structure first appear in the middle of the segment B3. When the load is 1778.4 kN, multiple cracks penetrate through the inner and outer arc surfaces of the arch crown and arch bottom, with the maximum widths of the cracks on the inner and outer arc surfaces reaching 4.5 mm and 11.5 mm, respectively. The steel bar of the segment B3 reaches the yield condition, and the segment structure is damaged due to local instability.

Seismic Design Method of Ultra-High Performance Concrete Prestressed Connection of Prefabricated Bridge Pier
ZHU Shengchun, CUI Bing, SONG Yingtong, WANG Kangkang, HAN Yongan, ZHAO Canhui
, Available online  , doi: 10.3969/j.issn.0258-2724.20250103
Abstract:

To enhance the industrialization level of the prefabricated bridge pier and develop connection forms compatible with the performance of the ultra-high performance concrete (UHPC) pier, an external prestressed connection joint and its tensioning process were proposed. The feasibility of this configuration and technique was validated through full-scale tests. Combined with finite element analysis, the influence of joint prestress and axial compression ratio on the performance of the UHPC pier was investigated, and a design method for prestressed joints based on the capacity design method was proposed. The research results have shown that the prefabricated pier with the external prestressed connection joint experiences typical bending failure, with the UHPC concrete of the pier body being crushed, but no interfacial separation occurs at the interface between the pier and the cap. The stress variation of the prestressed high-strength reinforcement is basically proportional to the horizontal load of the pier, with the maximum variation amplitude of only 9%, indicating that the connection performance of the external prestressed connection joint used in the specimen has reliable connection performance and good overall structural performance. When the tension force of the high-strength reinforcement in the pier is small, joint opening occurs at the interface between the pier and the cap, reducing the stiffness of the prefabricated bridge pier, but it has little effect on the peak load-bearing capacity of the specimen.

Intelligent Reflecting Surface-Assisted and Artificial Noise Enhancement-Based Beamforming Method in Covert High-Speed Rail Communications
LI Cuiran, SUN Shujing, ZHANG Zepeng, WANG Huiqin, XIE Jianli
, Available online  , doi: 10.3969/j.issn.0258-2724.20240424
Abstract:

To address the prevalent issues of low effective throughput and limited covertness in high-speed rail (HSR) wireless communication systems, an optimization problem was formulated to maximize the system’s effective throughput, subject to constraints on the covert requirement, transmit power of the maximum artificial noise (AN), and the unit modulus of the intelligent reflecting surface (IRS) phase shifts, and a beamforming method for covert HSR wireless communications based on IRS assistance and AN enhancement was designed. An alternating optimization strategy was adopted, decomposing the coupled optimization variables into three subproblems, including base station beamforming, IRS phase shift optimization, and AN transmit power optimization. The covert requirement constraint was mapped onto a complex circle manifold using the quadratic transform method from fractional programming. The conjugate gradient (CG) algorithm was employed to optimize the IRS phase shifts. The Dinkelbach algorithm was used to design the AN transmit power, and these steps were iterated alternately. Simulation results demonstrate that the proposed algorithm achieves a lower computational complexity. Under high-speed scenarios, it enhances the system’s effective throughput by 27.31% and improves the covert transmission performance, which is important for enhancing the security of information transmission in HSR wireless communication systems.

Study on Tool Condition-Integrated Online Optimization of Process Parameters
SUN Yi, GAO Hongli, SONG Hongliang, YOU Zhichao
, Available online  , doi: 10.3969/j.issn.0258-2724.20240578
Abstract:

As demands for manufacturing quality and production efficiency continue to rise in modern industry, tool wear has emerged as a critical constraint affecting surface roughness. Traditional tool condition monitoring and process parameter optimization methods are often based on empirical models or static optimization strategies, limiting their adaptability to complex, dynamically changing, multivariable environments. In response, this study proposes an innovative approach integrating multi-scale distribution ratio (MSDR) with Bayesian multi-armed bandit (BMAB) for process parameter online optimization, incorporating real-time tool condition data into the optimization framework. Additionally, by combining Bayesian optimization and multi-armed bandit strategies, this method enables real-time adjustments to process parameters in dynamic manufacturing environments, effectively balancing exploration and exploitation to maximize machining efficiency. Compared to mainstream methods, MSDR demonstrates exceptional precision and stability in tool condition monitoring, achieving MAE, SMER, and RMSE values of 0.145, 0.258, and 0.194, respectively. BMAB also performs exceptionally in optimizing cutting efficiency and computational effectiveness, achieving 2305 mm3/min and a runtime of 2.92 seconds, respectively. Therefore, tool state-aware online optimization of process parameters presents a novel and promising technical pathway for high-precision manufacturing.

Research and Application of Electrical Twin with Acoustic Metastructures for Vehicle NVH
DENG Shiqi, WU Yudong, HU Hao, QI Zhicheng, DING Weiping
, Available online  , doi: 10.3969/j.issn.0258-2724.20240473
Abstract:

Due to their exceptional wave manipulation characteristics, acoustic metastructures have attracted substantial attention in vehicle noise, vibration, and harshness (NVH). However, the further development and application of passive acoustic metastructures are limited by the narrow and non-tunable low-frequency bandgaps and bandwidth. To address this challenge, an electrically controlled acoustic metastructure was proposed to enable flexible bandgap tuning, and the corresponding electrical twin theory was established. According to the classical electromechanical analogy, a two-dimensional electrical twin circuit for a Kirchhoff-Love thin plate was established by the finite difference method. Then, an inductance-capacitance-resistance (LCR) resonant circuit was connected in series to form the twin circuit of a metastructure unit, with a tunable capacitor introduced to achieve electrical bandgap tuning. Finally, a spiral-shaped electrically controlled metastructure was derived from the twin circuit and verified through simulations and experiments. The results confirm that the twin circuit constitutes an exact electrical-domain mapping of the metastructure. The equivalent stiffness of the metastructure can be adjusted by the electrical control, thereby facilitating bandgap tuning. The resulting tuning law can be efficiently predicted and analyzed by the twin circuit. The designed spiral-shaped electrically controlled metastructure exhibits a significant order-tracking noise reduction effect for electric seats, with an average sound pressure level reduction of approximately 7.4 dB(A) in the wide frequency range of 200–460 Hz. The proposed twin circuit contributes to the electromechanical integrated design of electrically controlled metastructures, and it provides a theoretical paradigm for the study of different types of electrically controlled metastructures.

Metamodel-Driven Flexible Job Shop Embodied Agent and Its Scheduling System Construction
HU Mingzhu, ZHANG Weiwei, ZHANG Jian, ZHANG Haizhu
, Available online  , doi: 10.3969/j.issn.0258-2724.20240562
Abstract:
Objective

Flexible job shop scheduling optimization is an important research topic in digital manufacturing science, especially the random occurrence of abnormal disturbances such as machine failures and order changes, which disrupt the original production plan, causing problems such as unreasonable resource allocation, delayed order delivery, and increased production costs. In recent years, distributed multi-agent scheduling methods have been considered to be one of the most effective ways to improve the response speed of manufacturing system disturbances and reduce the negative impact of uncertain disturbances in the production process. In the context of job shop scheduling, designing an embodied scheduling agent that integrates the dynamic behavior of physical entities enables real-time environmental perception and autonomous decision-making based on behavioral feedback during disturbances. This ensures the efficient operation of the production system.

Method

To develop a multi-agent scheduling method for flexible job shops based on embodied agents and to enhance the stability and responsiveness of the system during disruptive events, a metamodel-driven approach for constructing embodied agents in flexible job shops was proposed. By instantiating this model, a scheduling system with a unified structure of embodied agents was achieved. First, to enable agents to perform autonomous decision-making and real-time dynamic adjustments, the concept of embodied intelligence was applied. Based on the existing interaction layer, decision layer, and adaptation layer of job shop scheduling agents, their instruction sets, behavior spaces, and sensory signal sets were associated and encapsulated to form embodied agents with physical bodies and behavioral spaces. Based on the resource composition characteristics of flexible job shops, the elements, relationships, and attributes of embodied agents were analyzed and abstracted. A metamodel for embodied agents in flexible job shops was proposed, enabling the unified modeling of embodied scheduling agents and providing a foundational model for their collaborative scheduling. Second, through instantiation operations such as inheritance, composition, aggregation, dependency, and association applied to the metamodel, a distributed multi-agent scheduling system with a unified structure and self-organizing collaborative operation capability was developed. Finally, a set of distributed multi-agent scheduling strategies was designed based on the different functions of different agents and the different information they can obtain. By integrating these strategies with a Q-game negotiation mechanism, collaborative scheduling among multiple agents was realized, thereby improving the stability of the scheduling method and enhancing its responsiveness to disruptions. This scheduling system, based on embodied agents, enabled the adjustment of scheduling strategies at the individual level when disruptive events occur. This approach effectively reduced the number of information exchanges during the scheduling process, improving the stability of the multi-agent system and enhancing its scheduling optimization capabilities in the face of disruptions.

Result

To validate the advantages of the proposed embodied scheduling agent modeling method and the multi-agent scheduling system, two small-scale manufacturing workshops producing structural components were used as case studies. The proposed method was compared with existing approaches in three aspects: embodied agent modeling, collaborative operation, and scheduling optimization. Experimental results demonstrate that the embodied agent modeling method proposed in this paper ensures model structure consistency, guaranteeing that the model adheres to predefined specifications and rules, thereby providing a unified modeling foundation for the collaborative scheduling of multi-agents. In the embodied multi-agent scheduling system, each agent generates a complete set of feasible individual strategies after evaluating all possible actions. Negotiation and interaction among agents are conducted based on these strategy sets. The number of interactions remains independent of the number of actions selected by the agents, resulting in an average reduction of 60.4% in communication volume and a 32.78% average decrease in computational response time. In terms of scheduling optimization performance, agents enhance the diversity of scheduling strategies during the negotiation process by adjusting their individual scheduling strategies, thereby improving the system’s global optimization capability. Compared to existing methods, the proposed approach achieves an improvement of more than 22.6%.

Conclusion

In summary, the proposed multi-agent scheduling method based on embodied agents not only improves system responsiveness and stability in handling disruptive events but also significantly enhances global optimization performance. This provides an efficient and robust solution for scheduling optimization in flexible job shops and offers new research directions and technical support for intelligent scheduling and collaboration in future manufacturing systems. In the future, further research could focus on applying the proposed method to larger and more complex manufacturing scenarios, integrating real-time sensing technologies and advanced machine learning algorithms to further enhance the adaptability and intelligence of the scheduling system.

Study on Hydraulic Interconnected Damper and Dynamic Performance of High-Speed Electric Multiple Units
DUAN Liang, SHI Huailong, SONG Chunyuan, LIN Jiazhi, CHEN Longfei, ZHANG Yaoxun
, Available online  , doi: 10.3969/j.issn.0258-2724.20240511
Abstract:

To address the inadequate anti-roll stiffness of the axle box in-board bogies of high-speed electric multiple units (EMUs), a primary suspension configuration was proposed to replace the traditional hydraulic damper with hydraulic interconnected units. The configuration could increase the anti-roll stiffness without increasing the vertical stiffness. Firstly, the equilibrium equations of oil pressure, flow rate, and output force were derived. A nonlinear dynamic model of the vehicle system was established using SIMPACK, and a simulation model of the hydraulic interconnected units was created in MATLAB/Simulink to facilitate co-simulation of the vehicle-hydraulic interconnected unit coupling system. Subsequently, the accuracy of the simulation model was validated based on the quasi-static characteristic test of the hydraulic interconnected units and the dynamic tests of the roller rig of the entire vehicle. The simulation analysis was conducted to ascertain the impact of pivotal parameters associated with the hydraulic interconnected units on the roll angle of the car body, derailment coefficient, and riding index for various operation conditions of vehicles. Finally, the field dynamic tests were conducted to verify the improvement in the dynamic performance of vehicles during curve negotiation. The results have shown that the roll stiffness of the interconnection unit is significantly greater than that of the traditional hydraulic dampers. The roll angle of the car body can be reduced by more than 0.5°, which is conducive to narrowing the dynamic limit and ensuring overturning safety. The field test results demonstrate that the dynamic indexes of hydraulic interconnected units are comparable to those of traditional oil pressure dampers. It is viable to address the issue of inadequate anti-roll capability of axle box in-board bogies by adopting hydraulic interconnected units.

Track–Bridge Longitudinal Dynamic Interaction during High-Speed Train Braking Process
LI Qi, LAI Yuchen, ZHANG Di, SHI Long, LI Kebing
, Available online  , doi: 10.3969/j.issn.0258-2724.20240409
Abstract:

To investigate the influence of the dynamic braking process of high-speed trains on the longitudinal force at the pier tops of simply supported beam bridges, the rail-level braking force time-history curve was first calculated and obtained via multibody system dynamics simulation. Longitudinal resistance tests were then conducted on WJ-8 type low-resistance fasteners to reveal the law governing the influence of loading frequency and vertical load on the longitudinal resistance properties of the fasteners. Finally, a finite element model for longitudinal track–bridge interaction in multi-span simply-supported girder bridges was established. In this model, the wheel’s vertical and longitudinal forces were applied as moving concentrated loads on the rail, taking into account the uneven distribution of dynamic vertical force among the fasteners and their corresponding vertical-load-dependent longitudinal resistances. The influence of braking stop position and number of spans on the dynamic response of the track–bridge system was analyzed using the dynamic time-history method, and the results were compared with those from static analysis. The findings indicate that the longitudinal resistance of the fasteners is not significantly influenced by the loading frequency but is sensitive to the vertical load they carry. The longitudinal forces in the rail and at the pier tops are maximized when the train stops braking at the abutment of the final span. While these forces increase with the number of spans, they stabilize beyond eight spans. A discrepancy is observed between the dynamic and static analysis results for the maximum rail stress and displacement, yielding a dynamic amplification factor of approximately 1.05. Furthermore, the dynamic amplification factor is about 1.07 for the pier experiencing the greatest braking force, but can reach up to 1.93 for piers subjected to smaller forces.

Complex Nonlinear Behavior of Parabolic Two-Hinged Arches Subjected to a Midspan Concentrated Force
HU Changfu, ZHU Shunshun, LV Jiabiao
, Available online  , doi: 10.3969/j.issn.0258-2724.20240363
Abstract:

To investigate the complex nonlinear behavior of parabolic two-hinged arches subjected to a midspan concentrated force, a theoretical method was proposed to reveal its rule. Based on the nonlinear strain–displacement relationship of arches in the Cartesian right-angled coordinate system, nonlinear equilibrium differential equations of parabolic two-hinged arches subjected to a midspan concentrated force were derived, as well as the corresponding high-precision approximate analytical solutions of these nonlinear equations. The common rules of complex nonlinear behavior of parabolic two-hinged arches subjected to a midspan concentrated force were investigated by the limitation analysis of these high-precision approximate analytical solutions in discontinuous points: 1) If and only if the modified slenderness ratio is greater than or equal to the limit-pattern critical slenderness ratio, limit-pattern nonlinear behavior occurs in parabolic two-hinged arches subjected to a midspan concentrated force. Moreover, multiple extreme points appear on the limit-pattern nonlinear equilibrium path, and the number of extreme points is positively correlated with the parameter k. 2) When limit-pattern nonlinear behavior occurs in parabolic two-hinged arches subjected to a midspan concentrated force, the limit-pattern nonlinear equilibrium path passes through specific points. The coordinates of these points are fixed and do not change with variations in the modified slenderness ratio. 3) If and only if the modified slenderness ratio is greater than or equal to the bifurcation-pattern critical slenderness ratio, bifurcation-pattern nonlinear behavior occurs in parabolic two-hinged arches subjected to a midspan concentrated force. This bifurcation-pattern nonlinear behavior exhibits multiple equilibrium paths. Comparisons against nonlinear finite element results demonstrate that the proposed approximate analytical solutions of nonlinear equilibrium of parabolic two-hinged arches subjected to a vertical midspan concentrated force have sufficient accuracy, and the rules of complex nonlinear behavior of parabolic two-hinged arches subjected to a midspan concentrated force agree well with nonlinear finite element results. The maximum relative error is 9.05%, which meets the needs of engineering accuracy.

Uncertainty Quantification for Seismic Vulnerability of Bridge Based on Bootstrap Method
CHEN Zhiqiang, ZENG Yongping, CHEN Zhiwei, DING Zihao, ZHANG Jin
, Available online  , doi: 10.3969/j.issn.0258-2724.20240630
Abstract:

To investigate the influence of ground motion uncertainty on the seismic demand and vulnerability of bridge structures and to clarify the propagation of this uncertainty in seismic vulnerability analysis, a quantitative method based on the Bootstrap method was proposed for assessing the uncertainty in the seismic vulnerability of bridges. Firstly, the relationship between the ground motion intensity index and the seismic demand of bridge structures was determined through probabilistic seismic demand analysis. Secondly, by considering the influence of ground motion sample size on the seismic demand model and vulnerability of bridge structures, the Bootstrap method was used to simulate the uncertainties in both probabilistic seismic demand model parameters and vulnerability curves. Finally, by taking a three-span simply supported beam bridge as an example, seismic vulnerability analyses were conducted using 50, 100, and 300 seismic records to quantify the variability of probabilistic seismic demand models and vulnerability under different ground motion samples. The results indicate that the seismic demand and vulnerability of bridge structures are subject to significant uncertainties under the ground motion. When 100 seismic records are used, the variability in failure probability of the bridge under various damage states exceeds 10%, and that under severe damage states reaches up to 30%. In seismic vulnerability analysis of bridges, it is advisable to represent the failure probability of bridge structures under different ground motion intensities as interval random variables to account for variability in seismic vulnerability due to seismic record samples. The Bootstrap method can effectively simulate the uncertainty in the seismic demand and vulnerability of bridge structures, providing an effective approach for statistical uncertainty simulation and seismic vulnerability analysis of probabilistic seismic demand models of bridge structures under small sample sizes.

Joint Dispatch of Cross-Regional Emergency Supplies Considering Differential Disaster Severity
WU Jiani, MENG Zejia
, Available online  , doi: 10.3969/j.issn.0258-2724.20250136
Abstract:

To enhance the response efficiency of cross-regional emergency rescue under major natural disasters, considering differential disaster severity in affected areas, the optimization of cross-regional emergency supplies dispatching with combined transport was conducted. Firstly, a differentiated disaster classification strategy and a comprehensive evaluation system were proposed. The CRITIC-TOPSIS method was employed to determine the risk level of each region. Then, a bi-level programming model was developed, in which the upper level minimizes the total emergency response time and the lower level maximizes fairness. The upper level incorporated the Beetle Antennae Search to improve the Particle Swarm Optimization algorithm for finding solutions, thereby determining the shortest time and the volume of supplies transported from supply points to distribution centers. This provides basic data and time constraints for the lower level. The lower level uses the NSGA-III to solve the supplies allocation problem, where its results influence the distribution of supplies to affected areas in the upper-level model. This interdependence may lead to adjustments in the upper-level transportation scheme, further optimizing the overall objective. Finanly, taking 5•12 Wenchuan Earthquake as a case study for the simulation, the results indicate that, in terms of the total emergency response time, the scheme considering disaster severity classification is 2.53% shorter than that without considering disaster severity classification. Regarding fairness, the scheme under disaster severity classification shows a positive correlation between the satisfaction rate of emergency supplies demands and the disaster severity level at different disaster-affected points, thereby better reflecting the differentiated strategies based on disaster severity levels and the fairness of emergency supplies dispatch.

Self-Sensing Performance of Ultra-High Performance Fiber-Reinforced Concrete Under Cyclic Loading
ZHOU Zhongyi, ZHAO Hongsheng, LIU Yan, CHEN Jianwei, ZHANG Wenming, HONG Yani
, Available online  , doi: 10.3969/j.issn.0258-2724.20240649
Abstract:

To study the self-sensing performance of ultra-high performance concrete (UHPC) mixed with steel fibers and multi-walled carbon nanotubes (MWCNTs) under different cyclic stress amplitudes, experimental studies were conducted on UHPC specimens with a steel fiber volume content of 2% and varying MWCNT contents. The results show that the initial resistivity of UHPC increases first and then decreases with the increase in MWCNT content, and the addition of 0.15% MWCNTs improves the conductivity of UHPC. When the MWCNT content is 0.15%, the sample exhibits optimal repeatability, with a repeatability coefficient of 0.019, and the linearity change of alternating current (AC) resistance presents a strong linear relationship with stress, with a linearity of 0.97. The stress sensitivity and strain sensitivity of the samples UHPC0 and UHPC0.05 first increase and then decrease with the increase in stress, while the stress sensitivity and strain sensitivity of samples UHPC0.1 and UHPC0.15 show a gradually decreasing trend. The maximum strain sensitivity and stress sensitivity of UHPC0.15 are 71.6% and 0.16%/MPa under different cyclic stress amplitudes, both appearing at a stress of 10 MPa. When the content of MWCNTs is 0.15%, UHPC exhibits the best self-sensing performance.

Influence of Axial Compression Ratio on Hysteretic Properties of Steel Shell-Concrete Pylon
LIANG Huanwei, XU Chunrong, LIN Yu, WU Jianli, XIA Fuyou, YAN Pengfan, ZHAO Canhui
, Available online  , doi: 10.3969/j.issn.0258-2724.20240559
Abstract:

To investigate the influence of the axial compression ratio on the hysteretic properties of the steel shell-concrete composite pylon, based on the composite pylon structure without longitudinal rebars, three hysteretic specimens were designed with the axial compression ratio as the research parameter. Through testing, the hysteresis curves, failure characteristics, and strain development of each specimen were obtained, and the mechanical behavior under large eccentric failure was analyzed. A finite element model was then established using ABAQUS for further analysis, and the boundary failure conditions of the pylon section were determined. Then, calculation formulas for the axial compression and bending moment of the section under boundary failure were proposed, and the effect of the steel ratio and concrete strength on the axial compression ratio under boundary failure was discussed. The research results indicate that under large eccentric failure, the section stiffness, peak bearing capacity, and energy dissipation capacity increase with the axial compression ratio. When the axial compression ratio increases from 0.056 to 0.166, the stiffness and the flexural capacity of the specimen improve by 20%. The boundary failure condition of the composite pylon section is defined by the yielding of the tensile-side steel shell and crushing of the compressive-side concrete. Under boundary failure, the section achieves its highest flexural capacity and stiffness. The proposed calculation formulas provide an accurate assessment of the axial compression ratio and flexural capacity under boundary failure. Both an increase in steel ratio and concrete strength lead to a reduction in the axial compression ratio at boundary failure. The axial compression ratio under boundary failure in the composite pylon section falls within the range of 0.44–0.56, making it well-suited for long-span suspension bridge towers with higher axial compression ratios.

Compression Performance of Cold-Formed Steel T-Shaped Composite Edge Columns with Web Stiffeners
CHEN Ming, HU Yunlong, HU Fangqi, LI Hengkai, LI Bushuan
, Available online  , doi: 10.3969/j.issn.0258-2724.20240481
Abstract:

To investigate the compression performance of cold-formed thin-walled steel T-shaped composite edge columns with web stiffeners, axial and eccentric compression tests were conducted on eight groups of specimens. The influence of “V”-shaped longitudinal stiffening ribs on the failure modes and bearing capacity of the components were revealed through finite element model validation and parameter analysis, and an improved calculation method for bearing capacity was proposed. The results indicate that under axial compression, local buckling first appears in the web of the T-shaped composite edge column without stiffening ribs, ultimately leading to overall crushing failure. After adding “V”-shaped stiffening ribs, the stiffness of the single-limb C-shaped steel web is enhanced; the local buckling mode of the T-shaped composite edge column is improved, and the bearing capacity increases by approximately 15%. As the eccentricity increases, the failure modes of the specimens remain similar, and the ultimate bearing capacity shows a decreasing trend. The bearing capacities under axial and eccentric compression predicted by the effective width method are conservative. Both the finite element results and the test results are greater than the calculated results, with the average ratios being 1.238 and 1.143, respectively. After modification, the ratio of the results predicted by the effective width method to the simulated values ranges from 1.000 to 1.074, indicating high prediction accuracy.

Experimental Study on Seismic Performance of Concrete Frame Structures Reinforced with High-Strength Steel Bars
ZHAO Hua, YUAN Weiguang, WEI Chengjin, LENG Donghang, CHEN Peng
, Available online  , doi: 10.3969/j.issn.0258-2724.20250036
Abstract:

To realize the requirements of minor post-earthquake damage, rapid repair, and functional recovery of reinforced concrete frame structures, three 1/2 scaled concrete frames were designed. One was an ordinary reinforced concrete frame, and two were concrete frames reinforced with high-strength steel bars (HG bars). Quasi-static tests were carried out to study the failure modes of the frames under cyclic loading. The effects of beams and columns with HG bars on seismic performance indexes including hysteretic curves, skeleton curves, residual deformation, repairability, and self-centering ability, were discussed. The results show that the use of HG bars in beams and columns effectively improves the overall bearing capacity and deformation capacity of the frame. Compared with the ordinary reinforced concrete frame, specimens NHGS2.5A15 and HGHGS2.5A15 show good displacement hardening effects. Their ultimate bearing capacities increase by 23% and 57%, respectively, and the displacement corresponding to ultimate load increases by 50% and 60%, respectively. These specimens have smaller residual deformations and higher reparability, and they demonstrate good self-centering capability and repairability performance.

Experimental Study on Dynamic Strength of Subgrade Loess under Continuous and Intermittent Loads
ZHUANG Xinshan, YANG Duan, LI Tong, LI Xiaofei
, Available online  , doi: 10.3969/j.issn.0258-2724.20240462
Abstract:

The cyclic dynamic stress generated during train operation presents a significant challenge to the dynamic strength of subgrade fill materials. Existing research has mostly simulated train loads using continuous loading methods, which fails to fully reflect the intermittency of these loads. To investigate the differences in dynamic strength of loess subgrade under continuous and intermittent loading, a series of consolidated undrained tests under continuous and intermittent loading conditions were conducted using a GDS dynamic triaxial apparatus. The influences of confining pressure and dynamic stress amplitude on the dynamic strength of the soil were examined. The effects of different loading methods on the dynamic strength and strength parameters of the subgrade loess were compared. The experimental results indicate that the dynamic strength of the loess subgrade increases with higher confining pressure, but the growth rate diminishes gradually. Both dynamic cohesion (cd) and dynamic friction angle ($ {\varphi _{\text{d}}} $) decrease with the increase in the failure cycles (lg Nf), showing an overall linear relationship. Under intermittent loading, the soil exhibits a marked increase in cd and$ {\varphi _{\text{d}}} $compared to continuous loading, with cd increasing by 2.18%–5.09% and $ {\varphi _{\text{d}}} $ by 4.03%–13.78%. By normalizing the dynamic strength using the static triaxial shear strength, an empirical formula for the dynamic strength of the loess subgrade based on static strength is proposed, which provides a critical basis for assessing the stability of the subgrade under dynamic loads.

Research Status and Prospects of Computer Vision-Based Crack Detection of Concrete Structure
XIE Mingzhi, FAN Dingmeng, JIANG Zhipeng, DENG Fei, WANG Kun, HAN Chen, YANG Yongqing
, Available online  , doi: 10.3969/j.issn.0258-2724.20240115
Abstract:

As one of the important contents of health monitoring of concrete structure, crack detection reflects the stress and damage state of the structure, and the detection and evaluation is the core technology to ensure structure safety for service. The traditional detection methods have limited coverage in time and space and are greatly affected by environmental and altitude factors, so the detection efficiency and accuracy are relatively low. Additionally, they are dependent on subjective judgment, which is easy to cause missed detection and false detection. The detection method based on computer vision is equipped with digital imaging equipment for data acquisition, input, and image processing to automatically analyze and identify the concrete surface, which has the advantages of high efficiency, accuracy, and objectivity and is widely used in the field of intelligent crack detection of concrete structures. The principle, method, and application of concrete crack detection based on computer vision were described in detail from four aspects: image acquisition, image processing, recognition algorithm, and structure evaluation. Besides, the application of crack image acquisition equipment and various image preprocessing methods in digital imaging technology was reviewed comprehensively, and the advantages, disadvantages, and applicability of different recognition algorithms were analyzed. At the same time, the shortcomings of current research were summarized, and the challenges and problems faced by the application of computer vision technology for equipment intelligence and lightweight network were analyzed. Then the corresponding solutions were proposed. Prospects are also presented from the aspects of multi-source data fusion and utilization, lightweight intelligent equipment, digital imaging and crack mapping, and high-efficiency and real-time structure evaluation.

Optimization of Automobile Firewall Acoustic Package for Multi-level Goals
HUANG Haibo, ZHENG Zhiwei, ZHANG Siwen, WU Yudong, YANG Mingliang, DING Weiping
, Available online  , doi: 10.3969/j.issn.0258-2724.20211086
Abstract:

To study the influence of automotive acoustic package design parameters on its multi-performance objectives, firstly, the traditional DBNs (deep belief networks) method was modified, and the SVR-DBNs (support vector regression-deep belief networks) model was proposed to improve the accuracy of model mapping. Secondly, from the perspective of vehicle noise transfer relationship and hierarchical target decomposition, a multi-level target prediction and analysis method was proposed. Finally, the proposed method was applied to the multi-objective prediction and optimization analysis of the MTL (mean transmission loss), weight and cost of the acoustic package for a real vehicle.The results show that the accuracy of SVR-DBNs method for the MTL, weight and cost target prediction of the acoustic package is higher than 0.975, which is better than that of the traditional BPNN(back propagation neural network), SVR and DBNs models. The optimization results based on the SVR-DBNs model are appropriate to the measured results, the comprehensive relative error of the predicted and tested targets is 1.09% (the absolute values of the relative errors of MTL, weight and cost are 1.44%, 1.04% and 0.71%, respectively). Compared with the original status, the MTL, weight and cost of the acoustic package have increased by 5.51%, 9.01% and 4.40%, respectively.

Anti-Disturbance Performance of Maglev Rotor Using Model Assisted Extended State Observer
JIN Chaowu, CAO Yingqing, ZHOU Jin, YE Zhoucheng, XIN Yu
, Available online  , doi: 10.3969/j.issn.0258-2724.20220803
Abstract:

With the increase in sinusoidal disturbance frequency, the performance of extended state observers (ESOs) will decrease. In order to improve the disturbance suppression ability of the ESO in the maglev rotor system, firstly, the mathematical model of a one-degree-of-freedom (1-DOF) maglev bearing rotor system was built. Secondly, ESO was designed, and the reasons for its reduced disturbance suppression effects were analyzed. On this basis, a model assisted ESO (MESO) was proposed to improve the bandwidth configuration and enhance the disturbance suppression effects. Then, the stability of the active disturbance rejection controller based on MESO was analyzed in the frequency domain. The effectiveness of the proposed observer was finally verified through simulation and experiments. The research results indicate that an increase in bandwidth amplifies the impact of system noises and increases the control voltage of the system. As the disturbance frequency increases, the suppression effect of MESO on high-frequency sinusoidal disturbance will decrease, but it can still reduce the modal amplitude of the rotor. After applying fundamental harmonic disturbance of 10 Hz−2 mm and fundamental impulse disturbance of 1g to the rotor at a rotating frequency of 50 Hz respectively, the rotor displacement under MESO control is reduced by 16.3% and 22.6%, respectively compared with that under ESO control, and the control voltage is reduced by about 14%.

, Available online  
Abstract:
Analysis of Limit Support Pressure Due to Shield Tunnelling with a Shallow Overburden Under Seepage
CAO Liqiang, ZHANG Dingli, LI Xinyu, LI Ao, SUN Zhengyu
, Available online  
Abstract:
Based on the semi-contained water model, the permeability of shield-crossing soil and overburden layers when the shield passed through the permeable soil was comprehensively analysed. The analytical solution of the head distribution along the tunnelling direction in the shield-crossing soil layer was derived, and the analytical solution of the two-dimensional seepage field was extended to the corresponding three-dimensional approximate solution. The active and passive failure modes of shallow-buried soil under steady-state seepage were determined using numerical modelling, and a corresponding cylinder-arc-corner-shaped model was established. Subsequently, the formulas of the two-limit support pressures at the tunnel face were obtained by introducing the above-mentioned three-dimensional seepage solution, and the newly calculated results were compared with those of the existing model. The results derived from the developed model were closer to the numerical solution. In addition, disturbance of the seepage field in front of the tunnel face was limited to three times the tunnel diameter and the values of the active and passive limit support pressure increased linearly with increasing head difference. The shield diameter and head difference were found to be the two main factors affecting the active limit pressure. The overburden thickness and shield diameter were the two major factors affecting the passive limit pressure. During tunnelling, the support pressure should be as close as possible to the in-situ transverse earth pressure using the approach that separately calculated the values for soil and water and should slightly fluctuate in the vicinity (preferably above it). The fluctuation range should be determined according to the deformation control standard.