利用ANSYS Solid 65单元分析复杂应力条件下的混凝土结构

简介： 在大型有限元软件ANSYS中，为混凝土材料专门定义了一种单元：Solid 65。它不但可以模拟混凝土材料特有的开裂，压碎等力学现象，而且预先定义好了混凝土的破坏准则,为使用者提供了很多方便。在本文中，作者结合自己的科研工作，通过4个科研项目实例，介绍了Solid 65单元在框架、叠合梁、约束混凝土柱、组合结构节点等问题中的应用。详细讨论了Solid 65单元和钢筋、型钢、钢管、碳纤维等不同材料共同工作时建模的方法及相互连接的处理。这些计算实例表明，在准确建立有限元模型并合理选择各种材料组合方法的前提下，利用ANSYS可以对各种混凝土结构进行准确的模拟，从而可以减小试验工作量，节省开支。但是，如果有限元模型建立不当，也会造成错误的结果，本文中对这些错误进行归类分析，提出了解决方法。
关键字：有限元 混凝土 ANSYS

一、引言

　　混凝土是目前应用最为广泛的建筑材料之一。它可以很方便地与钢筋、型钢、钢管以及其他一些结构材料组合成多种多样的构件和节点，应用于各种结构形式中。为了了解混凝土结构的详细受力机理和破坏过程，往往需要利用三维实体单元进行非线性有限元分析。而混凝土本身同时具有开裂、压碎、塑性等诸多复杂力学行为，在三维条件下这些力学行为更加难以确定，给实际应用带来了较大的困难。为了便于使用者应用，ANSYS软件内部设定了专门面对混凝土材料的三维实体单元形式SOLID 65。并建立了三维情况下混凝土的破坏准则，提供了很多缺省参数，从而为使用者提供了很大的方便。此外，ANSYS软件本身所拥有的大量单元形式，可以很方便的让使用者建立混凝土和其他材料之间的共同工作模型，因此在很多实际问题中都取得了成功应用。本文将结合作者实际参加的一些研究项目来说明其具体应用。

　　二、SOLID 65单元的使用方法

　　SOLID 65单元本身包含两部分。一是和一般的8节点空间实体单元SOLID 45相同的实体单元模型，但是加入了混凝土的三维强度准则。二是由弥散钢筋单元组成的整体式钢筋模型，它可以在三维空间的不同方向分别设定钢筋的位置，角度，配筋率等参数。在实际应用中，一般需要为SOLID 65 单元提供以下数据：

　　1）实参数real constants；在实参数中给定SOLID 65单元在三维空间各个方向的钢筋材料编号，位置，角度和配筋率。对于墙、板等钢筋分布比较密集而又均匀的构件形式，一般使用这种整体式钢筋混凝土模型。

　　2）材料模型Material Model；在这里设定混凝土和钢筋材料的弹性模量，泊松比，密度。

　　3）数据表 Data Table；在这里给定钢筋和混凝土的本构关系；对于钢筋材料，一般需要给定一个应力应变关系的Data Table，譬如双折线等强硬化或随动硬化模型等。而对于混凝土模型，则需要两个Data Table。一个是本构关系的Data Table，比如使用Multilinear kinematic hardening plasticity模型或者Drucker-Prager plasticity 模型等，用来定义混凝土的应力应变关系。另一个则是SOLID 65特有的Concrete element data，用于定义混凝土的强度准则，譬如单向和多向拉压强度等等。

　　本文附录中给出了定义混凝土材料的log文件流。

　　三、混凝土与其他材料的组合

　　一）与钢筋的组合

　　混凝土与钢筋组合是最常见的一种组合方式，一般说来，可供选择的方法有以下三种：

　　1）整体式模型

　　直接利用SOLID 65提供的实参数建模，其优点是建模方便，分析效率高，但是缺点是不适用于钢筋分布较不均匀的区域，且得到钢筋内力比较困难。主要用于有大量钢筋且钢筋分别较均匀的构件中，譬如剪力墙或楼板结构。

　　2）分离式模型，位移协调

　　利用空间杆单元link 8建立钢筋模型，和混凝土单元共用节点。其优点是建模比较方便，可以任意布置钢筋并可直观获得钢筋的内力。缺点是建模比整体式模型要复杂，需要考虑共用节点的位置，且容易出现应力集中拉坏混凝土的问题。

　　3）分离式模型，界面单元

　　前两种混凝土和钢筋组合方法假设钢筋和混凝土之间位移完全协调，没有考虑钢筋和混凝土之间的滑移，而通过加入界面单元的方法，可以进一步提高分析的精度。同样利用空间杆单元link 8建立钢筋模型。不同的是混凝土单元和钢筋单元之间利用弹簧模型来建立连接。不过，由于一般钢筋混凝土结构中钢筋和混凝土之间都有比较良好的锚固，钢筋和混凝土之间滑移带来的问题不是很严重，一般不必考虑。

　　二）与型钢及其他材料的组合

　　一般当混凝土和型钢相组合时，由于型钢截面尺寸大，界面上剪力很大，且型钢界面一般较光滑，因此，型钢和混凝土之间位移不协调情况比较显著，需要考虑。在ANSYS，有以下建模方法

　　1）利用非线性弹簧单元Combin 39

　　Combin 39 弹簧单元可以很方便的设定界面上的剪力－滑移关系。并且有多种滑移模型可供选择，是一种很实用的界面单元形式。

　　2）利用接触单元

　　当混凝土和型钢表面可能发生脱离，界面上又没有设置剪力连接件的情况下，混凝土和型钢之间的接触就需要考虑。ANSYS软件内部提供了点与点接触，点与面接触，面与面接触等多种接触分析模型。在综合对比了各种接触模型后，作者推荐使用空间点对点接触分析。因为在实际混凝土结构中，型钢和混凝土之间的脱开缝隙比较小，相对滑移也比较小，用点对点接触分析完全可以模拟二者之间的相互关系。而其效率和收敛速度都要高于其他的接触形式。

　　四、实例分析

　　4.1、异型柱框架结构

　　某带斜支撑异型柱框架如图1所示，在本例中，因为需要了解异型柱框架在外力作用下的整个荷载位移关系，及混凝土和钢筋应力变化及分布情况，因而需要建立混凝土的空间实体模型。同时考虑到异型柱中钢筋分布都沿着混凝土的外边缘，且钢筋应力是我们关心的重点，因此，我们选择位移协调的分离式钢筋模型。用link 8单元代表结构中的钢筋。有限元分析正确模拟了斜支撑对框架的刚度和承载力的影响，说明了斜支撑在改变框架受力机理中的作用，并提出了设计建议。最后得到框架柱和斜支撑中钢筋应力分布如图2所示。

　　4.2、预应力叠合梁结构

　　某预应力型钢－混凝土叠合梁如图3所示，其中型钢和混凝土面板是通过栓钉相连接。考虑到上部混凝土板内钢筋分布较均匀，因此使用了整体式钢筋混凝土模型，直接在实参数中定义配筋率，而连接型钢和混凝土的栓钉由于抗剪刚度比较小，因此应该考虑其滑移的影响。因此我们选用了非线性弹簧单元Combin39来模拟栓钉的影响，并设定栓钉的轴向刚度为，剪切－滑移关系为，这里，为栓钉弹模，为栓钉截面积，为栓钉长度，为滑移量，，为混凝土的弹性模量和抗压强度，为栓钉强度。计算得到叠合梁的荷载位移曲线如图4所示。得到剪力连接件的滑移度对整个构件的刚度的影响。