广厦建筑结构弹塑性静力和动力分析软件GSNAP说明书--第05章.docVIP

弹塑性静力推覆分析计算原理 钢筋、钢材与混凝土的应力应变方程 图5-1 钢材料模型 图5-2 混凝土材料Saenz曲线 钢与钢筋采用同一种材料模式，即理想弹塑性增强模式，图4-1所示 混凝土的单向应力应变关系的Saenz曲线 钢筋混凝土中的混凝土模型称为延性混凝土模式。程序中使用了单向应力应变关系，见图5-3。其中上升段仍然采用Saenz曲线，下降段采用折线。一般情况下，有加密箍筋的地方应该使用延性混凝土模式，比如抗震地区的梁和柱。 图5-4 混凝土在双向应力下强度变化 (5-3) 图5-3延性混凝土材料 混凝土在双向应力下强度会发生变化，直接影响了公式中的参数。本程序选择了Kupfer强度包络线，图5-4。各标号混凝土参数取值如下： 标号 E0 (107kN/m2) σt (103kN/m2) σ0 (103kN/m2) ε0 (m) σu (103kN/m2) εu (m) 15 2.20 1.27 -10.0 -0.00127 -1.00 -0.018 20 2.55 1.54 -13.4 -0.00134 -1.34 -0.018 25 2.80 1.78 -16.7 -0.00140 -1.67 -0.018 30 3.00 2.01 -20.1 -0.00147 -2.01 -0.018 35 3.15 2.20 -23.4 -0.00153 -2.34 -0.018 40 3.25 2.39 -26.8 -0.00159 -2.68 -0.018 45 3.35 2.51 -29.6 -0.00163 -2.96 -0.018 50 3.45 2.64 -32.4 -0.00168 -3.24 -0.018 55 3.55 2.74 -35.5 -0.00173 -3.55 -0.018 60 3.60 2.85 -38.5 -0.00177 -3.85 -0.018 65 3.65 2.93 -41.5 -0.00181 -4.15 -0.018 70 3.70 2.99 -44.5 -0.00184 -4.45 -0.018 75 3.75 3.05 -47.4 -0.00188 -4.74 -0.018 80 3.80 3.11 -50.2 -0.00192 -5.02 -0.018 不考虑剪切变形的弹塑性梁柱单元 图5-5 杆单元坐标系 梁柱单元的位移模式 图5-5 梁柱满足平截面假设，所以位移是Y、Z的一次函数。位移沿X轴采用三次Hermit插值。位移模式的变化主要体现在中轴线的位移。对于一般的杆件，中轴线的位移模式： 其中： 令u，v，w分别表示单元任意一点沿x，y，z轴的位移，则 一般杆单元中都忽略了Y和Z方向的正应力和剪应力，只考虑X方向的正应力，使问题由三维降为一维问题。本程序也采用同样的假设。所以只有X方向的正应变 其中 梁柱单元的切线刚度和等效力计算 设任意时刻的切线模量则单元切线刚度矩阵为： 上式采用高斯数值积分完成。高斯积分点的密度对积分精度的影响很大，因为积分点控制着塑性区域的变化。对于实心截面，每个单元在X，Y，Z的方向采用8×8×8＝512个高斯积分点。 其他类型的截面则按实际情况进行积分。例如，工字形截面分别沿两个翼缘和腹板各采用8个高斯点积分。 切线刚度在逐步积分过程中是不断更新的。 为了减少数值积分工作量，将上式记为如下形式： 其中，为应变矩阵，为弹（塑）性矩阵 等效结点力为： 其中： 为了减少数值积分工作量，将式记为如下形式： 其中： 为按照面积积分后得到的截面内力 单元的扭转刚度和扭转等效结点力按照弹性方法得到。 钢筋混凝土梁柱是由混凝土实心元和钢筋面积两部分叠成。程序采用了同一种位移模式，各自不同的材料特性和积分过程。钢筋看作连续的板，沿着板数值积分。这种方法称作整体模式。 正应力可以用切线增量公式表示。设任意位置的切线模量和截距，几何意义参见上图，则正应力可以表示为 有限元的平衡方程采用瞬间增量平衡方程的形式，正应力需要用上式的切线增量表示。设任意