纳米尺度气体滑动轴承的DSMC模拟.pdfVIP

第 28卷第 5期 摩 擦 学 学 报 V01．28 No．5 2008年 9月 Tribology Sept．2008 纳米尺度气体滑动轴承的DSMC模拟 兰旭东，李志信 (清华大学 航天航空学院，北京 100084) 摘要：依据硬盘驱动系统中读写头与盘面的构造，结合现有的楔型和阶跃型气体滑动轴承的特点，提出了一种气体滑 动轴承结构 ，并采用直接模拟 MonteCarlo(DSMC)法数值研究了滑动轴承的最小间距、移动速度和系统温度对轴承承 载能力的影响．计算结果表明，所设计轴承的承载能力随移动速度的增大而增加，随最小间距的增加和系统温度的升 高而减小；轴承的承载能力对轴承最小间距的变化更敏感． 关键词 ：纳米尺度；气体滑动轴承；DSMC 中图分类号： 文献标志码 ：A 文章编号：1004-0595(2008)05-0411-05 随着计算机技术的飞速发展，硬盘存储密度不 一 种改进的读写头／盘面模型结构，并使用DSMC法 断提高，其读写头与盘面问距由几十纳米减小到十 数值研究了最小间距、读写头移动速度以及系统温 几纳米 ¨。J，此时必须考虑由于尺度减小而产生的稀 度对气体滑动轴承性能的影响，探讨了轴承承载能 薄气体效应．研究微纳尺度气体滑动轴承有利于改 力． 进硬盘设计、提高硬盘存储能力．硬盘读写头、盘面 1 直接模拟MonteCarlo(DSMC)法 及之间的气膜构成气体滑动轴承，可用 Reynolds方 程描述 J．常规尺度Reynolds方程基于无滑移边界 DSMC法 1【是一种基于分子运动和统计规律、 N．S方程导出，在微纳尺度下需要进行修正．许多学 用来模拟稀薄气体流动的数值方法．它通过计算成 者采用理论和实验 的方法研究 了气体滑动轴 千上万模拟粒子的运动和碰撞再现气体分子的运动 承 ，导出了修正的Reynolds方程，其无量纲形式 过程．在模拟过程中，每个模拟粒子代表大量的真实 为 气体分子，当计算时间步长小于物理碰撞时间时，分 [朋dP]=A (朋) (1) 子的运动与碰撞可以解耦，同时在粒子运动与碰撞 其中，X：x／L，P=p／p0， =h／h0，A= UoL／p0ho 时保持动量和能量守恒． 为轴承数，是动力粘度，h为轴承间距， 是无量 在模拟开始阶段，所有的粒子都被均匀的随机 纲Poiseuille流率．也有学者直接采用数值计算的方 分布到网格中，在每个时问步中，粒子首先以各 自的 法对气体滑动轴承进行研究 ¨”，但这些研究多数 速度平移一段距离，并完成与边界的可能碰撞，然后 采用楔型结构 (倾角 1。)．2002年，Fukui和 Ya— 粒子按一定顺序标号．在每个网格中，根据无时间计 mane 建立了阶跃型模型结构，并采用 DSMC法数 数器 (NTC)法选取 “碰撞对”，并计算粒子碰撞后的 值研究了滑动轴承的压力和剪切力．事实上，当尺度 速度分量．流动的宏观特性通过统计采样计算得到． 小到十几甚至几个纳米时，需考虑分子间的作用力 重复上述的步骤直到统计误差足够小． (范德华力)…，再考虑到表面摩擦力和压力作用及 2 模型结构的改进 温度的影响，读写头在盘面上运动时可能产生一定 形变．因此，有必要依据实际结构、构造新的模型结 楔型结构是根据硬盘驱动系统而