结构力学数值方法：矩阵位移法：有限元法入门.pdf

结构力学数值方法：矩阵位移法：有限元法入门

1绪论

1.1有限元法的历史和发展

有限元法（FiniteElementMethod,FEM）的起源可以追溯到20世纪40年

代末，由工程师们在解决结构工程问题时首次提出。这一方法的早期发展与航

空工业紧密相关，特别是在飞机结构分析中。随着计算机技术的飞速进步，有

限元法在50年代末至60年代初得到了显著的推动，成为解决复杂工程问题的

有效工具。

1.1.1历史背景

20世纪40年代：有限元法的概念初步形成，主要用于解决结构

力学中的问题。

1956年：工程师Clough在《美国土木工程师学会期刊》上发表

了一篇关于有限元法的文章，标志着这一方法的正式提出。

1960年代：随着计算机的普及，有限元法开始广泛应用于各种工

程领域，包括土木、机械、航空和汽车工业。

1.1.2发展趋势

数值精度提升：从最初的线性单元到高阶单元，有限元法的精度

不断提高。

应用范围扩展：从结构力学扩展到流体力学、热力学、电磁学等

多个物理领域。

软件商业化：如ANSYS、ABAQUS、NASTRAN等专业有限元分析软

件的出现，使得有限元法的应用更加普及和便捷。

1.2有限元法的基本概念和应用范围

有限元法是一种数值分析方法，用于求解偏微分方程，特别适用于工程结

构的分析。它将连续体离散化为有限数量的单元，每个单元用一组节点来表示，

通过在节点上求解未知量来近似整个结构的行为。

1.2.1基本概念

离散化：将连续的结构或系统分解为有限数量的单元。

节点：单元的边界点，用于定义单元的几何形状和物理特性。

单元：结构的最小分析单元，可以是线、面或体。

1

插值函数：用于在单元内部从节点值推导出任意点的值。

刚度矩阵：描述单元在不同载荷下变形特性的矩阵。

载荷向量：表示作用在结构上的外力和边界条件。

1.2.2应用范围

结构力学：分析结构的应力、应变和位移。

热力学：模拟热传导、热对流和热辐射。

流体力学：研究流体的流动特性，如压力、速度和温度分布。

电磁学：分析电磁场的分布和电磁波的传播。

生物医学工程：模拟人体组织的力学行为，如骨骼、肌肉和血管

的力学分析。

1.3示例：简单梁的有限元分析

假设我们有一根简单的梁，长度为4米，两端固定，中间受到1000牛顿的

垂直载荷。我们将使用Python和NumPy库来演示如何进行有限元分析。

importnumpyasnp

#定义梁的属性

length=4.0#梁的长度

E=200e9#弹性模量，单位：帕斯卡

I=0.05**4/12#惯性矩，单位：米^4

num_elements=4#单元数量

num_nodes=num_elements+1#节点数量

#定义节点坐标

node_coords=np.linspace(0,length,num_nodes)

#定义单元连接

element_connectivity=np.array([[0,1],[1,2],[2,3],[3,4]])

#定义载荷

loads=np.zeros(num_nodes)

loads[2]=-1000#中间节点的垂直载荷

#定义边界条件

boundary_conditions=np.zeros(num_nodes)

boundary_conditions[0]=1#第一个节点固定

boundary_conditions[-1]=1#最后一个节点固定

#计算刚度矩阵

defstiffness_matrix(E,I,L):

计算单个梁单元的刚度矩阵

2

k=E*I