基于6DOF模型及动网格的动静压轴承刚度阻尼-动力学与控制学报.PDFVIP

第 １６卷第６期２０１８年１２月 动力 学 与控 制 学报 Ｖｏｌ．１６ Ｎｏ．６ １６７２⁃６５５３／ ２０１８／ １６⑹／ ４９０⁃６ ＪＯＵＲＮＡＬ ＯＦ ＤＹＮＡＭＩＣＳ ＡＮＤ ＣＯＮＴＲＯＬ Ｄｅｃ．２０１８ 基于６ＤＯＦ模型及动网格的动静压轴承刚度阻尼 数值计算∗ † 王攀　 刘保国　 冯伟　 赵耿 （河南工业大学 机电工程学院，郑州　 ４５０００１） 摘要　 针对传统雷诺方程求解三维油膜流场特性的局限性，提出基于６ＤＯＦ模型及动网格的动静压轴承刚 度阻尼计算方法．以具有典型结构形式的液体动静压轴承为例，通过加载６ＤＯＦ 自定义程序，采用非线性迭 代方法计算外载荷作用下轴心轨迹的瞬态变化过程，得到轴颈在外载荷作用下的静平衡位置；通过嵌入 ＵＤＦ宏程序以动网格更新方法实现对轴颈在静平衡位置的扰动，求解 Ｎａｖｉｅｒ⁃Ｓｔｏｋｅｓ方程得到轴颈扰动前后 位置变化后的瞬态油膜力，利用差分法求得动静压轴承油膜刚度和阻尼，并分析了不同转速下轴承刚度和 阻尼的变化规律． 关键词　 动静压轴承，　 ６ＤＯＦ，　 动网格，　 刚度，　 阻尼 ＤＯＩ：　 １０．６０５２／ １６７２⁃６５５３⁃２０１８⁃０６３ ［８］ 等 用ＣＦＤ程序对液体动静压轴承的静态和动态 引言 特性进行了计算，通过与标准润滑理论数值计算结 液体动静压轴承具有回转精度高、动态刚性和 果对比证明了该方法的有效性．为了反映油膜网格 ［９］ 阻尼减振性能好、使用寿命长等优点，在超高速精 在外界干扰下的扭曲和变形，熊万里 等提出基于 ［１］ 密磨削电主轴领域有着广泛的应用前景 ．电主轴 动网格模型的刚度阻尼计算方法，使油膜力计算结 系统的高速稳定性与轴承动力学特性紧密相关，因 果与高速时的工况更为吻合．采用动网格法求解轴 此对于轴承动态刚度和阻尼的精确计算成为研究 承刚度和阻尼时，首先需要准确计算轴颈在外载荷 的重点． 作用时的静平衡位置，以便对平衡位置施加位移扰 目前油膜轴承的数值计算模型主要有两种，一 动和速度扰动．本文基于６ＤＯＦ模型计算轴颈平衡 种是运用流体润滑理论并通过条件假设建立经典 位置，采用动网格更新方法实现轴颈在该位置的扰 ［２］ 的雷诺方程，如Ｒｏｗｅ等 基于经典雷诺方程的有 动，通过求解Ｎａｖｉｅｒ⁃Ｓｔｏｋｅｓ方程得到轴承刚度和阻 限差分法，用有限扰动法得到了非线性化的动静压 尼． ［３］ 轴承动力特性系数．贺玉岭等 采用类似的方法， １　 油膜刚度阻尼计算方法 建立了油膜在静平衡位置做微小扰动的计算模型， 进而得到动静压轴承的刚度和阻尼系数． 然而，该 采用 ６ＤＯＦ 模型及动网格的计算方法是在 方法不能很好描述扩散效应、挤压效应、惯性效应、 ＣＦＤ软件ＦＬＵＥＮＴ基础上，通过加载６ＤＯＦ