Decomposition of systems of equations for continuum mechanics 1. Elasticity, thermoelasticity, poroelasticity

Abstract


The work is devoted to the development of decomposition methods for systems of linear partial differential equations that arise in continuum mechanics, in particular, in the theory of elasticity and thermoelasticity and poro-elasticity. These methods are based on the decomposition (splitting) of systems of coupled equations into several independent equations. The decomposition significantly simplifies the qualitative study and interpretation of the most important physical properties related to three-dimensional equations and allows an effective study of their wave and dissipative properties. Moreover in certain cases the decomposition makes it possible to find exact analytical solutions of the corresponding boundary and initial-boundary value problems, and greatly simplifies the application of numerical methods, allowing us to use the appropriate routines for simpler equations and independent subsystems. In the first part of the work various systems of equations, including equations of elasticity theory in the form of Tedone and in the form of Beltrami-Donati-Michell are given, their dynamic generalizations are proposed, and various forms of the equations of classical and hyperbolic thermoelasticity as well as the equations poroelasticity are described. A number of historical facts, which are directly related with the considered questions and weakly reflected in Russian literature,are presented. Various types of decomposition and their generalizations are described. The representation of solutions of dynamical systems of equations resulting from the toroidal-poloidal decomposition, as well as the decompositions of Green-Lamé, Cauchy-Kovalevski-Somigliana, Naghdi-Hsu-Chandrasekharaiah, and Teodorescu types are discussed in details. Special attention is given to their static analogues. A generalization of the representation of Savin for the dynamic equations of elasticity is obtained. The representations in curvilinear coordinates, in particular, the representations of Boussinesq, Timpe, Love, Michell, and Muki types are given. The bibliographical references to the original papers are listed.

Full Text

Введение Различные способы декомпозиции систем уравнений механики деформируемого твердого тела и механики жидкости, приводящие их к более простым уравнениям (в частности, несвязанным), предлагались с середины XIX века. Ряд таких преобразований приведен в классических монографиях, например [1-6], однако их обоснование, как правило, основано на интуитивно понятных частных свойствах рассматриваемых систем. Кроме того, основноe внимание уделялось уравнениям линейной теории упругости. В настоящей работе предлагается систематический подход к декомпозиции систем уравнений различных классов трехмерных линейных (и модельных нелинейных) систем уравнений, которые используются в теории упругости, термоупругости и термовязкоупругости, в механике вязких и вязкоупругих несжимаемых жидкостей и сжимаемых баротропных газов. Статья состоит из двух частей*. Первая часть имеет обзорный характер. В ней приведены и систематизированы известные методы декомпозиции уравнений упругости, термоупругости и пороупругости. Первая часть также содержит материал исторического характера: указан ряд фактов, которые непосредственно связаны с рассматриваемыми в статье представлениями решений и слабо отражены в русскоязычной литературе. Отметим, что исторические аспекты довольно полно изложены в многотомном справочном издании Handbuch der Physik (Encyclopedia of Physics), издаваемом под редакцией S. Flügge с 1926 по 1988 год, в частности, в статьях [2, 7, 8]. 1. Теорема разложения Гельмгольца и ее обобщения Большинство преобразований опирается на представление достаточно гладкого векторного поля u в форме разложения на безвихревую u1 и соленоидальную u2 части: (1.1) Это разложение может быть получено из разложения Stokes-Helmholtz (1.2) Скалярная функция φ носит название скалярного потенциала, векторная функция ψ - векторного потенциала, ϑ - гармонической части разложения. Компоненты φ, ψ, ϑ представления (1.2) находятся из выражений (1.3) где операторы и , определяются следующим образом: Здесь a - скалярная функция; u - векторная функция; n¢ - внешняя единичная нормаль к границе области, в которой определяется разложение (1.2). Часто бывает удобным включить гармоническую часть ϑ в слагаемые, отвечающие скалярному и векторному потенциалам. Тогда представление (1.2) принимает сокращенный вид: (1.4) Здесь , - скалярный и векторный потенциалы, которые определяются операторами : (1.5) Очевидно, что первое слагаемое представления (1.4) определяет безвихревую часть разложения (1.1) u1, а второе - соленоидальную u2, т.е. Далее не будем различать обозначения потенциалов (1.3), (1.5) и опустим знак тильды, полагая, что их определение понятно из контекста. Формулировка и обоснование разложения (1.2) в декартовых координатах для ограниченной области евклидова пространства впервые была приведена в работе Stokes в 1849 г. [9, p. 9]. Заметим, что эта работа посвящена теории дифракции, которая разрабатывалась Stokes’ом на основе уравнений малых колебаний упругой среды. Чуть позже, в 1858 г., в работе [10, p. 38] Helmholtz использовал аналогичное представление (в декартовых координатах) для поиска некоторых точных решений линеаризованных уравнений гидродинамики. Впоследствии разложение Stokes-Helmholtz было сформулировано для произвольных криволинейных координат, в частности, в прямой бескоординатной форме (Gibbs, 1901 [11, p. 234]). Отдельно рассматриваются случаи, когда область оказывается неограниченной. При этом, в силу известной теоремы Лиувилля ([12, с. 270] или [13, с. 111]), гармоническая функция ϑ, при условии ее ограниченности и знакопостоянности, оказывается константой. В связи с этим разложение Stokes-Helmholtz (1.2) определяется с точностью до постоянной, которая должна обращаться в нуль, если дополнительно указываются условия для поля на бесконечности. Кроме того, для существования интегралов, определяющих скалярный и векторный потенциалы, необходимо указать асимптотические условия для поля u на бесконечности. В современных руководствах по математической физике (например, [14, p. 97]) приводится классическое доказательство для неограниченной области с условиями (1.6) Условия (1.6) существенно ограничивают применимость разложений, в частности, в задачах механики континуума, так как асимптотическое поведение, например, полей перемещений на бесконечности может иметь меньший порядок убывания. Кроме того, условие (1.6) исключает жесткое перемещение среды. Условия (1.6) могут быть уточнены, а именно, может быть указан асимптотический порядок для «деформационной» компоненты представляемого поля. В частности, Gurtin [15] в 1961 г. показал, что представление Stokes-Helmholtz (1.2) может быть использовано для более широкого класса полей, удовлетворяющих следующей асимптотической оценке: Здесь c - постоянное векторное поле, определяющее жесткое перемещение среды. Примеры использования представления Stokes-Helmholtz (1.2) как в классической интерпретации, т.е. для представления достаточно гладких векторных полей, так и в интерпретации, адаптированной к теории операторов в гильбертовых пространствах, в задачах гидромеханики приведены, например, в работах О.А. Ладыженской [16], А.Д. Полянина [17-19], Morino [20], О.В. Воинова [21] и др. Cоленоидальную часть (1.2), определяемую ротором векторного потенциала ψ, можно представить в виде двух слагаемых, определяемых ротором и повторным ротором произведений достаточно гладких скалярных функций ϕ, χ и радиус-вектора r, что приводит к следующему выражению: (1.7) Представление (1.7) впервые было введено в теории магнитного поля Земли и имеет название полоидально-тороидальной декомпозиции (Backus [22]). Доказательство того, что такое разложение может представить произвольную достаточно гладкую векторную функцию, определенную в области, гомеоморфной шару, получено в ряде работ, в частности Kaiser [23], Padmavati, Amaranath [24], Chadwick, Trowbridge [25]. Отметим еще один способ представления гладкого векторного поля через три скалярные функции φ, ϕ, χ [26]: (1.8) Преимущества представлений (1.7) и (1.8) состоит в том, что они не избыточны: три компоненты поля представляются через три скалярные функции, в то время как представление Stokes-Helmholtz (1.2) содержит четыре функции. 2. Различные формы уравнений линейной теории упругости Уравнения движения Navier-Cauchy для однородной изотропной среды имеют вид [2, 5, 6, 27] (2.1) Здесь - векторное поле перемещений; - вектор места; - упорядоченная тройка декартовых координат; t - время; λ, μ - упругие модули Ламе (μ - модуль сдвига); ρ - плотность массы. Точками обозначены производные по времени t. Символом Δ обозначен оператор Лапласа, (2.2) символами - градиент дивергенции, (2.3) Векторное поле внешних (объемных) сил обозначается символом f. Вместо упругих модулей Ламе часто используются «технические» модули (соотношения между различными модулями (см. например [5, с. 114] или [2, с. 79]), а именно модуль упругости E и коэффициент Пуассона ν: Из термодинамических соображений, сводящихся к условию того, что при любой бесконечно малой деформации из естественной конфигурации производится положительная работа, следуют неравенства [2, 5] Из этих неравенств вытекает, что, вообще говоря, коэффициент Пуассона ν может быть отрицательным (вещества с такими свойствами - ауксетики рассматриваются, например, в [28]). Уравнения (2.1) были впервые получены Navier (1821) из представлений о молекулярной микроструктуре сплошной среды, которая соответствовала так называемой одноконстантной модели (μ = λ) [29]. Общая форма (для μ ≠ λ) была получена Cauchy (1822, 1828) [30]. Исторические аспекты получения и исследования этих уравнений приведены в [31, 32]. Для анализа уравнения (2.1) целесообразно ввести безразмерные переменные где R - некоторый характерный линейный размер тела (например, диаметр шара, описанного около тела). В безразмерной форме уравнение (2.1) имеет вид (2.4) Здесь введены обозначения Безразмерная форма уравнений (2.4) является наиболее компактной и систематически используется в настоящей работе. Различие в обозначениях размерных и безразмерных величин (тильда) в (2.4) и далее опускается. Уравнение Navier-Cauchy (2.4) можно записать в формальной операторной форме: Эта форма записи будет использована для построения общих представлений решений. Построение аналитических решений краевых задач, определяемых уравнениями (2.4), осложняется тем, что они одновременно содержат оператор Лапласа (2.2) и градиент дивергенции (2.3). Различные способы представления решения в форме аналитических комбинаций решений более простых уравнений, таких как уравнения Лапласа, Пуассона, волнового уравнения и их степеней, связаны с различными формами уравнения (2.4). Кратко остановимся на некоторых из них. 1. С учетом тождества уравнение (3.4) может быть преобразовано к виду (2.5) Для формулировки уравнений удобно использовать волновые операторы (операторы Д’Аламбера) , Тогда уравнение (2.4) может быть записано в форме Уравнение (2.5) в размерном виде может быть записано следующим образом: где - скорость продольных волн (p-волн, первых волн, primary waves, т.е. первых волн, регистрируемых в геофизике первыми); - скорость поперечных волн (s-волн, secondary waves). В статическом приближении уравнения (2.4), (2.5) приводятся к виду Вскоре после формулировки уравнений движения (2.1) была отмечена их связь с волновыми уравнениями, а именно при b = 0. На первое из этих уравнений указал Lamé (1852) [27], а на второе - Cauchy (1840) [33]. 2. Вычисляя дивергенцию и ротор левой части уравнения (2.5), приходим к скалярному и векторному уравнениям: Последнее уравнение с учетом соотношения преобразуется к виду Если ввести новые искомые функции первая из которых имеет смысл дилатации (относительного изменения объема), то приходим к двум независимым уравнениям (2.6) Из последнего уравнения вытекает, что если и поле b удовлетворяет условиям то поле u в статическом приближении является бигармоническим, т.е. удовлетворяет уравнению На это обстоятельство впервые указал Cauchy (1828) (см. также [2]). 3. Из первого соотношения (2.6) и тождества вытекает выражение для градиента дилатации Подставляя результат в уравнения Navier-Cauchy, приходим к уравнениям в форме В статическом приближении приходим к двум независимым уравнениям (Tedone [34], 1903 г.) где введены следующие векторные переменные: 4. В цилиндрических координатах соответствующие нормированные (физические) локальные базисы определяются соотношениями Векторная форма уравнений Navier-Cauchy (2.4) остается инвариантной, однако явные выражения для разложений полей и дифференциальных операторов записываются в терминах разложений по локальным базисам, т.е. В осесимметричном случае уравнения Navier-Cauchy эквивалентны системе двух уравнений где Взаимосвязь различных форм уравнения Navier-Cauchy (2.4) показана на рис. 1. Еще один класс уравнений, связанный с уравнениями Navier-Cauchy (2.4), использует дифференциальный оператор несовместности Ink, который характеризует несовместность свободных деформаций ε (под свободными деформациями понимается гладкое тензорное поле, которое, вообще говоря, не является градиентом какого-либо гладкого векторного поля [35]): (2.7) ¯ ¯ ¯ ¯ ¯ ¯ (преобразование Тедоне) ¯ Рис. 1. Взаимосвязь уравнений линейной изотропной теории упругости Действие оператора Ink на деформации, порожденные некоторым векторным полем перемещений u, дает нуль, т.е. Здесь - оператор, вычисляющий симметрическую часть градиента u («деформация» векторного поля). Тензорное поле А характеризует несовместность свободных деформаций, которая является источником внутренних напряжений. В теории дефектов эта величина может быть ассоциирована с плотностью дислокаций [35]. Будем полагать, что тензорное поле А симметрично, и, следовательно, симметрично поле ε. Для произвольного симметричного поля ε имеет место тождество [2, c. 41], [26] (tr ε - след тензора ε): (2.8) Вычисление следа левой и правой части (2.7) с учетом этого тождества приводит к равенству (2.9) Полагая, что рассматривается линейно-упругая изотропная среда, воспользуемся законом Гука для выражения деформаций ε через напряжения T (2.10) Подстановка (2.10) в (2.7) с учетом тождеств (2.8), (2.9) приводит к уравнению (2.11) Полученное уравнение определяет связь между мерой несовместности свободных деформаций, определяемых симметричным тензорным полем A и симметричным полем напряжений T, т.е. систему шести уравнений. Эта система может быть преобразована к иному виду, если воспользоваться соотношением, связывающим след левой части и след правой части уравнения (2.11), т.е. (2.12) Подставляя выражение для из (2.12) в (2.11), приходим к уравнению (2.13) Уравнения (2.11), (2.13) не связаны с уравнениями равновесия и определяют несбалансированное поле напряжений, реализуемое, вообще говоря, только при некотором специально подобранном поле внешних сил. Для получения уравнений, характеризующих сбалансированные относительно заданного поля внешних сил f напряжения, используем уравнения баланса импульса (2.14) Исключая из системы уравнений (2.13), (2.14) след тензора напряжений, получаем уравнение (2.15) Уравнение (2.15) характеризует сбалансированное (относительно заданного поля внешних сил f) поле напряжений, которое возникает в теле вследствие двух причин: внешних f и внутренних A. Исключение напряжений из уравнения (2.15), которые могут быть выражены через перемещения u из напряженной отсчетной конфигурации в актуальную согласно закону состояния приводит к уравнениям относительно u: Здесь T0 - напряжения в отсчетной конфигурации. Левая часть полученного уравнения допускает факторизацию (2.16) Как и следовало ожидать, приходим к уравнению Navier-Cauchy (2.1) с дополнительным неоднородным членом fA («фиктивной нагрузкой») который, в свою очередь, удовлетворяет уравнению Исключение из (2.13) ускорений с помощью соотношений приводит к уравнению в форме (здесь ) (2.17) В частном случае, когда внутренние источники напряжений отсутствуют уравнение (2.17) принимает вид (Teodorescu [36]) Уравнение (2.15) при переходит в динамический вариант уравнения Beltrami-Michell: (2.18) В статическом приближении приходим к уравнениям (2.19) Однородное уравнение (2.19) (для f = 0) было получено в 1892 г. Beltrami [37], а неоднородное в 1894 г. Donati [38] и, независимо, в 1899 г. Michell [39]. Правильнее было бы называть уравнением Beltrami-Donati-Michell, а не Beltrami-Michell, как это принято в литературе. В безразмерных переменных уравнения преобразуются к виду, представленному в табл. 1. В записи, как и ранее, знак тильды опущен, а безразмерные коэффициенты определяются соотношениями Таблица 1 Безразмерные формы уравнений типа Beltrami-Michell Номер уравнения Уравнение в безразмерной форме (2.11) (2.13) (2.15) (2.16) (2.17) (2.19) 3. Представления решений динамических уравнений линейной теории упругости 1. Для построения решений неоднородного уравнения движения исторически первым было предложено представление Green-Lamé [27]. В рамках этого представления предполагается, что рассматривается конечная часть среды (регулярная в смысле Келлога [40]), в которой массовые силы могут быть представлены разложением Гельмгольца (1.4): (3.1) Искомое векторное поле перемещений u, которое предполагается достаточно гладким, может быть представлено разложением (3.2) Здесь - скалярное и векторное поля (скалярный и векторный потенциалы перемещений). Представление, аналогичное представлению Green-Lamé для одноконстантной теории, впервые было указано Poisson (теорема Пуассона, 1829) [41]. Подстановка представления (3.2) в уравнения Navier-Cauchy (2.4) преобразует их к виду Поочередное вычисление дивергенции и ротора левой и правой частей этого уравнения приводит к двум независимым уравнениям которые с учетом разложения (3.1) могут быть записаны следующим образом: Частные решения этих уравнений определяются решениями уравнений Ясно, что потенциалы, получаемые как решения этих уравнений, определяют решение, удовлетворяющее уравнению движения (т.е. представляют некоторое частное решение). Возможность представления любого решения уравнения Navier-Cauchy (2.4) определяется полнотой представления (3.2). Clebsch (1863) первым поставил вопрос о полноте представлений Green-Lamé и предложил вариант доказательства. Значительно позже Sternberg (1960) показал необходимость доработки этого доказательства (в оригинальном доказательстве Clebsch предполагал, что граница тела представляет собой одну гладкую замкнутую поверхность). Альтернативное, но все же неполное доказательство было дано Kelvin (1884). Доказательство в завершенном виде было получено Somigliana (1892). Следует отметить, что различные варианты доказательства позже предлагались Tedone (1897), Duhem (1898), Sternberg and Gurtin (1960, 1962) и другими [2]. 2. Решение однородного векторного волнового уравнения может быть выражено через решения скалярного волнового уравнения посредством тороидальной-полоидальной декомпозиции (1.7). В отсутствие объемных сил поле перемещений может быть представлено через три скалярных функции, удовлетворяющих волновым уравнениям где r - радиус-вектор точки, а φ, ζ, ω - скалярные функции, удовлетворяющие однородным волновым уравнениям Для произвольного поля скалярные функции φ, ζ, ω являются решениями неоднородных уравнений причем ϕ определяется как скалярный потенциал разложения векторного поля b, а скалярные функции ζ, ω удовлетворяют следующим уравнениям Здесь Δ0 - трансверсальная часть оператора Лапласа относительно радиального направления сферических координат, Математическое представление такой декомпозиции известно давно, однако его приложение к представлениям в теории упругости связывают с работой Chadwick, Trowbridge (1967) [25]. Авторы работы [25] доказали полноту этого представления для области, ограниченной концентрическими сферами (см. также [15, 23, 42]). 3. Другой класс представлений определяется так называемым решением Cauchy-Kovalevski-Somigliana: где g и определяется уравнением Подобное представление решения было впервые предложено Cauchy и впоследствии (независимо) С. Ковалевской [43]; см также Iacovache [44]. Доказательство полноты осуществлено Somigliana (1889) [45] и в других вариантах - Sternberg & Eubanks (1957) , Sternberg & Gurtin (1962) [2]. Следует отметить, что взаимосвязь этих представлений вытекает из теоремы Боджио (T. Boggio) [2]. Пусть g - решение однородного уравнения Тогда g может быть представлена как сумма где Основываясь на этой теореме, Sternberg показал, что при решения Green-Lamé и Cauchy-Kovalevski-Somigliana связаны соотношениями 4. Четвертую группу образуют представления Boussinesq-Papkovitch-Neuber, которые имеют вид где ψ и ϕ определяются уравнениями Возникает вопрос, является ли решение, найденное при помощи только трех функций Папковича-Нейбера, полным. Пионерские исследования этого вопроса принадлежат Слободянскому [46, 47]. Отметим работы Н.И. Остросаблина [48-52], в которых вопросы полноты исследовались методами анализа симметрий дифференциальных уравнений. См. также [53, 54]. 5. Еще один класс представлений образуют решения Naghdi-Hsu-Chandrasekharaiah [55-57]: позволяющие выразить поле перемещений через вектор-функцию ψ, определяемую уравнением 6. Решение динамических уравнений Beltrami-Michell (2.18) может быть представлено в форме [36], которая в безразмерных переменных принимает вид Здесь функции Ω, Ф определяются уравнениями Соответствующие перемещения находятся по формуле Следует отметить некоторые обобщения классических представлений решений уравнений (2.4), приведенные в работах [58, 59]. Соотношения, определяющие представления различных типов, сведены в табл. 2. Таблица 2 Представление решений динамической теории упругости № Представление перемещений Уравнения для элементов представления Вспомогательные уравнения 1 2 Отсутствуют 3 Отсутствуют 4 Отсутствуют 5 6 Отсутствуют 4. Представления решений статических уравнений линейной теории упругости 1. Первые способы декомпозиции искомой вектор-функции u для частных классов задач принадлежат Boussinesq [60] (1885) и Somigliana [45] (1889). В монографии [60] Boussinesq привел систему частных решений («базовую систему решений» по терминологии [1]), представленную семью гармоническими функциями : Считается, что решения u5, u6, u7 независимо от Boussinesq были повторно найдены П.Ф. Папковичем (1932) [61] и Neuber (1934) [62]. Вместе с тем следует отметить, что в своей работе П.Ф. Папкович ссылается на монографию Pigeaud [63], в которой цитируется работа [60] и воспроизводится усеченная система решений для случая осевой симметрии. Ясно, что решения Boussinesq не являются независимыми. Путем непосредственного вычисления легко убедиться, что имеют место следующие соотношения [64]: (4.1) (4.2) Заметим, что если воспользоваться переменными то соотношения (4.1), (4.2) могут быть записаны в более компактном виде: Таким образом, существует избыточное число «базовых решений», что позволяет представить поле перемещений неединственным образом. Действительно, легко убедиться в том, что, например, 2. Статическим аналогом решения Cauchy-Kovalevski-Somigliana считается решение Галеркина, которое имеет вид В результате подстановки этого соотношения в уравнения Navier-Cauchy (2.4) приходим к неоднородному бигармоническому уравнению относительно вспомогательной вектор-функции w - так называемого вектора Галеркина Б.Г. Галеркин, опубликовав в 1930 г. полученное им решение дифференциальных уравнений теории упругости [65] и считая свое решение общим, не привел доказательства общности полученного им решения, кроме того лишь факта, что ряд известных решений может быть получен из найденного. См. также [66]. Представление Галеркина может быть получено из операторного формализма. Рассмотрим пару операторов, первый из которых А определяется дифференциальными уравнениями Navier-Cauchy (2.4), а второй имеет подобную структуру, т.е. Без потери общности будем полагать, что область определения этих операторов определяется однородными краевыми условиями на границе Число k¢ выбирается из условия наиболее простого представления композиции этих операторов. Поскольку то наиболее простая форма определяется бигармоническим оператором ΔΔ, для чего следует положить Полагаем, что обратим. При этом число k¢ не может быть равно -1 и ∞ (ср. рассуждения Папковича). Тогда тождественный оператор I может быть представлен выражением а произвольная функция u, определенная в может быть представлена в виде разложения Обоснование представления Галеркина осуществил П.Ф. Папкович в 1932 г. [61], который предложил более удобную форму для представления перемещений: (4.3) где функции B, ϕ определяются из уравнений Аналогичное представление, найденное иным путем, дал Neuber (1934 г.) [62]. Следует также отметить, что в 1928 г. представление, аналогичное представлению Папковича для случая отсутствия объемных сил, получил Г.Д. Гродский. Обобщение на случай пространства произвольной размерности осуществлено Westergaard (1935 г.) [67]. Доказательство полноты представления (4.3) дал П.Ф. Папкович, основываясь на интуитивном предположении о существовании решения исходной задачи. Более строгое доказательство дано Sternberg и Gurtin (1961 г.). Считается, что вывод представления Папковича, основанный на разложении Гельмгольца (1.2), впервые был дан Mindlin (1936 г.) [68]. Миндлин приходит к бигармоническому уравнению Галеркина. Следует отметить, что именно этот путь был указан в статье Папковича ещё в 1932 г. [61]. Представление Папковича также может быть получено из операторного формализма. Вначале намеренно сократим функциональные возможности представления, положив скалярный потенциал ϕ равным нулю. Пусть оператор A, как и прежде, определяется уравнениями Navier-Cauchy (2.4) и однородными краевыми условиями на границе а оператор задается выражением где α - некоторое число, а - оператор, обратный к оператору Лапласа Δ, область определения которого, без потери общности, может быть определена однородными краевыми условиями. Композиция определенных таким образом операторов может быть представлена в виде Если положить, что k удовлетворяет уравнению чему соответствует следующее значение числа α, то полученное выражение существенно упрощается и приводится к виду Таким образом, композиция операторов A, представляет тождественный оператор на множестве функций, которые удовлетворяют условию Дополним конструкцию оператора . Пусть тогда При условиях имеем Таким образом, приходим к представлению (произвольного) поля: Это и есть представление Папковича: 3. Ещё один класс представлений образуют решения Naghdi-Hsu [69]: которые позволяют выразить поле перемещений через вектор-функцию ψ, определяемую уравнением Представление Naghdi-Hsu также может быть получено из операторного формализма. Пусть тогда имеет место следующее соотношение: Если положить, что k удовлетворяет уравнению т.е. параметр α связан с k соотношением то имеем Приходим к представлению Это и есть представление Naghdi-Hsu [69]: Заметим, что, если положить приходим к представлению Галеркина 4. Представление Аржаных-Слободянского [5] соответствует следующему оператору: Имеем (учитывая соотношение ) если Таким образом, представление Слободянского может быть использовано для полей, удовлетворяющих условию 5. Решение статических уравнений Beltrami-Michell (2.18) может быть представлено в форме Shaefer [70], которая в безразмерных переменных принимает вид Здесь функции Ω, Ф определяются уравнениями 6. Существуют представления, приводящие к уравнениям пятого порядка (ультрагармоническим по терминологии [61]). Одно из таких представлений было предложено С.А. Савиным в 1923 г. Представление Савина имеет вид разложения Helmholtz: в котором скалярный потенциал определяется произведением векторного потенциала на постоянный вектор а: Таким образом, перемещения представляются в виде (4.4) Подстановка в уравнение для объемной деформации приводит к следующему выражению: Вычисляя второй ротор левой и правой частей уравнения Navier-Cauchy (2.4), имеем В результате подстановки представления (4.4) Приходим к четырем уравнениям. Одно из уравнений можно исключить. Так как то имеет место равенство Учитывая уравнение для объемной деформации, приходим к системе трех уравнений Эти уравнения могут быть преобразованы к виду В статическом приближении приходим к уравнениям пятого порядка (4.5) В компонентах уравнения (4.5) записываются так: Заметим, что в [61] эта система уравнений написана с многочисленными ошибками. Кроме того, решения этой системы могут содержать лишние относительно исходной системы, что связано с повышением порядка уравнений в процессе их преобразований. Обобщения приведенных выше представлений на задачи для анизотропной упругой среды, а также для реологических соотношений, включающих вторые градиенты, представлены в работах [71-73]. 5. Представления решений статических уравнений линейной теории упругости в криволинейных координатах Представление решений в криволинейных координатах, в частности цилиндрических, позволяет построить ряд новых представлений, не сводящихся к обсуждаемым выше. Вначале перечислим специальные представления решений, не предполагающих какой-либо симметрии, в том числе осевой. 1. Решение однородных уравнений теории упругости, найденные Burgatti в цилиндрических координатах: Компоненты представления определяются тремя гармоническими функциями 2. Обобщенное (неосесимметричное) решение Boussinesq 3. Обобщенное решение Love 4. Расширенное обобщенное решение Love-Muki К решениям, обладающим осевой симметрией, относятся - решение Boussinesq [60] - решение Timpe [74] - решение Love - решение Michell [75] Приведенные решения определяют три компоненты векторного поля через меньшее количество независимых скалярных функций. Возможность такого представления зависит от формы границы тела. Исследованиям полноты подобных представлений посвящены работы [53, 54, 76-78]. 6. Декомпозиция уравнений связанной термоупругости и пороупругости Различные модели пороупругих и многофазных сред определяют линейные дифференциальные уравнения и системы уравнений, структура которых близка к обсуждаемым в настоящей работе. Общие принципы построения таких моделей изложены, например, в [79, 80]. Полные системы уравнений, определяющие совместное деформирование нескольких фаз и соответствующие аналитические представления фундаментальных решений приведены в [81]. Методы редукции системы, позволяющие при дополнительных условиях сократить количество независимых функций и, соответственно, порядок системы уравнений, изложены в [82]. Ниже рассмотрим редуцированные системы, содержащие сокращенное число уравнений. В наиболее простой постановке линейные уравнения связанной термоупругости и пороупругости для изотропного тела могут быть представлены в общем виде: (6.1) Так же как и в уравнениях, приведенных выше, поле u представляет собой перемещения среды: f - объемная плотность внешних сил; ρ - плотность массы; λ, μ - упругие модули Ламе. Физический смысл коэффициентов a1…a8 и полей θ, r зависит от интерпретации теоретической модели. Приведем некоторые из них. 1. В рамках моделей классической связной термоупругости [6] следует положить Здесь αT - коэффициент линейного температурного расширения; κ - коэффициент теплопроводности; c - теплоемкость при постоянной деформации; r - мощность источников тепла; T0 - отсчетная (характерная) температура, θ = T - T0; T - температура в деформированном состоянии. При этом связанные уравнения термоупругого движения и теплопроводности принимают вид 2. В моделях гиперболической теплопроводности (Cattaneo, Jeffreys) [83, 84] - коэффициент, характеризующий период релаксации внутренней энергии в нестационарных процессах; . Связные уравнения принимают вид 3. В бездиссипативной связной термоупругости (Green, Naghdi), которая моделирует эффект второго звука [85], следует положить Здесь - коэффициент, характеризующий волновой характер теплопроводности. Cвязанные уравнения бездиссипативного термоупругого движения и теплопроводности записываются следующим образом: 4. В рамках теории пороупругой среды, заполненной жидкостью (Biot) [82], скалярное поле θ определяет давление жидкости в порах, Уравнения принимают вид 5. В теории пористой среды с пустыми порами (voids) [86] скалярное поле θ определяет пористость (отношение объема пор к объему скелета). Уравнения могут быть записаны в следующей форме: Уравнения (6.1) могут быть представлены в безразмерной форме. Для этого вводятся безразмерные переменные где R - некоторый характерный линейный размер тела (например, диаметр шара, описанного около тела); θ0 - характерное значение скалярного поля (например, характерная температура). Как уже было принято, различие в обозначениях размерных и безразмерных величин (тильда) далее опускается. Безразмерная форма уравнений (6.1) имеет вид (6.2) Здесь Уравнения (6.2) удобно записать в операторной форме (6.3) где 1. Представление общих решений уравнений (6.1) может быть дано в форме Green-Lamé (3.2) [87]: (6.4) Подстановка представлений (6.4) в уравнения (6.3) приводит к следующим уравнениям: Таким образом, функции ϕ, ψ могут быть найдены как решения уравнений Здесь D1, D2 - дифференциальные операторы, задаваемые выражениями 2. Представление решений (6.1) в форме Boussinesq-Papkovitch-Neuber имеют вид (см. также [87]) (6.5) где α - постоянная. Подстановка представлений (6.5) в уравнения (6.3) приводит к уравнениям относительно функции ϕ, ψ: Последнее уравнение может быть преобразовано к виду Это уравнение становится более простым, если положить Имеем 3. Представление решений уравнений (6.1) в форме Cauchy-Kovalevski-Somiglianа могут быть записаны следующим образом (см. также [87]): где функции g, ϕ представляют решение уравнений Основные формулы для представлений сведены в табл. 3. Таблица 3 Представление решений динамической теории термо(поро)упругости № Представление перемещений и температуры Уравнения для элементов представления Вспомогательные уравнения 1 2 3 Отсутствуют 4 Отсутствуют

About the authors

S A Lychev

Moscow State University of Civil Engineering

A D Polyanin

National Research Nuclear University MEPhI (Moscow Engineering Physics Institute)

A L Levitin

Institute for Problems in Mechanics, Russian Academy of Sciences

References

  1. Green A.E. Zerna W. Theoretical Elasticity. - Dover Publications, 2002. - 457 p.
  2. Gurtin M.E. The Linear Theory of Elasticity. Handbuch der Physik (ed. S. Flügge). Vol. VI a/2 // Mechanics of Solids II. Ed. C. Truesdell. - Springer, 1972. - 295 p.
  3. Ильюшин А.А. Механика сплошной среды. - М.: Изд-во МГУ, 1990. - 310 c.
  4. Седов Л.И. Механика сплошной среды. - СПб.: Лань, 2004. - Т. 2. - 560 с.
  5. Лурье А.И. Теория упругости. - М.: Наука, 1970. - 940 с.
  6. Новацкий В. Теория упругости. - М.: Мир, 1975. - 872 c.
  7. Truesdell C., Noll W. The Non-Linear Field Theories of Mechanics. 3rd ed. - Springer, 2004. - 602 p.
  8. Truesdell C., Toupin R.A. The Classical Field Theories. Handbuch der Physik (ed. S. Flügge). Bd. III/1. - Springer, 1960. - P. 226-793.
  9. Stokes G.G. On the dynamical theory of diffraction // Trans. Cambridge Phil. Soc. - 1856. - Vol. 9. - P. 1-62.
  10. Helmholtz Н. Über Integrale der hydrodynamischen Gleichungen, welche den Wirbelbewegungen entsprechen // J. Reine und Angewandte Mathematik - 1858. - Vol. 55. - P. 25-55.
  11. Gibbs J.W. Vector analysis. - Yale Bicentennial Publications, 1901. - 436 p.
  12. Курант Р. Уравнения с частными производными. - М.: Мир, 1964. - 830 с.
  13. Положий Г.Н. Уравнения математической физики. - М.: Высш. шк., 1964. - 560 c.
  14. Arfken G.B., Weber H.J. Mathematical Methods for Physicists. Sixth Edition. - Elsiver Academic Press, 2005. - 1182 p.
  15. Gurtin M.E. On Helmboltz’s theorem and the completeness of the Papkovich-Neuber stress funaions for infinite domains // Arc. Rat. Mech. Anal. - 1962. - Vol. 9. - No. 1. - P. 225-233.
  16. Ладыженская О.А. Смешанная задача для гиперболического уравнения. - М.: Гостехиздат,1953. - 282 с.
  17. Polyanin A.D., Zhurov A.I. Integration of linear and some model nonlinear equations of motion of incompressible fluids // Int. J. Non-Linear Mechanics. - 2011. - Vol. 48. - P. 77-83.
  18. Полянин А.Д., Вязьмин А.В. Декомпозиция трехмерных линеаризованных уравнений вязкоупругих жидкостей Максвелла, Олдройда и их обобщений // Теоретические основы химической технологии. - 2012. - Т. 47, № 4. - С. 386-394.
  19. Липатов И.И., Полянин А.Д. Декомпозиция и точные решения уравнений вязкой слабосжимаемой баротропной жидкости // Доклады академии наук. - 2013. - Т. 449, № 3. - C. 290-294.
  20. Morino L. Helmholtz decomposition revised: Vorticity generation and trailing edge condition // Computational Mechanics. - 1986. - Vol. 1. - P. 65-90.
  21. Воинов О.В. Симметрии медленных движений и равновесий сплошных сред // ПММ, 2007. - Т. 71. - Вып. 5. - С. 788-793.
  22. Backus G.E. A class of self-sustaining dissipative spherical dynamos // Ann. Phys. - 1958. - Vol.4. - P. 372-447.
  23. Kaiser R. A Toroidal Magnetic Field Theorem // Commun. Math. Phys. - 2009. - Vol. 290. - P. 633-649.
  24. B. Sri Padmavati, T. Amaranath. Note on Decomposition of Solenoidal Fields // Applied Mathematics Letters. - 2002. - Vol. 15, 803-805 pp.
  25. Chadwick P., Trowbridge E.A. Elastic wave fields generated by scalar wave functions // Proc. Cambridge Phil. Soc. - 1967. - Vol. 63. - P. 1177-1187.
  26. Кочин Н.Е. Векторное исчисление и начала тензорного исчисления. - М.: Наука, 1965. - 427 с.
  27. Lamé G. Lections sur la théorie mathématique de l'élasticité des corps solides. - Paris: Bachelier, 1866. - 384 p.
  28. Гольдштейн Р.В., Городцов В.А., Лисовенко Д.С. Ауксетическая механика кристаллических материалов // Известия РАН, МТТ. - 2010. - № 4. - C. 43-62.
  29. Navier C.-L.-M.-H. Sur les lois de l’equilibre et du mouvement des corps solides elastiques // Bull. Soc. Philomath. - 1821. - P. 177-181.
  30. Cauchy A.-L. Sur les equations qui expriment les conditions d’equilibre ou les lois du mouvement interieur d’un corps solide, elastique ou non elastique // Ex. de Math. - 1822. - Vol. 3. - P. 160-187.
  31. Capecchi D., Ruta G., Trovalusci P. From classical to Voigt’s molecular models in elasticity // Arch. Hist. Exact. Sci. - 2009. - Vol. 64. - No. 5. - P. 525-558.
  32. Darrigol O. Between Hydrodynamics and Elasticity Theory: The First Five Births of the Navier-Stokes Equation // Arch. Hist. Exact. Sci. - 2002. - Vol. 56. - No. 2. - P. 95-150.
  33. Cauchy A.-L. Exercices d'analyse et de physique mathematique. - Paris: Bachelier, 1840. - 442 p.
  34. Tedone O. Saggio di una teoria generale delle equazioni dell'equilibrio elastico per un corpo isotropo // Annali di Matematica Pura ed Applicata. Pt. 1. - 1903. - Vol. 8. - No. 1. - P. 129-180; Pt. 2. - 1904. - Vol. 9. - No. 1. - P. 13-64; Pt. 3. - 1909. - Vol. 16. - No. 1. P. 141-177.
  35. Kröner E. Allgemeine Kontinuumstheorie der Versetzungen und Eigenspannungen // Arch. Rat. Mech. Anal. - 1959/60 (1967). - No. 4(1). - P. 273-334.
  36. Teodorescu P.P. Stress Functions in Three-Dimensional Elastodynamics // Acta Mechanica. - 1972. - Vol. 13. - P. 103-117.
  37. Beltrami E. Osservazioni sulla Nota precedente // Att. Accad. Nazl. Lincei. - 1888. - No. 1. - P. 141-142.
  38. Donati L. Ulteriori osservazioni intorno al. teorema del Menabrea // Memorie della Accademia delle Scienze dell'Istituto di Bologna. - 1894. - No. 4. - P. 449-474.
  39. Michell J.H. On the direct determination of stress in an elastic solid with application to the theory of plates // Proc. London Math. Soc. - 1899. - No. 31. - P. 100-124.
  40. Kellogg O.D. Foundation of Potential Theory. - Berlin, 1928. - 384 p.
  41. Poisson S.-D. Mémoire sur l'équilibre et le mouvement des corps élastiques // Mém Acad. Sci. Inst France. - 1828. - No. (2) 8. - P. 357-570.
  42. Eringen A.C., Suhubi E.S. Elastodynamics. Vol. I. Finite Motions. Vol. II. Linear Theory. - New York: Academic Press, 1975. - 1003 p.
  43. Kovalevski S. Über die Brechung des Lichtes in cristallinischen Mitteln // Acta Math. - 1885. - Vol. 6. - P. 249-304.
  44. Iacovache M. O extindere a metodei lui Galerkin pentru sistemul ecuaiilor elasticitatii // Bui. St. Acad. R. P. Romane. - 1948. - Vol. 1. - P. 593-596.
  45. Somigliana C. Sulle equazioni della elasticita // Ann. Mat. - 1889. - Vol. (2) 17. - P. 37-64.
  46. Слободянский М.Г. Функции напряжений для пространственной задачи теории упругости // Учен. зап. Моск. гос. ун-та. Вып. XXIV. Механика. - 1938 - Кн. 2. - С. 181-190.
  47. Слободянский М.Г. Об общих и полных формах решений уравнений упругости // ПММ. - 1958 - Т. 23. - Вып. 3. - С. 468-482.
  48. Остросаблин Н.И. Общие решения и приведенные системы уравнений линейной теории упругости к диагональному виду // ПМТФ - 1993 - Т. 33, № 5. - С. 112-122.
  49. Остросаблин Н.И. Операторы симметрии и общие решения уравнений линейной теории упругости // ПМТФ. - 1995. - Т. 36, № 5. - С. 98-104.
  50. Остросаблин Н.И. Упругий анизотропный материал с чисто продольными и поперечными волнами // ПМТФ. - 2003. - Т. 44, № 2. - С. 143-151.
  51. Остросаблин Н.И. Общее решение и приведение системы уравнений линейной изотропной упругости к диагональному виду // Сиб. журн. индустр. матем. - 2009. - Т. 12, № 2. - С. 79-83.
  52. Остросаблин Н.И. Диагонализация системы статических уравнений Ламе линейной изотропной упругости // Сиб. журн. индустр. матем. - 2012. - Т. 15, № 3. - С. 87-98.
  53. Tran Cong T., Steven G.P. On the representation of elastic displacement fields in terms of three harmonic functions // J. of Elasticity. - 1978 - Vol. 8 (3). - P. 325-333.
  54. Hackl K., Zastrow U. On the existence, uniqueness and completeness of displacements and stress functions in linear elasticity // J. of Elasticity. - 1988 - Vol. 17 - P. 3-23.
  55. Chandrasechariah D.S. Naghdi-Hsu Type Solution in Elastodynamics // Acta Mechanica. - 1988. - Vol. 76. - P. 235-241.
  56. Chandrasekharaiah D.S. A complete solution in elastodynamics // Acta Mechanica. - 1990. - Vol. 84. - No. 1-4. - P. 185-190.
  57. Wang M.Z. Note On the Naghdi-Hsu type solution in elastodynamics // Acta Mechanica. - 1991. - Vol. 86. - P. 225-226.
  58. Крутков Ю.А. Тензор функций напряжений и общие решения в статике теории упругости. - М.-Л.: Изд-во АН СССР, 1949. - 199 с.
  59. Деев В.М. О формах общего решения пространственной задачи теории упругости, выраженных при помощи гармонических функций // ПММ. - 1959. - Т. 23. - Вып. 6. - C. 1132-1133.
  60. Boussinesq J. Application des potentiels a l’etude de l’equilibre et des mouvements des solides elastiques. - Paris: Gauthier-Villars, 1885. - 734 p.
  61. Папкович П.Ф. Выражение общего интеграла основных уравнений теории упругости через гармонические функции // Известия Академии наук СССР. - 1932. - № 9. - С. 1425-1435.
  62. Neuber H. Ein neuer Ansatz zur Losung raumlicher Probleme der Elastizitatstheorie. Der Hohlkegel unter Einzellast als Beispiel // Z. Angew. Math. Mech. - 1934. - Vol. 14. - P. 203-211.
  63. Pigeaud G. Resistance des materiaux et Elasticite. - Paris: Gauthier-Villars, 1920. - 482 p.
  64. Benthem J.P. Note on the Boussinesq-Papkovich stress-functions // J. of Elasticity. - 1978 - Vol. 8. - No. 2 - P. 201-206.
  65. Galerkin B. Contribution a la solution ge ne rale du proble me de la the orie de l’e lasticite dans le cas de trios dimensions // Comptes Rendus. - 1930. - Vol. 190. - P. 1047-1048.
  66. Moisil G. Asupra formulelor lui Galerkin in teoria elasticitatii // Bui. St. Acad. R. P. Romane. - 1948. - Vol. 1. - P. 587-592.
  67. Westergaard H.M. General solution of the problem of elastostatics on n-dimensional homogeneous isotropic solid in an n-dimensional space // Bull. Amer. Math. Soc. - 1935. - Vol. 41. - No. 10. - P. 695-698.
  68. Mindlin R.D. Note on the galerkin and papkovitch stress functions // Bull. Amer. Math. Soc. - 1936. - Vol. 42. - No 6. - P. 373-376.
  69. Naghdi P.M., Hsu C.S. On a representation of displacements in linear elasticity in terms of three stress functions // J. Math. Mech. - 1961. - Vol. 9. - P. 233-245.
  70. Shaefer H. Die Spannungsfunctionen des dreidimensionalen Kontinuums und des elastischen Korpers // Z. Angew. Math. Mech. - 1953. - Vol. 33. - P. 356-362.
  71. Bickford W.B. A Generalization of the Boussinesq-Papkovich Functions to an Anisotropic Medium // ZAMP. - 1977. - Vol. 28. - P. 723-728.
  72. Eskandari-Ghadi M. A Complete Solution of the Wave Equations for Transversely Isotropic Media // J. of Elasticity - 2005. - Vol. 81. - No. 1. - P. 1-17.
  73. Goodman W.H., Naghdi P.H. The use of displacement potentials in second order elasticity // J. of Elasticity. - 1988. - Vol. 22. - P. 25-44.
  74. Timpe A. Achsensymmetrische Deformation von Umdrehungskörpern // Z. Angew. Math. Mech. - 1924. - Vol. 4. - No. 5. - P. 361-376.
  75. Michell J.H. The Uniform Torsion and Flexure of Incomplete Tores, with application to Helical Springs // Proc. London Math Soc. - 1900. - Vol. 31. - P. 130-146.
  76. Tran-Cong T. On the completeness and uniqueness of the Papkovich-Neuber and the non-axisymmetric Boussinesq, Love, and Burgatti solutions in general cylindrical coordinates // J. of Elasticity. - 1995. - Vol. 36. - P. 227-255.
  77. Wang M.Z. On the completeness of solutions of Boussinesq, Timpe, Love, and Michell in axisymmetric elasticity // J. of Elasticity. - 1988. - Vol. 17. - P. 85-92.
  78. Carlson D.E. A note on the solutions of Boussinesq, Love, and Marguerre in axisymmetric elasticity // J. of Elasticity. - 1983. - Vol. 12. - P. 345-348.
  79. Нигматулин Р.И. Динамика многофазных сред. - М: Наука, 1987. - Т. 1. - 464 с.; Т. 2. - 360 с.
  80. Николаевский В.Н. Механика пористых и трещиноватых сред. - М.: Недра, 1984. - 323 с.
  81. Burridge R., Vargas C.A. The fundamental solution in dynamic poroelasticity // Geophys. J. R. astr. Soc. - 1978. - Vol. 58. - P. 61-90.
  82. Ковтун А. Об уравнениях модели Био и их модификациях // Вопросы геофизики. - Вып. 44. - СПб., 2010. - С. 3-26.
  83. Cattaneo C. A form of heat conduction equation which eliminates the paradox of instantaneous propagation // Comptes Rendus. - 1958. - Vol. 247. - P. 431-433.
  84. Joseph D.D., Preziosi L. Heat waves // Reviews of Modern Physics. - 1988. - Vol. 61. - P. 41-73.
  85. Green A.E., Naghdi P.M. Thermoelasticity without energy dissipation // J. of Elasticity - 1993. - Vol. 31. - No. 3. - Р. 189-208.
  86. Cowin S.C., Nunziato J.W. Linear elastic materials with voids // J. of Elasticity. - 1983. - Vol. 11 (2). -P. 125-147.
  87. Chandrasechariah D.S., Cowin S.C. Unified complete solutions for the theories of thermoelasticity and poroelasticity // J. of Elasticity - 1988. - Vol. 21. - P. 121-126.

Statistics

Views

Abstract - 213

PDF (Russian) - 105

Cited-By


PlumX


Copyright (c) 2015 Lychev S.A., Polyanin A.D., Levitin A.L.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies