一篇文章入门多物理场有限元（上）.pdfVIP

一篇文章入门多物理场有限元（上）

邓子平

导读

本文主要写给从事CAD/CAE/EDA/CFD等相关的软件研发测试人

员，也可以作为有限元入门文章，对于从事仿真分析的工程师也可以

对软件实现细节有所了解。鉴于篇幅原因，文章分为上下篇。

前言

从《越狱》里的胡克定理说起

看过美剧《越狱》的朋友对其中一个桥段应该印象很深刻，男主

让狱友将事先画好的一张图铺在墙上，然后用打蛋器在固定点打孔，

最后将墙敲穿推到。狱友问这是什么原理，男主说：这是利用了胡克

定律。

高中物理课本对胡克定律的定义是：

弹簧在发生弹性形变时,弹簧的弹力F和弹簧的伸长量(或压缩量)x

成正比。

很显然，如果仅仅用胡克定律无法解释几个小孔就能把墙凿穿。

事实上男主角应该是事先拿到了墙的结构图，算好了砖块的距离，在

砖头的缝隙处打孔。至于说原理的话，一张白纸很难从中间撕开，但

如果纸中间有破损，沿着破损处就很容易撕开，非要扯上一个定律原

理的话，孔边应力集中原理可以很好地解释这种现象。所以男主所说

的胡克定律只不过是搪塞狱友的一个借口~

广义胡克定律

在材料的线弹性范围内，固体的单向拉伸变形与所受的外力成正

比。

高中物理课本对胡克定律的定义可视为广义胡克定律的简化。

从刚体到非刚体

一根柔软的钢尺一端固定，一端悬空，悬空的一端会发生明显变

形。但是高中物理只能告诉我们固定的一端受到向上的力以及力矩，

无法告诉悬空一端的变形大小，这是因为其研究对象的属于刚体范畴，

即研究的对象本身不会发生变形，而现实中绝大部分物体受力时自身

会发生变形，工程力学中使用应力，应变，位移来度量对象的受力情

况，而材料需要引入弹性模量和泊松比。有限元分析的几何对象是物

体本身。

从一维到二维

弹簧是典型的一维单元，扩展到二维需要引入偏微分方程PDE。

拉普拉斯方程和泊松方程是最简单的二维偏微分方程。这两个方程可

以描述静电场，导热，二维波动，以及天体运动等。简单讲，PDE

（PartialDifferentialEquation）是描述这个世界运行规律的一种方

式。从量子原子微观运动，到宏观热，声音，电磁，以及天体万有引

力规律，都有各自的方程来描述。声场对应赫姆霍兹方程，电磁波对

应麦克斯韦方程，流体对应纳维斯托克方程，热传导对应泊松方程，

量子力学对应薛定谔方程等等，而这个方程基本都是PDE。我们所学

的各种自然定律规律，比如牛顿三定律，浮力定律等只不过是常微分

方程的简化，而常微分方程则是PDE的简化。

有限元方法求解的对象是PDE，有限元方法的本质是将求解对象

离散成小的单元，在每个小的单元上用基函数和形函数表征其坐标和

物理量，使函数在其自由度所在的网格单元几何上满足偏微分方程的

解（之所以是几何不是顶点，是因为自由度不一定在顶点）

有限“单元”篇网格

网格是个很广泛的称谓，为了避免歧义，FEM中的网格做如下规

定：

1、网格的英文名为Mesh，非Grid

Grid英文翻译也是网格，但是基本意思为格子，通常指结构化网

格，即四四方方非常规整的网格。

2、网格是用于有限元计算的输入数据

网格对应于有限元的单元数据；线单元，面单元，实体单元都称

之为网格（对，线单元的一条直线也是网格）

3、几何离散成的三角形称之为“面片”，非网格

大部分显示引擎底层都使用三角形来渲染对象，因此2维、3维几

何都需要三角化（称之为“面片化”），比如一个长方体共有6个面，

每个面需要离散成两个三角形，总共12个三角形。这12个三角形我

们称之为12个“面片”，而非网格！一般情况下“面片”质量很差。

“面片”有两个功能：一是用来做显示渲染数据；二是可以作为网格

划分的输入数据，用来生成面网格。

4.网格计算泛指“分布式计算”

“网格计算”（GridComputing）是分布式计算的一种，它研究

如何把一个需要非常巨大的计算能力才能解决的问题分成许多小的部

分，然后把这些部分分配给许多计算机进行处理。这里的网格翻译也

是Gr