有限单元法

我早就想写点有限单元法的东西了,可是不知道该不该在这里写,今天下定了决心想写一下.

其实有限元不是计算机的知识.但是,利用有限元,我们可以让计算机做更多的事,而且,有限元对于软件工程思想的介入也是迫不及待的,至少我这么想.

首先,我先说一下有限元可以做什么吧.有限元是一种数值计算思想.它的理论基础是力学的虚功方程.利用有限单元法,我们能够计算各种建筑在动静荷载下的响应.现在,国外已经有了ANSYS,ABQUES等通用的有限元软件,现在在我们的工程及科研领域应用非常广泛.如果将有限单元法用于诸如三维动画,三维RPG游戏之中,我们将提高动画的质量和动画制作的效率.例如游戏中火焰效果,在地震做用下周围场景的响应情况等等.

其次,我浅谈一下有限单元法的基本原理.有限单元法,就是将连续体如板,壳,体等离散成有限个有简单形状的单元,如三角形,四边形,四面体等,然后对单元进行力学分析,最终得出单元刚度距阵[Ke],然后将每个单元的单元刚度距阵[Ke]整装到整体刚度距阵[K]中;根据连续体受到得荷载情况生成整体荷载向量[P].最后解线性方程组[K][δ]=[P]得出每个单元的位移情况[δ],这样我们就得到了整个连续体的变形和位移情况.这就是有限单元法的整个过程.在这个过程中,各种形状单元的单元刚度距阵[Ke]的推倒不是我们做软件的人的工作,而是力学研究者的工作,我们所要做的就是将这种方法形成我们需要的程序,然后将程序嵌入到我们软件所需要的地方.

最后,我要说一下有限单元法存在的问题.首先,应该是运行速度的问题,这里的运行时间主要消耗在线性方程组[K][δ]=[P]的求解上,不过解决这个问题的方法也是很简单,就是我们在离散连续体时,将我们的单元面积(体积)增大,这样就会大大减小整体刚度距阵[K]的结束从而提高计算时间,当然,这样做会减小计算的精度,但是这个精度可能已经远远超过了我们的要求.其次,应该就是有限单元法的精度的问题.当前有限元计算时都是假设连续体所用的材料在当前的应力状态下处于线弹性状态,并且也没有考虑大变形对刚度距阵的影响,不过随着力学的发展,新的考虑材料非线性和几何非线性的有限元方法日渐成熟,因此,有限元的应用会有一个非常好的空间.