弹塑性时程分析程序设计与实例-kinco hmiware使用手册

8.3 弹塑性时程分析程序设计与实例 实际上，结构在较大的地震过程中往往表现出非线性，即结构的刚度和阻尼是 随时间变化的参数。对于钢结构一般采用双线性模型，对于钢筋混凝土结构一般采 用三线型刚度退化模型。本节的程序设计采用复杂的三线型刚度退化模型。现按 Wilson-θ法编制结构的动力反应分析程序，由于该程序非常复杂，为使编程思路更 加清晰，将相应程序框图绘出，如图 8.8 所示。开始输入结构参数和地震参数，调 用子程序计算结构的割线刚度，即 ，并赋正副 大割线刚度卸载力（初值为正负 屈服力）。给程序各变量赋初值，从地震波的第一时刻起，依据结构的恢复力模型判 断结构各层的刚度，同时判断该地震时刻是否需要进行拐点处理，如果需要进行拐 点处理运算，将所计算的地震反应代替刚刚所求的地震反应数值（即上一时刻），然 后进行该时刻的地震反应求解运算，如果不需要进行拐点处理运算，则直接进行地 震反应求解运算。地震反应的求解过程按照 8.1 节的各公式进行，得出该时刻的地 震反应后，形成位移、速度、加速度和地震剪力的数值矩阵，如此反复循环，直至 地震的全过程计算完成。然后调用绘图程序，绘制地震反应时程曲线。 4k 例 8.2 某三层钢筋混凝土结构，结构的各层特性参数为：各层质量 =[1.690，1.581，1.412] ×10m 4kg，各层刚度 =[1.512，2.012，2.012] ×10k 7N/m，各 层第一刚度折减系数为 0.4，第二刚度折减系数为 0.1，各层的开裂位移为 ＝[6.3， 4.9，4.2]mm，各层的屈服位移为 ＝[21.8，18.9，17.2]mm，地震波采用 400gal El Centro波，采样周期为 0.02 秒。地震波采用El Centro波，结构的恢复力模型采用三 线型刚度退化模型。 cx yx 主程序 %%%%%%%%% MAIN PROGRAM %%%%%%%%% (made in July-October, 2000 by Xu Zhao-dong) %(revised in January 6--16, 2002)% % this program is three-poly stiffness degeneration model % this program is the reinforced concrete structure % structure parameters m0=[1.690 1.581 1.412]*1e+4; k0=[1.512 2.012 2.012]*1e+7; 12