• <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
    丰满乱子伦码专区| 伊人久久国产一区二区| 小蜜桃在线观看免费完整版高清| 午夜福利网站1000一区二区三区| 丰满人妻一区二区三区视频av| 天美传媒精品一区二区| 99热这里只有是精品在线观看| 精品亚洲成国产av| 伦理电影大哥的女人| 免费人妻精品一区二区三区视频| 成人特级av手机在线观看| 久久久色成人| 秋霞伦理黄片| 国产精品一区二区性色av| 日本vs欧美在线观看视频 | 嫩草影院入口| 久久久久久久国产电影| 黄色日韩在线| 亚洲av.av天堂| 观看美女的网站| 国产精品国产三级国产专区5o| 精品一品国产午夜福利视频| 久久精品人妻少妇| 国产亚洲最大av| 国产片特级美女逼逼视频| 特大巨黑吊av在线直播| 18禁动态无遮挡网站| 又粗又硬又长又爽又黄的视频| 国产成人免费观看mmmm| 免费观看a级毛片全部| 国产精品人妻久久久影院| 国产精品女同一区二区软件| 久久97久久精品| 最后的刺客免费高清国语| 国产精品女同一区二区软件| 成人18禁高潮啪啪吃奶动态图 | 狂野欧美激情性xxxx在线观看| 欧美国产精品一级二级三级 | 国产成人a∨麻豆精品| 国产精品嫩草影院av在线观看| 性色av一级| 老熟女久久久| 午夜免费观看性视频| 日韩国内少妇激情av| 永久网站在线| 日本av免费视频播放| 国产精品一及| 青春草亚洲视频在线观看| 亚洲成人手机| 欧美日韩国产mv在线观看视频 | 日韩三级伦理在线观看| 精品久久久久久电影网| 国产视频首页在线观看| 热99国产精品久久久久久7| 成年人午夜在线观看视频| 日韩 亚洲 欧美在线| 午夜激情福利司机影院| 久久精品久久久久久久性| 99久久人妻综合| 久久久精品免费免费高清| 少妇高潮的动态图| 三级国产精品欧美在线观看| 五月伊人婷婷丁香| 啦啦啦在线观看免费高清www| 日韩成人av中文字幕在线观看| 性色avwww在线观看| 三级国产精品片| 最近最新中文字幕大全电影3| 亚洲性久久影院| 久久久久久久久久人人人人人人| 在线观看av片永久免费下载| 久热久热在线精品观看| 国产乱人偷精品视频| 亚洲欧美一区二区三区国产| 色婷婷久久久亚洲欧美| tube8黄色片| 国产精品99久久久久久久久| 国产欧美日韩精品一区二区| 精品少妇黑人巨大在线播放| 国产永久视频网站| 亚洲精品乱久久久久久| 国产精品女同一区二区软件| 国产精品久久久久久精品电影小说 | 久久影院123| av在线播放精品| 久久久久网色| 久久久久久人妻| 午夜免费观看性视频| 亚洲经典国产精华液单| 日韩av在线免费看完整版不卡| 亚洲国产av新网站| 久久ye,这里只有精品| 亚洲欧美清纯卡通| 久久精品久久久久久久性| 一级爰片在线观看| 性高湖久久久久久久久免费观看| 欧美亚洲 丝袜 人妻 在线| 亚洲av福利一区| 欧美成人a在线观看| 国产乱人偷精品视频| 国产精品久久久久久久电影| 丰满乱子伦码专区| 久久久久久人妻| 99视频精品全部免费 在线| 国产精品人妻久久久影院| 亚洲国产色片| 97超碰精品成人国产| 欧美精品一区二区免费开放| 王馨瑶露胸无遮挡在线观看| 美女视频免费永久观看网站| 80岁老熟妇乱子伦牲交| 免费看日本二区| 久久久成人免费电影| 久久ye,这里只有精品| 亚洲欧美中文字幕日韩二区| 免费少妇av软件| 大陆偷拍与自拍| 久久99热这里只有精品18| 啦啦啦中文免费视频观看日本| 妹子高潮喷水视频| 人妻 亚洲 视频| 成年美女黄网站色视频大全免费 | 免费观看a级毛片全部| 亚洲精品自拍成人| 午夜福利在线在线| 国产精品国产三级国产av玫瑰| 最后的刺客免费高清国语| 国产欧美日韩一区二区三区在线 | 国产免费一区二区三区四区乱码| 美女高潮的动态| 精品视频人人做人人爽| 色综合色国产| 亚洲精品,欧美精品| 日韩av免费高清视频| 我要看黄色一级片免费的| 蜜桃久久精品国产亚洲av| 亚洲性久久影院| 精品国产露脸久久av麻豆| 国产精品久久久久成人av| 精品人妻熟女av久视频| 色5月婷婷丁香| 在线观看av片永久免费下载| 免费看av在线观看网站| 日韩一区二区视频免费看| 精品一区二区三卡| 国产在视频线精品| 亚洲av综合色区一区| 日日摸夜夜添夜夜添av毛片| 1000部很黄的大片| 精品一区二区免费观看| 天堂中文最新版在线下载| 赤兔流量卡办理| 国产伦在线观看视频一区| 99国产精品免费福利视频| 男女啪啪激烈高潮av片| 伦理电影免费视频| a级毛色黄片| 国产伦理片在线播放av一区| 亚洲精品亚洲一区二区| 欧美日韩视频精品一区| 国产精品女同一区二区软件| 少妇 在线观看| 国产色婷婷99| 免费看光身美女| 在线观看av片永久免费下载| 国产精品国产三级国产av玫瑰| 一级毛片久久久久久久久女| 99久久精品国产国产毛片| 免费黄网站久久成人精品| 黄色配什么色好看| 精品国产乱码久久久久久小说| 日韩成人av中文字幕在线观看| 美女高潮的动态| 国产一区二区三区综合在线观看 | 网址你懂的国产日韩在线| 只有这里有精品99| 国产亚洲最大av| 久久99热6这里只有精品| 尾随美女入室| 中文欧美无线码| av在线播放精品| 欧美3d第一页| 久久鲁丝午夜福利片| 哪个播放器可以免费观看大片| 日韩强制内射视频| 一本一本综合久久| 国产在视频线精品| tube8黄色片| 国产成人精品福利久久| 精品人妻视频免费看| 国产精品三级大全| 成年人午夜在线观看视频| 免费黄网站久久成人精品| 亚洲最大成人中文| 一区在线观看完整版| 视频中文字幕在线观看| 国产精品伦人一区二区| 色5月婷婷丁香| 午夜福利高清视频| 久久久久视频综合| 天美传媒精品一区二区| 熟女人妻精品中文字幕| 在线观看人妻少妇| 亚洲国产高清在线一区二区三| 久久久久久久久大av| 日韩欧美精品免费久久| 亚洲av不卡在线观看| 国产在线免费精品| 最新中文字幕久久久久| 久久97久久精品| 女的被弄到高潮叫床怎么办| 九草在线视频观看| av免费观看日本| 亚洲,一卡二卡三卡| 亚洲精品456在线播放app| 一个人看视频在线观看www免费| 精品久久久久久久久亚洲| 国产免费视频播放在线视频| 中文字幕免费在线视频6| a级一级毛片免费在线观看| 五月天丁香电影| 亚洲人成网站高清观看| 久久这里有精品视频免费| 国产av精品麻豆| 国产亚洲5aaaaa淫片| 日韩大片免费观看网站| 亚洲欧美日韩卡通动漫| av国产精品久久久久影院| 日韩一区二区三区影片| 91精品一卡2卡3卡4卡| 欧美人与善性xxx| 99热国产这里只有精品6| 一级毛片久久久久久久久女| 日本黄色片子视频| 国产av一区二区精品久久 | 日韩av不卡免费在线播放| 最近最新中文字幕免费大全7| 精品久久久久久久末码| 97在线人人人人妻| 日本-黄色视频高清免费观看| xxx大片免费视频| 国产片特级美女逼逼视频| 80岁老熟妇乱子伦牲交| 国产极品天堂在线| 亚洲av电影在线观看一区二区三区| 天堂俺去俺来也www色官网| 不卡视频在线观看欧美| 看非洲黑人一级黄片| 精品久久久久久久末码| 又黄又爽又刺激的免费视频.| 肉色欧美久久久久久久蜜桃| 精品午夜福利在线看| av女优亚洲男人天堂| 尾随美女入室| 久久久久精品久久久久真实原创| 最后的刺客免费高清国语| 亚洲成人一二三区av| 最近手机中文字幕大全| 国产精品久久久久久久久免| 久久久久国产网址| 少妇丰满av| 亚洲av电影在线观看一区二区三区| 人妻系列 视频| 久久精品熟女亚洲av麻豆精品| 观看免费一级毛片| 人人妻人人看人人澡| 色网站视频免费| 久久国产乱子免费精品| 久久国产乱子免费精品| 久久影院123| 欧美zozozo另类| 男人狂女人下面高潮的视频| tube8黄色片| 黄色日韩在线| 亚洲一区二区三区欧美精品| 少妇人妻久久综合中文| 插阴视频在线观看视频| 国产乱人偷精品视频| 中文资源天堂在线| 亚洲性久久影院| 国产黄色免费在线视频| 亚洲自偷自拍三级| 嫩草影院入口| 亚洲精品乱码久久久v下载方式| 成人亚洲精品一区在线观看 | 视频中文字幕在线观看| 99热这里只有精品一区| 色5月婷婷丁香| 国产成人一区二区在线| 亚洲av.av天堂| 韩国av在线不卡| 国产精品一区二区三区四区免费观看| 777米奇影视久久| 精品人妻熟女av久视频| 性高湖久久久久久久久免费观看| av又黄又爽大尺度在线免费看| 51国产日韩欧美| 伦理电影大哥的女人| 高清av免费在线| 香蕉精品网在线| 日本-黄色视频高清免费观看| 国产成人免费观看mmmm| 交换朋友夫妻互换小说| 欧美少妇被猛烈插入视频| 女的被弄到高潮叫床怎么办| 日本免费在线观看一区| 国产高清国产精品国产三级 | 爱豆传媒免费全集在线观看| 一本一本综合久久| 80岁老熟妇乱子伦牲交| 亚洲三级黄色毛片| 人妻 亚洲 视频| 久久精品国产自在天天线| 春色校园在线视频观看| 不卡视频在线观看欧美| 一边亲一边摸免费视频| 晚上一个人看的免费电影| 日本黄大片高清| 嫩草影院新地址| 日本欧美视频一区| 欧美日韩在线观看h| 日韩免费高清中文字幕av| 免费看不卡的av| 久久99热这里只有精品18| 99久国产av精品国产电影| 欧美老熟妇乱子伦牲交| 97超视频在线观看视频| 极品教师在线视频| 人妻夜夜爽99麻豆av| 搡老乐熟女国产| 国产日韩欧美亚洲二区| 日本-黄色视频高清免费观看| 精品一区在线观看国产| 亚洲美女搞黄在线观看| 欧美日韩一区二区视频在线观看视频在线| 偷拍熟女少妇极品色| 最近最新中文字幕大全电影3| 精品国产露脸久久av麻豆| 日本爱情动作片www.在线观看| 久久久精品免费免费高清| 国内揄拍国产精品人妻在线| 婷婷色综合www| 少妇的逼水好多| 少妇人妻 视频| 简卡轻食公司| 九草在线视频观看| 日韩中文字幕视频在线看片 | 欧美精品一区二区免费开放| 亚洲国产最新在线播放| 日韩在线高清观看一区二区三区| 一区二区av电影网| 亚洲欧美一区二区三区国产| 一个人看的www免费观看视频| 多毛熟女@视频| 高清毛片免费看| 中文字幕av成人在线电影| 在线观看免费高清a一片| 人妻 亚洲 视频| 亚洲综合色惰| 少妇熟女欧美另类| av网站免费在线观看视频| 国产精品熟女久久久久浪| 天堂中文最新版在线下载| 免费观看的影片在线观看| 亚洲国产高清在线一区二区三| 最后的刺客免费高清国语| 精品国产露脸久久av麻豆| 国语对白做爰xxxⅹ性视频网站| 国产精品一区二区三区四区免费观看| 伦理电影免费视频| 亚洲国产毛片av蜜桃av| 五月开心婷婷网| 男女啪啪激烈高潮av片| 亚洲国产最新在线播放| 欧美日本视频| 亚洲自偷自拍三级| 久久99热6这里只有精品| 在线观看av片永久免费下载| av国产久精品久网站免费入址| 精品一品国产午夜福利视频| 国产精品国产三级国产av玫瑰| 欧美高清性xxxxhd video| 91午夜精品亚洲一区二区三区| 18禁裸乳无遮挡免费网站照片| 五月开心婷婷网| 免费大片18禁| 青春草国产在线视频| 在线看a的网站| 日本一二三区视频观看| 夫妻性生交免费视频一级片| 亚洲欧美日韩卡通动漫| 日韩,欧美,国产一区二区三区| 一级二级三级毛片免费看| 亚洲丝袜综合中文字幕| 少妇人妻久久综合中文| 久久久久久久大尺度免费视频| 亚洲精品自拍成人| av国产精品久久久久影院| 尾随美女入室| 激情五月婷婷亚洲| 又黄又爽又刺激的免费视频.| 精华霜和精华液先用哪个| 99久久综合免费| 久久久久久久久久成人| 国产精品欧美亚洲77777| 男人爽女人下面视频在线观看| 女人久久www免费人成看片| 成人特级av手机在线观看| h视频一区二区三区| 高清日韩中文字幕在线| 丝瓜视频免费看黄片| 黄片wwwwww| 中文字幕精品免费在线观看视频 | 欧美日韩综合久久久久久| 成人影院久久| 国产久久久一区二区三区| 亚洲精品视频女| 男人和女人高潮做爰伦理| 国产精品一区二区三区四区免费观看| av免费在线看不卡| 妹子高潮喷水视频| 国产免费福利视频在线观看| 91精品国产国语对白视频| 亚洲精华国产精华液的使用体验| 日韩中文字幕视频在线看片 | 亚洲欧美成人综合另类久久久| 午夜精品国产一区二区电影| 欧美日韩精品成人综合77777| 欧美3d第一页| 欧美xxⅹ黑人| 一区二区三区免费毛片| 精品少妇黑人巨大在线播放| 久久久国产一区二区| 精品人妻一区二区三区麻豆| 少妇人妻一区二区三区视频| 久久精品夜色国产| 身体一侧抽搐| av国产久精品久网站免费入址| 日本黄色日本黄色录像| 亚洲精品国产色婷婷电影| 免费人成在线观看视频色| 亚洲伊人久久精品综合| 国产大屁股一区二区在线视频| 一本色道久久久久久精品综合| 国产淫语在线视频| 超碰av人人做人人爽久久| 久久久久国产网址| 一级毛片我不卡| 两个人的视频大全免费| 美女高潮的动态| 视频区图区小说| 久久久久久人妻| 菩萨蛮人人尽说江南好唐韦庄| 丝袜喷水一区| 丰满乱子伦码专区| 黄色配什么色好看| 国内精品宾馆在线| 在线观看免费日韩欧美大片 | 亚洲欧美一区二区三区国产| 午夜日本视频在线| 男人爽女人下面视频在线观看| 亚洲第一av免费看| av国产久精品久网站免费入址| 国产精品国产av在线观看| 午夜激情久久久久久久| 97在线人人人人妻| 日本av免费视频播放| 简卡轻食公司| 秋霞伦理黄片| 菩萨蛮人人尽说江南好唐韦庄| 欧美成人精品欧美一级黄| 久久久午夜欧美精品| 国产精品熟女久久久久浪| 大陆偷拍与自拍| 一二三四中文在线观看免费高清| 日日摸夜夜添夜夜添av毛片| 国产精品欧美亚洲77777| 联通29元200g的流量卡| 亚州av有码| 香蕉精品网在线| 日本av免费视频播放| 国产视频内射| 精品人妻偷拍中文字幕| 亚洲精品日本国产第一区| 在线天堂最新版资源| 天天躁日日操中文字幕| 欧美激情极品国产一区二区三区 | 午夜老司机福利剧场| 亚洲经典国产精华液单| 欧美日韩视频高清一区二区三区二| 人妻一区二区av| 亚洲精品国产成人久久av| 少妇的逼水好多| 日本爱情动作片www.在线观看| 熟女电影av网| 麻豆精品久久久久久蜜桃| 亚洲国产av新网站| 久久久精品免费免费高清| 久久6这里有精品| av.在线天堂| 建设人人有责人人尽责人人享有的 | 国产精品一区二区性色av| 2022亚洲国产成人精品| 99久久中文字幕三级久久日本| 日韩一区二区三区影片| 男人舔奶头视频| 精品久久国产蜜桃| 日日撸夜夜添| 精品99又大又爽又粗少妇毛片| 自拍偷自拍亚洲精品老妇| 国产精品嫩草影院av在线观看| 美女脱内裤让男人舔精品视频| 亚洲成人手机| 九九在线视频观看精品| av不卡在线播放| 在线观看av片永久免费下载| 欧美精品亚洲一区二区| 在线观看免费日韩欧美大片 | 亚洲精品日韩av片在线观看| 日本欧美国产在线视频| 久久久久国产网址| 精品久久久久久久末码| 国产av码专区亚洲av| 精品一区二区三区视频在线| 欧美精品人与动牲交sv欧美| 国语对白做爰xxxⅹ性视频网站| 久久久久国产网址| 18禁在线无遮挡免费观看视频| 99视频精品全部免费 在线| 日韩免费高清中文字幕av| 高清黄色对白视频在线免费看 | 丝袜脚勾引网站| 两个人的视频大全免费| 下体分泌物呈黄色| 国产有黄有色有爽视频| 久久 成人 亚洲| 99re6热这里在线精品视频| 欧美国产精品一级二级三级 | 中文字幕免费在线视频6| av播播在线观看一区| 国语对白做爰xxxⅹ性视频网站| 五月开心婷婷网| 多毛熟女@视频| 国产日韩欧美在线精品| 国产国拍精品亚洲av在线观看| 久久精品熟女亚洲av麻豆精品| 国产精品免费大片| 久久这里有精品视频免费| 多毛熟女@视频| 中文资源天堂在线| 色哟哟·www| 国产精品国产av在线观看| 久久ye,这里只有精品| 国产伦精品一区二区三区四那| 国产精品爽爽va在线观看网站| kizo精华| 在线观看国产h片| 免费看光身美女| 国精品久久久久久国模美| 一个人看的www免费观看视频| 欧美xxxx黑人xx丫x性爽| 人体艺术视频欧美日本| 久久久色成人| 午夜日本视频在线| 男人狂女人下面高潮的视频| 精品国产三级普通话版| 精品亚洲成a人片在线观看 | 久久精品熟女亚洲av麻豆精品| 亚洲精品国产成人久久av| 国产精品99久久久久久久久| 在线观看免费高清a一片| 一级毛片黄色毛片免费观看视频| 国产淫片久久久久久久久| 永久免费av网站大全| 啦啦啦啦在线视频资源| 午夜激情久久久久久久| 成人亚洲精品一区在线观看 | 边亲边吃奶的免费视频| 亚洲美女搞黄在线观看| 亚洲av电影在线观看一区二区三区| 精品熟女少妇av免费看| 在线亚洲精品国产二区图片欧美 | 久久久久视频综合| 欧美激情极品国产一区二区三区 | 久久精品久久久久久噜噜老黄| 少妇的逼水好多| 婷婷色综合www| 精品久久久久久久久亚洲| av视频免费观看在线观看| 亚洲va在线va天堂va国产| 亚洲熟女精品中文字幕| 人妻 亚洲 视频| 插阴视频在线观看视频| 婷婷色综合大香蕉| 久久久精品免费免费高清| 一区二区三区乱码不卡18| 欧美日韩视频精品一区| 人妻 亚洲 视频| 最近手机中文字幕大全| 黄片wwwwww| 蜜桃在线观看..| 国产精品熟女久久久久浪| 亚洲av成人精品一二三区| 亚洲欧美成人精品一区二区| 日韩亚洲欧美综合| 免费播放大片免费观看视频在线观看| 欧美老熟妇乱子伦牲交| 在线 av 中文字幕| 少妇人妻久久综合中文| 精品少妇黑人巨大在线播放| 日韩 亚洲 欧美在线| 97超视频在线观看视频| 伦理电影大哥的女人| 成年免费大片在线观看| 有码 亚洲区|