thelix_a.m

利用matalab编制的转子动力学计算软件。其中模型利用有限元建模

function [psi,ang_at_node,rx,ry]=thelix(R,v,psi)
%  [psi]=thelix(R,r)
%
% Transform a response vector r into helical co-ordinatr system
% such that one of the axes coincides some given angles
%
% Assume that the rotor is undaped and that the co-ordinates system
% is as described in C. W. Lee
%


if isfield(R,'SPRINGS'),
   SPRINGS=R.SPRINGS;  [nd md]=size(SPRINGS); 
else
   SPRINGS=[]; nd=0;
end
ang_node=[];
for q=1:length(psi.angle)
      ang_node=[ang_node;psi.node{q}];
end % for q

if nargin>1, % r provided
   
   Bcord=R.Bcord; nodes=R.NODES;
   
   if length(Bcord)>0,
      N=length(v)+length(Bcord);
      v2=zeros(N,1);
      z=1:N; z(Bcord)=[];  
      v2(z)=v;
   else,
      v2=v(:);
   end
   N=length(v2);
   vy=v2(1:N/2); vz=v2(N/2+1:N);
   L=diff(nodes); L=L(:); 
   Nnodes=length(nodes); Nelem=Nnodes-1;
   
   na=length(psi.angle);
   
   
   b0=psi.angle{1};  % This is the initial angle
%   b1=psi.angle{2};
   ang_at_node=ones(size(nodes))*b0;
   for q=1:na-1, 
      i1=ang_node(q);  i2=ang_node(q+1); 
      b0=psi.angle{q};  % This is the initial angle
      b1=psi.angle{q+1};

      ang_at_node(i1:i2)=([i1:i2]-i1)*(b1-b0)/(i2-i1)+b0;
    end
    
    if i2<Nnodes
         ang_at_node(i2:Nnodes)=ones(size(ang_at_node(i2:Nnodes)))* ang_at_node(i2);
      end
      
      % ====> Transform r into helical co-ordinates
      rx=zeros(Nnodes,1); ry=zeros(Nnodes,1);
      vz=-i*vz;
      for q=1:Nnodes
         a=-ang_at_node(q); T=[cos(a) sin(a); -sin(a) cos(a)];
         x=vy(2*q-1);    y=vz(2*q-1);
         xya=T*[x;y];  rx(q)=xya(1); ry(q)=xya(2);
      end
      
end