第一章非线性数值分析的基础.docVIP

非线性有限元法课程内容 第一部分 非线性方程组的数值方法 第二部分 材料非线性问题，以弹塑性问题为例 第三部分 几何非线性，以接触问题为例 第一部分非线性方程组的数值方法 对于数理方程问题的数值求解，线性问题最终都归结为现象方程组问题的求解，非线性问题最终都归结为非线性方程组的求解。当然，导致非线性方程组的问题不仅限于数理方程问题。 第一章 非线性数值分析的基础 1-1 非线性方程组 求n个实数，，使得 (1.1) 其中为给定的在n维Euclidian空间的区域D上的n元实值函数。 引进向量记号，令 , , 则(1.1)可以写成 (1.2) 这里f表示定义在而取值于Rn的非线性映象，也称为向量值函数。简记为 注意，目前在论文写作过程中，大都遵循如下的公式书写约定：大写黑体代表矩阵，小写黑体代表向量，其它量不用黑体；至于是否用斜体来书写数学符号则取决于个人爱好。 对于方程(1.2)要讨论的问题有很多，包括 方程是否有解(即解的存在性)； 解是否唯一(即解的存在性)，如果不唯一，有多少解； 如何求出求出(1.2)的解(包括多解情况)； 算法的收敛性与数值稳定性问题； 算法的复杂性问题。 目前，关于(1)、(2)两个问题，尚处在研究阶段。作为工科学生，我们主要关心问题(3)和(4)。 导出非线性方程组的例子。 例1-1 非线性两点边值问题 非线性两点边值问题的一般形式是 (1.3) 其中关于u是非线性的。 当时，(1.3)描述了单摆的强迫振动。 对于(1.3)，若假定f在集合 上二次连续可微，则我们可用差分法计算(1.3)的近似解。取n+1个等步长 用中心差商近似，，代表内部节点号 (1.4) 这时得到(1.3)的离散化方 , (1.5) 其中是近似式(1.4)时产生的误差，当u满足一定的光滑性条件时，我们有： 舍去(1.5)中的误差项，以表示的近似值，则得(1.5)的近似方程组 , ， (1.6) 如果引进n阶矩阵 ， 且设映象由下式定义： 则方程(1.6)可以写成 (1.7) 方程(1.7)就是方程(1.2)的形式，其中。 1-2 非线性映象微分学简介 1-2.1 几个基本概念 n维向量也被称为中的一个点。用表示中的向量的某一泛数，泛数是实数点R上的绝对值在中推广。常用的三种泛数是： (-泛数：； 1-泛数：； 2-泛数： 2-泛数又称为Euclid泛数和谱泛数。 对于任意的这三个泛数有如下关系 和 , , 设A = (aij)是一个m ( n阶矩阵，则将A的泛数定义为 也称为由向量诱导的矩阵泛数。如， 对于给定的，集合 称为点的邻域或开球，用 表示相应的闭球。 点称为D的内点，若存在，使。由所有的内点组成的集合称为D的内部，记为int(D)。若D = int(D)，则称D为开集。 点称为D的聚点或极限点，若对于任意的， D的聚点的全体称为D的闭包，记为。 点称为D的聚点或极限点，若对于任意的， D的聚点的全体称为D的闭包，记为。 设是中的m个向量，由这m个向量的所有线性组合所形成的集合是中的一个子空间，称为的张成空间，记为Span，即 Span = 若对于任何及，有 则称D为凸集。 一般地，若，，称为向量组的凸组合。用归纳法可证，若D是凸集，则D内任意个向量的凸组合仍在D内。 1-2.2 连续性与可微性 一般情形的多元函数(向量值函数)为，即 这里为的第i个分量函数。m = n为其特例。若各分量函数是线性的，即 ， 此时记 , 则可表示为 (1.8) 记所有由线性映象的集合为，则上述矩阵，此时称由(1.8) 定义的映象为仿射映象。 定义1-1 设，如果对任何固定的，恒有 则称f在是半连续的。若有 则称f在是连续的。 有定义可知连续蕴涵着半连续。 f在连续的(((定义：若对于任何，存在，使得对于任何，有 若映象f在D内每一点都连续，则称f在D上连续。 显然，f在(或D上)连续的充要条件是f的每一个分量函数在(或D上)连续。 定义1-2 称映象于D的内点x处G-可导(Gauteaux可导)，如果存在线性映象，使对任何，，有 (1.16) 此时，称A为f在x点处的G-导数，记为。 注意G-导数是按照直线方向来定义的，因此，f在x处G-可导并不能保证f在x处连续，但能保证f按任意直线方向连续，即在x处半连续。 设于处处G-可导，这时，对每一，存在，从而定义了由到的映象，称此映象为f的导映象。 对于任意固定的x，我们来求(1.16) 中的线性映象的表达式，记 ，取，这里是第j个坐标向量，即。由(1.16)显然有 可见 ; , (1.17) 于是，为在点x的Jacob