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

航天器電子產品抗隨機振動環境設計方法研究

2016-09-08 03:20:30劉天雄姜萬杰范燕平北京空間飛行器總體設計部北京100094
航天器工程 2016年3期
關鍵詞:振動產品

劉 晨 劉天雄 姜萬杰 范燕平(北京空間飛行器總體設計部,北京 100094)

航天器電子產品抗隨機振動環境設計方法研究

劉晨劉天雄姜萬杰范燕平
(北京空間飛行器總體設計部,北京 100094)

航天器電子產品在航天器發射過程中經歷惡劣、復雜的隨機振動環境,可能會引起電子產品失效。根據電子產品內部印制電路板(PCB)隨機振動原理和特點,文章研究梳理了航天器電子產品隨機振動失效模式及抗隨機振動設計的5條一般原則,總結出航天器電子產品進行抗隨機振動能力分析評估的兩種方法和風險判據,利用有限元仿真分析證明了上述分析方法和判據正確可行。針對失效模式,從兩個方面提出了5種降低產品隨機振動動力學響應的可行方案,通過計算仿真及試驗表明:產品內部PCB位移響應和變形極限大幅降低,減振方案有效。以上分析方法、判據及減振方案,可以作為同類電子產品抗隨機振動能力評估的參考,為隨機振動失效問題提供解決的思路和途徑。

航天器電子產品;隨機振動;分析方法;抗力學環境設計;減振方案

1 引言

航天器在發射過程中必須經歷運載火箭起飛引起的聲環境和運載器發動機工作時燃燒不穩定引起的振動環境激勵,此類聲環境呈現一定的隨機性,頻帶很寬,其能量通常集中在50~10 000 Hz[1]。隨機振動環境引起產品內部印制電路板組件(PCBA)產生劇烈的動力學響應,由此產生的變形和應力可導致模塊或電子器件失效、線路短路或斷路、焊點松動、開裂等問題,對航天器電子產品的穩定性和可靠性影響較大。根據美國空軍航空電子系統完整性項目(Avionics Integrity Program)的統計數據,振動環境因素所引起的電子產品失效所占比重已經接近30%[2-3]。

目前,國內航天器電子產品開展的抗隨機振動環境設計存在兩方面薄弱環節:①缺乏系統層面的統籌考慮,如產品安裝源頭減振、機箱內部PCBA結構優化、多塊PCB基板布局及集成模塊(如電源模塊及其他功能模塊)布局等;②缺少內部細節的計算仿真,如PCBA動力學響應、元器件引腳的焊接應力、元器件在PCBA的布局位置和安裝方式的優化等。如果產品未通過力學環境試驗考核,則產品需要重新設計并再次開展試驗驗證,造成時間和成本的巨大浪費。

針對上述問題,本文對電子產品抗隨機振動環境設計工作進行了系統的研究,根據隨機振動環境特點,梳理出航天器電子產品典型失效模式,給出2種可行的、可量化的電子產品PCBA動力學響應分析方法和失效風險判據。工程實例仿真分析和試驗數據表明:方法和判據正確,可提高產品設計人員抗隨機振動設計能力;此外,從2個方面探討了5種針對產品內部動力學響應過大問題的減振方案,試驗表明上述方案可行、有效,可為解決類似問題提供思路和途徑。

2 航天器電子產品隨機振動失效模式及設計原則

2.1常見失效模式

航天器電子產品隨機振動環境下常見失效模式如下:

(1)產品機箱結構、PCBA結構設計不合理,使PCBA與箱體頻率耦合,PCBA共振導致劇烈的動力學響應,使PCB基板上敏感元器件響應超出其抗振能力而失效;

(2)元器件與PCBA變形不匹配,導致引腳及連接處焊點局部應力水平超出焊接材料屈服極限,使引腳焊點開裂或斷裂而失效;

(3)PCBA布局不合理,兩個甚至多個板間距不足,在振動環境下振幅較大,導致板間發生接觸,板上元器件由于接觸碰撞引起短路失效;

(4)電連接器接口處未采取防脫落措施,導致在振動劇烈時脫插或脫落。

PCBA與其上器件連接處應力和敏感器件的振動響應超限,是造成電子產品隨機振動環境下失效的主要原因,因此,應充分考慮PCBA基頻與結構頻率錯開,以降低PCBA的振幅及響應,或采取減振措施從源頭降低PCBA振動幅值。

2.2設計原則

在抗隨機振動設計過程中,從PCBA的設計到產品的抗隨機振動環境設計,一般應遵循以下幾個原則:

(1)避免PCBA與機箱結構之間發生頻率耦合:在隨機振動環境下,機箱結構和PCBA的固有頻率都會受到激勵,在滿足剛度要求的前提下,應盡可能使內部PCBA基頻與機箱結構錯開,錯開頻率一般應滿足兩倍頻要求。

(2)使用減振器減小PCBA的位移和應力:為了減少與動態位移相關的PCBA失效問題,使用減振器可以減小PCBA的彎曲位移,進而減小焊接應力等。減振器安裝位置優選PCBA中心位置,軟橡膠減振器多用于固有頻率低于50Hz的PCBA上(應對在軌空間環境對橡膠材料的影響進行預先評估);而固有頻率高于100Hz的PCBA建議使用玻璃纖維環氧減振器。

(3)一般元器件應均勻布局,抗振性能較差、較重器件應靠近緊固點布局,避免由于質量效應引起大的振動響應,必要時可布局在機箱側壁或箱底位置。

(4)PCBA形狀應力求簡單,一般為長寬比例不大的長方形,慎用異型電路板并優化緊固點數量。

(5)對于產品內部有振動敏感的器件、模塊,可以在產品安裝面采取減振措施,如橡膠墊、金屬橡膠等。

3 航天器電子產品抗隨機振動環境設計分析方法

開展航天器電子產品抗隨機振動環境分析時,需要重點對PCBA上的振動響應、應力應變等數據進行計算,量化評估抗隨機振動能力。航天電子產品內部PCBA與箱體一般采取四周支撐或約束的連接方式,可以將它們簡化為單自由度系統,如圖1所示,m為單自由度系統質量塊,k為彈簧剛度,c為阻尼。

圖1 單自由度系統對隨機振動輸入響應的位移曲線Fig.1 Random vibration displacement response curve for single DOF system

通過積分隨機振動響應功率譜密度函數可以獲得單自由度振動系統響應均方根加速度水平,即

式中:GRMS=Gout為響應均方根加速度;Pout為響應功率譜密度;f為振動頻率;f1,f2為頻率范圍。

響應功率譜密度可以表示為

式中:Q為系統傳輸率;P為輸入功率譜密度[4-7]。

考慮一個諧振力P0cos(Ωt)作用在單自由度阻尼系統上,則根據振動微分方程有

式中:P0cos(Ωt)為諧振力函數;Ω為角頻率,將特解Y=Y0cos(Ωt—θ)帶入式(3),θ為初始相位角(rad),則最大位移表示為

將式(6)代入式(1)(2),有

若近似考慮弱阻尼系統頻率諧振處的傳輸率,即f/fn=1,其中fn為系統諧振頻率,則式(6)可近似為c/cc=1/2Q的關系;由積分公式(7)得到PCBA隨機振動均方根加速度響應:

將c/cc=1/2Q代入(8)得到對于寬帶隨機振動輸入的均方根響應:

根據頻率和加速度與位移的關系,單振幅位移可表示為[7]

隨機振動條件下,單自由度系統位移可用總均方根響應來評估,并使用3σ的加速度水平,得到:

式中:P為PCBA基頻處的輸入加速度功率譜密度,gn2/Hz;fn為PCBA的諧振頻率,Hz;Q為PCBA諧振頻率下的傳遞率,與諧振頻率存在如下關系(a取0.5~2的系數,在系統阻尼未知的情況下近似估算可取1.0[7-11])

(1)分析方法1:當我們需要評估PCBA隨機振動載荷下最大位移響應時,可通過獲得PCBA安裝點在PCBA基頻處的功率譜密度值P,將P帶入式(11)來得到預估位移響應。

(2)分析方法2:Lauetal和Dave.S.Steinberg曾做了大量研究和試驗,在積累大量試驗數據基礎上,提出了PCBA焊接引腳疲勞斷裂的失效判據[4,7](文獻單位為英寸):

當PCBA最大振幅

則PCBA焊接引腳振動疲勞壽命會低于2×107次振動周期,使PCBA早期斷裂失效。式(12)中:B為平行于元器件的PCB邊緣長度;L為元器件長度;h為PCB厚度;C為不同類型元器件安裝方式因子(標準雙列直插式封裝取1.0);r為元器件在PCB上相對位置因子(PCBA中心處取1.0)[7],其他安裝方式參數可參考文獻[7]。

(3)分析方法3:航天業內參考歐洲Alcatel ETCA公司標準,在承受力學載荷時,PCBA的彎曲變形可能會造成元器件及其引腳的損傷與破壞,彎曲變形的量度以基板的最大撓度除以基板的寬度,正則化為單位長度的變形來表達,由此規定的彎曲變形極限為0.006。

因此,在進行電子產品抗隨機振動分析時,可以針對敏感器件建立具體、詳細的有限元模型,計算敏感器件在PCB基板安裝處的隨機振動響應及焊點應力水平等數據,根據上述3個分析方法評估產品PCB基板振動幅值,評估其上元器件的抗振能力,從設計源頭避免出現諧振、動態位移及響應過大的問題。

4 某航天器電子產品隨機振動失效分析

某航天器電子產品采用鋁制箱式機箱,外形尺寸280mm×285mm×220mm,內部有9塊PCBA垂直X軸插接在底面板上,總質量12.6kg。產品側壁內表面安裝有3塊定制電源模塊,底板及側壁材料為鋁合金2A12,壁厚為2.5mm,內表面有截面為2.5mm× 3mm的加強筋,外表面平整,如圖2所示。

圖2 某航天器電子產品構型及電源模塊布局位置Fig.2 Electronic device structure and DC/DC position

產品在進行X向鑒定級隨機振動試驗過程中(100~2 000Hz,總均方根為12.6gn),發現產品電源模塊1#無輸出。該電源模塊為定制產品,上下分布兩塊PCBA(厚度均為1mm,采用通用玻璃環氧覆銅箔板FR-4),其中上PCBA尺寸80mm×90mm,通過4個鋁制桿件和M3螺釘安裝在電源模塊底板上,如圖3所示。試驗中止產品開蓋檢查,發現上板表面的過流保護三極管(三極管長度L=5mm)完全從焊盤上脫落且粘附在上蓋板導熱墊相鄰的導熱硅脂上。拆解模塊3#(經歷相同的力學環境試驗的同型號電源模塊)發現,在相同位置的同型號三極管引腳焊點部位也發現較大裂紋,存在失效隱患。

圖3 電源模塊內部結構圖及器件焊點斷裂脫落Fig.3 Inner structure of 1#DC/DC and elements broken off

4.1 建立電子產品有限元模型

某航天器電子產品由機箱殼體(含蓋板、內部支撐件等)、內部電路板組件、獨立安裝的模塊、元器件等部分組成,采用MSC.PATRAN商用FEM軟件建立電子產品有限元模型:機箱的側壁、底板、蓋板等結構采用板單元模擬,加強筋采用偏置梁單元模擬;對于相對產品面積、質量大的組件,采用剛性單元(RBE2)、梁單元模擬器件間的連接;印制板電路上分布許多小元器件,則以均布質量分布在印制板上;膠接、螺接、鉚接和焊接一般理想化為剛性連接或根據情況建立梁單元(模型單位為kg·mm·s),材料參數表見表1。

3塊電源模塊安裝在產品側板上,將電源模塊內部簡化為厚度h=1mm的PCBA(材料為玻璃環氧覆銅箔板FR4)并用板單元模擬,失效元器件通過4梁單元模擬引腳(含焊錫)連接在PCBA上,如圖4所示,模塊上其他附件質量配于側板安裝面上,模塊通過梁單元模擬緊固件安裝在產品側板上。

表1 材料參數表Table1 Parameters of material used in model

圖4 某航天器電子產品及電源模塊簡化模型Fig.4 Simplified model of electronic device and DC/DC

根據低量級正弦掃描試驗實測響應(測點位置如圖5所示),S7測點試驗曲線如圖6(a)所示,其中機箱側板一階基頻為295Hz,修正模型后機箱側板縱向一階基頻285Hz,仿真分析曲線與試驗曲線頻率、峰值、譜型基本一致,如圖6(b)所示,數據對比如表2所示。由此可知產品有限元模型較好地模擬了實際狀態,可以作為計算響應分析工作的模型。

圖5 模塊內上PCBA板測點及產品測點位置圖Fig.5 Measurement points location in PCBA and whole module

圖6 航天器電子產品S7測點試驗及分析響應曲線Fig.6 0.2gnsine vibration test response curve and simulation response curve at S7

表2 航天器電子產品正弦掃頻試驗測點數據與仿真結果主要頻點對比Table2 Main frequency data comparison between experiment and simulation

4.2隨機振動響應分析

產品按照鑒定級隨機振動條件進行加載,電源模塊1#的PCBA板在鑒定級隨機振動條件加載后,輸出PCBA緊固點及板上典型位置的響應曲線(如圖7(a)所示)、器件焊點應力(如圖7(b)所示)及結果數據(如表3所示)。

圖7 PCBA板響應輸出點響應曲線及Von Mises應力圖Fig.7 PCBA response point response curve and Von Mises stress

表3 PCBA隨機振動響應結果Table3 PCBA response results of random vibration

按照分析方法1,PCBA中心點振動響應均方根達到257.5gn,在PCBA基頻處輸入功率譜密度取P=150gn2/Hz,fn=303 Hz,系數a取1.0,代入公式(11),獲得估算的PCBA中心位移Y=9 mm。

按照分析方法2,將B=90mm、h=1mm、L= 5mm換算英寸單位帶入式(12)得到DC/DC模塊PCBA的最大位移Amax=1.2mm,即當PCBA中心位移大于1.2mm則焊接引腳振動疲勞壽命會低于2×107次振動周期。由方法1估算PCBA中心位移9mm超過該狀態的要求,因此PCBA存在早期斷裂失效的風險。

按照分析方法3,彎曲變形的量度以板的最大撓度除以板的寬度,則彎曲變形極限=9mm/90mm=0.1,超過了0.006的歐洲Alcatel ETCA公司的標準要求。

應力分析表明:PCBA緊固點錫焊連接處最大應力達到53.1MPa,如圖7(b)所示,超過錫焊接點材料(Sn3.5 Ag/Cu6Sn5)屈服極限35MPa,使得焊點處開裂破壞,與圖3中實際破壞形式是相符的,當焊點處有缺口損傷時會產生應力集中,其實際值會遠大于計算值。

因此,電源模塊失效的原因為隨機振動過程中模塊內PCBA板與安裝側板頻率耦合(機箱側板縱向基頻295Hz,PCBA縱向基頻303Hz),在隨機振動過程中,PCB基板共振產生劇烈彎曲變形,超出彎曲變形極限,三極管引腳焊點應力超出焊料許用應力導致開裂,使板上元器件發生引腳斷裂脫落。

4.3減振方案及驗證

根據第2節電子產品的失效模式,可從提高PCBA剛度來實現PCBA與機箱結構解耦和增加減振阻尼措施兩大方面,來降低PCB基板的位移幅值和加速度響應,具體給出5種減振方案。

(1)提高局部剛度實現解耦:方案1將失效元器件所在PCB板厚度h由1mm增加到1.6mm來提高PCBA剛度;方案2在失效元器件所在PCB板中心位置增加一個多點約束(MPC)剛性單元與產品側板連接來提高PCBA剛度;方案3在機箱側板外增加截面尺寸為5mm×11mm加強筋來提高機箱側板剛度。

(2)增加減振阻尼措施降低響應:方案4在模塊4個安裝腳引入阻尼器單元(阻尼0.2,垂直PCBA方向)來降低PCBA響應;方案5在產品安裝耳片處引入阻尼器單元(阻尼0.2,沿隨機載荷方向)來降低整機響應。

按照上述減振方案修改4.1節模型參數計算產品內部隨機振動響應,各方案輸出響應曲線如圖8所示,響應數據如表4所示。方案1和2通過不同方法改變了PCBA的基頻,方案3加強機箱側板的剛度實現單板與產品結構基頻錯開,在限制板上響應水平和位移方面效果明顯;方案4利用PCBA安裝處引入減振阻尼降低PCBA上的響應水平;方案5通過整機產品減振阻尼達到降低PCBA響應的目的。

圖8 方案1~5隨機振動響應曲線與原狀態對比Fig.8 State 1~5 response curves compared with that of original state

由5種方案分析估算中心位移均小于1.2 mm,滿足疲勞壽命的要求;變形極限除了方案5不滿足0.006的歐洲Alcatel ETCA公司的標準要求外,其他措施均滿足要求,證明5種減振措施均在降低PCBA振動響應和振幅方面有效,如表4所示。

表4 各方案隨機振動分析結果Table4 PCBA response results of random vibration for various states

綜上幾個方案,考慮空間電子產品安裝面散熱的需求,在盡量減少電源模塊內部改動的條件下,優選了方案1。按照方案1更改狀態后的產品,順利通過鑒定級隨機振動試驗,證明上述減振方案有效。

5 結束語

本文對航天器電子產品隨機振動失效模式及抗隨機振動設計分析方法進行了研究,給出了一種對PCBA元器件振動失效進行定量分析的方法和判據,通過工程案例仿真分析,表明上述分析方法和判據正確可行;為提高產品抗隨機振動能力,本文從兩個方面給出5種降低產品內部PCBA隨機振動動力學響應的方案,通過有限元分析模型(包括機箱結構、PCBA、連接件、元器件等),計算結果表明:產品內部PCB板位移響應和變形極限大幅降低80%以上,滿足判據要求,產品經試驗驗證了減振措施的有效性。

(References)

[1]馬興瑞,韓增堯.衛星與運載火箭力學環境分析方法及試驗技術[M].北京:科學出版社,2014 Ma Xingrui,Han Zengyao.Mechanics analysis method and testing technology for satellite and launcher[M]. Beijing:Science Press,2014(in Chinese)

[2]H D Solomon.Low cycle fatigue of surface-mounted chip-carrier/printed wiring board joints[J].IEEE Transactions on Components and Packaging Technologies,1989,12(4):473-479

[3]Reza Ghaffarian,Namsoo P Kim.Reliability and failure analysis of thermally cycled ball grid array assemblies [J].IEEE Transactions on Components and Packaging Technologies,2000,23(3):528-534

[4]Dave S Steinberg.Vibration analysis for electronic equipment[M].Second Edition.New York:A Wiley-Interscience Publication,John Wiley,1989

[5]Pitarrsi J M,Di Edwardo A V.A design approach for the systematic improvement of support locations for vibrating circuit cards[J].ASME J Electron Package,1993,115(1):118-12

[6]Pitarresi J M,Primavera A A.Comparison of modeling techniques for the vibration analysis of printed circuit cards[J].ASME J Electron Package,1992,114(4):378-38

[7]Dave S.Steinberg.電子設備振動分析[M].王建剛,譯.北京:航空工業出版社,2012 Vibration Analysis for Electronic Equipment[M].Wang Jiangang,translated.Beijing:Aviation and Industry Press,2012(in Chinese)

[8]劉孝保,杜平安.一種面向動態分析的PCB板等效建模方法[J].儀器儀表學報,2011,32(4):864-869 Liu Xiaobao,Du Ping'an.Dynamic property analysisoriented PCB equivalent modeling method[J].Chinese Jouranl of Scientific Instrument,2011,32(4):864-869 (in Chinese)

[9]范文杰.星載電子設備寬頻隨機振動響應分析[J].電子機械工程,2010,26(4):5-17 Fan Wenjie.High-frequency random vibration response analysis of spaceborne electronic equipment[J].Electro-Mechanical Engineering,2010,26(4):5-17(in Chinese)

[10]安彤,秦飛.跌落沖擊載荷下錫焊接點金屬間化合物層的動態開裂[J].固體力學學報,2013,32(2):117-124 An Tong,Qin Fei.Dynamic fracturing of the intermetallic compound layer in solder joints under impact drop loading[J].Chinese Journal of Solid Mechanics,2013,34(2):117-124(in Chinese)

[11]郭強,趙玫.隨機振動條件下SMT焊點半經驗疲勞壽命積累模型[J].振動與沖擊,2005,24(2):24-26 Guo Qiang,Zhao Mei.Fatigue life model for SMT solder joints under random vibration[J].Journal of Vibration and Shock,2005,24(2):24-26(in Chinese)

(編輯:張小琳)

Design to Resist Random Vibration Environment for Space Electronic Device

LIU Chen LIU Tianxiong JIANG Wanjie FAN Yanping
(Beijing Institute of Spacecraft System Engineering,Beijing 100094,China)

Space electronic devices suffers extreme and complicated random vibration in the process of rocket lifting off and ascending,which causes failure of device.According to the principle and characteristic of random vibration on PCB,the failure mechanism and five common principles for resisting random vibration environment are investigated,and two analysis methods and risk criteria are summarized,which are proved correct by FEM analysis.Then it is found out that five measures to decrease dynamics response are proved available with the displacement on PCB reduced obviously.The analysis methods,estimation criteria and the five measures can be as reference for estimating the capability of resisting random vibration for the similar electronic device. Key words:space electronic device;random vibration;analysis method;design to resist random vibration;vibration damping

V476.3

A

10.3969/j.issn.1673-8748.2016.03.013

2016-01-12;

2016-04-28

國家重大航天工程

劉晨,男,工程師,從事航天器總體設計工作。Email:jackychenjob@163.com。

猜你喜歡
振動產品
振動的思考
科學大眾(2023年17期)2023-10-26 07:39:14
好產品,可持續
現代裝飾(2022年4期)2022-08-31 01:39:32
從靈感出發,邂逅好產品
現代裝飾(2022年3期)2022-07-05 05:55:06
噴水推進高速艇尾部振動響應分析
This “Singing Highway”plays music
振動攪拌 震動創新
中國公路(2017年18期)2018-01-23 03:00:38
中立型Emden-Fowler微分方程的振動性
2015產品LOOKBOOK直擊
Coco薇(2015年1期)2015-08-13 02:23:50
UF6振動激發態分子的振動-振動馳豫
計算物理(2014年2期)2014-03-11 17:01:44
新產品
玩具(2009年10期)2009-11-04 02:33:14
主站蜘蛛池模板: 中文字幕亚洲另类天堂| 91视频99| 伊在人亚洲香蕉精品播放| 国产日本视频91| 国产亚洲现在一区二区中文| 亚洲va视频| 六月婷婷精品视频在线观看| 欧美午夜小视频| 亚洲男人在线天堂| 午夜在线不卡| 国产精品久久久久久久久kt| 国产毛片久久国产| 99热精品久久| 国产精品极品美女自在线网站| 精品无码国产自产野外拍在线| 国产96在线 | 日本免费精品| 高清久久精品亚洲日韩Av| 激情视频综合网| 国产久操视频| 黄色网页在线播放| 999国内精品视频免费| 中国毛片网| 中国一级特黄视频| 日韩在线播放欧美字幕| 国产精品蜜臀| 国产精品欧美在线观看| 亚洲男人的天堂网| 激情综合网址| 欧洲av毛片| 亚洲无码高清一区二区| 中文字幕av无码不卡免费| 无码中文字幕精品推荐| 18禁黄无遮挡网站| 欧美日本激情| 欧美另类视频一区二区三区| 538国产视频| 色噜噜狠狠色综合网图区| 久久国产免费观看| 国产va在线| 91综合色区亚洲熟妇p| 欧美日本一区二区三区免费| 国产成人无码Av在线播放无广告| 国产chinese男男gay视频网| 国产精品香蕉在线观看不卡| 成年人免费国产视频| 天天摸夜夜操| 黄色网址免费在线| 青青操视频在线| 日本福利视频网站| 黄色一级视频欧美| 丰满人妻中出白浆| 国产h视频免费观看| 思思热精品在线8| 天天摸夜夜操| 国产成人综合亚洲欧洲色就色| 国产精品成人AⅤ在线一二三四| 免费看美女自慰的网站| 国产精品亚洲片在线va| 国产精品第一区| 国产情侣一区| 夜夜爽免费视频| 91精品国产情侣高潮露脸| 国产白浆在线观看| 在线观看国产精品日本不卡网| 免费毛片a| 丁香综合在线| 91热爆在线| 免费看黄片一区二区三区| 亚洲精品午夜无码电影网| 亚洲国产精品日韩av专区| 国产主播在线观看| 欧美性精品| 色偷偷综合网| 欧美一区福利| 色偷偷男人的天堂亚洲av| 成年A级毛片| 亚洲国产看片基地久久1024| 狼友av永久网站免费观看| 人妻一本久道久久综合久久鬼色| 国内老司机精品视频在线播出| 国产噜噜在线视频观看|