基于LSDYNA刚性钢框架受撞击整体性能有限元分析.docVIP

基于LSDYNA刚性钢框架受撞击整体性能有限元分析 　　摘要：基于LS-DYNA软件,以钢梁吊装过程中坠落撞击钢框架的横梁为例,演示了数值仿真计算的过程。重点介绍了结构模型化,获得并分析了钢框架所受碰撞能量转换、以及框架受冲击响应的一般规律和特点。从而为钢框架设计、维护和损伤评估等提供理论上的支持。 　　关键词：撞击；钢框架；响应；动力分析 　　 　　1引言 　　 　　随着经济的发展、钢产量的日益提高，钢以其优异的材料性能越来越多地被应用到建筑结构中。然而，钢结构的现场施工过程大多采用吊装到位、组装。吊装或组装的施工过程中不可避免地要出现吊装或组装构件坠落对底部已装好结构产生撞击影响，研究钢框架在其受撞击中的响应及损伤有重要的意义。经过多年研究，对单个杆件的撞击理论研究较多，如文献[1]的脉冲响应函数和模态叠加法，文献[2]的应力波解析式方法，文献[3]以微分法为基础得到的撞击问题解析解。然而这些理论广泛应用到实际工程中有很大的难度，以至于现行各国规范中分析结构撞击问题所采用的方法都是基于有限试验条件下总结的经验公式[4]，而相对于价格高昂的试验，有限元模拟的优势越来越明显。本文以一榀钢框架受一根钢梁的坠落撞击为例，基于LS-DYNA软件进行有限元分析。 　　1.1 显式有限元动力分析的理论基础 　　钢框架受冲击是一个瞬态的复杂物理过程，它包含以大位移、大应变为特征的几何非线性和以材料弹塑性变形为典型特征的材料非线性，这些非线性物理现象的综合作用结果使框架碰撞过程的精确描述和求解十分困难。所以钢框架的撞击问题一般都采用显式有限元算法。 　　1.2 LS-DYNA软件简介 　　LS-DYNA是功能齐全的非线性分析程序包，可以模拟各种复杂的非线性动态过程，求解各种几何非线性、材料非线性和接触非线性问题，其显式算法特别适合于分析各种非线性冲击动力学问题，如冲击、爆炸、结构撞击等[5]。 　　 　　2钢框架受冲击有限元模型 　　 　　2.1接触算法 　　接触构件之间的相互作用通过接触算法来完成，本文采用对称罚函数法,其基本的分析原理为,每一时步先检查各从节点是否穿透主表面,没有穿透则对该从节点不作任何处理。如果穿透,则在该从节点与被穿透主表面之间引入一个较大的接触力,其大小与穿透深度、主表面刚度成正比,称为罚函数值。本文采用表面-表面接触类型(ASTS)。 　　2.2材料模型 　　本文采用塑性随动模型,此模型与应变率相关，可考虑失效。通过在0（仅随动硬化）和1（仅各向同性硬化）间调整硬化参数 来选择各向同性硬化或随动硬化。该模型采用Cowper-Symonds本构方程[6]，用与应变率有关的因数表示屈服应力。钢材密度?籽取为7800kg/m3；泊松比?自取为0.3；失效是由是失效应变来控制，失效应变取为0.34。 　　2.3框架结构有限元模型 　　本框架柱净高3米，梁净跨3米。梁柱均采用焊接H型钢，柱子截面为H300×300×10×16，梁截面为H200×300×10×16，节点采用刚性节点，柱强轴方向与梁连接，具体做法为梁翼缘与柱翼缘采用对接焊缝连接，梁腹板与柱翼缘处的连接板采用高强度螺栓摩擦型单剪连接。柱子与基础刚接。 　　有限元模型采用八节点六面体单元，该单元在默认情况下，采用单点积分方法，但是在本模型中由于是单点受力的撞击问题，计算结果会出现很严重的沙漏问题，为了防止沙漏问题影响计算结果的精度，本文采用全积分单元算法。 　　由于LS-DYNA软件的单元类型有限，所以在有限元建模时做了一定的简化：省略螺栓孔，由于梁腹板和柱子的焊接板之间用高强螺栓连接，所以把柱翼缘处连接板与梁腹板之间的连接视为完全连接，耦合所有自由度。整个框架以及框架节点有限元模型如图1所示。坠落梁几何尺寸与框架梁几何尺寸相同，假定梁从3m高度处坠落。 　　 　　3计算结果 　　 　　3.1碰撞能量分析 　　撞击过程中动能向应变能转换，应变能有坠落梁和整个框架承担。其中柱子和梁分担大部分的应变能，节点分担的应变能较少，可见刚性框架梁在受撞击作用下节点并不是主要耗能装置。 　　3.2框架梁受撞击处及框架节点塑性应变分析 　　节点的塑性区域随着时间一直在发展，并没有因为撞击的结束而停止发展。而且塑性应变的最大值一直位于框架节点的下面焊缝处，节点的上下焊缝附近是塑性应变最严重的区域。所以如果在撞击发生之前节点已经施焊，可考虑重新焊接或更换新梁重新施焊。 　　在与坠落梁两个翼缘接触的区域，框架梁腹板塑性发展比较严重，从接触开始 内，塑性区发展比较快，而 后塑性应变值略有增加而塑性区域发展较小。对于在以后主要承担弯矩的梁构件，梁中是弯矩最大的部分，所以对于框架梁上翼缘和梁腹板上部发展的塑性区域，应该在施工中引起重