我与有限元 by caoer.docVIP

我与有限元??by caoer7 J W+ [ K: c6 N7 ` I+ ~; V. R8 j2 f---------------------------------------某个周末的一天，电脑前坐了12小时完成了此文。/ a6 p- l6 [+ B. v4 ?: p, [/ d2 _+ ^+ f, Z. Q2 j E z) eCATALOG 4 x T9 x; R9 u z% ]# p% M6 X0 J0 a7 G( ] h# tPREFACE_写本文的初衷 L- ( ~0 u7 b$ G0 [% Z j. tCH1_什么是有限元2 Y0 I6 L6 C, {; U ECH2_有限元软件的介绍与比较* Y; x! F$ a1 z! \/ C: @) E7 }0 VCH3_有限元的数值方法: Y# F( D??u2 H3 nCH4_常用固体单元$ ^3 @2 I0 R# x8 E8 a/ V8 yCH5_有限元的求解器 {/ {: d3 f* V5 H$ XCH6_有限元的未来% U3 z+ \; }$ n \6 P* x5 u/ I! ~: W7 }4 Q- ~, oCONCLUDSION_我能为fem做什么PREFACE_写本文的初衷 by caoer3 _ g1 J, F% W??z @/ S4 _2 P. [5 R: G接触有限元时间不算常也不短，从03年到现在吧，中间断断续续的学习有限元，大大小小软件也用了不少，牛人的帖子大师的著作读了一些，最近常来simwe逛逛，发现还有很多学弟学妹初学fem，提出的问题五花八门但大多不对正路，不竟回想起了当年自己学习有限元时情景，孤立无助到处瞎碰，通宵熬夜啃手册也得不出一个所以然来，现在回头想想如果有几个良师益友会是多么幸福的事情，写此文仅仅只为了让后生少走些当年自己走过的弯路，有更多的精力去研究更深层次的问题，而不是在门槛上打转，同时也算是报答simwer多年来陪伴。$ Y% \0 }: t4 U9 P D/ G8 D) d) U: B$ {* C v }本文如有疏漏的地方，还请各位高手补充指正。CH1_什么是有限元　by caoer- X. V9 M4 c% J* U1 q) k; m7 m5 v9 ]% ^1.1 PDE5 J E9 [# d ~3 ^有限元是一种求解问题的数值方法，求解什么问题呢？－－求解PDE(偏微分方程).那么PDE是做什么用的呢？－－描述客观物理世界。我想如果这两个问题搞清楚了也就明白了为什么要用fem,fem可以做那些东西。 PDE可以描述很多物理现象，电磁，流体，换热，diffusion,力学，河床变迁，物种竞争，股票金融，等等等等。。。。乃至整个宇宙，当然也不是所有的物理现象都可以用PDE描述，如微观世界分子原子的运动等等，所以我从来都不建议用有限元方法仿真微观物质现象的原因，这是题外话就不多说了。除了PDE以外，ODE同样也可以描述客观世界，但ODE多用于控制系统，很多线性PDE的解法也都是将PDE转化为ODE来解的，如separating variables方法。4 s) O5 s* J; y3 E. } Z( m y$ X _1.2 求解PDE% N2 E) X7 I2 k5 k# d4 M7 ~. c8 ?3 _/ T1 z. \有了PDE以后，问题是如何求解并得到结果，首先要说明的是不是所有的PDE都有解的，往往有解的PDE才有实际工程意义。对于数值解法，常用的是有限差分，有限元和谱方法，还有蒙特卡罗法。有限差分出现的较早，计算精度相对较高，但是费时，且模型形状必须规则，边界条件处理困难。有限元方法效率高且满足精度要求，边界条件容易处理，得到了广大的应用。谱方法是基于FFT方法的解法，精度高，收敛快，必须使用periodic boundary condition，适合微观尺度的PDE解。蒙特卡罗法不是基于弱解形式的，多用于金融分析，这里还是着重有限元解PDE，顾名思义，有限元将整个计算几何模型划分为很多小的单元（element）,每个单元的含有一定数量的节点(node)，具体单个单元有多少节点，有对应的不同算法与差值方程，拿一个简单的线性4节点平面单元来说，每个单元包含4个节点，每个节点有对应的variable值，比如简单固体力学问题，每个节点就有对应的位移值，热力学问题每个节点就有对应的温度值，等等。然后单元内部的variables就通过差值方法计算得出。* S8 |/ E$ Z# [4 z: M C3 m5 B4 s9 U G. r??}( V7 I1.3 Galerkin approximation 与 weak formation4 O$ B/ f, C R