Математические модели при автоматизированном проектировании технологических процессов

 

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

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

Мv= {Tи, m, t, s, d, B, cv, kv, xv, yv, zv, rv}, (17.1)

где Ти - стойкость инструмента, мин; m - показатель относительной стойкости инструмента; t - глубина резания, мм; s - подача, мм/об (мм/зуб мм/дв.ход, мм/мин); d - диаметр обрабатываемой поверхности или диаметр инструмента; В - ширина обрабатываемой поверхности, мм; сv - коэффициент, характеризующий условия обработки; kv - поправочный коэффициент на скорость резания; хv, yv, zv,rv- показатели степени.

Логическое соотношение между приведенными выше параметрами и скоростью резания v имеют вид

v = Ти Ù m Ù сv Ù kv [(t Ù xv) Ú (s Ù yv) Ú (d Ù z) Ú (B Ù rv)] (17.2)

причем Ти , m, сv, kv всегда истинны, а истинные значения других переменных зависят от метода обработки резанием.

Формулы количественных соотношений между параметрами с учетом истинности их логических значений имеют вид:

при наружном точении

(17.3)

 

при сверлении -

(17.4)

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

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

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

у(х)= ВТf (х) + x(х), х, х Îc, (17.5)

где х = (х1, х2, ..., х1)Т- вектор входных (контролируемых) переменных;

В = (b1, b2,..., bm)Т - вектор неизвестных параметров;

f(х) = (f1(х), f2(x),...,fm(x))Т - вектор известных (регрессионных) функций;

x (х) - помеха (например, ошибка измерения);

c - область измерения входных переменных х, в которой справедлива модель.

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

на дискретном множестве точек {x1} ( ) по ri раз в точке xi.

Результаты измерений при этом могут быть записаны в виде

yik = yk (xi ) = ВТf (хi ) + xk (xi ) = ВТ f (хi) + xik , (17.6)

xi Î X,

где yik и xik - соответственно значения величины y(x) и помехи x(х) при k-м изменении в точке хi ( ).

В теории ПЭ совокупность величин

y11, y12, ..., yir1; ...; yn1, yn2, ..., ynrn

r1; ...; rn (17.7)

x1; ..., xn

обычно называют экспериментом ЭN. Здесь N - общее число изменений

( ) совокупность величин

p1, p2, ... , pn, , (17.8)

x1, x2, ... , xn

называют планом eN эксперимента ЭN. (Иногда план вида (17.8) называют нормированным планом эксперимента.) Совокупность точек {x1} ( ) называется спектром плана eN. И, наконец, часть области Х, в которой могут размещаться измерения, называется областью планирования Х. В частном случае область Х может совпадать с c. Предположим, что целью проведения эксперимента является нахождение наилучших в статическом смысле (т.е. имеющих минимальные дисперсии и несмещенных) линейных оценок неизвестных параметров Впо результатам измерений выходной величины (по данным эксперимента). Статистическая точность оценок (т.е. статическая близость их к истинным значениям, обычно и определяющая понятие «наилучшие оценки»), которую можно получить в результате проведения эксперимента ЭN по плану eN , описывается их дисперсионной (ковариационной) матрицей

, (17.9)

где - вектор несмещенных оценок параметров B.

Задача проведения и анализа результатов эксперимента естественным образом распадается на две самостоятельные задачи.

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

Б. Выбор оптимального в некотором смысле плана эксперимента eN при заданном способе обработки его результатов. Предметом теории ПЭ в основном является задача Б. Выбор плана эксперимента, соответствующего наилучшей дисперсионной матрице, предусматривает сравнение дисперсионных матриц различных планов. Обычно при этом сравниваются величины некоторого выпуклого функционала (скалярные меры) от дисперсионной матрицы оценок y[D(.)]. Вид функционала y[.] определяется выбранным критерием оптимальности плана. В настоящее время в литературе существует большое число различных критериев оптимальности планов и соответствующих этим критериям функционалов. Нами в дальнейшем будут использоваться следующие функционалы y[D] (здесь и ниже аргумент для краткости опущен).

1. Определитель дисперсионной матрицы: y[D] = |D|.

Планы, минимизирующие этот определитель, принято называть D - оптимальными. При этом мимизируется обобщённая дисперсия (эллипсоид рассеяния оценок).

2. Сумма диагональных элементов дисперсионной матрицы:

y[D] = Sp[D]. Планы минимизирующие Sp[D], называются А- оптимальными. При этом минимизируется средняя дисперсия оценок.

3. Дисперсия оценки линейной комбинации qТВ параметров В: y[D] = qТDq, где q = (q1, q2, ..., qm)T - заданный вектор. Важным частным случаем этого функционала является случай, когда q = F(x0). При этом величина d(x0) =f T(x0) Df (x0) есть дисперсия оценки так называемой поверхности отклика ВТ f(x) в точке х0. В общем случае точка х0Îc может не принадлежать области ХÌc.

4. Максимальное в области Х0 значение дисперсии оценки поверхности отклика: sup d(x). Планы, минимизирующие указанную величину, xÎc0 называются минимаксными. Если Х0 совпадает с Х, то соответствующий оптимальный план принято называть G - оптимальным.

Обычно различают три случая соотношения областей Х0 и Х:

область Х0 совпадает с областью Х;

область Х0 является частью области Х;

область Х0 не имеет общих точек с областью Х (задача экстраполяции) Но во всех трех случаях должно выполнятся Х0 Ìc и Х Ì c (c - область, в которой справедлива модель (17.5). При традиционной постановке задачи ПЭ считают, что помехи xik в (17.6) удовлетворяют следующим предположениям:

Е xik (Хi ) = 0; Е xk (Хi ) xs (Xj) = s2 (хi) при i=j и k=s, (17.10)

0 в остальных случаях;

Функцию l(х), обратную функции s2(х); где l(х)= s -2(х), принято называть функцией эффективности. Таким образом, из (17.10) видно, что в традиционном ПЭ результаты всех измерений уi считается статистически независимыми (точнее говоря, некоррелированными) независимо от того, есть ли среди точек хi совпадающие или нет.

Решением задачи А для модели (17.6), (17.10) являются оценки метода наименьших квадратов(МНК) Во (ЭN) (МНК-оценки). Точность МНК- оценок характеризуется дисперсионной матрицей

D(Bo(ЭN))= D (Э0N) = N-1D(eN),

где D(eN)= [S pi l(xi) f(xi) f Т (хi )]-1 (17.11)

Выражение (17.11) показывает, что элементы матрицы D(Э0N) зависят от расположения точек хi и частоты повторения в них измерений рi (т.е. от плана эксперимента) и не зависит от результатов измерений выходной величины в этих точках. Это означает, что, во-первых, здесь можно поставить задачу выбора наилучшего (оптимального) в некотором смысле плана эксперимента и, во-вторых, план эксперимента можно выбирать априори.

Обычно под оптимальным планом в ПЭ понимается план e*N1 для которого выполняем условие

y[D(e*N)] = inf y[D(eN)] (17.12)

где минимум (нижняя грань) берется по множеству всех допустимых планов. В общем случае поиск плана e*N, удовлетворяющего (17.12) представляет громоздкую экстремальную задачу. На практике различают три более простых случая.

а) Поиск оптимального плана при заданном числе измерений N. Область поиска оптимального плана ограничивается множеством всех планов с числом измерений, равным N.

б) Поиск оптимального плана при заданной точности определения оценок исходных параметров. ПЭ в этом случае заключается в поиске минимума y[D(eN)] по плану ex и числу измерений N при заданном y0. Практически задача сводится к последовательному поиску минимумов (17.12), который продолжается до тех пор, пока не будет найден план eN+1 для которого

y[D(eN)] > y0 > y[D(eN+1)] (17.13)

При этом эксперимент проводится по плану eN+1.

в) Наиболее сильные результаты удаётся получить, переходя к понятию непрерывного плана e путём устремления N к бесконечности. Множества дискретных планов тогда будут подмножествами непрерывных планов e соответствующими таким pi, для которых ri= piN - целые числа. Одно из важнейших свойств непрерывного плана, непосредственно следующие из его определения, заключается в том, что величина inf y[D(e)] не зависит от N. Поэтому планирование в этом случае сводится к отысканию плана e*, минимизирующего y[D] в классе всех допустимых планов e.

y[D(e*)] = y[D(e)] (17.14)

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

М(e)= òх М(х) р(х) d(x) = å pi М(хi), (17.15)

где М(х) = l(х) f(x) fТ(х) - информационная матрица плана, сосредоточенного в единственной точке х (здесь n < m(m+1)/ 2+1). Это свойство дает возможность рассматривать только планы с дискретным спектром.

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

 


Глава 18






Дата добавления: 2017-11-21; просмотров: 867; ЗАКАЗАТЬ НАПИСАНИЕ РАБОТЫ


Поиск по сайту:

Воспользовавшись поиском можно найти нужную информацию на сайте.

Поделитесь с друзьями:

Считаете данную информацию полезной, тогда расскажите друзьям в соц. сетях.
Poznayka.org - Познайка.Орг - 2016-2022 год. Материал предоставляется для ознакомительных и учебных целей.
Генерация страницы за: 0.044 сек.