Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Математическое моделирование процессов в машиностроении..pdf
Скачиваний:
46
Добавлен:
15.11.2022
Размер:
14.87 Mб
Скачать

3.2. АНАЛОГОВОЕ МОДЕЛИРОВАНИЕ ФИЗИЧЕСКИХ ПОЛЕЙ В ЭЛЕМЕНТАХ КОНСТРУКЦИИ

ТЕХНИЧЕСКИХ ОБЪЕКТОВ

3.2.1. Общие принципы моделирования полей

Стационарные физические поля описываются однород­ ными или неоднородными эллиптическими дифференциальны­

ми уравнениями в частных производных:

 

 

 

 

 

V (co (O V O O = 0,

 

(3.1)

 

 

V(co(^) VO(O)

Ф О -

 

(3.2)

Уравнения (3.1) и (3.2) включают в себя потенциал -

ска­

лярную функцию Ф(40> зависящую от координат £,г (где г -

ин­

декс,

равный

1, 2, 3). Для прямоугольных

координат

 

- х ,

=у,

= z,

для цилиндрических

= z, ^

= г, £,3 = а

и т. д.

Также уравнения (3.1) и (3.2) включают характеристику среды со(£,г) и мощность распределенных источников поляД^г, Ф(^г)).

При co(^r) = const уравнения

переходят в хорошо

известные

уравнения Лапласа и Пуассона:

 

 

> е

оои

(3.2)

АФ(^) = ~ / ( ^ ,Ф ( ^ ) ) .

(3.4)

СО

 

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

производные по времени:

 

 

V(co(£r, 1)'ЧФ&\1)) = а (£ , 0 9Ф(^

 

(3.5)

д2Ф(Ъг

t)

(3.6)

V(co(^, 0 V ® ( ^ ,0 ) = ^ r, 0

Произвольное поле описывается общим уравнением вто­ рого порядка:

V(co($r, /,Ф)УФ@г,0 ) = М ,

(,Ф) +а ^ г,1,Ф)

дФ(1г t)

- ; +

 

 

dt

 

 

(3.7)

+ а д .,,,Ф)^

£ ) .

 

Модель должна включать в себя также необходимые крае­ вые условия. Во-первых, должна быть задана область D, ограни­ ченная поверхностью (кривой) Г, в которой определяется реше­ ние, во-вторых, должны быть заданы условия на границе Г этой области. Могут быть заданы граничные условия первого, второ­ го и третьего рода:

 

Ф И /,(Г );Ф И /,(Г ,/),

 

(3.8)

 

дФ

дФ

 

 

(3.9)

 

= Л (Г ),^ Ч 0 =/ 2(Г,0,

 

дп

 

дп

 

 

 

 

Л

/.5Ф

= /з (П ,

 

(ЗЛО)

 

аФ +Ь

 

 

 

 

дп

 

 

 

дФ

(

 

 

дФ

= f 3{ r , t ) . ( 3.11)

аФ +Ь

= /э(Г)=> a(l)0(t) + b(t)— (t)

дп

I

 

 

дп

 

 

Для нестационарных полей, кроме того, должны быть за­

даны начальные условия,

характеризующие

состояние

поля

в начальный момент времени:

 

 

 

 

 

 

 

 

= № ) ■

 

(3.12)

/=о Некоторые уравнения полей с указанием конкретных фи­

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

76

Аналогичные уравнения полей различной физической природы

Область физики

Теплопередача (поле температур в твердом теле)

Гидродинамика (течение идеаль­ ной жидкости)

Электростатика

(электростатиче­ ское поле)

Электрическое поле в проводни­ ках

Уравнение

(закон)

Уравнение

теплопро­

водности

Закон

Бернулли

Теорема

Гаусса

Закон

Кирхгофа

Форма уравнения

V(X(^,7)-V7) =

=Т)

V(p(S', Ur)-VUr) =

= <7ЙГ.ОД

V(e(^, U)-VU) =

= g ( ^ U )

V(a(^,<p)-V(p) =

=<p)

Условпые обозначения

Г-температура;

X - теплопроводность; соинтенсивность объемноготепловыделения Ur - гидродинамический потенциал скоростей; р - плотность среды;

q ~ мощность распреде­ ленного источника жид­ кости

U - электростатический потенциал; е - диэлектрическая по­

стоянная среды;

q - плотность распределе­ ния объемных зарядов (р - потенциал электри­ ческого поля; ст - электропроводность среды;

i - распределенные токи, вводимые в среду

3.2.2. Особенности построения моделей

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

Каждое из уравнений, приведенных в табл. 3.1, включает в себя четыре величины, занимающие в этих уравнениях одина­ ковые позиции. Это потенциальная функция (Т, U, Ur, £/„, (р...),

характеристика среды (X, е, р, ц, а..)» мощность распределен­ ных источников поля (w, q , q r, и, наконец, геометрический масштаб модели (L, Г). Величины, занимающие в уравнениях одинаковые позиции, называются сходственными величинами, а для характеристики их отношений вводятся коэффициенты аналогии:

м л

a

IL

6

а

ф

ф

ф

 

I

I I

(3.13)

 

 

С помощью

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

аналогии можно выразить

и заменить величины изучаемого (например, температурного) поля величинами модельного (электрического) поля:

Т = М<рф, со = М^, L = M\l; X = Л/Стст,

(3.14)

дТ _ M v Эф

дТ _ Му дц>

дТ _ М 1( дер

 

дх М,

дх ’

ду М, ду

& М, dz

 

 

 

М

 

(3.16)

 

 

У7’ = ^ У ф .

 

 

 

М,

 

 

Тогда уравнение для стационарного температурного поля

(см. табл. 3.1) можно записать в виде

 

 

М М

 

 

(3.17)

ф

сг

 

 

м,м, У(а(^,ф)Уф) = /(^,ф ).

Важно подчеркнуть, что уравнение по-прежнему остается

уравнением теплопроводности, ведь изменилась только форма записи тепловых величин (Т, К, со, L). Вместе с тем уравнение теплопроводности в приведенной форме полностью совпадает

суравнением электрического поля при условии, что

ММ а

фО

м м ?

Данное соотношение выражает условие аналогии натур­ ного и модельного (в рассматриваемом примере температурного и электрического) полей и носит название индикатора аналогии.

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

<йг,7), /(^г,ф). Целесообразно заменить эту величину некоторой сосредоточенной величиной (W, I), характеризующей мощность источников (тепло- и токовыделение) в некотором элементар­ ном объеме (V, v):

W= (oV, f= iv.

(3.19)

Отношение величин W и / характеризуется коэффициен­

тами аналогии Mi и М,:

 

 

W

mV

(3.20)

М ,= — =— = М]М , ,

I

iv

 

 

М

 

 

м,

(3-21)

м *= т к -

Индикатор аналогии может быть представлен в виде

М щМ,Мв

(3.22)

М,

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

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

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

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

моделирования (рис. 3.2):

 

 

 

d = -^-rD \ г = — Г ; v =-X rV

(3 .20)

м?

м,

м)

 

Рис. 3.2. Задание области при изучении поля:

а- натурная область,

б- модельная область

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

(г = — Г). Граничные уело-

м,

вия первого рода для мо­ дельной функции задаются соотношением

г

Л

Ф|г = /( Г ) = > Ф|1._г.

j

М,

При этом на модели в каждой точке границы должен быть задан электрический потенциал <р. Точно также граничные усло­

вие второго рода задаются соотношением

(

дФ

= Mj_

дп т

М

(3.22)

/ 2

 

М,

 

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

границы. В силу очевидных соотношений:

 

(3.23)

; 1 = А = - а 5 — ,

дп

дп

 

гдеj - плотность тока, 7 - ток, можно получить уравнение

 

 

г т л

(3.24)

/ = А = - Ч ^ - Л

К**,;

К

 

где знак «-» учитывает направление тока.

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

быть записаны так:

 

 

 

 

 

 

 

Л идФ

 

 

.

 

,. бф

 

Г Л

аФ +Ь----

= /э(Г):

а

ш + Ь

 

= /э

• (3.25)

дп ,

 

 

дп

 

М,

 

 

 

 

 

 

■=.!_

V / /

 

 

 

 

 

 

 

М,

 

Условие (3.25) в силу (3.23) может быть записано в виде

f

 

и \

 

 

1

 

Г ^

(3.26)

ф -

аЛ7,а5г j

 

 

аМ. -/з

М,

 

 

и.

 

 

 

 

 

 

 

 

Если положить

 

 

 

 

 

 

1? =

Ь

 

 

1

f

Г Л

 

аМ,о5г

^>0

аМ,

3

 

 

 

 

 

 

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