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

集中式驅動純電動車振動特性試驗研究

2015-06-04 13:03:26王曉華
振動與沖擊 2015年14期
關鍵詞:振動

于 蓬,王曉華,章 桐,3,孫 玲,郭 榮

(1.同濟大學 新能源汽車工程中心,上海 201804;2.同濟大學 汽車學院,上海 201804;3.同濟大學 中德學院,上海 201804)

通常認為電動車振動、噪聲小于傳統車。由于驅動源內燃機被替代,原有“掩蔽效應”隨之消失,一些車輛運行中不易察覺的振動、噪聲問題凸顯,如懸置系統瞬態振動、傳動系統扭轉振動、電機嘯叫、齒輪嘯叫、電磁噪聲等,構成電動車NVH的新特性,使汽車NVH研究領域向外延展,給汽車設計帶來新挑戰[1]。

近年來,已有對新能源汽車的振動、噪聲問題進行的試驗研究[2-14]。而研究對象有混合動力汽車、純電動汽車、燃料電池汽車等,目標基本相同,即針對新能源車的NVH問題試驗、分析,闡明振噪特性。因此,量產電動車的NVH問題已成關注重點。

本文以某集中驅動式純電動車為研究對象,進行全面振動及模態試驗,據試驗結果總結其有別于傳統車的新特性,為解決相關工程問題提供試驗支持。

1 試驗用車及動力傳動系統

設整車坐標系為:車輛前進的反方向為X軸正向,駕駛員右手側為Y軸正向,豎直向上為Z軸正向。該車采用前置前驅式形式,驅動力由電機輸出,經減速器、差速器輸出到左右半軸驅動車輛前進。電機動力總成由電機與減/差速器組成,并通過3個懸置與副車架相連。選整車坐標系OXYZ為參考坐標系,在電機動力總成質心處建立質心坐標系OcXYZ,坐標軸分別與參考坐標系OXYZ坐標軸平行、同向。該車電機動力總成結構及質心坐標系見圖1,動力傳動系統模型(電機轉子到車輪)見圖2。

2 試驗方案

2.1 試驗目的

圖1 電機動力總成三維模型Fig.1 3D model of the motor driven power train

圖2 傳動系統模型Fig.2 Drive train model

據純電動車特征,采集各工況(穩態、瞬態、反拖、怠速)振動信號,記錄電機電流、轉速、轉矩及車速等信息。通過對試驗結果處理全面獲取電動車的振動信號。信號采集由5個布設于動力總成表面、6個布設于懸置系統支架主被動端、1個布設于座椅導軌處、1個布設于懸掛上止點及1個布設于方向盤處共14個三向加速度傳感器完成;電機轉矩、轉速采集由CANCASE獲取,電流信號由示波器獲得。

2.2 試驗設備

試驗在消聲轉鼓試驗室進行,試驗現場見圖3。設備為:① 純電動車樣車;② 5個三向加速度傳感器、27個單向加速度傳感器;③ 記錄、分析、處理振噪信號設備LMS Testlab;④ 記錄轉矩轉速車速信號設備CANCASE/CANnape;⑤ 記錄驅動電機電流信號設備LeC-roy HDO4034高分辨率示波器;⑥ 記錄試驗過程照相機。

圖3 整車聲振試驗現場Fig.3 Vibration and noise test

2.3 試驗工況

試驗過程按國家標準GB/T14365-93《機動車輛噪聲測量方法》布置振動加速度傳感器,主要測試工況有:①“怠速”工況:接通電源電機不驅動,用于配合查找振動源;② 穩速工況:在10~80 km/h穩速巡航條件下每間隔10 km/h測量一次;③ 瞬態工況:選取0~10/20/30/40(POT)/50/(POT)60(POT)/70/(WOT)80(WOT)km/h振動樣本,便于全面分析整個加速過程中振動特性,記錄電機轉矩轉速車速在車輛突加突減情況下測量信號,用于分析傳動系統扭轉振動;④轉鼓拖動工況:轉鼓帶動車輪速度從10~80 km/h巡航,每間隔10 km/h測量一次,用于配合查找振動源。另外,試驗前充滿電,全過程記錄SOC值,避免因電池性能大幅波動造成測試誤差。

3 試驗結果

3.1 激勵階次及傳遞路徑分析

3.1.1 電機殼體測點處階次

進行0~80 km/h瞬態加速工況下加速度時域響應獲取,并對測點信號進行階次分析。以電機處某測點數據處理結果為例,見圖4。由圖4看出,轉子基頻4、8、10、12、16、20、24、40、48 及 58 倍 10 個階次振動幅值較大。其中,4、8、12、16、20、24、48 倍對應驅動電機(電機極對數 p=4)變頻器供電頻率的 1、2、3、4、5、6、12倍[15]:引起1倍電流基頻原因為電流檢測時直流偏移檢測誤差、恒定電流成分、有限編碼器分辨率限制等;引起2、4倍頻原因為電流檢測時增益失配、相位偏差及諧波電流影響等;引起3、5倍頻原因為偶數倍諧波電流;引起6、12倍頻原因為諧波電流、定轉子氣隙磁場分布不均勻及PWM逆變器開關頻率影響等。諸多偏差反饋到傳感器中造成轉矩計算偏差,從而產生輸出轉矩階次波動。而10、20、40、58倍階次振動與減/差速器結構及速比i有關。其中20階次與電機供電電流諧波、減差速器均相關。

為進一步揭示由電流因素引起的電機轉矩波動,取車速v=40 km/h穩態工況,某向電流時域變化曲線及頻譜分析結果見圖5,此工況對應電流基頻f=178 Hz(f=vip/(7.2πr),r為車輪半徑)。此外在 356 Hz、712 Hz、892 Hz、1247 Hz、2317 Hz等電流諧波頻率處也存在峰值,分別為電流的 1、2、4、5、7、13 倍頻率,將分別產生轉子基頻 4、8、16、20、28、42 倍的轉矩波動。結合圖4幅值可知,1247 Hz、2317 Hz處電流幅值雖較大,但由于階次較高,產生的轉矩波動幅值不及低階次電流諧波影響及減差速器結構因素引發階次振動幅值大。

圖4 電機某測點加速度階次圖Fig.4 Acceleration order of motor point

3.1.2 減/差速器處階次圖

仍以獲得0~80 km/h瞬態加速工況減差速器處某測點數據處理結果,見圖6。由圖6看出,10、20、29、58四個階次振動幅值較大,其中58、29階次幅值最大,10、20階次次之。29階顯著原因為一級齒輪副小齒輪齒數為29;58階顯著原因為兩級齒輪副均存在一定偏心,而偏心會引發主要振動階次2倍頻階次振動;10階次出現對應一級小齒輪數與一級減速比比值;20階出現因差速器中半軸輪齒數與行星輪齒數比值為2,屬于10階次的2倍頻。幅值最大為58階次,可初步診斷減速器階次振動原因有一級齒輪副嚙合誤差、輸入軸齒輪磨損不均勻造成質量偏心、輸入軸軸承裝配不對中造成動不平衡等。

3.1.3 座椅導軌處某向階次圖

以上 4、8、10、12、16、20、24、29、40、48、58 振動階次均由激勵源動力總成產生的階次振動,4、8、12、16、24、48階次主要與電機(包括電機氣隙磁場、電路、電控等因素)相關,10、20、29、40、58 階次主要與減差速器結構及減速比有關。并非所有階次振動均能引發車身關注點的振動。駕駛員座椅導軌處階次振動見圖7,由此可對激勵的傳遞路徑進行初步探討。圖中幅值較大的4、8、12、24階次主要由電機因素(尤其電流諧波及磁場氣隙非正弦)引起,29、58兩階次在減速器表面較顯著而車身上并未體現,而10、20、40三個階次由于本身幅值較小,在車身測點的振動不明顯。由此可初步判斷動力總成產生的各階次通過不同路徑傳輸到駕駛員。

圖5 穩態工況下某向電流Fig.5 A phase current under steady condition

圖6 減差速器某測點加速度階次圖Fig.6 Acceleration order of reducer point

圖7 座椅導軌處加速度階次圖Fig.7 Order acceleration at the seat rail

圖8 動力總成右懸置主被動端階次切片Fig.8 Order slices of right power train mount

為更深入討論,將懸置系統右懸置點的主被動端進行“階次切片”處理,見圖8。由圖8看出,切片圖能基本反映動力總成的各主要階次,如4、8、12、24階由電機轉矩波動引起,20、29由齒輪嚙合等結構引起;而在被動端,4、8、12等轉矩波動階次并不存在,說明此類激勵被衰減或經其它路徑傳遞。在座椅導軌處明顯出現此類階次,可判斷由電機控制因素引起的激勵中,低階次(12階以下)部分通過傳動系統(電機轉子-齒輪系-半軸-輪胎 -車身)傳遞;被動端出現18、21、25階次分別為20階次、24階次邊頻,為因結構因素引起[16],可對應20、24階次進行探討,結合主動端及座椅導軌處階次圖知,由電機、減差速器共同產生的20階次通過懸置系統被衰減后經副車架進行二級衰減,在車身的作用較小;高階次的29、58幾乎由懸置系統完全衰減,不會作用于車身;由電機控制因素引起的激勵中,高階次(24階)部分主要通過懸置系統傳遞至車身。因此,懸置系統對20階次以下激勵有一定衰減作用,但此類激勵主要通過傳動系統傳輸,25階次以上激勵主要通過懸置系統傳遞,但懸置系統在此范圍內具有較好的衰減作用,會避免向車身的傳遞。而對20-25階次范圍內的振動則分別通過傳動系統、懸置系統傳遞。由圖8可知,懸置系統在起步、中低速(電機轉速1000 r/min、2000 r/min 對應車速 15 km/h、30 km/h)時產生的沖擊、振動具有較好的衰減作用,但在整個頻段傳遞率并不理想,可按需要進一步優化設計。

3.2 瞬態問題及整車振動舒適性評價

車輛起步加、減速時均會發生整車抖動,主要由電機瞬態轉矩特性及傳動系統固有特性引起。加、減速工況下電機轉矩及轉速波動見圖9。由圖9看出,突踩踏板時電機轉矩從零突變,驅動力矩及負載的相互作用引起傳動系沖擊、振動。沖擊表現為轉矩波動、電機轉速震蕩、車速波動及整車縱向抖動;振動表現為在沖擊未完全消除時由電機各階次轉矩波動引起的整個傳動系統共振,共振頻率較低,恰在人體對水平前后方向最敏感的0.5~2 Hz范圍內;能量回饋終止標定值造成1000 r/min處波動。此振動與沖擊均會影響駕乘的舒適性。

據ISO2631-1:1997(E)標準規定,峰值系數>9時用輔助評價方法即振動劑量(VDV)評價,可更好估計偶爾遇過大脈沖所致高峰值系數振動對人體影響。處理座椅導軌處加速度瞬態信號,獲得振動計量結果見表1。其中aw(t)為加權加速度時間歷程,aw為加權加速度均方根值。由表1看出,突加速時隨目標車速提高,舒適性下降;突減速時隨穩定初始車速提高,舒適性下降;突減速振動遠大于突加速,極不舒適;豎直X向振動大于水平(前后)Y向。由于緩慢加減速時振動不突出,水平(左右)Y向加速度值較小,表1未給出。對瞬態加減油門的振動問題,可考慮通過整車控制器控制目標電機轉矩或在傳動系統加裝扭轉減震器進行初步衰減,而優化懸置系統、懸掛系統、傳動系統等子系統可進一步減小整車振動及懸置系統的動反力幅值。

圖9 急加急減工況信號Fig.9 Tip-in and tip-out signal

表1 瞬態工況舒適性指標計算Tab.1 Comfort index calculation under transient conditions

3.3 相關模態試驗結果

進行相關系統及部件的模態試驗,包括傳動系統、懸置系統、懸掛系統、輪胎、轉向系統等,限于篇幅,試驗過程不再贅述,試驗結果見圖10。該圖全面反映出電動車的模態特性,將其與激勵階次特性結合,則可分析傳動系統引起的振動及聲輻射噪聲問題。圖中縱坐標電機轉速750 r/min對應車速11.5 km/h,375 r/min對應車速5.8 km/h,125 r/min對應車速2 km/h。由圖10看出,① 該車輛傳動系統機械結構從模態分布角度較理想,各系統部件間均有一定間隔,無嚴重的模態耦合現象。因此在解決整車出現的聲振問題時,可考慮多從整車控制、電機控制角度(如優化瞬態轉矩輸出曲線)進行優化,而添加機械部件或改變系統結構的傳統減振降噪方法可以作為第二選擇。② 該款電機動力總成懸置系統固有頻率較高,主要分布在30~100 Hz,普遍高于傳統內燃機汽車動力總成懸置系統的5~17 Hz,此為車輛設計初期為避免電機瞬態特性造成過大位移而進行的有意嘗試。結果表明,提高懸置剛度具有較好的限制瞬態工況動力總成位移過大效果,且在優化電動車動力總成模態分布時具有更多的模態布置空間,不受類似傳統車怠速工況限制;但中高頻振動傳遞率會下降,與主被動端切片分析結果一致,有需要可進行進一步折衷、優化。③ 考慮電磁剛度的傳動系統“零階”固有頻率為電傳動系統及內燃機傳動系統的區別之一;該車采用框架式車身,二級減振結構,在傳統車輛中較少使用;動力總成殼體模態特性與傳統內燃機動力源也有區別;而輪胎、轉向系、懸架系統則區別不大。

圖10 模態頻率排布圖Fig.10 Modal frequency alignment

4 電動車振動特性與傳統車異同

4.1 動力傳動系統振動現象匯總

車輛按系統級別可劃分為整車、系統、子系統、零部件4個級別,振動評價從整車評價、系統評價、部件評價三方面進行。電動車動力/傳動系統的振噪權重在整車總權重比例次于胎噪、風噪[17],電動車振動現象較內燃機有所減少,引發振動的各部件具有顯著特點,從部件劃分角度進行匯總能獲得全面認識[18-20]。將動力/傳動系統分為電動機、傳動系、減差速器、總成殼體、懸置系統5部分,由部件角度匯總電動車動力總成傳動系統振動現象,見表2。

表2 電動車動力總成傳動系統振動現象匯總Tab.2 Summary of vibration phenomenon of electric power train

4.2 激勵階次特性圖

綜合動力總成的階次特性,兼顧傳動系統靜、動偏心,繪制動力總成階次特性圖見圖11。圖中實線框實線為傳統內燃機常用階次;虛線框實線為傳統內燃機車輛可拓展階次(研究小幅高頻振動時用);無框線實線為純電動車驅動電機激勵階次,包括正向驅動(0~9000 r/min)及反向拖動(9000~0 r/min)兩部分;虛線為減速器嘯叫階次。

由圖11看出,電動機階次特性與內燃機車輛存在明顯差別。如發動機常用工作轉速范圍為800 ~6000 r/min,而電動機常用轉速范圍遠高于發動機;發動機主要階次由質量力、質量力矩及氣體扭矩組成,且幅值大于機械不平衡等影響,電動機轉矩波動階次主要為極對數倍頻;減/差速器階次振動,由于“無掩蓋”效應,也應引起關注。

4.3 振動幅值特性圖

某內燃機車與電動車動力傳動系統各子系統振動噪聲水平比較[21-23]見圖 12。圖中 EV(Electric Vehicle)代表電動車,IC(Internal Combustion)代表內燃機車輛,f0為固有頻率(懸置、傳動系統或總成殼體固有頻率)、A為振動或噪聲幅值。對比發現,新能源車輛優點為:振動激勵源數量明顯減少,各傳遞路徑間相互影響減少,利于降低整車振動水平;但結構簡化、輕量化要求提高及振動噪聲傳遞路徑變化給新能源車減振降噪提出新挑戰。

圖12 動力傳動系統振動噪聲分布Fig.12 Vibration and noise distribution of driveline

結合已有研究,將新能源汽車振動問題總結為5大特點,即① 無掩蔽效應:由于無類似內燃機氣體燃燒噪聲、慣性質量力/力矩產生振動等掩蔽,原在內燃機車輛中不易被察覺的次要振動、噪聲將成為必須關注的問題。② 瞬態振動更劇烈:更迅速的電動機動態響應、更小的動力總成系統質量/慣量、無扭轉減震器的常嚙合傳動系統,將產生更大的軸系扭轉振動及動力總成懸置系統扭轉振動問題。③ 存在機電耦合現象:新能源汽車傳動系統為典型的機電一體化系統,其動力傳動系具有一般電傳動系統普遍存在的機電耦合特性,會影響電動車傳動系統各頻段扭轉振動特性。④ 電磁振動突出:尤其表現為高速時電機嘯叫,電磁力波與總成殼體的磁固耦合振動問題突出。⑤ 齒輪撞擊、嘯叫:電機動態響應快、齒輪系統結構簡化、輪齒及軸系轉速范圍擴大將引發齒輪系統振噪問題的新特性。

5 結論

(1)通過試驗找到動力總成振動的主要階次及相關傳遞路徑,即動力總成產生轉子基頻4、8、10、12、16、20、24、29、40、48、58 倍階次振動,4、8、12、16、24、48 階次主要與電機(包括電機氣隙磁場、電路、電控等因素)相關,10、20、29、40、58 階次主要與減差速器結構、減速比有關。座椅導軌處4、8、12、24階次主要通過傳動系統而非懸置系統傳遞到駕駛員,對駕乘舒適性具有較顯著影響。

(2)評價該車動力總成懸置系統性能,即對20階次以下及25階次以上激勵具有較好的衰減作用,避免其向車身傳遞;而對20~25階次范圍內的振動隔振效果不理想。對車輛起步(0~5 km/h)、中低速(15 km/h,30 km/h)時產生的沖擊及傳動系統共振具有較好的衰減作用,但在整個頻段傳遞率并不理想。

(3)電機輸出瞬態特性造成的整車傳動系統沖擊且突減轉矩(Tip-out)大于突加轉矩(Tip-in),均可通過優化整車控制器目標轉矩得到緩解,該方法已在實車上驗證其可行性。

(4)由于無類似發動機等強烈振動源掩蓋,減/差速器的階次振動突顯,傳動系統靜、動偏心產生的主要階次二倍階次較顯著,造成電動車較嚴重嘯叫。可通過減小一、二級齒輪副小齒輪對傳遞誤差的靈敏度、輪齒齒面微觀修形及提高兩級齒輪副裝配精度進行改進。

(5)模態排布結果便于對集中驅動式電動車傳動系統結構進行全面把握,整車各子系統模態排布較合理;振動現象匯總、激勵階次特性、幅值特性分析闡明該電動車與傳統內燃機車輛相關特性區別,為進一步診斷、懸置/懸掛等子部件振動控制及優化提供試驗支持。

[1]于蓬,賀立釗,章桐,等.集中電機驅動車輛動力傳動系統NVH性能研究現狀與展望[J].機械設計,2014,31(3):1-5.YU Peng, HE Li-zhao, ZHANG Tong, et al. NVH performance study of centralized motor driven vehicle power train[J].Journal of Machine Design,2014,31(3):1 -5.

[2]Kuang M L.An investigation of engine start-stop NVH in a power split powertrain hybrid electric vehicle[R].SAE Technical Paper,2006.

[3]Humbert L,Pellerey P,Cristaudo S.Electromagnetic and structural coupled simulation to investigate NVH behavior of an electrical automotive powertrain[R].SAE Technical Paper,2012.

[4]Freeman T, Thom B, Smith S. Noise andvibration development for adapting a conventional vehicle platform for an electric powertrain[R].SAE Technical Paper,2013.

[5]Govindswamy K, Wellmann T, Eisele G. Aspects of NVHintegration in hybrid vehicles[R].SAE Technical Paper,2009.

[6]Rabofsky A,Walther C,Rothdeutsch J,et al.Electrical andhybrid drive lines:new demands for N&V integration[R].SAE Technical Paper,2010.

[7]La C,Poggi M,Murphy P,et al.NVH considerations for zero emissions vehicle driveline design[R].SAE Technical Paper,2011.

[8]Mammetti M,Arroyos M R.Thestudy of torque control characteristics for the optimization of the NVH of an electric vehicle[R].SAE Technical Paper,2014.

[9]閆明剛,張勇,侯之超,等.混合動力客車傳動系扭振響應及其影響因素分析[A].第22屆全國結構工程學術會議論文集第Ⅲ冊,2013:53-57.

[10]童煒,侯之超.輪轂驅動電動汽車垂向特性與電機振動分析[J].汽車工程,2014,36(4):398-403.TONG Wei,HOU Zhi-chao. Analyses on the vertical characteristics and motor vibration of an electric vehicle with motor-in-wheel drive[J].Automotive Engineering,2014,36(4):398-403.

[11]黃孝慈,唐小林,于海生.功率分流式混合動力轎車傳動系噪聲分析[J].機械傳動,2013,11:130-133.HUANG Xiao-ci,TANG Xiao-lin,YU Hai-sheng.Noise analysis of a power-split hybrid electric vehicle driveline[J].Journal of Mechanical Transmission,2013,11:130 -133.

[12]朱宇.純電動汽車車內聲品質分析評價研究[D].吉林:吉林大學,2013.

[13]郭榮,萬鋼,左曙光,等.燃料電池轎車主要噪聲源識別的試驗研究[J].汽車工程,2007,29(5):377-380.GUO Rong,WANG Gang,ZUO Shu-guang,et al.An experimental study on noise sources identification for fuel cell vehicle[J].Automotive Engineering,2007,29(5):377-380.

[14]相龍洋,左曙光,何呂昌,等.電動汽車及其驅動永磁電機聲振特性試驗[J].振動.測試與診斷,2014,34(1):90-95.XIANG Long-yang,ZUO Shu-guang,HE Lü-chang,et al.Experimental study on acoustic and vibration characteristics of electric vehicle and its permanent drive motor[J].Journal of Vibration,Measurement& Diagnosis,2014,34(1):90 -95.

[15]于蓬,章桐,郭榮.車輛傳動系扭轉振動主動控制研究現狀及趨勢[J].機電一體化,2014,6:3-10.YU Peng,ZHANG Tong,GUO Rong.The current state and trends of automotive drive train torsional vibration active control[J].Mechatronics,2014,6:3-10.

[16]Dupont J B,Aydoun R,Bouvet P.Simulation of the noise radiated by an automotive electric motor:influence of the motor defects[R].SAE Technical Paper,2014.

[17]Goetchius G.Leading the charge-the future of electric vehicle noise control[J].Sound & Vibration,2011,4:5 -7.

[18]Yoshioka T, Sugita H. Noise and vibration reduction technology in hybrid vehicle development[J].Fuel,2001,2013:4-4.

[19]Komada M,Yoshioka T.Noise and vibration reduction technology in new generation hybrid vehicle development[J].SAE Transactions,2005,114(6):2648 -2655.

[20]Eisele G,Wolff K,Wittler M,et al.Acoustics of hybrid vehicles[J].SAE paper,2010:1401-1402.

[21]于蓬,章桐,陳霏霏,等.集中式驅動電動車動力總成懸置系統模態影響因素分析[J].中國工程機械學報,2014,6(1):471-477.YU Peng,ZHANG Tong,CHEN Fei-fei,et al.Investigation on modeling method of powertrain mounting system of central driven electric vehicle[J].Chinese Journal of Construction Machinery,2014,6(1):471 -477.

[22]于蓬,陳霏霏,章桐,等.集中驅動式純電動車動力總成噪聲特性分析[J].機電一體化,2015,1:7-12.YU Peng,CHEN Fei-fei,ZHANG Tong,et al. Acoustic characteristic analysis of central driven power train of electric vehicle[J].Mechatronic,2015,1:7 -12.

[23]于蓬,陳霏霏,章桐,等.集中驅動式電動車動力總成系統振動特性分析[J].振動與沖擊,2015,34(1):44 -48.YU Peng,CHEN Fei-fei,ZHANG Tong,et al.Vibration characteristics analysis of a central-driven electric vehicle powertrain[J].Journal of Vibration and Shock,2015,34(1):44-48.

猜你喜歡
振動
振動的思考
科學大眾(2023年17期)2023-10-26 07:39:14
某調相機振動異常診斷分析與處理
大電機技術(2022年5期)2022-11-17 08:12:48
振動與頻率
天天愛科學(2020年6期)2020-09-10 07:22:44
This “Singing Highway”plays music
具非線性中立項的廣義Emden-Fowler微分方程的振動性
中立型Emden-Fowler微分方程的振動性
基于ANSYS的高速艇艉軸架軸系振動響應分析
船海工程(2015年4期)2016-01-05 15:53:26
主回路泵致聲振動分析
UF6振動激發態分子的振動-振動馳豫
計算物理(2014年2期)2014-03-11 17:01:44
帶有強迫項的高階差分方程解的振動性
主站蜘蛛池模板: 成人国内精品久久久久影院| 国产精品久久精品| 亚洲国产精品一区二区高清无码久久 | 色亚洲成人| 亚洲三级网站| 91精品专区国产盗摄| 欧美在线一二区| 欧美精品v欧洲精品| 国产免费久久精品99re不卡| 青青草a国产免费观看| 黑人巨大精品欧美一区二区区| 欧洲亚洲欧美国产日本高清| 国产精品亚洲综合久久小说| 久久亚洲黄色视频| 国产小视频免费| 亚洲精品另类| 网友自拍视频精品区| 日韩人妻少妇一区二区| 亚洲欧美人成电影在线观看| 色综合激情网| 99re热精品视频国产免费| 自拍亚洲欧美精品| 波多野结衣一区二区三区四区| 九九免费观看全部免费视频| 99热国产这里只有精品无卡顿"| 亚洲福利一区二区三区| 三上悠亚在线精品二区| 看av免费毛片手机播放| 亚洲中久无码永久在线观看软件 | 四虎成人免费毛片| 韩国v欧美v亚洲v日本v| 国产理论最新国产精品视频| 免费国产一级 片内射老| 欧美激情综合| 亚洲国语自产一区第二页| 萌白酱国产一区二区| 国产午夜一级毛片| 久久99国产乱子伦精品免| 自拍偷拍一区| 精品福利视频网| 青草精品视频| 亚洲精品色AV无码看| 国产一区二区三区免费观看| 国产三级视频网站| 亚洲精品麻豆| 91丝袜在线观看| 国产另类乱子伦精品免费女| 97狠狠操| 亚洲天堂视频在线免费观看| 日日噜噜夜夜狠狠视频| 国产欧美日韩18| 99精品国产自在现线观看| 欧美日本二区| 中文字幕在线永久在线视频2020| 在线精品视频成人网| 成人av手机在线观看| 国产剧情国内精品原创| 一本大道香蕉中文日本不卡高清二区 | 久久亚洲国产最新网站| 欧美劲爆第一页| 自拍偷拍欧美日韩| 精品人妻AV区| 久久国产精品影院| 欧美成人一级| 日本国产精品一区久久久| 在线国产91| 亚洲日本韩在线观看| 欧美v在线| 亚洲毛片一级带毛片基地| 色综合热无码热国产| 欧美不卡二区| 久久国产精品国产自线拍| 色综合网址| 天堂网亚洲系列亚洲系列| 青青草综合网| 日韩欧美国产成人| 五月天香蕉视频国产亚| 免费jjzz在在线播放国产| 午夜视频在线观看免费网站| 亚洲美女一区二区三区| 欧美日韩国产高清一区二区三区| 亚洲区欧美区|