原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
[航空航天]故障诊断作业
发动机状态监控与故障诊断大作业
第一题：发动机正常态数学模型
已知参数：
=1.40 =1.33 =1004.71 =1158.31
=0.99 =0.9500 =0.9800
=0.26397 =288.15 =101325 K=0.0397
X0=; X1=; X2=; X3=; X4=; X5=; X6=; X7=
源代码： 思路：
q1=sym(0.001116*x1^2-0.03825*x1+x0=0.5617); %压气机效率特性，控制条件压气机转数n=常数
q2=sym(0.3611*x1^2-4.500*x1+x6=98.00); %压气机流量特性
q3=sym(x3-288.15*(x1^ ((1.4-1)/1.4)-1)/x0=288.15); %压气机进出口截面面积参数关系
q4=sym(x6-(0.0397*0.9500*101325*0.072325*x1)/sqrt(x4)=0); %压气机与涡轮之间的连续方程
q5=sym(1004.71*(x3-288.15)-1158.31*(x4-x5)*0.99=0); %压气机与涡轮之间的功率平衡
q6=sym(x4-x5-0.9100*x4*(1-x2^(-(1.33-1)/1.33))=0); %涡轮进出口截面参数关系
q7=sym(1081.55*x6*(x4-x3)-0.9800x7=0); %燃烧室能量平衡
q8=sym(0.072325*x2*sqrt(x5)-0.26397*sqrt(x4)=0); %涡轮与尾喷管之间的连续方程
[x0,x1,x2,x3,x4,x5,x6,x7]=solve(q1,q2,q3,q4,q5,q6,q7,q8)% solve解方程组
计算结果：
x0 =0.85999701373329763609296449279423
x1 =12.000088412874692615393372511433
x2 =4.2892009508159647410202861471357
x3 =634.58164842164877933153995078615
x4 =1100.0097318689014719523661532824
x5 =796.48206413413851742309549262696
x6 =100.00123163377622105906413930964
x7 =1.1945655991850663939253290226252
第二题：假设压气机发生故障
源代码 思路：
q1=sym(0.01147*x1-x0+0.7124=0); %压气机故障实际特性方程,取控制条件n=常数，n为压气机转速
q2=sym(0.3611*x1^2-4.500*x1+x6=98.00); %压气机流量特性
q3=sym(x3-288.15*(x1((1.4-1)/1.4)-1)/x0=288.15); %压气机进出口截面参数关系
q4=sym(x6-(0.0397*0.9500*101325*0.072325*x1)/sqrt(x4)=0);%压气机与涡轮之间连续方程
q5=sym(1004.71*(x3-288.15)-1158.31*(x4-x5)*0.99=0); %压气机与涡轮之间功率平衡
q6=sym(x4-x5-0.9100*x4*(1-x2^(-(1.33-1)/1.33))=0);%涡轮进出口截面参数关系
q7=sym(1081.55*x6*(x4-x3)-0.9800x7=0);%燃烧室能量平衡
q8=sym(0.072325*x2*sqrt(x5)-0.26397*sqrt(x4)=0);%涡轮与尾喷管之间连续方程
[x0,x1,x2,x3,x4,x5,x6,x7]=solve(q1,q2,q3,q4,q5,q6,q7,q8);%运用solve解方程组并赋值给x0~x7
计算结果：
x0 =0.85065681543814275878270079735312
x1 =12.053776411346360835457785296697
x2 =4.2892009508159647410202861471357
x3 =639.26475751291967142
文档评论（0）