Full Text

Введение 1о. В настоящее время интенсивно развиваются технологии, позволяющие создавать слоистые структуры микронного масштаба с предварительно напряженными (самонапряженными) слоями. Такие структуры используются, например, для изготовления сферических микрокапсул, транспортирующих химические вещества в биосистемах [1-7]. Как правило, предварительные напряжения являются неизбежными нежелательными факторами и возникают вследствие физико-химических процессов, протекающих при создании или присоединении слоев. К ним, в частности, относятся: локальный нагрев при эпитаксиальном осаждении [8, 9], перераспределение ионных и вандерваальсовых связей на поверхностях соединяемых тонких пленок при LbL-процессинге [10, 11] и молекулярных агрегатов при самоассемблировании [12], усадки в результате фотохимической солидификации в ходе стереолитографии [13] и пр. Однако в некоторых случаях предварительные напряжения могут быть созданы намеренно для придания желаемой неоднородности и анизотропии материалам слоев. Ярким примером является технология Strained Silicon, позволяющая за счет деформирования тонкого кремниевого слоя на подложке создать заданную анизотропию проводимости. В результате тепловые потери вентилей снижаются и появляется возможность увеличить их плотность на кристалле [14]. Таким образом, помимо задачи моделирования, практически важными оказываются задачи оптимизации предварительно напряженных слоистых структур относительно целевой функции, характеризующей приближение к желаемому распределению внутренних напряжений в теле. При постановке перечисленных выше задач, как правило, предполагается, что: 1) деформации не могут считаться малыми; 2) для всей структуры в целом отсутствует свободная от напряжений отсчетная форма; 3) число слоев может быть очень велико. Математические методы моделирования механики и физики твердого тела, учитывающие эти факторы, далеки от завершения и представляют собой интенсивно развивающуюся область нелинейной механики континуума [15-20]. Поиск свободной от напряжений отсчетной формы для всей слоистой структуры в целом приводит к идее о погружении тела в пространство, геометрические свойства которого более общие, чем у евклидова. Это позволяет, используя дополнительные геометрические параметры, совместить без зазоров и наложений слои, не создавая в них напряжений. При этом параметры будут иметь конечные разрывы на границах слоев, но при большом количестве слоев они, согласно принятой в механике сплошной среды идее континуализации, могут быть аппроксимированы непрерывными функциями. В результате слоистой структуре ставится в соответствие гладкое тело (в смысле Нолла, [21]), отсчетная ненапряженная форма которого определяется в пространстве с неевклидовой непрерывной (аффинной) связностью. Эта связность называется материальной [15, 21], а определенное таким образом гладкое тело - неевклидовым телом, которое эквивалентно некоторому телу с непрерывно распределенными дефектами [22]. Эффективность геометрического подхода показана во многих работах [22-30]. Вместе с тем континуализация слоев, т.е. переход от дискретной слоистой структуры к непрерывной, может быть осуществлен различными способами, и вопросы об их строгом обосновании пока остаются открытыми. По этой причине сопоставление конкретных задач с дискретной структурой и соответствующих им непрерывных задач имеет как прикладной, так и теоретический интерес. Настоящая работа продолжает серию статей, в которых рассматриваются модельные задачи об определении напряженно-деформированного состояния многослойных структур при конечных деформациях и различных условиях нагружения [31-34]. Используются идеи, развитые в теории тел переменного материального состава (растущих тел) [20, 35, 36]. В отличие от предшествующих работ, впервые рассматривается случай сжимаемого гиперупругого материала. Цель настоящей работы - исследование модельной нелинейной задачи для самонапряженного слоистого полого шара, его континуализованного аналога и сопоставление результатов вычислительного моделирования. 2o. В начале работы приведем некоторые общие определения. Под телом [21] (гладким телом) будем понимать континуальное множество материальных точек, снабженное структурой гладкого трехмерного многообразия [37]. Под конфигурацией понимается гладкое вложение [37] тела в физическое пространство , которое, в свою очередь, является гладким трехмерным многообразием. Образ конфигурации, , представляет собой точечное множество в физическом пространстве - форму тела. Одно и то же тело может иметь сколь угодно много форм. Композиция[4] двух конфигураций, , , определяет деформацию - отображение одной формы, , на другую, . Среди всех форм выделяются особые - натуральные - формы, свободные от напряжений. Тело может быть результатом соединения некоторого конечного количества слоев. В этом случае тело будем называть сборкой (из указанных слоев). Сборка, в общем случае, не имеет натуральной формы, но каждый из слоев в отдельности - имеет. Работа имеет следующую структуру. В первом разделе исследуется задача о центрально-симметричных деформациях сборки из конечного числа тонкостенных полых шаровых слоев. Натуральные формы слоев геометрически несовместны: любое их соединение без деформации имеет либо зазоры, либо наложения. Величины этих зазоров или наложений определяются рекуррентно из эволюционной задачи, моделирующей технологический процесс. Аналитическое решение нелинейной краевой задачи для несжимаемого гиперупругого материала слоев (предполагается, что отклик определяется двухконстантным потенциалом Муни - Ривлина) построено в подразд. 1.4. Решение представлено в виде конечных алгебраических выражений, которые содержат корни трансцендентных уравнений, определяемых численно. Начальные приближения для них определяются из решений линейной задачи. В подразд. 1.5 нелинейная задача для сжимаемого материала, отклик которого определяется трехконстантным потенциалом Муни - Ривлина, приводится к системе обыкновенных квазилинейных уравнений и соответствующей двухточечной краевой задаче. Она решается методом пристрелки, в котором начальные значения искомых функций и их производных определяются из задачи для несжимаемого материала. Математически свойство самонапряженности выражается тем, что в отсутствие внешних силовых полей меры деформаций (более точно - локальных деформаций из ненапряженного состояния) отличны от нуля и определены кусочно-непрерывными функциями с разрывами на границе слоев. Вычисления, приведенные в подразд. 1.5, показывают, что при уменьшении толщины слоя и увеличении их числа величины разрывов уменьшаются. Это наблюдение подсказывает идею аппроксимации компонентов мер деформаций непрерывными функциями. В результате такой аппроксимации тело становится неевклидовым и приобретает новую физическую характеристику - меру несовместности локальных деформаций, или, в терминах континуальной теории дефектов, плотность дефектов. Вычислительное сопоставление слоистых структур и соответствующих им неевклидовых тел при различном числе слоев приведено в подразд. 2. 1. Эволюция дискретной структуры 1.1. Тела-слои и тела-сборки 3o. Рассмотрим сборки, состоящие из конечного числа шаровых слоев. Полагаем, что для каждого слоя в отдельности существует конфигурация, образ которой - форма слоя - свободен от напряжений. В общем случае эти формы геометрически несовместны, т.е. их нельзя совместить без зазоров или наложений. Совокупность отсчетных форм слоев сборки играет ту же роль, что и натуральная форма тела в классической теории упругости: относительно этих форм вычисляются меры деформаций, входящие как аргументы в функцию запасенной упругой энергии. В наиболее простых постановках задач для самонапряженных слоистых структур отсчетные формы слоев полагаются заранее известными (в [32, 33] им соответствует случай заданной отсчетной геометрии). Однако, если задача предполагает моделирование некоторого процесса послойного создания трехмерного тела, отсчетные формы его слоев заранее неизвестны и должны быть определены из некоторых дополнительных условий, характеризующих технологию процесса. В рамках настоящей работы будем полагать, что эти условия определяют геометрические соотношения для формы слоя в составе сборки на момент его присоединения и формы того же слоя в ненапряженном состоянии. В частности, такие соотношения характеризуют усадку, происходящую при фотохимической солидификации полимера, температурные деформации при охлаждении осажденного слоя и т.п. Следуя [34, 38, 39], будем называть краевую задачу с неизвестными отсчетными формами слоев, но дополненную условиями для их определения эволюционной задачей. 4o. Для формализации эволюционной задачи рассмотрим семейство гладких трехмерных тел, упорядоченных относительно вложения, т.е. Кроме того, будем полагать, что - подмногообразие , т.е. тело как топологическое пространство, обладает топологией, индуцированной . Образно говоря, это условие означает, что все тела из топологически эквивалентны и могут быть получены непрерывным отображением некоторого «модельного» тела. Для задачи, исследуемой в настоящей работе, это тело - полый шар, скажем, с единичным внешним радиусом и радиусом полости, равным 1/2. Поскольку в этом разделе рассматриваются сборки, состоящие из конечного числа слоев, то множество индексов имеет конечную мощность и без ограничения общности его можно считать конечным отрезком натурального ряда, т.е. , . Элементы семейства , т.е. тела , будем называть промежуточными сборками, - начальным телом, - финальным телом [40]. 5o. Каждая сборка, т.е. тело , является результатом соединения конечного числа слоев. Несмотря на то что представление тела как соединение его частей и, наоборот, разбиение его на конечные части часто используется в механике континуума, применение стандартных теоретико-множественных операций для определения таких разложений затрудняется следующим обстоятельством. Тела , будучи гладкими подмногообразиями финального тела , являются открытыми подмножествами [37]. На первый взгляд, слой сборки может быть определен как теоретико-множественная разность двух смежных тел, и , однако она не является открытым множеством, и поэтому не является гладким телом. Кроме того, теоретико-множественное объединение тела - промежуточной сборки и тела - слоя не будет содержать их общую границу и, следовательно, не может представлять следующую по порядку сборку. Эти сложности могут быть преодолены за счет введения особой операции над множествами, представляющими тела, - операции соединения : (1) Здесь - тело-слой, дополняющее промежуточную сборку до следующей по порядку сборки , - внутренность множества в топологии , - граница множества в топологии (напомним, что - подмножество и его топология индуцируется из ). Формула (1) имеет ясный геометрический смысл: тело-сборка является объединением трех множеств: предшествующей сборки , слоя и общей части границы и . Таким образом, операция соединения формализует процесс «склеивания» двух тел с общей границей, точки которой включаются в состав тела и изменяют свой статус с граничных на внутренние. Разложение промежуточного тела-сборки на составляющие его слои может быть записано в виде [39] (2) Физическое пространство, вмещающее формы тел-сборок и тел-слоев, полагаем трехмерным аффинно-евклидовым, т.е. тройкой , где - множество, элементами которого являются места материальных точек, - трансляционное векторное пространство со скалярным произведением Будем полагать, что формы тел-сборок и тел-слоев подчинены условиям: А. Актуальные формы тел , так же как и составляющих их слоев , - полые шары в , т.е. множества мест и соответственно: Здесь индекс обозначает номер сборки, а индекс - номер слоя. В. Натуральные формы слоев - полые шары в , т.е. С. Внутренний, , и внешний, , радиусы ненапряженной формы тела-слоя связаны с внешним радиусом актуальной формы тела-сборки соотношениями где , - заданные толщины слоев в ненапряженном состоянии; - величины, определяющие зазоры (при ) или наложения (при ) ненапряженных слоев и характеризующие технологическую усадку при солидификации слоя[5]. D. Заданы внутренний и внешний радиусы, , , начального тела (т.е. слоя с номером 0). Здесь и в дальнейшем используются следующие обозначения: одинарный нижний индекс обозначает номер соответствующего тела-слоя или тела-сборки, а двойной нижний индекс, элементы которого разделены запятой, - номер сборки и номер слоя в составе этой сборки. С учетом принятых обозначений и формулы (2) последовательность форм тел-сборок может быть представлена диаграммой[6] (рис. 1). Рис. 1. Последовательность форм тел-сборок Fig. 1. The sequence of assemblies Из предположений A-D следует, что для любого тела-сборки существуют конфигурации , , образы которых - полые шары в , причем части этих шаров, соответствующие слоям с номером , свободны от напряжений. Например, для тела-сборки с номером 2 существует три таких конфигурации; они показаны на диаграмме (рис. 2). Рис. 2. Конфигурации для тела-сборки с номером 2 (выделение темным цветом соответствует форме, свободной от напряжений) Fig. 2. Configurations for assembly with number 2 (stress free shapes are highlighted in dark) 1.2. Координаты и векторные базисы 6o. Поскольку, по предположению, физическое пространство евклидово, то в нем может быть определена декартова система координат с началом в точке и ортонормированным базисом : , ; , где - символ Кронекера. Так как формы всех сборок и тел-слоев обладают центральной симметрией, то, наряду с декартовыми координатами , в удобно использовать сферические координаты , связанные с декартовыми формулами перехода: и отсчетные координаты , связанные с координатами соотношениями (3) Здесь - продолжение функции , являющееся гомеоморфизмом на свой образ. Функция[7] задает расстояние от начала координат до актуальных позиций материальных точек формы тела-слоя в составе сборки в зависимости от их отсчетных позиций в натуральной форме . Каждая из таких координатных систем ассоциирована с формами слоев и играет роль «вмороженных» координат [41]. Таким образом, построен расширенный атлас физического пространства, содержащий, помимо декартовой карты, карты, определяемые различными вариантами сферических координат. Такое расширение позволяет упростить дальнейшие вычисления, используя различные координаты по мере необходимости. Каждая координатная система порождает поле локальных базисов (векторных базисов, по терминологии [42]), т.е. троек векторов и , касательных к соответствующим координатным линиям и определяемым формулами: Дуальные к ним базисы и определяются соотношениями биортогональности: Здесь - символы Кронекера, принимающие значение 1 при совпадении s и c и 0 - в противном случае. Разложения используемых далее векторов локальных базисов по элементам декартова базиса приведены в табл. 1. Таблица 1 Представление локальных базисов Table 1 Representation of local bases 0 0 0 0 1.3. Меры деформаций 7o. Согласно деформация каждого слоя в составе сборки из натуральной формы центральносимметрична (радиальные смещения задаются функцией ). В связи с этим градиент деформации, транспонированный градиент деформации и тензоры Коши - Грина могут быть представлены диагональными диадными разложениями: (4) (5) Главные инварианты тензоров определяются соотношениями: 1.4. Несжимаемый материал 8o. Полагаем, что материал тел-слоев несжимаем, однороден и изотропен[8], а плотность упругой энергии, запасаемой им при деформировании в составе тела-сборки, определяется двухконстантным потенциалом Муни - Ривлина: (6) Константы (упругие модули) , могут быть выражены через другую пару модулей b и μ следующим образом: , . Использование модулей и позволяет выражения, получаемые далее, записать в более лаконичной форме, в то время как модули μ и b указывают на взаимосвязь нелинейной гиперупругой модели материала и линейно-упругой модели: при малых деформациях из ненапряженного состояния μ совпадает с модулем сдвига линейно-упругого несжимаемого материала. Для потенциала (6) формула Ривлина [42, 44] определяет напряжения Коши: (7) где I - единичный тензор; - гидростатическая компонента, определяемая из уравнения равновесия и заданных краевых условий. Краевая задача для центрально-симметрично деформированного тела-сборки имеет вид (8) (9) (10) Краевые условия (9) характеризуют равномерное гидростатическое нагружение сферических границ сборок, а (10) - условия идеального контакта между слоями. 9o. Решения уравнения (8) могут быть определены следующим образом. Поскольку входящие в него функции являются функциями радиальной координаты, то оно является обыкновенным дифференциальным уравнением. Интегрируя его, имеем [42, 44] где - постоянная интегрирования. С учетом соотношения (7) и формул (4), (5) для разности интегрирование дает (11) где - параметр деформации[9]. Полученная формула (11) позволяет определить функцию . Тогда, используя (7), получаем (12) 10o. Постоянные , должны быть выбраны так, чтобы удовлетворялись краевые условия (9) и (10). Подстановка выражений (11) и (12) в них приводит к системам нелинейных алгебраических уравнений. Для начального тела ( ) система состоит из двух уравнений относительно двух неизвестных, , : (13) где . Системы для остальных тел-сборок , (с учетом равенств , ) состоят из уравнений относительно неизвестных , , : (14) Условие С дает уравнение, дополняющее систему (14): (15) которое преобразует недоопределенные по отдельности уравнения (14) в рекуррентную последовательность систем уравнений. Система (13) является базой рекурсии. Решая эту систему, определяем и . По известному значению из (15) определяем . Первый шаг рекурсии - система (14) для - система четырех уравнений относительно четырех неизвестных , . По известному из (15) определяем . Второй шаг рекурсии - система (14) для - система шести уравнений относительно , . Решая последовательно системы нелинейных уравнений, определяем все постоянные интегрирования в (11) и (12). 1.5. Сжимаемый материал 11o. Для сжимаемого материала плотность упругой энергии слоя задается трехконстантным потенциалом Муни - Ривлина в котором , постоянная d соответствует модулю объемного расширения. В этом случае выражение для тензора напряжений Коши примет вид [44] (16) где , а , , - скалярные функции реакции, Элементарные вычисления приводят к следующим выражениям для напряжений Коши (16): (17) Для центральной симметрии уравнение равновесия сводится к одному скалярному уравнению Подставляя выражение (17) в последнее уравнение, приходим к обыкновенному дифференциальному уравнению второго порядка относительно неизвестной функции : Краевые условия представлены выражениями (9), (10). Таким образом, для фиксированного α краевая задача относительно функций имеет вид (18) В качестве начального приближения для решения краевой задачи, соответствующей сжимаемому материалу, используется решение для несжимаемого материала. В свою очередь, это решение получено с использованием решения линейной задачи как начального приближения. Рис. 3, 4 иллюстрируют различия напряженных состояний, полученных для одного полого шара с использованием различных уравнений состояния для линейно-упругого, гиперупругого несжимаемого и гиперупругого сжимаемого материалов. Постоянная гидростатическая нагрузка полагалась равной и На графиках, приведенных на рис. 3, показаны распределения радиальных и окружных компонент напряжений Коши, соответствующие линейной и нелинейной моделям сжимаемого материала (здесь и далее - внутренний радиус начального тела). При этом компоненты тензора напряжений, соответствующие линейной модели, имеют вид [45, 46] Здесь q - гидростатическая нагрузка, - отсчетные радиусы полого шара. На графиках, приведенных на рис. 4, показаны распределения разностей напряжений, соответствующих сжимаемой и несжимаемой нелинейным моделям. Приведенные графики иллюстрируют, что напряженное состояние в сжимаемой модели может существенно отличаться от соответствующих состояний в линейной и нелинейной моделях. Это показывает, что влияние нелинейности отклика и сжимаемости материала может оказаться весьма существенным. Результаты вычислительного моделирования системы полых шаров, ненапряженные формы которых определяются из рекуррентных соотношений (18) при заданной послойной усадке S, иллюстрируются графиками на рис. 5, 6. Рассматривались системы, состоящие из 15 сферических слоев. Параметры гиперупругого материала брались теми же, что и в предыдущем расчете. Усадка S варьировалась в диапазоне (отрицательное значение S соответствует объемному расширению). Распределения радиальных и окружных компонент напряжений Коши, возникающих в самонапряженной структуре в отсутствие внешних силовых полей, показаны на рис. 5. В ходе решения рекуррентной последовательности краевых задач (18) при заданном значении S определяются внешние и внутренние радиусы ненапряженных шаровых форм. Эти формы геометрически несовместны: при их соединении возникают зазоры или наложения. Иллюстрация этого факта при различных значениях S приведена на рис. 6. На нем показаны экваториальные сечения ненапряженных шаровых слоев, которые расположены вплотную друг к другу вдоль радиуса, указанного жирной линией. Вдоль других направлений сечения располагаются с наложениями или зазорами. Радиальные компоненты напряжений Коши Окружные компоненты напряжений Коши q = 0,2μ q = 0,7μ Рис. 3. Сравнительные графики напряжений, полученных из линейно-упругой (пунктирная линия) и гиперупругой моделей (сплошная линия) Fig. 3. Comparative graphs of stresses for linear-elastic (dotted line) and hyperelastic models (solid line) Радиальные компоненты напряжений Коши Окружные компоненты напряжений Коши q = 0,2μ q = 0,7μ Рис. 4. Сравнительные графики напряжений, полученные из гиперупругих моделей для несжимаемого и сжимаемого материалов Fig. 4. Comparative graphs of stresses obtained for hyperelastic models of incompressible and compressible materials Радиальные компоненты напряжений Коши Окружные компоненты напряжений Коши S = -0,4 S = 0,1 S = 0,6 S = 1,6 Рис. 5. Распределения напряжений в самонапряженной структуре в отсутствие внешних силовых полей Fig. 5. The distribution of stresses in self-stressed structure with an absence of external force fields S = -0,4 S = 0,1 Отсчетная форма S = 0,6 S = 1,6 Отсчетная форма Рис. 6. Наложения ненапряженных отсчетных форм (изображены экваториальные сечения ненапряженных шаровых слоев) Fig. 6. Overlaps of the stress-free reference shapes (equatorial sections of non-stressed spherical layers are shown) 2. Непрерывные неевклидовы структуры 12o. Для определения напряженно-деформированного состояния слоистых структур приходится решать системы нелинейных уравнений (18), что для большого числа слоев становится вычислительно сложной задачей. В то же время «измельчение» слоев, т.е. увеличение их числа и уменьшение толщины при фиксированной суммарной толщине и достаточно медленном изменении параметров процесса послойной сборки, , , , приводит к распределениям для напряжений и деформаций, близким к непрерывным. В связи с этим представляется целесообразным для анализа слоистых структур с большим количеством слоев предварительно осуществлять их континуализацию, позволяющую заменить разрывные распределения для компонент меры деформаций Коши - Грина непрерывными аппроксимациями, близкими к исходным в некотором, определенном заранее, смысле. В рамках настоящей работы в качестве таких аппроксимаций используются полиномы, коэффициенты которых определяются из условий наименьших квадратичных отклонений от разрывных распределений. Аппроксимируются диагональные компоненты тензорного поля характеризующего послойно несовместные деформации финального тела: (19) где и - кусочно-непрерывные функции актуальной формы: (20) Функцию (19) можно рассматривать как «склейку» полей градиентов деформации каждого слоя из ненапряженного состояния, записанных в переменных актуальной формы. Поскольку в этих переменных области определения «склеиваемых» полей - смежные подобласти актуальной формы , то результатом будет непрерывное тензорное поле, определенное на отрезке Обозначим через А и а аппроксимации функций и , заданных формулой (20). Им соответствует тензор F - аппроксимация (21) Для сравнения рассматриваются тела-сборки, состоящие из 5, 10, 20, 100 и 200 слоев, суммарная толщина которых в ненапряженном состоянии равна Упругие модули принимались равными , ненапряженные формы слоев совпадают с ненапряженной формой подложки, Аппроксимирующие многочлены приведены в табл. 2. Сопоставление графиков разрывных функций и аппроксимаций для 5, 10, 30 слоев показано на рис. 7. Для большего числа слоев графики этих функций типографически неразличимы. 13o. В результате перехода к непрерывным функциям , финальное тело-сборка становится неевклидовым телом в том смысле, что его натуральная (свободная от напряжений) форма может быть определена только в пространстве с неевклидовой связностью. Таблица 2 Аппроксимирующие многочлены Table 2 Approximating polynomials α A a 5 10 30 50 100 200 α = 5 α = 10 α = 30 А а Рис. 7. Сопоставление дискретных и непрерывных распределений Fig. 7. Comparison of discrete and continuous distributions Согласно терминологии [21] такая связность называется материальной и может быть определена различными способами. Один из них - синтез материальной метрики [21, 35, 36] - состоит в том, что каждой материальной точке финального тела сопоставляется метрика (соответствующего касательного пространства), определяемая как где K - тензорное поле, определенное на , называемое полем локальных конфигураций [21]. Поле K «склеено» из градиентов конфигураций , принадлежащих континуальному семейству , т.е. ; здесь Т - операция взятия дифференциала [37]. Каждый элемент семейства, , переводит инфинитезимальную окрестность точки в ненапряженное состояние. Соотношение между локальной конфигурацией в точке и отображением иллюстрируется диаграммой Актуальная форма тела , полученного в результате континуализации, совпадает с актуальной формой тела , т.е. с . Пусть - актуальная конфигурация. Соотношения между телом , его актуальной формой и одной из его отсчетных форм иллюстрирует диаграмма На диаграмме: - отображение вложения, действующее по правилу . Отображение - деформация, переводящая свободную от напряжений окрестность точки в ее актуальную форму. Замечание 1. В результате континуализации свободные от напряжений формы слоев, составляющих финальное тело , вырождаются в сферы. Таким образом, форма обладает следующим свойством: не только точка , но и точки, лежащие вместе с ней на некоторой сфере, свободны от напряжений. 14o. Для получения координатного представления полей, заданных на теле , выберем в нем систему локальных координат. Используя картрирование актуальной формы сферическими координатами и актуальную конфигурацию , определим координатное отображение на теле . Здесь - координатное отображение, соответствующее сферическим координатам. Обозначим координатный репер на теле , соответствующий , через Это векторное поле на ; в каждой точке тела упорядоченная тройка является базисом касательного пространства . Через обозначим корепер на , связанный с координатным репером равенствами В паре выбранных систем координат на и актуальная конфигурация и ее градиент имеют представления в виде тождественных отображений 15o. Поскольку, в силу (21), то метрика G имеет диадное разложение: (22) Построенная таким образом метрика на всем теле превращает его в риманово пространство, связность которого определяется по формуле Леви-Чивита: Для рассматриваемой задачи отличные от нуля коэффициенты связности определяются следующим образом: (23) Отличие этой связности от евклидовой характеризуется скалярной кривизной (инвариантом Риччи): (24) На рис. 8 показано распределение скалярной кривизны (24) для аппроксимационного многочлена, построенного на основе дискретного распределения при Рис. 8. Скалярная кривизна для случая α = 200 Fig. 8. Scalar curvature for the case α = 200 16o. Тензор напряжений Коши , соответствующий непрерывному телу , задан на актуальной форме и является аппроксимацией тензора напряжений Коши Т, компоненты которого вычисляются по формулам (17). Таким образом, (25) где На рис. 9 приведены распределения физических компонент тензора напряжений на отрезке радиальных координат актуальной формы. Эти распределения сопоставляются с соответствующими распределениями компонент тензора напряжений (17) в дискретных структурах при разных значениях N. Рис. 9 иллюстрирует, что с увеличением N распределения стремятся к . 17o. Определение отсчетной ненапряженной формы финального тела с заданной на нем связностью (пусть даже неевклидовой) позволяет определить напряжения Пиолы [16, 19]. Для того чтобы получить представления в «отсчетных» координатах для неевклидова тела , предполагается использовать формализм геометрической механики континуума [16, 19, 22]. При этом материальное тело и физическое пространство рассматриваются как однотипные объекты - гладкие многообразия с определенной на них аффинной связностью. Связность физического многообразия евклидова, а связность на многообразии (материальная связность) - неевклидова и определяется коэффициентами связности (23). a б Рис. 9. Радиальные (a) и окружные (б) физические компоненты (1 - 5 слоев, 2 - 10 слоев, 3 - 30 слоев, 4 (непрерывная линия) - непрерывная неевклидова структура) Fig. 9. Radial (a) and circumferential (b) physical components of (1 is 5 layers, 2 is 10 layers, 3 is 30 layers, 4 (solid line) is non-Euclidean solid) Напряжения Коши Т интерпретируются как ковекторнозначные 2-формы [19], а напряжения Пиолы Р - как pullback этих форм относительно конфигурации : (26) Замечание 2. Уравнения баланса импульса в переменных карт на многообразии имеют вид [19] Здесь - оператор Картана, являющийся обобщением операции дивергенции; В - ковекторное поле, соответствующее плотности массовых сил; V - скорость движения частиц тела в материальном описании; - форма объема, порождаемая метрикой тела G; - музыкальный изоморфизм [37], соответствующий метрике G, переводящий ковекторные поля в векторные поля. Важно отметить, что уравнения баланса определены на всем теле , а не послойно, что позволяет осуществлять эффективные вычисления для большого числа слоев. 18o. Формула (26) позволяет по заданному тензору напряжений Коши определить тензор напряжений Пиолы. Именно, пусть - локальные координаты на теле , а - криволинейные координаты в . Тогда компоненты тензора напряжений Пиолы, определяются формулой [39] Здесь , и . Функции - ковариантные компоненты тензора напряжений Коши (25). Тензор g представляет метрику на , а тензор G - материальную метрику на . Тензор - транспонированный к обратному тензору , определяемый из равенства Для вычисления компонент тензора напряжений Пиолы выберем в физическом пространстве сферические координаты , а на теле - координаты согласно п. 140. Тогда где определены формулами (25). Метрика g имеет представление а G представлена формулой (22). В этой связи , где учтено, что . Поскольку матрица компонент единичная, то и мы приходим к окончательным выражениям для компонент тензора P: Заключение 19o. При большом числе слоев самонапряженной слоистой структуры кусочно-непрерывные параметры деформации могут быть аппроксимированы непрерывными функциями. При такой аппроксимации отсчетная форма может быть определена в пространстве с аффинной (неевклидовой) связностью. Получаемая при этом математическая модель самонапряженного тела эквивалентна нелинейной модели тела с распределенными дефектами. В рамках этой модели методами геометрической механики континуума могут быть определены обобщенные представления для напряжений Коши и Пиолы. Вычисления, проводимые для дискретной структуры и тела с неевклидовой отсчетной формой, определяемой аппроксимацией параметров деформации, численно иллюстрируют сходимость решений дискретной модели к непрерывной при увеличении числа слоев с неизменной суммарной их толщиной.

About the authors

S A Lychev

Ishlinsky Institute for Problems in Mechanics RAS; Bauman Moscow State Technical University

T N Lycheva

Ishlinsky Institute for Problems in Mechanics RAS

K G Koifman

Bauman Moscow State Technical University

References

Statistics

Views

Abstract - 605

PDF (Russian) - 454