第十五章 网格划分方法建立几何模型和选择单元类型以后，就应基于几何模型进行分网分网的工 作量大，需要考虑的问题很多，网格形式直接影响结果精度和模型规模，因此分 网是建模过程中最为关键的环节本节首先介绍网格划分的一般原则，然后介绍半自动和自动两种分网方法， 并介绍自适应分网的基本概念和过程第一节 网格划分原则划分网格时一般应考虑以下原则一、 网格数量网格数量又称绝对网格密度，它通过网格的整体和局部尺寸控制网格数量 的多少主要影响以下两个因素：1. 结果精度网格数量增加，结果精度一般会随之提高这是因为：⑴ 网格边界能够更好地逼近几何模型的曲线或曲面边界；⑵ 单元插值函数能够更好地逼近实际函数；⑶ 在应力梯度较大的部位，能够更好地反映应力值的变化但应注意，当网格数量太大时，数值计算的累计误差反而会降低计算精度2. 计算规模网格数量增加，将主要增加以下几个方面的计算时间⑴ 单元形成时间 这部分时间与单元数量直接相关当单元为高阶单元时， 由于计算单元刚度矩阵要进行高斯积分，所以单元形成要占相当大的比例⑵ 求解方程时间 网格数量增加，节点数量会增加，有限元方程的数量增 加，求解方程组的时间将大大增加。
⑶ 网格划分时间 网格数量增加时，无论采用半自动还是自动方法，都会 使网格划分更多的时间由于网格数量增加对结果精度和计算规模都将提高，所以应权衡两个因素综 合考虑一般原则是：首先保证精度要求，当结构不太复杂时尽可能选用适当多 的网格而当结构非常复杂时，为了不时计算精度而又不使网格太多，因采用其 他措施降低模型规模，如子结构法、分布计算法等图 15-1 中的实线表示结构位移随网格数量收敛的一般曲线，虚线代表时间 随网格数量的变化曲线可以看出，当网格数量较少时，增加网格数量可明显提 高精度，而计算时间不会明显增加当网格数量增加到一定程度后（例如点 P ），继续增加网格对精度提高甚微，而计算时间却大幅度增加因此并不是网格分得 越多越好，应该考虑网格增加的经济性实际应用时并不知道划分多少网格最合理、即不能事先确定点P的位置，这 时可先试算一次，然后适当增加网格，再进行计算比较两次计算的结果，如果 结果相差较大，则应继续增加网格；如果结果相差很小，则没必要继续增加结果精度与网格数量的关系因具体分析结构而异一些简单的结构在简单载 荷作用下，变形非常简单，则用少量网格就可以得到很高的精度例如受集中载 荷的等截面悬臂梁，即使用一个梁单元也可得到非常精确的结果。
但对一些复杂 工况下的复杂结构，由于内部位移场分布很复杂，即使采用较多网格，也不一定 能得到满意的结果在选择网格数量时还应考虑分析数据的类型和特点，一般可以遵循以下原 则⑴静力分析时，如果仅仅是计算变形，则网格可以取得较少如果需要计 算应力或应变，若保持相同精度，则应取相对多的网格例如，图15-2中的悬臂梁分别用图(a)所示的三种网格离散，计算出的最大 应力和最大变形如图(b)所示可以看出在网格数量相同的条件下，位移计算精度高于应力精度，或者说在精度相当的条件下，应力计算的网格应该比位移计算 多闯胳】•总力计轩姑应⑵在分析固有特性时，如果仅仅计算少数低阶模态，可以选择较少的网格 如果需要计算高阶模态，由于高阶振型更复杂，所以应选择较多的网格计算的模态阶次越高，要求模型越详细此外，选择网格数量时更应考虑质量矩阵的形 式由于一致质量矩阵的计算精度高于集中质量矩阵，所以在采用一致质量矩阵 计算时可以划分较少的网格，而采用集中质量矩阵时则应选择相对较多的网格图15-3是一梁的固有频率随网格数量的变化，其中实线表示一致质量矩阵, 虚线表示集中质量矩阵可以看出：①采用一致质量矩阵和集中质量矩阵时， 随着网管数量增加，固有频率分别从精确值的上方和下方收敛。
②在相同网格 数量的条件下，采用一致质量矩阵的计算精度高于集中质量③ 网格数量相同 时，高阶频率的计算精度要比低阶频率低或者说，如果要保持各阶固有频率具 有相同的计算精度，那么计算高阶频率时应选择更多的网格1,1*二肚二阶BFJT上 "小一 ■k — L/ Z -f r(J.ES i ///■-f f 二站三惭\SI所L.415a图口一孑PI克鞭率与糾i&sdt曲塔⑶在结构的响应分析中，如果仅仅是计算某些位置的位移响应，则网格数 量可以少一些如果需要计算应力响应，则应选择相对较多的网格⑷ 在热传导分析中，结构内部的温度梯度趋于常数，不需要大量的内部单 元，所以可以划分较少的网格但是如果热应变形和热应力计算采用同热传导分 析的模型，则应根据应力和位移计算的特点选择网管数量二、网格疏密网格疏密是指结构不同部位采用不同大小的网格，又称相对网格密度实际 应力场很少有均匀分布的，或多或少存在不同程度的应力集中为了反映应力场 的局部特性和准确计算最大应力值，应力集中区域就应采用叫密集的网格而在 其他非应力集中区域，由于应力变化梯度小，为减少网格数量，则可采用较稀疏 的网格因此整个结构显示出疏密不同的网格划分。
图15-4是一中心带圆孔的方形板的1/4模型，其网格形式反映了上述原则 即小孔附近存在应力集中，采用了较密的网格，而板的四周应力梯度小，网格相 对较稀其中图(a)网格疏密相差较小，模型共有132个单元而图(b)网格疏 密相差较大，只有84个单元但是通过两者计算出孔缘最大应力分别为 300.60MPa和296.36MPa (理论值为300MPa),误差仅相差1 %,而计算时间可 减少36%•可见，采用疏密不同的网格划分，既可保持相当的精度，又可使网格 数量减少该例说明，计算精度并不随着网格数量增加而绝对增加，网格数量应该增加到结构的关键部位，在次要部位增加网格是不经济的图15-4谢孔方扳的1/4模型采用不同密度的网格划分时，应注意疏密网格之间的额过渡过渡的一般原 则是使网格尺寸突然变化为最小，以避免出现畸形或者质量较差的网格过渡的 常见方式有以下几种：1. 单元过渡单元过渡是用三角形网格过渡不同大小的四边形网格，或用四面体和五面体 网格过渡不同大小的六面体网格图15-5是三角形网格过渡四边形网格情况， 其中图(a)中的四边形网格相差不大，只用了一层三角形网格过渡而图(b)中2. 强制过渡强制过度是指用约束条件保持大小网格之间的位移连续性，这时大小网格的 节点不可能完全重合，大小网格之间具有明显的界面，如图15-6 (a )所示。
为保 证过渡界面位移的连续性，必须对小网格在大网格边上的节点(即节点2、4、6、 8)进行位移约束约束方式有两种：一是利用多点约束等式；一是：通过约束 单元例如节点2,其位移u、u应满足约束等式2 2u + u和u +uu = —1 3u 二 ■2 22 23. 自然过渡单元自然过渡和强制过度均适合于半自动分网方法，而对于自动分网，目前 还不能完全按照人的意愿划分出图15-5或15-6(a)所示的网格形式，尽管这些形 式并不复杂对于一定大小的平面区域或空间体积，目前的自动分网算法只能划 分出具有相同形状和平滑过渡的疏密网格，如图15-6 (b )所示，这种大小网格之间的平滑过渡称为自然过渡自然过渡将引起网格变形，从而降低网格质量网 格尺寸越悬殊，过度距离越近，网格质量影响越严重15-&强制过渡利口憑过渡划分疏密不同的网格主要用于应力分析，包括静应力和动应力而在固有特 性分析时，则应采用比较均匀的网格，这时因为固有频率和振型仅与结构的质量 分布和刚度分布有关，不存在类似应力集中的现象其次均匀网格可使刚度矩阵 和质量矩阵的元素大小不致相差太大，可以减小数值计算误差例如，对于图 15-7 (a)所示的卡子，在计算其应力和固有频率时，就可以采用图15-7 (b)和图 15-7 (c)两种不同的网格形式。
此外，计算温度场时也趋于采用均匀的网格15-7 4 一铝钩制小同网命期廿实际建模中，有时并不能事先知道结构哪些部位存在应力集中，特别是热应 力分布这时可以先用均匀的网格进行初算，在根据得到的应力等值线图来调整 网格疏密等值线密集的部位说明应力梯度大，应增加网格密度等值线稀疏的 部位应力梯度小，可以保持或减少网格密度，然后再作精确计算这种过程还可 以重复进行，以逐步提高计算精度目前一些分析软件可以根据计算结果自动调 整网格疏密三、单元阶次很多单元都有低阶和高阶形式，采用高阶单元的目的是为了提高精度，这主 要基于两点考虑一是利用高阶单元的曲线或曲面边界更好地逼近结构的边界曲 线或曲面；二是利用高次插值函数更好地逼近复杂的实际函数但高阶单元的节 点较多使用时也应权衡精度和规模综合考虑增加网格数量和单元阶次都可以提高结果精度，但在节点总数相同的情况 下，增加阶次的效果更理想图15-8 (a )中的两种单元总是自由度相同，根据式 (12-1)可知，线性单元的应力计算误差为0(h'2),二次单元的应力计算误差为 0 (h )，显然后者比前者的计算误差小而对于图15-8 (b )所示的网格，二次单元 不仅误差小，且节点数还比线性单元少一个。
Lj-E纯性单元和二次单兀的比较图15-9是一悬臂梁分别用线性和二次三角形单元离散时，其顶端位移随网 格数量的收敛情况，从该图不难看出以下几点：⑴当网格数量较少时，两种阶次单元的精度 相差很大，这时采用低阶单元划分网格是不合适 的團巧-9不同阶次单.元■的收敛情况⑵当网格数量较多时，两种阶次单元的精度 相差很小，这时再采用高阶单元并不经济例如在 离散结构细节时，由于细节处的网格分得很密，这 时采用高阶单元的意义就不大了，采用密集的线性 单元可能比稀松的高阶单元效果更好[fl 15-10采用画种阶汰单元的齿轮■模皑⑶在精度一定时，需要的高阶单元数要远远少于线性单元数因此在使用 高阶单元时要选择适当的网格数量，太多的网格并不能带来明显的效益例如在 计算齿轮根部应力时，采用了图15-10所示的二次单元图(a)和线性单元图(b)， 尽管两者网格数量相差近115倍，但计算结果却非常相近有时为了兼顾进度和计算量，可以在不同部位采用不同借此的单元，即精度 要求高的重要部位采用高阶单元，而精度要求低的次要部位采用低阶单元不同 阶次单元之间应注意正确连接，以保证位移连续四、网格质量网格质量是指网格几何形状的合理性。
质量的好坏将直接影响结果的精度， 质量太差的网格甚至会终止计算过程直观上看，若网格各边和各个内角相差不 大，网格表面不过分扭曲，边角点位于边界等分点附近，则这类网格的质量较好网格质量可用一些具体指标定量表示网格划分之后，特别是自动划分的网格，应进行网格质量检查，并对质量差的网格(特别是重要部位的网格)进行修 改，以保证计算精度和使数值计算结果顺利完成在有限元模型中，图15-11 所示的几种网格是不允许的，它们将导致单元刚阵为零或负值，数值计算时将出 现致命错误而中断，这些网格称为畸形网格其中图(a )所示的网格节点交叉编 号，节点必须按顺时针或逆时针统一编号图(b)所示网格的内角大于或等于 180°，图(c )所示网格的两队节点重合，导致网格面积为零五、网格分界面和分界点划分网格时，结构中的一些特殊界面和特殊点应划分为网格边界或节点常 见特殊界面和特殊点有以下几种⑴不同材料的分界面；⑵几何尺寸的突变面，如板壳结构和平面应力结构不同厚度的分界面，杆 件结构不同截面的分界面；。