基于ansys的砌体结构非线性有限元分析.docxVIP

基于ansys的砌体结构非线性有限元分析 现代水泥结构注重高层建筑和地震的发展。为了适应这种发展需要,应加快对现代砌体结构理论的研究。毋庸质疑,有限元分析已经成为研究的重要途径之一。近年来随着一些大型有限元分析软件的开发应用,砌体结构的有限元分析得到越来越广泛的重视,不少研究工作者开展了用ANSYS进行砌体结构有限元分析的研究,获得了一些有益的成果。综观这些研究可以发现,由于砌体结构非线性有限元模拟的复杂性,目前在模型选取、材料本构、参数确定、迭代算法、收敛准则等有限元分析的关键问题方面尚未总结出一般性的规律和结论可供参考。为此,对ANSYS应用于砌体结构非线性有限元分析的适用性展开研究,讨论上述问题,并结合算例分析给出相应建议,供相关研究参考。 1 砌体分离式材料模型 砌体是由砌块和砂浆组成的二相复合材料,对其进行有限元分析常有两类模型:分离模型和整体连续体模型。分离模型即是将砌块和砂浆分别建模,可采用两种处理方式:一种是不考虑砖与砂浆之间的粘结滑移,将砌块与砂浆接触面的所有节点自由度耦合在一起;另一种是考虑砖与砂浆之间的粘结滑移,将砌块与砂浆通过接触单元或非线性弹簧单元联系在一起。由于砖与砂浆之间的粘结滑移关系曲线目前研究尚不成熟,接触面的水平粘附强度常常难以得到,目前使用前一种方式的较多且较易实现。整体连续体模型即是将砂浆和砌块作为一个整体来考虑,连续体单元的材料参数可通过以下三种方式获得:一是试验实测;二是按规范取值,规范无规定时则根据现有的砌体本构关系和材料性质的有关研究选用;三是根据由复合材料力学发展起来的均质化理论,将所有组成材料的几何和本构信息融入到一个代表性体积单元(RVE)中,按计算确定RVE的等效材料属性。前两种均以试验为基础,较符合实际,且容易获得,应用方便;第三种计算RVE等效材料参数的过程相当繁琐,且需已知砌块和砂浆的力学参数,其推广应用受到了限制。 总体来说,分离模型可以模拟砌块与砂浆之间的作用和砌体破坏机理,适用于模拟小型试验砌体的破坏行为,但计算量大,建模繁琐;整体连续体模型适于分析大规模的墙体,但对于详细的应力分析和理解砌体多样的失效机理却显得不足。采用哪种模型一般应视分析目的而定。 对于砌体组合结构,其中还包含有钢筋混凝土构件,其有限元模型可采用分离钢筋模式和整体分布钢筋模式。前者可考虑粘结滑移。 ANSYS分析时,砌体分离模型的砌块和砂浆可用三维实体单元SOLID45或SOLID65模拟,二者具有处理塑性、徐变、膨胀、应力刚化、大变形和大应变的能力,SOLID65还可以模拟开裂和压碎以及可以考虑三个方向的配筋;接触单元可采用CONTACT174等;非线性弹簧单元可用COMBIN39等。砌体整体连续体模型可用SOLID45或SOLID65,也可用三维实体壳单元SHELL63模拟。钢筋混凝土分离钢筋模式可用SOLID65+LINK8单元,整体分布钢筋模式可用带配筋的SOLID65单元。 2 砌体单轴应力应变曲线的确定 砌体是一种非均质的、各向异性材料,砌缝是其薄弱环节,对其本构关系及破坏准则进行研究非常困难。国内外已开展的大量研究中大多限于基于试验结果的平面受力砌体破坏准则研究,至今仍没有一个被广泛认可的砌体应力应变关系及其破坏准则。进行非线性有限元分析,合理选取材料的弹塑性本构关系是至关重要的,这正是进行砌体结构非线性有限元分析的难题。目前ANSYS尚没有能够反映砌体受砌缝影响的复杂的剪压破坏模型,只有参照类似材料的常用破坏准则,通过参数的适当选取来最大限度地模拟砌体的破坏。 弹塑性本构关系通常由三部分组成:初始屈服准则、流动法则和强化规律。初始屈服准则定义了多轴应力状态下的弹性极限。砌体、混凝土和岩石等材料为颗粒状,受压强度远大于受拉强度,且屈服应力受静水应力影响明显。砌体材料的屈服准则可以考虑使用岩土工程中广泛使用的Drucker -Prager(DP)准则或ANSYS专门为钢筋混凝土单元SOLID65开发的CONCRETE材料破坏准则。 DP屈服准则是修正的Von Mises准则,用以逼近Mohr-Coulomb准则,如图1所示。DP准则假设材料为理想弹塑性,无强化规律;流动法则可以使用相关流动法则,也可以使用不相关流动法则。需输入的三个参数即粘聚力C、内摩擦角φ和膨胀角φf,一般由试验确定或根据经验选取。对于砌块、砂浆和砌体而言,C、φ值的选取目前还没有相关的试验资料和经验取值,但可以根据其与材料的单轴受拉屈服应力σt和受压屈服应力σc之间的关系来确定: φ=sin-1[3√3β2+√3β]?β=σc-σt√3(σc+σt)(1)C=σy√3(3-sinφ)6cosφ?σy=2σcσt√3(σc+σt)(2) 式中:σt和σc可以取试验实测值或根据其强度等级按规范选取。 破坏面是