汪家瓊，姜國(guó)濤，王凱，錢文飛

(江蘇大學(xué)國(guó)家水泵及系統(tǒng)工程技術(shù)研究中心，江蘇 鎮(zhèn)江 212013)

磁力聯(lián)軸器[1]應(yīng)用于高速旋轉(zhuǎn)機(jī)械中時(shí)，對(duì)其結(jié)構(gòu)尺寸與質(zhì)量大小有嚴(yán)格要求，因此在滿足額定轉(zhuǎn)矩要求情況下，通常會(huì)采用減小隔離套及磁塊厚度與寬度的方式來縮小體積和重量.但在高速運(yùn)轉(zhuǎn)中，磁轉(zhuǎn)子與隔離套在多種因素影響下發(fā)熱變形，可能造成隔離套與磁轉(zhuǎn)子摩擦碰撞等危險(xiǎn)情況，使得磁力傳動(dòng)失效.因此對(duì)于高速運(yùn)轉(zhuǎn)的磁力聯(lián)軸器，基于多場(chǎng)耦合對(duì)隔離套這類薄壁件以及整體磁轉(zhuǎn)子，進(jìn)行強(qiáng)度分析顯得尤為重要.

近年來，多場(chǎng)耦合分析方法被廣泛應(yīng)用于電動(dòng)機(jī)和聯(lián)軸器的研究中.JIANG等[2]運(yùn)用ANSYS軟件分別對(duì)電動(dòng)機(jī)銅損、鐵損以及渦流損耗的影響因素進(jìn)行探究.XU等[3]運(yùn)用有限元方法對(duì)永磁耦合器簡(jiǎn)化的二維模型進(jìn)行計(jì)算，得到隔離套的渦流及溫度場(chǎng)分布.KONG等[4]模擬了不同轉(zhuǎn)速的磁力驅(qū)動(dòng)泵的內(nèi)流場(chǎng)，討論了不同轉(zhuǎn)速的磁力驅(qū)動(dòng)泵功率損失原因.ZHAO等[5]研究了磁熱單元順序耦合算法在永磁發(fā)電機(jī)溫度分析中的應(yīng)用.以往對(duì)熱磁耦合分析的研究多集中于渦流損耗，針對(duì)磁力聯(lián)軸器結(jié)構(gòu)強(qiáng)度與溫度分布的分析不足.

文中使用CFX對(duì)多轉(zhuǎn)速下磁力聯(lián)軸器的循環(huán)回路流場(chǎng)進(jìn)行計(jì)算，并結(jié)合Maxwell磁場(chǎng)計(jì)算結(jié)果，通過ANSYS Workbench平臺(tái)，對(duì)磁力聯(lián)軸器進(jìn)行磁-熱-流-固多場(chǎng)耦合分析，探究轉(zhuǎn)速對(duì)磁轉(zhuǎn)子體及隔離套的溫升與應(yīng)力、形變的影響.

1 結(jié)構(gòu)參數(shù)和理論基礎(chǔ)

1.1 結(jié)構(gòu)參數(shù)

現(xiàn)有高速“濕式”磁力傳動(dòng)測(cè)試試驗(yàn)臺(tái)[6]的結(jié)構(gòu)設(shè)計(jì)如圖1所示.圖中箭頭所示為冷卻介質(zhì)流通方向，可看出整個(gè)內(nèi)磁轉(zhuǎn)子體完全浸沒在冷卻介質(zhì)中，以保證冷卻液與轉(zhuǎn)子間能進(jìn)行充分熱交換.隔離套為靜止部件，內(nèi)表面與冷卻介質(zhì)充分接觸，一方面實(shí)現(xiàn)對(duì)隔離套的冷卻，另一方面保證介質(zhì)不外泄.磁力聯(lián)軸器主要由內(nèi)外磁轉(zhuǎn)子和隔離套組成，其樣機(jī)結(jié)構(gòu)如圖2所示，尺寸參數(shù)如表1所示.

圖1 試驗(yàn)臺(tái)結(jié)構(gòu)設(shè)計(jì)圖

圖2 圓筒型磁力聯(lián)軸器結(jié)構(gòu)設(shè)計(jì)圖

表1 磁力聯(lián)軸器主要結(jié)構(gòu)參數(shù)

該聯(lián)軸器磁極對(duì)數(shù)為16對(duì)，磁鋼材料為N38SH(釹鐵硼)，磁鋼表貼式嵌入導(dǎo)磁體中，充磁方式為徑向充磁.此外導(dǎo)磁體材料為Q235，隔離套為304材料，基體材料為2Cr13，且內(nèi)磁轉(zhuǎn)子由薄金屬包裹以防工作介質(zhì)對(duì)磁鋼產(chǎn)生銹蝕.N38SH的磁性能如下：退磁溫度大于等于150 ℃，剩磁Br為1.235 T，矯頑力HCB為900 kA/m，內(nèi)稟矯頑力HCJ大于等于1.595 kA/m，最大磁能積BHmax為300 kJ/m3.

1.2 多場(chǎng)耦合理論

首先對(duì)于磁生熱的物理過程，構(gòu)建數(shù)學(xué)模型建立控制方程確定每個(gè)變量間的耦合關(guān)系，將產(chǎn)生的焦耳熱通過Maxwell軟件與Workbench耦合接口映射到固體傳熱模塊，添加熱傳導(dǎo)以及對(duì)流換熱的邊界條件.基于磁熱耦合結(jié)果，順序加載到結(jié)構(gòu)力學(xué)模塊.根據(jù)材料的熱膨脹屬性運(yùn)用軟件模塊進(jìn)行熱力學(xué)形變計(jì)算，結(jié)果滿足控制方程[7].

熱應(yīng)力控制方程為

F=Eα(T-Tref)，

(1)

熱應(yīng)變控制方程為

ε=αT(T-Tref)，

(2)

上述式中：F為熱應(yīng)力，N/m2；ε為熱應(yīng)變；E為楊氏模量，GPa；α為熱膨脹系數(shù)，℃-1；Tref為相對(duì)參考溫度，℃.

其次以磁-熱-固耦合結(jié)果為基礎(chǔ)加入循環(huán)冷卻水體，針對(duì)以流場(chǎng)為主導(dǎo)的物理場(chǎng)展開對(duì)磁轉(zhuǎn)子各部件的溫度分布與結(jié)構(gòu)形變分析.根據(jù)質(zhì)量守恒方程、動(dòng)量守恒方程、能量守恒方程等得到耦合控制方程[8]為

(3)

式中：下標(biāo)f為流體；下標(biāo)s為固體；r為固體位移，m；τ為所受應(yīng)力，N/m2；q為熱流量，W；T為溫度，℃.

2 數(shù)值模擬

2.1 冷卻循環(huán)回路流場(chǎng)CFX數(shù)值模擬

2.1.1 模型建立及網(wǎng)格劃分

建立冷卻循環(huán)回路流場(chǎng)水體模型如圖3所示，水體與轉(zhuǎn)子結(jié)構(gòu)要完全貼合以保證多場(chǎng)耦合計(jì)算中載荷能準(zhǔn)確施加.其中冷卻循環(huán)回路中進(jìn)出口直徑均為12 mm，軸孔孔徑為6 mm.運(yùn)用ICEM軟件完成水體網(wǎng)格劃分如圖3b所示.其中在網(wǎng)格劃分過程中為保證網(wǎng)格質(zhì)量與計(jì)算結(jié)果精度，分別對(duì)隔離套、軸心回流孔等細(xì)小尺寸結(jié)構(gòu)處進(jìn)行面網(wǎng)格邊界層加密.

圖3 水體模型與網(wǎng)格劃分

2.1.2 邊界條件

由于RNGk-ε模型對(duì)剪切運(yùn)動(dòng)導(dǎo)致的湍流作用計(jì)算有優(yōu)勢(shì)，選用RNGk-ε模型[9].采用壓力進(jìn)口、質(zhì)量流量出口邊界條件.進(jìn)口壓力均設(shè)置為0.08 MPa，由增壓泵施加，出口流量均設(shè)置為1.67 kg/s.運(yùn)用CFX軟件計(jì)算不同轉(zhuǎn)速(3 000～9 000 r/min)下冷卻循環(huán)介質(zhì)在結(jié)構(gòu)中的流動(dòng)狀況[10].

2.2 磁力聯(lián)軸器Maxwell瞬態(tài)場(chǎng)數(shù)值模擬

2.2.1 模型建立及網(wǎng)格劃分

3D模型建立與網(wǎng)格劃分如圖4所示.模型網(wǎng)格為手動(dòng)劃分，渦流區(qū)的部件采用基于趨膚效應(yīng)的深入剖分法，并對(duì)模型交界處與細(xì)小處進(jìn)行加密.

圖4 3D模型與網(wǎng)格圖

2.2.2 邊界條件及控制參數(shù)

選擇能同時(shí)滿足散磁計(jì)算與無窮遠(yuǎn)處磁場(chǎng)計(jì)算的氣球邊界條件[11-12].此外對(duì)模型內(nèi)部繼續(xù)進(jìn)行邊界域細(xì)分，單獨(dú)將2個(gè)磁轉(zhuǎn)子用空氣層包裹，并定義旋轉(zhuǎn)域，其他區(qū)域定義為靜止域.瞬態(tài)磁場(chǎng)求解中激勵(lì)源的設(shè)置必不可少.在內(nèi)外磁轉(zhuǎn)子旋轉(zhuǎn)過程中，隔離套切割磁感線而產(chǎn)生感應(yīng)電流，將整個(gè)隔離套假想成一個(gè)單線圈導(dǎo)體，對(duì)隔離套施加一個(gè)初始電流源為0 A的激勵(lì).計(jì)算時(shí)長(zhǎng)t、計(jì)算步長(zhǎng)s分別由式(4)計(jì)算可得到.運(yùn)用Maxwell對(duì)不同轉(zhuǎn)速下的磁轉(zhuǎn)子進(jìn)行瞬態(tài)場(chǎng)數(shù)值模擬計(jì)算，計(jì)算公式為

(4)

式中：n為轉(zhuǎn)速，r/min；p為磁極對(duì)數(shù)，p=16.

2.3 磁-流-熱-固耦合數(shù)值計(jì)算

將不同轉(zhuǎn)速工況下Maxwell-3D的磁損結(jié)果與CFX的流場(chǎng)結(jié)果分別耦合至穩(wěn)態(tài)熱分析模塊與結(jié)構(gòu)靜力學(xué)分析模塊.對(duì)隔離套進(jìn)行溫度場(chǎng)分析、熱應(yīng)力形變分析以及磁轉(zhuǎn)子整體結(jié)構(gòu)強(qiáng)度分析[13-14].

2.3.1 模型建立及網(wǎng)格劃分

運(yùn)用CREO軟件構(gòu)建相應(yīng)尺寸的3D計(jì)算模型，保證磁轉(zhuǎn)子實(shí)體模型與CFX流場(chǎng)水體計(jì)算模型完全貼合.簡(jiǎn)化隔離套遠(yuǎn)端法蘭盤接口部分，將模型導(dǎo)入Workbench各計(jì)算模塊并設(shè)置各部件材料屬性，材料屬性如表2所示.表中：ρ為材料密度；γ為泊松比；Rd為拉伸屈服強(qiáng)度；Rp為壓縮屈服強(qiáng)度；h為熱傳導(dǎo)系數(shù)；Cp為熱容.在Mesh中對(duì)模型進(jìn)行網(wǎng)格劃分，并對(duì)結(jié)構(gòu)交界處與細(xì)小尺寸處進(jìn)行邊界層加密處理.模型網(wǎng)格如圖5所示.

表2 主要部件材料常溫下物理性能

圖5 隔離套與磁轉(zhuǎn)子體模型網(wǎng)格劃分

2.3.2 邊界條件及載荷

在穩(wěn)態(tài)熱分析模塊中，設(shè)定初始溫度為20 ℃，加載的熱源來自耦合的磁損結(jié)果，即隔離套的磁渦損耗會(huì)轉(zhuǎn)變成熱量，通過Maxwell軟件與ANSYS Workbench耦合接口將產(chǎn)生的焦耳熱映射到穩(wěn)態(tài)溫度場(chǎng)的模型中.考慮磁轉(zhuǎn)子在水冷時(shí)主要散熱方式為熱對(duì)流，因此選擇材料與流體介質(zhì)接觸面進(jìn)行換熱設(shè)置.換熱面主要包括隔離套內(nèi)壁、底部以及內(nèi)磁轉(zhuǎn)子包封套的外表面與兩端面.對(duì)流換熱系數(shù)可由努謝爾特準(zhǔn)則[15]進(jìn)行計(jì)算，計(jì)算公式為

(5)

式中:Nu為努謝爾特?cái)?shù)；D為當(dāng)量直徑，m；λf為流體的熱導(dǎo)率，W/(m·K)；Re為流體雷諾數(shù)；v為流速，m/s；μ為流體運(yùn)動(dòng)黏度，m2/s；Pr為普朗特?cái)?shù)；η為動(dòng)力黏度.

經(jīng)計(jì)算，冷卻水循環(huán)流量為1.67 kg/s時(shí)，以對(duì)流換熱系數(shù)12 000 W/(m2·K)進(jìn)行分析.

在靜力學(xué)結(jié)構(gòu)模塊計(jì)算中，分別定義整體轉(zhuǎn)子結(jié)構(gòu)上的重力、轉(zhuǎn)速以及旋轉(zhuǎn)軸與工作方向.對(duì)于載荷的施加，可分為2部分，溫度載荷：耦合數(shù)據(jù)源于穩(wěn)態(tài)熱模塊計(jì)算結(jié)果；壓力載荷：即磁場(chǎng)力與流場(chǎng)力，其耦合數(shù)據(jù)源于磁場(chǎng)與流場(chǎng)計(jì)算結(jié)果.

3 計(jì)算結(jié)果分析

3.1 隔離套穩(wěn)態(tài)熱分析

隔離套作為主要產(chǎn)熱部件，所產(chǎn)生形變大小是判定磁轉(zhuǎn)子能否正常運(yùn)轉(zhuǎn)的關(guān)鍵.對(duì)其單獨(dú)進(jìn)行穩(wěn)態(tài)熱分析，探究隔離套上的溫度場(chǎng)分布.圖6為不同轉(zhuǎn)速下隔離套上的溫度分布.

圖6 隔離套溫度分布云圖

從圖6中可以看出，高溫區(qū)域集中在筒壁中間，與磁鋼位置(50 mm)相對(duì)應(yīng)，因?yàn)檗D(zhuǎn)子由于磁耦合運(yùn)轉(zhuǎn)時(shí)，渦損主要產(chǎn)生于與磁鋼對(duì)應(yīng)長(zhǎng)度處，因此此段長(zhǎng)度內(nèi)筒壁所產(chǎn)生的渦損熱較大.

隨轉(zhuǎn)速的增加，隔離套的溫升會(huì)迅速增加，不同轉(zhuǎn)速下隔離套的最大溫升Tl如圖7所示.可以看出，在3 000 r/min時(shí)筒壁上最高溫度點(diǎn)為26.757 ℃，溫升為4.757 ℃.隨著轉(zhuǎn)速的增加，溫升迅速增加，在9 000 r/min時(shí)，最高溫度點(diǎn)為57.325 ℃，溫升達(dá)到35.325 ℃，且區(qū)域分布擴(kuò)大.高轉(zhuǎn)速工作下溫升幅度很大的主要原因是隨轉(zhuǎn)速的變大，磁渦損呈指數(shù)增加，因此溫升同樣增大.

3.2 隔離套與磁轉(zhuǎn)子結(jié)構(gòu)的應(yīng)力分析

對(duì)不同轉(zhuǎn)速工況下隔離套的形變Df與應(yīng)力F分布進(jìn)行分析，圖8為多場(chǎng)共同作用下的形變?cè)茍D與應(yīng)力分布云圖.從圖中可以看出，隔離套表面應(yīng)力主要集中在與磁鋼對(duì)應(yīng)長(zhǎng)度范圍內(nèi)，且此段長(zhǎng)度范圍內(nèi)形變量較大，形變大小均勻，呈向外膨脹狀態(tài).因?yàn)樵趯?duì)應(yīng)磁鋼范圍內(nèi)是磁渦損產(chǎn)生主要段，會(huì)產(chǎn)生較大熱應(yīng)力，影響套筒形變.另外隔離套是單獨(dú)的非旋轉(zhuǎn)體，只在頂端法蘭盤處存在固定約束，也易使套筒產(chǎn)生形變.而流場(chǎng)對(duì)隔離套形變產(chǎn)生的作用較小，冷卻循環(huán)系統(tǒng)進(jìn)口壓力為0.080 MPa，介質(zhì)循環(huán)中在經(jīng)冷卻回路壓降后，到達(dá)隔離套底部時(shí)的壓力已很小，從圖中即可看出，在隔離套底端的應(yīng)力大小是低于0.010 MPa的.

隨轉(zhuǎn)速增加，隔離套上的形變量與應(yīng)力值呈倍數(shù)增加.對(duì)比3 000 r/min轉(zhuǎn)速下形變(0.0036 mm)與應(yīng)力(5.060 MPa)大小，在9 000 r/min時(shí)增長(zhǎng)了7倍左右.圖9為不同轉(zhuǎn)速下隔離套的形變率δ，在9 000 r/min時(shí)最大形變率只有1.66%，相對(duì)于隔離套厚度1.5 mm，轉(zhuǎn)速對(duì)其形變率的影響較小.

對(duì)不同轉(zhuǎn)速工況下的磁轉(zhuǎn)子結(jié)構(gòu)體強(qiáng)度進(jìn)行分析，圖10為其結(jié)構(gòu)應(yīng)力分布云圖.從圖10中可以看出，應(yīng)力在磁轉(zhuǎn)子結(jié)構(gòu)上的分布較為規(guī)律，在周向上大小分布均勻，在軸向與徑向上遞減梯度明顯.在內(nèi)磁轉(zhuǎn)子體上，應(yīng)力主要集中在包封套表面，并呈梯度沿轉(zhuǎn)子內(nèi)部遞減，而且整個(gè)內(nèi)磁轉(zhuǎn)子形變量較小.對(duì)于外磁轉(zhuǎn)子體，應(yīng)力主要集中在磁鋼對(duì)應(yīng)位置與轉(zhuǎn)子內(nèi)表面上，由內(nèi)向外呈現(xiàn)遞減趨勢(shì)且有明顯分層現(xiàn)象，但結(jié)構(gòu)的形變量同樣不大.主要原因可能是因磁鋼被整個(gè)外導(dǎo)磁體包裹，在結(jié)構(gòu)上相對(duì)于隔離套而言更為穩(wěn)定，因此雖承受應(yīng)力較大，但產(chǎn)生的形變較小.

圖8 隔離套不同部位形變與應(yīng)力分布云圖

圖9 不同轉(zhuǎn)速下隔離套最大形變率

隨著轉(zhuǎn)速的增加，磁轉(zhuǎn)子體結(jié)構(gòu)所承受的應(yīng)力也越來越大.在3 000 r/min時(shí)應(yīng)力最大為6.35 MPa，但在轉(zhuǎn)速為9 000 r/min下應(yīng)力最大達(dá)到37.300 MPa.對(duì)比隔離套與磁轉(zhuǎn)子結(jié)構(gòu)上的應(yīng)力大小，如圖11所示.在不同轉(zhuǎn)速下，轉(zhuǎn)子上的應(yīng)力均大于隔離套上的應(yīng)力，這可能是因?yàn)橄啾雀綦x套所產(chǎn)生的熱應(yīng)力，磁鋼所產(chǎn)生的磁場(chǎng)力更大，因此在轉(zhuǎn)子結(jié)構(gòu)體上會(huì)出現(xiàn)更大應(yīng)力.

圖10 磁轉(zhuǎn)子結(jié)構(gòu)應(yīng)力分布云圖

圖11 隔離套與磁轉(zhuǎn)子結(jié)構(gòu)最大應(yīng)力對(duì)比圖

4 結(jié) 論

1) 分析在磁場(chǎng)與流場(chǎng)共同作用下轉(zhuǎn)速對(duì)隔離套溫度分布以及結(jié)構(gòu)形變的影響，結(jié)論如下：① 在隔離套圓筒壁上，高溫區(qū)域?qū)?yīng)于磁鋼長(zhǎng)度范圍并向圓筒壁兩端遞減，出現(xiàn)溫度分層現(xiàn)象.② 轉(zhuǎn)速對(duì)隔離套溫度分布的影響：隨轉(zhuǎn)速的增加，隔離套的溫升值呈線性增長(zhǎng).在3 000 r/min時(shí)溫升為4.757 ℃而當(dāng)轉(zhuǎn)速增加到9 000 r/min時(shí)溫升值達(dá)到35.325 ℃.③ 多場(chǎng)作用下對(duì)隔離套形變的影響：形變主要產(chǎn)生在磁鋼長(zhǎng)度范圍內(nèi)，與高溫集中區(qū)域一致.隔離套在周向上形變均勻，徑向上向外膨脹，軸向上在套筒兩端形變較小.磁熱效應(yīng)和特殊結(jié)構(gòu)是隔離套產(chǎn)生大形變的主要原因.④ 轉(zhuǎn)速對(duì)隔離套形變的影響：隨轉(zhuǎn)速的增加，圓筒壁最大形變量與應(yīng)力呈倍數(shù)增長(zhǎng)，對(duì)比3 000 r/min轉(zhuǎn)速下其形變(0.003 6 mm)與應(yīng)力(5.060 MPa)大小，在9 000 r/min時(shí)增長(zhǎng)了7倍.同時(shí)相對(duì)于隔離套厚度1.5 mm而言，在最大轉(zhuǎn)速9 000 r/min時(shí)其最大形變率只有1.66%左右，遠(yuǎn)小于隔離套間隙.此隔離套厚度滿足高速工作下的強(qiáng)度要求.

2) 對(duì)不同轉(zhuǎn)速下內(nèi)外磁轉(zhuǎn)子體進(jìn)行多場(chǎng)耦合分析，結(jié)論如下：① 應(yīng)力在磁轉(zhuǎn)子結(jié)構(gòu)上分布較規(guī)律，在周向上分布均勻，在軸向與徑向上遞減梯度明顯.內(nèi)磁轉(zhuǎn)子體上，應(yīng)力主要集中在包封套表面，產(chǎn)生形變較小.外轉(zhuǎn)子體上，應(yīng)力主要集中在與磁鋼交接面上，產(chǎn)生形變較小.② 轉(zhuǎn)速對(duì)磁轉(zhuǎn)子結(jié)構(gòu)的影響：隨轉(zhuǎn)速的增大，應(yīng)力與形變量均呈指數(shù)增加，在3 000 r/min時(shí)最大應(yīng)力為6.350 MPa，當(dāng)轉(zhuǎn)速增加到9 000 r/min時(shí)最大應(yīng)力達(dá)到37.300 MPa.而且由于內(nèi)外磁轉(zhuǎn)子在結(jié)構(gòu)上相對(duì)于隔離套而言更為穩(wěn)定，因此雖然承受應(yīng)力加大，但所產(chǎn)生的形變較小，故也是滿足運(yùn)轉(zhuǎn)時(shí)強(qiáng)度要求的.