非线性有限元 第1章绪论.pdfVIP

第一章 绪论 1.1 引言 有限单元法已成为一种强有力的数值解法来解决工程中遇到的大量问题，其应用范围从 固体到流体，从静力到动力，从力学问题到非力学问题。事实上，有限单元法已经成为在已 知边界条件和初始条件下求解偏微分方程组的一般数值方法。 对于有限元的线性分析，我们已经比较熟悉，它已作为设计工具被广泛采用。但绝大多 数实际问题属于非线性问题，根据产生非线性的原因，非线性问题主要有三种类型： 1．材料非线性(简称材料非线性或物理非线性) 其特点是应力 与应变之间为非线性关系，通常与加载历史有关，加载和卸载不是 同一途径（图 1.1），因而其物理方程 σ Dε 中的弹性矩阵 D 是应变ε 的函数。但材料非线性问题仍属于小变形问题，位移和应变是微 量，其几何方程是线性的。土、岩石、混凝土等具有典型的材料非线性性质，所以，土坝、 岩土地基的稳定性和加固，地下洞室和边坡的稳定性都应当按材料非线性问题处理。 2 ．几何非线性 几何非线性属于大变形问题，位移和应变或者它们中一 个是有限量。可能会有三种情况:大位移(包括线位移和角位 移) 、小应变，小位移、大应变和大位移、大应变。此时反 映应变和位移关系的几何方程是非线性方程，例如，正应变  可表示为 z u u v w  1  2  2  2 x  [( ) ( ) ( ) ]  x x x x  2    剪应变 表示为 xy v u u u v v w w    xy      x y x y x y x y 如果应力和应变之间的关系也是非线性的，就变成了更复杂的双重非线性问题。不过，在几 何非线性问题中一般都认为应力在弹性范围内， 与之间呈线性关系。工程中的实体结 构和板壳结构都存在几何非线性问题，例如弹性薄壳的大挠度分析，压杆或板壳在弹性屈曲 后的稳定性问题。 在采用有限元方法分析非线性问题时，以上两类都表现为结构的整体劲度矩阵 K 不再 是常量矩阵，而是结点位移δ 的函数，还有一类问题是结点荷载 R 与δ有关，这就是边界 非线性问题，又称接触非线性。 3 ．接触非线性 由于接触体的变形和接触边界的摩擦作用，使得部分边界条件随加载过程而变化，且不 可恢复。这种由边界条件的可变性和不可逆性产生的非线性问题，称为接触非线性。工程中 有许多接触非线性的问题，如混凝土坝纵缝和横缝缝面的接触，面板堆石坝中钢筋混凝土面 板与垫层之间的接触，岩体节理面或裂隙面的工作状态等。 目前，研究工程非线性问题比较有效的工具是非线性有限元方法。要使这一方法实用化， 有两个问题必须解决。第一，由于非线性问题的数值计算工作量大大增加，需要有相当大计 1 算能力的工具。近 10 年来，高速计算机的发展已基本上满足了这一需要，同时计算费用也 在继续减小。第二，非线性求解方法的精度和收敛性必须被验证。由于发展和改进了多种类 型的单元，以更好地模拟结构的工作，找到更有效的非线性求解方法，并且积累了许多经验 可应用于实际工程问题，现在已经能够比较有把握地完成非线性有限元分析。非线性有限元 方法正在成为一种强有力的计算工具被研究人员和工程人员所采用。 1.2 线弹性有限元方法的回顾 一般来说，非线性有限元法可归结为一系列线弹性问题。因而，线弹性有限元是非线性 有限元法的基础。两者不仅在分析方法和求解步骤上有相似之处，而且后者要不断调用前者 的结果。为此，本节将扼要地