Научная статья на тему 'Моделирование параметров технологического процесса обработки с использованием метода Монте-Карло'

Моделирование параметров технологического процесса обработки с использованием метода Монте-Карло Текст научной статьи по специальности «Математика»

CC BY
195
57
i Надоели баннеры? Вы всегда можете отключить рекламу.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Погонин Анатолий Алексеевич, Чепчуров Михаил Сергеевич, Белянкина Ольга Владимировна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Моделирование параметров технологического процесса обработки с использованием метода Монте-Карло»

© A.A. Погонин, М.С. Чепчуров, О.В. Белянкина, 2003

УЛК 681.06 (07)

A.A. Погонин, М.С. Чепчуров, О.В. Белянкина

МOЛEЛИPOBAHИE ПAPAМETPOB TEХHOЛOГИЧECKOГO ПPOUECCA OБPAБOTKИ С ИCПOЛЬЗOBAHИEМ МETOЛA МOHTE-KAPЛO

При решении оптимизационных задач постоянно возникает необходимость поиска максимума многоэкстремальной функции, простое решение дифференциального уравнения, то есть нахождение экстремума не приводит к положительным результатам, как это показано в примере из [4]. Для поиска глобального оптимума можно воспользоваться различными методами, одним из таких методов является "Метод Монте-Карло", т.е. статистический метод случайного поиска. Этот метод позволяет вычислить с относительной точностью искомое значение функции. Вычисление функции производится путем многократного ввода случайных значений, лежащих в заданном интервале. Такие значения создаются специальной функцией генерации случайного числа. Чем больше генерируется значений, тем точнее результат.

Определение оптимальных значений параметров технологического процесса при обработке неста-

R:= SOO

P := 5

f(a ,ф) := (r + p) • cos

a

2 • 9O

x := 6OO І

tan

i := І.. 3O k := І.. 45 aj := 3 • i фк := k

1 9O у

M(i,k) := f(i>fe)

а0 :=0 ф :=0 :=0 Начальные значения

ОИЮІЬЬ 1 п := 100000 і := 1.. п

аі :=гп^90) ф := т^45)

Ц := і( аьф) Максимум=та:(Б) Максимум= 0.25

Рис. 1. Определение максимума методом Монте-Карло

ционарным станочным модулем задача, связанная со значительными временными материальными затратами. Проведение эксперимента по обработке крупногабаритных изделий требует наличия оборудования, в конструкцию которого входят эти изделия, сам эксперимент проводится в тяжелейших условиях, для его осуществления расходуется много режущего инструмента. Аналитическое решение уравнения с несколькими неизвестными, требует применения мощных средств вычислительной техники и больших затрат времени. Допустим, требуется определить оптимальное расположение линии установки инструмента в плоскости опорного ролика? Согласно уравнения приведенного в [3], очевидно, что проведение натурных измерений или перебор возможных вариантов, связан с затратами времени и средств. Для реализации предлагаемого способа, воспользуемся средствами MathCAD [1, 2].

Как видно из графика, приведенного на рис. 1, искомая функция с двумя аргументами имеет множество экстремумов, то есть максимальных значений. В приведенном примере перед началом поиска производится задание индекса первого элемента массива с вычисляемыми значениями функций, затем происходит задание начальных значений аргументов, числа случайных значений аргументов. В качестве аргумента функции генерации случайного числа задается максимальное значение аргумента (то есть генерируется число в диапазоне от 0 до max). Затем происходит вычисление значений вектора функции, нахождение его максимального значения, это значение и будет являться глобальным оптимумом. В принципе при поиске оптимума при одном аргументе его значение можно найти, решив исходное уравнение. Для решения уравнения с двумя неизвестными приходится пользоваться численным метом. Вариант реализации алгоритма решения уравнения с двумя неизвестными численным методом с заданной точностью (TOL) показан на рис.2.

Решение уравнения задано в функции пользователя, передаваемыми параметрами являются максимальные значение аргументов, найденное макси-

Создание функции пользователя по нахождению корней уравнения

Корень(н, Xmax, YmaxR, p)

TOL ^ O1 X ^ O

while X < Xmax

Y ^ І

while Y < Ymax V2 ^ Y V, ^ X

return V if (R + p) • cos

Корень (Максимум , 90,45, Я, р)

Рис. 2. Численный метод решения уравнения

-x

X • п

н

-x

2 • !8O

tan Y

!8O

V

мальное значение искомой функции. Возвращаемым значением является вектор, содержащий значения аргументов. После получения результатов желательно выполнить проверку, как это показано в примере. Программу по приближенному нахождению корней уравнения можно использовать и в качестве самостоятельного приложения, количество неизвестных в уравнении можно увеличить, увеличив количество циклов в программе, описывающей функцию пользователя. Предложенный способ расчета глобального оптимума методом Монте-Карло

подходит не только для решения аналитических задач, но и логических, например, поиска оптимального управленческого решения. Методом можно пользоваться при моделировании протекания различных технологических процессов.

Как было отмечено выше, рассмотренный метод относится к статическим, то есть существует вероятность того, что решение будет неточным. Это заключение справедливо при малом начальном значении испытаний, при росте количества испытаний вероятность получения точного решения возрастает.

СПИСО K ЛИTEPATУPЫ

1. Погонпн А.А., Чепчуров М.С.

Инженерные расчеты в MathCAD 7.0 Prof. Учебное пособие. - Белгород Изд-во БелГТАСМ, 2001,-96 с.

2. Очков В.Ф. Mathcad 7 Pro для студентов и инженеров. -М.: КомпьютерПресс, 1998. -384 с. -ил.

3. Погонпн А.А. Восстановление точности крупногабаритных де-

талей машин без демонтажа в процессе эксплуатации. Горный информационно-аналитический бюллетень № 11, 2001, 0236-1493. -М. Из-

дательство МГГУ.

KOPOTKO ОБ ABTOPAX

Погонпн Анатолий Алексеевич — профессор, кандидат технических наук, Белгородская Государственная Технологическая Академия Строительных Материалов

Чепчуров Михаил Сергеевич — старший преподаватель, Белгородская Государственная Технологическая Академия Строительных Материалов

Белянкина Ольга Владимировна - ассистент, Московский государственный горный университет

© A.A. Погонин, М.С. Чепчуров, 2003

УЛK 681.06 (07)

A.A. Погонин, М.С. Чепчуров

BЫБOP ПAPAМETPOB УПPABЛEHИЯ TEXHOЛOГИЧECKИМ ПPOUECCOМ ПPИ OБPAБOTKE HECTAUИOHAPHЫМИ CTAHOЧHЫМИ МОЛУЛЯМИ

бравотка крупногабаритных изделий нестационар-■ымивтаночными модулями в настоящее время ста-ла^всновным способом восстановления их работоспособности. В последнее время разработаны новые совершенные конструкции встраиваемых и приставных станков, таких, как, например, описываемый в [4], определены параметры технологического процесса обработки.

Обработка изделия нестационарным станочным модуле требует значительных временных затрат, измеряемых сутками, или рабочими сменами. Все это время в зоне обработки должен находиться обслуживающий персонал, осуществляющий контроль за ходом технологического процесса, состоянием оборудования и режущего инструмента. Следовательно в себестоимости обработки изделия доля его заработной платы значительна [3]. А учитывая совре-

менные тенденции в развитии экономики процент себестоимости учитывающий заработную плату и отчисления на социальные нужды будет возрастать, увеличивая тем самым себестоимость. Решение этой проблемы видится в автоматизации технологического процесса, а также в автоматическом контроле за стоянием оборудования и режущего инструмента. Подобное техническое решение позволит также, вывести персонал из неблагоприятных условий работы.

Процесс перевода нестационарного станочного модуля на автоматический контроль и управление связан с дополнительными затратами на проектирование системы управления. Но все же главным является выбор параметров системы управления и оценка ее эффективности.

Началом работы на проектом является определение способа обработки, при обработка бандажей и опорных роликов цементных печей таковым является точение.

Прежде всего следует определиться с целью автоматизации процесса обработки, согласно рекомендациям приведенным в [1]: повышение производительности труда, снижение себестоимости, повышение точности, повышение качества поверхности,

i Надоели баннеры? Вы всегда можете отключить рекламу.