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

    光學(xué)元件撓性支撐結(jié)構(gòu)廣義建模及優(yōu)化設(shè)計(jì)

    2016-12-19 02:31:07曹玉巖王志臣韓西達(dá)張耀祖
    光學(xué)精密工程 2016年11期
    關(guān)鍵詞:面形撓性透鏡

    曹玉巖,王志臣,周 超,范 磊,韓西達(dá),張耀祖

    (中國(guó)科學(xué)院 長(zhǎng)春光學(xué)精密機(jī)械與物理研究所,吉林 長(zhǎng)春 130033)

    ?

    光學(xué)元件撓性支撐結(jié)構(gòu)廣義建模及優(yōu)化設(shè)計(jì)

    曹玉巖*,王志臣,周 超,范 磊,韓西達(dá),張耀祖

    (中國(guó)科學(xué)院 長(zhǎng)春光學(xué)精密機(jī)械與物理研究所,吉林 長(zhǎng)春 130033)

    建立了撓性支撐結(jié)構(gòu)的力學(xué)模型及優(yōu)化設(shè)計(jì)模型,以使光學(xué)元件撓性支撐結(jié)構(gòu)同時(shí)滿(mǎn)足元件定位的剛度要求和溫度適應(yīng)性的柔度要求,同時(shí)給出了相應(yīng)的建模方法??紤]撓性支撐結(jié)構(gòu)是由圓周對(duì)稱(chēng)分布的撓性單元組成的,故將撓性單元簡(jiǎn)化為超靜定梁結(jié)構(gòu),應(yīng)用虛功原理推導(dǎo)了撓性單元的徑向及切向剛度。然后,假設(shè)光學(xué)元件為剛性體,根據(jù)力平衡條件及變形協(xié)調(diào)條件,推導(dǎo)了撓性支撐結(jié)構(gòu)的整體剛度,并引入修正因子補(bǔ)償了剛體假設(shè)帶來(lái)的誤差。最后,以撓性支撐結(jié)構(gòu)總變形能為目標(biāo)函數(shù),推導(dǎo)了同時(shí)考慮撓性支撐結(jié)構(gòu)幾何構(gòu)形及參數(shù)的協(xié)同優(yōu)化設(shè)計(jì)模型,通過(guò)引入了整型變量將結(jié)構(gòu)整體剛度簡(jiǎn)化為整型變量和離散剛度的線(xiàn)性組合,從而使優(yōu)化模型中不含有諧波函數(shù)項(xiàng)。基于數(shù)值仿真和實(shí)驗(yàn)對(duì)結(jié)構(gòu)剛度模型進(jìn)行了驗(yàn)證,結(jié)果顯示:實(shí)驗(yàn)、仿真和理論計(jì)算結(jié)果一致。此外,以透鏡支撐為例,驗(yàn)證了撓性支撐結(jié)構(gòu)的優(yōu)化設(shè)計(jì)方法,有限元分析結(jié)果表明,透鏡面形精度較初始設(shè)計(jì)提高了23%。

    光學(xué)元件;撓性支撐結(jié)構(gòu);廣義模型;超靜定;變形能;優(yōu)化設(shè)計(jì)

    1 引 言

    撓性結(jié)構(gòu)已廣泛應(yīng)用于多個(gè)領(lǐng)域,如光學(xué)設(shè)備[1],陀螺儀[2],傳感器及驅(qū)動(dòng)器[3],微位移定位平臺(tái)[4]。撓性結(jié)構(gòu)得以廣泛應(yīng)用主要由于其具有以下優(yōu)點(diǎn):(1)具有很高的重復(fù)精度;(2)低摩擦性;(3)容易加工制造且成本相對(duì)較低。

    光學(xué)元件的面形精度決定了光學(xué)系統(tǒng)性能,而光學(xué)元件的安裝和定位結(jié)構(gòu)又是保證面形精度的關(guān)鍵因素。在光學(xué)設(shè)備中,撓性結(jié)構(gòu)主要用于光學(xué)元件的支撐與定位結(jié)構(gòu)[5-7]。光學(xué)元件,尤其是大口徑的元件,需要非常復(fù)雜的撓性支撐定位結(jié)構(gòu)來(lái)保持元件的位置及面形精度[8-9]。這是因?yàn)閾闲越Y(jié)構(gòu)在某些自由度上具有一定的柔度,容許該自由度上有少量可控的相對(duì)位移,在其它方向上則具有較大剛度。在環(huán)境溫度變化時(shí),光學(xué)元件及其支撐組件將產(chǎn)生變形,熱變形系數(shù)的不匹配將導(dǎo)致出現(xiàn)接觸應(yīng)力,接觸應(yīng)力將直接影響光學(xué)元件的面形,甚至造成破壞,而撓性結(jié)構(gòu)的柔度可以有效地減小這部分接觸應(yīng)力[10-13]。

    盡管撓性結(jié)構(gòu)具有以上優(yōu)點(diǎn),但撓性支撐結(jié)構(gòu)的建模非常復(fù)雜,尤其是多個(gè)撓性元件同時(shí)作用時(shí)。趙磊等[10]設(shè)計(jì)了光刻投影物鏡中透鏡的支撐結(jié)構(gòu),軸向采用的是多點(diǎn)撓性支撐結(jié)構(gòu),支撐中的彈性片簡(jiǎn)化為懸臂梁模型,但未深入分析彈性片數(shù)目對(duì)面形精度的影響。李宇軒等[14]針對(duì)空間遙感器大口徑主反射鏡提出了Cartwheel型雙軸柔鉸的三點(diǎn)撓性支撐結(jié)構(gòu),柔鉸的參數(shù)采用無(wú)量綱設(shè)計(jì)方法確定。盡管通過(guò)此方法確定的參數(shù)可以滿(mǎn)足設(shè)計(jì)要求,但無(wú)法確定此參數(shù)是否為最佳參數(shù)。曹乃亮等[15]提出將NiTi記憶合金絲復(fù)合于柔性件周邊,借助預(yù)應(yīng)變產(chǎn)生的拉應(yīng)力來(lái)降低柔性件根部的應(yīng)力集中問(wèn)題,將柔性件切口等效為彎扭組合的梁結(jié)構(gòu),根據(jù)NiTi合金絲的本構(gòu)方程設(shè)計(jì)了兩種復(fù)合方案,采用有限元軟件進(jìn)行分析。

    撓性支撐結(jié)構(gòu)功能包括兩個(gè)方面:(1)定位,即要求支撐結(jié)構(gòu)具有一定的剛度,使得光學(xué)元件處在不同姿態(tài)時(shí),能夠保持其光軸位置不變;(2)協(xié)調(diào)變形,即當(dāng)存在溫度梯度時(shí),支撐結(jié)構(gòu)要有一定撓性來(lái)協(xié)調(diào)因熱變形系數(shù)不匹配帶來(lái)的變形偏差。撓性結(jié)構(gòu)的這兩方面功能是矛盾的,因?yàn)槿绻3治恢貌蛔?,撓性件剛度?yīng)盡可能大,若為協(xié)調(diào)變形,則撓性件柔度應(yīng)盡可能大。為平衡上述矛盾,需要根據(jù)撓性支撐結(jié)構(gòu)的數(shù)學(xué)模型及使用要求,對(duì)各個(gè)參數(shù)進(jìn)行優(yōu)化,進(jìn)而得到合理方案。為此,本文提出了一種撓性支撐結(jié)構(gòu)的廣義模型,該模型適合于多種類(lèi)型的撓性支撐結(jié)構(gòu)。在力學(xué)模型的基礎(chǔ)上,提出了一種用于同時(shí)優(yōu)化撓性支撐結(jié)構(gòu)幾何構(gòu)形和參數(shù)的模型。模型中通過(guò)引入整型變量將撓性支撐結(jié)構(gòu)整體剛度表達(dá)為整型變量與離散剛度的線(xiàn)性組合,從而簡(jiǎn)化模型優(yōu)化的復(fù)雜程度。從數(shù)值仿真和實(shí)驗(yàn)角度對(duì)整體剛度模型進(jìn)行驗(yàn)證,此外,以420 mm口徑的透鏡為例,進(jìn)行了撓性支撐結(jié)構(gòu)的優(yōu)化設(shè)計(jì),并應(yīng)用有限元軟件分析了優(yōu)化設(shè)計(jì)前后透鏡表面面形精度。

    2 撓性支撐結(jié)構(gòu)力學(xué)建模

    撓性支撐結(jié)構(gòu)如圖1所示。由圖1可知,光學(xué)元件的邊緣被粘接在多個(gè)相同的均勻分布的撓性單元上。相比傳統(tǒng)剛性支撐結(jié)構(gòu),撓性支撐結(jié)構(gòu)具有很多優(yōu)點(diǎn):(1)在受到?jīng)_擊、振動(dòng)及溫度應(yīng)力等影響后,光學(xué)元件的位置及共軸精度不會(huì)受到影響,因?yàn)閾闲栽哂凶冃慰苫謴?fù)的特點(diǎn);(2)當(dāng)溫度變化導(dǎo)致光學(xué)元件與其它結(jié)構(gòu)變形不一致時(shí),撓性結(jié)構(gòu)可明顯減小溫度應(yīng)力,因?yàn)閾闲詷?gòu)件在光學(xué)元件熱變形方向上具有較大柔度。以Φ110 mm透鏡為例,溫度變化為30 ℃時(shí),傳統(tǒng)剛性支撐與撓性支撐兩種情況下,透鏡上表面的變形如圖2所示,其中剛性支撐情況下透鏡下面形峰-谷值(PV)和均方根面形(RMS)分別為16.34 nm和4.12 nm,撓性支撐情況下透鏡面形的PV和RMS分別為3.23 nm和0.82 nm。從圖2中可以看出撓性支撐情況下透鏡面形明顯優(yōu)于剛性支撐。此外,還可根據(jù)光學(xué)元件尺寸將該撓性支撐結(jié)構(gòu)設(shè)計(jì)為不同的幾何構(gòu)形,以滿(mǎn)足接觸點(diǎn)處的最大應(yīng)力要求,如圖3所示。

    盡管撓性支撐結(jié)構(gòu)具有上述優(yōu)點(diǎn),但由于其存在撓性環(huán)節(jié),使得光軸處于水平位置時(shí)會(huì)產(chǎn)生一定的偏移,且整個(gè)支撐組件的諧振頻率要低于剛性支撐。為了合理設(shè)計(jì)該結(jié)構(gòu),需要研究光軸偏差、諧振頻率與結(jié)構(gòu)參數(shù)及幾何構(gòu)形的關(guān)系,進(jìn)而使設(shè)計(jì)滿(mǎn)足性能指標(biāo)要求。

    圖1 撓性支撐結(jié)構(gòu)概念

    (a)剛性支撐

    (b)撓性支撐

    Fig.2 Deformation of top surface of lens under rigid and flexure mounting respectively

    (a)三撓性單元

    (b)四撓性單元

    (c)五撓性單元

    (d)六撓性單元

    2.1 撓性單元?jiǎng)偠确治?/p>

    圖1所示的撓性支撐結(jié)構(gòu)中,關(guān)鍵部分為對(duì)稱(chēng)分布的撓性單元,如圖4(a)所示。圖4(b)為其等效模型。對(duì)于圖1所示撓性支撐結(jié)構(gòu),Yolder[5]將其簡(jiǎn)化為彈簧結(jié)構(gòu),如圖5所示,每個(gè)撓性單元具有徑向和切向剛度。將撓性單元簡(jiǎn)化為兩個(gè)方向上的彈簧雖是合理的,但在撓性單元與光學(xué)元件結(jié)合點(diǎn)處存在約束不足的問(wèn)題,使得其等效剛度要小于實(shí)際模型剛度。為此,本文采用懸臂梁等效,如圖4(b)所示,在結(jié)合點(diǎn)處增加了轉(zhuǎn)動(dòng)約束ROTz=0,使得結(jié)構(gòu)成為超靜定形式。撓性單元結(jié)構(gòu)參數(shù)包括:有效長(zhǎng)度l,厚度t,面內(nèi)寬度b。

    (a)撓性單元

    (b)等效模型

    圖5 撓性支撐的彈簧模型

    解除圖4(b)所示超靜定梁上的連接點(diǎn)B和C的約束,其等效靜定結(jié)構(gòu)如圖6所示,其中Xi(i=1,2,3,4)為未知廣義力,F(xiàn)r為徑向力,F(xiàn)t為切向力。

    撓性單元的切向剛度Kt可以直接由下式得到:

    (1)

    式中:δt為C點(diǎn)處的切向位移,E為彈性模量,A為截面面積。

    圖6 超靜定梁的等效靜定結(jié)構(gòu)

    由圖6可知,由徑向力Fr和未知力Xi引起的懸臂梁截面上的彎矩為:

    (2)

    (3)

    (4)

    (5)

    (6)

    (7)

    根據(jù)虛功原理,在C點(diǎn)處由Fr和Xi引起的位移分別為:

    (8)

    (9)

    由B點(diǎn)和C點(diǎn)處位移邊界條件,結(jié)合式(8)和(9)可得:

    (10)

    (11)

    撓性單元的最大彎曲應(yīng)力為:

    (12)

    2.2 撓性支撐結(jié)構(gòu)整體剛度分析

    按撓性單元的懸臂梁等效模型,將圖3所示的撓性支撐結(jié)構(gòu)等效為如圖7所示的模型,若其安裝時(shí)存在夾角α,其受力分析和位移分析結(jié)果如圖8所示。模型中考慮了幾何構(gòu)形的因素,即撓性單元數(shù)量對(duì)剛度的影響。圖8中Fi表示連接點(diǎn)Pi處撓性件受到的外力,可將其分解為徑向力和切向力Fir和Fit,F(xiàn)iy為Fi的豎向分力,δir和δit分別表示連接點(diǎn)Pi處的徑向位移和切向位移,δiy表示豎向位移,δix表示橫向位移,n為撓性單元的數(shù)目i=1,2,…,n。

    圖7 撓性支撐結(jié)構(gòu)等效模型

    (a)力分析

    (b)位移分析

    根據(jù)幾何關(guān)系,豎向分力Fiy(i=1,2,…,n)可表示為:

    (13)

    根據(jù)豎向力平衡關(guān)系,即:

    (14)

    式中G為光學(xué)元件的重力。

    豎向位移δiy可表示為:

    (15)

    由于重力方向豎直向下,撓性支撐結(jié)構(gòu)連接點(diǎn)Pi處將僅有豎向位移,即:

    (16)

    根據(jù)式(16)知,δir和δit的關(guān)系可表示為:

    (17)

    與光學(xué)元件相比,撓性支撐整體具有一定的柔度,因此在推導(dǎo)光學(xué)元件整體剛度時(shí),可假設(shè)光學(xué)元件為一剛性體,則連接點(diǎn)Pi處的位移δiy均相等,即:

    (18)

    由式(16)和(17)得:

    (19)

    由式(1)知,每個(gè)連接點(diǎn)Pi的徑向位移和切向位移同樣滿(mǎn)足胡克定律,即:

    Fir=Krδir,F(xiàn)it=Ktδit.

    (20)

    由式(14),(18)~(20)得:

    (21)

    撓性支撐結(jié)構(gòu)的整體剛度為:

    (22)

    將式(1)和(11)代入式(22)得:

    (23)

    從式(23)可看出,該撓性支撐結(jié)構(gòu)的整體剛度是彈性模量E,厚度t,寬度b,長(zhǎng)度l,撓性單元數(shù)目n及安裝角度α的函數(shù),即:

    Kn=F(E,b,t,l,n,α).

    (24)

    這里要說(shuō)明,盡管剛度表達(dá)式(23)是按豎向位移推導(dǎo)而來(lái),但其具有一般性,即其它方向的剛度可以通過(guò)改變?chǔ)林档玫健?/p>

    此外,若考慮到光學(xué)元件或其他構(gòu)件自身的變形,則撓性支撐結(jié)構(gòu)的整體剛度小于由式(24)表示的剛度。為此,在整體剛度模型,即式(24)中引入了修正因子λ,即:

    Kn=λ·F(E,b,t,l,n,α),

    (25)

    式中:λ大小由與撓性元件相連構(gòu)件的剛度決定,完全剛性時(shí)λ=1。

    作為特例,b=8,t=0.5,l=20,λ=1,材料為65 Mn,整體剛度Kn隨安裝角度α及撓性單元數(shù)目n的變化曲線(xiàn)如圖9所示。

    從圖9中可以看出,隨安裝角度α變化撓性支撐的整體剛度保持不變。而保持參數(shù)不變情況下,隨著連接點(diǎn)數(shù)目增多,整體剛度增大。

    圖9 整體剛度與安裝角度及連接點(diǎn)數(shù)目關(guān)系

    Fig.9 Relationship between integral stiffness with mounting angle and the number of connecting points

    3 考慮幾何參數(shù)及構(gòu)形的撓性支撐結(jié)構(gòu)優(yōu)化設(shè)計(jì)

    對(duì)于一個(gè)光學(xué)元件,如透鏡或反射鏡,設(shè)計(jì)撓性支撐結(jié)構(gòu),涉及選擇合適的材料,確定其最佳構(gòu)形及幾何參數(shù),以滿(mǎn)足光學(xué)指標(biāo)要求。

    以往的研究主要是先假定一種初始構(gòu)形,然后優(yōu)化幾何參數(shù),若得到的幾何參數(shù)不符合工程應(yīng)用實(shí)際,則改變其構(gòu)形,重復(fù)上述過(guò)程。然而,這樣的優(yōu)化并不能保證幾何構(gòu)形及參數(shù)的最佳匹配。為此,本文提出了一種幾何構(gòu)形及參數(shù)協(xié)同優(yōu)化的設(shè)計(jì)方法。通常,撓性支撐結(jié)構(gòu)材料選擇熱脹系數(shù)較小的銦鋼,其設(shè)計(jì)一般需要滿(mǎn)足以下條件:光軸水平放置時(shí),光軸的偏移量要盡可能小,每個(gè)撓性單元的最大應(yīng)力要滿(mǎn)足其強(qiáng)度條件,每個(gè)連接點(diǎn)處的應(yīng)力要滿(mǎn)足光學(xué)材料的強(qiáng)度條件,此外,面形精度要盡可能高。

    根據(jù)式(23)和上述設(shè)計(jì)要求知,該撓性支撐結(jié)構(gòu)的一般模型可表達(dá)為:

    minΦ(E,b,t,l,n)

    (26)

    式中:Φ為目標(biāo)函數(shù),可根據(jù)需要選擇不同指標(biāo)作為目標(biāo)函數(shù),σf和σfmax分別為撓性單元的應(yīng)力和其材料的最大許用應(yīng)力,σg和σgmax分別為連接點(diǎn)處應(yīng)力及光學(xué)材料的許用應(yīng)力,b0,t0,l0,bmax,tmax,lmax分別為b,t,l的最小值與最大值。

    由于撓性單元數(shù)目n無(wú)法加工為任意數(shù)目,且非整數(shù)的連接點(diǎn)數(shù)也無(wú)意義,因此n僅能取整數(shù)值。由式(23)表達(dá)的剛度模型是非常復(fù)雜的諧波函數(shù),此類(lèi)函數(shù)的優(yōu)化非常困難。仔細(xì)觀察發(fā)現(xiàn),若n取為整數(shù)值,則剛度表達(dá)式(23)也變?yōu)榕cn相對(duì)應(yīng)的離散值,為此,在優(yōu)化模型中引入0-1型整數(shù)變量xi(i=n0,n1,…,nmax)來(lái)簡(jiǎn)化剛度表示式(23)。

    首先引入徑向剛度系數(shù)Krn和切向剛度系數(shù)Ktn,將式(23)表達(dá)為:

    (27)

    式中:Krn,Ktn為由n確定的常數(shù)。

    然后利用設(shè)計(jì)變量xi(i=n0,n1,…,nmax),將撓性支撐結(jié)構(gòu)整體剛度表示為:

    xn1Kn1+…+xnmaxKnmax,

    (28)

    式中:設(shè)計(jì)變量xi取值為0或1,有且僅有一個(gè)xi取值為1,其余為0。

    利用式(28),優(yōu)化模型(26)可重新表達(dá)為:

    minΦ(E,b,t,l,xi)

    (29)

    式(29)中不含有諧波函數(shù),簡(jiǎn)化了模型的復(fù)雜性,使求解過(guò)程變得簡(jiǎn)單。本文采用Lingo軟件整數(shù)和實(shí)數(shù)混合規(guī)劃程序完成上述模型的求解計(jì)算。

    4 數(shù)值模擬與實(shí)驗(yàn)測(cè)試

    為了驗(yàn)證提出的撓性支撐結(jié)構(gòu)的力學(xué)模型及設(shè)計(jì)撓性支撐結(jié)構(gòu)幾何構(gòu)形及參數(shù)的協(xié)同優(yōu)化算法,進(jìn)行了如下數(shù)值仿真及實(shí)驗(yàn)。

    4.1 撓性支撐結(jié)構(gòu)整體剛度仿真與實(shí)驗(yàn)

    以3個(gè)撓性單元情況為例,對(duì)撓性支撐結(jié)構(gòu)進(jìn)行仿真和實(shí)驗(yàn)測(cè)試。撓性支撐結(jié)構(gòu)試驗(yàn)件和有限元模型如圖10所示,結(jié)構(gòu)參數(shù)如表1所示。這里,撓性支撐與實(shí)際中用于支撐光學(xué)元件的結(jié)構(gòu)有一些差別,主要為了方便測(cè)試其剛度。測(cè)試內(nèi)容包括:(1)在結(jié)構(gòu)一側(cè)施加拉力,拉力大小由力傳感器監(jiān)視,采用直線(xiàn)光柵尺測(cè)量輸出位移,根據(jù)位移隨力的變化曲線(xiàn)得到結(jié)構(gòu)整體剛度,圖11所示x軸方向;(2)轉(zhuǎn)動(dòng)撓性支撐結(jié)構(gòu),使其具有一定的安裝角度,測(cè)試其它方向上的剛度,驗(yàn)證結(jié)構(gòu)剛度與安裝角度的關(guān)系,圖11所示y向和Ⅰ及Ⅱ方向。實(shí)驗(yàn)測(cè)試裝置如圖12所示,其中測(cè)力傳感器精度為0.1 N,光柵尺的測(cè)量精度為0.1 μm。

    圖10 撓性支撐結(jié)構(gòu)件及其有限元模型

    材料參數(shù)幾何參數(shù)彈性模量E/GPa210撓性單元數(shù)目n3泊松比ν0.3撓性單元厚度t/mm0.5密度ρ/(kg·m-3)7800撓性單元長(zhǎng)度l/mm35撓性單元寬度b/mm4

    圖11 剛度測(cè)試方向

    圖12 實(shí)驗(yàn)測(cè)試裝置

    分別對(duì)撓性結(jié)構(gòu)件x軸,y軸和Ⅰ方向施加0~400 N的作用力(Ⅱ方向與Ⅰ方向反對(duì)稱(chēng)),間隔為20 N,并在有限元模型及本文推導(dǎo)的理論模型上施加同樣的力,測(cè)量和計(jì)算得到的數(shù)據(jù)如表2~4所示。

    表2 試驗(yàn)件x向剛度測(cè)試數(shù)據(jù)

    表3 試驗(yàn)件y向剛度測(cè)試數(shù)據(jù)

    表4 試驗(yàn)件Ⅰ向剛度測(cè)試數(shù)據(jù)

    應(yīng)用最小二乘法分別對(duì)3個(gè)方向上剛度測(cè)試得到的離散位移數(shù)據(jù)進(jìn)行線(xiàn)性擬合,得到的擬合結(jié)果、數(shù)值仿真及理論計(jì)算結(jié)果如圖13~15所示。

    從圖13~15可以看出,3個(gè)方向上位移測(cè)試數(shù)據(jù)與有限元仿真數(shù)據(jù)吻合的非常好。然而,從圖13~15中還可以看出,由整體剛度模型(23)(即修正因子λ=1)計(jì)算得到的結(jié)果明顯小于實(shí)測(cè)及有限元仿真結(jié)果,其原因在于推導(dǎo)式(23)時(shí),假設(shè)除撓性單元外其余部分為剛體,而試驗(yàn)件的外框架并非理想剛體,因此實(shí)際撓性支撐結(jié)構(gòu)的整體剛度低于理論計(jì)算結(jié)果。通過(guò)比較發(fā)現(xiàn),當(dāng)修正因子λ=0.65時(shí),3個(gè)方向上的理論結(jié)果與實(shí)驗(yàn)及仿真結(jié)果完全吻合。

    此外,從3個(gè)方向的測(cè)試數(shù)據(jù)可看出,3個(gè)方向上的位移數(shù)據(jù)基本一致,證明此撓性件在各個(gè)方向上的剛度是一致的,與本文推導(dǎo)的理論模型結(jié)論一致。

    圖13 x方向數(shù)據(jù)結(jié)果

    圖14 y方向數(shù)據(jù)結(jié)果

    圖15 Ⅰ方向數(shù)據(jù)結(jié)果

    4.2 撓性支撐結(jié)構(gòu)的優(yōu)化設(shè)計(jì)

    以透鏡支撐為例,透鏡的結(jié)構(gòu)參數(shù)如表5所示,設(shè)計(jì)該透鏡的撓性支撐結(jié)構(gòu),使其面形精度達(dá)到最優(yōu),且在不同仰角位置光軸晃動(dòng)量小于0.5 μm。由于透鏡面形精度很難用解析式表達(dá),通常是采用有限元方法計(jì)算其變形,然后通過(guò)曲面擬合得到面形精度值[16-17],這使得優(yōu)化模型時(shí)無(wú)法直接利用面形精度作為目標(biāo)函數(shù)或約束條件,而只能采用間接的方法。已有研究表明,面形精度隨透鏡與撓性單元作用力的增大而變差,但并未考慮撓性單元數(shù)目的影響。為此,本文以撓性支撐結(jié)構(gòu)的總變形能作為目標(biāo)函數(shù),這樣能夠同時(shí)考慮撓性單元與透鏡的相互作用力及撓性單元數(shù)目對(duì)面形精度的影響,材料選取銦鋼4J32,參數(shù)如表6所示,將光軸位移作為約束條件,則優(yōu)化模型可表達(dá)為:

    (30)

    式中:δT是溫差為ΔT是透鏡的徑向變形,撓性件幾何參數(shù)約束由加工能力決定,這里為簡(jiǎn)化模型取l=35 mm。

    表5 透鏡參數(shù)

    表6 銦鋼材料參數(shù)

    以3個(gè)撓性連接點(diǎn)情況為初始設(shè)計(jì),首先調(diào)整撓性支撐結(jié)構(gòu)參數(shù),即厚度t和寬度b使得光軸水平時(shí),光軸偏移滿(mǎn)足要求,并以此為初始設(shè)計(jì)參數(shù)。應(yīng)用本文提出的優(yōu)化模型即式(30)對(duì)透鏡撓性支撐結(jié)構(gòu)進(jìn)行優(yōu)化,確定最佳的撓性單元數(shù)目及結(jié)構(gòu)參數(shù)。

    針對(duì)初始和最優(yōu)情況下的結(jié)構(gòu)參數(shù),分別建立相應(yīng)的有限元模型,并分析在環(huán)境溫差為50 ℃時(shí),透鏡上表面的變形,其中有限元分析在PATRAN/NASTRAN軟件中實(shí)現(xiàn)。然后,分別對(duì)得到的變形數(shù)據(jù)進(jìn)行擬合,除去剛體位移項(xiàng),計(jì)算透鏡表面的面形誤差RMS值。由于理想的溫差環(huán)境及在此環(huán)境下檢測(cè)透鏡面形需要的實(shí)驗(yàn)條件非常復(fù)雜,在此僅對(duì)其進(jìn)行了數(shù)值仿真分析。

    初始設(shè)計(jì)參數(shù)與最佳設(shè)計(jì)參數(shù)如表7所示,其中面形誤差RMS值為透鏡上表面面形,初始和最優(yōu)情況下透鏡的變形如圖16(彩圖見(jiàn)期刊電子版)所示。

    表7 初始及最佳參數(shù)

    圖16 初始情況和最佳情況下透鏡變形云圖

    模型的優(yōu)化目標(biāo)函數(shù)為環(huán)境存在溫差ΔT時(shí)撓性件的總變形能,從表7可以看出,最優(yōu)情況下的變形能及面形誤差均小于初始3個(gè)連接點(diǎn)情況,表明撓性件變形能可以反映透鏡的面形誤差,以此為目標(biāo)函數(shù)來(lái)簡(jiǎn)化優(yōu)化模型是合理的。此外,通過(guò)有限元分析得到的透鏡面形誤差RMS值可以看出,通過(guò)優(yōu)化透鏡的面形精度提高了23%。

    5 結(jié) 論

    本文針對(duì)光學(xué)元件撓性支撐結(jié)構(gòu)設(shè)計(jì)中剛度和柔度要求之間的矛盾,建立了光學(xué)元件撓性支撐結(jié)構(gòu)的力學(xué)模型及優(yōu)化設(shè)計(jì)模型。根據(jù)撓性支撐結(jié)構(gòu)的特性,將各個(gè)撓性單元簡(jiǎn)化為超靜定的梁,并推導(dǎo)了撓性單元的徑向及切向剛度。假設(shè)光學(xué)元件為剛性體,根據(jù)力平衡條件及變形協(xié)調(diào)條件,推導(dǎo)了撓性支撐結(jié)構(gòu)的整體剛度,并引入了修正因子補(bǔ)償由剛體假設(shè)帶來(lái)的誤差。以撓性結(jié)構(gòu)總變形能為目標(biāo)函數(shù),推導(dǎo)了同時(shí)考慮撓性支撐結(jié)構(gòu)幾何構(gòu)形及參數(shù)的協(xié)同優(yōu)化設(shè)計(jì)模型,模型中通過(guò)引入了整型變形將結(jié)構(gòu)整體剛度簡(jiǎn)化為整型變形和離散剛度的線(xiàn)性組合,從而使優(yōu)化模型中不含有諧波函數(shù)項(xiàng)。從數(shù)值仿真和實(shí)驗(yàn)角度對(duì)整體剛度模型進(jìn)行驗(yàn)證,實(shí)驗(yàn)結(jié)果與仿真及理論計(jì)算結(jié)果吻合的非常好。此外,透鏡撓性支撐結(jié)構(gòu)優(yōu)化設(shè)計(jì)結(jié)果表明,通過(guò)優(yōu)化,透鏡面形精度較初始設(shè)計(jì)提高了23%。

    [1] SCHWERTZ K, BURGE J.FieldGuidetoOptomechanicalDesignandAnalysis[M]. Washington USA: SPIE Press, 2012.

    [2] PAROS J M, WEISBORD L. How to design flexure hinge [J].MachineDesign, 1965, 37(1):151-157.

    [3] CAI K H, TIAN Y L, WANG F J,etal.. Development of a piezo-driven 3-DOF stage with T-shape flexible hinge mechanism [J].RoboticsandComputer-IntegratedManufacturing, 2016, 37(2):125-138.

    [4] XU Q. Design, testing and precision control of a novel long-stroke flexure micropositioning system [J].MechanismandMachineTheory, 2013, 70:209-224.

    [5] YOLDER P R.Opto-MechanicalSystemsDesign[M]. Third ed. Washington: CRC Press, 2006.

    [6] YOLDER P R.MountingOpticsinOpticalInstruments[M]. Second ed,Washington:SPIE Press,2008.

    [7] FATA R, KRADINOV V, FABRICANT D. Flexure mounts for high-performance astronomical lenses [J].Proc.SPIE, 2006,6269,doi:101117/12.672634.

    [8] 范磊, 張景旭, 吳小霞, 等. 大口徑輕量化主鏡邊緣側(cè)向支撐的優(yōu)化設(shè)計(jì)[J]. 光學(xué) 精密工程,2012,20(10):2207-2213. FAN L, ZHANG J X, WU X X,etal.. Optimum design of edge lateral support for large aperture lightweight primary mirror [J].Opt.PrecisionEng., 2012, 20(10): 2207-2213.(in Chinese)

    [9] 陳夫林, 張景旭, 吳小霞,等. 620 mm薄鏡面的主動(dòng)支撐結(jié)構(gòu)及面形校正[J]. 光學(xué) 精密工程,2011,19(5):1022-1029. CHEN F L, ZHANG J X, WU X X,etal.. Supporting structure of 620 mm thin primary mirror and its active surface correction [J].Opt.PrecisionEng., 2011, 19(5): 1022-1029.(in Chinese)

    [10] 趙磊,鞏巖. 光刻物鏡中透鏡高精度支撐結(jié)構(gòu)的設(shè)計(jì)及分析[J]. 光學(xué)學(xué)報(bào),2012,32(9):0922001-0922009. ZHAO L, GONG Y. Design and analysis for the high-precision lens support structure of ojective lens for lithography [J].ActaOpticaSinica, 2012, 32(9): 0922001-0922009.(in Chinese)

    [11] 倪明陽(yáng),鞏巖. 光刻投影物鏡光學(xué)元件運(yùn)動(dòng)學(xué)支撐結(jié)構(gòu)的設(shè)計(jì)與分析[J]. 中國(guó)光學(xué),2012,5(5):476-484. NI M Y, GONG Y. Design and analysis of kinematic lens positioning structure in lithographic projection objective [J].ChineseOptics, 2012, 5(5): 476-484.(in Chinese)

    [12] 馬磊,曹佃生,劉承志. 大口徑透鏡多點(diǎn)柔性支撐結(jié)構(gòu)設(shè)計(jì)與分析[J]. 光電工程,2015,42(5):88-94. MA L, CAO D SH, LIU CH ZH. Design and analysis on multi-points flexible supporting structure of large-aperture lens [J].Opto-ElectronicEngineering, 2015, 42(5): 88-94.(in Chinese)

    [13] 趙磊, 彭海峰, 于新峰,等. 具有整體徑向撓性的光學(xué)元件多點(diǎn)均勻支撐結(jié)構(gòu)[J]. 光電工程,2015,42(10):43-48. ZHAO L, PENG H F, YU X F,etal.. Lens support structure of multi point equal supporting with wholly radical freedom [J].Opto-ElectronicEngineering, 2015, 42(10): 43-48.(in Chinese)

    [14] 李宇軒, 陳雪, 張雷,等. 大口徑空間反射鏡Cartwheel型柔性鉸鏈支撐設(shè)計(jì)[J]. 光學(xué)學(xué)報(bào),2014,34(6):0622003-0622009. LI Y X, CHEN X, ZHANG L,etal.. Design of Cartwheel flexural support for a large aperture space mirror [J].ActaOpticaSinica, 2014, 34(6): 0622003-0622009.(in Chinese)

    [15] 曹乃亮, 徐宏, 辛宏偉, 等. 基于NiTi合金絲的反射鏡柔性支撐結(jié)構(gòu)的應(yīng)力補(bǔ)償[J]. 光學(xué) 精密工程,2012,20(10):2161-2169. CAO N L, XU H, XIN H W,etal.. Stress compensation of flexible supporting structures for mirrors using NiTi shape memory alloy [J].Opt.PrecisionEng., 2012, 20(10): 2161-2169.(in Chinese)

    [16] 曹玉巖, 王志, 周超,等. 壓電智能反射面靜態(tài)形狀控制與作動(dòng)器位置優(yōu)化[J]. 航空學(xué)報(bào),2015,36(2):527-537. CAO Y Y, WANG ZH, ZHOU CH,etal.. Static shape control of piezoelectric smart reflector and optimization of actuators placement [J].ActaAeronauticaetAstronauticaSinica, 2015, 36(2): 527-537.(in Chinese)

    [17] 王志, 曹玉巖, 周超,等. 柔性壓電智能反射面的靜態(tài)形狀控制[J]. 光學(xué) 精密工程,2014,22(10):2715-2724. WANG ZH, CAO Y Y, ZHOU CH,etal.. Static shape control of flexible piezoelectric smart reflectors [J].Opt.PrecisionEng., 2014, 22(10): 2715-2724.(in Chinese)

    曹玉巖(1986-),男,吉林大安人,助理研究員,2009年、2012年于西安電子科技大學(xué)分別獲學(xué)士、碩士學(xué)位,主要從事結(jié)構(gòu)有限元理論、結(jié)構(gòu)振動(dòng)控制技術(shù)研究。E-mail:yuyan_cao@126.com

    王志臣(1980-),男,黑龍江大興安嶺人,副研究員,主要從事望遠(yuǎn)鏡光機(jī)結(jié)構(gòu)設(shè)計(jì)研究。E-mail:zcwang911@163.com

    (本欄目編輯:馬 健)

    (版權(quán)所有 未經(jīng)許可 不得轉(zhuǎn)載)

    General modeling and optimal design of flexure supporting structure for optical components

    CAO Yu-yan*, WANG Zhi-chen, ZHOU Chao, FAN Lei, HAN Xi-da, ZHANG Yao-zu

    (Changchun Institute of Optics, Fine Mechanics and Physics,ChineseAcademyofSciences,Changchun130033,China)

    A mechanical and a parameter optimization model for flexure support structure of optical components were proposed to allow the flexure structure meet simultaneously the requirements of the stiffness for optical component position and the compliance for temperature adaptability。Meanwhile, the corresponding modeling method was investigated. As this flexure structure was consisted of several identical flexure parts, it was simplified into an indeterminate beam structure, and the radical stiffness and tangential stiffness were derived using the virtual work principle. Then, by assuming optical components for the rigid body, the whole stiffness of the flexure support structure was derived based on the force equilibrium and its compatible deformation, and the correction factor was introduced to compensate the error caused by the rigid assumption. Finally, the total strain energy of the flexure structure was taken as the objective function, and the collaborative optimization model was derived considering the geometrical pattern and parameters simultaneously. By introducing the integral variables, the whole stiffness of the structure was simplified into a linear combination of the integral variable and discrete stiffness, and the harmonic terms were eliminated. The whole stiffness model was verified by the simulation and experiment, and the experiment results are highly in agreement with the simulation results. A lens mounting was taken for an example, the optimization method of the flexure mounting structure was verified. The finite element simulation results show that the surface precision of the lens has been improved by 23%.

    optical component; flexure mounting structure; general model; indetermination; strain energy; optimization design

    2015-12-23;

    2016-02-04.

    國(guó)家自然科學(xué)基金資助項(xiàng)目(No.11403023)

    1004-924X(2016)11-2792-12

    TH703;TP273

    A

    10.3788/OPE.20162411.2792

    *Correspondingauthor,E-mail:yuyan_cao@126.com

    猜你喜歡
    面形撓性透鏡
    “透鏡及其應(yīng)用”知識(shí)延伸
    反射鏡面形三坐標(biāo)白光掃描檢測(cè)精度研究
    “透鏡”知識(shí)鞏固
    “透鏡及其應(yīng)用”知識(shí)拓展
    “透鏡”知識(shí)鞏固
    使用最小二乘迭代相移方法測(cè)量透明元件*
    超光譜儀光柵面形對(duì)光譜性能影響的仿真分析
    什么是撓性轉(zhuǎn)子軸承
    剛性轉(zhuǎn)子和撓性轉(zhuǎn)子的區(qū)別
    采樣點(diǎn)分布對(duì)基于面形斜率徑向基模型的自由曲面擬合精度的影響
    91精品一卡2卡3卡4卡| 激情五月婷婷亚洲| 日日摸夜夜添夜夜爱| 97精品久久久久久久久久精品| 国产亚洲av片在线观看秒播厂| 国产一区二区三区综合在线观看 | av国产精品久久久久影院| 国产人妻一区二区三区在| 精品亚洲成a人片在线观看 | 国产中年淑女户外野战色| 日韩视频在线欧美| 亚洲av欧美aⅴ国产| 国产欧美另类精品又又久久亚洲欧美| 一区二区三区乱码不卡18| 久久久久久久国产电影| 国产在线免费精品| 少妇丰满av| 人人妻人人添人人爽欧美一区卜 | 大话2 男鬼变身卡| 日韩三级伦理在线观看| 国产亚洲欧美精品永久| 丰满迷人的少妇在线观看| 多毛熟女@视频| 亚洲美女视频黄频| 美女xxoo啪啪120秒动态图| 狂野欧美白嫩少妇大欣赏| 国产免费又黄又爽又色| 最近中文字幕2019免费版| 亚洲欧洲国产日韩| 国模一区二区三区四区视频| 最近最新中文字幕大全电影3| 久久女婷五月综合色啪小说| 亚洲av男天堂| 大香蕉97超碰在线| 亚洲丝袜综合中文字幕| 91久久精品电影网| 搡老乐熟女国产| 久久久国产一区二区| 一级a做视频免费观看| 亚洲av免费高清在线观看| 中文字幕人妻熟人妻熟丝袜美| 女性生殖器流出的白浆| 菩萨蛮人人尽说江南好唐韦庄| 91精品国产国语对白视频| 国产精品蜜桃在线观看| av天堂中文字幕网| 久久久久久伊人网av| 一区二区三区四区激情视频| 国产深夜福利视频在线观看| 欧美日韩一区二区视频在线观看视频在线| 国产男人的电影天堂91| 国产中年淑女户外野战色| 另类亚洲欧美激情| 日本av免费视频播放| 少妇的逼好多水| 舔av片在线| 国产精品一区www在线观看| 三级国产精品片| 日韩欧美一区视频在线观看 | 日本wwww免费看| 国产成人a∨麻豆精品| 纯流量卡能插随身wifi吗| 不卡视频在线观看欧美| 一区二区av电影网| 男人和女人高潮做爰伦理| 一个人免费看片子| 99精国产麻豆久久婷婷| 成人黄色视频免费在线看| 大话2 男鬼变身卡| 久久亚洲国产成人精品v| 大话2 男鬼变身卡| 99热网站在线观看| 久久青草综合色| 大香蕉97超碰在线| 国产精品国产三级国产专区5o| 在线观看人妻少妇| 夫妻性生交免费视频一级片| 一区二区三区精品91| 日韩伦理黄色片| 最新中文字幕久久久久| 亚洲自偷自拍三级| 黄色怎么调成土黄色| 精品久久久久久久末码| 亚洲欧美日韩另类电影网站 | 秋霞伦理黄片| 久久国产精品大桥未久av | 亚洲熟女精品中文字幕| 国产成人精品婷婷| 男男h啪啪无遮挡| 午夜激情福利司机影院| av在线app专区| 在线精品无人区一区二区三 | a 毛片基地| 80岁老熟妇乱子伦牲交| 91aial.com中文字幕在线观看| 久久精品久久精品一区二区三区| 最近最新中文字幕大全电影3| av在线蜜桃| 免费久久久久久久精品成人欧美视频 | 亚洲精品国产av蜜桃| 亚洲中文av在线| 久久99热这里只有精品18| 一区二区三区乱码不卡18| 夜夜爽夜夜爽视频| 欧美亚洲 丝袜 人妻 在线| 人人妻人人添人人爽欧美一区卜 | 只有这里有精品99| 18禁在线无遮挡免费观看视频| videos熟女内射| 精品熟女少妇av免费看| 青春草亚洲视频在线观看| 国产精品久久久久久精品古装| 伦理电影大哥的女人| 精品久久久久久电影网| 水蜜桃什么品种好| 免费大片黄手机在线观看| 国产精品av视频在线免费观看| 18禁裸乳无遮挡免费网站照片| 在线观看一区二区三区激情| av播播在线观看一区| 一区二区三区四区激情视频| 大又大粗又爽又黄少妇毛片口| 亚洲精品视频女| 偷拍熟女少妇极品色| 女的被弄到高潮叫床怎么办| 男女免费视频国产| 我要看日韩黄色一级片| 精品亚洲成a人片在线观看 | 高清在线视频一区二区三区| 中文字幕av成人在线电影| 丰满迷人的少妇在线观看| 在线亚洲精品国产二区图片欧美 | 国产精品一区二区三区四区免费观看| 国产男人的电影天堂91| 女的被弄到高潮叫床怎么办| 中文天堂在线官网| 91久久精品国产一区二区成人| 极品少妇高潮喷水抽搐| 男人舔奶头视频| 97在线视频观看| 国产av国产精品国产| 亚洲精品aⅴ在线观看| 夜夜看夜夜爽夜夜摸| 日韩av在线免费看完整版不卡| 亚洲精品乱久久久久久| 高清不卡的av网站| a级毛色黄片| 男女无遮挡免费网站观看| 欧美日韩综合久久久久久| 美女主播在线视频| 久久影院123| 51国产日韩欧美| 日本欧美视频一区| 精品酒店卫生间| 国产亚洲5aaaaa淫片| 啦啦啦在线观看免费高清www| 好男人视频免费观看在线| 国产91av在线免费观看| 国产伦在线观看视频一区| 最新中文字幕久久久久| 在线免费观看不下载黄p国产| 尾随美女入室| 日韩不卡一区二区三区视频在线| 国产亚洲精品久久久com| 少妇的逼水好多| 精华霜和精华液先用哪个| av国产久精品久网站免费入址| 青春草视频在线免费观看| 免费人妻精品一区二区三区视频| 王馨瑶露胸无遮挡在线观看| 国产色婷婷99| 亚洲精品日韩av片在线观看| 热99国产精品久久久久久7| 欧美高清成人免费视频www| 久久久久精品久久久久真实原创| 国产美女午夜福利| 亚洲成色77777| 在线观看三级黄色| 91精品国产九色| 国产男女超爽视频在线观看| 精品久久久久久久久av| 直男gayav资源| 身体一侧抽搐| 一二三四中文在线观看免费高清| 成人午夜精彩视频在线观看| 最近最新中文字幕免费大全7| 久久人人爽人人片av| 亚洲成色77777| 国产精品三级大全| 极品少妇高潮喷水抽搐| 欧美日韩视频精品一区| 乱系列少妇在线播放| 国产精品一区二区三区四区免费观看| 狂野欧美激情性xxxx在线观看| 特大巨黑吊av在线直播| 91午夜精品亚洲一区二区三区| 欧美zozozo另类| 免费人成在线观看视频色| 中文欧美无线码| 亚洲精品国产av成人精品| 人人妻人人看人人澡| 久久99热这里只有精品18| 久久久久国产精品人妻一区二区| 99视频精品全部免费 在线| 亚洲精品国产成人久久av| 狠狠精品人妻久久久久久综合| 在线看a的网站| 综合色丁香网| 秋霞在线观看毛片| 熟女电影av网| 啦啦啦中文免费视频观看日本| tube8黄色片| 国产精品偷伦视频观看了| 欧美日韩国产mv在线观看视频 | 久久久国产一区二区| 国产精品久久久久久久久免| 日韩av在线免费看完整版不卡| 99久久综合免费| 午夜视频国产福利| 一本色道久久久久久精品综合| 国产精品国产三级国产专区5o| 九九在线视频观看精品| 啦啦啦视频在线资源免费观看| 中文字幕久久专区| 国产成人91sexporn| 秋霞在线观看毛片| 永久免费av网站大全| 最近中文字幕高清免费大全6| 亚洲三级黄色毛片| 国产伦精品一区二区三区视频9| 熟女人妻精品中文字幕| 国产av精品麻豆| 成人影院久久| 国产 一区 欧美 日韩| 中文欧美无线码| 大香蕉97超碰在线| 亚洲成人av在线免费| 一级毛片电影观看| 日日啪夜夜爽| 欧美日韩精品成人综合77777| 黄色视频在线播放观看不卡| 欧美三级亚洲精品| 毛片一级片免费看久久久久| 久久久久久久久久成人| 日韩强制内射视频| 精品一区二区三卡| 最近最新中文字幕免费大全7| 超碰97精品在线观看| av在线蜜桃| 91精品伊人久久大香线蕉| 国产在线男女| 亚洲内射少妇av| 久久精品国产亚洲av天美| 国产极品天堂在线| 一级毛片我不卡| 国产成人freesex在线| 日本猛色少妇xxxxx猛交久久| 特大巨黑吊av在线直播| 日产精品乱码卡一卡2卡三| 插阴视频在线观看视频| 亚洲中文av在线| 男女边吃奶边做爰视频| 亚洲三级黄色毛片| av又黄又爽大尺度在线免费看| 如何舔出高潮| 自拍偷自拍亚洲精品老妇| 国产老妇伦熟女老妇高清| 伦理电影大哥的女人| 舔av片在线| 欧美97在线视频| 国产精品av视频在线免费观看| 欧美丝袜亚洲另类| 国产乱人视频| 成年人午夜在线观看视频| 亚洲成色77777| 超碰av人人做人人爽久久| 青春草视频在线免费观看| 少妇人妻久久综合中文| 如何舔出高潮| 国产亚洲5aaaaa淫片| 国产在视频线精品| 亚洲欧美日韩东京热| 欧美精品国产亚洲| 国产日韩欧美亚洲二区| 久久精品国产自在天天线| 国产美女午夜福利| 永久网站在线| 99久久综合免费| 在线播放无遮挡| 久久久久视频综合| 国产爽快片一区二区三区| 亚洲精品视频女| 国产高清有码在线观看视频| 中文资源天堂在线| 亚洲综合色惰| 青春草国产在线视频| 人人妻人人爽人人添夜夜欢视频 | 热99国产精品久久久久久7| 亚洲三级黄色毛片| 高清日韩中文字幕在线| 少妇人妻精品综合一区二区| 99精国产麻豆久久婷婷| 久久 成人 亚洲| 欧美高清成人免费视频www| 亚洲一级一片aⅴ在线观看| 在线 av 中文字幕| 国产av码专区亚洲av| 一级毛片久久久久久久久女| 久热这里只有精品99| 久久99热这里只有精品18| 一本一本综合久久| 色综合色国产| 久久ye,这里只有精品| av在线老鸭窝| www.色视频.com| 热re99久久精品国产66热6| 亚洲av日韩在线播放| 日日摸夜夜添夜夜爱| 久久亚洲国产成人精品v| 亚洲欧洲日产国产| 亚洲国产毛片av蜜桃av| 爱豆传媒免费全集在线观看| 在线观看人妻少妇| 亚洲精品乱久久久久久| 成人美女网站在线观看视频| 大码成人一级视频| 国产日韩欧美在线精品| 一级片'在线观看视频| 亚洲国产毛片av蜜桃av| 精品午夜福利在线看| 99久久人妻综合| 国产一区二区在线观看日韩| 久久久久久人妻| 哪个播放器可以免费观看大片| 日韩一区二区三区影片| 国产精品成人在线| 韩国高清视频一区二区三区| 亚洲欧美精品专区久久| 91精品国产九色| 国产真实伦视频高清在线观看| 99热这里只有是精品在线观看| 亚洲av欧美aⅴ国产| 大香蕉97超碰在线| 成人毛片60女人毛片免费| 午夜福利影视在线免费观看| 亚洲第一区二区三区不卡| 亚洲aⅴ乱码一区二区在线播放| 午夜免费男女啪啪视频观看| 日本猛色少妇xxxxx猛交久久| 日韩视频在线欧美| 亚洲欧洲日产国产| 蜜桃在线观看..| 国产亚洲午夜精品一区二区久久| a级一级毛片免费在线观看| 免费久久久久久久精品成人欧美视频 | 黄色配什么色好看| 女人十人毛片免费观看3o分钟| 成人综合一区亚洲| 两个人的视频大全免费| 交换朋友夫妻互换小说| 18禁裸乳无遮挡动漫免费视频| 国产老妇伦熟女老妇高清| 在线精品无人区一区二区三 | 国产精品久久久久久精品电影小说 | 国产乱人视频| 在线观看免费高清a一片| 免费人成在线观看视频色| 久久久午夜欧美精品| 国产精品女同一区二区软件| av天堂中文字幕网| 我要看日韩黄色一级片| 久久久久久人妻| 亚洲精品日韩在线中文字幕| videossex国产| 中国美白少妇内射xxxbb| 各种免费的搞黄视频| 欧美3d第一页| 成年美女黄网站色视频大全免费 | 永久免费av网站大全| 久久 成人 亚洲| 国产精品国产av在线观看| www.色视频.com| 国产成人aa在线观看| 99国产精品免费福利视频| 又粗又硬又长又爽又黄的视频| 天堂俺去俺来也www色官网| av福利片在线观看| 在线观看国产h片| 丝袜喷水一区| 亚洲美女黄色视频免费看| 一本一本综合久久| 大香蕉久久网| 丰满乱子伦码专区| 国产精品蜜桃在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲国产最新在线播放| 欧美老熟妇乱子伦牲交| 国产精品国产三级国产av玫瑰| 丝袜脚勾引网站| 国产精品.久久久| 我要看日韩黄色一级片| kizo精华| 一级二级三级毛片免费看| 赤兔流量卡办理| 亚洲精品国产色婷婷电影| 久久久欧美国产精品| 18禁在线播放成人免费| 日韩,欧美,国产一区二区三区| 久久人人爽av亚洲精品天堂 | 视频中文字幕在线观看| 女人久久www免费人成看片| 中文精品一卡2卡3卡4更新| 成人美女网站在线观看视频| 日本wwww免费看| 蜜桃久久精品国产亚洲av| 色婷婷av一区二区三区视频| 联通29元200g的流量卡| 2018国产大陆天天弄谢| 在线观看免费视频网站a站| 国国产精品蜜臀av免费| 欧美激情国产日韩精品一区| 亚洲成人手机| av在线蜜桃| 午夜激情福利司机影院| 欧美日韩精品成人综合77777| 老熟女久久久| 国产成人精品婷婷| 日韩成人伦理影院| 日日摸夜夜添夜夜爱| 99精国产麻豆久久婷婷| 国产 精品1| 国产精品偷伦视频观看了| 欧美日韩亚洲高清精品| 丰满人妻一区二区三区视频av| 精品久久久噜噜| 免费av中文字幕在线| 久久精品国产鲁丝片午夜精品| 欧美日韩精品成人综合77777| 国产在线一区二区三区精| 各种免费的搞黄视频| 欧美97在线视频| 成人特级av手机在线观看| 熟女av电影| 亚洲精品久久午夜乱码| 九草在线视频观看| 一级a做视频免费观看| 一个人免费看片子| 草草在线视频免费看| 男人狂女人下面高潮的视频| 91午夜精品亚洲一区二区三区| 国产精品嫩草影院av在线观看| 免费av不卡在线播放| 青青草视频在线视频观看| 80岁老熟妇乱子伦牲交| 国产一区亚洲一区在线观看| 免费黄网站久久成人精品| 深夜a级毛片| 国产中年淑女户外野战色| 身体一侧抽搐| 夜夜骑夜夜射夜夜干| 黄色日韩在线| 日本av手机在线免费观看| 欧美日本视频| 丰满乱子伦码专区| 久热这里只有精品99| 免费久久久久久久精品成人欧美视频 | 三级经典国产精品| 少妇猛男粗大的猛烈进出视频| 亚洲精品第二区| 国产免费一级a男人的天堂| 丰满少妇做爰视频| 国产精品一区二区在线观看99| 观看美女的网站| 伦精品一区二区三区| 深爱激情五月婷婷| 在线播放无遮挡| 伦理电影免费视频| 亚洲国产精品国产精品| 国产午夜精品久久久久久一区二区三区| 97精品久久久久久久久久精品| 欧美精品一区二区大全| 丝袜脚勾引网站| 男女啪啪激烈高潮av片| 久久久久国产精品人妻一区二区| 日韩av在线免费看完整版不卡| 99热这里只有是精品在线观看| 一区二区三区精品91| 国产亚洲一区二区精品| 日本色播在线视频| 欧美高清性xxxxhd video| 亚洲av男天堂| 亚洲国产毛片av蜜桃av| 久久人人爽人人爽人人片va| av专区在线播放| 亚洲内射少妇av| 成人午夜精彩视频在线观看| 2018国产大陆天天弄谢| 国产高清三级在线| 午夜激情久久久久久久| 极品少妇高潮喷水抽搐| 草草在线视频免费看| 中文字幕人妻熟人妻熟丝袜美| 亚洲自偷自拍三级| 欧美一级a爱片免费观看看| 色哟哟·www| 久热这里只有精品99| 少妇的逼好多水| 免费看不卡的av| 午夜免费鲁丝| 国产黄片美女视频| 99久久精品国产国产毛片| 三级国产精品欧美在线观看| 最近2019中文字幕mv第一页| 亚洲精品,欧美精品| 下体分泌物呈黄色| 亚洲精品亚洲一区二区| 黄色视频在线播放观看不卡| 亚洲精品亚洲一区二区| 亚洲av综合色区一区| 国产一区二区在线观看日韩| 亚洲精品国产av蜜桃| 亚洲av男天堂| 内地一区二区视频在线| 欧美日韩精品成人综合77777| 晚上一个人看的免费电影| 少妇人妻久久综合中文| 99热这里只有精品一区| 一本色道久久久久久精品综合| 大又大粗又爽又黄少妇毛片口| a级毛片免费高清观看在线播放| 在线播放无遮挡| 在线观看美女被高潮喷水网站| 九草在线视频观看| 国产免费一区二区三区四区乱码| 一级毛片黄色毛片免费观看视频| 高清欧美精品videossex| 成人毛片a级毛片在线播放| 在线免费观看不下载黄p国产| 国产一区二区三区av在线| 99久久人妻综合| 亚洲国产欧美在线一区| 少妇的逼水好多| 观看美女的网站| 国产精品一区二区三区四区免费观看| 日韩欧美精品免费久久| 欧美精品一区二区大全| 最近2019中文字幕mv第一页| 身体一侧抽搐| 狂野欧美白嫩少妇大欣赏| 国产精品成人在线| 中文字幕制服av| 久久久成人免费电影| 免费少妇av软件| 夫妻午夜视频| 国产成人一区二区在线| 国产精品女同一区二区软件| 黄色配什么色好看| av卡一久久| 成人亚洲欧美一区二区av| 大香蕉97超碰在线| 22中文网久久字幕| 在线免费观看不下载黄p国产| 2022亚洲国产成人精品| 久久精品久久久久久久性| 五月伊人婷婷丁香| 成年av动漫网址| 亚洲一级一片aⅴ在线观看| 男女无遮挡免费网站观看| 一级毛片电影观看| 欧美+日韩+精品| 久久韩国三级中文字幕| 欧美日韩一区二区视频在线观看视频在线| 五月天丁香电影| 亚洲一区二区三区欧美精品| 妹子高潮喷水视频| 男女下面进入的视频免费午夜| 直男gayav资源| 亚洲国产欧美在线一区| 80岁老熟妇乱子伦牲交| 午夜免费鲁丝| 黑丝袜美女国产一区| 亚洲精品日本国产第一区| 国产精品国产三级专区第一集| 欧美成人精品欧美一级黄| 亚洲怡红院男人天堂| 国产色婷婷99| 我要看日韩黄色一级片| 伦理电影大哥的女人| 国产精品熟女久久久久浪| 国产免费一级a男人的天堂| 国产一级毛片在线| 夫妻午夜视频| 国产有黄有色有爽视频| 欧美日韩亚洲高清精品| 精品一区在线观看国产| 最近中文字幕2019免费版| 免费看日本二区| 免费av中文字幕在线| 国产精品久久久久久久电影| 中文字幕亚洲精品专区| 噜噜噜噜噜久久久久久91| 亚洲欧美一区二区三区黑人 | 联通29元200g的流量卡| 18禁在线无遮挡免费观看视频| 久久99蜜桃精品久久| 国产色婷婷99| 精品久久久久久久末码| 深夜a级毛片| 80岁老熟妇乱子伦牲交| 免费观看的影片在线观看| 99热国产这里只有精品6| 久久毛片免费看一区二区三区| 最近手机中文字幕大全| 婷婷色综合大香蕉| 夜夜爽夜夜爽视频|