MDAMP

[fancysim 2]

Я еще относительный новичок в ansys. Провожу динамические расчеты в воркбенч (механика).

 

У меня есть результаты модальных испытаний моего объекта, а именно собственные частоты и значения добротностей (коэффициентов демпфирования на этих частотах). Знаю, что в ансис есть команда mdamp .

В итоге в ветке проекта HARMONIC RESPONSE я добавил ветку СOMMANDS APDL и далее вcтавил следующее:

 

MDAMP, 1, 0.035, 0.035, 0.035, 0.035, 0.035, 0.035
MDAMP, 2, 0.035, 0.035, 0.035, 0.035, 0.035, 0.035
MDAMP, 3, 0.035, 0.035, 0.005, 0.035, 0.035, 0.035
MDAMP, 4, 0.035, 0.035, 0.035, 0.035, 0.035, 0.035
MDAMP, 5, 0.030, 0.030, 0.030, 0.030, 0.030, 0.030
MDAMP, 6, 0.030, 0.030, 0.030, 0.030, 0.030, 0.010
MDAMP, 7, 0.010, 0.010, 0.010, 0.005, 0.008, 0.008
MDAMP, 8, 0.008, 0.004, 0.004, 0.008, 0.008, 0.008
MDAMP, 9, 0.008, 0.004, 0.004, 0.008, 0.008, 0.008
MDAMP, 10, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 11, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 12, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 13, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 14, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 15, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 16, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 17, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 18, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046

 

 - это мои экспериментальные данные 

 

провожу расчет в диапазоне от 0 до 2000 Гц (суперпозициями).

 

Так вот , хотел сказать , что на выходе я смотрю Frequency response в некоторых точках (деталях) модели . Модель представляет собой сборку из нескольких деталей . 

Получаю некие графики.

Далее я комментирую приведенные строки из ветки СOMMANDS APDL (например  комментирую их). Провожу снова гармонический расчет. В итоге графики не меняются в тех же точках конструкции.

Поэтому у меня склад. впечатление что MDAMP не работает в ворбенче , либо я что - то не правильно делаю.

 

Я извиняюсь , быть может такая тема ранее и обсуждалась...

 

 

 

[soklakov 1 475]

7 минут назад, fancysim сказал:

провожу расчет в диапазоне от 0 до 2000 Гц (суперпозициями).

Выполните полным методом. И повторите свои изыскания. Расскажите, что получится.

 

[fancysim 2]

10 минут назад, soklakov сказал:

Выполните полным методом. И повторите свои изыскания. Расскажите, что получится.

 

Спасибо , попробую. Но насколько я знаю (может и ошибаюсь) команда MDAMP работает только для метода суперпозиции. Так написано в хелпе. (Во всяком случае про фул метод в контексте mdampа инфы нет).

Тогда (если это так), то тогда я посчитаю фул методом без MDAMP  и cо СВОИМИ вставками MDAMP методом суперпозиции  и сравню результаты.

Изменено 8 июня 2017 пользователем fancysim

[soklakov 1 475]

10 минут назад, fancysim сказал:

Спасибо , попробую. Но насколько я знаю (может и ошибаюсь) команда MDAMP работает только для метода суперпозиции.

Вот она где собака-то зарыта... Вместо того, чтобы попробовать, мы думаем, что знаем и не пробуем. А если бы попробовали, то может и вопроса не было бы.

 

Я не уверен, что Ваш MDAMP сработает, не проверял синтаксис и даже не спрашиваю, что это за таблицы чисел. Но методом суперпозиции корректно учесть демпфирование - это большая задача. Особенно с учетом того, что Вы будете складывать формы, полученные из модального анализа, в котором Вы, почти наверняка, пренебрегли демпфированием. Начинать стоит с полного метода. А там разберетесь. Если будете делать, а не думать. ;)

 

[fancysim 2]

3 минуты назад, soklakov сказал:

Я не уверен, что Ваш MDAMP сработает, не проверял синтаксис и даже не спрашиваю, что это за таблицы чисел. Но методом суперпозиции корректно учесть демпфирование - это большая задача. Особенно с учетом того, что Вы будете складывать формы, полученные из модального анализа, в котором Вы, почти наверняка, пренебрегли демпфированием. Начинать стоит с полного метода. А там разберетесь. Если будете делать, а не думать. ;)

 

Ок , посмотрю, но

  синтаксис MDAMP  ТАКОВ: 

 

MDAMP, STLOC, V1, V2, V3, V4, V5, V6
Defines the damping ratios as a function of mode.

SOLUTION: Dynamic Options

MP ME ST PR PRN <> <> <> PP EME MFS

 

 

STLOC

Starting location in table for entering data. For example, if STLOC = 1, data input in the V1 field applies to the first constant in the table. If STLOC = 7, data input in the V1 field applies to the seventh constant in the table, etc. Defaults to the last location filled + 1.

V1, V2, V3, . . . , V6

Data assigned to six locations starting with STLOC. If a value is already in this location, it will be redefined. Blank values for V2 to V6 leave the corresponding previous value unchanged.

Notes

Defines the damping ratios as a function of mode. Table position corresponds to mode number. These ratios are added to the DMPRAT value, if defined. Use STAT command to list current values. Applies to the mode-superposition harmonic (ANTYPE,HARMIC), the mode-superposition linear transient dynamic (ANTYPE,TRANS), and the spectrum (ANTYPE,SPECTR) analyses. Repeat MDAMP command for additional constants (10000 maximum).

MDAMP can also be defined in a substructure analysis using component mode synthesis with fixed-interface method (ANTYPE,SUBSTR with CMSOPT,FIX and SEOPT,,,3). The damping ratios are added to the diagonal of the reduced damping matrix as explained inComponent Mode Synthesis (CMS).

This command is also valid in PREP7.

 

А значит в каждой строке я прописываю 6 модов  , а v1 ... v6 в каждой строке  - это значение коэффициента демпфирования для каждой моды.

у меня 18 строк x6 =108 мод . Такое ~ количество мод у меня в дипазоне от 0... до 2000.

 

Как-то так.

 

  

[soklakov 1 475]

4 минуты назад, fancysim сказал:

но

  синтаксис MDAMP  ТАКОВ: 

Я не смотрел, потому что мне лень) а не потому что у меня нет справки.

9 минут назад, fancysim сказал:

Как-то так.

 

Вообще, похожу на правду. Но ведь не работает)) Посмотрите логи. (solution information - solution output). Может быть там есть какой-то фидбэк, по поводу использованных команд. Например, сообщения вида "это команда из преп7 и будет проигнорирована".

1 час назад, fancysim сказал:

а именно собственные частоты и значения добротностей (коэффициентов демпфирования на этих частотах)

С недавних версий поддерживается ввод демпфирования в табличном виде, в зависимости от частоты. Без команд. Может стоит использовать?

[hr4d 87]

@fancysim А если не секрет, каким образом находили демпфирование по частотам?

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

 

[fancysim 2]

17 часов назад, hr4d сказал:

@fancysim А если не секрет, каким образом находили демпфирование по частотам?

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

 

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

22 часа назад, soklakov сказал:

С недавних версий поддерживается ввод демпфирования в табличном виде, в зависимости от частоты. Без команд. Может стоит использовать?

А если не сложно подскажите  , где это в ворбенче есть ? В каком разделе модального или гармонического анализа ? (пользуюсь последними версиями воркбенча).

[fancysim 2]

В 08.06.2017 в 20:32, hr4d сказал:

@fancysim А если не секрет, каким образом находили демпфирование по частотам?

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

 

@hr4d , а где об этом почитать ? Я бы сравнил с экспериментом. Кстати я делал простое исследование : простукивание подвешенной стальной пластины (сталь 20, экспериментальный модальный анализ), получал собственные частоты . Проводил для этой же пластины модальный анализ в воркбенче.  Точность (относительная погрешность) моделирования по основным собственным частотам в диапазоне от 0 до 2000 Гц   лежит в пределах от 0.1 до 3 %. 

Изменено 13 июня 2017 пользователем fancysim

[soklakov 1 475]

В 09.06.2017 в 14:23, fancysim сказал:

где это в ворбенче есть ?

Кажется, я что-то напутал) А может это была какая-нибудь ACTшка. В любом случае сейчас ничего такого найти не могу, ложная тревога.

В 13.06.2017 в 09:36, fancysim сказал:

 Точность (относительная погрешность) моделирования по основным собственным частотам в диапазоне от 0 до 2000 Гц   лежит в пределах от 0.1 до 3 %. 

Это замечательно. А точность определения добротности вы оценивали? Интересно, какие там цифры, и что именно делали, что с чем сравнивали.

 

[AlexKaz 1 821]

В 08.06.2017 в 22:32, hr4d сказал:

@fancysim А если не секрет, каким образом находили демпфирование по частотам?

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

 

а также

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

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

- способ акустическо/радиотехнический - по "колоколу" (добротности) - для собственной на АЧХ берётся шапка 3 дБ х 2 sigma; более подробно хотелось бы узнать о нём самому; скорее всего софт @fancysim'а использует этот метод; порой достаточно миллиметровой линейки и графика АЧХ;

- теоретический - по предложению товарища Рэлея демпфирование линейно растёт с частотой;

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

 

Изменено 15 июня 2017 пользователем AlexKaz

[hr4d 87]

@AlexKaz Я примерно о таком способе

123.docx

3234.djvu

[AlexKaz 1 821]

@hr4d , построил отклики во времени по графикам с помощью обратного многомерного Фурье-преобразования . Если смотреть на огибающую одной из шести функций, у меня получилось демпфирование в районе 23-25 1/c, в скане же - диапазон 10-18 1/c. Методику в скане вижу впервые, хотя "слышать слышал".
 

Скрытый текст

 

 

 

Изменено 18 июня 2017 пользователем AlexKaz

[fancysim 2]

@hr4d  Скажите пожалуйста , а откуда у Вас файл 3234.djvu ? Вы знакомы с С.И.Досько ? Я с ним знаком. И утверждаю - этот файл 3234.djvu я лично сканировал в ГПНТБ в этом году. Я сравнил размеры файла , что Вы приложили и которого сканировал - они одинаково размера!!! Правда файл на моем компе с именем : "ХОМЯКОВ_ДОСЬКО_Параметрическая_идентификация_динамических_систем_станков.djvu"

 

Изменено 19 июня 2017 пользователем fancysim

[hr4d 87]

Познакомился недавно. Вот собираемся вместе попробовать одну вещь.

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

 

[fancysim 2]

Ошибка у меня была в вставке MDAMP, у меня было ранее :

 

MDAMP, 1, 0.035, 0.035, 0.035, 0.035, 0.035, 0.035
MDAMP, 2, 0.035, 0.035, 0.035, 0.035, 0.035, 0.035
MDAMP, 3, 0.035, 0.035, 0.005, 0.035, 0.035, 0.035
MDAMP, 4, 0.035, 0.035, 0.035, 0.035, 0.035, 0.035
MDAMP, 5, 0.030, 0.030, 0.030, 0.030, 0.030, 0.030
MDAMP, 6, 0.030, 0.030, 0.030, 0.030, 0.030, 0.010
MDAMP, 7, 0.010, 0.010, 0.010, 0.005, 0.008, 0.008
MDAMP, 8, 0.008, 0.004, 0.004, 0.008, 0.008, 0.008
MDAMP, 9, 0.008, 0.004, 0.004, 0.008, 0.008, 0.008
MDAMP, 10, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 11, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 12, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 13, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 14, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 15, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 16, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 17, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 18, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046

 

 

 Но оказывается , первый параметр команды MDAMP – это стартовый номер, номер моды (собственной формы) !!!. Т.е. если, написано MDAMP, 55, 0.1, 0.2, 0.3, то для 55-й собственной формы будет задано демпфирование 0.1, для 56-й будет 0.2, а для 57-й будет 0.3. Если после этого выполнить MDAMP, 57, 0.4, 0.5, то теперь для 57-й будет 0.4, для 58-й 0.5.

 Соответственно, мне необходимо поправить команды, чтобы номера шли через 6, т.е. 1,7,13,… Что я и сделал.

 

Поэтому, получилось примерно так  :

 

MDAMP, 1, 0.035, 0.035, 0.035, 0.035, 0.035, 0.035
MDAMP, 7, 0.035, 0.035, 0.035, 0.035, 0.035, 0.035
MDAMP, 13, 0.035, 0.035, 0.005, 0.035, 0.035, 0.035
MDAMP, 19, 0.035, 0.035, 0.035, 0.035, 0.035, 0.035
MDAMP, 25, 0.0030, 0.0030, 0.0030, 0.0030, 0.0030, 0.0030
MDAMP, 31, 0.0036, 0.0036, 0.0036, 0.0036, 0.0036, 0.0036
MDAMP, 37, 0.010, 0.010, 0.010, 0.005, 0.008, 0.008
MDAMP, 43, 0.008, 0.004, 0.004, 0.008, 0.008, 0.0008
MDAMP, 49, 0.0008, 0.004, 0.004, 0.008, 0.008, 0.008
MDAMP, 55, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046, 0.0046
MDAMP, 61, 0.0046, 0.46, 0.46, 0.46, 0.46, 0.46
MDAMP, 67, 0.46, 0.046, 0.046, 0.046, 0.46, 0.46
MDAMP, 73, 0.046, 0.46, 0.46, 0.46, 0.46, 0.046
MDAMP, 79, 0.46, 0.46, 0.46, 0.46, 0.046, 0.046
MDAMP, 85, 0.46, 0.46, 0.46, 0.46, 0.46, 0.46
MDAMP, 91, 0.46, 0.46, 0.46, 0.46, 0.46, 0.46
MDAMP, 97, 0.46, 0.46, 0.46, 0.46, 0.46, 0.46
MDAMP, 103, 0.46, 0.46, 0.46, 0.46, 0.46, 0.46

 

Вот такая вставка работает и меняет выходные данные НARMONIC RESPONSE . Т.е. то что мне надо работает. Поэтому согласно данным экспериментального модального анализа я смогу откорректировать за счет модального демпфирования выходные амплитуды на каждой рез. частоте или на участках частот. Соответственно я прописываю модальное демпфирование для 109 мод (собственных форм).

  

Изменено 20 июня 2017 пользователем fancysim