网格和单元的基本概念.docx

网格和单元的基本概念 前记：首先说明，和一般的有限元或者计算力学的教材不一样，本人也不打算去抄袭别人的著作，下面的连载是一 个阶段的学习或者专业感悟集大成，可以说深入浅出，也可以说浅薄之极——如果你认为浅薄，很好，说明我理解透了， 也祝贺你理解透了！好了，废话少说，书归正传。 无论是 CSD（计算结构力学）、 CTD（计算热力学）还是 CFD（计算流体动力学）——我们统一称 之为工程物理数值计算技术。支撑这个体系的 4 大要素就是：材料本构、网格、边界和荷载（荷载问题可以理解为 数学物理方程的初值问题），当然，如果把求解技术也看作一个要素，则也可以称之为 5 大要素。网格是一门复杂 的边缘学科，是几何拓补学和力学的杂交问题，也是支撑数值计算的前提保证。本番连载不做任何网格理论的探讨（网格理论是纯粹的数学理论），仅限于尽量简单化的应用技术揭秘。 网格出现的思想源于离散化求解思想，离散化把连续求解域离散为若干有限的子区域，分别求解各个子区域的物理变量，各个子区域相邻连续与协调，从而达到整个变量场的协调与连续。离散网格仅仅是物理量的一个“表征符号”，网格是有形的，但被离散对象既可以是有形的（各类固体），也可以是无形的（热传导、气体），最关键的核心在于网格背后隐藏的数学物理列式，因此，简单点说，看得见的网格离散是形式，而看不见的物理量离散才是本质核心。 对计算结构力学问题，网格剖分主要包含几个内容：杆系单元剖分（梁、杆、索、弹簧等）、二 维板壳剖分（曲面或者平面单元）、三维实体剖分（非结构化全六面体网格、四面体网格、金字塔网格、结构化六 面体网格、混合网格等），计算热力学和计算流体动力学的网格绝大部分是三维问题。对于 CAE工程师而言，任何 复杂问题域最终均直接表现为网格的堆砌， 工程师的任务等同于上帝造人的过程， 网格是一个机体， 承载着灵魂 （材 料本构、网格、边界和荷载），求解技术则是一个思维过程。 网格基本要素是由最基本的节点（ node）、单元线（ edge）、单元面（ face ）、单元体（ body） 构成，实质上，线、面、体只不过是为了让网格看起来更加直观，在分析求解过程中，线、面、体本质上并没有起 多大的作用，数值离散的落脚点在节点（ node）上，所有的物理变量均转化为节点变量实现连续和传递。在所有的 CAE环境下，网格的基本要素均可以直接构成，但对于复杂问题而言，这是一个在操作上很难实现的事情，因此，基于几何要素的网格划分技术成为现代网格剖分应用的支点，和网格基本要素完全相同，对应的几何要素分别称之为点（ point ）、线（ curve ）、面（ surface ）和实体（ solid ）。 数值离散求解器是不能识别几何元素的，要对其添加“饲料”，工程师必须对几何元素进行“精 加工”，因此，从这个意义上来说，网格剖分的本质就是把几何要素转换为若干离散的元素组，这些元素组堆砌成 形态上近似逼近原有几何域的简单网格集合体。因此，这里说明了一个网格“加工”质量的基本判别标准——和几 何元素的拟合逼近程度，理论上，越逼近几何元素的网格质量越好，当然，几何逼近只是一个基本的判别标准，网 格质量判别有一系列复杂的标准，后文详细阐述。 本篇将专门解释几个基本概念：点网格；一维线网格；二维三角形面网格、二维四边形面网格； 三维四面体网格 （ tetrahedra ）、三维金字塔单元 （ pyramid ）、五面体单元 （prism ）、三维六面体单元 （ hexahedra ）； 1 结构化网格（ structural grid ）、非结构化网格（ nonstructural grid ）、混合网格（ blend grid ）。需要专门 说明的是，网格（ grid net ）不等同于单元（ element ），单元是基于有限元思想的一个专用名词，而网格则是 网格理论体系下的专用称呼，当网格用于有限元分析的时候，便可以转而称为线单元、面单元、实体单元，因此两 者的本质差别在于，单元具有物理意义，一般具有特定的物理列式，而网格只是纯粹几何意义上的基本元素。 点网格主要针对 CSD中的质量单元和 CTD中的点状热源，通过对几何 point 直接 mesh 生成，应用 较为简单， 属性也较为简单， 一般仅仅包含质量特性或者温度特性。 值得一提的是， 在考虑行波效应的振动分析中， 质量单元常被妙用，作为基底无限大质量块，巧妙地将加速度激励转变为力激励，从而达到可以多点施加激励的作 用。 一维线网格主要针对计算结构力学问题，主要针对基本的桁架（ truss ）、梁（ frame or beam ）、 索（ cable ）、连杆 (link) 和弹簧（ spring ）等工程单元。需要特别说明的是，