matlab有限元编程求解实例_基于FreeFem++的有限元模拟

1、引言

有限元方法是２０世 纪５０年代伴随电子计算机的诞生，在计算数学和计算工程领域里诞生的一种高效而灵活的计算方法，它将古典变分法与分片多项式插值相结合，易于处理复杂的边值问题，具有有限差分法无可比拟的优越性，广泛应用于求解热传导、电磁场、流体力学等相关问题，已成为当今工程分析和计算中不可缺少的最重要的工具之一。有限元方法的“化整为零、积零为整”的基本思想与并行处理技术的基本原则“分而治之”基本一致，因而具有高度的内在并行性。

有限元基于分片插值函数，相比于有限差分方法，边界更接近于真实情况，误差相对较小。目前主流有限元软件众多，如 Ansys、MSC、Abaqus、COMSOL等,这些商业闭源软件在业界具有良好的口碑靠性较好。但这些软件由于相关模块被封装，难以理解有限元的计算原理,较难对相关模块作进一步的科研开发，限制了教学和科研的使用。除上述软件外,也有不少开源的有限元软件,如OpenFoam, Coder Aster、OOFEM、Emer、 ParaFEM和 Freefem这些软件不但对公众开放源代码,而且可靠性较高,相关模块的说明文档也较丰富,非常适合于科研。其中, Freefem++语句非常简洁,只需要将原方程转换成弱解形式,使用Freefem++很容易编写有限元程序,是最接近有限元的“计算语言”。

2、FreeFem++简介

FreeFem++是一款免费的、开放源码的二维、三维偏微分方程有限元计算软件，由巴黎第六大学研究人员开发，它集成网格生成器、线性方程组