О построении расчетных моделей динамических систем по результатам испытаний Текст научной статьи по специальности «Физика»

СПИСОК ЛИТЕРАТУРЫ

1. Гихман И.И., Скороход А.В., Ядренко М.И. Теория вероятностей и математическая статистика. - Киев: Вища школа, 1979. - 408 с.

2. Золотарев В.М. Одномерные устойчивые распределения. - М.: Наука, 1983. - 304 с.

3. Тимофеев В.С. Оценивание параметров регрессионных зависимостей на основе характеристической функции // Научный вестник НГТУ. - 2010. - № 2 (39). - С. 43-52.

4. Гнеденко Б.В. Курс теории вероятностей. - М.: Едиториал УРСС, 2001. - 320 с.

5. Feuerverger A., Mureika R.A. The empirical characteristic function and its applications // The annals of statistics. - 1977. - V. 5. -№1.- P. 88-97.

6. Paulson A.S., Holcomb E.W., Leitch R.A. The estimation of parameters of the stable laws // Biometrica. - 1975. - № 62. -P. 163-170.

7. Koutrouvelis I.A. Regression-type estimation of the parameters of stable laws // Journal of the American statistical association. 1980. -V. 75. - № 372. - P. 918-928.

8. Дрейпер Н., Смит Н. Прикладной регрессионный анализ. -М.: Статистика, 1973. - 392 с.

9. Оппенгейм А.В., Шафер Р.В. Цифровая обработка сигналов. -М.: Связь, 1979. - 416 с.

10. Ивахненко А.Г., Степашко В.С. Помехоустойчивость моделирования. - Киев: Наукова думка, 1985. - 216 с.

Поступила 04.10.2010 г.

УДК 629.7.018.4:534.1.015.1

О ПОСТРОЕНИИ РАСЧЕТНЫХ МОДЕЛЕЙ ДИНАМИЧЕСКИХ СИСТЕМ ПО РЕЗУЛЬТАТАМ ИСПЫТАНИЙ

В.А. Бернс

ФГУП «Сибирский научно-исследовательский институт авиации им. С.А. Чаплыгина», г. Новосибирск

E-mail: v.berns@yandex.ru

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

Ключевые слова:

Расчетная динамическая модель, обобщенные динамические характеристики, резонансные испытания. Key words:

Calculated dynamic model, generalized dynamic characteristics, resonance testing.

Введение

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

Методика построения расчетной динамической

модели по результатам испытаний

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

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

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

ента демпфирования /у-го собственного тона колебаний используются следующие формулы [2, 3]:

а=у/-Ц)/[(1+//)(р/-Щ^*2], /=у (ц+/)/[(1+/' (1)

Здесь и и V векторы, компонентами которых являются действительные и мнимые составляющие перемещений точек системы; Ц, Ц - векторы синфазной и квадратурной составляющих сил возбуждения, обеспечивающие выделение -го тона колебаний на частоте щ; / -котангенс сдвига фаз между перемещениями и синфазной составляющей возбуждения; р - собственная частота тона (частота фазового резонанса), определяемая по переходу / через 0; V'*- мнимая составляющая перемещения точки нормирования формы тона; _/=1,2,...,Д N -число исследуемых собственных тонов. Если в (1) положить /=0 при (условие «фиктивного» фазового резонанса), то получим известную формулу определения обобщенных масс введением квадратурной составляющей возбуждения [2].

Оценка погрешностей представления динамической

системы расчетной моделью

Рассмотрим построение математической модели с N степенями свободы системы, имеющей М степеней свободы, причем ЖМ. Воспользуемся методикой идентификации, разработанной в [3]: для анализа свойств динамической системы и определения характеристик ее собственных тонов используются свойства вынужденных монофазных колебаний. Вынужденные колебания динамических систем являются монофазными (аналогично и монофазное возбуждение), если различия в фазах колеблющихся точек равны 0 или л [5].

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

ир=/{Щ уот

Компонентами векторов ир и являются, соответственно, величины действительных и мнимых составляющиху'-х монофазных колебаний в N точках возбуждения колебаний. Индекс N означает, что данный параметр определяет монофазные колебания в этих N точках (далее везде верхним индексом N будем отмечать величины, относящиеся к модели с N степенями свободы). Исходная система имеет М степеней свободы. Действительные и мнимые составляющие вынужденных колебаний системы при гармоническом возбуждении определяются уравнениями

—ЩА" Vм) +Н") У") +СМ) Л" =ЕМ) ЩА<") У" +Н(М) иМ)—СМ) У") =ЦМ) поэтому величина /(Л) есть

//"> = у<МТ(-ЩЛ(М)+С(М)) у(М)/[ у(М)ТН(М)у(М)], (2)

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

Для упрощения выкладок будем считать, что первые N компонентов вектора у(М) составляют вектор уч; а<"), С(М), Н" - матрицы инерции, жесткости и демпфирования исходной системы.

Преобразуем (2), выделив при этом точки возбуждения колебаний

( =

I (-

2 Л М) ^ ЛМ Ьл,(Ю ),,(Ю )

щ а

)() у(Ю) +

I (-

2„(М) )\4МК1М)

щ а,

А "А

^( М),,(N )< N)

"¿"г + I нм)Vм)"(М)

у—1 / )=ю+1

"а ур У]1 + I "и ( "А

(3)

Кроме того, представим вектор у(М) в виде разложения по собственным векторам системы без демпфирования

у{М)=Жц.

Здесь Ж - матрица собственных векторов, ду -коэффициенты разложения.

С учетом ортогональности собственных векторов и условия их нормировки

1=1,2.,N.

Получим

I а(М)( рМ )2-щ2)^ +

/—1

М

+ I а(М)(р(М)2

/Ю) —

N М

IНМ"8р8]1 + I НМ)^

А

/ ,1

/1 —N +1

(4)

Здесь /;;(М) — элементы матрицы демпфирования в главных координатах. Первыми N столбцами матрицы Жявляются собственные векторы учитываемых тонов.

Формулы (3) и (4) отражают тот очевидный факт, что при N=М величина /т совпадает с точным значением. При ЖМ точность вычисления //щ зависит от числа точек возбуждения и регистрации колебаний, способа распределения сило-возбудителей по конструкции и определяется соотношением между величинами собственных частот, распределением масс, жесткостей и демпфирования. Так как распределение упруго-массовых характеристик и характеристик демпфирования неизвестно до проведения резонансных испытаний, то относительно размещения силовозбудите-лей можно дать лишь общие рекомендации, которые можно найти, например, в работе [6].

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

УМ=Щ+Г, (5)

причем компоненты вектора Гу равны нулю при /=1,2,., N. В этом случае величина использует-

ся для вычисления обобщенных масс и обобщенных коэффициентов демпфирования по формулам (1). Кроме того, по переходу Я^ через нуль определяются собственные частоты системы в предположении, что при 1/^=0 монофазные колебания являются собственными.

Выражение (3) с учетом (5) можно представить в виде

Я(ю =

Л М )(„(М )2

]

( р(М )2-о2) 8 2- +

I (-

2 Л М )

со а,

„(м)

ЬМ 8- + I НМ )(2Г-1у/! + У, у )

(6)

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

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

Для немонофазного возбуждения выражение (6) принимает следующий вид:

где к!М) - обобщенный коэффициент демпфирования у-го тона;

grh¡M)ЩТE¡M)/[a¡M2(p¡M2-CУ+h¡M2]. (7)

Первыми N компонентами вектора Е<М) являются компоненты вектора

Е (Ю) = -,

а остальные равны нулю. Вектор определяется из задачи о собственных значениях [3]

(ф -XV )С = 0.

Здесь и и V - есть матрицы, столбцами которых являются действительные и мнимые составляющие перемещений точек системы, соответствующие N линейно независимым силам возбуждения (матрица Е).

На основании выражения (6) можно сделать выводы о том, что точность определения ЯЯт, а значит и обобщенной массы и обобщенного коэффициента демпфирования -го тона, зависит от соотношения между частотой вынужденных колебаний и собственной частотой тона. При отсутствии близких частот она повышается с приближением о кр£т. С изменением числа степеней свободы модели изменяется точность расчета Я/4, собственной частоты, обобщенной массы и коэффициента демпфирования у-го тона. Если величины компонент вектора Гу определяются влиянием высших тонов, то Ц>\Я1М\ при ю<р(М) и Я\<\Я(М)\ при ет>ру(М). Это означает, что собственная частота тона р^, определяемая из условия Я(Л,=0, будет выше соответствующего точного значения, т. е. при N^■M-p¡N) стремится кру(М) сверху. Величина обобщенного коэффициента демпфирования ниже ЦМ) при оз<р М и выше ЦМ) при о>ру(М). Если на исследуемый тон преобладающее влияние оказывают низшие тона, то приведенные соотношения меняются на противоположные.

На точность определения обобщенных динамических характеристик оказывает влияние уровень демпфирования в системе. Изменение демпфирования одного порядка по всем тонам оказывает наибольшее влияние, как следует из (6) и (7), на величины Я тех тонов, собственные частоты которых находятся в исследуемом диапазоне частот. С уве-

Я( ю) =

]

,(м)(-)2-с2) 82

У(ы)

1«+

+ ! (-®2аМ + СМ^уу

I ЬМ^ -)

(8)

где у, /=1,2,..., N - компоненты вектора мнимой составляющей возбуждения колебаний по -му тону на частоте со.

Величина ЯЯМ) зависит от соотношения уровней действительной и мнимой составляющих возбуждения и может быть выбрана любой. В действительности реализуемым значением ЯЯМ) является Яр, определяемая по (8), которая отличается от Я/М) вследствие несоответствия степеней свободы системы и ее математической модели. Относительно соотношений между Я/4 и ЯМ можно сделать те же выводы, что и для случая монофазного возбуждения, если вместо собственной частоты тона оперировать частотой фиктивного фазового резонанса.

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

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

СПИСОК ЛИТЕРАТУРЫ

1. Галкин М.С. Вопросы идентификации в аэроупругости // Ча-плыгинские чтения: сб. докл. / под ред. В.Г. Сувернева. - М.: ЦНТИ Волна, 1983. - С. 84-97.

2. Смыслов В.И. Некоторые вопросы методики многоточечного возбуждения при экспериментальном исследовании колебаний упругих конструкций // Ученые записки ЦАГИ им. Н.Е. Жуковского. - 1972. - Т. 3. - № 5. - С. 110-118.

3. Бернс В.А. Модальная идентификация динамических систем на основе монофазных колебаний // Научный вестник НГТУ. - 2010. - № 3 (40). - С. 99-109.

4. Хейлен В., Ламменс С., Сас П. Модальный анализ: теория и испытания. - М.: ООО «Новатест», 2010. - 319 с.

5. Кононенко В.О., Плахтиенко Н.П. Методы идентификации механических нелинейных колебательных систем. - Киев: Наукова думка, 1976. - 114 с.

6. Смыслов В.И. Об экспериментальных способах исследования колебаний летательных аппаратов // Труды ЦАГИ им. Н.Е. Жуковского. - 1970. - Вып. 1217. - С. 3-63.

7. Бабаков И.М. Теория колебаний. - М.: Наука, 1968. - 560 с.

Поступила 22.11.2010 г.

УДК 621.313:534

ДЕСТРУКЦИЯ ТЕЛ СКОЛЬЗЯЩЕГО КОНТАКТА

М.С. Слободян, С.М. Слободян

Томский политехнический университет E-mail: sms1@tpu.ru; s.m.s@ngs.ru

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

Ключевые слова:

Стохастическая динамика, электрический контакт, эволюция, деструкция. Key words:

Stochastic dynamics, electrical contact, evolution, destruction.

Введение

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

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

Для процесса скольжения характерно то, что области структур поверхностных слоёв тел, находя-

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

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

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