matlab叶片强度计算程序

《matlab叶片强度计算程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、%叶片强度，离心拉伸应力，气动弯矩，离心力弯矩，合成弯矩，ABC三点的弯曲应力和总应力midu=8.2*10^3;%密度zhuansu=4700*pi/30;%转速rad/sw=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.0002

2、32;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

3、)=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);Lixi

4、nlashenyingli(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

5、/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

6、(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)*(

7、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)*(Zm(2)-Z(5))+Pxa(3)*(Zm(3)-Z(5))+Pxa(4)*(Zm

8、(4)-Z(5))+Pxa(5)*(Zm(5)-Z(5));MxaMyafigure;plotyy(Z,Mxa,Z,Mya);%求出气动弯矩%下面求离心弯矩dV=Am.*dZ/2;dFy=w*dV.*Ym;dFz=w*dV.*Zm;Mxc=zeros(5,1);Myc=zeros(5,1);forj=1:1:5%1,