论文部分内容阅读
网格结构有限元模型的选择会对计算精度产生很大的影响。如果将网格结构中的每一根杆件和每一个节点都用细小的实体/壳体单元来离散,该有限元模型的计算精度是很高的。但这种方法会造成巨大的计算量,这是现有的计算机水平所无法接受的。通常将网格结构中的每一根杆件采用一个梁/杆单元来离散,并忽略节点的影响。然而,数值分析和模型试验表明,这种有限元模型往往高估网格结构的承载力和延性,精度较差。因此,本文对网格结构的精细化有限元分析方法进行了深入的研究,该方法可以根据网格结构的受力情况选择单元的类型、尺寸、以及离散区域,在保证计算精度的前提下尽可能降低计算量。当网格结构的载荷较小,杆件和节点都没有发生塑性屈服时,网格结构可以采用梁单元模型来进行分析。为了判断网格结构有限元模型中梁单元的长度和插值函数是否合理并对此进行相应的调整,首先推导了梁在受拉、受压、纯弯三种情况下的挠曲微分方程,以有限元试算得到的梁单元两端节点力为边界条件,求出了梁单元广义力分布场的解析解,然后根据Zienkiewicz-Zhu后验误差估计理论,以该解析解为广义力相对精确解,推导了广义力有限元解和广义力相对精确解的能量范数以确定梁单元的相对误差。在试算过程中,如果网格结构中每个梁单元的有限元解相对误差满足都精度要求,则终止试算过程,否则调整梁单元的插值函数或长度后再进行试算。在节点发生塑性屈服前,节点的刚度是可测的。为了分析节点刚度和体积对结构分析精度的影响,采用与节点半径等长的刚臂来模拟节点的体积,在刚臂与杆件之间设置弹簧单元来模拟节点的刚度,根据经典梁理论,推导了考虑刚臂和节点刚度的梁的单元刚度矩阵;确定了影响节点刚度的主要因素;然后,以焊接球节点为例,根据节点刚度的回归公式,计算出符合施工构造要求的焊接球节点对杆件轴向刚度和弯曲刚度的影响范围。在此基础上,推导了节点体积和刚度对位移解影响的误差估计公式。随着载荷的增大,结构的节点和杆件可能会发生塑性屈服和弯曲破坏,无法继续采用全梁单元有限元模型来分析结构的受力,‘需要采用多种不同类型的单元来离散网格结构。为了考虑杆件弯曲破坏对网壳性能的影响,将杆件中可能发生弯曲破坏的部分作为微观尺度模型采用壳单元离散,杆件的其余部分作为宏观尺度模型采用梁单元离散,给出了梁-壳单元的耦合方法。在进行梁单元模型自适应分析的基础上,根据梁单元的内力分布来确定弯曲破坏可能发生的位置,并以塑性功的增量为指标来调整微观尺度模型的长度。为了分析节点塑性屈服以后网格结构的受力,将可能发生塑性屈服的焊接球节点以及它所连接的杆件长为S的一小段采用作为微观尺度模型采用实体单元离散,结构的其余部分作为宏观尺度模型采用梁单元离散。根据经典梁理论的平截面假定推导了两种尺度模型界面上的位移增量约束方程,给出了梁-实体单元的耦合方法。同时,给出了S的估计值,并采用Zienkiewicz-Zhu后验误差估计理论来检验它是否合理。理论推导和数值算例表明,精细化有限元分析方法可以根据结构的载荷状况,来调整有限元模型,它的计算精度和全实体/壳体单元有限元模型基本相似,但可以大幅度降低计算量,具有一定的工程实用价值。最后,论文还给出网格结构精细化有限元分析方法进一步研究工作的建议。