磁流变液微观结构数值模拟.pdfVIP

第39卷第3期 河 海 大 学 学 报 (自然 科 学 版 ) VoJ．39No．3 2011年 5月 JournalofHohaiUzfiversity(NaturalSciences) Mav2011 DOI：10．3876／j．issn．1000—1980．2011．03．014 磁流变液微观结构的数值模拟 李杰如 ，杜成斌2 (1．河海大学文天学院，安徽 马鞍山 243031；2．河海大学力学与材料学院，江苏 南京 210098) 摘要：通过分析磁流变液中固体颗粒的受力，利用分子动力学中的速度 Verlet算法，对磁流变液链 化过程和剪切过程进行了二维数值模拟，并对影响链化和剪切的主要因素进行了研究．研究结果表 明：粒子的运动最终会在模拟区域 内形成若干条几乎平行的链；颗粒链化时的变化过程先快后慢； 整个链化过程所耗费的时间在 2～3ins之间；颗粒体积百分率越大，模拟时所成通链越多；磁流变 液受到剪切荷载作用时，低剪切速率下，旧链 的断裂和新链 的生成基本上达到动态平衡的状态，高 剪切速率下，粒子向剪切方向的倾斜程度变大 ，但是通链减少，新链 的生成速度明显跟不上旧链的 断裂速度；单链的剪切强度先增大后减小；磁流变液的剪切屈服应力随磁感应强度和颗粒体积百分 率的增大而增大． 关键词 ：磁流变液 ；微观结构；链化；剪切 ；屈服应力 ；数值模拟 中图分类号：333381 文献标志码 ：A 文章编号：1000—1980(2011)03-031l-06 磁流变液是一种非常具有发展前途和工程应用价值的新兴智能材料，被认为是最具前途的智能材料之 一 … ． 磁流变液是 由非胶体的细小颗粒分散溶于绝缘基液 中形成的悬浮液，在外磁场作用下能产生明显的磁 流变效应，其力学、热学、流变学性能等会迅速发生变化，响应时间很短，为毫秒量级，且这种液态和固态之问 快速转化是可逆的．在该转化过程中，磁流变液的黏度保持连续、无级变化，固化强弱可 由外加磁场控制，能 耗极小，可实现实时的主动控制_2J． 1948年 RaJ)in0WJ最早发明了磁流变液及其应用装置离合器 ．Ginder等 7』认为在外磁场作用下，磁流变 液中的铁磁性固体颗粒形成的是一条一条的单链，并从磁性固体颗粒的非线性磁化、磁饱和以及相邻颗粒接 触位置处的局部磁场等方面进行了有限元计算，得出了磁流变液的剪切屈服应力的极值与磁性固体颗粒的 饱和磁化强度成正比的结论，该结论也被后来的实验所证实．Rosensweig~副设计出一种平均场模型——层模 型，由于该模型忽略了相邻颗粒问的磁场集 中效应，其结果与实验值相 比往往偏低．Bossis等 mJ将 Klingenberg等 】l¨利用电流变液推导出的恢复力公式应用于磁流变液，但由于没有考虑磁性粒子的非线性磁 化过程 ，其结果往往偏大 ． 近年来，国内武汉理工大学的阮中尉[z2J建立了BCT结构(体心四方结构)模型和物理材料模型，得到了 磁流变液受小角度剪切变形情况下其剪切应力的理论计算公式；中国矿业大学的赵四海等L13j研究了磁流变 液中颗粒的相互作用，分析了磁流变液产生链状结构的原因，推导出磁流变液在外磁场作用下的屈服应力公 式；广东工业大学陈伟俊 ¨J在对各种磁流变液的力学模型对比分析的基础上，根据磁场的库仑定律，提出了 一 种能够简洁表示各种因素对磁流变液力学性能影响的力学模型． _l 基于微结构的建模及力学分析只能描述磁流变液的宏观特性，不能描述微结构的形成机制及其演化过 程．影响磁流变液性能的因素很多，单纯从实验或宏观角度很难系统分析各因素的影响及效果，而基于微观 角度的数值模拟方法可以有效弥补这些缺陷，有利于优化磁流变材料的性能；设计高性能磁流变液．本文从 数值模拟角度出发，利用分子动力学中的速度 Verlet算法，基于磁流变液中磁性颗粒的实际受力，建立磁流 变液链化和剪切的分析模型，采用Matlab语言编程，分析粒子的运动规律，探讨磁流变液微观结构的链化形 成过程和剪切变形过程 ． ’ 收稿 日期：2010—06—23 基金项