Введение

При разработке проблемы безопасности атомных электростанций (АЭС) одним из направлений работы является обеспечение безопасности АЭС при возможном ударе по оболочке полуоторвавшимся или свободным телом [1].

Задача в окончательном виде заключается в том, чтобы создать алгоритмы определения параметров удара по результатам измерений колебаний оболочек.  Однако для окончательного решения обратной задачи необходимо решить прямую задачу о колебаниях оболочки после удара материальной точкой (телом).

В научной литературе имеются решения задач об упругом ударе точкой по цилиндрической и сферической оболочкам [2]. Нет ясности с решением подобных задач при упруго–пластическом ударе. Необходимость же решения обратных задач актуальна, но такие решения авторам неизвестны.

Весьма актуальна проверка достоверности решения подобных задач с помощью эксперимента, так как отклик оболочки на упругое воздействие является сложным сигналом.

В работе была поставлена задача о создании компактного измерительного комплекса с использованием АЦП (аналого-цифровых преобразователей) на базе персонального компьютера.

Подобраны соответствующие схемы и АЦП, проведена отладка схем, создан пакет программ для записи сигналов (колебаний) и их обработки.

Весь комплекс состоит из персонального компьютера с платой АЦП,  соединительных проводов и испытуемого образца [3].

1. Описание экспериментальной установки

Экспериментальная установка представляет собой программно-аппаратный комплекс, предназначенный для записи свободных колебаний оболочки после удара телом или материальной точкой (рис. 1).

Для фиксации ускорений используются пьезоэлектрические датчики KD‑91. Ниже приведены индивидуальные технические характеристики использованного набора датчиков:

·               Коэффициент передачи B_ua (125 Гц) = 0,52 мВ/мс^-2.

·               Емкость в сочетании с кабелем =0,58 нФ

·               Емкость кабеля C_k=0,15 нФ (длина 1,5 метра)

·               Резонансная частота (для массы 200 г) f_R =50,0 кГц

·               Коэффициент направленности Г_90max=10,0%

·               Сопротивление изоляции R_is > 1000 МОм

Рис. 1. Структурная схема экспериментальной установки

Формируемый электрический сигнал через специальный кабель, входящий в комплект датчиков, подается на вход платы АЦП. При создании установки был использован модуль аналогового ввода РСД-812 (производство компании «Демос»). Модуль предназначен для преобразования значений аналоговых сигналов напряжения, поступающих с внешнего разъема по шестнадцати каналам, в цифровой код. Конструктивно модуль представляет собой печатную плату для персонального компьютера типа IBM PC, выполненную для шины ISA.

Модуль содержит 16-тиканальный мультиплексор, устройство выборки-хранения и АЦП. Кроме этого, в модуле имеются два цифро-аналоговых преобразователя (ЦАП) которые не используются в данной работе. Выбор канала и запуск АЦП осуществляется программным способом. Имеется защита аналоговых входов от перегрузки.

Максимально допустимое напряжение: постоянное U_пост < ±25 В; импульсное U_имп < ±200 В при t_имп < 10 мкс.

Диапазон входных сигналов ±5В. Разрешение – 10 разрядов (старший разряд знаковый). Время преобразования – не более 50 мкс, таким образом,  максимальная частота дискретизации сигналов составляет 20 кГц. Ток потребления: +5В – 550 мА, ±12В – 45 мА.

Программный компонент предназначен для управления платой АЦП, программируемым контроллером таймера, записи и последующего отображения оцифрованных сигналов с датчиков.

Программа для записи результатов эксперимента осуществляет следующие функции:

·        Калибровка таймера

·        Управление контроллером АЦП

·        Замер времени

·        Сохранение полученной серии данных

·        Графическое воспроизведение серии данных

Под серией данных понимаются данные, снимаемые с датчиков экспериментальной установки в результате одного опыта. Серия состоит из некоторого фиксированного числа замеров, производимых через равные промежутки времени. В проводимой серии экспериментов принят размер серии в 2048 замеров, но он может быть легко изменен. Временной интервал между замерами обусловлен техническими ограничениями используемого АЦП и составляет величину порядка 50 мкс (время преобразования данных на АЦП согласно его техническим данным).

Структурная схема программы приведена на рис. 2.

Рис. 2. Структурная схема программы

Программа реализована на языке Pascal (версия Borland Pascal 7.0) с применением графической библиотеки Graph и пакета Object Professional 1.0 (© Turbo Power Software).

Модуль ADC предназначен для управления платой АЦП на нижнем уровне. Модуль реализован на языке Assembler 80х86 (компилятор Turbo Assembler 1.0 © Borland), программный интерфейс представлен в виде модуля unit на языке Pascal.

Модуль Timer реализует управление встроенным таймером компьютера. Для отсчета времени используется второй канал таймера. Таймер работает на частоте 1.193 МГц.

Главный модуль программы (FIX) выполняет две основные операции. Первая – запись непрерывного потока данных заданной длины с одновременным отображением их на экране и сохранение полученных данных в файле на диске. Вторая – графическое воспроизведение сохраненных данных.

Каждая серия данных идентифицируется целочисленным номером и сохраняется в файле текстового формата под именем HITnnnnn.DAT, где nnnnn – номер серии данных.

Первая строка файла является комментарием. Вторая строка является заголовком и содержит в себе количество замеров данных, сделанных в этом эксперименте. Далее идут строки, содержащие данные. Первый столбец содержит время в мкс, второй столбец – данные, считанные с АЦП, в условных единицах. Величине 512 соответствует измеренный потенциал +5В, величине –512 соответствует измеренный потенциал –5В. Такой формат файла позволяет легко экспортировать данные для дальнейшего анализа в другие программы, например, Microsoft Excel.

В программе также предусмотрен отладочный режим калибровки по времени. Он предназначен для точного определения интервала времени между двумя измерениями.  Определение интервала  времени производится при помощи внутреннего таймера компьютера при выполнении 10000 замеров. Абсолютная погрешность определения времени при этом составляет величину порядка 10^-8 c.

2. Работа установки

Пьезоэлектрический датчик, установленный на оболочке, фиксирует колебания оболочки и генерирует электрический сигнал, пропорциональный ускорению той точки оболочки, где он установлен. Сигнал поступает на вход платы АЦП, установленной в персональном компьютере.

Для определения начала колебательного процесса применена следующая программная схема. При запуске установки в течение 0.6 секунд (10000 измерений) происходит считывание показаний датчика и определение уровня помех (максимального и минимального значения, считанного с датчика за этот интервал времени). В качестве окончательного уровня помех берется модуль наибольшего из этих двух чисел.

После чего происходит непрерывное считывание состояния датчика. Как только будет считан сигнал, превышающий максимальный уровень помех в 1.5 раза, происходит запуск эксперимента. Счетчик времени устанавливается в ноль и выполняется серия из 2048 замеров. Интервал между замерами составляет 60.6 мкс. Из них 50 мкс уходит на выполнение преобразования сигнала в цифровую форму, а 10.6 мкс уходит на обработку данных процессором компьютера (использовался IBM-совместимый персональный компьютер с процессором Intel 386DX, тактовая частота 40 МГц). При использовании более быстродействующего компьютера время обработки можно было бы сократить, но это незначительно увеличит общее быстродействие схемы. Итоговая частота дискретизации составляет около 16.5 кГц при максимально возможной 20 кГц в идеальном случае.

Таким образом, алгоритм работы программы в режиме записи следующий:

1.   запрашивается номер первой серии данных,

2.   запрашивается количество серий данных,

3.   производится  определение  уровня   шумов,

4.   производится ожидание начала ударного процесса,

5.   производится серия из 2048 измерений,

6.   производится отображение серии на экране в заданном виде,

7.   производится вывод запроса на сохранение серии в файле на диске. Серия может быть сохранена, измерение может быть повторено или может быть произведен выход из программы.

3. Программа обработки экспериментальных данных

Для статистической обработки массивов данных, получаемых в результате проведения экспериментов на описанной выше экспериментальной установке, создана отдельная программа на основе Microsoft Excel 97 и Visual Basic for Application (версия 5.0). При обработке определяется средний коэффициент пересчета значений, измеренных датчиком-акселерометром, в величины ускорения. Передаточная характеристика платы АЦП при этом принимается линейной. В качестве эталонных данных для получения передаточного коэффициента используются данные, полученные программой расчета, описанной ниже.

Кроме этого, программа строит зависимости средних и максимальных (пиковых) показателей ускорения от параметров эксперимента – скорости соударения шарика и оболочки, расстояния от датчика до места удара (разность линейных и угловых координат), массы шарика.

4. Расчет зависимостей перемещений, скоростей и
 ускорений от времени в различных точках полого
цилиндра после удара по нему материальной точкой

Для проведения теоретического анализа вышеописанного эксперимента авторами решена задача и разработана программа расчета зависимостей перемещений, скоростей  и ускорений от времени материальной точки и точки оболочки во время и после удара. Расчет производится в два этапа. Вначале происходит расчет ударного процесса в точке удара. После этого происходит расчет параметров колебаний в заданной точке полого цилиндра.

Рассмотрим удар материальной точки массой M с начальной скоростью V₀ по нормали к образующей полого цилиндра (цилиндрической оболочки). Датчики ускорения установлены так, что фиксируют радиальные ускорения оболочки. Уравнение движения материальной точки при ударе имеет вид

                                                                                         (1)

Здесь  – нормальная сила взаимодействия материальной точки и оболочки в месте контакта,  – неизвестная функция. Координата s фиксирует движение точки по нормали к оболочке в процессе удара.

Начальные условия для движения точки при ударе: при  , .

Обозначим радиальное перемещение оболочки в точке, где производится удар, через w. Перемещение материальной точки s складывается из радиального перемещения точки оболочки w и величины «местного смятия» α оболочки и материальной точки:

                                                                                               (2)

Местные деформации предполагаются упругими. Для зависимости упру­гого смятия материалов в месте контакта шарика и оболочки от величины силы удара используется закон Герца упругого смятия [4]. Для малых скоростей соударения использование закона Герца для упругого контакта допустимо, что подтверждается физическими тенденциями в эксперименте.

Для упругих деформаций величина нормальной контактной силы принимается равной

                                           ,                                                 (3)  где коэффициент K зависит от параметров соударяющихся тел:

                                .                                     (4)

Величина q_k является функцией отношения A/B: . Для случая внешнего контакта цилиндрической оболочки радиуса R со сферой (материальной точкой) радиуса R₁ коэффициенты A и B примут вид:

                             , .                                   (5)

Величины δ₁ и δ₂ характеризуют упругие свойства материалов соударяющихся тел (материальной точки и оболочки соответственно):

                                           ,                                                (6)

где μ_i – коэффициент Пуассона, E_i – модуль упругости материала (i=1 для материальной точки, i=2 для оболочки) [4].

Рис. 3. Обозначения переменных

Рис. 4. Оси координат и составляющие перемещения точек оболочки

Величина w определяется из решения задачи о колебаниях цилиндричес­кой оболочки под действием неизвестной контактной силы .

Решение будем искать в виде рядов по собственным функциям, для чего сначала определим собственные частоты оболочки. Не учитывая волновой характер распространения деформаций, и пренебрегая инерционными силами в направлениях u и v (по оси x и по касательной к оболочке φ, см. рис. 4) и нелинейными членами, получим дифференциальные уравнения движения цилиндрической оболочки в виде:

                                        (7)

где , , , , , ρ, E, μ – плотность, модуль упругости и коэффициент Пуассона материала оболочки; h – толщина оболочки [5].

Из системы (7) можно получить:

                                        (8)

где .

Решив первое уравнение для w, затем из второго и третьего уравнений получим u и v.

Рассмотрим цилиндрическую оболочку с граничными условиями Навье: ,  при ,  (где l – длина цилиндрической оболочки). Решение ищем в рядах, где члены ряда имеют вид:

                                                  (9)

, ,  .

Определим собственные частоты колебаний оболочки при заданных граничных условиях, так как они далее используются в решении для радиальных перемещений w под действием силы P(t).

Подставив третье выражение (9) в первое уравнение (8), получим решение для собственных частот оболочки

,

откуда получаем значения ω_n,m:

                       .                            (10)

Минимальная частота получается для  и .

Далее ищем решение уравнения движения оболочки в виде

      .          (11)

Здесь через q учитывается неизвестная нагрузка (сила), которая возникает при ударе. Применяя к уравнению (11) интегральное преобразование Лапласа, получим при нулевых начальных условиях:

       .            (12)

Здесь p – параметр преобразования Лапласа, черта над функциями  и  означает, что в уравнения входят изображения функций:

, .

Функции ,  разложим в ряды по собственным функциям :

, .

Пусть функции  ортогональны, тогда

                               (13)

.

Умножив уравнение (12) на  и проинтегрировав по поверх­ности оболочки, а также учитывая, что  отвечают уравнению

,

получим

                               .                                   (14)

Выполнив обратное преобразование Лапласа, получим

                                     (15)

При допущении о том, что сила давления при ударе в точке с координа­тами  сосредоточена и зависит от времени, имеем

                       .                            (16)

Подставляя (16) в , получим

                                  (17)

Решение для уравнения (8) получено для условий Навье, т.е.

                                                 (18)

Учитывая (18), для D_n,m получим выражение

                                          .                                              (19)

Поэтому выражение (17) запишется в виде:

               (20)

По этому выражению рассчитывается прогиб в любой точке цилиндрической оболочки, если известен закон P(t). В частности, в точке удара  величина w примет вид:

                                   (21)

Расчет ударной силы осуществляется по шагам. На первом шаге расчета перемещение оболочки в месте удара положим равным 0, а величину местной деформации определим из скорости движения материальной точки V₀ (рассматривается нормальный удар). Тогда:

                        , , ,                            (22)

где Δt – временной интервал первого шага расчета, a – ускорение материальной точки, M – масса материальной точки.

На всех последующих шагах расчета перемещения материальной точки и оболочки определяются по следующим формулам (см. (21)):

                                             ,

                              ,

,     (23)

                            , .

Результатом данного этапа расчета являются зависимости s(t), w(t), P(t). Кроме этого определяются длительность ударного процесса τ₁ и максимальная величина нормальной контактной силы P_max. Эти данные используются на втором этапе расчета, который производит определение перемещений, скоростей и ускорений заданных точек цилиндрической оболочки.

Для радиального перемещения точки цилиндрической оболочки с координатами ξ, φ при граничных условиях Навье используется формула (20). Для упрощения расчета в случае упругого удара зависимость P(t) аппроксимирована синусоидой:

                                  .                                       (24)

При этом действие ударной силы ограничено временем удара τ₁ и прекращается при .

Подставив (24) в (20) и взяв интеграл по времени, получим формулы для расчета величин перемещений. Дифференцируя, получаем выражения для скоростей и ускорений заданной точки оболочки:

       ,

  (25)

                (26)

               (27)

В (25), (26), (27) использовано обозначение .

Из выражений для  видно, что ряды для  являются сходящимися, а ряд для  – условно сходящимся. Поэтому расчеты рядов  стабилизируются при гораздо меньшем количестве членов (по n и m), чем расчеты ускорений . При точности расчетов для перемещений и скоростей точек оболочки в 1%, точность расчетов по ускорению ~10% (при 1000–1500 членов ряда). В расчетах использовались 2550 (m=1…50, n=0…50) членов ряда для контактной задачи и 1020 (m=1…20, n=0…50) членов ряда для расчета перемещений, скоростей и ускорений точек оболочки.

Рассмотренное решение было доложено на нескольких научных конференциях [6–8].

В данную программу введена также возможность расчета упруго–пластического удара. При этом используются экспериментальные данные, а именно α_r – величина остаточного смятия после разделения соударяющихся тел. Закон изменения нормальной контактной силы при нагрузке имеет вид:

                                           ,                                                (28)

где D – коэффициент, определяемый экспериментально (в отличие от расчетного коэффициента K – см. (4)).

В момент, когда скорость материальной точки V обращается в 0, определяется величина максимального смятия α_m и максимальное значение нормальной контактной силы P_max. После этого наступает разгрузка, и закон изменения нормальной контактной силы приобретает вид:

                                .                                     (29)

Программа имеет оконный интерфейс, а также режим пакетного счета, при котором возможно обработать несколько наборов исходных данных.

5. Эксперимент и сравнение с расчетом в задаче об
ударе материальной точкой по цилиндрической
оболочке

Проведение и обработка эксперимента

В эксперименте использовалась в качестве цилиндрической оболочки стальная труба длины L=1520 мм, радиуса R=51 мм и толщиной h=2 мм. Схема установки датчиков и места замеров указаны на рис. 5.

Труба в экспериментах была свободно оперта. В эксперименте в качестве материальных точек использовались стальные шарики диаметром 2r₁=28,6 мм, 2r₂=22,3 мм, 2r₃=12,2 мм, 2r₄=9,6 мм. Их массы соответственно равны m₁=95,5 г, m₂=45,3 г, m₃=7,4 г, m₄=3,6 г. Высота падения шарика обозначалась H и изменялась от 5 см до 20 см, что соответствовало изменению относительной скорости соударения шарика и трубы от 0,99 м/с до 1,98 м/с. Удар шариком осуществлялся под углом 90º к образующей трубы (нормальный удар). Основные эксперименты проводились при фиксации сигналов датчиком Д2 (см. Рис. 5). Место удара шариком изменялось фиксированно, т. е. замер сигнала производился в точках (позициях) l₁=100 мм, l₂=300мм, l₃=500мм, l₄=600мм, l₅=650мм, l₆=700мм.  Замер сигнала у датчика (позиция 0, l₀=18мм) производился на некотором расстоянии от него. Это расстояние определялось из конструктивных возможностей.

Рис. 5. Установка датчиков, позиции точек замера

Из визуального анализа экспериментальных кривых видно, что используемая частота дискретизации недостаточна для того, чтобы оценить форму колебаний и построить спектральные характеристики сигнала по частоте. Возможно оценить только их амплитуду. Попытка частотного анализа (на основе быстрого преобразования Фурье) пока не дала отчетливых результатов. Следует повторить эксперимент, используя более совершенную плату АЦП с частотой дискретизации не менее 100 кГц на канал при поддержке не менее 2 каналов.

Анализ экспериментальных данных был построен на базе рассмотрения средних значений модуля амплитуды сигнала за выбранный отрезок времени. В качестве временного отрезка выбирались 1, 1.5, 2 мс. Наилучшие результаты получены при рассмотрении отрезка 1 мс, так как колебательный процесс в основном сосредоточен в этом диапазоне. В 1 мс укладывается 18 измеренных показаний датчика. Экспериментальные данные были перенесены в Microsoft Excel, построены таблицы абсолютных значений показаний датчика и проведено их усреднение на задаваемом временном отрезке. На графиках (Рис. 7) приведены полученные результаты, вместе с ними построены сглаживающие кривые. Сглаживание проведено полиномом второго порядка по методу наименьших квадратов (использование встроенных возможностей Excel).

Величины среднего значения показаний датчика на графиках показаны как |W|_i. Так как в эксперименте использовалось 4 шарика различной массы, в качестве индекса выступает номер шарика.

Полученные экспериментальные данные нуждались в градуировке, так как точный передаточный коэффициент пьезоэлектрического датчика неизвестен. Для этого на основании имеющейся математической модели рассчитаны колебания оболочки при тех же начальных и граничных условиях и данных, какие были использованы при проведении эксперимента. Расчет проведен на временном интервале 1 мс с дискретностью 1 мкс. Следствием этого является большое количество точек к усреднению.

Определение средних значений модуля ускорения проведено аналогич­ным образом. Градуировка осуществлена на всем множестве эксперимен­тальных данных путем линейной интерполяции зависимости |w| (|W|)  (то есть расчетного среднего ускорения от экспериментального среднего ускорения; если бы измеряемые в эксперименте значения были бы непосредственно ускорениями, то эта зависимость должна была бы иметь вид |w| = |W|). В результате получен единый для всего массива экспериментальных данных передаточный коэффициент К_w . С использованием передаточного коэффициента можно сравнивать экспериментальные и расчетные данные.

Экспериментальные данные сравнивались с расчетными, которые полу­чены в результате решения задачи об ударе материальной точкой по цилинд­рической оболочке. Оси координат выбраны были следующим образом – вдоль оси трубы x, перпендикулярно образующей трубы z и по касательной к образующей y (, где φ – угловая координата). Для того, чтобы можно было представить, каковы параметры максимальной силы при ударе точкой по оболочке, на рис. 6 представлены зависимость максимальной силы удара от массы материальной точки при H=10см (V=1,4 м/с).

Рис. 6. Зависимость максимальной силы удара от массы точки

,м/с2

Рис. 7. Зависимость ускорения точки оболочки от времени

На рис. 7 представлена зависимость ускорения точки оболочки от времени ( ). Здесь рассматривается составляющая ускорения вдоль оси z. Именно эта составляющая фиксируется датчиками ускорений (инерционны­ми датчиками). На теоретическую кривую нанесены экспериментальные точки – ускорения, замеренные в эксперименте в соответствующие моменты времени. Указаны значения по пяти одинаковым экспериментам. Разброс в эксперименте при одинаковых условиях значителен, поэтому общую картину наложения экспериментальных данных на расчетную кривую следует признать удовлетворительной.

На рис. 8–9 представлены расчетная и экспериментальная зависимости ускорения точки оболочки от времени. Экспериментальная зависимость представлена отдельными точками, соединенными отрезками прямых. Зависимости построены при одинаковых параметрах удара в сходных пози­циях. Следует отметить одинаковые тенденции сдвига максимума ускорений вправо по оси времени с удалением точки удара от датчика, снижение уровня максимума ускорений в одинаковых пропорциях. Если проследить за смещением максимума по времени, то можно установить величину средней скорости распространения возмущения по трубе, равную ~1000 м/с.

Проведя обработку экспериментальных данных, определили среднее значение коэффициента пересчета экспериментальных данных в напряжени­ях в числовые значения с размерностью ускорения. После обработки экспе­риментальных и расчетных данных и после определения среднего значения коэффициента пересчета для экспериментальных данных удалось построить ряд зависимостей и сравнить расчетные и экспериментальные данные. На рис. 10 представлены зависимости средних значений ускорений в точке за­мера на временном отрезке 1 миллисекунда от расстояния, которое занимает точка удара от датчика Д2 – точки замера. На рис. 11 представлены зависи­мости средних значений ускорений в точке замера в зависимости от скорости удара для шариков двух различных масс. На рис. 12 приведены значения среднего ускорения по окружности трубы. Удар производился в позиции l₂=300 при переменном угле φ. На рис. 13 представлены зависимости средних значений ускорения в точке замера от массы ударяющей материальной точки. Следует отметить весьма сложные механизмы замера и расчета средних ускорений. Поэтому представленные экспериментальные результаты можно считать удовлетворительно совпадающими с расчетными данными.

Рис. 8. Теоретическая зависимость

Рис. 9. Экспериментальная зависимость  (5 замеров на позицию)

Рис. 10. Зависимость среднего ускорения от расстояния до датчика

Рис. 10 (продолжение). Зависимость среднего ускорения от расстояния до датчика

Рис. 11. Зависимость значений среднего ускорения от скорости удара

Рис. 12. Распределение значений среднего ускорения по окружности трубы

Рис. 13. Зависимость значений среднего ускорения от массы шарика

Рис. 13 (продолжение). Зависимость значений среднего ускорения от массы шарика

Заключение

Аналитически исследована задача об упругом нормальном ударе точкой по цилиндрической оболочке и упругих колебаниях оболочки после удара. Местные деформации учитываются по модели Герца. Представлен алгоритм расчета зависимостей ускорений точек оболочки от времени при различных параметрах удара (массе и скорости материальной точки, месте удара и т.д.). Результаты численных расчетов позволили построить теоретические зависимости ускорений в функции времени и интегральные теоретические зависимости средних ускорений точек оболочки от параметров удара.

Создана экспериментальная установка, оснащенная программно-аппаратным комплексом, позволяющим записывать и обрабатывать характеристики свободных колебаний оболочки после удара по оболочке телом (материальной точкой). Полученные теоретические зависимости удовлетворительно подтверждаются экспериментами и являются основой для решения обратной задачи.

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

Литература

1. B. J. Olma. Source location and mass estimation in loose parts monitoring of LWR’s. Progress in Nuclear Energy, 1985, Vol. 15, pp.583–594. Printed in Great Britain.

2. И. Г. Кадомцев, М. Р. Фрейгейт (Ростов на Дону). Удар массивного тела по шарнирно опертому сферическому куполу. Проблемы прочности, 1991, №9, стр. 57–59.

3. В. В. Дубинин, С. Н. Банников, Е. М. Пашовкин. Разработка экспериментального метода определения параметров колебаний оболочки после удара материальной точкой. МГТУ им. Н.Э. Баумана. М. 1999 - 28 с. Деп. в ВИНИТИ, № 2251-В99 от 09.07.99 г., БО-121, № 9, 1999 г.

4. В. Гольдсмит. Удар. Теория и физические свойства соударяемых тел. Изд-во литературы по строительству. М. 1965, 448 с.

5. В. Д. Кубенко, П. С. Ковальчук, Н. П. Подчасов. Нелинейные колебания цилиндрических оболочек. Головное изд-во издат. объединения «Выща школа». К. 1989, 208 с.

6. В. В. Дубинин, С. Н. Банников. Разработка экспериментального метода определения параметров колебаний оболочки после удара материальной точкой. Научно-методическая конференция, посвященная 35-летию факультета «Фундаментальные науки» МГТУ им. Н. Э. Баумана. Тезисы докладов. 20 декабря 1999 г. – М.: Изд-во МГТУ им. Н. Э. Баумана, 1999. 140 с. (с.с. 66–68).

7. В. В. Дубинин, С. Н. Банников. Разработка экспериментального метода определения параметров колебаний оболочки после удара телом. Тезисы докладов научно-технической конференции. 21–23 ноября 2000 г. В 2 ч. Ч. 2. – М.: Изд-во МГТУ им. Н. Э. Баумана, 2000. – 160 с. (с. 38).

8. С. Н. Банников, В. В. Дубинин, Ю. Н. Жигулевцев. Развитие технологий обучения: обработка физического эксперимента и дистанционное образование. Восьмой всероссийский съезд  по теоретической и прикладной механике, Пермь, 23–29 августа 2001 г. Аннотации докладов. Екатеринбург. УрО РАН, 2001. (с. 72).