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

熱循環與隨機振動下焊點可靠性的研究綜述

2022-07-10 09:51:12徐鵬博呂衛民劉陵順孫晨峰
兵器裝備工程學報 2022年6期
關鍵詞:振動模型

徐鵬博,呂衛民,劉陵順,孫晨峰

(海軍航空大學, 山東 煙臺 264001)

1 引言

如今電子信息技術日益提高有目共睹,微電子技術不斷地發展創新對其產生了巨大作用。同時,微電子技術也在受微電子封裝技術的影響,封裝技術的水平決定了電子產品的發展趨勢以及工作壽命。早在20世紀70年代之前,主流的微電子封裝技術是雙列直插技術(dual in-line package,DIP),但是這種封裝技術的芯片面積與封裝面積的比值很小,這反映出DIP技術封裝效率十分低并且大量的面積沒有得到有效利用。80年代出現了薄型小尺寸封裝技術(thin small outline package,TSOP),它的特點是將封裝芯片四周的引腳與印制電路板(printed circuit board,PCB)貼合。改進后的TSOP技術如今仍在被廣泛使用。在90年代,為了處理功耗急劇增大問題,適應集成技術的進步和工業生產需求,球柵陣列封裝技術(ball grid array package,BGA)應運而生,它的散熱性良好,質量與厚度相較之前的封裝技術降低了許多,同時使用BGA封裝技術可以延長產品使用周期。

根據美國空軍的相關研究顯示:統計應用在導彈上的電子元器件各類失效原因,40%的失效原因是電連接器部件出現了問題,29%的失效原因是電子元器件自身失效導致無法正常工作。焊接技術是指通過高溫手段將金屬填充材料融化,待其重新冷卻凝固后將金屬或者其他具有熱塑特性的材料連接在一起的技術。而連接這2種材料的一個個金屬點即為焊點。焊點作為媒介將封裝芯片與PCB板連接起來。焊點既承擔了電力連接的任務,使得電路元件和電子器件之間可以傳遞電信號;又承擔了機械連接的任務,可以起到應力緩沖的作用,將封裝芯片固定在PCB板上使其無法脫落;也提供了空間使封裝芯片形成散熱通路,在一定程度上降低產生的熱量對電子元件的影響。一旦焊點出現問題,勢必會造成電路板上的其他組件無法正常工作。焊點的失效原因和可靠性問題一直都是關注重點。在進行電子元器件可靠性研究時,必須要考慮到工作環境產生的影響。如圖1所示,電子元器件工作環境的調查研究表明:50%的故障是由溫度引起的,20%的故障是由振動沖擊導致的,其余導致故障的環境因素則是鹽堿濕度、灰塵、沖擊、低氣壓等等。所以重點研究溫度和振動對焊點可靠性的影響十分必要。

圖1 各種環境引起的失效率占比圖

2 熱循環研究現狀

電子元件在服役期間,會在工作狀態和停機狀態不斷切換,并且工作環境的溫度也會不斷變化,這種工作特性導致了電子元件會受到周期性熱循環加載。在對電子元器件的可靠性進行研究時,為了能夠使仿真結果更加真實具有參考意義,通常采用加速壽命試驗。將電子元器件放進高低溫循環箱中,溫度范圍設置在-55 ℃~150 ℃,一次熱循環時間通常設定為15~120 min。

2.1 熱循環下焊點失效原因

由于電子元器件內部各個部分的材質不盡相同,導致它們的熱膨脹系數(coefficient of thermal expansion,CTE)差異很大,例如常用芯片材料硅的熱膨脹系數為2.8×10/℃,常用的錫鉛焊料Sn63Pb37的熱膨脹系數為24.5×10/℃,基板中各項材料熱膨脹系數為217×10/℃~240×10/℃。對比該系數發現,材料之間存在失配情況,在受到周期性應力的過程中焊點內部會持續累積應變能,從而形成焊點從萌生裂紋到裂紋擴展最后失效的過程。封裝芯片和電路板也會因溫度升高出現熱膨脹和翹曲現象。而在焊接過程中,焊點與引腳之間的尺寸縫隙和焊點內部空洞也會使焊點受到的應力增加,加速焊點失效。在熱循環加載下,材料產生疲勞破壞所需要的熱應力循環數即為熱疲勞壽命。

2.2 熱循環下危險焊點的研究

對電子封裝結構焊點進行疲勞失效分析時,首先要找到最容易破壞的焊點及其裂紋開始產生的位置。這樣可以提高電子顯微儀器的檢測效率,檢測結果更加準確、直觀;在使用計算機軟件建模過程中,可以精確繪制危險焊點的三維模型并對其他位置的焊點進行適當的簡化,從而提高仿真速度;最重要的是上述學者的研究成果可以提高焊點的制造工藝水平,對延長焊點的疲勞壽命起到重要作用。

黃大巍對表面貼裝技術(surface mounted technology,SMT)的可靠性進行研究。分析了無引線陶瓷封裝載體(leadless ceramic chip carrier,LCCC)模型焊點在熱循環條件加載下的應力應變情況,總結了非彈性形變和蠕變隨熱循環的變化趨勢,有限元仿真結果和實驗確定了圓角邊緣處為危險位置。Song對細間距球柵陣列(Fine-pitch ball grid array,FBGA)焊點進行分析,建立了三維的全局模型和子模型。發現8行6列共48個焊點當中最外圈的焊點疲勞壽命最短,同時PCB板厚度增加會導致自身受到更大的應力,導致疲勞壽命減少。吳玉秀找出了方型扁平式封裝(plastic quad flat package,QFP)元器件中J形引線焊點最易發生破壞的位置,并對QFP翼形引線和引腳進行設計,根據數值模擬仿真和針對引腳的抗拉力試驗確定一種兼備可靠性與經濟性的三維結構。趙新新進一步討論了熱循環中各項參數與BGA焊點疲勞壽命的關系,控制單因素變量,將焊點高度、焊點半徑、PCB板厚度和尺寸與焊點疲勞壽命變化趨勢建立起數學關系。并將這4種尺寸因素與焊點所受應力、塑性形變、壽命等結合起來確定了可靠的二階回歸方程,為今后設計和優化電子封裝元件提供了數學模型。

2.3 焊點材料的選取與改進

隨著現代電子制造工藝的發展,焊料的選取標準一直在提高。將傳統的鉛錫焊點與新型焊點進行對比分析出各自的優缺點也是專家們關注的重點。

樊強對BGA上下兩基板的材料選擇進行研究,設計了多個材料搭配組合,得出上下基板均為陶瓷材質時消除了CTEC差異性,錫鉛焊點的疲勞壽命最大。隨后樊又使用控制變量法,在材料搭配、熱循環加載、焊點模型等條件相同情況下,對比無鉛焊點與傳統Sn60Pb40焊點的可靠性。結果表明無鉛焊點熱疲勞壽命高于Sn60Pb40焊點四倍,證明了無鉛焊點更適合應用在電子封裝技術上。李躍以PCB混合組裝為研究對象,在PCB基板上安裝BGA等5種常用電子器件,將再流焊和熱循環這2種方式結合起來,通過分析電子元件焊接在PCB基板上的殘余應力以及熱循環載荷下的應力應變計算出電子元件疲勞壽命。C/R元件焊點是5種元件中疲勞壽命最高的,混裝釬料作為焊點疲勞壽命也高于Sn62Pb36Ag2。這表明無鉛焊料具有更為廣闊的研究意義,正在逐步取代傳統的鉛錫材料;同時在選取焊料時也會考慮混裝配制各類合金材料從而提高焊點的可靠性。

2.4 熱循環加載條件研究

不論采取什么樣的方法預測焊點的熱疲勞壽命,選取符合裝備服役環境規律的熱循環加載方式是非常重要的。李長庚就針對PBGA有限元模型采用了五種不同的熱循環加載條件,發現高溫加載時間維持得越長塑性應變累積得越多,焊點更易破壞;高低溫溫差也是焊點熱疲勞壽命的一個重要影響因素。Ladani根據能量分配損傷原理設計出能更加準確預測粘塑性材料疲勞壽命的模型——EPDE。并根據該模型中彈性、塑性性能的改變采用了更具有參考性的熱循環加載。Ramesham對陶瓷柱柵陣列封裝(ceramic column grid array packaging,CCGA)電子元件采用-185~125 ℃,溫度升高和降低的速率為5 ℃/min,每次循環高溫和低溫保持時間均為15 min。在對電子元件加載熱循環1 258次后,發現63.2%的元件在經歷664次循環后菊花鏈已經遭受破壞。趙鑫比較Sn63Pb37和Sn3.0Ag0.5Cu兩種焊料作為陶瓷球柵陣列封裝(ceramic ball grid array package,CBGA)焊點時在常規熱循環-55~125 ℃和大溫差循環-150~150 ℃的差異性。有限元計算結果表明2種熱循環載荷下Sn63Pb37疲勞壽命都遠遠高于Sn3.0Ag0.5Cu,同時溫度極值過大會縮短焊點疲勞壽命。邵陳希研究了電子元件的傳熱特性對疲勞壽命的影響,在建立有限元模型施加熱循環載荷時考慮了氣流流動、元件散熱導致環境溫度升高、熱應力帶來的耦合效果等因素。發現除了環境溫度,電子元件的功率損耗也是影響元件熱疲勞壽命的重要因素。

以上研究成果表明:在熱循環加載條件下,溫度幅值和溫度變化速率的改變會對焊點的疲勞壽命產生很大影響。針對惡劣極端的服役環境,電子封裝產品更加要注重焊點的熱疲勞壽命,確保其在服役期間能保證較高的可靠性。

3 隨機振動研究現狀

在考慮振動沖擊對電子元器件的影響時,可以采用隨機振動這一理論對元件進行分析。事實上目前裝備工作環境越來越復雜化,用單一線性、理想化的方式研究振動問題是不準確的。在裝備工作的過程中,會有許多不確定的因素影響著裝備的使用壽命。例如:飛行器和氣流因為流固耦合作用產生振動;車輛在不平整路面顛簸前進時會受到振動載荷;艦船在海上隨著海浪起伏搖擺的同時受到海水拍打也會受到振動沖擊。這種無法用某一確定數學模型表達出來的振動形式,可以采用隨機振動方法進行處理。

3.1 隨機振動下焊點失效原因

除了溫度交變循環產生的熱應力對電子元器件的壽命造成影響,周期性的機械振動疲勞也是電子元器件失效的一個主要因素。電子產品為了滿足工業需求,一直在向微型化方向發展,電子組件的尺寸、組裝間距越來越小。因此連接基板和PCB板的焊點尺寸在減少,數量在提高。電子元器件在工作期間總是會不可避免受到振動沖擊。振動沖擊帶來的應力會使PCB板和電子元件的基板發生動態的彎曲變形,隨著彎曲撓度增大,二者之間的焊點必然會受到應力并產生諸如塑性變形等一系列影響,加劇裂紋的萌生和擴展。焊點結構上的破壞將導致電信號不能傳輸,整個電子產品徹底失效。

3.2 隨機振動下危險焊點的研究

Jih根據斷裂力學相關理論,分別建立了整體模型和包含電子元件引腳、焊點的局部精確模型。研究了隨機振動載荷對焊點裂紋擴展的影響。王紅芳以SMT焊點作為研究對象,將焊點形狀和焊點工作環境作為試驗變量。根據試驗結果發現,裂紋的出現是從焊點底部的圓形區域開始,但是焊點頂部的裂紋擴展速率要遠遠高于其他位置。最終焊點失效斷裂的位置是在焊點頂部,而最早出現裂紋的底部擴展得并不顯著。嚴煥斌以電子混裝組件為對象進行隨機振動分析,根據有限元分析結果找到組件薄弱易損壞的位置和最大翹曲位置。同時嚴在分析危險焊點時選擇子模型法,在焊點的邊界條件處理上選擇用等效載荷來替代實際受到的應力載荷。在保證較高精度的情況下對BGA、翼型及J型引腳焊點進行比較分析,找到了各自的危險位置并預測出3種焊點的振動疲勞壽命。李春洋利用有限元方法對PCB板進行模態分析,得到PCB板的前五階固有頻率和振型。通過模態振型找出PCB板上元件容易損壞的位置,并在PCB板布局、PCB板固定形式、PCB板厚度、PCB板材料屬性上做出優化調整,提高PCB板的固有頻率,從而達到避免共振現象提高可靠性的目的。張龍利用Patran仿真軟件對疊層焊點進行頻率響應分析,在利用模型計算焊點疲勞壽命時首先將有限元仿真得出的隨機振動下焊點應力應變結果從頻域轉換到時域上,在處理時域信息時利用雨流計數法得出應變分布函數。雨流計數法能夠直觀地呈現出電子元件經歷的載荷歷程,使得計算危險焊點的疲勞壽命結果精準高效。黃春躍以光互連模塊作為研究對象,在隨機振動條件下以發射激光器VCSEL和耦合元件在水平、軸向、垂直這3個方向的偏移量作為數據參考。分別對焊點高度、體積等4種因素選取4種水平進行16組正交試驗,對這些三維模型得出的有限元結果進行方差分析,發現陶瓷基板焊點的高度對偏移量的影響最大。

3.3 功率譜密度加載方式研究

功率譜密度(power spectral density,PSD)作為隨機振動理論中的重要部分,是一個用作表示信號的功率能量與頻率之間關系的物理量。改變功率譜密度加載條件會對電子元件的疲勞壽命產生很大影響。王文在隨機振動試驗中采取對不同位置進行激勵的方法并比較試驗結果,同時考慮了試驗過程中邊界條件預緊力設置對結果的影響。作者在對SMT焊點進行壽命預測時采用2種方法:第一種方法是利用Basquin方程,以焊點上受到的最大拉應力作為參量,建立壽命預測模型;第二種方法根據統計學模型得出不同加速度PSD與焊點壽命的分布關系,建立了以加速度PSD幅值為參量的預測模型。劉文杰根據有限元法和Dirlik模型對如何篩選隨機振動剖面做出了分析。以加速度均方根、激勵頻率范圍、激勵方向3個參數作為變量,研究3種參數的變化給電子元件的疲勞壽命帶來的影響,結果表明:為了節約時間提高試驗效率可以提高加速度均方根,降低疲勞壽命;選取頻率范圍時需要考慮構件一階固有振動頻率;相比于垂直電路板方向施加激勵造成的影響,其余兩個方向所施加激勵可以忽略不計,在試驗時僅僅考慮垂直方向即可。

3.4 焊點的疲勞壽命預測模型

不論是熱循環加載還是隨機振動加載,焊點的損壞都是因為周期性應力不停地作用在焊點上致使其疲勞失效的。在分析焊點的疲勞壽命時,會采用實驗法和有限元仿真的方法。實驗法有云紋干涉儀、X射線、貼裝電阻應變片等方法,此類方法雖然直觀有效,能夠精準的計算焊點疲勞壽命,但也有著測試成本高,測試時間長,部分測試技術具有針對性不能適用于多數電子產品等缺點。隨著電子封裝技術向著元件尺寸更加微小,元件之間空隙逐漸縮小的方向發展,利用實驗儀器觀測電子元件預測疲勞壽命變得愈發困難。現如今利用有限元仿真軟件和經驗公式預測得出的結果更加接近真實壽命,根據仿真軟件得出的力學參量類型選擇合適的壽命預測模型。在熱疲勞壽命預測方面,Manson等總結了應力和疲勞壽命之間的關系以及裂紋產生的規律,分析在室溫和高溫情況下低周期和高周期循環產生的疲勞損傷結果,先后提出了Manson-Coffin模型、Ostergren模型。在機械振動疲勞壽命預測方面,Steinberg等則基于電子元件相對于PCB板是剛性部件、隨機振動載荷下的應力應變峰值均服從Rayleigh分布等假設提出了各自的疲勞壽命預測模型。如圖2所示,經過大量的研究和總結,壽命預測模型一般可分為4種:

圖2 經典疲勞壽命預測模型框圖

1) 塑性變形疲勞預測模型:該模型是以塑性變形量這一力學參量作為研究重點,塑性變形具有與時間無正相關關系的特性,將剪切變量和相關系數代入相對應的公式中得到元件最終失效所經歷的周期循環次數。典型的模型有Manson-Coffin模型、Ostergren模型、Miller模型、Engelmaier模型等。

2) 蠕變變形疲勞預測模型:當電子元件和焊點在經歷低周期高低溫交變循環過程時,會出現蠕變這一特性。蠕變的產生和累積過程十分復雜,并且此類模型不考慮塑性變形也無法完整表達出整個蠕變過程,應用并不廣泛。典型的模型有:Syed模型和Kencht-Fox模型。

3) 斷裂參量疲勞預測模型:根據斷裂力學原理,將斷裂參量作為研究對象進行壽命預測。在電子元件和焊點經歷循環載荷時,將裂紋從產生、逐漸累積最終到破壞的過程與疲勞壽命建立聯系。此類模型考慮到了材料不同帶來的影響,但需要對裂紋產生和發展趨勢做出預先假設同時無法預測出裂紋產生前的壽命,導致預測結果小于真實的疲勞壽命值。典型的模型有:Paris模型和J積分模型。

4) 能量疲勞預測模型:該模型是以電子元件和焊點內部應力應變引發的遲滯效應作為研究基礎,溫度載荷每一次循環都會產生相對應的應變能。當應變能累積到一定程度就會引發元件斷裂致使其失效。這一模型利用能量與疲勞推導出壽命經驗公式。典型的模型有:Akay模型和Darveaux模型。

4 熱振耦合研究現狀

隨著國內外學者對焊點可靠性的深入研究,對于單一載荷產生的影響已經有了很深刻的認識。事實上,電子產品不可能只在理想環境下受到一種載荷作用,為了更加準確地分析其可靠性,就必須模擬設計出最接近裝備真實服役工作的環境。而溫度和振動正是影響電子元件壽命的兩大最主要因素,所以要將熱循環加載和隨機振動加載放在一起進行研究。

事實上,正如前文所介紹的熱循環載荷和隨機振動載荷下的疲勞壽命研究已經十分成熟。無論是選取公式模型、試驗設計和后續驗證等環節,還是在深度學習、遺傳算法等領域都有很大進展。但是熱振耦合作用下的焊點可靠性研究還是不夠深入,載荷加載順序和疲勞壽命的計算方法值得深入研究。

4.1 線性損傷疊加法

計算熱振耦合條件下焊點的疲勞壽命時,在工程上應用較為廣泛的是線性損傷疊加原理,該理論認為各類載荷施加下造成的構件疲勞損傷都是獨立存在的,互不影響。這些獨立的疲勞損傷線性累加在一起構成了總損傷疲勞值。一旦總損傷疲勞值達到了預定值時,就會導致構件損壞。

Barker在研討熱循環和振動加載的綜合效應問題時,認為不能將材料僅僅看作簡單地彈性應變或是非彈性應變,需要在材料分析上考慮總應變原理。同時結合廣義上的Manson-Coffin模型和Miner準則來預測材料的疲勞壽命。Ping Yang認為以往的研究是將熱應力作為關鍵載荷,將熱循環和隨機振動建立起聯系,但是這樣忽略了蠕變損傷的作用,從而造成不可忽視的誤差。Yang利用有限元軟件對線性疊加損傷法進行改進,結果表明疲勞損傷程度與振幅幅值成正比,而蠕變損傷在逐漸趨于穩定時會經歷一個后減小的過程。Anthony Konoza在熱振載荷加載條件下對焊點材料進行分析。在振動方面,錫鉛材料可靠性要更優于無鉛合金材料SAC;在熱循環載荷方面,SAC材料明顯更加具有優勢。但是在組合加載條件下,多組數據顯示錫鉛材料可靠性和使用價值更高,這表明振動載荷在組合載荷中占據更主要的地位,振動載荷會累積更多的損傷。

4.2 遞增損傷疊加法

遞增損傷疊加法是K.Upadhyayula和A.Dasgupta在熱振耦合問題上提出一種新的壽命預測方法,在添加振動加載前對材料施以平均應力,以此來代替熱循環帶來的載荷作用。該方法認為:各類載荷施加在構件上時并非單獨作用,耦合作用下各類疲勞損傷值會發生改變。但K.Upadhyayula這種處理方式僅僅是考慮了溫度對振動的作用,忽略振動對溫度的影響。2種因素之間只有單一影響卻沒有相互作用,這在邏輯上并沒有形成閉環。

Haiyu Qi在遞增損傷疊加方法繼續深入研究,結果表明在將熱循環和隨機振動結合起來對PBGA焊點進行加載,其損傷程度要遠高于單獨加載所造成的。這一結論與K.Upadhyayula持有不同的意見,Haiyu Qi通過對損傷公式的檢查并分3種情況討論證明了自身觀點的正確性。Haiyu Qi又在此前研究的基礎上提出了一種預測疲勞壽命的快速方法,將熱應力和隨機振動產生的機械應力通過應變計算結合在一起。隨后將該方法應用在PBGA上進行驗證。Tilman Ecker對SAC305焊點材料進行組合載荷試驗,根據試驗數據和環境工況進行建模,并比較了線性損傷疊加法和遞增損傷疊加法。在2種方法沒有明顯優勢的情況下,Ecker認為用線性損傷疊加計算疲勞壽命更加簡單,更適合使用。

3.3 構建耦合作用下的疲勞壽命模型

上述總結的2種方法本質上還是將熱循環和振動載荷造成的疲勞損傷進行疊加處理。而構建耦合作用下的疲勞壽命模型這種方法是對構件在熱循環和隨機振動加載下的試驗數據進行分析和處理,建立出焊點在2種不同載荷耦合作用下的疲勞壽命數學模型。

Basaran Cemal分析和測量了熱循環和振動試驗下的焊點變形,由于尺寸結構不同所引發的尺寸效應,表明已知的線性疊加方法并不完全適用于各種材料。同時,Cemal創立了以本構模型為基礎的耦合分析方法,他將該方法與Miner法則進行比對,發現針對組合載荷Miner法則所得到疲勞壽命要高于材料的真實值。Andrew Eugene Perkins針對電子元件面對復雜的溫度、振動以及動力循環環境,將諸多不同類型的試驗數據采集起來,依靠多值預測方程和統一建模提出了非線性累積法。湯巍利用正交試驗對多組溫度和振動的耦合情況進行分析,并根據多項式擬合的方法建立預測模型,利用正交試驗法、單一時間因子傳遞熵建立了焊點的疲勞壽命預測模型和非經驗壽命預測模型。路良坤在對芯片尺寸封裝(chip scale package,CSP)焊點模型進行優化設計時,以隨機振動下的應力應變和電磁仿真得出的回波損耗作為參考。對于這種雙目標進行優化時作者采用了響應面法,先進行各結構參數的組合試驗,利用泰勒展開式建立了雙目標數學模型。最終通過遺傳算法對模型參數進一步優化,設計出模型最優參數組合。

5 結論

國內外研究學者對于焊點的失效機理仍在不斷建立和改進預測疲勞壽命的模型和理論方法。但研究成果絕大部分還是以單一載荷環境為前提。針對同時進行2種及2種以上載荷的疲勞失效研究比較薄弱。焊點在熱循環和隨機振動加載下的疲勞壽命僅僅根據損傷疊加法計算還不精確,載荷加載順序和元件疲勞壽命的計算方法值得深入研究。在2種類型的損傷疊加、耦合作用下的疲勞壽命模型應用并不廣泛仍需更多驗證的情況下,預測焊點的疲勞壽命應該更多地將有限元仿真方法與實驗法結合。因此,有限元仿真方法要提高計算機計算能力,建立更加真實的模型、減少簡化過程;實驗法也需要提高儀器的精度,使測試結果更加精準。這2種方法作為焊點可靠性研究中的重要方法,應當共同發展,不斷提高電子元件在實際工程應用中的可靠性。

猜你喜歡
振動模型
一半模型
振動的思考
科學大眾(2023年17期)2023-10-26 07:39:14
噴水推進高速艇尾部振動響應分析
重要模型『一線三等角』
This “Singing Highway”plays music
重尾非線性自回歸模型自加權M-估計的漸近分布
振動攪拌 震動創新
中國公路(2017年18期)2018-01-23 03:00:38
中立型Emden-Fowler微分方程的振動性
3D打印中的模型分割與打包
FLUKA幾何模型到CAD幾何模型轉換方法初步研究
主站蜘蛛池模板: 国产午夜在线观看视频| 国产精品一区二区久久精品无码| 日本五区在线不卡精品| 99热这里只有精品在线播放| 老色鬼久久亚洲AV综合| 国产免费a级片| AV片亚洲国产男人的天堂| 欧美成人A视频| 一区二区午夜| 色偷偷综合网| 操操操综合网| 无码一区二区三区视频在线播放| 精品天海翼一区二区| 国产人人射| 欧美专区在线观看| 无码中文AⅤ在线观看| aaa国产一级毛片| 久久香蕉国产线看观看式| 尤物视频一区| 国产真实乱了在线播放| 99成人在线观看| 国产精品毛片一区| 亚洲区视频在线观看| 秋霞午夜国产精品成人片| 亚洲欧美激情小说另类| 91九色最新地址| 国产真实二区一区在线亚洲| 精品偷拍一区二区| 又黄又湿又爽的视频| 欧美日在线观看| vvvv98国产成人综合青青| 免费无码AV片在线观看国产| 国产又粗又爽视频| 国产高清毛片| 麻豆AV网站免费进入| 中文无码伦av中文字幕| 国产第一页屁屁影院| 日韩高清欧美| 999国内精品视频免费| 亚洲一区二区三区在线视频| 国产成人高清亚洲一区久久| 国产成人高清在线精品| 91亚洲视频下载| 国产精品99久久久| 国产成人av一区二区三区| 国产精品久久自在自线观看| 国产99视频免费精品是看6| 五月婷婷综合在线视频| 男人的天堂久久精品激情| 欧美午夜网| 欧美日韩精品综合在线一区| 亚洲一级毛片在线观| 亚洲天堂区| 99色亚洲国产精品11p| 亚洲一级毛片免费观看| 久久综合伊人 六十路| 在线播放国产一区| 国产精品成人观看视频国产| 再看日本中文字幕在线观看| 欧美日韩国产精品va| 国产免费久久精品99re不卡 | 日韩欧美中文字幕一本| 日本精品αv中文字幕| 5388国产亚洲欧美在线观看| 免费在线不卡视频| 日本不卡在线播放| 国产精品欧美在线观看| 日韩欧美中文在线| 人妻丰满熟妇AV无码区| 成人国产精品2021| 欧美日韩精品一区二区视频| 在线免费a视频| 国产激情无码一区二区免费| а∨天堂一区中文字幕| 性色一区| 天天爽免费视频| 日韩麻豆小视频| 中文字幕在线观| 国产成人精品一区二区秒拍1o| 亚洲欧美精品一中文字幕| 91欧美在线| 国产免费黄|