资源简介
这个小程序是利用matlab进行运算的，是通过纤维截面来求解梁截面的弯矩-曲率的关系，混凝土和钢筋的本构均采用规范取值。
代码片段和文件信息
%混凝土柱尺寸为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=
评论
共有 条评论