永磁体在磁流体中悬浮平衡动态数值模拟.pdfVIP

第 16卷5期 应用基础与工程科学学报 V01．16．No．5 2008年l0月 JOURNALOFBASICSCIENCEANDENGINEERING October2o08 文章编号：1005-0930(2008)05-0679-07 中图分类号：TH701 文献标识码：A 永磁体在磁流体 中悬浮平衡动态数值模拟 刘桂雄，徐 晨，曹 东，蒲尧萍 (华南理工大学机械q-程学院，广东广州5~0640) 摘要：永磁体在磁流体中达到临界平衡状态时可以自悬浮．而以解析方法求解圆 柱形永磁体浸入磁流体中受到的磁性力难度很大．基于磁流体中磁场分布模型 及力学公式，利用有限元分析方法，采用新型移动边界控制方式，建立虚拟环境 和模型并进行磁场变化动态仿真．描绘永磁体位移与作用力关系曲线，精确求解 获得悬浮高度．试验结果表明，tkl0×20永磁体存在唯一悬浮平衡高度9．6mm， 数值计算数据与试验数据吻合，证明悬浮平衡位置的存在性和唯一性，同时验证 了动态数值分析结果真实有效性． 关键词 ：磁流体传感器；悬浮平衡；数值分析 磁流体是一种新型兼有磁性物体和流体特性的纳米功能材料．永磁体不施加任何外 力，可以自悬浮于磁流体中 ．永磁体浸入磁流体具有的这种近似真空中的零重力条件 下自悬浮能力，对于磁流体惯性传感机理研究具有巨大价值．基于该惯性传感原理研制的 新型传感器，在传感结构上可以避免直接的机械接触，不发生弹性变形，没有机械阻碍或 摩擦．永磁体在磁流体中运动下最终能否悬浮平衡，其受力临界状态是研究的关键．准确 描述临界状态和计算永磁体受力临界点，对建立永磁体在磁流体中的悬浮平衡理论有重 大意义 J．以磁流体动力学理论为基础，Naletova等 采用标量磁势法，研究球形永磁体 在充满磁流体的球形容器中静态悬浮问题．Vislovich等 假设给定永磁体一个微小位移， 计算出磁流体在方形容器中磁场分布，建立了静态悬浮力学模型．将解析式应用于实际受 力状态对整个过程求解 ，受函数复杂性制约非常难以精确求解受力临界点，目前在这方面 研究还没有太大进展．本文利用有限元分析方法，基于理想流体和磁体模型，将连续性力 学方程离散后逐点处理计算 引，采用数值模拟瞬态分析方法，动态描述出永磁体运动过 程磁场变化，从而研究永磁体在磁流体中的悬浮平衡问题． 1 磁流体中永磁体受力模型 对于浸没在磁流体 中的物体，根据磁流体中的应力张量，其受到磁流体的作用 收稿 日期：2007-07-．09；修订日期：2008-07-25 基金项 目：国家自然科学基金资助项 目；广东省自然科学基金资助项 目；高等学校博士学科 点专项科研基金(20060561003)；中国博士后科学基金 (2005037579)；华南理工大学 自然科学青年基金 项 目资助 作者简介：刘桂雄(1968一)，男，教授，博士生导师 应用基础与工程科学学报 力为 J： F= (p+ 。)n+HB小s (1) 式中，S为磁铁外表面； 为外表面S上的磁感应强度；rt为磁铁表面的外法线方向；日为 磁场强度，并且： H = +H,t (2) 式中， 为法向磁场强度； 为切向磁场强度． 设固定点的基准高度为h。、压强为P。，对于基准高度h的磁流体外表面层的任意点， P表达式 加]： P =一pgh+lXo(MH)+Po+pgho (3) 式中，P为磁流体密度；为平均磁化强度；。为真空磁导率． 将式(2)和式(3)代人式(1)，可得磁流体对浸没物作用力F表达式： F=[( f脚)n+H~B．tdS— V …