Инженерная эконометрика в задачах системного анализа Текст научной статьи по специальности «Математика»

УДК 681.5:330.43

ИНЖЕНЕРНАЯ ЭКОНОМЕТРИКА В ЗАДАЧАХ СИСТЕМНОГО АНАЛИЗА

В.А. Кочегуров, Л.И. Константинова, В.Г. Гальченко

Томский политехнический универсистет E-mail: gal@am.tpu.ru

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

Ключевые слова:

Эконометрика, система, математическая модель, обобщенные координаты, обобщенные векторы Лагранжа, системные единицы.

Key words:

Econometrics, system, mathematical model, generalized coordinates, generalized vectors of Lagrange, system unit.

Введение

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

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

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

Объективно существующие системы

и способы их представления

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

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

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

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

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

Исследование систем

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

Рисунок. Стратифицированное отображение системы

Модель поведения системы «М» реального объекта «О» является тогда, когда выполняются условия:

• реальный объект воспринимается так, что в нем выделяется существенное и несущественное;

• существенным величинам объекта «О» сопоставлены соответствующие величины объекта «М»;

• определены соотношения между существенными величинами конкретного объекта «О» и соответствующие соотношениям между образцами этих величин объекта «М».

Модели поведения должны адекватно отражать поведение реальных объектов в смысле соответствия между причинами и следствиями. Имеются классические задачи исследования систем:

• моделирование (анализ);

• прогнозирование;

• идентификация;

• управление.

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

Введение основных единиц измерения

системных величин

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

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

В качестве примера рассмотрим уравнение состояния автономных систем, широко используемое в теории управления [4]:

х (?) = р (X (0) . (1)

Здесь ДО - переменные состояния системы, единицы измерения:

? = [?я ], X = [VI ], X = Щ.

•• дР • дР

Из уравнения X(?) = ~^ X(?), где д— - матрица Якоби, введем меру инерционного замедления: Мг = [я1г, я2г,...,ят ]т

и соответствующие силы, определяющие покоординатные движения на интервале &:

Q(t ) = MrX (t ),

где

т1г 0 ... 0

Мг = 0 т2г 0

0 0 тпг_

и m = — 1 m 1 — . т_

dfT f

dxdx

dfT f

dxdx

dxdx

Введем единицу измерения замедления - inde-cel [idc]. Тогда обобщенная сила Q(t) будет иметь

. [idc][F?]

единицу измерения: Q(t) = —[—-2—, которую на-

[tm ]

зовем tensor [tns].

Количество движения вдоль каждой координаты (импульс) определиться:

Pi(t ) = Щг xi(t X

P2(t) = m2,x2{t),..., pn{t) = mmxn(t).

Тогда системная сила Q равна

Q(t ) =

dP(t ) dt

= CP ( X (t )),

где

C =

dfi df2 т1г df1m1r

dx1 dx2m2z dxnmrn

df2 т2г df2 df2m7.

dx1m1z dx2 dxnmnr n m

dfnmrn dfn

dxnm1T dxn

Матрица С характеризует упругие свойства си-

^ 1

стемы, ее единица измерения: С =-------.

[?я ]

Работа, совершаемая системной силой, равна: Ж =| дт (?)сК (?) = | (Ыг X (? )тХ (? )Ш? =

í (dX)TMTrX (t) = 2 X (t )T MTX (t ).

(2)

Здесь единицы измерения энергии:

-|2

V?

W = [idc]

VI

t„

Из приведенных системных единиц могут быть получены любые производные.

Системные переменные состояния

и обобщенные векторы Лагранжа

В механике движения для описания процессов широко используются обобщенные векторы Лагранжа д{дьд2,...,д„}и #{#ь#2,...,4,}. Для них введено понятие функции Лагранжа ЬЦдд), определяющей энергетическое состояние объекта как разность между кинетической энергией Жк(д,д), зависящей в общем случае от д и д,и потенциальной энергией Щ,(д), зависящей только от д:

Щ, д) = Ж ^, д) - Ж„ (д).

Частная производная кинетической энергии по вектору скорости называется обобщенным импульсом:

Рд ч) = дЖ (ч ч) =дь(ч ч).

дд дч

Частная производная потенциальной энергии Щ,(д) по обобщенному вектору д называется обобщенной силой:

дЖп (д)

Ô(q) =-

dq

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

Переход к обобщенным координатам

Примем, что новые координаты д не зависят от производных старых координат и переход к новым координатам определяется нелинейным преобразованием Д0=Дд,д).

Отсюда X(?) = дХ_ ^ = X(д, д).

дд ш

Тогда кинетическая и потенциальная энергия определяются как:

Ж (X(д, д)) = Жк (д, д), Шп (X(д)) = Шп (д).

Вычислим

дЖк (q, д) „ дЖк ^ д)

дд

дд

дЖк (X(д, д) дЖкт ¿X

дд дXт ¿Жк дд ¿X

дX ¿д

Ш дX^ Ж дд

Р( X),

дЖк (X(д, д))

дд

д_

дд

дд ¿X дд

Р( X) =

дX (?) дд

дXi

Обозначим Ои

дЖп

Р (X) =—Р (X).

дд

дXT ШР(д, д)

= Он + Оа ■

дд дц1 Ш

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

ш дЖк д д) дЖк (д, д)

Ж дд1 дд1

Ш дЬ(д, д) дЬ(д, д)

■ = Ож

= О,и

г = 1,..., п.

ддг %

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

ШР1 = дН Шдг = дН Ш ддг А др ’

где И=Жк+ЖП - полная энергия системы.

Равновесные сбалансированные модели

функционирования динамических систем

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

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

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

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

В соответствии с (1) поведение динамической системы, или характер изменения во времени переменных состояния X в пространстве состояний за счет внутренних и внешних источников энергии определяется некоторой силовой функцией ДХ(0).

Обозначим функциональное состояние системы при ее нормальном функционировании ^(Х(/))=^(Х0(/)). Тогда при малых отклонениях значений Х(0 от Х0(О уравнение (1) можно переписать в виде:

XX ) = ддХ± • X (?).

дх

(3)

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

хг + А?) - хг )

Кп =:

А?х (?)

и

При устремлении At к нулю имеем хі (t + At) X (t) Atxi (t) x (t)

limAt^0 KTi = limAt0

= К (t).

Отсюда х(1)=Ктх() или Х(0=КтХ(0 где Кт - диагональная матрица с элементами Кт.

Подставляя в уравнение (3):

XX ) = ддх) (?).

дх

Введем обобще.нные переменные состояния Р(0=КтХ(0 и О( 0=Х(0, их можно еще назвать фазовыми переменными. Тогда п дифференциальных уравнений второго порядка можно записать в виде 2п дифференциальных уравнений первого порядка:

ШР

■ = КтО(?),

dt dQ(t)

dF (x)

P(t).

Ш дх

В случае устойчивых систем эти уравнения описывают системы, состоящие из п взаимосвязанных гармонических осцилляторов. Собственные числа

матриц Кт и др(х1

дх

определяют частотные свойства

систем и являются мерой, определяющей темпы изменения переменных состояния.

Можно ввести скалярную величину, характеризующую энергию системы в пространстве обобщенных координат:

дЕ (х)

P (t)-

дх , Q (t)KTQ

W(P(t), Q(t) = - 2 2

Тогда уравнения (2) можно записать dP(t) dW (P (t), Q (t)

dt dQ(t)

dQ

dW (P (t), Q (t)

- = KTQ (t)

dF (x)

P(t).

Ш дР дх

Эти уравнения можно объединить, если ввести переменную

7 (?) = {Р(?), О(?)} =

= {2Х(?),2г($),..., 2п(?), 7п+1(?),..., 7п(?)};

7 (?) = Ы21 (?)

где

Ы2 =

0

dF (х)

дх

К

0

При кососимметричной матрице Mz описываемая система является самосопряженной, и для нее выполняется условие

(ZT (t )Z (t)) = | |Z (t )||2 = const.

По определению квадрат нормы характеризует энергию системы.

Условие кососимметричной матрицы Mz выпол-dF (x)

няется при Кт = -

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

Рассмотрим одиночный гармонический осциллятор. Для него

х (?) + дМ. х (?) = 0,

дх

х(?) = ктх(?).

Введем обобщенные переменные р(0=ктх(/), д(0=х(0.

Тогда

Ш.=1гМ1), ш.=-дхР(,У

ш ш дх

Обозначим 11( 0=р(0, г2(0=д( 0.

Тогда

¿(Г) = Mzz(t),

где Mz =

При kT =

0

df (х) дх df (х) дх

система становится самосопря-

женнои и

(zT (t) z (t) = zf(t) + z^(t) = p 2(t) + q 2(t) = const = c. Вернемся к исходным переменным:

(4)

, 2 2Ґ,\ , -2/.\ 2 х (t) , х і k, х (t) + х (t) = c , —— + — = I,

a

где a = -

< = c.

Уравнение (4) является уравнением эллипса с полуосями а и b.

Площадь эллипса равна

с2

Sm = nab = n—.

kT

Она остается неизменной со временем, если взаимосвязь между x(t) и Х(t) не нарушается, и является энергетической характеристикой осциллятора. Отсюда

k =—— ЯТ 2

ПС

= const.

дх

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

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

T

Выводы

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

процессов, не вдаваясь в их конкретное физическое содержание.

На основе понятия «гомеостаз» получены условия и критерии оценки сбалансированного равновесного динамического функционирования систем.

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

1. Лопатников Л. Экономико-математический словарь. - М.: Изд-во ДВУ, 1996. - 700 с.

2. Перегудов Ф.И., Тарасенко Ф.П. Введение в системный анализ. - М.: Высшая школа, 1989. - 367 с.

3. Системный анализ и принятие решений: словарь-справочник / под общ. ред. В.Н. Волковой, В.Н. Козловой. - М.: Высшая школа, 2004. - 612 с.

4. Сю Д., Майер А. Современная теория автоматического управления и ее применение. - М.: Машиностроение, 1972. -530 с.

5. Николис Дж. Динамика иерархических систем. Эволюционное представление.- М.: Мир, 1989. - 474 с.

Поступила 25.03.2011 г.

УДК 621.01:519.242

МЕТОДИКА ОПТИМИЗАЦИИ ФОРМЫ ЭЛЕМЕНТОВ МЕХАНИЧЕСКИХ СИСТЕМ

А.В. Янушевскис, А.Г. Мельников, И.Д. Кактабулис

Рижский технический университет, Латвия E-mail: janush@latnet.lv

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

Ключевые слова:

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

Key words:

Shape optimization, non-uniform rational B-splines, design of experiments, metamodel, tensometric wheel pair of wagon.

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

Р1 ^0;шзх^0

Г V- ’

г,шах г

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

Метод особенно успешно применяется [3 и др.] для оптимизации оболочечных конструкций. Однако его применение может быть затруднительным, так как количество параметров оптимизации может быть огромным (сотни тысяч и даже более миллиона) и для решения задачи требуется мощное программное обеспечение и вычислительные ресурсы. В случае оптимизации твёрдых тел полученные оптимальные формы часто трудно реализуемы или даже принципиально нереализуемы технологически. Поэтому актуальна разработка альтернативных подходов.

Предлагаемая альтернативная методика оптимизации формы элементов механических систем предусматривает применение высококачественных метамоделей и использование на соответствующих этапах программное обеспечение геометрического моделирования и метода конечных элементов (МКЭ), а также программ для оптимизации и аппроксимации. Разработанная методика основывается на последовательном использовании CAD/CAE