Full Text

Введение При численном решении разнообразных задач механики сплошной среды основными достаточно универсальными и широко распространенными являются методы конечных и граничных элементов в различных вариантах и сочетаниях с другими методами. Отличительная особенность этих методов состоит в том, что дискретизации подвергается не сама среда на этапе её моделирования, а уже созданная континуальная модель. Однако в ряде случаев использование дискретных подходов на более ранних этапах разработки и совершенствования моделей механики деформируемого твердого тела может привести к созданию новых эффективных численных схем и алгоритмов. Такие возможности открывает метод анализа полей деформаций и напряжений, использующий в качестве дискретных моделей сплошного тела ориентированные графы. Предлагаемый в работе подход основан на сочетании дискретных и энергетических представлений при моделировании процесса деформирования упругих сред. Дискретная модель при этом строится как первичная модель исследования, а способ перехода от континуума к системе с конечным числом степеней свободы основан прежде всего на физических соображениях. Основой для построения моделей в виде графов служит объективный характер операций измерения выбранных независимых переменных. Как известно, целью любого измерения является установление зависимости между значением переменной и показаниями прибора. В то же время измерение связано как с точками системы, между которыми оно осуществляется, так и с ориентацией измерителя. Эти свойства операции измерения можно представить направленным отрезком, т.е. дугой графа. С точки зрения операции измерения все переменные вне зависимости от их физической природы могут быть подразделены на два типа - параллельные и последовательные [1, 2]. Параллельные переменные характерны тем, что измеритель подключается к двум точкам системы параллельно потоку энергии. Для измерения последовательных переменных прибор должен быть установлен в разрыв энергетического потока. К параллельным переменным относятся линейные и угловые разности перемещений, относительные скорости и ускорения, деформации, разности давлений, температур, электрических потенциалов и т.д. К последовательным переменным относятся силы, моменты сил, электрический ток, тепловой поток, поток жидкости или газа и другие. Выбор связной пары переменных определяется тем, что произведение последовательной и параллельной переменных должно давать скаляр с размерностью мощности или работы [2, 3]. Анализ системы на основе графового подхода сводится к тому, что 1) среда делится на части, имеющие известное математическое описание (в рассматриваемом ниже случае - закон Гука), 2) для каждой части строится подграф (элементарная ячейка), являющийся моделью этой части среды, 3) элементарные ячейки объединяются в граф - модель анализируемого тела, после чего с помощью матриц, характеризующих структуру графа, и уравнений, описывающих элементарные ячейки, получают уравнения системы в целом. Вывод определяющей системы уравнений основан на использовании вершинного и контурного законов Кирхгофа [4]. Известно, что граф является моделью физической или технической системы только в том случае, если выбранные при его конструировании переменные удовлетворяют вершинному и контурному законам [2]. Эти законы были установлены Кирхгофом для электрических цепей применительно к токам и напряжениям. Впоследствии оказалось, что они носят фундаментальный характер и применимы ко всем последовательным и параллельным переменным, вне зависимости от их физической природы. Согласно вершинному закону алгебраическая сумма последовательных переменных, инцидентных любой вершине графа, равна нулю, а в соответствии с контурным законом алгебраическая сумма параллельных переменных на любом замкнутом контуре равна нулю [4]. Дискретная модель деформируемого тела разрабатывалась под влиянием работ Г. Крона [5-8]. Крон, применяя аналоговое моделирование, предложил использовать электрические сети для представления самых разнообразных физических и технических систем. Он разработал эквивалентные электрические сети, предназначенные для моделирования задач гидродинамики, квантовой механики, теории упругости и пластичности, теории электрического и магнитного поля. Основное внимание в своих работах Крон обращал на универсальность возможностей при использовании уравнений, записанных в тензорной форме. Однако в работах не описан алгоритм конструирования эквивалентных электрических сетей и способы определения параметров элементарных ячеек. Подход Крона не получил распространения применительно к задачам механики сплошной среды. В работе [9] Е.Г. Кузовков обратил внимание на то, что при моделировании напряженно-деформированного состояния упругих тел можно обойтись без электромеханической аналогии, если эквивалентную электрическую сеть трактовать как ориентированный граф. Способ конструирования графовой модели, ее конфигурация, применение специальным образом сконструированных матриц для вывода определяющей системы уравнений применительно к плоской и осесимметричной задачам теории упругости предложены Е.Г. Кузовковым и подробно изложены в [9-14]. В работе [11] на основе графовой модели построен двумерный сингулярный элемент для расчета напряженно-деформированного состояния в окрестности особых точек разреза. Использованию графового метода для анизотропных и неоднородных сред посвящены работы [16, 17]. Графовая модель упругой среды в полярной системе координат рассматривалась в [18]. В данной работе построена графовая модель упругой среды для решения пространственной задачи теории упругости. 1. Конструирование элементарной ячейки Способ конструирования графа тела связан с процессом измерения полного и независимого комплекта переменных, которые однозначно характеризуют деформированное состояние элементов, полученных в результате декомпозиции. При построении графа исследуемого тела за исходные удобно взять те переменные, которые измеряются установкой прибора непосредственно на исследуемый объект, т.е. параллельные переменные. Элементарной ячейкой будем называть подграф, соответствующий одному элементу, полученному при разбиении исходной области на мелкие части. При определении конфигурации элементарной ячейки в прямоугольной декартовой системе координат область, занятую телом, покрываем сетью координатных плоскостей x = const, y = const, z = const, между узлами которых устанавливаем гипотетические измерители, определяющие абсолютные нормальные деформации δхх, δуу, δzz и деформации δух, δху, δzх, δхz, δzу, δyz, обусловленные сдвигом. Рассмотрим шестигранный прямоугольный элемент со сторонами расположенными вдоль осей ox, oy, oz соответственно. Комплект переменных, определяющих деформацию элемента, состоит из трех групп. Первая группа включает в себя δхх, δxу, δxz, вторая - δух, δуу, δуz, третья - δzх, δzу, δzz. В каждой группе одна переменная измеряется между парой точек, лежащих на определяющей оси группы, а две другие - между парой точек, находящихся на осях, перпендикулярных к определяющей. Так, например, δхх измеряется между двумя точками оси ох, δxу - между парой точек, лежащих на оси оу, а δxz - между двумя точками оси oz. При этом δхх представляет собой нормальную абсолютную деформацию, а δху, δхz - относительные перемещения, вызванные поворотом грани со сторонами Измерители принимаем в качестве дуг графа, изображаем их направленными отрезками линий, причем ориентация дуг совпадает с ориентацией осей координат. Точки, между которыми проводим измерения, считаем вершинами графа. Здесь подразумевается, что используются идеализированные приборы, т.е. такие, что ножки нескольких из них могут быть одновременно установлены в одну и ту же точку тела. Поскольку одни и те же точки среды участвуют в разных группах измерений относительно осей х, у, z, то при построении графа они представляются разными вершинами. В результате получаем элементарную ячейку прямоугольного шестигранного элемента, имеющую 36 дуг и 24 вершины и состоящую из трех компонент (рис. 1). Рис. 1. Элементарная ячейка, состоящая из трех компонент Fig. 1. An elementary cell consisting of three components Вершины с номерами 22, 23, 24 - корни подграфов соответствующих компонент элементарной ячейки. Каждая из дуг графа отображает одну из следующих пар параллельных и последовательных переменных: ( , ), ( , ), ( , ), ( , ), ( , ), ( , ),( , ),( , ),( , ), где , , и , , , , , - нормальные и тангенциальные внутренние силы (обобщенные напряжения). Таким образом, дугам элементарной ячейки соответствует следующая пара векторов: (1) где символ t обозначает операцию транспонирования, а индекс с указывает на принадлежность переменных отдельной элементарной ячейке. 2. Матрица жесткости прямоугольного шестигранного элемента с линейным полем деформаций Связь элементарных векторов и , а также зависимость их от напряжений и деформаций упругого тела устанавливаем, принимая в качестве инварианта при переходе к дискретной модели энергию деформации произвольного элемента среды объемом V: (2) где и определяются формулой (1). В результате энергию элемента в виде графа можно представить, с одной стороны, в виде с другой стороны, энергия деформации элемента сплошной среды есть (3) Здесь напряжения связаны с деформациями законом Гука (4) или в матричной форме (5) где а матрица упругости представима в блочном виде (6) Матрицы и в (6) таковы: где - упругие постоянные Ламе. При определении энергии деформации (3) требуется интегрирование выражения, содержащего неизвестные деформации и напряжения. Для этого, учитывая дифференциальные зависимости Коши [19], аппроксимируем неизвестные деформации в пределах элемента следующими выражениями: (7) В результате соответствующую аппроксимацию получают Заметим, что поскольку аппроксимируем деформации, а не перемещения, исключается ряд проблем, возникающих при смещениях элементов как твердого тела. Кроме того, значительно сокращается размер матриц, использующихся в расчетах. В силу закона Гука (4) напряжения принимают вид (8) Рассматривая шестигранный элемент в форме прямоугольного параллелепипеда, поместим начало координат в центре элемента. Тогда в пределах элемента Выразим теперь неизвестные коэффициенты в (7) через деформации сторон шестигранного элемента. Как следует из (7) (9) (10) (11) Полагая в (9) поочередно получим при соответствующих комбинациях знаков (12) Из (12) находим, что (13) Аналогично из (10),(11) получаем (14) На основании контурного закона, алгебраическая сумма параллельных переменных на любом замкнутом контуре Ux - компоненты элементарной ячейки равна нулю. На рис. 2 в качестве примера представлены четыре из шести контуров, образованных дугами этой компоненты при описании деформации сторон одной из граней элемента. Рис. 2. Схема выбора контуров при записи уравнений Кирхгофа для параллельных переменных Ux - компоненты ячейки 1 Fig. 2. Scheme of cycles for Kirchhoff’s equations for arcs variables of cell 1 (Ux Component) Записывая все шесть соответствующих уравнений Кирхгофа, получим откуда (15) (16) (17) Складывая левые и правые части уравнений (15), а также (16) и (17), получим соответственно (18) Равенства (18) с учетом (13), (14) дают (19) Тот же результат следует из равенства смешанных частных производных: Интегрированием вдоль соответствующих сторон элемента выражений (7), связанных с uy, uz, аналогично предыдущему, найдем еще 18 коэффициентов: Использование контурных законов для дуг Uy- и Uz-компонент позволяет получить равенства (20) Подставляя напряжения (8) в уравнения равновесия: (21) и пользуясь (19), (20), получим три уравнения относительно оставшихся пока не найденными коэффициентов a1, b2, c3, d2, e1, f3, g1, h3, i2: (22) С целью их определения дополнительно предположим, что (23) Интерпретация этих условий состоит в следующем. Представление правой части в уравнениях (22) в виде суммы двух равных слагаемых позволяет произвести группировку коэффициентов, связанных с нормальными и тангенциальными составляющими тензора напряжений, после чего приравнять каждую группу к нулю. Кроме того, предполагая «равноправность» переменных введем условие на изменение компонент деформаций сдвига. Так, например, условие из (23) означает, что - угол поворота к оси Ох линейного элемента, параллельного оси Оу, изменяется по переменной точно так же, как изменяется по z , т.е. угол поворота к оси Ох линейного элемента, параллельного оси Оz. В этом случае В результате найдем значения оставшихся коэффициентов: где Кроме того, используя (19), (20) имеем С помощью найденных коэффициентов деформации внутри шестигранного элемента можно представить в виде , (24) где (25) Элементы матрицы таковы: (26) Элементы матрицы имеют следующий вид: (27) Подставляя из (24) в (2), получим Отсюда после сокращения на и транспонирования имеем (28) Используя теперь закон Гука ( 5), получим Таким образом, уравнение элементарной ячейки приобретает следующий вид: (29) где матрица жесткости элементарной ячейки определяется выражением Поскольку матрица из (25) транспонированная , а также матрица из (6) блочные, то уравнение (29) можно представить в виде где связь между нормальными составляющими определяется матрицей (30) а вторая часть между тангенциальными составляющими - матрицей (31) Используя матричное умножение и последующее интегрирование, находим элементы матриц (30) и (31). Матрица и имеют следующую структуру: Элементы матриц , ..., получены аналитически в явном виде, но здесь не приводятся из-за их громоздкости. 3. Сходимость аппроксимации графового метода Важно отметить, что континуальную модель сплошной среды удается полностью восстановить из ее дискретной модели в виде ориентированного графа. Покажем, что при уменьшении размеров сетки вершинный и контурный законы приводят к выполнению уравнений равновесия и совместности деформаций теории упругости, что влечет за собой сходимость приближенного решения к точному. Рис. 3. Тело, состоящее из восьми элементов Fig. 3. A solid consisting of eight elements Рассмотрим граф тела, состоящего из восьми элементов (рис. 3). Он состоит из трех компонент Ux, Uy и Uz. Для наглядности каждая компонента графа тела изображена на рис. 4 в виде двух частей. Компоненты Ux, Uy и Uz для элементов верхнего слоя представлены на рис. 4, а, б, в, для элементов нижнего слоя - на рис. 4, г, д, е. Поэтому узлы нижней грани для первой составляющей части и верхней грани для второй составляющей части общие. В соответствии с нумерацией центральная вершина графа тела получила для соответствующей компоненты номер 40, 41, 42 соответственно. Согласно вершинному закону алгебраическая сумма последовательных переменных на каждой вершине графа равна нулю. Если дуга выходит из вершины, то силу, действующую на соответствующую точку элемента, принимаем положительной, если дуга входит в вершину - направление считаем отрицательным. В соответствии с этим для центральной вершины 40 Ux-компоненты получим (32) Здесь нижний индекс показывает номер ячейки, к которой относилась дуга. Разложим значения внутренних сил элемента относительно их значений в смежных элементах в ряд Тейлора. Сохраняя два члена в разложении, получим (33) Рис. 4. Компоненты графа Ux, Uy, Uz: (а), (б), (в) - для ячеек элементов 1, 2, 5 и 6; (г), (д), (е) - для ячеек элементов 3, 4, 7 и 8 Fig. 4. Components of a graph Ux, Uy, Uz: а, б, в for cells 1, 2, 5, 6: г, д, е for cells 3, 4, 7, 8 Подставляя (33) в (32), имеем (34) Пользуясь теперь формулой (28) и матрицей (25) и выражая внутренние силы через напряжения, уравнение (34) после изменения порядка операций дифференцирования и интегрирования представим в виде Подставляя в последнее равенство значения (26), (27), получим откуда при произвольной области V следует первое уравнение из (21). Точно так же уравнения вершинного закона для вершины 41 Uy - компоненты и вершины 42 Uz - компоненты приводят к двум другим уравнениям равновесия (21). Покажем теперь, что контурный закон приводит к уравнению совместности деформаций. Рассмотрим сначала контуры, образованные дугами ячеек 1 и 3 Ux - компоненты графа (рис. 5, а): (35) Здесь и ниже верхний индекс показывает номер ячейки, к которой относится дуга. Поскольку (36) из (35) получим (37) Рис. 5. Схема к записи контурных уравнений: а - для дуг Ux - компоненты ячеек 1 и 3; б - для дуг Ux - компоненты ячеек 1 и 5; в - для дуг Uz - компоненты ячеек 1 и 3 Fig. 5. Scheme of cycles for Kirchhoff’s equations: а - for arcs of cell 1 and 3 (Ux Component); б - for arcs of cell 1 and 5 (Ux Component); в - for arcs of cell 1 and 3 (Uz Component). Разложим деформации в ряд Тейлора, сохраняя в разложении по два члена: Отсюда, используя (36), получим Следовательно, Но Поэтому (38) Разлагая в ряд Тейлора деформации получим Подставляя найденные отсюда выражения для а также (38) в (37), получим (39) Учитывая, что из (39) имеем (40) Рассматривая разность контуров, образованных дугами ячеек 1 и 2 Uy - компоненты графа (рис. 6, а), после аналогичных преобразований получим (41) Рис. 6. Схема к записи контурных уравнений: а - для дуг Uy - компоненты ячеек 1 и 2; б - для дуг Uz - компоненты ячеек 1 и 2 Fig. 6. Scheme of cycles for Kirchhoff’s equations: а - for arcs of cell 1 and 2 (Uy Component); б - for arcs of cell 1 and 2 (Uz Component) Из уравнения (40) и (41), следует или первое из уравнений совместности деформаций в традиционной записи [19] Совершенно аналогично получим: 1) разности контуров ячеек 1 и 5 Ux - компоненты (см. рис. 5, б) приводят к уравнению (42) 2) разности контуров ячеек 1 и 2 Uz - компоненты (см. рис. 6, б) к уравнению (43) 3) разности контуров ячеек 1 и 3 Uz - компоненты (рис. 5, в) к уравнению (44) 4) разности контуров ячеек 1 и 5 Uy - компоненты (рис. 7, а) к уравнению (45) Рис. 7. Схема к записи контурных уравнений Uy - компоненты: а - для ячеек 1 и 5; б - для ячеек 1 и 2 Fig. 7. Scheme of cycles for Kirchhoff’s equations: а - for arcs of cell 1 and 5 (Uy Component); б - for arcs of cell 1 and 2 (Uy Component). Складывая уравнения (42) и (43), а также (44) и (45), получим еще два уравнения совместности деформаций Покажем, что выполняются также оставшиеся три уравнения совместности деформаций [19]. Рассмотрим , например, уравнение (46) С учетом дифференциальных зависимостей Коши запишем его в виде (47) Используя контурный закон, покажем, что (48) (49) (50) (51) Действительно, разность контуров ячеек 1 и 2 Uy - компоненты (рис. 7, б) дает (52) С помощью разложений в ряды Тейлора и равенства получим (53) Пользуясь (53), из (52) получим (54) Так как (55) то после подстановки (55) в (54) и сокращения на приходим к (48). Точно так же уравнение (49) может быть получено из контурных законов, примененных к дугам ячеек 1 и 2 Uz - компоненты, а уравнения (50) и (51) из тех же законов для контуров Ux - компоненты ячейки 1, изображенных на рис. 3. В результате из справедливости равенства (47) следует уравнение совместности деформаций (46). Аналогично доказывается выполнение двух последних уравнений совместности деформаций Таким образом, уравнения равновесия и совместности деформаций, входящие в традиционную постановку задач теории упругости, на графовой модели обеспечиваются автоматически, как следствие присущих графу фундаментальных зависимостей - вершинного и контурного законов. В связи с этим меняется роль уравнений равновесия и совместности деформаций - из основных соотношений при определении напряженно-деформированного состояния они превращаются во вспомогательный инструмент, позволяющий проконтролировать правильность построения графовой модели. Выводы 1. С помощью графового подхода построена матрица жесткости для восьмиузлового шестигранного упругого элемента с 24 степенями свободы при линейной аппроксимации деформаций. В методе конечных элементов (МКЭ) для такой аппроксимации требуется элемент, имеющий 20 узлов (60 степеней свободы) [20]. В результате определяющая система уравнений графового метода содержит уравнений примерно в 3 раза меньше по сравнению с системой, получаемой традиционным способом МКЭ. 2. Показано, что фундаментальные законы теории графов (вершинный и контурный) обеспечивают выполнение условий равновесия и совместности деформаций для любого элемента, полученного при разбиении исходной области на мелкие части.

About the authors

A A Tyrymov

Volgograd State Technical University

Email: tyrymov2010@yandex.ru

References

Statistics

Views

Abstract - 341

PDF (Russian) - 53