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

    橢圓極化波對(duì)大鼠電磁暴露的數(shù)值模擬計(jì)算

    2020-04-28 09:59:40李孟達(dá)王向暉齊紅新
    關(guān)鍵詞:線極化波源電磁

    李孟達(dá) 王向暉 張 杰 齊紅新

    (華東師范大學(xué)物理與電子科學(xué)學(xué)院生物物理實(shí)驗(yàn)室 上海200241)

    電磁暴露引起的生物效應(yīng)與受輻照動(dòng)物在電磁場(chǎng)中的電磁暴露劑量有關(guān),電磁暴露劑量數(shù)值模擬計(jì)算和實(shí)驗(yàn)探究是電磁生物效應(yīng)研究中重要的一環(huán)。在動(dòng)物電磁暴露實(shí)驗(yàn)中,電磁波本身的不同特性(如頻率、入射方式、電場(chǎng)極化方向等)以及電磁波與受輻照體之間相對(duì)位置的變化,會(huì)使被輻照體所受到的電磁暴露劑量大?。ū任章剩⊿AR))發(fā)生改變,所以,考察波源特性在動(dòng)物電磁暴露實(shí)驗(yàn)中對(duì)SAR(包括局部SAR 和全身平均比吸收率(WBASAR))值的影響,對(duì)電磁暴露劑量的評(píng)估十分必要。同時(shí),隨著電磁暴露數(shù)值模擬方面的快速發(fā)展,如在算法運(yùn)算效率以及提升精度等方面研究[1-2]、對(duì)電磁暴露劑量防護(hù)的模擬計(jì)算[3-5],以及將電磁暴露數(shù)值模擬應(yīng)用在高鐵、純電動(dòng)汽車等[6-8]實(shí)際情景的輻射劑量探究,促使數(shù)值模擬方法在探究動(dòng)物電磁暴露劑量方面做出了許多工作。Wu 等[9]和Bamba 等[10]關(guān)于20 MHz~6 GHz頻段線極化波的SAR計(jì)算結(jié)果表明,波源的頻率是影響電磁暴露劑量大小的重要因素;Conil等[11]通過(guò)改變線極化波源與計(jì)算模型之間的方位角,分析不同入射角度的WBASAR和局部SAR 結(jié)果,發(fā)現(xiàn)線極化波入射角度同樣會(huì)對(duì)電磁暴露劑量產(chǎn)生影響。Chen 等[12]和Yelkenci等[13]綜合考慮線極化波源的頻率、入射角度以及極化方向3種因素進(jìn)行大量的數(shù)值計(jì)算,其結(jié)果進(jìn)一步說(shuō)明波源入射方式是影響電磁暴露劑量的重要因素。

    以上研究均以線極化波作為波源,而非線性極化(圓極化)電磁波同樣在無(wú)線訊通等實(shí)際工程中有著廣泛的應(yīng)用,如有些通訊設(shè)備被專門設(shè)計(jì)用于接收和發(fā)射圓極化、橢圓極化等非線極化波。因此,對(duì)非線極化波條件下的電磁暴露劑量評(píng)估同樣非常重要。Guy 和Chou 等[14-15]以圓形波導(dǎo)管中產(chǎn)生的圓極化波作為波源對(duì)大鼠進(jìn)行電磁輻照,首先開展了非線極化波對(duì)動(dòng)物電磁暴露研究;Wasoontarajaroen等[16]以圓形波導(dǎo)管中產(chǎn)生的圓極化波對(duì)大鼠進(jìn)行更進(jìn)一步的電磁輻照實(shí)驗(yàn),其結(jié)果表明圓極化波源的WBASAR值同樣與波源入射方向存在一定關(guān)系;此外,不同特性的電磁波對(duì)于細(xì)胞的生存能力也會(huì)產(chǎn)生不同的影響[17];Shckorbatov 等[18-19]從細(xì)胞的生物效應(yīng)角度證明了不同旋向的橢圓極化波會(huì)對(duì)人類成纖維細(xì)胞的遺傳物質(zhì)結(jié)構(gòu)產(chǎn)生不同的影響,結(jié)果還證明了旋向不同的圓極化波對(duì)人類細(xì)胞的損傷差異性極小。綜上,不同特性的非線極化波在不同入射條件下,會(huì)對(duì)細(xì)胞的生物效應(yīng)產(chǎn)生不同的影響。

    對(duì)于動(dòng)物電磁暴露劑量的評(píng)估,以非線極化電磁波或線極化波作為輻照波源時(shí)會(huì)呈現(xiàn)何種分布規(guī)律?對(duì)動(dòng)物電磁暴露劑量的評(píng)估有何不同?以及線極化波源是否足夠涵蓋不同極化方式的電磁波所產(chǎn)生的電磁輻照影響?目前在這些方面的研究還較少。本文通過(guò)自主建立雙波源的電磁暴露數(shù)值模擬系統(tǒng),以不同形態(tài)的橢圓極化波作為波源,在不同的入射方向上對(duì)大鼠進(jìn)行電磁暴露數(shù)值模擬計(jì)算,并與線極化波的WBASAR結(jié)果進(jìn)行比較,探討橢圓極化波電磁暴露的特性,為動(dòng)物電磁暴露的劑量評(píng)估提供方法和數(shù)據(jù)上的補(bǔ)充。

    1 計(jì)算方法與原理

    1.1 數(shù)值方法

    數(shù)值模擬計(jì)算基于時(shí)域有限差分法(Finitedifference time-domain,F(xiàn)DTD)[20-22],F(xiàn)DTD 方法是求解Maxwell微分方程的直接時(shí)域法,通過(guò)對(duì)電場(chǎng)分量E和磁場(chǎng)分量H在空間和時(shí)間上交替抽樣的離散,將Maxwell方程中的微分項(xiàng)進(jìn)行差分處理得到電場(chǎng)和磁場(chǎng)在空間上的時(shí)域遞推公式。模擬電磁場(chǎng)隨時(shí)間演化的過(guò)程可以很好地處理復(fù)雜物體的電磁散射和輻射問題。Maxwell 旋度方程見式(1)。

    式中:E、H 分別為電場(chǎng)和磁場(chǎng)強(qiáng)度;D、B 分別為電通量密度和磁通量密度;J、Jm分別為電流密度和磁流密度。文中所有變量的單位均為國(guó)際標(biāo)準(zhǔn)單位。在三維直角坐標(biāo)系的計(jì)算空間中,電場(chǎng)和磁場(chǎng)各自有3 個(gè)方向上的分量,故可以將Maxwell方程式寫為式(2)、(3)。

    式中的本構(gòu)關(guān)系可表示為式(4)。

    式中:ε、μ分別表示介電常數(shù)和磁導(dǎo)系數(shù),σ、σm分別表示電導(dǎo)率和導(dǎo)磁率。

    對(duì)式(2)、(3)進(jìn)行FDTD差分離散,即將式中的微分項(xiàng)利用差分項(xiàng)代替,以式(2)中第一分式為例[23],見式(5)。

    式中:Hy、Hz分別為磁場(chǎng)的y 方向分量和z 方向分量,對(duì)Maxwell方程中的其他分量的遞推式也做類似處理。此外,在計(jì)算空間邊界處利用各向異性完美匹配層(UPML)設(shè)置吸收邊界[23]模擬電磁波在無(wú)限大空間中的傳播。

    1.2 波源

    1.2.1 線極化波源

    選取不同極化方向的線極化波和不同形態(tài)的橢圓極化波作為入射波源進(jìn)行模擬計(jì)算。輻照波源為單波源線極化波時(shí),選取水平極化和垂直極化兩種線極化波作為波源,即極化角α分別為0和π/2,如圖1(a)所示。

    以+y 方向?yàn)槔?,用EzHxKy來(lái)代表沿著+y 方向入射,電場(chǎng)極化方向?yàn)閦 軸方向,而磁場(chǎng)分量沿x軸方向振動(dòng)的線極化波;同樣,以ExHzKy代表沿+y方向入射的線極化波,極化方向?yàn)閤軸。兩種線極化波的波矢方向以及電磁場(chǎng)分量的方向關(guān)系示意圖見圖1(b)、(c)。同理,以EyHzKx、EzHyKx表示沿x 軸方向入射的兩線極化波波源,以ExHyKz、EyHxKz表示沿z軸方向入射的兩線極化波波源。

    1.2.2 橢圓極化波源

    文中使用的橢圓極化波源利用兩個(gè)線極化波合成[24]。將極化方向相互垂直且頻率都相同的兩列線極化電磁波傳播至計(jì)算空間中,初始相位差?φ 從0~2π 每次間隔π/6 進(jìn)行變化。?φ 為0 和π 時(shí)得到兩極化角α 分別為π/4 和3π/4 的線極化波,其電場(chǎng)振動(dòng)方向如圖2(a)、(n)、(g)、(h)所示。?φ為π/2 和3π/2 時(shí)得到圓極化波分別如圖2(d)和(k)所示,其他情況則得到不同形態(tài)的橢圓極化波。

    在比較單波源線極化波的電磁暴露結(jié)果與不同形態(tài)橢圓極化波的結(jié)果時(shí),保持兩種電磁波的入射功率一致。線極化波源入射波電場(chǎng)幅值設(shè)置為1 V/m,橢圓極化波的每個(gè)單波源的振幅均設(shè)置為0.707 V/m[25],即兩者的輻射功率密度皆為0.001 33 W/m2。

    1.3 模型場(chǎng)景

    數(shù)值模擬計(jì)算中使用的電磁輻照對(duì)象是198 g大鼠的三維立體模型,購(gòu)買自IT’IS4公司。此大鼠模型根據(jù)解剖學(xué)結(jié)構(gòu)由60 種不同的組織構(gòu)建而成,其分辨率為1 mm,各個(gè)組織的電磁參數(shù)取自聯(lián)邦通信委員會(huì)的大鼠組織介電參數(shù)數(shù)據(jù)庫(kù)[26],表1列出了在輻照波源為2.45 GHz條件下重要組織的電導(dǎo)率與介電常數(shù)。

    在計(jì)算空間中,大鼠模型的位置關(guān)系如圖3(a)所示,波源入射的方向與計(jì)算空間坐標(biāo)系的關(guān)系如圖3(b)所示,其中φ表示波矢量與x軸的夾角,θ為波矢量與z軸的夾角。

    考慮以2.45 GHz的電磁波分別從頭部、尾部、背部和腹部以及側(cè)面等方向輻照大鼠模型,如圖3所示。大鼠模型的長(zhǎng)軸與z軸平行,體寬與體高分布分別沿x 軸和y 軸方向,頭部、尾部、背部和腹部入射方向?qū)?yīng)計(jì)算空間坐標(biāo)系的+z、-z、+y和-y 方向,側(cè)面入射則用+x/-x 表示。為了保證計(jì)算的精確度,綜合考慮入射波的頻率以及大鼠模型的電磁參數(shù)等因素,設(shè)定計(jì)算空間中FDTD元胞的網(wǎng)格邊長(zhǎng)為1 mm。

    在電磁暴露的研究中,常以SAR 值作為評(píng)價(jià)電磁暴露劑量大小的物理量[27-29]。SAR 值包括局部SAR 和WBASAR 計(jì)算方式,本文以大鼠的WBASAR 值作為電磁暴露水平的衡量標(biāo)準(zhǔn),觀察大鼠整體的能量吸收大小,其計(jì)算見公式(7)[30]。

    表1 大鼠模型中重要組織在2.45 GHz條件下對(duì)應(yīng)的電導(dǎo)率與介電常數(shù)Table 1 Electrical conductivity and permittivity of important tissues in a rat model under 2.45 GHz electromagnetic wave irradiation

    1.4 雙波源數(shù)值仿真的實(shí)現(xiàn)方法

    對(duì)于生物電磁劑量相關(guān)的數(shù)值仿真計(jì)算,Remcom公司開發(fā)的XFdtd 是該領(lǐng)域最常用的商用計(jì)算軟件,其高效穩(wěn)定的特點(diǎn)在生物電磁劑量研究中得到廣泛應(yīng)用。其不足之處是僅能在單一波源輻照條件下進(jìn)行電磁暴露劑量的計(jì)算,無(wú)法滿足橢圓極化波的產(chǎn)生條件,即需要極化方向相互垂直,且頻率都相同的兩列線極化電磁波以一定的相位差傳入計(jì)算空間。為了解決這一關(guān)鍵性難題,本文基于FDTD理論自主建立了可支持多波源入射的電磁暴露數(shù)值仿真系統(tǒng),以此來(lái)模擬在橢圓極化波等輻照條件下的電磁暴露場(chǎng)景。

    實(shí)現(xiàn)多波源入射電磁暴露數(shù)值仿真系統(tǒng)的核心在于能夠?qū)⒍鄠€(gè)波源按照設(shè)定的條件一同加入到FDTD計(jì)算空間中,本文根據(jù)FDTD方法中總場(chǎng)散射場(chǎng)邊界(TF/SF)加入波源的方法[23],將一維平面波按照不同的入射條件投影到三維計(jì)算空間的TF/SF 邊界上,在投影邊界上各電磁場(chǎng)分量見式(8)。式中:E0表示一維平面波的電場(chǎng)量,Z0表示真空中波阻抗,α、θ、φ 同為圖3 所示的波源入射的方位角。通過(guò)此方式,能夠?qū)蓚€(gè)波源一同加入到計(jì)算空間中,并分別體現(xiàn)在3個(gè)坐標(biāo)軸方向上的電場(chǎng)、磁場(chǎng)分量上,完成雙波源的加入。在此基礎(chǔ)上,通過(guò)與真空中解析解以及商業(yè)軟件存在介質(zhì)場(chǎng)景的結(jié)果進(jìn)行對(duì)比和分析,得到了與正確結(jié)果完全符合的數(shù)據(jù),由此說(shuō)明本文所建立的多波源電磁仿真系統(tǒng)在處理電磁輻射問題上是可靠的。

    2 結(jié)果與分析

    2.1 入射波垂直大鼠體長(zhǎng)方向

    首先考慮入射方向垂直大鼠體長(zhǎng)方向的情況。

    2.1.1 側(cè)面入射

    由+x 方向,即波從大鼠模型的左側(cè)面向右側(cè)面入射的數(shù)值模擬結(jié)果如圖4所示。圖4中曲線表示W(wǎng)BASAR值隨橢圓極化波源橢偏率的變化,上下兩條虛線分別表示EzHyKx、EyHzKx兩線極化波源的WBASAR結(jié)果。

    從圖4可以看出,在側(cè)面入射條件下,橢圓極化波輻照大鼠的WBASAR值隨著相位差的改變而變化,即表明橢圓極化波源的橢偏率是影響電磁暴露結(jié)果的重要因素之一。同時(shí),在橢圓極化波橢偏率的一個(gè)變化周期中,WBASAR 也呈現(xiàn)出一定的對(duì)稱性分布。

    觀察圖4中橢圓極化波的WBASAR變化曲線,以?φ 取π/6 和11π/6 為例進(jìn)行分析,如圖2(b)、(m)所示的橢圓極化波形態(tài)。這兩種情況所對(duì)應(yīng)波源的橢偏率相同,但橢圓極化波的旋向相反,相位差為π/6時(shí)對(duì)應(yīng)的是右旋橢圓極化波,11π/6對(duì)應(yīng)左旋橢圓極化波。對(duì)比這兩者的WBASAR結(jié)果可以發(fā)現(xiàn),橢偏率相同但旋向相反的橢圓極化波源得到的結(jié)果不同。這兩種情況下的波源表達(dá)見式(9)。

    兩種情況下的相位差分別為π/6和11π/6,且兩相位差之和為2π,故可以得到第二種橢圓極化波電場(chǎng)y方向分量的簡(jiǎn)化關(guān)系式見式(10)。

    比較Ey1和Ey2的表達(dá)式可以看出ωt 不斷變化時(shí),兩種情況下的橢圓極化波源的變化規(guī)律是不同的,故而產(chǎn)生了不同的電磁輻照結(jié)果,這與Shckorbatov[18]的研究發(fā)現(xiàn)是相吻合的。觀察其他橢偏率相同而旋向相反的波源所對(duì)應(yīng)的WBASAR結(jié)果,同樣有上述結(jié)論。但當(dāng)相位差為π/2 和3π/2時(shí),即得到圓極化波源的y 方向電場(chǎng)表達(dá)式如式(11)。

    從Ey1和Ey2的表達(dá)式中可以看出兩者隨時(shí)間的變化規(guī)律恰好等大反向,從而可知兩旋向相反橢圓極化波的電場(chǎng)量變化規(guī)律同樣呈現(xiàn)出大小相等、方向相反的規(guī)律,故兩種旋向相反的圓極化波源會(huì)得到相同的WBASAR 結(jié)果, 這同樣與Shckorbatov 關(guān)于不同旋向的圓極化波源對(duì)人類細(xì)胞影響的研究發(fā)現(xiàn)是相吻合的[19]。故可得出結(jié)論:當(dāng)橢圓極化波源以側(cè)面輻照大鼠時(shí),波源的旋向也是影響電磁暴露結(jié)果WBASAR大小的重要因素之一。

    此外,比較橢圓極化波與EzHyKx和EyHzKx兩種線極化波的WBASAR 結(jié)果值,可以發(fā)現(xiàn)前者的WBASAR 變化范圍始終被包含在兩線極化波的結(jié)果中。即可以得出結(jié)論:當(dāng)從側(cè)面方向輻照大鼠時(shí),與入射波能量相同條件下的橢圓極化波WBASAR 結(jié)果相比,線極化波EzHyKx的結(jié)果總大于不同橢偏率的橢圓極化波輻照結(jié)果,而線極化波EyHzKx的結(jié)果則總小于橢圓極化波。波源的電場(chǎng)振動(dòng)方向與大鼠體態(tài)的長(zhǎng)軸和短軸的關(guān)系是產(chǎn)生這一現(xiàn)象的重要原因,如圖5所示。

    線極化波EzHyKx的電場(chǎng)振動(dòng)方向與大鼠的長(zhǎng)軸方向一致,EyHzKx的電場(chǎng)振動(dòng)方向沿大鼠的體寬短軸方向,而橢圓極化波的電場(chǎng)振動(dòng)方向則是隨著時(shí)間空間不斷地旋轉(zhuǎn),即在大鼠的長(zhǎng)軸和短軸之間呈周期性分布,從而使得同樣條件下,橢圓極化波經(jīng)過(guò)大鼠組織時(shí)產(chǎn)生的能量吸收也處于兩線極化波之間,故表現(xiàn)為WBASAR的結(jié)果出現(xiàn)上述分布特點(diǎn)。

    2.1.2 背部入射

    為進(jìn)一步考察在入射方向改變時(shí)是否仍然能夠得到上述結(jié)果,改變電磁波為+y 方向入射,即從大鼠背部入射進(jìn)行數(shù)值模擬計(jì)算,得到WBASAR結(jié)果如圖6所示。

    整體觀察背部入射時(shí)的橢圓極化波電磁暴露數(shù)值模擬結(jié)果,可以發(fā)現(xiàn)WBASAR值隨著波源橢偏率的改變也呈現(xiàn)出對(duì)稱性的周期變化趨勢(shì),且其對(duì)稱性明顯高于從側(cè)面入射時(shí)橢圓極化波的WBASAR 變化曲線,這是由于大鼠的組織分布沿著背部入射時(shí)的截面較側(cè)面入射時(shí)呈現(xiàn)出較高的對(duì)稱性,如圖7所示。橢圓極化波橢偏率變化的高度對(duì)稱性和大鼠組織分布的對(duì)稱性共同作用,使得背部入射時(shí)的WBASAR結(jié)果呈現(xiàn)出高度的對(duì)稱性周期分布。

    為觀察在背部入射時(shí)旋向?qū)BASAR 結(jié)果的影響,分別對(duì)比相位差為π/6和11π/6、2π/6 和10π/6、4π/6 和8π/6、5π/6 和7π/6 等4 種橢偏率相同但旋向相反的橢圓極化波源,其WBASAR結(jié)果的差異平均在11.2%左右。這說(shuō)明在背部入射時(shí),橢圓極化波源的旋向依然是影響電磁暴露WBASAR結(jié)果的重要因素。對(duì)于圓極化波源,左旋和右旋的WBASAR結(jié)果分別為0.032 81 mW/kg和0.032 89 mW/kg,較側(cè)面入射時(shí)更為接近,即仍然可以認(rèn)為旋向?qū)A極化電磁波的電磁暴露結(jié)果幾乎不產(chǎn)生影響。

    此外,將橢圓極化波在背部入射所得WBASAR 結(jié)果,與此入射方向上兩線極化波EzHxKy和ExHzKy結(jié)果進(jìn)行比較,可以得到與側(cè)面入射時(shí)類似的結(jié)論,即電場(chǎng)極化方向平行于大鼠長(zhǎng)軸的WBASASR 結(jié)果總大于相同波源能量條件下的橢圓極化波結(jié)果,而電場(chǎng)極化方向平行于大鼠體寬短軸的線極化波的結(jié)果總小于橢圓極化波。

    2.1.3 其他入射波垂直大鼠體長(zhǎng)的方向

    取同樣垂直大鼠長(zhǎng)軸的腹部(-y)和側(cè)背部(xOy)方向進(jìn)行電磁輻照,側(cè)背部入射時(shí)的波源波矢方向沿計(jì)算空間的x 軸與y 軸角平分線方向,以xOy 表示,所得結(jié)果如圖8 所示。從圖8 可以看到,所得WBASAR 的變化規(guī)律與上述基本相同,進(jìn)一步論證了波源沿垂直大鼠長(zhǎng)軸方向進(jìn)行輻照時(shí)所得結(jié)論。

    2.2 入射波平行大鼠體長(zhǎng)方向

    與上述入射方向相比,平行于大鼠體長(zhǎng)方向是另一類較為特殊的電磁波入射方向,根據(jù)大鼠沿頭部方向的體態(tài)分布為例進(jìn)行分析,結(jié)果如圖9所示。大鼠在此方向上僅僅存在體寬短軸(x軸方向)和體高短軸(y 軸方向),即不存在明確的體長(zhǎng)軸與短軸的分布特點(diǎn),且大鼠在體寬短軸方向上的組織分布要略大于體高短軸。同時(shí)在此方向上大鼠的身體截面多呈近似圓形或近似橢圓的形狀,即表現(xiàn)出較為微小的組織分布變化,這與上述入射方向上的大鼠體態(tài)分布特點(diǎn)相差較大。

    2.2.1 頭部入射

    如圖10 所示,整體來(lái)看頭部入射的結(jié)果,無(wú)論是橢圓極化波源橢偏率的改變,還是線極化波源極化方向的改變都對(duì)WBASAR 值的影響較小,整體變化范圍僅在7.7%以內(nèi),這與大鼠沿體長(zhǎng)軸方向截面的組織分布變化微小的特點(diǎn)是相符合的。同時(shí),橢圓極化波的WBASAR結(jié)果表現(xiàn)出更強(qiáng)的對(duì)稱性分布特點(diǎn),且觀察橢偏率相同的WBASAR值幾乎相等,即表明頭部入射時(shí)旋向?qū)BASAR值的影響已經(jīng)可以忽略。橢圓極化波在頭部入射時(shí)表現(xiàn)出與垂直大鼠體長(zhǎng)方向入射時(shí)的變化規(guī)律明顯不同,這與此兩種情況下大鼠組織分布特點(diǎn)的顯著差異是相符合的。

    對(duì)比線極化波和橢圓極化波的WBASAR 結(jié)果可以發(fā)現(xiàn),有部分形態(tài)橢圓極化波的結(jié)果值比大鼠體高方向短軸極化的線極化波(EyHxKz)結(jié)果更小,即橢圓極化波的WBASAR 變化范圍超出了ExHyKz和EyHxKz對(duì)應(yīng)的WBASAR 值之外,意味著在頭部入射的情況下,線極化波的WBASAR值已經(jīng)不能夠代表橢圓極化波結(jié)果的變化范圍。

    2.2.2 尾部入射

    波源沿-z(尾部)方向入射同樣平行于大鼠體長(zhǎng)方向,其電磁暴露數(shù)值模擬的WBASAR結(jié)果如圖11 所示。從圖11 可以看出,尾部入射的WBASAR 結(jié)果整體變化范圍仍很小,且同樣表現(xiàn)出較高的對(duì)稱性分布,同時(shí)旋向產(chǎn)生的影響可以忽略。此外,橢圓極化波的變化范圍同樣超出ExHyKz和EyHxKz兩線極化波范圍之外。由此表明,尾部入射時(shí)橢圓極化波的WBASAR變化規(guī)律與頭部入射時(shí)幾乎一致,即論證了波源沿平行大鼠體長(zhǎng)方向入射時(shí)所得結(jié)論的可靠性。

    2.3 重要組織的平均SAR

    對(duì)于電磁暴露劑量的評(píng)估,大鼠重要組織的平均SAR 同樣是非常重要的參考依據(jù),在本文的數(shù)值仿真中,計(jì)算大鼠的大腦、心臟、肝臟、睪丸等重要組織的平均SAR,不同入射方向上的結(jié)果如圖12所示。圖12中的曲線代表不同組織的平均SAR以及此入射方向上WBASAR隨橢偏率的變化曲線,每種組織平均SAR 曲線對(duì)應(yīng)顏色相同的兩條虛線,他們分別代表與圖4、6、10 和11 中描述的線極化波相同條件下所計(jì)算得到的對(duì)應(yīng)組織的平均SAR,圖12中(e)為相位差?φ取π/6橢圓極化波源側(cè)面入射時(shí),沿坐標(biāo)軸方向某一截面上的SAR分布圖。

    2.3.1 橢圓極化波條件下組織平均SAR的特點(diǎn)

    從圖12 可以看出,大鼠重要組織的平均SAR隨橢圓極化波橢偏率的變化呈現(xiàn)與WBASAR相類似的變化規(guī)律,但不同的組織在不同入射條件下的結(jié)果又有所不同,大腦和睪丸等位于體表組織的平均SAR總大于此條件下的WBASAR,這也是大腦、睪丸等組織作為電磁暴露研究重要對(duì)象的一個(gè)原因[31-35];而肝臟等位于大鼠內(nèi)部靠中心位置,其組織平均SAR 相比其他組織,整體變化范圍很小,且始終小于WBASAR;對(duì)于心臟等偏向于大鼠身體一側(cè)的組織,其平均SAR與WBASAR的大小關(guān)系隨著入射方向的改變產(chǎn)生較大的變化,當(dāng)波源沿垂直大鼠長(zhǎng)軸方向入射時(shí),心臟的組織平均SAR幾乎完全小于WBASAR,而當(dāng)波源沿平行大鼠長(zhǎng)軸方向入射時(shí),心臟的組織平均SAR 與WBASAR 的大小關(guān)系隨橢偏率的改變產(chǎn)生較大差異。

    此外,比較圖12(c)和(d)中兩種輻照條件下大腦和睪丸的組織平均SAR 可以發(fā)現(xiàn),頭部入射時(shí)橢圓極化波的大腦平均SAR 隨著橢偏率的變化有較大范圍的波動(dòng),而睪丸組織平均SAR 的變化則相對(duì)平穩(wěn),然而當(dāng)波源沿尾部入射時(shí),兩者表現(xiàn)出相反的變化規(guī)律。這說(shuō)明波源與組織之間的相對(duì)位置關(guān)系對(duì)橢圓極化波條件下的組織平均SAR 也有較大的影響,即靠近波源的組織的平均SAR 會(huì)隨橢圓極化波橢偏率的改變產(chǎn)生較大的變化,而距離波源相對(duì)較遠(yuǎn)的組織,其平均SAR隨橢偏率的變化相對(duì)較小。

    2.3.2 橢圓極化波與線極化波對(duì)應(yīng)組織平均SAR之間的差異

    與WBASAR 類似,橢圓極化波源與線極化波源條件下的組織平均SAR 在分布范圍上亦有所不同。當(dāng)波源從側(cè)面或背部等垂直大鼠長(zhǎng)軸的方向入射時(shí),橢圓極化波的WBASAR整體分布在同方向入射的兩線極化波范圍之內(nèi),但是觀察波源從側(cè)面入射時(shí)的結(jié)果如圖12(a)所示,大腦區(qū)域的平均SAR明顯分布在EzHyKx和EyHzKx兩線極化波的大腦平均SAR 范圍之外,分別超出EzHyKx的19.9%,并低于EyHzKx的34.2%;然而,當(dāng)波源從背部入射時(shí),橢圓極化波條件下大腦組織的平均SAR的大小完全分布在線極化波SAR的范圍之內(nèi)。同樣,與線極化波條件下的組織平均SAR 相比,盡管睪丸組織的平均SAR大于WBASAR,但是無(wú)論波源從側(cè)面入射還是背部入射,其平均SAR 的大小都分布在兩線極化波條件下睪丸組織平均SAR的范圍之內(nèi)。

    當(dāng)波源平行大鼠長(zhǎng)軸方向入射時(shí),如圖12(c)、(d)所示,兩種波源入射條件下的心臟等組織的平均SAR在WBASAR的上下不斷浮動(dòng),不同之處在于波源從頭部入射時(shí),橢圓極化波條件下的心臟平均SAR分布范圍大于ExHyKz和EyHxKz對(duì)應(yīng)的線極化波結(jié)果,即部分橢偏率對(duì)應(yīng)橢圓極化波條件下的心臟平均SAR高出EyHxKz的10.7%,并低于ExHyKz的22.2%;而尾部入射時(shí)橢圓極化波的心臟平均SAR 的分布被完全涵蓋在線極化波中。對(duì)于肝臟組織,橢圓極化波沿頭部入射和尾部入射時(shí)的肝臟組織平均SAR 主要分布在線極化波肝臟組織平均SAR 的范圍之內(nèi),僅有部分超出了兩線極化波的范圍;相比之下,波源沿垂直大鼠長(zhǎng)軸方向入射時(shí),不同橢偏率條件下的肝臟組織平均SAR 幾乎都在兩線極化波范圍之外,且大部分結(jié)果比線極化波EyHzKx/ExHzKy對(duì)應(yīng)的肝臟平均SAR更小,這表明波源入射方向雖然對(duì)橢圓極化波輻照下肝臟等組織的平均SAR 整體變化范圍的影響較小,但對(duì)于橢圓極化波與線極化波結(jié)果的大小關(guān)系會(huì)產(chǎn)生較大的影響。

    3 結(jié)論

    在橢圓極化波的電磁暴露數(shù)值模擬計(jì)算中,波源的橢偏率是影響電磁輻照結(jié)果的重要因素,其WBASAR結(jié)果會(huì)隨著波源橢偏率的對(duì)稱性周期性變化,表現(xiàn)出一定的對(duì)稱性周期分布,且此對(duì)稱性隨波源入射方向上大鼠的組織分布對(duì)稱性的提高而升高。與線極化波類似,橢圓極化波源的入射方向同樣是影響動(dòng)物電磁暴露劑量評(píng)估的重要因素,表現(xiàn)為波源入射方向與大鼠體長(zhǎng)方向相互垂直和平行時(shí)分別表現(xiàn)出明顯不同的變化規(guī)律。

    橢圓極化波的旋向同樣會(huì)對(duì)電磁暴露的結(jié)果產(chǎn)生重要影響,當(dāng)波源垂直大鼠體長(zhǎng)方向入射時(shí),旋向?qū)E圓極化波的WBASAR值的影響平均達(dá)到11.0%以上,但波源以大鼠體長(zhǎng)方向入射時(shí),旋向產(chǎn)生的影響幾乎可以忽略。而對(duì)于橢圓極化波中的特殊情況——圓極化電磁波來(lái)說(shuō),無(wú)論入射方向如何變化,旋向幾乎不對(duì)電磁暴露結(jié)果產(chǎn)生影響。

    對(duì)比各個(gè)入射方向上相同能量條件下的線極化波和橢圓極化波電磁暴露結(jié)果可以得出結(jié)論:線極化電磁波作為波源對(duì)動(dòng)物電磁暴露劑量大小的評(píng)估不能夠完全代表非極化波電磁波,即不能夠完全涵蓋橢圓極化波的變化范圍。表現(xiàn)為當(dāng)波源入射方向垂直大鼠體長(zhǎng)方向時(shí),橢圓極化波的WBASAR 變化范圍總處在線極化波結(jié)果之間;但入射方向與大鼠體長(zhǎng)方向平行時(shí),線極化波的結(jié)果僅能包含部分形態(tài)的橢圓極化波,即此時(shí)不能夠代表橢圓極化波源對(duì)動(dòng)物電磁暴露劑量大小的評(píng)估。

    在兩種波源條件下對(duì)組織內(nèi)部電磁暴露劑量的評(píng)估同樣存在差異,橢圓極化波WBASAR的分布是否處于同能量條件下相同入射方向上“水平”或“垂直”極化的線極化波WBASAR的極值范圍內(nèi),并不能決定各組織的平均SAR 也處在此時(shí)線極化波條件下對(duì)應(yīng)組織的平均SAR 范圍之內(nèi),且橢圓極化波條件下的組織平均SAR 還受組織在大鼠體內(nèi)的分布位置、與體表的距離,以及與波源之間的相對(duì)位置關(guān)系等因素的影響,并在不同入射方向上呈現(xiàn)出不同的分布規(guī)律。

    所以,對(duì)動(dòng)物電磁暴露劑量的評(píng)估不僅要參考線極化波的結(jié)果,同時(shí)也需要綜合考慮非線極化波源的電磁暴露結(jié)果,從而得到更加全面的結(jié)論。

    猜你喜歡
    線極化波源電磁
    雙波源同向交替入射下大鼠的電磁暴露仿真
    基于回折線極化器技術(shù)的VICTS天線全極化可調(diào)實(shí)現(xiàn)
    微波無(wú)線功率傳輸中高效率雙線極化整流天線設(shè)計(jì)①
    平面電磁波極化場(chǎng)矢量的時(shí)空機(jī)理
    三維多孔電磁復(fù)合支架構(gòu)建與理化表征
    掌握基礎(chǔ)知識(shí) 不懼電磁偏轉(zhuǎn)
    衛(wèi)星通信地球站天線的線極化角計(jì)算研究*
    低阻抗沖擊電壓標(biāo)準(zhǔn)波源
    電磁換向閥應(yīng)用探討
    河南科技(2014年16期)2014-02-27 14:13:21
    瞬變電磁法在煤礦防治水中的應(yīng)用
    河南科技(2014年8期)2014-02-27 14:07:41
    中文在线观看免费www的网站| 中文字幕熟女人妻在线| 我的女老师完整版在线观看| 亚洲成人中文字幕在线播放| 亚洲国产欧美人成| 成人性生交大片免费视频hd| 国产综合懂色| 久久久久国产网址| a级一级毛片免费在线观看| 天堂av国产一区二区熟女人妻| 亚洲欧美成人精品一区二区| 白带黄色成豆腐渣| 麻豆成人av视频| 久久人妻av系列| 精品99又大又爽又粗少妇毛片| 亚洲精品456在线播放app| 成人鲁丝片一二三区免费| 国产69精品久久久久777片| 又爽又黄无遮挡网站| 免费在线观看成人毛片| 少妇熟女欧美另类| 国产日韩欧美在线精品| 国产精品无大码| 日日啪夜夜撸| 汤姆久久久久久久影院中文字幕 | 99国产精品一区二区蜜桃av| 免费大片18禁| 噜噜噜噜噜久久久久久91| 一个人观看的视频www高清免费观看| 高清av免费在线| 国产一区二区在线av高清观看| 高清av免费在线| 18禁在线播放成人免费| 国产高清有码在线观看视频| 男人和女人高潮做爰伦理| av线在线观看网站| 白带黄色成豆腐渣| 亚洲欧美成人综合另类久久久 | 免费看光身美女| 精品人妻偷拍中文字幕| 国产精品三级大全| 中文乱码字字幕精品一区二区三区 | 久久午夜福利片| 韩国高清视频一区二区三区| 中国美白少妇内射xxxbb| 国产爱豆传媒在线观看| 国产一区亚洲一区在线观看| 99久久精品国产国产毛片| 欧美成人午夜免费资源| 人人妻人人澡人人爽人人夜夜 | 高清日韩中文字幕在线| 欧美成人a在线观看| 欧美成人免费av一区二区三区| 免费av观看视频| 国产一级毛片七仙女欲春2| 国产精品av视频在线免费观看| 观看免费一级毛片| 欧美日本视频| 亚洲欧美日韩无卡精品| 久久精品久久精品一区二区三区| 久久99热6这里只有精品| 天堂av国产一区二区熟女人妻| 看免费成人av毛片| 日韩国内少妇激情av| 成人综合一区亚洲| 国产精品国产高清国产av| 三级国产精品欧美在线观看| 国产亚洲5aaaaa淫片| 日韩欧美精品v在线| 免费播放大片免费观看视频在线观看 | 国产乱人视频| 中文天堂在线官网| 少妇人妻一区二区三区视频| 日韩大片免费观看网站 | 亚洲一区高清亚洲精品| 欧美高清性xxxxhd video| 天堂av国产一区二区熟女人妻| 白带黄色成豆腐渣| 亚洲中文字幕日韩| 男的添女的下面高潮视频| 26uuu在线亚洲综合色| 亚洲欧洲日产国产| 亚洲精品影视一区二区三区av| 国产真实乱freesex| 干丝袜人妻中文字幕| 最近最新中文字幕大全电影3| 国产午夜精品一二区理论片| 国产精品国产高清国产av| av福利片在线观看| 国产亚洲5aaaaa淫片| 神马国产精品三级电影在线观看| 国产午夜精品久久久久久一区二区三区| 国产精品,欧美在线| 亚洲欧洲国产日韩| 国产高清视频在线观看网站| 激情 狠狠 欧美| 毛片女人毛片| 免费看av在线观看网站| 久久精品熟女亚洲av麻豆精品 | 亚洲成av人片在线播放无| 国产久久久一区二区三区| 久热久热在线精品观看| 精品99又大又爽又粗少妇毛片| videos熟女内射| 少妇的逼水好多| 国产又色又爽无遮挡免| 成人高潮视频无遮挡免费网站| 不卡视频在线观看欧美| 日本熟妇午夜| 你懂的网址亚洲精品在线观看 | 国产片特级美女逼逼视频| 国产精品伦人一区二区| 国产亚洲av片在线观看秒播厂 | 免费看a级黄色片| 日本免费一区二区三区高清不卡| 久久亚洲精品不卡| 久久精品综合一区二区三区| 欧美激情在线99| av在线播放精品| 一级毛片我不卡| 日本爱情动作片www.在线观看| 欧美日本亚洲视频在线播放| av黄色大香蕉| 国产欧美日韩精品一区二区| 中文字幕av在线有码专区| 国产免费又黄又爽又色| 亚洲av中文字字幕乱码综合| 国产伦在线观看视频一区| 欧美日韩一区二区视频在线观看视频在线 | 亚洲av电影不卡..在线观看| 欧美另类亚洲清纯唯美| 亚洲av成人精品一区久久| 一个人看视频在线观看www免费| 亚洲第一区二区三区不卡| 日韩人妻高清精品专区| 五月伊人婷婷丁香| 亚洲内射少妇av| 嫩草影院入口| 建设人人有责人人尽责人人享有的 | 男女边吃奶边做爰视频| 99久久无色码亚洲精品果冻| 亚洲人与动物交配视频| 亚洲国产精品国产精品| av在线亚洲专区| 亚洲欧美精品自产自拍| 欧美日韩精品成人综合77777| 高清av免费在线| 一二三四中文在线观看免费高清| 国产黄色视频一区二区在线观看 | 亚洲成人中文字幕在线播放| 中文精品一卡2卡3卡4更新| 91久久精品国产一区二区成人| 三级毛片av免费| 国产大屁股一区二区在线视频| 亚洲丝袜综合中文字幕| 在线观看66精品国产| 村上凉子中文字幕在线| 精品午夜福利在线看| 国产精品久久久久久av不卡| 亚洲成人久久爱视频| 女人被狂操c到高潮| 精品国内亚洲2022精品成人| 网址你懂的国产日韩在线| 汤姆久久久久久久影院中文字幕 | 日韩人妻高清精品专区| 久久久久久久久大av| 免费观看性生交大片5| 国产淫片久久久久久久久| 亚洲国产成人一精品久久久| 免费观看精品视频网站| 91精品一卡2卡3卡4卡| 日本-黄色视频高清免费观看| 精品久久久久久久久亚洲| 一级二级三级毛片免费看| 黑人高潮一二区| 床上黄色一级片| 国产成人a∨麻豆精品| 91在线精品国自产拍蜜月| 成人性生交大片免费视频hd| 极品教师在线视频| 国产伦一二天堂av在线观看| 亚洲精品aⅴ在线观看| 亚洲国产精品专区欧美| 青春草亚洲视频在线观看| 国模一区二区三区四区视频| 午夜精品国产一区二区电影 | 国产伦一二天堂av在线观看| 日本黄大片高清| 亚洲av成人精品一二三区| 男女边吃奶边做爰视频| 亚洲四区av| 国产一级毛片在线| 国产午夜精品一二区理论片| 在线观看66精品国产| 国产 一区精品| 欧美一区二区国产精品久久精品| 春色校园在线视频观看| 午夜精品国产一区二区电影 | 亚洲欧洲国产日韩| 国产精品女同一区二区软件| 亚洲欧美日韩高清专用| 真实男女啪啪啪动态图| 亚洲最大成人中文| 身体一侧抽搐| 色综合亚洲欧美另类图片| 丰满乱子伦码专区| 亚洲婷婷狠狠爱综合网| 国产精品久久久久久久久免| 看片在线看免费视频| 99视频精品全部免费 在线| 成人三级黄色视频| 久久久久久大精品| 国产精品熟女久久久久浪| 99久久人妻综合| 国产成人freesex在线| 久久久午夜欧美精品| АⅤ资源中文在线天堂| 1024手机看黄色片| 联通29元200g的流量卡| 欧美日韩一区二区视频在线观看视频在线 | 高清午夜精品一区二区三区| 免费观看人在逋| 又爽又黄a免费视频| 国产精品综合久久久久久久免费| 青春草亚洲视频在线观看| 国产精品永久免费网站| 欧美日本亚洲视频在线播放| 在线免费十八禁| 看免费成人av毛片| АⅤ资源中文在线天堂| 国产精品永久免费网站| 久久国产乱子免费精品| 三级毛片av免费| 亚洲经典国产精华液单| 九色成人免费人妻av| 亚洲av电影不卡..在线观看| 亚洲真实伦在线观看| 插阴视频在线观看视频| 大香蕉久久网| 国产精品女同一区二区软件| 秋霞伦理黄片| 黑人高潮一二区| 蜜桃亚洲精品一区二区三区| ponron亚洲| 国产精品女同一区二区软件| 国产精品野战在线观看| 亚洲精品色激情综合| 成年免费大片在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产伦理片在线播放av一区| 九色成人免费人妻av| 一卡2卡三卡四卡精品乱码亚洲| 1000部很黄的大片| 免费搜索国产男女视频| 久久亚洲国产成人精品v| 国产精品,欧美在线| 女人被狂操c到高潮| 九九在线视频观看精品| 69av精品久久久久久| 亚洲av中文字字幕乱码综合| 国产一区亚洲一区在线观看| 欧美成人免费av一区二区三区| 国产精品综合久久久久久久免费| 国产精华一区二区三区| 欧美最新免费一区二区三区| 综合色丁香网| 麻豆一二三区av精品| 国内揄拍国产精品人妻在线| 免费在线观看成人毛片| 国产精品综合久久久久久久免费| 99久久精品国产国产毛片| 日韩国内少妇激情av| 亚洲av福利一区| 国产91av在线免费观看| 日韩av在线免费看完整版不卡| 天天躁日日操中文字幕| 久久精品国产亚洲av涩爱| 亚洲内射少妇av| 亚洲国产精品国产精品| 男人舔女人下体高潮全视频| 国产亚洲一区二区精品| 男的添女的下面高潮视频| 日韩av在线大香蕉| 高清视频免费观看一区二区 | 亚洲精品国产av成人精品| 内射极品少妇av片p| 亚洲国产欧美人成| 青春草视频在线免费观看| 少妇丰满av| 国产人妻一区二区三区在| 国产精品一区二区在线观看99 | 日日干狠狠操夜夜爽| 欧美变态另类bdsm刘玥| 精品午夜福利在线看| 九九爱精品视频在线观看| 亚洲激情五月婷婷啪啪| 精品欧美国产一区二区三| 少妇高潮的动态图| 91在线精品国自产拍蜜月| 秋霞伦理黄片| 国产av不卡久久| 国产91av在线免费观看| 婷婷色综合大香蕉| 免费看日本二区| 日韩精品青青久久久久久| 毛片女人毛片| 夜夜看夜夜爽夜夜摸| 91久久精品国产一区二区成人| 亚洲精华国产精华液的使用体验| 国产国拍精品亚洲av在线观看| 精品久久久久久成人av| 日本色播在线视频| 成人鲁丝片一二三区免费| 成人毛片a级毛片在线播放| 久久韩国三级中文字幕| 高清视频免费观看一区二区 | 能在线免费观看的黄片| 天堂中文最新版在线下载 | 亚洲成人中文字幕在线播放| 亚洲一区高清亚洲精品| 只有这里有精品99| 三级毛片av免费| 三级国产精品欧美在线观看| 日本欧美国产在线视频| 国产高清不卡午夜福利| 亚洲婷婷狠狠爱综合网| 亚洲成色77777| 久久99热这里只频精品6学生 | 亚州av有码| 中文字幕免费在线视频6| 久久99蜜桃精品久久| 国产黄色视频一区二区在线观看 | 国产视频内射| 国产极品天堂在线| 51国产日韩欧美| 国产精品一区二区性色av| 欧美成人a在线观看| 一个人看视频在线观看www免费| 综合色av麻豆| kizo精华| 国产一区二区亚洲精品在线观看| 可以在线观看毛片的网站| 久久久久久久久大av| 精品久久久久久成人av| 日韩一本色道免费dvd| 国产黄色小视频在线观看| 国产乱人偷精品视频| 亚洲无线观看免费| 精华霜和精华液先用哪个| 两个人视频免费观看高清| 免费黄网站久久成人精品| 国产高清不卡午夜福利| 国产欧美另类精品又又久久亚洲欧美| 国产真实乱freesex| 麻豆成人av视频| 精品国产三级普通话版| 日本一二三区视频观看| av专区在线播放| 婷婷色综合大香蕉| 2021少妇久久久久久久久久久| 天堂av国产一区二区熟女人妻| 欧美bdsm另类| 一本一本综合久久| 国产伦一二天堂av在线观看| 日韩大片免费观看网站 | 欧美日韩国产亚洲二区| 免费av毛片视频| 一级黄片播放器| 麻豆久久精品国产亚洲av| 一区二区三区乱码不卡18| 成年版毛片免费区| 久久午夜福利片| 一个人看的www免费观看视频| 丝袜美腿在线中文| 网址你懂的国产日韩在线| 一级二级三级毛片免费看| 国产午夜精品一二区理论片| 我要看日韩黄色一级片| 亚洲精品乱码久久久v下载方式| 纵有疾风起免费观看全集完整版 | av播播在线观看一区| 午夜福利在线观看吧| 赤兔流量卡办理| 男女下面进入的视频免费午夜| 毛片女人毛片| 欧美激情国产日韩精品一区| 美女黄网站色视频| 99热这里只有是精品在线观看| 人体艺术视频欧美日本| 青春草亚洲视频在线观看| videos熟女内射| 91aial.com中文字幕在线观看| 一区二区三区乱码不卡18| 国产伦在线观看视频一区| 国产精品人妻久久久影院| 久久99热6这里只有精品| 亚洲精品亚洲一区二区| 国产麻豆成人av免费视频| 亚洲精品aⅴ在线观看| 久久久色成人| 国产免费一级a男人的天堂| 亚洲最大成人手机在线| 久久久久久九九精品二区国产| 好男人在线观看高清免费视频| 国产免费福利视频在线观看| 1024手机看黄色片| 亚洲精品亚洲一区二区| 欧美3d第一页| 久久久久精品久久久久真实原创| 蜜桃久久精品国产亚洲av| 色噜噜av男人的天堂激情| 亚洲国产成人一精品久久久| 久久久成人免费电影| 中文字幕人妻熟人妻熟丝袜美| 国产成人免费观看mmmm| 中文在线观看免费www的网站| 国产精品国产三级国产av玫瑰| 老师上课跳d突然被开到最大视频| 99久久中文字幕三级久久日本| 真实男女啪啪啪动态图| 日本猛色少妇xxxxx猛交久久| 久久这里有精品视频免费| 男女国产视频网站| 内射极品少妇av片p| 久久6这里有精品| 岛国毛片在线播放| 美女内射精品一级片tv| 久久精品91蜜桃| 中文乱码字字幕精品一区二区三区 | 国产老妇女一区| 午夜老司机福利剧场| 国产精品久久久久久久电影| 国产人妻一区二区三区在| 亚洲精品日韩在线中文字幕| 免费大片18禁| 免费看av在线观看网站| 国产一区二区在线观看日韩| 天堂影院成人在线观看| 亚洲一级一片aⅴ在线观看| 国产成人精品一,二区| 国产成人福利小说| 18+在线观看网站| 亚洲av一区综合| 少妇高潮的动态图| 99久久精品国产国产毛片| 亚洲精品aⅴ在线观看| 最近的中文字幕免费完整| 久久99精品国语久久久| 99久国产av精品| 九九久久精品国产亚洲av麻豆| 久久久a久久爽久久v久久| 九色成人免费人妻av| kizo精华| 91久久精品国产一区二区成人| 国产高清不卡午夜福利| 麻豆成人午夜福利视频| 欧美日本亚洲视频在线播放| 亚洲美女搞黄在线观看| 国产精品一及| 国产淫语在线视频| 久久久精品94久久精品| 麻豆成人av视频| 99久久无色码亚洲精品果冻| 色网站视频免费| 欧美色视频一区免费| 精品人妻偷拍中文字幕| 久久草成人影院| www.色视频.com| 国产一区有黄有色的免费视频 | 免费看美女性在线毛片视频| 免费无遮挡裸体视频| 亚洲综合精品二区| 男的添女的下面高潮视频| 国产精品一二三区在线看| 美女大奶头视频| 色尼玛亚洲综合影院| 又粗又硬又长又爽又黄的视频| 18+在线观看网站| 两个人视频免费观看高清| 九九在线视频观看精品| 最近的中文字幕免费完整| 哪个播放器可以免费观看大片| 国产精品永久免费网站| 有码 亚洲区| 久久久久久久久久成人| 性色avwww在线观看| 人妻夜夜爽99麻豆av| 青春草国产在线视频| 男人舔女人下体高潮全视频| 国产成人a区在线观看| 中文字幕熟女人妻在线| 一个人看的www免费观看视频| 一本一本综合久久| 黄片无遮挡物在线观看| 日韩欧美国产在线观看| 久久久久久大精品| 国产白丝娇喘喷水9色精品| 国产精华一区二区三区| 日本五十路高清| 亚洲精品乱码久久久v下载方式| 国产高清视频在线观看网站| 亚洲性久久影院| 免费搜索国产男女视频| 亚洲欧美精品自产自拍| 韩国高清视频一区二区三区| 久久亚洲精品不卡| 成人欧美大片| 女人十人毛片免费观看3o分钟| 免费观看在线日韩| 久久精品国产亚洲网站| 精品久久久久久电影网 | 免费看美女性在线毛片视频| 欧美日韩综合久久久久久| 99在线视频只有这里精品首页| 成人毛片60女人毛片免费| 久久国内精品自在自线图片| av.在线天堂| 欧美成人一区二区免费高清观看| 免费无遮挡裸体视频| 卡戴珊不雅视频在线播放| 欧美日韩国产亚洲二区| 我要搜黄色片| 欧美日韩国产亚洲二区| 久久久久久大精品| 国产成人精品一,二区| 日韩欧美三级三区| 欧美高清性xxxxhd video| 欧美日韩精品成人综合77777| 免费看av在线观看网站| 日本与韩国留学比较| 亚洲欧美中文字幕日韩二区| 插逼视频在线观看| 午夜福利在线观看吧| 国语自产精品视频在线第100页| 嫩草影院新地址| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品精品国产色婷婷| 少妇的逼水好多| 亚洲精品国产成人久久av| 天堂av国产一区二区熟女人妻| 亚洲国产欧美人成| 国产午夜精品一二区理论片| 久久久久久久午夜电影| 精品久久久久久久人妻蜜臀av| 欧美性猛交黑人性爽| 亚洲欧美成人精品一区二区| 舔av片在线| 中国美白少妇内射xxxbb| 日本免费a在线| 伦精品一区二区三区| 国产高清有码在线观看视频| 亚洲欧美日韩东京热| 最近中文字幕高清免费大全6| 亚洲国产精品合色在线| 一个人观看的视频www高清免费观看| 麻豆乱淫一区二区| 亚洲四区av| 国产精品综合久久久久久久免费| 国产亚洲一区二区精品| 91av网一区二区| 欧美成人午夜免费资源| av免费在线看不卡| 亚洲在线观看片| 九九热线精品视视频播放| 国产伦一二天堂av在线观看| 国产精品久久视频播放| 日韩欧美 国产精品| 桃色一区二区三区在线观看| 久久亚洲精品不卡| 午夜精品一区二区三区免费看| 日韩强制内射视频| 少妇熟女aⅴ在线视频| 插阴视频在线观看视频| 大话2 男鬼变身卡| 乱系列少妇在线播放| 1000部很黄的大片| 久久精品国产亚洲网站| 少妇人妻一区二区三区视频| 欧美另类亚洲清纯唯美| 国产一区二区三区av在线| 草草在线视频免费看| 国产精品久久久久久精品电影| 99热精品在线国产| 日本色播在线视频| 色尼玛亚洲综合影院| 国产亚洲91精品色在线| 久久精品熟女亚洲av麻豆精品 | 久久精品国产亚洲av涩爱| 色播亚洲综合网| 成年女人看的毛片在线观看| 最近2019中文字幕mv第一页| 国产三级中文精品| 亚洲国产成人一精品久久久| 国产精品久久久久久精品电影| 大又大粗又爽又黄少妇毛片口| 欧美性猛交黑人性爽| 亚洲欧美成人精品一区二区| 国产成人a∨麻豆精品| 日韩欧美 国产精品| 亚洲精品456在线播放app| 丰满人妻一区二区三区视频av| 大又大粗又爽又黄少妇毛片口| 欧美性感艳星| 三级经典国产精品| 亚洲精品国产av成人精品| 永久免费av网站大全| 精品一区二区三区视频在线| 亚洲乱码一区二区免费版| 欧美成人午夜免费资源| 久久人人爽人人爽人人片va| 午夜免费男女啪啪视频观看| 最近中文字幕高清免费大全6| 老司机影院毛片| 国产成人aa在线观看| 看非洲黑人一级黄片|