Дискретная модель процесса распространения импульса смещения в упругом стержне постоянного сечения при торцевом ударе

А. В. Затылкин, Г. В. Таньков, Д. В. Ольхов

ДИСКРЕТНАЯ МОДЕЛЬ ПРОЦЕССА РАСПРОСТРАНЕНИЯ ИМПУЛЬСА СМЕЩЕНИЯ В УПРУГОМ СТЕРЖНЕ ПОСТОЯННОГО СЕЧЕНИЯ ПРИ ТОРЦЕВОМ УДАРЕ

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

Ключевые слова: дискретная модель, импульс смещения, упругий стержень, торцевой удар.

Введение

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

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

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

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

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

Разработка дискретной модели распространения импульса смещения

в упругом стержне

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

мент Ь = о действует кратковременный удар, создающий импульс сжатия. Продольные колебания, возникающие при этом в стержне, описываются волновым уравнением

Е Э2и д2ы

Е-----= р---,

Эх2 Эх2

где и(х, Ь) - смещение текущего сечения стержня вдоль его оси х; Е - модуль Юнга; р - плотность материала.

Полагая, что с = Е - скорость распространения продольных волн в стержне, по-V Р

лучим:

2 Э2и Э2и , .

с2----=-----.

Эх2 ЭЬ2

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

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

X

-<рЛгОЛюЛ<>ЛгНк>Лн

I

Рис. 1. Геометрическая дискретная модель стержня

Заменив производную в левой части (2) ее разностным аналогом, получим

с2 д2и

—

к д

где их - смещение текущего узла вдоль оси х. Учитывая, что вторая производная от пере-

д2и

мещения по времени есть ускорение а узла, запишем (3) в виде а =-------- и, заменив про-

дЬ 2

изводные по времени разностным аналогом, получим:

т2а = их

где т - шаг дискретизации по времени.

Преобразуем (4) к виду явного разностного уравнения

т2а + 2их

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

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

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

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

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

Задача 1. Рассматривается процесс распространения по стержню импульса смещений от действия мгновенного ударного импульса. Правый крайний узел модели (см. рис. 1) закреплен -их=1 = о; в левый крайний узел (свободный конец стержня) в момент Ь = о

(рис. 2,а) задается единичное смещение и их=0 = 1, соответствующее действию мгновен-

( л ^

ного ударного импульса. Дойдя до закрепленного конца I Ь = — I, импульс смещения от-

I 4)

ражается с изменением знака смещения (рис. 2,6) и возвращается к левому свободному концу , при отражении от которого знак импульса не меняется (см. рис. 2,6). Дой-

дя до закрепленного правого конца , импульс снова отразится с изменением зна-

ка смещений (см. рис. 2,6) и вернется к свободному концу (Ь = 2л) с тем знаком, с каким началось его движение в момент Ь = о (см. рис. 2,а), т.е. картина полностью повторится.

Проведение вычислительных экспериментов

1=0

> 1

а)

1

1

1

б)

Рис. 2. Распространение одиночного импульса вдоль стержня

Полученная в численном решении картина процесса полностью соответствует нормальным колебаниям сплошного стержня, полученным в механике аналитическим путем [4].

Задача 2. В левый свободный узел модели задаются смещения, соответствующие действию прямоугольного ударного импульса (рис. з,а) и полусинусоидального импульса (рис. 4,а) в момент Ь = о. Длительность импульса в общих случаях выбирается равной 15 мкс. На графиках рис. 3,6 и 4,6 показаны картины распространения импульсов конечной длины, особенности взаимодействия падающего и отраженного импульсов на торцах

л л 3л

стержня (в моменты времени Ь = —, Ь = —, Ь =—, Ь = 2л) и моменты смены «знака» им-

424

пульсов. Это численное решение также подтверждается аналитической картиной процесса механики сплошных сред [5, 6].

1=0

а)

1=п

4

1=3п 1 4

1=п

1=3%

1=п

_3п

1=4

1=п

1=4

1=п

1=2п

1=2п

1=2п

1=2п

1=2п

б)

Рис. 3. Отражение границ стержня прямоугольного импульса

4

1

п

2

б)

Рис. 4. Отражение от границ стержня полусинусоидального импульса

Рис. 5. Отражение от границ стержня синусоидального импульса

83

Задача 3. Рассматривается действие на свободный край модели (рис. 5,а) синусоидального импульса длительностью 30 мкс (фактически короткая бегущая гармоническая волна). Из механики известно [5, 7], что в этом случае от закрепленного конца стержня волна смещений отражается с поворотом фазы на п; в случае же свободного конца стержня она отражается без изменения фазы. Эта аналитическая картина также соответствует результатам численного решения, представленного на графиках рис. 5,6.

Задача 4. Рассматривается задача получения в модели стержня стоячей волны. Из механики известно [5, 8], что если при распространении и отражении волны не происходит их затухания, то обе волны - падающая и отраженная - будут иметь одинаковую амплитуду, но фазы в какой-либо точке будут различны.

Эта картина получается и при численном решении на дискретной модели стержня. На рис. 6,а показана графическая картина распределения амплитуд смещений вдоль дискретной модели стержня, когда правый граничный узел закреплен, а в свободный левый граничный узел задается гармоническое колебание с частотой третьей гармоники и амплитудой 1 мкм (рис. 6,6).

и

т

6) г)

Рис. 6. Сложение падающей и отражающей волн

І

І

Из рис. 6,а видно, что в модели установилась стоячая волна смещений, когда на обоих концах модели получились узлы смещений, а по ее длине укладывается целое число полуволн. Но в этом случае амплитуда ближайшей к левому концу стержня пучности должна быть значительно большей по отношению к заданной амплитуде смещений левого конца стержня [5]. Это подтверждается в расчетах: седьмой узел модели (пучность) имеет амплитуду колебаний 8 мкм (в масштабе графика рис. 6,в). Узел, где синусоида проходит через нуль (пятнадцатый узел модели, см. рис. 6,а) имеет в численном решении амплитуду смещения, близкую к нулю (рис. 6,г).

Выводы

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

Список литературы

1. Слепян, Л. И. Нестационарные упругие волны / Л. И. Слепян. - Л. : Судостроение, 1972. -376 с.

2. Тартаковский, А. М. Математическое моделирование в конструировании РЭС : моногр. /

А. М. Тартаковский. - Пенза : Изд-во Пенз. гос. техн. ун-та, 1995. - 112 с.

3. Абрайтис, Л. Б. Автоматизация проектирования ЭВМ / Л. Б. Абрайтис, Р. И. Шейнаускас,

В. А. Жилевичюс ; под ред. Л. Б. Абрайтиса. - М. : Сов. радио, 1978. - 272 с.

4. Рындин, Д. А. Система генерации тестового сигнала для исследования динамических характеристик элементов конструкций РЭС / Д. А. Рындин, Г. В. Таньков, А. В. Затылкин // Цифровые модели в проектировании и производстве РЭС : межвуз. сб. науч. тр. / под ред. проф. Н. К. Юркова. - Пенза : Изд-во ПГУ, 2012. - Вып. 17. - 290 с.

5. Хайкин, С. Э. Физические основы механики / С. Э. Хайкин. - М. : Физматгиз, 1962. - 770 с.

6. Леонов, А. Г. Управление исследованиями моделей радиотехнических устройств на этапе проектирования / А. Г. Леонов, А. В. Затылкин, Н. К. Юрков // Прикаспийский журнал: управление и высокие технологии. - 2012. - № 1 (17). - С. 138-142.

7. Ольхов, Д. В. Система обработки экспериментальной информации в проектных исследованиях радиотехнических устройств / Д. В. Ольхов, А. В. Затылкин, Н. К. Юрков // Известия Южного федерального университета. Технические науки. - 2012. - № 5. - С. 94-99.

8. Затылкин, А. В. Моделирование изгибных колебаний в стержневых конструкциях РЭС / А. В. Затылкин, Г. В. Таньков // Надежность и качество : тр. Междунар. симп. / под ред. Н. К. Юркова - Пенза : Изд-во ПГУ, 2006. - С. 320-323.

Затылкин Александр Валентинович

кандидат технических наук, доцент, кафедра конструирования и производства радиоаппаратуры, Пензенский государственный университет E-mail: al.zatylkin@yandex.ru

Таньков Георгий Васильевич

кандидат технических наук, доцент, кафедра конструирования и производства радиоаппаратуры, Пензенский государственный университет E-mail: oldalez@yandex.ru

Ольхов Даниил Вадимович

студент,

Пензенский государственный университет E-mail: kipra@yandex.ru

Zatylkin Alexander Valentinovich

candidate of technical sciences, associate professor, sub-department of construction and production of radio equipment,

Penza State University

Tankov Georgiy Vasilievich

candidate of technical sciences, associate professor, sub-department of construction and production of radio equipment Penza State University

Ol'khov Daniil Vadimovich

student,

Penza State University

УДК 368.3.068 Затылкин, А. В.

Дискретная модель процесса распространения импульса смещения в упругом стержне постоянного сечения при торцевом ударе / А. В. Затылкин, Г. В. Таньков, Д. В. Ольхов // Вестник Пензенского государственного университета. - 2013. - № 4. - С. 79-85.