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

    同軸螺線管身管膛內(nèi)電磁感應(yīng)特性仿真與優(yōu)化

    2018-06-12 08:00:10毛保全蘭圖鄧威
    現(xiàn)代電子技術(shù) 2018年11期
    關(guān)鍵詞:螺線管身管電磁感應(yīng)

    毛保全 蘭圖 鄧威

    摘 要: 針對(duì)密繞螺線管在高頻外部激勵(lì)下感應(yīng)磁場難以穿透金屬圓筒壁的問題,基于麥克斯韋方程組推導(dǎo)不同激勵(lì)類型的通電螺線管金屬圓筒電磁感應(yīng)特性,并以某型火炮身管為仿真算例,建立同軸螺線管身管電磁仿真模型,對(duì)不同激勵(lì)頻率下膛內(nèi)電磁感應(yīng)特性進(jìn)行仿真分析。為進(jìn)一步提高交流激勵(lì)下膛內(nèi)磁感應(yīng)強(qiáng)度,以線圈參數(shù)、身管參數(shù)為設(shè)計(jì)變量,對(duì)同軸螺線管身管膛內(nèi)電磁感應(yīng)特性進(jìn)行優(yōu)化。仿真結(jié)果表明,膛內(nèi)磁感應(yīng)強(qiáng)度隨著激勵(lì)頻率的提高呈指數(shù)型降低;在相同激勵(lì)條件下,可通過優(yōu)化線圈、身管參數(shù)降低金屬身管對(duì)外激勵(lì)源的屏蔽效應(yīng),從而提高膛內(nèi)等效磁場強(qiáng)度。

    關(guān)鍵詞: 螺線管; 身管; 電磁感應(yīng); 優(yōu)化設(shè)計(jì); 激勵(lì)頻率; 等效磁場強(qiáng)度

    中圖分類號(hào): TN03?34; TP391.9 文獻(xiàn)標(biāo)識(shí)碼: A 文章編號(hào): 1004?373X(2018)11?0077?05

    Simulation and optimization of electromagnetic induction characteristic

    in chamber of coaxial solenoid′s barrel

    MAO Baoquan, LAN Tu, DENG Wei

    (Department of Weapon Engineering, Academy of Armored Force Engineering, Beijing 100072, China)

    Abstract: Since the induction magnetic field of close?coiled solenoid is difficult to penetrate the metal cylinder wall under the high frequency external excitation, the electromagnetic induction characteristic of the metal cylinder of electrified solenoid with different excitation types is deduced on the basis of Maxwell′s equations, and the electromagnetic simulation model of the barrel of coaxial solenoid is established by taking a certain gun barrel as the simulation example to perform the simulation analysis of the electromagnetic induction characteristic in the chamber at different excitation frequencies. The coil parameter and barrel parameter are taken as the design variables to optimize the electromagnetic induction characteristic in the chamber of coaxial solenoid′s barrel, so as to further improve the magnetic induction intensity in the chamber at AC excitation. The simulation results show that the magnetic induction intensity is exponentially decreased with the increase of the excitation frequency; under the same excitation conditions, the optimization of coil and barrel parameter can reduce the shielding effect of the metal barrel on external excitation source, so as to improve the equivalent magnetic field strength in the chamber.

    Keywords: solenoid; barrel; electromagnetic induction; optimization design; excitation frequency; equivalent magnetic field

    0 引 言

    施加磁場的圓筒腔體在軍事領(lǐng)域有著廣泛的應(yīng)用,如勵(lì)磁線圈火炮身管腔體產(chǎn)生的電磁場不僅能傳遞能量、提高身管發(fā)射彈丸的效能[1?2],而且能傳遞信息、提高身管的可檢測性[3],目前,已在火炮身管缺陷檢測、引信膛內(nèi)感應(yīng)儲(chǔ)能、電磁線圈炮發(fā)射等領(lǐng)域得到成功應(yīng)用。

    身管的激勵(lì)磁場可通過在身管上添加通電密繞螺線管實(shí)現(xiàn)。目前,關(guān)于螺線管的研究成果大部分集中在兩方面:一方面是僅針對(duì)密繞螺線管的單獨(dú)研究,主要研究螺線管的空間磁場分布特性[4?6];另一方面是將螺線管應(yīng)用于無損探傷領(lǐng)域,通過漏磁信號(hào)檢測金屬表面的缺陷[7?8]。然而,將螺線管與金屬身管結(jié)合并著眼于膛內(nèi)電磁感應(yīng)特性的研究還鮮有報(bào)道??紤]到炮鋼材質(zhì)的火炮身管在電磁場中可等效為金屬圓筒,在交變電流的激勵(lì)下其表面形成的渦流場勢(shì)必會(huì)削弱外激勵(lì)磁場,從而對(duì)膛內(nèi)等效磁場分布產(chǎn)生較大影響,因此,有必要對(duì)同軸螺線管身管膛內(nèi)電磁感應(yīng)特性進(jìn)行研究。本文基于麥克斯韋基本理論推導(dǎo)不同激勵(lì)類型下同軸螺線管身管膛內(nèi)磁場分布數(shù)學(xué)模型,以某型身管為算例,利用Ansoft Maxwell建立密繞螺線管身管的電磁場數(shù)值仿真模型,并對(duì)不同外激勵(lì)頻率下身管內(nèi)感應(yīng)磁場進(jìn)行仿真分析,其后固定激勵(lì)條件不變,通過優(yōu)化線圈參數(shù)及身管參數(shù)提高膛內(nèi)磁感應(yīng)強(qiáng)度。

    1 螺線管金屬圓筒電磁感應(yīng)特性

    1.1 通直流電螺線管金屬筒內(nèi)電磁感應(yīng)特性

    如圖1所示為螺線管線圈結(jié)構(gòu)簡圖,根據(jù)畢奧?薩伐爾定律,其單個(gè)載流圓環(huán)在空間任意點(diǎn)[P]產(chǎn)生的磁場為[9]:

    [B=μ04πIdl×RR3] (1)

    式中:[Idl]為圓環(huán)電流元;[R]為從電流元指向[P]點(diǎn)的矢量;[R]為電流元與[P]點(diǎn)間的距離;[ μ0]為真空磁導(dǎo)率;[B]為磁感應(yīng)強(qiáng)度;[I]為線圈中的電流。[dl],[R],[R]分別由下列公式計(jì)算得出:

    [dl=(-asinφi+acosφj)dφ] (2)

    [R=(rsinθ-acosφ)i-asinφj+rcosθk] (3)

    [R2=a2+r2-2arsinθcosφ] (4)

    式中:[a]為載流圓環(huán)半徑;[r]為[P]點(diǎn)到圓環(huán)圓心的距離;[θ]為[r]與[z]軸的夾角;[φ]為電流元與[x]軸的夾角。

    當(dāng)忽略電流的螺旋性以及線間距離時(shí),可視為多個(gè)載流圓環(huán)的疊加。設(shè)螺線管長度為[2l],線圈的匝數(shù)為[n]匝,則螺線管在[P]點(diǎn)產(chǎn)生的磁場為:

    [B=nμ04π-lldzIdl×RR3] (5)

    展開得:

    當(dāng)[xP=0]時(shí),可得螺線管軸線處的磁場為:

    [BzP=μ0nIa4π-lldz02πadφ(a2+(zP-z)2)3] (11)

    螺線管軸向、徑向磁場變化曲線分別如圖2,圖3所示。

    螺旋管軸線處的磁感應(yīng)強(qiáng)度在螺線管長度范圍內(nèi)緩慢變化,且關(guān)于[z=0]平面對(duì)稱,在[z=0]處最強(qiáng),在螺線管長度范圍外磁感應(yīng)強(qiáng)度迅速衰減;對(duì)于徑向磁場變化情況,磁感應(yīng)強(qiáng)度以螺線管中軸線為中心呈輻射對(duì)稱分布,即磁場接近軸線處較強(qiáng),距離軸線越遠(yuǎn)磁場強(qiáng)度越弱,在螺線管外部的磁場接近于零。

    1.2 通交流電螺線管金屬筒內(nèi)電磁感應(yīng)特性

    當(dāng)激勵(lì)源為正弦電流時(shí),螺線管中產(chǎn)生正弦變化的磁場,而變化的磁場將產(chǎn)生變化的渦流場,從而抵消原激勵(lì)磁場的變化,對(duì)整體磁場的分布造成影響,此時(shí)磁場[B]和電場[E]可表示為:

    [E=E0exp(-jωt)] (12)

    [B=B0exp(-jωt)] (13)

    考慮螺線管中有與螺線管半徑相同的金屬圓筒的情況,則金屬導(dǎo)體會(huì)產(chǎn)生渦旋電流。由于渦旋電流遠(yuǎn)大于位移電流,故位移電流可忽略不計(jì),可得:

    [?2H0-jωγμH0=0] (14)

    設(shè)螺線管線圈中交變電流角頻率為[ω],取柱坐標(biāo),設(shè)[k2=-jωγμ],則[x=ky],有:

    [d2Hdx2+1xdHdx+H=0] (15)

    式(5)為零階Bessel方程,其解為:

    [H=aJ0(kx)+bY0(kx)] (16)

    式中:[J0(kx)]為第一類Bessel函數(shù);[Y0(kx)]為第二類Bessel函數(shù)。

    由于[H0]在[x=0]處為有限值,而[Y0(0)=∞],因此,[b=0],故:

    [H0=aJ0(x)+aJ0(kr)] (17)

    該問題的邊界條件為:

    [H(x)=H0R, r=R] (18)

    [H0=H0RJ0(kR)J0(kr)] (19)

    [k]為復(fù)數(shù),在[J0(kr)]中展開,得:

    [Her(v)≡ReJ0(i-12v)] (20)

    [Hei(v)≡ImJ0(i-12v)] (21)

    [J0(i-12r)≡Her(v)+jHei(v)] (22)

    應(yīng)用上述關(guān)系式,則:

    [H(r)=H0Ber(rωμγ)+jBei(rωμγ)Ber(R2ωμγ)+jBei(R2ωμγ)] (23)

    令[HH0=HH0<θ],則有:

    [H(r)H0=B2er(rrωμγ)+B2ei(rrωμγ)B2er(R2rωμγ)+B2ei(R2rωμγ)12] (24)

    輻射角為:

    [θ=arctanBer(rωμγ)Bei(R2ωμγ)-Ber(rωμγ)Bei(R2ωμγ)Ber(R2ωμγ)Bei(rωμγ)+Ber(R2ωμγ)Bei(rωμγ)] (25)

    [HH0]為金屬圓筒內(nèi)磁場幅值與外激勵(lì)磁場幅值的比值,其在不同[ω]條件下隨半徑變化的關(guān)系曲線如圖4所示。

    由圖4可知,對(duì)于圓筒內(nèi)部,[ω]越小,圓筒內(nèi)部的磁場越接近均勻分布,隨著[ω]的增大,金屬圓筒內(nèi)部的渦流場對(duì)外激勵(lì)磁場的抵消作用越明顯。通過不同[ω]條件下磁感應(yīng)仿真曲線可以得出,當(dāng)[ω]=3 000 rad/s時(shí),在圓筒中心的磁場強(qiáng)度接近于零,此時(shí)[ω]外激勵(lì)磁場的諧振頻率為573 Hz。

    2 同軸螺線管身管膛內(nèi)電磁感應(yīng)特性

    2.1 仿真模型的建立

    以某型火炮身管為仿真算例,通過Maxwell構(gòu)建火炮身管模型及密繞線圈模型,以研究火炮身管加裝線圈后的電磁效應(yīng)。為提高仿真效率,火炮身管通過金屬圓筒進(jìn)行模擬,圓筒壁外密繞線圈形成螺線管,在Maxwell中建立外套密繞螺線管的火炮身管簡化模型,如圖5所示。

    圖5中,外部環(huán)狀為勵(lì)磁線圈,內(nèi)部圓筒為身管等效模型。模型中圓筒長度為600 mm,外徑及內(nèi)徑分別為40 mm,30 mm,;螺線管通過winding功能添加于身管外壁,位于孔長度的中間位置,其內(nèi)徑為35 mm;線圈外徑為50 mm,內(nèi)徑為40 mm,長度為300 mm;身管材料選用PCrNi3MoVA鋼,螺線管材料設(shè)置為銅;激勵(lì)設(shè)置于身管橫截面上。

    2.2 結(jié)果分析

    激勵(lì)和邊界條件設(shè)置完畢后進(jìn)行求解計(jì)算,得到模型的磁場強(qiáng)度矢量圖、磁力線圖及磁場強(qiáng)度分布圖,如圖6,圖7所示。

    由圖6,圖7可以看出,螺線管激勵(lì)磁場主要集中在螺線管內(nèi)部,在螺線管外部的磁感應(yīng)強(qiáng)度較低;外激勵(lì)在圓筒引發(fā)的渦流場反作用于外激勵(lì),大幅度削弱了螺線管內(nèi)部的磁感應(yīng)強(qiáng)度,導(dǎo)致磁場僅存在于金屬圓筒壁。以上結(jié)論表明,金屬圓筒的存在較大程度地改變了螺線管的空間磁場構(gòu)型,大大削弱了螺線管內(nèi)部的磁場強(qiáng)度,上述仿真結(jié)果進(jìn)一步驗(yàn)證了前述小節(jié)理論分析的正確性。

    2.3 不同激勵(lì)頻率對(duì)膛內(nèi)電磁感應(yīng)特性的影響

    為進(jìn)一步掌握不同外激勵(lì)頻率下螺線管金屬圓筒的磁感應(yīng)特性,分別以0,100 Hz,200 Hz,300 Hz,500 Hz等不同激勵(lì)頻率進(jìn)行仿真,結(jié)果如圖8所示。

    由圖8可知,當(dāng)外激勵(lì)頻率為0時(shí),圓筒腔體內(nèi)部感應(yīng)磁場沿軸線均勻分布,在圓筒端部迅速衰減;當(dāng)外激勵(lì)頻率不為0時(shí),圓筒內(nèi)部等效磁場隨外激勵(lì)頻率的提高而減弱,當(dāng)外激勵(lì)頻率為100 Hz時(shí),穿透金屬圓筒至軸線處的磁場強(qiáng)度為激勵(lì)磁場強(qiáng)度的61.9%,而激勵(lì)頻率達(dá)到500 Hz時(shí),穿透后的磁場僅為激勵(lì)磁場強(qiáng)度的3.4%。

    3 螺線管身管膛內(nèi)電磁感應(yīng)特性優(yōu)化

    為進(jìn)一步提高螺線管在身管內(nèi)部的激勵(lì)磁場強(qiáng)度,以線圈參數(shù)及身管參數(shù)為設(shè)計(jì)變量,對(duì)同軸螺線管金屬筒內(nèi)電磁感應(yīng)特性進(jìn)行優(yōu)化。線圈參數(shù)包括線圈密度、匝數(shù)及厚度,身管參數(shù)包括磁導(dǎo)率及電導(dǎo)率。通過優(yōu)化設(shè)計(jì)變量來提高軸線處磁場強(qiáng)度,優(yōu)化效果則可通過其與外激勵(lì)磁場強(qiáng)度的比值與基準(zhǔn)值的比較進(jìn)行量化表征。

    3.1 優(yōu)化模型

    根據(jù)上述仿真結(jié)果,以交變電流激勵(lì)下穿透金屬圓筒至軸線處的磁場強(qiáng)度為激勵(lì)磁場強(qiáng)度的50%作為基準(zhǔn)穿透比例,其對(duì)應(yīng)的線圈、身管參數(shù)作為初始條件,設(shè)計(jì)變量及取值范圍如表1所示。

    同時(shí),選擇螺線管激勵(lì)磁場與身管軸線處磁場的比值與基準(zhǔn)值(50%)的變化幅度作為表征電磁感應(yīng)特性優(yōu)化效果的特征量。

    [F=HmH0-0.50.5×100%] (26)

    3.2 優(yōu)化算法

    考慮到遺傳算法具有發(fā)展成熟、收斂速度快、不易陷入局部最優(yōu)解等優(yōu)點(diǎn)[10],因此,選擇其作為本文的優(yōu)化算法。遺傳算法將每個(gè)可行解等同于生物學(xué)個(gè)體,優(yōu)化前系統(tǒng)隨機(jī)將[N]個(gè)個(gè)體組成初始種群,然后人為確定適應(yīng)度函數(shù),系統(tǒng)則根據(jù)適應(yīng)度高低對(duì)個(gè)體進(jìn)行排序,并將其中高適應(yīng)度的個(gè)體進(jìn)行遺傳操作(選擇、交叉及變異),對(duì)低適應(yīng)度的個(gè)體進(jìn)行淘汰。該過程反復(fù)進(jìn)行,直至最終適應(yīng)度最好的個(gè)體被篩選出來,即為全局最優(yōu)解。尋優(yōu)過程如圖9所示。

    在本文仿真分析中,種群數(shù)量設(shè)置為44,交叉概率及變異概率分別為0.65,0.04,收斂閾值為0.01。

    3.3 優(yōu)化計(jì)算與結(jié)果分析

    對(duì)上述同軸螺線管身管膛內(nèi)電磁感應(yīng)特性優(yōu)化設(shè)計(jì)進(jìn)行仿真計(jì)算,結(jié)果如圖10,表2所示。

    由圖10可知,優(yōu)化過程經(jīng)過732次迭代后收斂到全局最優(yōu)解。由表2可知,全局最優(yōu)解為匝數(shù)1 769匝,密度0.03匝/mm2,厚度5.04 mm,相對(duì)磁導(dǎo)率712,電導(dǎo)率1.44×106 S/m,介電常數(shù)7.25×10-12 F/m,表明在該最優(yōu)參數(shù)組合下,同軸螺線管身管膛內(nèi)軸線處磁場強(qiáng)度為外激勵(lì)磁場強(qiáng)度的63%,相比初始值優(yōu)化幅度為26%,優(yōu)化效果較為顯著。

    4 結(jié) 論

    通過本文研究可得出如下結(jié)論:

    1) 金屬圓筒在交變電流的激勵(lì)下會(huì)產(chǎn)生渦流,削弱外激勵(lì)穿透至金屬圓筒的磁場強(qiáng)度,從而對(duì)外部磁場產(chǎn)生屏蔽效應(yīng)。

    2) 金屬身管對(duì)同軸螺線管激勵(lì)磁場的屏蔽效應(yīng)隨著激勵(lì)頻率的增加而劇烈,穿透后的磁場強(qiáng)度隨頻率的提高而呈指數(shù)規(guī)律衰減。仿真結(jié)果表明,若使穿透后的磁場強(qiáng)度不小于激勵(lì)磁場強(qiáng)度的61.9%,需將外激勵(lì)頻率控制在100 Hz以下。

    3) 在相同激勵(lì)條件下,通過優(yōu)化線圈參數(shù)及身管參數(shù),能夠使同軸螺線管身管膛內(nèi)磁場強(qiáng)度得到提高,對(duì)于軸線磁場強(qiáng)度為外激勵(lì)磁場強(qiáng)度50%的原有激勵(lì)條件,通過調(diào)節(jié)身管參數(shù)及線圈參數(shù),能夠使該穿透比例優(yōu)化至63%。

    參考文獻(xiàn)

    [1] 李軍,嚴(yán)萍,袁偉群.電磁軌道炮發(fā)射技術(shù)的發(fā)展與現(xiàn)狀[J].高電壓技術(shù),2014,40(4):1052?1064.

    LI Jun, YAN Ping, YUAN Weiqun. Electromagnetic gun technology and its development [J]. High voltage engineering, 2014, 40(4): 1052?1064.

    [2] 吳志亮,張合.小口徑引信膛內(nèi)感應(yīng)儲(chǔ)能磁場穿透特性[J].兵工學(xué)報(bào),2010,31(10):1310?1315.

    WU Zhiliang, ZHANG He. The magnetic field penetration characteristics in bore of small caliber fuze using inductive energy [J]. Acta armamentarii, 2010, 31(10): 1310?1315.

    [3] 郭希玲.交變漏磁檢測技術(shù)及螺紋缺陷識(shí)別的研究[D].長沙:國防科學(xué)技術(shù)大學(xué),2005.

    GUO Xiling. Research on technology of alternative current magnetic field leakage and flaw identification to the whorl [D]. Changsha: National University of Defense Technology, 2005.

    [4] 馮濤.有限長螺線管場分布特性分析[D].西安:西安電子科技大學(xué),2013.

    FENG Tao. Analysis on field distribution characteristics of limited solenoid [D]. Xian: Xidian University, 2013.

    [5] 廖斌,鄧春鳳,吳先映,等.利用Matlab計(jì)算螺線管內(nèi)磁場分布研究[J].北京師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2010,46(6):688?690.

    LIAO Bin, DENG Chunfeng, WU Xianying, et al. Inner magnetic field distribution of solenoid studied by Matlab [J]. Journal of Beijing Normal University (natural science), 2010, 46(6): 688?690.

    [6] 蔡偉,伍樊成,楊志勇,等.基于麥克斯韋方程的交變電流長螺線管磁場[J].強(qiáng)激光與粒子束,2015,27(12):162?167.

    CAI Wei, WU Fancheng, YANG Zhiyong, et al. Analysis of solenoid magnetic field based on Maxwell equation [J]. High power laser and particle beams, 2015, 27(12): 162?167.

    [7] 黃平.低頻電磁技術(shù)在儲(chǔ)罐罐底缺陷檢測中的應(yīng)用研究[D].沈陽:沈陽工業(yè)大學(xué),2016.

    HUANG Ping. Application research of low frequency electromagnetic technology in defect detection of tank bottom [D]. Shenyang: Shenyang University of Technology, 2016.

    [8] 吳德會(huì),游德海,柳振涼,等.交流漏磁檢測法趨膚深度的機(jī)理與實(shí)驗(yàn)研究[J].儀器儀表學(xué)報(bào), 2014,35(2):327?336.

    WU Dehui, YOU Dehai, LIU Zhenliang, et al. Mechanism and experimental research on skin depth in AC magnetic flux leakage testing [J]. Chinese journal of scientific instrument, 2014, 35(2): 327?336.

    [9] 蔡旭紅,李邵輝.有限長通電螺線管內(nèi)部空間磁場的模擬[J].汕頭大學(xué)學(xué)報(bào)(自然科學(xué)版),2004(2):28?31.

    CAI Xuhong, LI Shaohui. Simulation for spatial magnetic field distribution in a finite solenoid [J]. Journal of Shantou University (natural science), 2004(2): 28?31.

    [10] 梁傳建,楊國來,王曉鋒.基于神經(jīng)網(wǎng)絡(luò)和遺傳算法的火炮結(jié)構(gòu)動(dòng)力學(xué)優(yōu)化[J].兵工學(xué)報(bào),2015,36(5):789?794.

    LIANG Chuanjian, YANG Guolai, WANG Xiaofeng. Structural dynamics optimization of gun based on neural networks and genetic algorithms [J]. Acta armamentarii, 2015, 36(5): 789?794.

    猜你喜歡
    螺線管身管電磁感應(yīng)
    火炮身管壽命評(píng)估方法及其標(biāo)準(zhǔn)研究
    如何判斷通電螺線管的磁極
    千姿百態(tài)說電磁 歷久彌新話感應(yīng)——遵循“三步法”,搞定電磁感應(yīng)綜合題
    基于垂向固有振動(dòng)仿真的身管口徑倍數(shù)設(shè)計(jì)研究
    聚焦動(dòng)量觀點(diǎn)在電磁感應(yīng)中的應(yīng)用
    電磁感應(yīng)在生活中的應(yīng)用
    電子制作(2018年14期)2018-08-21 01:38:38
    基于數(shù)值仿真的四種類型身管強(qiáng)度比較
    無限長密繞通電螺線管磁場簡明計(jì)算
    安培定則之我見
    抓“四量”破電磁感應(yīng)綜合題お
    老司机深夜福利视频在线观看| 老熟妇仑乱视频hdxx| 天天躁日日操中文字幕| 综合色av麻豆| 国产精品日韩av在线免费观看| 国产三级黄色录像| 淫秽高清视频在线观看| 69人妻影院| 中亚洲国语对白在线视频| 午夜福利在线观看吧| 在线观看午夜福利视频| 午夜精品久久久久久毛片777| 首页视频小说图片口味搜索| 中文在线观看免费www的网站| 校园春色视频在线观看| 亚洲欧美激情综合另类| 亚洲色图av天堂| 一级毛片高清免费大全| 国产激情偷乱视频一区二区| 免费在线观看日本一区| 18+在线观看网站| 国产黄片美女视频| 1000部很黄的大片| 精品乱码久久久久久99久播| 欧美激情在线99| 欧美zozozo另类| 亚洲av成人av| 精品乱码久久久久久99久播| 19禁男女啪啪无遮挡网站| 叶爱在线成人免费视频播放| av视频在线观看入口| 欧美乱妇无乱码| 中文在线观看免费www的网站| 国产精品98久久久久久宅男小说| 亚洲自拍偷在线| 精品不卡国产一区二区三区| 色老头精品视频在线观看| 好看av亚洲va欧美ⅴa在| 国产黄片美女视频| 一进一出抽搐gif免费好疼| 久久久久亚洲av毛片大全| 日本 av在线| 亚洲精品国产精品久久久不卡| 岛国在线免费视频观看| 18禁黄网站禁片午夜丰满| 99久久综合精品五月天人人| 欧美日韩亚洲国产一区二区在线观看| 欧美色视频一区免费| 亚洲人成电影免费在线| 国产成人aa在线观看| 岛国在线免费视频观看| 亚洲av成人不卡在线观看播放网| 久久精品夜夜夜夜夜久久蜜豆| 91麻豆精品激情在线观看国产| 亚洲18禁久久av| 亚洲狠狠婷婷综合久久图片| 美女cb高潮喷水在线观看| 99热这里只有精品一区| 午夜日韩欧美国产| 精品国产亚洲在线| 亚洲国产精品999在线| 婷婷六月久久综合丁香| 一级a爱片免费观看的视频| 99热这里只有是精品50| 人妻丰满熟妇av一区二区三区| 国产av一区在线观看免费| 日本在线视频免费播放| 丰满人妻熟妇乱又伦精品不卡| 啦啦啦观看免费观看视频高清| 欧美乱码精品一区二区三区| 亚洲一区二区三区色噜噜| 国产一区二区在线观看日韩 | 香蕉久久夜色| 天堂网av新在线| 国产精品野战在线观看| 日韩国内少妇激情av| 婷婷精品国产亚洲av| 午夜福利在线在线| www.色视频.com| 国产精品爽爽va在线观看网站| 免费观看人在逋| 午夜福利在线观看吧| 亚洲欧美激情综合另类| 一本久久中文字幕| 国产视频内射| 在线播放国产精品三级| 9191精品国产免费久久| 神马国产精品三级电影在线观看| 一级黄片播放器| 国产视频内射| 欧美黄色淫秽网站| 免费看日本二区| 国产高清videossex| 国产v大片淫在线免费观看| 成人18禁在线播放| 日韩欧美在线二视频| 久久中文看片网| 国产亚洲欧美在线一区二区| 在线免费观看的www视频| 国产精品影院久久| 老鸭窝网址在线观看| 岛国在线免费视频观看| 日本撒尿小便嘘嘘汇集6| 国内精品久久久久精免费| 日韩精品中文字幕看吧| 久久久久久九九精品二区国产| 色综合亚洲欧美另类图片| 叶爱在线成人免费视频播放| 午夜福利视频1000在线观看| 在线看三级毛片| 亚洲久久久久久中文字幕| 丁香六月欧美| 免费av不卡在线播放| 淫秽高清视频在线观看| 母亲3免费完整高清在线观看| 99热6这里只有精品| 欧美国产日韩亚洲一区| 真人做人爱边吃奶动态| ponron亚洲| 亚洲国产欧洲综合997久久,| 欧美黄色淫秽网站| 日韩欧美在线乱码| 午夜两性在线视频| 女人高潮潮喷娇喘18禁视频| 十八禁网站免费在线| www.999成人在线观看| 亚洲,欧美精品.| 男女视频在线观看网站免费| 特大巨黑吊av在线直播| 国产极品精品免费视频能看的| 国产av不卡久久| 亚洲人成伊人成综合网2020| 又黄又粗又硬又大视频| 欧美午夜高清在线| 亚洲国产色片| 欧美激情久久久久久爽电影| 色综合亚洲欧美另类图片| 久久久色成人| 岛国在线观看网站| 麻豆国产97在线/欧美| aaaaa片日本免费| 成年免费大片在线观看| 亚洲欧美日韩东京热| 色综合亚洲欧美另类图片| 嫁个100分男人电影在线观看| 亚洲人成网站在线播| 色综合婷婷激情| ponron亚洲| 久久久久久久午夜电影| 午夜激情福利司机影院| 夜夜躁狠狠躁天天躁| 国产免费男女视频| 日本免费a在线| 免费人成视频x8x8入口观看| 十八禁人妻一区二区| 麻豆成人午夜福利视频| 国产精品永久免费网站| 老熟妇乱子伦视频在线观看| 亚洲国产精品999在线| 国产欧美日韩精品亚洲av| 国产麻豆成人av免费视频| 国产探花极品一区二区| 国产午夜精品久久久久久一区二区三区 | 日韩欧美 国产精品| 国产成人福利小说| 内地一区二区视频在线| 久久精品91蜜桃| 精品乱码久久久久久99久播| 18禁美女被吸乳视频| 在线看三级毛片| 日韩精品中文字幕看吧| 欧美日韩中文字幕国产精品一区二区三区| 国产精品永久免费网站| 亚洲成人精品中文字幕电影| 色精品久久人妻99蜜桃| 在线播放无遮挡| 精华霜和精华液先用哪个| 免费看十八禁软件| 国产精品亚洲美女久久久| 久久99热这里只有精品18| 一本一本综合久久| h日本视频在线播放| 九色成人免费人妻av| 又爽又黄无遮挡网站| 久久国产精品影院| 日韩欧美免费精品| 窝窝影院91人妻| 搡老岳熟女国产| 日韩欧美精品免费久久 | 欧美三级亚洲精品| 亚洲激情在线av| 亚洲精华国产精华精| 久久久国产成人精品二区| 国产伦精品一区二区三区四那| 亚洲欧美日韩无卡精品| www.www免费av| 亚洲午夜理论影院| 天堂动漫精品| 亚洲不卡免费看| 中文字幕av成人在线电影| 欧美zozozo另类| 国产探花在线观看一区二区| 黄色女人牲交| 国产单亲对白刺激| 日本免费一区二区三区高清不卡| 日韩欧美精品免费久久 | 波多野结衣高清无吗| 色播亚洲综合网| 国产精品野战在线观看| 色吧在线观看| 校园春色视频在线观看| 免费看a级黄色片| 18禁裸乳无遮挡免费网站照片| 亚洲午夜理论影院| 亚洲,欧美精品.| 国产97色在线日韩免费| 日韩欧美 国产精品| 亚洲aⅴ乱码一区二区在线播放| 国产色婷婷99| 国产精品久久久人人做人人爽| 中文字幕精品亚洲无线码一区| 波野结衣二区三区在线 | 黑人欧美特级aaaaaa片| 在线视频色国产色| 亚洲人与动物交配视频| 熟女电影av网| av在线天堂中文字幕| 久久国产精品人妻蜜桃| 在线观看一区二区三区| 日本撒尿小便嘘嘘汇集6| 国产老妇女一区| 亚洲狠狠婷婷综合久久图片| 国内精品久久久久精免费| 国产一区二区在线观看日韩 | 一级黄色大片毛片| 亚洲七黄色美女视频| 日韩大尺度精品在线看网址| 黄色女人牲交| 精品日产1卡2卡| 每晚都被弄得嗷嗷叫到高潮| 俺也久久电影网| 久久精品亚洲精品国产色婷小说| 国产成年人精品一区二区| 国产单亲对白刺激| 每晚都被弄得嗷嗷叫到高潮| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品一区二区www| 丰满人妻一区二区三区视频av | 成人av一区二区三区在线看| 欧美成人一区二区免费高清观看| 欧美精品啪啪一区二区三区| 久久香蕉国产精品| 国产久久久一区二区三区| 久久久久亚洲av毛片大全| 十八禁人妻一区二区| 欧美+日韩+精品| 变态另类丝袜制服| 精品久久久久久久久久免费视频| 日韩人妻高清精品专区| 99久久精品热视频| 国产黄片美女视频| 亚洲欧美日韩无卡精品| 久久6这里有精品| 首页视频小说图片口味搜索| 国产亚洲精品久久久久久毛片| 99在线人妻在线中文字幕| 桃色一区二区三区在线观看| 天天一区二区日本电影三级| 欧美黄色片欧美黄色片| 免费在线观看成人毛片| 看黄色毛片网站| av国产免费在线观看| 久久久色成人| 无限看片的www在线观看| 91av网一区二区| 欧美zozozo另类| 免费在线观看日本一区| 精品人妻1区二区| 午夜亚洲福利在线播放| 亚洲精品乱码久久久v下载方式 | 亚洲欧美日韩卡通动漫| 欧美极品一区二区三区四区| 久久精品亚洲精品国产色婷小说| 亚洲欧美日韩无卡精品| 51午夜福利影视在线观看| 欧美极品一区二区三区四区| 亚洲av二区三区四区| 免费电影在线观看免费观看| 久99久视频精品免费| 国产乱人伦免费视频| 免费看美女性在线毛片视频| 91在线观看av| 一个人免费在线观看的高清视频| 成人精品一区二区免费| 国产成人福利小说| 一级毛片女人18水好多| 日本精品一区二区三区蜜桃| 在线播放国产精品三级| 搡女人真爽免费视频火全软件 | av在线蜜桃| 成人无遮挡网站| 国内精品一区二区在线观看| 亚洲av第一区精品v没综合| 看黄色毛片网站| АⅤ资源中文在线天堂| 亚洲人成网站高清观看| 日韩欧美三级三区| 国产精品一区二区免费欧美| 国产伦精品一区二区三区视频9 | 欧美av亚洲av综合av国产av| 久久性视频一级片| 搡老熟女国产l中国老女人| 看免费av毛片| 久久久久亚洲av毛片大全| 成人特级av手机在线观看| 成人性生交大片免费视频hd| 怎么达到女性高潮| 日韩免费av在线播放| 深爱激情五月婷婷| 成熟少妇高潮喷水视频| 小说图片视频综合网站| 免费电影在线观看免费观看| 国产中年淑女户外野战色| 身体一侧抽搐| 午夜精品在线福利| 我的老师免费观看完整版| 国产成人影院久久av| 成年版毛片免费区| 最新在线观看一区二区三区| 欧美区成人在线视频| 欧美在线一区亚洲| 欧美高清成人免费视频www| 床上黄色一级片| 亚洲人成网站高清观看| 亚洲国产欧洲综合997久久,| 国产色婷婷99| 亚洲,欧美精品.| 久久伊人香网站| 久久精品国产自在天天线| 又爽又黄无遮挡网站| 最好的美女福利视频网| 中文字幕精品亚洲无线码一区| 日韩欧美 国产精品| 久久九九热精品免费| 九九在线视频观看精品| 久久性视频一级片| 身体一侧抽搐| 男人和女人高潮做爰伦理| 精品国产亚洲在线| 久久久色成人| 老熟妇仑乱视频hdxx| 亚洲欧美日韩高清专用| 免费在线观看影片大全网站| 欧美日韩福利视频一区二区| 日韩人妻高清精品专区| 少妇人妻精品综合一区二区 | 女同久久另类99精品国产91| 久久天躁狠狠躁夜夜2o2o| 免费无遮挡裸体视频| 亚洲av熟女| 日本熟妇午夜| 91久久精品国产一区二区成人 | 久久久久久国产a免费观看| 国产探花在线观看一区二区| 女生性感内裤真人,穿戴方法视频| 啦啦啦观看免费观看视频高清| 99精品欧美一区二区三区四区| 欧美大码av| 国产91精品成人一区二区三区| 国产av一区在线观看免费| 久久婷婷人人爽人人干人人爱| 欧洲精品卡2卡3卡4卡5卡区| 黄色女人牲交| 国产 一区 欧美 日韩| 一进一出好大好爽视频| 九色国产91popny在线| 91麻豆精品激情在线观看国产| 天堂影院成人在线观看| 深夜精品福利| 国产欧美日韩精品一区二区| or卡值多少钱| 国产精品亚洲美女久久久| 制服丝袜大香蕉在线| 好看av亚洲va欧美ⅴa在| 国产精品女同一区二区软件 | 啦啦啦观看免费观看视频高清| 中文字幕人妻丝袜一区二区| 怎么达到女性高潮| 成人特级av手机在线观看| 中文字幕av在线有码专区| 淫秽高清视频在线观看| 欧美+日韩+精品| 亚洲七黄色美女视频| 国产毛片a区久久久久| 色播亚洲综合网| 婷婷精品国产亚洲av| 欧美xxxx黑人xx丫x性爽| 又紧又爽又黄一区二区| 久久九九热精品免费| 亚洲一区二区三区色噜噜| tocl精华| 日韩欧美国产一区二区入口| 国产97色在线日韩免费| avwww免费| 国内久久婷婷六月综合欲色啪| 免费av毛片视频| 国产综合懂色| 长腿黑丝高跟| 在线播放无遮挡| 少妇裸体淫交视频免费看高清| 内射极品少妇av片p| 亚洲欧美日韩卡通动漫| h日本视频在线播放| 中亚洲国语对白在线视频| 国产精品亚洲一级av第二区| 黄片大片在线免费观看| 成人av一区二区三区在线看| 国产日本99.免费观看| 亚洲欧美日韩无卡精品| 亚洲精品乱码久久久v下载方式 | 亚洲午夜理论影院| 狠狠狠狠99中文字幕| 成人国产综合亚洲| 高清毛片免费观看视频网站| 91久久精品电影网| 成年免费大片在线观看| 啦啦啦免费观看视频1| 国产久久久一区二区三区| 91麻豆精品激情在线观看国产| 久久久久久久亚洲中文字幕 | 国产精品美女特级片免费视频播放器| av天堂在线播放| 狂野欧美激情性xxxx| 少妇丰满av| 在线播放国产精品三级| 亚洲内射少妇av| 中文字幕av在线有码专区| 免费看美女性在线毛片视频| 久久久久国产精品人妻aⅴ院| 欧美黑人巨大hd| 夜夜看夜夜爽夜夜摸| 一本一本综合久久| 亚洲精品在线观看二区| 成熟少妇高潮喷水视频| 欧美黄色片欧美黄色片| 久久久久久久久中文| 一个人免费在线观看电影| 成人国产综合亚洲| 俺也久久电影网| 一区二区三区高清视频在线| 国产69精品久久久久777片| 欧美bdsm另类| 亚洲av成人精品一区久久| 久久香蕉国产精品| 日韩国内少妇激情av| 亚洲五月天丁香| 999久久久精品免费观看国产| 国产成人aa在线观看| 色在线成人网| 精品久久久久久久久久久久久| 十八禁网站免费在线| 欧美高清成人免费视频www| 18禁美女被吸乳视频| 网址你懂的国产日韩在线| 免费看十八禁软件| 久久精品91蜜桃| 国产美女午夜福利| 日韩欧美精品v在线| 久久伊人香网站| 波多野结衣高清作品| 久久精品国产综合久久久| av视频在线观看入口| 亚洲 欧美 日韩 在线 免费| 操出白浆在线播放| 亚洲国产精品sss在线观看| 成人18禁在线播放| 床上黄色一级片| 国产免费男女视频| 国内揄拍国产精品人妻在线| 老司机福利观看| 丰满的人妻完整版| 美女被艹到高潮喷水动态| 香蕉丝袜av| 麻豆成人av在线观看| 99精品久久久久人妻精品| av中文乱码字幕在线| 亚洲国产中文字幕在线视频| 欧美黄色片欧美黄色片| 在线观看舔阴道视频| 嫩草影院入口| 12—13女人毛片做爰片一| 国产69精品久久久久777片| 99久久精品热视频| 波多野结衣巨乳人妻| 国产高清有码在线观看视频| 国产极品精品免费视频能看的| 97超级碰碰碰精品色视频在线观看| 国内揄拍国产精品人妻在线| 色视频www国产| 久久精品夜夜夜夜夜久久蜜豆| 国产av一区在线观看免费| 国产精品久久久久久精品电影| 日本三级黄在线观看| 亚洲成人精品中文字幕电影| 久久久久性生活片| 国产高清视频在线播放一区| 免费看十八禁软件| 九九在线视频观看精品| 性色avwww在线观看| 亚洲国产高清在线一区二区三| 久久九九热精品免费| 国产极品精品免费视频能看的| 高清毛片免费观看视频网站| 在线观看免费视频日本深夜| 成人无遮挡网站| 久久香蕉精品热| 一边摸一边抽搐一进一小说| 亚洲国产精品合色在线| 在线观看美女被高潮喷水网站 | 两人在一起打扑克的视频| 女人高潮潮喷娇喘18禁视频| 色综合欧美亚洲国产小说| 男女床上黄色一级片免费看| 国产精品99久久久久久久久| 国产亚洲av嫩草精品影院| 人人妻人人澡欧美一区二区| 中文在线观看免费www的网站| 中文字幕熟女人妻在线| 小说图片视频综合网站| 午夜激情福利司机影院| 久久99热这里只有精品18| 真人做人爱边吃奶动态| 欧美黑人巨大hd| 在线观看av片永久免费下载| 国产久久久一区二区三区| 亚洲国产欧美人成| 一个人看的www免费观看视频| 久99久视频精品免费| 国产不卡一卡二| 成人一区二区视频在线观看| 两个人视频免费观看高清| 成人av一区二区三区在线看| 欧美极品一区二区三区四区| 国产精品三级大全| 亚洲av二区三区四区| 狂野欧美白嫩少妇大欣赏| 搡老岳熟女国产| 国产精品 欧美亚洲| 亚洲,欧美精品.| 最近视频中文字幕2019在线8| tocl精华| 亚洲精品在线观看二区| 免费在线观看亚洲国产| 日本与韩国留学比较| 99久国产av精品| 午夜福利欧美成人| 亚洲第一电影网av| 久久久久久久久久黄片| 欧美大码av| 女生性感内裤真人,穿戴方法视频| 免费一级毛片在线播放高清视频| 国产一区二区亚洲精品在线观看| 欧美高清成人免费视频www| 免费看美女性在线毛片视频| 欧美激情久久久久久爽电影| 国产精品久久视频播放| 黄色成人免费大全| 亚洲五月婷婷丁香| 精品一区二区三区人妻视频| 99久久精品国产亚洲精品| 国内揄拍国产精品人妻在线| 国产成人a区在线观看| 免费电影在线观看免费观看| 日日夜夜操网爽| 久久久精品欧美日韩精品| 国产午夜精品论理片| 国产黄色小视频在线观看| 亚洲国产精品合色在线| 夜夜爽天天搞| 麻豆久久精品国产亚洲av| e午夜精品久久久久久久| 村上凉子中文字幕在线| av片东京热男人的天堂| 啦啦啦免费观看视频1| 日本三级黄在线观看| 免费电影在线观看免费观看| 国内精品久久久久精免费| 国产黄a三级三级三级人| 精品久久久久久久人妻蜜臀av| 精品乱码久久久久久99久播| 看黄色毛片网站| 国产精品 国内视频| www日本在线高清视频| 久久久久国内视频| 三级国产精品欧美在线观看| 在线免费观看不下载黄p国产 | 老司机午夜十八禁免费视频| 久久久久国产精品人妻aⅴ院| 丝袜美腿在线中文| 国产午夜福利久久久久久| 亚洲成人久久爱视频| 亚洲成人精品中文字幕电影| 非洲黑人性xxxx精品又粗又长| 日韩 欧美 亚洲 中文字幕| 一级毛片女人18水好多| 精品免费久久久久久久清纯| 国产精品1区2区在线观看.| 久久午夜亚洲精品久久| 久久人人精品亚洲av| 又爽又黄无遮挡网站| 亚洲国产精品999在线| 男女那种视频在线观看| 欧美性感艳星| 国产成+人综合+亚洲专区| 欧美在线一区亚洲|