kent-park模型_上叶侬

http://blog.sina.com.cn/u/1098546544

首页博文目录关于我

个人资料

微博

加好友发纸条

写留言加关注

博客等级：
博客积分：

博客访问：
关注人气：
获赠金笔：0支
赠出金笔：0支
荣誉徽章：

正文字体大小：大中小

kent-park模型

(2016-12-19 20:25:08)

标签：

abaqus二次开发

分类：我的学习生活

http://s7/mw690/001clo8Ezy77jDpjhSmc6&690

http://s6/mw690/001clo8Ezy77jDvkYzH45&690

%Kent_Park应力-应变全曲线
clear all;
clc;
fcR=[];%混凝土的圆柱体抗压强度；
Ps=[];%横向箍筋对核心混凝土（去箍筋外皮以内）的体积率；
bb=[];%从箍筋外皮
s=[];%箍筋间距
fcR=input('please input the conpressive strength of round column:');
ep=0:0.0001:0.006;
epp=0.002;
len=length(ep);
for i=1:len;
    epc(i)=ep(i)/0.002;
    if epc(i)<=1;
       fc(i)=fcR*(2*epc(i)-epc(i)^2);
       fcmax=max(fc);
    else
        ep5=(20.67+2*fcR)/(fcR-6.89)*1e-3;
        fc(i)=(0.5*fcmax-fcmax)/(ep5-1)*(epc(i)-1)+fcR;
        �(i)=0.5*fcR*(epp-epc(i))/(ep5-epp)+fcR
        if fc(i)<</span>=0.2*fcR;
           fc(i)=0.2*fcR;
        end
    end
end
%绘图；
plot(ep,fc);
grid on;
title('Kent Park modle');
xlabel('应变（单位：10^-3）');
ylabel('应力（单位：N/mm^2）');

阅读┊ 收藏 ┊ 喜欢 ▼ ┊打印┊举报/Report

前一篇：[转载]Fortran中的常用函数

后一篇：[转载]关于纤维模型的一些想法

新浪BLOG意见反馈留言板　欢迎批评指正