Многокритериальная оптимизация геометрии щелевого фильтра для очистки жидкостей Текст научной статьи по специальности «Компьютерные и информационные науки»

НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА

НАУКА и ОБРАЗОВАНИЕ

Эл № ФС77 • 48211. Государственная регистрация №0421200025. ISSN 1994-0408

электронный научно-технический журнал

Многокритериальная оптимизация геометрии щелевого фильтра для очистки жидкостей # 02, февраль 2013 Б01: 10.7463/0213.0539055

Карпенко А. П., Мухлисуллина Д. Т., Цветков А. А.

УДК 519.6

Россия, МГТУ им. Н.Э. Баумана [email protected] [email protected], [email protected]

Введение

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

В МГТУ им. Н.Э. Баумана под руководством профессора Н.Н. Зубкова разработан новый метод механической обработки - деформирующее резание. Метод позволяет с высокой производительностью получать полимерные щелевые фильтрующие структуры с регулируемой шириной щелей от 10 мкм

до нескольких мм [2]. Исследованию различных аспектов изготовления фильтрующих структур с помощью метода деформирующего резания посвящены работы [1-4].

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

- живое сечение (отношение суммарной площади щелей к площади поверхности нарезанной части фильтра);

- жесткость фильтра, обусловленная раскрытием щелей вследствие осевой нагрузки на фильтроэлемент;

- прочность фильтра, определяемая напряжениями, возникающими в материале фильтра вследствие радиальной нагрузки на него.

Целью работы является трехкритериальная оптимизация геометрии щелевого фильтра, изготовленного методом деформирующего резания.

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

В качестве метода решения задачи многокритериальной оптимизации (МКО-задачи) использован развиваемый авторами интерактивный метод решения МКО-задач, называемый методом РЯЕЕ (РЯЕЕвгепсе) [5]. Каждая итерация метода включает в себя фазу анализа, выполняемую лицом, принимающим решения (ЛПР), и фазу расчетов, выполняемую системой многокритериальной оптимизации (МКО-системой).

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

непосредственно выполнять оценку предлагаемых МКО-системой альтернатив (например, в терминах «лучше», «хуже», «одинаково»). Указанные оценки МКО-система интерпретирует как значения, так называемой, функции предпочтений ЛПР [5]. Метод требует аппроксимации данной функции, для чего в работе использован подход на основе нейронных сетей [5].

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

1. Постановка задачи

Математическую модель щелевого фильтра определяют следующие его параметры (рисунок 1):

• О - внешний диаметр щелевого фильтра;

• Ь - длина фильтра;

• И - толщина стенки фильтра;

• О - длина необработанной части фильтра;

• N - число винтовых рядов щелей;

• Ь - ширина щелей;

• О - угол наклона винтовых рядов щелей (рисунок 1);

• Яи - радиус траектории инструмента, формирующего сквозную щель (рисунок 2);

• 8о - осевой шаг щелей (рисунок 1).

Рисунок 1 - К модели фильтра

Взаимосвязи указанных параметров фильтра определяет следующая последовательность действий.

1) Задав численные значения параметров с, Яи, 8о, рассчитываем шаг

'С 180

витков спирали £ = пБ

/ V180у

2) Вычисляем высоту спирали (длина участка трубной заготовки со сквозными щелями) Ис = Ь - 2G.

3) Определяем число щелей на одной спирали пщ = И с/ 8о.

4) Вычисляем требуемое расстояние А между осью вращения заготовки и осью вращения инструмента по формуле А = Б/2 + Яи - г, где г = к +1, мм -максимальная глубина внедрения инструмента, заданная из технологических соображений [1].

5) Определяем длину прорезаемой щели на внутренней поверхности заготовки 1внщ = 2рг, где г = (Б/2 - к) - внутренний радиус заготовки;

гг2 + А2 - Д2Л

р = агееоБ

2 Аг

центральный угол, на который опирается дуга /

вн.щ

(рисунок 2). Это сделано в допущении о независимости длины щели от режима резания (встречное или попутное) [6].

Рисунок 2 - К расчету длины щели на внутренней поверхности фильтра: 1 - профиль трубной заготовки фильтра, 2 - траектория инструмента

В качестве компонентов вектора варьируемых параметров X рассматриваем - = о - угол наклона винтовой линии рядов щелей, град;

х2 = Бо - шаг щелей, мм; х3 = Яи - радиус вращения инструмента, мм. Допустимая область варьируемых параметров, полагаем, имеет вид Ох = {Х 130 <-<80; 1 <х2<3; 20 <х3 <50}.

Используем три критерия оптимальности: живое сечение ф1( X); жесткость ф2( X), мм; прочность ф3( X) МПа.

1) Живое сечение равно отношению суммарной площади щелей а1 к площади нарезанной части поверхности фильтра сг2:

Говоря более строго, живое сечение подлежит максимизации. Однако для единообразия записи в формуле (1) используем знак минус и решаем задачу минимизации функции ф1 (X).

2) Жесткость фильтра определяем относительным увеличением ширины щелей (их раскрытием) в результате растяжения фильтра осевой силой F, приложенной к одному из торцов фильтроэлемента при закреплении другого торца:

Ф2(X) = AL^ min. (2)

Hc X

Здесь AL - удлинение фильтра, вызванное силой F = 100 Н.

3) Прочность фильтра. В качестве данного критерия используем максимальное напряжение, возникающее в материале фильтра при перепаде давления, равном AP:

ф3(X) ^ min

X

(3)

Ставим трехкритериальную задачу оптимизации фильтра: найти значения его варьируемых параметров X е Эх, которые, в пределе, обеспечивают выполнение всех условий (1) - (3).

Полагаем, что фильтроэлемент изготовлен из полиэтилена низкого давления марки ПЭ80 (ГОСТ 18599-2001). Основные свойства материала ПЭ80 представлены в таблице 1 [1].

Таблица 1 - Основные свойства материала ПЭ80

Свойство Значение

Модуль упругости, МПа 800

Коэффициент Пуассона 0,43

Предел прочности при растяжении, МПа 30

2. Используемый метод решения задачи многокритериальной

оптимизации

Пусть X = (x1,x2,...,xn) - вектор варьируемых параметров задачи. Область допустимых значений вектора X образует множество DX (которое, в общем случае, формирует некоторый набор ограничивающих функций). Векторный критерий оптимальности Ф(X) = (ф1 (X), ф2 (X),..., фт (X)) определен на множестве DX и ЛПР стремится минимизировать значение каждого из частных критериев оптимальности ф1 (X), i е [1: т], что условно записываем в виде

min Ф(Х) = Ф(Х*), Х е DX. (4)

X

В нашем случае n = т = 3.

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

Обозначим р(X, Л) операцию скалярной свертки частных критериев оптимальности, где Л = Л2,..., Äm) еDЛ - вектор весовых множителей; DA -множество допустимых значений этого вектора. При каждом фиксированном векторе Л е Dk метод скалярной свертки сводит решение задачи (4) к решению однокритериальной задачи глобальной условной оптимизации (ОКО-задачи)

min р(Х, Л) = р(Х*, Л), Х е DX. (5)

X

Обозначим D(b множество достижимости задачи, то есть множество, в которое векторная целевая функция Ф^) отображает область DX. Если при каждом Л е D^^ решение задачи (5) единственно, то это решение ставит в http://technomag.bmstu.ru/doc/539055.html 227

соответствие вектору Л единственный вектор X [5]. Данное обстоятельство позволяет полагать, что функция у предпочтений ЛПР определена не на множестве DX, а на множестве . В результате МКО-задача сводится к задаче выбора вектора Л* е БЛ такого, что

тах у (Л) = у(Л ), Л е БЛ.

Л

Величину у полагаем лингвистической переменной со значениями от «очень - очень плохо» до «отлично». Ядро нечеткой переменной у

о

обозначим у, так что ее значению «очень - очень плохо» соответствует

о о

у = 1, а значению «отлично» - у = 9.

Таким образом, МКО-задачу сводим к задаче отыскания вектора Л* е DЛ, обеспечивающего максимальное значение целочисленной функции

у (Л):

0 0 ° таху(Л) = у(Л*) = у\ ЛеDЛ. (6)

Л

Детальная схема используемого метода решения задачи (6) представлена, например, в работе [5].

Предварительное исследование задачи, поставленной в п. 1, показало, что в условиях постановки задачи множество БФ является выпуклым. Это дает возможность в качестве функции ((X, Л) использовать аддитивную

m

свертку (X). Известно, что решения, доставляющие минимум этой

г=1

свертки, принадлежат множеству Парето задачи.

3. Программное обеспечение 3.1. Приложение FILTER-SW. Параметризованная модель щелевого фильтра (п.1) разработана в среде программной системы SolidWorks. С целью

исследования жесткости и прочности фильтра разработано программное приложение FILTER-SW, позволяющее программно управлять этой моделью. Разработка выполнена с использованием интегрированной среды разработки VisualStudio 2010 на языке программирования С++. Объем программного кода составляет ~1300 строк. Приложение использует интерфейс прикладного программирования (Application Programming Interface, API) системы SolidWorks, предназначенный для управления набором функциональных средств SolidWorks.

В качестве конечных элементов (КЭ) конечно-элементной сетки используем параболические тетраэдральные элементы, которые в отличие от линейных позволяют получить более точную аппроксимацию модели [5]. Программное приложение FILTER-SW позволяет не задавать максимальный размер ребра конечного элемента smax В этом случае приложение использует

рассчитанное им умолчательное значение smax = s'^^, которое определяется в

пакете CosmosWorks, исходя из геометрических параметров модели. Пакет CosmosWorks входит в состав SolidWorks и предназначен для выполнения конечно-элементного анализа.

3.2. Приложение PREF-FILTER-SW. Рассмотренный в п.2 метод решения МКО-задачи реализован в виде программного приложения PREF-FILTER-SW. Реализация выполнена с использованием интегрированной среды разработки VisualStudio 2010 на языке программирования C++. Объем программного кода составляет ~2900 строк.

Для аппроксимации функции предпочтений в программе PREF-FILTER-SW использован многослойный персептрон с семью нейронами в скрытом слое. Многослойный персептрон, по сути, строит указанную аппроксимирующую функцию в виде

i=0 V J=0

\\

где ai, Д j - параметры, вычисляемые в процессе обучения нейронной сети; /i("), ЛО) " сигмоидальная и линейная функции соответственно; H - число нейронов в скрытом слое нейронной сети (в нашем случае H=7). Программная реализация нейронных сетей выполнена на основе свободно распространяемого пакета FANN (Fast Artificial Neural Network Library) [8].

Для решения ОКО-задач, возникающих в методе PREF, использован свободно распространяемый пакет NLopt (NonLinear Optimization), поддерживающий широкий класс методов многомерной локальной и глобальной оптимизации [9]. В частности, в качестве метода глобальной условной оптимизации для решения задач (5) и (6) нами использован метод мультистарта в комбинации с прямым методом многомерной локальной условной оптимизации COBYLA (Constrained Optimization BY Linear Approximations), реализация которого имеется в указанном пакете.

4. Вычислительные эксперименты

Вычислительные эксперименты выполнены для модели реального фильтра, использованного в работе [1]. Параметры модели представлены в таблице 2. Эксперименты выполнены на персональном компьютере, имеющем следующие значения основных параметров: процессор - Intel Core(TM)2 Quad CPU Q9550; 2,83 ГГц; оперативная память - 4,00 Гбайт.

Таблица 2 - Параметры модели

Параметр Значение

Наружный диаметр заготовки (D), мм 50

Толщина стенки заготовки (К), мм 4

Число рядов щелей (N) 6

Радиус вращения инструмента (Яи), мм 95

Ширина щелей (Ъ), мм 0,1

Шаг щелей (8о), мм 1

Угол наклона винтовой линии рядов щелей (О), градус 43

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

Стандартной инженерной практикой при выборе размера КЭ является исследование энергии деформации при варьировании размера КЭ. При этом размеры КЭ уменьшают до тех пор, пока энергия деформации практически престанет меняется. Результаты такого сорта экспериментов, выполненных с помощью конечно-элементной SolidWorks-модели фильтра, представлены на рисунках 3-5. Исследование выполнено при варьировании максимального размера £тах ребра КЭ в диапазоне от 5 мм до 0,5 мм.

£ 4

® 3.5

ГО

о.

х 10

3"

ГО ^

о. о -8-да ч сх

о.

О)

2.5 2 1.5 1

0.5 0

! !

Т 1 • : | т*

т* | < 41 • 1 • •

и

*

* гТт •

• •

¥ • тт и».......

0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 Максимальный размер конечного элемента, мм

Рисунок 3 - Максимальная энергия деформации фильтра: SolidWorks-модель; осевое растяжение фильтра силой F=100 Н

.х 10

ф ш ш

го 7 .

с т

=г

го

О -&

ч 4 от

&3

п; го п: _о с

J 1 а

•

Т?

•

• • • II •

т* • •

1 • * •

•т •

гтТТ1 Г?!

0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 Максимальный размер конечного элемента, мм

Рисунок 4 - Максимальная энергия деформации фильтра: SolidWorks-модель;

перепад давления АР = 150 кПа

CL

о; ф

Cl

m

10'

10"

-1

0.5

1.5 2 2.5 3 3.5 4 4.5 5 Максимальный размер конечного элемента, мм

Рисунок 5 - Время расчета: SolidWorks-модель; осевое растяжение фильтра силой F=100 H; перепад давления AP = 150 кПа

Представленные результаты исследования позволяют сделать вывод о том, что при размерах конечных элементов smax, меньших 1 мм, энергия

деформации, обусловленная, как осевым растяжением фильтра, так и перепадом давления, меняется незначительно. В то же время, при такой величине smax вычислительные затраты оказываются приемлемыми. На этом основании, а также исходя из принципа гарантированного результата, дальнейшие исследования выполнены при smax = 1 мм.

Из экспериментов, результаты которых представлены выше, вытекает, что максимальные напряжения при указанных деформациях фильтра возникают у оснований щелей (рисунок 6).

Рисунок 6 - Распределение напряжений в фильтре

4.2. Трехкритериальная оптимизация фильтра. В качестве лица, принимающего решения, выступал профессор МГТУ им. Н.Э. Баумана Н.Н. Зубков. Число «разгонных» решений, которые предоставляются ЛПР для оценки на первой итерации, было принято равным пяти [5]. На рисунке 7 в трех разных проекциях изображены решения, предлагаемые ЛПР для оценки. Задача была решена за девять итераций диалога с ЛПР. Полученному решению соответствует точка с номером 13 на рисунке 10. Соответствующие значения варьируемых параметров и критериев оптимальности равны: угол наклона винтовой линии рядов щелей x1 = ю = 59 град; шаг щелей x2 = S0 = 1 мм; радиус вращения инструмента x3 = Ru = 49 мм; живое сечение

ф1 = - 0,093; жесткость фильтра ф2 = 0,012 мм; прочность фильтра ф3 = 3,675 МПа.

6.5 6 5.5 5 4.5 4 3.5 3 2.5

i

i 03

О6....................

о® OS 2 .......л.

11 л 12

13 • u Cu 9 ö

Ol С 1

О4

i COO

1 -009 -0 08 -007 -006 -005 -004 -0 03 -0 02 -001 fl

-0.01

-0.02

-0.03

-0.04

-0.05

-0.06

-0.07

-0.08

-0.09

-0.1

а)

7

............и

01Û

012

6 _..............iff )3.................

o2 О"

...............^ 5 3)11..............

0801

• 13

0.005

0.01

0.015

0.02

0.025

0.03

2 5-

MM

б)

в)

а) - проекция в осях ф1 и ф3; б) - проекция в осях ф2 и ф1; в) - проекция в осях ф2 и ф3

Рисунок 7 - Аппроксимация фронта Парето для трехкритериальной задачи

оптимизации геометрии фильтра

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

Заключение

Основные результаты работы заключаются в следующем.

1) Разработана математическая модель щелевого фильтра, изготавливаемого методом деформирующего резания.

2) Поставлена трехкритериальная задача оптимизации геометрии этого фильтра.

3) Разработана параметризованная SolidWorks-модель фильтра.

4) С использованием интегрированной среды разработки VisualStudio 2010 и языка программирования C++ разработано два программных приложения. Приложение FILTER-SW позволяет программно управлять SolidWorks-моделью фильтра. Приложение PREF-FILTER-SW реализует развиваемый авторами работы интерактивный метод решения многокритериальных задач оптимизации, основанный на нейросетевой аппроксимации функции предпочтений ЛПР.

5) C использованием SolidWorks-модели фильтра и указанного программного обеспечения выполнена трехкритериальная оптимизация геометрических параметров фильтра. Опыт решения задачи показал эффективность принятых алгоритмических и программных решений.

В развитие работы планируются эксперименты по многокритериальной оптимизации геометрии фильтра с использованием других развиваемых авторами методов многокритериальной оптимизации [10 ,11]. Кроме того, для повышения скорости вычислений авторы планируют параллельную реализацию рассмотренного и других методов многокритериальной оптимизации.

Авторы выражают глубокую благодарность Зубкову Н.Н. и Овчинникову В.А. за ценные советы и помощь в проведении вычислительного эксперимента.

Список литературы

1. Слепцов А.Д. Разработка и исследование метода получения щелевых фильтрующих структур: дис. ... канд. тех. наук: 05.02.07, 05.02.08 / Слепцов Аркадий Дмитриевич. МГТУ им. Н.Э. Баумана, 2011. 179 с.

2. Зубков Н.Н. Совмещение процессов резания и обработки давлением в новом методе формообразования развитых поверхностных структур // Кузнечно-штамповочное производство. Обработка материалов давлением. 2002. № 10. С.17-20, С. 29-34.

3. Зубков Н.Н., Слепцов А.Д. Получение микросеток и проницаемых щелевых труб механической обработкой // Известия ВУЗов. Машиностроение. 2007. № 3. С. 56-60.

4. Зубков Н.Н., Слепцов А.Д. Получение полимерных щелевых фильтрующих труб методом деформирующего резания // Вестник машиностроения. 2010. № 12. С. 51-53.

5. Карпенко А.П., Мухлисуллина Д.Т., Овчинников В.А. Нейросетевая аппроксимация функции предпочтений лица, принимающего решения, в задаче многокритериальной оптимизации // Информационные технологии. 2010. № 10. С. 2-9.

6. Иванов В.С., Зубков Н.Н., Иванов Д.В. Параметрический синтез технологической схемы при получении щелевых структур вращающимся лезвийным инструментом // Технология машиностроения. 2007. № 12. С. 16-19.

7. Справка по SolidWorks . Режим доступа: http://help.solidworks.com/ (дата обращения 10.10.2012).

8. FANN - Fast Artificial Neural Network Library. Режим доступа: http://leenissen.dk/fann/wp/ (дата обращения 10.10.2012).

9. NLopt. NonLinear Optimization. Режим доступа: http://ab-initio.mit.edu/wiki/index.php/NLopt (дата обращения 10.10.2012).

10.Карпенко А.П., Моор Д.А., Мухлисуллина Д.Т. Многокритериальная оптимизация на основе нечеткой аппроксимации функции предпочтений лица, принимающего решения // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2010. № 1. Режим доступа: http : //technomag. edu. ru/doc/135375. html (дата обращения 10.10.2012).

11.Карпенко А.П., Моор Д.А., Мухлисуллина Д.Т. Многокритериальная оптимизация на основе нейро-нечеткой аппроксимации функции предпочтений лица, принимающего решения // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2010. № 6. Режим доступа: http://technomag.edu.ru/doc/143964.html (дата обращения 10.10.2012).

SCIENTIFIC PERIODICAL OF THE RAIJMAN MS TU

SCIENCE and EDUCATION

EL № FS77 - 48211. №0421200025. ISSN 1994-0408

electronic scientific and technical journal

Multiobjective optimization of water slotted filter geometry

# 02, February 2013

DOI: 10.7463/0213.0539055

Karpenko A.P., D.T. Mukhlisullina, Cvetkov A.A.

Bauman Moscow State Technical University, 105005, Moscow, Russian Federation

[email protected] [email protected], [email protected]

The article deals with a three-objective optimization of slotted filter geometry produced by an innovative method of deforming cutting. This method was developed by Professor N.N. Zoubkov from BMSTU. The multiobjective problem is solved by using the interactive PREF method proposed by the authors. This method is based on the assumption that a decision maker (DM) is able to provide, in some form, preference information related to various solutions proposed by multiobjective optimization software. The DM's utility function is approximated by a neural network. This work presents a three-objective optimization problem statement of filter design, a short description of an interactive method and its software implementation along with results of computational experiments.

Publications with keywords: multiobjective optimization, slot filter, decision maker's utility function

Publications with words: multiobjective optimization, slot filter, decision maker's utility function

References

1. Sleptsov A.D. Razrabotka i issledovanie metodapolucheniia shchelevykh fil'truiushchikh struktur. Kand. diss. [Development and research of the method of obtaining of slotted filter structures. Cand. diss.]. Moscow, Bauman MSTU, 2011. 179 p.

2. Zubkov N.N. Sovmeshchenie protsessov rezaniia i obrabotki davleniem v novom metode formoobrazovaniia razvitykh poverkhnostnykh struktur [Combining processes of cutting and processing pressure in the new method of forming of developed surface structures]. Kuznechno-

shtampovoe proizvodstvo. Obrabotka materialov davleniem [Forging and stamping production. Processing of materials by pressure], 2002, no. 10, pp.17-20, pp. 29-34.

3. Zubkov N.N., Sleptsov A.D. Poluchenie mikrosetok i pronitsaemykh shchelevykh trub mekhanicheskoi obrabotkoi [Getting micro grids and permeable slotted pipe by machining]. Izvestiia VUZov. Mashinostroenie [Proceedings of Higher Educational Institutions. Machine Building], 2007, no. 3, pp. 56-60.

4. Zubkov N.N., Sleptsov A.D. Poluchenie polimernykh shchelevykh fil'truiushchikh trub metodom deformiruiushchego rezaniia [Production of slotted polymer filter tubes by deformational cutting]. Vestnikmashinostroeniia, 2010, no. 12, pp. 51-53. (Trans. version: Russian Engineering Research, 2010, vol. 30, no. 12, pp. 1231-1233.).

5. Karpenko A.P., Mukhlisullina D.T., Ovchinnikov V.A. Neirosetevaia approksimatsiia funktsii predpochtenii litsa, prinimaiushchego resheniia, v zadache mnogokriterial'noi optimizatsii [Neural network approximation of decisions maker's utility function in multicriteria optimization problem]. Informatsionnye tekhnologii, 2010, no. 10, pp. 2-9.

6. Ivanov V.S., Zubkov N.N., Ivanov D.V. Parametricheskii sintez tekhnologicheskoi skhemy pri poluchenii shchelevykh struktur vrashchaiushchimsia lezviinym instrumentom [Parametric synthesis of the technological scheme at obtaining slotted structure using rotating blade tool]. Tekhnologiia mashinostroeniia, 2007, no. 12, pp. 16-19.

7. Spravkapo SolidWorks [Help for SolidWorks]. Available at: http://help.solidworks.com/ , accessed 10.10.2012.

8. FANN - Fast Artificial Neural Network Library. Available at: http://leenissen.dk/fann/wp/ , accessed 10.10.2012.

9. NLopt. NonLinear Optimization. Available at: http://ab-initio.mit.edu/wiki/index.php/NLopt , accessed 10.10.2012.

10. Karpenko A.P., Moor D.A., Mukhlisullina D.T. Mnogokriterial'naia optimizatsiia na osnove nechetkoi approksimatsii funktsii predpochtenii litsa, prinimaiushchego resheniia [Multi-criteria optimization based on fuzzy approximation of the preferences function of a decision maker]. Nauka i obrazovanie MGTUim. N.E. Baumana [Science and Education of the Bauman MSTU], 2010, no. 1. Available at: http://technomag.edu.ru/doc/135375.html , accessed 10.10.2012.

11. Karpenko A.P., Moor D.A., Mukhlisullina D.T. Mnogokriterial'naia optimizatsiia na osnove neiro-nechetkoi approksimatsii funktsii predpochtenii litsa, prinimaiushchego resheniia [Multicriteria optimization based on neuro-fuzzy approximation of decision maker's utility function]. Nauka i obrazovanie MGTU im. N.E. Baumana [Science and Education of the Bauman MSTU], 2010, no. 6. Available at: http://technomag.edu.ru/doc/143964.html , accessed 10.10.2012.