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

    高強(qiáng)化柴油機(jī)鑄鐵缸蓋承載特性研究

    2017-03-08 05:22:01景國(guó)璽張樹勇付文清許春光文洋李鵬衛(wèi)軍朝
    車用發(fā)動(dòng)機(jī) 2017年1期
    關(guān)鍵詞:缸蓋熱應(yīng)力火力

    景國(guó)璽, 張樹勇, 付文清, 許春光, 文洋, 李鵬, 衛(wèi)軍朝

    (1. 中國(guó)北方發(fā)動(dòng)機(jī)研究所(天津), 天津 300400; 2. 中國(guó)兵器科學(xué)研究院, 北京 100089)

    高強(qiáng)化柴油機(jī)鑄鐵缸蓋承載特性研究

    景國(guó)璽1, 張樹勇1, 付文清2, 許春光1, 文洋1, 李鵬1, 衛(wèi)軍朝1

    (1. 中國(guó)北方發(fā)動(dòng)機(jī)研究所(天津), 天津 300400; 2. 中國(guó)兵器科學(xué)研究院, 北京 100089)

    研究了某柴油機(jī)鑄鐵缸蓋結(jié)構(gòu)溫度和應(yīng)力隨時(shí)間的變化行為,明確熱-機(jī)械載荷作用下缸蓋的承載規(guī)律,為缸蓋壽命預(yù)測(cè)模型建立提供依據(jù)。首先建立了缸蓋有限元分析模型,利用實(shí)測(cè)溫度和應(yīng)力數(shù)據(jù)對(duì)模型進(jìn)行了標(biāo)定,進(jìn)而基于該模型計(jì)算了標(biāo)定工況及怠速與標(biāo)定工況交替變化條件下缸蓋溫度和應(yīng)力的變化規(guī)律。仿真結(jié)果表明:缸蓋最高溫度和最大熱應(yīng)力出現(xiàn)在鼻梁區(qū)域;在工作循環(huán)內(nèi),缸蓋火力面溫度波動(dòng)幅值在30 ℃以內(nèi),由此引起的波動(dòng)熱應(yīng)力相對(duì)定常熱應(yīng)力較小,但相對(duì)高頻氣體應(yīng)力較大,故在火力面高周疲勞校核過程中必須考慮其帶來的影響;在怠速和標(biāo)定轉(zhuǎn)速交替變化工況下,鼻梁區(qū)載荷變化最明顯,溫度與應(yīng)力呈現(xiàn)反相狀態(tài),應(yīng)力幅值較高,該區(qū)域易發(fā)生低周疲勞損傷。

    柴油機(jī); 氣缸蓋; 仿真; 承載特性; 溫度; 應(yīng)力

    缸蓋是發(fā)動(dòng)機(jī)中結(jié)構(gòu)最復(fù)雜、承受機(jī)械載荷和熱載荷最嚴(yán)重的部件之一,其承載特性是影響缸蓋疲勞強(qiáng)度可靠性的重要因素[1-3]。缸蓋不僅承受螺栓預(yù)緊力的強(qiáng)約束作用,而且承受交變氣體力和熱應(yīng)力的耦合作用,同時(shí)由于發(fā)動(dòng)機(jī)服役工況范圍廣且多變,缸蓋承載規(guī)律極其復(fù)雜。研究表明,穩(wěn)態(tài)工況下螺栓預(yù)緊力、高頻溫度載荷與氣體力作用易導(dǎo)致缸蓋發(fā)生高周疲勞失效,失效部位主要出現(xiàn)在水腔隔板、氣道壁等區(qū)域。Mario Metzger的研究表明,工況大幅波動(dòng)(如停車—起動(dòng)—停車工況)易產(chǎn)生低頻波動(dòng)熱應(yīng)力,往往導(dǎo)致火力面鼻梁區(qū)發(fā)生低周疲勞失效[3]。目前,隨著發(fā)動(dòng)機(jī)功率密度的不斷提升,缸蓋承受的熱機(jī)械負(fù)荷增加,疲勞強(qiáng)度是當(dāng)前設(shè)計(jì)關(guān)注的焦點(diǎn)[4-5],因此,本研究重點(diǎn)探討了缸蓋溫度和應(yīng)力隨時(shí)間和載荷的變化規(guī)律,為缸蓋疲勞壽命預(yù)測(cè)模型建立提供依據(jù),為材料研制提供約束邊界。

    1 研究思路

    缸蓋必須同時(shí)滿足高周疲勞和低周疲勞的設(shè)計(jì)要求。高周疲勞強(qiáng)度按照無限壽命設(shè)計(jì),低周疲勞壽命與發(fā)動(dòng)機(jī)服役工況、服役壽命及用途相關(guān)。目前,主要基于發(fā)動(dòng)機(jī)臺(tái)架耐久性試驗(yàn)來考核發(fā)動(dòng)機(jī)結(jié)構(gòu)件疲勞特性,圖1示出發(fā)動(dòng)機(jī)典型臺(tái)架耐久規(guī)范的簡(jiǎn)化或等效示意,試驗(yàn)工況主要由標(biāo)定轉(zhuǎn)速穩(wěn)態(tài)工況及怠速與標(biāo)定交變工況兩部分組成,分別考核結(jié)構(gòu)件高周疲勞(HCF)和低周疲勞(LCF)性能。

    在圖1試驗(yàn)工況下,缸蓋主要承受以下作用力。1)預(yù)裝配應(yīng)力σf:主要由螺栓預(yù)緊力和氣門導(dǎo)管過盈引起;2)氣體應(yīng)力σg:主要由缸內(nèi)高頻波動(dòng)的氣體壓力引起;3)定常熱應(yīng)力σst:發(fā)動(dòng)機(jī)未工作前,結(jié)構(gòu)溫度較低且均勻,在高周疲勞考核穩(wěn)定工況(如標(biāo)定轉(zhuǎn)速工況)下,結(jié)構(gòu)溫度升高,由于熱慣性效應(yīng)溫度基本保持不變,由此形成的熱應(yīng)力即為定常熱應(yīng)力,通常該應(yīng)力已考慮預(yù)裝配應(yīng)力σf;4)高頻波動(dòng)熱應(yīng)力σht:在穩(wěn)態(tài)工況一個(gè)工作循環(huán)內(nèi),燃燒室氣體溫度、壓力和換熱邊界隨曲軸轉(zhuǎn)角高頻波動(dòng),引起缸蓋火力面表面溫度瞬態(tài)波動(dòng),并由此產(chǎn)生了高頻波動(dòng)的熱應(yīng)力;5)低頻波動(dòng)熱應(yīng)力σlt:在低周熱機(jī)疲勞考核工況下,缸蓋溫度隨著工況變化而發(fā)生大幅度波動(dòng),缸蓋處于急劇加熱和冷卻狀態(tài),由此產(chǎn)生了較大的低頻熱應(yīng)力。

    本研究主要開展以下工作:1)研究標(biāo)定轉(zhuǎn)速點(diǎn)穩(wěn)態(tài)工況下預(yù)裝配應(yīng)力σf、氣體應(yīng)力σg、定常熱應(yīng)力σst和高頻波動(dòng)熱應(yīng)力σht的分布特性和對(duì)缸蓋綜合應(yīng)力的影響規(guī)律;2)研究變工況條件下低頻波動(dòng)熱應(yīng)力σlt、溫度隨時(shí)間的變化關(guān)系;3)基于溫度和應(yīng)力分析結(jié)果,并結(jié)合發(fā)動(dòng)機(jī)臺(tái)架耐久規(guī)范進(jìn)一步構(gòu)建缸蓋疲勞壽命預(yù)測(cè)模型。

    圖1 發(fā)動(dòng)機(jī)臺(tái)架疲勞耐久試驗(yàn)考核工況示意

    主要研究思路見圖2。首先建立缸蓋有限元分析模型,利用實(shí)測(cè)溫度和應(yīng)力數(shù)據(jù)對(duì)模型進(jìn)行標(biāo)定,然后基于該模型計(jì)算標(biāo)定轉(zhuǎn)速工況穩(wěn)態(tài)條件下缸蓋溫度、應(yīng)力應(yīng)變行為,在此基礎(chǔ)上進(jìn)一步研究怠速工況與標(biāo)定工況交替變化條件下缸蓋載荷的變化特性。有限元網(wǎng)格模型見圖3,有限元模型包括缸蓋、進(jìn)排氣門、氣門導(dǎo)管、缸蓋螺栓、氣缸墊和氣缸體等部件,其中對(duì)缸蓋火力面等關(guān)鍵區(qū)域進(jìn)行了網(wǎng)格細(xì)化,模型共包含792 653個(gè)節(jié)點(diǎn),432 675個(gè)單元。

    圖2 總體研究思路

    圖3 有限元網(wǎng)格模型

    2 缸蓋溫度變化特性研究

    2.1 換熱邊界條件確定

    2.1.1 火力面熱邊界條件

    氣缸蓋火力面與高溫燃?xì)庵苯咏佑|,換熱過程以穩(wěn)態(tài)對(duì)流換熱為主。本研究基于AVLBOOST軟件采用Woshini模型計(jì)算得到了一個(gè)工作循環(huán)中的瞬時(shí)傳熱系數(shù)和燃?xì)饨橘|(zhì)溫度曲線,并基于式(2)得到了時(shí)間均化的當(dāng)量傳熱系數(shù)和燃?xì)鉁囟?,分別為1 703W/(m2·K)和840 ℃。

    ,

    (1)

    (2)

    式中:hgm和Tgm分別為時(shí)間均化的當(dāng)量傳熱系數(shù)和燃?xì)鉁囟取?/p>

    結(jié)合4氣門燃燒室空間結(jié)構(gòu),將火力面劃分為ININ,EXEX和EXIN3個(gè)區(qū)域,分別定義換熱邊界。各區(qū)傳熱系數(shù)沿著半徑方向的分布規(guī)律見圖4,圖中橫坐標(biāo)為實(shí)際位置相對(duì)缸徑的比值,縱坐標(biāo)為局部傳熱系數(shù)與當(dāng)量平均傳熱系數(shù)的比值。其中,ININ區(qū)域傳熱系數(shù)最小,EXEX區(qū)域傳熱系數(shù)最大,火力面熱邊界條件滿足式(3)要求。最后通過軟件編程實(shí)現(xiàn)了對(duì)火力面網(wǎng)格熱邊界條件的映射。

    圖4 火力面熱邊界條件分布規(guī)律

    (3)

    式中:r為氣缸半徑;b(r)為半徑r處的當(dāng)量傳熱系數(shù)。

    2.1.2 冷卻水腔表面熱邊界

    由于該發(fā)動(dòng)機(jī)采用并聯(lián)冷卻水腔設(shè)計(jì),各氣缸和單體式氣缸蓋采用單獨(dú)的進(jìn)出水口,因此僅對(duì)單缸冷卻水腔進(jìn)行分析。采用K-e模型對(duì)冷卻水進(jìn)行三維流動(dòng)數(shù)值模擬,入口邊界采用實(shí)測(cè)質(zhì)量流量(0.831 3kg/s),出口壓力為100kPa,壁面采用固定壁面邊界。通過計(jì)算獲得了缸體和缸蓋冷卻水腔表面?zhèn)鳠嵯禂?shù),并對(duì)有限元網(wǎng)格進(jìn)行了換熱邊界映射,其中缸蓋冷卻水腔映射前后傳熱系數(shù)見圖5。

    圖5 水腔表面?zhèn)鳠嵯禂?shù)

    2.1.3 其他區(qū)域熱邊界條件

    進(jìn)氣門與氣門座圈之間、排氣門與氣門座圈之間及氣門桿與氣門導(dǎo)管之間建立了接觸模型,接觸面間隙傳熱系數(shù)分別選取為10 000W/(m2·K),5 500W/(m2·K)和400W/(m2·K)[4]。其他區(qū)域換熱邊界采用經(jīng)驗(yàn)值。

    2.2 溫度場(chǎng)試驗(yàn)與模擬分析

    采用熱電偶方法進(jìn)行了缸蓋溫度場(chǎng)測(cè)試,熱電偶通過火力面鉆孔引出,測(cè)點(diǎn)布置見圖6。由于火力面輔助引線鉆孔直徑較大,阻礙火力面熱量向水側(cè)傳遞,對(duì)測(cè)溫結(jié)果有很大影響,因此,先基于試驗(yàn)測(cè)試用缸蓋結(jié)構(gòu)按照試驗(yàn)條件邊界進(jìn)行溫度場(chǎng)分析并與實(shí)測(cè)結(jié)果對(duì)比,確定火力面?zhèn)鳠嵯禂?shù)分布規(guī)律,進(jìn)而用于實(shí)際缸蓋溫度場(chǎng)預(yù)測(cè)。發(fā)動(dòng)機(jī)臺(tái)架試驗(yàn)工況條件下缸蓋溫度場(chǎng)測(cè)試和預(yù)測(cè)結(jié)果對(duì)比見圖7,由圖可知,實(shí)測(cè)最高溫度點(diǎn)位于兩個(gè)排氣門鼻梁區(qū),基于試驗(yàn)用缸蓋的多數(shù)測(cè)點(diǎn)實(shí)測(cè)值和計(jì)算值偏差在10 ℃以內(nèi),在相同熱邊界條件下實(shí)際缸蓋溫度偏低,峰值點(diǎn)溫度降低約40 ℃。

    圖6 熱電偶測(cè)點(diǎn)布置

    圖7 溫度場(chǎng)模擬結(jié)果與仿真結(jié)果對(duì)比

    在不考慮缸內(nèi)燃?xì)鉁囟群蛡鳠嵯禂?shù)隨時(shí)間變化對(duì)缸蓋整體溫度場(chǎng)影響的情況下,標(biāo)定工況點(diǎn)缸蓋溫度場(chǎng)模擬結(jié)果見圖8,高于300 ℃的區(qū)域主要位于與燃?xì)饨佑|的缸蓋底板區(qū)域,最高溫度為490 ℃,出現(xiàn)在兩排氣門鼻梁區(qū)域靠近氣門座孔邊緣處,缸蓋熱負(fù)荷極高。

    圖8 缸蓋溫度場(chǎng)分布

    2.3 工作循環(huán)內(nèi)溫度波動(dòng)分析

    溫度循環(huán)變化的燃?xì)饧捌渑c缸蓋的傳熱作用引起缸蓋火力面溫度的波動(dòng),研究表明[5]該溫度波動(dòng)向缸蓋內(nèi)部傳播,到達(dá)1~2mm時(shí)波動(dòng)幅度趨于零, 形成瞬變溫度層。三維瞬態(tài)溫度場(chǎng)可由下列微分方程式表達(dá):

    (4)

    初始條件

    (5)

    邊界條件

    (6)

    式中:T為溫度;τ為時(shí)間;k為導(dǎo)熱系數(shù);ρ為密度;c為比熱容;TW為壁面溫度;Tf為燃?xì)猱?dāng)量溫度;a為對(duì)流傳熱系數(shù)。

    矩陣常微分方程可描述為

    (7)

    本研究將前面穩(wěn)態(tài)溫度場(chǎng)計(jì)算結(jié)果作為瞬態(tài)熱分析的近似初始條件,進(jìn)而按照瞬態(tài)邊界計(jì)算了缸蓋在一個(gè)工作循環(huán)內(nèi)的溫度波動(dòng)。計(jì)算得到的缸蓋溫度波動(dòng)幅度見圖9,溫度波動(dòng)出現(xiàn)在靠近火力面約2mm的薄層內(nèi),標(biāo)定工況下溫度最大波動(dòng)幅度為26 ℃。計(jì)算得到的兩個(gè)排氣門之間特征區(qū)域溫度隨深度的變化關(guān)系見圖10,由圖可知,隨著深度增加,溫度波動(dòng)幅度逐漸減小,達(dá)到1mm后基本無波動(dòng)。

    圖9 缸蓋溫度場(chǎng)波動(dòng)幅

    圖10 火力面考察點(diǎn)溫度與深度的變化情況

    3 穩(wěn)態(tài)工況下載荷特性分析

    3.1 熱-機(jī)械耦合應(yīng)力分析

    為研究熱-機(jī)械載荷對(duì)缸蓋應(yīng)力的影響,利用有限元軟件計(jì)算了以下幾種工況:1)預(yù)裝配載荷(F),主要考慮氣門導(dǎo)管過盈裝配,加載缸蓋螺栓預(yù)緊力(130kN)和噴油器壓塊處螺栓預(yù)緊力(6 250N);2)預(yù)裝配載荷+熱載荷(F+T);3)預(yù)裝配載荷+最高燃燒壓力載荷(F+P);4)預(yù)裝配載荷+熱載荷+最高燃燒壓力載荷(F+T+P)。

    為驗(yàn)證仿真分析模型,在預(yù)裝配載荷條件下對(duì)缸蓋頂部和靠近螺栓側(cè)壁等區(qū)域進(jìn)行了應(yīng)力測(cè)試。與仿真值對(duì)比見圖11,由圖可見,多數(shù)測(cè)點(diǎn)應(yīng)力偏差較小,有限元仿真模型合理可信。

    圖11 預(yù)裝配載荷下應(yīng)力仿真與實(shí)測(cè)對(duì)比

    通過有限元應(yīng)力計(jì)算得到了預(yù)裝配載荷、熱載荷和氣體力載荷單獨(dú)作用及耦合作用下缸蓋應(yīng)力場(chǎng),結(jié)合溫度、應(yīng)力及疲勞計(jì)算結(jié)果,從缸蓋頂板、水腔、進(jìn)排氣道壁、底板冷卻水鉆孔和底板分別選取若干個(gè)特征點(diǎn)或危險(xiǎn)點(diǎn)作為應(yīng)力考察點(diǎn)(見圖12)。圖13示出各考察點(diǎn)應(yīng)力對(duì)比情況,圖14示出預(yù)裝配載荷、熱載荷、氣體力載荷對(duì)3種耦合作用下Mises應(yīng)力的貢獻(xiàn)系數(shù)。由圖13和圖14可知,預(yù)裝配載荷F對(duì)頂板、水腔和進(jìn)排氣道區(qū)域多數(shù)考察點(diǎn)的應(yīng)力貢獻(xiàn)較大,并占主導(dǎo)地位;熱載荷T作用下,由于底板溫度較高,底板向燃燒室方向膨脹,對(duì)頂部考察點(diǎn)耦合應(yīng)力有釋放緩解作用,而對(duì)底板冷卻鉆孔產(chǎn)生了較高的拉應(yīng)力,對(duì)鼻梁區(qū)產(chǎn)生了較高的壓應(yīng)力,熱載荷對(duì)底板考察點(diǎn)耦合應(yīng)力占主導(dǎo)作用。

    圖12 缸蓋應(yīng)力分析考察點(diǎn)

    圖13 不同載荷下考察點(diǎn)應(yīng)力

    圖14 不同載荷對(duì)耦合應(yīng)力的影響

    對(duì)于缸蓋結(jié)構(gòu)高周疲勞來說,預(yù)裝配載荷和熱載荷主要影響平均應(yīng)力,而氣體力載荷主要影響應(yīng)力幅。圖15示出考慮材料彈塑性影響修正后各考察點(diǎn)的應(yīng)力疲勞極限圖,由圖可見,載荷呈明顯非對(duì)稱分布,其中底板冷卻水鉆孔處部分考察點(diǎn)負(fù)載較高。

    圖15 各考察點(diǎn)應(yīng)力疲勞極限圖

    3.2 溫度波動(dòng)對(duì)熱機(jī)耦合應(yīng)力的影響

    在預(yù)裝配載荷和穩(wěn)態(tài)定常熱載荷的基礎(chǔ)上,考慮氣體力、高頻波動(dòng)熱載荷對(duì)缸蓋耦合應(yīng)力的影響,通過計(jì)算得到了其隨曲軸轉(zhuǎn)角的變化關(guān)系。圖16示出工作循環(huán)內(nèi)瞬態(tài)溫度波動(dòng)引起的高頻熱應(yīng)力的最大波動(dòng)幅值,由圖可知,應(yīng)力波動(dòng)區(qū)域與溫度瞬態(tài)變化直接相關(guān),高頻波動(dòng)熱應(yīng)力僅存在于火力面薄層區(qū)域內(nèi)[6]。取兩個(gè)排氣門之間鼻梁區(qū)D3考察點(diǎn)為研究對(duì)象,沿徑向正應(yīng)力的變化情況見圖17。由圖可見,氣體力作用導(dǎo)致壓應(yīng)力幅值進(jìn)一步增加,高頻波動(dòng)熱應(yīng)力相位與氣體力應(yīng)力明顯不同且數(shù)值較大,綜合應(yīng)力波動(dòng)幅為19.4MPa,是只考慮氣體力作用情況下應(yīng)力的1.9倍,可見,高頻波動(dòng)熱應(yīng)力對(duì)火力面區(qū)域結(jié)構(gòu)表面高周疲勞影響較大。

    圖16 熱應(yīng)力波動(dòng)幅

    圖17 考察點(diǎn)D3應(yīng)力波動(dòng)曲線

    4 變工況載荷特性分析

    按照發(fā)動(dòng)機(jī)熱沖擊規(guī)范,發(fā)動(dòng)機(jī)在圖1所示怠速和標(biāo)定轉(zhuǎn)速工況交替運(yùn)行,各工況運(yùn)行時(shí)間為3min,缸蓋溫度和應(yīng)力發(fā)生大幅度變化。在前文穩(wěn)態(tài)溫度場(chǎng)計(jì)算的基礎(chǔ)上,進(jìn)一步考慮材料溫度非線性特性,計(jì)算了變工況下缸蓋溫度和應(yīng)力隨時(shí)間的變化歷程(見圖18)。由圖18可知,在加熱過程中缸蓋底板特征點(diǎn)溫度急速增加,至60s后趨于穩(wěn)定,在加熱過程前60s內(nèi)平均加熱速率約為7.1 ℃/s,在冷卻過程前60s內(nèi)平均冷卻速率約為6.6 ℃/s;鼻梁區(qū)觀察點(diǎn)D2處應(yīng)力與溫度變化基本呈現(xiàn)反相狀態(tài),即最高溫度對(duì)應(yīng)最大壓應(yīng)力,而遠(yuǎn)離火力面的其他考察點(diǎn)溫度較低,應(yīng)力也較小。由此可見,除火力面鼻梁之外的其他區(qū)域溫度水平很低,最大應(yīng)力遠(yuǎn)小于材料屈服強(qiáng)度,而火力面鼻梁區(qū)結(jié)構(gòu)負(fù)載最高,為低周疲勞重點(diǎn)關(guān)注部位。

    圖18 變工況下溫度和應(yīng)力隨時(shí)間的變化曲線

    5 疲勞損傷與壽命預(yù)測(cè)模型建立

    在圖1所示試驗(yàn)工況下,缸蓋部件產(chǎn)生的總損傷可以表示為

    (8)

    式中:DTotal為材料發(fā)生的總損傷;DHCF為高周疲勞試驗(yàn)考核階段產(chǎn)生的總損傷;DLCF為低周疲勞試驗(yàn)考核階段所產(chǎn)生的總損傷。

    基于前文分析可知,式(8)中高周疲勞試驗(yàn)考核階段產(chǎn)生的總損傷DHCF主要由高頻波動(dòng)的氣體力和溫度載荷引起,可以描述為

    (9)

    式中:DF,HCF為高周高溫機(jī)械疲勞損傷;DC,HCF為蠕變損傷;kt為工作循環(huán)內(nèi)溫度高頻波動(dòng)熱應(yīng)力影響系數(shù),與缸蓋空間位置有關(guān),遠(yuǎn)離火力面溫度瞬態(tài)波動(dòng)層區(qū)域該系數(shù)為1,火力面區(qū)域可通過計(jì)算確定。

    在低周疲勞試驗(yàn)考核階段,缸蓋主要承受高低周疲勞載荷復(fù)合作用,材料產(chǎn)生的總損傷可以認(rèn)為是由于溫度大幅波動(dòng)而產(chǎn)生的低周熱機(jī)疲勞(TMF)損傷、高頻氣體力及熱載荷波動(dòng)引起高周疲勞損傷的疊加。其中,根據(jù)Sehitoglu理論[7],考慮疲勞與蠕變、疲勞與氧化之間的交互作用,熱機(jī)疲勞損傷DTMF主要包括疲勞損傷DF,TMF、氧化損傷DO,TMF和蠕變損傷DC,TMF三部分。因此,該階段產(chǎn)生的總損傷和低周熱機(jī)疲勞損傷可分別描述如下:

    (10)

    DTMF=DF,TMF+DO,TMF+DC,TMF。

    (11)

    由式(8)至式(11)可得,試驗(yàn)工況下缸蓋部件產(chǎn)生總損傷為

    (12)

    依照線性損傷法則,在該服役工況下缸蓋構(gòu)件的疲勞壽命方程[7]描述如下:

    (13)

    6 結(jié)論

    a) 在穩(wěn)態(tài)轉(zhuǎn)速工況下,在缸蓋火力面約1mm的薄層內(nèi)溫度存在明顯的瞬態(tài)波動(dòng),而其他區(qū)域溫度基本不變;工作循環(huán)內(nèi)瞬態(tài)溫度波動(dòng)引起的高頻熱應(yīng)力與溫度瞬態(tài)變化直接相關(guān),高頻波動(dòng)熱應(yīng)力僅存在于火力面薄層區(qū),氣體力作用導(dǎo)致該區(qū)域壓應(yīng)力進(jìn)一步增加,故高周疲勞強(qiáng)度評(píng)估時(shí)必須考慮高頻熱應(yīng)力的影響;

    b) 在怠速和標(biāo)定轉(zhuǎn)速交替變工況運(yùn)行下,缸蓋底板區(qū)域溫度發(fā)生大幅度變化,在加熱過程中缸蓋底板特征點(diǎn)溫度急速增加至60s后趨于穩(wěn)定,在加熱過程前60s內(nèi)平均加熱速率約為7.1 ℃/s,在冷卻過程前60s內(nèi)平均冷卻速率約為6.6 ℃/s;鼻梁區(qū)應(yīng)力與溫度變化基本呈現(xiàn)反相狀態(tài),而遠(yuǎn)離火力面的其他考察點(diǎn)溫度和應(yīng)力水平較低,鼻梁區(qū)為低周疲勞重點(diǎn)關(guān)注部位;

    c) 基于對(duì)缸蓋溫度和應(yīng)力變化特征的研究,提出了基于臺(tái)架服役工況的缸蓋疲勞損傷壽命預(yù)測(cè)模型,為后續(xù)壽命預(yù)測(cè)工作奠定了基礎(chǔ)。

    [1]ThomasHamm,MartinRebbert,Hermann-JosefEcker,etal.CylinderHeadDesignforHighPeakFiringPressures[C].SAEPaper2008-01-1196.

    [2]HazratMdAli,MasjukiHassan,KalamMd,etal.DeterminationofTimeVariant1D-3DTemperatureandHeatTransferDistri-butionInsidetheCoolingJacketofaSIEngineCoolingSystemafterKey-Off[C].SAEPaper2012-01-0642.

    [3]MarioMetzger,MarcLeidenfrost,EwaldWerner,etal.LifetimePredictionofEN-GJV450CastIronCylinderHeads[C].SAEPaper2014-01-9047.

    [4]HalfpennyA,AndersonR,LinX.IsothermalandThermo-MechanicalFatigueofAutomotiveComponents[C].SAEPaper2015-01-0548.

    [5]MetzgerM,SeifertT.Amechanism-basedmodelforLCF/HCFandTMF/HCFlifeprediction:Multiaxialformulation,finite-elementimplementationandapplicationtocastiron[J].TechnischeMechanik,2012,32:435-445.

    [6] 原彥鵬, 張衛(wèi)正, 程曉果,等.高強(qiáng)化內(nèi)燃機(jī)活塞瞬態(tài)溫度場(chǎng)分布規(guī)律研究[J].內(nèi)燃機(jī)工程,2005,26(4):35-38.

    [7]SehitogluH.,BoismierDA.Thermo-mechanicalfatigueofMar-M247:Part2Lifeprediction[J].TRANS.ASMEJ.Engng.MaterTechnol.,1990,112:80-89.

    [編輯: 李建新]

    Loading Characteristics of Cast Iron Cylinder Head for High-strengthened Diesel Engine

    JING Guoxi1, ZHANG Shuyong1, FU Wenqing2, XU Chunguang1, WEN Yang1, LI Peng1, WEI Junchao1

    (1. China North Engine Research Institute(Tianjin), Tianjin 300400, China;2. Ordnance Science and Research Academy of China, Beijing 100089, China)

    The changing behavior of temperature and stress for a cast iron cylinder head was investigated. Then the loading characteristics under thermal and mechanical loads were determined, which can provide the reference for building the life prediction model of cylinder head. The head finite element model was first built and calibrated with temperature and stress data and was further used for predicting the temperature and stress varying characteristics under steady working condition and altering condition between idle and rated speed. The simulation results show that the maximum temperature and thermal stress locate near the exhaust valve bridge region on the firing deck. The amplitude of temperature fluctuation in the firing deck region is less than 30℃, but it will still lead to thermal stress. The produced thermal stress is negligible compared with the steady thermal stress, but important to have to be considered in high cycle fatigue evaluation compared with high frequency gas stress. The load varies obviously in firing deck region under the altering condition between idle and rated speed. The temperature and stress are in out of phase, the stress amplitude is very high, and therefore the firing deck region generates the low cycle fatigue easily.

    diesel engine; cylinder head; simulation; loading characteristic; temperature; stress

    2016-01-26;

    2016-11-02

    景國(guó)璽(1984—),男,副研究員,博士,研究方向?yàn)閮?nèi)燃機(jī)現(xiàn)代設(shè)計(jì)方法、疲勞強(qiáng)度設(shè)計(jì)等;okjgx@163.com。

    10.3969/j.issn.1001-2222.2017.01.008

    TK413.2

    B

    1001-2222(2017)01-0042-07

    猜你喜歡
    缸蓋熱應(yīng)力火力
    WNS型鍋爐煙管管端熱應(yīng)力裂紋原因分析
    火力全開
    火力全開! 廣州上半年20條村改造,投入超800億!
    模擬缸蓋工藝在當(dāng)代內(nèi)燃機(jī)制造業(yè)中的應(yīng)用
    高性能連體缸蓋關(guān)鍵鑄造工藝研究及應(yīng)用
    《火力與指揮控制》投稿須知
    采用單元基光滑點(diǎn)插值法的高溫管道熱應(yīng)力分析
    缸蓋氣門導(dǎo)管座圈的加工及改進(jìn)
    鑄造缺陷對(duì)A356 鋁合金缸蓋疲勞性能的影響
    上海金屬(2015年5期)2015-11-29 01:13:43
    基于流熱固耦合的核電蒸汽發(fā)生器傳熱管熱應(yīng)力數(shù)值模擬
    在线观看午夜福利视频| 国产亚洲精品久久久久久毛片| 日韩欧美国产一区二区入口| 88av欧美| www.www免费av| 又爽又黄无遮挡网站| 亚洲av中文字字幕乱码综合| 亚洲中文字幕日韩| 国产又色又爽无遮挡免费看| 母亲3免费完整高清在线观看| 又爽又黄无遮挡网站| 黑人操中国人逼视频| 国产一区二区三区视频了| 精品国产美女av久久久久小说| 成人午夜高清在线视频| 久久亚洲真实| 亚洲熟妇熟女久久| 老司机午夜福利在线观看视频| a在线观看视频网站| a级毛片在线看网站| 后天国语完整版免费观看| 欧美丝袜亚洲另类 | 长腿黑丝高跟| 视频区欧美日本亚洲| av有码第一页| 听说在线观看完整版免费高清| 久久人妻福利社区极品人妻图片| 国模一区二区三区四区视频 | www.www免费av| 欧美人与性动交α欧美精品济南到| 亚洲国产欧美一区二区综合| 欧美黑人巨大hd| 最近最新中文字幕大全免费视频| 少妇裸体淫交视频免费看高清 | 少妇裸体淫交视频免费看高清 | 1024手机看黄色片| 99精品久久久久人妻精品| 久久久久久亚洲精品国产蜜桃av| 免费在线观看成人毛片| 国产日本99.免费观看| 操出白浆在线播放| 嫁个100分男人电影在线观看| 无人区码免费观看不卡| 欧美乱妇无乱码| 在线观看免费午夜福利视频| 欧美国产日韩亚洲一区| 日日爽夜夜爽网站| 欧美日韩国产亚洲二区| 老司机在亚洲福利影院| 妹子高潮喷水视频| 亚洲自拍偷在线| 亚洲成人国产一区在线观看| 一个人免费在线观看电影 | 亚洲五月婷婷丁香| 又大又爽又粗| 欧美人与性动交α欧美精品济南到| 可以在线观看的亚洲视频| 久久久国产成人精品二区| 高清在线国产一区| 亚洲av中文字字幕乱码综合| 午夜日韩欧美国产| 天堂av国产一区二区熟女人妻 | 男女午夜视频在线观看| 国产91精品成人一区二区三区| or卡值多少钱| 欧美成人性av电影在线观看| 国产精华一区二区三区| 国产日本99.免费观看| 91大片在线观看| 一级毛片精品| 少妇熟女aⅴ在线视频| 999久久久国产精品视频| 亚洲精品在线观看二区| 成人三级黄色视频| 色综合亚洲欧美另类图片| 露出奶头的视频| 18禁美女被吸乳视频| 1024视频免费在线观看| 这个男人来自地球电影免费观看| 十八禁人妻一区二区| 黄色片一级片一级黄色片| 欧美另类亚洲清纯唯美| 91老司机精品| 天天一区二区日本电影三级| 精品国产亚洲在线| 97碰自拍视频| 国产精品综合久久久久久久免费| 巨乳人妻的诱惑在线观看| 免费看a级黄色片| 日韩成人在线观看一区二区三区| 亚洲18禁久久av| 精品免费久久久久久久清纯| 国产av又大| 欧美国产日韩亚洲一区| 久久精品国产综合久久久| 动漫黄色视频在线观看| 老司机午夜十八禁免费视频| 极品教师在线免费播放| 免费看十八禁软件| 99久久精品国产亚洲精品| 美女扒开内裤让男人捅视频| 99热这里只有精品一区 | 国产精品精品国产色婷婷| 一个人观看的视频www高清免费观看 | 91九色精品人成在线观看| 亚洲精品中文字幕在线视频| 国产av一区二区精品久久| 色播亚洲综合网| 久久精品91无色码中文字幕| 日本一区二区免费在线视频| 欧美极品一区二区三区四区| 窝窝影院91人妻| 国产亚洲精品一区二区www| 亚洲av五月六月丁香网| 我的老师免费观看完整版| 黄色 视频免费看| 久久欧美精品欧美久久欧美| 久久久久久人人人人人| 丁香欧美五月| 丰满人妻一区二区三区视频av | 国产精品久久电影中文字幕| 一级毛片精品| 欧美精品亚洲一区二区| 精品电影一区二区在线| 在线观看一区二区三区| 天堂√8在线中文| 69av精品久久久久久| 欧美日韩福利视频一区二区| 日本三级黄在线观看| 欧美乱色亚洲激情| 亚洲 国产 在线| 波多野结衣高清作品| 黄色片一级片一级黄色片| 熟妇人妻久久中文字幕3abv| 99热只有精品国产| 最新在线观看一区二区三区| 欧美三级亚洲精品| 一夜夜www| 久久久久精品国产欧美久久久| 国产99白浆流出| 亚洲av日韩精品久久久久久密| 性欧美人与动物交配| 国产精品1区2区在线观看.| 99精品欧美一区二区三区四区| 午夜激情福利司机影院| 免费无遮挡裸体视频| 88av欧美| 97超级碰碰碰精品色视频在线观看| 国产精品国产高清国产av| 91麻豆精品激情在线观看国产| 亚洲色图 男人天堂 中文字幕| 看黄色毛片网站| 丁香欧美五月| 精品熟女少妇八av免费久了| 国内精品一区二区在线观看| 国产黄a三级三级三级人| 女人被狂操c到高潮| 变态另类成人亚洲欧美熟女| 亚洲狠狠婷婷综合久久图片| 久久人妻av系列| 日韩欧美国产一区二区入口| 午夜福利在线观看吧| 窝窝影院91人妻| 欧美av亚洲av综合av国产av| 欧美乱妇无乱码| 91九色精品人成在线观看| 久久亚洲精品不卡| 免费在线观看成人毛片| 亚洲精品一卡2卡三卡4卡5卡| 国产高清视频在线观看网站| 亚洲av电影不卡..在线观看| 男人舔女人的私密视频| 男女视频在线观看网站免费 | 我要搜黄色片| 观看免费一级毛片| 国产真人三级小视频在线观看| 一本久久中文字幕| 欧美性长视频在线观看| 久久香蕉国产精品| 级片在线观看| 亚洲国产精品999在线| 午夜福利在线观看吧| 亚洲av中文字字幕乱码综合| 亚洲成人国产一区在线观看| 婷婷丁香在线五月| 国产精品久久久久久亚洲av鲁大| 国产亚洲av嫩草精品影院| 一区二区三区国产精品乱码| 91麻豆精品激情在线观看国产| 精品无人区乱码1区二区| 手机成人av网站| 波多野结衣高清作品| 少妇的丰满在线观看| 亚洲国产高清在线一区二区三| 精品久久久久久久毛片微露脸| www.精华液| 99久久久亚洲精品蜜臀av| 男女之事视频高清在线观看| 丁香六月欧美| 18禁国产床啪视频网站| 欧美一级毛片孕妇| 亚洲av成人不卡在线观看播放网| 男人的好看免费观看在线视频 | 久久精品91无色码中文字幕| 真人做人爱边吃奶动态| 在线永久观看黄色视频| 又爽又黄无遮挡网站| 亚洲精品在线美女| 免费搜索国产男女视频| 真人做人爱边吃奶动态| 五月玫瑰六月丁香| 国产精品永久免费网站| 男人舔奶头视频| 精品久久久久久久末码| 又黄又粗又硬又大视频| www.999成人在线观看| 久久久国产成人精品二区| 精品一区二区三区视频在线观看免费| cao死你这个sao货| 一本综合久久免费| 亚洲国产欧美人成| 久久热在线av| 亚洲片人在线观看| 久久久国产成人精品二区| 亚洲一区二区三区色噜噜| 俄罗斯特黄特色一大片| 真人做人爱边吃奶动态| 波多野结衣高清作品| 人妻久久中文字幕网| 国产爱豆传媒在线观看 | 亚洲专区中文字幕在线| 香蕉国产在线看| 亚洲最大成人中文| 亚洲国产精品合色在线| 久久中文字幕一级| 久久久久久久久中文| 久久99热这里只有精品18| 午夜激情福利司机影院| 99久久久亚洲精品蜜臀av| 夜夜爽天天搞| 99国产精品一区二区三区| АⅤ资源中文在线天堂| 两个人看的免费小视频| 亚洲美女视频黄频| 精品午夜福利视频在线观看一区| 男人舔女人的私密视频| 久久久久九九精品影院| 国产精品一区二区三区四区免费观看 | 国产午夜福利久久久久久| 狠狠狠狠99中文字幕| 日日摸夜夜添夜夜添小说| 两个人免费观看高清视频| bbb黄色大片| 精品少妇一区二区三区视频日本电影| а√天堂www在线а√下载| 一级毛片精品| 国产欧美日韩精品亚洲av| 国产激情久久老熟女| 欧美又色又爽又黄视频| 欧美日韩福利视频一区二区| 老汉色∧v一级毛片| 夜夜爽天天搞| 又黄又爽又免费观看的视频| 欧美日韩亚洲综合一区二区三区_| 国产精品一区二区免费欧美| 一a级毛片在线观看| 亚洲人成伊人成综合网2020| 欧美一区二区精品小视频在线| 1024视频免费在线观看| 老司机在亚洲福利影院| 大型av网站在线播放| 亚洲午夜理论影院| videosex国产| 免费一级毛片在线播放高清视频| 天堂动漫精品| 日本一本二区三区精品| 一进一出好大好爽视频| 高清在线国产一区| 久久精品国产亚洲av香蕉五月| 欧美3d第一页| 亚洲熟女毛片儿| 亚洲精品粉嫩美女一区| 又紧又爽又黄一区二区| 1024香蕉在线观看| 操出白浆在线播放| 亚洲天堂国产精品一区在线| 999久久久精品免费观看国产| 亚洲,欧美精品.| 好男人在线观看高清免费视频| 中文资源天堂在线| 一区福利在线观看| 99久久无色码亚洲精品果冻| www.精华液| 中文亚洲av片在线观看爽| 一进一出好大好爽视频| 成人手机av| 脱女人内裤的视频| 又粗又爽又猛毛片免费看| 亚洲精品国产精品久久久不卡| 国产欧美日韩一区二区精品| 久久 成人 亚洲| 亚洲美女黄片视频| 国产精品久久久久久人妻精品电影| 午夜久久久久精精品| 少妇粗大呻吟视频| 久久人妻福利社区极品人妻图片| 精品久久久久久久久久久久久| 又粗又爽又猛毛片免费看| 国产成人精品无人区| 级片在线观看| 伊人久久大香线蕉亚洲五| 51午夜福利影视在线观看| 国产aⅴ精品一区二区三区波| 一级a爱片免费观看的视频| 热99re8久久精品国产| 久久婷婷成人综合色麻豆| 日韩大码丰满熟妇| 午夜两性在线视频| 在线观看午夜福利视频| 欧美成人午夜精品| 免费搜索国产男女视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲七黄色美女视频| 欧美乱码精品一区二区三区| 亚洲九九香蕉| 国产午夜精品论理片| 一个人免费在线观看的高清视频| 99久久99久久久精品蜜桃| 琪琪午夜伦伦电影理论片6080| 中文字幕久久专区| 麻豆国产97在线/欧美 | 日本免费一区二区三区高清不卡| 免费在线观看黄色视频的| 国产一级毛片七仙女欲春2| www.www免费av| 黄色成人免费大全| 午夜福利在线观看吧| 白带黄色成豆腐渣| 中亚洲国语对白在线视频| 欧美激情久久久久久爽电影| 国产成人精品无人区| 日本熟妇午夜| 免费搜索国产男女视频| av福利片在线观看| 久久精品成人免费网站| 国产视频一区二区在线看| 日本 欧美在线| 成年版毛片免费区| 久久精品成人免费网站| 亚洲电影在线观看av| 免费看十八禁软件| 亚洲 欧美 日韩 在线 免费| 亚洲五月婷婷丁香| 黄色 视频免费看| 真人一进一出gif抽搐免费| 国产精品久久久久久久电影 | 国产成人精品无人区| 亚洲最大成人中文| 男女视频在线观看网站免费 | 国产真实乱freesex| 黄色 视频免费看| 1024手机看黄色片| 亚洲 国产 在线| 国产成人av教育| 日韩精品青青久久久久久| 可以在线观看的亚洲视频| 午夜福利高清视频| 国产日本99.免费观看| 婷婷精品国产亚洲av在线| 午夜激情福利司机影院| 午夜福利18| 精品无人区乱码1区二区| 神马国产精品三级电影在线观看 | 日韩精品免费视频一区二区三区| 日本在线视频免费播放| 精品乱码久久久久久99久播| 久久国产精品人妻蜜桃| 一二三四社区在线视频社区8| 制服人妻中文乱码| 岛国视频午夜一区免费看| 夜夜躁狠狠躁天天躁| 老汉色av国产亚洲站长工具| 麻豆一二三区av精品| 亚洲av日韩精品久久久久久密| 男男h啪啪无遮挡| 亚洲av电影在线进入| 欧美日本视频| 男女之事视频高清在线观看| 国产黄片美女视频| 嫁个100分男人电影在线观看| 亚洲电影在线观看av| 好男人在线观看高清免费视频| 国产高清视频在线播放一区| 亚洲人与动物交配视频| 2021天堂中文幕一二区在线观| 久久久国产成人精品二区| 18禁黄网站禁片免费观看直播| 精品欧美一区二区三区在线| 国产精品av视频在线免费观看| 国产高清videossex| 国产精品久久久久久精品电影| 久久午夜综合久久蜜桃| 亚洲,欧美精品.| 啦啦啦观看免费观看视频高清| 国产精品九九99| 99国产精品99久久久久| √禁漫天堂资源中文www| 99久久精品热视频| 人妻丰满熟妇av一区二区三区| 床上黄色一级片| 国产99久久九九免费精品| 中文字幕人妻丝袜一区二区| 男男h啪啪无遮挡| 国产精品av久久久久免费| 亚洲va日本ⅴa欧美va伊人久久| 国产成人精品无人区| 亚洲欧美精品综合一区二区三区| 制服诱惑二区| 99久久久亚洲精品蜜臀av| 国产成人av教育| 久久久国产欧美日韩av| 午夜福利成人在线免费观看| 国产私拍福利视频在线观看| 黄色丝袜av网址大全| 在线免费观看的www视频| 夜夜躁狠狠躁天天躁| 国产精品日韩av在线免费观看| 99久久99久久久精品蜜桃| 日本免费a在线| 欧美日韩精品网址| 制服诱惑二区| 日本五十路高清| 老司机午夜福利在线观看视频| 国产主播在线观看一区二区| 91老司机精品| 舔av片在线| 超碰成人久久| 亚洲精华国产精华精| 久久久久久免费高清国产稀缺| 日本一区二区免费在线视频| 成人三级黄色视频| 正在播放国产对白刺激| 国产一区二区激情短视频| 香蕉久久夜色| 亚洲五月天丁香| 亚洲一区二区三区色噜噜| 国产成人一区二区三区免费视频网站| 五月伊人婷婷丁香| 亚洲天堂国产精品一区在线| 国产精品日韩av在线免费观看| 久久精品亚洲精品国产色婷小说| 大型av网站在线播放| 亚洲欧美日韩东京热| 国产精品 国内视频| 夜夜看夜夜爽夜夜摸| 欧美色视频一区免费| 欧美又色又爽又黄视频| 99热这里只有精品一区 | 听说在线观看完整版免费高清| 又黄又爽又免费观看的视频| 国产亚洲av嫩草精品影院| 嫩草影视91久久| 国产成人欧美在线观看| 国产精品久久久久久精品电影| xxx96com| 久久久久久久久久黄片| 成人18禁在线播放| 一区福利在线观看| 亚洲人成伊人成综合网2020| 嫩草影视91久久| 很黄的视频免费| 又粗又爽又猛毛片免费看| 可以在线观看的亚洲视频| 12—13女人毛片做爰片一| 午夜免费激情av| 欧美激情久久久久久爽电影| 国产高清有码在线观看视频 | 97人妻精品一区二区三区麻豆| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲成人免费电影在线观看| 中亚洲国语对白在线视频| 成年版毛片免费区| 亚洲av第一区精品v没综合| 好看av亚洲va欧美ⅴa在| 久久久久久久午夜电影| 怎么达到女性高潮| 五月玫瑰六月丁香| 精品久久久久久久久久免费视频| 久久婷婷人人爽人人干人人爱| 亚洲av电影不卡..在线观看| 国产一区在线观看成人免费| 又爽又黄无遮挡网站| 国产av在哪里看| 午夜激情福利司机影院| 欧美一区二区精品小视频在线| 丝袜人妻中文字幕| 日本成人三级电影网站| 最近最新中文字幕大全免费视频| 国产亚洲精品一区二区www| 成在线人永久免费视频| 午夜日韩欧美国产| 波多野结衣高清作品| 九色国产91popny在线| 一进一出好大好爽视频| 国产单亲对白刺激| 日韩欧美精品v在线| 久久久久久亚洲精品国产蜜桃av| 大型黄色视频在线免费观看| 后天国语完整版免费观看| 精品一区二区三区视频在线观看免费| 午夜福利在线观看吧| 在线观看免费视频日本深夜| 国产成人av教育| 久久亚洲精品不卡| 麻豆一二三区av精品| 99久久无色码亚洲精品果冻| 国内少妇人妻偷人精品xxx网站 | 日韩精品中文字幕看吧| 麻豆成人午夜福利视频| 美女午夜性视频免费| 欧美一级a爱片免费观看看 | 亚洲中文字幕日韩| 久久久久久亚洲精品国产蜜桃av| 18禁国产床啪视频网站| 好看av亚洲va欧美ⅴa在| 日韩高清综合在线| 婷婷精品国产亚洲av在线| 高清在线国产一区| 国产精品自产拍在线观看55亚洲| 一边摸一边做爽爽视频免费| 两人在一起打扑克的视频| 国产免费av片在线观看野外av| 中文字幕av在线有码专区| 亚洲黑人精品在线| 我的老师免费观看完整版| 欧美成人免费av一区二区三区| 国产精品免费一区二区三区在线| 国产熟女xx| 日本一二三区视频观看| 日韩成人在线观看一区二区三区| 国产成人av激情在线播放| 精品免费久久久久久久清纯| 黄色 视频免费看| 色av中文字幕| 国产精品日韩av在线免费观看| 午夜激情av网站| 亚洲人与动物交配视频| 国产精品久久久久久精品电影| 91在线观看av| 久久性视频一级片| 搡老岳熟女国产| 国产成人av激情在线播放| 日本a在线网址| 99久久久亚洲精品蜜臀av| 熟女少妇亚洲综合色aaa.| 床上黄色一级片| 欧美精品亚洲一区二区| 很黄的视频免费| 美女黄网站色视频| 国产免费男女视频| 国产欧美日韩一区二区精品| 国产精品香港三级国产av潘金莲| 国产激情偷乱视频一区二区| 久久亚洲真实| av天堂在线播放| 精品第一国产精品| 男女之事视频高清在线观看| 日韩欧美 国产精品| 国产成人啪精品午夜网站| or卡值多少钱| 国产精品久久久久久久电影 | 在线看三级毛片| 久久中文看片网| 一本大道久久a久久精品| 最近在线观看免费完整版| 亚洲精品久久成人aⅴ小说| 亚洲av片天天在线观看| 男插女下体视频免费在线播放| 人人妻人人澡欧美一区二区| 成人三级做爰电影| 在线免费观看的www视频| 精品一区二区三区av网在线观看| 麻豆国产97在线/欧美 | 深夜精品福利| 欧美久久黑人一区二区| 国产黄a三级三级三级人| 国内久久婷婷六月综合欲色啪| 国产免费av片在线观看野外av| 淫妇啪啪啪对白视频| 一个人免费在线观看的高清视频| 欧美日本视频| 日本a在线网址| 精品少妇一区二区三区视频日本电影| 成年人黄色毛片网站| 久久精品人妻少妇| 午夜精品在线福利| www.自偷自拍.com| 久久久久国产一级毛片高清牌| 免费av毛片视频| 日日爽夜夜爽网站| 国产精品av久久久久免费| 欧美乱色亚洲激情| 99re在线观看精品视频| 人妻久久中文字幕网| 久久亚洲精品不卡| 欧美久久黑人一区二区| 久久精品综合一区二区三区| 亚洲专区字幕在线| 夜夜爽天天搞| 亚洲精品在线美女| 久久久精品国产亚洲av高清涩受| 亚洲国产精品999在线| 午夜福利成人在线免费观看| 在线永久观看黄色视频| 观看免费一级毛片|