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

    壓氣機(jī)葉片輥軋模具型腔前滑補(bǔ)償方法

    2017-11-23 05:58:09靳淇超汪文虎蔣睿嵩趙德中崔康熊一峰
    航空學(xué)報(bào) 2017年1期
    關(guān)鍵詞:模腔型腔軋輥

    靳淇超,汪文虎,蔣睿嵩,趙德中,崔康,熊一峰

    壓氣機(jī)葉片輥軋模具型腔前滑補(bǔ)償方法

    靳淇超,汪文虎,蔣睿嵩*,趙德中,崔康,熊一峰

    西北工業(yè)大學(xué) 現(xiàn)代設(shè)計(jì)與集成制造技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,西安 710072

    輥軋模腔是葉片無(wú)余量輥軋成形的關(guān)鍵,然而壓氣機(jī)葉片輥軋成形過(guò)程中存在前滑現(xiàn)象,造成葉片沿輥軋方向精度偏差,需要在模具型腔設(shè)計(jì)過(guò)程中補(bǔ)償前滑量,以實(shí)現(xiàn)葉片沿積疊高度的精確成形。針對(duì)葉片輥軋前滑現(xiàn)象,本文提出并研究了基于輥軋前滑補(bǔ)償?shù)膲簹鈾C(jī)葉片輥軋模具型腔優(yōu)化設(shè)計(jì)方法。首先,在分析前滑成因及葉片截面前滑表征模型的基礎(chǔ)上,研究了前滑補(bǔ)償機(jī)理并建立了前滑補(bǔ)償模型,即型腔中心角與葉片積疊高度的對(duì)應(yīng)關(guān)系。在此基礎(chǔ)上,通過(guò)提取表征葉片的工藝模型截面線族,順次計(jì)算截面前滑值并基于前滑補(bǔ)償模型對(duì)輥軋模具型腔中心角進(jìn)行修正。然后,基于修正后的型腔中心角,建立葉片工藝模型截面線族到型腔截面線族的空間扇態(tài)映射法則并進(jìn)行截面線映射變換,進(jìn)而基于型腔截面線重構(gòu)了基于前滑補(bǔ)償?shù)娜~片輥軋模具型腔。最后,通過(guò)高置信度數(shù)值計(jì)算方法比較了前滑補(bǔ)償模腔和直接空間幾何映射模腔輥軋成形葉片的積疊高度。結(jié)果表明,優(yōu)化后的模具型腔能夠有效提高葉片積疊軸方向上的成形精度。

    壓氣機(jī)葉片;輥軋;前滑補(bǔ)償;模腔優(yōu)化設(shè)計(jì);積疊高度

    航空發(fā)動(dòng)機(jī)壓氣機(jī)的主要功能是將進(jìn)入發(fā)動(dòng)機(jī)的空氣逐級(jí)加壓并將壓縮空氣提供給燃燒室。葉片是進(jìn)行空氣壓縮的關(guān)鍵部件,其有著厚度薄、型面扭曲、品類繁多、數(shù)量巨大,制造精度要求極高等特點(diǎn)。無(wú)余量輥軋成形由于具有生產(chǎn)效率高、成形精度穩(wěn)定、在成形過(guò)程中產(chǎn)生均勻的加工硬化可以提高葉片的疲勞性能和服役性能[1]等優(yōu)點(diǎn),已成為高壓壓氣機(jī)葉片成形主要方法之一[2-4]。

    輥軋成形屬于模壓成形,模腔精度是保證輥軋成形精度和質(zhì)量的關(guān)鍵,其設(shè)計(jì)精度直接影響產(chǎn)品的尺寸精度、裝配性能和服役性能。特別是針對(duì)壓氣機(jī)葉片輥軋成形,要求成形精度高、質(zhì)量穩(wěn)定,因此精確穩(wěn)定的模具型腔是高精度輥軋葉片成形的基礎(chǔ),這就要求模具設(shè)計(jì)過(guò)程中綜合考

    慮葉身截面形狀、對(duì)應(yīng)截面間距和輥軋成形過(guò)程中前滑回彈等因素。21世紀(jì)初,為實(shí)現(xiàn)太行發(fā)動(dòng)機(jī)高壓壓氣機(jī)葉片成形要求,國(guó)內(nèi)引進(jìn)了無(wú)余量葉片輥軋成形技術(shù),并在輥軋成形結(jié)構(gòu)件的組織性能、力學(xué)性能等方面開(kāi)展了研究[5-6]。隨著輥軋葉片廣泛應(yīng)用于航空發(fā)動(dòng)機(jī),研究工作拓展到輥軋?jiān)O(shè)備開(kāi) 發(fā) 驗(yàn) 證[7-8]、輥 軋 工 藝 優(yōu) 化[9-10]、輥 軋 有限元模擬[11-12]、輥軋成形力學(xué)分析[13]等方面。然而葉片輥軋過(guò)程中存在前滑現(xiàn)象,在輥軋成形出口處,成形葉片截面速度大于模具型腔截面速度。以葉片截面線空間幾何映射變換獲得的型腔未考慮前滑對(duì)輥軋葉片成形的誤差補(bǔ)償,導(dǎo)致截面線變化的葉片相鄰截面間距擴(kuò)大,引起輥軋方向精度偏差。雖然采用修模-試模-再修模的方法,可以滿足生產(chǎn)需求,但是這種粗放的工藝延長(zhǎng)了模具定型周期。目前針對(duì)壓氣機(jī)葉片輥軋模具型腔設(shè)計(jì)和型腔優(yōu)化的研究極少。Sedighi和 Mahmoodi[14]通過(guò)實(shí)驗(yàn)和有限元方法進(jìn)行壓氣機(jī)輥軋葉片厚度誤差補(bǔ)償和優(yōu)化,提出軋輥軸心距調(diào)整算法,實(shí)現(xiàn)了葉片厚度方向精確成形。陳雅和周紅梅[15]將葉片截面理論數(shù)據(jù)、冷軋模截面測(cè)量數(shù)據(jù)、冷軋模截面理論數(shù)據(jù)進(jìn)行對(duì)比分析,對(duì)冷軋模的毛邊面和過(guò)渡面進(jìn)行了優(yōu)化,使其曲率接近一致,同時(shí)加大過(guò)渡面的寬度,徹底消除了葉片臨近前后緣處型面沿弦向皺紋,并縮短30%的修模時(shí)間。李文平等[16]提出一種三維鋁合金板料沖壓工件的沖壓模腔優(yōu)化方法,采用位移補(bǔ)償方法修正模腔型面,縮短了模具定型周期。在葉片輥軋成形中,模具型腔結(jié)構(gòu)和精度是影響成形精度最重要的因素,同時(shí)工藝參數(shù)、裝配精度和毛坯狀態(tài)等也是影響成形精度的重要因素[17]。裴景東[18]在分析輥軋成形工藝的基礎(chǔ)上,建立了葉片輥軋工藝模型,并研究了工藝模型型面映射到模腔型面的變換過(guò)程,實(shí)現(xiàn)了壓氣機(jī)葉片輥軋模腔CAD建模。王淵彬等[19]基于葉片輥軋成形過(guò)程提出了輥軋模腔數(shù)字化建模方法。以葉片中弧面曲率變化為基準(zhǔn),規(guī)劃葉片截面線位置和數(shù)量,對(duì)截面線族進(jìn)行空間映射變換進(jìn)而構(gòu)建型腔曲面,并基于CAD軟件開(kāi)發(fā)程序?qū)崿F(xiàn)葉片模具型腔自動(dòng)化建模。

    葉片輥軋成形時(shí)軋輥?zhàn)鐾鶑?fù)性回轉(zhuǎn)運(yùn)動(dòng),毛坯在一對(duì)逆向旋轉(zhuǎn)的軋輥擠壓作用下發(fā)生塑性變形,獲得給定尺寸形狀的葉片,最后模腔與構(gòu)件分離并返回初始位置。因此模具型腔不是簡(jiǎn)單的葉片型面映射到型腔曲面的空間幾何扇態(tài)分布,需要對(duì)輥軋成形誤差進(jìn)行補(bǔ)償,才能建立更精確的模具型腔模型。由于對(duì)成形誤差機(jī)理的研究匱乏,傳統(tǒng)的模腔設(shè)計(jì)優(yōu)化是基于誤差(實(shí)驗(yàn)測(cè)量的誤差或有限元分析的誤差)的反變形優(yōu)化和修模,具有迭代次數(shù)多,模腔定型周期長(zhǎng)等缺點(diǎn)。本文針對(duì)輥軋前滑引起的成形葉片對(duì)應(yīng)截面線間距擴(kuò)大現(xiàn)象,分析輥軋前滑成因并建立輥軋模腔前滑補(bǔ)償模型,即軋輥型腔截面線中心角與葉片截面線積疊高度對(duì)應(yīng)關(guān)系,精確計(jì)算模腔截面中心角;基于前滑補(bǔ)償模型建立葉片型面點(diǎn)到型腔型面點(diǎn)的映射法則,進(jìn)行葉片截面線族到型腔截面線族的影射變換,實(shí)現(xiàn)截面線族空間扇態(tài)分布并通過(guò)截面線族建立模具型腔,既保證出口截面形狀,又補(bǔ)償輥軋過(guò)程中的前滑,確保葉片截面尺寸精度同時(shí)提高了葉片積疊軸向精度,實(shí)現(xiàn)無(wú)余量輥軋葉片精確成形。

    1 葉片輥軋模具型腔前滑補(bǔ)償方法

    無(wú)余量輥軋成形要求葉片截面線與對(duì)應(yīng)的模具型腔截面同時(shí)出現(xiàn)在輥軋變形區(qū)出口截面,然而輥軋前滑使得對(duì)應(yīng)截面線在輥軋出口處沿輥軋方向速度存在差異,導(dǎo)致葉片積疊軸方向出現(xiàn)成形誤差,因此建立截面線積疊高度與模具型腔中心角對(duì)應(yīng)關(guān)系是進(jìn)行型腔設(shè)計(jì)前滑補(bǔ)償?shù)年P(guān)鍵。

    1.1 截面前滑表征

    毛坯沿輥軋方向運(yùn)動(dòng),在輥軋變形區(qū)受軋輥擠壓發(fā)生塑性變形,最終在接觸分離截面處獲得目標(biāo)構(gòu)件截面形狀,輥軋變形區(qū)內(nèi)塑性變形是輥軋成形的關(guān)鍵,而前滑區(qū)是保證成形精度核心區(qū)域。軋制過(guò)程中在毛坯厚度方向受到擠壓的金屬,一部分沿縱向流動(dòng)形成延伸;而另一部分金屬沿橫向流動(dòng)形成展寬。軋件的延伸是被壓下的金屬沿軋輥入口和出口兩個(gè)方向流動(dòng)的結(jié)果。根據(jù)輥軋基礎(chǔ)理論,在輥軋接觸弧內(nèi),軋件水平速度與軋輥水平速度相等的平面稱為中性面,輥軋變形區(qū)被中性面分為前滑區(qū)和后滑區(qū)兩部分,如圖1所示,圖中α為接觸角,γ為中性角,R為軋輥半徑,v為軋輥圓周速度,D為上下軋輥軸心距,H為毛坯厚度,vH為毛坯在輥軋入口處的速度,h為構(gòu)件厚度,vh為軋件在輥軋出口處的速度,ω為軋輥旋轉(zhuǎn)角速度。在前滑區(qū)內(nèi),沿著輥軋方向毛坯和軋輥的接觸質(zhì)點(diǎn)的速度越來(lái)越快,且毛坯質(zhì)點(diǎn)的速度增長(zhǎng)更大,因此構(gòu)件在軋輥出口處的速度vh大于軋輥圓周速度v,這種現(xiàn)象叫做前滑,其大小用前滑值Sh表示[20],其表達(dá)式為

    如果將構(gòu)件截面線空間幾何映射建立的模腔,未進(jìn)行前滑補(bǔ)償,單位時(shí)間內(nèi)軋輥轉(zhuǎn)角為ω,型腔截面間距為v,構(gòu)件成形長(zhǎng)度為vh,成形構(gòu)件截面單位時(shí)間內(nèi)沿著輥軋方向精度誤差為vh-v,如果構(gòu)件截面無(wú)變化,經(jīng)過(guò)裁剪對(duì)其使用性能無(wú)影響。然而截面變化構(gòu)件,將擴(kuò)大其相鄰截面間距,導(dǎo)致沿輥軋方向上的成形偏差,因此輥軋模腔設(shè)計(jì)過(guò)程中必須進(jìn)行前滑補(bǔ)償。

    對(duì)于截面變化類構(gòu)件,不同截面形狀存在差異,因此前滑值各不相同,建立表征構(gòu)件型面變化的截面線族,通過(guò)分析求解各個(gè)截面的前滑值,進(jìn)而建立工藝模型截面積疊高度與模腔截面的中心角對(duì)應(yīng)關(guān)系,實(shí)現(xiàn)前滑補(bǔ)償。當(dāng)輥軋工藝固化后,壓下率、摩擦系數(shù)和軋輥直徑都確定,前滑值Sh是構(gòu)件厚度h的函數(shù),即Sh=F(h),因此葉片截面前滑值為[21]

    1.2 葉片輥軋前滑補(bǔ)償模型

    型腔截面線與上下軋輥軸線所在平面的夾角稱為中心角,葉片截面線高度與模腔截面線的中心角一一對(duì)應(yīng),計(jì)算截面線對(duì)應(yīng)的中心角θk是保證模腔精度的關(guān)鍵。前滑補(bǔ)償要求葉片工藝模型截面線與對(duì)應(yīng)的模具型腔截面線同時(shí)運(yùn)動(dòng)到輥軋變形出口處,從而保證葉片縱向和橫向尺寸,前滑補(bǔ)償算法示意圖如圖2所示,Z0為葉片葉根截面線,對(duì)應(yīng)的型腔截面線為θ0;葉片總高度為Z,對(duì)應(yīng)的截面線為Zn,對(duì)應(yīng)型腔截面線中心角為θn;第k-1條截面線葉片積疊高度為Zk-1,映射變換后模具型腔截面線的中心角為θk-1,假設(shè)其同時(shí)出現(xiàn)在輥軋出口;第k條截面線葉片積疊高度為Zk,映射變換后模具型腔截面線的中心角為θk,其中θ′k為第k條截面線無(wú)前滑補(bǔ)償僅空間幾何變換位置,其中心角θ′k大于θk,是造成成形葉片積疊高度誤差的原因,因此準(zhǔn)確求解葉片截面線對(duì)應(yīng)的型腔截面線中心角是提高葉片成形積疊高度的核心問(wèn)題。

    葉片輥軋成形初始,葉片工藝模型第一條截面線與模具型腔第一條截面線在輥軋出口處重合。假設(shè)葉片第k-1條截面線Zk-1對(duì)應(yīng)的型腔中心角θk-1已知,第k-1條截面線與其對(duì)應(yīng)的型腔截面線同時(shí)出現(xiàn)在輥軋出口。則葉片第k條截面線與對(duì)應(yīng)型腔截面線同時(shí)出現(xiàn)在輥軋出口處的條件如式(3)所示,第k條截面線與第k-1條截面線間距為Zk-Zk-1,第k條截面線在變形出口處的速度為,第k條截面線對(duì)應(yīng)型腔中心角θk與第k-1條截面線對(duì)應(yīng)型腔中心角θk-1差值為θk-θk-1,為保證成形精度,相鄰截面成形時(shí)間差與軋輥對(duì)應(yīng)相鄰截面轉(zhuǎn)動(dòng)到輥軋出口處的時(shí)間差應(yīng)相等,即葉片和模具型腔相鄰截面線出現(xiàn)在輥軋出口處時(shí)間差為

    基于式(1)將式(3)進(jìn)行變換,則輥軋模具型腔圓心角的遞推關(guān)系為

    因此,輥軋葉片第k條截面線積疊軸高度Zk對(duì)應(yīng)的模具型腔圓心角θk為

    因此對(duì)于截面變化類構(gòu)件,如壓氣機(jī)葉片,可順次求解表征工藝模型的截面線的前滑值,進(jìn)而建立對(duì)應(yīng)截面線積疊高度與模腔截面中心角關(guān)系,精確控制型腔中心角,從而實(shí)現(xiàn)模具設(shè)計(jì)過(guò)程中的前滑補(bǔ)償和葉片積疊高度精確成形。

    2 考慮前滑的模具型腔優(yōu)化設(shè)計(jì)

    模具型腔設(shè)計(jì)是葉片輥軋成形的關(guān)鍵,其設(shè)計(jì)精度直接影響葉片成形精度,必須以葉片工藝模型、輥軋毛坯、輥軋工藝等信息為設(shè)計(jì)依據(jù)。本文針對(duì)葉片工藝模型,建立表征工藝模型的截面線族并按照式(2)順次求每一條截面線前滑量,以式(5)求工藝模型截面線對(duì)應(yīng)的型腔截面線中心角,進(jìn)而基于截面線映射法則實(shí)現(xiàn)工藝模型截面線族到模腔截面線族的映射變換,通過(guò)型腔截面線族建立葉片輥軋模具型腔曲面。

    2.1 工藝模型截面線族提取

    1)工藝模型建模

    葉片在輥軋成形過(guò)程中,需要在葉片前后緣頭添加溢料邊和在葉根葉尖添加工藝余量。為保證葉片前后緣頭尺寸精度并減少加工余量,雙向延伸中弧線l,沿葉片厚度方向分別偏置延伸后中弧線的端點(diǎn)t,t為前后緣曲線擬合圓的半徑,連接葉盆葉背曲線端點(diǎn)與偏置點(diǎn)并光順,如圖3(b)所示,沿著弦向延伸偏置點(diǎn)形成溢料邊。葉根截面線直紋外推lb,葉尖截面線直紋外推lt,針對(duì)本文葉片模型,lb為7mm,lt為5mm,建立葉根葉尖工藝余量。葉片模型、前后緣頭工藝添加和工藝模型如圖3所示。

    2)提取截面線族

    根據(jù)葉片CAD設(shè)計(jì)模型截面線,建立截面線中弧線,通過(guò)中弧線族建立葉片中弧面,中弧面能夠反映葉片在縱向的彎曲程度和積疊軸向的扭曲情況。由于壓氣機(jī)葉片葉身扭曲程度不均勻,局部區(qū)域扭曲程度較大,沿積疊軸方向重新規(guī)劃截面線位置,建立能夠表征中弧面的截面線族,如圖4(a)所示,并在相同積疊高度截取工藝模型葉盆葉背截面線,建立工藝模型截面線族,如圖4(b)所示。

    2.2 模腔建模映射法則

    將工藝模型截面線族分別繞軋輥軸線旋轉(zhuǎn),使其呈空間扇態(tài)分布,實(shí)現(xiàn)截面線族在軋輥表面的重新分布,映射法則是截面線族扇態(tài)分布的關(guān)鍵。葉片沿著葉身方向有扭曲,為防止輥軋出口處構(gòu)件與軋輥干涉,影響葉片輥軋成形精度,軋輥半徑小于軋輥在該處的曲率半徑。然而軋輥半徑較小時(shí),模腔中心角較大,成形精度不穩(wěn)定。針對(duì)本文葉片CAD模型和輥軋?jiān)O(shè)備參數(shù),軋輥中心距D取136mm。對(duì)工藝模型截面線族逐一進(jìn)行點(diǎn)離散,相鄰離散點(diǎn)對(duì)應(yīng)的曲線弧高控制在閥值Δ范圍內(nèi)(Δ=0.05mm),對(duì)離散的點(diǎn)進(jìn)行坐標(biāo)映射變換并重構(gòu)樣條曲線,獲得模具型腔截面線。建立截面線映射變換坐標(biāo)系,以工藝模型葉根截面最大厚度圓圓心為坐標(biāo)原點(diǎn),葉根截面法矢為Z軸,X軸與葉根截面線前后緣頭曲線中心連線平行,由葉片截面的進(jìn)氣邊指向排氣邊,Y軸由笛卡兒坐標(biāo)系定則確定。上模型腔映射變換公式為式(6),下模型腔映射變換公式為式(7),映射變換如圖5所示,其中:S1為工藝模型第1條截面線;Sk為工藝模型第k條截面線;Sk+1為工藝模型第k+1條截面線;S′1模腔第1條截面線;S′k模腔第k 條截面線;S′k+1模腔第k+1條截面線。

    式中:(x,y,z)為工藝模型截面線族上離散點(diǎn)坐標(biāo);(x′,y′,z′)為模腔截面線族上點(diǎn)坐標(biāo)。

    2.3 模腔建模

    葉片輥軋模具型腔不同于鑄造、鍛造等模具型腔建模方法,其模腔位于軋輥扇面上,因此映射變換獲取模腔截面線,并通過(guò)截面線族建立型腔面片體,進(jìn)而建立模具型腔體,并設(shè)計(jì)定位特征。輥軋模具安裝軸心距D為136mm,模具腔體與模架相連,模架安裝在軋輥軸上。

    1)模腔曲面建模

    對(duì)工藝模型截面線族分別進(jìn)行幾何空間映射變換和基于補(bǔ)償?shù)挠成渥儞Q,建立輥軋模腔截面線族,并通過(guò)型腔截面線族生成型腔曲面,模腔建模如圖6所示。從圖中可以看出,針對(duì)工藝模型截面線族,進(jìn)行幾何空間映射變換,模腔對(duì)應(yīng)中心角θ′為42.3°;而在求解截面前滑值的基礎(chǔ)上,進(jìn)行前滑補(bǔ)償?shù)挠成渥儞Q,模腔對(duì)應(yīng)中心角θ為39.0°。

    2)輥軋模具模架設(shè)計(jì)

    為進(jìn)行輥軋成形分析,需要建立模具CAD模型。將X軸向上向下分別偏置10mm,長(zhǎng)度為葉片工藝模型寬度。將閉合的偏置曲線分別繞軋輥軸線旋轉(zhuǎn)模腔中心角θ,建立扇形柱體,使用型腔面片體分別對(duì)應(yīng)分割建立上下模腔體?;谀>吲c軋輥裝備關(guān)系,建立模架并設(shè)計(jì)定位,輥軋模具安裝。建立模架體截面,模架內(nèi)環(huán)半徑為軋輥半徑,外半徑與模腔體相連,模架寬度為模具型腔寬度兩倍,模架中心角120°,如圖7所示。

    3 補(bǔ)償方法驗(yàn)證與討論

    根據(jù)模具型腔模型和毛坯模型,建立有限元分析模型并進(jìn)行葉片輥軋成形仿真分析,提取輥軋成形的葉片模型并分析輥軋成形誤差,校驗(yàn)?zāi)>咝颓粠缀尉忍匦?。隨著數(shù)值計(jì)算方法的成熟,科學(xué)計(jì)算已與理論分析、實(shí)驗(yàn)驗(yàn)證并列成為科學(xué)研究的三大支柱。已有研究表明,輥軋成形過(guò)程的數(shù)值計(jì)算結(jié)果具有很高的精度和置信度,可實(shí)現(xiàn)葉片成形過(guò)程仿真及接觸應(yīng)力分析[14,22-24],建立模具與毛坯、毛坯內(nèi)部力和熱載荷的作用過(guò)程,其結(jié)果用于指導(dǎo)生產(chǎn),縮減物理實(shí)驗(yàn)成本。為此,在已有研究成果的基礎(chǔ)上,采用高置信度數(shù)值計(jì)算方法對(duì)葉片輥軋成形過(guò)程進(jìn)行分析及驗(yàn)證。

    3.1 有限元分析與工藝參數(shù)

    基于有限元軟件DEFORM-3D,建立輥軋分析模型并進(jìn)行輥軋成形分析。輥軋成形過(guò)程中,模具彈性變形較小,因此采用剛性體。所選葉片為某型航空發(fā)動(dòng)機(jī)壓氣機(jī)某級(jí)靜子葉片,其材料為GH4169合金。材料本構(gòu)模型使用Jonson-Cook模型[25],其參數(shù)如表1所示,其中:A 為室溫條件參考應(yīng)變力下屈服極限;B為溫度的函數(shù);C為應(yīng)變率對(duì)屈服極限的影響系數(shù);n為溫度變化率對(duì)屈服極限的影響系數(shù);m為溫度對(duì)屈服極限的影響系數(shù)。根據(jù)調(diào)研葉片輥軋成形參數(shù)[21],上下模腔分別繞各自軋輥軸線逆向旋轉(zhuǎn),角速度ω為0.5rad/s;為更準(zhǔn)確的模擬葉片輥軋成形,同時(shí)降低輥軋壓力和變形能耗,軋輥沿著輥軋方向做平動(dòng),速度與模腔旋轉(zhuǎn)線速度相同,葉片平均 厚 度 為 1.76mm,軋 輥 平 動(dòng) 速 度 v-為0.033 65m/s;為防止模腔與毛坯劇烈摩擦,采用潤(rùn)滑油潤(rùn)滑,毛坯與模腔切向接觸摩擦系數(shù)設(shè)為0.1,熱傳導(dǎo)系數(shù)為50W/(m·K);毛坯采用壓入進(jìn)給方式,毛坯末端面沿著輥軋方向載荷qH為500MPa。

    表1 材料本構(gòu)模型參數(shù)Table 1 Parameters of material constitutive model

    根據(jù)本文設(shè)計(jì)方法建立壓氣機(jī)葉片模具型腔和模架,壓氣機(jī)葉片輥軋成形分析模型如圖8(a)所示。通過(guò)有限元網(wǎng)格劃分工具對(duì)毛坯進(jìn)行三角面片體網(wǎng)格劃分,由于葉片弦向厚度變化劇烈,前后緣頭厚度不超過(guò)0.4mm,因此劃分需要細(xì)密,然而網(wǎng)格數(shù)目增大將增加運(yùn)算時(shí)間和影響收斂性,因此本文設(shè)置有限元分析模型中所有毛坯網(wǎng)格最大尺寸為0.1mm。完成有限元分析建立的葉片輥軋成形模型如圖8(b)所示,在輥軋成形過(guò)程中,最大應(yīng)力出現(xiàn)在截面線扭曲嚴(yán)重的葉身與葉尖扭曲過(guò)渡區(qū)域,葉根與葉身過(guò)渡段應(yīng)力次之,這是因?yàn)槿~尖截面線最大厚度圓直徑1.48mm,小于葉根截面最大厚度圓直徑1.90mm,而毛坯截面線相同,因此靠近葉尖截面壓下率大,同時(shí)截面扭曲嚴(yán)重,該結(jié)果與輥軋塑性變形過(guò)程分析吻合。

    3.2 輥軋成形前滑和葉片積疊高度校核

    基于前滑計(jì)算模型和型腔映射法則,建立輥軋模具型腔。通過(guò)數(shù)值計(jì)算獲取型腔前滑補(bǔ)償模具和幾何映射變換型具成形的葉片毛坯并對(duì)其進(jìn)行精度分析。葉片CAD模型積疊高度為38.00mm,工藝模型積疊高度為50.00mm,測(cè)繪比較不同模腔條件下葉片成形積疊高度,前滑補(bǔ)償模腔輥軋成形的葉片毛坯高度為50.39mm,然而采用幾何變換模腔輥軋成形的葉片毛坯高度為54.16mm,如圖9所示。進(jìn)行前滑補(bǔ)償輥軋成形葉片積疊高度誤差為+0.78%,幾何映射變換模腔輥軋成形葉片積疊高度誤差為+8.3%,進(jìn)行前滑補(bǔ)償建立葉片輥軋模具型腔,葉片積疊軸向精度提高了近11倍。

    4 結(jié) 論

    1)基于輥軋成形原理和葉片模具型腔設(shè)計(jì)過(guò)程,建立工藝模型截面線積疊高度與模腔截面線中心角關(guān)系,提出輥軋前滑補(bǔ)償方法,即精確控制型腔中心角,能夠提高輥軋葉片積疊成形精度。

    2)基于前滑補(bǔ)償模型,建立工藝模型截面線點(diǎn)到型腔截面線點(diǎn)的映射法則,對(duì)葉盆葉背截面線族進(jìn)行基于前滑補(bǔ)償?shù)挠成渥儞Q,建立模腔截面線,進(jìn)而建立輥軋模具型腔及輥軋模架。

    3)基于補(bǔ)償?shù)哪承蛪簹鈾C(jī)葉片輥軋模具型腔中心角為39.0°,通過(guò)有限元方法校核葉片輥軋成形精度,成形葉片積疊高度誤差為0.78%,有效的保證了葉片積疊軸向精度。

    論文在模腔設(shè)計(jì)中針對(duì)葉片輥軋前滑引起的積疊軸偏差進(jìn)行了補(bǔ)償,葉片積疊軸方向精度得到了有效控制;但是在研究過(guò)程中發(fā)現(xiàn),葉片塑性回彈引起截面前后緣頭的位置誤差較大,截面成形精度有待提高。型腔優(yōu)化設(shè)計(jì)是前滑補(bǔ)償和回彈補(bǔ)償?shù)鸟詈?,前滑補(bǔ)償控制葉片積疊線高度,回彈補(bǔ)償控制截面形狀尺寸。為此,下一步工作擬圍繞輥軋葉片彎曲回彈引起的前后緣尺寸處位置精度不足等問(wèn)題展開(kāi)研究,以期實(shí)現(xiàn)葉片在弦向的成形精度控制。

    [1] 李榮斌,姚枚,劉文昌,等.冷軋對(duì)GH4169合金組織與性能的影響[J].金屬熱處理,2002,27(7):12-15.LI R B,YAO M,LIU W C,et al.Effects of cold rolling on microstructure and performance of GH4169alloy[J].Metal Heat Treatment,2002,27(7):12-15 (in Chinese).

    [2] 喬思佳,李深亮.某轉(zhuǎn)子葉片輥軋工藝[J].中國(guó)新技術(shù)新產(chǎn)品,2015,33(22):38.QIAO S J,LI S L.Roll rolling process for a rotor blade[J].New Technology & New Products of China,2015,33(22):38(in Chinese).

    [3] 王輝,吳寶海,李小強(qiáng).新一代商用航空發(fā)動(dòng)機(jī)葉片的先進(jìn)加工技術(shù)[J].航空制造技術(shù),2014,20(20):26-31.WANG H,WU B H,LI X Q.Advanced machining technology of new generation commercial aero engine blade[J].Aeronautical Manufacturing Technology,2014,20(20):26-31(in Chinese).

    [4] 趙升噸,趙承偉,邵中魁,等.現(xiàn)代葉片成形工藝的探討[J].機(jī)床與液壓,2012,40(21):167-170.ZHAO S D,ZHAO C W,SHAO Z K,et al.Discussion about the modern forming process of leaves[J].Machine Tool and Hydraulics,2012,40(21):167-170 (in Chinese).

    [5] HEDAYATI A,NAJAFIZADEH A,KERMANPUR A,et al.The effect of cold rolling regime on microstructure and mechanical properties of AISI 304L[J].Stainless Steel Journal of Materials Processing Technology,2010,210(8):1017-1022.

    [6] 馮瑩娟.鍛造-冷輥軋GH4169合金葉片組織性能研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2012:16-30.FENG Y J.Study of microstructure and mechanical property of forged and cold rolling GH4169alloy blade[D].Harbin:Harbin Institute of Technology,2012:16-30(in Chinese).

    [7] 于建民.葉片溫輥軋成型工藝與裝備研究[D].太原:中北大學(xué),2006:34-42.YU J M.The research on the technology of the warm rolling and equipment of blade[D].Taiyuan:North University of China,2006:34-42(in Chinese).

    [8] 裴潤(rùn)奇,張治民,李保成.葉片專用輥軋機(jī)軋輥的強(qiáng)度校核[J].特種成形,2006,2(3):73-75.PEI R Q,ZHANG Z M,LI B C.Checking the roller strength in rolling machine for vane[J].Special Forming,2006,2(3):73-75(in Chinese).

    [9] 于建民,張治民.葉片輥軋工藝的計(jì)算機(jī)模擬[J].鍛壓裝備與制造技術(shù),2005,40(3):833-836.YU J M,ZHANG Z M.The computer simulation of blade rolling technology[J].China Metal Forming Equipment &Manufacturing Technology,2005,40(3):833-836 (in Chinese).

    [10] 劉慧,王國(guó)棟,劉相華.軋制工藝參數(shù)對(duì)鋼板平面形狀的影響[J].東北大學(xué)學(xué)報(bào)(自然科學(xué)版),2004,25(7):674-677.LIU H,WANG G D,LIU X H.Effect of rolling parameters on plane shape of plates[J].Journal of Northeastern University(Natural Science),2004,25(7):674-677(in Chinese).

    [11] 董連超.變厚度軋制金屬塑性流動(dòng)規(guī)律[D].秦皇島:燕山大學(xué),2013:20-30.DONG L C.Metal flow law of longitudinally profiled flat steel[D].Qinhuangdao:Yanshan University,2013:20-30(in Chinese).

    [12] 周道.航空葉片冷輥軋過(guò)程仿真分析[D].沈陽(yáng):東北大學(xué),2010:41-69.ZHOU D.Simulation and analysis of blade cold roll forming process[D].Shenyang:Northeastern University,2010:41-69(in Chinese).

    [13] MAMALIS A G,JOHNSON W,HAWKYARD J B.Pressure distribution,roll force and torque in cold ring rolling[J].International Journal of Mechanical Sciences,1976,4(18):196-209.

    [14] SEDIGHI M,MAHMOODI M.An approach to simulate cold roll-forging of turbo-engine thin compressor blade[J].Aircraft Engineering and Aerospace Technology:An International Journal,2009,81(3):191-198.

    [15] 陳雅,周紅梅.冷軋成形整流葉片毛邊問(wèn)題分析及改進(jìn)方案[J].模具技術(shù),2014,5(12):23-30.CHEN Y,ZHOU H M.Analysis and improvement plan for the straightener blade burrs in cold rolling process[J].Mould Technology,2014,5(12):23-30(in Chinese).

    [16] 李文平,邊文德,聶紹珉,等.補(bǔ)償回彈的沖壓件模具設(shè)計(jì)方法[J].鍛壓技術(shù),2007,32(5):86-91.LI W P,BIAN W D,NIE S M,et al.Die design method for sheet metal stamping work-piece with compensating spring-back[J].Forging & Stamping Technology,2007,32(5):86-91(in Chinese).

    [17] 陳洪月,孟輝,毛君,等.模具安裝誤差對(duì)葉片輥軋精度的影響[J].熱加工工藝,2014,43(5):120-125.CHEN H Y,MENG H,MAO J,et al.Influence of die installation error on blade rolling accuracy[J].Hot Working Technology,2014,43(5):120-125(in Chinese).

    [18] 裴景東.壓氣機(jī)葉片輥軋模具型腔自動(dòng)化建模技術(shù)研究[D].西安:西北工業(yè)大學(xué),2014:33-60.PEI J D.Research on automatic modeling technology about the cavity of mould of compressor roll-forming blade[D].Xi’an:Northwestern Polytechnical University,2014:33-60(in Chinese).

    [19] 王淵彬,汪文虎,張艷,等.壓氣機(jī)葉片輥軋模具型腔快速建模技術(shù)[J].航空學(xué)報(bào),2014,35(11):3190-3198.WANG Y B,WANG W H,ZHAN Y,et al.Rapid modeling technology of rolling mold cavity of aero-engine com-pressor blade[J].Acta Aeronautica et Astronautica Sinica,2014,35(11):3190-3198(in Chinese).

    [20] 趙志業(yè).金屬塑性變形與輥軋理論[M].北京:冶金工業(yè)出版社,2012:264-269.ZHAO Z Y.Metal plastic deformation and rolling theory[M].Beijing:Metallurgical Industry Press,2012:264-269(in Chinese).

    [21] 靳淇超,汪文虎,蔣睿嵩,等.一種改進(jìn)的壓氣機(jī)葉片輥軋成形前滑計(jì)算模型[J].航空學(xué)報(bào),2016,37(10):3178-3185.JIN Q C,WANG W H,JIANG R S,et al.An improved calculation model for forward slip in rolling compressor blade[J].Acta Aeronautica et Astronautica Sinica,2016,37(10):3178-3185(in Chinese).

    [22] LU B,OU H,ARMSTRONG C G,et al.3Ddie shape optimization for net-shape forging of aero-foil blades[J].Materials and Design,2009,30(7):2490-2500.

    [23] OU H,LAN J,ARMSTRONG C G,et al.An FE simulation and optimization approach for the forging of aero-engine components[J].Journal of Materials Processing Technology,2004,151(1-3):208-216.

    [24] SEDIGHI M,MAHMOODI M.Pressure distribution in cold rolling of turbo-engine thin compressor blades[J].Materials and Manufacturing Processes,2012,27(4):401-405.

    [25] 王濤,陳國(guó)定,巨江濤.GH4169高溫合金高應(yīng)變率本構(gòu)關(guān)系試驗(yàn)研究[J].航空學(xué)報(bào),2013,34(4):946-953.WANG T,CHEN G D,JU J T.Experimental study of constitutive relationship of super-alloy GH4169under high strain rates[J].Acta Aeronautica et Astronautica Sinica,2013,34(4):946-953(in Chinese).

    Forward slip compensation method for designing rolling mold cavity of compressor blade

    JIN Qichao,WANG Wenhu,JIANG Ruisong*,ZHAO Dezhong,CUI Kang,XIONG Yifeng
    Key Laboratory of Contemporary Design and Integrated Manufacturing Technology,Ministry of Education,Northwestern Polytechnical University,Xi’an 710072,China

    Rolling cavity is the key of rolling allowance free blade,F(xiàn)orward slip occurs in rolling process of high-pressure compressor blade,which cause a deviation of formed blade along rolling direction,so forward slip compensation need to be considered in designing mold cavity to improve the stacking height precision of blade.Contraposing forward slip phenomenon in rolling blade,the paper proposes and researches the optimum designing method which based on forward slip compensation for mold cavity of rolling compressor blade.Firstly,based on the analysis of cause and forward slip representation of rolling blade,the forward slip compensation mechanism is studied and a compensation model,namely,relationship between roller’s central angle and stack height of section is proposed.After that,a set of sections which characterized the process mold is drawn,forward slip value of every section is calculated,and relationship between roller’s central angle and stack height of section is set up.Then,based on the adjusted central angle,sector mapping algorithm in space for section curves is proposed and section curves are mapped;the rolling cavity surface through mapped section curves is reconstituted.Finally,an effective numerical calculation method is used to compare the stacking height of blade formed by forward slip compensated cavity and geometric cavity.The results show that optimized design of mold cavity can accurately control the blade forming accuracy of compressor blade in stack axis direction.

    compressor blade;roll forming;forward slip compensation;cavity optimum design;stacking height

    2016-03-01;Revised:2016-03-21;Accepted:2016-04-08;Published online:2014-04-18 14:58

    URL:www.cnki.net/kcms/detail/11.1929.V.20160418.1458.004.html

    s:National Natural Science Foundation of China (51475374,51505387);the Fundamental Research Funds for the Central Universities(3102015ZY087)

    V263.1

    A

    1000-6893(2017)01-420173-10

    http:/hkxb.buaa.edu.cn hkxb@buaa.edu.cn

    10.7527/S1000-6893.2016.0117

    2016-03-01;退修日期:2016-03-21;錄用日期:2016-04-08;網(wǎng)絡(luò)出版時(shí)間:2016-04-18 14:58

    www.cnki.net/kcms/detail/11.1929.V.20160418.1458.004.html

    國(guó)家自然科學(xué)基金 (51475374,51505387);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金 (3102015ZY087)

    *通訊作者 .E-mail:jiangrs@nwpu.edu.cn

    靳淇超,汪文虎,蔣睿嵩,等.壓氣機(jī)葉片輥軋模具型腔前滑補(bǔ)償方法[J].航空學(xué)報(bào),2017,38(1):420173.JIN Q C,WANG W H,JIANG R S,et al.Forward slip compensation method for designing rolling mold cavity of compressor blade[J].Acta Aeronautica et Astronautica Sinica,2017,38(1):420173.

    (責(zé)任編輯:李世秋)

    *Corresponding author.E-mail:jiangrs@nwpu.edu.cn

    猜你喜歡
    模腔型腔軋輥
    梅鋼1780產(chǎn)線高速鋼軋輥使用維護(hù)
    注塑工藝參數(shù)對(duì)模腔壓力分布的影響
    注塑工藝參數(shù)對(duì)模腔壓力的影響研究
    模腔壓力曲線與制品質(zhì)量的關(guān)系研究
    汽車(chē)內(nèi)飾件組合型腔注塑模設(shè)計(jì)
    數(shù)控軋輥磨床主輔系統(tǒng)分離的探討
    基于STEP-NC型腔特征識(shí)別方法的研究
    基于Mastercam的復(fù)雜型腔加工方法及其參數(shù)研究
    軋輥用012Al鋼熱處理試驗(yàn)研究
    上海金屬(2014年1期)2014-12-18 06:52:03
    采用模腔壓力曲線的注塑成型實(shí)驗(yàn)
    国产一级毛片七仙女欲春2| 国产av在哪里看| 91av网一区二区| 久久久久久久久久黄片| 给我免费播放毛片高清在线观看| av黄色大香蕉| h日本视频在线播放| 日韩欧美精品v在线| 国产午夜精品久久久久久| 亚洲国产欧美一区二区综合| 国产精品久久电影中文字幕| 国产精品98久久久久久宅男小说| 日韩欧美一区二区三区在线观看| 亚洲精品久久国产高清桃花| 亚洲欧洲精品一区二区精品久久久| 成人高潮视频无遮挡免费网站| 最好的美女福利视频网| 久久久精品欧美日韩精品| 欧美乱妇无乱码| 日韩欧美国产在线观看| 日韩中文字幕欧美一区二区| 宅男免费午夜| 男插女下体视频免费在线播放| 欧美国产日韩亚洲一区| 激情在线观看视频在线高清| 99久久精品热视频| 嫁个100分男人电影在线观看| 国产久久久一区二区三区| 日韩国内少妇激情av| 人人妻人人澡欧美一区二区| 男女之事视频高清在线观看| 两人在一起打扑克的视频| 国产av麻豆久久久久久久| 亚洲人成网站在线播放欧美日韩| 久久精品国产清高在天天线| 青草久久国产| 高潮久久久久久久久久久不卡| 亚洲国产高清在线一区二区三| 国产高清激情床上av| 69av精品久久久久久| 老司机午夜福利在线观看视频| 国产熟女xx| 伊人久久大香线蕉亚洲五| 波多野结衣高清作品| 中国美女看黄片| 婷婷丁香在线五月| 亚洲专区国产一区二区| 亚洲av免费在线观看| 成人无遮挡网站| 成人18禁在线播放| 亚洲熟女毛片儿| 丰满的人妻完整版| tocl精华| 中文字幕高清在线视频| 久久人人精品亚洲av| 99国产综合亚洲精品| 日本撒尿小便嘘嘘汇集6| 99国产精品99久久久久| 成人三级黄色视频| 国产成人精品无人区| 老司机福利观看| 亚洲最大成人中文| 国产欧美日韩一区二区三| 色精品久久人妻99蜜桃| 真实男女啪啪啪动态图| 久久久久久人人人人人| 亚洲九九香蕉| 久久久国产成人精品二区| 欧美乱码精品一区二区三区| 午夜影院日韩av| 亚洲九九香蕉| 国产av一区在线观看免费| h日本视频在线播放| 男女做爰动态图高潮gif福利片| 男女视频在线观看网站免费| 精品久久久久久久久久免费视频| 两个人的视频大全免费| 床上黄色一级片| 国产欧美日韩一区二区精品| 9191精品国产免费久久| 巨乳人妻的诱惑在线观看| av视频在线观看入口| 欧美成人免费av一区二区三区| www日本在线高清视频| 亚洲成人中文字幕在线播放| 久久久久国内视频| 免费在线观看日本一区| 久久久久久九九精品二区国产| 国产三级黄色录像| 国产成人欧美在线观看| 日韩精品中文字幕看吧| 亚洲,欧美精品.| 精品一区二区三区视频在线 | 欧美高清成人免费视频www| 成人国产一区最新在线观看| 国产成人福利小说| 国产成人av教育| 亚洲国产日韩欧美精品在线观看 | 黄色视频,在线免费观看| e午夜精品久久久久久久| 香蕉久久夜色| 1024香蕉在线观看| 国产麻豆成人av免费视频| 亚洲国产日韩欧美精品在线观看 | 两个人视频免费观看高清| 精品99又大又爽又粗少妇毛片 | 黄片大片在线免费观看| 国产成人精品久久二区二区免费| 亚洲国产精品sss在线观看| 亚洲欧美日韩卡通动漫| 精品一区二区三区av网在线观看| 亚洲精品在线观看二区| 精品国产乱码久久久久久男人| 精品久久蜜臀av无| 亚洲国产高清在线一区二区三| 亚洲精品美女久久久久99蜜臀| 国产精品美女特级片免费视频播放器 | 一级黄色大片毛片| 免费高清视频大片| 久久精品国产99精品国产亚洲性色| 天天一区二区日本电影三级| 深夜精品福利| 最新在线观看一区二区三区| 夜夜看夜夜爽夜夜摸| 麻豆国产97在线/欧美| 在线免费观看不下载黄p国产 | 久久久久九九精品影院| 中文字幕av在线有码专区| 2021天堂中文幕一二区在线观| 亚洲精品中文字幕一二三四区| 精品99又大又爽又粗少妇毛片 | 99久久精品国产亚洲精品| 无限看片的www在线观看| 国产又色又爽无遮挡免费看| 亚洲色图 男人天堂 中文字幕| 免费一级毛片在线播放高清视频| 色噜噜av男人的天堂激情| 国模一区二区三区四区视频 | 91av网站免费观看| 法律面前人人平等表现在哪些方面| 欧美日韩综合久久久久久 | 国产精品自产拍在线观看55亚洲| 精品免费久久久久久久清纯| 中国美女看黄片| 久久性视频一级片| 在线观看免费午夜福利视频| 在线免费观看不下载黄p国产 | 精品一区二区三区视频在线 | 欧美日韩福利视频一区二区| 国产高清videossex| 日本熟妇午夜| 国产av在哪里看| 成人特级av手机在线观看| 精品欧美国产一区二区三| 日韩欧美国产在线观看| 欧美激情久久久久久爽电影| 男女午夜视频在线观看| 国产午夜精品久久久久久| 成年人黄色毛片网站| 国产av一区在线观看免费| 老司机在亚洲福利影院| 搡老岳熟女国产| 国产精品自产拍在线观看55亚洲| 1024手机看黄色片| 欧美精品啪啪一区二区三区| 国产三级黄色录像| 99久久成人亚洲精品观看| ponron亚洲| 午夜亚洲福利在线播放| 国产探花在线观看一区二区| 日本五十路高清| 69av精品久久久久久| 国产在线精品亚洲第一网站| 亚洲国产精品999在线| 麻豆av在线久日| 精品国内亚洲2022精品成人| 91av网一区二区| 日韩高清综合在线| 国产av麻豆久久久久久久| 精品久久久久久,| 国产成人影院久久av| 小蜜桃在线观看免费完整版高清| 久久中文看片网| 国产成人一区二区三区免费视频网站| 黑人巨大精品欧美一区二区mp4| 国产午夜精品久久久久久| 国产v大片淫在线免费观看| 精品国产亚洲在线| 亚洲欧美一区二区三区黑人| 亚洲 欧美一区二区三区| 在线视频色国产色| 熟女人妻精品中文字幕| 亚洲专区国产一区二区| 午夜亚洲福利在线播放| 日韩大尺度精品在线看网址| 午夜a级毛片| 在线观看午夜福利视频| 国产精品99久久久久久久久| 国产亚洲欧美在线一区二区| 日本三级黄在线观看| 成人av一区二区三区在线看| 亚洲精华国产精华精| 波多野结衣巨乳人妻| avwww免费| 欧美色视频一区免费| 国产成+人综合+亚洲专区| 国产私拍福利视频在线观看| 狂野欧美激情性xxxx| 国产精品久久久久久久电影 | av片东京热男人的天堂| 久久精品91蜜桃| 偷拍熟女少妇极品色| 中文在线观看免费www的网站| 校园春色视频在线观看| 在线观看66精品国产| 国产精品爽爽va在线观看网站| 欧美中文综合在线视频| 国产av在哪里看| 性色avwww在线观看| 搡老熟女国产l中国老女人| 男女下面进入的视频免费午夜| 美女cb高潮喷水在线观看 | 日韩高清综合在线| 国产精品国产高清国产av| 亚洲av成人av| 中文在线观看免费www的网站| 2021天堂中文幕一二区在线观| 高清毛片免费观看视频网站| 国产精品乱码一区二三区的特点| 久久中文看片网| 女人被狂操c到高潮| 不卡一级毛片| 欧美黑人巨大hd| 看黄色毛片网站| 国产私拍福利视频在线观看| 国产激情久久老熟女| 日本免费一区二区三区高清不卡| 两性夫妻黄色片| 欧美日韩瑟瑟在线播放| 久久中文字幕一级| 看片在线看免费视频| 国产免费av片在线观看野外av| 一进一出抽搐动态| 亚洲一区二区三区色噜噜| 男女之事视频高清在线观看| 亚洲av片天天在线观看| 日韩三级视频一区二区三区| 国产人伦9x9x在线观看| 啦啦啦免费观看视频1| 99精品久久久久人妻精品| 日韩精品中文字幕看吧| 免费看美女性在线毛片视频| 在线国产一区二区在线| 国产三级黄色录像| 怎么达到女性高潮| 亚洲激情在线av| 免费看日本二区| 12—13女人毛片做爰片一| 99久久无色码亚洲精品果冻| 欧美绝顶高潮抽搐喷水| a级毛片a级免费在线| 黑人欧美特级aaaaaa片| 国产伦精品一区二区三区视频9 | xxx96com| 国产精品av视频在线免费观看| 亚洲人成伊人成综合网2020| 日韩欧美 国产精品| 亚洲五月天丁香| 中文资源天堂在线| 国产亚洲欧美在线一区二区| 啦啦啦韩国在线观看视频| 日韩欧美一区二区三区在线观看| 老司机午夜福利在线观看视频| 18禁观看日本| 久9热在线精品视频| 国产精品乱码一区二三区的特点| 国产精品av久久久久免费| 欧美一级a爱片免费观看看| 久久精品aⅴ一区二区三区四区| 国产日本99.免费观看| 人妻夜夜爽99麻豆av| 色精品久久人妻99蜜桃| 欧美日韩精品网址| 91在线精品国自产拍蜜月 | 国产精品影院久久| 别揉我奶头~嗯~啊~动态视频| 日本成人三级电影网站| 嫩草影院入口| 国产精品一区二区精品视频观看| 日本a在线网址| 99国产精品99久久久久| 亚洲精品粉嫩美女一区| 久久午夜综合久久蜜桃| 久久精品国产清高在天天线| 精品久久久久久久久久久久久| 无遮挡黄片免费观看| bbb黄色大片| 精品国产亚洲在线| 久久性视频一级片| 成年女人看的毛片在线观看| 久久久久九九精品影院| 熟女电影av网| 18禁裸乳无遮挡免费网站照片| 成人国产综合亚洲| 国产黄色小视频在线观看| 成人精品一区二区免费| 国产亚洲欧美在线一区二区| 丁香六月欧美| 床上黄色一级片| 国产高清三级在线| 久久久久久久久久黄片| 亚洲人成网站高清观看| 中文字幕人妻丝袜一区二区| 91老司机精品| 亚洲熟妇中文字幕五十中出| 亚洲在线自拍视频| 亚洲国产欧洲综合997久久,| 国产成人av激情在线播放| 久久精品aⅴ一区二区三区四区| 狂野欧美激情性xxxx| 国产av在哪里看| 亚洲人与动物交配视频| 日韩精品青青久久久久久| 麻豆成人午夜福利视频| 亚洲电影在线观看av| 欧美黑人巨大hd| 岛国在线观看网站| 亚洲国产精品成人综合色| 搞女人的毛片| 亚洲欧美一区二区三区黑人| 757午夜福利合集在线观看| 精品国产乱码久久久久久男人| 国产成人精品久久二区二区免费| 亚洲av电影在线进入| 美女被艹到高潮喷水动态| 亚洲精品色激情综合| 久久伊人香网站| 一级a爱片免费观看的视频| 日本撒尿小便嘘嘘汇集6| www.自偷自拍.com| 女警被强在线播放| 日韩精品青青久久久久久| 国产精品美女特级片免费视频播放器 | 国产伦人伦偷精品视频| 亚洲自偷自拍图片 自拍| 欧美性猛交╳xxx乱大交人| 一级毛片高清免费大全| 国产一区二区三区在线臀色熟女| 国产熟女xx| 国产精品自产拍在线观看55亚洲| 欧美日韩国产亚洲二区| 国产精品女同一区二区软件 | а√天堂www在线а√下载| 欧美乱码精品一区二区三区| 制服丝袜大香蕉在线| 国产伦一二天堂av在线观看| 亚洲七黄色美女视频| 国产伦一二天堂av在线观看| а√天堂www在线а√下载| 一个人看视频在线观看www免费 | 欧美成人免费av一区二区三区| 免费看a级黄色片| 两人在一起打扑克的视频| av中文乱码字幕在线| 97超级碰碰碰精品色视频在线观看| 国产精品久久久久久亚洲av鲁大| 成年免费大片在线观看| 18禁黄网站禁片免费观看直播| 欧美丝袜亚洲另类 | 免费看日本二区| 在线视频色国产色| 欧美日韩福利视频一区二区| 午夜视频精品福利| 美女 人体艺术 gogo| 国产乱人视频| 精品一区二区三区四区五区乱码| 午夜免费成人在线视频| 色视频www国产| 亚洲成人精品中文字幕电影| 老司机午夜福利在线观看视频| 美女大奶头视频| 国产精品亚洲av一区麻豆| 日本与韩国留学比较| aaaaa片日本免费| 国产三级在线视频| 精品久久蜜臀av无| 久久这里只有精品中国| 俺也久久电影网| 国产视频一区二区在线看| a在线观看视频网站| 夜夜看夜夜爽夜夜摸| 香蕉国产在线看| 悠悠久久av| 午夜精品久久久久久毛片777| 午夜视频精品福利| 欧美色欧美亚洲另类二区| 熟妇人妻久久中文字幕3abv| 免费高清视频大片| 中文字幕av在线有码专区| 国产一区二区三区在线臀色熟女| 欧美绝顶高潮抽搐喷水| a级毛片在线看网站| 亚洲欧美精品综合一区二区三区| 久久性视频一级片| 男人舔女人的私密视频| 国产97色在线日韩免费| 最新中文字幕久久久久 | 欧美日本亚洲视频在线播放| 99久久综合精品五月天人人| 免费av毛片视频| 99热这里只有是精品50| 日韩三级视频一区二区三区| 久久精品国产综合久久久| 国产精品av久久久久免费| 99国产精品一区二区三区| 国产高清有码在线观看视频| 黑人巨大精品欧美一区二区mp4| 桃红色精品国产亚洲av| 午夜福利高清视频| 日本黄色片子视频| 一a级毛片在线观看| 日本三级黄在线观看| 他把我摸到了高潮在线观看| 免费在线观看日本一区| 国产精品精品国产色婷婷| 老司机深夜福利视频在线观看| 色播亚洲综合网| 久久香蕉精品热| 精品一区二区三区视频在线 | 最近视频中文字幕2019在线8| 特级一级黄色大片| 51午夜福利影视在线观看| 久久精品亚洲精品国产色婷小说| cao死你这个sao货| 国产精品久久久久久久电影 | 国产1区2区3区精品| 99国产综合亚洲精品| 757午夜福利合集在线观看| 久久精品91无色码中文字幕| 免费在线观看影片大全网站| 日本撒尿小便嘘嘘汇集6| www日本黄色视频网| 国产av一区在线观看免费| 亚洲18禁久久av| 一区二区三区国产精品乱码| 日本免费a在线| 亚洲av五月六月丁香网| 亚洲五月婷婷丁香| 嫩草影视91久久| 亚洲精品456在线播放app | 日日夜夜操网爽| 成人特级av手机在线观看| 国产免费男女视频| 麻豆av在线久日| 黄色 视频免费看| 午夜激情福利司机影院| 午夜福利欧美成人| 制服人妻中文乱码| 久久久久久大精品| 老熟妇乱子伦视频在线观看| 男人和女人高潮做爰伦理| 国产精品久久久久久亚洲av鲁大| 亚洲av电影不卡..在线观看| 精品乱码久久久久久99久播| 日韩欧美在线乱码| 成人午夜高清在线视频| 日韩中文字幕欧美一区二区| 精品无人区乱码1区二区| 精华霜和精华液先用哪个| 欧美国产日韩亚洲一区| 麻豆国产97在线/欧美| 一本综合久久免费| 一边摸一边抽搐一进一小说| 91字幕亚洲| 国产v大片淫在线免费观看| 午夜免费观看网址| 看黄色毛片网站| 最近最新中文字幕大全免费视频| 久久精品91蜜桃| 午夜激情福利司机影院| 精品国产三级普通话版| 亚洲 欧美一区二区三区| 我要搜黄色片| 亚洲熟女毛片儿| 热99re8久久精品国产| 啦啦啦免费观看视频1| 国产aⅴ精品一区二区三区波| 国产精品九九99| 国产精品乱码一区二三区的特点| 免费观看的影片在线观看| 人人妻,人人澡人人爽秒播| 婷婷丁香在线五月| 亚洲国产欧洲综合997久久,| 又黄又粗又硬又大视频| 男插女下体视频免费在线播放| www.自偷自拍.com| 欧洲精品卡2卡3卡4卡5卡区| 久久香蕉国产精品| 日韩精品中文字幕看吧| 国产欧美日韩一区二区精品| 美女免费视频网站| 淫妇啪啪啪对白视频| 亚洲狠狠婷婷综合久久图片| 精品福利观看| 国产成人aa在线观看| 日韩 欧美 亚洲 中文字幕| 精品一区二区三区av网在线观看| 免费看美女性在线毛片视频| 全区人妻精品视频| 两人在一起打扑克的视频| 在线永久观看黄色视频| 亚洲av成人一区二区三| 午夜福利在线观看吧| 久久久久久人人人人人| 精品99又大又爽又粗少妇毛片 | 日韩国内少妇激情av| 又大又爽又粗| 亚洲真实伦在线观看| 欧美国产日韩亚洲一区| 午夜福利免费观看在线| 变态另类成人亚洲欧美熟女| 非洲黑人性xxxx精品又粗又长| 热99在线观看视频| 99国产精品99久久久久| 亚洲五月天丁香| 欧美+亚洲+日韩+国产| www.精华液| 日韩三级视频一区二区三区| 国产精品一区二区三区四区免费观看 | aaaaa片日本免费| 丰满的人妻完整版| 黑人操中国人逼视频| 久久精品91无色码中文字幕| 欧美又色又爽又黄视频| 午夜久久久久精精品| 听说在线观看完整版免费高清| 免费av不卡在线播放| 一进一出抽搐动态| 九九久久精品国产亚洲av麻豆 | 欧美黑人欧美精品刺激| 国产免费男女视频| 欧美色视频一区免费| 国产一区在线观看成人免费| 看免费av毛片| 亚洲国产精品成人综合色| 日本免费一区二区三区高清不卡| 午夜成年电影在线免费观看| 欧美日韩中文字幕国产精品一区二区三区| 黄色日韩在线| 欧美色欧美亚洲另类二区| 美女黄网站色视频| 国产午夜精品论理片| 毛片女人毛片| 美女午夜性视频免费| 美女被艹到高潮喷水动态| 啪啪无遮挡十八禁网站| 波多野结衣高清无吗| 色综合欧美亚洲国产小说| 九色成人免费人妻av| 免费无遮挡裸体视频| 操出白浆在线播放| 日韩欧美一区二区三区在线观看| 在线观看舔阴道视频| 国产成人影院久久av| 久久久国产欧美日韩av| 久久久精品大字幕| 91在线精品国自产拍蜜月 | 露出奶头的视频| 三级男女做爰猛烈吃奶摸视频| 免费在线观看亚洲国产| 岛国视频午夜一区免费看| 亚洲欧美日韩东京热| 一级作爱视频免费观看| 婷婷六月久久综合丁香| 久久精品国产综合久久久| 青草久久国产| 午夜福利欧美成人| 成人精品一区二区免费| 十八禁人妻一区二区| 亚洲国产精品999在线| 成人欧美大片| 成人无遮挡网站| 精品久久久久久,| 亚洲 欧美一区二区三区| 日本 欧美在线| 国产三级在线视频| 中文字幕久久专区| 亚洲黑人精品在线| av天堂在线播放| 天堂网av新在线| 老鸭窝网址在线观看| 国产男靠女视频免费网站| 淫妇啪啪啪对白视频| 免费人成视频x8x8入口观看| 99国产极品粉嫩在线观看| 亚洲在线自拍视频| 夜夜爽天天搞| 亚洲美女黄片视频| 亚洲美女视频黄频| 亚洲成av人片在线播放无| 久久久色成人| 国产伦精品一区二区三区视频9 | 女同久久另类99精品国产91| www国产在线视频色| 国产一区二区在线av高清观看| 小说图片视频综合网站| 亚洲第一欧美日韩一区二区三区| 九色成人免费人妻av| 99国产精品一区二区三区| 狂野欧美白嫩少妇大欣赏| 97超视频在线观看视频| 手机成人av网站| 国产97色在线日韩免费| 国产黄a三级三级三级人| 在线观看66精品国产|