- •Министерство образования и науки Российской Федерации
- •Содержание
- •Введение
- •1. Современное состояние проблемы моделирования систем
- •1.1. Моделирование как метод научного познания. Философские аспекты моделирования
- •1.2. Использование моделирования при исследовании и проектировании систем
- •1.2.1. Особенности разработки систем
- •1.2.2. Особенности использования моделей
- •1.2.3. Перспективы развития методов и средств моделирования систем
- •2. Основные понятия теории моделирования систем
- •2.1. Принцип системного подхода в моделировании систем
- •2.1.1. Структура системы – совокупность связей между элементами системы
- •2.1.2. Экспериментальные исследования систем
- •2.2. Стадии разработки моделей
- •2.3. Понятие подобия
- •2.3.1. Общие положения
- •2.3.2. Основные понятия теории размерности
- •2.3.3. Примеры подобия
- •2.4. Общая характеристика проблемы моделирования систем
- •2.4.1. Объект моделирования.
- •2.4.2. Характеристики моделей систем
- •2.4.3. Цели моделирования систем
- •2.5. Классификация видов и методов моделирования систем
- •2.5.1. Классификационные признаки
- •2.5.2. Математическое моделирование.
- •2.6. Построение модели
- •2.7. Разработка вычислительного метода
- •2.8. Проверка (тестирование) модели
- •3. Математическое моделирование
- •3.1. Задачи и цели исследования математических моделей
- •3.2. Методология математического моделирования. Системный анализ
- •3.2.1. Понятие системы
- •3.2.2. Этапы системного анализа и декомпозиция
- •3.2.3. Экспертные оценки
- •3.3. Классификация математических моделей
- •3.4. Методы формализованного описания системы
- •3.4.1. Математическая модель по “входу-выходу”
- •3.4.2. Математическая модель в пространстве состояний
- •3.4.3. Описание линейных систем в пространстве состояний
- •3.4.4. Реализация систем в пространстве состояний
- •3.5. Методы построения математических моделей и их применение в сапр
- •3.5.1. Методы построения математических моделей
- •3.5.2. Математические модели с точки зрения сапр
- •3.5.4. Методика составления уравнений динамики элементов сау
- •3.6. Математические модели системы управления. Понятие об оптимальном управлении
- •4. Экспериментальное определение динамических характеристик объектов моделирования
- •4.1. Понятие о динамических характеристиках объектов
- •4.2. Определение динамических характеристик элементов систем по временным характеристикам
- •4.2.1. Определение статических характеристик
- •4.2.2. Определение динамических характеристик объектов с помощью периодических воздействий
- •4.4.1. Временные характеристики и их свойства
- •4.4.2. Определение характеристик апериодического звена
- •4.4.3. Определение характеристик колебательного звена
- •4.3. Формы описания динамических свойств объектов
- •4.4. Синтез пассивных двухполюсников и четырехполюсников
- •4.3.1. Разложение передаточной функции активного четырехполюсника
- •4.3.2. Способы синтеза двухполюсников
- •4.5. Экспериментальная отработка характеристик системы управления движущимся объектом
- •4.5.1. Общие положения
- •4.5.2. Алгоритмы обработки внешнетраекторных измерений
- •5. Динамические свойства воспринимающих элементов и датчиков
- •5.1. Основные определения и понятия
- •5.1.1. Понятие датчика
- •5.1.2. Классификация датчиков
- •5.2. Основные характеристики датчиков
- •5.2.1. Погрешности измерений
- •5.2.2. Чувствительность датчиков
- •5.2.3. Быстродействие датчика
- •5.3. Схемы формирования сигналов пассивных датчиков
- •5.3.1. Общие характеристики
- •5.4. Оптические датчики
- •5.4.1. Определения и основные зависимости
- •5.4.2. Фоторезисторы
- •5.4.3. Фотодиоды
- •5.4.4. Тепловые приемники излучения
- •5.4.5. Датчики изображения
- •5.4.6. Волоконная оптика
- •5.5. Датчики температуры
- •5.5.1. Методы измерения температуры
- •5.6. Датчики положения и перемещения
- •5.6.1. Методы определения положения и перемещения
- •5.6.2. Резисторные потенциометры
- •5.6.3. Индуктивные датчики
- •5.6.4. Емкостные датчики
- •5.6.5. Цифровые датчики
- •5.6.6. Датчики близости
- •5.7. Датчики деформации
- •5.7.1. Основные определения
- •5.7.2. Основные положения
- •5.8. Тахометрические датчики
- •5.8.1. Электродинамическая тахометрия
- •5.8.2. Импульсная тахометрия
- •5.8.3. Гирометры
- •5.9. Датчики ускорения, вибрации и удара
- •5.9.1. Общие положения
- •5.9.2. Принцип действия сейсмических датчиков
- •5.10. Датчики скорости, расхода и уровня жидкости
- •5.10.1. Элементарные понятия
- •5.10.2 Датчики и методы измерения скорости жидкости
- •5.10.3. Измерение расхода жидкости
- •5.10.4. Измерение и указание уровня жидкости
- •5.11. Датчики влажности
- •5.11.1. Определения
- •5.11.2. Гигрометры
- •5.12. Акустические датчики
- •5.12.1. Распространение плоской волны
- •5.12.2. Распространение трехмерной волны
- •5.12.3. Микрофоны
- •5.12.4. Измерение интенсивности
- •6. Основы технологии имитационного моделирования
- •6.1. Основные определения и понятия
- •6.2. Область применения и классификация имитационных моделей
- •6.3. Описание поведения системы
- •6.3.1. Общие положения.
- •6.3.2. Методика моделирования случайных факторов
- •6.3.3. Два подхода к моделированию случайных чисел
- •6.4. Оценка качества псевдослучайных чисел
- •6.5. Оценка качества имитационного моделирования
- •7. Методы испытаний систем управления и их применение в системах автоматизированного проектирования (сапр)
- •7.1. Полунатурное моделирование
- •7.1.1. Общие положения
- •7.1.2. Автоматизация испытаний на основе полунатурного моделирования
- •8. Анализ систем управления с эвм
- •8.1. Основные задачи
- •8.2. Особенности систем управления с эвм
- •8.2. Основные положения из теории дискретных линейных систем
- •8.2.1. Последовательности
- •8.2.2. Линейные системы с постоянными параметрами
- •8.2.3. Разностные уравнения
- •8.2.3.1. Решение разностных уравнений методом прямой подстановки
- •8.3. Расчет цифровых фильтров по фильтрам непрерывного времени
- •8.3.1 Методика синтеза цифровых фильтров. Общие положения
- •8.3.2 Методы дискретизации аналоговых фильтров
- •8.3.3. Геометрическая интерпретация методов расчета цифровых фильтров по фильтрам непрерывного времени
- •9. Моделирование свойств объектов с помощью системыMatLab
- •9.1. Введение
- •9.2. MatLab как научный калькулятор
- •9.2.1. Командное окно
- •9.2.2. Операции с числами
- •9.2.3. Простейшие операции с векторами и матрицами
- •9.2.4. Некоторые функции прикладной численной математики
- •9.2.5. Построение простейших графиков
- •9.3. Исследование линейных стационарных систем (лсс)
- •9.3.1. Классы пакета control.L
- •9.3.2. Ввод и преобразование моделей
- •Пример создания модели
- •9.3.3. Анализ системы
- •9.4. Моделирование динамических процессов с помощью подсистемы MatLab simulink
- •9.4.1. Краткие сведения о подсистеме MatLab simulink
- •9.4.2. Запуск подсистемы simulink
- •9.4.3. Создание модели
- •9.4.4. Некоторые основные приемы подготовки и редактирования модели
- •9.4.5. Установка параметров моделирования и его выполнение
- •9.2.2. Результат составления модели
- •Приложения п1. Динамические характеристики объектов моделирования
- •П2. Примеры составление функциональной и структурной схемы динамической системы
- •П2.1. Система управления угловой скорости вращения ротора двигателя при условии действия постоянного возмущения
- •П2.2. Система сопровождения цели
- •П2.3. Система автоматического наведения летательного аппарата на объект
- •П2.4. Система управления уровнем жидкости
- •П2.5. Система управления экономическими параметрами
- •Использованные источники
- •Основы теории и практики моделирования динамических систем
5.12.2. Распространение трехмерной волны
Волны в этом случае не являются плоскими и происходит уменьшение амплитуды с увеличением расстояния от источника до точки наблюдения (что следует из закона сохранения энергии).
Рассмотрим случай сферического источника с пульсирующим потоком энергии .
В этом случае имеют место следующие соотношения:
- давление
; (5.149)
- скорость
, (5.150)
- интенсивность
, (5.151)
Из этих выражений видно, что радиальная составляющая интенсивности убывает обратно пропорционально квадрату расстояния от источника при любых .
5.12.3. Микрофоны
Из названных в предыдущих разделах механических величин, характеризующих звуковые волны можно измерить только давление, а скорости частиц, соответствующие звуковым волнам, слишком малы и, следовательно, их нельзя измерить обычными анемометрами, используемых при исследовании жидкостей.
Однако компоненты скорости звука можно определить косвенным образом с использованием уравнения динамики, связывающего градиент давления и производную скорости по времени (см. (5.143))
.
Эту операцию можно произвести, используя датчик, применяемый для измерения разности давлений. Этот датчик, основанный на измерении градиента давления, называется микрофоном.
Другой часто используемый метод состоит в обработке с помощью электроники разности сигналов с двух микрофонов.
5.12.3.1. Принцип действия микрофона
Микрофон представляет собой датчик, преобразующий акустический сигнал в электрический.
Классическим можно считать микрофон, применяемый в телефонной трубке.
Основан на необратимом явлении, под действием давления происходит изменение сопротивления контакта частиц угля, заполняющих объем, закрытый мембраной.
Но микрофон этого типа не приспособлен для физических измерений, т.к. имеет ограниченную полосу пропускания.
Для создания измерительных микрофонов предпочтительнее методы преобразования, основанные на обратимых явлениях электромагнитного, электростатического, пьезоэлектрического типов.
5.12.3.2. Основные типы микрофонов
Выделяются группы микрофонов, предназначенных для измерения:
- давления;
- градиента давления.
Первые осуществляют физические измерения давления, вторые - физические измерения скорости звука.
Важное преимущество последних – их чувствительность к направлению распространения.
Для получения более выраженной направленности можно использовать микрофоны, чувствительные ко второму градиенту.
5.12.3.3. Электро-механо-акустические аналоги
Аналогия “течение – ток”
Сопоставляя уравнения для электрических цепей и уравнения акустических явлений, можно обнаружить некоторые аналогии.
Рассмотрим акустическую схему – простой акустический контур, представленный на рисунке 5.44.
На рисунке:
капиллярный канал;
полость объема V;
трубка с площадью сечения A;
u – скорость, p – давление.
Поскольку равенство давлений по разные стороны разграничивающих плоскостей означает равенство плотностей, из закона сохранения массы вытекает, что расход, вычисляемый как
, (5.152)
при переходе через эти плоскости сохраняется (т.е. течение сохраняется). Это течение можно рассматривать как аналог электрического тока, т.е. можно утверждать
.
Аналогия акустической массы
Находящийся в трубе газ массой перемещается как целое под действием разности давлений. Тогда справедливы соотношения
,
откуда
(5.153)
где - расход газа.
Приведем уравнение основного закона механики (динамики) имеет вид
, (5.154)
где - сила, действующее на твердое тело с массой, под действием которой оно движется со скоростью.
Также известно выражение для закона Ома
, (5.155)
в котором - индуктивность,- разность потенциалов,- ток.
Как следует из рассмотрения полученных соотношений, уравнение для газа (5.153) близко по форме к уравнению механики (5.154) и выражению закона Ома. Формально сомножитель в правой части выражения (5.153) можно интерпретировать какакустическую массу.
Таким образом, существует аналогия между механической массой , индуктивностьюи акустической массой. При этом перемещение массы играет роль тока
, (5.156)
а разность давлений - разность потенциалов
. (5.157)
Акустическое сопротивление
Падение давления в капилляре, вызванное вязкостью, пропорционально скорости (при ламинарном течении).
Таким образом, получаем
, (5.158)
где - аналог гидравлического сопротивления, которое называетсяакустическим сопротивлением.
Акустическая упругость
Из-за малых размеров сжатие газа в полости происходит без смещения его центра массы. Процесс сжатия описывается законами термодинамики.
Предполагая сжатие адиабатическим (изменение количества теплоты, сообщенное в процессе, отсутствует, при котором можно принять , где- удельная теплоемкость при постоянном давлении и объеме) получаем
. (5.159)
Справка.
Выражение (5.159) можно получить из следующих выкладок.
Введем обозначение
Из данного требования вытекает, что приращение переменной .
Найдем полный дифференциал .
.
Отсюда имеем .
Пусть
.
Тогда
. (5.160)
Сомножитель служит аналогом емкости в электричестве. Назовем егоакустической емкостью.
Уравнение акустического канала
Для представленной на рисунке схемы получим дифференциальное уравнение, связывающее давление с потоком.
Имеем
.
В правой части слагаемые заменим соответствующими выражениями из (5.153), (5.158), (5.160) и получим
. (5.161)
Полученное уравнение подобно закону Ома для последовательно включенных активного сопротивления , индуктивностии емкости
. (5.162)
и механическому уравнению, описывающему затухающее движение тела с массой M с одной степенью свободы
. (5.163)
Сопоставление последних трех уравнений указывает на наличие аналогии с активными и реактивными сопротивлениями, соответствующие обозначения для которых приведено в таблице 5.5.
Таблица 5.5. Электро-механо-акустические аналоги
Электрическая величина |
Механическая величина |
Акустическая величина |
Графическое обозначение |
Э.д.с., e |
Механическая сила F |
Акустическое давление, p |
- |
Ток, i |
Линейная скорость, V |
Объемный расход, Au |
|
Индуктивность, L |
Масса, M |
Акустическая масса, |
|
Электрическое сопротивление, |
Механическое сопротивление, |
Акустическое сопротивление, |
|
Емкость, c |
Податливость (обратная величина жесткости), 1/k |
Акустическая упругость, емкость,
|
|
Эквивалентная схема микрофона
Все электроакустические системы можно разделить более или менее условно на три части:
- акустическую;
- механическую;
- электрическую.
Руководствуясь изложенными выше принципами, можно составить электрическую схему, эквивалентную исходной акустической.
При этом все элементы схемы, через которые проходит один и тот же поток, размещаются последовательно, а элементы, соответствующие частям общего потока, - параллельно.
Заметим, что в большинстве случаев сложение потоков проявляется в емкостях.
5.12.3.4. Микрофон, чувствительный к давлению
В микрофонах этого типа флуктуации давления приводят к смещению диафрагмы. А изменение атмосферного давления не вызывает какого-либо эффекта (деформации) из-за наличия капиллярного канала для выравнивания давления.
Подобные микрофоны содержат три элемента:
- полость;
- диафрагму;
- канал для подхода воздуха (капилляр),
которые можно представить соответственно с помощью емкости (ей соответствует импедансZс), импеданса реактивного типа ZD и активного сопротивления .
Получить выражение для определения давления, действующего на диафрагму.
Эквивалентная акустическая схема, построенная на основе аналогии с активным и реактивным сопротивлением, представлена на рисунке 5.45.
Предположим, что давление, действующее на диафрагму и капилляр одинаковы. Далее выполняем расчет по эквивалентной электрической схеме (рисунок 5.45,в)).
Используем законы Кирхгофа.
Согласно первому закону (следствие закона сохранения зряда, согласно которому в любом узле заряд, одного знака не может ни накапливаться, ни убывать: алгебраическая сума токов ветвей, сходящихся в узле, равна нулю) имеем
;
.
Согласно второму закону Кирхгофа (алгебраическая сумма напряжений всех участков замкнутого контура равна нулю; при составлении уравнений принято следующее правило: направление обхода участков контура – по часовой стрелке; со знаком плюс принять те слагаемые, в которых токи и э.д.с. имеют направление, совпадающее с направлением обхода) можно записать:
(5.164)
Решение системы (5.164) дает выражение
,
в котором (активное сопротивление мало по сравнению с импедансом). Откуда можно записать
. (5.165)
Учитывая, что , из (5.165) получим
(5.166)
В полученном выражении тока в соответствии с (5.156) является эквивалентом расхода газа (скорсти потока) , а напряжениев соответствии с (5.157) – эквивалентом давления . Следовательно, выражение (5.166) можно переписать для искомых величин
(5.167)
Определим силу, действующую на диафрагму.
,
где - площадь диафрагмы микрофона;
- давление на диафрагму.
В соответствии со схемой (рисунок 5.25, б), в)) падение напряжения на импедансе (давление на диафрагму) есть разница потенциалов на двух концах.
Это падение напряжения (давление) определяется по зависимости:
, (5.168)
в которой первое слагаемое - есть давление на левом конце, а второе слагаемое - на правом конце импеданса. Тогда сила будет равна
.
Последнее выражение с учетом (5.167) дает
. (5.168)
Значение в (5.168) можно подобрать таким образом, чтобы всегда имело место, а постоянная времени.
Этими параметрами определяется нижняя граница частот.
5.12.3.5. Микрофон, чувствительный к градиенту давления
Практически градиент можно определить, используя приближение: ,
где ,-давление в двух точках, расстояние между которыми.
Пусть плоская волна направлена вдоль оси OX0, образуя угол с осьюOX (рисунок 5.46). Уравнение волны по давлению, направленному по оси OX0, согласно (5.142) имеет вид,
(5.169)
При из (5.169) имеем
(5.170)
Сила, действующая на диафрагму площадью A, равна
(5.171)
Для плоской волны имеет место
. (5.172)
Подставив отсюда давление в (5.171) получим
. (5.172)
Из полученного выражения следует, что амплитуда силы, действующей на диафрагму, пропорциональна произведению линейной скорости на круговую частоту.
Этот вывод не зависит от формы волны. В самом деле, основное уравнение динамики имеет вид (что следует из закона сохранения количества движения, см. (5.143))
.
Отсюда следует
. (5.173)
Напротив, чувствительность к давлению зависит от формы волны.
Для плоской волны с углом падения имеем
и .
Используя выражения (5.172) и (5.169), получим
,
откуда вытекает
. (5.174)
Для сферической волны радиусом и с центом на осиOX0, удовлетворяющей условию
,
имеем
.
С учетом очевидных соотношений
получим
.
Отсюда имеем
. (5.175)
Из сравнения выражений (5.174) и (5.175) очевидно, что чувствительность к давлению действительнозависит от формы волны.
5.12.3.6. Микрофоны смешанного типа (чувствительные как к давлению, так и к градиенту давления)
Рассмотрим микрофон (рисунок 5.47), полость которогозакрыта, а на противоположной стороне имеется окно, регулируемое акустическим сопротивлением (на рисунке: 1 – полость; 2 – сопротивление; 3 – механическая связь с преобразователем; 4 – диафрагма).
Эквивалентная расчетная электрическая схема микрофона представлена ниже на рисунке 5.48.
Вследствие существованния зазора Δl между двумя поверхностями давления p1 и p2, приложенные к этим поверхностям, различаются.
Для случая плоской волны, образующей с осью микрофона угол θ, получим
;
;
.
Отсюда
. (5.176)
Отметим, что акустическое сопаротивление R значительно меньше, чем в случае, рассмотренном выше, поскольку относится к окну, а не к капиллярному отверстию.
Используя акустическую аналогию, можно записать (что следует из законов Кирхгофа)
, (5.177)
. (5.178)
Очевидно, что даление на диафрагму есть разность потенциалов на двух концах.
Отсюда имеем:
.
Выражение для давления получим их трех последних уравнениий (5.176) - (5.178)
,
откуда следует выражение для силы , действующей на диафрагму:
. (5.179)
При условии имеем
.
Умножая числитель и знаменатель последнего выражения на , получим
.
Вводя обозначение , из последнего уравнения получим
. (5.180)
Полученное выражеие для общего слусая.
Частный случай, когда можно принять , соответствует микрофону, чувствительному кдавлению (см. разд. 5.12.3.4, соотношение (5.168)).
При имеем микрофон, чувствительный кградиенту.
Пример.
Пусть ,.
Используя зависимость (5.180), получим линии постоянных давлений (диаграмму направленности) микрофона.
Для этого находим модуль давления , определяемого по соотношению (5.180), извыражения
На рисунке 5.49 приведены такие линии для двух значений и.