分数阶导数振动模型中阻尼力与速度关系的分析.pdfVIP

分数阶导数振动模型中阻尼力 与速度关系的分析 胡帅 张晓棣 陈文* （河海大学工程力学系，软物质力学研究所，南京 210098） 摘要：分数阶导数广泛用于描述复杂介质中振动系统的频率依赖阻尼问题。但至今我们仍未 见到有关这类模型中阻尼力与速度关系的详细研究。本文利用有限差分法模拟了分数阶导数 振动模型中阻尼力与速度的关系，通过数值实验结果揭示了一些新的现象。 关键词：分数阶导数；频率依赖阻尼；振动；阻尼力；速度 振动在工业、交通、航空等很多领域都广泛存在。由于振动所带来了各种各样的不良影 响如材料老化、结构破坏等，都可能造成不同程度的损失。因此，有关减振的研究和利用在 工业以及建筑行业具有重要的意义。影响振动状态的阻尼力，尤其是阻尼力与振动速度的关 系是各种抗振和振动控制技术中的关键问题。至今许多学者在该领域已做了大量的研究工 作，提出了多个结构动力系统振动阻尼的模型，其中主要有负常数模量模型（也被称为结构 [1] 阻尼）、标准流变学模型、分数阶导数模型、分数指数模型及微振子模型 。 研究人员和工程师在复杂材料的研究中经常遇到所谓的频率依赖振动阻尼[2-5] 。例如高 速铁路路基的振动，磁流变抗震阻尼器等。这类阻尼的实验测量或现场观测结果通常拟合为 一个以频率为变量的幂函数c ωα [6,7] 。这里α∈[0,2]可为任意实数，ω代表角频率，参数 c0 0 和α则由测量数据拟合得到，是依赖于介质和结构的常数。整数阶导数的经典理论适用于描 述理想固体和液体材料，而频率依赖阻尼材料经常为介于理想固体和液体之间的多相软物质 介质，其性质与理想固体和液体有很大的差别。大量理论和实验分析表明，经典理论模型不 能够准确地描述复杂材料的频率依赖阻尼行为。因此，必须寻找其他的模型来替代。 目前，国内外广泛采用分数阶导数模型模拟频率依赖阻尼。理论和实验研究证实这类模 型的误差相对较小，物理意义清晰。分数阶导数模型最早在上世纪三十年代由 Gement 提出， 应用于描述粘弹性材料的本构关系。我们可以用s = E e 来表示弹性固体材料的本构关系， 牛顿流体的本构关系则表示为s = cde t / dt 。由此，Gement 认为介于理想弹性固体和牛 ( ) 顿流体之间的粘弹性材料本构关系可用分数阶导数表达为s = qd a e (t ) / dta 。至今已经出现 了多种分数阶导数模型。如 1/2 阶非线性分数阶微分方程模型和 Bagley-Torvik 方程模型等。 虽然对这类方程已有较多的研究，但现有工作的注意力是建立一个更好地描述粘弹性阻尼振 动现象的分数阶导数模型，而至今未对其中阻尼力项进行深入的探讨和分析。 阻尼力是实际工程振动控制中要考虑的关键因素，决定了速度衰减以及振幅的衰减都直 接。因此，阻尼力的变化规律及阻尼力与可测量速度的关系对减振和抗振分析与设计有重要 的指导作用。本文利用有限差分法模拟了分数阶导数振动模型中阻尼力与速度的关系，并尝 试对此类模型中阻尼力与速度关系进行理论解释。 本文结构如下：第一节主要介绍分数阶阻尼模型；第二节分析匀加速运动分数阶阻尼力 （广义速度）与速度的关系；第三节通过单自由度分数阶阻尼振动算例，分析了阻尼力与速 度的关系；第四节是结论和展望。 *基金项目：2007 年河海大学大学生创新训练计划项目（N2007027 ），国家自然科学基金面上项目 作者简介：胡帅（1987- ）本科生 hswyy@hhu.edu.cn 张晓棣（1985- ）硕士研究生zhangxd1985@hhu.edu.cn 陈文（1967- ）通讯作者，博士生导师chenwen@hhu.edu.cn 第 758 页 1. 分数阶阻尼模型简介 对于弹性固体由广义胡克定律可得其本构方程为：