320飞行控制从系统可靠性.rar_可靠性_可靠性分析_可靠性计算_可靠指标_客机

%%%%%%%%%%左升降舵的可靠性指标计算 底事件有：舵面损伤、ELAC1失效、ELAC2失效、B液压源失效、B液压传动装置失效、G液压源失效、G液压传动装置失效、 clc; clear; clear all; close all; %%%%%%%%%%%%%%%%%%%%%%%%%%1%%%%%%2%%%%%%%%3%%%%%%4%%%%%%5%%%%%%6%%%%%%%%%%7%%-----8----------9--------10------11--------12--- ---13---------14------15-------16---------17----18---- -19--------20-----21----22-----23------24-----25--------26------27-----28-----29 %--------------- ----%左副翼--左副翼B-左副翼G-右副翼—右副翼G-右副翼B---左升降舵--左升降舵B-左升降舵G-右升降舵--右升降舵Y--右升降舵B-方向舵-----方向舵G- 方向舵Y--方向舵B-- -马达1--马达2---ELAC1---ELAC2--SEC1---SEC2---SEC3---液压源B-液压源G---液压源Y- -脚蹬-----FAC1--FAC2 A320_fail_rate=[ 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 ];%每个底事件失效率 %%左副翼 1 19 24 2/ 20 25 3 %%右副翼 4 19 25 5/ 20 24 6 %%左升降舵 7 20 25 9/ 19 24 8/ 22 25 9/21 24 8 %%右升降舵 10 20 26 11/19 24 12/22 26 11/21 24 12 %%方向舵 13 25 14/24 16/26 15/28 17/29 18/27 A320_num_d=29;%底事件的个数 N=10000;%总的仿真次数 T=3000;%仿真时间，单位小时 M=3000;%把仿真时间分成多少段 delt_t=T/M;%平均每段长度 for i=1:M l_fuyi_f_num(i)=0;%计数器 l_fuyi_R(i)=0;%可靠度 l_fuyi_F(i)=0;%不可靠度 l_fuyi_P(i)=0;%系统失效概率分布 l_fuyi_f(i)=0;%系统失效概率密度函数 l_fuyi_y(i)=0;%系统失效率 r_fuyi_f_num(i)=0;%计数器 r_fuyi_R(i)=0;%可靠度 r_fuyi_F(i)=0;%不可靠度 r_fuyi_P(i)=0;%系统失效概率分布 r_fuyi_f(i)=0;%系统失效概率密度函数 r_fuyi_y(i)=0;%系统失效率 r_shenjiang_f_num(i)=0;%计数器 r_shenjiang_R(i)=0;%可靠度 r_shenjiang_F(i)=0;%不可靠度 r_shenjiang_P(i)=0;%系统失效概率分布 r_shenjiang_f(i)=0;%系统失效概率密度函数 r_shenjiang_y(i)=0;%系统失效率 l_shenjiang_f_num(i)=0;%计数器 l_shenjiang_R(i)=0;%可靠度 l_shenjiang_F(i)=0;%不可靠度 l_shenjiang_P(i)=0;%系统失效概率分布 l_shenjiang_f(i)=0;%系统失效概率密度函数 l_shenjiang_y(i)=0;%系统失效率 l_direction_f_num(i)=0;%计数器 l_direction_R(i)=0;%可靠度 l_direction_F(i)=0;%不可靠度 l_direction_P(i)=0;%系统失效概率分布 l_direction_f(i)=0;%系统失效概率密度函数 l_direction_y(i)=0;%系统失效率 A320_f_num(i)=0;%控制系统失效计数器 A320_R(i)=0;%可靠度 A320_F(i)=0;%不可靠度 A320_P(i)=0;%系统失效概率分布 A320_f(i)=0;%系统失效概率密度函数 A320_y(i)=0;%系统失效率 A320_f_num1(i)=0;%控制系统失效计数器 A320_R1(i)=0;%可靠度 A320_F1(i)=0;%不可靠度 A320_P1(i)=0;%系统失效概率分布 A320_f1(i)=0;%系统失效概率密度函数 A320_y1(i)=0;%系统失效率 A320_f_num2(i)=0;%控制系统失效计数器 A320_R2(i)=0;%可靠度 A320_F2(i)=0;%不可靠度 A320_P2(i)=0;%系统失效概率分布 A320_f2(i)=0;%系统失效概率密度函数 A320_y2(i)=0;%系统失效率 end for i=1:N l_shenjiang_n_mttf(i)=0;%记录每次仿真的左升降舵失效时间 r_shenjiang_n_mttf(i)=0;%记录每次仿真的右升降舵失效时间 l_direction_n_mttf(i)=0;%记录每次仿真的方向舵失效时间 l_fuyi_n_mttf(i)=0;%记录每次仿真的左副翼失效时间 r_fuyi_n_mttf(i)=0;%记录每次仿真的右副翼失效时间 A320_n_mttf(i)=0;%记录每次仿真的飞控系统失效时间 A320_n_mttf1(i)=0;%记录每次仿真的飞控系统失效时间 A320_n_mttf2(i)=0;%记录每次仿真的飞控系统失效时间 end l_shenjiang_s_t=0; r_shenjiang_s_t=0; l_fuyi_s_t=0; r_fuyi_s_t=0; l_direction_s_t=0; %%%判断标志位 flg_l_fuyi=0; flg_r_fuyi=0; flg_l_shenjiang=0; flg_r_shenjiang=0; flg_direction=0; % lshenjiang_mttf=0;%左升降舵平均寿命 for i=1:N A320_r_temp=rand(1,A320_num_d); for j=1:A320_num_d A320_t(j)=-1/A320_fail_rate(j)*log(1-A320_r_temp(j)); end l_fuyi_t_temp_1=min([A320_t(19) A320_t(24) A320_t(2)]);%%左副翼失效时间 l_fuyi_t_temp_2=min([A320_t(20) A320_t(25) A320_t(3)]); l_fuyi_t_temp_3=min(A320_t(1),max(l_fuyi_t_temp_1,l_fuyi_t_temp_2)); l_fuyi_s_t=l_fuyi_t_temp_3; if(A320_t(2)<A320_t(3)) if(l_fuyi_t_temp_3==A320_t(3)) if(rand(1)>0.5) flg_l_fuyi=1; end end else if(l_fuyi_t_temp_3==A320_t(2)) if(rand(1)>0.5) flg_l_fuyi=1; end end end if(l_fuyi_s_t<=T) l_fuyi_f_num(ceil(l_fuyi_s_t/delt_t))=l_fuyi_f_num(ceil(l_fuyi_s_t/delt_t))+1; l_fuyi_n_mttf(i)=l_fuyi_s_t; else l_fuyi_n_mttf(i)=T; end r_fuyi_t_temp_1=min([A320_t(19) A320_t(25) A320_t(5)]);%%右副翼失效时间 r_fuyi_t_temp_2=min([A320_t(20) A320_t(24) A320_t(6)]); r_fuyi_t_temp_3=min(A320_t(1),max(r_fuyi_t_temp_1,r_fuyi_t_temp_2)); r_fuyi_s_t=r_fuyi_t_temp_3; if(A320_t(5)<A320_t(6)) if(r_fuyi_s_t==A320_t(6)) if(rand(1)>0.5) flg_r_fuyi=1; end end else if(r_fuyi_s_t==A320_t(5)) if(rand(1)>0.5) flg_r_fuyi=1; end end end if(r_fuyi_s_t<=T) r_fuyi_f_num(ceil(r_fuyi_s_t/delt_t))=r_fuyi_f_num(ceil(r_fuyi_s_t/delt_t))+1; r_fuyi_n_mttf(i)=r_fuyi_s_t; else r_fuyi_n_mttf(i)=T; end l_shenjiang_t_temp_1=min([A320_t(20) A320_t(25) A320_t(9)]);%左升降舵失效时间 l_shenjiang_t_temp_2=min([A320_t(19) A320_t(24) A320_t(8)]); l_shenjiang_t_temp_3=min([A320_t(22) A320_t(26) A320_t(11)]); l_shenjiang_t_temp_4=min([A320_t(21) A320_t(24) A320_t(8)]); l_shenjiang_t_temp_5=max([l_shenjiang_t_temp_1 l_shenjiang_t_temp_2 l_shenjiang_t_temp_3 l_shenjiang_t_temp_4]); l_shenjiang_s_t=min(l_shenjiang_t_temp_5,A320_t(7)); if(A320_t(8)<A320_t(9)) if(l_shenjiang_s_t==A320_t(9)) if(rand(1)>0.5) flg_l_shenjiang=1; end end else if(l_shenjiang_s_t==A320_t(8)) if(rand(1)>0.5) flg_l_shenjiang=1; end end end if(l_shenjiang_s_t<=T) l_shenjiang_f_num(ceil(l_shenjiang_s_t/delt_t))=l_shenjiang_f_num(ceil(l_shenjiang_s_t/delt_t))+1; l_shenjiang_n_mttf(i)=l_shenjiang_s_t; else l_shenjiang_n_mttf(i)=T; end r_shenjiang_t_temp_1=min([A320_t(20) A320_t(26) A320_t(11)]);%右升降舵失效时间 r_shenjiang_t_temp_2=min([A320_t(19) A320_t(24) A320_t(12)]); r_shenjiang_t_temp_3=min([A320_t(22) A320_t(25) A320_t(9)]); r_shenjiang_t_temp_4=min([A320_t(21) A320_t(24) A320_t(12)]); r_shenjiang_t_temp_5=max([r_shenjiang_t_temp_1 r_shenjiang_t_temp_2 r_shenjiang_t_temp_3 r_shenjiang_t_temp_4]); r_shenjiang_s_t=min(r_shenjiang_t_temp_5,A320_t(10)); if(A320_t(11)<A320_t(12)) if(r_shenjiang_s_t==A320_t(12)) if(rand(1)>0.5) flg_r_shenjiang=1; end end else if(r_shenjiang_s_t==A320_t(11)) if(rand(1)>0.5) flg_r_shenjiang=1; end end end if(r_shenjiang_s_t<=T) r_shenjiang_f_num(ceil(r_shenjiang_s_t/delt_t))=r_shenjiang_f_num(ceil(r_shenjiang_s_t/delt_t))+1; r_shenjiang_n_mttf(i)=r_shenjiang_s_t; else r_shenjiang_n_mttf(i)=T; end direction_t_temp_1=min([A320_t(25) A320_t(14)]);%%方向舵 direction_t_temp_2=min([A320_t(24) A320_t(16)]); direction_t_temp_3=min([A320_t(26) A320_t(15)]); direction_t_temp_4=max([direction_t_temp_1 direction_t_temp_2 direction_t_temp_3]); direction_t_temp_5=min([A320_t(28) A320_t(17)]); direction_t_temp_6=min([A320_t(29) A320_t(18)]); direction_t_temp_7=direction_t_temp_5+direction_t_temp_6; direction_t_temp_8=max(direction_t_temp_7,A320_t(27)); l_direction_s_t=min([dire