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

    磁強(qiáng)計陣列測量一致性校正

    2022-05-11 07:51:28雷怡琴孫兆龍洪澤宏劉琪鐘佩娜
    關(guān)鍵詞:磁強(qiáng)計系統(tǒng)誤差磁通

    雷怡琴, 孫兆龍, 洪澤宏, 劉琪, 鐘佩娜

    (1.海軍工程大學(xué) 電氣工程學(xué)院, 湖北 武漢 430033; 2.中國船舶重工集團(tuán)公司第七一五研究所, 浙江 杭州 310012)

    磁強(qiáng)計自身具有一定的物理局限性,在許多較大的工程應(yīng)用場景中單個磁強(qiáng)計并不能滿足測量要求[1],尤其在需要大平面陣列進(jìn)行測量的情況中,單個磁強(qiáng)計受到技術(shù)原理以及加工工藝水平限制,短期內(nèi)并不能依托于提高性能實現(xiàn)工程要求[2-5],因此對磁強(qiáng)計陣列技術(shù)的準(zhǔn)確性研究尤為重要。磁強(qiáng)計陣列技術(shù)指的是在待測磁場空間或者待測磁源周圍布放多個磁強(qiáng)計,同時進(jìn)行磁場測量[6]。影響陣列系統(tǒng)測量準(zhǔn)確性的主要原因有磁強(qiáng)計自身的非正交、刻度因子、零偏等系統(tǒng)誤差因素[7-8],以及在磁強(qiáng)計陣列的搭建安裝過程中由于安裝工藝水平的限制,使各磁強(qiáng)計之間軸向姿態(tài)不一致引起的旋轉(zhuǎn)誤差因素。

    研究表明,當(dāng)磁強(qiáng)計的坐標(biāo)軸有1°偏差時,測量誤差將達(dá)到0.87%[9-10]。目前,針對磁強(qiáng)計陣列校正的研究較少,朱興樂[11]采用DE算法對磁強(qiáng)計的水平誤差角進(jìn)行校正,但是該校正方法的前提是各磁強(qiáng)計垂直分量軸需保持一致,屬于二維平面校正。磁強(qiáng)計三維校正方法可借鑒磁梯度儀的誤差校正,現(xiàn)有磁梯度儀的非對準(zhǔn)誤差的主要解決方法有2種:①基于機(jī)械結(jié)構(gòu)的磁強(qiáng)計姿態(tài)調(diào)整;②通過軟件的手段進(jìn)行校正。基于機(jī)械結(jié)構(gòu)的非對準(zhǔn)誤差校正方法不僅對機(jī)械結(jié)構(gòu)要求高,且陣列中磁強(qiáng)計數(shù)量較多時,校正不易實現(xiàn)。關(guān)于軟件校對方法,于振濤等對梯度儀系統(tǒng)誤差采用兩步法進(jìn)行補(bǔ)償校準(zhǔn),張光等對張量系統(tǒng)進(jìn)行線性化校準(zhǔn),具有一定的校正效果,但是兩者均忽略了實際測量輸出的高階小量,影響了測量精度[12-15]。

    本文提出了一種針對磁強(qiáng)計陣列的一致性校正模型,在不忽略高階小量的基礎(chǔ)上構(gòu)建九參數(shù)磁通門磁強(qiáng)計系統(tǒng)誤差校正模型,對構(gòu)成測磁陣列的磁強(qiáng)計的非正交、刻度因子、零偏誤差進(jìn)行校正得到理想正交三分量磁場測量值。在此基礎(chǔ)上,以陣列中一個磁強(qiáng)計的軸向為參考,構(gòu)建三維九參數(shù)歐拉旋轉(zhuǎn)矩陣校正模型,對其他磁強(qiáng)計進(jìn)行俯仰、橫傾、方位3種姿態(tài)轉(zhuǎn)換,實現(xiàn)了對磁強(qiáng)計陣列的軸向一致性校正。

    1 磁通門磁強(qiáng)計系統(tǒng)誤差校正

    磁通門磁強(qiáng)計陣列校正是對磁強(qiáng)計之間坐標(biāo)軸的一致性進(jìn)行校正,單個磁通門磁強(qiáng)計的測量精度直接影響了磁強(qiáng)計陣列的測量精度[16-18]。為了保證陣列一致性校正的精確性首先需要對磁強(qiáng)計的系統(tǒng)誤差進(jìn)行校正,保證每個磁強(qiáng)計均是理想正交磁強(qiáng)計。

    建立磁場測量實際空間坐標(biāo)O-X1Y1Z1與理想坐標(biāo)O-XYZ,如圖1所示。設(shè)磁場空間坐標(biāo)系的原點與理想坐標(biāo)系的原點重合,再令實際坐標(biāo)系垂向軸OZ1與理想坐標(biāo)系的垂向軸OZ同軸,Z1OY1與ZOY共面。設(shè)OY與OY1之間的夾角為ψ,OX1與XOY的夾角為φ;OX1在XOY面上的投影與OX的夾角為θ。

    圖1 磁通門三軸示意圖

    三軸磁通門的測量誤差來源主要有非正交誤差、零位誤差及刻度因子誤差,這些誤差的具體表達(dá)式如下所示:

    1) 非正交誤差

    在不考慮其他誤差的情況下對圖1進(jìn)行矢量分解,得到實際測量磁場值B1與理想測量磁場值B之間的關(guān)系為

    (1)

    2) 刻度因子誤差

    刻度因子誤差是不同通道的采集電路無法做到完全一致導(dǎo)致的,設(shè)刻度因子誤差矩陣為K,表示為

    (2)

    3) 零偏誤差

    零偏誤差是由磁強(qiáng)計內(nèi)部剩磁引起的,表達(dá)式為

    (3)

    綜上,可構(gòu)建實際輸出B1與理想輸出B之間的關(guān)系為

    (4)

    對(4)式進(jìn)行轉(zhuǎn)換,可將三軸磁通門磁強(qiáng)計存在誤差干擾時測量得到的校正值當(dāng)做測量磁場的真實值,得到輸出校正值為

    (5)

    磁強(qiáng)計的誤差參數(shù)校正模型為

    (6)

    利用最小二乘法可對kx,ky,kz,θ,φ,ψ,bx,by,bz進(jìn)行求解,則可完成對磁強(qiáng)計測量值的補(bǔ)償校正。

    2 測磁陣列一致性校正

    勻強(qiáng)靜磁場環(huán)境中,測磁陣列中的各磁強(qiáng)計經(jīng)過系統(tǒng)誤差校正后,不同位置上的磁強(qiáng)計測量輸出值均收斂于環(huán)境磁場值。由于安裝具有一定偏差,經(jīng)過校正后的磁強(qiáng)計各理想正交軸朝向錯亂,導(dǎo)致不同位置上的磁強(qiáng)計存在旋轉(zhuǎn)非對準(zhǔn)誤差。如圖2所示,本文選定測磁陣列中磁強(qiáng)計1的坐標(biāo)軸作為參考坐標(biāo)系,將測磁系統(tǒng)中的其他磁強(qiáng)計轉(zhuǎn)化到磁強(qiáng)計1的坐標(biāo)軸下。

    圖2 磁強(qiáng)計陣列布放模擬

    三維一致性校正在空間上進(jìn)行,空間中任意2個正交坐標(biāo)系可通過包含3個非對準(zhǔn)角度的旋轉(zhuǎn)矩陣進(jìn)行轉(zhuǎn)換。在均勻場環(huán)境中,將一個坐標(biāo)軸轉(zhuǎn)換到參考坐標(biāo)系中,雖然有不同的旋轉(zhuǎn)方法,但是旋轉(zhuǎn)后得到的磁場測量值是固定的,故通過建立歐拉旋轉(zhuǎn)矩陣可對磁強(qiáng)計進(jìn)行一致性校正。

    2.1 橫傾轉(zhuǎn)換模型

    本文分別從x,y,z3個坐標(biāo)系建立轉(zhuǎn)換模型,對測磁陣列進(jìn)行一致性校正。設(shè)待校正磁強(qiáng)計的實際坐標(biāo)軸為O-X2Y2Z2,實際磁場測量值為(Bx,By,Bz),本文采用模擬旋轉(zhuǎn)磁強(qiáng)計坐標(biāo)軸的方法,對磁強(qiáng)計進(jìn)行一致性校正。建立橫傾正交轉(zhuǎn)換模型,如圖3所示,令待校正磁強(qiáng)計與基準(zhǔn)磁強(qiáng)計的x軸重合,再繞x軸旋轉(zhuǎn),使By與XOY共面。

    圖3 橫傾轉(zhuǎn)換模型

    (7)

    2.2 俯仰轉(zhuǎn)換模型

    圖4 俯仰轉(zhuǎn)換模型

    (8)

    2.3 方位轉(zhuǎn)換模型

    圖5 方位轉(zhuǎn)換模型

    (9)

    2.4 三維一致性轉(zhuǎn)換模型

    磁強(qiáng)計安裝完畢后的角度固定,因此旋轉(zhuǎn)角度λ,θ,δ為固定值,不會隨著時間變化。磁強(qiáng)計三軸變換的歐拉旋轉(zhuǎn)矩陣表達(dá)式為

    (10)

    式中:Tλ為待校正磁強(qiáng)計的橫傾旋轉(zhuǎn)矩陣;Tθ為待校正磁強(qiáng)計的俯仰旋轉(zhuǎn)矩陣;Tδ為待校正磁強(qiáng)計的俯仰旋轉(zhuǎn)矩陣。綜合上述3次旋轉(zhuǎn)變換,空間上任意姿態(tài)的2個正交坐標(biāo)系,都可以通過橫傾、俯仰、方位的三維變換轉(zhuǎn)化到同一坐標(biāo)系下。由(10)式可以看出,矩陣的乘積順序與繞軸旋轉(zhuǎn)的順序相匹配,本文采取橫傾、俯仰、方位的旋轉(zhuǎn)順序進(jìn)行一致性校正,待校正磁強(qiáng)計與參考磁強(qiáng)計的轉(zhuǎn)換關(guān)系表示為

    (11)

    (12)

    向量形式的表達(dá)方式如(13)式所示。

    (13)

    (12)式可以看出,每一種旋轉(zhuǎn)次序?qū)?yīng)一種旋轉(zhuǎn)矩陣,λ,θ,δ互相匹配。將(13)式用增廣矩陣表示成3個線性方程組

    (14)

    (15)

    由(11)式可以看出,每一次旋轉(zhuǎn)校正是在上一次旋轉(zhuǎn)的基礎(chǔ)上進(jìn)行的,當(dāng)旋轉(zhuǎn)順序固定時,可計算出一致性校正的唯一角度,聯(lián)立方程(12)~(15)可以求解得到待校正磁強(qiáng)計的非對準(zhǔn)誤差唯一角度,此時磁強(qiáng)計陣列的非線性誤差校正完成。

    3 仿真分析

    3.1 磁通門磁強(qiáng)計校正仿真分析

    利用MATLAB對上述校正方法進(jìn)行性能分析,設(shè)待校正磁強(qiáng)計所處的環(huán)境磁場大小為48 750 nT。為了得到單個磁通門磁強(qiáng)計的校正參數(shù),需要得到磁強(qiáng)計在空間上任意方位的磁場測量數(shù)據(jù)。把磁強(qiáng)計陣列中的每個磁強(qiáng)計置于如圖6所示的三維旋轉(zhuǎn)臺上,從初始位置開始繞三維空間軸進(jìn)行旋轉(zhuǎn),進(jìn)而獲得全空間方位上的磁場測量矢量,在x軸方向上進(jìn)行以60°為采樣間隔的旋轉(zhuǎn),y軸方向上進(jìn)行以90°為采樣間隔的旋轉(zhuǎn),z軸方向上進(jìn)行以60°為采樣間隔的旋轉(zhuǎn),并記錄數(shù)據(jù)。

    圖6 三維轉(zhuǎn)臺

    為了方便觀察,對磁場測量數(shù)據(jù)進(jìn)行插值,圖7中紅色圓點為補(bǔ)償前的實測點,藍(lán)色點為插值擬合點。

    圖7 空間磁場測量值三維示意圖

    由圖7可以看出,在未對磁強(qiáng)計進(jìn)行補(bǔ)償校正的前提下,磁強(qiáng)計的實際空間輸出為橢球形磁場,通過對方程(6)的求解,得到每個磁強(qiáng)計校正的參數(shù)矩陣,把校正參數(shù)代入到實際測量結(jié)果中,得到磁強(qiáng)計的理想測量值。

    對4個磁通門磁強(qiáng)計分別進(jìn)行仿真校正,計算4個磁強(qiáng)計的誤差矩陣,分別為磁強(qiáng)計系統(tǒng)誤差矩陣M,以及零偏噪聲誤差矩陣b計算結(jié)果如表1所示。

    表1 磁強(qiáng)計補(bǔ)償校正系數(shù)矩陣

    把表1中的校正參數(shù)系數(shù)矩陣代入到磁強(qiáng)計的測量結(jié)果中,利用最小二乘法對磁場測量數(shù)據(jù)進(jìn)行后處理,對陣列中所有磁強(qiáng)計的測量誤差進(jìn)行補(bǔ)償,消除單個磁強(qiáng)計的橢球化影響。補(bǔ)償校正后的磁強(qiáng)計矢量三維空間磁場測量結(jié)果如圖8所示,圖中藍(lán)色球體是以空間磁場圓點為圓心,背景磁場為半徑的理想球體,橙色點為補(bǔ)償后的磁場測量值。

    圖8 空間磁場校正結(jié)果三維示意圖

    由圖8可以看出,對磁強(qiáng)計進(jìn)行校正之后,測量值均落在了理想球面上,磁場測量橢球化得到了很好消除。通過第一節(jié)中的校正方法對陣列中的磁通門磁強(qiáng)計測量結(jié)果進(jìn)行校正,反映到陣列磁場測量總場值中的校正效果如圖9所示。結(jié)果表明,進(jìn)行非正交補(bǔ)償后的磁場測量值收斂于空間磁場值,達(dá)到了理想效果。

    圖9 磁場總場值校正結(jié)果

    利用磁強(qiáng)計校正前后的輸出相對方差值對磁場測量校正效果進(jìn)行表征,以此來定量反映磁強(qiáng)計的系統(tǒng)誤差以及非正交誤差的校正效果。由于地磁環(huán)境時刻變化,因此磁強(qiáng)計進(jìn)行測量時總會存在波動,為了更明顯地展示校正效果,本文利用每次磁場測量最大值進(jìn)行精度表征,輸出相對方差值計算公式為

    (16)

    式中:N為測量姿態(tài)數(shù);|Bi|max為第i個磁強(qiáng)計最大測量絕對值;Br為磁場測量標(biāo)準(zhǔn)值。

    表2 測量相對方差值

    校正前后每個磁強(qiáng)計的相對方差誤差對比結(jié)果列于表2。由實驗結(jié)果可以看出,磁強(qiáng)計校正前的相對方差誤差在百級以上,校正后磁強(qiáng)計的均方根誤差在10 nT以內(nèi),校正效果明顯。

    3.2 磁強(qiáng)計陣列一致性校正仿真分析

    對磁強(qiáng)計測量結(jié)果進(jìn)行校正是對磁強(qiáng)計陣列進(jìn)行一致性校正的前提,在對磁強(qiáng)計測量結(jié)果校正的基礎(chǔ)上,對磁強(qiáng)計陣列進(jìn)行一致性校正。當(dāng)磁強(qiáng)計陣列位于勻強(qiáng)磁場的環(huán)境中時,磁強(qiáng)計陣列中每個磁強(qiáng)計的磁場測量軌跡應(yīng)該重合,利用上述磁強(qiáng)計姿態(tài)校正方法對各磁強(qiáng)計進(jìn)行校正,將1號磁強(qiáng)計作為參考對其他磁強(qiáng)計依次進(jìn)行一致性校正,計算得到各磁強(qiáng)計的一致性校正系數(shù)如表3所示。

    表3 一致性校正系數(shù)

    將表3中的一致性校正系數(shù)代入每個磁強(qiáng)計的磁場測量結(jié)果中,得到各個磁強(qiáng)計相對誤差校正結(jié)果如圖10所示。

    圖10 補(bǔ)償校正前后磁場測量結(jié)果對比

    由仿真結(jié)果可以看出,本文提出的磁強(qiáng)計陣列一致性校正方法可以有效補(bǔ)償磁強(qiáng)計陣列輸出,補(bǔ)償后磁強(qiáng)計陣列磁場測量輸出誤差小于±20nT,校正效果良好。

    4 實驗結(jié)果分析

    把參考磁強(qiáng)計與待校正磁強(qiáng)計置于轉(zhuǎn)臺上進(jìn)行實驗驗證,保證各磁強(qiáng)計到轉(zhuǎn)臺中心位置的距離相同,使測磁系統(tǒng)繞磁中心旋轉(zhuǎn)一周。本文主要針對由4個磁強(qiáng)計構(gòu)成的磁強(qiáng)計陣列進(jìn)行一致性校正,搭建一致性校正實驗平臺如圖11所示,由于實驗條件有限,將磁強(qiáng)計兩兩一組置于無磁轉(zhuǎn)臺上,通過無磁轉(zhuǎn)臺改變測磁系統(tǒng)姿態(tài)。磁強(qiáng)計間距離為0.5m,工作溫度范圍是-40~70℃,環(huán)境溫度為12℃。實驗在外界磁場穩(wěn)定且干擾較小的環(huán)境中進(jìn)行,根據(jù)世界地磁模型WMM2020查詢得知,實驗地點的地磁場總量Br約為48 500nT,以此作為磁強(qiáng)計陣列測量的標(biāo)準(zhǔn)磁場輸出值。

    圖11 無磁轉(zhuǎn)臺

    令無磁轉(zhuǎn)臺繞z軸進(jìn)行旋轉(zhuǎn),以此改變磁強(qiáng)計的空間姿態(tài),保持磁強(qiáng)計的初始姿態(tài)水平,此時x軸與y軸的旋轉(zhuǎn)角度均為0°,轉(zhuǎn)臺每轉(zhuǎn)10°記錄1次磁場測量值,共旋轉(zhuǎn)1周,每個磁強(qiáng)計有36組姿態(tài)數(shù)據(jù),該磁強(qiáng)計陣列共有4×36組姿態(tài)數(shù)據(jù)。

    4.1 磁通門磁強(qiáng)計非正交補(bǔ)償實驗結(jié)果分析

    以Br作為參考輸出值,利用上述方法對磁強(qiáng)計陣列中的4個磁強(qiáng)計進(jìn)行校正,得到磁強(qiáng)計的測量輸出校正結(jié)果如圖12所示。用各個磁強(qiáng)計的測量輸出標(biāo)準(zhǔn)差表征磁強(qiáng)計的補(bǔ)償校正性能,校正結(jié)果如表4所示。

    圖12 測量輸出校正結(jié)果

    表4 測量輸出標(biāo)準(zhǔn)差

    由圖12中對磁強(qiáng)計補(bǔ)償校正的實驗結(jié)果可知,經(jīng)過補(bǔ)償校正后磁強(qiáng)計的磁場測量輸出值穩(wěn)定在48 500nT左右,與標(biāo)準(zhǔn)輸出值擬合程度較高。由表4可知,經(jīng)過補(bǔ)償校正后的各磁強(qiáng)計測量輸出的誤差均方根值在10nT以內(nèi),磁強(qiáng)計的輸出達(dá)到了較高的精度。

    4.2 磁強(qiáng)計陣列一致性校正結(jié)果分析

    以磁強(qiáng)計陣列中選定的磁強(qiáng)計坐標(biāo)軸為參考,對待校正磁強(qiáng)計進(jìn)行一致性校正,采用最小二乘法進(jìn)行匹配估計求解,收斂速度如圖13所示,可以看出該方法收斂速度快,計算精度高。

    圖13 計算精度

    圖14以磁場值的三維形式直觀地顯示了校正結(jié)果。當(dāng)各磁強(qiáng)計坐標(biāo)軸之間具有較高的一致性時,在三維圖中表現(xiàn)為磁場測量結(jié)果重合度較高。其中圖14a)為三維轉(zhuǎn)臺上實際磁場測量值,圖14b)為對磁強(qiáng)計進(jìn)行系統(tǒng)誤差校正后的磁場測量值,圖14c)為經(jīng)過一致性校正后的磁場測量值。實驗結(jié)果表明,若對磁強(qiáng)計陣列只進(jìn)行系統(tǒng)誤差校正而不進(jìn)行一致性校正時,磁強(qiáng)計陣列測量不僅重合度未得到改善,反而使非重合度增加,進(jìn)行一致性校正之后,各個磁強(qiáng)計在均勻的磁場環(huán)境中測量得到的磁場值重合度達(dá)到了理想水平。

    圖14 三維轉(zhuǎn)臺實驗結(jié)果

    為了更精確地顯示一致性校正效果,表5列出了一致性校正前后陣列中的各待校正磁強(qiáng)計與標(biāo)準(zhǔn)參考磁強(qiáng)計之間的角度差。一致性校正前待校正磁強(qiáng)計與參考磁強(qiáng)計各坐標(biāo)軸之間的角度差在1°左右,甚至高達(dá)2°,這對弱磁測量結(jié)果會產(chǎn)生很大影響。由實驗結(jié)果可以看出,經(jīng)過一致性校正之后各坐標(biāo)軸之間的誤差減小到了0.01°以下,大大提高了磁強(qiáng)計陣列的測量精度。

    表5 待校正磁強(qiáng)計與標(biāo)準(zhǔn)參考磁強(qiáng)計的角度差

    5 結(jié) 論

    與單個磁強(qiáng)計相比,磁強(qiáng)計陣列能夠直接獲得多個位置的磁場信息,同時獲得不同位置的磁場變化規(guī)律,具有更高的可靠性。

    隨著磁強(qiáng)計陣列的廣泛使用,本文針對其存在的姿態(tài)不一致問題進(jìn)行了建模校正,首先針對磁通門磁強(qiáng)計構(gòu)建了非正交、刻度因子系統(tǒng)誤差校正模型。在不忽略高階小量的情況下,得到九參數(shù)測量補(bǔ)償校正系數(shù)矩陣。在對磁強(qiáng)計進(jìn)行測量補(bǔ)償校正的基礎(chǔ)上,通過歐拉旋轉(zhuǎn)矩陣搭建三維磁強(qiáng)計陣列一致性校正模型,分別進(jìn)行俯仰、橫傾、方位3種姿態(tài)變換。利用最小二乘法進(jìn)行求解,收斂速度快,計算精度高。實驗結(jié)果表明,經(jīng)過磁強(qiáng)計系統(tǒng)誤差校正后,磁場測量均方根誤差在10nT以內(nèi),有效降低了磁通門磁強(qiáng)計的測量誤差。經(jīng)過一致性校正后,各軸測量角度誤差在0.01°以下,換算到磁場中為4.22nT以下,說明經(jīng)過校正后的磁強(qiáng)計陣列具有良好的一致性,有效解決了磁強(qiáng)計陣列姿態(tài)不一致導(dǎo)致的測量誤差。

    本文建立的磁強(qiáng)計陣列的校正模型為大平面測磁陣列的一致性校準(zhǔn)提供了思路和方法,有利于提高磁強(qiáng)計陣列的測量精度。該模型對測磁陣列系統(tǒng)無特殊要求,誤差校正模型基于單磁通門磁強(qiáng)計而建立,因此適用于任何由三軸磁通門磁強(qiáng)計或加速度計構(gòu)成的陣列,具有很高的實用性。

    猜你喜歡
    磁強(qiáng)計系統(tǒng)誤差磁通
    軸向磁通電勵磁雙凸極電機(jī)及容錯運(yùn)行控制策略
    三軸磁強(qiáng)計及梯度計校準(zhǔn)的綜合系數(shù)法
    基于矢量磁強(qiáng)計的磁場梯度張量儀誤差校正方法
    組合導(dǎo)航中磁強(qiáng)計干擾估計與補(bǔ)償方法
    基于LabVIEW的微型磁通門磁強(qiáng)計測試系統(tǒng)搭建
    基于磁通門原理的零磁通交直流電流傳感器
    基于ADS-B的航空器測高系統(tǒng)誤差評估方法
    基于Bagging模型的慣導(dǎo)系統(tǒng)誤差抑制方法
    存在系統(tǒng)誤差下交叉定位系統(tǒng)最優(yōu)交會角研究
    基于FPGA的數(shù)字磁通計設(shè)計
    電測與儀表(2015年3期)2015-04-09 11:37:52
    内地一区二区视频在线| 日韩一区二区三区影片| 亚洲精品国产av蜜桃| 91精品一卡2卡3卡4卡| 亚洲欧美一区二区三区黑人 | 中文乱码字字幕精品一区二区三区 | 久久精品国产鲁丝片午夜精品| 男女那种视频在线观看| 三级经典国产精品| eeuss影院久久| 国产精品日韩av在线免费观看| 两个人视频免费观看高清| 高清欧美精品videossex| 女人被狂操c到高潮| 中文字幕av在线有码专区| 大话2 男鬼变身卡| 免费不卡的大黄色大毛片视频在线观看 | 99热这里只有精品一区| 一二三四中文在线观看免费高清| 女人十人毛片免费观看3o分钟| 久久久久精品久久久久真实原创| 亚洲精品日韩av片在线观看| 亚洲国产高清在线一区二区三| 99久久九九国产精品国产免费| 欧美不卡视频在线免费观看| 22中文网久久字幕| 久久这里只有精品中国| 亚洲av成人av| 国产成人免费观看mmmm| 男女边吃奶边做爰视频| 91精品伊人久久大香线蕉| 亚洲婷婷狠狠爱综合网| 嫩草影院新地址| 超碰av人人做人人爽久久| 欧美性感艳星| 色哟哟·www| 国内揄拍国产精品人妻在线| 免费看a级黄色片| 亚洲欧洲国产日韩| 亚洲最大成人av| 亚洲婷婷狠狠爱综合网| 一二三四中文在线观看免费高清| 91精品伊人久久大香线蕉| 亚洲精品日本国产第一区| 晚上一个人看的免费电影| 观看免费一级毛片| 欧美 日韩 精品 国产| 亚洲成人久久爱视频| 最近2019中文字幕mv第一页| 国产乱来视频区| 久久精品国产亚洲av天美| 免费人成在线观看视频色| 搞女人的毛片| 亚洲国产高清在线一区二区三| 在线观看av片永久免费下载| 国产 一区精品| 亚洲av免费在线观看| 超碰av人人做人人爽久久| 国产精品嫩草影院av在线观看| 久久人人爽人人片av| 蜜桃亚洲精品一区二区三区| 99久国产av精品| 高清毛片免费看| 小蜜桃在线观看免费完整版高清| 欧美一级a爱片免费观看看| 九草在线视频观看| 3wmmmm亚洲av在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品一二三区在线看| 在线播放无遮挡| 国产精品国产三级专区第一集| 亚洲国产欧美人成| 最近2019中文字幕mv第一页| 欧美日本视频| 街头女战士在线观看网站| 综合色av麻豆| 韩国高清视频一区二区三区| 搞女人的毛片| 最新中文字幕久久久久| 亚洲aⅴ乱码一区二区在线播放| 美女高潮的动态| 高清视频免费观看一区二区 | 中文乱码字字幕精品一区二区三区 | 男女边摸边吃奶| 久久久久久久午夜电影| 建设人人有责人人尽责人人享有的 | 舔av片在线| 国产精品国产三级国产专区5o| 久久久精品免费免费高清| 国产乱人视频| 久久久精品欧美日韩精品| 青春草国产在线视频| 亚洲无线观看免费| 日本午夜av视频| 国产免费视频播放在线视频 | 国产亚洲91精品色在线| 能在线免费观看的黄片| 日韩在线高清观看一区二区三区| 亚洲精品自拍成人| 99九九线精品视频在线观看视频| 水蜜桃什么品种好| 99热这里只有是精品50| 边亲边吃奶的免费视频| 街头女战士在线观看网站| 一夜夜www| 久久久久久久国产电影| 国产精品1区2区在线观看.| 99热网站在线观看| 亚洲国产色片| 亚洲欧美成人精品一区二区| 精品人妻视频免费看| 两个人视频免费观看高清| av线在线观看网站| 欧美三级亚洲精品| 一级a做视频免费观看| 非洲黑人性xxxx精品又粗又长| 午夜视频国产福利| 久久99蜜桃精品久久| 免费观看性生交大片5| 亚洲美女搞黄在线观看| 亚洲第一区二区三区不卡| 高清av免费在线| 久久久久性生活片| 99re6热这里在线精品视频| 看非洲黑人一级黄片| 午夜福利在线观看吧| 国产成人精品一,二区| 男女边摸边吃奶| 日本免费a在线| 精品99又大又爽又粗少妇毛片| 国产精品一区二区三区四区免费观看| 在线观看一区二区三区| 日韩欧美精品免费久久| 亚洲av成人精品一区久久| 一级片'在线观看视频| 日韩av不卡免费在线播放| 国产精品精品国产色婷婷| 欧美成人精品欧美一级黄| 老司机影院毛片| 亚洲av福利一区| 免费黄色在线免费观看| 伊人久久精品亚洲午夜| 日韩制服骚丝袜av| 99久久精品国产国产毛片| 亚洲丝袜综合中文字幕| 国产一区二区三区av在线| 少妇猛男粗大的猛烈进出视频 | 免费大片黄手机在线观看| 日本wwww免费看| 国产老妇伦熟女老妇高清| 欧美日本视频| eeuss影院久久| 国产在视频线精品| 国产午夜精品一二区理论片| 日韩中字成人| 丰满少妇做爰视频| 日韩欧美一区视频在线观看 | 精品久久久久久成人av| 国内精品美女久久久久久| 丰满乱子伦码专区| 非洲黑人性xxxx精品又粗又长| 国产在线男女| 欧美日韩精品成人综合77777| 国产综合精华液| 欧美高清性xxxxhd video| 高清午夜精品一区二区三区| 一边亲一边摸免费视频| 国产单亲对白刺激| 全区人妻精品视频| 亚洲一级一片aⅴ在线观看| 亚洲欧美一区二区三区国产| 日韩国内少妇激情av| 免费观看a级毛片全部| 日韩欧美精品免费久久| 赤兔流量卡办理| 亚洲精品,欧美精品| 69av精品久久久久久| 真实男女啪啪啪动态图| 国产精品国产三级专区第一集| 综合色丁香网| 尤物成人国产欧美一区二区三区| 亚洲色图av天堂| 国产 亚洲一区二区三区 | 美女内射精品一级片tv| 国产v大片淫在线免费观看| 欧美区成人在线视频| 丰满乱子伦码专区| 亚洲精品一区蜜桃| 婷婷色麻豆天堂久久| 天天躁日日操中文字幕| 亚洲丝袜综合中文字幕| 日本午夜av视频| 成人无遮挡网站| videos熟女内射| 99re6热这里在线精品视频| 全区人妻精品视频| 国产精品国产三级专区第一集| 欧美bdsm另类| 亚洲精品久久午夜乱码| 超碰97精品在线观看| 亚洲欧美日韩东京热| 日产精品乱码卡一卡2卡三| 欧美性猛交╳xxx乱大交人| 乱系列少妇在线播放| 婷婷色综合www| 免费看不卡的av| 日韩欧美精品v在线| 91在线精品国自产拍蜜月| 久久久久精品性色| 亚洲在线自拍视频| 亚洲三级黄色毛片| 看十八女毛片水多多多| 日韩伦理黄色片| 欧美xxxx黑人xx丫x性爽| 成人午夜精彩视频在线观看| 日韩中字成人| 国产黄色视频一区二区在线观看| 欧美激情久久久久久爽电影| 亚洲内射少妇av| 欧美日韩亚洲高清精品| 国产高潮美女av| 一区二区三区四区激情视频| 午夜免费观看性视频| 精品国产三级普通话版| 麻豆乱淫一区二区| 欧美+日韩+精品| 婷婷色综合大香蕉| 噜噜噜噜噜久久久久久91| 精品一区二区三区视频在线| 一级毛片aaaaaa免费看小| 亚洲国产日韩欧美精品在线观看| 神马国产精品三级电影在线观看| 国产精品国产三级国产av玫瑰| 午夜福利视频精品| 亚洲aⅴ乱码一区二区在线播放| a级毛色黄片| 欧美+日韩+精品| 亚洲欧美日韩无卡精品| 日韩伦理黄色片| 能在线免费看毛片的网站| 尾随美女入室| 国产av不卡久久| 日韩一区二区三区影片| 欧美日韩一区二区视频在线观看视频在线 | 女人久久www免费人成看片| 成人特级av手机在线观看| 国产在视频线在精品| 亚洲乱码一区二区免费版| 春色校园在线视频观看| 免费大片黄手机在线观看| 丰满乱子伦码专区| 久久久久久久亚洲中文字幕| 黄片wwwwww| 亚洲精品第二区| 欧美区成人在线视频| 视频中文字幕在线观看| 国产在线男女| 亚洲av福利一区| 高清毛片免费看| 自拍偷自拍亚洲精品老妇| 在现免费观看毛片| 欧美激情久久久久久爽电影| 一级二级三级毛片免费看| 日韩视频在线欧美| 日日撸夜夜添| 内地一区二区视频在线| 深夜a级毛片| 国产美女午夜福利| 永久网站在线| 午夜老司机福利剧场| 高清视频免费观看一区二区 | 欧美变态另类bdsm刘玥| 国产 亚洲一区二区三区 | 亚洲精品第二区| 亚洲第一区二区三区不卡| 欧美xxxx黑人xx丫x性爽| 亚洲av中文字字幕乱码综合| 国产午夜精品久久久久久一区二区三区| 色吧在线观看| 男女视频在线观看网站免费| 国产精品麻豆人妻色哟哟久久 | 亚洲18禁久久av| 五月天丁香电影| 欧美日本视频| 亚洲精品视频女| 精品久久久久久久久av| 少妇人妻精品综合一区二区| 国产伦精品一区二区三区四那| 成人亚洲精品一区在线观看 | 精品久久久久久久久av| 观看美女的网站| 亚洲精品国产av蜜桃| 美女脱内裤让男人舔精品视频| 国产成人91sexporn| 91精品伊人久久大香线蕉| 丝瓜视频免费看黄片| 欧美区成人在线视频| 久久精品综合一区二区三区| 成年免费大片在线观看| 成人美女网站在线观看视频| 精品人妻偷拍中文字幕| 欧美性感艳星| 亚洲图色成人| 欧美精品一区二区大全| 少妇裸体淫交视频免费看高清| 中文字幕亚洲精品专区| 欧美日本视频| 亚洲综合色惰| 亚洲精品日本国产第一区| 亚洲在久久综合| 最近最新中文字幕免费大全7| 午夜老司机福利剧场| 免费看不卡的av| 久久久国产一区二区| 国产 一区 欧美 日韩| 看非洲黑人一级黄片| 啦啦啦啦在线视频资源| 婷婷色麻豆天堂久久| 2021少妇久久久久久久久久久| 人妻一区二区av| 麻豆国产97在线/欧美| 99热这里只有是精品50| 美女黄网站色视频| 国产亚洲av片在线观看秒播厂 | 日本一二三区视频观看| 成人亚洲欧美一区二区av| 综合色丁香网| 国产欧美日韩精品一区二区| 久久99热这里只有精品18| 秋霞在线观看毛片| 久久精品人妻少妇| 亚洲成人久久爱视频| 我的老师免费观看完整版| 婷婷色综合www| 色网站视频免费| 欧美zozozo另类| 十八禁网站网址无遮挡 | 精品酒店卫生间| 一级二级三级毛片免费看| 亚洲人成网站在线观看播放| 一级a做视频免费观看| 亚洲三级黄色毛片| 777米奇影视久久| 国产高清三级在线| 丰满人妻一区二区三区视频av| 啦啦啦啦在线视频资源| 蜜臀久久99精品久久宅男| 国产黄色视频一区二区在线观看| 欧美日韩一区二区视频在线观看视频在线 | 色哟哟·www| 联通29元200g的流量卡| 日本猛色少妇xxxxx猛交久久| 高清在线视频一区二区三区| 色综合亚洲欧美另类图片| 亚洲精品乱码久久久v下载方式| 国产 一区精品| 亚洲一区高清亚洲精品| 97人妻精品一区二区三区麻豆| 国产麻豆成人av免费视频| 亚洲精品aⅴ在线观看| 成年人午夜在线观看视频 | 91在线精品国自产拍蜜月| 久久精品国产亚洲av涩爱| 校园人妻丝袜中文字幕| 免费看光身美女| 免费看美女性在线毛片视频| 91精品一卡2卡3卡4卡| 亚洲精品国产成人久久av| 午夜视频国产福利| 听说在线观看完整版免费高清| 精品久久久久久久人妻蜜臀av| 国产精品一区二区三区四区久久| 男女那种视频在线观看| 777米奇影视久久| 欧美97在线视频| 高清av免费在线| 美女被艹到高潮喷水动态| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产成人一精品久久久| 久久国产乱子免费精品| 99热这里只有精品一区| 久久久久久久久大av| 国产综合懂色| 五月玫瑰六月丁香| 国产探花在线观看一区二区| 男女边吃奶边做爰视频| 少妇的逼水好多| 色哟哟·www| 赤兔流量卡办理| 国产高清国产精品国产三级 | 精品久久久久久成人av| 色综合亚洲欧美另类图片| 国产在视频线在精品| 97超视频在线观看视频| 联通29元200g的流量卡| or卡值多少钱| 亚洲在久久综合| 欧美zozozo另类| av福利片在线观看| 一本久久精品| 男的添女的下面高潮视频| 三级毛片av免费| 91午夜精品亚洲一区二区三区| 久久久久精品性色| 国产亚洲5aaaaa淫片| 麻豆乱淫一区二区| 日韩中字成人| 国产在线一区二区三区精| 亚洲av男天堂| 久久精品熟女亚洲av麻豆精品 | 午夜精品国产一区二区电影 | 草草在线视频免费看| 内地一区二区视频在线| 爱豆传媒免费全集在线观看| 午夜免费观看性视频| 搡老乐熟女国产| 伦精品一区二区三区| 91在线精品国自产拍蜜月| 成人亚洲精品一区在线观看 | 日韩伦理黄色片| 久久精品国产亚洲av天美| 一级毛片 在线播放| 国产成人a区在线观看| 国产亚洲91精品色在线| 国产亚洲精品av在线| 深爱激情五月婷婷| 网址你懂的国产日韩在线| 能在线免费看毛片的网站| 日本黄大片高清| 欧美xxxx性猛交bbbb| 亚洲精品,欧美精品| 女的被弄到高潮叫床怎么办| 99久久中文字幕三级久久日本| 99久国产av精品| 精品久久久久久电影网| 国产探花在线观看一区二区| 2018国产大陆天天弄谢| 国产成人一区二区在线| 性插视频无遮挡在线免费观看| 国产精品1区2区在线观看.| 最近手机中文字幕大全| 91午夜精品亚洲一区二区三区| 亚洲精品第二区| 日本色播在线视频| 亚洲精品一区蜜桃| 少妇的逼好多水| 人人妻人人看人人澡| 美女xxoo啪啪120秒动态图| 精品少妇黑人巨大在线播放| 性插视频无遮挡在线免费观看| 久久久久久久久久人人人人人人| 亚洲va在线va天堂va国产| 一区二区三区高清视频在线| 亚洲一区高清亚洲精品| freevideosex欧美| 国产单亲对白刺激| 色吧在线观看| 18+在线观看网站| 91在线精品国自产拍蜜月| 色网站视频免费| 欧美成人精品欧美一级黄| 久久精品久久久久久噜噜老黄| 超碰97精品在线观看| 国产真实伦视频高清在线观看| 免费黄色在线免费观看| 亚洲国产精品sss在线观看| 熟妇人妻不卡中文字幕| 日韩大片免费观看网站| 亚洲av在线观看美女高潮| 美女国产视频在线观看| 成人午夜精彩视频在线观看| 女人十人毛片免费观看3o分钟| 搞女人的毛片| 一区二区三区免费毛片| 高清午夜精品一区二区三区| 久久久久网色| 91aial.com中文字幕在线观看| 国产v大片淫在线免费观看| 美女国产视频在线观看| 中文字幕av在线有码专区| 听说在线观看完整版免费高清| 亚洲av成人精品一区久久| av一本久久久久| 中文资源天堂在线| 国产精品久久久久久av不卡| 欧美精品一区二区大全| 国产亚洲一区二区精品| 亚洲国产精品专区欧美| 日韩欧美国产在线观看| 2018国产大陆天天弄谢| 少妇熟女aⅴ在线视频| 禁无遮挡网站| 中文欧美无线码| 国产在视频线在精品| 久久精品综合一区二区三区| 伦理电影大哥的女人| 五月天丁香电影| 中文欧美无线码| 亚洲天堂国产精品一区在线| 国产视频首页在线观看| 国产黄片美女视频| 国产人妻一区二区三区在| 亚洲成人中文字幕在线播放| 精品久久久噜噜| 成人漫画全彩无遮挡| h日本视频在线播放| 国产淫片久久久久久久久| 日韩欧美精品v在线| 乱码一卡2卡4卡精品| 久久久精品欧美日韩精品| 波野结衣二区三区在线| 少妇人妻精品综合一区二区| 亚洲怡红院男人天堂| 国产伦精品一区二区三区四那| 男女边吃奶边做爰视频| 女人久久www免费人成看片| 在线 av 中文字幕| 国产精品国产三级专区第一集| 欧美成人精品欧美一级黄| 国产一区亚洲一区在线观看| 高清视频免费观看一区二区 | www.av在线官网国产| 日本欧美国产在线视频| 哪个播放器可以免费观看大片| 久久精品国产亚洲av天美| 久久精品综合一区二区三区| 网址你懂的国产日韩在线| 国产人妻一区二区三区在| 看非洲黑人一级黄片| 日韩电影二区| 国产免费视频播放在线视频 | 又爽又黄无遮挡网站| 最近最新中文字幕免费大全7| 亚洲av中文字字幕乱码综合| 国产精品久久久久久av不卡| 日韩精品青青久久久久久| 亚洲国产精品国产精品| or卡值多少钱| 国产有黄有色有爽视频| 女人被狂操c到高潮| 五月玫瑰六月丁香| 看免费成人av毛片| 国产人妻一区二区三区在| 国产一区二区在线观看日韩| 国产欧美日韩精品一区二区| 国产成人a∨麻豆精品| 男女边摸边吃奶| 精品亚洲乱码少妇综合久久| or卡值多少钱| 你懂的网址亚洲精品在线观看| 秋霞在线观看毛片| 国产 亚洲一区二区三区 | 秋霞伦理黄片| 国产精品久久久久久av不卡| 国产黄色免费在线视频| 国产大屁股一区二区在线视频| 2018国产大陆天天弄谢| 久久久久久久久久久免费av| 一级黄片播放器| 国产 亚洲一区二区三区 | 青春草视频在线免费观看| 免费观看无遮挡的男女| 国产精品一区二区性色av| 欧美极品一区二区三区四区| 久久久久久国产a免费观看| 网址你懂的国产日韩在线| 久久人人爽人人片av| 亚洲一级一片aⅴ在线观看| 亚洲在久久综合| 在线免费十八禁| 日韩精品有码人妻一区| 97超碰精品成人国产| 亚洲综合精品二区| 国产真实伦视频高清在线观看| 十八禁国产超污无遮挡网站| 91在线精品国自产拍蜜月| 在线观看av片永久免费下载| 精品久久久久久久久久久久久| 夫妻午夜视频| 夫妻性生交免费视频一级片| 国产69精品久久久久777片| 国产一区亚洲一区在线观看| 国产成人精品福利久久| 久久精品夜色国产| 国产精品国产三级国产av玫瑰| 久久久欧美国产精品| 最近手机中文字幕大全| 美女主播在线视频| 中文字幕av在线有码专区| 亚洲综合精品二区| 五月伊人婷婷丁香| 少妇的逼水好多| 春色校园在线视频观看| 黄片wwwwww| 国产成人午夜福利电影在线观看| 亚洲国产最新在线播放| 精品人妻视频免费看| 午夜免费观看性视频| 亚洲va在线va天堂va国产| 五月伊人婷婷丁香| 赤兔流量卡办理| 免费观看的影片在线观看| 久久精品久久久久久噜噜老黄| 联通29元200g的流量卡| av国产久精品久网站免费入址| 国产极品天堂在线| 午夜精品一区二区三区免费看| 嫩草影院新地址| 亚洲精品视频女| 国产成人a区在线观看| 毛片一级片免费看久久久久| 亚洲欧美精品专区久久| 久久这里有精品视频免费| 久久午夜福利片| 亚洲精品视频女|