Управление процессом моделирования движения объекта с ударами о преграду

число услуг за счет передачи коротких сообщений в разных форматах от одного абонента другому. Кроме того, как правило, имеются шлюзы между однотипными сетями различных операторов, что позволяет организовать связь для их абонентов. Остается открытым вопрос сопряжения разнотипных сетей абонентского доступа, принадлежащих различным операторам. В этом случае шлюзовое оборудование должно принадлежать государству (федеральным или местным властям) и, конечно, его стоимость не должна быть слишком высокой. В качестве целевой функции здесь выступает либо минимизация приведенных затрат на это оборудование, либо среднее время передачи сообщения от одного абонента другому. Таким образом, можно воспользоваться методом решения задач размещения пунктов обслуживания в потоковых сетях [5] . Так как в основе подобных сетей лежат, как правило, разнотипные сети, к тому же принадлежащие различным собственникам, в качестве математических моделей в общем случае должны выступать нестационарные S-гиперсети.

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

На рисунке приведена блок-схема основного метода решения поставленных задач.

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

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

Литература

1. Портнов Э.Л. Принципы построения первичных сетей и оптические кабельные линии связи. М.: Горячая линия-Телеком, 2009. 543 с.

2. Соколов Н.А. Задачи планирования сетей электросвязи. СПб: Техника связи, 2012. 432 с.

3. Попков В.К. Математические модели связности. Новосибирск: ИВМиМГ СО РАН, 2006. 490 с.

4. Попков В.К. Применение теории S-гиперсетей для моделирования систем сетевой структуры // Проблемы информатики. 2010. № 4. С. 17-23.

5. Попков Г.В., Попков В.К. Об одном подходе размещения пунктов обслуживания в транспортных сетях // Проблемы информатики. 2012. № 3. С. 20-28.

References

1. Portnov E.L., Principy postroeniya pervichnykh setey i op-ticheskie linii svyazi, Moscow, Goryachaya liniya-Telekom, 2009.

2. Sokolov N.A., Zadachi planirovaniya setey elektrosvyazi, SPb, Tekhnika svyazi, 2012.

3. Popkov V.K., Matematicheskie modeli svyaznosti, Novosibirsk, Inst. of Сomputational Math. and Mathematical Geophys. of RAS (SB), 2006.

4. Popkov V.K., Problemy informatiki, 2010, № 4.

5. Popkov G.V., Popkov V.K., Problemy informatiki, 2012, № 3, pp. 34-39.

УДК 004.942

УПРАВЛЕНИЕ ПРОЦЕССОМ МОДЕЛИРОВАНИЯ ДВИЖЕНИЯ ОБЪЕКТА С УДАРАМИ О ПРЕГРАДУ

(Работа выполнена в рамках ФЦП «Научные и научно-педагогические кадры инновационной России»

на 2009-2011 гг., ГК№ П1122)

В.К. Манжосов, д.т.н., профессор, зав. кафедрой; Д.А. Новиков, ассистент кафедры

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

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

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

Предложены алгоритм управления процессом моделирования и элементы его реализации в проблемно-ориентированном программном комплексе при вычислительном эксперименте

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

SOFTWARE COMPLEX FOR MODELING OBJECT MOTION WITH IMPACT BARRIER Manzhosov V.K., Ph.D., Professor, Head of Chair; Novikov D.A., Teaching Assistant

Аbstract. The paper gives a review of the object motion model with a periodic force of the relay type and impact with the rigid barrier. Software complex provides a solution to the equations of motion taking into account the impacts, force discontinuous functions and chattering phenomenon. The mathematical modeling results generate numeric values base and they are reproduced in the form of the corresponding diagrams during the process of simulation. The motion animation sequence are implemented, the diagrams of displacement, velocity and object acceleration and motion phase diagram are reproduced. It is possible to implement turn-based a modeling process.

The ability to analyze the transitions and the motion limit cycle is provided. Statistical analysis of the numerical sequence pretonic speeds is performed to estimate the shock system to steady state of motion in addition to visual assessments. The analysis allows to implement a rigorous statement of the system output at steady state to a given level of deviation of numeric values that do not exceed a small value.

The implementation of the parallel computing of experiment and output of simulation results in the graphs of displacement, velocity and acceleration is provided to conduct a comparative analysis of simulation results of various experiments in the software product. This procedure is effectively used in the analysis of the motion stability.

Modeling management algorithm and its implementation elements in the problem-oriented software package for computational experiment are proposed.

Keywords: a program complex, modeling, impact, vibro-impact system.

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

Движение рассматриваемой системы опишем уравнениями

m

•х = р(t), х(i0) = x , X(i0) = x0

P (t) =

P1, (n-1 )T < t < ( n-l) T +11, P2, ( n-l) T +11 < t < n-T,

п=1, 2, ..., при х=хс, если х >0, то х+ =-Я х , где Р\, Р2 - силы, действующие на массу т соответственно в направлении ограничителя и в обратном направлении; Т - период силового воздействия; и - длительность действия силы Р\; п - номер цикла силового воздействия; х0 - координата массы в начальный момент времени при /=/0; х0 -скорость ударной массы в начальный момент времени; хс - координата ограничителя; х- - скорость ударной массы перед столкновением с ограничителем; х+ - скорость ударной массы после столкновения с ограничителем; Я - коэффициент восстановления скорости при ударе (0<Я< 1).

Методом припасовывания решений дифференциальных уравнений, описывающих смежные интервалы движения, разделенные моментом удара,

построены аналитические зависимости, определяющие параметры движения на различных интервалах [2]. Время нанесения первого удара ^ и

предударная скорость (х^))- могут быть найдены по зависимостям

t = t +-

Х0 + VX> 2A1

A

если tc < t,,

4 1 -

-x(t,)x2 (t,) + 24 (xc - x(t,)) tCi = t +--A---- , если tci >t,,

())- = 02 - 2А1 (Х0 - X ) , если К, < t1 ,

(х^))- = ^X2) + 2Л (хс -х«1)), если ^ > ^,

где А\=Р\/т, А2=Р2/т, х(^) и ) - координата и скорость ударной массы в момент переключения силы.

Если удар происходит до момента переключения силы t < ^, то удар тела об ограничитель может многократно повторяться на интервале I < ^ < . Если в первой фазе действия силы в момент времени ^ = tc происходит /-й удар, то на интервале t < t <

= (Х ('с-, ))+ + А "(' - 'с-) , Х() = Хс + (X (^ ))+  ( - ^ ) + -2 А . (* - )2 , где - время нанесения (/-1)-го удара;

- скорость ударной мас-

сы после нанесения (г'-1)-го удара, являющаяся начальной скоростью ее движения на интервале

t < ? < I

С-1 с

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

2R ,., чч-

t = t +-

4 А1 (1 - R) ^ .

Если 1С </ь то многократный ударный режим

движения ударной массы заканчивается еще в первой фазе действия силы Р(/) и следует переходить к рассмотрению движения во второй фазе действия силы при начальных условиях х(/1)=хс,

) = 0 .

Если неравенство 1С < /1 не выполняется, возникает необходимость определения конечного числа ударов j массы об ограничитель до начала второй фазы действия силы Р(/):

t - tc ln[1 -J—^ (R -1)]

J =-

ln R

-+1, a, =1- R • (xc(tC]))- .

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

причем (х(4))- = Я'-1 •(х))-. Так как Д<1, то Я'— при г'^да стремится к нулю. Если рассмотреть

отношение

Х- (О'Х- (',) =

то отношение

скорости удара на '-м соударении к предударной скорости первого удара при достаточно большом числе соударений становится малой величиной. Если задаться величиной этой малости е и полагать, что

если R' 1 >е,

10, если R' <s,

то можно определить минимальное число ударов /т1П, после которого следует считать, что ударная масса находится в покое у ограничителя, пока

/шт = 1п е /1п Я + 1.

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

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

Для численного решения дифференциальных уравнений движения с учетом соударений, разрывных функций силового воздействия на объект, явления дребезга разработан программный продукт (Свид. о гос. регистр. № 2010614035 от 22.06.10, авторы: Манжосов В.К., Новиков Д.А., Корняков Д.Е., Муромцев И.И.).

Блок-схема программного комплекса представлена на рисунке 1.

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

В работу программного комплекса заложена возможность его использования с несколькими типами БД - автономными и клиент-серверными. Обеспечиваются оба типа соединения посредством доступа к БД через Microsoft ActiveX Data Objects (ADO), связь с данными в котором осуществляется посредством технологии OLE DB.

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

Данная возможность реализована для построения рабочих моделей эксперимента вне зависимости от используемой аппаратной платформы как на локальных ПЭВМ, так и на объединенных в локальную сеть с выделенным сервером БД.

В случае использования клиент-серверной технологии регистрации экспериментов можно отойти от проблемы синхронизации наработанных БД. Используя автономную БД, представленную в виде БД Microsoft Access 2003, можно использовать программный продукт без предварительной подготовки ПЭВМ (за исключением случаев отсутствия на ЭВМ программного комплекса Microsoft Office, что встречается довольно редко).

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

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

а

Рис. 1

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

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

При активировании кнопки «Старт» происходит запуск таймера и реализуется процедура численного расчета дифференциальных уравнений движения. Задержка анимации позволяет регулировать скорость отображения процесса моделирования. Обеспечивается возможность изменения величины временного шага при расчете.

Отображение окна фазовой диаграммы происходит при выборе пункта «Отображать». Фазовая диаграмма представляет собой геометрическое место точек в координатной плоскости «скорость-координата». Масштабы оси абсцисс и ординат настраиваются через вкладку «Фазовая диаграмма». Здесь также можно настроить размеры отображаемого окна.

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

Важным моментом анализа функционирования системы является определение максимальной ско-

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

Блок отображения диаграмм координаты, скорости и ускорения ударника предусматривает возможность графического отображения параметров эксперимента.

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

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

Для сравнительного анализа результатов моделирования различных экспериментов в программном продукте существует вкладка «Фоновая модель», отвечающая за параллельный расчет эксперимента и вывод данных результатов моделирования для воспроизведения в виде диаграмм координаты, скорости и ускорения.

В программном продукте разработан блок табличного отображения числовых значений результатов моделирования. Настройки на вкладке «Таблица» основного экрана программы позволя-

ют выбирать необходимые к снятию параметры, а также отображать форму «Таблица параметров».

На самой форме «Таблица параметров» присутствуют дополнительные настройки записи значений параметров в табличную часть формы. Настроечное поле «Пропуск записей» позволяет регулировать частоту снятия записей в режимах разгона и торможения ударника, а «Критическая координата» указывает, при достижении какого значения координаты ударника начнется запись параметров на каждом расчетном шаге. Форма поддерживает копирование данных в OLE Container и вставку данных в любом из приложений, поддерживающих данный стандарт.

Общий вид интерфейса программного комплекса для моделирования процесса движения ударной системы при столкновении с преградой представлен на рисунке 3.

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

Непосредственно под набором вкладок, регулирующих настройки моделирования, находится панель управления БД экспериментов. Она позволяет добавлять, удалять, сохранять изменения и производить навигацию по хранящейся базе проведенных экспериментов. Внизу панели управления расположены кнопки «Старт», «Стоп», «Очистка».

Активируя курсором кнопку «Старт», пользователь осуществляет запуск процесса моделирования. На экране монитора реализуется анимационный процесс движения ударной массы, воспроизводятся диаграммы положения х, скорости v и

ускорения А (по желанию пользователя) ударной массы, а также фазовая диаграмма у=у(х) движения ударной массы.

Возможна пошаговая реализация процесса моделирования при помощи кнопки «Пошагово».

Отдельное внимание следует уделить вкладке «Статистика». Она используется для построения модели силового воздействия на ударник на основе экспериментальных данных, что позволяет оценивать достоверность полученной модели по реальным данным эксперимента.

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

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

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

Устойчивость процесса движения анализируется путем оценки сходимости ряда в виде числовой последовательности разницы предыдущего и текущего значений предударных скоростей с использованием критерия Коши. Сформулируем этот критерий для рассматриваемого процесса:

п

чтобы ряд суммы чисел (^ (V - У1_1) + ...) сходил-

0

ся, необходимо и достаточно, чтобы последовательность его частичных сумм ^, 53, ..., обладала следующим свойством: каково бы ни было е>0, существует такое п, что при любом т>0

К+» - <е.

При вычислительном эксперименте делаются выборка чисел (V-V«) и расчет частичных сумм для пе(50, 100, 200, 300, 500, 1000, 5000 и т.д.), где п - число циклов. Чем меньше значение малой величины е>0, тем большее число циклов необходимо для гарантированной констатации устойчивости процесса. Моделирование показало, что в

большинстве случаев при е<0,001 можно ограничиться выборкой для 100<n<500.

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

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

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

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

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

Литература

1. Крупенин В.Л. Ударные и виброударные машины и устройства // Вестн. науч.-технич. развития. Национальная технологическая группа. 2009. № 4 (20). С. 3-32.

2. Манжосов В.К., Новиков Д.А. Моделирование режимов движения виброударной системы при периодическом силовом воздействии // Изв. Саратов. ун-та: Сер. Математика. Механика. Информатика. 2010. Вып. 4. Т. 10. С. 65-71.

References

1. Krnpenin V.L., Vestnik nauchno-tekhnicheskogo razvitiya, 2009, no. 4 (20), pp. 3-32.

2. Manzhosov V.K., Novikov D.A., Izvestiya Saratovskogo Univ., 2010, vyp. 4, Vol. 10, pp. 65-71.