Full Text

Введение Особенности гистерезиса (форма петлеобразных траекторий, асимптотическое подобие или симметрия прямого и обратного процессов и др.) определяются природой процесса [1-11]. Поэтому для математического описания гистерезиса авторы, как правило, предлагают модели, опираясь на конкретные представления о физическом явлении. Однако, несмотря на специфику моделей, ряд из них может претендовать на определённую общность, выходящую за область специализированного исследования. Феноменологический подход является весьма эффективным инструментом к созданию обобщенных математических моделей. В современной литературе имеется достаточно большое число примеров, демонстрирующих «миграцию» идеи метода из одной области исследования в другую. Это свидетельствует не только о «глубине» математических моделей, но и, в большей степени, об общих чертах различных по природе гистерезисных процессов. Поэтому представляется целесообразным кратко остановиться на нескольких известных феноменологических моделях, которые успешно применяются для описания гистерезиса в различных областях науки и техники. В динамике сложных механических систем, составные элементы которых взаимодействуют друг с другом силами различной природы, построение теоретических моделей с идентификацией их параметров представляет собой непростую задачу. В таких случаях механическую систему можно рассматривать как «черный ящик» с известными из эксперимента значениями входных и выходных параметров. Взаимосвязи между этими параметрами устанавливаются на основе феноменологических моделей, параметры которых идентифицируются с использованием экспериментальных данных [12-14]. В этом направлении особое место занимают модели, которые строятся с использованием спектральных разложений по релейным нелинейностям. Такой подход был предложен в 1935 г. немецким физиком Ф. Прейзахом в работах по магнетизму [15-17], в которых процесс намагничивания рассматривается как статистический результат перемагничивания отдельных элементарных областей (доменов). Считается, что каждая такая область может находиться только в состоянии насыщения с направлением намагниченности вдоль или против действия внешнего поля. Соответственно этому намагниченность каждой области описывается с помощью функций-переключателей, определяющих петлю гистерезиса в виде прямоугольника. Важной составляющей модели Ф. Прейзаха является функция распределения ориентаций доменов, с помощью которой определяются значения намагниченности в произвольном поле. В настоящее время идеи Ф. Прейзаха превращены в строгий математический аппарат и существенно развиты в работах М.А. Красносельского, А.В. Покровского и их последователей [18]. Схожие феноменологические представления предложены и развиваются в различных областях механики и физики [19-22]. Однако для идентификации параметров таких моделей часто требуются сложные экспериментальные исследования и интерпретации полученных данных. В 1967 г. Бук предложил способ решения задачи о вынужденных колебаниях механической системы с гистерезисом зависимости восстанавливающей силы от перемещения [23]. Траектория гистерезиса описывалась с помощью нелинейного обыкновенного дифференциального уравнения первого порядка, для которого идентификация коэффициентов осуществлялась с использованием опорных экспериментальных данных. В 1971 г. Бук представил модель гистерезиса уже для абстрактной физической системы, рассматривая её как «чёрный ящик» с известными данными на её входе и выходе [24]. В 1976 г. и далее модель была обобщена в работах Вена [25, 26], и с тех пор она известна как модель Бук-Вена [12]. 1. Модель Бук-Вена В настоящее время дифференциальная модель Бук-Вена привлекает повышенный интерес благодаря возможности аналитического описания разнообразных по форме гистерезисных петель, возникающих в различных системах с демпфированием [27]. В частности, модель успешно была использована для моделирования гистерезиса пьезоэлектрических элементов [28], магнитореологических амортизаторов [29], деревянных соединений [30], изоляции фундаментов зданий и сооружений [31] и пр. В современной литературе модель Бук-Вена используется, как правило, в рамках подхода «черного ящика» [12, 32]. Рассмотрим систему с гистерезисом, преобразующую входной сигнал q, зависящий от времени t, в сигнал f на выходе. В соответствии с моделью Бук-Вена отображение устанавливается с использованием нелинейного обыкновенного дифференциального уравнения, в состав которого входят функции-переключатели направления процесса. Общей формой записи модели Бук-Вена является дифференциальная связь между входным и выходным сигналами в виде [32] , (1) где g - выбранная кусочно-гладкая функция, идентификация которой осуществляется по экспериментальным данным опорного сигнала. Уравнение (1) включается в общую систему уравнений динамики системы, содержащих и в качестве неизвестных функций. Классическим примером использования (1) является уравнение Бук-Вена [25, 26], описывающее динамику механической системы с гистерезисом диссипации энергии. Для одномерного осциллятора эти уравнения имеют вид [33, 34] (2) Здесь сумма представляет собой силу демпфирования , где c - коэффициент вязкого демпфирования; k - коэффициент жесткости; z - переменная, описывающая гистерезисную траекторию; - постоянная составляющая . Значения , , , , , , и идентифицируются по экспериментальным данным. В литературе описаны различные модификации и усложнения (2), позволяющие описывать с определённой точностью гистерезис разнообразных по природе физических процессов [12, 27-37]. В результате интегрирования уравнений динамики и (1) строится зависимость , определяющая кусочно-гладкую непрерывную траекторию гистерезиса. Значения q, когда , формируют последовательность точек , где k - порядковый номер при возрастании t от начала развития процесса. При переходе через эти точки производная последовательно меняет знак и происходит смена ветви гистерезисного процесса, как это показано на рис. 1. Жирные стрелки во втором и четвертом квадрантах показывают соответственно прямое и обратное направления гистерезисного процесса, соответствующие росту или убыванию , когда или . Применение методов идентификации позволяет определить параметры модели так, чтобы ошибка (расхождение) между выходными данными, полученными из эксперимента, и вычислениями по алгоритмам модели была достаточно мала. Расчёты осуществляют для заданного (опорного) входного сигнала. После этого модель используется для моделирования гистерезиса при других входных сигналах. Однако известны примеры, когда найденные параметры модели Бук-Вена не обеспечивают соответствие результатов вычислений данным, полученным из экспериментов для иных входных сигналов. Такие примеры говорят о неоднозначности идентификации, которая может приводить к неустойчивости модели относительно входного сигнала. В литературе сформулированы условия, которым должна удовлетворять модель Бук-Вена [12, 32-34]. Основными являются адекватность математической модели физическому процессу и её устойчивость. Модель считается устойчивой по входу, если из ограниченности сигнала на входе системы следует ограниченность сигнала на её выходе при любых начальных условиях (bounded input, bounded output (BIBO) stability). Условия устойчивости накладывают ограничения на диапазоны изменения параметров модели. Однако выбор параметров из области устойчивости не всегда обеспечивает адекватность модели Бук-Вена реальному физическому процессу. Математические конструкции модели, подобные (2), также не всегда гарантирует такую адекватность. Рис. 1. Гистерезисная траектория внутри объемлющего цикла Fig. 1. Hysteresis trajectory inside the surrounding cycle Настоящая статья является развитием кинематического подхода, предложенного в работах [35-37]. В соответствии с этим подходом физическая связь между параметрами системы, описывающими процесс гистерезиса, устанавливается с использованием дифференциального уравнения (1) с правой частью в виде полинома от двух переменных и . Неизвестные коэффициенты полинома определяются по экспериментальным данным для кривых объемлющего цикла, область которого содержит любую возможную траекторию гистерезиса. На рис. 1 объемлющий цикл показан двумя пунктирными кривыми и , соответствующими прямому и обратному процессам. Вычислительные эксперименты показали, что предложенная модель обладает асимптотической устойчивостью и позволяет описывать сложные гистерезисные траектории адекватно реальному процессу. 2. Свойства гистерезисного процесса Как известно, для динамических систем с гистерезисом зависимости между силовыми и кинематическими параметрами имеют циклический характер. На диаграммах гистерезиса траектория каждого цикла имеет форму петли, которая образуется двумя кривыми (ветвями), соответствующими росту или убыванию параметра процесса . Начальная точка каждой ветви определяется предысторией нестационарного процесса. Она может находиться в любой точке пространства изменений исследуемых параметров. В условиях неустановившихся колебаний гистерезисные петли могут быть замкнутыми и незамкнутыми, отличаться друг от друга как по форме, так и относительному расположению. Опыты показывают, что для широкого круга гистерезисных систем однотипные ветви локальных циклов асимптотически приближаются к соответствующим кривым объемлющего цикла при монотонном изменении параметра процесса. В настоящей работе для описания подобных гистерезисных процессов в рамках дифференциального подхода в виде (1) используется физическая модель, удовлетворяющая нижеследующим свойствам. Свойство 1. Считается, что любая физически возможная траектория гистерезиса принадлежит области объемлющего цикла, который строится по экспериментальным данным для максимально допустимого диапазона изменения параметра процесса . Свойство 2. Гистерезисный процесс рассматривается как частотно-независимый. В этом случае гистерезисные кривые объемлющего цикла можно получить на основе квазистатических испытаний. Свойство 3. Считается, что все локальные кривые прямых процессов асимптотически подобны кривой , так же как и все локальные кривые обратных процессов асимптотически подобны . Это означает, что при монотонном изменении параметра q каждая локальная кривая стремится к соответствующей кривой объемлющего цикла. Например, на рис. 1 кривая с начальной точкой стремится с уменьшением q к кривой , кривая с начальной точкой стремится с возрастанием q к кривой и т.д. Причём все локальные кривые, продолженные при или , сходятся в точках и согласно направлению процесса. Предлагаемый подход заключается в использовании обыкновенного дифференциального уравнения первого порядка в виде (1), устанавливающего зависимость для каждой ветви гистерезисной траектории. Правая часть уравнения подбирается из класса функций, обеспечивающих асимптотическое приближение решения к соответствующим кривым или объемлющего цикла. В этом случае удаётся описать бесконечное множество ветвей прямого и обратного процессов, имеющих разные точки «старта», но приближающиеся к кривым или в зависимости от направления процесса. 3. Дифференциальное уравнение гистерезиса Для моделирования ветвей гистерезисной траектории используется обыкновенное дифференциальное уравнение первого порядка с правой частью в виде полинома от двух переменных и : . (3) Коэффициенты в (3) определяются методами приближения с использованием экспериментальных данных для двух кривых (ветвей) объемлющего цикла. Поэтому матрица имеет два набора значений. Будем обозначать их соответственно как и . Выбор между ними осуществляется по правилу (4) Уравнение (3) должно быть присоединено к уравнениям движения системы для интегрирования по времени при заданных начальных условиях. Для этого производная преобразуется к виду и подставляется в левую часть (3). Последующее умножение на приводит к дифференциальному соотношению , (5) которым дополняются уравнения движения системы. Знак определяет направление гистерезисного процесса, а также набор констант в соответствии с (4). В силу (4) прямой и обратный процессы описываются разными уравнениями. В точках , когда , происходит смена одного уравнения на другое, однако непрерывность не нарушается. Эта смена чередуется в процессе интегрирования по времени. Уравнения (5), (4) можно объединить, если воспользоваться функцией-переключателем , где - скорость. Тогда вместо (5) и (4) можно записать одно уравнение , (6) которым моделируются как прямой, так и обратный процессы. Уравнение (6) является частной формой записи уравнения (1). Во многих случаях можно считать, что прямой и обратный процессы протекают одинаково, т.е. отличаются лишь направлением изменения параметра процесса . В этом случае кривые объемлющего цикла симметричны относительно начала координат, как показано на рис. 2. На диаграмме произвольной точке с координатами кривой прямого процесса отвечает точка с координатами кривой обратного процесса и наоборот. Следовательно, обратный процесс описывается уравнением прямого процесса -fk f q Рис. 2. Ветви цикла, симметричные относительно начала координат Fig. 2. Cycle branches symmetrical with regard to the origin (7) с переменными , Обратный переход к q и M вместо (7) даёт (8) Уравнения прямого (3) и обратного (8) процессов можно объединить и записать аналогично (6) в виде . (9) Способ определения коэффициентов описан в работах [11, 35]. Идентификация коэффициентов осуществляется с использованием только экспериментальных данных, соответствующих кривым прямого и обратного процессов объемлющего цикла для установившегося процесса. Таким образом, решается проблема выбора опорной гистерезисной траектории для идентификации параметров модели. Вычислительные эксперименты показали асимптотическую устойчивость модели и её адекватность реальным процессам со сложной формой кривых объемлющего цикла [11, 35-37]. В качестве примера на рис. 4 показаны результаты моделирования переходного процесса колебаний гасителя маятникового типа (рис. 3) при внешнем гармоническом воздействии с амплитудой . Конструкция демпферного узла гасителя состоит из одного или нескольких ведущих дисков 1, одного или нескольких ведомых дисков 2, скрепленных с маятником. Диски имеют центральную ось, позволяющую им вращаться друг относительно друга. Диски имеют специальные профильные углубления, между которыми помещаются и частично зажимаются эластомерные элементы 3, контактирующие с ведущими и ведомыми дисками. Движение конструкции через систему креплений передаётся к ведущим дискам гасителя. В силу инерционности маятника происходит прокручивание ведущих дисков относительно ведомых. При этом происходят перемещения эластомерных элементов по каналам, образуемым криволинейными углублениями на поверхностях дисков. При этом перемещение шариков сопровождается их значительным деформированием и трением о поверхности каналов, что является причиной гистерезиса. Подробные описания эксперимента, вариантов конструкции и способа идентификации параметров модели изложены в работах [14, 35]. б а б а 1 2 31 G I Рис. 3. Гасители колебаний маятникового типа: а - схема работы; б - варианты конструкции (I - инерционный, G - гравитационный) Fig. 3. Pendulum-type vibration suppressors: a - operational scheme; b - structure variants (I is inertial and G is gravitational ones) Траектории на рис. 4 отвечают соответственно колебаниям «инерциального» гасителя при внешнем гармоническом воздействии с частотой (рис. 4, а) и колебаниям «гравитационного» гасителя при воздействии с частотой (рис. 4, б). Красным цветом отмечены кривые прямого и обратного процессов объемлющего цикла, в пределах которого строится траектория гистерезиса. б а -150 -100 -50 0 50 100 150 -0,4 0 0,4 -0,4 0 0,4 -80 -40 0 40 80 Рис. 4. Гистерезис зависимостей момента сопротивления М от угла поворота штанги гасителя колебаний маятникового типа Fig. 4. Hysteresisof dependences between modulus of resistance М and angle of rotation of the rod of pendulum-type vibration suppressor Заключение В работе развивается кинематический подход, предложенный ранее авторами для описания гистерезиса различных нестационарных процессов. Определена связь подхода с известным методом Бук-Вена. Так же как и модель Бук-Вена, математическая модель строится на основе обыкновенного дифференциального уравнения первого порядка. Однако в отличие от модели Бук-Вена правая часть этого уравнения выбирается в виде полинома от двух переменных - параметра гистерезисного процесса, зависящего от времени, и функции этого параметра. Такими переменными являются, например, перемещение (угол поворота) и соответствующая ему сила (момент). Коэффициенты полиномиального представления отыскиваются в результате аналитической аппроксимации кривых объемлющего (включающего) цикла в условиях установившихся колебаний. Предполагается, что такое представление обеспечивает асимптотическую устойчивость решения, это подтверждается численными экспериментами и сравнениями с экспериментальными данными, опубликованными авторами ранее. Полиномиальные коэффициенты правой части определяются методами приближения, минимизируя невязку аналитического представления к множеству значений, полученных в экспериментах. Использование кривых объемлющего цикла позволяет преодолеть две проблемы, свойственные подходу Бук-Вена. Одна из них связана с однозначным выбором опорного сигнала (опорных кривых) для идентификации параметров модели, вторая - с адекватностью модели реальному процессу. Предложенная модель физических зависимостей гистерезисного типа является аналитической, что является важным для её включения в модели, описывающих нелинейное поведение механических и иных систем. Разработанный подход может быть использован при решении различных задач о нестационарных колебаниях различных конструкций и механизмов с гистерезисным характером рассеяния энергии. Подход может оказаться полезным и при решении иных задач механики, в том числе об упругопластическом циклическом деформировании различных материалов и конструкций.

About the authors

A N Danilin

Institute of Applied Mechanics RAS; Moscow Aviation Institute

E L Kuznetsova

Moscow Aviation Institute

N N Kurdumov

Institute of Applied Mechanics RAS; Moscow Aviation Institute

L N Rabinsky

Moscow Aviation Institute

S S Tarasov

Moscow Aviation Institute

References

Statistics

Views

Abstract - 242

PDF (Russian) - 98