About geometrically nonlinear constitutive relations for elastic material
- Authors: Trusov PV1, Kondratev NS2, Shveykin AI1
- Affiliations:
- Perm National Research Polytechnic University
- Research Institute of Mechanics Lobachevsky Nizhniy Novgorod State University
- Issue: No 3 (2015)
- Pages: 182-200
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/270
- DOI: https://doi.org/10.15593/perm.mech/2015.3.13
- Cite item
Abstract
In the field of solid mechanics it is often necessary to use constitutive (physical) relations in the rate form, for example, for formulation of the boundary value problem in the rate form with contact conditions when a contact area is a priori unknown and changes in a deformation process. The paper considers some questions of formulating geometrically nonlinear constitutive relations of elastic material in the rate form and the relationship of these equations and the constitutive relations in the finite form. In many existing elastic and elasto-plastic models Hooke's law (written in terms of the actual configuration) is used as the constitutive relation. As a rule, the rate measure of the strain is deformation rate tensor (the symmetrical part of velocity gradient) and the stress rate measure is some objective derivative (convective or corotational) of weighted Kirchhoff stress tensor. The use of the convective derivative leads to certain difficulties, for example, analysis of evolution stress state in a deformable basis is complicated; consequently, convective derivatives are excluded from consideration. Using a stress tensor corotational derivative instead of the material derivative (derivative by time) allows satisfying the principle of material indifference (the principle of independence of constitutive relation from the choice of reference frame), however a choice of the derivative type can be implemented by many ways. An arbitrary selection of stress corotational derivative leads to undesirable effects: stress oscillations for the simple shear monotonic loading (for example, for Jaumann derivative), “not closed” stress trajectories and nonzero stress work for a closed deformation trajectory. In papers of A. Meyers, H. Xiao, O. Bruhns corotational derivative (with logarithmic spin or logspin) was proposed.The derivative of the right Hencky's strain tensor is exactly equal to deformation rate tensor. When using this derivative in the constitutive relation the described effects are missing, on the basis of this remarked logspin authors state the logspin exclusivity and recommend it only for using in rate form constitutive relations. The article shows full compliance between Hooke's law in the rate and finite form under conditions: existence of the material basis in which properties of the body (in this case - elastic) remain unchanged and the use of the same type corotational derivative for stress and strain measures. Computational experiments for illustration of the proved assertion about the equivalence of different Hooke's law forms were conducted: various measures of strain state and their corotational derivatives are considered; a kinematic loading in a closed cycle (in the space of deformations) is applied. It is shown for isotropic elastic materials that stress trajectories are closed. When using work-conjugacy stress and strain measures the dissipation energy is absent. Thus it is possible to question the exclusive choice of logarithmic spin and the above-mentioned commonly used measures of stress-strain state. In formulating and solving the problems of solid mechanics there should be an opportunity to use different stress-strain measures and their objective derivatives, the choice of measures and constitutive relations should be justified from the point of physical analysis of the process.
Full Text
Введение и основные понятия Проблема построения определяющих соотношений для случая геометрической нелинейности в последние десятилетия является одним из весьма значимых направлений развития механики деформируемого твердого тела (МДТТ). Это вызвано как необходимостью дальнейшего развития теории механики деформируемого твердого тела, так и потребностями практики: - в промышленности постоянно повышаются требования к качеству материалов и готовой продукции, следовательно, к качеству математических моделей технологических процессов, для корректного описания многих из которых них необходим учет геометрической нелинейности. Обзоры литературы по данной проблеме можно найти в монографиях [1-3]. По мнению авторов, при построении определяющих соотношений в терминах отсчетной конфигурации (такие соотношения описываются, например, в работах [3, 4]) возникают некоторые проблемы с физическим анализом как самих соотношений, так и входящих в них мер деформированного и особенно напряженного состояния. При построении определяющих соотношений в терминах актуальной конфигурации эти вопросы частично снимаются, по крайней мере, по отношению к мере напряженного состояния: тензор напряжений имеет ясно выраженный физический смысл. Для ряда задач МДТТ имеется насущная потребность в формулировки формулировке определяющих соотношений в терминах актуальной конфигурации в скоростях. Например, это необходимо при постановке контактных краевых задач, когда изменяющаяся область контакта определяется в процессе решения задачи [3]. Можно также отметить, что постановки краевых задач в скоростной форме, более удобны для применения численных методов, например, метода конечных элементов [5, 6], - когда решение задачи осуществляется пошагово, в приращениях [3]. Поэтому в последние десятилетия весьма популярными являются определяющие соотношения (ОС) в скоростной форме, сформулированные с использованием мер напряженно-деформированного состояния (НДС), определёенных в текущей конфигурации. При построении таких соотношений, помимо вопроса выбора мер НДС (скоростей мер НДС), актуальной является проблема выбора не зависящей от системы отсчета производной [3, 7-9]. Действительно, в отличие от случая использования в ОС материальной производной, применение любой не зависящей от выбора системы отсчета производной позволяет удовлетворить принципу материальной индифферентности (принципу независимости определяющего соотношения от выбора системы отсчета [10, 11, 3]), однако различный выбор вида производной приводит к существенным (зачастую - качественным) отличиям получаемых НДС. В общем виде не зависящая от выбора системы отсчета производная тензора второго ранга определяется выражением , где - соответственно векторы основного и сопряженного базисов подвижной системы координат (СК). Образно говоря, подвижный наблюдатель в каждый момент времени фиксирует мгновенные изменения компонент тензора и относит их к своему базису, считая его неизменным. Таким образом, использование независимых от выбора системы отсчета производных подразумевает разделении разделение движения на две составляющие: квазитвердое движение подвижной СК с базисом и собственно деформационное движение (относительно этой подвижной СК), мера скорости которого связывается в ОС со скоростью изменения напряженного состояния. Очевидно, что выбор подвижной СК и, соответственно, коротационной производной можно осуществить множеством способов. Наиболее известные и часто используемые в работах по механике сплошной среды не зависимые зависящие от выбора системы отсчета производные рассмотрены в [4, 12-14, 3 и др.]. Базисы для конвективных производных являются деформируемыми (например, для производной Олдройда [12] и Коттер - Ривлина [13] подвижный базис - текущий лагранжев), что, по мнению авторов, приводит к неясности в определении тензора напряжений и невозможности анализа изменения длин материальных отрезков и углов между ними (которые наблюдаются из неподвижной лабораторной системы отсчета, но не «ощущаются» подвижным наблюдателем), равно как к сложностям отделения изменения компонент мер напряженного состояния, определяемых процессом деформирования, от их изменения из-за деформирования базисных векторов. Поэтому далее рассматриваются только коротационные производные (с ортонормированным базисом подвижной СК). Известно, что произвольный выбор коротационной производной меры напряжений приводит к нежелательным эффектам: осцилляциям напряжений при монотонной деформации простого сдвига (например, для производной Зарембы-Яуманна [3]), «незамкнутости» траекторий напряжений и отличию от нуля работы напряжений по замкнутой траектории упругих деформаций. В работах [15, 16] предложена коротационная производная (ее спин назван логарифмическим), обладающая следующим свойством: эта производная от правого тензора деформации Генки в точности равна тензору деформации скорости. На основании того, что при использовании этой производной в определяющем соотношении упругого (важно обратить внимание) изотропного материала описанные эффекты отсутствуют (соответствующие результаты приводятся, например, в работе [17]), авторы логарифмического спина рекомендуют только его к использованию в определяющих соотношениях в скоростной форме, причем не только для упругого материала - в работе [18] логарифмический спин предлагается использовать в модели упругопластичности. Отметим, что построение геометрически нелинейных ОС особенно затруднительно для физически нелинейных материалов, например, упругопластических и упруговязкопластических [1-3]. Следует отметить, что при значительных деформациях в поликристаллических металлах и сплавах происходит процесс текстурообразования, что приводит к появлению анизотропии свойств, в связи с чем обоснованность применения логарифмического спина в соотношениях макроуровня вызывает определенные сомнения. В последние десятилетия при формулировке геометрически и физически нелинейных соотношений все большую популярность приобретает многоуровневый подход, основанный на введении внутренних переменных [19-22]. В моделях авторов предлагаемой статьи, построенных в рамках этого подхода с явным описанием внутренней структуры и механизмов деформирования, формулируются соотношения на низшем масштабном уровне (на уровне кристаллитов используется закон Гука в скоростной форме [23]), соотношения макромасштабного уровня получаются путем проведения процедуры согласования ОС масштабных уровней [23]. Таким образом, при применении многоуровневого моделирования одной из важных задач является корректное описание ОС на низших масштабных уровнях, в частности, для двухуровневых моделей материалов [23] - корректная формулировка геометрически нелинейного закона Гука. В предлагаемой статье устанавливается связь между законом Гука для анизотропного материала в скоростной и конечных формах (раздел 2) при условиях существования материального базиса, в котором компоненты тензора упругих свойств тела остаются неизменными, и использования однотипных коротационных производных от мер деформаций и напряжений. В разделе 3Также рассмотрены некоторые меры деформаций и их не зависимые зависящие от выбора системы отсчета скорости изменения, в том числе - получаемые с применением логарифмического спина. Приведены Ррезультаты вычислительных экспериментов для замкнутых циклов деформаций, приведены в разделе 4, которые показывают, что при соответствующем выборе мер НДС и мер их изменения для изотропных упругих материалов при любом выборе коротационной производной наблюдается замкнутость траектории напряжений; если при этом на меры НДС накладывается условие энергетической сопряженности, то наблюдается и отсутствие диссипации энергии. Полученные результаты свидетельствует свидетельствуют о недостаточной аргументации вышеприведенного утверждения об исключительной необходимости использования логарифмического спина и соответствующих скоростей мер НДС. По мнению авторов, при постановке и решении задач механики твердого тела должна иметься возможность использования различных мер НДС и их скоростей, при этом выбор мер и ОС должен быть обоснован с позиций физического анализа рассматриваемого процесса. 1. Меры деформации и их не зависимые зависящие от выбора системы отсчета скорости изменения В настоящее время в нелинейной механике используется широкий спектр мер деформированного состояния. Распространение получило полярное разложение градиента места, с использованием которого вводится семейство мер Сетха [4]. В силу неособенности градиента места его можно представить в виде скалярного произведения симметричного и ортогонального тензоров: , (1) , где - набла-оператор, определенный в отсчетной конфигурации K0, ; r - радиус-вектор материальной частицы в актуальной конфигурации Kt,; U, V - соответственно левый и правый тензоры искажения, ; R - ортогональный тензор, сопровождающий деформацию. Поскольку U и V симметричные положительно определенные тензоры, то для них существует тройка положительных действительных собственных чисел λi (кратность удлинения) и тройки ортогональных векторов тензора U и тензора V. В главном базисе U и V тензоры, входящие в (1), примут вид: (2) , С использованием введенных тензоров U и V можно ввести широкий класс мер деформаций и [4, 3]: (3) , где - монотонно гладкие функции удовлетворяющие условиям: . (4) , Введенное семейство мер включает в себя семейство мер деформаций Сетха: , (5) , здесь где I - единичный тензор. Определяющие соотношения удобно записать в актуальной конфигурации, где компоненты мер напряжений и деформаций имеют ясный физический и геометрический смысл. Вследствие этого для примера будут рассмотрены случаи, когда в качестве мер деформаций используются тензоры Альманси A и правый Генки , а меры напряжений - тензор Кирхгоффа κ. Все компоненты упомянутые упомянутых тензоры тензоров имеют явный геометрический и физический смысл компонент. Следует отметить, что в отличие от пары тензоров κ и , пара κ и A не являются является энергетически сопряженной. Тензоры A и являются тензорами из семейства Сетха при n = -1 и n = 0, соответственно: (6) , где - набла-оператор в актуальной конфигурации Kt, ; - мера Альманси. Используя диадное представление градиента места ( , - взаимный базис лагранжевой системы координат в отсчетной и основной в актуальной конфигурации), можно убедится, что , R0 - радиус-вектор частицы в K0. Правый тензор Генки получается из семейства Сетха при n = 0. При раскрытии неопределенности типа ноль на ноль используется правило Лопиталя: (7) , Для определения производной по параметру, входящему в степень, правый тензор искажения V записывается в главном базисе с собственными числами λi, далее, учитывая независимость базисных векторов от степени, определялась производная компонент тензора V2n в этом базисе. В общем случае независимые от системы отсчета производные тензоров 2-го ранга (конвективные и коротационные), описывающие некоторые характеристики исследуемой частицы с ее ε-окрестностью, определяются соотношениями [3]: (8) , где qi, qi - основной и сопряженный базисы подвижной системы координат (ПСК). Множественность способов выбора этого базиса определяет многообразие объективных производных. Отметим, что в настоящей работе рассматривается жесткая ортогональная система координат, связанная с материалом, поэтому различие между векторами основного и сопряженного базиса пропадает и соотношения (8)1 и (8)2 являются эквивалентными. Вопрос о выборе ПСК и определении скорости ее вращения ω является актуальной задачей, возможный подход к ее решению предложен в [24]. Не конкретизируя способ определения тензора ω ( ), но предполагая его существование, запишем в общий общем виде коротационной коротационную производной производную индифферентной тензорзначной (второго ранга) функции времени B: , (9) , Далее рассматриваются три типа коротационных производных: 1). Зарембы-Яуманна [25, 26], ассоциированную ассоциированная с тензором вихря : , (10) , 2. ) Грина-Нагхди [27], ассоциированную ассоциированная со спином Ω над ортогональным тензором R из полярного разложения ( ),: , (11) , 3). определяемую определяемая с помощью логарифмического спина Ωlog [15, 16]:: . (12) , Первые две из перечисленных производных получили широкое распространение в механике сплошных сред, и способ их введения не вызывает вопросов, в отличие от последней производной. Вследствие этого кратко остановимся на выводе соотношения для определения логарифмического спина, используя теорию, предложенную в [28, 29]. Исходной задачей является определение тензора логарифмического спина Ωlog, такого, что коротационная производная правого тензора Генки , ассоциированная с этим спином, дает в точности тензор деформации скорости D: . (13) , Из соотношения (13) можно определить искомый тензор Ωlog, для этого (13) записывается в базисе правого трехгранника : , (14) , где - компоненты тензора скорости вращения базисных векторов относительно ЛСК в базисе этих же векторов; Dij, - компоненты тензора деформации скорости и правого тензора Генки в базисе . Поскольку тензор Ωlog является кососимметричным, то в любом базисе имеет шесть ненулевых, из которых только три независимые ( ). Искомые недиагональные компоненты ( ) определяются из (14): , (15) , или при неравенстве главных значений можно записать . (16) , Следует заметить, что на этом можно было бы остановиться и при определении логарифмического спина пользоваться (16). Однако при решении практических задач удобнее работать с компонентами тензоров деформации скорости и вихря (чем, например, с компонентами спина ΩV); компоненты тензора ΩV через компоненты W и D в базисе можно определить следующими соотношениями [3]: , (17) , где Wij - компоненты тензора вихря в базисе . С учетом (28) имеем: , (18) , или , (19) , где - гиперболический тангенс. Выше использована связь компонент тензора Генки с главными значениями тензора V: . Физический смысл компонент логарифмического спина, в отличие от смысла компонент тензора вихря и «материального» спина , непрозрачен. К сожалению, авторы цитируемых статей его не обсуждают. 2. Связь определяющих соотношений закона Гука в скоростной и конечной форме Покажем, что для класса материалов с подвижным базисом, в котором компоненты тензора упругих характеристик материала не меняются (для изотропных материалов это выполняется в любой системе отсчета) в процессе деформирования, формы закона Гука в скоростном и конечном виде полностью эквивалентны. Это, в свою очередь, означает, что в скоростном виде для «замкнутости» траекторий напряжений при замкнутых траекториях деформации достаточно взять от меры деформаций и напряжений одинаковые коротационные производные. Вследствие этого появляются определенные сомнения по поводу предпочтительности выбора в качестве скоростных мер тензора деформации скорости и логарифмической коротационной производной тензора напряжений Коши. Помимо сложности вычисления логарифмического спина, существует более значимый недостаток этой производной - отсутствие явной связи с материалом. Рассмотрим закон Гука в конечной форме для монокристаллического материала в следующем виде [10]: ,. (20) где κ - взвешенный тензор напряжений Кирхгоффа;, ζ - некоторый тензор деформации, определенный в актуальной конфигурации, ; п - тензор упругих свойств монокристалла. Несколько «забегая вперед», отметим, что тензор напряжений Кирхгоффа использован здесь скорее как дань традиции, сложившейся при формулировке ОС гипер- и/или гипоупругих материалов. В общем случае в качестве меры напряжений может быть введен некоторый индифферентный тензор напряжений, определенный в базисе актуальной конфигурации; в частности, на определение этого тензора могут быть наложены дополнительные ограничения, например энергетическая сопряженность с коротационной производной, используемой в скоростной форме ОС. Последнее требование вводится при использовании термодинамического подхода к построению определяющих соотношений; в настоящей работе, ориентированной на применение макрофеноменологического подхода, указанное ограничение не вводится. Найдем полную (материальную) производную от левой и правой части соотношения (20), при этом в общем случае тензор упругих свойств зависит от времени, п(t) (даже если его компоненты в том или ином подвижном базисе являются константами): . (21) Введем две системы отсчета: лабораторную условно неподвижную систему координат (ЛСК) с ортонормированным базисом и жесткая жесткую подвижная подвижную система систему координат (ПСК) с ортонормированным базисом ( ); здесь и далее черта сверху обозначает величины, определенные в ПСК. Без ограничения общности можно считать, что обе системы в начальный момент деформирования совпадают. Полагается, что ПСК связана с материалом, поэтому компоненты тензора упругих свойств в этом базисе остаются неизменными. Связь векторов базиса ПСК и ЛСК описывается ортогональным тензором o(t): , откуда получаем: , - спин ПСК относительно ЛСК. Рассмотрим первое слагаемое правой части (21): (22) Для преобразования двух последних слагаемых следует напомнить, что используется ортонормированный базис, т.е. справедливо ( - дельта Кронекера). Дифференцируя последнее соотношение, имеем: , или . Тогда (22) примет вид: Подставляя полученное выражение в (21), с учетом (20) получаем: , или . (23) Тогда окончательно соотношение запишется в скоростной форме (в ряде источников определяющее соотношение в этой форме называется гипоупругим законом): , (24) где верхний индекс «cr» обозначает ту или иную коротационную производную соответствующего тензора. Далее проведем обратную процедуру - проинтегрируем в ПСК соотношение закона Гука в скоростной форме (24) и покажем его соответствие соотношению в конечной форме (20). Для этого рассмотрим коротационную производную тензорной величины a, определенной в базисе подвижной системы координат , которая вращается со скоростью ω: (25) Полученное соотношение непосредственно следует из определения коротационной производной [3]. Тогда соотношение (24) в базисе ПСК имеет вид: , (26) , или . , Из последнего выражения следует равенство, справедливое для любого момента времени: , (27) , Интегрируя (27), получим: , , откуда с учетом неизменности базиса ЛСК в отсчетной конфигурации и компонент тензора упругих модулей в этой системе координат следует: , (28) , или, полагая отсчетную конфигурацию естественной (свободной от напряжений и деформаций), . (29) , Умножая последнее соотношение слева на o, а справа на oT, имеем: , (30) , что совпадает с соотношением в конечной форме (20). Таким образом, доказана эквивалентность упругих соотношений в конечной и скоростной формах (20) и (24) при выполнении следующих условий: компоненты тензора упругих характеристик в базисе ПСК остаются постоянными, коротационные производные мер напряжений и деформаций идентичны и определены с использованием спина, характеризующего скорость ротации введенной ПСК. 3. Результаты численного моделирования Следует отметить, что рассуждения, проведенные для доказательства эквивалентности законов Гука в скоростной и конечной форме, справедливы для анизотропного случая, однако ниже в качестве примера будет рассмотрен изотропный материал (развитие работы для анизотропных материалов будет отражено в последующих публикациях авторов). В этом случае компоненты тензора п будут неизменны в любой системе координат. Тензор упругих свойств может быть представлен с использованием линейной комбинации трех изотропных тензоров четвертого ранга: , (31) , где Тогда определяющее соотношение в скоростной форме имеет вид: , (32) , с применением свойств изотропных тензоров получается: , (33) , где - первый инвариант. Поскольку рассматриваемые меры деформации и их скорости являются симметричными, последнее соотношение можно переписать в виде: , (34) , где введено обозначение . Соотношение Гука в конечной форме для изотропного материала будет иметь вид: , (35) , Далее в качестве мер деформационного состояния выступают тензор Альманси (6) и правый тензор Генки (7), а используемые коротационные производные в (34) - производные Зарембы-Яуманна (10), Грина-Нагхди (11) и логарифмического спина (12). Важной характеристикой процесса деформирования является работа напряжений, которая при деформировании по замкнутому циклу для упругого материала должна быть в точности равна нулю. Известно, что локальная скорость изменения внутренней энергии u определяется только контактными взаимодействиями, является величиной, независимой от выбора системы отсчета, и равна определяется следующим образом [3]: , (36) где σ - тензор напряжений Коши. Последнее соотношение получается из закона сохранения механической энергии и уравнения движения. В случае симметричного тензора напряжений соотношение (36) можно переписать в виде: (37) , Плотность на единицу массы внутренней энергии u определяется соотношением: (38) , где - плотность в отсчетной и текущей конфигурациях, ; J -якобиан. В вычислительных экспериментах задавалось движение по замкнутой траектории, а отклик материала определялся двумя способами: с использованием (34) и (35). Компоненты тензора напряжений, приведенные ниже, определены в неподвижной лабораторной системе координат (ЛСК). Согласно приведенным выше аналитическим рассуждениям, напряженное состояние не должно быть отличным при применении различных форм упругого определяющего соотношения; другим вопросом является установление отсутствия диссипации энергии - в конце цикла нагружения значение внутренней энергии (38) должно равняться нулю. Для решения тестовых задач использованы два лагранжевых закона движения по замкнутой траектории, предложенные в работе [17]: , (39) (40) где , u - вектор перемещений, ; ki - базис неподвижной ЛСК; qi - лагранжевы координаты, фиксированные в отсчетной конфигурации, ; φ - угол (параметр, по которому ведется кинематическое нагружение), ); h - длина ребра материального прямоугольного в отсчетной конфигурации параллелепипеда с квадратом в поперечном сечении (на рис. . 1 показана проекция параллелепипеда на плоскость поперечного сечения Ox1x2), rd - второй параметр, условный радиус «окружности», по траектории которой частица совершает движение. а б Рис. 1. Траектория движения материальной точки в плоскости Ox1x2 с лагранжевыми координатами q1 = 0, q2 = h, при rd/h = 0,.3 , , (а) - лагранжев закон движения (39),; (б) - лагранжев закон движения (40) В расчетах использовались параметры материала, значения которых близки к характерным для некоторых резин: λ = 15 МПа, μ = 2 МПа, плотность в отсчетной конфигурации равнялась 1200 1200 кг/м3; при интегрировании скоростных соотношений использовался метод Эйлера с шагом интегрирования 10--3. На рис. . 2 показана эволюция компоненты κ12 тензора напряжений Кирхгоффа в лабораторной системе координат для нагружения (39) и (40) соответственно. Остальные компоненты не приводятся, поскольку также совпадают. Следует отметить, что на одном графике накладываются друг на друга четыре кривых, три из них - это зависимости, полученные при интегрировании в ЛСК упругого соотношения в скоростной форме и использовании коротационных производных Зарембы--Яумана, Грина--Нагхди и логарифмического спина, еще одна зависимость - результат расчетов с использованием закона Гука в конечной форме. В конце цикла значения компонент тензора напряжений обращаются в нуль с точностью до вычислительной погрешности, что является прямым следствием эквивалентности упругого закона в форме (24) определяющему соотношению в форме (20). а б б б в г Рис. 2. Эволюция компоненты κ12 тензора напряжений в лабораторной системе координат от параметра нагружения φ при численном интегрирования закона Гука в скоростной форме с мерой деформации тензора Альманси (а, в) и правым тензором Генки (б, г) при использовании производных Зарембы--Яумана, Грина--Нагхди, логарифмического спина, а также закона Гука в конечной форме; (а, б) - закон движения (39), (в, г) - закон движения (40) Анализируя приведенные на рис. 2 зависимости, можно отметить близость компонент тензора напряжений (при использовании ОС с мерой Альманси, - с одной стороны, с мерой Генки - с другой стороны) в области малых градиентов перемещений и их отличие при значительных градиентах перемещений. На рис. . 3 представлена зависимость массовой плотности внутренней энергии u от параметра нагружения φ при разных мерах деформирования - Альманси и Генки. а б Рис. 3. Эволюция массовой плотности внутренней энергии u от параметра нагружения φ:; сплошная линия - мера деформаций тензора Альманси, штриховая - Генки;, (а) - лагранжев закон движения (39);, (б) - лагранжев закон движения (40) Плотность внутренней энергии u определяется согласно формуле (38). Эквивалентность скоростной и конечной формы закона Гука позволяет сразу указать пару энергетически сопряженных тензоров - правый тензор Генки (его коротационная производная с логарифмическим спином равна тензору деформации скорости) и тензор Кирхгоффа. Для указанной пары тензоров при упругом деформировании работа в конце замкнутого цикла должна быть в точности нулем, что подтверждается проведенными расчетами, представленными на рис. 3. Тензоры Альманси и Кирхгоффа не являются энергетически сопряженными, поэтому внутренняя энергия, определяемая по соотношениям (37)-(38) для закона (24) (в изотропном случае - (34)), где ζ - тензор деформации Альманси, в конце цикла далека от нулевого значения (рис. 3, б). Ключевой результат, подтверждающий аналитические выводы о соответствии соотношений в конечной и скоростной форме для изотропных упругих материалов: использование любых (численно подтверждено для 3 вариантов) коротационных производных, одинаковых и для меры напряжений, и для меры деформаций, приводит к одинаковым результатам. Если при этом используются энергетически сопряженные меры НДС, то и диссипация энергии отсутствует. Заключение При построении геометрически нелинейных определяющих соотношений в скоростной форме важным является вопрос выбора мер напряженно-деформированного состояния и их не зависящих от выбора системы отсчета производных. Использование вместо материальной производной (по времени) тензора напряжений его коротационной производной и тензора деформации скорости как меры скорости изменения деформированного состояния позволяет удовлетворить принципу независимости определяющего соотношения от выбора системы отсчета, однако конкретизация вида производной может быть осуществлена множеством способов. Произвольный выбор коротационной производной меры напряжений в этом случае приводит к нежелательным эффектам: осцилляциям напряжений при монотонной деформации простого сдвига (например, для производной Зарембы-Яуманна), «незамкнутости» траектории напряжений и отличию от нуля работы напряжений при замкнутой траектории деформаций. Для упругого изотропного материала в работах [15-17] показано, что использование основанной на логарифмическом спине коротационной производной позволяет избежать этих эффектов, на основании чего авторы цитируемых работ заключают, что именно логарифмический спин должен использоваться в скоростных определяющих соотношениях. В работе установлена связь между законом Гука для анизотропного материала в скоростной и конечной формах при условиях существования подвижного материального базиса, в котором упругие характеристики тела остаются неизменными, и использования однотипных коротационных производных от мер деформаций и напряжений. Для иллюстрации проведены численные эксперименты с использованием некоторых мер деформации и их не зависящих от выбора системы отсчета производных, в том числе получаемых с применением логарифмического спина. При рассмотрении замкнутых циклов деформаций показано, что при одинаковом выборе коротационных производных для мер деформаций и напряжений в определяющем соотношении в скоростной форме, для изотропных упругих материалов при любом выборе коротационной производной выполняется замкнутость траектории напряжений, к тому же, если используемые меры напряжений и деформаций энергетически сопряжены, то наблюдается и отсутствие диссипации энергии. Таким образом, полученные результаты свидетельствует о неоднозначности утверждения об исключительной необходимости использования логарифмического спина при формулировке геометрически нелинейных определяющих соотношений в скоростной форме. По мнению авторов, при постановке и решении задач механики твердого тела выбор мер НДС и коротационных производных должен быть обоснован с позиций физического анализа рассматриваемого процесса.About the authors
P V Trusov
Perm National Research Polytechnic University
N S Kondratev
Research Institute of Mechanics Lobachevsky Nizhniy Novgorod State University
A I Shveykin
Perm National Research Polytechnic University
References
- Коробейников С.Н. Нелинейное деформирование твердых тел. - Новосибирск: Изд-во СО РАН, 2000. - 262 с.
- Левитас В.И. Большие упругопластические деформации материалов при высоком давлении. - Киев: Наукова думка, 1987. - 232 с.
- Поздеев А.А., Трусов П.В., Няшин Ю.И. Большие упругопластические деформации: теория, алгоритмы, приложения. - М.: Наука, 1986. - 232 с.
- Лурье А.И. Нелинейная теория упругости. - М.: Наука, 1980. - 512 с.
- Зенкевич О. Метод конечных элементов в технике. - М.: Мир, 1975. - 542 с.
- Зенкевич О., Морган К. Конечные элементы и аппроксимация. - М.: Мир, 1986. - 318 с.
- Новокшанов Р.С., Роговой А.А. Эволюционные определяющие соотношения для конечных вязкоупругих деформаций // Известия РАН. Механика твердого тела. - 2005. - № 4. - С. 122-140.
- Rogovoy A.A. Formalized approach to construction of the state equations for complex media under finite deformations // Continuum Mechanics and Thermodynamics. - 2012. - Vol. 24. - Р. 81-114. doi: 10.1007/s00161-011-0220-y
- Маркин А.А., Толоконников Л.А. Меры и определяющие соотношения конечного упругопластического деформирования // Прикладные проблемы прочности и пластичности: всесоюзн. межвуз. сб. - Горький: Изд-во Горьк. ун-та, 1987. - С. 32-37.
- Truesdell C., Noll W. The non-linear field theories of mechanics. - Berlin: Springer, 1965. - 602 p.
- Трусделл К. Первоначальный курс рациональной механики сплошных сред. - М.: Мир, 1975. - 592 с.
- Oldroid J.G. On the formulation of reological equations of state // Proc. Roy. Soc. London A. - 1950. - Vol. 200. - Р. 523-541.
- Cotter B.A., Rivlin R.S. Tensors associated with time-dependent stress// Quart. Appl. Math. - 1955. - Vol. 13. - No. 2. - Р. 177-188.
- Прагер В. Введение в механику сплошных сред. - М.: Изд-во иностр. лит., 1963. -312 с.
- Xiao H., Bruhns O.T., Meyers A. Hypo-elasticity model based upon the logarithmic stress rate // J. Elasticity. - 1997. - Vol. 47. - P. 51-68.
- Xiao H., Bruhns O.T., Meyers A. Logarithmic strain, logarithmic spin and logarithmic rate // Acta Mechanica. - 1997. - Vol. 124. - P. 89-105.
- Meyers A., Xiao H. and Bruhns O. Elastic stress ratchetting and corotational stress rates// Technische mechanik. - 2003. - Vol. 23. - P. 92-102.
- Bruhns O.T., Xiao H., Meyers A. Large simple shear and torsion problems in kinematic hardening elasto-plasticity with logarithmic rate// Int. J. of Solids and Structures. - 2001. - Vol. 38. - P. 8701-8722.
- McDowell D. L. A perspective on trends in multiscale plasticity // Int. J. Plasticity. - 2010. - Vol. 26. - Р. 1280-1309. doi: 10.1016/j.ijplas.2010. 02.008
- Trusov P.V., Ashikhmin V.N., Volegov P.S., Shveykin A.I. Constitutive relations and their application to the description of microstructure evolution // Physical Mesomechanics. - 2010. - Vol. 13. - Iss. 1-2. - P. 38-46.
- Trusov P.V., Shveykin A.I. Multilevel crystal plasticity models of single- and polycrystals. Statistical models // Physical Mesomechanics. - 2013. - Vol. 16. - No. 1. - P. 23-33.
- Trusov P.V., Shveykin A.I. Multilevel crystal plasticity models of single- and polycrystals. Direct models // Physical Mesomechanics. - 2013. - Vol. 16. - No. 2. - P. 99-124.
- Trusov P.V., Shveykin A.I., Nechaeva E.S., Volegov P.S. Multilevel models of inelastic deformation of materials and their application for description of internal structure evolution // Physical Mesomechanics. - 2012. -Vol. 15. - Iss. 3-4. - P. 155-175.
- Трусов П.В., Нечаева Е.С., Швейкин А.И. Применение несимметричных мер напряженного и деформированного состояния при построении конститутивных моделей материалов // Физическая мезомеханика. - 2013. - Т. 16, № 2. - С. 15-31.
- Zaremba S. Sur une forme perfectionnée de la théorie de la relaxation // Bull. Int. Acad. Sci. Cracovie. - 1903. - Р. 595-614.
- Jaumann G. Geschlossenes System physikalischer und chemischer Differential-gesetze // Sitzber. Akad. Wiss. Wien, Abt. IIa. - 1911. - В. 120. - S. 385-530.
- Green A. E., Naghdi P. M. A general theory of an elasto-plastic continuum // Arch. Rat. Mech. Anal. - 1965. - Vol. 18. - Р. 251-281.
- Reinhardt W.D., Dubey R.N. Eulerian strain-rate as a rate of logarithmic strain // Mechanics Research Communications. - 1995. -Vol. 22. - P. 165-170.
- Reinhardt W.D., Dubey R.N. Coordinate-independent representation of spins in continuum mechanics // Journal of Elasticity. - 1996. - Vol. 42. - P. 133-144.