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

    Research on Optimal Preload Method of Controllable Rolling Bearing Based on Multisensor Fusion

    2024-03-23 08:17:44KuoshengJiangChengruiHanandYashengChang

    Kuosheng Jiang ,Chengrui Han and Yasheng Chang

    1School of Mechanical and Electrical Engineering,Anhui University of Science and Technology,Huainan,232001,China

    2School of Optical and Electronic Information,Suzhou City University&Suzhou Key Laboratory of Biophotonics,Suzhou,215104,China

    3State Key Laboratory for Manufacturing Systems Engineering,Xi’an Jiaotong University,Xi’an,710054,China

    ABSTRACT Angular contact ball bearings have been widely used in machine tool spindles,and the bearing preload plays an important role in the performance of the spindle.In order to solve the problems of the traditional optimal preload prediction method limited by actual conditions and uncertainties,a roller bearing preload test method based on the improved D-S evidence theory multi-sensor fusion method was proposed.First,a novel controllable preload system is proposed and evaluated.Subsequently,multiple sensors are employed to collect data on the bearing parameters during preload application.Finally,a multisensor fusion algorithm is used to make predictions,and a neural network is used to optimize the fitting of the preload data.The limitations of conventional preload testing methods are identified,and the integration of complementary information from multiple sensors is used to achieve accurate predictions,offering valuable insights into the optimal preload force.Experimental results demonstrate that the multi-sensor fusion approach outperforms traditional methods in accurately measuring the optimal preload for rolling bearings.

    KEYWORDS Multi-sensor;information fusion;neural network;preload force

    Nomenclature

    dmPitch circle diameter

    D Ball or roller diameter

    nmBall or roller track speed

    nrThe speed at which a ball or roller rotates around its own axis

    MgGyroscopic moment

    PdRadial clearance

    q Roller/raceway unit length load

    Q Ball or roller normal load

    QaAxial load on a ball or roller

    QrRadial load on a ball or roller

    QiaBall thrust load

    αContact angle

    βBall stance angle

    δDisplacement or contact deformation

    J Moment of inertia

    ωAngular velocity of rotation

    ωRBall or roller track speed

    ωmThe speed at which the ball or roller rotates around its own axis

    ΔψAngular spacing between rolling elements

    1 Introduction

    As a basic industrial component,bearings reflect industrial development [1–3].The bearing preload is a crucial mechanical connection technique that establishes a robust link between the shaft and bearing by applying an appropriate pressure or tension.It effectively minimises clearance,enhances bearing rigidity and precision,reduces operational vibrations and noise,improves the load capacity of the machine,and ensures stable machinery operation[4,5].

    Bearing preloads can be applied in various ways,such as electromagnetic control,hydraulic control,and centrifugal force control [6–9].The mechanical control method offers the benefits of automation and stability,enabling more precise torque and position control while providing a consistent control torque.For instance,Tian et al.[10]used a fixture equipped with a pressure sensor to study the influence of a preload on the nonlinear dynamic response of typical bolt cantilever structures through finite element analysis.The hydraulic preload control method is primarily employed for large structures because of its high accuracy and ability to perform a wide range of torque adjustments.Xu et al.[11] studied the temperature rise in the outer ring of a bearing under varying preload and speed conditions by implementing hydraulic preload control.However,notably,this method entails higher equipment costs and requires specific operational and maintenance requirements.The centrifugal force preload control method offers advantages such as no additional energy consumption and uniform loading.However,it is dependent on the rotating equipment and has limitations in terms of the preload torque that can be achieved.Choi et al.[12] used the centrifugal force and rubber pressure to enhance the preload force and proposed an elbow joint mechanism for accommodating variable preloads.Considering the test conditions and requirements,an enhanced mechanical preloading method was employed for the experiments.In contrast to conventional mechanical preloading methods,which typically involve hardware replacement to modify the preload size,this test allows for the direct real-time adjustment of the preload magnitude while facilitating simultaneous observation and recording.Notably,hardware improvements and the integration of multisensor fusion and fieldprogrammable gate array(FPGA)contribute to improved precision and accuracy of the test[13,14].

    A data-driven approach is one that analyzes and leverages large amounts of data to drive problem-solving and decision-making,and deep learning and neural networks are typical data-driven approaches.With the widespread application of deep learning technology and the development of transfer learning and self-supervised learning,data-driven methods have become more and more important in industrial applications[15].However,it also has limitations such as data bias and limited data generalization ability.

    Multi-sensor fusion is widely used in industrial system monitoring,fault diagnosis,spatial positioning,logistics management and other fields.By fusing data from multiple sensors,the accuracy and reliability of information can be improved,which in turn can achieve more accurate measurement and sensing results,providing important support for efficient data processing and decision-making[16,17].Table 1 describes the mainstream methods of multi-sensor fusion and their advantages and disadvantages.

    Table 1:Mainstream multi-sensor fusion methods

    Lee et al.[18] proposed a centralised fusion algorithm for sensor systems along with a filter algorithm particularly designed for distributed sensor systems.The study concluded that this approach yielded higher effectiveness and accuracy than using a single sensor.Yu et al.[19]proposed a singlesensor network and a multisensor network and conducted a comparative analysis to validate the superiority of the multisensor network.Furthermore,they extensively discussed the optimal placement of multiple sensors.Peng et al.[20]proposed an improved Dempster-Shafer(D-S)synthesis formula based on the concept of evidence information and focus element generated by it and applied it to multisensor information fusion technology.Zhang et al.[21]proposed an asynchronous track fusion algorithm with information feedback,which was combined with a target quality management model for weight allocation to improve the performance of an asynchronous multisensor fusion system.These studies illustrated the feasibility and advantages of multisensor fusion.Therefore,a multisensor fusion method was selected to analyse the bearing preloads.

    In this article,a variety of sensors are used to capture signals such as sound,vibration,temperature,and torque.When the system works after receiving the influence of preload,compared with normal operation,the generated parameters such as sound vibration and temperature rise will change due to the change of force.Torque is different from other signals,the change of torque will change the initial state of the system,and the preload force can be applied to the system in different situations,as a variable group to obtain more data,so that the experimental accuracy and stability are improved.In this paper,the improved D-S evidence theory method is used to process the data,and the results are more accurate and stable than those of the traditional methods.

    To the best of our knowledge,predicting the bearing preload using multisensor fusion is a relatively new method,and the contributions of the current study are as follows:

    1) A new type of precision controllable mechanical preload test bench was established,and the mechanical preload data was measured in various aspects and scales.

    2) An improved multisensor fusion method is proposed to analyse and predict the measured bearing preload and propose a more accurate optimal bearing preload.

    The remainder of this paper is organised as follows:The second section introduces the principle of bearing preload and the influence of mechanical preload,introduces the new controllable preload test bench used in the third section,analyses and predicts the experimental data,and uses deep learning to fit the data to ensure the accuracy of the investigation.Finally,the paper is summarised,and future research is discussed.

    2 Experimental Principle

    2.1 Rolling Bearing Load Analysis

    Various factors generate a dynamic load between the rolling elements and bearing raceways.The dynamic loads are relatively small at medium and low speeds compared with the external forces acting on the balls or rollers.However,at high speeds,dynamic loads such as centrifugal forces and gyroscopic moments significantly influence the load distribution and create dynamic loading conditions [22–24].The centrifugal force Fc is produced when the balls or rollers rotate around the bearing axis.When rolling bearings operate at high speeds,the centrifugal force exerted by balls or rollers becomes substantial.In particular,it can be expressed as:

    Assuming that the preload force on the ball is Qia,and the centrifugal force Fcis applied,the force applied by the ball is as shown in Fig.1.

    When external forces cause rolling ball bearings to rotate around shafts misaligned with the inner or outer ring,the presence of balls inside the bearing induces a gyroscopic effect.This effect gives rise to a gyroscopic moment,denoted as Mg,which acts in a direction opposite to the rotation of the bearing.Consequently,increased friction and wear are observed.The magnitude of gyroscopic moment Mg can be described as follows:

    In high-speed operation,the ball load caused by the combination of centrifugal force and gyroscopic moment when the preload of the ball is Qiais shown in Fig.2.As the ball rolls on the bearing,it is subjected to a balanced force.These include the load,inertia,friction,and elastic deformation forces.At equilibrium,these forces cancel each other out,that is,the sum is zero.The ball force balance equation for rolling bearings is as follows:

    Figure 1:Ball under thrust load and centrifugal force

    Figure 2:Ball load under thrust

    When the centrifugal force acts on the ball,the line of action between the center of curvature of the inner and outer channels becomes a broken line because the contact angles between the ball and the inner and outer channels are different,such as Fig.3.At any azimuth j,the distance between the fixed center of curvature of the outer channel and the final position of the center of the sphere is:

    Because:

    whereδojandδijare the normal contact deformations of the inner and outer raceways,respectively.According to the relative axial displacementδaand the relative angular displacementθof the inner and outer rings,the axial distance and radial distance between the center trajectories of curvature of the inner and outer channels at the free kick position are as follows:

    when a rolling bearing works under a combined load of the axial force Fx,two radial forces Fyand Fz,and two moments Mxand Mz,a relative displacement occurs between the outer and inner rings,which can be expressed as axial displacementδxand radial displacementδy.Assuming that the rolling bearing is not affected by factors such as centrifugal force and gyroscopic moment when the roller is subjected to the maximum load,the total elastic deformation of its inner and outer rings and rollers isδψ,and the formula is as follows.The center angle of the roller isψandα,as shown in Fig.3.

    Figure 3:Displacement caused by axial preload and radial load

    As shown in Eq.(11),the elastomeric variable reaches its maximum when the centre angle of the roller is 0°,and the formula is as follows:

    Combining and simplifying Eqs.(11)and(12)yields Eqs.(13)and(14).

    Maximum ball loads can be obtained as:

    2.2 Optimal Preload Analysis Based on Bearing Performance

    During operation,rolling bearings experience friction and deformation,generating heat.If the heat generated from these sources exceeds the heat dissipation capacity of the bearing,the bearing temperature increases[25–27].Rolling bearings typically have an internal clearance,and applying an axial preload allows for adjusting this internal clearance,thereby enhancing bearing performance.The temperature of the bearing exhibits an approximately linear relationship with both speed and applied axial preload.However,the rate of temperature change decreases to a certain extent,as shown in Fig.4.

    Figure 4:Temperature rise under different rotating speed and preload

    Typically,the maximum allowable service temperature for the surface of a rolling bearing does not exceed 80°C.In addition,the temperature rise should not exceed 40°C above the ambient temperature.By continuously monitoring the maximum surface temperature and temperature rise of the bearing in real-time,we can assess the optimal axial bearing preload using vibration and noise analyses at various speeds.Subsequently,a variable-axial-preload device was used to adjust the preload to its optimal value.This approach enables the efficient control and adjustment of the bearing preload,thereby ensuring optimal performance.

    Bearing vibration refers to the deviation of bearing elements from their ideal positions.If a sufficient preload level is not achieved,the balls may undergo intermittent slipping and rolling,leading to instability.In such cases,the vibration level can be significantly higher (by one or two orders of magnitude) than typical bearing vibrations.In addition,the axial preload influences the damping characteristics of the spindle assembly,further affecting its performance.

    The vibration characteristics of rolling bearings vary at different speeds and are not directly correlated with the temperature increase.Bearing vibrations originate from the inherent factors within the bearing.It does not exhibit a linear increase with speed;instead,it reaches its peak vibration level at a specific critical speed [28–30].However,the vibrations induced by the axial preload were less pronounced within certain speed ranges,as shown in Fig.5.Therefore,reducing the axial preload is necessary when the bearing speed approaches the critical speed to minimise bearing and system vibrations.

    Figure 5:Vibration level under different rotating speeds and preload

    Conversely,when the influence of the axial preload is insignificant,increasing the preload is essential to ensure the stiffness of the bearing and the overall system.Bearing noise is generated by mechanical waves that cause air vibrations due to the vibration of the bearings.Even if individual bearing components are perfectly machined,the relative movements between the inner and outer rings,cages,and balls cannot be eliminated entirely during operation.

    In general,the axial preload induced a slight change in the bearing noise,as listed in Table 2.Applying an axial preload to the bearing increases its stiffness,thereby significantly influencing the bearing and entire system[31].Optimal axial preload applications can effectively mitigate issues such as acoustic noise.

    Table 2:Noise level at different speeds(dB)

    The key advantage of multisensor fusion technology is its ability to enhance the accuracy and perception capabilities of sensing systems.By combining the observation information acquired from multiple sensors at different levels and perspectives,we can extract more valuable insights.This,in turn,improves the decision-making reliability,confidence,and environmental adaptability of the system.Through applying multisensor fusion,we gained a deeper understanding of how the bearing preload is influenced by multiple factors,allowing us to determine the optimal preload value comprehensively.

    3 Controlled Preload Test Bench Design

    Preload is the most important part of an experiment,and the accuracy of the data must be guaranteed.This test bench was designed as a closed test bench that can be adjusted by a servomotor,and the test used a multisensor fusion method to measure the controllable preload.In the experiment,the sensor and FPGA data interacted with each other so that the test could obtain higher precision,and the data could also be fed back to the host computer in real time.

    3.1 Test Bench Design

    An appropriate preload is important for system operation and bearing life,among other factors.Determining the optimal preload force involves considering various factors such as noise,vibration,and temperature.This study developed a comprehensive test bench to measure and analyse all the relevant factors and provide real-time feedback.Herein,the test bench consisted of several key components,including a host computer,motor control system,torque loading unit,preload control unit,and data acquisition module(Fig.6).AC motors,servomotors,and stepper motors were used in this setup.

    Figure 6:Test bench for optimum preload of controllable bearing

    In particular,AC motors provided the drive to the system,servomotors were used to control the preload force,and stepper motors generated torque by acting on the shaft.The FPGA served as the core controller for the entire system.It facilitated the synchronisation of the eight-channel analogueto-digital(AD)conversion,enabling the conversion of changing analogue voltages from specified IO ports into numerical values.This process was an integral part of the signal sampling procedure.Using this systematic configuration and employing advanced technologies,the experiment ensured accurate measurement,analysis,and control of the preload,thereby enhancing our understanding of its impact on the system performance,bearing life,and other relevant aspects.

    3.2 Closed System Design

    A closed test bench,also known as a power flow closed test bench is primarily used to evaluate the performance of mechanical transmission components and their systems in an enclosed environment.Its purpose is to determine the service life of each component,assess the work efficiency,and identify the impact of factors such as materials and loads.A closed test bench offers several advantages,including a simple structure and low cost,low energy consumption,strong economic applicability,and a wide range of applications.Closed test benches can be classified based on the closed methods employed,with the main types being electrically closed,hydraulically closed,and mechanically closed.A mechanically closed method was adopted for the test bench.

    The test bench shown in Fig.7 consists of two pairs of gears.These gears were meshed together with z1,z4,z2,and z3,which were connected by shafts.The loading device was positioned on the right side of gear z4,and the torque was transmitted through gap matching between the gears,shafts,and system.Shafts 1 and 2 in the entire system,spanning from z3 to z4,received the torque retained within the system.This torque represents the load added to the system,known as the closed torque.Different closed powers were obtained by adjusting the torque and rotational speed.

    Figure 7:Closed system structure

    The loading device in this test bench primarily consisted of a loading gear and a motor.The gear and shaft were clearance-matched,facilitating the convenient adjustment of the applied load on the shaft.Although the system is closed,when functioning,it will cause unequal energy loss due to friction and other factors.These losses will be compensated by the drive device;therefore,the power output by the drive device is only used to make the system run,which can also compensate for the power lost by friction in the system.This part of the power is generally small;hence,the energy output required by the drive device is not significant.

    3.3 Dynamic Torque Loading Design

    The dynamic torque transducer employed in this test bench enabled torque measurement in both the positive and negative directions through noncontact power supply and signal output methods.The test bench incorporated a designed dynamic torque loading system,where the torque served as a control variable owing to its influence on the optimal preload.Before the commencement of the test bench system operation,the loading system applied torque to the entire system.Through closed transmission within the entire system,the power-flow direction was determined by the rotation direction of the loading and driving devices.A significant portion of the power of the active gear was transmitted to the driven gear drive system,which generated a force that sequentially propagated from the active gear to the driven gear.Eventually,z3 and z4 produced torques of the same magnitude but in opposite directions owing to their opposing rotation directions.Consequently,the system rotated under the influence of this torque.This validates the operation of the entire system,including the presence of torque,impact of preload,and other factors.

    Various closed powers were achieved on the test bench by manipulating the torque and rotational speed.The loading device primarily consisted of a loading gear and motor.The gear and shaft were designed with appropriate clearance,facilitating the convenient adjustment of the applied load on the shaft.

    A stepper motor was employed to control the system and apply torque to facilitate the experimentation.This allows for the real-time adjustment of the torque within the system.Furthermore,the motor was equipped with a self-locking mechanism that ensured the torque remained locked and did not change during operation.This feature enhanced the stability and consistency of the system during testing.

    3.4 Controlled Preload Design

    On this test bench,the application of a preload was achieved indirectly by rotating a threaded nut driven by a bevel gearbox.The motor drove the bevel gearbox to rotate at a specific angle via coupling.Within the bevel gearbox,the bevel gears meshed with each other,causing the output shaft to rotate.The front side of the output shaft was threaded and equipped with a suitably sized nut.This nut was secured on both sides by ferrules,preventing it from rotating circumferentially(Fig.8).

    Figure 8:Three-dimensional diagram of the test bench and explosion diagram of preload

    When the test shaft was rotated,the threaded portion at the front of the shaft caused the fixed nut to move horizontally along the threaded shaft.This horizontal movement of the nut occurred in the direction of the measured bearing.Consequently,the nut transmitted a force to the bearing through the top sleeve and pressure sensor.The current preload was accurately determined by collecting the force received by the pressure sensor.

    This method of obtaining preload offered both stability and controllability.The threaded shaft remained fixed owing to motor fixation and did not rotate during the measurement process.After the nut moved forward,it remained fixed by a ferrule.Consequently,both the circumference of the nut and threaded shaft remain fixed,ensuring that the position of the nut and magnitude of the preload remain constant.

    3.5 Test Bench Measurement and Control System Design

    The measurement and control systems of the controllable bearing optimal preload test bench were divided into hardware and software designs.The hardware design primarily included the core controller FPGA,various sensors,and a host computer.This part primarily realised a controllable preload and data acquisition and transmission of each sensor.The software design displayed the measured value and optimal preload of each sensor online.

    The hardware design of the measurement and control system primarily comprised a core controller field-programmable gate array(FPGA),pressure sensor,temperature rise sensor,noise sensor,vibration acceleration sensor,torque speed sensor,host computer,motor,and controller.A multisensor fusion architecture with two inputs and three outputs was constructed to determine the optimal preload for the best bearing performance(Fig.9).

    Figure 9:Hardware architecture of measurement and control system

    The hardware of the measurement and control system can be divided into two modules: the multisensor fusion module and optimal preload analysis and control module.The multisensor fusion module was primarily based on the AN706 eight-channel synchronous AD acquisition module.This module enabled the simultaneous collection of measurements from five different sensors.

    The sampling frequency of the system can reach 200 MHz with parallel control of a highprecision,low-timing FPGA.This allowed for high-speed and low-error voltage acquisition,ensuring synchronisation of the data from all five sensors.The architecture of the multisensor fusion module incorporated the pressure and torque speed sensors as two inputs,whereas the temperature rise,noise,and vibration acceleration sensors served as three outputs.

    All five sensors,either through their own voltage output or conversion by the corresponding transmitter,produced a voltage signal within the 0–5 V range.Therefore,multichannel AD acquisition was highly suitable for this setup.Among these sensors,the pressure,temperature rise,noise,and torque speed sensors all provided a voltage output within the 0–5 V range according to their respective measurement ranges.The voltage output of the vibration accelerometer increased with the vibration intensity but did not exceed 5 V.The data collected from the five sensors and FPGA interacted using an SPI bus,ensuring synchronisation and efficient data transfer.

    In the optimal preload analysis and control module,the field-programmable gate array(FPGA)served as the core controller responsible for data collection and real-time monitoring of preload values using pressure sensors.This allowed the determination of the current preload force through multisensor fusion.If the current preload force was not optimal,an iterative process was initiated to update the preload force gradually.

    The preload force was controlled by the angle of rotation of the servomotor.A larger angle corresponded to a greater preload force.The required angle of the servomotor was determined by the number of pulses generated by the servo controller.By writing a Verilog program in advance,the FPGA incrementally increased the pulse output of the corresponding pin to control the iterative update of the preload force.By continuously adjusting the pulse output,the FPGA enabled the servomotor to rotate incrementally,thereby increasing the preload force in small increments.This process facilitated the collection of new data and the subsequent evaluation of the optimal preload force.An iterative update of the preload value ensured that the system could dynamically adjust and maintain optimal preload conditions.

    4 Experimental Setup and Result Analysis

    4.1 Feature Extraction

    Feature extraction was performed after data preprocessing,which is the process of converting the original signal into a meaningful feature vector.This process eliminated unnecessary noise and interference as well as performed data normalisation and standardisation,enabling more effective learning and classification by neural networks.

    The vibration signal was captured using a vibration acceleration sensor that encompassed time,frequency,and time-frequency domain characteristics.The time-domain feature is a fundamental characteristic that provides insight into the overall trend and distribution of the signal.The frequencydomain characteristics reveal the frequency and energy distributions of the signal.Time-frequency domain characteristics depict both the time and frequency domain features of the signal.In this study,the time-frequency domain analysis method was employed to extract features such as the frequency range,amplitude,and time range,as depicted in Fig.10.It is evident that the application of preload results in a significant increase in vibration levels.

    The noise signal was captured using a noise sensor microphone,encompassing various characteristics,such as the sound pressure level,frequency characteristics,spectral analysis,and timedomain characteristics.The sound pressure level serves as a fundamental feature that reflects noise intensity.The frequency characteristics depict the energy distributions of the different frequency components within the noise signal.Spectral analysis enables the acquisition of power spectral density and frequency distribution of noisy signals.The time-domain characteristics reveal the waveform and statistical features of the noisy signal.

    Herein,a combination of frequency-and time-domain analyses was employed for the feature extraction of noise signals.Particularly,the time-domain characteristics of the noise signal were obtained through a time-domain analysis.Subsequently,a spectral analysis was conducted to obtain the frequency-domain characteristics of the noise signal.Fig.11 illustrates the noise signal and its frequency-domain representation under no preload force and with a preload force of 100 N.The application of a 100 N preload force significantly reduced the noise level.Unlike vibration and noise signals,temperature rise signals primarily reflect the operational and thermal equilibrium states of rolling bearings.The characteristics of the temperature rise signals include the steady-state temperature,temperature change rate,and temperature gradient.The feature extraction of temperature rise signals typically involves a combination of traditional statistical methods and machine learning techniques.

    Figure 10:Vibration signals and frequency domain signals under no preload and 100 N preload

    This study employed statistical functions for feature extraction of the temperature rise signals.The statistical features of the temperature increase signal were extracted,and Fig.12 shows the temperature increase signal under no preload force and with a preload force of 100 N.When no preload force was present,the change rate,mean,maximum,and minimum values were 0.005,21.935,21.301,and 21.612,respectively.Conversely,when a preload force of 100 N was applied,the change rate,mean,maximum value,and minimum value were 0.003,22.683,22.296,and 22.462,respectively.Therefore,applying the preload force significantly inhibited the rate of change in the temperature increase.

    4.2 D–S Evidence Theory

    All evidence does not have the same confidence level in the fusion of information and decisionmaking processes.Therefore,weighted evidence can be obtained by introducing the credibility of evidence,weighted by the credibility of evidence,weighted by each evidence frame,and a weighted sum to obtain the weighted evidence,thereby solving this problem.The calculation formula is as follows:

    wherem(A)is the weighted evidence obtained,wiis the feasibility of each piece of sensor evidence,m(ai)is the original evidence,and n is the number of sensors used.Ifm(A)>0,A is the focal element of the evidence,and the set of all focal elements is called the nucleus.

    Figure 11:Noise signal and frequency domain signal under no preload and 100 N preload

    Figure 12:Temperature rise signal under no preload and 100 N preload

    The comprehensive rule of the D–S evidence theory is used to integrate the three sets of weighted evidence,and the formula is:

    when dealing with conflicting information in evidence synthesis using synthetic rules,the Dempster–Shafer (D–S) evidence theory can introduce paradoxes and result in fusion failure.To address the challenge of conflicting evidence,evidence fusion often involves assigning different weights to individual sensors.However,determining the appropriate weights for each sensor can be a complex task,requiring substantial prior knowledge or experience.

    Therefore,this study proposes an enhanced multisensor fusion weighting algorithm to significantly enhance the accuracy of fault diagnosis.Using this algorithm,the fusion process can effectively handle evidence conflicts and improve the reliability of diagnostic outcomes.

    4.3 Improved Multisensor Fusion Weighting Algorithm

    In practical scenarios,different sensors exhibit varying levels of reliability,thereby influencing confidence in evidence synthesis.For an information fusion system comprising n sensors,the weight assigned to each sensor can be determined through reliability testing and simulation analyses.These weights are denoted as,where the superscript“pri”indicates the prior weight.

    The prior weight represents the relative reliability of each sensor based on the aforementioned analyses.However,the actual working environment introduces uncertainties,necessitating the calculation of evidence similarity to obtain real-time sensor reliability,referred to as the posterior weight.A novel method for evidence synthesis can be derived by considering both prior and posterior weights[32].

    By incorporating prior weights that reflect the initial reliability assessment and posterior weights obtained through evidence similarity calculations,this new methodology facilitates more accurate and adaptive evidence synthesis in dynamic environments.

    Let there be n incompatible pieces of evidence,where the basic trust distribution functions of any two pieces of evidenceEiandEjaremiandmj,and the focal elements areAiandAj,respectively.The distance between evidencemiandmjis expressed as:

    whereMi=,D is a matrix of 2n×2n,and the elements in the matrix are as follows:

    Combining Eqs.(20)and(21)yields,and the distance betweenmJcan be expressed as:

    Define the body of evidencemi,the similarity coefficient betweenmjis:

    Each evidence available supports theEias:

    Normalization is available as:

    Eq.(26) reflects the confidence of evidenceEiand satisfies=1,therefore,it can be a posterior weight.If the evidence is similar,they are highly supportive of each other;therefore,the greater the confidence,the greater the posterior weight.Considering both the a priori and posterior weights,the composite weights of the computational evidenceEiare:

    Corresponding signal processing methods were employed to extract the features from these sensors.To minimise the errors,the characteristic signals from each sensor were averaged.Subsequently,a multisensor-weighted fusion approach was employed using the aforementioned algorithm.This fusion technique measured the mean,variance,and standard deviation of the weighted signal to determine the optimal preload force required to achieve the maximum bearing performance(Fig.13).The smaller the parameters obtained after weighted fusion,the closer it is proved to be to the optimal preload of the bearing.We can get a more correct bearing preload through experiments.Compared with the traditional method of bearing preload determination,the method used in this paper uses multiple sensors and signals,and the results obtained are more stable and reliable,and the robustness is stronger,and will not have a serious impact on the whole system due to problems caused by a certain set of data.At the same time,sensors have different response capabilities to target information,and may give incomplete or indeterminate information.In order to improve the reliability and credibility of the information fusion results of multi-sensor systems,this paper obtains the posteriori weights of the evidence by calculating the similarity matrix between the evidences,and combines the prior weights pre-assigned according to the inherent reliability of the sensors to obtain the composite weights of the evidence body,and then weights the weighted average of the evidence,and finally uses the D-S synthesis formula to fuse the evidence.

    Figure 13:Sensor weighted fusion

    In the controllable preload test bench developed in this study,the temperature increase was monitored within a range of 25°C through orthogonal experiments involving various speeds and torques.Three sets of signals,namely,noise,vibration,and temperature rise,were collected from multiple sensors.These three sets of signals will reflect the changes of parameters such as preload torque when the system is working,and these three sets of data will be analyzed,and the signals before and after the initial conditions of the system are changed,and they will be fused.From the results shown in Fig.14,it can be observed that,under low-speed operating conditions for rolling bearings,applying a large torque is favourable for enhancing the bearing stiffness.Conversely,when operating at high speeds,a smaller torque is more suitable for reducing the performance degradation caused by a temperature increase.By applying a certain torque,the preload force was reduced to a certain extent.

    Compared with the latest articles in the field of data fusion [33,34],Pan uses an improved stochastic weighting algorithm to effectively deal with rolling bearing fault diagnosis,but its use has certain limitations,due to the large difference in vibration frequency between different bearing faults,its effectiveness will be limited when dealing with mixed faults or other faults with small frequency differences.Moreover,the stochastic weighting algorithm itself still has problems such as high requirements for samples and parameters,easy to be disturbed by abnormal data,and unable to deal with non-stationary environments,which causes difficulties in application.The method used in this paper is more flexible,adaptable and robust,which can better deal with uncertainty and conflicting information,and provide more accurate and reliable fusion results.It also better handles high-conflict evidence and converges faster.The results obtained in this paper are shown in Fig.14 and have been verified to be more than 99%accurate.In summary,the multi-sensor fusion method proposed in this paper has certain advantages over the traditional method,and can obtain more stable and accurate results.

    Figure 14:Optimum preload at different speeds and torques

    4.4 Neural Network Fitting

    Fitting a neural network is the process of adjusting the parameters of a neural network model so that it can learn from the input data and approximate the desired output [35].Neural networks are composed of multiple levels of nonlinear units that are capable of learning and expressing complex nonlinear patterns,so they can fit a variety of complex data patterns well.The fitting ability of a neural network allows it to be used to predict data,classify and identify.The accuracy of the results obtained in this paper is determined by fitting the best preload obtained in this paper.

    The neural network fitting process involves the iterative optimisation of the activation function coefficients for each neuron,creating a progressive improvement in performance.In this study,a Simple Neural Network (SNN) model was used because it offers optimal performance with a straightforward structure.The SNN model comprised a hidden layer containing multiple neurones and an output layer consisting of a single neuron.The topology of the complete SNN model is shown in Fig.15.

    For each neuron in the hidden layer of the SNN model,the output can be expressed as:

    In the SNN model,the weights and biases of each neuron are denoted asωiand bi,respectively.The input dataxicomprise ofx1,x2,...,andxn,whereasyirepresents the output neuron of the i-th neuron,and N denotes the total number of neurones in the hidden layer.The input signal xi undergoes nonlinear transformation through the neurones in the hidden layer to generate the output signalyi.This transformation is achieved using the activation functionσ(x),which is carefully chosen to ensure stable signal propagation and ease of differentiation.This can be mathematically expressed as follows:

    According to the topology of the SNN model,we can obtain the output of the SNN model as:

    Figure 15:Neural network topology

    The weights and biases of the output layer are denoted W and B,respectively.The output signal yifrom each neuron in the hidden layer is passed to the output layer.Through a linear transformation,the SNN model produced an output value of Y.The activation function of the output layer was linear.During training,mean squared error(MSE)determined the convergence of the function.The minimum MSE corresponded to the optimal curve fit of the dataset.To achieve this,the weights and biases of the neurones in the hidden and output layers were adjusted.The mean squared error(MSE)between the output value Y of the SNN model and test data was compared,and iterative training using the error descent method was performed along the gradient direction.The procedure was designed to determine the model parameters that resulted in minimal MSE.

    The optimal preload data obtained from the measurements were selected for fitting using deeplearning techniques.The fitting results are shown in Fig.16,which illustrates the regression between the measured and output data of the SNN model.An approximately linear relationship exists between the measured and output data,with an R-squared value of 0.99.This high R-squared value indicates an excellent fit and substantiates the effectiveness of the SNN model.

    Figure 16:Fit regression

    5 Conclusion

    This study proposes a novel model for controlling the preload force in bearings.This includes developing a controllable bearing preload experimental platform,measurement and control system based on an FPGA,mechanical power closed-loop design,and dynamic torque loading design.Moreover,a weighted fusion analysis method for determining the optimal preload force was introduced,which applied the D–S evidence theory in a multisensor setting.The optimal preload under different conditions obtained by algorithm fusion is shown in Fig.14.The accuracy of the experimental results was validated by fitting the results to neural networks.The method used in this paper is widely used and has no obvious shortcomings.It can better handle high-conflict evidence,converge faster,and has distinct advantages over traditional methods,resulting in more stable and accurate results.

    Acknowledgement:The authors thank the support from the Key Project of National Natural Science Foundation of China,the Open Project of State Key Laboratory of Mining Response and Disaster Prevention and Control in Deep Coal Mines,the National Key Research and Development Program of China,the Open Project of Key Laboratory of Conveyance and Equipment,the Science and Technology Research Project of Jiangxi Provincial Department of Education,the Supply and Demand Linking Employment Education Project of the Ministry of Education,Suzhou Basic Research Project,and the Open Project of State Key Laboratory for Manufacturing Systems Engineering.

    Funding Statement:This research is supported by: The Key Project of National Natural Science Foundation of China(U21A20125);The Open Project of State Key Laboratory of Mining Response and Disaster Prevention and Control in Deep Coal Mines (SKLMRDPC21KF03);The National Key Research and Development Program of China (2020YFB1314203,2020YFB1314103);The Open Project of Key Laboratory of Conveyance and Equipment (KLCE2021-05);The Science and Technology Research Project of Jiangxi Provincial Department of Education(GJJ210639);The Supply and Demand Linking Employment Education Project of the Ministry of Education (20220100621);The Open Project of State Key Laboratory for Manufacturing Systems Engineering(sklms2023009);The Suzhou Basic Research Project(SJC2023003)(Corresponding author:Yasheng Chang.).

    Author Contributions:The authors confirm contribution to the paper as follows:study conception and design:Kuosheng Jiang,Yasheng Chang;data collection:Chengrui Han;analysis and interpretation of results: Kuosheng Jiang,Chengrui Han,Yasheng Chang;draft manuscript preparation: Chengrui Han,Yasheng Chang.All authors reviewed the results and approved the final version of the manuscript.

    Availability of Data and Materials:The datasets used or analysed during the current study are available from the corresponding author on reasonable request.

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    午夜激情久久久久久久| 国产极品天堂在线| 天天操日日干夜夜撸| 黄频高清免费视频| 欧美在线一区亚洲| 肉色欧美久久久久久久蜜桃| av.在线天堂| 国产精品 国内视频| 亚洲精品自拍成人| 国产黄频视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 日本wwww免费看| 桃花免费在线播放| 国产亚洲最大av| 丰满少妇做爰视频| 九九爱精品视频在线观看| 熟女av电影| 美女大奶头黄色视频| 少妇人妻久久综合中文| av电影中文网址| 熟妇人妻不卡中文字幕| 国产伦理片在线播放av一区| 最新的欧美精品一区二区| 一本色道久久久久久精品综合| 91精品伊人久久大香线蕉| 飞空精品影院首页| 黄色视频不卡| 久久精品久久久久久噜噜老黄| 久久久久久久久免费视频了| 晚上一个人看的免费电影| 在线精品无人区一区二区三| 69精品国产乱码久久久| 黄色怎么调成土黄色| 亚洲欧美日韩另类电影网站| 一本一本久久a久久精品综合妖精| 国产精品 国内视频| 熟女少妇亚洲综合色aaa.| 欧美国产精品va在线观看不卡| 午夜福利视频在线观看免费| 丝袜美腿诱惑在线| 成人三级做爰电影| 亚洲欧美日韩另类电影网站| 热99久久久久精品小说推荐| 街头女战士在线观看网站| 国产日韩一区二区三区精品不卡| 久久人人爽人人片av| 亚洲人成网站在线观看播放| 91aial.com中文字幕在线观看| 久久影院123| 亚洲欧美日韩另类电影网站| 国产在线视频一区二区| 制服诱惑二区| 丁香六月欧美| 丰满乱子伦码专区| 夜夜骑夜夜射夜夜干| 精品亚洲乱码少妇综合久久| 大话2 男鬼变身卡| √禁漫天堂资源中文www| 亚洲欧美中文字幕日韩二区| 亚洲国产看品久久| 国产女主播在线喷水免费视频网站| 国产一区亚洲一区在线观看| 五月天丁香电影| 久久国产精品男人的天堂亚洲| 亚洲第一区二区三区不卡| 一本—道久久a久久精品蜜桃钙片| 99久久99久久久精品蜜桃| 伊人久久大香线蕉亚洲五| 午夜福利,免费看| 久久99精品国语久久久| 91精品伊人久久大香线蕉| 日韩制服丝袜自拍偷拍| 日韩熟女老妇一区二区性免费视频| 国产精品一国产av| 啦啦啦在线免费观看视频4| 亚洲精品国产区一区二| 黄色视频不卡| 可以免费在线观看a视频的电影网站 | 久久久国产一区二区| 中文乱码字字幕精品一区二区三区| tube8黄色片| 国产成人精品久久久久久| 可以免费在线观看a视频的电影网站 | 精品人妻熟女毛片av久久网站| 国产精品一国产av| 男男h啪啪无遮挡| 国产在线一区二区三区精| 国产男女超爽视频在线观看| 色吧在线观看| 一级a爱视频在线免费观看| 80岁老熟妇乱子伦牲交| 一区在线观看完整版| 一区二区三区四区激情视频| 国产在线一区二区三区精| 欧美日韩亚洲综合一区二区三区_| 亚洲一级一片aⅴ在线观看| 99久国产av精品国产电影| 日韩电影二区| 黄片无遮挡物在线观看| 久久精品国产综合久久久| 日本91视频免费播放| 国产有黄有色有爽视频| av又黄又爽大尺度在线免费看| 在现免费观看毛片| 午夜福利一区二区在线看| 十八禁人妻一区二区| 色视频在线一区二区三区| 欧美日韩亚洲高清精品| 国产av精品麻豆| 狂野欧美激情性xxxx| 老司机亚洲免费影院| 久久婷婷青草| 免费观看a级毛片全部| 丁香六月欧美| avwww免费| 亚洲国产欧美一区二区综合| 日日摸夜夜添夜夜爱| 亚洲伊人色综图| 亚洲精品一二三| 在线精品无人区一区二区三| 精品酒店卫生间| 伊人亚洲综合成人网| 日本wwww免费看| av一本久久久久| 国产免费一区二区三区四区乱码| 最近的中文字幕免费完整| 国产深夜福利视频在线观看| 捣出白浆h1v1| 男的添女的下面高潮视频| 日本av免费视频播放| 日韩大片免费观看网站| 制服人妻中文乱码| 亚洲自偷自拍图片 自拍| 国产精品一二三区在线看| 免费黄频网站在线观看国产| 国产毛片在线视频| 黄频高清免费视频| 亚洲欧洲国产日韩| 黄片播放在线免费| 国产 精品1| h视频一区二区三区| 男女免费视频国产| 国产精品人妻久久久影院| 免费观看性生交大片5| a级毛片黄视频| 久久人人97超碰香蕉20202| 亚洲欧美色中文字幕在线| 岛国毛片在线播放| 亚洲一区二区三区欧美精品| 性少妇av在线| 91老司机精品| 九色亚洲精品在线播放| 女人高潮潮喷娇喘18禁视频| 亚洲国产精品999| 国产成人精品无人区| 欧美人与性动交α欧美精品济南到| 丰满乱子伦码专区| 老司机影院成人| 亚洲国产欧美日韩在线播放| 亚洲精华国产精华液的使用体验| 久久99一区二区三区| 一本色道久久久久久精品综合| 国产亚洲最大av| 中文天堂在线官网| 又大又黄又爽视频免费| 国产男女内射视频| 国产在视频线精品| 国产精品蜜桃在线观看| 亚洲人成网站在线观看播放| 国产精品欧美亚洲77777| 久久青草综合色| 欧美乱码精品一区二区三区| 在线免费观看不下载黄p国产| 极品少妇高潮喷水抽搐| 亚洲成av片中文字幕在线观看| 18在线观看网站| 亚洲欧洲精品一区二区精品久久久 | 欧美亚洲 丝袜 人妻 在线| 女人精品久久久久毛片| 久久久精品94久久精品| 王馨瑶露胸无遮挡在线观看| 国产精品久久久久久精品古装| avwww免费| 国产成人精品无人区| 狂野欧美激情性xxxx| 99久久综合免费| 一边摸一边抽搐一进一出视频| 老汉色∧v一级毛片| 日韩一区二区视频免费看| 免费高清在线观看日韩| 一区在线观看完整版| 久久午夜综合久久蜜桃| 777久久人妻少妇嫩草av网站| 丁香六月天网| 夫妻性生交免费视频一级片| 亚洲天堂av无毛| 黑人欧美特级aaaaaa片| 97人妻天天添夜夜摸| 免费少妇av软件| 欧美激情高清一区二区三区 | 国产精品免费视频内射| 亚洲精品成人av观看孕妇| 国产又爽黄色视频| 亚洲精品一二三| 美女大奶头黄色视频| xxxhd国产人妻xxx| 久热这里只有精品99| 免费在线观看黄色视频的| 亚洲精品日本国产第一区| 亚洲在久久综合| 精品午夜福利在线看| 啦啦啦啦在线视频资源| av网站免费在线观看视频| 91精品国产国语对白视频| 精品一品国产午夜福利视频| 亚洲欧美成人综合另类久久久| 美女国产高潮福利片在线看| 日韩熟女老妇一区二区性免费视频| 国产精品 欧美亚洲| 熟妇人妻不卡中文字幕| 男女边摸边吃奶| 男男h啪啪无遮挡| 一区二区日韩欧美中文字幕| 少妇的丰满在线观看| 一区二区三区激情视频| 波野结衣二区三区在线| 观看美女的网站| 天天躁狠狠躁夜夜躁狠狠躁| 最黄视频免费看| 9热在线视频观看99| 黄色一级大片看看| av天堂久久9| 十八禁高潮呻吟视频| 18在线观看网站| 久久韩国三级中文字幕| 如何舔出高潮| 久久久精品94久久精品| 91精品伊人久久大香线蕉| 亚洲国产精品成人久久小说| 欧美精品高潮呻吟av久久| 纵有疾风起免费观看全集完整版| 自拍欧美九色日韩亚洲蝌蚪91| 老司机影院成人| 国产又色又爽无遮挡免| 18在线观看网站| 久久青草综合色| av在线app专区| 制服丝袜香蕉在线| 18禁裸乳无遮挡动漫免费视频| 国产免费视频播放在线视频| 老熟女久久久| 中文字幕另类日韩欧美亚洲嫩草| 久久久精品国产亚洲av高清涩受| 亚洲情色 制服丝袜| 亚洲欧美激情在线| 在线精品无人区一区二区三| 亚洲天堂av无毛| 天美传媒精品一区二区| 韩国av在线不卡| 亚洲人成网站在线观看播放| 永久免费av网站大全| 亚洲综合色网址| 国产精品熟女久久久久浪| 人人妻,人人澡人人爽秒播 | 欧美国产精品va在线观看不卡| 欧美97在线视频| 国产免费福利视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 男女午夜视频在线观看| 国产熟女欧美一区二区| 亚洲一级一片aⅴ在线观看| 黄色 视频免费看| 最近最新中文字幕大全免费视频 | 国产成人啪精品午夜网站| 国产精品久久久久久精品电影小说| 美女中出高潮动态图| 国产成人精品无人区| 操美女的视频在线观看| 国产 精品1| 欧美黄色片欧美黄色片| 汤姆久久久久久久影院中文字幕| 青青草视频在线视频观看| 悠悠久久av| 午夜福利在线免费观看网站| 99久久精品国产亚洲精品| 我要看黄色一级片免费的| 国产精品久久久av美女十八| 男的添女的下面高潮视频| 国产男女内射视频| 秋霞在线观看毛片| 深夜精品福利| 免费在线观看黄色视频的| 亚洲国产最新在线播放| 天天影视国产精品| 色精品久久人妻99蜜桃| 青春草亚洲视频在线观看| 亚洲国产欧美日韩在线播放| 人体艺术视频欧美日本| 中文字幕色久视频| 波多野结衣一区麻豆| www.精华液| 丝袜脚勾引网站| 国产成人一区二区在线| 国产1区2区3区精品| 精品人妻一区二区三区麻豆| 国产精品久久久久久人妻精品电影 | 青春草视频在线免费观看| 男女下面插进去视频免费观看| 精品第一国产精品| 超碰97精品在线观看| 啦啦啦 在线观看视频| 桃花免费在线播放| 五月天丁香电影| 宅男免费午夜| 久久久久精品久久久久真实原创| 无限看片的www在线观看| 日韩视频在线欧美| 国产精品国产三级专区第一集| 我的亚洲天堂| av在线app专区| 1024视频免费在线观看| 国产日韩一区二区三区精品不卡| 久久久久久人妻| 在线天堂最新版资源| 久久久久国产一级毛片高清牌| 性色av一级| 亚洲中文av在线| 日日撸夜夜添| 免费日韩欧美在线观看| 午夜免费男女啪啪视频观看| 精品一区二区免费观看| 国产乱人偷精品视频| 美女主播在线视频| 中文字幕制服av| 波多野结衣一区麻豆| 激情视频va一区二区三区| 最新的欧美精品一区二区| 国产黄色视频一区二区在线观看| 伊人亚洲综合成人网| 亚洲精品日本国产第一区| 黄频高清免费视频| 五月天丁香电影| 亚洲欧美精品自产自拍| 亚洲欧美一区二区三区久久| 美女脱内裤让男人舔精品视频| 热re99久久国产66热| 国产探花极品一区二区| 精品少妇久久久久久888优播| kizo精华| 国产人伦9x9x在线观看| 亚洲精华国产精华液的使用体验| www.熟女人妻精品国产| 热99久久久久精品小说推荐| 久久久久人妻精品一区果冻| 赤兔流量卡办理| 亚洲欧洲日产国产| 一区福利在线观看| 777米奇影视久久| 亚洲国产毛片av蜜桃av| 一级毛片 在线播放| 九九爱精品视频在线观看| 日本av免费视频播放| 麻豆精品久久久久久蜜桃| 色精品久久人妻99蜜桃| 免费黄色在线免费观看| 汤姆久久久久久久影院中文字幕| 亚洲五月色婷婷综合| 母亲3免费完整高清在线观看| 2021少妇久久久久久久久久久| 国产成人精品在线电影| 亚洲欧美一区二区三区久久| 天天添夜夜摸| 亚洲欧美清纯卡通| 人体艺术视频欧美日本| 亚洲欧美清纯卡通| 国产免费福利视频在线观看| 久久久久网色| 国产成人精品久久二区二区91 | 国产精品国产av在线观看| 国产精品三级大全| 秋霞伦理黄片| 日韩制服骚丝袜av| 国产亚洲av高清不卡| 一个人免费看片子| 熟女av电影| 欧美人与性动交α欧美精品济南到| 一区二区三区乱码不卡18| 日韩视频在线欧美| 欧美人与性动交α欧美精品济南到| 中文字幕精品免费在线观看视频| 亚洲,欧美精品.| 操美女的视频在线观看| 日韩中文字幕欧美一区二区 | 天天躁狠狠躁夜夜躁狠狠躁| 麻豆av在线久日| 9色porny在线观看| 国产乱人偷精品视频| 男女边摸边吃奶| 亚洲精品av麻豆狂野| 亚洲成人免费av在线播放| 视频在线观看一区二区三区| 色播在线永久视频| 两个人免费观看高清视频| 又大又黄又爽视频免费| 老司机影院成人| 十八禁网站网址无遮挡| 成人国产av品久久久| 欧美另类一区| 97在线人人人人妻| 久久性视频一级片| 一本大道久久a久久精品| 丝袜美腿诱惑在线| 在线 av 中文字幕| 美女大奶头黄色视频| 最近2019中文字幕mv第一页| 麻豆av在线久日| 一边摸一边做爽爽视频免费| 香蕉国产在线看| 蜜桃国产av成人99| 亚洲天堂av无毛| 纯流量卡能插随身wifi吗| 少妇被粗大猛烈的视频| 亚洲伊人色综图| 女人精品久久久久毛片| 青春草视频在线免费观看| 国产亚洲欧美精品永久| 久久婷婷青草| 久久人人爽av亚洲精品天堂| 欧美日本中文国产一区发布| 男男h啪啪无遮挡| 人人妻人人澡人人爽人人夜夜| 一级毛片我不卡| 黄频高清免费视频| 国产 精品1| 这个男人来自地球电影免费观看 | 交换朋友夫妻互换小说| 97精品久久久久久久久久精品| 一级,二级,三级黄色视频| 在线 av 中文字幕| 色精品久久人妻99蜜桃| 久久久久久免费高清国产稀缺| 最黄视频免费看| 91成人精品电影| av网站在线播放免费| 久久久精品免费免费高清| 在线观看国产h片| 日本黄色日本黄色录像| 国产在视频线精品| 满18在线观看网站| 97精品久久久久久久久久精品| 91aial.com中文字幕在线观看| 免费黄网站久久成人精品| 亚洲成国产人片在线观看| 久热这里只有精品99| 欧美日韩视频高清一区二区三区二| 啦啦啦 在线观看视频| 欧美乱码精品一区二区三区| a级毛片在线看网站| 婷婷成人精品国产| 日韩制服骚丝袜av| 日韩一区二区视频免费看| 国产日韩欧美视频二区| 两个人免费观看高清视频| 久久精品国产亚洲av高清一级| 美女福利国产在线| 男女床上黄色一级片免费看| www.熟女人妻精品国产| 精品视频人人做人人爽| 性少妇av在线| 嫩草影视91久久| 亚洲天堂av无毛| 韩国av在线不卡| 亚洲av欧美aⅴ国产| 无限看片的www在线观看| 欧美最新免费一区二区三区| 一区二区三区四区激情视频| 亚洲成av片中文字幕在线观看| 久久久久人妻精品一区果冻| 80岁老熟妇乱子伦牲交| 国产成人啪精品午夜网站| 国产又爽黄色视频| 各种免费的搞黄视频| 男女午夜视频在线观看| 国产精品秋霞免费鲁丝片| 亚洲专区中文字幕在线 | 亚洲国产欧美网| 啦啦啦 在线观看视频| 久久99热这里只频精品6学生| 18禁国产床啪视频网站| av天堂久久9| 国产精品国产av在线观看| 黄片播放在线免费| 叶爱在线成人免费视频播放| 欧美成人午夜精品| 精品一区二区免费观看| 热99久久久久精品小说推荐| 中文字幕高清在线视频| 男的添女的下面高潮视频| 悠悠久久av| 99香蕉大伊视频| 中文字幕制服av| 涩涩av久久男人的天堂| 中文天堂在线官网| 啦啦啦视频在线资源免费观看| 另类亚洲欧美激情| 成人三级做爰电影| 成年av动漫网址| 男女下面插进去视频免费观看| 天天躁日日躁夜夜躁夜夜| 十八禁网站网址无遮挡| 国产精品女同一区二区软件| 久久久欧美国产精品| 欧美中文综合在线视频| 欧美黑人精品巨大| 精品午夜福利在线看| 美女中出高潮动态图| 夫妻午夜视频| 久久久精品免费免费高清| 久久午夜综合久久蜜桃| 人妻一区二区av| 国产一区二区 视频在线| 国产成人午夜福利电影在线观看| 久久久久人妻精品一区果冻| 久久毛片免费看一区二区三区| 2018国产大陆天天弄谢| 久久国产精品男人的天堂亚洲| 成年女人毛片免费观看观看9 | 亚洲一卡2卡3卡4卡5卡精品中文| 老司机靠b影院| 日韩伦理黄色片| 亚洲美女搞黄在线观看| 久久这里只有精品19| 精品亚洲成a人片在线观看| e午夜精品久久久久久久| 国产野战对白在线观看| 18禁裸乳无遮挡动漫免费视频| 日本欧美视频一区| 亚洲国产中文字幕在线视频| 亚洲精品久久成人aⅴ小说| 国产片内射在线| 国产精品国产av在线观看| 欧美久久黑人一区二区| 久久久久精品人妻al黑| 国产免费视频播放在线视频| 久久女婷五月综合色啪小说| 高清欧美精品videossex| 91国产中文字幕| 国产精品麻豆人妻色哟哟久久| 国产精品女同一区二区软件| 51午夜福利影视在线观看| 日韩中文字幕欧美一区二区 | 亚洲综合精品二区| 免费人妻精品一区二区三区视频| tube8黄色片| xxxhd国产人妻xxx| 亚洲av欧美aⅴ国产| 嫩草影视91久久| 中文字幕最新亚洲高清| 欧美日韩一区二区视频在线观看视频在线| 麻豆精品久久久久久蜜桃| 一本—道久久a久久精品蜜桃钙片| 亚洲成人国产一区在线观看 | 国产深夜福利视频在线观看| 激情五月婷婷亚洲| 成人亚洲欧美一区二区av| 久久久久网色| 免费观看a级毛片全部| 国产一区有黄有色的免费视频| 观看美女的网站| 一级毛片黄色毛片免费观看视频| 国产一级毛片在线| 亚洲欧洲日产国产| 色网站视频免费| 免费看av在线观看网站| 一边亲一边摸免费视频| 久久久久久久久久久免费av| 日韩一区二区三区影片| 日韩中文字幕视频在线看片| 欧美黑人欧美精品刺激| 国产精品熟女久久久久浪| 精品国产乱码久久久久久小说| 考比视频在线观看| 亚洲综合精品二区| 欧美人与性动交α欧美精品济南到| 亚洲精品久久成人aⅴ小说| 9191精品国产免费久久| 国产精品无大码| 日韩精品免费视频一区二区三区| 在线观看免费高清a一片| 亚洲人成电影观看| 欧美av亚洲av综合av国产av | 香蕉丝袜av| 蜜桃国产av成人99| 国产精品欧美亚洲77777| 性高湖久久久久久久久免费观看| 欧美精品亚洲一区二区| 又大又黄又爽视频免费| 亚洲,欧美精品.| 超碰成人久久| 成人国产麻豆网| 久久av网站| 老司机影院成人| 菩萨蛮人人尽说江南好唐韦庄| 高清视频免费观看一区二区| 女性被躁到高潮视频| 一区二区av电影网| 久久av网站| a 毛片基地| 成人国语在线视频| 亚洲久久久国产精品| 欧美 亚洲 国产 日韩一| 视频区图区小说| 欧美激情 高清一区二区三区| 久久久久网色| 操美女的视频在线观看|