| 查看: 401 | 回复: 0 | ||
1607210637金虫 (正式写手)
|
[求助]
求大神matlab程序帮助
|
|
uD1=input(' 法兰截面内径的均值/mm μD1=' ); uD11=input(' 法兰截面内径的标准差/mm σD1='); uD0=input(' 法兰截面中径的均值/mm μD0=' ); uD00=input(' 法兰截面中径的标准差/mm σD0='); uh=input(' 法兰截面厚度的均值/mm μh='); bh=input(' 法兰截面厚度的标准差/mm σh=' ); uF=input(' 作用在法兰上载荷的均值/N μF=' ); bF=input(' 作用在法兰上载荷的标准差/N σF='); ur=input(' 法兰拉伸疲劳极限的均值/MPa ur=' ); br=input(' 法兰拉伸疲劳极限的标准差/MPa σr='); R=input(' 输入可靠度 R=' ); b=input(' 法兰截面尺寸的相关系数 ρ=' ); p=norminv(1-R); F=normcdf(p); fprintf(1,'可靠性指标 β=%3.6f\n',p); f=normpdf(p); A=3*uF*(uD0-uD1)/3.1415926*uD1+3*uF*uD0*uD11^2/3.1415926*uD1^3+9*uF*(uD0-uD1)*(0.005)^2/3.1415926*uD1; B=(-3*uF/3.1415926*uD1)^2*uD00^2+(3*uF*uD0/3.1415926*uD1^2)^2*uD11^2+[6*uF*(uD0-uD1)/3.1415926*uD1]^2*(0.005)^2+ [-3*uF*(uD0-uD1)/3.1415926*uD1]^2*p^2; c4=ur^2-p^2*br^2; fprintf(1,'联结方程多项式中四次项的系数 c4=%3.4f\n',c4) c2=-2*A*ur; fprintf(1,'联结方程多项式中二次项的系数 c2=%3.4f\n',c2) c0=A^2-p^2*B; fprintf(1,'联结方程多项式中常数项 c0=%3.4f\n',c0) y=[c4 0 0 0 c2 0 0 0 c0]; h=roots(y); fprintf(1,' h1=%3.3fmm\n',h(1)) fprintf(1,' h2=%3.3fmm\n',h(2)) fprintf(1,' h3=%3.3fmm\n',h(3)) fprintf(1,' h4=%3.3fmm\n',h(4)) RR=[0.5 0.6 0.7 0.8 0.9 0.99 0.999 0.9999 0.99999 0.999999]; pp=norminv(1-RR); figure(1); plot(pp,RR) title('\bf可靠性指标β与可靠度R的关系曲线\rmR=\int e^{-p^{2}/2}dp') xlabel('可靠性指标β') ylabel('可靠度R') grid; A=3*uF*(uD0-uD1)/3.1415926*uD1+3*uF*uD0*uD11^2/3.1415926*uD1^3+9*uF*(uD0-uD1)*(0.005)^2/3.1415926*uD1; B=(-3*uF/3.1415926*uD1)^2*uD00^2+(3*uF*uD0/3.1415926*uD1^2)^2*uD11^2+[6*uF*(uD0-uD1)/3.1415926*uD1]^2*(0.005)^2+ [-3*uF*(uD0-uD1)/3.1415926*uD1]^2*p^2; cc2=-2*A*ur; cc0=A^2-pp^2*B; y1=[cc4(1) 0 0 0 cc2(1) 0 0 0 cc0(1)]; hr1=roots(y1); fprintf(1,'可靠度R=0.5000 时的可靠性指标 β1=%3.6f\n',pp(1)); fprintf(1,' 法兰的厚度 hr1=%3.3fmm\n',hr1(4)); y2=[cc4(2) 0 0 0 cc2(2) 0 0 0 cc0(2)]; hr2=roots(y2); fprintf(1,'可靠度R=0.6000 时的可靠性指标 β2=%3.6f\n',pp(2)); fprintf(1,' 法兰的厚度 hr2=%3.3fmm\n',hr2(4)); y3=[cc4(3) 0 0 0 cc2(3) 0 0 0 cc0(3)]; hr3=roots(y3); fprintf(1,'可靠度R=0.7000 时的可靠性指标 β3=%3.6f\n',pp(3)); fprintf(1,' 法兰的厚度 hr3=%3.3fmm\n',hr3(4)); y4=[cc4(4) 0 0 0 cc2(4) 0 0 0 cc0(4)]; hr4=roots(y4); fprintf(1,'可靠度R=0.8000 时的可靠性指标 β4=%3.6f\n',pp(4)); fprintf(1,' 法兰的厚度 hr4=%3.3fmm\n',hr4(4)); y5=[cc4(5) 0 0 0 cc2(5) 0 0 0 cc0(5)]; hr5=roots(y5); fprintf(1,'可靠度R=0.9000 时的可靠性指标 β5=%3.6f\n',pp(5)); fprintf(1,' 法兰的厚度 hr5=%3.3fmm\n',hr5(4)); y6=[cc4(6) 0 0 0 cc2(6) 0 0 0 cc0(6)]; hr6=roots(y6); fprintf(1,'可靠度R=0.9900 时的可靠性指标 β6=%3.6f\n',pp(6)); fprintf(1,' 法兰的厚度 hr6=%3.3fmm\n',hr6(4)); y7=[cc4(7) 0 0 0 cc2(7) 0 0 0 cc0(7)]; hr7=roots(y7); fprintf(1,'可靠度R=0.9990 时的可靠性指标 β7=%3.6f\n',pp(7)); fprintf(1,' 法兰的厚度 hr7=%3.3fmm\n',hr7(4)); y8=[cc4(8) 0 0 0 cc2(8) 0 0 0 cc0(8)]; hr8=roots(y8); fprintf(1,'可靠度R=0.9999 时的可靠性指标 β8=%3.6f\n',pp(8)); fprintf(1,' 法兰的厚度 hr8=%3.3fmm\n',hr8(4)); y9=[cc4(9) 0 0 0 cc2(9) 0 0 0 cc0(9)]; hr9=roots(y9); fprintf(1,'可靠度R=0.99999 时的可靠性指标 β9=%3.6f\n',pp(9)); fprintf(1,' 法兰的厚度 hr9=%3.3fmm\n',hr9(4)); y10=[cc4(10) 0 0 0 cc2(10) 0 0 0 cc0(10)]; hr10=roots(y10); fprintf(1,'可靠度R=0.999999 时的可靠性指标 β10=%3.6f\n',pp(10)); fprintf(1,' 法兰的厚度 hr10=%3.3fmm\n',hr10(4)); dp=[hr1(4)*hr2(4)*hr3(4)*hr4(4)*hr5(4)*hr6(4)*hr7(4)*hr8(4)*hr9(4)*hr10(4)]; figure(2); plot(dp,RR) title('\bf 法兰的厚度h与可靠度R的关系曲线\it') xlabel('法兰的厚度h/mm') ylabel('可靠度R') grid; （图片1、2是我参照的。图片3、4是我要验证的，输入可靠度0.9999和答案对不上）文本程序是我改编的不知道哪里有问题？  1.jpg  2.jpg  3.jpg  4.jpg |
回复此楼» 猜你喜欢
限项规定
已经有7人回复
-
青椒八年已不青，大家都被折磨成啥样了？
已经有4人回复
-
救命帖
已经有8人回复
-
招博士
已经有4人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
-
英文综述是否需要润色及查重
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
最失望的一年
已经有18人回复
