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

馬鈴薯超聲強化遠紅外輻射干燥特性及神經(jīng)網(wǎng)絡(luò)模型研究

2019-04-24 09:18:56席慧涵劉云宏李慧洋李功偉趙蘭馨
食品與機械 2019年2期
關(guān)鍵詞:模型

席慧涵 劉云宏,2 王 琦 李慧洋 李功偉 趙蘭馨

(1. 河南科技大學(xué)食品與生物工程學(xué)院,河南 洛陽 471023;2. 農(nóng)產(chǎn)品干燥裝備河南省工程技術(shù)研究中心,河南 洛陽 471023)

馬鈴薯是老少皆宜的蔬菜之一[1]。馬鈴薯含有豐富的多酚類、維生素、礦物質(zhì)等營養(yǎng)物質(zhì),具有抗衰老、抗氧化、補氣養(yǎng)血、健腦益智等功效[2-3]。然而,馬鈴薯儲存不當(dāng)會發(fā)生褐變甚至腐敗,從而造成營養(yǎng)成分大量損失,因此,將馬鈴薯進行干燥是目前提高其經(jīng)濟價值的措施之一[4]。熱風(fēng)干燥是較常見的馬鈴薯干燥方法[5],該方法雖然成本低廉、操作方便,但有干燥時間過長、熱效率低、營養(yǎng)成分損失嚴(yán)重等缺點。由于馬鈴薯內(nèi)部組織結(jié)構(gòu)致密、水分?jǐn)U散較難,干燥脫水速率主要取決于內(nèi)部傳熱傳質(zhì)過程,因此,采取有效措施來促進內(nèi)部傳質(zhì)和傳熱,減小水分?jǐn)U散阻力并改善熱量傳遞狀況,將有利于提高脫水速率、縮短馬鈴薯的干燥時間。

近年來,遠紅外輻射因具有顯著的熱效應(yīng),相關(guān)干燥技術(shù)得到了廣泛應(yīng)用。遠紅外輻射具有一定的穿透性,可實現(xiàn)內(nèi)部加熱,從而加快物料水分從內(nèi)部向外擴散與蒸發(fā)。有關(guān)紫薯[6]、紅棗片[7]、胡蘿卜[8]等遠紅外輻射干燥的研究表明,與傳統(tǒng)干燥技術(shù)比較,遠紅外輻射干燥技術(shù)具有加熱均勻、干燥速度快、產(chǎn)品品質(zhì)好等優(yōu)點[9-10]。此外,由于超聲波明顯的傳質(zhì)強化效應(yīng),超聲強化干燥技術(shù)逐漸成為果蔬干燥研究的熱點[11-13],其中,直觸超聲強化干燥技術(shù)尤為引人關(guān)注。直觸式超聲是將超聲波及其能量直接傳入物料內(nèi)部,從而促進水分傳遞及提高能量效率。Liu 等[14]、Sch?ssler等[15]探討了直觸式超聲在干燥過程中的強化作用,結(jié)果表明,接觸式超聲可有效縮短物料干燥時間,提高干燥速率。綜上可知,遠紅外輻射干燥技術(shù)具有改善物料內(nèi)部傳熱的優(yōu)點,直觸式超聲技術(shù)能夠促進內(nèi)部傳質(zhì)進程,因此,將超聲強化與遠紅外輻射干燥相結(jié)合,理論上可改善物料的傳熱與傳質(zhì)狀況,從而有效提高干燥速率。但目前未見直觸超聲強化遠紅外干燥馬鈴薯的研究報道。

物料含水率是體現(xiàn)干燥進程的重要參數(shù),也是調(diào)節(jié)干燥參數(shù)的主要依據(jù)。但因干燥過程中含水率具有非線性和時變性的特點,并且常規(guī)的模型預(yù)測精度不高,存在局限性[16-17]。近些年,神經(jīng)網(wǎng)絡(luò)因其具有較高的非線性映射能力、很強的適應(yīng)性[18],已廣泛應(yīng)用于農(nóng)產(chǎn)品干燥領(lǐng)域。Nadian等[19]建立神經(jīng)網(wǎng)絡(luò)模型來預(yù)測蘋果片在干燥中的色澤、水分變化。白竣文等[20]將神經(jīng)網(wǎng)絡(luò)應(yīng)用于南瓜的真空脈動干燥,預(yù)測干燥過程中的含水率。但目前未見將神經(jīng)網(wǎng)絡(luò)用于超聲強化遠紅外輻射干燥的研究報道。

本試驗擬將直觸式超聲強化技術(shù)與遠紅外輻射干燥聯(lián)合應(yīng)用,研究超聲功率及遠紅外輻射溫度對馬鈴薯片干燥特性及有效水分?jǐn)U散系數(shù)的影響,并采用Elman、RBF、BP 3種神經(jīng)網(wǎng)絡(luò)來建立馬鈴薯超聲強化遠紅外輻射干燥過程的數(shù)學(xué)模型,確定最佳模型構(gòu)建方法,以期為超聲強化遠紅外輻射干燥技術(shù)及馬鈴薯干燥加工研究提供理論參考。

1 材料與方法

1.1 材料

馬鈴薯:購于河南洛陽。采用105 ℃[14]烘箱法測量并計算得到新鮮馬鈴薯的初始干基含水率為(567±8)%。

1.2 儀器與設(shè)備

超聲強化遠紅外輻射干燥裝置:河南科技大學(xué)自制,具體結(jié)構(gòu)和參數(shù)見文獻[21];

切片機:SHQ-1型,德州天馬機械有限公司;

電子天平:Scout SE 型,美國OHAUS公司。

1.3 干燥方法

干燥開始前,將馬鈴薯洗凈去皮,切成直徑40 mm、厚度3 mm的薄片,每次取10片[(55±1) g ]進行干燥試驗。為有效抑制干燥時的酶促褐變,在干燥前,先進行沸水漂燙2 min的鈍酶殺青處理,取出瀝水后,平鋪在超聲振動盤上并放入超聲聯(lián)合遠紅外輻射干燥機中進行干燥。每隔30 min進行稱重,直至干燥結(jié)束。

在研究過程中,將干燥介質(zhì)溫度、風(fēng)速、輻射板距物料的距離分別固定為30 ℃、1 m/s、10 cm,并維持干燥初始物料量與物料厚度不變,分別設(shè)定紅外輻射板溫度為200,240,280,320 ℃,超聲功率為0,30,60 W,研究遠紅外輻射溫度和超聲功率對馬鈴薯干燥特性的影響。每組干燥試驗均重復(fù)3次。

1.4 干燥指標(biāo)計算

1.4.1 干基含水率 參照文獻[21]。

1.4.2 干燥速率 參照文獻[22]。

1.4.3 有效水分?jǐn)U散系數(shù) 參照文獻[23]。

1.5 神經(jīng)網(wǎng)絡(luò)模型建立的數(shù)據(jù)采集和預(yù)處理

對于本試驗的干燥過程,其主要參數(shù)是超聲功率、遠紅外輻射溫度和干燥時間,物料水分比(MR)從初始到干燥結(jié)束共得到162組數(shù)據(jù),隨機選取140組作為神經(jīng)網(wǎng)絡(luò)訓(xùn)練數(shù)據(jù),剩余22組作為預(yù)測數(shù)據(jù)。分別采用Elman、RBF、BP 3種神經(jīng)網(wǎng)絡(luò)來建立馬鈴薯超聲強化遠紅外輻射干燥過程的數(shù)學(xué)模型。在進行網(wǎng)絡(luò)訓(xùn)練前需先對數(shù)據(jù)進行歸一化處理,將其控制在[-1,1]之內(nèi),這是為了避免因數(shù)據(jù)量綱和范圍的差異影響擬合精度。

1.6 數(shù)據(jù)處理與分析

采用Origin 8.5統(tǒng)計軟件進行數(shù)據(jù)處理,利用MATLAB(R2014a)進行數(shù)學(xué)模型計算。選用決定系數(shù)R2和均方根誤差 RMSE作為擬合評價指標(biāo),R2越大、 RMSE越小,則擬合精度越高[20]。

2 結(jié)果與討論

2.1 馬鈴薯超聲強化遠紅外輻射干燥特性

在輻射溫度分別為200,240,280,320 ℃及超聲功率為0,30,60 W的條件下,進行馬鈴薯片干燥試驗,對應(yīng)的干燥曲線及速率曲線如圖1所示。

由圖1可知,在所有的輻射溫度下,提高超聲功率均能縮短干燥時間、提高干燥速率。以240 ℃輻射溫度為例,與未施加超聲的干燥時間(540 min)相比,施加超聲功率為30,60 W后的干燥時間分別縮短至420,330 min,減少幅度分別達到22.2%和38.9%,對應(yīng)的平均干燥速率分別為1.34%/min和1.70%/min,與0 W時的1.04%/min相比,提高幅度分別為28.8%和63.5%。可見,在馬鈴薯遠紅外輻射干燥過程施加超聲,能夠獲得顯著的強化效果,并且超聲功率越大,強化作用越顯著。由于超聲的高頻振動導(dǎo)致內(nèi)部組織快速、反復(fù)壓縮和擴張,使其產(chǎn)生大量的微細(xì)空泡,氣泡瞬間爆裂后產(chǎn)生強大的動能和壓縮能[24],這些能量使得原本與微細(xì)管壁緊密結(jié)合的水分子附著力降低,增強水分子的流動性并加快其擴散與遷移。另一方面,超聲也會影響物料內(nèi)部孔隙結(jié)構(gòu),增大微孔道的尺寸和連通度,可顯著提高水分轉(zhuǎn)移速率[25]。Ortuo等[26]研究了超聲技術(shù)在橘子皮對流干燥中的應(yīng)用,發(fā)現(xiàn)超聲可顯著提高干燥傳質(zhì)系數(shù)和水分?jǐn)U散能力,這與本研究所得的超聲可有效提升干燥速率的結(jié)論相符。超聲功率增加,物料內(nèi)部振動加強,超聲波能夠產(chǎn)生更強的機械效應(yīng)與空化效應(yīng),有利于水分湍動與擴散,從而能夠快速完成干燥,提高干燥速率。

圖1 不同超聲功率和輻射溫度下馬鈴薯片干燥的干燥曲線與干燥速率曲線Figure 1 Drying curves and drying rate curves of potato slices dried at different ultrasound powers and radiation temperatures

此外,由圖1的干燥速率曲線還可知,隨著水分含量的不斷減少,干燥速率呈下降趨勢,表明馬鈴薯干燥屬于降速干燥過程。隨著干燥的進行,不同超聲功率下的干燥速率差異在不斷縮小,在干燥末期,所有干燥速率曲線幾乎重合。上述現(xiàn)象表明,超聲波的強化效應(yīng)與水分含量密切相關(guān),且隨含水率的降低而減弱。這是由于在馬鈴薯干燥的初始階段,物料含水率較高,大量的自由水有利于超聲的傳遞,超聲波的衰減系數(shù)較小,內(nèi)部傳遞阻力較小,有利于超聲波穿入物料內(nèi)部并產(chǎn)生較好的強化效應(yīng)[14],并且超聲功率增大,相應(yīng)產(chǎn)生的空化效應(yīng)和機械效應(yīng)更強,進一步加快水分脫除速率。但隨著干燥的進行,物料內(nèi)部水分含量不斷減少,超聲波衰減加劇,用于強化傳質(zhì)的超聲能量利用率降低。

圖1還表明了遠紅外輻射溫度對馬鈴薯干燥時間及干燥速率的影響。在相同超聲功率下,提高遠紅外輻射溫度也能夠明顯地減少干燥時間、增大干燥速率。以超聲功率30 W 時為例,輻射溫度分別為200,240,280,320 ℃ 時,對應(yīng)的干燥時間分別為510,420,360,300 min,對應(yīng)的平均干燥速率為1.10%/min,1.34%/min,1.55%/min,1.87%/min。在遠紅外輻射干燥過程中,采用熱電偶測量馬鈴薯片表面溫度,在遠紅外輻射溫度為200,240,280,320 ℃時,最終產(chǎn)品溫度分別達到(51.0±1.5),(60.0±1.0),(66.0±2.0),(71.0±0.8) ℃。這表明提高輻射溫度會導(dǎo)致物料表面溫度的升高,有利于物料內(nèi)部水分子運動及水分蒸發(fā)。遠紅外輻射板發(fā)射的能量傳至物料內(nèi)部并轉(zhuǎn)化為熱能,致使水分子運動加劇,提高水分子的能量及水分流動性,從而促進水分?jǐn)U散。提高輻射板溫度會產(chǎn)生更多的輻射能量,物料對能量的吸收率隨之提高。因此增加遠紅外輻射溫度,有利于加快內(nèi)部的水分?jǐn)U散與汽化,顯著縮短干燥時間,提高產(chǎn)品速率。

從圖1中還可看出,超聲功率在不同輻射溫度條件下對馬鈴薯干燥速率影響不同,在未施加超聲時,輻射溫度分別為200,240,280,320 ℃對應(yīng)的平均干燥速率為0.89%/min,1.04%/min,1.25%/min,1.56%/min,當(dāng)施加60 W的超聲功率時,平均干燥速率分別提高到1.44%/min,1.70%/min,2.08%/min,2.34%/min,與未施加超聲時相比,提高幅度分別為61.80%,63.46%,66.40%,50.00%。上述現(xiàn)象表明,提高遠紅外輻射溫度在一定程度上可增強超聲的強化作用,但過高的加熱溫度又會削弱超聲強化作用。這是由于隨著遠紅外輻射溫度的升高,馬鈴薯片吸收的能量隨之增加,物料內(nèi)部水分蒸發(fā)擴散能力也隨之增強,并且超聲波發(fā)射的能量直接傳到物料內(nèi)部之后,能夠產(chǎn)生較強的空化效應(yīng)和機械效應(yīng),進一步加劇水分子運動,提升其擴散速率。然而過高的輻射溫度又會使物料中水分溫度過高,相應(yīng)的蒸汽分壓升高到一定程度[21],將會增加超聲在物料內(nèi)部產(chǎn)生微泡的難度,微泡爆破力也會有所下降。因此,過高的遠紅外輻射溫度雖然能夠?qū)崿F(xiàn)良好的傳熱,但在一定程度會削弱超聲的傳質(zhì)強化效果。

2.2 馬鈴薯干燥的有效水分?jǐn)U散系數(shù)(Deff值)

馬鈴薯在不同超聲功率和輻射溫度下的有效水分?jǐn)U散系數(shù)變化規(guī)律如圖2所示,Deff值為1.15×10-10~1.80×10-9m2/s。圖2中表明,Deff值隨遠紅外輻射溫度的升高而增大,以超聲功率為30 W時為例,輻射溫度200 ℃時,Deff值為1.36×10-10m2/s,升高輻射溫度分別到240,280,320 ℃ 時,對應(yīng)的Deff值分別為1.66×10-10,1.96×10-10,2.43×10-10m2/s,提升幅度分別為22.06%,44.12%,78.68%。由上述現(xiàn)象可知,遠紅外輻射溫度升高,產(chǎn)生熱效應(yīng)更強烈,因此加快物料表面和內(nèi)部的微細(xì)管道擴張并增加其數(shù)量,水分子擴散能力隨之增強,相應(yīng)的Deff值就會提高。這一現(xiàn)象與Sharma等[27]的研究結(jié)果相同。此外,從圖2中還可看出,Deff值隨超聲功率的增加而增大,以輻射溫度為200 ℃為例,施加30,60 W的超聲功率時,對應(yīng)的Deff值分別為1.36×10-10,1.96×10-10m2/s,與未施加超聲時相比,增加比率分別為18.26%,70.43%。這是由于增加超聲功率可顯著增強其產(chǎn)生的空化效應(yīng)和機械效應(yīng),2種效應(yīng)共同作用能夠有效降低傳質(zhì)阻力,增強物料內(nèi)部水分湍動并減弱物料組織結(jié)構(gòu)對水分的吸附力和束縛力,有利于促進水分流動與擴散,從而提高有效水分?jǐn)U散系數(shù)。

不同字母表示在95%的置信水平有顯著差異;小寫字母表示超聲功率,大寫字母表示輻射溫度

圖2 馬鈴薯超聲強化遠紅外輻射干燥的Deff值

Figure 2 TheDeffvalues of ultrasound strengthened far-infrared radiation drying of potato

2.3 馬鈴薯干燥過程的神經(jīng)網(wǎng)絡(luò)模型建立

2.3.1 Elman神經(jīng)網(wǎng)絡(luò)模型 Elman神經(jīng)網(wǎng)絡(luò)是一種反饋型網(wǎng)絡(luò),結(jié)構(gòu)一般為4層:輸入層、隱含層、承接層和輸出層[28]。將干燥過程中的3個重要參數(shù)作為輸入層的神經(jīng)元,物料水分比(MR)作為輸出層的神經(jīng)元。從理論上講,Elman神經(jīng)網(wǎng)絡(luò)隱含層神經(jīng)元數(shù)目可以任意選定,但是神經(jīng)元數(shù)目越多,可以解決更復(fù)雜的問題,相應(yīng)的精度和速度也會提高。經(jīng)過多次嘗試,隱層神經(jīng)節(jié)點數(shù)為30時有最佳效果。進行權(quán)值修正時采用BP算法,該算法訓(xùn)練函數(shù)為默認(rèn)值traingdx,選擇tansig作為輸入層與隱含層之間的傳輸函數(shù),purelin作為隱含層與輸出層之間的傳輸函數(shù),mse作為性能函數(shù)。按照上述結(jié)構(gòu)建立網(wǎng)絡(luò)模型,訓(xùn)練結(jié)果如圖3所示,網(wǎng)絡(luò)的目標(biāo)誤差達到0.000 5時,訓(xùn)練步數(shù)為205。

圖3 Elman神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練誤差曲線Figure 3 Training error curve of Elman neural network model

2.3.2 RBF神經(jīng)網(wǎng)絡(luò)模型 RBF神經(jīng)網(wǎng)絡(luò)(徑向基網(wǎng)絡(luò))是一種3層前向網(wǎng)絡(luò)。與其他前向型神經(jīng)網(wǎng)絡(luò)最大的區(qū)別在于,隱含層是采用局部響應(yīng)的高斯函數(shù)作為激勵函數(shù),正因如此,RBF網(wǎng)絡(luò)具有訓(xùn)練時間短,能夠逼近任意非線性函數(shù)的優(yōu)勢。RBF網(wǎng)絡(luò)中需要設(shè)定的參數(shù)只有分度密度Spread,其值越大,函數(shù)越平滑,但太大又會造成傳遞函數(shù)的作用域擴大到全局,誤差較大,喪失原本網(wǎng)絡(luò)局部收斂的優(yōu)勢。經(jīng)過多次嘗試和誤差對比,當(dāng)Spread=3時,網(wǎng)絡(luò)預(yù)測性能最好。

2.3.3 BP神經(jīng)網(wǎng)絡(luò)模型 BP神經(jīng)網(wǎng)絡(luò)是一種典型的多層前向型網(wǎng)絡(luò),網(wǎng)絡(luò)結(jié)構(gòu)與RBF網(wǎng)絡(luò)相同。BP網(wǎng)絡(luò)隱含層神經(jīng)元個數(shù)不是固定的,需采用經(jīng)驗公式和多次訓(xùn)練選取最優(yōu)值,根據(jù)Kolmogorov定理[29]和不斷檢驗,隱含層節(jié)點數(shù)為7時有最佳的訓(xùn)練效果。模型選用優(yōu)化的Levenberg-Marquardt算法[30],訓(xùn)練函數(shù)為trainlm,BP網(wǎng)絡(luò)模型中其他參數(shù)與Elman網(wǎng)絡(luò)保持一致。按照上述算法建立3-7-1結(jié)構(gòu)的BP神經(jīng)網(wǎng)絡(luò)模型,訓(xùn)練誤差曲線如圖4所示,網(wǎng)絡(luò)的目標(biāo)誤差達到0.000 5時,訓(xùn)練步數(shù)為151。

圖4 BP神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練誤差曲線Figure 4 Training error curve of BP neural

2.3.4 模型預(yù)測結(jié)果分析 訓(xùn)練完成后,將剩余的22組數(shù)據(jù)用于網(wǎng)絡(luò)預(yù)測,與試驗所得的水分比(MR)作對比,結(jié)果如圖5所示,可明顯看出BP神經(jīng)網(wǎng)絡(luò)的預(yù)測值與真實值更為接近,預(yù)測性能高于其他2個神經(jīng)網(wǎng)絡(luò)模型。采用決定系數(shù)(R2)和均方根誤差(RMSE)作為評價指標(biāo),對比結(jié)果如表1所示,3種神經(jīng)網(wǎng)絡(luò)的R2均在99%以上,預(yù)測性能較好,但使用優(yōu)化算法的BP神經(jīng)網(wǎng)絡(luò)擬合精度最高,R2=0.998 6,RMSE=0.006 5。因此,BP神經(jīng)網(wǎng)絡(luò)逼近能力更強,誤差更小,能夠很好地擬合馬鈴薯直觸超聲強化遠紅外輻射干燥過程。

圖5 不同神經(jīng)網(wǎng)絡(luò)預(yù)測值對比Figure 5 Comparison of prediction results of different neural networks

模型R2RMSERBF0.991 70.021 2Elman0.992 40.016 3BP0.998 60.006 5

3 結(jié)論

在馬鈴薯片超聲強化遠紅外輻射干燥過程中,超聲功率增加、遠紅外輻射溫度升高均可有效促進物料內(nèi)部的傳質(zhì)與傳熱,從而縮短干燥時間,加快干燥速率。隨著干燥的進行及水分含量逐漸減少,超聲波產(chǎn)生的強化效應(yīng)不斷減弱,表明水分減少會導(dǎo)致超聲傳播衰減增大。在遠紅外輻射干燥過程中,超聲強化效果隨溫度升高呈先增加后下降的趨勢,表明過高的遠紅外輻射溫度會減弱超聲的傳質(zhì)強化效果。提高超聲功率及遠紅外輻射溫度可增加Deff值,表明兩者均能促進干燥過程中的水分流動與擴散。針對不同超聲功率、遠紅外輻射溫度、干燥時間對馬鈴薯含水率的影響,分別構(gòu)建了Elman、RBF、BP 3種神經(jīng)網(wǎng)絡(luò)模型。3種神經(jīng)網(wǎng)絡(luò)模型均有較好的預(yù)測性能,其中使用優(yōu)化算法的BP網(wǎng)絡(luò)模型預(yù)測值與真實值最為接近,能夠快速準(zhǔn)確地進行馬鈴薯超聲強化遠紅外輻射干燥過程中含水率的預(yù)測。

本試驗研究結(jié)果證實了將超聲的傳質(zhì)強化和遠紅外輻射的良好傳熱性能相結(jié)合,可以實現(xiàn)優(yōu)勢互補,從而顯著加快干燥進程、縮短干燥時間。然而,超聲對遠紅外輻射干燥過程的傳質(zhì)強化機理與品質(zhì)保護機制需要進一步的研究。此外,雖然本研究利用神經(jīng)網(wǎng)絡(luò)實現(xiàn)了馬鈴薯超聲強化遠紅外輻射干燥過程的模擬,但所建模型并未從微觀角度來深入揭示干燥過程的質(zhì)熱傳遞機理,因此,超聲強化遠紅外輻射干燥的質(zhì)熱傳遞解析也是未來研究的重要方向。

猜你喜歡
模型
一半模型
一種去中心化的域名服務(wù)本地化模型
適用于BDS-3 PPP的隨機模型
提煉模型 突破難點
函數(shù)模型及應(yīng)用
p150Glued在帕金森病模型中的表達及分布
函數(shù)模型及應(yīng)用
重要模型『一線三等角』
重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
3D打印中的模型分割與打包
主站蜘蛛池模板: 亚洲无限乱码一二三四区| 欧美一级高清片久久99| 天天色天天操综合网| 日本亚洲欧美在线| 高清久久精品亚洲日韩Av| 国产福利影院在线观看| 97国产在线视频| 四虎在线观看视频高清无码| 国产亚洲第一页| 国产性爱网站| 国产剧情无码视频在线观看| 中文字幕在线一区二区在线| 日韩在线视频网站| 亚洲天堂网2014| 国产91九色在线播放| 在线免费无码视频| 一本大道无码日韩精品影视| 精品视频在线观看你懂的一区| 国产91丝袜| 丝袜国产一区| 午夜国产在线观看| 国产成人精品免费视频大全五级| 99在线视频网站| 亚洲精品久综合蜜| 亚洲综合经典在线一区二区| 国产成人精品一区二区三区| 国产经典免费播放视频| 久久国产精品77777| 国产污视频在线观看| 国产亚洲精品自在久久不卡| 99热免费在线| 国产在线一二三区| 一级福利视频| 国产主播喷水| 美女高潮全身流白浆福利区| 久久久久亚洲AV成人网站软件| 在线a网站| 精品少妇人妻一区二区| 国产91在线|中文| 亚洲资源站av无码网址| 高清乱码精品福利在线视频| 91精选国产大片| 国产精品自拍露脸视频| 欧美在线导航| 亚洲精品国产乱码不卡| 精品亚洲麻豆1区2区3区| 亚洲日本韩在线观看| 亚洲 欧美 偷自乱 图片| 国产欧美在线视频免费| 国产乱子伦手机在线| 国产内射在线观看| 亚洲精品图区| 狠狠操夜夜爽| 91欧美在线| 欧美成人午夜在线全部免费| 国产色网站| 日韩视频免费| 亚洲国产理论片在线播放| 伊人久久福利中文字幕| 高清无码不卡视频| 波多野结衣的av一区二区三区| 亚洲中文精品人人永久免费| 亚洲成A人V欧美综合| 久久精品无码一区二区日韩免费| 欧美精品成人一区二区视频一| 亚洲精品另类| 无码AV日韩一二三区| 精品成人一区二区| 亚洲一区二区三区国产精品 | 免费国产高清精品一区在线| 亚洲精品成人7777在线观看| 国产日韩欧美精品区性色| 欧美精品啪啪| 国产乱子伦一区二区=| 无码综合天天久久综合网| 久久人体视频| 国产打屁股免费区网站| 午夜人性色福利无码视频在线观看| 欧美日韩专区| 婷婷99视频精品全部在线观看| 91麻豆国产在线| 日韩午夜福利在线观看|