内燃机设计 程序2.0

源程序

#include <cstdio>
#include <math.h>
#include<iostream>
using namespace std;
#define PI 3.14159
int main()
{
   
double a[72]={
   0.0857,0.0894,0.0833,0.0838,0.0822,0.0818,0.0841,0.0845,
0.0834,0.0858,0.0862,0.0855,0.0892,0.0919,0.0903,0.0929,
0.0947,0.0928,0.0980,0.0976,0.0968,
0.1033,0.1073,0.1167,0.1308,0.1552,0.1783,
0.2191,0.2740,0.3596,0.4990,0.7298,
1.1332,1.8783,3.1401,4.6997,6.7025,7.0849
,5.1156,3.5106,2.4167,1.6797,1.2474,
0.9562,0.7655,0.6405,0.5507,0.4842,0.4413
,0.4053,0.3800,0.3659,0.3423,0.3150,
0.2093,0.1557,0.1379,0.1143,0.0948,0.0821
,0.0915,0.0845,0.0888,0.0974,0.1031,
0.0946,0.0894,0.0842,0.0854,0.0930,0.0967,0.0936};//定义数组a
int i;
double b,pg1;
    double zj=0;
for(i=0;i<=71;i++)
{
   
pg1=a[i];
printf("曲柄转角(度)b=%5.2f\n",zj);
printf("缸内气体压力(MPa)pg1=%5.4f\n",pg1);
b=zj*PI/180;
zj+=10;
double function(double z);
double ma=1.38,mb=1.267,m1=0.543,mk1=0.563,mk2=1.15,p2=0.0415,
mk3=2.3,p3=0.0688,r=0.0575,l=0.21,n=2000,e=0.007;
double m2,h,s,w,B,w1,q1,q2,mk;
m2=m1+mb;//连杆质量
h=e/r;//偏心率
s=r/l;//连杆比
w=PI*n/30;//曲柄角速度
B=asin(s*(sin(b)-h));//连杆摆动角位移
w1=s*cos(b)/sqrt(1-s*s*(sin(b)-h)*(sin(b)-h));//连杆角速度
mk=mk1+2*p2*mk2/r-2*p3*mk3/r;
printf("连杆摆动角位移B=%5.2f\n",B);
printf("连杆摆动角速度w1=%5.2f\n",w1);
double x,v,g;
x=r*((1-cos(b))+0.5*s*sin(b)*sin(b)-s*h*sin(b));
v=r*w*(sin(b)+0.5*s*sin(2*b)-s*h*cos(b));
g=r*w*w*(cos(b)+s*cos(2*b)+s*h*sin(b));
printf("活塞位移x=%5.4f\n",x);
printf("活塞速度v=%5.4f\n",v);
printf("活塞加速度g=%5.4f\n",g);
double mj,A,P0,Pg,Pi,P,p,Pn,Pt,t,k;
mj=ma+m1;//往复运动质量
int D=95;
A=PI*D*D/4;
P0=0.101;
Pg=pg1-P0;
Pi=(-1)*mj*g/A;
P=Pg+Pi;
Pn=P*tan(B);
Pt=P/cos(B);
t=Pt*sin(b+B);
k=Pt*cos(b+B);
printf("缸内气压引起的作用力Pg=%5.4f\n",Pg);
printf("往复惯性力Pi=%5.4f\n",Pi);
printf("活塞上的总作用力P=%5.4f\n",P);
printf("活塞侧压力Pn=%5.4f\n",Pn);
printf("连杆力Pt=%5.4f\n",Pt);
printf("切向力t=%5.4f\n",t);
printf("径向力k=%5.4f\n",k);
double kr1,Pqx,Pqy,Pq,Pqx1,Pqy1,aq1,aq;
kr1=mb*r*w*w/A;
Pqx=kr1-k;
Pqy=t;
Pq=sqrt(Pqx*Pqx+Pqy*Pqy);
Pqx1=function(Pqx);
Pqy1=function(Pqy);
aq1=atan(fabs(Pqy/Pqx));
aq=(1-Pqx1)*90+Pqx1*Pqy1*aq1*180/PI;
printf("连杆轴颈正反切角度及负荷:aq=%5.4f,Pq=%5.4f\n",aq,Pq);
double Pp,ap;
Pp=Pq;//连杆轴承负荷
ap=aq+180+b*180/PI+B*180/PI;
printf("连杆轴承正反切角度及负荷:aq=%5.4f,Pp=%5.4f\n",ap,Pp);
double Pc,ac;
double krk,Pzx,Pzy,Pz,Pzx1,Pzy1,az1,az;
krk=mk*r*w*w/A;
Pzx=-0.5*(Pqx+krk);
Pzy=-0.5*t;
Pz=sqrt(Pzx*Pzx+Pzy*Pzy);
Pzx1=function(Pzx);
Pzy1=function(Pzy);
az1=atan(fabs(Pzy/Pzx));
az=(1-Pzx1)*90+Pzx1*Pzy1*az1*180/PI;
printf("第一主轴颈角度及负荷:az=%5.4f,Pz=%5.4f\n",az,Pz);
Pc=Pz;
ac=az+180+b*180/PI;
printf("第一主轴承角度及负荷:ac=%5.4f,Pc=%5.4f\n",ac,Pc);
printf("\n");
}
return 0;
}
double function(double z)//定义符号函数
{
   
double N,M;
N=fabs(z);
M=z/N;
return M;
}

结果如下

曲柄转角(度)b= 0.00
缸内气体压力(MPa)pg1=0.0857
连杆摆动角位移B=-0.03
连杆摆动角速度w1= 0.27
活塞位移x=0.0000
活塞速度v=-0.4014
活塞加速度g=3212.8380
缸内气压引起的作用力Pg=-0.0153
往复惯性力Pi=-0.8716
活塞上的总作用力P=-0.8869
活塞侧压力Pn=0.0296
连杆力Pt=-0.8874
切向力t=0.0296
径向力k=-0.8869
连杆轴颈正反切角度及负荷:aq=1.2667,Pq=1.3381
连杆轴承正反切角度及负荷:aq=179.3565,Pp=1.3381
第一主轴颈角度及负荷:az=189.8551,Pz=0.0864
第一主轴承角度及负荷:ac=369.8551,Pc=0.0864

曲柄转角(度)b=10.00
缸内气体压力(MPa)pg1=0.0894
连杆摆动角位移B= 0.01
连杆摆动角速度w1= 0.27
活塞位移x=0.0008
活塞速度v=2.2598
活塞加速度g=3147.4703
缸内气压引起的作用力Pg=-0.0116
往复惯性力Pi=-0.8539
活塞上的总作用力P=-0.8655
活塞侧压力Pn=-0.0123
连杆力Pt=-0.8656
切向力t=-0.1624
径向力k=-0.8502
连杆轴颈正反切角度及负荷:aq=-7.1153,Pq=1.3111
连杆轴承正反切角度及负荷:aq=183.6991,Pp=1.3111
第一主轴颈角度及负荷:az=129.4329,Pz=0.1051
第一主轴承角度及负荷:ac=319.4329,Pc=0.1051

曲柄转角(度)b=20.00
缸内气体压力(MPa)pg1=0.0833
连杆摆动角位移B= 0.06
连杆摆动角速度w1= 0.26
活塞位移x=0.0037
活塞速度v=4.8014
活塞加速度g=2927.9126
缸内气压引起的作用力Pg=-0.0177
往复惯性力Pi=-0.7943
活塞上的总作用力P=-0.8120
活塞侧压力Pn=-0.0491
连杆力Pt=-0.8135
切向力t=-0.3238
径向力k=-0.7463
连杆轴颈正反切角度及负荷:aq=-15.1371,Pq=1.2401
连杆轴承正反切角度及负荷:aq=188.3208,Pp=1.2401
第一主轴颈角度及负荷:az=95.2271,Pz=0.1626
第一主轴承角度及负荷:ac=295.2271,Pc=0.1626

曲柄转角(度)b=30.00
缸内气体压力(MPa)pg1=0.0838
连杆摆动角位移B= 0.10
连杆摆动角速度w1= 0.24
活塞位移x=0.0087
活塞速度v=7.1016
活塞加速度g=2571.6567
缸内气压引起的作用力Pg=-0.0172
往复惯性力Pi=-0.6977
活塞上的总作用力P=-0.7149
活塞侧压力Pn=-0.0744
连杆力Pt=-0.7187
切向力t=-0.4219
径向力k=-0.5819
连杆轴颈正反切角度及负荷:aq=-22.2219,Pq=1.1156
连杆轴承正反切角度及负荷:aq=193.7230,Pp=1.1156
第一主轴颈角度及负荷:az=72.2853,Pz=0.2215
第一主轴承角度及负荷:ac=282.2853,Pc=0.2215

曲柄转角(度)b=40.00
缸内气体压力(MPa)pg1=0.0822
连杆摆动角位移B= 0.14
连杆摆动角速度w1= 0.21
活塞位移x=0.0155
活塞速度v=9.0571
活塞加速度g=2106.1055
缸内气压引起的作用力Pg=-0.0188
往复惯性力Pi=-0.5714
活塞上的总作用力P=-0.5902
活塞侧压力Pn=-0.0851
连杆力Pt=