考虑楼板作用的钢筋混凝土框架有限元模型及并行计算效率
- 格式:pdf
- 大小:566.28 KB
- 文档页数:7
一级注册结构工程师专业部分试卷-总则、钢筋混凝土结构(二)-1(总分100,考试时间90分钟)单项选择题(每题的四个备选项中只有一个符合题意)1. 某钢筋混凝土次梁,下部纵向钢筋配置为,fy=300N/mm2,混凝土强度等级为C25,ft=1.27N/mm2。
在施工现场检查时,发现某处采用绑扎搭接接头,其接头方式如图所示。
试问,钢筋最小搭接长度l1(mm),应与下列______项数值最为接近。
A. 793B. 992C. 1100D. 1283某抗震设计烈度为9度,抗震等级为一级的现浇钢筋混凝土多层框架结构房屋,梁柱混凝土强度等级为C30,纵筋均采用HRB400级热轧钢筋。
框架中间楼层某端节点平面及节点配筋如图所示。
2. 该节点上、下楼层的层高均为4.5m,上柱的上、下端设计值分为别=440kN?m,=380kN?m;下柱的上、下端弯矩设计值分别为=440kN?m,=500kN?m;柱上除节点外无水平荷载作用。
试问,上、下柱反弯点之间的距离Hc(m),应与下列______项数值最为接近。
A. 4.2B. 4.6C. 4.8D. 5.53. 假定框架梁KL1在考虑x方向地震作用组合时的梁端最大负弯矩设计值Mb=600kN?m;梁端上部和下部配筋均为(As=A"s=2945mm2),as=a"s=40mm;该节点上柱和下柱反弯点之间的距离为4.5m。
试问,在x方向进行节点验算时,该节点核心区的剪力设计值Vj(kN)应与下列______项数值最为接近。
A. 1000B. 1100C. 1220D. 15054. 假定框架梁柱节点核心区的剪力设计值Vj=1200kN,箍筋采用HRB335级钢筋,箍筋间距s=100mm,节点核心区箍筋的最小体积配箍率ρvmin=0.67%;as=a"s=40mm。
试问,在节点核心区,下列______项箍筋的配置较为合适。
A.8@100B.10@100C.12@100D.14@100某多层民用建筑,采用现浇钢筋混凝土框架结构,建筑平面形状为矩形,抗弯刚度较大,属规则框架,抗震等级为二级;梁、柱混凝土强度等级均为C30。
钢筋混凝土结构设计计算书一、工程概况本工程为具体工程名称,位于工程地点。
建筑总高度为高度数值米,总建筑面积为面积数值平方米。
结构形式为钢筋混凝土框架结构,设计使用年限为使用年限数值年,抗震设防烈度为烈度数值度。
二、设计依据1、相关规范和标准《混凝土结构设计规范》(GB 50010-2010)(2015 年版)《建筑抗震设计规范》(GB 50011-2010)(2016 年版)《建筑结构荷载规范》(GB 50009-2012)2、工程地质勘察报告提供了场地的地质条件、土层分布和地基承载力等参数。
3、建筑设计方案包括建筑平面布置、层高、功能分区等。
三、荷载计算1、恒载结构自重:根据构件尺寸和材料容重计算。
楼面恒载:包括楼板面层、找平层、吊顶等的重量,取值为具体数值kN/m²。
2、活载楼面活载:根据不同功能房间的使用要求,取值如下:客厅:具体数值kN/m²卧室:具体数值kN/m²厨房:具体数值kN/m²卫生间:具体数值kN/m²屋面活载:取值为具体数值kN/m²3、风荷载基本风压:根据当地气象资料,取值为具体数值kN/m²。
风荷载体型系数:根据建筑物的形状和尺寸确定。
4、地震作用根据抗震设防烈度和设计分组,计算水平地震影响系数最大值和特征周期。
四、构件尺寸初步确定1、框架柱根据轴压比限值和预估的竖向荷载,初步确定框架柱的截面尺寸。
2、框架梁考虑跨度、荷载和构造要求,初步确定框架梁的截面尺寸。
五、内力计算1、竖向荷载作用下的内力计算采用分层法或弯矩二次分配法计算框架梁柱的内力。
2、水平荷载作用下的内力计算采用 D 值法计算框架在风荷载和地震作用下的内力。
六、配筋计算1、框架柱配筋根据柱的内力组合,计算正截面受压承载力和斜截面受剪承载力,确定柱的纵筋和箍筋。
2、框架梁配筋计算梁的正截面受弯承载力和斜截面受剪承载力,确定梁的纵筋和箍筋。
3、楼板配筋按照单向板或双向板计算楼板的内力,确定板的受力钢筋和分布钢筋。
结构动力学有限元混合分层并行计算方法结构动力学是研究结构在外界载荷作用下的响应及其稳定性的一门学科。
有限元方法是结构动力学分析中广泛使用的一种数值方法。
为了提高计算效率和精度,混合分层并行计算方法应运而生。
混合分层并行计算方法是指将有限元方法与分层并行计算相结合的一种计算方法。
在结构动力学中,混合分层并行计算方法被广泛应用于解决大型结构的复杂动力学问题。
它通过将结构进行分层划分,将计算任务分配给不同的处理器进行并行计算,从而大幅提高计算速度和效率。
混合分层并行计算方法的基本思想是将结构分为多个子结构,并将每个子结构分配给一个处理器进行计算。
每个处理器独立地计算与其对应的子结构,然后通过通信机制将计算结果交换,并进行整体求解。
这种并行计算方法充分利用了计算机集群的计算能力,提高了计算效率。
在混合分层并行计算方法中,有限元方法被用于对每个子结构进行离散化,并建立相应的有限元模型。
有限元模型中的自由度数目较少,计算量相对较小,可以降低计算复杂度。
同时,分层并行计算策略使得计算任务可以被同时执行,加速了计算速度。
混合分层并行计算方法的应用范围广泛。
例如,在工程领域中,可以用于模拟大型桥梁、高层建筑等结构的动力学响应;在航空航天领域中,可以用于模拟飞机、卫星等复杂结构的动力学特性;在地震工程中,可以用于模拟地震对建筑物的影响等。
混合分层并行计算方法可以准确预测结构的振动特性、动态响应和破坏过程,为结构设计和分析提供了有力的工具。
总之,结构动力学有限元混合分层并行计算方法是一种高效、准确的计算方法。
它通过将结构进行划分和并行计算,充分利用计算机集群的计算能力,实现了大规模结构动力学分析的快速求解。
混合分层并行计算方法在工程领域中的应用潜力巨大,有着广阔的发展前景。
多层混凝土框架结构设计中的关键问题及对策摘要:混凝土框架结构中的墙体不作为承重构件承受竖向和水平荷载而只是起到填充分隔作用,所以此结构自重较轻分隔灵活且节省材料,在实际中的建设上运用广泛。
本文即结合具体工程案例详细阐述了多层混凝土框架结构设计中的关键问题及对策。
关键词:多层混凝土;框架结构;1.多层混凝土框架结构的特点框架结构是最常见的竖向承重结构,具有以下优点:(1)结构轻巧,便于布置;(2)整体性比砖混结构和内框架承重结构好;(3)可形成大的使用空间;(4)施工较方便;(5)较为经济。
框架结构是由横梁和立柱组成的杆件体系,梁-柱通过节点构成承重骨架,刚性连接,共同承受水平荷载,竖向荷载。
混凝土框架结构中的墙体不作为承重构件承受竖向和水平荷载而只是起到填充分隔作用,所以此结构自重较轻分隔灵活且节省材料,在实际中的建设上运用广泛。
但框架结构也存在一定的缺点,在力学上则会表现出抵抗水平荷载能力差、侧向刚度小、侧移大、应力集中、当框架结构房屋的高宽比较大,则水平荷载作用下的侧移也比较大,引起的颠覆作用也是比较严重,受地基的不均匀沉降也比较严重,以及受外界环境影响较大的特点。
二、工程概况在设计工业厂房建筑结构期间,某多层厂房实际占地面积有 3000 平方米,厂房建筑以钢筋混凝土为主,防火等级设计为一级。
厂房共四层,总高度为23m总面积为12000平方米,抗震等级为三级,以工业厂房实际生产需求为基础,开展建筑结构优化设计。
三、多层混凝土框架结构设计中的关键问题及对策(一)结构设计要点1、结构荷载设计为了保证工业厂房在实际应用过程中的质量,同时还要为其日常生产作业的顺利开展提供安全的环境作为支持,必须要意识到钢筋混凝土框架结构设计在其中的重要性。
在结构设计工作的具体开展过程中,设计人员要与实际情况进行结合,同时还要对各种不同类型的荷载进行综合分析,其中包括竖向、水平等荷载类型。
只有对荷载进行深入分析,才能够根据现有的荷载实现对结构的有效设计,为厂房结构的稳定性和安全性提供有效保证。
PKPM软件中钢筋混凝土结构设计的参数设置PKPM软件广泛应用于土建工程,作为设计人员不应满足于会用该软件来计算和辅助绘图,而应弄清楚重要参数的含义。
在计算模型和荷载输入正确的情况下,关键参数的错误会导致结果错误,参数的正确设置具有更重要的意义。
下面是我结合规范谈谈在实际工作中易忽略的参数如何设置,以供设计人员参考和交流。
一、合理使用软件目前,PKPM程序拥有的空间计算程序有三个,即TAT、SATWE、PMSAP。
1、TAT它是一个空间杆件程序,对柱、墙、梁都是采用杆件模型来模拟的,特殊的就是剪力墙是采用薄壁柱原理来计算的。
因此,在用TAT程序计算框剪结构、剪力墙结构等含钢筋混凝土剪力墙的结构都要对剪力墙的洞口、节点做合理的简化,有点让实际工程来适应我们的计算程序的味道。
当然,在作结构方案时,对结构作这样的调整对建筑结构方案的简洁、合理有很大的好处。
它的楼盖是作为平面内无限刚、平面外刚度为零的假设。
在新版的TAT程序中,允许增设弹性节点,这种弹性节点允许在楼层平面内有相对位移,且能承担相应的水平力。
增加了这种弹性节点来加大TAT程序的适用范围,使得TAT程序可以计算空旷、错层结构。
2、SATWE空间组合结构有限元程序,与TAT的区别在于墙和楼板的模型不同。
SATWE对剪力墙采用的是在壳元的基础上凝聚而成的墙元模型。
采用墙元模型,在我们的工程建摸中,就不需要象TAT程序那样做那么多的简化,只需要按实际情况输入即可。
对于楼盖,SATWE程序采用多种模式来模拟。
有刚性楼板和弹性楼板两种。
其中弹性楼板又分为弹性板6、弹性板3和弹性膜。
SATWE程序主要是在这两个方面与TAT程序不同。
3、PMSAP是一个结构分析通用程序。
当然,它是偏向于建筑的,但它是一个发展方向。
现在的比较著名的通用计算程序有:SAP84、SAP91、SAP2000、ANSYS、ETABS等程序,这些程序各有特长。
二、重要和易忽略的参数设置1、SATWE中的刚性板与弹性板刚性板------平面内刚度无限大,平面外刚度为零,通过梁刚度放大系数来变相的考虑楼板的平面外实际刚度。
一、荷载取值1、面荷载1)二层楼面【建筑做法:采用钢筋混凝土楼板随捣随抹平,上刷耐磨地面硬化剂。
】砼折合板厚125mm 0.125x25=3.125 kN/㎡压型钢板自重0.15 kN/㎡各专业管线设备吊挂荷载0.6 kN/㎡恒载合计 3.9 kN/㎡活荷载标准值为7.5 kN/㎡,货架区域为10 kN/㎡(按业主要求)2)不上人屋面彩色钢板岩棉夹芯板0.25 kN/㎡檩条等0.15 kN/㎡各专业管线设备吊挂荷载0.5 kN/㎡恒载合计0.9 kN/㎡屋面活荷载标准值为 0.5 kN/㎡檩条计算时,风荷载按FM要求按100年取值为0.6 kN/㎡。
4.2m标高雨篷(6m跨&10m跨),输入模型中,荷载见简图22~23轴间有1t电动葫芦,在14.225标高。
2.墙体外墙彩色钢板两层 0.25-0.5*0.12=0.19 kN/m2檩条 0.06 kN/m2玻璃纤维保温棉填充200厚 0.20x1.0 kN/m3=0.2 kN/m2∑ 0.45 kN/ m2二层外墙 0.45*7.8=3.51 kN/m 取为3.6 kN/mKFZ自重 0.1kN/m*10.45=11 kN 取为15 kN一层外墙R轴及32轴:0.45*6.1+1.5*5.0+0.5*4.2=13 kN/m B轴: 0.45*6.1+1.5*5.0+(25*0.25+0.02*20*2)*1.5=26 kN/m 内墙材料统计:CSR(甲方指定)板材二层防火墙,4轴山墙尖,报警阀间,及建筑指定的房间(common room)板材150mm厚,0.3 kN /㎡(根据厂家资料)蒸压砂加气混凝土砌块(200厚)一层防火墙,楼梯间,卫生间,及建筑指定的房间15厚水泥砂浆抹面两层 0.015*20*2=0.6 kN/㎡加气混凝土砌块200厚 0.2x5.5 kN/m3=1.1kN/㎡∑ 1.7 kN/ ㎡1)二层防火墙(CSR) 0.3*(7.8-0.6+2.65/2)=2.6 kN/m 取为10 kN/m(考虑防火卷帘及构造钢柱钢梁等)2)4轴山墙尖(CSR) 0.3*3.2=0.96 kN/m 取为2.5 kN/m3)报警阀间等(CSR) 0.3*3.0=0.9 kN/m (二层3m高)0.3*7.6=2.3 kN/m (二层7.6m高)报警阀间板材墙体自重很小,按装修考虑,未做基础及构造钢梁钢柱。
多层钢筋混凝土框架结构设计中应注意的问题随着我国现代建筑造型与功能的日趋多样化,无论是民用建筑还是工业建筑,在结构设计中总会遇到各类难题,尤其在框架结构设计中还必须注意一些看似简单,却容易忽视的问题。
本文通过分析多层钢筋混凝土框架结构的设计要求,针对计算过程中应注意的问题进行探讨,以期通过本文的阐述使我们对框架结构梁、柱、板以及结构体系中的一些注意点有清晰的认识,使设计的工程既经济又合理。
【标签】多层钢筋混凝土;框架结构;设计;工程1 多层钢筋混凝土框架结构的设计要求1.1 强柱弱梁节点设计为了实现在罕遇地震作用下,让梁端形成塑形铰,柱端处于非弹性工作状态,而没有屈服,但节点还处于弹性工作阶段。
强柱弱梁措施的强弱,也就是相对于梁端截面实际抗弯能力而言柱端截面抗弯能力增强幅度的大小,是决定由强震引起柱端截面屈服后塑性转动能否不超过其塑性转动能力,而且不致形成“层侧移机构”,从而使柱不被压溃的关键控制措施。
柱强于梁的幅度大小取决于梁端纵筋不可避免的构造超配程度的大小,以及结构在梁、柱端塑性铰逐步形成过程中的塑性内力重分布和动力特征的相应变化。
因此,当建筑许可时,尽可能将柱的截面尺寸做得大些,使柱的线刚度与梁的线刚度的比值尽可能大于1,并控制柱的轴压比满足规范要求,以增加延性。
验算截面承载力时,人为地将柱的设计弯矩按强柱弱梁原则调整放大,加强柱的配筋构造。
梁端纵向受拉钢筋的配筋不得过高,以免在罕遇地震中进入屈服阶段不能形成塑性铰或塑性铰转移到立柱上。
注意节点构造,让塑性铰向梁跨内移。
1.2 强剪弱弯剪力墙设计为了提高抗震墙的变形能力,避免发生剪切破坏,对于一道截面较长的抗震墙,应该利用洞口设置弱连梁,使墙体分为小开口墙、多肢墙或单肢墙,并使每个墙段的高宽比不小于2。
所谓弱连梁,是指在地震作用下各层连梁的总约束弯矩不大于该墙段总地震弯矩的20%;连梁不能太强,以免水平地震作用下某个墙肢出现全截面受拉,这是比较危险的。
大震弹塑性分析软件PKPM-SAUSAGE技术特点刘春明,张宏(广州建研数力建筑科技有限公司,北京100013)提要:PKPM-SAUSAGE是广州建研数力新推出的一款大震下动力弹塑性分析程序,采用了考虑塑性损伤的混凝土本构模型以及高效GPU并行显式求解技术。
本文介绍了PKPM-SAUSAGE软件的技术特点以及应用方面的一些技术细节,对工程师处理实际工程动力弹塑性分析问题具有一定的指导意义。
关键词:PKPM-SAUSAGE,弹塑性分析近年来,我国超限项目迅猛发展,每年全国都有上千个超限项目通过专家审查。
在完成超限报告的过程中,弹塑性动力时程分析结果是最费时最不好把握的内容之一。
针对结构进行动力弹塑性分析的目的是了解结构在大震下是否还具备保护人身安全的能力(大震不倒),结构的抗震体系是否合理。
这通常是根据结构在大震下的变形情况以及结构破坏的部位、程度和次序等因素来判断的。
本文通过剖析PKPM-SAUSAGE的技术特点,展示软件的技术细节,希望对结构工程师有所启发。
1 弹塑性分析软件介绍我国建筑设计规范对高层建筑结构在大震下的位移以及变形能力等性能指标提出了严格要求。
因此建筑结构工程中陆续引进各种分析软件用于结构弹塑性动力时程分析。
由于软、硬件发展的局限性,这些软件都一定程度地满足了当时人们的需要。
随着结构理论以及计算分析技术的发展,国内的弹塑性分析也经历了相应的从简化分析方法到静力弹塑性推覆再到动力弹塑性的一个发展过程。
SAP2000、Perform3D引入了Pushover分析和非线性动力分析概念,使用框架的杆单元模型进行弹塑性分析,但是由于存在建模繁琐、不能直接使用设计配筋以及计算收敛性差等缺点限制了软件的工程应用。
Midas/Building将这种塑性铰推广到使用框架铰模拟框架和纤维铰模拟剪力墙,加入了丰富的滞回模型,并可以考虑施工加载,自动引入钢筋进行结构静力推覆或动力弹塑性分析。
EPDA/PUSH引入纤维铰考虑框架、非线性分层壳单元模拟混凝土剪力墙,首先在程序中使用渐变模型模拟弹塑性发展过程。
第29卷第1期2009年2月 地 震 工 程 与 工 程 振 动JOURNAL OF E ARTHQUAKE E NGI N EER I N G AND E NGI N EER I N G V I B RATI O NVol .29No .1 Feb .2009 收稿日期:2008-03-10; 修订日期:2008-06-04. 基金项目:国家自然科学基金项目(90815014;50608024);教育部减震控制与结构安全重点实验室基金项目(0808) 作者简介:王洪涛(1973-),男,博士研究生,主要从事结构抗震工程研究。
E 2mail:hit w ht@文章编号:1000-1301(2009)01-0063-07考虑楼板作用的钢筋混凝土框架有限元模型及并行计算效率王洪涛1,谢礼立1,2(1.哈尔滨工业大学土木工程学院,黑龙江哈尔滨150090; 2.中国地震局工程力学研究所,黑龙江哈尔滨150080)摘 要:对于钢筋混凝土框架,考虑楼板与梁的协同作用可以更充分地体现结构静动力性能。
本文首先采用分层壳单元模拟楼板,通过刚臂耦合梁板,建立了可以较好模拟梁板协同工作的框架结构有限元模型,并与其它3种采用简化假定的模型进行了比较。
为了控制这种模型的分析误差,必须将梁、板划分为较小单元,所以需要耗费大量的计算机内存和时间。
为此,比较分析了采用集群并行计算与普通串行计算的求解时间,结果说明对于中等到较大的框架模型,集群计算均可以大量减少计算时间。
关键词:并行计算;集群;梁板;钢筋混凝土框架中图分类号:T U352.11 文献标志码:AF i n ite elem en t m odel of RC fram e structures and eff i c i ency of para llel com put a ti onWANG Hongtao 1,X I E L ili1,2(1.School of Civil Engineering,Harbin I nstitute of Technol ogy,Harbin 150090,China; 2.I nstitute ofEngineeringMechanics,China Earthquake Adm inistrati on,Harbin 150080,China )Abstract:For RC frame structures,models including fl oor slabs can si m ulate the static and dyna m ic behavi ors more accurately .A finite ele ment model considering cooperati on bet w een slabs and bea m s are p r oposed at first,where slabs are modelled by layered shell ele ments,and rigid links are used t o coup le the bea m s and slabs .Another three models using si m p lified assu mp ti ons or methods are als o established t o compare their perf or mances with the p r o 2posed model .I n order t o m ini m ize the err ors caused by discrete cooperati on,bea m s and slabs should be divided in 2t o s mall size ele ments,which require substantial computati on res ource .Theref ore,the computati on ti m es t o s olve the p r oposed model,using a 42node cluster and an ordinary PC,are compared and analysed,and the results sho w that the computati on ti m e can be reduced substantially f or moderate and large models .Key words:parallel computati on;cluster;fl oor slab;fra me structures引言钢筋混凝土框架一般采用现浇形式,楼板既把梁作为支座,又自然形成梁的翼缘,与梁协同工作。
楼板除了将竖向荷载传递给梁,在框架结构中还起到很多作用,比如:联系整个楼层使结构抗侧力构件协同受力,与梁形成T 型截面增强梁的抗弯能力[1~4],约束梁的扭转避免梁平面外失稳等。
在结构平面规则且板没有46地 震 工 程 与 工 程 振 动 第29卷较大开洞的情况下,结构在抗震计算时一般假定楼板在平面内为刚性,平面外刚度为零,也就是忽略了板的抗弯刚度;当结构平面较复杂时,采用弹性薄壳单元模拟楼板的弹性变形[5],但是忽略梁板的T型梁效应。
实际上由于现浇楼板参与受弯工作,提高了梁的抗弯刚度和强度,导致梁抗弯能力被低估,有可能使破坏机制由梁铰机构改变为柱铰机构[2~4],对结构抗震非常不利。
石化部规划设计院对唐山地震中48幢框架结构调查统计,凡无楼板的空旷框架,裂缝集中在梁上,且形成梁端塑性铰;凡具有整浇楼板的框架结构,破坏均发生在柱中。
可见现浇楼板参与工作后,框架柱与梁承载力比值降低,导致难以实现“强柱弱梁”,可以认为这是造成上述震害的重要因素。
另外一些计算模型将梁侧一定宽度的板作为梁的翼缘,仍采用板平面内刚性假定,这种模型考虑了板对梁抗弯能力的加强。
但是影响板有效宽度的因素很多,目前如何取值尚存争议[1]。
此外这种模型也无法模拟板的其它作用,比如梁与板在扭矩作用下相互影响[6],梁竖向变形对板内弯矩的影响[7]等等。
所以,建立一种能够准确反映梁板协同工作的有限元模型对于考察框架结构静动力特性是必要的。
钢筋混凝土楼板有限元模型一般有两种形式:离散式三维实体元和采用层模型的板壳元[8]。
其中三维实体元模拟混凝土结合桁架单元模拟钢筋,可以较真实地模拟楼板。
但是随着板尺寸的增大单元数急剧增加,不适宜做整体框架分析。
采用层模型的壳元将混凝土分为若干层,同时将钢筋网弥散为钢筋层,由此可以模拟板在厚度方向的应力变化。
这种模型相对简单,并且具有较好的精度,单元数目相对实体元要少得多,适合于做结构整体分析。
钢筋混凝土梁柱一般采用弹塑性梁单元模型,根据刚度沿杆长的变化规律,可以分为集中刚度模型、分布刚度模型和平均刚度模型。
集中刚度模型是将塑性变形集中于单元端部的一点建立刚度矩阵,适用于带塑性铰的梁模型,主要有单分量模型,双分量模型,多弹簧模型等。
分布刚度模型更细致地考虑刚度沿杆长的变化规律,如分段变刚度杆单元模型。
平均刚度模型假设刚度沿长度方向不变,取平均值计算刚度矩阵,适用于细分单元的模型。
通过刚臂连接板与梁,是一种简单可行的模拟梁板协同工作的方法。
文献[9]把梁截面形心和板截面形心连线看作刚臂,实现两形心位移分量的变形协调。
但是,梁板在实际工作状态中是连续耦合在一起的,而刚臂单元只在单元节点处联系梁板,所以如果梁板单元尺寸过大,会造成较大的误差[4]。
以前的研究工作也表明,采用刚臂连接梁板模拟梁板协同工作,必须将梁板划分为较小单元才能减少误差。
由于板单元数量与长度方向划分份数的平方成正比,所以,当结构尺寸增大时,模型的单元数量急剧增加,对计算机内存及CP U时间的需求也急剧增加。
文献[4]采用子结构法和超级单元法对这种有限元模型进行分析,获得了很高的计算效率。
本文首先根据不同假定建立了四种有限元模型,比较了各模型特点及动力性能差别;进而为解决梁板协同工作框架模型的大规模有限元计算问题,尝试采用集群并行计算手段,与普通串行计算求解时间进行了比较分析,并给出了求解此类模型并行效率的规律。
1 梁板有限元模型介绍为研究梁板协同工作的钢筋混凝土框架有限元建模,本文采用4种方法,针对图1所示钢筋混凝土框架结构建立了有限元模型,并比较了各种模型的特点和动力特性差别。
图1 五层框架厂房简图Fig11 Sketch of a five-st ory fra me 图2为应用大型有限元程序ANSYS 建立梁板模型的4种有限元模型简图,包括上排显示梁元和壳元尺寸的截面图和下排的几何模型图。
楼板均采用SHELL91壳单元分层模拟钢筋混凝土和钢筋网,梁柱采用BEAM189梁单元,根据实际配筋情况定义纤维模型截面。
其中模型A 采用刚臂将梁轴线处节点与板节点连接,模拟梁板协同工作。
模型B 利用ANSYS 梁单元的轴线偏移功能将梁的轴线上移至距板顶一半板厚处,梁单元与板壳元共用节点,也实现了对T 型梁效应的模拟。
但是,梁轴线偏移后,框架系统分配给梁的轴力将使梁产生附加弯矩,造成一定程度的误差。
模型C 不采用刚臂,也不偏移梁轴线,板壳元与梁共用节点,这种模型考虑了板本身的抗弯贡献,但是忽略了T 型梁效应。
模型D 忽略了板的抗弯贡献,采用平面内刚性假定,是普通结构分析经常采用的模型。
根据四种模型的特点,可以认为模型A 与实际工作状态最为接近,只要网格尺寸划分足够细,也可以将误差控制在一定范围内。
本文分别采用四种模型建立图1所示五层钢筋混凝土框架有限元模型,板厚分别为120mm 和200mm 。
为了突出本文对四种建模方法的比较,四种模型除以上所述建模方法上的区别外,单元尺寸、材料特性等其他参数均相同。
图3显示采用模型A 所建立五层框架有限元模型,所有梁板结合处均采用刚臂连接,并使板顶与梁顶齐平。
对于模型A,为达到一定的精度,必须将梁板划分为较小的单元。
表1列出框架结构中6m 跨框架梁与板划分不同份数时重力作用下跨中和梁端绕强轴弯矩值,并给出相对于划分为20份时的误差百分比。
可见,要想达到较高精度必须划分很细的有限元网格。
当网格尺寸划分更小时,应力集中会导致收敛越来越困难。
图2 梁板有限元模型Fig 12 Finite ele ment model of bea m s and fl oorslabs图3 五层框架结构有限元模型Fig 13 Finite element model of a five -st ory frame表1 梁板单元尺寸引起误差Table 1 Err ors due t o size of bea m and slab elements 梁板划分份数跨中弯矩/(kN ・m )梁端弯矩/(kN ・m )跨中误差%梁端误差%210.3530.038-2.30254-99.6385410.686 4.4910.83986-57.2774610.6 6.20.02831-41.0198810.5867.293-0.1038-30.62211010.598.106-0.06606-22.88811210.5938.751-0.03775-16.75231410.5949.286-0.02831-11.66291610.5959.744-0.01887-7.305941810.59610.151-0.00944-3.434172010.59710.5122 四种框架模型动力特性比较对框架结构四种有限元模型分别进行了模态分析,表2列出板厚分别为120mm 和200mm 四种模型的前五阶自振频率。