论文部分内容阅读
网格生成是有限元分析计算中的一个重要研究领域。常规有限元计算时,计算者通常根据自己对计算对象的理解和经验,或根据已有的计算结果进行网格划分。这种人工划分网格的方式,一方面使计算结果难以反映求解问题的真实状态,而且数据准备工作量大容易出错;另一方面由于不同计算者对计算对象的理解不同,所划分的网格也不同,造成计算结果因人而异。自适应有限元因其前处理过程高效,对网格划分的客观控制性强,可以有效地减少网格离散化误差,因而得到广泛使用。本文在系统学习自适应有限元理论的基础上,研究了新型协同转动3节点三边形壳单元的自适应网格生成方法,并使用有限元后验误差估计法和h型自适应策略对三个发生大变形和大转角的板壳结构进行了非线性分析。论文主要内容如下:(1)梳理了平面网格的自动生成技术,针对三角形网格计算简单,发展成熟的特点专门介绍了两类Delaunay三角网格的自动生成算法,分别是Bowyer-Watson逐点插入算法和Lawson换边算法。然后又介绍了两类网格质量的优化方法,分别是Laplacian光顺算法和网格拓扑优化。通过对算法的分析对比,最终采用Bowyer-Watson逐点插入算法和是Laplacian光顺算法生成初始网格。(2)介绍了新型协同转动法3节点三边形平板壳单元,并定义了随单元共同旋转的局部坐标系,在由整体节点变量计算局部节点变量时,扣除了单元刚体转动成分,从而简化了局部坐标系下的单元计算公式。将每个节点的两个中面法向量分量作为矢量型转动变量,从而在非线性增量分析过程中可以采用简单加法直接相加。由于在计算单元应变能对局部节点变量的二阶偏微分时,节点变量的微分次序具有交换性,使得在局部和整体坐标系下都得到了对称的单元切线刚度矩阵。采用MacNeal的线积分法计算假定应变,以此减轻闭锁现象对单元计算的不利影响,同时在剪切应变柔度矩阵中引入了剪切刚度的修正系数,进一步提高了单元计算精度。(3)阐述了有限元后验误差估计理论,重点介绍Zienkiewicz-Zhu(简称ZZ)后验误差估计方法的核心过程,包括误差范数计算公式和构造光滑应力,并拓展了基于超收敛点的SPR误差估计方法。最后采用绕点平均法计算修匀应力,采用应力误差L2范数计算单元与整体的相对误差。(4)提出了一套基于Delaunay三角化网格的h型自适应网格加密方案:对通过Delaunay三角化生成的初始网格进行误差估计,根据结果对误差超限单元进行“外心插点”的操作,之后对目标区域进行局部网格重生成和Laplace光滑,形成新网格。(5)将自适应有限元技术与新型协同转动3节点三边形壳单元相结合,通过对多个算例进行非线性分析,并与均匀加密后的网格进行比较,进一步验证了自适应网格加密后的结果更为精确。