非线性分析：粘弹性材料分析_（7）.非线性动力学分析.docxVIP

PAGE1

PAGE1

非线性动力学分析

在前一节中，我们已经探讨了粘弹性材料的基本特性和分析方法。接下来，我们将深入讨论非线性动力学分析，这是工程力学中一个重要的分支，特别是在涉及粘弹性材料时。非线性动力学分析主要用于研究系统在外力作用下的动力响应，特别是在非线性行为显著的情况下。这一节将详细介绍非线性动力学的基本原理、数值方法以及实际应用中的案例分析。

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

非线性动力学是指研究系统在受到外力作用时，其动力响应表现出非线性特征的学科。非线性特征通常包括几何非线性、材料非线性和边界条件非线性。对于粘弹性材料，材料非线性是主要的研究对象，因为其应力-应变关系通常不是线性的，且具有时间依赖性。

1.1几何非线性

几何非线性主要发生在大变形和大位移的情况下，当结构的变形足够大时，初始几何形状和变形后的几何形状之间的差异显著，导致线性假设不再适用。在这种情况下，需要考虑结构变形对自身几何形状的影响。

1.2材料非线性

材料非线性是指材料的应力-应变关系不是线性的。对于粘弹性材料，其应力-应变关系不仅依赖于应变，还依赖于应变率和时间历史。常见的粘弹性模型包括Maxwell模型、Kelvin-Voigt模型和标准线性固体模型。

1.3边界条件非线性

边界条件非线性是指系统边界条件随时间或状态的变化而变化。例如，接触问题中的接触力和接触面的变化，或者温度场的变化导致的热应力变化。

2.非线性动力学分析的数学基础

非线性动力学分析的数学基础主要包括动力学方程、非线性方程的求解方法以及数值积分技术。

2.1动力学方程

非线性动力学方程通常可以表示为：

M

其中：-M是质量矩阵-C是阻尼矩阵-Ku是非线性刚度矩阵，依赖于位移u-Ft

2.2非线性方程的求解方法

非线性方程的求解方法主要有牛顿-拉夫森法（Newton-RaphsonMethod）和弧长法（Arc-LengthMethod）。

2.2.1牛顿-拉夫森法

牛顿-拉夫森法是一种迭代方法，用于求解非线性方程组。其基本步骤如下：1.选择一个初始猜测值u0。2.计算残差向量Ru=Mu+Cu+Ku?F

u

5.检查收敛条件，如果满足则停止迭代，否则返回步骤2。

2.2.2弧长法

弧长法是一种路径跟踪方法，用于解决非线性方程在某些情况下难以收敛的问题。其基本思想是通过引入一个虚拟参数λ来控制迭代路径，使得每一步的位移和外力增量之和保持在一定的弧长范围内。

2.3数值积分技术

数值积分技术用于求解动力学方程中的时间导数项。常见的数值积分方法包括中心差分法（CentralDifferenceMethod）、Newmark法和Wilson-θ法。

2.3.1中心差分法

中心差分法是一种显式时间积分方法，适用于求解高频率问题。其基本公式如下：

u

u

2.3.2Newmark法

Newmark法是一种隐式时间积分方法，适用于求解低频率问题。其基本公式如下：

u

u

2.3.3Wilson-θ法

Wilson-θ法是一种改进的隐式时间积分方法，适用于求解强非线性问题。其基本公式如下：

u

u

3.非线性动力学分析的数值方法

非线性动力学分析的数值方法主要包括有限元法（FiniteElementMethod,FEM）、有限差分法（FiniteDifferenceMethod,FDM）和边界元法（BoundaryElementMethod,BEM）。

3.1有限元法

有限元法是一种广泛应用于非线性动力学分析的方法。其基本步骤如下：1.将连续体离散化为有限个单元。2.建立单元的运动方程。3.组合单元方程，形成整体结构的运动方程。4.求解整体结构的运动方程。

3.1.1单元运动方程

对于一个粘弹性材料单元，其运动方程可以表示为：

M

其中：-Me是单元质量矩阵-Ce是单元阻尼矩阵-Keue是单元非线性刚度矩阵

3.2有限差分法

有限差分法是一种基于差分逼近的数值方法，适用于求解偏微分方程。其基本步骤如下：1.将连续体离散化为网格点。2.在每个网格点上建立差分方程。3.求解差分方程组。

3.2.1差分方程

对于一维粘弹性材料的非线性动力学方程，可以使用中心差分法建立差分方程：

M

3.3边界元法

边界元法是一种基于边界积分方程的数值方法，适用于求解边界条件复杂的非线性动力学问题。其基本步骤如下：1.将连续体的边界离散化为有限个边界元。2.在每个边界元上建立边界积分方程。3.求解边界积分方程组。

3.3.1边界积分方程

对于一个粘弹性材料的非线性动力学问题，边界积分方程可以表示为：

Γ

其中：-u*是基本解