如何绘制单轴循环加载（0→(ε=0.02↔ε=-0.005)×20次）应力应变图

一、循环加载应力应变可视化方法 #

针对0→(ε=0.02拉伸 ↔ ε=-0.005压缩)×20次的加载路径，按以下步骤实现可视化：

• 坐标轴设置：横轴为应变(ε)，建议范围设为[-0.006, 0.021]（覆盖压缩极值并留余量）；纵轴为应力(σ)，根据550级钢筋接头的力学性能设置合理范围（比如0到600MPa）。
• 曲线绘制规则：
  1. 按加载顺序逐次绘制每一轮循环的应力应变轨迹，可通过颜色渐变或线型区分循环次数（比如前几次用实线，后续用虚线，或用从浅到深的蓝色系区分）。
  2. 标注关键节点：在图中明确标出起始点(0,0)、拉伸极值点(0.02, σ_t)、压缩极值点(-0.005, σ_c)，并在第5、10、15、20次循环旁标注循环次数，方便观察循环衰减规律。
• 代码实现示例（Python Matplotlib）：

import matplotlib.pyplot as plt
import numpy as np

# 模拟550级钢筋接头的循环应力应变数据（需替换为实际试验数据）
total_cycles = 20
# 初始加载段（0→2%→-0.5%）
init_strain = [0, 0.02, -0.005]
init_stress = [0, 520, -130]  # 模拟初始应力值

plt.figure(figsize=(8, 6))
# 绘制初始循环
plt.plot(init_strain, init_stress, color='#1f77b4', label='Cycle 1')
# 绘制后续19次循环
for cycle in range(2, total_cycles+1):
    # 模拟循环应力衰减（实际需用真实试验数据）
    current_strain = [init_strain[-1], 0.02, -0.005]
    current_stress = [init_stress[-1], 520 - (cycle-1)*3, -130 + (cycle-1)*1.5]
    line_color = plt.cm.Blues(0.2 + (cycle-1)*0.04)
    plt.plot(current_strain, current_stress, color=line_color, 
             label=f'Cycle {cycle}' if cycle in [5,10,15,20] else "")

plt.xlabel('Strain (ε)')
plt.ylabel('Stress (MPa)')
plt.title('Stress-Strain Hysteresis Curve for 550-grade Rebar Mechanical Splices')
plt.xlim(-0.006, 0.021)
plt.grid(True, linestyle='--', alpha=0.7)
plt.legend(loc='upper right')
plt.tight_layout()
plt.show()

二、550级钢筋机械接头试验要求 #

• 加载循环规范：严格执行拉伸2%应变与压缩0.5%应变之间的非弹性循环，共完成20次循环，加载起始状态为无应力、无应变的自然状态。
• 应变测量要求：必须在紧邻但不包含机械拼接装置的钢筋段上布置应变采集元件（如应变片），测量点需避开接头本体及接头附近的应力集中区域，确保采集到的是钢筋母材的真实应变数据，避免接头变形干扰测量结果。
• 数据记录要求：每一次循环的拉伸、压缩极值点的应力应变数据需完整记录，用于绘制滞回曲线并分析循环耐久性。

内容的提问来源于stack exchange，提问作者Nurul Alma