Температурні режими дискового екструдера

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

Повний опис

Збережено в:
Бібліографічні деталі
Дата: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 Ukraine
id 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