999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

預處理對無核白葡萄熱風干燥特性的影響

2014-01-19 03:50:20關志強吳寶川
食品科學 2014年5期
關鍵詞:擴散系數模型

康 彥,關志強*,李 敏,吳寶川

(廣東省水產品加工與安全重點實驗室,廣東普通高等學校水產品深加工重點實驗室,廣東海洋大學食品科技學院,廣東 湛江 524088)

預處理對無核白葡萄熱風干燥特性的影響

康 彥,關志強*,李 敏,吳寶川

(廣東省水產品加工與安全重點實驗室,廣東普通高等學校水產品深加工重點實驗室,廣東海洋大學食品科技學院,廣東 湛江 524088)

為研究預處理對無核白葡萄熱風干燥特性的影響,用5 g/100 mL碳酸鉀和0.6 g/100 mL橄欖油混合液浸泡無核白葡萄,進行不同干燥溫度條件下的熱風干燥實驗。然后對4種常見的農產品薄層干燥模型進行非線性擬合,并比較評價決定系數(R2)、卡方(χ2)值和標準誤差(RMSE)。結果表明:Parabolic模型更加適合描述無核白葡萄熱風干燥 水分比與干燥時間之間的關系。在干燥溫度分別為50、60、70 ℃時,預處理組的干燥時間比對照組的干燥時間 分別縮短33.3%、19.2%、15.4%。預處理組的無核白葡萄的有效水分擴散系數分別提高47.4%、32.0%、32.4%;預處理組的平均活化能比對照組降低16.1%。預處理組葡萄干制品的復水率大于對照組。干燥前預處理可以縮短干燥時間、提高物料有效水分擴散系數、降低干燥活化能和保持干制品良好品質。

葡萄;預處理;干燥模型;水分擴散率;活化能

葡萄是一種世界性水果,是加工比例最高、國際貿易量最大的水果之一[1]。葡萄是一種季節性水果,除鮮食外,主要加工制品為葡萄酒、葡萄干、葡萄汁。葡萄干的傳統制作方法是使用晾房設備自然干燥,這種方式耗時長,不易控制且衛生狀況不能保證。熱風干燥簡便,干燥效率高,控制性良好。在葡萄干制前預處理的研究上,國內外同行作了大量的工作,干燥預處理的添加物質主要包括堿性物質、脂肪酸類和脂類,常見的預處理試劑包括:氫氧化鈉、氫氧化鉀、硫酸鉀、碳酸鉀、甲酯或乙酯的堿性乳濁液的混合物[2-7]。碳酸鉀可以增加葡萄表皮對水分的滲透性,橄欖油是一種具有美容等多重保健功效的天然植物油,這兩種添加劑對提高葡萄干燥速率和保持葡萄干的品質有很好的作用[6-7]。而碳酸鉀和橄欖油復配用于干制黑葡萄、蘋果渣等產品,使用較少且對其作用探討不夠深入。本實驗采用熱風干燥方式對無核白葡萄進行干燥,預處理中加入一定比例的碳酸鉀和橄欖油混合物,探討預處理對熱風干燥特性的影響,確定擬合度較好的干燥模型,為葡萄干制品的工業化生產提供理論參考依據。

1 材料與方法

1.1 材料與試劑

無核白葡萄購置于廣東省湛江霞山農貿市場,平均直徑為(1.56±0.2)cm。直接干燥法(GB 5009.3—2010《食品中水分的測定》)測得無核白葡萄平均含水率為(78.5±1.0)%。

預處理試劑碳酸鉀和橄欖油均為分析純和食品級。

1.2 實驗步驟

配制5 g/100 mL和0.6 g/100 mL的碳酸鉀和一定量橄欖油的混合預處理溶液,將配好的預處理溶液放入40 ℃水浴鍋進行預熱。稱取無核白葡萄(50.0±2.0)g,浸入預處理溶液,40 ℃浸泡3 min,取出,吸干表面水分后晾干(約1 h),放入熱風干燥裝置進行干燥。熱風干燥條件為:熱風溫度設定為50、60、70 ℃,風速2.0 m/s[6,8-9]。放入樣品前,先對熱風干燥設備進行調試,約2 h達到穩定的設定條件后再進行實驗。從放入樣品開始計時,每隔1 h用電子天平記錄一次無核白葡萄的質量,計算干基含水率。干燥到無核白葡萄的干基含水率為0.25±0.30時終止實驗。每組溫度的實驗重復兩次取平均值。冷卻后稱取適量干燥樣品測量復水率。剩余樣品放入聚乙烯樣品袋常溫密封保存。

1.3 干燥模型

無核白葡萄的干燥水分比由公式(1)計算。

式中:MR為水分比;Mt為任意時刻無核白葡萄的干基含水率;Me為終止干基含水率;M0為初始干基含水率。

選用4種常用的農產品干燥模型:

1)Lewis模型,用于描述大麥[10]和葡萄[11]的干制過程:

2)Page模型,成功地用于描述一些農產品[12-15]的干燥特性:

3)Parabolic模型,用于描述葡萄[6]的干燥特性:

4)Logarithmic模型,一般適用于描述桃子[3]、杏[4]和土豆[16]的干燥特性:

式中:t為干燥時間/h。

1.4 數據分析

用JMP7.0軟件對干燥實驗數據進行回歸分析,由R2、χ2和RMSE來評價模型擬合度。R2越大、χ2和RMSE越小,說明擬合度越好[17-19]。χ2和RMSE由以式(6)、(7)計算。

式中:N指總項數;z指第z項。

1.5 有效水分擴散系數

農產品的干燥特性可由菲克擴散方程來描述。假設:非穩態擴散在球面坐標中的水分遷移過程中忽略物質收縮性,擴散系數為常數,且溫度恒定[20],由菲克擴散方程推導得:

式中:Deff是有效水分擴散系數/(m2/s);r是物質半徑/m。對于較長時間的干燥過程,方程(8)還可以進一步簡化,簡化后方程為:

有效水分擴散系數可由lnMR對時間的曲線獲得。由方程(9)可以看出,lnMR-t是一條直線,斜率K為:

在不同的干燥條件下,用實驗數據擬合lnMR-t直線方程,根據直線方程的斜率計算Deff。

1.6 活化能

與溫度有關的活化能由阿倫尼烏斯方程[21]計算:

式中:D0為阿倫尼烏斯方程前置因子/(m2/s);Ea為活化能/(kJ/mol);R為氣體常數/(kJ/(mol·K));T為溫度/℃。

在不同的干燥溫度下,將方程兩邊取對數,對lnDeff與1/T進行線性擬合,由擬合直線的斜率計算出活化能Ea。

1.7 復水率

稱取一定質量的干燥樣品(m1)放入燒杯中加500 mL水,40 ℃水浴5 h后取出,用濾紙吸干表面水分,在電子天平上稱質量(m2),計算公式為:

2 結果與分析

2.1 干燥數學模型的擬合

把干燥實驗的干基含水率轉換成水分比,對前述的4個模型(式(2)~(5))進行非線性回歸。4個模型的擬合度對比由R2、χ2和RMSE進行比較,數據分析結果見表1~3。

表1 50 ℃薄層干燥模型擬合結果Table 1 Statistical results from various thin-layer drying models at 50 ℃

表2 60 ℃薄層干燥模型擬合結果Table 2 Statistical results from various t hin-layer drying models at 60 ℃Table 2 Statistical results from various t hin-layer drying models at 60 ℃

表3 70 ℃薄層干燥模型擬合結果Table 3 Statistical results from various t hin-layer drying models at 70 ℃Table 3 Statistical results from various t hin-layer drying models at 70 ℃

由表1~3可知,對于對照組和預處理組,4個模型中,Parabolic模型的R2最大、χ2和RMSE相對較小,R2范圍為0.803 39~0.999 11、χ2范圍為0.000 05~0.008 36、RMSE范圍為0.010 55~0.028 91。這說明無核白葡萄的干燥特性與二次多項式Parabolic模型最符合,分析結果與大蒜[22]和無核小葡萄[6]的模型分析相似,但是與黑葡萄[7]的最適模型Page模型不同,這可能與葡萄種類和干燥條件有關。

2.2 不同預處理對干燥時間的影響

圖1 預處理對無核白葡萄熱風干燥時間的影響Fig.1 Effect of pretreatment on drying time of white grapes

由圖1可知,隨著干燥過程的進行,無核白葡萄含水率不斷下降,且預處理組比對照組下降更快,干燥時間明顯縮短。要達到預期的最終含水率,當干燥溫度為50 ℃,對照組需要42 h,預處理組只需28 h;干燥溫度為60 ℃時,對照組需要26 h,預處理組只需21 h;干燥溫度為70 ℃時,對照組需要13 h,預處理組只需11 h。干燥溫度為50、60、70 ℃時,經過預處理的無核白葡萄干燥時間分別縮短33.3%、19.2%、15.4%,這與文獻[7, 23-26]的干燥預處理規律相似。干燥溫度相對較低時,干燥時間縮短的比例相對較大,說明預處理對縮短干燥時間的貢獻率隨著干燥溫度的降低而增大,即在一定干燥溫度范圍內,溫度越低,預處理對縮短干燥時間的作用越明顯。由此可見,采用預處理,可以顯著提高干燥速率,縮短干燥時間。

2.3 預處理對有效水分擴散系數的影響

有效水分擴散系數越大,則水分擴散能力越強,干燥速率越快。一般食品材料的有效水分擴散系數范圍為10-10~10-8m2/s。將實驗數據轉換為lnMR-t,并進行線性擬合,擬合方程的斜率為K,根據式(10)計算得Deff。不同干燥條件下無核白葡萄的有效水分擴散系數如表4所示,有效水分擴散系數范圍為1.9×10-10~4.9×10-10m2/s,這一結果與其他干燥實驗的數據相似[27-30]。

表4 不同干燥條件下無核白葡萄的有效水分擴散系數Table 4 The effective diffusivity of white grapes at different drying condittiioonnss

從圖2、表4可知,有效水分擴散系數隨著溫度的升高而增加,預處理組有效水分擴散系數大于對照組有效水分擴散系數。50 ℃預處理組和對照組水分擴散系數分別為2.8×10-10m2/s和1.9×10-10m2/s;60 ℃預處理組和對照組水分擴散系數分別為3.3×10-10m2/s和2.5×10-10m2/s;70 ℃預處理組和對照組水分擴散系數分別為4.9×10-10m2/s和3.7×10-10m2/s,預處理組比對照組的有效水分擴散系數分別提高了47.4%、32.0%、32.4%。從結果可以看出,干燥預處理可以顯著提高有效水分擴散系數,促進水分向物料表面擴散,從而提高干燥速率。預處理添加劑的作用就是增加干燥物料內外的滲透壓,從而加快水分擴散。一般有效水分擴散系數隨著溫度的升高而增大,上述結果顯示,干燥溫度為50 ℃時有效水分擴散系數提高的幅度最大,干燥溫度為60 ℃和70 ℃時有效水分擴散系數提高的幅度相當。這可能因為,當干燥溫度升高到60~70 ℃時,干燥表面傳入的能量大于水分去除所需能量時,物料微觀組織結構發生了變化,改變了原有水分脫除的通道,少量的糖分外滲使物料表面黏性變大,對內外水分擴散也產生較大的影響,這些復雜的因素使得有效水分擴散系數提高的幅度與干燥溫度的關系不呈單調變化規律。

圖2 不同干燥溫度的有效水分擴散系數Fig.2 Variation of effective diffusivity with drying temperature

2.4 預處理對平均活化能的影響

將式(11)兩邊取對數,然后將lnDeff與1/(T+273.15)的曲線進行線性擬合,根據擬合直線的斜率計算出無核白葡萄的活化能Ea。式(13)和(14)分別表示預處理組和對照組溫度的阿倫尼烏斯關系。一般食品材料干燥過程的活化能范圍為12.7~110 kJ/mol[27]。

由圖3可知,本實驗中預處理組和對照組的平均活化能分別為23.26 kJ/mol和27.71 kJ/mol,預處理組的平均活化能比對照組的平均活化能降低16.1%。活化能越低,水分越容易除去。由此可見,預處理對無核白葡萄干燥過程具有促進作用。

圖3 有效水分擴散系數與絕對溫度倒數關系Fig.3 Effective diffusivity vs. reciprocal absolute temperature

2.5 預處理對復水率的影響

復水率是評價干燥品品質的一個重要指標,它可以反應干燥過程或預處理對干燥產品的破壞程度,較低的復水率可能是由于樣品組織結構破壞程度較大造成的。

表5 不同干燥條件下無核白葡萄的復水率Table 5 Rehydration rates of white grapes at different drying conditions %

由表5可知,在相同溫度下,預處理組的復水率都比對照組高,且70 ℃預處理組的復水率最高,這與柿子片[8]、甜櫻桃[9]、檳榔葉[31]等干燥產品復水率變化的結論一致。

3 結 論

3.1 比較評價4種常用農產品薄層干燥模型非線性擬合的R2、χ2和RMSE,Parabolic模型的模擬效果最好,更加適合描述無核白葡萄熱風干燥水分比與干燥時間之間的關系。

3.2 無核白葡萄經5 g/100 mL碳酸鉀和0.6 g/100 mL橄欖油混合液進行干燥預處理后,能明顯改善干燥效果。通過比較預處理組和對照組的干燥特性,得出如下結論:干燥前預處理可以縮短干燥時間、提高物料有效水分擴散系數和降低干燥活化能,達到降低干燥能耗的目的。

3.3 預處理得到的干制品,復水率比對照組高,說明本實驗設計的預處理方法不僅可以提高干燥速率,而且對保持干制品的品質有較好的效果。

[1] 喬憲生. 中國葡萄及其加工品貿易現狀與發展方向[J]. 中國食物與營養, 2004(3): 27-31.

[2] SARAVACOS G D, MAROUSIS S N, RAOUZEOS G S. Effect of ethyl oleate on the rate of air-drying of foods[J]. Journal of Food Engineering, 1988, 7(4): 263-270.

[3] DOYMAZ I. Effect of pre-treatments using potassium metabisulphide and alkaline ethyl oleate on the drying kinetics of apricots[J]. Biosystems Engineering, 2004, 89(3): 281-287.

[4] KINGSLY R P, GOYAL R K, MANIKANTAN M R, et al. Effects of pretreatments and drying air temperature on drying behavior of peach slice[J]. International Journal of Food Science and Technology, 2007, 42(1): 65-69.

[5] KASSEM A S, SHOKR A Z, EL-MAHDY A R, et al. Comparison of drying characteristics of ‘Thompson’ seedless grapes using combined microwave oven and hot air drying[J]. Journal of the Saudi Society of Agricultural Sciences, 2011, 10(1): 33-40.

[6] DOYMAZ ?, ALTINER P. Effect of pretreatment solution on drying and color characteristics of seedless grapes[J]. Food Science and Biotechnology, 2012, 21(1): 43-49.

[7] DOYMAZ ?. Drying kinetics of black grapes treated with different solutions [J]. Journal of Food Engineering, 2006, 76(2): 212-217.

[8] DOYMAZ ?. Evalution of some thin-layer drying models of persimmon slices (Diospyros kaki L.)[J]. Energy Conversion and Management, 2012, 56: 199-205.

[9] DOYMAZ ?, ?SMAIL O. Drying characteristics of sweet cherry[J]. Food and Bioproducts Processing, 2011, 89(1): 31-38.

[10] BRUCE D M. Exposed-layer barley drying, three models fi tted to new data up to 150 ℃ [J]. Journal of Agricultural Engineering Research, 1985, 32(4): 337-347.

[11] ROBERTS J S, KIDD D R, PADILLA-ZAKOUR O. Drying kinetics of grape seeds[J]. Journal of Food Engineering, 2008, 89(4): 460-465.

[12] da SILVA M A, PINEDO R A, KIECKBUSCH T G. Ascorbic acid thermal degradation during hot air drying of Camu-Camu (Myrciaria dubia [H.B.K.] Mc Vaugh) slices at different air temperatures[J]. Drying Technology, 2005, 23(9/11): 2277-2287.

[13] HASSAN-BEYGI S R, AGHBASHLO M, KIANMEHR M H, et al. Drying characteristics of walnut (Juglans regia L.) during convection drying [J]. International Agrophysics, 2009, 23(2): 129-135.

[14] SINGH S, RAINA C S, BAWA A S, et al. Effect of pretreatments on drying and rehydration kinetics and color of sweet potato slices[J]. Drying Technology, 2006, 24(10): 1487-1494.

[15] JANGAM S V, JOSHI V S, MUJUMDAR A S, et al. Studies of dehydration of sapota (Achras zapota) [J]. Drying Technology, 2008, 26(3): 369-377.

[16] KHAZAEI J, CHEGINI G R, MARYAM B. A novel alternative method for modeling the effects of air temperature and slice thickness on quality and drying kinetics of tomato slices: superposition technique[J]. Drying Technology, 2008, 26(6): 759-775.

[17] SACILIK K, ELICIN A K, UNAL G. Drying kinetics of Ueryani plum in a convective hot-air dryer[J]. Journal of Food Engineering, 2006, 76(3): 362-368.

[18] HASSAN-BEYGI S R, AGHBASHLO M, KIANMEHR M H, et al. Drying characteristics of walnut (Juglans regia L.) during convection drying[J]. International Agrophysics, 2009, 23(2): 129-135.

[19] SOBUKOLA O. Effect of pre-treatment on the drying characteristics and kinetics of okra (Abelmoschus esculetus(L.) Moench) slices[J]. International Journal of Food Engineering, 2009, 5(2), doi: 10.2202/1556-3758.1191.

[20] CRANK J. The mathematics of diffusion[M]. Oxford : Clarendon Press, 1975.

[21] SIMAL S, FEMENIA A, C?RCEL J A, et al. Mathematical modeling of the drying curves of kiwi fruits: influence of the ripening stage[J]. Journal of the Science of Food and Agriculture, 2005, 85: 425-432.

[22] SHARMA G P, PRASAD S. Effective moisture diffusivity of garlic cloves undergoing microwave-convective drying[J]. Journal of Food Engineering, 2004, 65(4): 609-617.

[23] DOYMAZ I. Effect of pre-treatments using potassium metabisulphide and alkaline ethyl oleate on the drying kinetics of apricots[J]. Biosystems Engineering, 2004, 89(3): 281-287.

[24] TARHAN S, ERGUNES G, TASER O F. Selection of chemical and thermal pretreatment combination to reduce the dehydration time of sour cherry (Prunus cerasus L.)[J]. Journal of Food Process Engineering, 2006, 29(6): 651-663.

[25] EISSEN W, MüHLBAUER W, KUTZBACH H D. Solar drying of grapes[J]. Drying Technology, 1985, 3(1): 63-74.

[26] VEGA-G?LVEZ A, MIRANDA M, BILBAO-S?INZ C, et al. Empirical modeling of drying process for apple (Cv. Granny Smith) slices at different air temperatures[J]. Journal of Food Processing & Preservation, 2008, 32(6): 972-986.

[27] ZOGZAS N P, MAROULIS Z B, MARINOS-KOURIS D. Moisture diffusivity data compilation in foodstuffs[J]. Drying Technology, 1996, 14(10): 2225-2253.

[28] TOGRUL I T. Modelling of heat and moisture transport during drying black grapes[J]. International Journal of Food Science and Technology, 2010, 14(6): 23-32.

[29] XIAO Hongwei, PANG Changle, WANG Lihong, et al. Drying kinetics and quality of ‘Monukka’ seedless grapes dried in an air-impingement jet dryer[J]. Biosystems Engineering, 2010, 105(2): 233-240.

[30] V?ZQUEZ G, CHENLO F, MOREIRA R, et al. Effects of various treatments on the drying kinetics of ‘Muscatel’ grapes[J]. Drying Technology, 2000, 18(9): 2131-2144.

[31] BALASUBRAMANIAN S, SHARMA R, GUPTA R K, et al. Validation of drying models and rehydration characteristics of betel (Piper betel L.) leaves [J]. Journal of Food Science and Technology, 2011, 48(6): 685-691.

Effect of Pretreatment on Hot Air Drying Characteristics of White Seedless Grapes

KANG Yan, GUAN Zhi-qiang*, LI Min, WU Bao-chuan
(Guangdong Provincial Key Laboratory of Aquatic Product Process and Safety, Key Laboratory of Advanced Processing of Aquatic Products of Guangdong Higher Education Institution, College of Food Science and Technology, Guangdong Ocean University, Zhanjiang 524088, China)

The effects of pretreatment by soaking in a mixed solvent system containing 5 g/100 mL potassium carbonate and 0.6 g/100 mL olive oil and hot-air temperature (50, 60 and 70 ℃) on drying characteristics of white grapes were studied using a hot air dryer. Non-linear regression analysis was used to develop four common agricultural thin layer-drying models based on the data obtained and these models were evaluated by coeffi ci ent of determination (R2), chi-square value (χ2) and root mean squared error (RMSE). It turned out that the Parabolic model was more appropriate to describe the relationship between moisture ratio and drying time. Compared with the control group, the drying time for the pretreatment group was shorten by 33.3%, 19.2% and 15.4% at 50, 60 and 70 ℃, respectively, the effective diffusivity of white grapes was increased by 47.4%, 32.0% and 32.4%, respectively. In addition, the activation energy was decreased by 16.1%, and the rehydration of dried grapes from the pretreatment group was higher. This pretreatment is not only good for shortening the drying time, improving effective moisture diffusivity and reducing activation energy, but also help preserve the quality of dried products.

grape; pre-treatment; drying mathematical models; moisture diffusivity; activation energy

TS255.4

A

1002-6630(2014)05-0119-05

10.7506/spkx1002-6630-201405024

2013-02-04

康彥(1987—),女,碩士研究生,主要從事食品加工與貯藏研究。E-mail:kangyan@126.com

*通信作者:關志強(1956—),男,教授,碩士,主要從事食品冷凍與干燥工程技術研究。E-mail:mmcgzq@163.com

猜你喜歡
擴散系數模型
一半模型
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
一類具有變擴散系數的非局部反應-擴散方程解的爆破分析
3D打印中的模型分割與打包
基于Sauer-Freise 方法的Co- Mn 體系fcc 相互擴散系數的研究
上海金屬(2015年5期)2015-11-29 01:13:59
FCC Ni-Cu 及Ni-Mn 合金互擴散系數測定
上海金屬(2015年6期)2015-11-29 01:09:09
FLUKA幾何模型到CAD幾何模型轉換方法初步研究
非時齊擴散模型中擴散系數的局部估計
Ni-Te 系統的擴散激活能和擴散系數研究
上海金屬(2013年4期)2013-12-20 07:57:07
主站蜘蛛池模板: 欧美亚洲一区二区三区在线| 亚洲五月激情网| 欧美日韩国产综合视频在线观看| 欧美日韩中文字幕二区三区| 少妇被粗大的猛烈进出免费视频| 成人永久免费A∨一级在线播放| 欧美福利在线| 狼友视频一区二区三区| 亚洲精品国产自在现线最新| 国产主播喷水| 精品欧美日韩国产日漫一区不卡| 亚洲人妖在线| 99久久性生片| 婷婷色中文| 日本久久免费| 99热这里只有精品国产99| 五月天久久婷婷| 欧美亚洲国产日韩电影在线| 成人字幕网视频在线观看| 精品亚洲欧美中文字幕在线看| 亚洲精品成人福利在线电影| 国产视频你懂得| 精品久久人人爽人人玩人人妻| 久久a毛片| 在线观看精品自拍视频| 久久国产免费观看| yjizz国产在线视频网| 亚洲综合精品香蕉久久网| 日本精品视频一区二区| 亚洲无码视频一区二区三区 | 国产老女人精品免费视频| 国产成人福利在线| 国产日韩av在线播放| 国产精品香蕉| 国产成人麻豆精品| 免费看av在线网站网址| 91九色国产在线| 就去色综合| 99热精品久久| 久久中文电影| 久草网视频在线| 久久精品中文字幕免费| 亚洲色成人www在线观看| 国产剧情一区二区| 日本久久网站| 黄色在线不卡| 国产第八页| 精品自拍视频在线观看| 曰韩人妻一区二区三区| 日韩国产 在线| 国产综合网站| 久操线在视频在线观看| 999国内精品久久免费视频| 国产成人一区| 国语少妇高潮| aa级毛片毛片免费观看久| 高清色本在线www| 国产在线91在线电影| 自慰网址在线观看| 亚洲天堂网2014| 国产成人AV综合久久| 欧美天堂在线| 囯产av无码片毛片一级| 国产精品不卡片视频免费观看| 国产黄色免费看| 妇女自拍偷自拍亚洲精品| 69视频国产| 国产白丝av| 国产麻豆精品手机在线观看| 九九九九热精品视频| 2021最新国产精品网站| 欧美日韩成人在线观看 | 乱系列中文字幕在线视频| 精品综合久久久久久97| 欧美影院久久| 国产迷奸在线看| 免费人成视网站在线不卡 | 2020精品极品国产色在线观看| 国产96在线 | 青青草原国产免费av观看| A级全黄试看30分钟小视频| 97se亚洲综合在线天天|