纤维模型的描述(第二稿).docVIP

1、纤维模型简介 纤维模型的主要思路是沿单元纵向将各控制截面离散化为若干个纤维（包括混凝土纤维和钢筋纤维），同时忽略剪切变形和钢筋粘结滑移的影响，认为整个截面符合平截面假定，在假定截面上每根纤维应变分布均匀、处于单轴应力应变状态的基础上，可根据相应纤维材料的单轴应力应变关系来计算整个截面的力与变形的关系。与直接基于截面的恢复力模型相比，根据各纤维材料的应力应变关系来确定整个截面的力与变形关系的做法，能更为客观、真实地模拟截面的实际受力性能，特别是对模拟变化轴力、弯矩相互作用下柱子的非线性受力，以及双向弯曲和变轴力共同作用下柱子的空间地震反应方面具有明显的优势。此外，通过记录和追踪截面关键点处纤维材料的应力应变关系，可以对构件和结构的非线性反应规律有更为深入的了解和把握。纤维模型的计算精度与截面划分的纤维数目和纤维的划分方式紧密相关。一般来说，可采用较为均匀的划分方式，常见的截面纤维划分如图1所示。由图2.1 可见，截面主要离散化为核心区混凝土纤维（按有约束混凝土材料考虑）、外围混凝土纤维（按无约束混凝土材料考虑）和钢筋纤维。对于纤维划分数量的问题：当截面划分的纤维达到一定数量时，数值积分方法产生的误差将不再显著。 图 1 纤维模型示意图 每一根纤维的轴向应变与截面应变之间的关系如下： —表示第根纤维的应变 —表示第根纤维在构件截面的局部坐标 —表示纤维截面在构件中的位置 —分别表示截面绕轴的曲率 —表示截面的轴向应变 2、考虑剪切效应的模型化方法的提出 如上所述，虽然纤维模型可以模拟构件的弯曲和轴向效应，从理论上完全解决了双向弯曲和轴向效应的耦合，但它不能模拟构件的剪切效应。为了模拟构件的剪切变形，OpenSees 允许通过另外定义一种专门用来模拟构件截面剪切应变的材料，然后将此材料组合到纤维截面中。这一方法的特点是直接在截面层次上加入该剪切材料的，其形象的示意可参见图2，（b）中采用三折线描述应力—应变（或力—位移）关系，当第三点的纵坐标值大于第二点的纵坐标值时，在第三点后，保持切线刚度不变；若第三点的纵坐标值小于或等于第二点的纵坐标值，则在第三点后，材料切线刚度为零。 图 2 考虑剪切效应的模型化方法示意图 Fig 2 Discretization of RC section into fibers (a) 截面组合示意图 (b)剪切材料本构关系示意图 图 3 组合截面及剪切材料本构关系示意图 Fig 3 Aggregated section and the constitutive relation of shear material 如上所述，截面弯曲和轴向效应互相耦合，截面弯矩、轴力均为轴应变、曲率的一个函数。而本文定义的截面剪力与剪应变关系忽略了截面弯曲和轴向效应对其的影响，只与截面剪应变有关，故截面剪力仅是剪应变的一个函数。如式(1)示。 （1） 截面力与应变的向量如式(2)示，式中为弯矩，为轴力，为剪力。 （2） 截面的刚度矩阵可以写成： （3） 其中，受单轴循环本构关系影响，由截面层次上的混凝土、钢筋纤维刚度积分得到。代表剪切刚度，它是剪应变的一个函数。剩余其它参数为 0，从截面的刚度矩阵定义中也可以看出，模型未考虑弯曲、轴向效应与剪切效应的耦合。 弯曲、轴向刚度可以用传统方式得出，截面被分为个纤维，其中是截面中性轴到第个纤维中性轴的距离，是纤维面积，参数详细示意参见图4，是材料的切线模量。刚度计算式参见式4。 （4） 图 4 纤维模型示意图 Fig 4 Discretization of fibre model 为剪切滞回曲线的斜率，计算式参见式5。 （5） 3、有限元模型的建立 3.1混凝土本构模型 本文所介绍的混凝土本构关系模型如图5所示，混凝土受压时的骨架曲线采用Scott等人扩展后的Kent-Park模型。修正的Kent-Park模型是在简化与精确之间的一种较好的平衡，该模型通过修改图中混凝土受压骨架的峰值应力应变以及软化段斜率来考虑横向箍筋的约束影响。其受压骨架曲线分以下三段加以描述： 图5 混凝土单调加载骨架曲线 Fig 5 Monotonic envelopes of concrete （6） （7） 其中 （8） （9） —代表考虑约束所引起的混凝土强