How to rotating a stl object

Question

TOM SALIC 2022-9-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1807310-how-to-rotating-a-stl-object

评论： TOM SALIC 2022-9-19

采纳的回答： Fabio Freschi

Rotor.zip

在 MATLAB Online 中打开

Hi everyone,

I know this question has already been asked but here I am trying to spin the rotor of a wind turbine. I have the stl file and I applied the rotation matrix but the object is rotating but also in translation. My objective is simply to make it turn on the spot. Thanks for your advices.

[p,t,tnorm]=import_stl_fast('Rotor.stl',1,2);
TR=triangulation(t,p);
%trimesh(TR)
 
   ax=linspace(0,2*pi,10);
  
   az=pi/4;
   ay=pi/4;
   
   Ry=[cos(ay) 0 sin(ay); 0 1 0; -sin(ay) 0 cos(ay)];
   Rz=[cos(az) -sin(az) 0; sin(az) cos(az) 0;  0 0 1];
   filename='test.gif';
   figure;
  
    h1=trisurf(t,p(:,1),p(:,2),p(:,3), ...
       'FaceColor','black','EdgeColor','none','LineWidth',0.25);
   
for idx=1:length(ax)
   Rx=[ 1 0 0; 0 cos(idx) -sin(idx); 0 sin(idx) cos(idx)];
   PointR=p*Rx;
   
   
  
   h1=trisurf(t,PointR(:,1),PointR(:,2),PointR(:,3), ...
       'FaceColor','black','EdgeColor','none','LineWidth',0.25);
   %view(90,0)
        frame = getframe(1);
      im = frame2im(frame);
      [imind,cm] = rgb2ind(im,256);
      if idx == 1
          imwrite(imind,cm,filename,'gif', 'Loopcount',inf);
      else
          imwrite(imind,cm,filename,'gif','WriteMode','append');
      end
      pause(0.1)
end

I have attached the stl file

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

Jan 2022-9-18

If you explain, what "doesn't work" means, it is possible to help you. Do you get an error message? Does the result differ from your expectations?

TOM SALIC 2022-9-19

编辑：TOM SALIC 2022-9-19

在 MATLAB Online 中打开

test.gif

Sorry.

No I have no error but the result is the same as before. I specified the axis of rotation in the "AxelRot" function. But the stl object keeps moving off its axis of rotation.

 ax=linspace(0,2*pi,10);
 [p,t,tnorm]=import_stl_fast('Rotor.stl',1,2);
TR=triangulation(t,p);
filename='test.gif';
figure;
% x=[-13.30 100]; y=[0 0]; z=[89.35 89.35];
% line(x,y,z);
for idx=1:length(ax)
[PointR, R, t1] = AxelRot(p(:,:)',ax(idx)*(180/pi), [1 0 0],[-13.30 0 89.35]);
%  bhold on;
h1=trisurf(t,PointR(1,:),PointR(2,:),PointR(3,:), ...
       'FaceColor','black','EdgeColor','none','LineWidth',0.25);
   view(90,180)
        frame = getframe(1);
      im = frame2im(frame);
      [imind,cm] = rgb2ind(im,256);
      if idx == 1
          imwrite(imind,cm,filename,'gif', 'Loopcount',inf);
      else
          imwrite(imind,cm,filename,'gif','WriteMode','append');
      end
      pause(0.1)
end