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

苦瓜微波干燥特性及動力學(xué)模型

2017-12-02 14:15:42田華韓艷婷
食品研究與開發(fā) 2017年23期
關(guān)鍵詞:模型研究

田華,韓艷婷

(1.信陽師范學(xué)院生命科學(xué)學(xué)院,河南信陽464000;2.信陽師范學(xué)院大別山農(nóng)業(yè)生物資源保護(hù)與利用研究院,河南信陽464000)

苦瓜微波干燥特性及動力學(xué)模型

田華1,2,韓艷婷1,2

(1.信陽師范學(xué)院生命科學(xué)學(xué)院,河南信陽464000;2.信陽師范學(xué)院大別山農(nóng)業(yè)生物資源保護(hù)與利用研究院,河南信陽464000)

采用單因素試驗(yàn)法,研究切片厚度、微波功率和裝載量對苦瓜微波干燥特性的影響,并建立苦瓜微波干燥動力學(xué)模型。試驗(yàn)結(jié)果表明:微波功率對苦瓜干燥影響最大,其次是裝載量,最后是切片厚度;苦瓜微波干燥分加速干燥階段和降速干燥階段。對7種常用的薄層干燥動力學(xué)數(shù)學(xué)模型擬合,通過比較相關(guān)系數(shù)R2、殘基平方和RSS和卡方χ2得出,tian model最適于描述苦瓜片微波薄層干燥過程,其模型系數(shù)在0.99以上。

薄層干燥;微波干燥;干燥特性;苦瓜;動力學(xué)模型

苦瓜,葫蘆科苦瓜屬植物,苦瓜中除含有水分、蛋白質(zhì)、脂肪、碳水化合物、粗纖維和各種維生素外,還含有苦瓜甙、類蛋白活性物質(zhì)及多種氨基酸和多種人體必需礦物元素,具有清熱解毒、養(yǎng)血益氣、補(bǔ)腎保脾、養(yǎng)肝明目的功效,是夏季上火、發(fā)熱、痢疾、腫瘡等疾病的食療良品[1-2]。新鮮的苦瓜含水量高達(dá)90%以上,不易于常溫儲藏,苦瓜的干燥是對苦瓜鮮品進(jìn)行加工貯藏的主要手段之一。目前對苦瓜進(jìn)行干燥的技術(shù)主要有

熱風(fēng)干燥[3-10]、微波干燥[5-7,10-11]、熱泵干燥[7]、熱風(fēng)結(jié)合微波干燥[5,10]、真空微波干燥[7]、冷凍干燥[8,12-13]、遠(yuǎn)紅外干燥[9]、噴霧干燥[14]。而微波干燥是一項(xiàng)節(jié)能、干燥均勻、干燥效率高、清潔生產(chǎn)、產(chǎn)品質(zhì)量高的干燥技術(shù)[15]。國內(nèi)外很早就有關(guān)于微波干燥的研究[16-20],目前微波干燥技術(shù)已在水果、番薯、蔬菜、茶葉等農(nóng)產(chǎn)品干燥領(lǐng)域中廣泛應(yīng)用[21-32]。目前苦瓜微波干燥主要集中在干燥工藝的優(yōu)化及干燥品質(zhì)的研究,干燥動力學(xué)模型構(gòu)建方面的研究還沒有報(bào)道。因此,本文利用微波技術(shù)手段,開展苦瓜傳熱傳質(zhì)特性研究,探討水分在苦瓜中的移動規(guī)律,建立苦瓜薄層微波干燥的動力學(xué)模型,為實(shí)現(xiàn)苦瓜干燥過程自動化、改善干燥質(zhì)量、提高干燥效率和降低干燥能耗奠定基礎(chǔ),以期為苦瓜薄層微波干燥工藝和生產(chǎn)控制提供理論依據(jù)。

1 材料與方法

1.1 材料與儀器

苦瓜:市區(qū)超市。

美的EM7KCG4-NR微波爐:廣東美的廚房電器制造有限公司;奧豪斯AR2202CN電子天平:美國奧豪斯儀器(上海)有限公司。

1.2 方法

1.2.1 苦瓜切片制作

挑選沒蟲害、無腐爛、外形整齊、大小相近的苦瓜,去蒂,清洗干凈,根據(jù)試驗(yàn)要求的厚度切片,護(hù)色,瀝干水分,均勻平鋪成一層在盛物皿上,分別在不同的條件下進(jìn)行微波干燥,每隔1 min測一次質(zhì)量,直至恒重。每組試驗(yàn)重復(fù)3次,取平均值。

1.2.2 微波干燥

將苦瓜片置于盛物器皿上,放在微波爐中心位置,對苦瓜片進(jìn)行干燥處理。干燥過程中每隔1 min用電子天平稱量苦瓜片的質(zhì)量,每次測量時(shí)間控制在25 s之內(nèi),干燥苦瓜片至恒重時(shí)為止,每次平行測定3組樣品。

1.2.2.1 微波干燥功率的選擇

用電子天平稱取3份苦瓜切片,每份50 g,切片厚度為0.5 cm,然后平鋪成一層放在盛物器皿上,放入微波爐中,分別在700、560、350 W功率下進(jìn)行微波干燥,每隔1 min取出測一次質(zhì)量,直至恒重。

1.2.2.2 微波干燥切片厚度的選擇

用電子天平稱取3份苦瓜切片,每份50 g,切片厚度分別為0.3、0.5、0.7 cm,然后平鋪成一層放在盛物器皿置微波爐中,然后在700 W功率下進(jìn)行微波干燥,每隔1 min取出測一次質(zhì)量,直至恒重。

1.2.2.3 微波干燥干燥量的選擇

用天平稱取3份苦瓜切片,每份質(zhì)量為30、50、70 g,切片厚度為0.5 cm,然后平鋪成一層放在盛物器皿上,放入微波爐中,然后在700 W功率下進(jìn)行微波干燥,每隔1 min取出測一次質(zhì)量,直至恒重。

1.2.3 干基含水率Mt、干燥速率Dr、水分比MR、有效水分?jǐn)U散系數(shù)(Deff)的測定。

干基含水率Mt的計(jì)算方法見公式(1):

式中:Mt為試樣干燥至t時(shí)刻的干基含水率,g水/g干料;mt、mg分別為試樣干燥至t時(shí)刻與干燥恒重時(shí)的質(zhì)量,g。

干燥速率Dr的計(jì)算方法見公式(2):

式中:Dr為干燥速率,g/(g·h);Mt和Mt+Δt分別為干燥至t和t+Δt時(shí)刻干基含水率,g水/g干料;Δt干燥間隔時(shí)間,h。

水分比(MR)用于表示一定干燥條件下,t時(shí)刻樣品的水分含量與初始樣品水分含量的比值,可以用來反映物料干燥速率的快慢,計(jì)算方法見公式(3):

式中:Mt為試樣干燥至t時(shí)刻的的干基含水率,g水/g干料;M0為初始水分含量,g水/g干料;Me為平衡水分含量,g水/g干料,通常Me可以忽略。

有效水分?jǐn)U散系數(shù)(Deff),反映物料在一定干燥條件下的脫水能力,對深入分析物料內(nèi)部水分?jǐn)U散過程及優(yōu)化干燥工藝具有重要意義,計(jì)算方法見公式(4):

式中:Deff為有效水分?jǐn)U散系數(shù),m2/s;L為苦瓜切片厚度,m;t為干燥時(shí)間,s。當(dāng)物料進(jìn)行長時(shí)間干燥時(shí),公式(4)可簡化為:

1.3 薄層干燥模型

物料干燥是一個(gè)復(fù)雜的傳熱傳質(zhì)過程,本試驗(yàn)采用7種常用的薄層干燥模型對苦瓜片的干燥過程進(jìn)行描述,其表達(dá)式如表1所示。對數(shù)學(xué)模型①-⑥進(jìn)行線性處理后,數(shù)學(xué)模型①Newton、⑤Henderson and Pabis、⑥Logaritlunic顯示lnMR-t呈線性,②page模型顯示ln(-lnMR)-lnt呈線性,③Two term、⑦tian model顯示MR-t呈指數(shù)方程,④Wang and singh模型顯示MR-t呈二次多項(xiàng)式回歸,其中⑦tian model由本文擬合構(gòu)建。

表1 常用的薄層干燥數(shù)學(xué)模型Table 1 Commonly used mathematical models of thin-layer drying

2 結(jié)果與討論

2.1 切片厚度、微波功率、裝載量對苦瓜干燥速率Dr的影響

圖1~3分別是切片厚度、微波功率、裝載量對苦瓜干燥速率的影響。

圖1 切片厚度對Dr的影響(700 W、50 g)Fig.1 Effect of slice thickness on drying rate

圖2 微波功率對Dr的影響(0.5 cm、50 g)Fig.2 Effect of microwave power on drying rate

圖3 裝載量對Dr的影響(700 W、0.5 cm)Fig.3 Effect of loadage on drying rate

由圖1~圖3可以看出,在不同的切片厚度下,苦瓜微波干燥速率先升后降,干燥過程中最大干燥速率分別為208.085 g/(g·h)(0.3 cm)、178.957 g/(g·h)(0.5 cm)、162.087 g/(g·h)(0.7 cm)。在不同微波功率下,苦瓜微波干燥速率先升后降,干燥過程中最大干燥速率分別為178.957 g/(g·h)(700 W)、162.078 g/(g·h)(560 W)、87.173 g/(g·h)(350 W)。在不同的裝載量下,苦瓜微波干燥速率表現(xiàn)不一致,30、50 g裝載量下干燥速率先升后降,70 g裝載量下干燥速率持續(xù)下降,干燥過程中最大干燥速率分別為264.808 g/(g·h)(30 g)、178.957 g/(g·h)(50 g)、108.529 g/(g·h)(70 g)。由以上分析可以看出,微波功率對苦瓜干燥影響最大,其次是裝載量,最后是切片厚度。苦瓜微波干燥分加速干燥階段和降速干燥階段。

2.2 切片厚度、微波功率、裝載量對苦瓜對Deff的影響

圖4~6分別是切片厚度、微波功率、裝載量對苦瓜對Deff的影響。

圖4 切片厚度對Deff的影響(700 W、50 g)Fig.4 Effect of slice thickness on effective water diffusivity

圖5 微波功率對Deff的影響(0.5 cm、50 g)Fig.5 Effect of microwave power on effective water diffusivity

圖6 裝載量對Deff的影響(700 W、0.5 cm)Fig.6 Effect of loadage on effective water diffusivity

由圖4~圖6可以看出,隨著干燥時(shí)間的延長,Deff逐漸增加,微波功率和裝載量對苦瓜干燥Deff影響較大,切片厚度對苦瓜干燥Deff影響較小,這與微波功率、裝載量對苦瓜微波干燥速率的影響表現(xiàn)一致,0.3 cm或0.5 cm、30 g、700 W是最佳微波干燥條件。

3 苦瓜微波干燥動力學(xué)模型構(gòu)建

建立干燥模型對研究干燥規(guī)律、預(yù)測干燥工藝參數(shù)有重要作用。目前,用來描述薄層干燥過程的模型一般有7種,見表1。鑒于圖1~6,不同切片厚度、微波功率、裝載量對苦瓜微波干燥速率Dr、Deff的影響,0.3 cm或0.5 cm、30 g、700 W是最佳微波干燥條件。于是,選取切片厚度0.5 cm、裝載量50 g、微波功率700 W作為苦瓜微波干燥動力學(xué)研究的參數(shù)條件。其中裝載量選擇50 g而不是30 g,這是因?yàn)檠b載量為30 g時(shí)微波干燥到恒重的時(shí)間為6 min,參數(shù)較少,影響擬合方程的準(zhǔn)確性。因此,分別制作切片厚度為0.5 cm、載物量50 g、微波功率為700 W苦瓜微波干燥的 MR-t曲線、lnMR-t曲線、ln(-lnMR)-lnt曲線。

圖7是苦瓜微波干燥MR-t二次擬合曲線,回歸方程R2=0.997 5,符合Wang and singh模型。圖8是苦瓜微波干燥lnMR-t關(guān)系曲線,回歸方程R2=0.958 3,不符合表1所列薄層干燥數(shù)學(xué)模型;圖9是苦瓜微波干燥 ln(-lnMR)-lnt關(guān)系曲線,回歸方程 R2=0.997 7,符合Page模型。圖10是苦瓜微波干燥MR-t指數(shù)方程擬合,回歸方程R2=0.998 11,符合Two term和tian model模型。故選擇Two term方程、Page方程、Wang and singh方程和tian model作為苦瓜微波干燥動力學(xué)模型。為了確定最佳模型,對origin得到的統(tǒng)計(jì)數(shù)據(jù)進(jìn)行對比,如表2所示。

圖7 MR-t二次擬合曲線Fig.7 Two times fitting curve of MR and t

為得到適于微波干燥苦瓜片的數(shù)學(xué)模型,將MR-t關(guān)系圖中各曲線分別帶入模型中,處理數(shù)據(jù)得到苦瓜干燥的R2值。一般來說R越高,χ2值和RSS值越小,說明擬合度越高,準(zhǔn)確性也越高。由表2可以看出,tian model方程擬合得到方程的R2為0.999,χ2最小,RSS最小,說明tian model方程最適于描述苦瓜片微波薄層干燥過程。

圖8 lnMR-t關(guān)系曲線Fig.8 Relationship curve of lnMR and t

圖 9 ln(-lnMR)-lnt關(guān)系曲線Fig.9 Relationship curve of ln(-lnMR)and ln t

圖10 MR-t指數(shù)擬合曲線Fig.10 Exponential fitting curve of MR and t

4 結(jié)論

以苦瓜片為原料,研究了樣品經(jīng)過不同微波功率(700、560、350 W)、不同切片厚度(0.3、0.5、0.7 cm)、和不同裝載量(30、50、70 g)干燥過程的動力學(xué)變化,為微波干燥技術(shù)在苦瓜片的深加工應(yīng)用技術(shù)研究提供基礎(chǔ)。通過考察苦瓜片在不同微波功率下的干燥特性,建立微波干燥的水分含量比MR與干燥時(shí)間關(guān)系的數(shù)學(xué)模型,得出以下結(jié)論:苦瓜干燥分加速干燥階段和降速干燥階段。微波功率對苦瓜干燥影響最大,其次是裝載量,最后是切片厚度。對7種常用的薄層干燥動力學(xué)數(shù)學(xué)模型擬合,通過比較相關(guān)系數(shù)R2、殘基平方和RSS和卡方χ2得出,tian model最適于描述苦瓜片微波薄層干燥過程,其模型系數(shù)在0.99以上。

表2 苦瓜微波薄層干燥模型統(tǒng)計(jì)結(jié)果Table 2 Statistical results obtained from microwave thin-layer drying models of balsam pear

[1]董義,王斯慧,曾里,等.六種保健食品原料的輔助降血糖功效研究[J].食品研究與開發(fā),2014,35(23):121-124

[2]李越剛.苦瓜養(yǎng)生之我見[J].中國果菜,2016,36(5):78-79

[3]王華紅.苦瓜切片熱風(fēng)干燥特性及含水率預(yù)測模型研究[D].呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學(xué),2013

[4]王華紅,趙士杰,張博.苦瓜切片熱風(fēng)干燥試驗(yàn)研究[J].農(nóng)機(jī)化研究,2013,35(7):179-182

[5]聶凌鴻,張偉.苦瓜切片干制及其復(fù)水性的研究[J].食品工業(yè),2015(2):185-187

[6]聶凌鴻,張偉.苦瓜切片干燥工藝的研究[J].食品研究與開發(fā),2016,37(2):48-52

[7]鄒宇曉,徐玉娟,張名位,等.苦瓜的干制工藝研究[J].廣東農(nóng)業(yè)科學(xué),2006(11):62-63

[8]江英,胡小松,生吉萍.不同干燥方法對苦瓜粉品質(zhì)的影響[J].食品研究與開發(fā),2002,23(2):53-54

[9]朱香燕,張珺,何義雁,等.熱風(fēng)與遠(yuǎn)紅外干燥溫度對苦瓜全粉品質(zhì)的影響[J].現(xiàn)代食品科技,2015,31(7):265-269,325

[10]唐小俊,池建偉,張名位,等.熱風(fēng)結(jié)合微波干燥苦瓜片的工藝研究[J].食品科技,2008,34(3):76-81

[11]唐小俊,池建偉,張名位,等.苦瓜微波干燥工藝優(yōu)化[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2008,39(1):68-73

[12]吳茂玉,葛邦國,朱風(fēng)濤,等.苦瓜凍干粉生產(chǎn)技術(shù)研究[J].食品工業(yè)科技,2008(10):156-159

[13]祝瑩.苦瓜凍干超微粉對代謝綜合征的干預(yù)作用研究[D].鎮(zhèn)江:江蘇大學(xué),2013

[14]張敏,何俊萍,趙永會.噴霧干燥生產(chǎn)苦瓜復(fù)合粉工藝研究[J].飲料工業(yè),2011,12(6):34-37

[15]田加勝.微波技術(shù)基礎(chǔ)[M].武漢:華中科技大學(xué)出版社,2011:1-2

[16]Gomes A M,Leonhardt G F,Torloni M,et al.Microwave drying of microorganisms:Influence of the microwaveenergy and of the sample thicknessonthedryingofyeast[J].Journal of Microwave Power,1975,10(3):265-270

[17]Mothiron J C,Vialard-Goudou A,Maupas P.Proceedings:Use of an industrial prototype machine for drying and sterilizing pharmaceutical ampoules by microwaves at 2.45 GHz[J].The Journal of microwave power,1976,11(2):200-201

[18]M Tsang,Mui Chung,Sheldon C.Furutani.Microwave drying of macadamianuts[J].AmericanSocietyofAgriculturalEngineers,1989,5(4):565-567

[19]李耘庭,王嘉文,張鋼,等.八一糖廠應(yīng)用微波干燥方糖的初步試驗(yàn)[J].甜菜糖業(yè),1978(8):1-9

[20]李耘庭,王嘉文,張鋼,等.方塊糖生產(chǎn)應(yīng)用微波干燥初步試驗(yàn)[J].甘蔗糖業(yè),1979(3):10-15

[21]朱德泉,王繼先,朱德文.玉米微波干燥特性及其對品質(zhì)的影響[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2006,37(2):72-75

[22]孔慶新,祝冬青.微波干燥法生產(chǎn)南瓜粉的研究[J].食品工業(yè)科技,2006(9):150-152

[23]呂為喬,韓清華,李樹君,等.微波干燥姜片模型建立與去水機(jī)理分析[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2015,46(4):233-237

[24]段振華,汪菊蘭.微波干燥技術(shù)在食品工業(yè)中的應(yīng)用研究[J].食品研究與開發(fā),2007,28(1):155-158

[25]趙超,陳建,邱兵,等.花椒微波干燥特性試驗(yàn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2007,38(3):99-101

[26]蔣玉萍,王俊.番薯片微波干燥特性及干燥模型[J].浙江農(nóng)業(yè)學(xué)報(bào),2009,21(4):407-410

[27]余煉,顏棟美,侯金東.牡蠣微波干燥特性及動力學(xué)研究[J].食品科學(xué),2012,33(11):111-115

[28]朱國旺.金銀花微波干燥工藝研究及環(huán)形隧道腔式微波干燥設(shè)備的研制[D].鄭州:河南農(nóng)業(yè)大學(xué),2007

[29]許綽微,張璐,趙思明,等.發(fā)芽糙米熱風(fēng)和微波干燥特性及品質(zhì)研究[J].現(xiàn)代食品科技,2016,32(1):256-260

[30]林向陽,張麗晶,Roger Ruan,等.茶葉真空微波干燥特性及動力學(xué)模型[J].中國農(nóng)學(xué)通報(bào),2010,26(22):65-70

[31]任廣躍,化春光,段續(xù),等.鮮切懷山藥片微波干燥特性及品質(zhì)研究[J].食品科學(xué),2010,31(22):203-206

[32]Ernest Ekow Abano.熱風(fēng)干燥、真空微波干燥和紅外干燥對番茄片干燥動力學(xué)和品質(zhì)的影響[D].鎮(zhèn)江:江蘇大學(xué),2013

Microwave Drying Characteristics and Dynamic Model of Balsam Pear

TIAN Hua1,2,HAN Yan-ting1,2
(1.College of Life Science,Xinyang Normal University,Xinyang 464000,Henan,China;2.Key Laboratory of Simulation and Control for Dabie Mountains Population Ecology,Xinyang Normal University,Xinyang 464000,Henan,China)

In order to to study the influence of slice thickness,microwave power and loadage on the microwave drying characteristics and dynamic model of balsam pear,the single factor experiment was investigated.Results showed that microwave power had the greatest influence on the dryness of balsam pear,followed by loadage,and finally the thickness.Balsam pear drying had acceleration and reduction stages.The results of comparing the values of the correlation coefficient(R2),residual sum of squares(RSS)and reduced chi-square(χ2)showed that tian model was the best one to describe the relationship of moisture content and drying time during microwave drying,and the coefficient was above 0.99.

thin-layer drying;microwave drying;drying characteristics;balsam pear;dynamic model

10.3969/j.issn.1005-6521.2017.23.023

信陽師范學(xué)院大別山農(nóng)業(yè)生物資源保護(hù)與利用研究院開放課題

田華(1979—),女(漢),副教授,博士,研究方向:食品保藏及營養(yǎng)評價(jià)。

2017-04-19

猜你喜歡
模型研究
一半模型
FMS與YBT相關(guān)性的實(shí)證研究
2020年國內(nèi)翻譯研究述評
遼代千人邑研究述論
重要模型『一線三等角』
重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
視錯(cuò)覺在平面設(shè)計(jì)中的應(yīng)用與研究
科技傳播(2019年22期)2020-01-14 03:06:54
EMA伺服控制系統(tǒng)研究
新版C-NCAP側(cè)面碰撞假人損傷研究
3D打印中的模型分割與打包
主站蜘蛛池模板: 一级一级特黄女人精品毛片| 中文精品久久久久国产网址 | 日本精品影院| 免费一极毛片| 日本日韩欧美| 自拍中文字幕| 国产在线98福利播放视频免费| 熟妇丰满人妻av无码区| 亚洲一区二区精品无码久久久| 亚洲中文字幕97久久精品少妇| 色婷婷成人| 国产精品伦视频观看免费| 一级黄色网站在线免费看| 国产亚洲男人的天堂在线观看| 免费人成视频在线观看网站| 日本精品αv中文字幕| av一区二区无码在线| 高h视频在线| 日本欧美午夜| 亚洲午夜综合网| 欧美成人看片一区二区三区| 国产综合无码一区二区色蜜蜜| 91久久性奴调教国产免费| 成人综合网址| 精品久久久久久久久久久| 亚洲成a∧人片在线观看无码| 国产av无码日韩av无码网站| 久热中文字幕在线| 亚洲精品免费网站| 四虎永久在线| 国产在线视频欧美亚综合| 人人看人人鲁狠狠高清| 色久综合在线| 亚洲色图在线观看| 97se亚洲综合在线韩国专区福利| 无套av在线| 色网站在线视频| 亚洲第一成年网| 超清人妻系列无码专区| 午夜福利视频一区| 国产主播一区二区三区| 国产精品白浆在线播放| 香蕉综合在线视频91| 在线看片中文字幕| 欧美成一级| 91口爆吞精国产对白第三集| 456亚洲人成高清在线| 99久久精品视香蕉蕉| 国产永久无码观看在线| 免费va国产在线观看| 亚洲成人精品| 国产精品区视频中文字幕| 亚洲国产天堂在线观看| 国产永久在线观看| 欧美精品影院| 亚洲人成人无码www| 国产91导航| 四虎影视国产精品| 99热国产在线精品99| 免费jjzz在在线播放国产| 免费无码又爽又黄又刺激网站 | 91毛片网| 精品无码国产一区二区三区AV| 日本伊人色综合网| 精品无码国产一区二区三区AV| 日本不卡免费高清视频| 99精品伊人久久久大香线蕉| 日本伊人色综合网| 免费亚洲成人| YW尤物AV无码国产在线观看| 亚洲成在线观看| 国产精品美女网站| 国产免费观看av大片的网站| 国产小视频a在线观看| 激情综合网址| 67194亚洲无码| 色爽网免费视频| 午夜影院a级片| 久久人体视频| 囯产av无码片毛片一级| 久久夜色撩人精品国产| 夜精品a一区二区三区|