[hungry_Duck]

Добрый день!
В попытке реализовать программно расчет внецентренного сжатия по НДМ столкнулся с рядом проблем.
Предыдущую тему "Откуда пошла..." проштудировал и взял за основу алгоритм пользователя palexxvlad.

Хотелось бы спросить у тех, кто уже реализовывал этот алгоритм, сразу ли он заработал или пришлось над чем-то шаманить? Есть ли там какие-то подводные камни?
У меня пока получается так, что решение сходится к неверным результатам, из-за чего проверка по уравнениям равновесия не проходит. При этом коэффициенты Nu почти во всех участках равны либо 0 либо 1. 0 в растянутой зоне, где напряжения принимаю равными 0, 1 - в зоне, где материал работает линейно. "Нелинейную зону" увидеть пока толком не получается .

Банальные ошибки устранил, и теперь в раздумьях куда дальше двигаться и на что смотреть.
Как правильно поступать с напряжением Sigma_bi, если отн. деформации превышают значение Epsilon_b2 - крайнюю точку на диаграмме? Сейчас оно у меня принимается = 0.

Пишу на Python3.0. Разбиение сечения на элементарные участки провожу сторонней программой GMSH.

Вобщем, буду рад любым комментариям/советам по-существу.

[CEBEP]

https://dwg.ru/dnl/13011 - бери да смотри.там всё видать

[frostyfrost]

Сам писал на третьем Питоне+Numpy+PIL веб-приложение (https://civileng.ru/check/rc/ndm-rect-sym/ и другие). Алгоритм лучше брать представленный в работе Симбиркина с меняющимся приведенным модулем упругости в зависимости от величины напряжений и деформаций (подставляя нужные кривые деформирования). Для крошения в зависимости от НДС деформации могут ограничиваться менее eps_b2. Для себя учитывал диаграмму работы бетона на растяжение. Если ушли за eps_b2 и eps_bt2, то напряжения обнуляете (крошение и окончание участка разупрочнения), но не деформации (модуль упругости в этом случае лучше снизить на несколько порядков, не обнуляя его). Также рекомендую разбивать на 5-10 шагов нагрузку, чтобы экономить на итерациях. В сложных случаях у меня, если сечение целиком не разрушается, получается не более 10-30 итераций.

[tutanhamon]

Сообщение от frostyfrost Сам писал на третьем Питоне+Numpy+PIL веб-приложение (https://civileng.ru/check/rc/ndm-rect-sym/ и другие). Алгоритм лучше брать представленный в работе Симбиркина с меняющимся приведенным модулем упругости в зависимости от величины напряжений и деформаций (подставляя нужные кривые деформирования). Для крошения в зависимости от НДС деформации могут ограничиваться менее eps_b2. Для себя учитывал диаграмму работы бетона на растяжение. Если ушли за eps_b2 и eps_bt2, то напряжения обнуляете (крошение и окончание участка разупрочнения), но не деформации (модуль упругости в этом случае лучше снизить на несколько порядков, не обнуляя его). Также рекомендую разбивать на 5-10 шагов нагрузку, чтобы экономить на итерациях. В сложных случаях у меня, если сечение целиком не разрушается, получается не более 10-30 итераций.

Интересно. Вы не думали про реализацию произвольного сечения?

[frostyfrost]

У меня все для этого уже сделано. Я сразу делал с прицелом на любые сечения. На вход в расчете подается "файл" сетки для бетона и арматуры (диаметр и координаты, что позволяет потом считать трещины). Единственное ограничение, которое есть пока, я не нашел нормальную библиотеку мешера, чтобы в условиях сервера выйти на приемлемые затраты времени. Все сетки пока генерирую предварительно, т.к. PIL съедает кучу времени на генерацию графики, поэтому, кстати, пока отказался от визуализации деформаций (закомментировал эти блоки на будущее).

С сечениями по НДМ другие проблемы. Есть сечения типа треугольника, которые по НДМ проходят, а в объемном исполнении с учетом поверхностей разрушения для объемного НДС не проходят.

[tutanhamon]

У меня все для этого уже сделано. Я сразу делал с прицелом на любые сечения. На вход в расчете подается "файл" сетки для бетона и арматуры (диаметр и координаты, что позволяет потом считать трещины). Единственное ограничение, которое есть пока, я не нашел нормальную библиотеку мешера, чтобы в условиях сервера выйти на приемлемые затраты времени. Все сетки пока генерирую предварительно, т.к. PIL съедает кучу времени на генерацию графики, поэтому, кстати, пока отказался от визуализации деформаций (закомментировал эти блоки на будущее).

Не специалист в этой области, может быть есть какие-то проблемы в области реализации взаимодействия с сервером, но ведь для пользователя же не так критичны затраты по времени? То есть подождать секунд 10 или 30 - большой разницы для пользователя не будет, пока генерируется сетка, передается в "решатель", происходит итерационный расчет и т.п.
Offtop: Не исключаю, что написал глупость тут, ибо знания мои в области Web-программирования весьма скудны

Сообщение от frostyfrost С сечениями по НДМ другие проблемы. Есть сечения типа треугольника, которые по НДМ проходят, а в объемном исполнении с учетом поверхностей разрушения для объемного НДС не проходят.

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

[frostyfrost]

С испытаниями не сталкивались, к сожалению, но при моделировании с применением различных моделей бетона Willam-Warnke, Ottosen и др. получались близкие результаты по картине разрушения с откалыванием острой части треугольника. Я сейчас в отпуске, когда буду снова в Мск, то возможно статью в блог на dwg кину.

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

[tutanhamon]

Сообщение от frostyfrost С испытаниями не сталкивались, к сожалению, но при моделировании с применением различных моделей бетона Willam-Warnke, Ottosen и др. получались близкие результаты по картине разрушения с откалыванием острой части треугольника. Я сейчас в отпуске, когда буду снова в Мск, то возможно статью в блог на dwg кину.

Спасибо! Я думаю, многие будут Вам за это благодарны.

[hungry_Duck]

Да! Было бы здорово почитать!

----- добавлено через ~3 ч. -----

Сообщение от frostyfrost Алгоритм лучше брать представленный в работе Симбиркина

А что за работа? Я к сожалению, такой не знаю.
tutanhamon, не пойму такой момент в алгоритме palexxvlad: итеративный процесс организован за счет вычисления поправочных коэффициентов Nu_bi и Nu_sj. Но, пока мы находимся в зоне упругой работы материала, эти коэффициенты равны 1 и процесс сходимости останавливается. Как быть?

[hungry_Duck]

И еще один вопрос к tutanhamon, по поводу макроса.
При двухлинейной диаграмме работы бетона, при высчитывании коэффициента Nu_b вы делите на значение Eb, а в пункте 8.1.23 говорится, что надо брать E_b,red. Для чего у вас так сделано?

[Akim_1989]

hungry_Duck, давно писал посмотри как пример.Там конечно не без косяков)). Я использовал scipy и numpy. Диаграмму бетона и арматуры представлял в виде объекта spline из библиотеки scipy. А дальше там просто перемножение массивов в цикле.
https://yadi.sk/d/2BgVi0UpubKrM

[hungry_Duck]

Да, спасибо!
Дело с мертвой точки сдвинулось, работаю над визуализацией результатов.

Offtop: Как проекту надо дать "отлежаться", так и на код порой полезно взглянуть свежими глазами.

[frostyfrost]

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

[hungry_Duck]

frostyfrost, Спасибо, попробую.
Сейчас наваяю блок визуализации, с ним дело пойдет интереснее.

[Akim_1989]

Для визуализации удобная библиотека matplotlib

[hungry_Duck]

Я привык работать с VTK, она тяжелая и сложновата, но возможности практически неограниченные.
matplotlib посмотрю, что за зверь такой.

[румата]

У Хасана Мамаева есть код на пайтоне, реализующий "алгоритм palexxvladа" для прямоугольного сечения. У них была задумка сделать обсчет колонн в ревите через динамо на основе НДМ и робот.

----- добавлено через ~2 мин. -----
Если разрешат, выложу.

[eilukha]

Сообщение от frostyfrost в работе Симбиркина

- подробнее можно?

[frostyfrost]

Сообщение от eilukha - подробнее можно?

см. вложение.

Для питона хорошая графическая библиотека PIL. Из минусов для меня, слишком долго рендерит много текста в условиях сайта.

[Boxa]

Если поможет, то вот тут http://www.effect-project.ru/article...oi-modeli.html разобрано очень подробно и наглядно, и очень помогает на начальном этапе... ну а дальше уже накрутить можно любое сечение...
Мне эта статья, в свое время, очень помогла быстро включится.