基于MATLAB的高层框架剪力墙结构剪力墙位置变化对模态及时程分析结果的影响.docVIP

基于MATLAB的高层框架剪力墙结构剪力墙位置变化对模态及时程分析结果的影响 1 总体过程介绍 由于本文的计算结果是用MATLAB程序计算的，而对框架剪力墙结构无论是进行模态分析还是进行地震作用下的模态分析，首先要知道结构的刚度矩阵和质量矩阵。为了得到结构的刚度矩阵，采用有限元分析软件ANSYS来完成。 ANSYS软件得到结构的刚度矩阵采用了刚度法，为求某一层的层间刚度，锁住对应的上下层，加载单位力后，利用ansys强大的有限元分析能力，求出此层对应的最大位移，之后通过计算可以得到层间刚度，进而得到整个结构的刚度矩阵。而质量矩阵的求法，采用简单的手算来完成，考虑在ANSYS建模中所有的构件，乘以混凝土的容重后，得到结构每层的质量，之后得到结构的质量矩阵。 得到刚度矩阵和质量矩阵之后，在MATLAB中写出模态分析的程序，导入求得的刚度和质量，求得结构的模态分析结构。之后编写出时程分析的程序，导入结构的刚度、质量、及地震波情况，求得结构的时程分析结果。 2 ANSYS前处理过程 2.1 计算参数确定 整个框架剪力墙结构有限元模型采用的材料为混凝土,且在整个模型中选用的混凝土种类相同且其强度等级为C30,弹性模量为3.0xl010pa,泊松比是PRXY=0.2,密度DENS=2500kg/m3。此外,该模型的建立只考虑混凝土忽略钢筋各参数,在结构中,梁、板和柱节点处均为刚接,使建模简化。其他荷载工况如下: 抗震设防类别:乙类; 抗震等级:二级框架,一级剪力墙; 抗震设防烈度:8度,0.20g; 场地类别:m类 设计地震分组:第一组 根据《建筑抗震设计规范》(GB50011-2010)查表可知,水平地震影响系数最大值 =0.16,特征周期Tg=0.45s。 底层层高6m，2-20层层高为3m，建筑结构总长度36m，总宽度20m，总高度63m， 立柱—1-20层立柱截面采用500mm×500mm。 梁—1-20层梁截面采用300mm×600mm。 楼板厚度取100mm。 剪力墙厚度取300mm。 利用ANSYS有限元分析软件,在建立框架剪力墙结构的有限元模型过程中,整个模型采用两种单元:BEAM4和SHELL63。框架梁和框架柱采用的单元是BEAM188,而剪力墙和楼板采用的单元是SHELL63。 该框架剪力墙结构平面图如图1所示，剪力墙位置变化后的结构平面图如图2所示。图中剪力墙用加深黑色粗线标出。 图1：建筑结构平面图（剪力墙位置不变） 图2：建筑结构平面图（剪力墙位置变化） 2.2 边界条件 在高层建筑结构中,由于结构下部的基础埋置的深度较大,所以可以忽略地面处结构的侧移和转动,在建模过程中是用固支将结构模型的底部约束住,将其作为基础。从结构受力模型上来说,整个结构为一承受侧向载荷的悬臂梁,其底部嵌固在基础上。 2.3 ANSYS建模情况 高层建筑结构体系庞杂,形式多种多样,如若直接对实体建模,那么其工作量也将非常繁杂。因此,本论文为使建模工作简化,采用了直接生成法,APDL语言中的循环语句定义，利用命令流输入的方式,建立了20层框架-剪力墙结构有限元模型。框架剪力墙结构有限元模型如图3所示。 图3：框架剪力墙结构有限元模型 3 ANSYS后处理过程 3.1 ANSYS导出刚度具体过程 采用刚度法，导出每层的层间刚度，以求i层的层间刚度为例，介绍刚度法的步骤。步骤如下： （1）约束住i层上下两层。 （2）在i层上施加单位力 （3）求出i层的层间位移 （4）分析层间位移结果，得到i层的层间刚度 3.2 ANSYS刚度法的难点 在求i层层间刚度时，对i层楼板面的28个点加载力，力的大小为1/28N.加载力总大小为1N，求出对应位移Si 。难点有两点： （1）针对28个点加载力是否合理。 （2）对顶层加载和对中间层加载时梁板作用不同。 针对问题（1），在梁板刚度EX取值大小和立柱一样的情况下，对每层的28个节点和其中2个节点加载力，力的总和相同，但在ANSYS中计算结果相差较大（相差10倍）。由此可知，针对28个点加载力的合理性值得推敲。因为如果针对其中的2点加力，由于梁板刚度有限，会有局部的应力集中，导致最大位移变大，使得计算结果不准确。解决方案是扩大梁板的弹性模量，在ANSYS建模中体现为对材料属性的定义上。在对梁板弹性模量进行扩大的几次尝试中，依次将梁板的弹性模量扩大为原来的10倍、10