有限元分析可能会出现什么错误.docVIP

有限元分析可能会出现什么错误？ 当有限元分析扩散向可能没有正式数字化程序培训的设计者的时候，专业人员必须问“最适当的方法是否被采用？这些方法是否产生了精确的结果？” 在今天的设计领域，有限元方法被广泛的应用，其中包括各种各样的通用商业软件和适合某专业领域的专业软件。这些方法日益增长的被用，在帮助确定好的新的设计的同时改良了设计性能和成本。 图一 应变能随单元尺寸变化图 考虑到有限元方法在各个设计领域起着重要的作用，专业人员需要问他们自己“他们的设计程序是否是可获得的最适合的技术？这些方法是否会产生精确的结果？”。这些问题是非常重要的，因为越来越多的设计人员不见得受过数字化程序培训，但是他们却在他们的工作中应用有限元方法。 右图是应变能量随单元尺寸减小的变化示意图。从图中可以看出应变能会随着单元尺寸的减小而收敛。 当这些有限元方法被向越来越多且越来越广的群体广泛的应用的时候，用户必须问有限元分析会出现什么样的错误。本文目的不是在广义上解决这个问题，更恰当的说，我们必须集中焦点于有限元方法的可靠性和准确应用方面。为了便于说明，我们考虑线弹性问题，假设有限元的代数方程精确地被求解。对于复杂的分析，考虑这些条件的同时，还有一些额外的要求也有必要得到。 数学模型 首先，设计人员应该记住有限元方法是为了求解数学模型，这数学模型是实际物理问题的理想化结果。数学模型是建立在考虑几何、材料特性、加载条件和位移边界条件等假设的基础上的。数学模型的指导方程是考虑到边界条件的偏微分方程。这些方程不能用封闭的解析方式求解，因此，设计人员要借助有限元方法获得一个数值解。 例如，考虑一个几何和载荷为轴对称的阀套。在这种条件下，考虑轴对称分析条件是合理的。分析的数学模型可以通过指定几何尺寸、支撑条件、材料常数和加载条件来获得。 虽然通常情况下设计人员不能用封闭方程的数学模型的精确解，但是这个数学模型的精确解是存在的，且是唯一的。高精度的精确解的近似解可以用有限元方法求得。 为了充分理解这些观察到的，必???要有收敛概念。这里E表示应变能的数学模型的精确解（未知的），Eh表示对应于单元尺寸h的应变能的有限元解。那么收敛表示为： 上述图一中的示意图表明了收敛是如何达到。当单元网格趋于精确时（这就意味着单元尺寸h的减少），应变能Eh将趋向于数值E。E和Eh之间的误差值的减少速度视解答的题目，也依靠采用单元类型和网格质量。明显地，网格细化过程中高阶单元减少误差率比低阶单元快。 可靠性问题 有限元方法的可靠性是指，在提出很好的数学模型求解时，有限元程序有两个特性。第一，在任何材料特性、位移边界条件和加载条件下，当单元尺寸h趋向0时，有限元的结果会向着数学模型的精确解收敛，合理的有限元网格可以得到一个合理的有限元答案。此外，当材料特性（或板壳的厚度）变化时，有限元结果的质量不会发生彻底地剧烈变化。 上述条件时非常重要的。如果违背了第一个收敛性条件，这变细的网格所得到的有限元分析渐进的结果也不是数学模型的精确结果。这样的误差结果可能导致错误设计，随之带来不幸的后果。当然，违背收敛性的有限元方法也不应该被采用。 为了考虑第二个条件，即用合理的网格得到合理的结果，假设阀套采用的是钢铁材料（杨氏模量是2E5Mpa，泊松比是0.30）。用合理的网格能得到一个能够接受的结果，即误差|E—Eh|是可以接受的小。假如我们现在把材料变为塑料的，材料的泊松比是0.49（接近不能压缩情况0.50）。材料情况的变化会导致精确解相对小的变化，那么它也应该导致相应的有限元结果的一个小的变化。不幸的是，现用的有限元公式违背了第二个条件，这使得当泊松比变化到0.49时，有限元结果会产生一个非常大的误差。 这解答现象是遵守“以位移为基础的有限元”。大误差的存在是因为当泊松比v接近0.50时，单元的刚度太大，，所以当V＝0.50时，不能被用的。这个问题出现的力学原因是：考虑到应力p=Kev，p为压力，K为体积模量，ev为体积应变，当v趋近于0.50时，K变得非常的大，当V＝0.50时，它是无限大。再者，在精确解中，是当n趋近于0.50时，ev变得非常的小，当V＝0.50时，它为零。因此，在一个几乎不能压缩的分析中，压力是被赋予一个大的数（体积模量）和一个很小的数（体积应变）相乘，应力必须被精确计算来平衡这种表面力。当体积模量变的非常大时，压力一直是个有限数，因此当v趋于0.50时它通常的变化是不大。 作为“位移为基础的有限元方法”补救，“减缩积分”被采用。这意味着在单元刚度矩阵的数值积分中，精确的矩阵是不能被计算的。这个方法是简单地用于规划，它为了用少的计算时间建立单元刚度矩阵，常常通过经验来获得可以接受的结果。然而，这个方法也会导致非常大的误差。 考虑后续章节中提到的支架的频率分析。在这个分析中，采用九