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

    Influence of ground effect on flow field structure and aerodynamic noise of high-speed trains

    2024-02-05 01:02:54XiaomingTANLinliGONGXiaohongZHANGZhigangYANG
    關(guān)鍵詞:影響

    XiaomingTAN,LinliGONG,XiaohongZHANG,ZhigangYANG

    Research Article

    Influence of ground effect on flow field structure and aerodynamic noise of high-speed trains

    1Key Laboratory of Intelligent Manufacturing and Service Performance Optimization of Laser and Grinding in Mechanical Industry, Hunan Institute of Science and Technology, Yueyang 414000, China2Key Laboratory of Traffic Safety on Track (Central South University), Ministry of Education, School of Traffic & Transportation Engineering, Central South University, Changsha 410075, China

    The simulation of the ground effect has always been a technical difficulty in wind tunnel tests of high-speed trains. In this paper, large eddy simulation and the curl acoustic integral equation were used to simulate the flow-acoustic field results of high-speed trains under four ground simulation systems (GSSs): “moving ground+rotating wheel”, “stationary ground+rotating wheel”, “moving ground+stationary wheel”, and “stationary ground+stationary wheel”. By comparing the fluid-acoustic field results of the four GSSs, the influence laws of different GSSs on the flow field structure, aero-acoustic source, and far-field radiation noise characteristics were investigated, providing guidance for the acoustic wind tunnel testing of high-speed trains. The calculation results of the aerodynamic noise of a 350 km/h high-speed train show that the moving ground and rotating wheel affect mainly the aero-acoustic performance under the train bottom. The influence of the rotating wheel on the equivalent sound source power of the whole vehicle was not more than 5%, but that of the moving ground slip was more than 15%. The average influence of the rotating wheel on the sound pressure level radiated by the whole vehicle was 0.3 dBA, while that of the moving ground was 1.8 dBA.

    High-speed train; Aero-acoustics; Flow field structure; Large eddy simulation; Moving ground condition; Rotating wheel

    1 Introduction

    As high-speed ground transportation vehicles, high-speed trains inevitably encounter ground effect problems, that is, disturbance caused by the ground boundary layer to the structure of the train bottom flow fields. The accuracy of simulation of the ground effect phenomenon has a great effect on the accuracy of predictions of high-speed train aerodynamics and aerodynamic noise.

    Prediction models are widely used at present to understand the nature of railway ground waves and their associated characteristics (Kouroussis et al., 2021). Using boundary layer correction theory to correct the simulation results of aerodynamic force and aerodynamic noise of 3D complex structures such as high-speed trains is difficult. The main techniques available to eliminate the negative effects of ground effect in wind tunnel tests include the moving ground method, tangent blow-suction method, and suction method. The moving ground method is the most satisfactory for its capacity to simulate the relative motion between the train, ground, and air. Using the moving ground method, Tyll et al. (1996) estimated the aerodynamic drag coefficients of maglev to be 0.24 with moving rails and 0.17 without rails. Moreover, they pointed out that a moving track is very important to the accuracy of aerodynamic testing of a maglev train. Yi et al. (1997) investigated the ground effect of a high-speed train using the suction method. They showed that the elimination of the boundary layer can significantly increase the aerodynamic drag force coefficient, the increment of which depends on the shape of the train bottom. By comparing the test results from the tangent blow-suction method and the moving ground method, Kwon et al. (2001) concluded that the tangential blowing method was an effective substitute for the moving ground method for testing aerodynamics.

    Xia et al. (2016) compared train wind simulation results between the stationary ground and the moving ground conditions and concluded that the stationary ground result was larger than that of the moving ground. In the following year, they constructed seven ground simulation systems (GSSs) and pointed out that the elimination of the boundary layer can significantly increase the aerodynamic drag coefficient, and that elevation of the model cannot effectively eliminate the adverse effects of the ground effect (Xia et al., 2017). Zhang et al. (2016) numerically simulated the train aerodynamic simulation results of three GSSs: "stationary ground+stationary wheel", "moving ground+stationary wheel", and "moving ground+rotating wheel", and concluded that the moving ground could significantly increase the aerodynamic drag coefficient, while the rotating wheel could not. Paz et al. (2017) evaluated the influence of sleepers on the aerodynamic drag coefficient on the ground effect through the dynamic layering method. They found that the existence of sleepers had a great impact on the aerodynamic performance of the whole train, increasing the aerodynamic drag coefficient by about 15%, but had little impact on the lift coefficient. Zhu et al. (2017) numerically calculated the aerodynamic noise radiated by a wheelset of two GSS conditions: moving ground and no ground. They found that the result of moving ground was about 7 dB higher than that of no ground. The simulation results of aerodynamic noise of high-speed trains by Liu et al. (2013) showed that the result of moving ground was 4?–?6 dBA larger than that of stationary ground. Wen et al. (2019) found that the existence of bogies on the bottom of the train, especially the last bogie, not only enhanced the wake flow but also introduced large perturbances into the wake flow. The bogie effects were revealed through a systematic comparison of flow structures, slipstream characteristics, and aerodynamic forces between two generic train configurations (Shibo et al., 2018).

    To sum up, we now have a more comprehensive and accurate understanding of how the ground effect influences the prediction accuracy of aerodynamic force, but how the ground effect influences the prediction accuracy of aerodynamic noise is still unclear. The reason is that the current methods of eliminating the ground boundary layer cannot account for either changes in the fluctuating flow field at the train bottom or the generation of new noise. The ground is assumed to be stationary in current acoustic wind tunnel tests, and it is impossible to explore the influence laws of different GSSs on aero-acoustic tests. On the other hand, the current geometric models for numerical simulation are simplified. In this study, we constructed numerical models of four GSSs to explore the influence laws of different GSSs on the flow field structure and aerodynamic noise, and to provide guidance for the acoustic wind tunnel testing of high-speed trains.

    2 Numerical simulation model

    The geometric model was a 1/8-scale high-speed train model with three carriages, bogies and without a pantograph (Tan et al., 2018). The model was 7960 mm long, 4080 mm high, and 3360 mm wide. The geometric model was divided into 32 parts, i.e., Components #1?–?#32. The naming rules of the first 20 parts are shown in Fig. S1 of the electronic supplementary materials (ESM), and the last 12 parts correspond to the 12 wheelsets. From upstream to downstream, the wheelset was named sequentially from Component #21. For example, Component #21 corresponds to the front wheelset of Bogie 1, and Component #22 corresponds to the back wheelset of Bogie 1.

    The simulated incoming velocity was 350 km/h. The surface of the train model, excluding the 12 wheelsets, was set as a non-slip boundary condition with friction. All wheelsets were set to a speed of 1817.2 or 0 rad/s. The wheel speed is derived from Section S1 of the ESM. The ground moving velocity was 350 or 0 km/h. These could be combined freely to form four GSSs: "moving ground+rotating wheel", "stationary ground+rotating wheel", "moving ground+stationary wheel", and "stationary ground+stationary wheel". For the convenience of discussion, they were named Case 1, Case 2, Case 3, and Case 4, respectively.

    We used the Ansys Fluent software from the Wuxi Computing Centre (China) to carry out numerical simulations of the aerodynamic noise of the high-speed train. The unsteady calculation time step was 5×10-5s, and iterated 30 steps in each time step. A total of 10000 time steps were calculated. The computation of the first 2500 time steps was to ensure full development of the turbulence flow field, and that of the remaining 7500 time steps was to extract the noise source information. Sound source data in each time step were stored, with a time span of 0.375 s and a frequency resolution of 2.7 Hz. The number of body grids was about 110 million.

    Fig. 1 Grid distribution: (a) longitudinal symmetry plane; (b) on the same elevation as the tip of the nose; (c) a cross section of the tail-flow type shoulder; (d) streamlined area; (e) bogie

    3 Difference analysis of flow field structure characteristics

    Through the three physical quantities of velocity amplitude, vorticity amplitude, andvalue, the differences among the four GSSs in terms of the flow field structure around the train were revealed. For convenience of description, Fig. 2 shows all sections in this study. Section a is a-equivalent plane running through the calculation domain and the wheel pair, with avalue of 93 mm. Section b is theiso-surface penetrating the calculation domain, with avalue of 50 mm. Line 1 runs along thedirection and penetrates the calculation region, with avalue of 93 mm and avalue of 12.5 mm.

    Fig. 2 Sketch map of the section

    3.1 Flow velocity magnitude

    From Figs. 3 and 4, the influence of the four GSSs on the spatial distribution of the velocity amplitude around the high-speed train can be obtained. The four cases had no influence on the velocity amplitude distribution in Region I, the upper part of the train. In Region II, near the nose of the head car, whether the wheel was rotating or not had little effect on the form of the stagnation zone, and the volume of the stagnation zone caused by "stationary ground" was larger than that caused by "moving ground". In Region III, the area from the discharger of the head car to Bogie 1, the length of the acceleration zone under the "moving ground" condition was larger than that under the "stationary ground" condition. When the wheel was rotating, the velocity amplitude inside the cavity was increased. In Region IV, the area from beneath Bogie 2 to Bogie 3, the velocity amplitude under "moving ground" was larger than that under "stationary ground", and so was the velocity amplitude inside the cavity. In Region V, the wake area near the discharger of the tail car, the low-speed zone in "moving ground", was smaller than that in "stationary ground". Nevertheless, whether the wheel was rotating or not hardly affected the velocity amplitude distribution in this region. Note that the velocity amplitude distribution of the ground area under "stationary ground" was closely related to the grid scale and ground boundary layer effect.

    Fig. 3 Instantaneous velocity amplitude contour of Section a: (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4

    Fig. 4 Instantaneous velocity amplitude contour of Section b: (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4

    To quantitatively analyze the velocity amplitude distribution formed in the four GSSs, Table 1 shows the mean value of the velocity amplitude of Line 1.

    Table 1 Mean value of the velocity amplitude

    Whether the wheel was rotating or not hardly affected the velocity amplitude of Line 1, while the velocity amplitude of Line 1 on moving ground was 23% larger than that on stationary ground (Table 1).

    3.2 Vorticity magnitude

    From Figs. 5 and 6, we can obtain the effects of the four GSSs on the spatial distribution of vorticity amplitude around high-speed trains. The four cases had no effect on the vorticity amplitude distribution in Region I, the upper part of the train. In Region II, near the nose of the head car, the state of the wheel had little effect on the vorticity distribution. However, the vorticity distribution on "stationary ground" was uneven, but relatively uniform on "moving ground". In Region III, the area from the discharger of the head car to Bogie 1, the "moving ground" decreased the vorticity amplitude between the train and the ballast compared with the "stationary ground". In Region IV, the area from beneath Bogie 2 to Bogie 3, whether the ground was moving or not had little effect on the vorticity distribution in that region. In Region V, the wake area near the discharger of the tail car, both "stationary ground" and "moving ground" conditions could form strong vorticity separation flow at the shoulder of the tail car streamline and form a strong vortex mixing flow in the tail car cowcatcher bottom and downstream area. But the separation flow formed in the "stationary ground" condition was closer to the streamlined car body downstream, and the mixed flow formed by "stationary ground" showed large-scale rotating flow intensity and strong swing strength on both sides. Besides, whether the wheel was rotating or not hardly affected the vorticity distribution of this region. Note that the amplitude distribution of vorticity was closely related to the grid scale and the boundary layer effect. The larger the grid scale, i.?e, the further away from the high-speed train area, the more uneven the cloud map of vorticity amplitude distribution, and the larger the voracity amplitude. Nonetheless, these phenomena did not appear under the "moving ground" condition.

    Fig. 5 Instantaneous contour plot of vorticity amplitude at Section a: (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4

    Fig. 6 Instantaneous contour plot of vorticity amplitude at Section b: (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4

    The vorticity amplitude distribution formed in the four GSSs was analyzed. Table 2 shows the mean value of the vorticity amplitude of Line 1. Table 2 shows that whether the wheel was rotating or not hardly affected the vorticity amplitude of Line 1, while the vorticity amplitude of Line 1 on stationary ground was up to 631.3% larger than that on moving ground.

    Table 2 Mean value of the vorticity amplitude

    3.3 Q value

    As shown in Figs. 7 and 8, we can observe the influence of the four GSSs on the spatial distribution of the scale/strength of the vortex structure around the high-speed train. The four cases had no influence on the spatial distribution of the scale/strength of vortex structure in Region I, the upper part of the train. In Region II, near the nose of the head car, the four GSSs hardly affected the vortex structure scale/strength spatial distribution. In Region III, the area from the discharger of the head car to Bogie 1, the vortex structure formed in the "moving ground" condition was larger than that in the "stationary ground" condition, and the moving speed of the vortex structure was faster along both sides. Compared with the "stationary wheel" case, the "rotating wheel" case slightly increased the scale and strength of the vortex structure inside the cavity, but decreased the number of the vortexes inside the cavity and propagated a larger group of vortexes downstream. In Region IV, the area from beneath Bogie 2 to Bogie 3, the scale of the vortex structure under "moving ground" was larger than that under "stationary ground". There was also a stronger tendency of the vortex to move to both sides. In Region V, the wake area near the cowcatcher of the tail car, the vortex group structure rotated around the center in the "stationary ground" condition and formed a larger vortex structure on the upper part of the tail streamline than in the "moving ground" condition. However, whether the wheel was rotating or not hardly affected the formation/revolution of the vortex structure in this region. Note that the "moving ground" conveyed the vortex group structure formed in the wake region farther downstream than did the "stationary ground" condition.

    Fig. 7 Instantaneous contour plot of Q value at Section a: (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4

    Fig. 8 Instantaneous contour plot of Q value at Section b: (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4

    In summary, compared with "stationary ground", the "moving ground" cases increased the velocity amplitude of the flow between the train and the ballast, the scale of the vortex structure, and the vortex groups transporting downstream. Furthermore, it enhanced the tendency of the vortex group under the train to move to the sides. The "rotating wheel" increased the velocity amplitude inside the cavity, the vortex amplitude and structure scale/strength, and the scale of the vortex group conveyed downstream. Also, it enhanced the shear flow release strength of the bogie guides and decreased the amount of vortex groups inside the cavity. Note that the distributions of velocity and vorticity amplitude near the ground in the "stationary ground" case were linked to the grid scale and ground boundary layer effect. However, those phenomena did not appear in the "moving ground" condition. In addition, compared with the "moving ground" condition, the vortex structure rotating around the center in the wake mixing zone under the "stationary ground" condition was more identifiable, and a large scale of vortex structure was formed above the tail streamline.

    The reasons for the difference in the simulated ground flow field between the four GSSs are as follows. A boundary layer forms in the "stationary ground" and becomes thicker along the train, while the velocity inside the boundary layer is lower than the main flow, thereby preventing more air from flowing into the head car bottom and out of the tail car. Therefore, air accumulates at the head car upstream and the tail car downstream. The velocity gradient of the boundary layer is large resulting in the large vortex amplitude. But the vortex amplitude of the "moving ground" is small and conveys a strong vortex flow generated by the bogie cabin to the downstream. Therefore, the velocity amplitude of the upstream train bottom space under the "moving ground" is small, but the velocity amplitude of the downstream train bottom space is almost the same as that of the "stationary ground".

    Among the four GSSs, the "moving ground+rotating wheel" had the closest emulation of the true flow field, followed by the "moving ground+stationary wheel", while the other two GSSs hardly met the requirements for refined simulation of high-speed train aerodynamic noise.

    4 Aerodynamic noise source

    whereis the fluctuating pressure on the high-speed train, Pa;is the physical quantity of the Fourier transform for, Pa;is the frequency, Hz;is the fluctuating force, N;is the area of the noise source, m2; the symbol ′ represents the time derivative.

    4.1 Intensity characteristics

    The root mean square of the high-speed train surface fluctuating pressure change rate was defined as a function by

    whereis the total calculation time.

    Eq. (3) can be used to characterize the average effect of sound source intensity during the sampling time. Fig. 9 shows the dipole noise source intensity distribution cloud map of the train in each of the four GSSs.

    Fig. 9 shows that the strength distribution law of the vehicle dipole noise source was consistent in the four GSSs. For instance, the areas with strong dipole sources included the lower sides of each bogie cabin and the bottom of the headstock and the bogies. The strength of the dipole sources on Bogie 1 was larger than that on the other bogies. However, there was little difference in the magnitude of the train body dipole noise source intensity calculated by the four GSSs. To the train bottom downstream of the head car cowcatcher, whether the wheel was rotating or not hardly influenced the dipole noise source distribution in that region, and the area of strong dipole noise source distribution under the "moving ground" was significantly larger than that under the "stationary ground". With regard to the six-bogie region, the strong dipole noise source distribution area was larger when the wheel was rotating. For the five-bogie region upstream, the area of strong dipole noise source distribution under the "moving ground" was significantly larger than that under the "stationary ground". However, at the six-bogie domains, it was just the opposite. To quantitatively analyze the dipole noise sources calculated by the four GSSs, since the flow field restored by Case 1 was the most realistic, we counted the sound power of the other three GSSs based on the equivalent sound power of each component or the whole train calculated by Case 1. The component number is in Fig. S1 of the ESM.

    Fig. 9 Distribution cloud maps of on the surface of the train: (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4

    4.2 Frequency characteristics

    The effects of the four GSSs on aerodynamic noise sources were discussed above in terms of aerodynamic noise source strength. We focused only on the details of the effects of the four cases on the whole car, the bottom of the head car streamlined, the bottom of the tail car streamlined, Bogie 1, Bogie 6, and the upper part of the tail car streamlined from the sound source frequency characteristics.

    For the components discussed above and the whole car, we calculated the relative percentages of the other three GSSs of the sound power equivalent to Case 1 within the corresponding frequency range based on the equivalent sound power within every frequency band (1/3 octave) (Fig. 10).

    With regard to the equivalent sound power of the whole car, the result calculated by Case 3 at all frequency bands was almost the same as that of Case 1, showing a maximum difference value of 1.0 dB with the corresponding frequency band center of 40 Hz. At the frequency range of [125, 10000] Hz, the results of Case 2 and Case 4 were slightly smaller than that of Case 1, with a maximum difference of 1.9 dB with the corresponding frequency band center of 800 Hz. In the range of [20, 100] Hz, the calculation results of Case 2 and Case 4 were bigger than that of Case 1. The biggest difference was 3.2 dB and the corresponding frequency band center frequency was 50 Hz. Fig. 11 shows the percentage of the equivalent sound power result calculated by Case 2, Case 3, and Case 4 at 50 Hz relative to Case 1. The equivalent sound power at 50 Hz calculated by Case 2 and Case 4 apparently differed from that of Case 1 in the upper part of the middle car, the upper part of windshield 2, the upper part of the tail carriage and the tail car streamlined, the bottom of windshield 2, the bottom of the tail carriage and the tail car streamlined, and Bogies 4 to 6 and their wheel pairs at the downstream of the head car.

    As shown in Fig. 10, the equivalent sound power of the head car streamlined bottom calculated by Case 3 was almost the same as that of Case 1 in most of the frequency bands, especially in the high-frequency range, but differed greatly in several low-frequency bands. For example, in the frequency band with a center frequency of 20 Hz the difference was 3.5 dB. The results of Case 2 and Case 4 were smaller than that of Case 1 in all frequency bands and the maximum difference was about 2.4 dB with the frequency band center frequency being 12.5 Hz.

    The equivalent sound power of the tail car streamlined bottom calculated by Case 3 was almost the same as that of Case 1 in most frequency bands, but differed greatly in several low-frequency bands. For example, when the frequency band center frequency was 40 Hz the difference was up to 2.4 dB. In the frequency range where the center frequency was 12.5 Hz, the results of Case 2 and Case 4 were about 2.8 and 1.0 dB smaller than that of Case 1, respectively. At other frequency ranges, the results of Case 2 and Case 4 were bigger than that of Case 1, with a maximum of about 5.8 dB and a corresponding center frequency of 20 Hz.

    The equivalent sound power of Bogie 6 calculated by Case 3 was almost the same as that of Case 1 in most frequency bands, but the results were variable in several low-frequency bands. For instance, the difference was 1.8 dB when the center frequency was 16 Hz. In most frequency bands, the results of Case 2 and Case 4 were apparently larger than that of Case 1. The largest difference was about 5.3 dB with the corresponding center frequency being 20 Hz. The result of Case 2 was slightly smaller than that of Case 1 only when the frequency band center frequency was 12.5 Hz.

    In general, whether the wheel was rotating and whether the ground was moving affected the equivalent sound power mainly in the frequency range below 100 Hz. Whether the ground was moving had a more significant influence on the equivalent sound power than whether the wheel was rotating. In most frequency bands, the equivalent sound power of the head car downstream calculated by the "stationary ground" was larger than that of the "moving ground", while the components in the head car calculated by the "stationary ground" were smaller than those of the "moving ground". We adopted the same method to analyze the percentage of the equivalent sound source, and the results are shown in Fig. S2.

    5 Far field radiation noise

    Combined with the mirror imaging principle, the Curl acoustic integral equation can be used to calculate the far field radiation noise considering the ground effect. The calculation process is given in Section S2 of the ESM.

    To study the characteristics of high-speed train radiated noise, we set 16 measuring points every 5 m along the train. These points were 25 m from the train's central axis and 3.5 m above the ground. The first point was located in the 2.5 m-position downstream of the nose of the head car. The last point was located in the 2.5 m-position upstream of the nose of the tail car. The locations of these points are shown in Fig. 12. Considering the ground reflection, we set another 16 points which were symmetric about the ground according to the specular reflection principle.

    Fig. 12 Locations of the measuring points

    5.1 Intensity characteristics

    To study the aerodynamic noise sound pressure level (SPL), we compared the measuring point test results of the four GSSs (Fig. 13a). The difference distribution curve of whether the wheel was rotating or not is shown in Fig. 13b, and the curve of "Case 3-Case 1" represents the difference distribution curve of whether the wheel was rotating or not in the "moving ground" condition. Fig. 13c gives the difference distribution curve of whether the ground was moving or not, in which the curve marked with "Case 1-Case 2" represents the difference distribution curve in the "rotating wheel" condition and the curve "Case 3?-?Case 4" represents the difference distribution curve in the "stationary wheel" condition.

    Fig. 13 shows that the distribution of aerodynamic noise sound pressure levels calculated by the four GSSs was consistent. But the aerodynamic noise pressure level results of the four GSSs apparently differed in their simulation values. Compared with the "stationary wheel", the result of the "rotating wheel" was smaller in the upstream, while the result tested by the downstream points was the opposite. The mean difference value of whether the wheel was rotating or not at 16 measuring points was 0.3 dBA. Compared with "stationary ground", the result for the "moving ground" of each point was larger, especially the points in the upstream, while in the downstream the difference was smaller. The average difference between whether the ground was moving or not at the 16 points was 1.8 dBA. Whether the ground was moving or not had a more significant influence on the aerodynamic noise sound pressure level radiated by the whole train.

    Fig. 13 Aerodynamic noise sound pressure level radiated by the vehicle to each measuring point: (a) results of the four GSSs; (b) difference between whether the wheel pair rotated or not; (c) difference between whether the ground was moving or not. SPL represents the A-weighted sound pressure level (dBA) and the reference quantity is 2×10-5 Pa

    Fig. 14 Histogram of the difference in between the "rotating wheel" and "stationary wheel"

    Fig. 15 Histogram of the difference in between the "moving ground" and "stationary ground"

    5.2 Frequency characteristics

    In this section, we discuss the influence of the four GSSs on the characteristics of the aerodynamic noise radiated by the whole car and its components. We discuss the influence of the aerodynamic noise spectrum characteristics in detail for the whole car, the bottom of the head/tail car streamlined.

    Fig. 16 spectrum curve radiated by the whole car: (a) result of the four GSSs; (b) value of the difference between "rotating wheel" and "stationary wheel"; (c) value of the difference between "moving ground" and "stationary ground"

    According to Figs. 17 and 18, the spectrum curves of the above-mentioned components calculated by the four GSSs were identical, showing the distribution laws as small at both ends and large in the middle. However, the radiation of the components calculated by the four GSSs in each frequency band was quite different.

    Fig. 17 spectral curve radiated by the bottom of the head car streamlined: (a) result of the four GSSs; (b) value of the difference between "rotating wheel" and "stationary wheel"; (c) value of the difference between "moving ground" and "stationary ground"

    Fig. 18 spectral curve radiated by the bottom of the tail car streamlined: (a) result of the four GSSs; (b) value of the difference between "rotating wheel" and "stationary wheel"; (c) value of the difference between "moving ground" and "stationary ground"

    The "moving ground" condition increased the velocity amplitude of the flow beneath the train, the scale of the vortex, the number of vortex groups transported downstream, and the tendency of the vortex group beneath the train to move to both sides of the train. This phenomenon did not appear in the "moving ground" condition, which is why the aerodynamic noise sources of the high-speed train and the noise intensity and frequency characteristics of the far field radiation were visibly different in the "stationary ground" and "moving ground" conditions.

    6 Conclusions

    In this study, we investigated the influence of four GSSs on the flow structure of a high-speed train, the sources of noise, and the far field radiation noise. Whether the wheel rotated or not and whether the ground moved or not hardly affected the flow field above the roof of the train. They affected mainly the flow field at the bottom of the train, where they had little effect on the aero-acoustic performance of the upper parts of the train, but had a significant effect on the aero-acoustic performance of the lower parts of the train. The effect of the wheels on the equivalent sound power of the whole train was less than 5%, but the effect of the ground movement condition on the equivalent sound power of the train was more than 15%. The average level of influence of rotation or non-rotation of the wheelset on the sound pressure level of aerodynamic noise was 0.3 dBA, while the level of influence of ground movement or not was 1.8 dBA. Wheelset rotation and ground movement mainly affected the equivalent sound power of the sound source and the sound pressure level below 100 Hz.

    Considering that "wheelset rotation" increased the instability of numerical calculation, we suggest that only "ground slip" should be considered in the numerical calculation of the aerodynamic noise of high-speed trains. Considering that the ground was stationary in the acoustic wind tunnel experiment, we suggest that the experimental results of the high-speed train should be corrected appropriately in an acoustic wind tunnel.

    This work is supported by the National Natural Science Foundation of China (No. 52272363) and the Foundation of the Key Laboratory of Aerodynamic Noise Control (No. ANCL20200302), China.

    Xiaoming TAN designed the research. Xiaohong ZHANG and Zhigang YANG processed the corresponding data. Xiaoming TAN wrote the first draft of the manuscript. Linli GONG helped to organize the manuscript. Linli GONG revised and edited the final version.

    Xiaoming TAN, Linli GONG, Xiaohong ZHANG, and Zhigang YANG declare that they have no conflict of interest.

    Kouroussis G, Zhu SY, Vogiatzis K, 2021. Noise and vibration from transportation., 22(1):1-5. https://doi.org/10.1631/jzus.A20NVT01

    Kwon HB, Park YW, Lee DH, et al., 2001. Wind tunnel experiments on Korean high-speed trains using various ground simulation techniques., 89(13):1179-1195. https://doi.org/10.1016/S0167-6105(01)00107-6

    Liang XF, Liu HF, Dong TY, et al., 2020. Aerodynamic noise characteristics of high-speed train foremost bogie section., 27(6):1802-1813. https://doi.org/10.1007/s11771-020-4409-8

    Liu JL, Zhang JY, Zhang WH, 2013. Study of computational method of far-field aerodynamic noise of a high-speed train considering ground effect., 30(1):94-100 (in Chinese). https://doi.org/10.7511/jslx201301016

    Liu W, Guo D, Zhang Z, et al., 2019. Effects of bogies on the wake flow of a high-speed train., 9(4):759. https://doi.org/10.3390/app9040759

    Paz C, Suárez E, Gil C, 2017. Numerical methodology for evaluating the effect of sleepers in the underbody flow of a high-speed train., 167:140-147. https://doi.org/10.1016/j.jweia.2017.04.017

    Wang SB, Burton D, Herbst A, et al., 2018. The effect of bogies on high-speed train slipstream and wake., 83:471-489. https://doi.org/10.1016/j.jfluidstructs.2018.03.013

    Tan XM, Liu HF, Yang ZG, et al., 2018. Characteristics and mechanism analysis of aerodynamic noise sources for high-speed train in tunnel., 2018:5858415. https://doi.org/10.1155/2018/5858415

    Tyll JS, Liu D, Schetz JA, et al., 1996. Experimental studies of magnetic levitation train aerodynamics., 34(12):2465-2470. https://doi.org/10.2514/3.13425

    Xia C, Shan XZ, Yang ZG, 2016. Influence of ground configurations in wind tunnels on the slipstream of a high-speed train. The 8th International Colloquium on Bluff Body Aerodynamics and Applications.

    Xia C, Shan XZ, Yang ZG, 2017. Comparison of different ground simulation systems on the flow around a high-speed train., 231(2):135-147. https://doi.org/10.1177/0954409715626191

    Yi SH, Zou JJ, Wu GF, et al., 1997. Experimental investigation for ground effects of the high speed train models on a plate with uniform boundary layer suction., 11(2):95-100 (in Chinese).

    Zhang J, Li JJ, Tian HQ, et al., 2016. Impact of ground and wheel boundary conditions on numerical simulation of the high-speed train aerodynamic performance., 61:249-261. https://doi.org/10.1016/j.jfluidstructs.2015.10.006

    Zhu JY, Hu ZW, Thompson DJ, 2017. The effect of a moving ground on the flow and aerodynamic noise behaviour of a simplified high-speed train bogie., 5(2):110-125. https://doi.org/10.1080/23248378.2016.1212677

    Electronic supplementary materials

    Sections S1–S3, Figs. S1–S5, and Eqs. (S1)–(S3)

    地面效應(yīng)對高速列車流場結(jié)構(gòu)及氣動(dòng)噪聲的影響

    譚曉明1,2,龔林立1,張曉紅1,楊志剛2

    1湖南理工學(xué)院,機(jī)械工業(yè)激光磨削復(fù)合智能制造與服役性能優(yōu)化重點(diǎn)實(shí)驗(yàn)室,中國岳陽,414000;2中南大學(xué),軌道交通安全教育部重點(diǎn)實(shí)驗(yàn)室,中國長沙,410075

    高速列車作為高速地面交通工具,不可避免地會遇到地面效應(yīng)問題。地面效應(yīng)模擬一直是高速列車風(fēng)洞試驗(yàn)的技術(shù)難點(diǎn)。地面效應(yīng)現(xiàn)象的準(zhǔn)確模擬對高速列車空氣動(dòng)力學(xué)和氣動(dòng)噪聲的預(yù)測精度有很大的影響。通過對比4種地面模擬系統(tǒng)(GSS)的流聲場結(jié)果,研究不同GSS對流場結(jié)構(gòu)、氣動(dòng)聲源和遠(yuǎn)場輻射噪聲特性的影響規(guī)律,為高速列車聲學(xué)風(fēng)洞試驗(yàn)提供指導(dǎo)。

    1. 搭建高速列車地面模擬系統(tǒng),模擬不同邊界條件;2. 明確輪對旋轉(zhuǎn)與地面滑移對高速列車氣動(dòng)噪聲幅值的相對增量及影響頻率范圍。

    1. 在仿真系統(tǒng)中建立“移動(dòng)地面+旋轉(zhuǎn)輪對”、“靜止地面+旋轉(zhuǎn)輪對”、“移動(dòng)地面+靜止輪對”和“靜止地面+靜止輪對”四種地面模擬系統(tǒng);2. 采用大渦模擬和旋度聲學(xué)積分方程,對高速列車的流聲場結(jié)果進(jìn)行模擬;3. 通過對比4種GSS的流聲場結(jié)果,研究不同GSS對流場結(jié)構(gòu)、氣動(dòng)聲源和遠(yuǎn)場輻射噪聲特性的影響規(guī)律。

    1. 移動(dòng)地面和旋轉(zhuǎn)輪對是影響列車底部氣動(dòng)聲學(xué)性能的主要因素;2. 旋轉(zhuǎn)輪對對整車等效聲源功率的影響不大于5%,且移動(dòng)地面對整車等效聲源功率的影響大于15%;3. 旋轉(zhuǎn)輪對對整車輻射聲壓級的平均影響為0.3 dBA,且運(yùn)動(dòng)地面對整車輻射聲壓級的平均影響為1.8 dBA;它們主要影響100 Hz以下的氣動(dòng)聲學(xué)性能。

    高速列車;氣動(dòng)聲學(xué);流場結(jié)構(gòu);大渦模擬;移動(dòng)地面;旋轉(zhuǎn)輪對

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

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

    Zhigang YANG, yangzg1976@163.com

    Zhigang YANG, https://orcid.org/0000-0001-7255-7324

    Received Jan. 16, 2023;

    Revision accepted May 27, 2023;

    Crosschecked Jan. 5, 2024

    ? Zhejiang University Press 2024

    猜你喜歡
    影響
    美食網(wǎng)紅如何影響我們吃什么
    英語文摘(2022年4期)2022-06-05 07:45:18
    是什么影響了滑動(dòng)摩擦力的大小
    哪些顧慮影響擔(dān)當(dāng)?
    影響大師
    沒錯(cuò),痛經(jīng)有時(shí)也會影響懷孕
    媽媽寶寶(2017年3期)2017-02-21 01:22:28
    擴(kuò)鏈劑聯(lián)用對PETG擴(kuò)鏈反應(yīng)與流變性能的影響
    中國塑料(2016年3期)2016-06-15 20:30:00
    基于Simulink的跟蹤干擾對跳頻通信的影響
    如何影響他人
    APRIL siRNA對SW480裸鼠移植瘤的影響
    啦啦啦在线观看免费高清www| 亚洲第一av免费看| 国产精品嫩草影院av在线观看| 亚洲欧洲精品一区二区精品久久久 | 国产在线免费精品| 欧美日韩亚洲高清精品| 久久精品久久精品一区二区三区| 国产精品一国产av| 成年美女黄网站色视频大全免费| 一区福利在线观看| 香蕉精品网在线| 国产成人a∨麻豆精品| 丰满迷人的少妇在线观看| 国产精品成人在线| 老司机影院毛片| 精品国产一区二区三区四区第35| 晚上一个人看的免费电影| 中文乱码字字幕精品一区二区三区| 天美传媒精品一区二区| 日日撸夜夜添| 亚洲天堂av无毛| 午夜福利在线观看免费完整高清在| 国产xxxxx性猛交| 少妇人妻精品综合一区二区| 欧美中文综合在线视频| 久久久精品国产亚洲av高清涩受| 电影成人av| 亚洲一区二区三区欧美精品| 久久人妻熟女aⅴ| 在线观看人妻少妇| 在线免费观看不下载黄p国产| 高清欧美精品videossex| 91久久精品国产一区二区三区| 哪个播放器可以免费观看大片| 纵有疾风起免费观看全集完整版| 老汉色av国产亚洲站长工具| 另类亚洲欧美激情| 十八禁高潮呻吟视频| 另类精品久久| 2021少妇久久久久久久久久久| 日韩av免费高清视频| 精品国产乱码久久久久久男人| 宅男免费午夜| 久久久精品免费免费高清| 纯流量卡能插随身wifi吗| 曰老女人黄片| av国产久精品久网站免费入址| 日本爱情动作片www.在线观看| 97在线视频观看| 色婷婷av一区二区三区视频| 久久久久精品久久久久真实原创| 母亲3免费完整高清在线观看 | 国产成人精品久久二区二区91 | 一区二区三区激情视频| 午夜精品国产一区二区电影| 99久国产av精品国产电影| 久久精品国产综合久久久| 日本爱情动作片www.在线观看| 一级毛片电影观看| 亚洲天堂av无毛| 在线天堂最新版资源| 美女高潮到喷水免费观看| 欧美激情 高清一区二区三区| 精品亚洲乱码少妇综合久久| 性色av一级| 成年人午夜在线观看视频| av国产久精品久网站免费入址| 国产精品一区二区在线不卡| 亚洲美女视频黄频| 亚洲人成电影观看| 制服丝袜香蕉在线| 91精品三级在线观看| 黄色配什么色好看| 久久av网站| 国产av国产精品国产| 精品少妇内射三级| 中文字幕精品免费在线观看视频| 午夜日韩欧美国产| 少妇 在线观看| 啦啦啦啦在线视频资源| 欧美 日韩 精品 国产| 久久这里有精品视频免费| 精品国产一区二区三区久久久樱花| 欧美人与性动交α欧美精品济南到 | 亚洲美女黄色视频免费看| 国产熟女午夜一区二区三区| 亚洲欧美成人综合另类久久久| 搡老乐熟女国产| 91在线精品国自产拍蜜月| 我要看黄色一级片免费的| 建设人人有责人人尽责人人享有的| 丝瓜视频免费看黄片| 狠狠精品人妻久久久久久综合| 三级国产精品片| 热re99久久精品国产66热6| 成年女人在线观看亚洲视频| 成人亚洲欧美一区二区av| 久久精品国产自在天天线| 男女边吃奶边做爰视频| 亚洲精品国产av成人精品| 国产成人欧美| 宅男免费午夜| 精品福利永久在线观看| 免费不卡的大黄色大毛片视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 亚洲色图 男人天堂 中文字幕| 国产精品免费大片| 国产精品久久久久久av不卡| 亚洲国产欧美日韩在线播放| 欧美精品高潮呻吟av久久| 免费大片黄手机在线观看| 国产成人精品久久久久久| 亚洲色图 男人天堂 中文字幕| 久久久久国产一级毛片高清牌| 久久午夜福利片| 曰老女人黄片| av在线app专区| 精品国产一区二区三区四区第35| 日本欧美视频一区| 十八禁网站网址无遮挡| 欧美人与性动交α欧美精品济南到 | av卡一久久| 九九爱精品视频在线观看| 国产免费现黄频在线看| 亚洲伊人久久精品综合| 久久国产亚洲av麻豆专区| 国产一级毛片在线| 亚洲精华国产精华液的使用体验| 国产不卡av网站在线观看| 99re6热这里在线精品视频| 亚洲视频免费观看视频| 不卡视频在线观看欧美| 成人毛片a级毛片在线播放| av福利片在线| 久久精品国产自在天天线| 欧美老熟妇乱子伦牲交| 精品一品国产午夜福利视频| 欧美激情 高清一区二区三区| 成人毛片60女人毛片免费| 亚洲欧美成人综合另类久久久| 久久ye,这里只有精品| 久久午夜福利片| 97在线人人人人妻| 丁香六月天网| 国产成人精品在线电影| 国产极品粉嫩免费观看在线| 成人漫画全彩无遮挡| 国产精品国产三级国产专区5o| 黄色视频在线播放观看不卡| 丝袜美腿诱惑在线| 国产成人精品无人区| 国产 一区精品| 日本av手机在线免费观看| 久久久久网色| 最近2019中文字幕mv第一页| 欧美日韩亚洲高清精品| 欧美成人午夜免费资源| 91久久精品国产一区二区三区| 哪个播放器可以免费观看大片| 国产精品亚洲av一区麻豆 | 最黄视频免费看| 一级毛片 在线播放| 中国三级夫妇交换| 欧美日韩视频高清一区二区三区二| av卡一久久| 最近2019中文字幕mv第一页| 亚洲国产毛片av蜜桃av| 亚洲,欧美精品.| 美女大奶头黄色视频| 中国三级夫妇交换| 校园人妻丝袜中文字幕| 99热全是精品| 久久婷婷青草| 国产成人精品久久久久久| 日韩中文字幕视频在线看片| 国产在视频线精品| 九草在线视频观看| 国产精品无大码| 久久久精品94久久精品| 亚洲国产精品成人久久小说| 日本欧美视频一区| 高清黄色对白视频在线免费看| 在线观看三级黄色| 亚洲欧美清纯卡通| 亚洲精品一区蜜桃| 亚洲国产精品999| 午夜福利在线观看免费完整高清在| 制服诱惑二区| 黑人猛操日本美女一级片| 我要看黄色一级片免费的| 伊人久久大香线蕉亚洲五| 亚洲国产日韩一区二区| 亚洲精品视频女| av又黄又爽大尺度在线免费看| 91在线精品国自产拍蜜月| 热re99久久精品国产66热6| 人妻 亚洲 视频| 在线精品无人区一区二区三| 只有这里有精品99| 国产片内射在线| 如何舔出高潮| 亚洲激情五月婷婷啪啪| 亚洲婷婷狠狠爱综合网| 91精品三级在线观看| 日韩中文字幕视频在线看片| 国产深夜福利视频在线观看| 看十八女毛片水多多多| 欧美精品av麻豆av| 免费黄频网站在线观看国产| 麻豆av在线久日| 在线观看美女被高潮喷水网站| 麻豆精品久久久久久蜜桃| 免费在线观看视频国产中文字幕亚洲 | 18在线观看网站| 精品第一国产精品| 亚洲精品中文字幕在线视频| 国产一区二区激情短视频 | 中国三级夫妇交换| 2018国产大陆天天弄谢| 天美传媒精品一区二区| 777久久人妻少妇嫩草av网站| 九九爱精品视频在线观看| a级片在线免费高清观看视频| 成年美女黄网站色视频大全免费| 国产精品秋霞免费鲁丝片| 我的亚洲天堂| 久久精品久久久久久久性| 日韩电影二区| 亚洲在久久综合| 视频在线观看一区二区三区| 多毛熟女@视频| 一本久久精品| 国产精品麻豆人妻色哟哟久久| 亚洲人成77777在线视频| 精品一区二区免费观看| 亚洲国产最新在线播放| 亚洲av综合色区一区| 深夜精品福利| 午夜av观看不卡| 在线观看美女被高潮喷水网站| 婷婷成人精品国产| 女人久久www免费人成看片| 欧美日韩视频精品一区| 美国免费a级毛片| 欧美+日韩+精品| 在线观看三级黄色| 亚洲图色成人| 婷婷色综合大香蕉| 最近2019中文字幕mv第一页| 满18在线观看网站| 美女福利国产在线| 十八禁网站网址无遮挡| 色哟哟·www| 极品少妇高潮喷水抽搐| 青草久久国产| 狂野欧美激情性bbbbbb| 国产精品一区二区在线观看99| 国产精品欧美亚洲77777| 亚洲精品国产av蜜桃| 纵有疾风起免费观看全集完整版| 99九九在线精品视频| 国产精品二区激情视频| 多毛熟女@视频| 亚洲精品中文字幕在线视频| 免费高清在线观看视频在线观看| 性色avwww在线观看| 国产精品二区激情视频| 亚洲国产精品一区三区| 精品一品国产午夜福利视频| 国产精品亚洲av一区麻豆 | www.av在线官网国产| www.熟女人妻精品国产| 日韩精品免费视频一区二区三区| 女人精品久久久久毛片| 欧美亚洲日本最大视频资源| 国产极品天堂在线| 亚洲欧美精品综合一区二区三区 | 看免费av毛片| 亚洲成av片中文字幕在线观看 | 日本91视频免费播放| 一区在线观看完整版| 亚洲欧美清纯卡通| 国产精品av久久久久免费| 久久97久久精品| 成年动漫av网址| 美女国产视频在线观看| 综合色丁香网| 深夜精品福利| 亚洲情色 制服丝袜| 精品一品国产午夜福利视频| 国产精品秋霞免费鲁丝片| 少妇 在线观看| 国产成人精品在线电影| 高清不卡的av网站| 亚洲第一青青草原| 成人影院久久| 在线观看人妻少妇| 中文欧美无线码| 亚洲精品日本国产第一区| 久久精品国产亚洲av涩爱| 中文精品一卡2卡3卡4更新| 久久久久久人妻| 满18在线观看网站| videossex国产| 人妻一区二区av| av又黄又爽大尺度在线免费看| 美女视频免费永久观看网站| 18在线观看网站| 下体分泌物呈黄色| 亚洲欧美精品自产自拍| 毛片一级片免费看久久久久| 久久久国产精品麻豆| xxxhd国产人妻xxx| 日本黄色日本黄色录像| 久久热在线av| 亚洲国产精品999| 欧美日韩亚洲国产一区二区在线观看 | 99久久中文字幕三级久久日本| 欧美日韩亚洲高清精品| 满18在线观看网站| 亚洲欧美一区二区三区久久| 国产成人精品在线电影| 美女福利国产在线| 中文字幕人妻丝袜一区二区 | 不卡av一区二区三区| 亚洲成色77777| 欧美bdsm另类| 狠狠精品人妻久久久久久综合| av在线观看视频网站免费| 国产片特级美女逼逼视频| 色网站视频免费| 91久久精品国产一区二区三区| 天天影视国产精品| 久久久国产一区二区| av线在线观看网站| 久久99蜜桃精品久久| 国产黄色视频一区二区在线观看| 精品国产一区二区三区久久久樱花| 日日摸夜夜添夜夜爱| 国产在线视频一区二区| 青春草亚洲视频在线观看| 欧美在线黄色| 亚洲av福利一区| 欧美成人午夜免费资源| 成人毛片a级毛片在线播放| 99九九在线精品视频| 黄色配什么色好看| 毛片一级片免费看久久久久| 久久久久久久久久人人人人人人| 国产亚洲最大av| 亚洲国产看品久久| 亚洲欧洲日产国产| 亚洲国产看品久久| 日本猛色少妇xxxxx猛交久久| 成年人午夜在线观看视频| 国产熟女欧美一区二区| 咕卡用的链子| 国产av码专区亚洲av| 久久久久久久久免费视频了| 人人妻人人爽人人添夜夜欢视频| 最近中文字幕高清免费大全6| 熟女av电影| 69精品国产乱码久久久| 婷婷色av中文字幕| 中文字幕人妻熟女乱码| 天天躁夜夜躁狠狠躁躁| 精品一区二区免费观看| 波多野结衣av一区二区av| 中文字幕人妻熟女乱码| 欧美成人精品欧美一级黄| 婷婷色综合大香蕉| 老熟女久久久| 街头女战士在线观看网站| 亚洲欧美一区二区三区国产| 国产精品 国内视频| 丝袜脚勾引网站| 寂寞人妻少妇视频99o| 99久国产av精品国产电影| 久久精品国产亚洲av高清一级| 2021少妇久久久久久久久久久| 成年人免费黄色播放视频| 一区二区三区四区激情视频| 国产欧美亚洲国产| 成人18禁高潮啪啪吃奶动态图| 人妻 亚洲 视频| av又黄又爽大尺度在线免费看| 午夜福利视频在线观看免费| 777米奇影视久久| 久久久国产欧美日韩av| xxx大片免费视频| 天天躁夜夜躁狠狠躁躁| 女人精品久久久久毛片| 美女中出高潮动态图| 亚洲一级一片aⅴ在线观看| 建设人人有责人人尽责人人享有的| 在线天堂最新版资源| 99久久精品国产国产毛片| 国产精品一区二区在线观看99| 18禁裸乳无遮挡动漫免费视频| 一边摸一边做爽爽视频免费| 成人黄色视频免费在线看| 一二三四在线观看免费中文在| 久久久亚洲精品成人影院| 国产精品久久久久成人av| 欧美激情高清一区二区三区 | 久久久a久久爽久久v久久| 久久久久国产一级毛片高清牌| 青春草视频在线免费观看| 97在线视频观看| 国产又色又爽无遮挡免| 国产精品欧美亚洲77777| 秋霞伦理黄片| 日本猛色少妇xxxxx猛交久久| 日本欧美视频一区| 午夜久久久在线观看| 18+在线观看网站| 黑人欧美特级aaaaaa片| 女人久久www免费人成看片| 在线观看免费高清a一片| 国产在线免费精品| 亚洲国产欧美网| av.在线天堂| av女优亚洲男人天堂| 视频在线观看一区二区三区| 国产精品久久久久久久久免| 亚洲综合色惰| 亚洲精品久久成人aⅴ小说| 男女国产视频网站| 人妻一区二区av| 色视频在线一区二区三区| 欧美+日韩+精品| 黄色毛片三级朝国网站| 免费观看av网站的网址| 男人爽女人下面视频在线观看| 国产成人精品无人区| 亚洲五月色婷婷综合| 青草久久国产| 天天躁夜夜躁狠狠躁躁| av又黄又爽大尺度在线免费看| 国产色婷婷99| 精品第一国产精品| 久久精品人人爽人人爽视色| av不卡在线播放| 国产成人精品无人区| 亚洲欧美日韩另类电影网站| 日本色播在线视频| 制服人妻中文乱码| 久久精品国产亚洲av涩爱| 黄片播放在线免费| 国产精品.久久久| 春色校园在线视频观看| 亚洲av中文av极速乱| 国产精品久久久久久久久免| 99久久中文字幕三级久久日本| 黑人巨大精品欧美一区二区蜜桃| 国产精品人妻久久久影院| 少妇的逼水好多| 亚洲精品久久午夜乱码| 伦理电影大哥的女人| 春色校园在线视频观看| 18在线观看网站| 在线亚洲精品国产二区图片欧美| 亚洲一级一片aⅴ在线观看| 日韩欧美一区视频在线观看| 我的亚洲天堂| 在线天堂最新版资源| 国产黄色视频一区二区在线观看| 国产精品国产三级专区第一集| 国产黄色视频一区二区在线观看| 黄片无遮挡物在线观看| 极品少妇高潮喷水抽搐| 亚洲精品国产一区二区精华液| 熟女av电影| 亚洲色图综合在线观看| 99国产综合亚洲精品| 男人爽女人下面视频在线观看| 亚洲精品视频女| 啦啦啦啦在线视频资源| 99热网站在线观看| 国产午夜精品一二区理论片| 久久国产亚洲av麻豆专区| 免费高清在线观看日韩| 日韩 亚洲 欧美在线| 五月伊人婷婷丁香| 日本欧美视频一区| 在线观看国产h片| 99国产精品免费福利视频| 国产精品久久久久久av不卡| 欧美精品一区二区大全| 欧美日韩亚洲高清精品| av免费在线看不卡| 日韩在线高清观看一区二区三区| 国产成人a∨麻豆精品| 午夜福利网站1000一区二区三区| 亚洲精品成人av观看孕妇| 制服诱惑二区| 国产精品麻豆人妻色哟哟久久| 男女边吃奶边做爰视频| av有码第一页| 人妻少妇偷人精品九色| 日韩中字成人| 男人添女人高潮全过程视频| 国产日韩一区二区三区精品不卡| 黑丝袜美女国产一区| 男人爽女人下面视频在线观看| 最新的欧美精品一区二区| 九九爱精品视频在线观看| 日韩av免费高清视频| 国产毛片在线视频| 好男人视频免费观看在线| 美女大奶头黄色视频| 18禁动态无遮挡网站| 最近最新中文字幕免费大全7| 日韩一区二区视频免费看| 亚洲 欧美一区二区三区| 久久人妻熟女aⅴ| 侵犯人妻中文字幕一二三四区| 少妇被粗大猛烈的视频| 精品国产一区二区三区四区第35| 久久精品久久精品一区二区三区| 丝袜美腿诱惑在线| 韩国高清视频一区二区三区| 亚洲国产av新网站| 亚洲精品久久成人aⅴ小说| 日日爽夜夜爽网站| av网站免费在线观看视频| 免费看av在线观看网站| 精品午夜福利在线看| 久久久久精品人妻al黑| 女性被躁到高潮视频| 在线观看免费日韩欧美大片| 国产日韩欧美亚洲二区| 国产片内射在线| 夫妻性生交免费视频一级片| 久久人人爽人人片av| 亚洲综合色惰| 国产伦理片在线播放av一区| 亚洲国产欧美网| 亚洲av国产av综合av卡| 国产精品久久久久久精品古装| 日韩一卡2卡3卡4卡2021年| 免费看不卡的av| 久久 成人 亚洲| 最近最新中文字幕大全免费视频 | 国产片特级美女逼逼视频| 婷婷色综合大香蕉| 日韩一区二区视频免费看| 国产成人a∨麻豆精品| 男女边摸边吃奶| 丝袜在线中文字幕| 中国三级夫妇交换| 久久亚洲国产成人精品v| 亚洲综合精品二区| 天堂中文最新版在线下载| 国产精品一二三区在线看| 亚洲欧洲国产日韩| 亚洲视频免费观看视频| 日韩免费高清中文字幕av| 1024视频免费在线观看| 伊人久久大香线蕉亚洲五| 国产片内射在线| 精品国产露脸久久av麻豆| 99热全是精品| 少妇的丰满在线观看| 国产探花极品一区二区| 老司机影院成人| 久久精品久久久久久噜噜老黄| 美女中出高潮动态图| 最近中文字幕高清免费大全6| 中文天堂在线官网| 晚上一个人看的免费电影| 亚洲欧洲日产国产| 欧美日韩综合久久久久久| 啦啦啦啦在线视频资源| 大片电影免费在线观看免费| 日韩,欧美,国产一区二区三区| 欧美日韩国产mv在线观看视频| 免费人妻精品一区二区三区视频| av线在线观看网站| 夫妻午夜视频| 国产精品一区二区在线观看99| 最近的中文字幕免费完整| 一级毛片 在线播放| 精品人妻在线不人妻| 咕卡用的链子| 制服诱惑二区| 久久这里有精品视频免费| 咕卡用的链子| 国产一区二区三区av在线| 熟妇人妻不卡中文字幕| 欧美国产精品一级二级三级| 叶爱在线成人免费视频播放| 丝袜在线中文字幕| 亚洲内射少妇av| 大陆偷拍与自拍| 少妇熟女欧美另类| 久久久久国产一级毛片高清牌| 伦理电影大哥的女人| 亚洲三区欧美一区| av网站在线播放免费| 国产男女超爽视频在线观看| 韩国av在线不卡| 色吧在线观看| 黄色一级大片看看| 国产精品免费视频内射| 女人精品久久久久毛片| 成年人免费黄色播放视频| 亚洲精品第二区| 精品一区二区免费观看| 91精品三级在线观看| 亚洲欧美精品综合一区二区三区 | 午夜精品国产一区二区电影| 亚洲av福利一区| 女人久久www免费人成看片| 免费黄色在线免费观看|