基于ANSYS的钢筋混凝土结构实例分析
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
解放军理工大学工程兵工程学院课程论文
基于ANSYS的钢筋混凝土结构实例分析课程名称:高等军桥结构分析
专业:桥梁与隧道工程
学生姓名:马森
学号: S201103062
指导教师:郭志昆教授
陈万祥讲师
时间: 2012年7月14日
基于ANSYS的钢筋混凝土结构实例分析
马 森
(解放军理工大学工程兵工程学院,江苏,南京,210007)
摘要:讨论了在大型有限元软件ANSYS中用于混凝土材料的Solid65单元, 介绍了混凝土和钢筋共同工作时的建模方法及相互连接的处理, 利用Solid65单元对钢筋混凝土板进行了分析。分析表明, 用Solid65单元模拟钢筋混凝土材料所得到的结果, 能较好地反映钢筋混凝土板的非线性力学特征。
关键词:Solid65单元;ANSYS;钢筋混凝土结构
中图分类号:TU37511; O242121文献标识码: A
钢筋混凝土是土木工程中应用最广泛的材料,钢筋混凝土相关力学问题的分析是进行土木工程分析的基础。对于性质复杂的钢筋混凝土结构, 材料非线性与几何非线性常同时存在, 用传统的方法来分析和描述难度非常大,用有限元等工程软件进行钢筋混凝土的力学行为的模拟分析,对于结构设计的合理性与经济性非常有意义。
1 Solid65单元
1.1 单元简介
通常钢筋混凝土结构有限元分析的单元分为两种:杆系单元和实体单元。前者着重分析单元力(包括力和弯矩)与位移(包括位移和转角)之间的关系,而后者着重分析单元的应力-应变关系。单元类型的选取应兼顾计算规模、材料模型的精度等多方面的因素。对于全结构规模较大,可将结构离散成杆系单元进行分析。对于复杂区域(梁柱节点)或重要的构件等可将杆系结构计算的力和位移施加到实体单元模型上,分析局部应力和应变。在结构分析中应尽可能多地采用三维实体单元模型,力求最大程度地真实模拟实际结构构件。
Solid65单元用于含钢筋或不含钢筋的三维实体
模型。该实体模型可具有拉裂与压碎的性能。在混凝土的应用方面,如用单元的实体性能来模拟混凝土,而用加筋性能来模拟钢筋的作用。当然该单元也可用于其它方面,如加筋复合材料(如玻璃纤维)及地质材料(如岩石)。该单元具有八个节点,每个节点有三个自由度,即x,y,z三个方向的线位移;还可对三个方向的含筋情况进行定义。Solid65单元最多可以定义3种不同的加固材料,即此单元允许同时拥有4种不同的材料。混凝土材料具有开裂、压碎、塑性变形和蠕变的能力;加强材料则只能受拉压,不能承受剪切力。
图 1 Solid65单元
1.2 几点假设
(1)只允许在每个积分点正交的方向开裂;
(2)积分点上出现裂缝之后,将通过调整材料属性来模拟开裂,裂缝的处理方式采用分布模
型而非离散模型;
(3)混凝土材料初始时是各向同性的;
(4)除了开裂和压碎, 混凝土也会塑性变形, 常采Drucker- Prager屈服面模型模拟其塑性
行为的应力应变关系。在这种情况下, 一般
在假设开裂和压碎之前,塑性变形已经完成。
1.3 Solid65单元的应用
在实际应用中一般要为solid65单元提供以下数据:实参数,即设定Solid65单元在三维空间各个方向的钢筋材料编号、位置、角度和配筋率等;材料模型,设定混凝土和钢筋材料的弹性模量、泊松比、密度;数据表,给定钢筋和混凝土的本构关系和破坏准则。关于单元几何图形、节点位置、单元坐标系请见图1。单元性质为八节点各向同性材料,单元包括一种实体材料和三种钢筋材料,用命令MAT输入对混凝土材料的定义,而有关钢筋的细则需在实常数中定义,包括材料号、体积率、方向角(THETA, PHI),钢筋的方向角可通过命令/Eshape以图示方式校验。
体积率是指钢筋的体积与整个单元体积的比,钢筋的方向通过单元坐标系中的两个角度(度制)来定义。当钢筋的材料号为0或等于单元的材料号时则不考虑它的作用。
另外,有关混凝土的材料定义,如剪切传递系数,拉应力,压应力都应在数据表中给出,详细描述见表“Solid65混凝土材料数据表”。通常剪力传递系数为0~1.0,0表示平滑的裂缝(完全丧失剪力传递作用),1表示粗糙的裂缝(几乎没有失去剪力传递作用)。这就有利于对裂缝开裂与闭合进行描述。
有关单元荷载的描述见“节点单元荷载”(ANSYS 帮助中专有一节)。压力作为面荷载作用在单元表面如“Solid65的几何模型图”中带圈数字所示。主动力作用在单元内。温度和影响可在节点上作为单元体荷载输入。节点I的温度T(I)默认为TUNIF,如其它节点温度没有被指定,则它们默认为T(I)。对于其它的输入模型未指定温度时默认值都为TUNIF。对影响的设定除用0取代TUNIF外与温度的设定是相同的。
用命令TREF和BETAD分别用来设定整体的基准温度和阻尼值。用MAT命令指定与单元相关的基准温度值(MP,REFT)或阻尼值(MP,DAMP),但不能对钢筋的材料号进行以上定义。
KEYOPT(1)用于设定是否考虑大变形,KEYOPT (5)和KEYOPT(6)则提供是多种单元输出选项(详见单元解答)。
KEYOPT(7)是与是否考虑应力松弛相关的项,当KEYOPT=1时表示考虑,目的是加速裂缝即将开裂时计算的收敛(在混凝土材料数据表的第9个系数中的输入值即为拉伸应力松弛的折减系数)。应力松弛并不能反应因为次生裂缝的产生而引起的应力应变关系的变化。松弛系数在裂缝处为零,因此,相应的开裂面上的刚度也是零。
在几何非线性分析时可用“SOLCONTROL,,,INCP”命令设定考虑抗压刚度的影响。抗压刚度的影响在线性屈曲分析中会被自动考虑。
2 混凝土与钢筋的组合
整体式模型,直接利用Solid65提供的实参数建立钢筋模型;分离式模型,位移协调使用杆件单元模拟钢筋,混凝土和钢筋共用节点;分离式模型,界面单元在杆件单元和混凝土单元之间加入界面单元进行模拟。
在条件允许的情况下,优先使用整体式模型滑移影响可以通过折减钢筋弹性模量加以模拟,推荐折减量:65~80%;在类似节点受往复荷载作用的问题,由于滑移严重,必须使用界面单元;预应力考虑预应力损失,必须使用界面单元。
3 ANSYS建模中需要注意的问题3.1 单元尺寸大小
基于最大开裂应力准则,单元越细,应力集中越严重,开裂出现越早。
解决方法:使用半脆性裂缝模型,减小单元尺寸影响;控制网格大小,单元尺寸不宜小于5cm;控制网格划分,在容易出现应力集中的部位要避免过小的单元出现。
3.2 施加支座
支座是个非常严重的应力集中部位,尽量避免把约束直接施加在支座上。