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

    基于正交試驗的光電立靶光幕陣列結構參數(shù)優(yōu)化方法

    2017-12-01 09:35:58陳瑞倪晉平
    兵工學報 2017年11期
    關鍵詞:光幕紙板彈丸

    陳瑞, 倪晉平

    (西安工業(yè)大學 光電工程學院, 陜西 西安 710032)

    基于正交試驗的光電立靶光幕陣列結構參數(shù)優(yōu)化方法

    陳瑞, 倪晉平

    (西安工業(yè)大學 光電工程學院, 陜西 西安 710032)

    針對光電立靶光幕陣列結構參數(shù)標定誤差較大的問題,提出一種光幕陣列結構參數(shù)反演優(yōu)化的方法。以測量坐標與紙板靶坐標差值的平方根構建目標函數(shù),借助正交試驗在允許誤差范圍內產生多組不同參數(shù)量值的組合,設計了光幕陣列結構參數(shù)優(yōu)化方法。以雙N形天幕靶為例,在Matlab中進行仿真,優(yōu)化后的光幕陣列結構參數(shù)更接近給定真值,且結構參數(shù)的優(yōu)化程度與標定的初始值和誤差范圍相關。用實彈射擊試驗對設計的方法進行驗證,結果表明,優(yōu)化后的彈丸飛行參數(shù)測量結果更接近實際真值,設計的方法有效可行。該優(yōu)化方法為提高光電立靶測量設備的測量精度提供了一種新思路。

    兵器科學與技術; 光幕陣列; 多目標優(yōu)化; 反演; 結構參數(shù)

    0 引言

    彈丸飛行參數(shù)是衡量身管武器性能和武器效能評估的重要指標,主要包括彈丸空間坐標、飛行速度及速度方位角和俯仰角等[1-2]?,F(xiàn)有的彈丸飛行參數(shù)測量方法主要有:紙板靶法[3-4]、聲學法[5]、電荷耦合器件法[6-7]、光幕陣列法[8-9]。其中,光幕陣列法以其適用范圍廣、可靠性高、成本低廉等優(yōu)點在靶場外彈道測量領域得到越來越廣泛的應用。若干個光幕在空間以特定的幾何結構排列,彈丸依次穿過各個光幕面,通過記錄彈丸到達每個光幕面的時刻,配合光幕陣列結構參數(shù)信息,從而解算出彈丸的飛行參數(shù)。光幕陣列空間幾何結構多種多樣,雙光幕陣列以平行結構排列,可以通過定距測時的方法測量彈丸飛行速度;四光幕陣列以雙N形結構排列,可以測量垂直入射彈丸的速度和坐標[10];六光幕陣列以雙V形或雙N形結構排列,可以測量任意角度飛行的彈丸坐標、速度及速度方位角和俯仰角[11-12]等。各種類型光幕陣列的空間結構一般通過光幕之間的夾角進行描述,它們在光幕陣列測量設備成型后通過標定的方法得到,但由于形成的光幕在空間是肉眼不可見的,無法直接測量,常常導致標定的角度誤差較大,從而影響測量精度。本文以雙N形六光幕陣列天幕立靶為例提出一種光幕陣列結構參數(shù)反演優(yōu)化方法,將彈丸實際著靶坐標作為已知條件,以彈丸測量坐標與實際坐標差值的均方根構建目標函數(shù),采用正交試驗在多組不同結構參數(shù)量值組合下尋求目標函數(shù)的最小值,最終實現(xiàn)光幕陣列結構參數(shù)優(yōu)化。

    1 雙N形光幕陣列測量原理

    雙N形六光幕陣列天幕立靶的一個靶體可在空間形成一個N形光幕,兩個靶體(分別稱為啟動靶和停止靶)沿彈道線方向布置,可在空間形成雙N形排列的探測光幕。由于彈丸穿過光幕面時探測電路輸出的彈丸過幕時刻為彈丸到達每個幕面固定位置的時刻,在雙N形六光幕陣列測量原理中空間光幕可用理想平面表示[13-14]。圖1給出了雙N形六光幕陣列的測量原理,圖2是空間光幕在坐標系平面內的投影。

    圖1 雙N形六光幕陣列示意圖Fig.1 Structure of double-N shaped six-light-screen-array

    圖2 空間光幕在坐標系平面投影Fig.2 Projection view of light-screen structure

    測量坐標系如圖1所示,s1為起始靶與停止靶靶距、h為高度差,統(tǒng)稱為布靶參數(shù)。圖2中實線光幕表示圖1中光幕在坐標系內的真實投影,虛線光幕表示圖1中光幕在坐標系內投影的截面,則光幕陣列空間結構參數(shù)可表示為豎直夾角α=[α1,α2,α3,α4]和水平夾角β=[β1,β2],記彈丸到達每個光幕的時刻為t1~t6,將6個光幕面G1~G6用空間平面方程表示,在靶距較短的情況下彈丸在光幕陣列內所受的空氣阻力和重力作用可忽略不計,則此過程中可假定彈丸做勻速直線運動,彈道線與空間光幕的交點以矩陣形式[14]表示為

    M·X=N,

    (1)

    通過解算得到的彈丸在光幕G1內的空間坐標和飛行速度矢量為

    X=M-1·N.

    (2)

    從(2)式可見,彈丸的空間坐標測量結果與光幕陣列結構參數(shù)密切相關。現(xiàn)有的光幕結構參數(shù)獲取方法通常需要對天幕立靶進行拆卸,再由人工標定得到光幕之間的夾角,不僅標定過程繁瑣,且得到的光幕陣列結構參數(shù)往往也存在較大誤差,若有一種方法可對標定的結果進行優(yōu)化,將有效的提高光幕陣列設備的測量精度。

    2 優(yōu)化方法及目標函數(shù)構建

    由于光幕陣列在測量過程中對彈丸的飛行軌跡不產生任何影響,在其后垂直于彈道方向再放置一個紙板靶,則可通過彈孔直接測量彈丸在紙板靶上的坐標(視為準真值),以雙N形六光幕天幕立靶為例表示上述過程,如圖3所示。

    圖3 雙N形光幕陣列優(yōu)化方法模型Fig.3 Optimization model of double-N shaped six-light-screen-array

    圖3中S2為停止靶到紙板靶的距離,(2)式中的測量結果可從G1平面換算至紙板靶平面得到光幕陣列測量的紙板靶處彈丸坐標[15],理論上在紙板靶處光幕陣列測量的坐標與紙板靶坐標相等,但由于光幕陣列結構參數(shù)和其他誤差因素的影響,兩組坐標必然存在一定偏差,該偏差表征了光幕陣列測量結果與準真值的接近程度,在其他誤差因素不變的情況下,該偏差可以用來表征光幕陣列結構參數(shù)與結構參數(shù)真值的接近程度。

    在圖3所示的模型中,將雙N形光幕陣列和紙板靶測定的兩組坐標結果分別記為(z,y)和(z*,y*),構建兩個目標函數(shù)δy和δz分別表征光幕陣列測量結果與紙板靶測量結果縱坐標和橫坐標的偏差,即

    (3)

    式中:n是參與優(yōu)化的彈丸總數(shù);i是彈丸序號。

    由于光幕陣列測量的彈丸坐標結果是光幕陣列幾何結構的函數(shù),其值包含了光幕陣列的結構參數(shù)信息,且在雙N形六光幕陣列天幕立靶測量模型中,空間光幕夾角一般在設備成型后固定不變。因此在其他誤差因素不變的情況下,測量公式中選取的光幕陣列結構參數(shù)越接近實際真值則兩組坐標測量結果的偏差越小,反之則越大,構建的目標函數(shù)δy和δz可分別度量結構參數(shù)在豎直和水平方向的準確度。

    3 基于正交試驗設計的光幕陣列參數(shù)優(yōu)化方法

    將雙N形光幕陣列結構參數(shù)作為影響試驗結果的因素,基于正交試驗產生多組不同量值的結構參數(shù)組合,能充分反應全面試驗的要求。依據(jù)測量原理,坐標測量結果中縱坐標y僅受光幕豎直面夾角α的影響,而橫坐標z測量結果主要受光幕水平面夾角β的影響,因此本文依據(jù)(3)式設計兩組連續(xù)試驗先后對α=[α1,α2,α3,α4]和β=[β1,β2]進行優(yōu)化。

    首先將[α1,α2,α3,α4]作為影響試驗結果δy的4個因素即待反演優(yōu)化參數(shù),將標定的結果[α1o,α2o,α3o,α4o]作為初始值,在初始值附近選擇標定誤差±ε作為實際光幕夾角可能出現(xiàn)的范圍生成5水平,對應的因素水平表如表1所示。

    選取正交表L25(56)對表1中的4個參數(shù)設計4因素5水平試驗,過程中由于β對立靶縱坐標測量結果沒有影響,其值代入標定的初始值[β1o,β2o]不變,在α標定的初始值[α1o,α2o,α3o,α4o]附近產生25種不同形式的正交組合,不同的參數(shù)組合下對應的彈丸縱坐標y隨之變化,進而引起目標函數(shù)δy的變化,假設試驗中目標函數(shù)的最小值出現(xiàn)在第j組,則此時光幕夾角α與實際情況最接近,其對應的光幕陣列豎直方向夾角[α1j,α2j,α3j,α4j]即為α優(yōu)化后的結果。

    表1 豎直面內夾角因素水平表

    之后對水平面內夾角β進行優(yōu)化,由于[β1,β2]決定彈丸著靶橫坐標z的測量結果,此時目標函數(shù)變?yōu)棣膠,待反演優(yōu)化的參數(shù)只有兩個,根據(jù)正交規(guī)則β1的每一個水平和β2的每一個水平都需產生組合,以標定結果[β1o,β2o]作為初始值,在初始值附近選擇標定誤差±ε作為實際光幕夾角可能出現(xiàn)的范圍生成5水平,對應的因素水平表如表2所示。

    表2 水平內夾角因素水平表

    與α的優(yōu)化過程類似將表2中的2個5水平因素填入正交表L25(56)的前2列產生25種正交組合,取α為之前優(yōu)化得到的結果[α1j,α2j,α3j,α4j]不變,則不同的參數(shù)β組合下目標函數(shù)δz取值發(fā)生顯著變化,假設試驗中目標函數(shù)的最小值出現(xiàn)在第k次試驗,則對應的[β1k,β2k]即為水平夾角優(yōu)化后的結果。至此通過兩次連續(xù)試驗可對雙N形六光幕天幕立靶的6個光幕空間結構參數(shù)完成反演優(yōu)化。

    4 仿真分析及優(yōu)化

    在Matlab中建立雙N形六光幕陣列模型,假定光幕陣列結構參數(shù)的真值為α1=α2=α3=α4=25°,β1=β2=24°,并在模型中精確給定,在上述真值附近±0.6°隨機生成一組數(shù)字作為光幕陣列結構參數(shù)的標定結果:α1o=25.4°,α2o=25.5°,α3o=24.9°,α4o=24.5°;β1o=23.5°,β2o=23.6°;構建500 mm×500 mm靶面并生成n=15發(fā)彈丸穿過靶面,生成彈丸的飛行參數(shù)為彈速v=(720±30)m/s,速度俯仰角γ=0°±4°,速度方向角θ=0°±4°.

    用第3節(jié)設計的方法,選擇ε=0.6°作為各夾角真值可能出現(xiàn)的范圍,對生成的標定結果進行仿真優(yōu)化,得到所有組合下目標函數(shù)的最小值δy,min=1.859 4,δz,min=8.562 6對應的結構參數(shù)優(yōu)化結果見表3.

    表3 雙N形光幕陣列仿真結果

    從表3可見,仿真優(yōu)化后各結構參數(shù)的取值更接近于給定真值。

    為研究不同參數(shù)初始值和參數(shù)范圍ε對本文方法的影響,取ε=0.3°和ε=0.9°在真值附近重新隨機生成兩組天幕立靶初始結構參數(shù)并用第3節(jié)設計的方法對其進行仿真優(yōu)化,兩次仿真得到所有組合下目標函數(shù)的最小值分別為δy,min=1.502 3,δz,min=4.511 4和δy,min=2.134 7,δz,min=12.396 1,對應的結構參數(shù)初始值和優(yōu)化后的結果分別見表4和表5所示。

    表4 雙N形光幕陣列仿真結果(ε=0.3°)

    表5 雙N形光幕陣列仿真結果(ε=0.9°)

    從表4和表5中可以看出,本文方法對光幕陣列結構參數(shù)的優(yōu)化程度與光幕陣列結構參數(shù)的標定方法相關,標定精度高的方法初始值相對準確、參數(shù)誤差范圍(步長)ε較小,優(yōu)化結果也越好。另外還需要注意的是,理論上為了減小隨機誤差,此方法在實際工程應用時選取的彈丸總數(shù)不宜過少,為了盡可能地使試驗簡化,在射擊過程中可使參與優(yōu)化的彈丸在靶面內均勻分布,在減小射擊次數(shù)的同時反映空間不同位置處的光幕陣列結構信息,從而保證反演優(yōu)化結果的精度。

    5 實際測試驗證

    依據(jù)圖3所示的模型,用雙N形光幕陣列天幕立靶工程樣機進行試驗,試驗現(xiàn)場如圖4所示。

    圖4 試驗現(xiàn)場Fig.4 Photo of testing site

    采用木板貼坐標紙的方式構成500 mm×500 mm紙板靶后方放置收彈箱,圖4中激光水平儀、投線儀等用于保證立靶與紙板靶間相對關系,布靶過程中起始靶與停止靶高度差h=0 mm、距離s1=4 724 mm,停止靶與紙板靶距離s2=1 400 mm. 用氣槍進行20次射擊使彈丸均勻分布在整個靶面,分別記錄每發(fā)彈丸穿過光幕陣列的過靶時刻序列和彈丸在紙板靶上的彈孔坐標,用本文方法將彈丸過靶時刻序列代入不同試驗序號下對應的光幕陣列結構,解算天幕立靶著靶坐標,并換算至紙板靶所在的平面與紙板靶坐標進行比對,計算得到對應的目標函數(shù)值。

    以所選雙N形光幕陣列標定值作為初始值,精度取ε=0.6°,優(yōu)化前后的結構參數(shù)如表6所示。

    表6 雙N形光幕陣列實彈試驗結果

    圖5 優(yōu)化前后測量結果比對Fig.5 Comparison of measured results of structure parameters before and after optimization

    為驗證本文優(yōu)化結果的正確性,將優(yōu)化后的雙N型光幕陣列天幕立靶在室外重新進行15次隨機射擊,分別用優(yōu)化前后的兩組結構參數(shù)對著靶坐標進行測量,對優(yōu)化前后結構參數(shù)下測量結果與紙板靶的偏差進行統(tǒng)計,結果如圖5所示。

    從圖5可以看出,兩組試驗結果與實際彈孔坐標的測量誤差在橫坐標、縱坐標方向分布一致,均在0 mm的上下波動,但在初始結構參數(shù)下測量結果偏差的誤差絕對值及波動較大,優(yōu)化后的結構參數(shù)下測量結果偏差的誤差絕對值及波動較小。試驗驗證了本文方法可以實現(xiàn)雙N形天幕立靶光幕陣列結構參數(shù)的優(yōu)化,優(yōu)化后的光幕陣列測量結果與彈丸實際著靶坐標相比誤差更小。

    6 結論

    光幕陣列方法測量彈丸飛行參數(shù),結構參數(shù)的準確性直接影響測量精度,以天幕立靶為例,本文提出了一種光電立靶光幕陣列結構參數(shù)優(yōu)化方法,用于提高其準確性。得到主要結論如下:

    1) 將彈丸的彈孔坐標作為準真值,以測量坐標與彈孔坐標差值平方根最小為原則構建目標函數(shù),借助正交試驗表設計反演方法,實現(xiàn)了光幕陣列結構參數(shù)優(yōu)化。

    2) 在Matlab中進行仿真,結果顯示用設計的方法優(yōu)化后的雙N形光幕陣列結構參數(shù)更接近給定真值,且結構參數(shù)的優(yōu)化程度與結構參數(shù)標定的初始值和誤差范圍相關。

    3) 用實彈射擊試驗對設計的方法進行驗證,結果表明,優(yōu)化后的彈丸飛行參數(shù)測量結果更接近實際真值,設計的方法有效可行。

    本文方法可減小光電立靶測量設備結構參數(shù)的標定誤差,為提高該類設備的測量精度提供了一種新思路。

    References)

    [1] 王昌明. 實用彈道學[M]. 北京:兵器工業(yè)出版社, 1994.

    WANG Chang-ming. Utility ballistics[M]. Beijing: Publishing House of Ordnance Industry, 1994. (in Chinese)

    [2] Stepanov E, Vavilov S A. The main problem of external ballistics[J]. Computers & Mathematics with Applications, 1997, 33(5):95-101.

    [3] David C, Smith P E. Line of sight beyond line of sight advanced technology demonstrator[C]∥40th Annual Armament Systems: Gun-Ammunition-Rockets-Missiles Conference & Exhibition. New Orleans, LA, US: Department of Defense, 2005.

    [4] Corriveau D. System modeling of a 40 mm automatic grenade launcher[C]∥42nd Gun and Missile Systems Conference & Exhibition. Charlotte, NC, US:NDIA, 2007.

    [5] Canistraro H A, Jordan E H. Projectile-impact-location determination: an acoustic triangulation method[J]. Measurement Science and Technology, 1996, 7(12):1755-1760.

    [6] 劉澤慶, 張玉榮, 趙建新, 等. 基于數(shù)字攝影測量的靶場高速攝影測速方法[J]. 彈道學報, 2015, 27(4):47-51.

    LIU Ze-qing, ZHANG Yu-rong, ZHAO Jian-xin, et al. High-speed photography velocity measurement in range based on digital photogrammetry[J]. Journal of Ballistics, 2015, 27(4):47-51.(in Chinese)

    [7] Ma W H, Dong T, Tian H, et al. Line-scan CCD camera calibration in 2D coordinate measurement[J]. Optik, 2014, 125(17): 4795-4798.

    [8] Gauthier Jr L R., Drewry Jr D G, Brunner L. Optical sensor and method for detecting projectile impact location and velocity vector: US, US6931166[P]. 2005-08-26.

    [9] Kalonia R C, Chhachhia D P, Bajpai P P, et al. Multiple laser-based high-speed digital shadowgraphy system for small caliber projectile-target interaction studies[J]. Optical Engineering, 2014, 53(3):034104.

    [10] HPI.B572 optical target system[EB/OL]. [2015-03-08]. http:∥www.hpi-gmbh.com/redesign/products_10.htm.

    [11] Li H S, Gao J C, Wang Z M. Object location fire precision test technology by using intersecting photoelectric detection target[J]. Optik, 2014, 125(3): 1325-1329.

    [12] 田會, 倪晉平, 焦明星. 拋物線彈道彈丸飛行參數(shù)測量模型與精度分析[J]. 儀器儀表學報, 2016, 37(1):67-74

    TIAN Hui, NI Jin-ping, JIAO Ming-xing. Projectile flying parameter measurement model of parabolic trajectory and precision analysis[J]. Chinese Journal of Scientific Instrument, 2016, 37(1):67-74.(in Chinese)

    [13] 張斌, 趙冬娥, 劉吉, 等. 基于重心法的激光光幕彈丸過把信號特征點提取算法[J]. 火炮發(fā)射與控制學報, 2015, 36(3):68-71.

    ZHANG Bin, ZHAO Dong-e, LIU Ji, et al. Algorithm of extracting the moment when projectile passes through a laser screen using center of gravity method[J]. Journal of Gun Launch & Control, 2015, 36(3):68-71.(in Chinese)

    [14] 倪晉平, 田會. 用于光幕測試的時刻信息提取方法研究[J]. 光學技術, 2008, 34(1):141-148.

    NI Jin-ping, TIAN Hui. A study on method for acquiring moment on which a projectile going through a light screens[J]. Optical Technique, 2008, 34(1):141-148.(in Chinese)

    [15] 陳瑞, 倪晉平, 侯颯. 三種著靶坐標自動測量設備精度評價方法[J]. 測試技術學報, 2017, 31(1):70-75.

    CHEN Rui,NI Jin-ping, HOU Sa. Three accuracy assessment methods of automatic projectile impact coordinates measurement system[J]. Journal of Test and Measurement Technology, 2017, 31(1):70-75.(in Chinese)

    OptimizationMethodofLight-screen-arrayStructureParametersofPhotoelectricTargetBasedonOrthogonalTest

    CHEN Rui, NI Jin-ping

    (School of Optoelectronic Engineering, Xi’an Technological University, Xi’an 710032, Shaanxi, China)

    An optimization method for inversion of light-screen-array structure parameters is presented for their large calibrated errors. An objective function is established by using the square root of difference between measuring coordinates and actual coordinates, many groups of structure parameters are generated within the allowable error range based on orthogonal test, and the optimization method of light-screen-array structure parameters is designed. The simulated results of double-N shaped sky screen target in Matlab shows that the optimized results are affected by the initial value and the allowable error range, and the ball cartridge test results show that the proposed method is feasible and effective. The optimization method provides a new idea for improving the measuring accuracy of photoelectric target measuring equipment.

    ordnance science and technology; light-screen-array; multi-objective optimization; inversion; structure parameter

    TJ012.3+7

    A

    1000-1093(2017)11-2234-06

    10.3969/j.issn.1000-1093.2017.11.021

    2017-04-12

    國家自然科學基金項目(61471289); 陜西省自然科學基礎研究計劃項目(2015JM5189); 陜西省教育廳重點實驗室科研計劃項目(14JS035)

    陳瑞(1987—), 男, 博士研究生。 E-mail: chenrui_xatu@163.com

    倪晉平(1965—), 男, 教授, 博士生導師。 E-mail: jpni@263.net

    猜你喜歡
    光幕紙板彈丸
    超高速撞擊下球形彈丸破碎特性仿真研究
    紙板填數(shù)
    神秘的『彈丸』
    紙板俄羅斯方塊拼圖
    針對復雜光線環(huán)境而生,相伴九城巡展 ACME艾可美|常焦抗光幕
    非理想狀態(tài)下斜幕面對四光幕陣列精度靶測量誤差的影響*
    2013年芬蘭紙與紙板凈出口量高達947萬噸
    激光復合光幕靶測速系統(tǒng)設計
    彈丸對預開孔混凝土靶體侵徹的實驗研究
    一種電梯光幕保護系統(tǒng)的設計
    電腦與電信(2014年3期)2014-03-16 03:38:44
    av在线观看视频网站免费| 亚洲国产看品久久| 老女人水多毛片| 高清视频免费观看一区二区| 如日韩欧美国产精品一区二区三区| 国产成人精品婷婷| 18禁在线无遮挡免费观看视频| 欧美成人午夜精品| 巨乳人妻的诱惑在线观看| 99久久人妻综合| 中文欧美无线码| 精品熟女少妇av免费看| 大陆偷拍与自拍| 成人毛片60女人毛片免费| 国产精品成人在线| 视频区图区小说| 精品久久蜜臀av无| 久久热在线av| 欧美日韩视频高清一区二区三区二| 女性生殖器流出的白浆| 免费高清在线观看日韩| 人人妻人人澡人人爽人人夜夜| 国产在线免费精品| 狠狠婷婷综合久久久久久88av| 深夜精品福利| 好男人视频免费观看在线| 在线观看三级黄色| 久久ye,这里只有精品| 国产精品久久久久久av不卡| 国产成人a∨麻豆精品| 青青草视频在线视频观看| 免费观看无遮挡的男女| 咕卡用的链子| 97人妻天天添夜夜摸| 超碰97精品在线观看| 久久毛片免费看一区二区三区| 久久99一区二区三区| 成人国产av品久久久| 亚洲国产最新在线播放| 亚洲伊人色综图| 黄色一级大片看看| 日本-黄色视频高清免费观看| 香蕉国产在线看| av免费在线看不卡| 久久ye,这里只有精品| 亚洲人与动物交配视频| 热99久久久久精品小说推荐| 少妇精品久久久久久久| 人妻一区二区av| 2022亚洲国产成人精品| 人人妻人人添人人爽欧美一区卜| 欧美丝袜亚洲另类| 国产在线视频一区二区| 国产 一区精品| 免费大片18禁| 18禁动态无遮挡网站| 亚洲,欧美,日韩| 两性夫妻黄色片 | 日日啪夜夜爽| 精品熟女少妇av免费看| 黑人猛操日本美女一级片| 水蜜桃什么品种好| 中国美白少妇内射xxxbb| 国产黄色免费在线视频| 日本猛色少妇xxxxx猛交久久| 在现免费观看毛片| 考比视频在线观看| 国产 一区精品| 最近中文字幕高清免费大全6| 黄网站色视频无遮挡免费观看| 国产精品人妻久久久久久| 男女啪啪激烈高潮av片| 日韩成人av中文字幕在线观看| 欧美3d第一页| 黄色配什么色好看| 欧美人与性动交α欧美精品济南到 | 婷婷色麻豆天堂久久| 国产片内射在线| 麻豆乱淫一区二区| 欧美日本中文国产一区发布| 国产一区二区三区av在线| 久久热在线av| 成人二区视频| 久久久久精品人妻al黑| 久久这里只有精品19| 久久精品国产综合久久久 | 国产男女内射视频| 韩国高清视频一区二区三区| 黄色一级大片看看| 国产成人91sexporn| 咕卡用的链子| 80岁老熟妇乱子伦牲交| 国产黄频视频在线观看| 国产成人免费无遮挡视频| 久久久久久人人人人人| 全区人妻精品视频| 国产欧美另类精品又又久久亚洲欧美| 欧美xxxx性猛交bbbb| 久久久精品免费免费高清| 日韩人妻精品一区2区三区| 丰满乱子伦码专区| 性色av一级| 91精品三级在线观看| 哪个播放器可以免费观看大片| 一区二区三区四区激情视频| 久久久a久久爽久久v久久| 制服丝袜香蕉在线| 国产精品欧美亚洲77777| 免费黄色在线免费观看| 欧美日本中文国产一区发布| 免费不卡的大黄色大毛片视频在线观看| 2021少妇久久久久久久久久久| 国产日韩欧美在线精品| 日韩av不卡免费在线播放| 久久ye,这里只有精品| 啦啦啦啦在线视频资源| 亚洲人成77777在线视频| 女人精品久久久久毛片| 秋霞在线观看毛片| 女人久久www免费人成看片| 久久久久精品久久久久真实原创| 欧美国产精品va在线观看不卡| 精品一区二区三卡| 亚洲色图 男人天堂 中文字幕 | 亚洲人成网站在线观看播放| 男人操女人黄网站| 亚洲国产精品一区三区| 91精品国产国语对白视频| 香蕉国产在线看| 我的女老师完整版在线观看| 国产男人的电影天堂91| 综合色丁香网| 下体分泌物呈黄色| 精品亚洲成a人片在线观看| 国产欧美亚洲国产| 亚洲精品乱码久久久久久按摩| 啦啦啦在线观看免费高清www| 久久 成人 亚洲| 久久久久久久国产电影| 国产福利在线免费观看视频| 欧美少妇被猛烈插入视频| 丰满饥渴人妻一区二区三| 大香蕉久久网| 欧美精品国产亚洲| 韩国av在线不卡| 国产在线视频一区二区| 十分钟在线观看高清视频www| 曰老女人黄片| 大香蕉久久网| 97在线视频观看| 超碰97精品在线观看| 久久久国产精品麻豆| 一级毛片 在线播放| 国产精品久久久久久精品电影小说| 免费看不卡的av| 欧美日韩精品成人综合77777| 丝袜美足系列| 久久精品aⅴ一区二区三区四区 | 亚洲熟女精品中文字幕| 亚洲av免费高清在线观看| 午夜免费男女啪啪视频观看| 精品久久久久久电影网| 亚洲精品日本国产第一区| 欧美激情 高清一区二区三区| 26uuu在线亚洲综合色| 最近手机中文字幕大全| 99国产精品免费福利视频| 精品久久久久久电影网| 黄网站色视频无遮挡免费观看| 亚洲中文av在线| 少妇的丰满在线观看| 熟女人妻精品中文字幕| 国产免费现黄频在线看| 91久久精品国产一区二区三区| 如何舔出高潮| 久久久久久久久久久久大奶| 毛片一级片免费看久久久久| 久久毛片免费看一区二区三区| 欧美bdsm另类| 老司机亚洲免费影院| 午夜免费男女啪啪视频观看| 最近最新中文字幕免费大全7| 人人妻人人爽人人添夜夜欢视频| 亚洲一级一片aⅴ在线观看| 少妇人妻精品综合一区二区| 久久国产精品大桥未久av| 九色成人免费人妻av| 超碰97精品在线观看| 九九在线视频观看精品| 久久综合国产亚洲精品| 国产免费又黄又爽又色| 91国产中文字幕| 国产免费一区二区三区四区乱码| 国产又色又爽无遮挡免| 少妇人妻 视频| 在线 av 中文字幕| 成人亚洲精品一区在线观看| 国产精品三级大全| 久久久久久久国产电影| 日产精品乱码卡一卡2卡三| 国产精品久久久久久av不卡| 色5月婷婷丁香| 极品人妻少妇av视频| 丝袜脚勾引网站| 久久鲁丝午夜福利片| 一区在线观看完整版| 春色校园在线视频观看| 九草在线视频观看| 一边亲一边摸免费视频| 中文字幕亚洲精品专区| 最新的欧美精品一区二区| 嫩草影院入口| 在线精品无人区一区二区三| 国产激情久久老熟女| 黄色 视频免费看| 天堂中文最新版在线下载| 青青草视频在线视频观看| 91成人精品电影| 亚洲精品,欧美精品| 18禁在线无遮挡免费观看视频| av在线app专区| 国产免费一区二区三区四区乱码| 美国免费a级毛片| 欧美亚洲日本最大视频资源| 国产日韩欧美在线精品| 插逼视频在线观看| 日产精品乱码卡一卡2卡三| 久久人人爽人人爽人人片va| 日韩免费高清中文字幕av| 欧美3d第一页| 黄色怎么调成土黄色| 一本大道久久a久久精品| 精品99又大又爽又粗少妇毛片| 热99久久久久精品小说推荐| 久久人妻熟女aⅴ| 自拍欧美九色日韩亚洲蝌蚪91| 成年人午夜在线观看视频| 男女无遮挡免费网站观看| 母亲3免费完整高清在线观看 | 夫妻性生交免费视频一级片| 国产精品国产三级国产av玫瑰| 99久国产av精品国产电影| 黑人巨大精品欧美一区二区蜜桃 | 王馨瑶露胸无遮挡在线观看| 男女午夜视频在线观看 | www日本在线高清视频| 久久这里有精品视频免费| 国产精品偷伦视频观看了| av国产精品久久久久影院| 中文字幕制服av| 最近中文字幕2019免费版| 黑人高潮一二区| 一级毛片我不卡| 老女人水多毛片| 熟女电影av网| 国产极品天堂在线| √禁漫天堂资源中文www| 亚洲,一卡二卡三卡| 26uuu在线亚洲综合色| 十分钟在线观看高清视频www| 亚洲色图综合在线观看| 国产成人欧美| 免费人成在线观看视频色| 日韩欧美精品免费久久| 女的被弄到高潮叫床怎么办| 国产精品国产av在线观看| 日本爱情动作片www.在线观看| 色婷婷久久久亚洲欧美| 乱人伦中国视频| 久久婷婷青草| 丰满饥渴人妻一区二区三| 男人添女人高潮全过程视频| 国产免费福利视频在线观看| 七月丁香在线播放| 亚洲国产色片| 成人综合一区亚洲| 亚洲av日韩在线播放| 国产麻豆69| 18禁在线无遮挡免费观看视频| 国产精品麻豆人妻色哟哟久久| 天天操日日干夜夜撸| 夫妻午夜视频| 亚洲av综合色区一区| 久久久国产精品麻豆| av又黄又爽大尺度在线免费看| 水蜜桃什么品种好| 欧美老熟妇乱子伦牲交| 国产一区二区在线观看日韩| 毛片一级片免费看久久久久| 精品亚洲成国产av| 国产成人精品久久久久久| 国产亚洲午夜精品一区二区久久| 99热全是精品| www.av在线官网国产| 亚洲av男天堂| 在线 av 中文字幕| 18禁国产床啪视频网站| av免费观看日本| 王馨瑶露胸无遮挡在线观看| 18禁动态无遮挡网站| 久久精品国产综合久久久 | 亚洲人成网站在线观看播放| 这个男人来自地球电影免费观看 | 午夜免费观看性视频| 午夜日本视频在线| 精品少妇久久久久久888优播| 自线自在国产av| 美女脱内裤让男人舔精品视频| 2018国产大陆天天弄谢| a级毛色黄片| 飞空精品影院首页| 夫妻性生交免费视频一级片| 免费高清在线观看视频在线观看| 黑人猛操日本美女一级片| 国产 精品1| 久久精品国产鲁丝片午夜精品| 国产熟女欧美一区二区| 精品久久久久久电影网| 亚洲婷婷狠狠爱综合网| 亚洲图色成人| 日本免费在线观看一区| 国产男人的电影天堂91| 人妻系列 视频| 伦理电影大哥的女人| 久久久久久伊人网av| 最近2019中文字幕mv第一页| 伦精品一区二区三区| 最近最新中文字幕免费大全7| 成年av动漫网址| 91aial.com中文字幕在线观看| 少妇 在线观看| 成人毛片60女人毛片免费| xxxhd国产人妻xxx| 国产色婷婷99| 亚洲欧美精品自产自拍| 人妻一区二区av| 欧美bdsm另类| av电影中文网址| 成人午夜精彩视频在线观看| 国产一区二区三区av在线| 人妻人人澡人人爽人人| 欧美变态另类bdsm刘玥| 日本午夜av视频| 亚洲精品国产av成人精品| 午夜福利,免费看| 欧美国产精品一级二级三级| 80岁老熟妇乱子伦牲交| 男女国产视频网站| 午夜福利视频在线观看免费| 你懂的网址亚洲精品在线观看| 欧美人与善性xxx| 男人操女人黄网站| 性色avwww在线观看| 人妻 亚洲 视频| 一级,二级,三级黄色视频| 日本爱情动作片www.在线观看| 精品亚洲乱码少妇综合久久| 久久久亚洲精品成人影院| 亚洲欧美日韩卡通动漫| 天天影视国产精品| 国产黄频视频在线观看| 亚洲欧美清纯卡通| 午夜福利,免费看| 又粗又硬又长又爽又黄的视频| 国产亚洲精品第一综合不卡 | 亚洲精品中文字幕在线视频| 国产 一区精品| 成人二区视频| 欧美bdsm另类| 亚洲av在线观看美女高潮| 久久精品久久久久久噜噜老黄| 女的被弄到高潮叫床怎么办| 亚洲欧洲国产日韩| 十八禁高潮呻吟视频| 草草在线视频免费看| 亚洲国产成人一精品久久久| 精品国产一区二区久久| 高清在线视频一区二区三区| 亚洲精品色激情综合| 国产色爽女视频免费观看| 一本—道久久a久久精品蜜桃钙片| 日日啪夜夜爽| 一级毛片 在线播放| 涩涩av久久男人的天堂| 国产一区亚洲一区在线观看| 久久婷婷青草| 久久午夜综合久久蜜桃| 国产精品人妻久久久久久| 国产视频首页在线观看| 国产又色又爽无遮挡免| 视频中文字幕在线观看| 欧美激情极品国产一区二区三区 | 日本黄大片高清| 在线 av 中文字幕| 亚洲精品乱码久久久久久按摩| 久久久久国产精品人妻一区二区| 亚洲精品,欧美精品| 啦啦啦啦在线视频资源| 少妇被粗大猛烈的视频| 亚洲图色成人| 亚洲成国产人片在线观看| 国产精品久久久av美女十八| 视频区图区小说| 一二三四在线观看免费中文在 | 亚洲精品久久久久久婷婷小说| 精品国产国语对白av| 男人添女人高潮全过程视频| 在线观看免费日韩欧美大片| 一区二区三区乱码不卡18| 岛国毛片在线播放| 青春草亚洲视频在线观看| 欧美成人午夜免费资源| 免费高清在线观看视频在线观看| 亚洲综合精品二区| 18+在线观看网站| 日韩精品有码人妻一区| 九色亚洲精品在线播放| 国产探花极品一区二区| 国产国拍精品亚洲av在线观看| 人妻系列 视频| 老司机影院毛片| 亚洲精品国产色婷婷电影| 91成人精品电影| 亚洲天堂av无毛| 9热在线视频观看99| 波野结衣二区三区在线| 久久午夜综合久久蜜桃| 一级爰片在线观看| 久久久久精品久久久久真实原创| 国产精品久久久久成人av| 久热这里只有精品99| 波多野结衣一区麻豆| av女优亚洲男人天堂| 国产伦理片在线播放av一区| 777米奇影视久久| 欧美日韩亚洲高清精品| 美女大奶头黄色视频| 王馨瑶露胸无遮挡在线观看| 欧美人与善性xxx| 又粗又硬又长又爽又黄的视频| 中文字幕制服av| 一级,二级,三级黄色视频| 麻豆精品久久久久久蜜桃| 成人黄色视频免费在线看| 天天躁夜夜躁狠狠躁躁| 国产一区亚洲一区在线观看| 久久毛片免费看一区二区三区| 免费大片18禁| 久久女婷五月综合色啪小说| www日本在线高清视频| 久久久亚洲精品成人影院| 激情视频va一区二区三区| 亚洲精品乱码久久久久久按摩| 中文乱码字字幕精品一区二区三区| 日韩人妻精品一区2区三区| 丰满饥渴人妻一区二区三| 精品国产乱码久久久久久小说| 久久精品国产a三级三级三级| 亚洲成国产人片在线观看| 在线观看免费视频网站a站| 亚洲国产精品一区三区| 国产成人精品在线电影| 老司机影院成人| 免费高清在线观看日韩| 中文乱码字字幕精品一区二区三区| 婷婷色麻豆天堂久久| 丝袜在线中文字幕| 亚洲一级一片aⅴ在线观看| 日韩欧美一区视频在线观看| 日韩成人av中文字幕在线观看| 波野结衣二区三区在线| 久久久久精品性色| 香蕉国产在线看| 老熟女久久久| 色哟哟·www| 高清视频免费观看一区二区| 王馨瑶露胸无遮挡在线观看| 国产女主播在线喷水免费视频网站| 99热全是精品| 草草在线视频免费看| 精品一区二区三卡| 观看美女的网站| 久久久久久久久久人人人人人人| 97人妻天天添夜夜摸| 午夜av观看不卡| av黄色大香蕉| 男女高潮啪啪啪动态图| 深夜精品福利| 国产综合精华液| 尾随美女入室| 亚洲成国产人片在线观看| 国产精品蜜桃在线观看| 边亲边吃奶的免费视频| 午夜福利视频精品| 岛国毛片在线播放| 伦理电影免费视频| 久久这里有精品视频免费| av又黄又爽大尺度在线免费看| 亚洲少妇的诱惑av| 日本vs欧美在线观看视频| 久久精品国产亚洲av天美| 69精品国产乱码久久久| 亚洲av综合色区一区| 欧美日韩视频精品一区| 久久久精品区二区三区| 这个男人来自地球电影免费观看 | 啦啦啦视频在线资源免费观看| 久久这里有精品视频免费| 99热网站在线观看| 亚洲av免费高清在线观看| 免费看光身美女| 一级a做视频免费观看| 国产在线一区二区三区精| 国产精品一国产av| 久久人人爽人人片av| 天天操日日干夜夜撸| 又黄又粗又硬又大视频| videosex国产| 校园人妻丝袜中文字幕| 黑丝袜美女国产一区| 99久久综合免费| 天堂中文最新版在线下载| 少妇 在线观看| 日本色播在线视频| av黄色大香蕉| 多毛熟女@视频| 精品少妇黑人巨大在线播放| 一级毛片我不卡| 美女大奶头黄色视频| 另类精品久久| 极品人妻少妇av视频| 久久狼人影院| 久久人人爽av亚洲精品天堂| 老司机影院毛片| 亚洲国产精品成人久久小说| 91国产中文字幕| 久久精品国产a三级三级三级| 一级a做视频免费观看| av有码第一页| 一二三四在线观看免费中文在 | 免费黄色在线免费观看| 极品人妻少妇av视频| 免费看不卡的av| 精品国产国语对白av| 亚洲av国产av综合av卡| 国产日韩欧美亚洲二区| 久久午夜综合久久蜜桃| 国产精品久久久久久久电影| 热re99久久精品国产66热6| 日韩大片免费观看网站| 99国产精品免费福利视频| 亚洲四区av| 69精品国产乱码久久久| 精品一品国产午夜福利视频| 国产精品久久久久久av不卡| 亚洲av中文av极速乱| 高清在线视频一区二区三区| 欧美性感艳星| 大话2 男鬼变身卡| 18禁动态无遮挡网站| 亚洲欧洲日产国产| 国产男女超爽视频在线观看| 日韩一区二区三区影片| 七月丁香在线播放| 在线精品无人区一区二区三| 国产片特级美女逼逼视频| 最新中文字幕久久久久| 午夜免费男女啪啪视频观看| 美国免费a级毛片| 国产精品无大码| 国产又爽黄色视频| 啦啦啦视频在线资源免费观看| 一级片免费观看大全| 内地一区二区视频在线| 成人影院久久| 免费观看在线日韩| 中文字幕亚洲精品专区| 一级毛片电影观看| 国产av码专区亚洲av| 99热这里只有是精品在线观看| 黑人高潮一二区| 岛国毛片在线播放| 亚洲五月色婷婷综合| 中国三级夫妇交换| 久久久久精品人妻al黑| 国产成人aa在线观看| 国产成人免费观看mmmm| 欧美+日韩+精品| 曰老女人黄片| 亚洲综合色惰| 久久精品国产亚洲av涩爱| a 毛片基地| 菩萨蛮人人尽说江南好唐韦庄| 久久久欧美国产精品| 欧美人与性动交α欧美软件 | 国产又爽黄色视频| 久久韩国三级中文字幕| av在线观看视频网站免费| 国产又爽黄色视频| 亚洲成人手机| 两性夫妻黄色片 | 成人国产av品久久久| 三上悠亚av全集在线观看| 黑人猛操日本美女一级片| 香蕉丝袜av| 欧美日韩亚洲高清精品| 伦理电影免费视频| 亚洲欧美一区二区三区黑人 | 国产成人精品在线电影| 久久精品久久久久久久性| 免费久久久久久久精品成人欧美视频 | 国产精品偷伦视频观看了| 男人添女人高潮全过程视频| 精品卡一卡二卡四卡免费| 亚洲熟女精品中文字幕|