MATLAB-能量法计算骨架曲线性能点

以下内容均来自网络，如若侵权，请联系作者删除！！！

当我们对钢筋混凝土结构力学分析时，往往需要对骨架曲线进行分析，获得相应的抗震性能指标点（屈服位移、屈服荷载、峰值位移、峰值荷载、极限位移、极限荷载，延性系数）。常用的方法有几何作图法、Park法以及能量法。能量法的核心思想是将骨架曲线等效为理想弹塑性曲线进行计算，示意图如下图所示：

（侵权联系作者删除！！！）

折线OBC为理想弹塑性荷载-位移曲线，与实际荷载-位移曲线相交于A点，当曲边三角形的面积SABC = SOA时，过B点作垂直于x轴的直线，与实际荷载-位移曲线相交于D点，D点则为等效屈服点，与D点对于的实际位移则为屈服位移，与D点对应的实际荷载则为屈服荷载。极限点对应于荷载下降至0.85倍峰值荷载。

MATLAB程序基于“能量法”可获得。骨架曲线Skeletion的等效屈服位移Dy、等效屈服荷载Fy、峰值位移Dp、峰值荷载Fp、极限位移Du、极限荷载Fu（0.85×Fp）、延性系数（Du/Dy）Lamda以及计算误差Err。计算误差Err可通过增大插值点数Divided降低。

运行程序后会产生计算结果，便于查看计算是否正确。

结果在工作区的变量中显示。

需注意：增大插值点数Divided越大越好。当大到一定程度，计算结果稳定。自行调试；骨架曲线的极限荷载应小于等于0.85峰值荷载，否则会报错。