Ploting temperature evolution in a multi-layer tyre

2023 10 24
1 个回答
回答已采纳
更新时间：2024 2 8
6 次查看（30 天）

0 个投票

Hello,
I'm looking for a graphical representation of temperature changes in a tyre. My simulink gives me the real time temperatures of my 300 layers in a spreadsheet called "temperatureData" but I'm now trying to represent the 300 layers of the tyre according to their temperature using a colour scale as shown in the image above. I've tried several methods but I can't get anywhere. If anyone has an idea or a code proposal that might be suitable I'd be happy to discuss it and see what they come up with !
I'm also copying an excel of the temperature as an example (a shorter version).
Thank you in advance for your help.

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

 采纳的回答

Mathieu NOE
Mathieu NOE 2023-10-24

2 个投票

maybe this ?
I assumed your data has columns = layers and the rows spans the angular position between 0 and 2pi
as a result , a nice donut but does it represents really a tire?
A =readmatrix('temperatureData.xlsx');
[m,layers] = size(A); % rows = angular position (0 to 2pi), columns = layers
da = 2*pi/m;
angl = (0:da:2*pi-da)'*ones(1,layers);
% tire defined by inner radius ri and outer radius ro (units ??)
ri = 100;
ro = 200;
r = ones(m,1)*linspace(ri,ro,layers);
% create X,Y meshes
X = r.*cos(angl);
Y = r.*sin(angl);
% plot
surf(X,Y,A);
shading interp
view(2)
axis square

4 个评论
显示 2更早的评论 隐藏 2更早的评论

Jan Soller
Jan Soller 2023-10-24
Hello Mathieu,
First of all, thanks a lot for your answer, this donut looks just awesome !! :D
I'm sorry if I didn't explain the data properly. The columns in Excel represent the 12 layers of the tyre, you'r right. However, the rows are the temperature values over time. This means that line 235, for example, gives the temperatures of the 12 tyre layers after 235 seconds.
I imagine that this would change a large part of the programme but the result that I imagined is at least in this form.
Would it be possible to have a programme that shows the temperatures of each layer as the simulation progresses (as time goes by)?
Thanks again
Mathieu NOE
Mathieu NOE 2023-10-30
编辑：Mathieu NOE 2024-2-8
hello Jan
try this, it will generate a MP4 video file (demo.mp4)
as your time axis is very long the video file generation can take a while, be patient !
A =readmatrix('temperatureData.xlsx');
[time,layers] = size(A); % rows = time stamps (seconds) , columns = layers
% tire defined by inner radius ri and outer radius ro (units ??)
ri = 100;
ro = 200;
nb = 100; % number or radii to plot
r = ones(nb,1)*linspace(ri,ro,layers);
da = 2*pi/(nb-1);
angl = (0:da:2*pi)'*ones(1,layers);
% create X,Y mesh
X = r.*cos(angl);
Y = r.*sin(angl);
% fix colorbar min/max values
minA = min(A,[],'all');
maxA = max(A,[],'all');
%==============================================================================================
% Prepare the video file.
vidObj = VideoWriter('demo.mp4','MPEG-4');
open(vidObj);
numberOfFrames = time; % maybe consider decimation / resampling to speed up the video process
for ci = 1:numberOfFramesa
cla reset;
s = surf(X,Y,ones(nb,1)*A(ci,:));
s.EdgeColor = 'none'; % hide the edges
title(['Time index = ' num2str(ci), ' seconds '] , 'FontSize', 15);
colorbar('vert');
caxis([minA-5, maxA+5]);
shading interp
view(2)
axis square
drawnow;
% Write each frame to the file.
currFrame = getframe(gcf);
writeVideo(vidObj,currFrame);
end
% Close the video file.
close(vidObj);
Jan Soller
Jan Soller 2023-10-30
编辑：Jan Soller 2024-2-6
Hello Mathieu,
Thanks a lot for your help, it is exactly what I was searching for and the result looks really impressive!
Adding the video recording is also just perfect, as I can accelerate it in post-processing.
Thanks again and I wish you all the best,

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Foundation and Custom Domains 的更多信息

产品

版本

R2022b

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by