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

Нелинейный грунт EDP - проблема сходимости


Chingachguk

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

Вечер добрый!

Моделирую работу стальной сваи из круглой трубы с ребром в грунте. Собрал тестовую задачу - четверть.

Материал для стальных элементов нелинейный, билинейная диаграмма + кинематическое упрочнение.

Для грунта на первых этапах решил взять Extended Drucker-Prager с формулировками

TB,EDP,matid,,,LYFUN
....
TB,EDP,matid,,,LFPOT

Между трубой и грунтом контакт Frictional с коэффициентом 0,4. Собственно модель простая и сетка относительно ровная, НО задача не решается. Даже если задать жесткий контакт то вроде первые две три итерации решаются, но потом все равно процесс расходится. С линейным грунтом все нормально.

Выдает ошибку

An unknown error occurred during solution.  Check the Solver Output on the Solution Information object for possible causes.

В отчете предупреждение

EDP plasticity: plastic algorithm does not converge for element 47472 and material 102.

Кто нибудь сталкивался с похожим? может настройки решателя... или как-то EDP настроен не правильно? хотя на простых примерах типа осадки штампа решается нормально.

ссылка на *.wbpz проекта.

https://drive.google.com/open?id=1pwEE7wX0gRjO1E1EzjfRzRVi6LimbthU

Спасиб.

 

 

01.PNG

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


9 часов назад, Chingachguk сказал:

Кто нибудь сталкивался с похожим?

Да постоянно. Вообще не очень понятно, класс каких задач можно решать с использованием DP/EPD.

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

EDP plasticity: plastic algorithm does not converge for element 47472 and material 102.

Кто нибудь сталкивался с похожим? может настройки решателя...

Элементы надо брать квадратичными.

Настройки шага-подшага сильно влияют и прочие опции могут сильно влиять.

А по мне лучше перенести проект в APDL в классику - там веселее отлавливать глюки, настраивать решатели, да и их выбор больше.

Если совсем припрёт - переходите в LS-DYNA. Много головной боли исчезнет, т.к. схема Ньюмарка, имеющаяся в LS-DYNA, в стандартных солверах ANSYS не реализована.

 

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

Элементы надо брать квадратичными.

иногда напротив, стоит перейти на линейные в грунте, чтобы улучшить сходимость. в LS-dyna часто квадратичные используете?

 

11 час назад, Chingachguk сказал:

Для грунта на первых этапах решил взять Extended Drucker-Prager с формулировками

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

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

в LS-dyna часто квадратичные используете?

Каким образом связаны квадраты в LS-DYNA с квадратами в ANSYS? Насколько мне известно - никаким. Но всегда интересно узнать что-то новое.

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

Каким образом связаны квадраты в LS-DYNA с квадратами в ANSYS?

а что если посмотреть на вопрос вне контекста? :smile:

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

Солвер решает больше квадратов, о чём я и написал - go to "схема Ньюмарка википедия".

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

go to "схема Ньюмарка википедия"

то ли русской статьи нет, то ли запрос нужен другой.

интегрирование по времени по Ньюмарку в Ansys стоит по дефолту уже не знаю сколько лет. про запас есть вариант HHT.

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

если это вдруг секрет какой, то вовсе необязательно давать ссылки на схему Ньюмарка. секрет, так секрет.)

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

часто ли Вы использовали квадратичные элементы в LS-Dyna?

Ни разу. Я выбираю жизнь без наркотиков double float.

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

интегрирование по времени по Ньюмарку в Ansys стоит по дефолту уже не знаю сколько лет

Что ж тогда так хреново сходится? Что-то тут не так...

Статья https://en.wikipedia.org/wiki/Newmark-beta_method

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

Что ж тогда так хреново сходится?

схема неявная. чего уж.

с этого января можно попросить, чтобы когда плохо сходится, автоматически переключался на явную решалку, тогда сходится ловчее.

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

Если бы только "явная-неявная", сходится отлично - шаг требует менять. ANSYS до 16-й версии здесь тупит - после уменьшения шага и сходимости - шаг обратно не увеличивается. В итоге процесс решается ну оооочень долго. Как сейчас - не в курсе. A в дайне всё шустренько.

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

ANSYS до 16-й версии здесь тупит

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

а сегодня уже Ansys 2019 R1.

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

тогда сходится ловчее.

я вот тут неправильно выразился. явные алгоритмы в принципе не сходятся. поэтому у них и нет проблем со сходимостью.

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

A в дайне всё шустренько.

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

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

явные алгоритмы в принципе не сходятся. поэтому у них и нет проблем со сходимостью.

Эээ, вот почему нада хоть раз пописАть свой код - Ньюмарк де-факто явная схема. Хоть и обзывается иначе, хоть и быстрее явных в миллион раз. Понятия сходимости в этой схеме отсутствует. Каким макаром прогеры ансиса её переврали до аналога Ньютона-Рафсона - история до меня не дошла.

Предположу, что после шага по схеме Ньюмарка пытаются уравновесить суммы "силы + реакции" к нулю каким-нибудь методом, тем же Н-Р, а значит приходится в цикле вычислять Ньюмарка. Получается гибрид.

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

Эээ, вот почему нада хоть раз пописАть свой код - Ньюмарк де-факто явная схема.

Методе Ньюмарка известный неявный метод

Для этого нужно хоть раз написать свой код.)

 

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

Каким макаром прогеры ансиса её переврали до аналога Ньютона-Рафсона - история до меня не дошла.

ньютон рафсон - это не метод интегрирования по времени. все в кашу смешали.

 

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

Предположу, что после шага по схеме Ньюмарка пытаются уравновесить суммы "силы + реакции" к нулю каким-нибудь методом, тем же Н-Р, а значит приходится в цикле вычислять Ньюмарка.

а вот это похоже на правду

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

Ньюмарк де-факто явная схема.

в ваших терминах "явная" - это  там где время есть?

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

Элементы надо брать квадратичными.

Пробовал... не помогло, только время естественно увеличивается. В металлах, у меня это играло значение, особенно если считать закритику...

Цитата

 

А по мне лучше перенести проект в APDL в классику - там веселее отлавливать глюки, настраивать решатели, да и их выбор больше.

 

Собственно геометрию в этом проекте делаю в макросе, перешел в WB - потестировать мешер = хорошо сделали. Да и просто посмотреть как работает WB. А решатель я все равно вставками кода настраиваю. Но думаю придется вернуться к классике. Кстати чтоб два раза не вставать, в версии 19.2 в WB в свитке Component Systems не нашел Finite Element Modeler. Видел как через него подгружают проекты сделанные в классике на апдл. Его удалили? Кто-нибудь в курсе, как сейчас можно подгружать макросы с созданием геометрии?

LS-DYNA начал изучать, думаю в ней тоже попробовать, но позже.

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

попробуйте сначала кинематический с упрочнением(да побольше). потом без упрочнения. если с упрочнением работать будет, а без него - валиться, то задайте упрочнение в друкера-прагера, да побольше

Под упрочнением к EDP вы имеете ввиду это?

https://www.sharcnet.ca/Software/Ansys/16.2.3/en-us/help/ans_str/Hlp_G_STR8_3.html

/prep7
mp,ex,1,2.1e4			! Elastic Properties
mp,nuxy,1,0.1

ys=7.894657
sl=1000.0

tb,edp,1,,,LYFUN
tbdata,1,2.2526,ys

tb,edp,1,,,LFPOT
tbdata,1,0.566206

tb,plas,1,1,2,miso
tbpt,defi,0.0,7.894
tbpt,defi,1,1007.894

 

5 часов назад, soklakov сказал:

с этого января можно попросить, чтобы когда плохо сходится, автоматически переключался на явную решалку, тогда сходится ловчее.

Имеете ввиду версию 19R1 ?

PS: Поставил Mohr-Coulomb и решение сошлось, правда пришлось задирать угол дилатансии, но есть одно НО для глинистых грунтов он 0 или около...

 

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

в ваших терминах "явная" - это  там где время есть?

В моих - когда проверка сходимости не нужна, и не требуется искать корни уравнения итерациями.

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

Его удалили?

нет. спрятали.

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

Кто-нибудь в курсе, как сейчас можно подгружать макросы с созданием геометрии?

вместо него теперь external model

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

Имеете ввиду версию 19R1 ?

2019 R1

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

@soklakov

А есть ли возможность задать коэффициент трения хотяб линейно изменяющийся в процессе расчёта? В идеале конечно зависимостью от контактных напряжений....

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

@Chingachguk  может и есть. 

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

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

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • karlf
      Посмотрю в понедельник. Ну а какова вообще вероятность того, что дело в ПО, а не в железе?
    • vladimir.songin
      Так между прочим, кто-то на форуме фотографии очень похожего преса показывал. У них в поле энкодера значение энкодера.   Вся электросхема пожалуста.  GLADIATOR E3A 115-33 schema elettrico.pdf  
    • v-david
      Вы выбрали конфигурацию, но не активировали ее:       Set swCustPropMgr = CmpDoc.Extension.CustomPropertyManager(CompConfig) '------------------------------------------------------------------------------- lRetVal = CmpDoc.ShowConfiguration2(CompConfig) '-------------------------------------------------------------------------------     lRetVal = swCustPropMgr.Get4("Наименование", False, Val, ValOut)     If ValOut = Empty Then         Set swCustPropMgr = CmpDoc.Extension.CustomPropertyManager("")         lRetVal = swCustPropMgr.Get4("Наименование", False, Val, ValOut)
    • hlibhlib
      При двойном клике на имеющийся Размер для редактирования и последующем подтверждении действия слетают Настройки CAD на дефолт. ===> Бесит! Этого не было в предыдущих версиях. Может какие-то настойки нужны?  Сталкивался-ли кто с такой ситуацией?  Как выход из положения настроил линии и дуги в настройках Атрибутов. Вроде работает, но в меню Главная Цвет и ширина каркасной геометрии Дефолтная.  
    • Viktor2004
      Импортозамещение  
    • vad0000
      Нет, просто отключаете, и забываете о ней   Не надо ничего делать   Откуда приходит сигнал 24 вольта на 7 вход сервоусилителя? Фото другой страницы схемы
    • maxx2000
    • vladimir.songin
      То есть отключаю аналоговую обратную связь от аналогового входа, и подключаю вместо энкодера? Так не делал но сомнения у меня: Вход энкодера может  быть 12 или 5 вольт. У нас 12. Два пина питание, два пина Каналы A B соответственно, один пин REF (непонятно что он означает и как работает). Остальные соеденены вместе. На какие из них подавать аналаговый вход обратной связи дефлектора.  сигналы от энкодера приходят, проверял осцилографом на самом разьёме непосредственно на входе DM02, и некое значение уже могло бы отобразиться в поле энкодера при отключенном аналаговом входе.  Если отключить аналоговый вход, то видим единичку в обоих полях, но уже писал об этом.  
    • Алексей 1977
      Кто знает подскажите как отключить этот ненужный набор букв и символов в готовой УП? Я так думаю надо редактировать постпроцессор? Заранее спасибо ( Общая длина: 130.0) ( Заготовка:) ( MIN X: -10.970) ( MIN Y: -10.970) ( MIN Z: -6.500) ( MAX X: 10.970) ( MAX Y: 10.970) ( MAX Z: 0.000) ( COORDINATE SYSTEM: Глобальная СК) ( Кончик инструмента:) (   X: -0.000) (   Y: 0.000) (   Z: 10.000) ( Рекомендованная длина: 50.000) ( Количество кромок: 4) ( Инструмент:   Концевая фреза) ( DIAMETER: 10.000) ( Безопасность:) ( Рабочие ходы инструмента: Безопасная БЕЗ зарезов) ( Подводы инструмента: Безопасная БЕЗ зарезов) ( Переходы инструмента: Безопасная БЕЗ зарезов) ( Рабочие ходы патрона: Столкновения НЕ проверялись) ( Подводы патрона: Столкновения НЕ проверялись) ( Переходы патрона: Столкновения НЕ проверялись) ( Траектория: Шаблон) ( STEPOVER: 5.000) ( ДОПУСК:0.100) ( THICKNESS:0.000) ( Статистика:) ( LENGTH: 95.318)( LIFTS: ( TIME: 0/00/05) 1) G0X0Y0 G43Z10.H13 X4.75Y-8.227 Z5. G1Z0F500 X9.5Y-5.485F1000 Y5.485 X0Y10.97 X-9.5Y5.485 Y-5.485 X0Y-10.97 X4.75Y-8.227 G0Z10.
    • gudstartup
      считывание происходит при помощи вх\вых сигналов контроллера plc 
×
×
  • Создать...