对如下可靠性问题进行仿真.doc

《对如下可靠性问题进行仿真.doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、对如下可靠性问题进行仿真：一设备上三个相同的轴承，每个轴承正常工作寿命为随机变量，其概率分布如表所示．寿命h1000110012001300140015001600170018001900概率0.100.130.250.130.090.120.020.060.050.05在任何一个轴承损坏都可以使设备停止工作，从有轴承损坏，设备停止工作，到检修工到达开始更换部件为止，称为一个延迟时间．延迟时间也是随机变量，其概率分布如下表所示．延时时间min51015概率0.60.30.1设备停工时每分钟损失5元,检修工每小时工时费12元，轴承每个成本16元．更换一个轴承需要

2、20min，同时更换两个轴承需要30min，同时更换三个轴承需要40min．现在有两种方案：方案一是损坏一个更换一个；方案二是一旦有轴承损坏就全部更换．试通过计算机仿真对这两种方案做出评价．方案1%事件步长法clearcloseclct=0;%仿真时钟,成本一成本二cost1=0;cost2=0;%产生初始化事件表p1=rand(1,3);%产生三个轴承坏的时间概率p2=rand(1,3);%产生三个延迟时间概率fori=1:3ifp1(i)<0.1ta(i)=1000;elseifp1(i)<0.23ta(i)=1100;elseifp1(i)<0.48ta

3、(i)=1200;elseifp1(i)<0.61ta(i)=1300;elseifp1(i)<0.7ta(i)=1400;elseifp1(i)<0.82ta(i)=1500;elseifp1(i)<0.84ta(i)=1600;elseifp1(i)<0.9ta(i)=1700;elseifp1(i)<0.95ta(i)=1800;elseta(i)=1900;endifp2(i)<0.6tyc(i)=5;elseifp2(i)<0.9tyc(i)=10;elsetyc(i)=15;endtsj(i)=ta(i)+tyc(i);endt=min(tsj);

4、forj=1:3iftsj(j)==ti=j;break;endendcost1=cost1+(tyc(i)+20)*5+4+16;%仿真whilet<=100000p1(i)=rand(1);p2(i)=rand(1);ifp1(i)<0.1ta(i)=1000;elseifp1(i)<0.23ta(i)=1100;elseifp1(i)<0.48ta(i)=1200;elseifp1(i)<0.61ta(i)=1300;elseifp1(i)<0.7ta(i)=1400;elseifp1(i)<0.82ta(i)=1500;elseifp1(i)<0.84

5、ta(i)=1600;elseifp1(i)<0.9ta(i)=1700;elseifp1(i)<0.95ta(i)=1800;elseta(i)=1900;endifp2(i)<0.6tyc(i)=5;elseifp2(i)<0.9tyc(i)=10;elsetyc(i)=15;endtsj(i)=tsj(i)+ta(i)+tyc(i);t=min(tsj);forj=1:3iftsj(j)==ti=j;break;endendcost1=cost1+(tyc(i)+20)*5+4+16;end[t,cost1]ans=10053034450方案2%事件步长

6、法clearcloseclct=0;%仿真时钟,成本一成本二cost2=0;tsj=[000];whilet<=100000p1=rand(1,3);%产生三个轴承坏的时间概率p2=rand(1,3);%产生三个延迟时间概率fori=1:3ifp1(i)<0.1ta(i)=1000;elseifp1(i)<0.23ta(i)=1100;elseifp1(i)<0.48ta(i)=1200;elseifp1(i)<0.61ta(i)=1300;elseifp1(i)<0.7ta(i)=1400;elseifp1(i)<0.82ta(i)=1500;elseifp

7、1(i)<0.84ta(i)=1600;elseifp1(i)<0.9ta(i)=1700;elseifp1(i)<0.95ta(i)=1800;elseta(i)=1900;endifp2(i)<0.6tyc(i)=5;elseifp2(i)<0.9tyc(i)=10;elsetyc(i)=15;endtsj(i)=tsj(i)+ta(i)+tyc(i);endt=min(tsj);forj=1:3iftsj(j)==ti=j;break;endendcost2=cost2+(tyc(i)+40)*5+8+3*16;end[t,cost2]ans=10108

8、022531第二种方案好实验四蒙特卡洛