用matlab进行有限元法编程.docVIP

用matlab进行有限元法编程 Jack Chessa 2002年10月3日 1 引言 本文的目的在于给用户在使用matlab编写有限元代码方面遇到的一些问题作一个简短的概述和指引。估计读者对于有限元方法的理论知识已经有了基本的认识，因此，本文主要关注有限元方法的实施方法。本文将使用一个用mablab编写的用于分析静态线性弹性力学问题的有限元代码示例，以给大家说明如何进行有限元法的编程。 1.1 符号 为清楚起见，我们在本文中采用以下符号；粗体斜体字体V表示一个空间向量，其等同于在一个点上的位移或速度的空间维度，粗斜体字体d表示一个向量或矩阵的维数是未知数在离散系统即像STI内斯矩阵的系统矩阵，大写标表示节点编号，而小写标一般指沿笛卡尔单位向量的向量组成部分。所以，如果d是节点未知数的系统向量，uI是节点I的一个位移向量，那么 uIi是节点I在方向i上的位移向量即uI*ei。MATLAB的语法通常与数学符号混合使用，本文对此增加了增加了明确的解释。本文用打字机字体，字体指示MATLAB的语法被使用。 2 有关编写matlab程序 matlab编程语言在说明如何进行有限元法编程方便十分有用，因为其允许我们进行非常迅速的编写数值计算，并且其具有广阔的prede NED数学库。另外，矩阵（疏与密），矩阵和许多线性代数工具极影呗定义了，开发人员可以完全集中在没有定义这些数