Перейти к публикации

Динамика - падение груза на консоль, явное решение против теории


Рекомендованные сообщения

11 минуту назад, imt сказал:

взять реализацию, где ANSYS дребезжит, что АЧХ покажет?

Амплитуды вырастут, частоты слегка сдвинутся. Какая частота и куда уйдёт - сказать заранее нельзя. Но первые несколько частот практически не изменятся.

Добавить демпфирование - амплитуды упадут, высокие частоты исчезнут. Чем выше демпфирование, тем меньше всплесков будет на АЧХ.

Изменено пользователем AlexKaz
Ссылка на сообщение
Поделиться на других сайтах


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

Какой можно сделать вывод

Это я как бы намекаю @NNK , что пора бы уже демпфирование добавить ...

 

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

только их нет смысла сравнивать

... Вы хотите сказать, что на АЧХ в низкочастотной области будет только один пик на 9 Гц?

Глядя на реализацию @soklakov , не смогу с Вами согласиться.

25 минут назад, AlexKaz сказал:

Но первые несколько частот практически не изменятся

Я как раз про первые первые частоты ...

На АЧХ от @NNK нет 42 Гц.

@soklakov говорит про 42 Гц - это что?

Собственная частота в поперечной плоскости? Поэтому ее нет у @NNK ?

Изменено пользователем imt
Ссылка на сообщение
Поделиться на других сайтах
18 минут назад, imt сказал:

Это я как бы намекаю @NNK , что пора бы уже демпфирование добавить ...

 

... Вы хотите сказать, что на АЧХ в низкочастотной области будет только один пик на 9 Гц?

Глядя на реализацию, не смогу с Вами согласиться.

Я как раз про первые первые частоты ...

На АЧХ от @NNK нет 42 Гц.

@soklakov говорит про 42 Гц - это что?

Собственная частота в поперечной плоскости? Поэтому ее нет у @NNK ?

Если совсем точно, то первая собственная частота в моем решении 41.826 Гц

Прогоню с демпфированием, посмотрим что получится

Ссылка на сообщение
Поделиться на других сайтах
6 минут назад, NNK сказал:

Если совсем точно, то первая собственная частота в моем решении 41.826 Гц

Точно. Ведь 9 Гц это с прилипшим грузом...

Ссылка на сообщение
Поделиться на других сайтах
6 минут назад, NNK сказал:

Если совсем точно, то первая собственная частота в моем решении 41.826 Гц

Точно. Ведь 9 Гц это с прилипшим грузом...

Ссылка на сообщение
Поделиться на других сайтах

в обед решил ту же задачку, но без прилипания массы к балке

 

получилось куда как "лучше".....

 

НО - как я уже писал - хорошо бы увидеть ВЫВОД той формулы - про к-т динамичности - да бы понять - в какой мере вообще правомерно рассматривать ТУ ФОРМУЛУ - как эталон?

 

 

 

просьба к автору этой темы: где это можно прочитать - подскажите пожалуйста?

 

1.png

Изменено пользователем Victor_M
Ссылка на сообщение
Поделиться на других сайтах
34 минуты назад, Victor_M сказал:

в обед решил ту же задачку, но без прилипания массы к балке

 

получилось куда как "лучше".....

 

НО - как я уже писал - хорошо бы увидеть ВЫВОД той формулы - про к-т динамичности - да бы понять - в какой мере вообще правомерно рассматривать ТУ ФОРМУЛУ - как эталон?

 

нужно понять - ЧТО и с ЧЕМ сравнивает автор?

 

просьба к автору этой темы: где это можно прочитать - подскажите пожалуйста?

 

1.png

 

1.png

Изменено пользователем Victor_M
Ссылка на сообщение
Поделиться на других сайтах

@Victor_M , http://rk5.msk.ru/books/Teor_Kolebaniy/Panovko.pdf - Пановко, Основы прикладной теории колебаний и удара, стр. 303-304.

Ссылка на сообщение
Поделиться на других сайтах
11 минуту назад, AlexKaz сказал:

@Victor_M , http://rk5.msk.ru/books/Teor_Kolebaniy/Panovko.pdf - Пановко, Основы прикладной теории колебаний и удара, стр. 303-304.

Да, например. Или Беляев - Сопротивление материалов, там, может, чуть более обширное описание.

Ссылка на сообщение
Поделиться на других сайтах
1 минуту назад, NNK сказал:

Да, например. Или Беляев - Сопротивление материалов, там, может, чуть более обширное описание.

пожалуйста, укажите страницы?

ну не искать же по всей книге.....

нужно точно понимать чем Вы пользуетесь....

14 минуты назад, AlexKaz сказал:

@Victor_M , http://rk5.msk.ru/books/Teor_Kolebaniy/Panovko.pdf - Пановко, Основы прикладной теории колебаний и удара, стр. 303-304.

здесь нет вывода именно тех формул.... - там общие слова - типа "бла-бла-бла.... как всем известно...."

нужен ВЫВОД - и те условия - та постановка задачи - для которой эти формулы выводились.....

Ссылка на сообщение
Поделиться на других сайтах

Ну ок. Давайте подключим элементарную математику. Заглянем в сопромат, чтобы не забыть что есть напряжения: %sigma=M/W=(E*J/W)*diff(y,t,2) - все это помнят.

Рассматриваем случай, перемещения описываются синусоидой y(t)=A*sin(w*t), где A - амплитуда перемещений, w - круговая частота колебаний. Заметим про себя, что амплитуда перемещений в нашем случае не зависит от частоты. А что насчёт амплитуды напряжений? => %sigma = (E*J/W)*(- A *w^2) *sin(w*t). Все теперь заметили, что в амплитуды напряжений входит круговая частота w? 

Дальше - больше. Рассматриваем случай двухмассовой системы  y(t)=A1*sin(w1*t) + A2*sin(w2*t), где w2>w1. Догадайтесь, что теперь произойдёт с графиком напряжений.

А теперь возьмите 10-20-милльён масс и продифференцируйте уравнение перемещений...

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

 

Изменено пользователем AlexKaz
Ссылка на сообщение
Поделиться на других сайтах

Добавил 2% демпфирование в задачу с мгновенным приложением силы:

8a2d742e9a7ef143d63580a8c0bbf0c2.jpg

Максимальные перемещения - "почти", напряжения - тоже "почти", а реакция - в 2 раза больше. Энергия деформации изогнутой консоли пропорциональна квадрату силы. Энергия в статической задаче составляет 1.5Дж, в динамической - 6, т. е. в 4 раза больше. То есть сила должна увеличиться в 2 раза. У меня же сила 4 кН вместо 2. Где-то собака зарыта...

6 минут назад, Victor_M сказал:

пожалуйста, укажите страницы?

ну не искать же по всей книге.....

нужно точно понимать чем Вы пользуетесь....

здесь нет вывода именно тех формул.... - там общие слова - типа "бла-бла-бла.... как всем известно...."

нужен ВЫВОД - и те условия - та постановка задачи - для которой эти формулы выводились.....

Стр. 515-517.

Изменено пользователем NNK
Ссылка на сообщение
Поделиться на других сайтах

в данном случае все "проще":

приводится формула и она сравнивается с численным решением...

многое мне не понятно, в частности - ЧТО является физической моделью для вывода ТОЙ формулы и на сколько та модель - похожа на ту - которую автор численно моделировал ?

Ссылка на сообщение
Поделиться на других сайтах
7 минут назад, Victor_M сказал:

в данном случае все "проще":

приводится формула и она сравнивается с численным решением...

многое мне не понятно, в частности - ЧТО является физической моделью для вывода ТОЙ формулы и на сколько та модель - похожа на ту - которую автор численно моделировал ?

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

Ссылка на сообщение
Поделиться на других сайтах
38 минут назад, AlexKaz сказал:

Ну ок. Давайте подключим элементарную математику. Заглянем в сопромат, чтобы не забыть что есть напряжения: %sigma=M/W=(E*J/W)*diff(y,t,2) - все это помнят.

Рассматриваем случай, перемещения описываются синусоидой y(t)=A*sin(w*t), где A - амплитуда перемещений, w - круговая частота колебаний. Заметим про себя, что амплитуда перемещений в нашем случае не зависит от частоты. А что насчёт амплитуды напряжений? => %sigma = (E*J/W)*(- A *w^2) *sin(w*t). Все теперь заметили, что в амплитуды напряжений входит круговая частота w? 

Дальше - больше. Рассматриваем случай двухмассовой системы  y(t)=A1*sin(w1*t) + A2*sin(w2*t), где w2>w1. Догадайтесь, что теперь произойдёт с графиком напряжений.

А теперь возьмите 10-20-милльён масс и продифференцируйте уравнение перемещений...

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

 

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

Учет демпфирования, как видно выше, не привел к гашению высоких частот.

Я решил ту же задачу с мгновенной силой, но всего с одним элементом по длине балки. Результат почти совпал с теорией:

f9725d26f631edbcabf0374443809e7a.jpg

Напряжения в 2 раза меньше, но это правильно - они вычисляются в середине элемента.

 

 

Ссылка на сообщение
Поделиться на других сайтах
44 минуты назад, NNK сказал:

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

Конечно. Формулам ВЕРЮ!

Если их применять, то применять точно так и никак иначе.

Какие сомнения после прочтения Пановко (стр. 303-304) (Спасибо @AlexKaz ):

1. Масса балки почти 20 кг. Можно, конечно, и пренебречь по сравнению со 100 кг. Я бы не стал.

2. У вас консольная балка, т.е. маятник с прилипшим грузом. Реакция у маятника как вычисляется? Я не знаю, поэтому сделал бы расчетную схему, как нарисовано  у Пановко.

3. Уменьшил бы размерность задачи, как убедительно продемонстрировал @AlexKaz

Ссылка на сообщение
Поделиться на других сайтах
7 часов назад, NNK сказал:

Я решил ту же задачу с мгновенной силой, но всего с одним элементом по длине балки. Результат почти совпал с теорией

Напряжения в 2 раза меньше, но это правильно - они вычисляются в середине элемента.

т.е. если взять два КЭ - то эта математическая модель даст уже "неверные" результаты?

т.е. они будут отличаться от "эталона" - от формулы?

 

получается более чем "странная" ситуация....

 

м.б. все же забыть про формулу и не рассматривать её как "эталон"?

Изменено пользователем Victor_M
Ссылка на сообщение
Поделиться на других сайтах

кажется я понял в чем причина "нефизичности" результатов численного моделирования:

 

в описании постановки задачи - для вывода формулы к-та динамичности есть "гипотеза":

контакт/столкновение происходит "быстро", но за КОНЕЧНОЕ время - это в данном случае играет ключевую роль

 

в данной мат.модели - время "контакта" (время изменения энергии) -  в точке прилипания массы к балке в точке касания - определяет шаг интегрирования по времени, который в свою очередь зависит от модуля упругости и от минимального РАЗМЕРА КЭ 

т.е. в данном случае - чем больше КЭ задается - тем самым меньше размер КЭ - тем более нефизичный результат получаем.

*** учитывая что сила = масса умноженная на ускорение, а ускорение = скорость деленная на время (на шаг интегрирования по времени) - поскольку заданное "действие" происходит мгновенно - за один шаг интегрирования.

 

косвенным подтверждением правильности моего предположения - является результат численного моделирования в моем сообщении ID: 48 , , где я задал возможность отскока массы от балки (убрал эффект прилипания), таким образом программа получили возможность "регулировать" время контакта, учитывая все параметры движения и "начальные" условия.

 

- в мат.модели с возможностью отскока - увеличение кол-ва КЭ (в разумных пределах) ведет к "повышению" точности решения - результат получается все более физичным.

 

 

 

РЕЗЮМЕ:

 

1. поскольку "эталон" - имеет неявную зависимость от времени контакта - эту задачу не стоило использовать в качестве "эталонной"

 

2. Математическая модель - сделана "неудачно" - ибо она с точки зрения физики - лишает алгоритм решения "выяснять" реальное время "контакта"

Изменено пользователем Victor_M
Ссылка на сообщение
Поделиться на других сайтах
13 часа назад, imt сказал:

консольная балка, т.е. маятник с прилипшим грузом

В данном случае не важно, т.к. при r=1 m=J v=w

Ссылка на сообщение
Поделиться на других сайтах

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

Ссылка на сообщение
Поделиться на других сайтах

Присоединяйтесь к обсуждению

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

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставить изображения напрямую. Загрузите или вставьте изображения по ссылке.

  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу.




  • Сообщения

    • gudstartup
      так он так и позиционируется по m19 pmc выдает g70.6 а чпу отвечает f45.7 но ориентацию и смещение в 4077 он отрабатывает нормально шпиндель встает ровно и смена происходит хорошо. вопрос почему после ввода команды управления скоростью он все еще продолжает контролировать число импульсов между нуль метками хотя в принципе уже должен отменить позиционный контроль и просто считать обороты по 0 метке как он это делает без М19? это все понятно но почему оно продолжает проверять это после завершения ориентации мне непонятно
    • Александр1979
      SP9047 SSPA:47 ILLEGAL SIGNAL OF POSITION CODER "The relationship between the A/B phase and 1-rotation signal is incorrect (Pulse interval mismatch)." "Неправильное значение счетчика импульсов сигнала на энкодере ALPHAi. На фазах A и B энкодера за один оборот шпинделя насчитывается 4096 импульсов обратной связи. Программное обеспечение по управлению шпинделем проверяет количество импульсов на фазах A и B, соответствующее энкодеру, при каждой генерации сигнала одного оборота. Данный аварийный сигнал срабатывает, если регистрируется число импульсов, нарушающее заданный диапазон."
    • gudstartup
      после ввода 789 в edit mod надо перейти в режим ручное упр. нажать mod а потом помощь и откроется wzw.hlp  с возможностью выполнения команд управления atc. если не появляется содержимое wzw.hlp то возможно что у вас другой код хотя 789 'это стандарт для dmg в этом случае надо смотреть ваш plc раздел у меня пароль прописан в файле wzw_rett.mod
    • Anat2015
      SERVO TURRET AMP ALARM - переводить не пробовали, учитывая, что AMP - это и есть сервопривод? Что тут не ясного, язык не знаком?
    • Viktor2004
      я бы попробовал позиционировать шпиндель сигналом G70.6 Ответный сигнал, шпиндель сориентирован F45.7
    • gudstartup
      Здравствуйте уважаемые участники форума попрошу вас поделится мнениями. Имеется фрезерный шпиндель с позиционным кодером (зубчатое колесо и считывающая головка). Ошибку 9047 выдает после того отрабатывает M19 c последующим включением оборотов Sxxx M3 причем обороты могут быть любыми от 100-макс. Если задавать обороты без предварительной ориентации то 9047 не возникает и шпиндель нормально вращается на любой скорости при этом обороты не гуляют точность поддержания +-1 об\мин. Кабель заменен заменить сам датчик нет возможности. Мой вопрос в следующем почему шпиндель после M19 и включении управления скоростью Sxxx M3 остается в режиме позиционирования а не сразу переходит врежим управления скоростью и даже явная отмена ориентации кодом M20 ничего не меняет ведь режим позиционного управления должен быть отменен или он остается и отменяется не сразу после подачи команды M3 и кода скорости? Станок достаточно старый 2009 года плату управления шпинделем меняли на соседний ничего не изменилось. Это может быть логика PMC или скорее всего это датчик?   PMC-SB.001
    • Killerchik
      А мы-то думали, что пропорциональный гидравлический.....
    • Grig75
    • EnKlL
      Я написал полностью ошибку, серво привод
    • gudstartup
      https://www.aggsoft.com/cnc-dnc/maho-cnc-432-1000.htm
×
×
  • Создать...