Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Семененко В.А. Вычислительная техника в инженерных и экономических расчетах учеб. пособие для студентов всех специальностей

.pdf
Скачиваний:
5
Добавлен:
23.10.2023
Размер:
8.7 Mб
Скачать

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

Рис. 1—54. Установка начальных условий с при­ менением дополнительного суммирующего усили­ теля

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

В АВМ предусмотрена также установка тре'буёмых коэф­ фициентов передачи решающих элементов. Обычно такая установка производится на линейных .решающих элементах, так как нелинейные имеют, как правило, фиксированный коэффициент передачи. В линейных решающих элементах предусматривается возможность изменения коэффициента передачи в широких пределах (от 0,001 до 100). При этом ступенчатое изменение коэффициента передачи достигается с помощью многодѳкадных делителей, включенных на входе и в цепи обратной связи, а плавное регулирование — с по­ мощью прецизионных многооборотных потенциометров.

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

70

імгом f—i Ifincp

■ О - * ^-f----Il-

Рис. 1—55. Способы установки коэффициента передачи решающего элемента

§ 5—4. Методы задания функций возмущения

При решении на АВМ неоднородных дифференциальных уравнений вида

dnx

dn~lx

 

d tn

+ O-t 1 d t n~1 + • • • + d t

+aüx-bQy{t)

необходимо представить в виде напряжения U(t) функцию возмущения y(t). Задание функции y(t) может быть осуще­ ствлено двумя методами:

1.Представлением функции y(t) в виде решения соответ­ ствующего дифференциального уравнения.

2.С помощью функционального элемента.

Рассмотрим первый метод на нескольких примерах.

а) у — sin at.

 

Очевидно, что гармоническая функция sin ш t

является

решением дифференциального уравнения •

 

~ ^ - + * 2У = 0.

(1-48)

аР

 

Структурная схема, составленная по уравнению

(1— 48),

приведена на рис. 1—56. Устанрвив на интеграторах требуе­

71

мые

начальные

условия (при

t — О, Уо = 0, Уо

= 1), получим

на выходе схемы напряжение, моделирующее

в определен­

ном

масштабе

функцию у =

sin со t.

 

 

 

Рис.

1—(56. Структурная

схема

 

 

набора дифференциального

урав­

 

 

нения у" +

= О

 

б) y — ct,

tт.

е. возмущение

носит линейный характер.

Так как

j

с- dt =

с - 1, то

функцию

y = c-t можно

 

о

 

 

(рис. 1—57), на вход кото­

получить с помощью интегратора

рого включен источник постоянного напряжения.

в) y(t) =

Сі • t\

 

 

 

+L

Рис. 1—57. Схема реализации функции у = Сі t

Заданная функция может быть

реализована, очевидно, с

 

t

помощью схемы рис. 1—58, так как

Jctdt = схі2.

 

о

‘ 1

Рис: 1—58. Схема реализации функции

У = с і • Р

Рассмотрим теперь метод задания функции возмущения с использованием функционального элемента. Напряжение U(t) = f(t) можно получить на выходе функционального эле­ мента, если на его вход поступает напряжение t/BXs ^

72 ч

(рис. 1—59). В некоторых случаях возмущение, представля­ ющее заданную функцию времени, получают с помощью

" Г

Рис. 1—59.

Схема получения функции возмуще­

 

 

ния f(l)

с

помощью

функционального

элемента

 

введения,

так

называемого,

вариатора

коэффициентов

(рис. 1—60). В этом случае y(t)

целесообразно представить

в виде y(t)

= утах -а(/),

где у тах — максимальное

значение,

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

y(t)

на

рассматриваемом

отрезке

времени.

 

 

 

 

 

 

 

 

§ 5—5. Оценка погрешности решения задач на АВМ

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

 

у" + ш=у=Ц

(1—49)

с начальными условиями при

t = 0; у = С; у'

= 0.

Структурная схема решения данного уравнения приведена на рис.

1—61. Пусть передаточные функции отдельных

элементов схемы равны

а\(р), 0г(р)> °з(р) ■Тогда для

структурной схемы можно записать уравне­

ние

Y(p)

 

а^(р)-а3(р) а,(р) Y (р) =

 

ИЛИ

 

 

[1 — аг{р)а2 (р)а3 (p)]Y (р) = 0,

(1—50)

где Y — машинная переменная, связанная с

величиной у

в уравнении

определенным масштабным коэффициентом.

 

 

.

,

. 73

Для идеальных решающих элементов, не вносящих погрешности при выполнении математических операций,

(р) = — — ;

а2(рУ= -

; а3{р) = — а3,

Р

Р

Рнс. 1—61. Структурная

схема для получения

гармонических

колебаний

и уравнение (1—50) запишется следующем виде:

где

 

СРг + o,fl2a3)y(p) = О,

(1-51)

1

 

1

 

Ro

а, =

а. =

а3=

 

RXCо

 

R3Cо

 

R3

Точное решение заданного уравнения на АВМ можно записать в виде

 

 

y 1 =

csinQ<M

 

(1—52)

где 0. = а 1а г а 3\

tM— машинное

время..

 

 

Рассмотрим, каково будет отклонение от точного решения с учетом погрешностей, вносимых интеграторами из-за конечного значения коэффи­ циента усиления Ку УПТ. В этом случае передаточные функции интегра­ торов будут иметь вид (§ 2—5):

 

К у

 

 

 

 

 

Ку

(1 -53)

ах{р) = — ---------- Tr— Г Ѵ

I а2ІР) =

 

 

K y + l

 

Ky+l

 

 

 

1 + p

 

1 + P:

a,

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставив (1—53) в (1—51), получим

 

 

 

 

 

ІЧ-

 

Ау я7за

 

 

 

У(Р) = 0.

(1-54)

К у +

1

,

К

у +

1

 

1 + Р

 

 

al

I

1 + р -

у

 

 

 

 

\

 

а2

 

 

 

Учитывая, что Ку > 1 ,

уравнение

(1—54)

можно записать в виде

 

У

Ц]

 

 

 

 

 

 

р 2 + р

+ axa3a^jY(p) = 0

 

 

 

 

 

Ку

74

или

 

 

 

 

 

 

 

 

 

где

 

Y" + aY +

й К = О

 

 

(1-55)

 

 

 

 

CLо

 

 

 

 

 

 

a =

ß

j

 

 

 

 

 

 

----------— .

 

 

 

 

 

 

 

Ky

 

 

 

 

Таким образом, реально машиной моделируется уравнение (1-55)

вместо (1—49). Представим решение этого уравнения

в виде

 

 

Y =

Yx + aq,

 

 

(1-56)

где

a-q— величина,

определяющая

отклонение

от точного решения

У, =

csin й tM.

 

 

 

 

 

 

 

 

 

Подставив (1—56) в (1—55), получаем

 

 

 

 

 

У," +

ад" + а У", + o*q' + Й2^

+ ай2? = 0.

(1—57)

 

Пренебрегая в уравнении (1—57) членом

а2?'

(ввиду

малости а),

получим уравнение

 

 

 

 

 

 

 

 

 

Y'\ + й*У, + aq" +

ай2? +

а У / =

0.

(1—58)

 

Функцию q, определяющую

погрешность,

можно

найти

путем реше­

ния уравнения

 

 

 

 

 

 

 

 

 

,

?" +

й2? +

У'1 =

0.

 

 

 

 

Учитывая, что

У', = с Й cos й tK,

имеем

 

 

 

 

 

 

?" + й2 q = с й cos ч) .

 

(1—59)

q' =

Решив уравнение (1—59) с начальными условиями при

fM= 0:? = 0,

0, получим:

 

 

 

n /

 

 

 

 

 

q = — с

 

 

 

 

 

sin й f ,

 

 

 

и погрешность в машинном решении равна

 

 

 

 

 

 

аі + а2

J M

 

 

 

 

 

aq = — с

Ку

2

• sin й fM.

 

 

Так как точное решение должно быть равным Уі = с sin й /м, то отно­ сительная погрешность

0,2

■и

(1-60)

2Ку

 

Из выражения (1—60) следует, что

для уменьшения погрешности

желательно уменьшение коэффициентов передачи интеграторов аі, а2 и уменьшение времени решения tM.

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

§ 5—6. Состав и компоновка блоков АВМ. Порядок набора и решения задач

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

75

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

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

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

Перед началом решения задачи на АВМ следует осуще­ ствить проверку 'Нулевого уровня решающих усилителей. Настройка нулей производится с помощью Специальных по­ тенциометров, включенных в цепь связи между первым и вто­ рым каскадами УПТ. В процессе настройки добиваются минимума напряжения на выходе УПТ при отсутствии сиг­ нала на ого входе. Наряду с индивидуальной установкой ну­ левого уровня усилителей часто используется групповая установка, которая обеапечивает автоматизацию этого про­ цесса. 1

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

76

Важным достоинством аналоговых вычислительных машин является наглядность получаемого решения, так как АВМ обычно 'Снабжается различными средствам« для наблюдения и записи 'процесса решения. Для визуального наблюдения решения применяются катодные осциллографы (типов И—4, И—5М, И—40 и др.). Для регистрации медленных процессов с частотой, не превышающей 1 -г 2 гц, используются различ­ ные самопишущие приборы, например, типа регистрирующих автопотенциометров (ЭПП, MC и др.). При более высоких частотах используются шлейфовые осциллографы, например, типов МП02, К—'105 и др., позволяющие воспроизводить на фотобумаге одновременно до 12 кривых.

ГЛАВА VI

ЭЛЕКТРИЧЕСКИЕ МОДЕЛИ—АНАЛОГИ

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

турного

типа

или аналоговых

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

машин

в моделях-аналогах задача решается в ее физической

поста­

новке на

основе

систем аналогий,

существующих

между

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

§ 6—1. Подобие явлений. Константы и критерии подобия

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

77

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

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

- ^

= N m, 4 - = ^ >

4 - = ^

(1~ 61)

т2

 

h

 

где т — масса,

I — линейное перемещение; t — время.

Можно показать, что силы F\ и F2, действующие на мате­ риальные точки т\ и т2 и вызывающие их подобные движе­ ния, также подобны, то есть отношение этих сил может быть выражено константой подобия.

По второму закону Ньютона имеем

 

F2 = m,

 

 

(1-62)

Подставляя выражения (1—61) в

уравнения

(1—62),

получим

т,ЫтЫ[ dH,

 

 

 

 

 

 

N )

dt\

 

 

откуда

 

 

 

 

^2

N mN t

const

" N F ,

(1-63)

N i

 

 

 

что и требовалось показать.

 

 

 

Перепишем равенство (1—63) в виде

 

 

 

Ngi NI

_ I

 

(1-64)

 

 

 

 

N't N F

Зависимость (1—64) между константами подобия назы­ вается индикатором подобия. Она указывает на то, что между константами подобия, как и между физическими вели­ чинами, характеризующими протекание процесса, существует

78

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

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

Если подставить выражения (1—61) в равенство ^1—64), то получим для рассмотренного ранее примера соотношения между величинами масс, сил и перемещений

_ TtXyJ'l F2t\

которые являются критерием механичеокого подобия.

§ 6—2. Электрические модели механических систем

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

Рассмотрим эту теорию применительно к электрическому моделированию процессов в упругих механических системах.

На рис. 1—62 изображена механическая система, состоя­ щая из груза с массой т, .подвешенного с помощью пружины к неподвижной опоре. Груз может перемещаться лишь по вертикали вдоль ограничителей, создающих трение. Для

данной механической системы можно составить

уравнение

динамического равновесия Далам'бера

 

Qm Qs + Qe — Q(0>

(1—65)

 

 

 

Рис. 1—62.

Упругая

ме­

 

 

 

ханическая

система

с

 

 

 

одной степенью свободы

где Qm = m

d?x

 

■сила инерции;

 

 

dt*

 

 

 

 

Q - S

-

-

'Сила трения;

 

Ws

dt

-

 

 

 

79

Соседние файлы в папке книги из ГПНТБ