混凝土结构M-N关系图计算机数值分析C++编程
钢筋混凝土全截面M-φ曲线matlab程序原代码

Xyy(2)=850;
Gas(1)=400*12;
Gas(2)=400*12; %工字钢翼缘面积
gstrain(dd)=(x(j)-Xyy(dd))*s(i)/x(j); %工字钢翼缘应变
cstrain(kk)=(x(j)-Xc(kk))*s(i)/x(j); %混凝土第k层应变
cstress(kk)=fun2(cstrain(kk)); %混凝土第 k 层应力。fun2 混凝土受压时的应力应变关系
end
f=f1+f2+f3+f4
if f<=1&&f>=0
v(i)=s(i)/x(j); %曲率Φ
m(i)=r1+r2+r3+r4; %整个截面的力对中性轴取矩 M
for j=3: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;
ystress(ee)=fun3(ystrain(ee)); %工字钢腹板应力。fun3 为自定义函数,钢筋的应力应变关系
f4=f4+ ystress(ee)*ghs*12; %工字钢腹板合力
r4=r4+ ystress(ee)*ghs*12*(x(j)-Xfb(ee)); %工字钢腹板合力对中性轴的力矩
f1=f1+ stress(aa)*As(aa); %钢筋合力
r1=r1+ stress(aa)*As(aa)*(x(j)-Xbar(aa)); %钢筋合力对中性轴的力矩
end
混凝土及钢筋混凝土工程量计算资料

4.现浇混凝土墙按设计图示尺寸以体积计算。墙垛及突出
墙面部分并入墙体体积计算内。 •有梁的计至梁底,与墙同厚的梁,其工程量并入墙计算, 没有梁的计至板面。 •有地下室的从地下室底板面计起,没有地下室的从基础 面计起,楼层从板面计起。
91
5.现浇混凝土板按设计图示尺寸以体积计算。
•有梁板(包括主、次梁与板)按梁、板体积之和计算 •无梁板按板和柱帽体积之和计算 •平板按板体积计算。 •各类板伸入墙内的板头并入板体积内计算 •薄壳板的肋、基梁并入薄壳体积内计算。 •板伸出墙外150mm以内的仍按板计算
3
4
5-1
5 局部放大
6
第一节 现浇混凝土
一、定额说明 1 现浇混凝土工程分为混凝土制作和现浇混凝土浇
捣两部分。现浇混凝土浇捣工程量,按混凝土结构 体积计算。混凝土制作工程量,按现浇混凝土浇捣 子目中的混凝土含量(已包括损耗率)计算。 2 毛石基础子目,按毛石占毛石混凝土体积的20%考 虑。如设计不同,可以换算。 3 小型构件指每件体积在0.05m3以内的未列出子目 的构件。
100
92
补:现浇混凝土板项目的划分:
(1) 有梁板,指带肋的混凝土板,梁板一起整体浇注。
肋形板:板中梁为单向;梁中距>1m;
有梁板 密肋板:板中梁为单向;梁中距≤1m;
井字板:板中梁为双向交叉呈“井”字形布置。
左图: 梁中距 L>1m,为肋形板 梁中距 L≤1m,为密肋板
LL L
注意:圈梁、以柱为支座的连系梁、框架 梁,不得并入板内按有梁板计算。
22
5、桩基础。多数用于高层建筑或土质不好的情况下 ,由 若干桩来支承一个平台,然后由这个平台托住整个建筑物 ,叫做桩承台
23
(完整版)混凝土结构施工图平面整体表示方法制图规则和构造详图

平法讲座########公司2014年7月目录第一部分基本概念 (1)第二部分柱平面表示法 (3)一、注写 (3)二、抗震KZ钢筋构造 (5)三、非抗震KZ钢筋构造 (8)四、柱钢筋制作安装质量通病 (9)第三部分梁平面表示法 (10)一、注写 (10)二、抗震KL、WKL纵向钢筋构造 (12)三、梁配筋构造 (15)四、非抗震KL同抗震KL的异同点 (16)五、梁钢筋制作安装质量通病 (16)第四部分剪力墙平面表示法 (17)一、墙柱 (17)二、墙身 (20)三、墙梁 (22)四、洞口平面表示法 (24)混凝土结构施工图平面整体表示方法制图规则和构造详图(03G101-1)基本概念第一部分基本概念1、接头百分率(1)排扎接头连接区段长度为1.3L,区段内接头面积百分率梁板墙不宜>25%,柱不宜>50%,确有必要增加时梁不应>50%,其他构件可放宽。
(2)机械连接及焊接接头连接区段为35d(d为较大直径)且不小于500mm,凡接头中心点位于该区段内的接头均属于内一区段,区段内接头百分率,受拉区不宜>50%,接头不宜设在抗震框架梁、柱端加密区。
直接承受动力荷载的不宜采用焊接接头。
2、HPB235钢筋:180°弯钩,内直径≥2.5d,平直段长度不应小于3d,箍筋弯钩内直径不小于受力钢筋,一般结构弯90°,平直部分不宜小于5d,抗震结构应为135°,平直部分不应小于10d且不应小于75。
3、抗震等级钢筋混凝土房屋根据烈度、设防类型、结构类型和房屋高度,场地类别采用不同的抗震等级并应符合相应计算和构造措施要求。
4、锚固长度L 1、L1E纵向受拉钢筋绑扎搭接长度l1E、l1纵向受拉钢筋搭接长度修正数注:当不同直径的钢筋搭接时,其l1E与l1值按较小的直径计算;在任何情况下l1不得小于300mm;式中ζ为搭接长度修正数。
基本概念混凝土结构施工图平面整体表示方法制图规则和构造详图(03G101-1)L a 、LaE纵向受拉筋抗震锚固长度LaE受拉钢筋的最水锚固长度La注:1、四级抗震等级L aE=L a;2、当弯锚时,有些部位的锚固长度为≥0.4L aE+15d,见各类构件的标准构造祥图。
混凝土结构设计图表

混凝土结构设计图表目录1.基本设计规定 (1)1.1建筑结构的安全等级 (1)1.2受弯构件的挠度限值 (1)1.3结构构件的裂缝控制等级及最大裂缝宽度限值 (2)1.4耐久性规定 (2)1.4.1混凝土结构的环境类别 (2)1.4.2 结构混凝土耐久性的基本要求 (3)2 材料 (3)2.1混凝土 (3)2.2钢筋 (4)2.2.1普通钢筋强度标准值 (4)2.2.2预应力钢筋强度标准值 (5)2.2.3普通钢筋强度设计值 (5)2.2.4预应力钢筋强度设计值 (6)2.2.5钢筋弹性摸量 (6)2.2.6普通钢筋疲劳应力幅限值 (7)2.2.7预应力钢筋疲劳应力幅限值 (8)2.2.8钢筋公称截面面积、计算截面面积及理论重量 (8)2.2.9钢绞线公称截面面积、计算截面面积及理论重量 (9)2.2.10钢丝公称截面面积、计算截面面积及理论重量 (9)3 预应力混凝土结构构件计算要求 (9)3.1一般规定 (9)3.1.1 拉控制应力限值 (9)3.1.2预应力钢筋及非预应力钢筋合力位置图 (10)3.1.3预应力混凝土施工阶段验算图 (10)3.2预应力损失值计算 (11)3.2.1预应力损失值表 (11)3.2.2锚具变形与钢筋缩值 (12)3.2.3摩擦系数 (12)3.2.4预应力损失计算图 (12)3.2.5各阶段预应力损失值的组合表 (12)4 承载能力极限状态计算 (13)4.1正截面受弯承载能力计算 (13)4.1.1工形截面受弯构件受压区高度位置图 (13)4.1.2 T形、及工字形及倒L形截面受弯构件翼缘 (14)4.2正截面受压承载能力计算 (14)4.2.1 钢筋混凝土轴心受压构件的稳定系数 (15)4.2.2 刚性屋盖单层房屋排架柱、露天吊车柱和栈桥柱的计算长度 (15)4.2.3 刚性屋盖单层房屋排架柱、露天吊车柱和栈桥柱的计算长度 (16)5 正常使用极限状态验算 (16)5.1裂缝控制验算 (16)5.1.1 构件受力特征系数 (16)5.1.2 钢筋的相对粘结特性系数 (16)5.2受弯构件挠度验算 (17)5.2.1 截面抵抗矩塑性影响系数基本值 (17)6 构造规定 (17)6.1钢筋混凝土结构伸缩缝最大间距 (17)6.2 纵向受力钢筋的混凝土保护层最小厚度 (18)6.3 钢筋的外形系数 (18)6.4 纵向受拉钢筋搭接长度修正系数 (18)6.5 纵向受力钢筋的最小配筋百分率 (19)7 结构构件的基本规定 (19)7.1现浇混凝土板的最小厚度 (19)7.2 梁中箍筋最大间距 (20)8 混凝土结构构件抗震设计 (20)8.1一般规定 (20)8.1.1现浇混凝土房屋适用的最大高度 (20)8.1.2混凝土结构的抗震等级 (21)8.1.3承载力抗震调整系数 (21)8.2 框架梁 (22)8.2.1框架梁纵向受拉钢筋的最小配筋百分率 (22)8.2.2框架梁梁端加密区的构造要求 (22)8.3 框架柱及框支柱 (22)8.3.1柱全部纵向受力钢筋最小配筋百分率 (22)8.3.2柱全部纵向受力钢筋最小配筋百分率 (23)8.3.3柱全部纵向受力钢筋最小配筋百分率 (23)8.3.4柱全部纵向受力钢筋最小配筋百分率 (24)8.4铰接排架柱箍筋加密区的箍筋最小直径 (24)8.5剪力墙 (25)8.5.1 墙肢轴压比限值 (25)8.5.2 剪力墙设置构造边缘构件的最大轴压比 (25)8.5.3 约束边缘构件沿墙肢的长度及其配箍特征值 (25)8.5.4 约束边缘构件沿墙肢的长度及其配箍特征值 (26)9 混凝土多轴强度和本构关系 (27)9.1单轴应力应变曲线 (27)9.2多轴应力应变曲线 (28)10 其他 (29)注:《混凝土结构设计规》GB50010-2002在本文中简称规。
钢筋混凝土工程量计算程序公式

钢筋混凝土工程量计算程序公式:1、圆孔板工程量计算程序公式:预算定额中关于圆孔板的定额项目分别列有圆孔板预制(包括钢筋加工和绑扎)、圆孔板安装。
若在预制构件厂制作或购买,尚需计算圆孔板构件的蒸汽养护费、从预制厂至工地的运输费。
因此一般需要计算圆孔板的制作、蒸汽养护、运输、安装四项费用,即计算四项工程量。
按圆孔板计算出来的为安装工程量,安装工程量再增加1%的安装损耗是制作、养护、运输的工程量,钢筋数量也要计算出来。
工程量计算程序及方法如下:首先根据机构平面图统计出每中规格的圆孔板的数量,然后查《圆孔板标准图》,从经济指标栏中查得每块板的混凝土和钢筋数量,再计算出总的制作、安装工程量和钢筋数量。
2、现浇板缝工程量计算程序公式:首先统计板缝数量列于表中(表如下例)板缝数量统计定额规定,板缝宽20~200mm,套板缝定额;板缝宽200mm以上,套现浇板定额。
板缝计算公式:a、混凝土工程量:开间宽度X板缝宽度X圆孔板厚度(0.13)X板缝条数=?立方b、钢筋工程量:主筋(①XX):[开间宽度+两端弯钩长度]X钢筋根数X板缝条数X型钢理论重量=?kg箍筋(①XX):缝宽X开间宽度三箍筋间距(@)X型钢理论重量X板缝条数=?kg附:钢筋均未计算施工损耗,待工程量算完后,再统一按钢筋施工损耗表计算。
3、现浇板工程量计算程序公式:定额规定,现浇板分平板(板四周支承在砖墙或圈梁上)和有梁板(板和现浇梁浇制在一起),应分别计算,并分别套平板、后梁板定额。
平板厚度在100mm以内和以上时,也要分别计算和分别套相应定额。
板上开洞超过0.05时应扣除,但留洞口的工料应另列项目计算,计算单位以洞口个数计。
工程量计算公式如下:(1)平板100mm厚以内的:B1板=长乂宽X厚X块数=?立方+B2板=长乂宽X厚X块数=?立方+B3板=长乂宽X厚X块数=?立方+……(扣板上的洞)一工(洞口面积X板厚)=?立方(2)平板100mm以上的:B5板=长乂宽X厚X块数=?立方+B6板=长乂宽X厚X块数=?立方+B7板=长乂宽X厚X块数=?立方+……(扣板上的洞)一工(洞口面积X板厚)=?立方(3)有梁板:B1板=长乂宽X厚X块数=?立方+B2板=长乂宽X厚X块数=?立方+B3板=长乂宽X厚X块数=?立方+……(扣板上的洞)一工(洞口面积X板厚)=?立方附1:有梁板体积与板中梁体积合并,套有梁板定额。
15钢筋混凝土结构图详解

5
配筋图的绘制
2 钢筋混凝土结构的图示方法
2.2 配筋图
配筋平面图
对于水平放置,
纵、横向尺寸都比 较大的构件,通常 用平面图表示其配 筋情况。
2 钢筋混凝土结构的图示方法
2 钢筋混凝土结构的图示方法
2.2 配筋图
配配筋筋立平面面图图和断面图
对于细而长的构
件(如梁、柱)的 钢筋,常用配筋立 面图和若干断面图 表达。
2 钢筋混凝土结构的图示方法
2.1 模板图
模板图主要用来标注构件的外形尺寸及预埋件 的位置及编号,供制作模板使用,同时还可用 来计算混凝土的方量。 模板图一般就是构件的外形图。 形状简单的构件如矩形断面的梁、板、柱,可 不必单独绘制模板图,而与配筋图合并绘制。
2 钢筋混凝土结构的图示方法
2.2 配筋图
Ⅳ级钢筋强度太高,不适宜作为钢筋混凝土构件中 的配筋,一般冷拉后作预应力筋
1 钢筋混凝土结构的基本知识
1.6 钢筋骨架
受力钢筋:主要用来 承受外力。
架立钢筋:主要用来 固定箍筋和受力钢筋 的位置。主要用于钢 筋混凝土梁中。
箍筋:主要用来固定 受力钢筋的位置,也 承受一部分外力。
1 钢筋混凝土结构的基本知识
对于型号、直径、长 度和间距都相同的钢 筋,可以只画出第一 根和最末一根全长, 用标注的方法表明其 根数、规格和间距。
也可用粗实线画出其 中一根,同时用一两 端带斜短划线的横穿 细线表示其余钢筋及 其起止范围。
对于型号、直径和长 度都相同,而间距不 相同的钢筋,可只画 出第一根和最末一根 全长,中间用短粗线 表示其余钢筋的位置, 并用标注的方法表明 钢筋的根数、直径和 间距。
其它钢筋:如构造钢 筋、吊钩等。
中南院《混凝土结构计算图表》excel简化版

850~1000 1050~1300 >1300 1.07 1.11 1.06 1.08 1.04 1.06
数据输入
混凝土强度等级 (N/mm2) 底筋排数 钢筋等级 梁宽 (mm) 梁高 (mm) 钢筋数量 钢筋规格 (mm)
数据输出
混凝土抗压强度 (N/mm2) 保护层厚 (mm) 钢筋抗拉强度 (N/mm2) 钢筋弹性模量 (N/mm2) 相对极限受压区高度 极限配筋面积 (mm2) 底筋面0! 0
1.2 1.17 -
b
1
1 f
y
/(E s cu )
As bh0b f c / f y
M 1 f c bh 0 b (1 0 . 5 b )
2
制表公式: [M ] f
y
As (h0
f
y
As
c
2 bf
)
配三排钢筋 h-95 h-90
650 700~800 1.1 1.18 1.09 1.16
结论
抵抗矩 (kNm) 是否超筋 #DIV/0! #DIV/0! 极限弯矩M1 (kNm) 配筋率 (%) 0.00 #DIV/0!
输入底筋面积 (mm2) 重新计算抵抗矩 (kNm) #DIV/0!
表1 梁截面有效高度h0 混凝土强度等级 C20 C25,C30,C35,C40 配一排钢筋 h-40 h-35 表2 修改钢筋排数时的系数 梁高h 一排改二排 二排改三排 - 400 450 500 1.14 - - 550 1.13 - 600 1.11 配二排钢筋 h-65 h-60 配三排钢筋 h-95 h-90
钢筋混凝土构件基本知识和图示方法

§4.2 钢筋混凝土构件基本知识和图示方法混凝土是由水泥、砂子、石子和水按一定的比例拌和而成,灌入定型模板,经过振捣密实和养护凝固后坚硬如石,这种材料受压性能好,但受拉能力差,约为抗压强度的1/9-1/20。
钢筋混凝土是在混凝土受拉区域或相应部位加入适量钢筋,钢筋不但具有良好的抗拉性能,同时其热膨胀系数与混凝土接近,将两种材料黏结成整体,共同承受外力,这种材料解决了混凝土材料抗拉能力差的弊端。
用钢筋混凝土制成的梁、板、柱、基础等构件,称为钢筋混凝土构件。
图4-2为素混凝土与钢筋混凝土梁受力简图。
图4-2 素混凝土与钢筋混凝土梁受力比较一、混凝土的强度等级钢筋混凝土构件,有在工地现场浇制的,称为现浇钢筋混凝土构件。
也有在工厂或工地以外预先把构件制作好,然后运到工地安装的,称为预制钢筋混凝土构件。
此外还有的构件,制作时对混凝土预加一定的压力以提高构件的强度和抗裂性能,称为预应力钢筋混凝土构件。
混凝土的抗压强度,我国《规范》规定,用边长150mm标准立方体试块,在温度20°±3°C,相对湿度大于90%的环境中,养护28天后所测得的平均压力值为混凝土的抗压强度。
分为C7.5、C10、C15、C20、C25、C30、C35、C40、C45、C50、C55、C60,共12个等级,数字越大,表示混凝土抗压强度越高。
不同工程或用于不同部位的混凝土,对其强度标号的要求不一样。
一、结构说明1、结构形式(结构材料及类型;结构材料及规格、强度等级);2、地基与基础(包括地基土的地耐力等);3、施工技术要求及注意事项;4、选用的标准图集等。
二、结构布置平面图1、基础平面;2、楼层结构平面布置图;3、屋面结构平面布置图。
三、构件详图1、梁、板、柱、基础结构详图;2、楼梯结构详图;3、屋架(屋面)结构详图;4、其它详图,天沟、雨蓬、圈梁、过梁、门窗过梁、阳台、管道井、烟道井等。
特点:1、沿房屋防潮层的水平剖切表示基础平面图,沿每层楼板面水平剖切表示各层楼层结构平面图,沿屋面承重层的水平剖切表示屋面结构平面图。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
练习2:一、题目钢筋混凝土矩形截面:b=300mm ,h=600mm ,h 0=560mm ,a s ’=25mm ,a s =40mm ,A s ’=157mm2,A s =804mm2,f y ’=280MPa ,f y =280MPa ,E s =200GPa ,E c =25.5GPa ,f c =13.4MPa ,f t =1.54MPa ,ε0=0.002,εcu =0.0038,εs u ≤10%=0.10。
.利用数值方法计算截面的M~N 关系,并附简化计算结果N u 。
b=300mm2Φ10h=600mm4Φ16二、简单分析:本次作业是在上一次作业的基础上继续进行其他计算。
主要任务是利用计算机软件来计算特定截面偏心受压情况下的纵向压力与弯矩的关系。
可参考第一次作业的程序,做适当的修改即可。
混凝土应力—应变曲线采用的是R üsch 建议的曲线。
曲线的上升段采用抛物线形式,下降段为斜直线。
R üsch 建议的曲线:当0εε≤时,'2002[()]c cf εεσεε=-当0cu εεε<≤时,'c c f σ=根据《高等钢筋混凝土结构学》提供的公式: 由平衡条件:0N =∑,0M =∑可知,''0()()d x c ci s s s sN b y y A A σεσσ=+-⎰0'''00000()()()(+)d ()x s c ci i s s s N e y b y h x y y A h a σεσ+=-+-⎰对于离纵向力较远钢筋应力的取值可参照以下情况:平衡破坏:即受压边缘混凝土应变恰好达到极限应变时,受拉钢筋刚好达到屈服强度280MPa 。
受拉破坏:荷载偏心较大时,钢筋先屈服(达到280MPa ),经过一个过程后,混凝土达到极限压应变。
受压破坏:荷载偏心较小时,构件产生受压破坏。
受压破坏是指受压较大一侧的混凝土达到极限压应变,而离纵向力较远一侧的钢筋可能受拉或者受压但都不屈服。
此时钢筋应力可用以下代替:(1)s s s cu sh E E x σεε==-轴心受压时,e0=0,全截面受压且破坏时压应变均为0.0038,两侧钢筋均受压屈服,代入上式中可以得到'''13.43006002801572808072681.92c s y s y N b h A f A f kNσ=⨯⨯+⨯+⨯=⨯⨯+⨯+⨯='''00()()2242.60c s y s s h b h h A f h a y Nmmσ⨯⨯⨯-+⨯⨯-==. 由c++最终计算的数据所得到的N-M 图如下图所示。
计算出来的纯弯弯矩为近似等于121.731kN·m与手算结果120.89kN·m结果吻合。
实际上c++计算时,得到的轴心受压极限纵向力以及纯弯的极限弯矩其分别对应的弯矩和纵向力并不为0,只是值很小接近0。
所以这种计算只是近似计算。
x的步长越小所得的结果就越精确。
平衡破坏时,弯矩为266.485kN·m,纵向压力为1068.62kN。
轴压破坏时极限压力为2683.49kN与理论计算结果2687.92kN吻合。
三、程序设计#include<iostream>#include<math.h>#include<fstream>using namespace std;int main(){cout<<"设计中As=804mm2,As'=157mm2,fy=280MPa,fy'=280MPa,Es=200GPa,Ec=25.5GPa"<<endl;cout<<endl;cout<<"fc=13.4MPa,ft=1.54MPa,ε0=0.002,εcu=0.0038,εsu<=0.1"<<endl;cout<<endl;//给出题目的基本信息double b,h,as2,as1,x0,c,t,N,h0,x1,q,delta_epsilon_cc,delta_epsilon_cs,p;doublek,epsilon_ss,epsilon_sc,epsilon_cu,sigma_ss,sigma_sc,M,epsilon_c,m0, epsilon_c1,epsilon_cs,epsilon_c2;ofstream outfile;b=300;h=600;as1=40;as2=25; h0=h-as1; epsilon_cu=0.0038;//给出题目相关参数,as1为受拉区的钢筋保护层厚度,as2为受压区钢筋保护层厚度。
outfile.open("data.txt");//建立数据输出文件for(x0=0.0;x0<=10000;x0+=1){epsilon_sc=epsilon_cu/x0*(x0-as2);epsilon_ss=epsilon_cu/x0*(h0-x0 );//根据几何关系求出受压区和受拉区钢筋的应变。
下面可得相应应力sigma_sc=200*epsilon_sc*1000;sigma_ss=200*epsilon_ss*1000;if(sigma_ss>280)sigma_ss=280;if(sigma_ss<-280)sigma_ss=-280;if(sigma_sc>280)sigma_sc=280;p=0.0;m0=0.0;k=0.0;N=0;M=0;if(x0>600)//受压区高度超过梁高情况下应力分布,以及纵向力和弯矩求解{ for(int i=0;i<=1000;i++){ epsilon_c1=600*epsilon_cu/x0;epsilon_c=epsilon_c1*i/1000+epsilon_cu-epsilon_c1;if(epsilon_c<=0.002){k=13.4*(1000*epsilon_c-epsilon_c*epsilon_c/0.000004); }else k=13.4; //求出不同应变时受压区各个位置的压应力大小p=p+k*b*h/1000;m0=m0+k*b*h/1000*(h*(i+0.5)/1000-as1);}//受压区混凝土对离纵向受拉区钢筋取弯矩,不过此时受拉区钢筋也受压N=p+sigma_sc*157-sigma_ss*804;M=m0+sigma_sc*157*(h0-as2)-N*242.60;if(N>=0&&M>=0)//对N和M进行筛选,并存入txt文件里outfile<<N/1000<<" "<<M/1e+06<<endl;//统一成kN,m单位}epsilon_c2=(h-x0)*epsilon_cu/x0;epsilon_cs=1.54/25.5/1000;//受拉混凝土弹性极限应变if(x0<=600)//受压区高度未超过梁高情况应力分布,以及纵向力和弯矩求解{if(epsilon_c2>epsilon_cs)//用来判断是否开裂,此时开裂。
{x1=epsilon_cs*x0/epsilon_cu;for(int i=0;i<=1000;i++){ delta_epsilon_cc=epsilon_cu*(i+0.5)/1000;if(delta_epsilon_cc<=0.002){k=13.4*(1000*delta_epsilon_cc-delta_epsilon_cc*delta_epsilon_cc/0.000004);} else {k=13.4;}//求出不同应变时受压区各个位置的压应力大小delta_epsilon_cs=epsilon_cs*(i+0.5)/1000;q=delta_epsilon_cs*25.5*1000;p=p+k*b*x0/1000-q*b*x1/1000;m0=m0+k*b*x0/1000*(x0*(i+0.5)/1000+h0-x0)-q*b*x1/1000*(-(i+0.5)*x1/ 1000)+h0-x0;}//对纵向受拉钢筋求矩N=p+sigma_sc*157-sigma_ss*804;M=m0+sigma_sc*157*(h0-as2)-N*242.60;if(N>=0&&M>=0)//对N和M进行筛选,并存入txt文件里outfile<<N/1000<<" "<<M/1e+06<<endl;}//统一成kN,m单位if(epsilon_c2<epsilon_cs)//未开裂{for(int i=0;i<=1000;i++){ delta_epsilon_cc=epsilon_cu*(i+0.5)/1000;if(delta_epsilon_cc<=0.002){k=13.4*(1000*delta_epsilon_cc-delta_epsilon_cc*delta_epsilon_cc/0.000004);} else {k=13.4;}//求出不同应变时受压区各个位置的压应力大小delta_epsilon_cs=epsilon_c2*(i+0.5)/1000;q=delta_epsilon_cs*25.5*1000;p=p+k*b*x0/1000-q*b*(h-x0)/1000;m0=m0+k*b*x0/1000*(x0*(i+0.5)/1000+h0-x0)-q*b*(h-x0)/1000*(-(h-x0)* (i+0.5)/1000+h0-x0);}//对受压区受拉区应力零点求弯矩N=p+sigma_sc*157-sigma_ss*804;M=m0+sigma_sc*157*(h0-as2)-N*242.60;}if(N>=0&&M>=0)//对N和M进行筛选,并存入txt文件里outfile<<N/1000<<" "<<M/1e+06<<endl;//统一成kN,m单位}}outfile.close();system("pause");return 0;//程序结束}。