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

    核電汽輪機轉(zhuǎn)子日歷總壽命的計算方法

    2023-09-26 04:22:48史進淵謝岳生江路毅李汪繁范雪飛
    動力工程學(xué)報 2023年9期
    關(guān)鍵詞:汽輪機裂紋

    史進淵, 謝岳生, 江路毅, 李汪繁, 范雪飛

    (上海發(fā)電設(shè)備成套設(shè)計研究院有限責(zé)任公司,上海 200240)

    核能發(fā)電沒有二氧化碳排放,是實現(xiàn)“碳中和”的重要發(fā)電方式之一。單臺核電機組功率可達1 000~1 755 MW,相當(dāng)于500~877.5臺2 MW風(fēng)電機組的裝機容量;核電機組的利用小時數(shù)長達7 400~8 600 h,相當(dāng)于風(fēng)電機組利用小時數(shù)(約2 000 h)的4倍。發(fā)展核電是替代化石能源、減少能源領(lǐng)域碳排放的重要技術(shù)途徑。太陽能、風(fēng)能等可再生能源發(fā)電存在間歇性、隨機性和波動性等問題,而核電機組具有大容量和長利用小時數(shù)的特征,在新型電力系統(tǒng)中提供穩(wěn)定的基本負(fù)荷,可以為電力供應(yīng)起到兜底保障的“頂梁柱”和“壓艙石”的作用。

    壓水堆核電汽輪機的進汽為飽和蒸汽,其特點是進汽參數(shù)較低、進汽流量大,結(jié)構(gòu)尺寸大、半速多缸設(shè)計,中間再熱去濕、防水蝕設(shè)計[1]。3缸4排汽核電汽輪機有3根轉(zhuǎn)子,4缸6排汽核電汽輪機有4根轉(zhuǎn)子。1 000 MW等級核電汽輪機低壓轉(zhuǎn)子重300多噸,約為相同功率火電汽輪機低壓轉(zhuǎn)子質(zhì)量(75 t)的4倍。超大尺寸核電汽輪機轉(zhuǎn)子是確定核電汽輪機壽命的關(guān)鍵部件,核電汽輪機的使用壽命與轉(zhuǎn)子壽命密切相關(guān)。在核電汽輪機起動、停機與超速試驗等瞬態(tài)過程中轉(zhuǎn)子面臨低周疲勞損傷問題,而在核電汽輪機帶負(fù)荷穩(wěn)態(tài)運行工況下轉(zhuǎn)子面臨高周疲勞損傷問題。

    核電汽輪機轉(zhuǎn)子低周疲勞壽命的單位是疲勞循環(huán)次數(shù),核電汽輪機轉(zhuǎn)子高周疲勞壽命的單位是與運行小時數(shù)有關(guān)的疲勞循環(huán)次數(shù)。電站業(yè)主要求的第三代核電汽輪機的設(shè)計壽命是60 a[2],工程上急需以使用年數(shù)表示的核電汽輪機轉(zhuǎn)子日歷總壽命。文獻[3]中介紹了核電汽輪機轉(zhuǎn)子在高周疲勞與低周疲勞損傷機理交互作用下裂紋擴展壽命的計算方法,但研究交互作用下裂紋萌生壽命與日歷總壽命(CL)的文獻還比較少。因此,筆者研究了此交互作用下核電汽輪機轉(zhuǎn)子裂紋萌生日歷壽命與日歷總壽命的計算方法,有助于保障核電汽輪機的安全服役。

    1 轉(zhuǎn)子高周疲勞裂紋萌生壽命

    1.1 轉(zhuǎn)子高周疲勞損傷機理

    核電汽輪機帶負(fù)荷穩(wěn)態(tài)運行工況下,受離心力載荷與熱載荷的作用,核電汽輪機轉(zhuǎn)子中將產(chǎn)生穩(wěn)定的平均應(yīng)力。核電汽輪機轉(zhuǎn)子質(zhì)量為90~300 t,并采用大跨距軸承支撐。在動葉片與轉(zhuǎn)子自重的作用下,轉(zhuǎn)子重心以上部位將產(chǎn)生壓應(yīng)力,轉(zhuǎn)子重心以下部位將產(chǎn)生拉應(yīng)力。

    核電汽輪機帶負(fù)荷穩(wěn)態(tài)運行工況下,轉(zhuǎn)子外表面某一點除了承受因離心力載荷與熱載荷的作用產(chǎn)生的穩(wěn)定的平均應(yīng)力以外,當(dāng)該點旋轉(zhuǎn)到轉(zhuǎn)子上部90°位置時,在動葉片與轉(zhuǎn)子自重單獨作用下該點將產(chǎn)生壓應(yīng)力;當(dāng)該點旋轉(zhuǎn)到轉(zhuǎn)子下部270°位置時,在動葉片與轉(zhuǎn)子自重單獨作用下該點將產(chǎn)生拉應(yīng)力。

    核電汽輪機轉(zhuǎn)子外表面某一點旋轉(zhuǎn)到不同位置時的交變應(yīng)力是由動葉片與轉(zhuǎn)子自重引起的,每轉(zhuǎn)一周,核電汽輪機轉(zhuǎn)子產(chǎn)生一次疲勞循環(huán)。對于半轉(zhuǎn)速核電汽輪機,每秒疲勞循環(huán)25次,若每年運行7 000 h,則60 a高周疲勞循環(huán)3.78×1010次。

    核電汽輪機在帶負(fù)荷穩(wěn)態(tài)運行工況下每轉(zhuǎn)一周,轉(zhuǎn)子高周疲勞損傷一次,轉(zhuǎn)子在平均應(yīng)力σmH與高周疲勞應(yīng)力幅σaH的作用下也會產(chǎn)生裂紋萌生和裂紋擴展。鑒于核電汽輪機轉(zhuǎn)子高周疲勞壽命薄弱部位的平均應(yīng)力σmH遠(yuǎn)遠(yuǎn)大于高周疲勞應(yīng)力幅σaH,雖然核電汽輪機每轉(zhuǎn)一周時轉(zhuǎn)子的高周疲勞壽命損耗非常小,但是在核電汽輪機60 a服役期內(nèi),轉(zhuǎn)子高周疲勞循環(huán)次數(shù)的數(shù)量級非常大,因此仍需要考慮轉(zhuǎn)子累計高周疲勞壽命損傷引起的裂紋萌生和裂紋擴展。

    1.2 高周疲勞裂紋萌生壽命計算方法

    1.2.1 高周疲勞裂紋萌生的力學(xué)量

    根據(jù)文獻[3]和文獻[4],在核電汽輪機帶額定負(fù)荷穩(wěn)態(tài)運行工況下,計算核電汽輪機轉(zhuǎn)子某一部位高周疲勞的應(yīng)力幅σaH、平均應(yīng)力σmH、應(yīng)力范圍ΔσH、應(yīng)力比RH:

    (1)

    (2)

    ΔσH=σmaxH-σminH

    (3)

    (4)

    式中:σmaxH、σminH分別為核電汽輪機轉(zhuǎn)子某一部位高周疲勞的最大應(yīng)力、最小應(yīng)力。

    1.2.2 高周疲勞應(yīng)變幅

    在核電汽輪機帶負(fù)荷穩(wěn)態(tài)運行工況下,轉(zhuǎn)子的高周疲勞應(yīng)變幅εaH為:

    (5)

    式中:E為材料的彈性模量;μ為材料的泊松比。

    1.2.3 核電汽輪機轉(zhuǎn)子高周疲勞裂紋萌生壽命

    根據(jù)文獻[5],核電汽輪機轉(zhuǎn)子高周疲勞裂紋萌生壽命NH的計算公式為:

    (6)

    式中:σf為材料的疲勞強度系數(shù);εf為材料的疲勞延性系數(shù);bH為材料的高周疲勞強度指數(shù);c為材料的疲勞延性指數(shù)。

    2 轉(zhuǎn)子低周疲勞裂紋萌生壽命

    2.1 轉(zhuǎn)子低周疲勞損傷機理

    在核電汽輪機起動、停機與超速試驗等瞬態(tài)過程中,由于核電汽輪機超大尺寸轉(zhuǎn)子外表面與轉(zhuǎn)子中心存在較大的溫度差,轉(zhuǎn)子中會產(chǎn)生較大的熱應(yīng)力。熱應(yīng)力與由轉(zhuǎn)子離心力產(chǎn)生的機械應(yīng)力的合成應(yīng)力有可能為拉應(yīng)力,也有可能為壓應(yīng)力。核電汽輪機轉(zhuǎn)子的拉應(yīng)力與壓應(yīng)力均會引起轉(zhuǎn)子裂紋萌生,但壓應(yīng)力不會引起轉(zhuǎn)子裂紋疲勞擴展,只有拉應(yīng)力才會引起轉(zhuǎn)子裂紋疲勞擴展。

    通常在核電汽輪機起動過程中,轉(zhuǎn)子外表面溫度高,轉(zhuǎn)子內(nèi)部與內(nèi)表面溫度低,轉(zhuǎn)子外表面材料受壓,轉(zhuǎn)子內(nèi)部與內(nèi)表面材料受拉。在核電汽輪機起動過程中,轉(zhuǎn)子外表面的熱應(yīng)力為壓應(yīng)力,由熱應(yīng)力與轉(zhuǎn)子離心力產(chǎn)生的機械應(yīng)力合成的綜合應(yīng)力較小或有可能為壓應(yīng)力;轉(zhuǎn)子內(nèi)部與內(nèi)表面的熱應(yīng)力為拉應(yīng)力,由熱應(yīng)力與轉(zhuǎn)子離心力產(chǎn)生的機械應(yīng)力合成的綜合應(yīng)力為拉應(yīng)力。

    在核電汽輪機的停機過程中,轉(zhuǎn)子外表面溫度低,轉(zhuǎn)子內(nèi)部與內(nèi)表面溫度高,轉(zhuǎn)子外表面材料受拉,轉(zhuǎn)子內(nèi)部與內(nèi)表面材料受壓。轉(zhuǎn)子外表面的熱應(yīng)力為拉應(yīng)力,由熱應(yīng)力與轉(zhuǎn)子離心力產(chǎn)生的機械應(yīng)力合成的綜合應(yīng)力為拉應(yīng)力;轉(zhuǎn)子內(nèi)部與內(nèi)表面的熱應(yīng)力有可能為壓應(yīng)力,由熱應(yīng)力與轉(zhuǎn)子離心力產(chǎn)生的機械應(yīng)力合成的綜合應(yīng)力較小或有可能為壓應(yīng)力。

    核電汽輪機起動一次與停機一次,核電汽輪機轉(zhuǎn)子疲勞循環(huán)一次。在核電汽輪機60 a服役期內(nèi),核電汽輪機的起動次數(shù)不超過10 000次,核電汽輪機的起動與停機疲勞循環(huán)屬于低周疲勞。核電汽輪機快速起動時,轉(zhuǎn)子熱應(yīng)力有可能遠(yuǎn)超材料的屈服極限,轉(zhuǎn)子低周疲勞損傷較大,核電汽輪機起停一次,轉(zhuǎn)子的壽命損傷也較大。核電汽輪機轉(zhuǎn)子在起動與停機等瞬態(tài)過程的低周疲勞損傷機理作用下發(fā)生裂紋萌生與裂紋擴展。

    2.2 低周疲勞裂紋萌生壽命計算方法

    在核電汽輪機起動與停機過程中,拉應(yīng)力與壓應(yīng)力均會引起核電汽輪機轉(zhuǎn)子裂紋萌生。由于核電汽輪機起動與停機過程產(chǎn)生的轉(zhuǎn)子熱應(yīng)力并不相同,核電汽輪機轉(zhuǎn)子的低周疲勞循環(huán)屬于非對稱低周疲勞循環(huán)。核電汽輪機在冷態(tài)起動、溫態(tài)起動與熱態(tài)起動過程中的熱應(yīng)力也不相同,由一次冷態(tài)起動與一次正常停機構(gòu)成的一個低周疲勞循環(huán)、一次溫態(tài)起動與一次正常停機構(gòu)成的一個低周疲勞循環(huán)以及一次熱態(tài)起動與一次正常停機構(gòu)成的一個低周疲勞循環(huán),這三種低周疲勞循環(huán)的應(yīng)變幅并不相同,可見核電汽輪機轉(zhuǎn)子的低周疲勞屬于變幅載荷的低周疲勞循環(huán)。

    (1) 確定帶有正負(fù)號的等效應(yīng)變。區(qū)別于核電汽輪機轉(zhuǎn)子的高周疲勞循環(huán),在核電汽輪機轉(zhuǎn)子低周疲勞裂紋萌生壽命的計算中,低周疲勞循環(huán)的最大應(yīng)變εmax、最小應(yīng)變εmin、應(yīng)力幅σa、應(yīng)變幅εa和應(yīng)力比R等力學(xué)量不加下標(biāo)“L”。核電汽輪機轉(zhuǎn)子的低周疲勞裂紋萌生壽命計算中采用的低周疲勞循環(huán)的最大應(yīng)變εmax和最小應(yīng)變εmin均為等效應(yīng)變。在進行核電汽輪機轉(zhuǎn)子有限元數(shù)值計算時按照馮·米塞斯公式得出的等效應(yīng)變?yōu)檎?而在構(gòu)建核電汽輪機轉(zhuǎn)子非對稱低周疲勞循環(huán)裂紋萌生壽命計算力學(xué)模型時,先采用有限元分析方法計算轉(zhuǎn)子等效應(yīng)變εeq,再按照文獻[6]和文獻[7]中提出的方法確定帶有正負(fù)號的等效應(yīng)變εei:

    sign(ε1+ε3)=εeqsign(ε1+ε3)

    (7)

    式中:ε1、ε2、ε3為主應(yīng)變;sign(x)為符號函數(shù),若x>0,則sign(x)=1,若x<0,則sign(x)=-1。

    按照式(7)確定等效應(yīng)變的正負(fù)號,其物理意義表示等效應(yīng)變的正負(fù)號與3個主應(yīng)變ε1、ε2、ε3絕對值最大者的正負(fù)號相同。采用帶有正負(fù)號的等效應(yīng)變εei來計算核電汽輪機轉(zhuǎn)子的低周疲勞循環(huán)應(yīng)變幅εa,相當(dāng)于認(rèn)為最大應(yīng)變εmax與最小應(yīng)變εmin的方向相同;事實上,最大應(yīng)變εmax與最小應(yīng)變εmin的方向并不相同,采用上述處理為方向相同的方法來預(yù)測核電汽輪機轉(zhuǎn)子的低周疲勞裂紋萌生壽命,經(jīng)超臨界汽輪機高壓內(nèi)缸開裂事故驗證,這種計算方法計算結(jié)果偏于安全[3,7],且計算精度高[6]。采用sign(ε1+ε3)識別等效應(yīng)變的正負(fù)號,物理意義明確,數(shù)學(xué)公式簡單,工程上使用方便。

    (2) 計算最大應(yīng)變和最小應(yīng)變。在核電汽輪機一次起動與一次停機的低周疲勞循環(huán)中,轉(zhuǎn)子低周疲勞循環(huán)的最大應(yīng)變εmax和最小應(yīng)變εmin分別為:

    εmax=max{εe1,εe2,…,εei,…,εem}

    (8)

    εmin=min{εe1,εe2,…,εei,…,εem}

    (9)

    式中:εe1,εe2,…,εei,…,εem為一次低周疲勞循環(huán)中不同時刻的瞬態(tài)等效應(yīng)變。

    核電汽輪機轉(zhuǎn)子非對稱低周疲勞循環(huán)模型中分以下3種情況計算應(yīng)變幅εa:

    ① 若εmax>0且εmin<0,則

    (10)

    ② 若εmax>0且εmin>0,則

    (11)

    ③ 若εmax<0且εmin<0,則

    (12)

    式中:Δε為疲勞循環(huán)的應(yīng)變范圍或全應(yīng)變。

    對于εmax和εmin均大于0的情況,式(11)中低周疲勞循環(huán)最小應(yīng)變εmin取核電汽輪機停機后轉(zhuǎn)子處于靜止?fàn)顟B(tài)的0。對于εmax和εmin均小于0的情況,式(12)中低周疲勞循環(huán)最大應(yīng)變εmax取核電汽輪機停機后轉(zhuǎn)子處于靜止?fàn)顟B(tài)的0。

    (3) 計算應(yīng)力比和非對稱循環(huán)的疲勞極限??紤]到核電汽輪機轉(zhuǎn)子材料的不可壓縮與各向同性,核電汽輪機轉(zhuǎn)子非對稱低周疲勞循環(huán)的理論應(yīng)力幅σa的計算公式為:

    (13)

    核電汽輪機轉(zhuǎn)子低周疲勞循環(huán)的最大實際應(yīng)力σmax為:

    (14)

    核電汽輪機轉(zhuǎn)子低周疲勞循環(huán)的應(yīng)力比R為:

    (15)

    核電汽輪機轉(zhuǎn)子材料非對稱循環(huán)的疲勞極限σN為:

    (16)

    式中:σ-1、σb分別為在低周疲勞裂紋萌生壽命曲線試驗溫度θ下材料的疲勞極限和抗拉強度。

    (4) 計算考慮轉(zhuǎn)子材料抗拉強度因素的修正值

    (17)

    式中:σb(θ1)、σb(θ2)分別為最大應(yīng)變對應(yīng)工作溫度θ1、最小應(yīng)變對應(yīng)工作溫度θ2下材料的抗拉強度。

    (5) 確定材料缺口敏感性修正系數(shù)Cq:

    (18)

    Kf=1+q(Kt-1)

    (19)

    式中:Kf為有效應(yīng)力集中系數(shù);Kt為理論應(yīng)力集中系數(shù);q為缺口敏感系數(shù)。

    (6) 計算變幅載荷非對稱低周疲勞循環(huán)模型的等效應(yīng)變幅。根據(jù)文獻[6]~文獻[8],考慮了轉(zhuǎn)子應(yīng)力比、低周疲勞循環(huán)非對稱性、材料缺口敏感性、疲勞極限和抗拉強度等影響因素并進行修正后,核電汽輪機一次起動與一次停機的轉(zhuǎn)子變幅載荷下非對稱低周疲勞循環(huán)模型中低周疲勞應(yīng)變幅εai為:

    (20)

    (7) 計算轉(zhuǎn)子低周疲勞裂紋萌生壽命Ni。依據(jù)核電汽輪機第i種起動與停機工況下轉(zhuǎn)子低周疲勞應(yīng)變幅εai和材料疲勞力學(xué)性能數(shù)據(jù),計算此工況下轉(zhuǎn)子低周疲勞裂紋萌生壽命Ni。核電汽輪機第i種工況下轉(zhuǎn)子低周疲勞應(yīng)變幅εai包括核電汽輪機轉(zhuǎn)子的冷態(tài)起停低周疲勞應(yīng)變幅εac、溫態(tài)起停低周疲勞應(yīng)變幅εaw、熱態(tài)起停低周疲勞應(yīng)變幅εah、110%超速試驗低周疲勞應(yīng)變幅εa110和120%超速運行低周疲勞應(yīng)變幅εa120。核電汽輪機第i種工況下轉(zhuǎn)子低周疲勞裂紋萌生壽命Ni包括核電汽輪機轉(zhuǎn)子冷態(tài)起停低周疲勞裂紋萌生壽命Nc、溫態(tài)起停低周疲勞裂紋萌生壽命Nw、熱態(tài)起停低周疲勞裂紋萌生壽命Nh、110%超速試驗低周疲勞裂紋萌生壽命N110和120%超速運行低周疲勞裂紋萌生壽命N120。核電汽輪機第i種工況下轉(zhuǎn)子低周疲勞裂紋萌生壽命Ni為:

    (21)

    式中:bL為材料低周疲勞強度指數(shù)。

    3 轉(zhuǎn)子高周疲勞裂紋擴展壽命

    3.1 高周疲勞裂紋擴展材料試驗常數(shù)

    (22)

    式中:ΔKI為應(yīng)力強度因子的范圍。

    對于核電汽輪機的轉(zhuǎn)子材料,不同高周疲勞的應(yīng)力比RH對應(yīng)于不同的材料試驗常數(shù)C0H(RH)和m0H(RH)。通過試驗可以確定核電汽輪機轉(zhuǎn)子材料的各種應(yīng)力比RH對應(yīng)的材料試驗常數(shù)C0H(RH)和m0H(RH),但試驗費用較高。根據(jù)文獻[7],將C0H(RH)和m0H(RH)表示為RH的函數(shù)的4個經(jīng)驗公式分別如下:

    (23)

    (24)

    (25)

    (26)

    式中:C0(0)和m0(0)均為RH=0時疲勞裂紋擴展Paris公式的材料試驗常數(shù);n為待定材料試驗常數(shù)。

    對于核電汽輪機轉(zhuǎn)子的CrNiMoV材料,由試驗得出的C0H(RH)和m0H(RH)見表1。對于表1數(shù)據(jù),采用非線性回歸分析方法,按式(23)和式(24)擬合得出的材料試驗常數(shù)C0(0)、n與回歸分析的相關(guān)指數(shù)r2見表2,相關(guān)指數(shù)r2的計算方法詳見參考文獻[11]。按式(23)和式(24)得出的擬合曲線分別如圖1和圖2所示。按式(25)和式(26)擬合得出的材料試驗常數(shù)m0(0)、n與回歸分析的相關(guān)指數(shù)r2見表3,擬合曲線分別如圖3和圖4所示。

    圖1 CrNiMoV轉(zhuǎn)子材料C0H(RH)按式(23)的變化曲線

    圖2 CrNiMoV轉(zhuǎn)子材料C0H(RH)按式(24)的變化曲線

    圖3 CrNiMoV轉(zhuǎn)子材料m0H(RH)按式(25)的變化曲線

    圖4 CrNiMoV轉(zhuǎn)子材料m0H(RH)按式(26)的變化曲線

    表1 CrNiMoV轉(zhuǎn)子材料C0H(RH)和m0H(RH)的試驗結(jié)果

    表2 由式(23)和式(24)得到的材料試驗常數(shù)C0(0)與n的回歸分析結(jié)果

    表3 由式(25)和式(26)得到的材料試驗常數(shù)m0(0)與n的回歸分析結(jié)果

    從表2和表3可知,采用式(24)和式(26)得到的相關(guān)指數(shù)r2比采用式(23)和式(25)得到的相關(guān)指數(shù)r2大,表明采用式(24)和式(26)擬合表1中試驗數(shù)據(jù)的精度更高。通常核電汽輪機轉(zhuǎn)子高周疲勞的應(yīng)力比RH不小于0.9。從圖1~圖4可知,在核電汽輪機轉(zhuǎn)子壽命設(shè)計中使用的高周疲勞的應(yīng)力比RH大于0.9時,這4個公式擬合試驗數(shù)據(jù)的精度均較高,可以用來確定核電汽輪機轉(zhuǎn)子材料高周疲勞裂紋擴展材料試驗常數(shù)。

    3.2 高周疲勞裂紋擴展壽命計算方法

    3.2.1 轉(zhuǎn)子高周疲勞裂紋擴展壽命計算的力學(xué)量

    根據(jù)文獻[3]和文獻[4],在核電汽輪機帶額定負(fù)荷穩(wěn)態(tài)運行工況下,計算核電汽輪機轉(zhuǎn)子某一部位高周疲勞的最大應(yīng)力σmaxH、最小應(yīng)力σminH和應(yīng)力范圍ΔσH:

    σmaxH=σν,270°+σre

    (27)

    σminH=σν,90°+σre

    (28)

    式中:σν,270°為轉(zhuǎn)子某一部位旋轉(zhuǎn)到轉(zhuǎn)子下部270°位置時的最大主應(yīng)力,其中ν為最大主應(yīng)力的方向;σν,90°為轉(zhuǎn)子的某一點旋轉(zhuǎn)到轉(zhuǎn)子上部90°位置時在ν方向的正應(yīng)力;σre為轉(zhuǎn)子由于鍛造、焊接與熱處理等制造工藝產(chǎn)生的殘余應(yīng)力。

    3.2.2 轉(zhuǎn)子高周疲勞裂紋擴展的裂紋尺寸界限值

    核電汽輪機轉(zhuǎn)子高周疲勞裂紋擴展的裂紋尺寸界限值ath為:

    (29)

    ath的物理意義是核電汽輪機轉(zhuǎn)子在高周疲勞循環(huán)的應(yīng)力幅σaH作用下不發(fā)生裂紋疲勞擴展的門檻值,表示當(dāng)核電汽輪機轉(zhuǎn)子的初始裂紋尺寸不大于ath時,在高周疲勞循環(huán)的應(yīng)力幅σaH作用下轉(zhuǎn)子的初始裂紋不會發(fā)生高周疲勞擴展。

    3.2.3 轉(zhuǎn)子高周疲勞臨界裂紋尺寸

    核電汽輪機轉(zhuǎn)子高周疲勞臨界裂紋尺寸acH為:

    (30)

    式中:KIC為轉(zhuǎn)子材料的斷裂韌性。

    3.2.4 轉(zhuǎn)子高周疲勞裂紋擴展壽命

    核電汽輪機轉(zhuǎn)子裂紋從高周疲勞裂紋擴展的裂紋尺寸界限值ath至高周疲勞臨界裂紋尺寸acH的過程中高周疲勞裂紋擴展壽命NfH為:

    (31)

    4 轉(zhuǎn)子低周疲勞裂紋擴展壽命

    4.1 低周疲勞裂紋擴展材料試驗常數(shù)

    在核電汽輪機轉(zhuǎn)子低周疲勞裂紋擴展壽命的計算中,低周疲勞的最大應(yīng)力、最小應(yīng)力、應(yīng)力比和應(yīng)力范圍等力學(xué)量分別用σmaxL、σminL、RL和ΔσL表示。在核電汽輪機停機后轉(zhuǎn)子處于靜止?fàn)顟B(tài)的應(yīng)力可以視為0 MPa,即低周疲勞的最小應(yīng)力σminL=0 MPa。由于RL=0,核電汽輪機轉(zhuǎn)子低周疲勞裂紋擴展的低周疲勞循環(huán)為脈動循環(huán)。核電汽輪機轉(zhuǎn)子低周疲勞裂紋擴展的材料試驗常數(shù)是式(22)的C0和m0。

    (32)

    (33)

    4.2 低周疲勞裂紋擴展壽命計算方法

    4.2.1 轉(zhuǎn)子低周疲勞裂紋擴展壽命計算的力學(xué)量

    根據(jù)文獻[4]和文獻[7],考慮到核電汽輪機停機時轉(zhuǎn)子的應(yīng)力為0 MPa,即轉(zhuǎn)子低周疲勞裂紋擴展的最小應(yīng)力σminL=0 MPa,在核電汽輪機起動、停機、超速試驗與超速運行等瞬態(tài)過程中,計算核電汽輪機轉(zhuǎn)子低周疲勞裂紋擴展的最大應(yīng)力σmaxL、應(yīng)力比RL和應(yīng)力范圍ΔσL:

    σmaxL=σ1max+σre

    (34)

    (35)

    ΔσL=σmaxL-σminL=σmaxL

    (36)

    式中:σ1max為低周疲勞裂紋擴展最大主應(yīng)力。

    4.2.2 轉(zhuǎn)子低周疲勞臨界裂紋尺寸

    核電汽輪機轉(zhuǎn)子低周疲勞臨界裂紋尺寸acL為:

    (37)

    4.2.3 第一階段低周疲勞裂紋擴展壽命

    在低周疲勞損傷機理作用下,核電汽輪機轉(zhuǎn)子裂紋從初始裂紋尺寸ai至高周疲勞裂紋擴展的裂紋尺寸界限值ath的第一階段低周疲勞裂紋擴展壽命Nfi,1為:

    (38)

    根據(jù)文獻[7],在探傷沒有發(fā)現(xiàn)缺陷的情況下,假定初始裂紋尺寸ai=2 mm,計算低周疲勞裂紋擴展壽命。

    4.2.4 第二階段低周疲勞裂紋擴展壽命

    在低周疲勞損傷機理作用下,核電汽輪機轉(zhuǎn)子裂紋從高周疲勞裂紋擴展的裂紋尺寸界限值ath至低周疲勞臨界裂紋尺寸acL的第二階段低周疲勞裂紋擴展壽命Nfi,2為:

    (39)

    5 轉(zhuǎn)子日歷壽命

    5.1 日歷壽命的物理量

    核電汽輪機轉(zhuǎn)子的日歷壽命計算中需要考慮汽輪機帶負(fù)荷穩(wěn)態(tài)運行工況的高周疲勞損傷和起停工況的低周疲勞損傷。在核電汽輪機帶負(fù)荷穩(wěn)態(tài)運行工況下轉(zhuǎn)子產(chǎn)生高周疲勞損傷;在核電汽輪機起動、停機與超速試驗等瞬態(tài)過程中轉(zhuǎn)子產(chǎn)生低周疲勞損傷。在核電汽輪機的服役期內(nèi),其帶負(fù)荷穩(wěn)態(tài)運行工況與起動、停機、超速試驗等瞬態(tài)工況交替出現(xiàn),則核電汽輪機轉(zhuǎn)子在高周疲勞與低周疲勞損傷機理交互作用下發(fā)生裂紋萌生與裂紋擴展。

    核電汽輪機部件日歷壽命是以日歷年數(shù)為單位表示的核電汽輪機關(guān)鍵部件的設(shè)計壽命??紤]到核電汽輪機轉(zhuǎn)子尺寸大且高速旋轉(zhuǎn),是核電汽輪機的關(guān)鍵部件,其日歷壽命可以用來表征核電汽輪機的日歷壽命。核電汽輪機轉(zhuǎn)子的裂紋萌生日歷壽命和裂紋擴展日歷壽命也可以采用日歷年數(shù)表示,兩者之和稱為日歷總壽命。

    在高周疲勞與低周疲勞損傷機理交互作用下,核電汽輪機轉(zhuǎn)子的日歷壽命等于其年均壽命損耗的倒數(shù)。年均壽命損耗也稱年均壽命損傷,其物理意義為在高周疲勞與低周疲勞損傷機理交互作用下核電汽輪機轉(zhuǎn)子每年的平均壽命損耗。依據(jù)核電汽輪機年均起動次數(shù)、年均停機次數(shù)、年均超速試驗次數(shù)等,以及核電汽輪機轉(zhuǎn)子的高周疲勞裂紋萌生壽命、低周疲勞裂紋萌生壽命、高周疲勞裂紋擴展壽命、低周疲勞裂紋擴展壽命等,可以確定核電汽輪機轉(zhuǎn)子的年均壽命損耗與日歷壽命。

    在計算核電汽輪機轉(zhuǎn)子日歷壽命的過程中,采用以下2項假設(shè):第一項假設(shè)是核電汽輪機冷態(tài)起動、穩(wěn)態(tài)起動、熱態(tài)起動、正常停機、110%超速試驗、120%超速運行等瞬態(tài)過程中與帶負(fù)荷穩(wěn)態(tài)運行工況下轉(zhuǎn)子的等效應(yīng)力或拉應(yīng)力的方向相同,實際上這些應(yīng)力的方向并不相同,在壽命計算中假定方向相同,計算結(jié)果偏于安全[7];第二項假設(shè)是在高周疲勞與低周疲勞損傷機理交互作用下,核電汽輪機轉(zhuǎn)子裂紋萌生或裂紋擴展的壽命損耗可以使用Palmgern-Miner壽命損耗線性疊加法則[4,7,16-17],這是工程上常用的一種簡化處理方法。

    5.2 裂紋萌生日歷壽命的計算方法

    拉應(yīng)力(正應(yīng)變)與壓應(yīng)力(負(fù)應(yīng)變)均會引起核電汽輪機轉(zhuǎn)子裂紋萌生,同樣高周疲勞與低周疲勞均會導(dǎo)致核電汽輪機轉(zhuǎn)子裂紋萌生。在高周疲勞與低周疲勞損傷機理交互作用下,核電汽輪機轉(zhuǎn)子部件裂紋萌生年均壽命損耗ey0和裂紋萌生日歷壽命τCL0分別為:

    (40)

    (41)

    式中:yc為核電汽輪機年均冷態(tài)起動次數(shù);yw為核電汽輪機年均溫態(tài)起動次數(shù);yh為核電汽輪機年均熱態(tài)起動次數(shù);y110為核電汽輪機年均110%超速試驗次數(shù);y120為核電汽輪機年均120%超速運行次數(shù);yH為核電汽輪機轉(zhuǎn)子年均高周疲勞次數(shù);NH為高周疲勞裂紋萌生壽命。

    5.3 第一階段裂紋擴展日歷壽命的計算方法

    在核電汽輪機轉(zhuǎn)子裂紋擴展的第一階段,由于初始裂紋尺寸ai小于轉(zhuǎn)子高周疲勞裂紋擴展的裂紋尺寸界限值ath,轉(zhuǎn)子只在低周疲勞損傷機理作用下發(fā)生裂紋擴展。由于壓應(yīng)力不會引起疲勞裂紋擴展,核電汽輪機轉(zhuǎn)子產(chǎn)生疲勞裂紋擴展的條件之一是承受比較大的拉應(yīng)力。在核電汽輪機轉(zhuǎn)子的外表面,在停機等瞬態(tài)過程中會出現(xiàn)較大的拉應(yīng)力;在核電汽輪機轉(zhuǎn)子內(nèi)部與內(nèi)表面,在起動等瞬態(tài)過程中會出現(xiàn)較大的拉應(yīng)力。

    (1) 轉(zhuǎn)子外表面第一階段裂紋擴展日歷壽命。對于核電汽輪機轉(zhuǎn)子外表面,在核電汽輪機正常停機、110%超速試驗和120%超速運行等瞬態(tài)過程中會出現(xiàn)較大的拉應(yīng)力,從而引起轉(zhuǎn)子低周疲勞裂紋擴展。核電汽輪機轉(zhuǎn)子外表面的第一階段裂紋擴展年均壽命損耗ey1o和日歷壽命τCL1o分別為:

    (42)

    (43)

    式中:yn為核電汽輪機年均正常停機次數(shù);Nfn,1為核電汽輪機正常停機過程中轉(zhuǎn)子第一階段低周疲勞裂紋擴展壽命;Nf110,1為核電汽輪機110%超速試驗過程中轉(zhuǎn)子第一階段低周疲勞裂紋擴展壽命;Nf120,1為核電汽輪機120%超速運行過程中轉(zhuǎn)子第一階段低周疲勞裂紋擴展壽命。

    (2) 轉(zhuǎn)子內(nèi)部與內(nèi)表面第一階段裂紋擴展日歷壽命。對于核電汽輪機轉(zhuǎn)子內(nèi)部與內(nèi)表面,在核電汽輪機冷態(tài)起動、溫態(tài)起動、熱態(tài)起動、110%超速試驗和120%超速運行等瞬態(tài)過程中會出現(xiàn)較大的拉應(yīng)力,從而引起轉(zhuǎn)子低周疲勞裂紋擴展。核電汽輪機轉(zhuǎn)子內(nèi)部與內(nèi)表面的第一階段裂紋擴展年均壽命損耗ey1i和日歷壽命τCL1i分別為:

    (44)

    (45)

    式中:Nfc,1為核電汽輪機冷態(tài)起動過程中轉(zhuǎn)子第一階段低周疲勞裂紋擴展壽命;Nfw,1為核電汽輪機溫態(tài)起動過程中轉(zhuǎn)子第一階段低周疲勞裂紋擴展壽命;Nfh,1為核電汽輪機熱態(tài)起動過程中轉(zhuǎn)子第一階段低周疲勞裂紋擴展壽命。

    5.4 第二階段裂紋擴展日歷壽命的計算方法

    在核電汽輪機轉(zhuǎn)子裂紋擴展的第二階段,由于初始裂紋尺寸ai大于轉(zhuǎn)子高周疲勞裂紋擴展的裂紋尺寸界限值ath,轉(zhuǎn)子在低周疲勞與高周疲勞損傷機理交互作用下發(fā)生裂紋擴展。

    (1) 轉(zhuǎn)子外表面第二階段裂紋擴展日歷壽命。核電汽輪機轉(zhuǎn)子外表面的第二階段裂紋擴展年均壽命損耗ey2o和日歷壽命τCL2o分別為:

    (46)

    (47)

    式中:Nfn,2為核電汽輪機正常停機過程中轉(zhuǎn)子第二階段低周疲勞裂紋擴展壽命;Nf110,2為核電汽輪機110%超速試驗過程中轉(zhuǎn)子第二階段低周疲勞裂紋擴展壽命;Nf120,2為核電汽輪機120%超速運行過程中轉(zhuǎn)子第二階段低周疲勞裂紋擴展壽命;NfH為轉(zhuǎn)子高周疲勞裂紋擴展壽命。

    (2) 轉(zhuǎn)子內(nèi)部與內(nèi)表面第二階段裂紋擴展日歷壽命。核電汽輪機轉(zhuǎn)子內(nèi)部與內(nèi)表面的第二階段裂紋擴展年均壽命損耗ey2i和日歷壽命τCL2i分別為:

    (48)

    (49)

    式中:Nfc,2為核電汽輪機冷態(tài)起動過程中轉(zhuǎn)子第二階段低周疲勞裂紋擴展壽命;Nfw,2為核電汽輪機溫態(tài)起動過程中轉(zhuǎn)子第二階段低周疲勞裂紋擴展壽命;Nfh,2為核電汽輪機熱態(tài)起動過程中轉(zhuǎn)子第二階段低周疲勞裂紋擴展壽命。

    6 轉(zhuǎn)子日歷總壽命

    6.1 日歷總壽命的計算方法

    核電汽輪機轉(zhuǎn)子從首次投入使用到其裂紋尺寸擴展達到臨界裂紋尺寸為止,轉(zhuǎn)子所經(jīng)歷的時間或循環(huán)次數(shù)稱為總壽命,所經(jīng)歷的日歷年數(shù)稱為日歷總壽命。

    在高周疲勞與低周疲勞損傷機理交互作用下,核電汽輪機轉(zhuǎn)子外表面日歷總壽命τCLto、內(nèi)部與內(nèi)表面日歷總壽命τCLti及轉(zhuǎn)子日歷總壽命τCLr分別為:

    τCLto=τCL0+τCL1o+τCL2o

    (50)

    τCLti=τCL0+τCL1i+τCL2i

    (51)

    τCLr=min{τCLto,τCLti}

    (52)

    6.2 日歷總壽命的評價方法

    依據(jù)核電站業(yè)主的要求,第三代壓水堆核電汽輪機日歷總壽命的判據(jù)值為使用壽命τ0=60 a。根據(jù)核電汽輪機轉(zhuǎn)子日歷總壽命的計算結(jié)果和核電汽輪機日歷總壽命的判據(jù)值,按照以下判據(jù)對核電汽輪機轉(zhuǎn)子壽命設(shè)計結(jié)果進行評價:(1) 若τCLr≥τ0,在高周疲勞與低周疲勞損傷機理交互作用下,核電汽輪機轉(zhuǎn)子日歷總壽命設(shè)計合格;(2) 若τCLr<τ0,在高周疲勞與低周疲勞損傷機理交互作用下,核電汽輪機轉(zhuǎn)子日歷總壽命設(shè)計不合格,表明在設(shè)計階段需要對核電汽輪機轉(zhuǎn)子的結(jié)構(gòu)進行優(yōu)化改進,直到τCLr≥τ0為止。

    7 應(yīng)用實例

    以某1 000 MW等級半轉(zhuǎn)速核電汽輪機為研究對象,焊接低壓轉(zhuǎn)子材料為CrNiMoV鋼,低壓轉(zhuǎn)子有5條焊縫,并采用對稱結(jié)構(gòu)。核電汽輪機焊接低壓轉(zhuǎn)子的示意圖如圖5所示。通常核電汽輪機的起停工況只有冷態(tài)起動、溫態(tài)起動、熱態(tài)起動和正常停機4種工況。表4給出了電站業(yè)主要求的核電汽輪機的壽命指標(biāo)值[18]。

    圖5 核電汽輪機焊接低壓轉(zhuǎn)子的示意圖

    表4 核電汽輪機的壽命指標(biāo)值

    在完成該核電汽輪機焊接低壓轉(zhuǎn)子溫度場與應(yīng)力場有限元計算分析后,確定焊低壓轉(zhuǎn)子壽命薄弱部位是圖5中的低壓轉(zhuǎn)子母材外表面部位A。

    采用變幅載荷與非對稱循環(huán)模型計算該焊接低壓轉(zhuǎn)子母材外表面部位A的高周疲勞裂紋萌生壽命與低周疲勞裂紋萌生壽命,計算結(jié)果列于表5。

    表5 焊接低壓轉(zhuǎn)子裂紋萌生壽命的計算結(jié)果

    該核電汽輪機焊接低壓轉(zhuǎn)子母材外表面部位A的高周疲勞裂紋擴展壽命和低周疲勞裂紋擴展壽命的計算結(jié)果分別列于表6和表7。

    表6 焊接低壓轉(zhuǎn)子高周疲勞裂紋擴展壽命的計算結(jié)果

    表7 焊接低壓轉(zhuǎn)子低周疲勞裂紋擴展壽命的計算結(jié)果

    該核電汽輪機焊接低壓轉(zhuǎn)子母材外表面部位A的年均壽命損耗與日歷壽命的計算結(jié)果列于表8。從表8可知,該焊接低壓轉(zhuǎn)子母材外表面部位A的疲勞裂紋擴展日歷總壽命超過了40 a,但還未達到核電站業(yè)主要求的60 a。由于τCLr<τ0,該核電汽輪機焊接低壓轉(zhuǎn)子原設(shè)計結(jié)構(gòu)在高周疲勞與低周疲勞損傷機理交互作用下裂紋萌生與裂紋擴展日歷總壽命監(jiān)控不合格,表明在設(shè)計階段需要對該核電汽輪機焊接低壓轉(zhuǎn)子的結(jié)構(gòu)進行改進。在該轉(zhuǎn)子結(jié)構(gòu)改進設(shè)計中,采用增大焊接低壓轉(zhuǎn)子表面部位A的圓角半徑的結(jié)構(gòu)改進措施,將該焊接低壓轉(zhuǎn)子部位A處的圓弧半徑增大50%,則該部位在核電汽輪機正常停機過程中的最大應(yīng)力σmaxL可從458.027 MPa減少到423.311 MPa。采取結(jié)構(gòu)改進措施后,該部位的日歷壽命計算結(jié)果列于表9。從表9可知,該核電汽輪機焊接低壓轉(zhuǎn)子結(jié)構(gòu)改進后τCLr>τ0,表明該核電汽輪機焊接低壓轉(zhuǎn)子在高周疲勞與低周疲勞損傷機理交互作用下裂紋萌生與裂紋擴展日歷總壽命設(shè)計合格。

    表8 原結(jié)構(gòu)焊接低壓轉(zhuǎn)子日歷壽命計算結(jié)果

    表9 改進后焊接低壓轉(zhuǎn)子日歷壽命計算結(jié)果

    8 分析與討論

    文中建立的核電汽輪機轉(zhuǎn)子日歷壽命計算方法中不僅考慮了核電汽輪機在帶負(fù)荷穩(wěn)態(tài)運行工況下轉(zhuǎn)子的高周疲勞壽命損傷,還考慮了核電汽輪機冷態(tài)起動、溫態(tài)起動、熱態(tài)起動、正常停機、110%超速試驗、120%超速運行等瞬態(tài)過程中轉(zhuǎn)子的低周疲勞壽命損傷。

    核電汽輪機轉(zhuǎn)子的裂紋萌生日歷壽命計算方法中采用一次起動過程與一次停機過程構(gòu)成的非對稱低周疲勞循環(huán)與變幅載荷的計算模型,在考慮了低周疲勞裂紋萌生的應(yīng)力比、非對稱疲勞循環(huán)疲勞極限、抗拉強度與應(yīng)力集中影響的同時,還考慮了高周疲勞與低周疲勞損傷機理的交互作用。

    核電汽輪機轉(zhuǎn)子的裂紋擴展日歷壽命計算方法中采用冷態(tài)起動、溫態(tài)起動、熱態(tài)起動、正常停機、110%超速試驗、120%超速運行等變幅載荷模型來計算轉(zhuǎn)子低周疲勞裂紋擴展壽命損傷,另外還考慮了核電汽輪機帶負(fù)荷穩(wěn)態(tài)運行工況的轉(zhuǎn)子高周疲勞裂紋擴展壽命損傷。

    9 結(jié) 論

    (1) 基于年均壽命損耗構(gòu)建的核電汽輪機轉(zhuǎn)子日歷總壽命計算的理論模型,突破了核電汽輪機轉(zhuǎn)子在高周疲勞與低周疲勞損傷機理交互作用下考慮裂紋萌生與裂紋擴展2個階段的日歷總壽命計算的難題,為核電汽輪機長壽命設(shè)計奠定了理論基礎(chǔ)。

    (2) 建立的核電汽輪機轉(zhuǎn)子在高周疲勞與低周疲勞損傷機理交互作用下日歷總壽命的計算方法中,考慮了核電汽輪機帶負(fù)荷穩(wěn)態(tài)運行工況引起的轉(zhuǎn)子高周疲勞壽命損傷,以及核電汽輪機起動與停機等瞬態(tài)工況引起的轉(zhuǎn)子低周疲勞壽命損傷,日歷總壽命計算方法在工程上可行。

    (3) 應(yīng)用文中給出的核電汽輪機轉(zhuǎn)子在高周疲勞與低周疲勞損傷機理交互作用下日歷總壽命的計算方法,實例結(jié)果表明,該核電汽輪機焊接低壓轉(zhuǎn)子經(jīng)結(jié)構(gòu)優(yōu)化改進后,日歷總壽命由原設(shè)計的55.92 a增長到74.46 a,為保證核電汽輪機的安全服役提供了理論方法和技術(shù)依據(jù)。

    猜你喜歡
    汽輪機裂紋
    裂紋長度對焊接接頭裂紋擴展驅(qū)動力的影響
    一種基于微帶天線的金屬表面裂紋的檢測
    東汽百萬等級汽輪機低壓軸承偏載治理研究
    能源工程(2020年5期)2021-01-04 01:29:00
    Epidermal growth factor receptor rs17337023 polymorphism in hypertensive gestational diabetic women: A pilot study
    微裂紋區(qū)對主裂紋擴展的影響
    淺析給水泵汽輪機跳閘回路改造
    廣西電力(2016年4期)2016-07-10 10:23:38
    汽輪機排汽缸噴水量計算
    預(yù)裂紋混凝土拉壓疲勞荷載下裂紋擴展速率
    汽輪機高壓噴嘴組加工
    汽輪機供熱系統(tǒng)的技改創(chuàng)新
    河南科技(2014年24期)2014-02-27 14:19:40
    国产一区二区在线观看av| 中文字幕人妻丝袜一区二区| 精品一品国产午夜福利视频| 最近中文字幕2019免费版| 91精品三级在线观看| 中国美女看黄片| 成人国语在线视频| 久久精品国产综合久久久| 在线亚洲精品国产二区图片欧美| 国产在线免费精品| 日本色播在线视频| 久久综合国产亚洲精品| av不卡在线播放| 亚洲精品国产av蜜桃| 男女无遮挡免费网站观看| 精品国产一区二区三区久久久樱花| 日本午夜av视频| 亚洲国产精品一区三区| 啦啦啦啦在线视频资源| 日本黄色日本黄色录像| 中文字幕人妻熟女乱码| 亚洲精品av麻豆狂野| 亚洲精品国产区一区二| 九草在线视频观看| 国产国语露脸激情在线看| 脱女人内裤的视频| 成年av动漫网址| 捣出白浆h1v1| 日本色播在线视频| 精品国产国语对白av| 好男人电影高清在线观看| 精品少妇黑人巨大在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 精品一区二区三区av网在线观看 | 在线观看免费日韩欧美大片| 91老司机精品| 国产欧美日韩一区二区三区在线| 亚洲人成77777在线视频| 国产极品粉嫩免费观看在线| 久久久久久久精品精品| 自线自在国产av| 日本黄色日本黄色录像| 老司机午夜十八禁免费视频| 一本大道久久a久久精品| 国产日韩欧美亚洲二区| 在线看a的网站| 日韩中文字幕视频在线看片| 黄色视频在线播放观看不卡| 日韩欧美一区视频在线观看| 久久天躁狠狠躁夜夜2o2o | 亚洲成av片中文字幕在线观看| 国产黄频视频在线观看| 日日夜夜操网爽| 每晚都被弄得嗷嗷叫到高潮| 看十八女毛片水多多多| 亚洲欧洲国产日韩| 欧美人与性动交α欧美软件| 日韩欧美一区视频在线观看| 国产一区二区三区av在线| 黄片播放在线免费| www.自偷自拍.com| 亚洲欧美精品自产自拍| 一级毛片电影观看| 久久久久久免费高清国产稀缺| 国产熟女欧美一区二区| 天天添夜夜摸| 亚洲国产av影院在线观看| 国产成人免费无遮挡视频| 中文字幕制服av| 亚洲精品美女久久av网站| 成人黄色视频免费在线看| 日本av手机在线免费观看| 丝袜在线中文字幕| 99精国产麻豆久久婷婷| 女人精品久久久久毛片| 国产精品欧美亚洲77777| 看免费av毛片| 在线观看免费高清a一片| 国产老妇伦熟女老妇高清| 日韩中文字幕欧美一区二区 | 满18在线观看网站| 亚洲色图 男人天堂 中文字幕| 亚洲av男天堂| 老司机在亚洲福利影院| 精品人妻熟女毛片av久久网站| 一区在线观看完整版| 免费在线观看日本一区| 精品亚洲乱码少妇综合久久| 交换朋友夫妻互换小说| 久久综合国产亚洲精品| 狠狠精品人妻久久久久久综合| 日韩熟女老妇一区二区性免费视频| 一区二区av电影网| 精品福利观看| 久久鲁丝午夜福利片| 精品国产乱码久久久久久小说| 色播在线永久视频| 久久狼人影院| 男人操女人黄网站| 少妇猛男粗大的猛烈进出视频| 丝瓜视频免费看黄片| 精品国产一区二区三区久久久樱花| 各种免费的搞黄视频| 国产老妇伦熟女老妇高清| 考比视频在线观看| 久久性视频一级片| 一本综合久久免费| 日韩 欧美 亚洲 中文字幕| 久久久亚洲精品成人影院| 国产精品一区二区在线不卡| 久久精品国产a三级三级三级| 国产有黄有色有爽视频| 高清av免费在线| 少妇 在线观看| 国产伦理片在线播放av一区| 国产国语露脸激情在线看| 国语对白做爰xxxⅹ性视频网站| 巨乳人妻的诱惑在线观看| 老鸭窝网址在线观看| av天堂久久9| 久久精品久久久久久久性| 国产一区二区三区综合在线观看| 一本色道久久久久久精品综合| 国产黄色视频一区二区在线观看| 国产xxxxx性猛交| 天天躁夜夜躁狠狠久久av| 青春草亚洲视频在线观看| 蜜桃国产av成人99| 99国产综合亚洲精品| xxx大片免费视频| 亚洲欧美精品综合一区二区三区| 美女高潮到喷水免费观看| e午夜精品久久久久久久| 亚洲欧美一区二区三区久久| 国产日韩欧美在线精品| 成年人黄色毛片网站| 国产成人免费无遮挡视频| 久久久国产欧美日韩av| 日韩人妻精品一区2区三区| 亚洲国产精品一区三区| 亚洲国产毛片av蜜桃av| 汤姆久久久久久久影院中文字幕| 午夜日韩欧美国产| 国产欧美日韩一区二区三 | 国产一区二区三区av在线| 精品少妇内射三级| 欧美变态另类bdsm刘玥| bbb黄色大片| 免费在线观看视频国产中文字幕亚洲 | 久久久久久免费高清国产稀缺| 久久精品人人爽人人爽视色| 国产片内射在线| 欧美成人精品欧美一级黄| 欧美日韩综合久久久久久| 一级片'在线观看视频| 国产视频一区二区在线看| 国产一区二区三区综合在线观看| 亚洲国产中文字幕在线视频| 亚洲av综合色区一区| 亚洲人成77777在线视频| 欧美日韩亚洲国产一区二区在线观看 | 宅男免费午夜| 在线av久久热| 晚上一个人看的免费电影| 欧美激情高清一区二区三区| 黑人欧美特级aaaaaa片| 一个人免费看片子| 国产成人影院久久av| 午夜福利视频精品| 国产欧美日韩一区二区三 | svipshipincom国产片| 国产无遮挡羞羞视频在线观看| 国产男女超爽视频在线观看| 9色porny在线观看| 天天影视国产精品| 精品熟女少妇八av免费久了| 一级a爱视频在线免费观看| 九草在线视频观看| 国产片内射在线| 伊人久久大香线蕉亚洲五| 大型av网站在线播放| 黑人欧美特级aaaaaa片| 国产淫语在线视频| 国产极品粉嫩免费观看在线| 国产精品国产三级国产专区5o| videosex国产| 亚洲av电影在线观看一区二区三区| 嫩草影视91久久| 纵有疾风起免费观看全集完整版| 国产成人免费无遮挡视频| 男女无遮挡免费网站观看| 波多野结衣一区麻豆| 高清欧美精品videossex| 另类亚洲欧美激情| 国产精品亚洲av一区麻豆| 国产日韩欧美视频二区| 成人亚洲精品一区在线观看| 久久久国产一区二区| 午夜视频精品福利| 国产精品 国内视频| av国产精品久久久久影院| 亚洲第一av免费看| 国产国语露脸激情在线看| 欧美日韩av久久| 五月开心婷婷网| 国产一区亚洲一区在线观看| 午夜精品国产一区二区电影| 纵有疾风起免费观看全集完整版| 欧美日韩国产mv在线观看视频| 人妻 亚洲 视频| 久9热在线精品视频| 亚洲精品久久午夜乱码| 中文字幕亚洲精品专区| 欧美黄色片欧美黄色片| 亚洲精品自拍成人| 久久久精品94久久精品| 久久久久网色| 国产又爽黄色视频| 成年美女黄网站色视频大全免费| 国产97色在线日韩免费| 亚洲av片天天在线观看| 自线自在国产av| 国产男女超爽视频在线观看| 日韩伦理黄色片| av国产精品久久久久影院| 女人被躁到高潮嗷嗷叫费观| 亚洲国产精品999| 汤姆久久久久久久影院中文字幕| 欧美大码av| 一级黄片播放器| 1024香蕉在线观看| 超色免费av| 国产成人av激情在线播放| 亚洲熟女毛片儿| 大陆偷拍与自拍| 久久午夜综合久久蜜桃| 波多野结衣av一区二区av| 亚洲天堂av无毛| 国产精品熟女久久久久浪| 黄色视频不卡| 中国国产av一级| 十八禁高潮呻吟视频| 超色免费av| 一本大道久久a久久精品| 精品人妻熟女毛片av久久网站| 国产精品人妻久久久影院| 欧美成人精品欧美一级黄| 精品福利永久在线观看| 亚洲欧美清纯卡通| 亚洲国产日韩一区二区| 中国国产av一级| 波多野结衣av一区二区av| 亚洲天堂av无毛| 国产精品熟女久久久久浪| 性高湖久久久久久久久免费观看| 伦理电影免费视频| 国产精品香港三级国产av潘金莲 | 欧美精品高潮呻吟av久久| 欧美精品一区二区免费开放| 亚洲专区国产一区二区| 亚洲七黄色美女视频| 真人做人爱边吃奶动态| 国产伦理片在线播放av一区| 久久久久精品国产欧美久久久 | 免费在线观看影片大全网站 | 亚洲国产最新在线播放| 大香蕉久久网| 国产一区二区 视频在线| 精品国产一区二区三区四区第35| 性少妇av在线| 亚洲专区中文字幕在线| 中文欧美无线码| 亚洲一区中文字幕在线| 精品国产一区二区三区久久久樱花| 高潮久久久久久久久久久不卡| 美女午夜性视频免费| 亚洲伊人久久精品综合| 黄色 视频免费看| 亚洲精品日韩在线中文字幕| a级毛片黄视频| 成在线人永久免费视频| 看免费av毛片| 十八禁网站网址无遮挡| 中文字幕最新亚洲高清| 精品视频人人做人人爽| 中文字幕色久视频| 免费黄频网站在线观看国产| 极品人妻少妇av视频| 超碰97精品在线观看| av网站免费在线观看视频| 亚洲精品国产av成人精品| 一二三四社区在线视频社区8| 午夜激情久久久久久久| 在线观看人妻少妇| 国产成人91sexporn| 亚洲精品中文字幕在线视频| 青草久久国产| 又紧又爽又黄一区二区| 亚洲国产欧美一区二区综合| 99国产精品一区二区三区| 亚洲av电影在线进入| 我的亚洲天堂| 在线精品无人区一区二区三| 两个人看的免费小视频| 国产不卡av网站在线观看| 久久久久视频综合| 欧美av亚洲av综合av国产av| 亚洲免费av在线视频| 狠狠精品人妻久久久久久综合| 亚洲人成电影观看| 麻豆乱淫一区二区| 18在线观看网站| 午夜激情av网站| 狠狠婷婷综合久久久久久88av| 国产成人精品久久二区二区免费| 亚洲中文av在线| 50天的宝宝边吃奶边哭怎么回事| 国产一卡二卡三卡精品| 黑人猛操日本美女一级片| 国产在线观看jvid| 日韩中文字幕视频在线看片| 一个人免费看片子| xxx大片免费视频| 极品少妇高潮喷水抽搐| 一级毛片 在线播放| 三上悠亚av全集在线观看| 91国产中文字幕| 欧美精品亚洲一区二区| 天天躁夜夜躁狠狠久久av| 欧美性长视频在线观看| 日本a在线网址| 久久精品国产综合久久久| 看十八女毛片水多多多| 91精品伊人久久大香线蕉| 黄色视频不卡| 啦啦啦在线免费观看视频4| 男女床上黄色一级片免费看| 天堂8中文在线网| 黄片小视频在线播放| 男女边吃奶边做爰视频| 国产精品国产av在线观看| 日韩av免费高清视频| 欧美精品高潮呻吟av久久| 少妇人妻久久综合中文| 国产在线一区二区三区精| 日韩伦理黄色片| 国产高清视频在线播放一区 | 久久综合国产亚洲精品| 国产成人一区二区在线| 日本av手机在线免费观看| 欧美激情 高清一区二区三区| 欧美精品av麻豆av| 宅男免费午夜| 成人18禁高潮啪啪吃奶动态图| 成年动漫av网址| 大片电影免费在线观看免费| 久久国产精品男人的天堂亚洲| 成人国语在线视频| 国产在线观看jvid| 亚洲国产中文字幕在线视频| 日韩欧美一区视频在线观看| svipshipincom国产片| 青春草视频在线免费观看| 啦啦啦在线观看免费高清www| 新久久久久国产一级毛片| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久人人人人人| bbb黄色大片| 纯流量卡能插随身wifi吗| 一区在线观看完整版| 亚洲精品国产区一区二| 操出白浆在线播放| 国产深夜福利视频在线观看| 女人精品久久久久毛片| 亚洲国产毛片av蜜桃av| 午夜久久久在线观看| 成年动漫av网址| 一边亲一边摸免费视频| 免费看不卡的av| 欧美另类一区| 国产精品久久久人人做人人爽| 久久久国产一区二区| 成年动漫av网址| 日本猛色少妇xxxxx猛交久久| 涩涩av久久男人的天堂| 亚洲色图 男人天堂 中文字幕| 人妻人人澡人人爽人人| 亚洲成人免费av在线播放| 免费人妻精品一区二区三区视频| 在线观看免费高清a一片| 久久久精品免费免费高清| 麻豆国产av国片精品| 亚洲欧美一区二区三区久久| 国产不卡av网站在线观看| 亚洲人成电影免费在线| a级毛片黄视频| av有码第一页| 国产精品 欧美亚洲| 午夜免费成人在线视频| 久久精品成人免费网站| 久久精品国产综合久久久| 欧美人与善性xxx| 各种免费的搞黄视频| 欧美 日韩 精品 国产| 母亲3免费完整高清在线观看| 亚洲国产精品成人久久小说| 一级,二级,三级黄色视频| 午夜两性在线视频| 日本午夜av视频| 欧美+亚洲+日韩+国产| 亚洲国产av影院在线观看| 日本五十路高清| 国产97色在线日韩免费| 午夜福利乱码中文字幕| 欧美久久黑人一区二区| www.999成人在线观看| 啦啦啦中文免费视频观看日本| www.av在线官网国产| 好男人视频免费观看在线| 最近中文字幕2019免费版| 日本欧美国产在线视频| 久久久亚洲精品成人影院| 亚洲成人免费电影在线观看 | 极品少妇高潮喷水抽搐| 久久久久久久国产电影| 女警被强在线播放| av在线app专区| 777米奇影视久久| 男女无遮挡免费网站观看| 精品一区在线观看国产| www.999成人在线观看| 午夜福利视频精品| 中文字幕人妻丝袜制服| 黄色视频在线播放观看不卡| 制服诱惑二区| 欧美日韩亚洲高清精品| e午夜精品久久久久久久| 国产在线观看jvid| 国产精品一区二区在线不卡| 久久鲁丝午夜福利片| 9191精品国产免费久久| 永久免费av网站大全| 欧美精品亚洲一区二区| 久9热在线精品视频| 自线自在国产av| 国产成人欧美在线观看 | 久久久久久久精品精品| 人人妻人人澡人人爽人人夜夜| 美国免费a级毛片| 国产精品人妻久久久影院| 国产精品香港三级国产av潘金莲 | 亚洲人成77777在线视频| 欧美性长视频在线观看| 免费在线观看完整版高清| 成人18禁高潮啪啪吃奶动态图| 精品人妻熟女毛片av久久网站| 男男h啪啪无遮挡| 久久人妻福利社区极品人妻图片 | 欧美日韩一级在线毛片| 精品亚洲成国产av| 亚洲精品国产一区二区精华液| 久久av网站| 啦啦啦在线免费观看视频4| 精品国产乱码久久久久久男人| 亚洲国产最新在线播放| 久久久亚洲精品成人影院| 欧美亚洲日本最大视频资源| 精品一区二区三区四区五区乱码 | 老司机午夜十八禁免费视频| 精品一区二区三卡| 中国国产av一级| 看十八女毛片水多多多| 老司机亚洲免费影院| 亚洲精品日韩在线中文字幕| 久久人人97超碰香蕉20202| 成人国语在线视频| 欧美老熟妇乱子伦牲交| 久久 成人 亚洲| 亚洲九九香蕉| 欧美成人午夜精品| 精品一区二区三区四区五区乱码 | 国产一区二区在线观看av| 国产日韩欧美在线精品| 国产一区二区 视频在线| 亚洲中文av在线| 老司机亚洲免费影院| 大香蕉久久成人网| www.自偷自拍.com| 多毛熟女@视频| 一二三四社区在线视频社区8| 一区二区三区乱码不卡18| 精品福利永久在线观看| 久久久久久人人人人人| 午夜激情久久久久久久| 国产精品二区激情视频| 一本综合久久免费| 亚洲精品自拍成人| 黑人巨大精品欧美一区二区蜜桃| 男女边吃奶边做爰视频| 91精品三级在线观看| 国产在视频线精品| 这个男人来自地球电影免费观看| 亚洲图色成人| 欧美国产精品一级二级三级| 国产欧美亚洲国产| av不卡在线播放| 欧美亚洲 丝袜 人妻 在线| 亚洲av成人不卡在线观看播放网 | 久久人人爽人人片av| 国产熟女午夜一区二区三区| 国产又爽黄色视频| xxxhd国产人妻xxx| 久久女婷五月综合色啪小说| 美女国产高潮福利片在线看| 制服诱惑二区| 日本欧美视频一区| av国产精品久久久久影院| 狠狠婷婷综合久久久久久88av| 建设人人有责人人尽责人人享有的| 在线观看免费午夜福利视频| av欧美777| 亚洲欧美日韩高清在线视频 | 波野结衣二区三区在线| 水蜜桃什么品种好| 亚洲av成人精品一二三区| 国产主播在线观看一区二区 | 国产av精品麻豆| 亚洲精品一卡2卡三卡4卡5卡 | 夫妻午夜视频| 国产三级黄色录像| 亚洲午夜精品一区,二区,三区| 婷婷色综合大香蕉| 欧美97在线视频| 久久青草综合色| 亚洲九九香蕉| 国产精品一二三区在线看| 少妇人妻久久综合中文| 久久久欧美国产精品| 丝袜美足系列| 国产黄色免费在线视频| www.自偷自拍.com| 国产精品一区二区在线不卡| 天天添夜夜摸| 天堂俺去俺来也www色官网| 色婷婷av一区二区三区视频| 久久久欧美国产精品| 欧美日韩国产mv在线观看视频| 1024视频免费在线观看| av在线app专区| 9色porny在线观看| 免费观看人在逋| 久久精品久久久久久久性| 少妇的丰满在线观看| 人妻一区二区av| 老司机深夜福利视频在线观看 | 免费在线观看日本一区| 在线看a的网站| 纯流量卡能插随身wifi吗| 日韩av免费高清视频| 看十八女毛片水多多多| 亚洲精品成人av观看孕妇| 国产成人av教育| 国产午夜精品一二区理论片| 王馨瑶露胸无遮挡在线观看| 久久久亚洲精品成人影院| 精品人妻1区二区| 久久午夜综合久久蜜桃| 免费女性裸体啪啪无遮挡网站| 老汉色∧v一级毛片| 中文精品一卡2卡3卡4更新| 国产精品九九99| 老司机在亚洲福利影院| 少妇人妻久久综合中文| 亚洲,欧美,日韩| 欧美人与性动交α欧美软件| 美女主播在线视频| 久久亚洲精品不卡| 国产精品一区二区在线不卡| 一二三四社区在线视频社区8| 国产免费视频播放在线视频| 国产深夜福利视频在线观看| 国产福利在线免费观看视频| www日本在线高清视频| 欧美在线一区亚洲| 一级a爱视频在线免费观看| 国产精品偷伦视频观看了| 丰满饥渴人妻一区二区三| 美女主播在线视频| 韩国精品一区二区三区| 久久久久国产一级毛片高清牌| 成人午夜精彩视频在线观看| 啦啦啦中文免费视频观看日本| 最近中文字幕2019免费版| 欧美精品啪啪一区二区三区 | 考比视频在线观看| 啦啦啦中文免费视频观看日本| 午夜福利在线免费观看网站| 老司机亚洲免费影院| 国产黄色免费在线视频| 三上悠亚av全集在线观看| av片东京热男人的天堂| 人妻 亚洲 视频| 精品少妇久久久久久888优播| 大陆偷拍与自拍| 一边摸一边抽搐一进一出视频| 两性夫妻黄色片| 人人妻,人人澡人人爽秒播 | 老司机影院成人| 亚洲美女黄色视频免费看| 男人舔女人的私密视频| 十分钟在线观看高清视频www| 日韩视频在线欧美| 我的亚洲天堂| 一级,二级,三级黄色视频| 国产日韩欧美在线精品|