[发明专利]RC圆柱截面曲率极限状态快速分析方法

权利要求书

1.一种RC圆柱截面曲率极限状态快速分析方法，其特征在于，包括：

根据截面设计参数，建立OpenSees截面曲率有限元分析模型；

定义OpenSees软件输出的数据文件名称；

编写OpenSees计算不收敛的处理方案；

通过MATLAB软件，读取OpenSees软件输出的数据文件，并把数据写入到Excel表格中的指定位置。

2.如权利要求1所述的RC圆柱截面曲率极限状态快速分析方法，其特征在于，“根据截面设计参数，建立OpenSees截面曲率有限元分析模型；”具体包括：

将截面设计参数，轴向荷载、箍筋间距、纵筋直径形成各自的参数列表；

混凝土材料分为保护层混凝土和箍筋约束后的核心混凝土两种，有限元模型中采用Concrete02单轴材料；

箍筋约束混凝土本构关系采用Mander模型，在OpenSees命令流中编写Mander模型方程计算式，根据输入的截面设计参数，自动计算并输出箍筋约束后的混凝土本构参数；

有限元模型中钢筋材料采用考虑低周疲劳损伤的ReinforcingSteel单轴材料，疲劳参数采用OpenSees官网推荐参数；

zeroLengthSection单元建立截面曲率分析模型。

3.如权利要求1所述的RC圆柱截面曲率极限状态快速分析方法，其特征在于，“定义OpenSees软件输出的数据文件名称；”具体包括：

截面极限状态由边缘混凝土受压应变或纵筋受拉应变决定，OpenSees计算结果输出文件包括，截面曲率、边缘混凝土受压应变和纵筋受拉应变输出文件；

计算结果输出文件由输入的截面设计参数命名，即轴压比-箍筋率-纵筋率。

4.如权利要求1所述的RC圆柱截面曲率极限状态快速分析方法，其特征在于，“编写OpenSees计算不收敛的处理方案；”具体包括：

计算不收敛时自动切换积分方法：(KrylovNewton)-(Newton)-(Newton-initial)-(Broyden)-(NewtonLineSearch)-(ModifiedNewton)-(Linear)；

如切换积分方法仍不收敛时，增大判断收敛的容差：(1.0e-8)-(1.0e-4)；

在上述步骤改变的情况下仍不收敛时，改变判断收敛的条件，由位移增量判断收敛，改成能量增量判断收敛：(testNormDispIncr)-(test EnergyIncr)；在执行上述方案的情况下，仍不收敛，退出计算程序，检查模型。

5.如权利要求1所述的RC圆柱截面曲率极限状态快速分析方法，其特征在于，“通过MATLAB软件，读取OpenSees软件输出的数据文件，并把数据写入到Excel表格中的指定位置。”具体包括：

建立Excel表格，确定每列将要写入的数据内容；

建立MATLAB软件脚本，sprintf命令读取OpenSees计算结果输出文件中的数据，放入计算机内存中，index命令寻找给定的截面曲率极限状态时的混凝土受压应变或钢筋受拉应变的位置，即行数，然后通过xlswrite命令将所在行数位置的曲率，写入到Excel表格中的指定位置。

6.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1到5任一项所述方法的步骤。

7.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现权利要求1到5任一项所述方法的步骤。

8.一种处理器，其特征在于，所述处理器用于运行程序，其中，所述程序运行时执行权利要求1到5任一项所述的方法。