Torsional Analysis Holzer
Uploaded by
arsharmTorsional Analysis Holzer
Uploaded by
arsharmclear % Program for doing Torsional Analysis by Holzer's Method % Written By: Arun Sharma 04/25/07 C=1; M=1;
while C<=99 % Problem defined for a N rotor system N=input('Please enter number of rotors'); for i=1:N J(i)=input('Enter Polar MI for rotors 1, 2 and 3 respectively in lb-in.sec^2'); % Polar moment of inertia of rotor in lb-in.sec^2 end for j=1:N-1 K(j)=input('Enter Stiffness for shafts 1 and 2 respectively in lb-in/rad'); end % Choose starting natural frequency omega omega1=input('Please enter assumed value for starting natural frequency in rad/sec'); theta=1; % Assume angular displacement of rotor at one end % Calculate Jomega square and inertia torque jomegasq= J.*omega1^2; theta=1; Ti=jomegasq.*theta; text=zeros(1,N); text(1,1)=Ti(1,1); for i=2:N theta=theta-(text(1,i-1)/K(i-1)); a1=jomegasq(1,i); text(1,i)=text(1,i-1)+a1*theta; end ExternalTorque(M) = text(1,N) omega(M)=omega1 C=input('If u want to choose another omega enter 1 else enter 100'); M=M+1; end % Now printing the T vs Omega plot figure(1); Te=ExternalTorque; plot(omega,Te);xlabel('Natural Frequency \omega');ylabel('External Torque T'); Title('Torsional Analysis by Holzers Method');
