有限元分析中的一些问题.pdfVIP

有限元分析的一些基本考虑 单元形状对于计算精度的影响 笔者发现，在分析复杂问题时，我们所可能出现的错误，竟然是一些很根本的错误， 这些根本错误是由于对有限元的基本理论理解不清晰而造成的。 鉴于这个原因，笔者决定对一些基本问题（例如单元形状问题，单元大小问题，应力 集中问题等）展开调查，从而形成了一系列文章，本篇文章是这些系列文章中的第一篇。 本篇文章先考虑有限元分析中的第一个基本问题：单元形状问题。 我们知道，单元形状对于有限元分析的结果精度有着重要影响，而对单元形状的衡量 又有着诸多指标，为便于探讨， 这里首先只讨论第一个最基本的指标： 长宽比（四边形单元 的最长尺度与最短尺度之比），而且仅考虑平面单元的长宽比对于计算精度的影响。 为此，我们给出一个成熟的算例。该算例是一根悬臂梁，在其端面施加竖直向下的抛 物线分布载荷，我们现在考察用不同尺度的单元划分该梁时，对于 A 点位移的影响。 这五种不同的划分方式，都使用矩形单元，只不过各单元的长宽比不同。 例如第一种（ 1）AR=1.1 ，就是长宽比接近 1； 第二种（ 2）AR=1.5, 就是长宽比是 1.5.其它类推。 第五种（ 5）AR=24 ，此时单元的长度是宽度的 24 倍。 现在我们看看按照这五种单元划分方式对于 A 点位移的影响， 顺便我们也算出了 B 点 的位移，结果见下表。 我们现在仔细查看一下上表，并分析其含义。 我们先考虑第一行，它是第一种单元划分情况，此时每个单元的长宽比是 1.1，由此我 们计算出 A 点， B 点的垂直位移，可以看到， A 点的竖直位移是 -1.093 英寸，而 B 点的竖 直位移是 -0.346 英寸。 而这两点我们都是可以用弹性力学的方式得到精确解的， 其精确解分 别是 -1.152 以及-0.360. 这样，我们可以得到此时 A 点位移误差的百分比是 [(-1.093)-(-1.152)]/1.152 = 5.2%. 对于其它情况，也采用类似的方式得到 A 点位移误差的百分比。 从上表可以看出来，随着长宽比的增加，位移误差越来越大，竟然大到 56% 。因此， 如果我们是用长宽比为 24 的单元进行划分的话，那么我们的结果可以说是完全错误的。 下面按照上表绘制出一张图，该图从形象的角度表达了上表的含义。 由此可见，长宽比越接近于 1，那么结算结果越精确，越远离 1，则误差越大。 因此我们在进行有限元分析时， 应该尽量保证划分的单元长宽比接近 1，这意味着， 如 果我们使用了四边形单元， 则最好是正方形单元； 如果使用了三角形单元， 则最好是等边三 角形。 当然，对于一个复杂的零件而言，我们很难保证每个单元都满足这些要求，但是，我 们一定要确保， 在我们所关注的地方， 例如应力最大的地方， 单元形状要接近这一点， 否则， 我们得到的解就是不可相信的。 但是上述结果也告诉我们，即便是最好形状的单元（情况 1，长宽比为 1.1 ），结果的 计算精度也不容乐观，其误差达到 5.2%,那么，我们可以得到更高精度的解答吗？ 可以。这需要单元的细分，下一篇博文中将会详细说明这一点。 有限元分析的一些基本考虑 单元大小对于计算精度的影响 有限元分析一定可以得到问题的精确解吗？ 理论上可以证明，如果插值函数使用了