用matlab进行有限元法编程Jack Chessa2002年10月3日1 引言本文的目的在于给用户在使用matlab编写有限元代码方面遇到的一些问题作一个简短的概述和指引。估计读者对于有限元方法的理论知识已经有了基本的认识,因此,本文主要关注有限元方法的实施方法。本文将使用一个用mablab编写的用于分析静态线性弹性力学问题的有限元代码示例,以给大家说明如何进行有限元法的编程。1.1 符号 为清楚起见,我们在本文中采用以下符号;粗体斜体字体V表示一个空间向量,其等同于在一个点上的位移或速度的空间维度,粗斜体字体d表示一个向量或矩阵的维数是未知数在离散系统即像STI内斯矩阵的系统矩阵,大写标表示节点编号,而小写标一般指沿笛卡尔单位向量的向量组成部分。所以,如果d是节点未知数的系统向量,uI是节点I的一个位移向量,那么 uIi是节点I在方向i上的位移向量即uI*ei。MATLAB的语法通常与数学符号混合使用,本文对此增加了增加了明确的解释。本文用打字机字体,字体指示MATLAB的语法被使用。2 有关编写matlab程序 ma