结构力学数值方法：解析法：有限元方法原理.pdf

结构力学数值方法：解析法：有限元方法原理

1绪论

1.1有限元方法的历史与发展

有限元方法（FiniteElementMethod,FEM）起源于20世纪40年代末，最

初由工程师们在解决结构力学问题时提出。1943年，R.Courant在解决弹性力

学问题时，首次提出了使用分片多项式逼近函数的概念，这被认为是有限元方

法的雏形。然而，直到1956年，当O.C.Zienkiewicz和Y.K.Cheung在《工程计

算》杂志上发表了一篇关于有限元方法的文章后，这一方法才开始被广泛认识

和应用。

1960年代，随着计算机技术的发展，有限元方法得到了迅速的推广和应用，

特别是在航空航天、汽车、土木工程等领域。1970年代，有限元软件开始商业

化，如NASTRAN、ANSYS等，这使得有限元方法的应用更加普及。进入21世纪，

有限元方法不仅在工程领域，还在生物医学、材料科学、环境科学等众多领域

发挥着重要作用。

1.2有限元方法的基本概念

有限元方法是一种数值求解偏微分方程的通用技术，尤其适用于解决复杂

的工程问题。其基本思想是将连续的结构或系统离散化为有限个单元的集合，

每个单元用简单的函数（如线性或二次函数）来近似描述其行为。通过在每个

单元上应用基本的物理定律（如牛顿第二定律、能量守恒定律等），可以建立整

个结构或系统的数学模型。

1.2.1单元与节点

在有限元分析中，结构被划分为多个小的、可管理的单元，这些单元通过

节点连接。节点是单元的边界点，单元的变形和应力状态通过节点的位移和力

来描述。

1.2.2形函数

形函数是用于在单元内部插值节点位移的函数。例如，对于一个线性单元，

形函数可以是线性函数，而对于一个二次单元，形函数可以是二次函数。

1.2.3刚度矩阵与载荷向量

每个单元的物理行为可以用刚度矩阵和载荷向量来描述。刚度矩阵反映了

单元抵抗变形的能力，而载荷向量则表示作用在单元上的外力。通过将所有单

1

元的刚度矩阵和载荷向量组合，可以得到整个结构的总刚度矩阵和总载荷向量。

1.2.4求解过程

有限元方法的求解过程通常包括以下步骤：1.结构离散化：将结构划分为

多个单元。2.选择形函数：确定每个单元的形函数。3.建立单元方程：根据物

理定律，建立每个单元的刚度矩阵和载荷向量。4.组装总方程：将所有单元的

方程组合成一个总方程。5.施加边界条件：根据问题的实际情况，施加边界条

件。6.求解未知数：使用数值方法（如直接求解法、迭代法等）求解未知的节

点位移。7.后处理：分析求解结果，如应力、应变、位移等。

1.2.5示例：一维杆件的有限元分析

假设我们有一个长度为1米的一维杆件，两端固定，受到均匀分布的载荷

作用。我们使用有限元方法来分析杆件的位移和应力。

1.2.5.1离散化

将杆件离散化为两个线性单元，每个单元长度为0.5米。

1.2.5.2形函数

对于线性单元，形函数为：

#形函数定义

defshape_function(x,xi):

线性单元的形函数

:paramx:单元内部的坐标

节点位置（或）

:paramxi:01

:return:形函数值

ifxi==0:

return1-2*x

elifxi==1:

return2*x-1

1.2.5.3刚度矩阵

对于一维杆件，刚度矩阵为：

#刚度矩阵计算

defstiffness_matrix(E,A,L):

计算一维杆件的刚度矩阵

:paramE:材料的弹性模量

2

:paramA:杆件的截面积

:paramL:单元的长度

:return:刚度矩阵

k=E*A/L

returnnp.array([[k,-k],[-k,k]])

1.2.5.4载荷向量

对于均匀分布的载荷，载荷向量为：

#载荷向量计算

defload_vector(q,L):

计算一维杆件的载荷向量

:paramq:均匀分布的载荷

:paramL:单元的