В июле 2018 г. вышла версия Штуцер-МКЭ 2.15 [1], в которой произошло несколько серьезных изменений. Мы полностью переписали ядро программы, чтобы оно отвечало таким современным требованиям, как быстродействие, расширяемость, поддержка иерархических моделей на логическом уровне и многое другое.
Структура ядра состоит из трех основных компонент (рис. 1), каждая из которых может быть использована независимо друг от друга. Объектно-иерархическая модель представляет собой конструкцию в виде дерева объектов (например, штуцер является дочерним объектом несущей обечайки). В этой модели задаются только те геометрические данные, которые описывают сам объект, а также необходимые для привязки дочерних объектов.
Из объектно-иерархической генерируется геометрическая модель, которая является топологической, то есть содержит информацию о кривых (curves), поверхностях (surfaces), объемных объектах (solids), а также о связях между ними. Этот компонент использует лицензированное геометрическое ядро C3D [4]. На данный момент в ядре Штуцер-МКЭ применяется оболочечная топологическая модель, которая строится по серединной поверхности обечаек.
На основе геометрической формируется конечно-элементная модель (КЭ-модель), состоящая из оболочечных конечных элементов. Из нее непосредственно составляется система линейных алгебраических уравнений, соответствующая условиям равновесия деформируемого твердого тела.
В настоящее время ядро работает на 32- и 64-разрядных платформах. Для расширения возможностей конечно-элементный решатель был переработан под 64-разрядные платформы. Поэтому, начиная с версии 2.15, программа устанавливается только на 64-разрядные версии Windows OS.
В версии 2.15 были уточнены методики базовых расчетов прочности, допускаемых нагрузок на штуцер, жесткостей, а также коэффициентов интенсификации напряжений. При расчете жесткостей теперь используются два подхода: расчет жесткости системы «несущая обечайка — штуцер» (общая жесткость) и расчет жесткости за исключением балочной (местная жесткость в месте врезки штуцера). Второй подход применяется при совместной работе с программами, использующими балочную конечно-элементную модель (такими как СТАРТ [3]) с некоторым моделированием оборудования (например, обечайка сосуда или тройник в СТАРТ [3]). В этом случае Штуцер-МКЭ строит «за кадром» эквивалентную балочную модель, а местная «оболочечная» жесткость врезки определяется по разности перемещений оболочечной и балочной моделей:
где
- расчетное значение i-й степени свободы на конце штуцера в балочной конечно-элементной модели,
а 
- аналогичное значение в оболочечной модели.
Для уточнения методики расчета напряженно-деформированного состояния была реализована процедура экстраполяции напряжений, позволяющая решать следующие задачи:
Процедура экстраполяции напряжений реализована методом Hot Spot Stress (HSS) [5−8]. Суть этого метода заключается в том, что расчетные напряжения в сечении, где определяются их максимальные значения (рис. 2), вычисляются по линейной экстраполяции в точках t1 и t2 (Linear Surface Extrapolation — LSE), выбор которых зависит от толщины обечайки s. Такой подход повышает сходимость результатов: делает устойчивым процесс определения напряжений в местах концентраторов и позволяет получить близкие результаты при разных уровнях разбивки КЭ-модели.
Приведем сравнение результатов расчетов с использованием процедуры экстраполяции и без нее на примере врезки штуцера в цилиндрическую обечайку. Пусть на штуцер действует только осевая сила 100 кН. В этом случае при предыдущем подходе максимальные условно-упругие напряжения могут не определиться даже на 5-м уровне сетки (рис. 3б). Здесь суммарные напряжения (мембранные и изгибные) на внешней поверхности на 1-м уровне составили 64,1 МПа (рис. 3а), а с учетом коэффициента разбивки Km=1,30 расчетные будут равны 83,3 МПа. В то же время напряжения на 5-м уровне составили 83,9 МПа (рис. 3б), а с учетом коэффициента разбивки Km=1,05 расчетные будут равны 88,1 МПа. Отсюда видно, что максимальное значение условно-упругих напряжений на 1-м уровне не достигает значения на 5-м уровне. Ниже будет показано, что максимальные пиковые значения не получены и на 5-м уровне.
|
|
Процедура экстраполяции напряжений, примененная в данном примере (рис. 4), дала для 1-го и 5-го уровней разбивки схожие результаты, которые при этом превысили результаты, полученные без использования экстраполяции напряжений.
|
|
В версии программы Штуцер-МКЭ 2.16 предложен новый подход к типу расчетов «в месте врезки». Однако пользователи нередко не понимали, как рассчитываются напряжения в околошовной зоне. Вариант «в месте врезки» обычно используется, когда подводящий трубопровод заканчивается в районе врезки в обечайку (точка «i» на рис. 5). При этом расчетная схема трубопровода, как правило, заканчивается на внешней поверхности обечайки, а расчетные нагрузки на штуцер получены в данной точке.
В предыдущих версиях для построения конечно-элементной модели принималась минимально возможная расчетная длина штуцера, что приводило к более жесткой врезке трубопровода в обечайку, некорректному определению напряжений в районе врезки и, соответственно, к риску проявления краевых эффектов от места приложения нагрузки на штуцер. Однако, начиная с версии 2.16, длина штуцера определяется из условия достаточной гибкости:
Эта длина, полученная на основе многочисленных натурных испытаний [9], позволяет учитывать овализацию сечения в районе врезки, что соответствует «гибкому» соединению трубопровода с обечайкой (рис. 6).
|
|
|
На графике изменения допускаемой осевой силы Fy в зависимости от длины штуцера L1 (рис. 7) видно, что в районе врезки на расстоянии меньше 100 мм проявляются краевые эффекты от места приложения нагрузки. Результаты расчета представлены для косой врезки штуцера с внутренним диаметром 203 мм и толщиной 6 мм в цилиндрическую обечайку с внутренним диаметром 2000 мм и толщиной 12 мм.
При расчете на прочность программный комплекс Штуцер-МКЭ автоматически переводит нагрузки, заданные пользователем в месте врезки, в статически эквивалентные, соответствующие сечению с длиной L1 = Lmin (рис. 5).
Переработка ядра, а также использование современного геометрического ядра твердотельного моделирования C3D позволили добавить новые типы врезок: с отбортовкой и вварной торообразной вставкой (рис. 8). Теперь в Штуцер-МКЭ можно моделировать штампованные и штампосварные узлы врезок.
| В отличие от ПАССАТ [2], в Штуцер-МКЭ долгое время не было конического перехода со смещением второго конца (эксцентрический переход) (рис. 9), что было обусловлено сложностью реализации этой возможности в рамках старого ядра программы. Многие пользователи просили добавить в Штуцер МКЭ плоские прямоугольные крышки, с помощью которых можно моделировать врезки в крышки АВО (аппаратов воздушного охлаждения), стенки кубических емкостей |
|
Для повышения корректности соответствия моделируемых нагрузок реальным была добавлена возможность задания гидростатического давления, необходимого, в первую очередь, при расчете вертикальных сосудов и резервуаров, хотя эту нагрузку можно задавать и для других типов обечаек (рис. 11). Гидростатическое давление сейчас задается линейно вдоль вертикальной оси, где максимальное значение определяется на уровне самой нижней точки модели. Для вертикальных сосудов (резервуаров) это соответствует уровню днища.
Востребованным и необходимым дополнением к программе, реализованным в версии 2.16, стал учет температурных деформаций (расширение материалов элементов конструкции под действием температуры, а также дополнительные напряжения, возникающие при стесненных температурных деформациях — разные коэффициенты линейного расширения, разные температуры
На рис. 12а приведено распределение эквивалентных напряжений в эллиптическом днище, в котором действует только внутреннее давление 0,5 МПа. Максимальное значение напряжений составило 86,6 МПа. Для примера зададим материалы на обечайке и штуцере с разными коэффициентами линейного температурного расширения, а также разные температуры на днище (200 °С) и штуцере (220 °С). В результате мы получим стесненность температурных деформаций и в районе врезки штуцера напряжения резко возрастут (рис. 12б), при этом максимальное значение напряжений в данном случае составило 118,5 МПа.
|
|
Кроме того, учет температурного воздействия позволил обеспечить более полное соответствие требованиям нормативных документов. Например, в ASME BPVC.VIII.2 при проверке эквивалентных напряжений категории «Secondary Membrane plus Membrane» необходимо учитывать термическое расширение материалов. А в
или в проверке по формуле (13) — расчетные температуры, при которых допускаемые напряжения устанавливают по пределам длительной прочности или ползучести:
При разработке версии нового ядра и его доработке был также усовершенствован модуль проверок на прочность в соответствии с различными нормативными документами, что позволило достаточно гибко и в то же время быстро актуализировать как зарубежные, так и отечественные стандарты до версий ASME BPVC.VIII.1, 2 2017, EN 13445 3:2014 (E) Issue 4 (2017 07),
В перспективе мы представляем Штуцер-МКЭ, с одной стороны, как интерактивный конструктор сосудов и аппаратов, а с другой — как встраиваемое решение. Поэтому первоочередной задачей является переработка интерфейса программы. В результате, например, Штуцер-МКЭ сможет «бесшовно» интегрироваться c программами ПАССАТ и СТАРТ для расчета врезок непосредственно в их среде.
Не исключено, что в итоге получатся и будут одновременно существовать два программных продукта: один — в виде отдельной программы с расширенным набором средств задания топологических моделей, а другой — в качестве аналога текущей версии для создания топологических моделей с возможностями встраивания в другие продукты. При этом оба решения будут обладать большей, чем сейчас, функциональностью и интерактивностью.
Литература
