
- •Введение
- •1. Запуск математического пакета maple, структура окна
- •2. Описание основных команд–процедур пакета maple
- •2.1. Арифметические операции, типы данных в Maple
- •2.2. Структура команд Maple
- •Команда (параметр_1, параметр_2, …, параметр_n);
- •2.3. Преобразование математических выражений
- •2.3.1. Выделение частей арифметических выражений
- •2.3.2. Тождественные преобразования арифметических выражений
- •Collect(выражение, имя переменной относительно которой приводятся подобные);
- •3. Функции maple, решение уравнений и неравенств
- •3.1. Способы задания функций и замена переменных
- •Piecewise(cond_1,f1,cond_2,f2, …).
- •3.2. Операции оценивания
- •3.3. Решение уравнений и их систем
- •3.4. Решение неравенств
- •4. Графические и геометрические возможности пакета maple
- •4.1. Графические возможности пакета Maple
- •4.1.1. Двумерная графика
- •Polarplot ([переменная, выражение,параметры переменной],опции);
- •Pointplot(data_list,опции);
- •Sphereplot ((выражение), параметры_, параметры_, опции);
- •Sphereplot ([r_выражение, _выражение, _выражение], параметр_1, параметр_2, опции);
- •4.1.2. Трехмерная графика
- •4.2. Анимация
- •4.3. Геометрические пакеты Maple
- •4.3.1. Стереометрия
- •4.3.2. Планиметрия
- •5. Дифференциальное и интегральное исчисление функции одной переменной
- •5.1. Вычисление пределов в Maple для некоторых математических операций существует по две команды: прямого и отложенного исполнения. Имена команд состоят из одинаковых букв за исключением первой:
- •5.2. Дифференцирование
- •5.2.1. Вычисление производных
- •5.2.2. Дифференциальный оператор
- •5.3. Исследование функций
- •5.3.1. Непрерывность функции и точки разрыва
- •5.3.2. Экстремумы. Наибольшее и наименьшее значение функции
- •5.3.3. Исследование функции по общей схеме
- •{Arctan(0)}
- •6. Интегрирование
- •6.1. Аналитическое и численное интегрирование
- •6.2. Определенные интегралы
- •Is assumed to be: RealRange(Open(-1),3)
- •6.3. Основные методы интегрирования
- •7. Дифференциальные уравнения
- •7.1. Аналитическое решение дифференциальных уравнений
- •7.1.1. Общее решение дифференциальных уравнений
- •7.1.2. Фундаментальная (базисная) система решений
- •7.1.3. Решение задачи Коши или краевой задачи
- •7.1.4. Системы дифференциальных уравнений
- •7.1.5. Приближенное решение дифференциальных уравнений с помощью степенных рядов
- •7.2. Численное решение дифференциальных уравнений
- •7.2.1. Формат команд dsolve и odeplot для нахождения численного решения дифференциальных уравнений в графическом виде
- •Italic, 12]):
- •Italic, 12]):
- •7.2.2. Пакет графического представления решений дифференциальных уравнений dEtools
- •7.2.3. Построение фазовых портретов систем дифференциальных уравнений
- •8. Функции многих переменных, векторный анализ, ряды, интегральные преобразования
- •8.1. Дифференциальное исчисление функций многих переменных
- •8.1.1. Частные производные
- •8.1.2. Локальные и условные экстремумы функций многих переменных
- •Maximize(f,{x1,…,xn},range);
- •Minimize(f,{x1,…,xn}, range);
- •8.2. Интегральное исчисление функций многих переменных
- •8.3. Векторный анализ
- •8.4. Ряды и произведения
- •8.4.1. Вычисление суммы ряда и произведений
- •8.4.2. Разложение функции в степенной ряд и ряд Тейлора
- •8.5. Интегральные преобразования
- •8.5.1. Преобразование Фурье
- •8.5.2. Преобразование Лапласа
- •9. Линейная алгебра
- •9.1. Векторная алгебра
- •Warning, the protected names norm and trace have been redefined and unprotected
- •9.2. Действия с матрицами
- •9.2.1. Определение матрицы
- •9.2.2. Арифметические операции с матрицами
- •9.2.3. Определители, миноры и алгебраические дополнения, ранг и след матрицы
- •9.2.4. Обратная и транспонированная матрицы
- •Evalm(1/a);
- •Inverse(a).
- •9.2.5. Определение типа матрицы
- •1*Sqrt(3)/2,-1/2]]);
- •9.2.6. Функции от матриц
- •9.3. Спектральный анализ матрицы
- •9.4. Системы линейных уравнений, матричные уравнения
- •For имя переменной цикла in выражение 1 do выражение 2 od;
- •10.3. Процедуры-функции
- •10.4. Процедуры
- •Writeto (“имя файла”)
- •Appendto (“имя файла”)
- •Print(список Maple-выражений, перечисляемых через запятую);
- •Lprint(список Maple-выражений, перечисляемых через запятую);
- •Readdata(“имя файла”, тип переменной(integer/float – последний тип устанавливается по умолчанию),счетчик чисел);
- •Readlib(readdata):
- •12. Использование математического пакета maple для научных исследований
- •12.1. Исследование влияния изменяемых параметров плоской помольной камеры мельницы противоточного действия на скорость энергоносителя
- •12.1.1. Постановка задачи
- •12.1.2. Решение задачи
- •12.2. Определение коэффициента полезного заполнения материалом деформируемых бочкообразных камер устройств для тонкого и свертонкого помола материалов
- •12.2.1. Расчетная схема
- •12.2.2. Решение задачи
- •12.3. Исследование влияния основных конструктивных и технологических параметров на режимы работы трубной мельницы
- •Formula(1.8)
- •Formula(1.9)
- •Formula(2.1)
- •Formula(2.2)
- •Formula(2.3)
- •Formula(3.1)
- •Formula(3.2)
- •Formula(3.3)
- •13. Формат наиболее используемых команд аналитического пакета maple
- •Заключение
- •Библиографический список
- •Оглавление
12.1. Исследование влияния изменяемых параметров плоской помольной камеры мельницы противоточного действия на скорость энергоносителя
12.1.1. Постановка задачи
Струйные мельницы являются разновидностью ударных измельчителей и состоят из разгонного аппарата (одного или нескольких), в котором струя газа-энергоносителя сообщает, скорость частицам обрабатываемого материала, и камеры, в которой происходит взаимодействие потоков материала между собой и(или) со специальными отбойными поверхностями. В качестве энергоносителя в струйных мельницах чаще всего применяется воздух, реже – инертный газ, водяной пар, продукты сгорания.
Струйный помол дает возможность сочетания помола и разделения со смешением, сушкой и другими технологическими процессами. А работа в замкнутом цикле обеспечивает минимальное выделение пыли в окружающую среду.
Любой струйный аппарат включает в себя эжектор, представляющий собой узел, в котором происходит смешение и обмен энергией двух потоков (основного и эжектируемого) и помольную камеру, в которой взаимодействуют смешанные потоки. Ускоренные энергоносителем в разгонных трубках эжекторов частицы попадают в помольную камеру, а затем в зону встречи струй (рис. 12.1.).
Струя, выходящая из разгонной трубки, не сразу заполняет все поперечное сечение помолной камеры, струя в месте входа в нее отрывается от стенок и дальше движется в виде свободной струи, отделенной от остально среды поверхностью раздела. Поверхность раздела неустойчива, на ней возникают вихри, в результате чего струя перемешивается с окружающей средой.
При истечении струи из разгонной трубки скорости потока в ее выходном сечении 1-1 во всех точках сечения равны между собой. На протяжении длины– начальном участке, осевая скорость постоянна по величине и равна скорости на срезе разгонной трубки V0. В области треугольника АВС (рис. 12.1.) во всех точках струи скорости энергоносителя равны между собой и также равны V0 - эта область образует так называемое ядро струи. Далее осевая скорость постепенно уменьшается и на основном участке длинной lосн осевая скорость VОС < V0.
Рис. 12.1. Схема струи в помольной камере
Известно, что скорость энергоносителя от среза разгонной трубки до плоскости соударения струй изменяется по закону
,
(12.1)
где Vz – скорость энергоносителя с помольной камере на расстоянии z от среза разгонной трубки, м/с;
V0 – скорость энергоносителя на срезе разгонной трубки, м/с;
z0 – расстояние от среза разгонной трубки до плоскости встречи струй, м.
При определении изменения кинетической энергии конечного объема сплошной среды, необходимо знать работу сил межкомпонентного взаимодействия частиц измельченного материала и энергоносителя. Эта работа зависит от вектора силы динамического воздействия энергоносителя на частицу, которая вычисляется следующим образом
,
(12.2)
где R – вектор силы динамического воздействия воздуха на частицу, Н;
Fm – площадь сечения частицы, м2;
,
(12,3)
где dэ – эквивалентный диаметр частицы, м;
- коэффициент
лобового сопротивления частиц;
V – вектор скорости энергоносителя, м/с;
U – вектор скорости частицы материала, м/с;
-
плотность энергоносителя, кг/м3.
Изменение скорости частицы после вылета из разгонной трубки и движение вдоль камеры помола до оси Y – границы помольной камеры определяется уравнением динамики одиночной частицы
,
(12.4)
или в проекции на ось Z
,
.
(12.5)
Введем обозначение
.
(12.6)
Получим
.
(12.7)
Обозначим
,
(12.8)
где m – масса частицы измельчаемого материала, кг.
,
(12.9)
где
- плотность частиц измельчаемого
материала, кг/м.
Выражение (12.7) примет вид
.
(12.10)
Полученное уравнение может быть использовано для определения изменения скорости частиц, измельчаемого материала в помольной камере на участве от среза разгонных трубок до области взаимодействия встречных потоков.
Система дифференциальных уравнений, описывающих процесс изменения скорости частиц и энергоносителя в помольной камере от среза разгонной трубки до области соударения встречных потоков
.
(12.11)
Расстояние lстр – между срезом разгонной трубки и серединной плоскость в помольной камере выбрано из условия
,
(12.12)
где dтр = 18 диаметр разгонной трубки, мм.