Ц) Check for updates
УДК 624.04 Научная статья
https://doi.org/10.23947/2949-1835-2022-l-3-43-52
О влиянии выбора основной системы метода сил на точность расчетов при решении плоской задачи теории упругости в напряжениях
Г. Т. Гагалаян , А. С. Фоминых
Донской государственный технический университет, Российская Федерация, г. Ростов-на-Дону, пл. Гагарина, 1 Н ggagalayan@mail.ru
Аннотация
Введение. Известно, что решение задачи теории упругости возможно в перемещениях, напряжениях, а также смешанным методом. Наиболее распространен метод конечных элементов (МКЭ) в форме метода перемещений, однако он имеет свои недостатки. Целью настоящей статьи является изучение особенностей численного решения плоской задачи теории упругости в напряжениях с использованием метода конечных разностей (МКР) для возможности его использования в качестве альтернативы МКЭ. Рассматриваются варианты выбора основной системы метода сил для получения граничных условий при расчете прямоугольной балки-стенки методом сеток и их влияние на конечный результат при вычислении напряжений.
Материалы и методы. Расчет прямоугольной балки-стенки реализован в программе БаШа^внка, разработанной в ДГТУ. Выполняется сравнение решения на основе метода конечных разностей с решением методом конечных элементов в форме метода перемещений в программном комплексе «ЛИРА-САПР». Точность расчетов определяется путем выполнения статических проверок равновесия отсеченной части конструкции. Для рассматриваемой конструкции проанализированы пять вариантов выбора основной системы метода сил (ОСМС).
Результаты исследования. Впервые установлено, что, в отличие от расчета стержневых систем, выбор ОСМС несущественно влияет на окончательные величины напряжений при решении плоской задачи теории упругости.
Обсуждение и заключение. Сравнение результатов расчета при помощи метода конечных элементов в форме метода перемещений и метода конечных разностей показало некоторые преимущества последнего. Метод конечных разностей может выступить альтернативой методу конечных элементов в форме метода перемещений для решения плоской задачи теории упругости при условии автоматизации процесса построения эпюр внутренних усилий в контуре балки-стенки.
Ключевые слова: метод конечных элементов, метод конечных разностей, плоская задача теории упругости, балка-стенка, плоское напряженное состояние, рамная аналогия.
Для цитирования. Гагалаян, Г. Т. О влиянии выбора основной системы метода сил на точность расчетов при д
о
решении плоской задачи теории упругости в напряжениях / Г. Т. Гагалаян, А. С. Фоминых // Современные
тенденции в строительстве, градостроительстве и планировке территорий. — 2022. — Т. 1, № 3. — С. 43-52. https://doi.org/10.23947/2949-1835-2022-1-3-43-52
-а
ад
и
а £ л
© Г. Т. Гагалаян, А. С. Фоминых, 2022
Original article
About Influence of Choosing the Force Method Principal System on the Accuracy of Calculations in Solving a Stress-Based Plane Elasticity Problem
Gurgen T. Gagalayan ED, Anna S. Fominykh
Don State Technical University, Gagarin sq., 1, Rostov-on-Don, Russian Federation H ggagalayan@mail.ru
Abstract
Introduction. It is known that the solution of a problem in elasticity theory is possible based on displacements, stresses, as well as mixed method. The most common is the displacement-based finite element method (FEM), however, it has its imperfections. The aim of this article is to study the features of numerical solution of a stress-based plane elasticity problem using the finite difference method (FDM) as a possible alternative to the FEM. The options for choosing the force method principal system are considered aiming to obtain the boundary conditions in calculating a rectangular beam-wall by the grid method along with the influence thereof on the final result during stresses calculation. Materials and methods. The calculation of a rectangular beam-wall is made in Balka-Stenka software developed in DSTU. Comparison between a solution made by the finite difference method and a solution made by the displacement-based finite element method is made in the LIRA-SAPR software. The accuracy of calculations is determined by performing the static balance tests of the cut-off part of a structure. For the structure under consideration, the choice of five options of the force method principal system (FMPS) is analysed.
Results. For the first time the conclusion was made that in contrast to the calculation of bar systems, the choice of FMPS does not significantly affect the final stress values in solving a plane elasticity problem.
Discussion and conclusion. The comparison of calculation results made by the displacement-based finite element method and by the finite difference method showed some advantages of the latter. The finite difference method can act as an alternative to the displacement-based finite element method in solving a plane elasticity problem on condition of automation of the process of building the internal forces epures in a wall beam-wall contour.
Keywords: finite element method, finite difference method, plane stress problem, beam-wall, plane stress state, frame analogy.
For citation. G. T. Gagalayan, A. S. Fominykh. About Influence of Choosing the Force Method Principal System on the Accuracy of Calculations in Solving a Stress-Based Plane Elasticity Problem. Modern Trends in Construction, Urban and Territorial Planning, 2022, vol. 1, no. 3, pp. 43-52. https://doi.org/10.23947/2949-1835-2022-1-3-43-52
Введение. В качестве основного метода расчета строительных конструкций в настоящее время выступает
метод конечных элементов (МКЭ) в форме метода перемещений [1-4].
5 Недостатками МКЭ в форме метода перемещений являются [5-7]:
13
1. Пониженная, по сравнению с перемещениями, точность вычисления усилий (напряжений).
о
тз 2. Сложность учета разрывов значений усилий (напряжений) на межэлементных границах.
ад
3. Неиспользование граничных условий, выраженных в усилиях (напряжениях).
4. «Запирание» при расчете с учетом элементов, обладающих большой сдвиговой жесткостью, когда энергия деформации от сдвига значительно превышает энергию деформации от изгиба.
5. Содержащиеся в кинематических граничных условиях смещения конечного элемента, как жесткого целого, не связанные с работой внутренних сил. Они приводят к ухудшению обусловленности матрицы жесткости
и
всей конструкции при сгущении конечно-элементной сетки, когда необходимо повысить точность решения в опасных зонах или при итерационном решении нелинейных задач.
6. Сложности при включении в конструкцию недеформируемых элементов и при расчетах геометрически изменяемых систем в виде кинематической цепи.
Альтернативой МКЭ при решении плоской задачи теории упругости (ТУ) выступает метод конечных разностей. При этом в качестве искомой функции выступает функция напряжений.
Классическим объектом, для которого имеет место плоская задача ТУ, является балка -стенка — тонкая пластина, подверженная действию нагрузок, приложенных в ее срединной плоскости.
В качестве граничных условий при расчете балки-стенки выступает рамная аналогия. Если рассмотреть контур балки-стенки как раму, то функция напряжений на контуре будет равна изгибающему моменту, а ее производная по нормали - продольной силе:
р = М;
^ = , (1) дУ
Замкнутый контур балки-стенки при рассмотрении его как рамы является трижды статически неопределимым, однако эпюры М и N можно строить в основной системе метода сил (ОСМС). При этом в зависимости от выбора ОСМС эпюры будут различными, а значит и будет отличаться функция напряжений <р. Однако окончательные значения напряжений должны быть одинаковы, поскольку все функции напряжений, получаемые на основе различных ОСМС, отличаются на линейную функцию Ах + Ву + С. Известно, что для стержневых систем при решении методом сил точность результатов зависит от выбора ОСМС [8-9].
Целью настоящей работы выступает исследование влияния выбора ОСМС на точность определения напряжений в балке-стенке. В рамках поставленной цели решаются следующие задачи:
1. Расчет балки-стенки методом конечных разностей при различных вариантах выбора ОСМС.
2. Расчет той же конструкции при помощи метода конечных элементов в форме метода перемещений.
3. Сопоставление результатов расчета при помощи МКР в различных ОСМС, а также МКР и МКЭ. Материалы и методы. Задача определения напряжений деформированного состояния балки-стенки сводится к однородному дифференциальному уравнению четвертого порядка в частных производных [10-11]:
= 0, (2)
где ф — функция напряжений (функция Эри);
д4 д4 д4
V4 = —■ + 2 _ _ + ■— - двойной оператор Лапласа.
дх4 дх2дг2 дг4
Разностная аппроксимация уравнения (2) имеет вид:
(6 а2 + 8 а + 6)(рк — 4(1 + а)(ра + асрь + <рс + асра) + +2а(сре + фг + фд + фь) + ф1 + а2р1 + рт + а2рп = 0,
где а = йг2/йх2.
2
Разностная схема (3) имеет второй порядок точности. ^
Схема узлов к формуле (3) показана на рис. 1. д
о
Уравнение (3) составляется для всех внутренних узлов сетки. При составлении уравнений для смежных к контуру узлов возникает необходимость определения функции напряжений в «законтурных» точках. Для этого используется второе условие в (1).
Нормальные и касательные напряжения при отсутствии объемных сил вычисляются по формулам:
_ д2р _ д2р _ д2р (4)
д х д
ТЗ
ад
и и
а £ л
Или в форме метода конечных разностей:
°z,fc =
°х,к = '
ФЬ - 2Фк + Vd Ах2 '
Фа - 2<рк + фс 4z2 ;
(<Рд + <Ре) - (Vf + <Ph)
Txz,k ' 4AxAz
(5)
Рис. 1. Схема узлов для составления уравнений МКР (рисунок авторов)
В данной работе рассматривается прямоугольная балка-стенка (тонкая пластина) из изотропного материала, соответствующего закону Гука, загруженная нормальной к контуру поверхностной нагрузкой. При этом сосредоточенная сила ^ рассматривается как предельный случай равномерно распределенной нагрузки на малом участке балки-стенки. Толщина балки-стенки принимается равной единице.
Требуется определить напряженное состояние пластины. Расчетная схема приведена на рис. 2. Расчет
кН
выполняется при а = 9,9 м; Ь = 13,5 м; а= 30 —; Р= 130 кН; Дг = 1,1 м; Дг = 0,9 м. Конструкция разбивается на
м
9 интервалов по горизонтали и 15 — по вертикали.
р-
й о
ТЗ
ад
и
и CP
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
. 114, 115, . 116, . 117, . 118, ■ 119, , 120, - 121 , 122,
1 2 3 4 5 6 7 8 123 2Р
9 10 11 12 13 14 15 16 124
17 18 19 20 21 22 23 24 125 о о II 3
25 26 27 28 29 30 31 32 126
33 34 35 36 37 38 39 40 127 Ь= 13,5
41 42 43 44 45 46 47 48 128
49 50 51 52 53 54 55 56 129
57 58 59 60 61 62 63 64 130
65 66 67 68 69 70 71 72 1.11
73 74 75 76 77 78 79 80 132
81 82 83 84 85 86 87 88 133 2q
89 90 91 92 93 94 95 96 134
97 9В 99 100 101 102 103 104 135
105 106 107 108 109 110 111 112 136
153 154 155 156 157 158 159 160 137
■ р AX=1.l| ^ а = 9,9 г
Рис. 2. Исходная схема (рисунок авторов)
В качестве ОСМС рассматриваются варианты, представленные на рис. 3. Опоры в узлах 113 и 152 в ОСМС остаются нетронутыми, в вариантах (а) и (в) врезается три промежуточных шарнира в угловых точках, а в остальных вариантах делается разрез на одной из сторон.
Рис. 3. Варианты ОСМС (рисунок авторов)
Расчет реализован в программе ВаШа^епка, разработанной на кафедре «Сопротивление материалов» Донского государственного технического университета. Общее количество уравнений и неизвестных для рассматриваемой схемы — 112.
Рис. 4. Статическая проверка (рисунок авторов)
и Й О
-о
ад
и
СР £ Л
После определения напряжений выполняются две статические проверки равновесия отсеченной части. На рис. 4 показано равновесие верхней и правой отсеченной части при расчете на основе ОСМС варианта (а).
Также выполняется расчет при помощи МКЭ в форме метода перемещений в программном комплексе «ЛИРА-САПР» при том же количестве узлов сетки. При расчете в ПК «ЛИРА» используются плоские прямоугольные конечные элементы балки-стенки (КЭ 21) с двумя степенями свободы в узле (линейные перемещения и и w вдоль осей х и г).
Результаты исследования. В таблицах 1-2 представлены полученные в результате расчета методом конечных разностей с использованием различных ОСМС значения напряжений по линиям 144-129 (горизонтальное сечение) и 117 - 156 (вертикальное сечение). Также в указанных таблицах приведены результаты на основе МКЭ.
Таблица 1
Значения напряжений а2 и тХ2 по линии 144-129
й о -о
ад
и
и
а
Вариант Номер узла 144 49 50 51 52
а Oz -24,924 -37,231 -46,807 -50,975 -50,023
Txz 0 -11,445 -17,068 -17,384 -14,047
б Oz -24,946 -37,251 -46,818 -50,977 -50,019
Txz 0 -11,451 -17,076 -17,393 -14,054
в Oz -24,948 -37,251 -46,817 -50,977 -50,019
Txz 0 -11,451 -17,077 -17,393 -14,054
г Oz -25,17 -37,324 -46,783 -50,883 -49,918
Txz 0 -11,501 -17,133 -17,424 -14,05
д Oz -25,051 -37,284 -46,793 -50,931 -49,977
Txz 0 -11,487 -17,104 -17,399 -14,043
ЛИРА-САПР Oz -30,7 -36,1 -45,4 -50,5 -50,2
Txz -6,33 -10,8 -17,1 -17,9 -14,4
Вариант Номер узла 53 54 55 56 129
а Oz -44,513 -34,072 -17,462 6,209 34,671
Txz -8,575 -2,431 2,363 3,261 0
б Oz -44,506 -34,064 -17,453 6,219 34,685
Txz -8,58 -2,435 2,361 3,26 0
в Oz -44,506 -34,064 -17,453 6,219 34,685
Txz -8,58 -2,435 2,361 3,26 0
г Oz -44,433 -34,034 -17,471 6,152 34,559
Txz -8,549 -2,392 2,402 3,287 0
д Oz -44,48 -34,056 -17,463 6,191 34,638
Txz -8,561 -2,414 2,379 3,271 0
ЛИРА-САПР Oz -44,4 -32,9 -15 8,83 22
Txz -8,33 -1,81 2,49 2,62 1,54
Таблица 2
Значения напряжений ах и тх2 по линии 117-156
Вариант Номер узла 117 4 12 20
а Ох -125,231 -99,051 -64,507 -34,823
0 -1,958 -5,954 -9,754
б Ох -125,212 -99,038 -64,499 -34,821
0 -1,964 -5,962 -9,763
в Ох -125,213 -99,038 -64,499 -34,821
0 -1,963 -5,962 -9,763
г Ох -125,251 -99,067 -64,526 -34,849
0 -1,965 -5,965 -9,766
д Ох -125,225 -99,049 -64,51 -34,834
0 -1,964 -5,963 -9,764
ЛИРА-САПР Ох -114 -99,8 -88,3 -37,3
-0,582 -1,86 -5,22 -9,33
Вариант Номер узла 28 36 44 52
а Ох -14,398 -2,907 1,454 0,406
-12,224 -13,406 -13,821 -14,047
б Ох -14,398 -2,911 1,448 0,397
-12,233 -13,414 -13,829 -14,054
в Ох -14,398 -2,911 1,448 0,397
-12,232 -13,414 -13,829 -14,054
г Ох -14,429 -2,942 1,422 0,386
-12,235 -13,415 -13,827 -14,05
д Ох -14,414 -2,929 1,429 0,383
-12,233 -13,413 -13,824 -14,043
ЛИРА-САПР Ох -14,6 -1,88 1,76 0,61
-12,6 -13,9 -14,3 -14,4
Вариант Номер узла 60 68 76 84
а Ох -4,885 -13,717 -25,506 -39,38
-14,564 -15,675 -17,342 -18,897
б Ох -4,897 -13,732 -25,523 -39,396
-14,57 -15,679 -17,342 -18,89
в Ох -4,897 -13,732 -25,523 -39,396
-14,57 -15,679 -17,342 -18,89
г Ох -4,877 -13,667 -25,408 -39,247
-14,568 -15,688 -17,379 -18,962
д Ох -4,898 -13,706 -25,455 -39,283
-14,551 -15,653 -17,318 -18,892
ЛИРА-САПР Ох -4,14 -13,4 -25,8 -40,6
-15 -16,4 -18,5 -20
Вариант Номер узла 92 100 108 156
а Ох -53,872 -66,502 -73,431 -70,302
-18,924 -15,761 -8,761 0
б Ох -53,881 -66,496 -73,404 -70,264
-18,91 -15,741 -8,744 0
в Ох -53,881 -66,496 -73,404 -70,264
-18,91 -15,741 -8,744 0
г Ох -53,735 -66,382 -73,386 -71,111
-18,983 -15,723 -8,633 0
д Ох -53,759 -66,469 -73,571 -70,422
-18,968 -15,842 -8,732 0
ЛИРА-САПР Ох -55,4 -86,9 -71,8 -72,2
-18,8 -14,1 -7,1 -3,15
В таблицах 3-4 представлены полученные погрешности статических проверок равновесия левой и верхней отсеченной части при различных вариантах ОСМС, а также при расчете с помощью МКЭ в программном комплексе «ЛИРА-САПР».
Таблица 3
Погрешности статической проверки равновесия левой отсеченной части
Погрешность IZ ЕМ
а) 0,313 % 0,478 % 0,0805 %
б) 0,001 % 0,476 % 0,133 %
в) 0,001 % 0,476 % 0,133 %
г) 0,0036 % 0,461 % 0,0827 %
д) 0,00103 % 0,465 % 0,16 %
ЛИРА-САПР 0,0405 % 0,611 % 0,522 %
Таблица 4
Погрешности статической проверки равновесия верхней отсеченной части
Погрешность EZ ЕМ
а) 0,979 % 0,0017 % 0,454 %
б) 0,99 % 0,73 % 0,464 %
в) 0,99 % 0,73 % 0,464 %
г) 1,007 % 0,0056 % 0,453 %
д) 0,988 % 0,0017 % 0,462 %
ЛИРА-САПР 0,942 % 0,0069 % 0,542 %
и SP
Из таблиц 1-2 видно, что напряжения при расчете методом конечных разностей при выборе основной системы метода сил по схемам «а-д» отличаются незначительно. Среднее расхождение результатов, полученных при помощи МКР, с решением на основе метода конечных элементов составляет 13,9 % по напряжениям ах, 7,2 % по напряжениям а2 и 7,7 % по напряжениям тХ2 (без учета крайних узлов). Следует также обратить внимание на то, что в лежащих на контуре балки-стенки узлах 144, 129, 117, 156 при расчете методом конечных элементов касательные напряжения тХ2 оказываются отличными от нуля, что противоречит условиям на поверхности:
Xv = ax-l + Txz- п, zv — Tzx ' I + Gz ' П,
(6)
где Ху и составляющие поверхностной нагрузки по осям х и г; 1 = cos(x,v); п = cos(z,v) - направляющие косинусы нормали к поверхности балки-стенки V.
Невыполнение условий на поверхности по касательным напряжениям при расчете в ПК «ЛИРА» можно объяснить тем, что напряжения в данном программном комплексе вычисляются только в центрах конечных элементов, и выведенные значения для узлов 144, 129, 117, 156 на самом деле соответствуют точкам, смещенным на Ах/2 или А212.
В то же время эта особенность никак не отразилась на погрешностях статических проверок, они сопоставимы с результатами, полученными при помощи метода конечных разностей региона [9].
й о
^ Обсуждения и заключения. Полученные результаты позволяют сделать следующие выводы:
1. Выбор основной системы метода сил практически не влияет на результат по напряжениям.
2. Точность решения при помощи метода конечных элементов сопоставима с точностью метода конечных разностей, но при решении методом конечных элементов при одинаковой густоте сетки число неизвестных в узле в 2 раза больше (перемещения и и w).
3. Недостатком МКЭ по сравнению с МКР является невозможность точного определения напряжений непосредственно на поверхности конструкции.
4. МКР может служить альтернативой МКЭ в форме метода перемещений при решении плоской задачи теории упругости в случае автоматизации процесса построения эпюр изгибающего момента и продольной силы в контуре балки-стенки, как в раме.
Библиографический список
1. Фокин, В. Г. Метод конечных элементов в механике деформируемого твёрдого тела / В. Г. Фокин. — Самара: Самарский гос. технический ун-т, 2010. — С. 22. URL: http ://meh. samgtu.ru/sites/meh. samgtu.ru/files/fokin.pdf
2. Языев, Б. М. Расчёт трёхслойной пластинки методом конечных элементов с учётом ползучести среднего слоя / Б. М. Языев и др. // Вестник Дагестанского государственного технического университета. Технические науки. — 2014. — Т. 33, №. 2. — С. 47-55.
3. Chepurnenko, A. Finite element modeling of the creep of shells of revolution under axisymmetric loading / A. Chepurnenko, N. Neumerzhitskaya, M. Turko // Energy Management of Municipal Transportation Facilities and Transport. — Springer, Cham, 2017. — P. 808-817. DOI: https://doi.org/10.1007/978-3-319-70987-1 86.
4. Chepurnenko, A. Calculation of a three-layer plate by the finite element method taking into account the creep of the filler / A. Chepurnenko, A. Savchenko, S. Yazyeva // MATEC Web of Conferences. — EDP Sciences, 2017. — Vol. 129. — P. 05008. DOI: https://doi.org/10.1051/matecconf/201712905008.
5. Игнатьев, А. В. Развитие метода конечных элементов в форме классического смешанного метода строительной механики / А. В. Игнатьев // Дисс. д-ра техн.технических техн. 543 наук: 05.23.17. — Волгоград, 2019.
— С. 294. — URL: https://www.vstu.ru/upload/iblock/65a/65a405ea74517bf03e889045eda7ad00.pdf.
6. Игнатьев, А. В. Применение метода конечных элементов в форме классического смешанного метода к расчету систем с односторонними связями / А. В. Игнатьев, В. А. Игнатьев, М. И. Бочков // Строительная механика и расчет сооружений. — 2017. — № 2. — С. 52-61.
7. Игнатьев, А. В. Расчет тонких пластин по методу конечных элементов в форме классического смешанного метода с исключением перемещений конечных элементов как жесткого целого / А. В. Игнатьев, В. А. Игнатьев, Е. А. Гамзатова // Известия высших учебных заведений. Строительство. — 2018. — №. 3. — С. 5-13.
8. Радкевич, Р. А. О влиянии выбора варианта основной системы метода сил для бесшарнирной арки на устойчивость решения канонических уравнений / Р. А. Радкевич // Строительные конструкции. — Новополоцк: ПГУ, 2015.
— С. 77. URL: https://elib.psu.bv:8080/bitstгeam/123456789/15391/5/Трvлы%20мололых 2015%20Вып.9 СТРОИ-ТЕЛЬНЫЕ%20КОНС'ТРУКЦИИ.pdГ#paцe=9
9. Романова, М. В. О влиянии выбора варианта основной системы метода сил на достоверность расчета рамных конструкций / М. В. Романова // Строительные конструкции. — Новополоцк: ПГУ, 2014. — P. 52. — URL: https://elib.psu.by/bitstream/123456789/17465/5/3 74 %20Строительные %20констрvкции.pdf#page=4 g
10. Kozelskaya, M. Stress strain behavior research of triangular dam using analytical and numerical methods / «
й
M. Kozelskaya, D. Donskova, V. Ulianskaya, P. Shvetsov // IPICSE-2018. — 2018. — No 251. — P. 8.
DOI: https://doi.org/10.1051/matecconf/201825104043.
11. Донскова, Д. О. Особенности расчета треугольной плотины численным и аналитическим методом / Д. О. Донскова // Молодой исследователь Дона. — 2019. — № 2 (17). — С. 9. URL: https://cvberleninka.ru/arti-cle/n/osobennosti-rascheta-treugolnoi-plotiny-chislennym-i-analiticheskim-metodom.
ад
и и
а £ л
Поступила в редакцию 18.11.2022
Поступила после рецензирования 25.11.2022
Принята к публикации 30.11.2022
Об авторах:
Гагалаян Гурген Татеосович — студент факультета «Промышленное и гражданское строительство» Донского государственного технического университета (344003, РФ, г. Ростов-на-Дону, пл. Гагарина, 1), ORC1D. ggagalavan@mail.ru, тел.: +7 (950) 865-43-26.
Фоминых Анна Сергеевна — студент факультета «Промышленное и гражданское строительство» Донского государственного технического университета (344003, РФ, г. Ростов-на-Дону, пл. Гагарина, 1). ORC1D. fomanva03@gmail.com, тел.: +7 (908) 517-41-66.
Заявленный вклад авторов:
А. С. Фоминых — проведение расчетов в программе Balka-Stenka, подготовка текста статьи. Г. Т. Гагалаян — проведение расчетов в ПК «ЛИРА», подготовка графической части, редактирование текста статьи.
Конфликт интересов
Авторы заявляют об отсутствии конфликта интересов.
Все авторы прочитали и одобрили окончательный вариант рукописи.
й о -о
ад
и
и
а