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

    2D Finite Element Analysis of Asynchronous Machine Influenced Under Power Quality Perturbations

    2022-03-14 09:26:26JasminPamelaSaranyaIndragandhiRajaSinghSubramaniyaswamyYuvarajaTeekaramanShabanaUroojandNorahAlwadai
    Computers Materials&Continua 2022年3期

    Jasmin Pamela S,R.Saranya,V.Indragandhi,R.Raja Singh,V.Subramaniyaswamy,Yuvaraja Teekaraman,Shabana Uroojand Norah Alwadai

    1School of Electrical Engineering,Vellore Institute of Technology,Vellore,India

    2School of Computing,SASTRA Deemed University,Thanjavur,India

    3Mobility,Logistics,and Automotive Technology Research Centre,Faculty of Engineering,ETEC,Department of Electrical Engineering and Energy Technology,Vrije Universiteit Brussel,Brussel,1050,Belgium

    4Department of Electrical Engineering,College of Engineering,Princess Nourah bint Abdulrahman University,Riyadh,84428,Saudi Arabia

    5Department of Physics,College of Sciences,Princess Nourah bint Abdulrahman University,Riyadh,11671,Saudi Arabia

    Abstract: Asynchronous machines are predominantly preferred in industrial sectors for its reliability.Power quality perturbations have a greater impact on industries;among the different power quality events,voltage fluctuations are the most common and that may cause adverse effect on machine’s operation since they are longer enduring.The article discusses a numerical technique for evaluating asynchronous motors while taking into account magnetic saturation, losses, leakage flux, and voltage drop.A 2D linear analysis involving a multi-slice time stepping finite element model is used to predict the end effects.As an outcome,the magnetic saturation and losses are estimated using a modified 2D nonlinear time-stepping finite element formulation.The method takes the electromagnetic fields at the ends of the motor into account using limited computer resources.The proposed method will greatly reduce computation time with limited computer resources for analyzing the machine’s performance with high precision.The analyzed findings assist in preventing voltagevariance issues in the power network system and provide suggestions for developing a robust system.

    Keywords: Asynchronous machine; finite element model analysis; voltage variations; power quality; multi-slice model

    1 Introduction

    It is well-known fact that more than 50 percent of global electrical energy is utilized by electrical machines, especially by asynchronous machines.The induction machine is known as asynchronous machine, actually it is a more precise designation.Asynchronous machines are engaged with various industrial applications including, blowers, air conditioning systems, pumps,refrigeration, materials processing and manufacturing sectors.Three-phase asynchronous machines(AM) are widely preferred as they are reliable, robust, self-starting, and require less maintenance.Nevertheless, the power quality problems in grid are inevitable, even in advanced power networks.Hence, most of the industries are influenced by power quality perturbations.IEEE 1159-2019 defined the power quality problems as transients, sag, swell, interruptions, voltage fluctuations,waveform distortion, voltage variations, voltage imbalance and power frequency variations [1].Among the power quality events, the voltage variations are most frequent, which may cause adverse effect on machine’s operation.Although AMs are highly efficient at full load operating conditions, they are highly sensitive to supply voltage variations.When the supply voltage is increased, the magnetic components of the AM are pushed into saturation.These variations in supply voltage could slow down the operation, degrade the machine, even burnout the machine leading to interruption in the process industries.Hence it is vital to investigate the voltage variations on AM [2].As per the IEEE 1159-2019 Std., Overvoltage and undervoltage are described as a fluctuate in the RMS value of voltage by more than 110 V or less than 90 V for more than one minute at the power frequency.The voltage changes into two categories based on their form and period: under voltage and over voltage differences [3].These voltage variations, especially over voltage are caused due to incorrect tap settings on transformers, single line to ground fault,and abrupt load deduction.Similarly, the under-voltage is due to sudden loading, system faults,starting of large motors, transformer energization, power system equipment malfunction, severe weather, pollution, and capacitor bank switching off [4].The three-phase grid voltage variation of a low voltage network is shown in Fig.1.It can be observed that, the maximum of 465.62 voltage and minimum of 380.56 voltage is recorded in a duration of 7 days.

    Figure 1: Weekly average three-phase grid voltage profiles

    Thus, the performance study of AM under power perturbations becomes more vital.Various types of analysis including magneto-static analysis, time-harmonic analysis, electrostatic analysis,current flow analysis, and thermal stress analysis are performed on AM [5,6].Automated Finite Element Analysis (FEA) is an efficient tool for analyzing various factors of the machine such as, harmonics of flux linkages, skewed rotor geometries and winding layout.This paper employed finite element model using Motor analysis tool in Matlab environment for detailed analysis of AM.Motor analysis tool is an efficient algorithm for designing and analyzing the AM.The performance and behavior of three phase asynchronous machine under 10 percent over-voltage and under-voltage are analyzed.The parameters such as air gap flux, airgap MMF, magnetic field intensity, iron loss, joule loss, relative permeability, torque, efficiency, power factor, and total losses are studied.Conclusively, a detailed comparative analysis is carried out with abnormal and normal operating conditions.Additionally, through sensitivity analysis the change in critical parameters of the system during voltage variation is identified.

    1.1 Problem Overview and Its Significance in the Present Study

    The machine’s performance is degraded by voltage variations in the electrical grid, and it could trip the system, causing economic loss to industries due to outages.The voltage fluctuation has the potential to affect machine characteristics such as overheating, noise, torque oscillation,vibrations, excessive slip losses, and machine efficiency [1].The issues with voltage fluctuation are well outlined in this study for developing a dependable system and enhancing driving characteristics.Infusion craters, severe pitting and fluting damage to machine bearings are caused by frequent voltage variations, resulting in a shorter lifespan [7-9].The machine’s thermal stress is increased by the uneven flux distribution in the airgap, which causes overheating of the rotor cage,as well as rotor material abrasion and dynamic stress on the shaft [10].Asymmetrical operation,pulsing torque, reduced torque, and a decline in machine efficiency result from these losses.The difficulties mentioned above have a greater influence on large and medium industrial machineries.To address these concerns, asynchronous machines are tested using the IEEE 1159-2019 standard,and machine sensitivity is determined using simulation and analytic findings.The test criteria and equipment needed to develop a dependable drive system are clearly outlined in this research.The research contributions are highlighted as follows;

    ? Performance of asynchronous machines under voltage fluctuation using FEM analysis.

    ? Two-dimensional FEM analysis of magnetic flux density, iron loss density, joule loss density,and relative permeability as a function of voltage.

    ? Sensitivity study on asynchronous machines under various voltage circumstances.

    ? The losses of a textile spinning mill motor subjected to a 10% voltage variation were assessed in a case study.

    1.2 Paper Organization

    The article is constructed as; Section 2 deals about generalized magnetic circuit modelling of three phase asynchronous machine with the design parameters.Section 3 describes about the finite element model analysis of 9.2 kW asynchronous machine to explore the magnetic saturation.This section explains how the FEM algorithm works in detail.Section 4 interprets the finite element analysis for a 9.2 kW asynchronous machine and determines its performance under voltage fluctuation.This section also includes a two-dimensional FEM study of magnetic flux density, iron loss density, joule loss density, and relative permeability as a function of voltage.The sensitivity study of asynchronous machines under various voltage circumstances is presented in Section 5.A case study for a textile spinning mill with a 10% voltage fluctuation was developed in Section 6,and their losses were evaluated.Finally, the paper is concluded in Section 7.

    2 Magnetic Circuit Modelling of Three-Phase Asynchronous Machine

    For the design of asynchronous devices, this session presents a simplified magnetic circuit model.Traditional methods to asynchronous machine design have been used, but they may be inadequately accurate or simplified while taking pole-slot counts into account.As a consequence,the early design stage’s layout precision will suffer.Furthermore, since it is hard to ascertain the flux in individual stator teeth, magnetic saturation is often skipped or balanced by correction factors in generalized models.Using the established FEM model, the flux provided by stator winding currents can be measured efficiently and consistently while taking saturation into consideration.The stator side windings are similar with sinusoidally distributed type windings, placed 120?apart,where Nsand Rsare the equivalent turns and resistance of the stator, respectively [6,11].A typical squirrel-cage asynchronous rotor has vertical bars that are made using copper conductor.To permit current to flow through the system, the bars are fixed to the top and bottom conductive rings.The differing currents in the stator windings generate the spinning magnetic field radially across the rotor.Eddy currents in the vertical bars of the squirrel cage cause it to rotate due to the spinning magnetic field [4].The rotor and magnetic field rotate at different speeds.Other assumptions for deriving the equations that describe with the uniform airgap and friction is neglected.The flux linkages, as noticed in the voltage equations, can be stated in terms of inductances.Lsris the magnitude of the mutual inductances between the windings.

    The flux linkages can now be stated as

    The voltage equations in terms of machine variables referred to the stator winding is given by,

    The quantitative accuracy of the air-gap magnetic flux density distribution, however, becomes inadequate when the degree of magnetic saturation by the leakage flux is high [6].When the motor is driven with a high slip by a grid voltage source, the line current rises many times that of the rated state.Therefore, during high slip operation, the leakage flux should be considered for estimating the motor characteristics.For an instance, motor starting, pull-out torque operation,the large amount of flux leakage occurs in AM.Magnetic saturation due to leakage flux is depicted in Fig.2 as a magnetic model.This circuit model will account for magnetic saturation caused by leakage flux.This magnetic model is used to measure the constants on the T-type equivalent circuit.When designing induction motors, FEA and algebraic equations research based on the permeance model are used to study the harmonic components of the air-gap flux density distribution [12,13].

    The flux density (B) in air gap is given as

    Here, F-MMF and P-Permeance-.The total F is expressed as

    Here,FS(1,2,3)is the stator MMF, andFR(1,2,3)is the rotor MMF.The gap reactance Rg is expressed as

    Figure 2: Magnetic circuit modelling of three-phase asynchronous machine

    Here, n1-stator turns per phase, fw1-stator fundamental winding factor, g-air-gap length, and kc-carter’s coefficient.P-Number of poles, and m-number of phases.Dscis the inside diameter of stator core, c is the axial length of rotor core andμ0is the permeability.The torque equation for a three-phase AM is widely known and is given by

    An electric motor consists of the electrical and mechanical system and requires an equation that couples both.

    Here, the number of poles is denoted by P, a moment of inertia is denoted by J, torque linked to the shaft is denoted by TL,and the angular velocity of the rotor is denoted byωr.The specifications and design results of 3-phase asynchronous machine of the AM are given in Tab.1.

    3 Finite Element Analysis of Asynchronous Machine

    Recently, Finite Element Analysis (FEA) has been mostly used in the machine’s analysis and design.To explore the magnetic saturation of 3-phase AM through analyzing the air gap flux component with respect to voltage variation, the finite element model is generated.The voltage variations hardly cause saturation of the entire magnetic circuit; besides, it could produce local saturation in the stator tooth of the motor [10].Later, two FEA procedures were used to make a direct analysis of the induction motor, and their results were compared to establish the better of the two procedures [6].The work was carried out in steady-state conditions with constant supply voltage.The behavior of a three-phase AM with unbalance voltage through phase frame analysis is analyzed in [14].In recent equivalent circuit analysis of AM was carried out using a combination of two fine element models to calculate the equivalent circuit parameters under high magnetic saturation is discussed in [15].Many other attempts to analyze the AM in 2D using FEM, such as various frequency conditions, the influence of different soft magnetic materials, motor with three different modelling, and flux density over the 2D view of the machine,according to a literature review.However, no clear study has been done on motor efficiency in over voltage and under voltage situations.When opposed to 3D analysis, researchers believe that the smaller the electric machine measurements, the greater the error in output assessment.However, other than knowing the degree of error resulting from the use of 2D analysis, there is no concrete evidence for this assertion [13,16].As a result, in order to provide a more accurate response to this claim, electrostatic FEM analysis of a selected AM was performed in 2D analysis in this investigation.The degree of error arising from the use of 2D regression is seen to be important.As a result, it is advised that no additional computing workload be incurred as a result of using 3D analysis.Furthermore, owing to assembly defects, AM’s production systems suffer from manufacturing tolerances.During undervoltage and overvoltage conditions, this results in undesirable variations in the machine’s output and behavior the equivalent circuit parameters are affected [5].As a result, the authenticity of the electrical specifications of a reference system for prototype machines is questioned.In this paper, the FEA of three-phase asynchronous machine under power quality disturbances is carried out using two different approaches, the first approach is based on the motor’s analogous circuit.To evaluate the parameters of the equivalent circuit,the field solutions replicate the indirect (no load and locked rotor) tests.The direct (load) tests are replicated using the second process, which is based solely on the field solution.It blends the field and circuital solutions, directly imposing the motor voltage.Both methods are used on a four-pole 9.2 kW asynchronous computer with a two-dimensional discretized domain.

    Table 1: Ratings and design results of 3-phase asynchronous machine

    3.1 Equivalent Circuit Analysis

    This analysis consists of a no-load test and load test is based on the calculation of equivalent circuit parameters of the AM.The equivalent circuit is shown in Fig.6 for the specified operating point.Once parameters of the equivalent circuit are found, the user can examine how motor parameters depend on different operating conditions.The no-load test is based on the nonlinear magneto-static FEA assuming zero current in rotor bars during the test.The accuracy of the FEA solution is defined by the convergence tolerance of FEA solution field; the default value is 1e-005.To run the no-load test, at first specify stator current values used during the test in the stator current values field.Since it is assumed that there is no rotor current during the no-load test(I2= 0), the stator current is basically equal to the magnetizing current.The load test is based on the nonlinear time-harmonic FEA.The accuracy of the FEA solution is defined by convergence tolerance of the FEA solution field; the default value is 1e-005.The RMS phase current, supply frequency, and rotor speed (or slip) is defined.By default, the speed is set to zero, so the standard locked rotor test can be performed.During the load test, the parameters of the equivalent circuit is computed.The motor parameters are computed assuming that L1, L2, R1, R2, and Rmdo not change while Lmchanges depending on the magnetizing current according to no-load test data.In reality, L1, L2, R2, and Rmvalues more or less depending on rotor speed and saturation.AC Analysis is based on the time-harmonic FEA and the equivalent circuit but provides more accurate results comparing with equivalent circuit analysis since it can take into account variation of equivalent circuit parameters at different operating conditions.AC Analysis is performed for the specified operating points defined by the power supply (3-phase voltage source or 3-phase current source), RMS phase voltage (RMS phase current), and supply frequency fields.The speed can be defined whether by the number of rotations per minute or by slip.

    3.2 Finite Element Method

    Time harmonic finite element method is centered on the assumption of the sinusoidal timedependence of the field, which allows the computation time to be reduced radically comparing with the stepping FEM [17].This method provides only steady-state motor parameters and cannot be used for the analysis of transient regimes.Time-stepping finite element method is the representation of the time derivative of the magnetic vector potential for each time step and can be used for the analysis of both transient and steady-state [18].In both methods, the discretization of the problem over the plane of the machine’s cross-section using FEM results in the system of linear equations written form as follows;

    Here K is the stiffness matrix, X is the vector of unknowns including values of the magnetic vector potential in mesh node, rotor and stator currents, and voltages, F is right side vector.The residual vector of the gauss ?Xis defined as;

    Here, tol-convergence is the tolerance setting the desired accuracy of the solution, the generalized algorithm used in motor analysis for a nonlinear instance of the time-stepping process is shown in Fig.3.An outer loop of the algorithm represents integration over time with stepΔt.The magnetic field issue is addressed on each integration step in the nonlinear solver loop.

    Figure 3: Equivalent circuit of the AM

    Gauss-newton iteration continues until the inequality is satisfied [19].In motor analysis,the first-order triangular finite element mesh is used.Every time while the time-stepping FEA simulation is in progress when the rotor position is changed, the mesh should also be changed.It is implemented by rotation of the rotor mesh connected to the fixed stator mesh via a sliding layer located in the air gap.In motor analysis, the air gap always has an odd number of mesh layers(3, 5, 7, or 9), and the sliding layer is always located at the center of the air gap, as shown in Fig.4.To reduce the number of finite elements, and as a consequence reduce computation time,periodic/antiperiodic boundary conditions can be applied.Using this type of boundary conditions is based on the periodicity of the magnetic field of the motor [5].The mesh model is shown in Figs.5a and 5b.

    Figure 4: Time-stepping FEM algorithm

    Figure 5: Finite element mesh model.(a) 2D FEM model of 9.2 kW AM; (b) multi-slice FEM model

    Figure 6: Parameter variations with respect to slip; (a) 10% undervoltage (b) 10% overvoltage

    4 Interpretation of Finite Element Analysis

    The paper analysis the three-phase asynchronous motor, fed by a symmetrical three phase voltage sources with 10 percent variation.The two different FEA approaches, showing their main characteristics; 1.FEA method for the indirect test simulation2.FEA method for the load tests simulation.Both methods use fixed mesh and eddy current models, they do not consider the field harmonics and they are applied to the same two-dimensional (2D), domain ‘D’.

    The field components are supposed varying sinusoidally in time.Thus in a Cartesian coordinate system (x, y, z) the values of the generic quantity G(P, t), in the point P = (x, y, z) and at the instant t, can be written as;

    Here, i = x, y, z;ωis the frequency in rad/s andθdenotes the electric phase difference in rad.

    It is supposed that the current density vector J and the magnetic potential vector A have only the z-axis components, i.e., J = (0, 0, Jz) and A = (0, 0, Az).Thus, magnetic field strength H and flux density B have only components, on the (x, y) plane perpendicular to the z-axis, i.e., H= (Hx, Hy, 0) and B = (Bx, By, 0).In this way the induction motor analysis may be realised by a 2D finite element analysis, simpler than a three-dimensional (3D) one, without a significant loss of accuracy.The 3D effects have been by suitable circuit elements outside the field solution.The sensitivity of the variables to change in voltage is studied in the sensitivity analysis part of the paper.

    4.1 Performance of Asynchronous Machine Under Voltage Variations

    In this first case, the induction motor is modelized by a simplified equivalent circuit Fig.3.FEA method is used to reproduce the usual indirect tests, that are no load and locked rotor tests.From the field solution the non-linear parameters of the simplified equivalent circuit, essential to evaluate the motor performance, are, calculated.In order to evaluate the total losses and the motor efficiency, R0 the stator resistance R1 is evaluated and presented in Figs.3 and 6.Similarly,in order to analyse the load operation, the non-linearity of magnetic materials, the current induced in the rotor bars and the rotor movement have to be considered simultaneously.The field problem is a nonlinear and time-invariant one.It presents both motional and induced current density term.The simulation results for the standard asynchronous motor design with 10 percent undervoltage and 10 percent overvoltage is presented in Tab.2 with 6 percent slip (rated load).As well, the parameter profiles during no load to full load variation with 10 percent undervoltage and 10 percent overvoltage symmetrical supply is shown in Figs.6a and 6b respectively.It is clearly observed that there is an increase in losses as the supply voltage increase and a reduction in losses when voltage is lowered 10 percent below the nominal voltage.The reduction in losses results in increased efficiency, whereas the efficiency decreases by 10 percent OV.Upon comparison, it can be mentioned that the performance variables of stator winding losses and efficiency is highly affected by the change in voltages.Whereas, torque, input power, and output power are moderately affected by voltage variations.

    Table 2: Effect of voltage variations on asynchronous machines

    (Continued)

    Table 2: Continued

    4.2 Two-Dimensional FEM Analysis

    The objective of the analysis is to determine the extent to which flux can be routed in the asynchronous machine to achieve an electronically adjustable output voltage.The flux distribution will be observed around the windings and in each slot.The method chosen for analysis will be simulation using the Finite Element Method.The FEM software package Motor Analysis-IM is used.Ideally a method is required that allows analysis of nonlinear magnetic circuits with three voltage level.Airgap distribution plot allows the view of the distribution of the different parameters such as flux and MMF over the machine air gap beside the air gap harmonic components.The results are obtained for a selected power supply, supply frequency, and slip/speed.Fig.7 shows the airgap flux distribution, airgap MMF distribution for 10 percent under voltage, 0 percent voltage variation and 10 percent over voltage respectively.It is clearly noticed that the magnitude of the airgap flux and MMF are slightly affected by variations in voltage.The corresponding harmonic component magnitude is also observable.Undervoltage leads to a reduction of airgap flux, and overvoltage leads to a mild increase in flux magnitude, but the flux distribution pattern remains the same.In this paper, more than 800 simulations were performed over 120 degrees of rotor motion.A list of flux linkagevs.time was numerically differentiated with respect to time to yield line-to-neutral voltage.The difference between line-to-neutral voltages was taken to yield line-to-line voltage.Both a Mathematica notebook and a Matlab script performing these calculations.

    Figure 7: Airgap flux and MMF distribution with harmonic order (a) flux distribution with 10%under voltage; (b) MMF distribution with 10% under voltage; (c) flux distribution with 0% voltage variation; (d) MMF distribution with 0% voltage variation; (e) flux distribution with 10% over voltage; (f) MMF distribution with 10% over voltage

    The full cross-section view has been zoomed to show the pattern of distribution of flux lines in the machine.Flux lines are plotted so that the closeness of the flux lines indicates the magnitude of the flux density x-axis and y-axis limits will be determined by the size of the machine’s cross-section.Z-axis limits field sets the color scale limits to a specified minimum and maximum values.Values between z-axis limits are linearly mapped to the used color scale (color map).Data values less or greater than specified z-axis limits are mapped to the minimum limit or to the maximum limit, respectively.It is clearly visible from the cross-section plot of magnetic flux density that the magnetic flux density increases with an increase in voltage and decreases under voltage.This is evident from the color scaling.During 10 percent under voltage, the flux density is less in the stator and rotor when compared to the 10 percent over voltage, where the flux density value is higher viewable from the color change of area with pale green to bright yellow under 10 percent over voltage.The increase in flux density also results in an increase in iron losses.The iron losses increase and decrease with overvoltage and under-voltage, respectively, as evident from the color change of the iron loss plot from light blue to greenish-blue in Fig.8.A possible solution for reducing iron losses will be the choice of good magnetic material with lower specific iron losses.

    Figure 8: 2D color scale analysis of magnetic flux density-Z-axis of AM (a) 10% under voltage;(b) 0% voltage variation; (c) 10% over voltage

    The joule loss density is the copper losses that occur in the windings of the IM.It is also called the I2R loss and is the main source of heat generation in the AM.The joule loss is very sensitive to variation in supply voltage.There is a sharp rise in copper losses under overvoltage and a reduction in losses under 10 percent UV.This can be clearly observed from Fig.9, where the color change in the stator is evident with fewer losses in 10 percent UV and heavy losses in 10 percent OV.Also, the same can be verified with values from Tab.2, where the stator winding losses reduces to 332 W from the nominal value of 635.5 W under 10 percent UV and increases to 1210.77 W under 10 percent OV.The huge increase and decrease in losses lead to lower efficiency under over-voltage and higher efficiency at under voltage.From Fig.9,it can be noticed that the relative permeability increases with over-voltage and reduces in 10 percent UV.High relative permeability can degrade the performance of the machine.For relative permeability of the machine can be altered by using different slot wedges or by using different magnetic materials.

    Figure 9: 2D color scale analysis of machine loss densities and relative permeabilities Iron loss density on (a) 10% under voltage; (b) 0% voltage variation; (c) 10% over voltage; Joule loss density on (d) 10% under voltage; (e) 0% voltage variation; (f) 10% over voltage; relative permeability on(g) 10% under voltage; (h) 0% voltage variation; (i) 10% over voltage

    5 Sensitivity Analysis

    A sensitivity analysis is performed in this paper under various voltage conditions.This research is carried out to see how design variables (such as the air gap) affect the electrical parameters [5,20].Hence, it is proceeded with the analysis, observing how an induction motor performance variable are affected by 10% under-voltage (UV), 10% overvoltage (OV), and 0%voltage variations (VV).Further, sensitivity analysis on the various performance variables is carried out, and the results have been used to identify whether 10% OV or 10% UV impacts the AM most.Sensitivity analysis, also known as simulation analysis, is used to measure the impact of fluctuations on the performance of the system.The analysis involves the measurement of the performance variables of the IM under variations in the input parameter viz.10% OV and 10% UV to establish the parameter that affects the performance of the AM most.In industries,necessities related to the quality and control of AM are essential.In order to predict the influence of variation in motor parameters on the machine performance, sensitivity analysis is employed by machine engineers and designers.Performance variables such as torque, current, input power,efficiency, total loss, rotor cage loss, stator winding loss, output power, power factor, and iron core loss have been investigated by means of the Motor analysis tool in MATLAB.The sensitivity index of the various performance variables of the AM can be calculated using the relation

    where,Nx-performance parameter when theαis increased or lowered by 10%, and Nr is the rated value of the performance parameter.The sensitivity index of performance variables for 6% slip and 50 Hz supply frequency is tabulated in Tab.3.

    Figure 10: Performance of induction motor during transient period (a) magnetic Flux density plot;(b) energy density plot; (c) Air gap distribution plot

    The supply voltage variation has a considerable effect on almost all performance variables.The performance variables like stator winding loss and total loss are highly sensitive with respect to 10% OV and moderately sensitive to 10% UV.Furthermore, the sensitivity of phase current,power factor, and input power are moderately sensitive with respect to 10% UV and 10% OV.Torque, output power, efficiency, rotor cage loss, and iron loss are less sensitive with respect to 10% UV and 10% OV.From the above observation, it is clear that 10% OV has a maximum effect on the performance variables.Hence, we can conclude that the AM is prominently affected more by overvoltage than under voltage.Transient Analysis is based on the time-stepping FEA and has the highest accuracy among other analysis methods allowing to take into account rotation of the rotor, cogging torque, higher harmonics, PWM switching and etc.However, this type of analysis takes considerably more time comparing with previous two analysis types (equivalent circuit and AC analysis).Transient Analysis uses Simulation script file—MATLAB-function which is called on each simulation time step and allows the user to change all simulation settings.The magnetic flux density and energy density distribution plots are shown in Figs.10a and 10b respectively.The air gap distribution of the MMF at time point of 1 s in the first slice’s cross-section is shown in Fig.10c.

    Table 3: Sensitivity index of performance variables

    Figure 11: Textile spinning frame with the load diagram (a) spinning frame with drive motor (b)load diagram of spinning drive motor

    6 Case Study on Textile Spinning Mill with 10 Percent Voltage Variation

    Here, ENomis total energy losses under nominal condition (0% VV) in kWh, EUVis total energy losses in under voltage condition (10% UV) in kWh, EOVis total energy losses in under voltage condition (10% OV) in kWh, n-number of samples,x-process time in hours,yNom-power losses under nominal condition (0% VV) in watts,yUVis power losses in under voltage condition(10% UV) in watts,YOVis power losses in over voltage condition (10% OV) in watts (Fig.12).

    Figure 12: Total losses under different voltages

    In this part, a spinning machine in the textile industry is considered and studied in terms of energy conservation analysis.The textile mill uses a variety of processes such as combing,carding, spinning, and so on, with the spinning process accounting for 50% of overall energy consumption [6].The spinning process, depicted in Fig.11a, is a mechanized method of obtaining thread from cotton fibers.As shown in Fig.11b, the amount of thread on the spindle increases slightly from zero to the maximum permissible limit during this operation.As a result, the load on the spinning drive motor grows exponentially from very light to rated value.

    Taking the process cycle of the spinning frame in a textile manufacturing driven with 9.2 kW.From the results, noted that the machine with 10 percent over voltage creates 18,93,672 kW losses per year, which is almost 101.4 percent losses increased from the nominal value.Whereas, the AM with 10 percent under voltage creates 4,80,217 kW losses per year, which is almost 48.9 percent losses decreased from the nominal value.

    7 Conclusion

    A detailed analysis of the AM under the effect of power quality perturbations of 10%overvoltage and 10% under voltage was done.The effect of variation in voltage on the various performance variables and the distribution of magneto-static quantities in the airgap and cross-section have been studied in detail.Further, sensitivity function was used to determine the sensitivity analysis of the performance parameters like input power, efficiency, current, power factor, torque,power factor, output power, cage, iron, stator winding, and total loss of AM corresponding to voltage variations by using FEA simulations with the help of motor analysis tool in MATLAB.By sensitivity analysis, the impact of parameter variations in IM is analyzed.The 10% UV and 10% OV have a considerable effect on almost all performance variables.The AM is worst affected by 10% OV than 10% UV.An overall view of the AM behavior was acquired, and the results of the analysis are very useful in the design of an AM that can withstand the expected variation of the parameter, which in turn can lead to giving the best performance, long machine life, less machine damage to AM and satisfy the load for various application.

    Acknowledgement:This research was funded by the Deanship of Scientific Research at Princess Nourah bint Abdulrahman University through the Fast-track Research Funding Program.

    Funding Statement:This research was funded by the Deanship of Scientific Research at Princess Nourah bint Abdulrahman University through the Fast-track Research Funding Program.

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

    亚洲av成人av| 18禁在线播放成人免费| 久久久色成人| 最近最新中文字幕大全电影3| 黄色一级大片看看| 国产精品福利在线免费观看| 免费高清视频大片| 一本久久中文字幕| 在线播放国产精品三级| 色5月婷婷丁香| 一进一出抽搐gif免费好疼| 他把我摸到了高潮在线观看| 亚洲黑人精品在线| 校园春色视频在线观看| 亚洲欧美日韩东京热| 日韩欧美国产一区二区入口| 天堂√8在线中文| 免费搜索国产男女视频| 亚洲男人的天堂狠狠| 老司机深夜福利视频在线观看| 亚洲av中文av极速乱 | 日日夜夜操网爽| 亚洲电影在线观看av| 99热网站在线观看| 欧美激情在线99| 51国产日韩欧美| 国内精品久久久久精免费| 亚洲成人免费电影在线观看| 天堂动漫精品| 精品一区二区三区视频在线| 中文字幕久久专区| 少妇高潮的动态图| 精品一区二区三区人妻视频| 国内精品美女久久久久久| 一级黄色大片毛片| 亚洲中文字幕一区二区三区有码在线看| 成人毛片a级毛片在线播放| 少妇人妻一区二区三区视频| 丰满人妻一区二区三区视频av| 91久久精品国产一区二区三区| 69av精品久久久久久| 亚洲av免费高清在线观看| 久久久久性生活片| 内地一区二区视频在线| 欧美最黄视频在线播放免费| 欧美+日韩+精品| 国产精品久久久久久久电影| 精品久久久久久,| АⅤ资源中文在线天堂| 美女xxoo啪啪120秒动态图| 少妇猛男粗大的猛烈进出视频 | 亚洲av二区三区四区| 国产熟女欧美一区二区| 我的老师免费观看完整版| 国产三级在线视频| 男女做爰动态图高潮gif福利片| 亚洲av免费在线观看| 午夜亚洲福利在线播放| 免费观看的影片在线观看| 国产视频内射| 国产在视频线在精品| 中文字幕av在线有码专区| 看十八女毛片水多多多| ponron亚洲| 婷婷六月久久综合丁香| 久久久精品欧美日韩精品| 亚洲欧美日韩卡通动漫| 很黄的视频免费| 99热这里只有精品一区| 可以在线观看毛片的网站| 看免费成人av毛片| 两个人的视频大全免费| 日韩中字成人| 国国产精品蜜臀av免费| 黄色日韩在线| 99久久久亚洲精品蜜臀av| 亚洲精品成人久久久久久| 国产男人的电影天堂91| а√天堂www在线а√下载| 春色校园在线视频观看| 国产v大片淫在线免费观看| 日韩大尺度精品在线看网址| 久久久久国产精品人妻aⅴ院| 国产老妇女一区| av中文乱码字幕在线| 99热只有精品国产| 校园人妻丝袜中文字幕| 国产人妻一区二区三区在| 亚洲avbb在线观看| 日日摸夜夜添夜夜添小说| 亚洲av五月六月丁香网| 一进一出抽搐gif免费好疼| 国产蜜桃级精品一区二区三区| 91久久精品国产一区二区三区| 99riav亚洲国产免费| 婷婷丁香在线五月| 中文字幕av在线有码专区| 国产午夜福利久久久久久| 一a级毛片在线观看| 欧美性猛交黑人性爽| 1024手机看黄色片| 精品一区二区三区av网在线观看| 亚洲经典国产精华液单| 日韩中字成人| 22中文网久久字幕| 国产午夜精品论理片| 午夜福利在线在线| 欧美日韩黄片免| 深夜精品福利| 欧美成人a在线观看| 免费在线观看成人毛片| 内地一区二区视频在线| 成人欧美大片| 一进一出抽搐gif免费好疼| 久久精品久久久久久噜噜老黄 | 欧美精品国产亚洲| 日韩欧美精品免费久久| 51国产日韩欧美| 国产aⅴ精品一区二区三区波| 看片在线看免费视频| 国产黄片美女视频| 他把我摸到了高潮在线观看| 搡老妇女老女人老熟妇| 亚洲成人免费电影在线观看| 国产又黄又爽又无遮挡在线| 免费看光身美女| 亚洲一区二区三区色噜噜| 国内揄拍国产精品人妻在线| 国产高清三级在线| 亚洲三级黄色毛片| 亚洲三级黄色毛片| 亚洲精华国产精华精| 一个人观看的视频www高清免费观看| 国产淫片久久久久久久久| av黄色大香蕉| 两性午夜刺激爽爽歪歪视频在线观看| 无遮挡黄片免费观看| 我的女老师完整版在线观看| 久久6这里有精品| 国产麻豆成人av免费视频| 久久久久久久久久久丰满 | 不卡一级毛片| 我的女老师完整版在线观看| 国产精品久久久久久av不卡| 日韩一区二区视频免费看| 人人妻,人人澡人人爽秒播| av视频在线观看入口| 亚洲国产精品久久男人天堂| 国内精品美女久久久久久| 国产一区二区三区在线臀色熟女| 成熟少妇高潮喷水视频| 亚洲成人久久性| 欧美一区二区亚洲| 观看免费一级毛片| 亚洲专区中文字幕在线| 亚洲av第一区精品v没综合| 尤物成人国产欧美一区二区三区| 99国产极品粉嫩在线观看| 亚洲国产精品成人综合色| 亚洲第一电影网av| 波多野结衣高清无吗| 日本a在线网址| 美女高潮的动态| 不卡视频在线观看欧美| 91麻豆精品激情在线观看国产| 简卡轻食公司| 色综合色国产| 国产成人福利小说| 蜜桃久久精品国产亚洲av| 国产一区二区三区视频了| 99九九线精品视频在线观看视频| 久久久久精品国产欧美久久久| 午夜福利高清视频| 亚洲aⅴ乱码一区二区在线播放| 观看免费一级毛片| 91av网一区二区| 日韩欧美精品免费久久| 国产精品人妻久久久久久| 男女下面进入的视频免费午夜| 精品欧美国产一区二区三| 欧美日本视频| 国产精品久久久久久精品电影| 国产真实乱freesex| 两性午夜刺激爽爽歪歪视频在线观看| 日本精品一区二区三区蜜桃| 性欧美人与动物交配| 小说图片视频综合网站| 一夜夜www| 成人美女网站在线观看视频| 一级a爱片免费观看的视频| 国产亚洲精品久久久com| 午夜视频国产福利| 黄色女人牲交| 久久久国产成人免费| 一个人看视频在线观看www免费| 色吧在线观看| 一级av片app| 999久久久精品免费观看国产| 18禁在线播放成人免费| 窝窝影院91人妻| 亚洲精品成人久久久久久| 一本一本综合久久| 白带黄色成豆腐渣| 国产精品日韩av在线免费观看| 婷婷丁香在线五月| 国产乱人视频| 久久九九热精品免费| 老司机深夜福利视频在线观看| 国产精品人妻久久久影院| 国产单亲对白刺激| 22中文网久久字幕| a在线观看视频网站| 91狼人影院| 88av欧美| 欧美日韩亚洲国产一区二区在线观看| 中文字幕av成人在线电影| 久久亚洲真实| 国产亚洲精品av在线| 婷婷六月久久综合丁香| 桃色一区二区三区在线观看| 国内精品美女久久久久久| 狠狠狠狠99中文字幕| 免费人成在线观看视频色| 久久精品国产自在天天线| 一本一本综合久久| 国产精品野战在线观看| 日日啪夜夜撸| 高清毛片免费观看视频网站| 级片在线观看| 欧美日本亚洲视频在线播放| 麻豆国产av国片精品| 欧美日韩瑟瑟在线播放| 久久精品综合一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 国产精品久久久久久亚洲av鲁大| 男人和女人高潮做爰伦理| 国产探花极品一区二区| 成年女人永久免费观看视频| 十八禁国产超污无遮挡网站| 日本 av在线| 很黄的视频免费| 动漫黄色视频在线观看| 69人妻影院| 日日夜夜操网爽| 久久这里只有精品中国| 免费在线观看成人毛片| 一进一出抽搐gif免费好疼| 性插视频无遮挡在线免费观看| 免费av观看视频| 在线免费观看的www视频| 成人鲁丝片一二三区免费| 两个人的视频大全免费| 99久久精品国产国产毛片| 一个人看视频在线观看www免费| 长腿黑丝高跟| av黄色大香蕉| 精品国产三级普通话版| 国产中年淑女户外野战色| 亚洲av二区三区四区| 美女高潮喷水抽搐中文字幕| 男女之事视频高清在线观看| 婷婷色综合大香蕉| 长腿黑丝高跟| 欧美日韩瑟瑟在线播放| 男女下面进入的视频免费午夜| 日本a在线网址| 男人狂女人下面高潮的视频| 亚洲av二区三区四区| 国产成年人精品一区二区| 国产精品嫩草影院av在线观看 | 欧美日韩瑟瑟在线播放| 99热网站在线观看| 精品一区二区三区视频在线| 18+在线观看网站| 无遮挡黄片免费观看| 天天躁日日操中文字幕| 久久国产精品人妻蜜桃| 免费在线观看影片大全网站| 中文字幕熟女人妻在线| 人妻久久中文字幕网| 午夜福利在线观看免费完整高清在 | 国产成人a区在线观看| 精华霜和精华液先用哪个| 51国产日韩欧美| 精品一区二区三区视频在线观看免费| 国产成人影院久久av| 又黄又爽又免费观看的视频| 天美传媒精品一区二区| 久久久久久久久久成人| 美女xxoo啪啪120秒动态图| 窝窝影院91人妻| 久久中文看片网| 亚洲人成网站在线播| 最好的美女福利视频网| 国内精品久久久久久久电影| 国内少妇人妻偷人精品xxx网站| 桃红色精品国产亚洲av| 一级黄色大片毛片| 中国美女看黄片| 黄色配什么色好看| 最近中文字幕高清免费大全6 | 久久人人爽人人爽人人片va| 男女那种视频在线观看| 国产私拍福利视频在线观看| 尤物成人国产欧美一区二区三区| 九九久久精品国产亚洲av麻豆| 99久国产av精品| av在线观看视频网站免费| 国产成人aa在线观看| 不卡一级毛片| 最近在线观看免费完整版| 欧美区成人在线视频| 99久久无色码亚洲精品果冻| 国产精品国产高清国产av| 制服丝袜大香蕉在线| 色综合婷婷激情| 有码 亚洲区| 精品久久久噜噜| 久久久久久伊人网av| 中文字幕免费在线视频6| 久久亚洲精品不卡| 99热6这里只有精品| 精品一区二区三区人妻视频| 麻豆国产av国片精品| 日韩中文字幕欧美一区二区| 高清毛片免费观看视频网站| 日本黄色视频三级网站网址| 禁无遮挡网站| 亚洲国产高清在线一区二区三| 一区二区三区免费毛片| 在线a可以看的网站| 日日干狠狠操夜夜爽| 国产精品99久久久久久久久| 免费大片18禁| 内射极品少妇av片p| av在线天堂中文字幕| 日韩一区二区视频免费看| 久久人人爽人人爽人人片va| 国产欧美日韩精品一区二区| 有码 亚洲区| 大又大粗又爽又黄少妇毛片口| 免费不卡的大黄色大毛片视频在线观看 | 国产在线精品亚洲第一网站| 在线免费观看的www视频| 精品人妻一区二区三区麻豆 | 99在线人妻在线中文字幕| 级片在线观看| 免费av毛片视频| 亚洲国产精品久久男人天堂| 亚洲综合色惰| 99热这里只有是精品50| 欧美精品国产亚洲| 国产伦精品一区二区三区四那| 亚洲av成人av| 亚洲av中文字字幕乱码综合| 国产一区二区三区在线臀色熟女| 色5月婷婷丁香| 极品教师在线免费播放| 两性午夜刺激爽爽歪歪视频在线观看| 午夜亚洲福利在线播放| 久久久久国产精品人妻aⅴ院| 亚洲七黄色美女视频| 亚洲va日本ⅴa欧美va伊人久久| 欧美一区二区亚洲| 成人高潮视频无遮挡免费网站| 69人妻影院| 欧美精品啪啪一区二区三区| 国产精品久久久久久av不卡| 亚洲人成网站在线播放欧美日韩| 美女cb高潮喷水在线观看| 国产视频一区二区在线看| 久久天躁狠狠躁夜夜2o2o| 三级男女做爰猛烈吃奶摸视频| 啪啪无遮挡十八禁网站| av福利片在线观看| 毛片女人毛片| 亚洲成人久久性| 国产精品一区www在线观看 | 亚洲国产精品久久男人天堂| 免费电影在线观看免费观看| 免费在线观看日本一区| 九九热线精品视视频播放| 亚洲中文日韩欧美视频| 一级毛片久久久久久久久女| 婷婷亚洲欧美| 在线观看午夜福利视频| 成年女人看的毛片在线观看| 简卡轻食公司| 免费看a级黄色片| 女人被狂操c到高潮| 亚洲精品成人久久久久久| 欧美日韩亚洲国产一区二区在线观看| www.色视频.com| 国产色爽女视频免费观看| 丰满人妻一区二区三区视频av| 亚洲人成网站高清观看| 久久人人爽人人爽人人片va| 亚洲精华国产精华液的使用体验 | 亚洲熟妇熟女久久| 真人做人爱边吃奶动态| 少妇裸体淫交视频免费看高清| 欧美绝顶高潮抽搐喷水| 嫩草影院入口| 日韩中文字幕欧美一区二区| 国产精品一区www在线观看 | 精品人妻一区二区三区麻豆 | 国产精品女同一区二区软件 | 午夜福利视频1000在线观看| 麻豆久久精品国产亚洲av| 亚洲国产色片| 国产精品福利在线免费观看| 欧美成人免费av一区二区三区| 精华霜和精华液先用哪个| 亚洲av不卡在线观看| 最新中文字幕久久久久| 在线播放国产精品三级| 久久人人爽人人爽人人片va| 久久久久久久久中文| 最近视频中文字幕2019在线8| 91在线观看av| 日本爱情动作片www.在线观看 | 九九久久精品国产亚洲av麻豆| 国产大屁股一区二区在线视频| 99热只有精品国产| 成人精品一区二区免费| 嫩草影院新地址| 欧美不卡视频在线免费观看| 桃红色精品国产亚洲av| 男女啪啪激烈高潮av片| 国产欧美日韩一区二区精品| .国产精品久久| 久久人人爽人人爽人人片va| 少妇人妻精品综合一区二区 | 十八禁网站免费在线| 亚洲av不卡在线观看| 人人妻,人人澡人人爽秒播| 99在线视频只有这里精品首页| 国产高清不卡午夜福利| 亚洲av第一区精品v没综合| 国产亚洲精品综合一区在线观看| 日本精品一区二区三区蜜桃| 亚洲人成网站在线播| 无遮挡黄片免费观看| 成年人黄色毛片网站| 国产老妇女一区| 啪啪无遮挡十八禁网站| 欧美绝顶高潮抽搐喷水| 噜噜噜噜噜久久久久久91| 美女被艹到高潮喷水动态| 亚洲自偷自拍三级| 乱码一卡2卡4卡精品| 国产蜜桃级精品一区二区三区| 国产主播在线观看一区二区| 日韩欧美国产在线观看| 亚洲国产欧美人成| 少妇熟女aⅴ在线视频| 午夜老司机福利剧场| 身体一侧抽搐| 色视频www国产| 国产高清激情床上av| 精品一区二区三区av网在线观看| 久久久精品欧美日韩精品| 日韩精品中文字幕看吧| 亚洲四区av| 亚洲精品在线观看二区| 日韩欧美在线二视频| 88av欧美| 色5月婷婷丁香| 欧美黑人巨大hd| 国产精品一区二区性色av| 亚洲欧美激情综合另类| 婷婷六月久久综合丁香| 国产欧美日韩精品一区二区| 91在线精品国自产拍蜜月| 欧美日韩综合久久久久久 | 在线天堂最新版资源| 一进一出抽搐动态| 色噜噜av男人的天堂激情| 十八禁国产超污无遮挡网站| 舔av片在线| 婷婷亚洲欧美| 精品一区二区三区视频在线| 深夜a级毛片| 中文在线观看免费www的网站| 琪琪午夜伦伦电影理论片6080| 不卡视频在线观看欧美| 亚洲最大成人手机在线| 村上凉子中文字幕在线| 老女人水多毛片| 麻豆国产97在线/欧美| 内射极品少妇av片p| 国产成人影院久久av| 九九在线视频观看精品| 国产成人av教育| 婷婷色综合大香蕉| 国产色婷婷99| 一级毛片久久久久久久久女| 熟女电影av网| 一级黄片播放器| 热99在线观看视频| 免费人成在线观看视频色| 观看美女的网站| 国产熟女欧美一区二区| xxxwww97欧美| 99热网站在线观看| 午夜免费男女啪啪视频观看 | 国产av一区在线观看免费| 国产一区二区激情短视频| 日本熟妇午夜| 欧美色视频一区免费| 深夜a级毛片| 九色国产91popny在线| 丰满乱子伦码专区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产色片| 在线观看一区二区三区| 日日摸夜夜添夜夜添小说| 久久6这里有精品| 悠悠久久av| 99热这里只有是精品在线观看| 久久精品国产亚洲av香蕉五月| 国产伦精品一区二区三区四那| 99国产精品一区二区蜜桃av| 国产精品综合久久久久久久免费| 欧美另类亚洲清纯唯美| 男女边吃奶边做爰视频| 亚洲色图av天堂| 国产淫片久久久久久久久| www日本黄色视频网| 国内精品久久久久精免费| 日本爱情动作片www.在线观看 | 精品久久久久久久久av| 欧美激情久久久久久爽电影| 午夜爱爱视频在线播放| 日本 av在线| 精品久久久久久久人妻蜜臀av| 一个人观看的视频www高清免费观看| 国产精品综合久久久久久久免费| 亚洲久久久久久中文字幕| 久久久久久国产a免费观看| 亚洲av五月六月丁香网| 欧美xxxx性猛交bbbb| avwww免费| 99久久精品一区二区三区| 女人十人毛片免费观看3o分钟| 身体一侧抽搐| 一卡2卡三卡四卡精品乱码亚洲| 久久人人爽人人爽人人片va| 少妇高潮的动态图| bbb黄色大片| 久久精品夜夜夜夜夜久久蜜豆| 真人做人爱边吃奶动态| 一边摸一边抽搐一进一小说| 婷婷精品国产亚洲av在线| 国产 一区 欧美 日韩| 欧美日韩乱码在线| 国产免费一级a男人的天堂| 亚洲av电影不卡..在线观看| 亚洲精品粉嫩美女一区| 国产白丝娇喘喷水9色精品| 国产麻豆成人av免费视频| 久久久成人免费电影| 亚洲国产日韩欧美精品在线观看| 美女cb高潮喷水在线观看| 97人妻精品一区二区三区麻豆| 日韩大尺度精品在线看网址| 少妇猛男粗大的猛烈进出视频 | 国产高潮美女av| 色av中文字幕| 国产一区二区亚洲精品在线观看| 99久久九九国产精品国产免费| 51国产日韩欧美| 久久久久久伊人网av| 欧美+日韩+精品| 国产淫片久久久久久久久| 亚洲美女视频黄频| 中文资源天堂在线| 久久午夜亚洲精品久久| 精品久久久久久久久亚洲 | 亚洲黑人精品在线| 国产精品久久久久久亚洲av鲁大| 国产色爽女视频免费观看| 亚洲七黄色美女视频| 男女啪啪激烈高潮av片| 国产大屁股一区二区在线视频| 男女之事视频高清在线观看| 午夜福利高清视频| 级片在线观看| 最近中文字幕高清免费大全6 | 午夜福利视频1000在线观看| 婷婷色综合大香蕉| 看黄色毛片网站| 俺也久久电影网| 精品久久久久久久久亚洲 | 欧美成人a在线观看| av中文乱码字幕在线| 国内精品久久久久精免费| 性插视频无遮挡在线免费观看| 亚洲乱码一区二区免费版| 联通29元200g的流量卡| videossex国产| 国产一区二区三区av在线 | 99热这里只有是精品50| 69av精品久久久久久| 国产av在哪里看| 亚洲三级黄色毛片| 欧美xxxx黑人xx丫x性爽| 国产午夜精品久久久久久一区二区三区 | 中文字幕免费在线视频6| 我要看日韩黄色一级片| 免费在线观看影片大全网站| 亚洲第一电影网av| 日本爱情动作片www.在线观看 | 欧美一区二区亚洲| 一区二区三区高清视频在线|