- •Введение
- •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.2.2. Решение задачи
Расчитаем коэффициент полезного заполнения материалом бочкообразных рабочих камер для трех типоразмеров:
1 – RC = L = 300 мм, D = 150 мм;
2 – RC = L = 500 мм, D = 250 мм;
3 – RC = L = 1000 мм, D = 500 мм.
Размер мелющих тел изменяется в пределах rШ = 4 - 18.
[> restart;
Используем пакет plot для вывода графиков на экран.
[> with(plots):
Warning, the name changecoords has been redefined
Вводим исходные данные для первого типоразмера рабочей камеры.
Колическтво секций у камеры
[> nn:=1;
Радиус кривизны бочкообразной рабочей камеры.
[> Rc:=300;
Длина бочкообразной камеры.
[> Lc:=300;
Диаметр рабочей камеры
[> D1:=150;
Мимальный радиус мелющих тел r=4 мм. Для организации цикла для расчета коэффициента полезного заполнения материалом обозначим радиус мелющих тел символом j, первый шаг цикла j:=j+1, поэтому
[> j:=3;
[>while j<18 do j:=j+1; m:=0; n:=0; s:=0; k1:=Lc/4/j-trunc(Lc/4/j); k:=trunc(Lc/4/j); yn:=(sqrt(Rc^2-((2*n-1)*j)^2)-Rc+D1/2); m:=0; n:=0; s:=0; while n<k do n:=n+1;l:=trunc(yn/j/2); s:=s+Pi*yn/j*l-Pi*l^2 od; n2:=2*nn*evalf(s); if (k1>=0.5) then k:=k+1 fi ;m:=0; n:=0; s:=0; yn:=(sqrt(Rc^2-((2*n-1)*j)^2)-Rc+D1/2); while n<k do n:=n+1;l:=trunc(yn/j/2); s:=s+Pi*yn/j*l-Pi*l^2 od;n1:=2*nn*evalf(s); n:=trunc((n2+n1)/2); y(x):=sqrt(Rc^2-x^2)-Rc+D1/2; Ic:=int(y(x)^2 ,x=0..Lc/2); re1:=Pi*Ic*2*nn; m1:=n*4*Pi*j^3/3;oo:=evalf(re1); ot[j]:=evalf(1-m1/re1) od;
ot4 := .4579582763
………………
ot18 := .4579582763
В связи с громоздкостью представления результатов вычислений предсталяем только окончательные и исходные данные.
Для второго типоразмера рабочей камеры RC = L = 500 мм,
D = 250 мм;
ot4 := .4667585
…………………….
ot18 := .5847463
Для третьего типоразмера рабочей камеры RC = L = 1000 мм, D = 500 мм;
ot4 := .4579582763
…………………….
ot18 := .4579582763
Построим график зависимости коэффициента полезного заполнения от размера мелющих тел в зависимости от геометрических парметров рабочей камеры по координатам точек (рис. 12.13).
[>date_list1:=[[8,0.45],[13,0.46],[18,0.45], [23,0.46],[28,0.47],[33,0.46],[36,0.47]]:
[>date_list2:=[[8,0.46],[13,0.47],[18,0.46], [23,0.48],[28,0.47],[33,0.49],[36,0.5]]:
[>date_list3:=[[8,0.49],[13,0.5],[18,0.49], [23,0.52],[28,0.51],[33,0.53],[36,0.55]]:
[>pointplot(date_list1, date_list2, date_list3, style=line, color=bluek, thickness=3);
КПЗ
1
2
3
Выводы:
Характер зависимостей однотипен, близок к линейному на интервале значений мелющих тел от 8 до 20 мм с одним угловым коэффициентом, а на интервале от 20 до 32 мм с другим угловым коэффициентом. Линейный характер зависимостей показывает, что увеличение размера мелющих тел не оказывает существенного влияния на плотность упаковки, а следовательно мало влияет на количество точечных контактов, а значит и на интенсивность процесса помола. Коэффициент полезного заполнения выше у бочкообразных камер меньшего типоразмера (кривая 1), как и в случае рабочей камеры цилиндрической формы, это объясняется тем, что в пристенной зоне пустоты имеют больший объем из-за того, что радиус кривизны камеры меньше. Следует отметить, что для камер, имеющих одинаковый объем, внутренняя поверхность камер бочкообразной формы больше, чем у цилиндрических. За счет того, что у камеры бочкообразной формы имеется кривизна поверхности не только в поперечном сечении, но и в продольном, что способствует увеличению объема пустот в пристенной зоне. И как следствие коэффициент полезного заполнения выше почти в 2 раза. На основании изложенного можно считать, что бочкообразная форма более технологична для процесса помола.
