到底什么是有限单元法？

由于个人的原因最近一直在回顾有限元的一些内容，再加上之前一直觉得自己写作能力有很大的问题，也算是自己给自己一个总结，就想非常认真的写一篇有限元的科普文章，为广大结构工作者能在使用有限元程序的时候能更加好的理解软件到底在算些什么，顺便看看自己的写作能力还有没有救，也希望广泛地接受各位大神给我提点建议或者说一起讨论讨论。

在我学习有限元的过程中，也是参考了各行各业或者不同网站的资源，首先我觉得何晓明教授上传到B站上的《有限元基础编程》是非常不错的一个系列课程，但是需要你有一点有限元和Matlab以及编程的基础。

知乎上的Dr. Stein大神的文章也是非常不错的一个有限元资料，但是他讲的非常数学也非常专业，需要多读读。

再就是如果能翻墙的话Allan Bower 的网站也是非常好的资料并且给了不同版本的代码，Brown University毕业的学生跟ABAQUS软件有千丝万缕的联系（当然也是听说），去年上课的时候教授是Brown University毕业的大神，他说比较早期的ABAQUS是Brown University的学生参与了早期的设计，因此编程思路上跟Allan Bower教授网站上给的代码非常相似，包括形成总刚度矩阵的方式等。

如果真的想好好理解有限单元法，最好再去看看数值分析或者微分方程相关的基础知识。后面如果有什么想起来比较好的资源在补充。

首先，从本质上来说，有限元是一种求近似解的思路，最简单的明了的了解就是对于定义域在[0,π]上，f(x)=sin(x)是一条曲线，如果想求f(x)在定义域上的面积，最简单的办法就是积分。

但是如果我不知道具体的函数表达式，仅仅知道若干个点在定义域上所对应的函数值（f(0)=0, …,f(π/2) =1, …, f(π)=0）, 只要知道的点的个数足够多，那么我也可以把不同的点连成直线然后求面积，也可以得到一个近似的积分值。

那么，有限元到底是怎么应用到结构分析上的，在我看来首先要理解一个应力（stress），应变（strain），位移（displacement用u来表示）以及位移的导数

之间的关系。

用最简单的例子来解释，假设一个轴向受力单元（1D-bar element），材料的弹性模量是E，截面面积是A，受作用力P之前的位置是x1，x2，受力以后的位置是x1’，x2’，如图1：

△  图1

待解的未知量u1=x1’-x1，u2=x2’-x2，由于作用了一个力在单元上，单元会有一个拉长Δu = (x2’-x2)- (x1’-x1)，再带入材料力学（我隐约记得是第二章）给出的应变的表达式：

假设这个单元非常非常小，也就是对等式的右边取极限，就可以将应变转换为导数的形式：

在一维线弹性单元，应力应变之间遵循胡克定律：

再将上式稍微调整一下位置，再把P换成一般常用的字母小f来表示，Govern equation（实在不知道应该怎么翻译）可以表达为：