Конечноэлементное моделирование гармонических электромагнитных полей тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат технических наук Рояк, Светлана Хаимовна

• Рояк, Светлана Хаимовна
• кандидат технических наук
• 2000

Оглавление диссертации кандидат технических наук Рояк, Светлана Хаимовна

Введение.

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

1.1. Электромагнитное поле в двумерной среде, возбуждаемое аксиальным источником тока.

1.2. Осесимметричное электромагнитное поле, возбуждаемое гармоническим круговым током в однородной по магнитной проницаемости среде.,.

1.3. Осесимметричное электромагнитное поле, возбуждаемое удаленным от проводящей среды низкочастотным гармоническим током в круговой петле.

1.4. Трехмерное электромагнитное поле, возбуждаемое низкочастотным гармоническим источником в однородной по магнитной проницаемости среде.

1.4.1. Математическая модель.

1.4.2. Вычисление матриц и векторов правых частей конечноэлементной СЛАУ.

1.4.3. Учет условий симметрии.

1.5. Форматы хранения конечноэлементных матриц.

2. Итерационные методы решения конечноэлементных СЛАУ и ускорение их сходимости на основе неполной факторизации матриц.

2.1. Структура матриц конечноэлементных СЛАУ, аппроксимирующих трехмерные задачи.

2.2. Использование неполной факторизации для предобусловливания СЛАУ.

2.2.1. Алгоритм построения ¿[/(^-факторизации.

2.2.2. Алгоритм построения неполной блочной Ш^ц)-факторизации.

2.3. Учет симметрии матрицы конечноэлементной СЛАУ при построении матрицы неполного разложения.

2.4. Проблема существования неполного ¿[/(^-разложения.

2.5. Блочно-диагональное Ь и (б д)-разложение.

2.6. Итерационные методы решения СЛАУ с несимметричными матрицами.

2.6.1. ОМШ-Э.

2.6.2. Локально оптимальная схема.

2.6.3. Метод сопряженных градиентов для симметризованной СЛАУ.

2.6.4. Метод бисопряженных градиентов.

2.6.5. Метод бисопряженных градиентов стабилизированный.

2.6.6. Метод последовательной верхней релаксации.

3. Использование предложенных методик моделирования гармонических электромагнитных полей при решении практических задач.

3.1. Математическое моделирование электромагнитного поля в трехфазном кабеле.

3.2. Исследование электромагнитного поля в пазах двухклеточного ротора асинхронного короткозамкнутого двигателя.

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

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

3.4.1. Сравнение с результатами, полученными разностным методом и экспериментально.

3.4.2. Трехмерные задачи каротажа.

3.4.3. Сравнительный анализ сходимости различных итерационных методов решения СЛАУ в задачах индукционного каротажа.

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

3.5.1. Методика моделирования.

3.5.2. Сравнительный анализ сходимости различных итерационных методов решения СЛАУ в задачах аэроразведки.

3.5.3. Результаты исследования, проведенного с помощью пакета прикладных программ HAREM.

4. Применение разработанного подхода в программном комплексе TELMA.

4.1. Структура программного комплекса TELMA.

4.2. Структуры данных программного комплекса TELMA.

4.3. Программные модули для генерации и решения конечноэлементных СЛАУ при моделировании гармонических электромагнитных полей.

4.4. Подсистема HAREM программного комплекса TELMA.

4.4.1. Возможности препроцессоров комплекса HAREM.

4.4.2. Принципы конечноэлементной аппроксимации.

4.4.3. Управляющая оболочка программного комплекса

HAREM.

4.4.4. Расчетные модули программного комплекса HAREM.

4.4.5. Представление результатов.

Введение диссертации (часть автореферата) на тему «Конечноэлементное моделирование гармонических электромагнитных полей»

Математическое моделирование электромагнитных полей является мощным инструментом теоретических исследований и широко используется в различных областях науки и техники. Среди всех задач моделирования электромагнитных полей можно выделить широкий спектр задач, в которых электромагнитное поле возбуждается гармоническими источниками тока. Математическому моделированию таких полей посвящено большое количество работ [1, 16, 21, 22, 26, 33-35, 41, 42, 53, 63, 64, 81, 84, 109, 110, 112, 118-120].

Довольно часто для вычисления характеристик гармонических электромагнитных полей применяют аналитические методы, использующие некоторые упрощения математической модели поля [26, 53, 63, 109, 119, 120]. Основное достоинство таких методов состоит в том, что предлагаются достаточно простые с точки зрения вычислений формулы, описывающие требуемые характеристики поля. Недостатком этих методов является невозможность учета всех закономерностей поведения электромагнитного поля и необходимость в каждом частном случае разрабатывать новые подходы. Методы численного моделирования, рассматриваемые в дальнейшем, являются с этой точки зрения более универсальными, хотя в тех частных случаях, когда аналитические методы применимы, могут проигрывать им из-за более высокой сложности получения конечного результата.

На основе изучения электромагнитных полей, возбуждаемых гармоническими источниками, разрабатываются многие индукционные методы электроразведки [1, 4, 27, 34, 38, 63, 69, 84, 105, 109, 117, 118-120]. При использовании гармонических источников могут измеряться либо абсолютные характеристики поля (амплитуды и фазы декартовых компонент поля), либо относительные (отношения амплитуд и сдвиги фаз одноименных компонент поля в двух точках пространства или разноименных компонент в одной точке). Вместо амплитуды и фазы поля могут измеряться его действительная и мнимая части (квадратурное или синфазное с током в возбудителе напряжения, снимаемые с приемных рамок) [117]. Кроме того, могут изучаться интегральные характеристики (э.д.с. в контуре, совмещенном с возбуждающей поле петлей или отнесенным от нее на некоторое расстояние). Методы гармонических полей применяются в наземном, скважинном и аэровариантах. Характерной особенностью методов аэроэлектроразведки является то, что в силу большой удаленности приемника поля от возбудителя (по сравнению с высотой полета) величины наблюдаемых в этих методах аномальных эффектов от проводящей среды становятся сравнимы с аномалиями, регистрируемыми при наземных исследованиях.

Первых заметных успехов при решении дву- и трехмерных задач геоэлектромагнетизма удалось достичь исследователям, использовавшим для проведения численных расчетов метод интегральных уравнений (МИУ) в различных его вариантах [1, 21, 33, 41, 42, 64, 81, 84, 112, 117]. Однако один из главных недостатков МИУ, заключающийся в необходимости решения получающейся в результате численной аппроксимации интегрального уравнения системы линейных алгебраических уравнений (СЛАУ) с плотной матрицей, очень часто не позволяет выполнить дискретизации на уровне, необходимом для получения численного решения заданной точности [33, 42]. Из-за заполненности матрицы СЛАУ, аппроксимирующей интегральное уравнение, вычислительные затраты на решение задачи растут с измельчением сетки для МИУ столь быстро, что даже для компьютеров большой мощности становится актуальной проблема нехватки вычислительных ресурсов. Подтверждение этого можно найти в приведенном в работе [33] сравнительном анализе различных методов на примере решения трехмерных задач геоэлектромагнетизма в проекте СОММЕМ1. Представленные результаты доказывают, что различия в аномальных составляющих поля практически для всех трехмерных моделей во многих точках съема достигают 50% и более даже для тех программ, которые реализованы на основе МИУ. При этом отмечается, что расчеты проводились на сетках с таким числом ячеек дискретизации, что получаемые СЛАУ требовали объема вычислительных ресурсов на пределе возможностей использованных ЭВМ.

Некоторым развитием МИУ можно считать подход, изложенный в [1, 41, 81] и основанный на построении решения интегрального уравнения в виде функционального ряда Неймана. Этот подход позволяет заметно сократить (по сравнению с классической реализацией МИУ) требуемые вычислительные ресурсы, в результате чего появляется возможность использовать гораздо более подробные сетки. Но оценить реальные возможности этого метода при решении трехмерных задач пока довольно трудно, так как приведенные в [1, 41, 81] результаты несколько ограничены и не позволяют получить более или менее полного представления о точности и надежности метода при решении трехмерных задач. Также из [1, 41, 81] остается неясной зависимость скорости сходимости аппроксимации решения к точным значениям от величины ячейки дискретизации области. Кроме того, и МИУ, и разрабатываемые на его основе подходы при наличии в задаче негоризонтальных слоев, наклонных контактов и других трудностей, связанных с усложнением геометрии, а также в случае использования сеток с ячейками неодинакового объема становятся слишком затратными по памяти.

Все эти трудности, связанные с применением МИУ и близких к нему методов, заставляют искать новые подходы к решению трехмерных задач геоэлектромагнетизма, основанные на таких универсальных методах решения дифференциальных уравнений в частных производных, как методы конечных разностей (МКР) и конечных элементов (МКЭ) [4, б, 7, 14, 18, 28, 31, 39, 40, 43, 54, 59, 60, 65, 66, 70-78, 82, 85-87, 89-91, 94-98, 100-104, 107, 108, 113, 115]. МКЭ, позволяющий использовать симплициидальные (т.е. треугольные в двумерном случае и тетраэдральные в трехмерном) сетки, имеет значительные преимущества перед МКР благодаря возможностям построения нерегулярных сгущающихся или разреживающихся сеток с минимальным количеством так называемых «лишних» узлов (т.е. узлов, удаление которых из сетки не приводит к увеличению погрешности приближенного решения). МКР также серьезно уступает МКЭ по возможностям учета сложной трехмерной геометрии исследуемой среды: негоризонтальных или выклинивающихся слоев, наклонных контактов и т.д.

Однако применение стандартных вычислительных конечноэлемент-ных схем для решения трехмерных задач геоэлектромагнетизма также требует очень больших вычислительных затрат [16, 33]. Это приводит либо к большим погрешностям решения из-за использования недостаточно подробной сетки, либо к чрезмерно большой стоимости получения решения важных практических задач. Предложенные в [14, 39, 71-75, 96-98, 100-103] подходы к конечноэлементному моделированию с выделением двумерной части поля позволяют при . решении многих трехмерных задач геоэлектромагнетизма во много раз снизить вычислительные затраты и тем самым сделать эти задачи доступными для решения с очень высокой точностью при относительно небольших вычислительных затратах. В предлагаемой диссертационной работе будут построены основанные на таких подходах конечноэлементные схемы моделирования электромагнитных полей с гармоническими источниками.

Получаемые при конечноэлементной аппроксимации СЛАУ имеют разреженные матрицы. Методы решения таких СЛАУ рассматривались и исследовались в [2, 3, 5, 8-12, 15, 17, 20, 29, 32, 36, 43-51, 56-58, 79, 80, 114]. Но при этом разрабатываемые методы, как правило, не учитывают специфические особенности матриц СЛАУ, заключающиеся в существенных различиях между объемом памяти, необходимом для хранения этих матриц при стандартной их упаковке (например, при использовании разреженного строчного или диагонального форматах хранения), и объемом памяти, требуемым для хранения этих матриц при специальной упаковке. Заметим, что для некоторых классов задач применение специальных упаковок может на порядок сократить объем памяти под хранение матриц даже по сравнению с такими довольно эффективными универсальными форматами, как разреженный строчный или диагональный. К такому классу задач относятся и задачи конечноэлементного моделирования трехмерных электромагнитных полей с гармоническими источниками.

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

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

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

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

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

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

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

2. Предложены специальные структуры данных для экономичного хранения матриц конечноэлементных СЛАУ и матриц предобусловлива-ния.

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

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

Практическая ценность работы и реализация результатов.

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

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

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

Апробация работы. Основные результаты работы были представлены на IV Международной конференции «Актуальные проблемы электронного приборостроения» АПЭП-98 (Новосибирск, 1998.); Третьем сибирском ■ конгрессе по прикладной и индустриальной математике, посвященным памяти С.Л.Соболева (Новосибирск, 1998.); The Third Russian-Korean International Symposium on Science and Technology KORUS'99 (Новосибирск, 1999.), семинарах НГТУ. Результаты проведенных исследований включались в отчеты по НИР НГТУ.

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

Структура работы. Диссертация состоит из введения, четырех глав, заключения, списка использованных источников (120 наименований) и 4 приложений. Работа изложена на 206 страницах, включая 37 иллюстраций и 35 страниц приложений.

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Рояк, Светлана Хаимовна

ЗАКЛЮЧЕНИЕ

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

1. Построена конечноэлементная схема для моделирования трехмерного электромагнитного поля, возбуждаемого гармоническим током в круговой петле. Эта схема базируется на методике разделения искомого поля ~ на основное (осесимметричное) и добавочное (поле влияния трехмерных объектов). Такой подход дает возможность решать задачи моделирования электромагнитного поля в средах с трехмерными неоднородностями различной формы при достаточно небольших вычислительных затратах.

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

3. Предложены структуры данных, позволившие существенно сократить затраты памяти на хранение матриц конечноэлементных СЛАУ и матриц . предобусловливания.

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

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

6. Проведено исследование влияния токов смещения на примере решения осесимметричных задач электроразведки и показано, при каких частотах влияние токов смещения становится несущественным.

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

161 подсистему HAREM программного комплекса TELMA, которая применялась при решении практических задач геофизики.

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

Список литературы диссертационного исследования кандидат технических наук Рояк, Светлана Хаимовна, 2000 год

1. Avdeev D.B., Kuvshinov A.V., Pankratov O.V., Newman G.A. HighPerformance Three-Dimensional Electromagnetic Modelling Using Modeled Newman Series. Wide-Band Numerical Solution and Examples / /J. Geomag. Geoelectr. 1997. - V.49, №11-12. - P. 1519-1539.

2. Fletcher R. Conjugate gradient methods for indefinite systems / / Lecture Notes in Mathematics. 1976. - V.506. - P.73-89.

3. Lanczos C. Solution of systems of linear equations by minimized itera-tione / /J. Res. Nat. Bur. Standards 1952. - V.49, №1. - P. 33-53.

4. Peters A. Non-symmetric CG-like schemes and the finite element solution of the advection-dispersion equation // International Journal for Numerical Methods in Fluids. 1993. - Vol.17. - P.955-974.

5. Saad Y. Varietions of Arnoldi's method for computing eigenelements of large unsymmetric matrices// Linear Algebra Appl. 1980. - V.34. -P.269-295.

6. Saad Y., Iterative methods for sparse linear systems. NY: Intern. Thomson Publishing, 1998.

7. Saad Y., Krylov subspace methods for computing eigenelements of large unsymmetric matrices// Linear Algebra Appl. 1980. - V.34. - P.269-295.

8. Saad Y., Krylov subspace methods on supercomputers// SIAM J. Sci. Stat. Comput. 1989. - V.10. - P. 1200-1232.

9. Saad Y., Shultz M.H. GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems / / SIAM J.Sci.Stat.Comput. -1986. V.7. - P.856-869.

10. Soloveichik Y.G. Iterative method for solving finite element systems of algebraic equations / / Computers Math. Applic. V.33, N6, 1997. -P.87-90.

11. Van der Vorst H.A. Bi-CGStab: A fast and smoothly converging variant of Bi-CG for the solution nonsymmetric linear systems: Preprint 633. -Univ. Utrecht, 1990.

12. Weidelt P. Electromagnetic induction in three-dimensional structure //J. Geophys., 1975, V.41, P.85-109

13. Xu J., Cai X.-C. A preconditioned GMRES method for nonsymmetric or indefinite problems //Mathematics of computation. 1992. - V.59, №200. - P.311-319.

14. Андреев В.Б. Лекции по методу конечных элементов: Учеб. Пособие.- М.: Изд-во МГУ, 1997. 178 с.

15. Асинхронные двигатели общего назначения /Е.П.Бойко, Ю.В.Гаинцев, Ю.М.Ковалев и др./ под ред. В.М.Петрова и А.Э.Кравчина. М.: Энергия, 1980. - 488 с.

16. Блатов И.А. Об оценках ¿(/-разложения разреженных матриц и их приложения к методам неполной факторизации //ЖВМ и МФ -1997/ Т.З, №3. - С.259-276.

17. Вержбицкий В.В. Прямые задачи электрокаротажа //Изв. РАН, Сер. Физика Земли. 1997. - №3. - С.71-74.

18. Вешев А.В. Электропрофилирование на постоянном и переменном токе. Л.: Недра, 1980. - 391с.

19. Воеводин В.В. О методах сопряженных направлений //ЖВМ и МФ.- 1979. № 15. С.1313-1317.

20. Воеводин В.В. Численные методы алгебры. М.:Наука, 1966.

21. Воеводин В.В., Кузнецов Ю.А. Матрицы и вычисления. М.:Наука, 1984. - 319 с.

22. Войтович P.A., Кадомская К.П. Продольные погонные параметры и потери в оболочках в кабелях трехфазного исполнения/ / Научный вестник НГТУ. Новосибирск, НГТУ, 1995г., №1, - с59-67.

23. Вычислительные математика и техника в разведочной геофизике: Справочник геофизика / Под ред. В.И.Дмитриева. М.: Недра, 1990. - 498 с.

24. Галлагер Р. Метод конечных элементов. Основы. М., 1984. - 428 с.

25. Голуб Дж., Ван Лоун Ч. Матричные вычисления: Пер. с англ. М.: Мир, 1999. - 548 с.

26. Горелик Л.В., Дарьин С.Г., Голубков В.П. Анализ существующих конструкций зубцовой зоны ротора асинхронного двигателя и методов расчета его параметров в пусковых режимах / / Информэлектро. -Владимир, 1983. 38 с.

27. Деклу Ж. Метод конечных элементов. М., 1976. - 95 с.

28. Джордж А., Лю Дж. Численное решение больших разреженных систем уравнений. М.: Мир, 1984. 333 с.

29. Дмитриев В.И. Методы моделирования электромагнитных полей.: Материалы междунар. проекта СОМЕМ1/ М.С. Жданов, И.М. Варенцов, К.Г. Голубев, В.А. Крылов. М.: Наука, 1990. - 198 с.

30. Дмитриев В.И., Барышникова И.А., Захаров Е.В. Аномальные электромагнитные поля пластовых тел. Л.: Недра, 1977. - 168 с.

31. Дмитриев В.И., Захаров Е.В. Интегральные уравнения в краевых задачах электродинамики.: Учеб. пособие. М.: Изд-во МГУ, 1987. -167 с.

32. Доледенок O.A., Ильин В.П. О. скорости сходимости метода неполной факторизации для диффузионно-конвективных уравнений //Тр. Вычислительного центра СО РАН. Сер. Вычислительная математика. -Новосибирск, 1995. Вып. 3. - С.41-51.

33. Дьяков Е.Г. Минимизация вычислительной работы. М.: Наука, 1989.

34. Жданов М.С. Электроразведка: Учебник для Вузов. М.: Недра, 1986. - 316 с.

35. Зенкевич О., Морган К. Конечные элементы и аппроксимация. М., 1986.- 318 с.

36. Зингер Б.Ш., Панкратов О.В., Файнберг Э.Б. Пленочное моделирование поверхностных и глубинных неоднородностей //Изв. РАН, Сер. Физика Земли. 1992. - №10. - С.93-108.

37. Зингер Б.Ш., Файнберг Э.Б. Электромагнитная индукция в неоднородных тонких слоях. М.: ИЗМИРАН, - 1985. - 234 с.

38. Ильин В.П, Юдин А.Н. Модификация факторизованной матрицы системы трехмерных разностных уравнений/ /Вычислительные методы и технология решения задач математической физики. Новосибирск, 1993. - С.125-135.

39. Ильин В.П. Методы неполной факторизации для решения алгебраических систем. М.:Физматлит, 1995. - 288 с.

40. Ильин В.П. Численные методы решения задач электрофизики. М.: Наука, 1985. - 336 с.

41. Ильин В.П., Ицкович Е.А. Некоторые варианты неявных методов неполной факторизации //Тр. Вычислительного центра СО РАН. Сер. Вычислительная математика, Новосибирск, 1995. - Вып. 5. - С.93-113.

42. Ильин В.П., Ицкович Е.А. Трехсеточный метод решения двумерных краевых задач //Тр. Вычислительного центра СО РАН. Сер. Вычислительная математика. Новосибирск, 1996. - Вып. 5. - С.44-68.

43. Ильин В.П., Ицкович Е.А. Экспериментальный анализ явных методов неполной факторизации / / Численные методы и математическое моделирование. Новосибирск, 1990. - С.85-94.

44. Ильин В.П., Карначук В.И., Ларин М.Р. Древовидный подход к организации структуры данных для разложения Холесского / / ЖВМ и МФ 1994. - Т.34, №12. - С.1747-1756.

45. Ильин В.П., Ларин М.П. Многоуровневый итерационный метод неполной факторизации решения пятиточечных систем уравнений / / Тр. Вычислительного центра СО РАН. Сер. Вычислительная математика. Новосибирск, 1996. - Вып. 5,— С.69-89.

46. Ильин В.П., Юдин А.Н. Решение трехмерных разностных уравнений методом Булеева с сопряженными градиентами / / Технология моделирования задач математической физики/Сб. Научн. тр.; Под ред. В.П.Ильина. Новосибирск: ВЦ СОРАН, 1989. - С.152-165.

47. Инкин А.И., Темлякова З.С. Метод расчета комплексного сопротивления зубцового деления ротора с учетом насыщения зубцов / / Электричество. 1997. - №7. - С.37-42.

48. Инкин А.H., Рейхердт A.A. Математическая модель для расчета электромагнитных процессов в трехфазных кабелях с проводящей оболочкой // Электричество. 1999. - №5. - С.28-34.

49. Каледин В.О., Ластовецкий В.П. Решение прямой задачи электроразведки постоянным током методом конечных элементов / / Изв. РАН, Сер. Физика Земли. 1988. - №12. - С.31-38.

50. Каменецкий Ф.М. Электромагнитные геофизические исследования методом переходных процессов. М.: ГЕОС, 1997. - 162 с.

51. Капорин И.Е. О предобусловливании метода сопряженных градиентов при решении дискретных аналогов дифференциальных задач / /Дифференциальные уравнения. 1990. - Т.26, №7. - С. 12251236.

52. Коротков Д.Ю. Методы приближенной факторизации для решения уравнений эллиптического типа. М.: ВЦ АН СССР, 1989.

53. Кузнецов Ю.А. Метод сопряженных градиентов, его обобщения и применения. Вычислительные процессы и системы. М.: Наука, 1983. -Вып. 1.

54. Кулон Ж.-Л., Сабоннадьер Ж.-К. САПР в электротехнике: Пер. с франц. М.: Мир, 1988. - 208 с.

55. Лаевский Ю.М. Метод конечных элементов решения многомерных параболических уравнений. Новосибирск, 1993. - 103 с.

56. Ландау Л.Д., Лившиц Е.М. Теоретическая физика. М.: Наука, 1992.- Т.8. Электродинамика сплошных сред. 664с.

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

58. Мартаков C.B., Эпов М.И.Прямые задачи электромагнитного каротажа // Геология и геофизика. 1999. - Т.40, №2. - С.249-254.

59. Мартышко П.С. Об интегральных преобразованиях электромагнитных полей // Изв. РАН. Сер. Физика Земли. 1997. - №2,- С.69-70.

60. Марчук Г.И. Методы вычислительной математики. М.: Наука, 1989.- 608 с.

61. Марчук Г.И., Агошков В.И. Введение в проекционно-сеточные методы.- М.: Наука, 1981. 416 с.

62. Марчук Г.И., Кузнецов Ю.А. Теория и применение обобщенного метода сопряженных градиентов. Новосибирск: ВЦ СО АН СССР, 1980,

63. Матвеев Б.К. Электроразведка: Учеб. пособие для вузов. М.: Недра, 1990. - 368 с.

64. Митчел Э., Уэйт Р. Метод конечных элементов для уравнений с частными производными. М., 1981. - 216 с.

65. Моисеев B.C., Соловейчик Ю.Г., Рояк М.Э. Математическое моделирование сложнопостроенных сред / /Сб. рефератов №2 Международной геофизической конференции и выставки по разведочной геофизике SEG-EAGO. М., 1993. - С.15.

66. Молчанов И.Н., Николаенко Л.Д. Основы метода конечных элементов. Киев: Наук, думка, 1989. - 272 с.

67. Норри Д., де Фриз Ж. Введение в метод конечных элементов. М.: Мир, 1981. - 304 с.

68. Обэн Ж.-П. Приближенное решение эллиптических краевых задач. -М.: Мир, 1977. 383 с.

69. Оганесян Л.А., Руховец Л.А. Вариационно-разностные методы решения эллиптических уравнений. Ереван: Изд-во АН Арм.ССР, 1979.- 192 с.

70. Ортега Дж. Введение в параллельные и векторные методы решения линейных систем. М.: Мир, 1991. - 367 с.

71. Панкратов О.В., Авдеев Д.Б., Кувшинов А.В. Рассеяние электромагнитного поля в неоднородной земле. Решение прямой задачи. //Изв. РАН. Сер. Физика Земли. 1995. - №3.- С. 17-25.

72. Пересветов В. В. Математическое моделирование методом конечных элементов магнитотеллурических полей в двумерных средах с криволинейными границами. Препр. - Владивосток: ДВО АН СССР. -1990. - 33.

73. Писсанецки С. Технология разреженных матриц. М.: Мир, 1988. -410 с.

74. Прямые и обратные задачи геоэлектрики /Ин-т земного магнетизма, ионосферы и распространения радиоволн АН СССР. М.: Наука, 1990. - 101с.

75. Рояк М.Э., Соловейчик Ю.Г. Алгоритмы построения нерегулярных треугольных и тетраэдральных сеток // Сб. научн. тр. НГТУ. Новосибирск: НГТУ, 1996. - №2(4). - С.39-46.

76. Рояк М.Э., Соловейчик Ю.Г., Шурина Э.П. Сеточные методы решения краевых задач математической физики: Учеб. пособие. Новосибирск: Изд-во НГТУ, 1998. - 120 с.

77. Сабоннадьер Ж.-К., Кулон Ж.-Л. Метод конечных элементов и САПР: Пер. с фр. М.: Мир, 1989. - 190 с.

78. Самарский А.А, Николаев Е.С. Методы решения сеточных уравнений. М.: Наука, 1978. - 560 с.

79. Самарский A.A. Введение в численные методы. М.: Наука, 1997. -239 с.

80. Сегерлинд Л. Применение метода конечных элементов. М.: Мир,- . 1979. - 392 с.

81. Сильвестер П., Феррари Р. Метод конечных элементов для радиоинженеров и инженеров-электриков.- М.: Мир, 1986. 229 с.

82. Соловейчик Ю.Г. Вычислительные схемы МКЭ-моделирования трехмерных электромагнитных и тепловых полей в сложных областях: Ав-тореф. дис. . докт. техн. наук. Новосибирск, НГТУ, 1997.

83. Соловейчик Ю.Г., Рояк М.Э. Расчет трехмерного нестационарного электромагнитного поля с учетом вихревых токов / /Сб. научн. тр. НГТУ. Новосибирск: НГТУ, 1996. - №3(5). - С.71-80.

84. Соловейчик Ю.Г., Рояк М.Э., Моисеев B.C., Васильев A.B. Математическое моделирование на базе метода конечных элементов трехмерных электрических полей в задачах электроразведки //Изв. РАН. Сер. Физика Земли. 1997. - №9. - С.67-71.

85. Соловейчик Ю.Г., Рояк М.Э., Моисеев B.C., Тригубович Г.М. Моделирование нестационарных электромагнитных полей в трехмерных средах методом конечных элементов //Изв. РАН. Сер. Физика Земли. 1998. - №10. - С.78-83

86. Соловейчик Ю.Г., Рояк М.Э., Рояк С.Х., Тригубович Г.М. Применение МКЭ для расчета трехмерных гармонических электромагнитных полей в задачах каротажа и аэроразведки полезных ископаемых //Науч. вестник НГТУ. Новосибирск: НГТУ, 1998. - №1. - С.146-160

87. Юб.Сохранов Н.Н., Аксельрод С.М. Обработка и интерпретация с помощью ЭВМ результатов геофизических исследований нефтяных и газовых скважин. М.: Недра, 1984. - 255 с.

88. Справочник по электрическим машинам: в 2 т. / Под общ. ред. И.В.Копылова и Б.К.Клокова. М.: Энергоатомиздат, 1988,- Т.1. -456 с.

89. Стренг Г., Фикс Дж. Теория метода конечных элементов. М.: Мир, 1977. - 350 с.

90. Сьярле Ф. Метод конечных элементов для эллиптических задач. -М.: Мир, 1980. 512 с.

91. Табаровский Л.А., Эпов М.И. Электромагнитные поля гармонических источников в слоистых анизотропных средах / / Геология и геофизика. 1977. - №1. - С.101-109.

92. ПО.Уэйт Дж. Р. Геоэлектромашитизм: Пер. с англ./Ред. М.Н. Бердичев-ский. М.: Недра, 1987. - 235 с.

93. Ш.Фаддеев Д.К., Фаддеева В.Н. Вычислительные методы линейной алгебры. М., Л.: Физматгиз, 1963.

94. Файнберг Э.Б. Глобальное и региональное магнитовариационное зондирование Земли. Автореф. дис. . докт. техн. наук.- Троицк. Моск. область, ИЗМИРАН, 1983г.

95. ПЗ.Флетчер К. Численные методы на основе метода Галеркина. М.: Мир, 1988. - 352 с.

96. Хейгеман Л., Янг Д. Прикладные итерационные методы. М.: Мир, 1986.

97. Шайдуров В.В. Многосеточные методы конечных элементов. М.: Наука, 1989. - 288 с.

98. Шуйский В.П. Расчет электрических машин. Л.: Энергия, 1968. -731с.171

99. Электроразведка: Справочник геофизика. В 2-х кн. /Под ред. В.К. Хмелевского и В.М. Бондаренко. М.: Недра, 1989. - Кн. 1. - 438 с.

100. Эпов М.И. Электромагнитные методы исследования скважин. Новосибирск: Наука, 1979. - 104 с.

101. Эпов М.И., Никитенко М.И. Система одномерной интерпретации данных высокочастотных индукционных каротажных зондирований //Геология и геофизика. 1993. - Т.34, №2. - С.124-130.

102. Эпов М.И., Сухорукова К.В.,Никитенко М.Н., Антонов Ю.Н. Особенности высокочастотных индукционных зондирований в скважинах с горизонтальным завершением //Геология и геофизика. 1998. -Т.39, №5. - С.649-656.