有限元分析报告

西安市新城区某公司科研办公楼结构设计有限元分析报告撰 写 人： 王 平 班 级： 工程力学1203 学 号： 121010321 指导教师： 张卫喜 2016年6月15日目录1 工程概况22 分析依据33 荷载与计算工况43.1荷载简化及荷载组合43.2 边界条件43.3 工况54 有限元模型54.1 基本假定54.2 力学模型64.3 主要物理参数取值64.4单元选取74.5分网与有限元模型85 静力分析95.1模态结果95.2静力分析结果125.3 强度校核156基于ANSYS、PKPM、手算的误差分析176.1计算原理的不同176.2 研究对象的复杂性181 工程概况工程名称：西安市新城区某公司科研办公楼；建筑所在地：西安市；建设规模：总建筑面积约4700m2，主体结构6层，无地下室。结构总高度22.5m，底层结构高度4.5m，其余层结构高度为3.6m，几何模型图如图1所示；抗震设防烈度：抗震设防烈度为8度，设计基本地震加速度值0.2g,第一组。场地类别为类，特征周期为0.35s。周期折减系数为0.75。建筑设计使用年限：50年。结构重要性等级：二级。图1 框架几何模型图2 分析依据框架结构是由梁、板、柱以刚接相连接而成，构成承重体系的结构，即由梁、板、柱组成框架共同抵抗使用过程中出现的水平荷载和竖直荷载。本设计报告采用ANSYS有限元软件分析。根据框架结构体系特点，本结构分析主要依据以下国家规范：1国家标准：建筑结构荷载规范（GB50009-2012).北京：中国建筑工业出版社.2012；2国家标准：建筑抗震设计规范（GB50011-2010).北京：中国建筑工业出版社.2010；3国家标准：混凝土结构设计规范（GB50010-2010).北京：中国建筑工业出版社.2010；4建筑、勘察等技术文件。3 荷载与计算工况3.1荷载简化及荷载组合1）恒载：含柱、梁、板自重；2）梁间活载以换算线荷载施加（含墙、门、窗等），如图2所示；3）板面活荷载按建筑功能分类取值。图2 框架结构梁间荷载示意图3.2 边界条件1）柱底固结；2）板四周与梁固结，如图3所示。图3 框架结构板边界示意图3.3 工况工况一：模态分析；工况一：分析结构可变荷载组合下的刚度与强度；（注:由于永久荷载组合与可变荷载组合相比较小，故在此不进行校核）工况二：对结构在水平地震力和恒载组合作用下进行分析。3.3.1模态分析给模型加约束之后，进行前6阶模态分析。3.3.2 静力分析在可变荷载效应控制的荷载作用下，对框架进行刚度、强度校核。4 有限元模型4.1 基本假定1）假设混凝土为线弹性材料；2）钢筋混凝土取整体分析，计钢筋贡献；3）节点假定为刚接，底层柱与基础刚接，楼板与梁四边刚接。4.2 力学模型1）主要受力构件为柱、梁、板；2）抗侧力构件为柱；3）水平联系构件为框架梁（考虑板的贡献）；4）柱为压弯构件，框架梁含弯曲、轴向变形、扭转变形。4.3 主要物理参数取值表1 主要物理参数楼层混凝土等级弹性模量（N/m2）泊松比板厚（mm）柱（bh）（mm）横梁(bh)纵梁(bh)（mm）次梁（bh）（mm）边跨梁（bh）（mm）中跨梁（bh）（mm）1C30310100.2100600600300600300400300600250500255005006120查阅相关资料得：钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种，即分离式、组合式、整体式。分离式模型把钢筋和混凝土作为不同的单元处理，即混凝土和钢筋各自被划分为足够小的单元。考虑到钢筋是一细长材料，通常可以忽略其抗剪，一般钢筋可作为线单元来处理。混凝土则作为实体单元来处理。分离式模型中在钢筋和混凝土之间可以插入连接单元来模拟钢筋与混凝土的粘结和滑移。当钢筋和混凝土粘结较好，可以认为两者无滑移时，采用组合式模型。组合式模型中最常用两种方式第一种是分层式，即在横截面上分成许多混凝土层和若干钢筋层，并对截面的应变做出某些假定，主要适用于杆件系统的分层组合式；另一种组合方式是采用等参数单元，使用带钢筋膜的单元来进行模拟，主要适用于二维或三维结构的带膜组合式或带钢筋的组合式。整体式模型也是在考虑钢筋与混凝土之间无滑移的情况下使用的。假定钢筋弥散于整个单元中，并把单元视为连续均匀材料，由此求出的是综合了混凝土与钢筋两种材料的单元刚度矩阵。在ANSYS中可以考虑的钢筋混凝土模型一般是两种，即分离式模型和整体式模型。分离式模型虽然在理论上可行，但在进行复杂结构的分析时，对计算机的性能要求较高，计算效率较低。因此本文空间结构计算模型采用的是整体式有限单元模型来处理钢筋混凝土结构。4.4单元选取4.4.1 梁、柱单元的选取梁、柱采用空间Timoshenko单元，拟合其弯曲、拉压及扭转受力状态。Beam 188是一个二节点的三维线性梁, 适用于分析细长的梁, 其元素是基于Timoshenko 梁理论的，具有扭切变形效果。Beam 188 在每个节点上有6或7个自由度，（自由度）数目的变化是由KEYOPT(1)来控制的。当 KEYOPT(1) = 0时 (默认), 每节点有6个自由度。 分别是沿x,y,z的位移及绕其的转动。 当 KEYOPT(1) = 1时,会添加第七个自由度 (翘曲量) ，Beam 188单元示意图如图4所示。图4 Timoshenko梁单元示意图4.4.2 板单元的选取板采用弹性壳单元。Shell 181适用于薄到中等厚度的壳结构。该单元有四个节点,单元每个节点有六个自由度,分别为沿节点X,Y,Z方向的平动及绕节点X,Y,Z轴的转动。退化的三角形选项用于网格生成的过渡单元。Shell 181单元具有应力刚化及大变形功能。该单元有强大的非线性功能,并有截面数据定义,分析,可视化等功能,还能定义复合材料多层壳；Shell 181单元示意图如图5所示。图5 弹性壳单元示意图4.5分网与有限元模型4.5.1 建模思路经查阅ANSYS相关书籍得知Beam188 单元适合于分析从细长到中等粗短的梁结构，该单元基于铁木新科梁结构理论，并考虑了剪切变形的影响。有大量研究也表明该单元可以很好的模拟型钢梁与柱，计算结果与实际比较接近，那是因为钢材的材质比较均匀。但框架混凝土梁与柱材质并不均匀，混凝土内部含各种钢筋，故采用整体式模型，并把单元视为连续均匀材料，那么Beam188便可以很好的模拟框架混凝土梁与柱了。Shell 181适用于薄到中等厚度的壳结构.该单元有四个节点,单元每个节点有六个自由度,分别为沿节点X,Y,Z方向的平动及绕节点X,Y,Z轴的转动。在此处刚好可以模拟板结构。综上所述，得出以下建模思路：用Beam188、Shell181单元建立空间框架结构模型，与PKPM计算结果进行比较；4.5.2 分网对梁和柱均设置为按份数进行划分，份数取为20；对板进行划分时，考虑到板的边界与梁重合，故分网大小与梁相同，有限元模型图如图4所示。图6 有限元模型图5 静力分析 5.1模态结果对框架结构进行模态分析，提取前六阶振型图，并统计前六阶的频率，周期和振型特点。表2 ANSYS分析结构动力参数阶数周期（s)频率（Hz）振幅特点10.8121.2309平动，横向20.7561.3219平动，纵向30.7121.4046扭转40.2603.8421平动，横向50.2444.0930平动，纵向60.2294.3756扭转 表3 PKPM计算结构动力参数阶数周期（s)频率（Hz）振幅特点10.9051.1045平动，横向20.8511.1758扭转30.8421.1884平动，纵向40.2953.3944平动，横向50.2773.6088平动，扭转60.2763.6271平动，扭转对表2及表3的数据进行分析可知，对于PKPM计算结果，结构在二阶模态时发生扭转变形；对于ANSYS计算结果，结构在三阶模态时出现扭转变形。ANSYS计算各阶模态位移图如图813所示。MNMXXYZ 0.131E-03.261E-03.392E-03.522E-03.653E-03.783E-03.914E-03.001045.001175NODAL SOLUTIONSTEP=1SUB =1FREQ=1.2309USUM (AVG)RSYS=0DMX =.001175SMX =.001175图7 第一阶振形MNMXXYZ 0.928E-04.186E-03.278E-03.371E-03.464E-03.557E-03.649E-03.742E-03.835E-03NODAL SOLUTIONSTEP=1SUB =2FREQ=1.32186USUM (AVG)RSYS=0DMX =.835E-03SMX =.835E-03图8 第二阶振形MNMXXYZ 0.162E-03.324E-03.486E-03.648E-03.811E-03.973E-03.001135.001297.001459NODAL SOLUTIONSTEP=1SUB =3FREQ=1.40461USUM (AVG)RSYS=0DMX =.001459SMX =.001459图9 第三阶振形MNMXXYZ 0.124E-03.248E-03.372E-03.496E-03.620E-03.744E-03.868E-03.992E-03.001116NODAL SOLUTIONSTEP=1SUB =4FREQ=3.84214USUM (AVG)RSYS=0DMX =.001116SMX =.001116图10 第四阶振形MNMXXYZ 0.910E-04.182E-03.273E-03.364E-03.455E-03.546E-03.637E-03.728E-03.819E-03NODAL SOLUTIONSTEP=1SUB =5FREQ=4.09304USUM (AVG)RSYS=0DMX =.819E-03SMX =.819E-03图11 第五阶振形MNMXXYZ 0.156E-03.312E-03.467E-03.623E-03.779E-03.935E-03.00109.001246.001402NODAL SOLUTIONSTEP=1SUB =6FREQ=4.37558USUM (AVG)RSYS=0DMX =.001402SMX =.001402图12 第六阶振形5.2静力分析结果5.2.1 工况二MNMXXYZ 0.001286.002573.003859.005145.006432.007718.009004.010291.011577NODAL SOLUTIONSTEP=1SUB =1TIME=1USUM (AVG)RSYS=0DMX =.011577SMX =.011577图13 可变荷载组合下框架结构变形图由图13可得，在可变荷载组合作用下框架结构最大变形发生在屋面楼板，最大挠度为11.5mm。5.2.2 工况三MNMXXYZ 0.004579.009159.013738.018317.022896.027476.032055.036634.041214NODAL SOLUTIONSTEP=2SUB =1TIME=2UZ (AVG)RSYS=0DMX =.042232SMX =.041214图14 整体框架位移图MNMXXYZ 0.004579.009159.013738.018317.022896.027476.032055.036634.041214NODAL SOLUTIONSTEP=2SUB =1TIME=2UZ (AVG)RSYS=0DMX =.042232SMX =.041214图15 单榀框架位移图由上图14和图15可知，在只有地震力作用下，框架水平最大位移为36.634mm，PKPM电算结果为27.80mm。5.2.3层间位移角由ANSYS计算的各层位移及层间位移角与手算、PKPM电算的各层位移及层间位移角如表4、表5和表6所示。表4 ANSYS计算各层位移、层间位移角层数绝对位移相对位移层间位移角636.6344.9891/721531.6455.2691/683426.3765.0591/711321.3175.7791/622215.5386.3791/56419.1596.8691/655表5 手算各层位移、层间位移角层数绝对位移相对位移层间位移角627.4242.4141/1491525.0103.6621/983421.3484.6721/771316.6765.4451/661211.2315.9811/60215.2505.2501/857表6 PKPM计算各层位移、层间位移角层数绝对位移相对位移层间位移角627.802.501/1798525.303.301/1011422.004.651/744317.355.751/617211.606.451/55615.155.151/874由表4表6可知，ANSYS计算、PKPM电算以及手算结果的层间位移角均符合规范要求且相差不大。5.3 强度校核5.3.1 工况二可变荷载组合下框架结构的弯矩图如图17所示。XYZ -208747-166835-124923-83010.6-41098.5813.62142725.784637.8126550168462LINE STRESSSTEP=1SUB =1TIME=1SMIS2 SMIS2MIN =-208747ELEM=840MAX =168462ELEM=40418图16 可变荷载组合下弯矩图由图16可知，在可变荷载组合下最大弯矩发生在梁跨中，为168.462kNm。5.3.2 工况三水平地震力下框架结构的弯矩图如图18所示。XYZ -321605-243949-166293-88636.4-1098066676.3144333221989299645377302LINE STRESSSTEP=1SUB =1TIME=1SMIS2 SMIS2MIN =-321605ELEM=107871MAX =377302ELEM=2841图17 水平地震力下弯矩图由图可知，在水平地震荷载作用下，最大弯矩发生在一层柱底，为377.302kNm。6 基于ANSYS、PKPM、手算的误差分析6.1计算原理的不同6.1.1手算框架结构是多次超静定结构，纵横向的构件相互制约、相互影响。如果要精确计算结构内力，只能采用空间结构计算模型，而不能采用平面结构计算模型，所列的力法、位移法方程是很繁琐的。为了简化工作量，手算只计算一榀框架结构模型，而忽略其他方向，这样计算是建立在结构规则的前提下的。建筑抗震设计规范里提出，建筑形体规则的参考指所有平面规则和竖向规则。平面规则有扭转规则、凹凸规则；竖向规则则有侧向刚度规则和楼层承载力突变不能太大。只有满足这些条件，才可以简化成平面结构。因此，可以在计算框架结构承受水平荷载作用下的内力、位移时采用底部剪力分配法；在计算框架结构承受竖向荷载作用下的内力、位移时采用弯矩分配法。6.1.2 ANSYS与PKPM两者都是有限元软件，在有限元基本理论上都是一致的，表现如下：第一步：问题及求解域定义。根据实际问题近似确定求解域的物理性质与几何区域；第二步：求解域离散化。将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域，习惯性称为有限元网格划分；第三步：确定状态变量及控制方法。一个具体的物理问题通常可以用一组包含问题状态变量边界条件的微分方程表示，为适合有限元求解，通常将微分方程化为等价的泛函形式；第四步：单元推导。对单元构造一个适合的近似解，即推导有限单元的形式，其中包括选择合理的单元坐标系，建立单元试函数，以某种方法给出单元各状态变量的离散关系，从而形成单元矩阵；第五步：总装求解。将单元总装形成离散域的总矩阵方程（联合方程组），反映对近似求解域的要求，即单元函数的连续性要满足一定的连续条件；第六步：联立方程组求解和结果解释。有限元法最终导致联立方程组。联立方程组的求解可以用直接法、迭代法和随机法。求解结果是单元节点处状态变量的近似值。不同的是，两者在单元库、建模、后处理等子块程序有较大差异。ANSYS拥有丰富的单元库，有杆单元、梁单元、壳单元、实体单元、质量单元、接触单元、管单元等，每个单元又划分有许多单元，如solid45、solid65、solid66等。而PKPM中SATWE是采用空间杆单元模拟梁、柱，在壳单元基础上凝聚成的墙单元模拟剪力墙。综上所述，由于各种计算原理的不同，所以计算结果可能存在较大差异。6.2 研究对象的复杂性人们用ANSYS建模进行有限元分析的目的就是要模拟结构真实的内力和位移，从而校正结构是否安全。因此建模首要问题是如何选取合适的单元来模拟混凝土结构。人们之所以用ANSYS中的Beam单元分析钢结构并取得了较高的精度，是因为钢结构的材质比较均匀、连续，力学性能能用Beam单元来模拟。然而对于钢筋混凝土，人们对于其力学性能还未完全掌握。因为钢筋混凝土有钢筋、水泥、水、沙子、石子及各种掺和料或外加剂混合硬化而成的建筑材料，与一般连续均匀介质力学中的有限元方法相比，对钢筋混凝土结构进行有限元分析还存在不少困难。这些困难主要表现在：1）钢筋混凝土是由钢筋与混凝土两种力学性能很不相同的材料组成的；2）混凝土材料性质复杂。它不仅成分多样，硬化后留有空隙和自由水分，甚至还有未水化的水泥颗粒，形成好多微观裂缝。因此混凝土的应力应变关系是高度非线性的，且受其组成成型工艺及使用环境的严重影响，特别是在复杂的应力状态和加载历史下，混凝土的本构关系还有许多问题有待研究；（3）在荷载的作用下，一般钢筋混凝土结构是带裂缝工作的，而这些裂缝是随着荷载的增减和时间的推移而发生变化的；（4）混凝土的变形与时间有关，如收缩、徐变等，其规律还有待深入研究。