matlab叶片强度计算程序.pdfVIP

下载本文档

31
0
约 3页
2015-09-11 发布于重庆
举报
版权申诉

matlab叶片强度计算程序.pdf

1、本文档共3页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

matlab叶片强度计算程序,潮流计算matlab程序,matlab能带计算程序,matlab计算定积分程序,matlab叶片识别,matlab并行计算,精通matlab最优化计算,精通matlab金融计算,matlab数值计算,matlab计算相关系数

%叶片强度，离心拉伸应力，气动弯矩，离心力弯矩，合成弯矩，ABC 三点的弯曲应力和总应力 midu=8.2*10^3; %密度 zhuansu=4700*pi/30; %转速rad/s w=midu*zhuansu*zhuansu; X=[0.0053;0.0041;0.0041;0.0040;0.0024;0.0012]; Y=[-0.0041;-0.0038;-0.0030;-0.0019;-0.0011;-0.0002]; Z=[0.628;0.591;0.56;0.53;0.494;0.458]; A=[0.00018;0.000232;0.000312;0.00041;0.000548;0.000705]; XX=X; YY=Y; ZZ=Z; AA=A; X(1)=[]; Y(1)=[]; Z(1)=[]; A(1)=[]; %后五行 XX(6)=[]; YY(6)=[]; ZZ(6)=[]; AA(6)=[]; %前五行 Xm=(XX+X)/2; % 中间插值 Ym=(YY+Y)/2; % 中间插值 Zm=(ZZ+Z)/2; % 中间插值 Am=(AA+A)/2; % 中间插值 dZ=ZZ-Z; %逐差 %离心拉伸应力 dPc=w.*Am.*Zm.*dZ; dFc=0; Fc=zeros(5,1); Fc(1)= dPc(1); Fc(2)= dPc(1) +dPc(2); Fc(3)= dPc(1) +dPc(2) +dPc(3); Fc(4)= dPc(1) +dPc(2) +dPc(3) +dPc(4); Fc(5)= dPc(1) +dPc(2) +dPc(3) +dPc(4) +dPc(5); Lixinlashenyingli=zeros(5,1); Lixinlashenyingli(1)=Fc(1)/Am( 1); Lixinlashenyingli(2)=Fc(2)/Am(2); Lixinlashenyingli(3)=Fc(3)/Am(3); Lixinlashenyingli(4)=Fc(4)/Am(4); Lixinlashenyingli(5)=Fc(5)/Am(5); Lixinlashenyingli %计算出离心拉伸应力 figure; plot(Z,Lixinlashenyingli,k-*); %求出离心拉伸应力Lixinlashenyingli %下面求气动力弯矩 c1a=297;c1u=-410;midu1=0.894;p1=0.222*10^6; c2a=313;c2u=38;midu2=0.75;p2=0.178*10^6; Q=68; %输入数据 kx=2*pi/Q*(midu1*c1a^2-midu2*c2a^2+p1-p2); ky=2*pi/Q*(midu1*c1a*c1u-midu2*c2a*c2u);%计算数值积分系数 Pxa=kx.*dZ.*Zm; Pya=ky.*dZ.*Zm; Mxa=zeros(5,1);Mya=zeros(5,1); Mxa(1)= -Pya(1)*(Zm(1)-Z(1)); Mxa(2)= -Pya(1)*(Zm(1)-Z(2)) -Pya(2)*(Zm(2)-Z(2)); Mxa(3)= -Pya(1)*(Zm(1)-Z(3)) -Pya(2)*(Zm(2)-Z(3))- Pya(3)*(Zm(3)-Z(3)); Mxa(4)= -Pya(1)*(Zm(1)-Z(4)) -Pya(2)*(Zm(2)-Z(4))- Pya(3)*(Zm(3)-Z(4)) -Pya(4)*(Zm(4)-Z(4)); Mxa(5)= -Pya(1)*(Zm(1)-Z(5)) -Pya(2)*(Zm(2)-Z(5)) -Pya(3)*(Zm(3)-Z(5)) -Pya(4)*(Zm(4)-Z(5)) - Pya(5)*(Zm(5)-Z(5)); Mya(1)= Pxa(1)*(Zm(1)-Z(1)); Mya(2)= Pxa(1)*(Zm(1)-Z(2)) +Pxa(2)*(Zm(2)-Z(2)); Mya(3)= Pxa(1)*(Zm(1)-Z(3)) +Pxa(2)*(Zm(2)-Z(3))+Pxa(3)*(Zm(3)-Z(3)); Mya(4)= Pxa(1)*(Zm(1)-Z(4)) +Pxa(2)*(Zm(2)-Z(4))+Pxa(3)*(Zm(3)-Z(4)) +Pxa(4)*(Zm(4)-Z(4)); Mya(5)=Pxa(1)*(Zm(1)-Z(5))+Pxa(2)