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

    電磁層析成像系統(tǒng)中標量磁勢的數(shù)值解法

    2014-06-05 15:30:34郝建娜尹武良
    關鍵詞:磁勢層析成像標量

    趙 倩,郝建娜,尹武良

    (天津大學電氣與自動化工程學院,天津 300072)

    電磁層析成像系統(tǒng)中標量磁勢的數(shù)值解法

    趙 倩,郝建娜,尹武良

    (天津大學電氣與自動化工程學院,天津 300072)

    電磁層析成像技術是一種基于電磁感應定律的工業(yè)過程成像技術,激勵線圈產(chǎn)生的交變磁場在目標物體中產(chǎn)生渦流,進而產(chǎn)生二次磁場.接收線圈檢測到感應電壓后利用重建算法可以得到物場的分布信息.邊界元法以積分方程為數(shù)學基礎,同時采用了與有限元法相似的劃分單元離散技術,將邊界積分方程離散為代數(shù)方程組后用數(shù)值方法求解.邊界元法在電磁層析成像技術中已有一定的應用,而基于標量磁勢的邊界元法使得求解過程的速度和效率顯著提高.本文針對一個簡化的電磁層析成像系統(tǒng)模型,利用3種數(shù)值方法求解目標物體表面的標量磁勢,并利用Matlab編程得到仿真結果.對結果進行對比分析后,選出最優(yōu)解法.

    電磁層析成像;邊界元法;標量磁勢;數(shù)值仿真

    電磁層析成像(electromagnetic tomography,EMT)技術是20世紀90年代發(fā)展起來的一種新型電學層析成像(electrical tomography,ET)技術,EMT技術基于電磁感應原理,載流線圈在目標物體中產(chǎn)生渦流,進而激勵出二次磁場,這樣接收線圈可以測量到感應電壓.最后采用一定的圖像重建算法得到被測場域內的分布信息.它具有非接觸、非侵入、響應快速、信息量豐富等特點,因而在工業(yè)生產(chǎn)過程和生物醫(yī)學領域都有著廣泛的應用前景[1-7].

    邊界元法(boundary element method,BEM)具有邊界積分方程的深厚的數(shù)學根基.它以邊界積分方程為數(shù)學基礎,通過單元離散技術將邊界離散為邊界元,將邊界積分方程離散為代數(shù)方程組.再用數(shù)值方法求解代數(shù)方程組,從而得到原問題邊界積分方程的解.邊界元法的優(yōu)點為可以降低求解問題的維數(shù),提高了求解效率,并且由于求解過程中采用了解析基本解,所以具有更高的精度[8-14].

    BEM在EMT中的應用已經(jīng)有了許多的研究[15].在邊界元法求解電磁問題的過程中,標量磁勢作為基本變量具有重要意義.在求解外部磁場分布之前,需要先求出目標物體表面的標量磁勢,以此為基礎再進行后續(xù)計算.目前對于標量磁勢的求解方法一般表示為對要求點處求立體角[16],該方法同樣涉及如何求解積分的問題.對于復雜不規(guī)則物體,快速多極算法是一種快速有效的計算方法[17].但是該算法的步驟較為繁瑣,對于一些簡單問題反倒不適用.

    筆者通過建立一個簡單的EMT系統(tǒng)模型,推導出標量磁勢的最終表達式.通過3種不同的數(shù)值方法——曲面積分定義法、高斯-勒讓德積分公式法和積分坐標變換法,分別仿真得到不同的結果.對這些數(shù)據(jù)進行分析可以得出,在保證其他條件不變和相同求解精度的前提下,高斯-勒讓德積分公式法是最簡單快速的方法.

    1 建立模型并推導公式

    假設一個簡單的EMT系統(tǒng),只含有1個激勵線圈和1個接收線圈.激勵線圈和接收線圈的中心分別位于(0,0,4,m)和(0,0,-4,m)處,半徑均為1,m.目標球體為一般金屬,中心位于原點,半徑為1,m.

    對于恒定磁場,由磁場環(huán)路定理可得

    式中L為曲面S的邊界.這種情況下如果想引入標量磁勢,其前提條件是,對于求解區(qū)域內的任何閉合回路,都有

    根據(jù)上面推導出的公式,這里采用了3種方法來求解球體表面離散點處的標量磁勢.之后利用Matlab仿真,對這3種方法的求解速度和所得數(shù)據(jù)精度進行了分析對比.

    圖1 一個簡單的EMT系統(tǒng)Fig.1 A simple EMT system

    在實際應用中,對于單連通,不存在自由電流的區(qū)域,全區(qū)域皆可以引入標量磁勢描述磁場.如果區(qū)域中存在自由電流,則挖去電流及電流線圍著的一個曲面S,剩余空間里磁場是無旋的[18-19],可以假設為標量磁勢,A.

    圖1所示的EMT系統(tǒng)中,載流線圈在物體表面產(chǎn)生的磁場可由畢奧-薩伐兒定理求出[19],即

    式中:上標pr表示外部激勵產(chǎn)生的場量;B為磁感應強度;r=r1-r2,r1為球面上的點到原點的矢徑,r2為激勵線圈所圍平面上的點到原點的矢徑.

    式中:S為由激勵線圈L圍成的平面;n0為S的單位法向量.

    對場域中的球體表面進行三角形網(wǎng)格劃分,如圖2所示,可以得到38個離散點和72個網(wǎng)格.劃分時,縱軸z的范圍為1,m到-1,m,每隔0.2,m取一組,對應11個不同的z值.在每個z確定的圓上,角度范圍為0到2π,間隔為π/2.

    圖2 EMT系統(tǒng)中球體表面的網(wǎng)格劃分Fig.2 Mesh of the target surface in the EMT system

    若激勵線圈上載有單位電流,即I=| I |=1,A,則球體表面由外加磁場產(chǎn)生的標量磁勢為

    2 Matlab數(shù)值仿真

    2.1 由曲面積分的定義求解(方法1)

    由曲面積分的定義可得[20],通過將積分區(qū)域劃分成正方形網(wǎng)格,可將積分離散化,有

    式中:iSΔ為第i個正方形網(wǎng)格的面積;2ir為激勵線圈所圍平面上的點到原點的矢徑.在式(6)中,積分平面為激勵線圈所包圍的圓面,將該平面分成N個小正方形,r2i為矢徑取第i個正方形的中心到原點的矢徑.首先對r2i的位置進行判斷,若是在圓內則將其帶入積分公式;若在圓外,則其對積分的作用為0;再將所有的積分進行求和,最終得到在球體表面每個離散點處的標量磁勢.

    由仿真結果可得,隨著劃分網(wǎng)格數(shù)目的增多,積分結果成收斂趨勢.當網(wǎng)格數(shù)為6,400時,積分基本不再變,計算時間為t=0.556,122,s.38個離散點對應11個z值,通過仿真可知,每個z值相同的點具有相等的標量磁勢,故可用11個不同z值上的標量磁勢來表示所有離散點處的標量磁勢.

    表1 方法1得到的標量磁勢Tab.1 Magnetic scalar potential obtained by the first method

    以表1第1行為例,可以畫出隨著網(wǎng)格數(shù)增加標量磁勢逐漸收斂的仿真圖形,如圖3所示.

    由于每個確定的z值對應的圓上的4個點處具有相同的標量磁勢,所以下面的表格中可以只列出不同z值處的標量磁勢.圖4為x=0處的球體切面上的等勢線,可以看出,球體從上往下標量磁勢逐漸減小,并且在每個確定z值處具有相同的標量磁勢.

    圖3 隨著網(wǎng)格數(shù)標量磁勢數(shù)值逐漸收斂Fig.3Convergence of the magnetic scalar potential as the increase of mesh

    圖4 x=0處的切面上等勢線分布Fig.4 Isopotential lines on the section where x=0

    2.2 利用N段高斯-勒讓德積分公式來求解(方法2)

    首先將面積分轉換成二重積分[21]為

    則原積分可以簡寫為

    積分變量θ和r′分別表示角度和半徑,積分函數(shù)為

    對積分區(qū)間進行標準化得

    變換方程為

    式中:n為節(jié)點數(shù);k為插值節(jié)點數(shù);Ai和Aj(i,j=1,2,…,R)均為系數(shù)矩陣A中第i個和第j個值.具體的節(jié)點和對應的系數(shù)矩陣如表2所示.為了使精度更高,在這里先將積分區(qū)域劃分為N段小區(qū)域,然后在每段小區(qū)域上進行插值積分運算.運行程序得到的結果表明,當N=2,即分成2段,然后進行5點插值可得高精度的解,用時為t=0.109,799,s.

    表2 高斯-勒讓德積分公式的節(jié)點及系數(shù)矩陣Tab.2 Nodes and coefficient matrices of Gaussian-Legendre integral formula

    由表3可知,對于N(N=1,2,4)段積分公式而言,所求得的標量磁勢基本與N無關.只有在z=0時,不同N對應的求解結果才有細微差別.故可知,利用N段高斯-勒讓德積分公式來求解標量磁勢時,段數(shù)N對結果精度的影響甚微,從而為了節(jié)約計算時間,可以直接使用1段5點積分公式來求解.

    表3 方法2得到的標量磁勢Tab.3 Magnetic scalar potential obtained by the second method

    2.3 利用積分坐標變換法來求解(方法3)

    這種方法首先將積分平面劃分為三角形網(wǎng)格,然后利用坐標系轉換得到的積分公式來求解[22-24].

    劃分網(wǎng)格時,先將圓面上的點進行離散,然后按照逆時針順序將點依次連接為三角形.根據(jù)其中心點的位置,可以判斷出三角形是否在圓內.為了方便求解,在這里將全局坐標系轉換為局部坐標系.對每個網(wǎng)格而言,新坐標系平行于該網(wǎng)格所在的平面,并且以其中的一邊為其橫坐標.坐標系的各參數(shù)項如圖5所示.

    圖5 局部坐標系與全局坐標系的關系及局部坐標系中的參數(shù)Fig.5 Relationship between the local coordinate system and the global coordinate system and the parameters of the local coordinate system

    采用與文獻[22]中類似的參數(shù)和函數(shù),這里以一個三角形網(wǎng)格為例.

    V1、V2和V3為第e個三角形網(wǎng)格eΓ的3個頂點.iS?為eΓ的第i條邊.li、si和mi(i=1,2,3)分別是iS?的長度、單位切向量和單位法向量.則局部坐標系可以定義為

    根據(jù)以上所列參數(shù),則原積分可化為

    其中n0=(0,0,1),ξ′=0.故有

    式中Ne為三角形網(wǎng)格的個數(shù).

    由文獻[22]可得eΓ上的積分為

    數(shù)值計算表明,當網(wǎng)格數(shù)為12,800時,標量磁勢可以收斂,但用時較長,t=157.126,357,s.

    表4 方法3得到的標量磁勢Tab.4 Magnetic scalar potential obtained by the third method

    2.4 利用Matlab進行數(shù)據(jù)對比分析

    利用Matlab軟件對3種方法得到的數(shù)據(jù)進行分析,如圖6所示.

    從圖6可以看出,3種方法得到的數(shù)值解具有很好的一致性,即3種方法都可以得到精確的解,從程序運行時間上看,高斯-勒讓德積分公式用時最短,t=0.109,799,s,積分坐標變換法用時最長,t= 157.126,357,s.

    圖6 3種方法得到的數(shù)據(jù)曲線Fig.6 Results obtained by the three methods

    3 結 語

    快速求解標量磁勢為后續(xù)計算磁場分布提供了良好基礎.本文采用了3種數(shù)值方法來求解物體表面的標量磁勢,通過以上計算可以看出,以上3種方法都能有效求解目標物體表面的標量磁勢.若綜合考慮的話,高斯-勒讓德積分公式是最好的計算方法,不僅精度高,而且速度快.

    [1] Griffiths H. Magnetic induction tomography[J]. Meas Sci Technol,2001,12:1126-1131.

    [2] Watson S,Wee H C,Griffiths H,et al. A highly phase-stable differential detector amplifier for magnetic induction tomography[J]. Physiol Meas,2011,32:917-926.

    [3] Watson S,Williams R J,Gough W A,et al. A magnetic induction tomography system for samples with conductivities less than 10 S·m-1[J]. Meas Sci Technol,2008,19:1-11.

    [4] Yin W L,Peyton A J. A planar EMT system for the detection of faults on thin metallic plates[J]. Meas Sci Technol,2006,17:2130-2135.

    [5] Soleimani M,Lionheart W R B,Riedel C H,et al. Forward problem in 3,D magnetic induction tomography(MIT)[C]//Proc 3rd World Congr Industrial Process Tomography. Banff,Canada,2003:275-280.

    [6] Yin W L,Peyton A J. Sensitivity formulation including velocity effects for electromagnetic induction systems[J]. IEEE Transactions on Magnetics,2010,46(5):1172-1176.

    [7] Yin W L,Peyton A J,Zysko G,et al. Simultaneous noncontact measurement of water-level and conductivity [J]. IEEE Transactions on Instrumentation and Measurement,2008,57:2665-2669.

    [8] Brebbia C A,Telles J C F,Wrobel L C. Boundary Element Techniques:Theory and Applications in Engineering[M]. New York:Springer -Verlag,1984.

    [9] 姚振漢,王海濤. 邊界元法[M]. 北京:高等教育出版社,2010.

    Yao Zhenhan,Wang Haitao. The Boundary Element Method[M]. Beijing:Higher Education Press,2010(in Chinese).

    [10] 劉永健,姚振漢. 三維接觸邊界元法的一種誤差直接估計[J]. 清華大學學報:自然科學版,2003,43(11):1499-1502.

    Liu Yongjian,Yao Zhenhan. Direct error estimate for the boundary element method for a 3-D contact problem[J]. Journal of Tsinghua Univ:Sci & Tech,2003,43(11):1499-1502(in Chinese).

    [11] 單 磊,董天臨. 應用邊界元法對電磁場計算的分析及優(yōu)化[J]. 信息通信,2007,33(1):31-33.

    Shan Lei,Dong Tianlin. Analysis and optimization of boundary element method to irregular electromagnetic field[J]. Information & Communications,2007,33(1):31-33(in Chinese).

    [12] Harrington R F. Field Computation by Moment Methods[M]. New York:Acmillan Company,1968.

    [13] Banerjee P K. The Boundary Element Methods in Engineering[M]. New York:McGraw-Hill Book Company,1994.

    [14] Chen W,F(xiàn)u Z. A novel numerical method for infinite domain potential problems[J]. Chinese Science Bulletin,2010,55(16):1598-1603.

    [15] Pham M H,Peyton A J. A model for the forward problem in magnetic induction tomography using boundary integral equations[J]. IEEE Transactions on Magnetics,2008,44(10):2262 -2267.

    [16] 張澤明,鄧亦仁. 運用畢-沙定律導出磁標量位表達式的新方法[J]. 重慶建筑工程學院學報,1985(4):99-103.

    Zhang Zeming,Deng Yiren. A new method using Biot-Savart Law to derive magnetic scalar potential notation[J]. Journal of Chongqing Institute of Civil Engineering and Architecture,1985(4):99-103(in Chinese).

    [17] Hafla W,Buchau A,Rucker W M. Efficient computation of source magnetic scalar potential [J]. Adv Radio Sci,2006,4:59-63.

    [18] 阮放鳴. 有限區(qū)域里的磁標勢和磁多極展開[J]. 貴州師范大學學報:自然科學版,1997,15(3):102-107.

    Ruan Fangming. Magnetic scalar potentials in finite regious and the multipole expansion in material media[J]. Journal of Guizhou Normal University:Natural Science,1997,15(3):102-107(in Chinese).

    [19] 謝處方,饒克謹. 電磁場與電磁波[M]. 北京:高等教育出版社,1987.

    Xie Chufang,Rao Kejin. The Electromagnetic Field and Waves[M]. Beijing:Higher Education Press,1987(in Chinese).

    [20] 同濟大學數(shù)學教研室. 高等數(shù)學 [M]. 4版. 北京:高等教育出版社,1996.

    The Mathematical Teaching Team of Tongji University. Advanced Mathematics[M]. 4th ed. Beijing:Higher Education Press,1996(in Chinese).

    [21] 張 軍. 數(shù)值計算[M]. 北京:清華大學出版社,2008.

    Zhang Jun. Numerical Computation[M]. Beijing:Tsinghua University Press,2008(in Chinese).

    [22] Graglia R D. On the numerical integration of the linear shape functions times the 3-D Green's function or its gradient on a plane triangle [J]. IEEE Transactions on Antennas and Propagation,1993,41(10):1448-1455.

    [23] Morrison J A. Integral equations for electromagnetic scattering by perfect conductors with two-dimensional geometry [J]. Bell Syst Tech J,1979,58(2):409-425.

    [24] Graglia R D. Static and dynamic potential integrals for linearly varying source distributions in two- and threedimensional problems [J]. IEEE Transactions on Antennas and Propagation,1987,AP-35(6):662-669.

    (責任編輯:孫立華)

    Numerical Approach for the Magnetic Scalar Potential of Electromagnetic Tomography System

    Zhao Qian,Hao Jianna,Yin Wuliang
    (School of Electrical Engineering and Automation,Tianjin University,Tianjin 300072,China)

    Electromagnetic tomography (EMT) is an industrial process tomography based on the Faraday law of electromagnetic induction. In EMT, current-carrying coils are used to induce eddy currents in the object which in turn radiates a scattered field. The induced voltages are sensed with receiving coils and image reconstruction algorithms are used to solve for the spatial distribution of the complex conductivity in the target. Boundary element method (BEM) based on integral formulations is an effective way to analyze the problems of EMT system. By adopting the similar mesh technique as FEM, the boundary integral equations can be transformed into linear equations. Then numerical method is used to solve the linear equations and the solution of the original integral equations can be obtained. The magnetic scalar potential was the important part of the application of BEM into EMT which improves efficiency and speed. Three numerical approaches were used to calculate the magnetic scalar potential on the surface of the target for a simple EMT model. With the aid of Matlab, the optimal method was chosen by comparing the results obtained by the three methods.

    electromagnetic tomography;boundary element method;magnetic scalar potential;numerical simulation

    O241.8

    A

    0493-2137(2014)07-0613-06

    10.11784/tdxbz201210003

    2012-10-06;

    2013-01-23.

    國家自然科學基金國際重大合作資助項目(60910001);教育部博士點基金資助項目(20090032110062).

    趙 倩(1987— ),女,博士研究生,shmilyshenzhen@163.com.

    尹武良,wuliang.yin@gmail.com.

    猜你喜歡
    磁勢層析成像標量
    變壓器和異步機激磁磁勢研究
    防爆電機(2022年2期)2022-04-26 13:57:02
    基于大數(shù)據(jù)量的初至層析成像算法優(yōu)化
    基于快速行進法地震層析成像研究
    同步機和異步機之磁勢和電勢的比較
    防爆電機(2020年4期)2020-12-14 03:11:02
    一種高效的橢圓曲線密碼標量乘算法及其實現(xiàn)
    一種靈活的橢圓曲線密碼并行化方法
    交、直流電機繞組磁動勢求取方法之比較
    防爆電機(2016年6期)2016-12-18 13:35:18
    基于分布式無線網(wǎng)絡的無線電層析成像方法與實驗研究
    雷達學報(2014年4期)2014-04-23 07:43:22
    單調Minkowski泛函與Henig真有效性的標量化
    基于多級小波域變換的時域擴散熒光層析成像方法
    天天添夜夜摸| 两个人看的免费小视频| 性高湖久久久久久久久免费观看| 免费黄色在线免费观看| 美女福利国产在线| 成人三级做爰电影| 欧美日韩福利视频一区二区| 欧美精品一区二区大全| 国产野战对白在线观看| 男人爽女人下面视频在线观看| 91精品国产国语对白视频| 亚洲美女视频黄频| 1024视频免费在线观看| 少妇精品久久久久久久| 9热在线视频观看99| 久久国产亚洲av麻豆专区| 99九九在线精品视频| 亚洲精品视频女| 久久av网站| 亚洲精品国产区一区二| 欧美日韩av久久| 国产99久久九九免费精品| 最近最新中文字幕免费大全7| 男女边摸边吃奶| 国产女主播在线喷水免费视频网站| 亚洲成色77777| h视频一区二区三区| 欧美国产精品va在线观看不卡| 嫩草影院入口| 搡老乐熟女国产| 久久久精品94久久精品| 波多野结衣一区麻豆| 国产一区二区激情短视频 | 午夜免费观看性视频| 亚洲av综合色区一区| 老司机影院成人| 精品少妇内射三级| 久久久久人妻精品一区果冻| 啦啦啦在线观看免费高清www| 日韩一卡2卡3卡4卡2021年| 国产精品国产三级专区第一集| 51午夜福利影视在线观看| www.av在线官网国产| 久久久欧美国产精品| 亚洲第一av免费看| 黄色 视频免费看| 国产黄色视频一区二区在线观看| 午夜日韩欧美国产| 制服人妻中文乱码| 亚洲天堂av无毛| 亚洲av日韩在线播放| 国产av一区二区精品久久| 热re99久久精品国产66热6| 中文字幕精品免费在线观看视频| 亚洲精品国产色婷婷电影| 美女大奶头黄色视频| 19禁男女啪啪无遮挡网站| 91精品三级在线观看| 日韩电影二区| 美女午夜性视频免费| 亚洲国产成人一精品久久久| 亚洲国产成人一精品久久久| 亚洲欧美一区二区三区久久| 国产99久久九九免费精品| 久久久久精品国产欧美久久久 | 男女高潮啪啪啪动态图| 日本午夜av视频| 母亲3免费完整高清在线观看| av有码第一页| 亚洲成av片中文字幕在线观看| 少妇被粗大的猛进出69影院| 在线亚洲精品国产二区图片欧美| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品自拍成人| 亚洲成色77777| 看免费成人av毛片| 欧美黄色片欧美黄色片| 97在线人人人人妻| 亚洲伊人色综图| 99久久人妻综合| 亚洲成人一二三区av| 亚洲欧美一区二区三区国产| www日本在线高清视频| 国产一区二区 视频在线| 高清av免费在线| 免费观看a级毛片全部| 日本av免费视频播放| av在线观看视频网站免费| 久久久久久久久久久久大奶| 最近中文字幕2019免费版| 久久人妻熟女aⅴ| 九草在线视频观看| 精品国产露脸久久av麻豆| 国产av精品麻豆| 制服丝袜香蕉在线| 最近最新中文字幕大全免费视频 | 国产探花极品一区二区| 国产精品二区激情视频| 亚洲男人天堂网一区| 中文字幕人妻丝袜一区二区 | 性高湖久久久久久久久免费观看| 久久久国产精品麻豆| 国产av码专区亚洲av| 99国产综合亚洲精品| 亚洲欧美一区二区三区黑人| 日韩一本色道免费dvd| 母亲3免费完整高清在线观看| 18禁国产床啪视频网站| 亚洲欧美色中文字幕在线| av不卡在线播放| 无限看片的www在线观看| 中文字幕精品免费在线观看视频| 亚洲欧美清纯卡通| 欧美日韩精品网址| 亚洲欧美日韩另类电影网站| 日韩 欧美 亚洲 中文字幕| 岛国毛片在线播放| 日韩精品免费视频一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 亚洲三区欧美一区| 婷婷色av中文字幕| 久久久久精品性色| 日韩成人av中文字幕在线观看| 一级毛片电影观看| 久久久久精品国产欧美久久久 | 美女国产高潮福利片在线看| 国产男女内射视频| 亚洲美女视频黄频| 99久久精品国产亚洲精品| 性少妇av在线| 成人国产av品久久久| 国产淫语在线视频| 欧美精品高潮呻吟av久久| 国产深夜福利视频在线观看| 韩国精品一区二区三区| 极品少妇高潮喷水抽搐| 飞空精品影院首页| 嫩草影视91久久| 成年av动漫网址| 男女午夜视频在线观看| 国产日韩欧美在线精品| 国产av一区二区精品久久| 久久久欧美国产精品| 嫩草影视91久久| 一个人免费看片子| 韩国精品一区二区三区| 国产野战对白在线观看| 美女中出高潮动态图| 大香蕉久久网| 亚洲七黄色美女视频| 久久久亚洲精品成人影院| 青草久久国产| 极品人妻少妇av视频| 十分钟在线观看高清视频www| 夫妻午夜视频| 最近中文字幕高清免费大全6| 丰满少妇做爰视频| 一本大道久久a久久精品| 亚洲欧美日韩另类电影网站| 熟妇人妻不卡中文字幕| 午夜影院在线不卡| 亚洲久久久国产精品| 亚洲欧美一区二区三区黑人| 18在线观看网站| 国产精品成人在线| 亚洲美女视频黄频| 久久久国产一区二区| 久久人人爽av亚洲精品天堂| 夜夜骑夜夜射夜夜干| 另类精品久久| 精品亚洲成国产av| 老汉色av国产亚洲站长工具| 日本黄色日本黄色录像| 国产精品久久久久久人妻精品电影 | avwww免费| 亚洲第一av免费看| tube8黄色片| 久久精品国产亚洲av涩爱| 女性生殖器流出的白浆| 日韩大码丰满熟妇| 一级片免费观看大全| 伦理电影大哥的女人| 免费观看av网站的网址| 亚洲一区二区三区欧美精品| 日韩不卡一区二区三区视频在线| 亚洲一级一片aⅴ在线观看| 国产人伦9x9x在线观看| 午夜福利视频精品| 久久人人爽人人片av| 极品少妇高潮喷水抽搐| 国产伦理片在线播放av一区| 黑丝袜美女国产一区| 丝袜人妻中文字幕| 女人精品久久久久毛片| 大片电影免费在线观看免费| 中文字幕高清在线视频| 精品少妇一区二区三区视频日本电影 | 久久 成人 亚洲| 女性被躁到高潮视频| 少妇 在线观看| 18禁观看日本| 巨乳人妻的诱惑在线观看| 国产欧美亚洲国产| 电影成人av| 国产国语露脸激情在线看| 欧美中文综合在线视频| 丰满少妇做爰视频| 狂野欧美激情性bbbbbb| 久久狼人影院| 国产在线免费精品| 国产精品 国内视频| 少妇人妻久久综合中文| 久久综合国产亚洲精品| 一级毛片我不卡| 韩国高清视频一区二区三区| 新久久久久国产一级毛片| 国产一卡二卡三卡精品 | 一本色道久久久久久精品综合| 精品少妇久久久久久888优播| 亚洲成人手机| 久久久久久免费高清国产稀缺| 亚洲人成网站在线观看播放| 国产av精品麻豆| 国产成人精品久久二区二区91 | 一边亲一边摸免费视频| 男男h啪啪无遮挡| 婷婷色av中文字幕| 国产av一区二区精品久久| 国产精品一区二区在线观看99| 免费女性裸体啪啪无遮挡网站| 欧美日韩成人在线一区二区| 又粗又硬又长又爽又黄的视频| 久久韩国三级中文字幕| 少妇 在线观看| 熟妇人妻不卡中文字幕| 蜜桃在线观看..| 中文字幕另类日韩欧美亚洲嫩草| 国产精品一国产av| av又黄又爽大尺度在线免费看| xxxhd国产人妻xxx| 一二三四在线观看免费中文在| 日日啪夜夜爽| 日韩 欧美 亚洲 中文字幕| 亚洲成色77777| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品国产一区二区精华液| 亚洲 欧美一区二区三区| 最近手机中文字幕大全| 亚洲国产成人一精品久久久| 考比视频在线观看| 国产黄色免费在线视频| 成人18禁高潮啪啪吃奶动态图| 亚洲国产最新在线播放| 三上悠亚av全集在线观看| 中文字幕人妻丝袜一区二区 | 我要看黄色一级片免费的| 国产欧美日韩一区二区三区在线| 成年动漫av网址| 免费高清在线观看日韩| 91国产中文字幕| 欧美精品av麻豆av| 亚洲精品国产av成人精品| 亚洲国产欧美日韩在线播放| av视频免费观看在线观看| 亚洲av日韩精品久久久久久密 | 欧美黑人精品巨大| 亚洲av电影在线进入| 欧美亚洲 丝袜 人妻 在线| 女人精品久久久久毛片| 91精品国产国语对白视频| 在现免费观看毛片| 悠悠久久av| 国产淫语在线视频| 叶爱在线成人免费视频播放| 国产成人精品久久二区二区91 | 丰满乱子伦码专区| 青青草视频在线视频观看| 日日摸夜夜添夜夜爱| 99热全是精品| 亚洲国产看品久久| 人体艺术视频欧美日本| 国产精品无大码| 一级毛片电影观看| 欧美国产精品va在线观看不卡| 熟妇人妻不卡中文字幕| 大香蕉久久网| 哪个播放器可以免费观看大片| 国产亚洲最大av| 9热在线视频观看99| 久久久国产欧美日韩av| 叶爱在线成人免费视频播放| 性少妇av在线| 大片免费播放器 马上看| 亚洲av综合色区一区| 久久久欧美国产精品| 欧美精品人与动牲交sv欧美| 一级片'在线观看视频| 亚洲精品国产av成人精品| 国产成人系列免费观看| 男的添女的下面高潮视频| 搡老乐熟女国产| 老司机影院成人| 国产精品香港三级国产av潘金莲 | 国产免费现黄频在线看| 欧美日韩一级在线毛片| 美国免费a级毛片| 在线观看人妻少妇| 久久女婷五月综合色啪小说| av福利片在线| 国产99久久九九免费精品| 悠悠久久av| 美女高潮到喷水免费观看| 国产黄频视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 91精品伊人久久大香线蕉| 国产爽快片一区二区三区| 亚洲精品美女久久av网站| 精品国产乱码久久久久久小说| 国产成人免费无遮挡视频| 伊人久久大香线蕉亚洲五| 黄色一级大片看看| 欧美人与性动交α欧美精品济南到| 啦啦啦视频在线资源免费观看| 超碰97精品在线观看| 99国产精品免费福利视频| 国产欧美日韩综合在线一区二区| 免费观看性生交大片5| 中文字幕亚洲精品专区| 国产精品二区激情视频| 一区在线观看完整版| 最近中文字幕高清免费大全6| 侵犯人妻中文字幕一二三四区| 午夜福利一区二区在线看| 欧美日韩成人在线一区二区| 日韩中文字幕视频在线看片| a 毛片基地| 亚洲国产精品999| 国产精品三级大全| 一级毛片我不卡| 久久狼人影院| 午夜av观看不卡| 一级毛片我不卡| 一区福利在线观看| 男女下面插进去视频免费观看| 大码成人一级视频| 在线 av 中文字幕| 亚洲国产毛片av蜜桃av| 亚洲第一区二区三区不卡| 色婷婷av一区二区三区视频| 一二三四中文在线观看免费高清| 麻豆乱淫一区二区| 宅男免费午夜| 久久天躁狠狠躁夜夜2o2o | 国产成人免费观看mmmm| av电影中文网址| 国产精品香港三级国产av潘金莲 | 日韩一本色道免费dvd| 欧美国产精品一级二级三级| 日韩中文字幕欧美一区二区 | 成年人免费黄色播放视频| 亚洲美女视频黄频| 一本久久精品| 成年动漫av网址| 国产精品成人在线| 啦啦啦在线免费观看视频4| 涩涩av久久男人的天堂| 51午夜福利影视在线观看| 欧美另类一区| 免费黄色在线免费观看| 五月开心婷婷网| 久久久精品国产亚洲av高清涩受| 午夜91福利影院| 精品国产一区二区三区四区第35| 1024视频免费在线观看| 欧美日韩视频高清一区二区三区二| 国产成人欧美| 亚洲久久久国产精品| 午夜老司机福利片| 亚洲欧美清纯卡通| 美女中出高潮动态图| 国产一区二区三区av在线| 久久国产亚洲av麻豆专区| 自线自在国产av| 国产av精品麻豆| 亚洲av福利一区| 亚洲欧美清纯卡通| 国产精品久久久av美女十八| 街头女战士在线观看网站| 91精品国产国语对白视频| 精品国产一区二区三区四区第35| 亚洲人成电影观看| 黄色怎么调成土黄色| av线在线观看网站| 大香蕉久久网| 爱豆传媒免费全集在线观看| 国产 精品1| 国产在线视频一区二区| 亚洲免费av在线视频| a 毛片基地| 久久久国产精品麻豆| 777米奇影视久久| 韩国高清视频一区二区三区| 国产在线视频一区二区| 少妇人妻久久综合中文| 我要看黄色一级片免费的| 少妇被粗大的猛进出69影院| 最黄视频免费看| 一本久久精品| 国产午夜精品一二区理论片| 亚洲美女黄色视频免费看| netflix在线观看网站| 午夜福利视频在线观看免费| 高清黄色对白视频在线免费看| 人人妻人人添人人爽欧美一区卜| 在线亚洲精品国产二区图片欧美| 大陆偷拍与自拍| 国产激情久久老熟女| 啦啦啦啦在线视频资源| 夫妻性生交免费视频一级片| 2018国产大陆天天弄谢| 免费黄色在线免费观看| 黄色一级大片看看| 亚洲精品乱久久久久久| 亚洲五月色婷婷综合| 国产在线免费精品| 国产精品人妻久久久影院| 18禁动态无遮挡网站| 你懂的网址亚洲精品在线观看| 搡老乐熟女国产| 高清在线视频一区二区三区| 在线看a的网站| 久久午夜综合久久蜜桃| 国产深夜福利视频在线观看| 欧美激情高清一区二区三区 | 国产野战对白在线观看| 视频在线观看一区二区三区| 欧美日韩国产mv在线观看视频| 少妇被粗大的猛进出69影院| 日韩不卡一区二区三区视频在线| 国产黄色视频一区二区在线观看| 女人久久www免费人成看片| 午夜久久久在线观看| 久久人人97超碰香蕉20202| av网站免费在线观看视频| 亚洲婷婷狠狠爱综合网| 欧美激情 高清一区二区三区| 啦啦啦在线观看免费高清www| 热re99久久精品国产66热6| 天堂中文最新版在线下载| 免费人妻精品一区二区三区视频| av在线观看视频网站免费| 超碰成人久久| 国产精品.久久久| 青春草亚洲视频在线观看| 免费黄频网站在线观看国产| 1024视频免费在线观看| 国产精品久久久久久精品电影小说| 亚洲国产日韩一区二区| 在线观看三级黄色| 超碰97精品在线观看| 蜜桃国产av成人99| 黄片小视频在线播放| 亚洲精品视频女| 国产精品欧美亚洲77777| 美女国产高潮福利片在线看| 国产成人系列免费观看| 男女下面插进去视频免费观看| 免费久久久久久久精品成人欧美视频| 午夜福利影视在线免费观看| 亚洲第一av免费看| 国产成人精品无人区| 天天操日日干夜夜撸| 又粗又硬又长又爽又黄的视频| av电影中文网址| 少妇精品久久久久久久| 丝袜美腿诱惑在线| 日韩中文字幕视频在线看片| 黄色视频在线播放观看不卡| 日韩欧美一区视频在线观看| 亚洲人成电影观看| 一区二区三区四区激情视频| 亚洲国产成人一精品久久久| videosex国产| 久久久久国产一级毛片高清牌| 晚上一个人看的免费电影| 美女大奶头黄色视频| 久久韩国三级中文字幕| 久久久欧美国产精品| 国产男人的电影天堂91| 黑人巨大精品欧美一区二区蜜桃| 天天操日日干夜夜撸| 精品少妇内射三级| 亚洲成av片中文字幕在线观看| 欧美精品av麻豆av| 美女中出高潮动态图| 日韩一卡2卡3卡4卡2021年| 久久久久久久大尺度免费视频| 一区二区三区精品91| 国产精品久久久久久精品电影小说| 国产精品av久久久久免费| 亚洲 欧美一区二区三区| 国产成人午夜福利电影在线观看| 国产男人的电影天堂91| 在线亚洲精品国产二区图片欧美| 国产精品无大码| 免费在线观看视频国产中文字幕亚洲 | 欧美日韩一区二区视频在线观看视频在线| 色播在线永久视频| 大香蕉久久网| 国产亚洲最大av| 亚洲精品久久成人aⅴ小说| 桃花免费在线播放| 久久人人97超碰香蕉20202| 日本午夜av视频| 日韩一区二区三区影片| 最近最新中文字幕大全免费视频 | 精品国产国语对白av| 国产欧美亚洲国产| 观看美女的网站| 国产不卡av网站在线观看| 久久久国产一区二区| 一个人免费看片子| 国产成人免费观看mmmm| 久久久久久人妻| 最黄视频免费看| 999精品在线视频| 91国产中文字幕| 侵犯人妻中文字幕一二三四区| 男女下面插进去视频免费观看| 国产在视频线精品| 天天添夜夜摸| 啦啦啦视频在线资源免费观看| 99re6热这里在线精品视频| 欧美日韩视频高清一区二区三区二| 尾随美女入室| 久久午夜综合久久蜜桃| 热99久久久久精品小说推荐| 国产片内射在线| h视频一区二区三区| av片东京热男人的天堂| av线在线观看网站| 日韩不卡一区二区三区视频在线| 亚洲国产欧美日韩在线播放| 亚洲美女黄色视频免费看| 在线观看www视频免费| 国产黄频视频在线观看| 日韩一区二区三区影片| 色精品久久人妻99蜜桃| 久久人人爽人人片av| 最新的欧美精品一区二区| 韩国av在线不卡| 最新的欧美精品一区二区| 亚洲精品第二区| 麻豆av在线久日| 一二三四在线观看免费中文在| 国产精品久久久久久人妻精品电影 | 欧美另类一区| 男人添女人高潮全过程视频| 国产成人欧美在线观看 | 一二三四中文在线观看免费高清| 19禁男女啪啪无遮挡网站| 最近中文字幕高清免费大全6| 欧美激情极品国产一区二区三区| 久久久久久久精品精品| 看免费成人av毛片| 91精品国产国语对白视频| www.av在线官网国产| 亚洲成人手机| 国产高清不卡午夜福利| 午夜91福利影院| 成年美女黄网站色视频大全免费| 一本大道久久a久久精品| 狠狠婷婷综合久久久久久88av| 美女福利国产在线| 最近中文字幕2019免费版| 亚洲欧美一区二区三区久久| www.av在线官网国产| 久久韩国三级中文字幕| 国产av精品麻豆| 亚洲精品国产av成人精品| 香蕉丝袜av| 亚洲综合精品二区| 亚洲七黄色美女视频| 黄色视频在线播放观看不卡| 男女边摸边吃奶| 韩国av在线不卡| 国产精品偷伦视频观看了| 亚洲av综合色区一区| 中文欧美无线码| 亚洲精品中文字幕在线视频| 国产 一区精品| 午夜精品国产一区二区电影| 亚洲国产av影院在线观看| 国产乱人偷精品视频| 操美女的视频在线观看| 日韩 欧美 亚洲 中文字幕| 中文字幕最新亚洲高清| 不卡视频在线观看欧美| 欧美日韩视频高清一区二区三区二| 高清视频免费观看一区二区| 亚洲情色 制服丝袜| 看免费成人av毛片| 天天躁日日躁夜夜躁夜夜| 亚洲国产精品一区二区三区在线| 久久久久国产精品人妻一区二区| 大香蕉久久成人网| 女的被弄到高潮叫床怎么办| 日韩av不卡免费在线播放| 韩国av在线不卡| 日本午夜av视频| 嫩草影视91久久| 可以免费在线观看a视频的电影网站 | 日韩一区二区视频免费看| 一本—道久久a久久精品蜜桃钙片| 久久国产精品男人的天堂亚洲|