• <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
    一级毛片电影观看| 国产精品一二三区在线看| 久久国产精品男人的天堂亚洲 | 一级毛片电影观看| 秋霞在线观看毛片| 80岁老熟妇乱子伦牲交| 久久久午夜欧美精品| 免费久久久久久久精品成人欧美视频 | 亚洲精品第二区| 国产大屁股一区二区在线视频| 黄色配什么色好看| 亚洲美女视频黄频| 一二三四中文在线观看免费高清| 亚洲美女搞黄在线观看| 免费大片黄手机在线观看| 免费高清在线观看视频在线观看| 大片电影免费在线观看免费| 黄色视频在线播放观看不卡| 日韩一区二区视频免费看| 国产成人午夜福利电影在线观看| 99精国产麻豆久久婷婷| 亚洲av男天堂| 热99国产精品久久久久久7| 乱码一卡2卡4卡精品| 黄色日韩在线| 99久久中文字幕三级久久日本| 99久久综合免费| 在线亚洲精品国产二区图片欧美 | 亚洲av中文字字幕乱码综合| 亚洲精品日韩在线中文字幕| 欧美少妇被猛烈插入视频| 老司机影院成人| 亚洲人成网站高清观看| 欧美精品一区二区大全| 亚洲精品国产成人久久av| 丝袜脚勾引网站| 精品一区二区三卡| 中文字幕av成人在线电影| 久久久国产一区二区| 99久国产av精品国产电影| 国产色爽女视频免费观看| 色网站视频免费| 一区二区av电影网| 一区二区三区免费毛片| 免费少妇av软件| 天堂中文最新版在线下载| 久久久久久久久大av| 精品久久久久久久末码| 久久精品国产a三级三级三级| 国产成人精品福利久久| 七月丁香在线播放| 99久国产av精品国产电影| 成人美女网站在线观看视频| 人妻夜夜爽99麻豆av| 各种免费的搞黄视频| 春色校园在线视频观看| 日韩欧美一区视频在线观看 | 国产视频首页在线观看| 欧美丝袜亚洲另类| 亚洲av在线观看美女高潮| 97精品久久久久久久久久精品| 久久久久精品久久久久真实原创| 久久久久久久久大av| 精华霜和精华液先用哪个| 爱豆传媒免费全集在线观看| 久久久久人妻精品一区果冻| 久久综合国产亚洲精品| 在线观看三级黄色| 亚洲国产精品一区三区| 一区二区三区四区激情视频| 国产淫片久久久久久久久| 中文字幕制服av| 尾随美女入室| 看十八女毛片水多多多| 交换朋友夫妻互换小说| 哪个播放器可以免费观看大片| 国产高清国产精品国产三级 | 免费av不卡在线播放| 亚洲天堂av无毛| 一区二区av电影网| 国产av国产精品国产| 国产淫片久久久久久久久| 国产日韩欧美在线精品| www.av在线官网国产| 不卡视频在线观看欧美| 国产成人a∨麻豆精品| 久久人人爽人人爽人人片va| 色视频在线一区二区三区| 午夜福利视频精品| 国产在视频线精品| 晚上一个人看的免费电影| 秋霞伦理黄片| 久久久久久久久久人人人人人人| 永久网站在线| 搡女人真爽免费视频火全软件| 高清av免费在线| 久久99精品国语久久久| 精华霜和精华液先用哪个| 一级毛片 在线播放| 亚洲欧洲日产国产| 好男人视频免费观看在线| 久久人人爽人人片av| 国产无遮挡羞羞视频在线观看| 国产免费又黄又爽又色| 日韩成人伦理影院| 久久久久久久久久人人人人人人| 国产精品成人在线| av免费在线看不卡| 成人漫画全彩无遮挡| 大片免费播放器 马上看| 成人二区视频| 秋霞伦理黄片| 1000部很黄的大片| 成人亚洲精品一区在线观看 | 精品人妻视频免费看| 一级毛片 在线播放| 日日摸夜夜添夜夜爱| 下体分泌物呈黄色| 欧美极品一区二区三区四区| 亚洲三级黄色毛片| 91精品一卡2卡3卡4卡| 国产免费福利视频在线观看| 全区人妻精品视频| 在线观看一区二区三区| 欧美一区二区亚洲| 国产老妇伦熟女老妇高清| 久久精品国产鲁丝片午夜精品| 黄色配什么色好看| 中文欧美无线码| 自拍偷自拍亚洲精品老妇| 精品久久久久久久久av| 成年人午夜在线观看视频| 老师上课跳d突然被开到最大视频| 美女福利国产在线 | 国产欧美日韩一区二区三区在线 | 久久国产亚洲av麻豆专区| 亚洲精品日本国产第一区| 亚洲图色成人| 蜜桃在线观看..| 国产探花极品一区二区| 国产av码专区亚洲av| 国产精品福利在线免费观看| 亚洲精品色激情综合| 久久97久久精品| 91精品国产九色| 在线精品无人区一区二区三 | 一区二区三区免费毛片| 99热这里只有是精品在线观看| 国产伦精品一区二区三区四那| 人人妻人人澡人人爽人人夜夜| 欧美日韩在线观看h| 国产亚洲午夜精品一区二区久久| 插阴视频在线观看视频| 九九爱精品视频在线观看| 日韩在线高清观看一区二区三区| 欧美精品人与动牲交sv欧美| 99热网站在线观看| 亚洲,一卡二卡三卡| 欧美成人精品欧美一级黄| 18+在线观看网站| 蜜臀久久99精品久久宅男| 噜噜噜噜噜久久久久久91| 一级毛片我不卡| 久久久久久久大尺度免费视频| 精品亚洲成国产av| 91精品一卡2卡3卡4卡| 大陆偷拍与自拍| 久久99热这里只有精品18| 少妇裸体淫交视频免费看高清| 色婷婷久久久亚洲欧美| av线在线观看网站| 亚洲中文av在线| 久久精品人妻少妇| 亚洲国产精品一区三区| 亚洲精品日本国产第一区| 国产乱人偷精品视频| 午夜视频国产福利| 最近中文字幕2019免费版| 观看av在线不卡| 久久99热这里只频精品6学生| 婷婷色av中文字幕| 亚洲最大成人中文| 日韩人妻高清精品专区| 色哟哟·www| 午夜激情福利司机影院| 国产一区有黄有色的免费视频| 性高湖久久久久久久久免费观看| 亚洲精品久久久久久婷婷小说| 亚洲国产欧美人成| 亚洲va在线va天堂va国产| 国产在线一区二区三区精| 亚洲,欧美,日韩| 中文字幕av成人在线电影| xxx大片免费视频| 水蜜桃什么品种好| 成人综合一区亚洲| 成年美女黄网站色视频大全免费 | 亚洲不卡免费看| 中国美白少妇内射xxxbb| 老司机影院毛片| 国产精品无大码| 免费观看性生交大片5| 久久久久久久久久久丰满| 国产探花极品一区二区| 免费观看a级毛片全部| 色视频在线一区二区三区| av女优亚洲男人天堂| 欧美一级a爱片免费观看看| 国产免费一级a男人的天堂| 丰满乱子伦码专区| 亚洲人与动物交配视频| 久久久精品免费免费高清| 久久精品人妻少妇| 青青草视频在线视频观看| 欧美xxxx性猛交bbbb| 亚洲精华国产精华液的使用体验| 日本色播在线视频| 卡戴珊不雅视频在线播放| 欧美日韩亚洲高清精品| 国产精品一区二区在线观看99| 国产真实伦视频高清在线观看| 国产精品女同一区二区软件| 水蜜桃什么品种好| 久久精品熟女亚洲av麻豆精品| 日产精品乱码卡一卡2卡三| 久久97久久精品| 一级av片app| 一级毛片aaaaaa免费看小| 中文天堂在线官网| 搡老乐熟女国产| 97精品久久久久久久久久精品| 国产午夜精品一二区理论片| 亚洲精品aⅴ在线观看| 麻豆国产97在线/欧美| 美女高潮的动态| 尾随美女入室| 日本色播在线视频| 精品久久久精品久久久| 老熟女久久久| 大话2 男鬼变身卡| 一级毛片电影观看| 国产69精品久久久久777片| 亚洲精品一二三| 日韩av在线免费看完整版不卡| 18禁在线播放成人免费| 麻豆精品久久久久久蜜桃| 日韩成人伦理影院| 久久精品人妻少妇| 成人18禁高潮啪啪吃奶动态图 | 成人高潮视频无遮挡免费网站| 婷婷色综合大香蕉| 春色校园在线视频观看| 免费人成在线观看视频色| 国产女主播在线喷水免费视频网站| 午夜福利网站1000一区二区三区| 五月玫瑰六月丁香| 亚洲欧美清纯卡通| 99热这里只有是精品50| 美女脱内裤让男人舔精品视频| 久久久久久久久大av| 日韩免费高清中文字幕av| 精品亚洲乱码少妇综合久久| 亚洲电影在线观看av| 国产成人精品婷婷| 中文字幕精品免费在线观看视频 | 成人高潮视频无遮挡免费网站| av在线观看视频网站免费| 久久ye,这里只有精品| 少妇 在线观看| 日韩成人伦理影院| 亚洲av.av天堂| 性色avwww在线观看| 涩涩av久久男人的天堂| 插逼视频在线观看| 一本色道久久久久久精品综合| 乱系列少妇在线播放| 亚洲av中文字字幕乱码综合| 成人毛片a级毛片在线播放| 国产亚洲欧美精品永久| 少妇 在线观看| 欧美+日韩+精品| www.色视频.com| 日韩亚洲欧美综合| 日韩欧美一区视频在线观看 | 亚洲精品久久久久久婷婷小说| 亚洲精品久久午夜乱码| 熟妇人妻不卡中文字幕| 国产熟女欧美一区二区| 欧美精品国产亚洲| 黄色视频在线播放观看不卡| 国产精品不卡视频一区二区| 美女主播在线视频| 欧美日韩国产mv在线观看视频 | 在线观看三级黄色| 久久精品夜色国产| 人妻制服诱惑在线中文字幕| 国产精品一及| 久久久久国产精品人妻一区二区| 国产精品一区二区三区四区免费观看| 亚洲精品色激情综合| 麻豆精品久久久久久蜜桃| 三级经典国产精品| 精品久久久久久久久av| 欧美 日韩 精品 国产| 高清午夜精品一区二区三区| 国产精品爽爽va在线观看网站| 国产精品蜜桃在线观看| 在线看a的网站| 国产高清有码在线观看视频| 日韩欧美一区视频在线观看 | 免费av中文字幕在线| 成年av动漫网址| 亚洲精品aⅴ在线观看| 亚洲色图综合在线观看| 国产大屁股一区二区在线视频| 精品国产乱码久久久久久小说| 色视频www国产| 久久热精品热| av不卡在线播放| 国产精品一区二区在线观看99| 91在线精品国自产拍蜜月| 亚洲欧美成人精品一区二区| 97在线视频观看| 亚洲第一av免费看| 国产精品一及| 久久久久久久久久久丰满| 日韩中字成人| 天堂俺去俺来也www色官网| 亚洲精品aⅴ在线观看| 久久国内精品自在自线图片| h日本视频在线播放| 一区二区三区四区激情视频| 亚洲久久久国产精品| 亚洲av电影在线观看一区二区三区| 日韩电影二区| 黑人高潮一二区| 国产日韩欧美亚洲二区| 日韩不卡一区二区三区视频在线| 女的被弄到高潮叫床怎么办| 深夜a级毛片| 人妻系列 视频| 亚洲av国产av综合av卡| 乱系列少妇在线播放| 成人毛片60女人毛片免费| 大香蕉97超碰在线| 国产精品久久久久久久电影| 国产免费一区二区三区四区乱码| 国产黄片美女视频| 精品少妇久久久久久888优播| 国产成人一区二区在线| 久久ye,这里只有精品| 国产在线男女| 免费av中文字幕在线| 99久久精品一区二区三区| 又粗又硬又长又爽又黄的视频| 国产高清三级在线| 亚洲一级一片aⅴ在线观看| 国产日韩欧美亚洲二区| 美女国产视频在线观看| 国产乱人偷精品视频| 日本wwww免费看| 国产欧美日韩一区二区三区在线 | 国产亚洲5aaaaa淫片| 99热6这里只有精品| 22中文网久久字幕| 国产精品久久久久久精品电影小说 | 搡老乐熟女国产| 18+在线观看网站| 久久ye,这里只有精品| 国产又色又爽无遮挡免| 搡老乐熟女国产| 国产精品不卡视频一区二区| 在线观看国产h片| 久久亚洲国产成人精品v| 亚洲av电影在线观看一区二区三区| 亚洲精品一二三| 亚洲av电影在线观看一区二区三区| 在线播放无遮挡| 欧美性感艳星| 91狼人影院| 综合色丁香网| 国产一区二区三区综合在线观看 | 大话2 男鬼变身卡| 在线 av 中文字幕| 中文天堂在线官网| 狂野欧美激情性bbbbbb| 亚洲精品亚洲一区二区| 日韩一区二区三区影片| 国产精品成人在线| 蜜桃亚洲精品一区二区三区| 欧美日韩精品成人综合77777| 亚洲国产成人一精品久久久| 国内揄拍国产精品人妻在线| 亚洲不卡免费看| 精品一区二区免费观看| 亚洲精品一区蜜桃| 亚洲av.av天堂| 亚洲av成人精品一区久久| 欧美精品国产亚洲| 人妻 亚洲 视频| 国产精品爽爽va在线观看网站| 熟女av电影| 一级二级三级毛片免费看| 欧美激情国产日韩精品一区| 久久青草综合色| 国产在线免费精品| 18禁在线无遮挡免费观看视频| 成人亚洲欧美一区二区av| 欧美xxxx黑人xx丫x性爽| 你懂的网址亚洲精品在线观看| 毛片女人毛片| 寂寞人妻少妇视频99o| 在线 av 中文字幕| 久久精品人妻少妇| 久久99热6这里只有精品| 五月伊人婷婷丁香| 免费黄网站久久成人精品| 精品久久久久久电影网| 又黄又爽又刺激的免费视频.| 中文精品一卡2卡3卡4更新| 夜夜爽夜夜爽视频| 黑人猛操日本美女一级片| 欧美日本视频| 亚洲精品久久午夜乱码| 一个人看的www免费观看视频| 插阴视频在线观看视频| 成人毛片60女人毛片免费| 免费黄色在线免费观看| 国产精品无大码| 久久久久久久国产电影| 亚洲av成人精品一二三区| 精品亚洲成国产av| 一边亲一边摸免费视频| 欧美一级a爱片免费观看看| 亚洲欧洲日产国产| 久久av网站| av国产精品久久久久影院| 夫妻性生交免费视频一级片| 毛片女人毛片| 亚洲精品自拍成人| 免费看日本二区| 亚洲怡红院男人天堂| 国产一区有黄有色的免费视频| 精品久久国产蜜桃| 亚洲国产精品999| 久久综合国产亚洲精品| 国产精品不卡视频一区二区| 少妇的逼好多水| 亚洲国产精品国产精品| 欧美精品一区二区免费开放| 一级毛片我不卡| 久久久欧美国产精品| 国产视频首页在线观看| 人妻 亚洲 视频| 七月丁香在线播放| 中文字幕制服av| 欧美日韩视频精品一区| 国产精品国产三级国产专区5o| 成人毛片60女人毛片免费| 天美传媒精品一区二区| 国产精品人妻久久久久久| 国国产精品蜜臀av免费| 夜夜爽夜夜爽视频| 久久久久视频综合| 高清在线视频一区二区三区| 精品国产乱码久久久久久小说| 啦啦啦啦在线视频资源| 国产精品不卡视频一区二区| 22中文网久久字幕| 美女cb高潮喷水在线观看| 日韩国内少妇激情av| 夫妻午夜视频| 久久精品国产a三级三级三级| 欧美日韩在线观看h| 一级黄片播放器| 亚洲精品aⅴ在线观看| 大香蕉97超碰在线| 99精国产麻豆久久婷婷| 六月丁香七月| 黄色一级大片看看| 久久久久久久精品精品| 欧美老熟妇乱子伦牲交| 伦理电影免费视频| 亚洲第一av免费看| 精品午夜福利在线看| av在线老鸭窝| 最近最新中文字幕免费大全7| 亚洲av在线观看美女高潮| 国内揄拍国产精品人妻在线| 纯流量卡能插随身wifi吗| 日本色播在线视频| 国产精品久久久久久精品古装| 国产黄频视频在线观看| 青春草国产在线视频| 亚洲欧美日韩无卡精品| 涩涩av久久男人的天堂| 搡老乐熟女国产| h日本视频在线播放| av女优亚洲男人天堂| 亚洲熟女精品中文字幕| 亚洲av成人精品一区久久| 久久久久人妻精品一区果冻| 成人亚洲欧美一区二区av| 亚洲人与动物交配视频| 亚洲欧美成人综合另类久久久| 日韩电影二区| 国产av国产精品国产| 国产精品一区二区性色av| 超碰97精品在线观看| 最近最新中文字幕免费大全7| 欧美高清成人免费视频www| 欧美成人a在线观看| 22中文网久久字幕| 熟妇人妻不卡中文字幕| 多毛熟女@视频| 精品酒店卫生间| 国产淫片久久久久久久久| 成人黄色视频免费在线看| 一个人免费看片子| 青青草视频在线视频观看| 全区人妻精品视频| 色视频www国产| 十分钟在线观看高清视频www | 亚洲在久久综合| 一本久久精品| 黄片无遮挡物在线观看| 涩涩av久久男人的天堂| 成人漫画全彩无遮挡| 人体艺术视频欧美日本| 免费黄色在线免费观看| 国产成人91sexporn| 久热这里只有精品99| 国产亚洲最大av| 国产av一区二区精品久久 | 日韩欧美一区视频在线观看 | 美女cb高潮喷水在线观看| 自拍偷自拍亚洲精品老妇| 最近最新中文字幕大全电影3| 18禁在线播放成人免费| 肉色欧美久久久久久久蜜桃| 国产无遮挡羞羞视频在线观看| 亚洲精品色激情综合| 精品少妇黑人巨大在线播放| 我要看黄色一级片免费的| 国产午夜精品久久久久久一区二区三区| 最近最新中文字幕大全电影3| 久久人人爽av亚洲精品天堂 | 人妻系列 视频| 高清日韩中文字幕在线| 人人妻人人爽人人添夜夜欢视频 | 国产69精品久久久久777片| 免费不卡的大黄色大毛片视频在线观看| 国产成人一区二区在线| 热99国产精品久久久久久7| 一区在线观看完整版| 成人一区二区视频在线观看| av卡一久久| 91午夜精品亚洲一区二区三区| 亚洲av在线观看美女高潮| 激情 狠狠 欧美| 国产精品久久久久久久电影| 免费看日本二区| 国产男女内射视频| 亚洲成色77777| 麻豆乱淫一区二区| 黑人高潮一二区| 黄色配什么色好看| 狂野欧美激情性bbbbbb| 亚洲av男天堂| 99re6热这里在线精品视频| 中文乱码字字幕精品一区二区三区| 少妇丰满av| 国产成人午夜福利电影在线观看| 日韩免费高清中文字幕av| 深夜a级毛片| 久久婷婷青草| 欧美成人一区二区免费高清观看| 国产黄频视频在线观看| 国产免费又黄又爽又色| 在线观看免费高清a一片| 日日啪夜夜撸| 成人美女网站在线观看视频| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品乱码久久久久久按摩| 乱码一卡2卡4卡精品| 久久久久久久精品精品| 精品人妻偷拍中文字幕| 国产日韩欧美亚洲二区| 久久99热这里只频精品6学生| 国内少妇人妻偷人精品xxx网站| 亚洲精品乱码久久久v下载方式| 青春草亚洲视频在线观看| 国产免费福利视频在线观看| 蜜桃久久精品国产亚洲av| 我要看日韩黄色一级片| 亚洲欧美一区二区三区黑人 | 你懂的网址亚洲精品在线观看| 美女视频免费永久观看网站| 男女边吃奶边做爰视频| 99视频精品全部免费 在线| 嘟嘟电影网在线观看| 欧美另类一区| 亚洲电影在线观看av| 久久久久久久国产电影| 久久影院123| 两个人的视频大全免费| 老熟女久久久| 午夜免费鲁丝| 天天躁夜夜躁狠狠久久av| 老师上课跳d突然被开到最大视频| 久久青草综合色| 亚洲一区二区三区欧美精品| 久久久久久久亚洲中文字幕| 国产成人免费无遮挡视频| 久久99精品国语久久久|