基于向量式有限元的大跨度钢结构施工力学分析方法.docVIP

基于向量式有限元的大跨度钢结构施工力学分析方法 　　摘 要：针对大跨度钢结构施工容易产生刚体位移等强烈非线性过程的特点，采用向量式有限元分析方法进行施工过程模拟.基于向量式有限元的基本理论，引入张拉索单元，通过控制张拉索单元的原长，实现大跨度预应力钢结构张拉全过程模拟；引入千斤顶单元，实现大跨度钢结构的脱架模拟.在此基础上，利用MATLAB语言编制了基于向量式有限元理论的施工过程分析程序，并通过悬臂梁拼装、索桁架张拉成形以及大跨度张弦桁架的张拉施工过程模拟分析，与传统分析方法进行比较，验证理论推导和自编程序的正确性与有效性.结果表明，采用基于向量式有限元的施工力学分析方法能够准确模拟结构拼装、预应力张拉成形和结构脱架等施工过程. 中国论文网 /4/view-7257380.htm 　　关键词：大跨度钢结构；施工过程；向量式有限元；张拉索单元；千斤顶单元 　　中图分类号：TU311.4；TU393.3 文献标识码：A 　　文章编号：1674-2974（2016）03-0048-07 　　早期的施工力学问题主要存在于桥梁[1-3]和高层建筑[4]中，随着大跨空间结构、复杂结构的蓬勃发展，结构施工的周期和复杂性都大大增加，而且施工过程与结构最终成型状态关系更加密切，施工力学问题在大跨度钢结构中受到了充分的重视，但国外在大跨度钢结构施工力学问题方面公开发表的文献较少[5-6].国内对施工力学的研究则主要基于时变力学理论[7]，将施工过程离散为若干施工阶段进行分析，常采用生死单元法和分步建模法[8]，将连续的施工过程进行离散化求解.生死单元技术采用一次性建模，然后按照实际施工步骤逐步“杀死”或“激活”单元来模拟整个施工过程结构的受力及变形状态，避免了单元网格的重新划分，只需建立一次整体模型，但其缺点是单元被激活后可能发生漂移而与实际的安装位形不符，出现较大偏差甚至求解不能收敛；分步建模法是按照施工步骤边建模边求解，可精确控制施工过程中构件的安装位形，不存在生死单元技术由于“死”单元的“漂移”而导致刚度矩阵病态的问题，其缺点是每个施工步骤都需导入上个施工步分析的应力状态作为本次分析的初始应力状态，重复建模.而且传统有限元方法在大变形、大变位等这类施工过程中经常涉及到的非线性问题求解方面往往存在较大困难. 　　向量式有限元[9-11]是一种基于动力学求解的数值方法，它从传统的牛顿力学出发，建立起一套完整的理论.此方法可以应用于所有符合牛顿定律的力学问题求解，不需求解联立方程组，不存在非线性求解的收敛问题，尤其适合于动力问题.国内已有部分学者将其引入到结构分析中[12-16]，可以完成诸如大变形、大位移，甚至是刚体位移等一系列非线性分析.本文利用向量式有限元理论计算与时间的依存性，进行大跨度钢结构施工力学分析，为大跨度钢结构的施工力学分析提供了一种新的手段. 　　1 向量式有限元概述 　　向量式有限元的理论构架不同于经典结构力学，选择了一组不同的概念描述和简化假设.在向量式有限元基本理论中推翻了传统结构力学中的一些简化假设，例如刚性杆件、运动和变形的分解以及路径独立的过程和静态解，杆件的变位量和变形量是没有限制的，而把时间也作为分析的一个变量来考虑.因此，向量式有限元能够考虑运动进行的全部过程，处理作用力和操作环境持续变化的真实状况.同时，向量式有限元引入了数值计算方法，避免了多层次的迭代计算，求解过程中不形成刚度矩阵，因此不仅能够方便地处理大变形、大变位等几何非线性问题，也能够处理材料非线性和状态非线性等不连续行为. 　　1.1 求解过程 　　根据牛顿第二定律，对于每个质点有： 　　2.2 千斤顶单元 　　大跨度钢结构在安装过程中采用支撑胎架，为便于卸载，一般使用千斤顶作为临时支撑与结构之间的连接，千斤顶在卸载施工中有较大的承载能力，且便于控制.基于千斤顶工作中受压而不受拉的特点，可采用与张拉索单元类似的模拟方法，建立千斤顶单元的内力计算公式.不同的是，千斤顶单元只能受压不能受拉，因此，当fA2B20时，E0=0. 　　3 大跨度钢结构施工力学分析 　　施工力学分析方法主要包括有限单元法、时变单元法和拓扑变化法等.时变单元法是指离散网格不变，通过单元大小的变化来实现求解区域的变化，但存在数值积分稳定性问题.拓扑变化法应用拓扑学原理用数值手段实现求解区域的变化，但要求时变次数不能太多，否则计算效率不高.有限单元法因为理论成熟，易于程序化，得到了广泛的应用.但对于大变形、大变位甚至刚体位移等非线性过程的求解往往很难收敛.本文采用向量式有限元方法，可根据实际施工顺序通过确定新增单元或节点，直接建立新增构件加入初始模型进行分析.由于向量式有限元求解本身即为动态求解过程，因此不需调整参数，真实模拟实际施工顺序，跟踪受力和变形