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

    行波管中高溫平板結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)預(yù)報(bào)

    2017-10-23 07:25:56田艷妮王建民秦朝紅華宏星
    噪聲與振動(dòng)控制 2017年5期
    關(guān)鍵詞:振型固有頻率平板

    田艷妮,王建民,秦朝紅,程 昊,華宏星

    (1.上海交通大學(xué) 振動(dòng)、沖擊、噪聲研究所 機(jī)械系統(tǒng)與振動(dòng)國家重點(diǎn)實(shí)驗(yàn)室,上海 200240;2.北京環(huán)境強(qiáng)度研究所 可靠性與環(huán)境工程技術(shù)重點(diǎn)實(shí)驗(yàn)室,北京 100076)

    行波管中高溫平板結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)預(yù)報(bào)

    田艷妮1,王建民2,秦朝紅2,程 昊2,華宏星1

    (1.上海交通大學(xué) 振動(dòng)、沖擊、噪聲研究所 機(jī)械系統(tǒng)與振動(dòng)國家重點(diǎn)實(shí)驗(yàn)室,上海 200240;2.北京環(huán)境強(qiáng)度研究所 可靠性與環(huán)境工程技術(shù)重點(diǎn)實(shí)驗(yàn)室,北京 100076)

    研究熱噪聲復(fù)合載荷作用下結(jié)構(gòu)的動(dòng)力學(xué)響應(yīng)對高超聲速飛行器結(jié)構(gòu)的設(shè)計(jì)具有重要的指導(dǎo)意義,實(shí)驗(yàn)室多采用高溫行波管進(jìn)行熱噪聲復(fù)合試驗(yàn)研究,本文針對行波管中高溫平板結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)預(yù)示方法進(jìn)行研究。選取四邊固支C/SiC材料平板為研究對象,首先研究熱效應(yīng)對平板結(jié)構(gòu)固有模態(tài)特性的影響,然后在平板表面施加一維行波聲載荷模擬行波管中聲場,利用解析方法求解結(jié)構(gòu)在熱噪聲復(fù)合載荷下的動(dòng)力學(xué)響應(yīng)。結(jié)果表明,在加熱過程中平板固有頻率先減小至極小值,而后逐漸回升。在固有頻率下降階段,材料彈性模量減小和熱應(yīng)力對于結(jié)構(gòu)剛度的軟化效應(yīng)強(qiáng)于熱變形的硬化效應(yīng);在固有頻率回升階段,熱變形的硬化效用起決定性作用。線性響應(yīng)范圍內(nèi),平板在熱噪聲復(fù)合載荷下的加速度響應(yīng)特性取決于其熱模態(tài)。

    振動(dòng)與波;行波管;熱模態(tài);熱屈曲;熱噪聲復(fù)合載荷

    高超聲速飛行器服役時(shí),處于氣動(dòng)載荷、熱載荷、噪聲載荷和機(jī)械載荷的復(fù)雜環(huán)境中。一方面,由于熱對流和熱傳導(dǎo)作用,熱環(huán)境和溫度場隨時(shí)間變化,改變結(jié)構(gòu)的材料性能,同時(shí)使結(jié)構(gòu)產(chǎn)生熱變形,誘發(fā)熱應(yīng)力,進(jìn)而改變結(jié)構(gòu)的固有特性;另一方面,噪聲載荷本質(zhì)上是一種具有一定空間分布和頻率分布特征的隨機(jī)壓力載荷,飛行器結(jié)構(gòu)通常采用的復(fù)合材料壁板結(jié)構(gòu)對噪聲載荷極為敏感,在某些頻率上易產(chǎn)生足夠大的應(yīng)力,使得結(jié)構(gòu)發(fā)生疲勞失效[1]。

    近年來,國內(nèi)外學(xué)者分別從熱噪聲試驗(yàn)[2]和響應(yīng)預(yù)示方法[1]的建立兩個(gè)方面對薄壁結(jié)構(gòu)在熱噪聲復(fù)合載荷下的響應(yīng)和疲勞壽命問題進(jìn)行了諸多研究。文獻(xiàn)[3–5]均利用商用軟件及其二次平臺進(jìn)行熱噪聲復(fù)合環(huán)境下動(dòng)響應(yīng)預(yù)示方法的開發(fā),研究對象涉及梁、板和復(fù)雜結(jié)構(gòu),考慮真實(shí)服役熱噪聲復(fù)合載荷。沈陽航空航天大學(xué)先后對聲載荷和恒定或帶有溫度梯度的熱載荷[6]復(fù)合作用下的各向同性四邊簡支薄板、加筋板[7]、復(fù)合材料薄板[8]進(jìn)行響應(yīng)計(jì)算,結(jié)果表明,熱載荷通過改變結(jié)構(gòu)剛度特性曲線的形狀影響響應(yīng)的非線性特性,噪聲載荷使結(jié)構(gòu)工作在剛度曲線的不同區(qū)域。吳振強(qiáng)等研究了熱環(huán)境對固有模態(tài)的影響,并對比分析了均勻和非均勻溫度場的影響效果[9]。賀爾銘等基于時(shí)域分析法研究了金屬薄壁結(jié)構(gòu)在熱噪聲復(fù)合載荷下的非線性振動(dòng)響應(yīng)特性,采用四種應(yīng)力模型預(yù)測了薄壁梁的熱噪聲疲勞壽命[10]。

    采用頻率域、波數(shù)域的表征方法既能更好地揭示隨機(jī)載荷特征與結(jié)構(gòu)固有特性的相互作用機(jī)制,同時(shí)相對于時(shí)間域、空間域分析方法,也能夠減小計(jì)算成本,提高計(jì)算速度。因此,目前對隨機(jī)脈動(dòng)壓力作用下加筋板殼結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)問題的求解,普遍采用在波數(shù)-頻率域的求解方法[11–15]。其中,隨機(jī)脈動(dòng)壓力采用功率譜密度矩陣的形式輸入,加筋板殼結(jié)構(gòu)采用波數(shù)或者模態(tài)波數(shù)進(jìn)行表征,對波數(shù)進(jìn)行截?cái)?,通過波數(shù)域積分求得結(jié)構(gòu)的動(dòng)力學(xué)響應(yīng)。Moosrainer和Neubiberg在Sysnoise中,采用了耦合的有限元/邊界元方法,研究了薄壁結(jié)構(gòu)在高強(qiáng)度散射聲場中的響應(yīng)[16],其中聲載荷采用PSD矩陣描述,隨機(jī)響應(yīng)由輸入的PSD矩陣和耦合系統(tǒng)的復(fù)頻響函數(shù)矩陣獲得。

    目前,國內(nèi)外的熱噪聲試驗(yàn)設(shè)備多采用高溫行波管裝置。針對此現(xiàn)狀,本文借助有限元求解方法,研究了熱效應(yīng)對四邊固支平板固有模態(tài)特性的影響機(jī)理,并以此為基礎(chǔ),建立求解行波管中高溫薄壁結(jié)構(gòu)的動(dòng)力學(xué)響應(yīng)的預(yù)示方法,并從模態(tài)特征、頻率變化、幅值變化三個(gè)方面對響應(yīng)結(jié)果進(jìn)行詳細(xì)分析。

    1 結(jié)構(gòu)及載荷的有限元建模

    1.1 考慮熱環(huán)境影響的各向異性平板動(dòng)力學(xué)建模

    綜合考慮熱效應(yīng)對平板動(dòng)力學(xué)特性的影響以及熱噪聲復(fù)合載荷的作用,本文采用有限元方法推導(dǎo)得到正交各項(xiàng)異性平板的動(dòng)力學(xué)響應(yīng)控制方程。設(shè)常溫下有限元方程中的結(jié)構(gòu)剛度矩陣只包含傳統(tǒng)剛度矩陣K;高溫條件引起的熱應(yīng)力導(dǎo)致附加應(yīng)力剛度矩陣Kσ產(chǎn)生;氣動(dòng)加熱條件異常嚴(yán)酷,必須考慮薄板結(jié)構(gòu)的大位移變形因素,由此產(chǎn)生附加的非線性剛度矩陣KNL。忽略阻尼效應(yīng)的情況下,含預(yù)應(yīng)力并考慮幾何非線性變形的動(dòng)力學(xué)響應(yīng)分析的控制方程表示為[17]

    其中M為質(zhì)量矩陣,U為位移向量,P在本文中為聲載荷激勵(lì)向量,K為傳統(tǒng)的剛度矩陣,Kσ為應(yīng)力剛度矩陣,表示為

    其中G為應(yīng)變-位移矩陣,Nx,Ny,Nxy,Nyx為薄膜力,Ai為第i個(gè)單元的面積。KNL為考慮非線性時(shí)的結(jié)構(gòu)剛度矩陣,表示為

    式中BL為結(jié)構(gòu)變形矩陣的線性部分,BNL為非線性部分。

    1.2 一維行波聲場建模

    假設(shè)行波管中的一維行波聲壓載荷表達(dá)式為

    其中ω為圓頻率,為波數(shù),v為聲波傳播速度,A為聲壓幅值。將行波聲壓轉(zhuǎn)化為作用在結(jié)構(gòu)網(wǎng)格節(jié)點(diǎn)上的集中載荷,則兩不同節(jié)點(diǎn)m與n的壓力互功率譜密度表示為

    而某一節(jié)點(diǎn)聲壓的自功率譜密度為

    2 熱環(huán)境對結(jié)構(gòu)模態(tài)特性影響分析

    2.1 數(shù)值模型

    本文以尺寸為380 mm×260 mm,厚度為1.5 mm的矩形平板為研究對象,邊界條件為四周固支,實(shí)際暴露于噪聲載荷的面積為320 mm×200 mm。使用Abaqus建立有限元模型,采用殼單元模擬平板結(jié)構(gòu),網(wǎng)格尺寸為5 mm,共計(jì)3 952個(gè)單元。平板采用C/SiC材料,密度為2 100 kg/m3。

    2.2 均勻溫度載荷下平板結(jié)構(gòu)熱振動(dòng)模態(tài)分析

    由式(1)可知,熱效應(yīng)對結(jié)構(gòu)剛度的影響主要體現(xiàn)在三個(gè)方面:

    (1)材料的彈性模量隨溫度改變而發(fā)生變化,引起結(jié)構(gòu)剛度特性的變化;

    (2)結(jié)構(gòu)在高溫環(huán)境或約束狀態(tài)下,有熱應(yīng)力的產(chǎn)生,并引起剛度特性的改變;

    (3)結(jié)構(gòu)非線性幾何變形對剛度的影響。因此,本文借助有限元軟件Abaqus分別就以上3個(gè)因素研究熱效應(yīng)對結(jié)構(gòu)固有特性的影響。

    只考慮材料特性對結(jié)構(gòu)固有模態(tài)影響時(shí),可以得到圖1所示的固有頻率-溫度曲線,可見隨著溫度升高,平板固有頻率緩慢降低。這一現(xiàn)象可由表1中彈性模量隨溫度升高而減小解釋。其中,第6、第7階固有頻率相差無幾,頻率曲線幾近重合。

    圖1 只考慮材料特性,對結(jié)構(gòu)固有頻率

    由式(2)可知,結(jié)構(gòu)溫度變化導(dǎo)致的熱預(yù)應(yīng)力會產(chǎn)生附加剛度。在本文中,先通過靜力學(xué)分析獲得該熱預(yù)應(yīng)力,再對高溫平板進(jìn)行含預(yù)應(yīng)力的動(dòng)力學(xué)分析。同時(shí)考慮材料特性變化和熱預(yù)應(yīng)力對固有頻率的影響,可得到圖2所示的固有頻率-溫度曲線。

    圖2 考慮材料特性和熱應(yīng)力,結(jié)構(gòu)固有頻率

    對比圖1可知,熱預(yù)應(yīng)力使得平板各階模態(tài)頻率隨溫度增加而減小的速率顯著增加,由此可知,結(jié)構(gòu)受熱引發(fā)的材料特性變化和熱預(yù)應(yīng)力均使結(jié)構(gòu)的固有頻率降低。隨著溫度升高,固有頻率自低階至高階依次減小至零。第1階固有頻率減少到零時(shí)的對應(yīng)溫度即為平板的1階臨界屈曲溫度。

    由式(3)可知,幾何非線性變形影響結(jié)構(gòu)的剛度特性,進(jìn)而影響結(jié)構(gòu)固有模態(tài)。同時(shí)考慮材料特性、熱預(yù)應(yīng)力、非線性幾何變形對結(jié)構(gòu)固有頻率的影響,取T0=300 K作為平板初始溫度,首先通過計(jì)算獲得1階臨界屈曲溫度Tc=376 K,并定義屈曲系數(shù)S=(T-T0)/(Tc-T0)。根據(jù)1、2階屈曲振型和屈曲系數(shù)S模擬平板結(jié)構(gòu)的幾何變形,計(jì)算得到不同屈曲系數(shù)S下平板的熱模態(tài)特性,如圖3所示。

    圖3 考慮材料和幾何非線性時(shí),固有頻率隨溫度變化曲線

    可以看出,對于低階模態(tài),平板固有頻率的整個(gè)變化過程可以分為兩個(gè)階段,下降階段和上升階段;在第1個(gè)階段中,固有頻率隨著溫度的升高而減小,在第2個(gè)階段中,固有頻率隨溫度升高而增大。

    根據(jù)上述分析以及式(2)可知,熱效應(yīng)導(dǎo)致材料彈性模量減小、固支板內(nèi)產(chǎn)生的殘余壓應(yīng)力均會產(chǎn)生負(fù)的附加剛度,進(jìn)而降低結(jié)構(gòu)的固有頻率。同時(shí),平板在受熱過程中逐漸變形,且隨著溫度增加,平衡位置的中面撓度越來越大,如式(3)所示,非線性幾何形變將導(dǎo)致正的附加結(jié)構(gòu)剛度產(chǎn)生,使固有頻率增大。根據(jù)式(1),最終剛度由常溫下結(jié)構(gòu)傳統(tǒng)剛度和以上三個(gè)因素共同決定。

    在下降階段,熱變形并不明顯,此時(shí)熱應(yīng)力及材料特性的變化對結(jié)構(gòu)剛度的影響較大,因此隨著溫度升高,固有頻率逐漸降低;溫度繼續(xù)上升,熱變形產(chǎn)生的剛度硬化逐漸增強(qiáng),直至與熱應(yīng)力和材料特性變化產(chǎn)生的剛度軟化效果相平衡,固有頻率第一次達(dá)到極小值,下文將此極小值點(diǎn)定義為平衡點(diǎn),對應(yīng)溫度定義為平衡點(diǎn)溫度。此后,非線性變形產(chǎn)生的剛度硬化起主要作用,固有頻率曲線回升。

    本文通過在Abaqus軟件中添加初始缺陷的方式引入平板的熱變形,根據(jù)文獻(xiàn)[17]中論證可知,初始缺陷使得結(jié)構(gòu)固有頻率不會隨著溫度升高而減小至零,同時(shí)平衡點(diǎn)溫度低于其臨界屈曲溫度,亦即圖3中固有頻率最低點(diǎn)對應(yīng)的屈曲系數(shù)小于1。

    進(jìn)一步分析圖3可以發(fā)現(xiàn),(1,1)與(2,1)階振型、(3,1)與(2,2)階振型、(1,2)與(2,2)階振型、(3,1)與(3,2)階振型、(1,1)與(1,2)階振型、(2,1)與(2,2)階振型之間先后發(fā)生了模態(tài)順序互換現(xiàn)象。

    屈曲后的薄板模態(tài)還會發(fā)生圖4所示的振型漸變現(xiàn)象。以第1階模態(tài)為例,在S=0時(shí),是(1,1)階模態(tài);當(dāng)S=1.2時(shí),原來的模態(tài)被分割為2個(gè)同相振動(dòng)的模態(tài)(與第2階模態(tài)存在差別),即這階模態(tài)融合了第1階模態(tài)和第2階模態(tài)的特征。在溫度繼續(xù)升高時(shí),第3階模態(tài)(1,2)同樣變成2個(gè)同相的振動(dòng)模態(tài)。

    2.3 熱噪聲復(fù)合載荷下平板的動(dòng)力學(xué)響應(yīng)分析

    對不同溫度下的均勻受熱矩形平板施加一維行波聲載荷,其聲壓幅值不隨頻率變化,如圖5所示。

    圖4 模態(tài)振型漸變現(xiàn)象

    圖5 噪聲載荷功率譜

    當(dāng)W值取14 087.94 Pa2/Hz時(shí),總聲壓級為162 dB。提取薄板中點(diǎn)和受載區(qū)域長四等分點(diǎn)的加速度響應(yīng)分別如圖6、圖7所示。

    結(jié)合圖4所示的相應(yīng)屈曲系數(shù)下平板的固有模態(tài)可以發(fā)現(xiàn),在本節(jié)所研究的溫度范圍內(nèi),平板中點(diǎn)的加速度響應(yīng)峰值均落在第(1,1)、(3,1)階振型對應(yīng)固有頻率處;長四等分點(diǎn)的加速度響應(yīng)峰值均落在(1,1)、(2,1)、(3、1)階振型對應(yīng)的固有頻率處。除了長四等分點(diǎn)在屈曲系數(shù)S等于1.4時(shí)的1、2階固有頻率非常接近(分別為340 Hz和344 Hz),響應(yīng)曲線最大值對應(yīng)的是第2階固有頻率以外,其余響應(yīng)曲線均在(1,1)階振型對應(yīng)頻率處取得最大值。由上述分析可知,平板在熱噪聲復(fù)合載荷下的加速度響應(yīng)主要受其模態(tài)特性的影響。

    進(jìn)一步研究圖6、圖7,可以發(fā)現(xiàn),達(dá)到平衡點(diǎn)溫度之前,平板的動(dòng)力學(xué)響應(yīng)隨著溫度的升高逐漸向低頻方向移動(dòng)。達(dá)到平衡點(diǎn)溫度之后,結(jié)構(gòu)響應(yīng)隨溫度增加向高頻方向移動(dòng),這一變化趨勢主要是由熱環(huán)境下的固有頻率特性決定的。同時(shí)可以發(fā)現(xiàn),達(dá)到平衡點(diǎn)溫度之后,相鄰共振峰之間的間隔隨溫度升高而逐漸減小,共振峰向中間頻率靠近。

    當(dāng)平板熱膨脹受限制時(shí),結(jié)構(gòu)內(nèi)部將有壓力產(chǎn)生。達(dá)到平衡點(diǎn)溫度之前,熱應(yīng)力導(dǎo)致的軟化效應(yīng)起決定性作用,結(jié)構(gòu)剛度減小,位移響應(yīng)幅值增加,同時(shí),溫度升高使得1階共振峰出現(xiàn)的頻率降低。對于穩(wěn)定狀態(tài),加速度響應(yīng)幅值等于位移幅值乘以響應(yīng)頻率的平方。因此,響應(yīng)頻率和位移幅值的變化使得加速度響應(yīng)幅值有所波動(dòng)。達(dá)到平衡點(diǎn)溫度之后,應(yīng)變硬化起決定性作用,結(jié)構(gòu)剛度增大,位移幅值減小,同時(shí),1階共振峰出現(xiàn)的頻率升高,響應(yīng)頻率和位移幅值的變化同樣引起結(jié)構(gòu)加速度響應(yīng)幅值的波動(dòng)。

    圖6 162 dB聲壓下受熱薄板中點(diǎn)加速度響應(yīng)功率譜密度

    圖7 162 dB聲壓下受熱薄板長四等分點(diǎn)加速度響應(yīng)功率譜密度

    3 結(jié)語

    首先,本文基于薄板理論,從材料特性、熱預(yù)應(yīng)力、幾何非線性變形三個(gè)方面研究了熱效應(yīng)對平板結(jié)構(gòu)固有模態(tài)特性的影響。之后,以熱模態(tài)計(jì)算結(jié)果為基礎(chǔ),采用解析方法計(jì)算了熱環(huán)境中四邊固支C/SiC復(fù)合材料平板在一維行波聲場作用下的動(dòng)力學(xué)響應(yīng),并從模態(tài)特征、頻率變化、幅值變化等方面進(jìn)行比較分析。獲得結(jié)論如下:

    (1)結(jié)構(gòu)溫度的變化會引起材料特性的變化、導(dǎo)致殘余應(yīng)力以及幾何變形,這些因素都會較大地影響四端固支平板的剛度特性,因此在研究四端固支平板的動(dòng)力學(xué)特性時(shí),需要綜合考慮這些因素的影響。

    (2)達(dá)到平衡點(diǎn)溫度之前,材料彈性模量的減小和熱應(yīng)力對于結(jié)構(gòu)剛度的軟化效應(yīng)強(qiáng)于非線性幾何變形的硬化效應(yīng),因此,隨著溫度升高,結(jié)構(gòu)固有頻率逐漸降低;達(dá)到平衡點(diǎn)溫度之后,熱變形的硬化效用強(qiáng)于熱應(yīng)力的軟化效用,因此,隨著溫度升高,固有頻率逐漸回升。

    (3)隨著溫度升高,兩階模態(tài)振型之間出現(xiàn)順序互換現(xiàn)象;當(dāng)溫度高于臨界屈曲溫度后,平板的部分模態(tài)出現(xiàn)振型漸變現(xiàn)象。

    (4)平板在熱噪聲復(fù)合載荷下的加速度響應(yīng)主要受其模態(tài)振型的影響。達(dá)到平衡點(diǎn)溫度之前,平板的動(dòng)力學(xué)響應(yīng)隨著溫度的升高逐漸向低頻方向移動(dòng);達(dá)到平衡點(diǎn)溫度之后,結(jié)構(gòu)響應(yīng)峰值隨溫度增加向高頻方向移動(dòng),且共振峰向中間頻率靠近。

    [1]郭靜,吳振強(qiáng),張偉,等.熱噪聲復(fù)合環(huán)境下飛行器結(jié)構(gòu)動(dòng)響應(yīng)預(yù)示技術(shù)研究進(jìn)展[J].強(qiáng)度與環(huán)境,2014,41(6):1-10.

    [2]吳振強(qiáng),任方,張偉,等.飛行器結(jié)構(gòu)熱噪聲試驗(yàn)的研究進(jìn)展[J].導(dǎo)彈與航天運(yùn)載技術(shù),2010(2):24-30.

    [3]GORDON R W,HOLLKAMP J J.Reduced-order models for acoustic response prediction[R].AFRL-RB-WP-TR-2011-3040,2011.

    [4]LIGUORE S L,PITT D M,WHITE E V.Air vehicle technology integration program(AVTIP)Delivery order 0086:applied nonlinear low order response prediction methods evaluation[R].The Boeing Company,2009.

    [5]LIGUORE S L,THOMAS M J,PITT D M.Application and demonstration of nonlinear reduced order modeling(NLROM) for thermal/acoustic response[C].AIAA 2012-1626,53stAIAA/ASME/ASCE/AHS/ASCstructures,structural dynamics and materials conference,2012.

    [6]沙云東,朱林,欒孝馳,等.帶有溫度梯度的熱載荷與聲載荷作用下薄板動(dòng)態(tài)響應(yīng)[J].振動(dòng)與沖擊,2014,33(18):102-109.

    [7]鐘皓杰,沙云東,魏靜.加筋板在熱噪聲載荷作用下的非線性動(dòng)態(tài)響應(yīng)分析[J].沈陽航空航天大學(xué)學(xué)報(bào),2013,30(6):1-7.

    [8]鮑冬冬,沙云東,蔣娜娜.復(fù)合材料薄壁結(jié)構(gòu)在熱噪聲載荷作用下的非線性動(dòng)態(tài)響應(yīng)特性分析[J].沈陽航空航天大學(xué)學(xué)報(bào),2013,30(1):39-42.

    [9]吳振強(qiáng),程昊,張偉,等.熱環(huán)境對飛行器壁板結(jié)構(gòu)動(dòng)特性的影響[J].航空學(xué)報(bào),2013,34(2):334-342.

    [10]賀爾銘,劉峰,胡亞琪,等.熱噪聲載荷下薄壁結(jié)構(gòu)非線性振動(dòng)響應(yīng)分析及疲勞壽命預(yù)測[J].振動(dòng)與沖擊,2013,32(24):135-139.

    [11]A YU GOLUBEV.Experimental estimate of wave spectra of wall pressure fluctuations of the turbulent boundary layer in the subconvective region[J].Acoustical Physics,2012,58(4):396-403.

    [12]MAXIT L.Wavenumber space and physical space responses of a periodically ribbed plate to a point drive:A discrete approach[J].Applied Acoustics,2009,133(1):146-160.

    [13]MAXIT L,DENNIS V.Prediction of low induced sound and vibration of periodically stiffened plate[J].The Journal of the Acoustical Society of America,2013,133(1):146-160.

    [14]BIRGERSSON F,FINNVEDEN S,NILSSON C-M.A spectral super element for modelling of plate vibration.Part1:generaltheory[J].JournalofSound and Vibration,2005(287):297-314.

    [15]BIRGERSSON F,FINNVEDEN S.A spectral super elementformodelling ofplate vibration.Part2:turbulence excitation[J]. Journal of Sound and Vibration,2005,287:315-328.

    [16]MOOSRAINER M.Acoustic excitation of structures[A].InCAD-FEM UsersMeetingInternationaleFEMTechnologietage,Fredrichshafen[C].2000.

    [17]GENG Q,WANG D,LIU Y,et al.Experimental and numericalinvestigations on dynamic and acoustic responses of a thermal post-buckle plate[J].Sci China Tech Sci,2015,58:1414-1424.

    Thermal-acoustic Response Prediction of High-temperature Flat Plates in Traveling-wave Tubes

    TIAN Yan-ni1,WANG Jian-min2,QIN Zhao-hong2,CHENG Hao2,HUA Hong-xing1
    (1.Institute of Vibration,Shock and Noise Research,State Key Laboratory of Mechanical System and Vibration,Shanghai Jiaotong University,Shanghai 200240,China;2.Science and Technology on Reliability and Environment Engineering Laboratory,Beijing Institute of Structure and Environment Engineering,Beijing 100076,China)

    Research on the dynamic and acoustic responses of flat plates in thermal environments is significant to the optimal design of hypersonic aircrafts.For laboratory tests,the high-temperature traveling-wave tube is usually employed to simulate the thermal and acoustic environment to investigate the thermal-acoustic responses of high-temperature structures.In this paper,the thermal-acoustic responses of high-temperature plates in traveling-wave tubes are studied.Based on thinplate theory and finite element method,the influences of high temperature on the inherent characteristics of flat plates are investigated.The dynamic responses of a simply supported heated C/SiC plates under random acoustic excitation are calculated and analyzed.First of all,the influence of thermal effect on the intrinsic modal characteristics of the flat plate is analyzed.Then,the acoustic load of a one-dimensional traveling wave is added on the surface of the plate.Finally,the analytical method is used to solve dynamic response of the plate under the thermal-acoustic combined loading.The results show that in the heating process,the intrinsic frequency of the plate decreases initially to a minimum value and then increases gradually.In the stage of intrinsic frequency decreasing,the softening effect due to elastic modulus reduction and thermal stresses is stronger than the hardening effect due to the thermal deformation.While in the stage of the intrinsic frequency re-increasing,the hardening effect due to the thermal deformation is dominant.In the range of linear response,the acceleration response of the plate to the thermal-acoustic combined loading depends on its thermal modes.

    vibration and wave;traveling-wave tube;thermal mode;thermal buckling;a combination of thermal and acoustic environments

    V414.3+3

    A

    10.3969/j.issn.1006-1355.2017.05.033

    1006-1355(2017)05-0160-06

    2017-03-30

    973資助項(xiàng)目(111402028)

    田艷妮(1991-),女,山東省威海市人,碩士研究生,目前從事熱噪聲載荷下時(shí)變結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)分析方法研究。

    E-mail:tian_yan_ni@163.com

    華宏星,男,博士生導(dǎo)師。

    Email:hhx@sjtu.edu.cn

    猜你喜歡
    振型固有頻率平板
    關(guān)于模態(tài)綜合法的注記
    縱向激勵(lì)下大跨鋼桁拱橋高階振型效應(yīng)分析
    現(xiàn)場測定大型水輪發(fā)電機(jī)組軸系的固有頻率
    屬于你的平板電腦
    出彩的立體聲及豐富的畫面層次 華為|平板M6
    塔腿加過渡段輸電塔動(dòng)力特性分析
    結(jié)構(gòu)振型幾何辨識及應(yīng)用研究
    山西建筑(2015年14期)2015-06-05 09:37:07
    總溫總壓測頭模態(tài)振型變化規(guī)律研究
    A novel functional electrical stimulation-control system for restoring motor function of post-stroke hemiplegic patients
    八寸新標(biāo)桿四核皓麗H8平板發(fā)布
    大话2 男鬼变身卡| 国产成人精品福利久久| 黄色一级大片看看| 亚洲色图 男人天堂 中文字幕| 久久久精品国产亚洲av高清涩受| 在线观看国产h片| 91精品国产国语对白视频| 男男h啪啪无遮挡| 久久精品久久精品一区二区三区| 我的亚洲天堂| 国产 精品1| 欧美av亚洲av综合av国产av | 国产一区二区三区av在线| 宅男免费午夜| 伊人久久国产一区二区| 如何舔出高潮| 无限看片的www在线观看| 99热网站在线观看| 纵有疾风起免费观看全集完整版| 久久久久久久久久久免费av| 欧美激情高清一区二区三区 | 九九爱精品视频在线观看| 国产精品久久久人人做人人爽| 亚洲欧美一区二区三区黑人| 中文字幕人妻丝袜制服| 岛国毛片在线播放| 亚洲美女视频黄频| 只有这里有精品99| 国产成人精品久久二区二区91 | 天天添夜夜摸| 亚洲一码二码三码区别大吗| 交换朋友夫妻互换小说| 国产一区亚洲一区在线观看| 国产老妇伦熟女老妇高清| 亚洲国产精品国产精品| 99精国产麻豆久久婷婷| 成人国语在线视频| 大片免费播放器 马上看| 女人被躁到高潮嗷嗷叫费观| 伦理电影大哥的女人| 国产 一区精品| 丰满迷人的少妇在线观看| 侵犯人妻中文字幕一二三四区| 日韩 欧美 亚洲 中文字幕| 欧美激情高清一区二区三区 | 国产又色又爽无遮挡免| 青草久久国产| 国精品久久久久久国模美| 婷婷色av中文字幕| 亚洲在久久综合| 国产黄频视频在线观看| av在线播放精品| 亚洲精品中文字幕在线视频| 久久久久久久久久久免费av| 一区二区av电影网| 一级毛片我不卡| 免费黄网站久久成人精品| 国产欧美亚洲国产| 国产成人91sexporn| 美女福利国产在线| 国产乱人偷精品视频| 成人亚洲精品一区在线观看| 最近最新中文字幕免费大全7| 9191精品国产免费久久| 亚洲精品国产av蜜桃| 深夜精品福利| 99精品久久久久人妻精品| 久久人人爽av亚洲精品天堂| 日日撸夜夜添| 丰满乱子伦码专区| 亚洲三区欧美一区| 老汉色av国产亚洲站长工具| 日本wwww免费看| 在线观看免费高清a一片| 国产欧美日韩一区二区三区在线| 成人亚洲欧美一区二区av| 极品人妻少妇av视频| 好男人视频免费观看在线| 嫩草影视91久久| 亚洲av日韩精品久久久久久密 | 亚洲综合色网址| 不卡av一区二区三区| 国产成人系列免费观看| 日韩av在线免费看完整版不卡| svipshipincom国产片| 国产黄色视频一区二区在线观看| 午夜免费男女啪啪视频观看| 午夜av观看不卡| 校园人妻丝袜中文字幕| 亚洲图色成人| 久久精品国产亚洲av高清一级| 精品午夜福利在线看| 老司机影院成人| 美女中出高潮动态图| 最近最新中文字幕大全免费视频 | 亚洲一区中文字幕在线| 亚洲精品一二三| 亚洲精品久久久久久婷婷小说| 中国国产av一级| 国产无遮挡羞羞视频在线观看| 91精品三级在线观看| 美女主播在线视频| 我要看黄色一级片免费的| 日本91视频免费播放| 欧美最新免费一区二区三区| 可以免费在线观看a视频的电影网站 | 美女国产高潮福利片在线看| 精品人妻在线不人妻| 欧美亚洲 丝袜 人妻 在线| 制服人妻中文乱码| 国产成人欧美| 国产男女超爽视频在线观看| 亚洲精品国产一区二区精华液| 在线免费观看不下载黄p国产| 亚洲第一青青草原| bbb黄色大片| 在线观看免费高清a一片| 精品亚洲成国产av| 久久99精品国语久久久| 高清不卡的av网站| 久久国产精品男人的天堂亚洲| 国产女主播在线喷水免费视频网站| 久久久久精品国产欧美久久久 | 久久人人爽av亚洲精品天堂| 999久久久国产精品视频| 午夜福利在线免费观看网站| 色吧在线观看| 老司机在亚洲福利影院| 精品人妻熟女毛片av久久网站| 亚洲熟女精品中文字幕| 亚洲四区av| 亚洲精品久久成人aⅴ小说| 最新的欧美精品一区二区| 丰满迷人的少妇在线观看| 18禁动态无遮挡网站| 黑丝袜美女国产一区| 久久久国产一区二区| 国产乱来视频区| 国产一区二区三区综合在线观看| 爱豆传媒免费全集在线观看| 久热爱精品视频在线9| 人人妻人人爽人人添夜夜欢视频| 女人精品久久久久毛片| a 毛片基地| 18禁国产床啪视频网站| 精品一区二区免费观看| 亚洲精品第二区| 在线精品无人区一区二区三| 十分钟在线观看高清视频www| 欧美日韩国产mv在线观看视频| 女人被躁到高潮嗷嗷叫费观| 国产av一区二区精品久久| 赤兔流量卡办理| 久久久精品免费免费高清| 捣出白浆h1v1| 色婷婷久久久亚洲欧美| 日本欧美国产在线视频| 欧美成人精品欧美一级黄| 熟妇人妻不卡中文字幕| 老汉色∧v一级毛片| 天天躁夜夜躁狠狠躁躁| 国产精品二区激情视频| 亚洲欧洲国产日韩| 九色亚洲精品在线播放| 99re6热这里在线精品视频| 日韩一卡2卡3卡4卡2021年| 老司机亚洲免费影院| 另类亚洲欧美激情| 久久性视频一级片| 中文字幕精品免费在线观看视频| 80岁老熟妇乱子伦牲交| 日韩不卡一区二区三区视频在线| 九色亚洲精品在线播放| 99国产综合亚洲精品| 欧美日韩亚洲国产一区二区在线观看 | 一本—道久久a久久精品蜜桃钙片| 国产 精品1| 久久鲁丝午夜福利片| 国产欧美日韩综合在线一区二区| 久久久久精品性色| 亚洲av日韩精品久久久久久密 | 欧美日韩视频高清一区二区三区二| 人人妻,人人澡人人爽秒播 | 精品一区二区三区四区五区乱码 | 久久久久久久大尺度免费视频| 亚洲成色77777| 午夜激情久久久久久久| 久久久久久人人人人人| 国产精品二区激情视频| 自拍欧美九色日韩亚洲蝌蚪91| 免费黄色在线免费观看| 免费观看av网站的网址| 精品少妇久久久久久888优播| av线在线观看网站| 成人三级做爰电影| 侵犯人妻中文字幕一二三四区| 看十八女毛片水多多多| 国产视频首页在线观看| 午夜福利视频在线观看免费| 汤姆久久久久久久影院中文字幕| 午夜福利影视在线免费观看| 岛国毛片在线播放| 久久久国产一区二区| 晚上一个人看的免费电影| 成人亚洲欧美一区二区av| 国产片内射在线| 亚洲精品中文字幕在线视频| 男女边摸边吃奶| 黄色毛片三级朝国网站| 亚洲国产精品999| 精品一区二区免费观看| 亚洲国产欧美日韩在线播放| 日韩人妻精品一区2区三区| 亚洲,一卡二卡三卡| 涩涩av久久男人的天堂| 亚洲精品aⅴ在线观看| 精品少妇久久久久久888优播| 国产又爽黄色视频| 久久久久精品久久久久真实原创| 亚洲五月色婷婷综合| 成人亚洲精品一区在线观看| www.熟女人妻精品国产| 丁香六月天网| 99精品久久久久人妻精品| 曰老女人黄片| 国产av精品麻豆| 欧美精品av麻豆av| 免费少妇av软件| av福利片在线| 青春草国产在线视频| 美女脱内裤让男人舔精品视频| 国产亚洲最大av| 亚洲图色成人| 婷婷色综合大香蕉| 精品少妇一区二区三区视频日本电影 | 国产精品无大码| 老司机影院成人| 捣出白浆h1v1| 9色porny在线观看| 亚洲国产看品久久| 日日撸夜夜添| 嫩草影视91久久| 亚洲精品久久久久久婷婷小说| 精品久久久精品久久久| 久久久久久久久久久免费av| 精品少妇久久久久久888优播| 国产淫语在线视频| 日本wwww免费看| 国产亚洲一区二区精品| av网站免费在线观看视频| 国产成人免费无遮挡视频| 午夜免费观看性视频| 少妇人妻 视频| 丝袜脚勾引网站| 色94色欧美一区二区| 欧美精品av麻豆av| 欧美黑人精品巨大| 啦啦啦在线观看免费高清www| 如何舔出高潮| 两性夫妻黄色片| a级毛片黄视频| 国产在线免费精品| 黄色怎么调成土黄色| 国产成人系列免费观看| 亚洲欧美成人综合另类久久久| 最近的中文字幕免费完整| 在线天堂最新版资源| 9色porny在线观看| 1024香蕉在线观看| 免费久久久久久久精品成人欧美视频| 国产精品久久久av美女十八| 亚洲美女视频黄频| 天天躁夜夜躁狠狠躁躁| 别揉我奶头~嗯~啊~动态视频 | 一本久久精品| 久久精品国产亚洲av涩爱| 色综合欧美亚洲国产小说| 国产欧美日韩一区二区三区在线| 男人操女人黄网站| 久久人人爽人人片av| 国产高清不卡午夜福利| 日日摸夜夜添夜夜爱| 1024香蕉在线观看| 国产视频首页在线观看| 日本欧美视频一区| 可以免费在线观看a视频的电影网站 | 亚洲精品日韩在线中文字幕| 欧美精品亚洲一区二区| 男的添女的下面高潮视频| 亚洲国产欧美网| 一本一本久久a久久精品综合妖精| 国产男人的电影天堂91| 一级爰片在线观看| 成人国产av品久久久| 丝瓜视频免费看黄片| 国产在线免费精品| 极品人妻少妇av视频| 日日摸夜夜添夜夜爱| 丝袜喷水一区| 18在线观看网站| 亚洲国产av新网站| 看免费av毛片| 老熟女久久久| 一区福利在线观看| 国产免费现黄频在线看| 久久99精品国语久久久| 亚洲国产中文字幕在线视频| 9热在线视频观看99| 欧美 日韩 精品 国产| 精品视频人人做人人爽| 午夜福利网站1000一区二区三区| 亚洲精品视频女| 天堂8中文在线网| 成年动漫av网址| 18禁裸乳无遮挡动漫免费视频| 午夜精品国产一区二区电影| 亚洲天堂av无毛| 一边摸一边抽搐一进一出视频| 天天添夜夜摸| 午夜91福利影院| 九草在线视频观看| 18禁观看日本| 99精品久久久久人妻精品| 国产亚洲一区二区精品| 91国产中文字幕| 曰老女人黄片| 美女大奶头黄色视频| 国产精品嫩草影院av在线观看| 国产男女超爽视频在线观看| 亚洲国产最新在线播放| 亚洲av电影在线进入| 久久99热这里只频精品6学生| 两个人免费观看高清视频| 亚洲综合色网址| 国产精品麻豆人妻色哟哟久久| 亚洲美女黄色视频免费看| 一二三四在线观看免费中文在| 人成视频在线观看免费观看| 777久久人妻少妇嫩草av网站| 男女免费视频国产| 亚洲国产精品成人久久小说| 免费观看av网站的网址| 国产精品熟女久久久久浪| 一级毛片黄色毛片免费观看视频| 欧美乱码精品一区二区三区| 亚洲成人av在线免费| 日日爽夜夜爽网站| 久久精品人人爽人人爽视色| 热re99久久国产66热| 大香蕉久久网| 亚洲国产欧美网| 色94色欧美一区二区| 国产免费视频播放在线视频| 久久久久久久久久久免费av| 亚洲国产欧美网| 久久久久久久久久久免费av| 久久精品国产综合久久久| 亚洲av电影在线进入| 91精品伊人久久大香线蕉| 欧美人与善性xxx| 国产精品一区二区在线不卡| 欧美老熟妇乱子伦牲交| 精品一区二区三区av网在线观看 | 日本vs欧美在线观看视频| 日韩制服丝袜自拍偷拍| 国产精品秋霞免费鲁丝片| www.自偷自拍.com| 久久这里只有精品19| 中文字幕亚洲精品专区| 久久99一区二区三区| 久久 成人 亚洲| 丰满乱子伦码专区| 久久久久国产精品人妻一区二区| 久久精品人人爽人人爽视色| 天堂8中文在线网| 精品人妻一区二区三区麻豆| 哪个播放器可以免费观看大片| 久久精品国产综合久久久| 午夜福利一区二区在线看| 成人国产av品久久久| 日本91视频免费播放| 操出白浆在线播放| 中文乱码字字幕精品一区二区三区| 免费看不卡的av| √禁漫天堂资源中文www| 涩涩av久久男人的天堂| 天天添夜夜摸| 久久青草综合色| 久久久精品免费免费高清| 18禁裸乳无遮挡动漫免费视频| 国产爽快片一区二区三区| 日韩av不卡免费在线播放| 黄色毛片三级朝国网站| 午夜福利乱码中文字幕| 成人国产麻豆网| 国产在视频线精品| 国产精品国产av在线观看| 国产精品秋霞免费鲁丝片| 9热在线视频观看99| 一本一本久久a久久精品综合妖精| 国产精品亚洲av一区麻豆 | 国产福利在线免费观看视频| 我的亚洲天堂| 精品一区二区三卡| 免费在线观看黄色视频的| 国产精品99久久99久久久不卡 | 9191精品国产免费久久| 国产探花极品一区二区| 丰满乱子伦码专区| 午夜福利视频在线观看免费| 操出白浆在线播放| 天堂俺去俺来也www色官网| 男的添女的下面高潮视频| 亚洲国产最新在线播放| 高清视频免费观看一区二区| 最近中文字幕高清免费大全6| 国产成人精品福利久久| 毛片一级片免费看久久久久| 免费女性裸体啪啪无遮挡网站| 高清欧美精品videossex| 成人毛片60女人毛片免费| a级毛片黄视频| 亚洲精品日韩在线中文字幕| 亚洲国产欧美日韩在线播放| 99久久精品国产亚洲精品| 精品亚洲成国产av| 国产亚洲午夜精品一区二区久久| 人人妻人人爽人人添夜夜欢视频| 91精品伊人久久大香线蕉| 叶爱在线成人免费视频播放| 欧美 亚洲 国产 日韩一| 男女免费视频国产| 国产精品国产三级国产专区5o| 天天躁狠狠躁夜夜躁狠狠躁| 又大又黄又爽视频免费| 婷婷色av中文字幕| 丝袜人妻中文字幕| 亚洲 欧美一区二区三区| 国产高清国产精品国产三级| 九九爱精品视频在线观看| 欧美精品人与动牲交sv欧美| 99久久综合免费| 国产亚洲精品第一综合不卡| 两性夫妻黄色片| 国产精品久久久久久精品电影小说| 老司机影院毛片| 美女扒开内裤让男人捅视频| 中文字幕制服av| 欧美亚洲 丝袜 人妻 在线| 9191精品国产免费久久| 久久亚洲国产成人精品v| 国产成人欧美在线观看 | 国产精品一区二区在线不卡| 免费少妇av软件| 少妇人妻 视频| 97精品久久久久久久久久精品| 国产成人精品久久二区二区91 | 青春草视频在线免费观看| xxx大片免费视频| 亚洲成人免费av在线播放| 精品国产一区二区久久| 青春草视频在线免费观看| 成人影院久久| 国产欧美亚洲国产| 欧美 亚洲 国产 日韩一| 久久天躁狠狠躁夜夜2o2o | 看免费成人av毛片| 日韩欧美精品免费久久| 满18在线观看网站| 啦啦啦啦在线视频资源| 在线观看免费午夜福利视频| 免费少妇av软件| 国产精品一区二区在线不卡| 久久久国产精品麻豆| av不卡在线播放| 又大又爽又粗| 三上悠亚av全集在线观看| 七月丁香在线播放| 成年动漫av网址| 色视频在线一区二区三区| 国产高清不卡午夜福利| 最近2019中文字幕mv第一页| 日韩伦理黄色片| 在线观看一区二区三区激情| 精品第一国产精品| 久久精品人人爽人人爽视色| 曰老女人黄片| 别揉我奶头~嗯~啊~动态视频 | 三上悠亚av全集在线观看| 精品久久久久久电影网| 国产又爽黄色视频| 观看美女的网站| 精品国产一区二区三区久久久樱花| 午夜激情av网站| av免费观看日本| 涩涩av久久男人的天堂| 国产人伦9x9x在线观看| 欧美最新免费一区二区三区| 亚洲欧美精品综合一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| av在线app专区| 男女国产视频网站| 成人影院久久| 久久久国产精品麻豆| 久久久久精品性色| 欧美日韩一级在线毛片| 亚洲成人免费av在线播放| 欧美精品av麻豆av| 亚洲人成77777在线视频| 水蜜桃什么品种好| 日本一区二区免费在线视频| 女人高潮潮喷娇喘18禁视频| 亚洲国产中文字幕在线视频| 美女高潮到喷水免费观看| 久久久久网色| 最黄视频免费看| 精品国产一区二区三区四区第35| 黑人欧美特级aaaaaa片| 九草在线视频观看| 曰老女人黄片| 国产极品粉嫩免费观看在线| 三上悠亚av全集在线观看| 一级片'在线观看视频| 不卡视频在线观看欧美| 午夜福利视频精品| 天天躁夜夜躁狠狠久久av| 狠狠精品人妻久久久久久综合| e午夜精品久久久久久久| 久久人人97超碰香蕉20202| 又大又黄又爽视频免费| 伊人久久大香线蕉亚洲五| 亚洲伊人久久精品综合| 十八禁高潮呻吟视频| 国产有黄有色有爽视频| 麻豆av在线久日| 国产亚洲一区二区精品| 国产日韩欧美在线精品| 亚洲国产av新网站| 涩涩av久久男人的天堂| 一边亲一边摸免费视频| 亚洲欧美成人精品一区二区| 在线观看国产h片| 久久久久精品人妻al黑| 99久久99久久久精品蜜桃| 免费观看av网站的网址| a 毛片基地| 国产精品成人在线| 久久久精品区二区三区| 亚洲精品视频女| 欧美亚洲日本最大视频资源| 国产麻豆69| 亚洲第一av免费看| 国产亚洲欧美精品永久| 久久亚洲国产成人精品v| 久久 成人 亚洲| 免费女性裸体啪啪无遮挡网站| 丝袜脚勾引网站| 国产精品偷伦视频观看了| 在线亚洲精品国产二区图片欧美| 日日撸夜夜添| 国产精品香港三级国产av潘金莲 | 国产成人精品无人区| av卡一久久| 亚洲av电影在线进入| 欧美人与性动交α欧美软件| 精品福利永久在线观看| 亚洲熟女毛片儿| 青草久久国产| 免费黄网站久久成人精品| 国产成人免费观看mmmm| 久久精品熟女亚洲av麻豆精品| www.熟女人妻精品国产| 亚洲成av片中文字幕在线观看| 中文字幕高清在线视频| 最近中文字幕2019免费版| 色婷婷久久久亚洲欧美| 亚洲av国产av综合av卡| 女人被躁到高潮嗷嗷叫费观| 免费不卡黄色视频| 国产免费又黄又爽又色| 亚洲国产中文字幕在线视频| 熟女av电影| 国产免费视频播放在线视频| 国产有黄有色有爽视频| 国产精品久久久久久精品电影小说| 1024视频免费在线观看| 国产成人欧美| 少妇人妻久久综合中文| 欧美av亚洲av综合av国产av | 天天躁狠狠躁夜夜躁狠狠躁| 一边亲一边摸免费视频| 黄片小视频在线播放| 国精品久久久久久国模美| 精品国产一区二区三区久久久樱花| 中文乱码字字幕精品一区二区三区| 最新的欧美精品一区二区| 纯流量卡能插随身wifi吗| 日韩制服骚丝袜av| 美女中出高潮动态图| 青春草亚洲视频在线观看| 成年人免费黄色播放视频| 精品少妇内射三级| 丝袜美腿诱惑在线| 美女国产高潮福利片在线看| 汤姆久久久久久久影院中文字幕| 97精品久久久久久久久久精品| 69精品国产乱码久久久| 国产亚洲精品第一综合不卡| www.精华液| 国产在线视频一区二区| 免费在线观看视频国产中文字幕亚洲 | 亚洲欧洲国产日韩| 在线观看一区二区三区激情| 中文字幕最新亚洲高清| 国产精品.久久久|