ANSYS软件在砌体结构非线性分析中的应用.docVIP

安阳工学院学报 第六期 JournalofAnyangInstituteofTechnology 69 ANSYS软件在砌体结构非线性分析中的应用 宋扬1，高均昭2，赵军1 （1.安阳工学院建筑工程系，河南安阳455000；2.许昌学院城市与环境学院，河南许昌461000） 摘 要：利用ANSYS软件建立空间模型，对轴心受压时设置构造柱的砖砌体结构进行了模拟，得到的结果与试验结果 较吻合，并优于以往平面有限元的计算结果。说明可以利用ANSYS软件强大的模拟功能对砌体结构进行有限元分析，为目前日益增多的各类砌体结构提供非线性分析依据。 关键词：ANSYS非线性有限元砌体结构中图分类号：TU312 文献标识码：A 文章编号：1673-2928（2009）06-0069-03 为了环境保护，可持续发展，建筑节能等，我国鼓励使用蒸压蒸养砖（灰砂砖、粉煤灰砖和炉渣砖等）和砌块（普通混凝土小型砌块，轻骨料混凝土空心砌块和加气混凝土砌块等），2010年底，所有城市城区将禁止使用实心粘土砖，全国实心粘土砖产量控制在4000亿块以下，为适应这种要求，各种新型砌体材料应运而生。在已知物理、力学性质的前提下如何较准确、方便地对其进行承载力以及非线性分析，是工程结构研究、设计面临的一个问题。本文利用ANSYS通用有限元软件对湖南大学所做试验中的模型进行模拟，并将计算结果同试验结果和原平面有限元分析结果进行比较。 2）首先将荷载步（Loadstup）时间定为1，大概 计算出破坏荷载，再将其分为若干个荷载步逐级加载，这样就可以较准确地找到破坏荷载。但采用弧长方法时子步（Substupe）和时间数（Time）不是按比例发展，所以不能用它们计算对应荷载。 图1应力－应变曲线控制点 3）选用Newton-Raphson迭代法（FULLNR）。 如果结合弧长法，计算时间会增加，但通过比较证实结果要更精确许多。需注意的是如果不选择弧长法，也不要人为地固定子步数，否则程序只能按指定的子步计算，Newton-Raphson的优点不能发挥，结果误差较大。为了增大计算精度，本文采用了FULLNR方法并结合了弧长法。 1用ANSYS软件求解结构非线性问题 ANSYS软件强大的功能已经在结构分析中得 到广泛应用，但将其应用到砌体结构中的却几乎没有，本文尝试利用ANSYS软件对砌体结构进行模拟。首先我们可以认为当程序计算到刚刚不收敛时上一子步所对应的荷载为破坏荷载，然后利用子步比例及对应的应力（或应变）来判断开裂荷载。当然，这需要多次尝试，步步逼近。其中最关键的问题是如何对非线性问题（包括材料非线性和接触非线性）的求解进行控制。经过数百次的计算比较，发现如果计算参数或计算方法选择不当，结果会产生较大的偏差，而同时有些选择只会占用过多内存，却对计算精度影响不大。关键的控制点有： 4）本文根据文献[1]的建议，以节点力作为收 敛的控制标准，并选用范数（L2norm）来判断。至于收敛允许值，通过比较发现取0.1％（程序默认）比1％所耗的时间多两倍以上，但计算结果差别很小，建议取1％即可。 5）利用柔体－柔体的面－面接触单元模拟了混 凝土与砌体之间的连接，利用接触单元的各种属性设置代替实际结构中的马牙槎、拉筋，其中摩擦系数按文献[2]取为0.6。 1）材料非线性分析的前提是要有较准确的应 力－应变曲线，但如图1所示，控制收敛的关键是上升段和应力峰值及对应应变。 收稿日期：2008-12-23 6）判断开裂荷载前， 重新定义时间控制命令 作者简介：宋扬（1972-），男，河南省安阳市人，安阳工学院建筑工程系讲师，工学硕士，主要从事建筑结构与岩土工程的 教学与研究， 70 安阳工学院学报 JournalofAnyangInstituteofTechnology 2009年 （TimeContral），固定子步数（Numberofsubstupes）或者时间增量（Timeincrement），这样即可准确地计算出每一个子步所对应的荷载。然后建立一个数据文本文件，用来存贮依据破坏准则编制的函数，从而判断在每一子步对应荷载作用下，构件是否开裂。 大、下部小，且中下部分布相对均匀、差别不大。它沿墙体水平方向是中间大、两端小。其应力峰值位于墙体上部跨中处，它随构造柱间距的减小而减小。图6为单片墙体中上部从外侧到内侧的竖向压应力分布曲线，可以明显地看出中间大、两边小，而且靠近中柱处压应力比靠近边柱处要大一些。中部较均匀分布应力之外的水平距离可视为是构造柱对墙体竖向压应力分布的影响长度。 2设置钢筋混凝土构造柱砌体结构的轴心受压非 线性分析 湖南大学曾对这种结构轴心受压时进行了试验研究和理论计算分析[4][5]，但其计算方法是按平面有限元问题考虑。为了更接近实际情况，本文建立空间模型，为了验证A