Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
u_course_instrument_making.pdf
Скачиваний:
93
Добавлен:
12.02.2015
Размер:
3.48 Mб
Скачать

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

Окончание табл. 2.2

1

2

3

4

5

6

Консерватив-

 

 

 

 

 

ный

 

Индуктив-

 

Коэффици-

Тепловая

параметр I ро-

Емкость

Масса

ент податли-

ность

емкость

да

 

 

вости

 

 

 

 

( -параметр)

 

 

 

 

 

 

 

 

 

 

 

Консерватив-

 

 

 

 

 

ный

Индуктив-

 

Коэффициент

 

Тепловая

параметр

Емкость

Масса

индуктив-

ность

податливости

II рода

 

 

 

 

ность

( -параметр)

 

 

 

 

 

2.2. Функции параметрической чувствительности

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

Чувствительностью к изменению внутренних параметров обладают практически все выходные характеристики технического процесса. В общей теории чувствительности систем [11; 16] широко используется количественный показатель параметрической чувствительности в виде частной производной выходной характеристики yj по соответствующему, например, k-му параметру q, который называют функцией параметрической чувствительности (ФПЧ):

Ay j

( ) =

yj

 

 

 

 

,

(2.13)

 

 

 

 

Q

 

 

qk

 

q

k

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Q0 – вектор расчетных значений параметров, при которых вычисляется

частная производная.

Знак численного значения ФПЧ определяет направление изменения yj-й характеристики при изменении qk-го параметра в сторону увеличения.

Смысл введения количественных показателей ФПЧ выходных характеристик схемы, конструкции или технологического процесса ПС в виде функции (2.13) состоит в том, что в инженерной практике в большинстве случаев при

малых вариациях параметров Q

2

 

k 2

Ji y j

3 ,

 

 

Lj

Компьютерные технологии в приборостроении. Учеб. пособие

-56-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

Q = q1, q2 , ..., qL T ,

вариации Y = y1, с y2, ..., yN T выходных характеристик Y точностью можно выразить через абсолютные ФПЧ (2.13):

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

Y

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ay1

Ay1

...

Ay1

 

 

 

 

 

 

 

 

 

 

 

 

q1

 

q2

 

 

 

qL

 

 

 

 

 

 

 

 

 

 

 

 

Ay2

Ay2

...

Ay2

 

 

 

 

 

 

 

Y

 

 

где

A

 

 

 

q1

q2

 

 

 

qL

.

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

... ... ... ...

 

 

 

 

 

 

 

 

 

 

 

AyN

AyN

...

AyN

 

 

 

 

 

 

 

 

 

 

 

 

q1

q2

 

 

 

qL

 

 

 

 

 

абсолютная

матрица

 

чувствительности,

причем

Y

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с приемлемой

(2.14)

(2.15)

Q Q Q0 ,

а Y Y Y 0 .

При решении многих задач проектирования ПС, когда по результатам расчетов необходимо количественно сравнивать ФПЧ, применяют относи-

тельную форму представления ФПЧ:

 

 

 

 

S y j ( )

= Ay j ( )

qk

.

(2.16)

 

q

q

y j

( )

 

k

k

 

 

 

 

 

При этом для проектных процедур «ПРИНЯТИЕ РЕШЕНИЯ» в методиках автоматизированного проектирования используют матрицу относительных ФПЧ:

 

 

 

 

y

y

...

y

 

 

 

 

 

Sq 1

Sq21

SqL1

 

 

 

 

 

 

1

 

 

 

 

 

SY

 

Sqy22

 

SqyL2

 

Sqy12

...

.

(2.17)

Q

 

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yN

yN

...

yN

 

 

 

 

Sq1

Sq2

SqL

 

 

Информация о функциях чувствительности позволяет выполнить операции настройки ПС. Так, если перейти от дифференциала к приращению, то можно записать:

S

y

j

y j

 

q

k

или

y j

S

y

j

q

k .

 

qk

 

qk

 

y j

 

y j

 

qk

qk

 

Компьютерные технологии в приборостроении. Учеб. пособие

-57-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

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

и при этом SRU4вых = – 1,25, то увеличение, допустим, R4 на 1 % ( R4 R4 = 0,01)

приведет к уменьшению Uвых на 1,25 %, т. е. Uвых 1, 25 0,01 0,0125 . Если также принять Uвых = 100 В, то Uвых 1, 25 100 0,0125 ,

т. е. можно ожидать, что Uвых изменяется от 100 до 98,75 В.

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

2.2.1. Методаналитическогодифференцирования

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

ных характеристик Y к изменению внутренних параметров Q путем непо-

средственного дифференцирования явных (y = f(x, q, ) или неявных (F(x, y, q, ) = 0) функциональных зависимостей, выступающих как наиболее общая форма аналитических расчетных моделей. Рассмотрим такой подход на простых примерах.

Пример 1. Вынужденное перемещение y = амортизированного блока ПС при движении объекта установки блока с линейным ускорением x = a описывается выражением

 

 

a

 

2

1

2

 

 

 

или в унифицированных обозначениях:

y =

 

x

,

(2.18)

q2

1+ q2

1

2

 

 

где q1 = – собственная частота колебаний системы «блок – амортизатор»; q2 = – коэффициент демпфирования.

Допустим, надо определить S или в унифицированных обозначениях Sqy1 . Для этого необходимо использовать правила вычисления производных

(табл. 2.3) и частных производных, а также выражение (2.16). Учитывая это, получим

Компьютерные технологии в приборостроении. Учеб. пособие

-58-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

Sqy

=

y

 

q

=

y

 

q3

1+ q2

=

 

2x

 

q3

1+ q2

 

 

1

 

 

1

2

 

 

 

1

2

.

q1

q1

 

x

q3

1+ q2

 

 

1

 

 

y

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

В результате сокращений получим Sqy

2 .

 

 

 

 

1

 

 

Пример 2. Резонансная частота y = 0

параллельного резонансного

контура (рис. 2.9) определяется, как

 

 

1

, или в унифицированном ви-

 

 

 

 

 

 

0

 

LC

 

 

1

 

 

 

 

 

де, как y =

 

.

 

 

 

 

q q

 

 

 

 

 

1

2

 

 

 

 

 

C

L

R

Рис. 2.9. Схема параллельного резонансного контура

Для вычисления SLω0 SCω0

по аналогии с примером 1 осуществляем не-

обходимые операции:

 

 

 

 

 

 

 

 

 

 

 

Sqy

=

y

q1 =

0

 

L LC

=

1

;

Sqy

= Sqy

=

1 .

q

1

2

1

 

y

L

 

 

 

2

1

 

2

 

1

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.3 Некоторые правила вычисления производных

Функция

Производная

 

 

1

2

C (const)

0

 

 

X

1

 

 

Компьютерные технологии в приборостроении. Учеб. пособие

-59-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

 

X n

 

 

 

 

nX n 1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

X 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

X n

 

 

 

X n 1

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n x

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

n

n

 

x

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ex

 

 

 

 

 

 

 

 

 

 

ex

 

 

 

 

 

 

 

ax

 

 

 

 

 

axlnx

 

 

 

 

 

lnx

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

loga x

 

 

1 logae

 

 

 

1

 

 

xlna

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

lgx

 

 

1 lge

 

0,4343

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

x

sin x

 

 

 

 

 

 

 

 

cos x

 

 

 

 

 

 

cos x

 

 

 

 

 

sin x

 

 

 

 

 

 

 

 

 

 

 

 

Окончание табл. 2.3

1

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

tg x

 

 

1

 

 

 

 

 

 

sec2x

 

 

 

cos2x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ctg x

 

1

 

 

 

 

 

cosec2x

sin2x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arcsin x

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arccos x

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arctg x

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arcctg x

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1 x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ch x

 

 

 

 

 

 

 

 

sh x

 

 

 

 

 

 

 

th x

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ch2x

 

 

 

 

 

 

Компьютерные технологии в приборостроении. Учеб. пособие

-60-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

1 2

 

 

1

2

 

 

 

 

 

 

1 2

 

1 2 1 2

 

 

 

 

 

 

 

 

 

2

 

 

1 2

1

 

1

2

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

( (t))

 

 

 

 

 

 

( (t)) (t)

Для автоматизации проектных работ, прежде всего задач синтеза, необходимо унифицировать представления исходных расчетных моделей. Это означает, что аналитические расчетные модели должны быть записаны не в общепринятых обозначениях физических величин, разнообразных для электрических, электромагнитных, механических и тепловых процессов, а в унифицированных обозначениях – x, y, q, . Тогда один раз разработанный программный модуль аналитического дифференцирования может быть встроен в любое программное обеспечение или алгоритм расчета системы автоматизированного проектирования. При этом унификация аналитических моделей в программном модуле заключается в алгоритмизации перехода моделей, представленных в матричном виде ([C]. [Y]=[X]; [C] – матрица параметров) или в виде системы дифференциально-алгебраических уравнений (составленных, например, по методу переменных состояний (yj = Фj (x, y, q, t), где j – первая производная по времени) [16; 21], к моделям в форме y = (x, q, ) и их исследовании методом аналитического дифференцирования.

2.2.2.Таблицачувствительности

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

ренним параметром qk. Для этого проводят дифференцирование по параметру, умножение на qk0/y0 и упрощение полученных выражений. Пример полученных таким образом выражений приведен в табл. 2.4. Для повышения практической ценности таблицы выходная характеристика y рассмотрена в

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

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

Компьютерные технологии в приборостроении. Учеб. пособие

-61-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

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

 

 

 

 

 

 

 

Таблица чувствительности

 

 

 

 

Таблица 2.4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналитическая модель

 

 

Правило вычисления

 

п/п

 

 

 

относительной ФПЧ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общие случаи

 

 

 

 

 

 

 

 

1

 

 

 

y f (q )

 

 

 

 

 

 

 

 

 

S y 0

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

qk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Продолжение табл. 2.4

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

y f1(h1); h1 f2(h2), ...

 

 

Sqy Shy Shh1 ...

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

1

 

 

2

 

 

..., h

 

f

 

(h ), h

f

 

(q )

 

 

...Shn 1 Shn

 

 

 

 

n

 

 

 

 

 

 

n 1

 

 

n

n

 

 

 

n 1 k

 

 

 

 

hn

 

 

qk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y f h1 qk ;

h2 qk ;

 

 

 

 

hn qk

 

 

 

 

n

Shy

h

 

 

3

 

...;

 

 

Sqy

Sqi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

i 1

i

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

n

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

Sqyk hi Sqhki

hi

 

 

4

 

 

y h (q )

 

 

 

 

i

1

 

 

 

i 1

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1 i k

 

 

 

 

 

Sqhki

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hi

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

1

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

5

 

 

y hi

(qk )

 

 

 

 

Sqyk Sqhki

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

i

1

 

 

 

6

 

 

y h1(qk ) h2(qk )

 

 

 

Sqy

Sqh1

Sqh2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

k

 

k

 

 

7

 

y f (q ); h f

2

(q )

 

 

S y

S y

Sh

 

 

 

 

 

1

 

 

k

 

 

 

 

k

 

 

 

h

qk

 

qk

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

h

 

 

y f1 h(qk ) ; h

f2 y(qk )

 

 

 

Sh 1 Sy

 

Компьютерные технологии в приборостроении. Учеб. пособие

-62-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

9

10

11

12

13

14

15

1

16

17

18

19

20

F(y, y ) 0

 

 

 

Sqy

 

 

SqF

 

SyF

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

k

 

 

 

 

 

Частные случаи

 

 

 

 

 

S y

Sh

 

1

 

y h(q ) q

 

 

 

 

 

 

 

 

 

k

k

 

 

 

 

 

 

 

 

q

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

k

 

 

 

y c h(q ), c f (q )

 

 

 

 

 

 

 

 

S y Sh

 

 

 

k

k

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

Sqy

k

 

 

 

 

k

 

 

y hn (q ),

n const

 

 

 

 

 

 

n Sqh

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

k

1

 

 

 

 

k

 

y n h(q ),

n const

 

 

 

 

 

 

S y

n

Sh

 

 

k

 

 

 

 

 

 

 

 

 

qk

 

 

 

 

 

qk

 

y аh(qk ), а const

 

 

 

Sqy

 

 

h ln a Sqh

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

k

 

y logаh(qk ), а const

 

 

 

Sqyk

 

1ln h Sqqk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончание табл. 2.4

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

y h (q ) h2 (qk )

 

Sqy

k

h2Sqh1

h2 ln h1 Sqh2

 

1

k

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

k

 

y = sinh(q )

 

 

 

 

Sqy

 

h ctgh Sqh

 

 

 

k

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

k

 

y cosh(q )

 

 

 

Sqy

 

h tgh Sqh

 

 

 

k

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

y tgh(q )

 

 

 

S y

 

 

 

 

2h

 

 

 

Sh

 

 

 

 

 

 

sin2h

 

 

 

k

 

 

 

 

q k

 

 

 

 

 

q k

 

 

 

 

 

S y

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

Sh

 

 

 

 

 

 

1 h2 arcsin h

y аrcsin h(qk )

 

q k

 

 

qk

 

 

 

 

 

 

 

 

 

 

 

1

 

 

Sh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y ctg y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qk

 

 

Рассмотрим пример, иллюстрирующий применение метода таблицы чувствительности (табл. 2.4).

Пример. Используем функциональную зависимость из примера 1, при-

веденного в п. 2.2.1. Допустим, необходимо получить ФПЧ Sqy и

Sqy Для оп-

 

 

Sqy

 

1

2

ределения

, применяя п. 11 и 12 табл. 2.4, можем записать, что Sqy 2.

 

Sqy

1

 

 

1

ФПЧ

находится по табл. 2.4 следующим образом. Согласно п. 11 имеем

 

2

 

 

 

 

Sqy Sqh1 ,

где

h1 1

1 q22 . Далее, подставив h1 h2 1 2 , где h

1 q2 , по

2

2

 

 

2

2

 

 

 

 

Компьютерные технологии в приборостроении. Учеб. пособие

-63-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

 

Sqy

 

1

h

 

 

1

 

h3Sqh3

п. 12 получим

 

2

Sq

2

. Затем применим п. 4:

Sqy

 

 

2

, где

2

 

 

2

 

 

2

 

2

 

1 h3

 

 

 

 

 

 

 

 

 

h q2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sqh3

 

 

Но

 

при

 

этом,

 

согласно

п.

12,

2 .

Следовательно,

 

 

 

 

 

q2

2

 

 

 

 

 

 

 

 

2

 

 

y

 

1

 

 

 

 

 

q2

 

 

 

 

 

 

 

Sq

 

 

 

 

2

 

 

 

 

2

 

.

 

 

 

 

 

2

1

 

2

1

 

2

 

 

 

 

 

2

 

 

q2

 

q2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.2.3. Методприращений

Наиболее простым методом получения ФПЧ является метод приращений (метод вариаций) [8; 16], который применим к любым видам моделей и сводится к поочередной вариации (малому изменению, обычно в практических расчетах изменение принимается 1–5 % от значения параметра) каждого параметра q и регистрации соответствующих вариаций выходных характеристик y. Так, абсолютная ФПЧ рассчитывается приближенно следующим образом:

y j

 

y j

 

y j q1, q2 , ..., qk qk , ..., qi

 

y j q1, q2 , ..., qk , ..., qi

, (2.19)

qk

qk

qk qk qk

qk qk qk

а относительная ФПЧ:

y j

y j

qk

 

.

(2.20)

Sqk

Aqk

 

 

y j q1 , ..., qN ,...,

qL

Рассмотрим, например, расчет вектора относительных ФПЧ (строка матрицы (2.14) резонансной частоты f0 печатного узла (крепление в 4 точках) к изменению его геометрических и физико-механических параметров. Резонансная частота для такого способа крепления печатного узла рассчитывается следующим образом [11; 25]:

f

 

 

 

 

 

а2

 

Eh3 аb

 

(2.21)

 

 

1

 

 

 

,

0

2а2

b2

12 1 2 М

 

 

 

 

 

 

 

где а и b – длина и ширина печатной платы; Е – модуль Юнга;

h – толщина печатной платы; – коэффициент Пуассона; М – масса печатной платы с радиоэлементами.

Тогда, например, Sаf0 и Shf0 будут рассчитываться следующим образом:

Компьютерные технологии в приборостроении. Учеб. пособие

-64-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

f

 

 

f0 а а,b,h, E, ,M f0 а,b,h, E, ,M

 

 

 

а

 

Sа

0

 

 

 

 

 

 

 

,

а

f0

а,b,h, E,μ,M

 

 

 

 

 

 

f

 

 

f0 а,b,h h, E, ,M f0 а,b,h, E, ,M

 

 

 

h

 

Sh

0

 

 

 

 

 

 

 

.

h

 

f0

а,b,h, E, ,M

 

 

 

 

 

 

Таким образом, вычисляя остальные значения ФПЧ, можно получить следующий вектор относительных ФПЧ (строка матрицы S QY – выражение

(2.17):

SQf0 Sаf0 , Sbf0 , Shf0 , SEf0 , S f0 , SMf0 .

Сравнивая между собой составляющие полученного вектора, можно выбрать параметр, в наибольшей степени влияющий на резонансную частоту анализируемого печатного узла. Изменяя выбранный параметр, можно добиться, например, повышения резонансной частоты печатного узла при данном способе крепления. Можно также проранжировать полученные ФПЧ по степени влияния параметров (а, b, h, ..., М) на f0. В дальнейшем, выбирая варьируемые параметры из полученного ряда по конструктивнотехнологическим соображениям, можно изменять резонансную частоту проектируемого печатного узла.

Рассмотренный метод можно представить в виде алгоритма, приведенного на рис. 2.10.

Метод приращений имеет ряд существенных недостатков. Во-первых, ФПЧ определяются с ошибкой, так как функции y(q) в (2.16), как правило, нелинейные. Во-вторых, для вычисления ФПЧ yi по j параметрам требуется (j+1) раз рассчитывать модель, что связано со значительными временными затратами в случае применения моделей, детально описывающих исследуемые технические процессы. В практике автоматизированного проектирования, как правило, используются наиболее эффективные методы, такие как метод преобразованной модели [11], метод присоединенной цепи (смежной модели), метод транспонированной матрицы [11; 16]. Указанные методы получения ФПЧ хорошо отработаны для топологических моделей в унифицированной форме ненаправленных графов, позволяющих на основе единого программного обеспечения исследовать ФПЧ как электрических, так и тепловых и механических характеристик ПС.

В основе топологических методов получения ФПЧ лежат определенные принципы, например: принцип построения дополнительной модели и принцип автономного анализа двух моделей. На основе этих принципов построены такие методы получения ФПЧ, как метод преобразованной модели,

Компьютерные технологии в приборостроении. Учеб. пособие

-65-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

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

Начало

Решение математической модели

Цикл по исследуемым

А выходным характеристикам

0

yi , i 1, N

Цикл по варьируемым пара-

метрам qj ,

*

q j q j q j

Решение математической модели с использованием q* j

Вычисление

* 0

Aqyij *yi yi и Sqyij

q j qj

*

qj q j qj

Результат решения:

0

0

0

 

y1

, y2

,..., yN

вектор выходных характеристик

Вывод

S QY

Конец

А

Результат

*

решения: yi

Рис. 2.10. Блок-схема алгоритма вычисления матрицы S QY методом вариации параметров

Компьютерные технологии в приборостроении. Учеб. пособие

-66-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

2.2.4.ТопологическиеметодыполученияФПЧ

2.2.4.1.Методпреобразованноймодели

На рис. 2.11 приведена схема принципа построения дополнительной модели.

 

 

 

 

Исходная

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

модель

 

 

 

 

 

 

 

A

 

 

 

y1

 

y2

 

yN

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

 

,

,...,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qk

 

q

k

 

q

k

 

q

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дополни-

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тельная

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

модель

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 2.11. Принцип построения дополнительной модели

Принцип можно сформулировать следующим образом: любую ФПЧ можно представить как входную величину некоторой расчетной модели, дополнительно построенной и присоединенной к исходной модели так, что входные X, выходные Y и некоторые промежуточные U переменные исходной модели служат входными переменными дополнительной модели. Дополнительная модель может быть построена путем дифференцирования исходной модели по соответствующему параметру [11; 16]. В методе преобразованной модели (ПМ) это приводит к необходимости преобразования пассивных и активных зависимых компонентов (источников тока и напряжения, управляемых током и напряжением и т. п. [11; 16]), т. е. в дифференцировании исходной модели, представленной в топологическом виде (рис. 2.7), по интересующему параметру:

 

 

ijt hijt ( i

j ) hijt ij ,

 

 

(2.22)

 

ijt

kijt

 

;

ijt

kijt

 

mn

,

 

 

 

 

m

m

 

m

 

 

 

 

 

ijt

kijt

 

;

ijt

kijt

 

mn

,

(2.23)

 

m

m

 

 

m

 

 

 

 

где kimjt – коэффициенты пропорциональности.

В результате дифференцирования выражения показывают, какие преобразования в исходной модели надо выполнить, чтобы она превратилась в дополнительную, т. е. преобразованную модель (cм. табл. 2.5). Продифференцируем выражение (2.22) для двух случаев.

1-й случай. Параметр ветви hijt не зависит от первичного параметра qk, т. е. чувствительность к которому исследуется hijt (qk) (см. п. 1 табл. 2.5).

Компьютерные технологии в приборостроении. Учеб. пособие

-67-

ij

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

Aq ijt hijt

ij

hijt Aq ij .

(2.24)

qk

k

k

 

 

 

2-й случай. Параметр ветви hijt зависит от qk, т. е. hijt= (qk). В этом случае правая часть (2.22) дифференцируется как произведение двух функций,

зависящих от одного аргумента (см. п. 2 табл. 2.4):

ijt

 

ij

 

hijt

ij

 

 

 

Aqk

hijt

 

 

 

 

ij hijt Aqk

hijt

ij .

(2.25)

q

q

k

 

 

k

 

 

 

 

 

 

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

п

 

(2.26)

iju hijt ij ,

где ij относится к исходной модели, а знак «п» обозначает преобразованную модель.

Рассмотрим также два случая для нелинейной пассивной ветви в исходной модели.

1-й случай. Ветвь описывается зависимостью ijt= ( ij), которая не

содержит первичный параметр qk. Дифференцирование этого уравнения по параметру qk дает

A ijt f '(

)

ij

 

f ( ij )

A ij .

(2.27)

 

ij

qk

ij

 

qk

qk

 

Так как потенциальная переменная величина является функцией независимого аргумента (см. выражение (1.1) и результатом решения исходной аналитической модели, то, сравнивая выражения (2.27) и (2.24), получаем, что нелинейная ветвь исходной топологической модели переходит в линей-

ную ветвь с переменным параметром hijt f ( ij ) (см. табл. 2.5, п. 3).

Компьютерные технологии в приборостроении. Учеб. пособие

-68-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

Таблица 2.5 Правила перехода от исходной топологической модели к преобразованной

Ветвь исходной

Ветвь преобразованной

п/п

модели

модели

1

2

3

 

 

 

 

 

 

h

 

 

f(q

)

 

j

 

Aqki

 

 

 

 

 

h

 

 

 

 

 

 

 

Aqkj

1

 

 

 

i

 

 

 

 

ijt

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A ijt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Aqki

 

 

h

 

 

 

 

 

 

Aqkj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijuп

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

ij

i

 

j

 

j

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

ijt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

h'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iju

 

 

ijt

 

 

ij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

 

п

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

ijt

 

 

 

iju

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt f ( ij )

 

j

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A j

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Aq

 

 

 

 

 

hijt f ( ij )

 

 

qk

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

3

 

 

ijt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A ijt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A i

 

 

hijt fq ( ij ,qk )

 

 

A j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

i

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

f ( ,q )

 

 

 

 

 

 

 

 

ijtп

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ij

 

k

 

 

 

 

j

 

 

 

 

 

 

 

 

 

iju

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

i

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п f

 

' ( , q

);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iju

 

q

 

ij k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A ijt

 

 

k

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

ijt

 

 

 

iju

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Компьютерные технологии в приборостроении. Учеб. пособие

-69-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

Продолжение табл. 2.5

1

 

 

 

2

 

 

 

3

 

 

 

i

 

 

j

 

j

 

i

 

ijt

 

Aqki

Aqk

Aqk

 

 

 

 

 

j

 

ij

 

 

i

 

 

i,j

0

 

 

 

 

5

 

 

 

 

Aqk

 

 

ijt

 

j

A i

 

j

i

 

 

 

q

 

A

i

 

j

k

A ijt 0

qk

 

 

i

j

6

ijt ijt

 

 

q

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

A i

 

ijtп

 

 

 

 

 

hijt

 

 

A j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

h

f (q

)

 

 

k

 

 

 

 

 

XУП

 

 

 

 

 

 

 

 

k

 

 

 

 

ijt

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

k

 

j

 

 

 

ijtп

 

 

 

 

 

 

iju

 

 

 

 

 

 

i

 

 

 

XУП

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

ijt j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A l

 

 

 

 

 

 

 

 

 

 

Aq m

 

 

l

hlmu f (qk )

m

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qk

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lmu

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lmu

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

A lm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

h

(

ij

 

ijt

),

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

 

 

 

iju

 

 

ijt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A ijt п

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

ijt

 

 

 

 

iju

 

 

 

 

 

 

 

 

 

 

а)

ijt

kijt

,

 

 

а) п

kijt A lm ;

 

 

 

 

 

 

 

 

 

 

 

lm

 

 

 

lm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kijt

f (q );

 

 

ijt

 

 

 

lm

 

qk

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

lm

 

 

k

 

 

 

 

 

 

A lm A l A m ;

 

 

б) ijt

klm lm ,

 

 

 

 

 

 

 

 

 

k

 

 

 

 

k

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

 

 

б)

п

kijt

 

A lm

 

 

 

 

 

 

 

 

 

 

 

 

klm

f (qk );

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(kijt

)'

,

в) ijt

klm

lm ,

 

 

в) п

kijt A lm

 

 

 

 

lm

lm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

 

lm

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

klm f (qk )

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Компьютерные технологии в приборостроении. Учеб. пособие

-70-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

Окончание табл. 2.5

1

2

3

 

 

 

 

 

 

ijt

 

 

 

j

 

 

A i

 

 

 

 

 

п

 

 

 

 

 

 

A j

 

 

i

 

 

 

 

 

 

 

 

 

qk

 

 

 

 

 

 

 

 

 

 

 

 

qk

 

 

 

 

 

 

 

 

 

j

 

 

i

 

 

 

 

 

 

 

ijt

 

 

 

 

 

j

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hlmu

f (qk )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

ijv

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

A l

 

 

 

 

 

 

 

 

 

 

 

 

A m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

lmu

 

 

 

 

 

 

 

 

m

 

 

q

 

 

 

 

 

h

 

 

 

 

 

 

 

 

qk

8

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

lmu

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lmwп

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

lmu hlmu ( l m ) hlmu lmu

 

 

 

 

 

lmu

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п h'

 

lm

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lmu

 

lmu

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A lmu

п

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

lmu

 

 

 

lmw

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а) п

kijt

A lm ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

 

lm

 

q

 

 

 

 

 

 

 

 

 

 

а) ijt

klmijt lm ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijvп

0;

 

 

 

 

 

 

 

klmijt

f (qk );

 

 

 

 

 

 

 

 

 

 

 

 

б) п

 

kijt A lm

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

 

lm

q

 

 

 

 

 

 

 

 

 

б) ijt

klmijt

 

 

lm ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijuп

(klmijt )' lm;

 

kijt

f (q );

 

 

 

 

 

 

 

 

 

 

 

в) ijtп

klmijt

hlmu Aq lm ,

9

lm

 

 

k ijt

 

 

 

 

 

 

 

 

 

 

 

 

в)

ijt klm lm ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

f (qk );

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

kijt

h'

;

 

klm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijν

 

lm

lm

lm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г)

ijt

kijt

lm

,

 

 

 

 

 

 

г) ijtп

klmijt

hlmu Aqk lm ,

 

 

 

lm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijпν

klmijt

hlm'

lm

 

klmijt

f (qk )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(klmijt )' lmu

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2-й

случай.

 

Ветвь

 

описывается

зависимостью

ijt= ( ij,qk)

(см. табл. 2.5, п. 4). В этом случае вместо (2.27) получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

A ijt

f '

( ,q )

ij

f '

 

( , q )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

ij

k

 

qk

 

 

 

q

 

 

 

 

ij

 

k

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

ij

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.28)

 

f ( ij ,qk )

A ij

f ( ij ,qk )

.

 

 

 

ij

q

qk

 

k

 

 

При этом частная производная f ( ij ,qk ) ветви ijt исходной модели,

qk

Компьютерные технологии в приборостроении. Учеб. пособие

-71-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

как и второе слагаемое (hijt ij ) в выражении (2.25), играет роль зависимого

потокового активного компонента ветви iju, который вводится в преобразованную модель параллельно ветви ijt:

 

 

п

f '

( , q ).

(2.29)

 

 

iju

 

q

 

ij

k

 

 

 

 

 

k

 

 

 

 

Введенный активный компонент ветви ijt, как и пассивный параметр,

который по (2.28) равен

h

f '

 

( ,q )

, зависит от потенциальной пере-

 

ijt

 

 

ij

k

 

 

 

 

 

ij

 

 

 

 

 

менной величины ij исходной модели, но не зависит от переменных величин данной преобразованной модели. В связи с этим преобразованная модель становится линейной с переменными параметрами моделью даже в случае нелинейности исходной модели.

Теперь рассмотрим случай (табл. 2.5, п. 5, 6), когда в исходной модели находится только один активный потенциальный или потоковый независимый компонент (например независимый источник напряжения Xijt = E или тока ijt = I , см. рис. 2.7, б), моделирующий внешнее воздействие. В данном случае в преобразованной модели соответствующая ветвь исклю-

чается, а узлы i и j объединяются (для Xijt) или остаются необъединенными

(для ijt).

Это объясняется

тем,

что дифференцирование

выражений

(см. рис. 2.7, б) ij=Xijt (qk)

и ijt= ijt (qk) по параметру qk дает нулевой

 

 

 

 

 

 

ijt

 

 

 

 

результат, т. е.

 

ij

0

и

 

0 .

 

 

 

q

qk

 

 

 

 

 

k

 

 

 

 

Если

в

модели

присутствуют

зависимые активные

компоненты

(см. табл. 2.5, п. 7, 8), например, источник потокового воздействия, управ-

ляемый разностью потенциалов, т. е.

 

kijt

 

(см. рис. 2.7, в), и при

 

 

ijt

lm

ij

 

этом kijt

f (q ), то в результате дифференцирование исходного выражения

lm

k

 

 

 

 

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

п

kijt

A lm

и п

0 .

(2.30)

ijt

lm

qk

ij

 

 

Таким образом, проведя вышеописанные операции, дифференцирование

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

Используя описанный подход, можно составить систему уравнений по методу узловых потенциалов для преобразованной модели ([c]п[y]п = [x]п), в которой переменные величины исходной модели [y], входящие в зависимые активные компоненты преобразованной модели (см. выражения (2.26), (2.29),

(2.30)), будут играть роль воздействий, т. е. войдут в вектор [x]п. Решение по-

Компьютерные технологии в приборостроении. Учеб. пособие

-72-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

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

 

 

 

 

 

 

 

 

 

Т

 

 

j , ..., A

 

. На основе полученного вектора решения осуществ-

п A i , A

 

N

 

 

 

qk

qk

qk

 

 

ляются все необходимые действия [11], включая переход к относительной форме ФПЧ (см. выражение (2.16)).

В ряде случаев в процессе синтеза требуется получить ФПЧ разности потенциалов mn к параметрам различных ветвей модели. Допустим, что исходная модель будет линейной без зависимых компонентов. В этом случае дополнительная модель будет иметь входное воздействие (2.26). При вычислении ФПЧ, с принятыми допущениями, число дополнительных моделей будет равно числу исследуемых характеристик, для которых анализируется

ФПЧ Aqk mn . При этом дополнительные модели будут отличаться друг от дру-

га только значением и расположением входного воздействия (2.26), а искомые ФПЧ в каждой модели будут представлять собой разность потенциалов между одними и теми же узлами m и n. В силу линейности модели при решении каждой дополнительной модели можно принять значение входного воз-

действия равным единице, т. е. ijuп 1, тогда

mn

д

 

д

п

д

(2.31)

Aq

mn hijt ij mn iju mn ,

k

 

 

 

 

 

 

где дmn – разность потенциалов между узлами m и n дополнительной (преобразованной) модели при ijuп 1.

Таким образом, метод преобразованной модели позволяет получить за два решения системы уравнений вектор ФПЧ всех выходных характеристик к изменению одного параметра. Как правило, в этом случае в качестве параметра выступает входное воздействие (активный компонент топологической модели).

Для примера рассмотрим алгоритм получения вектора ФПЧ выходных характеристик модели тепловых процессов герметичного блока ПС к изменению одного параметра входного воздействия. Конструкция блока и его тепловая модель представлены на рис. 2.12. Тепловая модель (рис. 2.12, б) описывается системой уравнений (2.32).

Компьютерные технологии в приборостроении. Учеб. пособие

-73-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

Корпус транзистора

Окружающая среда

Кристалл Транзистор

Радиатор

Печатный узел

Корпус блока

Рис. 2.12. Эскиз конструкции герметичного блока ПС (а) его тепловая модель (б)

Компьютерные технологии в приборостроении. Учеб. пособие

-74-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

σ11

σ12

σ31

0

0

 

σ

σ

22

0

0

0

 

12

 

 

 

 

 

σ31

0

σ33

0

35

 

0

0

0

44

45

 

 

0

0

σ35

45

55

 

σ16

0

36

0

56

 

(2.32)

16 036

0

5666

t1t2t3t4t5t6

0

 

 

 

 

P

 

 

 

 

пу

 

 

 

 

0

 

n!

,

 

 

 

 

 

 

PVT

r! n r !

 

0

 

 

 

 

 

 

 

 

 

Tос

д

 

 

 

 

 

 

 

 

где

11 12 31 16 ;

22 12 ;

33 31 35 36 ;

...

;

 

66

σ

 

36

 

56

 

Ä

;

Д

– дополнительная проводимость

(103–104),

 

16

 

 

 

 

 

 

 

 

 

включаемая параллельно источнику Тос при его преобразовании в источник тепловой мощности; Pпу, PVT – мощности тепловыделения в печатном узле и транзисторе.

Для получения ФПЧ всех температур (t1, t2, , t6) к изменению, например, Pпу необходимо [8]:

1.Согласно п. 5 табл. 2.5 закоротить в схеме тепловой модели (уравнение (2.32) все источники потенциального воздействия (источник Тос на рис. 2.12, б).

2.Согласно п. 6 табл. 2.5 обнулить все источники потокового воздействия. Для уравнения (2.32) Pi = 0 / i = 1,6.

3.В интересующий узел подключить единичный источник потокового воздействия. Для рассматриваемого примера P2 = 1.

4.Решить систему (2.32) с новым вектором Pп [0,1,0,0,0,0]т . Результатом

решения будет вектор температур T п [t1п, t2п, t3п, t4п, t5п, t6п]т («п» знак преобразованной модели). При этом вектор решения будет вектором ФПЧ, т. е.

tп = At1

, tп = At2

, tп = At3

и т. д.

1

Pпу

2 Pпу

3 Pпу

 

5. Для получения SPТПУ необходимо использовать выражение (2.16).

Описанный выше метод преобразованной модели позволяет вести анализ ФПЧ для стационарных задач (анализ по постоянному току) в частотной и временной областях.

2.2.4.2. Методсопряженноймодели

Как отмечалось в п. 2.2.4, в основе метода сопряженной модели (присоединенной цепи [8; 16]) лежит принцип автономного анализа двух моделей, который представлен на рис. 2.13. Возможность построения сопряженной модели основывается на теореме Теллегена [16] или в более простом виде на теореме взаимности [11]. Теорема взаимности в общефизической трактовке может быть сформулирована следующим образом: «Если в элементе a сложной схемы дей-

ствует возбуждение F, вызывающее в другом элементе b этой системы отклик (реакцию) H, то, если перенести возбуждение F в элемент b, оно вызовет

Компьютерные технологии в приборостроении. Учеб. пособие

-75-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

вэлементе a тот же отклик H». В нашем случае под элементом a можно по-

нимать узлы i и j, а под элементом b – узлы m и n дополнительной модели (см. выражение (2.31).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 (единичное

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

воздействие)

 

X

Исходная

 

Y

 

 

 

 

Сопряженная

 

 

 

 

 

 

 

 

 

 

 

 

 

 

модель

 

 

 

 

 

 

 

модель

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

y

j a

1

 

2

Ay j , Ayj ,..., Ay j

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

q1 q2

qL

Рис. 2.13. Принцип автономного анализа двух моделей (а – вектор весовых коэффициентов)

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

Таким образом, если подключение к узлам i и j входного воздействияijuП =1 даст на узлах m и n разность потенциалов дmn , то перенос единичного

воздействия на узлы m и n (т. е. задание ijuс =1, где «с» – знак сопряженной модели) дает на узлах i и j то же значение разности потенциалов ijс дmn .

Полученную модель с единичным воздействием на выходных узлах m и n назовем сопряженной. Заменив в выражении (2.31) дmn на ijc , получим

окончательное выражение для расчета ФПЧ на основе линейной модели без зависимых компонентов:

mn

 

c

(2.33)

Aq

hijt ij ij ,

k

 

 

 

где hijt элементы матрицы параметров [С].

Таким образом, для получения ФПЧ одной выходной характеристики ( ) к первичным параметрам, содержащимся в различных ветвях топологи-

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

значения ( 1 и 2) перемножить в соответствии с (2.33). Производную hijt по первичному параметру qk легко найти путем непосредственного дифференцирования зависимости hijt (qk ) .

В табл. 2.6 приведены правила перехода от исходной модели к сопряженной для наиболее часто встречающихся случаев. Рассмотрим выражение

Компьютерные технологии в приборостроении. Учеб. пособие

-76-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

для вычисления ФПЧ относительно некоторых частных моментов, представленных в табл. 2.6.

Таблица 2.6 Правила перехода от исходной топологической модели

к сопряженной (присоединенной)

 

Ветвь исходной модели

 

 

 

 

Ветвь сопряженной модели

п/п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

 

 

3

 

1

ic

hijt

cj

ic

hijt

cj

 

 

 

i

 

 

 

j

 

i

 

 

 

j

 

 

2

i

 

hijt f(qk )

 

 

j

ic

hijt

 

 

cj

 

i

ijt f ij

 

 

 

j

i

 

 

 

j

 

i

 

j

ic h

f

 

 

cj

3

i

 

 

 

 

 

 

j

ijt

 

ij

 

 

 

 

 

 

 

i

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xijt

 

 

 

j

 

ic cj

 

 

i

 

 

 

 

j

 

 

 

 

 

4

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i, j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

 

j

c

 

 

 

c

5

i

 

 

 

 

 

 

 

i

 

 

 

j

i

 

 

 

 

 

 

j

i

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

ijt

 

h

 

j

ci

h

 

 

cj

 

i

 

 

ijt

 

 

 

ijt

 

 

 

6

i

 

 

 

 

 

 

j

i

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончание табл. 2.6

1

 

 

 

2

 

 

 

 

 

 

3

 

 

 

Компьютерные технологии в приборостроении. Учеб. пособие

-77-

 

 

 

 

 

 

 

 

2. МАТЕМАТИЧЕСКИЕ МОДЕЛИ

 

 

 

 

 

 

 

 

2.2. Функции параметрической чувствительности

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

 

c

 

 

h

ijt

cj

 

 

 

Xijt = klv

l

j

i

i

 

c

 

j

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X УП

 

 

ijt

h

 

 

7

 

 

 

 

 

h

j

 

l

 

 

lv u

v

 

 

 

 

 

h

ijt

v

l

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lvu

 

v

 

 

 

 

> >

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lcv u k lijv t

ijс t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijt

 

 

ic

 

 

hijt

 

 

 

cj

 

 

 

Xijt = klν lvu

j

i

ψijtc

 

 

 

 

 

j

 

 

 

X

 

 

 

 

 

 

 

 

 

i

 

 

j

 

 

 

 

 

 

 

 

8

 

i

 

 

 

 

 

УП hijt

 

Xlc

klijt

ψijtc

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

v

 

 

l

 

 

 

 

 

lvu

 

 

 

XУП

 

 

 

 

l

lvu

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

kijt

 

ic

 

 

 

 

 

cj

 

 

 

ijt

 

i

 

 

 

 

 

 

j

 

 

i

 

 

 

 

l

lv

j

 

 

hlνu

 

 

 

 

 

 

 

 

 

lc

 

 

 

 

vc

 

 

i

 

 

 

j

 

 

 

 

9

 

 

 

 

 

 

 

 

l

 

 

 

 

 

v

 

l

 

 

hlvu

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

v

 

c

= kijt c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lvu

 

lv

 

ij

 

 

 

 

 

 

 

ijt

ψlvu

 

ic

 

 

 

 

 

 

cj

 

 

ijt klv

 

 

i

 

 

 

 

 

 

j

 

 

i

 

 

 

 

j

 

X= klvijt ijc

 

 

i

 

 

 

j

 

 

 

 

10

 

 

 

 

 

 

 

lc

vc

l

hlvu

 

 

 

v

 

 

 

X УП

 

 

 

l

ψlvu

 

 

 

 

 

v

 

l

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

При анализе электромагнитных тепловых и механических моделей час-

то встречается ситуация, когда от первичного параметра qk одновременно за-

висит несколько параметров ветвей (топологической модели) с пассивными

компонентами:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hij

= hijt qk , ..., hplv = hplv qk и т. д.

 

 

 

(2.35)

 

 

Компьютерные технологии в приборостроении. Учеб. пособие

 

-78-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

Тогда в выражении (2.33) появятся новые слагаемые:

A mn = h'

ijt

 

ij

c

+ ...+ h

'

 

pl

c .

(2.36)

qk

 

ij

 

plv

 

pl

 

При практической реализации рассмотренного метода можно воспользоваться следующим обобщенным алгоритмом получения ФПЧ [16]:

Шаг 1. Решается исходная модель:

 

 

 

 

 

(2.37)

C Y

= X .

Шаг 2. Решается сопряженная (присоединенная) модель:

C Т Y

 

 

с ,

(2.38)

с = X

где X с формируется на основе правил, изложенных в табл. 2.5 (п. 4–6), а также с учетом номера узла, к потенциалу которого рассчитываются ФПЧ

(для номера j узла модели, к потенциалу которого вычисляются ФПЧ, задается единичное воздействие, т. е. xj = 1).

Шаг 3. На основе вектора Y с , используя выражения (2.33) и (2.36), вычисляются абсолютные ФПЧ.

Шаг 4. Переходят к относительной форме представления ФПЧ, используя выражение (2.16).

В заключение рассмотрим несколько примеров.

Пример 1. Рассмотрим алгоритм получения ФПЧ температуры t1 к изменению тепловых проводимостей 12 , 31 и т. д. (см. рис. 2.12, б). Алгоритм будем

рассматриватьвобобщенномвиде. Дляполучения At1ij необходимо:

1.Закоротить в модели источники потенциального воздействия (см. п. 4 табл. 2.5). Длятепловоймодели(рис. 2.12, б) закоротитьисточникТос.

2.Согласно п. 5 табл. 2.6 обнулить все источники потокового воздействия. Для системы (2.32): Pi = 0 / i = 1,6.

3.В интересующий узел модели подключить единичный источник пото-

кового воздействия. В рассматриваемом примере P1 = 1.

4. Осуществить анализ модели с новым вектором в правой части системы уравнений. Для нашего случая решить систему (2.32) с новым вектором

Pc =[1,0,0,0,0,0]т, используя выражение (2.38), т. е. предварительно транс-

понировав матрицу тепловых проводимостей и представив Pc в виде Pc =[ 1,0,0,0,0,0]т . Результатом решения будет вектор T с .

5. Согласно выражению (2.33) рассчитать абсолютные ФПЧ:

Компьютерные технологии в приборостроении. Учеб. пособие

-79-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

At112 = t1 t2 t1c t2c ; Aσt131 = t3 t1 t3c t1c ;

At1 = t t tc tc .

σ45 4 5 4 5

6. Для перехода к относительной форме ФПЧ следует использовать выражение (2.16).

Рассмотренные в примере ФПЧ At112 , At131 , ..., At145 не являются достаточ-

но информативными, так как при этом разработчик не может оценить влияние на выходную характеристику первичных параметров Q = q1, q2 , ..., qL ,

которые в рассмотренном примере могут играть роль коэффициентов теплопроводности материала радиатора λp : коэффициентов черноты поверхностей

блока или радиатора kч , геометрических параметров (размеры корпусаx , y , z ; площадь теплоотдающей поверхности ПУ – Snj; параметры оребрения: р – длина ребра; dp – толщина ребра и т. д.) и т. п. В этом случае следует полученные в примере ФПЧ умножить на составляющую hijt . При этом следует принять hijt = ij , продифференцировать все значения ij по ин-

тересующему параметру qk = lp , kч, x , y , z , p , dp , S и т. д., домножить полученные ранее ФПЧ на значение hijt согласно выражению (2.33) или (2.36).

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

ключить, что только σ12 f (Sп) . Однако 12 представляет собой две парал-

лельно соединенные тепловые проводимости 121 и 122.

Причем 121 = aл Sп , а σ122 = aк Sп , где 121 – моделирует тепловое излучение с поверхности печатного узла на корпус блока, 122 – моделирует конвективный перенос тепловой энергии с поверхности печатного узла на кор-

пус блока; αл – коэффициент лучистой теплоотдачи; αk

– коэффициент кон-

вективной теплоотдачи.

 

 

 

 

 

 

 

 

Учитывая вышеизложенное, применяем для вычисления, например,

ФПЧ

At1

выражение (2.36):

 

 

 

 

 

 

 

 

Sп

= t1 t2 t1c t2c

121 + t1 t2 t1c t2c

122

 

 

AS 10

= ASt1

= At1

aл At1 aк;

п

п

Sп

 

 

 

Sп

12

12

 

 

 

St1

= At1

Sп

.

 

 

 

 

 

 

 

 

 

 

 

 

 

Sп

Sп

t

 

 

 

 

 

 

 

 

1

 

 

 

 

Компьютерные технологии в приборостроении. Учеб. пособие

-80-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

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

ложим S = 2.

1

C1

= 1 Ф

2

 

 

I1 = 1 A

Uвых = U2

1 = 1 см

 

2 = 1 см

0

Рис. 2.14. Электрическая цепь

Система уравнений для представленной на рис. 2.14 цепи будет иметь

вид

 

 

 

 

 

 

 

sC

sC

U

 

 

 

1

 

 

 

 

 

 

 

 

 

1

1

 

 

1

 

 

 

1

 

 

.

 

 

 

 

 

 

 

 

 

sC1

2 sC1 U2

 

0

 

 

 

Подставив числовые значения, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

2

U

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

2

3

U2

0

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

0,6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

1

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U2

0,4

 

 

 

 

 

 

 

 

 

Анализ

сопряженной

 

модели,

 

согласно

 

выражению (2.38), при

 

 

c 0

1 т :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

2

U с

0

 

 

 

 

 

 

 

 

 

 

U с

=

 

0,4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1с

 

дает результат Y с =

 

 

1с

 

 

0,6

.

 

 

 

3

U2

1

 

 

 

 

 

 

 

 

 

U2

 

 

 

 

На основе векторов Y и Y с вычисляем ФПЧ:

Компьютерные технологии в приборостроении. Учеб. пособие

-81-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

AU вых = U2 =U1 U1c 0,24;

1 1

AU вых =U2 U2c 0,24;

 

2

 

 

 

 

 

 

ACU1вых = s U12 U12c

= s U1 U2 U1с U2c 0,08;

SUвых 0,24

1

= 0,24

1

 

= 0,6;

 

 

1

Uвых

0,4

 

 

 

 

 

 

 

 

SUвых = 0,24

2 = 0,6;

1

 

 

U2

 

SUвых = 0,08

C2 = 0,2 .

 

C1

 

 

U2

 

Пример 3. Рассчитать чувствительности выходного напряжения цепи, рассмотренной в примере 2, к температурам элементов. Допустим, что

температурный коэффициент емкости С1 чрезвычайно мал. Положим также

S = 2.

Кроме этого примем

1 =

1

 

,

R1 = R01 1+ aR tR1 t0 ;

2 = 1

, R2 = R02 1+ aR tR2 t0 ,

 

 

R1

 

где aR

температурный коэффициент

 

R2

 

 

 

 

 

сопротивления; R0 – значение сопротивления при нормальной температуре резистора t0; tR – рабочая температура резистора.

Исходная и сопряженная системы совпадают с системами, приводимыми в примере 2.

При вычислении ФПЧ будем применять выражение (2.33). С учетом вышеизложенного получаем

AU2

AU2

 

s1

;

AU2

AU2

 

σ2

.

tR1

 

 

tR1

t

R2

σ

 

 

tR2

 

1

 

 

 

2

 

Согласно приведенным выражениям вычисляем ФПЧ. При этом значения ФПЧ AU12 и AU22 возьмем из примера 2. В результате получаем:

AU2

= 0,24

 

 

R01 aR

 

 

;

 

 

 

 

 

 

t

 

 

 

 

 

2

 

 

 

 

 

R1

 

 

 

R1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

R

a

R

 

 

 

A

2

= 0,24

 

02

 

 

;

 

 

2

 

 

tR2

 

 

 

R

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

Компьютерные технологии в приборостроении. Учеб. пособие

-82-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

SU 2 = AU 2

tR1

;

U2

tR1

tR1

 

 

SU2

= AU 2

tR2 .

tR 2

tR2

 

U2

 

 

 

 

Если в рассмотренном примере вместо температур tR1 и tR2 использовать

температуру окружающей среды tос, то чувствительность U2 цепи (рис. 2.14) к tос будет вычисляться на основе выражений (2.33) и (2.36):

AU2

= AU2

 

1 + AU2

 

2 .

tоc

1

 

tоc

2

 

tоc

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

числения положим S = 2, а также будем учитывать, что при формировании системы уравнений паразитная проводимость [16] появляется в матрице [С] в следующих позициях:

 

i

… j

 

 

 

i

 

= 0

 

 

… = 0 .

 

 

 

 

j

 

= 0

… = 0

 

 

 

 

 

 

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

Учитывая данное обстоятельство, используем векторы Y и Y с из примера 2, на основе которых вычисляем ФПЧ:

AU2 s U

U с = 2 0,24

= 0,48;

1

1

 

Cп2

 

 

AU2 s U2 U2с = 0,48;

Cп3

AσUп22 U1 U2 U1с U2c = 0,04 .

Компьютерные технологии в приборостроении. Учеб. пособие

-83-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

п3 = 0

U1

1

C1 = 1Ф

U2 2

I1=1A

Cп2

= 0

Uвых = U2

 

 

 

 

 

 

Cп3 = 0

1 = 1См

 

2 = 1См

 

 

 

0

 

 

Рис. 2.15. Электрическая цепь с паразитными параметрами

Пример 5. Рассчитать ФПЧ выходного напряжения U2 электрической цепи (рис. 2.16) к параметрам C1, 2 и k 20 .

В дальнейших расчетах примем S = 2, а также будем учитывать, что источник потокового воздействия, управляемый разностью потенциалов

(см. п. 9 табл. 2.5) вносит компоненты в матрицу параметров [C] в соответст-

вии с рис. 2.17 [11; 16].

 

 

 

1

 

C1

= 1 Ф

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 = 1 см

 

 

 

 

 

 

 

 

 

 

Uвых = U2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Е1 = 1 В

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I2 = k20 U1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k20 = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

= 1 см

 

 

 

 

3 = 1 см

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

Рис. 2.16. Электрическая цепь с зависимым источником тока

Компьютерные технологии в приборостроении. Учеб. пособие

-84-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

l

l

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ij

= kij ( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

В случае j = = 0

Рис. 2.17. Источник потокового воздействия, управляемый разностью потенциалов (а), и вносимые им в структуру матрицы проводимостей [С] компоненты (б, в)

 

 

… l …

 

 

 

 

 

 

 

 

 

 

 

 

 

 

… l

 

i

 

+k

ij

k

ij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= C

i

+kij

 

= C

 

 

 

 

 

 

 

 

 

 

 

j

k

ij

+k

ij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б

 

 

 

 

в

 

 

Рис. 2.17. Окончание

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

+

 

+ sC

C

 

U

 

E s

 

 

1

2

 

20

1

1

 

1

 

 

1 1

.

 

σC1 + k

 

 

σ3 + σC1

U2

 

0

 

Подставив числовые значения, получим

4

2

U

1

 

1

 

1

3

 

 

 

.

 

 

U2

0

Компьютерные технологии в приборостроении. Учеб. пособие

-85-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

Решение представленной системы дает результат Y = 0,30,1 .

 

 

Решение

сопряженной модели

в соответствии с выражением (2.38)

 

4

1

U c

 

0

 

 

 

 

 

 

 

0,1

 

 

 

 

 

 

 

 

3

 

1c

 

 

дает результат Y c =

.

2

 

U2

 

1

 

 

 

 

 

 

0,4

 

 

На

 

 

 

 

с вычисляем ФПЧ:

 

 

основе векторов Y

и Y

AσU22 =U1 U1с 0,03;

ACU2 = s U1 U

2 U1с U2с 0,12 ;

 

1

 

AkU202 =U1 U2с 0,12.

Для перехода к относительной форме ФПЧ можно использовать выражение (2.16).

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

Полученная на основе вышеописанных методов матрица ФПЧ S QY по-

зволяет в процессе разработки прибора выбрать, например, наиболее влияющие на j-ю характеристику параметров, вариация которых позволяет обеспечить необходимые количественные показатели j-й выходной характеристики, качественно не изменяя при этом остальные выходные характеристики. Например, технический процесс имеет выходные характеристики,y1, y2, y3, y4 и изменяемые первичные параметры h1, h2, ..., h5, q1, q2 и qk = f(h) Предположим, что необходимо выявить параметры, в наибольшей степени влияющие на выходную характеристику y3. При этом одним из представленных выше способов была получена матрица относительных ФПЧ, имеющая следующий вид:

 

 

 

h1

h2

h3

h4

h5

 

 

 

 

0,09

0,14

0,18

0,47

0,30

y1

 

 

 

0,15

0,50

0,15

0,41

0,30

y2.

SY

 

 

 

 

 

 

 

 

y3

Q

0,05

0,76

1,25

2,1

0,12

 

 

 

0,01

0,45

0,10

0,52

0,22

y4

Анализ матрицы показывает, что наибольшими значениями ФПЧ будут –2,10; 1,25; 0,76, т. е. выходная характеристика y3 наиболее чувствительна к параметрам h2, h3, h4. Однако при изменении параметра h2 могут значительно

Компьютерные технологии в приборостроении. Учеб. пособие

-86-

2.МАТЕМАТИЧЕСКИЕ МОДЕЛИ

2.2.Функции параметрической чувствительности

изменяться характеристики y2 и y4. Поэтому целесообразно варьировать параметры h3 и h4. При изменении h3 характеристика y4 меняется незначительно, а y1 и y2 меняются меньше, чем в том случае, если бы варьировались другие параметры.

Компьютерные технологии в приборостроении. Учеб. пособие

-87-

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]