WWW.DISS.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА
(Авторефераты, диссертации, методички, учебные программы, монографии)

 

Pages:   || 2 |

«ЧИСЛЕННО-АНАЛИТИЧЕСКИЕ МЕТОДЫ СТАНДАРТНЫХ ЭЛЕМЕНТОВ ДЛЯ МОДЕЛИРОВАНИЯ СТАЦИОНАРНЫХ ФИЗИЧЕСКИХ ПОЛЕЙ В ЛИНЕЙНЫХ КУСОЧНООДНОРОДНЫХ И НЕЛИНЕЙНЫХ СРЕДАХ ...»

-- [ Страница 1 ] --

На правах рукописи

Пашковский Александр Владимирович

ЧИСЛЕННО-АНАЛИТИЧЕСКИЕ МЕТОДЫ СТАНДАРТНЫХ

ЭЛЕМЕНТОВ ДЛЯ МОДЕЛИРОВАНИЯ СТАЦИОНАРНЫХ

ФИЗИЧЕСКИХ ПОЛЕЙ В ЛИНЕЙНЫХ КУСОЧНООДНОРОДНЫХ И НЕЛИНЕЙНЫХ СРЕДАХ

05.13.18 – Математическое моделирование,

численные методы и комплексы программ

Автореферат диссертации на соискание ученой степени доктора технических наук

Новочеркасск – 2014 2

Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «ЮжноРоссийский государственный политехнический университет (НПИ) имени М.И.

Платова».

Научный консультант: доктор технических наук, профессор Ткачев Александр Николаевич

Официальные оппоненты: Князев Сергей Юрьевич, доктор технических наук, профессор, ФГБОУ ВПО «Донской государственный технический университет», профессор кафедры «Математика»

Наседкин Андрей Викторович, доктор физико-математических наук, доцент, ФГАОУ ВПО «Южный федеральный университет», профессор кафедры «Математическое моделирование»

Сипливый Борис Николаевич, доктор технических наук, профессор, ФГАОУ ВПО «Волгоградский государственный университет», профессор кафедры «Теоретическая физика и волновые процессы»

Ведущая организация: Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Воронежский государственный технический университет»

Защита состоится « 4 » июля 2014 г. в 10-00 на заседании диссертационного совета Д 212.304.02, созданного на базе Южно-Российского государственного политехнического университета (НПИ) имени М.И. Платова, в 149 ауд. главного корпуса университета по адресу: 346428, Ростовская область, г. Новочеркасск, ул.

Просвещения, 132.

С диссертацией можно ознакомиться в научно-технической библиотеке ФГБОУ ВПО ЮРГПУ (НПИ) имени М.И. Платова и на сайте www.npi-tu.ru Автореферат разослан “” _ 2014 г.

Ученый секретарь диссертационного совета Д 212.304. к.т.н., профессор А.Н. Иванченко

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы исследований.




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

Это приводит к увеличению размерности решаемых систем линейных алгебраических уравнений (СЛАУ) и, как следствие, к значительному росту вычислительных затрат, времени расчетов, особенно в тех случаях, когда эти расчеты являются многовариантными.

Подобные особенности делают проблематичным использование вышеперечисленных методов для решения следующих часто возникающих на практике задач:

– комплексных, когда необходимо выполнить одновременный расчет нескольких физических полей;

– непрямого контроля параметров полей и параметрического управления;

– расчета физических полей, выполняемого в условиях многовариантности, например, связанной с изменением геометрии (характерных размеров устройств и взаимного положения их элементов);

– оптимизации, связанной с выбором оптимальных параметров устройств;

– обратных задач, например, при:

a) непрямом измерении характеристик материалов;

б) большом числе угловых точек, в которых решение имеет особенности;

в) возникновении осцилляции решения;

г) наличии мили- и микровключений, вносящих большие погрешности и приводящих к неустойчивости численного решения;

д) расчете трехмерных физических полей в нелинейных кусочно-однородных средах.

Это обосновывает актуальность разработки новых эффективных методов, обеспечивающих снижение необходимых вычислительных ресурсов, сокращение времени расчетов и их устойчивость, повышение точности расчетов физических полей в кусочно-однородных линейных и нелинейных средах, в том числе, когда решение в угловых точках области имеет особенность, при наличии тонких включений, возникновении осцилляции решения. В диссертационной работе построены новые численно-аналитические методы стандартных элементов для расчета физических полей, в основе которых лежит представление расчетной области в виде объединения стандартных элементов (СЭ), а также нахождение приближенного решения с использованием аналитических решений краевых задач в СЭ, принадлежащих пространству Соболева W22 и задание операторов «связи» следов решения с нормальными производными на границах СЭ. Разработанные методы являются эффективными, обеспечивают снижение необходимых вычислительных ресурсов, сокращение времени расчетов, устойчивость и высокую точность. Это достигается с помощью:





– снижения объема необходимых вычислений, которое осуществляется в результате увеличения размеров элементов разбиения расчетной области при введении СЭ и блоков из них без потери точности;

– повышения гладкости приближенного решения, достигаемого при одновременном уменьшении числа СЭ в расчетной области, например, по сравнению с числом конечных элементов в МКЭ;

– использования точных (аналитических) решений краевых задач в СЭ;

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

– повышения точности нахождения нормальных производных и точности выполнения условий их склейки на границах СЭ;

– использования коэффициентов рядов Фурье для следов решения на границах СЭ при «склейке» граничных условий в отличие от общепринятой «склейки» с использованием узловых значений;

– снижения размерности матриц решаемых СЛАУ и их блочно-ленточной структуры;

– учета особенности решения в окрестности угловых точек;

– возможности использования при реализации математических моделей параллельного программирования и многоядерных ПК.

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

Для достижения цели были поставлены и решены следующие задачи:

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

2. Разработаны теоретические основы новых численно-аналитических комбинированных методов: метода СЭ с использованием рядов Фурье и конечных элементов (КМСФиКЭ), метода вспомогательных функций и конечных элементов (КМВФиКЭ) для расчета поля в нелинейных средах, а также методики их численной реализации.

3. Разработаны методики «склейки» двух и трехмерных СЭ с использованием граничных узловых значений решения и коэффициентов Фурье, которые реализованы в комбинированном МСЭФ с МКЭ.

4. Разработана новая методика «склейки» СЭ при использовании вспомогательных функций, обеспечивающая возможность построения блоков СЭ, комбинирование МВФ и МКЭ, и позволяющая повысить точность расчетов.

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

– дифференциальным оператором, соответствующим решаемой в СЭ задаче;

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

6. В результате проведенных численных экспериментов выполнены сравнительные оценки точности МКЭ, МГЭ, МТИ, МСЭФ, МВФ, КМГиКЭ, КМСФиКЭ, КМВФиКЭ.

7. Разработан алгоритм линеаризации нелинейных характеристик среды, заполняющей СЭ при использовании МСЭФ и МВФ для моделирования физических полей в нелинейных средах.

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

9. Получены и проанализированы результаты решений двух и трехмерных модельных краевых задач при условии варьирования свойств расчетных областей (однородная, кусочно-однородная, с источниками и без них).

10. Определена и проанализирована структура блочно-ленточной СЛАУ, возникающей при решении исходных краевых задач разработанными методами.

11. Проведены численные эксперименты с целью сравнения разработанного МВФ и МКЭ при различных размерах конечных элементов и варьировании свойств расчетных областей.

12. Оценена точность МВФ и КМВФиКЭ при наличии осцилляции решения в расчетной области, а также значительной диспропорции размеров СЭ в различных направлениях.

13. Выполнена теоретическая оценка точности представлений следов решения краевой задачи на границах СЭ при использовании МВФ и МСЭФ.

14. Выполнена оценка точности МВФ, КМСФиКЭ, КМВФиКЭ при использовании экспериментальных данных, результатов расчетов, выполненных МКЭ, МГЭ, МТИ, КМГиКЭ с помощью пакетов прикладных программ применительно к расчетам:

а) магнитного поля линейного двигателя с постоянными магнитами;

б) температурного поля якоря тягового электродвигателя;

в) электростатического поля каталитической системы;

г) магнитного поля трех- и двухмерного подъемного модуля и создаваемой им подъемной силы.

Методы исследования. Для решения поставленных в диссертационной работе задач использованы численные методы решения линейных и нелинейных дифференциальных уравнений в частных производных, численные методы решения СЛАУ, методы математической физики, теории вариационного исчисления, интегрального и дифференциального исчислений, методы математического моделирования, методы вычислительной математики, специализированные программные среды, выполнены численные и физические эксперименты.

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

Научная новизна результатов исследования состоит в следующем:

– в МСЭФ и КМСФиКЭ предложена вариационная постановка решаемых краевых задач, отличающаяся использованием новых “основных” функционалов, найденных в результате преобразования интеграла Дирихле, минимизация которых в области, полученной в результате объединения СЭ Gi, осуществляется в классе функций, принадлежащих пространству Соболева W22 Gi. Это обеспечивает существование нормальных производных решения на границах СЭ, а также существование разложений следов решений на границах СЭ в ряды Фурье;

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

а) с использованием узловых значений решения;

б) с использованием коэффициентов Фурье следов решения на границах СЭ.

– предложены новые подходы к решению краевых задач (МВФ и КМВФиКЭ), отличающиeся тем, что они не предполагают использование вариационного принципа для решения краевых задач, а построение блочно-ленточных СЛАУ при их реализации осуществляют с помощью вспомогательных функций, которые определяются с учетом решаемых в СЭ краевых задач;

– разработан новый алгоритм построения блочно-ленточных СЛАУ, к решению которой сводится решение краевой задачи, отличающийся от МКЭ тем, что он не требует минимизации функционала по узловым значениям решения, а в качестве неизвестных в СЛАУ служат коэффициенты Фурье следов решений на границах СЭ или узловые значения решений и коэффициенты Фурье одновременно;

– предложена методика оценки погрешности математической модели, построенной с использованием СЭ, отличающаяся тем, что позволяет обосновать выбор СЭ для дискретизации расчетной области с учетом требуемой точности решения и необходимых вычислительных ресурсов;

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

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

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

– снизить степень дискретизации (размеры элементов разбиения) расчетных областей, уменьшить размерность решаемой СЛАУ, повысить скорость выполнения многовариантных расчетов поля за счет введения СЭ;

– уменьшить погрешность выполненных расчетов за счет более точного (с использованием аналитических решений) представления приближенного решения в СЭ;

– получать непрерывно дифференцируемые решения внутри СЭ, что обеспечивает снижение погрешности вычисления векторов физических полей по их потенциалам и высокую точность нахождения их интегральных характеристик;

– обеспечить высокую точность расчетов поля при наличии особенностей решения, его осцилляции и наличии тонких включений в расчетной области;

– учитывать асимптотику решения и наличие особенности в окрестности угловых точек, решать внешние краевые задачи с использованием аналитических решений в рассмотренных неограниченных СЭ;

– использовать разработанные методы для решения комплексных задач неразрушающего и непрямого контроля, параметрического управления, задач оптимизации и других многовариантных задач;

– решать задачи расчета поля в расчетных областях, заполненных неоднородными средами, в том числе с нелинейными характеристиками, при диспропорции размеров СЭ, при наличии мили- и микровключений;

– служить в качестве быстродействующих эффективных решателей в комплексах программ расчёта физических полей.

Основные положения, выносимые на защиту:

1. Численно-аналитические методы МСЭФ, МВФ и алгоритмы их реализации, позволяющие с повышенной точностью решать полевые задачи в двухмерных и трехмерных расчетных областях, в том числе при наличии нелинейных включений.

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

3. Вариационные постановки краевых задач в МСЭФ и КМСФиКЭ и новые «основные» функционалы, позволяющие повысить гладкость решения, точность нахождения нормальных производных потенциалов, обеспечивающие существование разложений следов решений на границах СЭ в ряды Фурье.

4. Новые методики построения уравнения «связи» между следами нормальной производной и решения на границах СЭ, согласно которым вместо функционала в виде интеграла Дирихле используются «основные» функционалы, содержащие интегралы по границам СЭ.

5. Численные реализации МСЭФ, КМСФиКЭ на основе двух разработанных методик «склейки» СЭ, обеспечивающие возможность комбинирования МСЭФ с МКЭ с использованием:

– узловых значений решения на границах СЭ;

– коэффициентов Фурье следов решения на границах СЭ.

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

7. Новые методики построения уравнений «связи» между следами нормальной производной и решения на границах СЭ на базе МВФ и КМВФиКЭ, позволяющие формировать блоки СЭ, которым соответствуют блочные СЛАУ, что обеспечивает повышение точности расчетов.

8. Новые способы формирования блочно-ленточных СЛАУ согласно КМСФиКЭ и КМВФиКЭ, с использованием коэффициентов Фурье следов решения на границах СЭ и узловых значений решения в нестандартных подобластях.

9. Методики построения «основных» функционалов для краевых задач в случаях двух и трехмерных типовых СЭ, отличающихся от интеграла Дирихле, которые образуют библиотеку «основных» функционалов.

10.Алгоритмы приближенного учета нелинейных материальных характеристик среды в СЭ на каждом шаге итерации, построенные на базе МСЭФ и МВФ, позволяющие выполнять расчет двух и трехмерных физических полей путем сведения нелинейной среды к кусочно-однородной линейной.

Апробация работы. Основные результаты работы докладывались на международных конференциях:

– Международная научно-техническая конференция «Молодая наука новому тысячелетию». – Набережные Челны: КамПИ, 1996;

– Международная научно-техническая конференция «Компьютерные технологии в науке, производстве, социальных и экономических процессах». – Новочеркасск: ЮРГТУ (НПИ), 2000;

– Международная научно-практическая конференция «Теория, методы и средства измерений, контроля и диагностики». – Новочеркасск: ЮРГТУ(НПИ), 2000;

– XIV Международная научная конференция “Математические методы в технике и технологиях». –Смоленск, 2001;

– III Международная научно-техническая конференция «Компьютерное моделирование 2002». – Санкт-Петербург, 2002;

– XVII Международная научная конференция «Математические методы в технике и технологиях». – Кострома: КГТУ, 2004;

– XXVI Международная научная конференция «Математические методы в технике и технологиях». – Саратов: СГТУ, 2013;

Российских конференциях:

– I Российская национальная конференция по теплообмену. – Москва, 1994;

– II Российская национальная конференция по теплообмену. – Москва, 1998;

– IV Всероссийский симпозиум «Математическое моделирование и компьютерные технологии». – Кисловодск, 2000;

– 13 Школа-семинар молодых ученых и специалистов «Физические основы экспериментального и математического моделирования процессов газодинамики и тепломассобмена в энергетических установках». – Санкт-Петербург, 2001;

– II Общероссийская научно-техническая конференция «Новые технологии в азотной промышленности». – Невинномысск: НТИ, 2007;

– Общероссийская научно-практическая конференция «Математическое моделирование, компьютерные и информационные технологии в технике, экономике и образовании». – Невинномысск: СевКавГТУ, 2009;

– XII Всероссийский симпозиум «Прикладная и промышленная математика». – Сочи-Адлер: СГУ, 2011.

Публикации. Основное содержание диссертационной работы отражено в научных публикациях, среди которых 17 статей в научных журналах, рекомендованных ВАК, 1 авторское свидетельство на изобретение, 2 свидетельства о государственной регистрации программы для ЭВМ.

Реализация и внедрение результатов работы. Тема диссертационной работы соответствует «Основным направлениям НИР ЮРГПУ (НПИ)» теме «Интеллектуальные электромагнитные устройства, системы и комплексы», разделу «Математическое моделирование физических полей и процессов в технических устройствах и системах». Результаты диссертационной работы использовались при выполнении НИР в соответствии с тематическим планом НИР Невинномысского технологического института (филиал) ФГАОУ ВПО «СКФУ» «Математическое моделирование и информационные технологии в решении задач машиностроения»

(ГР 01200614265), а также НИР в рамках программы «Участник молодежного научно-инновационного конкурса» («У.М.Н.И.К») (государственный контракт №:6642р/9109 от 02.03.09). Разработанные в работе численные методы, а также алгоритмы, методики, пакеты программ внедрены на следующих предприятиях:

ОАО «Всероссийский научно-исследовательский и проектно-конструкторский институт электровозостроения» (ОАО «ВЭлНИИ», г. Новочеркасск), ООО «ХимПроект» (г. Невинномысск), ООО «Экоцентр» (г. Невинномысск), в учебном процессе НТИ (филиал) ФГАОУ ВПО «СКФУ» (г. Невинномысск). Программы для ЭВМ зарегистрированы в Объединенном фонде электронных ресурсов «Наука и образование» (№ 15811 от 03 июня 2010 г.) и реестре программ Федеральной службы по интеллектуальной собственности, патентам и товарным знакам (№ 2010615042 от 05 августа 2010 г.). Внедрение результатов работы подтверждено соответствующими актами.

Структура и объем диссертации. Диссертационная работа состоит из введения, пяти глав, заключения, списка используемой литературы из 156 наименований и приложений. Работа изложена на 364 страницах, включая 38 страниц приложений и 137 иллюстраций.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе «Численные методы моделирования физических полей»

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

Например, для МКР и МКЭ, в развитие которых большой вклад внесли, в том числе работы ученых Ю.А. Бахвалова, Г.А. Гринберга, К.С. Демирчяна, Ю.А.

Иосселя, О.Ф. Ковалева, Г.И. Марчука, А.В. Наседкина, И.И. Пеккера, А.А. Самарского, А.Н. Тихонова и др., а также работы зарубежных ученых И. Бабушки, Р.

Галлагера, Ж. Деклу, О. Зенкевич, Э. Митчелла, Дж. Одена, Ж. Сабоннадьер, Л.

Сегерлинда, П. Сильвестера, Г. Стренга, Ф. Сьярле, Дж. Фикса и др., характерны:

– линейная или полиномиальная (невысокого порядка) аппроксимация решения;

– невысокая точность выполнения условий склейки нормальных производных решения на границе сред с различными свойствами материалов;

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

– рост погрешности полевых расчетов при наличии осцилляции решения и его особенности в окрестности угловых точек;

– необходимость использования специальных методов формирования, хранения и обработки матриц решаемых СЛАУ;

– значительные затраты вычислительных ресурсов, возникающие при решении трехмерных задач;

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

МГЭ и комбинированные методы расчета, на нем основанные, получившие свое развитие, в том числе в работах В.И. Астахова, В.Г. Баженова, Ю.А. Бахвалова, П. Бенерджи, К. Бребия, Э.В. Колесникова, С. Крауча, А.Г. Никитенко, А.Н.

Ткачева, О.В. Тозони и др., хотя и позволяют сократить степень дискретизации в расчетной области, по сравнению с МКЭ, однако:

– требуют решения интегральных уравнений с интегралами, содержащими функции с особенностями;

– приводят к формированию полностью заполненных матриц СЛАУ, часто имеющих большие числа обусловленности;

– требуют разработки и использования специальных приемов формирования матриц, особых методов решения СЛАУ для обеспечения приемлемой скорости сходимости итерационных процессов;

– методики использования некоторых комбинированных методов оказываются неприменимыми для решения трехмерных задач.

Проведенный анализ показал, что при расчете физических полей в многосвязных кусочно-однородных средах с линейными и нелинейными характеристиками, возможности перечисленных методов ограниченны. Их использование, особенно в условиях наличия особенностей решения, его осцилляции, наличии включений, многофазности требует увеличения необходимого для подготовки исходных данных времени, высокой степени дискретизации расчетных областей, увеличения размерности решаемых СЛАУ, роста вычислительных затрат и времени расчета.

Альтернативными методами расчета физических полей в этих условиях становятся предложенные в работе численно-аналитические методы стандартных элементов, которые, как следует из анализа проведенных исследований, могут позволить:

– сократить степень дискретизации областей благодаря введению СЭ больших размеров;

– получать непрерывно дифференцируемые решения внутри каждого СЭ;

– повысить точность расчетов и сократить вычислительные затраты при наличии в расчетной области особенностей решения, осцилляции, узких включений.

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

где u uM ; D DM, u ; M ; u uM – подлежащая определению функция;

– область (двухмерная или трехмерная), в которой ищется распределение поля;

M – точка области ; D DM, u – функция, задающая материальные свойства среды в области ; f f M – заданная функция, определяемая источниками поля.

На границе области должны быть заданы граничные условия, в качестве которых могут выступать условия первого рода (Дирихле), второго рода (Неймана) или третьего рода (смешанные).

При расчете физических полей в линейных кусочно-однородных средах расчетная область представима в виде k. При этом функция D DM в каждой области k задается равенством: DM Dk const ; M k, а на границе раздела сред должны выполняться условия непрерывности (условия связи), характерные для каждого конкретного вида физического поля.

Для решения краевых задач в постановке (1) разработанными методами в главе 1 введено понятие СЭ, которое предлагает рассматривать СЭ как математический объект, определенный геометрией, граничными условиями на частях его границ, дифференциальным оператором, определенным в СЭ, а также оператором, определяющим «связи» между решениями в смежных (соседних) СЭ. С учетом введенного определения, а также анализа типовых постановок краевых задач, решаемых при расчете физических полей, сформирована библиотека из 19 различных СЭ, представленных в таблице 1.

Таблица 1 Библиотека СЭ Постоянный магнит Сектор двойного двумя границами Параллелепипед Для обозначения СЭ в соответствии с его геометрической формой и граничными условиями на частях границы использовались следующие записи. Например, для прямоугольного СЭ с условиями Дирихле на сторонах, указывается: RECD (REC-сокращение английского перевода слова «прямоугольник», а D1234 – задание условия Дирихле на всех сторонах СЭ). Аналогично, для СЭ в виде параллелепипеда с условиями Неймана на 1, 2, 3 гранях и условиями Дирихле на 4, 5, гранях указывается: PARN123D456.

В главе сформулированы основные задачи исследований, определены пути их решения.

Во второй главе «Численно-аналитические методы стандартных элементов для моделирования физических полей в линейных средах» для моделирования физических полей в линейных кусочно-однородных средах разработаны два базовых метода: МСЭФ и МВФ. Сформулированы основные виды краевых задач, к решению которых сводится расчет физических полей (магнитных, электрических, тепловых и т.д.). Изложены теоретические основы МСЭФ и МВФ, суть которых состоит в следующем.

Согласно МСЭФ предлагается следующая вариационная постановка краевых задач. Предполагается, что решение, определенное в СЭ Gi принадлежит пространству Соболева W22(Gi). Используются два новых, так называемых, «основных» функционала. Первый из них имеет вид:

где функции Fi ui определены через узловые значения следов решения ui на границах СЭ Si. Второй функционал представим в виде:

где ai, a j – наборы коэффициентов Фурье следов решения на границах СЭ Si, а сij – найденные арифметические выражения, зависящие от размеров СЭ и его положения в пространстве.

Согласно МСЭФ операторы «связи» выражаются через следы решения на границах СЭ в виде одномерного ряда Фурье. Для нахождения решения в СЭ предложено две методики: с помощью узловых значений решения на границах СЭ и с помощью коэффициентов рядов Фурье следов решения на границах СЭ. Данные методики соответствуют двум реализованным в МСЭФ способам «склейки»

СЭ и решений в них. При этом используются уравнения «связи» между решением ui и его нормальной производной на границе СЭ, а также полученные «основные»

функционалы из библиотеки СЭ.

Для каждого СЭ заданной геометрии из построенной библиотеки получены «основные» функционалы МСЭФ вида (2), (3), соответствующие типовым краевым задачам в СЭ при возможных условиях Дирихле и Неймана на границах. Найденная в соответствии с МСЭФ совокупность функционалов образует библиотеку «основных» функционалов, которая применяется при выполнении разбиения расчетных областей. При этом формирование функционала, соответствующего краевой задаче, сводится к простому суммированию «основных» функционалов.

Показано, что решения ui в СЭ Gi и их нормальные производные ui n на границах СЭ, входящие в каждый интеграл функционала (2) и (3), выражаются либо через неизвестные значения решения в узлах, расположенных на границах СЭ и лежащих внутри расчетной области, либо через неизвестные коэффициенты рядов Фурье следов решений на границах СЭ. Формирование блочно-ленточной СЛАУ, соответствующей решаемой краевой задаче, выполняется в результате минимизации функционала общей задачи по неизвестным узловым значениям или коэффициентам Фурье. Решение исходной краевой задачи во всех точках СЭ определяется по найденным в результате решения СЛАУ неизвестным. Для нахождения решения в отдельных точках СЭ или линиях внутри него, вводится дополнительный СЭ, границы которого содержат эти точки или совпадают с линиями. Это позволяет, не находя вид решения внутри СЭ, вычислять искомые значения с использованием следов решения на рассматриваемых границах дополнительных СЭ.

Разработанный МВФ также требует введения СЭ в расчетной области. Однако в нем реализован иной подход для обеспечения «склейки» СЭ и формирования СЛАУ, соответствующей решаемой задаче. Согласно МВФ решение в каждом СЭ Gi ищется в классе функций из пространства W22 (Gi ). При этом в каждом СЭ вводятся специальным образом определенные в нем вспомоRECD1234 гательные функции. При решении в СЭ задачи Дирихле, Если на частях границы СЭ должно выполняться условие Неймана, то на данных частях границы в нуль должна обращаться нормальная производная вспомогательных функций. Такая методика позволяет получить интегродифференциальное выражение, связывающее решение и нормальную производную на границе СЭ. Это позволяет:

– исключить в уравнениях «склейки» неизвестные значения нормальных производных решения и использовать для выполнения условий «склейки» лишь неизвестные значения решения на границах СЭ;

– получить на каждой стороне «склейки» СЭ предельные значения интегродифференциальных выражений, приравнивание которых позволяет сформировать уравнения блочно-ленточной СЛАУ, соответствующей решаемой краевой задаче.

Реализации МВФ не требует использования вариационного подхода, что отличает его от МСЭ, МСЭФ и КМСФиКЭ. Так для пары смежных четырехугольных СЭ G1, G2, решение в которых задается функциями u1, u2 при введении вспомогательных функций V12, V21 интегро-дифференциальное уравнение «связи» для следов решений на общей границе «склеенных» СЭ имеет вид:

где ij – следы функций u1, u2 на границах СЭ. Уравнения (4) записываются на каждой стороне «склейки» для всех пар СЭ. Так, при «склейке», например, бесконечного СЭ с соседними СЭ в качестве неизвестных в уравнениях (4) являются коэффициенты Фурье следов решения на сторонах конечных и бесконечных СЭ. В результате формируется блочно-ленточная СЛАУ, соответствующая исходной краевой задаче. После решения СЛАУ находятся коэффициенты разложения искомой функции u в ряд Фурье на всех сторонах СЭ. Далее, решение краевой задачи в области СЭ определяется с использованием найденных на его границах коэффициентов Фурье. При необходимости нахождения решения в отдельных точках или на линиях, лежащих внутри СЭ, как и в МСЭФ, вводятся дополнительные СЭ, границы которых содержат эти точки или совпадают с линиями. Далее, искомое решения определяется через следы решения на границах дополнительных СЭ. Уравнения «связи» типа (4) конкретизированы для всех СЭ сформированной библиотеки.

В соответствии с алгоритмом МВФ, построенным для случаев двух-, трехмерных СЭ и бесконечного СЭ получены уравнения «связи» для краевой задачи Дирихле в области, образованной объединением СЭ. Сформированы соответствующие СЛАУ, проанализирована их структура. Анализ показал, что при выполнении требований, сформулированных в алгоритме МВФ к разбиению области, СЛАУ является разреженной и имеет блочно-ленточную структуру.

Показано, что при распределении в области СЭ источников поля при решении краевой задачи для уравнения Пуассона могут быть использованы уравнения «связи», полученные ранее МСЭФ и МВФ применительно к решению краевых задач для уравнения Лапласа. При этом уравнения «связи» сохраняют свою структуру, однако, содержат дополнительные слагаемые, учитывающие наличие источников поля.

Разработаны алгоритмы использования МСЭФ и МВФ для решения задач расчета магнитного поля, создаваемого постоянными магнитами, наличие которых приводит к необходимости решения краевой задачи для уравнения Пуассона. Показано, что введение СЭ, покрывающих область, заполненную постоянными магнитами, позволяет учесть намагниченность М условиями «склейки» смежных с постоянным магнитом СЭ. Показано, что в этом случае СЛАУ может быть сформирована с использованием уравнений «связи», полученных ранее для решения краевых задач для уравнения Лапласа. При этом коррекция каждого уравнения «связи» СЛАУ для соответствующей части границы двух- или трехмерного СЭ, покрывающего постоянный магнит, сводится к изменению правой части уравнения. Дополнительное слагаемое для уравнения Пуассона выражается через намагниченность М.

Рассмотрен также случай, когда источниками поля в МС выступают катушки с током. Моделирование магнитного поля сводится к формированию и решению СЛАУ относительно проекций векторного потенциала. При формировании СЛАУ учитываются условия на границах «смежных» СЭ, которые могут быть выражены через производные проекций векторного потенциала. Дополнительное слагаемое может быть получено с использованием частных решений уравнения Пуассона, либо интегрированием выражения, содержащего заданную в СЭ функцию плотности тока, распределенной по его области.

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

В третьей главе «Численно-аналитические методы стандартных элементов для моделирования физических полей в нелинейных средах» разработанные ранее для кусочно-однородных линейных сред алгоритмы МСЭФ и МВФ обобщены применительно к решению задач расчета физических полей в нелинейных средах.

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

Разработана методика расчета на каждом итерационном шаге среднего значения индукции B в СЭ, которое используется для линеаризации материального уравнения. Найденные в результате решения блочно-ленточной СЛАУ коэффициенты рядов Фурье следов потенциала А на границах СЭ используются для определения проекций индукции на границах СЭ, а также на расположенных внутри СЭ вспомогательных линиях для последующего вычисления среднего значения индукции B в СЭ.

Например, для прямоугольного СЭ (рис.2) проекции индукции находятся на границах s1, s2, s3, s4 и на вспомогательных линиях l1, l2. Среднее значение B в используются значения коэффициентов Фурье, найденРис.2 К линеаризации материального оператора ные для следов решения на граничных si и вспомогательных отрезках li, а искомое значение B определяется усреднением всех найденных значений B r при r 1,9.

Такой способ определения индукции на каждом шаге не требует нахождения потенциала поля в виде отрезка ряда Фурье на границе или в области СЭ. После нахождения среднего значения и линеаризации материального уравнения выполняется расчет. По его результатам аналогичным образом находятся средние значения индукции Bk в частях СЭ (рис.2), ограниченных сторонами si СЭ и отрезками прямых lj, которые сравниваются со значением B r. Если погрешность не превосходит допустимую, осуществляется переход к вычислению нового значения магнитной проницаемости, а затем к следующему шагу итерационного процесса. В противном случае выполняется дополнительное разбиение СЭ. Процедура выполнения дополнительного разбиения автоматизирована и не требует корректировки исходных данных. Соответствующее значение магнитной проницаемости в СЭ, расположенном в части расчетной области с нелинейными магнитными характеристиками, на каждом r-м шаге (r = 0, 1, 2…) определяется по кривой намагничивания.

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

Первый из них реализован в виде двух модификаций, применение которых проиллюстрировано применительно к решению краевых задач расчета магнитного поля в области, содержащей включения с нелинейными характеристиками. В частности, реализация КМСФиКЭ подробно расG2 G3 G смотрена на примере расчета физического поля, потенциал которого в расчетной двухмерной обG новке:

где D кусочно-постоянная функция в одних часиспользованием узловых значений тях области (G1, G2,..), нелинейная в других (1, 2) и uW21(G). Показано, что с использованием КМСФиКЭ задача (5) сводится к минимизации функционала:

который представим в виде суммы интегралов по границам СЭ, размещенным в частях Gi, и по областям k. Методика, реализованная в комбинированном КМСФиКЭ, состоит в том, что функционал J(u) выражается через совокупность узловых значений решения на границах СЭ и в подобластях k (рис. 3). При этом, в результате минимизации функционала (6) по полученной совокупности неизвестных узловых значений формируется СЛАУ, соответствующая исходной краевой задаче. Решением полученной СЛАУ служат узловые значения искомой функции на границах СЭ и в подобластях k. После их нахождения определяется решение исходной краевой задачи в каждом СЭ и КЭ. В случае необходимости нахождения решения в точках или на линиях внутри СЭ вводятся дополнительные СЭ, границы которых включают эти точки или линии.

При решении задачи (5) КМСФиКЭ с использованием коэффициентов Фурье сумма интегралов в (6), принимает вид:

Функционал (7) выражается через неизвестные коэффициенты Фурье следов решения на границах СЭ и узловые значения решения на конечно-элементной сетке в области k (рис. 3). Минимизация функционала (7) по данным переменным позволяет сформировать СЛАУ, соответствующую краевой задаче. Отметим, что построенная библиотека «основных» функционалов МСЭФ может быть непосредственно использована для реализации КМСФиКЭ. При этом формирование «основного» функционала в (6) и (7) сводится к простому суммированию «основных»

функционалов в СЭ, выделенных в части расчетной области, в которой материальные свойства среды линейны.

Для оценки возможности применения разработанных методов с помощью КМСФиКЭ был выполнен расчет магнитного поля, создаваемого высококоэрцитивными постоянными магнитами (область 2) с намагниченностью М, размещенными в области 0 с проницаемостью 0, в присутствии ферромагнитных тел i (рис. 4). Задача сводится к минимизации функционала вида:

Рис. 4 Расчетная область цам СЭ части области, в которой размещены постоянные магниты. Функционал (8) выражается через коэффициенты Фурье разложения векторного потенциала по границам СЭ, не примыкающим к нелинейной части области 1 и узловые значения потенциала в нелинейной части области 1, включая ее границу. В результате минимизации функционала по совокупности коэффициентов Фурье и узловым значениям формируется система СЛАУ относительно них. Разработана методика вычисления интегралов, входящих в функционал, в зависимости от свойств сред в СЭ, указаны алгоритмы формирования матриц СЛАУ.

Рассмотренный в этой же главе КМВФиКЭ предполагает использование другого подхода к «склейке» СЭ и формированию СЛАУ, соответствующей краевой задаче. Поиск решения в СЭ в классе функций W22(Gi), а также введение вспомогательных функций позволяет использовать иной, отличный от вариационного подход в части расчетной области, заполненной СЭ. Он основан на том, что при «склейке» СЭ, смежных друг с другом и имеющих общие границы, в КМВФиКЭ используются принципы МВФ. При этом оказываются справедливыми уравнения «связи» вида (3). С помощью данных уравнений формируется СЛАУ относительно коэффициентов разложения следов решения в ряд Фурье на границах «склеиваемых» СЭ.

Вывод уравнений СЛАУ выполняется с учетом условия «склейки» СЭ, граничащих с подобластями k, в которых вводится конечно-элементная сетка. При этом уравнения «склейки» связывают коэффициенты Фурье на общих границах областей Gi и k и узловые значения решения на этих же границах.

Третья часть уравнений СЛАУ формируется с использованием введенной конечно-элементной сетки в областях k. Для решения краевой задачи (6) выполняется минимизация функционала (7) или (8) по введенным в областях i, заполненным ферромагнетиком, узловым значениям решения Ai. Решением сформированной СЛАУ являются узловые значения потенциала на границах СЭ в областях k, также коэффициенты Фурье следов решения на границах СЭ. Решение исходной краевой задачи в каждом СЭ определяется с использованием найденных коэффициентов Фурье следов решения и узловых значений на границах СЭ. В случае необходимости вычисления значений решения в отдельных точках или на линиях в СЭ, можно ввести в СЭ дополнительные СЭ, границы которых содержат эти точки или совпадают с линиями.

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

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

Снижение точности решения таких задач связано с ростом погрешности полиномиальной аппроксимации. Известные способы повышения точности МКЭ предполагают использование в окрестности угловых точек изопараметрических сингулярных элементов или включения в базис сингулярных функций (Митчелл Э., Марчук Г.И.), что практически выполнить сложно. Другой подход (Ю.А. Бахвалова) предполагает использование нерегулярных конечных элементов (НКЭ). Это позволяет сократить порядок решаемой СЛАУ за счет замены большого количества конечных элементов вблизи угловой точки одним элементом, размеры которого определяются размером области. Однако, при большом количестве угловых зон в среде процедура выделения и учета НКЭ является достаточно трудоемкой, а вопрос разработки новых подходов к решению подобных задач в пространственном случае вообще остается открытым.

Выполненные расчеты показали, что использование МСЭФ и МВФ для решения краевых задач при наличии в расчетных областях угловых точек и источников поля, расположенных вблизи них, обеспечивает высокую точность решения.

Для оценки точности МВФ в области, содержащей угловую точку, с особенностью на границах (рис. 5). Расчетная область была заполнена тремя прямоугольными СЭ, причем отношение длин сторон двух из них равно 1/10 (рис. 5).

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

Таблица 2 Значения потенциала в контрольных точках на оси DF нерегулярный КЭ с согласованными КЭ, uНС 10,0149 30,0526 50,1248 70,2806 90, нерегулярный КЭ с несогласованными КЭ, uНН 9,9729 29,8919 49,7714 69,5647 88, Таблица 3 Погрешность расчета потенциала (%) Таблица 4 Погрешность расчета модуля градиента потенциала (%) Оценка точности МВФ была также выполнена в случае, когда в области имеются угловые зоны с источниками поля, размешенные вблизи угловых точек (рис. 6).

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

При решении двухмерной задачи МКЭ было использовано 3784 конечных элемента. Решения КМВФиКЭ находилось в случае, когда расчетная КОС заполнялась 42 СЭ и 160 конечными элементами. Конечно-элементной сеткой был покрыт лишь горизонтальный слой наиболее вытянутых СЭ. Это позволило оценить погрешность, вносимую в расчет МКЭ. Результаты расчетов поля в области с указанной структурой неоднородности, а также содержащей узкие разнородные включения (отношение длин сторон 1/20), осцилляцию решения и особенности в угловых точках, приводятся в табл. 5.

Таблица 5 Погрешности методов расчета для двухмерной области (%) Метод расчета потенциала При решении задачи расчета поля в аналогичным образом заполненной КОС трехмерной области, в соответствии с ее неоднородной структурой, были введены вертикальные слои СЭ в виде параллелепипедов (рис. 7). Результаты расчетов для точек угловой зоны СЭ, смежных с узкими включениями, приведены в табл. 6.

Таблица 6 Погрешность значений потенциала в трехмерной КОС (%) Анализ точности разработанных методов был выполнен также путем сравнения результатов расчетов, полученных МВФ, а также МКЭ, при моделировании областью П-образного вида (рис. 8), форма которой характерна для задач электротехники. Решалась следующая для задачи с точным решением ( n k ) (рис. 9). Рассмотрены случаи отсутствия в области источников ( f 0 ) и их наличии ( f 0 ), а также ситуация, когда в области имеются узкие включения (отношение сторон СЭ 1/20), магнитные проницаемости Рис. 8П-образная область и ее покрытие СЭ полученное МКЭ в программной среде Matlab при очень методы сравнивались с МКЭ по степени дискретизации области, размерностям матриц решаемых СЛАУ, быстродействию, точности и т.д.

Также проведен ряд следующих теоретических исследований:

– для МСЭФ обоснован переход к «основным» функционалам нового типа;

– с использованием теоремы Римана-Лебега проведена оценка сходимости в МСЭФ рядов, аппроксимирующих следы решения на границах СЭ;

– на основании теоремы Римана-Лебега проведена оценка сходимости в МВФ рядов, аппроксимирующих следы решения на границах СЭ.

Отмечается, что коэффициенты ряда Фурье следа функции uk на рассматриваемом участке границе Гk убывают не медленнее, чем n -2, т. е.:

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

В результате численных экспериментов установлено, что при реализации МСЭФ на основе «склейки» СЭ с использованием узловых граничных значений положительный эффект с точки зрения повышения точности достигается при отступе узлов, в которых вычисляются значения искомой функции, от концов сторон «склейки» СЭ. При реализации МВФ более высокая точность расчетов обеспечивается при разложении следов решения в ряд по косинусам (используется чётное продолжение следов решения на виртуальное продолжение соответствующих сторон СЭ). При этом, для нормы остатка ряда, представляющего след решения uk на рассматриваемом участке границы Гk также оказывается справедливой оценка (9).

Результаты численных экспериментов, выполненных при расчете поля в рассмотренных расчетных областях, имеющих аналитическое решение в случаях двух- и трехмерного полей показали, что МВФ и КМВФиКЭ позволяют:

– сократить число элементов разбиения и размерность дискретной модели поля;

– сократить размерность решаемой СЛАУ, соответствующей поставленной краевой задаче;

– обеспечить при нахождении потенциала поля значительно более высокую точность по сравнению с МКЭ даже при наличии в области осцилляции решения, особенностей в угловых зонах, очень узких СЭ (соотношение длин сторон прямоугольных СЭ равно 1/10);

– повысить точность вычисления градиентов потенциала и интегральных характеристик поля за счет получения непрерывно дифференцируемых решений внутри СЭ.

Результаты численных экспериментов, выполненных в условиях специально промоделированной неоднородности среды, при наличии особенности в условиях узких включений, резком изменении материальных характеристик и осцилляции решения показали, что – МВФ и КМВФиКЭ позволяют обеспечить высокую точность решения полевых задач;

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

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

При использовании МВФ для расчета магнитного поля (МП) линейного двигателя (ЛД) с постоянными магнитами, для учета намагниченности при построении СЛАУ на сторонах «склейки» СЭG А (Вб/м) Рис. 11Распределение результаты расчета МП ЛД выполненные с использованием программного пакета FEMM при разбиении расчетной магнитного поля в контрольных точках воздушного зазора между полюсами двигателя иллюстрирует рис. 11.

Рис. 12 – Расчетный виям выполненного физического эксперимента, формулируется так:

где f – заданная функция плотности потерь; Т – функция решения (температура); С4 – температура, поддерживаемая индуктором на границе.

В качестве СЭ были выбраны сектор кольца, сектор кольца с нижней границей. Для разбиения пазовой рованное заполнение расчетной области G конечными руемым полем индуктора. Распределение потерь в стали шихтованного пакета якоря находилось по эмпирическим формулами, полученным экспериментально. Результаты расчета значений температуры КМСФиКЭ в узлах конечно-элементной сетки в стали и обмотке приведены на рис. 14. Отличия экспериментальных (Тмеди, Тстали) и расчетных значений T меди эксперим.

температуры составило 1,9 % в зоне зубца и 1,3 % в Для оценки эффективности разработанных методов расчета поля в средах с мили и микро включениями, для которых точный расчет физических по- Рис. 14 Распределение лей существующими методами затруднен, был вы- температуры в узлах сетки содержит регулярно расположенные вертикальные варьируемыми величинами. В рассматриваемом ный потенциал приближении, природа сложных процессов форми- боковой рования ионов платины (испарение, электролиз и поверхно- сти т.д.) в устройстве не рассматривалась, а качество аноде, боковых дискретных поверхностях устройст- Рис. 15 Катализатор и ва (рис. 15) и подложке катализатора, которое необ- устройство осаждения ходимо было оценить в результате расчета. Сравнительная оценка точности МВФ и МКЭ выполнена в результате решения краевой задачи Дирихле вида:

При отсутствии внутренних источников поля оценка точности методов выполнялась при изменении условий на внешней границе КОС и варьировании геометрических размеров катализатора. Для решения поставленной краевой задачи МВФ в соответствии с кусочно-однородной структурой среды было введено прямоугольных СЭ, которые «склеивались» по внутренним сторонам. Сравнительная оценка точности МВФ и МКЭ выполнена в контрольных точках в окрестностях угловых зон СЭ, смежных с узкими включениями. Для решения поставленной полевой задачи МКЭ использовался пакет pde программной среды Matlab при разбиении расчетной области на 6208 конечных элемента. В качестве эталона принималось решение, полученное при разбиении расчетной области на 397312 КЭ. Результаты расчета потенциала, полученные МВФ, МКЭ и «эталонным» МКЭ для 7и контрольных точек в угловой зоне элемента подложки приведены в табл. 7.

Таблица 7 – Значения потенциала, полученные МКЭ и МВФ

МКЭ МВФ

№ точки «Эталонный» МКЭ (397312КЭ) МКЭ (6208 КЭ) МВФ (28 СЭ) Для оценки возможности сокращения вычислительных затрат и повышения точности МВФ и КМВФиКЭ был выполнен расчет трех- и двухмерных физических полей (индукции В) и подъемной силы магнитного подвеса. Расчеты проводились с использованием:

– разработанного пакета программ SEMF-3D, реализующего МВФ при 72 СЭ в виде параллелепипедов (12540 уравнений СЛАУ);

– пакета FEMM, реализующего МКЭ и разработанного пакета программ SEMF-2D, реализующего МВФ при нелинейных характеристиках частей МС;

Для учета остаточной намагниченности M в уравнение «связи», согласно алгоритму построения матриц СЛАУ краевой задачи, при рассмотрении граней «склейки» СЭ, расположенных в области постоянных магнитов, введено дополнительное слагаемое. Наличие СЭ одинаковых размеров дает возможность многократно использовать найденные элементы матрицы СЛАУ. Это позволяет сократить объем используемых вычислительных ресурсов. Расчетные значения индукции в нескольких контрольных точках, расположенных в воздушном зазоре, приведены в табл. 8.

Таблица 8 Магнитная индукция, найденная МВФ и МКЭ в зазоре L NКЭ=415000; Nуравнений=209173 NСЭ=120; Nуравнений= Из табл. 8 следует, что результаты расчета F(н) трехмерного магнитного поля, выполненные но меньше. Значения подъемной силы при раз- ги- риги личных зазорах, найденные с использованием МВФ, оказались также близкими к «эталонным», мен найденным МКЭ (рис. 17), и в два раза точнее реа пример, методом точечных источников. ные Также был выполнен расчет магнитного МС поля в области, заполненной КОС, с линейными иподъемной силы от величины зазора и нелинейными характеристиками. Известно, СФ что в нелинейных средах насыщение поля приводит уЭ к выравниванию распределения магнитной индукции, ее «вытеснению» из насыщенных зон в ненасыщенные. Это позволяет ограничиться применением МВФ или МСЭФ со СЭ допустимо малых разме- шаге расчета характеристики материала можно считать постоянными. При таком допущении была ре- 1 шена задача расчета стационарного магнитного поля П-образного подъемного модуля (рис. 18). Зависи- Рис. 18 Разбиение области СЭ мость между индукцией и напряженностью В(Н) в и контрольные точки магнитных системах модуля считалась нелинейной.

Намагничивающая сила iW менялась от 0 до значения, обеспечивающего глубокое насыщение магнитопровода. Величина зазора принималась равной 0,001м. Задача расчета магнитного поля подъемного модуля рассматривалась как плоскопараллельная краевая задача относительно векторного потенциала А в кусочнооднородной среде. Расчетная область, не смотря на простую геометрию, содержит угловые зоны, тонкий воздушный зазор, а также области с линейными и нелинейными магнитными характеристиками. Оценка точности расчета выполнялась путем сравнения результатов, полученных МВФ и МКЭ с использованием программного комплекса FEMM. При решении поставленной задачи МКЭ рассматривалось два варианта разбиения расчетной области, при «эталонном»

(16 007 КЭ, 8183 узла) и для сравнения менее точном (4019 КЭ, 2104 узла). Результаты оценки используемых вычислительных ресурсов и точности расчета поля, выполненного с использованием МКЭ и МВФ, в сечении воздушного зазора (рис.18) в контрольных точках приводятся в табл. 9, 10.

Таблица 9 Сравнительная характеристика примененных методов решения Таблица 10 – Значения потенциала в контрольных точках сечения МКЭ (эталон) -0,0160 -0,0156 -0,0150 -0,0135 -0,0037 -0,0029 -0,0023 -0, С учетом нелинейных характеристик части расчетной области была также решена двумерная задача расчета силовых характеристик подъемного модуля КМВФиКЭ и выполнено сравнение рассчитанных значений со значениями, полученными другими авторами с использованием КМГиКЭ, а также с экспериментальными данными. Намагничивающая сила iW обмотки варьировалась в пределах от 6440 до 12880 А. Для разбиении использовалось 25 СЭ и 864 конечных элемента (рис. 19). Расчетные значения подъемной силы Fy, полученные КМВФиКЭ и КМГиКЭ, и ее экспериментальные значения приводятся в табл. 11.

Таблица 11 Расчетные и экспериментальные значения подъемной силы Fy iW (A) Сравнение результатов расчетов температурных, электростатических, магнитных полей и силовых параметров рассмотренных устройств и систем, выполненных разработанными в работе методами, с результатами расчетов этих полей КМГиКЭ, МТИ, МКЭ, «эталонным» МКЭ и экспериментальными данными, показало, что разработанные методы обеспечивают: у 21 22 23 24 элементов разбиения расчетных областей, уменьшение размерности блочно-линейной СЛАУ и необходимых вычислительных ресурсов при сохранении вы- 13 ции типовых размеров СЭ, наличии особенностей в угловых точках КОС, что подтверждает возможность их использования для решения задач расчета физиче- Рис. 19. Разбиение области ских полей в кусочно-однородных линейных и нели- на КЭ и СЭ нейных средах;

– возможность решения широкого круга задач с использованием созданных программных пакетов, реализующих МВФ, МСЭФ, КМСФиКЭ и КМВФиКЭ.

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

В приложениях приводится вывод «основных» функционалов для СЭ из библиотеки СЭ, акты внедрения результатов диссертационных исследований и их использования в учебном процессе СКФУ, НТИ (филиал СКФУ), свидетельства о регистрации алгоритмов и программ.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ И ВЫВОДЫ

1. Разработан численно-аналитический метод МСЭФ, отличающийся новой вариационной постановкой краевых задач, в которой вместо функционала в виде интеграла Дирихле используются «основные» функционалы по границам СЭ, а вместо приближенной полиномиальной аппроксимации решений краевых задач в СЭ используются точные (аналитические) решения в СЭ, позволяющие с повышенной точностью проводить расчеты физических полей в кусочно-однородных линейных и нелинейных средах, с одновременным сокращением необходимых вычислительных ресурсов.

2. Разработан численно-аналитический метод МВФ, отличающийся тем, что решение в СЭ Gi ищется в пространстве W22(Gi), в каждом СЭ вводятся специальные вспомогательные функции, формирование блочно-ленточной СЛАУ, соответствующей решаемой краевой задаче, выполняется без использования вариационного подхода, что обеспечивает повышенную точность расчетов физических полей в кусочно-однородных линейных и нелинейных средах при одновременном сокращении необходимых вычислительных ресурсов.

3. Разработан численно-аналитический комбинированный метод КМСФиКЭ, отличающийся новой вариационной постановкой краевых задач, использованием аналитических решений краевых задач в СЭ, новыми способами формирования блочно-ленточных СЛАУ, позволяющий сочетать достоинства МСЭФ при расчете поля в областях с линейными свойствами среды с достоинствами МКЭ при расчете поля в нелинейных средах, обеспечивающий повышенную точность расчетов физических полей в линейных и нелинейных средах.

4. Разработан численно-аналитический комбинированный метод КМВФиКЭ, отличающийся тем, что поиск решений в СЭ осуществляется в классе функций из W22(Gi), в каждом СЭ вводятся определенные специальным образом вспомогательные функции, обеспечивающие возможность нахождения решения только на границе СЭ в результате решения СЛАУ, соответствующей краевой задаче. Метод сочетает достоинства МВФ при расчете поля в областях с линейными свойствами среды с достоинствами МКЭ при расчете поля в нелинейных средах, что приводит к снижению погрешности расчета физических полей в кусочно-однородных линейных и нелинейных средах.

5. Впервые создана библиотека двух- и трехмерных СЭ, которые задаются геометрией и типом решаемых краевых задач, определены соответствующие им вспомогательные функции, найдены в аналитическом виде уравнения «связи», получены «основные» функционалы, что позволяет формализовать процессы разбиения расчетной области, построения в области блоков из СЭ, нахождения «основного» функционала, соответствующего решаемой краевой задаче и формирования СЛАУ с блочно-ленточной структурой.

6. Разработаны новые методики «склейки» и задания уравнений «связи» для:

– МСЭФ с использованием на границах СЭ следов решения;

– МВФ с использованием вспомогательных функций;

– КМСФиКЭ и КМВФиКЭ с использованием «склейки» СЭ с конечными элементами.

Методики обеспечивают более точную «склейку» нормальных производных на общих границах СЭ, а также повышение гладкости решения.

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

8. Разработана методика линеаризации нелинейных характеристик среды в СЭ на каждом шаге итерации и коррекции разбиения области на СЭ, которая не требует выполнения операции численного дифференцирования для определения векторов поля и обеспечивает возможность использования МСЭФ и МВФ для моделирования физических полей в нелинейных средах без их комбинирования с МКЭ.

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

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

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

12. Основные результаты диссертационной работы отражены в 48 публикациях, в том числе в 17 публикациях в ведущих научных журналах из списка ВАК.

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

ОСНОВНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

Публикации в ведущих научных журналах из списка ВАК 1. Пашковский, А.В. Метод (p,q) связанных операторов в двухфазной модели тепловых процессов электрических машин / А.В. Пашковский // Изв. вузов.

Электромеханика. – 1985. – № 5. – С. 17-23.

2. Пашковский, А.В. Двухфазная модель в идентификации процессов теплообмена в элементах электрических машин и аппаратов / А.В. Пашковский // Изв.

вузов. Электромеханика. – 1988. – № 7. – С. 30-33.

3. Пашковский, А.В. Метод (p,q) связанных операторов в решении дифференциальных уравнений третьего порядка и обратных задач / А.В. Пашковский // Дифференц. уравнения. – 1988. –Т. 24. –№ 8. – С. 1407-1410.

4. Пашковский, А.В. Комбинированный метод конечных элементов для расчета температурных полей электрических машин и его программная реализация / А.В. Пашковский // Изв. вузов. Электромеханика. – 1988. – № 8. – C. 10-15.

5. Пашковский, А.В. Прямоугольный стандартный элемент в моделировании температурных и электромагнитных полей в кусочно-однородных средах / А.В. Пашковский, И.В. Пашковская // Изв. вузов. Электромеханика. – 2003. – № 3. – C. 9Пашковский, А.В. Прямоугольный стандартный элемент с условиями Дирихле на границе / А.В. Пашковский, И.В. Пашковская // Известия вузов. Электромеханика. – 2003. – № 3. – C. 71-73.

7. Пашковский, А.В. «Склеенные» прямоугольные стандартные элементы в решении модельной полевой задачи / А.В. Пашковский, И.В. Пашковская // Изв.

вузов. Электромеханика. – 2007. – № 1. – C. 78-80.

8. Пашковский, А.В. Метод стандартных элементов в решении задач магнитостатики при особенностях в окрестностях угловых точек / А.В. Пашковский // Научно-технические ведомости СПбГПУ. – 2010. – № 1. – C. 103-106.

9. Пашковский, А.В. Решение тестовых полевых задач в кусочнооднородной области методом стандартных элементов / А.В. Пашковский // Научно-технические ведомости СПбГПУ. – 2009. – № 6. – C. 147-151.

10. Пашковский, А.В. Метод стандартных элементов в моделировании стационарного поля в области с П-образным сердечником / А.В. Пашковский, И. В.

Пашковская // Изв. вузов. Электромеханика. –2009. – № 2. – C. 10-12.

11. Пашковский, А.В. Метод стандартных элементов в расчете магнитного поля линейного двигателя с постоянными магнитами / А.В. Пашковский // Научно-технические ведомости СПбГПУ. – 2010. – № 2. – C. 130-136.

12. Пашковский, А.В. Метод стандартных элементов в управлении процессами термообработки электромеханических устройств / А.В. Пашковский // Известия ЮФУ. Технические науки. – 2010. – № 5. – C. 74-79.

13. Пашковский, А.В. Метод стандартных элементов в решении задач электростатики при особенностях в окрестности угловых точек / А.В. Пашковский // Вестник СевКавГТУ. – 2010. – № 3. – C. 125-129.

14. Пашковский, А.В. Повышение точности расчетов магнитного поля и силовых характеристик электромеханических устройств комбинированным методом стандартных и конечных элементов / А.В. Пашковский // Научно-технические ведомости СПбГПУ. – 2010. – № 3. – C. 140-145.

15. Пашковский, А.В. Комбинированный метод стандартных и конечных элементов в расчете магнитного поля и силовых характеристик асинхронного тягово-подъемного модуля / А.В. Пашковский // Информационные системы и технологии. – 2010. – № 4(60). – C. 33-42.

16. Пашковский, А.В. Блочные численно-аналитические методы и новые математические модели в расчете силовых взаимодействий наночастиц / А.В.

Пашковский, В.И. Пашковский // Научно-технические ведомости СПбГПУ. – 2012. – № 4. – C. 39-44.

17. Пашковский, А.В. Блочный численно-аналитический метод вспомогательных функций в расчете магнитного поля в нелинейных средах / А.В. Пашковский, А.Н. Ткачев // Известия вузов. Электромеханика. – 2013. – № 3. – C. 3-7.

Свидетельство о регистрации программ для ПЭВМ 18. Пашковский, А.В. Программа «SEMF» расчета электромагнитных и температурных полей в кусочно-однородных средах на основе метода стандартных элементов [Электронный ресурс] / А.В. Пашковский // Хроники Объединенного Фонда Электронных Ресурсов Наука и образование. – 2010. – № 7. Режим доступа: http://ofernio.ru/portal/newspaper/ofernio/2010/7.doc.

19. Программа расчета электромагнитных полей, создаваемых постоянными магнитами в присутствии ферромагнитных тел, на основе метода стандартных элементов: cвидетельство о гос. регистр. программы для ЭВМ №2010615042 :

Пашковский А.В., –заявка № 2010613485; заявл. 16.06.2010; зарег. в Реестре программ для ЭВМ Федеральной службы по интеллектуальной собственности, патентам и товарным знакам 05.08.2010.

20. А. с. 1527673 CCCР, МКИ Н01F41/06 Способ изготовления электрических катушек / А.В. Пашковский (CCCР). – № 4325891/24-07; заявлено 10.11.1987; опубл. 07.12.89, Бюл. № 45. – 1 с.

21. Пашковский, А.В. Метод стандартных элементов в идентификации процессов теплообмена в пористых средах с невысокой степенью неоднородности / А.В. Пашковский // Материалы 1-ой Российской национальной конференции по теплообмену. – М., 1994. – С. 1087-1088.

22. Пашковский, А.В. Некоторые аспекты повышения надежности машин при их конструировании / А.В. Пашковский // Молодая наука – новому тысячелетию: матер. Междунар. науч.-техн. конф. – Набережные Челны: КамПИ, 1996. – С.

91-92.

23. Пашковский, А.В. Метод стандартных элементов в краевой задаче расчета температурного поля неоднородного тела / А.В. Пашковский, И.В. Пашковская // Математическое моделирование и компьютерные технологии: материалы 4го Всероссийского симпозиума. – Кисловодск, 2000. – С. 54-55.

24. Пашковский, А.В. Теория расчета коэффициентов межфазного теплообмена многофазных диэлектрических материалов / А.В. Пашковский // Теория, методы и средства измерений, контроля и диагностики: материалы Международной науч.-практ. конф. –Новочеркасск, 2000. – С. 71-72.

25. Пашковский, А.В. Математическое и экспериментальное моделирование нестационарных температурных полей / А.В. Пашковский // Физические основы экспериментального и математического моделирования процессов газодинамики и тепломассобмена в энергетических установках: материалы 13 Школысеминара молодых ученых и специалистов. – Санкт-Петербург, 2001. – С. 191-193.

26. Пашковский, А.В. Компьютерное моделирование температурных и электромагнитных полей в кусочно-однородных средах / А.В. Пашковский, И.В.

Пашковская // Сб. тр. 3-й Междун. науч.-техн. конф. – Санкт-Петербург, 2002. –С.

69-73.

27. Пашковский, А.В. Прямоугольные стандартные элементы в моделировании электромагнитных полей и их склейка / А.В. Пашковский // Математические методы в технике и технологиях: материалы трудов XVII Междун. науч.

конф. – Кострома: КГТУ, 2004. – С. 136-138.

28. Токарева, Д. В. Высокоточные методы расчета в создании катализаторов и микрореакторов в экологической направленности катализа / Д.В. Токарева, И.В. Пашковская, А.В. Пашковский // Новые технологии в азотной промышленности: материалы 2-й Общеросс. науч.-техн. конф. – Невинномысск: СевКавГТУ, 2007. – С. 39-42.

29. Пашковский, А.В. Метод стандартных элементов в моделировании полей наложения системы «катализатор-среда-электрод» / А.В. Пашковский // Математическое моделирование, компьютерные и информационные технологии в технике, экономике и образовании: материалы общеросс. научн.-прак. конф. – Невинномысск: СевКавГТУ, 2009. – С. 94-98.

30. Ткачев, А.Н. Программный комплекс на основе комбинированного метода вспомогательных функций и конечных элементов для компьютерного моделирования потенциальных физических полей / А.Н. Ткачев, И.В. Шкуропадский, А.В.

Пашковский // Обозрение прикладной и промышленной математики. – 2011. – Т.18. – Вып. 4. – C. 470-471.



Pages:   || 2 |
 
Похожие работы:

«ОЛЕНЦЕВИЧ Виктория Александровна МЕТОДИЧЕСКОЕ И ПРОГРАММНОЕ ОБЕСПЕЧЕНИЕ АНАЛИЗА И ПРОГНОЗИРОВАНИЯ БЕЗОПАСНОСТИ ФУНКЦИОНИРОВАНИЯ ПОДСИСТЕМ ЖЕЛЕЗНОДОРОЖНОГО ТРАНСПОРТА Специальность 05.13.01 – Системный анализ, управление и обработка информации (промышленность) АВТОРЕФЕРАТ диссертации на соискание учёной степени кандидата технических наук Иркутск – 2014 Работа выполнена в ФГБОУ ВПО Иркутский государственный университет путей сообщения Научный руководитель : доктор технических...»

«Филиппов Алексей Александрович ФОРМИРОВАНИЕ НАВИГАЦИОННОЙ СТРУКТУРЫ ЭЛЕКТРОННОГО АРХИВА ТЕХНИЧЕСКИХ ДОКУМЕНТОВ НА ОСНОВЕ ОНТОЛОГИЧЕСКИХ МОДЕЛЕЙ Специальность 05.13.12 – Системы автоматизации проектирования (промышленность) Автореферат диссертации на соискание ученой степени кандидата технических наук Ульяновск – 2013 Работа выполнена на кафедре Информационные системы в Ульяновском государственном техническом университете. Научный руководитель : кандидат технических наук,...»

«Капустин Дмитрий Сергеевич МОДЕЛИ И АЛГОРИТМЫ ПАРАЛЛЕЛЬНЫХ ВЫЧИСЛЕНИЙ НА ГРАФИЧЕСКИХ ПРОЦЕССОРАХ И ИХ ПРИМЕНЕНИЕ В ПРОГРАММНЫХ СРЕДСТВАХ АВТОМАТИЧЕСКОГО ТЕСТИРОВАНИЯ ГРАФИЧЕСКИХ ПРИЛОЖЕНИЙ Специальность 05.13.11 – Математическое и программное обеспечение вычислительных машин, комплексов и компьютерных сетей АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Санкт-Петербург 2013 2 Работа выполнена на кафедре Автоматика и вычислительная техника в...»

«ЗАГРЕБНЕВА Анна Дмитриевна СТРУКТУРООБРАЗОВАНИЕ В ПОПУЛЯЦИОННЫХ СИСТЕМАХ, ОБУСЛОВЛЕННОЕ ЯВЛЕНИЕМ ТАКСИСА 05.13.18 – математическое моделирование, численные методы и комплексы программ АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук Ростов-на-Дону 2010 Работа выполнена в отделе математических методов в экономике и экологии НИИ механики и прикладной математики им. Воровича И.И. Южного федерального университета, г. Ростов-на-Дону Научный...»

«ЛИБМАН МИХАИЛ СЕРГЕЕВИЧ АЛГОРИТМИЧЕСКОЕ И МЕТОДИЧЕСКОЕ ОБЕСПЕЧЕНИЕ ПОВЫШЕНИЯ ОПЕРАТИВНОСТИ ПОИСКА ДАННЫХ В КОРПОРАТИВНЫХ ИНФОРМАЦИОННЫХ СИСТЕМАХ Специальность: 05.13.01 – Системный анализ, управление и обработка информации (в наук е и промышленности) Автореферат диссертации на соискание ученой степени кандидата технических наук Калуга - 2013 Работа выполнена в Калужском филиале Московского государственного технического университета им. Н.Э. Баумана. кандидат технических наук,...»

«УСОВ СЕРГЕЙ ВЛАДИМИРОВИЧ МЕТОДИКА ПРОВЕРКИ НАЛИЧИЯ ВОЗМОЖНОСТИ НЕСАНКЦИОНИРОВАННОГО ДОСТУПА В ОБЪЕКТНООРИЕНТИРОВАННЫХ СИСТЕМАХ Специальность: 05.13.19 – Методы и системы защиты информации, информационная безопасность АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Омск-2012 Работа выполнена в ФГБОУ ВПО ОмГУ им. Ф.М.Достоевского. Научный руководитель : доктор физико-математических наук, доцент Белим Сергей Викторович Официальные оппоненты :...»

«ГОДОВНИКОВ Евгений Александрович Автоматизированная система исследования алгоритмов идентификации и прогнозирования аварийных состояний в импульсных системах преобразования энергии. Специальность 05.13.06 – Автоматизация и управление технологическими процессами и производствами (промышленность) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Ханты-Мансийск – 2011 2 Работа выполнена в Федеральном государственном бюджетном образовательном...»

«Сачкова Елена Федоровна Методы, алгоритмы и программы приближенного решения задачи управления 05.13.11 Математическое и программное обеспечение вычислительных машин, комплексов и компьютерных сетей 05.13.01 Системный анализ, управление и обработка информации (технические наук и) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Переславль-Залесский 2009 г....»

«Жериков Андрей Валерьевич ПРИМЕНЕНИЕ КВАЗИГИДРОДИНАМИЧЕСКИХ УРАВНЕНИЙ ДЛЯ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ ТЕЧЕНИЙ ВЯЗКОЙ НЕСЖИМАЕМОЙ ЖИДКОСТИ 05.13.18 – Математические моделирование, численные методы и комплексы программ АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук. Москва, Работа выполнена на...»

«БУБНОВ ДМИТРИЙ ВЛАДИМИРОВИЧ РАЗРАБОТКА МЕТОДОВ И АЛГОРИТМОВ ПРИНЯТИЯ УПРАВЛЕНЧЕСКИХ РЕШЕНИЙ В АВТОМАТИЗИРОВАННЫХ МАШИНОСТРОИТЕЛЬНЫХ СИСТЕМАХ 05.13.06 – Автоматизация и управление технологическими процессами и производствами (технические системы) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Москва 2012 г. Работа выполнена в ФГБОУ ВПО Московском государственном технологическом университете СТАНКИН. Научный руководитель : доктор технических...»

«ГУЛИН Артур Игоревич АНАЛИЗ И СИНТЕЗ ЦЕПНЫХ ТРЕХПОЛЮСНЫХ СТРУКТУР ПРЕОБРАЗОВАТЕЛЕЙ СИСТЕМ УПРАВЛЕНИЯ И УСТРОЙСТВ ВЫЧИСЛИТЕЛЬНОЙ ТЕХНИКИ МЕТОДОМ ФУНКЦИЙ ПРЕОБРАЗОВАНИЯ Специальность 05.13.05 – Элементы и устройства вычислительной техники и систем управления АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора технических наук Уфа 2014 Работа выполнена на кафедре телекоммуникационных систем ФГБОУ Уфимский государственный авиационный технический университет Научный...»

«Нигматулин Равиль Михайлович УСТОЙЧИВОСТЬ СТАЦИОНАРНОГО УРОВНЯ ЧИСЛЕННОСТИ ПОПУЛЯЦИИ В ДИСКРЕТНОЙ МОДЕЛИ ПИЕЛОУ С ДВУМЯ ЗАПАЗДЫВАНИЯМИ 05.13.18 математическое моделирование, численные методы и комплексы программ АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук ЧЕЛЯБИНСК – 2008 Работа выполнена на кафедре математического анализа ГОУ ВПО Челябинский государственный педагогический университет. Научный руководитель : доктор...»

«ЗЯЗИН СЕРГЕЙ НИКОЛАЕВИЧ РАЗРАБОТКА РЕШЕНИЙ ПО ИНТЕГРАЦИИ ТЕРРИТОРИАЛЬНО-РАСПРЕДЕЛЕННЫХ СИСТЕМ ГИБДД И СТРАХОВЩИКОВ Специальность: 05.13.13 – Телекоммуникационные системы и компьютерные сети АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Москва – 2008 Работа выполнена в Московском государственном институте электроники и математики на кафедре Научный руководитель : доктор технических наук, профессор Черкасов Александр Сергеевич Официальные...»

«Трифонов Сергей Владимирович ОПТИМИЗАЦИЯ РАБОТЫ МАЛОМОЩНОЙ БЕСПРОВОДНОЙ СЕНСОРНОЙ СЕТИ НА БАЗЕ ЕЁ ИМИТАЦИОННОЙ МОДЕЛИ Специальность 05.13.18 – Математическое моделирование, численные методы и комплексы программ. АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук Москва – 2013 2 Работа выполнена на кафедре вычислительной математики Московского физико-технического института (государственного университета) Научный руководитель : кандидат...»

«Портнов Игорь Сергеевич РАЗРАБОТКА ИНФОРМАЦИОННОЙ СИСТЕМЫ УПРАВЛЕНИЯ ПОТРЕБЛЕНИЕМ ТОПЛИВНОЭНЕРГЕТИЧЕСКИХ РЕСУРСОВ Специальность: 05.13.01– Системный анализ, управление и обработка информации (промышленность) АВТОРЕФЕРАТ Диссертация на соискание ученой степени кандидата технических наук Владикавказ 2009 Работа выполнена в ГОУ ВПО Северо-Кавказский горнометаллургический институт (государственный технологический университет) Научный руководитель : доктор технических наук, доцент...»

«Достовалов Дмитрий Николаевич СПЕЦИФИКАЦИЯ И ИНТЕРПРЕТАЦИЯ МОДЕЛЕЙ ПЕРЕХОДНЫХ ПРОЦЕССОВ В СИСТЕМАХ ЭЛЕКТРОЭНЕРГЕТИКИ 05.13.11 – Математическое и программное обеспечение вычислительных машин, комплексов и компьютерных сетей АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Новосибирск – 2014 Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Новосибирский государственный...»

«ИВАЩУК ОЛЬГА АЛЕКСАНДРОВНА ТЕОРЕТИЧЕСКИЕ ОСНОВЫ ПОСТРОЕНИЯ АВТОМАТИЗИРОВАННЫХ СИСТЕМ УПРАВЛЕНИЯ ЭКОЛОГИЧЕСКОЙ БЕЗОПАСНОСТЬЮ ПРОМЫШЛЕННОТРАНСПОРТНОГО КОМПЛЕКСА Специальность: 05.13.06 – Автоматизация и управление технологическими процессами и производствами (промышленность) АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора технических наук Орел 2009 2 Работа выполнена на кафедре Информационные системы Государственного образовательного учреждения высшего...»

«Грибанова Екатерина Борисовна АЛГОРИТМЫ И КОМПЛЕКС ПРОГРАММ ДЛЯ РЕШЕНИЯ ЗАДАЧ ИМИТАЦИОННОГО МОДЕЛИРОВАНИЯ ОБЪЕКТОВ ПРИКЛАДНОЙ ЭКОНОМИКИ Специальность 05.13.18 Математическое моделирование, численные методы и комплексы программ АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Томск – D Работа выполнена в Томском государственном университете систем управления и радиоэлектроники. Научный руководитель : доктор технических наук, профессор Мицель Артур...»

«Максаков Алексей Владимирович ПОВЫШЕНИЕ РЕЛЕВАНТНОСТИ ПЕРИОДИЧЕСКОГО ТЕМАТИЧЕСКОГО ПОИСКА ИНФОРМАЦИИ В WEB Специальность 05.13.11 – математическое и программное обеспечение вычислительных машин, комплексов и компьютерных сетей АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук МОСКВА 2007 Работа выполнена на кафедре автоматизации...»

«Скворцова Мария Ивановна МАТЕМАТИЧЕСКИЕ МОДЕЛИ И АЛГОРИТМЫ В ИССЛЕДОВАНИЯХ СВЯЗИ МЕЖДУ СТРУКТУРОЙ И СВОЙСТВАМИ ОРГАНИЧЕСКИХ СОЕДИНЕНИЙ 05.13.18 – математическое моделирование, численные методы и комплексы программ АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора физико-математических наук Москва – 2007 1 Работа выполнена в Московской государственной академии тонкой химической технологии (МИТХТ) им. М. В. Ломоносова ОФИЦИАЛЬНЫЕ ОППОНЕНТЫ: доктор...»






 
© 2013 www.diss.seluk.ru - «Бесплатная электронная библиотека - Авторефераты, Диссертации, Монографии, Методички, учебные программы»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.