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

    Transfer learning: A new aerodynamic force identification network based on adaptive EMD and soft thresholding in hypersonic wind tunnel

    2023-09-02 10:18:20YiSUNShichoLIHongliGAOXioqingZHANGJinzhouLVWeixiongLIUYingchunWU
    CHINESE JOURNAL OF AERONAUTICS 2023年8期

    Yi SUN, Shicho LI,*, Hongli GAO, Xioqing ZHANG, Jinzhou LV,Weixiong LIU, Yingchun WU

    a School of Mechanical Engineering, Southwest Jiaotong University, Chengdu 610031, China

    b China Aerodynamics Research and Development Center, Mianyang 621000, China

    KEYWORDS Aerodynamic intelligent identification model;Transfer learning;Force measurement system;Residual attention block with soft threshold;Dense block with adaptive EMD

    Abstract The aerodynamic test in the pulse combustion wind tunnel is very important for the design, evaluation and optimization of aerodynamic characteristics of the hypersonic aircraft.The test accuracy even affects the success or failure of hypersonic aircraft development.In the aerodynamic test of pulse combustion wind tunnel, the aerodynamic signal is disturbed by the inertial force signal, which seriously affects the test accuracy of aerodynamic force.Aiming at the above problems, this paper innovatively proposes an aerodynamic intelligent identification method, that is the transfer learning network based on adaptive Empirical Modal Decomposition (EMD) and Soft Thresholding (TLN-AE&ST).Compared with the existing aerodynamic intelligent identification model based on deep learning technology, this study introduces the transfer learning idea into the aerodynamic intelligent identification model for the first time.The TLN-AE&ST effectively alleviates the problem of scarcity of training samples for intelligent models due to the high cost of wind tunnel tests,and provides a new idea for further implementation of deep learning technology in the field of wind tunnel aerodynamic testing.And this study designed residual attention block with soft threshold and dense block with adaptive EMD in TLN-AE&ST model.Residual attention block with soft threshold module can more effectively suppress the influence of instrument noise signal on model training effect.Dense block with adaptive EMD makes the deep learning model no longer a black box to a certain extent,and has certain physical significance.Finally,a series of wind tunnel tests were carried out in the φ = 2.4 m pulse combustion wind tunnel of China Aerodynamic Research and Development Center to verify the effectiveness of TLN-AE&ST.

    1.Introduction

    Hypersonic wind tunnel aerodynamic test is very important for the design and optimization of aerodynamic shape of a Hypersonic Aircraft(HA),and even determines the success or failure of HA development.At present, the pulse combustion wind tunnel is a mainstream hypersonic wind tunnel used to generate high-temperature and high-speed airflow to simulate airflow under hypersonic flight conditions.The high-speed airflow (greater than 5 times the speed of sound) produces a transient impact on the aircraft test model, which causes the transient vibration of the aircraft test model.However, due to the short effective test time of the pulse combustion wind tunnel (only 500 ms), the vibration of the aircraft test model is difficult to attenuate rapidly in such a short time.Thus,the inertial force is introduced into the output of the wind tunnel Force Measurement System(FMS),which causes the aerodynamic force to be overwhelmed by inertial forces.1–3At the same time, instrument noise is introduced to the output signals.The above phenomenon resulting in unable to directly distinguish the real dynamic aerodynamic forces regularity,which seriously affects the accuracy of aerodynamic test.The structure of FMS in this paper and the noise source of the signal are shown in the top part of Fig.1.

    Moreover, as the aircraft test model gradually develops to the full-scale (the full-size model is 7 m-8 m in length and several tons in weight, compared with the scaled-size model less than 2 m in length and 1ton in weight), its stiffness further decreases,which will further reduce the identification accuracy of the FMS output signals.

    Many researchers in the field of hypersonic wind tunnel testing have carried out many valuable studies on the interference of inertial force on aerodynamic force,which can be summarized into two aspects, namely develop new FMS and develop effective aerodynamic load identification technology.For developing new FMS, such as the accelerometer balance,4–6the stress-wave balance,7–9the compensation balance10–13and the pulse strain balance14–16have been proposed in recent years.Although the above balance can effectively reduce the proportion of inertial force component in FMS output signals, it cannot fundamentally eliminate the influence of inertial force on aerodynamic force.For the aerodynamic load identification technology, researchers try to extract the hidden aerodynamic force from the oscillation signal output by FMS.Therefore, time–frequency transform methods such as filtering and Fast Fourier Transform (FFT)are introduced into FMS.These methods can effectively eliminate fixed frequency and noise interference, but also lead to some loss of the original signal component17near the fixed frequency component is inevitably eliminated.And the identification accuracy of FFT depends on the effective test time of wind tunnel.Because the effective test time of pulse combustion wind tunnel is only 200 ms, the highest frequency resolution of FFT is only 5 Hz, which will affect the identification accuracy of FFT.Therefore, FFT is often used in conventional wind tunnels with test time greater than 1 s.Wavelet transform can achieve 1 Hz frequency resolution, but its identification accuracy depends on the choice of a wavelet basis function.In addition, wavelet transform has Heisenberg uncertainty defect, that is, frequency-domain and time-domain cannot achieve high resolution at the same time.Compared with the above methods, EMD is adaptive, but it is prone to modal aliasing, which will reduce the identification accuracy of aerodynamic force, as shown in the bottom right part of Fig.1.18

    To solve the shortcomings, experts have tried to introduce Deep Learning (DL) techniques into the field of aerodynamic identification.19,20Deep learning, like dynamic calibration of FMS is to obtain the input-to-output model of FMS under dynamic force.Therefore, DL based aerodynamic intelligent identification techniques can eliminate the interference of inertial force on aerodynamic identification and has higher accuracy, as shown in the bottom left part of Fig.1.Compared to traditional aerodynamic identification techniques, aerodynamic intelligent identification techniques based on DL have the following advantages:

    (1) This kind of method does not need manual participation in the aerodynamic identification process.Therefore,the influence of professional knowledge on identification effect is reduced.Moreover, the identification accuracy is not affected by the test time of wind tunnel and frequency resolution.

    (2) The aerodynamic intelligent identification method based on deep learning consists of stacking multiple convolutional layers.Each channel of the convolution layer can be regarded as a trainable band-pass filter, which can suppress the specific noise signal in the FMS output signals adaptively.In contrast, traditional aerodynamic identification methods process the signal by manually selecting a filtering method.This mode has no universality for signals of different aircrafts and different FMSs,and needs to be adjusted according to the frequency of noise.However,the accuracy of the aerodynamic intelligent recognition model based on DL depends on the sample size used to train the model.The larger the sample size, the higher the identification accuracy.Moreover, the large-scale aircraft test model needs to be tested in the larger hypersonic wind tunnel.Highquality and sufficient data at the same Mach number are difficult to collect due to the excessively expensive single test cost.Therefore, when training the aerodynamic intelligent identification model, the scarcity of sample size is the key to restricting the identification accuracy of the aerodynamic intelligent identification model.This is also the direct reason for the low identification effect of intelligent methods based on Convolution Neural Network (CNN),20–23Recurrent Neural Network (RNN)24or Deep Belief Network (DBN).25

    Fig.1 Noise sources and pre-processing of signal.The VMD’s calculation process and parameters in the lower right corner refer to the paper.18

    Based on the above challenges, this paper proposes a new aerodynamic force identification method,transfer learning network with adaptive Empirical Mode Decomposition (EMD)and Soft Thresholding (TLN-AE&ST), from the engineering point of view.Firstly, the dynamic calibration test of FMS(i.e.applying the known excitation force and measuring the output of FMS) is carried out to obtain a large amount of training data, and then the pre-training of the model is completed.Subsequently, a small number of wind tunnel aerodynamic test samples were used to carry out the real aerodynamic test,and the feature distribution was transferred,so that the model can accurately identify the aerodynamic force.The contributions of this paper are summarized as follows:

    (1) A new intelligent aerodynamic force identification model based on deep learning is proposed in this paper.The model is used to identify and filter out the inertial force signal and instrument noise signal from the FMS output signal in the wind tunnel test, so as to obtain‘‘pure”aerodynamic force.The model will greatly reduce the design difficulty of strain gauge balance with high precision, low cost and mature development, and improve the accuracy of pulse combustion wind tunnel aerodynamic force test.

    (2) Compared with the existing aerodynamic intelligent identification model based on deep learning technology,this study introduces the transfer learning idea into the aerodynamic intelligent identification model for the first time.The TLN-AE&ST effectively alleviates the problem of scarcity of training samples for intelligent models due to the high cost of wind tunnel tests.

    (3) To improve the identification accuracy of the aerodynamic force intelligent identification model, the residual attention block with soft threshold is introduced into the aerodynamic force intelligent identification model to filter the instrument noise signal, and a Dense block with adaptive EMD is innovatively proposed, which greatly filters the inertial force noise introduced by the aircraft test model.

    The subsequent sections of this paper are structured as follows.Section 2,the principle of the related work is introduced.The details of the architecture proposed in this paper are described in Section 3.Section 4, the experimental set-up is described,and the performance of our architecture is analyzed.We close the paper with conclusions in Section 5.

    2.Principle

    2.1.Analysis of FMS output signals

    Fig.2 Test stage of FMS and components of its output signals.

    The upper left quarter of Fig.2 shows an aerodynamic test in the φ = 2.4 m impulse combustion wind tunnel with a Mach number of 6.5.The wind tunnel is a large hypersonic wind tunnel of the China Aerodynamics Research and Development Center (CARDC).The blue dotted line is the identified lift force;the red dotted line is the total pressure signal in the wind tunnel test section.Usually, a pulse-combustion wind tunnel aerodynamic test consists of three main stages: preparation stage, start-up stage, and effective experimental stage, and the signal components in the FMS output signal are different in different stages.Therefore, this section analyzes the signal components of the FMS output signals at different test stage,clarifies the characteristics of aerodynamic force signal and noise signal components, and provides theoretical support for gradually removing the noise signal component.Then,according to the characteristics of each noise signal component, the matching processing method is proposed.The components of the FMS output signals and their characteristics are as follows:

    (1) Aerodynamic force.As shown in the upper left quarter of Fig.2.In the preparation stage, the wind tunnel has not yet started, and the aircraft test model will not be affected by aerodynamic force.Then the wind tunnel test enters the start-up stage, and the flow field around the aircraft test model is in the establishment process.The aircraft test model will be subjected to nonstationary aerodynamic forces at this stage.Finally, in the effective experimental stage of the wind tunnel, the flow field around the aircraft test model becomes stable,and the aircraft test model will be subjected to stable aerodynamic forces at this stage.

    (2) Inertia force.In the wind tunnel start-up stage,the transient airflow is applied to the aircraft test model.The strong aerodynamic force causes the FMS to vibrate,and the effective test time of the pulse combustion wind tunnel is short,only 500 ms.The vibration of the aircraft test model is difficult to rapidly attenuate,which in turn introduces the inertial force signal into the output signals of FMS.Because of the high-amplitude periodic components of inertia force with different frequencies,the aerodynamic force signal is submerged by the inertial force signal, as shown in the right middle position of Fig.2.

    (3) Quantization noise.The strain sensor outputs analog signals.ADC is needed to convert the aerodynamic force from the hypersonic wind tunnel test into a digital signal.The quantization bits of ADC are limited, so ADC can be regarded as a nonlinear mapping from continuous amplitude input to discrete amplitude output,which causes quantization noise.Quantization noise is a random error in output signals of FMS.Usually, the frequency of the quantization noise is high, but the amplitude of the quantization noise is small, as shown in the lower right quarter of Fig.2.

    2.2.Mathematical description of FMS output signals

    Based on the above analysis,it is not difficult to find that in the effective experimental stage of the pulse combustion wind tunnel aerodynamic test, the output signals of FMS mainly include the following three kinds: (A) the stationary aerodynamic force component with almost zero frequency; (B) periodic components of inertia force with different frequencies;(C)low amplitude and high frequency quantization noise component.The time-displacement of the balance signal could be modeled as follows:

    3.Proposed method

    The TLN-AE&ST method proposed in this study mainly includes the following two parts:

    (1) To improve the identification accuracy of the aerodynamic force intelligent identification model, the residual attention block with soft threshold is introduced into the aerodynamic force intelligent identification model to filter the instrument noise signal, and a Dense block with adaptive EMD is innovatively proposed, which greatly filters the inertial force noise introduced by the largescale test model.

    (2) In view of the scarcity of wind tunnel test samples, the idea of transfer learning is proposed to achieve the proliferation of wind tunnel test samples.Firstly, the FMS dynamic calibration test bench is used to simulate the loading of aerodynamic force on FMS in wind tunnel test,and a large number of training samples for aerodynamic intelligent recognition model are obtained.Subsequently, the feature similarity between the training samples generated by the test bench and the training samples generated by the wind tunnel test is evaluated,and the training samples generated by the test bench are updated to generate more training samples similar to the wind tunnel test.The above improvements can provide more abundant training sample sets for the aerodynamic intelligent identification model and improve the identification accuracy of the aerodynamic intelligent identification model.

    3.1.Residual attention block with soft threshold

    The residual attention block with soft threshold in this paper is a variant of the residual block.The module is mainly used to filter out the instrument noise component in the FMS output signal,and avoid the effect of instrument noise on the learning and identification of inertial force characteristics by the intelligent model.

    Fig.3(a)is the residual attention block,and Fig.3(b)is the residual attention block with soft threshold.Different from the residual block, the residual attention block with soft threshold has an attention mechanism module to calculate the soft threshold for each channel.In the attention mechanism module, the absolute value of feature map is pooled globally and reduced to 1-dimensional vector.In the attention mechanism module,the absolute value of the feature map is reduced to onedimensional vector by a Global Average Pooling (GAP) operation.Subsequently,one-dimensional vector is processed into scaling parameter through two fully connected layers.The two Fully Connected (FC) layers with batch normalization and leaky ReLU activation function accelerate network training by reducing internal covariate shifts.Finally,the output of the fully connected layer is scaled to the range of (0,1) using the sigmoid activation function, and the process can be expressed as:

    Fig.3 Structure of residual attention block and residual attention block with a soft threshold.

    where zcand αcare the feature and scaling parameter of the c th channel, respectively.The threshold γ of the soft threshold function is the average value of the scaling parameter αcmultiplied by the data of the c th channel.The purpose of this operation is to determine an adaptive dynamic threshold.Thus, the threshold γ of soft threshold function can be obtained as follows:

    where i, j and c are the width index, height index and channel index of the feature graph x,respectively;γcis the threshold of the c th channel of the feature graph x.It can be seen from Eq.(3), the threshold γcmust be positive and cannot reach the maximum value of the data in the c th channel.At this time,the threshold γccan be maintained in a reasonable range to prevent the output features from being completely zero.

    Multiple residual blocks with soft threshold are stacked so that the model can learn features through nonlinear changes.In this case, soft thresholding is used as a shrinkage function to eliminate noise-related features.The thresholds are adaptively updated in the model architecture and do not need to be set manually by researchers with specialized knowledge.

    3.2.Dense block with adaptive EMD

    EMD and mutual information entropy can effectively identify the trend characteristics of aerodynamics and suppress lowfrequency noise caused by inertial force, which can effectively improve the identification accuracy of the model.The structure of dense block with adaptive EMD is shown in Fig.4.In dense block with adaptive EMD, the dense connection between layers is processed by adaptive EMD.After adaptive EMD processing, the trend feature of aerodynamic force F(t) is obvious, but the resonance of inertia force Fi(t) is suppressed.Instead, the data in the dense block is directly connected.The dense block not only introduced the resonance caused by inertial force Fi(t) into the data processed by the convolutional neural network, it will also make the model difficult to train,because the data contains more high-frequency noise.

    Fig.4 Structure of dense block with adaptive EMD.

    The input–output relationship of convolutional neural network is shown in the following equation:

    where l is the layer index of convolutional neural network;Hl(?)is a non-linear transformation of the l th layer in the network;tland tl-1are the output of the l th layer and the l-1 th layer,respectively.The dense block improves the flow of information between layers.Because, it introduces a direct connection from any layer to all subsequent layers.That is, the l th layer has the ability to integrate the feature maps of all preceding layers:

    Fig.5 Process of TLN-AE&ST.

    where [t0,t1,...,tl-1] is the integrate of the feature maps generated from the 1st layer to the l-1 th layer.The dense block with adaptive EMD optimizes the connection mode of feature maps, and the structural are shown in Fig.5.The connection method is improved from direct connection to pre-processing the feature maps from 1st layer to the l-1 th layer by adaptive EMD, as follows:

    where HEMD(?) is the adaptive EMD.It is not difficult to find that based on the above structure, the convolution layer can correct the output of adaptive EMD layer by layer.This improvement makes the model’s output close to the real aerodynamic force while retaining the rich high-frequency impact.

    3.3.Procedure of TLN-AE&ST

    Fig.5 shows the process of TLN-AE&ST.The main steps of TLN-AE&ST are as follows: high-frequency noise suppression,feature distribution alignment,aerodynamic force correction, and aerodynamic source identification.

    Step 1.Suppression of high-frequency noise

    The main purpose of Step 1 is to reduce the high-frequency noise resulting from measuring aerodynamic force.Therefore,the feature extraction part of the model is composed of stacked residual attention blocks with soft threshold.Also, this step uses examples from the source and target domains to train the model.

    Methods such as autoregressive moving average inverse model, Kalman filter, recursive least square, L2 norm have been used by researchers to suppress high-frequency noise in the signal by the aerodynamic force measurement.26It is not difficult to find that the essence of these methods is to parse the signal into real aerodynamic force through a transfer matrix H, as follows:

    However, the matrix H is usually an ill-conditioned matrix due to the increase in the dimension of solving the inverse problem, that is, the matrix H is a non-singular matrix with a large condition number.Therefore,the matrix H is extremely sensitive to the noise in the output F(t) of the FMS, and it is difficult to output real aerodynamic force.

    To solve the above problem, this paper uses a residual attention block with a soft threshold for characterization learning of the data.This block is a deep learning structure with multiple convolutional layers, which can adaptively filter instrument noise.Further, the method of stacking residual blocks can form a more abstract high-level representation of low-level features to completely describe the aerodynamic features.

    Step 2.Transfer of features

    The main purpose of Step 2 is to make network learning source domain and target domain signal common features.The feature similarity between the training samples generated by the test bench and the training samples generated by the wind tunnel test is evaluated, and the training samples generated by the test bench are updated to generate more training samples similar to the wind tunnel test.The output signal of hypersonic FMS is called target domain signal.The signal of FMS test bench is called source domain signal,which also conforms to the above feature.The similarity of trends between the two signals makes the introduction of transfer learning possible.Using these signals for training can make the network quickly understand the trend features of aerodynamic force.

    The output of an input signal in the neural network should be consistent with its target value,and the consistency evaluation function is called the objective function.In this paper,the cross-entropy loss function is used to evaluate the consistency between the source domain signal and the target domain signal.To make the cross-entropy loss function focus on the difference of signal trend,rather than the difference of amplitude,the softmax activation function is used in the output of Step 1 to convert the aerodynamic trend into a probability distribution with sum of 1.The process is as follows:

    where m is the size of the mini-batch.

    Step 3.Correction of aerodynamic force

    The main purpose of Step 3 is to suppress the inertial force Fi(t)in the signal.After processing Step 1 and Step 2,the highfrequency noise caused by sensor noise of the aerodynamic force measurement have been suppressed,but the inertial force Fi(t) still causes interference to the real aerodynamic force.Therefore, Dense block with adaptive EMD is adopted to reduce the influence of inertial force Fi(t) on the identification accuracy of aerodynamic loads.

    Through EMD reconstruction, the network can quickly identify the trend of aerodynamic components in the signal.However, EMD has no training parameters, for a fixed input,the reconstructed signal is the same.Therefore, the impact component of the aerodynamic force is suppressed as highfrequency noise by the EMD.To address the above situation,this paper uses convolutional neural networks to dynamically compensate the immersion impact component.The Step 3 can improve the accuracy of pneumatic identification, and the process is as follows:

    Step 4.Learning of source domain features

    The main purpose of Step 4 is to make the network learn the trend features of aerodynamic force.MMD, a nonparametric method, is adopted in this paper to measure the difference between the network output signals and the real aerodynamic signals in the source domain.The principle of MMD: (A) the aerodynamic trend features in the source and target domains are mapped to the Regenerated Kernel Hilbert Space (RKHS) through a high-dimensional mapping function f:X →S.(B) MMD is the supremum of the expected difference between two domain features in RHKS, as follows:

    Step 5.Optimization object

    In the above part,the function of each module of the TLNAE&ST and its corresponding loss function are introduced respectively.The optimization objective of the model is established by comprehensively considering each loss function, and the TLN-AE&ST is trained by random gradient descent method.Lcrsand LMMDare the difference of domain feature and the difference of source domain trend feature extraction,respectively.The optimization objective of TLN-AE&ST can be expressed as:

    4.Experiment

    4.1.Experimental setup

    Fig.6 Three-dimensional model of suspension force measuring system.

    In the field of hypersonic wind tunnel aerodynamic test, the aircraft test model is developing towards large-scale and heavy-duty (the large-size aircraft test model is 7 m-8 m in length and several tons in weight, compared with the scaledsize aircraft test model less than 2 min length and 1 ton in weight).When the aerodynamic test of hypersonic wind tunnel is carried out on the large-scale and heavy-duty aircraft test model,it is required that the FMS has high stiffness,small disturbance to the flow field, and can carry out the airframe/propulsion integration test at the same time.The traditional force measurement technology is difficult to meet the above requirements at the same time.In view of the above problems,our research group proposed a new FMS, namely the suspension FMS,as shown in Fig.6.The device has been successfully applied in φ = 2.4 m pulse combustion wind tunnel of China Aerodynamics Research and Development Center(CARDC).27Therefore,the suspension FMS is selected to verify the effectiveness of TLN-AE&ST.The working principle of the suspension force measuring system is shown in Fig.7.During the hypersonic wind tunnel test,the aerodynamic force acting on the aircraft model is transmitted to the force measurement sensors through the force measurement pull rods, and the voltage signal is output to realize the aerodynamic force test.

    The specific process of verifying the effectiveness of the aerodynamic intelligent identification model (TLN-AE&ST)is as follows:

    (1) The sample data used for TLN-AE&ST model training is very important, which will directly affect the training process and results.To obtain high-quality training sample data, this study designed a force loading device that can apply a step force to an aircraft test model to simulate the aerodynamic forces experienced by the aircraft test model in hypersonic wind tunnel tests, as shown in the left upper part of Fig.8.The step load is achieved by suddenly releasing weights of known mass.The loading point is fixed on the surface of the aircraft test model.The step force transfer path from the loading point to the weight is shown in the middle left part of Fig.8.The principle of generating step force is as follows: first place a known weight on the weight plate,then send a control command to make the electromagnetic sucker power off and lose suction, thus generating step force.In the test,the force sensor measures the step force signal, and the suspension FMS outputs the response signal.The multi-component dynamic response signals were collected by LMS SCADAS Mobile SCM 202 acquisition instrument.In this study, 126 sets of training samples were generated based on the dynamic calibration test bench, i.e.source domain samples.

    (2) The suspension FMS was tested in the φ= 2.4 m pulse combustion wind tunnel of CARDC,and its output signals were used as the target domain sample,as shown in Fig.9.This study obtained a total of 12 sets of training sample sets(real wind tunnel test signals),namely target domain samples.The φ=2.4 m pulse combustion wind tunnel is mainly used to carry out the performance test research of large-scale scramjet and the performance test research of integrated hypersonic vehicle with power.The test contents and parameters of wind tunnel test are shown in Table 1.The test conditions are Ma=5.5,6.5,7.0.The aircraft test model angle of attack and sideslip are both 0°.The suspension FMS used in the wind tunnel test in this paper is shown in the right middle part in Fig.8.In the figure,Fx,Fy,Fz,Mx,Myand Mzrepresent the output signals of six degrees of freedom (normal,axial, yaw, pitch, roll and nod) of the force measuring system.

    Fig.7 Principle diagram of suspension force measuring system.

    (3) 126 source domain samples and 12 target domain samples were used to train and verify the accuracy and effectiveness of the TLN-AE&ST model.

    In this section, various traditional aerodynamic identification methods are compared with TLN-AE&ST.The methods involved in this section are shown in Table 219,22,28–30.

    Frequency domain model: Considering the complexity of wind tunnel FMS signal, the time–frequency aerodynamic force identification model is used to reduce noise signal and process the wind tunnel FMS signal of the pulse wave of the aircraft test model, and output reliable aerodynamic force results.In this paper, three mainstream aerodynamic force identification models are used as comparison,which are Wavelet Transform (WT),19EMD28and Extremum field Mean Mode Decomposition (EMMD).29

    Intelligence model:Since 2017,the dynamic self-calibration model based on deep learning technology has been applied to the FMS of wind tunnel (especially hypersonic high-enthalpy wind tunnel).The intelligent model obtains the ideal aerodynamic force by dynamically correcting the force measurement results of inertial vibration interference.In this paper, CNN and RNN are selected as the comparison models.22,30

    4.2.Evaluation indicators and results

    The model is evaluated by Mean Absolute Error (MAE),Mean Square Error (MSE) and Root Mean Square Error(RMSE).MAE reflects the average absolute error between real and predicted values.MSE reflects the error between real and predicted values expectation,that is,the stability of the model.RMSE reflects the square root of the mean of the squared differences between true and predicted values, and is more sensitive to outliers than MSE.The formulas for the different evaluation indicators are as follows:

    where y and ^y are the true and predicted value, respectively.n is the number of data points.t is the number of channels output by the FMS.The real aerodynamic force in the target domain is obtained by filtering mean method.The TLNAE&ST is compared with models in Table 2 by the above evaluation indexes.The results are shown in Fig.10, and the detailed data are shown in Table 3.The following conclusions can be drawn from the analysis:

    Fig.8 Structure and output signals of rod-type suspension force measuring device and its scaled-down force measurement device.

    (1) Compared with the mainstream frequency domain models and intelligent models in the field of FMS, TLNAE&ST achieves the best results.

    (2) Among the frequency domain models, the MEA results of the three models are close, but the MSE results of EMD and EEMD are significantly smaller than WT.This is because the WT can not effectively suppress the pulse signal caused by inertial force in large scale balance.On the contrary, EMD and EEMD models have stronger anti-interference ability.

    (3) Interestingly, in the scenario of aerodynamic feature transfer, the identification accuracy of the intelligent models is lower than that of the traditional frequency domain models.The parameters of CNN and RNN are fixed after training.Therefore, they only have good identification effect for a specific FMS (that is, both the training and the test set come from the source domain), and lack sufficient robustness.TLN-AE&ST can learn the common features of source and target domain based on Eq.(14), and still maintain good aerodynamic identification effect in the target area.

    Table 1 Test contents and parameters of φ = 2.4 m pulse combustion wind tunnel.

    Fig.11 Performance of each model under different sample sizes.

    Table 2 Task settings for each model.

    Fig.10 Performance of each model.

    Table 3 Detailed performance of each model.

    4.3.Model performance analysis

    4.3.1.Sample size analysis

    The results of intelligent models are dependent on the number of samples due to their excessive trainable parameters.In fields such as face recognition, natural language processing, and fault diagnosis,intelligent models are established based on sufficient data.However, in the FMS of hypersonic wind tunnel,due to the expensive experimental cost, there are not enough samples to train the model.Therefore, it is necessary to studythe performance of intelligent models in the scarcity of samples.In this section, four datasets with different specifications are selected as the source domain data, namely 10,25,50 and 100.The MSE was selected to evaluate the performance of models.and the results are shown in Fig.11, and the detailed data are shown in Table 4.

    Table 4 Detailed performance of each model under different sample sizes.

    The frequency domain methods are not affected by the number of samples,all of which can maintain good discrimination accuracy.When the number of samples in the source domain is less than 25, the TLN-AE&ST is underfitted and its effect is slightly worse than EEMD.When the source domain samples are more than 50, the TLN-AE&ST is fully trained,and its identification effect is the best among all models.On the contrary, the identification accuracy of intelligent models such as CNN and RNN decreases significantly when the source domain samples are scarce (less than 50).When the source domain samples are 100, these models can only maintain the similar identification accuracy of WT.It can be seen that such intelligent models have a certain distance in the large-scale application of FMS.

    Fig.12 Performance of TLN-AE&ST under different sample sizes.

    For each dataset, the circle on the left is the MSE for each experiment, and the violin plot on the right is the logarithmic normal fitting to the data from 50 replicate experiments.

    To further analyze the stability of TLN-AE&ST, 50 replicate experiments were performed with different sample sizes,and the results are shown in Fig.12.It is not difficult to find that when the sample size is 10, the stability of the model is poor, and the average and the standard deviation of the 50 experiments are 0.108 and 0.017, respectively.When the number of samples is 50 and 100,the identification accuracy of the model is slightly improved, but its stability is significantly improved, and its average is increased from 0.085 to 0.082.

    4.3.2.Model complexity analysis

    Real-time aerodynamic force identification is a challenge for most intelligent models at present.This is because, the FMS is a typical industrial environment with sensitive and private output signals that cannot be processed using online highperformance servers.Real-time identification of aerodynamic force on the limited computing resources of mobile devices is also one of the problems solved in this study.In order to evaluate the usability of intelligent models in FMS, it is necessary to validate the model complexity in addition to the aerodynamic identification accuracy.The evaluation indexes of model complexity include model parameters and the multiplyaccumulate operations(MADDs)time.Fig.13 shows the comparison between TLN-AE&ST and the current mainstream intelligent models.The radius of the circle represents the number of model parameters.

    The experimental platform is Jetson TX2.To prevent the interference of other tasks on the platform,20 repeated experiments were carried out and the minimum value was taken.As shown in the figure, the RNN model has achieved the fastest operation time of 265 ms, which fully meets the real-time requirements of aerodynamic identification on the mobile platform.The parameter number of TLN-AE&ST is 1.3 M,which has a significant reduction compared with CNN and RNN.The number of model parameters also explains that TLNAE&ST can be fully trained when there are only 50 samples in the source domain, while CNN and RNN are underfitting.

    4.3.3.Ablation experiment

    (1) TLN-AE&ST w/ or w/o cross entropy function

    Fig.13 Complexity of each model.

    Fig.14 Effect comparison of models in source domain and target domain.

    To clearly determine the role of the cross-entropy function in STEP-2, TLN-AE&ST and TLN-AE&ST without crossentropy function are plotted in Fig.14 for both source and target domain effects.In the source domain, both methods achieve ideal results with MES of 0.06.However,in the target domain, TLN-AE&ST without cross-entropy function cannot effectively identify aerodynamic signal, and its performance declines sharply.Because TLN-AE&ST without crossentropy function is unable to learn the aerodynamic signal features during the training process.TLN-AE&ST, due to its cross-entropy function, can learn the common features of the aerodynamic force in the source and the target domain, so it maintains a good identification effect.

    (2) TLN-AE&ST w/ or w/o Residual attention block

    To verify the importance of the attention mechanism module on the accuracy of aerodynamic force identification, this study uses the experimental data of the φ=2.4 m pulse combustion wind tunnel to analyze the influence of TLN-AE&ST and TLN-AE&ST without residual attention block on the accuracy of aerodynamic identification.The test conditions are Ma=6.5.The aircraft test model angle of attack and sideslip are both 0°.The identification results are shown in Fig.15.The blue line in Fig.15 is the output signals of TLN-AE&ST without Residual attention block.It is not difficult to find that the output signals contain a large amount of burr noise,which shows that the network cannot deal with medium–high frequency noise well.The output signals of TLN-AE&ST are smoother,and the MAE and MSE are much higher than those of TLN-AE&ST without residual attention block.In conclusion,residual attention block with soft threshold is an effective method to suppress medium–high frequency noise.And TLNAE&ST can combine the advantages of convolutional neural network and soft threshold to adaptively reduce the medium–high frequency noise component in the signal.

    Fig.15 Output signals of TLN-AE&ST w/ or w/o residual attention block.

    Fig.16 Output signals of TLN-AE&ST w/ or w/o dense block.

    (3) TLN-AE&ST w/ or w/o without dense block

    From a certain point of view, the implementation process of dense block can be regarded as the intelligent filtering process of inertial force, which can be more clearly understood from the following spectrum comparison.From Fig.16,FFT is performed on the output signal of FMS in φ=2.4 m pulse combustion wind tunnel,the signal processed by the TLN-AE&ST model and TLN-AE&ST without dense block.It can be clearly seen that the main frequencies of the inertial force in the wind tunnel test are completely filtered out.The upper right corner of Fig.16 is the spectrogram of the output signals of TLN-AE&ST without dense block.As shown by the yellow line, the spectrogram has the selfoscillation frequency near 100 Hz, which can be inferred that the signals are disturbed by the inertial force caused by the intrinsic frequency of the FMS.The lower right corner of Fig.16 is the spectrogram of the output signals of TLNAE&ST.As shown by the red dotted line, the spectrogram has the self-oscillation frequency near 20 Hz.It can be seen that TLN-AE&ST can effectively eliminate the inertial force interference caused by the intrinsic frequency.

    5.Conclusions and future work

    At present, the research on the intelligence of FMS based on deep learning is relatively new, and the related methods and technologies are still in the development stage.In this paper,the transfer learning is introduced into the hypersonic FMS for the first time, which greatly reduces the requirements of sample quantity and quality for intelligent model.Moreover,the residual attention block with soft threshold and dense block with adaptive EMD are proposed, which can not only accurately identify the aerodynamic signal, but also effectively suppress other interference signals.TLN-AE&ST has application prospect,which provides key technology and data support for accurate evaluation of hypersonic vehicle aerodynamic characteristics with high enthalpy.

    Although aerodynamic identification results of cross-model size are promising, some limitations of this work should be pointed out.First, in this study, the FMS involved in the test is single.The aerodynamic identification network with good robustness is not considered (by using the dataset of the existing full-scale FMS)and used to guide the design of new FMS.Second,time–frequency transform method is widely studied in the field of FMS.By incorporating the latest research results,enhancing the identification capabilities of TLN-AE&ST is a task that should be pursued in the future.

    Declaration of Competing Interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgements

    This study was co-supported by the National Natural Science Foundation of China (52105562), the Fundamental Research Funds for the Central Universities, China (XJ2021KJZK037),the Fundamental Research Funds for the Central Universities,China (2682022CX058).

    国产日韩欧美视频二区| 黄色怎么调成土黄色| 一区二区av电影网| 高清黄色对白视频在线免费看| 久久av网站| videos熟女内射| 黄色视频在线播放观看不卡| 一二三四中文在线观看免费高清| 18在线观看网站| 成人午夜精彩视频在线观看| 免费黄色在线免费观看| 国产99久久九九免费精品| 一本色道久久久久久精品综合| 日韩制服骚丝袜av| 一级片免费观看大全| 80岁老熟妇乱子伦牲交| 无遮挡黄片免费观看| 老司机在亚洲福利影院| 国产精品欧美亚洲77777| 国产精品一国产av| 欧美日韩福利视频一区二区| 18在线观看网站| 午夜福利网站1000一区二区三区| 黄片小视频在线播放| 国产精品香港三级国产av潘金莲 | 午夜免费男女啪啪视频观看| 久久久久网色| 色播在线永久视频| 午夜福利视频精品| 日日爽夜夜爽网站| 欧美中文综合在线视频| av在线观看视频网站免费| 99国产综合亚洲精品| 深夜精品福利| 男女无遮挡免费网站观看| 国产乱人偷精品视频| 中文字幕精品免费在线观看视频| 黄色怎么调成土黄色| 操美女的视频在线观看| 欧美成人午夜精品| 亚洲色图综合在线观看| 久久免费观看电影| 狠狠婷婷综合久久久久久88av| 美女脱内裤让男人舔精品视频| 黄色一级大片看看| 蜜桃国产av成人99| 国产免费视频播放在线视频| 一本一本久久a久久精品综合妖精| 国产精品久久久久成人av| 国产人伦9x9x在线观看| 大片电影免费在线观看免费| 色综合欧美亚洲国产小说| 欧美人与性动交α欧美精品济南到| 超色免费av| 高清黄色对白视频在线免费看| videos熟女内射| 欧美日韩av久久| 亚洲精品国产av成人精品| 中文乱码字字幕精品一区二区三区| 美女国产高潮福利片在线看| 亚洲精品第二区| 亚洲欧美精品综合一区二区三区| 99久国产av精品国产电影| 成年人午夜在线观看视频| 亚洲欧洲日产国产| 一级毛片电影观看| 男女国产视频网站| 男女国产视频网站| 9色porny在线观看| 中国国产av一级| 欧美精品av麻豆av| 日本av手机在线免费观看| 亚洲自偷自拍图片 自拍| 又黄又粗又硬又大视频| 99久国产av精品国产电影| 亚洲 欧美一区二区三区| 精品卡一卡二卡四卡免费| 一本大道久久a久久精品| 下体分泌物呈黄色| 妹子高潮喷水视频| 久久久久久久大尺度免费视频| 久久精品国产a三级三级三级| 乱人伦中国视频| 成年人午夜在线观看视频| 国产精品嫩草影院av在线观看| 色婷婷久久久亚洲欧美| 老司机影院成人| 丝袜脚勾引网站| 国产亚洲av高清不卡| a级毛片黄视频| 国产日韩欧美视频二区| 晚上一个人看的免费电影| 亚洲精品国产av蜜桃| av线在线观看网站| 天天躁夜夜躁狠狠躁躁| 日日啪夜夜爽| 国产亚洲精品第一综合不卡| 国产精品香港三级国产av潘金莲 | 国产国语露脸激情在线看| 国产97色在线日韩免费| 欧美少妇被猛烈插入视频| 国产高清国产精品国产三级| 久久久久精品国产欧美久久久 | 欧美国产精品一级二级三级| 中文字幕最新亚洲高清| √禁漫天堂资源中文www| 又粗又硬又长又爽又黄的视频| 自拍欧美九色日韩亚洲蝌蚪91| 操出白浆在线播放| 亚洲av中文av极速乱| 19禁男女啪啪无遮挡网站| 成人免费观看视频高清| 91成人精品电影| 老司机在亚洲福利影院| 久久久精品国产亚洲av高清涩受| 香蕉国产在线看| 无限看片的www在线观看| 国产激情久久老熟女| 日韩大码丰满熟妇| 国产熟女午夜一区二区三区| 99热网站在线观看| av免费观看日本| 高清黄色对白视频在线免费看| 岛国毛片在线播放| 不卡av一区二区三区| 高清不卡的av网站| 国产激情久久老熟女| 亚洲美女黄色视频免费看| 成人毛片60女人毛片免费| 精品少妇黑人巨大在线播放| 热99久久久久精品小说推荐| 欧美变态另类bdsm刘玥| 中文乱码字字幕精品一区二区三区| 2021少妇久久久久久久久久久| 在线观看国产h片| 国产精品人妻久久久影院| 成年女人毛片免费观看观看9 | 国产精品麻豆人妻色哟哟久久| 国产精品国产三级国产专区5o| www日本在线高清视频| 制服人妻中文乱码| 天天躁夜夜躁狠狠久久av| 久久女婷五月综合色啪小说| 欧美激情高清一区二区三区 | 久久久精品免费免费高清| 韩国精品一区二区三区| 熟妇人妻不卡中文字幕| 中文字幕色久视频| 国产福利在线免费观看视频| 亚洲成色77777| 国产淫语在线视频| 老司机在亚洲福利影院| 天美传媒精品一区二区| 日韩av在线免费看完整版不卡| 国产又色又爽无遮挡免| 无限看片的www在线观看| 99热国产这里只有精品6| 国产av一区二区精品久久| 国产乱人偷精品视频| 精品久久久精品久久久| 日韩制服骚丝袜av| 日本黄色日本黄色录像| 精品人妻一区二区三区麻豆| 99久久99久久久精品蜜桃| 久久精品亚洲av国产电影网| 久热这里只有精品99| 一区二区三区乱码不卡18| 一边亲一边摸免费视频| 最近最新中文字幕大全免费视频 | 精品一区二区三区四区五区乱码 | 丝袜美腿诱惑在线| 青春草视频在线免费观看| 男女午夜视频在线观看| 黄网站色视频无遮挡免费观看| 麻豆精品久久久久久蜜桃| 国产精品一区二区精品视频观看| 精品久久久精品久久久| www.熟女人妻精品国产| 日韩不卡一区二区三区视频在线| 十八禁人妻一区二区| 成年美女黄网站色视频大全免费| 国产高清不卡午夜福利| 久久久久精品人妻al黑| 少妇精品久久久久久久| 亚洲精品国产av成人精品| 爱豆传媒免费全集在线观看| 亚洲国产精品国产精品| 韩国高清视频一区二区三区| 一区二区三区四区激情视频| 国产精品免费大片| 日韩,欧美,国产一区二区三区| 高清欧美精品videossex| 免费在线观看黄色视频的| 精品久久蜜臀av无| 黑丝袜美女国产一区| 国产爽快片一区二区三区| 好男人视频免费观看在线| 波多野结衣一区麻豆| 久久人妻熟女aⅴ| 久久狼人影院| 丝瓜视频免费看黄片| 最近2019中文字幕mv第一页| av网站在线播放免费| 欧美亚洲 丝袜 人妻 在线| 新久久久久国产一级毛片| netflix在线观看网站| 亚洲成国产人片在线观看| 日韩一本色道免费dvd| 国产男女超爽视频在线观看| 免费黄频网站在线观看国产| 综合色丁香网| 一级片'在线观看视频| 免费观看a级毛片全部| 欧美黑人欧美精品刺激| 天天躁夜夜躁狠狠躁躁| 国产免费一区二区三区四区乱码| 精品一区在线观看国产| 爱豆传媒免费全集在线观看| 国产男女超爽视频在线观看| 岛国毛片在线播放| 热99久久久久精品小说推荐| 亚洲少妇的诱惑av| 亚洲欧美日韩另类电影网站| 日日摸夜夜添夜夜爱| 欧美精品av麻豆av| 最新在线观看一区二区三区 | 日韩电影二区| 国产日韩欧美亚洲二区| 亚洲专区中文字幕在线 | 欧美黑人欧美精品刺激| 精品国产一区二区三区四区第35| www.自偷自拍.com| 亚洲国产成人一精品久久久| 亚洲熟女毛片儿| 亚洲国产欧美日韩在线播放| 亚洲精品日本国产第一区| 亚洲精品国产区一区二| 国产精品蜜桃在线观看| 女人被躁到高潮嗷嗷叫费观| 韩国精品一区二区三区| 97在线人人人人妻| 国产色婷婷99| 另类亚洲欧美激情| 欧美激情 高清一区二区三区| 色网站视频免费| 香蕉丝袜av| 亚洲精品自拍成人| 超色免费av| 亚洲精品美女久久久久99蜜臀 | 国产免费福利视频在线观看| 亚洲中文av在线| 国产欧美亚洲国产| 男女之事视频高清在线观看 | 下体分泌物呈黄色| 久久精品亚洲熟妇少妇任你| 国产精品 欧美亚洲| 亚洲综合色网址| 色综合欧美亚洲国产小说| 国产成人欧美在线观看 | 别揉我奶头~嗯~啊~动态视频 | 国产精品麻豆人妻色哟哟久久| 亚洲少妇的诱惑av| 免费黄频网站在线观看国产| 九色亚洲精品在线播放| 热99国产精品久久久久久7| 亚洲自偷自拍图片 自拍| xxx大片免费视频| 久久午夜综合久久蜜桃| 热re99久久国产66热| 美国免费a级毛片| 久久久久精品性色| 国产精品国产三级国产专区5o| 国产成人一区二区在线| 亚洲av在线观看美女高潮| 免费少妇av软件| 99re6热这里在线精品视频| 女性生殖器流出的白浆| 深夜精品福利| 黄频高清免费视频| 成年人免费黄色播放视频| 日韩视频在线欧美| 国产一区有黄有色的免费视频| 久热这里只有精品99| 亚洲第一av免费看| 久久久久久人妻| 老汉色∧v一级毛片| 狠狠精品人妻久久久久久综合| 亚洲专区中文字幕在线 | 丝袜在线中文字幕| 美女中出高潮动态图| 亚洲国产精品999| 亚洲人成77777在线视频| 香蕉丝袜av| 捣出白浆h1v1| 大陆偷拍与自拍| 亚洲第一青青草原| 免费高清在线观看日韩| 制服丝袜香蕉在线| 女人久久www免费人成看片| 国产精品二区激情视频| 欧美av亚洲av综合av国产av | 曰老女人黄片| 成年av动漫网址| 考比视频在线观看| 99热全是精品| 91成人精品电影| 成人亚洲欧美一区二区av| 国产精品偷伦视频观看了| 亚洲成国产人片在线观看| 嫩草影院入口| 亚洲精品国产av成人精品| 久久久久久人人人人人| 日韩 欧美 亚洲 中文字幕| 亚洲一码二码三码区别大吗| 国产成人免费观看mmmm| e午夜精品久久久久久久| 女人精品久久久久毛片| a级毛片在线看网站| 中文字幕另类日韩欧美亚洲嫩草| 精品国产乱码久久久久久男人| 97在线人人人人妻| 在线观看人妻少妇| 亚洲精品一区蜜桃| 久久国产精品男人的天堂亚洲| 涩涩av久久男人的天堂| 波多野结衣av一区二区av| 欧美av亚洲av综合av国产av | 免费观看a级毛片全部| 一区二区av电影网| 国产精品免费大片| 国产免费视频播放在线视频| 国产黄色视频一区二区在线观看| 999精品在线视频| 巨乳人妻的诱惑在线观看| 久久精品aⅴ一区二区三区四区| 欧美激情 高清一区二区三区| 亚洲在久久综合| 精品免费久久久久久久清纯 | 一级片'在线观看视频| 成人漫画全彩无遮挡| 又黄又粗又硬又大视频| 在线亚洲精品国产二区图片欧美| 久久精品久久久久久久性| 精品亚洲乱码少妇综合久久| 国产不卡av网站在线观看| 国产成人精品在线电影| 18在线观看网站| 亚洲av男天堂| 久久精品久久精品一区二区三区| 国产精品久久久久久精品电影小说| 久久久久久久大尺度免费视频| 在线观看免费视频网站a站| 国产在线视频一区二区| 亚洲欧美一区二区三区黑人| 欧美精品一区二区免费开放| 亚洲美女黄色视频免费看| av又黄又爽大尺度在线免费看| 亚洲色图综合在线观看| 国产精品久久久久久精品电影小说| 亚洲图色成人| 国产精品国产av在线观看| 精品国产一区二区三区久久久樱花| 丝袜人妻中文字幕| 免费在线观看完整版高清| 一区二区三区精品91| 校园人妻丝袜中文字幕| 久久久久久久国产电影| 一级,二级,三级黄色视频| 一区二区三区四区激情视频| 97在线人人人人妻| 精品福利永久在线观看| 少妇人妻 视频| 人人澡人人妻人| 永久免费av网站大全| 午夜91福利影院| 99热全是精品| 精品少妇一区二区三区视频日本电影 | 在线观看www视频免费| 国产爽快片一区二区三区| 亚洲第一青青草原| 亚洲国产精品成人久久小说| 欧美日韩av久久| 国产日韩欧美亚洲二区| 亚洲伊人色综图| 国产精品久久久久成人av| 男女边摸边吃奶| 在线观看免费日韩欧美大片| 狂野欧美激情性xxxx| 丰满迷人的少妇在线观看| 国产乱来视频区| 男人舔女人的私密视频| 亚洲第一青青草原| 国产有黄有色有爽视频| 久久性视频一级片| 一级毛片 在线播放| 久久久久精品国产欧美久久久 | 人人妻人人澡人人看| 91精品伊人久久大香线蕉| 国产精品蜜桃在线观看| 精品少妇黑人巨大在线播放| 中文字幕人妻熟女乱码| 亚洲欧洲精品一区二区精品久久久 | 又粗又硬又长又爽又黄的视频| 热re99久久国产66热| 国产免费一区二区三区四区乱码| 国产亚洲av片在线观看秒播厂| 午夜福利乱码中文字幕| 欧美xxⅹ黑人| 免费在线观看黄色视频的| 天堂8中文在线网| av女优亚洲男人天堂| av在线播放精品| 精品亚洲乱码少妇综合久久| 99国产综合亚洲精品| 亚洲国产av新网站| www.精华液| 亚洲五月色婷婷综合| 久久亚洲国产成人精品v| 久久天躁狠狠躁夜夜2o2o | 成年美女黄网站色视频大全免费| 一区二区av电影网| 精品第一国产精品| 日韩一区二区三区影片| 一区福利在线观看| 高清黄色对白视频在线免费看| 男女午夜视频在线观看| 久久天堂一区二区三区四区| 丝袜美腿诱惑在线| 国产精品亚洲av一区麻豆 | 国产色婷婷99| 中文欧美无线码| 国产免费视频播放在线视频| 成年动漫av网址| av有码第一页| 国产男女内射视频| 国产精品一二三区在线看| 涩涩av久久男人的天堂| 精品国产一区二区三区久久久樱花| 亚洲色图 男人天堂 中文字幕| 亚洲av欧美aⅴ国产| 十八禁人妻一区二区| 高清不卡的av网站| 久热这里只有精品99| 国产 一区精品| av国产久精品久网站免费入址| 精品一区二区三卡| 亚洲人成电影观看| 五月天丁香电影| 国产在线免费精品| av又黄又爽大尺度在线免费看| 亚洲精品在线美女| 国产精品av久久久久免费| 久久久久精品国产欧美久久久 | 欧美激情 高清一区二区三区| 一区二区av电影网| 亚洲精品久久午夜乱码| 一边摸一边抽搐一进一出视频| 青春草视频在线免费观看| 日本色播在线视频| 午夜91福利影院| 欧美变态另类bdsm刘玥| 国产97色在线日韩免费| 侵犯人妻中文字幕一二三四区| 你懂的网址亚洲精品在线观看| 久久久久久久久免费视频了| av线在线观看网站| 叶爱在线成人免费视频播放| 欧美久久黑人一区二区| 老熟女久久久| 国产成人精品在线电影| 爱豆传媒免费全集在线观看| 麻豆乱淫一区二区| 亚洲欧美一区二区三区黑人| av又黄又爽大尺度在线免费看| 日韩制服丝袜自拍偷拍| 欧美日韩视频精品一区| 丰满饥渴人妻一区二区三| 大陆偷拍与自拍| 在线观看三级黄色| 国产精品一二三区在线看| 日韩熟女老妇一区二区性免费视频| 国产精品欧美亚洲77777| 伊人久久国产一区二区| 日韩大码丰满熟妇| 欧美乱码精品一区二区三区| 亚洲国产精品999| 国产精品亚洲av一区麻豆 | 国产在线视频一区二区| 视频区图区小说| 777久久人妻少妇嫩草av网站| 日韩,欧美,国产一区二区三区| 国产亚洲欧美精品永久| 欧美亚洲日本最大视频资源| 狠狠婷婷综合久久久久久88av| 国产片内射在线| 国产日韩欧美亚洲二区| 大香蕉久久成人网| 亚洲国产最新在线播放| 99re6热这里在线精品视频| 中文字幕色久视频| 久久久精品区二区三区| 男女边摸边吃奶| 又黄又粗又硬又大视频| 成人免费观看视频高清| 国产精品一区二区在线不卡| 涩涩av久久男人的天堂| 欧美日韩视频精品一区| 国产乱来视频区| 亚洲国产精品999| 在线观看三级黄色| 久久久精品区二区三区| 久久女婷五月综合色啪小说| 9色porny在线观看| 久久久国产欧美日韩av| 亚洲第一青青草原| 亚洲熟女精品中文字幕| netflix在线观看网站| 高清在线视频一区二区三区| 纵有疾风起免费观看全集完整版| 免费观看性生交大片5| 日韩欧美一区视频在线观看| 日日撸夜夜添| 精品国产超薄肉色丝袜足j| 久久天躁狠狠躁夜夜2o2o | 国产成人免费观看mmmm| av女优亚洲男人天堂| 巨乳人妻的诱惑在线观看| 老鸭窝网址在线观看| 一二三四中文在线观看免费高清| 国产成人系列免费观看| av在线老鸭窝| 免费高清在线观看日韩| 制服人妻中文乱码| 国产xxxxx性猛交| 亚洲色图 男人天堂 中文字幕| av天堂久久9| 18禁观看日本| 青草久久国产| 久久久久网色| 免费观看av网站的网址| 黄色怎么调成土黄色| 亚洲精品一二三| 国产精品三级大全| 亚洲国产毛片av蜜桃av| 99久久精品国产亚洲精品| 精品国产超薄肉色丝袜足j| a级毛片在线看网站| 亚洲欧美清纯卡通| 大香蕉久久网| 成人国产麻豆网| 久久久久精品性色| 国产成人精品久久二区二区91 | 国产免费一区二区三区四区乱码| 在线观看一区二区三区激情| 无限看片的www在线观看| 晚上一个人看的免费电影| 亚洲精品视频女| 欧美国产精品va在线观看不卡| 天天添夜夜摸| 卡戴珊不雅视频在线播放| 精品酒店卫生间| 欧美国产精品一级二级三级| 嫩草影视91久久| 18在线观看网站| 日日摸夜夜添夜夜爱| 亚洲精华国产精华液的使用体验| 一本—道久久a久久精品蜜桃钙片| 美女视频免费永久观看网站| 18在线观看网站| 97精品久久久久久久久久精品| 欧美精品人与动牲交sv欧美| 超碰成人久久| 欧美xxⅹ黑人| 热99国产精品久久久久久7| 久久久久久久久久久免费av| av在线播放精品| 十八禁网站网址无遮挡| 午夜免费鲁丝| 宅男免费午夜| 欧美精品av麻豆av| 午夜福利一区二区在线看| 免费久久久久久久精品成人欧美视频| 亚洲第一av免费看| 亚洲成国产人片在线观看| 女人高潮潮喷娇喘18禁视频| 日韩精品有码人妻一区| 国产成人欧美| 国产熟女午夜一区二区三区| 晚上一个人看的免费电影| 中文字幕人妻熟女乱码| 老熟女久久久| 极品人妻少妇av视频| 91成人精品电影| 国产激情久久老熟女| 免费少妇av软件| 99久久精品国产亚洲精品| 国产亚洲精品第一综合不卡| 亚洲一级一片aⅴ在线观看| 免费黄频网站在线观看国产| av网站在线播放免费| 久久毛片免费看一区二区三区| 欧美久久黑人一区二区| 欧美最新免费一区二区三区| 咕卡用的链子| 亚洲av日韩在线播放| 秋霞在线观看毛片| 久久99精品国语久久久| 国产成人欧美| 国产欧美日韩综合在线一区二区| 黄色 视频免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av综合色区一区| 免费观看a级毛片全部| 晚上一个人看的免费电影| 亚洲欧洲国产日韩|