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

    有機(jī)溶劑對(duì)液體活性染料分散和水解穩(wěn)定性影響的模擬研究

    2021-03-06 09:37:46徐保律毛志平
    紡織學(xué)報(bào) 2021年2期
    關(guān)鍵詞:靜電勢(shì)吡咯烷酮活性染料

    徐保律, 吳 偉, 鐘 毅,3, 徐 紅,3, 毛志平,3,4

    (1. 東華大學(xué) 生態(tài)紡織教育部重點(diǎn)實(shí)驗(yàn)室, 上海 201620; 2. 東華大學(xué) 化學(xué)化工與生物工程學(xué)院,上海 201620; 3. 東華大學(xué) 紡織科技創(chuàng)新中心, 上海 201620; 4. 東華大學(xué) 國家染整工程技術(shù)研究中心, 上海 201620)

    液體染料比粉狀染料在產(chǎn)業(yè)化生產(chǎn)時(shí)所用添加劑的量少,有助于降低企業(yè)生產(chǎn)成本。隨著電子配色和染色自動(dòng)化工藝的發(fā)展,以及企業(yè)對(duì)清潔生產(chǎn)的要求,使用過程中無粉塵污染,對(duì)自動(dòng)化染色工藝適用性較好的液體染料在噴墨印花等領(lǐng)域得到廣泛的應(yīng)用,具有良好的發(fā)展前景,代表著現(xiàn)代染料工業(yè)的發(fā)展方向[1-2]。液體染料作為將原粉狀商品染料溶解在某種溶劑中使其液體化的一種染料,需要有較好的溶解度及穩(wěn)定性[3]。對(duì)于液體活性染料來說:一是如何提高染料的分散穩(wěn)定性,增加染料溶解度;二是如何抑制染料在水相條件下的水解[4]。作為在液體活性染料的配制中添加量通常僅次于水的物質(zhì),有機(jī)溶劑對(duì)提高液體活性染料的穩(wěn)定性具有重要的意義[5]。

    液體染料體系中,染料的大量聚集通常意味著沉淀的產(chǎn)生,伴隨著染料分散穩(wěn)定性的顯著變化。近年來,關(guān)于染料聚集的研究已經(jīng)很多,對(duì)有機(jī)溶劑影響液體染料的分散穩(wěn)定性也有大量的文獻(xiàn)報(bào)告。如將紫外-可見吸收光譜、流變學(xué)和表面張力等測(cè)試手段相結(jié)合,研究C.I.活性紅218與內(nèi)酰胺類化合物(2-吡咯烷酮和N-甲基吡咯烷酮)和聚乙二醇的相互作用,發(fā)現(xiàn)有機(jī)溶劑分子可以通過疏水作用和分離作用將C.I.活性紅218分子的聚集體轉(zhuǎn)變?yōu)檩^小的聚集基團(tuán)[6-7]。Park等[8-10]通過理論計(jì)算染料與水溶性聚合物聚乙烯吡咯烷酮(PVP)的結(jié)合常數(shù)和染料在聚合物水溶液中的聚集常數(shù)證明了PVP存在時(shí)對(duì)染料的解聚作用。目前文獻(xiàn)所報(bào)道關(guān)于有機(jī)溶劑對(duì)活性染料分散穩(wěn)定性的提高,主要是由于有機(jī)溶劑分子可以打破染料疏水基團(tuán)周圍的冰狀結(jié)構(gòu)水,從而靠近染料疏水基團(tuán),使染料疏水部分親水以防止染料分子聚集,但活性染料的高濃度使染料分子在溶液中處于更加復(fù)雜的狀態(tài)[11],實(shí)驗(yàn)也只能以不同的角度從側(cè)面反映染料混合溶液的部分性能,無法直接解釋有機(jī)溶劑提高染料分散穩(wěn)定性的機(jī)制。并且針對(duì)液體活性染料的水解穩(wěn)定性,大都為額外添加pH值緩沖溶液,而對(duì)有機(jī)溶劑影響活性染料水解穩(wěn)定性的詳細(xì)報(bào)道并不多。

    近年來,分子動(dòng)力學(xué)(MD)模擬作為一種重要的研究方法應(yīng)用于科學(xué)研究領(lǐng)域,是連接物質(zhì)的微觀結(jié)構(gòu)與宏觀現(xiàn)象的重要橋梁。MD模擬手段已經(jīng)被廣泛應(yīng)用于研究聚合物在水中的行為(溶解性、溫敏性、自組裝等),但對(duì)于MD模擬在活性染料應(yīng)用方面的研究報(bào)道相對(duì)較少。Wu等[12]通過分子動(dòng)力學(xué)模擬研究了3種主要活性染料在水中的聚集行為,確定了系統(tǒng)中染料的最大聚集尺寸不大于五聚體,以圖形方式顯示了染料-染料和染料-水分子間相互作用的位置、強(qiáng)度和類型。

    液體活性染料體系中,分子間相互作用通過實(shí)驗(yàn)手段難以檢測(cè)和量化,而通過分子動(dòng)力學(xué)模擬可以系統(tǒng)地進(jìn)行計(jì)算,因此,本文的重點(diǎn)是在參考以往研究的基礎(chǔ)上,采用全原子動(dòng)力學(xué)模擬的手段,結(jié)合實(shí)驗(yàn)測(cè)試,探究有機(jī)溶劑改善液體活性染料分散及水解穩(wěn)定性的機(jī)制,為配制液體活性染料更好地選擇合適的有機(jī)溶劑提供理論支持。

    1 實(shí)驗(yàn)部分

    1.1 實(shí)驗(yàn)材料

    C.I.活性藍(lán)176(商品染料,江蘇德美科化工有限公司);N-甲基吡咯烷酮(NMP)、己內(nèi)酰胺(CPL)、一縮二乙二醇(DEG)、二甲基亞砜(DMSO)、N,N-二甲基甲酰胺(DMF)、三氯甲烷,分析純,均為上海泰坦科技股份有限公司。

    1.2 實(shí)驗(yàn)方法

    1.2.1 染料的提純

    稱取定量商品染料置于燒杯中,加入DMF將其完全溶解,過濾。向?yàn)V液中加入適量三氯甲烷使染料結(jié)晶析出,靜置沉淀12 h,過濾,烘干即得提純?nèi)玖瞎腆w。

    1.2.2 活性染料溶液的配制

    用MS7-H550-Pro型恒溫磁力攪拌器在25 ℃下攪拌,將活性染料在去離子水和有機(jī)溶劑混合溶液中溶解2 h。其中染料和有機(jī)溶劑的質(zhì)量分?jǐn)?shù)分別為10%和15%。

    1.2.3 有機(jī)溶劑影響染料穩(wěn)定性的測(cè)定

    1.2.3.1紫外-可見吸收光譜的測(cè)量 使用U-3310型紫外-可見分光光度計(jì)(日本日立公司)測(cè)定不同體系的染料溶液在300~700 nm波長范圍內(nèi)的吸收光譜。

    1.2.3.2染料粒徑測(cè)定 將測(cè)試樣品經(jīng)砂芯漏斗過濾除雜,而后使用Nano-ZS型納米粒度儀(英國馬爾文公司)測(cè)試染液顆粒粒徑,樣品粒徑數(shù)據(jù)均采用平均粒徑值,所有樣品的數(shù)據(jù)選取3次測(cè)試結(jié)果的平均值。

    1.2.3.3染料的沉淀質(zhì)量 配制的染料溶液在25 ℃ 恒溫條件下靜置72 h后,經(jīng)離心分離機(jī)在3 000 r/min 條件下分離30 min后過濾,沉淀烘干后稱取質(zhì)量。

    1.2.3.4染料水解的高效液相色譜測(cè)定 配制的不同體系的染料溶液在25 ℃恒溫振蕩水浴鍋中水解反應(yīng)12 h,應(yīng)用反相離子對(duì)高效液相色譜對(duì)染液樣品進(jìn)行分析。色譜柱:Hypersil ODS,5 μm, 416 mm×250 mm的C18柱;流動(dòng)相:離子對(duì)緩沖溶液(含0.002 mol四丁基溴化胺和 0.01 mol 醋酸銨)與乙腈的混合液(體積比為40∶60);流速為1 mL/min; 檢測(cè)波長為600 nm;溫度為35 ℃。

    圖1 染料及有機(jī)溶劑的化學(xué)結(jié)構(gòu)Fig.1 Chemical structures of dye and organic solvents

    2 分子動(dòng)力學(xué)模擬

    活性藍(lán)176染料(RB176)與有機(jī)溶劑分子的化學(xué)結(jié)構(gòu)如圖1所示。使用GROMACS程序(版本2018.2)[13]對(duì)活性染料有機(jī)溶劑-水溶液進(jìn)行了模擬和分析。詳細(xì)的模擬設(shè)置如表1所示。包括溶液的混合比例、染料分子個(gè)數(shù)和三維周期系統(tǒng)的尺寸,其中不同的有機(jī)溶劑體系具有相同的染料質(zhì)量百分?jǐn)?shù)、總?cè)玖蠑?shù)和體系尺寸。

    表1 體系的模擬設(shè)定Tab.1 Simulation setting of systems

    染料陰離子的力場是GAFF力場[14],通常用于研究有機(jī)小分子,還使用了與之兼容的TIP3P水模型[15]。為確保模擬的準(zhǔn)確性,參考Kashefolgheta等[16]的研究,使用了Na+的Lennard-Jones(LJ)參數(shù)以及Na+與染料的磺酸鹽或硫酸鹽的氧原子之間的非鍵相互作用參數(shù)。對(duì)5種染料溶液體系進(jìn)行預(yù)平衡,并在恒溫恒壓下模擬200 ns。用Nose-hoover耦合方法[17]將溫度控制在298.15 K,并使用Parrinello-Rahman方法[18]將壓力耦合到101 kPa。所有的鍵約束都是通過LINCS算法[19]實(shí)現(xiàn),計(jì)算范德華相互作用的截?cái)喟霃綖?.2 nm,靜電相互作用采用PME方法[20]計(jì)算。牛頓方程采用蛙跳式Verlet算法[21]進(jìn)行積分,時(shí)間步長為2 fs。在所有方向上使用周期性邊界條件。模擬完成后,使用GROMACS程序和VMD程序(版本1.9.3)[22]進(jìn)行分析和可視化演示。

    2.1 分子表面靜電勢(shì)

    化學(xué)中的一個(gè)基本概念是:大多數(shù)分子具有相對(duì)富含電子(即電子層)部分,而其他部分則相應(yīng)地為電子貧乏(即原子核)部分。在原子核附近,由于離核較近,靜電勢(shì)為正值,而在分子表面,原子核與電子的貢獻(xiàn)可以相互平衡,不均勻的電子密度分布會(huì)導(dǎo)致分子表面靜電勢(shì)正負(fù)不一。當(dāng)在一個(gè)分子的外部區(qū)域評(píng)估靜電勢(shì)時(shí),表明該分子是如何被接近的反應(yīng)物“看見”的,因此對(duì)于該分子的行為,特別是在非共價(jià)相互作用中,分子表面靜電勢(shì)的計(jì)算很有意義。

    分子表面靜電勢(shì)(ESP)的計(jì)算公式為

    式中:V(r)為分子表面位置r處的靜電勢(shì), kJ/mol;r′為r附近位置;ZA為原子核A的電荷,C;RA為原子核A的位置;|RA-r|代表其到r的距離;ρ(r′)為r′位置處的電子密度。

    2.2 徑向分布函數(shù)

    徑向分布函數(shù)(RDF)的計(jì)算公式為

    RDF量化了在給定的A原子的r距離的球?qū)犹幈唤y(tǒng)計(jì)的B原子的數(shù)密度<ρB(r) >與平均數(shù)密度<ρB>local的比值,也即從A原子到一個(gè)給定的距離r找到B原子的概率。NA、NB分別為體系中A、B原子的總個(gè)數(shù);δ(rij-r)為接近r附近的微小厚度。

    2.3 空間分布函數(shù)

    在RDF的基礎(chǔ)上,可以將分布函數(shù)的概念推廣到三維,三維空間的分布函數(shù)也即空間分布函數(shù)(SDF)是粒子在體系中各個(gè)位置的局部數(shù)密度與平均數(shù)密度之比:

    通過SDF可以量化染料分子周圍的有機(jī)溶劑分子堆積,從而可視化系統(tǒng)內(nèi)染料基團(tuán)周圍有機(jī)溶劑分子的分離度,擴(kuò)展了對(duì)體系結(jié)構(gòu)以及染料與溶劑相互作用的常規(guī)分析。

    3 模擬結(jié)果分析

    3.1 染料陰離子的溶劑可及表面積分析

    溶劑可及表面積(SASA)的概念是由Lee和Richards在1971年首次提出,此后常被用于判斷有機(jī)小分子(如染料陰離子)模擬系統(tǒng)的平衡,更適用于估計(jì)聚集-分解平衡[12]。圖2示出不同體系中染料陰離子的SASA變化。在初始狀態(tài)下,72個(gè)未聚集的染料陰離子在所有體系中的SASA值都在700 nm2左右,因此,本文認(rèn)為當(dāng)染料陰離子SASA在一段時(shí)間內(nèi)的變化不超過35 nm2時(shí),體系便達(dá)到了相對(duì)動(dòng)態(tài)平衡,此條件下大多數(shù)染料陰離子的聚集狀態(tài)沒有發(fā)生變化[12]。

    圖2 模擬過程中不同體系染料陰離子的溶劑可及表面積變化Fig.2 SASA changes of dye anions in different systems during simulation process

    圖2中5個(gè)體系染料陰離子的SASA均隨時(shí)間變化而逐漸降低,并趨于平緩,取最后50 ns作為相對(duì)平衡時(shí)段。以下分析都是在此時(shí)間段(150~200 ns) 內(nèi)進(jìn)行的。

    對(duì)比5個(gè)模擬體系達(dá)到平衡狀態(tài)時(shí)對(duì)應(yīng)的SASA值發(fā)現(xiàn),相對(duì)于體系V,加入有機(jī)溶劑后染料陰離子的SASA值更高,即染料分子與溶劑分子間的接觸面積大,因此染料分子間的聚集程度小,分散性好。而根據(jù)最后50 ns的SASA值,在4種有機(jī)溶劑中,NMP的加入更有利于C.I.活性藍(lán)176的分散。

    3.2 溶液中染料聚集分析

    溶質(zhì)分子團(tuán)簇的尺寸分布可更直接地表征溶質(zhì)分子間聚集情況,當(dāng)2個(gè)染料陰離子中任意2個(gè)原子之間的距離小于截止值時(shí),就判定2個(gè)分子為聚集態(tài)。截止值設(shè)為0.34 nm,是碳原子范德華半徑的2倍[12],因?yàn)槿玖戏肿又械拇蠖鄶?shù)原子是碳原子。為確定染料在溶液中的聚集情況,在相對(duì)平衡期間計(jì)算了具有不同染料聚集大小的簇的平均數(shù)目,如圖3所示。

    圖3 不同體系中平衡期間具有不同聚集尺寸的染料簇平均數(shù)目Fig.3 Average number of dye clusters with different aggregation sizes during equilibrium period in different systems

    圖3中的染料水溶液體系Ⅴ,染料單分子與二聚體的平均數(shù)目都在13左右,更傾向于以單分子形式存在,表現(xiàn)出活性染料良好的水溶性。而在加入有機(jī)溶劑的Ⅰ~Ⅳ體系中,大多數(shù)染料陰離子以單分子形式存在,主要的二聚體與三聚體的平均數(shù)目明顯降低。其中添加NMP的體系Ⅰ,大約三分之二的染料以單體形式存在溶液中,染料分散性最好。該系統(tǒng)的聚集趨勢(shì)與平衡狀態(tài)下SASA值的趨勢(shì)相同。對(duì)C.I.活性藍(lán)176分散穩(wěn)定性效果遵循以下順序:N-甲基吡咯烷酮>己內(nèi)酰胺>一縮二乙二醇>二甲基亞砜。

    3.3 染料與有機(jī)溶劑相互作用分析

    3.3.1 分子表面靜電勢(shì)

    染料與有機(jī)溶劑分子表面靜電勢(shì)等值面如圖4所示。染料分子的表面靜電勢(shì)數(shù)ESP值整體為負(fù),在紅色的磺酸根氧原子附近具有更多ESP值為負(fù)的位點(diǎn),吸引親電試劑水的能力強(qiáng),這與活性染料良好的水溶性一致。而對(duì)于其他4種有機(jī)溶劑來說,除典型的氧原子附近的ESP值多為負(fù)值外,其他部分的分子表面靜電勢(shì)皆為正值。染料分子與有機(jī)溶劑分子會(huì)因各自對(duì)其外部(范德華表面及更遠(yuǎn)區(qū)域)產(chǎn)生的靜電效應(yīng)造成的二者間的靜電吸引而相互接近,從而發(fā)生疏水相互作用。

    注:ESP的局部最小值和最大值分別表示為青色和橙色的球體,并分別用深藍(lán)色和棕紅色文本標(biāo)記。圖4 染料及有機(jī)溶劑分子的表面靜電勢(shì)Fig.4 Surface electrostatic potential of dye and organic solvent molecular

    對(duì)于分子體系來說,靜電勢(shì)是由體系電荷分布決定的,分子中電荷的分布會(huì)體現(xiàn)在靜電勢(shì)的分布上,因此通過分子表面靜電勢(shì)分析可衡量分子的極性。利用Multiwfn程序[23]還可計(jì)算分子的極性和非極性表面分布。ESP絕對(duì)值小于41.8 kJ/mol的面為非極性面積,反之為極性表面。圖5示出染料分子非極性表面的分布??芍玖系姆菢O性表面主要分布在均三嗪環(huán)及鄰近的苯環(huán)區(qū)域。結(jié)合染料及有機(jī)溶劑分子表面靜電勢(shì)的分析,在染料的非極性區(qū)域二者最有可能發(fā)生相互作用。

    圖5 染料分子非極性部分的靜電勢(shì)分布Fig.5 Electrostatic potential distribution of non-polar part of dye molecule

    3.3.2 空間分布函數(shù)

    利用空間分布函數(shù)可以可視化染料分子周圍有機(jī)溶劑分子及水分子在體系三維空間的分布,如圖6所示。由圖可知,與3.3.1分子表面靜電勢(shì)分析結(jié)果相一致,水分子因受磺酸根氧原子的吸引,主要分布在染料磺酸根附近,而有機(jī)溶劑分子則集中于染料的非極性表面。此外還觀察到,有機(jī)溶劑分子在染料分子的萘環(huán)處也有較多分布,這可能是由于萘環(huán)具有良好的共平面性,而這也是染料聚集的影響因素之一。在模擬體系中有機(jī)溶劑分子對(duì)活性染料分子非極性以及萘環(huán)區(qū)域的“包裹”是有機(jī)溶劑減少染料聚集,提高染料分散穩(wěn)定性的原因。而在同一等值面下, 4種有機(jī)溶劑分子在染料分子周圍分布程度不一,N-甲基吡咯烷酮的包裹更完整,而二甲基亞砜的分布更零散,這可能是對(duì)染料分散效果產(chǎn)生差異的原因,還需要進(jìn)一步的分析。

    注:紅色為水分子,其他不同顏色分別代表4個(gè)有機(jī)溶劑分子。圖6 同等值面下染料周圍水和有機(jī)溶劑分子的空間分布函數(shù)Fig.6 SDF of water and organic solvent molecules around dye under same isosurface

    3.3.3 徑向分布函數(shù)

    為進(jìn)一步闡明染料與有機(jī)溶劑分子的相互作用以及對(duì)染料分散穩(wěn)定性差異的原因,量化了染料周圍的有機(jī)溶劑分子堆積,計(jì)算了染料周圍有機(jī)溶劑分子的徑向分布函數(shù)(RDF)。RDF峰的位置越接近原點(diǎn),峰的強(qiáng)度越高,表明該有機(jī)溶劑分子與染料基團(tuán)相互作用越強(qiáng)。

    在3.3.1和3.3.2節(jié)分析的基礎(chǔ)上,染料分子部分原子周圍有機(jī)溶劑分子的徑向分布函數(shù)如圖7所示。圖7(a)、(b)為染料分子非極性部分與有機(jī)溶劑分子的RDF圖。從圖7(a)可以看出:無論是RDF峰的位置還是峰值,N-甲基吡咯烷酮分子與染料疏水苯環(huán)的相互作用最強(qiáng),且存在雙峰;其次是己內(nèi)酰胺,二甲基亞砜的作用較弱。而對(duì)于圖7(b)中的染料均三嗪環(huán)在其附近,CPL的RDF曲線具有最大的峰值,也即己內(nèi)酰胺分子出現(xiàn)在均三嗪環(huán)周圍的概率更高,但在較短的分離距離(r<0.5 nm)下NMP分子更易接近染料的均三嗪環(huán)。在SDF基礎(chǔ)上分析染料萘環(huán)附近有機(jī)溶劑分子的RDF如圖7(c) 所示,NMP分子在r=0.4 nm處有效堆積,且峰強(qiáng)最高;CPL分子在此處也有概率出現(xiàn),雖然峰的強(qiáng)度不及DEG與DMSO,但后二者RDF中的峰位置(r=0.5 nm)距萘環(huán)表面的距離比己內(nèi)酰胺RDF中的更遠(yuǎn)。

    圖7 不同體系中染料基團(tuán)周圍有機(jī)溶劑分子的徑向分布函數(shù)Fig.7 RDFs of organic solvent molecules around groups of dyes in different systems. (a) Benzene ring;(b) Triazine;(c) Naphthalene nucleus

    綜合上述模擬分析,4種有機(jī)溶劑中,N-甲基吡咯烷酮分子與染料分子相互作用最強(qiáng),對(duì)染料聚集有效區(qū)域的“包裹”最緊密,從而阻礙了染料分子的聚集,提高了液體活性染料的分散穩(wěn)定性。

    在染料分散穩(wěn)定性分析的基礎(chǔ)上,對(duì)于一氯均三嗪型活性藍(lán)176染料來說(見圖8),染料因分子中不穩(wěn)定的氯原子被羥基所取代進(jìn)而發(fā)生水解。有機(jī)溶劑分子在氯原子周圍的分布差異會(huì)影響染料的水解。圖9示出染料氯原子周圍有機(jī)溶劑分子的徑向分布函數(shù)。由圖可知,4種有機(jī)溶劑分子RDF峰的位置都在0.4 nm左右,但峰強(qiáng)不同,其中CPL分子在染料氯原子周圍有較多分布,且存在雙峰,對(duì)氯原子的“包裹”阻礙了羥基的取代,從而緩解了染料的水解。DEG與DMSO分子在氯原子周圍分布較少,緩解染料水解效果相對(duì)較差,對(duì)抑制C.I.活性藍(lán)176 水解效果為:己內(nèi)酰胺>N-甲基吡咯烷酮>一縮二乙二醇>二甲基亞砜。

    注:D代表母體染料。圖8 染料的水解示意圖Fig.8 Schematic diagram of dye hydrolysis

    圖9 不同體系中染料氯原子周圍有機(jī)溶劑分子的徑向分布函數(shù)Fig.9 RDFs of organic solvent molecules around chlorine atoms of dyes in different systems.

    4 實(shí)驗(yàn)結(jié)果驗(yàn)證

    4.1 染料的分散穩(wěn)定性

    以不同有機(jī)溶劑按照表1比例配制體系Ⅰ~Ⅴ,比較靜置72 h后生成的沉淀質(zhì)量、染料的平均粒徑以及溶液紫外-可見吸收光譜來評(píng)價(jià)染料的分散穩(wěn)定性。

    表2示出不同體系生成的沉淀質(zhì)量及染料平均粒徑??梢钥闯觯袡C(jī)溶劑的添加減少了染料的沉淀,染料的粒徑變小,抑制了染料的聚集。向染料溶液中添加有機(jī)溶劑會(huì)導(dǎo)致610 nm左右的染料吸收峰高明顯增加,4種有機(jī)溶劑均可增加染料溶液的吸收,但添加NMP的染料吸收峰比其他有機(jī)溶劑高,也即NMP表現(xiàn)出更強(qiáng)大的增強(qiáng)染料溶液吸收的能力。這意味著NMP的加入可幫助染料溶液釋放更多的單染料分子。實(shí)驗(yàn)結(jié)果表明,N-甲基吡咯烷酮對(duì)C.I.活性藍(lán)176的分散穩(wěn)定性效果最好,其次是己內(nèi)酰胺,與模擬分析結(jié)果具有很好的一致性,如圖10所示。值得注意的是,有機(jī)溶劑N-甲基吡咯烷酮被國際環(huán)保紡織協(xié)會(huì)列入“溶劑殘留物”目錄, 在2020年新標(biāo)準(zhǔn)中限量值為0.05%(質(zhì)量分?jǐn)?shù))。在實(shí)際使用中應(yīng)盡量減少NMP的使用。

    表2 不同體系生成的沉淀質(zhì)量及染料平均粒徑Tab.2 Precipitation quality and average particle size of dyes produced in different systems

    4.2 染料的水解穩(wěn)定性

    C.I.活性藍(lán)176水解的高效液相譜圖如圖11所示。

    圖10 不同體系的染料溶液的紫外-可見吸收光譜變化Fig.10 UV-Vis absorption spectra change of dye solutions in different systems

    圖11 不同體系中C.I.活性藍(lán)176水解的高效液相譜Fig.11 HPLC spectra of C.I. Reactive Blue 176 hydrolyzed in different systems. (a) Peak of dyes; (b) Peak of hydrolyzed dyes

    從圖11(a)可以觀察到,染料在水解過程中有2種物質(zhì)存在,對(duì)于水解時(shí)間0 h的染料溶液,只有保留時(shí)間為1.7 min左右的峰,而在水解12 h后在保留時(shí)間2 min附近又出現(xiàn)峰。由此可以確定,保留時(shí)間1.7 min左右的是RB176,保留時(shí)間2 min左右的是水解染料。結(jié)合圖11(b)可知,水解12 h后,相對(duì)于體系Ⅴ,添加了有機(jī)溶劑的都在一定程度上緩解了染料的水解。其中加入CPL的體系Ⅱ在水解12 h后,染料峰的降低及生成的水解染料都較少,對(duì)染料的保護(hù)更好。其次是添加N-甲基吡咯烷酮的體系Ⅰ,染料水解也較少。符合模擬的分析結(jié)果。

    5 結(jié) 論

    本文基于分子動(dòng)力學(xué)模擬方法,分別研究了液體活性染料體系中4種常用的有機(jī)溶劑對(duì)C.I.活性藍(lán)176分散及水解穩(wěn)定性的影響,探索活性染料在不同溶劑條件下分散與水解效果差異的原因。主要結(jié)論如下:

    1)模擬結(jié)果表明,有機(jī)溶劑分子對(duì)染料疏水及萘環(huán)區(qū)域的“包裹”阻礙了染料分子的聚集,而與染料相互作用的強(qiáng)弱決定了對(duì)染料的分散穩(wěn)定性。4種有機(jī)溶劑中,N-甲基吡咯烷酮分子與染料疏水區(qū)域相互作用最強(qiáng),周圍分布最多,所以對(duì)染料的分散穩(wěn)定性效果最好。與之類似,己內(nèi)酰胺在一氯均三嗪染料分子中不穩(wěn)定的氯原子周圍分布最多,相互作用最強(qiáng),減小羥基的作用力,從而抑制了染料的水解。

    2)實(shí)驗(yàn)結(jié)果驗(yàn)證了模擬分析,有機(jī)溶劑的添加都在一定程度上提高了染料的分散穩(wěn)定性,緩解了染料的水解。考慮到N-甲基吡咯烷酮使用過程中在織物上的殘留可能對(duì)人體存在的潛在危害,并且對(duì)于染料的分散穩(wěn)定性,添加N-甲基吡咯烷酮與己內(nèi)酰胺的效果差異不大,而己內(nèi)酰胺的加入對(duì)抑制活性染料的水解效果更好,所以綜合對(duì)比,4種有機(jī)溶劑中,優(yōu)選己內(nèi)酰胺作為提高C.I.活性藍(lán)176儲(chǔ)存穩(wěn)定性的有機(jī)溶劑。

    猜你喜歡
    靜電勢(shì)吡咯烷酮活性染料
    N-甲基吡咯烷酮降解菌株的篩選鑒定及應(yīng)用
    表面活性劑十八烷基磺酸鈉的Multiwfn研究*
    廣州化工(2022年19期)2022-11-09 11:30:46
    一種水處理劑:氨基三亞甲基膦酸的Multiwfn研究*
    廣州化工(2022年18期)2022-10-22 10:27:00
    聚乙烯吡咯烷酮分子三級(jí)中紅外光譜研究
    菠蘿纖維活性染料染色工藝研究
    毛紡科技(2021年8期)2021-10-14 06:51:20
    水分子在高嶺石(001)面吸附的密度泛函計(jì)算
    歐盟重新評(píng)估聚乙烯吡咯烷酮(E1201)和聚乙烯聚吡咯烷酮(E1202)作為食品添加劑的安全性
    纖維素纖維負(fù)載鈷酞菁對(duì)活性染料X-3B的降解
    聚乙烯吡咯烷酮對(duì)Pt/Ru雙金屬納米簇的結(jié)構(gòu)和催化性能影響
    L型活性染料對(duì)蠶絲的低溫染色研究
    絲綢(2014年7期)2014-02-28 14:55:24
    欧美xxxx性猛交bbbb| 国产精品影院久久| 久久精品国产亚洲av香蕉五月| 成人特级av手机在线观看| 男人和女人高潮做爰伦理| 老司机午夜福利在线观看视频| 午夜日韩欧美国产| 国内少妇人妻偷人精品xxx网站| 在线免费观看的www视频| 欧美黑人巨大hd| 亚洲第一欧美日韩一区二区三区| 亚洲欧美日韩东京热| 日本五十路高清| 欧美黑人欧美精品刺激| 麻豆国产av国片精品| 天堂影院成人在线观看| 午夜两性在线视频| 亚洲av美国av| 日韩亚洲欧美综合| 欧美潮喷喷水| 此物有八面人人有两片| 精品久久国产蜜桃| 在现免费观看毛片| 可以在线观看毛片的网站| 动漫黄色视频在线观看| 一二三四社区在线视频社区8| 日韩高清综合在线| 日本一本二区三区精品| 精品一区二区三区人妻视频| 亚洲国产欧洲综合997久久,| 久久久久九九精品影院| 午夜精品久久久久久毛片777| 直男gayav资源| 亚洲男人的天堂狠狠| 一级av片app| 久久香蕉精品热| 精品一区二区三区人妻视频| 少妇丰满av| 极品教师在线免费播放| 欧美乱妇无乱码| 久久伊人香网站| 搞女人的毛片| 国产精品野战在线观看| 国产av麻豆久久久久久久| 亚洲av二区三区四区| 午夜福利高清视频| 亚洲综合色惰| 亚洲内射少妇av| 亚洲av电影在线进入| 日本 欧美在线| 1024手机看黄色片| 国产人妻一区二区三区在| 国产日本99.免费观看| 一卡2卡三卡四卡精品乱码亚洲| 九色国产91popny在线| a级一级毛片免费在线观看| 韩国av一区二区三区四区| h日本视频在线播放| 国产一级毛片七仙女欲春2| 97碰自拍视频| 18禁黄网站禁片免费观看直播| 美女高潮喷水抽搐中文字幕| av在线天堂中文字幕| 亚洲成av人片免费观看| 一卡2卡三卡四卡精品乱码亚洲| 国内精品一区二区在线观看| 高清日韩中文字幕在线| 成人av在线播放网站| 日日夜夜操网爽| 每晚都被弄得嗷嗷叫到高潮| 午夜福利免费观看在线| 少妇人妻精品综合一区二区 | 亚洲一区二区三区不卡视频| 91av网一区二区| 狠狠狠狠99中文字幕| 如何舔出高潮| 夜夜看夜夜爽夜夜摸| 在线免费观看的www视频| 午夜精品在线福利| 99久久99久久久精品蜜桃| 波多野结衣高清无吗| 久久久精品欧美日韩精品| 狂野欧美白嫩少妇大欣赏| 丁香六月欧美| 亚洲欧美日韩无卡精品| 久久天躁狠狠躁夜夜2o2o| 黄色丝袜av网址大全| 欧美区成人在线视频| 国产精品一及| 日韩 亚洲 欧美在线| 成人精品一区二区免费| 最好的美女福利视频网| 99久久无色码亚洲精品果冻| 99久久成人亚洲精品观看| 又紧又爽又黄一区二区| 国产精品一区二区三区四区久久| 少妇人妻精品综合一区二区 | 亚洲精品粉嫩美女一区| 亚洲精品456在线播放app | 男人舔女人下体高潮全视频| 亚洲五月婷婷丁香| 午夜免费男女啪啪视频观看 | 日韩人妻高清精品专区| 男女之事视频高清在线观看| 国产大屁股一区二区在线视频| 久久精品91蜜桃| 亚洲熟妇中文字幕五十中出| 国产日本99.免费观看| 内射极品少妇av片p| 美女免费视频网站| 最近最新免费中文字幕在线| 免费一级毛片在线播放高清视频| 2021天堂中文幕一二区在线观| 免费大片18禁| 国产色爽女视频免费观看| 悠悠久久av| 欧美日韩中文字幕国产精品一区二区三区| 欧美成人免费av一区二区三区| 伦理电影大哥的女人| 色综合亚洲欧美另类图片| 禁无遮挡网站| 啦啦啦韩国在线观看视频| 日韩精品青青久久久久久| 欧美一区二区亚洲| 成人亚洲精品av一区二区| 欧美日韩综合久久久久久 | 国产精品自产拍在线观看55亚洲| 亚洲不卡免费看| 成年女人看的毛片在线观看| av在线天堂中文字幕| 精品久久久久久成人av| 国产爱豆传媒在线观看| 亚洲国产精品成人综合色| 在线播放国产精品三级| 此物有八面人人有两片| 中国美女看黄片| 亚洲欧美日韩高清在线视频| 亚洲精品在线观看二区| 日韩高清综合在线| 色5月婷婷丁香| 久久精品夜夜夜夜夜久久蜜豆| 18禁裸乳无遮挡免费网站照片| 夜夜躁狠狠躁天天躁| 中文字幕av成人在线电影| 欧美日韩乱码在线| 久久人妻av系列| 日韩人妻高清精品专区| 18美女黄网站色大片免费观看| 亚洲av成人精品一区久久| 成人国产一区最新在线观看| 毛片一级片免费看久久久久 | 色5月婷婷丁香| 久久亚洲精品不卡| 欧美xxxx性猛交bbbb| 天天一区二区日本电影三级| 九色成人免费人妻av| 亚洲一区高清亚洲精品| 国内精品美女久久久久久| 90打野战视频偷拍视频| 国产大屁股一区二区在线视频| ponron亚洲| 国产黄色小视频在线观看| 国产精品永久免费网站| 亚洲乱码一区二区免费版| 精品久久久久久久末码| 国产激情偷乱视频一区二区| 国产精品av视频在线免费观看| 午夜精品在线福利| 伊人久久精品亚洲午夜| 色哟哟·www| 又爽又黄无遮挡网站| 亚洲在线自拍视频| 国产乱人视频| 日韩高清综合在线| 小蜜桃在线观看免费完整版高清| 婷婷六月久久综合丁香| 很黄的视频免费| 网址你懂的国产日韩在线| 欧美精品国产亚洲| 亚洲av电影不卡..在线观看| 日本成人三级电影网站| 我要看日韩黄色一级片| 天堂网av新在线| 精品人妻1区二区| 色精品久久人妻99蜜桃| 欧美日本视频| 国产成年人精品一区二区| aaaaa片日本免费| 国产乱人伦免费视频| 色吧在线观看| 99久国产av精品| 久久香蕉精品热| 嫩草影院精品99| 国产久久久一区二区三区| 两个人的视频大全免费| 日韩欧美免费精品| 成人一区二区视频在线观看| av欧美777| 99在线视频只有这里精品首页| 一本综合久久免费| 国产国拍精品亚洲av在线观看| 如何舔出高潮| 中文字幕av成人在线电影| 久久久久久久精品吃奶| 午夜视频国产福利| 久久婷婷人人爽人人干人人爱| 日韩有码中文字幕| 亚洲人成网站在线播放欧美日韩| 超碰av人人做人人爽久久| 国产色爽女视频免费观看| 色在线成人网| 简卡轻食公司| 国产精品久久久久久久电影| 国产精品久久久久久人妻精品电影| 怎么达到女性高潮| eeuss影院久久| 简卡轻食公司| 久久精品91蜜桃| 一夜夜www| 怎么达到女性高潮| 高清日韩中文字幕在线| 亚洲avbb在线观看| 人人妻人人澡欧美一区二区| 嫁个100分男人电影在线观看| 国产精品不卡视频一区二区 | 99久久99久久久精品蜜桃| 丰满的人妻完整版| 99久国产av精品| 国产精品综合久久久久久久免费| 国产成人影院久久av| 国产欧美日韩精品一区二区| 久久婷婷人人爽人人干人人爱| 中文字幕精品亚洲无线码一区| 国产精品1区2区在线观看.| 国产一区二区激情短视频| 99热这里只有精品一区| 亚洲国产精品999在线| 久久久色成人| 在线a可以看的网站| 露出奶头的视频| 精品久久久久久久久久久久久| 欧美又色又爽又黄视频| 精品一区二区三区视频在线| 日本成人三级电影网站| 亚洲国产精品成人综合色| 中文字幕av在线有码专区| 我的女老师完整版在线观看| 中文字幕免费在线视频6| 日韩人妻高清精品专区| 欧美色视频一区免费| 白带黄色成豆腐渣| 日本免费a在线| 中文字幕久久专区| 搡女人真爽免费视频火全软件 | 乱人视频在线观看| 国内精品久久久久精免费| 日韩中字成人| 首页视频小说图片口味搜索| 男人舔奶头视频| 免费观看的影片在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产三级黄色录像| 亚洲人成网站高清观看| avwww免费| 99热6这里只有精品| 精品国产亚洲在线| 亚洲五月天丁香| 国产精品99久久久久久久久| 长腿黑丝高跟| 又黄又爽又刺激的免费视频.| 搡老岳熟女国产| 少妇人妻精品综合一区二区 | 97热精品久久久久久| 日日摸夜夜添夜夜添av毛片 | 精品人妻熟女av久视频| 欧美三级亚洲精品| 国模一区二区三区四区视频| 欧美不卡视频在线免费观看| 国产黄色小视频在线观看| 看片在线看免费视频| 国产91精品成人一区二区三区| 又爽又黄无遮挡网站| 中文字幕av成人在线电影| 午夜福利视频1000在线观看| 一级黄色大片毛片| 久久天躁狠狠躁夜夜2o2o| 老司机深夜福利视频在线观看| 精品久久久久久久久av| 脱女人内裤的视频| 亚洲18禁久久av| 色尼玛亚洲综合影院| 国产精品久久久久久亚洲av鲁大| 欧美成狂野欧美在线观看| 亚洲经典国产精华液单 | 亚洲电影在线观看av| 免费人成视频x8x8入口观看| 两性午夜刺激爽爽歪歪视频在线观看| 少妇裸体淫交视频免费看高清| 日日夜夜操网爽| 亚洲av免费在线观看| 亚洲国产精品久久男人天堂| 精品国产三级普通话版| 免费黄网站久久成人精品 | 黄色丝袜av网址大全| 18禁黄网站禁片免费观看直播| x7x7x7水蜜桃| 亚洲精品在线观看二区| 日韩 亚洲 欧美在线| 亚洲一区二区三区色噜噜| 久久久久国内视频| 熟女电影av网| 欧美成人一区二区免费高清观看| 国产乱人伦免费视频| 亚洲七黄色美女视频| 91在线观看av| or卡值多少钱| 91九色精品人成在线观看| 国产一区二区激情短视频| 天堂网av新在线| 免费电影在线观看免费观看| 天堂影院成人在线观看| 宅男免费午夜| 尤物成人国产欧美一区二区三区| 亚洲激情在线av| 欧美极品一区二区三区四区| 色噜噜av男人的天堂激情| 老司机午夜福利在线观看视频| 国产男靠女视频免费网站| 国产真实乱freesex| 不卡一级毛片| 麻豆国产av国片精品| 成年女人看的毛片在线观看| 男人和女人高潮做爰伦理| 熟女电影av网| 毛片女人毛片| 午夜福利在线观看免费完整高清在 | 亚洲精品粉嫩美女一区| 观看免费一级毛片| 精品午夜福利在线看| 尤物成人国产欧美一区二区三区| 国产私拍福利视频在线观看| 中文在线观看免费www的网站| 国产不卡一卡二| 一级av片app| 国产午夜精品久久久久久一区二区三区 | 深夜a级毛片| 国产亚洲欧美98| 国产av在哪里看| 亚洲人成伊人成综合网2020| 国产三级在线视频| 欧美丝袜亚洲另类 | 亚洲av电影在线进入| 波多野结衣巨乳人妻| 精品久久久久久久久av| 亚洲中文日韩欧美视频| 俄罗斯特黄特色一大片| 一个人看视频在线观看www免费| 超碰av人人做人人爽久久| 亚洲无线观看免费| 嫩草影院新地址| 欧美+日韩+精品| 亚洲性夜色夜夜综合| 看片在线看免费视频| 女生性感内裤真人,穿戴方法视频| 国产精品电影一区二区三区| 亚洲无线在线观看| 欧美日韩国产亚洲二区| 天美传媒精品一区二区| 亚洲av五月六月丁香网| 午夜福利在线观看吧| 网址你懂的国产日韩在线| 毛片一级片免费看久久久久 | 日韩人妻高清精品专区| 很黄的视频免费| av在线观看视频网站免费| 两人在一起打扑克的视频| 久久欧美精品欧美久久欧美| 最新在线观看一区二区三区| 亚洲国产精品999在线| 中出人妻视频一区二区| 色噜噜av男人的天堂激情| 99热这里只有精品一区| 高清日韩中文字幕在线| 啦啦啦观看免费观看视频高清| 亚洲av成人av| av专区在线播放| 有码 亚洲区| 午夜视频国产福利| 一级毛片久久久久久久久女| 久久久久国产精品人妻aⅴ院| 欧美在线一区亚洲| 熟女电影av网| 国产真实乱freesex| 大型黄色视频在线免费观看| 麻豆av噜噜一区二区三区| 深夜a级毛片| 亚洲av不卡在线观看| 搡女人真爽免费视频火全软件 | 亚洲欧美清纯卡通| 精品一区二区三区人妻视频| 中文字幕精品亚洲无线码一区| 日本五十路高清| 国产成人欧美在线观看| av欧美777| 又爽又黄a免费视频| 91麻豆精品激情在线观看国产| 国产精品一区二区免费欧美| 亚洲内射少妇av| 少妇的逼水好多| 中文字幕熟女人妻在线| 久久精品国产亚洲av涩爱 | 精品一区二区三区人妻视频| 国产成人av教育| 亚洲精品在线观看二区| 亚洲av不卡在线观看| 国产伦精品一区二区三区视频9| 午夜久久久久精精品| 美女黄网站色视频| 啦啦啦韩国在线观看视频| 欧美一区二区国产精品久久精品| 国产精品99久久久久久久久| 性插视频无遮挡在线免费观看| av福利片在线观看| 亚洲最大成人手机在线| 久久欧美精品欧美久久欧美| 亚洲欧美精品综合久久99| 在线观看一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 男女床上黄色一级片免费看| 日韩欧美在线乱码| 欧美激情在线99| 小说图片视频综合网站| 亚洲国产色片| av中文乱码字幕在线| 性欧美人与动物交配| 少妇丰满av| 老女人水多毛片| 少妇的逼好多水| 在线国产一区二区在线| 久久6这里有精品| 波野结衣二区三区在线| 一本一本综合久久| 亚洲精品一区av在线观看| 亚洲av不卡在线观看| 亚洲激情在线av| 啦啦啦韩国在线观看视频| 亚洲av电影在线进入| 亚洲最大成人中文| 男人的好看免费观看在线视频| 人妻夜夜爽99麻豆av| 久久久久性生活片| 97人妻精品一区二区三区麻豆| 亚洲人成网站在线播| 97超视频在线观看视频| 欧美高清成人免费视频www| 69人妻影院| 一级黄色大片毛片| 日韩 亚洲 欧美在线| 日韩欧美精品v在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲18禁久久av| 在线播放无遮挡| 丝袜美腿在线中文| 午夜两性在线视频| 很黄的视频免费| 在线观看一区二区三区| 亚洲欧美日韩高清在线视频| 偷拍熟女少妇极品色| 国产黄片美女视频| 久9热在线精品视频| 免费在线观看亚洲国产| 成人三级黄色视频| 国产色爽女视频免费观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲最大成人中文| 欧美成狂野欧美在线观看| 亚洲精品在线观看二区| 直男gayav资源| 99在线人妻在线中文字幕| 色5月婷婷丁香| 欧美xxxx黑人xx丫x性爽| av在线蜜桃| 99久久精品国产亚洲精品| 久久这里只有精品中国| 国产欧美日韩精品一区二区| 岛国在线免费视频观看| 身体一侧抽搐| 日本黄色片子视频| 国产精品一及| netflix在线观看网站| 老鸭窝网址在线观看| 国产一级毛片七仙女欲春2| 国产视频内射| 国产精品久久久久久精品电影| 丁香六月欧美| 日韩av在线大香蕉| 国产精品一区二区性色av| 欧美黄色淫秽网站| 精品乱码久久久久久99久播| 亚洲无线观看免费| 桃红色精品国产亚洲av| 2021天堂中文幕一二区在线观| 美女高潮的动态| 亚洲精品日韩av片在线观看| 好男人在线观看高清免费视频| 一级av片app| 熟女电影av网| 亚洲黑人精品在线| 人人妻人人看人人澡| 国产高潮美女av| 亚洲欧美日韩高清专用| av天堂中文字幕网| 最近中文字幕高清免费大全6 | 亚洲黑人精品在线| 欧美不卡视频在线免费观看| 丝袜美腿在线中文| 久久国产精品影院| 国产毛片a区久久久久| 午夜亚洲福利在线播放| 欧美绝顶高潮抽搐喷水| a级毛片a级免费在线| 中文字幕av在线有码专区| 90打野战视频偷拍视频| 在线国产一区二区在线| 国产精品爽爽va在线观看网站| 欧美日韩国产亚洲二区| 我要搜黄色片| 99久久九九国产精品国产免费| 国产精品久久视频播放| 亚洲国产高清在线一区二区三| 午夜老司机福利剧场| 又紧又爽又黄一区二区| 最近最新免费中文字幕在线| 精品国产亚洲在线| 午夜激情福利司机影院| 精品人妻1区二区| 中文字幕精品亚洲无线码一区| 成人三级黄色视频| 久久精品国产亚洲av天美| 亚洲国产欧洲综合997久久,| 看片在线看免费视频| 国产亚洲欧美在线一区二区| 国产伦精品一区二区三区四那| 久久人人爽人人爽人人片va | 成人国产一区最新在线观看| 首页视频小说图片口味搜索| 99久久精品热视频| 免费人成在线观看视频色| 我要看日韩黄色一级片| 日韩精品中文字幕看吧| 午夜免费激情av| 亚洲专区中文字幕在线| 女人被狂操c到高潮| 国产伦在线观看视频一区| 国产亚洲精品久久久久久毛片| 欧美xxxx黑人xx丫x性爽| 日本黄色视频三级网站网址| 欧美日韩福利视频一区二区| 波多野结衣高清无吗| 成年免费大片在线观看| 热99在线观看视频| 天堂影院成人在线观看| 亚洲精品色激情综合| 看片在线看免费视频| 女人十人毛片免费观看3o分钟| 国产极品精品免费视频能看的| 97超级碰碰碰精品色视频在线观看| 18禁黄网站禁片免费观看直播| 国产色爽女视频免费观看| 久久久精品欧美日韩精品| 国产乱人伦免费视频| 精品欧美国产一区二区三| 97碰自拍视频| 免费在线观看亚洲国产| 久久精品久久久久久噜噜老黄 | 91在线精品国自产拍蜜月| 国产老妇女一区| 亚洲国产色片| av天堂在线播放| 精品日产1卡2卡| 国产精品av视频在线免费观看| 国内精品久久久久精免费| 亚洲av电影在线进入| 国产午夜精品论理片| 99国产精品一区二区蜜桃av| 999久久久精品免费观看国产| 欧美色欧美亚洲另类二区| 日本 av在线| 男女那种视频在线观看| 国产av不卡久久| 脱女人内裤的视频| 亚洲无线在线观看| 亚洲,欧美,日韩| 中国美女看黄片| 99久久久亚洲精品蜜臀av| 婷婷丁香在线五月| 男女床上黄色一级片免费看| 人妻丰满熟妇av一区二区三区| 黄片小视频在线播放| 久久久精品大字幕| 亚洲狠狠婷婷综合久久图片| 天堂av国产一区二区熟女人妻| 88av欧美| 欧美一级a爱片免费观看看| 成年女人永久免费观看视频| 特级一级黄色大片| 亚洲精品色激情综合| 性色avwww在线观看| 国内精品美女久久久久久| 色视频www国产| 九九久久精品国产亚洲av麻豆| 亚洲av成人精品一区久久| 听说在线观看完整版免费高清| 老熟妇乱子伦视频在线观看| 男人的好看免费观看在线视频|