MATLAB编写的预应力CFRP加固程序，请各位指正完善!

shihonglei

先对混凝土梁加载，然后粘贴预应力CFRP加固，计算全过程曲线。本人编写调试了N多次，总是难如意：
1、加载到初始弯矩M0时，想得到预应力CFRP放张前后的曲线突变，可是实现不了。求助!
如果第一点是因为编程水平低的原故，那么第二点更令我困惑!
2、弯矩--曲率关系曲线在M0后应该是单调上升的，可是中间怎么会出现波动性的上升？？
请各位XDJM改改调试调试!!


我是根据条带划分法来编程的，具体步聚如下：
整个计算程序按以上分析分为三大步：
第一步：分级加变形计算 ，直到弯矩达到初始弯矩 ；
第二步：预应力碳纤维加固后截面的 计算；
第三步：继续增加变形计算 ，直到满足截面破坏条件而停止。
具体程序步骤：
（1）取曲率： ；
（2）假定梁截面中心轴处混凝土的应变 ；
（3）求各混凝土条带和钢筋、碳纤维布的应变；
（4）按它们各自的应力－应变关系求与应变相对应的应力值；
（5）把各条带的内力求和，判断是否满足平衡条件；
（6）若不满足平衡条件，则需要调整中心轴处的应变 值，重复（3）－（5）；
（7）满足平衡条件后，求得与 相对应的弯矩值 ；
（8）判断截面弯矩 是否达到加固前的初始弯矩 ，若 ＜M0 ，重复（1）－（7）；若 >M0 ，转入第（9）步；
（9）每次取曲率 ；
（10）重复（2）－（7）步，判断是否满足预应力碳纤维布放张后的平衡条件，若不满足，转入第（9）；若满足，转入第（11）步；
（11）重复（1）－（7）步，即继续增加变形计算，直到满足截面破坏条件而停止。
        
        你可以对照程序看，程序后也注解 的。

       你也可以参考 吕西林编的<钢筋混凝土结构非线性有限元理论与应用>，书中有用C语言编写的条带划分法计算混凝土梁的全过程曲线!

shihonglei

大家帮帮忙，谢谢了!

H.Q.Yang

感谢分享