钢筋混凝土柱的Mphi曲线matlab代码

以下是一个简单的钢筋混凝土柱的M-phi曲线的Matlab代码：

% 定义参数
fc = 30; % 混凝土抗压强度
fy = 400; % 钢筋抗拉强度
b = 250; % 截面宽度
h = 500; % 截面高度
As = 4*25^2; % 钢筋面积
d = 450; % 钢筋所在位置到截面底部的距离
eps_c = -0.0035; % 混凝土应变
eps_s = -0.0025; % 钢筋应变

% 定义函数
phi = @(eps_s) eps_s*fy*As - (fc*b*(h-d) + As*fy*d)*eps_c;

% 绘制M-phi曲线
eps_s = linspace(-0.005,0.005,1000);
M = phi(eps_s);
plot(M, eps_s)
xlabel('弯矩（kN·m）')
ylabel('钢筋应变')
title('M-phi曲线')

这个代码中假设柱子是矩形截面，有4根直径为25mm的钢筋，钢筋距离截面底部450mm，混凝土抗压强度为30MPa，钢筋抗拉强度为400MPa。通过计算得到M-phi曲线，其中phi是一个关于钢筋应变的函数。最终绘制出M-phi曲线图。