marlab绘制混凝土梁弯矩-曲率曲线

本文利用MATLAB实现了绘制混凝土梁弯矩与曲率曲线关系。

clear 
clc 
format long 
h=1000;b=600;h0=h-40;                          %长度单位均为 mm 
fc=26.8;ft=2.39;Es=200000;Ec=32500;
s=0.00001:0.00004:0.008;                      %自变量为受压区边缘混凝土压应变，0-0.008，步长 0.00004
n=100;q=50;                                    %混凝土划分为n 层工字钢q层
m=zeros（1200）;                                %定义数组
for  i=1:200 
    x（1）=0;x（2）=h0;x（3）=0.5*h0;    
for  j=3:200                               %如果迭代次数超过 200 次，即认为迭代发散，终止迭代 
    f1=0;f2=0;f3=0;f4=0;r1=0;r2=0;r3=0;r4=0;
for  aa=1:6 
        Xbar（aa）=40+（aa-1）*184;
        As（1）=3.14*25*25*6/4;
        As（2）=3.14*25*25*2/4;
        As（3）=3.14*25*25*2/4;
        As（4）=3.14*25*25*2/4;
        As（5）=3.14*25*25*2/4;
        As（6）=3.14*25*25*6/4;
        strain（aa）=（x（j）-Xbar（aa））*s（i）/x（j）;           %钢筋应变 
    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         131  2018-11-05 10:29  fun1.m
     文件         236  2018-11-05 16:47  fun2.m
     文件         143  2018-10-30 17:23  fun3.m
     文件        3204  2018-11-03 22:00  fun.m