1.网格数量
网格数量的多少将影响计算结果的精度和计算规模的大小。一般来讲,网格数量增加,计算精度会有所提高,但同时计算规模也会增加,所以在确定网格数量时应权衡两个因素综合考虑。
图B-1中的曲线1表示结构中的位移u随网格数量收敛的一般曲线,曲线2代表计算时间随网格数量的变化。可以看出,网格较少时增加网格数量可以使计算精度明显提高,而计算时间不会有大的增加。当网格数量增加到一定程度后,再继续增加网格时精度提高甚微,而计算时间却有大幅度增加。所以应注意增加网格的经济性。实际应用时可以比较两种网格划分的计算结果,如果两次计算结果相差较大,可以继续增加网格,相反则停止计算。
图B-1 位移精度和计算时间随网格数量的变化
在决定网格数量时应考虑分析数据的类型。在静力分析时,如果仅仅是计算结构的变形,网格数量可以少一些。如果需要计算应力,则在精度要求相同的情况下应取相对较多的网格。同样在响应计算中,计算应力响应所取的网格数应比计算位移响应多。在计算结构固有动力特性时,若仅仅是计算少数低阶模态,可以选择较少的网格,如果计算的模态阶次较高,则应选择较多的网格。在热分析中,结构内部的温度梯度不大,不需要大量的内部单元,这时可划分较少的网格。
2.网格疏密
网格疏密是指在结构不同部位采用大小不同的网格,这是为了适应计算数据的分布特点。在计算数据变化梯度较大的部位(如应力集中处),为了较好地反映数据变化规律,需要采用比较密集的网格。而在计算数据变化梯度较小的部位,为减小模型规模,则应划分相对稀疏的网格。这样,整个结构便表现出疏密不同的网格划分形式。下面通过实例给出网格疏密对计算精度的影响。
图B-2是中心带圆孔方板的对称模型,其网格划分反映了疏密不同的划分原则。小圆孔附近存在应力集中,采用了比较密的网格。板的四周应力梯度较小,网格分得较稀。其中图B-2在缺口处网格划分较疏;而图B-4在缺口处的网格划分较密。其应力计算结果:图B-4在缺口处的计算精度高于图B-2中的有限元模型计算的结果,如图B-3及图B-5所示。由此可见,不同的地方应该采用不同的网格划分。因此,网格数量应增加到结构的关键部位,在次要部位增加网格是不必要的,也是不经济的。
图B-2 较粗网格的有限元模型
图B-3 网格对应得环向应力云图
图B-4 缺口处较细网格
图B-5 较密网格所得的环向应力云图
划分疏密不同的网格主要用于应力分析(包括静应力和动应力),而计算固有特性时则推荐采用较均匀的网格形式。这是因为固有频率和振型主要取决于结构质量分布和刚度分布,不存在类似应力集中的现象,采用均匀网格可使结构刚度矩阵和质量矩阵的元素不致相差太大,可减小数值计算误差。同样,在结构温度场计算中也推荐采用均匀网格。
3.单元阶次
许多单元都具有线性、二次和三次等形式,其中二次和三次形式的单元称为高阶单元。选用高阶单元可提高计算精度,因为高阶单元的曲线或曲面边界能够更好地逼近结构的曲线和曲面边界,且高次插值函数可更高精度地逼近复杂场函数,所以当结构形状不规则、应力分布或变形很复杂时可以选用高阶单元。但高阶单元的节点数较多,在网格数量相同的情况下由高阶单元组成的模型规模要大得多,因此在使用时应权衡考虑计算精度和时间。(www.daowen.com)
图B-6中的有限元模型采用了8节点的单元,图B-2中的单元采用了4节点的单元,从其计算结果中可以看出,高阶单元在应力集中处即使采用了较粗糙的网格划分,也可以计算得到较精确的应力值,如图B-7所示。因此,在有应力集中和刚度突变的地方,应该采用高阶单元来对其进行网格划分。
图B-6 高阶单元的有限元网格
图B-7 高阶单元的计算结果
增加网格数量和单元阶次都可以提高计算精度。因此在精度一定的情况下,用高阶单元离散结构时应选择适当的网格数量,太多的网格并不能明显提高计算精度,反而会使计算时间大大增加。为了兼顾计算精度和计算量,同一结构可以采用不同阶次的单元,即精度要求高的重要部位用高阶单元,精度要求低的次要部位用低阶单元。不同阶次单元之间或采用特殊的过渡单元连接,或采用多点约束等式连接。
4.网格质量
网格质量是指网格几何形状的合理性。质量好坏将影响计算精度。质量太差的网格甚至会中止计算。直观上看,网格各边或各个内角相差不大、网格面不过分扭曲、边节点位于边界等分点附近的网格质量较好。网格质量可用细长比、锥度比、内角、翘曲量、拉伸值、边节点位置偏差等指标度量。划分网格时一般要求网格质量能达到某些指标要求。在重点研究的结构关键部位,应保证划分高质量网格,即使是个别质量很差的网格也会引起很大的局部误差。而在结构次要部位,网格质量可适当降低。当模型中存在质量很差的网格(称为畸形网格)时,计算过程将无法进行。网格分界面和分界点,结构中的一些特殊界面和特殊点应分为网格边界或节点以便定义材料特性、物理特性、载荷和位移约束条件。即应使网格形式满足边界条件特点,而不应让边界条件来适应网格。常见的特殊界面和特殊点有材料分界面、几何尺寸突变面、分布载荷分界线(点)、集中载荷作用点和位移约束作用点等。
单元的质量和数量对求解结果和求解过程影响较大,如果结构单元全部由等边三角形、正方形、正四面体、立方六面体等单元构成,则求解精度可接近实际值,但由于这种理想情况在实际工程结构中很难做到。因此根据模型的不同特征,设计不同形状种类的网格,有助于改善网格的质量和求解精度。单元质量评价一般可采用以下几个指标。
1)单元的边长比、面积比或体积比以正三角形、正四面体、正六面体为参考基准。理想单元的边长比为1,可接受单元的边长比的范围是:线性单元长宽比小于3,二次单元小于10。对于同形态的单元,线性单元对边长比的敏感性较高阶单元高,非线性比线性分析更敏感。
2)扭曲度:单元面内的扭转和面外的翘曲程度。
3)疏密过渡:网格的疏密主要表现为应力梯度方向和横向过渡情况,应力集中的情况应妥善处理,而对于分析影响较小的局部特征应分析其情况,如外圆角的影响比内圆角的影响小得多。
4)节点编号排布:节点编号对于求解过程中的总体刚度矩阵的元素分布、分析耗时、内存及空间有一定的影响。合理的节点、单元编号有助于利用刚度矩阵对称、带状分布、稀疏矩阵等方法提高求解效率,同时要注意消除重复的节点和单元。
5.位移协调性
位移协调是指单元上的力和力矩能够通过节点传递相邻单元。为保证位移协调,一个单元的节点必须同时也是相邻单元的节点,而不应是内点或边界点。相邻单元的共有节点具有相同的自由度性质。否则,单元之间须用多点约束等式或约束单元进行约束处理。
免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。