qiang.cpp

北京航空航天大学能源与动力工程学院 强度与振动第一题大作业

#include<math.h>
#include<iostream.h>

double x[6]={0.0053,0.0041,0.0041,0.0040,0.0024,0.0012};
double y[6]={-0.0041,-0.0038,-0.003,-0.0019,-0.0011,-0.0002};
double z[6]={0.628,0.591,0.56,0.53,0.494,0.458};
double A[6]={1.8*pow(10,-4),2.32*pow(10,-4),3.12*pow(10,-4),4.1*pow(10,-4),5.48*pow(10,-4),7.05*pow(10,-4)};
double st[6],Am[6],zm[6],dz[6],St,Dct[6],Px[6],Py[6],PX,PY,Mxq[6],Myq[6],ym[6],Pz[6],dV[6],xm[6],PZ,Mxl[6],Myl[6];  
double c1am,c1um,d1m,p1m,c2am,c2um,d2m,p2m,Mx[6],My[6],m,n,Mvs[6],Mut[6];
double a[6]={31.667*0.0174,27.817*0.0174,25.317*0.0174,22.092*0.0174,16.95*0.0174,12.717*0.0174};	
double Da[6],Db[6],Dc[6],DA[6],DB[6],DC[6];
double Iu[6]={0,6.694*pow(10,-8),9.332*pow(10,-8),12.52*pow(10,-8),17.57*pow(10,-8),23.74*pow(10,-8)};
double Iv[6]={0,0.242*pow(10,-8),0.304*pow(10,-8),0.484*pow(10,-8),0.939*pow(10,-8),1.802*pow(10,-8)};
double Ua[6]={0,-0.02685,-0.02847,-0.02938,-0.02889,-0.02894},Va[6]={0,0.00797,0.00951,0.01094,0.01232,0.01319};	
double Ub[6]={0,-0.00084,-0.00205,-0.00303,-0.00219,-0.00302},Vb[6]={0,-0.00481,-0.00521,-0.00655,-0.00749,-0.01015};	
double Uc[6]={0,0.03728,0.03909,0.0406,0.04366,0.04597},Vc[6]={0,0.00773,0.00824,0.00840,0.0113,0.01305};
double pi;

void main()
{	
	int i,j;
	pi=3.14159,St=0,PX=0,PY=0;
	c1am=297,c1um=-410,d1m=0.894,p1m=0.222*pow(10,6);
	c2am=313,c2um=38,d2m=0.75,p2m=0.178*pow(10,6);
	for(i=1;i<6;i++)
	{		
		Am[i]=0.5*(A[i-1]+A[i]),dz[i]=z[i-1]-z[i],zm[i]=0.5*(z[i-1]+z[i]),
		st[i]=8200*pow((2*3.14159*4700/60),2)*Am[i]*zm[i]*dz[i],St=St+st[i],Dct[i]=St/A[i];   
		cout<<"离心应力Pc:";
		cout<<Dct[i]<<'\n';
	}
	cout<<'\n';

	for(i=1;i<6;i++)
	{
		zm[i]=0.5*(z[i-1]+z[i]);
		Px[i]=2*pi*zm[i]/68*(d1m*pow(c1am,2)-d2m*pow(c2am,2)+p1m-p2m)*(z[i-1]-z[i]);
		Py[i]=2*pi*zm[i]/68*(d1m*c1am*c1um-d2m*c2am*c2um)*(z[i-1]-z[i]);
	}
	
	for(j=1;j<6;j++)
	{		
		PX=0,PY=0;		
		for(i=1;i<=j;i++)
		{	
			PX=PX+Px[i]*(zm[i]-z[j]),PY=PY+Py[i]*(zm[i]-z[j]);
			}		
		Mxq[j]=-PY,Myq[j]=PX;
		cout<<"气动力弯矩:";
		cout<<"Mx:"<<Mxq[j]<<";My:"<<Myq[j]<<'\n';
	}
	cout<<'\n';

	PY=0,PZ=0,PX=0;  
	for(i=1;i<6;i++)		
	{
		dV[i]=0.5*(A[i-1]+A[i])*(z[i-1]-z[i]);
		xm[i]=0.5*(x[i-1]+x[i]),ym[i]=0.5*(y[i-1]+y[i]),zm[i]=0.5*(z[i-1]+z[i]);
		Py[i]=8200*pow((2*pi*4700/60),2)*dV[i]*ym[i],Pz[i]=8200*pow((2*pi*4700/60),2)*dV[i]*zm[i];
	}
	for(j=1;j<6;j++)
	{		
		PY=0,PZ=0,PX=0;
		for(i=1;i<=j;i++)
		{
			PY=PY+Py[i]*(zm[i]-z[j]),PZ=PZ+Pz[i]*(ym[i]-y[j]),PX=PX+Pz[i]*(xm[i]-x[j]);
		}
		Mxl[j]=PZ-PY,Myl[j]=-PX;
	}
	for(i=1;i<6;i++)
	{
		cout<<"离心弯矩:";
		cout<<"Mx1:"<<Mxl[i]<<";My1:"<<Myl[i]<<'\n';
		Mx[i]=Mxl[i]+Mxq[i],My[i]=Myl[i]+Myq[i];
		cout<<"合成弯矩:";
		cout<<"Mx:"<<Mx[i]<<";My:"<<My[i]<<'\n';
	}
	cout<<'\n';

	for(i=1;i<6;i++)
	{   
		Mvs[i]=Mx[i]*cos(a[i])+My[i]*sin(a[i]),Mut[i]=-Mx[i]*sin(a[i])+My[i]*cos(a[i]);
		Da[i]=Mvs[i]/Iv[i]*Va[i]-Mut[i]/Iu[i]*Ua[i];
		Db[i]=Mvs[i]/Iv[i]*Vb[i]-Mut[i]/Iu[i]*Ub[i];
		Dc[i]=Mvs[i]/Iv[i]*Vc[i]-Mut[i]/Iu[i]*Uc[i];
	}
	for(i=1;i<6;i++)
	{
		cout<<"ABC弯曲应力:";
		cout<<"Da:"<<Da[i]<<";Db:"<<Db[i]<<";Dc:"<<Dc[i]<<'\n';
		DA[i]=Dct[i]+Da[i],DB[i]=Dct[i]+Db[i],DC[i]=Dct[i]+Dc[i];
		cout<<"ABC总应力:";
		cout<<"DA:"<<DA[i]<<";DB:"<<DB[i]<<";DC:"<<DC[i]<<'\n';
	}
    cout<<'\n';
}