• <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
    亚洲18禁久久av| 麻豆av在线久日| 一二三四在线观看免费中文在| 国产激情欧美一区二区| 夜夜躁狠狠躁天天躁| 夜夜躁狠狠躁天天躁| 亚洲一区中文字幕在线| 亚洲美女视频黄频| 99国产精品一区二区蜜桃av| 校园春色视频在线观看| 国产成+人综合+亚洲专区| 免费在线观看影片大全网站| 亚洲美女黄片视频| 亚洲成人久久性| 搡老熟女国产l中国老女人| 可以免费在线观看a视频的电影网站| 国产av又大| 青草久久国产| 日韩欧美在线二视频| 一本大道久久a久久精品| 久久欧美精品欧美久久欧美| 亚洲片人在线观看| www.999成人在线观看| 丁香欧美五月| 久久久久久九九精品二区国产 | 国内少妇人妻偷人精品xxx网站 | 亚洲专区字幕在线| 亚洲欧美日韩无卡精品| 国产精品久久久久久人妻精品电影| 欧美黄色淫秽网站| 亚洲专区字幕在线| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲,欧美精品.| videosex国产| 90打野战视频偷拍视频| 亚洲欧美日韩无卡精品| 在线a可以看的网站| 91九色精品人成在线观看| 国产精品乱码一区二三区的特点| 一级黄色大片毛片| 欧美日韩黄片免| 亚洲激情在线av| 亚洲av电影在线进入| 欧美zozozo另类| 一个人免费在线观看的高清视频| 日韩精品免费视频一区二区三区| 午夜激情av网站| 熟女电影av网| 精品不卡国产一区二区三区| 午夜两性在线视频| 男女床上黄色一级片免费看| 午夜激情福利司机影院| 欧美 亚洲 国产 日韩一| 午夜影院日韩av| 成在线人永久免费视频| 精品午夜福利视频在线观看一区| 国产熟女xx| 午夜福利高清视频| 午夜免费成人在线视频| 夜夜看夜夜爽夜夜摸| 亚洲在线自拍视频| 国产视频内射| av福利片在线| 狂野欧美白嫩少妇大欣赏| 久久热在线av| 在线观看www视频免费| 国产成人系列免费观看| 亚洲精品久久成人aⅴ小说| 国产真人三级小视频在线观看| 亚洲精品一区av在线观看| 妹子高潮喷水视频| 三级国产精品欧美在线观看 | 亚洲国产日韩欧美精品在线观看 | 亚洲精品在线美女| videosex国产| 欧美成人免费av一区二区三区| 99国产精品一区二区三区| 亚洲精品一区av在线观看| 亚洲精品色激情综合| 一区福利在线观看| 97碰自拍视频| 不卡一级毛片| 老汉色∧v一级毛片| 午夜福利高清视频| 美女午夜性视频免费| 亚洲精品中文字幕一二三四区| 91麻豆av在线| 亚洲欧美精品综合一区二区三区| 丁香欧美五月| 色综合站精品国产| 成人午夜高清在线视频| 给我免费播放毛片高清在线观看| 国产v大片淫在线免费观看| 级片在线观看| 麻豆一二三区av精品| 99精品在免费线老司机午夜| 国产精品永久免费网站| 欧洲精品卡2卡3卡4卡5卡区| 淫妇啪啪啪对白视频| 色老头精品视频在线观看| 亚洲 欧美 日韩 在线 免费| 最好的美女福利视频网| 国产精品一区二区三区四区久久| 桃色一区二区三区在线观看| 国产精品99久久99久久久不卡| 巨乳人妻的诱惑在线观看| av视频在线观看入口| 午夜成年电影在线免费观看| 又紧又爽又黄一区二区| av国产免费在线观看| 亚洲无线在线观看| 欧美黑人欧美精品刺激| 亚洲,欧美精品.| 变态另类丝袜制服| 亚洲精品粉嫩美女一区| 亚洲中文av在线| av超薄肉色丝袜交足视频| 一卡2卡三卡四卡精品乱码亚洲| 老司机在亚洲福利影院| 精品国产美女av久久久久小说| 搡老妇女老女人老熟妇| 欧美日韩福利视频一区二区| 午夜老司机福利片| 18禁黄网站禁片午夜丰满| 高清毛片免费观看视频网站| 成人av在线播放网站| 国产视频内射| 久久午夜综合久久蜜桃| 全区人妻精品视频| a级毛片a级免费在线| 欧美3d第一页| 日本成人三级电影网站| 亚洲av电影在线进入| 可以免费在线观看a视频的电影网站| 超碰成人久久| 国产精品久久久久久亚洲av鲁大| 国产精品久久久人人做人人爽| 国产av一区二区精品久久| 久久婷婷人人爽人人干人人爱| 国产精品免费一区二区三区在线| 中文字幕人妻丝袜一区二区| 国产精品自产拍在线观看55亚洲| 国模一区二区三区四区视频 | 亚洲va日本ⅴa欧美va伊人久久| avwww免费| av视频在线观看入口| www国产在线视频色| 国产片内射在线| 午夜福利成人在线免费观看| 中出人妻视频一区二区| 国产久久久一区二区三区| 黄频高清免费视频| 欧美人与性动交α欧美精品济南到| 欧美另类亚洲清纯唯美| 国产99久久九九免费精品| 国产精品久久久久久久电影 | 国产高清激情床上av| 日韩高清综合在线| 国产成人啪精品午夜网站| 日日摸夜夜添夜夜添小说| 亚洲av片天天在线观看| 亚洲成av人片在线播放无| 亚洲午夜理论影院| 妹子高潮喷水视频| av国产免费在线观看| 欧美性猛交╳xxx乱大交人| 欧美成人免费av一区二区三区| 亚洲熟女毛片儿| 亚洲中文av在线| 亚洲第一欧美日韩一区二区三区| 老汉色∧v一级毛片| 免费看美女性在线毛片视频| 色综合欧美亚洲国产小说| 不卡av一区二区三区| 人人妻,人人澡人人爽秒播| 亚洲美女黄片视频| 日韩精品青青久久久久久| 日韩欧美国产一区二区入口| 久久婷婷人人爽人人干人人爱| 亚洲九九香蕉| 久久久久国内视频| 97超级碰碰碰精品色视频在线观看| 色老头精品视频在线观看| 中文字幕人妻丝袜一区二区| 99热这里只有精品一区 | 免费在线观看影片大全网站| 听说在线观看完整版免费高清| 18禁黄网站禁片午夜丰满| xxx96com| 好看av亚洲va欧美ⅴa在| 欧美高清成人免费视频www| 精品福利观看| 精品电影一区二区在线| 国产主播在线观看一区二区| 午夜福利免费观看在线| 国产区一区二久久| 久热爱精品视频在线9| 国产三级黄色录像| 婷婷丁香在线五月| 一本大道久久a久久精品| 亚洲精品在线观看二区| 国产欧美日韩一区二区三| 国产精品日韩av在线免费观看| 久久久久久久精品吃奶| xxxwww97欧美| 成人欧美大片| 三级男女做爰猛烈吃奶摸视频| 人人妻人人澡欧美一区二区| a级毛片a级免费在线| 国产成人一区二区三区免费视频网站| 国产99久久九九免费精品| 国产野战对白在线观看| 久久香蕉激情| 亚洲成人免费电影在线观看| 亚洲七黄色美女视频| 午夜精品一区二区三区免费看| 一区二区三区高清视频在线| 九色成人免费人妻av| 日本 欧美在线| 亚洲人成伊人成综合网2020| 精品久久久久久成人av| netflix在线观看网站| 精品国产超薄肉色丝袜足j| 五月伊人婷婷丁香| 极品教师在线免费播放| 国产日本99.免费观看| 国产视频一区二区在线看| 久久久精品欧美日韩精品| 日本在线视频免费播放| 亚洲 欧美 日韩 在线 免费| 亚洲精品色激情综合| 99国产精品一区二区蜜桃av| 丝袜人妻中文字幕| 日本撒尿小便嘘嘘汇集6| 日韩大尺度精品在线看网址| 日韩精品中文字幕看吧| 日韩欧美 国产精品| 91成年电影在线观看| 成年女人毛片免费观看观看9| 亚洲精品在线美女| 18美女黄网站色大片免费观看| 在线观看舔阴道视频| 国产成年人精品一区二区| 特大巨黑吊av在线直播| 亚洲乱码一区二区免费版| 久久久久九九精品影院| 脱女人内裤的视频| 久久精品91蜜桃| 国产精品av视频在线免费观看| 长腿黑丝高跟| 国产精品av久久久久免费| 精品欧美国产一区二区三| 蜜桃久久精品国产亚洲av| 老司机在亚洲福利影院| 嫩草影院精品99| 亚洲人成网站高清观看| 丰满的人妻完整版| 欧美成人免费av一区二区三区| 两个人的视频大全免费| 天天添夜夜摸| 首页视频小说图片口味搜索| 91老司机精品| 又黄又粗又硬又大视频| 妹子高潮喷水视频| 国产精品一及| 国产精品久久视频播放| 91麻豆av在线| 久久久久国产精品人妻aⅴ院| 香蕉av资源在线| 成人国语在线视频| av中文乱码字幕在线| 午夜福利视频1000在线观看| 久热爱精品视频在线9| 制服丝袜大香蕉在线| 啦啦啦韩国在线观看视频| 美女免费视频网站| 国产97色在线日韩免费| 午夜免费成人在线视频| 久久久精品国产亚洲av高清涩受| 悠悠久久av| 一边摸一边抽搐一进一小说| 国产精品98久久久久久宅男小说| 99国产精品一区二区蜜桃av| 久99久视频精品免费| 国产免费男女视频| 极品教师在线免费播放| 亚洲成av人片在线播放无| 日韩欧美三级三区| 亚洲人成伊人成综合网2020| 成人高潮视频无遮挡免费网站| 亚洲国产看品久久| 在线观看美女被高潮喷水网站 | 美女大奶头视频| 久久久水蜜桃国产精品网| 亚洲人与动物交配视频| 日韩欧美 国产精品| a在线观看视频网站| 少妇熟女aⅴ在线视频| 欧美性猛交╳xxx乱大交人| 亚洲在线自拍视频| 日韩精品青青久久久久久| 51午夜福利影视在线观看| 国产精品永久免费网站| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久久毛片微露脸| 国产亚洲欧美98| 精品久久久久久久人妻蜜臀av| 熟女电影av网| 五月玫瑰六月丁香| 黄片大片在线免费观看| 国产片内射在线| 亚洲欧美精品综合一区二区三区| 免费在线观看完整版高清| 国产视频内射| 亚洲精品一区av在线观看| 国产成人aa在线观看| 不卡一级毛片| 精品欧美国产一区二区三| 亚洲欧美日韩东京热| 久久精品国产清高在天天线| 亚洲精华国产精华精| 校园春色视频在线观看| 欧美日本视频| 三级国产精品欧美在线观看 | 在线永久观看黄色视频| 亚洲国产欧美一区二区综合| 50天的宝宝边吃奶边哭怎么回事| 日韩 欧美 亚洲 中文字幕| 久久久精品国产亚洲av高清涩受| 久久欧美精品欧美久久欧美| 国内精品久久久久久久电影| 欧美又色又爽又黄视频| 日本一区二区免费在线视频| 日韩精品青青久久久久久| 亚洲人成电影免费在线| 又粗又爽又猛毛片免费看| 一个人免费在线观看电影 | 国内久久婷婷六月综合欲色啪| 久久精品国产综合久久久| 欧美又色又爽又黄视频| 免费在线观看完整版高清| 亚洲一卡2卡3卡4卡5卡精品中文| 成人18禁高潮啪啪吃奶动态图| 成人av一区二区三区在线看| 国产精品精品国产色婷婷| 精品国内亚洲2022精品成人| 99精品欧美一区二区三区四区| 欧美丝袜亚洲另类 | 久久久精品国产亚洲av高清涩受| 不卡av一区二区三区| 国产不卡一卡二| 亚洲国产欧美一区二区综合| 色综合欧美亚洲国产小说| 亚洲av成人一区二区三| 亚洲国产精品999在线| 最近最新中文字幕大全免费视频| √禁漫天堂资源中文www| 国产伦一二天堂av在线观看| 久久99热这里只有精品18| 国产成人一区二区三区免费视频网站| 97碰自拍视频| 99热6这里只有精品| 国内少妇人妻偷人精品xxx网站 | 精品国内亚洲2022精品成人| 一二三四社区在线视频社区8| 黄色a级毛片大全视频| 午夜免费观看网址| 狂野欧美激情性xxxx| 999久久久精品免费观看国产| 亚洲精品美女久久av网站| 国产成人一区二区三区免费视频网站| 又黄又粗又硬又大视频| 美女免费视频网站| 在线观看午夜福利视频| 亚洲av片天天在线观看| 精品久久久久久,| 免费看十八禁软件| 在线国产一区二区在线| 国产精品爽爽va在线观看网站| 日日夜夜操网爽| 亚洲精品美女久久av网站| 99精品在免费线老司机午夜| 国产精品久久久人人做人人爽| 又爽又黄无遮挡网站| 亚洲中文字幕一区二区三区有码在线看 | 国产成人av教育| 中文字幕高清在线视频| 99re在线观看精品视频| 在线观看日韩欧美| 熟妇人妻久久中文字幕3abv| 久久99热这里只有精品18| av福利片在线| 精品电影一区二区在线| xxxwww97欧美| 午夜福利欧美成人| 天堂影院成人在线观看| 亚洲精品美女久久久久99蜜臀| 啦啦啦韩国在线观看视频| 国产一区二区在线av高清观看| 黄色毛片三级朝国网站| 国产精品精品国产色婷婷| 黄色毛片三级朝国网站| 18禁国产床啪视频网站| 免费在线观看亚洲国产| 一二三四在线观看免费中文在| 亚洲,欧美精品.| 欧美日韩福利视频一区二区| 看黄色毛片网站| 久久这里只有精品中国| 亚洲人成电影免费在线| 一级毛片女人18水好多| 久久精品aⅴ一区二区三区四区| 国产97色在线日韩免费| 国产不卡一卡二| svipshipincom国产片| 亚洲成人国产一区在线观看| 国产精品电影一区二区三区| 欧美黄色片欧美黄色片| 亚洲avbb在线观看| 少妇粗大呻吟视频| 亚洲国产欧洲综合997久久,| 给我免费播放毛片高清在线观看| 99精品在免费线老司机午夜| 欧美又色又爽又黄视频| 男人的好看免费观看在线视频 | 97人妻精品一区二区三区麻豆| 变态另类成人亚洲欧美熟女| 91成年电影在线观看| 男女那种视频在线观看| 国产高清激情床上av| 99国产综合亚洲精品| 可以免费在线观看a视频的电影网站| 日韩精品中文字幕看吧| 亚洲18禁久久av| 黄色片一级片一级黄色片| 又黄又爽又免费观看的视频| 亚洲精品色激情综合| 亚洲欧洲精品一区二区精品久久久| 午夜福利高清视频| 欧美3d第一页| 日韩欧美一区二区三区在线观看| svipshipincom国产片| 99re在线观看精品视频| 亚洲一区高清亚洲精品| 每晚都被弄得嗷嗷叫到高潮| 久久亚洲真实| netflix在线观看网站| 亚洲,欧美精品.| 熟女少妇亚洲综合色aaa.| 又粗又爽又猛毛片免费看| 日韩高清综合在线| 亚洲专区字幕在线| 亚洲熟妇熟女久久| 国产高清videossex| 国产一区二区在线观看日韩 | 日韩精品免费视频一区二区三区| 国产单亲对白刺激| 久久久久久九九精品二区国产 | 757午夜福利合集在线观看| 丰满人妻一区二区三区视频av | 日韩大尺度精品在线看网址| av欧美777| 久久久国产成人精品二区| 97碰自拍视频| 亚洲人成电影免费在线| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲天堂国产精品一区在线| 欧美一区二区国产精品久久精品 | 长腿黑丝高跟| 国产成+人综合+亚洲专区| 日本 欧美在线| 日本黄色视频三级网站网址| 亚洲精品久久国产高清桃花| 黄色成人免费大全| 国产精品久久久久久人妻精品电影| 中文字幕av在线有码专区| 国产又色又爽无遮挡免费看| 亚洲人与动物交配视频| 不卡一级毛片| 国模一区二区三区四区视频 | 搞女人的毛片| 国产一区二区在线av高清观看| 高潮久久久久久久久久久不卡| 国产精品久久久久久精品电影| 国产精品一区二区免费欧美| 国产精品爽爽va在线观看网站| 欧美激情久久久久久爽电影| 久久午夜亚洲精品久久| 亚洲成人精品中文字幕电影| 性欧美人与动物交配| 久久国产乱子伦精品免费另类| 无人区码免费观看不卡| 悠悠久久av| 男插女下体视频免费在线播放| 亚洲一区中文字幕在线| 最近最新免费中文字幕在线| 亚洲国产高清在线一区二区三| 国产av又大| 国产一区二区三区视频了| 香蕉国产在线看| 亚洲中文字幕日韩| 久久久久久久精品吃奶| 亚洲国产日韩欧美精品在线观看 | 亚洲欧洲精品一区二区精品久久久| 国产真人三级小视频在线观看| 夜夜夜夜夜久久久久| 亚洲专区中文字幕在线| av福利片在线| 18禁观看日本| 亚洲男人的天堂狠狠| 亚洲18禁久久av| 成人永久免费在线观看视频| 一个人免费在线观看电影 | 757午夜福利合集在线观看| 制服丝袜大香蕉在线| 亚洲av成人一区二区三| 1024香蕉在线观看| 亚洲av电影不卡..在线观看| 日韩欧美精品v在线| 欧美色欧美亚洲另类二区| 日韩免费av在线播放| 日本一区二区免费在线视频| 亚洲成人久久性| 亚洲一卡2卡3卡4卡5卡精品中文| 两个人的视频大全免费| 日韩有码中文字幕| 免费在线观看影片大全网站| 性欧美人与动物交配| 岛国视频午夜一区免费看| 欧美日韩国产亚洲二区| 久久性视频一级片| 亚洲国产精品999在线| 久久久久久久久免费视频了| 久久中文字幕一级| 欧美三级亚洲精品| www国产在线视频色| 五月玫瑰六月丁香| 亚洲av成人一区二区三| 欧美日韩乱码在线| 免费无遮挡裸体视频| 日本 欧美在线| 成人高潮视频无遮挡免费网站| 午夜亚洲福利在线播放| 色综合婷婷激情| 日韩三级视频一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 一进一出好大好爽视频| 国产久久久一区二区三区| 两人在一起打扑克的视频| 国产乱人伦免费视频| 亚洲免费av在线视频| 欧美精品亚洲一区二区| 免费人成视频x8x8入口观看| 一本综合久久免费| а√天堂www在线а√下载| 老司机深夜福利视频在线观看| 长腿黑丝高跟| 国产精品电影一区二区三区| 九色成人免费人妻av| 丁香六月欧美| 嫩草影院精品99| 国产伦一二天堂av在线观看| 久久中文字幕一级| 在线十欧美十亚洲十日本专区| 人成视频在线观看免费观看| 国产亚洲精品久久久久久毛片| 91字幕亚洲| av超薄肉色丝袜交足视频| 色哟哟哟哟哟哟| 国产日本99.免费观看| 亚洲真实伦在线观看| 久久久久国产一级毛片高清牌| 国产主播在线观看一区二区| 日本三级黄在线观看| 亚洲成人国产一区在线观看| 一级毛片精品| 性色av乱码一区二区三区2| 精品福利观看| 精品欧美国产一区二区三| 嫩草影视91久久| 成人三级做爰电影| 18禁国产床啪视频网站| 中文字幕高清在线视频| 久久天躁狠狠躁夜夜2o2o| 欧美日韩黄片免| 很黄的视频免费| 午夜福利欧美成人| 欧美久久黑人一区二区| 露出奶头的视频| 手机成人av网站| 两性夫妻黄色片| 亚洲av成人不卡在线观看播放网| 亚洲欧美精品综合久久99| 欧美午夜高清在线| 一级毛片精品| 51午夜福利影视在线观看| 国产精品亚洲一级av第二区| 一进一出抽搐动态| 国产成人系列免费观看| 亚洲第一电影网av| 亚洲熟妇熟女久久| 不卡av一区二区三区| 欧美成人一区二区免费高清观看 | 婷婷精品国产亚洲av| 欧美成狂野欧美在线观看| 神马国产精品三级电影在线观看 | 久久亚洲精品不卡| 中文字幕久久专区| 免费看美女性在线毛片视频| 又紧又爽又黄一区二区| 日韩有码中文字幕| av国产免费在线观看| 精品一区二区三区av网在线观看| 欧美性猛交黑人性爽| 无限看片的www在线观看|