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

    汽車零部件熱成形極限裕度云圖的構(gòu)建與應(yīng)用

    2017-03-01 09:27:18崔俊佳李光耀
    中國(guó)機(jī)械工程 2017年3期
    關(guān)鍵詞:壓邊裕度板料

    崔俊佳 胡 明 李光耀

    1.湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室,長(zhǎng)沙,4100822.智能型新能源汽車協(xié)同創(chuàng)新中心,上海,201804

    汽車零部件熱成形極限裕度云圖的構(gòu)建與應(yīng)用

    崔俊佳1,2胡 明1李光耀1,2

    1.湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室,長(zhǎng)沙,4100822.智能型新能源汽車協(xié)同創(chuàng)新中心,上海,201804

    在熱沖壓工藝中,熱沖壓鋼板處于高溫狀態(tài),不同位置溫度分布有明顯差異,傳統(tǒng)方法難以評(píng)價(jià)沖壓件的成形性,考慮溫度的三維成形極限圖能有效評(píng)價(jià)熱沖壓件的成形性能,但無法判斷熱沖壓件的安全裕度。針對(duì)上述問題,在考慮溫度的三維成形極限圖的基礎(chǔ)上,首先提出了成形裕度的計(jì)算方法;然后結(jié)合數(shù)值模擬和理論分析手段,構(gòu)建了適用于評(píng)價(jià)高強(qiáng)鋼板熱沖壓的安全裕度云圖;最后對(duì)某汽車B柱熱沖壓不同工藝過程進(jìn)行了仿真分析,并與試驗(yàn)進(jìn)行了對(duì)比。結(jié)果表明,在不同工藝參數(shù)和板料尺寸下,構(gòu)建的裕度云圖均能有效預(yù)測(cè)高強(qiáng)鋼板熱沖壓件的開裂部位和程度。

    高強(qiáng)鋼板;熱沖壓;成形極限圖;成形裕度云圖

    0 引言

    隨著汽車輕量化的發(fā)展,汽車車身中高強(qiáng)度鋼板的應(yīng)用越來越廣泛。高強(qiáng)度鋼板不僅能夠減輕車重,而且可以提高汽車的燃油經(jīng)濟(jì)性,減少尾氣排放。然而,高強(qiáng)度鋼板在室溫下的成形性較差,易產(chǎn)生較大的回彈,所需成形力也較大[1-2]。熱沖壓技術(shù)能夠有效解決上述問題,獲得幾何形狀復(fù)雜、尺寸精度高的高強(qiáng)鋼零件,且經(jīng)過熱沖壓技術(shù)后,零件抗拉強(qiáng)度可達(dá)1500 MPa[3],因而受到汽車行業(yè)的廣泛關(guān)注。

    熱沖壓件開裂缺陷的預(yù)測(cè)是其得以推廣和發(fā)展的必要基礎(chǔ)。目前國(guó)內(nèi)外學(xué)者對(duì)熱成形件的斷裂和成形性進(jìn)行了大量研究。LI等[4]建立了基于位錯(cuò)密度統(tǒng)一形式的超塑性斷裂本構(gòu)方程,并通過有限元模擬預(yù)測(cè)了熱成形過程中的缺陷。盈亮等[5]基于連續(xù)介質(zhì)損傷力學(xué)建立了耦合損傷的熱成形本構(gòu)方程,并利用有限元預(yù)測(cè)了某B柱熱成形中的破裂現(xiàn)象。借鑒室溫下成形極限圖(forming limit diagram,F(xiàn)LD)的建立和使用方法,SHI等[6]提出了熱成形極限圖(thermal forming limit diagraph,TFLD),將600 ℃、700 ℃、800 ℃下的成形極限曲線擬合成曲面,將位于曲面以下的應(yīng)變-溫度狀態(tài)點(diǎn)判斷為安全,落在曲面上判斷為頸縮,落在曲面以上判斷為破裂。在文獻(xiàn)[6]方法的基礎(chǔ)上,CUI等[7]提出了一種溫度范圍更廣且考慮了成形過程中相變影響的三維成形極限圖,通過有限元數(shù)值模擬,成功預(yù)測(cè)了某B柱的熱成形破裂缺陷。上述研究只能定性判定熱沖壓件是否破裂,并不能評(píng)價(jià)破裂的程度,或者未破裂區(qū)域的安全程度。NACEUR等[8]將成形極限曲線(forming limit curve,FLC)作為板料破裂的判斷標(biāo)準(zhǔn),用單元成形后的主應(yīng)變與FLC處對(duì)應(yīng)的主應(yīng)變之差來判定板料破裂的程度。WEI等[9]針對(duì)冷沖壓成形的優(yōu)化設(shè)計(jì),將應(yīng)變狀態(tài)點(diǎn)到FLC的距離定義為約束條件,用于評(píng)價(jià)試件數(shù)值模擬的沖壓性能。上述方法定量地評(píng)價(jià)了板材冷沖壓性能,但均沒有涉及熱沖壓件成形性的評(píng)價(jià)。

    目前國(guó)內(nèi)外對(duì)于使用FLC定量評(píng)價(jià)高強(qiáng)度鋼板熱沖壓性能的研究較少,對(duì)使用三維成形極限圖評(píng)價(jià)高強(qiáng)度鋼板熱沖壓成形缺陷也不夠系統(tǒng)?;诖?,本文首先提出了一種預(yù)測(cè)熱沖壓缺陷的方法;然后,在溫度相關(guān)的FLC基礎(chǔ)上,利用數(shù)值模擬通過彩色云圖表示成形裕度,結(jié)合裕度云圖和成形極限圖評(píng)價(jià)熱沖壓件的成形性;最后,將該方法用于某汽車B柱的熱沖壓數(shù)值模擬中,對(duì)比各種工藝參數(shù)下零件的沖壓性能,并進(jìn)行了試驗(yàn)驗(yàn)證。

    1 成形裕度云圖的建立方法

    FLC被廣泛用于預(yù)測(cè)金屬板材成形的破裂失效預(yù)測(cè)[10]。FLC是指沖壓件在不同應(yīng)變路徑下發(fā)生頸縮或破裂之前可以達(dá)到的極限應(yīng)變線,可以作為分析沖壓成形過程中材料是否失效的判據(jù)[11]。

    成形裕度是指材料在變形一定程度之后還具備的變形能力。在熱沖壓中,熱沖壓鋼板在不同溫度下的材料性能不同,導(dǎo)致其在不同溫度下成形極限曲線也有所不同。因此,判斷熱沖壓中某點(diǎn)的成形性能時(shí),要考慮該點(diǎn)的溫度狀態(tài),使用這一溫度下的成形極限對(duì)其進(jìn)行評(píng)價(jià)。

    熱成形極限圖是由不同溫度、等溫條件下的成形極限曲線組成的。若將這些曲線擬合成一空間曲面,即形成熱沖壓的成形極限面[7]。

    在空間曲面擬合時(shí),要明確因變量與自變量之間的關(guān)系,本文中極限主應(yīng)變是因變量,溫度和次應(yīng)變是自變量,自變量x、y與因變量z之間的函數(shù)關(guān)系可以寫成

    (1)

    (2)

    (3)

    式中,pij為多項(xiàng)式系數(shù)。

    圖1 熱成形極限曲面的擬合結(jié)果Fig.1 Result of fitting hot forming limit surface

    借鑒文獻(xiàn)[9]中冷沖壓成形性能的表示方法,將應(yīng)變-溫度狀態(tài)點(diǎn)至該單元溫度下成形極限曲線的最小距離表示為板料的成形裕度。設(shè)曲線外一點(diǎn)p0(x0, y0, z0),將溫度坐標(biāo)y0代入擬合曲面,得y0溫度下的FLC:z=z(x, y0),但y0溫度下的成形極限曲線次應(yīng)變的取值范圍無法直接確定,故采用Lagrange多項(xiàng)式插值法確定在此溫度下的次應(yīng)變的取值范圍,Lagrange插值多項(xiàng)式可以表示為[11]

    (4)

    其中,li(x)是Lagrange基函數(shù),其表達(dá)式為

    li(x)=

    (5)

    li(x)具有以下性質(zhì):

    (6)

    i,j=0,1,…,n

    應(yīng)用上述Lagrange多項(xiàng)式插值方法,由已知溫度下次應(yīng)變的范圍,求得對(duì)應(yīng)溫度下FLC的次應(yīng)變范圍[xd,xu]。在xz平面坐標(biāo)系下,設(shè)曲線上一點(diǎn)為p(x,z),對(duì)應(yīng)溫度p0(x0,z0)下成形極限曲線的斜率kp和直線pp0的斜率kpp0的乘積kpkpp0=-1,得

    (7)

    聯(lián)立下式

    z(x,y0)=z

    (8)

    可解得p(x, z)。將成形裕度M定義如下:

    (9)

    其中,d為p到p0點(diǎn)的距離。F(x0,y0,z0)=z(x0,y0)-z0>0表示應(yīng)變-溫度狀態(tài)點(diǎn)在FLC的下方,如圖2中p01點(diǎn)的位置所示,p1為FLC:z=z(x,y01)上離其最近的點(diǎn),直線距離p1p01即定義為該溫度下的成形裕度。F(x0,y0,z0)=z(x0,y0)-z0<0表示應(yīng)變-溫度狀態(tài)點(diǎn)在FLC

    圖2 單元成形裕度的計(jì)算方法Fig.2 Computing method of forming margin for single element

    的上方,如圖2中p02點(diǎn)的位置所示,p2為FLC:z=z(x, y02)上離其最近的點(diǎn),此時(shí)判定板料發(fā)生破裂。F(x0,y0,z0)=z(x0,y0)-z0=0表示應(yīng)變-溫度狀態(tài)點(diǎn)在FLC上,為臨界破裂。將上述方法編寫成MATLAB程序,讀取數(shù)值模型上所有單元的主次應(yīng)變和溫度信息計(jì)算成形裕度值,由PAM-STAMP2G后處理最終形成裕度云圖,成形裕度的生成程序流程圖見圖3。

    圖3 所有單元成形裕度計(jì)算程序流程圖Fig.3 Flow chart of forming margin computing program for all elements

    2 有限元模擬

    2.1 有限元模型的建立

    選用厚度為1.5 mm的熱沖壓硼鋼板BR1500HS作為研究對(duì)象,其主要元素含量見表1。

    表1 硼鋼板BR1500HS的化學(xué)成分(質(zhì)量分?jǐn)?shù))

    采用PAM-STAMP 2G進(jìn)行熱沖壓建模,板料使用殼單元進(jìn)行離散,網(wǎng)格的大小為13 mm,單元個(gè)數(shù)為4209,在數(shù)值模擬過程中對(duì)網(wǎng)格進(jìn)行自動(dòng)細(xì)化,凸凹模間隙設(shè)置為1.65 mm,沖壓工具設(shè)置為剛體。模具、板料之間傳熱系數(shù)與接觸壓力有關(guān)[12],其關(guān)系如圖4所示。可以看出,隨著板料與模具間壓力的增大,傳熱系數(shù)先增大后趨于穩(wěn)定。

    圖4 板料-工具間傳熱系數(shù)隨壓力的變化Fig.4 Pressure dependent heat transfer coefficient between blank and tools

    試驗(yàn)中使用石墨對(duì)模具表面進(jìn)行潤(rùn)滑,以減小板料和工具之間的摩擦,數(shù)值模擬中將摩擦因數(shù)設(shè)置為0.1。板材在應(yīng)變速率為0.1 s-1、不同溫度θ下的真實(shí)應(yīng)力-應(yīng)變曲線如圖5所示??梢钥闯觯瑴囟葘?duì)熱沖壓鋼板的力學(xué)性能有很大的影響,溫度越高,板料流動(dòng)性能越好。建立的有限元模型及各部件擺放的位置如圖6所示。

    圖5 應(yīng)變速率為0.1 s-1時(shí)的真實(shí)應(yīng)力應(yīng)變曲線Fig.5 True strain-stress curves at strain rate 0.1 s-1

    圖6 沖壓成形的有限元模型Fig.6 Finite element model of the stamping process

    熱沖壓工藝參數(shù)(如壓邊力、沖壓速度、板料初始溫度)對(duì)沖壓性能有較大的影響[13],因此,為了評(píng)價(jià)各工藝參數(shù)對(duì)熱沖壓性能的影響,對(duì)各沖壓參數(shù)下的數(shù)值模擬進(jìn)行研究,分析各工況下板材的成形裕度,工藝參數(shù)見表2。

    表2 沖壓成形工藝參數(shù)Tab.2 Process parameter of stamping

    2.2 有限元模擬結(jié)果

    2.2.1 板料形狀對(duì)熱沖壓性能的影響

    為了驗(yàn)證成形裕度云圖的準(zhǔn)確性,選取兩種板料形狀進(jìn)行數(shù)值模擬。板料初始溫度為900 ℃、沖壓速度為350 mm/s、壓邊力為2000 kN。板料調(diào)整前后的尺寸及形狀如圖7所示。

    (a)板料調(diào)整前

    (b)板料調(diào)整后圖7 板料調(diào)整前后的形狀Fig.7 Blank shape before and after adjusting

    圖8a和圖8b分別為板料形狀調(diào)整前后的數(shù)值模擬結(jié)果。圖8a的結(jié)果表明,B柱的側(cè)壁出現(xiàn)了嚴(yán)重的破裂現(xiàn)象,過多的壓料面限制了板料在B柱側(cè)壁和拐角處的流動(dòng)。減小B柱零件修邊線以外的坯料后(圖8b),由成形裕度云圖可以看出,應(yīng)變-溫度狀態(tài)點(diǎn)顯著下移,板料形狀調(diào)整前后成形裕度最小值分別為-1.73和-0.34,成形裕度為負(fù),說明板料依然存在破裂現(xiàn)象,經(jīng)過分析發(fā)現(xiàn)破裂部分在修邊線以外(圖9)。修邊線以外成形裕度均為正,說明在板料尺寸調(diào)整后,B柱側(cè)壁處破裂現(xiàn)象消失,板料的熱沖壓性能有所改善。

    (a)板料形狀調(diào)整前

    (b)板料形狀調(diào)整后圖8 板料形狀調(diào)整前后數(shù)值模擬結(jié)果Fig.8 Results of numerical modeling before and after blank adjusting

    圖9 B柱零件修邊線以外破裂情況Fig.9 Cracking Out of B-pillar Trimming Line

    綜上分析可以看出,使用成形極限圖僅能判斷成形件是否破裂,而成形裕度云圖不僅可以在數(shù)值模型上顯示破裂的區(qū)域,而且可以通過不同的顏色表示破裂的程度。另外,對(duì)于未破裂的區(qū)域,也可以根據(jù)顏色的分布來預(yù)測(cè)材料還可變形的程度。

    2.2.2 壓邊力對(duì)熱沖壓性能的影響

    選取板料的初始溫度為900 ℃,沖壓速度為250 mm/s,壓邊力為1500~3000 kN進(jìn)行數(shù)值模擬,分析不同壓邊力下B柱的沖壓性能。圖10~圖13所示為不同的壓邊力下的數(shù)值模擬結(jié)果。當(dāng)壓邊力在1500~2500 kN之間時(shí),B柱修邊線以內(nèi)部分單元對(duì)應(yīng)的應(yīng)變-溫度狀態(tài)點(diǎn)在成形極限曲面之下,說明未發(fā)生開裂(圖10~圖12)。當(dāng)壓邊力增大到3000 kN時(shí),B柱直壁發(fā)生了較嚴(yán)重的破裂,最小成形裕度達(dá)到-2.38,成形極限曲面之上的應(yīng)變-溫度狀態(tài)點(diǎn)明顯增多(圖13)。壓邊力越大,對(duì)板料的約束增加,板料的成形性能越差。

    圖10 壓邊力1500 kN下的數(shù)值模擬結(jié)果Fig.10 Numerical modeling results at 1500 kN pressure force conditions

    圖11 壓邊力2000 kN下的數(shù)值模擬結(jié)果Fig.11 Numerical modeling results at 2000 kN pressure force conditions

    圖12 壓邊力2500 kN下的數(shù)值模擬結(jié)果Fig.12 Numerical modeling results at 2500 kN pressure force conditions

    圖13 壓邊力3000 kN下的數(shù)值模擬結(jié)果Fig.13 Numerical modeling results at 3000 kN pressure force conditions

    2.2.3 沖壓速度對(duì)熱沖壓性能的影響

    選取板料初始溫度為900 ℃,壓邊力為2000 kN進(jìn)行數(shù)值模擬,分析沖壓速度對(duì)成形性能的影響。圖14~圖17所示為不同沖壓速度下的數(shù)值模擬結(jié)果。當(dāng)沖壓速度在250~400 mm/s之間變化時(shí),最小成形裕度隨沖壓速度的增大略有減小(表3),在400 mm/s沖壓速度時(shí)達(dá)到最小值,此時(shí)B柱成形性能較好。其原因是沖壓速度增大,導(dǎo)致成形時(shí)間變短,板料溫度高,材料流動(dòng)性好,但沖壓速度繼續(xù)增大會(huì)導(dǎo)致材料來不及流動(dòng),對(duì)變形較大區(qū)域難以補(bǔ)料。綜合上述原因,沖壓速度應(yīng)在低速范圍內(nèi)盡量大。

    圖14 沖壓速度250 mm/s下的數(shù)值模擬結(jié)果Fig.14 Numerical modeling results at 250 mm/s die velocity

    圖15 沖壓速度300 mm/s下的數(shù)值模擬結(jié)果Fig.15 Numerical modeling results at 300 mm/s die velocity

    圖16 沖壓速度350 mm/s下的數(shù)值模擬結(jié)果Fig.16 Numerical modeling results at 350 mm/s die velocity

    圖17 沖壓速度400 mm/s下的數(shù)值模擬結(jié)果Fig.17 Numerical modeling results at 400 mm/s die velocity

    沖壓速度(mm/s)250300350400成形裕度-0.40~0.39-0.38~0.40-0.34~0.40-0.30~0.40

    2.2.4 板料初始溫度對(duì)熱沖壓性能的影響

    選取修改后的板料,將沖壓速度和壓邊力分別設(shè)置為350 mm/s和2000 kN進(jìn)行數(shù)值模擬。圖18~圖21所示為不同初始溫度下的數(shù)值模擬結(jié)果??梢钥闯?,溫度從750 ℃升高到900 ℃時(shí),沖壓件的最小裕度由-1.92增大到-0.34(表4)。板料初始溫度為750 ℃時(shí),B柱拐角附近出現(xiàn)嚴(yán)重的開裂現(xiàn)象(圖18);當(dāng)初始溫度升高到800 ℃時(shí),開裂現(xiàn)象明顯改善(圖19);當(dāng)板料初始溫度繼續(xù)升高到850 ℃以上時(shí),B柱的成形性能有較大的改善,最小成形裕度分別提高到-0.36,成形極限曲面之上的應(yīng)變-溫度狀態(tài)點(diǎn)均靠近成形極限曲面(圖20)。與板料初始溫度850 ℃時(shí)相比,900 ℃時(shí)的成形裕度略有提高(圖21)。原因是溫度越高,材料流動(dòng)性越強(qiáng),成形極限也相應(yīng)提高。

    圖18 板料初始溫度750 ℃下的數(shù)值模擬結(jié)果Fig.18 Numerical modeling results at 750 ℃ blank initial temperature

    圖19 板料初始溫度800 ℃下的數(shù)值模擬結(jié)果Fig.19 Numerical modeling results at 800 ℃ blank initial temperature

    圖20 板料初始溫度850 ℃下的數(shù)值模擬結(jié)果Fig.20 Numerical modeling results at 850 ℃ blank initial temperature

    圖21 板料初始溫度900 ℃下的數(shù)值模擬結(jié)果Fig.21 Numerical modeling results at 900 ℃ blank initial temperature

    初始溫度(℃)750800850900最小成形裕度-1.92-1.39-0.36-0.34

    3 熱成形裕度云圖的應(yīng)用與評(píng)價(jià)

    為了對(duì)熱成形裕度云圖進(jìn)行應(yīng)用和評(píng)價(jià),同時(shí)驗(yàn)證數(shù)值模擬的正確性,進(jìn)行B柱熱沖壓試驗(yàn)。板料的初始溫度均為900 ℃,試驗(yàn)中板料的形狀及工藝參數(shù)見表5。

    表5 B柱熱沖壓工藝參數(shù)Tab.5 B-pillar hot stamping process parameter

    首先研究板料尺寸對(duì)熱沖壓成形性的影響。當(dāng)板料尺寸為初始尺寸,沖壓速度為350 mm/s、壓邊力為2000 kN時(shí),進(jìn)行數(shù)值模擬與試驗(yàn)結(jié)果的對(duì)比(圖22)。將B柱側(cè)壁分為A1、A2、A3三個(gè)區(qū)域,從試驗(yàn)結(jié)果來看,A1、A2和A3段均發(fā)生破裂,A3段破裂程度段小于A1和A2段破裂程度段。在圖8a中,側(cè)壁處區(qū)域裕度均為負(fù)值,數(shù)值見表6,A1段由右向左成形裕度有減小的趨勢(shì),A1至A3段成形裕度范圍的變化與試驗(yàn)中破裂程度基本吻合。

    (a)仿真結(jié)果

    (b)試驗(yàn)結(jié)果圖22 B柱側(cè)壁處有破裂情況對(duì)比Fig.22 Contrast of B-pillar side wall when cracking

    區(qū)域A1區(qū)域A2區(qū)域A3成形裕度-1.13~-0.52-0.82~-0.22-0.52~0.22

    圖23所示為工藝參數(shù)不變,板料修改后B柱側(cè)壁部分尺寸(修改尺寸1)的數(shù)值模擬與試驗(yàn)結(jié)果的對(duì)比。圖23a中B柱側(cè)壁處裕度在0.08~0.29之間,即仿真中沒有出現(xiàn)破裂,圖23b中對(duì)應(yīng)區(qū)域也沒有出現(xiàn)破裂,試驗(yàn)和仿真結(jié)果一致。

    (a)仿真結(jié)果

    (b)試驗(yàn)結(jié)果圖23 B柱側(cè)壁處無破裂情況對(duì)比Fig.23 Contrast of B-pillar side wall when no cracking

    圖24所示為沖壓速度為350 mm/s、壓邊力為2000 kN時(shí),B柱拐角處數(shù)值模擬與試驗(yàn)結(jié)果的對(duì)比。圖24a中拐角處區(qū)域的裕度值為-0.02~0.29,即仿真中沒有破裂的區(qū)域,圖24b試驗(yàn)中拐角處對(duì)應(yīng)區(qū)域亦沒有出現(xiàn)破裂,仿真和試驗(yàn)結(jié)果吻合較好。

    (a)仿真結(jié)果

    (b)試驗(yàn)結(jié)果圖24 B柱拐角處未裂的預(yù)測(cè)Fig.24 Forecast of no Cracking at B-pillar Corner

    圖25所示為板料修改尺寸(修改尺寸2)后,沖壓速度降低至250 mm/s、壓邊力增大至3000 kN時(shí),B柱拐角處數(shù)值模擬與試驗(yàn)的對(duì)比。圖25a成形裕度云圖中,B柱左右拐角處區(qū)域的成形裕度在-0.4~0之間,對(duì)應(yīng)的應(yīng)變-溫度狀態(tài)點(diǎn)位于成形極限曲面之上(圖13),判定為破裂。圖25b試驗(yàn)中對(duì)應(yīng)的這兩個(gè)位置也出現(xiàn)破裂的現(xiàn)象。

    (a)仿真結(jié)果

    (b)試驗(yàn)結(jié)果圖25 B柱拐角處破裂的預(yù)測(cè)Fig.25 Forecast of crack at B-pillar corner

    綜上所述,選用不同規(guī)格尺寸和工藝參數(shù),所建立的成形裕度云圖均能有效地預(yù)測(cè)開裂位置和程度。

    4 結(jié)論

    (1)本文基于三維成形極限圖,提出了高強(qiáng)度鋼板熱沖壓成形裕度的計(jì)算方法,并將該方法編寫成程序;結(jié)合PAM-TAMP 2G的后處理程序,構(gòu)建了熱成形裕度云圖。

    (2)板料的初始溫度、壓邊力和板料的幾何形狀對(duì)沖壓性能均有較大的影響。隨著溫度的升高,材料流動(dòng)性變好,板料成形極限提升;隨著壓邊力的增大,B柱熱沖壓性能變差。

    (3)改變板料尺寸、沖壓速度和壓邊力等參數(shù)對(duì)建立的裕度云圖進(jìn)行驗(yàn)證,發(fā)現(xiàn)B柱側(cè)壁、兩拐角處的仿真結(jié)果和試驗(yàn)結(jié)果基本吻合,板料破裂程度變化趨勢(shì)基本一致。

    [1] NADERI M, UTHAISANGSUK V, PRAHL U, et al.A Numerical and Experimental Investigation into Hot Stamping of Boron Alloyed Heat Treated Steels[J].Steel Research International,2008, 79(2): 77-84.

    [2] 王敏,張春,張宇,等.汽車高強(qiáng)鋼縱梁熱沖壓試驗(yàn)研究[J]. 中國(guó)機(jī)械工程,2016,27(3):391-397. WANG Min, ZHANG Chun, ZHANG Yu, et al. Experiments on Hot Stamping of Car Caring of High Strength Steel [J]. China Mechanical Engineering, 2016,27(3):391-397.

    [3] 莊百亮,單忠德,姜 超,等.汽車防撞梁超高強(qiáng)鋼熱成形工藝研究[J].中國(guó)機(jī)械工程,2012,23(2):225-228. ZHUANG Bailiang, SHAN Zhongde, JIANG Chao, et al.Research on Ultra-high-strength Hot Stamping Process of Automobive Anti-collision Beam[J]. China Mechanical Engineering, 2012, 23(2):225-228.

    [4] LI Nan, SUN Chaoyang, GUO Ning, et al. Damage Investigation of Boron Steel at Hot Stamping Conditions[C]//11th International Conference on Technology of Plasticity.Nagoya, 2014: 19-24.

    [5] 盈亮,史棟勇,胡平,等.基于連續(xù)介質(zhì)損傷力學(xué)的高強(qiáng)度鋼板熱成形性數(shù)值預(yù)測(cè)[J].機(jī)械工程學(xué)報(bào),2016,52(4):36-44. YING Liang, SHI Dongyong, HU Ping, et al. Numerical Prediction of the Formability during Hot Forming of High Strength Steel Based on Continuum Damage Mechanics[J]. Journal of Mechanical Engineering, 2016,52(4):36-44.

    [6] SHI D Y, YING L, HU P, et al. Experimental and Numerical Determination of Thermal Forming Limit Diagrams (TFLD) of High Strength Steel 22MnB5[C]// 11th International Conference on Numerical Methods in Industrial Forming Processes:Shenyang, 2013: 406-413.

    [7] CUI Junjia, SUN Guangyong, XU Junrui, et al. A Method to Evaluate the Formability of High-strength Steel in Hot Stamping[J]. Materials and Design, 2015, 77: 95-109.

    [8] NACEUR H, DELAMéZIERE A, BATOZ J L, et al. Some Improvements on the Optimum Process Design in Deep Drawing Using the Inverse Approach[J]. Journal of Materials Processing Technology, 2004, 146(2): 250-262.

    [9] WEI Donglai, CUI Zhenshan, CHEN Jun. Optimization and Tolerance Prediction of Sheet Metal Forming Process Using Response Surface Model[J].Computational Materials Science, 2008, 42(2): 228-233.

    [10] ZHUANG Weimin, ZHANG Mengxi, CHEN Yanhong. Error Analysis of FLC Experimental Data at Warm/Hot Stamping Conditions[J]. Chinese Journal of Mechanical Engineering, 2014, 27(4): 730-737.

    [11] ZHANG Huawei, LIU Xianghua, LIU Lizhong, et al. Study on Nonuniform Deformation of Tailor Rolled Blank during Uniaxial Tension[J]. Acta Metall. Sin. (Engl. Lett.), 2015, 28(9): 1198-1204.

    [12] CHANG Ying, TANG Xinghui, ZHAO Kunmin, et al. Investigation of the Factors Influencing the Interfacial Heat Transfer Coefficient in Hot Stamping [J]. Journal of Materials Processing Technology, 2016, 228: 25-33.

    [13] 高云凱,鄧有志,曹 偉.超高強(qiáng)度鋼車身B柱加強(qiáng)板熱成形工藝參數(shù)多目標(biāo)優(yōu)化[J].中國(guó)機(jī)械工程,2011,22(5):621-629. GAO Yunkai, DENG Youzhi, CAO Wei. Multi-objective Optimization for Ultra High Strength Steel B-pillar of Car Body Hot Forming Process Parameters[J]. China Mechanical Engineering, 2011,22(5):621-629.

    (編輯 陳 勇)

    作者簡(jiǎn)介:崔俊佳,男,1983年生。湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室副教授。主要研究方向?yàn)橄冗M(jìn)成形技術(shù)、異種金屬連接。E-mail: cuijunjia@hnu.edu.cn。胡 明,男,1991年生。湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室碩士研究生。李光耀(通信作者),男,1963年生。湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室教授、博士研究生導(dǎo)師。E-mail:gyli@hnu.edu.cn。

    Establishment and Applications of Thermal Forming Limit Margin Field Graph for Automobile Components

    CUI Junjia1,2HU Ming1LI Guangyao1,2

    1.State Key Laboratory of Advanced Design and Manufacture for Vehicle Body of Hunan University, Changsha,410082 2.Collaborative Innovation Center of Intelligent New Energy Vehicle, Shanghai,201804

    In hot stamping processes, the sheet maintained in a high temperature state and the temperature varied at different positions, for which traditional method was difficult to evaluate the sheet formability. Temperature dependent 3D forming limit graph was a new method to evaluate formability of hot stamping components, however, the prediction of forming margin remained a blank. To solve the problems discussed above, a way to calculate the forming margin was proposed herein based on the temperature dependent 3D forming limit graph. Then, combining numerical simulations and theoretical analyses, a safe margin field graph for evaluating high strength steel sheet in hot stamping was established. Using the safe margin field graph, different hot stamping processes of an automobile B pillar were analyzed by numerical simulation and were compared to experiments. Results show that the proposed forming margin field graph may effectively predict the fracture positions and cracked degrees of HSS components in hot stamping at different process parameters and blank shapes.

    high strength steel(HSS) sheet; hot stamping; forming limit diagram; forming margin field graph

    2016-03-31

    國(guó)家自然科學(xué)基金資助重點(diǎn)項(xiàng)目(61232014);國(guó)家自然科學(xué)基金資助項(xiàng)目(51405149)

    TG306DOI:10.3969/j.issn.1004-132X.2017.03.018

    陳清根,男,1990年生。南昌航空大學(xué)航空制造工程學(xué)院碩士研究生。主要研究方向?yàn)樾澢霃綇濐^的冷推彎成形。徐雪峰,男,1983年生。南昌航空大學(xué)航空制造工程學(xué)院講師。E-mail:xfwinzy@163.com。馬媛媛,女,1988年生。中航工業(yè)洪都航空工業(yè)集團(tuán)鈑金廠工程師。李玲玲,女,1974年生。中航工業(yè)洪都航空工業(yè)集團(tuán)制造工程部高級(jí)工程師。趙 爽,男,1982年生。中航工業(yè)洪都航空工業(yè)集團(tuán)部裝一廠高級(jí)工程師。徐 龍,男,1963年生。中航工業(yè)洪都航空工業(yè)集團(tuán)制造工程部研究員。

    猜你喜歡
    壓邊裕度板料
    沖床板料輸送機(jī)分離裝置的動(dòng)力學(xué)和運(yùn)動(dòng)學(xué)仿真分析
    拉延模具板料定位設(shè)計(jì)
    薄板成形電磁壓邊技術(shù)研究進(jìn)展
    漸進(jìn)成形A3003鋁板減薄帶分析及數(shù)值模擬研究
    重力對(duì)中臺(tái)在厚板沖壓自動(dòng)生產(chǎn)線中的應(yīng)用
    鍛造與沖壓(2018年6期)2018-03-20 02:10:11
    柔性壓邊多點(diǎn)成形技術(shù)
    基于DFIG可用無功裕度的風(fēng)電場(chǎng)無功電壓控制方法
    三環(huán)路核電廠的抗震裕度評(píng)價(jià)
    基于ANN模型的在線電壓穩(wěn)定裕度評(píng)估
    電流互感器磁飽和裕度及其試驗(yàn)裝置的探討
    久久久久久久久中文| 亚洲国产欧洲综合997久久,| 久久午夜福利片| 99久久精品国产亚洲精品| 99热这里只有是精品50| 亚洲精品在线美女| 日韩欧美一区二区三区在线观看| 国产精品,欧美在线| 日日干狠狠操夜夜爽| 一级黄片播放器| 亚洲在线观看片| 欧美日韩综合久久久久久 | 一a级毛片在线观看| 18禁在线播放成人免费| 老司机午夜福利在线观看视频| av专区在线播放| 亚洲电影在线观看av| 久久99热这里只有精品18| 亚洲欧美日韩卡通动漫| 国产精品久久久久久久久免 | 欧美黄色片欧美黄色片| 久久久国产成人免费| 男女那种视频在线观看| 两个人的视频大全免费| 国产三级黄色录像| 99精品久久久久人妻精品| 国产精品一区二区免费欧美| 我的老师免费观看完整版| 在线免费观看的www视频| 成人无遮挡网站| 亚洲美女搞黄在线观看 | 日日夜夜操网爽| 亚洲电影在线观看av| 国产精品一及| 国产欧美日韩一区二区精品| 麻豆一二三区av精品| 99久久久亚洲精品蜜臀av| 一进一出好大好爽视频| 在线观看午夜福利视频| 久久久精品大字幕| 免费在线观看影片大全网站| 我要看日韩黄色一级片| 男人狂女人下面高潮的视频| 天堂网av新在线| 亚洲国产高清在线一区二区三| 欧美午夜高清在线| 一本一本综合久久| 日韩大尺度精品在线看网址| 又粗又爽又猛毛片免费看| 久久精品影院6| 1024手机看黄色片| 人人妻人人看人人澡| 国产三级中文精品| 搡老妇女老女人老熟妇| 国产中年淑女户外野战色| 亚洲天堂国产精品一区在线| 在线播放无遮挡| 久久久色成人| 午夜精品在线福利| 日韩有码中文字幕| 又爽又黄无遮挡网站| 午夜精品一区二区三区免费看| 免费看光身美女| 久久欧美精品欧美久久欧美| 色综合亚洲欧美另类图片| 欧美一区二区亚洲| 精品一区二区三区视频在线| 桃色一区二区三区在线观看| 免费在线观看日本一区| 淫妇啪啪啪对白视频| 精品久久久久久久久久免费视频| 人妻丰满熟妇av一区二区三区| 国产伦一二天堂av在线观看| 九九在线视频观看精品| 97碰自拍视频| 国产精品久久久久久精品电影| 色哟哟·www| 亚洲aⅴ乱码一区二区在线播放| 亚洲五月婷婷丁香| 国产精品亚洲av一区麻豆| 国产av一区在线观看免费| 久久伊人香网站| 国产大屁股一区二区在线视频| 在线播放国产精品三级| 成人高潮视频无遮挡免费网站| 亚洲精品成人久久久久久| 欧美zozozo另类| 蜜桃久久精品国产亚洲av| 亚洲欧美日韩东京热| av欧美777| 国产精品久久视频播放| 天堂av国产一区二区熟女人妻| 日韩人妻高清精品专区| 亚洲aⅴ乱码一区二区在线播放| 蜜桃久久精品国产亚洲av| 在线观看免费视频日本深夜| 国产精品日韩av在线免费观看| x7x7x7水蜜桃| 国语自产精品视频在线第100页| 亚洲美女视频黄频| 久久久国产成人免费| 小蜜桃在线观看免费完整版高清| 婷婷色综合大香蕉| 给我免费播放毛片高清在线观看| 欧美潮喷喷水| 亚洲成人精品中文字幕电影| 天天躁日日操中文字幕| 日本一本二区三区精品| 亚洲美女黄片视频| 久久久久九九精品影院| 精品久久久久久久人妻蜜臀av| 国产精品影院久久| 中亚洲国语对白在线视频| 一本精品99久久精品77| 国产精品久久久久久亚洲av鲁大| 一级毛片久久久久久久久女| 十八禁国产超污无遮挡网站| 国产野战对白在线观看| 国产成人影院久久av| 国产成人a区在线观看| 一进一出抽搐gif免费好疼| 露出奶头的视频| 欧美3d第一页| 亚洲国产高清在线一区二区三| 床上黄色一级片| 免费在线观看成人毛片| 婷婷亚洲欧美| 99热只有精品国产| 欧美丝袜亚洲另类 | 亚洲第一区二区三区不卡| 有码 亚洲区| 日本撒尿小便嘘嘘汇集6| 国产探花极品一区二区| 亚洲熟妇中文字幕五十中出| 久久99热6这里只有精品| 91九色精品人成在线观看| www日本黄色视频网| 成人无遮挡网站| 成人欧美大片| 国产真实伦视频高清在线观看 | 亚洲美女搞黄在线观看 | 亚洲精品一卡2卡三卡4卡5卡| 国产高清激情床上av| 日本一二三区视频观看| 成人永久免费在线观看视频| 日韩欧美三级三区| 免费看日本二区| 亚洲专区国产一区二区| 日本精品一区二区三区蜜桃| 亚洲精品456在线播放app | 丰满的人妻完整版| 两人在一起打扑克的视频| 国产精品国产高清国产av| 88av欧美| 欧美另类亚洲清纯唯美| 亚洲精品在线美女| 欧美激情在线99| 国产成人a区在线观看| 午夜福利在线在线| 亚洲欧美清纯卡通| 搞女人的毛片| 亚洲综合色惰| 免费在线观看亚洲国产| 国产精品久久久久久久电影| 一个人免费在线观看电影| 99久久成人亚洲精品观看| 每晚都被弄得嗷嗷叫到高潮| 中亚洲国语对白在线视频| 国产探花在线观看一区二区| 亚洲成人精品中文字幕电影| 国产真实乱freesex| 内射极品少妇av片p| 免费av毛片视频| 亚洲人成电影免费在线| 日本成人三级电影网站| 亚洲色图av天堂| 亚洲激情在线av| 久9热在线精品视频| 赤兔流量卡办理| 最新中文字幕久久久久| 国产淫片久久久久久久久 | 精品久久国产蜜桃| 欧美成人免费av一区二区三区| 全区人妻精品视频| 性色av乱码一区二区三区2| 天堂影院成人在线观看| 亚洲成人精品中文字幕电影| 成人av在线播放网站| 内射极品少妇av片p| 亚洲av美国av| 老司机福利观看| 国产免费一级a男人的天堂| 欧美一区二区亚洲| av专区在线播放| 欧美bdsm另类| 很黄的视频免费| 一个人看的www免费观看视频| 十八禁网站免费在线| 免费看日本二区| 午夜福利在线观看免费完整高清在 | 一二三四社区在线视频社区8| 亚洲人成电影免费在线| 欧美日韩综合久久久久久 | 欧美成人a在线观看| 免费一级毛片在线播放高清视频| 中文字幕人成人乱码亚洲影| 国产人妻一区二区三区在| 免费看光身美女| 长腿黑丝高跟| 丰满人妻熟妇乱又伦精品不卡| 午夜福利高清视频| 久久午夜福利片| 国产亚洲精品av在线| 亚洲综合色惰| 欧美极品一区二区三区四区| 欧美又色又爽又黄视频| 国产淫片久久久久久久久 | 我要看日韩黄色一级片| 日韩大尺度精品在线看网址| www.色视频.com| 国内精品久久久久久久电影| 三级国产精品欧美在线观看| 色播亚洲综合网| 免费一级毛片在线播放高清视频| 日韩欧美三级三区| 日本在线视频免费播放| 自拍偷自拍亚洲精品老妇| 国产精华一区二区三区| 白带黄色成豆腐渣| 深爱激情五月婷婷| a级一级毛片免费在线观看| 国产高清激情床上av| 一级黄片播放器| 亚洲人成网站在线播放欧美日韩| 丝袜美腿在线中文| 搞女人的毛片| 欧美性猛交╳xxx乱大交人| 久久人妻av系列| 一级a爱片免费观看的视频| 在现免费观看毛片| 国产视频一区二区在线看| 天美传媒精品一区二区| 国产精品人妻久久久久久| 亚洲激情在线av| 性欧美人与动物交配| 色哟哟哟哟哟哟| 女人十人毛片免费观看3o分钟| 天堂av国产一区二区熟女人妻| 国产免费av片在线观看野外av| 免费在线观看亚洲国产| 精品人妻偷拍中文字幕| 亚洲五月婷婷丁香| 又紧又爽又黄一区二区| 亚洲国产精品成人综合色| 免费av观看视频| 宅男免费午夜| 91麻豆av在线| 亚洲乱码一区二区免费版| 亚洲欧美日韩高清专用| 国产视频一区二区在线看| 男女之事视频高清在线观看| 国产午夜精品论理片| 国产av不卡久久| 久久久久精品国产欧美久久久| av中文乱码字幕在线| 给我免费播放毛片高清在线观看| 欧美日韩综合久久久久久 | 桃红色精品国产亚洲av| .国产精品久久| 欧美一区二区亚洲| 在线观看舔阴道视频| 免费高清视频大片| 男女那种视频在线观看| 久久中文看片网| 婷婷丁香在线五月| 美女大奶头视频| bbb黄色大片| 动漫黄色视频在线观看| 99riav亚洲国产免费| 757午夜福利合集在线观看| 琪琪午夜伦伦电影理论片6080| 99久久成人亚洲精品观看| 婷婷色综合大香蕉| 国产精品久久久久久久久免 | 午夜福利视频1000在线观看| 亚洲专区国产一区二区| 天天躁日日操中文字幕| 伦理电影大哥的女人| 国产精品自产拍在线观看55亚洲| 看片在线看免费视频| 久久中文看片网| 亚洲美女搞黄在线观看 | 嫩草影视91久久| 美女cb高潮喷水在线观看| 天天一区二区日本电影三级| 51午夜福利影视在线观看| 看片在线看免费视频| 亚洲av免费在线观看| 校园春色视频在线观看| 3wmmmm亚洲av在线观看| 最新在线观看一区二区三区| 亚洲乱码一区二区免费版| 51国产日韩欧美| 999久久久精品免费观看国产| 色播亚洲综合网| 老司机深夜福利视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 免费观看精品视频网站| АⅤ资源中文在线天堂| 啦啦啦观看免费观看视频高清| 他把我摸到了高潮在线观看| 级片在线观看| 一进一出好大好爽视频| 国产综合懂色| 久久精品久久久久久噜噜老黄 | 久久久久久久久久成人| 综合色av麻豆| 中文字幕av在线有码专区| 少妇的逼好多水| 国产在线精品亚洲第一网站| 88av欧美| 免费观看精品视频网站| 麻豆国产97在线/欧美| 十八禁国产超污无遮挡网站| 深夜精品福利| 久久久久久久亚洲中文字幕 | 欧美潮喷喷水| 国产精品三级大全| 少妇的逼水好多| 国产一区二区在线av高清观看| 精品久久久久久久久久免费视频| 男女床上黄色一级片免费看| 中文字幕久久专区| 亚洲精品粉嫩美女一区| 丁香六月欧美| 日本a在线网址| 亚洲中文日韩欧美视频| 国产一区二区三区在线臀色熟女| 国产成人a区在线观看| 久久久久久久久大av| 国产精品久久电影中文字幕| 欧美bdsm另类| 色综合欧美亚洲国产小说| 久久国产乱子免费精品| 色av中文字幕| 国产精品av视频在线免费观看| 婷婷丁香在线五月| 国产高清激情床上av| 蜜桃久久精品国产亚洲av| 一a级毛片在线观看| 91麻豆av在线| 久久国产乱子伦精品免费另类| 一级黄片播放器| 男人舔奶头视频| 18禁黄网站禁片午夜丰满| 午夜福利高清视频| 国内精品久久久久久久电影| 国产午夜精品久久久久久一区二区三区 | 成人毛片a级毛片在线播放| 欧美成狂野欧美在线观看| 在线天堂最新版资源| 18美女黄网站色大片免费观看| 成年版毛片免费区| 亚洲aⅴ乱码一区二区在线播放| 在线十欧美十亚洲十日本专区| 在线观看美女被高潮喷水网站 | 亚洲熟妇中文字幕五十中出| av专区在线播放| 国产一区二区激情短视频| 窝窝影院91人妻| 麻豆国产97在线/欧美| 日日摸夜夜添夜夜添小说| 国产av不卡久久| 国产欧美日韩精品一区二区| 噜噜噜噜噜久久久久久91| 亚洲午夜理论影院| 亚洲av一区综合| 国产午夜精品论理片| 制服丝袜大香蕉在线| 一个人看视频在线观看www免费| 超碰av人人做人人爽久久| 久久久久免费精品人妻一区二区| 免费看美女性在线毛片视频| 伊人久久精品亚洲午夜| 日本 欧美在线| 小说图片视频综合网站| 精品欧美国产一区二区三| 日韩欧美 国产精品| bbb黄色大片| 婷婷精品国产亚洲av在线| 国产亚洲欧美98| 2021天堂中文幕一二区在线观| 免费在线观看成人毛片| 亚洲精品成人久久久久久| 亚洲avbb在线观看| 一进一出好大好爽视频| 免费观看精品视频网站| 舔av片在线| 亚洲无线观看免费| 日韩av在线大香蕉| 国产视频一区二区在线看| 久久九九热精品免费| 身体一侧抽搐| x7x7x7水蜜桃| 日韩中字成人| 国产aⅴ精品一区二区三区波| 亚洲,欧美精品.| 日韩大尺度精品在线看网址| 九九热线精品视视频播放| 久久久久免费精品人妻一区二区| 成人午夜高清在线视频| 国产视频一区二区在线看| 欧美乱妇无乱码| ponron亚洲| 少妇人妻精品综合一区二区 | 日韩有码中文字幕| 色在线成人网| 欧美高清性xxxxhd video| 精品一区二区免费观看| 中文字幕久久专区| 听说在线观看完整版免费高清| 欧美日韩国产亚洲二区| 男女视频在线观看网站免费| 国产午夜福利久久久久久| av女优亚洲男人天堂| 神马国产精品三级电影在线观看| 日本撒尿小便嘘嘘汇集6| 丁香欧美五月| 午夜久久久久精精品| 一二三四社区在线视频社区8| av专区在线播放| 国产精品自产拍在线观看55亚洲| 久久久久国产精品人妻aⅴ院| 性插视频无遮挡在线免费观看| 99视频精品全部免费 在线| 1024手机看黄色片| 神马国产精品三级电影在线观看| 欧美成人免费av一区二区三区| 久久精品国产亚洲av香蕉五月| 国产色爽女视频免费观看| АⅤ资源中文在线天堂| 一本综合久久免费| 毛片一级片免费看久久久久 | 亚洲av成人不卡在线观看播放网| 亚洲一区高清亚洲精品| 精品一区二区免费观看| 成人欧美大片| 嫩草影院新地址| 欧美日韩综合久久久久久 | 久久草成人影院| 久久婷婷人人爽人人干人人爱| 久久国产精品人妻蜜桃| 成人特级av手机在线观看| 内射极品少妇av片p| 久久天躁狠狠躁夜夜2o2o| 日本黄大片高清| 精品无人区乱码1区二区| 亚洲精品成人久久久久久| 国产在视频线在精品| 久久人人爽人人爽人人片va | 少妇的逼好多水| 搡老熟女国产l中国老女人| 免费人成视频x8x8入口观看| 日本黄大片高清| 俺也久久电影网| 国产免费一级a男人的天堂| 99久久精品国产亚洲精品| 男人舔女人下体高潮全视频| 男女那种视频在线观看| 亚洲性夜色夜夜综合| 黄色一级大片看看| 真人一进一出gif抽搐免费| 亚洲av成人不卡在线观看播放网| 婷婷丁香在线五月| 男女床上黄色一级片免费看| 国产午夜精品论理片| 亚洲精华国产精华精| 12—13女人毛片做爰片一| 国产69精品久久久久777片| 久久99热这里只有精品18| 99精品在免费线老司机午夜| 亚洲av一区综合| 麻豆成人av在线观看| 一本一本综合久久| 亚洲成人精品中文字幕电影| 国产aⅴ精品一区二区三区波| 女生性感内裤真人,穿戴方法视频| 国产亚洲精品av在线| 男女视频在线观看网站免费| 久久久久久九九精品二区国产| 欧美中文日本在线观看视频| 免费电影在线观看免费观看| 在线播放无遮挡| 最近最新免费中文字幕在线| 中文字幕av成人在线电影| 免费看日本二区| 丰满人妻一区二区三区视频av| 国产av麻豆久久久久久久| 亚洲七黄色美女视频| 亚洲 欧美 日韩 在线 免费| 国内精品久久久久久久电影| 久久久成人免费电影| 亚洲av中文字字幕乱码综合| 免费人成视频x8x8入口观看| 欧美一区二区精品小视频在线| 欧美黄色淫秽网站| 99热只有精品国产| 欧美+日韩+精品| 男人的好看免费观看在线视频| 神马国产精品三级电影在线观看| 亚洲精品乱码久久久v下载方式| АⅤ资源中文在线天堂| 人妻制服诱惑在线中文字幕| 青草久久国产| 成人三级黄色视频| 高潮久久久久久久久久久不卡| 深爱激情五月婷婷| 一区二区三区激情视频| 日本撒尿小便嘘嘘汇集6| 色综合亚洲欧美另类图片| 赤兔流量卡办理| 日本 欧美在线| 久久精品人妻少妇| 直男gayav资源| 久久久久久大精品| 国产毛片a区久久久久| 麻豆国产av国片精品| av国产免费在线观看| 亚洲在线自拍视频| 最近中文字幕高清免费大全6 | 午夜影院日韩av| 最近最新中文字幕大全电影3| 亚洲狠狠婷婷综合久久图片| 国产精品1区2区在线观看.| 无人区码免费观看不卡| 久久这里只有精品中国| 国产精品人妻久久久久久| 成人无遮挡网站| 午夜免费成人在线视频| 精品人妻一区二区三区麻豆 | eeuss影院久久| 国产在线精品亚洲第一网站| 一夜夜www| 我的老师免费观看完整版| 色综合欧美亚洲国产小说| 哪里可以看免费的av片| 少妇高潮的动态图| 51国产日韩欧美| 他把我摸到了高潮在线观看| 中文资源天堂在线| 国产亚洲精品久久久com| 亚洲精品日韩av片在线观看| 99精品在免费线老司机午夜| 老鸭窝网址在线观看| 国产精品久久久久久亚洲av鲁大| 国产在视频线在精品| 欧美一区二区亚洲| 国内精品美女久久久久久| 成熟少妇高潮喷水视频| 欧美一区二区国产精品久久精品| 亚洲色图av天堂| x7x7x7水蜜桃| 欧美黄色片欧美黄色片| 欧美日韩瑟瑟在线播放| 国产成人aa在线观看| 久久这里只有精品中国| 国产精品不卡视频一区二区 | 国产三级中文精品| 一本精品99久久精品77| 搡老妇女老女人老熟妇| 国内少妇人妻偷人精品xxx网站| 91麻豆av在线| 欧美日韩中文字幕国产精品一区二区三区| 欧美最新免费一区二区三区 | 国产免费一级a男人的天堂| 中文字幕免费在线视频6| 久久久久久大精品| 夜夜看夜夜爽夜夜摸| 特大巨黑吊av在线直播| 久久精品综合一区二区三区| 国产精品亚洲一级av第二区| 男人狂女人下面高潮的视频| 精品久久国产蜜桃| 国产老妇女一区| 成年女人毛片免费观看观看9| 九色成人免费人妻av| 变态另类成人亚洲欧美熟女| 女人被狂操c到高潮| 床上黄色一级片| 在线a可以看的网站| 又爽又黄无遮挡网站| 亚洲精品粉嫩美女一区| 国产老妇女一区| 蜜桃亚洲精品一区二区三区| 琪琪午夜伦伦电影理论片6080| 精品不卡国产一区二区三区| 性欧美人与动物交配| 成年女人永久免费观看视频| 在线a可以看的网站| 国产精品伦人一区二区| 精品久久久久久久久av| 亚洲专区国产一区二区| 国产成人欧美在线观看| 国产精品日韩av在线免费观看| 九九在线视频观看精品| 色综合欧美亚洲国产小说| a级毛片a级免费在线| 一区二区三区激情视频| 精品99又大又爽又粗少妇毛片 | 欧美黄色淫秽网站| 亚洲精品粉嫩美女一区| 在线天堂最新版资源| 校园春色视频在线观看| 18禁在线播放成人免费| ponron亚洲|