1、1钢筋混泥土梁单元双非线性分析的共旋坐标法（1. 湖南大学 土木工程学院，湖南 长沙 410082；2. 长沙理工大学 土木与建筑工程学院， 湖南 长沙 410076） 摘 要：针对已有的钢筋混凝土梁单元非线性分析模型采用较多的假定和近似从而导致计算量增大及计算精度下降的问题，基于共旋坐标法建立了考虑材料和几何非线性的任意截面钢筋混凝土梁的数值分析模型.首先利用虚功原理计算共旋坐标系下完全粘结钢筋混凝土梁考虑材料非线性的切线刚度矩阵，再通过结构坐标系与共旋坐标系下节点力之间及节点位移之间的总量关系及微分导出的增量关系，最终获得钢筋混凝土梁在结构坐标系中考虑几何与材料双重非线性的切线刚度矩阵.算
2、例结果表明，本文算法可减少计算量、不累积误差、精度高. 关键词：钢筋混凝土； 梁单元；共旋坐标法；双非线性；微分法；切线刚度矩阵 中图分类号：TU323.3 文献标识码：ACorotational Procedure for the Binonlinear Analysis of Reinforced Concrete Beam Element DENG Jihua1，2， SHAO Xudong1 （1. College of Civil Engineering， Hunan Univ， Changsha， Hunan 410082， China； 2. School of Civil En
3、gineering and Architecture， Changsha Univ of Science and 2Technology， Changsha， Hunan 410076， China） Abstract：Multiple assumptions and approximations in nonlinearity analysis models of existing reinforced concrete beam element result in low calculation efficiency and low calculation accuracy. In thi
4、s paper， based on corotational procedure， a numerical model for a given section considering material and geometrical nonlinear analysis of reinforced concrete beam element was developed. Firstly， by means of virtual work， a tangent stiffness matrix for the material nonlinearity of perfectlybonded re
5、inforced concrete beam element was derived in corotational coordinate system. Then， by building total and incremental relationships derived from differential equations of nodal displacements and forces between global coordinate system and corotational coordinate system， respectively， tangent stiffne
6、ss in global coordinate system reinforced concrete beam element was developed by considering geometric and material nonlinearity. A comparison between the results in this paper and those from existing references has demonstrated that the algorithm developed is highly efficient and accurate with many
7、 advantages， such as noncumulative calculation errors and reduction in computation. Key words：reinforced concrete； beam element；corotational 3procedure；binonlinear；variational method；tangent stiffness matrix - 在目前已有的钢筋混凝土杆系结构考虑几何与材料双重非线性的各种分析模型中，分层梁模型由于具有不受截面形状限制、不同梁层可采用不同材质、能模拟已开裂的梁结构等优点而应用较广1-3.但目
8、前的分层梁模型存在以下缺点：1）在分层梁模型中，将钢筋划分成截面的一层，认为截面每一层的应变沿梁轴向均匀分布.因此，为保证计算精度，单元长度需划分得很短，这在非线性计算中是非常不利的；2）将钢筋沿单元轴向理想化为平行于梁轴线的直线段，而在实际钢筋混凝土梁中，由于受力或构造的需要，钢筋并不总是完全平行于梁轴线；3）材料非线性分析中，一般通过假定截面形心处的应变和曲率得到各层的应变，由材料本构关系得到应力，再由截面平衡得到计算内力，将其与实际内力比较以确定单元刚度，这一过程往往需要多次反复迭代，甚至由于混凝土和钢筋的应力 - 应变曲线为分段曲线而不收敛4；4）几何非线性一般通过应变计算中计入高阶项
9、，基于总体拉格朗日法（T.L 法）或修正拉格朗日法（U.L 法）来考虑，这不仅导致单元切线刚度矩阵异常复杂，还具有 T.L法和 U.L 法本身所固有的缺点，如 T.L 法在非线性程度较高时所得计算结果精度很差5，以及 U.L 法为保证计算精度，需将荷载步取得较4小，导致计算量显著增加和误差累积的问题6. 为解决上述问题，须对现有非线性算法进行改进.在几何非线性分析方法中，共旋坐标法相对于 T.L 法和 U.L 法而言，具有列式简单、力学概念明确以及计算精度高的特点，因而成为研究的热点7-8.本文在上述文献的基础上，首先利用共旋坐标系下应变与扣除刚体位移后的变形呈线性关系的特点计算出完全粘结钢筋
10、混凝土梁考虑材料非线性的切线刚度矩阵，再基于静力平衡通过微分获得钢筋混凝土梁在结构坐标系中考虑几何与材料双非线性的切线刚度矩阵，多个算例表明本文算法是正确的. 湖南大学学报（自然科学版） 2013 年 第 8 期 邓继华等：钢筋混凝土梁单元双非线性分析的共旋坐标法 1 钢筋混凝土梁元切线刚度矩阵 图 1 所示为钢筋混凝土平面梁元，基于实际情况，可设钢筋在梁单元内为直线. 图 2 示意了在结构坐标系 XY 中初始时刻和计算时刻 t 钢筋混凝土梁元的几何参数及位移的即时变量；设钢筋混凝土梁元的共旋坐标系为xy，该坐标系是随单元变形而转动的，它始终以节点 i 为原点，以节点i 到 j 的连线方向为
11、x 轴，由 x 轴逆时针转 90为 y 轴. 图 1 钢筋混凝土梁元 Fig.1 Reinforced concrete beam element 图 2 变形前后梁单元 Fig.2 Beam element before and after deformation 51.1 共旋坐标系下混凝土梁元切线刚度矩阵 设初始时刻单元节点在结构坐标系里的坐标为（xi，yi）和（xj，yj） ，在计算时刻 t 结构坐标系中的位移向量为 d=uiviiujvj jT，在共旋坐标系中的位移向量为 dl=ui （1） 混凝土梁单元在共旋坐标系中的节点位移为： （2） 式中：0l 和 tl 分别为变形前、后梁单
12、元长度. 设单元任一截面形心轴处的应变及截面曲率为 0，在共旋坐标系下单元的应变 - 位移关系只需考虑线性项，有： （3） 式中：B0 为线性应变矩阵. 采用沿梁高分层的方法，对计算截面任意梁层 k，基于平截面假定，由 0 和 可求出该梁层的应变值 k 为： （4） 式中：yk 为梁层 k 到形心轴的距离. 对式（4）微分，有： （5） 通过混凝土的应力 -6- 应变关系由 k 可求得梁层 k 的应力 k 和切线模量 ETk，截面力P=（ncmc）由截面平衡条件求得为： （6） 式中： Ak 为梁层 k 的面积；n 为截面分层总数. 对式（6）微分，有： （7） 将式（5）代入式（7）并写成矩
13、阵形式有： （8） （9） 由虚功原理可建立单元的平衡方程组有： （10） 式中：fc=fxicfyicmicfxjcfyjcmjcT 为共旋坐标系下混凝土单元等效节点力列阵. 将式（9）代入式（10）可得： （11） 对式（11）微分，并考虑式（8）有： （12） 式中：kTc 为共旋坐标系下考虑材料非线性的混凝土元切线刚度矩阵. 1.2 共旋坐标系下钢筋单元切线刚度矩阵 参考图 1 所示，设钢筋上距离左节点所在截面为 x 的任意一点 A 在7竖向距离中性轴为 d（x） ，且 x 与 d（x）的值在变形过程中始终不变.对于钢筋只需考虑轴向拉压应变，与前面混凝土应变计算类似，在共旋坐标系下钢筋
14、的应变计算也只考虑线性项，由应变旋转公式9知钢筋的应变为： （13） 式中：c，A 为混凝土在 A 点的应变. 对式（13）微分，并联立式（3）有： （14） 式中：B0，s 为钢筋的应变矩阵值，只需将钢筋具体位置的值代入B0 就很容易求得. 通过钢筋的应力 - 应变关系由 s 可求得该点的应力 s 及切线模量 ETs，钢筋的轴向力 Ps 为： （15） 式中：As 为钢筋的面积. 对式（15）微分，并联立式（14）有： （16） 由虚功原理可建立共旋坐标系下钢筋单元的平衡方程有： （17） 式中：fs=fxisfyismisfxjsfyjsmjsT 为共旋坐标系下钢筋元等效节8点力列阵. 将
15、式（14）代入式（17）可得： （18） 对式（18）微分，并联立式（16）有： （19） 式中：kTs=cos 30l0BT0，sAsETsB0，sdx 即为共旋坐标系下考虑材料非线性的钢筋元切线刚度矩阵. 1.3 结构坐标系下钢筋混凝土梁切线刚度矩阵 从式（2）知，由于 ui，vi，vj 恒为 0，对式（2）的后 3 项微分，不难得到共旋坐标系下位移微分 dl 用结构坐标系下位移微分d 表达的形式： （20） 设钢筋混凝土梁单元在结构坐标系下的节点力向量F=FxiFyiMiFxjFyjMjT，由静力平衡可知： （21） 式中：t 为坐标转换矩阵. 将式（21）两边微分可得： （22） 对矩
16、阵 t 微分，为表述方便，将 tT（fc+fs）改写成： （23） 为得到 tT（fc+fs）用 d 表达的形式，联立（12） （19） （20）有： 9（24） 联立式（22） （23）和（24） ，可得到结构坐标系下钢筋混凝土梁单元考虑几何与材料双重非线性的单元切线刚度矩阵为： （25） 2 非线性分析流程 1）根据上一计算时刻单元 i，j 节点在结构坐标系下的总位移向量d，由式（2）求出共旋坐标系中的位移向量 dl，由式（11）和（12）求出混凝土梁元在共旋坐标系下的切线刚度矩阵 kTc 及等效杆端力 fc；再基于式（18）和（19）求出钢筋梁元在共旋坐标系下的切线刚度矩阵 kTs及等效
17、杆端力 fs. 2）通过式（25）得到钢筋混凝土单元在结构坐标系下的切线刚度矩阵 KT，基于式（21）得到结构坐标系下的等效杆端力 F. 3）重复 1）至 2）的步骤，生成结构切线刚度矩阵KT 和等效杆端力合力F. 4）计算不平衡力 R=L-F，其中 L 为到计算时刻 t 施加的总外荷载的等效节点力. 5）求解结构方程K?d=R，得到节点位移增量 d，将其叠加到总位移向量 d 中. 6）收敛条件判断，如收敛，转到 t+t 时刻计算，如不收敛，返回1）进行下一次迭代. 3 材料的应力 -10- 应变关系 受压区混凝土采用的应力 - 应变关系为： （26） 式中：fc 为混凝土轴心抗压强度；0 为与 fc 对应的应变，且有0=0.002；u 为极限压应变，且有 u=0.003 8. 受拉区混凝土应力 - 应变关系为： 普通钢筋受拉和受压时都采用理想弹塑性的应力 - 应变关系. 4 算例分析 例 1 如图 3 所示的肘式框架，两端嵌固，由 William 完成其试验分析工作，由于该结构的几何非线性十分突出，同时又具有相应的试验结果，因而成为众多研究者，如 Wood， Papadrakakis， Yang， Chan， Meek 等人，检验各自所建立理论的有效性与准确性的标准算例.杆件为一
