Сглаживание внутренних сил поверхности в программе RFEM

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

Пример: Пластина по Миндлину и Кирхгофу

Пример:

Сначала создадим двумерную модель пластины размером 4 м ⋅ 3 м. Пластина должна быть с обеих сторон шарнирно закреплена по более коротким граничным линиям. Более длинные граничные лини защемлены против кручения. Чтобы ограничить действие поперечной деформации, установим для материала нулевой коэффициент Пуассона. Для сравнения внутренних сил, определенных с помощью теории изгиба пластин по Кирхгофу либо по Миндлину, зададим толщину пластины 80 см. Таким образом, мы получим соотношение d/L = 0,2, которое находится в пределах действия обеих теорий пластин (см. https://www.dlubal.com/ru/podderzhka-i-obuchenije/podderzhka/faq/003158).

Для простоты шаг сетки КЭ задан равным 1 м. Таким образом мы получим 4 ⋅ 3 конечных элемента.

На пластину действует нагрузка p = 5 кН/м². Собственный вес в нагружении не учитывается.

Не сглаженные внутренние силы

Результаты, полученные решателем для внутренних сил, v_x и m_x относятся к продольному сечению, проходящему через центр пластины. Распределение не сглаженных внутренних сил по элементам можно отобразить с помощью функции «Не непрерывно». Распределение по теориям пластин Миндлина и Кирхгофа показано на рисунке 02.

Несглаженные результаты по Миндлину (слева) и Кирхгофу (справа)

Стандартное сглаживание для внутренних узлов

Для узлов сетки КЭ, лежащих внутри поверхности, сначала определяется среднее арифметическое из узловых значений смежных конечных элементов. В данном подходе элементы должны лежать на поверхности и на одной и той же стороне потенциальной внутренней линии. Узел не должен быть узлом поверхности, заданным пользователем.

При положении разрыва в точке x = 1,00 м данные внутренние значения являются следующими (результат данного шага недоступен для вывода).

Миндлин:

Кирхгоф:

Сглаживание для краевых узлов

На краях поверхности (функция сглаживания «Непрерывно внутри поверхности») или модели («Непрерывно в целом») нет смежных конечных элементов, которые можно было бы применить для усреднения узловых значений. Поэтому используется другой метод, который выполняется в два этапа.

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

Значения краевых узлов в x = 0,00 м являются следующими.

Миндлин:

Кирхгоф:

Силы сдвига

По Миндлину сдвигающие силы рассчитываются через первую производную прогиба.

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

В теории изгиба пластин по Кирхгофу сдвигающие силы рассчитываются через третью производную прогиба.

Такое определение значения с помощью третьей производной приводит к значительной потере точности. По этой причине сдвигающие силы, рассчитанные решателем, не используются, а определяются по улучшенному методу с помощью дифференцирования моментов.

Моменты m_x, m_y и m_xy в уравнениях представляют собой сглаженные значения, которые определяются по методам, описанным выше. В результате получены более точные значения сдвигающих сил, чем из расчета с помощью решателя.

Моменты

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

Данный подход соответствует указанным выше уравнениям для вычисления сдвигающих сил через производные моментов. Затем моменты определяются с помощью следующих уравнений.

Константы интегрирования m_x,0 и m_y,0 рассчитываются по условиям значения в центре стержня.

S представляет собой поверхность элемента.

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

По теории изгиба пластин Миндлина в нашем примере мы получим следующие значения.

Первый элемент:

Второй элемент:

Сглаженное распределение внутренних сил v-x и m-x по Миндлину

По теории изгиба пластин Кирхгофа мы получим следующие значения.

Первый элемент:

Второй элемент:

Сглаженное распределение внутренних сил v-x и m-x по Кирхгофу

В процессе сглаживания в программе последовательно выполняются описанные выше шаги. Сглаженные значения можно отобразить графически с помощью функции отображения «Непрерывно внутри поверхности».