分数阶微积分及其在粘弹性材料和控制理论中的应用

论文摘要

“数学是给予万物高深莫测的别名的艺术。而‘分数阶微积分’这一美丽的-乍听起来有些神秘的一名称正是这些别名中的一个。它同时也是数学理论的精华.”——Igor Podlubny,1999。到目前为止,分数阶微分算子和分数阶积分算子在粘弹性理论中得到了最广泛的应用。许多文献中提到应用分数阶微分作为粘弹性材料的数学模型是件很自然的事情。值得一提的是,分数阶理论之所以在粘弹性材料建模上得到如此大的发展的主要原因是分数阶材料在工程领域的广泛应用。并且,只要给予适当的假设,几乎所有的分数阶模型都能很好的归纳材料的变形特征,并对解释实际问题起着相当大的作用。基于这些理论和实践的研究,并且考虑到粘弹性材料的四大特征,松弛、蠕变、预调和滞后,徐明瑜教授和我提供了一种新的粘弹性材料特性的数学描述并将其作为本学术论文的第一大部分内容。这些理论可以直接指导各种材料试验并且丰富了现有的结果。这篇博士论文的另一部分阐述的是分数阶微积分在控制理论中的应用。目前,受许多物理现象的启发,世界范围内的学者们应用含有分数阶导数和积分的分数阶本构关系和分数阶系统方程在不同领域做出了越来越多的贡献。这些新模型比起经典的整数阶模型更胜一筹。另外,分数阶导数和积分为描述不同事物的记忆性和遗传性提供了一个有力的工具。这是分数阶模型最明显的特征之一。同时,也是分数阶粘弹性模型在实际中得到如此成功实践的原因之一。尽管如此,由于缺少相应的数学工具,分数阶动力系统的研究只占据了控制系统理论和应用很小的一部分空间。尽管在分数阶控制领域也有一些成功的尝试,但总的来说在时域中的讨论几乎被避开了。为了填补这项空白,我们（我的导师,合作者和我自己）发表了一系列这方面的学术论文并希望这些有关基础理论的文章有助于分数阶微积分在理论和应用方面的发展。首先,我们导出了分数阶Lyapunov直接方法,它是系统渐进稳定的充分条件。同时,我们还提出了Mittag-Leffler稳定的概念,它是指数稳定的广义形式。另外,我们还研究了分数阶通用自适应稳定理论。基于这些理论,我们还做了一系列的有很好结果的试验,但基于本文的连贯性和可读性,这一部分被省略掉了。第三,在对分数阶通用自适应稳定问题的研究中,我们回答了“什么时候（广义）Mittag-Leffler函数是Nussbaum函数?”。这一理论提供了一类全新的Nussbaum函数。通过应用这一类新的函数,我们得到了许多成功的试验和仿真结果。最后,在对上述问题的研究中,我们得到了许多灵感和启示,丰富了我们下一步的研究内容。总之,在未来的时间里将会涌现出更多的“分数阶”物理问题。我们引用下面的话来作为本文陈述部分的结束:“…所有的系统需要一个含有分数阶导数的方程来描述…系统对所有之前发生的事情都有记忆性。用发生过的事情来预测未来是有必要的。…这个结论是明显且不可避免的:‘逝去的事物仍有记忆’。换句话说,一切自然事物都是分数阶的。”——S.Westerlund,1991。本学位论文的正文部分共分六章:在第一章中,我们介绍了Riemann-Liouville和Caputo分数阶算子的基本概念和性质。另外,我们还列出了Mittag-Leffler函数,两参数Mittag-Leffler,广义Mittag-Leffler函数（三参数Mittag-Leffler函数）和四参数Mittag-Leffler函数的基本定义。这些函数是许多分数阶微分方程的基本解。此外,我们给出了Fox H-函数的定义和无穷项级数表示。这种函数包含了几乎所有在应用数学和统计学领域被广泛应用的特殊函数:Mittag-Leffler类型函数,广义超越几何函数,广义Bessel函数和Meijer’s G-函数,等等。在本章的最后,我们对分数阶微积分及其在粘弹性材料和控制理论中应用的发展作了精要概括。本章内容是一个分数阶微积分及其应用的基本介绍和精要概括。所有的这些定义和内容都在随后的各个章节中得到了应用,并且,这些定义和特殊函数的性质使本文更易于读者理解。在第二章中,我们首先证明了整数阶Kelvin模型在锯齿波加载下的滞后和预调。随后,我们证明了产生滞后和预调的模型和输入的充分条件。我们还讨论了滞后和预调的单调和极限性质。我们证明了滞后回环的宽度单调递减趋于非负值。同时我们还证明了,应用分数阶Kelvin模型和准线性理论以上结果仍然成立。我们是第一个给出分数阶粘弹性固体材料滞后和预调的数学描述的科研小组。这一工作给许多材料学试验提供了理论依据和支持。另外,基于各式各样真实的材料试验,我们可以通过调整模型参数来更好的描述和指导分数阶模型在更大范围内的应用和实践。在第三章中,我们讨论了粘弹性固体材料模型滞后回环和能量损耗的变化和趋势。我们的结论之一是:（3.4.15）式的正负是判断第（n+1）周期和第n个周期能量损耗变化的充要条件。我们还证明了Boltzmann叠加原理对于周期输入的适用性。应用分数阶Kelvin模型和准线性理论,本章中上述结论仍成立。在真实材料的运动中,能量损耗几乎是不可避免的现象。尤其是在周期加载下能量损耗的趋势和极限性质是一项非常有意义的工作,这也是本章的主要内容。基于整数阶模型和分数阶模型,叠加原理和拟线性原理,这一问题得到了较好的回答并希望其结论有助于指导实际材料试验。在第四章中,我们提出了Mittag-Leffler稳定的定义和分数阶里亚普诺夫直接方法。我们还证明了分数阶比较原理。另外,应用Caputo分数阶微积分算子,我们扩展了Riemann-Liouville分数阶系统的应用。为了支持这些理论,我们还提供了一个实例。Mittag-Leffler稳定这一新定义是经典的指数衰减的扩展。它还可以在更广泛的意义下来描述整数阶和分数阶系统的稳定性理论。在第五章中,应用分数阶通用自适应稳定理论,我们研究了三种分数阶标量系统的稳定性问题。这三种系统是:（Ⅰ）分数阶动力,整数阶控制;（Ⅱ）分数阶动力,分数阶控制;（Ⅲ）含有分数阶粘弹性材料和分数阶电磁力学模型的通用自适应稳定问题。文章证明了这三种情况下控制量的有界性。另外,我们还首次证明了对于Mittag-Leffler函数Eα（—λtα）当α∈(2,3]和λ＞0时是Nussbaum函数。几个仿真结果验证了以上的理论。通用自适应稳定理论看似一件神奇的事情,因为,它可以在不知道系统参数的情况下保证系统的稳定性.这一理论在整数阶系统中已经得到了完善的发展。但是,它在分数阶系统中的应用还是一个新课题。另外,我们已经证明了应用当Mittag-Leffler函数是Nussbaum函数时这一理论在理论和实验中都是非常有意义的工作。在第六章中,我们建立了广义Mittag-Leffler函数和Nussbaum函数之间的关系。另外,当广义Mittag-Leffler函数是Nussbaum函数时,这些广义Mittag-Leffler函数的导数仍然是Nussbaum函数。其中的Matlab图像应用了YangQuan Chen老师编辑的M-文件。最后,我们讨论了Hessenberg矩阵的行列式问题。这一工作意义非常,因为,当我们把“（广义）Mittag-Leffler函数是Nussbaum函数”应用到直流永磁电动机的通用自适应稳定实验中时,得到了非常好的结论。出于连贯性的考虑,这一部分并没有出现在本文中。

论文目录

• 中文部分
• 符号说明
• 中文摘要
• 英文摘要
• 第一章 基础知识
• §1.1 分数阶微积分概述
• §1.1.1 Riemann-Liouville分数阶导数
• §1.1.2 Caputo分数阶导数
• §1.2 Mittag-Leflter类型函数的概述
• §1.3 Fox H-函数概述
• §1.4 分数阶微积分在粘弹性材料和控制理论中的应用的概述
• §1.4.1 分数阶粘弹性材料及模型
• §1.4.2 控制理论中的分数阶微积分
• 第二章 粘弹性固体模型的滞后和预调
• §2.1 简介
• §2.2 Kelvin模型和锯齿波输入
• §2.3 滞后和预调
• §2.3.1 滞后
• §2.3.2 预调
• §2.4 关于滞后和预调的进一步讨论
• §2.5 分数阶Kelvin模型和锯齿波输入
• §2.6 更一般的应变输入
• §2.7 拟线性理论下滞后和预调的讨论
• §2.8 小结
• 第三章 粘弹性固体模型的滞后回环和能量损耗
• §3.1 简介
• §3.2 粘弹性固体模型基本理论
• §3.3 非负应变输入下的滞后回环和能量损耗
• §3.3.1 情况一
• §3.3.2 情况二
• §3.3.3 情况三
• §3.3.4 情况四
• §3.4 应变值正负相间的情况
• §3.5 分数阶Kelvin模型和拟线性理论的应用
• §3.6 结果和讨论
• 第四 章分数阶非线性系统的Mittag-Leffler稳定性
• §4.1 简介
• §4.2 分数阶非自治系统
• §4.2.1 Lipschitz条件和Caputo分数阶非自治系统
• §4.3 Mittag-Leffler稳定性
• §4.4 Lyapunov直接方法的分数阶拓展
• §4.5 引用κ-类函数的分数阶Lyapunov直接方法
• §4.6 例题
• §4.7 结论和展望
• 第五章 分数阶通用自适应稳定
• §5.1 简介
• §5.2 基础知识之Nussbaum函数
• §5.3 分数阶UAS的稳定性分析
• §5.3.1 分数阶系统整数阶控制
• §5.3.2 分数阶系统分数阶控制
• §5.3.3 粘弹性理论和电磁学理论在UAS中的应用
• §5.4 Mittag-Leffler函数为Nussbaum函数的情况
• §5.5 仿真
• §5.6 结论和展望
• 第六章 什么时候广义Mittag-Leffler函数是Nussbaum函数
• §6.1 简介
• §6.2 基础知识
• §6.2.1 Nussbaum函数及其性质
• §6.2.2 一个特殊Hessenberg矩阵的行列式
• §6.3 广义Mittag-Leffler函数的分解
• ρ,μ,γ（t）'>§6.3.1 有界项f_ρ,μ,γ（t）
• ρ,μ,γ（t）'>§6.3.2 震荡项g_ρ,μ,γ（t）
• §6.4 图形
• §6.5 小结
• 本文的主要贡献和未来的展望
• 参考文献
• 致谢
• 攻读博士学位期间完成论文情况
• 个人简历
• 学位论文评阅及答辩情况表
• 英文部分
• Notation
• Abstract
• 摘要
• 1 Elementary Knowledge
• 1.1 Synopsis of Fractional Calculus
• 1.1.1 The Riemann-Liouville Fractional Derivative
• 1.1.2 The Caputo Fractional Derivative
• 1.2 Synopsis of Mittag-Leffler Type Functions
• 1.3 Synopsis of Fox H-Function
• 1.4 A Survey of the Applications of Fractional Calculus to Viscoelastic Material and Control Theory
• 1.4.1 Fractional Viscoelastic Materials and Models
• 1.4.2 Fractional Calculus in Control Theory
• 2 Hysteresis and Precondition of Viscoelastic Solid Models
• 2.1 Introduction
• 2.2 Kelvin Model and Saw-Tooth Wave Input
• 2.3 Hysteresis and Precondition
• 2.3.1 Hysteresis
• 2.3.2 Precondition
• 2.4 Further Discussion About Hysteresis and Precondition
• 2.5 Fractional-Order Kelvin Model and Saw-Tooth Wave Input
• 2.6 More General Situations of Inputted Strains
• 2.7 Discussion of Hysteresis and Precondition with Quasi-Linear Theory
• 2.8 Concluding Remarks
• 3 Hysteresis Loop and Energy Dissipation of Viscoelastic Solid Models
• 3.1 Introduction
• 3.2 The Basic Theories of Viscoelastic Solid Models
• 3.3 Hysteresis Loop and Energy Dissipation Under Normegative Inputted Strain
• 3.3.1 Situation 1
• 3.3.2 Situation 2
• 3.3.3 Situation 3
• 3.3.4 Situation 4
• 3.4 Situation of Strain's Value Being Positive or Negative,Alternative
• 3.5 Applications of Fractional Kelvin Model and Quasi-Linear Theory
• 3.6 Results and Discussion
• 4 Mittag-Leffler Stability of Fractional Order Nonlinear Dynamic Systems
• 4.1 Introduction
• 4.2 Fractional Nonautonomous Systems
• 4.2.1 Lipschitz Condition and Caputo Fractional Nonautonomous Systems
• 4.3 Mittag-Leffler Stability
• 4.4 Fractional Order Extension of Lyapunov Direct Method
• 4.5 Fractional Lyapunov Direct Method by Using The Class-κ Functions
• 4.6 An Illustrative Example
• 4.7 Conclusion and Future Works
• 5 Fractional Order Universal Adaptive Stabilization
• 5.1 Introduction
• 5.2 Preliminary:Introduction of the Nussbaum Function
• 5.3 Stability Analysis of Fractional UAS
• 5.3.1 Fractional Dynamics with Integer-order Control Strategy
• 5.3.2 Fractional Dynamics with Fractional Control
• 5.3.3 Application of Viscoelastic and Electromagnetic Theories in UAS
• 5.4 Mittag-Leffler Function as Nussbaum Function
• 5.5 Simulations
• 5.6 Conclusion and Future Works
• 6 When Generalized Mittag-Leffler function is Nussbaum function
• 6.1 Introduction
• 6.2 Preliminaries
• 6.2.1 Nussbaum Function and its Properties
• 6.2.2 Determinant of A Special Hessenberg Matirx
• 6.3 Decomposition of Generalized Mittag-Leffler Functions
• ρ,μ,γ（t）'>6.3.1 The Bounded Term f_ρ,μ,γ（t）
• ρ,μ,γ（t）'>6.3.2 The oscillated term g_ρ,μ,γ（t）
• 6.4 Figures
• 6.5 Conclusions
• Contributions and Recommendations
• Bibliography
• Acknowledgement
• List of Publications during Study for the Doctorate
• Curriculum Vitae
• 学位论文评阅及答辩情况表