Mass scaling при расчете SPH

[rasta89 113]

Всем привет.

Опытные дайняне, подскажите пожалуйста. Нужен совет.

При расчете, например резания, когда обьект резания представлен частицами SPH и необходимо оценить состояние режущего инструмента, приходится увеличивать время расчета, например до 1 сек. Это связано с тем, что при расчете например 0,0001 сек инструмент испытывает сильное напряжение из-за большой скорости движения. Но если при расчете на 0,0001 сек можно было обойтись без массового масштабирования, то на 1 сек уже приходится его применять, иначе время расчета тысячи часов. Итак можно/нужно использовать масштабирование. Но тут опять проблема, чтобы значительно снизить время расчета приходится брать число-коэффициент поменьше, что приводит к увеличению массы на миллионы, после чего частицы реагируют уже неадекватно, просто рассыпаются в стороны. Если следовать правилам, чтобы увеличение массы не превышало 5%, тогда и время расчета не значительно снижается. Замкнутый круг. Как быть? Как сократить время расчета? Можно еще немного уменьшить тайм степ скейл, но это не сильно помогает.

[AlexKaz 1 822]

Для начала в *define_curve надо поменять время и амплитуду местами.

В *control_timestep можно задать dt2ms.

И не забывайте про hourglass и damping.

 

 

[rasta89 113]

2 часа назад, AlexKaz сказал:

Сначала выложите небольшой рабочий проект с SPH. Увидел файл. Дело может быть и в других настройках.

Там размер не изменён и контрол_датабаза дурацкая

Вложил нормальный файл

wdfiles.ru/2NE7e

2 часа назад, AlexKaz сказал:

Для начала в *define_curve надо поменять время и амплитуду местами.

В *control_timestep можно задать dt2ms.

И не забывайте про hourglass и damping.

 

 

dt2ms это же и есть масс скейл

А зачем хоруглас для бетона?

Там модель Джонсона Холмквиста с ЕОСом, этого вроде как достаточно.

А какое время расчета? Контрол_терминейшн

2 часа назад, AlexKaz сказал:

*define_curve

Косяк, но я бы заметил потом если бы перемещение пошло

[AlexKaz 1 822]

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

А зачем хоруглас для бетона?

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

26 минут назад, rasta89 сказал:

Там модель Джонсона Холмквиста с ЕОСом, этого вроде как достаточно.

Было задано слишком много критериев разрушения. Решатель останавливался при первом же контакте тел из-за разрушений.

27 минут назад, rasta89 сказал:

А какое время расчета? Контрол_терминейшн

Прежнее. На скрине просто форма тел через ~0.5-1 минуту работы решателя.

[rasta89 113]

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

В качестве дополнительного демпфера

А исходя из чего вы выбрали этот параметр? Стоковый? Или для горных пород свой - ни разу не встречал, почитаю об этом ещё

 

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

Было задано слишком много критериев разрушения. Решатель останавливался при первом же контакте тел из-за разрушений.

И я на простом песчанике пробовал, который в базе был, на 0.0001 и с холмквистом считало же

 

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

Прежнее. На скрине просто форма тел через ~0.5-1 минуту работы решателя

Прежнее это какое? 1мс или 1 с? Имеете в виду что в файле было записано. Попробуйте на 1 сек, в этом вся соль.

 

По итогу у нас напряжения ГПа, размеры в мм, время в мс, плотность СИ. Да ещё и массскейл. Что скажете по результатам они выглядят адекватно?

Куча вопросов сразу скопилось)

[AlexKaz 1 822]

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

Куча вопросов сразу скопилось)

Вы попробуйте сами, а там посмотрим.

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

[rasta89 113]

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

Вы попробуйте сами, а там посмотрим.

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

Так, не понял. Там же не скорость задна а перемещение.

Так а как ее уменьшить-то. Например инструмент должен пройти за 1 сек 1 мм.

Что-то я запутался в край

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

Вы попробуйте сами, а там посмотрим.

Так я поэтому тут пишу, что не выходит

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

46 минут назад, rasta89 сказал:

 

Вы попробуйте сами, а там посмотрим.

Ещё хуже стало) тело не режется, ни на какое время. Я ранее кэписал что массскейл не приводит к добру

 

Блин, может внесёте правки в файл а я посмотрю? Ну и вопросы потом задам?

Шансов мало, но я спросил)

[AlexKaz 1 822]

43 минуты назад, rasta89 сказал:

Например инструмент должен пройти за 1 сек 1 мм.

Посмотрите на размеры деталей.

Ну, скорость или перемещение. Просто уменьшил на три порядка.

Изменённый файл.

test_v1_2.k.bz2

[rasta89 113]

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

Посмотрите на размеры деталей.

Ну, скорость или перемещение. Просто уменьшил на три порядка.

Изменённый файл.

test_v1_2.k.bz2 441 \u041a\u0431 · 0 скачиваний

Тогда получается если мы плотность подняли до СИ, то перемещение стало в метрах? Я логику хочу понять

[AlexKaz 1 822]

Вы не туда смотрите. Сначала задают длину-массу-время. Из них формируют всё остальное, включая плотность. Но затем плотностью можно варьировать по своему усмотрению (поглядывая на влияние динамики на результаты), остальными величинами варьировать нельзя.

[rasta89 113]

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

Посмотрите на размеры деталей.

Ну, скорость или перемещение. Просто уменьшил на три порядка

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

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

Вы не туда смотрите. Сначала задают длину-массу-время. Из них формируют всё остальное, включая плотность. Но затем плотностью можно варьировать по своему усмотрению (поглядывая на влияние динамики на результаты), остальными величинами варьировать нельзя.

Все понятно, спасибо.

 

А вот это не понятно, поменял плотности, ригид тело сожрало деформируемую насадку

[rasta89 113]

Все проде предельно понятно, но ничего не понятно

Файл

wdfiles.ru/27GO0

Исравлен размер детелаей.

[AlexKaz 1 822]

33 минуты назад, rasta89 сказал:

Исравлен размер детелаей.

Если изменена размерность у длины/времени/массы, то как минимум менять надо и константы материалов под новые размерности. Таблички выше есть

18 часов назад, AlexKaz сказал:

Изменённый файл.

test_v1_2.k.bz2 441 \u041a\u0431 · 1 скачивание

Используйте этот файл, здесь ужё всё настроено, только массскейлинг чрезмерен.

Изменено 31 августа 2023 пользователем AlexKaz

[rasta89 113]

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

Если изменена размерность у длины/времени/массы, то как минимум менять надо и константы материалов под новые размерности. Таблички выше есть

Используйте этот файл, здесь ужё всё настроено, только массскейлинг чрезмерен.

Так я же писал в самом начале об этом. Может я не так понял тогда.

С самого начала я так и делал, менял размерности материала плотность и контрастны, пытаясь добиться того чтобы я мог решать задачу не только на 0.0001 с но и на 1 секунде (время расчета). Я сказал, что пытаюсь сократить время расчета массскейлом но он портит жизнь частицам спш. И тогда появились предложение выключить массскейл и корректировать плотность. Но, смысл в смене полной системы юнитов (плотность, констатны, размеры моделей и т д), не даёт результата по вопросу о котором я говорю с самого начала, при 0.0001 с решается, при 1 сек время решения заоблачное. Мы опять получается к этому вернулись, что нужно включать массскейл но на разных юнитах (СЕИ)?

[fzbm 44]

Я бы для SPH вообще не применял массовое масштабирование

[rasta89 113]

23 минуты назад, fzbm сказал:

Я бы для SPH вообще не применял массовое масштабирование

Расскажите свое видение? Как быть? Если нужно увеличить время, скажем резания, но чтобы время расчета не было при этом 10000 часов?

[AlexKaz 1 822]

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

плотность

Здрасте. Без массскейлинга с исходной плотностью расчётное время на компе 3 часа 40 минут, без масс-га с плотностью выше на порядок - 1 час с чем-то. Вы что-то не то делаете.

[rasta89 113]

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

Здрасте. Без массскейлинга с исходной плотностью расчётное время на компе 3 часа 40 минут, без масс-га с плотностью выше на порядок - 1 час с чем-то. Вы что-то не то делаете.

Вы точно ничего не путаете?

Вот файл с исходной плотностью.

Время 1 мм, без масскейла, все юнииты в согласованной системе, время расчета 6750 часов

wdfiles.ru/LLUc

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

Здрасте. Без массскейлинга с исходной плотностью расчётное время на компе 3 часа 40 минут, без масс-га с плотностью выше на порядок - 1 час с чем-то. Вы что-то не то делаете.

В файле который вы прикладввали ранее, в первый раз, плотность была по си, а констатны в ГПа.

 

---

 

 

55 минут назад, fzbm сказал:

Я бы для SPH вообще не применял массовое масштабирование

 

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

Здрасте. Без массскейлинга с исходной плотностью расчётное время на компе 3 часа 40 минут, без масс-га с плотностью выше на порядок - 1 час с чем-то. Вы что-то не то делаете.

Коллеги, извиняюсь за наглость.

Может быть сделаем так.

Сделайте пожалуйста два файла из исходного, один с решением при 0.0001 сек, второй при 1 сек. И все стразу станет ясно.

Изменено 31 августа 2023 пользователем rasta89

[AlexKaz 1 822]

24 минуты назад, rasta89 сказал:

точно ничего не путаете?

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

 

Ещё раз - берёте МОЙ файл в архиве bz2 и меняете в картах *MAT плотности у упругого и пластичного материалов, отключаете масскейлинг.

Изменено 31 августа 2023 пользователем AlexKaz

[rasta89 113]

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

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

 

Ещё раз - берёте МОЙ файл в архиве bz2 и меняете в картах *MAT плотности у упругого и пластичного материалов.

У вас там плотность с системе си, а константы в ГПа. Вы же говорили с плотностью нужно и константы остальные менять.

[AlexKaz 1 822]

Вы не выспались? :) Вы бы передохнули, потом ещё раз перечитали топик, потому как у Вас опять выводы неверные.