calulate the indicated mean pressure with a p-v diagramm

Question

Julian 2024-6-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2130706-calulate-the-indicated-mean-pressure-with-a-p-v-diagramm

编辑： Mathieu NOE 2024-6-28

Hello,

I have to calculate the indicated mean pressure. I have a p-v diagramm (in the attachements). I separated the two areas by finding the intersection. But know im not sure if the polyarea() function is the correct function to calculate the two areas.

The integral i have to implement is also in the attachements

(p: pressure in the cylinder; A: piston surface; V: volume of the cylinder; s: piston travel (path) as a function of the crank angle

Thanks for help!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Mathieu NOE 2024-6-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2130706-calulate-the-indicated-mean-pressure-with-a-p-v-diagramm#answer_1477261

编辑：Mathieu NOE 2024-6-28

在 MATLAB Online 中打开

data.mat

hello Julian

yes you can use polyarea, I believe you could also with trapz like in my example below (it uses this Fex submission : Fast and Robust Self-Intersections - File Exchange - MATLAB Central (mathworks.com))

results (two loops area) :

area = 71.3220 125.0467 (with trapz)

area2 = 72.5056 125.3691 (with polyarea)

code :

load('data.mat') 
[x0,y0,segments]=selfintersect(x,y); % fex : https://fr.mathworks.com/matlabcentral/fileexchange/13351-fast-and-robust-self-intersections
figure(1)
plot(x,y,'b',x0,y0,'dr','markersize',15);
axis square
hold on
% compute area for each loop 
for k = 1:numel(x0)
    ind = (segments(k,1):segments(k,2));
    x_tmp = x(ind);
    y_tmp = y(ind);
    % compute area 
    area(k) = trapz(x_tmp,y_tmp);
    area2(k) = polyarea(x_tmp,y_tmp);
    plot(x_tmp,y_tmp)
end
area
area2