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

    Design and fabrication of a full elastic submicron-Newton scale thrust measurement system for plasma micro thrusters

    2021-10-31 08:14:44ZhongkaiZHANG張仲愷GuanrongHANG杭觀榮JiayunQI齊佳運ZunZHANG張尊ZheZHANG章喆JiubinLIU劉久鑌WenjiangYANG楊文將andHaibinTANG湯海濱
    Plasma Science and Technology 2021年10期
    關(guān)鍵詞:海濱

    Zhongkai ZHANG(張仲愷),Guanrong HANG(杭觀榮),Jiayun QI(齊佳運),Zun ZHANG (張尊), Zhe ZHANG (章喆), Jiubin LIU (劉久鑌),Wenjiang YANG (楊文將) and Haibin TANG (湯海濱),5,6

    1 School of Space and Environment, Beihang University, Beijing 100191, People’s Republic of China

    2 Shanghai Institute of Space Propulsion, Shanghai Engineering Research Center of Space Engine,Shanghai 201112, People’s Republic of China

    3 School of Astronautics, Beihang University, Beijing 100191, People’s Republic of China

    4 The Institute of Space Systems, University of Stuttgart, Stuttgart D-70569, Germany

    5 Key Laboratory of Spacecraft Design Optimization & Dynamic Simulation Technologies, Ministry of Education, Beijing 100083, People’s Republic of China

    6 Laboratory of Space Environment Monitoring and Information Processing, Ministry of Industry and Information Technology, Beijing 100083, People’s Republic of China

    Abstract In this work, a force measurement system is proposed to measure the thrust of plasma microthruster with thrust magnitude ranging from sub-micro-Newtons to hundreds micro-Newtons.The thrust measurement system uses an elastic torsional pendulum structure with a capacitance sensor to measure the displacement,which can reflect the position change caused by the applied force perpendicular to the pendulum axis.In the open-loop mode, the steady-state thrust or the impulse of the plasma micro-thruster can be obtained from the swing of the pendulum,and in the closed-loop mode the steady-state thrust can be obtained from the feedback force that keeps the pendulum at a specific position.The thrust respond of the system was calibrated using an electrostatic weak force generation device.Experimental results show that the system can measure a thrust range from 0 to 200 μN in both open-loop mode and closed-loop mode with a thrust resolution of 0.1 μN, and the system can response to a pulse bit at the magnitude of 0.1μN s generated by a micro cathode arc thruster.The background noise of the closed-loop mode is lower than that of the open-loop mode, both less than 0.1/μN Hz in the range of 10 mHz to 5 Hz.

    Keywords: plasma thruster, micro-thrust measurement, torsional pendulum

    1.Introduction

    In recent years, more and more plasma propulsion devices(electric thrusters) have been widely used in various applications such as orbit keeping, drag force compensation, altitude control and adjustment,and active deorbiting[1–3].Due to the high specific impulse and the ability to precisely control micro thrust, electric thrusters, especially micro electric thrusters, have better adaptability than traditional chemical thrusters and cold gas micro thrusters, and can provide more options for diverse space missions[2].Based on the operating principle, electric thrusters can be classified as electrostatic,electromagnetic and electrothermal.Various specific micro plasma thrusters, such as low-power Hall thruster, radio-frequency or microwave ion thrusters, field emission electric propulsion (FEEP), micro cathode arc thruster (μCAT),pulsed plasma thruster (PPT) and other technologies either working in steady state or pulse mode [4–11] are under continuous research.In general, the propellant of a thruster can be gaseous,liquid or solid.When the propellant is heated,ionized and accelerated, it can gain energy to be discharged from the thruster while momentum exchange takes place.Since micro plasma thrusters produce thrusts ranging from sub-micro-Newtons to sub-milli-Newtons, ground-based thrust measurement systems with high thrust resolution and wide range are urgently needed for more accurate measurement of parameters and characteristics of plasma micro thrusters.

    There have been several studies on thrust measurement systems with different ranges and resolutions for thrust ranges in the magnitude of sub-micro-Newtons to mili-Newtons[12].From the classification of measurement principle and device structure,the measurement methods of weak forces of micro-Newton magnitude can be divided into torsional pendulum systems, torsional suspend wires system, single and double pendulum and cantilever type systems, etc to meet different needs of thrust measurement range and resolution.Andrewset al[13] developed a torsional pendulum thrust stand with a gas dynamic calibration device that was capable of resolving thrusts of tens of nano-Newtons to micro-Newtons.Manuelet al[14] used a torsional pendulum measurement device which detected the deflection of the torsion pendulum,and the device was able to measure the impulse thrust of plasma micro thrusters with typical values of a few tens of micro-Newtons, and its control of the system background noise provided the technical support for the Laser Interferometer Space Antenna (LISA) [15, 16] project.A similar torsional pendulum device to measure laser-ablative micro thrusters with the impulse bit of 1–10μN s has been developed by researcher from the University of Tokyo[17].Kurtet al[18]developed a vertical pendulum thrust measurement system and it was capable of measuring micro thrusts of 100μN or more, and its impressive pendulum structure was capable of carrying thrusters with the mass up to 125 kg.Beihang University [19, 20] developed a pulsed thrust measurement system with novel electromagnetic calibration techniques on the magnitude of sub-milli-Newtons earliest in China.Boccalettoet al[21]designed a thrust measurement device with a double-suspending wire structure and measured the thrust of a FEEP thruster, and Yanget al[22] measured the thrust characteristics of a PPT thruster using a similar suspending wire pendulum device in both open and closed operation.Zhouet al[23]developed a torsional pendulum structure with the resolution as low as 24.4 nN,and Soniet al[24]also used a similar torsional pendulum with low spring rate to achieve the resolution of 10 nN.Subhaet al[25] designed a novel compact thrust measurement device to verify a 50μN cold gas thruster using a suspended plate.The Brazilian researchers[26]used a unique electrostatic calibration device for static and dynamic calibration of the torsional pendulum system with tens of micro-Newton in recent years.

    At present,there are relatively few measurement systems applicable to plasma micro thrusters with thrust ranges from sub-micron-Newtons to hundreds micro-Newtons, and most of the existing measurement systems cannot achieve the submicron-Newton resolution and the minimum thrust response for thrusters working in either steady state or pulsed state with low background noise applicable to plasma micro thrusters.According to the requirement of the thrust range and resolution, an elastic micro-Newton thrust balance based on torsional pendulum structure is designed.A capacitance sensor is employed to measure the displacement, which can reflect the position change after the force is imposed.Compared with other types of micro-Newtonian thrust measuring devices,this system can operate in open-loop measurement mode or closed-loop measurement mode,where the measured thrust or pulse bit is related to the displacement signal or the feedback signal respectively.When operating in closed-loop mode,there is an additional advantage of eliminating the errors caused by the thruster’s cable, while greatly reducing the noise in the thrust measurement system.In addition, the application of an electrostatic weak force generator in the thrust measurement system calibrates the thrust response of the system by generating a previously known weak force.Furthermore, the pulsed response of a micro cathode arc thruster was evaluated by the thrust measurement system,which validated the practicality of the system in open-loop mode.In summary, the thrust measurement system was validated for the first time in an elastic torsional pendulum using both open-loop and closed-loop methods for steadystate and impulsive forces, achieving measurement ranges in the sub-micro-Newtons to hundreds of micro-Newtons with the ability of 0.1 micro-Newton magnitude and low background noise.

    In this work, the dual-modes full elastic thrust measurement system proposed shows credibility and practicability under sub-micro-Newton to hundreds of micro-Newton magnitude in section 2, the structure of the thrust measurement system and the weak force generation device will be presented in detail, while in section 3, the steady thrust response and characteristic of the system will be discussed with the dynamic response of a micro cathode arc thruster.Finally, the conclusion will be given in the final part.

    2.Construction and principle of the thrust measurement system

    The image and the schematic diagram of the thrust measurement system are shown in figure 1.The system is composed of the rigid frame,the torsional pendulum of the thrust stand, a displacement sensor and a weak force generation device.The torsional pendulum is connected by two vertically placed shafts and is connected with the rigid frame consequently.The plasma micro thruster is installed at one end of the torsional pendulum, and the counterweight is installed at the other end accordingly.The position limit frame is used to limit the rotation and pitch angle of the axis.The capacitive displacement sensor with the resolution of 40 nm is used to measure the position change of the torsional pendulum, and the micro-Newton weak force is generated by an electrostatic device attached to the torsional pendulum.A magnet sheet is fixed on the other side of the pendulum, and the current generated by the feedback coil interacts with it to pull or push the torsional pendulum in the closed-loop operation.The pipeline and cables of the plasma micro thruster is placed along the pendulum based on the actual demand.

    2.1.Open-loop methodology

    When the thrust measurement system is operating in openloop mode,the pendulum is a regular torsional thrust balance and the principle can be schematically described in figure 2.In figure 2,rexandrsensorare the distance from the axis to the external force point and the displacement sensor respectively, and it is assumed that the external forceFexis perpendicular to the pendulum for convenience.Δxis the displacement signal detected by the displacement sensor,corresponding to the rotation angleθof the pendulum.When the torsional pendulum deviates from the balance position and within a small angle, the restoring moment spontaneously generated by the axis is opposite to the direction of the moment caused by the external force.The theoretical description of the mechanical behavior of the torsional pendulum is well described in the literature [27, 28].Applying Newton’s second law, the result is expressed as [27]:

    whereJis the moment of inertia, λ is the total torsional damping factor,kis the torsional spring constant andFexis the external force mentioned above, which can be generated by either micro thrusters or the weak force generation device.

    When the torsional pendulum deflects within a small angle, the coefficientsJ, λ andkcan be considered as constants [27].As the deflection angle θ is relatively small, the relationship between θ and displacementΔxcan be written as:

    Therefore, the mechanical behavior equation of the torsional pendulum can be written as:

    When it is required to measure the steady thrust, the external torque and the restoring moment should be balanced after the torsional pendulum deflecting to a small angle, and the pendulum will finally stay static or slightly oscillate around the equilibrium position whereclose to zero, consequently the relationship betweenFexandΔxcan be simplified as [27]:

    The impulse of micro thruster can be measured according to the impulse exchange between the thruster and the torsional pendulum.When the impulse generated by thruster is transmitted to the torsional pendulum, the pendulum will oscillate and the amplitude of the oscillation is related to the torsional spring constant and the moment of inertia of the pendulum.The maximum amplitude of oscillationAmaxis defined as[28]:

    HereIbitis the impulse provide by the external force,Jis the moment of inertia andωdis the frequency of damped oscillation defined as follows [28]:

    Here the parameterζ[27] is the damping ratio andandωnis the natural frequency of oscillation of the open-loop system.

    For an oscillation system such as the torsional pendulum,the oscillation amplitude of the pendulum could be described as [28]:

    From the equation (7), the moment of inertia of the torsional pendulum can be acquired by the attenuation envelope of the oscillation curve.

    2.2.Closed-loop methodology

    When the thrust measurement system is operating in the closed-loop mode, an additional torque is applied to the pendulum to keep it in the same position while being affected by the external torque [29, 30].The detailed principle is schematically described in figure 3.

    The feedback closed-loop is composed of the control circuit, feedback electromagnetic coil and permanent magnet sheets.When the feedback loop is activated, the electrical signal received by the displacement sensor is continuously input into the control circuit, and compared with the preset reference voltage in the control circuit.The reference voltage here is the sensor output voltage corresponding to a preset distance between the pendulum and the sensor.When there is a voltage difference between the sensor signal and the reference signal,it means that the current position of the torsional pendulum deviates from the preset position, and the control circuit will output a specific current according to the proportional-integral-derivative (PID) control method [29].The current through the feedback coil will interact with the magnet sheet fixed on the pendulum and generate a certain force to push or pull the pendulum back to the preset position.

    The schematic diagram of the PID control circuit used here is shown in figure 1.The output current is composed of the proportional, integral and derivative operation of the difference voltage between the sensor signal voltage and the preset voltage.When there is an external torque caused by the forceFexacting on the torsional pendulum, the converse feedback torque generated by the PID circuit has to be equal to the external torque, which can counteract the effect of the external torque and make the torsional pendulum return to the preset position.In the experiment, by measuring the force of the feedback circuit, or the feedback voltage signal in the feedback circuit, the external force torque can be indirectly measured.

    The main advantage of using feedback method is that it can effectively avoid the effect of external influence such as mechanical friction,pipeline and cable restraint in the thruster measurement and eliminate the errors caused by external factors as far as possible.In addition,the closed-loop method can effectively suppress the natural oscillation of the torsional pendulum and eliminate the influence effectively, furthermore, it can reduce the thrust measurement noise and effectively improve the accuracy of steady-thrust measurement.

    2.3.The electrostatic weak force generation device

    The electrostatic weak force is generated by two electrostatic comb pairs, and the structure is shown in figure 4.

    The two electrostatic combs are made of aluminum and have 20 fins and 19 fins respectively.The length and thickness of the fins in this project are 10 and 1 mm and the space between two fins is 2 mm.Both combs are insulated from their holders by PTFE insulator, and an electrostatic comb is fixed on the torsional pendulum while another is installed on a two-dimensional displacement table, which can adjust the relative position of the combs precisely.When there is a voltage between the two combs and fins, the electrostatic attraction force will be generated.According to the literature,the approximate equation of electrostatic force between the two combs is described as follows [31, 32]

    HereNis the total number of fins on the combs,candgare the thickness of each fin and the separation between fins respectively, andx0is the engagement distance between fins[32].For the weak force generation device in this experiment,they are constants when the relative position of the combs is determined, and the electrostatic forceFcan be described as:

    Firstly, the weak force generated by the electrostatic comb is measured by a high-precision electronic balance with the resolution of 10μg(Sartorius CPA26P series)in advance.The electrostatic combs are powered by an IT6723G auto range DC stabilized power and it is able to provide voltages from 0 to 600 V with the rise time less than half a second.According to the balance reading and the acceleration of gravityg, the relationship between the force of the electrostatic combs and the voltage applied is obtained.For convenience, the force of the electrostatic comb has been transformed into an equivalent external force at the position of the thrusters under the same torque, and the images of the measurement process and the results of the electrostatic comb are shown in figure 5.

    In the experiment, the voltage range on the electrostatic combs is 0–550 V with the interval of 50 V, and the horizontal axis is the square of the voltage.The fitting line of the electrostatic force and the square of the voltage is obtained by the least-squares method [33] and the equation is:

    Here,a=1.1 ×10?3(μNV?2)andb= ? 6 .1 ×10?3(μN) .The linear goodness of fitR2of the slope value is greater than 0.9999 with the error of slopeδ=1.4 ×10?6,and it is suggested that the electrostatic force has good linearity for calibration.The enlarged part of the blue box in the figure 5 is the straight line fitting curve of the corresponding electrostatic force when the voltage ranges from 0 to 100 V.In this voltage range the linear goodness of fitR2is 0.9997 and the electrostatic force is close to the minimum resolution of the electronic balance which is easy to be disturbed by the environment.According to the results of several measurements of force of the electrostatic comb,the error range of the electrostatic force affected by the comb and external environment at different voltages was noted in figure 5.The range of electrostatic forces generated at higher voltages was larger, when less interference was applied from the ambient airflow and the error ratio was near stable according to the data in figure 5.Analysis of the error ratio showed that the calibration process produces an average error of 4.4%.

    3.Results and discussion

    3.1.Open-loop measurement

    The static calibration process can determine the torsional spring rate of the pendulum system.In the open-loop mode of the experiment, the known constant electrostatic force was applied on the torsional pendulum, and the actual deflection angle of the pendulum was determined according to the output signal of the displacement sensor.Therefore, the rotation angle of the axis or the signal voltage change of the displacement sensor should be proportional to the external force.In the experiment, the voltage applied on the electrostatic comb and the corresponding force were from the calibration chart shown in figure 5.When increasing the voltage on the comb pair step by step,the corresponding displacement could be obtained and recorded by a data recorder with a frequency of 1 kHz.

    Figure 6 exhibits the displacement data from the sensor when the electrostatic comb pair working voltage increased from 0 to 450 V with the increment of 50 V and then returned to 0 V.The thrusts corresponding to each static thrust step were noted in figure 6,and the thrust values were taken from the calibration chart.

    Figure 1.The image and schematic diagram of the thrust measurement system and other experimental setups.

    Figure 2.(a) Torsional pendulum at the rest position and (b)torsional pendulum deflected immediately after an impulse was applied to the arm.

    Figure 3.(a)Torsional pendulum deflected by the external force and(b) back to the equivalent position by the electromagnetic force of feedback coil.

    Figure 4.The image and schematic diagram of the structure of the electrostatic comb.

    Figure 5.The images of the measurement process and the linear fit results of the electrostatic force versus the square of applied voltage with error bars.

    Figure 6.The displacement data versus time for large spans of applied voltage on the electrostatic comb.

    Figure 7.The displacement data versus time for small spans of applied voltage on the electrostatic comb.

    Figure 8.The linear fit of the displacement versus square of voltage.

    Figure 9.The oscillatory decay curve of the torsional pendulum system and the fitted line of the envelope exponential model.

    Figure 10.FFT of natural oscillation of the torsional pendulum.

    As shown in figure 6,the high voltage on the electrostatic comb was maintained for 30 s at each step of external force.When the external forces changed, there were obvious and stable steps of the displacement signal after a few seconds.According to the raw data, when the voltage increased or decreased the torsional pendulum oscillated sinusoidally near the position due to the change in applied force, and then the amplitude gradually decreased because of the damping effect.Because the damping factor of the system was relatively small, the torsional pendulum would oscillate near the equilibrium position in a short period of time instead of being completely stationary.If the torsional pendulum was continuously oscillating, the final equilibrium position was the average value of oscillation according to equation (7).

    In addition, the output signal of the displacement sensor increased gradually when the external force increased.At the moment when the external force changed, large amplitude oscillation sometimes occurred,such as when the electrostatic comb voltage rose from 150 V to 200 V or 250 V to 300 V.On the contrary, there were small oscillations such as the slight oscillation at the point when the voltage of the comb increased from 300 to 350 V.This was probably due to that the torsional pendulum was not completely stationary at the moment when the external force was applied to the torsional pendulum.When the external force which had a same direction to the angular velocity of the pendulum was applied on the torsional pendulum,there would be an enhanced effect which made the torsional pendulum obtain a greater angular velocity, on the contrary, the direction of the external force was opposite to the angular velocity of the torsional pendulum,there would be an opposite effect and it would reduce the angular velocity of the pendulum after deflecting a certain angle, and the amplitude of the oscillation consequently was relatively small.Furthermore, when the voltage on electrostatic comb decreased step by step after 250 s in figure 6, a certain amplitude of oscillation occurred at the moment exactly when the external force decreased.Besides the above factors, this was also probably due to the release of elastic potential energy stored by the spring.When the torsional pendulum returned to the next equilibrium position with the decrease of the external force, the excess elastic potential energy was released and converted into the kinetic energy,and finally the oscillation of the pendulum was produced.

    When applied smaller electrostatic force where the comb voltage ranged from 0 to 100 V and increased with 10 V steps,the relationship between the displacement measured by the sensor and time is shown in figure 7.

    As shown in figure 7,the gray line in the figure is the raw data and the red line is the filtered data by Savitzky–Golay filter [34].When the electrostatic comb voltage increased from 0 to 100 V, the measurement point at the pendulum moved about 18 μm approximately due to the electrostatic force.It was suggested by the filtered data that the torsional pendulum could maintain stable at the equilibrium position with the electrostatic force increased and decreased.

    In order to evaluate the linearity of the deflection displacement, a linear fit of the displacement versus square of voltage from 0 to 450 V was performed and the results are shown in figure 8.

    As shown in figure 8, the red line was the fitting line of the torsional pendulum displacement and square of applied voltage when the voltage of the electrostatic comb ranged from 0 to 400 V.The fitting equation is:

    The linear goodness of fitR2of the slope value was greater than 0.9999 with the standard error of slopeδ=4.6 ×10?5,and it was suggested that the displacement had good linearity for calibration.According to equations (10) and (11), the relationship between the displacement and the applied force was obtained by eliminating the square term of voltage:

    Equation (12) can be considered as a proportional relationshipF=(k1/k2)Δxand the deviation of the final fitted force wasδF/F=δk1/k1Δx+k1/k2δΔx?δk2/k2Δx.From the relationship between the fitted force and the displacement of the pendulum, and the error was found to be approximately 1.2%.Considering the effect caused by the procedure of the calibration of the electrostatic comb,and the error of open-loop measurement was approximately 5.5%.In addition, according to the straight line fitting model of relationship between deflection displacement and the external force, and the effective force arm of 0.5 m, the torsional spring rate of the pendulum system could be obtained as:

    Similarly, when the electric comb voltage ranged from 0 to 100 V the measured displacement voltage straight line of the torsional pendulum fits the same straight line and the fitting slope was1.7 ×10?3(μm V?2)with a linear goodness of fit of 0.9999.It was suggested that torsional pendulum had a good linear ability with small force of thrusters.

    It is needed to acquire the dynamic characteristics and parameters of the torsional pendulum system when measuring the impulsive external force.According to equation (7),the oscillation of the torsional pendulum system satisfied the characteristic of exponential attenuation [28], in which the attenuation coefficient was related to the damping rate and the moment of inertia of the pendulum.Therefore, the relationship between the damping coefficient and the moment of inertia could be obtained by recording the change of the voltage signal corresponding to the displacement with time during the decay of the free oscillation of the torsion pendulum.The oscillation attenuation curve of the torsional pendulum system and the fitting line of an exponential model of the envelope line are shown in figure 9.

    As shown in figure 9, the horizontal axis was time and the vertical axis was the absolute displacement of the displacement sensor, and the black line was the displacement value of the corresponding point of the torsional pendulum,and the amplitude decreased gradually to the center with time.The blue line was the upper envelope of the displacement oscillation, the red line was the exponential model fitting of

    the envelope with the fitting equationy=A1·exp ( ?x/t1)+y0.The results showed that the goodness of fit was over 0.9999, indicating that thecurve fitted by the index was relatively accurate.According to the fitting equation,the termt1=119.63 ±0.31,and the damping termcan be obtained according to equation (7) asconsequently the damping ratio could be obtained as:

    Figure 10 shows the FFT used to determine the natural frequency of the torsional pendulum oscillation.

    Figure 11.The feedback signal versus time for large spans of applied voltage on the electrostatic comb.

    As shown in figure 10, the frequency of the damping oscillation of the torsional pendulum was 0.41 Hz, and the full width at half maximum (FWHM) was 0.014 Hz.According to equations(6)and(14),the natural frequency of the system could be obtained asfn=0.41 Hz.The natural frequency of the pendulum system was obtained from the FFT of the oscillation, and the natural frequency of the oscillation was used to determine the moment of inertia according to the equation [28, 35]:

    Consequently the moment of inertiaJcould be obtained according to equation (15) as 0.048 kg m .2

    3.2.Closed-loop measurement

    The closed-loop control method can be used to measure the static force, and its advantage is that it can reduce the influence of external factors such as cables.In the closed-loop mode,the feedback force which is balanced with the torque of external force is applied to the torsional pendulum to maintain the pendulum at a specific position(as shown in figure 3).At the moment when the external force changes, the feedbackcontrol coil applies a balanced force through the PID control circuit,and finally realizes the balance with the external force again.In this experiment, the input signal of the PID control loop was the position signal of the displacement sensor, and the output signal was the current on the feedback coil in figure 1.Therefore,the corresponding variable of the external force measurement was transformed from the position signal to the voltage signal on the feedback device.

    Figure 11 exhibits the feedback voltage of the coil when the voltage on the electrostatic comb pair increased from 0 to 450 V with the increment of 50 V and then returned to 0 V.The corresponding thrust of each static thrust step was noted and the enlarged image showed the raw data of electrostatic comb voltage from 0 to 50 V.

    Figure 12.The feedback signal versus time for small spans of applied voltage on the electrostatic comb.

    Figure 13.The linear fit of the feedback signal versus square of voltage from 0 to 400 V.

    Figure 14.Background noise spectrum density of the system in open-loop and closed-loop mode.

    Figure 15.The μCAT ignition test image(left)and the displacement signal of pendulum versus time (right).

    As shown in figure 11, the electrostatic comb was maintained at the each voltage step for 30 s,corresponding to the electrostatic force range from 0 to 171.1μN previously measured by the electronic balance in figure 4.From the raw data in figure 11, it could be suggested that the feedback voltage signal changed accordingly when the external force increased or decreased, and the electromagnetic force generated by the feedback coil made an equilibrium with the external force torque consequently and the corresponding feedback voltage recorded was about 1.5–8.5 V.In addition,it was suggested from the data line that when the external force changed, there were almost no sinusoidal oscillations produced by the feedback signal as similar as those shown in figure 6.On the contrary, the feedback signal reached the balance position relatively quickly.This was because the PID control circuit has the ability of rapid response, and its differential control branch (as shown in figure 1) could make a certain prediction of the input signal, so as to control the torsional pendulum under the combined effect of PID circuit and made it return to the equilibrium position as soon as possible.

    In figure 11,the enlarged image is the raw data line with time around the moment when the voltage of the electrostatic comb increased from 0 to 50 V, which produced an obvious step of the feedback voltage signal at the time when the force was applied.Similarly, when the electrostatic comb voltage ranged from 0 to 100 V and increased at 10 V intervals, the feedback voltage was measured under several relatively smaller external forces.The relationship between the feedback signal and time is shown in figure 12.

    As shown in figure 12, the gray line is the raw data and the red line is the filtered data by Savitzky–Golay filter.When the electrostatic comb voltage increased from 0 to 100 V,the feedback voltage signal ranged from 1.1 to 1.5 V approximately.The enlarged image is the corresponding data line when the voltage of electrostatic comb increased from 0 to 30 V.Although the raw data corresponding to the gray line had some signal noise, the filtered signal line still clearly distinguished the feedback voltage corresponding to the external force of submicron Newton magnitude.According to the previous calibration results of external force, when the voltage of electrostatic comb increased from 0 to 10 V, the corresponding increase of equivalent external force was 0.1μN, and the data curve after filtering showed that the voltage signal increased distinguishably, and the amplitude was approximately 5 mV.

    In order to evaluate the linearity of deflection displacement,the results of a linear fit of the displacement versus the square of the voltage are shown in figure 13.

    As shown in figure 13,the red line is the fitting line of the feedback voltage signal and square of applied voltage when the voltage of the electrostatic comb ranged from 0 to 400 V.The fitting equation is:

    The linear goodness of fitR2of the slope value was greater than 0.9999 with the standard error of the slopeδ=1.8 ×10?8,and it was suggested that the feedback voltage had good linearity with external force compared to that of the open-loop mode.According to equations (10) and(15), the relationship between the displacement and the applied force was obtained by eliminating the square term of voltage:

    Similarly, when the electric comb voltage was 0–100 V and the interval was 10 V, the measured feedback signal straight line of the torsional pendulum fit the straight line and the fitting relationship had a linear goodness of fit over than 0.9999.It was suggested that torsional pendulum had a good linear ability with the small force of thrusters at the closedloop mode.From the relationship between the fitted force and the feedback voltage,the error was found to be approximately 0.9%in the closed-loop.Considering the effect caused by the procedure of the calibration of the electrostatic comb,and the error of open-loop measurement was approximately 5.3%.

    3.3.Background noise analysis

    Plasma micro thrusters of sub-micro-Newton to hundreds of micro-Newton thrust usually exhibit low thrust noise,thus the background noise of the thrust measurement system should be low enough to better reflect the parameters and characteristics of the micro thrusters.In this experiment,the output signal of the thrust stand was measured with no external disturbance in the open-loop mode and closed-loop mode respectively after a period of resting until stabilization,and the signal was converted to the equivalent external force perpendicular to the torsional pendulum (Fexat the thruster mounting position in figure 2) by equations (12) and (16), respectively.The background noise spectrum of the system was obtained by analyzing the power spectral density of the signal, as shown in figure 14.

    As shown in figure 14, the gray curve is the background noise power spectral density of the system measured in openloop mode,and the red curve is the noise power spectral density in closed-loop mode.In the open-loop mode, the noise power spectral density had a peak near 0.41 Hz with an amplitude of about 0.4which was caused by the inherent oscillation frequency of the torsional system itself,and if the damping method, such as electromagnetic vortex damping, was applied,the background noise generated by the inherent frequency oscillation would be greatly reduced.In addition, this inherent frequency was known and digital filtering methods could be used to reduce the noise here, minimizing interference with the thrust measurements of the plasma micro thrusters.In other frequency bands with frequencies greater than 0.01 Hz, the frequency spectral density of the background noise in open-loop mode was less than 0.1in general, with a roughly trend of decreasing with increasing frequency.

    The power spectral density curve shown in the red curve in figure 14 shows that the noise power spectral density in the closed-loop mode generally tended to decrease with the increase of frequency.When the frequency was less than 0.01 Hz,the background noise in closed-loop mode and openloop mode was approximately the same or in one order of magnitude, while with the gradual increase of frequency, the background noise in closed-loop mode decreased significantly, and the amplitude of the noise spectrum in the interval greater than 0.01 Hz was basically less than 0.1and the amplitude in the frequency zone greater than 0.1 Hz was greatly reduced compared with that in the open-loop mode, and the difference was generally about one order of magnitude.It is worth noting that in the closedloop mode there was no peak near the intrinsic frequency of the system as described above, which indicated that the feedback loop could control the motion of the torsional pendulum and keep it as stable as possible in the equilibrium position well.

    3.4.Dynamic evaluation of a μCAT

    Micro cathode arc thruster (μCAT) is a kind of pulsed miniature plasma thruster based on vacuum arc ablation [9].The main principle is that the vacuum arc is formed between the cathode and anode under the applied electric field and the conductive film, and the metal vapor from the arc ablation is accelerated and ejected out of the thruster to generate thrust under the coupling effect of heat energy and electromagnetic field.Usually μCAT operates in the pulsed mode, and the pulse bit of the thruster is approximately in the magnitude of 0.01–10μN s depending on the size and discharge energy of the thruster,cathode and anode material,discharge parameters and magnetic field characteristics.

    In this experiment, the dynamic response of the thrust measurement system operated on the pulsed discharge mode of the μCAT with titanium anode was obtained, and the μCAT ignition test image taken with long exposure and the displacement signal output from the sensor versus time is shown in figure 15.

    As shown in figure 15, the μCAT was operating in pulsed mode injecting plasma at the exit,and the process was recorded where the gray curve represented for the raw data and the red curve represented for Savitzky–Golay filter processed data.The μCAT ignited once at the moment of 4.5 s as shown in figure 15, at that time a spike of signal interference was generated due to the transient high current of the thrusters,and then the torsional pendulum started to oscillate.The maximum amplitudeA2maxof 0.56 μm was obtained from the processed results, consequently the impulse value was calculated from equation (5) as 0.07μN s.The error in impulse bit was found to be 21.6% according to the raw data due to the limit of the data recorder, and 3.4% with respect to theJand other parameters.Because of the weak magnetic field of the μCAT in the experiment, the impulse of the μCAT was relatively small and only the magnitude can be evaluated.In the follow-up work the measurement system using more precise displacement sensor will be proposed to make it more suitable for sub-micro-Newton plasma thrusters.

    4.Summary and conclusion

    In this work, a device for measuring the thrust of plasma micro thrusters with sub-micro-Newton to hundreds micro-Newtons was proposed, and the structure of the device and the measurement principles in open-loop and closed-loop modes were introduced respectively, and a capacitive sensor was used to measure the displacement to reflect the oscillatory motion of the torsional pendulum and the feedback voltage was obtained to indirectly measure the equilibrium force.An electrostatic weak force generation device was used to simulate the external force generated at the thruster mounting position, and the linear fitting relationship between the external force and the output signal of the thrust measurement system could be obtained by the measurement results of the electronic balance previously, and it was suggested that the thrust measurement system had good fitting linearity in the range of 0–200μN approximately, with the ability to resolve external force variations of 0.1μN.Based on the relationship between the angular deflection of the torsional pendulum and the applied force, the torsional spring rate of the system was obtained, and the damping and inertia parameters were obtained consequently from the pendulum decay oscillation.In addition,the analysis results of the background noise suggested that the thrust measurement system was able to achieve a noise level of 0.1/μN Hz or less in open-loop and closed-loop modes in the frequency range from 10 mHz to 5 Hz approximately, which was advantageous for plasma micro thruster parameters and characteristics measurements.The actual micro cathode arc thruster impulse discharge test results showed that the thruster could produce impulses of 0.1μN s, which could be distinguished by the thrust measurement system.

    Acknowledgments

    This work was partly supported by the Shanghai Engineering Research Center of Space Engine (No.17DZ2280800).

    猜你喜歡
    海濱
    Fringe visibility and correlation in Mach–Zehnder interferometer with an asymmetric beam splitter
    夏日海濱
    Structural Reliability Analysis Based on Support Vector Machine and Dual Neural Network Direct Integration Method
    海濱的夏天
    孩子(2019年9期)2019-11-07 01:35:49
    古城里的海濱新居
    海濱書簡
    散文詩(2017年17期)2018-01-31 02:34:19
    海濱1
    海濱
    講一題通一類得一法
    海濱風(fēng)光掠影
    麻豆av噜噜一区二区三区| 亚洲成人久久性| 美女大奶头视频| 国产黄a三级三级三级人| 午夜老司机福利剧场| 性欧美人与动物交配| 成人一区二区视频在线观看| 在线a可以看的网站| 国产男靠女视频免费网站| 欧美黄色片欧美黄色片| 毛片女人毛片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩欧美国产一区二区入口| 丰满乱子伦码专区| 国产麻豆成人av免费视频| 淫妇啪啪啪对白视频| 亚洲经典国产精华液单 | 久久伊人香网站| 少妇高潮的动态图| 亚洲综合色惰| 日本一二三区视频观看| 黄色丝袜av网址大全| 欧美激情久久久久久爽电影| 久久精品夜夜夜夜夜久久蜜豆| 我要搜黄色片| 亚洲第一欧美日韩一区二区三区| 人妻制服诱惑在线中文字幕| 日韩中字成人| 在线观看美女被高潮喷水网站 | 亚洲国产高清在线一区二区三| 草草在线视频免费看| av视频在线观看入口| 国产69精品久久久久777片| 亚洲性夜色夜夜综合| 可以在线观看的亚洲视频| 高潮久久久久久久久久久不卡| 精品不卡国产一区二区三区| 久久精品国产清高在天天线| 亚洲成av人片免费观看| 日本一本二区三区精品| 悠悠久久av| 免费av毛片视频| 日韩av在线大香蕉| 级片在线观看| 亚洲中文字幕一区二区三区有码在线看| 色吧在线观看| 欧美zozozo另类| 成年女人毛片免费观看观看9| 9191精品国产免费久久| 少妇人妻一区二区三区视频| 亚洲aⅴ乱码一区二区在线播放| 性欧美人与动物交配| 免费黄网站久久成人精品 | 桃红色精品国产亚洲av| 国产aⅴ精品一区二区三区波| av福利片在线观看| 日本精品一区二区三区蜜桃| 熟女人妻精品中文字幕| av在线老鸭窝| 99国产精品一区二区三区| 中文字幕av成人在线电影| 亚洲乱码一区二区免费版| a级毛片a级免费在线| 国产人妻一区二区三区在| 午夜精品一区二区三区免费看| 精品国产亚洲在线| 久久精品影院6| 久久久色成人| 国产精品电影一区二区三区| 免费在线观看影片大全网站| 淫秽高清视频在线观看| 国产v大片淫在线免费观看| 国产精品国产高清国产av| 3wmmmm亚洲av在线观看| 亚洲黑人精品在线| 一本精品99久久精品77| 亚洲成人久久性| 国产一区二区激情短视频| 欧美日韩福利视频一区二区| 亚洲无线观看免费| 18禁黄网站禁片免费观看直播| 舔av片在线| 日本撒尿小便嘘嘘汇集6| 国产精品人妻久久久久久| 波多野结衣高清无吗| a级毛片免费高清观看在线播放| 亚洲成人久久爱视频| 国产伦精品一区二区三区四那| ponron亚洲| 国产精品日韩av在线免费观看| bbb黄色大片| 成人精品一区二区免费| 日韩成人在线观看一区二区三区| 中国美女看黄片| 国产在线精品亚洲第一网站| 好男人电影高清在线观看| 舔av片在线| 成年女人毛片免费观看观看9| 成人国产综合亚洲| 久久中文看片网| 国产亚洲精品久久久com| 久久久精品大字幕| 日韩人妻高清精品专区| 性色av乱码一区二区三区2| 美女高潮的动态| 18+在线观看网站| 亚洲专区国产一区二区| 国产视频一区二区在线看| 身体一侧抽搐| 永久网站在线| 国产视频一区二区在线看| 国内精品一区二区在线观看| 精品人妻偷拍中文字幕| 亚洲精品在线美女| 日韩欧美在线乱码| 99久久精品国产亚洲精品| 午夜精品在线福利| 精品免费久久久久久久清纯| 嫩草影院新地址| 超碰av人人做人人爽久久| 少妇裸体淫交视频免费看高清| 一a级毛片在线观看| 村上凉子中文字幕在线| 3wmmmm亚洲av在线观看| 久久欧美精品欧美久久欧美| 国产精品一区二区性色av| av福利片在线观看| 欧美日韩国产亚洲二区| 国产精品永久免费网站| 国产成人a区在线观看| 免费人成在线观看视频色| 亚洲,欧美,日韩| 黄色一级大片看看| 免费在线观看影片大全网站| 日韩中字成人| 国产三级黄色录像| 一夜夜www| 99在线视频只有这里精品首页| 成熟少妇高潮喷水视频| 日日摸夜夜添夜夜添小说| 欧美不卡视频在线免费观看| 中文字幕av成人在线电影| 高清日韩中文字幕在线| 我要搜黄色片| 国产麻豆成人av免费视频| 日韩欧美国产在线观看| 国产精品一区二区三区四区免费观看 | 天美传媒精品一区二区| 午夜精品一区二区三区免费看| 亚洲精品亚洲一区二区| 男女视频在线观看网站免费| 99久久精品一区二区三区| 亚洲精品一区av在线观看| 免费搜索国产男女视频| 国产亚洲欧美98| 老熟妇乱子伦视频在线观看| 综合色av麻豆| 国产av一区在线观看免费| 国产午夜福利久久久久久| 国产一区二区三区在线臀色熟女| 99久久成人亚洲精品观看| 国产精品一区二区性色av| 极品教师在线免费播放| 精品人妻一区二区三区麻豆 | av专区在线播放| av在线观看视频网站免费| 久久6这里有精品| 国产综合懂色| 亚洲自拍偷在线| 日韩欧美精品免费久久 | 狠狠狠狠99中文字幕| 又爽又黄a免费视频| 日韩人妻高清精品专区| 在线观看午夜福利视频| netflix在线观看网站| 色噜噜av男人的天堂激情| 高清毛片免费观看视频网站| 色5月婷婷丁香| www.色视频.com| 国产真实乱freesex| 国产成+人综合+亚洲专区| 久久这里只有精品中国| 非洲黑人性xxxx精品又粗又长| 亚洲欧美日韩东京热| 国产91精品成人一区二区三区| 可以在线观看的亚洲视频| 亚洲国产精品合色在线| 国产精品三级大全| 中文在线观看免费www的网站| 欧美一区二区精品小视频在线| 亚洲av一区综合| 99久久成人亚洲精品观看| 国产乱人视频| 啪啪无遮挡十八禁网站| 亚洲久久久久久中文字幕| 精品福利观看| 99国产极品粉嫩在线观看| 精品久久久久久久久亚洲 | 亚洲国产高清在线一区二区三| 18+在线观看网站| 亚洲无线观看免费| 精品国内亚洲2022精品成人| 黄色丝袜av网址大全| 日日摸夜夜添夜夜添小说| 变态另类成人亚洲欧美熟女| 日韩成人在线观看一区二区三区| 亚洲av电影在线进入| 真人一进一出gif抽搐免费| 国产av一区在线观看免费| 亚洲经典国产精华液单 | 91麻豆av在线| 欧美乱色亚洲激情| 最好的美女福利视频网| 女生性感内裤真人,穿戴方法视频| 亚洲欧美日韩高清专用| 午夜久久久久精精品| 国产精品98久久久久久宅男小说| 热99在线观看视频| 性插视频无遮挡在线免费观看| 日本三级黄在线观看| 在线观看66精品国产| 色综合站精品国产| 国产大屁股一区二区在线视频| 久久久久亚洲av毛片大全| 在线观看午夜福利视频| 日韩国内少妇激情av| 亚洲av美国av| 高清在线国产一区| 亚洲国产精品sss在线观看| 久久性视频一级片| 免费无遮挡裸体视频| 色哟哟·www| 亚洲avbb在线观看| 嫩草影视91久久| 欧美成人a在线观看| 看黄色毛片网站| www.色视频.com| 少妇人妻精品综合一区二区 | 尤物成人国产欧美一区二区三区| 欧美成人免费av一区二区三区| 97超级碰碰碰精品色视频在线观看| 亚洲人与动物交配视频| av女优亚洲男人天堂| 精品人妻熟女av久视频| 日本撒尿小便嘘嘘汇集6| 97超级碰碰碰精品色视频在线观看| 成人一区二区视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 色吧在线观看| 蜜桃久久精品国产亚洲av| 久久精品人妻少妇| 一本一本综合久久| 午夜两性在线视频| 精品久久久久久久末码| 伊人久久精品亚洲午夜| 露出奶头的视频| 国产大屁股一区二区在线视频| 少妇人妻一区二区三区视频| 日本撒尿小便嘘嘘汇集6| 精品一区二区三区视频在线观看免费| 久久亚洲真实| 亚洲va日本ⅴa欧美va伊人久久| АⅤ资源中文在线天堂| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 99精品久久久久人妻精品| 亚洲av二区三区四区| 综合色av麻豆| 国产伦精品一区二区三区四那| 欧美午夜高清在线| 色尼玛亚洲综合影院| 人人妻,人人澡人人爽秒播| 嫁个100分男人电影在线观看| 欧美乱色亚洲激情| 欧美性猛交╳xxx乱大交人| 高清在线国产一区| 久久热精品热| 亚洲熟妇熟女久久| 最近在线观看免费完整版| 午夜激情福利司机影院| 99精品久久久久人妻精品| 久99久视频精品免费| 亚洲人与动物交配视频| 国产人妻一区二区三区在| 国产亚洲精品久久久com| 精品人妻一区二区三区麻豆 | 成人av在线播放网站| 激情在线观看视频在线高清| av欧美777| 国产日本99.免费观看| 三级国产精品欧美在线观看| 久久亚洲真实| 成年版毛片免费区| 日日干狠狠操夜夜爽| 亚洲午夜理论影院| 亚洲最大成人av| 成人美女网站在线观看视频| 久久久久久久久久黄片| av在线老鸭窝| 欧美黄色片欧美黄色片| 最近最新中文字幕大全电影3| 精华霜和精华液先用哪个| 无人区码免费观看不卡| 成人毛片a级毛片在线播放| 日韩欧美精品免费久久 | 久久久国产成人免费| 99久久九九国产精品国产免费| 国产精品爽爽va在线观看网站| 免费观看精品视频网站| 丰满的人妻完整版| 欧洲精品卡2卡3卡4卡5卡区| 久久亚洲精品不卡| 色av中文字幕| 婷婷色综合大香蕉| 久久草成人影院| 国产精品影院久久| 久久久久性生活片| 97碰自拍视频| 热99re8久久精品国产| 12—13女人毛片做爰片一| 日本一二三区视频观看| 亚洲电影在线观看av| 欧美黑人巨大hd| 黄色配什么色好看| 日韩av在线大香蕉| 国语自产精品视频在线第100页| 午夜亚洲福利在线播放| 欧美+日韩+精品| a级毛片免费高清观看在线播放| 啦啦啦观看免费观看视频高清| 丰满的人妻完整版| 国内久久婷婷六月综合欲色啪| 91av网一区二区| 美女高潮的动态| 午夜免费成人在线视频| 国产亚洲欧美98| 亚洲精品在线观看二区| 久久99热6这里只有精品| 午夜福利高清视频| 国产精品美女特级片免费视频播放器| 99久久久亚洲精品蜜臀av| 欧美乱妇无乱码| 在线天堂最新版资源| 老鸭窝网址在线观看| 在线十欧美十亚洲十日本专区| 极品教师在线免费播放| 91狼人影院| 亚洲性夜色夜夜综合| 国产亚洲精品综合一区在线观看| 极品教师在线免费播放| 一个人看视频在线观看www免费| 亚洲av电影在线进入| 成人精品一区二区免费| 国产精品精品国产色婷婷| 国产亚洲精品综合一区在线观看| 欧美性猛交╳xxx乱大交人| 18+在线观看网站| 九九热线精品视视频播放| a级一级毛片免费在线观看| 精品久久久久久久久av| 看黄色毛片网站| 亚洲国产欧美人成| 波多野结衣巨乳人妻| 亚洲熟妇熟女久久| 中文字幕高清在线视频| 噜噜噜噜噜久久久久久91| 在线国产一区二区在线| 最近视频中文字幕2019在线8| 99久久99久久久精品蜜桃| 琪琪午夜伦伦电影理论片6080| 国产伦精品一区二区三区四那| 床上黄色一级片| 人妻丰满熟妇av一区二区三区| 欧美精品国产亚洲| 窝窝影院91人妻| 亚洲国产欧洲综合997久久,| 免费观看精品视频网站| 九九热线精品视视频播放| 亚洲中文日韩欧美视频| 超碰av人人做人人爽久久| 国产免费男女视频| 啪啪无遮挡十八禁网站| 亚洲人成网站在线播放欧美日韩| 午夜福利视频1000在线观看| 一区二区三区高清视频在线| 欧美xxxx性猛交bbbb| 午夜福利高清视频| 无人区码免费观看不卡| 亚洲自拍偷在线| 国产精品久久久久久久久免 | 国产视频内射| 婷婷色综合大香蕉| 高清在线国产一区| 国产精品一区二区三区四区久久| 亚洲在线自拍视频| 成人特级av手机在线观看| 日韩欧美在线乱码| 日本a在线网址| 欧美激情久久久久久爽电影| 久久午夜福利片| 婷婷精品国产亚洲av| 老鸭窝网址在线观看| 一级黄色大片毛片| 免费观看的影片在线观看| 日本黄大片高清| 日韩亚洲欧美综合| 久久久久亚洲av毛片大全| 色视频www国产| 无遮挡黄片免费观看| 国产成人欧美在线观看| 草草在线视频免费看| 日韩欧美国产在线观看| 99在线人妻在线中文字幕| a级一级毛片免费在线观看| 亚洲最大成人av| 久久久久久久久久成人| 欧美激情在线99| 欧美最新免费一区二区三区 | 成年版毛片免费区| 美女cb高潮喷水在线观看| 动漫黄色视频在线观看| 三级毛片av免费| 国产午夜精品久久久久久一区二区三区 | av天堂在线播放| 国产视频内射| 国内揄拍国产精品人妻在线| 国产精品一区二区免费欧美| 有码 亚洲区| 中文字幕高清在线视频| 日日干狠狠操夜夜爽| .国产精品久久| 午夜福利视频1000在线观看| 91字幕亚洲| 国产美女午夜福利| 午夜福利在线观看免费完整高清在 | 国产乱人伦免费视频| 精品久久久久久,| 亚洲av一区综合| 亚洲avbb在线观看| 白带黄色成豆腐渣| 日韩欧美精品v在线| 亚洲aⅴ乱码一区二区在线播放| 欧美3d第一页| 亚洲国产高清在线一区二区三| 亚洲国产欧洲综合997久久,| 国产精品永久免费网站| 午夜福利在线在线| 欧美xxxx黑人xx丫x性爽| 成人永久免费在线观看视频| 国产精品日韩av在线免费观看| 国产成人a区在线观看| 午夜精品久久久久久毛片777| 网址你懂的国产日韩在线| 两个人视频免费观看高清| 亚洲专区国产一区二区| 简卡轻食公司| 亚洲欧美日韩无卡精品| 中文字幕人妻熟人妻熟丝袜美| 免费高清视频大片| 欧美日本视频| 丝袜美腿在线中文| 永久网站在线| 最新在线观看一区二区三区| 很黄的视频免费| 精品久久久久久久久亚洲 | 国产亚洲精品久久久久久毛片| 成年人黄色毛片网站| 久久热精品热| 极品教师在线视频| 久久久久久久久久黄片| netflix在线观看网站| 69av精品久久久久久| 女同久久另类99精品国产91| 无人区码免费观看不卡| 久久精品久久久久久噜噜老黄 | 国产亚洲精品久久久久久毛片| 日韩亚洲欧美综合| 波多野结衣高清无吗| 欧美日韩亚洲国产一区二区在线观看| 国产精品一及| a级一级毛片免费在线观看| 观看美女的网站| 婷婷精品国产亚洲av| 好男人在线观看高清免费视频| 免费看a级黄色片| 中出人妻视频一区二区| 中文字幕人成人乱码亚洲影| 色av中文字幕| 国产一区二区在线观看日韩| 一级黄片播放器| 日本黄大片高清| 午夜精品久久久久久毛片777| 日韩国内少妇激情av| 亚洲狠狠婷婷综合久久图片| 精品99又大又爽又粗少妇毛片 | 人人妻,人人澡人人爽秒播| 高清毛片免费观看视频网站| 欧美日韩瑟瑟在线播放| 亚洲国产欧美人成| 欧美国产日韩亚洲一区| 男人舔奶头视频| 国产精品久久视频播放| 国产精品国产高清国产av| 性色av乱码一区二区三区2| 国产亚洲精品av在线| 天堂√8在线中文| 国产精品一区二区免费欧美| 色吧在线观看| 国产综合懂色| 美女黄网站色视频| 久久国产精品人妻蜜桃| 黄片小视频在线播放| 国产精品一区二区免费欧美| 亚洲av免费在线观看| 久久国产精品人妻蜜桃| 日韩欧美三级三区| 淫妇啪啪啪对白视频| 一二三四社区在线视频社区8| 男女那种视频在线观看| 我的女老师完整版在线观看| 最近最新中文字幕大全电影3| 亚洲精品色激情综合| 国产高清三级在线| 在线免费观看的www视频| 精品久久久久久久久久免费视频| 亚洲国产欧洲综合997久久,| a级一级毛片免费在线观看| 丝袜美腿在线中文| 在线观看av片永久免费下载| 亚洲不卡免费看| 丰满乱子伦码专区| 久久国产精品影院| 免费在线观看亚洲国产| 国产私拍福利视频在线观看| netflix在线观看网站| 日韩欧美精品免费久久 | 色尼玛亚洲综合影院| 亚洲无线在线观看| 色5月婷婷丁香| av视频在线观看入口| 老司机深夜福利视频在线观看| 嫩草影院精品99| 亚洲一区二区三区不卡视频| 九九热线精品视视频播放| 成人欧美大片| 长腿黑丝高跟| 久久国产精品人妻蜜桃| 搡女人真爽免费视频火全软件 | 婷婷亚洲欧美| 极品教师在线免费播放| av黄色大香蕉| 3wmmmm亚洲av在线观看| 亚洲av电影不卡..在线观看| 在线观看一区二区三区| 国产v大片淫在线免费观看| 日韩人妻高清精品专区| 老司机午夜十八禁免费视频| 国产男靠女视频免费网站| 国产精品伦人一区二区| 五月伊人婷婷丁香| 欧美国产日韩亚洲一区| 99热只有精品国产| 免费观看精品视频网站| 国产一区二区激情短视频| 欧美日本视频| 亚洲欧美日韩卡通动漫| 三级男女做爰猛烈吃奶摸视频| 久久性视频一级片| 日本免费一区二区三区高清不卡| 免费在线观看亚洲国产| 一级a爱片免费观看的视频| 欧洲精品卡2卡3卡4卡5卡区| 黄片小视频在线播放| 精品国内亚洲2022精品成人| 90打野战视频偷拍视频| 两性午夜刺激爽爽歪歪视频在线观看| 两个人视频免费观看高清| 成年版毛片免费区| 欧美乱妇无乱码| 波多野结衣高清无吗| 中出人妻视频一区二区| 97超视频在线观看视频| 人人妻人人看人人澡| 两个人视频免费观看高清| 欧美三级亚洲精品| 亚洲av电影不卡..在线观看| 综合色av麻豆| 精品一区二区三区视频在线观看免费| 老司机福利观看| 日韩欧美三级三区| 色综合婷婷激情| а√天堂www在线а√下载| 欧美不卡视频在线免费观看| 欧美成人一区二区免费高清观看| 欧美三级亚洲精品| 麻豆成人午夜福利视频| 亚洲精品亚洲一区二区| xxxwww97欧美| 日韩中字成人| 精品人妻一区二区三区麻豆 | 亚洲五月婷婷丁香| 人妻夜夜爽99麻豆av| 给我免费播放毛片高清在线观看| 国产午夜福利久久久久久| 日韩有码中文字幕| 国产一区二区三区视频了| 欧美成人a在线观看| 午夜福利在线观看免费完整高清在 | 亚洲黑人精品在线| 在线观看66精品国产| 国产欧美日韩一区二区精品| 国产精品影院久久| 国产精品久久电影中文字幕| 精品午夜福利在线看| 在线播放无遮挡| 天天躁日日操中文字幕|