基于改进的传递矩阵法识别电主轴系统滑动轴承油膜特性系数
|
摘要
针对高速电主轴系统在重载下滑动轴承油腔复杂的情况,以及油膜特性系数不便于直接获取的难题,提出了基于改进的传递矩阵法识别滑动轴承油膜特性系数的参数识别方法。该方法将参数识别问题转化为最优化问题,通过最小化电主轴系统不平衡响应的实验值和计算值的偏差来识别滑动轴承油膜特性系数。针对建立电主轴系统动力学模型时要考虑电机转子的自重和电磁转矩的问题,提出改进Riccati传递矩阵法,建立滑动轴承油膜特性系数和转子不平衡量参数识别中的流固耦合动力学模型,引入电机的自重和电磁转矩,扩大状态向量维数为十维,获得计算电主轴系统的不平衡响应的位移传递矩阵模型。并采用隔代映射遗传算法来寻优。最后通过数值算例对这一方法的效率和精度进行了验证。
Considering the bearing dynamic parameters cannot be gained by calculating directly when the oil cavity of sliding bearing is complex under the condition of high speed and heavy loads,a parameter identification method is presented to identify the bearing dynamic parameters of a sliding bearing-rotor system of high speed spindle based on the improved transfer matrix.In this method,the parameter identification problem is formulated as an optimization problem.The bearing dynamic parameters have been characterized through minimizing the error squared of the motorized spindle system unbalance response between the experiment results and the computational ones.In order to consider the weight and electromagnetic torque of the motor in building the fluid-solid coupling dynamics model of the motorized spindle system,the Riccati transfer matrix is improved.The status vector of typical unit is modified to ten for the introduction of the weight and torque of the motor.The displacement Riccati transfer matrix models are built to calculate the unbalance response of the motorized spindle system.The Intergeneration Projection Genetic Algorithm is used to minimize the error squared.The unbalance response in the motorized spindle system is gained by the improved transfer matrix method.Finally,the efficiency and robustness of the bearing dynamic parameters identification method are verified.
Considering the bearing dynamic parameters cannot be gained by calculating directly when the oil cavity of sliding bearing is complex under the condition of high speed and heavy loads,a parameter identification method is presented to identify the bearing dynamic parameters of a sliding bearing-rotor system of high speed spindle based on the improved transfer matrix.In this method,the parameter identification problem is formulated as an optimization problem.The bearing dynamic parameters have been characterized through minimizing the error squared of the motorized spindle system unbalance response between the experiment results and the computational ones.In order to consider the weight and electromagnetic torque of the motor in building the fluid-solid coupling dynamics model of the motorized spindle system,the Riccati transfer matrix is improved.The status vector of typical unit is modified to ten for the introduction of the weight and torque of the motor.The displacement Riccati transfer matrix models are built to calculate the unbalance response of the motorized spindle system.The Intergeneration Projection Genetic Algorithm is used to minimize the error squared.The unbalance response in the motorized spindle system is gained by the improved transfer matrix method.Finally,the efficiency and robustness of the bearing dynamic parameters identification method are verified.
引文
[1]陈双桥,蒋学武,吴新跃.采用边界元法测试轴承的动态刚度和阻尼[J].海军工程大学学报,2008,3(20):66-70.(CHEN Shuangqiao,JIANG Xuewu,WU Xinyue.Application of boundary element method to measuring dynamic stiffness and damping of bearing[J].Journal of Naval University of Engineering,2008,3(20):66-70(in Chinese)).
[2]DING Y,ZHAO B Y,WU B,et al.A condition assessment method for time-variant structures with incomplete measurements[J].Mechanical Systems and Signal Processing,2015,58:228-244.
[3]FRANULOVIC M,BASAN R,PREBIL I.Genetic algorithm in material model parameters’identification for low-cycle fatigue[J].Computational Materials Science,2009,45(2):505-510.
[4]AGUIR H,BELHADJSALAH H,HAMBLI R.Parameter identification of an elasto-plastic behavior using artificial neural networks-genetic algorithm method[J].Materials and Design,2011,32(1):48-53.
[5]LIEW A,FENG N S,HAHN E J.A non-linear transfer matrix technique for statically indeterminate rotor bearing systems[J].Proceedings of the Institution of Mechanical Engineers,Part C:Journal of Mechanical Engineering Science,2001,215(11):1343-1355.
[6]STANDER N,ROUX W,GEOL T,et al.Ls-opt user's manual[M].California:Livermore Software Technology Corporation,2009.
[7]伍乾坤,韩旭,胡德安.一种陶瓷脆性材料动态本构参数的计算反求方法[J].固体力学学报,2009,30(3):280-285.(WU Qiankun,HAN Xu,HU Dean.An inverse technique for identification of dynamic constitutive parameters of a ceramic brittle material[J].Chinese Journal of Solid Mechanics,2009,30(3):280-285(in Chinese)).
[8]Thakur M.A new genetic algorithm for global optimization of multimodal continuous functions[J].Journal of Computational Science,2014,5(2):298-311.
[9]闻邦椿,顾家柳,夏松波,等.高等转子动力学[M].北京:机械工业出版社,2000.(WEN Bangchun,GU Guliu,XIA Songbo,et al.Higher rotor dynamics[M].Beijing:Mechanical Industry Press,2000(in Chinese)).
[10]LUND J W.The stability of an elastic rotor in journal bearing with flexible damped supports[J].ASME Journal of Applied Mechanics,1965,32(4):911-920.
[11]HAN X,LIU G R.Computational inverse technique for material characterization of functionally graded materials[J].American Institute of Aeronautics and Astronautics,2003,41(2):288-295.
[12]LIU G R,HAN X.Computational inverse techniques in nondestructive evaluation[M].Raton:Boca,2003.
[2]DING Y,ZHAO B Y,WU B,et al.A condition assessment method for time-variant structures with incomplete measurements[J].Mechanical Systems and Signal Processing,2015,58:228-244.
[3]FRANULOVIC M,BASAN R,PREBIL I.Genetic algorithm in material model parameters’identification for low-cycle fatigue[J].Computational Materials Science,2009,45(2):505-510.
[4]AGUIR H,BELHADJSALAH H,HAMBLI R.Parameter identification of an elasto-plastic behavior using artificial neural networks-genetic algorithm method[J].Materials and Design,2011,32(1):48-53.
[5]LIEW A,FENG N S,HAHN E J.A non-linear transfer matrix technique for statically indeterminate rotor bearing systems[J].Proceedings of the Institution of Mechanical Engineers,Part C:Journal of Mechanical Engineering Science,2001,215(11):1343-1355.
[6]STANDER N,ROUX W,GEOL T,et al.Ls-opt user's manual[M].California:Livermore Software Technology Corporation,2009.
[7]伍乾坤,韩旭,胡德安.一种陶瓷脆性材料动态本构参数的计算反求方法[J].固体力学学报,2009,30(3):280-285.(WU Qiankun,HAN Xu,HU Dean.An inverse technique for identification of dynamic constitutive parameters of a ceramic brittle material[J].Chinese Journal of Solid Mechanics,2009,30(3):280-285(in Chinese)).
[8]Thakur M.A new genetic algorithm for global optimization of multimodal continuous functions[J].Journal of Computational Science,2014,5(2):298-311.
[9]闻邦椿,顾家柳,夏松波,等.高等转子动力学[M].北京:机械工业出版社,2000.(WEN Bangchun,GU Guliu,XIA Songbo,et al.Higher rotor dynamics[M].Beijing:Mechanical Industry Press,2000(in Chinese)).
[10]LUND J W.The stability of an elastic rotor in journal bearing with flexible damped supports[J].ASME Journal of Applied Mechanics,1965,32(4):911-920.
[11]HAN X,LIU G R.Computational inverse technique for material characterization of functionally graded materials[J].American Institute of Aeronautics and Astronautics,2003,41(2):288-295.
[12]LIU G R,HAN X.Computational inverse techniques in nondestructive evaluation[M].Raton:Boca,2003.