非线性分析：几何非线性_（6）.非线性动力学.docxVIP

PAGE1

PAGE1

非线性动力学

1.引言

非线性动力学是研究非线性系统在时间演化中的行为的学科。在工程力学和固体力学领域，非线性动力学的应用非常广泛，尤其是在处理大变形、材料非线性、接触问题和多体系统动力学等方面。本节将详细介绍非线性动力学的基本概念、数学模型和数值方法，并通过具体的应用实例来说明这些概念和技术的实现。

2.非线性动力学的基本概念

2.1非线性系统

非线性系统是指其行为不能用线性关系来描述的系统。与线性系统相比，非线性系统的行为更加复杂，可能表现出多稳态、混沌、分岔等现象。在固体力学中，常见的非线性系统包括大变形的结构、非线性材料（如橡胶、塑料等）和接触问题。

2.1.1非线性方程

非线性方程是指方程中包含非线性项的方程。例如，一个简单的非线性方程可以是：

u

其中，u是位移，α和β是线性阻尼和刚度系数，γ是非线性刚度系数，ft

2.2相空间和相轨迹

相空间是描述系统状态的多维空间，每个维度表示一个状态变量。相轨迹是系统在相空间中的时间演化路径。通过相空间和相轨迹，可以直观地分析非线性系统的动态行为。

2.2.1相空间的构建

假设我们有一个两自由度的非线性系统，其状态变量为u1和u2，以及它们的导数u1和

x

系统的动力学方程可以写成：

x

其中，f是一个非线性函数。

2.3稳定性和分岔

稳定性是指系统在受到小扰动后能否回到平衡状态。分岔是指系统在参数变化时出现新的行为模式。非线性系统中常见的分岔类型包括鞍节点分岔、叉形分岔和霍普夫分岔。

2.3.1鞍节点分岔

鞍节点分岔发生在系统的一个平衡点消失，同时出现两个新的平衡点。这种分岔通常导致系统的多稳态行为。

2.3.2叉形分岔

叉形分岔发生在系统的一个平衡点分裂成两个或多个新的平衡点。这种分岔导致系统的对称性破缺。

2.3.3霍普夫分岔

霍普夫分岔发生在系统的一个平衡点失去稳定性，同时出现一个周期解。这种分岔导致系统的振荡行为。

3.数学模型

3.1非线性振动方程

非线性振动方程是描述非线性系统振动行为的基本方程。常见的非线性振动方程包括Duffing方程和VanderPol方程。

3.1.1Duffing方程

Duffing方程是一个典型的非线性振动方程，形式为：

u

其中，δ是阻尼系数，α和β是刚度系数，γ是外力幅值，ω是外力频率。

3.1.2VanderPol方程

VanderPol方程是一个描述自激振荡的非线性振动方程，形式为：

u

其中，μ是非线性阻尼系数。

3.2非线性微分方程

非线性微分方程是描述非线性系统动态行为的数学工具。常见的非线性微分方程包括Lorenz方程和R?ssler方程。

3.2.1Lorenz方程

Lorenz方程是一个描述大气对流的非线性微分方程组，形式为：

x

其中，σ、ρ和β是系统参数。

3.2.2R?ssler方程

R?ssler方程是一个描述化学反应的非线性微分方程组，形式为：

x

其中，α、β和γ是系统参数。

4.数值方法

4.1显式和隐式时间积分方法

显式和隐式时间积分方法是求解非线性动力学问题的常用数值方法。显式方法计算速度快，但稳定性较差；隐式方法计算速度较慢，但稳定性较好。

4.1.1显式时间积分方法

显式时间积分方法中最常用的是Newmark-β方法。该方法的公式为：

u

其中，Δt是时间步长，β

4.1.2隐式时间积分方法

隐式时间积分方法中最常用的是中心差分法。该方法的公式为：

u

4.2Newmark-β方法

Newmark-β方法是一种广泛应用于结构动力学问题的隐式时间积分方法。该方法通过选择不同的参数α和β来平衡计算速度和稳定性。

4.2.1Newmark-β方法的原理

Newmark-β方法的基本原理是将位移、速度和加速度在时间步长Δt

u

u

其中，α和β是Newmark参数，通常取α=0.5和β

4.3中心差分法

中心差分法是一种显式时间积分方法，适用于线性问题和某些非线性问题。该方法通过在时间步长Δt

4.3.1中心差分法的原理

中心差分法的基本原理是：

u

u

4.4有限元方法

有限元方法（FEM）是一种将连续体离散化为有限个单元的方法，广泛应用于非线性动力学问题的数值求解。

4.4.1有限元方法的基本原理

有限元方法的基本原理是将结构划分为若干个单元，每个单元用节点位移来描述。通过节点位移的插值函数，可以将单元内的位移场表示为：

u

其中，Nix是插值函数，u

4.4.2有限元方法的应用

有限元方法可以应用于大变形、材料非线性和接触问题的数值模拟。例如，对于一个大变形的弹性梁，可以使用有限元方法来求解其非线性动力学方