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

    基于力約束的空心渦輪葉片陶芯定位方法

    2017-11-20 01:20:43崔康汪文虎蔣睿嵩趙德中靳淇超
    航空學報 2017年9期
    關(guān)鍵詞:位姿布局元件

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

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

    基于力約束的空心渦輪葉片陶芯定位方法

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

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

    精鑄蠟型作為空心渦輪葉片精鑄過程重要的前期工藝轉(zhuǎn)接件,其壁厚精度主要由蠟型模具型腔與內(nèi)部陶芯的位置匹配關(guān)系決定。由于陶芯在模具內(nèi)完全依靠定位元件實現(xiàn)空間定位,為減小由定位誤差引起的陶芯位姿漂移,提出了一種基于力平衡約束的空心渦輪葉片精鑄模具陶芯定位布局優(yōu)化方法。首先,通過建立陶芯定位誤差傳遞模型,揭示了定位誤差與陶芯空間位姿擾動量之間的映射關(guān)系;其次,根據(jù)力平衡原理構(gòu)建了基于力約束的陶芯定位布局優(yōu)化模型;之后,針對陶芯表面定位候選點的離散分布特性,結(jié)合遺傳算法給出了陶芯定位布局點的詳細求解策略。最后,仿真對比證明了利用本文所提方法獲得的陶芯定位方案可以在保證陶芯定位穩(wěn)定性的同時提高陶芯定位精度,此外,按照優(yōu)化后的定位方案壓制實際蠟型,壁厚檢測結(jié)果也進一步表明所提方法的有效性。

    空心渦輪葉片; 蠟型; 陶芯; 定位布局; 力平衡; 遺傳算法

    空心渦輪葉片是航空發(fā)動機的核心熱端部件,由于結(jié)構(gòu)復(fù)雜、精度及強度要求高,目前主要通過熔模精密鑄造方法制備,具體流程包括:制芯,模具設(shè)計、壓蠟、制殼、脫蠟、澆鑄等環(huán)節(jié)。精鑄蠟型作為重要前期工藝轉(zhuǎn)接件,其形狀精度是影響葉片成形質(zhì)量的重要因素[1]??招臏u輪葉片精鑄蠟型形狀精度主要包含外形輪廓精度以及壁厚精度。研究表明:通過提升蠟料性能[2]、預(yù)測型面誤差[3-5]、優(yōu)化鑄件結(jié)構(gòu)[6]及相關(guān)工藝參數(shù)[7-8]可以有效改善蠟型外形輪廓精度;而蠟型壁厚精度主要由蠟型模具型腔與內(nèi)部陶芯的位置匹配關(guān)系決定[9]。在蠟型模具中,陶芯完全依靠定位及夾緊元件限位固定,在模具型腔不變的情況下,由定位誤差引起的陶芯位姿漂移是導致蠟型壁厚偏差的決定性因素。因此,通過合理設(shè)計陶芯定位布局方案,減小定位誤差對陶芯定位精度的影響,可以有效提高精鑄蠟型壁厚精度。

    針對空心渦輪葉片精鑄蠟型模具陶芯定位問題,崔康等[10]提出了一種通過逆向調(diào)整定位元件尺寸進而提高蠟型壁厚精度的方法。另外,基于模型匹配算法,馮煒等[11]提出了一種陶芯定位元件尺寸計算方法。然而,上述文獻并未對陶芯定位布局方案設(shè)計展開探討。事實上,陶芯定位可近似為一般工件定位問題。Asnte[12]基于彈性接觸力學相關(guān)理論以及有限元仿真技術(shù)給出了工件加工過程中各定位夾緊元件載荷計算方法。Wu和Chan[13]以及Chou等[14]以定位夾緊元件載荷為優(yōu)化目標,分別基于遺傳算法和線性規(guī)劃理論給出了工件定位夾緊布局的相關(guān)計算方法。Liao[15]針對板類零件變形量提出了一種板類零件定位夾緊位置空間布局優(yōu)化方法。Li和Melkote[16]假定工件在加工過程中的彈性位移是由定位夾緊點處的工件彈性變形引起的,并基于此給出了這一種工件定位夾緊布局優(yōu)化方法。Kaya[17]、Prabhaharan等[18]以及Padmanaban等[19-20]將工件加工過程中表面型值點變形量作為優(yōu)化目標,分別對呈離散或連續(xù)分布狀態(tài)的定位夾緊候選點集給出了基于有限元理論及啟發(fā)式算法的定位夾緊布局優(yōu)化方法。Rex和Ravindran[21]則遍歷工件所有可能的定位夾緊方案,通過對加工過程進行數(shù)值仿真,得到了最優(yōu)定位夾緊方案。然而,上述方法均是以減小定位夾緊元件受力或工件變形量為主要優(yōu)化目標,并未考慮定位誤差對工件定位精度的影響。Wang等[22-23]和Xiong等[24]提出一種基于定位誤差模型的工件定位布局優(yōu)化方法,但由于沒有考慮工件自重對定位方案的影響,因此,優(yōu)化結(jié)果常無法滿足工件定位穩(wěn)定性要求;此外,由于該模型采用Greed算法和Interchange算法求解,因此,計算結(jié)果也常陷入局部最優(yōu)。

    為此,本文針對空心渦輪葉片精鑄蠟型模具陶芯定位問題,提出了一種基于力平衡約束的定位布局求解方法。首先通過建立定位誤差傳遞模型,獲得了陶芯定位誤差與空間位姿擾動量之間的映射關(guān)系。其次,考慮到陶芯自重問題,建立了基于力平衡約束的陶芯定位布局優(yōu)化模型,并結(jié)合遺傳算法給出了該模型的詳細求解步驟。最后,結(jié)合具體實例對本文及Wang等[22-23]所提方法進行了對比驗證,結(jié)果表明本文所提方法在陶芯定位穩(wěn)定性以及定位誤差控制方面具有一定優(yōu)勢,可為空心渦輪葉片蠟型壁厚控制提供基礎(chǔ)理論。

    1 定位誤差傳遞模型

    如圖1所示,在空心渦輪葉片蠟型壓制過程中,陶芯完全依靠定位及夾緊元件實現(xiàn)空間限位,定位元件決定了陶芯空間姿態(tài),夾緊元件保證了陶芯與定位元件始終接觸。由于定位元件存在定位誤差,陶芯實際位姿相較于理論位置必然存在一定擾動。為研究定位布局方案對陶芯定位精度的影響,首先需要建立定位誤差與陶芯空間位姿擾動量之間的函數(shù)關(guān)系。

    圖1 陶芯定位示意圖Fig.1 Illustration of ceramic core localization

    根據(jù)剛體運動學知識可知,剛體在三維空間o-xyz內(nèi)的任意運動可以分解為沿某方向λ的平移運動及繞某方向ω的旋轉(zhuǎn)運動。假設(shè)由定位誤差引起的陶芯空間位姿擾動量為

    Δq=[ΔbTΔθT]T

    (1)

    式中:Δb=[ΔbxΔbyΔbz]T為平移擾動量;Δθ=[ΔθxΔθyΔθz]T=Θ·τ為旋轉(zhuǎn)擾動量,且τ為沿旋轉(zhuǎn)軸ω方向的單位向量,Θ為轉(zhuǎn)動角。如圖2所示,陶芯表面上任意一點i的位置偏移量Δri也必然包含平移分量Δb以及另外一個繞軸ω的旋轉(zhuǎn)分量Δφi,即

    Δri=Δb+Δφi

    (2)

    此外,由圖2可知,當旋轉(zhuǎn)擾動量Δθ較小時,Δφi的表達式為

    Δφi≈(Θdi)·ci=Θ·(τ×ri)=Δθ×ri

    (3)

    式中:di為點i到旋轉(zhuǎn)軸ω的歐式距離;Θdi為旋轉(zhuǎn)弧長;ci為沿旋轉(zhuǎn)方向且與由旋轉(zhuǎn)軸ω及點i所

    圖2 陶芯定位誤差傳遞模型Fig.2 Model for locating error transfer of ceramic core

    構(gòu)成的平面垂直的單位向量;ri為點i的坐標向量。

    由此,聯(lián)立式(2)和式(3)可以得陶芯空間位姿擾動量Δq與陶芯表面點i位置偏移量Δri之間的映射關(guān)系為

    Δri=Δb+Δθ×ri

    (4)

    另外,由于點i坐標ri滿足陶芯型面方程:

    f(ri)=0

    (5)

    因此,對式(5)微分可以進一步求得陶芯表面點i處指向陶芯內(nèi)部的單位法向矢量為

    (6)

    事實上,若點i恰好是陶芯表面某定位點,且將偏移量Δri在點i法矢方向的投影近似為定位元件在該點處的定位誤差,那么,陶芯定位誤差與陶芯空間位姿擾動量之間的映射關(guān)系可簡化為

    (7)

    式中:ji=[(dri)T(ri×dri)T]T。

    當陶芯表面存在k個定位點,且定位誤差分別為Δe1、Δe2、Δe3、…、Δek,代入式(7)聯(lián)立可得定位誤差與陶芯空間位姿擾動量之間的函數(shù)關(guān)系:

    (8)

    式中:J為定位Jacobian矩陣。

    另外,考慮工件定位完整性,當陶芯表面各定位點不存在定位誤差:ΔE=0,陶芯的空間位姿應(yīng)唯一存在,即Δq=0應(yīng)是式(8)的唯一解。由此可知:陶芯表面至少需要6個定位點,且對應(yīng)定位Jacobian矩陣J滿秩,才能保證陶芯空間位姿唯一不變。

    2 定位布局優(yōu)化模型

    根據(jù)式(8)可知,陶芯空間位姿擾動量Δq不僅取決于定位誤差,同時還與定位Jacobian矩陣J有關(guān),前者是由定位元件制造誤差及裝配誤差造成,不可避免;后者取決于定位點在陶芯表面的空間分布。因此,在保證陶芯定位完整性的前提下,尋求一種最佳定位布局方案使得相同定位誤差條件下陶芯位姿擾動量最小是空心渦輪葉片精鑄蠟型模具陶芯定位布局設(shè)計的主要目標。

    2.1 定位布局求解基本模型

    對式(8)作如下變換:

    ΔE=JΔq?

    (9)

    同時令:

    M=JTJ

    (10)

    綜上所述,陶芯定位布局基本求解模型可以表述為:在陶芯表面的候選點集P_set中找出使接觸信息矩陣行列式det(M)最大的6個定位點,且該行列式最大值det(M)max≠0,即

    (11)

    2.2 重力平衡約束條件

    陶芯定位元件不僅決定了陶芯的空間位姿,其還應(yīng)該滿足陶芯定位穩(wěn)定性需求,即在定位點處由陶芯自重引起的定位元件反作用力應(yīng)與陶芯自重滿足靜力平衡關(guān)系,從而防止陶芯置放時滑移、滑落。

    陶芯在蠟型模具中一般采用葉盆朝下或者葉背朝下兩種置放形式,如圖3所示。將兩種置放形式下的陶芯重力載荷分別表示為Gbasin、Gback,且

    (12)

    圖3 陶芯定位方式Fig.3 Locating modes of ceramic core

    式中:g為葉盆朝下時的重力矢量;t為葉盆朝下時的重力矩矢量。在不考慮陶芯與定位元件之間摩擦力的情況下,上述兩種置放形式對應(yīng)的陶芯定位靜力平衡方程可表示為

    (13)

    (14)

    在陶芯定位過程中,定位元件只能提供沿接觸點法矢方向且指向陶芯內(nèi)部的正壓力。針對上述兩種置放形式,如果在任意一種置放形式下,若某種陶芯定位布局方案能使得定位點平衡力大于零,即正壓力,那么可以認為該定位布局方案滿足陶芯定位穩(wěn)定性需求。因此,基于重力平衡的陶芯定位布局求解約束條件可表示為

    (15)

    據(jù)此,基于重力約束的空心渦輪葉片精鑄蠟型模具陶芯定位布局求解問題可以表示為

    (16)

    為了便于求解,現(xiàn)將式(16)模型作出如下簡化:

    1) 針對約束條件一,如果優(yōu)化結(jié)果det(M)max=0,那么本次候選點集中不存在滿足陶芯定位唯一性的定位布局方案,因此需要更新候選點集。但是在模型優(yōu)化求解過程中,不需要考慮約束條件一對優(yōu)化結(jié)果的影響。

    2) 針對約束條件二,可以采用內(nèi)部懲罰法,將約束條件整合到優(yōu)化模型中。最終得到無約束的空心渦輪葉片蠟型模具陶芯定位布局求解模型:

    6 points→from:P_set

    (17)

    式中:γ、η為力平衡約束因子。

    3 模型優(yōu)化算法

    由于陶芯表面定位候選點集呈離散分布狀態(tài),且候選點與優(yōu)化模型不存在明顯的梯度關(guān)系,因此針對連續(xù)參數(shù)空間的傳統(tǒng)優(yōu)化算法并不適用于上述模型的求解。目前,針對離散空間的模型優(yōu)化問題,一般采用啟發(fā)式算法求解。遺傳算法就是一種針對非連續(xù)參數(shù)空間,通過模擬生物進化機制從而實現(xiàn)高效、并行、全局搜索的啟發(fā)式優(yōu)化算法。該方法不依賴優(yōu)化模型的梯度信息,采用概率化的尋優(yōu)方式,能夠跳出局部最優(yōu)解并自適應(yīng)的調(diào)整搜索方向,非常適合本文所提陶芯定位布局優(yōu)化模型。因此,本文基于遺傳算法,建立了空心渦輪葉片精鑄蠟型模具陶芯定位布局優(yōu)化模型求解方法,具體步驟如下:

    步驟1 候選點選擇。在陶芯表面適合定位的部位選取N個候選點,對其進行二進制編碼。

    步驟2 初始化。選取m種布局方案作為初始種群,且m為偶數(shù)。每種方案被稱之為種群個體,每個種群個體包含了6個從候選點集中隨機挑選的定位點。

    步驟3 個體評價。按照式(18)對初始種群進行個體評價,計算每個個體的適應(yīng)度Ti。

    i=1,2,…,m

    (18)

    步驟4 選擇運算。按照式(19)計算初值種群每個個體被挑選到的概率分布pi,并利用隨機數(shù)在初始種群中挑選m/2對個體。

    (19)

    步驟5 交叉及變異運算。基于每個個體在候選點中的二進制編碼,對步驟4中m/2對個體進行交叉及變異置換,并按照式(18)重新對交叉變異后的m個個體的適應(yīng)度進行求解。

    步驟6 種群迭代。在初始種群以及交叉變異后的種群中分別選取適應(yīng)度最大的2個和m-2個個體作為下次迭代的初始種群,并跳轉(zhuǎn)至步驟3。

    步驟7 停止條件。當?shù)螖?shù)超過預(yù)設(shè)最高迭代次數(shù)n時,計算終止。

    4 算法驗證

    結(jié)合具體空心渦輪葉片精鑄蠟型模具陶芯定位布局問題,對上述算法及Wang等[22-23]所提算法進行了對比分析,并通過壓制實際蠟型對優(yōu)化結(jié)果進行了實例驗證,分析及驗證過程如下。

    4.1 對比計算

    圖4為驗證用空心渦輪葉片精鑄蠟型陶芯三維模型。由陶芯定位原理可知,在蠟型壓制完成后蠟型表面會殘留定位孔,為保證蠟型葉身型面質(zhì)量,定位元件候選點應(yīng)分布在葉尖以及葉根位置,同時綜合考慮蠟型模具結(jié)構(gòu)及陶芯裝夾工藝特點,最終在陶芯表面選取了256個定位候選點,如圖5所示。

    另外,基于陶芯型面方程,獲取了各候選點空間坐標以及指向陶芯內(nèi)部的單位法向矢量,并將陶芯定位布局優(yōu)化模型及求解算法中的相關(guān)參數(shù)設(shè)置如表1所示。

    圖4 實驗用陶芯三維模型Fig.4 3D model of experimental ceramic core

    最后,利用遺傳算法對目標方程式(17)進行了迭代求解,結(jié)果表明每代種群的個體適應(yīng)度最大值Tmax呈收斂狀態(tài),如圖6所示。最終得到的陶芯定位布局方案如圖7所示。圖8為基于Wang等[22-23]所提方法求解得到的陶芯定位布局方案。

    為驗證上述方案是否滿足陶芯定位穩(wěn)定性要求,根據(jù)式(14)分別計算了兩種方案下為平衡陶芯重力各定位點所需力的大小,結(jié)果如表2所示。

    圖5 陶芯定位候選點Fig.5 Locating candidate points of experimental ceramic core

    表1 陶芯定位布局優(yōu)化模型相關(guān)參數(shù)值

    Table 1 Parameter values involved in optimization model for layout of ceramic core localization

    ParameterγηmnValue11100500

    圖6 各代種群適應(yīng)度最大值Fig.6 Maximum individual fitness in each generation

    圖7 陶芯定位方案1Fig.7 Locating strategy 1

    圖8 陶芯定位方案2Fig.8 Locating strategy 2

    表2 兩種定位方案各定位元件陶芯重力反作用力大小Table 2 Reacting forces required from locating rods to balance core gravity with respect to two locating strategies

    LocatingstrategyForcerequiredfromlocation/NRod1Rod2Rod3Rod4Rod5Rod6100.00150.6398 0.01670.2929 0.5507200.18620.7117-0.04930.8375-0.0237

    在方案1中,各定位點所需要提供的力均大于0,即:所有定位元件只需提供正壓力便可滿足陶芯靜力平衡。相反,在方案2中,定位元件需要在定位點4、6處提供拉力,明顯不滿足陶芯定位穩(wěn)定性需求。由此可見,本文所提算法能夠有效約束定位元件布局方案使其滿足陶芯靜力平衡需求。優(yōu)化后的方案可以在置放陶芯時為陶芯提供足夠支撐,防止其滑移或滑落。

    另外,為驗證本文所提方法在陶芯定位誤差控制方面的性能,本文設(shè)計如下對比實驗:

    步驟1 由于陶芯定位誤差具有統(tǒng)計意義,即各個定位點誤差均服從獨立的正態(tài)分布。為此,基于正態(tài)分布N(0,0.05)為上述兩種定位方案分別隨機生成了1 000組定位誤差。

    步驟2 根據(jù)式(8)分別計算每組定位誤差引起的陶芯空間位姿擾動量Δq。

    步驟3 陶芯空間位姿擾動量包含平移擾動量以及旋轉(zhuǎn)擾動量,由于平移和旋轉(zhuǎn)具有累加效應(yīng),因此無法直觀的通過位姿擾動量對陶芯位姿變化進行判定。為此在陶芯表面對蠟型壁厚影響最為敏感前后緣位置挑選了100個評估點(如圖9 所示),分別計算陶芯漂移前后位姿評估點的空間位移量,并將最大值作為陶芯位姿變化的衡量標準。

    計算結(jié)果如圖10和圖11所示。在定位方案1中,1 000組定位誤差引起的陶芯位姿評估點最大

    圖9 陶芯漂移評估點Fig.9 Evaluating points for shift of ceramic core

    圖10 基于定位方案1的陶芯最大偏移量 Fig.10 Maximum displacement of ceramic core with respect to locating strategy 1

    圖11 基于定位方案2的陶芯最大偏移量 Fig.11 Maximum displacement of ceramic core with respect to locating strategy 2

    位移量的平均值為0.068 7 mm。在定位方案2中,該值為0.086 3 mm。上述結(jié)果說明,相較于Wang等[22-23]的方法,本文所提算法在保證陶芯定位穩(wěn)定性的同時在陶芯定位誤差控制方面也具有一定優(yōu)越性。

    4.2 實例驗證

    按照圖7所示的陶芯定位布局方案及陶芯置放形式(葉盆朝下),作者設(shè)計了相應(yīng)的定位元件以及蠟型模具,如圖12所示。最終,基于表3所示的工藝參數(shù)壓制了實際蠟型5組,如圖13所示。對關(guān)鍵截面進行壁厚測量,結(jié)果表明蠟型壁厚滿足設(shè)計公差需求(±0.05 mm),如表4所示。該實驗證明了本文所提方法具有一定實際可行性。

    圖12 蠟型模具及壓蠟過程Fig.12 Wax pattern die and wax injection process

    表3 壓蠟工藝參數(shù)Table 3 Parameter values of wax injection experiment

    ParameterInjectionpressure/MPaHoldingpressure/MPaWaxtempera-ture/℃Holdingtime/sCoolingtime/sValue1.01.0602020

    圖13 精鑄蠟型Fig.13 Wax patterns

    表4 精鑄蠟型測量截面最大壁厚偏差

    Table 4 Maximum wall-thickness errors for measurement sections of wax patterns

    Numberofmea-surementsectionError/mmPattern1Pattern2Pattern3Pattern4Pattern51 0.048 0.039-0.026 0.039-0.0312-0.036-0.041 0.037-0.048 0.0383 0.042 0.047 0.046-0.041-0.039

    5 結(jié) 論

    1) 基于剛體定位誤差傳遞理論以及靜力平衡原理,提出了一種空心渦輪葉片精鑄蠟型模具陶芯定位布局優(yōu)化模型,并結(jié)合遺傳算法給出了該模型的求解方法。

    2) 仿真對比分析發(fā)現(xiàn):相較于Wang等[22-23]所提的方法,利用本文所提方法得到的陶芯定位布局方案可以在保證陶芯定位力學穩(wěn)定性的同時進一步減小由陶芯定位誤差引起的陶芯位姿漂移。

    3) 壓蠟實驗表明:利用本文所提方法設(shè)計的陶芯定位方案可以保證壓制出的蠟型壁厚精度滿足設(shè)計公差要求(小于±0.05 mm)。

    4) 值得指出的是:壓蠟工藝參數(shù)對蠟型壁厚精度也會產(chǎn)生一定影響,而本文未就工藝參數(shù)的選擇及優(yōu)化展開討論,該方面的研究將在后續(xù)論文中進行深入探討。

    [1] PETES F, VOIGT R, BLAIR M. Dimensional repeatability of investment castings[C]//9th World Conference on Investment Casting. Montvale, NJ: Investment Casting Institute, 1996: 22.

    [2] BEMBLAGE O, KARUNAKAR D B. A Study on the blended wax patterns in investment casting process[C]//Proceedings of the World Congress on engineering. The International Association of Engineers, 2011, 1: 6-8.

    [3] SINGH B, KUMAR P, MISHRA B K. Simulation of wax pattern dimensions for accuracy improvement in ceramic shell investment casting[J]. International Journal of Surface Engineering & Materials Technology, 2013, 3(1): 45-50.

    [4] SABAU A S, VISWANATHAN S. Prediction of wax pattern dimensions in investment casting[J]. Transactions-American Foundrymens Society, 2002, 1: 733-746.

    [5] SABAU A S, VISWANATHAN S. Material properties for predicting wax pattern dimensions in investment casting[J]. Materials Science and Engineering: A, 2003, 362(1): 125-134.

    [6] LIU C, JIN S, LAI X, et al. Influence of complex structure on the shrinkage of part in investment casting process[J]. The International Journal of Advanced Manufacturing Technology, 2015, 77(5-8): 1191-1203.

    [7] SINGH B, KUMAR P, MISHRA B K. Simulation of wax pattern dimensions for accuracy improvement in ceramic shell investment casting[J]. International Journal of Surface Engineering and Materials Technology, 2013, 3(1): 45-50.

    [8] PATTNAIK S, KARUNAKAR D B, JHA P K. Multi-characteristic optimization of wax patterns in the investment casting process using grey-fuzzy logic[J]. The International Journal of Advanced Manufacturing Technology, 2013, 67(5-8): 1577-1587.

    [9] JIANG R S, WANG W H, ZHANG D H, et al. Wall thickness monitoring method for wax pattern of hollow turbine blade[J]. The International Journal of Advanced Manufacturing Technology, 2016, 83(5): 949-960.

    [10] 崔康, 汪文虎, 蔣睿嵩, 等. 渦輪葉片精鑄模具陶芯定位元件逆向調(diào)整算法[J]. 航空學報, 2011, 32(10): 1924-1929.

    CUI K, WANG W H, JIANG R S, et al. Reverse adjustment algorithm of ceramic core locators in hollow turbine blade investment casting die[J]. Acta Aeronautica et Astronautica Sinica, 2011, 32(10): 1924-1929 (in Chinese).

    [11] 馮煒, 汪文虎, 王孝忠, 等. 空心渦輪葉片精鑄蠟型陶芯定位元件尺寸計算方法[J]. 航空學報, 2013, 34(1): 181-186.

    FENG W, WANG W H, WANG X Z, et al. Size calculation method of ceramic core locators for hollow turbine blade investment casting wax pattern[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(1): 181-186 (in Chinese).

    [12] ASNTE J N. A combined contact elasticity and finite element-based model for contact load and pressure distribution calculation in a frictional workpiece-fixture system[J]. The International Journal of Advanced Manufacturing Technology, 2008, 39(5-6): 578-588.

    [13] WU N H, CHAN K C. A genetic algorithm based approach to optimal fixture configuration[J]. Computers & Industrial Engineering, 1996, 31(3): 919-924.

    [14] CHOU Y C, CHANDRU V, BARASH M M. A mathematical approach to automatic configuration of machining fixtures: analysis and synthesis[J]. Journal of Engineering for Industry, 1989, 111(4): 299-306.

    [15] LIAO Y G. A genetic algorithm-based fixture locating positions and clamping schemes optimization[J]. Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, 2003, 217(8): 1075-1083.

    [16] LI B, MELKOTE S N. Improved workpiece location accuracy through fixture layout optimization[J]. International Journal of Machine Tools and Manufacture, 1999, 39(6): 871-883.

    [17] KAYA N. Machining fixture locating and clamping position optimization using genetic algorithms[J]. Computers in Industry, 2006, 57(2): 112-120.

    [18] PRABHAHARAN G, PADMANABAN K P, KRISHNAKUMAR R. Machining fixture layout optimization using FEM and evolutionary techniques[J]. The International Journal of Advanced Manufacturing Technology, 2007, 32(11-12): 1090-1103.

    [19] PADMANABAN K P, PRABHAHARAN G. Dynamic analysis on optimal placement of fixturing elements using evolutionary techniques[J]. International Journal of Production Research, 2008, 46(15): 4177-4214.

    [20] PADMANABAN K P, ARULSHRI K P, PRABHAHARAN G. Machining fixture layout design using ant colony algorithm based continuous optimization method[J]. The International Journal of Advanced Manufacturing Technology, 2009, 45(9-10): 922-934.

    [21] REX F M T, RAVINDRAN D. An integrated approach for optimal fixture layout design[J]. Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, 2017, 231(7): 1217-1228.

    [22] WANG M Y, PELINESCU D M. Optimizing fixture layout in a point-set domain[J]. IEEE Transactions on Robotics and Automation, 2001, 17(3): 312-323.

    [23] WANG M Y. An optimum design for 3-D fixture synthesis in a point set domain[J]. IEEE Transactions on Robotics and Automation, 2000, 16(6): 839-846.

    [24] XIONG Z, WANG M Y, LI Z. A near-optimal probing strategy for workpiece localization[J]. IEEE Transactions on Robotics, 2004, 20(4): 668-676.

    [25] ATKINSON A C, DONEV A N, TOBIAS R D. Optimum experimental designs, with SAS[M]. Oxford: Oxford University Press, 2007: 137-147.

    (責任編輯: 李世秋)

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

    Force-constraint method for localization of ceramic core of hollow turbine blade

    CUI Kang, WANG Wenhu, JIANG Ruisong*, ZHAO Dezhong, JIN Qichao

    KeyLaboratoryofContemporaryDesignandIntegratedManufacturingTechnologyofMinistryofEducation,NorthwesternPolytechnicalUniversity,Xi’an710072,China

    The wax pattern is always used as a dimension transfer component in near-net-shape casting process for a hollow turbine blade, and its wall-thickness accuracy entirely depends on the positional relationship between the die cavity of the wax pattern and the internal ceramic core. Generally, the ceramic core is located in the wax pattern die through a series of locating rods. In order to reduce the positional shift of the ceramic core caused by locating errors, a locating layout optimization method based on the force-balance constraint is proposed in this paper. An error transfer model, which formulates the mapping relationship between localization errors and perpetuation of the ceramic core, is established. According to the static equilibrium theory, an optimization model for locating the layout of the ceramic core is then proposed based on gravity constraint. Considering the discrete feature of locating the candidate point on the surface of the ceramic core, a solving strategy for the optimization model is given by utilizing the genetic algorithm. Comparisons of simulation results prove that the locating layout optimized with the method in this paper can improve the localization accuracy of the ceramic core, while guaranteeing the localization stability. Based on a wax injection experiment, the feasibility of the optimization result is also demonstrated.

    hollow turbine blade; wax pattern; ceramic core; localization layout; force balance; genetic algorithm

    2017-03-01; Revised: 2017-04-13; Accepted: 2017-04-21; Published online: 2017-05-19 09:27

    URL: www.cnki.net/kcms/detail/11.1929.V.20170519.0927.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)09-421209-09

    2017-03-01; 退修日期: 2017-04-13; 錄用日期: 2017-04-21; 網(wǎng)絡(luò)出版時間: 2017-05-19 09:27

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

    國家自然科學基金(51475374, 51505387); 中央高?;究蒲袠I(yè)務(wù)費專項資金(3102015ZY087)

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

    崔康, 汪文虎, 蔣睿嵩, 等. 基于力約束的空心渦輪葉片陶芯定位方法[J]. 航空學報, 2017, 38(9): 421209. CUI K, WANG W H, JIANG R S, et al. Force-constraint method for localization of ceramic core of hollow turbine blade[J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(9): 421209.

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

    10.7527/S1000-6893.2017.421209

    猜你喜歡
    位姿布局元件
    BP的可再生能源布局
    能源(2017年5期)2017-07-06 09:25:57
    基于共面直線迭代加權(quán)最小二乘的相機位姿估計
    基于CAD模型的單目六自由度位姿測量
    VR布局
    QFN元件的返工指南
    小型四旋翼飛行器位姿建模及其仿真
    2015 我們這樣布局在探索中尋找突破
    在新興產(chǎn)業(yè)看小元件如何發(fā)揮大作用
    寶馬i3高電壓元件介紹(上)
    Face++:布局刷臉生態(tài)
    伦精品一区二区三区| 少妇被粗大的猛进出69影院| 看免费成人av毛片| 制服诱惑二区| 91在线精品国自产拍蜜月| 一级a爱视频在线免费观看| 日本91视频免费播放| 国产欧美亚洲国产| 欧美日韩视频精品一区| 久久青草综合色| 黄网站色视频无遮挡免费观看| 美女中出高潮动态图| 中文字幕另类日韩欧美亚洲嫩草| 国产精品亚洲av一区麻豆 | 丰满少妇做爰视频| 熟妇人妻不卡中文字幕| 欧美最新免费一区二区三区| 国产黄频视频在线观看| 国产综合精华液| 久久人人爽人人片av| 国产精品久久久av美女十八| 午夜福利网站1000一区二区三区| 亚洲综合色网址| 亚洲综合精品二区| 亚洲美女搞黄在线观看| 精品国产一区二区三区久久久樱花| 视频区图区小说| 久久97久久精品| 免费看av在线观看网站| 高清视频免费观看一区二区| 少妇被粗大猛烈的视频| 少妇猛男粗大的猛烈进出视频| 老汉色∧v一级毛片| 日韩在线高清观看一区二区三区| 如日韩欧美国产精品一区二区三区| 国产高清国产精品国产三级| 亚洲av综合色区一区| 国产伦理片在线播放av一区| 国产熟女午夜一区二区三区| 精品少妇内射三级| 丁香六月天网| 成人亚洲精品一区在线观看| 亚洲人成77777在线视频| 亚洲精品av麻豆狂野| 美女国产视频在线观看| 国产麻豆69| 亚洲成国产人片在线观看| 99九九在线精品视频| 飞空精品影院首页| 国产男女超爽视频在线观看| 国产精品人妻久久久影院| 中文字幕精品免费在线观看视频| 成年女人在线观看亚洲视频| 纵有疾风起免费观看全集完整版| 亚洲欧洲精品一区二区精品久久久 | freevideosex欧美| 观看美女的网站| 亚洲成人一二三区av| 中文字幕人妻丝袜一区二区 | 久久人人爽av亚洲精品天堂| 久久久久久久大尺度免费视频| 国产探花极品一区二区| 亚洲av男天堂| a 毛片基地| 久热这里只有精品99| 最近最新中文字幕大全免费视频 | 极品人妻少妇av视频| 精品少妇内射三级| 日日啪夜夜爽| 两个人免费观看高清视频| 啦啦啦中文免费视频观看日本| 丰满乱子伦码专区| 午夜91福利影院| 在现免费观看毛片| 最近手机中文字幕大全| 亚洲欧美清纯卡通| 亚洲精品中文字幕在线视频| 国产成人精品一,二区| 尾随美女入室| 精品酒店卫生间| 亚洲av日韩在线播放| 成人国产av品久久久| 国产成人免费无遮挡视频| 亚洲美女搞黄在线观看| 国产福利在线免费观看视频| 久久热在线av| av线在线观看网站| 天堂中文最新版在线下载| 成年女人毛片免费观看观看9 | 亚洲欧洲日产国产| 热re99久久精品国产66热6| 人人妻人人爽人人添夜夜欢视频| 国产熟女午夜一区二区三区| 午夜福利视频精品| 国产精品香港三级国产av潘金莲 | 欧美激情高清一区二区三区 | 丝袜脚勾引网站| 人妻 亚洲 视频| 午夜日本视频在线| 精品少妇久久久久久888优播| 久热久热在线精品观看| 五月伊人婷婷丁香| www.精华液| 韩国av在线不卡| 国产精品不卡视频一区二区| 日韩精品有码人妻一区| 欧美 亚洲 国产 日韩一| 日韩视频在线欧美| 免费高清在线观看视频在线观看| 啦啦啦视频在线资源免费观看| 乱人伦中国视频| 熟女电影av网| 亚洲av福利一区| 欧美精品亚洲一区二区| 久久精品国产亚洲av涩爱| 少妇被粗大猛烈的视频| 免费日韩欧美在线观看| h视频一区二区三区| 丝袜人妻中文字幕| 91精品伊人久久大香线蕉| 新久久久久国产一级毛片| 欧美xxⅹ黑人| 久久久精品94久久精品| 午夜精品国产一区二区电影| 波野结衣二区三区在线| 国产精品 欧美亚洲| 亚洲三区欧美一区| 亚洲av电影在线进入| 亚洲精品乱久久久久久| 国产精品免费视频内射| 性高湖久久久久久久久免费观看| 欧美精品国产亚洲| 亚洲三区欧美一区| 国产有黄有色有爽视频| 成人影院久久| 女人高潮潮喷娇喘18禁视频| 国产极品天堂在线| 精品久久久精品久久久| 日韩在线高清观看一区二区三区| 99热网站在线观看| 在线免费观看不下载黄p国产| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 丰满少妇做爰视频| 人妻 亚洲 视频| 国产 一区精品| 在线观看免费高清a一片| 免费黄频网站在线观看国产| 男女边摸边吃奶| 色视频在线一区二区三区| 18禁观看日本| 国产成人精品无人区| 成人18禁高潮啪啪吃奶动态图| 亚洲精品国产av蜜桃| 香蕉丝袜av| 黑人猛操日本美女一级片| 黄色 视频免费看| 国产精品av久久久久免费| 精品午夜福利在线看| 亚洲三区欧美一区| 免费观看av网站的网址| 国产视频首页在线观看| 老熟女久久久| 亚洲精品自拍成人| videossex国产| 久久婷婷青草| 精品一区二区三区四区五区乱码 | 午夜久久久在线观看| 热re99久久精品国产66热6| 老司机影院毛片| 国产精品免费大片| 天天躁夜夜躁狠狠躁躁| 亚洲av日韩在线播放| 亚洲情色 制服丝袜| 日韩制服丝袜自拍偷拍| 午夜福利网站1000一区二区三区| 美女福利国产在线| 日韩在线高清观看一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品第二区| 黄色配什么色好看| 一级毛片电影观看| 久久久久久人妻| 久久精品国产亚洲av高清一级| 制服人妻中文乱码| 韩国精品一区二区三区| 久久精品亚洲av国产电影网| 老司机影院成人| 午夜福利,免费看| 国产欧美日韩一区二区三区在线| 一级片免费观看大全| 欧美日韩视频高清一区二区三区二| 亚洲精品一区蜜桃| 国产黄频视频在线观看| 国产片特级美女逼逼视频| 99re6热这里在线精品视频| 国产97色在线日韩免费| 久热久热在线精品观看| 麻豆av在线久日| 久久韩国三级中文字幕| 国产精品不卡视频一区二区| 亚洲综合色惰| 色婷婷久久久亚洲欧美| 亚洲三级黄色毛片| 亚洲欧美精品自产自拍| 亚洲精品视频女| 久久这里只有精品19| 国产一区二区 视频在线| a 毛片基地| 久久精品国产a三级三级三级| 亚洲第一青青草原| 国产女主播在线喷水免费视频网站| 亚洲精品乱久久久久久| av不卡在线播放| 在线观看人妻少妇| 国产成人精品在线电影| 街头女战士在线观看网站| 大码成人一级视频| 日韩免费高清中文字幕av| 男人舔女人的私密视频| 欧美日韩综合久久久久久| 免费在线观看黄色视频的| 性色avwww在线观看| 亚洲,欧美精品.| 久久久久视频综合| 中文字幕制服av| 亚洲av欧美aⅴ国产| 黄色毛片三级朝国网站| 欧美黄色片欧美黄色片| 国产精品 欧美亚洲| 乱人伦中国视频| 一区二区三区精品91| 亚洲欧美清纯卡通| 90打野战视频偷拍视频| 久久久国产一区二区| 最近的中文字幕免费完整| 国产精品久久久av美女十八| 夜夜骑夜夜射夜夜干| 国产精品久久久久久精品电影小说| 在线观看人妻少妇| 免费观看在线日韩| 日韩av免费高清视频| 人妻系列 视频| 人人妻人人澡人人爽人人夜夜| 亚洲伊人久久精品综合| 人人澡人人妻人| 国产精品熟女久久久久浪| 国产成人aa在线观看| 国产乱人偷精品视频| 一区二区日韩欧美中文字幕| 久久久国产精品麻豆| 黄片无遮挡物在线观看| 精品第一国产精品| 亚洲欧美日韩另类电影网站| 日韩在线高清观看一区二区三区| 男人操女人黄网站| 午夜福利视频精品| 中文字幕制服av| 久久久国产精品麻豆| 欧美日韩一级在线毛片| 亚洲成人手机| 午夜影院在线不卡| 久久免费观看电影| 中国国产av一级| 亚洲五月色婷婷综合| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久精品国产自在天天线| 国产老妇伦熟女老妇高清| 精品一品国产午夜福利视频| av在线观看视频网站免费| 午夜日韩欧美国产| 亚洲av在线观看美女高潮| 中文天堂在线官网| 久久人人爽av亚洲精品天堂| 男女午夜视频在线观看| 母亲3免费完整高清在线观看 | 国产日韩欧美视频二区| 性高湖久久久久久久久免费观看| 精品久久蜜臀av无| 人妻人人澡人人爽人人| 在线亚洲精品国产二区图片欧美| 国产成人午夜福利电影在线观看| 一级片免费观看大全| 男女高潮啪啪啪动态图| 精品国产超薄肉色丝袜足j| 老汉色av国产亚洲站长工具| 亚洲精品久久久久久婷婷小说| 欧美日韩视频高清一区二区三区二| 亚洲欧美一区二区三区国产| 亚洲成av片中文字幕在线观看 | 啦啦啦视频在线资源免费观看| 九九爱精品视频在线观看| 黑丝袜美女国产一区| 99香蕉大伊视频| 精品卡一卡二卡四卡免费| 国产色婷婷99| 婷婷色麻豆天堂久久| av国产久精品久网站免费入址| 捣出白浆h1v1| 黄色怎么调成土黄色| 久久精品熟女亚洲av麻豆精品| 精品一区二区免费观看| 九九爱精品视频在线观看| 91aial.com中文字幕在线观看| 九色亚洲精品在线播放| 精品福利永久在线观看| 丝袜脚勾引网站| 久久人妻熟女aⅴ| 99久久综合免费| 99国产精品免费福利视频| 激情五月婷婷亚洲| 少妇的逼水好多| 高清av免费在线| 国产一区有黄有色的免费视频| 精品久久久久久电影网| av线在线观看网站| 又大又黄又爽视频免费| 久久久久久久国产电影| 美女主播在线视频| 如何舔出高潮| 欧美日韩精品成人综合77777| 日韩中文字幕欧美一区二区 | 欧美日韩一区二区视频在线观看视频在线| 天天躁日日躁夜夜躁夜夜| 久久精品人人爽人人爽视色| 大片免费播放器 马上看| 免费大片黄手机在线观看| 亚洲国产精品一区三区| 亚洲欧美成人精品一区二区| 制服人妻中文乱码| 欧美激情 高清一区二区三区| 国产激情久久老熟女| 国产免费又黄又爽又色| 一本色道久久久久久精品综合| 天天躁夜夜躁狠狠躁躁| 亚洲av成人精品一二三区| 欧美 日韩 精品 国产| 青草久久国产| 一级黄片播放器| 日韩av不卡免费在线播放| av女优亚洲男人天堂| 大话2 男鬼变身卡| 欧美日韩av久久| 成人毛片60女人毛片免费| 亚洲四区av| 欧美激情高清一区二区三区 | 国产黄色视频一区二区在线观看| 久久精品久久久久久久性| 最近2019中文字幕mv第一页| 精品人妻在线不人妻| 欧美亚洲 丝袜 人妻 在线| 黄色怎么调成土黄色| 亚洲第一av免费看| 制服丝袜香蕉在线| 免费在线观看完整版高清| 叶爱在线成人免费视频播放| 国产一区二区在线观看av| 9191精品国产免费久久| 久久精品熟女亚洲av麻豆精品| 哪个播放器可以免费观看大片| 99久久综合免费| 午夜福利视频在线观看免费| 视频在线观看一区二区三区| xxxhd国产人妻xxx| 新久久久久国产一级毛片| 大话2 男鬼变身卡| 久久综合国产亚洲精品| 久久久精品区二区三区| 亚洲人成77777在线视频| 色哟哟·www| 精品99又大又爽又粗少妇毛片| 午夜日本视频在线| 丰满饥渴人妻一区二区三| 久久99精品国语久久久| 国产成人欧美| 中文字幕人妻熟女乱码| 伊人久久大香线蕉亚洲五| 成人国产麻豆网| 韩国精品一区二区三区| 女人精品久久久久毛片| 伊人久久大香线蕉亚洲五| 丝袜人妻中文字幕| 99精国产麻豆久久婷婷| 久久久国产一区二区| 久久人人爽人人片av| 亚洲在久久综合| 一区二区av电影网| 欧美国产精品va在线观看不卡| 亚洲人成77777在线视频| videossex国产| 叶爱在线成人免费视频播放| 五月天丁香电影| 纵有疾风起免费观看全集完整版| 午夜激情av网站| 亚洲精品视频女| 国产乱来视频区| 亚洲国产最新在线播放| 狠狠婷婷综合久久久久久88av| 日本欧美国产在线视频| 日韩av不卡免费在线播放| 欧美成人午夜精品| 五月伊人婷婷丁香| 精品亚洲乱码少妇综合久久| 中文精品一卡2卡3卡4更新| 日韩视频在线欧美| 搡老乐熟女国产| 欧美+日韩+精品| 免费观看性生交大片5| 黄色配什么色好看| av电影中文网址| 亚洲第一区二区三区不卡| 亚洲精品日本国产第一区| 精品人妻偷拍中文字幕| 国产精品一二三区在线看| 午夜av观看不卡| 男人舔女人的私密视频| 少妇的丰满在线观看| 美女中出高潮动态图| 精品一区在线观看国产| 欧美日本中文国产一区发布| 日韩欧美精品免费久久| 久久精品国产亚洲av高清一级| 人妻少妇偷人精品九色| 久久精品亚洲av国产电影网| 熟妇人妻不卡中文字幕| 久久久久国产一级毛片高清牌| 欧美av亚洲av综合av国产av | 国产av国产精品国产| 国产白丝娇喘喷水9色精品| 一级毛片我不卡| 亚洲欧美一区二区三区黑人 | 最黄视频免费看| 国产爽快片一区二区三区| 国产一区二区在线观看av| 欧美精品高潮呻吟av久久| www.精华液| 尾随美女入室| 亚洲成人手机| 国产女主播在线喷水免费视频网站| 少妇人妻精品综合一区二区| 哪个播放器可以免费观看大片| 亚洲图色成人| 自拍欧美九色日韩亚洲蝌蚪91| 成人午夜精彩视频在线观看| 国产成人av激情在线播放| 国产白丝娇喘喷水9色精品| 夜夜骑夜夜射夜夜干| 欧美精品人与动牲交sv欧美| 婷婷成人精品国产| 日日撸夜夜添| 两个人免费观看高清视频| 91久久精品国产一区二区三区| 精品少妇黑人巨大在线播放| 欧美人与性动交α欧美精品济南到 | 国产精品香港三级国产av潘金莲 | 免费久久久久久久精品成人欧美视频| 欧美老熟妇乱子伦牲交| 在线观看免费高清a一片| 日本免费在线观看一区| 亚洲欧美日韩另类电影网站| 日韩伦理黄色片| 国产深夜福利视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 日日爽夜夜爽网站| 黄频高清免费视频| 秋霞在线观看毛片| 免费播放大片免费观看视频在线观看| 巨乳人妻的诱惑在线观看| 搡老乐熟女国产| 国产一区二区激情短视频 | 久久99蜜桃精品久久| 自线自在国产av| 国产精品人妻久久久影院| 国产成人午夜福利电影在线观看| 久久 成人 亚洲| 免费观看在线日韩| 青春草国产在线视频| 26uuu在线亚洲综合色| 免费黄频网站在线观看国产| 搡老乐熟女国产| 女性生殖器流出的白浆| 久久久久精品性色| 哪个播放器可以免费观看大片| 涩涩av久久男人的天堂| 国产视频首页在线观看| 久久精品人人爽人人爽视色| 最近手机中文字幕大全| 在线天堂最新版资源| 好男人视频免费观看在线| 中文字幕人妻丝袜一区二区 | 97在线视频观看| 一级爰片在线观看| 99热网站在线观看| 久久精品国产a三级三级三级| 久久精品久久久久久噜噜老黄| 国产成人欧美| 男人操女人黄网站| 麻豆精品久久久久久蜜桃| 最新的欧美精品一区二区| 黄色视频在线播放观看不卡| 纵有疾风起免费观看全集完整版| 国产精品熟女久久久久浪| 男女免费视频国产| 女性被躁到高潮视频| 久久精品国产亚洲av高清一级| 高清欧美精品videossex| 香蕉丝袜av| 最近手机中文字幕大全| 欧美日韩一区二区视频在线观看视频在线| 天天躁日日躁夜夜躁夜夜| 99国产综合亚洲精品| 春色校园在线视频观看| 天天操日日干夜夜撸| 超碰97精品在线观看| 中国三级夫妇交换| 18禁国产床啪视频网站| 丰满少妇做爰视频| 大话2 男鬼变身卡| 亚洲国产色片| 亚洲成av片中文字幕在线观看 | 欧美日韩视频精品一区| 亚洲三区欧美一区| 各种免费的搞黄视频| 亚洲av福利一区| 精品国产一区二区三区四区第35| 女人精品久久久久毛片| 涩涩av久久男人的天堂| 又粗又硬又长又爽又黄的视频| 国产精品二区激情视频| 不卡av一区二区三区| 两性夫妻黄色片| 97在线人人人人妻| 亚洲精品国产色婷婷电影| 老司机影院毛片| 一边摸一边做爽爽视频免费| 亚洲av.av天堂| 亚洲精品一二三| 国语对白做爰xxxⅹ性视频网站| 18在线观看网站| av在线app专区| 亚洲国产精品一区三区| 免费观看在线日韩| 黄色 视频免费看| 一级片'在线观看视频| 黄色配什么色好看| 卡戴珊不雅视频在线播放| 亚洲,欧美,日韩| 在线天堂最新版资源| 日本av手机在线免费观看| 精品人妻熟女毛片av久久网站| 青春草亚洲视频在线观看| 欧美老熟妇乱子伦牲交| 久久午夜综合久久蜜桃| 尾随美女入室| 国产成人免费无遮挡视频| 母亲3免费完整高清在线观看 | 日本免费在线观看一区| 老汉色∧v一级毛片| 午夜激情久久久久久久| 亚洲中文av在线| 久久久久国产精品人妻一区二区| 极品人妻少妇av视频| 欧美激情 高清一区二区三区| 高清欧美精品videossex| 男人操女人黄网站| 一区福利在线观看| 亚洲婷婷狠狠爱综合网| 只有这里有精品99| 国产日韩欧美在线精品| 啦啦啦视频在线资源免费观看| 日韩人妻精品一区2区三区| 日日啪夜夜爽| 综合色丁香网| 中文天堂在线官网| 高清在线视频一区二区三区| 人人澡人人妻人| 国产成人精品久久二区二区91 | 一区福利在线观看| 亚洲国产欧美网| 久久99蜜桃精品久久| 美女午夜性视频免费| 啦啦啦视频在线资源免费观看| 亚洲精品自拍成人| 久久久a久久爽久久v久久| 日本av手机在线免费观看| 欧美精品人与动牲交sv欧美| 精品人妻一区二区三区麻豆| 观看美女的网站| 国产麻豆69| 欧美黄色片欧美黄色片| 成年人午夜在线观看视频| 国产国语露脸激情在线看| 宅男免费午夜| 五月天丁香电影| 亚洲国产av新网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 99re6热这里在线精品视频| a级毛片黄视频| 性少妇av在线| 老汉色av国产亚洲站长工具| 国产成人aa在线观看| 日韩不卡一区二区三区视频在线| 国产又色又爽无遮挡免| 99久国产av精品国产电影| 成人国产麻豆网| 国产成人一区二区在线| 丝袜脚勾引网站| 女人精品久久久久毛片| 成年人免费黄色播放视频| 一级黄片播放器| xxx大片免费视频| 久久精品人人爽人人爽视色| 国产亚洲一区二区精品| 日本猛色少妇xxxxx猛交久久| 国产av码专区亚洲av| 久热久热在线精品观看|