• <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)綜合題お
    嫩草影视91久久| 国产亚洲精品一区二区www| 亚洲精品中文字幕在线视频| 欧美黄色片欧美黄色片| 午夜福利在线观看吧| 极品教师在线免费播放| 国产极品粉嫩免费观看在线| 一级毛片高清免费大全| 黄网站色视频无遮挡免费观看| 18禁黄网站禁片免费观看直播| 欧美性猛交黑人性爽| 国产av又大| 久久伊人香网站| 亚洲自偷自拍图片 自拍| 亚洲专区中文字幕在线| 一个人免费在线观看的高清视频| 久久久精品国产亚洲av高清涩受| 黑人操中国人逼视频| 国产爱豆传媒在线观看 | 日日干狠狠操夜夜爽| aaaaa片日本免费| 99精品在免费线老司机午夜| 国产成人欧美| 天堂√8在线中文| 99riav亚洲国产免费| 午夜久久久在线观看| 色老头精品视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 日韩欧美免费精品| 国产成人av教育| 国产99久久九九免费精品| 不卡一级毛片| 亚洲久久久国产精品| 国产亚洲av嫩草精品影院| 日本在线视频免费播放| 午夜两性在线视频| 亚洲男人天堂网一区| 久久久久国内视频| 黄色丝袜av网址大全| 精品福利观看| 草草在线视频免费看| 人成视频在线观看免费观看| 国产久久久一区二区三区| 亚洲国产精品合色在线| 亚洲国产看品久久| 久久亚洲真实| 欧美日韩亚洲国产一区二区在线观看| 欧美绝顶高潮抽搐喷水| www日本黄色视频网| 黄频高清免费视频| 狂野欧美激情性xxxx| 又黄又粗又硬又大视频| 男人舔女人的私密视频| 久久久久久国产a免费观看| 最近在线观看免费完整版| 精品国产乱子伦一区二区三区| 中文亚洲av片在线观看爽| 日韩欧美一区二区三区在线观看| 午夜福利视频1000在线观看| 淫妇啪啪啪对白视频| 精品国产超薄肉色丝袜足j| 久久久久久久久久黄片| 99在线视频只有这里精品首页| 国内精品久久久久精免费| 亚洲精品中文字幕在线视频| av有码第一页| 久久这里只有精品19| 色播在线永久视频| 国产精品久久久av美女十八| 两个人视频免费观看高清| 老司机福利观看| 在线十欧美十亚洲十日本专区| 天天躁狠狠躁夜夜躁狠狠躁| 母亲3免费完整高清在线观看| 法律面前人人平等表现在哪些方面| 成人国产综合亚洲| 亚洲成人国产一区在线观看| 国产真实乱freesex| 免费在线观看日本一区| 久久久久久九九精品二区国产 | 久久香蕉激情| 国产1区2区3区精品| 亚洲专区中文字幕在线| 欧美激情高清一区二区三区| 校园春色视频在线观看| 午夜两性在线视频| 亚洲精品美女久久久久99蜜臀| 99在线视频只有这里精品首页| 99精品在免费线老司机午夜| 亚洲成人精品中文字幕电影| 国产成人精品久久二区二区91| www国产在线视频色| 日韩精品中文字幕看吧| 香蕉av资源在线| 日韩大码丰满熟妇| 黄色毛片三级朝国网站| 亚洲国产精品999在线| 一夜夜www| 99热这里只有精品一区 | 午夜成年电影在线免费观看| 国产午夜福利久久久久久| 亚洲精品美女久久久久99蜜臀| 午夜精品久久久久久毛片777| 麻豆一二三区av精品| 人人妻,人人澡人人爽秒播| 听说在线观看完整版免费高清| 国产成人欧美在线观看| 变态另类成人亚洲欧美熟女| 女生性感内裤真人,穿戴方法视频| 中文字幕精品亚洲无线码一区 | 国内揄拍国产精品人妻在线 | 在线观看日韩欧美| 好男人在线观看高清免费视频 | 亚洲精品国产区一区二| 国产av在哪里看| 亚洲熟妇熟女久久| 久久久久久亚洲精品国产蜜桃av| 久热爱精品视频在线9| 欧美黄色片欧美黄色片| 中文字幕精品亚洲无线码一区 | 亚洲男人天堂网一区| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看舔阴道视频| 久久人人精品亚洲av| 国产成人精品久久二区二区91| 日本在线视频免费播放| 亚洲专区国产一区二区| 国产亚洲精品第一综合不卡| av有码第一页| 欧美色视频一区免费| 高清在线国产一区| 在线十欧美十亚洲十日本专区| 亚洲专区中文字幕在线| 中文字幕久久专区| 夜夜看夜夜爽夜夜摸| 村上凉子中文字幕在线| 琪琪午夜伦伦电影理论片6080| 一区二区三区精品91| 中文字幕最新亚洲高清| 老司机在亚洲福利影院| 欧美又色又爽又黄视频| 亚洲男人天堂网一区| 大香蕉久久成人网| 一区二区三区国产精品乱码| 美女免费视频网站| 久久午夜亚洲精品久久| 日本熟妇午夜| 亚洲av成人av| 十分钟在线观看高清视频www| 国产乱人伦免费视频| 免费在线观看完整版高清| 好男人在线观看高清免费视频 | 精品免费久久久久久久清纯| 亚洲在线自拍视频| 久久久久亚洲av毛片大全| 欧美成人性av电影在线观看| 人妻丰满熟妇av一区二区三区| av天堂在线播放| 成人午夜高清在线视频 | 国产精品自产拍在线观看55亚洲| 国产精品永久免费网站| 亚洲一区高清亚洲精品| 手机成人av网站| 亚洲成人国产一区在线观看| 人人妻人人澡人人看| 欧美日韩瑟瑟在线播放| 国产一区二区激情短视频| 免费在线观看成人毛片| 人人妻人人看人人澡| 妹子高潮喷水视频| 久久狼人影院| 90打野战视频偷拍视频| 真人一进一出gif抽搐免费| 国产精品 欧美亚洲| 在线国产一区二区在线| 亚洲 国产 在线| 19禁男女啪啪无遮挡网站| 最近在线观看免费完整版| 香蕉国产在线看| 日韩欧美免费精品| 国产成人一区二区三区免费视频网站| 国产1区2区3区精品| 12—13女人毛片做爰片一| 久久婷婷人人爽人人干人人爱| 51午夜福利影视在线观看| 操出白浆在线播放| 最近最新中文字幕大全电影3 | 日本免费一区二区三区高清不卡| 免费在线观看完整版高清| 在线观看午夜福利视频| 午夜成年电影在线免费观看| 老司机在亚洲福利影院| 成人三级黄色视频| 国产精品爽爽va在线观看网站 | 国产aⅴ精品一区二区三区波| 天天一区二区日本电影三级| 精品人妻1区二区| videosex国产| 久久久国产成人免费| 精品一区二区三区av网在线观看| 久久国产精品影院| 91av网站免费观看| 久久国产亚洲av麻豆专区| 国产又色又爽无遮挡免费看| 中文字幕精品亚洲无线码一区 | 性色av乱码一区二区三区2| 最近最新中文字幕大全电影3 | 午夜免费激情av| 亚洲av成人一区二区三| 亚洲精品色激情综合| 在线av久久热| 老司机靠b影院| 嫩草影院精品99| 午夜激情福利司机影院| 高潮久久久久久久久久久不卡| 国产高清视频在线播放一区| 久久精品成人免费网站| 好看av亚洲va欧美ⅴa在| 宅男免费午夜| 亚洲av电影不卡..在线观看| 99国产极品粉嫩在线观看| www.www免费av| 在线观看舔阴道视频| 宅男免费午夜| 美女扒开内裤让男人捅视频| 亚洲专区中文字幕在线| 禁无遮挡网站| 国产免费男女视频| 亚洲无线在线观看| 国产人伦9x9x在线观看| 99re在线观看精品视频| 久久久精品国产亚洲av高清涩受| 国产精品亚洲一级av第二区| 国产成人精品久久二区二区91| 精品国内亚洲2022精品成人| 1024香蕉在线观看| 精品日产1卡2卡| 午夜福利高清视频| 久久久久国内视频| 女性被躁到高潮视频| 亚洲精品粉嫩美女一区| 91老司机精品| 国产av一区二区精品久久| 国产麻豆成人av免费视频| 麻豆成人av在线观看| 丝袜美腿诱惑在线| 中文字幕av电影在线播放| 一a级毛片在线观看| 人人妻人人澡人人看| 久久国产精品影院| 色播亚洲综合网| 色综合婷婷激情| 国产精品免费一区二区三区在线| 麻豆国产av国片精品| 啦啦啦韩国在线观看视频| 午夜福利在线在线| 麻豆国产av国片精品| 欧美日韩精品网址| 丝袜美腿诱惑在线| 免费在线观看日本一区| 在线免费观看的www视频| 欧美成人免费av一区二区三区| 国产爱豆传媒在线观看 | 色尼玛亚洲综合影院| 亚洲av成人不卡在线观看播放网| 制服诱惑二区| 国产亚洲精品久久久久久毛片| 国内毛片毛片毛片毛片毛片| 一区二区三区国产精品乱码| 久久精品影院6| 欧美激情久久久久久爽电影| 色综合站精品国产| 国产精品久久电影中文字幕| 巨乳人妻的诱惑在线观看| 国产极品粉嫩免费观看在线| 欧美精品亚洲一区二区| 少妇裸体淫交视频免费看高清 | 日韩国内少妇激情av| 久久久久久国产a免费观看| 欧美大码av| 中文字幕人成人乱码亚洲影| 亚洲成a人片在线一区二区| 丝袜在线中文字幕| 中亚洲国语对白在线视频| 中文字幕精品免费在线观看视频| 国产精品自产拍在线观看55亚洲| 一级毛片高清免费大全| 女生性感内裤真人,穿戴方法视频| 级片在线观看| 色综合亚洲欧美另类图片| 99久久综合精品五月天人人| 色在线成人网| 桃色一区二区三区在线观看| 999久久久国产精品视频| 男人的好看免费观看在线视频 | 成人三级黄色视频| 首页视频小说图片口味搜索| 国产男靠女视频免费网站| 亚洲一卡2卡3卡4卡5卡精品中文| 国产99白浆流出| 色精品久久人妻99蜜桃| 国产黄片美女视频| 母亲3免费完整高清在线观看| 国产成+人综合+亚洲专区| 久久久久久国产a免费观看| 国产又黄又爽又无遮挡在线| 亚洲激情在线av| 欧美日韩亚洲国产一区二区在线观看| 男女午夜视频在线观看| 两人在一起打扑克的视频| 久久久久久九九精品二区国产 | 高潮久久久久久久久久久不卡| 少妇被粗大的猛进出69影院| 亚洲国产中文字幕在线视频| 久久狼人影院| 色综合欧美亚洲国产小说| 国产三级在线视频| av片东京热男人的天堂| 嫩草影视91久久| 国产亚洲av嫩草精品影院| 一二三四在线观看免费中文在| 日本免费a在线| 一区二区三区国产精品乱码| 亚洲精品粉嫩美女一区| 一级毛片精品| 女性生殖器流出的白浆| 真人做人爱边吃奶动态| 色综合站精品国产| 午夜日韩欧美国产| 日韩视频一区二区在线观看| 欧美日本视频| 少妇裸体淫交视频免费看高清 | 国产黄片美女视频| 欧美黑人巨大hd| 99国产精品一区二区三区| 久久中文字幕人妻熟女| 国产亚洲精品一区二区www| a级毛片a级免费在线| 日日摸夜夜添夜夜添小说| 波多野结衣高清作品| 中文在线观看免费www的网站 | 中国美女看黄片| 国产熟女xx| 在线永久观看黄色视频| 欧美色欧美亚洲另类二区| 一边摸一边做爽爽视频免费| 听说在线观看完整版免费高清| 国产精品久久电影中文字幕| 黄网站色视频无遮挡免费观看| 国产午夜精品久久久久久| 女人被狂操c到高潮| 女同久久另类99精品国产91| videosex国产| www.熟女人妻精品国产| 97碰自拍视频| 俄罗斯特黄特色一大片| 精品久久久久久久久久久久久 | 99精品欧美一区二区三区四区| 亚洲欧美激情综合另类| 国产精品1区2区在线观看.| 女生性感内裤真人,穿戴方法视频| 国产精品亚洲av一区麻豆| 变态另类成人亚洲欧美熟女| 亚洲av电影在线进入| 亚洲av中文字字幕乱码综合 | 国产亚洲精品一区二区www| 亚洲性夜色夜夜综合| 99久久99久久久精品蜜桃| 午夜福利免费观看在线| 人人妻人人看人人澡| 成人亚洲精品一区在线观看| 男男h啪啪无遮挡| 成人亚洲精品一区在线观看| 亚洲av成人一区二区三| www.999成人在线观看| 国产亚洲欧美精品永久| 999久久久国产精品视频| 男女视频在线观看网站免费 | 热re99久久国产66热| 久久久久久国产a免费观看| 国产精品野战在线观看| 欧美性猛交黑人性爽| 精品久久久久久久人妻蜜臀av| 色综合站精品国产| 国产亚洲精品一区二区www| 久久中文字幕一级| 欧美绝顶高潮抽搐喷水| 人成视频在线观看免费观看| 曰老女人黄片| 又大又爽又粗| 亚洲精品在线美女| 后天国语完整版免费观看| 身体一侧抽搐| bbb黄色大片| 国产精品影院久久| 亚洲成国产人片在线观看| 亚洲国产精品999在线| 老司机福利观看| 日韩精品免费视频一区二区三区| 日本在线视频免费播放| 欧美黑人欧美精品刺激| 久久国产亚洲av麻豆专区| 午夜福利免费观看在线| 欧美绝顶高潮抽搐喷水| 男女视频在线观看网站免费 | 丰满的人妻完整版| 亚洲欧美日韩高清在线视频| 色综合站精品国产| 成人永久免费在线观看视频| 免费搜索国产男女视频| 亚洲成人免费电影在线观看| 麻豆国产av国片精品| 男人舔女人的私密视频| 成人国语在线视频| 久久亚洲精品不卡| 久久人妻av系列| 亚洲精品美女久久久久99蜜臀| 99精品在免费线老司机午夜| 久久久久久大精品| 亚洲va日本ⅴa欧美va伊人久久| 狠狠狠狠99中文字幕| 国产精品99久久99久久久不卡| 又大又爽又粗| 国产高清视频在线播放一区| x7x7x7水蜜桃| 99re在线观看精品视频| 正在播放国产对白刺激| 久久久水蜜桃国产精品网| 热re99久久国产66热| av在线天堂中文字幕| 色婷婷久久久亚洲欧美| 丝袜人妻中文字幕| 日本免费一区二区三区高清不卡| bbb黄色大片| av中文乱码字幕在线| 免费看日本二区| 日韩欧美一区视频在线观看| 最近最新中文字幕大全电影3 | 欧美中文综合在线视频| 99精品欧美一区二区三区四区| 日日夜夜操网爽| 精品久久久久久久毛片微露脸| 欧美绝顶高潮抽搐喷水| 中文字幕高清在线视频| 亚洲国产精品sss在线观看| 久99久视频精品免费| 欧美三级亚洲精品| 男女视频在线观看网站免费 | 麻豆成人av在线观看| 最近最新中文字幕大全免费视频| 亚洲av成人一区二区三| 欧美激情极品国产一区二区三区| 久久国产精品人妻蜜桃| 国产高清激情床上av| 久久久久国产一级毛片高清牌| 午夜免费成人在线视频| 日韩有码中文字幕| 久久婷婷成人综合色麻豆| 久久久久九九精品影院| av免费在线观看网站| 欧美成人性av电影在线观看| 免费在线观看视频国产中文字幕亚洲| 国产精品爽爽va在线观看网站 | 美女大奶头视频| 成人18禁在线播放| АⅤ资源中文在线天堂| 一本综合久久免费| 久久天躁狠狠躁夜夜2o2o| 麻豆久久精品国产亚洲av| 老司机午夜十八禁免费视频| 美女午夜性视频免费| 久久午夜综合久久蜜桃| 久久久久久大精品| 老司机在亚洲福利影院| 久9热在线精品视频| 波多野结衣高清无吗| 欧美成人免费av一区二区三区| 亚洲激情在线av| 久久 成人 亚洲| 制服丝袜大香蕉在线| 久热爱精品视频在线9| 99精品久久久久人妻精品| 精品电影一区二区在线| 香蕉国产在线看| 欧美亚洲日本最大视频资源| 久99久视频精品免费| 色在线成人网| 国产高清videossex| 国产一区二区在线av高清观看| 亚洲av成人一区二区三| 老鸭窝网址在线观看| 久久婷婷成人综合色麻豆| 在线观看免费午夜福利视频| 久久久久久人人人人人| 精品熟女少妇八av免费久了| 美女午夜性视频免费| www.熟女人妻精品国产| 亚洲激情在线av| 正在播放国产对白刺激| 国产精品日韩av在线免费观看| www日本在线高清视频| 大香蕉久久成人网| 国产精品久久久久久亚洲av鲁大| 欧美中文综合在线视频| 国产精品九九99| 天天一区二区日本电影三级| 禁无遮挡网站| 国产三级在线视频| 午夜久久久久精精品| 免费高清视频大片| 亚洲va日本ⅴa欧美va伊人久久| 高清毛片免费观看视频网站| 亚洲真实伦在线观看| www.熟女人妻精品国产| 久久久国产精品麻豆| 老司机在亚洲福利影院| 亚洲一卡2卡3卡4卡5卡精品中文| 97人妻精品一区二区三区麻豆 | 国产精品香港三级国产av潘金莲| 人妻丰满熟妇av一区二区三区| 免费看美女性在线毛片视频| 亚洲熟女毛片儿| 97超级碰碰碰精品色视频在线观看| 久久婷婷人人爽人人干人人爱| 一级毛片精品| 人成视频在线观看免费观看| 久久久久久人人人人人| 精品午夜福利视频在线观看一区| 免费高清视频大片| 757午夜福利合集在线观看| 久久午夜亚洲精品久久| 美女大奶头视频| 国产成人av激情在线播放| 香蕉丝袜av| 国产亚洲精品久久久久久毛片| 免费在线观看成人毛片| 久久天堂一区二区三区四区| 国产亚洲精品综合一区在线观看 | 欧美中文综合在线视频| 久久久久久国产a免费观看| 老司机在亚洲福利影院| 99久久国产精品久久久| 少妇被粗大的猛进出69影院| 午夜日韩欧美国产| 久久久久国内视频| 精品人妻1区二区| 国产人伦9x9x在线观看| 极品教师在线免费播放| 精品国产超薄肉色丝袜足j| 国产精品亚洲av一区麻豆| 黄片小视频在线播放| 国产伦一二天堂av在线观看| 男人的好看免费观看在线视频 | 久久精品aⅴ一区二区三区四区| 高清毛片免费观看视频网站| 欧美成狂野欧美在线观看| 日本免费一区二区三区高清不卡| 亚洲精华国产精华精| 香蕉久久夜色| 国产亚洲精品久久久久5区| 免费在线观看视频国产中文字幕亚洲| 人成视频在线观看免费观看| 黄网站色视频无遮挡免费观看| 嫁个100分男人电影在线观看| 中文亚洲av片在线观看爽| 99riav亚洲国产免费| 亚洲精品一卡2卡三卡4卡5卡| 黄片大片在线免费观看| 国产精品综合久久久久久久免费| 亚洲免费av在线视频| 亚洲精品久久成人aⅴ小说| 色av中文字幕| 亚洲第一青青草原| 国产精品野战在线观看| 日韩欧美免费精品| 国产精品一区二区精品视频观看| 成年人黄色毛片网站| 男人舔女人的私密视频| 日本撒尿小便嘘嘘汇集6| 美女大奶头视频| 欧美另类亚洲清纯唯美| 国产一卡二卡三卡精品| 国产精品野战在线观看| av有码第一页| 91在线观看av| a级毛片在线看网站| 国产男靠女视频免费网站| 亚洲在线自拍视频| 亚洲黑人精品在线| 日本五十路高清| 男女那种视频在线观看| 99国产精品一区二区蜜桃av| 97碰自拍视频| 黑人欧美特级aaaaaa片| 国产成人精品无人区| 日本免费一区二区三区高清不卡| 久久久久精品国产欧美久久久| 精品国产乱码久久久久久男人| 亚洲成人久久性| www.熟女人妻精品国产| 免费女性裸体啪啪无遮挡网站| 成人av一区二区三区在线看| 国产精品永久免费网站| 少妇裸体淫交视频免费看高清 | 少妇 在线观看| 国产单亲对白刺激| av有码第一页| 中文亚洲av片在线观看爽| 老司机福利观看| 精品久久久久久久毛片微露脸| 不卡一级毛片| 久久伊人香网站| 亚洲第一青青草原| 一区二区三区国产精品乱码| 国产黄色小视频在线观看|