CC BY
57
4. Eisenbarth E., Meyle J., Nachtigall W., Breme J. Influence of the surface structure of titanium materials on the adhesion of fibroblasts // Biomaterials. -1996. -V. 17. - P. 1399-1403.
5. Понтер В.Э., Дамбаев Г.Ц., Сысолятин П.Г и др. Медицинские материалы и имплантаты с памятью формы. — Томск: Изд-во Том. ун-та, 1998. -487 с.
6. Buchanan R.A., Rigney E.D., Williams S.M. Ion implantation of surgical Ti-6Al-4V for improved resistance to wear-accelerated corrosion // J. Biomed. Hater. Res. -1990. - № 21. - P. 355.
7. Wintermantel E., Mathey M., Mayer J. Biomaterials for Interbody Fusion. Final Report 1995-1999. Swiss Priority Program on Materials Research. Project 4.2B. - P. 68-69.
8. Hiroshi Y. Strength of biological materials / Ed. by F. Gaynor Evang. - Williams and Wilkins Company, Baltimore, 1970. - 297 p.
9. Collier J.P., Surprenant V.A., Jengen R.E. etc. Corrosion between the components of modular femoral hip prostheses // J. Bone St. Surg. -1992. -V. 7413. - P. 511-518.
10. Псахье С.Г., Коростелев С.Ю., Смолин А.Ю. и др. Метод подвижных клеточных автоматов как инструмент физической ме-зомеханики материалов // Физическая мезомеханика. — 1998. -Т. 1. — № 1. —С. 95—100.
11. Psakhie S.G., Horie Y., Korostelev S.Yu. etc. Movable cellular automata method as a new technique to simulate powder metallurgy materials // Proc. Intern. Conf. on Deformation and Fracture in Structural PM Materials, Stara Lesna, Slovakia, 1996. — P. 210—220.
12. Psakhie S.G., Horie Y., Ostermeyer G.P. etc. Movable cellular automata method for simulating materials with mesostructure // Theor. and Appl. Fract. Mech. —2001. —№ 1—3. —P. 311—333.
13. Физические величины. Справочник / Под ред. И.С. Григорьева, Е.З. Мейлихова. — М.: Энергоатомиздат, 1991. —1232 с.
14. Глазер Р. Очерк основ биомеханики. — М.: Мир, 1988. — 128 с.
УДК 621.317.727.1
РАСЧЕТ ПОГРЕШНОСТЕЙ МНОГОДЕКАДНОГО ИНДУКТИВНОГО ДЕЛИТЕЛЯ НАПРЯЖЕНИЯ
В.Л. Ким
Томский политехнический университет E-mail: sov@camsam.tpu.ru
Предложен аналитический способ расчета случайной погрешности коэффициента передачи многодекадного индуктивного делителя напряжения для области верхних частот с учетом стохастических свойств жгутов. Приведены интервалы погрешностей при доверительной вероятности 0,997.
Индуктивные делители напряжения (ИДН) относятся к одной из функциональных основных групп средств измерений - меры и цифро-аналоговые преобразователи, в том числе многозначные меры [1]. Для ИДН как многозначной меры отношения переменных напряжений в комплексе нормируемых метрологических характеристик важнейшей является динамическая характеристика.
Во многих измерительных системах ИДН работает при гармонических воздействиях. Для установившегося режима, где ИДН представляет собой линейное динамическое звено, наиболее приемлема амплитудно-частотная характеристика (АЧХ) делителя. Последняя есть зависимость модуля коэффициента передачи (Кп=ивш/Ц,х, где - выходное, а Цвх - входное напряжение) от круговой частоты а входного сигнала. Для генеральной совокупности ИДН, объединяемых понятием "тип ИДН", нормируют номинальную АЧХ и допускаемые отклонения от неё [1].
В общем случае аналитическое представление коэффициента передачи ИДН имеет вид [2]
кп( а) = кн[1+7» +
где Кн - номинальный коэффициент передачи; 7а(а) и 7ф(а) - амплитудная и фазовая составляющие относительной погрешности.
При пренебрежении 7ф(а) ввиду малости фазовых сдвигов [3] погрешность коэффициента передачи определяется вещественной составляющей, т.е. 7а(®).
В работе [4] проведен расчет амплитудной погрешности однодекадного жгутового ИДН, являющейся для указанного типа делителей случайной величиной. Однако наибольшее применение находят многокаскадные ИДН, поэтому определение их погрешностей является весьма актуальной задачей.
Цель данной статьи - получить аналитическое выражение амплитудной погрешности двухдекад-ного ИДН с учетом стохастических свойств жгутов.
Среди многокаскадных делителей значительный интерес представляют многодекадные ИДН, в которых каскады выполняются путем намотки двух делительных (декадных) обмоток на одном (общем) ферромагнитном тороидальном сердечнике ТУ1 (рис. 1). Отводы первой декадной обмотки ¿1 и второй ¿2 коммутируются многопозиционными переключателями 51 и £2. При таком способе выполнения и соединения декад необходимо, чтобы число витков обмотки ¿1 в десять раз превышало число витков обмотки ¿2. Очевидно, что из-за конструктивных усложнений размещение на общем сердечнике более двух декад нерационально.
&— и™
ТУ1
И 10
ивх
г
>2 10
_9_
_8_
_7_
_6_
_5_
_4_
_3_
_2_
_1_
0
82
ю
1.8, Ф 4,5 О
1-з, =Р12,5С1
1*, Ф 18,5 С]
I.,, ф 22,5 С,
10
1-32 -г 4,5 Сг
-Б 2
=Р 12,5 0.7
=Р 24,5 С,
1+1
1_з, =Р 24,5 О
1*1 =Р 22,5 С,
и, Ф 18,5 С,
Ц, "г12,5С1
""° 1 Ф 4,5 О
1-32 -Г" 18,5Сг
1_32 Ф 22,5 С2
1-32 Ф 24,5 Сг
1_зг Ф 24,5 Сг
1_зг Ф 22,5 Сг
Ьзг ~г 18,5 Сг
Ьвг -р 12,5 Сг
° 1
Ф 4,5 Сг
ивых
Рис. 1. Принципиальная схема двухдекадного ИДН
При преобразовании принципиальной схемы в эквивалентную (рис. 2) учтем, что из-за общего магнитного потока ЭДС второй обмотки равна напряжению секции первой обмотки, т.е. 0,Щх.
Расчет погрешностей двухдекадного ИДН для нормальных условий применения в областях верхних частот (20...200 кГц) проведем при следующих допущениях:
1. индуктивности и активные сопротивления соединительных проводов и коммутационных элементов равны нулю;
2. емкостные проводимости между обмотками и обмотками и сердечником равны нулю;
3. импеданс нагрузки
4. ЭДС секций в каждой делительной обмотке одинаковы.
Тогда в расчетной схеме (рис. 2) учитываем только индуктивности рассеяния секций первой ЬА и второй ЬЛ декад, а также межпроводные емкости жгутов С1, С2. Длина жгута второй обмотки в десять раз меньше длины жгута первой обмотки. При одинаковом диаметре проводов жгутов параметры второй декады также будут меньше в десять раз.
Рис. 2. Упрощенная эквивалентная схема ИДН
Следовательно, малая входная емкость второй декады, определяющая входной импеданс этой декады в области верхних частот, практически не оказывает существенного влияния на точность первой декады. В этом случае полагаем, что первая декада работает в ненагруженном режиме. Тогда для расчета напряжений на отводах первой декады можно использовать формулу [4]
Ь2 '
Ь1' + 12,'
1 + а2 Ь1,
Ь, + Ь2,
V 2,1 у
^ а2Ь Ь2.
— 1 Си + _ '•'_2-' • с2
к, + Ь2'
и„,
где /=(1,10) - номер отвода; Ьи, Ь2/ - эквивалентные индуктивности рассеяния и Си, С2/ - эквивалентные емкости секций, включенных соответственно между отводами /, 10 и 0, /.
Эквивалентные параметры в свою очередь определяются в виде
Ь1,' = X Ь.1т • Ь2, = XЬ 1т •
т='+1 т =1
10 ' С1,' = (10 -') - X С 1т • С2, =Г2 X Ст
(1)
где Ь1т и С1т - индуктивность рассеяния и шунтирующая емкость т-й секции; т=(1,10)).
Заметим, что параметры второй декады определяются по (1) заменой ЬАт-^-Ьйт, С1т^С2т, г^и, где «=(1,10) - номер отвода второй декады.
Схемы декад идентичны, поэтому выходное напряжение ИДН, как напряжение на отводе и второй декады относительно отвода 0 первой декады, определяется в виде
и„ттх = и„„ п = и„„, + ивых п =
ВЫГХ ВЫГХ 1,П ВЫГХ 1 вых п
= ив.
¿1,1 + 4,1-
1 + ю Ь
1,11 т + т 1 + 2,1
ю2 Ь к
-1 |х
С + Ь ' С21
т2,
4,1 + Ь2,1
1 + ю Ь
4,1 + Ь2,1
1 |х
уС + ю 4,1+1 4,1+1 . С
1, 1+1 + Т л. Т 2,1+1
1, 2,
Ь1,1 + Ь2,1
1 + ю2 Ьи
1 + ю Ь
2,
1,11 т + т
1, 2,
ю2 Ь к
- 1|х
уС I " 2,'' . С 1, 2,
4, п + 4,
Ь2,п , ^ , ю Ь1,пЬ2,п С
. С2,п
4, п + Ь
2,п
- 1| С1, п +
4, п + Ь
2,п
Из этого выражения после преобразований получаем формулу амплитудной погрешности двух-декадного делителя
у. = —ю2 у
< 1,п
(9 — 1 )2(10 — 1)-2 + 4,1+1 Ь2,,14, и
(4,1 + Ь2,1)[ Ь2,1 '(4п + Ь2,п )" (¿1,1 + "¿2, 1)]
уС1,1+1 +ю 2у
(1 + 1)2Г 2 г Ь Ь + Ь к ь Н, 1+12,1+12, п
(4, + Ь2,1)[ Ь2,1 '(4п + Ь2,п )" Ю,^ (¿1,1 + "¿2, 1)]
уС2,1+1 - ю 2у
(10 -1) -2 ¿2 Ь Ь + 1,1 2,1 1, п ^ "1 кЬ2,к, п
(4, + Ь2,1)[ Ь2,1 '(4п + Ь2,п )" Ю,^ (¿1,1 + "¿2, 1)]
уС.+1 - ю' У у
¿2 ¿2 (Ь - + 1,п 2, 1,1 -¿и)
10(4, п + 4, п)[ 4, Л4, п + 4, п) + 0,14, п (4, 1 + 4, ,)]
уС1,п +ю2 у
_А,„4„ (4,1 + 4,1)_
' 10(4, п + 4, п)[ 4, Л4, п + 4 п) + 0,14, п (4, 1 + к 1)]
у С...
(2)
и индуктивностей рассеяния влияет качество изготовления жгута. Несмотря на тщательное скручивание проводов в жгуте, не удается достичь полного выравнивания указанных выше параметров. Разброс параметров жгута и случайная последовательность соединения проводов делительной обмотки обуславливают случайный характер амплитудной погрешности. Заметим, что параметры жгутов отличаются высокой временной и температурной стабильностью. Поэтому полагаем амплитудную погрешность не случайным процессом, а многомерной случайной величиной, зависящей от двадцати индуктивностей рассеяния и двадцати шунтирующих емкостей секций.
Расчет амплитудной погрешности практически не выполним, так как требуется знание совместной функции распределения сорока случайных величин. Очевидно, что решение задачи возможно лишь при некоторых упрощающих допущениях. Воспользуемся, например, критерием существенности случайной составляющей погрешности [1].
Так амплитудную погрешность (2), обусловленную разбросом параметра секции полагаем пренебрежимо малой, если среднее квадратическое отклонение (СКО) параметра не превышает 10 % среднего значения.
В работе [5] приведены результаты исследований различных жгутов, показывающие, что СКО индуктивностей рассеяния не превышает 6 %, тогда как СКО межпроводных емкостей достигает 20 %. Эти данные подтверждены измерениями параметров более 100 жгутов, проведенными автором [4].
Тогда, если не учитывать влияние индуктивности рассеяния, случайная погрешность коэффициента передачи ИДН будет зависеть только от статистических свойств емкости между двумя проводами многопроводного жгута.
В выражении (2) суммирование случайных величин необходимо проводить с учетом степени их коррелированности. Емкости разных жгутов считаем независимыми случайными величинами. В работе [4] показано, что емкости в жгуте являются зависимыми случайными величинами с сильной обратной корреляцией, т.е. с коэффициентами корреляции Л[С(.)С(.)]=-1.
С учетом этого из (2) получаем следующее выражение для дисперсии амплитудной погрешности
Ъ[у п ]=А Я[СУ+1]+Д Я[С,1+1]+
Амплитудная погрешность, как следует из (1) и (2), зависит от индуктивностей рассеяния секций и межпроводных емкостей жгутов. Как показано в работах [2, 5] на значения межпроводных емкостей
+а2Щ.СМ] -2АД ^С^МС^ 1+1] -
,1+1
МС
1+1J
-2 Ва,р[С2, ,|/)|С. ,| + ЗДС^ ] + (3)
+т2/[С2п ] - 2 ]/[С2, п ],
где
А =-ю2 •
(9-1)2(10-1)-2 ¿2,Л,А п + 421+1 ¿2,1+14,п
(4,1 + 4,)[ 4,1 (4п + 4,п)+0,14,п (к + 4,)]г
+
у
В, = т ■
й = -®2' Б1 = -т2 ■ Т = т2 ■
(г +1)21 2/.. I:. I, „ + Ь1; +Ь\ , +Ь2 п (Ц. + ¿2, )[/г (/п + Ь2п ) + 0 ,1/2,п (А,г + ¿2, )] ' (10 - г)-2Цг/2гАпП + Г2А, А, А,п
(А,1 + ¿2,г )[¿2 ,г (А„ + ¿2п ) + 0 ,1/п С¿1., + ¿2, )] '
_¿1, ,¿2,, (А,г+ ¿2, г)_
10(^,п + 4,П )[¿2, (¿,,п + ^П ) + 0,1^,п (А,г + ¿г )] ;
А,пА,п (А,г + ¿2,,)
10(А,п + ¿2,п )[¿2, (¿1,п + ¿2п ) + 0,!, (¿1,, + ¿2, )]
ДСи+1], ДСМ], дсу, ДС,„] -дисперсии
величин С1>/+1, Су+1, С+1, С1,„, С2,„.
Далее определяем дисперсии емкостей ДСН]. Пусть для первой декады среднее значение межпроводных емкостей равно С1, а дисперсия - Б[С1] для второй декады - С2 и Б[ С2] соответственно. Тогда по методике [4] можно рассчитать дисперсии эквивалентных емкостей для отводов г, п декад. Например, для отвода г=2 из (1) при г:=г+1 находим
С = С =
1
= —(2 ■ 22,5 + 2 ■ 24,5 +18,5 +12,5 + 4,5) С1 = 2,64 С1;
С2, ,+1 = С23 = 32 (4, 5 +12 , 5 + 18, 5)С = 3 , 94СР С.+1 = С3 = 18,5С1. Дисперсии этих емкостей
£>[С13] = Д[2,64С1] = 7 Б[С1]; £>[С2,3] = Б[3,94С1] = 16 Б[С1]; В[С3] = Д[18,5С1] = 342 ЩС] Для п=1 второй декады из (1) при г:=1 получаем
С1,п = С1,1 = 912 х
х( 2 -12,5 + 2 18,5 + 2 ■ 22,5 + 2 ■ 24,5 + 4,5 ) С2 = 2 С2;
2 12
ДСЦ] = £[2С2] = 4 ДС2]; £>[С2Д] = £>[4,5С2] = 20 ЩС2].
Результаты расчетов для других отводов приведены в табл. 1, 2.
Таблица 1. Дисперсии эквивалентных емкостей первой декады
С = С =— 4 5С = 4 5С •
2,п 2,1 ..2 > 2
/ о[си ОД«] ОД+1]
0 40[д 20ОД 20ОД]
1 5 ОД 18 ОД] 156ОД]
2 7ОД] 16ОД] 3420[С1]
3 9ОД] 13ОД] 506ОД]
4 11 од] 11 ОД] 600ОД]
5 13ОД] 9ОД] 600ОД]
6 16ОД] 7ОД] 506ОД]
7 18ОД] 5ОД] 342ОД]
8 20ОД] 4ОД] 156ОД]
9 0 3ОД] 20ОД]
случаев: для первой декады (К,„>0,10) и второй декады (К,„<0,10). Примем также, что ХЙ=0,1ХЙ, С2=0,1С1. Тогда, например, для коэффициента передачи К,„=0,21, т. е. г=2, п=1 из (2), (3) находим характеристики основной погрешности - математическое ожидание и дисперсию
У2,1$ =-7,58, ДГ2д] = 760К 28 2,
где 8 =т2¿$1С1,V = •У/ДС1]/С1 - коэффициент вариации.
Таблица 2. Дисперсии эквивалентных емкостей второй декады
п ОДп] ОДп]
1 4ОД] 20ОД]
2 5ОД] 18 0[ С2]
3 7ОД] 160[С2]
4 9 ОД] 13ОД]
5 11ОД] ЦОД]
6 13ОД] 9 ОД]
7 16ОД] 7ОД]
8 18ОД] 5ОД]
9 20ОД] 4ОД]
Числовые характеристики случайной погрешности двухдекадного ИДН определяем для двух
Как показано в [4, 5], погрешность коэффициента передачи ИДН распределена по нормальному закону. Тогда при доверительной вероятности Рд = 0,997 получаем интервальную оценку амплитудной погрешности
у21 =-(7,5 ±16,5)8.
При расчете первой декады особое внимание уделяем точкам диапазона преобразования (табл. 3), находящимся слева и справа от номинальных значений коэффициента передачи первой декады на расстоянии минимального шага квантования (цены деления) двухдекадного ИДН, т.е. ±0,01.
Проведем анализ данных табл. 3. Во-первых, в точках квантования первой декады, т.е. К,„ = 0,10, 0,20, 0,30, ..., 0,90 границы интервальных оценок погрешностей совпадают с интервалами погрешностей на отводах однодекадного ИДН [4]. Во-вторых, монотонный характер зависимости систематической погрешности у,^ от коэффициента передачи К,„ не нарушается в точках квантования. В-третьих, изменение разброса амплитудной погрешности в пределах каждого шага квантования происходит скачкообразно, достигая максимального значения в начале и убывая в конце кванта.
Первое обстоятельство подтверждает правомочность предлагаемого подхода к определению случайной погрешности многодекадного ИДН, а второе объясняется меньшим весом второй декады в коэффициенте передачи К„. Однако, вторая декада, как и первая, имея наибольший разброс погрешностей на своих нижних отводах, трансформирует его в общую погрешность двухдекадного делителя.
Очевидно, что и первая декада влияет на вторую. При установке К,„<0,10 вторая декада подключается к нижним отводам 0 и 1 первой декады. Тогда для расчета погрешностей второй декады формулу (2) преобразуем к виду
Таблица 3. Характеристики случайной погрешности двухде-кадного ИДН
Из анализа данных табл. 4 следует важный вывод: интервалы амплитудных погрешностей на всех отводах второй декады одинаковы и полностью определяются характеристиками погрешности первого отвода первой декады, т.е. при К,п = 0,10.
Следует отметить практическую значимость результата расчета второй декады на верхних частотах: повышение точности многодекадных ИДН при малых уровнях выходного напряжения возможно путем уменьшения погрешности на первом отводе первой (старшей) декады, например, емкостной коррекцией, применением новых типов делителей, таких как ИДН с симметрирующей обмоткой [6, 7].
Таблица 4. Характеристики случайной погрешности ИДН при Кп< 0,10
Коэффици- Математи- Дисперсия, Интервал, x„
ент переда- ческое DU
чи, к Ожидание, ytris
0,10 -12,05 400 V252 -(12,0±12,0)5
0,11 -11,25 900 V 252 -(11,2±18,0)5
0,19 -8,15 470 V 2S2 -(8,1±13,0)5
0,20 -8,05 460 V 252 -(8,0±12,8)5
0,21 -7,55 760 V 2S2 -(7,5±16,5)5
0,29 -4,95 460 V 2S2 -(4,9±12,8)5
0,30 -4,75 440 V 252 -(4,7±12,6)5
0,31 -4,35 620 V 2S2 -(4,3±15,0)5
0,39 -2,25 380 V 2S2 -(2,2±11,7)5
0,40 -2,05 370 V 252 -(2,0±11,5)5
0,41 -1,85 470 V 252 -(1,8±13,0)5
0,49 -0,15 280 V 2S2 -(0,1±10,0)5
0,50 0 270 V252 ±9,85
0,51 0,15 310 V 252 (0,1±10,6)5
0,59 1,25 170 V 2S2 (1,2±7,8)5
0,60 1,35 165 V252 (1,3±7,7)5
0,61 1,45 180 V 2S2 (1,4±8,0)5
0,69 1,95 88 V 2S2 (1,9±5,6)5
0,70 2,05 81 V2S2 (2,0±5,4) 5
0,71 2,05 85 V 252 (2,0±5,5)5
0,79 2,05 33 V 2S2 (2,0±3,5)5
0,80 2,05 29 V 2S2 (2,0±3,2)5
0,81 2,05 29 V 2S2 (2,0±3,2)5
0,89 1,45 7 V2 52 (1,4±1,5)5
0,90 1,35 5 V 2S2 (1,3±1,3)5
0,91 1,25 4 V 252 (1,2±1,2)5
0,99 0,15 0,05 V 252 (0,1±0,1)5
Коэффициент передачи, Kn Математическое ожидание, Tos Дисперсия, DU Интервал, Y0n
0,01 -12,35 410 V252 -(12,3 + 12,1)5
0,02 -12,25 410 V252 -(12,2 + 12,1)5
0,03 -12,25 410 V252 -(12,2±12,1) 5
0,04 -12,25 410 V252 -(12,2±12,1)5
0,05 -12,15 410 V252 -(12,1±12,1)5
0,06 -12,15 410 V252 -(12,1±12,1)5
0,07 -12,15 410 V252 -(12,1±12,1)5
0,08 -12,15 410 V252 -(12,1±12,1)5
0,09 -12,15 410 V252 -(12,1±12,1)5
Yo, n = А С1,1 + B0 C2,1 + SnC1, n + TnC2,
где A = —10га
2 A.l T2,1 . D 1 A 2 A,! A.l
2 ; B0 = 10га 2
T T
•4,0 4,0
L2LL
о 2 1, n rr 2 1, n 2, n
S =—га • —-L-—; T = га
T1, n + T2,n
T,n + T2,n '
Дисперсия амплитудной погрешности второй декады
D[r0nп] = АДСЦ] + БЩС2Л] -
-2А»Б0Л1 ЩСи \D\C\i | + ^[С,я] +
+ГИ2Ъ[С2п] -25иГиЛ/Л[С1,и \D\C2. ,, ].
В табл. 4 приведены результаты расчетов погрешностей при Рд = 0,997 и V = 0,2.
Экспериментальные исследования ИДН, проведенные на этапах проектирования и подготовки серийного выпуска таких установок, как В1-20, К2-41, Ф7090, подтвердили достоверность полученных результатов.
Таким образом, предложенный в работе метод расчета амплитудной погрешности как случайной величины применим для расчета многокаскадных ИДН, делительные обмотки которых могут содержать любое число секций с известными распределениями шунтирующих емкостей секций и статистическими характеристиками параметров секций. Аналитические выражения амплитудных погрешностей позволяют выявить степень влияния отдельных каскадов на результирующую случайную погрешность многокаскадного ИДН, нормировать его номинальную АЧХ и допускаемые отклонения от нее.
СПИСОК ЛИТЕРАТУРЫ
1. Нормирование и использование метрологических характеристик средств измерений: Нормативно-технические документы. — М.: Изд-во стандартов, 1985. — 152 с.
2. Байков В.М. Трансформаторные делители напряжения. — М.: Машиностроение, 1984. —52 с.
3. Ройтман М.С., Калиниченко Н.П. Индуктивные делители напряжения // Измерения, контроль, автоматизация: Научно-техн. сб. обзоров ЦНИИТЭИ приборостроения. — М., 1978. — Вып. 2(14). —С. 24—32.
4. Ким В.Л. Расчет амплитудной погрешности индуктивного делителя напряжения // Измерительная техника. — 2004. — № 3. -С. 28—31.
5. Иыерс Р.Р. Модели и характеристики многопроводного жгута // Труды Таллинского политехнического института. — 1977. — № 432. —С. 77—88.
6. Ким В.Л. Новый способ расширения частотного диапазона индуктивных делителей напряжения // Электричество. —2003. — № 12. —С. 46—49.
7. Ким В.Л. Индуктивные делители напряжения с симметрирующей обмоткой // Известия Томского политехнического университета. — 2003. — Т. 306. — № 3. — С. 93—95.
