• <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真有效性的標量化
    基于多級小波域變換的時域擴散熒光層析成像方法
    69av精品久久久久久| 午夜福利高清视频| 18禁黄网站禁片免费观看直播| 色综合站精品国产| 日韩 亚洲 欧美在线| 欧美一级a爱片免费观看看| 免费看光身美女| 成人国产一区最新在线观看| 国产精品日韩av在线免费观看| 亚洲五月婷婷丁香| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品一区二区免费观看| 99国产精品一区二区蜜桃av| 国产精品久久久久久久久免 | 日本与韩国留学比较| 国产亚洲欧美在线一区二区| 国产精品嫩草影院av在线观看 | 真实男女啪啪啪动态图| 国产精品亚洲一级av第二区| 美女高潮喷水抽搐中文字幕| 亚洲国产色片| 免费观看精品视频网站| 欧美日韩黄片免| 丰满人妻一区二区三区视频av| 一个人观看的视频www高清免费观看| 国产成人欧美在线观看| 欧美日韩瑟瑟在线播放| 日本 av在线| 日本a在线网址| 国产视频内射| 99久久99久久久精品蜜桃| 日本免费一区二区三区高清不卡| 欧美一区二区国产精品久久精品| 精品熟女少妇八av免费久了| 欧美日本视频| 久久热精品热| 中文字幕熟女人妻在线| 又爽又黄无遮挡网站| 国产伦在线观看视频一区| 日本五十路高清| 亚洲av电影在线进入| 一进一出抽搐动态| 丁香六月欧美| 自拍偷自拍亚洲精品老妇| 亚洲狠狠婷婷综合久久图片| 成年女人毛片免费观看观看9| 一本综合久久免费| 一a级毛片在线观看| 日本a在线网址| 蜜桃亚洲精品一区二区三区| 丁香六月欧美| av天堂在线播放| 欧美日韩中文字幕国产精品一区二区三区| 很黄的视频免费| 18+在线观看网站| 身体一侧抽搐| 99在线人妻在线中文字幕| 亚洲国产色片| 色5月婷婷丁香| 免费看a级黄色片| 黄色视频,在线免费观看| 亚洲欧美激情综合另类| 一个人看视频在线观看www免费| 亚洲国产精品sss在线观看| 美女xxoo啪啪120秒动态图 | 国产精品一及| 国产真实乱freesex| 国产av麻豆久久久久久久| 九九在线视频观看精品| 黄色配什么色好看| 狠狠狠狠99中文字幕| 亚洲av电影在线进入| 亚洲精品影视一区二区三区av| 一进一出抽搐动态| 国产单亲对白刺激| 成人鲁丝片一二三区免费| 在线观看午夜福利视频| 国产麻豆成人av免费视频| 亚洲avbb在线观看| 黄色丝袜av网址大全| 人妻久久中文字幕网| 欧美乱色亚洲激情| 悠悠久久av| 少妇丰满av| 一进一出好大好爽视频| 嫩草影院精品99| 51午夜福利影视在线观看| 久久九九热精品免费| 精品久久久久久成人av| 国内精品一区二区在线观看| 色在线成人网| 亚洲狠狠婷婷综合久久图片| 真人一进一出gif抽搐免费| 国产单亲对白刺激| 精品久久久久久久人妻蜜臀av| 精品不卡国产一区二区三区| 精华霜和精华液先用哪个| 中文字幕熟女人妻在线| 性插视频无遮挡在线免费观看| 国产人妻一区二区三区在| 亚洲三级黄色毛片| 亚洲成人免费电影在线观看| 九九久久精品国产亚洲av麻豆| 人人妻人人澡欧美一区二区| www.999成人在线观看| 性插视频无遮挡在线免费观看| 桃红色精品国产亚洲av| 一级黄片播放器| 欧洲精品卡2卡3卡4卡5卡区| 日韩大尺度精品在线看网址| 在线播放无遮挡| 波野结衣二区三区在线| 两个人的视频大全免费| 搡老岳熟女国产| 国产成年人精品一区二区| or卡值多少钱| 美女cb高潮喷水在线观看| 中文字幕熟女人妻在线| 久久久久久久精品吃奶| 免费在线观看日本一区| 成人三级黄色视频| 18禁裸乳无遮挡免费网站照片| 亚洲在线观看片| 日本免费一区二区三区高清不卡| 深爱激情五月婷婷| 亚洲第一电影网av| 日本黄色视频三级网站网址| 国产精品日韩av在线免费观看| 久久99热6这里只有精品| 黄色一级大片看看| 亚洲内射少妇av| 亚洲人成伊人成综合网2020| 在线观看66精品国产| 深夜精品福利| 国内精品久久久久久久电影| 久久久久久久久久成人| 欧美黑人欧美精品刺激| 狂野欧美白嫩少妇大欣赏| 欧美激情国产日韩精品一区| 色吧在线观看| 日本撒尿小便嘘嘘汇集6| 99国产精品一区二区三区| 久久99热这里只有精品18| 久久九九热精品免费| 麻豆一二三区av精品| 国产真实伦视频高清在线观看 | 亚洲五月婷婷丁香| 亚洲一区二区三区不卡视频| 五月伊人婷婷丁香| 简卡轻食公司| 两性午夜刺激爽爽歪歪视频在线观看| 欧美黄色片欧美黄色片| 国产单亲对白刺激| 1000部很黄的大片| 亚洲黑人精品在线| 国产色婷婷99| 在线观看一区二区三区| 51国产日韩欧美| 免费看光身美女| 国产精品久久久久久精品电影| 欧美色视频一区免费| 18+在线观看网站| 中文字幕人成人乱码亚洲影| 成人av在线播放网站| 韩国av一区二区三区四区| 一a级毛片在线观看| 久久九九热精品免费| 亚洲午夜理论影院| 国产高清激情床上av| 国产欧美日韩一区二区精品| 伊人久久精品亚洲午夜| 国产一区二区在线av高清观看| 99热这里只有是精品在线观看 | 久久午夜亚洲精品久久| 特级一级黄色大片| 亚洲18禁久久av| 夜夜夜夜夜久久久久| 波野结衣二区三区在线| av在线天堂中文字幕| 国内久久婷婷六月综合欲色啪| 1000部很黄的大片| 久久亚洲真实| 99久久99久久久精品蜜桃| 精品国产三级普通话版| 免费在线观看亚洲国产| 欧美三级亚洲精品| 欧美成人性av电影在线观看| 亚洲av中文字字幕乱码综合| 国产主播在线观看一区二区| 日本撒尿小便嘘嘘汇集6| 黄色日韩在线| 久久精品综合一区二区三区| 别揉我奶头 嗯啊视频| 1000部很黄的大片| 国产精品久久视频播放| 国产伦在线观看视频一区| 高清日韩中文字幕在线| 国内精品久久久久久久电影| 亚洲第一区二区三区不卡| 极品教师在线免费播放| 免费观看精品视频网站| 3wmmmm亚洲av在线观看| 亚洲人成伊人成综合网2020| 成人特级黄色片久久久久久久| 久久久精品欧美日韩精品| 精品久久久久久久久久久久久| 成人国产一区最新在线观看| 久久性视频一级片| 网址你懂的国产日韩在线| 国产亚洲欧美在线一区二区| 国产成人啪精品午夜网站| 国产乱人视频| 午夜免费激情av| 一本综合久久免费| 长腿黑丝高跟| 国产aⅴ精品一区二区三区波| av天堂在线播放| 国模一区二区三区四区视频| 18+在线观看网站| 色综合亚洲欧美另类图片| 热99re8久久精品国产| 亚洲精品一卡2卡三卡4卡5卡| 麻豆国产97在线/欧美| 亚洲久久久久久中文字幕| 精品人妻偷拍中文字幕| 国模一区二区三区四区视频| 国产午夜精品论理片| 亚州av有码| 欧美丝袜亚洲另类 | 99久久精品热视频| 欧美3d第一页| 一进一出抽搐gif免费好疼| 岛国在线免费视频观看| 成人一区二区视频在线观看| 久久久久久久久大av| 亚洲成a人片在线一区二区| 小说图片视频综合网站| 听说在线观看完整版免费高清| 中文字幕人妻熟人妻熟丝袜美| 国产野战对白在线观看| 可以在线观看的亚洲视频| 国产黄片美女视频| 制服丝袜大香蕉在线| 午夜精品一区二区三区免费看| 精品久久久久久久久亚洲 | 国产精品野战在线观看| 99热只有精品国产| 亚洲av电影不卡..在线观看| 国产精品久久久久久亚洲av鲁大| 欧美成狂野欧美在线观看| 757午夜福利合集在线观看| 国产黄a三级三级三级人| 国产久久久一区二区三区| av天堂在线播放| 给我免费播放毛片高清在线观看| 国内精品一区二区在线观看| 国产精品爽爽va在线观看网站| 丰满人妻熟妇乱又伦精品不卡| av中文乱码字幕在线| 亚洲最大成人手机在线| 小说图片视频综合网站| 久久久久久久精品吃奶| 搞女人的毛片| 国产一区二区三区在线臀色熟女| 国产aⅴ精品一区二区三区波| 午夜福利18| av女优亚洲男人天堂| 国产精品1区2区在线观看.| 欧美性感艳星| 国产av麻豆久久久久久久| 亚洲综合色惰| 亚洲三级黄色毛片| 一a级毛片在线观看| 非洲黑人性xxxx精品又粗又长| 小说图片视频综合网站| 最近中文字幕高清免费大全6 | 国产熟女xx| 久久中文看片网| 精品午夜福利视频在线观看一区| 国产精品影院久久| 亚洲人与动物交配视频| 不卡一级毛片| 久久性视频一级片| 波野结衣二区三区在线| 欧美成人性av电影在线观看| 亚洲成av人片在线播放无| 欧美精品国产亚洲| 老女人水多毛片| 色噜噜av男人的天堂激情| 在线看三级毛片| av在线蜜桃| 成人无遮挡网站| 人妻制服诱惑在线中文字幕| 亚洲狠狠婷婷综合久久图片| 午夜精品一区二区三区免费看| 久久精品夜夜夜夜夜久久蜜豆| 99久久精品国产亚洲精品| 每晚都被弄得嗷嗷叫到高潮| 久久久色成人| 在线天堂最新版资源| 亚洲国产精品成人综合色| 亚洲无线观看免费| or卡值多少钱| 欧美高清成人免费视频www| 日本精品一区二区三区蜜桃| 色尼玛亚洲综合影院| 亚洲精品日韩av片在线观看| 成人鲁丝片一二三区免费| 久久精品人妻少妇| 欧美日韩福利视频一区二区| 麻豆久久精品国产亚洲av| 国模一区二区三区四区视频| .国产精品久久| 国产高清有码在线观看视频| 国产成+人综合+亚洲专区| 久久久久久久午夜电影| 成人美女网站在线观看视频| 免费电影在线观看免费观看| 亚洲国产精品久久男人天堂| 一级a爱片免费观看的视频| aaaaa片日本免费| 成人毛片a级毛片在线播放| 最后的刺客免费高清国语| 亚洲色图av天堂| а√天堂www在线а√下载| 色5月婷婷丁香| or卡值多少钱| 91九色精品人成在线观看| 久久久久久久午夜电影| 色综合婷婷激情| 国产精品三级大全| 少妇的逼水好多| 国产精品三级大全| 中文字幕av在线有码专区| 精品久久国产蜜桃| 午夜免费男女啪啪视频观看 | 久久久久久久久久黄片| 亚洲在线观看片| 久久久久久久久久黄片| 欧美一级a爱片免费观看看| 美女黄网站色视频| 怎么达到女性高潮| 国产精品三级大全| 国产精品人妻久久久久久| 欧美一区二区国产精品久久精品| 看黄色毛片网站| 国产成人福利小说| 中文字幕高清在线视频| 久久久精品欧美日韩精品| 中文字幕高清在线视频| 乱人视频在线观看| 内射极品少妇av片p| 亚洲av一区综合| 欧美+日韩+精品| 超碰av人人做人人爽久久| 色5月婷婷丁香| 18禁在线播放成人免费| 亚洲va日本ⅴa欧美va伊人久久| 少妇裸体淫交视频免费看高清| 欧美日韩国产亚洲二区| 深爱激情五月婷婷| 人妻制服诱惑在线中文字幕| 欧美成人a在线观看| 国语自产精品视频在线第100页| 免费搜索国产男女视频| 国产毛片a区久久久久| 欧美激情在线99| 在线十欧美十亚洲十日本专区| 搡老妇女老女人老熟妇| 国内精品一区二区在线观看| 中亚洲国语对白在线视频| av福利片在线观看| 日韩欧美精品免费久久 | 午夜激情欧美在线| 精品人妻1区二区| 国产视频内射| 五月伊人婷婷丁香| 欧美国产日韩亚洲一区| 久久久久久久久久黄片| 可以在线观看毛片的网站| 免费高清视频大片| 看片在线看免费视频| 女生性感内裤真人,穿戴方法视频| 国产伦精品一区二区三区四那| 美女黄网站色视频| 极品教师在线视频| 97碰自拍视频| 亚洲aⅴ乱码一区二区在线播放| 国产aⅴ精品一区二区三区波| 麻豆久久精品国产亚洲av| 中国美女看黄片| 久久久久久久久中文| av欧美777| 51国产日韩欧美| 12—13女人毛片做爰片一| 午夜免费激情av| 日日摸夜夜添夜夜添av毛片 | 制服丝袜大香蕉在线| 精品久久久久久久人妻蜜臀av| 亚洲av免费在线观看| 免费黄网站久久成人精品 | 国产精品亚洲av一区麻豆| 身体一侧抽搐| 国产成人影院久久av| 久久久久国产精品人妻aⅴ院| 噜噜噜噜噜久久久久久91| 国产精品野战在线观看| 搡女人真爽免费视频火全软件 | 欧美不卡视频在线免费观看| 色综合站精品国产| av黄色大香蕉| 国产真实伦视频高清在线观看 | 热99re8久久精品国产| 夜夜看夜夜爽夜夜摸| 欧美激情在线99| 9191精品国产免费久久| 日本精品一区二区三区蜜桃| 欧美另类亚洲清纯唯美| 大型黄色视频在线免费观看| 一本久久中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产欧美人成| www.熟女人妻精品国产| 久久人人精品亚洲av| 国产真实乱freesex| 12—13女人毛片做爰片一| 亚洲人成电影免费在线| 国内少妇人妻偷人精品xxx网站| 很黄的视频免费| 亚洲精品日韩av片在线观看| 色尼玛亚洲综合影院| 免费电影在线观看免费观看| 99久久精品一区二区三区| 国产精品爽爽va在线观看网站| 特级一级黄色大片| 18美女黄网站色大片免费观看| 久久香蕉精品热| 久久久国产成人免费| 中文字幕av成人在线电影| 色视频www国产| 3wmmmm亚洲av在线观看| 一本精品99久久精品77| 久久精品综合一区二区三区| 国产精品99久久久久久久久| 中国美女看黄片| 亚洲精品在线美女| aaaaa片日本免费| 亚洲国产精品999在线| 亚洲国产精品久久男人天堂| 日本精品一区二区三区蜜桃| 午夜免费男女啪啪视频观看 | 日韩欧美精品v在线| 91麻豆av在线| 国产不卡一卡二| 最近最新中文字幕大全电影3| 成年女人毛片免费观看观看9| av中文乱码字幕在线| 亚洲av不卡在线观看| 国产伦精品一区二区三区视频9| 国产又黄又爽又无遮挡在线| 久久久久久久午夜电影| 国产伦在线观看视频一区| 欧美黑人巨大hd| 精品国产亚洲在线| av女优亚洲男人天堂| 欧美日本亚洲视频在线播放| 91麻豆av在线| 免费电影在线观看免费观看| 免费观看精品视频网站| 国内毛片毛片毛片毛片毛片| 黄色一级大片看看| 男人的好看免费观看在线视频| 在线观看av片永久免费下载| 桃色一区二区三区在线观看| 男女做爰动态图高潮gif福利片| 久久国产精品人妻蜜桃| 国产精品日韩av在线免费观看| 欧美绝顶高潮抽搐喷水| 亚洲精品影视一区二区三区av| 精品一区二区三区人妻视频| 国产日本99.免费观看| 少妇人妻一区二区三区视频| 天堂影院成人在线观看| 亚洲熟妇熟女久久| 99国产精品一区二区三区| 日本 欧美在线| 国产高清视频在线播放一区| 99久久精品国产亚洲精品| 一级a爱片免费观看的视频| 一本精品99久久精品77| 成年女人永久免费观看视频| 久久久久国产精品人妻aⅴ院| 亚洲精品久久国产高清桃花| 日韩免费av在线播放| 亚洲第一欧美日韩一区二区三区| 91麻豆精品激情在线观看国产| 国产在线男女| 亚洲av熟女| 成人特级黄色片久久久久久久| 亚洲人与动物交配视频| 3wmmmm亚洲av在线观看| 日韩大尺度精品在线看网址| 欧美性猛交╳xxx乱大交人| .国产精品久久| 1024手机看黄色片| 18禁黄网站禁片免费观看直播| 国产大屁股一区二区在线视频| 久久精品综合一区二区三区| 白带黄色成豆腐渣| 日韩亚洲欧美综合| 97人妻精品一区二区三区麻豆| 国产69精品久久久久777片| 婷婷精品国产亚洲av在线| 老鸭窝网址在线观看| 午夜精品久久久久久毛片777| 网址你懂的国产日韩在线| 国内毛片毛片毛片毛片毛片| 国产精品电影一区二区三区| 亚州av有码| 又爽又黄a免费视频| 国产不卡一卡二| 国产日本99.免费观看| 1000部很黄的大片| av女优亚洲男人天堂| 动漫黄色视频在线观看| 97热精品久久久久久| 国产精品1区2区在线观看.| 女人十人毛片免费观看3o分钟| 国产精品不卡视频一区二区 | 欧美日韩亚洲国产一区二区在线观看| 午夜影院日韩av| 国产69精品久久久久777片| 欧美又色又爽又黄视频| 可以在线观看的亚洲视频| av天堂在线播放| 99热这里只有精品一区| 男女做爰动态图高潮gif福利片| 国内久久婷婷六月综合欲色啪| 精品久久久久久久久av| 国产探花极品一区二区| 99国产精品一区二区蜜桃av| 国产老妇女一区| 最近在线观看免费完整版| 丁香欧美五月| 国产探花在线观看一区二区| 国产伦人伦偷精品视频| 免费在线观看影片大全网站| 蜜桃亚洲精品一区二区三区| 亚洲内射少妇av| 国内精品一区二区在线观看| 可以在线观看的亚洲视频| 午夜福利欧美成人| 成人精品一区二区免费| 毛片一级片免费看久久久久 | 五月伊人婷婷丁香| 国产探花在线观看一区二区| 狠狠狠狠99中文字幕| ponron亚洲| 精品99又大又爽又粗少妇毛片 | 欧美成人a在线观看| 国产精品亚洲一级av第二区| 禁无遮挡网站| 老熟妇仑乱视频hdxx| 亚洲人与动物交配视频| 久久久久久九九精品二区国产| 一进一出抽搐动态| 国产精品影院久久| 免费在线观看成人毛片| 日本五十路高清| 首页视频小说图片口味搜索| 18禁黄网站禁片免费观看直播| 精品久久久久久成人av| 亚洲av.av天堂| 久久国产精品人妻蜜桃| 男女床上黄色一级片免费看| av专区在线播放| 亚洲狠狠婷婷综合久久图片| 国产伦在线观看视频一区| 欧美激情在线99| 又紧又爽又黄一区二区| 亚洲五月婷婷丁香| 成人毛片a级毛片在线播放| 精品无人区乱码1区二区| 精品人妻1区二区| 人妻久久中文字幕网| 97超视频在线观看视频| 欧美+日韩+精品| 给我免费播放毛片高清在线观看| 午夜福利欧美成人| 国产伦精品一区二区三区四那| 亚洲精品久久国产高清桃花| 人人妻人人澡欧美一区二区| 免费观看精品视频网站| 久久久久久久精品吃奶| 午夜免费成人在线视频| 真实男女啪啪啪动态图| 91久久精品国产一区二区成人| 中文字幕av成人在线电影| 老女人水多毛片| 国产黄色小视频在线观看| 大型黄色视频在线免费观看| 国产av一区在线观看免费| 一个人免费在线观看电影| 人妻夜夜爽99麻豆av| 在线播放国产精品三级| av在线天堂中文字幕| 成人av一区二区三区在线看| 蜜桃久久精品国产亚洲av| 免费大片18禁| 我的老师免费观看完整版| 尤物成人国产欧美一区二区三区| 18禁在线播放成人免费| ponron亚洲| av黄色大香蕉| 51午夜福利影视在线观看|