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

    紅外熱成像折轉(zhuǎn)光學(xué)系統(tǒng)的光軸靜態(tài)敏感度分析

    2023-08-28 09:31:30楊開宇董樹林胡健釧李曉君
    紅外技術(shù) 2023年8期
    關(guān)鍵詞:光軸蒙特卡羅公差

    楊 丹,金 寧,楊開宇,李 晶,董樹林,胡健釧,李曉君

    紅外熱成像折轉(zhuǎn)光學(xué)系統(tǒng)的光軸靜態(tài)敏感度分析

    楊 丹1,金 寧1,楊開宇1,李 晶1,董樹林1,胡健釧1,李曉君2

    (1. 昆明物理研究所,云南 昆明 650223;2. 空裝成都局駐昆明地區(qū)軍代表室,云南 昆明 650223)

    紅外熱成像折轉(zhuǎn)光學(xué)系統(tǒng)在復(fù)雜環(huán)境條件下,光軸容易因為光學(xué)元件的偏心或傾斜而發(fā)生漂移,影響系統(tǒng)對目標(biāo)的指示精度。在紅外熱成像系統(tǒng)設(shè)計之初對光學(xué)系統(tǒng)開展光軸靜態(tài)敏感度分析,能夠識別出光學(xué)系統(tǒng)的敏感點,為滿足光軸穩(wěn)定性的結(jié)構(gòu)優(yōu)化設(shè)計提供約束條件。通過基于旋轉(zhuǎn)矩陣的坐標(biāo)變換,建立了光學(xué)元件旋轉(zhuǎn)過程量和傾斜狀態(tài)量的轉(zhuǎn)換關(guān)系,從而實現(xiàn)了光學(xué)元件在任意方向傾斜的空間姿態(tài)模擬,確保了光軸敏感度公差分析中的蒙特卡羅采樣與結(jié)構(gòu)設(shè)計的約束條件相對應(yīng),并在此基礎(chǔ)上搭建了對紅外折轉(zhuǎn)光學(xué)系統(tǒng)光軸靜態(tài)敏感度分析的流程,編制了程序。用所編程序?qū)δ车湫图t外熱成像折轉(zhuǎn)光學(xué)系統(tǒng)進行實例分析,根據(jù)光軸穩(wěn)定性的指標(biāo)要求,依次對光學(xué)系統(tǒng)中各光學(xué)件的偏心量和傾斜量進行了光軸的靈敏度和反靈敏度分析,得出了初始公差限,再針對初始公差限數(shù)據(jù)進行了任意方向采樣的蒙特卡羅分析,最終得出了各光學(xué)元件能夠滿足光軸穩(wěn)定性指標(biāo)的偏心和傾斜公差限數(shù)據(jù),通過建立多重坐標(biāo)系的方法驗證了所得數(shù)據(jù)的準(zhǔn)確性,為指導(dǎo)光機熱優(yōu)化設(shè)計奠定了基礎(chǔ)。

    紅外折轉(zhuǎn)光學(xué)系統(tǒng);光軸漂移;靜態(tài)敏感度;公差分析;蒙特卡羅采樣

    0 引言

    紅外熱成像系統(tǒng)工作在較為嚴(yán)酷的使用環(huán)境中時,受安裝力、高低溫和沖擊振動等因素影響,內(nèi)部的機械支撐結(jié)構(gòu)產(chǎn)生微小的變形,使光學(xué)元件出現(xiàn)空間狀態(tài)的偏移(包含偏心或傾斜),導(dǎo)致系統(tǒng)的光軸發(fā)生漂移,影響系統(tǒng)對目標(biāo)的指示精度[1-6]。紅外折轉(zhuǎn)光學(xué)系統(tǒng)具有體積小、布局緊湊的優(yōu)勢,但在工程實踐中發(fā)現(xiàn),相較于直線形光學(xué)系統(tǒng),光軸在環(huán)境影響下的漂移量相對更大。為了滿足使用要求,光軸的最大漂移角度即光軸穩(wěn)定性必須限定在允許范圍之內(nèi)。因此,在系統(tǒng)設(shè)計階段,需嚴(yán)格分析各光學(xué)元件的偏移誤差對系統(tǒng)光軸帶來的影響,識別出光學(xué)系統(tǒng)在靜態(tài)下的敏感點,同時基于光軸穩(wěn)定性指標(biāo)分析出各個鏡片偏心和傾斜的合理公差限,從而能夠得到光學(xué)系統(tǒng)所容許各支撐結(jié)構(gòu)件最大的形變范圍,以作為結(jié)構(gòu)優(yōu)化的約束條件[7-14],為光軸穩(wěn)定性的設(shè)計提供指導(dǎo)。

    對于折轉(zhuǎn)光學(xué)系統(tǒng),由于其非旋轉(zhuǎn)對稱性,光學(xué)元件的偏移狀態(tài)需要考慮偏移量和偏移方向兩方面信息。也就是說,在分析某個元件允許的最大偏心量或傾斜量時,需要對所有可能的偏心方向或傾斜方向進行采樣,選取最嚴(yán)苛的情況作為結(jié)構(gòu)優(yōu)化的約束。國內(nèi)外有關(guān)光學(xué)系統(tǒng)公差靈敏度的研究,在考慮元件的偏心量或傾斜量時大多僅針對-軸方向或-軸方向進行采樣,而極少考慮其他傾斜方向,對折轉(zhuǎn)光學(xué)系統(tǒng)的適用性較低[15-17]。

    目前主流的光學(xué)軟件具有完備分析公差的流程,但其對光學(xué)元件傾斜狀態(tài)的建模是先后繞、、各個坐標(biāo)軸的旋轉(zhuǎn)形成,其旋轉(zhuǎn)值為過程量,不能代表旋轉(zhuǎn)后形成的空間狀態(tài)量[18-19]。利用現(xiàn)有軟件進行公差分析時,只能完全按照建模的輸入方式,對光學(xué)元件旋轉(zhuǎn)過程量采樣,無法對其傾斜的空間狀態(tài)量采樣;而對結(jié)構(gòu)優(yōu)化來說,所需的約束條件為狀態(tài)量,這就導(dǎo)致采樣結(jié)果與約束條件不能對應(yīng),對結(jié)構(gòu)設(shè)計的指導(dǎo)存在局限性。

    本文基于公差分析的思想,立足于紅外折轉(zhuǎn)光學(xué)系統(tǒng)元件的偏心、傾斜靜態(tài)公差,建立光學(xué)元件旋轉(zhuǎn)過程量與旋轉(zhuǎn)形成的空間狀態(tài)量(即傾斜量和傾斜方向)之間的轉(zhuǎn)換關(guān)系,實現(xiàn)對光學(xué)元件在任意方向傾斜狀態(tài)的準(zhǔn)確采樣;在此基礎(chǔ)上,建立光學(xué)元件偏移對光軸穩(wěn)定性的影響的靈敏度分析、反靈敏度分析及蒙特卡羅分析整個公差計算流程,以獲得折轉(zhuǎn)光學(xué)系統(tǒng)的靜態(tài)敏感點、以及結(jié)構(gòu)優(yōu)化設(shè)計所需的約束條件。

    1 光學(xué)元件旋轉(zhuǎn)過程量和傾斜狀態(tài)量的轉(zhuǎn)換關(guān)系建立

    蒙特卡羅公差分析是光軸靜態(tài)敏感度分析中最為重要的一步,是最終獲得光學(xué)系統(tǒng)中各鏡片偏心和傾斜合理公差限的必要過程。在紅外折轉(zhuǎn)光學(xué)系統(tǒng)光軸敏感度公差的蒙特卡羅分析中,為了和結(jié)構(gòu)優(yōu)化設(shè)計能夠?qū)?yīng),需在每個元件垂直于光軸的參考平面的360°所有方向上,對傾斜量進行蒙特卡羅采樣,然后根據(jù)采樣的狀態(tài)構(gòu)建出對應(yīng)的傾斜模型,從而能夠?qū)υ撃P瓦M行非順序光線追跡,計算出光軸的最終漂移情況。

    對傾斜模型的構(gòu)建,需要通過旋轉(zhuǎn)元件來完成。在幾何空間中,每個元件自帶一個右手局部坐標(biāo)系-,其中-軸沿光軸方向。元件的旋轉(zhuǎn)運動姿態(tài)采用歐拉角來表述,旋轉(zhuǎn)后的最終姿態(tài)除了與角度相關(guān),還與旋轉(zhuǎn)順序相關(guān)[20]。在對元件的傾斜進行建模時,假設(shè)旋轉(zhuǎn)按圖1所規(guī)定的順序,即元件先繞自身坐標(biāo)軸的-軸左旋角度,再繞旋轉(zhuǎn)后形成的坐標(biāo)軸的-軸左旋角度,再繞旋轉(zhuǎn)后形成的坐標(biāo)軸的-軸右旋角度,得到元件的新坐標(biāo)系-111,則空間中某點在原坐標(biāo)系-和新坐標(biāo)系-111中的坐標(biāo)描述對應(yīng)關(guān)系為:

    式中:[X,Y,Z]是點W在原坐標(biāo)系O-XYZ中的坐標(biāo),[X1,Y1,Z1]是點W在旋轉(zhuǎn)后坐標(biāo)系O-X1Y1Z1中的坐標(biāo)。

    由于-軸為元件的光軸方向,因此繞-軸的旋轉(zhuǎn)并不影響方向向量最后的朝向,所以由式(2)可以看到,方程組是一個與無關(guān)的等式。

    如果要實現(xiàn)任意方向傾斜采樣的蒙特卡羅分析,需要建立光學(xué)元件旋轉(zhuǎn)過程,和空間狀態(tài)量的,之間的轉(zhuǎn)換關(guān)系。因此聯(lián)立式(2)與式(3),得到了,與之間的對應(yīng)關(guān)系:

    這樣,通過旋轉(zhuǎn)矩陣的推導(dǎo),解出了光學(xué)元件傾斜狀態(tài)量(方向角、傾斜角)到傾斜建模過程量(旋轉(zhuǎn)角)的轉(zhuǎn)換對應(yīng)公式。在任意方向蒙特卡羅分析中,只需對方向角和傾斜角進行采樣,就可利用上述公式將采樣數(shù)據(jù)轉(zhuǎn)換建模需要的旋轉(zhuǎn)角,,從而構(gòu)建出元件傾斜的模型,進行非順序光線追跡。

    2 光軸靜態(tài)敏感度分析流程構(gòu)建

    2.1 基本思路

    對紅外折轉(zhuǎn)光學(xué)系統(tǒng)光軸靜態(tài)敏感度分析,本質(zhì)上是公差分析,需將光學(xué)元件的偏心、傾斜公差引入光學(xué)系統(tǒng)中,分析由元件偏移引起的光軸漂移與光軸穩(wěn)定性指標(biāo)的匹配情況。構(gòu)建紅外折轉(zhuǎn)光學(xué)系統(tǒng)的光軸敏感度公差分析流程的基本思路如圖3所示,輸入為光學(xué)系統(tǒng)的仿真模型以及系統(tǒng)的光軸穩(wěn)定性指標(biāo),經(jīng)過光軸敏感度分析以后,最終輸出的是光學(xué)系統(tǒng)中各個元件合理的公差限數(shù)據(jù)。該思路對光軸敏感度分析分解為三步:靈敏度分析-反靈敏度分析-蒙特卡羅分析[21],每一步的分析需要基于上一步分析的結(jié)果數(shù)據(jù)作為輸入,而輸出的數(shù)據(jù)將作為下一步分析的輸入。

    流程中首先開展靈敏度公差分析,根據(jù)光軸穩(wěn)定性指標(biāo)進行初步分解,將光學(xué)系統(tǒng)的公差大致設(shè)定在某個范圍,用非順序建模實現(xiàn)對光學(xué)元件空間姿態(tài)的精確模擬,并通過非順序光線追跡,仿真計算這些公差所造成的光軸漂移量數(shù)據(jù)(如圖4所示),此數(shù)據(jù)即為各元件偏心和傾斜的光軸靈敏度數(shù)據(jù),以此識別敏感元件。

    圖3 構(gòu)建紅外光學(xué)系統(tǒng)光軸敏感度公差分析流程的基本思路

    圖4 透鏡傾斜(a)和透鏡偏心(b)引起的光軸漂移

    反靈敏度分析則假設(shè)光學(xué)系統(tǒng)的評價標(biāo)準(zhǔn)(即光軸漂移量)變化特定的數(shù)值,評估造成相應(yīng)變化的各項公差范圍,稱之為光軸反靈敏度數(shù)據(jù),可根據(jù)上一步光軸靈敏度數(shù)據(jù)的多次計算迭代求得。

    光軸的反靈敏度數(shù)據(jù),可作為蒙特卡羅公差分析中的輸入,即各光學(xué)件偏心和傾斜的光差限初值,以公差限初值作為反射鏡各運動分量的蒙特卡洛采樣區(qū)間,在區(qū)間內(nèi)進行各元件偏心和傾斜隨機取值并疊加成組合誤差,進行了蒙特卡洛分析,得出光軸漂移量的初步統(tǒng)計評估結(jié)果,與系統(tǒng)的光軸穩(wěn)定性的指標(biāo)要求對比后,放寬或收緊采樣區(qū)間,進行多輪分析,最終得到合理的公差限。

    2.2 程序設(shè)計

    根據(jù)上述思路,對光學(xué)系統(tǒng)進行光軸敏感度分析的程序流程如圖5所示,基于Matlab編制計算程序,通過調(diào)用LightTools軟件實現(xiàn)光學(xué)元件偏移及非順序光線追跡,從而完成光軸漂移仿真數(shù)據(jù)計算。該流程以光學(xué)系統(tǒng)的光軸穩(wěn)定性指標(biāo)axis作為流程輸入,以光學(xué)系統(tǒng)中各鏡片的偏心和傾斜合理公差限為輸出。

    程序中,調(diào)取各透鏡的原始偏心值0和原始傾斜值0,分別疊加標(biāo)準(zhǔn)偏心值stan和標(biāo)準(zhǔn)傾斜值stan,通過對軸上主光線追跡計算光軸漂移,并進行公差分析,得到光學(xué)元件的靈敏度數(shù)據(jù)?;诠鈱W(xué)元件的靈敏度數(shù)據(jù)和系統(tǒng)的光軸穩(wěn)定性指標(biāo)axis,計算出光學(xué)元件的反靈敏度數(shù)據(jù)?;诠鈱W(xué)元件的光軸反靈敏度數(shù)據(jù),得到公差限初值0,即光學(xué)元件各運動分量的蒙特卡洛采樣區(qū)間。給定循環(huán)次數(shù)loop,在公差限初值范圍內(nèi)對每一光學(xué)件的偏心方向,偏心量,旋轉(zhuǎn)方向,旋轉(zhuǎn)量四個值分別產(chǎn)生loop組隨機數(shù)。偏心方向和旋轉(zhuǎn)方向在0~360°的范圍內(nèi)采樣,偏心量和旋轉(zhuǎn)量則在公差限初值的范圍內(nèi)采樣。由偏心方向,偏心量,旋轉(zhuǎn)方向,旋轉(zhuǎn)量四個值,通過上節(jié)中的數(shù)據(jù)變換,計算得出光學(xué)元件建模所需的過程量輸入,即沿-軸方向的平移量,沿方向的平移量,繞-軸的旋轉(zhuǎn)量,繞-軸的旋轉(zhuǎn)量。

    ,與,的對應(yīng)關(guān)系可由簡單的三角函數(shù)求得:

    式中:u,v與a,b的對應(yīng)關(guān)系由上節(jié)中推導(dǎo)得出的光學(xué)元件旋轉(zhuǎn)過程量和傾斜狀態(tài)量的轉(zhuǎn)換關(guān)系式(4)求得。

    基于計算得出的位置與旋轉(zhuǎn)角度輸入,進行光學(xué)元件偏移狀態(tài)的建模和光軸漂移的仿真分析,得到各光學(xué)元件偏心與傾斜隨機采樣數(shù)據(jù)搭配后的光線輸出數(shù)據(jù),對輸出數(shù)據(jù)進行統(tǒng)計分析,得出光軸漂移量的蒙特卡洛分析結(jié)果。

    將輸出的蒙特卡羅分析結(jié)果與光軸漂移量的指標(biāo)要求進行對比,得出新一輪蒙特卡羅分析的公差限調(diào)整系數(shù)P+1,如下式計算:

    P+1=axis/I(6)

    式中:I為本輪光軸漂移量蒙特卡洛分析結(jié)果。若P+1已接近1,則說明輸入的公差限已經(jīng)合理,無需再調(diào)整。否則,則需按照下式放寬或收緊公差限:

    +1=P+1×(7)

    式中:為本輪的公差限;+1為新一輪公差限。獲取新一輪公差限后,再次展開蒙特卡羅分析,重復(fù)迭代此過程,直至公差限調(diào)整系數(shù)接近1,即光軸漂移量的蒙特卡羅統(tǒng)計結(jié)果與系統(tǒng)光軸穩(wěn)定性指標(biāo)要求接近,則可將該輪的公差限作為光學(xué)元件偏移的約束值。

    3 光軸敏感度分析實例

    根據(jù)上述搭建的光軸敏感度分析流程和計算程序,選取了光軸較為敏感的某型紅外熱成像光學(xué)系統(tǒng)作為分析對象,光學(xué)系統(tǒng)的仿真模型如圖6所示。該系統(tǒng)為制冷型中波熱像儀,為減少體積,采用了兩個折轉(zhuǎn)反射鏡形成U型折轉(zhuǎn)光路。本例以系統(tǒng)的光軸穩(wěn)定性指標(biāo)axis≤0.25mrad為例,對系統(tǒng)中8個光學(xué)透鏡和2個折轉(zhuǎn)反射鏡的傾斜和偏心的公差限進行分析計算。

    圖6 光學(xué)系統(tǒng)的靜態(tài)仿真模型

    在傾斜和偏心的靈敏度分析中,本例的標(biāo)準(zhǔn)偏心值stan設(shè)置為0.01mm,標(biāo)準(zhǔn)傾斜值stan設(shè)置為0.01°。在傾斜和偏心的反靈敏度分析中,由于系統(tǒng)的光軸穩(wěn)定性指標(biāo)axis≤0.25mrad,而光學(xué)系統(tǒng)中包括反射鏡在內(nèi),共有10個光學(xué)件,且每個光學(xué)件包含兩個公差(偏心和傾斜),這樣一來,共由20個公差對光軸漂移產(chǎn)生貢獻,因此,將0.25mrad的指標(biāo)均分成20份,以光軸漂移量為0.0125mrad作為反靈敏度分析的標(biāo)準(zhǔn)。

    圖7(a)為傾斜靈敏度分析結(jié)果,橫坐標(biāo)是光學(xué)元件的序號,縱坐標(biāo)是光軸漂移量;圖7(b)為傾斜反靈敏度分析結(jié)果,橫坐標(biāo)是光學(xué)元件的序號,縱坐標(biāo)是光學(xué)件的傾斜量。從傾斜靈敏度折線圖中可以較為直觀地看出,當(dāng)每個透鏡同時傾斜0.01°,對光軸漂移量產(chǎn)生的貢獻存在很大差異。對比傾斜靈敏度折線圖與傾斜反靈敏度折線圖中可以發(fā)現(xiàn),傾斜反靈敏度的數(shù)據(jù)與傾斜靈敏度數(shù)據(jù)正好相反,第2反射鏡的傾斜靈敏度最高,傾斜反靈敏度最低,第5透鏡的傾斜靈敏度最低,傾斜反靈敏度最高。從圖中能夠識別出,第1反射鏡和第2反射鏡是對傾斜最敏感的關(guān)鍵元件。

    圖8(a)為偏心靈敏度分析結(jié)果,橫坐標(biāo)是光學(xué)件的序號,縱坐標(biāo)是光軸漂移量;圖8(b)為偏心反靈敏度分析結(jié)果,橫坐標(biāo)是光學(xué)件的序號,縱坐標(biāo)是光學(xué)件的偏心量。對比偏心靈敏度折線圖與偏心反靈敏度折線圖中可以發(fā)現(xiàn),偏心反靈敏度的數(shù)據(jù)與偏心靈敏度數(shù)據(jù)同樣正好相反,第2透鏡偏心靈敏度最高,偏心反靈敏度最低,第5透鏡的偏心靈敏度最低,偏心反靈敏度最高。從圖中能夠識別出,第2透鏡和第3透鏡是對偏心最敏感的關(guān)鍵元件。

    圖7 傾斜靈敏度(a)和傾斜反靈敏度(b)折線圖

    圖8 偏心靈敏度(a)和偏心反靈敏度(b)折線圖

    將光軸反靈敏度分析數(shù)據(jù)設(shè)置為各光學(xué)元件偏心和傾斜的初始公差限,在8個透鏡與2個反射鏡的偏心與傾斜初始公差限范圍內(nèi)進行蒙特卡羅分析,本例中循環(huán)次數(shù)loop設(shè)置為10000次。圖9(a)為蒙特卡羅分析后輸出的10000組光學(xué)系統(tǒng)光軸漂移量數(shù)據(jù)的分布區(qū)間直方圖,橫坐標(biāo)為光軸漂移量,縱坐標(biāo)為分布在相應(yīng)光軸漂移量區(qū)間的蒙特卡羅實驗次數(shù),由圖可知,分布在光軸漂移量為0.021mrad附近區(qū)間的實驗次數(shù)最多。圖9(b)為數(shù)據(jù)的累積分布函數(shù)圖,橫坐標(biāo)為光軸漂移量,縱坐標(biāo)為累計分布概率,可以看到,50%的概率光軸移動量在0.025mrad范圍內(nèi),70%的概率光軸移動量在0.033mrad范圍內(nèi),90%的概率光軸移動量在0.0456mrad范圍內(nèi)。

    本文中以90%的概率為例,因此基于初始公差限,光學(xué)系統(tǒng)的光軸移動量評估結(jié)果大概率分布在0.0456mrad范圍內(nèi),而對比本例設(shè)置的指標(biāo),系統(tǒng)光軸穩(wěn)定性axis≤0.25mrad。因此基于蒙特卡羅分析評估的光軸漂移量結(jié)果還遠低于光軸一致性指標(biāo),說明設(shè)置的初始公差限過于保守,因此按照式(6)計算,公差調(diào)整系數(shù)等于5.48,并按照式(7)放寬各光學(xué)件的偏心與傾斜公差限。

    根據(jù)調(diào)整后的公差限,進行蒙特卡羅隨機采樣。新一輪蒙特卡羅分析的輸出數(shù)據(jù)分布區(qū)間直方圖及累積分布函數(shù)圖如圖10(a)和圖10(b)所示。

    本輪輸出的光軸漂移量數(shù)據(jù)以90%的概率分布在0~0.256mrad范圍內(nèi),基本與光學(xué)系統(tǒng)的光軸一致性指標(biāo)相匹配。因此將本輪輸入的公差限作為最終各光學(xué)件的傾斜和偏心合理公差。

    表1為采用本文程序?qū)δ繕?biāo)光學(xué)系統(tǒng)依次進行了靈敏度分析、反靈敏度分析、蒙特卡羅分析之后,得出的各個光學(xué)元件的合理公差限。此公差限可作為結(jié)構(gòu)優(yōu)化設(shè)計的約束條件,只要由結(jié)構(gòu)件引起的光學(xué)元件偏心和傾斜不超過表1數(shù)據(jù)的范圍,即可保證光學(xué)系統(tǒng)的光軸漂移量不超出系統(tǒng)光軸穩(wěn)定性0.25mrad的指標(biāo)要求。從表1數(shù)據(jù)和前面靈敏度、反靈敏度分析結(jié)果可以看出,結(jié)構(gòu)優(yōu)化設(shè)計需要著重考慮高低溫或外力的作用下,第1、2反射鏡的傾斜、以及第2、3透鏡的偏心。因此,該分析識別出了光學(xué)系統(tǒng)的敏感點,同時能夠指導(dǎo)結(jié)構(gòu)優(yōu)化設(shè)計,以滿足光軸穩(wěn)定性指標(biāo)。

    圖9 光軸漂移量分布直方圖(a)和累積分布圖(b)

    圖10 調(diào)整公差限后的光軸漂移量分布直方圖(a)和光軸漂移量累積分布圖(b)

    表1 各光學(xué)件傾斜與偏心最終公差限

    4 數(shù)據(jù)準(zhǔn)確度驗證

    本文為實現(xiàn)任意角度傾斜的蒙特卡羅分析,采用數(shù)據(jù)轉(zhuǎn)換的思路,利用坐標(biāo)變換推導(dǎo)出了方向角和傾斜角到旋轉(zhuǎn)建模角度的轉(zhuǎn)換關(guān)系式。為了對文中角度轉(zhuǎn)換公式的準(zhǔn)確性進行驗證,利用上一章實例中輸出的公差限,并通過上文中提到的通過手動建立多重臨時坐標(biāo)系的替代方法,對實例中光學(xué)系統(tǒng)再次進行蒙特卡洛分析,觀察是否能得出一致的光軸移動量分布數(shù)據(jù)。

    在仿真軟件中,如圖11所示,根據(jù)每一個光學(xué)元件的原始坐標(biāo)系分別建立一個相同的臨時坐標(biāo)系,并使臨時坐標(biāo)系繞¢-軸左旋90°旋轉(zhuǎn),使臨時坐標(biāo)系的¢-軸朝向原內(nèi)層坐標(biāo)系的-軸,臨時坐標(biāo)系的¢-軸朝向原內(nèi)層坐標(biāo)系的-軸,臨時坐標(biāo)系的¢-軸朝向原內(nèi)層坐標(biāo)系的-軸負方向,將原始坐標(biāo)系嵌套在臨時坐標(biāo)系的內(nèi)層并綁定,讓外層坐標(biāo)系帶動光學(xué)元件一同旋轉(zhuǎn),實現(xiàn)讓光學(xué)件先繞光軸旋轉(zhuǎn),在繞與光軸垂直方向的坐標(biāo)軸旋轉(zhuǎn)。

    以實例中基于最后一輪公差限產(chǎn)生隨機采樣數(shù)據(jù),輸入到臨時坐標(biāo)系中,進行蒙特卡洛分析。得到的數(shù)據(jù)分布區(qū)間直方圖如圖12(a),與上一章中得出的數(shù)據(jù)分布區(qū)間直方圖如圖12(b)進行對比,可以看出,兩種方法得出的結(jié)果分布完全一致,證明了數(shù)據(jù)轉(zhuǎn)換方法的準(zhǔn)確性。

    此外,相較于手動建立多重臨時坐標(biāo)系的方法,本文建立的方法適用性更強、效率更高,省去了很多繁瑣過程,不需要為每一個新的光學(xué)系統(tǒng)的所有光學(xué)元件手動建立與其坐標(biāo)相對應(yīng)的兩重臨時坐標(biāo)系,避免了數(shù)據(jù)輸入或匹配錯誤。

    圖11 臨時坐標(biāo)系示意圖

    圖12 數(shù)據(jù)轉(zhuǎn)換方法得出(a)和臨時坐標(biāo)系方法得出(b)的分布直方圖

    5 結(jié)論

    在對紅外熱成像折轉(zhuǎn)光學(xué)系統(tǒng)進行光軸穩(wěn)定性的結(jié)構(gòu)優(yōu)化過程中,非常有必要針對光學(xué)系統(tǒng)開展光軸的靜態(tài)敏感度分析,基于系統(tǒng)的光軸穩(wěn)定性指標(biāo)分析出各個光學(xué)元件偏心和傾斜的合理公差限,從而指導(dǎo)結(jié)構(gòu)優(yōu)化。為解決蒙特卡羅分析中對折轉(zhuǎn)光學(xué)系統(tǒng)元件的任意方向傾斜采樣的難點,通過基于旋轉(zhuǎn)矩陣的坐標(biāo)變換,建立了光學(xué)元件旋轉(zhuǎn)過程量,和空間狀態(tài)量的,相互轉(zhuǎn)換的方法,使蒙特卡羅分析能夠模擬真實情況下的光學(xué)元件運動,計算效率高且具有通用性。在建立的數(shù)據(jù)轉(zhuǎn)換方法的基礎(chǔ)上,搭建了對紅外光學(xué)系統(tǒng)光軸靜態(tài)敏感度分析的流程,并運用Matlab對LightTools進行二次開發(fā),實現(xiàn)了程序化。與此同時,以一個光軸較為敏感的實際光學(xué)系統(tǒng)為例,運用建立的方法對其開展了光軸靜態(tài)敏感度分析,得出了光學(xué)系統(tǒng)的敏感點和各光學(xué)元件的合理公差限,并結(jié)合建立臨時外重坐標(biāo)系的方法,對數(shù)據(jù)準(zhǔn)確性進行了驗證,證實了公式推導(dǎo)以及程序的可靠性和準(zhǔn)確性。

    文中所建立的光軸靜態(tài)敏感度分析流程為實現(xiàn)完整的光機熱優(yōu)化設(shè)計奠定了基礎(chǔ),它雖然是針對光學(xué)系統(tǒng)的光軸分析所搭建的,但與此同時,文中提出的任意方向傾斜的蒙特卡羅采樣方法也為折轉(zhuǎn)光學(xué)系統(tǒng)像質(zhì)的公差分析提供了解決思路。后續(xù)將對如何優(yōu)化文中方法分析出的敏感元件開展研究,以期進一步提高光學(xué)系統(tǒng)的光軸穩(wěn)定性。

    [1] 郭夏銳, 廖志波, 王春雨, 等. 光軸一致性誤差對空間透射式系統(tǒng)像差和質(zhì)量的影響[J]. 紅外與激光工程, 2012, 41(2): 437-441. GUO Xiarui, LIAO Zhibo, WANG Chunyu, et al. Aberration and alignment quality investigation of space transmission optical system with centering error of optical axes[J]., 2012, 41(2): 437-441.

    [2] 馬宏川, 范宏波, 林宇, 等. 熱像儀光機熱集成分析綜述[J]. 紅外技術(shù), 2019, 41(2): 134-141. MA Hongchuan, FAN Hongbo, LIN Yu, et al. Review of thermal-structural-optical integration analysis of thermal imager[J]., 2019, 41(2): 134-141.

    [3] 王陽, 黃煜, 李占峰, 等. 利用恒星對天文觀測系統(tǒng)光軸平行性檢校[J].紅外與激光工程, 2017, 46(5): 0517003. WANG Yang, Huang Yu, Li Zhanfeng, et al. Calibration of optical axis parallelism by using star for astronomical observation system[J]., 2017, 46(5): 0517003.

    [4] 陳志斌, 肖文健, 馬東璽, 等. 大間距多光軸一致性野外在線檢測方法[J]. 光學(xué)學(xué)報, 2017, 37(1): 0112006.Chen Zhibin, Xiao Wenjian, Ma Dongxi, et al. A method for large distance multi-optical axis parallelism online detection[J]., 2017, 37(1): 0112006.

    [5] 楊耀山, 趙泉, 王建龍, 等. 紅外光學(xué)系統(tǒng)光軸穩(wěn)定性建模與分析[J]. 電光與控制, 2017, 24(12): 82-84.Yang Yaoshan, Zhao Quan, Wang Jianlong, et al. Optical axis stability modeling and analysis of infrared optical system[J]., 2017, 24(12): 82-84.

    [6] 隋杰, 程會艷, 余成武, 等. 星敏感器光軸熱穩(wěn)定性仿真分析方法[J]. 空間控制技術(shù)與應(yīng)用, 2017, 43(4): 38-41.Sui Jie, Cheng Huiyan, Yu Chengwu, et al. A thermal stability analysis and simulation method for boresight axis of star sensor[J]., 2017, 43(4): 38-41.

    [7] 董雪巖, 許方宇, 陳驥, 等. 1m紅外太陽望遠鏡鏡面位置誤差致光軸偏移分析[J]. 紅外技術(shù), 2016, 38(10): 870-876. DONG Xueyan, XU Fangyu, CHEN Ji, et al. Analysis on optic axis deviation induced by lens position disorder of 1m new vacuum solar telescope[J]., 2016, 38(10): 870-876.

    [8] Kalikivayi V, Kumar V C P, Kannan K, et al. Tolerance analysis of misalignment in an optical system using Shack–Hartmann wavefront sensor: experimental study[J]., 2015, 54(7): 075104.

    [9] 高磊. 基于蒙特卡洛法的公差分析及優(yōu)化設(shè)計方法研究[D]. 哈爾濱: 哈爾濱理工大學(xué), 2015.Gao Lei. Research on the Tolerance Analysis Based on Monte Carlo Method and Optimal Design Method[D]. Harbin: Harbin University of Science and Technology, 2015.

    [10] 劉琳, 張興德, 賀誼亮. 基于蒙特卡洛模擬法的紅外光學(xué)系統(tǒng)公差分析[J]. 激光與紅外, 2020, 40(5): 497-499.Liu Lin, Zhang Xingde, He Yiliang. Monte Carlo simulation and its application in the IR optical system[J]., 2020, 40(5): 497-499.

    [11] 袁貴星, 王平. 蒙特卡洛模擬及其在公差設(shè)計中的應(yīng)用[J]. 天津科技大學(xué)學(xué)報, 2008, 23(2): 61-64.Yuan Guixing, Wang Ping. Monte Carlo simulation and its application in tolerance design[J]., 2008, 23(2): 61-64.

    [12] Yang K Y, Jin N, Yang D, et al. Monte Carlo calculation on stray radiation scattered by the baffle with a heterotypic surface for uncooled infrared system[C]//, 2018, 10815: 108150R.

    [13] 蘇朋. 365nm光刻照明系統(tǒng)中變焦系統(tǒng)的設(shè)計及公差分析[J]. 紅外與激光工程, 2022, 51(7): 20210524.Su Peng. Design and tolerance analysis of the zoom system in 365nm UV lithography illumination system[J]., 2022, 51(7): 20210524.

    [14] 王川川, 岑兆豐, 李曉彤, 等. 基于線性模型的聚焦光學(xué)系統(tǒng)靜態(tài)公差優(yōu)化算法[J]. 光子學(xué)報, 2017, 46(11): 1122001. WANG Chuanchuan, CEN Zhaofeng, LI Xiaotong, et al. Optimization algorithm of static tolerance in focusing system based on linear model[J]., 2017, 46(11): 1122001.

    [15] 蔡東浩, 李雅燦, 魏立冬, 等. 高光譜成像儀中曲面棱鏡的裝調(diào)公差研究[J]. 光學(xué)學(xué)報, 2021, 41(6): 60-68.Cai Donghao, Li Yacan, Wei Lidong, et al. Fabricating tolerance of curved prisms in hyperspectral spectrometer[J]., 2021, 41(6): 60-68.

    [16] 李美萱, 董連和. 浸沒式光刻照明系統(tǒng)中會聚鏡的設(shè)計及公差分析[J]. 應(yīng)用光學(xué), 2019, 40(5): 876-881. LI Meixuan, DONG Lianhe. Design and tolerance analysis of converging lens in immersion Lithography illumination system[J]., 2019, 40(5): 876-881.

    [17] 曹宇澤, 馬文禮. 兩步式靈敏度矩陣法在卡塞格林望遠鏡裝調(diào)中的應(yīng)用[J]. 光電工程, 2020, 47(2): 5-12.Cao Yuze, Ma Wenli. Application of two step sensitivity matrix method in Cassegrain telescope alignment[J]., 2020, 47(2): 5-12.

    [18] Synopsys Inc.9.1[M]. Mountain View: Synopsys, 2020.

    [19] ZEMAX LLC.21.3Version September 2021[M]//Kirkland: ZEMAX, 2021.

    [20] 王增偉, 趙知誠, 楊溢, 等. 基于剛體運動完備方程的光機熱集成分析方法[J]. 紅外與激光工程, 2022, 51(6): 20210617.Wang Zengwei, Zhao Zhicheng, Yang Yi, et al. Thermal -structural-optical integrated analysis method based on the complete equations of rigid body motion[J]., 2022, 51(6): 20210617.

    [21] 張坤, 鐘興, 孟遙, 等. 大視場納型星敏光學(xué)系統(tǒng)公差靈敏度研究[J]. 光子學(xué)報, 2019, 48(5): 0522001. ZHANG Kun, ZHONG Xing, MENG Yao, et al. Tolerance sensitivity research of nano-star sensor optical system with large field[J]., 2019, 48(5): 0522001.

    Optical Axis Static Sensitivity Analysis for Infrared Thermal Imaging Folding Optical System

    YANG Dan1,JIN Ning1,YANG Kaiyu1,LI Jing1,DONG Shulin1,HU Jianchuan1,LI Xiaojun2

    (1. Kunming Institute of Physics, Kunming 650223, China;2. Military Representative Office of Chengdu Bureau of AIR Force Stationed in Kunming Area, Kunming 650223, China)

    The optical axis of an infrared thermal imaging folding optical system is prone to shift owing to decentering of the tilt of optical components under complex environmental conditions, which affects the indication accuracy of the system for the target. Static sensitivity analysis of the optical axis for the optical system at the beginning of the design of the infrared thermal imaging system is useful for identifying the sensitive points of the optical system and provides constraints for the structural optimization design to meet the stability of the optical axis. The conversion relationship between the rotation process and spatial state quantities of the optical components was established by coordinate transformation based on the rotation matrix to simulate the spatial attitude of the optical component tilted in any direction and ensure that the Monte Carlo sampling in the optical axis sensitivity analysis corresponds to the constraint conditions of the structural design. On this basis, the flow of the static sensitivity analysis of the optical axis of the infrared folding optical system was established, and a program was compiled. A typical infrared thermal imaging folding optical system was analyzed using this program. According to the index requirements of the optical axis stability, the optical axis sensitivity and inverse sensitivity of the decenter and tilt of each optical component in the optical system were analyzed, and the initial tolerance limit was obtained. Then, Monte Carlo analysis sampling could be performed in any direction according to the initial tolerance limit data; thus, the decenter and tilt tolerance limit data that meet the optical axis stability index could be obtained, and the accuracy of the obtained data was verified by establishing a multi-coordinate system. Static sensitivity analysis provides a foundation for guiding the design of optical–mechanical thermal optimization.

    infrared folding optical system, optical axis shift, static sensitivity, tolerance analysis, Monte Carlo sampling

    O435

    A

    1001-8891(2023)08-0828-09

    2022-12-27;

    2023-02-15.

    楊丹(1991-),女,云南昆明人,工程師,碩士,主要從事光學(xué)系統(tǒng)設(shè)計和仿真方面的研究。E-mail:naomi626391@126.com。

    楊開宇(1984-),女,云南昆明人,正高級工程師,博士,主要從事光學(xué)系統(tǒng)設(shè)計和仿真方面的研究。E-mail:yky20030634@126.com。

    國家重點研發(fā)計劃(2017YFA0701200)。

    猜你喜歡
    光軸蒙特卡羅公差
    非平行光軸三維立體檢測技術(shù)在毫山峪隧道中的應(yīng)用
    光軸的形態(tài)和應(yīng)用
    攝影之友(2020年12期)2020-12-23 09:39:39
    利用蒙特卡羅方法求解二重積分
    智富時代(2019年6期)2019-07-24 10:33:16
    帶有中心值的百分?jǐn)?shù)的公差表示方式
    公差里的“五朵金花”
    帶有中心值的百分?jǐn)?shù)的公差表示方式
    探討蒙特卡羅方法在解微分方程邊值問題中的應(yīng)用
    多視場電視觀瞄具的光軸調(diào)校技術(shù)
    復(fù)合型種子源125I-103Pd劑量場分布的蒙特卡羅模擬與實驗測定
    同位素(2014年2期)2014-04-16 04:57:20
    壓制觀瞄系統(tǒng)多光譜光軸平行性調(diào)校技術(shù)研究
    大陆偷拍与自拍| 在线观看免费午夜福利视频| 成人18禁高潮啪啪吃奶动态图| 无限看片的www在线观看| 成人特级黄色片久久久久久久 | 亚洲伊人色综图| 久久av网站| 欧美在线一区亚洲| 首页视频小说图片口味搜索| 亚洲欧美日韩高清在线视频 | 免费人妻精品一区二区三区视频| 大片电影免费在线观看免费| 一本综合久久免费| 国产无遮挡羞羞视频在线观看| 精品一品国产午夜福利视频| 欧美在线一区亚洲| 黑人猛操日本美女一级片| 欧美日韩国产mv在线观看视频| 欧美成人午夜精品| 一区福利在线观看| 久久人妻熟女aⅴ| 两性夫妻黄色片| 91麻豆av在线| 国产精品久久久久久精品电影小说| 人人妻,人人澡人人爽秒播| 亚洲精品美女久久av网站| 国产成人精品久久二区二区91| 18禁黄网站禁片午夜丰满| 91大片在线观看| 国产精品国产高清国产av | 成人亚洲精品一区在线观看| 国产极品粉嫩免费观看在线| 黄色怎么调成土黄色| 亚洲专区国产一区二区| 国产欧美日韩一区二区三| 亚洲国产精品一区二区三区在线| 欧美日韩一级在线毛片| 91老司机精品| 一区二区三区激情视频| 91麻豆精品激情在线观看国产 | 亚洲第一av免费看| 少妇被粗大的猛进出69影院| av超薄肉色丝袜交足视频| 色婷婷av一区二区三区视频| 亚洲中文日韩欧美视频| 国产欧美日韩一区二区三区在线| av有码第一页| 国产欧美日韩一区二区精品| 精品一区二区三区av网在线观看 | 黄色视频,在线免费观看| 两性夫妻黄色片| www.熟女人妻精品国产| 亚洲熟女毛片儿| 久久精品aⅴ一区二区三区四区| 黄频高清免费视频| 成年女人毛片免费观看观看9 | 操美女的视频在线观看| 丰满迷人的少妇在线观看| av电影中文网址| 黄网站色视频无遮挡免费观看| 男女无遮挡免费网站观看| 一本一本久久a久久精品综合妖精| 性色av乱码一区二区三区2| 嫁个100分男人电影在线观看| 欧美激情极品国产一区二区三区| 日韩视频一区二区在线观看| 免费观看av网站的网址| 日本a在线网址| 国产精品偷伦视频观看了| 日本wwww免费看| 精品午夜福利视频在线观看一区 | 欧美性长视频在线观看| 美女国产高潮福利片在线看| 国产亚洲一区二区精品| 好男人电影高清在线观看| 国产日韩欧美亚洲二区| 久久精品熟女亚洲av麻豆精品| 日韩欧美三级三区| 日本黄色日本黄色录像| 91av网站免费观看| 久久精品91无色码中文字幕| 日本一区二区免费在线视频| 首页视频小说图片口味搜索| 欧美久久黑人一区二区| 亚洲av成人一区二区三| 九色亚洲精品在线播放| 嫩草影视91久久| 国产三级黄色录像| 最新在线观看一区二区三区| 国产成+人综合+亚洲专区| 嫩草影视91久久| 露出奶头的视频| 亚洲国产av影院在线观看| 男男h啪啪无遮挡| 天天躁夜夜躁狠狠躁躁| 国内毛片毛片毛片毛片毛片| 高清av免费在线| 久久久久网色| 国产91精品成人一区二区三区 | 人人妻,人人澡人人爽秒播| 亚洲人成77777在线视频| 午夜老司机福利片| 热99国产精品久久久久久7| 丁香六月天网| 天天操日日干夜夜撸| 亚洲国产中文字幕在线视频| 999久久久精品免费观看国产| 亚洲成人国产一区在线观看| 99久久国产精品久久久| 欧美日韩国产mv在线观看视频| 男女午夜视频在线观看| 少妇裸体淫交视频免费看高清 | 久久精品人人爽人人爽视色| 精品久久蜜臀av无| 国产精品1区2区在线观看. | 99精国产麻豆久久婷婷| 色综合欧美亚洲国产小说| 女性被躁到高潮视频| 新久久久久国产一级毛片| 国内毛片毛片毛片毛片毛片| 欧美另类亚洲清纯唯美| 成年版毛片免费区| 人妻 亚洲 视频| 欧美精品人与动牲交sv欧美| 高潮久久久久久久久久久不卡| 亚洲国产欧美日韩在线播放| 久久精品aⅴ一区二区三区四区| 免费看十八禁软件| 国产一区二区 视频在线| 超碰97精品在线观看| 一二三四社区在线视频社区8| 免费人妻精品一区二区三区视频| 99精品在免费线老司机午夜| www.自偷自拍.com| 国产精品 国内视频| 国产男女内射视频| 麻豆成人av在线观看| 成年动漫av网址| 久久免费观看电影| 两个人看的免费小视频| 999久久久国产精品视频| 国产亚洲欧美在线一区二区| 精品午夜福利视频在线观看一区 | 香蕉丝袜av| av又黄又爽大尺度在线免费看| 老熟妇乱子伦视频在线观看| 日本精品一区二区三区蜜桃| bbb黄色大片| 欧美日韩国产mv在线观看视频| 日韩 欧美 亚洲 中文字幕| 免费看十八禁软件| 在线观看舔阴道视频| 美女午夜性视频免费| 18在线观看网站| 黄片播放在线免费| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品偷伦视频观看了| 免费一级毛片在线播放高清视频 | 日韩成人在线观看一区二区三区| 丝袜人妻中文字幕| 精品国产国语对白av| 国产精品影院久久| 1024视频免费在线观看| 午夜91福利影院| 国产xxxxx性猛交| 亚洲av欧美aⅴ国产| 最新在线观看一区二区三区| 日日夜夜操网爽| 免费黄频网站在线观看国产| 夜夜爽天天搞| 久久国产精品人妻蜜桃| 在线观看免费视频日本深夜| 成年动漫av网址| 香蕉国产在线看| 熟女少妇亚洲综合色aaa.| 欧美 亚洲 国产 日韩一| 啦啦啦免费观看视频1| 国产精品偷伦视频观看了| 中文字幕高清在线视频| 人妻 亚洲 视频| 久久香蕉激情| 成年版毛片免费区| 在线天堂中文资源库| 亚洲成a人片在线一区二区| 9热在线视频观看99| 男女床上黄色一级片免费看| 操美女的视频在线观看| 亚洲中文字幕日韩| 国产精品.久久久| 91麻豆av在线| 午夜两性在线视频| 久久精品国产亚洲av高清一级| 久久久精品94久久精品| av线在线观看网站| 亚洲av日韩在线播放| 最近最新中文字幕大全电影3 | 丁香六月天网| 国产黄频视频在线观看| 在线观看舔阴道视频| 考比视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 日日夜夜操网爽| 女同久久另类99精品国产91| 国产伦理片在线播放av一区| 又大又爽又粗| 国产aⅴ精品一区二区三区波| 亚洲色图av天堂| 国产淫语在线视频| 国产成人精品无人区| 久久久久视频综合| 精品福利观看| 午夜福利欧美成人| 啪啪无遮挡十八禁网站| 人人妻人人澡人人看| 黑人欧美特级aaaaaa片| 热99国产精品久久久久久7| 国产深夜福利视频在线观看| 亚洲男人天堂网一区| 亚洲精华国产精华精| 成年女人毛片免费观看观看9 | 日本wwww免费看| 久久久久久久久久久久大奶| 欧美精品一区二区免费开放| 法律面前人人平等表现在哪些方面| 99久久99久久久精品蜜桃| 激情视频va一区二区三区| 在线观看免费视频网站a站| 丝袜美腿诱惑在线| 老司机午夜十八禁免费视频| 国产精品一区二区免费欧美| 在线观看66精品国产| 欧美精品高潮呻吟av久久| 性高湖久久久久久久久免费观看| 亚洲天堂av无毛| 久久影院123| 国产精品亚洲一级av第二区| 巨乳人妻的诱惑在线观看| 男女免费视频国产| 操美女的视频在线观看| 真人做人爱边吃奶动态| 亚洲中文字幕日韩| 精品国产一区二区三区久久久樱花| 在线观看免费午夜福利视频| 99久久99久久久精品蜜桃| 亚洲国产欧美日韩在线播放| 日韩熟女老妇一区二区性免费视频| 大香蕉久久成人网| 欧美久久黑人一区二区| 九色亚洲精品在线播放| 亚洲视频免费观看视频| 亚洲一区二区三区欧美精品| 国产熟女午夜一区二区三区| 日日夜夜操网爽| 欧美日韩黄片免| 午夜福利欧美成人| 亚洲一区二区三区欧美精品| 精品第一国产精品| 飞空精品影院首页| 美女高潮到喷水免费观看| 欧美亚洲 丝袜 人妻 在线| 国产精品久久电影中文字幕 | 成年人黄色毛片网站| 久久ye,这里只有精品| 国产精品99久久99久久久不卡| 免费看a级黄色片| 国产免费现黄频在线看| 两人在一起打扑克的视频| 日日摸夜夜添夜夜添小说| 欧美亚洲日本最大视频资源| 国产一区二区三区综合在线观看| 亚洲第一欧美日韩一区二区三区 | 精品国产乱码久久久久久男人| 多毛熟女@视频| 亚洲 国产 在线| 国产亚洲精品久久久久5区| 50天的宝宝边吃奶边哭怎么回事| 欧美变态另类bdsm刘玥| 久久久久久亚洲精品国产蜜桃av| 久久久久久久大尺度免费视频| 久久中文字幕一级| 国产高清国产精品国产三级| 一个人免费看片子| 午夜老司机福利片| 欧美亚洲日本最大视频资源| 中文字幕另类日韩欧美亚洲嫩草| 正在播放国产对白刺激| 国产成人精品在线电影| 人人妻人人添人人爽欧美一区卜| 国产精品久久久人人做人人爽| 国产成人一区二区三区免费视频网站| 国产无遮挡羞羞视频在线观看| 狠狠婷婷综合久久久久久88av| 国产一区二区 视频在线| 国产精品麻豆人妻色哟哟久久| 两性夫妻黄色片| 女性被躁到高潮视频| 成人免费观看视频高清| 色视频在线一区二区三区| 国产又色又爽无遮挡免费看| 欧美黄色淫秽网站| 久久精品国产a三级三级三级| 最新的欧美精品一区二区| 热re99久久国产66热| 国产精品久久久人人做人人爽| 最新在线观看一区二区三区| 嫩草影视91久久| 超色免费av| 欧美人与性动交α欧美软件| 精品少妇黑人巨大在线播放| 欧美大码av| 国产精品久久久av美女十八| 亚洲国产av影院在线观看| 丰满饥渴人妻一区二区三| 男女免费视频国产| 午夜91福利影院| 国产欧美日韩精品亚洲av| 欧美中文综合在线视频| 新久久久久国产一级毛片| 中文欧美无线码| 老司机深夜福利视频在线观看| 大香蕉久久网| 黄色怎么调成土黄色| 欧美变态另类bdsm刘玥| 亚洲国产中文字幕在线视频| 亚洲国产毛片av蜜桃av| 欧美日韩视频精品一区| 黑人巨大精品欧美一区二区mp4| 少妇裸体淫交视频免费看高清 | 免费观看av网站的网址| 咕卡用的链子| 国产免费现黄频在线看| 黄色毛片三级朝国网站| 午夜福利,免费看| 国产区一区二久久| 97在线人人人人妻| 精品国产乱码久久久久久小说| 下体分泌物呈黄色| 91av网站免费观看| 午夜福利影视在线免费观看| 国产精品自产拍在线观看55亚洲 | 激情视频va一区二区三区| 在线看a的网站| 大片免费播放器 马上看| 日本精品一区二区三区蜜桃| 久久ye,这里只有精品| 日韩视频在线欧美| 国产单亲对白刺激| 亚洲国产毛片av蜜桃av| 欧美亚洲 丝袜 人妻 在线| 9热在线视频观看99| 狠狠婷婷综合久久久久久88av| 99精品在免费线老司机午夜| 欧美激情 高清一区二区三区| 精品久久蜜臀av无| 国产野战对白在线观看| 热99re8久久精品国产| 天堂俺去俺来也www色官网| 99国产综合亚洲精品| 91老司机精品| 热re99久久精品国产66热6| 午夜福利在线免费观看网站| 亚洲天堂av无毛| 老司机影院毛片| 满18在线观看网站| av又黄又爽大尺度在线免费看| 欧美国产精品va在线观看不卡| 天天影视国产精品| 欧美激情极品国产一区二区三区| 免费看十八禁软件| 中文字幕最新亚洲高清| 午夜精品国产一区二区电影| 欧美日韩亚洲综合一区二区三区_| 大片电影免费在线观看免费| 日本av手机在线免费观看| 极品人妻少妇av视频| 老司机在亚洲福利影院| 亚洲男人天堂网一区| 国产成人一区二区三区免费视频网站| 免费日韩欧美在线观看| 国产真人三级小视频在线观看| 这个男人来自地球电影免费观看| 国产99久久九九免费精品| 免费黄频网站在线观看国产| 这个男人来自地球电影免费观看| 久久国产亚洲av麻豆专区| 他把我摸到了高潮在线观看 | 国产精品欧美亚洲77777| 自线自在国产av| 国产福利在线免费观看视频| 欧美大码av| 天天影视国产精品| 精品福利永久在线观看| 91麻豆av在线| 国产精品免费视频内射| 亚洲精品粉嫩美女一区| 夜夜夜夜夜久久久久| 一级黄色大片毛片| 精品少妇黑人巨大在线播放| 精品久久蜜臀av无| 黄色成人免费大全| 极品教师在线免费播放| 国产精品一区二区免费欧美| 啦啦啦免费观看视频1| 69精品国产乱码久久久| 超碰成人久久| 亚洲精品美女久久久久99蜜臀| 精品欧美一区二区三区在线| 亚洲精品中文字幕在线视频| 亚洲精品中文字幕一二三四区 | 亚洲熟女毛片儿| 老熟女久久久| 9191精品国产免费久久| 久久久国产成人免费| 亚洲国产成人一精品久久久| 久久精品国产亚洲av香蕉五月 | 欧美一级毛片孕妇| 欧美精品一区二区免费开放| 黄色丝袜av网址大全| 久久久久精品国产欧美久久久| 久久精品亚洲精品国产色婷小说| 免费日韩欧美在线观看| 亚洲成国产人片在线观看| 人人妻人人澡人人看| av网站在线播放免费| 亚洲av片天天在线观看| 999精品在线视频| 日韩欧美国产一区二区入口| 久久狼人影院| 国产伦理片在线播放av一区| 色94色欧美一区二区| 大码成人一级视频| 丝袜美足系列| 亚洲精品在线美女| 大香蕉久久网| 黄片大片在线免费观看| 国产精品电影一区二区三区 | 最黄视频免费看| 看免费av毛片| av电影中文网址| www.精华液| 欧美性长视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲欧美激情在线| 亚洲少妇的诱惑av| 国产伦理片在线播放av一区| 少妇被粗大的猛进出69影院| 久久性视频一级片| 香蕉国产在线看| 亚洲综合色网址| 又紧又爽又黄一区二区| 在线观看人妻少妇| 1024视频免费在线观看| 亚洲人成电影免费在线| 91大片在线观看| 少妇粗大呻吟视频| 国产淫语在线视频| a级毛片在线看网站| 桃花免费在线播放| 97人妻天天添夜夜摸| 黑人巨大精品欧美一区二区蜜桃| 久久精品国产亚洲av香蕉五月 | 成人永久免费在线观看视频 | 欧美日本中文国产一区发布| 少妇裸体淫交视频免费看高清 | 黄频高清免费视频| 国产精品 国内视频| 日本av手机在线免费观看| 国产亚洲av高清不卡| 午夜福利乱码中文字幕| 少妇裸体淫交视频免费看高清 | 啪啪无遮挡十八禁网站| 成年人黄色毛片网站| 黑人操中国人逼视频| 国产高清视频在线播放一区| 91av网站免费观看| 香蕉国产在线看| 天堂8中文在线网| 国产精品影院久久| 国产精品免费一区二区三区在线 | 日本一区二区免费在线视频| 中文字幕最新亚洲高清| 不卡av一区二区三区| 久久精品aⅴ一区二区三区四区| 午夜福利在线免费观看网站| 国产欧美日韩一区二区精品| 99久久99久久久精品蜜桃| 国产日韩欧美视频二区| 午夜福利欧美成人| 男女免费视频国产| 欧美日韩中文字幕国产精品一区二区三区 | 午夜激情久久久久久久| 国产在线精品亚洲第一网站| 一区二区三区乱码不卡18| 久久 成人 亚洲| 99久久人妻综合| 亚洲精品美女久久av网站| 国产区一区二久久| 精品少妇内射三级| 亚洲专区中文字幕在线| 午夜福利在线免费观看网站| 2018国产大陆天天弄谢| 香蕉久久夜色| 亚洲精品一二三| 中文字幕人妻熟女乱码| av国产精品久久久久影院| 99九九在线精品视频| 午夜福利一区二区在线看| 亚洲国产精品一区二区三区在线| 性少妇av在线| 亚洲欧美精品综合一区二区三区| 国产精品久久久久久精品古装| 国产亚洲精品第一综合不卡| 免费观看av网站的网址| 首页视频小说图片口味搜索| 夜夜夜夜夜久久久久| 高清在线国产一区| 黑人操中国人逼视频| 99精品欧美一区二区三区四区| 欧美精品人与动牲交sv欧美| 91老司机精品| 国产精品1区2区在线观看. | 少妇裸体淫交视频免费看高清 | 免费看a级黄色片| 亚洲av国产av综合av卡| xxxhd国产人妻xxx| 999精品在线视频| 国产三级黄色录像| 人人妻人人澡人人看| 丰满人妻熟妇乱又伦精品不卡| 国产亚洲一区二区精品| 国产欧美日韩一区二区三| 99香蕉大伊视频| 亚洲中文日韩欧美视频| 日韩欧美一区二区三区在线观看 | 久久久久精品国产欧美久久久| 亚洲va日本ⅴa欧美va伊人久久| 国产一区有黄有色的免费视频| 99香蕉大伊视频| 怎么达到女性高潮| 9191精品国产免费久久| 久热这里只有精品99| av有码第一页| 69av精品久久久久久 | av网站在线播放免费| 国产精品久久电影中文字幕 | 国产精品1区2区在线观看. | 后天国语完整版免费观看| 黄色成人免费大全| 在线观看人妻少妇| 在线观看免费午夜福利视频| 国产一区二区 视频在线| 久久99热这里只频精品6学生| 亚洲精品一二三| 日本撒尿小便嘘嘘汇集6| 丝瓜视频免费看黄片| 国产精品二区激情视频| 欧美激情极品国产一区二区三区| 一进一出抽搐动态| 日韩大片免费观看网站| 亚洲情色 制服丝袜| 老司机靠b影院| 老司机影院毛片| 黄色怎么调成土黄色| 热99久久久久精品小说推荐| 亚洲欧洲日产国产| 国产成人免费无遮挡视频| 欧美激情 高清一区二区三区| 国产精品欧美亚洲77777| 男人舔女人的私密视频| 青草久久国产| 日日爽夜夜爽网站| 国产一区二区三区综合在线观看| 99热网站在线观看| 欧美精品高潮呻吟av久久| 满18在线观看网站| 91九色精品人成在线观看| 久久狼人影院| 自拍欧美九色日韩亚洲蝌蚪91| 麻豆乱淫一区二区| 满18在线观看网站| 窝窝影院91人妻| 十八禁网站网址无遮挡| 亚洲国产av影院在线观看| 欧美日韩亚洲高清精品| 中文字幕另类日韩欧美亚洲嫩草| 日韩中文字幕视频在线看片| 精品亚洲成a人片在线观看| 久久人妻熟女aⅴ| 久久精品国产99精品国产亚洲性色 | 久久久国产欧美日韩av| 欧美日韩视频精品一区| 热re99久久国产66热| 露出奶头的视频| 日本一区二区免费在线视频| 丁香六月欧美| 1024香蕉在线观看| av网站在线播放免费| 国产精品亚洲一级av第二区| 老鸭窝网址在线观看| 女人爽到高潮嗷嗷叫在线视频| 岛国毛片在线播放| 国产一区二区 视频在线| 99国产极品粉嫩在线观看| 黄片小视频在线播放| 国产精品欧美亚洲77777| 99国产极品粉嫩在线观看| aaaaa片日本免费| 国产成人系列免费观看| 中亚洲国语对白在线视频| 女同久久另类99精品国产91| 免费人妻精品一区二区三区视频| 黑人操中国人逼视频| 久久中文字幕人妻熟女| 精品一区二区三区av网在线观看 | 99久久人妻综合|