Температурні режими дискового екструдера
З метою вдосконалення управління температурними режимами у процесі виробництва полімерних композиційних матеріалів на лініях грануляції досліджено розподіл температури і розроблено математичний опис нестаціонарних процесів, що проходять в проміжку між диском та корпусом дискового екструдера....
Збережено в:
Дата: | 2003 |
---|---|
Автори: | , , , , , |
Формат: | Стаття |
Мова: | Ukrainian |
Опубліковано: |
Інститут бiоорганiчної хiмiї та нафтохiмiї НАН України
2003
|
Онлайн доступ: | http://dspace.nbuv.gov.ua/handle/123456789/3754 |
Теги: |
Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
|
Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Цитувати: | Температурні режими дискового екструдера / Л. Шкарапута, Л. Митрохіна, І. Морозова, Л. Тищенко, В. Даниленко, В. Абросімов // Катализ и нефтехимия. — 2003. — № 11. — С. 42-45. — Бібліогр.: 4 назв. — укp. |
Репозитарії
Digital Library of Periodicals of National Academy of Sciences of Ukraineid |
irk-123456789-3754 |
---|---|
record_format |
dspace |
spelling |
irk-123456789-37542009-07-10T12:00:32Z Температурні режими дискового екструдера Шкарапута, Л. Митрохіна, Л. Морозова, І. Тищенко, Л. Даниленко, В. Абросімов, В. З метою вдосконалення управління температурними режимами у процесі виробництва полімерних композиційних матеріалів на лініях грануляції досліджено розподіл температури і розроблено математичний опис нестаціонарних процесів, що проходять в проміжку між диском та корпусом дискового екструдера. С целью усовершенствования управления температурными режимами в процессе производстве полимерных композиционных материалов на линиях грануляции исследовано распределение температуры и разработано математическое описание нестационарных процессов, которые проходят в промежутке между диском и корпусом дискового экструдера. The allocation of temperature has been studied for the purpose of improving the control of temperature conditions while producing the polymeric composites on granulation lines of proceeding and the mathematical specification of non-stationary processes in a gap between disk and body of disk extruder has been developed. 2003 Article Температурні режими дискового екструдера / Л. Шкарапута, Л. Митрохіна, І. Морозова, Л. Тищенко, В. Даниленко, В. Абросімов // Катализ и нефтехимия. — 2003. — № 11. — С. 42-45. — Бібліогр.: 4 назв. — укp. http://dspace.nbuv.gov.ua/handle/123456789/3754 665.7.038.5:519.876.2 uk Інститут бiоорганiчної хiмiї та нафтохiмiї НАН України |
institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
collection |
DSpace DC |
language |
Ukrainian |
description |
З метою вдосконалення управління температурними режимами у процесі виробництва полімерних композиційних матеріалів на лініях грануляції досліджено розподіл температури і розроблено математичний опис нестаціонарних процесів, що проходять в проміжку між диском та корпусом дискового екструдера. |
format |
Article |
author |
Шкарапута, Л. Митрохіна, Л. Морозова, І. Тищенко, Л. Даниленко, В. Абросімов, В. |
spellingShingle |
Шкарапута, Л. Митрохіна, Л. Морозова, І. Тищенко, Л. Даниленко, В. Абросімов, В. Температурні режими дискового екструдера |
author_facet |
Шкарапута, Л. Митрохіна, Л. Морозова, І. Тищенко, Л. Даниленко, В. Абросімов, В. |
author_sort |
Шкарапута, Л. |
title |
Температурні режими дискового екструдера |
title_short |
Температурні режими дискового екструдера |
title_full |
Температурні режими дискового екструдера |
title_fullStr |
Температурні режими дискового екструдера |
title_full_unstemmed |
Температурні режими дискового екструдера |
title_sort |
температурні режими дискового екструдера |
publisher |
Інститут бiоорганiчної хiмiї та нафтохiмiї НАН України |
publishDate |
2003 |
url |
http://dspace.nbuv.gov.ua/handle/123456789/3754 |
citation_txt |
Температурні режими дискового екструдера / Л. Шкарапута, Л. Митрохіна, І. Морозова, Л. Тищенко, В. Даниленко, В. Абросімов // Катализ и нефтехимия. — 2003. — № 11. — С. 42-45. — Бібліогр.: 4 назв. — укp. |
work_keys_str_mv |
AT škaraputal temperaturnírežimidiskovogoekstrudera AT mitrohínal temperaturnírežimidiskovogoekstrudera AT morozovaí temperaturnírežimidiskovogoekstrudera AT tiŝenkol temperaturnírežimidiskovogoekstrudera AT danilenkov temperaturnírežimidiskovogoekstrudera AT abrosímovv temperaturnírežimidiskovogoekstrudera |
first_indexed |
2025-07-02T07:00:21Z |
last_indexed |
2025-07-02T07:00:21Z |
_version_ |
1836517552570761216 |
fulltext |
42 Катализ и нефтехимия, 2003, №11
УДК 665.7.038.5:519.876.2 © 2003
Температурні режими дискового екструдера
Л. Шкарапута, Л. Митрохіна, І. Морозова, Л. Тищенко, В. Даниленко, В. Абросімов
Інститут біоорганічної хімії та нафтохімії НАН України,
Україна, 02094 Київ, вул. Мурманська,1, факс: (044) 559-98-00
З метою вдосконалення управління температурними режимами у процесі виробництва полімерних компо-
зиційних матеріалів на лініях грануляції досліджено розподіл температури і розроблено математичний
опис нестаціонарних процесів, що проходять в проміжку між диском та корпусом дискового екструдера.
Виробництво сучасних композиційних полімерних
матеріалів для багатошарових антикорозійних покрит-
тів включає стадію змішування інгредієнтів у дисково-
му екструдері [1]. На якість матеріалів суттєво впли-
вають температурні умови в проміжку між диском і
корпусом екструдера (рис. 1). На жаль, безпосередній
контроль за розподілом температур, зазвичай, відсут-
ній. У промислових екструдерах контролюють лише
температуру корпусу. Виникають запитання, як вимі-
рювана величина співвідноситься з фактично існую-
чою в розплаві, яким є розподіл температур у зазорі,
яким чином одержати найдостовірнішу інформацію
про тепловий режим під час переробки матеріалу і,
найсуттєвіше, як можна здійснити автоматичне керу-
вання процесом. Останнє потребує знання динамічних
характеристик.
Розглянемо розплав багатокомпонентної полімер-
ної композиції під час виготовлення клею-розплаву [2]
в зазорі дискового екструдера, де відбувається диспер-
гування, гомогенізація і остаточна пластифікація мате-
ріалу, як плоску стінку ,0
z
t
y
t
2
2
2
2
=
∂
∂
=
∂
∂ (рис. 1).
Рис. 1. Розрахункова схема дискового екструдера
Рівняння теплопровідності матиме такий вигляд [3]:
λ
−=
∂
∂ q
x
t
2
2
(1)
де t – температура, °С; x – координата, м; q – інтенсив-
ність тепловиділення у розплаві, кВт /м3; λ – коефіцієнт
теплопровідності, Вт/(м⋅град).
Двічі інтегруючи рівняння (1), одержимо
1Cxq
dx
dt
+
λ
−= (2)
01
2
CxC
2
xqt ++⋅
λ
−= (3)
Для визначення констант С0 та С1 розглянемо
граничні умови ІІІ роду [3]:
( )
0x0
0x
tt
dx
dt
=
=
−α=λ− (4)
( )tt
dx
dt
x1
x
−α=λ−
δ=
δ=
(5)
де α та α1 – відповідні коефіцієнти тепловіддачі,
Вт/(м2⋅град) (рис. 1).
З рівняння (2) випливає 1
0x
C
dx
dt
=
=
, а з виразу (3)
00x
Ct =
= . Тоді гранична умова (4) матиме вигляд
.010 tCC =
α
λ
−
У результаті аналогічних перетворень одержали в
явному вигляді вираз і для другої граничної умови (5):
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
−+δ+
δ
⋅
λ
−α=⎟
⎠
⎞
⎜
⎝
⎛ +δ
λ
−λ− 101
2
11 tCC
2
qCq
Визначивши із системи останніх двох рівнянь зна-
чення сталих С0 та С1, знайшли вираз для розрахунку
розподілу температури за товщиною шару матеріалу:
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
⎟
⎠
⎞
⎜
⎝
⎛
δ
−⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
λ
δ
+
α
⎟
⎠
⎞
⎜
⎝
⎛
δ
+
αδ
λ
λ
δ
+
+⎥
⎦
⎤
⎢
⎣
⎡
⎟
⎠
⎞
⎜
⎝
⎛
α
+
λ
+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
α
+
λ
−δ
=
2
1
2
1
1
0
x
K2
1
2
1xKq
1xt1xtKt
(6)
Катализ и нефтехимия, 2003, №11 43
де
1
1
11K
−
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
α
+
λ
δ
+
α
=
Після експериментального визначення аргументів
останнього виразу (екструдер лінії грануляції ЛГП-
60): δ = 6⋅10-3 м; λ = 0,348 Вт / (м⋅град); α = α1 =
= 58 Вт/(м2 град); q = 1044 кВт/м3 було розраховано
розподіл температур в проміжку між диском та кор-
пусом за різних температур корпусу (t0) та диску (t1).
Встановлено, що за рахунок високої інтенсивності
тепловиділення локальні перегріви можуть досягати
багатьох десятків градусів (рис. 2) (вихідні t0 = 100 °С,
t1 = 80 °С). Це підтверджується експериментально: під
час переробки цієї композиції за зазначених умов
спостерігали значне газовиділення і підгоряння мате-
ріалу. Встановлено, що керування тепловим режимом
слід здійснювати не за температурою
корпусу екструдера (зона I), а за пока-
заннями термопари наступної зони II,
яка є менш інерційною і дає вірогід-
нішу інформацію про температуру
розплаву матеріалу.
Рис. 2. Розподіл температур (t) по довжині (х) в зазорі ди-
скового екструдера
Аналізуючи неусталені режими, доцільно перейти
до рівнянь у відхиленнях. Рівняння нестаціонарної
теплопровідності має такий вигляд:
τ∂
Δ∂
=
γ
+
∂
Δ∂ t
c
q
x
t
2
2
а (7)
де а – коефіцієнт температуропровідності,
a = λ/γc, м2/с; γ – густина, кг/м3;
с – теплоємність, кДж /(кг⋅град).
Виконавши перетворення за Лапласом
( ∫
∞
τ−
−
τΔ=
0
p
ii dett , де р – оператор Лапласа),
отримали рівняння теплопровідності у звичай-
них похідних:
з граничними умовами
(8)
.; 1
x
x
1
0
0x
0x
tt
dx
tdtt
x
td −
δ=
−
δ=
−
−
=
−
=
−
Δ=Δ+
Δ
α
λ
−Δ=Δ+
∂
Δ
α
λ
−
Розв’язка рівняння (8) у загальному вигляді має та-
кий вид [4]:
.expexp ⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
−+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
+
γ
=Δ
−
aa
pxCpxC
cp
qt 43
Використовуючи граничні умови, визначили С3 та
С4. Після перетворень одержали вираз, який описує по-
ведінку відхилення температури як функцію характери-
стик клею-розплаву та оператора Лапласа:
Передаточні функції по окремих каналах опису-
ються наступними виразами:
• канал
−−
→ tt 0 (вплив температури корпусу на роз-
поділ температур у розплаві композиції):
• канал
−−
→ tt1 (вплив температури диску на розподіл
температур у клеї-розплаві):
• канал
−−
→ tq (вплив інтенсивності тепловиділення
на розподіл температур у клеї-розплаві):
140
144
148
152
156
0 1 2 3 4 5 6x · 103, м
t, °C
( )
( ) ;
a
pexp
a
p1
a
p1
a
pexp
a
p1
a
p1
a
px1exp
a
p1
a
px1exp
a
p1
pt
ptW
2
1
2
1
2
1
2
1
0
0
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
δ−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λα
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
−−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
δ⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λα
+⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
+
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⎟
⎠
⎞
⎜
⎝
⎛
δ
−δ−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λα
−−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⎟
⎠
⎞
⎜
⎝
⎛
δ
−δ⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λα
+
=
Δ
Δ
=
−
.
a
pexp
a
p1
a
p1
a
pexp
a
p1
a
p1
a
px1exp
a
p1
a
px1exp
a
p1
a
pexp
a
p1
a
p1
a
pexp
a
p1
a
p1
a
pxexp
a
p1
a
pxexp
a
p1
1
p
a
q
tW
2
1
2
1
2
1
2
1
2
1
2
1
q
⎪
⎪
⎭
⎪
⎪
⎬
⎫
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
δ−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λα
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
−−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
δ⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λα
+⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
+
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
+⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λα
−−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λα
+
+
+
⎪
⎪
⎩
⎪
⎪
⎨
⎧
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
δ−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λα
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
−−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
δ⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λα
+⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
+
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
−−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
+
−
λ
=
Δ
Δ
=
−
−
(9)
(10)
( ) ( )
.
a
pexp
a
p1
a
p1
a
pexp
a
p1
a
p1
a
pxexp
a
p1
a
pxexp
a
p1
cp
qt
a
pexp
a
p1
a
p1
a
pexp
a
p1
a
p1
a
pxexp
a
p1
a
pxexp
a
p1
cp
qt
cp
qt
11
11
0
11
1
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
δ−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
−−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
δ⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
+⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
+
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
−δ−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
−−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
−δ⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
γ
−
+
+
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
δ−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
−−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
δ⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
+⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
+
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
−−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
α
λ
+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
γ
−
+
γ
=Δ
−
−
−
λ
−=−
−
−
qtp
dx
td
2
2
а
44 Катализ и нефтехимия, 2003, №1
Було проаналізовано зміну температури у найбільш
”гарячому” місці: х = 5⋅10-3 м при δ = 0,01 м, λ = 0,348
Вт/(м⋅град), с = 2,51 кДж/(кг⋅град), γ = 940 кг/м3, а =
= 1,47⋅10-7 м2/с. За таких умов амплітудно-фазова хара-
ктеристика (АФХ) по каналу
−−
→ tt 0 має такий вигляд:
( )
( ) ( ) ( ) ( )
( ) ( ) ( ) ( )ω−ω−−ωω+
ω−ω−−ωω+
=
=ω
j4340j26101j4340j26101
j2710j26101j2170j26101
0050jW
2
,exp,,exp,
,exp,,exp,
,;
де j = 1− ; ω – частота.
Результати розрахунку АФХ в координатах Im(ω)–
Re(ω) наведено на рис. 3.
Рис. 3. Амплітудно-фазова характеристика
Була розглянута поведінка замкнутої системи з ре-
гулятором (канал t0 → t). Виходячи із співвідношення
Wзам = Wр⋅W(1 + Wр⋅W)-1 було розраховано АФХ за-
мкнутої системи та її дійсну частотну характеристику
Re(ω), рис. 4.
Рис. 4. Залежність дійсної частотної характеристики
Для визначення перехідного процесу залежність
Re(ω) було апроксимоване відрізками прямих Re =
= Re2 + (Re1-Re2)(ω2 - ω)/(ω2 - ω1), де індекси 2 і 1 ви-
значають координати кінця і початку відрізка. В на-
шому випадку слід виділити чотири інтервали, п’ять
характерних частот. Нижче наведені координати ін-
тервалів:
Характерна частота ωi Дійсна частина АФХ Rei(ω)
0 0,333
30 0,076
60 -0,010
100 -0,031
300 -0,013
Зазначимо, що при ω = 300 Re300 < 5 % R0.
Використовуючи інтеграл Фур’є [4], одержали за-
гальний вираз для визначення динаміки зміни
температури:
Для ”характерного” інтервалу, де ω2>ω1 відповідна
складова Δt(τ) буде
Після визначення коефіцієнтів при ( )
ω
ωτsin та ( )ωτsin
одержали вираз для розрахунку перехідного процесу
по каналу t0 → t (рис. 5).
Рис. 5. Перехідний процес
-0,3
-0,2
-0,1
0,0
0,1
-0,1 0 0,1 0,2 0,3 0,4 0,5
Re (ω )jIm
( ω
)
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) .dsinRedsinRedsinRedsinRe2)(t
300
100
43
100
60
32
60
30
21
30
0
10
⎪⎭
⎪
⎬
⎫
ω
ω
ωτ
ω+ω
ω
ωτ
ω+ω
ω
ωτ
ω
⎪⎩
⎪
⎨
⎧
+ω
ω
ωτ
ω
π
=τΔ ∫∫∫∫
( ) ( ) ( ) ( ) ( ) ( )
( )
( )
( ) ( )
( ) ( )
( )
( )
( ) ( )
( ) ( ) .dsinReResinReReRe
dsinReResinReReRe2
dsinRedsinRe2dsinRe2
1
2
122
1
0 12
21
2
12
21
2
0 12
21
2
12
21
2
00
⎪⎭
⎪
⎬
⎫
ω
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
ωτ
ω−ω
−
−
ω
ωτ
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
ω
ω−ω
−
+−
⎪⎩
⎪
⎨
⎧
−ω
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
ωτ
ω−ω
−
−
ω
ωτ
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
ω
ω−ω
−
+
π
=
=
⎪⎭
⎪
⎬
⎫
⎪⎩
⎪
⎨
⎧
ω
ω
ωτ
ω−ϖ
ω
ωτ
ω
π
=ϖ
ω
ωτ
ω
π
∫
∫
∫∫∫
ϖ
=ϖ
ϖ
=ϖ
ϖ
=ϖ
ϖ
=ϖ
ϖ
ϖ
0
0,1
0,2
0,3
0,4
0 10 20 30 40 50 60τ, хв.
Δτ, °C
-0,1
0
0,1
0,2
0,3
0,4
0,5
0,6
0 200 400 600 800 1000
ω
Re (ω )
Катализ и нефтехимия, 2003, №11 45
Запропоновані співвідношення дають змогу розра-
хувати систему автоматичного управління температу-
рними режимами дискового екструдера.
Література
1. Бортников В.Г., Основы технологи переработки
пластических масс, Ленинград, Химия, 1983.
2. Композиція клею-розплаву адгезійного КРА-99,
ТУ У 24.1-03563790-046-2002.
3. Исаченко В.П., Осипова В.А., Сукомел А.С., Теп-
лопередача, Москва, Энергия, 1969.
4. Корн Г., Корн Т., Справочник по математике для
научых работников и инженеров, Москва, Наука, 1974.
Надійшла до редакції 26 травня 2003 р.
Температурные режимы дискового экструдера
Л. Шкарапута, Л. Митрохина, И. Морозова, Л. Тищенко, В. Даниленко, В. Абросимов
Институт биоорганической химии и нефтехимии НАН Украины,
Украина, 02094 Киев, ул. Мурманская,1; факс: (044) 573-25-52
С целью усовершенствования управления температурными режимами в процессе производстве поли-
мерных композиционных материалов на линиях грануляции исследовано распределение температуры
и разработано математическое описание нестационарных процессов, которые проходят в промежутке
между диском и корпусом дискового экструдера.
Temperature conditions of a disk extrude
L. Shkaraputa, L. Mitrohina, I. Morozova, L. Tishchenko, V. Danilenko, V. Abrosimov
Institute of Bioorganic Chemistry and Petrochemistry, National Academy of Sciences of Ukraine,
1, Murmanska Str., Kyiv, 02094 Ukraine, Fax.: (044) 573-25-52
The allocation of temperature has been studied for the purpose of improving the control of temperature condi-
tions while producing the polymeric composites on granulation lines of proceeding and the mathematical speci-
fication of non-stationary processes in a gap between disk and body of disk extruder has been developed.
( ) ( )[ ] ( )[ ]
( ) ( )[ ] ( ) ( )[ ]
( ) ( )[ ] ( )[
( )] ( ) ( )[ ]
( ) ( )[ ] .coscos,
,cos
cos,,
coscos,,
cos,,
⎭
⎬
⎫
τ−τ
τ
⋅
−
−τ−τ−τ−
−τ
τ
⋅
+τ−τ+
+τ−τ
τ
⋅
+τ−τ+
+
⎩
⎨
⎧
−τ
τ
⋅
+τ
π
=τΔ
−
−
−
−
10030010090
100Si300Si040060
10010525060Si100Si02150
30601087230Si60Si1620
1301057830Si33302t
3
3
3
3
|