3. Амплитудно-частотная и фазовая характеристики

Разделив каждый член уравнения (1) на массу m, получим следующие дифференциальное уравнение

(6)

где - частота собственных колебаний;

Из уравнения (6), представленного в операторном виде, подставляя , находим АЧХ и ФЧХ колебательной системы прибора [1]:

а) чувствительность механической системы к ускорению a

(7)

где - отношение частоты приложенной (вынуждающей) силы к частоте собственных колебаний;

б) чувствительность механической системы к скорости v

(8)

в) чувствительность механической системы к перемещению массы y

(9)

M - файл программы расчета АЧХ и ФЧХ механической колебательной системы инерционного прибора.

% Анализ АЧХ и ФЧХ механической системы

% Анализ чувствительности системы к ускорению

% Программу составил доцент Литвин А.В.

% кафедра "Приборостроение" ДГТУ май 2004

%

% Введем исходные параметры

% колебательной механической системы

% m - масса тела в [г];

% d - диаметр проволоки пружины;

% ni- число свободных витков пружины;

m=1e-3*input('Input mass [g] m=');

d=1e-3*input('Input spring wire diameter [mm] d=');

ds=1e-3*input('Input average spring diameter [mm] ds=');

ni=input('Input spring coils number ni=');

G=83670e6; % Модуль упругость второго рода [Па]

nw=ni+0.5; % Число рабочих витков пружины

k=G*d^4/(8*nw*ds^3); % Жесткость пружины

omega0=sqrt(k/m);

% Диапазон отношений частоты действующей силы

% к частоте собственных колебаний

eta=0:0.01:4;

etf1=0:0.01:0.99;

etf2=1.01:0.01:4;

eta=eta';

etf1=etf1';

etf2=etf2';

etf12=etf1.^2;

etf22=etf2.^2;

eta2=eta.^2;

% Диапазон коэффициентов демпфирования

dzeta=[0.125 0.25 0.6];

dzeta2=dzeta.^2;

% Вычисление матрицы чувствительности (АЧХ)

% к ускорению

% и матрицу ФЧХ

Sna=[]; Fna=[];

for i=1:3

Sna(:,i)=1./((1-eta2).^2+eta2.*dzeta2(i)*4).^0.5;

end

for i=1:3

Fna1(:,i)=-atan((etf1.*dzeta(i).*2)./(1-etf12));

Fna2(:,i)=-atan((etf2.*dzeta(i).*2)./(1-etf22));

end

Fna2=Fna2-pi;

Fna=[Fna1;Fna2];

etf=[etf1;etf2];

Fnaan=(Fna.*180)./pi;

% Графики зависимости чувствительности Sa(eta)

% к ускорению от eta=omega/omega0

figure(1)

plot(eta,Sna), grid

title('Amplitude-ratio vs. Frequency ratio eta')

xlabel('Frequency ratio eta=omega/omega0')

ylabel('Relative acceleration sensitivity')

pause

figure(2)

plot(etf, Fna), grid

title('Phase angle vs. Frequency ratio eta')

xlabel('Frequency ratio eta=omega/omega0')

ylabel('Phase angle, radian')

pause

figure(3)

plot(etf, Fnaan), grid

title('Phase angle vs. Frequency ratio eta')

xlabel('Frequency ratio eta=omega/omega0')

ylabel('Phase angle, degree')

pause

4. Задание

1. Ответьте на приведенные ниже контрольные вопросы.

2. Наберите тексты программ в системе MATLAB, выполните их отладку по тестовому примеру: m=50 г.; d=0.5 мм; ds=5 мм; ni=6.

3. Выполните анализ механической системы прибора по заданию преподавателя.

4. Составьте собственные программы для анализа АЧХ и ФЧХ для чувствительности механической системы прибора к скорости и перемещению массы.

5. Выполните анализ АЧХ и ФЧХ механической системы прибора по заданию преподавателя.

6. Проанализируйте полученные результаты и дайте заключение.

7. Составьте эквивалентную схему механической системы.

8. Составьте граф механической системы по эквивалентной схеме.

9. Приведите электрическую схему – аналог механической системы.

5. Контрольные вопросы

1. Что является математической моделью динамической системы?

2. Какими уравнениями описываются колебательные системы?

3. Какие свойства системы определяет ее постоянная времени?

4. Почему при определении жесткости пружины используется модуль упругости второго рода?

5. Чему соответствует сила Ньютона в электрических системах?

6. Что получается в результате решения дифференциального уравнения?

7. Каким методом решается система ОДУ в лабораторной работе?

8. При каких условиях будет наблюдаться затухание колебаний?

9. Как изменяется график движения массы с увеличением коэффициента демпфирования?

10. Какую систему можно назвать осциллятором?

11. Какую траекторию имеют затухающие колебания линейного осциллятора в фазовой плоскости?

12. Какие оси имеет фазовая плоскость?

13. Какое явление наблюдается в колебательной системе при равенстве частот собственных колебаний системы и вынуждающей силы?

14. В каких координатных осях построены графики АЧХ и ФЧХ?

6. Содержание отчета

1. Схема прибора с указанием параметров механической системы.

2. Эквивалентная схема и ее граф.

3. Математическая модель механической системы.

4. Алгоритм программы анализа динамики механической системы.

5. Результаты анализа: Графики переходного процесса; Графики АЧХ и ФЧХ; Заключение по результатам анализа.

6. Текст программы по пункту 4 задания.

7. Результаты анализа, выполненного по программе пункта 6.