CC BY
26
ЕСТЕСТВЕННЫЕ И ТОЧНЫЕ НАУКИ
УДК 539.3
ОБ ОБЪЕКТИВНОСТИ УРАВНЕНИЙ В НЕЛИНЕЙНЫХ ТЕОРИЯХ МЕХАНИКИ ДЕФОРМИРУЕМОГО ТВЕРДОГО ТЕЛА
Канд. техн. наук ШВЕД О. Л.
Объединенный институт проблем информатики НАН Беларуси
Тензорно-инвариантная запись физических уравнений среды гарантирует их независимость от системы координат. Для обеспечения объективности описания процесса (независимости от расположения наблюдателя, связанного с системой отсчета) вводится принцип материальной индифферентности (объективности), т. е. независимости материальных характеристик от выбора системы отсчета [1-3]. Его значение велико, так как «отказ от принципа материальной объективности означает отказ от признания существования объективных закономерностей природных явлений» [2]. Кроме того, принцип является важным инструментом при выборе структуры определяющих уравнений. Соответствующие примеры приведены в [1-3]. Принцип появился в механике несколько десятилетий назад, но на интуитивном уровне использовался с давних времен. Формально он совпадает с принципом инвариантности материальных характеристик при наложении жестких движений за одним исключением, на которое указано в [2, с. 249]. Определитель ортогонального тензора жестких движений не может быть равен -1, так как при замене системы отсчета, как показано в этой работе, ее ориентация не участвует в рассмотрении. В связи с этим следует уточнить формулировку в [3, с. 86].
При жестких движениях актуальной конфигурации входящие в уравнение скаляры, векторы, тензоры преобразуются. Будем обозначать преобразованные величины верхним индексом «звездочка». Тем самым изменяется и уравнение, в которое они входят. Если преобразованное уравнение оказывается эквивалентным исходному, то последнее называется объектив-
ным, т. е. удовлетворяющим принципу материальной объективности, либо необъективным в противном случае. В теориях, учитывающих нелинейность изменения геометрии деформируемого тела, проверка объективности уравнений обязательна, необъективные уравнения использовать нельзя.
Рассмотрим, например, уравнение вязко-упругой среды Максвелла
АТ + Т -|D = 0, (1)
где T, D - индифферентные тензоры напряжений Коши и деформации скорости; А, ц - индифферентные скаляры и точка над символом вверху или справа означает его материальную производную. Обозначим Ot - собственно ортогональный тензор жесткого движения, наложенного на актуальную конфигурацию, и его
спин кососимметричный тензор ю = OTt • Ot =
= -OT • Ot, где индекс T означает операцию транспонирования; точка - скалярное произведение.
Покажем, что уравнение (1) не удовлетворяет принципу материальной индифферентности. Учитывая определения индифферентных скаляра и тензора [1], имеем:
А* = А, ц* = ц, D* = OT • D• Ot, T* = OT • TOt, (T - |D)* = T* - |*D* = OT (T - |D) • Ot,
(T)* = (T*)D = (OT T • Ot)D = OT T-Ot+ (2)
+ ©-oTT-Ot-Ot • T• Ot ю.
□ Т
Уравнение ОГ(ХТ+ Т-цБ) -О, = 0, рав-
□ *
носильное уравнению (X Т + Т -цБ) = 0, получается в виде
П Т Т
X Т + Т-цБ + ш - О, -Т-О, - О, -Т-О,-ш = 0. (3)
Соотношения (1), (3) не являются эквивалентными, и, следовательно, уравнение (1) будет необъективным.
Все тензорные величины можно подразделить на три типа: индифферентные, инвариантные и все остальные. Понятия индифферентности и инвариантности для скаляров совпадают [3]. Инвариантная величина по определению не меняется при преобразовании «звездочка». Понятно, что для объективности уравнения достаточно (но не необходимо), чтобы все входящие в него величины были представителями либо первого, либо второго типов. В нашем случае в уравнение (1), кроме тензоров первого типа, входит тензор - материальная производная Т, который является представителем третьего типа. Материальная производная индифферентного тензора - не индифферентный тензор и не инвариантный. Универсальным средством восстановить индифферентность производной является ее замена на объективную производную, но такая процедура неоднозначна [1]. Существует множество объективных производных и возникает проблема выбора. Самой простой из них, вероятно, будет популярная яуманнская w □
производная Т = Т - W - Т + Т - W, где W -тензор вихря. Однако ее использование может привести к непредсказуемым последствиям, например обнаружена осцилляция напряжений при простом сдвиге [3]. Нас будет интересовать здесь только О-производная Т
Т = Т- П-Т+Т-П, (4)
где П = ОТ - О = - ОТ -О - спин собственно ортогонального тензора поворота О, сопровождающего общую деформацию. Уравнение (1) представляется в виде
хТ + Т-цБ = 0. (5)
Предполагаем, что на отсчетную конфигурацию также наложено жесткое движение с собственно ортогональным тензором От. Символ «звездочка» теперь означает преобразование обеих конфигураций. Попытаемся показать, что уравнение (5) удовлетворяет принципу материальной индифферентности. Имеем:
О* = ОТ - О-О(; О* = ОТ- (О + О- О, О,)-О,;
(ОТ)* = ОТ (ОТ + О, О)Т -ОТ)-От;
(П)* =оТ - П- О, + ш. (6)
Из (2), (4)-(6) следует, что О-производная Т является индифферентным тензором и
(X Т + Т - цБ)* = ОТ - (X Т + Т - цБ) - О, = 0 .
Вероятно, что и требовалось показать. Однако, просматривая доказательство объективности уравнения (5), замечаем, что оно теряет силу, если обе конфигурации совпадают [3, с. 115; 4, с. 83]. В самом деле, тогда имеем От = О, и в выражении для (О* появляет-
т
ся дополнительное слагаемое ш - О, - О - О,. Следовательно, в правой части (6) аддитивно
т
возникает еще тензор -О, РО,, где Р =
т т
= О - О, - ш - О, - О. Таким образом, уравнение
П
О, - (X Т + Т - цБ) - О[ = 0 получается в виде
п
ХТ + Т-цБ + Р-Т - Т-Р = 0. (7)
Уравнения (5), (7) оказываются неэквивалентными. В этом случае уравнение (5) будет необъективным.
Заметим, что объективность уравнений не требуется в геометрически линейной теории, вероятно, ввиду нестрогого обращения там с поворотами, поскольку при их строгом учете обнаруживается, что линейный тензор деформации при жестком движении становится ненулевым [1, с. 26], а это недопустимо. Относительное смещение частиц тела отсутствует и его сплошность не нарушается, Следовательно, по определению деформация не возникает.
Использование О-производной в активном процессе нагружения сыграло решающую роль
при обобщении классической модели нелинейной упругости [1] на нелинейную модель упру-гопластичности [5]. Объективная О-производ-ная, в выражении которой общий поворот заменен на упругий, позволяет однозначно решить отмеченную выше проблему выбора скорости напряжений. С ее помощью двойственное описание процесса в двух пространствах основных мер упругой деформации распространяется на пространства напряжений Т и О - Т - ОТ. Это дает возможность исключить влияние поворотов и определить упругий спин [6].
Построение было выполнено ценой неизбежных потерь. Вводимая промежуточная (разгрузочная) конфигурация является отсчетной для обобщенного упругого закона [4] и, как правило, оказывается неизвестной. Полностью разгрузить тело после появления необратимых деформаций обычно невозможно. Однако всегда могут быть определены основной и взаимный векторные базисы разгрузочной конфигурации. Утверждение о том, что промежуточная конфигурация определена с точностью до поворота [4], на наш взгляд, ошибочное. Разгрузочная конфигурация либо не определена, либо может быть найдена точно там, где все тело деформируется как один элемент.
Законченное определение пластической деформации получается с использованием мультипликативного разложения Крёнера - Ли транспонированного градиента общей дефор-0 0 мации УЫТ (() = Ее - Ер; Ер = (Ее)-1 - VЯТ {г). В условиях изотропии Ер - единичный тензор. Градиент общей деформации определяется известными соотношениями [1]. Строгое определение Ее через его полярное разложение обеспечивает возможность найти Ер. Тензоры Ее, Ер утрачивают смысл градиентов деформации, но оказываются невырожденными. Полярное разложение Ер дает тензоры пластических меры искажений и поворота. Относительно пластической деформации никаких предположений делать не требуется, если строго определить тензор Ее, используя не только кинематические понятия.
Модель нелинейной упругости достаточно разработана, и существуют подходящие изотропные потенциалы напряжений [1], в которые следует добавить анизотропные структуры,
по крайней мере, до второго порядка включительно. Трансформация обобщенного упругого закона (вторая сторона упрочнения по Генки) осуществляется по вводимым определяющим уравнениям для тензора остаточных напряжений и коэффициентов анизотропии [5, 6]. Надо обеспечить и не требующуюся в упругости неотрицательность такого потенциала по его смыслу запасенной энергии, поскольку при обобщении модели нужно постулировать определяющее уравнение для скорости изменения запасенной энергии. Это дает возможность получить дифференциальное уравнение для тензора меры упругих искажений с использованием О-производной.
В [5] дано представление тензора упругого спина, позволяющее найти упругий поворот. Более громоздкое соотношение существует и для трехмерного напряженно-деформированного состояния [6]. Тензор спина разлагается в сумму тензора вихря и некоторого индифферентного тензора. Поскольку яуманнская производная индифферентного тензора индифферентна [1], и О-производная Т будет индифферентна (4). Это обеспечивает объективность предложенных в [5, 6] определяющих уравнений. Рассмотренный выше случай совпадения отсчетной и актуальной конфигураций, когда теряется индифферентность О-производной Т, не реализуется в активном процессе. Тогда имеет место пассивный процесс нагружения, и значит, используется обобщенный закон упругости.
Отметим, что свойства инвариантности относительно отсчетной и актуальной конфигураций тензоров Ее и градиента общей деформации совпадают. Разложение Крёнера - Ли является объективным уравнением.
В Ы В О Д Ы
1. Рассмотрен принцип независимости материальных характеристик от выбора системы отсчета (принцип материальной объективности). Указано условие его совпадения с принципом инвариантности материальных характеристик при наложении жестких движений.
Принцип материальной объективности утверждает, что все материальные характеристики (включая и необъективные) должны вхо-
дить в физические законы специальным образом: при замене системы отсчета изменяется и определяющее уравнение, но полученное уравнение должно быть равносильным исходному.
2. Уравнения геометрически линейной теории часто оказываются непригодными для нелинейной вследствие нарушения принципа материальной индифферентности. Замена материальной производной тензора напряжений Коши на объективную производную позволяет сделать уравнение объективным. На примере классического уравнения Максвелла показаны некоторые тонкости такого преобразования при использовании О-производной (1), (5).
3. Объективная О-производная индифферентного тензора, в частности напряжений Ко-ши, не будет индифферентным тензором только при условии совпадения актуальной и от-счетной конфигураций. Показано, что при активном процессе в обсуждаемой модели упругопластической среды такой случай ис-
ключается. Это позволяет обеспечить обязательную объективность определяющих уравнений при моделировании конечных деформаций и поворотов.
Л И Т Е Р А Т У Р А
1. Лурье, А. И. Нелинейная теория упругости / А. И. Лурье. - М., 1980. - 512 с.
2. Жилин, П. А. Векторы и тензоры второго ранга в трехмерном пространстве / П. А. Жилин. - СПб.: Изд-во СПбГТУ, 2001. - 275 с.
3. Поздеев, А. А. Большие упругопластические деформации: теория, алгоритмы, приложения / А. А. Поздеев, П. В. Трусов, Ю. И. Няшин. - М., 1986. - 232 с.
4. Левитас, В. И. Большие упругопластические деформации материалов при высоком давлении / В. И. Ле-витас. - Киев, 1987. - 229 с.
5. Швед, О. Л. К теории упругопластичности при конечных упругих деформациях и поворотах / О. Л. Швед // Доклады НАН Беларуси. - 2005. - Т. 49, № 3. - С. 45-48.
6. Швед, О. Л. Двойственное описание упругопласти-ческого процесса / О. Л. Швед // Вестник БГУ. - 2007. -№ 2. - С. 88-93.
Поступила 25.10.2006
УДК 629.11.001.24:531.3
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ЗАДЕРЖЕК ПРОЦЕССОВ ПРОИЗВОДСТВА
Канд. техн. наук, доц. ЧЕПЕЛЕВА Т. И., студ. ЧЕПЕЛЕВ А. Н.
Белорусский национальный технический университет, Минский государственный медицинский университет
Одной из ведущих областей машиностроения является производство транспортных машин. Процессы, связанные с производством машин, можно рассматривать как сосредоточенными, так и распределенными. Процесс производства будем считать сосредоточенным, если один и тот же интеллектуальный исполнитель одновременно задействован для нескольких производственных подструктур и управляет ими. Распределенным же будем считать процесс, если несколько различных интеллектуальных исполнителей управляют определенным количеством подструктур.
Для поддержки принятия решений задействованы информационно-аналитические системы, которые предназначены для сопровождения процессов производства машин. При создании систем параллельной обработки данных необходимо знать минимальное время реализации рассматриваемых процессов, а также время взаимодействия параллельных процессов, связанных с производством машин. Реализация сложных параллельных процессов - это особый механизм, который ведет к росту накладных расходов, что значительно влияет на производство машин.
