磁流变液流变特性的数值模拟分析.pdfVIP

助 锨 财 抖 磁流变液流变特性的数值模拟分析+ 李海涛1，彭向和1，陈伟民2 (1．重庆大学资源及环境科学学院，重庆400044；2．重庆大学光电工程学院，重庆400044) 摘 要： 建立了微观颗粒的动力学模型，提出了相 有人触及，粒子链在剪切运动中断裂后如何重新成链 应的算法，模拟了二维状态下颗粒成链及其剪切流动 以保证磁流变液宏观应力不致下降，成为本文的研究 过程，表明剪切过程中存在旧链断裂和新链形成，得 重点。 到了流动时磁流变液的剪切应力，并模拟了高剪切速 本研究通过分析颗粒所受的磁力、接触时的排斥 率下微结构的演变过程，发现新链形成速度小于旧链 力和运动时的粘性阻力，建立了颗粒运动的动力学模 断裂速度，出现了剪应力下降现象。 型，模拟了磁场作用下磁流变液链状微结构的形成以 关键词： 磁流变液；数值模拟；流变特性 及剪切运动下的流变过程，得到了不同剪切速率下磁 中图分类号：TB303 文献标识码：A 流变液的剪应力。 文章编号：1001-9731(2006)05—0710-03 2动力学模型 1 引 言 颗粒在载液中运动，既要与载液中其它粒子发生 磁流变液的奇妙特性与其微观机理密切相关。无 磁相互作用，还存在颗粒间的碰撞，并要受到载液粘性 外加磁场时，磁流变液的铁磁性颗粒呈无序状悬浮于 阻力的影响。因其粒径为微米量级，在外加磁场作用 载液中。一旦受磁场作用，颗粒逐渐沿磁场方向排列， 下，布朗热运动对颗粒的影响非常微弱，可以忽略不 最后达到一种稳定的链状有序结构[1’2]。这种微观结计‘6|。 构的演化导致磁流变效应的发生，并决定磁流变液的 2．1 磁 力 流变性能。 假设磁流变液中的固体铁磁性颗粒为圆球状，大 基于微观角度的数值模拟方法是研究磁流变液特 小一致，其半径均为Ot，载液不可磁化。颗粒被磁化后 产生的磁矩为[7]： 性的一条途径。S．Melle等模拟了旋转磁场下单链的 长度和形状，指出其受磁化率影响显著n]。王瑞金模 ^ m=VM=要7r口3YH (i) ’’ 拟了旋转磁场下不同Mason数对应的磁流变液聚合 J 链形态uj，J．Martin等研究了剪切过程中的剪切稀化式中，V是单个颗粒的体积，M—xH为磁化强度， y为磁化率，H为磁场强度。 现象¨J，H．Ly等利用Maxwell应力张量模拟了颗粒 外加磁场不强时，忽略颗粒被磁化后产生的磁场 在单向磁场下的成链状态[6]，然而，在这些动力学模拟 中，颗粒在磁流变液剪切流动时的动态演化过程很少 使周围的颗粒磁化，颗粒i受到的磁场作用力为【8]： cmj--5mirmjk7mir／7／J ㈤ Fm。蚤隘az08 r)ro+鹄．4 hm)] 其中肌为真空磁导率，rij表示两颗粒i和J之间