999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

硅微陀螺正交誤差直流校正設(shè)計與分析

2014-07-19 10:13:44倪云舫李宏生黃麗斌趙立業(yè)
中國慣性技術(shù)學(xué)報 2014年1期
關(guān)鍵詞:信號檢測

倪云舫,李宏生,楊 波,黃麗斌,趙立業(yè),王 攀,劉 嘉

(東南大學(xué) 儀器科學(xué)與工程學(xué)院 微慣性儀表與先進(jìn)導(dǎo)航技術(shù)教育部重點(diǎn)實(shí)驗室,南京 210096)

硅微陀螺正交誤差直流校正設(shè)計與分析

倪云舫,李宏生,楊 波,黃麗斌,趙立業(yè),王 攀,劉 嘉

(東南大學(xué) 儀器科學(xué)與工程學(xué)院 微慣性儀表與先進(jìn)導(dǎo)航技術(shù)教育部重點(diǎn)實(shí)驗室,南京 210096)

為硅微陀螺設(shè)計了一種正交誤差直流校正方法,通過設(shè)計校正結(jié)構(gòu)和加載直流電壓實(shí)現(xiàn)正交校正。建立了存在正交誤差時敏感質(zhì)量的運(yùn)動微分方程,分析明確了正交響應(yīng)的成因和對驅(qū)動運(yùn)動軌跡的影響。為硅微陀螺設(shè)計了正交校正結(jié)構(gòu),實(shí)現(xiàn)構(gòu)建靜電耦合彈性系數(shù)并利用驅(qū)動運(yùn)動產(chǎn)生校正靜電力。校正力的頻率和相位無需電路控制,使該方法相較于傳統(tǒng)校正方法具備特殊優(yōu)勢。設(shè)計并分析了校正電壓的兩種加載方式,通過實(shí)驗測試驗證了正交響應(yīng)幅值隨校正電壓的變化規(guī)律。實(shí)測校正電壓接近理論值,證明了校正結(jié)構(gòu)設(shè)計的正確性。該方法在雙線振動式硅微陀螺中具有重要應(yīng)用價值。

硅微陀螺;正交誤差;正交校正;直流電壓;結(jié)構(gòu)設(shè)計

硅微陀螺基于硅微機(jī)械加工工藝制造,因其相較于傳統(tǒng)陀螺體積小、功耗低、可量產(chǎn)等優(yōu)勢,在商業(yè)和軍事領(lǐng)域均有廣泛應(yīng)用[1-3]。硅微陀螺為哥氏振動陀螺的一種,利用哥氏效應(yīng)使驅(qū)動振動耦合到檢測方向進(jìn)行角速度檢測。然而,由于加工誤差等非理想因素的存在,驅(qū)動方向和檢測方向的運(yùn)動之間還存在其他耦合機(jī)制,如彈性耦合、阻尼耦合、靜電耦合等,使得在沒有角速度輸入的情況下,硅微陀螺的檢測方向仍有振動信號輸出,直接影響其零偏穩(wěn)定性和零偏溫度穩(wěn)定性等關(guān)鍵性能。由于硅微陀螺的彈性梁對加工誤差十分敏感,彈性耦合一般最為嚴(yán)重,其造成的檢測振動與哥氏效應(yīng)引起的檢測振動在相位上相差90°,故被稱為正交誤差[4]。正交誤差通常遠(yuǎn)大于有用信號,其校正為硅微陀螺研制中的關(guān)鍵技術(shù)之一[5]。目前常用的校正方法有同步解調(diào)法和閉環(huán)抑制法兩種[6]。前者利用與有用信號同頻同相的電壓信號對輸出信號進(jìn)行解調(diào),從而消除輸出中與有用信號相位相差 90°的正交信號。后者通過設(shè)計閉環(huán)反饋控制回路,產(chǎn)生與引起正交信號的彈性耦合力同頻反相的反饋力抑制正交運(yùn)動,從而從源頭上消除正交信號。采用同步解調(diào)法時,正交信號仍然存在,其相對于有用信號而言很大的幅值增加了對前端電路動態(tài)范圍的要求,限制了陀螺的分辨率和信噪比。閉環(huán)抑制法從源頭上消除了正交運(yùn)動,對前端電路的限制放寬,但與同步解調(diào)法需要精確的解調(diào)電壓相位控制類似,閉環(huán)抑制法需要精確的反饋力相位控制。兩種方法都對電路的相位精度提出了很高的要求,尤其當(dāng)環(huán)境變化造成陀螺結(jié)構(gòu)參數(shù)和電路器件參數(shù)變化時更為困難。

針對以上問題,本文在硅微陀螺設(shè)計中采用了一種正交誤差直流校正方法[7-8],通過設(shè)計專用的校正結(jié)構(gòu)產(chǎn)生校正靜電力抵消彈性耦合力,從而抑制正交運(yùn)動。校正力的幅值由加載的直流電壓控制,頻率和相位信息由驅(qū)動振動直接提供,無需專門的控制電路。該法對前端電路動態(tài)范圍和電路相位精度均無特殊要求,具備實(shí)用價值。

1 正交誤差對敏感質(zhì)量運(yùn)動的影響分析

存在彈性耦合時,硅微陀螺敏感質(zhì)量的二維運(yùn)動微分方程可表示為:

式中,mc為敏感質(zhì)量,cx、cy、kx、ky分別為驅(qū)動方向(X向)和檢測方向(Y向)的阻尼、彈性系數(shù),kxy、kyx為驅(qū)動和檢測方向之間的耦合彈性系數(shù);Fd為簡諧靜電驅(qū)動力,幅值和角頻率分別為F0和ωd;Fc為哥氏力,Ωz為輸入角速度(Z向)。

一般而言,模態(tài)不匹配時,檢測振幅遠(yuǎn)小于驅(qū)動振幅,故kxy的影響可忽略不計。彈性耦合對敏感質(zhì)量運(yùn)動的影響主要表現(xiàn)為kyx與驅(qū)動位移x結(jié)合構(gòu)成作用于檢測方向的彈性耦合力Fq,其頻率和相位與x一致,故與哥氏力Fc同頻且相位相差90°,定義為正交力。

硅微陀螺工作時,敏感質(zhì)量沿驅(qū)動方向恒幅諧振,定義振幅為Ad,則x、Fq、Fc可表示為:

定義正交力Fq、哥氏力Fc引起的敏感質(zhì)量沿檢測方向的位移響應(yīng)分別為正交響應(yīng)yq和哥氏響應(yīng)yc,有:

式中,ωy、Qy分別為檢測模態(tài)的諧振頻率和品質(zhì)因數(shù)。可見,yq和yc相較于Fq和Fc存在同樣的相移φ,仍為同頻且相位相差90°。正交誤差等效角速度Ωq可表示為:

假設(shè)驅(qū)動頻率為4 kHz,Qy在幾千量級,Ωz=0 (°)/s,存在Ωq=100 (°)/s時,敏感質(zhì)量質(zhì)心的運(yùn)動軌跡如圖1所示。

圖1 存在正交誤差時敏感質(zhì)量質(zhì)心的運(yùn)動軌跡Fig.1 Motion trajectory of the center of proof mass when quadrature error exists

可見,正交誤差的存在使得敏感質(zhì)量的驅(qū)動運(yùn)動偏離理想的X向,并且由于yq和x之間相差φ的存在,呈現(xiàn)橢圓形運(yùn)動軌跡。耦合彈性系數(shù)kyx不變時,一定的驅(qū)動頻率條件下,檢測與驅(qū)動模態(tài)的頻差Δf越小,正交響應(yīng)yq的幅值越大,yq和x之間的相差φ也越大。yq的幅值由Y向檢測電極提取后,成為陀螺輸出中的正交信號??梢娬恍盘柵c有用信號類似,是一個動態(tài)信號,其幅值和相位同樣受到固有頻率、Q值等結(jié)構(gòu)參數(shù)的影響,因而在環(huán)境變化導(dǎo)致結(jié)構(gòu)參數(shù)波動時,正交信號也將發(fā)生變化,從而使相位控制成為應(yīng)用同步解調(diào)法和正交抑制法進(jìn)行正交校正時的一個難點(diǎn)。正交誤差直流校正方法可以緩解這一問題。

2 正交誤差直流校正結(jié)構(gòu)設(shè)計與分析

正交誤差直流校正的原理為,通過設(shè)計校正結(jié)構(gòu),構(gòu)建如下作用于敏感質(zhì)量的Y向校正靜電力:

式中,keyx為靜電耦合彈性系數(shù),其符號和大小由校正結(jié)構(gòu)與所加直流電壓共同決定。驅(qū)動位移x為機(jī)械信號,由敏感質(zhì)量的驅(qū)動運(yùn)動直接提供給校正結(jié)構(gòu),因此無需進(jìn)行頻率和相位控制。比較式(3)(6)可見,校正力Fey與正交力Fq的形式一致,當(dāng)構(gòu)建的靜電耦合彈性系數(shù)keyx與結(jié)構(gòu)耦合彈性系數(shù)kyx等值反號時,F(xiàn)ey實(shí)現(xiàn)對Fq的平衡,從而抑制正交運(yùn)動,實(shí)現(xiàn)正交校正。

本文在自主研制的某型雙質(zhì)量硅微陀螺[9]上設(shè)計了校正結(jié)構(gòu),并完成了第一回合的制版加工。包含校正結(jié)構(gòu)的硅微陀螺結(jié)構(gòu)示意圖及校正結(jié)構(gòu)的局部電鏡照片如圖2所示。

圖2 雙質(zhì)量硅微陀螺結(jié)構(gòu)示意圖及校正結(jié)構(gòu)局部照片F(xiàn)ig.2 Structure of the dual-mass silicon micro-gyroscope and a local SEM photo of the quadrature correction structure

該陀螺結(jié)構(gòu)為全解耦形式,包含驅(qū)動框架、敏感質(zhì)量、檢測框架三部分質(zhì)量,驅(qū)動框架和檢測框架分別只具備X向、Y向自由度,敏感質(zhì)量具備X、Y兩向自由度,實(shí)現(xiàn)利用哥氏效應(yīng)檢測Z向角速度。

校正結(jié)構(gòu)主要由布置于敏感質(zhì)量中軸線兩側(cè)的不等間隙梳齒構(gòu)成,初始的大、小間隙分別為pd0和d0??蓜邮猃X位于敏感質(zhì)量上,左右梳齒對稱布置。固定梳齒連接于兩組校正電極A和B上,分別施加有直流電壓VA和VB,左右梳齒根據(jù)大、小間隙的Y向位置反對稱布置。單組校正結(jié)構(gòu)由4個可動梳齒和相應(yīng)的固定梳齒構(gòu)成,此種組合方式的設(shè)計目的在于消除伴隨校正靜電力同時出現(xiàn)的附加干擾靜電力。敏感質(zhì)量的Y向位移很小(y<

式中,F(xiàn)ex、Fey分別為X、Y向靜電力,kexx、keyy分別為X、Y向靜電彈性系數(shù),kexy、keyx為兩個方向之間的靜電耦合彈性系數(shù)。

校正結(jié)構(gòu)對敏感質(zhì)量運(yùn)動的影響主要表現(xiàn)為構(gòu)建靜電耦合彈性系數(shù)keyx,結(jié)合位于敏感質(zhì)量上的可動梳齒的驅(qū)動運(yùn)動x,形成作用于Y向的校正靜電力Fey,用于平衡引起正交誤差的正交力Fq。keyx的大小由校正結(jié)構(gòu)參數(shù)n、h、d0、p和所加直流電壓VA、VB的平方差決定,正、負(fù)符號由VA、VB的相對大小決定。結(jié)合式(3)(5)(7)可知,校正等效角速度為Ωq的正交誤差所需的直流電壓加載條件為:

另外,由式(8)可見,校正結(jié)構(gòu)同時也產(chǎn)生Y向靜電彈性系數(shù)keyy,影響檢測模態(tài)的諧振頻率。根據(jù)目的不同,可對keyy的大小進(jìn)行相應(yīng)設(shè)計上的調(diào)節(jié)。對照keyx和keyy的表達(dá)式可見:若希望正交校正對檢測模態(tài)諧振頻率的影響最小化,可在設(shè)計中令校正結(jié)構(gòu)梳齒的交疊長度l0盡可能小,并在結(jié)構(gòu)模態(tài)設(shè)計中將正交校正引起的頻差變化納入考量;若希望利用keyy進(jìn)行檢測頻率調(diào)節(jié),則可適當(dāng)增大l0的設(shè)計值,并對直流電壓VA、VB的加載方式進(jìn)行相應(yīng)的優(yōu)化設(shè)計。

3 校正電壓加載方式設(shè)計與分析

正交校正過程中,校正電壓的調(diào)節(jié)依賴于正交響應(yīng)yq與所加直流電壓VA、VB的相對變化規(guī)律。校正結(jié)構(gòu)對敏感質(zhì)量檢測方向運(yùn)動微分方程的影響主要表現(xiàn)為引入靜電彈性系數(shù)keyx和keyy,即有:

求解該式可得此時的正交響應(yīng)yq為:

對于圖2所示的校正結(jié)構(gòu),若記:

則靜電彈性系數(shù)keyx和keyy可表示為:

校正電壓的加載方式主要有兩種:(a)只在單個電極上施加電壓,根據(jù)正交信號的相位,判斷校正所需keyx的正、負(fù)符號,在校正電極A、B中選擇一個施加直流電壓,調(diào)節(jié)該電壓使正交信號最小化。該加載方式實(shí)現(xiàn)簡單,但需要手動選擇校正電極,在模擬電路中僅適用于開環(huán)校正。(b)在兩個電極上施加電壓,分別為VA、VB,通過調(diào)節(jié)其平方差使正交信號最小化,一般將二者取為VA=VD+Vq,VB=VD-Vq,有:

式中,VD為預(yù)置直流電壓,Vq為調(diào)節(jié)電壓。VD一定時,調(diào)節(jié)Vq既可控制keyx的大小也可改變其正、負(fù)符號,因而該加載方式適用于對正交幅值進(jìn)行閉環(huán)控制。

選擇加載方式(a)時,假設(shè)校正所需的keyx為負(fù),正交響應(yīng)幅值|yq|與VA的關(guān)系為:

選擇加載方式(b)時,取VD為固定值,正交響應(yīng)幅值|yq|與Vq的關(guān)系為:

可見,采用加載方式(a)時,|yq|與VA為非線性關(guān)系,且在調(diào)節(jié)VA的過程中,檢測模態(tài)諧振頻率受到直接影響,變化較大。采用加載方式(b)時,若取VD>>Vq,則調(diào)節(jié)Vq對檢測諧振頻率的影響可以降至很小,且|yq|與Vq為近似線性關(guān)系。若希望正交校正對檢測諧振頻率的影響最小化,加載方式(b)具備優(yōu)勢。

無論哪種加載方式,使正交響應(yīng)幅值為零的電壓值為所需的正交校正電壓,電壓小于該值或大于該值時正交響應(yīng)均不為零,處于欠校正或過校正狀態(tài)。

4 實(shí)驗測試

加工得到帶正交校正結(jié)構(gòu)的雙質(zhì)量硅微陀螺芯片,以陶瓷管殼進(jìn)行真空封裝。安裝在PCB電路板上的硅微陀螺表頭如圖3所示。

圖3 硅微陀螺表頭及測試電路Fig.3 Vacuum packaged gyroscope on PCB with test circuits

圖4 無角速度輸入時的驅(qū)動速度信號與輸出信號Fig.4 Waveform comparison between the drive velocity signal and the zero-rate output one

選擇正交誤差較大的某硅微陀螺表頭進(jìn)行正交校正實(shí)驗。無角速度輸入條件下,令驅(qū)動回路正常工作,表頭輸出信號與驅(qū)動速度信號的波形比較如圖4所示,圖中上方波形為驅(qū)動速度信號,下方波形為輸出信號,可見輸出信號在相位上滯后驅(qū)動速度信號大約90°,主要為正交信號,且結(jié)構(gòu)耦合彈性系數(shù)kyx為正。由輸出信號幅值結(jié)合相關(guān)電路參數(shù)折算可得該表頭的正交誤差等效角速度約為354 (°)/s。采用兩種校正電壓加載方式對該表頭進(jìn)行正交校正,記錄輸出信號幅值隨校正電壓的變化情況。同時,由式(15)(16)計算正交響應(yīng)幅值與校正電壓相對變化關(guān)系的理論曲線。選擇加載方式(a)(b)時,理論曲線與實(shí)測數(shù)據(jù)的比較分別如圖5、圖6所示。

圖5 在單個電極上施加電壓時正交信號隨VA的變化Fig.5 Relative relationship between quadrature signal and correction voltage when only one electrode is in use

可見,兩種電壓加載方式下,實(shí)測數(shù)據(jù)的變化趨勢與理論曲線基本一致,證明了理論分析的正確性。設(shè)計的校正結(jié)構(gòu)參數(shù)為:n= 45,d0= 4,p= 2.4。敏感質(zhì)量mc≈ 0.5 mg,驅(qū)動頻率為3.8 kHz。由式(9)計算可得,理論上校正354 (°)/s正交誤差所需的VA、VB平方差約為119.5 V2。在單個電極上加載時,所需理論校正電壓VA約為10.9 V,實(shí)測值為12 V,略大于理論值。在兩個電極上加載時,若取VD為10 V,則所需理論校正電壓Vq約為2.9 V,實(shí)測值為3 V,與理論值吻合良好。

實(shí)際測試中,調(diào)節(jié)校正電壓使輸出信號幅值最小化時,仍然存在100 mV左右的殘余信號,將其與驅(qū)動速度信號對比可以發(fā)現(xiàn),殘余信號為與有用信號同頻同相的偏移誤差部分,正交信號已被完全消除。

圖6 在兩個電極上施加電壓時正交信號隨Vq的變化(VD= 10 V)Fig.6 Relative relationship between quadrature signal and correction voltage when both electrodes are in use(VD= 10 V)

5 結(jié) 論

本文設(shè)計并分析了硅微陀螺的一種正交誤差直流校正方法。對敏感質(zhì)量的運(yùn)動分析表明,正交誤差主要由結(jié)構(gòu)耦合彈性系數(shù)結(jié)合驅(qū)動運(yùn)動形成的正交力作用于檢測方向引起正交響應(yīng)產(chǎn)生。合理設(shè)計正交校正結(jié)構(gòu),可實(shí)現(xiàn)構(gòu)建靜電耦合彈性系數(shù),并利用敏感質(zhì)量的驅(qū)動運(yùn)動產(chǎn)生校正靜電力,用以平衡正交力。校正力的頻率和相位與正交力自然匹配,無需電路控制。校正電壓可由結(jié)構(gòu)參數(shù)和正交誤差等效角速度進(jìn)行估計。相較于單電極加載方式,校正電壓的雙電極加載方式適用于對正交幅值進(jìn)行閉環(huán)控制,并且可使校正電壓調(diào)節(jié)對檢測諧振頻率的影響最小化。實(shí)驗測試驗證了正交信號隨校正電壓的變化規(guī)律,實(shí)際校正電壓接近理論值,證明了校正結(jié)構(gòu)設(shè)計的正確性。該方法在雙線振動式硅微陀螺中具備通用性和實(shí)用價值。

(References):

[1]Shkel A M,Horowitz R,Seshia A A,et al.Dynamics and control of micromachined gyroscopes[C]// Proceedings of the American Control Conference.San Diego,USA,1999:2119-2124.

[2]Geiger W,Bartholomeyczik J,Breng U,et al.MEMS IMU for AHRS applications[C]//IEEE/ION Position Location and Navigation Symposium.Monterey,USA,2008:225-231.

[3]Lapadatu D,Blixhavn B,Holm R,et al.SAR500 - A high-precision high-stability butterfly gyroscope with north seeking capability[C]// IEEE/ION Position Location and Navigation Symposium.Indian Wells,USA,2010:6-13.

[4]YANG Bo,WANG Shou-rong,LI Hong-sheng,et al.The coupling error analysis of the decoupled silicon micro-gyroscope[C]// 5th IEEE International Conference on Nano/Micro Engineered and Molecular Systems.Xiamen,China,2010:356-361.

[5]羅兵,張輝,吳美平.硅微陀螺正交誤差及其對信號檢測的影響[J].中國慣性技術(shù)學(xué)報,2009,17(5):604-607, 613.LUO Bing,ZHANG Hui,WU Mei-ping.Quadrature signal of microgyroscope and its effect on signal detection [J].Journal of Chinese Inertial Technology,2009,17(5):604-607,613.

[6]Saukoski M,Aaltonen L,Halonen K A I.Zero-rate output and quadrature compensation in vibratory MEMS gyroscopes [J].IEEE Sensors Journal,2007,7(12):1639-1652.

[7]Sharma A,Zaman M F,Ayazi F.A sub-0.2°/hr bias drift micromechanical silicon gyroscope with automatic CMOS mode-matching[J].IEEE Journal of Solid-State Circuits,2009,44(5):1593-1608.

[8]Tatar E,Alper S E,Akin T.Effect of quadrature error on the performance of a fully-decoupled MEMS gyroscope [C]//IEEE 24th International Conference on Micro Electro Mechanical Systems.Cancun,Mexico,2011:569-572.

[9]殷勇,王壽榮,王存超,等.一種雙質(zhì)量硅微陀螺儀[J].中國慣性技術(shù)學(xué)報,2008,16(6):703-706,711.YIN Yong,WANG Shou-rong,WANG Cun-chao,et al.Dual-mass silicon micro-gyroscope[J].Journal of Chinese Inertial Technology,2008,16(6):703-706,711.

Design and analysis of a quadrature correction method using DC voltages for silicon micro-gyroscope

NI Yun-fang,LI Hong-sheng,YANG Bo,HUANG Li-bin,ZHAO Li-ye,WANG Pan,LIU Jia
(Key Laboratory of Micro-Inertial Instrument and Advanced Navigation Technology of Ministry of Education,School of Instrument Science and Engineering,Southeast University,Nanjing 210096,China)

A quadrature correction method using DC voltages for silicon micro-gyroscope was designed.The quadrature correction was accomplished by a specially designed correction structure in combination with DC voltages.The differential equations of the proof mass’s motion were established in consideration of quadrature error.The causes of the quadrature response and how it affects the drive motion trajectory were then derived.The quadrature correction structure was designed for silicon micro-gyroscope,through which a correction electrostatic force can be generated by using the created electrostatic coupling elastic coefficient and the drive motion.Comparied with traditional quadrature correction methods,this method has the advantage of eliminating the need of frequency and phase control.Two voltage loading methods were designed and analyzed,and the relative relationship between the quadrature response and the correction voltage was verified through experimental test.The actual values of correction voltage are close to the theoretical ones,proving the correctness of the structure design.This method has important application value in linear vibration silicon micro-gyroscope.

silicon micro-gyroscope;quadrature error;quadrature correction;DC voltage;structure design

TH824+.3

:A

1005-6734(2014)01-0104-05

10.13695/j.cnki.12-1222/o3.2014.01.021

2013-08-06;

:2013-11-10

國家自然科學(xué)基金項目(61101021);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項資金資助(3222003102)

倪云舫(1986—),女,博士研究生,從事MEMS微慣性儀表研究。E-mail:niyunfang@126.com

聯(lián) 系 人:李宏生(1964—),男,教授,博士生導(dǎo)師。E-mail:hsli@seu.edu.cn

猜你喜歡
信號檢測
“不等式”檢測題
“一元一次不等式”檢測題
“一元一次不等式組”檢測題
信號
鴨綠江(2021年35期)2021-04-19 12:24:18
“幾何圖形”檢測題
“角”檢測題
完形填空二則
孩子停止長個的信號
小波變換在PCB缺陷檢測中的應(yīng)用
基于LabVIEW的力加載信號采集與PID控制
主站蜘蛛池模板: 伊人查蕉在线观看国产精品| 欧美日本激情| 国产一级无码不卡视频| 午夜精品区| 青青青草国产| 无遮挡一级毛片呦女视频| 在线色国产| 亚洲综合经典在线一区二区| jizz亚洲高清在线观看| 在线综合亚洲欧美网站| 亚洲色图欧美激情| 2021亚洲精品不卡a| 国产黄色爱视频| 国产成人亚洲欧美激情| 国产欧美中文字幕| 国产SUV精品一区二区6| 波多野吉衣一区二区三区av| 国内丰满少妇猛烈精品播| 亚洲AV无码久久天堂| 国产精品毛片一区视频播| 国产小视频免费观看| 久久99国产精品成人欧美| 国产青榴视频在线观看网站| 亚洲国产中文精品va在线播放| 国产精品漂亮美女在线观看| 亚洲精品图区| 欧美日韩在线第一页| AV不卡在线永久免费观看| 久久精品国产亚洲AV忘忧草18| 亚洲欧美国产五月天综合| 国产无码制服丝袜| 日韩午夜伦| 久久精品人妻中文视频| 人妻少妇乱子伦精品无码专区毛片| 丰满人妻中出白浆| 美女免费黄网站| 欧美性久久久久| 国产真实二区一区在线亚洲| 草草影院国产第一页| 小说区 亚洲 自拍 另类| 国产精品视频第一专区| 在线观看亚洲精品福利片| 国产欧美在线| 久久性妇女精品免费| 亚洲欧美色中文字幕| 老色鬼欧美精品| 精品国产自| 91娇喘视频| 激情综合网址| 伊人色在线视频| 亚洲欧美日韩久久精品| 欧美 国产 人人视频| 国产乱子伦一区二区=| 国产自在线播放| 久久精品国产精品青草app| 国产微拍一区| 中文字幕欧美日韩高清| 亚洲高清中文字幕在线看不卡| 国产一区二区在线视频观看| 国产第一页免费浮力影院| 无码高潮喷水专区久久| 日本少妇又色又爽又高潮| 日韩色图区| 天堂岛国av无码免费无禁网站| 精品视频91| 99爱视频精品免视看| 无码视频国产精品一区二区| 白浆免费视频国产精品视频| 久久亚洲美女精品国产精品| 日本久久网站| 91福利在线观看视频| 黄色网页在线播放| 久操中文在线| 亚洲中文字幕在线一区播放| 欧美日韩v| 亚洲永久精品ww47国产| h视频在线观看网站| 奇米精品一区二区三区在线观看| 国产小视频在线高清播放| 国产精品乱偷免费视频| 免费va国产在线观看| 六月婷婷激情综合|