第二章 CAE技术的核心.docVIP

第二章 CAE技术的核心—有限元法：概念与原理 上一章简要地讨论了CAE的基本概念，并指出：凡是可以用于产品性能分析、评估与优化的数值仿真方法，都应归属于CAE技术。在当前被广泛应用的CAE方法中，有限元法最有代表性。有限元方法之所以获得广泛的应用，是因为它的数值仿真内涵十分丰富，例如： 静力仿真——其中包括材料非线性、几何非线性、边界条件非线性 动态仿真——其中包括模态、瞬态、频谱 屈曲和失稳仿真——其中包括线性屈曲、非线性屈曲、后屈曲 热仿真——其中包括稳态/瞬态热传导、幅射分析、金属相变分析 失效/破坏仿真——其中包括裂纹扩展、疲劳破坏、内部缺陷 耦合仿真——其中包括流—结构耦合、热—结构耦合、电—磁耦合 接触仿真——其中包括静态接触、动态接触 制造工艺仿真——其中包括锻造、挤压、拉伸、铸造、轧辗、冲压、超塑、蠕变、焊接 自从1943年Courant第一次正式发表他的有限元思想，到20世纪60年代Clough第一次公开提出“有限元单元法”之后的几十年里，有限元法的理论已经相当成熟。对有限元理论有兴趣的读者可以从本书所列的参考专著中深刻领会。在这一章里，讨论的重点不是较深的有限元理论，而是从工程应用的角度介绍有限元法最基本的概念与最核心的基本思想。之所以将重点放在这里，其基本出发点是： 有限元法的理论基础非常坚实，然而并不神秘； 把握住有限元法的基本思想是创建一个品质优良模型的必要条件，“黑箱”模式不可能创建出一个好的计算模型； 对有限元法的理解有助于消化CAE的其它方法。 第一节 平衡状态与微分方程 施加给一个结构的外力不随时间变化或基本不随时间变化时，所研究的对象是一个静平衡问题。反之，当外力随时间明显变化时，所研究的对象则是一个动力问题，然而，不论是静力问题，还是动力问题，从数学角度看，其力学模型的实质就是求解一个或一组微分或偏微分方程。 首先从一个简单的静力平衡的微分方程问题入手。 设有一个承受均匀载荷的小变形悬臂粱如图2-1所示。 图2-1 悬臂梁 已知材料，梁截面等参数。在载荷的作用下，根据材料力学知识，可以推导出它处于平衡状态下的变形方程： （2-1） 式中是待求的梁的变形曲线。（2-1）表明，所研究的梁的平衡状态等价于求解这个微分方程。考虑到悬臂梁的具体边界条件，求解微分方程（2-1）后，即得它处于平衡状态下的变 形 （2-2） 由（2-2），很容易进一步求得梁内任一点的内力与应力。 再举一例。 设有一均匀拉伸下的矩形薄板如图2—2所示，板中开一小孔，求小孔处应力分布，这是经典的小孔应力集中问题。 图2—2 薄板中小孔应力集中问题 根据弹性力学知识，上述问题可以等价于求解下面这样一个带有边界条件的偏微分方程组： （2-3） 式中，，或为板上任一点(x，y)的应力状态，为单位体积上的力，微分方程（2－3）实质也是一个平衡问题。 用求解平面问题的极坐标方式求解（2-3）后，可得处于平衡状态下的该板的应力分布，其中小孔处最大应力为外载荷强度q的3倍。 类似地，不论结构几何多么复杂，也不论载荷等边界条件种类多么繁多，其三维空间体的平衡问题，也归结为求解下面一组偏微分方程（纳维叶方程）： （2-4） 考虑到具体的边界条件后，可求出处于平衡状态的体内任一点应力状态如图2-3所示。 总之，结构受力后，如果不被破坏，必然处于平衡状态。尽管（2-2）、（2-3）、（2-4）那样的描写平衡状态的数学模型形式简洁，由于求解这些微分方程时，必须同时建立相应边界条件的数学模型，而工程实际问题的边界条件又很难数学模型化，所以通过直接求解这些微分方程去研究例如货车的摇枕、侧架，客车的转向架等等工程结构的平衡状态是不可行的。 图2-3 体内任一点的应力状态 那么，可以不可以绕过这个困难另寻办法呢？回答是肯定的，这个办法就是变分原理，以及基于变分原理的有限元法。 第二节 变分原理的基本思想 对工厂的许多设计人员来说，平衡状态易于想象和理解，而变分原理似乎太抽象，难以理解。其实，变分原理并不神秘，它是高等数学微分自变量的延伸。下面从一个经典的最速降线问题引进泛函和变分原理的基本概念。 已知空间中两点A和B，A高于B，要求在这两点间构造一条路径，使得重物m从A沿此路径自由滑下时，从A到B所需时间最小（忽略摩擦力）。由于需要求解的最小滑下时间取决于路径的形状，而路径形状本身即是一个数学上的函数，当取不同函数或路径曲线时，下滑所需时间即不同，因此下滑时间是一个以函数为