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

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


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

1 минуту назад, NNK сказал:

Я ответил на это Ваше сообщение выше

всегда допускайте, что Вас троллят) это ж форум. этот форум)

(правда, в этот раз это не так, судя по всему, но с такими допущениями жить проще)

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


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

всегда допускайте, что Вас троллят) это ж форум. этот форум)

(правда, в этот раз это не так, судя по всему, но с такими допущениями жить проще)

Да мне не жалко )

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

влияет РАЗМЕР КЭ - увеличиваете размер => получаете "лучшее" решение

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

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

Вот что получилось в НЕЯВНОМ решателе:

мелкая сетка

fb5724e8e1fd871f6c89d40d49768055.jpg

 

крупная сетка

1e995f848996408b1540836fdbcf66a1.jpg

 

Тут-то шаг интегрирования безусловен...

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

Тут-то шаг интегрирования безусловен...

Вы выполнили фильтрацию сеткой. 2% демпфирования впихните и у нас будут одинаковые картинки реакции)

А если не будут - будет интересно и продолжим.

@NNK Хотя я вот смотрю - они кажется уже одинаковые.

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

Вы выполнили фильтрацию сеткой. 2% демпфирования впихните и у нас будут одинаковые картинки реакции)

А если не будут - будет интересно и продолжим.

@NNK Хотя я вот смотрю - они кажется уже одинаковые.

Фильтрация сеткой означает неучет высших форм?

Вот при 2% демпфировании:

4497646ea8a5d31fb95093dad847a347.jpg

Видно, что к концу второго периода амплитуда стала меньше, но высокочастотные пики никуда не делись

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

Фильтрация сеткой означает неучет высших форм?

угу.

53 минуты назад, NNK сказал:

Вот при 2% демпфировании:

А вот это уже странно. Есть ли у Вас уверенность, что демпфирование действительно задано и именно 2%? К сожалению, не смогу проверить Ваш проект.

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

Не все варианты демпфирования испробованы =), в ANSYS есть возможность задать минимум два: пропорционально частоте и по Рэлею. В LS-Dyna есть понятие релаксации энергии. Вряд ли возможности abaqus скромнее. 

Victor_M выше написал про длину балочного элемента. На графиках сетка равномерная. А если в заделке будет более длинный элемент по сравнению с остальными - колебания напряжений в заделке сгладятся?

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

Кстати, хороший вопрос про демпфирование. Можно задать по Релею, коэффициентами альфа и бета. Каждый из них пропорционален собственной частоте и доле от критического. Если бы мы рассматривали модальную динамику, то можно было бы задать на каждую учитываемую форму свой коэффициент. Что делать в случае прямого интегрирования, мне не совсем понятно - там тоже задаются альфа и бета, но как их определить. Выше я использовал демпфирование, пропорциональное массе, альфа вычислена по первой частоте. Теперь это как-то сомнительно.

Я перешел на линейное модальное решение - гулять, так гулять. Тут хотя бы можно явно выбрать количество учитываемых форм. Так вот при учете только первой формы реакция составляет 3-3.1 кН (при разном разбиении,  от 4 до 200 элементов по длине). 2 кН получается только при 1 элементе. Истина где-то рядом...

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

Victor_M выше написал про длину балочного элемента. На графиках сетка равномерная. А если в заделке будет более длинный элемент по сравнению с остальными - колебания напряжений в заделке сгладятся?

расчет времени интегрирования ведется на основании НАИМЕНЬШЕГО КЭ в модели

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

Выше я использовал демпфирование, пропорциональное массе, альфа вычислена по первой частоте.

используйте бету

:smile:

damp.png

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

@NNK

подставьте beta= 1.5158e-004

34 минуты назад, NNK сказал:

Истина где-то рядом...

тема уже длинная, цифр много, графики, таблицы какие-то. попробуйте переначать) что сейчас Вас больше всего беспокоит?

ну кроме того, что ничего непонятно)

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

Попробуем переначать ) Тем более от первого сообщения я уже несколько отошел.

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

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

Бету подставил, задемпфировалось сильнее:

два решения - при альфа и при бета

e157ed285a3c254e1d3103e2da5bd530.jpg

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

Предлагается не переходить на силу, а продолжать работать с исходной задачей (с массой):

 

Модель в ANSYS Static Sructural

 1. Стальная балка (beam) сечением 50 х 50 х 1000 мм с уменьшенной плотностью 78.50  кг/мм3

 2. Груз (rigid) - шар d10 мм массой 100 кг

 3. Сила тяжести 9806.65 мм/сек2

 4. Малые перемещения Z стат = 3.15 мм

 5. Реакция Z стат = 982.65 Н

 6. По формуле (30.10) К дин = 3.94

 

Модель в ANSYS Transient:

 1. Начальная скорость груза при падении с высоты 12 мм - 485 мм/с

 2. Сила тяжести 9806.65 мм/сек2

 3. Малые перемещения

 4. Damping Ratio 0.02 на 9 Гц. Картинка красивая - высокочастотные колебания затухают до 0.005 сек

 5. Деформация Z дин = 12.02 мм

 6. Реакция Z дин =  3752.00 Н

 

Итак, К дин = 3752/983 = 3.82

 

 

 

 

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

"Правильное" значение реакции

Таки оно неправильное.

В 19.01.2018 в 16:31, soklakov сказал:

что сейчас Вас больше всего беспокоит?

 

Ссылка на сообщение
Поделиться на других сайтах
В ‎17‎.‎01‎.‎2018 в 15:31, AlexKaz сказал:

решение даёт погрешность менее 5%, вообще то

 

В ‎22‎.‎01‎.‎2018 в 11:39, imt сказал:

По формуле (30.10) К дин = 3.94

 

В ‎22‎.‎01‎.‎2018 в 11:39, imt сказал:

ANSYS Transient К дин = 3752/983 = 3.82

А может, кто выложит проект ANSYS Dynamic Explicit или Autodyn  с правильным ответом?...

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

с правильным ответом?

правильный ответ - это который?

Коэффициент динамичности - концепт, который позволяет приближенно отыскивать решение динамической задачи, не решая ее.

Своего рода калибровочный коэффициент. Явные решения тут при чем?

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

Явные решения тут при чем?

Захотелось на чужом горбу в рай въехать :)...

С явным решением у меня даже модель не работает, не говоря уже о правильном ответе.

Опять же настройки решения, демпфирования всякие ...

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • Алекале
      @esergey Подскажите а для того что бы создавать "интерактивное электронное технологическое руководство" с 22 по 26 минуту в вебинаре, для этого нужно создавать техпроцесс в системе тимцентр с электронной моделью ДСЕ, или же можно как то поверхностно этот путь пройти?
    • Nalad4eG
      Как этот "свист" сказывается на тех характеристиках? Если не существенно, то мой вас совет не трогайте будет только хуже.
    • lem_on
      Поздравляю. Но замена подшипников на таком шпинделе тот ещё квест. Нужно название шпинделя, там наклейка на борту и по нему искать чертеж. Подшипники скорее всего керамические. Сначала найти чертеж, а потом разбирать. А если нет опыта то лучше отдать тем кто умеет. Балансировку обязательно и стразу пакет тарельчатых пружин со всеми уплотнениями перетрясти. 
    • Ветерок
      Зачем это нажимать? Лучше нажать проецировать кривую. Как выяснилось :) такая функция есть. @mannul даже файл выложил - бери и смотри/повторяй.
    • Stanislav
      Относительно POP and PUP в фале постпроцессора вот такое нашел для POP прослеживается возростающая закономерность, а вот для  PUP не ясно почему прыгает значение. # POP and PUP functions for initialization. H parameters are used to  # # get the correct POP and PUP values. If the thickness of the sheet   # # is between the minimum and maximum thickness, the postprocessor     # # will use POP and PUP parameters to write them in the PS file.       # # Example: H000 0.0 1.0 30000 25000                                   # #                                                                     # #     Minimum      Maximum       POP        PUP                       # #    Thickness    Thickness     Value      Value                      # #---------------------------------------------------------------------# H001    0.0          0.5        14000      9000 H002    0.51         1.0        14500      9000 H003    1.01         1.5        15000      9000 H004    1.51         2.0        15500      9000 H005    2.01         2.5        16000      8500 H006    2.51         3.0        16500      8000 H007    3.01         3.5        17000      8000 H008    3.51         4.0        17500      8000 H009    4.01         4.5        18000      9000 H010    4.51         5.0        18500      8000 H011    5.01         5.5        19000      9000 H012    5.51         6.0        19500      8000
    • Krusnik
      Да там же ничего сложного.   Делаете вашу спираль.   Потом делаете эскиз. В эскизе нажимаете Преобразование объектов и тыкаете на созданную спираль.   Нажимаете поверхность вытянуть и выбираете только что созданный эскиз. В опциях вытягивание ставите "До поверхности". И вытягиваете спиральку до вашего купола.   И последний шаг - вырез по траектории. Думаю с этим вы без моих советов справитесь.   Лишнее скрываете и вуаля.  
    • IgP
      Это, конечно, да, не хорошо ... Но и конструктор не бог чтобы всё знать ... Для этого есть и д.б. этап "технологическая проработка", со всеми вытекающими.
    • arsenev
      Спасибо за ответ. Нет, в видео он говорит, что это такая стратегия. Пробовал использовать шаблон и обрабатывать его со смещением, получается не то (именно, когда нужно обработать цилиндрическое отверстие). Лучше отводится стружка и не мешает резанию. (см. видео, где-то с 18 минуты).
    • scrimline
      Если ещё нужна информация, fanuc oi-tf станок akira-seiki, М14 на опускание руки, для привязки и М15 для поднятия
    • Богоманшин Игорь
      Эта проблема от незнания некоторыми конструкторами технологических особенностей порошковой окраски. Максимум на чертеже научатся указывать РАЛ и толщину слоя. А необходимо предусмотреть места подвеса, можно туда шпильки или гайки запрессовать в нужных местах.
×
×
  • Создать...