Full Text

Введение Графы широко используются как структурные модели физических и технических систем, допускающих идеализированное представление в виде схем с сосредоточенными параметрами [1-6]. Сюда относятся электрические цепи, механизмы, стержневые конструкции, теплообменники, акустические системы, разнообразные комбинированные системы, такие как электромеханические, гидромеханические, механоакустические и др. Объясняется это рядом достоинств графовых моделей. Одно из обстоятельств, способствующих применению графов для описания различных технических систем, заключается в их связи с операциями измерения основных переменных, однозначно определяющих состояние этих систем. Физическая природа переменных изменяется от одного класса систем к другому, в то время как техническая суть операции по их измерению остается неизменной. Подобный подход позволяет представить систему ее переменными, изображенными в виде графа, а граф использовать при анализе этой системы. Однако помимо преимуществ, связанных с наглядностью представления структуры и связи основных переменных, графы позволяют целенаправленно в рамках определенного алгоритма провести сложные математические преобразования. Анализ физических и технических систем сводится к получению уравнений состояния системы на основании уравнений состояния ее частей. Теория графов позволяет разработать систематическую методику вывода этих уравнений. Объединение описаний частей в описание связной системы осуществляется с помощью матриц, характеризующих структуру графа. Анализ системы с сосредоточенными параметрами с использованием матричных методов возможен на основе сетевой модели Г. Крона. Подход Крона [7, 8] основан на объединении структурных свойств электрической цепи, представляемых схемой соединения ее компонентов, с функциональными свойствами, представленными уравнениями состояния компонентов. Применяя аналоговое моделирование, Крон предложил использовать электрические цепи для представления самых разнообразных физических, технических и экономических систем. Были разработаны эквивалентные электрические сети, предназначенные для моделирования задач гидродинамики, электрического поля, квантовой механики, линейного программирования [7-10]. При построении электрической сети, эквивалентной упругому телу, Крон принял, что механические напряжения моделируются электрическими токами, а деформации - разностями потенциалов. Однако необходимость учета характеристик реально существующих электрических компонентов и процессов привело к тому, что элементарные ячейки Крона отличаются сложностью конфигурации. Так, в элементарную ячейку для непосредственного физического моделирования плоского напряженного состояния упругого тела входит 12 компонентов со сложным их соединением [11]. Обращая основное внимание в своих работах на универсальность возможностей использования в качестве моделей эквивалентных электрических сетей в различных областях науки и техники, Крон не приводит методик по их конструированию и разработке алгоритмов при определении параметров элементарных ячеек. Подход Крона не получил распространения применительно к задачам механики сплошной среды. Опираясь на идеи Крона, Е.Г. Кузовков [12-14] предложил использовать электрические ячейки только для отображения связей и характеристик моделируемой упругой среды. Это позволило в компонентах ячеек не учитывать такие свойства реально существующих электрических элементов, как потери, утечки, взаимные наводки. Подобный подход дал возможность при описании упругого элемента применить ячейку, состоящую всего из четырех компонентов. Параметры ячейки определялись сравнением выражений, полученных на основе законов Ома и Гука. В работе [15] Е.Г. Кузовков обратил внимание на то, что при моделировании напряженно-деформированного состояния упругих тел можно обойтись без электромеханической аналогии, если эквивалентную электрическую сеть трактовать как ориентированный граф. Способ конструирования графовой модели, ее конфигурация, применение специальным образом сконструированных матриц для вывода определяющей системы уравнений применительно к плоской и осесимметричной задачам теории упругости предложены Е.Г.Кузовковым и подробно изложены в [15-20]. В работе [21] на основе графовой модели построен двумерный сингулярный элемент для расчета напряженно-деформированного состояния в окрестности особых точек разреза. Использованию графового метода для анизотропных и неоднородных сред посвящены работы [22, 23]. Графовая модель упругой среды в полярной системе координат рассматривалась в [24]. В работе [25] в качестве дискретной модели трехмерного упругого тела в декартовой системе координат используется ориентированный граф как идеализация гипотетических приборов, необходимых для измерения деформированного состояния тела. В этой работе приведено описание элементарной ячейки - подграфа, соответствующего одному элементу, полученному при разбиении исходной области на мелкие части. При определении конфигурации элементарной ячейки, так же как предложено в [25], «область, занятую телом, покрываем сетью координатных плоскостей x = const, y = const, z = const, между узлами которых устанавливаем гипотетические измерители, определяющие абсолютные нормальные деформации и деформации ,, обусловленные сдвигом. Полученную сеть измерителей преобразовываем в граф, представляя каждый из приборов дугой, а узлы, между которыми осуществлялось измерение, - вершинами графа. При этом вместо собственно измерителя будем ставить символ, отражающий природу измеряемой величины. В результате получаем элементарную ячейку прямоугольного шестигранного элемента, имеющую 36 дуг и 24 вершины и состоящую из трех компонент (рис. 1): Рис. 1. Элементарная ячейка, состоящая из трех компонент Fig. 1 An elementary cell consisting of three components Вершины с номерами 22, 23, 24 - корни подграфов соответствующих компонент элементарной ячейки. Каждая из дуг графа отображает одну из следующих пар параллельных и последовательных переменных: (,), (,), (,), (,), (,), (,), (,), (,), (,), где , , и , , , , , - нормальные и тангенциальные внутренние силы (обобщенные напряжения). Таким образом, дугам элементарной ячейки соответствует следующая пара векторов: (1) где символ t обозначает операцию транспонирования, а индекс с указывает на принадлежность переменных отдельной элементарной ячейке. Использование в качестве инварианта перехода к модели энергию деформации элемента среды позволило в [25] получить матрицу жесткости для восьмиузлового шестигранного прямоугольного упругого элемента с 24 степенями свободы при линейной аппроксимации деформаций. В данной работе рассматривается метод построения определяющей системы, описывающей связное пространственное тело. Он основан на объединении описаний отдельных элементов, полученных в [25], в описание связного тела, которое осуществляется с помощью матриц, характеризующих структуру графа. Решение полученной системы позволяет определить напряженно-деформированное состояние упругого тела. 1. Матричное представление структурных элементов графа В анализе напряженно-деформированного состояния упругого тела на основе графовой модели используются разнообразные структурные элементы. Поскольку терминология теории графов до конца не стандартизирована, дадим определения основных используемых ниже понятий [5, 26-28]. Маршрут - конечная последовательность необязательно различных дуг, у которых одна из вершин первой дуги является вершиной второй, оставшаяся вершина второй - вершиной третьей и т.д. Путь - незамкнутый ориентированный маршрут, не содержащий повторяющихся дуг. Граф называется связным, если в нем существует путь между каждой парой вершин. Контур графа - замкнутый ориентированный маршрут, не содержащий повторяющихся дуг. Дерево - связанный подграф, не содержащий контуров. Остовное дерево - дерево, содержащее все вершины исходного графа. Дуги дерева называются ветвями. Хорды - дуги графа, не принадлежащие дереву. Компонента графа - максимально связный подграф этого графа. Если исходный граф состоит из нескольких компонент, то деревья компонент образуют лес. Вершина является корнем дерева, если из нее существует ориентированный путь по всем остальным вершинам этого дерева. Базисный контур - контур, образованный хордой и теми ветвями дерева, которые составляют единственный путь между ее концевыми вершинами. Число линейно независимых базисных контуров равно числу хорд графа и определяется формулой где - цикломатическое число; m - число дуг; n - число вершин; k - число компонент графа. При изучении структурных связей, порождаемых графом, удобно использовать матричный аппарат. Рассмотрим матрицы, задающие отношения инцидентности между дугами, с одной стороны, и вершинами, путями, хордами и контурами - с другой [5,6]. Расположение строк и столбцов в этих матрицах определяется номером вершин и обозначением дуг графа. Матрица инцидентности. Матрица инцидентности [A] ориентированного графа, имеющего n вершин и m дуг, является матрицей размером n ´ m с элементами aij (i = 1, 2, …, n; j = 1, 2, …, m), определяемыми следующим образом: Ранг матрицы [A] равен n-k, где k - число компонент графа и совпадает с рангом графа rk. Строки матрицы инцидентности соответствуют вершинам, которые остаются изолированными при удалении дуг, перечисленных в соответствующих строках. Если оставить невырезанными корневые вершины, то получим усеченную матрицу инцидентности. Строки в этой матрице перечислены в порядке возрастания номеров вершин от 1 до 21, столбцы соответствуют дугам графа в порядке xx1, yy1, zz1, xx2, yy2, zz2, xx3, yy3, zz3, xx4, yy4, zz4, xy1, yx1, xz1, zx1, yz1, zy1, xy2, yx2, xz2, zx2, yz2, zy2, xy3, yx3, xz3, zx3, yz3, zy3, xy4, yx4, xz4, zx4, yz4, zy4. Матрица корневых путей. Матрицей корневых путей остовного дерева [P] ориентированного графа называется матрица, элементы которой pij определяются по следующему правилу: Матрица путей строится в соответствии с рис. 2, где изображены деревья элементарной ячейки. Рис. 2. Деревья элементарной ячейки Fig. 2.The trees of the elementary cell from Матрица путей имеет размер rk´m. Строки матрицы путей соответствуют вершинам, из которых начинается путь к корню, столбцы - дугам, образующим путь. За положительное направление пути принимаем направление от вершины к корню. Каждый путь обозначаем индексом вершины, которую данный путь связывает с корнем. Нумерация строк и столбцов точно такая, как и в матрице инцидентности. Матрица хорд. Матрицей хорд [B] называется матрица с элементами bij, где Матрица хорд имеет ранг, равный цикломатическому числу m(G) = m - n + k, ее размер - m(G)´m. В матрице хорд индексация строк и столбцов такая же, как у матрицы базисных контуров. Элементы матрицы [B], стоящей на пересечении строк и столбцов с нумерацией xy1, yx1, xz1, zx1, yz1, zy1, zx2, zy2, xy3, xz3, yz3, yx4, xz4, zx4, yz4, равны единице, а остальные равны нулю. Матрица базисных контуров. Матрицей базисных контуров [C] ориентированного графа называется матрица, элементы cij которой определяются так: Каждый базисный контур содержит только одну хорду и те ветви дерева, которые кратчайшим путем связывают ее вершины. Направление базисного контура задается направлением определяющей хорды. Ранг матрицы базисных контуров определяется линейно независимыми контурами, их число m(G) = m - n + k, т.е. совпадает с цикломатическим числом графа. Размер матрицы контуров [C] равен m(G)´m. Контуры матрицы [C] упорядочены согласно выбору задающих их хорд. Порядок расположения строк соответствует хордам xy1, yx1, xz1, zx1, yz1, zy1, zx2, zy2, xy3, xz3, yz3, yx4, xz4, zx4, yz4], столбцы - дугам графа с прежней индексацией. В трехмерной декартовой системе координат транспонированная матрица [C]t выглядит так: xy1 yx1 xz1 zx1 yz1 zy1 zx2 zy2 xy3 xz3 yz3 yx4 xz4 zx4 yz4 xx1 -1 -1 yy1 -1 -1 zz1 -1 -1 xx2 1 -1 yy2 1 -1 zz2 1 1 xx3 1 1 yy3 1 -1 zz3 1 -1 xx4 1 -1 yy4 1 1 zz4 1 -1 xy1 1 yx1 1 xz1 1 zx1 1 yz1 1 [C]t = zy1 1 xy2 -1 -1 -1 yx2 -1 -1 -1 xz2 -1 -1 -1 zx2 1 yz2 -1 -1 -1 zy2 1 xy3 1 yx3 1 1 -1 xz3 1 zx3 -1 -1 -1 yz3 1 zy3 -1 -1 -1 xy4 1 -1 1 yx4 1 xz4 1 zx4 1 yz4 1 zy4 1 -1 1 В матрице базисных контуров и ниже - в матрице линейных комбинаций базисных контуров - записаны только ненулевые элементы. Заметим, что в третьем, четвертом и пятом столбцах матрицы [C]t число отличных от нуля элементов больше четырех. С увеличением числа элементов, на которые рассечено тело, растет количество дуг графа, и появляются столбцы, содержащие большее число единиц (с тем или иным знаком). Это усложняет процесс программирования и, кроме того, ухудшает точность вычислений. В то же время число отличных от нуля элементов в столбцах матрицы [C]t можно сделать равным четырем при любой степени декомпозиции. Базисные контуры, которые содержат только одну хорду и три ветви, остаются без изменения. В остальных случаях следует перейти от базисных контуров к их линейным комбинациям. Например, из рис. 2, видно, что в базисный контур Cxz1 входят дуги xz1, xx4, xy4, -xz2, -xy2, -xx1, а в контур Cxz4 - дуги xz4, xy4, -xz2, -xy2. Разность этих базисных контуров даст контур , который включает в себя дуги xz1, xx4, -xz4, -xx1. Таким образом, контуры, полученные линейной комбинацией базисных, содержат две хорды и две ветви. Одна хорда имеет положительное направление и принимается в качестве определяющей данный контур, вторая хорда входит в контур с отрицательным знаком. В результате получим транспонированную матрицу линейных комбинаций базисных контуров: xy1 yx1 xz1 zx1 yz1 zy1 zx2 zy2 xy3 xz3 yz3 yx4 xz4 zx4 yz4 xx1 -1 -1 yy1 -1 -1 zz1 -1 -1 xx2 1 -1 yy2 1 -1 zz2 1 1 xx3 1 1 yy3 1 -1 zz3 1 -1 xx4 1 -1 yy4 1 1 zz4 1 -1 xy1 1 yx1 1 xz1 1 zx1 1 = yz1 1 zy1 1 xy2 -1 -1 -1 yx2 -1 -1 xz2 -1 -1 zx2 1 yz2 -1 zy2 1 xy3 1 -1 yx3 1 xz3 1 zx3 -1 -1 yz3 -1 1 zy3 -1 -1 xy4 -1 1 yx4 1 xz4 -1 1 zx4 -1 1 yz4 1 zy4 -1 1 Структурные матрицы для элементарной ячейки в двумерном случае представлены в [20, 23]. 2. Свойства топологических матриц. Построение несингулярных матриц преобразования Между топологическими матрицами существуют взаимные связи, выражающие свойства ортогональности или позволяющие по одной из них определить другую [6]. Известно, что матрицы инцидентности и контуров обладают свойством ортогональности: [A][C]t = 0, [C][A]t = 0. (2) Кроме того, из определения матрицы корневых путей остовного дерева и матрицы хорд следует [B][P]t = 0, [P][B]t = 0. (3) Действительно, в столбцах матрицы [P], соответствующих хордам, стоят нули, а единицы расположены в столбцах, соответствующих ветвям, в то время как в матрице [B]t в строках, соответствующих хордам, стоят единицы, а ветвям - нули. С другой стороны, произведение матриц инцидентности и путей основного дерева, а также хорд и контуров дает единичную матрицу: [A][P]t = [I], [P][A]t = [I]. (4) [B][C]t = [I], [C][B]t = [I]. (5) Ранг каждой из матриц [A] и [P] равен рангу графа rk, а ранг матриц [B] и [C] равен цикломатическому числу m. Это позволяет сконструировать квадратные матрицы размером m´m, так как rk + m = m, а число столбцов у матриц [A], [B], [C], [P] равно m. Для этого следует применить операцию сочленения двух матриц по вертикали [28], которая заключается в простом приписывании второй матрицы к первой снизу. Таким образом, сначала располагаются все строки первой матрицы, а затем второй - без нарушения их порядка в исходных матрицах. В результате можно построить две квадратные матрицы Матрицу [AB] можно назвать матрицей дополнений, так как она построена с помощью матрицы инцидентности [A], содержащей разрезы, дуги которых дополняют множество изолированных вершин до исходного графа, а матрица [B] включает в себя хорды, дополняющие ветви остовного дерева до исходного графа. Матрицу [PC] назовем матрицей маршрутов, поскольку она составлена из матрицы корневых путей [P] и матрицы контуров [C], т.е. с использованием матриц открытых и замкнутых маршрутов. Из правил оперирования с блочными матрицами вытекают соотношения [AB]t = [At, Bt], [PC]t = [Pt, Ct]. Учитывая свойства (2)-(5), получаем Отсюда следует, что матрицы [AB] и [PC]t взаимно обратные, т.е. [PC]t = [AB]-1, [AB]t = [PC]-1. (6) Заметим, что строки и столбцы матриц [AB] и [PC] имеют одинаковые обозначения. Поскольку связь (6) между этими матрицами получена не с помощью математического обращения, а путем искусственного построения, то при выполнении операции сочленения необходимо специально следить за соответствием распределения строк и столбцов в каждой из матриц, обеспечивая одинаковый порядок. 3. Уравнения графовой модели тела Тело, представленное в виде отдельных элементов, и соответствующая ему совокупность элементарных ячеек описывается уравнением (7) Здесь {f} и {d} - векторы внутренних сил и деформаций. причем определяются (1), - глобальная несвязная матрица жесткости тела, представленного в дискретном виде: [K] = diag[K1c, K2c,…,Knc], где - число ячеек, образующих граф, а определяется формулой (29) из работы [25]. Уравнение (7) связывает внутренние силы с деформациями элементов и не позволяет в общем случае определить напряженно-деформированное состояние тела, поскольку, как правило, заданными бывают внешние силы и перемещения. Для получения уравнений связного тела, следуя [18], используем преобразования обобщенных координат, описывающих разрезанное на элементы и связное тело. Вводим преобразования переменных с помощью квадратных матриц [AB] и [PC]: (8) (9) Поскольку матрицы дополнений и маршрутов несингулярны, то с учетом (6) возможны обратные преобразования, которые не требуют процедур численного обращения: (10) Подставив в правую часть (7) значения {δ} из (10), а также умножив это уравнение слева на [AB], получаем уравнение, описывающие связное тело: {F} = [AB][K][AB]t{Δ}, (11) где Следует отметить, что уравнение (11),описывающее напряженно-деформированное состояние дискретной модели, несет информацию о структуре графа и физическом законе Гука. Дифференциальные уравнения равновесия и совместности деформаций при выводе (11) не использовались. В то же время для построенной графовой модели, как показано в [25], уравнения равновесия и совместности деформаций обеспечиваются автоматически как следствие вершинного и контурного законов Кирхгофа [6, 27], которые имеют фундаментальный характер. Рассмотрим общесистемный вектор сил (8). Его член [A]{f} на основании вершинного закона [5,15] представляет собой вектор вершинных переменных {F}e, которые интерпретируются как заданные внешние силы. Член [B]{f} выделяет часть дуг графа {f}b, являющихся хордами, которые интерпретируются как неизвестные внутренние силы. Вектор общесистемных перемещений (9) включает в себя член [P]{δ}, который представляет собой перемещения узловых точек в глобальной системе координат {U}. На основании контурного закона [6,18] член [C]{δ} равен нулю. В результате уравнение (11) принимает вид откуда следует Решение последнего уравнения получаем в виде {U} = [Ks]-1{F}e. По известным перемещениям определяем деформации, а затем с помощью закона Гука находим напряжения, действующие в каждом элементе. 4. Численные результаты В качестве примера, иллюстрирующего предлагаемый метод, приведем задачу из [29]. В ней рассматривается расчет напряженного состояния для пространственной структуры, изображенной на рис. 3. Как указано в [29], эта расчетная схема может быть использована для моделирования шахтной выработки для поддерживания кровли при проведении горных работ. На рис. 3 показана расчетная схема и дискретизация структуры на шестигранные прямоугольные элементы. Изображена восьмая часть расчетной области. Плоскости xy, yz и zx являются плоскостями симметрии. Нормальные перемещения узлов, лежащих в плоскостях симметрии, задавались равными нулю. К нижней поверхности приложена равномерно распределенная нагрузка p = 0,1 МПа, на верхней поверхности исключены вертикальные перемещения. Другие грани конструкции свободны от нагрузки. Размеры, как и в [29], верхнего параллелепипеда 20×20×40 (см), нижнего - 60×60×80 (см). В [29] приведены решения рассматриваемой задачи методом граничных элементов (МГЭ), методом конечных элементов с использованием модифицированных гексаэдров Айронса-Зенкевича, а в работе [30] с помощью шестигранного прямоугольного конечного элемента с применением линейных базисных функций для аппроксимации напряжений. Анализируется концентрация напряжений σу/р, вдоль диагонали плоскости соединения кровли и крепи. В расчетах число граничных сегментов в зависимости от версии программ МГЭ составляло 16 либо 44, в МКЭ сетка разбиения на элементы редкая с использованием 274 степеней свободы для ее узлов [29]. Авторы отмечают, что дискретизация слишком грубая и требует более детального описания. Анализ напряжений методом КЭ в [30] проведен при разбиении нижней части конструкции на 196 элементов, верхней - на 64 элемента. В графовом подходе верхний параллелепипед разбивается на 150 элементов, нижний - 550 элементов. На рис. 4 показаны графики напряжений σу/р, полученные МКЭ, МГЭ и графовым методом. Заметим, что, как и в [29, 30], отсчет расстояния ведется от вершины угла, лежащего на ребре крепи. Рис. 3 Геометрия задачи Fig. 3. Geometry of the problem Рис. 4. Распределение напряжений σу/р, по диагонали сечения крепи Fig. 4. The σу/р stresses distribution along the diagonal cross section of the lining Как видно из графиков, результаты расчетов графовым методом согласуются с общей тенденцией распределения напряжений из [29, 30], позволяя при этом получить более точное решение.

About the authors

A A Tyrymov

Volgograd State Technical University

References

Statistics

Views

Abstract - 265

PDF (Russian) - 130