《叶片受力分析.doc》由会员分享，可在线阅读，更多相关《叶片受力分析.doc（6页珍藏版）》请在金锄头文库上搜索。
1、叶片截面1截面2截面3截面4截面5xm0.00470.00410.004050.00320.0018ym0.003950.00340.002450.00150.00065Am0.0002060.0002720.0003610.0004790.0006265xh0.001200.00010.00160.0012yh0.00030.00080.00110.00080.0009zh0.0370.0310.030.0360.036Pz9228.049639.2411724.417537.821325.4Py59.804456.947752.706251.380229.1208拉257777374205
2、413313402990354552Fx11.30668.944698.19749.241258.59148Fy245.431194.16177.939200.598186.493Mx_gas-4.540461.5309812.489227.404353.5My_gas0.209173-0.0705302-0.575361-1.26248-2.46467Mx_离0.2778252.69518-3.12148-25.5983-47.0795My_离-5.536825.536824.9506-6.0204658.1915Mx-4.262644.226169.36773 1.806 6.42054M
3、y -5.327655.46629 4.37524-7.2829455.7268M -1.37934 1.520057.046613.74348 -5.30602 M-6.68217 6.7402 7.56587-6.5030155.844弯a -1.86244 2.6988614.15225.98086-10.6914弯b 2.82543 -2.75315 -9.71932 -2.90496 -0.411531 弯c -8.12731 6.9434714.68312.888996.971 （以上个参量均使用国际单位制）使用C+编程，程序如下：#include#include#define N
4、 6#define M 5#define Q 68void average (double a,int n,double b) /平均int i=0;for(i;in-1;i+)bi=(ai+ai+1)/2;void altitude (double a,int n,double b) int i=0;for(i;in;i+)bi=ai-ai+1;void load(double a,double b,double c,double d,double f,double g) /离心力for(int i=0;iM;i+)ci=a*b*b*di*fi*gi;void sentrifugal_for
5、ce(double a,double b,double c) /离心拉应力double d;for(int i=0;iM;i+)d=0;for(int j=0;j=i;j+)d=d+bj; ai=d/ci+1;void main() double pi=3.1415926;double w=4700*2*pi/60; int i,j; double den=8200;double radianN=0;radian0=(30+40/360)*pi/180;radian1=(27+49/360)*pi/180;radian2=(25+19/360)*pi/180;radian3=(22+5/360
6、+30/360/360)*pi/180;radian4=(16+57/360)*pi/180;radian5=(12+43/360)*pi/180; double xN=0.0053,0.0041,0.0041,0.0040,0.0024,0.0012;double yN=0.0041,0.0038,0.0030,0.0019,0.0011,0.0002;double zN=0.628,0.591,0.56,0.53,0.494,0.458;double AN=0.00018,0.000232,0.000312,0.000410,0.000548,0.000705;double c1am=29
7、7; double c1um=410; double den1m=0.894; double p1m=222;double c2am=313; double c2um=-38; double den2m=0.75; double p2m=178;double Jat5=0.242,0.304,0.484,0.939,1.802;double Jkc5=6.694,9.332,12.52,17.57,23.74;double ata5=-2.685,-2.847,-2.938,-2.889,-2.894;double kca5=-0.797,-0.951,-1.094,-1.232,-1.319
8、;double atb5=-0.084,-0.205,-0.303,-0.219,-0.302;double kcb5=0.481,0.521,0.655,0.749,0.1015;double atc5=3.728,3.909,4.060,4.366,4.597;double kcc5=-0.773,-0.824,-0.840,-1.130,-1.305;double xm5=0; average(x,N,xm); coutxm ;for(i=0;iN-1;i+) coutxmit; coutn;double ym5=0; average(y,N,ym); coutym ;for(i=0;i
9、N-1;i+) coutymit; coutn;double zm5=0; average(z,N,zm); coutzm ;for(i=0;iN-1;i+) coutzmit; coutn; double Am5=0; average(A,N,Am); coutAm ;for(i=0;iN-1;i+) coutAmit; coutn; double xh5=0; altitude(x,M,xh);coutxh ;for(i=0;iN-1;i+) coutxhit; coutn;double yh5=0; altitude(y,M,yh);coutyh ;for(i=0;iN-1;i+) co
10、utyhit; coutn;double zh5=0; altitude(z,M,zh);coutzh ;for(i=0;iN-1;i+) coutzhit; coutn;double Pz5=0;load(den,w,Pz,Am,zm,zh);coutPz ;for(i=0;iN-1;i+) coutPzit; coutn;double Py5=0;load(den,w,Py,Am,ym,zh);coutPy ;for(i=0;iN-1;i+) coutPyit; coutn; coutn; double stress_sen5=0;sentrifugal_force(stress_sen,
11、Py,A);coutstress_sen ;for(i=0;iN-1;i+) coutstress_senit; coutn; coutn;double Fx5=0; double Fy5=0;double a=den1m*c1am*c1am-den2m*c2am*c2am+(p1m-p2m);double b=den1m*c1am*c1um-den2m*c2am*c2um;for(i=0;iM;i+) Fxi=2*pi*zmi*zhi*a/Q; Fyi=2*pi*zmi*zhi*b/Q; coutFx ;for(i=0;iN-1;i+) coutFxit; coutn;coutFy ;for(i=0;iN-1;i+) coutFyit; coutn; double Mx_gasM; double My_gasM; double Mx_senM; double My_senM; /弯矩double MxM;double MyM; double MkcM; double MatM;double stress_aM; double stress_bM;double stress_cM;double d,e,f,g;for(j=0;jM;j+)d=0;e=0;f=0;g=0; for(i=0;i=j;i+)d=d+Fyi*(z
