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

    雙金屬復(fù)合翅片管振動特性的研究

    2016-04-21 01:36:46郭建英馬騰飛劉生寶白艷艷太原理工大學(xué)礦業(yè)工程學(xué)院太原理工大學(xué)應(yīng)用力學(xué)與生物工程研究所太原000中石化寧波工程有限公司寧波50太原理工大學(xué)機械工程學(xué)院太原000
    振動與沖擊 2016年5期

    郭建英, 馬騰飛. 劉生寶, 白艷艷(.太原理工大學(xué) 礦業(yè)工程學(xué)院;. 太原理工大學(xué) 應(yīng)用力學(xué)與生物工程研究所,太原 000;.中石化寧波工程有限公司,寧波 50; .太原理工大學(xué) 機械工程學(xué)院,太原 000)

    ?

    雙金屬復(fù)合翅片管振動特性的研究

    郭建英1, 馬騰飛2. 劉生寶3, 白艷艷4(1.太原理工大學(xué) 礦業(yè)工程學(xué)院;2. 太原理工大學(xué) 應(yīng)用力學(xué)與生物工程研究所,太原030024;3.中石化寧波工程有限公司,寧波315103; 4.太原理工大學(xué) 機械工程學(xué)院,太原030024)

    摘要:雙金屬復(fù)合翅片管是一種高效傳熱元件,容易發(fā)生流體誘導(dǎo)振動破壞,對其進行振動模態(tài)理論分析具有重要的工程指導(dǎo)意義。針對雙金屬復(fù)合翅片管的結(jié)構(gòu)特征,將其簡化為串、并聯(lián)剛度系統(tǒng),采用組合截面等效彎曲剛度、等效扭轉(zhuǎn)剛度和等效抗拉壓剛度,并結(jié)合等效質(zhì)量和等效轉(zhuǎn)動慣量的方法,對其彎曲、扭轉(zhuǎn)和軸向振動模態(tài)進行理論解析。為了驗證理論分析方法的準確性,對雙金屬復(fù)合翅片管的振動模態(tài)進行了實驗測試和有限元分析。研究了翅片幾何參數(shù)對雙金屬復(fù)合翅片管振動頻率的影響規(guī)律。結(jié)果表明,對于矩形翅片形式的鋼鋁雙金屬翅片管,其彎曲、扭轉(zhuǎn)和軸向振動頻率均隨翅片高度和翅片厚度的增大而減小,隨翅片間距的增大而增大。

    關(guān)鍵詞:雙金屬翅片管;振動模態(tài);等效截面剛度;等效質(zhì)量;等效轉(zhuǎn)動慣量

    換熱器是國民經(jīng)濟和工業(yè)生產(chǎn)領(lǐng)域中應(yīng)用十分廣泛的熱量交換設(shè)備,約占設(shè)備總量的30%左右。為了提高換熱器的傳熱效率,常用各種強化換熱元件(如翅片管)代替?zhèn)鹘y(tǒng)光滑管[1]。翅片管有單金屬翅片管和雙金屬復(fù)合翅片管。目前,雙金屬整體螺旋復(fù)合翅片換熱管在國內(nèi)外的應(yīng)用非常廣泛,已遍及能源、動力及化工加工的各種領(lǐng)域[2-5]。雙金屬復(fù)合翅片管由基管和包裹在外的整體螺旋翅片管組成,基管的材質(zhì)可根據(jù)熱流體腐蝕情況和工作壓力選定,通常為碳鋼、不銹鋼或銅;外管可選用延展性好,抗大氣腐蝕能力強,傳熱性能優(yōu)良的金屬,如鋁或銅;將鋁管或銅管套在基管外面,然后用專用軋機將鋁或銅管擠壓成翅片,同時在基管外形成一層薄的底層將基管緊緊裹住。

    換熱器在使用過程中,極容易產(chǎn)生流體誘導(dǎo)振動[6-7],它會導(dǎo)致劇烈的噪聲及換熱管的斷裂損壞。有效預(yù)防和控制換熱器流體誘導(dǎo)振動破壞的首要措施是,將換熱管的固有振動頻率和流體干擾力頻率分隔開。因此,換熱器的設(shè)計必須進行振動計算,即準確計算換熱管的固有振動頻率[8]。雙金屬整體螺旋翅片管在結(jié)構(gòu)上屬于不同材料的雙層殼體復(fù)合結(jié)構(gòu),且外管由翅片組成,其截面剛度與質(zhì)量沿軸向分布不均勻。目前,雖然采用有限元法可以獲得雙金屬復(fù)合翅片管較精確的動力學(xué)特性,但不便于工程設(shè)計采用。建立準確的雙金屬復(fù)合翅片管振動頻率計算模型,對雙金屬翅片管換熱器的設(shè)計具有重要的理論指導(dǎo)和工程應(yīng)用價值。

    對于雙層殼體復(fù)合結(jié)構(gòu),謝志勇等[9-11]的研究發(fā)現(xiàn),內(nèi)外殼的耦合作用會影響其動態(tài)特性。而對于翅片管,工程上常采用經(jīng)驗公式將其簡化為有效直徑的光滑管,但計算結(jié)果誤差較大。

    本文擬將雙金屬翅片管看作串、并聯(lián)剛度組合系統(tǒng),采用組合截面等效剛度法并結(jié)合等效質(zhì)量和等效轉(zhuǎn)動慣量,推導(dǎo)雙金屬復(fù)合翅片管彎曲、扭轉(zhuǎn)及軸向振動模態(tài)特性的理論解析模型;并采用實驗測試和有限元模擬的方法對理論結(jié)果進行驗證;同時分析不同翅片幾何尺寸時雙金屬復(fù)合翅片管振動頻率的變化規(guī)律。

    1雙金屬復(fù)合翅片管幾何結(jié)構(gòu)

    雙金屬復(fù)合翅片管為兩種不同材料的組合管,內(nèi)外管之間牢固連接而無相對錯動;外管翅片形狀為矩形,且翅片高而薄,其結(jié)構(gòu)如圖1所示。圖1中,管子總長為L,內(nèi)管內(nèi)徑為D1,外徑為D2,內(nèi)管材料的楊氏彈性模量為E1,剪切彈性模量為G1,密度為ρ1。外翅片管內(nèi)徑為D2,翅根直徑為D3,翅頂直徑為D4,翅片高度為h,翅片厚度為a,翅片間距為b,外翅管材料的楊氏彈性模量為E2,剪切彈性模量為G2,密度為ρ2。

    圖1 雙金屬復(fù)合翅片管幾何結(jié)構(gòu)與參數(shù)Fig.1 The geometry and parameters of bimetallic finned tube

    2雙金屬復(fù)合翅片管彎曲振動模態(tài)

    2.1理論解析

    在橫向載荷激勵下,雙金屬復(fù)合翅片管的振動模態(tài)為彎曲振動。假設(shè)雙金屬翅片管為線彈性體,材料均勻連續(xù),各向同性,且振動過程中各截面仍保持為平面。由于雙金屬翅片管的長度遠大于其直徑,為細長結(jié)構(gòu),因此采用Bernoulli-Euler梁理論建立其橫向彎曲自由振動動力學(xué)方程:

    (1)

    式中,y(x,t)為距原點x處管子截面在t時刻的橫向位移。

    由圖1可知,雙金屬翅片管的剛度和質(zhì)量沿軸向(x方向)分布不均勻。為了便于解析,本文將雙金屬翅片管等效為等截面管,并將EeIe定義為雙金屬翅片管的等效截面彎曲剛度,me定義為單位長度的等效質(zhì)量。求解該微分方程,得雙金屬翅片管彎曲振動的各階固有頻率為:

    (2)

    式中,λwi(i=1, 2, 3,…)為各階彎曲振動頻率常數(shù),取決于邊界支撐條件。

    由式(2)可知,準確獲得雙金屬復(fù)合翅片管的等效截面彎曲剛度EeIe及其等效質(zhì)量me是精確解析其彎曲振動頻率的關(guān)鍵。為此,基于質(zhì)量等效原則,本文將圖1所示的雙金屬復(fù)合翅片管的螺旋結(jié)構(gòu)簡化為槽狀結(jié)構(gòu),再根據(jù)其翅片周期性分布的結(jié)構(gòu)特征,選取其中一個翅片單元,并利用組合截面法來推導(dǎo)其等效截面彎曲剛度和等效質(zhì)量。圖2(a)為簡化后的雙金屬復(fù)合翅片管槽狀結(jié)構(gòu),圖2(b) 為一個翅片單元的幾何參數(shù)示意圖。

    圖2 雙金屬翅片管槽狀結(jié)構(gòu)和一個翅片單元的幾何參數(shù)圖Fig.2 Simplified trough-like geometry of bimetallic finned tube and the geometry diagram of a fin unit

    從圖2(a)可以看出,雙金屬復(fù)合翅片管相當于內(nèi)管與外翅管的并聯(lián)剛度系統(tǒng),而外翅管可看做不同截面管段的串聯(lián)剛度系統(tǒng)。對圖2(b)所示的一個翅片單元管段(長度為a+b)而言:內(nèi)管(內(nèi)徑D1,外徑D2)的截面剛度沿軸向(x方向)均勻分布,定義其橫截面積為A1,其截面對中性軸的慣性矩為I1,其截面的極慣性矩為Ip1,該部分管段的轉(zhuǎn)動慣量為J1。相應(yīng)長度的外翅管相當于兩段階梯截面的串聯(lián)剛度系統(tǒng),其中含翅片部分(內(nèi)徑D2,外徑D4),定義其橫截面積為A21,該截面對中性軸的慣性矩為I21,截面的極慣性矩為Ip21,該部分轉(zhuǎn)動慣量為J21;不含翅片部分(內(nèi)徑D2,外徑D3),定義其橫截面積為A22,該截面對中性軸的慣性矩為I22,截面的極慣性矩為Ip22,該部分轉(zhuǎn)動慣量為J22。這些參數(shù)的值分別為:

    (3)

    (4)

    (5)

    J1=(a+b)ρ1Ip1,J21=aρ2Ip21,J22=bρ2Ip22

    (6)

    對于這樣的串、并聯(lián)剛度系統(tǒng)的組合截面,可以推導(dǎo)出其等效截面彎曲剛度為:

    (7)

    另外,可以推導(dǎo)出雙金屬復(fù)合翅片管單位長度的等效質(zhì)量為:

    (8)

    這樣,將式(7)、(8)代入式(2),可以解出雙金屬復(fù)合翅片管的各階彎曲振動固有頻率。

    2.2實驗測試

    為了驗證等效剛度法解析雙金屬復(fù)合翅片管彎曲振動模態(tài)方法的精確性,本文采用北京東方振動和噪聲技術(shù)研究所的INV3018A智能信號采集系統(tǒng)和DASP-V10振動信號分析系統(tǒng),并用捶擊法對鋼鋁雙金屬復(fù)合翅片管在兩端自由狀態(tài)下的橫向彎曲振動固有頻率進行了實驗測試。兩端自由支撐條件是通過采用兩條高彈性皮筋懸掛翅片管的方法實現(xiàn)的。實驗測試所用的雙金屬復(fù)合翅片管的內(nèi)基管為碳鋼管,外管為鋁翅片管,結(jié)構(gòu)幾何參數(shù)見表1。實驗測試獲得了兩端自由狀態(tài)下,鋼鋁雙金屬復(fù)合翅片管振動模態(tài)的自譜-FFT分析的幅值譜(PEAK)圖,如圖3所示。

    表1 鋼鋁雙金屬復(fù)合翅片管結(jié)構(gòu)

    圖3 實測幅頻特性曲線Fig.3 Measured amplitude-frequency characteristic curve

    2.2有限元分析

    本文采用有限元軟件ANSYS對鋼鋁雙金屬復(fù)合翅片管進行了模態(tài)分析。實驗測試用的翅片管為整體螺旋結(jié)構(gòu),在進行理論解析時將其簡化為槽狀翅片管結(jié)構(gòu)。在進行有限元分析時,分別建立了兩種模型,一種為槽狀翅片管模型,見圖4(a);另一種為螺旋翅片管模型,見圖4(b)。有限元模型的幾何尺寸與實驗用管完全相同,單元類型采用四面體10節(jié)點三維實體單元,內(nèi)外管之間無相對錯動的牢固連接通過GLUE操作來實現(xiàn)。

    圖4 槽狀翅片管模型及網(wǎng)格和螺旋翅片管模型及網(wǎng)格 Fig.4 Trough-like and Helix model and mesh of bimetallic finned tube

    2.3彎曲振動頻率結(jié)果分析

    表2顯示了兩端自由條件下,鋼鋁雙金屬復(fù)合翅片管前四階彎曲振動固有頻率的實測值、有限元模擬值及理論解析值的比較。在表2中,誤差1指的是理論解析值與螺旋狀模型有限元分析值之間的差別,誤差2指的是理論解析值與實測值之間的差別。表2數(shù)據(jù)表明,① 對雙金屬復(fù)合翅片管振動模態(tài)進行有限元分析時,槽狀模型和螺旋狀模型的分析結(jié)果非常接近,說明將雙金屬整體螺旋翅片管簡化為槽狀模型是可行的。② 采用等效截面彎曲剛度EeIe并結(jié)合等效質(zhì)量me,計算得到的雙金屬復(fù)合翅片管的前四階彎曲振動固有頻率,與實測值相比,一階頻率誤差為5.00%;四階頻率誤差為6.96%。通過分析發(fā)現(xiàn)該誤差產(chǎn)生的主要原因是,實測用的雙金屬翅片管的幾何尺寸加工不均勻。如實測的雙金屬翅片管的單位長度質(zhì)量me為1.386 kg/m,而按照實測的管子幾何參數(shù)計算得到的單位長度質(zhì)量me為1.302 kg/m,二者誤差達到6.06%。這種加工不均勻同樣會影響管子的彎曲剛度,最終使振動頻率的解析值與實測值之間產(chǎn)生大的誤差。③雙金屬翅片管前四階彎曲振動固有頻率的解析值與有限元模擬值相比,一階頻率誤差最小,僅為1.39%;四階頻率誤差最大,為3.63%;這表明采用等效剛度法并結(jié)合等效質(zhì)量解析雙金屬復(fù)合翅片管的彎曲振動模態(tài)是精確的。

    表2 雙金屬復(fù)合翅片管彎曲振動固有頻率

    3雙金屬復(fù)合翅片管扭轉(zhuǎn)與軸向振動模態(tài)

    3.1扭轉(zhuǎn)振動模態(tài)解析

    在扭矩激勵下,雙金屬復(fù)合翅片管表現(xiàn)為扭轉(zhuǎn)振動模態(tài)。扭轉(zhuǎn)振動為高階模態(tài),必須考慮轉(zhuǎn)動慣量的影響,因此其扭轉(zhuǎn)自由振動的動力學(xué)方程為:

    (9)

    式中,θ(x,t)為在距原點x處截面在t時刻的角位移。與彎曲振動解析方法相同,這里將GeIpe定義為雙金屬翅片管的等效截面扭轉(zhuǎn)剛度,Je定義為其單位長度的等效轉(zhuǎn)動慣量。求解得雙金屬翅片管扭轉(zhuǎn)振動的各階固有頻率為:

    (10)

    式中,λni(i=1, 2, 3…)為各階扭轉(zhuǎn)振動頻率常數(shù),取決于邊界支撐條件。同樣,對于雙金屬翅片管這樣的串、并聯(lián)剛度系統(tǒng)的組合截面,可以推導(dǎo)出其等效截面扭轉(zhuǎn)剛度為:

    (11)

    而其單位長度的等效轉(zhuǎn)動慣量為:

    (12)

    這樣,將式(11)、(12)代入式(10),可以解出雙金屬復(fù)合翅片管的各階扭轉(zhuǎn)振動固有頻率。

    3.2軸向振動模態(tài)解析

    在軸向載荷激勵下,雙金屬復(fù)合翅片管的軸向自由振動的動力學(xué)方程為:

    (13)

    式中,u(x,t)為距原點x處截面在t時刻的軸向位移。同樣,將EeAe定義為雙金屬復(fù)合翅片管的等效截面抗拉壓剛度,me為其單位長度等效質(zhì)量。求解可得雙金屬翅片管的各階軸向振動固有頻率為:

    (14)

    式中,λli(i=1, 2, 3,…)為各階軸向振動頻率常數(shù),取決于邊界支撐條件。

    同樣采用串、并聯(lián)剛度法,可以推導(dǎo)出雙金屬翅片管的等效截面抗拉壓剛度為:

    (15)

    將式(15)、(8)代入式(14) 可以獲得其各階軸向振動固有頻率。

    3.3扭轉(zhuǎn)和軸向振動頻率結(jié)果分析

    表3顯示了兩端自由條件下,鋼鋁雙金屬復(fù)合翅片管扭轉(zhuǎn)與軸向振動固有頻率的有限元模擬值及理論解析值的比較,其中管子幾何參數(shù)同表1。在表3中,誤差指的是理論解析值與螺旋狀模型有限元分析值之間的差別。從表3可以看出,① 槽狀模型與螺旋狀模型對于扭轉(zhuǎn)振動和軸向振動固有頻率的分析結(jié)果也非常接近,進一步表明對雙金屬螺旋翅片管振動模態(tài)進行分析時,將其簡化為槽狀模型是可行的。② 采用等效截面扭轉(zhuǎn)剛度GeIpe并結(jié)合等效轉(zhuǎn)動慣量Je,計算得到的雙金屬復(fù)合翅片管的前兩階扭轉(zhuǎn)振動固有頻率,與有限元值相比,一階頻率誤差為1.41%,二階頻率誤差為1.25%。這說明該解析方法具有較高精度。③ 采用等效截面抗拉壓剛度EeAe并結(jié)合等效質(zhì)量me,計算得到的雙金屬復(fù)合翅片管的軸向振動固有頻率,與有限元值相比,一階頻率誤差為3.87%,二階頻率誤差為31.4%。本文經(jīng)分析發(fā)現(xiàn),這個誤差產(chǎn)生的主要原因是,在管子軸向振動時,翅片會同時沿軸向產(chǎn)生彎曲振動(相當于環(huán)形薄板的彎曲振動),二者耦合作用會對管子的軸向振動頻率,尤其是高階頻率產(chǎn)生較大影響。這種影響可以從管子的軸向振動模態(tài)振型圖中得到進一步驗證,如圖5(h)和圖5(j)所示。

    圖5為雙金屬復(fù)合翅片管各階模態(tài)振型。從圖中可以看出,雙金屬翅片管的彎曲和扭轉(zhuǎn)模態(tài)振型與光滑管完全相同,但其軸向振動的模態(tài)振型卻有不同。由圖5 (h),(j)可以看出,在管子軸向振動時,翅片沿軸向發(fā)生倒伏現(xiàn)象,即翅片沿管子軸向發(fā)生彎曲振動,二者耦合作用會對翅片管的軸向振動頻率產(chǎn)生重要影響。因此,式(14)由于沒有考慮翅片沿軸向的彎曲振動而產(chǎn)生較大誤差。另外需要注意的是,管子軸向振動時,翅片會沿軸向出現(xiàn)雙向反復(fù)倒伏現(xiàn)象,這時翅片由于受到交變載荷作用容易產(chǎn)生疲勞破壞,即翅片會在根部斷裂而喪失傳熱作用。

    表3 雙金屬復(fù)合翅片管扭轉(zhuǎn)和軸向振動固有頻率

    圖5 雙金屬復(fù)合翅片管各階模態(tài)振型Fig.5 Modal shapes of bimetallic finned tube

    4不同翅片幾何參數(shù)時雙金屬復(fù)合翅片管的振動頻率

    換熱管設(shè)置翅片是為了增加傳熱面積。翅片幾何參數(shù),如翅片高度h、翅片厚度a和翅片間距b是影響傳熱面積的重要因素。為滿足不同的傳熱要求,工程中需設(shè)計不同的翅片幾何參數(shù), 而翅片幾何參數(shù)的變化會直接影響翅片管的剛度和質(zhì)量,從而改變整個管子的振動頻率。 表4~表6顯示了翅片幾何參數(shù)對鋼

    鋁雙金屬復(fù)合翅片管振動頻率的影響規(guī)律和影響程度,其中管子的基本幾何參數(shù)D1,D2,D3和D4見表1。

    由表4可知,① 對于鋼鋁雙金屬復(fù)合翅片管,當其翅片高度h由0 mm增大到12 mm時,在兩端自由邊界條件下,其一階彎曲固有頻率(以下均指螺旋狀模型的有限元分析值)由122.7 Hz降低到104.9 Hz,降低幅度為14.5%;其一階扭轉(zhuǎn)固有頻率由1 571.4 Hz降低到1 026.0 Hz,降低幅度為34.6%;其一階軸向固有頻率由2 545.1 Hz降低到1 794.2 Hz,降低幅度為29.5%。②對不同的翅片高度,一階彎曲和扭轉(zhuǎn)振動頻率的解析值與有限元值均非常吻合,最大誤差僅為1.57%。③ 一階軸向振動頻率的解析值與有限元值相比,其誤差隨翅片高度的增加而增大。當翅片高度h小于或等于9.0 mm時,誤差最大值為3.87%;當翅片高度h為12 mm時,誤差達到21.71%。這進一步表明翅片沿軸向的彎曲振動會對管子的軸向振動頻率產(chǎn)生大的影響,且翅片越高,影響越大。但在工程實際中,翅片太高容易在翅片間隙積塵,不利于傳熱,因此實際翅片高度會有所限制。

    由表5可知,① 當翅片厚度a由0.35 mm增大到0.80 mm時,其一階彎曲頻率的降低幅度為8.6%,一階扭轉(zhuǎn)頻率的降低幅度為15.9%,一階軸向頻率的降低幅度為6.4%。② 對于不同的翅片厚度,管子的一階彎曲、扭轉(zhuǎn)和軸向振動頻率的解析值均與有限元值較為吻合,一階軸向頻率的誤差最大,為3.87%。

    由表6可知,① 當翅片間距b由2.8 mm增大到4.75 mm時,其一階彎曲頻率增大4.2%,一階扭轉(zhuǎn)頻率增大9.9%,一階軸向頻率增大7.0%。② 對不同的翅片間距,管子的一階彎曲、扭轉(zhuǎn)和軸向振動頻率的解析值也與有限元值非常吻合,最大誤差仍為軸向頻率的誤差,為3.87%。

    另外,由式(2)、(10)和(14)可知,翅片幾何參數(shù)對雙金屬翅片管振動頻率的影響規(guī)律在不同邊界支撐條件下是相同的,在高階頻率下也是相同的。因此可以說,對于矩形翅片的鋼鋁雙金屬翅片管,在任意邊界條件下,管子的各階彎曲、扭轉(zhuǎn)和軸向振動頻率均隨翅片高度h和翅片厚度a的增大而減小,隨翅片間距b的增大而增大。

    表4 不同翅片高度時雙金屬翅片管的固有頻率 (a=0.35 mm,b=2.80 mm)

    表5 不同翅片厚度時雙金屬翅片管的固有頻率(h=9.0 mm,b=2.8 mm)

    表6 不同翅片間距時雙金屬翅片管的固有頻率(h=9.0 mm,a=0.35 mm)

    5結(jié)論

    (1) 將雙金屬復(fù)合翅片管結(jié)構(gòu)簡化為串、并聯(lián)剛度系統(tǒng),采用組合截面等效彎曲剛度、等效扭轉(zhuǎn)剛度和等效抗拉壓剛度,并結(jié)合等效質(zhì)量和等效轉(zhuǎn)動慣量的方法,對其彎曲、扭轉(zhuǎn)和軸向振動模態(tài)進行理論解析。其中,彎曲和扭轉(zhuǎn)振動頻率的理論解析值與有限元分析值相比,誤差較??;而軸向振動頻率的誤差相對較大,主要原因是沒有考慮翅片沿管子軸向的彎曲振動對管子軸向振動的耦合作用。

    (2) 翅片幾何參數(shù)對雙金屬復(fù)合翅片管的振動頻率有重要影響。對于矩形翅片的鋼鋁雙金屬翅片管,管子的各階彎曲、扭轉(zhuǎn)和軸向振動頻率均隨翅片高度h和翅片厚度a的增大而減小,隨翅片間距b的增大而增大。

    (3) 對雙金屬螺旋翅片管振動模態(tài)進行有限元模擬時,將其簡化為槽狀結(jié)構(gòu)產(chǎn)生的誤差很小。

    參 考 文 獻

    [ 1 ] Goldstein R J, Ibele W E, Patankar S V. Heat transfer—A review of 2004 literature [J]. International Journal of Heat and Mass Transfer, 2010, 53:4343-4396.

    [ 2 ] Kuntysh V B, Sankovich E S, Mulin V P. An energy-saving nondestructive tester for the thermal performance of bimetallic finned tubes [J]. Chemical and Petroleum Engineering, 2007, 43(3/4):227-233.

    [ 4 ] 楊小瓊, 王啟杰, 王曉軍.翅片管接觸熱阻的理論預(yù)測和試驗研究[J]. 化工學(xué)報, 1996, 47(4): 427-432.

    YANG Xiao-qiong, WANG Qi-jie, WANG Xiao-jun. heoretical prediction and experimental investigation of finned tubes thermal contact resistance [J]. Journal of Chemical Industry and Engineering, 1996, 47(4): 427-432.

    [ 5 ] 劉偉軍, 付國民, 張書華. Al—Cu翅片管的單管換熱及阻力特性[J]. 哈爾濱理工大學(xué)學(xué)報, 1996, 1(1): 17-20.

    LIU Wei-jun, FU Guo-ming, ZHANG Shu-hua. Convective transfer and resistance characteristics of single aluminum-copper fin-tube [J]. Journal of Harbin University of Science and Technology, 1996, 1(1): 17-20.

    [ 6 ] Longatte E, Baj F, Hoarau Y, et al. Advanced numerical methods for uncertainty reduction when predicting heat exchanger dynamic stability limits: Review and perspectives [J]. Nuclear Engineering and Design, 2013, 258:164-175.

    [ 7 ] Goyder H G D. Flow induced vibration in heat exchangers [J]. Institution of Chemical Engineers, 2001, 80, Part A: 226-232.

    [ 8 ] GB151-1999. 管殼式換熱器[S]. 北京: 中國標準出版社, 1999.

    [ 9 ] 謝志勇, 周其斗, 紀剛. 雙層柱殼的流固耦合模態(tài)計算與試驗研究[J]. 海軍工程大學(xué)學(xué)報, 2009, 21(2): 97-101.

    XIE Zhi-yong, ZHOU Qi-dou, JI Gang.Computation and measurement of double shell vibration mode with f luid load [J]. Journal of Naval University of Engineering, 2009, 21(2): 97-101.

    [10] 金廣文, 何琳, 姜榮俊. 流固耦合對雙層圓柱殼體振動特性的影響[J]. 武漢理工大學(xué)學(xué)報:交通科學(xué)與工程版, 2007, 31(5): 882-885.

    JIN Guang-wen,HE Lin,JIANG Rong-jun. Influence of fluid-structure coupling for vibration transmissibility of ribbed cylindrical double-shell [J]. Journal of Wuhan University of Technology:Transportation Science & Engineering, 2007, 31(5): 882-885.

    [11] 魏建輝, 陳美霞, 喬志李,等. 湍流激勵下單雙層圓柱殼振動特性對比[J]. 華中科技大學(xué)學(xué)報:自然科學(xué)版, 2012, 40(8): 44-48.

    WEI Jian-hui, CHEN Mei-xia, QIAO Zhi-li, et al. Comparing the vibration characteristics of single cylindrical shells and double ones excited by turbulent boundary layers [J]. Journal of Huazhong University of Science and Technology:Natural Science Edition, 2012, 40(8): 44-48.

    Vibration characteristics of bimetallic finned tubes

    GUOJian-ying1,MATeng-fei2,LIUSheng-bao3,BAIYan-yan4(1. College of Mining Engineering, Taiyuan University of Technology, Taiyuan 030024, China; 2. Institute of Applied Mechanics and Biomedical Engineering, Taiyuan University of Technology, Taiyuan 030024, China; 3. Sinopec Ningbo Engineering Co., Ltd., Ningbo 315103, China; 4. College of Mechanical Engineering, Taiyuan University of Technology, Taiyuan 030024, China)

    Abstract:Bimetallic finned tubes are highly efficient heat transfer components, and they are easy to have flow-induced vibration damages. Their modal analysis has an important engineering significance. A bimetallic finned tube was simplified as a system with series and parallel stiffnesses according to its structural characteristics. The bending, torsion and axial vibration modes of the bimetallic finned tube were analyzed theoretically by using the methods of equivalent cross-sectional bending stiffness, equivalent torsion stiffness and equivalent tension-compression stiffness combined with the equivalent mass and equivalent moment of inertia. The tests and finite element analysis for vibration modes of bimetallic finned tubes were conducted in order to verify the correctness of the theoretical methods. The effects of fin geometry on the natural frequencies of a steel-aluminum bimetal finned tube with rectangular fins were investigated. The results showed that its bending, torsional and axial vibration natural frequencies decrease with increase in fin height and fin thickness, and increase with increase in fin pitch.

    Key words:bimetallic finned tubes; vibration modes; equivalent cross-sectional stiffness; equivalent mass; equivalent moment of inertia

    中圖分類號:TB123

    文獻標志碼:A

    DOI:10.13465/j.cnki.jvs.2016.05.033

    收稿日期:2015-02-02修改稿收到日期:2015-03-15

    基金項目:山西省自然科學(xué)基金(2013011025-2)

    第一作者 郭建英 女,博士,副教授,1972年11月生

    日韩av不卡免费在线播放| 菩萨蛮人人尽说江南好唐韦庄| 菩萨蛮人人尽说江南好唐韦庄| 精品少妇久久久久久888优播| 亚洲国产色片| 国产日韩欧美在线精品| 国产精品久久久久久久电影| 久久精品国产自在天天线| 成人特级av手机在线观看| 亚洲av免费高清在线观看| 亚洲精品国产av蜜桃| 成人一区二区视频在线观看| av.在线天堂| 亚洲欧美成人综合另类久久久| 亚洲伊人久久精品综合| 高清不卡的av网站| 婷婷色av中文字幕| 国产精品99久久久久久久久| 国产精品国产三级专区第一集| 久久久国产一区二区| 插阴视频在线观看视频| 久久久a久久爽久久v久久| 国产欧美亚洲国产| 久久精品久久久久久噜噜老黄| 欧美日韩视频高清一区二区三区二| 一级黄片播放器| 男人添女人高潮全过程视频| 交换朋友夫妻互换小说| 丰满人妻一区二区三区视频av| 91精品伊人久久大香线蕉| 成人一区二区视频在线观看| 精品久久久噜噜| 久热这里只有精品99| 亚洲经典国产精华液单| 寂寞人妻少妇视频99o| 边亲边吃奶的免费视频| 亚洲成人av在线免费| 国产白丝娇喘喷水9色精品| 国产精品一区二区三区四区免费观看| 国产精品麻豆人妻色哟哟久久| videos熟女内射| 一级毛片aaaaaa免费看小| 赤兔流量卡办理| 少妇高潮的动态图| 国产av精品麻豆| 男的添女的下面高潮视频| 中文欧美无线码| 免费播放大片免费观看视频在线观看| 大片电影免费在线观看免费| 国产精品熟女久久久久浪| 久久精品国产a三级三级三级| 人妻制服诱惑在线中文字幕| 观看av在线不卡| 偷拍熟女少妇极品色| 亚洲精品成人av观看孕妇| 日韩一区二区视频免费看| 日本免费在线观看一区| 亚洲一级一片aⅴ在线观看| 久久午夜福利片| 亚洲美女搞黄在线观看| 亚洲成色77777| 亚洲av成人精品一二三区| 欧美zozozo另类| 一级黄片播放器| 国产精品伦人一区二区| 老女人水多毛片| 一级爰片在线观看| 久久久欧美国产精品| 久久久久久久久久久免费av| 我的老师免费观看完整版| 中文欧美无线码| 成人国产av品久久久| 亚洲精品成人av观看孕妇| 伊人久久国产一区二区| 精品国产露脸久久av麻豆| 99久久精品热视频| 毛片女人毛片| 久久这里有精品视频免费| 交换朋友夫妻互换小说| 亚洲精品国产色婷婷电影| 国产成人aa在线观看| 久久久精品免费免费高清| 亚洲av日韩在线播放| 自拍欧美九色日韩亚洲蝌蚪91 | 国产久久久一区二区三区| 黄色一级大片看看| www.色视频.com| 欧美成人一区二区免费高清观看| 国产大屁股一区二区在线视频| 欧美日韩综合久久久久久| 日本vs欧美在线观看视频 | 欧美高清性xxxxhd video| 亚洲欧美日韩东京热| 哪个播放器可以免费观看大片| 校园人妻丝袜中文字幕| 99久久人妻综合| 网址你懂的国产日韩在线| 视频区图区小说| 六月丁香七月| 观看免费一级毛片| 国产乱人视频| 日韩成人av中文字幕在线观看| 亚洲成人一二三区av| 99国产精品免费福利视频| 精品人妻视频免费看| 国产黄片视频在线免费观看| 国产v大片淫在线免费观看| 久久鲁丝午夜福利片| 最黄视频免费看| 欧美zozozo另类| 国产在线一区二区三区精| 人体艺术视频欧美日本| 亚洲欧洲日产国产| 超碰97精品在线观看| 不卡视频在线观看欧美| 久久久久精品久久久久真实原创| 最近手机中文字幕大全| 在线观看免费视频网站a站| 亚洲精品亚洲一区二区| 日韩强制内射视频| 99热这里只有精品一区| 午夜福利在线观看免费完整高清在| xxx大片免费视频| 国产精品久久久久久精品电影小说 | 成人免费观看视频高清| 色吧在线观看| 久久久色成人| h视频一区二区三区| 男人舔奶头视频| 丰满迷人的少妇在线观看| 亚洲精品中文字幕在线视频 | 美女脱内裤让男人舔精品视频| 五月玫瑰六月丁香| 性色av一级| 精品少妇黑人巨大在线播放| 国产精品国产三级国产av玫瑰| 久久久色成人| 蜜桃在线观看..| 国产真实伦视频高清在线观看| 高清不卡的av网站| h视频一区二区三区| 只有这里有精品99| videossex国产| 国产高清三级在线| 精品人妻视频免费看| 精品久久久久久久末码| 亚洲三级黄色毛片| 亚洲欧美日韩无卡精品| 亚洲精品久久午夜乱码| 久久久a久久爽久久v久久| 91久久精品电影网| 多毛熟女@视频| 久久人人爽人人片av| 色5月婷婷丁香| 最黄视频免费看| 国产在视频线精品| 成人午夜精彩视频在线观看| 国产一区二区在线观看日韩| 老司机影院成人| 国产精品国产av在线观看| 国产精品成人在线| 一级毛片我不卡| 大又大粗又爽又黄少妇毛片口| 日产精品乱码卡一卡2卡三| 欧美少妇被猛烈插入视频| 国产精品女同一区二区软件| 1000部很黄的大片| 国内少妇人妻偷人精品xxx网站| 国产精品一及| 狂野欧美激情性xxxx在线观看| 久久人人爽av亚洲精品天堂 | 如何舔出高潮| 尤物成人国产欧美一区二区三区| 久久亚洲国产成人精品v| 国产免费一级a男人的天堂| 精品一区二区免费观看| 色视频在线一区二区三区| 在线 av 中文字幕| 直男gayav资源| 久久鲁丝午夜福利片| 日本av免费视频播放| 激情五月婷婷亚洲| 日本一二三区视频观看| 久久精品国产a三级三级三级| 黄色配什么色好看| 亚洲激情五月婷婷啪啪| 国产亚洲欧美精品永久| 国产在视频线精品| 中文字幕人妻熟人妻熟丝袜美| 天堂中文最新版在线下载| 精品一区在线观看国产| 免费黄频网站在线观看国产| 国产精品福利在线免费观看| 亚洲va在线va天堂va国产| 免费看不卡的av| av国产免费在线观看| 国产精品嫩草影院av在线观看| 亚洲欧美成人精品一区二区| 18禁裸乳无遮挡动漫免费视频| 久久精品国产亚洲网站| 身体一侧抽搐| 一级二级三级毛片免费看| 极品少妇高潮喷水抽搐| 日本午夜av视频| 国产成人免费无遮挡视频| 日韩电影二区| 男女啪啪激烈高潮av片| 久久人人爽人人爽人人片va| 亚洲精品,欧美精品| 天堂中文最新版在线下载| 日韩欧美 国产精品| 午夜福利在线在线| 成人特级av手机在线观看| 久久久久性生活片| 国产精品不卡视频一区二区| 日产精品乱码卡一卡2卡三| 我要看日韩黄色一级片| 下体分泌物呈黄色| 一级爰片在线观看| 男女免费视频国产| 国产精品免费大片| 美女cb高潮喷水在线观看| 精品亚洲成a人片在线观看 | 亚洲国产欧美人成| 最近手机中文字幕大全| 国产色婷婷99| 高清黄色对白视频在线免费看 | 老司机影院成人| 日本黄色日本黄色录像| 久久久久人妻精品一区果冻| 中文在线观看免费www的网站| freevideosex欧美| 国产精品嫩草影院av在线观看| 十八禁网站网址无遮挡 | 免费观看的影片在线观看| 亚洲久久久国产精品| 久久久久久伊人网av| 国产乱来视频区| 大陆偷拍与自拍| a级毛片免费高清观看在线播放| 精品亚洲乱码少妇综合久久| 精品人妻视频免费看| 国产片特级美女逼逼视频| 大香蕉久久网| 国产精品免费大片| 在线观看一区二区三区激情| 国内揄拍国产精品人妻在线| 韩国高清视频一区二区三区| 欧美精品一区二区大全| 18禁裸乳无遮挡免费网站照片| 欧美一级a爱片免费观看看| 看非洲黑人一级黄片| a 毛片基地| 久久久久久久久久人人人人人人| 亚洲色图综合在线观看| 国产毛片在线视频| 日韩强制内射视频| 在线观看免费视频网站a站| 久久99热这里只有精品18| 97精品久久久久久久久久精品| 女人久久www免费人成看片| 夫妻午夜视频| 老师上课跳d突然被开到最大视频| 一区二区三区精品91| 亚洲精品国产av成人精品| 免费看光身美女| 一边亲一边摸免费视频| 极品少妇高潮喷水抽搐| 国产精品av视频在线免费观看| 久久久久国产精品人妻一区二区| 国产成人精品久久久久久| 亚洲色图综合在线观看| 欧美激情国产日韩精品一区| 国产精品一区二区在线不卡| 亚洲三级黄色毛片| 国产一区亚洲一区在线观看| 亚洲国产精品国产精品| 少妇人妻 视频| 性色av一级| 国产爱豆传媒在线观看| 久久久久网色| 看免费成人av毛片| 国产亚洲最大av| 国产色爽女视频免费观看| 日韩伦理黄色片| 欧美xxxx性猛交bbbb| 精品国产一区二区三区久久久樱花 | 午夜老司机福利剧场| 久久久久久久久久成人| 国产精品一区二区在线不卡| 国产在线男女| 亚洲av综合色区一区| 精品国产一区二区三区久久久樱花 | 五月天丁香电影| 人妻制服诱惑在线中文字幕| 精品少妇久久久久久888优播| av视频免费观看在线观看| 久久久久视频综合| 在线观看人妻少妇| 3wmmmm亚洲av在线观看| 欧美 日韩 精品 国产| 22中文网久久字幕| 国产 精品1| 精品国产露脸久久av麻豆| 高清不卡的av网站| 高清视频免费观看一区二区| 超碰av人人做人人爽久久| 亚洲丝袜综合中文字幕| 亚洲内射少妇av| tube8黄色片| 日韩制服骚丝袜av| 久久久久久久久久久丰满| 国产 一区精品| 亚洲人与动物交配视频| 日韩人妻高清精品专区| 久久久久精品性色| 在线看a的网站| 卡戴珊不雅视频在线播放| 国产免费又黄又爽又色| 黄色欧美视频在线观看| 久热久热在线精品观看| 亚洲人成网站高清观看| 亚洲欧美精品专区久久| 在现免费观看毛片| 免费大片黄手机在线观看| 九九久久精品国产亚洲av麻豆| 国产片特级美女逼逼视频| 久久婷婷青草| 蜜臀久久99精品久久宅男| 精品一区在线观看国产| 亚洲av综合色区一区| 亚洲色图av天堂| 国产精品人妻久久久影院| 国产精品99久久99久久久不卡 | 交换朋友夫妻互换小说| 在线观看美女被高潮喷水网站| 成年av动漫网址| 肉色欧美久久久久久久蜜桃| 黑人猛操日本美女一级片| 中文字幕亚洲精品专区| 成人免费观看视频高清| 又粗又硬又长又爽又黄的视频| 国产乱来视频区| 黄色日韩在线| 久久精品久久精品一区二区三区| 国产欧美亚洲国产| 91精品国产国语对白视频| 免费大片黄手机在线观看| 综合色丁香网| 97超碰精品成人国产| 国产男女内射视频| 大又大粗又爽又黄少妇毛片口| 亚洲精品色激情综合| 国产爽快片一区二区三区| 午夜日本视频在线| 国产一区亚洲一区在线观看| 欧美精品一区二区免费开放| 国产视频首页在线观看| 夜夜骑夜夜射夜夜干| 久久久久久久亚洲中文字幕| 久久国产精品大桥未久av | 日韩伦理黄色片| 在线 av 中文字幕| 亚洲欧美日韩东京热| 国产男人的电影天堂91| 国产精品一二三区在线看| 成人亚洲欧美一区二区av| 免费久久久久久久精品成人欧美视频 | a级毛片免费高清观看在线播放| 久久精品国产亚洲av天美| 日韩av在线免费看完整版不卡| 国产精品熟女久久久久浪| 免费观看在线日韩| 国产精品久久久久久精品电影小说 | 人人妻人人添人人爽欧美一区卜 | 一边亲一边摸免费视频| 80岁老熟妇乱子伦牲交| 青春草亚洲视频在线观看| 人人妻人人看人人澡| 狂野欧美激情性xxxx在线观看| 精品一品国产午夜福利视频| 久久久成人免费电影| 中文字幕精品免费在线观看视频 | 免费久久久久久久精品成人欧美视频 | 午夜福利网站1000一区二区三区| 少妇人妻一区二区三区视频| 国产精品女同一区二区软件| 2021少妇久久久久久久久久久| 尤物成人国产欧美一区二区三区| 精品酒店卫生间| 亚洲伊人久久精品综合| 国产精品嫩草影院av在线观看| 18+在线观看网站| 亚洲av成人精品一二三区| 偷拍熟女少妇极品色| 成年av动漫网址| 岛国毛片在线播放| 一本—道久久a久久精品蜜桃钙片| 建设人人有责人人尽责人人享有的 | 国产在线男女| 看非洲黑人一级黄片| 国产片特级美女逼逼视频| 成人免费观看视频高清| 丰满乱子伦码专区| 亚洲国产成人一精品久久久| 人妻夜夜爽99麻豆av| 91精品国产国语对白视频| 久久毛片免费看一区二区三区| 色婷婷av一区二区三区视频| 一级黄片播放器| 亚洲电影在线观看av| 日韩成人av中文字幕在线观看| 国产淫语在线视频| 亚洲av国产av综合av卡| 亚洲婷婷狠狠爱综合网| 国产成人一区二区在线| 亚洲国产欧美在线一区| av在线老鸭窝| 一边亲一边摸免费视频| 亚洲无线观看免费| 国产成人a∨麻豆精品| 精品国产一区二区三区久久久樱花 | 亚洲av不卡在线观看| 欧美bdsm另类| 欧美成人午夜免费资源| 国产成人精品久久久久久| 老司机影院成人| 日韩一本色道免费dvd| 亚洲精品日本国产第一区| 天天躁夜夜躁狠狠久久av| 狂野欧美激情性bbbbbb| 深夜a级毛片| 内地一区二区视频在线| 一级毛片 在线播放| 亚洲人成网站高清观看| 亚洲精品国产成人久久av| 麻豆成人av视频| 91在线精品国自产拍蜜月| 国产成人精品福利久久| 午夜福利视频精品| 日本免费在线观看一区| 国产淫片久久久久久久久| 亚洲av中文av极速乱| 国产精品久久久久久精品古装| 国产 一区精品| 亚洲不卡免费看| 99热网站在线观看| 国产免费又黄又爽又色| 97在线视频观看| 人人妻人人爽人人添夜夜欢视频 | 日本一二三区视频观看| 日本爱情动作片www.在线观看| 建设人人有责人人尽责人人享有的 | 一级毛片 在线播放| 秋霞在线观看毛片| 欧美最新免费一区二区三区| 免费观看性生交大片5| 亚洲精品一二三| 大码成人一级视频| 中文字幕人妻熟人妻熟丝袜美| 少妇人妻 视频| 久久久久久久大尺度免费视频| 成人毛片60女人毛片免费| 国产精品爽爽va在线观看网站| 在线观看一区二区三区激情| 少妇人妻一区二区三区视频| 97热精品久久久久久| 伊人久久精品亚洲午夜| 一级毛片电影观看| 秋霞伦理黄片| 亚洲一级一片aⅴ在线观看| 直男gayav资源| 少妇被粗大猛烈的视频| 国产69精品久久久久777片| 久久久久久九九精品二区国产| 日本欧美国产在线视频| .国产精品久久| 91精品国产国语对白视频| 人人妻人人爽人人添夜夜欢视频 | 欧美精品一区二区免费开放| 麻豆乱淫一区二区| 日韩一本色道免费dvd| 国产无遮挡羞羞视频在线观看| 亚洲欧美成人综合另类久久久| 成人美女网站在线观看视频| 久久ye,这里只有精品| 精品酒店卫生间| 国内精品宾馆在线| 中文欧美无线码| 久久精品国产鲁丝片午夜精品| 亚洲内射少妇av| 欧美老熟妇乱子伦牲交| 精品一区在线观看国产| 日韩一区二区视频免费看| 色网站视频免费| 日本vs欧美在线观看视频 | 久久 成人 亚洲| 亚洲丝袜综合中文字幕| 精品午夜福利在线看| 国产大屁股一区二区在线视频| 菩萨蛮人人尽说江南好唐韦庄| 在线亚洲精品国产二区图片欧美 | 免费黄频网站在线观看国产| 少妇被粗大猛烈的视频| 极品教师在线视频| av国产免费在线观看| 国产v大片淫在线免费观看| 国产精品无大码| 国产精品久久久久久精品古装| 黄片无遮挡物在线观看| 99热全是精品| 自拍偷自拍亚洲精品老妇| 成年免费大片在线观看| 欧美精品一区二区免费开放| 建设人人有责人人尽责人人享有的 | 欧美精品一区二区免费开放| 啦啦啦啦在线视频资源| 99久久中文字幕三级久久日本| 99久久精品一区二区三区| 亚洲,一卡二卡三卡| 国产精品一区二区性色av| 亚洲av综合色区一区| 亚洲精品乱久久久久久| 伦精品一区二区三区| 国产在视频线精品| 亚洲欧美精品专区久久| 国产高清三级在线| 午夜免费鲁丝| 免费观看av网站的网址| 欧美3d第一页| 精品午夜福利在线看| 久久精品国产亚洲av涩爱| 毛片女人毛片| 在线观看美女被高潮喷水网站| av天堂中文字幕网| 免费在线观看成人毛片| 国产精品一及| 高清视频免费观看一区二区| 毛片女人毛片| av免费观看日本| 亚洲三级黄色毛片| 亚洲aⅴ乱码一区二区在线播放| 交换朋友夫妻互换小说| 成年av动漫网址| 多毛熟女@视频| 三级经典国产精品| 一级黄片播放器| 日日撸夜夜添| 色5月婷婷丁香| 中文字幕精品免费在线观看视频 | 欧美成人精品欧美一级黄| 中文精品一卡2卡3卡4更新| 亚洲欧美日韩东京热| 国产淫片久久久久久久久| 精品国产乱码久久久久久小说| 精品一品国产午夜福利视频| 美女福利国产在线 | 国产成人a∨麻豆精品| 男女啪啪激烈高潮av片| 极品少妇高潮喷水抽搐| 美女xxoo啪啪120秒动态图| 欧美日韩综合久久久久久| 少妇人妻精品综合一区二区| 亚洲美女视频黄频| 97超碰精品成人国产| 国产在视频线精品| 成人黄色视频免费在线看| 青春草视频在线免费观看| 久久97久久精品| 日韩人妻高清精品专区| 极品少妇高潮喷水抽搐| 天堂8中文在线网| 在线天堂最新版资源| 男人狂女人下面高潮的视频| 午夜激情久久久久久久| 看十八女毛片水多多多| 韩国av在线不卡| 18禁裸乳无遮挡动漫免费视频| 十八禁网站网址无遮挡 | 联通29元200g的流量卡| 91精品伊人久久大香线蕉| 中文字幕av成人在线电影| 尤物成人国产欧美一区二区三区| 免费大片18禁| 交换朋友夫妻互换小说| 精品亚洲乱码少妇综合久久| 免费大片18禁| 日韩欧美 国产精品| 久久精品国产a三级三级三级| 成人18禁高潮啪啪吃奶动态图 | 日日啪夜夜爽| 少妇裸体淫交视频免费看高清| 亚洲第一区二区三区不卡| 中文字幕亚洲精品专区| 伦理电影免费视频| 尾随美女入室| 在线 av 中文字幕| 亚洲国产av新网站| 亚洲,欧美,日韩| av不卡在线播放| 精品久久久久久久久亚洲| 国产午夜精品一二区理论片| 日本一二三区视频观看| 日韩伦理黄色片| 日韩在线高清观看一区二区三区| 亚洲精品国产av成人精品| 亚洲精品自拍成人| 特大巨黑吊av在线直播| 日韩三级伦理在线观看| 人妻系列 视频| 久久精品国产亚洲av涩爱| 久久精品久久精品一区二区三区| 五月开心婷婷网| 亚洲综合精品二区| 高清黄色对白视频在线免费看 |