RC圆柱截面曲率极限状态快速分析方法技术

本发明专利技术公开了一种RC圆柱截面曲率极限状态快速分析方法。本发明专利技术一种RC圆柱截面曲率极限状态快速分析方法，包括：根据截面设计参数，建立OpenSees截面曲率有限元分析模型；定义OpenSees软件输出的数据文件名称；编写OpenSees计算不收敛的处理方案；通过MATLAB软件，读取OpenSees软件输出的数据文件，并把数据写入到Excel表格中的指定位置。本发明专利技术的有益效果：本发明专利技术提出了一种RC圆柱截面曲率极限状态快速分析方法，克服了在截面设计参数，如轴压比、箍筋率和纵筋率改变的情况下，需手动更改设计参数，无法批量进行曲率分析的难题，极大提高了计算效率，使科研人员把更多的精力花在参数分析上。

【技术实现步骤摘要】
RC圆柱截面曲率极限状态快速分析方法
本专利技术属于建筑或桥梁结构工程抗震分析领域，具体涉及一种土木工程专业中的RC圆柱截面曲率极限状态快速分析方法。
技术介绍
柱子是建筑和桥梁结构承受竖向荷载和抗侧向力的主要构件，在发生地震时，防止结构倒塌起着关键的作用。在基于性能的结构抗震设计方法中，对构件的变形进行量化是确定结构系统性能水平的依据，用材料应变极限状态和截面曲率极限状态作为构件弹塑性变形的量化指标得到了广泛认可。材料应变的极限状态由于受力复杂，难以检测，而截面曲率相对稳定，因此材料应变极限状态对应的截面曲率得到了更为普遍的应用。截面曲率极限状态的大小受截面设计参数的影响，故明确轴压比、箍筋率和纵筋率对截面曲率极限状态的作用，是准确量化构件变形量的前提。传统技术存在以下技术问题：目前计算截面曲率的方法包括手算和电算两种。手算：科研工作者在大量试验数据的基础上，通过线性回归总结出计算截面屈服曲率的经验公式，截面屈服曲率表示构件处于弹性状态，而地震作用下的构件大多进入了塑形发展阶段，截面屈服曲率已无法量化构件的变形，此外计算出的截面屈服曲率与材料应变也无法形成连续的关系。电算即分析软件，如：OpenSees、XTRCT和USC_RC分析软件，这三种软件可以根据使用者输入的轴压比、箍筋率和纵筋率等参数进行截面曲率分析，并可以实时输出材料应变和与之对应的截面曲率。但如果要研究某一变量对截面曲率极限状态的影响，需要手工记录上一次的分析数据，然后重新输入新的截面设计参数，再进行分析，以此类推。在此指出，OpenSees软件不可视化，分析数据不能自行呈现图表，进一步限制其应用范围。如果研究的变量多，每次变动得幅度小，像这样机械化操作，效率低下，消耗大量的时间和精力，根本无法大规模的进行分析，无法总结出每一变量对截面曲率极限状态的影响，更无法对实际的工程作出优化。
技术实现思路
本专利技术要解决的技术问题是提供一种RC圆柱截面曲率极限状态快速分析方法。为了解决上述技术问题，本专利技术提供了一种RC圆柱截面曲率极限状态快速分析方法，包括：根据截面设计参数，建立OpenSees截面曲率有限元分析模型；定义OpenSees软件输出的数据文件名称；编写OpenSees计算不收敛的处理方案；通过MATLAB软件，读取OpenSees软件输出的数据文件，并把数据写入到Excel表格中的指定位置。在其中一个实施例中，“根据截面设计参数，建立OpenSees截面曲率有限元分析模型；”具体包括：将截面设计参数，轴向荷载、箍筋间距、纵筋直径形成各自的参数列表；混凝土材料分为保护层混凝土和箍筋约束后的核心混凝土两种，有限元模型中采用Concrete02单轴材料；箍筋约束混凝土本构关系采用Mander模型，在OpenSees命令流中编写Mander模型方程计算式，根据输入的截面设计参数，自动计算并输出箍筋约束后的混凝土本构参数；有限元模型中钢筋材料采用考虑低周疲劳损伤的ReinforcingSteel单轴材料，疲劳参数采用OpenSees官网推荐参数；zeroLengthSection单元建立截面曲率分析模型。在其中一个实施例中，“定义OpenSees软件输出的数据文件名称；”具体包括：截面极限状态由边缘混凝土受压应变或纵筋受拉应变决定，OpenSees计算结果输出文件包括，截面曲率、边缘混凝土受压应变和纵筋受拉应变输出文件；计算结果输出文件由输入的截面设计参数命名，即轴压比-箍筋率-纵筋率。在其中一个实施例中，“编写OpenSees计算不收敛的处理方案；”具体包括：“编写OpenSees计算不收敛的处理方案；”具体包括：计算不收敛时自动切换积分方法：(KrylovNewton)-(Newton)-(Newton-initial)-(Broyden)-(NewtonLineSearch)-(ModifiedNewton)-(Linear)；如切换积分方法仍不收敛时，增大判断收敛的容差：(1.0e-8)-(1.0e-4)；在上述步骤改变的情况下仍不收敛时，改变判断收敛的条件，由位移增量判断收敛，改成能量增量判断收敛：(testNormDispIncr)-(testEnergyIncr)；在执行上述方案的情况下，仍不收敛，退出计算程序，检查模型。在其中一个实施例中，“通过MATLAB软件，读取OpenSees软件输出的数据文件，并把数据写入到Excel表格中的指定位置。”具体包括：建立Excel表格，确定每列将要写入的数据内容；建立MATLAB软件脚本，sprintf命令读取OpenSees计算结果输出文件中的数据，放入计算机内存中，index命令寻找给定的截面曲率极限状态时的混凝土受压应变或钢筋受拉应变的位置，即行数，然后通过xlswrite命令将所在行数位置的曲率，写入到Excel表格中的指定位置。一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现任一项所述方法的步骤。一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现任一项所述方法的步骤。一种处理器，所述处理器用于运行程序，其中，所述程序运行时执行任一项所述的方法。本专利技术的有益效果：本专利技术提出了一种RC圆柱截面曲率极限状态快速分析方法，克服了在截面设计参数，如轴压比、箍筋率和纵筋率改变的情况下，需手动更改设计参数，无法批量进行曲率分析的难题，极大提高了计算效率，使科研人员把更多的精力花在参数分析上。建立了MATLAB软件数据后处理脚本，可自行读取OpenSees计算结果输出文件中的数据，自动寻找给定的截面曲率极限状态时的混凝土受压应变或钢筋受拉应变的位置，即行数，然后将所在行数位置的曲率，写入到Excel表格中的指定位置。在此充分实现了批量计算，无需人工干预的优点；另外随着软件、测量仪器精度的增加，结构性能水平所对应的材料应变极限状态已有变化，本方法可以灵活定义材料应变极限状态，提取与之对应的截面曲率极限状态，可对既有的截面曲率极限状态公式做进一步的修正提供极大的便利。附图说明图1RC圆柱截面曲率极限状态快速分析方法的流程图。图2OpenSees截面曲率有限元分析模型。图3保护层未约束混凝土与Mander箍筋约束混凝土本构关系。图4材料骨架曲线：(a)Concrete02单轴材料；(b)ReinforcingSteel单轴材料。图5MATLAB读写输出数据流程图。图6混凝土受压应变0.004时的截面曲率极限状态分析图：(a)纵筋率为变量；(b)轴压比为变量。具体实施方式下面结合附图和具体实施例对本专利技术作进一步说明，以使本领域的技术人员可以更好地理解本专利技术并能予以实施，但所举实施例不作为对本专利技术的限定。一种RC圆柱截面曲率极限状态快速分析方法，包括：步骤1，根据截面设计参数，建立OpenSees截面曲率有限元分析模型；步骤2，定义OpenSees软件输出的数据文件名称；步骤3，编写OpenSees计算不收敛的处理方案；步骤4，通过MATLAB软件，读取OpenSees软件输出的数据文件，并把数据写入到Excel表格中的指定位置；本专利技术的有益效果在于：本专利技术提出了一种RC圆柱截面曲率极限状态快速分析方法，克服了在截本文档来自技高网...

【技术保护点】
1.一种RC圆柱截面曲率极限状态快速分析方法，其特征在于，包括：根据截面设计参数，建立OpenSees截面曲率有限元分析模型；定义OpenSees软件输出的数据文件名称；编写OpenSees计算不收敛的处理方案；通过MATLAB软件，读取OpenSees软件输出的数据文件，并把数据写入到Excel表格中的指定位置。

【技术特征摘要】
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计算不收敛的处理方案；”具体...

【专利技术属性】
技术研发人员：刘赛，顾冬生，熊能，
申请(专利权)人：江南大学，
类型：发明
国别省市：江苏,32