Научный отдел
МЕХАНИКА
Известия Саратовского университета. Новая серия. Серия: Математика. Механика. Информатика. 2022. Т. 22, вып. 1. С. 48-61 Izvestiya of Saratov University. Mathematics. Mechanics. Informatics, 2022, vol. 22, iss. 1, pp. 48-61 https://mmi.sgu.ru
https://doi.org/10.18500/1816-9791-2022-22-1-48-61
Научная статья УДК 539.3
Нелинейное деформирование осесимметрично нагруженной оболочки вращения на основе МКЭ при различных вариантах определяющих уравнений
А. Ш. Джабраилов10, А. П. Николаев1, Ю. В. Клочков1, Н. А. Гуреева2, Т. Р. Ищанов1
1 Волгоградский государственный аграрный университет, Россия, 400002, г. Волгоград, просп. Университетский, д. 26 2Финансовый университет при Правительстве Российской Федерации, Россия, 125993, г. Москва, Ленинградский просп., д. 49
Джабраилов Арсен Шахнавазович, кандидат технических наук, доцент, доцент кафедры высшей математики, arsen82@yandex.ru, https://orcid.org/0000-0001-6494-1377
Николаев Анатолий Петрович, доктор технических наук, профессор, профессор кафедры прикладной геодезии, природообу-стройства и водопользования, anpetr40@yandex.ru, https://orcid.org/ 0000-0002-7098-5998
Клочков Юрий Васильевич, доктор технических наук, профессор, заведующий кафедрой высшей математики, Klotchkov@bk.ru, https://orcid.org/0000-0002-1027-1811
Гуреева Наталья Анатольевна, доктор физико-математических наук, доцент, доцент департамента анализа данных, принятия решений и финансовых технологий, Natalya-gureeva@yandex.ru, https: //orcid.org/0000-0003-3496-2008
Ищанов Тлек Рахметолович, кандидат технических наук, доцент кафедры высшей математики, ishchanov.volgau@yandex.ru, https: //orcid.org/0000-0002-2137-4269
Аннотация. Использован криволинейный конечный элемент срединной линии осесимметрично нагруженной оболочки вращения с матрицей жесткости размером 8x8 при выборе узловых неизвестных в виде перемещений и их первых производных. Определяющие уравнения на шаге нагружения реализованы в двух вариантах. В первом варианте использованы соотношения деформационной теории пластичности, состоящие из выражений упругих и пластических частей. Соотношения между приращениями деформаций и приращениями напряжений определялись дифференцированием используемых уравнений. Во втором варианте гипотеза о разделении деформации
© Джабраилов А. Ш, Николаев А. П., Клочков Ю. В., Гуреева Н. А, Ищанов Т. Р., 2022
на упругую и пластическую части не использовалась. Разработанные авторами определяющие уравнения получены на основе предложенной гипотезы о пропорциональности компонент девиаторов приращений напряжений и компонент девиаторов приращений деформаций с коэффициентом пропорциональности в виде функции хордового модуля диаграммы деформирования. Представлен пример расчета, показывающий эффективность разработанного алгоритма.
Ключевые слова: оболочка вращения, тензор деформаций, вектор перемещения, конечный элемент, шаг нагружения, физическая нелинейность, напряжение
Благодарности: Исследование выполнено при финансовой поддержке РФФИ и администрации Волгоградской области (проект № 19-41-340002р_а).
Для цитирования: Джабраилов А. Ш., Николаев А. П., Клочков Ю. В., Гуреева Н. А., Ищанов Т. Р. Нелинейное деформирование осесимметрично нагруженной оболочки вращения на основе МКЭ при различных вариантах определяющих уравнений // Известия Саратовского университета. Новая серия. Серия: Математика. Механика. Информатика. 2022. Т. 22, вып. 1. С. 48-61. https://doi.org/10.18500/1816-9791-2022-22-1-48-61
Статья опубликована на условиях лицензии Creative Commons Attribution 4.0 International (CC-BY 4.0)
Article
Nonlinear deformation of axisymmetrically loaded rotation shell based on FEM with different variants of definitional equations
A. Sh. Dzhabrailov10, A. P. Nikolaev1, Yu. V. Klochkov1, N. A. Gureeva2, T. R. Ishchanov1
1Volgograd State Agricultural University, 26 Universitetskiy pr., Volgograd 400002, Russia 2Financial University under the Government of the Russian Federation, 49 Leningradsky Prospekt, Moscow 125993, Russia
Arsen Sh. Dzhabrailov, arsen82@yandex.ru, https://orcid.org/0000-0001-6494-1377 Anatoly P. Nikolaev, anpetr40@yandex.ru, https://orcid.org/0000-0002-7098-5998 Yuri V. Klochkov, Klotchkov@bk.ru, https://orcid.org/0000-0002-1027-1811
Natalia A. Gureeva, Natalya-gureeva@yandex.ru, https://orcid.org/0000-0003-3496-2008 Tlek R. Ishchanov, ishchanov.volgau@yandex.ru, https://orcid.org/0000-0002-2137-4269
Abstract. A curvilinear finite element of the median line of an axisymmetrically loaded shell of revolution with a stiffness matrix of 8x8 size is used when choosing nodal unknowns in the form of displacements and their first derivatives is used. The constitutive equations at the loading step are implemented in two versions. In the first version, the relations of the deformation theory of plasticity are used, which consist of expressions for the elastic and plastic parts. The relationships between strain increments and stress increments were determined by differentiating the equations used. In the second version, the hypothesis of separation of the deformation into elastic and plastic parts was not used. The constitutive equations developed by the authors are obtained on the basis of the hypothesis of the proportionality of the components of the deviators of the stress increments and the components of the deviators of the increments of deformations with the coefficient of proportionality as a function of the chord modulus of the deformation diagram. An example of calculation showing the effectiveness of the developed algorithm is presented.
Keywords: shell of revolution, strain tensor, displacement vector, finite element, loading step, physical nonlinearity, stress
Acknowledgements: This work was supported by the Russian Foundation for Basic Research
and the Administration of the Volgograd Region (project No. 19-41-340002).
For citation: Dzhabrailov A. Sh., Nikolaev A. P., Klochkov Yu. V., Gureeva N. A., Ishchanov T. R.
Nonlinear deformation of axisymmetrically loaded rotation shell based on FEM with different
variants of definitional equations. Izvestiya of Saratov University. Mathematics. Mechanics.
Informatics, 2022, vol. 22, iss. 1, pp. 48-61 (in Russian). https://doi.org/10.18500/1816-9791-
2022-22-1-48-61
This is an open access article distributed under the terms of Creative Commons Attribution 4.0 International License (CC-BY 4.0)
Введение
В настоящее время расчеты тонкостенных конструкций опираются на достаточно развитую теорию оболочек [1-7]. В современных технологических реалиях значительно возрастает необходимость наиболее полного использования прочностных свойств материалов, которые должны обеспечивать устойчивую работу оболочечной конструкции даже при достижении предельных нагрузок. Действие этих нагрузок носит ярко выраженный местный характер, что в свою очередь может привести к возникновению пластических деформаций. Учет физической нелинейности применяемого материала, несомненно, приводит к более экономичным решениями, позволяя эффективно оценить работу конструкции [8-13]. В расчетах за пределами упругости широко используется деформационная теория пластичности с разделением деформации на упругую и пластическую части. Соотношения между приращениями деформаций и приращениями напряжений на шаге нагружения определяются дифференцированием определяющих уравнений деформационной теории пластичности с принятием неизменности объема при пластическом деформировании. При выводе определяющих уравнений в настоящей работе предложена гипотеза о пропорциональности компонент девиаторов приращений деформаций компонентам девиаторов приращений напряжений, на основе которой коэффициент пропорциональности получен в виде функции хордового модуля диаграммы деформирования. В разработанном варианте получения определяющих уравнений на шаге нагружения исключена операция разделения деформации на упругую и пластическую части, что позволило значительно упростить алгоритм формирования определяющих уравнений.
1. Материалы и методы 1.1. Основные соотношения
В процессе деформирования рассматриваются три положения точки срединной линии продольного сечения оболочки вращения: исходное (точка Я0), после j шагов нагружения (точка Н, вектор перемещения v) и после (j + 1) шага (точка Н*, вектор перемещения Av). Соответствующие точки, отстоящие на расстоянии ( от срединной поверхности, обозначаются символами Н0^, Н^ и Н*^.
Положение точки Н0 оболочки вращения описывается радиус-вектором
R0 = xi + rk,
(1)
где х — осевая координата, г — радиус вращения, являющийся функцией от х, 1, к — орты декартовой системы координат.
Вектор локального базиса, касательный к срединной линии продольного сечения оболочки, можно получить дифференцированием (1) по криволинейной координате в:
в? = =(1 + г» ж>в. (2)
Орт нормали к срединной поверхности определяется векторным произведением
e0 = e° х j = -r,xxsi + Соотношения (2) и (3) можно представить в матричном виде
{e0} = [rn°]{i}, {i} =[m0]-1{e0},
2x1 2x2 2x1 2x1 2x2 2x1
(3)
(4)
где |e°}T = |e°e°}; |i}T = {ik}.
Дифференцированием (2) и (3) с учетом (4) можно сформировать матричное соотношение
{еП = [m°J[m°]-1{e°} = [n] {e°}. (5)
21
2x2 2x2 2x1 2x2 2x1
Рассматриваемые положения произвольной точки оболочки вращения при учете гипотезы прямой нормали можно описать следующими радиус-векторами (рис. 1):
R0Í = R0 + С e0, R = R0 + v, R* = R + Av,
Rí
R0Í + v + С(e - e0), R*c = Rc + Av + ((e* - e)
*c — -rC
(6)
где в и в* — орты нормалей к срединной поверхности оболочки в точках Н и Н* соответственно.
Рис. 1. Положение рассматриваемых точек при деформировании оболочки Fig. 1. The position of the points under consideration when the shell is being deformed
Векторы перемещения V и Ду, входящие в (6), можно определить выражениями
V = и?е? + уе0, Ду = Д^?е? + Д^е°. (7)
Для определения базисных векторов в точках , Н^ и Н*^ выполняется дифференцирование (6) по глобальной криволинейной координате в:
е
* = И?, е? = К< = + V, + ((е, - е°),
^ Л Г (8)
е?с = = е? + Ду, + ((е* - еД
Входящие в (8) производные векторов перемещений V и Ду определяются с учетом (5) соотношениями
ув = р?е? + р?е°, = р??е° + рп е°,
(9)
Ду,8 = /?е° + /?е°, Ду>вв = /1? е? + /??е°
Величины р?.../?? в соотношениях (9) являются функциями компонент векторов перемещений V и Ду, а также их производных.
Орты нормалей к срединной поверхности оболочки в точках Н и Н* определяются векторными произведениями
е = е? х j = (К° + V),, х j = (е? + у,3) х j, е* = е* х j = (К + Дv)>s х j = (е? + Дув) х j.
Для нахождения деформаций и их приращений в отстоящем на расстоянии ( от срединной линии слое оболочки вращения при осесимметричном деформировании используются соотношения механики сплошной среды [14]
£
(е?! - е°?) дсс _ (е*? - е?!)
И = о , Д£И = О
2 2 (11) Кк-Б°ск . Л к-Кк
= ____ Д^
Ь22 = тю^ , ДЬ22
И°Ск 22 Кк
где ковариантные компоненты метрических тензоров для исходного и деформированного состояний определяются произведениями соответствующих базисных векторов
е?? = е? • ^ е?? = • , е?? = • е?^. (12)
Выражения (11) с учетом (8)-(10) и (12) можно переписать в виде
£??? = г»? — ку + ((-ку? - к,3у? - ь>3в), £22 = + к?у? + ((-к?у,3 - кк?у?),
Де?? = Д^ - кДу + ((-кДу?>3 - к8Дь? - Д^), (13)
Де£2 = Дук2 + к?Ду? + ((-к?Ду,3 - кк?Ду?),
где к = к? = г,хх,3/г, к2 = х,8/г.
Деформации (13) представляются в более компактной форме следующим образом:
= е + , Д = Д£ар + СД Х«^, (14)
где £а/3 и Ха/з — деформации и искривления срединной поверхности в точке Н, Деа// и ДХа/з — приращения деформаций и искривлений срединной поверхности в точке Н*. Соотношения (14) можно представить в матричном виде
{е<} = [СМ = [ОДМ, {Д^} = [С](Дв} = [С][Ь]{АУ], (15)
2x1 2x44x1 2x44x22x1 2x1 2x4 4x1 2x44x2 2x1
где {е}т = {£11^22X11X22}, {Де}т = {ДецД^Дхп ДХ22}, Мт = {над}, {Д^}т = = {ДмДад}, [Ь] — матрица дифференциальных и алгебраических операторов.
1.2. Зависимости между приращениями деформаций и приращениями напряжений на шаге нагружения
В первом варианте соотношения между компонентами деформаций и напряжений на шаге нагружения устанавливались на основе положений деформационной теории пластичности [15]:
Ьа/3 = °а/3 + °а/3
где ^, 4 — упругие и пластические деформации соответственно.
Согласно [15] компоненты девиаторов упругих и пластических деформаций про порциональны компонентам девиатора напряжений:
4з = & + 4рз, (16)
3 £е 3 £Р
4а - ^ = о - (ааа - <Т8), ^ - ^ = ~ ~ (°аа - <Г8), « = 1, 2, (17)
2 Ог 2 а,-
где ее и е; — средняя упругая и пластическая деформации, ^ и £\ — интенсивности упругих и пластических деформаций соответственно, а3 — среднее напряжение.
В соотношениях (17) полагается, что при пластическом деформировании изменения объема не происходит (£р3 = £11+£|2+£33 = 0).
Таким образом, суммарные деформации на шаге нагружения можно представить соотношениями
4а - ^ = 3 — (^аа - a8), (18)
2 04
где £3 = ее + ^ = ^, £г = £е + .
Связь между средней деформацией и средним напряжением, входящими в (18), такая же, как и в пределах упругости [15]:
{Ч^т)*8 = К1°3) (19)
где V — коэффициент поперечной деформации, Е — модуль упругости материала оболочки.
Приращения деформаций на шаге нагружения определялись дифференцированием (18) и представлялись в общем виде следующим образом:
ДбСаа = ^ Д*11 + ТТ^ Д^22. (20)
о а 11 о а 22
На основе (20) формируется матричное соотношение
{Д^} = [^]{Да}, (21)
2 1 2 2 2 1
£
й
где (Дес}т = {Де^Де^}, {Д^}Т = {До-пДа22}, [Di] — матрица упругопластическо-го деформирования на шаге нагружения.
Во втором варианте определяющие уравнения записаны на основе предложенной гипотезы о пропорциональности компонент девиаторов приращений деформаций компонентам девиаторов приращений напряжений:
Де^ = , (22)
где Де^ = Д&у — Де3, Дв^ = Ди^ — Ди3; Де3, Ди3 — первые инварианты тензоров приращений деформаций и приращений напряжений.
Для определения коэффициента входящего в (22), используются соотношения [16]
ef = У^Де^Де^, af = Дзч, (23)
3Деч ДеИ, аЛ
где £л — интенсивность приращений деформаций, ил — интенсивность приращений напряжений.
Подставляя (22) в (23), можно получить выражения коэффициента у в виде
3 е,л 2 о,
На шаге нагружения принимается равенство
£_Л Де,
(г:
f =1 TT. (24)
f = Д„ • (25)
где Де, — приращение интенсивности деформации, Да, — приращение интенсивности напряжений.
Отношение приращений интенсивностей напряжений и деформаций обозначается хордовым модулем диаграммы деформирования (Еь), а так как на шаге нагружения он неизвестен, то хордовый модуль заменяется касательным модулем диаграммы деформирования (Ек).
При учете (24), (25) коэффициент у определяется выражением
= з Д, = з д = з д (26)
* 2 Да, 2 Ен 2 Ек. ()
На основании эксперимента при растяжении стержня между средними значениями устанавливается зависимость
3 '•{ = (27)
где Е3 — секущий модуль диаграммы деформирования.
На основе (26) и (27) из (22) определяются соотношения между приращениями деформаций и приращениями напряжений. Для осесимметрично нагруженной оболочки они запишутся выражениями
Деспп = —Дааа — Да3( 3— — К2) аа 2Ek * V2 Ек )
или в матричном виде
{д^ } = [адДа}. (28)
2 i 2 2 2 i
£
1.3. Конечный элемент и матрица жесткости на шаге нагружения
Для случая осесимметричного нагружения наиболее целесообразным является использование одномерного конечного элемента с узлами г и ] [17].
Компоненты вектора перемещения внутренней точки конечного элемента в исходном и деформированном состояниях определяются через узловые неизвестные следующим образом:
М =ти1у}, {Д^} = [А]{Аи1у} {и1у} = [ТВД}, (29)
2x1 2x8 8x1 2x1 2x8 8x1 8x1 8x8 8x1
1 т т
где матрица [А] содержит полиномы Эрмита третьей степени, {иу} , {иа} — строки
8x1 8x1
узловых неизвестных конечного элемента в локальной и глобальной системах координат, [Т] — матрица преобразования узловых неизвестных из локальной системы
8x8
координат в глобальную.
Матрица-столбец приращений деформаций на шаге нагружения при учете (15) и (29) может быть представлена в виде
{Д4 } = [С][Ь][А]{Диу} = [С][В ]{Д^У}. (30)
2x1 2x44x22x8 8x1 2x44x8 8x1
Для формирования матрицы жесткости конечного элемента используется функционал, определяющий равенство возможных работ внешних и внутренних сил на шаге нагружения:
Ф = / {Дес}Т ({а} + {ДаЛ ¿V - I {Д^}т ({Р} + {ДРЛ ¿Р, (31)
,/ 2x1 4 2x1 2x1 / ,/ 2x1 4 2x1 2x1 /
V Р
где V — объем элемента, Р — площадь приложения нагрузки, {Р}т = {РиР22}, {ДР}т = {ДРИДР22} — внешние нагрузки после ] шагов нагружения и ее приращения на (] + 1)-м шаге нагружения.
Подставляя в функционал (31) выражения (21), (28)-(30) и выполняя его минимизацию по глобальным узловым неизвестным [17], можно получить матрицу жесткости конечного элемента [К] и вектор его узловых усилий [/], связанных соотношением
у
8x8 8^ 8x1 8x1
где
[К]{Щ} = {/} - {Щ, (32)
[К] = [Т]т / [В]т[С]тттв^у[Т], {/} = [Т]т / [А]т{ДР
8x8 8x8 7 8x4 4x2 2x2 2x44x8 8x8 8x1 8x8 7 8x2 2x1 V Р
{Р} = [Т]т / [А]т{Р} (IР - [Т]т / [В]т [С] {Да} ¿V — невязка Ньютона - Рафсона на
8x1 8x8 Р 8x2 2x1 8x8 V 8x4 4x2 2x1
рассматриваемом шаге нагружения.
Интегралы, входящие в (32), определялись численно на основе формул Гаусса.
2. Результаты
В качестве примера была решена тестовая задача по определению напряженного состояния оболочки вращения в форме усеченного эллипсоида, изображенного на рис. 2.
Оболочка, нижний край которой жестко защемлен (т = 0, т3 = 0), загружена внутренним давлением интенсивности д. Были приняты следующие исходные данные:
внутреннее давление q = 5 МПа; коэффициент Пуассона v = 0.32; модуль упругости Е = 7 ■ 105 МПа; толщина оболочки t = 0.02 м; параметры эллипса: a = 1.3 м, b = 0.9 м; высота усеченной оболочки h = 1.2 м; интенсивность напряжения, соответствующая пределу текучести, &ít = 200 МПа, интенсивность деформации, соответствующая пределу текучести, егт = 0.0023. По условию сходимости вычислительного процесса на первом шаге нагружения при различных вариантах дискретизации оболочки было принято достаточным разбиение рассматриваемой конструкции на 32 элемента.
Было реализовано два варианта расчета. В первом варианте получены конечно-элементные решения на основании определяющих уравнений деформационной теории пластичности (21). Во втором варианте найдены конечно-элементные решения при использовании определяющих уравнений (28), полученных на основе предложенной гипотезы о пропорциональности компонент девиаторов приращений деформаций и девиаторов приращений напряжений.
Для представления диаграммы деформирования на участке упрочнения использовалась параболическая функция вида
а, = Ае2г + Вег + С, (33)
где А = -23461.55 МПа, В = 181201.17 МПа, С = 1574.3 МПа.
Результаты конечно-элементных решений в зависимости от числа шагов нагружения представлены в табл. 1, 2. В них приведены значения меридиональных напряже-
Таблица 1 / Table 1 Меридиональные напряжения при использовании предположения о неизменности объема в результате
пластических деформаций Meridional stresses using the assumption of volume invariance as a result of plastic deformations
b
Рис. 2. Усеченная оболочка вращения Fig. 2. Truncated shell of revolution
Координаты Напряжение, МПа
Coordinates Stress, MPa
Число шагов нагружения, п
X, см / cm (, см / cm Number of loading steps, n
30 50 70
-1.0 274.21 274.50 277.50
-0.5 235.20 229.50 238.20
0.0 0.0 119.00 121.40 119.50
0.5 -28.30 -28.30 -29.10
1.0 -177.50 -176.80 -178.30
-1.0 4.80 5.17 5.13
-0.5 3.70 2.65 2.68
120.0 0.0 0.15 0.19 0.17
0.5 -2.40 -2.36 -2.30
1.0 -4.80 -4.78 -4.79
Таблица 2 / Table 2
Меридиональные напряжения при использовании
второго варианта определяющих уравнений Meridional stresses when using the second version of the constitutive equations
Координаты Напряжение, МПа
Coordinates Stress, MPa
Число шагов нагружения, п
x, см / cm (, см/cm Number of loading steps, n
30 50 70
-1.0 273.05 275.50 276.10
-0.5 226.30 227.60 230.50
0.0 0.0 119.60 120.40 122.60
0.5 -25.20 -26.30 -29.50
1.0 -176.60 -177.80 -178.40
-1.0 4.90 5.12 5.12
-0.5 3.62 2.72 2.65
120.0 0.0 0.15 0.17 0.18
0.5 -2.30 -2.32 -2.30
1.0 -4.74 -4.75 -4.89
ний в зависимости от координаты ( в характерных точках конструкции: в опорном сечении и верхней границе.
Как видно из данных табл. 1, 2, картина напряженного состояния оболочки при нелинейном упрочнении диаграммы деформирования в обоих вариантах определяющих уравнений изменяется несущественно (отклонения в пределах 1%).
Для наглядности и удобства анализа полученных результатов при числе шагов п = 70 в опорном сечении оболочки были построены эпюры напряжений (рис. 3),
300 277.5 266.2
238.2
к
X
<и «
С й К
200 100 0 -100 -200
119.5
1.0
-29.1
-108.5
Толщина оболочки, см -178.3
а / a
300 276.1 265.2
230.5
Е| 200
g 100
X
CD
0 0
-200
-1.0
109.5
Толщина оболочки, см б / b
-178.4
Рис. 3. Эпюры меридиональных напряжений в опорном сечении: а — классический вариант
определяющих уравнений, б — предложенный вариант определяющих уравнений Fig. 3. Diagrams of meridional stresses in the reference section: a — the classical version of constitutive equations, b — the proposed version of the constitutive equations
а также зависимости напряжений вдоль координаты й на внутренней стороне оболочки (рис. 4).
á 300
щ 250
к
1 200
w а
Ü 150 и
2 100 к
й и о к
50
277.3
45 60 75 90 105 120 135 150 Координата S, см
а / a
сз 200
§§ 180
О 160
! 140
¡ 120106.
& 100
g 80
g 60 и
& 40
g 20
о
W 0
175.8
171.4
45 60 75 90 105 120 135 150 Кордината S, см
б / b
Рис. 4. Графики изменения меридиональных (а) и кольцевых (б) напряжений на внутренней
стороне оболочки в зависимости от координаты s Fig. 4. Graphs of changes in the meridional (a) and hoop (b) stresses on the inner side of the
shell depending on the coordinate s
0
0
3. Обсуждения
Для верификации полученных результатов производился анализ сходимости вычислительного процесса. При различных числах шагов нагружения сравнивались полученные значения конечно-элементных решений в рассматриваемых точках конструкции.
Как видно из данных табл. 1, 2, сходимость численных значений меридиональных напряжений в опорном сечении в зависимости от числа шагов нагружения выполняется. Значения контролируемых параметров напряженного состояния несущественно отличаются друг от друга (в пределах 1%).
Для рассматриваемой конструкции проверялось выполнение уравнения статики: сумма проекций всех сил на ось Ох равна нулю
Рг - = 0, (34)
где Рх — проекция равнодействующей всех заданных сил на ось Ох, — проекция равнодействующей всех внутренних усилий на ось Ох.
Численное значение Рх определялось по формуле
Рг = д(^2 - ) = 50(3.14(89)2 - 3.14(33.61)2) = 10667.13 кН. (35)
Чтобы определить внутреннее усилие в опорном сечении, необходимо разбить эпюру напряжений на элементарные геометрические фигуры и найти их площади, затем умножить полученные усилия на соответствующие данному слою длины окружностей. Суммированием полученных растягивающих и сжимающих усилий определяется значение силы в рассматриваемом сечении.
Суммарные внутренние усилия получились равными Pv = 10843.9 кН и Pv = = 10756.9 кН соответственно при использовании определяющих соотношений, полученных на основании (21) и на основании (28). Погрешности вычислений составили 8 = 1.64% и 8 = 0.84% соответственно.
Третьим критерием верификации вычислительного процесса был контроль значения меридионального напряжения на срединной линии поперечного сечения в верхнем крае оболочечной конструкции. В соответствии с физическим смыслом задачи значение напряжения должно стремиться к нулю, так как данный край не загружен.
Анализируя результаты из табл. 1, 2, можно с уверенностью сказать, что данный критерий выполняется.
Заключение
Исходя из вышеизложенного материала, можно заключить, что алгоритм получения определяющих уравнений осесимметрично нагруженных оболочек вращения при учете физической нелинейности применяемого материала на основе предложенной гипотезы о пропорциональности компонент девиатора приращений деформаций и компонент девиатора приращений напряжений является корректным и позволяет получать конечно-элементные решения достаточной степени точности. Использование предположения о неизменности объема в результате пластических деформаций является, по мнению авторов, не совсем корректным. Как показал сравнительный анализ, выполненный для различных вариантов формирования матрицы пластичности на шаге нагружения, связь между средним линейным напряжением и средней линейной деформацией можно осуществлять в виде (27). Это, по мнению авторов, более соответствует физическому смыслу процесса деформирования и значительно упрощает процесс формирования матрицы пластичности на шаге нагружения.
Список литературы
1. Амосов А. А. Техническая теория тонких упругих оболочек. Москва : АСВ, 2011. 304 с.
2. Петров В. В. Нелинейная инкрементальная строительная механика. Москва : Инфа-Инженерия, 2014. 480 с.
3. Cohen H, De Silva C. N. Nonlinear theory of elastic surfaces // Journal of Mathematical Physics. 1966. Vol. 7, iss. 2. P. 246-253. https://doi.org/10.1063/1.1704926
4. Kirillova I. V., Kossovich L. Y. Elliptic boundary layer in shells of revolution under normal edge shock loading // Multiscale Solid Mechanics / ed. by H. Altenbach, V. A. Eremeyev, L. A. Igumnov. Advanced Structured Materials. Vol. 141. Cham : Springer, 2021. P. 249-260. https://doi.org/10.1007/978-3-030-54928-2_19
5. Кабриц C. A., Михайловский Е. И., Товстик П. Е., Черных К. Ф., Шамина В. А. Общая нелинейная теория упругих оболочек. Санкт-Петербург : Изд-во С.-Петербургского университета, 2002. 388 с.
6. Kayumov R. A. Postbuckling behavior of compressed rods in an elastic medium // Mechanics of Solids. 2017. Vol. 52, iss. 5. P. 575-580. https://doi.org/10.3103/S0025654417050120
7. Badriev I. B., Paimushin V. N. Refined models of contact interaction of a thin plate with postioned on both sides deformable foundations // Lobachevskii Jurnal of Mathematics. 2017. Vol. 38, iss. 5. P. 779-793. https://doi.org/10.1134/S1995080217050055
8. Beirao da Veiga L, Lovadina C., Mora D. A virtual element method for elastic and inelastic problems on polytope meshes // Computer Methods in Applied Mechanics and Engineering. 2015. Vol. 295. P. 327-346. https://doi.org/10.1016/j-.cma.2015.07.013
9. Aldakheel F., Hudobivnik B., Wriggers P. Virtual element formulation for phase-field modeling of ductile fracture // International Journal for Multiscale Computational Engineering. 2019. Vol. 17, iss. 2. P. 181-200. https://doi.org/10.1615/IntJMultCompEng. 2018026804
10. Magisano D., Leonetti L., Garcea G. Koiter asymptotic analysis of multilayered composite structures using mixed solid-shell finite elements // Composite Structures. 2016. Vol. 154. P. 296-308. https://doi.org/10.1016/j-.compstruct.2016.07046
11. Ломакин Е. В., Минаев Н. Г. Осесимметричное поле напряжений вблизи кругового выреза в теле с зависящими от вида напряженного состояния пластическими свойствами // Известия Саратовского университета. Новая серия. Серия : Математика. Механика. Информатика. 2019. Т. 19, вып. 3. С. 317-325. https://doi.org/10.18500/ 1816-9791-2019-19-3-317-325
12. Karpov V. V., Ignatev O. V., Semenov A. A. The stress-strain state of ribbed shell structures // Magazine of Civil Engineering. 2017. Iss. 6 (74). P. 147-160. https://doi.org/ 10.18720/MCE.74.12
13. Dzhabrailov A. Sh., Klochkov Yu. V., Marchenko S. S., Nikolaev A. P. The finite element approximation of vector fields in curvilinear coordinates // Russian Aeronautics. 2007. Vol. 50, № 2. P. 115-120. https://doi.org/10.3103/S1068799807020018
14. Седов Л. И. Механика сплошной среды : в 2 т. Т. 1. Москва : Наука, 1976. 492 с.
15. Малинин Н. Н. Прикладная теория пластичности и ползучести. Москва : Машиностроение, 1975. 400 с.
16. Ильюшин А. А. Пластичность. Упруго-пластические деформации. Санкт-Петербург : Ленанд, 2018. 352 с.
17. Клочков Ю. В., Николаев А. П., Джабраилов А. Ш. Конечно-элементный анализ осесимметрично нагруженных оболочек вращения с ветвящимся меридианом при упруго-пластическом деформировании // Строительная механика инженерных конструкций и сооружений. 2013. № 3. С. 50-56.
References
1. Amosov A. A. Tekhnicheskaia teoriia tonkikh uprugikh obolochek [Technical Theory of Thin Elastic Shells]. Moscow, ASV, 2011. 304 p. (in Russian).
2. Petrov V. V. Nelineinaia inkremental'naia stroitel'naia mekhanika [Nonlinear Incremental Building Mechanics]. Moscow, Infa-Inzhenerija, 2014. 480 p. (in Russian).
3. Cohen H., De Silva C. N. Nonlinear theory of elastic surfaces. Journal of Mathematical Physics, 1966, vol. 7, iss. 2, pp. 246-253. https://doi.org/10.1063/U704926
4. Kirillova I. V., Kossovich L. Y. Elliptic boundary layer in shells of revolution under normal edge shock loading. In: H. Altenbach, V. A. Eremeyev, L. A. Igumnov, eds. Multiscale Solid Mechanics. Advanced Structured Materials, vol. 141. Springer, Cham, 2021, pp. 249-260. https://doi.org/10.1007/978-3-030-54928-2_19
5. Kabrits S. A., Mikhailovsky E. I., Tovstik P. E., Chernykh K. F., Shamina V. A. Obshchaia nelineinaia teoriia uprugikh obolochek [General Nonlinear Theory of Elastic Shells]. St. Petersburg, Izd-vo S.-Peterburgskogo universiteta, 2002. 388 p.
6. Kayumov R. A. Postbuckling behavior of compressed rods in an elastic medium. Mechanics of Solids, 2017, vol. 52, iss. 5, pp. 575-580. https://doi.org/10.3103/S0025654417050120
7. Badriev I. B., Paimushin V. N. Refined models of contact interaction of a thin plate with postioned on both sides deformable foundations. Lobachevskii Jurnal of Mathematics, 2017, vol. 38, iss. 5, pp. 779-793. https://doi.org/10.1134/S1995080217050055
8. Beirao da Veiga L., Lovadina C., Mora D. A virtual element method for elastic and inelastic problems on polytope meshes. Computer Methods in Applied Mechanics and Engineering, 2015, vol. 295, pp. 327-346. https://doi.org/10.1016/j-.cma.2015.07.013
9. Aldakheel F., Hudobivnik B., Wriggers P. Virtual element formulation for phase-field modeling of ductile fracture. International Journal for Multiscale Computational Engineering, 2019, vol. 17, iss. 2, pp. 181-200. https://doi.org/10.1615/IntJMultCompEng. 2018026804
10. Magisano D., Leonetti L., Garcea G. Koiter asymptotic analysis of multilayered composite structures using mixed solid-shell finite elements. Composite Structures, 2016, vol. 154, pp. 296-308. https://doi.org/10.1016/j-.compstruct.2016.07.046
11. Lomakin E. V., Minaev N. G. Axisymmetric stress field near a circular cut in a solid with stress state dependent plastic properties. Izvestiya of Saratov University. Mathematics. Mechanics. Informatics, 2019, vol. 19, iss. 3, pp. 317-325 (in Russian). https://doi.org/10. 18500/1816-9791-2019-19-3-317-325
12. Karpov V. V., Ignatev O. V., Semenov A. A. The stress-strain state of ribbed shell structures. Magazine of Civil Engineering, 2017, iss. 6 (74), pp. 147-160. https://doi.org/10.18720/ MCE.74.12
13. Dzhabrailov A. Sh., Klochkov Yu. V., Marchenko S. S., Nikolaev A. P. The finite element approximation of vector fields in curvilinear coordinates. Russian Aeronautics, 2007, vol. 50, no. 2, pp. 115-120. https://doi.org/10.3103/S1068799807020018
14. Sedov L. I. Mekhanika sploshnoi sredy [Continuum Mechanics: in 2 vols.]. Vol. 1. Moscow, Nauka, 1976. 492 p. (in Russian).
15. Malinin N. N. Prikladnaia teoriia plastichnosti i polzuchesti [Applied Theory of Plasticity and Creep]. Moscow, Mashinostroenie, 1975. 400 p. (in Russian).
16. Ilyushin A. A. Plastichnost'. Uprugo-plasticheskie deformatsii [Plastic. Elastic-plastic Deformation]. St. Petersburg, Lenand, 2018. 352 p. (in Russian).
17. Klochkov Yu. V., Nikolaev A. P., Dzhabrailov A. Sh. A finite element analysis of axisymmetric loaded shells of revolution with a branching meridian under elastic-plastic deforming. Structural Mechanics of Engineering Constructions and Buildings, 2013, no. 3, pp. 50-56 (in Russian).
Поступила в редакцию / Received 15.02.2021
Принята к публикации / Accepted 19.07.2021
Опубликована / Published 31.03.2022