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

    燃料電池車用超高速空壓機永磁體結構強度研究

    2021-01-07 11:23張智明武賽潘佳琪倪玥章桐
    湖南大學學報·自然科學版 2021年12期
    關鍵詞:空壓機

    張智明 武賽 潘佳琪 倪玥 章桐

    摘 ? 要:針對大功率燃料電池系統(tǒng)高速空壓機轉子永磁體離心拉應力破壞,采用解析法計算了額定轉速時護套與永磁體間過盈量的大小及應力,采用有限元方法建立二維軸對稱模型,計算了永磁體和護套內的應力分布,對比解析法與有限元方法的準確性. 基于有限元法研究了溫升和護套材料對永磁體應力的影響,結果表明,解析法和有限元法計算的各應力的相對誤差不超過2.5%,而影響因素中轉速和溫升均會造成永磁體應力的顯著提升,需要增大過盈量對轉子永磁體保護. 碳纖維護套相比合金護套可以至少降低42.1%的永磁體應力,且在高溫工況下的效果更為顯著. 降低轉速、增大裝配過盈量以及改善轉子冷卻均能有效優(yōu)化結構強度,而復合材料護套相比鋼護套更適用于保護永磁體.

    關鍵詞:燃料電池汽車;空壓機;永磁體;過盈配合;強度分析

    中圖分類號:U473.4 ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻標志碼:A

    Study on Mechanical Strength of Permanent Magnet

    for Super High-speed Air Compressor in Fuel Cell Vehicle

    ZHANG Zhiming?,WU Sai,PAN Jiaqi,NI Yue,ZHANG Tong

    (School of Automotive Studies,Tongji University,Shanghai 201804,China)

    Abstract:Aiming at the centrifugal tensile stress failure of permanent magnets on the rotor for super high-speed air compressor in high-power fuel cell system, the magnitude of interference between the sleeve and the permanent magnet as well as the stress at the rated speed are calculated by the analytical method. Then, a two-dimensional axisymmetric model was established by the finite element method to calculate the stress distribution, and the accuracy of the analytical method and the finite element method was compared. Furthermore, the influence of temperature rise and sleeve material on stress is studied based on the finite element method. The results show that the calculations of the two methods differ no more than 2.5%. Among the influencing factors, both the rotational speed and the temperature rise cause a significant increase in the stress of the permanent magnet, and it is necessary to increase the interference in order to provide protection for the permanent magnet. Additionally, the carbon fiber sleeve can reduce the stress of the permanent magnet by at least 42.1% when compared with the alloy sleeve, and the effect is particularly obvious under high temperature conditions. Therefore, reducing the rotating speed, increasing the assembly interference and improving rotor cooling can effectively optimize the structural strength, and the composite material sleeve is more suitable for protecting permanent magnets than the steel sleeve.

    Key words:fuel cell vehicle;air compressors;permanent magnet;interference fit;strength analysis

    燃料電池發(fā)動機系統(tǒng)研發(fā)時工作壓力應選擇在2×105 Pa附近[1],相關研究也表明電池性能會隨著工作壓力的增大而提升,但當壓力超過2 × 105 Pa后,電池性能的提升明顯減緩,同時濃差損失開始增加,目前最常用的燃料電池系統(tǒng)工作壓力在2 × 105 Pa左右[2]. 燃料電池供氣系統(tǒng)中空壓機的功耗占整個輔助系統(tǒng)的90%,約為整體輸出功率的13%,成本占系統(tǒng)總體的16%[3]. 空壓機作為供氣系統(tǒng)的核心部件是燃料電池系統(tǒng)除負載以外最大的附加能耗,其綜合性能對整個燃料電池發(fā)動機系統(tǒng)的輸出功率有著決定性影響. 因此研究可靠且低成本的空氣供應部件是當前燃料電池發(fā)動機研發(fā)中的一個重點.

    在多種空壓機類型中,離心式空壓機在功率密度、效率和噪聲等方面具有最佳的綜合表現,又易于與渦輪相匹配以實現排氣能量的回收,被認為是未來發(fā)展的主流方向[4]. 高速永磁電機由于其功率密度高、體積小、動態(tài)響應快、結構簡單等優(yōu)點,尤其適用于安裝空間有限的車用高速空壓機. 目前一般選用高速永磁電機作為空壓機的驅動電機[5],其中又以表貼式應用最為廣泛. 表貼式高速永磁電機轉子一般由轉軸、永磁體以及保護套三個部件裝配而成,電機的輸出功率在轉矩不變的情況下與轉速成正比,在保持空壓機體積不變的情況下,要實現高壓比、大流量的供氣以增加燃料電池的輸出功率就必須要提高轉子的轉速,即朝高速化發(fā)展[3]. 目前超高速轉子一般在5 × 104 ~ 1.5 × 105 r/min以滿足30~100 kW左右燃料電池發(fā)動機的空氣供應需求. 例如搭載于本田Clarity氫燃料電池汽車中的兩級電動渦輪增壓空壓機的最高轉速為1 × 105 r/min,能提供103 W的最大輸出功率. 高速永磁電機主要選用第三代稀土永磁合金釹鐵硼制造永磁體,該材料有優(yōu)異的磁性能且可以承受很大的壓應力(約為1 000 MPa),但抗拉性能較差,僅有80 MPa. 如果未采取有效保護措施,永磁體將無法承受轉子高速旋轉時所產生的巨大離心力,從而發(fā)生斷裂破壞[5],因此轉子的結構強度是車用燃料電池空壓機轉子設計中需要解決的關鍵問題之一.

    最常采用的方法是在永磁體外設置過盈配合護套以施加預壓應力與離心拉應力相抵消的方式來保證結構強度,護套的常用材料主要分為非導磁高溫合金[6]和碳纖維復合材料[7]兩類. 張鳳閣等[8-9]在兩項研究中指出碳纖維護套在強度要求相同時所需的厚度更小,且在散熱條件相同時的溫升更低,有著出色的綜合性能. 吳震宇等[10]的研究也表明碳纖維護套的機械強度和電磁性能優(yōu)于其他材料,動力學特性也更佳. 張超等[11]研究了轉軸材料熱膨脹系數在高溫時對永磁體應力巨大的影響,指出了熱態(tài)強度校核的必要性;并在另一項研究中發(fā)現合金護套和碳纖維護套的等效應力均隨溫升線性增大,而碳纖維護套的增長率更高[12]. 馬振杰[13]分別研究了轉速和溫升對轉子結構強度的影響,結果表明溫升對應力的影響遠高于轉速,是最需要重視的影響因素.

    從保證空壓機正常工作、結構強度可靠性和運行穩(wěn)定性的角度出發(fā),對高速空壓機轉子各部件應力及其分布的研究是十分必要的. 本文分別計算了轉子在靜止和額定工作轉速兩種工況時不同半徑位置的應力,校驗了永磁體和護套所受應力的極值是否在安全范圍內,以達成設計的強度要求. 由于永磁體抗壓能力很強但抗拉能力較差,永磁體和保護套之間可以采用過盈量較大的配合,使永磁體在靜態(tài)就承受一定的預壓應力[14],以部分或完全抵消高速離心力造成的拉應力影響. 除上述研究中提及的過盈量大小和護套材料外,離心力與轉子轉速的平方成正比,與轉子的結構強度有著密切關聯,因此本文也將在其后的章節(jié)中分析該因素的影響.

    1 ? 過盈量的確定

    永磁體和護套之間可以采用過盈量較大的配合,使永磁體在靜態(tài)時能承受一定的預壓應力,過盈量的大小可通過解析計算來確定. 根據材料力學,當電機轉子的護套與永磁體間存在過盈量時,其配合表面會產生使護套內徑向外膨脹、永磁體外徑向內收縮的壓應力. 在過盈量初值的理論計算方法中,忽略護套和永磁體在軸向的長度和應變,把結構視為兩個無限長厚壁圓筒間的過盈配合,并采用拉美方程進行分析. 轉子的永磁體與護套結構尺寸如圖1所示.

    圖1中,a為永磁體內半徑,b為護套的內半徑,也是永磁體的外半徑,c為護套外半徑.

    該燃料電池空壓機電機的設計額定轉速為8×104 r /min,其轉子的護套和永磁體的尺寸參數分別為a = 11.5 mm,b = 17.5 mm,c = 19.0 mm,長度l = 0.06 m. 護套和永磁體的材料屬性如表1所示.

    1.1 ? 過盈量的計算方法分析

    工程材料通??煞譃樗苄圆牧虾痛嘈圆牧蟽纱箢? 前者以屈服為主要失效形式,極限應力為屈服極限;后者的主要失效形式是斷裂,以強度極限為極限應力. 釹鐵硼永磁體為脆性材料,應采用最大拉應力理論和最大伸長線應變理論進行校核;碳纖維護套屬于塑性材料,應采用最大剪應力理論和形狀改變比能理論進行校核. 對永磁體主要分析徑向和切向應力,對護套主要分析等效應力[15].

    對于繞軸線旋轉的圓柱體,其應力分量σθ、σr、σz依次為切向、徑向和軸向這3個方向的主應力. 在高速永磁電機的轉子護套和永磁體的過盈量的計算中,應該首先考慮滿足的條件是永磁體的切向應力σθ小于其抗拉強度,在σθ滿足條件的情況下求出過盈配合量的值,然后根據過盈量值進行高速旋轉時永磁體內應力的計算和分析. 轉子護套中的等效應力可在永磁體強度校核完成后進行驗證. 護套與永磁體的理論法計算流程如圖2所示[16-17].

    圖2中,Δδ為過盈量的減少量,δs為靜態(tài)過盈量,δd為動態(tài)過盈量,σ′?θm和σ″?θm為永磁體靜態(tài)和動態(tài)時的切向應力,p為接觸面壓應力,σm為切向應力極值,[σ]為強度極限.

    1.2 ? 過盈量的計算結果

    根據丁鴻昌等[17]提出的理論計算方法、表1和圖2,可以得出過盈配合面半徑r = b 處轉子護套的徑向位移分量u″

    e = 3.492 48 × 10-3 mm;永磁體的徑向位移分量u″

    m ?= 8.790 26 × 10-3 mm;以及兩者間的差值,即過盈量的減小量Δδ = -5.3 × 10-3 mm.

    然后計算半徑r = a處旋轉的永磁體的最大切向應力分量,其值為:σ″

    θm max = 145.54 MPa.

    由于永磁體抗拉強度[σ]為80 MPa,要保證永磁體在高速工況時不損壞就需要滿足σθ ≤[σ]/k. 其中k為安全系數,取值1.2,r = a處的靜態(tài)切向應力分量σ′

    θm = -78.87 MPa,取整-80 MPa.

    根據靜態(tài)切向應力分量σ′

    θm可以求出相應的靜態(tài)過盈裝配壓力ps = 22.73 MPa. 再通過ps求出靜態(tài)過盈量,其值為δs = 2.915 × 10-2 mm. 考慮到加工精度的要求,取靜態(tài)過盈量δs = 3 × 10-2 mm,此時靜態(tài)過盈裝配壓力為ps = 23.39 MPa.

    將靜態(tài)過盈量與過盈量的減少量相減計算動態(tài)過盈量δd = 3.53 × 10-2 mm. 再由δd計算得到對應的動態(tài)過盈壓力pd = 27.52 MPa.

    綜上,解析計算所得的結果為:護套與永磁體的靜態(tài)過盈量為0.03 mm,靜態(tài)過盈裝配壓力約為23.39 MPa,而動態(tài)過盈壓力會提升至27.52 MPa. 根據初步計算結果可以求出高速旋轉時永磁體內不同半徑位置的應力,在工作轉速時永磁體內徑向應力為-27.52 ~ 0 MPa,永磁體內切向應力為15.92 ~ 48.66 MPa.

    2 ? 轉子強度有限元分析

    2.1 ? 建立有限元接觸力學模型

    永磁體和護套的結構簡圖如圖1所示,永磁體和護套的結構均為空心圓柱體,因此可作為軸對稱問題使用彈性力學中的厚壁圓筒理論進行解析法計算,但是在采用解析法計算時忽略了永磁體和護套長度這一屬性,未考慮兩者的軸向應變,簡化該模型就會導致一定的誤差. 本節(jié)應用ANSYS軟件進行有限元法分析,有限元模型包含了結構的長度,考慮永磁體和護套在軸向的應變,因此應力計算結果會更加精確.

    永磁體與護套的尺寸參數及材料屬性都與解析法計算時相同,另設兩者的長度均為60 mm. 建立的y軸對稱模型如圖3所示.

    永磁體和護套選用PLANE 183單元來建模,此次建模為二維軸對稱模型,在定義單元時需要將Element behavior設置為Axisymmetric. 而永磁體和護套間的過盈配合可按照有限元的接觸理論設置一對線-線接觸單元來模擬,通常是成對使用CONTA 172單元和TARGE 169單元. 故轉子護套和永磁體之間的過盈配合是通過在配合面對應的線上生成CONTA 172單元和TARGE 169單元來定義的,該方法可以較為準確地模擬軸對稱模型中的過盈配合[18]. ?設置過盈量值為前文中理論計算所得的0.03 mm,分配各截面對應的材料屬性并選用映射網格對模型進行劃分.

    劃分網格之后,對整個模型建立約束. 先約束永磁體內表面和護套外表面上所有節(jié)點在徑向上的位移,再固定兩邊中間節(jié)點在軸向上的位移,設置約束后的模型如圖4所示. 最后設置轉子繞y軸的轉速為8×104 r/min

    2.2 ? ANSYS仿真結果分析

    為了分析轉子在靜態(tài)和工作轉速工況下永磁體和護套內部應力以及接觸面接觸壓力的變化,設置過盈量為0.03 mm,對轉速為8×104 ?r/min和0 r/min時的動態(tài)和靜態(tài)工況進行仿真. 所得各應力云圖如圖5和圖6所示,整理數據結果匯總于表2中.

    在轉速為8 × 104 r/min時,從接觸應力以及永磁體徑向應力和切向應力的數值可知,仿真結果與前文中解析法計算所得到的結果相近,相對誤差不大于2.5%,證明了該有限元仿真模型和計算方法的正確性. 在轉子以8 × 104 r/min的高轉速工作時,永磁體中的最大切向應力為48.2 MPa,最大徑向應力為1.5 MPa,均遠低于永磁體的抗拉強度80 MPa,故永磁體不會出現破壞;護套中的最大等效應力為395 MPa,而制造護套的碳纖維材料其抗拉強度一般超過3 500 MPa,所以護套在工作時也能滿足其強度要求.

    對比0 r/min和8 × 104 r/min轉速下兩種工況的應力可以看出:在轉子高速旋轉時接觸面上的最大接觸應力由23.2 MPa提升至27.9 MPa,接觸應力增大;最大徑向應力由0.0 MPa變?yōu)?.5 MPa,最大切向應力由-58.2 MPa變?yōu)?8.2 MPa,切向應力和徑向應力也同樣增大,切向應力由壓應力變?yōu)榱死瓚? 根據前文理論計算時所得,在高速旋轉時永磁體與其碳纖維護套間的過盈量較靜態(tài)時有所增加,這也是導致接觸面上接觸壓力增大的原因. 雖然接觸面上接觸壓力的增大有助于維持過盈配合的狀態(tài),即有利于使永磁體受壓應力的作用,但仍然無法與高速離心力引起的拉應力相抵消,最終永磁體在切向方向從原來的受壓轉變?yōu)槌惺芤欢ǖ睦瓚? 鑒于永磁體的抗拉強度很小,所以設置的過盈量既不宜過小也不宜過大,存在一個適中的取值范圍.

    一般來說,當轉子高速旋轉受到很大離心力時,接觸面應該有分離的趨勢,會引起接觸面上的接觸應力有所減小,但表2中的數據顯示永磁體與其碳纖維護套間的接觸應力反而有所增加,這是由于碳纖維材料的低密度使得護套所受的離心力較小,而高彈性模量又使其抗變形能力較強. 由此可以看出,護套材料對護套和永磁體之間的接觸有很大影響,進而對永磁體的強度造成影響. 除了護套的材料之外,在轉子結構不改變的情況下,轉子轉速、裝配過盈量以及溫升等參數影響還需要進一步分析.

    2.3 ? 永磁體強度的關鍵因素分析

    分別設置轉速為6 × 104 r/min,過盈量為0.05 mm,護套材料為合金鋼,重新進行仿真分析. 鋼護套和碳纖維護套的材料屬性如表3所示,仿真結果如表4所示.

    對比轉速為6 × 104 r/min與8 × 104 r/min兩種工況下的應力可以發(fā)現,永磁體內的最大切向應力、最大徑向應力,護套內的等效應力以及接觸面的接觸應力均隨轉速的降低而有所減小,有利于永磁體和護套的保護,使其遠離各材料的強度極限. 但考慮到空壓機的工作方式是電機驅動轉子高速旋轉帶動葉輪對外做功,將自然空氣壓縮成高壓空氣并提供給燃料電池系統(tǒng),所以,對于空壓機,轉速對外的輸出性能有著至關重要的影響. 由上可知,降低轉速,理論上的確是降低永磁體所受拉應力的有效方式,但在實際應用中會降低輸出功率,故不宜采用.

    在轉速為8 × 104 r/min,對比過盈量分別為0.03 mm和0.05 mm兩種情況下的應力可以看到,隨著過盈量的增加,接觸面上的最大接觸應力也隨之增加,由27.9 MPa增加到43.2 MPa,由此所能提供的預壓應力也就增加了. 在高速旋轉時,永磁體的最大切向應力由原來的48.2 MPa拉應力轉變?yōu)?.93 MPa壓應力,即在轉速為8 × 104 r/min時,永磁體在切向方向依然承受的是壓應力;而在徑向方向雖然永磁體的最大徑向應力基本保持不變,但永磁體在整體上承受了更大的壓應力,這對于保護永磁體來說是有利的. 護套的最大等效應力也隨著過盈量的增加而呈現顯著提升,從原來的395 MPa增長為591 MPa,提高約49.62%.

    雖然護套的最大等效應力依然沒有超過它的抗拉強度,但可以看出護套的最大等效應力會隨著過盈量的增加而迅速增大. 由此可知,雖然增加過盈量對保護永磁體免受拉應力的影響很有幫助,但考慮到過盈量的增加會同時增加護套的最大等效應力,而且太大的過盈量也會使裝配流程更為困難,所以一味地增加過盈量的值也是不可取的. 對于過盈量,仍然需要通過適當的計算來選取最合適的值,這樣才能使整個轉子最安全有效地工作.

    對比0 r/min與8 × 104 r/min兩種工況時的應力可以看出,在使用鋼護套且轉速為8 × 104 r/min時,接觸面上的接觸應力要小于0 r/min時的情況,即永磁體和護套間的過盈配合有分離的趨勢,造成這個結果的原因是過盈量的減小,其值是由旋轉轉子護套和永磁體在徑向位移量的差值所決定的. 在永磁體位移量相同的情況下,護套的位移量就會決定過盈量是減小還是增大,即接觸面間是有分離趨勢還是壓緊趨勢. 當護套的位移量大于永磁體的位移量時,接觸面就會有分離的趨勢,反之則有壓緊的趨勢.

    從護套位移量的計算式中可以看到,泊松比對位移量的影響相對較小,護套材料的密度越大、彈性模量越小,護套的徑向位移量也就越大. 鋼護套的彈性模量值與碳纖維護套的彈性模量值很接近,而密度卻遠大于碳纖維護套,因此在護套材料為鋼時,在離心力作用下的位移量更大,也就造成了接觸面有分離的趨勢,接觸面上的接觸應力也相應減小.

    此外,由于鋼護套和碳纖維護套的彈性模量相差不大,因此在使用鋼護套時,永磁體的預壓應力比使用碳纖維護套時只是略有減小. 對比轉速同為8×104 r/min時鋼護套和碳纖維護套下永磁體的應力可以看出,由于在使用鋼護套時接觸面間有分離的趨勢,加之轉子高速旋轉時產生的巨大離心力作用,才會造成永磁體的切向應力與使用碳纖維護套時的巨大差別;永磁體切向的拉應力大幅增加,甚至已經超過了永磁體的抗拉極限,導致永磁體出現強度破壞.

    由此可見,在同等條件下,只要改變護套的材料就會對高速旋轉時的永磁體內部應力帶來顯著的變化. 正因為護套材料的影響如此之大,所以在選擇護套材料后必須重新進行校核,以免無法在高速運轉時有效保護永磁體. 理論上,選擇密度更小、彈性模量更大的材料更有利于保護永磁體免受拉應力的影響,但在實際選擇護套材料時還需考慮實際的情況選擇最合適的護套.

    2.4 ? 溫升對永磁體強度的影響

    除了過盈量、轉速和護套材料這三項因素外,溫升也會對轉子強度產生影響[19]. ?轉子在高速運行一段時間后,機械和電磁損耗均會使轉子的溫度升高. 由于護套和永磁體材料熱膨脹系數不同,接觸面間的過盈量有所減小,從而導致接觸壓力的減小,造成護套的等效應力和永磁體的拉應力均相應增大,最終可能由壓應力轉變?yōu)槔瓚?,甚至拉應力的極值超過永磁體的抗拉強度. 因此,在設置過盈量時,需考慮是否適當增加過盈量的大小以抵消溫度上升對過盈量和永磁體強度帶來的影響. 釹鐵硼的熱膨脹系數為6.5 × 10-6 ℃-1,碳纖維的熱膨脹系數為8.8 × 10-6 ℃-1,鋼的熱膨脹系數為10.5 × 10-6 ℃-1 . 對2.1節(jié)中的模型施加不同的溫度載荷,并設置參考溫度為20 ℃,過盈量為0.03 mm,轉速為8 × 104 r/min,護套材料分別選用碳纖維和鋼. 部分仿真結果如圖7和圖8所示,數據匯總于表5中. 將表5中永磁體最大拉應力和溫升的關系進行數值擬合,所得關系曲線如圖9所示.

    由表5可知,隨著溫升的提高,護套和永磁體間的接觸應力逐漸降低,這表明了裝配過盈量會因熱變形而減小. 由于過盈配合對永磁體的保護作用削弱,永磁體所受的最大拉應力逐漸增大. 溫升100 ℃時,永磁體搭配碳纖維護套和鋼護套所受的最大拉應力相比室溫時分別提升了30.0%和39.1%. 尤其是鋼護套保護的永磁體在溫度升高20 ℃時就已超過了材料的強度極限. 這是由于鋼材料的熱膨脹系數高于碳纖維和永磁體,所以裝配過盈量隨溫升的減小會更顯著. 溫升從0 ℃升至100 ℃時,改用碳纖維護套使永磁體應力的降幅從42.1%增至45.9%.

    溫升所引起的裝配過盈量減小值同樣可使用彈性力學理論計算:

    Δδt = (αst - αpm)rΔT ? ? ? ?(1)

    式中:αst為護套的熱膨脹系數(℃-1);αpm為永磁體的熱膨脹系數(℃-1);r為配合面半徑(m);ΔT為溫度升高量(℃).

    由式(1)可以計算出鋼護套由溫升導致的過盈量減小量為7×10-3 mm,相比碳纖維護套增長約73.91%. 因此,在護套材料的選取上應盡量選擇熱膨脹系數較低且與永磁體相近的材料,這樣可以有效抑制溫升對過盈量的影響,從而提升轉子的最高工作溫度.

    從圖9的關系曲線可以看出,兩種護套材料的永磁體最大拉應力均隨溫升呈近似線性增大,兩者間的數值關系為:

    σcf,t max = 0.29ΔT + 75.48 ? ? ? (2)

    σst,t max = 0.13ΔT + 43.90 ? ? ? (3)

    式中:σcf,t max和σst,t max為搭配碳纖維護套和鋼護套時的永磁體最大拉應力(MPa).

    從式(2)和式(3)中ΔT的斜率可以看出,搭配鋼護套時拉應力極值隨溫度的增長率是碳纖維護套時的2.23倍,這表明在使用鋼護套時結構強度受溫升的影響更大. 通過此數值關系也可以推算出在安全系數為1.2時,安裝碳纖維護套的轉子最高工作溫度可達176.5 ℃,遠高于電機轉子的正常工作溫度. 結果表明碳纖維護套材料在高溫工況時具有顯著優(yōu)勢,不僅在工作溫度范圍內無需進一步增大過盈量來加強保護作用,而且受熱變形的影響也更小.

    3 ? 結束語

    在對驅動燃料電池車用空壓機所使用的高速永磁電機轉子分別進行解析法和有限元法強度計算后,可得出以下主要結論:

    1)過盈量的初值可通過護套與永磁體間的接觸壓力確定,接觸應力的最小值為永磁體抗拉強度與安全系數相除所得;過盈量的減小量與轉子轉速的平方成正比,且受護套材料密度和彈性模量的影響.

    2)解析計算方法和有限元模型仿真方法的計算結果十分相符,各應力的相對誤差小于2.5%. 這是由于解析計算方法忽略了軸向應變的影響,而在有限元分析中考慮了軸向應變的影響. 鑒于2種方法計算的相對誤差較小,故可認為解析計算方法具有足夠的計算精度.

    3)轉子轉速、過盈量大小以及護套材料等都是永磁體和護套結構強度的重要影響因素,降低轉速或增加過盈量都能有效減小永磁體所受拉應力. 分析結果表明,選用密度更小、彈性模量更大的護套材料可有助于抑制工作轉速時實際過盈量的減小. 但在實際選擇護套材料時還需綜合考慮加工工藝、制造成本和裝配方式等多方面的因素.

    4)溫升對過盈量的減小有顯著影響,會使接觸面壓力降低,永磁體所受拉應力的極值隨溫升近似線性增大. 因此需要對轉子的熱態(tài)工況進行計算以驗證是否應適當增加靜態(tài)過盈量以抵消熱變形的影響. 此外,選擇熱膨脹系數與永磁體相近的護套材料更有助于抑制溫升對轉子強度的削弱作用. 結果表明碳纖維護套在高溫工況時對永磁體的保護效果優(yōu)于鋼護套.

    參考文獻

    [1] ? ?衛(wèi)國愛. 車用PEMFC空氣供給系統(tǒng)建模及控制策略研究[D]. 武漢:武漢理工大學,2010:36—37.

    WEI G A. Study on model and control strategy in air supply system of PEMFC in vehicle[D]. Wuhan :Wuhan University of Technology,2010:36—37. (In Chinese)

    [2] ? ?溫泉,盛苗苗,董天哥. 操作參數對質子交換膜燃料電池性能的影響[J]. 電子技術,2018,47(10):30—33.

    WEN Q,SHENG M M,DONG T G. Effects of operating parameters on performance of proton exchange membrane fuel cell[J]. Electronic Technology,2018,47(10):30—33. (In Chinese)

    [3] ? ?吳震宇. 車載燃料電池用空氣壓縮機工作性能研究及能效分析[D]. 杭州:浙江大學,2007:5—14.

    WU Z Y. Performance research and power analysis of air compressor in fuel cell vehicle[D]. Hangzhou:Zhejiang University,2007:5—14. ?(In Chinese)

    [4] ? ?鮑鵬龍,章道彪,許思傳,等. 燃料電池車用空氣壓縮機發(fā)展現狀及趨勢[J]. 電源技術,2016,40(8):1731—1734.

    BAO P L,ZHANG D B,XU S C,et al. Development status and trend of air compressor in fuel cells vehicle[J]. Chinese Journal of Power Sources,2016,40(8):1731—1734. (In Chinese)

    [5] ? ?王繼強,王鳳翔,鮑文博,等. 高速永磁電機轉子設計與強度分析[J]. 中國電機工程學報,2005,25(15):140—145.

    WANG J Q,WANG F X,BAO W B,et al. Rotor design and strength analysis of high speed permanent magnet machine[J]. Proceedings of the CSEE,2005,25(15):140—145. (In Chinese)

    [6] ? ?李振平,占彥. 高速永磁同步電機的轉子結構強度分析研究[J]. 機電工程,2016,33(7):900—903.

    LI Z P,ZHAN Y. Rotor strength analysis of high speed permanent magnet machine[J]. Journal of Mechanical & Electrical Engineering,2016,33(7):900—903. (In Chinese)

    [7] ? ?王保俊,畢劉新,陳亮亮,等. 碳纖維綁扎表貼式高速永磁電機轉子強度分析[J]. 浙江大學學報(工學版),2013,47(12):2101—2110.

    WANG B J,BI L X,CHEN L L,et al. Strength analysis of a surface mounted high speed permanent magnetic machine rotor with carbon fiber bandage[J]. Journal of Zhejiang University(Engineering Science),2013,47(12):2101—2110. (In Chinese)

    [8] ? ?張鳳閣,杜光輝,王天煜,等. 高速永磁電機轉子不同保護措施的強度分析[J]. 中國電機工程學報,2013,33(S1):195—202.

    ZHANG F G,DU G H,WANG T Y,et al. Rotor strength analysis of high-speed permanent magnet under different protection measures[J]. Proceedings of the CSEE,2013,33(S1):195—202. (In Chinese)

    [9] ? ?張鳳閣,杜光輝,王天煜,等. 基于多物理場的高速永磁電機轉子護套研究[J]. 電機與控制學報,2014,18(6):15—21.

    ZHANG F G,DU G H,WANG T Y,et al. Rotor containment sleeve study of high-speed PM machine based on multi-physics fields[J]. Electric Machines and Control,2014,18(6):15—21. (In Chinese)

    [10] ?吳震宇,曲榮海,李健,等. 表貼式高速永磁電機多場耦合轉子設計[J]. 電機與控制學報,2016,20(2):98—103.

    WU Z Y,QU R H,LI J,et al. Multi-field coupling rotor design for surface-mounted high-speed permanent magnet machine[J]. Electric Machines and Control, 2016, 20(2): 98—103. (In Chinese)

    [11] ?張超,朱建國,韓雪巖. 高速表貼式永磁電機轉子強度分析[J].中國電機工程學報,2016,36(17):4719—4728.

    ZHANG C,ZHU J G,HAN X Y. Rotor strength analysis of high-speed surface mounted permanent magnet rotors[J] Proceedings of the CSEE,2016,36(17):4719—4728. (In Chinese)

    [12] ?張超,陳麗香,于慎波,等. 不同保護型式下的高速表貼式永磁轉子應力與溫升分析[J]. 電工技術學報,2019,34(9):1815—1824.

    ZHANG C,CHEN L X,YU S B,et al. Stress and temperature rise of high speed surface-mounted permanent magnet rotor with different protection types[J]. Transactions of China Electrotechnical Society,2019,34(9):1815—1824. (In Chinese)

    [13] ?馬振杰. 高速永磁電機機械特性分析[D]. 沈陽:沈陽工程學院,2016:10—18.

    MA Z J. The analysis of mechanical properties about high speed permanent magnet motor[D]. Shenyang:Shenyang Insitute of Engineering,2016:10—18. (In Chinese)

    [14] ?徐秉業(yè),劉信聲. 應用彈塑性力學[M]. 北京:清華大學出版社,1995:45—46.

    XU B Y,LIU X S. Applied elastic-plastic mechanic[M]. Beijing:Tsinghua University Press,1995:45—46. (In Chinese)

    [15] ?方程. 高速永磁電機轉子特性分析及結構優(yōu)化設計[D]. 沈陽:沈陽工業(yè)大學,2010:38—39.

    FANG C. Properties analysis and optimization design of rotor for high speed PM machine[D]. Shenyang:Shenyang University of Technology,2010:38—39. (In Chinese)

    [16] ?王繼強,馬俊強,李連春. 高速永磁電機轉子護套強度分析[J]. 沈陽化工學院學報,2009,23(1):46—49.

    WANG J Q,MA J Q,LI L C. Enclosure strength analysis for the rotor of a high speed permanent magnet machine[J]. Journal of Shenyang University of Chemical Technology,2009,23(1):46—49. (In Chinese)

    [17] ?丁鴻昌,肖林京,張華宇,等. 高速永磁電機轉子護套過盈配合量計算及應力分析[J]. 機械設計與研究,2011,27(5):95—98.

    DING H C,XIAO L J,ZHANG H Y,et al. Interference fit calculation and stress analysis for rotor sleeve of high-speed permanent magnet electric machine[J]. Machine Design and Research,2011,27(5):95—98. (In Chinese)

    [18] ?陳家新,鄭傳銀,白洋,等. 高速永磁電機護套配合過盈量的計算及強度校核[J]. 微特電機,2015,43(1):21—24.

    CHEN J X,ZHENG C Y,BAI Y,et al. Interference fit calculation and strength verification for rotor sleeve of permanent magnet motor operating on high speed[J]. Small & Special Electrical Machines,2015,43(1):21—24. (In Chinese)

    [19] ?陳亮亮,祝長生,王萌. 碳纖維護套高速永磁電機熱態(tài)轉子強度[J]. 浙江大學學報(工學版),2015,49(1):162—172.

    CHEN L L,ZHU C S,WANG M. Strength analysis for thermal carbon-fiber retaining rotor in high-speed permanent magnet machine[J]. Journal of Zhejiang University(Engineering Science),2015,49(1):162—172. (In Chinese)

    收稿日期:2020-10-05

    基金項目:上海市自然科學基金資助項目(19ZR1460000),Natural Science Foundation of Shanghai(19ZR1460000)

    作者簡介:張智明(1979—),男,遼寧遼陽人,同濟大學講師,博士

    ?通信聯系人,E-mail:zhangzm@#edu.cn

    猜你喜歡
    空壓機
    空壓機曲軸箱進水故障實例分析
    淺談車用電動空壓機在新能源汽車上的應用
    紡織廠空壓機預處理吸氣參數對其能耗的影響
    噴油螺桿式空壓機組電耗高原因分析及對策
    解決3#空壓機排氣效率不高的難題
    基于S7—1200PLC的鋼鐵企業(yè)空壓機集中管控系統(tǒng)
    煤礦空壓機余熱利用系統(tǒng)研究及應用
    空壓機系統(tǒng)運行優(yōu)化研究
    淺析空壓機節(jié)能技術及其應用
    醫(yī)院空氣正壓系統(tǒng)設備的檢修及維護
    男人舔女人下体高潮全视频| www.av在线官网国产| 国产成人a区在线观看| 性色avwww在线观看| 久久久精品大字幕| 男插女下体视频免费在线播放| 男女啪啪激烈高潮av片| 国产精品久久久久久精品电影小说 | 免费观看在线日韩| 九九在线视频观看精品| 国产成人午夜福利电影在线观看| videos熟女内射| 国产91av在线免费观看| 欧美成人免费av一区二区三区| 一本一本综合久久| 免费观看人在逋| 日韩欧美精品v在线| 欧美日韩综合久久久久久| 色尼玛亚洲综合影院| 韩国av在线不卡| av在线亚洲专区| 欧美3d第一页| 国产精品一二三区在线看| 三级毛片av免费| 国产视频内射| 天天躁日日操中文字幕| 乱人视频在线观看| 我要看日韩黄色一级片| 天堂√8在线中文| 亚洲av免费在线观看| 亚洲18禁久久av| 日韩三级伦理在线观看| 美女内射精品一级片tv| 午夜福利成人在线免费观看| 亚洲在线观看片| 久久久久久久久久成人| 亚洲精品日韩在线中文字幕| av又黄又爽大尺度在线免费看 | 汤姆久久久久久久影院中文字幕 | 久久久精品94久久精品| 色吧在线观看| 3wmmmm亚洲av在线观看| 在线免费观看的www视频| av又黄又爽大尺度在线免费看 | 免费黄色在线免费观看| 久久草成人影院| 午夜激情欧美在线| 狂野欧美白嫩少妇大欣赏| 欧美xxxx性猛交bbbb| 日本一二三区视频观看| 亚洲国产精品成人久久小说| 国产人妻一区二区三区在| 嫩草影院新地址| 国产伦一二天堂av在线观看| 国产69精品久久久久777片| 韩国av在线不卡| 特大巨黑吊av在线直播| 一卡2卡三卡四卡精品乱码亚洲| 一夜夜www| 国产亚洲午夜精品一区二区久久 | 欧美成人a在线观看| 国产免费视频播放在线视频 | 国产高清视频在线观看网站| 日韩av在线大香蕉| 中文字幕精品亚洲无线码一区| 国模一区二区三区四区视频| 老司机影院毛片| 午夜福利成人在线免费观看| 2021少妇久久久久久久久久久| 国产精品久久久久久av不卡| 一区二区三区四区激情视频| 91aial.com中文字幕在线观看| 亚洲精品乱码久久久久久按摩| 一级黄色大片毛片| av福利片在线观看| 三级男女做爰猛烈吃奶摸视频| 三级经典国产精品| .国产精品久久| 人体艺术视频欧美日本| 99视频精品全部免费 在线| 看免费成人av毛片| 亚洲精品456在线播放app| 日韩欧美精品v在线| 日韩人妻高清精品专区| 三级国产精品欧美在线观看| 中文字幕av在线有码专区| 国产欧美日韩精品一区二区| 国产在视频线在精品| 高清在线视频一区二区三区 | 午夜福利在线观看吧| 丰满少妇做爰视频| 禁无遮挡网站| 欧美极品一区二区三区四区| 少妇熟女aⅴ在线视频| 欧美一区二区精品小视频在线| 汤姆久久久久久久影院中文字幕 | 亚洲精品,欧美精品| 高清av免费在线| 2021少妇久久久久久久久久久| 男人狂女人下面高潮的视频| 高清视频免费观看一区二区 | 两性午夜刺激爽爽歪歪视频在线观看| 国产高潮美女av| 国产精品人妻久久久影院| 亚洲婷婷狠狠爱综合网| 色尼玛亚洲综合影院| 特级一级黄色大片| .国产精品久久| 国产精品国产三级专区第一集| 日韩一本色道免费dvd| 水蜜桃什么品种好| 亚洲丝袜综合中文字幕| 亚洲精品乱码久久久久久按摩| 最新中文字幕久久久久| 男女视频在线观看网站免费| 极品教师在线视频| 国产激情偷乱视频一区二区| 边亲边吃奶的免费视频| 国产黄a三级三级三级人| 99九九线精品视频在线观看视频| av女优亚洲男人天堂| 一级爰片在线观看| 91久久精品电影网| 亚洲av免费在线观看| 中文乱码字字幕精品一区二区三区 | 日本爱情动作片www.在线观看| 能在线免费观看的黄片| 久久精品熟女亚洲av麻豆精品 | 国产成人午夜福利电影在线观看| 边亲边吃奶的免费视频| 男女国产视频网站| 色视频www国产| 免费av观看视频| ponron亚洲| 久久久久久伊人网av| 嫩草影院入口| 午夜精品国产一区二区电影 | 日韩三级伦理在线观看| 亚洲欧美清纯卡通| 十八禁国产超污无遮挡网站| 午夜福利在线观看免费完整高清在| 日韩中字成人| 好男人视频免费观看在线| 麻豆精品久久久久久蜜桃| 欧美高清性xxxxhd video| 人妻系列 视频| 中文天堂在线官网| 亚洲三级黄色毛片| 国产伦在线观看视频一区| 国产色爽女视频免费观看| 久久这里有精品视频免费| 天天一区二区日本电影三级| 黄色日韩在线| 成年版毛片免费区| 国产精品一区二区性色av| 中文欧美无线码| 国产精品久久久久久精品电影| 免费观看a级毛片全部| 黑人高潮一二区| 亚洲国产精品成人综合色| 男插女下体视频免费在线播放| 国产乱人偷精品视频| 波野结衣二区三区在线| av免费观看日本| 我要搜黄色片| 亚洲国产精品国产精品| 午夜福利网站1000一区二区三区| 亚洲伊人久久精品综合 | 少妇人妻精品综合一区二区| 欧美+日韩+精品| 3wmmmm亚洲av在线观看| 欧美日韩在线观看h| 丝袜美腿在线中文| 亚洲人成网站在线观看播放| 99久久中文字幕三级久久日本| 嫩草影院精品99| 亚洲精品456在线播放app| 人体艺术视频欧美日本| 日本av手机在线免费观看| 国产 一区 欧美 日韩| 婷婷色麻豆天堂久久 | 免费一级毛片在线播放高清视频| 日本午夜av视频| 久久久久九九精品影院| 国产伦理片在线播放av一区| 色吧在线观看| 亚洲国产成人一精品久久久| 亚洲激情五月婷婷啪啪| 性色avwww在线观看| 国产精品人妻久久久久久| 男人舔女人下体高潮全视频| 国内精品宾馆在线| 看十八女毛片水多多多| 最近视频中文字幕2019在线8| 一卡2卡三卡四卡精品乱码亚洲| 亚洲电影在线观看av| 亚洲欧洲日产国产| 久久久a久久爽久久v久久| av在线观看视频网站免费| 能在线免费看毛片的网站| 亚洲精华国产精华液的使用体验| 国产精品久久久久久精品电影| 亚洲四区av| 久久精品熟女亚洲av麻豆精品 | 国产精品国产高清国产av| 午夜免费男女啪啪视频观看| 人妻制服诱惑在线中文字幕| 国产精品国产三级专区第一集| 五月玫瑰六月丁香| 日本爱情动作片www.在线观看| 插阴视频在线观看视频| 别揉我奶头 嗯啊视频| 免费无遮挡裸体视频| av国产免费在线观看| 一夜夜www| 女的被弄到高潮叫床怎么办| 日韩国内少妇激情av| 男人舔奶头视频| 亚洲国产欧美人成| 国产麻豆成人av免费视频| 天堂av国产一区二区熟女人妻| 麻豆久久精品国产亚洲av| 亚洲在久久综合| 边亲边吃奶的免费视频| 久久99热这里只频精品6学生 | av播播在线观看一区| 久久国内精品自在自线图片| 三级经典国产精品| 综合色丁香网| 成人午夜精彩视频在线观看| 村上凉子中文字幕在线| 午夜a级毛片| 99久久精品一区二区三区| 亚洲精华国产精华液的使用体验| 国语自产精品视频在线第100页| 日韩人妻高清精品专区| 亚洲国产欧洲综合997久久,| 久久午夜福利片| 婷婷色av中文字幕| 欧美激情国产日韩精品一区| 成人三级黄色视频| 久久久精品欧美日韩精品| 中文字幕精品亚洲无线码一区| 99久久精品一区二区三区| 亚洲真实伦在线观看| 激情 狠狠 欧美| 国产精品伦人一区二区| 日韩,欧美,国产一区二区三区 | 天天躁日日操中文字幕| 国产在视频线在精品| 永久网站在线| a级毛片免费高清观看在线播放| 亚洲人与动物交配视频| 亚洲国产日韩欧美精品在线观看| 亚洲av男天堂| 女的被弄到高潮叫床怎么办| 亚洲成人精品中文字幕电影| 亚洲欧美中文字幕日韩二区| 亚洲中文字幕日韩| 男女视频在线观看网站免费| 亚洲国产精品合色在线| 国产精品日韩av在线免费观看| 在线免费观看不下载黄p国产| 国产91av在线免费观看| 国产在视频线精品| videos熟女内射| 夜夜看夜夜爽夜夜摸| 婷婷色综合大香蕉| 三级男女做爰猛烈吃奶摸视频| 久久久成人免费电影| 久久99热这里只频精品6学生 | 免费人成在线观看视频色| 男女啪啪激烈高潮av片| 视频中文字幕在线观看| 国产精品永久免费网站| 国产精品一区二区在线观看99 | 人人妻人人澡欧美一区二区| 国产精品女同一区二区软件| 极品教师在线视频| 欧美区成人在线视频| 日本猛色少妇xxxxx猛交久久| 亚洲成人久久爱视频| 小蜜桃在线观看免费完整版高清| 久久久亚洲精品成人影院| 欧美成人一区二区免费高清观看| 两性午夜刺激爽爽歪歪视频在线观看| 91久久精品国产一区二区三区| 精品人妻一区二区三区麻豆| 亚洲欧洲国产日韩| 欧美极品一区二区三区四区| 精品人妻视频免费看| 乱系列少妇在线播放| 国产精品久久久久久久电影| 69av精品久久久久久| 99热这里只有精品一区| 18禁在线无遮挡免费观看视频| 99国产精品一区二区蜜桃av| 六月丁香七月| 国产色婷婷99| 三级国产精品欧美在线观看| 欧美性感艳星| 蜜桃亚洲精品一区二区三区| 日韩欧美在线乱码| 亚洲乱码一区二区免费版| 日日摸夜夜添夜夜爱| 最近最新中文字幕大全电影3| 亚洲成人av在线免费| 麻豆国产97在线/欧美| 桃色一区二区三区在线观看| 内射极品少妇av片p| 久久久久久伊人网av| 乱系列少妇在线播放| 国产成人freesex在线| 高清日韩中文字幕在线| 国产美女午夜福利| 成人高潮视频无遮挡免费网站| 亚洲美女视频黄频| 熟女人妻精品中文字幕| 国产精品一区二区三区四区免费观看| 国产精品国产三级国产av玫瑰| 老司机福利观看| 欧美日韩在线观看h| 久久久久性生活片| 欧美日本亚洲视频在线播放| 秋霞在线观看毛片| 久久欧美精品欧美久久欧美| 国产精品美女特级片免费视频播放器| 国产午夜福利久久久久久| 久久韩国三级中文字幕| 国产白丝娇喘喷水9色精品| 一夜夜www| 欧美日韩国产亚洲二区| 美女xxoo啪啪120秒动态图| 国产激情偷乱视频一区二区| 少妇的逼好多水| 丝袜喷水一区| 国产一区亚洲一区在线观看| 一级爰片在线观看| 又粗又硬又长又爽又黄的视频| 又爽又黄a免费视频| 久久这里只有精品中国| 99久久成人亚洲精品观看| 3wmmmm亚洲av在线观看| 内地一区二区视频在线| 亚洲国产精品专区欧美| 国产精品久久电影中文字幕| 国产一区二区在线观看日韩| 亚洲美女搞黄在线观看| 亚洲怡红院男人天堂| 亚洲精品乱码久久久v下载方式| 亚洲国产精品sss在线观看| 国产黄色视频一区二区在线观看 | 最后的刺客免费高清国语| 国产黄片美女视频| 少妇猛男粗大的猛烈进出视频 | 国产色婷婷99| 国产一区二区在线av高清观看| 成年av动漫网址| 在线观看一区二区三区| 亚洲欧美精品综合久久99| 午夜福利在线在线| 最后的刺客免费高清国语| 国产高清不卡午夜福利| 国产亚洲精品久久久com| 插逼视频在线观看| 老司机影院成人| 久久精品熟女亚洲av麻豆精品 | 美女xxoo啪啪120秒动态图| 国产精品不卡视频一区二区| 少妇裸体淫交视频免费看高清| kizo精华| 一夜夜www| 日本熟妇午夜| 成人午夜高清在线视频| 夫妻性生交免费视频一级片| 国产在视频线精品| 亚洲精品aⅴ在线观看| 国产真实乱freesex| 身体一侧抽搐| 国产淫片久久久久久久久| 亚洲精品aⅴ在线观看| 午夜爱爱视频在线播放| 日本五十路高清| 久久精品国产亚洲网站| 国产av在哪里看| 午夜福利在线观看吧| 在线播放无遮挡| 一个人看视频在线观看www免费| 精品久久久久久久人妻蜜臀av| 国产午夜福利久久久久久| 亚洲精品乱久久久久久| 最后的刺客免费高清国语| 精品久久久噜噜| 久久99蜜桃精品久久| 一边摸一边抽搐一进一小说| videos熟女内射| 久久国产乱子免费精品| 男人的好看免费观看在线视频| 人人妻人人看人人澡| 亚洲欧洲国产日韩| 欧美xxxx性猛交bbbb| 高清午夜精品一区二区三区| 在线观看美女被高潮喷水网站| 观看免费一级毛片| 日本与韩国留学比较| 亚洲国产色片| 精品酒店卫生间| 免费大片18禁| 国产欧美日韩精品一区二区| 我的女老师完整版在线观看| 女人被狂操c到高潮| 少妇高潮的动态图| 我要看日韩黄色一级片| 99热精品在线国产| 国内揄拍国产精品人妻在线| 日本黄色片子视频| 男女下面进入的视频免费午夜| 国内精品宾馆在线| 亚洲色图av天堂| 精品国产一区二区三区久久久樱花 | 视频中文字幕在线观看| 免费人成在线观看视频色| 亚洲伊人久久精品综合 | 国产午夜精品一二区理论片| 我的女老师完整版在线观看| 亚洲av不卡在线观看| 午夜日本视频在线| 国产精品熟女久久久久浪| 国产精品1区2区在线观看.| 我要看日韩黄色一级片| 久久精品熟女亚洲av麻豆精品 | 午夜日本视频在线| 亚洲欧美日韩高清专用| 中文亚洲av片在线观看爽| 午夜亚洲福利在线播放| 99热精品在线国产| 日日撸夜夜添| 国产精品永久免费网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 美女内射精品一级片tv| 1000部很黄的大片| 桃色一区二区三区在线观看| 少妇被粗大猛烈的视频| 又爽又黄a免费视频| 国产精品一区二区三区四区久久| 看片在线看免费视频| 欧美一区二区国产精品久久精品| 日韩欧美精品免费久久| 亚洲成色77777| 亚洲av免费高清在线观看| 亚洲精品乱码久久久久久按摩| 精品无人区乱码1区二区| 精品人妻偷拍中文字幕| 国产v大片淫在线免费观看| 亚洲精品,欧美精品| 国产91av在线免费观看| 欧美精品国产亚洲| 亚洲性久久影院| 男女视频在线观看网站免费| av在线蜜桃| 99热全是精品| 热99在线观看视频| 啦啦啦韩国在线观看视频| 亚洲欧洲日产国产| 少妇熟女aⅴ在线视频| 性插视频无遮挡在线免费观看| 在线观看66精品国产| 日本欧美国产在线视频| 2021天堂中文幕一二区在线观| 成人二区视频| 国产成人免费观看mmmm| 小蜜桃在线观看免费完整版高清| 美女xxoo啪啪120秒动态图| 日韩国内少妇激情av| 欧美zozozo另类| 黄色日韩在线| 中文乱码字字幕精品一区二区三区 | 日产精品乱码卡一卡2卡三| 成人亚洲精品av一区二区| 欧美日韩在线观看h| 午夜激情欧美在线| 久久久久久久久中文| 国产精品.久久久| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美日韩高清专用| 免费av不卡在线播放| 国产成人a∨麻豆精品| 两性午夜刺激爽爽歪歪视频在线观看| 草草在线视频免费看| 少妇熟女欧美另类| 国产精品综合久久久久久久免费| 久久99热6这里只有精品| 国产不卡一卡二| 亚洲av成人精品一区久久| 亚洲精品国产成人久久av| 色综合色国产| av女优亚洲男人天堂| 亚洲精品日韩av片在线观看| 亚洲天堂国产精品一区在线| 久久婷婷人人爽人人干人人爱| 亚洲在线自拍视频| 久久精品影院6| 久久久久久久久中文| 中文字幕熟女人妻在线| 三级男女做爰猛烈吃奶摸视频| 亚洲综合精品二区| 亚洲激情五月婷婷啪啪| 久久久色成人| 国产乱人偷精品视频| 亚洲av成人精品一区久久| 亚洲精品亚洲一区二区| 亚洲婷婷狠狠爱综合网| 黄色配什么色好看| av在线天堂中文字幕| 十八禁国产超污无遮挡网站| 久久久成人免费电影| 午夜精品国产一区二区电影 | 91aial.com中文字幕在线观看| 久久久久性生活片| 日产精品乱码卡一卡2卡三| 中文字幕精品亚洲无线码一区| 亚洲av日韩在线播放| 国产真实乱freesex| 最近2019中文字幕mv第一页| 亚洲激情五月婷婷啪啪| 精品无人区乱码1区二区| 99久久人妻综合| 国产精品av视频在线免费观看| 一本一本综合久久| 欧美成人精品欧美一级黄| 婷婷色综合大香蕉| 欧美丝袜亚洲另类| 汤姆久久久久久久影院中文字幕 | 国产真实伦视频高清在线观看| 欧美xxxx性猛交bbbb| 国产精品国产三级专区第一集| 三级国产精品片| 欧美日韩在线观看h| 免费黄网站久久成人精品| 精品久久久久久久久亚洲| 高清毛片免费看| 少妇的逼水好多| 亚洲内射少妇av| 秋霞伦理黄片| 一级毛片久久久久久久久女| 国产伦精品一区二区三区四那| 婷婷六月久久综合丁香| 乱系列少妇在线播放| 国产高潮美女av| 久久精品国产99精品国产亚洲性色| 3wmmmm亚洲av在线观看| 内地一区二区视频在线| 亚洲伊人久久精品综合 | 中国美白少妇内射xxxbb| АⅤ资源中文在线天堂| 最近最新中文字幕大全电影3| 亚洲怡红院男人天堂| 国产 一区 欧美 日韩| 国产高清视频在线观看网站| 欧美变态另类bdsm刘玥| 视频中文字幕在线观看| 黄色一级大片看看| 亚洲精品国产av成人精品| 国产高潮美女av| 99热6这里只有精品| 人体艺术视频欧美日本| 国产精品久久视频播放| 麻豆成人av视频| 亚洲国产精品成人综合色| 亚洲在久久综合| 日日干狠狠操夜夜爽| 久久鲁丝午夜福利片| 一区二区三区四区激情视频| 免费看光身美女| 亚洲高清免费不卡视频| 日本黄色视频三级网站网址| АⅤ资源中文在线天堂| 午夜久久久久精精品| 我的女老师完整版在线观看| 黄色欧美视频在线观看| 国产精品三级大全| 搞女人的毛片| 99国产精品一区二区蜜桃av| 波野结衣二区三区在线| 国产精品嫩草影院av在线观看| 在线a可以看的网站| 午夜精品在线福利| 22中文网久久字幕| 国产视频内射| 亚洲精品成人久久久久久| 精品人妻一区二区三区麻豆| 九色成人免费人妻av| 91精品一卡2卡3卡4卡| av视频在线观看入口| 免费看光身美女| av女优亚洲男人天堂| 国产精品久久久久久精品电影| 国产激情偷乱视频一区二区| 日韩av在线大香蕉| 日韩av在线免费看完整版不卡| 高清毛片免费看| 在线观看美女被高潮喷水网站| 亚洲电影在线观看av| 嘟嘟电影网在线观看| a级毛片免费高清观看在线播放| 在线a可以看的网站| 亚洲精品456在线播放app| 久久久久久久午夜电影| 精品不卡国产一区二区三区| 久久精品久久久久久久性| АⅤ资源中文在线天堂| 国产白丝娇喘喷水9色精品| 麻豆av噜噜一区二区三区| 麻豆国产97在线/欧美| 中文亚洲av片在线观看爽|