книги из ГПНТБ / Железнов, Ю. Д. Статистические исследования точности тонколистовой прокатки
.pdfВозможность системы регулирования, необходимо провести ана лиз работы исполнительного механизма, включающего в себя трубопровод и гидроцилиндры.
Для нахождения полосы пропускания исполнительного меха низма надо исследовать динамику гидравлического привода — гидроцилиндров, помещенных между подушками опорных и ра бочих валков.
Дифференциальное уравнение движения сервомотора выво дится из условия равенства расхода рабочей жидкости в силовом цилиндре и в рабочем окне золотника; это условие вытекает из
неразрывности потока рабочей |
жидкости |
[42]. |
|
|
F 1U= |
( 0 |
V ~у( Р п “ |
Рсл)’ |
( 1 5 5 ) |
где F — эффективная |
площадь поршня; |
|
|
|
у — координата поршня; |
|
|
|
|
а — коэффициент |
расхода в рабочем окне золотника; |
|||
/ (t) — закон открытия проходного сечения |
рабочего окна |
|||
золотника; |
|
|
|
|
g — ускорение свободного падения; у — плотность рабочей жидкости; Рп— подводимое давление;
Рcjj — давление на сливе.
При учете массы поршня, действующих на него сил (нагрузки на поршень) и сжимаемости рабочей жидкости в полостях цилиндра процесс взаимодействия жидкости с твердым телом в поршневом исполнительном механизме описывается приведенными ниже уравнениями.
йРпр
Рис. 107. Принципиальная схема системы автоматического регулирования поперечной разнотолщинности
173
Уравнение |
неразрывности |
потока |
жидкости: |
,у dAP -г 2f |
- f - = 2а/ (0 у |
у [Рп - |
Рсл - 2Рс - АР sign /(<)]; (156) |
уравнение равновесия сил, действующих на поршень:
'"-S'+'K-f-' >)+Х‘ЧР-%-=ГЫ>. <‘<*>
где |
ky — коэффициент упругости столба рабочей жидко |
|
сти в полости цилиндра; |
т— приведенная масса поршня;
ф( й ’ — силы сопротивления, действующие на поршень;
АР — перепад давления в рабочих полостях силового цилиндра;
Рс — потери давления на преодоление местных со противлений на пути рабочего окна золотника до поршня [при турбулентном режиме течения
|
Рс = |
^(-gf" )2> гДе |
£ коэффициент пропорцио |
п |
нальности]; |
|
|
dy |
сухого трения |
на поршне. |
|
К sign |
----- сила |
||
Исключая из этих уравнений перепад давлений АР и полагая sign dyldt = sign f (/), получаем уравнение движения поршня гидравлического исполнительного механизма:
kyrn > d*y_ |
■ |
ky_'_d_(b( |
у\ + 2Р |
dy |
|
|
|||
F |
dt3 |
F |
dt |
\ dt ’ |
y) 1 |
dt |
|
|
|
- 2af (0 Y - f |
{я„ - |
Pc, - |
P- |
R ~ 2£ |
- |
|
d2y |
-> |
|
|
|
||||||||
F |
Hi2 |
|
|||||||
|
|
+ |
ф { ч г > |
i/)sign/(o]}. |
|
|
(158) |
||
|
|
|
|
|
|||||
Таким образом, при учете сжимаемости рабочей жидкости движение поршня описывается нелинейным дифференциальным уравнением третьего порядка. Уравнение (158), линеаризованное в возмущениях, исследовано в работе [43].
Учет сжимаемости рабочей жидкости все же не дает возмож ности полностью определить динамическую характеристику ги дравлического исполнительного механизма, так как в уравнении (158) не учитываются волновые процессы в трубопроводах, а жид кость рассматривается как безмассовая пружина. Волновые про цессы в соединительных трубопроводах и полостях механизмов гидроавтоматики существенно влияют на работу этих механизмов, особенно при длинных трубопроводах.
174
Волновые процессы в системе трубопровод— гидроцилиндр описываются [44 ] для одномерного течения невязкой сжимаемой жидкости относительно возмущения по давлению и скорости жидкости при постоянных установившихся давления Р 0 и скоро сти жидкости 1>0 системой
дР |
8яц |
|
dQ_ |
= |
0 ; |
||
дх |
^ |
Q+ 4 |
dt |
|
|||
|
|
Т |
|
|
|
|
(159) |
|
dQ |
FT . |
др |
n |
|
|
|
|
|
|
|
||||
|
дх |
рa2 |
dt |
|
’ |
|
|
где Q = FTv — расход рабочей жидкости; |
|
||||||
Fr — площадь |
поперечного |
сечения трубопровода; |
|||||
ц, р — динамическая вязкость |
|
и |
плотность; |
||||
а — скорость |
звука |
в данной |
среде. |
||||
Вводя обозначения R = вяц//^; |
L = |
|
рIFT; С = Frlpa2, полу |
||||
чим систему дифференциальных уравнений в частных производ
ных, |
описывающих |
|
процесс |
распространения волны давления: |
|||
|
|
|
|
др |
■RQ+L dQ |
= 0, |
|
|
|
|
|
дх |
|
dt |
(160) |
|
|
|
|
|
dQ + С дР |
||
|
|
|
|
|
0 . |
||
|
|
|
|
|
дх |
dt |
|
Начальные условия: |
Pl=o, |
P t=m = Рп\ граничные условия: |
|||||
Р х = |
о |
=Q x = l |
S= |
CЛ^ P г^ i . , |
|
уравнением данной линии |
|
Система |
уравнений (160) является |
||||||
как системы с распределенными параметрами. Эта система ана логична уравнению длинных электрических линий без утечки [10 ].
Здесь R — коэффициент, учитывающий сопротивление трубо провода течению рабочей жидкости из-за трения на стенке трубо провода; L —■коэффициент, учитывающий плотность протекаемой среды; С — коэффициент, учитывающий объем трубопровода.
Для решения этой системы используем интеграл Лапласа, который преобразует уравнение в частных производных в обыкно
венное |
дифференциальное уравнение: |
|||
|
|
dP |
|
RQ + LSQ = 0, |
|
|
dx |
|
|
|
|
|
(161) |
|
|
|
|
dQ |
|
|
|
|
4- CSP = 0, |
|
|
|
|
dx |
|
где P |
и |
Q — изображения |
соответствующих функций P (x , t) |
|
|
|
и Q (x , 0, |
t . |
e. |
|
|
P (x, |
S) = j p (x, t) e~st dt, |
|
|
|
|
|
о |
где 5 = |
a + jb — комплексное число, оператор преобразования. |
|||
175
Из системы уравнений (161) исключим Q, тогда |
|
||
d-J> |
_ |
( 1 6 2 ) |
|
- f i f l - C S ( L S + R ) P ~ 0 . |
|||
Характеристическое |
уравнение имеет следующий |
вид: |
|
где |
fe2— v2 = о, |
|
|
|
|
|
|
|
v = VCS (LS + R). |
|
|
Уравнение (162) решается так: |
|
||
Р (х, S) = |
Лх ch vx + А2sh vx. |
(163) |
|
Аналогично получим для расхода: |
|
||
Q(x, S) |
■ |
(/4Хsh vx + А2ch vx), |
(164) |
где |
|
R + L C |
|
|
- V : CS |
|
|
Произвольные константы А х и А 2 находятся с помощью гра ничных условий и уравнений (163) и (164). Окончательно получим:
р
{ ' ’ ~
Q {х, S) =
ch у (х — /) -|- ЗСр sh у (/ — х) |
- |
|
|
ch уI + SCp sh vI |
Нп' |
Р |
|
(165) |
_ sh у (/ — х) + SCp ch у (/ — х) |
||
р |
ch у/ -f- SCp sh уI |
|
Пренебрегая С в (165), так как он мал по отношению к емкости трубопровода, и считая изменение давления на входе трубопро вода скачкообразным, получим:
|
Р (х, S) |
Рп |
ch у (/ — х) |
|
(166) |
|
S |
ch у/ |
|
||
|
|
|
|
||
Для получения оригинала комплексной функции Р (х , S) |
|||||
используем вторую |
теорему |
разложения Хевисайда |
[10]: |
||
Р (х, t) |
У{0) |
|
U(Sk) |
eV, |
(167) |
|
W (0) |
|
|
|
|
где U = Pnchv (l — х), W = |
chx. |
|
|
||
Найдем значения S, превращающие в нуль знаменатель урав |
|||||
нения (167) |
, . . |
|
, ,2 k + l |
|
|
|
|
|
|
||
|
ch у/ = 0, |
у* / = / — g— я - |
|
|
|
где k — любое отрицательное, положительное целое число или нуль.
176
Зная, |
что v |
|
|/ |
(LS + R) CS, |
найдем |
|
|
|
|||||
|
|
|
|
|
LCS2+ |
RCS — v2 = О, |
|
|
_ Y -1 |
||||
|
|
S = — 1 JL |
Л Г R 2 |
|
LC |
( 2 k + i |
j |
||||||
|
|
|
2 ' |
L |
V |
4L2 |
|
\ |
2/ |
’ |
|||
|
|
|
|
|
|
Si = — a ± |
Pi; |
|
|
|
|
||
при S = |
0 |
U{0) = |
Pn, |
W(0) = 1 |
И |
|
|
|
|
|
|||
|
|
S = |
SkU (Sk) = cos ( |
|
1 ■ 1 |
я j |
, |
|
|||||
|
- |
d |
t„ |
, |
, |
, , dv |
, |
, |
, |
2SfcLC -{- RC |
-- |
||
|
7K— It7 |
(S/{) = |
/ sh vl -Jcr ~ 1sh vl ---- - |
|
|
||||||||
|
dSk |
|
|
|
dS |
|
2 V (L C + R ) CS k |
||||||
= i (— i) * |
2 t c ( s , + |
|
, |
, чА, 7-C ( — a - | - p * . - f - a ) |
|||||||||
|
|
2\k |
|
|
' |
|
2 * 4 - 1 |
Я |
|||||
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
2/ |
||
Итак |
|
|
= |
|
|
|
|
|
|
|
|
|
|
Подставив 7/ |
(0), |
W (0), |
U (Sk), |
W (Sk) |
в |
уравнение (167), |
|||||||
находим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P (x, |
t) = P„ |
1 + |
|
я |
e-a t |
X |
|
|
|
|
|
|
|
|
|
|
|
l2LC |
|
|
|
|
|
|
|
CD |
|
|
|
' 2 * + l |
/ — x |
»)<2t+ l> |
|
||||
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
COS |
' |
|
/ |
(168) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
2Pi (— a ± |
Pi) |
|
|
|||
|
|
|
|
|
|
|
|
|
|
||||
И Л И |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p ( X , i) = P„ 1 + |
|
|
e -at X |
|
|
||||
00
X s
k = 0
cos |
'2k + 1 |
l — x |
я) (2*+ 1) |
|
|
v 2 |
' / |
,±M |
(169) |
||
( - 1 Г |
2P* (— a ; |
|
|
|
|
12 JO- Д . Ж е л е з н о в |
177 |
С помощью формулы |
Эйлера |
для гиперболических |
синусов |
|
и косинусов |
получим |
|
|
|
|
Р (х, () = Рп |
Я |
|
|
|
Же e~ at X |
|
||
оо |
(2k + I \ ( ^ |
“h 1 l — |
X Jl^ |
|
X |
|
|||
K - * |
(-^-sh + ch P*^ |
• (170) |
||
иk=0 |
|
|
||
Последнее выражение является решением уравнения распро странения волны давления в трубопроводе.
Так как нас интересует конец трубопровода — вход гидро цилиндра, то, подставив в выражение (170) х = I, получим:
Р(1, 0 - Рп |
• (171) |
|
А=0 |
Для расчета переходного процесса были взяты k — 10, Р = 1 |
|
и формула |
(171) была запрограммирована на ЭЦВМ «Наири» |
по исходным данным гидропривода. По полученным результатам построены кривые переходного процесса (рис. 108), которые показывают, что в начале процесса имеется время транспортного запаздывания. Это связано с тем, что при изменении давления на входе вначале происходит сжатие жидкости, после чего пере дается волна давления. Для гидравлических систем время сжа тия очень мало в отличие от пневматических систем, где запазды вание составляет 1— 1,5 с.
Рис. 108. Переходные процессы движения поршня гидроцилиндра (а) и измене* ния давления на конце трубопровода (б):
------- расчетные кривые;----------------- |
аппроксимированные кривые |
178
Рис. 109. Структурная схема САР поперечной разнотолщинности
На рис. 108 пунктирной линией показаны аппроксимирован ные кривые, соответствующие инерционным звеньям первого порядка. Передаточная функция для гидроцилиндра
1 |
(172) |
Ф(Р) = 1 + О.ООЗЗр |
Для расчета гидропривода были использованы следующие
данные: |
150 |
кгс/см2— номинальное давление; |
|
р о = |
|||
D = |
100 |
мм — диаметр |
гидроцилиндра; |
Н = |
120 |
мм — высота |
гидроцилиндра; |
/ = |
50 |
мм — длина трубопровода; |
|
d — 18 мм — диаметр трубопровода.
Рабочей жидкостью служило масло индустриальное-20, обла дающее при 20° С плотностью (р) 0,88— 0,89 г/см3.
Трубопровод является звеном с распределенными параме трами. Как известно из теории, системы с распределенными параметрами можно аппроксимировать звеном с запаздыванием.
Тогда передаточная функция трубопровода |
будет: |
|||
WT(р) = |
1 |
е-0,002Р> |
(173) |
|
1 + 0,012р |
||||
|
||||
Как видно из структурной схемы САРТ (рис. 109), исполни тельный механизм состоит из двух звеньев, соединенных последо вательно. Передаточная функция исполнительного органа запи шется так:
(р) — ^гц (р) WT(р )-- |
1 |
-0,002р |
(174) |
(1 + 0,033р)(1 + 0,012р) |
|
||
Частотная характеристика |
исполнительного механизма: |
|
|
W(/fill = ______________ !______________ p~i (0,02-(-arctg О.ОЗЗш+arctg 0,012w)_
/( 1 + 0,001 1oj2) (1 + 0,000144co2)
(175)
12* |
179 |
Обозначим
______________ 1______________ |
: A ( I D ) . |
(176) |
|
\ 1 + 0,001 ко2) (1 + 0,000144со2) |
|||
|
|
По формуле (176) рассчитана и построена амплитудно-частот
ная характеристика исполнительного |
механизма (см. рис. 103): |
||||||
(о .............................. |
0 |
4 |
8 |
12 |
16 |
20 |
30 |
А (ы ).......................... |
1 |
0,98 |
0,92 |
0,86 |
0,80 |
0,66 |
0,44 |
Частота среза исполнительного механизма на уровне 0,7 составляет /ср = со/2я = 3,1 Гц.
Таким образом, исполнительный механизм пропускает с дву кратным запасом все возмущения поперечной разнотолщинности.
Для измерения этих возмущений, как видно из рис. 103, тре буется микрометр с постоянной времени Т = 0,1 с.
на |
Для более точного измерения на более высоком уровне, т. е. |
|
уровне 0 ,8Л (0), требуется микрометр с постоянной времени |
||
Т = 0,08 с, |
а на уровне 0,9Л (0) микрометр с постоянной времени |
|
не |
больше |
Т = 0,05 с. |
| Требования к точности измерителей поперечной разнотолщин ности зависят от того, для чего предназначена информация о поперечной разнотолщинности. Так, если нужно измерить попе речную разнотолщинность на входе и выходе из клети и в даль нейшем использовать эту информацию для контроля планшетности, то требования к точности измерителей должны быть очень высокими. Относительная разница вытяжек, которая приводит к искажению формы (волнистости, коробоватости и серповидности), исчисляется сотыми долями процента. Этими же величинами должна определяться относительная разница обжатий по участкам ширины, из-за которой возникает неравномерность продольной деформации. Отсюда следует, что контроль за планшетностью через измерение поперечного профиля подразумевает измерение поперечной разнотолщинности с точностью до сотых долей про цента; однако при нынешнем уровне развития техники измерения это невозможно (точность современных толщиномеров не превы шает 1%). Кроме того, поскольку самые последние исследования указывают на возможность поперечного перемещения металла даже при тонколистовой холодной прокатке, представляется сом нительной сама идея использовать соотношение профиля попе речного сечения полосы до и после прокатки для оценки неравно мерности продольной деформации.
Если же необходимо проконтролировать поперечную разно толщинность горячекатаной полосы, предназначенной в последу ющем для холодной прокатки, при оптимальных значениях чече вицеобразной горячекатаной полосы 0,02— 0,05 мм достаточной была бы точность измерения до 0,005 мм. Это, однако, означает, что по отношению к средней толщине полосы эта точность должна
180
составить 1% и практически является Достижимой в том случае, когда для измерения поперечной разнотолщинности используется один толщиномер с периодическим сканированием измерений по ширине.
Возможность измерять поперечную разнотолщинность с по мощью отдельно стоящих по ширине трех толщиномеров, име ющих точность более 1 %, является весьма сомнительной, по скольку погрешность отдельного толщиномера будет соизмерима с заданной точностью измерения и, кроме того, очень трудно до биться стабильной и абсолютно идентичной работы в течение дли тельного времени всех приборов. Все же не исключено, что при включении толщиномеров, расположенных у края и у середины в дифференциальной схеме, требования к итоговому сигналу, пропорциональному значению поперечной разнотолщинности, и в отношении измерения могут быть существенно снижены.
До настоящего времени вопрос об измерении поперечной раз нотолщинности в процессе холодной прокатки не ставился с до статочной остротой. Если же потребовались бы такие измерения для каких-то специфических видов продукции, то, очевидно, предпочтительной оказалась бы дифференциальная схема, при которой отдельные толщиномеры располагаются у краев и по середине полосы. Использование одного сканирующего по тол щине толщиномера при достаточно высоких скоростях холодной прокатки было бы весьма затруднительным.
3. Статистичесние характеристики неровностей тонколистового проката
а. Зависимость формы полос от продольной и поперечной разнотолщинности подката
Форма полосы на выходе из стана определяется рядом техно логических параметров: профилировкой валков, величиной дав ления металла на валки и усилия гидроизгиба валков, распре делением условий трения по длине бочки валков, профилем и формой полосы на входе в стан. В процессе прокатки технологи ческие параметры непрерывно изменяются, что может привести к изменению толщины и формы полосы.
Статистический анализ частот и амплитуд колебаний регу лируемого параметра произведен на основании опытных данных о характере поведения показателей геометрии в процессе прокатки. Экспериментальные данные могут быть получены:
а) осциллографированием давления или толщины с пересчетом их изменений на бр/р либо осциллографированием характери стики неравномерности вытяжек — распределения удельных на тяжений ба2 = бр/р£;
б) непосредственным контролем толщины и формы полосы при перемотках рулона или нарезанных листов.
181
Учитывая, что изменение толщины полосы приводит к колеба ниям давления прокатки и, следовательно, профиля прокатной щели и формы полосы, и что толщина является одним из опреде ляющих параметров холоднокатаного листа, провели экспери ментальное исследование зависимости формы полосы от ее про дольной и поперечной разнотолщинности. Исследовали прокатку полосы размером 0,5X 1030 мм из стали 08кп на дрессировочном стане. Рулон прокатывали при выключенной системе гидроизгиба валков с волной по всей длине полосы. После дрессировки рулон
нарезали на |
листы длиной |
2 |
м. Вручную |
замеряли |
толщину, |
|
поперечную |
разнотолщинность, |
амплитуду |
А (мм) |
и длину |
||
волны X (мм). |
|
|
|
|
|
|
По формулам [44 ] |
|
|
|
|
|
|
|
А/ |
|
/ |
яА \2 |
|
(177) |
|
1 |
~ |
\ |
2к ) ’ |
|
|
|
|
|
||||
|
да = — ^ - Е |
|
(178) |
|||
рассчитывали относительную разность длин и разницу удельных натяжений между серединой и краем полосы. На рис. 110 приве дены накопительные кривые расчетной разности удельных натя жений одного из участков рулона. Экспериментальные данные обрабатывали методами математической статистики. Вероятност
ные |
характеристики |
рассчитывали |
по |
формулам П 1 ), |
(20), |
|
(21), |
(48). |
расчета приведены |
в |
табл. 47 (rAaJl = |
0,08; |
|
Результаты |
||||||
^ест, б* |
== 0,12; |
/■§/,, * = |
0,05). |
|
|
|
f
Рис. ПО. Накопительные кривые разности |
удельного |
натяжения |
||||||
и |
относительной разности Др/р холоднокатаных полос: |
|
||||||
а |
— сортамент, мм; |
1 — 1,0Х 1250; |
2 — 1,5Х 1200; |
3 — 2,0Х 1000; |
||||
замеры А вручную |
и |
пересчет на |
Др/р; |
6 |
— непосредственное |
|||
осциллографирование |
параметра а; |
сортамент, |
мм; |
1 |
— 2Х 1420; |
|||
2 |
— 0,8Х 1030 |
|
|
|
|
|
|
|
182
