有限元分析报告

1、西安市新城区某公司科研办公楼结构设计有限元分析报告撰写人： 王平班 级：工程力学1203学 号：_J指导教师：张卫喜2016年6月15日错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签 错误!未定义书签目录1工程概况2分析依据3 荷载与计算工况 荷载简化及荷载组合边界条件工况4有限元模型基本假定力学模型主要物理参数取值单元选

2、取分网与有限元模型5静力分析模态结果静力分析结果强度校核6基于ANSYS PKPM手算的误差分析.计算原理的不同研究对象的复杂性1工程概况工程名称：西安市新城区某公司科研办公楼；建筑所在地：西安市；2建设规模：总建筑面积约4700m,主体结构6层，无地下室。结构总高度，底层结构高度，其余层结构高度为，几何模型图如图1所示；抗震设防烈度：抗震设防烈度为 8度，设计基本地震加速度值，第一组。场 地类别为U类，特征周期为。周期折减系数为。建筑设计使用年限：50年。结构重要性等级：二级。图1框架几何模型图2分析依据框架结构是由梁、板、柱以刚接相连接而成，构成承重体系的结构，即由梁、 板、柱组成框架共同

3、抵抗使用过程中出现的水平荷载和竖直荷载。 本设计报告采 用ANSYSt限元软件分析。根据框架结构体系特点，本结构分析主要依据以下国家规范：1 国家标准：建筑结构荷载规范(GB50009-2012).北京：中国建筑工业 出版社.2012 ;2 国家标准：建筑抗震设计规范(GB50011-2010).北京：中国建筑工业 出版社.2010 ;3 国家标准：混凝土结构设计规范(GB50010-2010).北京：中国建筑工 业出版社.2010 ;4建筑、勘察等技术文件3荷载与计算工况荷载简化及荷载组合1)恒载：含柱、梁、板自重;2)梁间活载以换算线荷载施加(含墙、门、窗等)，如图2所示;3)板面活荷载按

4、建筑功能分类取值AMCVQELEMENTSACELPRES-N!| HrROTNFORNMOMJUN 14 201618:35:44勺 alfnlz、x-9800-5933.33-7866.67-4000-2066.67-133.33318005666.673733.337600图2框架结构梁间荷载示意图边界条件1）柱底固结；2）板四周与梁固结，如图3所示。图3框架结构板边界示意图工况工况一：模态分析；工况一：分析结构可变荷载组合12乐1馄k下的刚度与强度；（注：由于永久荷载组合与可变荷载组合相比较小，故在此不进行校核）工况二：对结构在水平地震力和恒载组合作用下进行分析。模态分析给模型加约束之

5、后，进行前6阶模态分析静力分析在可变荷载效应控制的荷载作用下，对框架进行刚度、强度校核4有限元模型基本假定1）假设混凝土为线弹性材料；2）钢筋混凝土取整体分析，计钢筋贡献；3）节点假定为刚接，底层柱与基础刚接，楼板与梁四边刚接力学模型1）主要受力构件为柱、梁、板；2）抗侧力构件为柱；3）水平联系构件为框架梁（考虑板的贡献）；4）柱为压弯构件，框架梁含弯曲、轴向变形、扭转变形主要物理参数取值b x io:o表1主要物理参数混凝土弹性模量泊松板厚(bX h边跨梁中跨梁纵梁次梁楼层(b X h)(bX h )等级(N/ )比(mm )(bX h)(b X h)(mm(mm(mm(mm(mm柱横梁（b

6、xh ）1C30100600 X 600 300 X 600 300 X 400 300 X 600 250 X 50025500X5006 120查阅相关资料得：钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种，即分离式、组合式、整体式。分离式模型把钢筋和混凝土作为不同的单元处理， 即混凝土和钢筋各自被划 分为足够小的单元。考虑到钢筋是一细长材料，通常可以忽略其抗剪，一般钢筋 可作为线单元来处理。混凝土则作为实体单元来处理。分离式模型中在钢筋和混 凝土之间可以插入连接单元来模拟钢筋与混凝土的粘结和滑移。当钢筋和混凝土粘结较好，可以认为两者无滑移时，采用组合式模型。组合 式模型中最常用两种方

7、式第一种是分层式，即在横截面上分成许多混凝土层和若 干钢筋层，并对截面的应变做出某些假定，主要适用于杆件系统的分层组合式； 另一种组合方式是采用等参数单元， 使用带钢筋膜的单元来进行模拟，主要适用 于二维或三维结构的带膜组合式或带钢筋的组合式。整体式模型也是在考虑钢筋与混凝土之间无滑移的情况下使用的。假定钢筋弥散于整个单元中，并把单元视为连续均匀材料，由此求出的是综合了混凝土与 钢筋两种材料的单元刚度矩阵。在ANSY防可以考虑的钢筋混凝土模型一般是两种，即分离式模型和整体式 模型。分离式模型虽然在理论上可行，但在进行复杂结构的分析时，对计算机的 性能要求较高，计算效率较低。因此本文空间结构计算

8、模型采用的是整体式有限 单元模型来处理钢筋混凝土结构。单元选取梁、柱单元的选取梁、柱采用空间Timoshenko单元，拟合其弯曲、拉压及扭转受力状态。Beam 188是一个二节点的三维线性梁,适用于分析细长的梁,其元素是基 于Timoshenko梁理论的，具有扭切变形效果。Beam188在每个节点上有6或7个自由度，（自由度）数目的变化是由 KEY0PT（1来控制的。当KEY0PT（1） = 0 时（默认）,每节点有6个自由度。分别是沿x,y,z的位移及绕其的转动。 当KEY0PT（1） = 1时,会添加第七个自由度（翘曲量），Beam 188单元示意图如图 4所示。图4 Timoshenko

9、梁单元示意图板单元的选取板采用弹性壳单元。Shell 181适用于薄到中等厚度的壳结构。该单元有四个节点 ,单元每个节 点有六个自由度,分别为沿节点X,Y,Z方向的平动及绕节点X,Y,Z轴的转动。退 化的三角形选项用于网格生成的过渡单元。 Shell 181单元具有应力刚化及大变 形功能。该单元有强大的非线性功能，并有截面数据定义，分析,可视化等功能, 还能定义复合材料多层壳；Shell 181单元示意图如图5所示。图5弹性壳单元示意图分网与有限元模型建模思路经查阅ANSYS相关书籍得知Beam188单元适合于分析从细长到中等粗短的 梁结构，该单元基于铁木新科梁结构理论，并考虑了剪切变形的影响

10、。有大量研究也表明该单元可以很好的模拟型钢梁与柱，计算结果与实际比较接近，那是因为钢材的材质比较均匀。但框架混凝土梁与柱材质并不均匀，混凝土内部含各种 钢筋，故采用整体式模型，并把单元视为连续均匀材料，那么Beam188便可以很 好的模拟框架混凝土梁与柱了。Shell 181适用于薄到中等厚度的壳结构该单元有四个节点,单元每个节点 有六个自由度,分别为沿节点X,Y,Z方向的平动及绕节点X,Y,Z轴的转动。在此 处刚好可以模拟板结构。综上所述，得出以下建模思路：用Beam188 Shell181单元建立空间框架结构模型，与 PKPM计算结果进行比较；分网对梁和柱均设置为按份数进行划分，份数取为2

11、0；对板进行划分时，考虑到板的边界与梁重合，故分网大小与梁相同，有限元模型图如图4所示。图6有限元模型图5静力分析模态结果对框架结构进行模态分析，提取前六阶振型图，并统计前六阶的频率，周期和振型特点。表2 ANSYS分析结构动力参数阶数周期（s）频率（Hz）振幅特点1平动，横向2平动，纵向3扭转4平动，横向5平动，纵向6扭转表3 PKPM计算结构动力参数阶数周期（s）频率（Hz）振幅特点1平动，横向扭转2平动，纵向平动，横向平动，扭转对表2及表3的数据进行分析可知，对于 PKPM计算结果，结构在二阶模态时发生扭转变形；对于ANSYS计算结果，结构在三阶模态时出现扭转变形。ANSYS计算各阶模态

12、位移图如图813所示图7第一阶振形图8第二阶振形图9第三阶振形图10第四阶振形图11第五阶振形图12第六阶振形静力分析结果工况二图13可变荷载组合下框架结构变形图由图13可得，在可变荷载组合作用下框架结构最大变形发生在屋面楼板,最大挠度为工况三图14整体框架位移图图15单榀框架位移图由上图14和图15可知，在只有地震力作用下，框架水平最大位移为，PKPM 电算结果为。层间位移角由ANSYS+算的各层位移及层间位移角与手算、PKPM电算的各层位移及层 间位移角如表4、表5和表6所示。表4 ANSYS计算各层位移、层间位移角层数绝对位移相对位移层间位移角61/72151/68341/71131/6

13、2221/56411/655表5手算各层位移、层间位移角层数绝对位移相对位移层间位移角61/149151/98341/77131/66121/60211/857表6 PKPM计算各层位移、层间位移角层数绝对位移相对位移层间位移角61/179851/101141/74431/61721/55611/874由表4表6可知，ANSYS+算、PKPM电算以及手算结果的层间位移角均符 合规范要求且相差不大。强度校核工况二可变荷载组合下框架结构的弯矩图如图 17所示图16可变荷载组合下弯矩图由图16可知，在可变荷载组合下最大弯矩发生在梁跨中，为工况三水平地震力下框架结构的弯矩图如图18所示、:图17水平地

14、震力下弯矩图由图可知，在水平地震荷载作用下，最大弯矩发生在一层柱底，为6基于ANSYS PKPM手算的误差分析计算原理的不同手算框架结构是多次超静定结构，纵横向的构件相互制约、相互影响。如果要精 确计算结构内力，只能采用空间结构计算模型，而不能采用平面结构计算模型， 所列的力法、位移法方程是很繁琐的。为了简化工作量，手算只计算一榀框架结 构模型，而忽略其他方向，这样计算是建立在结构规则的前提下的。建筑抗震设计规范里提出，建筑形体规则的参考指所有平面规则和竖向规则。 平面规则 有扭转规则、凹凸规则；竖向规则则有侧向刚度规则和楼层承载力突变不能太大。 只有满足这些条件，才可以简化成平面结构。因此，

15、可以在计算框架结构承受水 平荷载作用下的内力、位移时采用底部剪力分配法；在计算框架结构承受竖向荷载作用下的内力、位移时采用弯矩分配法。ANSYS与 PKPM两者都是有限元软件，在有限元基本理论上都是一致的，表现如下：第一步：问题及求解域定义。根据实际问题近似确定求解域的物理性质与几何区域；第二步：求解域离散化。将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域，习惯性称为有限元网格划分；第三步：确定状态变量及控制方法。一个具体的物理问题通常可以用一组包 含问题状态变量边界条件的微分方程表示， 为适合有限元求解，通常将微分方程 化为等价的泛函形式；第四步：单元推导。对单元构造

16、一个适合的近似解，即推导有限单元的形式， 其中包括选择合理的单元坐标系，建立单元试函数，以某种方法给出单元各状态 变量的离散关系，从而形成单元矩阵；第五步：总装求解。将单元总装形成离散域的总矩阵方程（联合方程组），反映对近似求解域的要求，即单元函数的连续性要满足一定的连续条件；第六步：联立方程组求解和结果解释。有限元法最终导致联立方程组。 联立 方程组的求解可以用直接法、迭代法和随机法。求解结果是单元节点处状态变量 的近似值。不同的是，两者在单兀库、建模、后处理等子块程序有较大差异。ANSYS拥有丰富的单元库，有杆单元、梁单元、壳单元、实体单元、质量单元、接触单元、 管单元等，每个单元又划分有

17、许多单元，如 solid45、solid65、solid66等。而 PKPM中 SATWE是采用空间杆单元模拟梁、柱，在壳单元基础上凝聚成的墙单元 模拟剪力墙。综上所述，由于各种计算原理的不同，所以计算结果可能存在较大差异。研究对象的复杂性人们用ANSYS建模进行有限元分析的目的就是要模拟结构真实的内力和位 移，从而校正结构是否安全。因此建模首要问题是如何选取合适的单元来模拟混 凝土结构。人们之所以用ANSY防的Beam单元分析钢结构并取得了较高的精度， 是因为钢结构的材质比较均匀、连续，力学性能能用Beam单元来模拟。然而对于钢筋混凝土，人们对于其力学性能还未完全掌握。因为钢筋混凝土有钢筋、水 泥、水、沙子、石子及各种掺和料或外加剂混合硬化而成的建筑材料，与一般连 续均匀介质力学中的有限元方法相比，对钢筋混凝土结构进行有限元分析还存在