• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    A method for support stiffness failure identification in a steel spring floating slab track of urban railway: a case study in China

    2024-03-18 06:57:48JunyuanZHENGCaiyouZHAODuojiaSHIPingWANGJianWANGBolongJIANGXiSHENG
    關(guān)鍵詞:浮置轉(zhuǎn)角敏感性

    Junyuan ZHENG, Caiyou ZHAO, Duojia SHI, Ping WANG, Jian WANG, Bolong JIANG, Xi SHENG

    Research Article

    A method for support stiffness failure identification in a steel spring floating slab track of urban railway: a case study in China

    1Key Laboratory of High-speed Railway Engineering, Ministry of Education, Southwest Jiaotong University, Chengdu 610031, China2School of Civil Engineering, Southwest Jiaotong University, Chengdu 610031, China301 DYNA (Qingdao) Test & Data Co., Ltd., Qingdao 266000, China4Rail Transit Digital Construction and Measurement Technology National Engineering Research Center, China Railway Design Corporation, Tianjin 300308, China5Institute of Urban Smart Transportation & Safety Maintenance, Shenzhen University, Shenzhen 518060, China

    The extensive use of steel spring floating slab tracks has effectively addressed the challenge of alleviating the environmental vibrations induced by urban rail transit systems. However, under the combined action of train dynamic loads and complex environmental factors, problems, such as the fracture of steel spring vibration isolators and suspension vibrations induced by the uneven settlement of the base, often occur. The failure of isolator support stiffness is often hidden in its early stages and is challenging to identify by conventional detection methods. At the same time, it will aggravate the wheel?–?rail interaction, accelerate the deterioration of track structure, and even affect the driving safety. This study first establishes a detailed coupled train-floating slab track-foundation analytical model. Then the influence of the vibration isolator support stiffness failure on the dynamic indices of the floating slab track system response is analyzed. A set of defect identification methods that can detect the number of failed steel springs, severity of damage, and their location is proposed. Finally, an intelligent monitoring system for support stiffness of floating slab track is built by combining the density-based spatial clustering of applications with noise algorithm and statistical data analysis and is applied to a rail line in southern China. During a three-year monitoring campaign, a suspension failure and a fracture of a steel spring were each successfully detected and detailed failure information was obtained. Field investigation results were consistent with the damage identification results. After repair, the track structure dynamic response returned to the average pre-damage level and further deterioration had been arrested. The proposed damage identification methods and monitoring system provide an approach for intelligent identification of track structure support stiffness failures.

    Floating slab track; Support stiffness; Detailed analytical model; Failure identification; Monitoring system

    1 Introduction

    With the increase in train line length and train speeds, the environmental vibrations caused by urban rail transit systems are becoming increasingly serious. Vibration reduction measures for metro lines can be divided into three types according to their different implementation locations. At the vibration source, many measures, such as floating slab track bed (Tamagawa, 2021), damping fasteners (Cui et al., 2016), and elastic sleepers (Auersch, 2017), can be used to weaken the propagation of vibration energy at the wheel?–?rail contact location and thus reduce the vibration energy input to the soil. Along the propagation path, vibration isolation trenches (Sitharam et al., 2018) and rigid barriers can be used for vibration isolation. Measures such as foundation isolation systems (Talbot, 2016) can be introduced for vibration control. The control measures at the vibration source reduce the outward spreading of vibration energy from the source and are more widely used. The floating slab track bed has a better vibration reduction performance than other measures applied at the vibration source. In the floating slab track bed, instead of the traditional cast-in-place track bed and foundation, a track slab elastically supported by steel springs (Bashir and Akhtar, 2022) or rubber pads (Yuan et al., 2020) is adopted. The basic structure is shown in Fig. 1. The existence of elastic support makes the inertia of the floating slab absorb the dynamic train loads and reduce the vibration transmission from the track to the substructure (Wilson et al., 1983).

    Although steel spring floating slabs have the best vibration isolation effect, up to 30 dB (Zou et al., 2022), many studies have recently attempted to further improve their performance (Nelson et al., 2018; Wang et al., 2022). However, the steel springs not only damp vibrations but also carry all the static and dynamic loads of trains and the track. In addition to determining the vibration damping performance, their condition directly affects the running condition and even the running safety of trains. There is currently a lot of research on track support failure. Dersch et al. (2021) proposed reducing fastener failure through parameter optimization using a finite element method but did not conduct an experimental verification. Yu et al. (2021) studied the pull-out resistance of different materials of sleepers and proposed a prediction formula for pull-out resistance. Lam et al. (2012) studied the relationship between ballast damage and sleeper vibration, providing ideas for damage detection. Ramos et al. (2021) compared the settlement performance of ballast and slab railway track structures in the laboratory and obtained consistent numerical results.

    However, relatively little research has been conducted on the failure of steel springs in floating slab track. Xu et al. (2020) established a 2D coupled train-floating slab track model and analyzed the influence of the number of failed springs on soil vibrations in combination with a finite element model. However, their model did not consider the spatial distribution of failed springs. Zhao et al. (2019) established a 3D coupled train-floating slab-track dynamic model, studied the effect of spring failure location on the train track dynamic responses, and concluded that spring failures at the slab ends are more harmful to driving safety and environmental vibrations.

    Fig. 1 Comparison between traditional track bed and floating slab track bed

    The above studies mainly focused on the adverse effects of steel spring failures. These harmful effects can be avoided if the spring failure is detected early. Generally, structural damage detection methods analyze the dynamic responses of the structure and detect structural damage from changes in the natural frequencies (Dilena et al., 2015), modal shapes (Gomes and Giovani, 2022), or flexibility matrix (Wickramasinghe et al., 2020). However, these approaches have several shortcomings (Zhu et al., 2021). For example, damage usually only affects the local response of the structure, which requires proper arrangement of sensors to be detectable. Modal shape measurement requires high accuracy and is easily affected by noise and human error. The flexibility matrix method requires prior knowledge of the undamaged structure.

    In addition to these traditional monitoring methods, some novel detection methods have been studied and applied in rail transit. The identification of suspension has been extensively carried out in vehicles. Ankrah et al. (2020) used vehicle acceleration combined with a support vector machine to predict the fault state of suspension systems. Similarly, Hong et al. (2020) proposed a data-driven framework for identifying suspension system status through train vibration. Some scholars also use train vibration for track state evaluation (Cong et al., 2020), but separating the influence of suspension system and track state is a challenge which makes this method difficult for stiffness failure identification of steel springs. For the track structure, damage detection approaches frequently use visual, ultrasonic, or electromagnetic inspections for identifying rail damage (Loveday et al., 2020), fastener failure (Chandran et al., 2022), and track slab cracking (Li et al., 2020). However, these technologies require access to the inspected element surface, but the steel springs are wrapped by sleeves and are therefore inaccessible. The spring break indicator developed in recent years tries to address this challenge. If a steel spring is broken, an indicator light will switch on. However, the device has many shortcomings: each steel spring requires an individual indicator, the light will not turn on until a certain degree of fracture is reached, and the degree of suspension cannot be displayed. In addition, the light can only by detected by night patrols rather than monitored in real-time. Conventional vibration monitoring of track structure is easy to implement but signal anomaly detection is challenging. On the one hand, the conventional time-domain data change little, so it is difficult to identify failure whereas; on the other hand, the regular detection cycle is long and its time is short. Even if abnormal data are found, the failure may already be severe. If failure can be discovered early, the maintenance costs can be significantly reduced. Zhu et al. (2021) proposed a failure detection method based on a convolution neural network using the vibration data from a floating slab, but their approach required undamaged and damaged data sets, which is difficult to satisfy in the field.

    Fig. 2 Proposed technology structure. RMS: root mean square

    This study presents a more straightforward method to identify the number, degree, and location of steel spring failures. Firstly, a refined train-floating slab track-foundation analytical model is established and the sensitivity of various floating slab dynamic response indices to failure is analyzed. Methods for identifying the failure of steel springs are proposed based on the specific sensitive dynamic response indices. Based on the proposed method, an intelligent monitoring system for the support stiffness of the floating slab track (ISFST) system was developed and applied in the field, where it detected steel spring failures. An on-site inspection confirmed the failure identification results. Furthermore, the monitoring data returned to normal pre-damage values after repairs. The structure of the proposed technology is shown in Fig. 2. Although steel spring support failure is studied in this paper, the proposed method will be applicable to other types of supports as well.

    2 Detailed coupled train???floating slab track?foundation analytical model

    The developed coupled train?floating slab track?foundation dynamic model comprises three subsystems: the train, the track, and the foundation, as shown in Fig. 3.,, andrefer to the translational directions.,, andrefer to the rotational directions. The subscripts c, w, and t refer to the train body, wheel, and bogie, respectively. The subscripts s and p refer to the first and second suspensions, respectively. The train is a multi-rigid body structure, while the track and foundation are flexible bodies. The train consists of six A-type vehicles, each consisting of a train body, two bogies, four wheelsets (Chinese LMA tread), and eight axle boxes. The train body and bogies have three translational and three rotational degrees of freedom (DOFs), while the wheelsets have two translational DOFs (the translation in the-direction is realized through the specified train driving speed) and three rotational DOFs. The wheels and axle box are connected through primary springs and dashpots, and the bogies and train body are connected through secondary springs and dashpots. The train parameters are listed in Table 1. The track and foundation model comprise, from top to bottom, the rails, fasteners, floating slabs, steel springs, substrates, and soil. A Timoshenko beam was used to simulate the rail and the finite element method was adopted to model the floating slabs and foundation. The rail was connected to the floating slabs through the fasteners (6-DOF spring-dashpot elements), the floating slabs were connected to the foundation through steel springs (3-DOF spring-dashpot elements), and eventually the foundation was connected to the soil. The longitudinal spacing of fasteners was 0.6 m and the transverse spacing was 1.5 m. The longitudinal spacing of steel springs was 1.6 m, while the transverse spacing was 1.5 m. The slabs were connected by shear hinges (simulated by 3-DOF springs with a large stiffness) to ensure the consistency of slab end displacements.

    Fig. 3 Train structure

    Table 1 A-type train parameters

    Table 2 Track parameters

    3 Dynamic response sensitivity to damage

    Numerous physical quantities can describe the vertical motions of the floating slab, such as linear acceleration, velocity, displacement, angular displacement, angular velocity, and angular acceleration. Selecting physical quantities sensitive to steel spring failure can improve the damage identification accuracy. A failed spring was considered to have two states: complete failure, i.e., the steel spring is fully fractured, and partial failure, i.e., the steel spring is suspended. A partially failed, or suspended, steel spring may still provide some support to the floating slab. Both damage states are collectively referred to as failure in the following text. Slab No. 17 was the central one of all the floating slabs in the model and was the focus of analyses. Different slab No. 17 spring failure cases were considered and compared. All failures in this section are fracture for obvious comparison. The steel springs were numbered as shown in Fig. 4.

    Fig. 4 Steel spring failure cases

    The dynamic response results are shown in Fig. 5. In Figs. 5a and 5b, numerous signal bursts in acceleration and velocity can be seen, caused mainly by local impact, resulting in significant fluctuations in peak values. However, the RMS values are more stable. It can be seen that the fracture of a spring causes an increase in the RMS value, and the sensitivity of velocity to fracture is higher than that of acceleration. The displacement in Fig. 5c is smoother relative to acceleration and velocity, and the changes in peak and RMS values are more sensitive, evident, and stable. Figs. 5d and 5e show angular displacements with respect to the- and-axes, respectively. Positive angular displacement satisfies the right-hand screw rule and is not sensitive to the fracture of the steel springs. However, in Fig. 5d, the angular displacement corresponding to one spring failure is the largest, and the fracture position is located in the positive-direction, which may suggest a technique for spring fracture location identification. As in translational motion, the angular displacement stability and sensitivity were higher than those of the angular velocity and acceleration, which are not shown here.

    Fig. 5f shows the 1st-order vertical mode natural frequency of the floating slab. This vibration mode is the overall vertical resonance of the floating slab?–steel spring system. The theoretical natural frequency was obtained by finite element modal analysis. The experimental frequency identified in the figure was obtained through operational modal analysis (OMA) and the data source came from Fig. 5a. This method can obtain the natural frequency from the dynamic response of the slab under external load excitation, avoiding separate modal tests such as hammer tests (Rajaram and Nelson, 2019). The modal identification was divided into two steps. Firstly, the free attenuation signal was extracted from the vertical acceleration time-domain signal, and then the modal parameters were calculated using a time-domain modal parameter identification method. Fig. 6 shows the main steps of modal parameter identification and the identified 1st-order natural frequency without failure. The random decrement technique (RDT) (Vandiver et al., 1982) only uses the center point data, while the natural excitation technique (NExT) (James III et al., 1995) uses the center point and corner point data. It can be seen that the differences in data handling by the different methods are the main cause of the difference in frequency identification. RDT is not suitable because the signal has no apparent attenuation, which is mainly caused by poor stability of train excitation. However, NExT uses more data and the identification results are very good. Autoregressive moving average model (ARMA) (Bertha and Golinval, 2017) and stochastic subspace identification (SSI) are more accurate for the preprocessed data. However, ARMA has fewer calculation steps and higher efficiency and is more suitable for processing data from many train runs. As shown in Fig. 5f, the difference between the calculation results using NExT/ARMA and the theoretical results is within 0.3 Hz. The theoretical and experimental identification results show that the 1st-order vertical mode natural frequency of the floating slab is sensitive to spring fracture.

    Fig. 5 Sensitivity analysis of dynamic response indices of slab No. 17: (a) vertical acceleration; (b) vertical velocity; (c) vertical displacement; (d) X-rotation; (e) Y-rotation; (f) 1st-order vertical mode natural frequency. References to color refer to the online version of this figure

    4 Steel spring failure identification

    4.1 Failure number identification

    The 1st-order vertical mode natural frequency at different failure numbers, gaps, and positions is investigated in this section. The failure gap here refers to the vertical distance between the spring contact point and the floating slab after the system reaches static equilibrium under gravity. Different spring failure gaps for slab No. 17 were simulated by a nonlinear spring. Four failure gaps of 0, 2, 4, and 6 mm were considered, where 0 mm indicated that the failed spring was just in contact with the floating slab, and 6 mm meant the spring was broken because it could not contact the slab at all (Fig. 5c). The force–displacement relationships for steel springs with different failure gaps are illustrated in Fig. 7 and Table 3. The failure causes nonlinear behavior of the floating slab and it is challenging to solve the natural frequency by the finite element method. The modal identification method in the previous section can avoid this problem. This paper only studies the failure of one or two springs. Four failure cases are considered: spring No. 1 failed, spring No. 2 failed, springs Nos. 1 and 4 failed, and springs Nos. 2 and 4 failed. In the two-spring failure cases, e.g., springs Nos. 1 and 4 failed by 1 mm and 3 mm, respectively, the spring positions were noted as positions 14, and the gaps were abbreviated as 1+3.

    Fig. 6 Steps and results of natural frequency identification. ITD: Ibrahim time domain technique; STD: sparse time domain technique; QR: an orthonormal matrix Q and an upper triangular matrix R

    Table 3 1st-order vertical mode natural frequencies of floating slab subjected to different combined failures

    The 1st-order slab vertical mode natural frequencies of different failure cases are shown in Table 3. Each 1st-order frequency dropped significantly when the steel spring just failed (0 mm). After the failure gap exceeded 2 mm, the frequency gradually tended to the spring fracture frequency, which hardly changed with the failure gap. In addition, when the number of failed springs is the same, different failure gaps and positions had little effect on the frequency. This is because the failed spring can only be compressed and not stretched, equivalent to changing from a bidirectional spring to a unidirectional spring, which plays a significant role in the change of frequency. These characteristics show that the 1st-order vertical mode natural frequency can be used robustly to identify the number of failed springs.

    Fig. 7 Force?–?displacement relationships for steel springs with different failure gaps

    4.2 Failure severity identification

    The failure gap was used to quantify failure severity. Since the frequency was not sensitive to the failure gap, the maximum displacement was considered as an index to quantify the failure gap. Table 4 shows the peak displacements at various failure gaps. For one-spring failures, failure position No. 1 or No. 2 was selected. Before complete failure (gap<5 mm), there was a stable linear relationship between the displacement peak and the failure gap. For two-spring failure cases, the symmetrical case of positions 16 and asymmetrical case of positions 15 were selected. It can be seen that the different positions did not affect the maximum displacement in the case of two-spring failures. In addition, except for the case of 2+4, when the total failure gap of two springs was equal to the total failure gap of one spring, the maximum displacement of both springs was the same. This was because, in the case of 2+4, the two springs failed only partially, but in the case of one spring, 6 mm signified a complete failure. This indicates that the maximum dynamic displacement of the floating slab was only related to the total failure gap. However, when the failure gap changed, the displacement change rate was low. The on-site results will be further complicated because there are inevitable fluctuations in the actual track structure parameters. It was difficult to identify the failure gap accurately from displacements except when the failure gap was large.

    Table 4 Maximum vertical dynamic displacements of floating slab subjected to different failures

    4.3 Failure location identification

    Based on the results shown in Fig. 5d, failure position identification from the rotation angle was attempted. Since the fluctuation of peak values was too large, the RMS of rotation angle was used. Two cases of one-spring failure, Nos. 1 and 2, were considered. Fig. 8 shows the RMS of each rotation direction in these two cases, where the purple solid line is the RMS of the rotation angle when there is no defect. Note that there are two types of RMS in Fig. 8. One is the RMS of the positive rotation part, and the other is the RMS of the negative rotation part. For example, 1: RMS+in Fig. 8a indicates the RMS of positive rotation angle around the-direction when spring No. 1 failed, while 2: RMS?-?represents the RMS of the negative rotation angle around the-direction when spring No. 2 failed. It can be seen from Figs. 8a and 8b that 1: RMS->1: RMS+and 1:RMS?-?>1:RMS+in the case of spring No. 1 failure. This difference increased with the increase in the failure gap. In the case of spring No. 2 failure, 2: RMS?-?≈2: RMS+and 2: RMS?-?>2: RMS+. The above observations can be summarized as follows: if a spring in a specific direction fails, the rotation in this direction will increase. For other failure cases, such as spring No. 6 failure, a similar analysis can be conducted, which showed that 6:RMS-<6:RMS+and 6:RMS-<6:RMS+. The equivalence here assumes that the rotation angle is symmetric when there is no failure, which is confirmed by the purple line in Fig. 8a.

    Failure of two springs can be analyzed in a similar way by division into six independent cases, as shown in Fig. 8c. Damage positions 1 and 2, 1 and 3, and 2 and 5 can be identified from the- and-rotations. The remaining positions require special consideration. For positions 1 and 6, since the failure positions are centrosymmetric, the rotation should be considered as diagonal. Diagonal rotation can be calculated from the- and-rotations. Note that the directions of diagonal lines for damage cases 16 and 43 are denoted asand, respectively (Fig. 8c). Considering the complete spring fracture at positions 16, the calculated RMS±is 5.3×10-5rad, and RMS±is 4.4×10-5rad, indicating that failure occurred at positions 1 and 6 instead of 3 and 4. For damage positions 1 and 4 and 1 and 5, their rotations in the- and-directions are similar and should be judged in combination with the diagonal rotation. In case of complete failure, the RMS?-?and RMS-of positions 14 are both 6.6×10-4rad, and the RMS?-?and RMS?-?of positions 15 are 5.4×10-4and 4.8×10-4rad. Therefore, it can be discerned whether the failure is in positions 14 or positions 15. Fig. 8d shows the rotations of the three exceptional cases.

    Fig. 8 Effective value of rotation angle for different failure positions: (a) RMS of Y-rotation; (b) RMS of X-rotation; (c) six cases for two-spring failures; (d) effective values of diagonal corner rotations in special cases. References to color refer to the online version of this figure

    (1) For one-spring failure,

    1: RMS->RMS+, RMS->RMS+;

    2: RMS->RMS+, RMS-≈RMS+.

    (2) For two-spring failure,

    12: RMS->RMS+, RMS->RMS+;

    13: RMS->RMS+, RMS-≈RMS+;

    14: RMS-≈RMS+, RMS->RMS+, RMS-≈RMS-;

    次日清晨,我早早兒就爬起。心想今朝一定要見到狼剩兒,我都跟大梁說了,今朝要趕回去的。抬眼瞄賽十娘,她困得正香。我躡手躡腳走到門口,伸手拉門,外面還鎖著。輕輕搖了兩下,冇得人應(yīng)。

    15: RMS-≈RMS+, RMS->RMS+, RMS->RMS-;

    25: RMS-≈RMS+, RMS-≈RMS+;

    16: RMS->RMS-.

    5 Floating slab monitoring system

    Based on the above analysis, the ISFST was proposed. The system has five parts: monitoring module, monitoring terminal, cloud server, data analysis system, and user terminal, as shown in Figs. 9 and 10. The monitoring module comprises sensors and communication units. The sensors are three-axial acceleration and three-axial angular velocity sensors. The communication unit includes a converter and a transmission cable. The converter changes analog electrical signals into digital ones. The transmission cable connects the monitoring module and the monitoring terminal. The monitoring terminal receives the digital signal from the monitoring module and transmits the data to the cloud server through a wireless network. The data analysis system processes the vibration data on the cloud server and then sends the data processing results to the user terminal. If the data is abnormal, the system automatically sends an alarm message to the user.

    The system was successfully implemented on a section of a metro line in southern China. The radius of this rail line section is 390 m, the dimensions of the floating slab are 3.576 m×2.600 m×0.340 m, the gap between slabs is 12 mm, and the natural design frequency is 12.5 Hz. A monitoring module was arranged at the center of each slab. The monitoring area contains 16 boards, with a total length of about 60 m. The tunnel power supply was used to operate the entire monitoring system. Slightly different from the above theoretical analysis, the vibration isolators at the adjacent slab ends support the two adjacent floating slabs simultaneously and their stiffness is about twice that of the steel springs in the middle of the slab. Although this section is located on a curve, this does not affect the above identification rules. The functions of the monitoring system include more than steel spring condition monitoring. The training data can be obtained by analyzing the responses of the floating slabs caused by different train runs. The rail status can be monitored by comparing the responses of floating slabs at different times caused by the same train type. However, this study only focuses on the steel spring support stiffness monitoring.

    Fig. 9 Architecture of ISFST

    Fig. 10 On-site monitoring system: (a) monitoring module; (b) monitoring terminal; (c) power supply; (d) signal line

    Fig. 11 shows the clustering results of the data from 20 train runs at the beginning of the system operation. There were two types of data. Type I is all the data of slab No. 4, but a few data of slab No. 4 are classified as type II (outliers), indicating that a steel spring in slab No. 4 was damaged. As shown in Fig. 11b, although the planar projections of two types of data overlapped, the method achieved satisfactory results using three types of data (Fig. 11a), which shows that the method can reduce the impact of single-type data fluctuations on failure identification. This method was still very effective for data clustering from one train run, which shows that this method can effectively find failures.

    The failure characteristics were identified by statistical analysis. Fig. 12 shows the average value of all data on a particular day at the beginning of system operation, and the shaded area represents the standard deviation. It can be seen from Fig. 12a that the frequency of slab No. 4 was much lower than that of other slabs, indicating that the support stiffness was insufficient, which may be due to the suspension of the support or the fracture of the steel spring. Compared to the design frequency, the frequency of slab No. 4 has decreased by 1.4 Hz. According to Fig. 5f, there is a high probability of one failure. The dynamic displacement of slab No. 4 shown in Fig. 12b was high, and the difference from other slabs was almost 0.5 mm. Based on these observations, it was speculated that lifting might have occurred. Figs. 12c and 12d show the difference between each slab RMSs of positive and negative rotations. It is evident that for slab No. 4, RMS+>RMS?-?and RMS+

    Fig. 11 Cluster analysis based on sensitive response indices: (a) space distribution; (b) projection in XYZ. References to color refer to the online version of this figure

    After a field inspection performed by maintenance personnel, it was found that the shared isolator of slabs Nos. 3 and 4 had failed. The failure gap was about 2–3 mm, and a steel ruler with a thickness of 1.5 mm could be inserted. The isolator was repaired by adding vibration isolation gaskets, and the natural frequency of slab No. 4 was monitored the next day. It was found that the frequency had returned to the normal value, as shown in Fig. 13.

    In the following three years of long-term monitoring, the system also identified a steel spring fracture damage and cluster analysis quickly found that the abnormality was located in slab No. 7. Statistical analysis showed that steel spring No. 2 in slab No. 7 was broken. Since the process was similar to the above suspension failure identification, only the 1st-order vertical mode natural vibration frequency and vertical dynamic displacement are shown here. Based on the frequency difference, it could be determined that the number of failures was 1. Sensitivity analysis showed that the maximum dynamic displacement was the main difference between steel spring fracture and suspension data. The dynamic displacement of slab No. 7 in Fig. 14b was almost 0.5 mm higher than that of other slabs, so the possibility of fracture was high. The identification results were confirmed on-site, as shown in Fig. 15. After the damaged steel spring was removed and replaced with a new one, the natural frequency returned to the pre-damage value.

    Fig. 12 Damage indices of floating slab: (a) 1st-order vertical mode natural frequency; (b) vertical dynamic displacement; (c) difference in effective values of rotation angles in X-direction; (d) difference in effective values of rotation angles in Y-direction

    Fig. 13 On-site confirmation of steel spring hoisting

    Fig. 14 Damage indices of floating slab: (a) 1st-order vertical mode natural frequency of floating slab; (b) vertical dynamic displacement

    Fig. 15 On-site confirmation of steel spring fracture

    6 Conclusions

    This study was concerned with the problem of steel spring failures in steel spring floating slab tracks. It established a refined coupled train-floating slab track-foundation analytical model and analyzed the sensitivity of several dynamic response indices of a floating slab to failure. Adopting the damage-sensitive physical quantities, a method for identifying the number, severity, and location of steel spring failures was proposed. Based on this method, the ISFST was established and installed on a rail line in southern China. Cases of steel spring suspension and fracture were successfully identified and detailed failure information was obtained.

    The sensitivity analysis of the dynamic model responses showed that the 1st-order vertical mode natural frequency, the maximum vertical dynamic displacement, and the RMS of the rotation of the floating slab under the train loads are easily affected by steel spring failures. The specific features are as follows:

    1. The experimental 1st-order vertical mode natural frequency of the floating slab under train load was obtained using the NExT/ARMA algorithms. This index is most sensitive to the number of failed springs but is insensitive to the spring failure gap or failure location. Thus, it can be used to determine the number of failed steel springs.

    2. The maximum vertical dynamic displacement has an excellent linear relationship with the total failure gap, but is relatively unaffected by the specific distribution of the locations of failed springs. Due to the complexity of track parameters, it is difficult to identify the failure gap accurately but it can nevertheless be roughly estimated.

    3. The RMS of the response angles is susceptible to the number, severity, and location of failures. However, the magnitudes of the practical values of the rotation angles are very stable for the specific distribution of failure positions, enabling the failure position to be determined.

    According to the above observations, a set of identification methods for steel spring failure was established. Firstly, the number of failures was determined using the identified 1st-order vertical mode natural frequency and then the failure severity was determined from the maximum vertical dynamic displacement. The failure position was obtained from the RMS of the rotations in each direction. This study considered only the damage location identification of one or two springs.

    Using the above methods, the ISFST was established utilizing the DBSCAN algorithm and statistical analysis and used successfully to detect the failures of steel springs on-site, which were consistent with the field investigation results. After repair, the 1st-order vertical mode natural frequency returned to its pre-damage average level.

    In addition to the steel spring failure mentioned in this study, there are also shear hinge fractures and fastener failures in floating slab track. The identification of these defects is worth further research.

    This work is supported by the National Natural Science Foundation of China (Nos. 51978585 and 52008264), the Applied Basic Research Programs of Science and Technology Commission Foundation of Sichuan Province (No. 2020YJ0214), the Foundation of High-speed Rail Joint Fund Key Projects of Basic Research (No. U1734207), and the Foundation of National Engineering Laboratory for Digital Construction Evaluation Technology of Urban Rail Transit, China (No. 2023JZ01).

    Junyuan ZHENG designed the research. Caiyou ZHAO, Jian WANG, and Bolong JIANG processed the corresponding data. Junyuan ZHENG wrote the first draft of the manuscript. Duojia SHI, Ping WANG, and Xi SHENG helped to organize the manuscript. Junyuan ZHENG revised and edited the final version.

    Junyuan ZHENG, Caiyou ZHAO, Duojia SHI, Ping WANG, Jian WANG, Bolong JIANG, and Xi SHENG declare that they have no conflict of interest.

    Ankrah AA, Kimotho JK, Muvengei OM, 2020. Fusion of model-based and data driven based fault diagnostic methods for railway vehicle suspension., 12(3):51-81. https://doi.org/10.4236/jilsa.2020.123004

    Auersch L, 2017. Static and dynamic behaviours of isolated or unisolated ballast tracks using a fast wavenumber domain method., 87(3):555-574. https://doi.org/10.1007/s00419-016-1209-6

    Bashir S, Akhtar N, 2022. Development of low-frequency mass spring system for underground high-speed railways., 10(2):559-579. https://doi.org/10.1007/S42417-021-00392-W

    Bertha M, Golinval JC, 2017. Identification of non-stationary dynamical systems using multivariate ARMA models., 88:166-179. https://doi.org/10.1016/j.ymssp.2016.11.024

    Chandran P, Thiery F, Odelius J, et al., 2022. Unsupervised machine learning for missing clamp detection from an in-service train using differential eddy current sensor., 14(2):1035. https://doi.org/10.3390/su14021035

    Collier M, 2003. A micro-AGV for flexible manufacturing in small enterprises., 14(5):442-448. https://doi.org/10.1108/09576060310477843

    Cong JL, Gao MY, Wang Y, et al., 2020. Subway rail transit monitoring by built-in sensor platform of smartphone., 21(8):1226-1238. https://doi.org/10.1631/FITEE.1900242

    Cui XL, Chen GX, Yang HG, et al., 2016. Study on rail corrugation of a metro tangential track with Cologne-egg type fasteners., 54(3):353-369. https://doi.org/10.1080/00423114.2015.1137955

    Dersch MS, Khachaturian C, Edwards JR, 2021. Methods to mitigate railway premium fastening system spike fatigue failures using finite element analysis., 121:105160. https://doi.org/10.1016/j.engfailanal.2020.105160

    Dilena M, Limongelli MP, Morassi A, 2015. Damage localization in bridges via the FRF interpolation method., 52-53:162-180. https://doi.org/10.1016/j.ymssp.2014.08.014

    Gomes GF, Giovani RS, 2022. An efficient two-step damage identification method using sunflower optimization algorithm and mode shape curvature (MSDBI–SFO)., 38(2):1711-1730. https://doi.org/10.1007/s00366-020-01128-2

    Hong N, Li LS, Yao WR, et al., 2020. High-speed rail suspension system health monitoring using multi-location vibration data., 21(7):2943-2955. https://doi.org/10.1109/tits.2019.2921785

    James III GH, Carne TG, Laufer JP, 1995. The natural excitation technique (NExT) for modal parameter extraction from operating structures., 10(4):260-277.

    Lam HF, Wong MT, Yang YB, 2012. A feasibility study on railway ballast damage detection utilizing measured vibration of in situ concrete sleeper., 45:284-298. https://doi.org/10.1016/j.engstruct.2012.06.022

    Li SS, 2020. An improved DBSCAN algorithm based on the neighbor similarity and fast nearest neighbor query., 8:47468-47476. https://doi.org/10.1109/ACCESS.2020.2972034

    Li ZW, Liu XZ, Lu HY, et al., 2020. Surface crack detection in precasted slab track in high-speed rail via infrared thermography., 13(21):4837. https://doi.org/10.3390/ma13214837

    Loveday PW, Long CS, Ramatlo DA, 2020. Ultrasonic guided wave monitoring of an operational rail track., 19(6):1666-1684. https://doi.org/10.1177/1475921719893887

    Nelson JT, Watry DL, Amato MA, et al., 2018. Sound transit prototype high performance low frequency floating slab testing and evaluation.: Anderson D, Gautier PE, Iida M, et al. (Eds.), Noise and Vibration Mitigation for Rail Transportation Systems. Springer, Heidelberg, Germany, p.607-618. https://doi.org/10.1007/978-3-319-73411-8_48

    Rajaram S, Nelson JT, 2019. High-performance floating slab track: design and construction improvements based on lessons learned from prototype slabs., 2673(1):300-309. https://doi.org/10.1177/0361198118823004

    Ramos A, Correia AG, Cal?ada R, et al., 2021. Influence of track foundation on the performance of ballast and concrete slab tracks under cyclic loading: physical modelling and numerical model calibration., 277:122245. https://doi.org/10.1016/j.conbuildmat.2021.122245

    Rio G, Soive A, Grolleau V, 2005. Comparative study of numerical explicit time integration algorithms., 36(4):252-265. https://doi.org/10.1016/j.advengsoft.2004.10.011

    Shen ZY, Hedrick JK, Elkins JA, 1983. A comparison of alternative creep force models for rail vehicle dynamic analysis., 12(1-3):79-83. https://doi.org/10.1080/00423118308968725

    Sitharam TG, Sebastian R, Fazil F, 2018. Vibration isolation of buildings housed with sensitive equipment using open trenches–case study and numerical simulations., 115:344-351. https://doi.org/10.1016/j.soildyn.2018.08.033

    Talbot JP, 2016. Base-isolated buildings: towards performance-based design., 169(8):574-582. https://doi.org/10.1680/jstbu.15.00057

    Tamagawa S, 2021. Determination of load test conditions for rail fastenings of a floating slab track., 9(1):14-27. https://doi.org/10.2495/CMEM-V9-N1-14-27

    Vandiver JK, Dunwoody AB, Campbell RB, et al., 1982. A mathematical basis for the random decrement vibration signature analysis technique., 104(2):307-313. https://doi.org/10.1115/1.3256341

    Wang L, Zhao YN, Sang T, et al., 2022. Ultra-low frequency vibration control of urban rail transit: the general quasi-zero-stiffness vibration isolator., 60(5):1788-1805. https://doi.org/10.1080/00423114.2021.1874428

    Wickramasinghe WR, Thambiratnam DP, Chan THT, 2020. Damage detection in a suspension bridge using modal flexibility method., 107:104194. https://doi.org/10.1016/j.engfailanal.2019.104194

    Wilson GP, Saurenman HJ, Nelson JT, 1983. Control of ground-borne noise and vibration., 87(2):339-350. https://doi.org/10.1016/0022-460X(83)90573-4

    Xu FZ, Song XL, Yang JJ, 2020. Influence of steel spring failure of floating slab track on vibration characteristics of infrastructure.: Tutumluer E, Chen XB, Xiao YJ (Eds.), Advances in Environmental Vibration and Transportation Geodynamics. Springer, Singapore, p.987-998. https://doi.org/10.1007/978-981-15-2349-6_67

    Yu P, Manalo A, Ferdous W, et al., 2021. Failure analysis and the effect of material properties on the screw pull-out behaviour of polymer composite sleeper materials., 128:105577. https://doi.org/10.1016/j.engfailanal.2021.105577

    Yuan XC, Zhu SY, Xu L, et al., 2020. Investigation of the vibration isolation performance of floating slab track with rubber bearings using a stochastic fractional derivative model., 234(9):992-1004. https://doi.org/10.1177/0954409719883552

    Zhao CY, Liu DY, Zhang XM, et al., 2019. Influence of vibration isolator failure on vehicle operation performance and floating slab track structure vibration reduction effectiveness., 2019:8385310. https://doi.org/10.1155/2019/8385310

    Zhao CY, Zheng JY, Sang T, et al., 2021. Computational analysis of phononic crystal vibration isolators via FEM coupled with the acoustic black hole effect to attenuate railway-induced vibration., 283:122802. https://doi.org/10.1016/j.conbuildmat.2021.122802

    Zhu SY, Zhang QL, Zhai WM, et al., 2021. Sensor deploying for damage identification of vibration isolator in floating-slab track using deep residual network., 183:109801. https://doi.org/10.1016/j.measurement.2021.109801

    Zou JH, Du TF, Chen W, et al., 2022. Experimental study of concrete floating slab municipal road with steel spring isolators under vehicle loads., 315:125686. https://doi.org/10.1016/j.conbuildmat.2021.125686

    城市軌道交通鋼彈簧浮置板軌道支承剛度失效識(shí)別方法:一項(xiàng)中國的案例研究

    鄭鈞元1,2,趙才友1,2,師多佳1,2,王平1,2,王建3,姜博龍4,盛曦5

    1西南交通大學(xué),高速鐵路線路工程教育部重點(diǎn)實(shí)驗(yàn)室,中國成都,610031;2西南交通大學(xué),土木工程學(xué)院,中國成都,610031;3青島零一動(dòng)測(cè)數(shù)據(jù)科技有限公司,中國青島,266000;4中國鐵路設(shè)計(jì)集團(tuán)有限公司,城市軌道交通數(shù)字化建設(shè)與測(cè)評(píng)技術(shù)國家工程研究中心,中國天津,300308;5深圳大學(xué),城市智慧交通與安全運(yùn)維研究院,中國深圳,518060

    在列車往復(fù)動(dòng)力作用和復(fù)雜環(huán)境因素耦合作用下,常出現(xiàn)諸如鋼彈簧隔振器斷裂、基底不均勻沉降誘發(fā)吊空等隔振器支承剛度失效的問題。本文提出一種鋼彈簧浮置板軌道支承剛度失效的方法,能及時(shí)識(shí)別失效數(shù)量、失效程度、失效位置,以避免剛度失效引起的軌道結(jié)構(gòu)劣化和行車安全威脅。

    1. 分析隔振器支承剛度失效對(duì)浮置板軌道系統(tǒng)各動(dòng)力指標(biāo)的影響,并構(gòu)建一套浮置板軌道鋼彈簧失效識(shí)別方法;2. 結(jié)合具有噪聲的基于密度的聚類方法(DBSCAN)算法和統(tǒng)計(jì)分析構(gòu)建一套浮置板軌道支承剛度智能監(jiān)測(cè)系統(tǒng),并應(yīng)用于中國南部某線路。

    1. 通過動(dòng)力響應(yīng)敏感性分析,得出浮置板自振頻率、動(dòng)位移、轉(zhuǎn)角對(duì)支承剛度失效的敏感性較高(圖5和6);2. 通過分析多種剛度失效的情況,得出敏感性指標(biāo)和支承剛度失效的關(guān)系(表3和4,圖8);3. 基于上述成果,構(gòu)建智能監(jiān)測(cè)系統(tǒng)并應(yīng)用到現(xiàn)場,驗(yàn)證方法的有效性(圖13和15)。

    1. 動(dòng)力學(xué)模型的敏感性分析表明,列車荷載下的浮置板一階垂向自振頻率、最大垂向動(dòng)位移、轉(zhuǎn)角有效值易受剛度失效影響。2. 浮置板的一階垂向頻率主要跟失效數(shù)量有關(guān);浮置板最大垂向動(dòng)位移與失效高度具有較好的線性關(guān)系;浮置板轉(zhuǎn)角之間的大小關(guān)系對(duì)失效位置很敏感。3. 以上述方法為理論基礎(chǔ)構(gòu)建浮置板軌道支承剛度智能監(jiān)測(cè)系統(tǒng),成功識(shí)別出鋼彈簧失效。

    浮置板軌道;支承剛度;精細(xì)化模型;失效識(shí)別;監(jiān)測(cè)系統(tǒng)

    https://doi.org/10.1631/jzus.A2300085

    https://doi.org/10.1631/jzus.A2300085

    Duojia SHI, https://orcid.org/0000-0001-5081-7973

    Received Feb. 14, 2023;

    Revision accepted June 26, 2023;

    Crosschecked Jan. 3, 2024

    ? Zhejiang University Press 2024

    猜你喜歡
    浮置轉(zhuǎn)角敏感性
    玩轉(zhuǎn)角的平分線
    移動(dòng)荷載作用下離散支承浮置板軌道振動(dòng)響應(yīng)研究
    釔對(duì)Mg-Zn-Y-Zr合金熱裂敏感性影響
    三次“轉(zhuǎn)角”遇到愛
    解放軍健康(2017年5期)2017-08-01 06:27:42
    永春堂贏在轉(zhuǎn)角
    AH70DB鋼焊接熱影響區(qū)組織及其冷裂敏感性
    焊接(2016年1期)2016-02-27 12:55:37
    如何培養(yǎng)和提高新聞敏感性
    新聞傳播(2015年8期)2015-07-18 11:08:24
    預(yù)制浮置板軌道隔振性能研究
    下一個(gè)轉(zhuǎn)角:邁出去 開啟“智”造時(shí)代
    基于多島遺傳算法的橡膠浮置板彈性的優(yōu)化
    久久欧美精品欧美久久欧美| 黄频高清免费视频| 日本五十路高清| 女同久久另类99精品国产91| 国产成+人综合+亚洲专区| 一级毛片精品| 免费一级毛片在线播放高清视频 | 欧美日韩瑟瑟在线播放| 日日摸夜夜添夜夜添小说| 免费观看精品视频网站| xxx96com| 午夜久久久久精精品| 级片在线观看| 亚洲第一电影网av| 国产成人啪精品午夜网站| 精品国产一区二区久久| 一本久久中文字幕| av网站免费在线观看视频| 最近最新中文字幕大全免费视频| av在线天堂中文字幕| 色综合站精品国产| 亚洲人成77777在线视频| 成人亚洲精品一区在线观看| 在线天堂中文资源库| 欧美黑人欧美精品刺激| 在线永久观看黄色视频| avwww免费| 欧美午夜高清在线| 99国产精品免费福利视频| 免费高清视频大片| 亚洲第一青青草原| 中文字幕人成人乱码亚洲影| 丝袜在线中文字幕| 国产99久久九九免费精品| 国产精品秋霞免费鲁丝片| 亚洲精华国产精华精| 美女大奶头视频| 九色亚洲精品在线播放| 国产精品野战在线观看| 级片在线观看| 欧美中文综合在线视频| 久久久国产成人免费| 狠狠狠狠99中文字幕| 日韩欧美三级三区| av欧美777| 亚洲成人精品中文字幕电影| 亚洲专区国产一区二区| 久久国产精品人妻蜜桃| 日韩有码中文字幕| 免费看美女性在线毛片视频| 亚洲熟妇中文字幕五十中出| 国产一区在线观看成人免费| 久久欧美精品欧美久久欧美| 亚洲国产高清在线一区二区三 | 一级毛片精品| 日韩欧美在线二视频| 在线观看免费午夜福利视频| 大型av网站在线播放| 久久久久久大精品| 久久亚洲精品不卡| 亚洲自拍偷在线| 国产一级毛片七仙女欲春2 | 99久久综合精品五月天人人| 久9热在线精品视频| 国产精品日韩av在线免费观看 | 97人妻精品一区二区三区麻豆 | 精品免费久久久久久久清纯| 日韩 欧美 亚洲 中文字幕| 国产精品综合久久久久久久免费 | 激情在线观看视频在线高清| 看免费av毛片| 国产精品一区二区在线不卡| 少妇的丰满在线观看| 欧美色视频一区免费| 99riav亚洲国产免费| 波多野结衣一区麻豆| 一级作爱视频免费观看| 两人在一起打扑克的视频| 两个人视频免费观看高清| 国产成人精品久久二区二区91| 色播亚洲综合网| 亚洲九九香蕉| 波多野结衣av一区二区av| 欧美黑人欧美精品刺激| 在线观看免费视频日本深夜| 在线十欧美十亚洲十日本专区| 日韩欧美一区二区三区在线观看| 国产精品1区2区在线观看.| 美国免费a级毛片| 亚洲欧洲精品一区二区精品久久久| 亚洲国产中文字幕在线视频| 国产精品亚洲av一区麻豆| 两人在一起打扑克的视频| 国产成人欧美在线观看| 欧美黄色淫秽网站| 国产精品美女特级片免费视频播放器 | 国产亚洲av高清不卡| 岛国在线观看网站| 超碰成人久久| 美国免费a级毛片| 中文字幕精品免费在线观看视频| 久久久久久久久中文| 亚洲 欧美一区二区三区| 高清在线国产一区| 老司机午夜福利在线观看视频| 亚洲成国产人片在线观看| 欧美黄色淫秽网站| 国产麻豆69| 国产欧美日韩综合在线一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| tocl精华| 亚洲精品久久国产高清桃花| 亚洲人成77777在线视频| 国产精品自产拍在线观看55亚洲| 91麻豆精品激情在线观看国产| 一级a爱片免费观看的视频| 亚洲自拍偷在线| 男女午夜视频在线观看| av中文乱码字幕在线| 神马国产精品三级电影在线观看 | av中文乱码字幕在线| 久久国产精品影院| 老司机深夜福利视频在线观看| 搡老熟女国产l中国老女人| 黑人操中国人逼视频| 久久人人爽av亚洲精品天堂| 丰满人妻熟妇乱又伦精品不卡| 欧美另类亚洲清纯唯美| 国产亚洲欧美98| 麻豆av在线久日| 母亲3免费完整高清在线观看| 亚洲七黄色美女视频| 日韩欧美免费精品| 露出奶头的视频| 变态另类丝袜制服| 在线观看免费视频网站a站| 亚洲av成人不卡在线观看播放网| 久久久久久人人人人人| 亚洲成人久久性| 亚洲专区中文字幕在线| 九色亚洲精品在线播放| 欧美 亚洲 国产 日韩一| 不卡av一区二区三区| 久久精品亚洲熟妇少妇任你| 777久久人妻少妇嫩草av网站| 色婷婷久久久亚洲欧美| www.精华液| 夜夜爽天天搞| 国产成人精品在线电影| АⅤ资源中文在线天堂| 欧美日本亚洲视频在线播放| 人人妻人人澡欧美一区二区 | 婷婷精品国产亚洲av在线| 男女午夜视频在线观看| 亚洲国产精品999在线| 亚洲熟女毛片儿| 变态另类成人亚洲欧美熟女 | 一区二区三区高清视频在线| 麻豆成人av在线观看| 又紧又爽又黄一区二区| 亚洲精品一卡2卡三卡4卡5卡| 手机成人av网站| 激情在线观看视频在线高清| 欧美成人午夜精品| 国产精品99久久99久久久不卡| 黄色成人免费大全| 精品久久久精品久久久| 麻豆一二三区av精品| 夜夜躁狠狠躁天天躁| 成人欧美大片| 国产熟女xx| 一本综合久久免费| 9191精品国产免费久久| 欧美激情久久久久久爽电影 | 久久人人精品亚洲av| 精品人妻1区二区| 91麻豆精品激情在线观看国产| 午夜亚洲福利在线播放| tocl精华| 精品一区二区三区四区五区乱码| 99香蕉大伊视频| 国产97色在线日韩免费| 国产国语露脸激情在线看| 欧美黄色淫秽网站| 欧美成人免费av一区二区三区| 成年女人毛片免费观看观看9| 午夜老司机福利片| 极品人妻少妇av视频| 国产单亲对白刺激| 亚洲avbb在线观看| 麻豆一二三区av精品| 国产精品国产高清国产av| 亚洲一区中文字幕在线| 激情视频va一区二区三区| 中文亚洲av片在线观看爽| 老汉色∧v一级毛片| 久久久久久人人人人人| 日韩 欧美 亚洲 中文字幕| 精品人妻在线不人妻| 亚洲熟妇熟女久久| 国产高清激情床上av| 99国产精品免费福利视频| 亚洲成av片中文字幕在线观看| 成人免费观看视频高清| 人成视频在线观看免费观看| 国产精华一区二区三区| 狠狠狠狠99中文字幕| cao死你这个sao货| 嫁个100分男人电影在线观看| 高清黄色对白视频在线免费看| 久久久久久久久久久久大奶| 叶爱在线成人免费视频播放| 麻豆av在线久日| 美女 人体艺术 gogo| svipshipincom国产片| xxx96com| 亚洲av片天天在线观看| 亚洲自偷自拍图片 自拍| 97人妻天天添夜夜摸| 欧美成人午夜精品| 人人妻人人澡欧美一区二区 | 亚洲第一av免费看| 一区二区日韩欧美中文字幕| 变态另类丝袜制服| 亚洲人成77777在线视频| 国产高清有码在线观看视频 | 免费女性裸体啪啪无遮挡网站| 免费女性裸体啪啪无遮挡网站| 18美女黄网站色大片免费观看| 久久香蕉激情| 国内精品久久久久久久电影| 熟女少妇亚洲综合色aaa.| 久久午夜综合久久蜜桃| 色播在线永久视频| 超碰成人久久| 动漫黄色视频在线观看| 又紧又爽又黄一区二区| 伊人久久大香线蕉亚洲五| av中文乱码字幕在线| 少妇被粗大的猛进出69影院| 老汉色av国产亚洲站长工具| 日韩av在线大香蕉| 操出白浆在线播放| 精品午夜福利视频在线观看一区| 国产精品一区二区免费欧美| 夜夜躁狠狠躁天天躁| 韩国av一区二区三区四区| 侵犯人妻中文字幕一二三四区| 日韩免费av在线播放| 神马国产精品三级电影在线观看 | 一a级毛片在线观看| 午夜免费成人在线视频| 精品人妻在线不人妻| 亚洲av电影在线进入| 一本久久中文字幕| 操美女的视频在线观看| xxx96com| 又黄又粗又硬又大视频| 国产熟女xx| 午夜激情av网站| 男人舔女人的私密视频| 91成人精品电影| 午夜福利,免费看| 亚洲狠狠婷婷综合久久图片| 桃红色精品国产亚洲av| 首页视频小说图片口味搜索| 亚洲色图综合在线观看| 欧美另类亚洲清纯唯美| 久久人人爽av亚洲精品天堂| 成人欧美大片| 欧美激情久久久久久爽电影 | 亚洲av成人av| 女性生殖器流出的白浆| 黄色视频不卡| 中亚洲国语对白在线视频| 日韩大码丰满熟妇| 十分钟在线观看高清视频www| 午夜成年电影在线免费观看| 黄色毛片三级朝国网站| 国产精品一区二区精品视频观看| 免费在线观看影片大全网站| 黄色女人牲交| 欧美日韩亚洲综合一区二区三区_| 亚洲av熟女| 男人操女人黄网站| 免费久久久久久久精品成人欧美视频| 国产激情欧美一区二区| 变态另类成人亚洲欧美熟女 | 婷婷六月久久综合丁香| АⅤ资源中文在线天堂| 母亲3免费完整高清在线观看| 操出白浆在线播放| 日本在线视频免费播放| 国产亚洲精品av在线| 少妇粗大呻吟视频| 免费av毛片视频| 69av精品久久久久久| 热99re8久久精品国产| 制服人妻中文乱码| 日韩国内少妇激情av| 999精品在线视频| 日韩精品免费视频一区二区三区| 国产精品爽爽va在线观看网站 | 又黄又粗又硬又大视频| 夜夜看夜夜爽夜夜摸| 国产在线观看jvid| or卡值多少钱| 欧美乱妇无乱码| 久久久久国产一级毛片高清牌| 在线av久久热| 窝窝影院91人妻| 香蕉久久夜色| 黑人操中国人逼视频| 久久婷婷成人综合色麻豆| 两个人免费观看高清视频| 国产黄a三级三级三级人| 久久精品91无色码中文字幕| 免费久久久久久久精品成人欧美视频| 男男h啪啪无遮挡| 国产午夜福利久久久久久| 国产成人系列免费观看| 一区福利在线观看| 99国产综合亚洲精品| 久久九九热精品免费| 大香蕉久久成人网| 亚洲第一欧美日韩一区二区三区| 9191精品国产免费久久| 一进一出抽搐gif免费好疼| 一a级毛片在线观看| 国产av一区二区精品久久| 日本免费a在线| av片东京热男人的天堂| 日韩欧美国产在线观看| 午夜激情av网站| 国产主播在线观看一区二区| 黑丝袜美女国产一区| 一区二区三区精品91| 亚洲av电影不卡..在线观看| 国产精品一区二区在线不卡| 国产精品一区二区精品视频观看| 嫩草影视91久久| 美女高潮到喷水免费观看| 久9热在线精品视频| 国产精品一区二区在线不卡| 久久香蕉精品热| 久久亚洲精品不卡| 亚洲aⅴ乱码一区二区在线播放 | 精品免费久久久久久久清纯| 久99久视频精品免费| 久久伊人香网站| www日本在线高清视频| 啦啦啦免费观看视频1| 国产av精品麻豆| 黄色片一级片一级黄色片| 黄色a级毛片大全视频| 激情在线观看视频在线高清| 婷婷精品国产亚洲av在线| 欧美日本中文国产一区发布| 男人操女人黄网站| 色尼玛亚洲综合影院| 国产精品国产高清国产av| 身体一侧抽搐| 国产精品免费一区二区三区在线| 国产成人影院久久av| 亚洲欧美精品综合久久99| 欧美绝顶高潮抽搐喷水| 国产三级在线视频| 亚洲男人天堂网一区| 电影成人av| 亚洲国产精品久久男人天堂| 亚洲色图av天堂| 亚洲欧美日韩无卡精品| 嫩草影院精品99| 国产精品美女特级片免费视频播放器 | 国产亚洲av高清不卡| 久久精品国产亚洲av香蕉五月| 一边摸一边抽搐一进一出视频| 成人国产一区最新在线观看| 欧美黄色片欧美黄色片| 大香蕉久久成人网| 最近最新中文字幕大全电影3 | 欧美激情 高清一区二区三区| 搡老妇女老女人老熟妇| 欧美一级毛片孕妇| 国产精品乱码一区二三区的特点 | 亚洲欧美一区二区三区黑人| 无人区码免费观看不卡| 久久草成人影院| 又黄又爽又免费观看的视频| 日韩成人在线观看一区二区三区| 亚洲无线在线观看| 在线观看免费午夜福利视频| 国产亚洲欧美精品永久| 成人永久免费在线观看视频| 亚洲国产欧美日韩在线播放| 国产一区二区在线av高清观看| 日韩有码中文字幕| 久久青草综合色| 午夜福利一区二区在线看| 99久久国产精品久久久| 又紧又爽又黄一区二区| 国产精品九九99| 欧美中文日本在线观看视频| 激情视频va一区二区三区| 黄色视频,在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 999精品在线视频| 制服丝袜大香蕉在线| 人妻丰满熟妇av一区二区三区| 搡老岳熟女国产| 如日韩欧美国产精品一区二区三区| 亚洲国产日韩欧美精品在线观看 | 好看av亚洲va欧美ⅴa在| 精品国产美女av久久久久小说| 欧美日韩一级在线毛片| 午夜激情av网站| 日韩视频一区二区在线观看| www.精华液| 成人国产一区最新在线观看| 亚洲五月天丁香| 欧美乱码精品一区二区三区| 欧美老熟妇乱子伦牲交| 性欧美人与动物交配| 国产99白浆流出| 色精品久久人妻99蜜桃| 黑人操中国人逼视频| 一区福利在线观看| 亚洲色图 男人天堂 中文字幕| 亚洲人成电影观看| 欧美久久黑人一区二区| 国产成人欧美| 国产亚洲精品久久久久5区| 99精品久久久久人妻精品| 精品高清国产在线一区| 自线自在国产av| 两个人视频免费观看高清| 亚洲精品国产一区二区精华液| 侵犯人妻中文字幕一二三四区| 岛国视频午夜一区免费看| 老鸭窝网址在线观看| 国产精品一区二区在线不卡| 欧美国产精品va在线观看不卡| av欧美777| 久久久久国内视频| 国产精品一区二区免费欧美| 99久久久亚洲精品蜜臀av| 国产精品日韩av在线免费观看 | 女性被躁到高潮视频| 动漫黄色视频在线观看| 欧美日本中文国产一区发布| 热99re8久久精品国产| 日日摸夜夜添夜夜添小说| 美女大奶头视频| 岛国视频午夜一区免费看| 在线观看日韩欧美| 国产亚洲精品综合一区在线观看 | 91麻豆精品激情在线观看国产| 给我免费播放毛片高清在线观看| 欧美中文综合在线视频| 久热爱精品视频在线9| АⅤ资源中文在线天堂| aaaaa片日本免费| 美女国产高潮福利片在线看| 国产精品免费一区二区三区在线| av福利片在线| 天天一区二区日本电影三级 | 夜夜爽天天搞| 成人18禁在线播放| 每晚都被弄得嗷嗷叫到高潮| 亚洲熟妇熟女久久| 日本vs欧美在线观看视频| 亚洲国产精品合色在线| 美女午夜性视频免费| 免费观看人在逋| 欧美日韩福利视频一区二区| 性少妇av在线| 级片在线观看| 99精品在免费线老司机午夜| 91在线观看av| 又黄又粗又硬又大视频| 国语自产精品视频在线第100页| 国产一区在线观看成人免费| 欧美成人一区二区免费高清观看 | 黑人操中国人逼视频| 国产精品久久久久久亚洲av鲁大| 精品久久久久久久人妻蜜臀av | 99久久99久久久精品蜜桃| 亚洲第一av免费看| 在线观看午夜福利视频| 中文字幕高清在线视频| 亚洲成人精品中文字幕电影| 视频在线观看一区二区三区| 亚洲色图综合在线观看| 亚洲国产精品sss在线观看| 嫁个100分男人电影在线观看| 高清在线国产一区| 国产成人欧美| 国产精品精品国产色婷婷| 午夜久久久久精精品| 亚洲中文字幕日韩| 亚洲熟妇熟女久久| 久久久久久亚洲精品国产蜜桃av| 成年版毛片免费区| 亚洲精品国产精品久久久不卡| 精品久久久久久成人av| 精品一区二区三区视频在线观看免费| 亚洲性夜色夜夜综合| 级片在线观看| 亚洲 欧美 日韩 在线 免费| 欧美日韩亚洲综合一区二区三区_| 精品欧美一区二区三区在线| 国产亚洲精品久久久久久毛片| 亚洲视频免费观看视频| 老司机午夜十八禁免费视频| 国产精品一区二区精品视频观看| 美女高潮到喷水免费观看| 亚洲第一欧美日韩一区二区三区| 不卡一级毛片| 一个人免费在线观看的高清视频| 亚洲av电影不卡..在线观看| 不卡一级毛片| 午夜亚洲福利在线播放| 午夜久久久久精精品| 亚洲精品在线观看二区| 18美女黄网站色大片免费观看| 色播在线永久视频| 午夜a级毛片| 国产精品久久视频播放| 久久青草综合色| 啦啦啦免费观看视频1| 亚洲九九香蕉| 精品无人区乱码1区二区| 老司机深夜福利视频在线观看| 99热只有精品国产| a在线观看视频网站| 麻豆一二三区av精品| 亚洲片人在线观看| 成年女人毛片免费观看观看9| 18禁观看日本| 男女床上黄色一级片免费看| 国产精品久久电影中文字幕| 99re在线观看精品视频| 日本三级黄在线观看| 精品一品国产午夜福利视频| 色av中文字幕| 99re在线观看精品视频| 亚洲avbb在线观看| 国产不卡一卡二| 欧美日韩亚洲国产一区二区在线观看| 午夜福利欧美成人| 午夜久久久在线观看| 精品人妻在线不人妻| 国产人伦9x9x在线观看| 18禁国产床啪视频网站| 制服人妻中文乱码| 国产成年人精品一区二区| 午夜免费鲁丝| 国产一区二区激情短视频| 免费在线观看日本一区| 老汉色av国产亚洲站长工具| 看免费av毛片| 亚洲av日韩精品久久久久久密| 国产亚洲精品综合一区在线观看 | 一级毛片精品| 国产一区二区激情短视频| 男男h啪啪无遮挡| 国产成人精品在线电影| 久久精品国产综合久久久| 叶爱在线成人免费视频播放| 欧美乱色亚洲激情| 久久精品国产亚洲av高清一级| 老熟妇乱子伦视频在线观看| 亚洲精华国产精华精| 国产亚洲精品久久久久久毛片| 欧美成狂野欧美在线观看| 老司机靠b影院| 美女 人体艺术 gogo| 亚洲精品一区av在线观看| 欧美黄色淫秽网站| 久久国产精品男人的天堂亚洲| 亚洲国产中文字幕在线视频| 亚洲视频免费观看视频| 精品久久久久久久久久免费视频| 黄色视频不卡| 99久久99久久久精品蜜桃| 欧美日韩亚洲综合一区二区三区_| 天天一区二区日本电影三级 | 日韩国内少妇激情av| 亚洲av第一区精品v没综合| 激情在线观看视频在线高清| 18美女黄网站色大片免费观看| 国产一区二区三区视频了| 人人妻人人澡人人看| 涩涩av久久男人的天堂| 狂野欧美激情性xxxx| 一区二区日韩欧美中文字幕| 欧美久久黑人一区二区| 色综合站精品国产| 无限看片的www在线观看| 欧美在线一区亚洲| 国产亚洲av高清不卡| 麻豆国产av国片精品| 午夜免费成人在线视频| 婷婷丁香在线五月| 一区二区日韩欧美中文字幕| 国产在线精品亚洲第一网站| av有码第一页| 国产亚洲精品一区二区www| 亚洲七黄色美女视频| 黄色成人免费大全| 如日韩欧美国产精品一区二区三区| 天天躁夜夜躁狠狠躁躁| 首页视频小说图片口味搜索| 亚洲国产看品久久| 久久久久久久久免费视频了|