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

    Scaling procedures of cabin noise generated by turbulent boundary layer excitation

    2018-04-19 08:29:06XiojinZHAOWeiSONGJieZHOUBngchengAI
    CHINESE JOURNAL OF AERONAUTICS 2018年3期

    Xiojin ZHAO,Wei SONG,Jie ZHOU,Bngcheng AI

    aChina Academy of Aerospace Aerodynamics,Beijing 100074,China

    bUniversity of Southampton,Southampton,England SO17IBJ,United Kingdom

    1.Introduction

    Interior noise alleviation is an issue of primary interest in the design of modern aircraft.The acoustic nuisance inside the cabin may adversely influence passenger comfort,particularly in helicopter and propeller-driven aircraft.The typical sources for generating cabin noise include the fuselage boundary layer and airborne and structure-borne noise,which may cause unacceptable ride discomfort and impact the fatigue life of the structure.

    Investigations on the cabin noise induced by Turbulent Boundary Layers(TBLs)began early on.In Ref.1,the contribution of TBLs to cabin noise was studied,and the problem was formulated as a sequence of two linear couplings,namely,the TBL fluctuations exciting the fuselage skin in lateral vibrations and the skin vibrations inducing sound inside the fuselage.Moreover,it was assumed that the boundary layer was locally homogeneous and the fuselage skin was flat,and only outgoing waves were considered.Airplane interior broadband noise under the cruise condition was considered to be dominated by the TBL.For the prediction of cabin noise,a mathematical model2was derived from the Statistical Energy Analysis(SEA)technique to evaluate the interior TBL noise of a Boeing 737 airplane.A method for measuring the in- flight noise and vibration of an aircraft,developed by SAAB(Svenska Aeroplan AktieBolaget),3was shown to be a powerful tool for monitoring the quality of an aircraft.In Ref.4,ground and flight measurements were conducted on a pusher-propeller configuration with different interior furnishings.The main objective of the measurements was to identify the sources of interior noise on a specific airplane and revise the computer program so that it could be used to predict the noise effects of interior treatment modifications.A new wavenumberbased approach for predicting sound transmission through an aircraft fuselage was introduced in Ref.5.The proposed method was accurate at mid-frequencies,for which the method involved nearly no approximation for a large doubly periodic panel.A new stochastic model of simulating the surface pressure fluctuations that induced cabin noise was developed in Ref.6.The method proved to be efficient when a random particle mesh method with recursive filters was used.Ref.7investigated the cabin noise generated by engine noise,and the proposed methodology considered the engine as a spatially extended source and assumed that the jet-mixing source was a line source distributed on the engine axis.With a focus on the cabin noise induced by TBL excitation,the interior noise environment in a Large Civil TiltrotoR(LCTR2)was assessed in Ref.8.The contributions of different aero-acoustic sources to the cabin noise of an aircraft,including the TBL,jet noise and the air conditioning system,were compared in Ref.9.Extensive measurements with microphones mounted at various longitudinal positions were used to assess the cabin noise.Ref.10presented the results of flight measurements,including the field of pressure fluctuation on the fuselage and cabin noise,and identified the main sources of noise.

    The methods developed for assessing cabin noise typically only provide tools to support the noise control in the cabin of the aircraft.Ref.11reported the reduction of the engine rotor vibration-induced cabin noise of DC-9 due to several improvement designs.Ref.12investigated the control of lowfrequency cabin noise and proposed a concept based on intrinsic tuning and damping of fuselage structural elements.Ref.13discussed the performance of various structural and cabin sidewall treatments applied to reduce cabin noise,and the measurement and analysis were conducted on a DC-9 test section.Passive vibration absorbers are often ineffective for the analysis of frequencies below 500 Hz because the wavelength is large compared to the thickness of a damping layer.Thus,an active control system for cabin noise consisting of actuators and sensors was used in Refs.14–17.

    The precision of the simulated TBL excitation typically limits the cabin noise assessment.A wind tunnel test is the most effective method for simulating TBL excitations over an aircraft with a scaled model.Thus,this paper proposes to assess the cabin noise of a craft induced by TBL excitation in a wind tunnel.Moreover,scaling laws are developed to relate the cabin noise of a scaled model with the full-scale model.The remainder of the paper is organized as follows.Section 2 describes the cabin noise scaling procedures,which are derived from the theoretical models of cylindrical shell vibrations induced by TBL excitations.Furthermore,to validate the derived scaling procedures,theoretical results are presented in Section 3.In Section 4,the measurements are introduced,and the experimental results are discussed in Section 5.The last section provides conclusions concerning the work in this paper.

    2.A scaling procedure for cylindrical shell response

    The scaling procedures used to predict the structural response of a typical aircraft when subjected to TBL excitations were derived based on a curved cylindrical shell(Fig.1).In Fig.1,a and b indicate the length and the width of the cylindrical shell,respectively,and R is the curvature.The scaling laws of the structure response induced by TBL excitation introduced in Section 2.1 were developed in Refs.18,19Because the interior sound distribution is also an important factor that must be considered in the structural design of an actual aircraft,the other scaling procedure was developed in Section 2.2 to relate the cabin noise distribution of scaled aircraft structures with that of full-scale aircraft structures.The scaling proceduresdeveloped can provide theoreticalsupportfor predicting cabin noise in a wind tunnel.

    2.1.Scaling procedure of structure response

    The differential equations governing the vibration of a curved cylindrical shell(as presented in Refs.20,21)are as follows:

    where w is the transverse displacement of panel.

    μ= ρh,which indicates the area density;h is the thickness of the curved cylindrical shell and ρ is the density of structure.D=Eh3/12(1- ν2),where v is Poisson’s ratio of the shell,E is Young’s modulus;pi=exp[j(ωt-kxx-kyy-kzz)],which indicates an incident wave acting on the shell,kx,ky,kzare the wave number along x,y,z,respectively,prdenotes the acoustic pressure radiated by the shell,ζ is Airy’s stress function,ω is circular excitation frequency,t is the time,indicates the second-order operator,▽4indicates the forth-order operator,j is the imaginary unit.

    If the scaling coefficient is denoted as σ,each side of the scaled curved cylindrical shell can be expressed asˉa=σa,ˉb=σb,andˉh=σh.The parameters with‘ˉ’correspond to the scaled model.

    The relationship characterizing the velocity Power Spectral Density(PSD)response between the scaled model and the fullscale model is as follows:18,19

    Fig.1 Schematic of cylindrical shell.

    In addition,Svis the velocity power spectral density of the structure,which is circular frequencyandindicatesthe response frequency of the scaled model,V is free stream velocity,q∞is dynamic pressure of the flow.

    2.2.Scaling procedure for internal sound

    The internal sound assessment of the aircraft is based on a weak coupling approximation.The flow induces the structural vibration,and then the structural vibration determines the internal acoustic pressure;however,the counterforce of the sound pressure applied to the structure is ignored.The structural vibration of a curved cylindrical shell is more complex than the sound radiation of a flat plate.The internal sound pressure of an aircraft is typically determined based on two important sound source components:direct sound21and reverberation sound.The sound pressure at an arbitrary point in the cabin of an aircraft is expressed as

    When the aircraft is scaled,the direct sound can be expressed as

    According to the research presented in the paper discussed in Section 2.1,

    where K indicates the velocity response relationship between the scaled model and the original structure,which is determined by the coefficients of the derived scaling procedures presented in Section 2.1.Eq.(3)presents the velocity PSD response relationship between the scaled model and the original structure.Therefore,the coefficient K could be obtained by combining Eq.(5)with Eq.(3).

    Thus,

    From Eq.(6),a scaling relationship between the direct sound of the scaled model and that of a full-dimensional model is obtained if either=1 or G=is satisfied.The scaling laws can be expressed as

    In addition,the reverberation sound also satisfies a similar formula:

    Thus,the following expression can be obtained:

    According to the research in Section 2.1,is obtained,q is the dynamical pressure of air flow.

    3.Theoretical results

    To validate the scaling laws developed regarding the cabin noise induced by TBL excitation,the Finite Element Method(FEM)solver in the VA-One software was used to compute the cabin noise of a typical aircraft layout(Fig.2).The Corcos model used to simulate the TBL excitation can be expressed as20

    where Spis power spectral density of the wall pressure distribution;Ucis the convection velocity;ξxand ξyindicate the separation distances of two arbitrary points in the stream wise and span wise directions on the panel,respectively;αxand αyare the decay rates,Sppis the amplitude of the wall pressure fluctuation.

    Fig.2 Distribution of test points corresponding to original model.

    The cabin noises of both the scaled model and the full-scale model were computed for comparison.The thickness of the aircraft panel was 4 mm.The cabin noise is typically determined by the vibration responses of the structure investigated.Thus,two test points on the surface of the model,i.e.,P1 and P2,were used to measure the structural vibration.The cabin noises were measured at points P3,P4 and P5.The investigated full-scale model was made of aluminum(E=7.0×1010Pa;ρ=2.7×103kg/m3;ν=0.33),and a scaling coefficient of σ=0.25 was chosen.According to the scaling laws derived,the scaled models must also be made of aluminum to satisfy the scaling procedures.The TBL excitation was simulated using the Corcos model,where αx=0.16, αy=0.46,ρ0=1.21 kg/m3,Spp=1 Pa,and Uc=0.7U0,with U0being the velocity of the mean flow;here,U0=30 m/s.Usually,the Corcos model is only suitable for low-speed flow.However,the TBL excitation used in the scaling procedure is not related to the exact expression of TBL excitation.The Corcos model used here provides only a kind of expression for the TBL excitation,which does not weaken the rationality of the scaling procedure developed.The shell elements were the FEM elements used for computation,and the grid number for the structure was 4753,which ensured that at least 8 nodes existed for a sound wave.Hexahedron grids were applied to compute the internal sound,and the number of acoustic grids used was 1240.The same computation and the comparison of the structural responses were presented in Ref.19However,the works in this study focused on the cabin noise prediction.

    According to Eq.(4),the cabin noise is mainly determined based on the structure response induced by TBL excitation.However,the results of the structure response were previously given in Ref.19and thus will not be re-introduced here.Fig.3 presents the Sound Pressure Level(SPL)in the cabin of the full-scale model.One result was obtained using the FEM with the full-scale model;the other result was computed using the scaling procedure with the scaled model.Fig.3 shows that the mean sound levels computed using the two different methods are in good agreement.However,the local sound levels of the two methods,especially for low-frequency noise,have some difference.The reason for this difference is that the computation error usually induces a different cabin noise distribution,and the nonuniformity is more obvious for low-frequency noise because of its lower mode number in bands compared to that of high-frequency noise.However,the mean noise computation can weaken the nonuniform cabin noise distribution.Therefore,the mean noise offers a better comparison compared to the local noise.In Fig.3,the results based on the scaled model are compared with those obtained with the fullscale model,thus verifying the validity of the scaling procedure derived for the cabin noise prediction.

    4.Measurement

    The numerical methods typically only provide validation of the scaling equations in theory.For a real structure,many factors might affect the structural vibrations,including the design,manufacturing and installation,resulting in differences from the theoretical response.Thus,two types of experiments were designed to validate the scaling laws derived:a modal test and a wind tunnel test.The modal test was used to measure the modal characteristics of the structure,including the scaled model and the original model.A wind tunnel was used to simulate the TBL excitation of the structure on both the scaled model and the original model.

    4.1.Modal test

    A modal test was performed to determine the natural vibration of the structure and to provide theoretical support for the revision of the results of the wind tunnel test.For the modal test,a full-scale model with a thickness of 3 mm,which was the same as that of the wind tunnel test model,was made of steel(E=2.1×1011Pa;ρ=7.8×103kg/m3;ν=0.313).The scaling coefficient for the scaled model,which was made of aluminum(E=7.0×1010Pa;ρ=2.7×103kg/m3;ν=0.33),was 0.5.

    The test points used to measure the structural vibrations and noise distributions in the test model are shown in Fig.4.The dynamic response of the model was measured under the excitation of an external force generated by a vibration exciter(YMC,Sharker)at F1 and F2(Fig.4(a)).The excitation signal was generated by SIGLENT SDG1005 and magnified with a power amplifier(YMC,LA-100).The vibration was detected by accelerometers(PCB,352C33)mounted at P1 and P2 on the plate shown in Fig.4,and the internal noise was measured using the microphone(BSWA,MPA 416)located at P3.The signals for the structural vibration and internal noise were collected by a collector(VXI-16026A)and then transferred to a computer.

    4.2.Wind tunnel test

    Fig.3 Sound level prediction with a scaled model under Condition 1.

    Fig.4 Experimental test design and model.

    The experiment performed in the wind tunnel was designed to utilize the wind tunnel to simulate the TBL excitation more effectively than other methods.The TBL excitation for transonic flows22–26is typically more serious.Hence,the TBL excitation of transonic flows was simulated in the wind tunnel to verify the developed scaling procedure.The cabin noise of the full-scale model induced by over flow was compared with the results obtained using the scaling laws with a scaled model to validate the scaling procedures.

    The TBL excitation was simulated in the FD-12 wind tunnel,which was built at the China Academy of Aerospace Aerodynamics(CAAA).The FD-12 wind tunnel is a type of intermittent wind tunnel.The dimensions of the wind tunnel are 1.2 m×1.2 m.The range of the flow Mach numbers simulated in FD-12 is from Ma 0.4 to Ma 4.0.

    For the wind tunnel test,the models used were the same as those used for the modal test in Section 4.1.Because the internal noise of the models was easily influenced by external noise in the wind tunnel,the end of the test model was sealed using the special structure shown in Fig.5.The model design was presented in Ref.19A small hole was maintained in the center of the airpro of structure to allow signal transmission.Simultaneously,the support of the test model was designed with sufficient stiffness to avoid model vibrations due to the impact of the air flow in the wind tunnel.For the wind tunnel test,the model was fixed in the middle of the test section to ensure that the model was not affected by the wall of the wind tunnel.Fig.6 shows the test model fixed in the wind tunnel.

    Fig.5 Design scheme of test model.

    Fig.6 Test model fixed in wind tunnel.

    Fig.7 Distribution of test points corresponding to original model.

    The test points used for the measurement of structural vibrations and noise distributions in the test model are shown in Fig.7.The two accelerometers used for the structural vibration measurement were fixed at points PP1 and PP2.Because the structural vibration comparison of the wind tunnel test results had been previously presented in Ref.19,only the cabin noise data were provided to validate the theory.Fluctuating pressure sensors used to measure the internal noise were fixed at points PP3 and PP4.The structural vibration and internal noise were measured with accelerometers(PCB,333C33)and fluctuating pressure sensors(Kulite XCL-100),respectively.When the wind tunnel was running,the signals for the structural vibration and the cabin noise were collected by a collector(VXI-16026A)and then transferred to a PC.

    5.Comparison

    Figs.8 and 9 present the velocity PSD of the model derived from the modal test.The solid line is the test results of the full-scale model,and the dash-dot line was obtained using the scaling laws with the test result of a scaled model.The purpose of the modal test is twofold.The first purpose is to validate the scaling procedure of the cabin noise;the second purpose is to test the frequency offset due to the model manufacture and installation.The frequency offset could be obtained with the test results of the structure response,and the frequency offset should be introduced to compute the cabin noise for a full-scale structure with the scaling procedure.Figs.10 and 11 show the sound pressure of the cabin under excitation at F1 and F2,respectively.The frequency offsets are considered in Figs.10 and 11.

    Fig.10 One-third octave of sound pressure for P3 under excitation at F1.

    Figs.12 and 13 show the 1/3 octave of the cabin noise measured in the wind tunnel of the aircraft model for Ma=0.85 and Ma=0.9,respectively.The results of the scaled model were obtained using the scaling procedure with a scaled model measurement.The results of the original model were obtained via direct measurement of the full-scale model.The results shown,scaled with frequency revision,were obtained using the scaling procedure with the frequency offset considered.Figs.12 and 13 clearly show that a notable frequency offset occurred between the results of the original model and those of the scaled model when the frequency offset was not considered.

    Fig.8 One-third octave of velocity PSD under excitation at F1.

    Fig.9 One-third octave of velocity PSD under excitation at F2.

    Fig.11 One-third octave of sound pressure for P3 under excitation at F2.

    The frequency offsets are usually generated for different reasons;the effect of the material used to construct the scaled model on the frequency offset was discussed in Ref.18The frequency offset due to the use of different materials can be eliminated via the frequency offset correction shown in Ref.18Meanwhile,the errors introduced during the design,manufacture and installation of the structure are also important reasons for generating frequency offsets,particularly in the low-frequency range.However,these errors can be reduced via improvements in the experimental design.This was introduced in Ref.19In this paper,the modal test is used for frequency offset correction.In Figs.12 and 13,the comparisons clearly show that the cabin sound pressure of the aircraft model with the frequency offset considered has a better agreement with the original model results than the scaled model without frequency offset correction.

    However,for the cabin sound of the model,the SPL based on the scaled model is in good agreement with the original model results only at point PP3.The background noise may be a factor affecting the sound pressure prediction in the cabin of the models.Although the end of the model was sealed using a special structural design,background noise could still have diffused into the cabin through the hole that was maintained for signal transmission in the support.Because point PP4 is closer to the support than point PP3,it is more easily affected by the background noise of the wind tunnel.The comparison at point PP3 shows more consistent results than that at point PP4.

    Fig.12 One-third octave of cabin noise measured in wind tunnel for Ma=0.85.

    Fig.13 One-third octave of cabin noise measured in wind tunnel for Ma=0.9.

    6.Conclusions

    This paper presents a method for scaling cabin sound.The scaling procedure provides an academic foundation for the prediction of the cabin noise of the structure in a wind tunnel.In addition,a new test skill is proposed for investigating cabin noise in a wind tunnel.

    The numerical computation provides theoretical results to validate the scaling laws developed for the investigation of the cabin noise generated by TBL excitation on the structure.Because the cabin sound is determined by both the external excitation and the modal characteristic of the structure,a modal test was performed before the wind tunnel test.The modal test provides an effective method to obtain the frequency offset of the structure response due to the model manufacture and installation.

    A wind tunnel test was performed to simulate the TBL excitation.The cabin sound of the full-scale model was compared to the results obtained using the scaling procedure based on the scaled model.The comparisons of the measurement results show that the scaling procedures developed in this paper are effective for vibro-acoustic predictions in wind tunnels.

    Background noise was an important factor affecting the sound pressure in the cabin of the models for the wind tunnel test.The tail end of test model should be sealed to prevent background noise from entering the cabin of the model.Although the background noise could diffuse into the cabin through the hole that was maintained for signal transmission in the support,a better result might be obtained if a better test scheme was designed.

    ?

    1.Corcos GM,Liepmann HW.On the contribution of turbulent boundary layers to the noise inside a fuselage.Washington,D.C.:NACA;1956.Report No.:NACA-TM-1420.

    2.Alex LJS.Airplane interior noise modeling using Statistical Energy Analysis Approach.Reston:AIAA;1999.Report No.:AIAA-1999-1903.

    3.Joakim H.SAAB’S experience in aircraft interior noise and vibration measurements.Reston:AIAA;1999.Report No.:AIAA-1999-1837.

    4.Francesco M,Antonia S,Marco A.Interior noise sources identifications,in- flight measurements and numerical correlations of an advanced business aircraft.Reston:AIAA;2006.Report No.:AIAA-2006-2491.

    5.Blakemore M,Heron KH,Davis EB.A new wavenumber-based approach for predicting aircraft cabin noise.Reston:AIAA;2008.Report No.:AIAA-2008-3009.

    6.Malte S,Roland E,Olaf H,Oliver U.A synthetic wall pressure model for the efficient simulation of boundary layer induced cabin noise.Reston:AIAA;2010.Report No.:AIAA-2010-3760.

    7.Alessandro B,Sebastien G.Semi empirical jet noise modeling for cabin noise prediction-acoustic loads in the geometric near field.Reston:AIAA;2011.Report No.:AIAA-2011-2925.

    8.Ferdinand WG.Large Civil Tiltrotor(LCTR2)interior noise predictions due to turbulent boundary layer excita tion.Washington,D.C.:NASA;2013.Report No.:NASA/CR-2013-218005.

    9.Nan H,Heino B,Michaela H.Contributions of different aeroacoustic sources to aircraft cabin noise.Reston:AIAA;2013.Report No.:AIAA-2013-2030.

    10.Robert GA,Alexey JG.Identification of noise sources in the cabin and determination of the sound energy local passage through the board design based on the results of in- flight measurements of the super jet.Reston:AIAA;2015.Report No.:AIAA-2015-3114.

    11.Van Dyke JD,Schendel J,Gunderson C,Ballard M.Cabin noise reduction in the DC-9.Reston:AIAA;1967.Report No.:AIAA-1967-0401.

    12.SenGupta G.Methods of reducing low frequency cabin noise and sonically induced stresses,based on the intrinsic structural tuning concept.Reston:AIAA;1977.Report No.:AIAA-1977-0444.

    13.Simpson MA.Cabin noise control ground tests for UHB aircraft.Reston:AIAA;1989.Report No.:AIAA-1989-1121.

    14.Swanson A,Billoud GD.Aircraft cabin noise reduction through active vibration control.Reston:AIAA;1999.Report No.:AIAA-1999-1935.

    15.Maury C,Gardonio P,Elliott SJ.Active control of the flowinduced noise transmitted through a panel.Reston:AIAA;2000.Report No.:AIAA-2000-2042.

    16.Julian G,Delf S,Sten B.Global active noise control in aircraft cabins.Reston:AIAA;2009.Report No.:AIAA 2009-3242.

    17.Giovanni B,Claudio T,Marco MC.A comprehensive approach for the optimal control of tiltrotor cabin noise through activelydriven piezoelectric actuators.Reston:AIAA;2013.Report No.:AIAA-2013-2032.

    18.Zhao XJ,Ai BC,Liu Z,Li D.A scaling procedure for panel vibroacoustic response induced by turbulent boundary layer.J Sound Vibrat 2016;380:165–79.

    19.Zhao XJ,Ai BC.Predicting the structural response induced by turbulent boundary layer in wind tunnel.AIAA J 2017;55(4):1221–9.

    20.Liu BL,Feng LP,Anders NA.Sound transmission through curved aircraft pannels with stringer and ring frame attachments.J Sound Vibrat 2007;300:949–73.

    21.Liu BL.Noise radiation of aircraft panels subjected to boundary layer pressure fluctuation.J Sound Vibrat 2008;314:693–711.

    22.Ezra AA,Per-erson HC.Determination of design criteria for transonic buffeting forces acting on launch vehicles.ARS launch vehicles:structures and materials conference;Arizona,USA,April 1962.

    23.Robertson JE. Unsteady pressure phenomena for basic missle shapesat speeds.Reston:AIAA;1964. Report No.: AIAA-1964-0003.

    24.Riddle D.Investigation of surface pressure fluctuation association with buffet.Reston:AIAA;1975.Report No.:AIAA-1975-0067.

    25.Serpa JM,Lessmann RC,Hagist WM.Turbulent separated and reattached flow over a curved surface.Reston:AIAA;1986.Report No.:AIAA-1986-1064.

    26.Piatak DJ,Sekula MK,Rausch RD.Comparison of ares I-X wind-tunnel derived buffet environment with flight data.Reston:AIAA;2011.Report No.:AIAA-2011-3013.

    啦啦啦观看免费观看视频高清| 女人久久www免费人成看片 | 久久欧美精品欧美久久欧美| 成人国产麻豆网| 国产伦精品一区二区三区视频9| 最新中文字幕久久久久| 久久99热这里只有精品18| 亚洲自拍偷在线| 在线观看66精品国产| 亚洲性久久影院| 夫妻性生交免费视频一级片| 一个人看的www免费观看视频| 国产精品一区二区在线观看99 | 国产伦精品一区二区三区四那| 人人妻人人看人人澡| 亚洲天堂国产精品一区在线| av在线老鸭窝| a级一级毛片免费在线观看| 麻豆一二三区av精品| av在线亚洲专区| 大香蕉久久网| 三级经典国产精品| 国产精品一区二区性色av| 精品人妻偷拍中文字幕| 久久精品91蜜桃| 国产探花极品一区二区| 看免费成人av毛片| av福利片在线观看| 我的老师免费观看完整版| 汤姆久久久久久久影院中文字幕 | 麻豆久久精品国产亚洲av| 婷婷六月久久综合丁香| 97热精品久久久久久| 99久国产av精品| 精品少妇黑人巨大在线播放 | 高清毛片免费看| 伊人久久精品亚洲午夜| 在线观看一区二区三区| 嫩草影院入口| 亚洲人成网站高清观看| 22中文网久久字幕| 国产视频首页在线观看| 国产69精品久久久久777片| 欧美高清成人免费视频www| 精品人妻偷拍中文字幕| 欧美性感艳星| 久久精品久久久久久久性| 亚洲精华国产精华液的使用体验| 中文欧美无线码| 超碰av人人做人人爽久久| 麻豆乱淫一区二区| 精品少妇黑人巨大在线播放 | 国产一区有黄有色的免费视频 | 又粗又爽又猛毛片免费看| 国产成人福利小说| 99在线视频只有这里精品首页| 国产日韩欧美在线精品| 纵有疾风起免费观看全集完整版 | 亚洲av电影在线观看一区二区三区 | 成人毛片a级毛片在线播放| 丰满人妻一区二区三区视频av| 午夜福利在线观看吧| 久久久久久久久久久丰满| 大香蕉久久网| 97超视频在线观看视频| 人妻少妇偷人精品九色| 又黄又爽又刺激的免费视频.| 免费看光身美女| 国产乱人视频| 免费观看人在逋| 性插视频无遮挡在线免费观看| 国产高清国产精品国产三级 | 精品无人区乱码1区二区| 一级毛片aaaaaa免费看小| 国产视频内射| 女人被狂操c到高潮| 久久久久久久久久成人| 精华霜和精华液先用哪个| 九九热线精品视视频播放| 亚洲欧美日韩卡通动漫| 免费不卡的大黄色大毛片视频在线观看 | 亚洲,欧美,日韩| 天堂√8在线中文| 中文字幕av成人在线电影| 人人妻人人看人人澡| 久久久国产成人免费| 丰满乱子伦码专区| 日韩一本色道免费dvd| 天天一区二区日本电影三级| 亚洲色图av天堂| 亚洲欧美精品综合久久99| 日韩成人av中文字幕在线观看| 国产精品人妻久久久影院| 又粗又爽又猛毛片免费看| 免费无遮挡裸体视频| 亚洲最大成人中文| 可以在线观看毛片的网站| 国产乱来视频区| 国产亚洲一区二区精品| 国产成年人精品一区二区| 国产乱人视频| 国产午夜福利久久久久久| 国产探花在线观看一区二区| 欧美xxxx性猛交bbbb| 国产乱来视频区| 色尼玛亚洲综合影院| 简卡轻食公司| 国产精品嫩草影院av在线观看| 亚洲人成网站在线播| 精品久久久久久久人妻蜜臀av| 夫妻性生交免费视频一级片| 免费电影在线观看免费观看| 97热精品久久久久久| 十八禁国产超污无遮挡网站| 欧美日韩一区二区视频在线观看视频在线 | 久久久久精品久久久久真实原创| 一区二区三区四区激情视频| 亚洲一区高清亚洲精品| 国产伦一二天堂av在线观看| 日韩精品有码人妻一区| 一二三四中文在线观看免费高清| 久久久久久久国产电影| 精品熟女少妇av免费看| 日日撸夜夜添| 观看美女的网站| 国产v大片淫在线免费观看| 亚洲人成网站在线观看播放| 一区二区三区高清视频在线| 青春草亚洲视频在线观看| 午夜激情福利司机影院| 午夜免费男女啪啪视频观看| 精品久久久久久久久亚洲| 欧美性猛交╳xxx乱大交人| 久久韩国三级中文字幕| 午夜精品一区二区三区免费看| 欧美日韩综合久久久久久| 久久久久久九九精品二区国产| 精品久久久久久久久av| 少妇熟女欧美另类| 九九在线视频观看精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 村上凉子中文字幕在线| 一区二区三区四区激情视频| 又黄又爽又刺激的免费视频.| 久久热精品热| 国产老妇女一区| 天天一区二区日本电影三级| 九九在线视频观看精品| 亚洲av电影在线观看一区二区三区 | 在线天堂最新版资源| 婷婷色麻豆天堂久久 | kizo精华| 国产淫片久久久久久久久| 一夜夜www| 一个人观看的视频www高清免费观看| 免费黄色在线免费观看| 精品一区二区三区视频在线| 亚洲欧美中文字幕日韩二区| 国产av码专区亚洲av| 午夜福利成人在线免费观看| 久久99热这里只频精品6学生 | 精品一区二区免费观看| 久久久久精品久久久久真实原创| 欧美bdsm另类| 在线免费观看的www视频| 国产黄a三级三级三级人| 国产亚洲91精品色在线| 一级毛片久久久久久久久女| 在线免费观看不下载黄p国产| 99久国产av精品国产电影| 男插女下体视频免费在线播放| 最近2019中文字幕mv第一页| 久久久国产成人精品二区| 国产精华一区二区三区| 人人妻人人澡欧美一区二区| 日本熟妇午夜| 国产精品不卡视频一区二区| 啦啦啦观看免费观看视频高清| 变态另类丝袜制服| 中文精品一卡2卡3卡4更新| 午夜a级毛片| 日韩三级伦理在线观看| 中文欧美无线码| 免费电影在线观看免费观看| 噜噜噜噜噜久久久久久91| 女人十人毛片免费观看3o分钟| 免费观看的影片在线观看| 日韩一区二区三区影片| 老司机福利观看| 久久人妻av系列| av在线亚洲专区| 国产精品一区二区在线观看99 | 一个人免费在线观看电影| 久久久午夜欧美精品| 麻豆一二三区av精品| 国产精品不卡视频一区二区| 成人亚洲精品av一区二区| 99热这里只有精品一区| 国产又黄又爽又无遮挡在线| 国产成人午夜福利电影在线观看| 国产黄片美女视频| 高清午夜精品一区二区三区| 久久热精品热| 国产亚洲最大av| 丰满少妇做爰视频| 1024手机看黄色片| 三级毛片av免费| 亚洲四区av| 18禁裸乳无遮挡免费网站照片| 日韩高清综合在线| 一级av片app| 国产在视频线在精品| 97超碰精品成人国产| 国产黄片美女视频| 久久精品久久久久久噜噜老黄 | 麻豆国产97在线/欧美| 国产成人精品婷婷| 久久久久久久久久黄片| eeuss影院久久| 97在线视频观看| 国产黄片视频在线免费观看| 午夜精品在线福利| 日本黄色片子视频| 国产亚洲一区二区精品| 亚洲av中文字字幕乱码综合| 中文资源天堂在线| 欧美激情国产日韩精品一区| 91狼人影院| 少妇高潮的动态图| 91精品一卡2卡3卡4卡| 久久精品91蜜桃| 亚洲精品久久久久久婷婷小说 | 午夜a级毛片| 亚洲av一区综合| 午夜亚洲福利在线播放| 激情 狠狠 欧美| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 观看美女的网站| 午夜爱爱视频在线播放| 国产免费又黄又爽又色| av福利片在线观看| 日本五十路高清| 国内揄拍国产精品人妻在线| 日日摸夜夜添夜夜爱| 在线免费十八禁| 乱码一卡2卡4卡精品| 少妇的逼水好多| 色网站视频免费| 深夜a级毛片| 亚洲av福利一区| 2021天堂中文幕一二区在线观| 日本色播在线视频| 亚洲一级一片aⅴ在线观看| 我要搜黄色片| 亚洲最大成人手机在线| 国产一级毛片在线| 久久久精品大字幕| 国产 一区精品| 最近手机中文字幕大全| 国产91av在线免费观看| 欧美日韩国产亚洲二区| 国产av码专区亚洲av| 九九在线视频观看精品| 精品久久久久久久久亚洲| 久久人人爽人人爽人人片va| 我的老师免费观看完整版| 久久久久久久久中文| av.在线天堂| 91av网一区二区| av在线观看视频网站免费| 亚洲真实伦在线观看| 人妻夜夜爽99麻豆av| www日本黄色视频网| 国产精品一及| 精品一区二区三区视频在线| 在线a可以看的网站| 欧美一区二区亚洲| 亚洲人成网站在线观看播放| 国产一区二区在线av高清观看| 亚洲美女搞黄在线观看| 国产精品无大码| 能在线免费观看的黄片| 插逼视频在线观看| 99在线人妻在线中文字幕| 成人国产麻豆网| 国产精品99久久久久久久久| 国产午夜精品一二区理论片| 日本一本二区三区精品| 亚洲国产成人一精品久久久| 国产真实乱freesex| 在线免费十八禁| 亚洲四区av| av女优亚洲男人天堂| 日本一二三区视频观看| 女人十人毛片免费观看3o分钟| 91久久精品国产一区二区三区| 免费av毛片视频| 国产亚洲5aaaaa淫片| 1000部很黄的大片| 欧美潮喷喷水| 亚洲人成网站高清观看| 日本wwww免费看| 高清av免费在线| 亚洲婷婷狠狠爱综合网| 欧美成人一区二区免费高清观看| 一夜夜www| 亚洲自偷自拍三级| 国产精品国产三级国产av玫瑰| 欧美bdsm另类| 我的老师免费观看完整版| 小说图片视频综合网站| 精品酒店卫生间| 国产亚洲精品av在线| 伦精品一区二区三区| 亚洲人成网站高清观看| 99久久人妻综合| 男人狂女人下面高潮的视频| 国产亚洲精品久久久com| 91久久精品国产一区二区成人| 色综合亚洲欧美另类图片| 日韩中字成人| 97人妻精品一区二区三区麻豆| 午夜福利在线在线| 99视频精品全部免费 在线| 免费av不卡在线播放| 人人妻人人澡人人爽人人夜夜 | 欧美不卡视频在线免费观看| av在线亚洲专区| 久久久久久久久久成人| 在线播放无遮挡| 少妇的逼水好多| 欧美日本视频| 伦精品一区二区三区| 熟妇人妻久久中文字幕3abv| 人妻制服诱惑在线中文字幕| 大香蕉久久网| 网址你懂的国产日韩在线| 在线播放无遮挡| 欧美+日韩+精品| 高清视频免费观看一区二区 | eeuss影院久久| 热99在线观看视频| 日韩欧美三级三区| 亚洲综合色惰| 99久久精品热视频| 嫩草影院精品99| 国产黄片视频在线免费观看| 欧美激情在线99| 一级毛片久久久久久久久女| 一个人免费在线观看电影| 少妇裸体淫交视频免费看高清| 一个人免费在线观看电影| 91在线精品国自产拍蜜月| 91精品伊人久久大香线蕉| 亚洲aⅴ乱码一区二区在线播放| 校园人妻丝袜中文字幕| 国产精品久久视频播放| 亚洲精品国产av成人精品| 一本久久精品| 欧美色视频一区免费| 两个人视频免费观看高清| 久久婷婷人人爽人人干人人爱| 99热网站在线观看| 汤姆久久久久久久影院中文字幕 | 中文字幕av在线有码专区| 久久精品夜色国产| 精品无人区乱码1区二区| 国产亚洲一区二区精品| 麻豆精品久久久久久蜜桃| 亚洲不卡免费看| 国产探花在线观看一区二区| 免费观看精品视频网站| 久久久精品大字幕| 如何舔出高潮| 欧美成人午夜免费资源| 久久久久久久午夜电影| 欧美一区二区亚洲| 久久99热这里只频精品6学生 | 亚洲av男天堂| 变态另类丝袜制服| 欧美一区二区亚洲| 免费搜索国产男女视频| 性插视频无遮挡在线免费观看| videossex国产| 又爽又黄无遮挡网站| 久久久精品大字幕| 日韩欧美在线乱码| 久久精品人妻少妇| 少妇人妻一区二区三区视频| 久久精品夜夜夜夜夜久久蜜豆| 精品久久久久久电影网 | 婷婷色av中文字幕| 男人的好看免费观看在线视频| 99热6这里只有精品| 特大巨黑吊av在线直播| 国产av一区在线观看免费| 天堂av国产一区二区熟女人妻| 久久久久免费精品人妻一区二区| 国产毛片a区久久久久| 边亲边吃奶的免费视频| 精品人妻偷拍中文字幕| 久久久久精品久久久久真实原创| 一边亲一边摸免费视频| 一个人免费在线观看电影| 91久久精品国产一区二区三区| 亚洲欧美日韩卡通动漫| 亚洲欧美清纯卡通| 国产乱来视频区| 国产大屁股一区二区在线视频| 高清视频免费观看一区二区 | 国产亚洲av嫩草精品影院| 九九爱精品视频在线观看| 中文在线观看免费www的网站| 亚洲天堂国产精品一区在线| 尤物成人国产欧美一区二区三区| 国产一区有黄有色的免费视频 | 一本久久精品| 亚洲四区av| 久久精品夜色国产| 亚洲欧美日韩卡通动漫| 欧美三级亚洲精品| 日本免费a在线| 久久精品夜色国产| 精品人妻一区二区三区麻豆| 91av网一区二区| 日本猛色少妇xxxxx猛交久久| 老司机影院成人| 久久久精品大字幕| 国产一级毛片七仙女欲春2| 国产精品,欧美在线| 国产在视频线精品| 亚洲精品乱码久久久v下载方式| 噜噜噜噜噜久久久久久91| 免费观看性生交大片5| 黑人高潮一二区| 丰满乱子伦码专区| 欧美激情在线99| 偷拍熟女少妇极品色| 亚洲国产色片| 天堂影院成人在线观看| 国产v大片淫在线免费观看| 18禁动态无遮挡网站| 中文字幕亚洲精品专区| 久久精品国产鲁丝片午夜精品| 丝袜喷水一区| 七月丁香在线播放| www.色视频.com| 乱人视频在线观看| 国产亚洲精品久久久com| 亚洲一区高清亚洲精品| 青春草视频在线免费观看| 亚洲欧美日韩东京热| 边亲边吃奶的免费视频| 国产女主播在线喷水免费视频网站 | 丰满少妇做爰视频| 日韩高清综合在线| 国产91av在线免费观看| 国产精品人妻久久久久久| 中文亚洲av片在线观看爽| 91aial.com中文字幕在线观看| 97超视频在线观看视频| 搡老妇女老女人老熟妇| 成人性生交大片免费视频hd| 国产免费男女视频| 看非洲黑人一级黄片| av国产久精品久网站免费入址| 99热精品在线国产| 久久国产乱子免费精品| 午夜福利视频1000在线观看| 一本久久精品| 久久久久久久久中文| 精华霜和精华液先用哪个| 国模一区二区三区四区视频| 国产日韩欧美在线精品| 一区二区三区高清视频在线| a级毛片免费高清观看在线播放| 精品人妻偷拍中文字幕| av播播在线观看一区| 天美传媒精品一区二区| 久久精品国产99精品国产亚洲性色| 可以在线观看毛片的网站| 天堂中文最新版在线下载 | 免费不卡的大黄色大毛片视频在线观看 | 我要搜黄色片| 欧美极品一区二区三区四区| 久久精品国产亚洲av涩爱| 国产极品天堂在线| 久久久久久久午夜电影| av女优亚洲男人天堂| 一卡2卡三卡四卡精品乱码亚洲| 高清视频免费观看一区二区 | 我要看日韩黄色一级片| 久久久色成人| 搡女人真爽免费视频火全软件| 综合色丁香网| 国产 一区 欧美 日韩| 天天一区二区日本电影三级| 亚洲精品,欧美精品| 免费观看性生交大片5| 免费黄网站久久成人精品| 不卡视频在线观看欧美| 久久久久久伊人网av| 网址你懂的国产日韩在线| 亚洲18禁久久av| 欧美另类亚洲清纯唯美| 91午夜精品亚洲一区二区三区| 天天躁日日操中文字幕| 国产激情偷乱视频一区二区| 国产一区二区在线观看日韩| 亚洲精品国产av成人精品| 国产精品乱码一区二三区的特点| 国产av不卡久久| 91精品国产九色| 亚洲精品日韩在线中文字幕| 99在线人妻在线中文字幕| 亚洲人与动物交配视频| 国产免费男女视频| 国产精品国产三级国产av玫瑰| 亚洲国产精品成人综合色| 久久人妻av系列| av福利片在线观看| 九色成人免费人妻av| 免费黄色在线免费观看| 国产高清有码在线观看视频| 国产精品电影一区二区三区| 国产av在哪里看| 精品人妻熟女av久视频| 日韩av在线免费看完整版不卡| 国产麻豆成人av免费视频| 成人毛片60女人毛片免费| 一级黄片播放器| 亚洲成人精品中文字幕电影| 国产在视频线在精品| 亚洲一级一片aⅴ在线观看| 国产大屁股一区二区在线视频| 国产精品人妻久久久影院| 成人欧美大片| 免费观看人在逋| 精品国产露脸久久av麻豆 | 亚洲中文字幕一区二区三区有码在线看| 国产在视频线精品| 欧美不卡视频在线免费观看| 成人午夜高清在线视频| 成人亚洲欧美一区二区av| 好男人在线观看高清免费视频| 欧美又色又爽又黄视频| 老司机福利观看| 精品久久久久久电影网 | 日本wwww免费看| 黄色欧美视频在线观看| 噜噜噜噜噜久久久久久91| 青春草视频在线免费观看| 国产精品三级大全| 国产午夜精品论理片| 成人无遮挡网站| 男女啪啪激烈高潮av片| 成年av动漫网址| 免费不卡的大黄色大毛片视频在线观看 | 99久国产av精品| 丰满人妻一区二区三区视频av| 少妇的逼水好多| 好男人视频免费观看在线| 中文字幕av成人在线电影| 欧美成人午夜免费资源| 麻豆国产97在线/欧美| 欧美激情久久久久久爽电影| 简卡轻食公司| 成人三级黄色视频| 成人二区视频| 色综合色国产| 免费av毛片视频| 午夜福利视频1000在线观看| 亚洲图色成人| 高清午夜精品一区二区三区| 激情 狠狠 欧美| 国产精品美女特级片免费视频播放器| 国产在线一区二区三区精 | 欧美人与善性xxx| 欧美日本亚洲视频在线播放| 国产男人的电影天堂91| 18禁裸乳无遮挡免费网站照片| 久久久久免费精品人妻一区二区| 18禁裸乳无遮挡免费网站照片| 日本与韩国留学比较| 国产又黄又爽又无遮挡在线| 免费观看人在逋| 寂寞人妻少妇视频99o| 色视频www国产| 国产精品国产三级国产av玫瑰| 欧美人与善性xxx| 好男人在线观看高清免费视频| 26uuu在线亚洲综合色| 亚洲av中文字字幕乱码综合| 男女视频在线观看网站免费| 久久久亚洲精品成人影院| 久久精品91蜜桃| 精品久久久久久久人妻蜜臀av| 日本免费在线观看一区| 精品午夜福利在线看| 久久久久九九精品影院| 级片在线观看| av天堂中文字幕网| 久久精品久久久久久噜噜老黄 | 亚洲成人久久爱视频| 亚洲在久久综合| 精品一区二区免费观看| 一本一本综合久久| 伊人久久精品亚洲午夜| 高清av免费在线| 国产久久久一区二区三区| 欧美精品一区二区大全| 久久99热这里只频精品6学生 | 欧美一级a爱片免费观看看| 久久久亚洲精品成人影院| 亚洲aⅴ乱码一区二区在线播放|