（机械工程专业论文）往复轴磁性液体密封机理及在车门密封中的应用研究.pdf

北方交通大学工程硕士专业学位论文 y 8 7 8 8 3 6 摘要 随着国民经济的发展，人民生活水平的不断提高，人们对出行的重要交通 工具铁路客车，要求越来越高，今年铁路实施了第五次大提速，明年还将实旌了 第六次大提速，大部分铁路线路的时速将提高到2 0 0 公里每小时，随着速度 的提高，对客车车门密封要求也越来越高，而车门活动踏板传动轴密封，一直没 有得到很好地解决，而对往复传动轴密封的最好办法就是采用磁性液体密封。 磁性液体是伴随着高科技发展而产生的一种新型的功能材料，在许多 方面具有其它材料难以实现的功效，尤其是在密封应用领域。往复轴动密 封是密封的主要形式之一，随着人们对往复轴动密封性能要求的不断提高， 现有材料难以满足需要，于是本世纪八十年代初期，科学家们提出了一种 全新概念的动密封技术往复轴磁性液体密封技术。往复轴磁性液体密 封技术无论在理论上还是实验中都处在探索阶段。本论文的主要任务就是 对磁性液体在往复轴密封方面应用而进行的基础性理论探索和实验研究，并 且尝试在铁路客车上进行应用解决车门密封问题。 在第一章中，作者在搜集了大量国内外资料和现有文献的基础上，对 磁性液体的各项物理性能和应用优缺点作了概括性的论述，同时还论述了 往复轴磁性液体密封的发展及原理，使读者对磁性液体和磁性液体往复轴 密封的概念、发展现状有一个清晰的认识。 在第二章和第三章中，作者在对磁性液体密封结构进行了深入的理论 分析的基础上，对现有的计算机密封磁场计算程序进行了改编，并且为程 序增添了许多新的功能，使程序运行更快，使用更方便，更易于维护。 在第四章中，作者着重讲解了磁性液体往复轴实验台的设计和要求。在 设计中，提出了一些新的结构和设计方案，并且附上了实物照片和结构 简图及系统原理图。 在第五章中，作者根据车门活动踏板传动轴的具体结构，在实验台上进行 了实验研究，记录了实验数据和实 验现象，并且对实验数据和实验现象进行了深入的分析，总结出了可能造 成往复轴磁性液体密封件密封失效的主要因素。 总之，作者在往复轴磁性液体密封的理论和实验上做了的这些工作， 是对今后往复轴磁性液体密封在铁路客车上的设计和应用有指导意义。 关键词：车门密封、往复轴密封、磁性液体密封、 磁场强度计算、 往复轴密封实验台 北方交通大学工程硕士专业学位论文 a b s t r a c t a i o n 窖w i t hm ed 州e l o p 呲毗o ft h en a t i o n 叠l 盹o o m y t h el i v i n gl e v e l i np 印p i e i n c r e a s e sc o n d n u o u s l y ，t h er e q u 鹳to f p p i e2 r em o 阼8 n dm o r eh 培ht oi m p o r t a n l p i l e u pr a i l m a dp a 鹳e n 2 e rc a r ，t h es p e e do fr a i l m a dp u t si n c n a s e f o rt h en f t ht i m et h i s y e 疆r ，t h es p e e dw m t i i mt oi n c 件a s ef o rt h es i x t ht i m e 既ty e a r ，t h es p e e do fm o s t r a i l r o a d sc i k u i tw mj n c m 叠s et o2 0 0l 【i l o m e t e r so fp e rh o u r s ，a l o n gw 油t h et h es p e e do f r a i l r o a dp u t si n c 仲叠s e ，t h er e q u e s to fs e a it ot h ep 叠s 辨n g e rc a rd o o rt m o 仲a n dm o m h i g h ，b u tt h ec a rd rm o v a b i ep e d a ls p 弛a d st 0m o v e 恤es t a s e a i s m p l e t e l y ，h 矗s i n g b e e nh a v i n gn o 协g e t t os o l 他n i c e i y b u ts p 他盛dt ot h eb a c k8 n df o r t ht h es t a l kt h a t i n o v et h eb e s tw 4 ys 嘲lc o m p i e t e l yi st o _ d o p tt h em a g n e t i s mi i q u i ds e a i sc o m p l e t e l y m a 鼬l e t i cn u 柚i sah e ws o r to f f l l n c t i o a lm a t e r i a i ，w h i c hi sc o m ei t ob e i n gw i t ht h e d e v e l o p m e n to ft h eh i g h t 髓h nh a sm a n yg df - l n c t i o st h 矗to t h e rm 8 t e r i 叠l sd on o t h a v e ，e s p e c i a y i n t h e n e l do fs e a t h e d y n a m i cs e a l o f t h er e c i p r o c a t i n gs h a f t i sa m a i n f o r mo f s e a l b u tw n ht h ej n c 件8 s j 2c a p a b m t i e so f t h ed y n a m i cs e a lo f t h er e c i p m c a t i n g s h a f tt h a tm a nm q u i n d t i i em a t e r i a li ne x i s t 蛆c ec o u m n tm 钟t s oa tt h eb e g i n i i i go f t h e1 9 8 0 s ，恤es c i e n l i s t sp u tl b n y a r dan e wc o n c e p to f d y n a i n i cs e 丑lt e c h 一t h em a g n e t i c n u 甜s e a lo f t h er 盹i p m c a t i i i 譬s h _ f t n o wb o 山t h et h 呻忡t i ca n dt h ee x p e r i m e n t so f a p p l y i n gm em a g n e t i cn u i d 叠1w i t ht h e 愀i p m c a t i n gs h a f t 叠r es 蝴la t t h ee x p l o r i n g s t a ge t h em a i nt a s ko ft h i sp a p e ri se x p l o r i n gt h ef u n d a m 蛐t 4 lt h e o 仲t i co fa p p l y i n gt h e d y n a m i cs e a lo f 蜘er e c i p m c 盘t i n 2s h a f ta n de x p e r i m e n ti t a n dt r yt og of o r w a r da n a p p i i c a t i o ni nt h er a i l m 叠dp a s s e n g e rc a r r e s o l v e ，t h ec a rd o o rs e a i sc o m p i e t e i yt h e p r o b i e m i n t h ef l r s tc h a p t er ，t h ea u t h o rr e s u m p t i v e l yd i s c u s s e st h ep h y s j c a l p e r f o r m a n c 鹤a n d 恤ee x c e i l e n c eo fa p p h c a t i o n so ft h em 叠g n e t i cn u ma f t e rh e h e c t e d l o t so ft h ei f h r m 鼻t i o ni nh o m ea n da b i 畸a d 1 nt h em e a n 廿m e 。h ea l s od i s c u s s e st h e d e v e l o p m e n ta n df u n c t i o no ff h em 曩i f 皿e t 量cn u i ds e a l - s on 睡r e 8 d e 璐c a nk n o ww h a tj s m a 2 n e t i cn u i da n dm _ l r e t i cn u i ds e a i 1 nt h es e dc h a p t e r 矗n dt l l et h i r dc h a p t e r t h ea u t h o rc o m p n e st h en n i t ee l e m e n t 地t h o dp r o g r a mt oc a l c u l a t et h em a g n e t i co n 协eb 鼬eo fa n a l y z i n gt h es t r u c t u r eo ft h e m a 暨n e t i cs e a i a dm a n yn e wf u n c t i o n sa na d d e dt ot h ep r 0 窖r a mt om a k ei tr u n n i n g m o r eq u i c u y ，u s e dm o 比c o n v e n i e n t i y 矗n dm a i n t a i n e dm o r ee a s i l y i nt h ef o u r t hc 蛔p t e bt h ea u t h o re m p h a s i z e so nt h ed e s i g no ft h ee x p e r i m e n t a ir i g o fm a 2 n e t j cn u 试s e a lo far 凹i p r o c a t i n 2s h a f t t h e r ea r em a n yn e ws t r u c t u r e sb e i n g b r o u 曲tf o r w a r d t h ea u t l i o ra l s o 窟i v e sm a n yp h o t oa n ds k e t c ho ft h o s en e ws t r u c t u 代s t or e a d e r si nt h j sc h a p t e l i nt h e 矗f t hc h a p t e r ai o to fe x p e r i i n e n t so fm a 窖n e t i cn u i ds e a lh a v eb e e nd o n eo n t l i er j 2 a f t e r 叠n a l y z i n gt h ed _ t a 日n dp h 蛳o i 眦胍o ft h o 辩e x p e r i m 啪缸，t h ea u t h o r c o n c i u d e st h em a i nl 啪s o n st h 鱼tm a i t h em a 2 n e t i cn u i ds e a u n gd i s a b i e i ns h o r t ，w h 叠tt h ea u t h o rd o e si nt h i sp 鱼p e rc a nd i r 盹tt h ed 姻i g na p p c a t i o no f m a g n e t i c 玎u 坩s 姻lo far e c i p r 仳a t i n gs h a f ti nt h er a i l m a dp a 鹳e n g e rc a ri nt h ef u t u r e k e vw o r d s ： t h ec a rd o o rs e a l s ， r e c i p r o c a t i n gs h a f t ，m a g n e t i cn u d m a g n e t i c n u i ds e a k a n d ， 1 a g n e t j cc a i c u l a t i o n 第一章前言 第1 章前言 1 1 磁性液体 磁性材料科学是材料科学的一个分支，和整个科学一样，磁性材料近 几年；也有很大的发展。磁性材料从晶体结构来看可分为晶态材料和非晶 态材料。晶态材料又可细分为单晶材料和多晶材料。从化学成分来看又分 为金属材料( 合金或金属间化合物) 与氧化物材料，从形体来看，又可分为 块状与膜状。但是所有这些材料都是固体材料，而磁性液体却是液态的磁 性材料，是将众多微小的铁磁性或亚铁磁性微粒高度弥散于液态载液中而 构成的一种具有很高稳定性的胶体溶液。微粒与载液通过界面活性剂浑成 这种磁液，即使在重力场、电场、磁场的作用下也能长期稳定地存在，不 产生沉淀与分离，因此具有很高的实用性。如图1 _ l 所示。i l 刮 图卜l 磁性液体的组成 磁性液体有许多优越的物理性质。由于磁性液体是一种两相混合物， 它的基本组成部分是液体载体和悬浮着的固相颗粒，因此在研究磁性液体 的物理性质时，总是假定它是一种均匀的两相混合物，和普通的两相流体 一样，其密度和粘度受两相问的混合物控制。但不同的地方在于磁性液体 的粘度还受外磁场的影响，由于磁性液体是一种可以受磁场控制的流体， 所以它的磁化性能特别重要。 1 1 1 磁性液体的物理特性： a 磁性液体的磁特性f 1 5 l 郎之万的顺磁性经典理论很好地描述了磁性液体的磁学性质。因为液 体中的铁磁颗粒的体积约为原子体积的1 0 5 倍，颗粒的磁性也就强多了。这 些颗粒显现“超顺磁性”，也就是说，它们具有高的饱和磁化强度，没有剩 磁，没有矫顽力。图1 2 显示了某些磁性液体的典型磁化曲线。 北方交通大学工程硕士专业学位论文 糍柏凝墒h k 0 e ) 图卜2 蒹些磁性i 瘦体的磁化曲线 假定它是一个无相互作用的球形颗粒系统，液体的磁化强度就可写 作： f = 中胁( c o t h 偿一1 a ) 式中的口= 咖3 胁日，2 4 i ，而中是饱和磁化强度为m s 的磁性液体 的磁性颗粒的体积比。 但e l m o r e 的测量结果表明，实验观察结果同郎之万曲线符合得不 好，非把颗粒尺寸的分布纳入考虑不可：c h a i l t r e u 等证明，运用对数正 态型体积分布可以得到很好的符合。 在撒去外加磁场后，磁性液体磁化强度弛豫的机理有两个。第一个 磁化强度弛豫的机理是液体中颗粒的转动，这个机理可由b r o w n i a l l 转动 的弥散时间tb 来说明其特性，tb 是由下式给出的： 1b = 3 咋t 叩丁 ( 1 1 ) 式中的v p l 是颗粒的体积，n 是液体的粘度。但是，在冷冻状念或 在聚合物包覆了颗粒的状态下，这一机理将会被有效地拘制住。这时的 弛豫机理就是颗粒内磁化向量克服磁各向异性而转动，涅耳首先设想了 这一过程，同这一过程相关的另一个弛豫时间tw 。 对于单轴各向异性材料，这一弛豫时间是由下式给出的： 1n = 1o 2e x p x ( 1 2 ) 式中的x = k v p l k t = nd 3 k 6 k o t 。to 由下式给出： 2 北方交通大学工程硕士专业学位论文 o _ m “8 肌噬( 1 3 ) 式中的r 是旋磁比，o 是阻尼因子( 1 0 之) 。这里的k 是磁各向异性常数。 在液体下，小颗粒的磁化强度发生涅耳弛豫的可能性要大一些，而 大颗粒的磁化强度发生b m w n i a i l 转动弛豫的可能性要更大一些。驰豫过 程过渡的直径d s 可由式( 1 1 ) 、式( 1 2 ) 相等得出。在固态下则只可能发生 涅耳弛豫。 b 磁性液体的粘度【1 5 】 ( 1 ) 零磁场下的磁性液体的粘度在没有磁场时磁性液体显 示宏观各向同性，因为颗粒的布朗运动将使得磁晶各向异性和形状各向 异性都平均为零。 若磁性液体内的固体颗粒不仅尺寸足够小，而且固相的体积分量d 也很小，这时可以假定颗粒间没有相互作用，且圆球体大小相同，在这 样的情况下，磁性液体的动力粘性系数可以按e i n s t e i n 公式计算： = 仉( 1 + 三妒) ( 1 4 ) 上 式中仇 是磁性液体基载液体的动力粘性系数。 对于体积分量驴不很小的情况，e i n s t e i n 公式不能使用，r o s e n s w e i g 提出一种修正的办法，他假定磁性液体的动力粘性系数和固相体积分量 成 二次式的关系，即 吼一 町。1 + 口矽+ 6 2 ( 1 5 ) 式中a 、b 两个系数可由稀疏相、及固相球形颗粒的最大密度堆彻来 确 定，a 、b 的值分别为a 一2 5 和b = 1 5 5 。 在磁性液体的磁性固相颗粒的表面上都附着一层分散剂，若颗粒的 坐 径是r p ，而分散剂层的厚度为6 ，则其有效半径增加到r p + 6 ，于是每个 北方交通大学_ 下程硕士专业学位论文 颗粒的体积从v p l 外增加到( 1 + 兰) 3 ，相应地其体积分量也就成为 0 ( 1 + 兰) 3 庐，结果上式就成为： 0 旦生：1 + 2 5 ( 1 + 生) 3 一1 5 5 ( 1 + 鱼) s z( 1 6 ) 彳1 cr pr p ( 2 ) 有磁场时磁性液体的粘度磁性液体在外加一个磁场后会 产生磁各向异性，m c t a 印e 使用毛细管粘度计测量过包有聚合物的钴颗 粒在甲苯中的稀悬浮液置于磁场下的粘度，实验表明，在磁场平行于剪 切面时的粘度变化量大约是磁场垂直于流动方向时粘度变化量的2 倍。 s h l i o m a s 对稀悬浮液中的磁粘效应作过理论分析，得到了磁化强度 的 横向分量和纵向分量的弛豫时间同磁场强度的关系式；影响粘度的弛豫 时间是t ，是垂直于磁场的磁化强度的弛豫时间。 当mh k o t 时， 在这种情况f ，叩，会取极限值 1 刁，= 三叩庐s i n2 口 ( 1 一1 1 ) z s h i o m i s 证明，平行于磁场方向的粘度的饱和值( n ，) 和垂直于磁 场方向的粘度的饱值( n ) 分别由下式给出： 11 n 一2 叩( 1 + 吾矿) ，叩- 2 叩( 1 十 ) ( 卜1 2 ) 式中的n 可由e i n s t e i n 公式给出，因而( n ，n ) ( n rn ) = 2 这个 4 北方交通大学工程硕士专业学位论文 式子和m c t a g u e 所得结果是一致的。 r o s c n s w c i g 等研究过市售的铁磁性悬浮液的粘度和施加磁场的关 系，量纲分析预言，在磁场中的液体粘度是流体动压力与磁压力的比值 l 叩。m 日的函数，式中l 是切变率，叩。 是载液的粘度，m 是磁性液 体的平均磁化强度，而h 是施加磁场，其实验结果如图1 3 。 l q ”m 1 1 0 6 翻卜：j 磁场和切变率对含f e ；0 。颗粒抟磁性液体粘滞艘孵影响 r o s e n s w e i g 等研究过高浓度悬浮液的粘度和磁场的关系，有证据表 明，当加上磁场时在浓悬浮液中有凝聚成簇团的现象发生，所有的切平 面受到一个有效的旋转磁场的作用，这将使情况更加复杂。 c 磁性液体的密度【1 5 】 磁性液体的密度可用下列表达式表述： p f = q 一扪p w ：七母p n p 旷赢篙 ( 1 1 3 ) ( 1 一1 4 ) p n 。 式中p ，为磁性液体的密度：p 。为基载体的真实密度；p 为固相 颗 粒的真实密度：庐为固相颗粒的磁性液体的体积分量；p 为固相颗粒在 磁性液体的重量分量。 北方交通大学工程硕士专业学位论文 d 磁性液体的光学特性 1 5 铁磁氧体在水中或在其它非金属性液体中的悬浮液系统的光学性质 已被人们研究过了，光学研究一般都要观察在磁场中或在梯度场中簇团 的形成，其中b i b i k 等将液体透明度的变化和凝聚过程按下式联系了起 来 t - t o c x p ( 一h )( 1 - 1 5 ) 式中的h 是与时间有关的吸收系数，而t 0 是在零磁场下的透明度。 g o l 曲e r g 等所做的进一步研究是关于放置在磁场中的磁性液体的偏 光性质，他们认为，偏振同由颗粒凝聚成链所产生的液体光学各向异性 有关。 h a y e s 也观察过针状链的形成，磁场消除后，链会解体。m a n i n e t 及 d a y i e s 和l i e w e l l y i l 观察过将铁氧体的非电导性磁性液体置于磁场中时 所发生的强烈的双折射现象和线性的双色性现象。在敝o “8 0 k a 血。) 的磁 场下，f e 3 0 4 颗粒在煤油中的悬浮液的同磁场方向平行的折射率与其同 磁场方向相垂直的折射率的差的典型值为1 2 1 0 ，在磁场不高时( 5 m m ) ，尽管极靴的密封级数不同，但它们的耐压能力趋近于相等。 密封级数和泄漏率的关系是，随着密封级数的增加，泄漏率下降。”1 n 。3 酬辟能力p 、 n = 2 、 m 叫j 、。一 r 一、羞 。o a 幔茸鸵 t w h s 州同行程卜学科蛾瓤nj j m 中 = 械体常列制l 蝌j = - ；_ 能，o i 影榔 1 3 4 理论分析的成果 直线运动的密封情况较为复杂，因而磁性液体动密封的动力分析比 较困难。磁性液体被拖甩出这一现象是瞬时的，同时还要需要考虑到永 磁铁产生的磁场强度的作用，因此还未建市合适的数学物理模型对直线 型密封作出正确的理论分析。实验导出拖甩力与粘膜状态的关系可以表 北方交通大学工程硕士专业学位论文 不如f ： t = u v d 式中一一拖甩力；u 一磁性液体粘度；v 一是轴速；d 一是轴的直径。 由此可见，随着轴速的增加，拖用力随着增大，磁性液体泄漏量增多。 1 9 9 0 年，前苏联学者s i e v s 矾g 等人对他们的往复轴磁性液体 密封的研究进行了系统地总结，指出往复轴磁性液体密封失效是由于 密封间隙中磁性液体的变形以及磁性液体被带走而引起的，这两个因素 是引起磁性液体临界耐压能力下降和引起密封微漏的主要原因。他们在 实验基础上提出了磁性液体带定量的公式： h = e x p ( f 。；) 万万，f r = 矿2 撕 其中，v 是磁性液体的动力粘度；a 是磁性液体在磁场中运动的加速 度；b s 为磁性液体的宽度：和l 是由实验确定的无量纲系数；其数学 模型如图1 9 所示。【7 一。i ：，_ 、。m ? | 一j ? ：。 。、。_ 、! _。 k l 瓷蕊i l 囊聂i 疆j 蕊瓢摹蠢漶耍 图l9 磁性液体变形和携带量研究模型 该研究还分析了磁性液体变形时耐压情况，数学模型如图1 1 0 所示 得出如下公式： m = o 朋0 ( 日1 一日2 ) p c r ：= 。m s ( 日。一日2 ) 北方交通大学工程硕士专业学位论文 m ：= o m s ( h l 一日c ) ， j 力露b 、? ，? i j 。 0 ，扯 i 影。r 、 6 pj 一 、 _w 一一 融 、。 dfd 。 8( 王) 轴怒静d ：的， 一妯向着压差的方i 匈运动 剁u 蓿与af ) 襁反的乃向远曲 吲卜】( ) 撩。心撕f 密封的位置和磁场强魁 的分n ， 1 9 9 6 年，李德才对密封间隙中磁性液体的流动状态进行了显微境厂 的微观观察，其微观实验台如图1 1 l 所示，发现即使往复轴运动的速度再 低，磁性液体膜与齿的接触点，磁性液体膜与往复轴的接触点都是变化 的，在一定的往复轴运动速度下、磁性液体膜的形状以及静止位景是一 定的，在不同的速度下限的变形形状以及静止位置是不同的，而且磁性 液体膜在往复轴速度作用下，静止到某一确定位置，不会因振幅的不同 而不同。实验表明，当磁性液体量比较充足且往复轴运动速度不大时， 磁性液体膜在往复运动速度为零后，仍可以回到往复轴静止位置。当往 复轴运动速度很高时，磁性液体脱离密封齿的束缚，而被带到与运动速 度方向一致的最前端齿区及其附近，而往复轴的速度不是常数时，磁性 液体膜随着往复轴速度的变化而发生振动。这使我们很可能只根据轴速 便可确定磁性液体密封性能和使用寿命，并选择实现最佳密封效果的轴 速。经初步观察和研究发现，磁性液体膜的形状与密封齿f 磁力线的形 状在离密封也较远时，并不完全复合，磁性液体量的变化也可以引起磁 性液体膜厚和形状的变化，然后在相应速度下决定出膜变形后的形状。 另外，磁性液体在往复轴上粘附的膜厚度很薄，与磁性液体密封间隙等 特征尺寸相比，在进行磁性液体动力学的理论分析时可以忽略不计。 1 7 北方交通大学工程硕士专业学位论文 、r ” r l ： 定位坎 授鞭。水磁 掘靴1 定位= l | # 往复轴 图卜1 1 观察磁性液体微观状态的实验台 在文献中，作者在建立上述微观实验的基础上。结合如图1 1 2 所示 数学模型图，得出了往复轴磁性液体的耐压公式： p h p 乎 h ( x 。) 一h ( x b ) o m 。+ 11 嘶高彳c 一面面而也 式中p i 一高压侧的压力；p a - 一大气压力；x ( r 点的坐标：x b _ b 点的 坐标： h ( x c 广x c 点的磁场强度：h ( x b 卜- x b 点的磁场强度；uo 一空气中的磁导率；m f 一磁性液体的饱和磁化强度：n 一磁性液体的动 力粘度；v 一往复轴的速度；h ( x c ) _ x c 点对应的磁性液体的膜厚：d 一 密封齿间往复轴的间隙：o 一磁性液体的表面张力； cr8 h 嘲1 】2 推导锰复式磁性液 奉密封耐压套t 的数举强曾j 图1 1 2 中实线为轴静止时磁性液体的位嚣。虚线为磁性液体在一定 轴速作用下的位置。由此公式即可确定密封耐压能力的情况，并定性地 判定往复轴运动对耐压能力的影响。【1 6 ，i：oo一荔 第二章往复轴密封件有限元磁场分析 第2 章往复轴密封件有限元磁场分析 确定磁性液体密封装置中的磁场强度( 主要是指装置中极靴和密封齿形 成的密封间隙处的磁场强度) 的大小和方向，是确定磁性液体密封装置最大 耐压能力理论值的前提条件。由此才能确定该密封装置在理论上能否满足密 封要求，并还能在此基础上确定在磁性液体密封装置中加入磁性液体量与密 封最大理论耐压能力的关系。另外，确定了密封件内密封间隙处的磁力线分 布，就可以分析往复轴静止或运动时磁性液体膜形态变化或磁性液体膜的运 动状态与磁力线分布的关系。 磁场强度的求解方法一般有图解法、解析法和模拟法。这些万法曾有 效地解决了不少问题，但它们各自存在着许多不足之处，如：图解法要花 费大量的时间，且精度不高：解析法不能求解具有任意边界形状的实际问 题；模拟法则设备制作比较费事，精度也不能保证；此外，对于解决非线 性磁场或边界介质复杂磁场的问题，这些方法都很难奏效。 1 3 磁场测量方法对于确定一般构件内磁场强度是比较有效的，但这种方 法对被测量对象的空间尺寸有一定的要求，不能对密封间隙为0 2 m m 以下 的磁性液体密封件内的磁场强度进行有效测量。 随着电子计算机的发展，数值计算方法得到了广泛的应用。对于电磁 场问题的计算，尤其是非线性电磁场和边界复杂电磁场问题的计算，数值 计算方法的效果比较好。求解电磁场常用的数值解法有有限差分法和有限 元法，这两种方法都是基于网格划分的，但有限元法的网格剖分比差分法 有更大的灵活性，能较好地适应边界曲折多变的形状，因而，再一般情况 下，有限元法有较高的精度。本次研究根据本实验装置的自身结构特点， 从电磁场理论的最基本方程麦克斯韦方程组出发，并作出适当合理的 近似处理，利用矢量磁位的意义，采用有限元计算的方法来确定磁性液体 密封装胃中的磁场强度。 2 1 磁场分析的理论基础一麦克斯韦方程组 麦克斯韦方程组是电磁场理论的基本方程，是研究电磁场的理论基 础和进行电磁场数值分析的计算基础，也是本次有限单元计算的数学基础， 北方交通人学工程硕士专业学位论文 它的微分形式如下： r d f 万：了+ 丝( 2 1 ) 西 。f 罾：旦曼 西 d 如d = d ( 2 2 ) ( 2 3 ) 硪v 否：o( 2 4 ) e 述麦克斯韦方程组的微分形式中，各矢量符号分别表示的物理意义 如下： 耳：磁场强眯培米) 了：传导电流密度( 安培米2 ) 五：电位移( 库仑米2 ) e ：电场强度( 伏特米) 口：磁通密度( 特斯拉) p ：自由电荷的体密度( 库仑米3 ) 麦克斯韦方程组是麦克斯韦在提出位移电流的假设下，全面总结电生 磁和磁生电现象后提出来的。方程组的意义如下： 方程f 2 1 ) 表示传导电流和变化电场是磁场的“旋涡源”，它表示变化电 场能够产生磁场这个重要的客观现象： 方程( 2 2 ) 表示变化磁场产生电场这个重要的事实； 方程佗一3 1 表示电场是有通量源的场，即电场是由电荷产生的： 方程( 2 4 ) 表示磁场无“磁通量源”，即磁场不可能由磁荷产生。 麦克斯韦方程是宏观电磁现象的基本规律，能适用于恒定电磁场、似 稳定电磁场和高频交变电磁场等不同情况。但在具体计算时，将逃到许多 北方交通大学工程硕士专业学位论文 问题。其中之一就是方程式少而求解的未知量多。因为方程组中共有四个 方程式但却有8 个标量方程1 5 个独立的矢量分量，并且上述的8 个标量方 程并非全部独立，其中两个散度方程可以分别从两个旋度方程及电流连续 方程： 讲v 了：一望 a ( 2 5 ) 推导出来。 ”1 为了增加方程的数目，现引入媒质的构成关系式，这也是表 征在电磁场作用下媒质的宏观电磁特性的媒质方程： d = 占e b = u h ( 2 - 6 ) ( 2 7 ) ，= 盯e ( 2 8 ) 其中： e ：介电常数( 法米2 ) ； u ：磁导率( 亨，米) ； o ：电导率( 1 ，欧姆) 。 应当注意的是，只有在线性并且各向同性媒质的情况下，这三个特征 参数才是常数。对于非线性和各向异性的媒质，它们不是常数，随场强的 变化而变化。 麦克斯韦方程组的微分形式只适用于媒体的物理性质( 由u 、e 、o 等 参数来表征) 处处连续的空间。但是实际中经常要遇到媒质的性质在一个或 多个界面处有突变的情况。因而不同媒质的分界面( 如空气和铁的分界面) 处的场矢量也会有不连续的突变情况出现。所以对于这些交界面上各点来 说，麦克斯韦方程组的微分形式已失去意义而必须去考虑有限空间中场量 之间的关系，这种关系是由麦克斯韦方程组的积分形式所制约的。同时， 电磁场的边界条件即可由之而导出。麦克斯韦方程组的积分形式是： 伊如，= 肌；+ 詹蕊 牺撕：一盟：一旦f - d ； j a西 扣r d j = g = p 咖 ( 2 9 ) f 2 1 0 ) ( 2 1 1 ) 2 1 北方交通大学工程硕士专业学位论文 c r 否d ；= o 占 ( 2 1 2 ) 方程组中各标量符号分别表示的物理意义为： i ：传导电流( 安培) ； q ：区域内自由电荷总数( 库仑) p ：自由电荷体密度( 库仑米3 ) 巾：磁通量( 韦伯1 其中，在方程式( 2 9 ) 和( 2 1 0 ) 中，1 为电磁场中任意闭合曲线，s 为由l 所 界定的任意曲面；方程式( 2 1 1 ) 和( 2 1 2 ) 中s 为电磁场中任意闭合曲面，v 为 由s 所界定的体积。 方程式( 2 9 ) 反映了磁场与时变电场之间的关系，其中，表示s 面内 传导电流密度，= 兰表示s 面内位移电流密度。该式的物理意义是：在任意 班 电磁场内，磁场强度沿任意闭合曲线的线积分，等于穿过该闭合曲线所包 围面的全电流。 方程式( 2 1 0 ) 反映了电场与时变磁场之间的关系，其物理意义表示：在 任意电磁场内，电场强度沿任意闭合曲线的线积分，等于穿过该闭合曲线 所包围任意曲面的磁通量的时间变化率的负值。 方程式( 2 1 1 ) 反映了电场的基本性质：在任意电场内，穿过任意闭合蓝 面的电位移通量，等于闭合曲面所包围的自由电荷电量的代数和。 方程式( 2 1 2 ) 反映磁场的基本性质：磁感应线都是闭合的，在任意磁场 内，穿过任意闭合曲面的磁通量恒等于零。 综上所述，麦克斯韦方程组和媒质方程组组合在一起，就足以求解所 需要的未知量和求解边界条件。但从计算的观点来看，极少采用麦克斯韦 方程组去直接求解，因为方程中的待求量毕竟太多。如果以位函数来描写 场，待求量的数目可以减少，较为有利。 2 2 矢量位及其微分方程 在分析计算电磁场问题时，为了求出场量( e 或b ) 与场源( p 和j ) 之间 北方交通大学工程硕士专业学位论文 的关系，引用位函数( 或称势函数) 作为辅助量，可以减少未知数的个数，使 问题简化，有时也使物理概念更加清晰。在无旋场( 即旋度为零的场) 中可以 采用标量位函数，而在有旋场中，则必须采用矢量位函数。对于有电流的 区域，因为r o t 日：了，属于有旋场，因此必须引入矢量磁位j ，一般它是 空间坐标和时间的函数，包含三个空间分量。在国际单位制中，4 的单位 为韦米。 采用矢量磁位之后，应使麦克斯韦方程仍旧能满足，这只要丑和4 之 间满足 b = v 4 的关系，便能符合要求。因为d i v b = d i v “d t 4 ) 一0 ，即表示磁通连续的条 件永远满足。 对于没有电流存在的区域，矢量磁位爿同样能满足应用。 求解具体磁场问题时，根据电流分布决定矢量磁位4 ，再用( 2 - 1 3 ) 式求出曰，即得到解答。 为了从电流分布求爿，必须建立爿与电流之间的关系。 在磁性液体密封件中的磁场属于稳定磁场，因此进行磁场计算， 方程中不包含时变分量，故( 2 1 ) 式中，为常数，右端第二项为零。 将此式两边同乘以p 得：r o t 口= u _ ，。 将式( 2 1 3 ) 代入，得：r o t r o t 4 = uj 。 利用矢量恒等式m 廿o t 4 = 孕a d a i v 4 v 2 爿，再引入库仑条件v 爿= 0 北方交通大学工程硕士专业学位论文 则可以推导出下列矢量位函数微分方程： v 2 j = 一u 了 ( 2 1 4 ) 上式( 2 - 1 4 ) 称为线性泊松方程。 对于非线性介质，由于存在磁饱和效应，其磁导率u 不是常数， 而是磁通密度占的函数。如果再假设铁磁材料为各向同性，则根据静态 磁场基本方程、介质的构成关系式和向量磁位函数j 定义式( 2 1 3 ) 即可导出下列矢量位函数方程： v 。f 上v 。j 1 - 了 ( 2 1 5 ) l 方程( 2 1 5 ) 称为非线性泊松方程。 下面分析一下矢量磁位一的物理意义，我们把口对磁场内任意以 回路l 为边界的曲面s 积分，根据前面给出的斯托克斯方程和矢量磁位 彳的定义式( 2 1 3 ) ，得通过曲面s 的磁通量： 庐= p 正= f ( ，。f _ ) d ；= 伊斫 ( 2 - 1 6 ) 由旋度定义可得： d 爿d ， 一一 l i m 一= r d f 彳= 曰 m _ o 如 ( 2 1 7 ) 由此可以看出，矢量磁位j 沿任一闭合曲线的线积分等于通过以 此闭合曲线为边界的曲面的磁通量。当以该闭合曲线为边界的曲面限 定的面积趋于零时，其单位面积的线积分的极限值等于r o t j ，即等于磁 通密度否的值。 北方交通人学工程硕士专业学位论文 根据上述说明可知，由矢量磁位4 可以确定磁通密度b ，并进而 确定其它场量。在直角坐标系下，根据： b = 玩f + 占。，+ 酏七 ( 2 1 8 ) r o t j = ( 等一警卜+ ( 警一警) 丑 芸一等 _ c z 舶，i 砂玉jl 瑟缸。i 缸砂j 。 以及4 和b 的关系，可得口的各个分量与爿分量的关系： 耻等一警d v韶 b ，：娑一娑 。一翩，弛 曰：= 一孚 0 洲 ( 2 2 0 ) ( 2 - 2 1 ) ( 2 2 2 ) 针对本有限元计算问题的求解磁场为轴对称磁场，把上面直角坐 标系下的方程转化到圆柱坐标系下有： b = b ，n r + b o n e + b ：n = m t j = ( ；嚣一警卜，+ ( 警一等) ira 口出i 出西 9 + ；( 掣一澍五： 由此可得圆柱坐标下否的各个分量与j 分量的关系 ( 2 2 3 ) ( 2 2 4 ) 北方交通大学工程硕士专业学位论文 驴嚣一警r 口d 。 ( 2 2 5 ) ( 2 2 6 ) ( 2 - 2 7 ) 上述关系式( 2 - 2 5 ) 、( 2 2 6 ) 、( 2 2 7 ) 是本次磁场有限元计算的理论基础。 2 3 电磁矢量位的边界条件 在前面通过对场量和位函数的讨论，已经给出了各种状态下电磁场特 性的数学描述。然而，实际物理问题的数学模型构造还有另一方面，即定 解条件问题。 在前面已经给出了场量和位函数的方程，这些方程均属于一元二次线 性偏微分方程，这类偏微分方程定解问题的定解条件应包含给定待求场函 数u ( x ，y ，z ，t ) 的初始条件和边界条件，而对于本次计算的静态场函数u ( x ， v ，z ) 仅有边界条件。 所谓边界条件就是指与空间坐标( x ，y ，z ) 相联系，给出场域边界s 上 的待求场函数u ( x ，y ，z ) 的边值，通常有下列三种情况： i 给定的边界条件是整个场域边界s 上的场函数值： 此类边界条件称为第一类边界条件，与此类边界条件对应的边界问题 称为狄罩赫利问题； i i 给定的是场函数在整个场域边界s 上的法向导数值： 矗 詈i ；= 厂( 工，) ，z ) ( 2 2 9 ) 此类边界条件称为第二类边界条件，与此类边界条件对应的边界问题 称为聂依曼问题； i i i 给定的是场域边界s 上的场函数与其法向导数的线性组合： 丝加 一 丝加 盟， a。 坚西 丝良半 = = 瑟 北方交通火学工程硕士专业学位论文 i “+ g ( z ，y ，z ) 豢乳= ，( z ，y ，：) ( 2 3 0 ) lj 此类边界条件称为第三类边界条件，与此类边界条件对应的边界问题 称为洛平问题，通常把洛平问题化为前两种问题的和。 由前所述，本次有限元计算引入了磁矢量位，那么正确写出磁矢量位 在坐标系中不同媒质分界面上矢量位的边界条件，是利用矢量位解算电磁 场的前提条件。下面我们来推证用表示的边界条件的一般形式。 图2 1 j 的边界条件 如果我们紧贴两媒质的分界面取一个任意小回路，其高度h 为一个高 阶微量，如图2 1 所示，则依据斯托克斯定理，穿过此回路的通量为 妒= j 占d ；= l ( m f j ) d ；= j - d j ( 2 3 1 ) 因为回路高度h 为一高阶微量，故 s = h x v ，o = 盈= o 由此可得：a 1 t 刊k t ( 2 3 2 ) 即分界面两侧矢量位的切向分量相等，这是矢量的第一个边界条件。 应当注意，一般情况下j ，。j ：。，矢量磁位的法向分量并不相等，故 一i 42 。 北方交通大学工程硕士专业学位论文 其次，再沿上述任意回路计算磁场强度的积分回路，按照麦克斯韦方 程，有 f 西盈= 4 v 再- 玉= ( 了，+ 警) - d ； 在s - h d 1 一o 的条件下，如果界面上有面电流密度j 。，且对静磁场有 竽：o ，则上式可改写为 d r h 2 t h 1 t = j s ( 2 3 3 ) 式中，h 2 。与回路，方向相同取正，h n 与回路，的方向相反取负，回路， 的方向由所包围的电流j s 按右手定则确定。 如果用矢量来记述，可写为： h 2xh2 + h l 以l = ，； 式中：n 。：为媒质l 界面的外法线方向； h ：为媒质2 界面的外法线方向 西、；与了，间符合右手定则。 如果记；：；，：一二：与磁场夹角小于9 0 。的界面法线方向，并代入 万：! v 。j ，则式( 2 3 3 ) 可改写为 ；。 l ( v 。j ：) 一土( 甲。j ) 】：了。( 2 3 4 ) 2i 这就是矢量磁位的第二个边界条件。 北方交通大学工程硕士专业学位论文 2 4 磁场有限元分析 有了上面的电磁场理论基础，下面进行有限元分析。 2 4 1 有限元概论 有限元法( f i n i t ee l e m e n tm c m o d ) 这个名称是由c l o u 曲于1 9 6 0 年在其 著作中首先提出。三十多年来，以变分原理为基础建立起来的有限元法， 因其理论依据的普遍性，不仅广泛地被应用于各种结构工程，而且作为一 种声誉很高的数值分析方法已被普遍推广并成功地用来解决其他工程领域 中的问题，例如热传导、渗流、流体力学、空气动力学、土壤力学、机械 零件强度分析、电磁场工程问题等等。 传统的有限元法以变分原理为基础，即首先把所要求解的微分方程数 学模型边值问题转化为相应的变分向题，即泛函求极值问题：然后， 利用剖分插值，离散化变分问题为普通多元函数的极值问题，即最终归结 为一组多元代数方程组，解之即得待求边值问题的数值解。可以看出，有 限元法的核心在于剖分插值，它是将所研究的连续场分割为有限个单元， 然后用比较简单的插值函数