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

книги из ГПНТБ / Тепловые процессы при обработке металлов и сплавов давлением учеб. пособие для студентов металлург. спец. вузов

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

 

 

 

 

Po(p)exp

 

M2

 

M

 

 

 

 

 

 

 

 

• +

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mi

 

M

 

 

 

 

 

VMiAV T +

p ~ T + g

 

- f С exp

-

^

( /

f

+ ^ f )

]

+ « e x

p [

x ( / f +

, - f ) ]

 

 

 

 

 

 

 

 

 

 

(2.2.133)

Удовлетворяя

 

(2.2.133)

 

граничным

условиям

(2.2.128)-

(2.2.130), находим постоянные С и D:

 

 

 

С-

 

Кі(/>)

 

 

 

Po

(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

|

/

t

+

^

t

2

1 / т +

Ч ] / т

 

 

 

 

+

 

 

 

 

 

 

(p+Mg-g2)

Po (p)

ЛГ2

/ / M2

' M

/

Подставляя значения С и D в выражение (2.2.133), находим ре­ шение поставленной задачи в пространстве изображений

 

 

Ѵ(Х,

 

р) =

 

Ро(р)Х

 

 

ІГехр

- * (

і / т

 

' +

т

)

exp ( —

gX)

X

+

Mg-*2>(]/

 

 

 

 

/> + ЛІ£ — ^2

(P +

+ P +

 

 

 

Af2

 

Al

 

 

 

Kl (p)exp ! — X

 

 

+ p

M

 

 

 

 

+

(2.2.134)

 

 

VAf2

 

 

 

 

 

 

M

 

 

Перепишем

(2.2.134) следующим

образом:

 

где

V(X,

p) = Z(X,

p)+T(X,

p),

(2.2.135)

 

 

 

 

 

 

 

 

 

Z(X,

p) = P~o(p)Qx{X,

p);

 

70

 

g

exp

 

X (

, /~M2

 

M

\

 

 

exp ( — gX) .

 

 

 

 

\

Ѵ

^

+ р +

т

)

 

 

 

 

 

 

 

 

 

 

 

 

 

p + M g — #2

 

 

{P + Mg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.2.135a)

F(X,p)

определяется

формулой

(2.2.1136).

 

 

 

 

 

Обратное преобразование функции G{(X,

р)

 

получим с

по­

мощью таблиц операционного исчисления:

 

 

 

 

 

G:(X,

Ф) = \ {1 exp [ -

gX

-

(Mg-

g2) t] eric | " - ^ +

( | — g-) V't

+

 

g

exp [ -

(M - g) X-

(Mg -g2)t]

erîc

X

 

 

2Vt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'M

 

 

 

-exp[-gX-(Mg-g*)t]-

 

 

 

 

 

-iî-i)VT

 

 

 

 

 

 

 

 

 

M

 

erîc U - i / T

 

 

 

 

(2.2.136)

 

 

2 (M - g - )

2

 

 

 

или для неподвижной системы

координат:

 

 

 

 

 

 

О, (A-,, t ) =

j l

- l exp [ - gXx

+ (Л* ,g + £ * ) / ] erfc

 

M

 

+ 2 ( A f - * )

exp [ - (M-g) (Xi-Mit

l-+g)Vt

 

 

 

g

 

+ ^ ) ] e r î c

X

 

 

-exp[-gXl-\.(Mlg

1_2]Л

\

 

 

1

 

 

 

M

erîc

і / Л І Л .

 

 

2{M-g)

 

 

 

2 V T

 

 

 

 

Оригинал функции

Z(X, p)

получим, используя

дывания

 

 

 

 

+

+g2)t}-

(2.2.136a)

теорему запаз­

Z(X,

p) = Z(X,

Ф),

 

гд>е Z(X, ty) принимает значения:

 

 

для

(п + 1)-й

паузы

 

Z(X, Ф) = Ро2 2 ЮгІХ,

 

 

^-(т~\)хй1^-(т-\)хй1\~

m = l

 

 

 

— Ох (X, ф— т т 0 ) т] ( т т 0 ) } ;

(2.2.1366)

71

Z(X, ф)=Ро2

для

(п+1)-го обжатия

 

m=02 °ЛХ>

Ф— т т 0 t j ) ^(ф—mr0-г^-

 

п

 

 

 

— 2

^ 1

^ '

Ф - ^ о ^ ^ - ^ о )

(2.2.136B)

Принимая во внимание соотношения (2.2.115) и (2.2.115а), за­ пишем решение поставленной задачи в области действительной пе­ ременной:

 

 

 

для (п + 1)-й

паузы

 

 

v(X, t ) = P o 2 2

{Oi[Xv

 

t

- ( m - l

) t 0 - t j K i [ t - ( m -

l ) t 0 - f J -

 

m-l

 

 

 

 

 

 

 

 

~G1(XV

t - m t 0 ) ï i ( t - m t 0 ) } + K W A ' i , r ) f

 

(Ki2-Ki,)j2

{ ф І * і >

- « - ( / n - l J t o - t ^ T j l t - C m - l J t o - t J -

 

 

 

 

t - m T 0 ) Y ) ( t - m t 0 ) } ;

 

(2.2.137)

 

 

оля

(п+1)-го

обжатия

 

 

v(Xv

t ) = P o 2

Qi(Xv

X— mt0rjijfc

— т т 0 — t

j —

 

 

2

 

m2.-l O1 (Jr1 ,t-mt0 )ïi(t-mt0 )] + K i , ® ( A ' ï , t ) +

 

 

л

 

 

 

 

 

 

+ ( K i 2 - K i x ) 2 Ф ( * і , t — m t 0 - t , ) T ] ( t — m t o - t j -

 

 

л

 

 

 

 

 

 

 

 

2

Ф

^ ,

t - m t 0

) r i ( t - m t 0

)

(2.2.138)

Выражения (2.2.137)

и

(2.2.138) определяют

температурное по­

ле раската

в любом пропуске для случая, когда

начальное

распре­

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

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

72

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

В начальный момент времени температура по сечению тела име­ ет постоянное значение, а в течение времени t2 по его объему выде­ ляются равномерно распределенные тепловые источники мощно­ стью W2.

Если принять, что теплообмен между поверхностью тела и окру­ жающей средой происходит по закону Ньютона, а также если пре­ небречь компонентами температурного градиента по осям у и г, то для определения функции температуры тела необходимо решить дифференциальное уравнение теплопроводности (2.2.1) при краевых условиях (2.2.3), (2.2.4), (2.2.5), а также при следующем гранич­ ном условии:

 

}дТ(хь

 

p

 

 

« w i ^ ^ - ^ w ] -

 

( 2 - 2 - 1 3 9 )

 

 

 

 

 

 

 

где

 

 

 

 

 

Tc

(t) = Т, + {T2-Tx)

? (t);

a

= +

(a2

-

a,) <p (/);

cp(t) определяется выражениями

(1.6.3) и (1.6.4).

 

 

Перейдя

к подвижной

системе

координат и используя безраз­

мерные величины, можно поставленную задачу записать

следующим

образом:

 

 

 

 

 

 

 

 

 

 

 

д ѵ

д 2 Ѵ

+

м - ^ - - P o

(ф), ( 0 < Х < о о ) ;

 

(2.2.140)

 

 

дХ2

 

'

дХ

 

 

 

 

 

 

 

дѴ

 

 

= Ъ($)[Ѵх-о-Ѵе№->

 

(2.2.141)

 

 

дХ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ѵ(со,

ф)^оо;

 

(2.2.142)

 

 

 

 

 

дѴ

 

=

0;

 

(2.2.143)

 

 

 

 

 

дХ

 

 

 

 

 

 

 

Х=<х>

 

 

 

 

 

 

 

 

Ѵ(Х,

0) =

0,

 

(2.2.144)

где

 

 

 

 

 

 

 

 

 

 

остальные обозначения встречались раньше. Обозначим

q№ = b(Wx-o-VQW

(2.2.145)

73

и применим к системе уравнений (2.2.140) — (2.2.144) интегральное преобразование Лапласа — Карсона. После этого запишем

™ % Р )

где

^ М ^

Л

- р

Ѵ { Х ,

p)s=fTo(p),

 

dX

 

 

dV

(X,

p)

 

 

 

dx

 

X=Q

 

 

 

 

 

 

V

(со,

р)фсо,

 

dV (X,

p)

=

0,

 

 

 

dX X=oo

(2.2.146)

(2.2.147)

(2.2.148)

(2.2.149)

Общее решение дифференциального уравнения (2.2.146), удов­ летворяющее условиям (2.2.148) и (2.2.149),

Ѵ(Х,

р) =

А(р)ехр

+ / т + ' ) Ь ^

( 2 , 2 Л 5 0 )

Отсюда

находим

 

 

 

 

 

 

dV (X,

р)

 

 

 

 

 

dX

•х=о

\

2

 

Условие

(2.2.147) запишется так:

 

 

 

 

 

А(р) =

-Я(Р)

 

Решение, поставленной задачи в области изображений

получает

следующий вид:

 

 

 

 

 

 

е х р

X

Л12

 

 

 

 

 

.q(p)-MÛ.

(2.2.151)

Ѵ(Х,

р)=-

M + VМ2

 

 

 

 

 

Используя теорему о свертке оригиналов, запишем решение за­

дачи в области действительной переменной:

 

 

 

 

для

(п + 1)-й

паузы.

 

 

Ѵ(Х,

Ф ) = - J k + CV-Si)^

г{[Е-(т-\)х0-хЛ-

 

 

 

 

0 {

771 = 1

 

- 7 i ( / - m t o ^ 3 ( ^ Ф - 0 ^ ( 0 . 0 Л + | | 8 і + ( 8 2 - 8 і ) Х

74

X 2

• ' d i - ^ ~ ^ o - ^ i ] ~ ^ ( i - ' n t 0

) \ \ v l + ( V 2 ~ V l ) X

m = l

 

)

\

 

 

 

п

 

 

 

 

 

Х m2= 1 т( [/ (/и — 1>г0

— T J ч / и * 0

) Ф3(Х,

ф - / ) Л -

 

m =

l

1)тг0 —tjTjf^—(m— 1) тг0

— r j —

 

 

Р о 2 2

{[ф —(m—

(2.2.152)

 

 

— (Ф—отг0)т)(ф —/ит0 )};

 

оля ( я + 7)-го обжатия

Х Ф , ( * .

)-t)V(0,

t)dt+\\

\ + ( \ - \ )

У.

i j ^ - m t o - t i ) -

 

 

 

o l

 

 

L ^ o

 

m =

l

 

 

 

2

 

w c ö T i ) —

 

 

 

 

л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 ( ф - / я т 0 - ^ ) Х

/и=1

 

 

n

 

 

m=0

 

 

 

 

 

 

 

(2.2.153)

 

 

 

 

 

 

 

где

Х І ( Ф - ^ 0 — Т Х

) — 2

(Ц>

WTT0) TT] (ф-OTto

 

 

 

 

 

 

 

 

ф з ( * .

ф) =

Г / A f 2 ,

X2

Л*

\1

 

p [ - ( T ^ +

Яфl ? + T x ) J

 

 

e x

 

 

erfc

 

і/ф

(2.2.151)

 

 

2

2 К ф

 

 

 

 

 

 

 

 

 

Для перехода от функции Ѵ(Х,

г|з) к функции и ( Х ь т) необходи­

мо в выражениях

(2.2.152)

и (2.2.153)

заменить X на Хі—М/х и ф

на т. Функция Фз(А'і, т) будет иметь следующий вид:

 

 

 

 

г2

 

 

 

 

 

ехр (-4)

ü e r f c /

(2.2.155)

 

ф а № - *) =

 

 

 

 

 

Ят

 

2

\ 2

(/ x

 

 

 

 

~

В формулах (2.2.154) и (2.2.155) имеется неизвестная функция

1/(0, -ф). Дл я ее определения положим Х = 0.

Получаем:

для

(п + 1)-й паузы

 

1/(0, ф ) = - | | 8 1 + ( 8 а - 8 1

) 2

•nlt-im-^-X^—nd-mx^X

75

Х Ф 3 ( О , ф - о 1 / ( 0 , ю ^ н - J Ь і + ( * 2 — т <

і ) ѵ - * і і -

 

0 I

m = l

 

 

m = l

 

 

X Ф 8 ( 0 , *t-t)dt-Po22

{ [ Ф - ( w -

I J t o - t i ] il [ф— {m-

1 ) t o - t j -

- ( ф _ т т 0 ) т ] ( ф - т т 0 ) } ;

( 2 . 2 . 1 5 6 )

(Эля

(п + 1)-го

обжатия

 

Х Ф 3 ( 0 , ф - ^ ) 1 / ( 0 , / ) Л + |Ь1

+ ( 8 2 - 8 1 )

2

T i p - O T T o - t , ) -

 

 

 

О I

 

 

Lm=0

 

~2

l ^ - д а т о ) И^ + ^ - К , )

2 l ^ - w t o - f i ) -

-2

r,^ - /nt 0 )

} Ф З ( 0 , ф _ * ) Л - Р о 2

2

( t - m t o - t j X ,

 

 

 

J)

 

Lm=o

 

 

Х 7 1 ( ф - / п г 0 - г , ) - 2

-

'»To)

( ф — тт0)1,

( 2 . 2 . 1 5 7 )

где

 

 

m =

l

 

 

 

 

 

 

Af2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp

 

f e „

c

( f ^ ) .

 

 

Ф 3 (0,

ф) =

 

( 2 . 2 . 1 5 8 )

Выражения (2 .2 .156) и (2 . 2 . 157) являются интегральными урав ­

нениями Вольтерра

I I рода относительно температуры движущейся

поверхности V(0,

гр) или У ( М / Т , Т ) .

 

 

 

 

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

шить интегральные уравнения (2.2 .156)

и (2 . 2 . 157) . Подставляя за­

тем полученные значения функции Ѵ(0,

гр) в выражения

(2 . 2 . 152) и

( 2 . 2 . 1 5 3 ) ,

можно определить температуру в любой точке по сечению

раската.

 

 

 

 

 

 

 

 

Пример

Рассчитать температуру на поверхности раската * от момента выдачи из печи до конца второй паузы, если параметры теплообмена и прокатки имеют следую­

щие

значения:

 

 

 

1. Средняя температура по сечению заготовки

после вы­

 

дачи

из печи (Г0 )

.

1200° С

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

76

 

2.

Коэффициент

теплоотдачи

от

поверхности

раската

вт/м2

к окружающей

среде

в^, течение

пауз

(ai)

 

 

 

.

250

 

3.

Коэффициент теплоотдачи от поверхности раската

к по­

вт/м2°С

верхности валков

во

время

обжатия

(аг)

 

 

 

 

8-Ю3

 

4.

Длительность

паузы

(^і)

 

 

 

 

 

 

 

5 сек

 

5.

Длительность

пребывания

металла

в

очаге

деформа­

 

ции

(t2)

 

 

 

теплопроводности

прокатываемого

0,05 сек

 

 

6.

Коэффициент

метал­

 

ла

(X)

 

 

 

 

 

 

 

 

 

 

 

 

35

вт/м-град

 

7.

Коэффициент

температуропроводности

прокатываемого

 

металла

(а)

 

 

 

 

 

 

 

 

 

 

.

0,02 ж2

 

8.

Средняя

температура

окружающей

среды (Г]) . . . .

 

30° С

2 )

9.

Средняя

за

прокатку

температура

поверхности

валков

300° С

10. Средняя за прокатку скорость перемещения

поверхно­

 

 

сти

раската

 

 

 

 

 

 

 

 

 

 

 

20 м/ч

 

 

11. Среднее значение предела текучести

стали

(сг«) . . .

85 н/мм?

 

12.

Начальная

высота

раската

0 )

(#і )

 

 

 

740 мм

 

13.

Высота

раската после

обжатия

 

 

 

700 мм

 

14.

Скорость

конвективного

переноса

(si)

 

 

 

0,0

м/ч

По р я д о к р а с ч е т а

I.Определение безразмерных параметров теплообмена

Кр и т е р и й П о м е р а н ц е в а

 

 

 

 

 

Po (х) =

 

 

;

 

 

 

 

 

 

для

пауз

 

 

 

П71 =

0; Роі =

0,

 

 

 

 

 

 

 

для

обжатий

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я Л

 

 

 

 

 

 

740

 

 

 

 

 

 

 

 

as I n — ^

 

 

 

 

8,5-107.In

 

 

 

 

 

'

Po, =

 

Н\

 

=

 

 

 

 

700

 

 

=

п

 

 

 

 

I

8-103—0,25-103 \2

0,042.

 

1

*2ХА2 0 — ТЛ

 

 

 

 

 

г

- х 0

1 J

0,05-35 (

 

J

(1200 — 30)

 

 

 

Б е , з р а з м е р н а я с к о р о с т ь д в и ж е н и я п о в е р х н о с т и

р а с к а т а

 

 

 

 

s[

 

 

20

 

 

 

 

 

 

 

 

 

 

 

/И,

=

=

8-103-0,25-103

=

4,52.

 

 

 

 

 

 

1

 

ah

 

 

 

 

 

 

 

 

 

 

 

 

0

) 0 2

. _ _

_

 

 

 

 

 

 

 

 

Поскольку

M = M[ — Afj,

а Мх

= 0,

то M — М\ = 4,52.

 

 

 

 

П р о д о л ж и т е л ь н о с т ь

п а у з ы в б е з р а з м е р н о м

 

в и д е

 

 

 

 

 

 

 

/7,75-103 \2

5

 

 

 

 

 

 

 

 

т , = а А 2 Л = 0 , 0 2

 

 

• - —

=

1,360.

 

 

 

 

 

 

1

 

1

.

V

35

/

3600

 

 

 

 

 

 

П р о д о л ж и т е л ь н о с т ь

о б ж а т и я в б е з р а з м е р н о м

 

в и д е

 

 

•c2

=

ü;A2jf2 = 0,02

/ 7,75-103 \2 0,05

 

 

 

 

 

 

 

г —

- J да 0,014.

 

 

 

 

 

 

 

 

 

1

\

35

}

3600

 

 

 

 

 

 

Б е з р а з м е р н а я т е м п е р а т у р а

о к р у ж а ю щ е й с р е д ы

в т е ч е ­

н и е п а у з ы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^

= 4 F 4 r - = i,o.

 

 

 

 

 

 

 

 

 

 

 

 

То-Тг

 

 

 

 

 

 

 

 

 

77

Б е з р а з м е р н а я

т е м п е р а т у р а

п о в е р х н о с т и

в а л к о в

( с р е д -

я я я з а п р о к а т к у )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ѵ 2

=

Т0 7\

1200 — 300

=

0,770.

 

 

 

 

 

1200 — 30

 

 

 

 

 

 

 

Б е з р а з м е р н ы й

 

к о э ф ф и ц и е н т

 

т е п л о о т д а ч и

в

т е ч е н и е

н а у з ы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

et!

250

 

 

=

0,032.

 

 

 

 

 

 

 

0,2 0,1

8-103 — 250

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Б е з р а з м е р н ы й

 

к о э ф ф и ц и е н т

 

т е п л о о т д а ч и

в

т е ч е н и е

о б ж а т и я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а2

 

=

1,032.

 

 

 

 

 

 

 

 

 

а2—аг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

II. Расчет температуры поверхности

 

 

 

 

 

 

 

 

 

 

 

И н т е г р а л ь н о е

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

о п р е д е л е н и я

т е м п е р а т у ­

ры п о в е р х н о с т и

в

т е ч е н и е

п е р в о й

п а у з ы

получаем

из

общего

выражения

(2.2.156), положив п = 0 :

 

 

 

 

 

 

 

 

 

 

 

Ѵ(0, ф) =

—0,032 | " ф 3 ( 0 ,

Ф — r ) K ( 0 ,

t)dt

+

0,032

|*Ф3 (0,

<\> — t)dt.

(2.2.159)

 

6

 

 

 

 

 

ô

 

 

 

 

 

 

Необходимо подчеркнуть, что, поскольку определению подлежит

температура

новерхности раската в течение первой паузы,

под

1|) следует

понимать

время,

заключенное в интервале 0^г|)=£;ті, т. е. 0 ^ і | ) ^

1,360. Соответственно интеграль­

ное уравнение (2.2.159) позволяет определить

функцию Ѵ(0, а|)) только для вре­

мени первой паузы. Дл я получения данных о

температуре,

когда

1,360, сле­

дует учесть

в выражениях

(2.2.156) и (2.2.157)

изменение

числа

я, так что

общая

формула (вернее, интегральное уравнение) для определения температуры поверх­ ности будет иной, чем (2.2.159).

Принимая во внимание, что

I

ехр

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

=

 

 

 

 

 

 

2

 

l 2

 

 

 

] А г ( ф - г )

 

 

 

 

 

 

 

 

1

J M

 

 

 

 

 

 

Af

/ ф - ^

 

 

e r f ( ™ V l t - t ) + y ( Ф - г ) е г і с

2

 

 

меті[^

 

 

 

 

 

 

 

 

 

 

 

i r

 

Г

M 2

 

•t)

 

 

 

 

 

 

У Ф - ^ е х р

-

— (Ф

 

 

 

 

 

 

 

+ c.

 

(2.2.160)

 

 

 

 

 

У л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

запишем интегральное уравнение

(2.2.159) следующим образом:

 

 

 

 

ф) =

ср(ф) — 0,032

т

 

ф — О ^ ( 0 ,

t)dt.

 

 

 

Ѵ(0,

| " ф 3 ( 0 ,

(2.2.161)

где

 

 

 

 

 

 

 

 

 

 

Mi

, \ 1

 

 

 

 

 

 

 

 

 

V

Ф ехр

 

1

 

 

Af

 

 

 

Т

*

?(ф) = 0,032

erf —

 

 

У

 

 

У ф — — фerfc

 

 

 

 

 

M

2

Y

2

Y

I

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2 . 2 . 161а)

78

Уравнение (2.2.161) будем

решать

методом последовательных

приближений*.

Д л я

этого искомую функцию

 

Ѵ(0, -ф) представим

в виде

ряда [63, 65—67]

 

 

Ѵ(0,

Ф) =

2

 

г =

0,032).

 

(2.2.162)

 

 

 

 

 

 

й =0

 

 

 

 

 

 

 

Подставляя

ряд (2.2.162)

в

уравнение

(2.2.161) и

сравнивая

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

при

одинаковых

степенях

г,

находим

последовательные

приближения функции

Ѵ(0,

і|>):

 

 

 

 

 

 

f (Ф);

 

 

 

(2.2.163)

 

 

 

 

 

 

Л) (Ф) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F1 (ф) =

-

 

0,032 j

Ф3

(0, ф -

о

(/) dt;

(2.2.164)

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

Ф

 

 

 

 

 

 

^ 2 ( Ф ) = - | ф 3 ( 0 , Ф - О Л ( О Л .

(2.2,165)

ит. д.

Дл я £-го приближения имеем

 

т

 

 

 

 

 

 

 

П ( Ф ) = - | ф 3 ( 0 , Ф - О ^ - і С ) Л .

 

(2.2.165а)

 

о

 

 

 

 

 

 

 

После определения k функций

Fk(ty)

находим

функцию У(0, ф) в &-м при­

ближении:

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ѵ*{0,

ф) = 2^П(Ф).

 

(2.2.166)

Безразмерная длительность паузы

составляет

т і = 1,360. Дл я

расчета разби­

ваем этот интервал на 4 отрезка

(фі = 0,34;

ф>2=0,68; фз=1,02;

ір4 = 1,36).

Подставляя эти значения \|) в

формулу

(2.2.161а), а также

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

Л1 = 4,52, определяем нулевое приближение

функции

Ѵ(0, •ф):

 

 

Г 1

/ 4 , 5 2 , /

\

4,524,52

/ 4 , 5 2 , / — — Л

Т (фО = 0.032 { ^ e r f

l / 0 , 3 4 j -

0,34 erf с {-j-

V

0.34J +

]/o,34 exp

4,52\2

 

 

 

 

 

 

 

 

0,34

 

 

 

 

 

 

 

 

 

 

 

=

0,032 [0,221 erf (1,32)

— 0,767 erfc (1,32) + 0,330 exp( — 1,74)] = 0,032 (0,221-0,938 — 0.767 X X 0,062 + 0,330-0,176) = 0,007;

<p (ф2 ) =

0,032 [0,221 erf (1,86) — 1,534erfc (1,86) + 0,465 exp ( — 3,45)] =

=

0,032 (0,221-0,990— 1,534-0,010 +

0.465 0,032) = 0.007.

Аналогично определяем

ф(фз) и <р(ф.і):

 

 

 

f

(Фз) = 0,007; <р4) =

0,007.

 

Д л я уточнения значений функции ср(ф) в

интервале О^грsg0,34 рассчиты­

ваем ф ( ф 5 ) ,

ф(фб), ф(Ф?) (ф5=0,08; г|)6 =0,16; ф 7 = 0 , 2 4 ) .

Получаем:

 

9 5 ) = 0,005; ер 6 ) = 0,006; <р 7 ) =

0,007.

* Процесс последовательных приближений при решении интегральных урав­ нений Вольтерра I I рода является сходящимся.

79