comsol如何定义狄利克雷边界_有限元分析（FEA）是个什么东东

一、有限元能干什么

二、偏微分方程如何求解(有限元发明以前)

三、有限元法的基本思路

四、有限元法的数学基础——降维

五、如何获得“弱形式”的解

六、二维、三维有限元计算

一、有限元能干什么

通常能通过第二牛顿定律能得出这样的方程：

我们把这种方程称之为常微分方程，即自变量只和时间

有关系，和空间位置没关系，换句话说，常微分方程描述的是单质点的变化规律，如：某个物体在重力作用下做自由落体运动，下落距离随时间变化的规律；火箭在发动机推动下在空间飞行，飞行的轨道等等。常微分方程一般是把研究对象当成一个质点或者刚体，研究整体的运动规律。

值得高兴的是，线性常微分方程相对比较好解决，我们可以通过傅里叶变换(J Pan：傅里叶变换后面的到底有什么小秘密)和拉普拉斯变换(J Pan：从另一个角度看拉普拉斯变换)将微分方程变成代数方程，进而得出准确的解析解。

事实证明，对于繁复纷杂的大自然，只用常微分方程是不够的，因为有些研究对象不能简化成质点，一个典型的例子就是琴弦：

琴弦是一个柔性体，在弹拨的时候每个点振动都是不一样的，一根均匀的弦，其自由振动的数学方程是什么呢？假定

表示

点在

时刻的位移，如果取琴弦中一个微元进行力学分析，就可以得到琴弦应遵守的方程为：

可以发现，这个方程的自变量不仅仅是时间，还有空间坐标，我们把这种方程称之为偏微分方程。也就是说，偏微分方程能描述的连续体的各个点随时间连续变化的情况，本质上就是一种“场”的描述了。

严格来说，自然界的各种现象，都需要用偏微分方程来描述，只不过有些情况下为简单起见，我们忽略了位置项带来的影响，如火箭发射的时候，在强烈的振动下壳体会发生变形的，理论上要完整描述火箭的状态，是需要用偏微分的，但是实际上我们还是把它简化成一个质点来描述，因为这样可以简化很多工作。

当然，这个琴弦除了满足上面的方程，还有其它约束，比如说琴弦的两端要固定在支架上，琴弦内部与外界是通过这个边界联系起来的，边界的变化是会影响琴弦的波动的，我们把这些称之为边界条件，比对于琴弦，边界条件一般为：

；

；

除了边界，琴弦的波动显然还和它的“历史”有关。两根同样的弦，一根在薄刀背的敲击下发生的声音就比较刺耳，另一根在手指的弹拨下就比较和谐，两根弦由于初始时刻的振动情况不一样，后来振动下发出来的声音就不一样，我们把初始时刻的状态称之为初始条件：

；

；

通常边界条件(含边界条件和初始条件)有三类，一类是变量本身的约束，如

，称之为狄利克雷边界条件；一类是变量的导数有约束，如

，称之为诺依曼边界条件；还有第三类就是前两者的混合，成为罗宾边界条件。

遗憾的是，对于偏微分方程，除了少数极简单的情况下能计算出解析解(形式还很复杂)，绝大多数情况下我们还无法很好的求解——有限元就是在这种情况下应运而生，简要来说，有限元法就是为了求解偏微分方程！

二、偏微分方程如何求解(有限元发明以前)

前面我们说了，对于大多数偏微分方程而言，我们是没有办法求得解析解的，只能退而求其次，寻求个近似解，即便是这个近似解的过程也是相当不容易的，方法也很简单、粗暴——试凑法！比如我们现在我们假定一个位移函数，其中包含若干个待定系数：

当然这些函数也不是随意假定的，它需要满足位移边界条件，这个靠

来满足，

主要用来模拟非边界条件的点，它们在边界处取值应该为零，否则和

叠加后会破坏位移边界条件。

为待定系数。

系统在外力的作用下，发生位移，产生变形。位移可以是各种各样的，但必须满足位移的边界条件。但是，满足位移边界条件的位移位移有无穷多组，我们应该选哪一组呢？——那就是当总势能取极小值时的位移，这种现象，我们称之为最小作用量原理，这是自然界中一个普遍成立的原理。

实际我们是怎么做的呢？我们知道，不同的

组合就会获得不同的位移，这样我们就可以列出总势能的表达式，其中只包含待定系数

，由前面我们说的最小势能原理，

只有一组是真实的解，那就是当总势能变分

为零的时候，怎么做到的呢？很简单，将

分别对

求导取零就行了，这个时候偏微分方程组就变成线性方程组了，这种解偏微分的方法称之为利兹法。

如果我们所取的位移不仅满足位移边界条件，而且根据它们求得的应力还满足应力边界条件(不要求满足平衡方程)，这种方法称之为伽辽金方法，这种方法对位移函数的要求较高，但计算量小一些。

上面说的都太抽象，我们举个具体的例子，注意下面的例子不用看懂，只是为了演示传统求解偏微分方程是多么费劲！

我们举一个利兹法的例子，可能大家会更容易理解一些：两端简支的等截面梁，受均匀分布载荷

作用如图所示，试求解梁的挠度

。

首先构造位移试函数(说是构造，其实就是猜) ，

当然也不是瞎猜，这个试函数要满足位移边界条件

；

；

则总的势能为(边看形式复杂，就是按公式)：

整理得到