《钢筋混凝土结构开裂计算方案》由会员分享,可在线阅读,更多相关《钢筋混凝土结构开裂计算方案(7页珍藏版)》请在金锄头文库上搜索。
1、钢筋混凝土构造开裂计算方案1概述 开裂计算是工程中比较关心的问题,但一直是有限元分析的一个难点,涉及到材料本构、计算收敛性等诸多问题。ANSYSCivilFEM提供了钢筋混凝土构造开裂计算功能,其中土木专用模块CivilFEM提供的非线性混凝土计算适用于混凝土梁构造的非线性计算包括开裂,可以直接通过截面定义钢筋,从而模拟钢筋混凝土梁。但对于更一般的构造,用梁单元来模拟不一定适宜,需要采用更一般的单元,ANSYS提供了专用的钢筋混凝土实体单元SOLID65来模拟钢筋混凝土构造,该单元材料采用混凝土材料模型,可定义混凝土的开裂、压碎准那么。另外可以定义钢筋方向和体积率,可用来模拟钢筋混凝土的破坏。
2、本文将通过算例对ANSYSCivilFEM开裂计算的效果进展讨论,并针对一些计算难点提出初步的解决方案。2CivilFEM开裂计算 CivilFEM适宜于梁构造开裂分析,另外为了与后面SOLID65单元开裂计算结果进展比较,先讨论了CivilFEM的开裂计算。 CivilFEM开裂计算需要考虑的要点: 1、激活CivilFEM非线性模块CFACTIV,NLC,Y,这是CivilFEM非线性计算的前提。 2、即使事实上为小变形,也必须翻开几何非线性效应NLGEOM,ON,否那么无法激活非线性迭代。 3、通常应该关闭求解控制SOLCONTROL,OFF,由于CivilFEM非线性计算通过修改实常数
3、的等效方法,自动求解控制反而可能导致发散。 4、在收敛不好的情况下,可以增加子步数、翻开自动步长AUTOTS,ON或可以给定一个比较大的迭代数NEQIT,NUM,以改善收敛,线性搜索有时也可以改善收敛LNSRCH,ON。 5、有些情况下上述调整可能仍然无法保证收敛,这通常发生在一些开裂、受压区状态转换的临界点,尤其在动力分析中更易出现,可以结合两个方法抑制,一是放松收敛准那么CNVTOL,开裂分析状态变化剧烈,往往是接近收敛但出现振荡,放松收敛可以保证在较松的准那么下收敛,但可得到足以满足要求的结果。另一个方法是在未收敛情况下仍然继续下一步计算NCNV,0。不收敛往往发生在一些临界点,该命令可
4、以保证跳过这些点,而后续载荷步往往可以迅速收敛,只要构造事实上具有平衡状态,没有失效,那么这种处理不会影响到总体结果,后面的动力分析实例也可说明这一点。 以悬臂梁为例,该悬臂梁长10m,截面如图1所示,尺寸为0.6m0.5m,钢筋直径20mm,混凝土保护层厚40mm。混凝土参数国际单位制:弹模E=28.848E9, 泊松比m=0.2, 密度D=2600,钢筋参数:E=200E9, m=0.3, D=7800。2.1CivilFEM开裂静力分析图1 截面 计算采用的命令流为文件crack_static_cv.txt。图2为CivilFEM定义的截面,计算采用梁单元beam54在CivilFEM中
5、beam54、beam44梁单元可以进展非线性计算,CivilFEM在定义梁单元截面后自动计算beam54单元的实常数,无需用户输入。图3为实际形状显示的悬臂梁模型。图2 截面定义图3 有限元模型 梁一端固支,一端施加Y向力,考虑图4所示可变载荷,最大为1500N,最小为-1500N,计算变化载荷作用下梁的开裂,以验证CivilFEM开裂计算可以考虑这种交变载荷情况。 计算没有考虑混凝土抗拉强度,并进展静力分析。图4 可变载荷 图5为载荷到达1500N时固端截面混凝土部分的正应力,红色区域即为开裂区。图6为相应的钢筋正应力。混凝土不抗拉,受压区混凝土最大压应力为-0.74MPa,受压钢筋应力为
6、-3.60MPa,受拉钢筋应力为15.27MPa, 拉裂区应力全部由钢筋承担。图6 固端钢筋正应力载荷:1500N图5 固端混凝土正应力载荷:1500N 图7为载荷到达-1500N时固端截面混凝土部分的正应力,图8为相应的钢筋正应力。受压区混凝土最大压应力为-0.81MPa,受压钢筋应力为-3.41MPa,受拉钢筋应力为22.40MPa, 拉裂区应力全部由钢筋承担。图8 固端钢筋正应力载荷:-1500N图7 固端混凝土正应力载荷:-1500N 从计算结果来看,载荷为1500N时,梁下部拉裂,载荷为-1500N时,梁上部拉裂,由于上部配筋数少于下部,所以上部拉裂时的下部混凝土压力要大于下部拉裂时
7、上部混凝土压力。且开裂后混凝土仍然具有抗压才能。这说明CivilFEM开裂计算可以考虑交变载荷作用。 图9为载荷1500N 时开裂计算得到的Y向位移图,梁端最大为0.0063m,图10为不考虑开裂的Y向位移图,梁端为0.0017m不考虑开裂的命令流为static_cv.txt,可以看到开裂使得位移增加了很多。图10 不考虑开裂位移载荷:1500N图9 考虑开裂位移载荷:1500N2.2CivilFEM开裂动力分析 计算采用的命令流为文件crack_dynamic_cv.txt,悬臂梁载荷为图11所示的地基加速度历程,计算时间为1秒。 图12 为在这个动力载荷作用下梁端点的Y向位移历程,图13为
8、同一模型不考虑开裂情况下梁端点的Y向位移历程不考虑开裂的命令流为dynamic_cv.txt。图11 加速度历程 计算结果非常合理,开裂后不仅位移大幅增加,而且由于刚度变小,导致构造的振动周期延长。图13 梁端位移历程不考虑开裂图12 梁端位移历程考虑开裂3实体单元SOLID65开裂计算 对于不适于用梁单元模拟的实体或其它钢筋混凝土构造,需要采用ANSYS的钢筋混凝土单元SOLID65来模拟,SOLID65单元可以定义开裂、压碎强度,也可以定义分布钢筋,但分布钢筋由于将钢筋均质化于单元内,缺乏以反映实际情况,可以用梁单元或杆单元来模拟离散钢筋,用约束方程CEINTF来建立钢筋节点与混凝土单元的
9、位移协调关系。本例采用BEAM188来模拟钢筋。 SOLID65单元开裂计算需要考虑的要点: CivilFEM中考虑的一些事项仍然适用: 1、关闭求解控制SOLCONTROL,OFF,根据测试结果,由于刚度变化剧烈,自动求解控制情况下,收敛参数不易操纵,有时更易导致发散。 2、在收敛不好的情况下,可以增加子步数、翻开自动步长AUTOTS,ON或可以给定一个比较大的迭代数NEQIT,NUM,以改善收敛,线性搜索有时也可以改善收敛LNSRCH,ON。 3、放松收敛准那么CNVTOL,以及在未收敛情况下仍然继续下一步计算NCNV,0,这两个方法非常有效而且可以得到可信的结果,详细讨论见第2节。 此外
10、,SOLID65本身有一些参数可以增强收敛: 4、指定极小的分布钢筋体积率通过单元实常数,比方1e-6,这不会影响计算结果,但可以使得开裂后单元具有一个小刚度,不致奇异,从而增强收敛。 5、开裂单元指定适当的剪力传递系数通过混凝土材料,一般张开裂纹0.1,闭合裂纹1.0,可以极大地进步收敛性。 6、指定开裂起始刚度松弛因子为1.0通过混凝土材料,并设SOLID65单元的KEYOPT(7)=1,可以使得开裂后刚度逐渐减小至0,增强收敛。 7、给混凝土指定合理的抗拉强度通过混凝土材料,可以极大地进步收敛性。 8、为了得到一个好的结果,在构造主要受弯厚度方向单元不要太少,不要少于4层,8层以上比较好
11、如本例为8层。 用SOLID65模拟混凝土,BEAM188单元模拟钢筋,混凝土给定抗拉强度0.1MPa,对前述算例进展开裂分析,并进展比较。3.1SOLID65开裂静力分析 计算采用的命令流为文件crack_static_65.txt,图14为计算模型,钢筋节点与混凝土节点之间建立位移约束方程,梁端部建立刚性区来施加载荷。图15为钢筋分布图,混凝土采用透明来直观显示构造。图14 计算模型图15 钢筋分布 图16为载荷到达1500N时混凝土梁轴向正应力,固端红色区域中除了小部分为小于抗拉强度的受拉区外,大部分为开裂区。图17为相应的钢筋正应力。受压区混凝土最大压应力为-0.78MPa,受压钢筋应
12、力为-3.60MPa,受拉钢筋应力为15.7MPa, 拉裂区应力全部由钢筋承担。与前面CivilFEM的计算结果比较可以看出,二者是非常接近的,虽然本例中混凝土指定了0.1MPa的抗拉强度,但由于抗拉强度很小,所以计算结果变化不大,但抗拉强度的存在却可以大大进步SOLID65的收敛性而且更接近实际。图17 钢筋轴向应力载荷:1500N图16 混凝土轴向应力载荷:1500N 图18为载荷到达-1500N时混凝土梁轴向正应力,图19为相应的钢筋正应力。受压区混凝土最大压应力为-0.89MPa,受压钢筋应力为-3.62MPa,受拉钢筋应力为23.4MPa, 拉裂区应力由钢筋承担。计算结果同样与Civ
13、ilFEM非常接近。图19 钢筋轴向应力载荷:-1500N图18 混凝土轴向应力载荷:-1500N 图20为载荷1500N 时开裂计算得到的Y向位移图,梁端最大为0.0062m,图21为不考虑开裂的Y向位移图,梁端为0.0017m将抗拉强度设置为无穷大即可不考虑开裂,结果与CivilFEM一样,开裂使得位移增加了很多。图21 不考虑开裂Y向位移载荷:1500N图20 考虑开裂Y向位移载荷:1500N SOLID65单元可以绘图表示开裂区,图22、23为载荷1500N和-1500N 时开裂区分布。加载过程中开裂区域的扩展过程见动画crack_static.avi。图22 开裂区域载荷:1500N
14、图23 开裂区域载荷:-1500N3.2SOLID65开裂动力分析 对第2.2节同样的问题用SOLID65和BEAM188单元进展动力分析,计算采用的命令流为文件crack_dynamic_65.txt。 图24 为在动力载荷作用下梁端一点的Y向位移历程,图25为同一模型不考虑开裂情况下梁端一点的Y向位移历程将抗拉强度设置为无穷大即可不考虑开裂。 计算结果与前面CivilFEM结果很接近,只是由于考虑了抗拉强度,所以开裂后的峰值位移要稍小一些。图25 梁端一点位移历程不考虑开裂图24 梁端一点位移历程考虑开裂 图26为时间为1秒时开裂分布图,图27为开裂状态,灰色部分没有开裂,其余为开裂部分,不同颜色代表不同张开度,红颜色表示完全闭合,篮颜色表示最大张开。振动过程中开裂区域的扩展过程见动画crack_dynamic1.avi和crack_dynamic2.avi。图27 开裂状态时刻1秒图26 开裂区域时刻1秒4结论