НАВИГАЦИОННЫЕ И ГИРОСКОПИЧЕСКИЕ
СИСТЕМЫ
УДК 621.31(075.8)
С. А. Васюков
СИЛЫ И МОМЕНТЫ В ЦИЛИНДРИЧЕСКОМ ЭЛЕКТРОСТАТИЧЕСКОМ ПОДВЕСЕ
Рассмотрены вопросы силообразования в цилиндрическом электростатическом подвесе. Проанализированы статические силовые и моментные характеристики для подвесов, работающих на постоянном токе, как при малых, так и больших смещениях оси ротора. Рассмотрено влияние неидентичности зазоров левой и правой электростатических опор на силовые характеристики подвеса.
В акселерометрах, относящихся к широко применяемому классу поплавковых маятниковых приборов и работающих в компенсационном режиме, подвижная система (рис. 1) содержит два цилиндрических поплавка, которые служат для обеспечения гидростатической разгрузки. Поверхности поплавков могут быть использованы в качестве роторов электростатических подвесов (ЭСП), выполняющих функции дополнительных прецизионных элементов центрирования подвижной системы.
В отличие от вакуумных подвесов в поплавковых демпфирование реализуется посредством вязкой жидкости, тем самым легко достигается устойчивое взвешивание. Подвижная система акселерометра, охваченная обратной связью по измерительной оси, практически не совершает разворотов вокруг нее. Конечно, и в этом случае имеет место воздействие ЭСП на точностные характеристики прибора через возмущающий момент на измерительной оси. Однако возможность высококачественного изготовления цилиндрического ротора подвеса существенно снижает влияние момента и делает его вклад в уравнение движения вокруг измерительной оси пренебрежимо малым.
Рис. 1. Подвижная система ЭСП
Взвешенный элемент, окруженный системой электродов, схематично изображен на рис. 1 (здесь ДМ — магнитоэлектрический датчик момента).
Максимальные силы и моменты, развиваемые электростатическим подвесом, зависят от габаритных размеров поплавков и электродной системы. Эти размеры во многом определяются условиями, заданными при проектировании прибора. Действительно, заданная при проектировании маятниковость прибора Мм является исходной для определения его габаритных размеров 13 и 14, с которыми непосредственно связан максимальный радиус поплавкового узла атах < 14/2.
Габаритная длина прибора Ь, заданная при проектировании, складывается из габаритной длины подвижной системы, размеров сильфо-на, корпуса и др. Опыт проектирования показывает, что длина поплавка 1п < (0,7Ь - 13)/2.
При определенных габаритных размерах подвеса ключевыми параметрами, влияющими на значения максимальных сил и моментов, остаются зазоры электрод-ротор Н и 5, а также предельная напряженность электрического поля в жидкости Етах, с которой непосредственно связано значение максимального опорного напряжения на электродах подвеса.
В результате работы с поплавковыми приборами выявили, что для достижения оптимальных параметров подвеса зазоры между ротором и электродами следует выбирать из следующих неравенств:
30 мкм ^ Н ^ 70, мкм;
30 мкм ^ 5 ^ 70 мкм.
Выбор зазора менее 30 мкм нецелесообразен, поскольку чрезмерное увеличениие линейного и углового демпфирования ведет к затягиванию процесса всплытия и переходных процессов при действии перегрузок, а также из-за негативного влияния посторонних примесей в жидкости, остающихся после фильтрации, воздушных пузырьков и влаги.
Зазоры более 70 мкм служат причиной заниженных величин демпфирования, что создает трудности в обеспечении устойчивости при больших жесткостях подвеса. Важно также отметить, что большие зазоры приводят к необходимости повышения напряжения на электродах для поддержания необходимой перегрузочной способности. Это существенно усложняет схемотехнику подвеса.
Для уменьшения вероятности электрического пробоя в жидкости при начальном всплытии чувствительного элемента (ЧЭ) рекомендуется устанавливать в приборах специальные страховочные упоры, с
которых происходит всплытие ротора. Обычно в качестве упоров применяют камневые опоры (см. рис. 1). Общий ход оси подвижной системы в упорах в радиальном направлении ограничивается величиной ^уп = 2Ду, а в осевом направлении — ^уп = 2Д^. Типовое значение общего хода цапфы в камневой опоре лежит в диапазоне 5... 10мкм, т. е. составляет от 7 до 30% рекомендуемых зазоров электрод-ротор.
Необходимо учитывать, что из-за конструктивных особенностей приборов с цилиндрическим ЭСП (ЦЭСП) не удается обеспечить полного равенства зазоров левого и правого подвесов. Очевидно, что неидентичность зазоров отражается на силовых и моментных характеристиках. Электростатический подвес поплавкового маятникового акселерометра должен устойчиво функционировать при перегрузках и обеспечивать необходимый уровень жесткости.
Задачей исследования является: получение силовых и моментных характеристик подвеса и анализ их отдельных составляющих; определение глубины перекрестных связей, вызванных неидентичностью левого и правого подвесов; определение диапазонов действия восстанавливающей силы и момента при смещениях ротора в различных направлениях.
На рис. 2 показана расчетная схема подвеса. Взвешиваемое тело цилиндрической формы окружено системой электродов, представляющих собой тонкие металлизированные площадки, расположенные на внутренней поверхности корпуса. В принятой модели межэлектродные промежутки, располагающиеся на одной с электродами поверхности, также металлизированы и объединены в один электрод, который назовем внешним проводящим экраном. Экран изолирован от остальных электродов, и его потенциал принимаем равным нулю.
Присвоим ротору индекс г = 0, электродам г = 1,..., 10, экрану г = 11. Потенциал и заряд г-го проводника обозначим как и qi соответственно. Тогда для системы проводников, составляющих электростатический подвес, можно записать
n
qj = 2^ °з ^ ^ =0,...,щ (1)
г=0
где Оу — коэффициенты электростатической индукции, которые связывают заряды qi тел с их потенциалами Эти коэффициенты являются сложными функциями, зависящими как от геометрических размеров подвеса, так и от линейных и угловых смещений оси ротора.
Обозначим Д^ произвольное малое линейное смещение оси ротора. Тогда выражение для проекции силы на направление £ можно записать как [1]
1 „ „
2
f = 2 ЕЕ ^v• (2)
i=0 j=0
Рис. 2. Расчетная схема сил (а) и моментов (б) конкретного подвеса
Предположим, что в качестве малого смещения было принято угловое перемещение А9, тогда для момента сил электрического поля получили бы следующее выражение:
1 n n M = & V Cij. i=0 j=0
(3)
Выражения (2) и (3) позволяют установить зависимость центрирующей силы и момента подвеса от потенциалов на его электродах, которые меняются при регулировании, а также от параметров перемещений ротора, входящих в эти формулы посредством коэффициентов Сц.
Анализируя выражения (2) и (3), можно отметить, что их структура представляет собой квадратичную форму, которую удобно записать в виде произведения вектора-столбца (т = [(0, (,..., (п], квадратной матрицы градиентов коэффициентов электростатической индукции уС = у {Сц} (г, ] = 0 ... п) и вектора-строки (р:
F = 1 C<.
(4)
Выражение (4) удобно использовать для вычисления сил (моменты вычисляются аналогично) конкретного подвеса на основе сформированных полных матриц коэффициентов электростатической индукции.
В работе [2] были получены базовые коэффициенты электростатической индукции как функции нормированных линейных X, y, z и угловых а, в смещений ротора из центра подвеса. Если воспользоваться преобразованиями координат, то можно на основе базовых коэффициентов построить полные матрицы градиентов VC и определить проекции сил в ЦЭСП на координатные оси. Здесь следует отметить, что в силу симметрии подвеса достаточно найти выражения сил и моментов, действующих в какой-либо одной плоскости, например XOZ. Остальные проекции легко получить перестановкой индексов.
Введем коэффициенты неидентичности левого (ci) и правого (с2) подвесов. Если номинальный зазор при центральном положении ротора равен h, то с введенными коэффициентами зазоры левого и правого подвесов выражаются как
h-л = h/ci; h-п = h/c2. (5)
Тогда суммарную проекцию силы левого и правого подвесов на ось X можно записать следующим образом:
Fx =
2h
c1(<-<4) + c2(<7-<s) <0 + 2CX3+y c{(ipi + <>)(<3-<4) +
+ c2(<5 + <б)(<7 - <s)
x f 1
+ Cy
11
c1(<2 - + c2(<7 - <S)
+
+ hi 2(cl + c3^<0 + 2(o <0 [c1(<1 + <2) + c3(<5 + <6)
+
+ 2Cy/() <0 c1(<3 + <4) + c2(<7 + <s)
+
+ 2C
13
c1(<i + <2)(<3 + <4) + c3(<5 + <б)(<7 + <s)
+
+ CX
c3(<2 + <2) + c2(<5 + <2)
+ Cf
+ CX2
33
c-<1<2 + c3<5<6
+ 2C?2
c3(<2 + <4) + c2(<7 + <i)
cf<3<4 + c3<7 <8 > +
+
в
+ 2h ^ 2Св?0Ж<0 c1(<1 + <2) - c2(<5 + <6)
+ 2C1c0/<0 c3(<3 + <4) - c3(<7 + <s)
+ Cf
yll
+
c3(<2 + <2)-
- +
+ cay
11
+ - + d)
+
+ 2C
12
_ 2C
3 3
+ 2CäS
12
33
+ ^2)(^2 + (4) _ ^((5 + + (s)
+
1
+h(C^^+C
-уаж+/3;у —
13
а
c1(^1_^2 )(^2_^4 )+c2(^5_^6)(^7_^s) . (6)
Найдем теперь проекцию силы на ось Z.
Z 2h
2C90^0(^9 _ ^10) + CS^9(^2 _ ^10)
+
z +hl
C00 ^ + 2C9Z0 + ^10) + C99 (^2 + ^20) . (7)
Чтобы вычислить проекцию момента на ось У, необходимо предварительно сформировать матрицу градиентов коэффициентов электростатической индукции по углу в:
1
My = 2h 12ca [c? ((4 _ (2) + c2 (^r _ ы] ^0+ + 2ci2+e [c! (^ + ^2) (^4 _ ^2) +
+c2 (^5 + Ы (^r _ Ы] + cf1 [c1 (^ _ + c2 (^r _ p8)] \ +
+ h|2Ce?02(0 H (^1 + ^2) + c2 ((5 + Ы] +
+ 2C002(0 [c1 ((2 + (4) + c2 (^r + (s)] +
+ 2C02+/? [cf (^1 + + (4) + c2 ((5 + Ы (^r + (S)] +
2 ( 2 ) ( )] [ ( ) (
+С^1 [с? (( + (;) + с3 (( + (2)] +с1с\ [с1 (( + + с2 (( + (!)] + + С^^2 [с?(1(2 + с2(5(б] + 2С022 [с?(3(4 + с2(7(8^ + + {2С£)>0 [с1 ((1 + (2) - -с2 ((5 + (б)] -- 2С0о'со [с3 ((3 + (4) - с2 ((7 + (8)] + СЙ8 [с? ((2 + (2) -
-с2 ((5 + (2)] - СОТ [с3 ((3 + (2) - с2 ((2 + (8)] +
+ 2С{?2С [с1(1(2 - ¿^(б] + 2Са2 [с1(7(8 - ^(3(4] -
- 2С17+^ [с3 ((1 + (2) ((3 + (4) - с2 ((5 + (б) ((7 + (8)] } +
+ ; (с^ Й-^+^ур) [с3 ((1 -(2)((3-(4) +с2 ((б-(5)((7-(83 . (8)
Проекция момента на ось Z равна нулю, так как в рассматриваемой модели поверхность цилиндра предполагалась идеальной.
Соотношения (6)-(8) являются функциями потенциалов на электродах, смещений ротора из центра подвеса, а также функциями потенциала ротора.
Проанализируем силовые характеристики цилиндрического подвеса, работающего на постоянном токе. Управление в подвесах осуществляется по следующим законам:
V = - Д^г)Дг; Vj = (j + ; (9)
здесь i, j — номера пар диаметрально противоположных электродов.
Для цилиндрического подвеса номера пар (i, j) принимают значения (1, 2), (3, 4), (5, 6), (7, 8), (9,10). Коэффициенты Дг, Д.,- равны плюс или минус единица. Потенциал ^г0 является положительным и представляет собой начальную уставку или опорное напряжение на электродах. Приращение потенциала Д^г в статическом режиме пропорционально напряжению на выходе датчика перемещений идг и не превышает по абсолютной величине потенциала уставки:
{Vi0 при Мдг > Vio/куг;
кугМдг при |Мдг| < ^0/^ (10)
-^г0 при Мдг < -^го/куг.
здесь куг — коэффициент усиления следящей системы, подключенной к паре электродов (i, j).
Коэффициенты Дг и Д. подбираются таким образом, чтобы минимизировать потенциал ротора. В работе [3] рассмотрены несколько вариантов управления с точки зрения минимизации потенциала ротора. Рассмотрим вариант, для которого
Дг = 1, i = 1, 2, 7, 8, 9; Д. = -1, j = 3, 4, 5, 6,10. (11)
Проанализируем сначала случай малых смещений (подвес работает в линейной зоне), при котором потенциалы на электродах не выходят на режим насыщения.
Примем ^г0 = ^0э (i = 1,..., 10), т.е. по абсолютной величине опорное напряжение на всех электродах подвеса одинаковое. Подставляя уравнения (9), (10), (11) в выражение (8), получаем проекцию силы ЦЭСП на ось X:
5 2 /х
Fx=E F«; Fxi=A?(c^3 + c^7); A = — (С? - 2С1а;з+у)
г=1
FX2=A2x; A2 = 2^0э (c3 + c3)(Cfi + Cj? + Cf2 + Cg - 4Cl3+y 2) /h; Fxз=Аз/5; A3 = ^(c? - c3)(Cff - Cff + Cf2f - C2 + 4С0Гв1 )/h;
FX 4 — —
2x h
Cg2 _ Cg2) (c2A(2 + c2A(5) + + (Cg _ Cg) (c1 A(2 + c2A(?)
+
+h
(Cff _ Cgf) (c1 Atf _ c2A(2) _
_ (Cg _ Cg) (c2A(2 _ c2A(?) _ h (CgУ _ Cgа) (c2A(1 A(2 + c2A(5A(r);
(12)
F
4 h
_2(0 " h
X5—
Cf0(c2A(2_c2A(r)+2X(03(c2 _ c2) Cg _ Cg, +
+ /5(0з (c2 + c2) (Cg + Cg) + ^Cg2-
(0X
Проанализируем отдельные составляющие силовой характеристики.
Составляющая пропорциональна приращению потенциалов на электродах левого (Д(3) и правого (Д(7) подвесов, причем оси электродов, формирующих эти приращения, совпадают с осью X. Найдем выражения Д(3 и Д(7. Для этого воспользуемся соотношением (10), а также обобщенными выражениями для перемещений датчиков цилиндрического подвеса, полученными в работе [3]:
(13)
Ид1 = кд1ур - кд2«; Мд3 = кд1 р + кд2/5; «д5 = кд1ур + кд2<Р; Мд7 = кд1 р - А^/З; Мд9 = -кд3р.
Подставляя в уравнение (10) выражения (13), получаем следующие значения приращений потенциала:
Д(3 = с1 (^1 р + ^2/5) ; Д(7 = с2 (^1 р - ^2/3) ; (14)
здесь к1, к2 — коэффициенты усиления в канале регулирования, которые содержат составляющие от датчика (г = 1, 2) и от дополнительного усилителя ку1 радиального канала к1 = кд1ку1, к2 = кд2ку1. Используя выражения (14), запишем соотношение для составляющей 1 в статическом режиме:
Fx1 — A1[(c1 + c2) k1 x + (c1 _ c2) k2/3]
(15)
Видно, что имеет не только полезную составляющую, пропорциональную р, но и перекрестную связь по /, глубину которой можно
оценить коэффициентом, %
I с3 _ I
щ = • 100%.
При отклонении зазоров подвесов в разные стороны от номинального зазора на 5 % коэффициент перекрестной связи принимает значение 14,9%, если же подвесы сделаны некачественно и отклонение составляет 10%, то п1 = 26,6 %.
При с1 = с2 перекрестная связь исчезает и 1 можно представить
как
-к
= сх1 X; сх 1 = 2^*1 = -(С* - 2СХз+у) ^э,
где сх 1 < 0. Сила и перемещение здесь имеют разные знаки, поэтому можно говорить о том, что стремится ликвидировать смещение. В связи с этим составляющую 1 будем называть стабилизирующей.
В противоположность 1 составляющая всегда имеет положительную крутизну по X, т.е. является дестабилизирующей составляющей. Естественно, одним из условий устойчивости подвеса является преобладание стабилизирующей составляющей над дестабилизирующей, т.е. должно выполняться соотношение 2/^Х1| < 1. Если учесть, что *1 = Д^/X = ^Сэ/Хтах, где Хтах — линейная зона подвеса, на границе которой Д^ = ^Сэ, то при с1 = с2 получаем (в процентах)
П2 =
F
X2
F
X1
100 =
Xm
12 C11
+ Cj1 + Cf22 + Cj2 - 4C
12
12
1з
r<y onf+y
C11 - 2C13
100. (16)
Из выражения (16) видно, что дестабилизирующая составляющая увеличивается с ростом линейной зоны Хтах; с этой точки зрения наиболее выгодным является подвес с нулевой линейной зоной (релейное управление).
Коэффициенты электростатической индукции, входящие в выражение (16), сложным образом зависят от таких параметров подвеса, как зазор, линейные и угловые размеры электродов, их взаимного расположения, и проследить аналитическую зависимость п2 от этих параметров не представляется возможным. Учитывая это, оценим влияние зазора Н (как одного из основных параметров, определяющих максимальную силу подвеса) на п2 для конкретных типоразмеров опытного образца прибора с цилиндрическим подвесом (см. рис. 2):
Ь = 55 мм; /1 = 1 мм; /с = 11 мм; #с = 50°; а = 7,5 мм.
Воспользовавшись соотношениями для коэффициентов электростатической индукции, полученными в работе [2], определим, что при зазоре Н = 50 мкм и Хтах = 0,1, коэффициент п2 = 11,36 %. С ростом зазора коэффициент п2 возрастает (правда, незначительно): П2 = 11,43 % (Н = 200 мкм); п2 = 11,64 % (Н = 500 мкм); п2 = 13,41 % (Н = 3000 мкм).
Слагаемое 3 представляет собой перекрестную связь по углу в, которая возникает при неидентичности зазоров подвесов. Оценим ее вклад (в процентах) коэффициентом
Пз =
F
хз
• 100 =
FX1 ,
X f^ay . f^ßx ^-»öy . f^ay+/3X _C11 — C11 + C12 — C12 + 4C13
2(C?i - 2CX+y)[(c? + c2)x/Xmax + (c3 - c2)/?//5max]
(c3-c2)e
•100.
При С1 = 1,05, С2 = 0,95, Хтах = 0,1, /тах = 0,2, X = Хтах, в = втах, Н = 50 мкм коэффициент п3 = 0,65 %, т.е. довольно незначителен.
Составляющие ^х4 и 5 характеризуют собой нелинейности, нелинейные перекрестные связи и перекрестные связи через потенциал ротора. При малых смещениях ротора из центра подвеса их значения не превышают десятых долей процента.
Учитывая численные оценки, можно утверждать, что при работе подвеса в линейной зоне при смещениях оси ротора, составляющих не более 10% зазора, обосновано использование линеаризованных силовых характеристик, которые в статическом режиме имеют вид
^х = XXX + Хвх XX = (с? + с3) + А2;
(17)
Хв = (с3 - с2) А? + Аз,
где Хх и Хв — линейная и угловая жесткости.
Перейдем теперь к вычислению моментной характеристики Му. Так же как и силовая характеристика , она содержит пять характерных составляющих:
Му = ¿ Муг; Му 1 = В? (с?Д^з - с2Д^7);
г=1
2 г / _ - \
В = — (С- - 2СГз+в) ^сэ; Му2 = В2/Х; Му3 = -^Гй;
В = 2^0э (с3 + с2) (с-2 + С- + С«22 + С- - 4^) Г/Н;
^ [(СГ - С?) (с!д^2 + с3Д^2) +
MY 4 =
+
xL
+ XL
+ (c?í - C?2a) (c?A(2 + c2A(?)
Cff-Cf2x) (c?A(2 - c2A(2) - (Cff-Cfl) (c?A(2+c2A(?)
- 4L (cf/a - (cfA^A^ - c3A^A(7)
My5 = [-CÍ5, (c?A^3 - c2A(7) + 2^оэ (c3 - c2) (cg - Cg?) +
+ X^L (c? + c2) (cegf + cg)
(18)
Стабилизирующую составляющую в статическом режиме с учетом выражений (14) запишем в виде
Myi = Bi[(c? + c?) М + (c? - c?) kix],
она содержит перекрестную связь по X, исчезающую при c1 = c2. Соотношения между отдельными составляющими в (18) представляют собой величины того же порядка, что и в (12), поэтому правомерно применить линеаризованную статическую моментную характеристику
My = се в + сх X; Се = (c1 + c?) k2Bi + B2;
(19)
Сх = (c1 - c2) kiBi + A?L.
Рассмотрим случай больших смещений ротора. Составляющие перемещения по осям X и Y (см. рис. 2), выраженные через абсолютное значение смещения d и направление смещения можно записать как
x = dcos y = dsin, d = -\/x2 + y2.
Угол рассогласования между вектором силы и направлением смещения будет равен
Fx x + Fy y Fsd
ф = arccos F X^ YУ; F = л/F2 + F2.
Очевидно, что направление полезной составляющей силы противоположно направлению смещения. Действительно, для возвращения ротора в центр подвеса необходимо, чтобы сила и смещение имели разные знаки. Назовем полезную составляющую восстанавливающей силой и запишем выражение для нее: F = Fs cos ф. Сила будет восстанавливающей в том диапазоне смещений, где F < 0.
Проекции сил FX и FY рассчитывали по обобщенному выражению (6). Силовые характеристики цилиндрического подвеса в функ-
Рис. 3. Силовые характеристики цилиндрического подвеса (левая и правая опоры идентичны):
1 - Фо = 0о; 2 - Фо = 15°; 3 - Фо = 30о; 4 - Фо = 45°; а = в = 0°; г = 0
ции нормированного смещения ¿/к приведены на рис.3. Коэффициенты усиления по осям X и У выбирались так, чтобы линейная зона подвеса была в пределах 10 % зазора к. Анализ характеристик показывает равножесткость при смещениях в любом направлении. Назовем диапазоном действия восстанавливающей силы предельное смещение ротора, при котором сила и перемещение имеют противоположные знаки. Для рассматриваемого случая диапазон действия восстанавливающей силы (¿/к)тах = 0,43.
Рассмотрим, каким образом линейная зона подвеса влияет на вид силовой характеристики. На рис. 4, а приведены силовые характеристики при смещении = 45° с линейными зонами от 4 до 10% зазора. Максимальная сила при уменьшении линейной зоны растет, однако диапазон восстанавливающей силы не меняется.
И наконец, на рис. 4, б показано изменение силовых характеристик для трех значений зазора при сохранении одного и того же значения опорного напряжения на электродах. Также наблюдается изменение максимальной силы при сохранении диапазона (все характеристики пересекаются в одной точке при ^ = 0).
Полученные результаты позволяют надлежащим образом выбрать размер упоров ¿уп.
В результате анализа силовых и моментных характеристик выявили, что при выполнении надлежащих мер по минимизации потенциала ротора и при малых смещениях ротора из центра подвеса возможна линеаризация силовых характеристик, что значительно упрощает исследование подвесов. На этапе начального всплытия ротора с упоров и для расчета предельного диапазона действия силы при смещении в произвольном направлении необходимо использовать обобщенные
Рис. 4. Силовые характеристики цилиндрического подвеса при смещении ф0 = 45° (а) и при трех значениях зазора h (б)
силовые характеристики (6)-(8) и учитывать конкретный закон управления потенциалами.
СПИСОК ЛИТЕРАТУРЫ
1. Т а м м И. Е. Основы теории электричества. - М.: Наука, 1966.
2. В а с ю к о в С. А., Д р о б ы ш е в Г. Ф. Математическая модель цилиндрического электростатического подвеса как системы заряженных проводников // Вестник МГТУ им. Н.Э. Баумана. Сер. "Естественные науки". - 2007. - № 1. -С.72-87 .
3. Васюков С. А., Дробышев Г. Ф. Алгоритмы управления потенциалами на электродах электростатического подвеса // Вестник МГТУ им. Н.Э. Баумана. Сер. "Приборостроение". - 2007. - № 2. - С. 69-81.
Статья поступила в редакцию 11.05.2007
Сергей Александрович Васюков родился в 1958г., окончил в 1981г. МВТУ им. Н.Э. Баумана. Канд. техн. наук, доцент кафедры "Электротехника и промышленная электроника" МГТУ им. Н.Э. Баумана. Автор 32 научных работ в области электростатических подвесов и датчиков автосигнализаций.
S.A. Vasyukov (b. 1958) graduated from the Bauman Moscow Higher Technical School in 1981. Ph. D. (Eng.), assoc. professor of "Electrical Engineering and Industrial Electronics" department of the Bauman Moscow State Technical University. Author of 32 publications in the field of electrostatic suspensions and sensors of automatic signaling system.