Математический аппарат в моделях разных иерархических уровней

 

Кматематическому обеспечению (МО) САПР относятся:

· математические модели,

· численные методы,

· алгоритмы выполнения проектных процедур.

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

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

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

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

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

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

 

Требования к математическим моделям и численным методам в САПР

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

Адекватность имеет место, если модель отражает заданные свойства объекта сприемлемойточностью.Точность математической модели заключается в максимальном учете всех факторов, присутствующих в реальном физическом процессе, чтобы математическая модель позволяла моделировать этот процесс. Адекватность МО относительна, так как оценить ее мы можем только с позиций современного состояния средств измерения. Например в 16 веке Галлилео Галлилей открыл закон свободного падения, в котором описал влияние на падающее тело притяжения Земли. В 17 веке Исаак Ньютон сформулировал закон всемирного тяготения, где учел взамное влияние двух тел как силу прямопропорциональную их массам и обратнопропорциональную квадрату расстояния между ними. Этот закон уточнял результат Галилея как частный случай но носил качественный характер и опиралсяь на результаты экспериментов, выполненных на тогдашней лабораторной базе. Только в 19 веке эксперименты Кавендиша и Пуассона дали возможность вычислить в формуле Ньютона гравитационную постоянную G = 6,67428 10^-11 м³с^-2кг-¹ или Н м² кг^-2. На этом примере мы видим как со временем менялась и уточнялась математическая модель (формула), описывающая реально существующий физический процесс, так как совершенствовались средства измерения при проведении экспериментов. Соответственно менялось наше представление об этом физическом процессе, сам же процесс все время оставался одним и тем же и не зависит от точности его описания.

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

 

max|ε₀ | ≤ ε_доп,

 

где ε_ij и ε_доп — допущенная и предельно допустимая относительные погрешностимоделирования поверхности,максимумберетсяповсемкоординатам и контролируемым точкам;

ε_ij= (x _ijист – x_ijмод)/ x _ijист

x_ijист, x_ijмод — i -я координата j -той точки поверхности в объекте и моделисоответственно.

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

С областью адекватности тесно связана еще одна характеристика МО – чувствительность к изменению данных. При этом различают устойчивые и неустойчивые МО. Для устойчивых моделей можно гарантировать, что погрешность будет находиться в некоторых заданных пределах, если изменения параметров МО тоже будут происходить в ограниченном диапазоне значений. Математически устойчивость описывается как отношение изменения результата Df к изменению конкретного параметра Dx:

Для нескольких параметров x₁, x₂, … x_n устойчивость будет определяться как сумма устойчивостей по каждому параметру:

|¶f /¶x₁| + |¶f /¶x₂| + |¶f /¶x₃| +... + |¶f /¶x_n|.

Для устойчивых МО это отношение меньше 1, для неустойчивых – больше 1. Примером устойчивой МО в механике является игрушка-неваляшка, а неустойчивой – спичечный коробок, поставленный на ребро. В первом случае небольшое изменение положения центра тяжести вызовет возврат к положению покоя, а во втором к нарастанию этого смещения. В радиотехнике устойчивость работы усилительных каскадов приемника обеспечивается отрицательной обратной связью, в противном случае в этих каскадах возникает возбуждение (резкий свист). Наоборот для генератора нужно, чтобы состояние было неустойчивым. И электронная схема переключалась поочередно из одного поллярного состояния в другое.

Рассмотрим еще один пример МО, описываемой уравнением

x³ – 3x² + 3x – a = 0

где a – некоторый параметр.

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

Для этого несколько преобразуем уравнение к виду:

 

(x-1)³ + (1-a) = 0

Теперь определим корень при а=1. Несложно видеть, что x_k = 1.

Изменим a на 0,001 (0,1% от первоначального значения) до величины a = 1,001. Теперь корень уравнения также изменится и будет вычислен как x_k = 1,1. (изменится на 10% от первоначального значения). Налицо неустойчивая МО, так как ¶x_k /¶a = 10% / 0,01% = 1000.

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

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

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

Пример. Найти решение системы уравнений:

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

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