CAD和有限元分析.ppt

CAD与有限元分析 有限元法及其应用 有限元法（有限单元法、有限元素法） 求解复杂工程问题的一种近似数值分析方法 有限元网格生成技术 将物体离散成有限个简单单元的组合 是有限元法的主要瓶颈 自动和交互方法 有限元网格的要求 单元类型： 二维：三角形、四边形 三维：四面体、五面体、六面体 基本要求： 几何逼近 单元合法 密度合理 拓扑相容 结构化网格、非结构化网格 二维网格生成技术 拓扑分解法 节点连元法 网格模板法 映射法 几何分解法 Delaunay法 映射法 手工划分为简单区域 建立网格模板（参数空间相应的单位正子单元的网格） 网格模板到欧氏空间中子区域的映射 网格剖分、相容 优点：算法速度快；算法简单； 单元质量好； 密度可以控制； 可生成多种类型的单元； 可以生成结构化网格和非结构化网格。易于实现，网格规整 缺点：手动划分区域、不规则形状不适用 不能满足对网格疏密过渡的要求 拓扑分解法 根据实体元素的拓扑关系，通过连接相应的顶点将多边形分解成一系列的粗糙的三角形 细化 对于复杂形状，首先将区域的洞切开，形成单连通域；然后每次切下一个三角形；直到最后只有三个顶点 原理简单，对剖分实体的拓扑结构依赖强 节点连元法 首先在边界和区域内布节点，然后按照一定规则生成单元 节点生成： 自动 人工 节点生成： 随机法：边界等距生成；内部区域分成一系列的子区域，随机在子区域插入节点； 扫描线法： 单元生成 夹角最大 Delaunay法：最小角和最大；外接圆准则 基线法 Watson法 Lawson法 几何分解法 1）区域递归细分法 分解成凸区域 对每个区域，在边界插入节点 在区域的“最长轴”的中点处加一个分割线划分区域，在分割线上布节点 递归划分两个子区域，直到所得部分均为三角形 2）单元迭代移去法 跟踪算法 划分边界 从边界移去一个三角形 网格前沿法 划分边界 边界向内部移动，形成一组单元 波前法的优点： 节点和单元同时生成，从而可以控制单元形状和尺寸，取得高质量的网格单元，有利于网格的局部加密和网格过渡； 在新生成单元时，要进行大量的相交判断、包含判断，从而进一步的保证了网格单元的质量。 铺砌法 基本思想是：根据所要划分的网格的密度和大小形成初始化边界节点，这些节点形成一条双向链表。外铺砌边界的节点按逆时针方向排列，内铺砌边界的节点按顺时针方向排列，外边界向里一层层铺砌，内边界向外一层层铺砌，铺砌边界一层层更新，直到铺砌边界只剩六个节点，就对六节点封闭环进行封闭处理 铺砌法的优点： 所有网格基本上都是四边形； 边界敏感性。由于铺砌法是从边界开始铺砌的，能保证所生成的网格绝对通过边界的节点； 网格的尺寸转化容易。使用这种方法可以容易地改变网格单元尺寸的大小； 方向无关性。旋转和平移对于网格生成的结果没有影响。 基于栅格法 1）正则栅格法 用正则栅格覆盖区域，删除区域外单元，调整相交单元 2）四叉树法 正方形包围盒 四分包围盒，判断子区域与物体的关系，决定是否继续四分子包围盒 处理相交的小包围盒，逼近物体边界 优点： 适用于任何复杂的二维和三维域网格生成； 几何适应性强； 网格密度便于控制。 缺点： 生成的网格与所选择的初始栅格及其取向有关，网格单元密集的地方，其网格质量欠佳； 模型内部有少量大网格单元存在，内部网格和边界网格相容性的问题需要解决； 边界网格质量差。 三维曲面网格生成 与二维类似 曲面网格品质评价标准 网格的质量评价应该同时考虑结点、单元和整个曲面网格的质量。结点内角大于150小于2250；每个结点周围的单元数为3或4为宜，若大于6为不规则结点。单元网格质量直接反映单元质量的好坏程度，三角形的质量可以用变形因子值(面积与三边平方和之比)来衡量： 基于三角形的变形因子，可以定义四边形的形状参数。一个四边形沿着两条对角线可分为四个三角形，设这四个三角形分别对应四个值，1≥2≥3≥4,则四边形的变形因子可定义为： 值越大，表明四边形质量越好。凹四边形的值小于0，凸四边形的值在0与1之间，矩形的值为最大1，当四边形退化为三角形时为0。当一个三角形与其相邻的三角形合并为四边形时，按值最大原则进行。 锥度（Taper） 它反映四边形网格单元由两对角线形成的四个三角形面积的差异程度 翘曲度（Warpage） 翘曲度主要用翘曲因子和翘曲角度两种方法来表示，它反映了四边形网格单元的扭曲程度，如图2－4所示。翘曲因子为四边形网格单元对角线的最短距离与单元面积之比。翘曲角度为网格单元对角线所分割的两三角形垂直矢量的夹角。 细长比（Aspect Ratio） 细长比为网格单元最长边长与最短边长之比，它映外观边界的差异。对于理想网格单元，该值为1。 目前大体上有两类定义曲面网格质量的方法， 一类是将网格质量定义为所有单元