下面的代码实现了。4-1进气过程:V4,P4)到(V1,P1),此时P恒定不变,V逐渐增大;1-2压缩过程(V1,P1)到(V2,P2),服从过程方程,此时P增大,V减小;2-3排气过程从(V2,P2)到(V3,P3),此时P恒定不变,V逐渐减小;3-4膨胀过程从(V3,P3)到(V4,P4),服从过程方程,此时P减小,V增大至原点(V4,P4)。
format long;
% 给定参数
p_in = 4.937; % 进气压力,Mpa
T1 = 0 + 273; % 进气温度,K
p_out = 2.087; % 排气压力,Mpa
T2 = 59.92 + 273; % 排气温度,K
eta_v = 0.7703; % 容积效率
c = 0.2; % 头端余隙
m = 1.277; % 膨胀指数
lambda_v = 0.7861; % 容积系数
lambda_p = 1.0849; % 压力系数
lambda_T = 0.94695; % 温度系数
lambda_l = 0.953782; % 泄漏系数
A_h = 1; % 假设活塞面积为 1(需根据实际情况修改)
A_hg = 0.9; % 假设活塞杆面积为 0.9(需根据实际情况修改),这里是求盖侧的,这个变量