999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

Optimization and Performance Comparison of Hairpin-Winding PMSM for Electric Vehicles Under Drive Cycle

2021-11-23 08:58:54,,,

,,,

1.School of Electrical Engineering,Beijing Jiaotong University,Beijing 100044,P.R.China;

2.Electric Engineering Department,Suzhou Inovancel Technology Co.,Ltd,Suzhou 215100,P.R.China

Abstract:Aiming at the problem of large AC copper loss caused by skin effects and proximity effects,and low efficiency at high speed of the hairpin-winding permanent magnet synchronous motor(PMSM)for electric vehicles(EVs),this paper firstly established the electromagnetic analytical model of the hairpin winding to calculate AC resistance.And the finite element model(FEM)of the hairpin-winding driving motor is established to calculate the AC characteristic of the hairpin winding at different speeds and temperatures.Then,combining modified particle swarm optimization(MPSO)and FEM,a 60 kW hairpin-winding PMSM is optimized under driving cycle conditions,and the electromagnetic performance and heat dissipation performance are compared with that of the traditional strand-winding motor.Finally,a prototype is made and an experimental platform is built to test the efficiency Map and temperature rise of the hairpin-winding motor over the whole speed range and verify the accuracy of the proposed optimization design method.The results show that the hairpin-winding PMSM not only has higher slot filling rate,high-efficiency range and power density,but also has better heat dissipation performance,which is suitable for application in the field of electric vehicles.

Key words:electric vehicles(EVs);hairpin-winding;optimization design;drive cycle;performance comparison

0 Introduction

In the new global economy,the electric vehicles(EVs)play an important role in addressing the issue of increasingly environmental pollution.As the key parts of electric vehicles,the electromagnetic and heat dissipation performance of the driving motor directly affect the operating quality and cruising ability of electric vehicles.Therefore,the improvement of the driving motor’s power density and efficiency has received considerable critical attention[1-3].

Recently,hairpin windings,also known as form-wound windings have attracted the attention of motor researchers because of its high slot filling rate and better heat dissipation performance[4-6].Ref.[7]compared the efficiency and continuous powerspeed trajectory of the traditional strand-winding and the hairpin-winding induction motors used for EVs in the full speed range,and found that the hairpin-winding induction motor has a higher efficiency range than the strand-winding motor,due to the higher slot filling rate.Ref.[8]used the hairpinwindings motor in the new generation bolt EVs to improve system efficiency,and found that increasing the number of layers of hairpin windings in each slot can effectively reduce its AC eddy current loss.Aiming at the problem of AC eddy loss of hairpin windings,Ref.[9]proposed an asymmetric structure to reduce the AC loss of the hairpin winding at high speed and make the winding current evenly distributed,thereby reducing the hotspot temperature and increasing continuous output capability.But when the motor is at low speed and the DC resistance of each layer of the winding is different,the performance of the motor will be reduced.

Ref.[10]studied the influence of the arrangement mode of hairpin windings in the slots on the current circulation problem.It found that if the conductors in each parallel branch are placed unevenly on each layer,it will lead to uneven distribution of current in two or more branches,causing circulation to increase AC loss and decrease the motor efficiency.To determine the effects of the number of parallel turns and the size of the hairpin winding on the AC loss of the motor at high speed,Ref.[11]conducted numerical analysis on the winding of a section hairpin-winding stator motor.Ref.[12]studied the heat dissipation performance of an oil-cooled hairpin-winding motor,and found that the use of oil injection cooling instead of water jacket cooling can greatly increase the power density and output torque of the hairpin-winding motor.

However,affected by the skin effect and proximity effect[13],the AC resistance of the hairpin winding will increase sharply at high speed,thereby affecting the output capacity and the motor efficiency in the high-speed range.Therefore,the low DC copper loss at low speed and the high AC copper loss at high speed will complicate the optimization design process of the hairpin-winding motor.How to optimize the comprehensive loss of the hairpinwinding motor under the driving cycle to improve the endurance mileage of EVs has become particularly important.

In view of the above problems,this paper systematically analyzes the AC loss characteristics of hairpin winding,aiming to propose an optimal design method of hairpin-winding motor under certain drive cycle.Firstly,the AC resistance analytical expressions of the hairpin windings at different speeds are deduced and an FEM model is established to calculate the AC resistance characteristics of the hairpin windings.Then,combining the modified particle swarm optimization(MPSO)algorithm with the FEM model,the optimization design method of the hairpin-winding PMSM under the driving cycle is proposed,and a 60 k W hairpin-winding drive motor is optimized.Meanwhile,the electromagnetic and the heat dissipation performance of optimized hairpin-winding motor is compared with the traditional strand-winding motor.Finally,a hairpin-winding prototype is manufactured and an experimental platform is built to test the electromagnetic and heat dissipation performance of the hairpin-winding motor,so as to verify the feasibility and accuracy of the optimized design method.

1 Structure and Theoretical Analysis

1.1 Hairpin-winding structure

The structures of traditional strand winding and hairpin winding are shown in Fig.1.Because of the reduction of air gap area and insulation area between conductors in the slots,the hairpin windings can effectively improve the slot filling rate and reduce the DC resistance and DC copper loss of stator windings.However,as the area of each layer of winding conductor increases,the eddy current loss of hairpin windings will increase greatly at high speed due to skin effect and proximity effect.

Fig.1 Conventional strand winding and hairpin winding

1.2 Analytical expression of AC resistance

Fig.2 shows the equivalent magnetic structure of hairpin windings.The width and height of stator slots arebsandhs,respectively.There arezlayers of hairpin conductors in each slot of the stator,and the width and height of each conductor arebcandhc,respectively.

Fig.2 Equivalent structure diagram of hairpin winding

The DC resistanceRDCof thekth layer hairpin conductor is

whereγis the winding’s conductivity.According to Ref.[14],the AC resistanceRACof thekth layer hairpin conductor is

whereKPkis the increasing coefficient of AC resistance influenced by skin effect and proximity effect,and its calculation expression is[14]

whereξis the relative height of each conductor,φ(ξ)the AC resistance coefficient caused by skin effect,andψ(ξ)the AC resistance coefficient caused by proximity effect.βis the reciprocal of the equivalent penetration depth and its expression is

wherefis the electrical frequency of the winding current andμthe winding permeability.

Since there are a total ofzlayers hairpin conductors in each slot,it is assumed that all conductors in each slot belong to the same phase,then the average resistance coefficientKP1of hairpin winding is

Similarly,if the conductors in each slot belong to two phases respectively,the average resistance coefficientKP21of the lower winding is[14]

The average resistance coefficientKP22of the upper winding is[14]

wherem=z/2,θ=±60°.

Based on the above expression,the AC resistance of the hairpin winding can be calculated theoretically at different speeds or electrical frequencies.

2 AC Characteristics of Hairpin Winding

In order to accurately analyze the AC characteristics of the hairpin winding,a 2-D FEM of a 60 kW hairpin-winding PMSM used for EVs is established,as shown in Fig.3.The main performance requirements and parameters of motor are shown in Table 1.

Fig.3 2D-FEM model of hairpin-winding motor

Table 1 Main requirements and parameters of 60 k W PM SM

2.1 Influence of electrical frequency on AC resistance

Different from the strand winding,the AC resistance of the hairpin winding will change greatly with speed.Fig.4 presents the intercorrelations(RAC/RDC)between AC/DC resistance ratio and frequency of hairpin windings at 20℃and 120℃.Fig.4 shows that the AC resistance of the hairpin winding increases with the increase of motor speed and frequency due to the skin effect and proximity effect.The AC resistance of the winding under the highest speed and 120℃is about 4 times of the DC resistance.

Fig.4 AC/DC resistance ratio vs frequency of hairpin winding

Moreover,because the winding resistivity is greatly affected by temperature,the higher the temperature,the lower resistivity.And according to Eq.(5),it can be seen that the lower the resistivity,the greater the skin depth.Therefore,the growth rate of the AC/DC resistance ratio with the speed at 120℃is slower than that at 20℃.

Meanwhile,the theoretical results have a certain deviation from FEM,which is mainly caused by the fact that the theoretical analysis does not consider the core saturation.

Another interesting phenomenon can be found through Eq.(4),that is,each layer of conductor in each slot is affected differently by proximity effect.Fig.5(a)shows the electrical densityJdistribution diagram of each layer conductor of the hairpin winding at different frequencies by FEM.It can be seen that as the frequency increases,each conductor of the hairpin winding has obvious eddy current effects.And the eddy current effect of the conductor close to the air gap is the most obvious.That is,the eddy loss of the conductor close to the air gap is large,and the current distribution is the most uneven.

Fig.5 Layered characteristics of hairpin winding

Figs.5(b,c)show the change trend of the AC/DC resistance ratio of each layer of hairpin conductor with frequency and the proportion of AC copper loss of each layer of conductor at different frequencies.Similar to the above conclusions,it can be found that the AC resistance of the conductor close to the air gap has the fastest growth,and the AC copper loss accounts for the highest proportion,accounting for about 58%of the total AC copper loss at the highest speed.The AC resistance of the conductor near the stator yoke increases slowly,and only accounts for about 9%of the total AC copper loss at the highest speed.

This layered characteristic of the hairpin winding makes the AC copper loss distribution of each layer of conductor extremely uneven.It is not conducive to the uniform heat dissipation of the stator windings,which is prone to hot spots,and at the same time increases the risk of excessive temperature rise of air gap.

2.2 Influence of winding temperature on AC resistance

Fig.6 shows the change trend of the AC resistance of the hairpin winding with frequency at different winding temperatures.What stands out in Fig.6 is that the DC resistance and the AC resistance at low frequencies of the hairpin winding increase with the increase of temperature.But as the frequency increases,the tendency of AC resistance to increase with the increasing temperature slows down.The most surprising aspect of the data is in the high frequency(>600 Hz),the AC resistance of the hairpin winding decreases with the increase of winding temperature,which is completely contrary to the change trend of DC resistance.

Fig.6 Temperature characteristic of hairpin winding

This result may be explained by the fact that the skin depth of the hairpin windings increases with the increase of temperature and the decrease of electrical conductivity.This phenomenon is beneficial to improving the heat dissipation performance of hairpin-winding motor at high speed and high temperature.

2.3 Influence of winding size on AC resistance

Obviously,the DC resistance and AC resistance of the hairpin winding will be affected by the conductor sizes(height,width).Among them,the larger the cross-sectional area of the winding,the lower the DC resistance.However,AC resistance is not only affected by the winding size,but also by the electric frequency.Fig.7(a)shows the changing trend of AC resistance of the hairpin winding with the conductor widthbcand the widthhcat 6 000 r/min and 12 000 r/min(120℃).The change trend of AC resistance is different at different speeds.At 6 000 r/min,the AC resistance decreases first and then increases with the increase of the conductor widthbcand the heighthc.But at the high speed of 12 000 r/min,the minimum value of AC resistance appears when the conductor widthbcand the heighthcare the smallest,which is completely opposite to the change trend of the winding DC resistance.

In addition,the change of the winding size will inevitably cause the change of the stator slot size,which will change the torque of the motor.Fig.7(b)shows the change of the output torque of the hairpinwinding motor with the conductor widthbcand the heighthc.The trend shows that the smaller the conductor widthbcand the heighthc,the higher the motor output torque.

Fig.7 AC loss and torque vs size of hairpin winding

Then,how to choose the optimal hairpin winding size has become a complicated problem for multiple working conditions.It is not only necessary to consider the selection of DC resistance in the lowspeed range,but also the change in AC resistance in the high-speed range to maximize the overall efficiency of the drive motor.

2.4 AC resistance of the end-winding

Since there is no magnetic circuit constraint around the end winding,the AC resistance of the end winding changing with the frequency is different from the effective winding.A 3D FEM of the hairpin-winding PMSM is established to calculate the AC resistance of the end winding with different frequencies,as shown in Fig.8.From the graph we can see that because the end winding is less affected by the proximity effect,its AC resistance increases slowly,which is much lower than the effective winding.As a result,the growth trend of the entire winding’s AC resistance decreases.

Fig.8 AC resistance of end winding vs frequency

3 Optimization Design of Hairpin-Winding PMSM

3.1 Number of winding layers

Before optimizing the hairpin-winding PMSM,the number of conductor layers per slot of the hairpin winding must be selected.Fig.9 presents the structure,slot filling rate and AC resistance of hairpin winding when the motor takes different conductor layersNs.It can be seen from Fig.9 that as the number of conductor layers increases,the conductor insulation area in the slot increases,resulting in a decrease in the slot filling rate and a corresponding increase in the DC resistance.But at the same time,the structure of each conductor tends to be flatter,so its AC resistance decreases with the increase of the number of conductor layers at high speed.In order to reduce the DC resistance of the motor and reduce the manufacturing cost,the number of conductor layers of the hairpin-winding motor is selected as four layers in this paper.

Fig.9 Characteristics of hairpin with different conductor layers

3.2 Optimization process

Table 1 shows the constraints and design objectives of the PMSM designed in this paper.It has been previously observed that the main objectives of this optimization are to minimize the comprehensive loss of the motor under certain driving cycle,and to maximize the power density or torque density.

The working condition scatter distribution of the driving motor designed in this paper under the Urban Dynamometer Driving Schedule(UDDS)driving cycle is shown in Fig.10.From the scatter,it is apparent that the motor mainly works in lowspeed(climbing,starting)overload peak operating conditions,long-term medium-speed stable rated operating conditions(5 000 r/min@60 k W)and short-term high-speed operating conditions.Therefore,it is necessary to comprehensively consider the losses in these operating conditions in the motor design process to improve cruising range of EVs.

Fig.10 Working scatter of motor by UDDS cycle condition

However,in the actual optimization design process,it is unrealistic to calculate the loss of all cycle operating points to obtain the comprehensive loss of the motor.For the purpose of rapid optimization,the center points of the three operating conditions are defined as the rated operating point,the peak operating point and the highest speed operating point,and the equivalent comprehensive loss of the motor under UDDS is defined as

Equivalent Loss=wrPrated+wpPpeak+whPhigh(9)wherePrated,PpeakandPhighare the total losses of the drive motor at rated operating point,peak operating point and maximum speed point,respectively,including iron loss,AC copper loss,mechanical loss and additional loss.wr,wp,andwhare the weight proportions of the motor under the rated condition,peak condition and maximum speed condition with the UDDS driving cycle,respectively.wr,wp,andwhare obtained by accumulating the power of all operating points in the three regions in Fig.10 and by calculating the proportions.

Fig.11 shows the flow chart of the optimized design method for the hairpin-winding motor proposed in this paper,which mainly includes three steps:

Fig.11 Optimization flow chat

Step 1Carrying out the preliminary electromagnetic theoretical design of hairpin-winding motor.

First,analyze the desige requirements and design objectives,use the theoretical formulas and design experience to calculate the initial size of the motor,and then use sensitivity analysis to define the optimization variables.The optimization objectives in this paper is to increase the output torque of the motor and reduce the comprehensive loss of the motor.The optimization variables are the inner diameter of the motor stator,the width and the height of the hairpin-winding conductor.

Step 2Optimizing the torque and comprehensive loss of the motor.

At this step,MPSO[15]is used to initialize and update the particle swarm population parameters.In the iterative process,the FEM is used to calculate the iron loss,AC copper loss and peak output torque under different parameters,the weight coefficient is used to calculate the comprehensive loss of the motor under the three working conditions and different size parameters,and MPSO is used for iterative optimization until the design scheme converges.

Step 3Verification of temperature rise and structural strength.

When the electromagnetic performance of the optimized motor meets the requirements and reaches the optimum,the temperature rise and structural strength of the designed motor are verified.

3.3 Optimization r esults and compar ison

Through the method,the performance of the hairpin-windings motor is optimized for multiple iterations.Its comprehensive loss and output torque are optimized,as shown in Fig.12.

It is obvious from Fig.12 that as the optimization progresses,the population gradually approaches to the direction of high torque output and low comprehensive loss,and finally converges to a certain position in space,forming the optimal solution set of the Pareto front.

Fig.12 Pareto front and optimal solutions of optimization

Table 2 shows the parameter comparison of the hairpin-winding motor before and after optimization.It can be seen that although the winding DC resistance of the hairpin-winding motor increases after optimization,its comprehensive loss decreases from 5 699.7 W to 4 879.6 W and the peak torque is also improved from 286.87 N·m to 321.85 N·m,which verifies the effectiveness of the proposed optimization design method.

Table 2 Par ameter comparison of hair pin-winding motor before and after optimization

4 Comparison of Hairpin-Winding Motor with Strand-Winding Motor

In order to verify the superiority of the performance of the hairpin-winding motor in the field of EVs,this paper also optimizes a strand-winding PMSM under the same restrictive conditions according to the performance requirements of Table 1.The FEM models of the two winding type motors are shown in Fig.13.

Fig.13 Models of hairpin-winding and strand-winding motors

Here,the number of conductors per slot of the motor is determined by the current limit,the required output torque and the number of stator slots.Therefore,the number of conductors per slot of the two winding types of motors is the same.The difference is that the hairpin winding has only one conductor per turn.However,when the number of conductors is determined for the strand-winding motor,the appropriate wire diameter needs to be selected,and each turn of the coil is composed of multiple strand wires.

4.1 Electromagnetic performance comparison

Table 3 shows the comparison of the specific performance parameters of the two winding types of motors,and Fig.14 shows the comparison curves of the AC resistance,and the torque and power external trajectory of the two winding types of motors.

Fig.14 Comparison of electromagnetic performance

It is apparent from Table 3 that under the same restriction conditions,the hairpin-winding motor has a higher slot full rate and a lower DC resistance value than the strand-winding motor.Therefore,the AC resistance of the hairpin-winding motor in the low speed range(<6 000 r/min)is smaller.However,when the speed exceeds 6 000 r/min,the AC resistance of the hairpin-winding motor will exceed the strand-winding motor.

Table 3 Comparison of electromagnetic and heat dissipation

In addition,due to the decrease of the core saturation of the stator teeth,the output torque of the hairpin-winding motor increases and the active length decreases.The reduction of the axial length reduces the no-load back EMF value.Therefore,the flux weakening current of the hairpin-winding motor is small and the peak power is high in the high-speed range.The peak power of the strandwinding motor is 139.2 kW,while the peak power of the hairpin-winding motor is 148.4 kW,which is an increase of 6.6%.In addition,the active length of the hairpin-winding motor is reduced by 10.6%,so the effective power density of the hairpin-winding motor is increased by about 19.2%compared with the strand motor.

Fig.15 shows the efficiency map of the two winding types of motors.We can see that efficiency of hairpin-winding motor at low and medium speed range is superior to the traditional strand-winding motor due to the low copper loss and iron loss.However,the efficiency of the hairpin-winding motor is lower than that of the strand-winding motor at high speed range due to the extremely increased AC copper loss.Generally,the high efficiency range of hairpin-winding motor is significantly higher than that of traditional strand-winding motor,which can make the driving system run more efficiently.

Fig.15 Comparison of efficiency map at full speed range

4.2 Thermal per formance compar ison

In order to verify the heat dissipation performance of the hairpin-winding motor,this paper uses Fluent software to establish the fluid-temperature field simulation model and calculates the steadystate temperature rise of the two winding types of motors.

Table 4 presents all the material properties of motor parts.It is worth noting that due to the complex structure of the strand windings,it is difficult to establish an accurate model of winding in the fluid-temperature field.Therefore,the equivalent method given in Ref.[16]is used to make all strand coils,air gaps and film in each slot of strand winding equivalent to a conductor,and the end of the winding is equivalent to an end ring.

The cooling method of the motor studied in this paper is housing water jacket cooling,and the water channel is the axial Z type.Fig.16 shows the meshing of the fluid-temperature coupling model and the water channel structure.In the simulation,the inlet of the water channel adopts the boundary condition of the flow inlet,the inlet flow is 10 L/min,and the coolant temperature is 60℃.The outlet adopts the pressure outlet condition,as shown in Table 5.

Table 4 Material properties of motor parts

Fig.16 Meshing and cooling channel structure of motor

Table 5 Simulation boundar y condition

Firstly,the two winding types of motors can obtain the same temperature rise(80℃)by changing the copper loss(other losses are the same),as shown in Fig.17.At this time,the AC copper loss of the hairpin-winding motor is 1 384 W,while the AC copper loss of the strand-winding motor is 932 W.Therefore,the hairpin-winding motor can output higher power under the same steady-state temperature rise limit,which further improves the power density of the motor.

Fig.17 Same steady temperature rises under different losses

Similarly,the temperature rise distribution of the two motors at the same loss is shown in Fig.18.

Fig.18 Steady temperature rises under the same losses

It can be known that under the same loss and coolant temperature(60℃),the steady-state highest temperature rises of the hairpin-winding motor and the strand-winding motor are 124℃and 140℃(all appear in the windings),respectively.There is a significant difference between the two winding types of motors,which verifies that the hairpinwinding motor has better heat dissipation performance.

5 Experiment

In order to verify the accuracy of the optimization design method and understand the actual working performance of the hairpin-winding PMSM,a hairpin-winding prototype is manufactured according to the aforementioned design scheme,as shown in Fig.19.

Fig.19 Hairpin winding and prototype

Fig.20 shows the simulated and measured power,torque external trajectory,efficiency and AC copper loss comparison curves.In addition,the AC copper loss is calculated by testing the no-load mechanical loss and no-load iron loss,and then through the inverse calculation of the input power and efficiency under load conditions.

Fig.20 Comparison of simulation and tested characteristic curve

We can see from Fig.20 that the simulated torque and power are basically consistent with the measured results at low speed range.However,when the motor is running at high speed,the measured output power is obviously lower than the simulated one,which is mainly caused by the decrease in voltage utilization rate during high-speed testing.The simulated efficiency and AC copper loss are lower than the actual measured results at high speed,which is mainly due to the fact that the highfrequency carrier introduced by the controller is not considered during the simulation.

Fig.21 shows the comparison of simulated and measured steady temperature of the hairpin-winding motor under rated conditions at different positions.The simulation error is basically kept within 10%,which verifies the accuracy of the temperature field simulation in this paper.

Fig.21 Comparison of temperature of steady condition

6 Conclusions

This paper firstly analyzed the AC characteristics of the hairpin-winding motor through the analytical model and the FEM.Then,combining FEM and the MPSO algorithm,a 60 k W hairpin-winding PMSM is optimized under the UDDS driving cycle.And the optimized results are compared with traditional strand-winding motor in electromagnetic and heat dissipation performance.Finally,the prototype test is used to verified the accuracy of the optimized design method.The main conclusions are as follows:

(1)Affected by the skin effect and proximity effect,the AC resistance of the hairpin winding increases with the increase of the motor speed.And the AC resistance of the effective part of the winding reaches 4—6 times of the DC resistance at the highest speed.Meanwhile,the AC resistance of the hairpin winding will decrease with the increase of temperature at high speed.

(2)The increasing tendency of AC resistance at the hairpin end winding with the increase of speed is slower than that of the effective part,mainly because the end winding is less affected by the proximity effect.

(3)Due to the influence of the eddy effect,the AC resistances of hairpin winding at low speed,medium speed and high speed are completely different with the change of the winding size,resulting in the optimization of hairpin-winding motor becoming a complex multi-condition problem.

(4)Compared with traditional strand-winding motor,the high efficiency range of hairpin-winding motors in the low and medium speed ranges is obviously better than that of traditional strand motors,but the efficiency at high-speed range is lower due to the larger AC copper loss.In addition,the hairpinwinding motor also has higher peak output power and power density and better heat dissipation performance compared with the strand-winding motor,which is suitable for application in the field of EVs.

主站蜘蛛池模板: 亚洲人成网18禁| 无码精品国产dvd在线观看9久 | 亚洲精品久综合蜜| 啪啪永久免费av| 亚洲男人天堂2020| 四虎永久在线| 狠狠综合久久| 国内精品一区二区在线观看| 欲色天天综合网| 久久香蕉国产线看观看式| 欧美天堂久久| 精品久久久久成人码免费动漫| 欧美国产日韩在线观看| 国产97公开成人免费视频| 国产精品流白浆在线观看| 国产精品密蕾丝视频| 国产h视频免费观看| 欧美人与牲动交a欧美精品 | 色婷婷亚洲综合五月| 91在线精品麻豆欧美在线| 亚洲婷婷丁香| 国产综合另类小说色区色噜噜| 麻豆精品在线播放| 久久99国产精品成人欧美| 国产性生交xxxxx免费| 少妇精品网站| 欧美激情视频在线观看一区| 波多野结衣无码中文字幕在线观看一区二区 | 免费人成网站在线观看欧美| 国产麻豆91网在线看| 又黄又爽视频好爽视频| 国产欧美日韩资源在线观看| 国产精品妖精视频| 日韩欧美中文| 伊人久久久久久久| 国产99久久亚洲综合精品西瓜tv| 色综合成人| 欧美在线黄| 人妻一区二区三区无码精品一区| 精品1区2区3区| 欧美日本二区| 人妖无码第一页| 亚洲人成在线精品| 91精品亚洲| 久久伊人操| 国产欧美视频一区二区三区| 色综合久久88| 国产精品一区二区不卡的视频| 国产精品成人免费综合| 久久精品人妻中文系列| 亚洲中久无码永久在线观看软件| 美女国内精品自产拍在线播放 | 一区二区三区毛片无码| 国产欧美精品午夜在线播放| 国产网站一区二区三区| 亚洲永久色| 国产菊爆视频在线观看| 久久综合色视频| 亚洲国产精品不卡在线| 亚洲人在线| 国产无码制服丝袜| 国产精品免费露脸视频| 色亚洲成人| 亚洲天堂日本| 亚洲浓毛av| 国模沟沟一区二区三区| 中文字幕佐山爱一区二区免费| 亚洲欧美不卡中文字幕| 免费av一区二区三区在线| 一级毛片免费的| 四虎成人精品| 亚洲国产综合精品中文第一| 国产欧美亚洲精品第3页在线| 精品91在线| 手机在线国产精品| 国产乱子伦无码精品小说| 日本午夜网站| 国产精品99在线观看| 日韩a在线观看免费观看| 另类综合视频| 日本免费精品| 久热re国产手机在线观看|