现阶段结构力学分析的江湖几乎已被有限单元法一统江山。但是实际上，有限单元法只是数值方法求解微分方程的众多手法之一，在此文中，我将介绍现阶段几位一样能够参与结构分析华山论剑的绝顶高手。
盟主的身世——有限单元法的发展
有限单元法思想最早出现在1943年，Courant采用定义在一系列三角形区域的分片连续函数结合最小势能原理求解St.Venant扭转问题。由于当时的计算工具的限制，这种方法并没有在接下来的时间里得到较好的发展。直到1960年以后，电子计算机的广泛应用和发展，现在意义上的有限单元法才真正发展起来。
现代意义上的有限单元法第一个成功的尝试出现在1956年，这是由Tuner和Clough等人在分析飞机结构时提出的。他们使用有限单元法三角形单元得出了平面应力问题的正确解答。1960年Clough将其进一步推广到平面弹性问题，并第一次提出“有限单元法”这一名称。而在之后的1965年O. C. Zienkiewicz和Y. K. Cheung发现存在变分形式的所有场问题即满足变分原理的所有场问题，都可以用与固体力学有限元法相同的步骤求解。
但是并不是所有问题都能找到对应的变分原理，在不能找到变分原理的问题中应用有限单元法就必须找到能够给出有限元格式的方法。在1969年，Oden J T从所谓的能量平衡法出发，得出热弹性问题的有限元分析方程组。在1969年B. A. Szabo和G. C. Lee将加权余量法，特别是Galerkin法引入有限单元法导出标准的有限元过程来求解非结构问题。加权余量法引入有限单元法后使得有限单元法的应用范围得到了极大拓展。而且有限单元法同门中还有两位骨骼惊奇的奇侠，一般的有限单元法，形函数都是一阶线性的，对于复杂的几何边界条件比较吃力，所以有了等参元法和Wilson非协调元，这两位形函数中含二次项的选手。
20世纪70年代，有限单元法的发展已达到高峰，其基本理论和方法已趋成熟。
盟主的招式——有限单元法的基本步骤
盟主有限单元法在解决敌手时总有个固定的套路。
起手式，将求解域离散为一系列单元。就是现在各路软件中划分 格的过程。
紧接一招隔空猜物，就是在各个单元内用基于单元各结点值的单元插值函数来拟合此单元上的未知场函数，就我们常见的力学求解问题，一般都是先求解整体结构的位移函数。
然后一招万法归宗，根据力学原理求出各单元的刚度矩阵。再由单元刚度矩阵凝聚为总刚矩阵，配合边界条件解出数值解。
在一般的结构力学问题中，这几招下来基本无坚不摧。
盟主的罩门——有限单元法的瑕疵
盟主虽然已然一统江湖，但是其余高手仍未被除尽，可见盟主武功中也存在弱点。
首先就是单元的划分，虽然许多有限元程序有了自动划分单元的功能，但是至今还无法保证划分的 格一定合理，由于 格划分不合理导致计算精度差甚至不收敛也不能完全避免。
其次，现有的有限单元法大都采用一阶线性形函数，在拟合复杂问题时效率不高，而等参元法拟合的二次形函数并不完备，Wilson非协调元由于不符合变形协调条件，不能保证收敛。
再者，有限单元法在涉及大变形问题时，由于计算中 格会出现严重的扭曲，这样不但需要 格的重构，而且也会导致计算精度的下降。在处理裂纹发展问题时，由于裂纹发展的不确定性，在计算中随着裂纹的发展不得不不断地重新划分单元来模拟整个裂纹发展的动态过程。
二、隐侠—— 无 格法
隐侠的由来——无 格法的发展
针对盟主的罩门，另一位高手修炼了另一门绝学——无 格法，其计算与 格无关。这种无 格法的出现，在大变形和裂缝发展问题中，不像有限单元法一样需要重新划分 格，因此无 格法在处理这类问题时不仅能够提高计算精度，并且还能减小计算难度。
无 格法的出现最早可以追溯到20世纪70年代。当时人们开始研究对非规则 格有限元差分法的研究，这是无 格法的前身。1977年，Lucy提出了“Smoothed Particle Hydrodynamics”方法，此方法即被简称为SPH方法。在1995年Liu等人在其基础上提出了再生核粒子方法——“Reproducing Kernel Particle Method”简称为RKPM方法。RKPM方法和SPH方法一样都是基于核近似的无 格方法，它在SPH方法基础上引入修正函数施加再生条件并且采用高斯积分。RKPM方法不仅解决了SPH方法在边界条件上不一致性的问题，而且完全消除了SPH方法的张力不稳定。但是由于RKPM方法不满足Kronecker delta 函数性质，这使得这种方法施加边界条件比较困难。
除了基于核近似的无 格法以外，还有采用移动最小二乘法近似的移动最小二乘法——“Moving Least Square”简称为MLS方法。在1992年Nayroles等人将Galerkin法引入移动最小二乘法，从而提出了弥散单元法——“Diffuse Element Method”简称DEM方法。此后的1994年Belytschko等人对弥散单元法进行了改进，在形函数的导数中保留了被弥散单元法省去的项，并且引入Lagrange乘子法应用边界条件，这就是无 格Galerkin法——“The element-free Galerkin Method”简称为EFG方法。其后又有诸多如有限点法——“The Finite Point Method”简称FPM方法、无 格局部Petrov-Galerkin法——Meshless Local Petrov- Garlerkin Method简称MLPG方法、再生核阶谱单位分解方法等方法。
隐侠的招式——无 格法的基本思想
鄙人对无 格法涉猎不深，从基本思想上来讲无 格法是在求解域上布置有限数量的结点，在每个结点上定义覆盖结点周围一定区域的结点权函数或称结点核函数来拟合改结点附近某一区域内的物理量。每个结点所影响的区域称为其支撑域(支撑域可部分互相重合)，由所有结点的支撑域来覆盖整个求解域。这样通过各个结点的权函数的叠加就可以拟合待求解函数。
为何归隐——无 格法的缺点
隐侠无 格法一身奇绝武艺，甚至摆脱了有限单元法 格的桎枯，可以无视 格划分是否合理达到无招胜有招境界，可是为何现实中却无一款商业软件采用此法呢根结底是由于隐侠招法过于奇诡，难以驾驭。就以我们熟悉的结构分析为例，我们已经习惯于在模型的一些的节点上愉快的施加边界条件。给一个节点刚接约束节点就不能有位移的畅快体验在无 格法中几乎成为泡影。因为无 格法中节点的支撑域重合，通常不满足Kronecker函数特性，即一个节点上的位移并不仅由此节点的形函数确定，周围的节点都可能会对其有影响，所以就算要施加一个位移为0的约束也是千难万难，而且无 格法的积分计算较有限元法复杂，所以这位高手暂时只能在特定领域内栖身，江湖中往往仅仅只留下他的传说。
三、盟主的觊觎者—— 广义有限单元法
本部分是半弥散单元法展示实力的部分。
本节中的符 约定：
T3——表示传统三结点三角形单元的有限单元法。
本节以数值算例比较半弥散单元法与传统三结点三角形单元在稳态热传导问题计算过程中的优劣。
四边边界温度已知的矩形温度场问题
如下图所示一无限长矩形柱体的横断面，边长分别为L1与L2，材料性质均匀，且柱体内部无热源。由于柱体无限长，而边界上热交换条件与纵向坐标无关。因此问题简化为二维稳态热传导问题。
边界条件为
半弥散单元法和传统有限单元法在计算中采用同样的 格划分。 格划分如下图所示。
采用以上四种 格划分方法，分别采用半弥散单元法和传统有限单元法求解该问题，求得的考察点P温度值与用式(3)求得的考察点P温度值解析解列于下表。
如以自由度数作为比较标准。在表2中列出了半弥散单元法和传统有限单元法自由度数对比。
可以看出半弥散单元法无论是在相同 格还是在相同自由度数目下，在计算的收敛速度和精度方面都远远优于传统有限单元法。
四、乱点英雄谱
观今时结构分析手段之江湖，豪杰并起，百舸争流。而有限单元法凭借其成熟的理论，严整的招法，为各大软件开发商所喜，因之一统江湖，一时之间风光无两。但其严整的招法既有好处也带来了缺陷，对于一些大变形问题，由于对手招式过于诡奇，应付起来难免措手不及。而无 格法招式轻灵多变，应付有限单元法所不能解决的对手时长袖善舞，应付自如，但也因自身套路难以琢磨而难以为众人所用。广义有限单元法则有些集二者之长的意味，但作为新晋高手，至今仍未显大名于世。
本人学识粗鄙，天下英雄实不能尽识，如有错漏，还望各位看官海涵。
参考文献
[1] 王勖成.有限单元法.北京：清华大学出版 ，2003
[2] 陈晓珞. 二次广义有限元性质及其边界条件的处理方法：[硕士学位论文]. 上海：同济大学，2011
[3] 李刚. 半弥散单元法：[硕士学位论文]. 上海：同济大学，2007
相关资源：IC卡破解软件破解红门、捷顺等一切停车场软件-交通工具类资源…
声明：本站部分文章及图片源自用户投稿，如本站任何资料有侵权请您尽早请联系jinwei@zod.com.cn进行处理,非常感谢！