1、有限元网格划分的基本原则 杜平安 机械设计与制造 划分网格是建立有限元模型的一个重要环节， 它要求考虑的问题较多， 需要的工作量较大， 所划分的网格形式对计算精度和计算规模将产生直接影响。 为建立正确、 合理 的有限元模型， 这里介绍划分网格时应考虑的一些基本原则。 1 网格数量 网格数量的多少将影响计算结果的精度和计算规模的大小。一般来讲，网格数量增加，计算精度会有所提高，但同时计算规模也会增加，所以在确定网格数量时应权衡两个因数综合考虑。 图 1 中的曲线 1 表示结构 中的位移随网格数量收敛的一般曲线， 曲线 2 代表计算时间 随网格数量的变化。 可以看出， 网格较少时增加网格数量可以使
2、计算精度明显提高， 而计算时间不会有大的增加。 当网格数量增加到一定程度后， 再继续增加网格时精度提高甚微， 而计算时间却有大幅度增加。 所以应注意增加网格的经济性。 实际应用 时可以比较两种网格划分的计算结果，如果两次计算结果相差较大，可以继续增加网格，相反则停止计算。 图 1 位移精度和计算时间随网格数量的变化 在决定网格数量时应考虑分析数据的类型。 在静力分析时， 如果仅仅是计算结构的变形， 网格数量可以少一些。如果需要计算应力， 则在精度要求相同的情况下应取相对较多的网格。 同样在响应计算中， 计算应力响应所取的网格数应比计算位移响应多。 在计算结构固有动力特性时， 若仅仅是计算少数低
3、阶模态， 可以选择较少的网格， 如 果计算的模态阶次较高， 则应选择较多的网格。 在热分析中， 结构内部的温度梯度不大， 不需要大量的内部单元， 这时可划分较少的网格。 2 网格疏密 网格疏密是指在结构不同部位采用大小不同的网格， 这是为了适应计算数据的分布特点。 在计算数据变化梯度较大的部位 (如应力集中处 )，为了较好地反映数据变化规律，需要采用比较密集的网格。而在计算数据变化梯度较小的部位，为减小模型规模，则应划分相对稀疏的网格。这样，整个结构便表现出疏密不同的网格划分形式。 图 2 是中心带圆孔方板的 四分之一模型， 其网格反映了疏密不同的划分原则。 小圆孔附近存在应力集中， 采用了比
4、较密的网格。 板的四周应力梯度较小， 网格分得较稀。 其中图 b 中网格疏密相差更大 ， 它比图 a 中的网格少 48个，但计算出的孔缘最大应力相差 1%，而计算时间却减小了 36%。由此可见，采用疏密不同的网格划分，既可以保持相当的计算精度， 又可使网格数量减小。 因此， 网格数 量应增加到结构的关键部位， 在次要部位增加网格是不必要的，也是不经济的。 1图 2 带孔方板的四分之一模型 划分疏密不同的网格主要用于应力分析 (包括静应力和动应力 )，而计算固有特性时则趋于采用较均匀的钢格形式。 这是因为固有频率和振型主要取决于结构质量分布和刚度分布， 不存在类似应力集中的现象， 采用均匀网格可
5、使结构刚度矩阵和质量矩阵的元素不致相差太大， 可减小数值计算误差。 同样， 在结构温 度场计算中也趋于采用均匀网格。 3 单元阶次 许多单元都具有线性、 二次和三次等形式， 其中二次和三次形式的单元称为高阶单元。 选用高阶单元可提高计算精度， 因为高阶单元的曲线或曲面边界能够更好地逼近结构的曲线和曲面边界， 且高次插值函数可更高精度地逼近复杂场函数，所以当结构形状不规则、应力分布或变形很复杂时可以选用高阶单元。但高阶单元的节点数较多，在网格数量相同的情况下由高阶单元组成的模型规模要大得多，因此在使用时应权衡考虑计算精度和时间。 图 3 是一悬臂梁分别用线 性和二次三角形单元离散时， 其顶端位移
6、随网格数量的收敛情况。 可以看出， 但网格数量较少时， 两种单元的计算精度相差很大， 这时采用低阶单元是不合适的。 当网格数量较多时， 两种单元的精 度相差并不很大， 这时采用高阶单元并不经济。 例如在离散细节时， 由于细节尺寸限制， 要 求细节附近的网格划分很密，这时采用线性单元更合适。 图 3 不同阶次单元的收敛情况 增加网格数量和单元阶次都可以提高计算精度。 因此在精度一定的情况下， 用高阶单元离散结构时应选择适当的网格数量， 太多的网格并不能明显提高计算精度， 反而会使计算时间大大增加。 为了兼顾计算精度和计算量， 同一结构可以采用不同阶次的单元， 即精度要求高的重要部位用高阶单元，
7、精度要求低的次要部位用低阶单元。 不同阶次单元之间或采用特殊的过渡单元连接，或采用多点约束等式连接。 4 网格质量 网格质量是指网格几何形状的合理性。 质量好坏将影响计算精度。 质量太差的网格甚至会中止计算。 直观上看，网格各边或各个内角相差不大、 网格面不过分扭曲、 边节点 位于边界等份点附近的网格质量较好。 网格质量可用细长比、锥度比、内角、翘曲量、拉伸值、边节点位置偏差等指标度量。 划分网格时一般要求网格质量能达到某些指标要求。 在重点研究的结构关键部位， 应保证划分高质量网格， 即2使是个别质量很差的网格也会引起很大的局部误差。 而在结构次要部位， 网格质量可适当降低。 当模型中存在质
8、量很差的网格 (称为畸形网格 )时， 计算过程将无法进行。 图 4 是 三种常见的畸形网格， 其中 a 单元的节点交叉编号， b单元的内角大于 180， c 单元的两对节点重合，网格面积为零。 图 4 几种常见的畸形网格 5 网格分界面和分界点 结构中的一些特殊界面和特殊点应分为网格边界或节点以便定义材料特性、物理特性、载荷和位移约束条件。即应使网格形式满足边界条件特点， 而不应让边界条件来适应网格。 常见的特殊界面和特殊点有材料分界面、 几何尺寸突变面、分布载荷分界线 (点 )、集中载荷作用点和位移约束作用点等。图 5 是具有上述几种界面的结构及其网格划分形式。 图 5 特殊界面和特殊点网格
9、划分 6 位移协调性 位移协调是指单元上的力和力矩能够通过节点传递相邻单元。 为保证位移协调， 一个单元的节点必须同时也是相邻单元的节点， 而不应是内点或边界点。 相邻单元的共有节点具有相同的自由度性质。 否则， 单元之间须用多 点约束等式或约束单元进行约束处理。 图 6 是两种位移不协调 的网格划分， 图 a 中的节点 1 仅属于一个单元， 变形后会产生材料裂缝或重叠。图 b 中的平面单元和梁单元节点的自由度性质不同，粱单元的力矩无法传递到平面单元。 图 6 位移不协调的网格划分 37 网格布局 当结构形状对称时，其网格也应划分对称网格，以使模型表现出相应的对称特性 (如集中质矩阵对称 )。
10、不对称布局会引起一定误差，如在图 7 中，悬臂粱截面相对 y 轴对称，在对称载荷作用下，自由端两对称节点 1、 2 的挠度值本应相等。但若分图 b 所示的不对称网格，计算出的 y1=0.0346,y2=0.0350。若改用图 c 所示的网格，则 y1 和y2 完全相同。 图 7 网格布局对计算结果的影响 8 节点和单元编号 节点和单元的编号影响结构总刚矩阵的带宽和波前数， 因而影响计算时间和存储容量的大小， 因此合理的编号有利于提高计算速度。 但对复杂模型和自动分网而言， 人为确定合理的编号很困难， 目前许多有限元分析软件自带有优化器，网格划分后可进行带宽和波前优化，从而减轻人的劳动强度。 4
