一種空間飛網裝置參數優化設計系統和優化方法

            文檔序號:6546234閱讀:248來源:國知局
            一種空間飛網裝置參數優化設計系統和優化方法
            【專利摘要】一種空間飛網裝置參數優化設計系統,包括人機交互模塊M1,控制與計算模塊M2,數據庫模塊M3和后處理模塊M4;人機交互模塊M1將輸入的優化模型的參數、運行環境參數和指令信息記錄并傳遞給控制與計算模塊M2,控制與計算模塊M2根據這些信息從數據庫模塊M3選擇優化算法、模型函數,再通過控制與計算模塊M2的優化子模塊更新設計點,控制與計算模塊M2的仿真分析子模塊計算新設計點的仿真結果,通過循環計算獲得最優設計點;后處理模塊M4讀取最優設計和尋優過程信息并進行分析;用戶通過人機交互模塊M1查看后處理模塊M4分析與匯總的結果;一種空間飛網裝置參數優化方法,它有六大步驟。本發明在航天【技術領域】里具有應用前景。
            【專利說明】一種空間飛網裝置參數優化設計系統和優化方法
            【技術領域】
            [0001]本發明涉及一種空間飛網裝置參數優化設計系統和優化方法,該優化系統可以快速地完成空間飛網裝置參數的優化設計,屬于航天【技術領域】。
            技術背景
            [0002]隨著人類對空間資源的利用,空間碎片的數量逐漸增加,對在軌正常運行的航天器構成了嚴重的威脅,空間飛網裝置就是為了清除空間碎片而設計的一類裝置,其目的主要是通過直接拋射展開一張面積較大的柔性網來包覆較大的空間碎片,以實行空間碎片的離軌。空間飛網裝置的捕獲效果是有柔性網的展開效果決定,而柔性網的展開效果的關鍵因素是空間飛網裝置的相關設計參數,例如發射速度與角度。

            【發明內容】

            [0003]1.目的:本發明的目的是為了提供一種空間飛網裝置參數優化設計系統和優化方法,它可以快速地獲得空間飛網裝置最佳的參數組合,實現柔性網的最佳展開效果。
            [0004]2.技術方案:為了實現上述目的,本發明采用以下技術方案
            [0005](I)本發明一種空間飛網裝置參數優化設計系統,Optimization Platform for aSpace Net Capture System V1.0,簡稱 OPSNCS V1.0,包括以下 4 個部分:0PSNCS V1.0 人機交互模塊Ml,OPSNCS V1.0控制與計算模塊M2,OPSNCS V1.0數據庫模塊M3和OPSNCS V1.0后處理模塊M4。這四個模塊之間的關系是:
            [0006]OPSNCS V1.0人機交互模塊Ml將用戶輸入的優化模型的參數、運行環境參數和軟件運行控制指令等信息記錄并傳遞給OPSNCS V1.0控制與計算模塊M2,OPSNCS V1.0控制與計算模塊M2根據這些信息從OPSNCS V1.0數據庫模塊M3選擇優化算法、試驗設計方法、近似模型函數等,再通過OPSNCS V1.0控制與計算模塊M2的優化子模塊更新設計點,0PSNCSV1.0控制與計算模塊M2的仿真分析子模塊計算新設計點的仿真結果,通過這種循環計算獲得最優設計點。OPSNCS V1.0后處理模塊M4讀取最優設計點、最優設計和尋優過程的信息,并進行分析。用戶可通過OPSNCS V1.0人機交互模塊Ml查看OPSNCS V1.0后處理模塊M4分析與匯總的結果。
            [0007]下面詳細說明各部分的結構與功能:
            [0008]所述OPSNCS V1.0人機交互模塊Ml是用戶對軟件進行設置、控制與查看的模塊,包括以下三個子模塊:優化定義模塊M11,運行環境參數設置模塊M12,軟件運行控制工具M13。它們之間的關系是相互并列的。
            [0009]該優化定義模塊Mll是對空間飛網裝置參數優化中設計變量、優化目標和設計約束的定義。
            [0010]在優化定義模塊Mll中,定義了空間飛網裝置中牽引塊、艙蓋和空間網的發射參數等6個設計變量。數學表達如下:
            _] X = {mTm,,mh, vTm, Θ Tm, VJ ;[0012]其中,X為設計變量組,由五個變量組成:mTm為空間飛網裝置中牽引塊的質量,浮點變量,單位為kg ;mh為空間飛網裝置中網艙蓋的質量,浮點變量,單位為kg ;vTm為空間飛網裝置中牽引塊發射速度的大小,浮點變量,單位為m/s; θΤπ為空間飛網裝置中牽引塊發射速度的角度,浮點變量,單位為度;vh為空間飛網裝置中牽引塊發射速度的大小,浮點變量,單位為m/s。各變量的取值范圍需要根據空間飛網裝置的結構、材料、工況以及捕獲任務等設定。
            [0013]在優化定義模塊Mll中,定義了空間飛網裝置參數優化目標,也即通過優化設計變量,使得空間網在指定距離處展開到最大面積,數學表達如下:
            [0014]Max (A)
            [0015]其中,A為A為空間網在指定位置處展開的面積,單位為m2。
            [0016]在優化定義模塊Mll中,定義了設計約束,分別從空間網的強度和展開過程的效果出發,強度滿足空間網的材料性能,以空間網的最大內應力小于材料的許用應力作為約束,展開效果要求空間網在展開過程中不出現空間網面積收縮和空間網纏繞的情況,以空間網在到達指定位置前展開面積變化率大于零作為約束。其數學表達式如下:
            [0017]C= (C1,…,cm, ac” …,ack);
            [0018]其中,C為約束條件,由2部分組成:Cl,-,Cffl為強度約束,分別對應著空間網的各段網繩的最大應力小于網繩材料的許用應力;aCl,…,ack為空間網展開過程的效果約束,對應的是在空間網在到 達指定位置前的各個時刻空間網展開面積變化率應大于零。各約束具體取值與空間網的材料屬性、發射參數等有關。
            [0019]該運行環境參數設置模塊M12是對空間飛網裝置參數優化系統OPSNCS V1.0運行環境進行設置的模塊。包括求解器的選擇、文件存儲路徑的設置等。
            [0020]該軟件運行控制工具M13是對空間飛網裝置參數優化系統OPSNCS V1.0運動過程的檢查、開始、中止、終止過程的控制模塊。
            [0021]所述OPSNCS V1.0控制與計算模塊M2是空間飛網裝置參數優化系統的核心模塊,實現對優化過程中數據的管理與控制。OPSNCS V1.0控制與計算模塊M2包括以下三個子模塊:數據管理模塊M21,優化模塊M22以及仿真分析模塊M23。它們之間的關系是相互并列的。
            [0022]該數據管理模塊M21實現空間飛網裝置參數優化過程數據流的控制,以及訪問空間飛網裝置參數優化系統的數據庫模塊M3。
            [0023]該優化模塊M22是利用數據管理模塊M21選擇的優化算法計算出下一個設計方案的程序模塊。
            [0024]該仿真分析模塊M23是對設計方案進行空間網裝置展開過程的仿真與分析。
            [0025]所述OPSNCS V1.0數據庫模塊M3包括三個子模塊:優化算法庫M31,試驗設計方法庫M32,近似模型函數庫M33。它們之間的關系是相互并列的。
            [0026]該優化算法庫M31是包括常用的全局算法如多島遺傳算法和常用的梯度算法如SQP的算法數據庫。其中SQP算法是序列二次規劃算法的英文縮寫。
            [0027]該試驗設計方法庫M32是包括常用的試驗設計方法的方法數據庫,如正交試驗設計方法。
            [0028]該近似模型函數庫M33是包括數據擬合常用函數的函數數據庫,如多項式響應面函數。
            [0029]所述OPSNCS V1.0后處理模塊M4包括兩個子模塊:結果讀取模塊M41和結果分析模塊M42。它們之間的關系是相互并列的。
            [0030]該結果讀取模塊M41是在優化循環停止后,查詢優化和仿真結果文件并提取最優設計結果的程序模塊。
            [0031]該結果分析模塊M42是根據結果讀取模塊M41讀取到的信息,分析優化結果的可信度的程序模塊。
            [0032](2)本發明一種空間飛網裝置參數優化方法,該方法具體步驟如下:
            [0033]步驟一:空間飛網裝置參數優化模型定義SI ;
            [0034]其中,優化模型SI的定義包括優化目標函數S11,優化約束條件S12,優化設計變量S13,三者之間是并列關系。
            [0035]優化目標函數Sll是指空間飛網裝置參數優化的目標,也就是空間網在指定距離處展開到最大面積,其數學表達式如下:
            [0036]Max ⑷;
            [0037]其中,A為空間網在指定位置處展開的面積,單位為m2。
            [0038]優化約束條件S12是指空間飛網裝置參數優化的約束條件,分別從空間網的強度和展開過程的效果出發,強度滿足空間網的材料性能,以空間網的最大內應力小于材料的許用應力作為約束,展開效果要求空間網在展開過程中不出現空間網面積收縮和空間網纏繞的情況,以空間網在到達指定位置前展開面積變化率大于零作為約束。其數學表達式如下:
            [0039]C = (C1,…,cm, ac” …,ack);
            [0040]其中,C為約束條件,由2部分組成:Cl,-,Cffl為強度約束,分別對應著空間網的各段網繩的最大應力小于網繩材料的許用應力;aCl,…,ack為空間網展開過程的效果約束,對應的是在空間網在到達指定位置前的各個時刻空間網展開面積變化率應大于零。各約束具體取值與空間網的材料屬性、發射參數等有關。
            [0041]優化設計變量S13是指空間飛網裝置參數優化的各個設計變量,包括空間飛網裝置的要進行優化設計的發射速度參數和結構參數,其數學表達式如下:
            [0042]X= (XpX2WXj ;
            [0043]其中,X為設計變量組,由若干個變量組成:Xi(i = 1,…,η)的具體的物理含義根據不同的空間飛網裝置設計任務而不同,例如:Xl為空間飛網裝置中牽引塊的質量,單位為kg ;&為空間飛網裝置中網艙蓋質量,單位為kg ;x3為空間飛網裝置牽引塊的發射速度的大小,單位為m/s ;x4為空間飛網裝置牽引塊的發射角度,單位為度;x5為空間飛網裝置網艙蓋發射速度的大小,單位為m/s ;x6為空間飛網裝置中牽引塊和網艙蓋發射時間差,單位為S,具體使用時設計變量不限以上幾種。
            [0044]步驟二:空間飛網裝置展開過程的仿真分析S2 ;
            [0045]空間飛網裝置展開過程的仿真分析S2是對空間網發射展開過程進行仿真分析和提取優化所需要分析結果,包括六個子步驟:定義任務參數S21,定義結構尺寸S22,定義空間網的材料屬性S23,定義仿真單元屬性S24,定義載荷工況S25,分析求解S26和結果讀取S27。[0046]定義任務參數S21是指空間飛網任務的參數,包括捕獲的目標的最大展長和距離,捕獲目標的最大展長單位為m,取值范圍5m?30m ;捕獲目標的距離單位為m,取值范圍30m ?120mo
            [0047]定義結構尺寸S22是指空間網的形狀、大小、網格密度以及網繩的粗細,例如,正方形網的邊長,單位為m,取值為20m?40m;網格密度是指空間網網格疏密程度,通過網邊單位為個,取值為整數,取值范圍為20?100。
            [0048]定義空間網的材料屬性S23是實現空間飛網裝置材料屬性的定義,包括空間網網繩的材料定義、牽引塊的材料定義、網艙蓋的材料定義,其中材料定義是指指定材料的彈性模量、密度和泊松比等參數,例如:牽引塊和艙蓋可以選擇45號鋼或者其他合金鋼,網繩材料可選尼龍或者其他復合材料等。
            [0049]定義仿真單元屬性S24實現對仿真結構模型單元的定義。包括結構仿真模型各部分的單元定義。例如牽引塊和網艙蓋使用集中質點單元,網繩使用繩索單元。
            [0050]定義分析求解S25實現空間飛網裝置在定義的發射參數下空間網展開過程求解。我們選用求解器對空間網的展開過程求解。其中,求解器選用現有成熟軟件。
            [0051]定義結果讀取S26實現對分析結果進行處理。在空間飛網裝置分析求解完畢后,求解器會輸出從空間網發射開始到空間網到達指定距離所經歷的時間段內,空間網展開面積和繩索單元的應力等隨時間變化的信息。
            [0052]步驟三:空間飛網裝置設計變量的試驗設計S3
            [0053]空間飛網裝置設計變量的試驗設計S3是實現試驗點在求解空間的合理分布的數理統計技術,包括三個子步驟,即選擇試驗設計方法S31,設置試驗設計參數S32和制定試驗設計表格S33。
            [0054]定義選擇試驗設計方法S31是從數理統計學科中各種試驗設計方法中選擇其中之一,以明確試驗設計點在設計空間分布規律。例如,可選擇正交試驗設計,正交試驗設計是一種快速高效的試驗設計方法,是數理統計學科中的一門成熟技術。
            [0055]定義設置試驗設計參數S32是指設置試驗設計的因子數和設計變量的水平數,以確定試驗設計點的數量和分布。其中,因子數是指設計變量的個數,設計變量的水平數是指試驗設計中該設計變量在取值范圍選定試驗值的個數。選定實驗設計方法后,試驗設計點的個數取決于因子數和設計變量的水平數,例如正交試驗設計中,對于η個設計變量,因子數為η,若選定每一個設計變量的水平數為m,則可以選擇常用的正交試驗表格,確定樣本點的分布。樣本點越多,分布越均勻,越能反映出設計變量與優化目標函數、設計變量與設計約束之間的關系,但相應的計算效率也會越低。
            [0056]定義制定試驗設計表格S33是指根據選擇的試驗設計方法和試驗設計參數設置,指定出試驗設計表格。例如因子數為n,每個設計變量的水平數為m的正交試驗表格中包含每個設計變量的具體取值。
            [0057]在完成以上三個子步驟后,將試驗設計表格的設計點進行空間飛網裝置展開過程的仿真分析S2。空間飛網裝置展開過程的仿真分析S2結束后,獲得設計點的仿真結果。
            [0058]步驟四:構建近似數學函數模型S4
            [0059]構建近似數學函數模型S4是指根據試驗設計得到的設計變量與優化目標函數、設計變量與設計約束之間的關系,利用數學函數擬合,獲得可替代空間飛網展開過程分析的近似數學函數模型。構建近似數學函數模型S4包括2個子步驟:選擇構建近似數學函數模型的方法S41 ;擬合數據S42。
            [0060]定義選擇構建近似數學函數模型的類型S41是指選擇數據擬合的方法,實現設計變量與優化目標函數、設計變量與設計約束之間的關系的最佳。
            [0061]定義擬合數據S42是指利用S41選擇的方法,對試驗設計所有設計方案與結果進行擬合,構建出可替代空間飛網展開過程仿真分析的近似數學函數模型。
            [0062]步驟五:近似模型計算分析S5
            [0063]近似模型計算分析S5是指在獲得設計方案后,利用近似數學函數模型給出對應的優化目標函數值和設計約束值。
            [0064]步驟六:優化S6
            [0065]優化S6是實現優化循環與收斂判斷的過程,包括選擇優化路徑S61、選擇優化算法S62、更新設計方案S63和優化收斂判斷S64這四個子步驟。
            [0066]選擇優化路徑S61是指優化過程中從仿真分析和近似數學函數模型兩種中選擇一種路徑獲得設計方法的結果反饋。
            [0067]選擇優化算法S62是指針對選擇的優化路徑,從優化算法中選出一種或多種方法用于優化迭代計算中。
            [0068]更新設計方案S63是指根據選擇的優化算法,計算得出新的設計方案的過程。
            [0069]優化收斂判斷S64是指計算前后兩組設計方案所得到的優化目標函數值的差的絕對值,并判斷如果結果小于或等于指定收斂精度,則停止優化,最后設計方案為最優設計方案,否則更新設計方案后進行空間飛網裝置展開過程的仿真分析S2或者近似模型計算分析S5,獲得相應的結果反饋。其中,收斂精度一般選用0.1 %?10%,收斂精度越高,優化效果越好,但相應地會延長系統工作時間。
            [0070]3.本發明一種空間飛網裝置參數優化設計系統的優點:
            [0071]空間飛網裝置參數優化設計系統通過界面操作,實現便捷的優化過程操作與控制;
            [0072]空間飛網裝置參數優化設計系統可以高效、快捷的實現空間飛網裝置參數的優化設計;
            [0073]空間飛網裝置參數優化設計系統按功能的模塊進行設計,便于軟件升級;
            [0074]空間飛網裝置參數優化設計系統提供擴展的接口,便于實現軟件的拓展。
            [0075]本發明一種空間飛網裝置參數優化方法的優點:
            [0076]可以高效、快速地實現空間飛網裝置參數優化設計;
            [0077]各模塊之間在軟件結構上相互獨立,具有較強的靈活性。
            【專利附圖】

            【附圖說明】
            [0078]圖1為本發明一種空間飛網裝置參數優化設計系統OPSNCS V1.0參數優化設計對象不意圖;
            [0079]圖2為本發明空間飛網裝置參數優化設計系統OPSNCS V1.0的組成模塊框架圖;
            [0080]圖3本發明一種空間飛網裝置參數優化系統的分步結構示意圖;
            [0081]圖4為本發明空間飛網裝置參數優化方法的流程圖。[0082]圖中具體標號說明如下:
            [0083]I—空間飛網裝置牽引塊;2—空間飛網裝置網艙蓋;3—空間網;4—空間網儲箱。
            【具體實施方式】
            [0084]圖1為本發明實施例空間飛網裝置參數優化設計對象的示意圖。
            [0085]圖2為本發明空間飛網裝置參數優化系統OPSNCS V1.0的組成模塊框架圖。
            [0086](I)本發明一種空間飛網裝置參數優化設計系統,Optimization Platform for aSpace Net Capture System V1.0,簡稱 OPSNCS V1.0,包括以下 4 個部分:0PSNCS V1.0 人機交互模塊Ml,OPSNCS V1.0控制與計算模塊M2,OPSNCS V1.0數據庫模塊M3和OPSNCS V1.0后處理模塊M4。這四個模塊之間的關系是:
            [0087]OPSNCS V1.0人機交互模塊Ml將用戶輸入的優化模型的參數、運行環境參數和軟件運行控制指令等信息記錄并傳遞給OPSNCS V1.0控制與計算模塊M2,OPSNCS V1.0控制與計算模塊M2根據這些信息從OPSNCS V1.0數據庫模塊M3選擇優化算法、試驗設計方法、近似模型函數等,再通過OPSNCS V1.0控制與計算模塊M2的優化子模塊更新設計點,OPSNCSV1.0控制與計算模塊M2的仿真分析子模塊計算新設計點的仿真結果,通過這種循環計算獲得最優設計點。OPSNCS V1.0后處理模塊M4讀取最優設計點、最優設計和尋優過程的信息,并進行分析。用戶可通過OPSNCS V1.0人機交互模塊Ml查看OPSNCS V1.0后處理模塊M4分析與匯總的結果。
            [0088]下面詳細說明各部分的結構與功能:
            [0089]所述OPSNCS V1.0人機交互模塊Ml是用戶對軟件進行設置、控制與查看的模塊,包括以下三個子模塊:優化定義模塊M11,運行環境參數設置模塊M12,軟件運行控制工具M13。它們之間的關系式相互并列的。
            [0090]該優化定義模塊Mll是對空間飛網裝置參數優化中設計變量、優化目標和設計約束的定義。
            [0091]在優化定義模塊Mll中,定義了空間飛網裝置中牽引塊、艙蓋和空間網的發射參數等6個設計變量。數學表達如下:
            [0092]X = {mTm,,mh, vTm, Θ Tm, vJ ;
            [0093]其中,X為設計變量組,由六個變量組成:mTm為空間飛網裝置中牽引塊的質量,浮點變量,單位為kg ;mh為空間飛網裝置中網艙蓋的質量,浮點變量,單位為kg ;vTm為空間飛網裝置中牽引塊發射速度的大小,浮點變量,單位為m/s; θΤπ為空間飛網裝置中牽引塊發射速度的角度,浮點變量,單位為度;vh為空間飛網裝置中牽引塊發射速度的大小,浮點變量,單位為m/s; At為空間飛網裝置中牽引塊和網艙蓋發射時間間隔,浮點變量,單位為S。各變量的取值范圍需要根據空間飛網裝置的結構、材料、工況以及捕獲任務等設定。
            [0094]在優化定義模塊Mll中,定義了空間飛網裝置參數優化的目標,也即通過優化設計變量,使得空間網在指定距離處展開到最大面積,數學表達如下:
            [0095]Max (A)
            [0096]其中,A為A為空間網在指定位置處展開的面積,單位為m2。
            [0097]在優化定義模塊Mll中,定義了優化約束,分別從空間網的強度和展開過程的效果出發,強度滿足空間網的材料性能,以空間網的最大內應力小于材料的許用應力作為約束,展開效果要求空間網在展開過程中不出現空間網面積收縮和空間網纏繞的情況,以空間網在到達指定位置前展開面積變化率大于零作為約束。其數學表達式如下:
            [0098]C = (C1,…,cm, ac” …,ack);
            [0099]其中,C為約束條件,由2部分組成:Cl,-,Cffl為強度約束,分別對應著空間網的各段網繩的最大應力小于網繩材料的許用應力;aCl,…,ack為空間網展開過程的效果約束,對應的是在空間網在到達指定位置前的各個時刻空間網展開面積變化率應大于零。各約束具體取值與空間網的材料屬性、發射參數等有關。
            [0100]該運行環境參數設置模塊M12是對空間飛網裝置參數優化系統OPSNCS V1.0運行環境進行設置的模塊。包括求解器的選擇、文件存儲路徑的設置等。如設置程序運行系統為windows7professional,設置安裝目錄為”D: \Program Files\”,設置系統當前路徑為” C:\Users\”。
            [0101]定義軟件運行控制工具M13是對空間飛網裝置參數優化系統OPSNCS V1.0運動過程的檢查、開始、中止、終止過程的控制模塊。如檢查輸入的數據是否合法,初始參數是否在設計空間內,系統運行路徑是否正確;系統運行過程中需要中止或終止以修改參數等。
            [0102]所述OPSNCS V1.0控制與計算模塊M2是空間飛網裝置參數優化系統的核心模塊,實現對優化過程中數據的管理與控制。OPSNCS V1.0控制與計算模塊M2包括以下三個子模塊:數據管理模塊M21,優化模塊M22以及仿真分析模塊M23。
            [0103]該數據管理模塊M21實現空間飛網裝置參數優化過程數據流的控制,以及訪問空間飛網裝置參數優化系統的數據庫模塊M3。
            [0104]該優化模塊M22是利用數據管理模塊M21選擇的優化算法計算出下一個設計方案的程序模塊。
            [0105]該仿真分析模塊M23是對設計方案進行空間網裝置展開過程的仿真與分析
            [0106]所述OPSNCS V1.0數據庫模塊M3包括三個子模塊:優化算法庫M31,試驗設計方法庫M32,近似模型函數庫M33。
            [0107]該優化算法庫M31是包括常用的全局算法如多島遺傳算法和常用的梯度算法如SQP的算法數據庫。其中SQP算法是序列二次規劃算法的英文縮寫。
            [0108]該試驗設計方法庫M32是包括常用的試驗設計方法的方法數據庫,如正交試驗設計方法。
            [0109]該近似模型函數庫M33是包括數據擬合常用函數的函數數據庫,如多項式響應面函數。
            [0110]所述OPSNCS V1.0后處理模塊M4包括兩個子模塊:結果讀取模塊M41和結果分析模塊M42。
            [0111]該結果讀取模塊M41是在優化循環停止后,查詢優化和仿真結果文件并提取最優設計結果的程序模塊。
            [0112]該結果分析模塊M42是根據結果讀取模塊M41讀取到的信息,分析優化結果的可信度的程序模塊。
            [0113]圖3為本發明一種空間飛網裝置參數優化系統的分步結構示意圖。
            [0114](2)本發明一種空間飛網裝置參數優化方法,主要包括以下幾個步驟:[0115]步驟一:S1空間飛網裝置參數優化模型定義。
            [0116]空間飛網裝置參數優化模型定義SI步驟包括Sll優化目標函數,S12優化約束條件,S13優化設計變量,三者之間是并列關系。
            [0117]所述Sll優化目標函數是指空間飛網裝置參數優化的目標,也就是空間網在指定距離處展開到最大面積,其數學表達式如下:
            [0118]Max (A);
            [0119]其中,A為空間網在指定位置處展開的面積,單位為m2。
            [0120]所述S12優化約束條件是指空間飛網裝置參數優化的約束條件,分別從空間網的強度和展開過程的效果出發,強度滿足空間網的材料性能,以空間網的最大內應力小于材料的許用應力作為約束,展開效果要求空間網在展開過程中不出現空間網面積收縮和空間網纏繞的情況,以空間網在到達指定位置前展開面積變化率大于零作為約束。其數學表達式如下:
            [0121]C = (C1,…,cm, acl,…,ack);
            [0122]其中,C為約束條件,由2部分組成:Cl,-,Cffl為強度約束,分別對應著空間網的各段網繩的最大應力小于網繩材料的許用應力,空間網3使用的是聚酰胺材料,拉伸彈性模量約2000Mpa,許用拉應力約為200Mpa,牽引塊I和艙蓋2使用45號鋼,材料彈性模量約為196Gpa,泊松比約為0.3 ;aCl,…,ack為空間網展開過程的效果約束,對應的是在空間網在到達指定位置前的各 個時刻空間網展開面積變化率應大于零。
            [0123]所述S13優化設計變量是指空間飛網裝置參數優化的各個設計變量,主要包括空間飛網裝置的要進行優化設計的發射速度參數和結構參數,其數學表達式如下:
            [0124]X = {x1; X2, X3, X4, X5I ;
            [0125]其中,X為設計變量組,由五個變量組成=X1為空間飛網裝置中牽引塊的質量,浮點變量,單位為kg,取值范圍為0.5~5 ;x2為空間飛網裝置中網艙蓋質量,浮點變量,單位為kg,取值范圍為0.2~4 ;x3為空間飛網裝置牽引塊的發射速度的大小,浮點變量,單位為m/s,取值范圍為10~45 ;x4為空間飛網裝置牽引塊的發射角度,浮點變量,單位為度,取值范圍為10~60 ;x5為空間飛網裝置網艙蓋發射速度的大小,浮點變量,單位為m/s,取值范圍為5~20。
            [0126]步驟二:S2空間飛網裝置展開過程的仿真分析;
            [0127]所述S2空間飛網裝置展開過程的仿真分析是對空間網發射展開過程進行仿真分析和提取優化所需要分析結果,主要包括六個子步驟:S21定義任務參數,S22定義結構尺寸,S23定義空間網的材料屬性,S24定義仿真單元屬性,S25分析求解和S26結果讀取。
            [0128]所述S21任務參數主要是指空間飛網任務的參數,主要包括捕獲的目標的最大展長和距離,捕獲目標的最大展長單位為m,取值范圍5m~30m,可取目標的最大展長為15m ;捕獲目標的距離單位為m,取值范圍30m~120m,可取目標的距離為60m。
            [0129]所述S22結構尺寸主要是指空間網的形狀、大小、網格密度以及網繩的粗細,例如,正方形網的邊長,單位為m,取值為20m~40m ;網格密度是指空間網網格疏密程度,通過網邊單位為個,取值為整數,取值范圍為20~100。
            [0130]所述S23空間網的材料屬性主要是實現空間飛網裝置材料屬性的定義,主要包括空間網網繩的材料定義、牽引塊的材料定義、網艙蓋的材料定義,其中材料定義是指指定材料的彈性模量、密度和泊松比等參數,例如:牽引塊和艙蓋可以選擇45號鋼或者其他合金鋼,網繩材料可選尼龍或者其他復合材料等。
            [0131]所述S24仿真單元屬性主要實現對仿真結構模型單元的定義。主要包括結構仿真模型各部分的單元定義。例如牽引塊和網艙蓋使用集中質點單元,網繩使用繩索單元。
            [0132]分析求解S25主要實現空間飛網裝置在定義的發射參數下空間網展開過程求解。我們選用求解器對空間網的展開過程求解。其中,求解器選用現有成熟軟件。
            [0133]所述S26結果讀取主要實現對分析結果進行處理。在空間飛網裝置分析求解完畢后,求解器會輸出從空間網發射開始到空間網到達指定距離所經歷的時間段內,空間網展開面積和繩索單元的應力等隨時間變化的信息。
            [0134]步驟三:S3空間飛網裝置設計變量的試驗設計
            [0135]S3空間飛網裝置設計變量的試驗設計是實現試驗點在求解空間的合理分布的數理統計技術,包括三個子步驟,即選擇試驗設計方法S31,設置試驗設計參數S32和制定試驗設計表格S33。
            [0136]所述S31選擇試驗設計方法是從數理統計學科中各種試驗設計方法中選擇其中之一,以明確試驗設計點在設計空間分布規律。可選擇正交試驗設計,正交試驗設計是一種快速高效的試驗設計方法,是數理統計學科中的一門成熟技術。
            [0137]所述S32設置試驗設計參數是指設置試驗設計的因子數和設計變量的水平數,以確定試驗設計點的數量和分布。其中,因子數是指設計變量的個數,設計變量的水平數是指試驗設計中該設計變量在取值范圍選定試驗值的個數。選定實驗設計方法后,試驗設計點的個數取決于因子數和設計變量的水平數。根據所選擇的正交試驗設計,對于5個設計變量,因子數為5,選定每一個設計變量的水平數為3,可選取水平數為3的常用正交表,以確定試驗樣本點的個數。試驗樣本點越多,分布越均勻,越能反映出設計變量與優化目標函數、設計變量與設計約束之間的關系,但相應的計算效率也會越低。
            [0138]所述S33制定試驗設計表格是指根據選擇的試驗設計方法和試驗設計參數設置,指定出試驗設計表格。對于5因子3水平的正交試驗設計,其試驗表格可記為L27(3)因子數為n,每個設計變量的水平數為m的正交試驗表格中包含的每一個試驗設計點中各個設計變量的具體取值。
            [0139]在完成以上三個子步驟后,將試驗設計表格的設計點進行空間飛網裝置展開過程的仿真分析S2。空間飛網裝置展開過程的仿真分析S2結束后,獲得設計點的仿真結果。
            [0140]步驟四:構建近似數學函數模型S4
            [0141]S4構建近似數學函數模型是指根據試驗設計得到的設計變量與優化目標函數、設計變量與設計約束之間的關系,利用數學函數擬合,獲得可替代空間飛網展開過程分析的近似數學函數模型。構建近似數學函數模型S4包括2個子步驟:選擇構建近似數學函數模型的方法S41 ;擬合數據S42。
            [0142]所述S41選擇構建近似數學函數模型的類型是指選擇數據擬合的方法,實現設計變量與優化目標函數、設計變量與設計約束之間的關系的最佳,例如可選擇二階響應面函數近似模型。
            [0143]所述S42擬合數據是指利用S41選擇的方法,對試驗設計所有設計方案與結果進行擬合,構建出可替代空間飛網展開過程分析的近似數學函數模型。[0144]步驟五:S5近似模型計算分析
            [0145]S5近似模型計算分析是指在獲得設計方案后,利用近似數學函數模型給出對應的優化目標函數值和設計約束值。
            [0146]步驟六:S6優化
            [0147]S6優化主要是實現優化循環與收斂判斷的過程,包括S61選擇優化路徑、S62選擇優化算法、S63更新設計方案和S64優化收斂判斷這四個子步驟。
            [0148]S61選擇優化路徑是指優化過程中從仿真分析和近似數學函數模型兩種中選擇一種路徑獲得設計方法的結果反饋。
            [0149]S62選擇優化算法是指針對選擇的優化路徑,從優化算法中選出一種或多種方法用于優化迭代計算中。
            [0150]S63更新設計方案是指根據選擇的優化算法,計算得出新的設計方案的過程。
            [0151]S64優化收斂判斷是指計算前后兩組設計方案所得到的優化目標函數值的差的絕對值,并判斷如果結果小于或等于指定收斂精度,則停止優化,最后設計方案為最優設計方案,否則更新設計方案后進行S2空間飛網裝置展開過程的仿真分析或者S5近似模型計算分析,獲得相應的結果反饋。其中,收斂精度一般選用0.1%?10%,收斂精度越高,優化效果越好,但相應地會延長系統工作時間。
            [0152]參見圖4,本發明空間飛網裝置參數優化方法的流程圖:
            [0153]1、根據優化對象空間飛網裝置定義優化問題;
            [0154]2、根據初始設計方案進行仿真分析;
            [0155]3、選擇是否進行試驗設計,如果進行試驗設計,進入下一步,否則進入第17步;
            [0156]4、制定空間飛網裝置系統參數試驗設計方案表;
            [0157]5、對試驗設計方案表中的方案進行仿真分析;
            [0158]6、判斷是否所有的方案仿真都完成,若已完成則進入下一步,否則繼續回到第5步對下一個試驗設計方案進行仿真分析;
            [0159]7、選擇是否構建近似模型,如果是則進入下一步,否則進入第18步;
            [0160]8、構建近似模型;
            [0161]9、將試驗點中最佳結果點作為優化的初始點;
            [0162]10、選擇優化算法;
            [0163]11、根據所選擇的優化算法,更新設計點;
            [0164]12、將新的設計點的代入近似模型進行計算分析;
            [0165]13、判斷優化計算是否收斂,若優化不收斂則返回第11步,若優化已經收斂表明已經求解最優解,則進入下一步;
            [0166]14、將優化得到的最優點進行仿真分析;
            [0167]15、判斷最優點是否可信,若可信則結束,若不可信,則返回第5步;
            [0168]16、將試驗點中最佳結果點作為優化的初始點;
            [0169]17、選擇優化算法;
            [0170]18、根據所選擇的優化算法,更新設計點;
            [0171]19、對設計點進行仿真分析;
            [0172]20、判斷優化計算是否收斂,若優化不收斂則返回第18步,若優化已經收斂表明已經求解最優解,則結束優化。
            [0173]綜上所述,本發明一種空間飛網裝置參數優化方法包括以下六個步驟:S1空間飛網裝置參數優化模型定義;S2空間飛網裝置展開過程的仿真分析;S3空間飛網裝置設計變量的試驗設計;S4構建近似數學函數模型;S5近似模型計算分析;S6優化。
            【權利要求】
            1.一種空間飛網裝置參數優化設計系統,其特征在于:它包括以下4個部分=OPSNCSV1.0人機交互模塊Ml,OPSNCS V1.0控制與計算模塊M2,OPSNCS V1.0數據庫模塊M3和OPSNCS V1.0后處理模塊M4 ;OPSNCS V1.0人機交互模塊Ml將用戶輸入的優化模型的參數、運行環境參數和軟件運行控制指令信息記錄并傳遞給OPSNCS V1.0控制與計算模塊M2,OPSNCS V1.0控制與計算模塊M2根據這些信息從OPSNCS V1.0數據庫模塊M3選擇優化算法、試驗設計方法、近似模型函數,再通過OPSNCS V1.0控制與計算模塊M2的優化子模塊更新設計點,OPSNCS V1.0控制與計算模塊M2的仿真分析子模塊計算新設計點的仿真結果,通過這種循環計算獲得最優設計點;0PSNCS V1.0后處理模塊M4讀取最優設計點、最優設計和尋優過程的信息,并進行分析;用戶可通過OPSNCS V1.0人機交互模塊Ml查看OPSNCS V1.0后處理模塊M4分析與匯總的結果; 所述OPSNCS V1.0人機交互模塊Ml是用戶對軟件進行設置、控制與查看的模塊,包括以下三個子模塊:優化定義模塊M11,運行環境參數設置模塊M12,軟件運行控制工具M13,它們之間是相互并列的關系; 該優化定義模塊Mll是對空間飛網裝置參數優化中設計變量、優化目標和設計約束的定義;在優化定義模塊Mll中,定義了空間飛網裝置中牽引塊、艙蓋和空間網的發射參數共6個設計變量,數學表達如下:X = {mTm, , mh, vTm, Θ Tm, vh}; 其中,X為設計變量組,由五個變量組成:mTm為空間飛網裝置中牽引塊的質量,浮點變量,單位為kg ;mh為空間飛網裝置中網艙蓋的質量,浮點變量,單位為kg ;vTm為空間飛網裝置中牽引塊發射速度的大小,浮點變量,單位為m/s ; θΤπ為空間飛網裝置中牽引塊發射速度的角度,浮點變量,單位為度;vh為空間飛網裝置中牽引塊發射速度的大小,浮點變量,單位為m/s ;各變量的取值范圍需要根據空間飛網裝置的結構、材料、工況以及捕獲任務設定;在優化定義模塊Mll中,定義了空間飛網裝置參數優化目標,也即通過優化設計變量,使得空間網在指定距離 處展開到最大面積,數學表達如下:Max㈧ 其中,A為A為空間網在指定位置處展開的面積,單位為m2 ;在優化定義模塊MlI中,定義了設計約束,分別從空間網的強度和展開過程的效果出發,強度滿足空間網的材料性能,以空間網的最大內應力小于材料的許用應力作為約束,展開效果要求空間網在展開過程中不出現空間網面積收縮和空間網纏繞的情況,以空間網在到達指定位置前展開面積變化率大于零作為約束;其數學表達式如下:C = (C1,…,cm, ac^…,ack); 其中,C為約束條件,由2部分組成:Cl,-,Cffl為強度約束,分別對應著空間網的各段網繩的最大應力小于網繩材料的許用應力;aCl,…,ack為空間網展開過程的效果約束,對應的是在空間網在到達指定位置前的各個時刻空間網展開面積變化率應大于零;各約束具體取值與空間網的材料屬性、發射參數有關; 該運行環境參數設置模塊M12是對空間飛網裝置參數優化系統OPSNCS V1.0運行環境進行設置的模塊,包括求解器的選擇、文件存儲路徑的設置; 該軟件運行控制工具M13是對空間飛網裝置參數優化系統OPSNCS V1.0運動過程的檢查、開始、中止、終止過程的控制模塊; 所述OPSNCS V1.0控制與計算模塊M2是空間飛網裝置參數優化系統的核心模塊,實現對優化過程中數據的管理與控制;0PSNCS V1.0控制與計算模塊M2包括以下三個子模塊:數據管理模塊M21,優化模塊M22以及仿真分析模塊M23,它們之間是相互并列的關系;該數據管理模塊M21實現空間飛網裝置參數優化過程數據流的控制,以及訪問空間飛網裝置參數優化系統的數據庫模塊M3 ;該優化模塊M22是利用數據管理模塊M21選擇的優化算法計算出下一個設計方案的程序模塊;該仿真分析模塊M23是對設計方案進行空間網裝置展開過程的仿真與分析; 所述OPSNCS V1.0數據庫模塊M3包括三個子模塊:優化算法庫M31,試驗設計方法庫M32,近似模型函數庫M33 ;它們之間是相互并列的關系;該優化算法庫M31是包括常用的全局算法和常用的梯度算法;該試驗設計方法庫M32是包括常用的試驗設計方法的方法數據庫;該近似模型函數庫M33是包括數據擬合常用函數的函數數據庫; 所述OPSNCS V1.0后處理模塊M4包括兩個子模塊:結果讀取模塊M41和結果分析模塊M42 ;它們之間是相互并列的關系; 該結果讀取模塊M41是在優化循環停止后,查詢優化和仿真結果文件并提取最優設計結果的程序模塊;該結果分析模塊M42是根據結果讀取模塊M41讀取到的信息,分析優化結果的可信度的程序模塊。
            2.一種空間飛網裝置參數優化方法,其特征在于:該方法具體步驟如下: 步驟一:空間飛網裝置參數優化模型定義SI ; 其中,優化模型SI的定義包括優化目標函數S11,優化約束條件S12,優化設計變量S13,二者之間是并列關系; 優化目標函數Sll是指空間飛網裝置參數優化的目標,也就是空間網在指定距離處展開到最大面積,其數學表達式如下=Max(A); 其中,A為空間網在指定位置處展開的面積,單位為m2 ; 優化約束條件S12是指空間飛網裝置參數優化的約束條件,分別從空間網的強度和展開過程的效果出發,強度滿足空間網的材料性能,以空間網的最大內應力小于材料的許用應力作為約束,展開效果要求空間網在展開過程中不出現空間網面積收縮和空間網纏繞的情況,以空間網在到達指定位置前展開面積變化率大于零作為約束;其數學表達式如下:
            C (。1,...,Cm, SIC”...,BCj^); 其中,C為約束條件,由2部分組成:Cl,-,Cffl為強度約束,分別對應著空間網的各段網繩的最大應力小于網繩材料的許用應力;aCl,…,ack為空間網展開過程的效果約束,對應的是在空間網在到達指定位置前的各個時刻空間網展開面積變化率應大于零;各約束具體取值與空間網的材料屬性、發射參數有關; 優化設計變量S13是指空間飛網裝置參數優化的各個設計變量,包括空間飛網裝置的要進行優化設計的發射速度參數和結構參數,其數學表達式如下:
            X — {Xj, X2,X3...XrJ ; 其中,X為設計變量組,由復數個變量組成:Xi(i = 1,...,η)的具體的物理含義根據不同的空間飛網裝置設計任務而不同,Xl為空間飛網裝置中牽引塊的質量,單位為kg ;χ2為空間飛網裝置中網艙蓋質量,單位為kg ;x3為空間飛網裝置牽引塊的發射速度的大小,單位為m/s ;x4為空間飛網裝置牽引塊的發射角度,單位為度;x5為空間飛網裝置網艙蓋發射速度的大小,單位為m/s ;x6為空間飛網裝置中牽引塊和網艙蓋發射時間差,單位為S,具體使用時設計變量不限以上幾種; 步驟二:空間飛網裝置展開過程的仿真分析S2 ;空間飛網裝置展開過程的仿真分析S2是對空間網發射展開過程進行仿真分析和提取優化所需要分析結果,包括六個子步驟:定義任務參數S21,定義結構尺寸S22,定義空間網的材料屬性S23,定義仿真單元屬性S24,定義載荷工況S25,分析求解S26和結果讀取S27 ;定義任務參數S21是指空間飛網任務的參數,包括捕獲的目標的最大展長和距離,捕獲目標的最大展長單位為m,取值范圍5m~30m ;捕獲目標的距離單位為m,取值范圍30m~120m ; 定義結構尺寸S22是指空間網的形狀、大小、網格密度以及網繩的粗細,正方形網的邊長,單位為m,取值為20m~40m ;網格密度是指空間網網格疏密程度,通過網邊單位為個,取值為整數,取值范圍為20~100 ; 定義空間網的材料屬性S23是實現空間飛網裝置材料屬性的定義,包括空間網網繩的材料定義、牽引塊的材料定義、網艙蓋的材料定義,其中材料定義是指指定材料的彈性模量、密度和泊松比參數; 定義仿真單元屬性S24實現對仿真結構模型單元的定義,包括結構仿真模型各部分的單元定義;牽引塊和網艙蓋使用集中質點單元,網繩使用繩索單元; 定義分析求解S25實現空間飛網裝置 在定義的發射參數下空間網展開過程求解,選用求解器對空間網的展開過程求解,其中,求解器選用現有成熟軟件; 定義結果讀取S26實現對分析結果進行處理,在空間飛網裝置分析求解完畢后,求解器會輸出從空間網發射開始到空間網到達指定距離所經歷的時間段內,空間網展開面積和繩索單元的應力隨時間變化的信息; 步驟三:空間飛網裝置設計變量的試驗設計S3 空間飛網裝置設計變量的試驗設計S3是實現試驗點在求解空間的合理分布的數理統計技術,包括三個子步驟,即選擇試驗設計方法S31,設置試驗設計參數S32和制定試驗設計表格S33 ; 定義選擇試驗設計方法S31是從數理統計學科中各種試驗設計方法中選擇其中之一,以明確試驗設計點在設計空間分布規律; 定義設置試驗設計參數S32是指設置試驗設計的因子數和設計變量的水平數,以確定試驗設計點的數量和分布;其中,因子數是指設計變量的個數,設計變量的水平數是指試驗設計中該設計變量在取值范圍選定試驗值的個數;選定實驗設計方法后,試驗設計點的個數取決于因子數和設計變量的水平數; 定義制定試驗設計表格S33是指根據選擇的試驗設計方法和試驗設計參數設置,指定出試驗設計表格; 在完成以上三個子步驟后,將試驗設計表格的設計點進行空間飛網裝置展開過程的仿真分析S2,空間飛網裝置展開過程的仿真分析S2結束后,獲得設計點的仿真結果; 步驟四:構建近似數學函數模型S4 構建近似數學函數模型S4是指根據試驗設計得到的設計變量與優化目標函數、設計變量與設計約束之間的關系,利用數學函數擬合,獲得可替代空間飛網展開過程分析的近似數學函數模型;構建近似數學函數模型S4包括2個子步驟:選擇構建近似數學函數模型的方法S41 ;擬合數據S42 ; 定義選擇構建近似數學函數模型的類型S41是指選擇數據擬合的方法,實現設計變量與優化目標函數、設計變量與設計約束之間的關系的最佳; 定義擬合數據S42是指利用S41選擇的方法,對試驗設計所有設計方案與結果進行擬合,構建出可替代空間飛網展開過程仿真分析的近似數學函數模型; 步驟五:近似模型計算分析S5 近似模型計算分析S5是指在獲得設計方案后,利用近似數學函數模型給出對應的優化目標函數值和設計約束值; 步驟六:優化S6 優化S6是實現優化循環與收斂判斷的過程,包括選擇優化路徑S61、選擇優化算法S62、更新設計方案S63和優化收斂判斷S64這四個子步驟; 選擇優化路徑S61是指優化過程中從仿真分析和近似數學函數模型兩種中選擇一種路徑獲得設計方法的結果反饋; 選擇優化算法S62是指針對選擇的優化路徑,從優化算法中選出一種或多種方法用于優化迭代計算中; 更新設計方案S63是指根據選擇的優化算法,計算得出新的設計方案的過程; 優化收斂判斷S64是指計算前后兩組設計方案所得到的優化目標函數值的差的絕對值,并判斷如果結果小于或等于指定收斂精度,則停止優化,最后設計方案為最優設計方案,否則更新設計方案后 進行空間飛網裝置展開過程的仿真分析S2或者近似模型計算分析S5,獲得相應的結果反饋;其中,收斂精度選用0.1%~10%,收斂精度越高,優化效果越好,但相應地會延長系統工作時間。
            【文檔編號】G06F17/50GK104008233SQ201410197237
            【公開日】2014年8月27日 申請日期:2014年5月12日 優先權日:2014年5月12日
            【發明者】王曉慧, 萬長煌 申請人:北京航空航天大學
            網友詢問留言 已有0條留言
            • 還沒有人留言評論。精彩留言會獲得點贊!
            1
            婷婷六月激情在线综合激情,亚洲国产大片,久久中文字幕综合婷婷,精品久久久久久中文字幕,亚洲一区二区三区高清不卡,99国产精品热久久久久久夜夜嗨 ,欧美日韩亚洲综合在线一区二区,99国产精品电影,伊人精品线视天天综合,精品伊人久久久大香线蕉欧美
            亚洲精品1区 国产成人一级 91精品国产欧美一区二区 亚洲精品乱码久久久久久下载 国产精品久久久久久久伊一 九色国产 国产精品九九视频 伊人久久成人爱综合网 欧美日韩亚洲区久久综合 欧美日本一道免费一区三区 夜夜爽一区二区三区精品 欧美日韩高清一区二区三区 国产成人av在线 国产精品对白交换绿帽视频 国产视频亚洲 国产在线欧美精品 国产精品综合网 国产日韩精品欧美一区色 国产日韩精品欧美一区喷 欧美日韩在线观看区一二 国产区精品 欧美视频日韩视频 中文字幕天天躁日日躁狠狠躁97 视频一二三区 欧美高清在线精品一区二区不卡 国产精品揄拍一区二区久久 99久久综合狠狠综合久久aⅴ 亚洲乱码视频在线观看 日韩在线第二页 亚洲精品无码专区在线播放 成人亚洲网站www在线观看 欧美三级一区二区 99久久精品免费看国产高清 91麻豆国产在线观看 最新日韩欧美不卡一二三区 成人在线观看不卡 日韩国产在线 在线亚洲精品 亚洲午夜久久久久中文字幕 国产精品成人久久久久久久 精品国产一区二区在线观看 欧美精品国产一区二区三区 中文在线播放 亚洲第一页在线视频 国产午夜精品福利久久 九色国产 精品国产九九 国产永久视频 久久精品人人做人人综合试看 国产一区二区三区免费观看 亚洲精品国产电影 9999热视频 国产精品资源在线 麻豆久久婷婷国产综合五月 国产精品免费一级在线观看 亚洲国产一区二区三区青草影视 中文在线播放 国产成人综合在线 国产在线观看色 国产亚洲三级 国产片一区二区三区 久久99精品久久久久久牛牛影视 亚洲欧美日韩国产 四虎永久免费网站 国产一毛片 国产精品视频在 九九热在线精品 99精品福利视频 色婷婷色99国产综合精品 97成人精品视频在线播放 精品久久久久久中文字幕 亚洲欧美一区二区三区孕妇 亚洲欧美成人网 日韩高清在线二区 国产尤物在线观看 在线不卡一区二区 91网站在线看 韩国精品福利一区二区 欧美日韩国产成人精品 99热精品久久 国产精品免费视频一区 高清视频一区 精品九九久久 欧美日韩在线观看免费 91欧美激情一区二区三区成人 99福利视频 亚洲国产精品91 久热国产在线 精品久久久久久中文字幕女 国产精品久久久久久久久99热 成人自拍视频网 国产精品视频久久久久久 久久影院国产 国产玖玖在线观看 99精品在线免费 亚洲欧美一区二区三区导航 久久久久久久综合 国产欧美日韩精品高清二区综合区 国产精品视频自拍 亚洲一级片免费 久久久久久九九 国产欧美自拍视频 视频一区二区在线观看 欧美日韩一区二区三区久久 中文在线亚洲 伊人热人久久中文字幕 日韩欧美亚洲国产一区二区三区 欧美亚洲国产成人高清在线 欧美日韩国产码高清综合人成 国产性大片免费播放网站 亚洲午夜综合网 91精品久久一区二区三区 国产无套在线播放 国产精品视频网站 国产成人亚洲精品老王 91在线网站 国产视频97 欧美黑人欧美精品刺激 国产一区二区三区免费在线视频 久久久国产精品免费看 99re6久精品国产首页 久久精品91 国产成人一级 国产成人精品曰本亚洲 日本福利在线观看 伊人成综合网 久久综合一本 国产综合久久久久久 久久精品成人免费看 久久福利 91精品国产91久久久久久麻豆 亚洲精品成人在线 亚洲伊人久久精品 欧美日本二区 国产永久视频 国产一区二 一区二区福利 国产一毛片 亚洲精品1区 毛片一区二区三区 伊人久久大香线蕉综合影 国产欧美在线观看一区 亚洲国产欧洲综合997久久 国产一区二区免费视频 国产91精品对白露脸全集观看 久久亚洲国产伦理 欧美成人伊人久久综合网 亚洲性久久久影院 久久99国产精一区二区三区! 91精品国产欧美一区二区 欧美日韩亚洲区久久综合 日韩精品一二三区 久久久夜色精品国产噜噜 国产在线精品福利91香蕉 久久久久久久亚洲精品 97se色综合一区二区二区 91国语精品自产拍在线观看性色 91久久国产综合精品女同我 日韩中文字幕a 国产成人亚洲日本精品 久久国产精品-国产精品 久久国产经典视频 久久国产精品伦理 亚洲第一页在线视频 国产精品久久久久三级 日韩毛片网 久久免费高清视频 麻豆国产在线观看一区二区 91麻豆国产福利在线观看 国产成人精品男人的天堂538 一区二区三区中文字幕 免费在线视频一区 欧美日韩国产成人精品 国产综合网站 国产资源免费观看 亚洲精品亚洲人成在线播放 精品久久久久久中文字幕专区 亚洲人成人毛片无遮挡 国产一起色一起爱 国产香蕉精品视频在 九九热免费观看 日韩亚洲欧美一区 九九热精品在线观看 精品久久久久久中文字幕专区 亚洲欧美自拍偷拍 国产精品每日更新 久久久久国产一级毛片高清板 久久天天躁狠狠躁夜夜中文字幕 久久精品片 日韩在线毛片 国产成人精品本亚洲 国产成人精品一区二区三区 九九热在线观看 国产r级在线观看 国产欧美日韩精品高清二区综合区 韩国电影一区二区 国产精品毛片va一区二区三区 五月婷婷伊人网 久久一区二区三区免费 一本色道久久综合狠狠躁篇 亚洲综合色站 国产尤物在线观看 亚洲一区亚洲二区 免费在线视频一区 欧洲精品视频在线观看 日韩中文字幕a 中文字幕日本在线mv视频精品 91精品在线免费视频 精品国产免费人成在线观看 精品a级片 中文字幕日本在线mv视频精品 日韩在线精品视频 婷婷丁香色 91精品国产高清久久久久 国产成人精品日本亚洲直接 五月综合视频 欧美日韩在线亚洲国产人 精液呈暗黄色 亚洲乱码一区 久久精品中文字幕不卡一二区 亚洲天堂精品在线 激情婷婷综合 国产免费久久精品久久久 国产精品亚洲二区在线 久久免费播放视频 五月婷婷丁香综合 在线亚洲欧美日韩 久久免费精品高清麻豆 精品久久久久久中文字幕 亚洲一区网站 国产精品福利社 日韩中文字幕免费 亚洲综合丝袜 91精品在线播放 国产精品18 亚洲日日夜夜 伊人久久大香线蕉综合影 亚洲精品中文字幕乱码影院 亚洲一区二区黄色 亚洲第一页在线视频 一区二区在线观看视频 国产成人福利精品视频 亚洲高清二区 国内成人免费视频 精品亚洲性xxx久久久 国产精品合集一区二区三区 97av免费视频 国产一起色一起爱 国产区久久 国产资源免费观看 99精品视频免费 国产成人一级 国产精品九九免费视频 欧美91精品久久久久网免费 99热国产免费 久久精品色 98精品国产综合久久 久久精品播放 中文字幕视频免费 国产欧美日韩一区二区三区在线 精品久久蜜桃 国产小视频精品 一本色道久久综合狠狠躁篇 91在线免费观看 亚洲精品区 伊人成综合网 伊人热人久久中文字幕 伊人黄色片 99国产精品热久久久久久夜夜嗨 久久免费精品视频 亚洲一区二区三区高清不卡 久久久久国产一级毛片高清板 国产片一区二区三区 久久狠狠干 99久久婷婷国产综合精品电影 国产99区 国产精品成人久久久久 久久狠狠干 青青国产在线观看 亚洲高清国产拍精品影院 国产精品一区二区av 九九热在线免费视频 伊人久久国产 国产精品久久久久久久久久一区 在线观看免费视频一区 国产精品自在在线午夜区app 国产精品综合色区在线观看 国产毛片久久久久久国产毛片 97国产免费全部免费观看 国产精品每日更新 国产尤物视频在线 九九视频这里只有精品99 一本一道久久a久久精品综合 久久综合给会久久狠狠狠 国产成人精品男人的天堂538 欧美一区二区高清 毛片一区二区三区 国产欧美日韩在线观看一区二区三区 在线国产二区 欧美不卡网 91在线精品中文字幕 在线国产福利 国内精品91久久久久 91亚洲福利 日韩欧美国产中文字幕 91久久精品国产性色也91久久 亚洲性久久久影院 欧美精品1区 国产热re99久久6国产精品 九九热免费观看 国产精品欧美日韩 久久久久国产一级毛片高清板 久久国产经典视频 日韩欧美亚洲国产一区二区三区 欧美亚洲综合另类在线观看 国产精品自在在线午夜区app 97中文字幕在线观看 视频一二三区 精品国产一区在线观看 国产欧美日韩在线一区二区不卡 欧美一区二三区 伊人成人在线观看 国内精品91久久久久 97在线亚洲 国产在线不卡一区 久久久全免费全集一级全黄片 国产精品v欧美精品∨日韩 亚洲毛片网站 在线不卡一区二区 99re热在线视频 久久激情网 国产毛片一区二区三区精品 久久亚洲综合色 中文字幕视频免费 国产视频亚洲 婷婷伊人久久 国产一区二区免费播放 久久99国产精品成人欧美 99国产在线视频 国产成人免费视频精品一区二区 国产不卡一区二区三区免费视 国产码欧美日韩高清综合一区 久久精品国产主播一区二区 国产一区电影 久久精品国产夜色 国产精品国产三级国产 日韩一区二区三区在线 久久97久久97精品免视看 久久国产免费一区二区三区 伊人久久大香线蕉综合电影网 99re6久精品国产首页 久久激情网 亚洲成人高清在线 国产精品网址 国产成人精品男人的天堂538 香蕉国产综合久久猫咪 国产专区中文字幕 91麻豆精品国产高清在线 久久国产经典视频 国产精品成人va在线观看 国产精品爱啪在线线免费观看 日本精品久久久久久久久免费 亚洲综合一区二区三区 久久五月网 精品国产网红福利在线观看 久久综合亚洲伊人色 亚洲国产精品久久久久久网站 在线日韩国产 99国产精品热久久久久久夜夜嗨 国产综合精品在线 国产区福利 精品亚洲综合久久中文字幕 国产制服丝袜在线 毛片在线播放网站 在线观看免费视频一区 国产精品久久久精品三级 亚洲国产电影在线观看 最新日韩欧美不卡一二三区 狠狠综合久久综合鬼色 日本精品1在线区 国产日韩一区二区三区在线播放 欧美日韩精品在线播放 亚洲欧美日韩国产一区二区三区精品 久久综合久久网 婷婷六月激情在线综合激情 亚洲乱码一区 国产专区91 97av视频在线观看 精品久久久久久中文字幕 久久五月视频 国产成人福利精品视频 国产精品网址 中文字幕视频在线 精品一区二区三区免费视频 伊人手机在线视频 亚洲精品中文字幕乱码 国产在线视频www色 色噜噜国产精品视频一区二区 精品亚洲成a人在线观看 国产香蕉尹人综合在线 成人免费一区二区三区在线观看 国产不卡一区二区三区免费视 欧美精品久久天天躁 国产专区中文字幕 久久精品国产免费中文 久久精品国产免费一区 久久无码精品一区二区三区 国产欧美另类久久久精品免费 欧美精品久久天天躁 亚洲精品在线视频 国产视频91在线 91精品福利一区二区三区野战 日韩中文字幕免费 国产精品99一区二区三区 欧美成人高清性色生活 国产精品系列在线观看 亚洲国产福利精品一区二区 国产成人在线小视频 国产精品久久久久免费 99re热在线视频 久久久久久久综合 一区二区国产在线播放 成人国产在线视频 亚洲精品乱码久久久久 欧美日韩一区二区综合 精品久久久久免费极品大片 中文字幕视频二区 激情粉嫩精品国产尤物 国产成人精品一区二区视频 久久精品中文字幕首页 亚洲高清在线 国产精品亚洲一区二区三区 伊人久久艹 中文在线亚洲 国产精品一区二区在线播放 国产精品九九免费视频 亚洲二区在线播放 亚洲狠狠婷婷综合久久久久网站 亚洲欧美日韩网站 日韩成人精品 亚洲国产一区二区三区青草影视 91精品国产福利在线观看 国产精品久久久久久久久99热 国产一区二区精品尤物 久碰香蕉精品视频在线观看 亚洲日日夜夜 在线不卡一区二区 国产午夜亚洲精品 九九热在线视频观看这里只有精品 伊人手机在线视频 91免费国产精品 日韩欧美中字 91精品国产91久久久久 国产全黄三级播放 视频一区二区三区免费观看 国产开裆丝袜高跟在线观看 国产成人欧美 激情综合丝袜美女一区二区 国产成人亚洲综合无 欧美精品一区二区三区免费观看 欧美亚洲国产日韩 日韩亚州 国产欧美日韩精品高清二区综合区 亚洲午夜国产片在线观看 精品久久久久久中文字幕 欧美精品1区 久久伊人久久亚洲综合 亚洲欧美日韩精品 国产成人精品久久亚洲高清不卡 久久福利影视 国产精品99精品久久免费 久久久久免费精品视频 国产日产亚洲精品 亚洲国产午夜电影在线入口 精品无码一区在线观看 午夜国产精品视频 亚洲一级片免费 伊人久久大香线蕉综合影 国产精品久久影院 久碰香蕉精品视频在线观看 www.欧美精品 在线小视频国产 亚洲国产天堂久久综合图区 欧美一区二区三区不卡 日韩美女福利视频 九九精品免视频国产成人 不卡国产00高中生在线视频 亚洲第一页在线视频 欧美日韩在线播放成人 99re视频这里只有精品 国产精品91在线 精品乱码一区二区三区在线 国产区久久 91麻豆精品国产自产在线观看一区 日韩精品成人在线 九九热在线观看 国产精品久久不卡日韩美女 欧美一区二区三区综合色视频 欧美精品免费一区欧美久久优播 国产精品网址 国产专区中文字幕 国产精品欧美亚洲韩国日本久久 日韩美香港a一级毛片 久久精品123 欧美一区二区三区免费看 99r在线视频 亚洲精品国产字幕久久vr 国产综合激情在线亚洲第一页 91免费国产精品 日韩免费小视频 亚洲国产精品综合一区在线 国产亚洲第一伦理第一区 在线亚洲精品 国产精品一区二区制服丝袜 国产在线成人精品 九九精品免视频国产成人 亚洲国产网 欧美日韩亚洲一区二区三区在线观看 在线亚洲精品 欧美一区二区三区高清视频 国产成人精品男人的天堂538 欧美日韩在线观看区一二 亚洲欧美一区二区久久 久久精品中文字幕首页 日本高清www午夜视频 久久精品国产免费 久久999精品 亚洲国产精品欧美综合 88国产精品视频一区二区三区 91久久偷偷做嫩草影院免费看 国产精品夜色视频一区二区 欧美日韩导航 国产成人啪精品午夜在线播放 一区二区视频在线免费观看 99久久精品国产自免费 精液呈暗黄色 久久99国产精品 日本精品久久久久久久久免费 精品国产97在线观看 99re视频这里只有精品 国产视频91在线 999av视频 亚洲美女视频一区二区三区 久久97久久97精品免视看 亚洲国产成人久久三区 99久久亚洲国产高清观看 日韩毛片在线视频 综合激情在线 91福利一区二区在线观看 一区二区视频在线免费观看 激情粉嫩精品国产尤物 国产成人精品曰本亚洲78 国产成人精品本亚洲 国产精品成人免费视频 国产成人啪精品视频免费软件 久久精品国产亚洲妲己影院 国产精品成人久久久久久久 久久大香线蕉综合爱 欧美一区二区三区高清视频 99热国产免费 在线观看欧美国产 91精品视频在线播放 国产精品福利社 欧美精品一区二区三区免费观看 国产一区二区免费视频 国产午夜精品一区二区 精品视频在线观看97 91精品福利久久久 国产一区福利 国产综合激情在线亚洲第一页 国产精品久久久久久久久久久不卡 九色国产 在线日韩国产 黄网在线观看 亚洲一区小说区中文字幕 中文字幕丝袜 日本二区在线观看 日本国产一区在线观看 欧美日韩一区二区三区久久 欧美精品亚洲精品日韩专 国产日产亚洲精品 久久综合九色综合欧美播 亚洲国产欧美无圣光一区 欧美视频区 亚洲乱码视频在线观看 久久无码精品一区二区三区 九九热精品免费视频 久久99精品久久久久久牛牛影视 国产精品成久久久久三级 国产一区福利 午夜国产精品视频 日本二区在线观看 99久久网站 国产亚洲天堂 精品国产一区二区三区不卡 亚洲国产日韩在线一区 国产成人综合在线观看网站 久久免费高清视频 欧美在线导航 午夜精品久久久久久99热7777 欧美久久综合网 国产小视频精品 国产尤物在线观看 亚洲国产精品综合一区在线 欧美一区二区三区不卡视频 欧美黑人欧美精品刺激 日本福利在线观看 久久国产偷 国产手机精品一区二区 国产热re99久久6国产精品 国产高清啪啪 欧美亚洲国产成人高清在线 国产在线第三页 亚洲综合一区二区三区 99r在线视频 99精品久久久久久久婷婷 国产精品乱码免费一区二区 国产在线精品福利91香蕉 国产尤物视频在线 五月婷婷亚洲 中文字幕久久综合伊人 亚洲精品一级毛片 99国产精品电影 在线视频第一页 久久99国产精品成人欧美 国产白白视频在线观看2 成人精品一区二区www 亚洲成人网在线观看 麻豆91在线视频 色综合合久久天天综合绕视看 久久精品国产免费高清 国产不卡一区二区三区免费视 欧美国产中文 99精品欧美 九九在线精品 国产中文字幕在线免费观看 国产一区中文字幕在线观看 国产成人一级 国产精品一区二区制服丝袜 国产一起色一起爱 亚洲精品成人在线 亚洲欧美精品在线 国产欧美自拍视频 99精品久久久久久久婷婷 久99视频 国产热re99久久6国产精品 视频一区亚洲 国产精品视频分类 国产精品成在线观看 99re6久精品国产首页 亚洲在成人网在线看 亚洲国产日韩在线一区 久久国产三级 日韩国产欧美 欧美在线一区二区三区 国产精品美女一级在线观看 成人午夜免费福利视频 亚洲天堂精品在线 91精品国产手机 欧美日韩视频在线播放 狠狠综合久久综合鬼色 九一色视频 青青视频国产 亚洲欧美自拍一区 中文字幕天天躁日日躁狠狠躁97 日韩免费大片 996热视频 伊人成综合网 亚洲天堂欧美 日韩精品亚洲人成在线观看 久久综合给会久久狠狠狠 日韩精品亚洲人成在线观看 日韩国产欧美 亚洲成aⅴ人片在线影院八 亚洲精品1区 99久久精品免费 国产精品高清在线观看 国产精品久久久免费视频 在线亚洲欧美日韩 91在线看视频 国产精品96久久久久久久 欧美日韩国产成人精品 91在线亚洲 热久久亚洲 国产精品美女免费视频观看 日韩在线毛片 亚洲永久免费视频 九九免费在线视频 亚洲一区网站 日本高清二区视频久二区 精品国产美女福利在线 伊人久久艹 国产精品久久久久三级 欧美成人精品第一区二区三区 99久久精品国产自免费 在线观看日韩一区 国产中文字幕一区 成人免费午夜视频 欧美日韩另类在线 久久99国产精品成人欧美 色婷婷中文网 久久天天躁夜夜躁狠狠躁2020 欧美成人伊人久久综合网 国产精品福利资源在线 国产伦精品一区二区三区高清 国产精品亚洲综合色区韩国 亚洲一区欧美日韩 色综合视频 国语自产精品视频在线区 国产高清a 成人国内精品久久久久影 国产在线精品香蕉综合网一区 国产不卡在线看 国产成人精品精品欧美 国产欧美日韩综合精品一区二区三区 韩国电影一区二区 国产在线视频www色 91中文字幕在线一区 国产人成午夜免视频网站 亚洲综合一区二区三区 色综合视频一区二区观看 久久五月网 九九热精品在线观看 国产一区二区三区国产精品 99久热re在线精品996热视频 亚洲国产网 在线视频亚洲一区 日韩字幕一中文在线综合 国产高清一级毛片在线不卡 精品国产色在线 国产高清视频一区二区 精品日本久久久久久久久久 亚洲国产午夜精品乱码 成人免费国产gav视频在线 日韩欧美一区二区在线观看 欧美曰批人成在线观看 韩国电影一区二区 99re这里只有精品6 日韩精品一区二区三区视频 99re6久精品国产首页 亚洲欧美一区二区三区导航 欧美色图一区二区三区 午夜精品视频在线观看 欧美激情在线观看一区二区三区 亚洲热在线 成人国产精品一区二区网站 亚洲一级毛片在线播放 亚洲一区小说区中文字幕 亚洲午夜久久久久影院 国产自产v一区二区三区c 国产精品视频免费 久久调教视频 国产成人91激情在线播放 国产精品欧美亚洲韩国日本久久 久久亚洲日本不卡一区二区 91中文字幕网 成人国产在线视频 国产视频91在线 欧美成人精品第一区二区三区 国产精品福利在线 久久综合九色综合精品 欧美一区二区三区精品 久久国产综合尤物免费观看 久久99青青久久99久久 日韩精品免费 久久国产精品999 91亚洲视频在线观看 国产精品igao视频 色综合区 在线亚洲欧国产精品专区 国产一区二区三区在线观看视频 亚洲精品成人在线 一区二区国产在线播放 中文在线亚洲 亚洲精品第一国产综合野 国产一区二区精品久久 一区二区三区四区精品视频 99热精品久久 中文字幕视频二区 国产成人精品男人的天堂538 99精品影视 美女福利视频一区二区 久久午夜夜伦伦鲁鲁片 综合久久久久久久综合网 国产精品国产欧美综合一区 国产99视频在线观看 国产亚洲女在线精品 婷婷影院在线综合免费视频 国产亚洲3p一区二区三区 91成人爽a毛片一区二区 亚洲一区二区高清 国产欧美亚洲精品第二区首页 欧美日韩导航 亚洲高清二区 欧美激情观看一区二区久久 日韩毛片在线播放 亚洲欧美日韩高清中文在线 亚洲日本在线播放 国产精品一区二区制服丝袜 精品国产一区二区三区不卡 国产不卡在线看 国产欧美网站 四虎永久在线观看视频精品 国产黄色片在线观看 夜夜综合 一本色道久久综合狠狠躁篇 欧美亚洲综合另类在线观看 国产91在线看 伊人久久国产 欧美一区二区在线观看免费网站 国产精品久久久久三级 久久福利 日韩中文字幕a 亚洲午夜久久久久影院 91在线高清视频 国产亚洲一区二区三区啪 久久人精品 国产精品亚洲午夜一区二区三区 综合久久久久久 久久伊人一区二区三区四区 国产综合久久久久久 日韩一区精品视频在线看 国产精品日韩欧美制服 日本精品1在线区 99re视频 无码av免费一区二区三区试看 国产视频1区 日韩欧美中文字幕一区 日本高清中文字幕一区二区三区a 亚洲国产欧美无圣光一区 国产在线视频一区二区三区 欧美国产第一页 在线亚洲欧美日韩 日韩中文字幕第一页 在线不卡一区二区 伊人久久青青 国产精品一区二区在线播放 www.五月婷婷 麻豆久久婷婷国产综合五月 亚洲精品区 久久国产欧美另类久久久 99在线视频免费 伊人久久中文字幕久久cm 久久精品成人免费看 久久这里只有精品首页 88国产精品视频一区二区三区 中文字幕日本在线mv视频精品 国产在线精品成人一区二区三区 伊人精品线视天天综合 亚洲一区二区黄色 国产尤物视频在线 亚洲精品99久久久久中文字幕 国产一区二区三区免费观看 伊人久久大香线蕉综合电影网 国产成人精品区在线观看 日本精品一区二区三区视频 日韩高清在线二区 久久免费播放视频 一区二区成人国产精品 国产精品免费精品自在线观看 亚洲精品视频二区 麻豆国产精品有码在线观看 精品日本一区二区 亚洲欧洲久久 久久中文字幕综合婷婷 中文字幕视频在线 国产成人精品综合在线观看 91精品国产91久久久久福利 精液呈暗黄色 香蕉国产综合久久猫咪 国产专区精品 亚洲精品无码不卡 国产永久视频 亚洲成a人片在线播放观看国产 一区二区国产在线播放 亚洲一区二区黄色 欧美日韩在线观看视频 亚洲精品另类 久久国产综合尤物免费观看 国产一区二区三区国产精品 高清视频一区 国产精品igao视频 国产精品资源在线 久久综合精品国产一区二区三区 www.五月婷婷 精品色综合 99热国产免费 麻豆福利影院 亚洲伊人久久大香线蕉苏妲己 久久电影院久久国产 久久精品伊人 在线日韩理论午夜中文电影 亚洲国产欧洲综合997久久 伊人国产精品 久草国产精品 欧美一区精品二区三区 亚洲成人高清在线 91免费国产精品 日韩精品福利在线 国产一线在线观看 国产不卡在线看 久久99青青久久99久久 亚洲精品亚洲人成在线播放 99久久免费看国产精品 国产日本在线观看 青草国产在线视频 麻豆久久婷婷国产综合五月 国产中文字幕一区 91久久精品国产性色也91久久 国产一区a 国产欧美日韩成人 国产亚洲女在线精品 一区二区美女 中文字幕在线2021一区 在线小视频国产 久久这里只有精品首页 国产在线第三页 欧美日韩中文字幕 在线亚洲+欧美+日本专区 精品国产一区二区三区不卡 久久这里精品 欧美在线va在线播放 精液呈暗黄色 91精品国产手机 91在线免费播放 欧美视频亚洲色图 欧美国产日韩精品 日韩高清不卡在线 精品视频免费观看 欧美日韩一区二区三区四区 国产欧美亚洲精品第二区首页 亚洲韩精品欧美一区二区三区 国产精品视频免费 在线精品小视频 久久午夜夜伦伦鲁鲁片 国产无套在线播放 久热这里只精品99re8久 欧美久久久久 久久香蕉国产线看观看精品蕉 国产成人精品男人的天堂538 亚洲人成网站色7799在线观看 日韩在线第二页 一本色道久久综合狠狠躁篇 国产一区二区三区不卡在线观看 亚洲乱码在线 在线观看欧美国产 久久福利青草精品资源站免费 国产玖玖在线观看 在线亚洲精品 亚洲成aⅴ人在线观看 精品91在线 欧美一区二三区 日韩中文字幕视频在线 日本成人一区二区 日韩免费专区 国内精品在线观看视频 久久国产综合尤物免费观看 国产精品系列在线观看 一本一道久久a久久精品综合 亚洲免费播放 久久精品国产免费 久久人精品 亚洲毛片网站 亚洲成a人一区二区三区 韩国福利一区二区三区高清视频 亚洲精品天堂在线 一区二区三区中文字幕 亚洲国产色婷婷精品综合在线观看 亚洲国产成人久久笫一页 999国产视频 国产精品香港三级在线电影 欧美日韩一区二区三区四区 日韩国产欧美 国产精品99一区二区三区 午夜国产精品理论片久久影院 亚洲精品中文字幕麻豆 亚洲国产高清视频 久久免费手机视频 日韩a在线观看 五月婷婷亚洲 亚洲精品中文字幕麻豆 中文字幕丝袜 www国产精品 亚洲天堂精品在线 亚洲乱码一区 国产日韩欧美三级 久久999精品 伊人热人久久中文字幕 久热国产在线视频 国产欧美日韩在线观看一区二区三区 国产一二三区在线 日韩国产欧美 91精品国产91久久久久 亚洲一区小说区中文字幕 精品一区二区免费视频 国产精品视频免费 国产精品亚洲综合色区韩国 亚洲国产精品成人午夜在线观看 欧美国产日韩精品 中文字幕精品一区二区精品