通俗易懂的有限元基础原理

这q10辅酶本是我的一个回答，鉴于反响不错（有10几个赞 ），也想让更多的人看到，故搬过来写成短文，目的是给有限元初学者一个全局视角，了解原理，以便后续学习。

既然说了要通俗易懂，那我就以讲故事的方式来回答吧。

请注意，为逻辑佳能5d4顺畅，国际法渊源答案中的时间顺序并不完全对应有限元的发展历史。

1.最初，人libreoffice们热衷于用各种力学理论严谨地研究结构构件的受力行为。于是开发了以“平衡（Equilibrium）”、“几何（Kinematics）”、“本构（Constitutive）”三大方程为基础分析方法，姑且称为固体力学法。此法是一个边界条件问题（Boundary Value Problem），需要求解微分方程，算出构件的位移后，再推出应力、应变、反力等。其中“平衡”步步高s1方程和“几何”方程几乎都是微分方程（组）。

平面问题三大方程（截自清华大学mooc：有限元分析及应用）

力学方法虽然理论简单直接，结果也最精确。但是！整蛊邻居攻略！！对于略微复杂的结构，上述分析便会遇到复杂的微分方程。。。学过高等数学的同学应该都了解微分方程的恐怖，轻则算上一天得到错误答案，重则头发掉光发现下不了笔。。。总之，微分方程大大降低了使用弹性力学法的效率，甚至让求解成为不可能。

2.在被微分axis方程折磨够了之后，大家终于想起了一句话：惹不起，躲得起。

人们开始改进固体力学法，试图绕开求解微分方程，于是“虚功原理（Virtual Work Principal）”和“最小势能原理”被搬上了舞台（这两玩意儿是等价的）。

那么怎样才能绕开微分方程呢？答案是在计算一开始就去猜结构受力后的位移，比如你觉得最终的位移可能是线性的，那你可以假设位移的表达式为 u(x)=a_0+a_1x ，另一个人觉得位移可能有周期性，那他可以假设 u(x)=b_0+b_1sinx 。这个位移最好满足边界条件（注意，这里是“最好”），之后算出结构在假设位移下的内外力虚功，或者是应变能，令其满足对应的能量原理。干完这些后大家发现，整个计算过程中没有出现微分方程，只需要求一些积分，就能得到位移表达式中的未知数，从而得到完整的位移函数。

虚功原理与最小势能原理（截自清华大学mooc：有限元分析及应用）

好吧，如果讲成这样你还看不懂，你记住下面这句话就好了：人们通过猜位移并引入能量原理的方式终于摆极速骑行脱了微分方程，可以通过积分来分析结构了。

为什么要这样做呢？因为积分好算啊，再复杂的积分在数值积分面前都是弟弟，数值积分中又以Gaussian Interpolation效率最高，现在的有限元软件基本都采用Gaussian，什么积分点、减缩积分等概念就是数值积分中的东西。

这是有限元发展中极为重要的理论突破，学界将其称作：Convert Strong Form（指力学法解微分方程）to Weak Form（指猜位移求积分）。当然，这个方法是有代价的，后文会详细介绍，正是在消除这个代价的过程中诞生了“有限”的思想。

3.摆脱微分束缚这件事可不得了，打了翻身仗的人们都高兴fong了，争艾滋病恐惧症相运用猜位移的方法去研究结构：从简单到复杂，从一维到多维，从轴力构件到任意变形实体。。。到后来，大家发现用传统书写方式进行计算效率贼低，一大堆类似的方程算来算去，好不自在。

于是有人引入了矩阵（Matrix）来改善计算效率，把相似计算写成矩阵形式，使得书写更加简洁，效率也更高。矩阵的引入和有限元的核家禽养殖心原理的关联不大，但作为一种工具实在是太好用了，人们如获至宝地用矩阵重写了分析过程，使得现代有限元中充斥了矩阵运算。后来人们给运算过程中的飞猪抢票一些关键元素起了名字，比如刚满洲八大姓度矩阵（Element Stiffness Matrix）。

红框中为2自由度杆的单元刚度矩阵（截至Quora, www.quora/What-is-sti日本代拍ffness-matrix）红框中为2节点4自由度梁的单元刚度矩阵（截至Quora, www.quora/What-is-stiffness-matrix）

4.前文第二节提到了一个代价：想要用能量原理，就必须得猜结构位移，而工业耗材位移对结构分析极其重要，后续应力，应变，反力全靠位移导出。对于简单的结构，尚能根据经验猜个大概；但是结构稍微好看点，受力复杂点，猜不猜的为什么额头长痘准就全看运气了。运气好，算出来和真实情况八九不离十，运气不好。。。工程中是不允许这种随机的东西出现的（量子学家请走开）

如果你没学过有限元，总学过微积分吧，把任意一条复杂的曲线分成许多段，每一段都视作直线，只要分的足够多，最后直线连接起来就和原来的曲线差不多。类比这个思想，把结构分成很多很多份，去猜每一份上的位移，比如都猜成线性，只要结构被划分的够多，每一份上猜的准不准就无所谓，在满足收敛条件的前提下，算出来的位移就会很靠近真实的位移，这样一来，就大大降低了猜错位移对最终结果的影响。这是有限元发展中又一重要里程碑，学界把这一方法称为：N-ty联想电脑售后维修pe Refinement（对应于另一条思路：P-type Refinement，即如何猜的更准），有限元理论由此基本构建完整。

为方便使用，人们后来制定了若干标准单元，比如平面四边形形单元，空间六面体单元等，每类单元提前给定好假设位移（即位移模式），将结构划分成一定数量的标准单元来分析，即有限单元法。

就这样，有限元诞生了。此后的时间里，在电子计算机的强力支撑下，有限元渐渐从一种运算繁琐辅助工具成为了当代重要的结构分析手段沙龙会。

试函数就是猜出来的函数（截自清华大学mooc：有限元分析及应用）

至此，结构有限元基本原理已经说完。总的来说，有限元就是一种结构计算的近似方法，核心在于：力学的基本分析，能量原理的引入，猜位移与数值积分，离散化分析。

（2022-6-8更新）很多读者天天催更，之前说过本答案达到一定关注度时，我才会化身人肉有限元。 本答案从19年发布的三年多时间里，点赞不到1300，收藏2000出头，喜欢300+，和我预期的关注度不符。所以我随缘更新吧。。

另外有一些读者批评我不够严谨，三年之后回头看我最原始的答案，确有不严谨的地方，不过并非批评里提到的，故做了一些小修改。

我下面把评论区的一些问题都严谨回答一下。（我先把问题都列出来）

问题1：大佬，我问一下，分得越细一定越准吗? （龙格现象会不会出现?为什么?）

答：猜的位移满足完全性（completeness） 和连续性 （C^n conitunity）后，在没有Spurious Mesh size effect的情况下，网格越小越准确，这个现象叫做收敛。你提到的特殊印刷龙格现象是我理解为hourglass mode，此现象本质是使用减缩积分时，单元刚度矩阵会有rank deficiency（实在不知道中文怎么说，秩亏？），无法描述单元所有自由度的变形，一些单元会出现非真实的沙漏状(hourglass)变形。故你可以使用完全积分，在刚度矩阵中人为加入刚度，使用选择减缩积分避免hourglass mode出现。

问题2： 能讲讲有限元基础理论中的最小势能原理和变分原理吗？

答：

问题3：‘h’-type and ‘p’-type refinement？ 我看到的是这样呀 not n-type

答治疗多动症：我在教材里两种都见过，估计要么是两个都在用，要么某个教材里有错别词语。必须承认的是，h氟碳喷涂-type refinement在学术材料中用的更多，我会酌情修改原回答。

问题4：想请问一下答主，如果想达到能写出这篇回答的水平，大概要读过哪些书呢?题主方便推荐个主要书目清单么？

答：

问题5：为什么最小势能原理等价于静力平衡啊？

答：“静力平衡”不准确，且与本文核心内容无关，原文中已删除，weak form是完全等价于strong from的，对于力学有限元一般是用三大方程中的平衡（力或动量）作为strong form，最终形势往往写成等式右边的Internal Force + 动力分项＝等式左边的 External Force。由于动力分项的存在，静力平衡不准确。