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

    溶劑對丁腈橡膠溶解、尺寸、結(jié)構(gòu)和催化加氫的影響研究

    2021-01-30 08:10:28葛冰青陰義軒王亞溪張宏偉袁珮
    化工學(xué)報 2021年1期
    關(guān)鍵詞:雙鍵鏈段溶解度

    葛冰青,陰義軒,王亞溪,張宏偉,袁珮

    (福州大學(xué)石油化工學(xué)院,福建福州350108)

    引 言

    丁腈橡膠(NBR)是一種重要的合成橡膠,但由于分子鏈中碳碳雙鍵的存在,其耐溫性和耐老化性較差,應(yīng)用也受到一定限制。NBR 經(jīng)選擇性加氫得到的產(chǎn)物氫化丁腈橡膠(HNBR),不僅繼承了NBR的耐油、耐磨等性能,而且具有更優(yōu)異的耐熱、耐氧化、耐臭氧、耐化學(xué)品腐蝕等性能,同時還具有高強度、高撕裂性能、較高的抗壓縮永久變形性能,是綜合性能極為出色的橡膠之一,應(yīng)用領(lǐng)域得到極大的拓展[1-4]。

    NBR均相與非均相催化溶液加氫是制備HNBR的兩種主要方法[5-8]。其中,均相溶液加氫法的缺點是催化劑不易與產(chǎn)物分離,這不僅會消耗大量的貴金屬,造成制備成本較高,且產(chǎn)物中殘存的催化劑會影響HNBR 的使用性能,給實際生產(chǎn)和應(yīng)用帶來很大的隱患。相比之下,用非均相溶液加氫法制備HNBR 時,加氫催化劑容易分離和回收,不僅可以降低生產(chǎn)成本,而且還能保證產(chǎn)品質(zhì)量,展現(xiàn)出更好的工業(yè)應(yīng)用前景[9-10]。非均相溶液加氫中需要用到有機溶劑,目前的研究報道NBR 溶解和加氫時多用丙酮或四氯化碳作為溶劑,對于其他溶劑體系較少涉及,目前研究者多集中于高效催化劑的開發(fā)[11-15]。溶劑作為反應(yīng)媒介發(fā)揮著至關(guān)重要的作用[16-20],不僅能溶解反應(yīng)物,影響NBR 的構(gòu)型和尺寸,還會與反應(yīng)物、催化劑相互作用,影響反應(yīng)物的擴散以及在活性位點的吸附與活化,進而引起反應(yīng)速率、加氫活性和選擇性的變化。因此,研究溶劑對NBR 溶解、尺寸和結(jié)構(gòu)的影響并確定適合NBR 加氫反應(yīng)的溶劑體系至關(guān)重要。

    NBR 在溶劑中的溶解不同于其他簡單小分子,小分子化合物通常采用溶解度來描述其在溶劑中的溶解情況(如食鹽在水中的溶解),在溶解過程中溶劑分子與溶質(zhì)分子或者離子之間相互吸引,導(dǎo)致溶質(zhì)分子與溶質(zhì)分子互相分離,繼而擴散,最后溶質(zhì)均勻分布到溶劑中形成均相溶液。高分子化合物的溶解過程較為復(fù)雜,分為溶脹階段和溶解階段兩部分,其溶解表現(xiàn)為由表及內(nèi)、從外到內(nèi)逐漸溶劑化,最終完全溶解。溶劑對高分子化合物的溶解能力取決于溶劑分子與高分子之間的親和作用,且這種相互作用還會影響聚合物在溶劑中的尺寸和構(gòu)象。聚合物的溶解除了要遵循“溶解度參數(shù)相近”原則,還要考慮“溶劑化作用”,換句話說,聚合物和溶劑分子要滿足雙方分別為電子受體(親電性)或電子給體(親核性),才會產(chǎn)生溶劑化作用。尤為重要的是,NBR 在溶液中的尺寸和構(gòu)型都依賴于NBR 分子與溶劑間的相互作用,從而會進一步影響NBR 加氫速率和加氫度[21-22]。而溶劑對NBR 溶解、尺寸和結(jié)構(gòu)的影響問題并沒有深入的研究,因此,本文采用實驗表征結(jié)合密度泛函理論的方法系統(tǒng)地考察了不同溶劑對NBR 溶解、尺寸及結(jié)構(gòu)的影響,試圖探索溶劑性質(zhì)與NBR 分子尺寸和構(gòu)型的對應(yīng)關(guān)系,并研究了不同溶劑對NBR 催化加氫性能的影響,為構(gòu)建適合NBR 非均相加氫反應(yīng)的溶劑體系提供理論基礎(chǔ)。

    1 實驗材料和方法

    1.1 實驗材料

    根據(jù)溶劑的性質(zhì)不同,將其分為三類:強親電子溶劑、弱親電子溶劑及給電子溶劑。所選溶劑如下:

    (1) 強親電子溶劑:乙醇(EtOH, δ=26)、乙二醇(EG,δ=33.3)、異丙醇(IPA,δ=23)、水(water,δ=47.8)和乙腈(MeCN,δ=24.3);

    (2)弱親電子溶劑:甲苯(TOL,δ=18.2)、鄰二甲苯(OX, δ=18.4)、間二甲苯(MX, δ=18.0)、對二甲苯(PX,δ=17.9)、異丙苯(TIPB, δ=15.5)、氯苯(MCB, δ=19.8)、二氯甲烷(DCM, δ=19.8)、氯仿(TCM, δ=19.8)、1,2-二氯乙烷(EDC,δ=19.6);

    (3)給電子溶劑:丙酮(ACE,δ=20.2)、丁酮(MEK,δ=19.0)、4-甲基-2-戊酮(MIBK,δ=17)、環(huán)己酮(CYC,δ=20.2)、N-甲基吡咯烷酮(NMP, δ=23.0)、乙酸甲酯(MAC,δ=19.6)、乙酸乙酯(EA, δ=17)、乙酸異丁酯(IBAC, δ=16.4)、N,N-二甲基甲酰胺(DMF, δ=24.8); δ為溶解度參數(shù),(J/cm3)0.5。

    上述溶劑由阿拉丁提供。醋酸鈀[Pd(OAc)2,鈀含量為46.0%~48.0%(質(zhì)量分?jǐn)?shù))]、N-氨乙基-γ-氨丙基三甲氧基硅烷由阿拉丁提供。丁腈橡膠(NBR,重均分子量為3.6×105,丙烯腈質(zhì)量分?jǐn)?shù)為36.5%)由蘭州石化提供。

    1.2 實驗方案

    1.2.1 溶劑對NBR 的溶解 量取40 ml 上述一種強親電子溶劑、弱親電子溶劑或給電子溶劑置于廣口瓶中,將0.4 g 的NBR 放置到溶劑中,對比NBR 溶解前后的形態(tài),對NBR 在不同溶劑中的溶解行為進行研究。

    1.2.2 不同溶劑中NBR 的表征 不同溶劑中NBR的尺寸測定:配制濃度為1 g/L 的NBR 膠液,采用馬爾文公司的Zetasizer Nano ZS90 型納米粒度儀對NBR 的粒徑進行測定,測試溫度為25℃,光源為633 nm 的He-Ne 激光器,90°的光學(xué)散射原件。樣品制備方法:配制待測液,將待測液放置在1.5 ml的錐形離心管中,在轉(zhuǎn)速為10000 r/min 的條件下離心30 min,離心結(jié)束后,用吸管吸取離心后樣品的上清液,垂直加入到潔凈的比色皿中,關(guān)閉樣品倉門,在軟件中點擊開始等待測量結(jié)束。

    不同溶劑中NBR 膠液黏度測定:不同溶劑中NBR 膠液的黏度η 采用烏氏黏度計進行測定,調(diào)節(jié)恒溫槽的溫度范圍于(30±1)℃。測定5次(相差0.2~0.3 s),取平均值。聚合物稀溶液的黏度反映了溶液在流動時內(nèi)摩擦力的大小。在此動態(tài)過程中,通常將內(nèi)摩擦分為三類:純?nèi)軇┓肿娱g的內(nèi)摩擦,記作η0;聚合物分子與溶劑分子間的內(nèi)摩擦;聚合物分子間的內(nèi)摩擦。三種內(nèi)摩擦綜合的結(jié)果稱為聚合物溶液的黏度,記作η。事實表明,等溫條件下:η>η0。采用增比黏度(ηsp)進一步比較η和η0,ηsp反映了聚合物與聚合物之間、純?nèi)軇┡c聚合物之間的內(nèi)摩擦效應(yīng),由式(1)得到:

    通常情況下,溶液濃度C 增加,聚合物溶液ηsp隨之增大。作出如下定義便于比較:比濃黏度(單位濃度的增比黏度ηsp/C)、比濃對數(shù)黏度(lnηr/C)。

    當(dāng)溶液無限稀時又滿足:

    溶液被無限稀釋,聚合物分子之間距離較大,相互作用力可以被忽略。因此,通過將ηsp/C 對C 作圖,外推到C→0 的截距值,同時將ηr/C 對C 作圖,根據(jù)外推法,得到C→0的截距值,兩根線合于該點,該點縱坐標(biāo)的值即為所求特性黏度[η]。

    1.2.3 密度泛函理論(DFT)計算 密度泛函理論的計算在Gaussian 09 軟件上進行,結(jié)構(gòu)建模和最終結(jié)構(gòu)的可視化采用Gaussian View,選用包含了彌散函數(shù)和極化函數(shù)[23]的B3LYP/6-31G 基組水平對NBR分子的初始構(gòu)象進行優(yōu)化和計算[24]。在此基組水平上,采用自洽反應(yīng)理論中的SMD 模型[25]計算NBR 分子片段在不同溶劑氛圍中的偶極矩。結(jié)合SMD 溶劑化模型對NBR 分子在氣相和不同溶劑中的C C雙鍵、前線軌道和溶劑化能等做了統(tǒng)計分析。

    NBR分子建模:本實驗所用NBR的分子量為36萬,該聚合物是由千百個C、H 和N 原子以共價鍵形式連接成的高分子,其分布范圍較寬,因此對整個大分子進行模擬計算存在一定的困難。因而,根據(jù)對NBR 的H1NMR 譜圖(圖1)分析可知,NBR 分子內(nèi)鏈重復(fù)單元包括:1,4-丁二烯(x)、丙烯腈(y)和1,2-丁二烯(z),其中,H 化學(xué)位移δ =5.50 左右的峰對應(yīng)NBR 中1,4-C C—雙鍵上的H;在δ=4.9~5.10 范圍的峰對應(yīng)NBR 中1,2-C C—雙鍵上的H;δ =2.65左右的峰對應(yīng)NBR 中—CH(CN)—結(jié)構(gòu)上的H(表1)。經(jīng)過計算得出此重復(fù)結(jié)構(gòu)的近似摩爾比為x∶y∶z=5∶3∶1,以此為后續(xù)有關(guān)NBR 結(jié)構(gòu)設(shè)計和理論建模的基礎(chǔ)。

    1.2.4 NBR 催化加氫反應(yīng) 加氫催化劑采用實驗室自制的Pd/SiO2催化劑,載體的制備和表面修飾、負(fù)載型催化劑的制備方法與文獻[26]一致,不再贅述。

    圖1 NBR的1H NMR譜圖Fig.1 1H NMR spectrum of NBR

    表1 NBR中不同結(jié)構(gòu)對應(yīng)的H 化學(xué)位移Table 1 Chemical shifts of H protons in different microstructures of NBR

    將1 g NBR 溶解于80 g 溶劑中,攪拌均勻,得到NBR 膠液。所用溶劑包括:TOL、OX、TCM、MCB、EA、MAC、CYC、MEK、ACE。將自制催化劑與NBR膠液混合后移入到高壓反應(yīng)釜中密封;通入N2將反應(yīng)釜內(nèi)的空氣排除干凈,再通入H2(氫氣壓力為1 MPa)。在控制臺顯示器上設(shè)置如下反應(yīng)參數(shù):反應(yīng)溫度為60℃、反應(yīng)時間為2 h 及攪拌速率為800 r/min。反應(yīng)結(jié)束后,待反應(yīng)釜內(nèi)溫度降低至室溫,排空后開釜取樣。反應(yīng)后的膠液和催化劑離心分離,固體催化劑回收,膠液通過適量乙醇聚沉,得到固態(tài)HNBR,在85℃烘箱中干燥24 h,用于后續(xù)表征測試。NBR 及HNBR 的化學(xué)結(jié)構(gòu)表征在Thermo Fisher Scientific 公司制造的Nicolet iS50 型紅外光譜儀上獲取,通過FT-IR譜圖分析得到—CN基團、C C及C—C的特征峰強度,計算可求出氫化度[26]。

    2 實驗結(jié)果與討論

    2.1 不同溶劑對NBR溶解的影響

    圖2 NBR在強親電子溶劑中放置24 h前(a)、后(b)對比照片(從左至右依次為:EtOH、EG、IPA、水和MeCN)Fig.2 Pictures for NBR in strong electrophilic solvents before(a)and after(b)24 h(Solvents from left to right:EtOH,EG,IPA,water and MeCN)

    根據(jù)溶劑的給、親電子能力,將溶劑分為三類:強親電子溶劑、弱親電子溶劑及給電子溶劑,分別選擇幾種代表性的溶劑,研究溶劑性質(zhì)對NBR 溶解的影響。圖2為NBR在強親電子溶劑中放置24 h前后對比照片,可以看出NBR 在EtOH、EG、IPA 和水中放置24 h后,其體積和外形均未發(fā)生明顯變化[圖2(b)],說明在這些溶劑中NBR 既不能溶解也不能溶脹,原因可能在于上述四種溶劑均具有強親電子的官能團(—OH),而NBR 分子中也具有強親電子的官能團(—CN),根據(jù)“溶劑化原則”可知,NBR和上述強親電子溶劑不具備一方親電子和另一方給電子的條件,因而此類溶劑不能溶解NBR 分子。對于MeCN 溶劑,NBR 在放置24 h 后[圖2(b)]呈現(xiàn)發(fā)白、體積變大的現(xiàn)象,說明NBR 在MeCN 中只能溶脹不能溶解,這可能是由于MeCN 和NBR 分子中均含有—CN,MeCN 可以滲透到部分NBR 分子中,但兩者不滿足高分子溶解的“溶劑化原則”。此外,NBR 在乙醇、IPA 及MeCN 溶劑中置于溶劑底部,說明NBR 的密度大于上述三類溶劑;而在EG 和水中NBR 處于其上部,表明NBR 的密度小于EG 和水的密度。另外,強親電子溶劑與NBR 的溶解度(20.0 (J/cm3)0.5)參數(shù)相差較大,也是其不能溶解NBR的另一原因。

    圖3為NBR在弱親電子溶劑中放置24 h前后對比照片,可以看出NBR 在TOL、二甲苯(OX、MX、PX)、TIPB及EDC溶劑中放置24 h后,與初始階段相比,橡膠塊消失,溶液呈現(xiàn)較為渾濁的狀態(tài);對MCB、DCM 及TCM 溶劑而言,經(jīng)過24 h 后橡膠塊完全消失,溶液呈現(xiàn)清澈的狀態(tài),說明這3 種溶劑是NBR 的良溶劑。TOL、二甲苯及TIPB 為非極性溶劑(ε <2.8),且此類溶劑與NBR 的溶解度參數(shù)絕對差值(Δδ=|δp-δs|)均大于1.5,根據(jù)“溶解度參數(shù)相近”和“極性相近”的原則,這些溶劑對NBR 溶解度不高,因此膠液呈現(xiàn)渾濁狀態(tài);MCB、DCM 與TCM 為極性溶劑(ε >3.6),且其溶解度參數(shù)與NBR 的溶解度參數(shù)值非常接近,有利于NBR 的溶解,溶液最終呈現(xiàn)澄清透明狀態(tài)。

    圖3 NBR在弱親電子溶劑中放置24 h前(a)、后(b)對比照片(從左至右溶劑依次為:TOL、OX、MX、PX、TIPB、MCB、DCM、TCM和EDC)Fig.3 Pictures for NBR in weak electrophilic solvents before(a)and after(b)24 h(Solvents from left to right:TOL,OX,MX,PX,TIPB,MCB,DCM,TCM and EDC)

    圖4為NBR在給電子溶劑中放置24 h前后對比照片,可以看出經(jīng)過24 h 后NBR 橡膠塊在所有給電子溶劑中均消失,溶液澄清,說明給電子溶劑能夠很好地溶解NBR。而這些給電子溶劑的溶解度參數(shù)卻與NBR 溶解度參數(shù)大相徑庭,暗示了給電子溶劑與NBR 之間的溶解規(guī)律不能單從“溶解度參數(shù)相近”方面考慮,還需要考慮溶劑與NBR 之間的“溶劑化作用”,即NBR 和溶劑具備一方親電子和另一方給電子的條件,因而可以很好地溶解NBR。

    圖4 NBR在給電子溶劑中放置24 h前(a)、后(b)對比照片(從左至右溶劑依次為:ACE、MEK、MIBK、CYC、NMP、MAC、EA、IBAC和DMF)Fig.4 Pictures for NBR in electron donating solvents before(a)and after(b)24 h(Solvents from left to right:ACE,MEK,MIBK,CYC,NMP,MAC,EA,IBAC and DMF)

    給電子溶劑是NBR 溶解的良溶劑,實驗進一步選取了三種對NBR 有較好溶解能力的酮類溶劑(ACE、CYC 和MIBK),在相同溶解條件下研究NBR在不同溶劑中的溶解快慢。實驗結(jié)果表明NBR 在這三種酮類溶劑中的溶解速度呈現(xiàn)如下規(guī)律:ACE>MIBK>CYC。由于溶劑分子極性強弱會影響溶質(zhì)分子與溶劑分子間相互作用力,因而NBR 溶解速度快慢可能與溶劑的極性有關(guān),而分子的偶極矩數(shù)值可以反映分子極性的強弱,溶劑分子偶極矩排序為:丙酮(2.88 D,1 D=3.33564×10-30C·m) >4-甲基-2-戊酮(2.69 D) >環(huán)己酮(0.53 D),因此在一定范圍內(nèi),溶劑分子的偶極矩越大,極性越強,其對NBR 的溶解越快。

    2.2 不同溶劑中NBR的尺寸分布

    NBR 在弱親電子溶劑中的尺寸分布如圖5(a)所示(濃度為1 g/L),NBR 在DCM、TCM、MCB 及TOL 中粒徑尺寸與分布范圍相差不大(4~11 nm),在DCM、TCM 和MCB 中的分散系數(shù)(PDI,值越小越接近單分散)小于0.3,表明在這三種溶劑中NBR 的粒徑分布較為集中,該測試結(jié)果可信度較高;在TOL 中,NBR的粒徑雖然較小,但PDI>0.3,表明NBR 并非單分散,溶液中還有為數(shù)不多但粒徑較大的NBR 顆粒存在;在二甲苯(OX、MX、PX)溶劑中,NBR 的粒徑較大且分布范圍寬,表明NBR 并未完全溶解,以較大顆粒的形態(tài)存在;在EDC 和TIPB 溶劑中,NBR 粒徑及其分布范圍均不大,但兩者的PDI 較大(>0.3),說明NBR 在這兩種溶劑中非單分散,仍有大顆粒的NBR存在溶液體系中。由給電子溶劑中NBR 的尺寸分布情況可知[圖5(b)],NBR 的粒徑及其分布范圍非常接近,其粒徑約為10 nm。結(jié)合溶解實驗可知,NBR在良溶劑中時,溶解度高,膠液澄清,NBR 粒徑分布集中,粒徑較?。辉诓涣既軇┲袝r,膠液渾濁,NBR粒徑較大且分布寬。

    2.3 不同溶劑中NBR膠液的黏度

    計算所得不同溶劑中NBR的增比黏度(ηsp)及[η]列于表2。由表可知,在弱親電子溶劑中,NBR 在DCM 和TCM 溶劑中的[η]和ηsp均大于在其余溶劑中所得的兩種黏度,說明NBR 與DCM、TCM 之間的相互作用力較強,分子鏈段在該兩種溶劑中最為舒展;NBR 在二甲苯和TIPB 溶劑中的黏度較小,表明NBR 與這些溶劑相互作用力較小,因此溶解不充分,膠液呈現(xiàn)渾濁狀態(tài)(圖3)且粒徑較大(圖5)。對于給電子溶劑,NBR 在NMP 溶劑中的ηsp和[η]兩種黏度均最大,其次是CYC 溶劑,說明NMP 和CYC 與NBR 之間的相互作用力較大,NBR 在這兩種溶劑中更為舒展。

    圖5 弱親電子溶劑(a)和給電子溶劑(b)中NBR的尺寸分布(a1~a9分別為:DCM、TCM、EDC、MCB、TOL、OX、MX、PX和TIPB;b1~b9分別為:ACE、MEK、MIBK、CYC、NMP、MAC、EA、IBAC和DMF)Fig.5 The size distribution of NBR in weak electrophilic solvents(a)and electron donating solvents(b)(Solvents a1—a9:DCM,TCM,EDC,MCB,TOL,OX,MX,PX and TIPB;solvents b1—b9:ACE,MEK,MIBK,CYC,NMP,MAC,EA,IBAC and DMF)

    NBR 在溶劑中的[η]雖然能夠直觀地反映溶劑和NBR 之間相互作用力的強弱,但其數(shù)值較??;ηsp雖然數(shù)值較大,但同時包含了溶劑-NBR、溶劑-溶劑之間相互作用力的結(jié)果。將[η]與ηsp作圖,如圖6(a)和(b)所示,可知無論在弱親電子溶劑中還是在給電子溶劑中,NBR 的兩種黏度變化趨勢一致,且不同ηsp之間的差值更為直觀。從NBR 溶液增比黏度與溶解度參數(shù)的絕對差值的關(guān)系看,弱親電子溶劑與NBR 的Δδ 越小,溶液的黏度越大[圖6(c)],說明溶劑與NBR 之間的相互作用力越強,NBR 分子在對應(yīng)溶液中溶解性越好。但對給電子溶劑而言[圖6(d)],溶劑與NBR 的Δδ 與溶液體系的黏度無明顯線性變化的規(guī)律,進一步說明了給電子溶劑對NBR 的影響不能通過溶解度參數(shù)等宏觀物化角度進行表述,需要考慮微觀分子層面上的溶質(zhì)與溶劑間“溶劑化能”作用。

    表2 弱親電子(a1~a9)和給電子溶劑(b1~b9)中NBR溶液的黏度Table 2 Viscosity of NBR solutions in weak electrophilic solvents(a1—a9)and electron donating solvents(b1—b9)

    2.4 不同溶劑對NBR結(jié)構(gòu)的影響:DFT模擬計算

    圖6 弱親電子溶劑(a)和給電子溶劑(b)中NBR溶液特性黏度和增比黏度之間的關(guān)系;弱親電子溶劑(c)和給電子溶劑(d)中NBR溶液增比黏度與溶解度參數(shù)的絕對差值的關(guān)系(溶解度參數(shù)的絕對差值Δδ=|δNBR-δsolvent|)Fig.6 The relationship between the characteristic viscosity and the specific viscosity of NBR solutions in weak electrophilic solvents(a)and electron donating solvents(b);the relationship between specific viscosity of NBR solutions and absolute difference of solubility parameters in weak electrophilic solvents(c)and electron donating solvents(d)(absolute difference of solubility parameters Δδ=|δNBR-δsolvent|)

    2.4.1 不同溶劑對NBR 鏈段中C C 的影響 NBR分子的結(jié)構(gòu)根據(jù)核磁實驗(圖1)測得的1,4-丁二烯、丙烯腈和1,2-丁二烯三種結(jié)構(gòu)單元的近似比5∶3∶1作為建模的依據(jù),模擬了NBR 初始可能的4 種構(gòu)型(圖7),并根據(jù)頻率分析結(jié)果中無虛頻的出現(xiàn)來確定其為勢能面上極小點。通過計算得到了此4種構(gòu)型的能量(E),分析數(shù)據(jù)可知,在所有模擬計算的NBR結(jié)構(gòu)中,圖7(a)所示結(jié)構(gòu)的能量最低(-1450.0695 a.u.),說明圖7(a)所示結(jié)構(gòu)是這4種同分異構(gòu)體中相對最穩(wěn)定的構(gòu)型,因此選用圖7(a)所示結(jié)構(gòu)來研究和討論溶劑對NBR結(jié)構(gòu)的影響。

    NBR 的穩(wěn)定構(gòu)型[7(a)]中不同原子序號、雙鍵及腈基的分布情況如圖8所示,基于穩(wěn)定構(gòu)型[7(a)],使用SMD 溶劑化模型來研究溶劑效應(yīng),選擇代表性弱親電子溶劑(TOL、OX、TCM、MCB)和給電子溶劑(EA、MAC、CYC、MEK、ACE)共9 種不同性質(zhì)的溶劑進行了優(yōu)化和計算。對NBR 分子在氣相和不同溶劑中的鍵長、軌道能量、溶劑化能和偶極矩進行了統(tǒng)計分析。

    表3 列出了NBR 分子在氣相和不同溶劑中的1,4-C C—和1,2-C C—雙鍵鍵長,并根據(jù)數(shù)據(jù)得到NBR 鏈段不同位置的雙鍵鍵長隨著介電常數(shù)的變化曲線,如圖9所示。從氣相中的數(shù)據(jù)可以看出,不同位置的C C 雙鍵鍵長存在差異,3 位置的1,4-C C—鍵長最長,其次是1 位置的1,4-C C—鍵長,4,5 位置的1,4-C C—鍵長接近,最短的是2 位置的1,4-C C—鍵和1,2-C C—雙鍵的鍵長。對比不同溶劑中NBR 分子上C C 雙鍵的鍵長可知,無論是弱親電子還是給電子溶劑,1,4-C C—和1,2-C C—雙鍵鍵長均隨著溶劑的極性增加而呈現(xiàn)增長的趨勢,對比在氣相和不同溶劑中C C鍵長的計算結(jié)果(表3),給電子溶劑中雙鍵鍵長增長程度大于在弱親電子溶劑中,表明在具有較大極性的給電子溶劑中NBR 鏈段上C C 雙鍵的鍵能更小,更有利于雙鍵的活化,換言之,在該C C 上發(fā)生化學(xué)反應(yīng)的活性可能越大。

    圖7 氣相中NBR鏈段的幾何模型(白色:H,灰色:C,藍色:N)Fig.7 NBR geometric patterns in gas phase(White:H,gray:C,and blue:N)

    圖8 NBR鏈段中C C和C N雙鍵的分布情況Fig.8 The C C and C N bonds distribution in NBR chain

    表3 不同溶劑中NBR鏈段里不同位置C C雙鍵鍵長Table 3 The length of C C double bonds at different positions of NBR molecule in different solvents

    通過DFT 模擬得到不同溶劑中NBR 鏈段1,4-C C—紅外伸縮振動峰的峰位置繪制于圖10。氣相中1,4-C C—振動峰的出峰位置是974.04 cm-1,在弱親電子溶劑氛圍中[圖10(a)],對于NBR 溶解能力不好的弱親電子溶劑TOL 和OX,1,4-C C—伸縮振動峰的峰位置向高波數(shù)方向移動,而在對NBR溶解能力較好的弱親電子溶劑TCM 和MCB 及給電子溶劑氛圍中[圖10(b)],該峰的峰位置均往低波數(shù)方向移動,在給電子的酮類溶劑中移動最顯著,其峰位置均低于969.82 cm-1,說明酮類給電子溶劑更有利于1,4-C C—雙鍵的活化。

    圖9 弱親電子溶劑(a)和給電子溶劑(b)中NBR鏈段不同位置的C C雙鍵的鍵長Fig.9 The bond length of the C C double bonds at different positions of NBR molecule in weak electrophilic solvents (a)and electron donating solvents(b)

    圖10 弱親電子溶劑(a)和給電子溶劑(b)中NBR鏈段里1,4-C C—雙鍵紅外伸縮振動峰位置Fig.10 The stretching vibration peak position of the 1,4-C C—double bond of NBR chain in weak electrophilic solvents(a)and electron donating solvents(b)

    2.4.2 不同溶劑對NBR 鏈段前線軌道的影響 由Fukui等[27]提出的前線軌道理論屬于分子軌道理論,其中有電子占據(jù)的高能級軌道稱為HOMO 軌道,沒有電子占據(jù)的低能級軌道稱為LUMO軌道。前線軌道是確定系統(tǒng)穩(wěn)定性和化學(xué)反應(yīng)活性中心的關(guān)鍵。LUMO 和HOMO 軌道以及它們之間的能級差(能隙值)[28]決定和影響了分子的許多物理化學(xué)性質(zhì),溶劑效應(yīng)是影響NBR 分子前線軌道的一個因素,通過理論計算獲得了NBR 分子在不同溶劑中前線軌道的空間分布并得到了軌道能級和能級差,氣相中NBR鏈段的HOMO 和LUMO 軌道分布如圖11 所示,其HOMO 和LUMO 軌道主要集中在整個分子骨架的C C 雙鍵以及腈基上,其中C C 雙鍵對HOMO 和LUMO貢獻較大。

    圖11 氣相NBR鏈段的前線軌道圖Fig.11 Frontier orbital of NBR chain in gas phase

    圖12 為不同溶劑中NBR 鏈段的HOMO 軌道分布情況,可以看出在弱親電子溶劑中HOMO 軌道分布與氣相中分布情況一致[圖12(a)],未隨溶劑的改變發(fā)生變化,表明弱親電子溶劑并未對其化學(xué)反應(yīng)活性中心造成影響。而在給電子溶劑ACE、MEK 及CYC 中,NBR 鏈段的HOMO 軌道分布均發(fā)生了變化,由邊緣位置的(5)1,4-C C—雙鍵和1,2-C C—雙鍵轉(zhuǎn)到分子鏈內(nèi)部(3)1,4-C C—雙鍵和(4)1,4-C C—雙鍵的位置[圖12(b)],表明受溶劑效應(yīng)的影響其化學(xué)反應(yīng)活性中心由分子邊緣向分子內(nèi)部移動,對NBR 鏈段內(nèi)部C C雙鍵產(chǎn)生積極的影響;而酯類給電子溶劑并未對NBR 鏈段的HOMO 軌道分布造成影響,與在氣相中相似。不同溶劑中NBR 鏈段的LUMO 軌道分布情況與氣相一致,溶劑的改變并未對其LUMO軌道分布產(chǎn)生影響。

    圖12 NBR在弱親電子溶劑(a)和給電子溶劑(b)中HOMO軌道分布圖Fig.12 The HOMO orbital distribution of the NBR molecule in weak electrophilic solvents(a)and electron donating solvents(b)

    2.4.3 不同溶劑中NBR 的軌道能隙、溶劑化能及偶極矩 表4 是不同溶劑中NBR 分子的前線軌道能量、能隙差及溶劑化能(ΔG),可以看出不同溶劑對NBR 分子的HOMO 和LUMO 軌道能級能量會產(chǎn)生一定影響。整體來看,除OX 外,酮類給電子溶劑中NBR 的HOMO 能級相對能量最低,可能是由于溶劑對其軌道分布的影響所致[由分子邊緣向分子內(nèi)部移動,圖12(b)]。HOMO 和LUMO 軌道能隙差值(ΔEgap)的大小能反映分子所在系統(tǒng)的穩(wěn)定性,NBR的ΔEgap值越大,NBR 分子就會越穩(wěn)定。由表中可以看出,在氣相中ΔEgap最低,分子體系能量較高,不穩(wěn)定,在溶劑氛圍中ΔEgap均高于氣相中,NBR 在酮類溶劑中ΔEgap值較大(0.24180~0.24310 a.u.),分子體系能量較低且穩(wěn)定。NBR 分子在溶劑氛圍中受到靜電和非靜電的溶劑化作用,其作用結(jié)果用溶劑化能ΔG來描述,計算所得溶劑化能均為負(fù)值,說明NBR在溶劑中的能量相對其氣相均有所降低。對比NBR 在弱親電子溶劑和給電子溶劑中的|ΔG|還發(fā)現(xiàn):酮類給電子溶劑(除CYC 外)>帶有鹵素原子的弱親電子溶劑>酯類給電子溶劑>芳香烴類弱親電子溶劑,說明NBR 與酮類給電子溶劑之間的溶劑化作用最強,與芳香烴類弱親電子溶劑之間的相互作用力最弱。此外,NBR 在OX 中溶劑化能的值(|ΔG|=18.69 kcal/mol, 1 cal=4.18 J)最小,說明NBR與OX 之間相互作用力最小,這也解釋了為什么NBR 在OX 溶劑中溶解較差、黏度較小、膠液渾濁且顆粒較大的原因。整體而言,NBR 在酮類溶劑中的|ΔG|(22.20~26.31 kcal/mol)較大,說明NBR 與酮類溶劑之間的相互作用力較大,有利于NBR 的溶解及在溶劑中的穩(wěn)定存在。

    偶極距是反映分子中電荷分布的物理量[29]。分子電荷不均勻分布的變化反映在偶極距的波動上,由于分子極性的大小與分子內(nèi)正負(fù)電荷分布不均勻有關(guān),從而可以通過分子偶極矩的變化來判斷整個分子極性的變化。分子極性大小會影響溶質(zhì)分子與溶劑分子間相互作用力,在一定程度上決定了溶解度[30]。因此,分析不同溶劑中NBR 分子偶極距有助于理解溶劑效應(yīng)對整個分子的極性及溶解度的影響。如表4 所示,NBR 分子的極性隨著溶劑極性的增大而增大,到一定程度后趨于穩(wěn)定,可以看出酮類溶劑對NBR 分子的極性影響最大,其對NBR的溶解性也最優(yōu)。

    2.5 不同溶劑對NBR加氫的影響

    選用上述不同性質(zhì)的9種溶劑分別溶解NBR分子后,加入自制催化劑,在高壓反應(yīng)釜中進行催化加氫反應(yīng),得到產(chǎn)品HNBR,其紅外譜圖及加氫度如圖13 所示。紅外譜圖中出現(xiàn)的重要特征峰2236 cm-1為—CN 基團的特征吸收峰,從圖中可以看出此處的吸收峰強度在加氫前后基本不變,且加氫后在3250~3350 cm-1處均未出現(xiàn)胺基的特征吸收峰,說明無論在何種溶劑中—CN 基團在反應(yīng)過程中均未被加氫。970 cm-1和920 cm-1處的吸收峰歸屬于1,4反式-CH CH—和1,2-CH CH—的特征吸收峰,明顯看出在不同的溶劑中1,2-CH CH—的吸收峰基本消失,1,4 反式-CH CH—峰強度有不同程度地降低,說明1,2-CH CH—加氫更容易。723 cm-1處的吸收峰為加氫后生成的-[CH2]n(n>4)的變形振動吸收峰,其峰強度在不同溶劑中也有差異。根據(jù)紅外法加氫度計算的方法[26],通過紅外譜圖中出現(xiàn)的特征吸收峰,可得出不同溶劑中所得HNBR 的加氫度,結(jié)果如圖13(b)和(d)所示。在弱親電子溶劑(TCM、MCB、TOL、OX)中得到的HNBR 加氫度分別為:72%,69%,76%和79%;在給電子溶劑(ACE、MEK、CYC、MAC、EA)中得到的HNBR 加氫度分別為:93%,89%,94%,85%和90%。結(jié)果表明,NBR在弱親電子溶劑中的加氫度均低于80%,而在給電子溶劑中的加氫度均高于80%,即給電子溶劑更有利于NBR的加氫轉(zhuǎn)化。結(jié)合溶劑對NBR溶解、尺寸和結(jié)構(gòu)的影響結(jié)果,NBR 在溶解較好的給電子溶劑中分子尺寸粒徑更小,特別是在酮類溶劑中溶劑化能較大、C C 雙鍵更易活化,因而NBR 在給電子溶劑中的加氫效果更優(yōu)。

    表4 不同溶劑中NBR前線軌道能量、能隙值、溶劑化能及偶極矩Table 4 Frontier orbit energy,energy gap,solvation energy and dipole moments of NBR in different solvents

    3 結(jié) 論

    圖13 在弱親電子溶劑[(a),(b)]和給電子溶劑[(c),(d)]中得到的HNBR紅外譜圖及其加氫度Fig.13 FT-IR spectra and hydrogenation degree of HNBR obtained in weak electrophilic solvents [(a),(b)]and electron donating solvents[(c),(d)]

    本研究系統(tǒng)考察了不同性質(zhì)的溶劑對NBR 的溶解、尺寸、結(jié)構(gòu)和催化加氫的影響,結(jié)果表明能夠較好溶解NBR 的良溶劑需要符合“極性相近”及“溶劑化能”的規(guī)則,給電子溶劑和部分帶有鹵素原子的弱親電子溶劑是NBR 的良溶劑。NBR 在良溶劑中溶解度高、膠液澄清、粒徑小且分布集中。溶液特性黏度和增比黏度之間的關(guān)系說明弱親電子溶劑與NBR 的溶解度參數(shù)的絕對差值越小,溶液的黏度越大,NBR 分子在對應(yīng)溶劑中溶解性越好,而給電子溶劑需要考慮微觀分子層面上的溶質(zhì)與溶劑間“溶劑化能”作用。DFT 計算結(jié)果表明對NBR 溶解能力較好的弱親電子溶劑以及給電子酮類溶劑都會導(dǎo)致1,4-C C—伸縮振動峰的峰位置發(fā)生紅移,NBR 鏈段中雙鍵鍵長隨著溶劑的極性增加而呈增長的趨勢。前線軌道理論計算結(jié)果表明在酮類溶劑中,NBR 鏈段的HOMO 軌道分布發(fā)生明顯變化,由分子邊緣向分子內(nèi)部移動,對鏈段內(nèi)部C C雙鍵產(chǎn)生積極影響。催化加氫結(jié)果表明在給電子溶劑中比在弱親電子溶劑中更有利于NBR 的加氫轉(zhuǎn)化,雖然不同溶劑對NBR 加氫轉(zhuǎn)化率的影響各不相同,但對NBR 中C C 加氫選擇性均為100%。本研究對NBR 非均相溶液催化加氫制備高附加值、高性能HNBR 過程中溶劑的篩選具有重要的指導(dǎo)意義,也為研究其他聚合物在不同溶劑中的溶解行為、結(jié)構(gòu)影響提供一定的借鑒。

    猜你喜歡
    雙鍵鏈段溶解度
    ●一種服用阻燃聚酯纖維及其制備方法(申請?zhí)朇N202011602594.5 公開日 2021-12-17 申請人 江蘇恒力化纖股份有限公司)
    軟硬鏈段含量及相對分子質(zhì)量對PA 6基彈性體結(jié)構(gòu)的影響
    “溶解度曲線”考點擊破
    復(fù)雜斷塊調(diào)剖提效研究
    《有機化學(xué)》課程中不對稱烯烴親電加成反應(yīng)教學(xué)研究
    光化學(xué)反應(yīng)—串聯(lián)質(zhì)譜法鑒定細胞中不飽和卵磷脂雙鍵的位置
    溶解度曲線的理解與應(yīng)用例析
    CO2捕集的吸收溶解度計算和過程模擬
    丁腈橡膠中丙烯腈鏈段的分布特性對其與三元乙丙橡膠共混的影響
    溶解度計算錯誤種種和對策
    亚洲第一区二区三区不卡| 精品国产一区二区三区久久久樱花 | 国产高清有码在线观看视频| 国产麻豆成人av免费视频| 99久久人妻综合| 91狼人影院| 国产爱豆传媒在线观看| 97在线视频观看| 亚洲欧美日韩高清专用| 国产探花在线观看一区二区| av在线蜜桃| 高清日韩中文字幕在线| 国产精华一区二区三区| 亚洲天堂国产精品一区在线| 夜夜看夜夜爽夜夜摸| 色综合亚洲欧美另类图片| 少妇裸体淫交视频免费看高清| 三级毛片av免费| 免费观看精品视频网站| 日韩一本色道免费dvd| 国产成人精品婷婷| 春色校园在线视频观看| 大又大粗又爽又黄少妇毛片口| 亚洲中文字幕日韩| 免费电影在线观看免费观看| 禁无遮挡网站| 成人一区二区视频在线观看| 欧美一区二区精品小视频在线| 爱豆传媒免费全集在线观看| 亚洲欧美日韩高清专用| 可以在线观看毛片的网站| 国产极品精品免费视频能看的| 99热网站在线观看| 国产精品嫩草影院av在线观看| videos熟女内射| 亚洲成色77777| 午夜a级毛片| 两个人视频免费观看高清| 男人狂女人下面高潮的视频| 少妇高潮的动态图| 亚洲av.av天堂| 午夜福利在线观看吧| 高清毛片免费看| 九草在线视频观看| 亚洲av男天堂| 国产淫片久久久久久久久| 精品久久久噜噜| 日本黄色视频三级网站网址| 精品午夜福利在线看| 男插女下体视频免费在线播放| 久久精品国产鲁丝片午夜精品| 白带黄色成豆腐渣| 黑人高潮一二区| 99久久精品国产国产毛片| 老师上课跳d突然被开到最大视频| 国语自产精品视频在线第100页| 久久久久精品久久久久真实原创| 国产精品久久视频播放| 青春草视频在线免费观看| 国产极品精品免费视频能看的| 国产欧美另类精品又又久久亚洲欧美| 久久久色成人| 国产又色又爽无遮挡免| www日本黄色视频网| 校园人妻丝袜中文字幕| 久久精品国产亚洲av天美| 国产毛片a区久久久久| 欧美日韩精品成人综合77777| 一夜夜www| 如何舔出高潮| 欧美区成人在线视频| 精品国产三级普通话版| 男人的好看免费观看在线视频| 日韩欧美国产在线观看| 亚洲av日韩在线播放| 亚洲国产欧美在线一区| 亚洲欧美中文字幕日韩二区| 91精品一卡2卡3卡4卡| 亚洲精品久久久久久婷婷小说 | 亚洲三级黄色毛片| 蜜臀久久99精品久久宅男| 中国国产av一级| 久久精品国产亚洲网站| 亚洲无线观看免费| 热99在线观看视频| 国产亚洲精品久久久com| 亚洲人成网站高清观看| 又粗又爽又猛毛片免费看| 丰满乱子伦码专区| 中文字幕亚洲精品专区| 在线天堂最新版资源| 亚洲av.av天堂| 午夜精品在线福利| 99热精品在线国产| 在线播放无遮挡| 69人妻影院| 视频中文字幕在线观看| 久久综合国产亚洲精品| 久久久久久久午夜电影| 亚洲内射少妇av| 我的女老师完整版在线观看| 欧美日韩在线观看h| 久久这里只有精品中国| 亚洲精品日韩在线中文字幕| 男人舔奶头视频| 午夜日本视频在线| 亚洲精品久久久久久婷婷小说 | 亚洲精品,欧美精品| 成年版毛片免费区| 国产高清三级在线| 亚洲在久久综合| 三级国产精品欧美在线观看| 国产精华一区二区三区| 蜜桃亚洲精品一区二区三区| 亚洲精品一区蜜桃| 色5月婷婷丁香| 国产精品综合久久久久久久免费| 久久久成人免费电影| 床上黄色一级片| 舔av片在线| 美女xxoo啪啪120秒动态图| 一边摸一边抽搐一进一小说| 18禁动态无遮挡网站| 成人二区视频| 日韩一区二区视频免费看| .国产精品久久| 美女大奶头视频| 国产精品三级大全| 毛片女人毛片| 久久精品国产鲁丝片午夜精品| 亚洲最大成人中文| 国产午夜福利久久久久久| 在线观看一区二区三区| 我要看日韩黄色一级片| 亚洲av免费高清在线观看| 免费黄网站久久成人精品| 国产成人午夜福利电影在线观看| 亚洲欧美精品自产自拍| 最新中文字幕久久久久| 日韩成人av中文字幕在线观看| 国产美女午夜福利| av专区在线播放| av在线蜜桃| 狂野欧美激情性xxxx在线观看| 国产精品蜜桃在线观看| 国产女主播在线喷水免费视频网站 | 少妇的逼好多水| 少妇熟女欧美另类| 在线观看66精品国产| 五月玫瑰六月丁香| 桃色一区二区三区在线观看| 中文亚洲av片在线观看爽| 亚洲欧洲日产国产| 亚洲不卡免费看| 精品熟女少妇av免费看| 偷拍熟女少妇极品色| 成人亚洲精品av一区二区| 欧美区成人在线视频| 久久久久久国产a免费观看| 久久精品久久久久久久性| 黄色日韩在线| 我要看日韩黄色一级片| 国产在线一区二区三区精 | 狠狠狠狠99中文字幕| 日韩亚洲欧美综合| 非洲黑人性xxxx精品又粗又长| 大又大粗又爽又黄少妇毛片口| 午夜视频国产福利| 两个人的视频大全免费| 3wmmmm亚洲av在线观看| 亚洲综合色惰| 国产成年人精品一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 一个人免费在线观看电影| 免费大片18禁| 丝袜美腿在线中文| 午夜激情欧美在线| 18禁动态无遮挡网站| 国产精品永久免费网站| 两性午夜刺激爽爽歪歪视频在线观看| 汤姆久久久久久久影院中文字幕 | 精华霜和精华液先用哪个| 精品久久久久久久久亚洲| 国产高清国产精品国产三级 | 中文乱码字字幕精品一区二区三区 | 男人的好看免费观看在线视频| av线在线观看网站| 99久久九九国产精品国产免费| 国产精品福利在线免费观看| 亚洲精品国产成人久久av| 亚洲成人久久爱视频| 大又大粗又爽又黄少妇毛片口| 大话2 男鬼变身卡| 亚洲av免费在线观看| 亚洲欧美清纯卡通| 亚洲天堂国产精品一区在线| 久久精品91蜜桃| av又黄又爽大尺度在线免费看 | 色尼玛亚洲综合影院| 亚洲国产精品合色在线| 我的女老师完整版在线观看| 欧美成人精品欧美一级黄| 99视频精品全部免费 在线| 久久久精品94久久精品| 最近最新中文字幕大全电影3| 男人舔奶头视频| 美女内射精品一级片tv| 国产伦一二天堂av在线观看| 九色成人免费人妻av| 亚洲av不卡在线观看| 亚洲中文字幕一区二区三区有码在线看| 欧美性感艳星| 国产探花在线观看一区二区| 亚洲国产精品专区欧美| 国产精品伦人一区二区| 亚洲av中文av极速乱| 欧美xxxx黑人xx丫x性爽| 国产激情偷乱视频一区二区| 在线观看一区二区三区| 最近2019中文字幕mv第一页| 高清午夜精品一区二区三区| av在线天堂中文字幕| 尤物成人国产欧美一区二区三区| 国产一区二区三区av在线| 春色校园在线视频观看| h日本视频在线播放| 亚洲国产欧美人成| 亚洲av成人精品一区久久| 深夜a级毛片| 午夜精品在线福利| 国语自产精品视频在线第100页| 国产成年人精品一区二区| 麻豆成人午夜福利视频| 午夜精品国产一区二区电影 | 国产精华一区二区三区| 日本三级黄在线观看| 精品国产露脸久久av麻豆 | 又爽又黄a免费视频| 久久鲁丝午夜福利片| 26uuu在线亚洲综合色| 亚洲婷婷狠狠爱综合网| 亚洲人成网站高清观看| 久久久久久九九精品二区国产| 欧美一区二区国产精品久久精品| 日韩视频在线欧美| 大话2 男鬼变身卡| av免费在线看不卡| 超碰av人人做人人爽久久| 国产人妻一区二区三区在| 欧美+日韩+精品| 亚洲欧美成人综合另类久久久 | 国产精品久久久久久久电影| 国产亚洲最大av| 亚洲不卡免费看| 精品一区二区三区人妻视频| 我要搜黄色片| 看片在线看免费视频| 免费观看a级毛片全部| 1024手机看黄色片| 国产av码专区亚洲av| 嫩草影院精品99| 亚洲图色成人| 精品久久久久久久久亚洲| 99在线视频只有这里精品首页| 欧美人与善性xxx| 高清日韩中文字幕在线| 老司机福利观看| 国产女主播在线喷水免费视频网站 | 亚洲欧美精品专区久久| 国产亚洲av嫩草精品影院| 久久亚洲国产成人精品v| www.色视频.com| 2021少妇久久久久久久久久久| 在线播放无遮挡| 久久久久久久久大av| 午夜激情福利司机影院| 天堂中文最新版在线下载 | 欧美一区二区精品小视频在线| 欧美xxxx黑人xx丫x性爽| 特级一级黄色大片| 国产高清国产精品国产三级 | 日本爱情动作片www.在线观看| 日韩一区二区视频免费看| 最近最新中文字幕免费大全7| 国产精品久久久久久久久免| 国产亚洲5aaaaa淫片| 国产老妇伦熟女老妇高清| 三级男女做爰猛烈吃奶摸视频| 校园人妻丝袜中文字幕| 国产伦在线观看视频一区| 99在线人妻在线中文字幕| 亚洲精品一区蜜桃| 日韩制服骚丝袜av| 亚洲人成网站在线播| 2021天堂中文幕一二区在线观| 简卡轻食公司| ponron亚洲| 最近手机中文字幕大全| 国内精品一区二区在线观看| 国产黄色小视频在线观看| 免费一级毛片在线播放高清视频| 亚洲av成人av| 大又大粗又爽又黄少妇毛片口| 久久久久久伊人网av| 久久久久久久国产电影| 久久精品久久久久久噜噜老黄 | 99久久中文字幕三级久久日本| 欧美极品一区二区三区四区| 久久国产乱子免费精品| 国产精品久久久久久久久免| av播播在线观看一区| 在线免费十八禁| 亚洲最大成人中文| 国产av不卡久久| 高清毛片免费看| 久久欧美精品欧美久久欧美| av国产久精品久网站免费入址| 尾随美女入室| 超碰av人人做人人爽久久| 国产成人午夜福利电影在线观看| 有码 亚洲区| 久久久久久久亚洲中文字幕| 国产精品人妻久久久影院| 日本三级黄在线观看| 欧美bdsm另类| 精品99又大又爽又粗少妇毛片| 亚洲自拍偷在线| 午夜激情福利司机影院| 久久久午夜欧美精品| 久99久视频精品免费| 在线天堂最新版资源| 午夜福利在线观看吧| 三级经典国产精品| 久久午夜福利片| 黄色一级大片看看| 一个人观看的视频www高清免费观看| 亚洲无线观看免费| 国产黄色小视频在线观看| 国产伦理片在线播放av一区| 欧美一区二区亚洲| 两性午夜刺激爽爽歪歪视频在线观看| 久久精品久久精品一区二区三区| 亚洲国产精品合色在线| 99久久无色码亚洲精品果冻| av福利片在线观看| 精品久久久久久久末码| 在线a可以看的网站| 国产免费视频播放在线视频 | 国产免费福利视频在线观看| 国产免费视频播放在线视频 | 国产高清国产精品国产三级 | 国产精品乱码一区二三区的特点| 最近最新中文字幕免费大全7| 嫩草影院入口| 尾随美女入室| 99在线人妻在线中文字幕| 亚洲国产成人一精品久久久| 日韩欧美 国产精品| 激情 狠狠 欧美| 成人美女网站在线观看视频| 在线a可以看的网站| 69人妻影院| 久久精品熟女亚洲av麻豆精品 | 国内精品美女久久久久久| 综合色丁香网| 国产大屁股一区二区在线视频| 精品人妻视频免费看| 丝袜美腿在线中文| 成人欧美大片| 国产免费一级a男人的天堂| 成人午夜高清在线视频| 亚洲在线观看片| 亚洲成人中文字幕在线播放| 麻豆成人午夜福利视频| 国产白丝娇喘喷水9色精品| 97在线视频观看| 国产在视频线在精品| 男插女下体视频免费在线播放| 一级毛片aaaaaa免费看小| 久久精品国产自在天天线| 成人性生交大片免费视频hd| 三级经典国产精品| 嫩草影院精品99| 亚洲在线自拍视频| 亚洲精品成人久久久久久| 老司机影院成人| 国产精品一区二区三区四区免费观看| 精品久久久久久久久av| 欧美成人午夜免费资源| 久久久午夜欧美精品| 亚洲综合精品二区| 国产免费男女视频| 成年女人永久免费观看视频| 九九爱精品视频在线观看| 亚洲精品aⅴ在线观看| 国产成人a∨麻豆精品| 色5月婷婷丁香| 亚洲av免费在线观看| 精品99又大又爽又粗少妇毛片| 男人舔奶头视频| 一级毛片我不卡| 国产探花极品一区二区| 国产一区二区在线av高清观看| 欧美成人午夜免费资源| 神马国产精品三级电影在线观看| 久久久久九九精品影院| 国产毛片a区久久久久| 国产精品日韩av在线免费观看| 中文在线观看免费www的网站| 国产一区二区三区av在线| 成人性生交大片免费视频hd| av天堂中文字幕网| 少妇裸体淫交视频免费看高清| 免费人成在线观看视频色| 久久久久久伊人网av| 极品教师在线视频| 久久热精品热| 性色avwww在线观看| 少妇高潮的动态图| 我要搜黄色片| 国产又色又爽无遮挡免| 国产精品一区二区三区四区久久| 麻豆一二三区av精品| or卡值多少钱| 九草在线视频观看| 国产在视频线精品| 日韩人妻高清精品专区| 亚洲一区高清亚洲精品| 国内揄拍国产精品人妻在线| 中文在线观看免费www的网站| 国产精品一区二区三区四区久久| 国产又色又爽无遮挡免| 欧美激情在线99| 精品国产一区二区三区久久久樱花 | 高清毛片免费看| 免费电影在线观看免费观看| 国产淫语在线视频| 美女xxoo啪啪120秒动态图| 少妇熟女aⅴ在线视频| 国产精品久久视频播放| 国产高潮美女av| 国产精品国产三级国产专区5o | 夜夜爽夜夜爽视频| 午夜福利网站1000一区二区三区| 麻豆成人午夜福利视频| 黄片无遮挡物在线观看| 亚洲精品一区蜜桃| 国产成人福利小说| 内射极品少妇av片p| 性插视频无遮挡在线免费观看| 久久精品91蜜桃| 国产精品美女特级片免费视频播放器| 日韩欧美精品v在线| 99热这里只有精品一区| 亚洲熟妇中文字幕五十中出| 欧美高清成人免费视频www| 嫩草影院新地址| 国产精品1区2区在线观看.| 美女cb高潮喷水在线观看| 久久久亚洲精品成人影院| 婷婷色av中文字幕| 亚洲av.av天堂| av专区在线播放| 视频中文字幕在线观看| 国产老妇女一区| 国产人妻一区二区三区在| 黄片无遮挡物在线观看| av黄色大香蕉| 成年女人永久免费观看视频| 久久鲁丝午夜福利片| 插阴视频在线观看视频| 又粗又爽又猛毛片免费看| 极品教师在线视频| 国产成年人精品一区二区| 97超碰精品成人国产| 日韩av在线大香蕉| 两个人的视频大全免费| 国产一级毛片在线| 97超视频在线观看视频| 久久亚洲国产成人精品v| 亚洲国产精品成人综合色| 国产精品久久久久久久电影| 久久久久性生活片| 午夜a级毛片| 亚洲在线观看片| 欧美性猛交╳xxx乱大交人| 一个人看的www免费观看视频| 老司机影院成人| 欧美日本亚洲视频在线播放| 日韩欧美国产在线观看| 国产精品爽爽va在线观看网站| 免费看光身美女| 日韩欧美精品免费久久| 黄色一级大片看看| 日韩欧美精品免费久久| 波多野结衣高清无吗| 亚洲久久久久久中文字幕| 亚洲av成人av| 日韩,欧美,国产一区二区三区 | 午夜福利在线在线| 久久精品国产亚洲av天美| 国产成人福利小说| 国产毛片a区久久久久| 久久久色成人| 哪个播放器可以免费观看大片| 国产又黄又爽又无遮挡在线| 亚洲av电影在线观看一区二区三区 | 亚洲不卡免费看| 国产成人福利小说| 色综合色国产| 国产淫片久久久久久久久| 久久国产乱子免费精品| 精品久久久噜噜| 美女xxoo啪啪120秒动态图| 欧美日韩综合久久久久久| 免费黄色在线免费观看| 精品国产三级普通话版| 天堂av国产一区二区熟女人妻| 久久99热这里只频精品6学生 | 国产精品久久视频播放| 国产三级中文精品| 亚洲美女搞黄在线观看| 男女国产视频网站| 三级男女做爰猛烈吃奶摸视频| 最近的中文字幕免费完整| 国产不卡一卡二| 国产精品电影一区二区三区| 日韩av在线大香蕉| 日韩人妻高清精品专区| 久久精品国产亚洲av天美| 一级毛片电影观看 | 日韩成人av中文字幕在线观看| 欧美日韩一区二区视频在线观看视频在线 | 看非洲黑人一级黄片| 99热精品在线国产| 亚洲国产精品成人综合色| 22中文网久久字幕| 亚洲久久久久久中文字幕| 亚洲欧洲日产国产| 亚洲成色77777| 免费黄色在线免费观看| 亚洲av男天堂| 色5月婷婷丁香| 午夜福利视频1000在线观看| 精品人妻视频免费看| 黄色欧美视频在线观看| 国产一区二区在线av高清观看| 18+在线观看网站| 综合色丁香网| 欧美一级a爱片免费观看看| 精品一区二区三区人妻视频| 日韩 亚洲 欧美在线| 亚洲精品久久久久久婷婷小说 | 国产成人精品一,二区| 亚洲自偷自拍三级| 欧美潮喷喷水| 久久精品综合一区二区三区| 精品少妇黑人巨大在线播放 | 亚洲欧洲日产国产| 欧美97在线视频| 成人特级av手机在线观看| 看黄色毛片网站| 一区二区三区四区激情视频| 亚洲中文字幕日韩| 99热6这里只有精品| 精品久久久久久久久亚洲| 丝袜喷水一区| 国产成人精品婷婷| 亚洲,欧美,日韩| 成人无遮挡网站| 纵有疾风起免费观看全集完整版 | 国产精品一区www在线观看| 亚洲欧美精品专区久久| 别揉我奶头 嗯啊视频| 九色成人免费人妻av| 91狼人影院| 日韩欧美精品免费久久| 91久久精品电影网| 久久精品国产亚洲网站| 少妇丰满av| 亚洲av成人精品一区久久| 亚洲国产欧美人成| 久久久久久久久久成人| 91精品伊人久久大香线蕉| 寂寞人妻少妇视频99o| 久久精品国产亚洲网站| av免费观看日本| 亚洲精品乱码久久久v下载方式| 国产精品一及| 99久久精品一区二区三区| 乱人视频在线观看| 国产精品.久久久| 韩国av在线不卡| 日韩欧美三级三区| 成人高潮视频无遮挡免费网站| 国产淫片久久久久久久久| 久久久成人免费电影| 欧美激情国产日韩精品一区| 久99久视频精品免费| 国产私拍福利视频在线观看| 国产免费男女视频| 日韩 亚洲 欧美在线| 丰满人妻一区二区三区视频av| 不卡视频在线观看欧美| 欧美xxxx性猛交bbbb| 七月丁香在线播放| 不卡视频在线观看欧美| 狂野欧美白嫩少妇大欣赏| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产麻豆成人av免费视频| 亚洲一级一片aⅴ在线观看| 人妻制服诱惑在线中文字幕| 韩国av在线不卡| 国产国拍精品亚洲av在线观看| 亚洲av二区三区四区|