基于MCFT理论的割线刚度法介绍.docVIP

基于MCFT理论的割线刚度法介绍 F. J. Vecchio于1990年在ACI发表的论文中采用了基于MCFT理论的割线刚度法来进行钢筋混凝土的有限元分析。通过比较割线刚度法与基于切线刚度的Newton-Raphson法，可以得出以下几点结论： 1. 在结构的应力-应变曲线进入到下降段后，对于切线刚度法，刚度会出现负值，从而给计算处理带来难度。而割线刚度法则没有这种负刚度问题，但是割线刚度法会在卸载曲线穿越坐标轴时出现数值问题，这是由于此时割线刚度极小造成的。 2. 以往采用切线刚度法的程序，在更新应变时，对于一个单元采用了4个高斯积分点，每次迭代得到了一个新的结构位移，就需要在4个积分点上分别求出它们的应变，再以此计算出4个新的割线刚度，用以判断收敛与否。而F. J. Vecchio采用的割线刚度法只涉及了1个积分点，即中心积分点。 根据附表1可知，单元应变和位移的关系为： 其中： 令，，则可以得到割线刚度法中心积分点的应变计算式为： 每次循环迭代出的新位移，仅需要换算成中心积分点上的应变，通过一次计算便可以得出割线刚度用于收敛判断。相比于切线刚度法，割线刚度法简化了计算过程。 3. 切线刚度法的过程，需要通过外荷载节点力和结构恢复力求解出一个残余力，再以此计算出增量位移来判断是否满足收敛条件，若不收敛，便更新结构位移并重新计算结构的切线刚度矩阵。对于割线刚度法，并不需要求解残余力，而是直接利用外荷载节点力和结构割线刚度来求解出一个全量位移，之后通过计算出的应变来重新计算结构的割线刚度矩阵。二者相比起来，割线刚度法不需要启用一个新的数组来储存残余力值，降低了对程序的存储空间要求并提高了程序的运算速度。 从上面三点可以看出，切线刚度法和割线刚度法各有特点和不足，在某些方面，割线刚度法还体现出了一定的优势： 割线刚度法相比于切线刚度法的收敛准则更为简单。 割线刚度法相比于切线刚度法的有限元模型更为简单。 割线刚度法的稳定性较好。 所以对割线刚度法进行研究佐证是十分有必要的。 F. J. Vecchio在其论文中所述的割线刚度法的具体过程如下： 第1步：输入结构和材料的属性数据 第2步：计算外荷载节点力 第3步：计算材料的刚度矩阵： 1) 确定初始应变：，，； 2) 计算主方向应变，，，，以及主拉应变角； 3) 根据混凝土的本构关系，确定混凝土主拉和主压应力， 4) 根据钢筋的本构关系，确定钢筋的主拉和主压应力，； 5) 计算混凝土和钢筋的割线刚度元素值： (3-1) 第4步：计算混凝土和钢筋的材料割线刚度矩阵 1) 计算混凝土主方向上的材料割线刚度矩阵： (3-2) 2) 计算转换矩阵 3) 计算整体坐标系（X-Y坐标）下的混凝土割线刚度矩阵 (3-3) 4) 计算整体坐标系(X-Y坐标系)下的钢筋割线刚度矩阵： (3-4) 第5步：计算组合材料刚度： (3-5) 第6步：计算单元刚度矩阵： (3-6) 第7步：计算单元预应变产生的位移， 第8步：计算单元预应变产生的节点力， 第9步：计算单元节点力： 第10步：计算结构总刚矩阵： (3-7) 图 3.5 单元尺寸 第11步：计算结构总刚矩阵的逆矩阵： 第12步：计算节点位移 Fig3.5 Element dimension (3-8) 第13步：计算单元应变（取单元形心处的应变为平均应变）： (3-9) (3-10) (3-11) 注：，为相应节点的水平位移和竖向位移。 第14步：检查单元应力： 第15步：重新计算单元割线刚度值： 1) 根据求解出的位移，重新确定主方向上的应变值，，，；