第五章仿真系统构造方法及其关键技术.pptVIP

2.新网格的生成 塑性有限元对网格生成法的特别要求： （1）严格控制沿工件和模具接触面上网格的规整性和密度，以提高解的精度和减少网格重分次数； （2）使网格过渡区域最小化，以缩短模拟时间。 四边形四节点网格的生成方法(非结构化网格）: （1）区域递归细分法 基于目标域的几何描述，递归地生成有限元网格。 Loop单元划分法 首先把目标域边界按一定的网格密度离散化，形成由若干节点组成的环，环中节点按逆时针方向排列，其节点数为偶数，然后从环上选择两个合适的节点将它一分为二，同时在分割线上进行布点，所生成的两个子区域被递归分解下去，直到生成四节点或六节点单元。最后再把六节点单元划分成四节点单元，这样就在目标域中生成了完全四边形的网格。 四边形四节点网格的生成方法 （2）单元迭代生成法 这类方法是从实体边界向里一次移去一个单元，迭代执行下去直到剩下一个单元为止。最典型的是波文（paving）法，首先把初始几何边界离散为一系列逆时针方向排列的偶数个点，以这些点为基点生成第一层单元，然后再以新节点为基点向区域生成第二层单元，依次迭代下去，直到四边形网格布满整个区域。 四边形四节点网格的生成方法 （3）基于栅格法 基本思想是将一完全包含目标域的规则网格置于其上，除去落在目标域之外的栅格单元，并对与目标域相交的栅格单元进行调整或裁剪以满足边界条件的要求。 3. 网格质量优化 网格应具备以下条件： （1）所有单元都应该接近理想形状； （2）主要场量变化梯度较大的地方网格应较密； （3）疏密网格之间过渡应自然均匀。 两种优化方法： （1）网格光滑处理：保持网格的拓扑关系不变，将每个内部结合的位置摄动到由相邻节点组成的多边形的质心处，使得每个单元更接近理想形状。 缺点：实质上是进行了网格“局部平均化”，即牺牲一些较好质量的网格来提高一些较差单元的质量。 （2）自适应网格加密：首先用较稀疏的网格进行有限元分析，再根据计算结果通过某种误差指示器来判断哪些区域需要加密，经局部加密后再进行有限元分析。其中关键是如何确定单元的离散误差。 4. 节点编号优化 网格划分完后，编号混乱，尤其是在自适应加密后，带宽相当大。 （1）如果相连节点个数相同，则它们的编号之和随节点的编号增大而增大。如1、3、10和12相连的节点数相同，则它们的和值分别为12、16、34和40，以递增的次序排列； （2）平均值即相连节点的编号之和与个数之比，随节点编号的增大而增大。 （3）跨距值即相连节点的最大、最小节点编号的和，随节点编号的增大而增大。 三维网格重新划分 三维网格重新划分与二维网格重分相比，难度和复杂性大大增加，主要体现在： (1) 进行网格重新划分时变形体的几何重构； (2) 在重建的几何实体中建立新网格。 在三维网格划分中，通常采用四面体和六面体单元，其中六面体单元具有求解精度高，可以更好地逼近求解域的优点，因此在三维有限元模拟中大多采用这种单元。但对于具有复杂空间表面形状的变形体，几何重构有一定的难度，且网格划分不易实现自动化。 六面体网格划分方法 （1）由二维四边形有限元网格通过旋转、扫描和拉伸等几何变换形成六面体网格。这种方法难度较低，易实现。 （2）映射法。先把三维实体交互地分成几个区域，然后用形函数映射技术，把各个区域映射为很多细小的八节点六面体单元。该方法较易实现，但缺点是：当三维物体表面是十分复杂的自由曲面时，该方法逼近的精度不高，且人工分区十分复杂，难以实现自动化。 （3）基于网格填充法的六面体网格划分。这种方法是从实体内部开始，尽可能多地填充规则的长方体或正方体网格，在实体边界上，根据实体边界的具体特征更改网格形状和相互连接关系，使边界上的六面体尽可能准确地逼近物体的边界形状。该方法可实现完全自动化，但其最大的弱点是边界单元质量较差，单元较易畸变。 （4）采用波前法逐层由实体边界向实体内部生成六面体网格。该方法生成的六面体网格质量较高，但由于在三维实体内部各单元边与边、面与面之间的相互关系十分复杂，并且在实体表面布点必须满足一定条件才能生成完全的六面体网格，所以在算法的实现过程中难度较高。 三、新旧网格间的数据转换 有限元网格重新划分后，为了保证分析的连续性，必须把旧网格下与加载历史有关的场量数据转换到新网格上。数据分为两类，一类是单元场量，另一类是节点场量。 单元场量的转换分为以下两步： （1）先将该单元场量转换到旧网格节点上，即将旧网格下单元场量变为节点场量；将旧网格节点的场变量转换到新网格节点上； （2）把新网格节点数据通过插值到各个单元。 数据转换步骤： 旧网格系统内的转换 产生新网格节点值（包含测试、插值计算） 计算单元场量 5-6 初始速度场的确定 用牛顿-拉夫森法迭代求解刚塑性有限元法的非线性方程组时，需要设定一个初始速度场。初始速