A modifiedBouc-Wen model to describe the hysteresis of non-stationary processes

Abstract


A number of known phenomenological models are considered, which are used to describe a variety of hysteresis effects in nature. In this case, the system is considered as a "black box" with known experimental values of input and output parameters. Correlations between them are established by mathematical functions, whose parameters are identified using experimental data. Amongthe phenomenological models there are marked the Bouc-Wen model andits analogsthat have been successfully usedin variousscientific and technical fieldsdue to the possibilityof the analytical description ofvarious hysteresis loopsof non-stationary processes. The conditions are formulatedwhich must be satisfied by the Bouc-Wen model. The main ones are the model adequacy of the physical process and stability. To describe the hysteresis, a mathematical model is suggested, according to which the force and kinematic parameters are bound by a special differential equation of the first order. In contrast to the Bouc-Wen model, the right side of this equation is chosen in the form of a polynomial of two variables determining the trajectory of a hysteresis in the process diagram. It is stated that this presentation provides the asymptotic approximation of the solution to the curves of the comprehending (including) hysteresis cycle.This cycle is formed by curves of direct and reverse processes ("loading-unloading" processes), which are based on experimental data for the maximum possible or permissible intervals of parameter changes during the steady vibrations. Coefficients in the right part are determined from experimental data for the comprehending hysteresis cycle under conditions of steady-state oscillations. Approximation curves of the comprehending cycle are constructed using the methods of minimizing the discrepancy of analytical representations to the number of experimental points. The proposed approach allows by one differential equation to describe the trajectory of hysteresis with a random starting point within the area of the comprehending cycle.

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

  1. Мayergoyz I.D. Mathematical models of hysteresis and their applications: 2nd edn. (Electromagnetism). - Amsterdam: Elsevier, Academic Press, 2003. - 498 p.
  2. Bertotti G. Hysteresis in magnetism: for physicists, materials scientists, and engineers, Part 2. - Boston: Elsevier, Academic Press, 1998. - 558 p.
  3. Rieger M.O. Young measure solutions for nonconvex elastodynamics // SIAM Journal on Mathematical Analysis. - 2003. - Vol. 34. - No. 6. - P. 1380-1398. doi: 10.1137/S0036141001392141
  4. Rieger M.O. A model for hysteresis in mechanics using local minimizers of Young measures // Progress in Nonlinear Differential Equations and Their Applications. - 2005. - Vol. 63. - P. 403-414. doi: 10.1007/3-7643-7384-9_39
  5. Mielke A. Analysis of energetic models for rate-independent materials // Proceedings of the International Congress of Mathematicians, Beijing, China. - 2002. - Vol. 3. - P. 817-828.
  6. Mielke A., Roubícek T. A rate-independent model for inelastic behavior of shape-memory alloys // Multiscale Modeling and Simulation. - 2003. - Vol. 1. - No. 4. - P. 571-597. doi: 10.1137/S1540345903422860
  7. Мишустин И.В., Мовчан А.А. Моделирование фазовых и структурных превращений в сплавах с памятью формы, происходящих под действием немонотонно меняющихся напряжений // Известия РАН. МТТ. - 2014. - № 1. - С. 37-53.
  8. Мишустин И.В., Мовчан А.А. Аналог теории пластического течения для описания деформации мартенситной неупругости в сплавах с памятью формы // Известия РАН. МТТ. - 2015. - № 2. - С. 78-95.
  9. The investigation on the nonlinearity of plasticine-like magnetorheological material under oscillatory shear rheometry / X. Gong, Ya. Xu, S. Xuan, C. Guo, L. Zong // Journal of Rheology. - 2012. - Vol. 56. - No. 6. - P. 1375-1391.
  10. Large amplitude oscillatory shear rheology for nonlinear viscoelasticity in hectorite suspensions containing poly(ethylene glycol) / Z. Tong, W.X. Sun, Y.R. Yang, T. Wang, X.X. Liu, C.Y. Wang // Polymer. - 2011. - Vol. 52. - No. 6. - P. 1402-1409.
  11. Danilin A.N., Yanovsky Yu.G., Semenov N.A., Shalashilin A.D. Kinematic model of the rheological behavior of non-Newtonian fluids in conditions of nonstationary cyclic loading // Composites: Mechanics, Computations, Applications: An International Journal. - 2012. - Vol. 3. - No. 4. - P. 1-15.
  12. Visintin A. Differential Models of hysteresis (Applied Mathematical Sciences). - Berlin: Springer, 1994. - 418 p.
  13. Nová I., Zemánek I. Аnalytical model with flexible parameters for dynamic hysteresis loops modeling //Journal of Electrical Engineering. - 2010. - Vol. 61. - No. 7. - P. 46-49.
  14. Danilin A.N., Shalashilin V.I. A method to identify hysteresis by an example of an antigalloping device // International Applied Mechanics. - 2010. - Vol.46. - No. 5. - P. 588-595.
  15. Preisach F. Über die magnetische Nachwirkung // Zeitschrift für Physik. - 1935. - P. 277-302.
  16. Preisach modeling of magnetization changes in steel / S.F.H. Parker, C.A. Faunce, P.J. Grundy, M.G. Maylin, J.L.C. Ludlow, R. Lane // Journal of Magnetism and Magnetic Materials. - 1995. - Vol. 145. - P. 51-56.
  17. Torre E.D. A Preisach model for accommodation // IEEE Transactions on Magnetics. - 1994. - Vol. 30. - No. 5. - P. 2701-2707.
  18. Красносельский М.А., Покровский А.В. Системы с гистерезисом. - М.: Наука, 1983. - 271 с.
  19. Smith R.Smart material systems: model development. - Philadelphia: SIAM, 2005. - 525 p.
  20. Leenen R. The modeling and identification of an hysteretic system. The wire-rope as a nonlinear shock vibration isolator // Eindhoven University of Technology: Dept. Mechanical Engineering. - DCT 2002.72, 2002.
  21. Rosensweig R.E.Ferrohydrodynamics. - New York: Dover Publications, Inc., 2014. - 348 p.
  22. Berkovski B., Bashtovoy V. Magnetic Fluids and Applications. - New York: Begell House Inc. Publishers, 1996. - 350 p.
  23. Bouc R. Forced vibrations of a mechanical system with hysteresis // Proceedings of the Fourth Conference on Nonlinear Oscillations, Prague, Czechoslovakia, 1967. - P. 315-321.
  24. Bouc R. Modèle mathématique d’hystérésis (A mathematical model for hysteresis) // Acustica. - 1971. - Vol. 21. - P. 16-25.
  25. Wen Y.K. Method for random vibration of hysteretic systems // Journal of the Engineering Mechanics Division. - 1976. - Vol. 102 (EM2). - P. 246-263.
  26. Wen Y.K. Equivalent linearization for hysteretic systems under random excitation // Journal of Applied Mechanics. - 1980. - Vol. 47. - P. 150-154.
  27. Development of adaptive modeling techniques for non-linear hysteretic systems / A.W. Smyth, S.F. Masri, E.B. Kosmatopoulos, A.G. Chassiakos, T.K. Caughey // International Journal of Non-Linear Mechanics. - 2002. - Vol. 37. - P. 1435-1451.
  28. Low T., Guo W. Modelling of a three-layer piezoelectric bimorph beam with hysteresis // IEEE Journal of Microelectromechanical Systems. - 1995. - Vol. 4. - No. 4. - P. 230-237.
  29. Yoshioka H., Ramallo J.C., Spencer Jr. B.F. “Smart” base isolation strategies employing magnetorheological dampers // Journal of Engineering Mechanics. - 2002. - Vol. 128. - No. 5. - P. 540-551.
  30. Foliente G.C. Hysteresis modelling of wood joints and structural systems // ASCE Journal of Structural Engineering. - 1995. - Vol. 121. - No. 6. - P. 1013-1022.
  31. Nagarajaiah S., Xiaohong S. Response of base-isolated USC hospital building in Northridge earthquake // ASCE Journal of Structural Engineering. - 2000. - Vol. 126. - No. 10. - P. 1177-1186.
  32. Ismail M., Ikhouane F., Rodellar J. The hysteresis Bouc-Wen model, a survey // Archives of Computational Methods in Engineering. - 2009. - Vol. 16. - P. 161-188.
  33. Ikhouane F., Mañosa V., Rodellar J. Dynamic properties of the hysteretic Bouc-Wen model // Systems & Control Letters. - 2007. - Vol. 56. - P. 197-205.
  34. Ikhouane F., Rodellar J. On the hysteretic Bouc-Wen model // Nonlinear Dynamics. - 2005. - Vol. 42. - P. 63-78.
  35. Данилин А.Н., Кузнецова Е.Л., Рабинский Л.Н. Модель гистерезиса энергорассеяния при колебаниях механических систем // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2014. - № 4. - С. 45-67.
  36. Danilin A.N., Vinogradov A.A., Lilien J.-L. Kinematic model for hysteretic dissipation of vibration energy for Torsional Damper and Detuner // Proceedings of 7th International Symposium on Cable Dynamics. Vienna (Austria), 10-13 December 2007. - 2007. - P. 247-253.
  37. Виноградов А.А., Данилин А.Н., Рабинский Л.Н. Деформирование многослойных проволочных конструкций спирального типа. Математическое моделирование, примеры использования / под ред. А.Н. Данилина. - М.: Изд-во Моск. авиац. ин-та, 2014. - 168 с.

Statistics

Views

Abstract - 224

PDF (Russian) - 93

Cited-By


PlumX


Copyright (c) 2016 Danilin A.N., Kuznetsova E.L., Kurdumov N.N., Rabinsky L.N., Tarasov S.S.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies