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

    在加熱階段滾塑模具內(nèi)表面?zhèn)鳠嵯禂?shù)的兩種研究方法

    2023-02-07 07:36:10劉學(xué)軍
    中國塑料 2023年1期
    關(guān)鍵詞:粉料傳熱系數(shù)百分比

    劉學(xué)軍

    (北京工商大學(xué)人工智能學(xué)院,北京100048)

    0 前言

    滾塑工藝的主要缺點是在加熱階段的能耗比較大[1],早期對滾塑模具是采用明火直燒的加熱方式,但后來出于節(jié)能環(huán)保及安全方面的考慮,大部分滾塑模具采用在烘箱內(nèi)通過熱空氣循環(huán)對流的方式來進(jìn)行加熱。但由于烘箱的壁面和門、泄壓孔等都會向環(huán)境散熱,實際上烘箱加熱方式的能量利用率仍然很低[2]。為了進(jìn)一步降低能耗,近些年來一些研究者提出了對模具采用油加熱和電加熱的方式[3?4]。本文所研究的滾塑模具是一個繞雙軸轉(zhuǎn)動的電加熱鋁制模具,嵌在該模具外表面上的加熱絲通電加熱模具,用來生產(chǎn)氫氣瓶。所研究的加熱階段是指從通電加熱開始到粉料即將開始熔融的這一過程。

    與采用傳統(tǒng)加熱方式的滾塑模具一樣,電加熱模具的內(nèi)表面也是通過強迫對流換熱將熱量傳給模內(nèi)的粉料和空氣,同時粉料一邊隨模具翻滾一邊與空氣混合,因此也會從空氣中吸收熱量。無論滾塑工藝采用哪一種加熱方式,模具內(nèi)表面向模內(nèi)混合物的傳熱都是一個復(fù)雜的氣固兩相流問題。為了強化該傳熱過程,都希望模具內(nèi)表面的傳熱系數(shù)是越大越好。要想定量地評價滾塑模具向模內(nèi)混合物傳熱的強弱程度,就必須準(zhǔn)確地計算出模具內(nèi)表面上的傳熱系數(shù)。另外在建立傳熱理論模型對滾塑工藝進(jìn)行仿真計算時,模具的內(nèi)表面?zhèn)鳠嵯禂?shù)通常是作為事先已知的輸入?yún)?shù),而傳熱系數(shù)的準(zhǔn)確性直接影響了仿真的精度。因此探索模具內(nèi)表面?zhèn)鳠嵯禂?shù)的研究方法在滾塑研究中具有重要的理論意義。

    以往很少有文獻(xiàn)對模具的內(nèi)表面?zhèn)鳠嵯禂?shù)進(jìn)行理論或?qū)嶒灧矫娴难芯?,在用到該傳熱系?shù)時都是按經(jīng)驗給出一個估計值,并沒有相應(yīng)的理論或?qū)嶒炓罁?jù)。Nugent等[5]測量了一個未裝粉料的鋼制模具的外表面溫度和模內(nèi)溫度,然后通過傳熱模型得出加熱過程中其內(nèi)表面?zhèn)鳠嵯禂?shù)是1 W/(m2?K),但未裝料的模具與實際的滾塑模具有本質(zhì)的差別。在文獻(xiàn)[6]?[9]中采用的模具內(nèi)表面?zhèn)鳠嵯禂?shù)都是5 W/(m2?K)。Banerjee等[10]認(rèn)為模具內(nèi)表面與粉料間的傳熱系數(shù)為25 W/(m2?K),與模內(nèi)空氣間的傳熱系數(shù)為 2 W/(m2?K)。Abdullah等[11]認(rèn)為模具與粉料接觸處的內(nèi)表面?zhèn)鳠嵯禂?shù)為30 W/(m2?K),與模內(nèi)空氣接觸處的內(nèi)表面?zhèn)鳠嵯禂?shù)為0.55 W/(m2?K)??梢钥闯?,現(xiàn)有的研究對模具內(nèi)表面?zhèn)鳠嵯禂?shù)的取值差異很大。這是因為該傳熱系數(shù)的影響因素比較多,也難以預(yù)測。它除了與模具的結(jié)構(gòu)形狀、加熱方式和粉料的物性參數(shù)有關(guān)之外,還很大程度上受模內(nèi)混合物的流動狀態(tài)的影響,而模內(nèi)混合物的流動狀態(tài)又與模具的轉(zhuǎn)動方式和轉(zhuǎn)速、模內(nèi)粉料的體積百分比直接相關(guān)。Olinek等[12]通過實驗研究認(rèn)為當(dāng)單軸轉(zhuǎn)動模具的轉(zhuǎn)速在4~20 r/min的范圍內(nèi)模內(nèi)粉料的運動狀態(tài)為塌落或者滾落模式,實際商用的滾塑模具的轉(zhuǎn)速都在這個范圍內(nèi)。田昊等[13]通過實驗研究了二維回轉(zhuǎn)圓筒內(nèi)顆粒流的速度分布及顆粒溫度,并重點研究了滾落的流動模式。雖然回轉(zhuǎn)圓筒與實際的滾塑模具有一定的差異,但他們的研究結(jié)論有助于理解粉料的運動情況影響模內(nèi)傳熱系數(shù)的機理。文獻(xiàn)[14]通過對傳熱模型進(jìn)行計算來研究一個繞單軸轉(zhuǎn)動的圓筒形滾塑模具的內(nèi)表面?zhèn)鳠嵯禂?shù)隨模具轉(zhuǎn)速的變化規(guī)律,但沒有研究該傳熱系數(shù)與模內(nèi)粉料的體積百分比的關(guān)系。

    本文提出了兩種研究方法來獲得在加熱階段滾塑模具的內(nèi)表面?zhèn)鳠嵯禂?shù)。第一種方法是實驗與理論相結(jié)合的方法,即首先通過實驗測量了在4種情形下電加熱滾塑模具的外表面溫度和模內(nèi)溫度,然后應(yīng)用能量守恒原理建立一個描述該加熱階段的傳熱模型,并將實測溫度值的擬合多項式代入該傳熱模型的微分方程式中即可計算出模具的內(nèi)表面?zhèn)鳠嵯禂?shù)。第二種方法是將繞雙軸轉(zhuǎn)動的實際模具等效地簡化為一個繞單軸轉(zhuǎn)動的圓筒形模具,然后通過FLUENT軟件的多相流模塊中的Mixture模型對其進(jìn)行仿真計算來得出模具的內(nèi)表面?zhèn)鳠嵯禂?shù)。最后比較了這兩種方法所得的結(jié)果,并得出了該傳熱系數(shù)隨模內(nèi)粉料的體積百分比的變化規(guī)律。所得結(jié)論可為強化電加熱滾塑工藝的模內(nèi)傳熱過程提供理論依據(jù)。

    1 模內(nèi)溫度和模具外表面溫度的實驗結(jié)果

    本實驗所采用的粉料為線性低密度聚乙烯,其牌號為M3204RUP,粉料的表觀密度為317 kg/m3、真實密度為932 kg/m3、熔點為125 ℃、結(jié)晶點為110 ℃、熔體流動速率為4 g/10 min,生產(chǎn)廠家為泰國SCG ICO有限公司。

    電加熱絲的最大功率為700 W/m,滾塑機的內(nèi)、外軸轉(zhuǎn)速分別設(shè)定為11 r/min和1.2 r/min。所采用的滾塑模具如圖1所示,用來測量模具外表面溫度的熱電偶探頭分別安裝在1點和2點的位置,模具外表面上的凹槽用來嵌入電加熱絲。所制造的氫氣瓶的縱向剖面形狀和尺寸如圖2所示,δ為氫氣瓶的設(shè)計厚度。本文將在4種情形下測量模具的外表面溫度和模內(nèi)溫度。這4種情形下的氫氣瓶厚度、模內(nèi)裝填的粉料質(zhì)量、粉料在模具內(nèi)部空間所占的體積百分比如表1所示。

    表1 4種實驗情形的描述Tab.1 Description for four kinds of experimental cases

    圖1 電加熱的滾塑模具Fig.1 Rotational mold heated by electricity

    圖2 氫氣瓶的剖面形狀和尺寸Fig.2 Cross section and dimension of the hydrogen tank

    圖3表示了在情形1時實測的模內(nèi)溫度和模具外表面溫度隨時間的變化規(guī)律。在該圖中的A點是模具外表面的平均溫度上升到135 ℃時所對應(yīng)的時刻,這個溫度比本實驗所用粉料的熔點高10 ℃。本文認(rèn)為從A點開始模內(nèi)的粉料開始熔融,即從初始時刻到A點所對應(yīng)的時刻是本文所研究的加熱階段。以往的研究通常認(rèn)為當(dāng)模具的外表面溫度升高到粉料的熔點時粉料開始熔融[6?8]。但是由于模具沿厚度方向的導(dǎo)熱熱阻和模具內(nèi)表面與模內(nèi)粉料間的對流換熱熱阻都不等于0,因此當(dāng)模具的外表面溫度剛升到粉料的熔點時,模內(nèi)粉料的溫度肯定低于其熔點。從圖3可以看出,在A點時刻之前模具外表面上1點和2點間的溫差很小,即模具表面溫度分布的均勻性較好。但從A點時刻開始此二者間的溫差隨時間顯著增大,這個現(xiàn)象說明此時已有粉料變成了熔融的塑料層,且已經(jīng)覆蓋了部分的模具內(nèi)表面。這是因為熔融塑料層的流動性比熔融前粉料和空氣的混合物差很多,所以此時模具與模內(nèi)混合物的對流換熱的強度變得較弱,模具表面的溫度分布也就變得越來越不均勻。因此選擇當(dāng)模具外表面不同位置處的溫差開始顯著變大的時刻即A點作為粉料開始熔融的時刻是比較合理的。在圖3中B點是模具外表面上的電加熱絲斷電的時刻,即從B點開始模具的外表面溫度開始下降,但由于模具內(nèi)部積蓄的熱量還會繼續(xù)向模內(nèi)傳遞,因此模內(nèi)的溫度還會繼續(xù)升高一段時間。即模內(nèi)的最高溫度比模具外表面的最高溫度滯后一段時間。

    圖3 在情形1實測的模內(nèi)溫度和模具外表面溫度隨時間的變化Fig.3 Variation of tested temperature inside the mold and at outer surface of the mold with time in case 1

    在圖3中還有一個值得注意的現(xiàn)象是在整個加熱過程的前期,模內(nèi)的溫度隨時間呈現(xiàn)較大幅度的波動變化,這是因為在本實驗中伸入模具內(nèi)部的測溫探頭周期性地被塌落或滾落的粉料沖刷。如果把模內(nèi)的粉料和空氣當(dāng)成兩相混合物的話,那么隨著模具的運動,在不同的時刻測溫探頭會周期性地遇到不同體積百分比的兩相混合物,當(dāng)然也包括純空氣和純粉料,從而測量到不同的溫度值??梢哉J(rèn)為在圖3中模內(nèi)溫度的波峰值是當(dāng)測溫探頭遇到純空氣即模內(nèi)空氣的體積百分比為100 %、粉料的體積百分比為0的混合物所測得的溫度值;而模內(nèi)溫度的波谷值是當(dāng)測溫探頭遇到純粉料即空氣百分比為0、粉料百分比為100 %的混合物所測得的溫度值,在這兩者之間對應(yīng)的則是模內(nèi)混合物在其他濃度時的溫度值。在粉料開始熔融后,隨著時間的推移,模內(nèi)殘留的粉料越來越少,另外這些殘留的粉料與空氣也混合得越來越均勻,因此在整個加熱過程的后期,模內(nèi)溫度的波動也就變得越來越小。

    2 實驗數(shù)據(jù)的處理

    為了便于在之后的傳熱模型中應(yīng)用上述實測溫度值來計算模具內(nèi)表面的傳熱系數(shù),需要對這些實驗數(shù)據(jù)進(jìn)行處理。首先將圖3中在模具外表面1點和2點處的實測溫度值取其平均值,求得模具外表面的平均溫度,并將該圖中A點以后的溫度曲線截掉,從而得到圖4。在情形1,A點時刻即模具外表面的平均溫度上升到135℃所對應(yīng)的時刻為171s。

    圖4 在情形1的加熱階段實測的模內(nèi)溫度和模具外表面的平均溫度Fig.4 Temperature inside the mold and average temperature at outer surface of the mold tested for heating phase in case 1

    模內(nèi)粉料和空氣的混合平均溫度應(yīng)該介于圖4中模內(nèi)溫度實測值的波峰與波谷之間,即低于純空氣的溫度,又高于純粉料的溫度。如式(1)所示,對模內(nèi)溫度的實測值采用二階擬合多項式來表示模內(nèi)混合物的平均溫度隨時間的變化規(guī)律。這樣既能保證擬合的精度,同時又能保證擬合所得的溫度值隨時間總是單調(diào)上升。

    式中 a、b、c——擬合多項式的系數(shù)

    t——時間,s

    Tap——模內(nèi)混合物的平均溫度,℃

    如式(2)所示,對實測的模具外表面平均溫度隨時間的變化規(guī)律則采用三次多項式來擬合表示。

    式中 d、e、f、g——擬合多項式的系數(shù)

    Tm——模具外表面的平均溫度,℃

    圖5表示了在情形1~4按式(1)和(2)擬合所得的模內(nèi)混合物及模具外表面的平均溫度的擬合曲線,模具外表面平均溫度上升到135 ℃所對應(yīng)的時刻分別為186、162、147 s。擬合多項式中的所有系數(shù)由Matlab軟件計算給出,如表2所示。

    圖5 在情形1~4的加熱階段模內(nèi)混合物及模具外表面的平均溫度的擬合曲線Fig.5 Fitting curves of average temperature of mixture inside and at outer surface of the mold for heating phase in case 1~4

    表2 在4種情形下實測溫度的擬合多項式中的系數(shù)Tab.2 Coefficients of fitting polynomials of tested temperature in four cases

    3 傳熱模型的建立

    3.1 第一種方法的傳熱模型

    在第一種方法中,需要建立一個傳熱模型才能將實測的模具和模內(nèi)溫度轉(zhuǎn)換為模具內(nèi)表面的傳熱系數(shù)。在建立該模型時,先作出以下兩點假設(shè):(1)、忽略模具沿其壁厚方向的溫度梯度。由于鋁制模具的熱導(dǎo)率比較大,壁厚很薄,因此該假設(shè)與實際情況相差不大。(2)、將模內(nèi)粉料與空氣的混合物當(dāng)成一個整體來處理,在計算中采用模內(nèi)粉料與空氣的混合平均溫度,即采用圖5中所示的擬合之后的模內(nèi)平均溫度。

    根據(jù)以上假設(shè),模具溫度和模內(nèi)溫度都只是時間的函數(shù),而與空間位置坐標(biāo)無關(guān),因此所列出的控制方程為簡單的常微分方程。根據(jù)能量守恒原理可知模具內(nèi)表面與模內(nèi)混合物的對流換熱量等于模內(nèi)混合物內(nèi)能的增量,由此得到本傳熱模型的控制方程(3)。

    初始條件為t=0,Tm=Tm0,Tap=Tap0。

    式中 Ai——模具內(nèi)表面的面積,m2

    cap——模內(nèi)混合物的定壓比熱容,J/(kg·K),根據(jù)模內(nèi)空氣和粉料的質(zhì)量百分比按如式(4)計算所得。

    ca——模內(nèi)空氣的定壓比熱容,J/(kg·K)

    cp——模內(nèi)粉料的定壓比熱容,J/(kg·K)

    hi——模具內(nèi)表面瞬時的傳熱系數(shù),W/(m2·K)

    Ma——模內(nèi)空氣的質(zhì)量,kg

    Mp——模內(nèi)粉料的質(zhì)量,kg

    t——時間,s

    Tap——模內(nèi)混合物的平均溫度,℃

    Tap0——模內(nèi)混合物的初始溫度,℃

    Tm——模具的平均溫度,℃

    Tm0——模具初始溫度,℃

    參數(shù)設(shè)置:Ai=0.362 m2,ca=1 005 J/(kg·K)。根據(jù)文獻(xiàn)[15]的DSC實驗測得粉料在室溫下其比熱容為1 910.7 J/(kg·K),當(dāng)粉料溫度高于134 ℃時其比熱容為2 491.7 J/(kg·K),本文采用二者的平均值為cp=2 201.2 J/(kg·K)。在情形1~4中模內(nèi)粉料的質(zhì)量Mp如表1所示,模內(nèi)空氣質(zhì)量分別為Ma=0.015、0.014、0.013、0.012 kg。

    將上述參數(shù)及擬合多項式(1)和(2)代入式(3)可以分別計算出在4種情形下模具內(nèi)表面瞬時的傳熱系數(shù)hi。在實際應(yīng)用中平均的傳熱系數(shù)更有意義,因此將瞬時的傳熱系數(shù)代入式(5)可計算出在加熱階段模具內(nèi)表面平均的傳熱系數(shù)him,式中的Δt為所研究的加熱階段的時長。以后如不特別提出,本文中模具內(nèi)表面的傳熱系數(shù)指的均是在加熱階段平均的傳熱系數(shù)。

    3.2 第二種方法的傳熱模型

    采用FLUENT軟件的多相流模塊對模內(nèi)粉料與空氣的運動和傳熱進(jìn)行仿真計算,以獲得模具內(nèi)表面的傳熱系數(shù)。就是將模內(nèi)的空氣當(dāng)作兩相流中的主相,由于模內(nèi)粉料的顆粒尺寸非常小,可以近似看成是連續(xù)流體,因此將粉料當(dāng)作是該兩相流中的第二相。粉料流體的黏度系數(shù)沒有任何實驗數(shù)據(jù)可供查詢,也無法由精確的計算公式得出,本文采用式(6)[16]近似計算該黏度系數(shù)。

    式中 μa——模內(nèi)空氣的動力黏度,Pa?s

    μp——粉料流體的動力黏度,Pa?s

    ρa——模內(nèi)空氣的密度,kg/m3

    ρp——粉料的表觀密度,kg/m3

    FLUENT軟件的多相流模塊提供了VOF、Mix?ture、Eulerian 3種模型。VOF模型適用于兩相間不混溶的流體,這與模內(nèi)空氣和粉料的混合運動不相符,故不采用。Eulerian模型在兩相間的阻力定律已知的前提下比Mixture模型的計算結(jié)果更準(zhǔn)確,但目前滾塑模內(nèi)粉料與空氣之間相互作用力的規(guī)律是未知的,而且Eulerian模型求解的方程數(shù)更多,計算量更大,因此本文決定采用較簡單的Mixture模型來仿真求解模具內(nèi)表面的傳熱系數(shù)。

    為進(jìn)一步減小計算量,本文將圖1所示的實際模具簡化為橫截面形狀如圖6所示的一個圓筒,這樣三維問題就簡化為更簡單的二維問題。為保證流動與傳熱的等效性,在簡化時要求模具和模內(nèi)粉料的質(zhì)量都保持不變,模具的內(nèi)徑以及模內(nèi)粉料的體積百分比也都保持不變。 實際的模具繞雙軸轉(zhuǎn)動,在對簡化后的二維問題進(jìn)行仿真時,只考慮圓筒繞軸線以11 r/min轉(zhuǎn)動的情形,即忽略了實際模具以1.2 r/min繞外軸的轉(zhuǎn)動。由于繞外軸的轉(zhuǎn)速很低,因此忽略該轉(zhuǎn)動對仿真精度的影響應(yīng)該比較小。

    圖6 第二種方法的傳熱模型的示意圖Fig.6 Schematic of heat transfer model of the second method

    仿真的區(qū)域包括模具本身和模內(nèi)的粉料和空氣。將圖1所示的模具外表面的電加熱絲等效轉(zhuǎn)換為圖6所示的圓筒壁的內(nèi)熱源,該內(nèi)熱源的強度qv等于實際所用電加熱絲的總發(fā)熱功率除以該圓筒的體積,經(jīng)過計算可得qv=2 047.6 kW/m3。在電加熱絲的總發(fā)熱量中,有一部分熱量Qe被用于加熱模具本身和模內(nèi)的粉料和空氣,另一部分熱量Qc則傳給了外界環(huán)境。Qc包括模具外表面以對流換熱方式傳給環(huán)境空氣的熱量、以輻射換熱方式傳給環(huán)境壁面的熱量以及以導(dǎo)熱方式傳給與其接觸的模架等附件的熱量。如圖6所示,這部分熱量是通過模具的外表面?zhèn)鳠嵯禂?shù)ho來考慮的,也就是說ho是綜合考慮了模具外表面的對流換熱、輻射換熱和導(dǎo)熱的總傳熱系數(shù)。在ho的大小已知后,模具的外表面就可以被定義為第三類熱邊界條件。本文按如下方法計算ho,首先按式(7)計算Qe的大小。因為模內(nèi)空氣的質(zhì)量很小,因此計算Qe時忽略了空氣的吸熱量。

    式中 cm——模具的定壓比熱容,J/(kg·K)

    L——粉料的熔融潛熱,J/kg

    Mm——模具的質(zhì)量,kg

    Tm——在電加熱絲斷電時模具的溫度,℃

    Tp——在電加熱絲斷電時粉料的溫度,℃

    Tm0——模具的初始溫度,℃

    Tp0——粉料的初始溫度,℃

    參數(shù)設(shè)置:cm= 904 J/(kg·K),L = 180 000 J/kg[6],Mm=9.46 kg。在情形 1,通過實驗測得 Tm=230.5 ℃、Tm0=64.2 ℃、Tp=179.7 ℃、Tp0=41.1 ℃。另外通過實驗還可測得從加熱過程開始到電加熱絲斷電所花的時間為677 s,在這段時間模具外表面的平均溫度為180.5 ℃,以及扣除電機空轉(zhuǎn)后的耗電量即實際的加熱電能為4 377.6 kJ。該加熱電能減去由式(7)算出的Qe后得到的就是模具外表面?zhèn)鹘o外界環(huán)境的總熱量,然后通過牛頓冷卻公式[17]即可算出ho=33 W/(m2?K)。采用同樣的方法可算出在情形2,ho=30.1 W/(m2?K),可見不同情形的ho相差不大,而且模具內(nèi)表面的傳熱系數(shù)受外表面?zhèn)鳠嵯禂?shù)的影響很小,因此在所有情形的仿真計算中統(tǒng)一采用ho=33 W/(m2?K)。

    這是一個非穩(wěn)態(tài)的流動與傳熱問題,將模具外表面的對流換熱、模具的導(dǎo)熱與模內(nèi)空間的多相流換熱耦合在一起,將時間步長設(shè)為0.02 s。在每一個時間步長的仿真結(jié)束時,記錄模具和模內(nèi)混合物的平均溫度以及模具內(nèi)表面上的熱流密度,然后根據(jù)牛頓冷卻公式可計算出模具內(nèi)表面在該瞬時的傳熱系數(shù)。當(dāng)計算所得模具的平均溫度上升到135 ℃時,即認(rèn)為模內(nèi)的粉料即將開始熔融,仿真計算結(jié)束。然后將所有的瞬時傳熱系數(shù)代入式(5)即可計算出模具內(nèi)表面平均的傳熱系數(shù)。

    4 計算結(jié)果及分析

    首先按照這兩種方法分別計算出如表1所示的4種情形下的模具內(nèi)表面?zhèn)鳠嵯禂?shù),并進(jìn)行比較,其計算結(jié)果和相對誤差如表3所示。因為第一種方法是以實驗測得的溫度為基礎(chǔ)計算出來的,肯定更接近真實值,因此以第一種方法的結(jié)果作為參照值來計算第二種方法在各情形的相對誤差??梢钥闯鲈谇樾?、2、4下第二種方法對模具內(nèi)表面?zhèn)鳠嵯禂?shù)的計算結(jié)果具有較高的精度,但在情形3下誤差較大。

    表3 2種方法在4種情形下所得的模具內(nèi)表面?zhèn)鳠嵯禂?shù)及其相對誤差Tab.3 Heat transfer coefficients at inner surface of the mold and relative errors in four cases by two methods

    根據(jù)第一種方法在4種情形下所得的結(jié)果以及第二種方法在粉料的體積百分比為0~90 %范圍內(nèi)的仿真結(jié)果,可以得到如圖7所示的模具內(nèi)表面?zhèn)鳠嵯禂?shù)隨模內(nèi)粉料體積百分比的變化規(guī)律??梢钥闯鲈诜哿象w積百分比不超過53 %的時候,兩種方法的結(jié)果吻合得很好。這是因為模內(nèi)的空氣質(zhì)量很小,模具內(nèi)表面與模內(nèi)混合物間的換熱強度主要取決于其與模內(nèi)粉料間的換熱強度。當(dāng)模具轉(zhuǎn)動時,模具內(nèi)表面與粉料間的換熱強度受2個因素的影響,一個是模內(nèi)粉料與模具內(nèi)表面的接觸面積的大小,濃度越高的粉料與模具內(nèi)表面的接觸面積越大,則模內(nèi)傳熱系數(shù)就越大。另一個因素是模內(nèi)粉料在自由面或者說是粉料與空氣的分界面上的運動情況。如圖6所示,當(dāng)粉料處于滾落的運動狀態(tài)時,模具的內(nèi)表面通過摩擦力將溫度較低的粉料從粉料池里帶出到粉料的自由面上,然后沿著切線方向運動,一邊與模內(nèi)的空氣進(jìn)行對流換熱,一邊將熱量傳給粉料池里溫度較低的粉料,因此粉料在自由面上的切向運動可以提高模具的內(nèi)表面?zhèn)鳠嵯禂?shù)。一方面,由于Mixture模型不考慮兩相間的相互作用力,因此采用該模型無法準(zhǔn)確地仿真出粉料在自由面上的運動情況,從而導(dǎo)致第二種方法所得的模內(nèi)傳熱系數(shù)低于第一種方法。但當(dāng)粉料的體積百分比較小時,粉料在自由面上的運動比較弱,因此其對提高模內(nèi)傳熱系數(shù)的作用不是很大,即此時第二種方法的結(jié)果比第一種方法的結(jié)果低得不是很多。另一方面,仿真模型是一個二維的圓筒,沒有端部效應(yīng),而如圖1所示實際的滾塑模具是一個三維實體,其兩端的換熱強度肯定要低于中間的圓筒直段,因此仿真模型的內(nèi)表面?zhèn)鳠嵯禂?shù)又要高于實際滾塑模具的內(nèi)表面?zhèn)鳠嵯禂?shù),即這導(dǎo)致第二種方法的結(jié)果又要高于第一種方法。綜合來考慮,這兩個反向的影響大部分相互抵消了,因此在粉料的體積百分比不高的情形,兩種方法的結(jié)果吻合得很好,即此時采用第二種方法仿真所得的模具內(nèi)表面?zhèn)鳠嵯禂?shù)具有足夠的精度。

    圖7 模具內(nèi)表面的傳熱系數(shù)隨模內(nèi)粉料的體積百分比的變化Fig.7 Variation of heat transfer coefficient at inner surface of the mold with volume percentage of powder

    當(dāng)粉料的體積百分比超過53 %,如果繼續(xù)增大該體積百分比,粉料在自由面上的運動就會極大地強化粉料與空氣間的對流換熱,從而使得模內(nèi)傳熱系數(shù)快速升高。而第二種方法無法仿真出粉料在自由面上運動的強弱變化,所以此時第二種方法計算出的模內(nèi)傳熱系數(shù)比第一種方法低得較多。但是當(dāng)粉料的體積百分比超過70 %以后,由于粉料自由面的面積明顯變小,粉料在自由面的運動變?nèi)趿撕芏?。所以盡管此時隨著粉料體積百分比的提高,粉料與模具內(nèi)表面的接觸面積在增大,在一定程度上有強化模內(nèi)傳熱的趨勢,但還是抵消不了變?nèi)醯姆哿献杂擅孢\動對模內(nèi)傳熱的削弱作用,因此綜合來看模內(nèi)傳熱系數(shù)開始變小了,從而使得兩種方法的結(jié)果趨于接近,即第二種方法的仿真誤差也變小了。

    由圖7還可以看出,當(dāng)粉料的體積百分比達(dá)到62 %時,第二種方法所得的模內(nèi)傳熱系數(shù)隨體積百分比的增加而變化很小。圖8是采用第二種方法對該種情形進(jìn)行仿真所得的模內(nèi)粉料分別在1 s和100 s時的濃度分布云圖??梢钥闯鲇捎贛ixture模型忽略了兩相間的相互作用力,在較高體積百分比的情形下仿真所得的結(jié)果總是粉料較為均勻地分布并充滿整個模內(nèi)空間,也就是說此時繼續(xù)增大體積百分比,均勻濃度的粉料與模具內(nèi)表面的接觸面積總保持為最大值。而由于Mixture模型無法仿真粉料在自由面的運動情況,因此第二種方法所得的模內(nèi)傳熱系數(shù)只取決于粉料與模具內(nèi)表面接觸面積的大小,所以此時繼續(xù)增大粉料的體積百分比,第二種方法的仿真結(jié)果幾乎已經(jīng)達(dá)到了最大值而不再變化。

    圖8 采用第二種方法對體積百分比為62 %的情形仿真所得的模內(nèi)粉料濃度在不同時刻的分布情況Fig.8 Distribution of volume fraction of powder inside the mold at different moments by the second method when the volume percentage of the powder is 62 %

    綜合兩種方法的結(jié)果來看,模具內(nèi)表面的傳熱系數(shù)隨著粉料的體積百分比先是快速增大,然后從10 %開始增大的速率變慢,并在65.3 %時達(dá)到了峰值61.2 W/(m2·K),然后隨體積百分比開始減小。另外當(dāng)粉料的體積百分比不在58 %~74 %的范圍內(nèi),由第二種方法仿真所得的模具內(nèi)表面?zhèn)鳠嵯禂?shù)的相對誤差不超過10 %。考慮到Mixture模型的計算量較小,仿真計算的成本也低于實驗成本,因此在58 %~74 %范圍之外的體積百分比的情形,通過第二種方法來仿真計算模具的內(nèi)表面?zhèn)鳠嵯禂?shù)在滾塑工藝的研發(fā)方面還是具有實際的應(yīng)用價值。

    5 結(jié)論

    (1)為一個電加熱的滾塑模具在粉料開始熔融前的加熱階段提出了兩種研究模具內(nèi)表面?zhèn)鳠嵯禂?shù)的方法,第一種方法是通過傳熱模型將實測的模具和模內(nèi)溫度轉(zhuǎn)換為該傳熱系數(shù),第二種方法是應(yīng)用多相流的Mix?ture模型通過FLUENT軟件仿真計算出該傳熱系數(shù);

    (2)隨著粉料體積百分比的增加,模具內(nèi)表面?zhèn)鳠嵯禂?shù)先是快速增大,然后增大的速率變慢,當(dāng)體積百分比為65.3 %時達(dá)到最大值61.2 W/(m2·K),隨后開始減??;

    (3)兩種方法所得的模具內(nèi)表面?zhèn)鳠嵯禂?shù)在情形1、2、4都吻合得很好,當(dāng)粉料的體積百分比不在58 %~74 %的范圍內(nèi),由第二種方法仿真所得的模具內(nèi)表面?zhèn)鳠嵯禂?shù)的相對誤差不超過10 %。

    猜你喜歡
    粉料傳熱系數(shù)百分比
    近凈成型粉體制備技術(shù)的研究
    廣東建材(2023年10期)2023-10-17 08:25:42
    探析寒冷地區(qū)75%建筑節(jié)能框架下圍護(hù)結(jié)構(gòu)熱工性能的重組
    預(yù)拌混凝土粉料入庫智能防錯系統(tǒng)研究
    四川建筑(2019年6期)2019-07-20 09:39:32
    干法制粉與噴霧干燥制粉的粉料性能對比研究
    佛山陶瓷(2017年7期)2017-09-06 21:27:34
    普通照明用自鎮(zhèn)流LED燈閃爍百分比測量不確定度分析
    電子制作(2017年20期)2017-04-26 06:57:46
    新型鋁合金節(jié)能窗傳熱系數(shù)和簡化計算
    肝癌患者外周血Treg、Th17百分比及IL-17水平觀察
    聚乳酸吹膜過程中傳熱系數(shù)的研究
    中國塑料(2015年2期)2015-10-14 05:34:24
    一種粉料氣力輸送槽
    商貿(mào)投資最適宜國家和地區(qū)
    海外星云 (2012年24期)2012-04-29 00:30:09
    天天影视国产精品| 欧美成人午夜精品| 久久久久精品人妻al黑| av中文乱码字幕在线| 国产精品亚洲一级av第二区| 成人永久免费在线观看视频| 亚洲国产中文字幕在线视频| 国产国语露脸激情在线看| 老鸭窝网址在线观看| а√天堂www在线а√下载 | 国产精品久久视频播放| 亚洲精品中文字幕一二三四区| 亚洲av成人不卡在线观看播放网| 黄片播放在线免费| 中文欧美无线码| 国产精品久久久人人做人人爽| 黑人巨大精品欧美一区二区蜜桃| 一级a爱视频在线免费观看| 国产aⅴ精品一区二区三区波| 久久精品91无色码中文字幕| 少妇裸体淫交视频免费看高清 | 嫁个100分男人电影在线观看| 久久天堂一区二区三区四区| 国产成+人综合+亚洲专区| 欧美中文综合在线视频| 99香蕉大伊视频| 男女下面插进去视频免费观看| 日韩欧美三级三区| 亚洲成av片中文字幕在线观看| 日韩 欧美 亚洲 中文字幕| 美国免费a级毛片| 丰满迷人的少妇在线观看| 久久草成人影院| bbb黄色大片| 超色免费av| 啪啪无遮挡十八禁网站| 美女高潮到喷水免费观看| 午夜福利,免费看| 久久亚洲精品不卡| 亚洲少妇的诱惑av| 妹子高潮喷水视频| 亚洲成人手机| 精品久久蜜臀av无| 精品午夜福利视频在线观看一区| 午夜老司机福利片| 国内毛片毛片毛片毛片毛片| 99国产精品免费福利视频| 久久亚洲精品不卡| 一a级毛片在线观看| 国产精品久久久av美女十八| 无遮挡黄片免费观看| 在线视频色国产色| av中文乱码字幕在线| 成人av一区二区三区在线看| 欧美精品亚洲一区二区| 成人18禁高潮啪啪吃奶动态图| 丰满迷人的少妇在线观看| 亚洲在线自拍视频| 国产aⅴ精品一区二区三区波| 人妻久久中文字幕网| 三级毛片av免费| 日本wwww免费看| 午夜精品久久久久久毛片777| 最近最新免费中文字幕在线| 亚洲男人天堂网一区| 人妻 亚洲 视频| 亚洲片人在线观看| 国产成人欧美| 下体分泌物呈黄色| 一级毛片高清免费大全| 99久久99久久久精品蜜桃| 精品国产乱子伦一区二区三区| 正在播放国产对白刺激| 久久ye,这里只有精品| av片东京热男人的天堂| 精品一区二区三区四区五区乱码| 国产av又大| 免费在线观看黄色视频的| 色精品久久人妻99蜜桃| 亚洲av第一区精品v没综合| 久久国产精品大桥未久av| 久久人妻福利社区极品人妻图片| 国产欧美日韩一区二区三| 丰满的人妻完整版| 日韩 欧美 亚洲 中文字幕| 国产又爽黄色视频| 国产精品香港三级国产av潘金莲| 亚洲国产毛片av蜜桃av| av天堂在线播放| 久久国产精品男人的天堂亚洲| 天天影视国产精品| 免费女性裸体啪啪无遮挡网站| 国产精品秋霞免费鲁丝片| 色婷婷av一区二区三区视频| 在线av久久热| 熟女少妇亚洲综合色aaa.| www.熟女人妻精品国产| 欧美精品av麻豆av| 高清av免费在线| 色综合欧美亚洲国产小说| 日韩欧美在线二视频 | 亚洲午夜精品一区,二区,三区| 不卡av一区二区三区| 高清欧美精品videossex| 在线观看舔阴道视频| 国产亚洲精品一区二区www | 1024香蕉在线观看| 欧美午夜高清在线| 亚洲av美国av| 天天添夜夜摸| 国产成人精品久久二区二区91| 中亚洲国语对白在线视频| 成人手机av| 18禁黄网站禁片午夜丰满| 日本vs欧美在线观看视频| 十分钟在线观看高清视频www| 99久久综合精品五月天人人| 久久人人爽av亚洲精品天堂| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产中文字幕在线视频| 国产高清国产精品国产三级| 国产aⅴ精品一区二区三区波| 久久影院123| 亚洲av片天天在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 人人妻人人添人人爽欧美一区卜| 看免费av毛片| 一级a爱视频在线免费观看| 一边摸一边做爽爽视频免费| 大型黄色视频在线免费观看| 纯流量卡能插随身wifi吗| www.自偷自拍.com| 丁香欧美五月| 黄色女人牲交| 高清视频免费观看一区二区| 国产欧美日韩精品亚洲av| 色94色欧美一区二区| 精品国产乱码久久久久久男人| 久久天躁狠狠躁夜夜2o2o| 又大又爽又粗| 成人手机av| 久久久精品国产亚洲av高清涩受| aaaaa片日本免费| 精品欧美一区二区三区在线| av一本久久久久| 窝窝影院91人妻| 在线国产一区二区在线| 欧美人与性动交α欧美软件| 人妻丰满熟妇av一区二区三区 | 欧美日韩精品网址| 国产不卡一卡二| tube8黄色片| 午夜免费观看网址| 国产精品二区激情视频| 国产免费现黄频在线看| 亚洲一码二码三码区别大吗| 一边摸一边抽搐一进一小说 | 久久香蕉精品热| 悠悠久久av| 久久久国产欧美日韩av| 人人妻,人人澡人人爽秒播| 午夜福利乱码中文字幕| 80岁老熟妇乱子伦牲交| 精品乱码久久久久久99久播| 捣出白浆h1v1| 久久人妻熟女aⅴ| 在线观看日韩欧美| 免费在线观看黄色视频的| 久久久国产一区二区| 国产精品电影一区二区三区 | 亚洲av成人av| 国产有黄有色有爽视频| 一边摸一边抽搐一进一小说 | 亚洲精品中文字幕在线视频| 人人妻,人人澡人人爽秒播| 国产精品久久久人人做人人爽| 岛国毛片在线播放| 在线永久观看黄色视频| 久久中文看片网| 国产av精品麻豆| 久久精品91无色码中文字幕| 国产国语露脸激情在线看| xxx96com| 国产精品一区二区精品视频观看| av天堂在线播放| 亚洲精品美女久久av网站| 少妇裸体淫交视频免费看高清 | 9色porny在线观看| 首页视频小说图片口味搜索| 在线av久久热| 成人国语在线视频| 欧美日韩国产mv在线观看视频| 久久久久久久国产电影| 久久久久国内视频| 超碰成人久久| 色94色欧美一区二区| 国产成人啪精品午夜网站| 动漫黄色视频在线观看| 大型av网站在线播放| 91老司机精品| 亚洲成人免费av在线播放| 亚洲国产欧美一区二区综合| 亚洲欧美一区二区三区久久| 中文字幕人妻丝袜一区二区| 欧美中文综合在线视频| 天天躁夜夜躁狠狠躁躁| 国产成人影院久久av| 亚洲国产欧美网| 精品国产美女av久久久久小说| 无人区码免费观看不卡| 国产99白浆流出| 香蕉丝袜av| 最新美女视频免费是黄的| 日本五十路高清| 麻豆av在线久日| 黑人操中国人逼视频| 国产精品免费一区二区三区在线 | 999精品在线视频| 香蕉丝袜av| 国产aⅴ精品一区二区三区波| 国产深夜福利视频在线观看| 后天国语完整版免费观看| 欧美日韩国产mv在线观看视频| 国产91精品成人一区二区三区| 很黄的视频免费| 在线观看免费视频日本深夜| 亚洲精品中文字幕一二三四区| 免费在线观看日本一区| a级片在线免费高清观看视频| e午夜精品久久久久久久| 动漫黄色视频在线观看| 一级,二级,三级黄色视频| 91精品国产国语对白视频| 巨乳人妻的诱惑在线观看| av福利片在线| 国产99白浆流出| 满18在线观看网站| 王馨瑶露胸无遮挡在线观看| 亚洲熟妇熟女久久| 搡老乐熟女国产| 91成人精品电影| 久久精品亚洲av国产电影网| 777久久人妻少妇嫩草av网站| 91成年电影在线观看| 成人18禁高潮啪啪吃奶动态图| 国产野战对白在线观看| 热99re8久久精品国产| 亚洲五月色婷婷综合| 黄色视频,在线免费观看| 91大片在线观看| 91成人精品电影| 一区在线观看完整版| 欧美精品亚洲一区二区| 久久99一区二区三区| 亚洲精品av麻豆狂野| 日韩欧美一区二区三区在线观看 | 亚洲全国av大片| 天天操日日干夜夜撸| 美女福利国产在线| 亚洲熟妇中文字幕五十中出 | 国产精品99久久99久久久不卡| 免费高清在线观看日韩| 黄色a级毛片大全视频| 男人舔女人的私密视频| 亚洲久久久国产精品| 俄罗斯特黄特色一大片| 国产97色在线日韩免费| 欧美乱码精品一区二区三区| 国产男靠女视频免费网站| 狠狠婷婷综合久久久久久88av| 亚洲人成伊人成综合网2020| 搡老熟女国产l中国老女人| 99久久人妻综合| 51午夜福利影视在线观看| 国产男女内射视频| 在线观看66精品国产| 亚洲久久久国产精品| 亚洲人成伊人成综合网2020| 国产免费男女视频| 欧美日韩国产mv在线观看视频| 丁香欧美五月| 亚洲一区二区三区欧美精品| 高潮久久久久久久久久久不卡| 好男人电影高清在线观看| а√天堂www在线а√下载 | 国产高清激情床上av| 视频在线观看一区二区三区| 国产99久久九九免费精品| av线在线观看网站| 亚洲av日韩在线播放| 欧美精品av麻豆av| 久久婷婷成人综合色麻豆| 亚洲黑人精品在线| 国产精品乱码一区二三区的特点 | 91麻豆精品激情在线观看国产 | 曰老女人黄片| 丝袜美腿诱惑在线| 91九色精品人成在线观看| 亚洲av成人av| 欧美大码av| 午夜福利乱码中文字幕| 美女扒开内裤让男人捅视频| 国产xxxxx性猛交| 精品久久久久久电影网| 午夜精品国产一区二区电影| 在线观看日韩欧美| 999精品在线视频| 国产精品久久视频播放| 精品福利观看| 18禁美女被吸乳视频| 看片在线看免费视频| 99热只有精品国产| 精品第一国产精品| a级片在线免费高清观看视频| 大片电影免费在线观看免费| 久久国产精品男人的天堂亚洲| 免费不卡黄色视频| 成人国产一区最新在线观看| 热99re8久久精品国产| 黄色 视频免费看| 国产精品偷伦视频观看了| 国产精品98久久久久久宅男小说| 69av精品久久久久久| 亚洲va日本ⅴa欧美va伊人久久| 日韩视频一区二区在线观看| 黑人操中国人逼视频| 日韩熟女老妇一区二区性免费视频| 亚洲在线自拍视频| 热99久久久久精品小说推荐| 免费久久久久久久精品成人欧美视频| 亚洲片人在线观看| 国内久久婷婷六月综合欲色啪| 免费少妇av软件| 免费不卡黄色视频| 波多野结衣av一区二区av| 亚洲人成伊人成综合网2020| 日本五十路高清| 精品国产一区二区久久| 亚洲 国产 在线| 一区在线观看完整版| 999久久久精品免费观看国产| 这个男人来自地球电影免费观看| 亚洲三区欧美一区| 久久午夜亚洲精品久久| 婷婷成人精品国产| 免费在线观看黄色视频的| 欧美 日韩 精品 国产| 青草久久国产| 999久久久精品免费观看国产| 一区二区三区精品91| 日韩成人在线观看一区二区三区| 国产av又大| 日韩中文字幕欧美一区二区| 午夜精品在线福利| 一进一出抽搐gif免费好疼 | 人人妻人人添人人爽欧美一区卜| 捣出白浆h1v1| 男男h啪啪无遮挡| 久久香蕉国产精品| 亚洲av成人一区二区三| 香蕉丝袜av| 咕卡用的链子| 亚洲欧美一区二区三区黑人| 欧美不卡视频在线免费观看 | 正在播放国产对白刺激| 黄片小视频在线播放| 国产日韩一区二区三区精品不卡| 亚洲熟妇中文字幕五十中出 | 午夜精品国产一区二区电影| x7x7x7水蜜桃| 村上凉子中文字幕在线| 精品午夜福利视频在线观看一区| 一区二区日韩欧美中文字幕| 欧美中文综合在线视频| 宅男免费午夜| 看免费av毛片| 亚洲av电影在线进入| 国产激情欧美一区二区| 欧美激情 高清一区二区三区| 国产亚洲欧美98| 国产精品影院久久| 午夜老司机福利片| 超碰97精品在线观看| 亚洲黑人精品在线| 久久香蕉精品热| 国产精品偷伦视频观看了| 国产99白浆流出| 一区在线观看完整版| 少妇粗大呻吟视频| aaaaa片日本免费| 村上凉子中文字幕在线| 欧美另类亚洲清纯唯美| 午夜免费观看网址| 国产99久久九九免费精品| 亚洲精品国产精品久久久不卡| 99国产精品99久久久久| 中文字幕另类日韩欧美亚洲嫩草| 一边摸一边抽搐一进一出视频| 久久亚洲真实| 黄网站色视频无遮挡免费观看| 欧美精品一区二区免费开放| 亚洲专区中文字幕在线| aaaaa片日本免费| 色在线成人网| 日韩欧美在线二视频 | 在线国产一区二区在线| 久久人妻熟女aⅴ| 一边摸一边抽搐一进一出视频| 亚洲精华国产精华精| 黄色成人免费大全| 久久性视频一级片| 999久久久国产精品视频| 免费在线观看亚洲国产| 国产成人系列免费观看| 欧洲精品卡2卡3卡4卡5卡区| 一本一本久久a久久精品综合妖精| 国产99久久九九免费精品| 久久久久久久午夜电影 | 丁香欧美五月| 国产片内射在线| 高清视频免费观看一区二区| 热re99久久国产66热| 久久国产精品大桥未久av| 欧美人与性动交α欧美软件| 一级,二级,三级黄色视频| 亚洲三区欧美一区| 久久午夜亚洲精品久久| 久久国产精品影院| 欧美激情极品国产一区二区三区| 别揉我奶头~嗯~啊~动态视频| 欧美精品啪啪一区二区三区| 国产精品一区二区在线观看99| 亚洲国产欧美网| 国产有黄有色有爽视频| 亚洲欧美精品综合一区二区三区| 青草久久国产| 国产精品影院久久| 欧美另类亚洲清纯唯美| 高清视频免费观看一区二区| 亚洲av成人一区二区三| 欧美在线一区亚洲| 999久久久精品免费观看国产| 99精品久久久久人妻精品| 99国产极品粉嫩在线观看| 婷婷丁香在线五月| 宅男免费午夜| 热re99久久国产66热| 一夜夜www| 久久国产精品影院| 亚洲人成电影免费在线| 男女午夜视频在线观看| 黑人猛操日本美女一级片| 国产精品二区激情视频| 国产亚洲精品久久久久久毛片 | 亚洲av成人不卡在线观看播放网| 91麻豆av在线| 亚洲男人天堂网一区| 一区二区三区国产精品乱码| 在线看a的网站| 久久国产精品影院| 国产精品美女特级片免费视频播放器 | 夜夜躁狠狠躁天天躁| 国产免费现黄频在线看| 久久久久国产精品人妻aⅴ院 | av线在线观看网站| 亚洲avbb在线观看| 男女床上黄色一级片免费看| 精品卡一卡二卡四卡免费| 亚洲一码二码三码区别大吗| 亚洲三区欧美一区| 精品久久久久久电影网| 日韩欧美一区视频在线观看| 国产精品电影一区二区三区 | 日韩中文字幕欧美一区二区| 国产欧美日韩综合在线一区二区| 免费少妇av软件| 91大片在线观看| 久久精品人人爽人人爽视色| 欧美不卡视频在线免费观看 | 一级毛片精品| 久久99一区二区三区| 一区二区三区激情视频| 免费不卡黄色视频| 日韩欧美一区二区三区在线观看 | 日韩欧美三级三区| 久久精品国产综合久久久| av网站免费在线观看视频| 妹子高潮喷水视频| 欧美黄色淫秽网站| av视频免费观看在线观看| 日日夜夜操网爽| 一夜夜www| 激情视频va一区二区三区| aaaaa片日本免费| 久热这里只有精品99| 窝窝影院91人妻| 久久精品国产亚洲av高清一级| 成人特级黄色片久久久久久久| 18禁黄网站禁片午夜丰满| 国产不卡一卡二| 精品少妇久久久久久888优播| 久久中文字幕一级| 成人特级黄色片久久久久久久| 日本一区二区免费在线视频| 日韩精品免费视频一区二区三区| av国产精品久久久久影院| 一夜夜www| 在线观看www视频免费| 国产亚洲欧美98| 久久精品国产清高在天天线| 午夜91福利影院| 视频在线观看一区二区三区| 亚洲伊人色综图| 丰满人妻熟妇乱又伦精品不卡| 国产精品.久久久| 曰老女人黄片| 国产精品美女特级片免费视频播放器 | 老司机亚洲免费影院| 欧美成狂野欧美在线观看| 男女午夜视频在线观看| 久久国产精品人妻蜜桃| 99国产综合亚洲精品| 亚洲欧美激情在线| 极品人妻少妇av视频| 日韩欧美一区视频在线观看| 精品高清国产在线一区| 久久精品91无色码中文字幕| 国产在线一区二区三区精| 国内久久婷婷六月综合欲色啪| 日韩三级视频一区二区三区| 欧美精品啪啪一区二区三区| 午夜精品久久久久久毛片777| 国产一区在线观看成人免费| 国产激情欧美一区二区| av网站免费在线观看视频| 久久久久久久国产电影| 夫妻午夜视频| 久久国产精品大桥未久av| 久热爱精品视频在线9| 精品无人区乱码1区二区| 日韩欧美一区二区三区在线观看 | 99在线人妻在线中文字幕 | 后天国语完整版免费观看| 亚洲欧美色中文字幕在线| 人妻丰满熟妇av一区二区三区 | 搡老熟女国产l中国老女人| 18禁裸乳无遮挡免费网站照片 | 精品国产一区二区久久| 精品国产美女av久久久久小说| 9191精品国产免费久久| 国产黄色免费在线视频| 国产野战对白在线观看| 黄色a级毛片大全视频| 狠狠狠狠99中文字幕| 成人av一区二区三区在线看| www.精华液| 露出奶头的视频| 亚洲三区欧美一区| 亚洲第一欧美日韩一区二区三区| 日韩三级视频一区二区三区| 九色亚洲精品在线播放| 日本撒尿小便嘘嘘汇集6| 午夜免费观看网址| 高清欧美精品videossex| 欧美色视频一区免费| 精品视频人人做人人爽| 国产成人精品久久二区二区免费| a在线观看视频网站| 日韩成人在线观看一区二区三区| xxxhd国产人妻xxx| 老鸭窝网址在线观看| 99国产极品粉嫩在线观看| 国产日韩一区二区三区精品不卡| 男女之事视频高清在线观看| 老熟妇仑乱视频hdxx| 夫妻午夜视频| 日韩视频一区二区在线观看| 妹子高潮喷水视频| 国产精品免费大片| 久久精品人人爽人人爽视色| 亚洲成人免费电影在线观看| 欧美一级毛片孕妇| 女人久久www免费人成看片| 亚洲五月色婷婷综合| 久久精品国产99精品国产亚洲性色 | 9191精品国产免费久久| 国产色视频综合| 午夜福利欧美成人| 一区二区三区国产精品乱码| 亚洲五月婷婷丁香| 国产成人系列免费观看| 国产成人免费无遮挡视频| 18禁裸乳无遮挡免费网站照片 | 色婷婷av一区二区三区视频| 国产精品二区激情视频| 国产乱人伦免费视频| 老司机午夜十八禁免费视频| 国产精品二区激情视频| 国产乱人伦免费视频| e午夜精品久久久久久久| 捣出白浆h1v1| 一级黄色大片毛片| 久久草成人影院| 久久精品人人爽人人爽视色| 在线观看www视频免费| 美女扒开内裤让男人捅视频| 国产av精品麻豆| 欧美日韩亚洲高清精品| svipshipincom国产片| 国产亚洲精品久久久久5区| 天天躁夜夜躁狠狠躁躁| xxx96com| 久久久久久久久久久久大奶| 亚洲成人免费av在线播放| 黄色毛片三级朝国网站|