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

    一種提高飛行器結(jié)構(gòu)天地力學(xué)環(huán)境地面試驗(yàn)有效性的方法及其應(yīng)用

    2016-11-12 02:12:19閻桂榮董龍雷宋利強(qiáng)
    裝備環(huán)境工程 2016年5期
    關(guān)鍵詞:振型載荷有限元

    閻桂榮,董龍雷,宋利強(qiáng)

    (西安交通大學(xué) 機(jī)械結(jié)構(gòu)強(qiáng)度與振動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室,西安 710049)

    一種提高飛行器結(jié)構(gòu)天地力學(xué)環(huán)境地面試驗(yàn)有效性的方法及其應(yīng)用

    閻桂榮,董龍雷,宋利強(qiáng)

    (西安交通大學(xué) 機(jī)械結(jié)構(gòu)強(qiáng)度與振動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室,西安 710049)

    目的 研究提高飛行器結(jié)構(gòu)地面試驗(yàn)有效性的途徑。方法 計(jì)算同一被試件結(jié)構(gòu)在飛行狀態(tài)和地面試驗(yàn)狀態(tài)下的有限元模型,測(cè)量地面試驗(yàn)狀態(tài)下的模態(tài)以驗(yàn)證有限元模型的正確性;計(jì)算各特征點(diǎn)(也可以是遙測(cè)點(diǎn))在天地狀態(tài)下的響應(yīng),用機(jī)器學(xué)習(xí)法獲取各特征點(diǎn)的映射關(guān)系模型;基于該模型由飛行點(diǎn)響應(yīng)(或遙測(cè)數(shù)據(jù))確定出地面試驗(yàn)件對(duì)應(yīng)點(diǎn)的響應(yīng),并用載荷反求法得到它們的等效載荷;最終確定施加在試驗(yàn)系統(tǒng)上的載荷。結(jié)果 以細(xì)長(zhǎng)體結(jié)構(gòu)為例,所得到由其組成試驗(yàn)系統(tǒng)的有限元模型與實(shí)測(cè)模型之間的固有頻率最大相對(duì)誤差為 6.76%,利用映射關(guān)系模型預(yù)測(cè)出對(duì)應(yīng)點(diǎn)在飛行狀態(tài)下的振動(dòng)響應(yīng)。確定了飛行狀態(tài)下結(jié)構(gòu)響應(yīng)的特征點(diǎn),由地面試驗(yàn)系統(tǒng)所對(duì)應(yīng)的響應(yīng)點(diǎn)反推出應(yīng)施加的載荷為60 N。結(jié)論 利用天地?cái)?shù)值計(jì)算-地面試驗(yàn)驗(yàn)證聯(lián)合法,無(wú)需在地面試驗(yàn)狀態(tài)下刻意模擬飛行狀態(tài)的邊界條件,確定出所需要施加的載荷,從而提高了飛行器地面試驗(yàn)的有效性。

    振動(dòng)試驗(yàn);映射關(guān)系模型;有限元分析;試驗(yàn)有效性

    地面振動(dòng)試驗(yàn)是飛行器研制全過(guò)程中的重要環(huán)節(jié),地面振動(dòng)試驗(yàn)的有效性不僅直接關(guān)系到飛行器的工作可靠性及性能,而且對(duì)飛行器的精細(xì)化設(shè)計(jì)具有更重要的意義,但如何提高地面試驗(yàn)的有效性是一個(gè)值得深思的問(wèn)題。提高地面試驗(yàn)的有效性,需要解決兩個(gè)關(guān)鍵問(wèn)題:要制定正確的地面試驗(yàn)規(guī)范或條件;要能夠準(zhǔn)確地實(shí)現(xiàn)試驗(yàn)規(guī)范,即給出試驗(yàn)方法。對(duì)于地面試驗(yàn)規(guī)范的制定,最常用的方法有兩類:總體設(shè)計(jì)部門多用類比法、歸納法、包絡(luò)法、外推法等,隨著數(shù)值計(jì)算的日趨完善,有限元法也受到青睞;另一類的依據(jù)就是國(guó)內(nèi)外制定的環(huán)境振動(dòng)試驗(yàn)標(biāo)準(zhǔn),試驗(yàn)標(biāo)準(zhǔn)只是規(guī)定了試驗(yàn)條件,并未告知這些試驗(yàn)條件如何得到的。準(zhǔn)確地實(shí)現(xiàn)試驗(yàn)規(guī)范就是按照所制定的試驗(yàn)規(guī)范給出試驗(yàn)條件并實(shí)現(xiàn)它。當(dāng)一個(gè)結(jié)構(gòu)設(shè)計(jì)完成后,其固有特性就確定了,其結(jié)構(gòu)的振動(dòng)響應(yīng)還與邊界條件、載荷有關(guān)。載荷是外激勵(lì),它通過(guò)結(jié)構(gòu)的固有特性影響著振動(dòng)響應(yīng)。邊界條件一旦確定,就影響著系統(tǒng)的固有特性。邊界條件參數(shù)的影響卻難以模擬,特別在試驗(yàn)臺(tái)上安裝的零部件、分系統(tǒng),很難模擬出它們?cè)陲w行器整體結(jié)構(gòu)上的聯(lián)接方式。為了考核零部件、分系統(tǒng)在地面試驗(yàn)裝置上的響應(yīng),通常以結(jié)構(gòu)的局部響應(yīng)等效為依據(jù)。由于試驗(yàn)系統(tǒng)本身的動(dòng)力學(xué)特性、及試件在試驗(yàn)系統(tǒng)中的安裝條件所限,結(jié)構(gòu)局部響應(yīng)等效多是依靠試驗(yàn)系統(tǒng)控制器的負(fù)反饋實(shí)現(xiàn)。雖然被控點(diǎn)可以滿足規(guī)范要求,但未控點(diǎn)可能會(huì)出現(xiàn)過(guò)試驗(yàn)或欠試驗(yàn)狀況,特別是對(duì)大型結(jié)構(gòu),甚至?xí)霈F(xiàn)試驗(yàn)難以進(jìn)行而不得不修改試驗(yàn)條件的局面[1—3]。

    細(xì)長(zhǎng)體結(jié)構(gòu)是火箭、導(dǎo)彈等常用的結(jié)構(gòu),文中將以其為對(duì)象介紹一種飛行器結(jié)構(gòu)天地?cái)?shù)值計(jì)算-地面試驗(yàn)驗(yàn)證聯(lián)合法,以提高地面試驗(yàn)的有效性。

    1 理論基礎(chǔ)

    1.1 映射關(guān)系模型

    映射關(guān)系模型是應(yīng)用天地?cái)?shù)值計(jì)算-地面試驗(yàn)驗(yàn)證聯(lián)合法的基礎(chǔ)。映射關(guān)系就是將輸入空間的變量映射到高維的特征空間,在高維空間中尋求它們的關(guān)系,利用映射關(guān)系獲取的模型稱為映射關(guān)系模型[4]。

    文獻(xiàn)[4]中還引入了一個(gè)時(shí)域決策函數(shù) ψi(t),而 ψ(t)為時(shí)域決策函數(shù) ψi(t)的集合,ψ(t)描述了相同載荷下系統(tǒng)1和系統(tǒng)2各對(duì)應(yīng)點(diǎn)時(shí)域響應(yīng)之間的關(guān)系。

    根據(jù)文獻(xiàn)[4]的證明和對(duì)基于時(shí)域的第一映射關(guān)系模型和基于頻域的第二映射關(guān)系模型定義,發(fā)展了天地?cái)?shù)值計(jì)算-地面試驗(yàn)驗(yàn)證聯(lián)合法(簡(jiǎn)稱聯(lián)合法)。所謂“天”是指飛行器被試部件(或整體)結(jié)構(gòu)的飛行狀態(tài),所謂“地”是指飛行器被試部件(或整體)結(jié)構(gòu)在試驗(yàn)系統(tǒng)中的狀態(tài),需要指出的是無(wú)須對(duì)結(jié)構(gòu)在試驗(yàn)系統(tǒng)中的邊界條件提出模擬飛行狀態(tài)的要求。

    1.2 實(shí)施步驟

    以映射關(guān)系為理論基礎(chǔ),以機(jī)器學(xué)習(xí)為基本算法,實(shí)現(xiàn)天地?cái)?shù)值計(jì)算-地面試驗(yàn)驗(yàn)證聯(lián)合法的主要步驟如下。

    1) 建立被試部件結(jié)構(gòu)在飛行狀態(tài)下(簡(jiǎn)稱系統(tǒng)1)的有限元模型,建模時(shí)應(yīng)考慮總體結(jié)構(gòu)前后連接件對(duì)該限元模型的影響。

    2) 建立被試部件結(jié)構(gòu)在地面試驗(yàn)狀態(tài)下(簡(jiǎn)稱系統(tǒng)2)的有限元模型,此時(shí)也應(yīng)考慮試驗(yàn)系統(tǒng)通過(guò)各聯(lián)接件對(duì)被試件有限元模型的影響。

    3) 分別對(duì)系統(tǒng)1、系統(tǒng)2中被試部件結(jié)構(gòu)各點(diǎn)的響應(yīng)進(jìn)行聚類,以確定并關(guān)注其響應(yīng)特征點(diǎn)的位置。

    4) 建立被試部件結(jié)構(gòu)的天地響應(yīng)的映射關(guān)系模型。具體做法是給系統(tǒng)1和系統(tǒng)2的有限元模型施加相同的載荷,為了能反映被試件的特性,載荷的類型應(yīng)為平譜,且參加機(jī)器學(xué)習(xí)的訓(xùn)練樣本載荷的涵蓋范圍足夠?qū)挕?/p>

    5) 試驗(yàn)驗(yàn)證。其目的為驗(yàn)證結(jié)構(gòu)在系統(tǒng) 2中所建有限元模型的正確性;從系統(tǒng)2實(shí)測(cè)結(jié)構(gòu)的振動(dòng)響應(yīng)預(yù)測(cè)出結(jié)構(gòu)在系統(tǒng)1的響應(yīng)。

    6) 確定施加在系統(tǒng)2上的載荷。具體做法是:由步驟3)確定系統(tǒng)1聚類的特征點(diǎn)和步驟4)取得的與特征點(diǎn)對(duì)應(yīng)的系統(tǒng) 2上結(jié)構(gòu)的響應(yīng)樣本和模型,反求出一組被試件上的載荷(u1,u2,…,un)[5—8];通過(guò)映射關(guān)系模型和系統(tǒng)1測(cè)得的遙測(cè)信號(hào),確定相應(yīng)的載荷ui;利用該點(diǎn)的載荷ui與實(shí)際加載點(diǎn)之間的傳遞函數(shù)求出實(shí)際加載點(diǎn)的載荷。

    7) 假設(shè)取m個(gè)特征點(diǎn),可得到一組數(shù)據(jù)(um1,um2,…,umm),選取該組的最大值或程控法作為實(shí)際施加載荷,完成地面考核試驗(yàn)。

    2 試驗(yàn)系統(tǒng)研究

    以一個(gè)細(xì)長(zhǎng)體結(jié)構(gòu)為實(shí)例,介紹天地?cái)?shù)值計(jì)算-地面試驗(yàn)驗(yàn)證聯(lián)合法的具體應(yīng)用。

    2.1 試驗(yàn)系統(tǒng)設(shè)計(jì)

    圍繞實(shí)現(xiàn)天地?cái)?shù)值計(jì)算-地面試驗(yàn)驗(yàn)證聯(lián)合法的目標(biāo),設(shè)計(jì)了一個(gè)細(xì)長(zhǎng)體結(jié)構(gòu)試驗(yàn)件,該試驗(yàn)件由圓柱殼、截錐殼、二次曲面殼、后蓋板組成,后蓋板和圓柱殼、圓柱殼和截錐殼、截錐殼和二次曲面殼的連接方式都是內(nèi)外螺紋緊配合。細(xì)長(zhǎng)體結(jié)構(gòu)試驗(yàn)件借助于兩個(gè)支座安裝于底板上,形成地面試驗(yàn)系統(tǒng),如圖1所示。

    圖1 試驗(yàn)系統(tǒng)Fig.1 Test system

    試驗(yàn)件尾部安裝了一個(gè)振動(dòng)臺(tái)作為外激勵(lì)。系統(tǒng)還包括傳感器-測(cè)試分析系統(tǒng),用于完成試驗(yàn)建模、結(jié)構(gòu)響應(yīng)測(cè)試。

    2.2 系統(tǒng)建模

    利用子結(jié)構(gòu)法測(cè)取整體結(jié)構(gòu)的模態(tài),完成系統(tǒng)建模。測(cè)點(diǎn)布置、測(cè)點(diǎn)數(shù)量的選擇基于下述原則:能夠明確顯示在試驗(yàn)頻段內(nèi)所有模態(tài)的基本特征及相互間的區(qū)別,保證關(guān)心的結(jié)構(gòu)點(diǎn)都在所選的測(cè)點(diǎn)之中;為提高信噪比,測(cè)點(diǎn)不應(yīng)選在各階振型節(jié)點(diǎn)附近。根據(jù)以上原則,在試驗(yàn)對(duì)象上選擇了138個(gè)測(cè)試點(diǎn),為保證不遺漏重要模態(tài),在圓柱殼上選擇7個(gè)激勵(lì)點(diǎn),在每個(gè)支架上選擇一個(gè)激勵(lì)點(diǎn)。拾取各激勵(lì)點(diǎn)錘擊激勵(lì)下+x,+y,+z方向的響應(yīng)信號(hào),并計(jì)算出所有的頻率響應(yīng)函數(shù)FRF。用LSCF法進(jìn)行參數(shù)識(shí)別,模態(tài)識(shí)別穩(wěn)態(tài)圖如2所示,分析帶寬取為10~512 Hz。

    圖2 平臺(tái)總頻響的穩(wěn)態(tài)圖Fig.2 The steady-state frequency response diagram of the system

    結(jié)合穩(wěn)態(tài)圖、各敲擊點(diǎn)數(shù)據(jù)的對(duì)比和MAC表及振型分析,通過(guò)對(duì)參數(shù)識(shí)別的結(jié)果進(jìn)行比較可知,不同方向、不同點(diǎn)激勵(lì)以及不同方向頻率響應(yīng)分析結(jié)果基本一致,綜合分析后所得平臺(tái)模態(tài)描述見表1,模態(tài)振型如圖3所示。為保證試驗(yàn)?zāi)P偷挠行院驼_性,用互易性檢驗(yàn)、相干性檢驗(yàn)、MAC值檢驗(yàn)等方法對(duì)試驗(yàn)結(jié)果進(jìn)行了檢驗(yàn)。

    表1 平臺(tái)的固有頻率Table 1 Natural frequency of the system

    圖3 平臺(tái)整體前六階振型Fig.3 The preceding six mode shapes of the system

    3 數(shù)值仿真

    3.1 無(wú)邊界條件約束

    為模擬飛行狀態(tài),建立了被試件的有限元模型如圖4所示,模態(tài)描述見表2,模態(tài)振型如圖5所示。

    圖4 試驗(yàn)件有限元模型Fig.4 The finite element model of the specimen

    表2 自由狀態(tài)下試件的固有頻率Table 2 The natural frequency of the specimen in free state

    圖5 自由狀態(tài)下的模態(tài)振型Fig.5 The mode shapes of the sepcimen in free state

    3.2 地面試驗(yàn)系統(tǒng)

    利用試件、支架的有限元模型,采用節(jié)點(diǎn)對(duì)應(yīng)重合的連接方式建立了地面試驗(yàn)系統(tǒng)的有限元模型,建立試驗(yàn)臺(tái)整體的有限元模型如圖6所示,有限元計(jì)算與實(shí)測(cè)試驗(yàn)系統(tǒng)的固有頻率比較見表3,有限元分析得到的整體前四階振型如圖所示。

    由于試驗(yàn)中前兩階的剛體運(yùn)動(dòng)振型無(wú)法在有限元分析中得到,其余振型都能在有限元分析中得到,且有限元分析得到的固有頻率與實(shí)驗(yàn)值相比,相對(duì)誤差最大為6.76%,滿足工程要求。由于結(jié)構(gòu)較大,模態(tài)實(shí)驗(yàn)測(cè)點(diǎn)分布相對(duì)稀疏,實(shí)驗(yàn)只能得到前六階振型,若要得到更多振型,應(yīng)對(duì)測(cè)點(diǎn)布局加密。

    圖6 地面試驗(yàn)系統(tǒng)的有限元模型Fig.6 The finite element model of the ground test system

    表3 地面試驗(yàn)條件下試件實(shí)測(cè)與有限元計(jì)算固有頻率的比較Table 3 Comparison of the natural frequency of the test specimen and the finite element calculation under the ground test condition

    圖7 平臺(tái)有限元分析前四階振型Fig.7 The preceding four mode shapes of the system in finite element analysis

    4 響應(yīng)預(yù)示與試驗(yàn)載荷的獲取

    在模擬件尾端施加振動(dòng)激勵(lì),根據(jù)結(jié)構(gòu)響應(yīng)進(jìn)行聚類分析,確定了模擬件端部某點(diǎn)作為預(yù)示點(diǎn),如圖8所示。針對(duì)該模擬件結(jié)構(gòu),構(gòu)造振動(dòng)激勵(lì)數(shù)據(jù)集,某振動(dòng)激勵(lì)譜型如圖9所示,計(jì)算模擬件預(yù)示點(diǎn)在飛行狀態(tài)及地面試驗(yàn)狀態(tài)下的振動(dòng)加速度響應(yīng)。

    圖8 基于聚類分析確定的預(yù)示測(cè)點(diǎn)Fig.8 The prediction point based on clustering analysis

    圖9 振動(dòng)激勵(lì)力譜型Fig.9 The spectrum of the vibration excition force

    根據(jù)振動(dòng)響應(yīng)數(shù)據(jù)及構(gòu)造學(xué)習(xí)樣本,采用機(jī)器學(xué)習(xí)的方法建立映射關(guān)系模型。然后,依據(jù)地面試驗(yàn)設(shè)計(jì)所得到的振動(dòng)響應(yīng)作為輸入樣本,輸入到映射關(guān)系模型即可得到飛行狀態(tài)下的振動(dòng)響應(yīng),并與飛行狀態(tài)計(jì)算所得的振動(dòng)響應(yīng)進(jìn)行比較,實(shí)現(xiàn)振動(dòng)響應(yīng)映射預(yù)示方法的仿真驗(yàn)證。

    基于30組振動(dòng)激勵(lì)數(shù)據(jù),選取25組進(jìn)行映射關(guān)系模型的建立,其他未參與學(xué)習(xí)的5組數(shù)據(jù)進(jìn)行預(yù)示與驗(yàn)證。將試驗(yàn)狀態(tài)的響應(yīng)數(shù)據(jù)輸入到映射關(guān)系模型,得到飛行狀態(tài)下對(duì)應(yīng)測(cè)點(diǎn)的振動(dòng)響應(yīng)預(yù)示結(jié)果,并與計(jì)算得到的飛行狀態(tài)下的振動(dòng)響應(yīng)進(jìn)行對(duì)比,結(jié)果如圖10所示,預(yù)示結(jié)果與計(jì)算結(jié)果基本吻合。

    圖10 由試驗(yàn)狀態(tài)預(yù)示飛行狀態(tài)的響應(yīng)Fig.10 The response of the flight state is indicated by the test condition

    另外,若已知飛行狀態(tài)的遙測(cè)數(shù)據(jù),如圖 11所示,將該遙測(cè)振動(dòng)響應(yīng)作為工作樣本,輸入到映射關(guān)系模型,得到地面試驗(yàn)狀態(tài)的振動(dòng)響應(yīng),如圖12所示。那么,基于載荷反求可確定施加在模擬件尾部的激勵(lì)力譜型,如圖13所示。

    圖11 飛行狀態(tài)振動(dòng)響應(yīng)遙測(cè)數(shù)據(jù)Fig.11 The vibration response telemetry data in flight state

    圖12 根據(jù)天地映射確定的地面狀態(tài)振動(dòng)響應(yīng)數(shù)據(jù)Fig.12 The vbiration response data in ground state determined by mapping relation

    圖13 基于載荷反求確定的地面試驗(yàn)激勵(lì)力譜型Fig.13 The ground test excitation spectrum based on load identification

    映射關(guān)系模型代表了被研究結(jié)構(gòu)的固有特性,當(dāng)建立該模型時(shí)需要利用受兩種不同邊界約束的被研究結(jié)構(gòu)的響應(yīng),要求建立的模型正確,并給系統(tǒng)1和系統(tǒng)2施加相同的載荷。為了能充分反映系統(tǒng)1和系統(tǒng)2的特征,施加到系統(tǒng)1和系統(tǒng)2模型上的載荷的類型應(yīng)為平譜,且載荷的范圍足夠大。

    5 結(jié)語(yǔ)

    提出天地?cái)?shù)值計(jì)算-地面試驗(yàn)驗(yàn)證聯(lián)合法,通過(guò)映射關(guān)系模型將飛行和地面試驗(yàn)環(huán)境下結(jié)構(gòu)的振動(dòng)響應(yīng)有機(jī)地結(jié)合起來(lái),在不需要模擬飛行狀態(tài)邊界條件的情況下,完成地面模擬試驗(yàn)。以細(xì)長(zhǎng)體結(jié)構(gòu)為例,實(shí)現(xiàn)了計(jì)算和試驗(yàn)的全過(guò)程,證明了該方法的可實(shí)施性。

    [1] GJB 150A—2009, 軍用裝備實(shí)驗(yàn)室環(huán)境試驗(yàn)方法[S]. GJB 150A—2009, Laboratory Environmental Test Methods for Military Materiel[S].

    [2] MIL-STD-810G, Environmental Engineering Considera-tions and Laboratory Tests[S].

    [3] MIL-STD-1540C, Test Requirements Launch Upper-stage Space[S].

    [4] 閆桂榮, 董龍雷, 喻磊. 基于機(jī)器學(xué)習(xí)動(dòng)力學(xué)環(huán)境的預(yù)示方法[J]. 應(yīng)用力學(xué)學(xué)報(bào), 2013, 30(1): 13—18. YAN Gui-rong, DONG Long-lei, YU Lei. A New Dynamical Environment Predication Method Based on Machine Learning[J]. Chinese Journal of Applied Mechanics, 2013, 30(1): 13—18.

    [5] 毛文濤. 支持向量回歸機(jī)模型選擇研究及在綜合力學(xué)環(huán)境預(yù)示中的應(yīng)用[D]. 西安: 西安交通大學(xué), 2011. MAO Wen-tao. Research on Model Selection for Support Vector Regression and Its Application in Prediction of Combined Dynamics Environment[D]. Xi'an: Xi'an Jiaotong University, 2011.

    [6] MAO Wen-tao, HU Di-ke, YAN Gu-irong. A New SVM Regression Approach for Mechanical Load Identification[J]. International Journal of Applied Electromagnetics and Mechanics, 2010, 33(3/4): 1001—1008.

    [7] MAO Wen-tao, YAN Gui-rong, DONG Long-lei. Weighted Solution Path Algorithm of Support Vector Regression Based on Heuristic Weight-setting Optimization[J]. Neurocomputing, 2009, 73(1—3): 495—505.

    [8] 王成, 焦吉祥, 黃翔龍, 等. 基于一元線性回歸模型的多源頻域載荷識(shí)別[C]// 第八屆全國(guó)隨機(jī)振動(dòng)理論與應(yīng)用學(xué)術(shù)會(huì)議暨第一屆全國(guó)隨機(jī)動(dòng)力學(xué)學(xué)術(shù)會(huì)議論文集.成都: 西南交通大學(xué), 2012: 848—854. WANG Cheng, JIAO Ji-xiang, HUANG Xiang-long, et al. A Multi-sources Load Identification Algorithm in Frequency Domain Based on Linear Regression[C]// International Conference on Theory Application of Random Vibration. Chengdu: Southwest Jiaotong University, 2012: 848—854.

    [9] VAPNIK V N. Estimation of Dependences Based on Empirical Data[M]. Berlin: Springer, 1982.

    [10] VAPNIK V N. The Nature of Statistical Learning Theory[M]. Berlin: Springer, 1995.

    [11] CHERKASSKY V, MULIER F. Learning from Data: Concepts, Theory and Methods[M]. New York: John Viley& Sons, 1997.

    [12] MJOLSNESS E, DECOSTE D. Machine Learning for Science: State of the Art and Future Prospects[J]. Science, 2001, 293(14): 2051—2055.

    [13] 張學(xué)工. 關(guān)于統(tǒng)計(jì)學(xué)習(xí)理論與支持向量機(jī)[J]. 自動(dòng)化學(xué)報(bào), 2000, 26(1): 32—42. ZHANG Xue-gong. Introduction to Statistical Learning Theory and Support Vector Machine[J]. Chinese Journal of Automation, 2000, 26(1): 32—42.

    [14] 鄧乃揚(yáng), 田英杰. 數(shù)據(jù)挖掘中的新方法——支持向量機(jī)[M]. 北京: 科學(xué)出版社, 2006. DENG Nai-yang, TIAN Ying-jie. New Methods in Data Mining—Support Vector Maching[M]. Beijing: Science Press, 2006.

    [15] 王彥飛. 反演問(wèn)題的計(jì)算方法及其應(yīng)用[M]. 北京: 高等教育出版社, 2007. WANG Yan-fei. Computational Methods Inverse Problems and Their Application[M]. Beijing: Higher Education Press, 2007.

    A Method and Its Application for Improving the Validity of Ground Testing of Mechanical Environment's Effects on Aircraft Structure

    YAN Gui-rong, DONG Long-lei, SONG Li-qiang
    (State Key Laboratory for Strength and Vibration, Xi'an Jiaotong University, Xi'an 710049, China)

    Objective To study the way of improving the validity of the aircraft structure ground test. Methods The finite element models of the specimen structure in flight and ground test conditions were calculated. Then the test model of the specimen structure in ground testing was measured to verify the validity of the finite element model. The vibration response of feature points in flight and ground testing was calculated. The mapping relation model of the feature points was calculated by machine learning methods. Then the response of each point in ground testing was determined by flight point response (or telemetry data) and their equivalent load was obtained by load reverse method. Finally, the load imposed on the test system was determined. Results With a slender body structure as an example, the maximum relative error between the finite element model and the measured model of the test system was 6.76%. The vibration response of the feature points in flight state could be predicted via mapping relation model. The vibration response of the cor-responding ground testing points could be determined. The load to be applied was set at 60 N by reverse deliberation according to response points of ground test system. Conclusion By the unified method of the flight-ground numerical calculation and the ground testing, the load to be applied in the ground test can be determined without deliberately simulating the boundary conditions in flight state, which improves the validity of the aircraft ground testing.

    vibration test; mapping relation model; finite element analysis; test validity

    2016-08-12;Revised:2016-08-29

    10.7643/ issn.1672-9242.2016.05.002

    TS206

    A

    1672-9242(2016)05-0010-07

    2016-08-12;

    2016-08-29

    國(guó)家自然科學(xué)基金項(xiàng)目(10276032)

    fund program:National Natural Science Foundation of China(10276032)

    閻桂榮(1942—),女,陜西西安人,教授,主要研究方向?yàn)閺?fù)合動(dòng)力學(xué)環(huán)境仿真與控制、振動(dòng)控制、現(xiàn)代力學(xué)測(cè)控理論與技術(shù)等。

    Biography:YAN Gui-rong (1942—), Female, from Xi'an, Shanxi, Professor, Research focus: environmentdynamic simulation & control,theory andtechnology of vibration measure & control.

    猜你喜歡
    振型載荷有限元
    關(guān)于模態(tài)綜合法的注記
    交通運(yùn)輸部海事局“新一代衛(wèi)星AIS驗(yàn)證載荷”成功發(fā)射
    水上消防(2022年2期)2022-07-22 08:45:00
    縱向激勵(lì)下大跨鋼桁拱橋高階振型效應(yīng)分析
    塔腿加過(guò)渡段輸電塔動(dòng)力特性分析
    滾轉(zhuǎn)機(jī)動(dòng)載荷減緩風(fēng)洞試驗(yàn)
    結(jié)構(gòu)振型幾何辨識(shí)及應(yīng)用研究
    山西建筑(2015年14期)2015-06-05 09:37:07
    磨削淬硬殘余應(yīng)力的有限元分析
    一種基于白噪聲響應(yīng)的隨機(jī)載荷譜識(shí)別方法
    底排藥受力載荷及其分布規(guī)律
    基于SolidWorks的吸嘴支撐臂有限元分析
    黄色毛片三级朝国网站| 国产xxxxx性猛交| 免费在线观看完整版高清| 香蕉丝袜av| 亚洲一卡2卡3卡4卡5卡精品中文| 丰满迷人的少妇在线观看| 精品久久久精品久久久| 新久久久久国产一级毛片| 亚洲av美国av| 最近最新免费中文字幕在线| 多毛熟女@视频| 一本大道久久a久久精品| 久久久久久久国产电影| 黄片小视频在线播放| 美女脱内裤让男人舔精品视频| 成人黄色视频免费在线看| 欧美另类一区| 亚洲欧美日韩另类电影网站| 男人添女人高潮全过程视频| 欧美激情高清一区二区三区| 成人国产av品久久久| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av美国av| 后天国语完整版免费观看| 成人三级做爰电影| 一边摸一边抽搐一进一出视频| 成人国产一区最新在线观看| 国产一卡二卡三卡精品| 日韩大码丰满熟妇| 日本猛色少妇xxxxx猛交久久| 亚洲欧美一区二区三区久久| 亚洲成人免费av在线播放| 日韩,欧美,国产一区二区三区| 亚洲精品久久久久久婷婷小说| av网站在线播放免费| 精品国产国语对白av| 狠狠精品人妻久久久久久综合| 久久久久精品人妻al黑| 欧美日韩福利视频一区二区| 精品少妇内射三级| 日本五十路高清| 成人免费观看视频高清| 国产97色在线日韩免费| 我要看黄色一级片免费的| 好男人电影高清在线观看| www.av在线官网国产| 国产色视频综合| 成人影院久久| 欧美 日韩 精品 国产| 亚洲国产精品999| 最近最新免费中文字幕在线| 淫妇啪啪啪对白视频 | 男女之事视频高清在线观看| 中文字幕色久视频| 午夜激情av网站| 欧美av亚洲av综合av国产av| av网站免费在线观看视频| 日韩中文字幕视频在线看片| 亚洲国产精品成人久久小说| 1024香蕉在线观看| 久久久久国内视频| 久久99热这里只频精品6学生| 亚洲国产欧美日韩在线播放| 人人妻人人澡人人爽人人夜夜| 亚洲人成电影免费在线| 女人爽到高潮嗷嗷叫在线视频| 日韩欧美国产一区二区入口| 一个人免费在线观看的高清视频 | 一边摸一边做爽爽视频免费| 欧美激情 高清一区二区三区| 久久久久国内视频| 日本av手机在线免费观看| 999久久久精品免费观看国产| 亚洲视频免费观看视频| 最近最新免费中文字幕在线| 两个人看的免费小视频| 久久九九热精品免费| 两个人看的免费小视频| 亚洲专区中文字幕在线| 国产精品av久久久久免费| 最新在线观看一区二区三区| 老司机影院成人| 亚洲欧美色中文字幕在线| 亚洲欧美精品自产自拍| 高清在线国产一区| 另类精品久久| 岛国毛片在线播放| 人人妻人人澡人人看| 午夜福利视频精品| 亚洲成人免费av在线播放| 肉色欧美久久久久久久蜜桃| 国产国语露脸激情在线看| 精品少妇久久久久久888优播| 不卡一级毛片| 久久女婷五月综合色啪小说| 各种免费的搞黄视频| 两性午夜刺激爽爽歪歪视频在线观看 | 人成视频在线观看免费观看| 正在播放国产对白刺激| 午夜精品国产一区二区电影| 视频区欧美日本亚洲| 国产国语露脸激情在线看| 国产成人a∨麻豆精品| 在线亚洲精品国产二区图片欧美| 国产精品一区二区在线观看99| 人人妻人人澡人人看| netflix在线观看网站| 精品人妻1区二区| 丝袜美腿诱惑在线| 可以免费在线观看a视频的电影网站| 国产在视频线精品| 日韩精品免费视频一区二区三区| 大香蕉久久成人网| 亚洲成人手机| 国产精品麻豆人妻色哟哟久久| 亚洲九九香蕉| 亚洲avbb在线观看| 侵犯人妻中文字幕一二三四区| 精品人妻在线不人妻| 日韩人妻精品一区2区三区| 一级毛片电影观看| 麻豆乱淫一区二区| 亚洲精品美女久久久久99蜜臀| 在线观看免费午夜福利视频| av片东京热男人的天堂| 九色亚洲精品在线播放| 王馨瑶露胸无遮挡在线观看| 老司机影院成人| 在线十欧美十亚洲十日本专区| 丰满饥渴人妻一区二区三| videos熟女内射| 宅男免费午夜| 黄色视频不卡| 狠狠狠狠99中文字幕| 久久九九热精品免费| 亚洲成国产人片在线观看| 91老司机精品| 高清欧美精品videossex| 蜜桃国产av成人99| 久久狼人影院| 日韩三级视频一区二区三区| 国产成人精品久久二区二区91| 国产亚洲欧美精品永久| 丰满饥渴人妻一区二区三| 捣出白浆h1v1| 国产精品亚洲av一区麻豆| 超碰成人久久| 午夜福利乱码中文字幕| 超色免费av| 亚洲精品av麻豆狂野| 99国产极品粉嫩在线观看| 肉色欧美久久久久久久蜜桃| 男女高潮啪啪啪动态图| 18禁观看日本| 欧美在线一区亚洲| 国产免费视频播放在线视频| 搡老岳熟女国产| 91av网站免费观看| 亚洲国产欧美在线一区| 欧美日韩精品网址| 精品人妻在线不人妻| 成年美女黄网站色视频大全免费| 波多野结衣一区麻豆| 国产精品一区二区免费欧美 | 亚洲精品粉嫩美女一区| 国产一区二区三区av在线| 午夜福利视频精品| 涩涩av久久男人的天堂| 女性生殖器流出的白浆| 极品人妻少妇av视频| 热99国产精品久久久久久7| 丝瓜视频免费看黄片| 国产伦理片在线播放av一区| 亚洲一区二区三区欧美精品| 一本—道久久a久久精品蜜桃钙片| 最近最新中文字幕大全免费视频| 精品少妇黑人巨大在线播放| 亚洲色图 男人天堂 中文字幕| 亚洲成人手机| a级片在线免费高清观看视频| a级毛片在线看网站| 99香蕉大伊视频| 50天的宝宝边吃奶边哭怎么回事| 女性被躁到高潮视频| 久久人人97超碰香蕉20202| 国产在线视频一区二区| 国产欧美日韩综合在线一区二区| 国产99久久九九免费精品| netflix在线观看网站| 乱人伦中国视频| 天天操日日干夜夜撸| 可以免费在线观看a视频的电影网站| 王馨瑶露胸无遮挡在线观看| 成年av动漫网址| 一本一本久久a久久精品综合妖精| 国产成人精品无人区| 午夜免费观看性视频| 高清视频免费观看一区二区| 人人妻人人澡人人看| 亚洲精品久久午夜乱码| 动漫黄色视频在线观看| 欧美激情极品国产一区二区三区| 久久久国产成人免费| 日韩免费高清中文字幕av| 亚洲精品一区蜜桃| 精品国产乱码久久久久久男人| 99热国产这里只有精品6| 91精品国产国语对白视频| 超碰成人久久| 人妻一区二区av| 国产黄色免费在线视频| 亚洲情色 制服丝袜| 日韩一区二区三区影片| 亚洲成人免费av在线播放| 日韩大码丰满熟妇| av天堂久久9| 一区二区三区激情视频| 亚洲精品美女久久久久99蜜臀| 操美女的视频在线观看| 日韩免费高清中文字幕av| 在线观看一区二区三区激情| 中文字幕人妻丝袜制服| 亚洲免费av在线视频| 亚洲精品第二区| 国产又色又爽无遮挡免| 久久久国产成人免费| 真人做人爱边吃奶动态| 老汉色∧v一级毛片| 免费女性裸体啪啪无遮挡网站| 桃花免费在线播放| 夜夜骑夜夜射夜夜干| 巨乳人妻的诱惑在线观看| 俄罗斯特黄特色一大片| 50天的宝宝边吃奶边哭怎么回事| 久久久久久久国产电影| 国产亚洲精品久久久久5区| 婷婷丁香在线五月| 狂野欧美激情性bbbbbb| 亚洲av日韩精品久久久久久密| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜91福利影院| 精品亚洲乱码少妇综合久久| 最近最新免费中文字幕在线| 亚洲avbb在线观看| 午夜91福利影院| 精品少妇内射三级| 久久久久国产精品人妻一区二区| 老汉色∧v一级毛片| 亚洲av成人不卡在线观看播放网 | 午夜激情av网站| 欧美亚洲日本最大视频资源| 免费少妇av软件| 国产伦理片在线播放av一区| 黄色视频,在线免费观看| 1024视频免费在线观看| 精品国产一区二区久久| 午夜成年电影在线免费观看| 窝窝影院91人妻| 在线十欧美十亚洲十日本专区| 精品国产一区二区三区四区第35| 欧美在线一区亚洲| svipshipincom国产片| 黑人巨大精品欧美一区二区mp4| 亚洲国产看品久久| 亚洲精品国产精品久久久不卡| 性色av乱码一区二区三区2| 国产成人精品无人区| 如日韩欧美国产精品一区二区三区| 国产国语露脸激情在线看| 欧美国产精品一级二级三级| xxxhd国产人妻xxx| 下体分泌物呈黄色| 伊人久久大香线蕉亚洲五| 欧美黑人精品巨大| 人人妻人人澡人人爽人人夜夜| 日韩大码丰满熟妇| 大码成人一级视频| 亚洲av片天天在线观看| 国产精品国产三级国产专区5o| 日本91视频免费播放| 日韩熟女老妇一区二区性免费视频| 亚洲av日韩精品久久久久久密| 中国美女看黄片| 久久精品久久久久久噜噜老黄| 成年美女黄网站色视频大全免费| 国产男女内射视频| 亚洲欧美激情在线| 国产福利在线免费观看视频| 国产日韩一区二区三区精品不卡| 国产高清视频在线播放一区 | 国产熟女午夜一区二区三区| 97人妻天天添夜夜摸| 涩涩av久久男人的天堂| 啦啦啦 在线观看视频| 日韩制服骚丝袜av| 人人澡人人妻人| 国产成人欧美| 久久久国产精品麻豆| 日韩制服骚丝袜av| 黄色 视频免费看| 亚洲人成77777在线视频| 国产高清videossex| 国产精品久久久久久精品古装| 黑丝袜美女国产一区| 成人18禁高潮啪啪吃奶动态图| 老熟妇仑乱视频hdxx| 亚洲av电影在线进入| 18禁裸乳无遮挡动漫免费视频| 水蜜桃什么品种好| 女人久久www免费人成看片| 国产人伦9x9x在线观看| 国产欧美日韩一区二区三区在线| 天天操日日干夜夜撸| 男女国产视频网站| 久久精品熟女亚洲av麻豆精品| 黄色视频,在线免费观看| 国产熟女午夜一区二区三区| 日韩免费高清中文字幕av| 亚洲国产看品久久| 波多野结衣一区麻豆| 两个人看的免费小视频| 一区二区日韩欧美中文字幕| 免费观看a级毛片全部| 亚洲国产欧美网| 日韩 欧美 亚洲 中文字幕| 男女下面插进去视频免费观看| 丝袜人妻中文字幕| 丁香六月欧美| 精品视频人人做人人爽| 99香蕉大伊视频| 精品国产国语对白av| 啦啦啦啦在线视频资源| 久久中文看片网| 熟女少妇亚洲综合色aaa.| 妹子高潮喷水视频| 老司机影院毛片| www.熟女人妻精品国产| 亚洲欧美日韩另类电影网站| 国产精品.久久久| 老司机午夜福利在线观看视频 | 久9热在线精品视频| 亚洲成人免费电影在线观看| 免费观看人在逋| 欧美另类一区| 精品人妻1区二区| 伊人亚洲综合成人网| 国产成+人综合+亚洲专区| 久久毛片免费看一区二区三区| 欧美精品一区二区免费开放| 丝袜美足系列| 精品视频人人做人人爽| 欧美人与性动交α欧美精品济南到| 又紧又爽又黄一区二区| 亚洲精品粉嫩美女一区| 夫妻午夜视频| 国产野战对白在线观看| 99热网站在线观看| 欧美成人午夜精品| 久久国产精品大桥未久av| 午夜成年电影在线免费观看| 久久久久久久国产电影| 别揉我奶头~嗯~啊~动态视频 | 久久久精品免费免费高清| 如日韩欧美国产精品一区二区三区| 麻豆乱淫一区二区| 90打野战视频偷拍视频| 美女扒开内裤让男人捅视频| 老司机在亚洲福利影院| 啦啦啦 在线观看视频| 可以免费在线观看a视频的电影网站| 亚洲成人免费电影在线观看| 俄罗斯特黄特色一大片| 久久久久精品国产欧美久久久 | 国产成人精品久久二区二区91| 91老司机精品| 日韩精品免费视频一区二区三区| 亚洲欧美一区二区三区黑人| 欧美精品人与动牲交sv欧美| 国产免费av片在线观看野外av| 久久久国产精品麻豆| 欧美亚洲 丝袜 人妻 在线| 亚洲精品国产一区二区精华液| 91九色精品人成在线观看| 老司机深夜福利视频在线观看 | 日韩 亚洲 欧美在线| h视频一区二区三区| 欧美激情高清一区二区三区| 十分钟在线观看高清视频www| 在线观看免费高清a一片| 天天躁夜夜躁狠狠躁躁| 美女高潮到喷水免费观看| 欧美日韩一级在线毛片| 真人做人爱边吃奶动态| 亚洲 国产 在线| 国产一区二区激情短视频 | 99久久99久久久精品蜜桃| 国产真人三级小视频在线观看| 69av精品久久久久久 | 精品国产乱码久久久久久男人| 一区二区三区乱码不卡18| 亚洲成人国产一区在线观看| 两人在一起打扑克的视频| 亚洲色图综合在线观看| 久久久久精品人妻al黑| 法律面前人人平等表现在哪些方面 | 99久久人妻综合| 欧美激情极品国产一区二区三区| 18在线观看网站| 国产人伦9x9x在线观看| 国产成人系列免费观看| 欧美人与性动交α欧美精品济南到| www日本在线高清视频| 国产欧美日韩精品亚洲av| 国产精品一区二区免费欧美 | 国产av又大| 老司机深夜福利视频在线观看 | 大片免费播放器 马上看| 人妻 亚洲 视频| bbb黄色大片| 亚洲自偷自拍图片 自拍| 亚洲精品美女久久av网站| 亚洲欧洲日产国产| 极品少妇高潮喷水抽搐| 麻豆国产av国片精品| 久久久久久人人人人人| 亚洲精品久久久久久婷婷小说| 无限看片的www在线观看| 1024香蕉在线观看| 精品国内亚洲2022精品成人 | 99香蕉大伊视频| 777久久人妻少妇嫩草av网站| 国产在线视频一区二区| 97人妻天天添夜夜摸| 高清视频免费观看一区二区| 国产高清videossex| 人妻一区二区av| 国产激情久久老熟女| 在线观看一区二区三区激情| 18禁黄网站禁片午夜丰满| av天堂在线播放| 久久香蕉激情| 精品亚洲成a人片在线观看| 后天国语完整版免费观看| www日本在线高清视频| 日韩熟女老妇一区二区性免费视频| 精品熟女少妇八av免费久了| 女人高潮潮喷娇喘18禁视频| 国产成人av教育| 人人妻,人人澡人人爽秒播| 日韩有码中文字幕| 麻豆乱淫一区二区| 在线观看免费日韩欧美大片| a级毛片黄视频| 这个男人来自地球电影免费观看| 欧美性长视频在线观看| 成人免费观看视频高清| 国产成+人综合+亚洲专区| 成人亚洲精品一区在线观看| 丁香六月欧美| 亚洲视频免费观看视频| 一二三四在线观看免费中文在| 久久女婷五月综合色啪小说| 精品一区二区三区av网在线观看 | 视频区图区小说| 一区二区三区乱码不卡18| 日韩精品免费视频一区二区三区| 视频在线观看一区二区三区| av天堂久久9| 成人国产一区最新在线观看| 久久人妻福利社区极品人妻图片| 午夜久久久在线观看| 久久人人爽人人片av| 人妻 亚洲 视频| 两性夫妻黄色片| 精品少妇久久久久久888优播| 下体分泌物呈黄色| 亚洲,欧美精品.| 久久香蕉激情| 在线 av 中文字幕| 热99久久久久精品小说推荐| 99精品久久久久人妻精品| 国产免费一区二区三区四区乱码| 国产又爽黄色视频| 国产av又大| 国产亚洲av高清不卡| 国产成人精品久久二区二区91| 精品高清国产在线一区| 丝袜美腿诱惑在线| 国产成人精品久久二区二区免费| 欧美激情极品国产一区二区三区| 91老司机精品| 99精国产麻豆久久婷婷| 午夜福利视频精品| 欧美日本中文国产一区发布| 国产人伦9x9x在线观看| 大香蕉久久网| 9色porny在线观看| 国产一卡二卡三卡精品| 久久 成人 亚洲| 精品乱码久久久久久99久播| 精品一区在线观看国产| 日本精品一区二区三区蜜桃| 色视频在线一区二区三区| 他把我摸到了高潮在线观看 | 如日韩欧美国产精品一区二区三区| av国产精品久久久久影院| 精品人妻在线不人妻| 亚洲综合色网址| 亚洲精品第二区| 成人黄色视频免费在线看| 热re99久久精品国产66热6| 亚洲黑人精品在线| 久久久久久免费高清国产稀缺| 又大又爽又粗| 久久国产精品影院| 色94色欧美一区二区| 丰满少妇做爰视频| 精品久久久久久久毛片微露脸 | 免费在线观看视频国产中文字幕亚洲 | 国产成人精品无人区| 青春草亚洲视频在线观看| 国产视频一区二区在线看| 一二三四在线观看免费中文在| 男女床上黄色一级片免费看| 中文字幕人妻熟女乱码| 人妻人人澡人人爽人人| 两个人免费观看高清视频| 国产成人影院久久av| 色播在线永久视频| 少妇精品久久久久久久| 男女下面插进去视频免费观看| 两个人看的免费小视频| 亚洲精品一区蜜桃| 亚洲国产精品一区二区三区在线| 国产日韩一区二区三区精品不卡| 18禁观看日本| 欧美精品人与动牲交sv欧美| 国产1区2区3区精品| 日韩欧美一区二区三区在线观看 | 王馨瑶露胸无遮挡在线观看| 亚洲精品美女久久久久99蜜臀| 大片免费播放器 马上看| 人人妻人人澡人人爽人人夜夜| 久久国产亚洲av麻豆专区| 成年人免费黄色播放视频| 99久久99久久久精品蜜桃| 午夜福利在线免费观看网站| 亚洲国产欧美网| 国产免费现黄频在线看| 一边摸一边抽搐一进一出视频| 中文字幕人妻丝袜一区二区| 中文字幕最新亚洲高清| 精品亚洲成国产av| 最近最新免费中文字幕在线| 侵犯人妻中文字幕一二三四区| 69av精品久久久久久 | 色播在线永久视频| 母亲3免费完整高清在线观看| 老司机影院成人| 满18在线观看网站| 成年av动漫网址| 黄色片一级片一级黄色片| 久久久久久免费高清国产稀缺| 日本a在线网址| 老司机福利观看| 他把我摸到了高潮在线观看 | 亚洲精品国产一区二区精华液| 亚洲一区中文字幕在线| 亚洲成av片中文字幕在线观看| 国产精品 国内视频| 黄色毛片三级朝国网站| 成年动漫av网址| 国产精品秋霞免费鲁丝片| 一级黄色大片毛片| 久久国产精品大桥未久av| 天堂8中文在线网| 亚洲全国av大片| 51午夜福利影视在线观看| 亚洲欧美日韩高清在线视频 | 丝袜美腿诱惑在线| 极品人妻少妇av视频| 国产日韩欧美视频二区| 日韩熟女老妇一区二区性免费视频| 91精品国产国语对白视频| 久久国产精品大桥未久av| 亚洲三区欧美一区| 国产一区二区三区av在线| 美女视频免费永久观看网站| 亚洲欧洲日产国产| 国产不卡av网站在线观看| 两个人看的免费小视频| 国产成人精品久久二区二区91| 97人妻天天添夜夜摸| 亚洲精品日韩在线中文字幕| 国产日韩欧美亚洲二区| 中文字幕色久视频| av欧美777| 久久久久久久久免费视频了| 国产xxxxx性猛交| 国产又色又爽无遮挡免| 日韩三级视频一区二区三区| 免费在线观看黄色视频的| 亚洲欧洲精品一区二区精品久久久| 久久人人97超碰香蕉20202| 久久久国产精品麻豆| 国产激情久久老熟女| 国产亚洲av片在线观看秒播厂| 一区二区三区精品91| 啦啦啦啦在线视频资源| 久久精品人人爽人人爽视色| 高清黄色对白视频在线免费看| e午夜精品久久久久久久| 欧美大码av| 国产精品国产av在线观看| 欧美少妇被猛烈插入视频|