JC法可靠度计算程序

JC法可靠度计算程序，时变可靠度计算实例以及时变可靠度曲线绘制，使用寿命预测

clear;
clc;
fun=‘24726*x_1+45960*x_2-70686*x_3-70686*x_4‘;%%%%%%%%%功能函数形式
beta=[];
mark=[1;1;2;3];%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%随机变量类型
for t=1:1:50%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%所需求解时间点
    mu_1=1.3988*exp（-0.0195*（exp（0.75*t）-1）^2）*40;
    sigma_1=（0.318*0.75*t+0.9881）*5;
    mu_2=（1-0.000015*（0.75*t）^3）*400;
    sigma_2=（1+0.289*0.75*t）*20;
    mu_3=15;
    sigma_3=10;
    mu_4=40;
    sigma_4=5;
    muX=[mu_1;mu_2;mu_3;mu_4];
    sigmaX=[sigma_1;sigma_2;sigma_3;sigma_4];
    beta=[beta;JC_method_ty（muXsigmaXmarkfun）];
end
plot（1:1:50beta‘--‘‘linewidth‘2）;
hold on
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%以上为不加固工况下的时变可靠度
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%以下为加固工况下的时变可靠度
bbeta=[];
fun=‘24726*x_1+45960*x_2-70686*x_3-70686*x_4‘;%%%%%%%%%功能函数形式
t=15;
    mu_1=1.3988*exp（-0.0195*（exp（0.75*t）-1）^2）*40;
    sigma_1=（0.318*0.75*t+0.9881）*5;
    mu_2=（1-0.000015*（0.75*t）^3）*400;
    sigma_2=（1+0.289*0.75*t）*20;
    mu_3=15;
    sigma_3=10;
    mu_4=40;
    sigma_4=5;
    muX=[mu_1;mu_2;mu_3;mu_4];
    sigmaX=[sigma_1;sigma_2;sigma_3;sigma_4];
    bbeta=[bbeta;JC_method_ty（muXsigmaXmarkfun）];
fun=‘1.005*24726*x_1+1.005*45960*x_2-70686*x_3-70686*x_4‘;%%%%%%%%%功能函数形式
for t=15:1:55%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%加固点之后
    mu_1=1.3988*exp（-0.0195*（exp（0.75*（t-14））-1）^2）*40;
    sigma_1=（0.318*0.75*（t-14）+0.9881）*5;
    mu_2=（1-0.000015*（0.75*（t-14））^3）*400;
    sigma_2=（1+0.289*0.75*（t-14））*20;
    mu_3=15;
    sigma_3=10;
    mu_4=40;
    sigma_4=5;
    muX=[mu_1;mu_2;mu_3;mu_4];
    sigmaX=[sigma_1;sigma_2;sigma_3;sigma_4];
    bbeta=[bbeta;JC_method_ty（muXsigmaXmarkfun）];
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%绘图%%%%%%%%%%%%%%%%%%%%%%%%%
t=[1515:1:55];
plot（tbbeta‘r‘‘linewidth‘2）
xlabel（‘时间t（年）‘）
ylabel（‘可靠度\beta‘）
hold on
t=[0:1:60];
bbeta=2.5*ones（161）;
plot（tbbeta‘k‘）
legend（‘不加固短柱‘‘加固短柱‘‘允许可靠度指标‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2105  2020-09-22 05:50  JC_Calculate\beta_time.m
     文件        1274  2020-09-22 05:50  JC_Calculate\beta_time1.m
     文件        3516  2020-09-22 05:50  JC_Calculate\JC_method_ty.m
     文件        1382  2020-09-22 05:50  JC_Calculate\│╠╨≥╦╡├≈.txt
     文件      138240  2020-09-22 05:50  JC_Calculate\╩╣╙├╩┘├ⁿ╝╞╦π╩╡└².doc