超导磁悬浮动力运动阻尼特性分析.doc

超导磁悬浮动力运动阻尼特性分析( 苟晓凡(( 河海大学 工程力学系 南京 210098 摘要：在超导磁悬浮系统中，悬浮体在外界干扰下发生振动，其中阻尼特性是影响系统动力运动很重要的一个方面。为此，本文对超导磁悬浮动力运动过程的阻尼特性展开理论研究。针对超导磁悬浮动力运动的强非线性及多物理场耦特征，建立了描述悬浮体动力运动的数学模型。通过数值求解，定量研究了系统阻尼随时间的变化及其与系统主要物理参数的关系，并给出了系统阻尼对系统动力运动的定量影响结果。 引言 超导悬浮作为超导体的重要应用之一，与其它传统的物理悬浮相比较，具有十分突出的优点，广泛应用于超导磁悬浮列车及悬浮轴承等领域。在系统悬浮体动态响应过程中，其位移随时间均有不同程度的衰减[1]，由此可见系统中有阻尼存在。除过空气阻尼外，对系统阻尼的研究对全面把握系统动力运动特性是至关重要的。许多研究者仅认识到系统阻尼比较小，且在增加系统阻尼方面作了一些工作[2,3,4]，但系统动态响应过程中影响阻尼的因素以及作用规律也不完全清楚。目前，针对超导悬浮系统动态响应过程中阻尼特性的定量研究尚未见到。为此，本文将对系统阻尼展开定量研究。在充分考虑到系统强非线性和多物理场耦合特征基础上，建立系统悬浮体动力运动的理论分析框架，并对系统阻尼的作用规律及其影响因素展开研究，得到相应的定量结果。 基本方程 基本方程 最为典型的超导悬浮系统如图1所示。圆柱形的超导体（HTS）和永磁体（PM）轴线重合，超导体处于超导态，永磁体静止悬浮于静平衡位置。在外界初始干扰下，悬浮体（永磁体及附加质量）发生振动，描述其动力运动的微分方程为 , (1) 初始条件可写为 : . (2) 其中 ，为系统的静平衡位置（）。是超导体对永磁体的电磁悬浮力，它由Lorentz 公式 (3) 给出。这里是超导体内某一确定空间位置的电流密度；代表磁感应强度，其中和分别是相应位置的外加磁场以及诱导电流所引起磁场的磁感应强度。而超导体的电流密度和磁感应强度必须满足maxwell方程组 , (4) 电荷守恒定律 。其中为梯度算子，为超导体中的磁场矢量，为真空中的磁导率。 更为重要的是，超导体作为特殊的电磁材料，其电磁物理量必须满足超到电磁本构关系。对于超导体而言，当其处于交变磁场或悬浮系统动力运动过程中，磁通流动与蠕动模型是合适的，即 (5) 式中，为磁通蠕动阻抗， 临界电流密度，钉扎势，伯茨曼常数，绝对温度，磁通流动阻抗。 需要指出的是，式（5）所展示的超导本构关系是非线性的。此外，由(1)-(3)看出，欲求解系统动力运动方程，必先已知电磁悬浮力，其前提是已知超导体中的电磁场分布，即要联立求解Maxwell方程组和超导本构关系(5)。而要联立求解Maxwell方程组和超导本构关系，由必须已知超导体和永磁体的相对位置和速度。显然，求解系统动力运动方程和求解超导内部电磁场分布是相互耦合在一起的。因此，求解超导悬浮系统动力运动实际上是解决多物理场耦合和强非线性关系的问题。 阻尼系数 评价系统阻尼对振动系统的影响，一个较为直观的方法就是考察其动态响应曲线振幅衰减情况。为此，引入减幅系数的概念，它定义为相邻两个振幅的比值： (6) 其中和分别为第和第个周期的振幅。实际中为了方便常用下面的对数衰减率取代： (7) 对于小阻尼情形（相对阻尼系数），省略中间分析推导过程，直接给出相对阻尼系数为 (8) 这是由振幅衰减曲线的对数衰减率求相对阻尼系数的一种较为有效的方法[5,6]。 结果与讨论 为了清楚的展示系统阻尼的作用，这里暂不考虑空气阻尼的影响，即具体计算时令空气阻尼系数，并且以自由振动为例来研究。首先，定量模拟了系统在不同的超导材料参数、初始冷却高度以及悬浮高度处的动态响应曲线，并发现相应的曲线衰减程度不同，由此可见，在以上不同的条件下，系统动态响应过程中的阻尼是不同的。为此，我们将定量研究以上各个量对超导悬浮系统动态响应过程中阻尼的影响程度及其作用规律。图2给出了不同材料参数下相对阻尼系数随时间的变化。从图2中的四个图可明显看出，相对阻尼系数随时间均不断减小，在初始阶段，随时间下降较快，到稳定阶段后，其变化逐渐减缓。以下对每个图给予具体分析。由图2 (a)可看出，临界电流密度越大，相对阻尼系数反而越小。主要因为，临界电流密 (a) 临界电流密度 (b) (c)　钉扎势