变分法简介精要.pptVIP

* 结构分析是有限单元法最早、也是最广泛应用的领 域。 前面以弹性力学平面问题为例，阐释了有限单元法的 基本内容。这样的介绍具有直观性，但缺乏系统性和深刻 性。为加深对有限单元法的理解，本章将系统而深入地阐述 有限单元法的基本原理， 内容包括： (1)介绍定解问题的微分方程提法； (2)根据微分方程的等效积分形式，推导虚位移原理及势能变分原理，从而建立定解问题的泛函变分提法； (3)基于势能变分原理推导位移有限单元法的普遍公式，并对 位移有限元解的性质和收敛性作简要说明； 4.1微分方程提法 在物理或工程问题中，位移、应力、温度、电流等物 理量称为场变量，它们在一定区域内满足某些控制方程； 在域边界上满足给定的边界条件，有时还有初始条件，它们统 称为定解条件； 控制方程和定解条件构成所谓定解问题的微分方程或数学模型， 这种以微分方程形式提出问题的方法称为定解问题的微分方程提法。 为了获得数学模型，必须引入某些前提假设以建构几何模 型、物理模型或力学模型等，它们统称为分析模型。 4．1．1 结构分析模型 对任何复杂事物的研究，出发点都是对事物进行逼真而又可行的理想 化以建构分析模型；而结果的可靠性和实用价值主要取决于确立模型时对各 种控制条件和参数的正确反映。 何谓模型? 待分析的事物称为原型，其理想化的替代物就是模型。 任何模型都是为了某种特定目的而将原型的某些特征信息简缩、提炼而构造 出来的。 原型有各方面的因素和各种层次的特征，而模型只要求反映与某种目的有关 的那些因素和层次。模型成功的关键是必须反映原型事物的主要属性和特 征，而什么是主要属性和特征则与我们所关注的问题有关。 从本质上讲，有限单元法是求解微分方程的数值方法，即 在物理或工程问题的数学模型之基础上进行近似计算。因此， 有限元计算的精度并不意味着实际问题求解的精度。在采用有 限单元法解题时，必须时刻牢记：问题的分析模型具有根本的 重要性。 4．1．2 微分方程的形式 连续介质问题的分析方法是：首先从介质中取微元进行分析， 建立控制方程；然后结合具体的定解条件(边界条件和初始条件)求 解控制方程。显然，问题的物理实质不同，控制方程和定解条件也 就不同。然而，它们可被一般地表示为(图2．2) 在 内 在 内 本待求解的未知函数u可以是标量场(例如温度)，也可以是 若干变量组成的向量场(例如位移、应力)。A和B为对于独立变 量(例如空间坐标)的微分算子。 上述微分方程可以是单个方程，也可以是一组方程。 如直角坐标系下弹性静力问题的控制方程和边界条件，其建立 方法可参考弹性力学教科书。 4．2 泛函变分提法 4．2．1 等效积分形式 现在来研究微分方程的等效积分形式。由于控制方程在域内每一 点都必须为零，因此有: (a) 其中， 是函数向量，它是一组同微分方程个数相等的任意 可积(在 内)函数。 不难证明，式(a)是与微分方程组: 完全等效的积分形式，即式(a)成立，则A(u)式成立。假设A(u)在 内 不处处为零，则由于可以任意选择，我们选为处处与A(u)同号的函数， 于是 在域内恒正。可见，要满足式(a)，必须 在 内 同理，假如边界条件: 在边界上每点都得到满足，对于一组任意可积(在 上)函数 应当成立: 在 内 (b) 这样，就可得到与控制方程和边界条件等效的积分形式 : 统称为微分方程的等效积分形式。 4．2．1 微分方程弱形式 通常情况下，可以对式 进行分部积分，从而得到: 其中，C，D，E，F是微分算子。通常上式称为微分方程的弱形式(weakform)， 相对而言，定解问题的微分方程称为强形式(strong form)。 由于分部积分的缘故，场函数u的导数的阶次在弱形式中比在等效积分形 式中为低。这样，使用弱形式时对场函数便只要求较低阶的连续性。当然， 降低对u的连续性要求是以提高 和 的连续性要求为代价的。不过，由 于原来对 和 并无连续性要求，故适当提高其连续性并不困难。 4.3 泛函与变分的基本概念 在工程中常常遇到 z=f(x) 类型的函数，这时因变量z的值由自变量x的值来确定。 但有时我们还会遇到另外一种特殊类型的函数，它的因变量的值 是由一个函数y=y(x)或几个函数(y1(x)，y2(x)，…)来确定。 例如，求平面上任意给定的两点A和B之间曲线的长度L，由数学 分析知识，有 y=y(x)不同，曲线的长度