[发明专利]基于OpenSees与MATLAB的逐步增量动力分析与数据处理方法

权利要求书

1.一种基于OpenSees与MATLAB的逐步增量动力分析与数据处理方法，其特征在于，包括：

根据结构设计参数，建立OpenSees结构有限元分析模型；

地震波文件名称命名，编写地震波列表；

编写目标调幅峰值列表和地震动时间间隔列表；

编写OpenSees输出的数据文件名称命令流；

编写OpenSees加载地震波命令流；

编写在计算不收敛时自动切换算法命令流；

编写在倒塌临界状态分析终止命令流；

编写MATLAB软件读取OpenSees软件输出数据并写入到Excel表格中的命令流；

其中，“根据结构设计参数，建立OpenSees结构有限元分析模型”具体包括：

梁柱单元端部采用IMK转角弹簧单元，上部单元采用elasticBeamColumn Element与之串联；其中，IMK转角弹簧的具体参数根据截面的设计参数，由公式计算得来；由于构件是由IMK转角弹簧和elasticBeamColumn Element串联形成，所以，上部elasticBeamColumnElement需要对抗弯刚度进行修正，计算公式如下：

节点采用Joint2D Element；

采用rayleigh阻尼，利用region命令分别为弹性单元和有质量的节点赋予刚度比例系数和质量比例系数，计算过程中采用初始刚度。

2.如权利要求1所述的基于OpenSees与MATLAB的逐步增量动力分析与数据处理方法，其特征在于，“地震波文件名称命名，编写地震波列表”具体包括：

由于输入的地震波是数条甚至是数十条，因此将待分析的每条地震波进行命名，输出数据文件的名称由每条地震波名称决定，此名称亦是MATLAB后处理时寻找的文件名，为简便起见，以阿拉伯数字1,2,3......进行命名；

在OpenSees中编写地震波名称列表，记为AA，用于foreach循环。

3.如权利要求1所述的基于OpenSees与MATLAB的逐步增量动力分析与数据处理方法，其特征在于，“编写目标调幅峰值列表和地震动时间间隔列表”具体包括：

在OpenSees模型中编写目标调幅峰值列表记为BB，用于foreach循环；

在OpenSees模型中利用list命令，编写读取每条地震动时间间隔数组，记为dt_AA，用于识别每条地震波的时间间隔。

4.如权利要求1所述的基于OpenSees与MATLAB的逐步增量动力分析与数据处理方法，其特征在于，“编写OpenSees加载地震波命令流”具体包括：

首先，利用open命令打开地震动文件；利用read命令读取地震动文件中数据、lsort命令按递增顺序对地震动数据进行排序，然后，利用lindex命令获取地震动峰值，最后，llength命令获取地震动总步数。

5.如权利要求1所述的基于OpenSees与MATLAB的逐步增量动力分析与数据处理方法，其特征在于，“编写在倒塌临界状态分析终止命令流”具体包括：

已判定到达倒塌临界状态，则分析一步，使时程分析结束；利用open--r命令打开输出文件，puts命令写入0.20.2，覆盖掉分析不收敛前的记录，以此作为分析不收敛的标志，close关闭文件；

file mkdir FailRecord创建分析失败记录文件夹，open—w命令创建新的.txt文件，记录下此时分析不收敛地震波名称和目标幅值。

6.如权利要求1所述的基于OpenSees与MATLAB的逐步增量动力分析与数据处理方法，其特征在于，“编写MATLAB软件读取OpenSees软件输出数据并写入到Excel表格中的命令流”具体包括：

利用MATLAB中的sprintf命令读取输出的文件数据，max函数寻找损伤指标最大值，xlswrite函数，将地震动名称和目标峰值及损伤指标最大值写入Excel表格中的指定位置。