CC BY
38
УДК 519.1:539.3
ОСЕСИММЕТРИЧНАЯ ГРАФОВАЯ МОДЕЛЬ УПРУГОГО ТЕЛА С ПЕРЕМЕННЫМ МОДУЛЕМ УПРУГОСТИ
А. А. Тырымов
Волгоградский государственный технический университет,
400005, Россия, Волгоград, пр. В. И. Ленина, 28.
E-mail: [email protected]
Предлагается численный метод анализа упругой среды на основе дискретной модели в виде ориентированного графа. В процессе анализа на основе графового подхода тело рассекается на элементы и для каждого из них строится элементарная ячейка (подграф), являющаяся его моделью. С использованием матриц, представляющих структурные элементы графа, а также уравнений, описывающих разрезанное тело, получены уравнения связного тела. Приведён числовой пример, показывающий 'работоспособность предложенного подхода.
Ключевые слова: математическое моделирование, теория упругости, ориентированный граф, деформация, напряжения.
Введение. Графовая модель однородной упругой среды в осесимметричной постановке рассматривалась в работах Е. Г. Кузовкова и А. А. Тырымо-ва [1,2], где деформации в пределах элемента предполагались постоянными, и в [3], где использовалась линейная аппроксимация деформаций, однако при определении коэффициентов было сделано несколько упрощающих допущений. Автор настоящей работы показывает возможности графового подхода при определении напряжённо-деформированного состояния в упругих телах с непрерывной неоднородностью, причём при выводе матрицы жёсткости графового элемента он не делает никаких упрощающих предположений.
1. Основные уравнения осесимметричной теории упругости. В цилиндрической системе координат г, г, ср закон Гука для изотропного тела запишем, используя тензор плотности напряжений {сг/} , компоненты которого имеют вид а'10 = гагз. В результате получим
а/гг = (Д _|_ 2[л)гєгг + Xrezz +
\
(J — \vSrr (А £(рірі
А .А А -\-
(1)
<7/W — ~єгг + “£zz +
a,rz = (Inrz,
г
з Є‘Р‘Р1
где 7rz = erz + ezr, А и ц — упругие постоянные Ламе, которые в общем случае изотропного неоднородного тела являются функциями координат Г И Z. В дальнейшем будем рассматривать наиболее важный частный случай, когда коэффициент Пуассона постоянен (г/ = const), а модуль упругости изменяется по закону степенной функции Е(г) = Еога, где Eq — постоянный коэффици-
Александр Александрович Тырымов (к.ф.-м.н., доц.), доцент, каф. прикладной математики.
ент, а — вещественное число, г — расстояние от оси симметрии тела. Тогда ц{г) = 0,БЕога(1 + г/)-1, \(r) = коц(г), к0 = 2v(l-2v)~l. (2)
Уравнения равновесия с использованием тензорных плотностей имеют вид [4] г)п1тг f)nlrz f)nlzz f)nlzr
— + — - ro'w = 0, — + — = 0. (з)
О 1 О / О 1 О ' '
Or Oz Oz Or
2. Конструирование элементарной ячейки. Основой для построения моделей в виде графа служит объективный характер операций измерения выбранных независимых переменных. Как известно, целью любого измерения является установление зависимости между значениями переменной и показаниями прибора. В то же время измерение связано как с точками системы, между которыми оно осуществляется, так и с ориентацией измерителя. Эти свойства операции измерения можно представить направленным отрезком, то есть дугой графа.
С точки зрения операции измерения все переменные независимо от их физической природы могут быть подразделены на два типа — параллельные и последовательные [5]. Параллельные переменные характерны тем, что измеритель подключается к двум точкам системы параллельно потоку энергии. Для измерения последовательных переменных прибор должен быть установлен в разрыв энергетического потока. К параллельным переменным относятся линейные и угловые разности перемещений, относительные скорости и ускорения, деформации, разности давлений, температур, электрических потенциалов и так далее. К последовательным переменным относятся силы, моменты сил, электрический ток, тепловой поток, поток жидкости или газа и другие. Выбор связной пары переменных определяется тем, что произведение последовательной и параллельной переменных должно давать скаляр с размерностью мощности или работы. При построении графа исследуемого тела за исходные удобно взять те переменные, которые измеряются установкой прибора непосредственно на исследуемый объект, то есть параллельные переменные. Способ конструирования графа тела связан с процессом измерения полного и независимого комплекта переменных, которые однозначно характеризуют деформированное состояние элементов, полученных в результате декомпозиции. Элементарной ячейкой (графовым элементом) будем называть подграф, соответствующий одному элементу, полученному при разбиении исходной области на отдельные части. При построении конфигурации элементарной ячейки в осесимметричной постановке сечение тела плоскостью, проходящей через его ось, покрываем сеткой координатных линий г = const, z = const, между узлами которых устанавливаем гипотетические измерители. Полученную сеть измерителей преобразовываем в граф, представляя каждый из приборов дугой, а узлы, между которыми осуществлялось измерение, — вершинами графа. При этом вместо собственно измерителя будем ставить символ, отражающий природу измеряемой величины. Если вершины элемента обозначить А, В, С и D, то измерители необходимо установить так, чтобы измерять относительные перемещения точек А, В, С, D друг относительно друга (деформации 5rr, 5rzi 8ZZ, 5zr) и относительно глобальной системы координат (перемещения u^, uf, , и®). Поскольку одни
и те же точки участвуют в разных группах измерений относительно осей г и г, они представлены разными вершинами. В соответствии с описанной последовательностью получаем элементарную ячейку прямоугольного элемента в осесимметричной постановке (рис. 1). Дуги элементарной ячейки вводят два вектора: вектор деформаций {£}с и вектор внутренних сил {/}с:
ные перемещения, ииуелиьленные пивиритим етирин элемента, иг , иг , иг , иР — перемещения точек, соответствующих вершинам в глобальной системе координат; /^г, ГУГ, /^ — нормальные внутренние силы (обобщённые напряжения); /|г, /р, /„г, /^ — тангенциальные внутренние силы; /д, /£, /£, /£> — внутренние силы в глобальной системе координат. Индексы и, (1,1, г обозначают соответственно верхние, нижние, левые и правые дуги элементарной ячейки, а символ Т — операцию транспонирования.
Рис. 1. Элементарная ячейка в осесимметричной постановке: а) компонента 11г, б) компонента IIг; стрелками показано направление обхода контура; в кружках — верхний и нижний индексы переменных, представленных дугами ячейки
3. Построение матрицы жёсткости элементарной ячейки. Описание элементарной ячейки, т. е. связь векторов {/}с и {5}с, а также зависимость их от напряжений {<т} и деформаций {е} моделируемого тела устанавливаем, принимая в качестве инварианта перехода к модели энергию деформации элемента среды:
где {£}с и {/}с определяются формулами (4). В результате энергию элементарной ячейки можно представить, с одной стороны, в виде
а
б
(5)
с другой стороны, энергия деформации элемента сплошной среды есть
[ {(т'}т{е}(1У = [ (а,ггегг + а'ггехх + а'^е^ + а,ггЪх) (IV. (7)
При определении энергии деформации (7) требуется интегрирование выражения а,ггегг + (т1ххехх + (т'^е^ + (т'гх^ГХ1 содержащего неизвестные деформации и напряжения. В связи с этим аппроксимируем неизвестные деформации в пределах элемента полиномами, линейными относительно переменных г иг:
диг
£гг = “Л— = Оо + Я-іТ + (^2Z, дг
дих
Єгх = т— = Оо + 01Г + Ь2г, дг
і І
= Со + Сі г + с2,г, —1 = с?о + (кг + с?2,г, (8)
аг аг
ди ди
Ъг = + -7^7 = Со + (І0 + (сі + <іі)г + (с2 + С12)г, = гиг.
Выбирая систему координат так, чтобы ось г совпадала с осью симметрии тела, а ось г проходила через центр элемента, в пределах элемента имеем Ні ^ г ^ Ке, —Аг ^ 2г ^ Аг. Неизвестные коэффициенты полиномов в (8) выразим через деформации сторон элемента. С помощью интегрирования вдоль соответствующей грани элемента имеем
/ ) іі ^ г Ах/2 ‘
5„ = / ~^~йг = (ао + а\Кс + а2г)Аг, 5ГХ = / -тт-<1г = (со + сіг)Аг,
і Ні
дг
=
-Ах/2
Ке ди
дг
!-Ах/2
~—(Іг = (Ьо + Ьіг)Аг, Ьхг = / —— <1г = ((10 + й\Кс + й2-г)Дг,
дг
/Ні
дг
где Дг = Ее - Дг, ДС = (Дг + Де)/2. В результате получим
5й + 5Л ао + (цПс= " "
(І2
5^г - Щ:г
АгАг
Со + сіДс =
Заметим, что
2 Дг
6а - 5й
51 + 5Г 6о + ЬіДс= 22
АгАг ’ 5Гх ^Г-г 2Аг ’
сі
АгАг ’
2Дг
5Г -51
АгАг
(9)
5й +5Л
<*> + * =
(10)
а2 — С1, 61 — (12.
Это следует из контурных законов [6], применённых к иг и 11х компонентам элементарной ячейки 51гх + 6^г — Ьггх — Ь^г = 0, Ь1ХХ + 5ХГ — 5ГХХ — 5ХГ = 0 и соответствующих выражений из (9). Такой же результат следует из равенств смешанных производных д(диг/дг)/дг = д(диг/дг)/дг и д{дих/дг)/дг = = д{дих/дг)/дг.
В (9) найдена только часть неизвестных коэффициентов из (8). Чтобы найти оставшиеся, воспользуемся (1), (8) и, обеспечивая повышенную точность метода, удовлетворим в пределах элемента уравнениям равновесия (3). В результате получим
((к0 + 2)ц(г) - к0гц'(г))иг =
= г2ц'(г) ((ко + 2)(а0 + (1\Г + a2z) + /с0(Ь0 + hr + b2z)) +
+ (ко + 2 )гц(г)(ао + 2air + a2z) + r2/j,(r)(kobi + с2 + d2), (11)
r/j,(r)(k0ao + к0Ъ2 + 2b2) + ko/J.(r)(c0 + cir + c2z)+
+ r//(r) (c0 + d0 + (ci + di)r + (c2 + d2)z) +
+ //(r) (c0 + d0 + 2(ci + di)r + (c2 + d2)z) = 0. (12)
Из уравнения (11) для функции ц(г), представленной в (2), имеем
ur = (аг((к0 + 2 )(а0 + air + a2z) + fc0(6o + hr + Ь2г)) +
+ (fc0 + 2)(a0 + 2air + a2z)r + (k0h + c2 + d2)r2^(k0 + 2 - fc0a:)_1. (13)
Кроме этого, из (8) получим
ur = (а0 + 0,5a\r + a2z)r + f(z), (14)
где f(z)—произвольная функция. Приравнивая (13) и (14), получим
f(z) = r^2a(ko + 1)ао + 0,5(Зо;А:о + Зко + 4а + Q)a\r + 2a(fco + 1 )a2z+
+ afco(bo + hr + 62<г) + (к0Ъ\ + с2 + d2)r^(k0 + 2 - fcoa:)-1, а это возможно только при выполнении условий
2 (ко + l)a0 + kobo = 0,
2(к0 + 1)а2 + к0Ъ2 = 0, (15)
(3fco + 3fcoQ; + 4а + 6)ai + 2fco(a: + l)&i + 2(с2 + d2) = 0, (16)
Таким образом, f(z) = 0, а
ur = aor + 0,5air2 + a2rz. (17)
Из уравнения (12) следует
(к0с2 + с2 + d2)ii(r) + (с2 + d2)r/j,'(r) = 0,
(к0а0 + к0Ъ2 + 2 b2)r/j,(r) + fc0(c0 + cir)//(r) +
+ (с0 + d0 + (ci + di)r)rn'(r) + (c0 + do + 2(ci + di)r)n(r) = 0.
Отсюда, используя вид функции ц(г) из (2), получаем
кос2 + (а + 1)(с2 + d2) = 0, (18)
kod2 + (fco + 2)b2 + fc0ci + (ск + 2)(ci + di) = 0, (19)
к0с0 + (а + 1)(со + с?0) = 0.
В результате из (15) и (18) имеем Ь2 = к2а2, с2 = кі(і2, к2 = — 2(/с0+1), /с4 = -(а+1)(А:о+о;+1)_1. (20)
Из (16) с учётом (10) и (20) получим я. і = к\Ъ\,
= —2ко((си + I)2 + (ск + 1)Ао + і) х (21)
х ((ко о: -Ь 1) (З/со 3/ьо4а: -Ь 6))
Коэффициент с?і найдём из (19) с учётом (20) и равенства (10):
к = к3с\, к3 = (4к0 + 4 - к0а) [к0(а + 2)]-1. (22)
Используя (17), найдём перемещения вершин элемента АВСИ:
и^ = а0Пг + 0,Ба\К2 — 0,5а2КіАг, и^ = аоИе + 0,5аіД2 — 0,5а2КеАг,
и^ = аоПі + 0,5аіД2 + 0,Ба2ЯіАг, и^ = аоЯе + 0,5аіД2 + 0,Ба2ЯеАг.
С учётом этих зависимостей перемещение -иг представим в виде
иг
+ и^ + Пг +
Вводя в (8) для сокращения записи обозначения а = ао + аі-йс, Ь = Ъо + Ъ\КС1
с = Со + С\їісі (1 = (1о + (кКс имеем
єгг = а + аі(г - Дс) + а2-г, = Ь + &і(г - Дс) + Ъ2г,
Ъх = с + (1 + {сл+ к)(г - Не) + (с2 + <к)г, ^
о і + и1?
гиг = г
'Ч>Ч>
+ \аі{г-нс-ж)+а2~)
где коэффициенты представлены в (9), (20), (21) и (22). В результате деформациям (23) можно придать вид
_ Є + 5Л„ Ы2 - 4 , € ~ Є „
“ 2Аг + 1 АгАг [ с) + АгАг ’
г #гг+Ъг , С ~ 4 , , , Є Є,
“ 2Д,г + АгАг Кс) + АгАг
+ <Г + 5а 5Г - 51
Ъг = \ д г* + *г + (1 + к3) °гI °гг (г - Пс) +
2Аг 2А г АгАг , 24)
ха %и V /
<Г - 5Ъ
+<1 + ^1^дг
р —г —:_____________:_и
6^-Г \ Шс +
1. 5ггг- 51гг ( А г2 \ 5гГ - &гг
+ -к\ л г- Кс - —— + " Л ггг
2 АгАг V 4Кг) АгАг
Подставляя (24) в правую часть равенства (7) и пользуясь формулами (5), (6) и приравнивая выражения при соответствующих компонентах вектора {5}с, получим
я = X {(ш + ^У"++ ш^)л ■
=X ((& - аУ'*' -
&і(г - Дс) /гг ( 1 г - Дс\
_ І / _ '“ХУ /ГГ і / -1- _ -‘•’'С \ _
Д ~ 1у\ АгАг Ч2Дг АгАг/
_1 гЧг-йс-О&Аг^) |
2 АгАг 1 ’
І I к\(т Не) !ГГ / 1 V ПЛ /22
^ = Л ДгД2 ° Чгді+д^іГ +
IV
1 г (г-Д. -0 25Аг я; )а1гЛЛу (25)
2 АгАг у ’
Г г2
Га = Гв = Гс = Гв = ]
/г= I (-х- + х-^'гУгг^>
^ Уу \2Аг АгАг /
/Г = /у(2^-Ж(^-^))^У,
/Г = [
^ Уу \2Аг АгАг У
/Г= / (ттУ + Яс))<т,г^У.
^ .1У\2Аг АгАг V
Выполняя сначала подстановку (24) в формулы (1), а затем подстановку полученных выражений в (25), после интегрирования получим
Ше = Ие{Пе. (26)
Элементы симметрической матрицы [_К"]С получены аналитически в явном виде, но здесь не приводятся из-за их громоздкости.
4. Уравнения графовой модели тела. Граф всего тела строится по тому же принципу, что и графовый элемент, и является результатом операции объединения элементарных ячеек. На рис. 2 представлен граф тела, состоящего из четырёх элементов. Тело, расчлененное на отдельные элементы, и соответствующая ему совокупность элементарных ячеек описываются уравнением
{/} = [КШ- (27)
Здесь {/}т = {/£, /^,..., /Т.} и {6}т = {<^с, 6^,51с) - векторы внутренних сил и деформаций, причём {/}^, {$}7С определяются через (4), [К] = = diag{ [К] 1с, [К] 2с •••, [-К"]пс} — глобальная несвязная матрица жёсткостей тела, расчленённого на отдельные элементы, где [К]гС определены в (26),
Рис. 2. Граф тела из четырёх элементов
а п — число ячеек, образующих граф. Теория графов обеспечивает возможность искусственного конструирования квадратных матриц преобразования, с помощью которых можно исходя из уравнений, описывающих разрезанное на элементы тело, получить уравнение связного тела.
Структурные элементы графа — разрезы, контуры, пути и хорды — могут быть заданы матрицами [6]. Примеры этих матриц для элементарной ячейки приведены в [3]. После преобразования с помощью квадратных матриц дополнений [Дв],д и маршрутов [-Pels’, аналогично [3,7], имеем общесистемный вектор сил F и общесистемный вектор перемещений А:
{F}s = [AB]s{f}, {А}^ = [РсЫ^}, (28)
где индексом S обозначены переменные глобальной системы координат гра-
фа. Поскольку матрицы дополнений и маршрутов не сингулярны, а также в силу равенств [Рс]Т = [Лв]-1 и [Дв]т = [Рс]~1 получим
{/} = №1{fb = [-Pc]J{fb, И = [РоЙ1{АЬ = [Ав]Т{ДЬ- (29)
Тогда, используя (28), (27), (29), получим уравнения, описывающие связное тело:
{F}s = [ВДЛЬ, (30)
где [K]s = [Ab}s[K}[Ab}]i■
Заметим, что при выводе уравнения (30), связывающего общесистемный вектор перемещений и вектор внешних сил, использовалась информация только о структуре графа и физическом законе среды. Кроме этого, используемые в (29) квадратные матрицы преобразования получены не путём их численного обращения, а путём сочленения квадратных матриц, представляющих собой структурные элементы графа. Важно отметить, что уравнения равновесия и совместности деформаций на графовой модели обеспечиваются автоматически как следствие присущих графу фундаментальных зависимостей —
вершинному и контурному законам [6]. Эти законы были установлены Кирхгофом для электрических цепей применительно к токам и напряжениям. Впоследствии оказалось, что они носят фундаментальный характер и применимы ко всем переменным вне зависимости от их физической природы [5]. Согласно вершинному закону алгебраическая сумма последовательных переменных, представленных дугами графа, равна нулю, а в соответствии с контурным законом алгебраическая сумма параллельных переменных равна нулю на любом замкнутом контуре.
Покажем, что континуальную модель сплошной среды удаётся полностью восстановить из её дискретной модели в виде ориентированного графа.
Вершинный закон, если его применить к центральным вершинам (их номера 9 и 10) графа, изображенного на рис. 2, принимает вид
£ГГ I £ГГ _ £ГГ I £ГГ _ I
' Ы2 •/ иЗ ~г м “Г ]г3
ягг I ягг £г £г £г £г _ п /о 1Л
Н2 ' *14 Зт 2 1В3 1АА ~ и>
£гг I £гг £гг I £гг I I _ п
' Ы2 ^иЗ "Г ]иА М ]г3 Jl2 А- Jц — и.
Разложим значения внутренних сил в ряд Тейлора относительно их значений в смежных элементах, сохранив по два члена разложения:
Яг = ЯГ + (ЭДМ)Лг + ..., /11 = ГАз + (дКУдг)Аг + ...,
Я! = ГА + (дГЛ/дг) Ах + ..., /[* = /Ц + (д/%/дх)Ах + ..., , ,
Я2 = ЯГ + {дГЛ/дг)Аг + ..., /£ = Ги1 + (дГи1/дг)Аг + ..., М
ГА = ГА + (д/%/дх)Ах + ..., Щ + (д^/дх)Ах + ....
Подставив в (31) выражения из (32) с учётом (25), получим
(зз>
В силу произвольности выбора размеров ячейки из (33) следует, что подынтегральные выражения в них должны обращаться в ноль, а это влечёт выполнение уравнений равновесия (3).
С использованием связи между физическими компонентами тензора напряжений и контравариантными компонентами
7*7* 2 ^0^0 V X
О'(гг) ® 1 & (гг) ® 1 ^(фф) V (7 , (7(гх) СГ
и соотношений
а(гг) = о''гг/г, а{гг)=а'гг/г, ам=га^, о[гг)=о'™/г уравнения равновесия приобретают традиционную форму [8]:
да(гг) да(гх) 1 _ до(гх) да^ _
~дГ+ + г( (гг) “ ^ + +
Поступая, как в [3], покажем, что контурный закон графов, применённый к определённой комбинации контуров, приводит к уравнениям совместности деформаций. Рассмотрим разность контуров, образованных дугами ячеек с номерами 1 и 3 компоненты иг графа (ячейки левого ряда на рис. 2, а):
^ _ ^1 + §г1 _ 5из + 8аз + £13 _ §гз = 0_
Контурные законы для дуг ячеек с номерами 1 и 2 компоненты 11х (ячейки верхнего ряда на рис. 2, б) приводят к равенству
?и1 £(11 гИ | гг1 ги2 , ?с12 , г/2 сг2 __ п /ос:\
дгг огг + + дгг + 0. (35)
Дуги 5^1 и 5™г образуют замкнутый контур и поэтому 5^ = 6'^. Разложим
деформации 5гг,$гг и К\^ГЛ в РЯД Тейлора, сохраняя только два члена:
5% = С1 + (д5^/дх)Ах,
6$г = $гг + (д5гг/дх) Ах = 5™г + (д5^/дх) Ах =
= С3 + (^(с1 + {д5^/дх)Ах)/дх) Ах,
5Ц = 51г\ + {д51г\/дг)Аг, (36)
5Ц = 51г\ + (дёЦ/дх) Ах,
5ГГ1 = 61Г3Х + (д613х/дг)Аг = 513х + (д (61г\ + (д61г\/дх) Ах) / дг) Аг.
Подстановка (36) в (34) даёт
— (35^1 / дх) Ах + (0(5^ + (05^ / дх) Ах)/дх)Ах+
+ (д61г\/дг)Аг — (д(61г\ + (д61г\/дх)Ах)/дг)Аг =
= (д25'^/дх2)Ах2 - {с)251г\/{дгдх))АгАх = 0.
Принимая во внимание, что 5ГГ ~ (диг/дг)Аг, 5ГХ ~ (диг/дх)Ах, и сократив на АгАх2, получаем, что компонента графа 17г приводит к уравнению
д2(диг/дг)/дх2 — д2 (диг / дх) / (дгдх) = 0. (37)
После аналогичных преобразований уравнение (35) приобретает вид
д2(дих/дх)/дг2 — д2(дих/дг)/(дгдх) = 0. (38)
Складывая (37) и (38), получим известное уравнение совместности деформаций:
д2(диг/дг)/дх2 + д2(дих/дх)/дг2 — д2(диг/дх + дих/дг)/(дгдх) = 0.
Таким образом, при уменьшении размеров сетки вершинный и контурный законы приводят к выполнению уравнений равновесия и совместности деформаций теории упругости, что влечёт сходимость приближённого решения к точному.
Важным преимуществом графового метода является возможность аппроксимации поля деформаций четырёхузлового элемента линейными полиномами. В результате при таком же уровне аппроксимации графовый метод использует систему уравнений, содержащую примерно в 3 раза меньше уравнений по сравнению с системой, используемой в методе конечных элементов.
Специально был рассмотрен вопрос о влиянии размеров элементов на результаты расчётов. Тестирование показало, что при соотношении сторон прямоугольного элемента, равном 1/100, решение не искажается. Анализ также
показал, что в пределах допустимых соотношений сторон элементов координатные линии можно сгущать произвольно и после широких элементов можно сразу же размещать узкие.
5. Численные результаты. Проведено сравнение результатов, полученных по предложенной методике, с известным аналитическим решением [8] для толстостенной цилиндрической трубы с внутренним Г 1=3 мм и внешним г2 = 7 мм радиусами, находящейся под действием равномерного давления р = 1 МПа. При этом предполагается, что модуль Юнга материала является функцией вида Е = Еога, где Ео = 100 МПа, а коэффициент Пуассона V = = 0,3. В таблице приведены отношения аг/р и а^/р для различных точек в радиальном направлении в зависимости от показателя степени а.
аг/р °v/p
г а 3 4 5 6 7 3 4 5 6 7
0 -1,000 -0,991 -0,464 -0,467 -0,216 -0,217 -0,081 -0,082 0,000 -0,001 1,450 1,443 0,914 0,917 0,666 0,668 0,531 0,532 0,450 0,451
1 -1,000 -0,993 -0,542 -0,545 -0,284 -0,285 -0,117 -0,118 0,000 -0,001 0,904 0,900 0,779 0,781 0,727 0,728 0,706 0,707 0,701 0,701
2 -1,000 -0,996 -0,756 -0,755 -0,542 -0,540 -0,303 -0,300 0,000 -0,002 0,487 0,486 0,610 0,611 0,745 0,746 0,887 0,887 1,033 1,030
5 -1,000 -0,994 -0,613 -0,615 -0,355 -0,356 -0,159 -0,160 0,000 -0,002 -0,139 -0,138 0,114 0,113 0,555 0,561 1,286 1,290 2,427 2,418
Верхняя строчка каждого ряда соответствует аналитическому решению, а нижняя строчка — графовому методу. Видно, что результаты хорошо согласуются. Отметим, что выделяемый в расчётах фрагмент трубы покрывался не очень плотной (16 х 9) сеткой элементарных ячеек.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Кузовков Е. Г., Тырымов А. А. Графовая модель упругой среды в осесимметричной постановке / В сб.: Моделирование в механике: Сб. научных трудов. СО АН СССР: Новосибирск, 1990. С. 103-109. [Kuzovkov Е. G., Tyrymov A. A. Graph model of an elastic solid in the axisymmetric formulation / In: Modelling in Mechanics. SO AN SSSR: Novosibirsk, 1990. Pp. 103-109].
2. Тырымов А. А. Треугольный элемент графовой модели для осесимметричной задачи теории упругости / В сб.: Численные методы решения задач теории упругости и пластичности: Тр. XVIII Межресп. конф. (1-3 июля 2003 г.). Новосибирск: Нонпарель, 2003. С. 187-191. [Tyrymov A. A. Triangular element of the graph model for the axisymmetric problem of elasticity theory / In: Numerical methods for solving problems in the theory of elasticity and plasticity. Novosibirsk: Nonparel, 2003. Pp. 187-191].
3. Kuzovkov E. G. Axisymmetric graph model of an elastic solid // Strength of Materials, 1996. Vol. 28, no. 6. Pp. 470-485.
4. Крон Г. Исследование сложных систем по частям — диакоптика. М.: Наука, 1972. 542 с. [Kron G. Ivestlgatlon of Complex Systems by Parts (Diacoptics). Moscow: Nauka, 1972. 542 pp.]
5. Trent H. M. Isomorphisms between Oriented Linear Graphs and Lumped Physical Systems // J. Acoust. Soc. Am., 1955. Vol. 27, no. 3. Pp. 500-527.
6. Swamy M. N. S., Thulasiraman K. Graphs, Networks and Algorithms. New York: John Wiley, 1981. 592 pp.; русск. пер.: Свами М., Тхуласираман К. Графы, сети и алгоритмы. М.: Мир, 1984. 454 с.
7. Тырымов А. А. Сингулярный элемент графовой модели упругой среды в декартовой системе координат // Вычисл. мех. сплошн. сред, 2011. Т. 4, №4. С. 125-136. [Tyrymov А. А. A singular element of the graph model of an elastic medium in a cartesian coordinate system // Computational Continuum Mechanics, 2011. Vol. 4, no. 4. Pp. 125-136].
8. Демидов С. П. Теория упругости. М.: Высш. шк., 1979. 432 с. [Demidov S. P. Theory of Elasticity. Moscow: Vyssh. Shk., 1979. 432 pp.]
Поступила в редакцию 28/XII/2010; в окончательном варианте — 27/11/2012.
MSC: 74S30; 74В05
NUMERICAL SOLUTION OF AXISYMMETRIC PROBLEM OF THE THEORY OF ELASTICITY ON THE BASIS OF CONTINUUM GRAPH MODEL
A. A. Tyrymov
Volgograd State Technical University,
28, Lenin avenue, Volgograd, 400005, Russia.
E-mail: [email protected]
A numerical method for analysis of the stress strain state of elastic media based on a discrete model in form of directed graph is suggested. To analyze a deformable body using the graph approach, we partition a solid body on elements and replace each element by its model in the form of an elementary cell. The matrices, presenting several structure elements of the graph, and the equations, describing the elementary cells, contribute to deriving the constitutive equations of the intact body. Numerical examples are presented.
Key words: mathematical simulation, elasticity, directed graph, strain, stress.
Original article submitted 28/XII/2010; revision submitted 27/11/2012.
Alexander A. Tyrymov (Ph.D. (Phys. & Math.)), Dept, of Applied Mathematics.
