Математическая модель поведения пластификатора в неоднородно нагруженном эластомере
Рассмотрим меру упругой деформации Коши — Грина G.
Она представляет собой симметричный положительный тензор, т. е.
1) его правые собственные векторы gi являются одновременно и левыми:
2) собственные значения представляют собой неотрицательные числа (обозначим их квадратами некоторых параметров λi); 3) собственные векторы gi всегда могут быть представлены ортонормированной тройкой векторов. Обозначим символом О тензор поворота, переводящий базисные векторы прямоугольной Декартовой системы координат ei в ортонормированную тройку собственных векторов gi тензора G:
Нa основании сказанного, меру упругой деформации Коши — Грина G можно записать в виде выражения от ее собственных значений λi2, тензора поворота О и базисных векторов ei:
а квадраты упругих удлинений вдоль главных осей λi2 представить с помощью формул
В данной работе рассматриваются материалы, плотности свободных энергий компонент которых fi являются функциями температуры среды 0, характеристик ее обратимых деформаций λi и концентраций жидких компонент смеси φ1,...,φN. Полагаем, что они имеют следующий вид:
Под концентрациями жидких компонентов смеси подразумеваются выражения
В зависимости от физического смысла констант βj они могут иметь значения объемной доли жидкого компонента, массовой доли, мольной доли (отношение числа молекул данного вида к общему числу молекул в малом объеме) и т. п.
Изменение параметров состояния упругой среды
Рассмотрим, как меняются по времени характеристики обратимых деформаций λi. Продифференцируем равенство (35) по времени t. С помощью кососимметричпого тензора W
Правила свертки тензоров и кососимметричность тензора W позволяют упростить полученное выражение:
Таким образом, задача расшифровки закона изменения λi во времени свелась к анализу поведения производной от тензора G.
Еще раз отметим, что пашей целью является формулировка всех уравнений в координатах t, r*. Это касается и особенностей записи тензора G. Перейти к ней позволяет связь
Она справедлива, поскольку полагается возможным выражение векторов ra и r* друг через друга и оно не зависит от времени t:
В результате мера упругой деформации Коши — Грина G представляется выражением
От него нетрудно взять производную по времени t. Поскольку тензор QA не зависит от времени, искомая производная имеет вид
Изменим правую часть равенства (40), включив в слагаемые свертку с единичным тензором:
Связь (38) позволяет записать в новом, более коротком виде полученное выражение:
Для дальнейшего анализа потребуется более конкретная расшифровка смысла тензора Qа, Рассмотрим ее. В соответствии с теоремой о полярном разложении тензоров, справедливо равенство
в котором положительный симметричный тензор H и ортогональный тензор Oа определяются формулами
На основании представлений (32), (34) легко установить непосредственный вид тензора Н:
Поэтому искомый тензор Qа может быть представлен с помощью записи
в которой символом Оf обозначен тензор поворота
Легко проверить, что это тот самый ортогональный тензор, который получается при представлении собственных векторов hi меры упругой деформации Фингера F через базисные векторы ei прямоугольной Декартовой системы координат (подобно тому, как это было сделано выше при анализе вида меры деформации Коши — Грина). Отметим также, что собственными значениями тензора F выступают величины λi2, которые являются одновременно и собственными значениями тензора G:
Оно позволит нам переписать производную по времени от параметра λi (37) в новом виде:
Используя формулировку (26) и связи (33), (44), получаем окончательный вид искомого выражения:
Определяющие уравнения процессов
Вернемся к анализу термодинамического неравенства (31). Массовые плотности свободных энергий fi зависят только от параметров 0, φ1,...,φN, λ1, λ2, λ3 (согласно сделанному предположению (36)). В свою очередь величины φ1,...,φN выражаются через характеристики содержания жидких компонентов смеси Q1*,...,QN*. Они изменяются с течением времени (в противоположность неизменной характеристике Q0*). Поэтому в дальнейшем изложении всюду будем полагать в качестве изменяемых параметров состояния среды величины 0, Q1*,...,QN*, λ1, λ2, λ3, ξ1. Раскроем значения производных от плотностей fi:
Суммирование по индексу к начинается не с нуля, а с единицы, в соответствии с условием (18). Конкретный вид производной от параметра λi по времени t мы уже определили (45). Используя обозначение
перепишем термодинамическое неравенство (31) в более удобном для анализа виде:
Учитывая, что скорости vi* определены выражениями (7), справедливо равенство (19) и законы сохранения массы (17) могут быть сформулированы в виде
перепишем ограничение:
Таким образом, требование непротиворечивости математической модели первому и второму законам термодинамики свелось к удовлетворению неравенства диссипации, сформулированного в виде. Существует множество способов построить определяющие уравнения поведения материала, которые удовлетворяют этому неравенству. Ниже приводится одна из возможных физически разумных систем уравнений. Рассмотрим ее.
Пусть поведение континуума определяется связями
Причем для любого набора векторов аi (r = 1,..., N) должно выполняться неравенство
Первое равенство в системе уравнений (50) фиксирует известное в термодинамике утверждение. Массовая плотность энтропии среды si равна производной от массовой плотности свободной энергии fi по температуре 0, взятой с противоположным знаком.
Следующие выражения (47), (51) известны в механике. Величина главных значений равновесных компонент напряжений σei однозначно определяется массовыми равновесными плотностями свободной энергии среды fi. При этом их пространственная ориентация определяется собственными векторами hi меры деформации Фингера F.
Равенство (56), переписанное с учетом связи между операторами градиента места в актуальной r и рассматриваемой r* конфигурациях
представляет собой закон теплопроводности Фурье. Его физический смысл сводится к утверждению того, что передача тепла в континууме происходит всегда из более горячих областей материала в более холодные.
Рассмотрим теперь особенности моделирования массообменных процессов. Подставим для этого значения напряжений Ti, действующих на r-й континуум, и силы взаимодействия fi из формул (52)-(55) в уравнение движения (24):
Определим аналоги химических потенциалов жидких компонентов смеси с помощью выражений
Они отличаются от реальных химических потенциалов /л,- только сомножителем Qio:
где mi — масса одного моля r-го континуума. В результате закон движения (59), с учетом связи между операторами градиента места в актуальной r и рассматриваемой r* конфигурациях, переписывается в следующей формулировке:
Таким образом, движение i-го континуума происходит вследствие неоднородности химического потенциала μi (или его аналога μio) и поля температур 0. При этом движению противодействуют силы вязкого взаимодействия компонентов смеси ηQik Qk (vk — v0).
Последнее, в чем осталось убедиться, — это справедливость неравенства (49). Сразу следует отметить, что симметричность полного тензора напряжений T и кососимметричность тензора WD обращают в ноль результат их двойной свертки:
Используя выражения (50)-(56), (62), преобразуем ограничение (49):
Второе слагаемое в полученном выражении не может принимать отрицательных значений. Доказывается это с помощью следующего преобразования
Первое слагаемое неотрицательно в соответствии с условием (58).
Для моделирования проходящих в среде процессов необходимо сформулировать еще одно уравнение — теплопроводности. Для его получения необходимо подставить в закон сохранения энергии (27) с помощью связи (30) выражение массовой плотности еi через характеристики fi, 0, si:
В равенстве (64) необходимо раскрыть значения производных но времени от величин Это осуществляется с помощью приведенных ниже формул. Первая получается па основе зависимостей (45), (47), (51) и имеет вид
В итоге рассматриваемое выражение (64) при подстановке выписанных формул и зависимостей (50), (62) принимает вид
Эта связь и представляет собой закон теплопроводности. Изменение энтропии системы происходит за счет переноса ее потоками компонентов смеси, теплообмена среды с соседними областями, в результате выделения тепла в условиях вязкого взаимодействия компонентов смеси при движении друг относительно друга.
Краткая информация о модели
Для записи математических уравнений в отсчетной конфигурации используем основные характеристики состояния среды, перечисленные в табл. 4. Часть из них представляют собой известные физические параметры состояния среды, другая часть является удобным обозначением нужных нам математических выражений.
Все процессы в материале в рамках рассматриваемой модели описываются уравнениями, сформулированными в табл. 5.
Свойства материала задаются скалярными функциями fi, ηq, ηQik (табл. 6). Их необходимо определять экспериментально. При этом величины ηq, ηQik могут быть функциями параметров состояния материала 0, φ1,...,φN, λ1, λ2, λ3 и их производных по времени. Важно только, чтобы удовлетворялись неравенства (57), (58).
Тензор напряжений среды Т, плотности энтропии компонентов смеси si и аналоги химических потенциалов жидких составляющих материала μi0 однозначно выражаются через плотности свободных энергий fi.
- Основные термодинамические законы и их следствия в теории Грина и Надхи
- Теории массообменных процессов
- Структурные и эффективные характеристики поврежденных систем
- Структурные характеристики ансамблей со случайным расположением включений
- Эффективные и структурные характеристики ансамблей с регулярным расположением включений