拟静力试验做完，数据出来了，但真正的麻烦才刚刚开始。提取骨架曲线——到底用最外包还是相切环？计算延性系数——Park法、通用弯矩屈服法还是能量等效法？还有累积耗能面积、等效黏滞阻尼系数、割线刚度退化、残余变形占比……每一项都要手动处理，每一张图都要反复调整。
本系统旨在为结构工程研究人员及工程技术人员提供高精度、自动化的滞回曲线分析解决方案。通过对加载-位移时程数据的深度处理与算法辨识，实现材料及构件力学性能指标的量化评估。
一、 骨架曲线提取与数据预处理
骨架曲线是反映构件在反复荷载作用下最大位移与相应承载力的包络曲线，是结构非线性分析的基础。
1.1 骨架下降段辨识逻辑
本系统集成三种先进的辨识算法，以适应不同材料与加载制度下的性能退化特征：
• 优化的最外包法 (Optimized Outer Envelope): 基于几何包络的基础上，引入最小二乘或样条插值对退化段进行平滑处理，有效剔除因数据噪声导致的局部畸变。
• 最外包法 (Outer Envelope): 提取每一级加载循环中的最大正/负位移点的极值轨迹，直接构建外包络线。
• 相切所有环法 (Tangent to All Loops): 针对复杂滞回形状，通过搜索每一个滞回环的顶点（切点），强制曲线与每一循环的峰值点相切，确保能量守恒边界的准确性。
1.2 数据处理功能
• 数值骨架曲线： 将离散的时程数据转化为具备单调递增（加载）或递减（卸载）特征的坐标序列。
• 滞回环分解： 采用零点交叉检测算法，将连续的加载序列自动拆解为独立的闭合循环（Cycle）。
• 滞回曲线平滑： 采用 Savitzky-Golay 滤波器，通过多项式拟合去除高频传感噪声，同时保持曲线的非线性特征。
• 滞回曲线代数化： 利用切比雪夫多项式或自定义解析函数，将复杂的滞回数据拟合为显式代数表达式，便于进行数学推导与计算。
二、 构件延性与屈服点评估
延性系数 μ=Δu/Δy，其中 Δu为极限位移，Δy为屈服位移。核心难点在于屈服点怎么取。 软件集成了三种主流方法：
2.1 屈服点判定方法
• Park法 (Park's Method): 通过寻找骨架曲线切线刚度退化速率突变点，或利用能量分配法确定屈服位移。
操作步骤：
• 在骨架曲线上找到峰值荷载点 A
• 确定荷载为峰值荷载75%的点 B
• 连接原点 O 与 B
• 延长交过 A 点的水平线于 C
• 过 C 作垂线交骨架曲线于 D
• D 即为屈服点
几何表达：
优缺点：计算简便，但75%系数的选取有一定经验性。
• 通用弯矩屈服法，作图法 (General Moment-Yield Method): 基于规范或构件力学模型，通过能量等效原则求得理想弹塑性模型下的屈服载荷。
操作步骤：
• 过原点 O 作骨架曲线的切线 OA
• 交过峰值点的水平线于 A
• B 为 A点在骨架曲线上的垂直投影
• 连接 OB并延长，交过峰值点的水平线于 C
• C 点的垂直投影即为屈服点
特点：纯几何作图，不引入经验系数，结果更客观。
• 能量等效法 (Energy Equivalence): 令理想弹塑性曲线下的三角形/矩形面积等于骨架曲线实测的吸收能量面积，从而确定等效屈服点。
操作步骤：
• 过原点 O 作骨架曲线的割线交骨架曲线于 B，交过峰值点水平线于 C
• 当三角形 ABC面积与弓形 OB 面积相等时，B点即为屈服点
数学表达：
优缺点：物理意义最明确（能量等效），但需要迭代求解，计算较复杂。
极限位移的确定
通常取骨架曲线上承载力下降至峰值荷载85%对应的位移：
三、 结构性能参数与动力指标
针对各滞回环的详细动力学参数分析，揭示构件在反复作用下的退化规律。
3.1 累积耗能面积
滞回环所包围的面积 SH代表一个加载循环内耗散的能量。
数学上通过数值积分：
累积耗能即所有滞回环面积之和
3.2 等效黏滞阻尼系数
反映滞回环的“饱满程度”，值越大耗能能力越强：
其中 S△OAB和 S△OCD为滞回环正负向峰值点与原点围成的三角形面积。
3.3 割线刚度
单圈加载的等效刚度：
其中 +Pi、−Δi分别为第 i圈正负向峰值荷载和峰值位移。
3.4 割线刚度退化系数
以第一圈刚度为基准，反映刚度随循环的衰减程度，βi越小，说明刚度退化越严重：
3.5 残余变形及占比
取恢复力为零时对应的位移值，即滞回环与位移轴的交点：
残余变形占比 = Δresidual/Δmax，反映构件的复位能力。
3.6 滞回环峰值承载力系数
本圈峰值荷载与骨架曲线峰值荷载的比值，反映承载力退化程度。
四、关于“区分正负向加载”的意义
构件在正向推和负向拉两个方向的受力机制可能不对称（如偏心受压柱、非对称截面构件）。软件支持正负向独立计算上述所有参数，分别输出正向和负向的延性系数、耗能、刚度退化等指标，避免 “平均化处理” 掩盖真实力学行为。
注：由于软件不能准确判断什么时候开始位移加载，因此无法计算强度退化系数，但软件输出了滞回环的峰值，可以导出数据后自行根据公式计算
