网格划分的原则.pdfVIP

划分网格是建立有限元模型的一个重要环节，它要求考虑的问题较多，需要的工作量较大， 所划分的网格形式对计算精度和计算规模将产生直接影响。 为建立正确、 合理的有限元模型， 这里介绍划分网格时应考虑的一些基本原则。 1 网格数量 网格数量的多少将影响计算结果的精度和计算规模的大小。 一般来讲， 网格数量增加， 计算 精度会有所提高， 但同时计算规模也会增加， 所以在确定网格数量时应权衡两个因数综合考 虑。 图 1 中的曲线 1 表示结构中的位移随网格数量收敛的一般曲线，曲线 2 代表计算时间随网 格数量的变化。 可以看出， 网格较少时增加网格数量可以使计算精度明显提高， 而计算时间 不会有大的增加。 当网格数量增加到一定程度后， 再继续增加网格时精度提高甚微， 而计算 时间却有大幅度增加。 所以应注意增加网格的经济性。 实际应用时可以比较两种网格划分的 计算结果，如果两次计算结果相差较大，可以继续增加网格，相反则停止计算。 图 1 位移精度和计算时间随网格数量的变化 在决定网格数量时应考虑分析数据的类型。在静力分析时，如果仅仅是计算结构的变形， 网格数量可以少一些。如果需要计算应力，则在精度要求相同的情况下应取相对较多的网 格。同样在响应计算中，计算应力响应所取的网格数应比计算位移响应多。 在计算结构固 有动力特性时， 若仅仅是计算少数低阶模态， 可以选择较少的网格， 如果计算的模态阶次较 高，则应选择较多的网格。 在热分析中， 结构内部的温度梯度不大， 不需要大量的内部单元， 这时可划分较少的网格。 2 网格疏密 网格疏密是指在结构不同部位采用大小不同的网格， 这是为了适应计算数据的分布特点。 在 计算数据变化梯度较大的部位 (如应力集中处 )，为了较好地反映数据变化规律，需要采用 比较密集的网格。 而在计算数据变化梯度较小的部位， 为减小模型规模， 则应划分相对稀疏 的网格。这样，整个结构便表现出疏密不同的网格划分形式。 图 2 是中心带圆孔方板的四分之一模型，其网格反映了疏密不同的划分原则。小圆孔附近 存在应力集中，采用了比较密的网格。板的四周应力梯度较小，网格分得较稀。 其中图 b 中网格疏密相差更大，它比图 a 中的网格少 48 个，但计算出的孔缘最大应力相差 1% ，而 计算时间却减小了 36% 。由此可见，采用疏密不同的网格划分，既可以保持相当的计算精 度，又可使网格数量减小。因此， 网格数量应增加到结构的关键部位，在次要部位增加网 格是不必要的 ，也是不经济的。 图 2 带孔方板的四分之一模型 划分疏密不同的网格主要用于应力分析 (包括静应力和动应力 )，而计算固有特性时则趋于采 用较均匀的钢格形式。 这是因为固有频率和振型主要取决于结构质量分布和刚度分布， 不存 在类似应力集中的现象，采用均匀网格可使结构刚度矩阵和质量矩阵的元素不致相差太大， 可减小数值计算误差。同样，在结构温度场计算中也趋于采用均匀网格。 3 单元阶次 许多单元都具有线性、 二次和三次等形式， 其中二次和三次形式的单元称为高阶单元。 选用 高阶单元可提高计算精度， 因为高阶单元的曲线或曲面边界能够更好地逼近结构的曲线和曲 面边界， 且高次插值函数可更高精度地逼近复杂场函数， 所以 当结构形状不