混凝土Matlab程序分析
基于Matlab的 M-φ全过程分析43号 辛明明 研1504班 Matlab主程序Matlab子程序程序运行流程分析思路此程序采用的本构关系钢筋:理想弹塑性加硬化(三折线式)混凝土(压):德国H.Usch 混凝土(拉):斜直线图像成果问题回答:1、钢筋混凝土截面承载力分析中的基本自变量是受压区高度和钢筋面积,这里的Matlab程序中的基本自变量是什么?(即解释s(i=1~660)的做法)答:这里的自变量是s(受压区边缘的压应变),每次取混凝土受压区边缘的压应变s=s+▽s,然后假定受压区高度x0,再根据其他条件计算。2、试解释弯矩的组成,并分述各弯矩组成的计算实现。答:弯矩由受压区混凝土合力对中性轴的力矩、受拉区混凝土合力对中性轴的力矩、受压钢筋对中性轴力矩、受拉钢筋对中性轴力矩四种弯矩。其中受压(拉)区混凝土合力对中性轴的力矩是将受压(拉)区的混凝土平均分成若干层,然后根据混凝土受压(拉)时的应力应变关系求出第k层的受压(拉)区的混凝土应力,再求和得到受压(拉)区混凝土的合力,再对中性轴取矩的到;受压(拉)钢筋对中性轴力矩是由钢筋的应力应变关系求出受压(拉)钢筋的应力再对中性轴取矩得到。问题回答:3、二分法的基本原理?在本计算中的使用?答:基本原理——若 f ? C[a, b],且 f (a) · f (b) < 0,则 f 在 (a, b) 上必有一根。 具体做法就是通过二等分不断缩小有根区间的长度,直到满足精度为止。在本计算中的应用给定有根区间 [a, b] ( f(a) · f(b) < 0) (即混凝土分层的厚度)和 精度要求 ?1. 令 x = (a+b)/22. 如果 b – a <= 2?, 停止计算,输出 x ,否则执行第3步3. 如果 f (a) f (x) < 0 , 则令 b = x,否则令 a = x, 返回第1步问题回答:受弯钢筋截面面积 As=M/fy*rs*h0(M是正截面弯矩设计值、fy是纵向钢筋的抗拉强度设计值、rs内力臂系数、h0截面有效高度。4、如何把主程序编为子程序?——弯矩、曲率和钢筋面积(或者截面高度)的关系。答:将主程序制成.m文件放到工作目录上去,就可以直接运行(需要用的时候直接调用函数即可)。弯矩、曲率和钢筋面积(或者截面高度)的关系 THANKS 谢谢观看模板来自于