Научная статья на тему 'Модель прогнозирования трудоемкости изготовления мотогондолы авиационного гражданского ГТД в АСТПП'

Модель прогнозирования трудоемкости изготовления мотогондолы авиационного гражданского ГТД в АСТПП Текст научной статьи по специальности «Механика и машиностроение»

CC BY
356
179
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
КОНКУРЕНТОСПОСОБНОСТЬ / ТЕХНОЛОГИЧЕСКОЕ ПРОГНОЗИРОВАНИЕ / АСТПП / МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ / ТРУДОЕМКОСТЬ / ПОЛИМЕРНЫЕ КОМПОЗИЦИОННЫЕ МАТЕРИАЛЫ / COMPETITIVENESS / TECHNOLOGICAL FORECASTING / CAM / MATHEMATICAL MODELING / LABOR INPUT / POLYMER COMPOSITE MATERIALS

Аннотация научной статьи по механике и машиностроению, автор научной работы — Юрьев В. Л., Иркабаев И. М., Сухова Н. А., Кадырова Ж. М.

В рамках методики определения трудоемкости изготовления авиационных гражданских ГТД в условиях кооперации предлагается математическая модель определения трудоемкости изготовления деталей и сборочных единиц из полимерных композиционных материалов (ПКМ) мотогондолы ГТД.

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

Forecasting model the complexity of manufacturing nacelles civil Aviation gas turbine engine in CAM

As part of the development of a methodology of determining the complexity of producing Aviation gas turbine engine in cooperation mathematical model determining the complexity of manufacturing parts and assembly units of polymeric composite materials Aviation gas turbine engine nacelle.

Текст научной работы на тему «Модель прогнозирования трудоемкости изготовления мотогондолы авиационного гражданского ГТД в АСТПП»

ISSN 1992-6502 (P ri nt)_

2015. Т. 19, № 1 (67).С. 37-44

Ъыьмт QjrAQnQj

ISSN 2225-2789 (Online) http://journal.ugatu.ac.ru

УДК 004.65

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

вастпп

в. л. Юрьев\ и. м. Иркабаев2, н. а. Сухова3, ж. м. Кадырова4

4 prep08@mail.ru

1-4 ОАО «Институт технологии и организации производства» (НИИТ) 1,3,4 ФГБОУ ВПО «Уфимский государственный авиационный технический университет» (УГАТУ)

Поступила в редакцию 10 октября 2014 г.

Аннотация. В рамках методики определения трудоемкости изготовления авиационных гражданских газотурбинных двигателей (ГТД) в условиях кооперации предлагается математическая модель определения трудоемкости изготовления деталей и сборочных единиц из полимерных композиционных материалов (ПКМ) мотогондолы ГТД.

Ключевые слова: конкурентоспособность; технологическое прогнозирование; АСТПП; математическое моделирование; трудоемкость; полимерные композиционные материалы.

ВВЕДЕНИЕ

Актуальность данной публикации обусловлена государственной программой «Развитие авиационной промышленности на 2013-2025 годы», которая преследует цель создания высококонкурентной авиационной промышленности и закрепление ее позиции на мировом рынке в качестве третьего производителя по объемам выпуска авиационной техники. Данная программа имеет ряд подпрограмм, в т.ч. подпрограмму 1 «Самолетостроение» и подпрограмму 3 «Авиационное двигателестроение», в конечном итоге предполагающие выпуск конкурентоспособной продукции.

Цели и задачи. На стадии проектирования нового изделия для обеспечения в дальнейшем необходимого уровня его конкурентоспособности возникает необходимость определения основных показателей технологичности конструкции, а в частности трудоемкости изготовления. Согласно ГОСТ 14.201 «Обеспечение технологичности конструкции изделий. Общие требования» рекомендуемыми показателями технологичности конструкции изделия являются трудоемкость изготовления изделия, удельная материалоемкость изделия, себестоимость [1]. Для разработки и управления конкурентоспособными проектами в условиях автоматизированных систем технологической подготовки производ-

ства (АСТПП) целесообразно использовать актуальные на данный момент методы математического моделирования на основе функционального моделирования [2].

Научная новизна определяется построенными регрессионными моделями трудоемкости мотогондолы авиационного гражданского ГТД.

Практическая полезность обуславливается тем, что разработанная модель прогнозирования трудоемкости мотогондолы авиационного гражданского ГТД апробирована в рамках методики комплексной оценки технологичности изготовления ДСЕ узлов мотогондолы из полимерных композиционных материалов.

ТЕХНОЛОГИЧЕСКОЕ ПРОГНОЗИРОВАНИЕ

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

ках работы управления инвестиционными и инновационными проектами и технологического прогнозирования. Функции НТПП в системе

инновационной подготовки производства в ОПК и авиационной промышленности представлены на блок-схеме на рис. 1 [4].

Рис. 1. Блок-схема функции НТПП в системе инновационной подготовки производства в ОПК и авиационной промышленности

МОДЕЛЬ РАСЧЕТА ТРУДОЕМКОСТИ УЗЛОВ МОТОГОНДОЛЫ

Алгоритм построения модели прогнозирования трудоемкости изготовления заготовок ДСЕ мотогондолы авиационных гражданских ГТД из полимерных композиционных материалов (ПКМ) (рис. 2).

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

• без усиливающих и сотовых элементов (класс 1);

• с усиливающими и сотовыми элементами (класс 2).

Трудоемкость изготовления деталей из ПКМ для обоих классов может быть определена следующим образом:

Т.

Т „„. + Т.

Рис. 2. Алгоритм построения модели прогнозирования трудоемкости изготовления узлов мотогондолы авиационных гражданских ГТД

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

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

Т — т К К К

з.пкм база 1 ксл сл а'

где Тбаза;- базовая трудоемкость изготовления

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

Значение Кксл определялосьэкспертно в зависимости от группы конструктивной сложности детали из ПКМ (табл. 1, нумерация групп прямо пропорциональна сложности деталей, от наипростейших в группе под номером 1 и самых сложных в группе под номером 6).

Таблица 1 Характеристика групп конструктивной сложности деталей из ПКМ

№ Характеристика группы

гр.

1. Плоские монолитные детали из ПКМ с небольши-моребрением, выступами, с небольшой кривизной

2. Монолитные незамкнутые кольцевые (секторные) оболочки с поверхностями второго порядка

3. Монолитные детали сложной формы, представляющие набор оболочек второго порядка

4. Кольцевые сэндвич-панели, корпуса с поверхностями второго порядка, имеющие элементы повышения жесткости

5. Незамкнутые кольцевые сэндвич-панели, имеющие поверхности второго порядка

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

Коэффициент сложности собственной конструкции (Ксл) ПКМ определяется по формуле

К = К К К

сл ус. сот.эл. наш.эл.'

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

Коэффициент, отражающий усиление конструкции (Кус) , определяется по формуле

Кус =

гтвыкл дт . гтвыкл дт

осн.сл. осн.сл. доп.С доп.с

д ВЫЮ, N

осн.сл. осн.сл.

где 5дСнКСл ■ - площадь выкладки основных слоев; Л?0 сн сл - количество основных слоев; З'доп1^ -площадь выкладки дополнительных слоев; Л?Д0 п ■ сл ■ - количество дополнительных слоев. Все значения взяты из конструкторской документации ДСЕ узлов мотогондолы.

Коэффициент наличия сотопакетов определяется по формуле

К„

д + д N

сотопак сот.

д.

где - площадь сотопакетов; - коли-

чество рядов сотопакетов; - площадь детали. Все значения взяты из конструкторской документации ДСЕ узлов мотогондолы.

Коэффициент, учитывающий наличие накладных элементов, определяется по формуле

К

д + д N

_ дет. накл. накл.

где - площадь накладных элементов;

- количество накладных элементов. Все значения взяты из конструкторской документации ДСЕ узлов мотогондолы.

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

Моделированию предшествовала проверка качества информации. К исходным данным предъявлялись следующие требования: сопоставимость, однородность, устойчивость (табл. 2).

Таблица 2 Проверка качества информации

Требования Обеспечение требований

к информации

Сопоставимость Все показатели, использованные в моделях, принадлежат к однородным объектам, при расчетах использовался единый инструментарий (редактор электронных таблиц MicrosoftExcel).

Однородность (на отсутствие аномальных значений) Все значения рядов прошли проверку на принадлежность к интервалу (х—3 о,х+3 о) , где х - среднее, о - стандартное отклонение

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

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

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

4 осн.слоев7

ружной поверхности детали (5"дет).

Базовая трудоемкость 1 мм2площади выкладываемой поверхности с учетом количества основных слоев является отражением затрат труда на саму структуру конструкций ПКМ и определяется как

Х

Т

база I

д.

выкл

осн.слоев

Функциональная зависимость базовой трудоемкости 1 мм2площади выкладываемой поверхности с учетом количества основных слоев

(Хбазш) от площади наружной поверхности детали представленных классов конструкций заго-товки^дет) также определяет критерий новизны настоящей публикации:

Х . = I

база I

В рамках построения регрессии был проведен графический анализ (рис. 3, 4) Хбазш-и площади наружной поверхности ДСЕ узлов мотогондолы, выбранных для построения модели в зависимости от класса конструкции заготовки.

На рис. 3, 4 явно прослеживается гиперболическая связь, возможно функция полинома первой, второй степени, логарифмической и показательной функций. Для проверки значимости данной комбинации функций от переменной 5"дет, эта комбинация была изначально добавлена как факторы в регрессии двух классов узлов ДСЕ мотогондолы (табл. 3). Все факторы про-

шли проверку по /-критерий Стьюдента с 5-процентной вероятностью.

Процедура отбора факторов в уравнение представляла собой шаговую регрессию. Сущность метода шаговой регрессии заключается в последовательном исключении факторов из уравнения регрессии, если факторы при проверке их значимости оказывались незначимыми. Первыми из регрессии удаляются те переменные, абсолютные критическое значение ^Стью-дента которых намного меньше критического, и так по очереди все остальные. В результате применения метода пошаговой регрессии остались переменные, которые оказались значимыми. Таким образом, в табл. 4 представлены регрессионные модели, которые реализованы с помощью табличного процессора Мюго80ЙЕхсе1, имеющего в своем распоряжении необходимые функции и механизмы.

X

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

100 90 80 70 „ 60 1 50 40 30 20 10

35

30

25

2 20

х115 10

я,

0,039

0,041

0,065

0,072

0,082

0,125

0,744

0,794

1,5

Рис. 3. График базовой трудоемкости 1 мм площади выкладываемой поверхности с учетом количества основных слоев от площади наружной поверхности детали

ДСЕ 1 класса

Хфакт Храсч

сС? сЛл ^ Л0 П<ЪЛ & пль П<ЪЬ „<£ ^ Л> & Л> ¿¡> ъ? ^ й? Су? ^ V V ^ ^ V ы

0

5

Рис. 4. График базовой трудоемкости 1 мм площади выкладываемой поверхности с учетом количества основных слоев от площади наружной поверхности детали ДСЕ 2 класса

Таблица 3

Проверка значимости переменных регрессии

Показатели значи мости Переменные регрессий Комбинация функций до применения метода пошаговой регрессии Комбинация функций, после применения метода пошаговой регрессии

1 класс 2 класс 1 класс 2 класс

^Стьюдента при переменных регрессии 1 0,52 0,73 — 6,48

1 б 1 0,67 0,22 2,27 —

С б дет 0,95 0,74 2,23 —

б 1 -0,16 0,66 -2,21 2,102

0,51 0,51 2,29 2,101

1 5дет| -0,91 1,43 — 4,54

t критическое 2,26 2,10 2,26 2,10

Значимость связи была проверена на основе показателей качества подбора уравнения (в случае сравнения нормативно-расчетной с трудоемкостью, рассчитанной по экономико-математической модели) и на основе показателей тесноты связи (в случае сравнения экспертно-технологической с трудоемкостью, рассчитанной по экономико-математической модели) в сравнении с критическими значениями (табл. 5).

Для оценки качества (точности) зависимостей рассмотрен коэффициент детерминации В2. Достаточно тесная связь наблюдается при условии, что коэффициент детерминации больше 0,8 [5, 6].Фактическое значение ^-критерия Фишера сравнивается с критическим табличным значением ^крит. (а, к1, к2) при заданном уровне значимости а и степенях свободы к1=т и к2=п-т-1

(п - число наблюдений, т - число параметро при переменных регрессии). Фактическое значение F-критерия больше табличного ^факт>Ркрит. Таким образом, признается статистическая значимость зависимости в целом [5, 6].Все значения показателей качества связи больше критических, что свидетельствует о статистической значимости полученных зависимостей.

Качество зависимости нормативно-расчетной трудоемкости (ТНОрМ _расч.) и трудоемкости рассчитанной по модели (ТзпКМ), оценивалось с помощью средней ошибки абсолютной аппроксимации (А):

А =1Е

Пу

1 п Т - Т

1 ^^ | з.пкм норм-расч.

Т

х 100%.

Качество зависимости экспертно-техноло-гической трудоемкости (ТЭКст-технол.) и трудоемкости, рассчитанной по экономико-математической модели (ТзпКМ), оценивалось с помощью среднего абсолютного отклонения (А):

1 п Т

д =1Е ь

п Е

- Т

■ эксп.-технол.

Т

х 100%.

Объективность, погрешность расчетов, получаемых методами укрупненного определения трудоемкости, должны быть в допустимых пределах (см. табл. 6) [7]. Погрешности А и А не превышают 20%, таким образом, можно сделать вывод о достоверности применяемого метода расчета.

Таблица 4

Значения трудоемкости изготовления детали из ПКМ соответствующего класса

Название класса

Трудоемкость изготовления изделия, Тбазаь н.-ч.

конструкции без усиливающих и сотовых элементов (1=1)

^ + 34,63.5ДеТ -20,853.5^ + ^

5

V

дет

е

• Б

осн.слоев

У

{

конструкции с усиливающими и сотовыми элементами ^=2)

-у _ -у о выкл

база2 база2 осн.сслоев

0,743

^ТТРТ

+0,083 • б д2еТ + 5,061+_6328_

дет е " Ьп(Бдет)

• б в

осн.слоев

__£ Т выкл

база1 база1 осн.слоев

У

Таблица 5

Показатели качества выбора формы связи

Показатели трудоемкости Коэффициент детерминации, R2 Критерий Фишера, F Средняя ошибка аппроскимации4, % Среднее абсолютное отклонение Д, %

Расчетное Критическое Расчетное Критическое Расчетное Критическое Расчетное Критическое

ДСЕ без усиливающих и сотовых элементов

£Тнорм.-расч. 0,97 0,8 18,6 7,71 15,0 20,0 - -

£Тэксп. -технолог. 0,99 534 5,59 - - 17,2 20,0

ДСЕ с усиливающими и сотовыми элементами

£Тнорм.-расч. 0,96 0,8 87,5 3,18 17,5 20,0 - -

£Тэксп.-технолог. 0,97 479 3,01 - - 19,5 20,0

Таблица 6

Допускаемая погрешность укрупненного расчета трудоемкости в зависимости от типа производства

Предполагаемое Допускаемая

производство погрешность, (±)%

Массовое 10,0

Крупносерийное 15,0

Среднесерийное 20,0

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

ЗАКЛЮЧЕНИЕ

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

СПИСОК ЛИТЕРАТУРЫ

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

1. ГОСТ 14.201. Обеспечение технологичности конструкции изделий. Общие требования.URL: http://www.opengost.ru (дата обращения 22.04.2014). [[Providing technological effectiveness of design products. General requirements, (in Russian), Federal standard ]].

2. ГОСТ Р 50-54-86-88. Автоматизированная система технологической подготовки производства. Состав и порядок разработки. М.,1988.URL:http://www.opengost.ru(дата обращения 22.04.2014).[[ The automated system of technological preparation of production. Composition and development, Federal standard, (in Russian). Moscow, 1988. ]]

3. Селиванов С. Г., Гузаиров М. Б. Системотехника инновационной подготовки производства в машиностроении. М.: Машиностроение, 2012. 132 с. [[ S. G. Selivanov, M. B. Guzairov, Systems engineering innovative pre-production in mechanical engineering, (in Russian). Moscow: Mashinostroeniye, 2012. ]]

4. Кузнецова К. С., Селиванов С. Г. Экспертно-логический метод оптимизации предварительных проектов технологической документации // Вестник УГА-ТУ.2012.Т. 16, №1(46). С.81-87. [[ K. S.Kuznetsova, S. G. Selivanov, ''Expert-logical method optimization preliminary drafts of technical documentation,', (in Russian), in Vestnik UGATU, vol. 16, no. 1 (46), pp. 81-87, 2012.]]

5. Елисеева И. И., Курышева С. В., Костеева Т. В., Бабаева И. В., Михайлов Б. А. Эконометрика. М.: Проспект, 2010. [[ I. I.Yeliseyeva, S. V.Kurysheva, T. V. Kosteyeva, I. V. Babayeva, B. A. Mikhaylov, Econometrics, (in Russian). Moscow: Prospekt, 2010.]]

6. Елисеева И. И., Курышева С. В., Гордиенко Н. М., Костеева Т. В., Бабаева И. В., Михайлов Б. А. Практикум по эконометрике: учеб. пособие. М.: Финансы и статистика, 2003. 187 с. [[ I. I. Yeliseyeva, S. V. Kurysheva, N. M. Gordi-yenko, T. V. Kosteyeva, I. V. Babayeva, B. A. Mikhaylov, Workshop on Econometrics, (in Russian). Moscow: Finansy i statistika, 2003.]]

7. Амиров Ю. Д., Алферова Т. К. Технологичность конструкции изделия. М.: Машиностроение, 1990. 132 с. [[ Y. U. D. Amirov , T. K. Alferova, Technological effectiveness of product design, (in Russian). Мoscow: Mashinostroeniye, 1990.]]

ОБ АВТОРАХ

ЮРЬЕВ Виктор Леонидович, проф. каф. Технологии машиностроения. Дипл. инж. (УГАТУ, 1973). Д-р техн. наук (УГА-ТУ, 2000). Иссл. в обл. комплексной технологии сложнофа-сонных деталей ГТД.

ИРКАБАЕВ Ирик Мигранович, вед. инж.-технол. Дипл. инж. (УГАТУ, 1976).

СУХОВА Надежда Александровна, доц. каф. экономики предпринимательства. Дипл. экон.-менедж. (УГАТУ, 2009). Канд. техн. наук. (УГАТУ, 2005). Иссл. в обл. технологии ионно-плазменного насаждения и термической обработки многослойных покрытий.

КАДЫРОВА Жанна Минифануровна, асп. каф.технологии машиностроения. Дипл. экон.-мат. (БашГУ, 2009). Готовит дис. о моделях трудоемкости изготовления авиационных гражданских ГТД в условиях кооперации.

METADATA

Title: Forecasting model the complexity of manufacturing nacelles civil Aviation gas turbine engine in CAM.

Authors: 1 V. L. Yuriev, 2 I. M. Irkabaev, 3 N. A. Sukhova, 4 Zh. M. Kadyrova.

Affiliation:

1-4 Institute of Technology and the organization of production (NIIT), Russia.

1,3,4 Ufa State Aviation Technical University (UGATU), Russia.

Email: 4prep08@mail.ru.

Language: Russian.

Source: Vestnik UGATU (scientific journal of Ufa State Aviation Technical University), vol. 19, no. 1 (67), pp. 37-44,2015. ISSN 2225-2789 (Online), ISSN 1992-6502 (Print).

Abstract: As part of the development of a methodology of determining the complexity of producing Aviation gas turbine engine in cooperation mathematical model determining the complexity of manufacturing parts and assembly units of polymeric composite materials Aviation gas turbine engine nacelle.

Key words: competitiveness; technological forecasting; CAM; mathematical modeling; labor input; Polymer composite materials.

About authors:

YURIEV, Viktor Leonidovich, Prof., Dept. of Mechanical Engineering. Dr. of Tech. Sci. (UGATU, 2000). IRKABAEV, Irik Migranovich, leading engineer-tehnol. (UGATU, 1976).

SUKHOVA, Nadezhda Aleksandrovna, assis Prof., Dept. of Business Economics, candidate of Tech. Sci. (UGATU, 2005).

KADYROVA, Zhanna Minifanurovna, Postgrad. (PhD) Student, Dept. of Mechanical Engineering. Master of Mathematical Methods & Economics (BashGU, 2009).

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