论文部分内容阅读
薄壳结构在实际工程中的应用十分广泛,在外荷载作用下,薄壳结构往往会出现大变形大变位、弹塑性材料非线性以及碰撞接触、破裂破碎和穿透等严重不连续行为。传统有限元在处理这类复杂结构行为分析时需事先区分行为类型并采用专用的公式和程序来处理,对于结构行为预测和破坏全过程模拟往往难以得到满意的结果。本文基于向量式有限元这一新型分析方法,以向量力学为理论基础,通过理论分析、程序开发和数值模拟等手段实现金属薄壳结构的复杂力学行为分析。研究工作沿两条主线:一是单元发展,包括薄膜单元、薄板单元和薄壳单元;二是复杂结构行为实现,包括大变形大转动、屈曲、碰撞、断裂和穿透。第二章介绍了向量式有限元法的基本概念、假设、原理和推导思路,并对中央差分公式、共转坐标和逆向运动、材料(非)线性、静力和动力求解、等效质量和惯量矩阵以及误差分析和内力自平衡机制等问题进行了探讨;针对大变形大转动运动以及屈曲、碰撞、断裂和穿透等复杂结构行为问题及其求解思路进行了简单讨论;总结了向量式有限元法的分析步骤与流程。第三章建立了三角形CST常应变薄膜单元和四节点等参薄膜单元的向量式有限元基本公式,描述运动解析的原理及变形坐标系下单元节点内力的求解方法,同时对四节点膜单元的位置模式和内力计算的数值积分等问题提出处理方法;编制了相应的分析程序并通过算例进行验证。在复杂结构行为分析方面,针对膜材大变形大转动问题,采用本文方法跟踪获得其运动变形全过程;针对膜材的碰撞接触行为,采用“点-三角形”检测和基于中央差分式的罚接触力方法处理膜材的碰撞检测和碰撞响应问题,跟踪获得其碰撞接触全过程;针对膜材的破裂破碎行为,采用失效应力判断准则和质点分裂的向量式有限元处理,跟踪获得其破裂和破碎的全过程。第四章建立了三角形DKT薄板单元的向量式有限元基本理论,同时对质点的质量矩阵与惯量矩阵、应力计算的数值积分及插值方法、时间步长及阻尼参数的取值等问题提出处理方式;并通过平板结构的静、动力算例进行分析验证。第五章基于三角形CST薄膜单元和三角形DKT薄板单元的叠加组合,建立了三角形薄壳单元的向量式有限元基本理论,同时对质点位移和内力的合并和分离、单元应变和应力的合并和分离以及单元节点内力积分方案的处理等关键问题提出处理方法;基于C-S粘塑性材料本构模型和动态Mises屈服准则,推导获得C-S材料本构模型的弹塑性增量分析步骤,并将其引入向量式有限元薄壳单元的基本理论推导,实现金属薄壳结构同时考虑塑性硬化和应变率硬化效应的非线性分析。第六章研究了向量式有限元在薄壳结构复杂行为分析中的应用。针对薄壳结构的屈曲和后屈曲行为,采用位移和力控制处理方法来跟踪获得其变形全过程,并分析了两种控制方法的特点和应用范围;针对薄壳结构的碰撞接触行为,仍采用“点-三角形”检测和罚接触力法处理碰撞检测和碰撞响应问题,跟踪获得其碰撞接触的全过程;针对薄壳结构的破裂破碎行为,采用失效应变断裂判断准则和质点分裂的向量式有限元处理,跟踪获得其破裂和破碎的全过程;针对薄壳结构的穿透行为,结合碰撞检测、碰撞响应处理机制以及基于失效应变的断裂准则,实现刚体对薄壳结构的穿透过程。第七章综合本文基于向量式有限元的薄壳结构非线性分析基本理论,包括薄壳单元、材料本构、破裂破碎的理论公式与有效算法,实现了钢储罐结构在爆炸冲击荷载作用下的动力响应分析和破坏全过程分析,分析了内部三角形简化形式爆炸荷载作用下平顶罐、锥顶罐和拱顶罐的应力变形发展以及破坏模式。本文成果体现了向量式有限元分析方法在处理结构大变形、大变位、碰撞、断裂和穿透等强非线性和不连续力学行为中的独特优势,推动了向量式有限元的理论发展与工程应用,为薄膜结构、金属薄壳结构的复杂行为分析提供了一种新的有效方法。