%混凝土柱尺寸為1000x2000，混凝土等級為C30
b=1000
h=2000
fc=14.3
ec=0.00164
ft=1.43
et=65*10^（-6）*ft^0.54????????%過鎮海給出的公式
at=0.312*ft^2????%與抗拉有關的參數
n=50?????%n為混凝土纖維的高度
fy=360
Es=2*10^5
es=fy/Es

??r=zeros（1100）????????%將每一個曲率對應的彎矩儲存為向量形式，選取100個曲率
??for?k=1:1:100
?????o=k/10000000
?????x1=0x2=hdf=10
??while?df>=10??????????%循環終止條件，對于這個截面，總拉力和總壓力的差小于10kN已經有足夠的精度
????x0=（x1+x2）/2x=h-x0?
??%受壓區?
?s=zeros（1（round?（x0/n）-1））??t=zeros（1（round?（x0/n）-1））??
?for?p=1:1:（round?（x0/n）-1）
?????i=n*p
?????e1=-（i-n/2-x0）*o????????%混凝土的應變
????e2=-（i-x0）*o?????????????????%鋼筋的應變
??????if?i==50?|?i==1950?
????????????steelarea=8844
??????elseif?i==150?|?i==250?|?i==350?|?i==450?|?i==550?|?i==650?|?i==750?|?i==850?|?i==950?|?i==1050?|?i==1150?|?i==1250?|?i==1350?|?i==1450?|?i==1550?|?i=