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

    Design optimization for a launching system with novel structure

    2019-11-18 02:34:20RunduoCaoXiaobingZhang
    Defence Technology 2019年5期

    Run-duo Cao, Xiao-bing Zhang

    School of Energy and Power Engineering, Nanjing University of Science and Technology, Nanjing, 210094, Jiangsu, China

    Keywords:Interior ballistics Mobile chamber Design optimization

    ABSTRACT A launching system with novel structure using optimization method is investigated to improve the muzzle velocity of guns in this article. This system has two tandem chambers of which the front one is ignited first and the other is ignited after a while.The launching process of this novel system is modelled and a series of different schemes are simulated,to discover the interior ballistic performance of this novel launching system. In order to obtain the optimal loading conditions, an optimization model combined with the combustion model is established. The optimal schemes can improve the muzzle velocity by 20.6% without changing the parameters of barrel. It means that this novel launch system could improve the interior ballistics performance significantly and it still has considerable potential to be ameliorated.?2019 The Authors.Production and hosting by Elsevier B.V.on behalf of China Ordnance Society.This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

    1. Introduction

    Since guns were invented, nearly all guns applied a launching principle in which the projectile is pushed by high-pressure gas generated from burning propellants in a chamber. However, this launching principle has an inherent deficiency that, with the motion of the projectile in chamber,there is a pressure gradient along the gun barrel because the expansion of gas is lower than the projectile.Thus,the pressure behind the projectile is always lower than the pressure at the breech. For the interior ballistics of guns,the muzzle velocity of the projectile depends on the projectile base pressure,but the breech pressure is limited by the barrel design and materials,so that the guns could not get a higher muzzle velocity in the conventional launching method. So far, many basic methods such as enhancing the impetus of propellant,increasing the mass of propellant, extending the length of gun barrel, have been investigated but fail to solve this problem.

    In order to improve the muzzle velocity of guns, designers worldwide have invented and developed many novel methods by adopting some novel technologies.For example,a kind of traveling charge technique to maintain higher pressure have been tested for decades[1e3],in this method,the propellant is combusting behind the projectile constantly, which pushes the projectile more efficiently. And some electrothermal launchers are also proposed by some researchers [4e6], in which the propellant in the chamber reacts with the plasma generated by the pulse power and it releases much energy to give the projectile a high velocity. Comparatively,the electromagnetic launch is developed more deeply, and some prototypes have been developed in America by Dr. Fair [7e9] and some other countries [10]. This method can provide a sustaining acceleration during the complete interior ballistics process so that the projectile could get a high muzzle velocity. In addition, some other launching systems were also proposed, such as ram accelerator [11], rarefaction wave gun [12,13] and so on. These launching methods all applied advanced theory to improve the muzzle velocity, however, these methods also have disadvantages which are not able to be solved at present. For example, the traveling charge system requires liquid propellants or solid propellant with high burning rate, whose combustion is difficult to control. As for the electromagnetic launcher, it has a high requirement for electrical power system so it is not convenient to bring to service in the short term. In addition, the electrothermal launchers and rarefaction wave gun have been developed for many years but few achievements have been made.

    In this paper a mobile chamber launching system (MCLS) is proposed, of which the combustion chamber is divide into two chambers, separated by a piston, with one chamber fixed and the other moving. Although the similar structure has been proposed[14],the ignition order and motion of projectile for this method in this paper are both different.For this method,the propellants in the mobile chamber are ignited first and then generate gas to drive the projectile. After several milliseconds, the propellants in the fixed chamber are ignited and the piston starts moving under the pressure in both chambers. In the wake of the piston moving, the mobile chamber expands slowly relatively so that there are two pressure peaks in the mobile chamber and the pressure at the base of projectile drops more slowly than that in a conventional gun.Thus,the propellent gas does more work on the projectile along the barrel so the projectile can attain a higher muzzle velocity. Moreover, this launching system could be applied into existing guns because it avoids too much modification on the guns. Considering the difference in launching process between proposed method and conventional launching method, the modified physical and mathematical models are established in section 2. Because this novel launching system is an original structure for guns, the design experience for the conventional guns is for reference only but not able to obtain the best design scheme. Thus, the optimization design combines the interior ballistic model and optimization model is necessary. Some discussions about selecting variable and the design optimization applying a standard particle swarm optimization(SPSO)algorithm is described in section 3.The results and discussion about the launching process are expressed in section 4.This paper is for the purpose of providing a new approach to improve the muzzle velocity on the conventional guns. The simulation results indicate that this method is effect.

    2. Model establishment and simulations

    The schematic of the mobile chamber launching system is shown in Fig.1, which includes three main parts: a fixed chamber including the grain propellant and igniter tube, a mobile chamber consisting of a piston, stick propellant and grain propellant, and a projectile. For the proposed system, the igniter tube connects between the breech and the piston, which can avoid the piston moving earlier.As similar with the conventional launching system,a lumped parameter model proposed by Zhang [15], which is a useful and robust model for the interior ballistic cycle, is used to simulate this novel launching process.In this model,propellants in each chamber are assumed to be ignited simultaneously and uniformly, while their combustion is assumed to take place in a smoothly varying, well-stirred mixture, and the burning rate is determined by the instantaneous, space-mean chamber pressure.The state equation for gunpowder gas is Nobel-Abel equation given in Ref.15. The engraving process of projectile is ignored. In the process of combustion, the impetus of the propellant, the specific ratio and the covolume are constant.

    2.1. The model of launching process

    The launching process of this novel method can be divided into 4 stages and each stage is analyzed below. The schematic of launching process is shown in Fig. 3 (a)~(f). In the mathematical model of launching process, subscript “1” and “2” represent the parameters of fixed chamber and mobile chamber respectively,and subscript “c” and “p” represent the parameters of piston and projectile respectively.

    1. Combustion model of propellants

    For the propellants in the two chambers, a common parameter model is applied,which is described in Eq.(1).In this article,a 19-perforation rosette propellant is used, whose cross-section schematic is shown in Fig. 2.

    Fig.1. Schematic of the mobile chamber launching system.

    Fig. 2. 19-perforation rosette propellant.

    2. The first stage of launching process

    The first stage,which is shown in Fig.3(a),begins with ignition of the propellant in the mobile chamber and ends when the propellant in the fixed chamber is ignited. In this stage, the projectile moves ahead and other parts are immobilized. The state equation in mobile chamber is shown as Eq.(2)and the energy conservation equation in mobile chamber can be described as Eq. (3).

    According to Eqs. (2) and (3), the pressure in the mobile chamber can be obtained by Eq. (4). Besides, the motion of piston and projectile are described as Eq. (5) and pressure in fixed chamber remains unchanged in this stage as Eq. (6).

    2 The second stage

    Several milliseconds after the projectile starts moving,which is called delay time(td),this stage starts when the propellants in fixed chamber begin burning and ends when the piston begins to move,which is shown in Fig.3(b).In this stage,the pressure in the fixed chamber is less than that in the mobile chamber so that it is not able to push the piston forward. Because the projectile is still moving individually while the other parts are kept still, the expression of pressure in the mobile chamber and the motion of projectile are both same as previous as Eq.(4)and Eq.(5).However,the pressure in the fixed chamber changes as Eq. (3), which is similar to the combustion of propellants in a closed bomb vessel.

    3 The third stage

    Fig. 3. Schematic of the launching process.

    When the piston starts moving, which means the pressure in the fixed chamber is large enough to push the piston forward and the motion of piston can be described as Eq. (8). In this stage, the projectile is driven only by the pressure in the mobile chamber and its motion is described as Eq. (9).

    Meanwhile,based on the law of conversation of energy and state equation of gas,the state in fixed chamber and mobile chamber can be described as Eq. (10) and Eq. (11).

    Thus,the pressure in the fixed chamber and the mobile chamber can be obtained by Eq. (12) and Eq. (13).

    This stage finishes when the projectile leaves the gun barrel as shown in Fig.3(d).So far,the main process of launching is over and the muzzle velocity can be obtained.

    4 The last stage

    After the projectile has left the gun barrel, the main process of the interior ballistic cycle has been finished.However,because the piston is still moving, this stage ends when the piston leaves the gun barrel. In this stage, the pressures in the two chambers drive the piston and the outflow of gas in the mobile chamber is assumed to be adiabatic. It could be seen as another projectile launching from the gun. Therefore, the pressure in both chambers could be calculated as Eq.(14)and Eq.(15),which are given by Ref.[16],the motion of piston is described as Eq. (16) and Eq. (17).

    where,p2;outand toutrepresent the pressure in the mobile chamber and the time at the moment of the projectile leaving the gun barrel.

    Thus, the interior ballistic process of this method can be simulated by Eq. (1) to Eq. (17).

    Considering the pressure distribution in two chambers, the motion equation of piston (9) is calculated by the pressure at the piston base (pb1) and the pressure of the projectile base (pb2). In order to obtain them, we assume that the propellant gas in both chambers obey the Lagrange hypothesis[15].Because of the muzzle velocity is over 600 m/s, for which the classic The Lagrange pressure gradient model is not reasonable [17,18], a modified Lagrange pressure gradient model given by Ref.[19]is adopted.This model is modified by a series of experiment and the results show that it applied to for some guns of which the muzzle velocity is excess of 1700 m/s. Thus, this model is suit for the proposed method in this paper.The detailed process of this method is presented as follows:

    2.2. Numerical simulations

    After establishing the launch model, some numerical simulations are proceeded to investigate if this novel method could improve the interior ballistic performance.The original parameters of a 130-mm cannon and propellant are listed in Table 1.And some additional parameters for mobile chamber are listed in Table 2,which are inspired by the previous work and the conventional launching system.

    Table 1 The original parameters of a 130-mm cannon and propellant.

    Table 2 Additional parameters of propellant for MCLS.

    The simulation results are shown in Figs.4e7,Figs.4 and 5 show the pressure history and velocity history of different methods.Due to the limit of the pressure in fixed chamber,the pressure in mobile chamber cannot attain a high level, and vpis lower than v0before the piston starting moving. However, the piston starts moving while the p1exceeds p2so that the volume of mobile chamber expands slowly and p2can hold a high level and drop more slowly than p0. Figs. 6 and 7 is the pressure and velocity to the projectile displacement of different methods.Fig.6 shows that the maximum of p0appears while the piston moves about 1 m,and the maximum of p1appears at the breech of gun. The simulation results indicate that two methods can obtain similar muzzle velocity, even if p0is higher than p2in the initial stage.However, the parameters of the mobile chamber in this condition is just obtained by experience,the pressure in mobile chamber is too low to do enough work on the projectile. Thus,it is necessary to obtain the best scheme by using the design optimization to achieve full potential of this launching system adequately.

    3. The design optimization

    According to the model established above, the shape of the mobile chamber should be designed carefully to meet the requirements. For the mobile chamber, the most important parameters are its volume and charge mass,which is similar to the fixed chamber.Besides,the mass of piston and the delay time,which are new variables involved in interior ballistics process, also influence the performance. Therefore, some schemes are simulated to investigate the relationship between these parameters and interior ballistic preference. These simulations also provide us some experience on choosing value ranges of each variables.

    3.1. Analyses of some variables

    Fig. 4. Pressure history of different methods.

    Fig. 5. Velocity history of different methods.

    Fig. 6. Pressure to the projectile displacement of different methods.

    Fig. 7. Velocity to the projectile displacement of different methods.

    At first, this article discusses about effects of the piston mass(mh), delay time (td) and the propellants thickness (2e2) in the mobile chamber on muzzle velocity, and the results are shown in Figs. 8e10. It is important to note that while one variable is considered, some necessary changes of the propellants' shape are taken to ensure an equal maximum pressure in the mobile chamber every time, but the mass of charge is unchanged. In Fig. 8, muzzle velocity (vg) declines and the maximum pressure rises with the mass of piston increasing,because the lighter the piston is,it can be accelerated more easily, so that the pressure in mobile chamber could be maintained to work on the projectile more efficiently.However, with the mcincreasing, the piston needs a higher pressure to be push forward and more energy of propellant is converted to the velocity of the piston, so that the velocity of projectile declines but the maximum pressure in fixed chamber(P1max)rises.As for the influence of the delay time, Fig. 9 shows that both vgand P1maxdecline with the tdbeing extended. If the delay time were shorter, pressure peaks in two chambers would appear at the similar time so the piston was difficult to move to enlarge the fixed chamber and the P1maxwas at a very high level. If the delay time were too long, the propellant in the fixed chamber does not have enough time to burn out, so the projectile fails to obtain a high velocity. In Fig.10, vgand P1maxboth decline with 2e2increasing,because with thinner 2e2, the propellants generate gas more quickly, leading to larger P1max.

    Fig. 8. vg and P1max with changes of mc/mp.

    Fig. 9. vg and P1max with changes of td.

    Fig.10. vg and P1max with changes of 2e2.

    Fig.11. Pressure history with different mh/mp.

    Fig.12. Pressure history with different td.

    Fig.13. Pressure history with different 2e2.

    Moreover, in order to investigate the influences of these variables on pressure, more schemes are simulated to observe the changing process of pressure in the two chambers,and the results are shown in Figs.11e13. It is worthy to note that in each simulation, only the corresponding variable is changed and other parameters of the charge are unchanged. Fig. 11 shows that the maximum pressure in the fixed chamber rises along with increasing mh/mp, but the second pressure peak in the mobile chamber changes on the contrary.Because with larger mh,it needs a higher p2to give piston sufficient acceleration to shrink the mobile chamber relatively. As for p2, because before piston starts moving,the volume of mobile chamber does not change,so it does not change obviously.Similarly,Fig.12 shows the influence of tdon p1and p2,it is obvious that when the tdis short,the pressure peak in the fixed chamber appears early so that the piston moves early and p1declines slowly. When tdis too long, the piston starts moving later to shrink the mobile chamber and p1is nearly not affected.As for 2e2, p2rises more quickly while 2e2is thinner, which means it needs a larger p1to drive the piston so that pressures in both chambers are affected by 2e2,as it is shown in Fig.13.Considering the safety during the launching process, a too thin thickness of propellants in the mobile chamber should be avoided.

    The previous simulations and analyses have shown the significant influence on the interior ballistic performance due to the mass of piston, the delay time and thickness of extra propellants.Moreover, different combinations of them will lead to more variations. Because there is no similar design before, and this novel launching system involves many multidisciplinary variables of charges and canon structure in both chambers.It is hard to obtain a superior design scheme based on experience. Therefore, using an optimization method is an ideal way to predict the overall performance of an interior ballistic design. Our group has successfully applied some optimization algorithms into the charge design and other complex design optimizations [20,21]. Thus, similar approaches are used to solve the design optimization of the proposed launching system in this article as follows.

    3.2. Optimization mathematical model

    Generally, the optimization problem involves objective, design variables and contains.

    1) Design objective

    The mobile launching system is applied to improve the muzzle velocity,so the only design objective is muzzle velocity of projectile(vg).

    2) Design variables

    To satisfy the design requirements, it is important to choose some appropriate design variables.Considering the above analyses on the proposed method, it is known that if the structural parameters of gun are pre-determined, only the volume of mobile chamber(Vm),mass of piston(mc),and delay time of the ignition in fixed chamber (td) can be modified. In addition,the main characteristic parameters of the propellants in both chamber, which include web thickness(2e1(2)),perforation diameter(d1(2)),and the charge mass (u1(2)), can be redesigned. Therefore, all these nine parameters should be regarded as design variables in the design optimization. The structural parameters of guns and some design variables are shown in Fig.14.

    3) Constraints

    Considering the particularity of the proposed launching system,in order to insure the launching safety and tactical requirements,some special and necessary constraints are listed as follow, which include maximum pressure and charge density.

    At first, because of the limit of the gun barrel's strength, the pressure in both chambers(Pm1(2))must be lower than a maximum pressure(Pmax),of which the pressure distribution along the barrel have been considered. Furthermore, the charge density (D1e2T)which is calculated by Eq.(19)should be limited as follows,because it has a great effect on the launching safety.

    According to the previous analyses, an optimization mathematical model applied in this article could be described as follows:

    To achieve the optimization of this novel launching system, an enhanced SPSO algorithm [22] is applied in this article. Using a 130 mm cannon as an example, the combustion model for multiperforated propellant and the lumped parameter interior ballistic model for mobile chamber are firstly combined to establish the optimization model, and then SPSO is applied to find the optimal result. The flowchart of optimization design is shown in Fig.15.

    4. Optimization results

    Here,the first aim is to obtain the optimal scheme of the design variables and the second aim is to validate whether the mobile launching system improves the interior ballistic performance significantly. The values of the original parameters of a 130-mm cannon and a 19-perforation rosette propellant have been listed in Table 1. As for the optimization model, according to the result and discussion above, the lower bounds (LB) and upper bounds(UB) of each variable are given in Table 3. The values of each variable are inspired by the simulation results above,which ensure the maximum pressure in chamber not exceeding the design limit.

    After 500 iterations,a converged result is obtained and listed in Table 3,and the convergence process is illustrated in Fig.16(a)-(c).It is obvious that after hundreds of iterations, all variables and objective remain unchanged and the muzzle velocity changes from 975 m/s to 1026 m/s, which means the optimization algorithm is useful. Moreover, each variable has an obvious change at the objective's every improvement, which means the interior ballistics performance of this novel launching system depends on these variables and our choices about these variables are appropriate.

    Then the parameterized model is applied to simulate the launch process of the optimum scheme and results are listed in Table 4.The muzzle velocity of this novel method can attain 1026 m/s, which has an improvement of 20.6%. And the maximum pressure in two chambers are in the safe range. These results indicate that our proposed mobile launching system can improve the muzzle velocity significantly.

    In order to explain why this novel launching system could improve the performance of the interior ballistics system, Fig. 17(a)-(b) show the pressure and projectile velocity variation over time and projectile displacement. Compared with the parameters obtained by experience in Tables 1 and 2, 2e2and u2increases evidently leading that the energy in mobile chamber is higher and the propellant burns slowly so that it releases gas in a long time.Besides, the mass of piston is much lighter so that it can be accelerated more easily.

    Fig.14. Parameters of cannon and projectile.

    Fig.15. The flowchart of optimization design.

    Table 3 The Lower Bound, Upper Bound and Optimum results of each variable.

    Fig. 17 (a)e(b) show that when the propellant in the mobile chamber begins burning,p2raises quickly due to the charge density is pretty high. Considering the delay time, the p2is higher than p1even pass the first pressure peak. However, after the ignition in fixed chamber, the p1climb more fast and exceed the p2, then the piston starts moving and obtains a large acceleration. Fig. 17 (b)shows that the piston moves even faster than piston in a short period. Thus, the volume of the mobile chamber to expand slowly and even shrink so that the pressure in the mobile chamber can achieve the second pressure peak, which is higher than first one.For the piston adjusting two chambers dynamically in the effect of the p1and p2,Fig.17(a)shows that after about 8.2 ms,the pressure in both chambers are equal,which means the velocity of the piston stays unchanged as shown in Fig.17(b) so that energy of the propellant in two chambers both work on the projectile. Because the pressure in the mobile chamber remains at a high level, the projectile can obtain energy effectively. As shown in Fig.17(c), the p2does more work on projectile than p0after the projectile moving 2 m, so that the projectile can obtain a higher muzzle velocity by the proposed method. Fig.17 (b) show that the v1increases more quickly than v0in the most of launching process,the projectile can maintain a large acceleration in the whole time so that the proposed method can obtain a higher muzzle velocity. These results indicate that without increase the maximum pressure in chamber,the proposed method can improve the muzzle velocity obviously.

    Fig.16. The convergence process of each variable.

    Table 4 The original parameters of a 125 mm gun and propellant.

    Fig.17. The launch process of different launch method.

    5. Discussion

    Although our optimal design scheme could improve the muzzle velocity from 851 m/s to 1026 m/s,this novel launching system still has potential for perfection. For example, the simulation above is based on an existing cannon, whose fixed chamber and barrel are unchangeable, the propellants in two chambers are the same in shape and composition. If the fixed chamber or barrel could be redesigned and different propellant are applied considering the different requirement of two chambers, this novel launching system might achieve a much higher muzzle velocity. Besides, the mass of the mobile chamber also influences the acceleration in the initial state, and it may decrease by applying some new materials.

    For this proposed method, some methods are adopted to deal with new problems. Before the piston starting moving, a highstrength igniter tube is applied to prevent the piston moving earlier. Besides, the stick propellant and multi-point ignition are used to ensure the ignition in mobile chamber uniform as far as possible. In addition, there are also some problems existing with this system.For example,there are two pressure peaks in the barrel and their position is after that in the conventional method and some additional measures should be applied to ensure the launching safety.On the other hand,owing to the shorter effective movement and much higher velocity, the movement time is short so the proposed method requires a high-burning rate propellant.

    6. Conclusions

    In this article, a novel mobile chamber launching system is proposed, with its interior ballistics model being established and the launching process being simulated.Compared with other novel launching methods,such as electromagnetic gun,rarefaction wave gun and so on, the proposed method can be applied into the existing guns without too many modifications for guns, so it is a low-cost system and suitable for being equipped in a short time.

    Then a design optimization for the novel launching system using the SPSO is applied to obtain an optimum scheme. In order to evaluate its performance, it is compared with an original scheme and an optimum scheme in conventional method.The results show that the novel method could improve the muzzle velocity by 20.6%,while the maximum pressure in chamber has no significant change.Owing to much potential for the mobile launching system,it could be a very worthwhile technology to improve the interior ballistic performance. Although the result is satisfactory, this paper only provides this new approach and adopts a simplified interior ballistic model to simulate the launching process. For the purpose of describing launching process in more detail,some complex model can be applied in the future work.

    美女高潮的动态| 久久久久久久久久黄片| 国产精品美女特级片免费视频播放器| 国产男人的电影天堂91| 不卡一级毛片| 亚洲人成网站在线观看播放| 午夜激情福利司机影院| 我要看日韩黄色一级片| 久久人人精品亚洲av| 国产精品国产高清国产av| 别揉我奶头 嗯啊视频| 在线观看午夜福利视频| 波多野结衣高清作品| 久久久久久九九精品二区国产| 波多野结衣巨乳人妻| 91精品国产九色| 少妇裸体淫交视频免费看高清| 精品国产三级普通话版| 高清午夜精品一区二区三区 | 亚洲av一区综合| 久久99热这里只有精品18| 1000部很黄的大片| 欧美激情久久久久久爽电影| 午夜视频国产福利| 99热这里只有是精品在线观看| 嫩草影院入口| 日韩中字成人| 午夜福利在线在线| 欧美区成人在线视频| 欧美日韩一区二区视频在线观看视频在线 | 日本a在线网址| 久久久久久大精品| 亚洲精品日韩av片在线观看| 国产视频一区二区在线看| 日韩在线高清观看一区二区三区| 亚洲欧美精品综合久久99| 少妇人妻精品综合一区二区 | 亚洲乱码一区二区免费版| 悠悠久久av| 99久久九九国产精品国产免费| 99精品在免费线老司机午夜| 亚洲第一电影网av| 久久综合国产亚洲精品| 久久久成人免费电影| 精品一区二区三区人妻视频| 搡老熟女国产l中国老女人| 精品99又大又爽又粗少妇毛片| 免费观看的影片在线观看| 久久综合国产亚洲精品| 我的女老师完整版在线观看| 中文字幕av成人在线电影| 国产毛片a区久久久久| 国产色婷婷99| 国产精品伦人一区二区| 亚洲精品国产av成人精品 | 亚洲不卡免费看| 乱系列少妇在线播放| 狠狠狠狠99中文字幕| 波多野结衣高清作品| 五月伊人婷婷丁香| 不卡一级毛片| 国产白丝娇喘喷水9色精品| 国产色爽女视频免费观看| 老司机福利观看| 又爽又黄无遮挡网站| 午夜福利成人在线免费观看| 在线观看午夜福利视频| 精品久久久久久久末码| 一级毛片我不卡| 亚洲最大成人手机在线| 简卡轻食公司| 国产亚洲欧美98| 亚洲最大成人手机在线| 99久久无色码亚洲精品果冻| videossex国产| 18禁裸乳无遮挡免费网站照片| 日本一二三区视频观看| 日本爱情动作片www.在线观看 | 亚洲一区高清亚洲精品| 好男人在线观看高清免费视频| 久久久久九九精品影院| 国产成人a∨麻豆精品| 精品久久久久久久末码| 给我免费播放毛片高清在线观看| 亚洲av免费高清在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲婷婷狠狠爱综合网| 99久久精品国产国产毛片| 人妻少妇偷人精品九色| 久久人人爽人人片av| 中国美女看黄片| 国产精品女同一区二区软件| 99热只有精品国产| 国产真实伦视频高清在线观看| 国产一级毛片七仙女欲春2| 国产精品一区二区三区四区免费观看 | 波多野结衣巨乳人妻| 国产大屁股一区二区在线视频| 男人狂女人下面高潮的视频| 好男人在线观看高清免费视频| 99久久无色码亚洲精品果冻| 免费看美女性在线毛片视频| 欧美精品国产亚洲| 国产精品av视频在线免费观看| 蜜臀久久99精品久久宅男| 我的老师免费观看完整版| 国语自产精品视频在线第100页| 91在线观看av| 国产探花在线观看一区二区| 精品久久久久久久久久久久久| 欧美高清性xxxxhd video| 舔av片在线| 九九久久精品国产亚洲av麻豆| 久久亚洲精品不卡| 国产精品国产三级国产av玫瑰| 99在线人妻在线中文字幕| 噜噜噜噜噜久久久久久91| 日日干狠狠操夜夜爽| 亚洲婷婷狠狠爱综合网| 免费搜索国产男女视频| 久久韩国三级中文字幕| 给我免费播放毛片高清在线观看| 亚洲成人久久爱视频| 色视频www国产| 欧美xxxx性猛交bbbb| 亚洲av免费高清在线观看| 日韩av在线大香蕉| 1024手机看黄色片| av女优亚洲男人天堂| 精品欧美国产一区二区三| 国产片特级美女逼逼视频| 亚洲一区高清亚洲精品| 国产69精品久久久久777片| 美女 人体艺术 gogo| 亚洲精品456在线播放app| 午夜久久久久精精品| 国产亚洲精品综合一区在线观看| 亚洲av五月六月丁香网| 欧美绝顶高潮抽搐喷水| 色5月婷婷丁香| 久久久久久久久中文| 精品国内亚洲2022精品成人| 淫妇啪啪啪对白视频| 亚洲国产精品国产精品| 国产久久久一区二区三区| 极品教师在线视频| 干丝袜人妻中文字幕| 美女高潮的动态| 最新在线观看一区二区三区| 少妇人妻精品综合一区二区 | 97人妻精品一区二区三区麻豆| 日本五十路高清| 又粗又爽又猛毛片免费看| 日韩人妻高清精品专区| 欧美一区二区精品小视频在线| 夜夜看夜夜爽夜夜摸| 亚洲中文字幕日韩| 精品人妻偷拍中文字幕| avwww免费| 国产欧美日韩精品亚洲av| 国产精品久久久久久av不卡| 成人无遮挡网站| 国产精品嫩草影院av在线观看| 国产日本99.免费观看| 欧洲精品卡2卡3卡4卡5卡区| 成人美女网站在线观看视频| 欧美激情国产日韩精品一区| 国产伦在线观看视频一区| 国产黄片美女视频| 永久网站在线| 精品人妻视频免费看| 国内揄拍国产精品人妻在线| 久久这里只有精品中国| 香蕉av资源在线| 亚洲18禁久久av| 少妇的逼好多水| 成人午夜高清在线视频| 日韩亚洲欧美综合| 热99re8久久精品国产| 高清毛片免费看| 国产激情偷乱视频一区二区| 国产精品不卡视频一区二区| 少妇熟女aⅴ在线视频| 欧美日韩一区二区视频在线观看视频在线 | 日本一本二区三区精品| 大又大粗又爽又黄少妇毛片口| 联通29元200g的流量卡| 天堂av国产一区二区熟女人妻| 青春草视频在线免费观看| 欧美潮喷喷水| 欧美日本视频| 国产免费男女视频| 白带黄色成豆腐渣| 亚洲精品一卡2卡三卡4卡5卡| 看十八女毛片水多多多| av福利片在线观看| 日本-黄色视频高清免费观看| 国内精品美女久久久久久| 我的老师免费观看完整版| 亚洲乱码一区二区免费版| 搡老岳熟女国产| 91麻豆精品激情在线观看国产| a级毛片a级免费在线| 欧美成人一区二区免费高清观看| 日本撒尿小便嘘嘘汇集6| 亚洲欧美日韩高清在线视频| 亚洲性久久影院| 一进一出抽搐gif免费好疼| 精品久久久久久久久亚洲| 久久久国产成人精品二区| 午夜老司机福利剧场| 一边摸一边抽搐一进一小说| 最近2019中文字幕mv第一页| 男人狂女人下面高潮的视频| 在线播放国产精品三级| 在线播放无遮挡| 亚洲第一区二区三区不卡| av国产免费在线观看| 少妇的逼好多水| 久久久久九九精品影院| 别揉我奶头 嗯啊视频| 欧美日韩国产亚洲二区| 国产白丝娇喘喷水9色精品| 久久久午夜欧美精品| 亚洲精品色激情综合| 亚洲精品影视一区二区三区av| 成年免费大片在线观看| 欧美潮喷喷水| 日日摸夜夜添夜夜添av毛片| 亚洲最大成人中文| 亚洲国产高清在线一区二区三| 国产乱人视频| 少妇丰满av| 久久中文看片网| 亚洲精品国产成人久久av| 久久6这里有精品| 亚洲欧美清纯卡通| 99热全是精品| 亚洲久久久久久中文字幕| 精品久久久久久久久久免费视频| 成人精品一区二区免费| 日韩精品中文字幕看吧| 国产一级毛片七仙女欲春2| 神马国产精品三级电影在线观看| 男人舔女人下体高潮全视频| 精品乱码久久久久久99久播| 国产一区二区三区在线臀色熟女| 免费人成视频x8x8入口观看| 免费观看的影片在线观看| 国产精品一二三区在线看| 欧美国产日韩亚洲一区| 一级毛片aaaaaa免费看小| 久久久久国产网址| 此物有八面人人有两片| 久久精品国产自在天天线| 国产国拍精品亚洲av在线观看| 一区二区三区免费毛片| 老司机福利观看| 国产探花极品一区二区| 欧美日本亚洲视频在线播放| a级毛片免费高清观看在线播放| 亚洲国产精品成人久久小说 | 人妻少妇偷人精品九色| 免费一级毛片在线播放高清视频| 国产极品精品免费视频能看的| 成人毛片a级毛片在线播放| 丰满的人妻完整版| 欧美激情在线99| 欧美色欧美亚洲另类二区| 国产精品野战在线观看| 亚洲一区二区三区色噜噜| 搡老岳熟女国产| 91久久精品国产一区二区三区| 高清毛片免费观看视频网站| 亚洲第一电影网av| 三级经典国产精品| 国产一区二区在线观看日韩| 男女之事视频高清在线观看| 啦啦啦啦在线视频资源| 亚洲欧美成人精品一区二区| 色在线成人网| 久久综合国产亚洲精品| 特大巨黑吊av在线直播| 伦理电影大哥的女人| 最好的美女福利视频网| 国产一区二区激情短视频| 亚洲欧美中文字幕日韩二区| 国产又黄又爽又无遮挡在线| 欧美日本视频| av在线蜜桃| 丰满的人妻完整版| 校园人妻丝袜中文字幕| 午夜爱爱视频在线播放| 久久久久精品国产欧美久久久| 欧美激情国产日韩精品一区| 国产高潮美女av| 九九久久精品国产亚洲av麻豆| 美女xxoo啪啪120秒动态图| 久久久色成人| 国产精品永久免费网站| 国产精品99久久久久久久久| 成人毛片a级毛片在线播放| 国产精品久久久久久精品电影| 在线观看美女被高潮喷水网站| 好男人在线观看高清免费视频| 亚洲美女搞黄在线观看 | 美女黄网站色视频| 岛国在线免费视频观看| 在线看三级毛片| 精品一区二区三区av网在线观看| 白带黄色成豆腐渣| 日韩欧美精品v在线| 国国产精品蜜臀av免费| 亚洲精品国产av成人精品 | 国产精品福利在线免费观看| 国产黄a三级三级三级人| 美女cb高潮喷水在线观看| 欧美高清性xxxxhd video| 国产免费男女视频| 婷婷精品国产亚洲av在线| 亚洲熟妇熟女久久| 国产亚洲av嫩草精品影院| 成人国产麻豆网| 久久综合国产亚洲精品| 51国产日韩欧美| 黄色欧美视频在线观看| 小说图片视频综合网站| 亚洲欧美成人综合另类久久久 | 波多野结衣巨乳人妻| 三级国产精品欧美在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品一二三区在线看| 真人做人爱边吃奶动态| 亚洲欧美精品自产自拍| 久久精品国产清高在天天线| 国产精品久久视频播放| 校园春色视频在线观看| 成人永久免费在线观看视频| 此物有八面人人有两片| 精品久久久久久久久av| 欧美色视频一区免费| 国产色爽女视频免费观看| 欧美区成人在线视频| 国产蜜桃级精品一区二区三区| 成人午夜高清在线视频| 亚洲精品456在线播放app| 欧美日韩国产亚洲二区| 婷婷精品国产亚洲av在线| 亚洲国产欧美人成| 在现免费观看毛片| 乱码一卡2卡4卡精品| 欧美区成人在线视频| 久久精品国产清高在天天线| 精品久久久久久久末码| 亚洲美女搞黄在线观看 | 男人舔女人下体高潮全视频| 日本爱情动作片www.在线观看 | 国产三级中文精品| 国产一区二区亚洲精品在线观看| 国产精品嫩草影院av在线观看| 欧美一级a爱片免费观看看| 日本-黄色视频高清免费观看| 97超视频在线观看视频| 男插女下体视频免费在线播放| 直男gayav资源| 三级男女做爰猛烈吃奶摸视频| 婷婷六月久久综合丁香| 18禁在线无遮挡免费观看视频 | 少妇人妻一区二区三区视频| 日韩欧美免费精品| 成人二区视频| 天天一区二区日本电影三级| 国产真实伦视频高清在线观看| 观看美女的网站| 老司机午夜福利在线观看视频| 国产精品亚洲美女久久久| 自拍偷自拍亚洲精品老妇| 久久综合国产亚洲精品| 国产单亲对白刺激| 日韩欧美 国产精品| 国产精品三级大全| 美女高潮的动态| 国产伦一二天堂av在线观看| 午夜福利18| 日本a在线网址| 亚洲人成网站在线播| av福利片在线观看| 免费电影在线观看免费观看| 日本一二三区视频观看| 精品久久久久久久久久免费视频| 我的老师免费观看完整版| 少妇的逼好多水| 热99re8久久精品国产| 在线播放无遮挡| 亚洲精品日韩av片在线观看| 狠狠狠狠99中文字幕| 欧美高清性xxxxhd video| 免费av不卡在线播放| 国产私拍福利视频在线观看| 久久九九热精品免费| 亚洲丝袜综合中文字幕| 亚洲精品亚洲一区二区| 国产精品女同一区二区软件| 亚洲欧美精品综合久久99| 国产黄片美女视频| 日韩欧美三级三区| 国产精品一区www在线观看| 欧美日本视频| 色噜噜av男人的天堂激情| 在线观看午夜福利视频| 亚洲欧美日韩高清在线视频| 1000部很黄的大片| 九九热线精品视视频播放| 美女高潮的动态| 国产精品一区二区三区四区免费观看 | 国产成年人精品一区二区| 国产精品一区www在线观看| 欧美日本视频| 日韩精品有码人妻一区| 日本一二三区视频观看| 欧美潮喷喷水| 午夜福利18| 天天躁日日操中文字幕| 亚洲人成网站高清观看| 大型黄色视频在线免费观看| 日韩欧美免费精品| 日本三级黄在线观看| av在线蜜桃| 美女xxoo啪啪120秒动态图| 精品99又大又爽又粗少妇毛片| 国产亚洲精品久久久久久毛片| 有码 亚洲区| 日韩欧美国产在线观看| 日韩精品青青久久久久久| 欧美xxxx黑人xx丫x性爽| 麻豆一二三区av精品| 久久亚洲精品不卡| 国产美女午夜福利| 黄色一级大片看看| 色哟哟·www| 国产v大片淫在线免费观看| 色哟哟哟哟哟哟| 国产高清有码在线观看视频| 中文在线观看免费www的网站| 国内少妇人妻偷人精品xxx网站| 国产av在哪里看| 欧美日韩精品成人综合77777| 少妇猛男粗大的猛烈进出视频 | 又爽又黄无遮挡网站| 亚洲电影在线观看av| 亚洲精品日韩在线中文字幕 | 91麻豆精品激情在线观看国产| 久久国内精品自在自线图片| 人人妻人人看人人澡| 欧美又色又爽又黄视频| 日韩高清综合在线| 欧美人与善性xxx| 欧美精品国产亚洲| 欧洲精品卡2卡3卡4卡5卡区| 欧美日本亚洲视频在线播放| 给我免费播放毛片高清在线观看| 国产白丝娇喘喷水9色精品| av在线天堂中文字幕| 免费观看精品视频网站| 草草在线视频免费看| 能在线免费观看的黄片| 精品福利观看| 99九九线精品视频在线观看视频| 久久久精品大字幕| 99久久精品热视频| 欧美3d第一页| 日本免费一区二区三区高清不卡| av卡一久久| av在线亚洲专区| 少妇的逼好多水| 国产精品久久久久久久久免| 国产免费一级a男人的天堂| 国产精品久久视频播放| 伦精品一区二区三区| 久久人人爽人人爽人人片va| 亚洲国产精品国产精品| 少妇人妻一区二区三区视频| 少妇人妻精品综合一区二区 | 国产精品伦人一区二区| 成人特级黄色片久久久久久久| 麻豆国产97在线/欧美| 蜜桃亚洲精品一区二区三区| 国产精品一区二区三区四区免费观看 | 国产一区二区在线av高清观看| 日本黄色片子视频| 国产视频一区二区在线看| 床上黄色一级片| 色尼玛亚洲综合影院| 天天躁夜夜躁狠狠久久av| 日本爱情动作片www.在线观看 | 色噜噜av男人的天堂激情| 久久人人爽人人片av| 亚洲欧美精品自产自拍| 一夜夜www| 亚洲欧美日韩高清在线视频| 久久精品夜夜夜夜夜久久蜜豆| 99在线视频只有这里精品首页| 精品熟女少妇av免费看| 卡戴珊不雅视频在线播放| 亚洲国产日韩欧美精品在线观看| 一进一出好大好爽视频| 少妇熟女欧美另类| 国产成人精品久久久久久| 国模一区二区三区四区视频| 日韩,欧美,国产一区二区三区 | 欧洲精品卡2卡3卡4卡5卡区| 深夜a级毛片| 免费在线观看成人毛片| 国内精品美女久久久久久| 性色avwww在线观看| 波多野结衣巨乳人妻| 久久精品综合一区二区三区| 国产在线精品亚洲第一网站| 国产精品一及| 国产乱人偷精品视频| 国产午夜精品久久久久久一区二区三区 | 性色avwww在线观看| 亚洲av五月六月丁香网| 嫩草影院入口| 欧美最新免费一区二区三区| 欧美激情在线99| 婷婷精品国产亚洲av| 欧美一区二区亚洲| 国产高清三级在线| 男人舔女人下体高潮全视频| 成人鲁丝片一二三区免费| 12—13女人毛片做爰片一| 国内精品宾馆在线| 在线a可以看的网站| 天美传媒精品一区二区| 欧美高清成人免费视频www| 欧美精品国产亚洲| 天堂影院成人在线观看| 看非洲黑人一级黄片| 桃色一区二区三区在线观看| 国产爱豆传媒在线观看| 国产亚洲精品久久久com| 亚洲av五月六月丁香网| 校园人妻丝袜中文字幕| 欧美日本视频| 色综合站精品国产| 青春草视频在线免费观看| 九九久久精品国产亚洲av麻豆| 亚洲图色成人| av国产免费在线观看| 桃色一区二区三区在线观看| av天堂在线播放| 国产麻豆成人av免费视频| 1000部很黄的大片| 春色校园在线视频观看| 国产黄片美女视频| 国产一区二区亚洲精品在线观看| 一本精品99久久精品77| avwww免费| 国产久久久一区二区三区| 日日啪夜夜撸| 老司机影院成人| 九九在线视频观看精品| 欧美xxxx黑人xx丫x性爽| 色综合色国产| 日本爱情动作片www.在线观看 | 日韩欧美三级三区| 日本一本二区三区精品| 亚洲国产欧洲综合997久久,| 99热精品在线国产| 色哟哟·www| 赤兔流量卡办理| 激情 狠狠 欧美| 在线国产一区二区在线| 91久久精品国产一区二区三区| 18+在线观看网站| 国产精品精品国产色婷婷| 91狼人影院| 久久久久久国产a免费观看| 精品久久久久久久久久免费视频| 久久精品国产亚洲av香蕉五月| 欧美另类亚洲清纯唯美| 精品久久久久久久久亚洲| 久久草成人影院| 日产精品乱码卡一卡2卡三| 欧美3d第一页| eeuss影院久久| 亚洲国产日韩欧美精品在线观看| 少妇的逼好多水| 麻豆一二三区av精品| 淫秽高清视频在线观看| 亚洲成av人片在线播放无| 黄色视频,在线免费观看| 国产免费男女视频| 一个人观看的视频www高清免费观看| 欧美三级亚洲精品| 麻豆国产av国片精品| 97碰自拍视频| 成人二区视频| 日韩中字成人| 天堂影院成人在线观看| 日韩三级伦理在线观看| 日韩一区二区视频免费看| 黄色视频,在线免费观看| 内射极品少妇av片p| 两个人视频免费观看高清| 亚洲经典国产精华液单| 插阴视频在线观看视频| 欧美三级亚洲精品| 乱系列少妇在线播放| 99国产极品粉嫩在线观看| 国产精品伦人一区二区| 午夜a级毛片| 国产大屁股一区二区在线视频| 国产成年人精品一区二区|