cantbeam_ss_freq.m

MATLAB PROGRAMS FOR "VIBRATION SIMULATION USING MATLAB AND ANSYS" All the M-files which are liste

				 55347.    
				 65384.    
				 76254.    
				 87959.    
				0.10050E+06
				0.11387E+06
				0.12807E+06
				0.14311E+06
				0.15898E+06
				0.17569E+06
				0.19323E+06
				0.21160E+06
				0.23081E+06
				0.25085E+06
				0.27173E+06
				0.29345E+06
				0.31600E+06
				0.33940E+06
				0.36365E+06
				0.38874E+06
				0.41468E+06
				0.44147E+06
				0.46911E+06
				0.49762E+06
				0.52699E+06
				0.55722E+06
				0.58833E+06
				0.62032E+06
				0.65320E+06
				0.68696E+06
				0.72162E+06
				0.75718E+06
				0.79365E+06
				0.83104E+06
				0.86935E+06
				0.90859E+06
				0.94878E+06
				0.98990E+06
				0.10320E+07
				0.10750E+07
				0.11190E+07
				0.11639E+07
				0.12098E+07
				0.12566E+07
				0.13043E+07
				0.13528E+07
				0.14021E+07
				0.14519E+07
				0.15020E+07
				0.15517E+07
				0.15997E+07
				0.16432E+07
				0.16758E+07
				0.18739E+07
				0.19105E+07
				0.19598E+07
				0.20161E+07
				0.20764E+07
				0.21398E+07
				0.22054E+07
				0.22732E+07
				0.23428E+07
				0.24144E+07
				0.24877E+07
				0.25628E+07
				0.26398E+07
				0.27186E+07
				0.27992E+07
				0.28817E+07
				0.29661E+07
				0.30525E+07
				0.31408E+07
				0.32312E+07
				0.33236E+07
				0.34182E+07
				0.35149E+07
				0.36138E+07
				0.37149E+07
				0.38182E+07
				0.39238E+07
				0.40317E+07
				0.41420E+07
				0.42546E+07
				0.43695E+07
				0.44868E+07
				0.46064E+07
				0.47283E+07
				0.48524E+07
				0.49787E+07
				0.51071E+07
				0.52375E+07
				0.53697E+07
				0.55036E+07
				0.56389E+07
				0.57755E+07
				0.59130E+07
				0.60511E+07
				0.61894E+07
				0.63274E+07
				0.64648E+07
				0.66008E+07
				0.67349E+07
				0.68664E+07
				0.69945E+07
				0.71185E+07
				0.72376E+07
				0.73507E+07
				0.74569E+07
				0.75555E+07
				0.76452E+07
				0.77254E+07
				0.77950E+07
				0.78532E+07
				0.78993E+07
				0.79328E+07
				0.79530E+07
				0.88695E+07];


%	64red

	freq64red = [149.20    
				 935.00    
				 2618.0    
				 5130.0    
				 8479.8    
				 12667.    
				 17690.    
				 23550.    
				 30246.    
				 37778.    
				 46145.    
				 55348.    
				 65385.    
				 76256.    
				 87962.    
				0.10050E+06
				0.11388E+06
				0.12809E+06
				0.14313E+06
				0.15901E+06
				0.17573E+06
				0.19329E+06
				0.21169E+06
				0.23093E+06
				0.25102E+06
				0.27196E+06
				0.29377E+06
				0.31643E+06
				0.33998E+06
				0.36440E+06
				0.38972E+06
				0.41595E+06
				0.44311E+06
				0.47121E+06
				0.50028E+06
				0.53034E+06
				0.56143E+06
				0.59357E+06
				0.62680E+06
				0.66117E+06
				0.69672E+06
				0.73349E+06
				0.77154E+06
				0.81092E+06
				0.85168E+06
				0.89387E+06
				0.93753E+06
				0.98269E+06
				0.10293E+07
				0.10775E+07
				0.11270E+07
				0.11778E+07
				0.12297E+07
				0.12823E+07
				0.13353E+07
				0.13880E+07
				0.14397E+07
				0.14894E+07
				0.15363E+07
				0.15789E+07
				0.16162E+07
				0.16467E+07
				0.16695E+07
				0.16835E+07];


%	calculate the percent frequency differences between the frequencies of the 
%   10 element Block Lanczos model and the theoretical frequencies

	diff_16_10 = abs(freq10bl(1:10) - freq16(1:10));

	pc_diff_16_10 = 100*diff_16_10./freq10bl(1:10);

	semilogy(pc_diff_16_10,'k-')
	title('percent frequency differences between Lanczos 10 element and theoretical models')
	xlabel('mode number')
	ylabel('percent resonant frequency difference')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

	loglog(freq16(1:10),pc_diff_16_10,'k-')
	title('percent frequency differences between Lanczos 10 element and theoretical models')
	xlabel('resonant frequency, theoretical model')
	ylabel('percent resonant frequency difference, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

%	calculate the percent frequency differences between the frequencies of the 
%   10 element Reduced model and the theoretical frequencies

	diff_16_10 = abs(freq10red(1:10) - freq16(1:10));

	pc_diff_16_10 = 100*diff_16_10./freq10red(1:10);

	semilogy(pc_diff_16_10,'k-')
	title('percent frequency differences between Reduced 10 element and theoretical models')
	xlabel('mode number')
	ylabel('percent resonant frequency difference')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

	loglog(freq16(1:10),pc_diff_16_10,'k-')
	title('percent frequency differences between Reduced 10 element and theoretical models')
	xlabel('resonant frequency, theoretical model')
	ylabel('percent resonant frequency difference, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

	plot(freq2bl,'k-')
	hold on
	plot(freq4bl,'k-')
	plot(freq6bl,'k-')
	plot(freq8bl,'k-')
	plot(freq10bl,'k-')
	plot(freq12bl,'k-')
	plot(freq16bl,'k-')
	plot(freq32bl,'k-')
	plot(freq64bl,'k-')
	title('Block Lanczos resonances versus number of modes extracted')
	xlabel('number of modes, 2 times number of elements')
	ylabel('resonant frequency, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

	hold off

	loglog(freq2bl,'k-')
	hold on
	loglog(freq4bl,'k-')
	loglog(freq6bl,'k-')
	loglog(freq8bl,'k-')
	loglog(freq10bl,'k-')
	loglog(freq12bl,'k-')
	loglog(freq16bl,'k-')
	loglog(freq32bl,'k-')
	loglog(freq64bl,'k-')
	title('Block Lanczos resonances versus number of modes extracted')
	xlabel('number of modes, 2 times number of elements')
	ylabel('resonant frequency, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

	hold off

	axis([1 inf 100 100000])
	disp('execution paused to display figure, "enter" to continue'); pause


	plot(freq2red,'k-')
	hold on
	plot(freq4red,'k-')
	plot(freq6red,'k-')
	plot(freq8red,'k-')
	plot(freq10red,'k-')
	plot(freq12red,'k-')
	plot(freq16red,'k-')
	plot(freq32red,'k-')
	plot(freq64red,'k-')
	title('Reduced resonances versus number of modes extracted')
	xlabel('number of modes, number of elements')
	ylabel('resonant frequency, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

	hold off

	loglog(freq2red,'k-')
	hold on
	loglog(freq4red,'k-')
	loglog(freq6red,'k-')
	loglog(freq8red,'k-')
	loglog(freq10red,'k-')
	loglog(freq12red,'k-')
	loglog(freq16red,'k-')
	loglog(freq32red,'k-')
	loglog(freq64red,'k-')
	title('Reduced resonances versus number of modes extracted')
	xlabel('number of modes, number of elements')
	ylabel('resonant frequency, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

	hold off

	axis([1 inf 100 100000])
	disp('execution paused to display figure, "enter" to continue'); pause


%	calculate the frequency differences between the two 10 element models, 
%   reduced and Block Lanczos

	diff_10_10 = freq10bl(1:10) - freq10red;

	plot(diff_10_10,'k-')
	title('frequency differences between the reduced and Lanczos 10 element models')
	xlabel('mode number')
	ylabel('resonant frequency difference, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

%	calculate the percent frequency differences between the two 10 element models, 
%   reduced and Block Lanczos

	pc_diff_10_10 = 100*(freq10bl(1:10) - freq10red)./freq10bl(1:10);

	plot(diff_10_10,'k-')
	title('percent frequency differences between the reduced and Lanczos 10 element models')
	xlabel('mode number')
	ylabel('percent resonant frequency difference, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

	plot(freq10bl(1:10),pc_diff_10_10,'k-')
	title('percent frequency differences between the reduced and Lanczos 10 element models')
	xlabel('resonant frequency, 10 element Lanczos model')
	ylabel('percent resonant frequency difference, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

%	calculate the percent frequency differences between the first 16 frequencies of the 
%   64 element Block Lanczos model and the theoretical frequencies

	diff_16_64 = freq16 - freq64bl(1:16);

	pc_diff_16_64 = 100*diff_16_64./freq64bl(1:16);

	semilogy(pc_diff_16_64,'k-')
	title('percent frequency differences between first 16 modes of Lanczos 64 element and theoretical models')
	xlabel('mode number')
	ylabel('percent resonant frequency difference, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

	loglog(freq16,pc_diff_16_64,'k-')
	title('percent frequency differences between first 16 modes of Lanczos 64 element and theoretical models')
	xlabel('resonant frequency, theoretical model')
	ylabel('percent resonant frequency difference, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

%	calculate the percent frequency differences between the two 64 element models, 
%   reduced and Block Lanczos

	diff_64_64 = freq64bl(1:64) - freq64red;

	pc_diff_64_64 = 100*(freq64bl(1:64) - freq64red)./freq64bl(1:64);

	plot(pc_diff_64_64,'k-')
	title('percent frequency differences between the reduced and Lanczos 64 element models')
	xlabel('mode number')
	ylabel('percent resonant frequency difference, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

	plot(freq64bl(1:64),pc_diff_64_64,'k-')
	title('percent frequency differences between the reduced and Lanczos 64 element models')
	xlabel('resonant frequency, 64 element Lanczos model')
	ylabel('percent resonant frequency difference, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

%	calculate the frequency differences between the 10-element model and the
%   64 element model for the reduced models (only use the first 10 eigenvalues
%	from the 64 element model)

	diff_10_64 = freq64red(1:10) - freq10red;

	plot(diff_10_64,'k-')
	title('frequency differences between the reduced 10 and 64 element models')
	xlabel('mode number')
	ylabel('resonant frequency difference, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

%	calculate the % frequency differences between the 10-element model and the
%   64 element model for the reduced models (only use the first 10 eigenvalues
%	from the 64 element model)

	pc_diff_10_64 = 100*(freq64red(1:10) - freq10red)./freq64red(1:10);

	plot(pc_diff_10_64,'k-')
	title('percent frequency differences between the reduced 10 and 64 element models')
	xlabel('mode number')
	ylabel('percent resonant frequency difference, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause

	plot(freq64red(1:10),pc_diff_10_64,'k-')
	title('resonant frequency differences between the reduced 10 and 64 element models')
	xlabel('resonant frequency, 64 element model')
	ylabel('percent resonant frequency difference, hz')
	grid off
	disp('execution paused to display figure, "enter" to continue'); pause