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

等離子硅刻蝕及其工藝參數(shù)的多尺度優(yōu)化

2014-10-30 08:41:26閻軍等
計算機(jī)輔助工程 2014年5期
關(guān)鍵詞:優(yōu)化

閻軍等

摘要: 用CFDACE+和CFDTOPO分別對容性耦合等離子體反應(yīng)腔室放電和等離子硅刻蝕過程進(jìn)行仿真,討論不同射頻電壓和腔室條件對等離子體特性的影響.結(jié)果表明:隨著射頻電壓的升高,離子的通量增大;在低射頻電壓時,離子通量隨腔室壓強(qiáng)的升高而減小,而在高射頻電壓時趨勢則相反.用Kriging模型對影響刻蝕形貌的參數(shù)(腔室壓強(qiáng)和射頻電壓)進(jìn)行優(yōu)化,結(jié)果表明該優(yōu)化方法可以為工藝條件相近的刻蝕機(jī)設(shè)備的設(shè)計提供參考.

關(guān)鍵詞: 硅刻蝕; 容性耦合等離子體; 射頻電壓; 腔室壓強(qiáng); Kriging模型; 優(yōu)化

中圖分類號: TN405.98; TB115.1文獻(xiàn)標(biāo)志碼: B

Abstract: The discharge process in capacitively coupled plasma reaction chamber and the plasma silicon etch process are simulated by CFDACE+ and CFDTOPO. The effect of different radio frequency voltage and chamber conditions on plasma characteristics are discussed. The results show that, with the increase of radio frequency voltage, the flux of ion increases; the flux of ion decreases with the increase of chamber pressure when the radio frequency voltage is low, but the change trend shows in opposite direction while the radio frequency voltage is high. The Kriging model is used to optimize the parameters such as chamber pressure and ration frequency voltage, which have influences on the etching profile. The result shows that the optimization method is feasible, which can provide reference for the design of the etch device in similar process condition.

Key words: silicon etch; capacitively coupled plasma; radio frequency voltage; chamber press; Kriging model; optimization

引言

技術(shù)的不斷進(jìn)步對微電子機(jī)械系統(tǒng)工藝和器件的設(shè)計提出更高的要求[1],如元器件本身的尺寸減小:其表面特征尺寸已由原來的微米級上升到現(xiàn)在的納米級.傳統(tǒng)的濕法刻蝕由于各向異性差、均勻性差、不容易控制等缺點,逐漸被各向異性好、均勻性好、容易控制的干法刻蝕所取代.

在干法刻蝕中等離子體刻蝕應(yīng)用最廣泛,也是微納米加工能力最強(qiáng)的技術(shù),但是其目前主要依賴于刻蝕經(jīng)驗.由于試驗周期長、成本高[2],等離子體刻蝕的工藝仿真技術(shù)顯得尤為重要.工藝仿真技術(shù)可通過建立適當(dāng)?shù)哪P停玫侥繕?biāo)結(jié)構(gòu)的仿真效果,因此可以根據(jù)刻蝕機(jī)腔室的工藝相關(guān)參數(shù)預(yù)測刻蝕形貌,還能夠?qū)に噮?shù)進(jìn)行優(yōu)化,進(jìn)而得到理想的刻蝕形貌.仿真模擬不僅降低微電子機(jī)械系統(tǒng)加工工藝的試驗成本,提高生產(chǎn)效率,更為等離子體刻蝕的研究提供一定的參考.

1刻蝕形貌分析

1.1腔室模型分析

建立二維軸對稱腔室模型,見圖1,模擬最簡單的容性耦合等離子體放電:Cl2從入口進(jìn)入腔室內(nèi)部,被射頻電壓電離成總體呈電中性的等離子體.電場對離子加速,帶有一定能量和角度分布的等離子體轟擊刻蝕材料,既發(fā)生物理反應(yīng)又發(fā)生化學(xué)反應(yīng),在刻蝕材料表面形成所需的圖案.

2基于Kriging模型的優(yōu)化設(shè)計

通常,評價刻蝕結(jié)果質(zhì)量的參數(shù)有剖面的垂直度、刻蝕速率、刻蝕均勻性和刻蝕的選擇比等.如果刻蝕的形貌與理想形貌差別較大,即刻蝕垂直度不夠好,將直接影響由晶片組成的元器件的質(zhì)量,甚至不能使用.腔室壓強(qiáng)和射頻電壓是影響刻蝕形貌的重要參數(shù).以這2個參數(shù)作為變量,分析不同參數(shù)組合下的刻蝕形貌,并以刻蝕垂直度作為目標(biāo)函數(shù)進(jìn)行初步優(yōu)化.

2.1試驗點選取和優(yōu)化模型選擇

腔室壓強(qiáng)和射頻電壓是影響等離子體密度、能量和角度分布等特性的2個重要參數(shù),而等離子體的特性決定等離子體刻蝕的結(jié)果.以腔室壓強(qiáng)和射頻電壓作為變量,分析其對刻蝕垂直度的影響.

采用2個變量3水平的全因子設(shè)計方法,通過對不同參數(shù)組合下的數(shù)值模擬,得到9組參數(shù)組合下的刻蝕形貌(垂直度).射頻電壓的取值范圍為0~200 V,腔室壓強(qiáng)的范圍為0~5 Pa.通過調(diào)節(jié)不同宏觀參數(shù)得到微觀的理想刻蝕形貌,是具有大難度的多尺度優(yōu)化問題,所以采用代理模型的方法實現(xiàn)優(yōu)化目的.同時,該問題屬于非線性問題,可采用Kriging模型.Kriging模型是一種估計方差最小的無偏估計模型[811],該模型既可以用來解決各向同性問題,又可以用來解決各向異性問題,而且在解決非線性程度較高的問題時容易取得理想的擬合效果.

以圖5所示的刻蝕形貌與理想形貌的垂直度差作為目標(biāo)函數(shù),計算不同條件下得到的離子流量和目標(biāo)函數(shù),見表1.隨著射頻電壓的升高,離子的流量增大;在低射頻電壓時,離子流量隨腔室壓強(qiáng)的升高而減小,高射頻電壓下趨勢相反.

3結(jié)論

(1)等離子體刻蝕中Cl2刻蝕硅是比較簡單的一類,腔室壓強(qiáng)和射頻電壓不同程度地影響等離子體的流量和能量分布等特性.隨著射頻電壓的升高,離子的流量增大;在低射頻電壓時,離子流量隨腔室壓強(qiáng)的升高而減小,高射頻電壓下趨勢相反.

(2)總體來說,雖然優(yōu)化后的刻蝕形貌不是完全理想的直上直下的刻蝕形貌,但呈現(xiàn)各向異性很強(qiáng)、剖面光滑的形貌.將優(yōu)化后的參數(shù)組合再進(jìn)行刻蝕模擬,得到的目標(biāo)函數(shù)為65.5,并不是71.3,可能是由2個原因造成的:一是優(yōu)化模型選取不當(dāng),二是試驗點太少(個變量可以分為4個水平或更多).

通過對Cl2等離子體刻蝕硅的數(shù)值模擬與優(yōu)化,從科研應(yīng)用的角度提出一種等離子刻蝕機(jī)的優(yōu)化方法,可以為工藝條件相近的刻蝕機(jī)設(shè)備的研究和設(shè)計提供參考.參考文獻(xiàn):

[1]徐棟梁, 張鑒, 許曉琳. 微機(jī)械加工工藝的三維仿真模型及算法[J]. 電子科技, 2010, 23(8): 14.

XU Dongliang, ZHANG Jian, XU Xiaolin. Threedimensional models and algorithms for micromanufacture simulations[J]. Electronic Sci & Technol, 2010, 23(8): 14.

[2]王曉光, 朱曉明, 尹延昭. 等離子體刻蝕工藝的優(yōu)化研究[J]. 中國新技術(shù)新產(chǎn)品, 2010(14): 2224.

WANG Xiaoguang, ZHU Xiaoming, YIN Yanzhao. Optimization research of plasma etching process[J]. China New Technol and Products, 2010(14): 2224.

[3]朱炳金, 張國棟, 張向鋒. 電感耦合等離子體刻蝕InSb芯片工藝的研究[J]. 紅外技術(shù), 2009, 31(8): 467470.

ZHU Bingjin, ZHANG Guodong, ZHANG Xiangfeng. Inductive couple plasma etching processing of InSb wafer[J]. Infrared Technol, 2009, 31(8): 467470.

[4]陳曉南, 楊培林, 龐宣明, 等. 等離子體刻蝕工藝中工藝參數(shù)對刻蝕速率影響的研究[J]. 西安交通大學(xué)學(xué)報, 2004, 38(5): 546547.

CHEN Xiaonan, YANG Peilin, PANG Xuanming, et al. Influence of process parameters on the etching rate in inductively coupled plasma etcher[J]. J Xian Jiaotong Univ, 2004, 38(5): 546547.

[5]程嘉, 朱煜. ICP刻蝕機(jī)反應(yīng)腔室氣流仿真研究[J]. 半導(dǎo)體技術(shù), 2007, 32(1): 4346.

CHENG Jia, ZHU Yu. Gas flow simulation research on reaction chamber of ICP etcher[J]. Semiconductor Technol, 2007, 32(1): 4346.

[6]程嘉, 朱煜, 段廣洪, 等. 基于回歸正交設(shè)計的ICP刻蝕機(jī)工藝腔室流場特性分析[J]. 半導(dǎo)體學(xué)報, 2008, 29(4): 780784.

CHENG Jia, ZHU Yu, DUAN Guanghong, et al. Analysis of processing chamber flow field characteristics for an ICP etcher based on regression orthogonal design[J]. J Semiconductors, 2008, 29(4): 780784.

[7]SEO Seung T, LEE Yong H, LEE Kwang S, et al. Runtorun control of inductively coupled C2F6 plasma etching of SiO2: construction of a number process with a computational fluid dynamics code[EB/OL]. (20050705)[20130610]. http://2005.iccas.org/submission/paper/upload/ 12ICP_simulation.pdf.

[8]王凌, 吉利軍, 鄭大鐘. 基于代理模型和遺傳算法的仿真優(yōu)化研究[J]. 控制與決策, 2004, 19(6): 626630.

WANG Ling, JI Lijun, ZHENG Dazhong. Simulation optimization based on surrogate model and genetic algorithm[J]. Contr & Decision, 2004, 19(6): 626630.

[9]高月華, 王希誠. 基于Kriging代理模型的穩(wěn)健優(yōu)化設(shè)計[J]. 化工學(xué)報, 2010, 61(3): 676681.

GAO Yuehua, WANG Xicheng. Robust optimization based on Kriging surrogate model[J]. CIESC J, 2010, 61(3): 676681.

[10]謝延敏, 于滬平, 陳軍, 等. 基于Kriging模型的可靠度計算[J]. 上海交通大學(xué)學(xué)報, 2007, 41(2): 177193.

XIE Yanmin, YU Huping, CHEN Jun, et al. Reliability estimation based on Kriging model[J]. J Shanghai Jiaotong Univ, 2007, 41(2): 177193.

[11]隋允康, 宇慧平. 響應(yīng)面方法的改進(jìn)及其對工程優(yōu)化的應(yīng)用[M]. 北京: 科學(xué)出版社, 2011: 119.(編輯武曉英)

3結(jié)論

(1)等離子體刻蝕中Cl2刻蝕硅是比較簡單的一類,腔室壓強(qiáng)和射頻電壓不同程度地影響等離子體的流量和能量分布等特性.隨著射頻電壓的升高,離子的流量增大;在低射頻電壓時,離子流量隨腔室壓強(qiáng)的升高而減小,高射頻電壓下趨勢相反.

(2)總體來說,雖然優(yōu)化后的刻蝕形貌不是完全理想的直上直下的刻蝕形貌,但呈現(xiàn)各向異性很強(qiáng)、剖面光滑的形貌.將優(yōu)化后的參數(shù)組合再進(jìn)行刻蝕模擬,得到的目標(biāo)函數(shù)為65.5,并不是71.3,可能是由2個原因造成的:一是優(yōu)化模型選取不當(dāng),二是試驗點太少(個變量可以分為4個水平或更多).

通過對Cl2等離子體刻蝕硅的數(shù)值模擬與優(yōu)化,從科研應(yīng)用的角度提出一種等離子刻蝕機(jī)的優(yōu)化方法,可以為工藝條件相近的刻蝕機(jī)設(shè)備的研究和設(shè)計提供參考.參考文獻(xiàn):

[1]徐棟梁, 張鑒, 許曉琳. 微機(jī)械加工工藝的三維仿真模型及算法[J]. 電子科技, 2010, 23(8): 14.

XU Dongliang, ZHANG Jian, XU Xiaolin. Threedimensional models and algorithms for micromanufacture simulations[J]. Electronic Sci & Technol, 2010, 23(8): 14.

[2]王曉光, 朱曉明, 尹延昭. 等離子體刻蝕工藝的優(yōu)化研究[J]. 中國新技術(shù)新產(chǎn)品, 2010(14): 2224.

WANG Xiaoguang, ZHU Xiaoming, YIN Yanzhao. Optimization research of plasma etching process[J]. China New Technol and Products, 2010(14): 2224.

[3]朱炳金, 張國棟, 張向鋒. 電感耦合等離子體刻蝕InSb芯片工藝的研究[J]. 紅外技術(shù), 2009, 31(8): 467470.

ZHU Bingjin, ZHANG Guodong, ZHANG Xiangfeng. Inductive couple plasma etching processing of InSb wafer[J]. Infrared Technol, 2009, 31(8): 467470.

[4]陳曉南, 楊培林, 龐宣明, 等. 等離子體刻蝕工藝中工藝參數(shù)對刻蝕速率影響的研究[J]. 西安交通大學(xué)學(xué)報, 2004, 38(5): 546547.

CHEN Xiaonan, YANG Peilin, PANG Xuanming, et al. Influence of process parameters on the etching rate in inductively coupled plasma etcher[J]. J Xian Jiaotong Univ, 2004, 38(5): 546547.

[5]程嘉, 朱煜. ICP刻蝕機(jī)反應(yīng)腔室氣流仿真研究[J]. 半導(dǎo)體技術(shù), 2007, 32(1): 4346.

CHENG Jia, ZHU Yu. Gas flow simulation research on reaction chamber of ICP etcher[J]. Semiconductor Technol, 2007, 32(1): 4346.

[6]程嘉, 朱煜, 段廣洪, 等. 基于回歸正交設(shè)計的ICP刻蝕機(jī)工藝腔室流場特性分析[J]. 半導(dǎo)體學(xué)報, 2008, 29(4): 780784.

CHENG Jia, ZHU Yu, DUAN Guanghong, et al. Analysis of processing chamber flow field characteristics for an ICP etcher based on regression orthogonal design[J]. J Semiconductors, 2008, 29(4): 780784.

[7]SEO Seung T, LEE Yong H, LEE Kwang S, et al. Runtorun control of inductively coupled C2F6 plasma etching of SiO2: construction of a number process with a computational fluid dynamics code[EB/OL]. (20050705)[20130610]. http://2005.iccas.org/submission/paper/upload/ 12ICP_simulation.pdf.

[8]王凌, 吉利軍, 鄭大鐘. 基于代理模型和遺傳算法的仿真優(yōu)化研究[J]. 控制與決策, 2004, 19(6): 626630.

WANG Ling, JI Lijun, ZHENG Dazhong. Simulation optimization based on surrogate model and genetic algorithm[J]. Contr & Decision, 2004, 19(6): 626630.

[9]高月華, 王希誠. 基于Kriging代理模型的穩(wěn)健優(yōu)化設(shè)計[J]. 化工學(xué)報, 2010, 61(3): 676681.

GAO Yuehua, WANG Xicheng. Robust optimization based on Kriging surrogate model[J]. CIESC J, 2010, 61(3): 676681.

[10]謝延敏, 于滬平, 陳軍, 等. 基于Kriging模型的可靠度計算[J]. 上海交通大學(xué)學(xué)報, 2007, 41(2): 177193.

XIE Yanmin, YU Huping, CHEN Jun, et al. Reliability estimation based on Kriging model[J]. J Shanghai Jiaotong Univ, 2007, 41(2): 177193.

[11]隋允康, 宇慧平. 響應(yīng)面方法的改進(jìn)及其對工程優(yōu)化的應(yīng)用[M]. 北京: 科學(xué)出版社, 2011: 119.(編輯武曉英)

3結(jié)論

(1)等離子體刻蝕中Cl2刻蝕硅是比較簡單的一類,腔室壓強(qiáng)和射頻電壓不同程度地影響等離子體的流量和能量分布等特性.隨著射頻電壓的升高,離子的流量增大;在低射頻電壓時,離子流量隨腔室壓強(qiáng)的升高而減小,高射頻電壓下趨勢相反.

(2)總體來說,雖然優(yōu)化后的刻蝕形貌不是完全理想的直上直下的刻蝕形貌,但呈現(xiàn)各向異性很強(qiáng)、剖面光滑的形貌.將優(yōu)化后的參數(shù)組合再進(jìn)行刻蝕模擬,得到的目標(biāo)函數(shù)為65.5,并不是71.3,可能是由2個原因造成的:一是優(yōu)化模型選取不當(dāng),二是試驗點太少(個變量可以分為4個水平或更多).

通過對Cl2等離子體刻蝕硅的數(shù)值模擬與優(yōu)化,從科研應(yīng)用的角度提出一種等離子刻蝕機(jī)的優(yōu)化方法,可以為工藝條件相近的刻蝕機(jī)設(shè)備的研究和設(shè)計提供參考.參考文獻(xiàn):

[1]徐棟梁, 張鑒, 許曉琳. 微機(jī)械加工工藝的三維仿真模型及算法[J]. 電子科技, 2010, 23(8): 14.

XU Dongliang, ZHANG Jian, XU Xiaolin. Threedimensional models and algorithms for micromanufacture simulations[J]. Electronic Sci & Technol, 2010, 23(8): 14.

[2]王曉光, 朱曉明, 尹延昭. 等離子體刻蝕工藝的優(yōu)化研究[J]. 中國新技術(shù)新產(chǎn)品, 2010(14): 2224.

WANG Xiaoguang, ZHU Xiaoming, YIN Yanzhao. Optimization research of plasma etching process[J]. China New Technol and Products, 2010(14): 2224.

[3]朱炳金, 張國棟, 張向鋒. 電感耦合等離子體刻蝕InSb芯片工藝的研究[J]. 紅外技術(shù), 2009, 31(8): 467470.

ZHU Bingjin, ZHANG Guodong, ZHANG Xiangfeng. Inductive couple plasma etching processing of InSb wafer[J]. Infrared Technol, 2009, 31(8): 467470.

[4]陳曉南, 楊培林, 龐宣明, 等. 等離子體刻蝕工藝中工藝參數(shù)對刻蝕速率影響的研究[J]. 西安交通大學(xué)學(xué)報, 2004, 38(5): 546547.

CHEN Xiaonan, YANG Peilin, PANG Xuanming, et al. Influence of process parameters on the etching rate in inductively coupled plasma etcher[J]. J Xian Jiaotong Univ, 2004, 38(5): 546547.

[5]程嘉, 朱煜. ICP刻蝕機(jī)反應(yīng)腔室氣流仿真研究[J]. 半導(dǎo)體技術(shù), 2007, 32(1): 4346.

CHENG Jia, ZHU Yu. Gas flow simulation research on reaction chamber of ICP etcher[J]. Semiconductor Technol, 2007, 32(1): 4346.

[6]程嘉, 朱煜, 段廣洪, 等. 基于回歸正交設(shè)計的ICP刻蝕機(jī)工藝腔室流場特性分析[J]. 半導(dǎo)體學(xué)報, 2008, 29(4): 780784.

CHENG Jia, ZHU Yu, DUAN Guanghong, et al. Analysis of processing chamber flow field characteristics for an ICP etcher based on regression orthogonal design[J]. J Semiconductors, 2008, 29(4): 780784.

[7]SEO Seung T, LEE Yong H, LEE Kwang S, et al. Runtorun control of inductively coupled C2F6 plasma etching of SiO2: construction of a number process with a computational fluid dynamics code[EB/OL]. (20050705)[20130610]. http://2005.iccas.org/submission/paper/upload/ 12ICP_simulation.pdf.

[8]王凌, 吉利軍, 鄭大鐘. 基于代理模型和遺傳算法的仿真優(yōu)化研究[J]. 控制與決策, 2004, 19(6): 626630.

WANG Ling, JI Lijun, ZHENG Dazhong. Simulation optimization based on surrogate model and genetic algorithm[J]. Contr & Decision, 2004, 19(6): 626630.

[9]高月華, 王希誠. 基于Kriging代理模型的穩(wěn)健優(yōu)化設(shè)計[J]. 化工學(xué)報, 2010, 61(3): 676681.

GAO Yuehua, WANG Xicheng. Robust optimization based on Kriging surrogate model[J]. CIESC J, 2010, 61(3): 676681.

[10]謝延敏, 于滬平, 陳軍, 等. 基于Kriging模型的可靠度計算[J]. 上海交通大學(xué)學(xué)報, 2007, 41(2): 177193.

XIE Yanmin, YU Huping, CHEN Jun, et al. Reliability estimation based on Kriging model[J]. J Shanghai Jiaotong Univ, 2007, 41(2): 177193.

[11]隋允康, 宇慧平. 響應(yīng)面方法的改進(jìn)及其對工程優(yōu)化的應(yīng)用[M]. 北京: 科學(xué)出版社, 2011: 119.(編輯武曉英)

猜你喜歡
優(yōu)化
超限高層建筑結(jié)構(gòu)設(shè)計與優(yōu)化思考
PEMFC流道的多目標(biāo)優(yōu)化
能源工程(2022年1期)2022-03-29 01:06:28
民用建筑防煙排煙設(shè)計優(yōu)化探討
關(guān)于優(yōu)化消防安全告知承諾的一些思考
一道優(yōu)化題的幾何解法
由“形”啟“數(shù)”優(yōu)化運算——以2021年解析幾何高考題為例
圍繞“地、業(yè)、人”優(yōu)化產(chǎn)業(yè)扶貧
事業(yè)單位中固定資產(chǎn)會計處理的優(yōu)化
4K HDR性能大幅度優(yōu)化 JVC DLA-X8 18 BC
幾種常見的負(fù)載均衡算法的優(yōu)化
電子制作(2017年20期)2017-04-26 06:57:45
主站蜘蛛池模板: 狠狠干欧美| 久久青草热| 在线视频亚洲欧美| 亚洲第一国产综合| 99人妻碰碰碰久久久久禁片| 日韩人妻精品一区| 国产永久在线观看| 国产视频欧美| 日韩麻豆小视频| 少妇精品网站| 九色视频一区| 中文字幕第1页在线播| 国产精品制服| 亚洲欧洲自拍拍偷午夜色| 亚洲最黄视频| 白丝美女办公室高潮喷水视频| 日韩成人在线一区二区| 久久美女精品| 免费无码AV片在线观看中文| 精品国产三级在线观看| 在线观看91香蕉国产免费| 国产迷奸在线看| 91av成人日本不卡三区| 9999在线视频| 国产精品白浆无码流出在线看| 国产精品第一区| 国产一区二区三区在线无码| 欧美精品H在线播放| 91麻豆精品国产91久久久久| 欧美精品成人一区二区在线观看| 国内精品久久久久久久久久影视 | 欧洲高清无码在线| 精品色综合| 亚洲成网777777国产精品| 国产国语一级毛片在线视频| 国产成人亚洲精品无码电影| 国产精品性| 亚洲第一成年免费网站| 色偷偷综合网| 国产日韩欧美视频| 丝袜国产一区| 99在线观看视频免费| 99久久精品国产精品亚洲| av手机版在线播放| 91麻豆国产视频| 欧美成人午夜在线全部免费| 欧美在线精品怡红院| 亚洲一区色| 欧美在线精品怡红院| 色悠久久综合| 久久综合色播五月男人的天堂| 黄色网址免费在线| 国产精彩视频在线观看| 九九热精品在线视频| 91在线一9|永久视频在线| 日韩在线观看网站| 情侣午夜国产在线一区无码| 日韩av在线直播| 午夜老司机永久免费看片| 久久大香伊蕉在人线观看热2| 欧美日韩精品一区二区在线线| 国产呦精品一区二区三区网站| 国产人碰人摸人爱免费视频| 国产欧美日韩专区发布| 国产91视频免费观看| 天天做天天爱夜夜爽毛片毛片| 任我操在线视频| 亚洲精品午夜天堂网页| 婷婷午夜影院| 日韩在线成年视频人网站观看| 国产超薄肉色丝袜网站| 国产福利免费视频| 亚洲综合第一页| 成人国产精品一级毛片天堂 | 亚洲人精品亚洲人成在线| 久久久久人妻一区精品| 国产精品亚洲а∨天堂免下载| 欧美成人区| 国产精品欧美亚洲韩国日本不卡| 日本三级欧美三级| 国产黄色视频综合| 国产无码性爱一区二区三区|