Full Text

Введение В системах защиты от вибраций и ударов широкое распространение нашли пневматические элементы с резинокордными оболочками [1-3], а также резино-гидравлические виброопоры [4]. При проведении расчетов, как правило, материал оболочек полагается упругодеформируемым [5, 6]. В ряде случаев влиянием резины полностью пренебрегают по сравнению с несущей способностью армирующих нитей корда [7-9]. В то же время предварительные испытания показывают, что релаксация избыточного давления и усилия пневмоэлемента может достигать 25 % и более при выдержке постоянной высоты пневмоэлемента, в котором в качестве рабочей среды используется практически несжимаемая жидкость (например, вода). Причиной этого является вязкоупругость резинокордного композита, и прежде всего резины. Вследствие этого для повышения точности расчетов пневмоэлементов с резинокордными оболочками возникает необходимость в моделировании вязкоупругих свойств резины при конечных деформациях. При этом для симметрично нагруженных оболочек вращения достаточно ограничиться частным случаем, когда главные оси напряжений и деформаций совпадают и не изменяют своей ориентации относительно материальных линий (волокон) оболочки. Основополагающим математическим моделям, учитывающим влияние скорости деформации, посвящена обширная научная литература [10-17]. В обзоре [18] особо выделяется модель Бергстрёма-Бойс [19, 20] как одна из наиболее перспективных современных моделей гипервязкоупругости эластомеров. Среди многочисленных формулировок закона вязкоупругости для резиноподобных материалов модель Бергстрёма-Бойс наиболее популярна [21, 22], так как позволяет достаточно точно описывать экспериментальные данные, включая гистерезис при циклическом нагружении и зависимость диаграмм материала от скорости деформирования. Данная модель нашла широкое применение и признание, прошла экспериментальную проверку и получила дальнейшее развитие [23-28]. Настоящая статья посвящена обсуждению модели гипервязкоупругости Бергстрёма-Бойс, ее критическому анализу и феноменологическому обобщению с целью получения более достоверных результатов моделирования высокоэластичных оболочек вращения при симметричном нагружении. 1. Нелинейная гипервязкоупругая модель Бергстрёма-Бойс Йорген Бергстрем и Мэри Бойс разработали модель гипервязкоупругих материалов на основе соображений микроструктуры [19, 20]. Исходя из анализа экспериментальных данных было выработано положение, по которому механический отклик резиноподобных материалов может быть разложен на две части: реакцию равновесия и отклонение от равновесия, зависящее от времени. Это предполагает, что материал может быть смоделирован в виде двух полимерных сетей, действующих параллельно. Первая сеть фиксирует равновесный отклик материала и может моделироваться любой существующей моделью гиперэластичности. Вторая сеть моделируется идеальной сетью, подобной сети , и включенным последовательно элементом, который действует так, чтобы ослабить нагрузку на идеальную сеть со временем, позволяя описать характеристики, наблюдаемые в экспериментальном исследовании [19, 20]. Наглядно сформулированное положение иллюстрирует реологическая (механическая) модель Кельвина-Пойнтинга1 применительно к случаю одноосного растяжения / сжатия (рис. 1). В модели Кельвина-Пойнтинга упругий элемент параллельно соединен со звеном , называемым элементом Максвелла и состоящим из последовательно соединенных упругого и вязкого элементов. При медленном (квазистатическом) нагружении вязкий элемент не сопротивляется деформации и не воспринимает приложенные напряжения. Поэтому поведение материала определяется деформацией одного упругого элемента . В случае конечной скорости деформации работают все элементы модели, причем сопротивление вязкого элемента (демпфера) зависит от скорости деформации. При очень быстром динамическом (практически мгновенном) нагружении демпфер как бы «запирается», в нем не происходит деформации, и напряжение распределяется между упругими элементами и . В этом случае реакция системы является упругой, однако жесткость повышается по сравнению с медленным упругим деформированием. В начальный момент времени, при деформировании из ненагруженного состояния с произвольной по величине скоростью, напряжение принимает нулевое значение, что является характерной особенностью модели Кельвина-Пойнтинга. Как отмечалось Кельвином [15], упругость такой системы является «совершенной» в том смысле, что вся деформация полностью исчезает при разгрузке и является поэтому «упругой» (обратимой). Рис. 1. Реологическая модель Кельвина-Пойнтинга Fig. 1. Rheological model of Kelvin-Poynting При построении определяющих соотношений в [19] материал полагается однородным, изотропным и сжимаемым. Градиент полных деформаций считается одинаковым для обеих ветвей: (здесь и - радиусы-векторы материальной точки в отсчетной и актуальной конфигурациях соответственно). Градиент деформации представляется через градиент упругих деформаций и градиент вязких (неупругих) деформаций по формуле мультипликативного разложения (рис. 2): . (1) Тензор переводит локальную (бесконечно малую2) отсчетную конфигурацию материала в локальную промежуточную конфигурацию , получаемую из локальной актуальной конфигурации в результате полной виртуальной упругой разгрузки сети (см. рис. 1) в состояние без напряжений (рис. 2). По этой причине в [19] конфигурация называется разгруженной конфигурацией (relaxed configuration). По теореме Коши о полярном разложении градиент можно представить в виде , (2) где - собственный ортогональный тензор поворота, сопутствующий деформации; и - правый и левый тензоры растяжения, которые являются симметричными положительно определенными тензорами. Аналогичные разложения имеют место для градиентов и (см. рис. 2): , (3) . (4) Рис. 2. Мультипликативное разложение деформации: , , - бесконечно малая окрестность материальной точки среды в отсчетной, актуальной и промежуточной конфигурациях Fig. 2. Multiplicative strain decomposition: , , - is the infinitely small neighborhood of the material point of the medium in the current, current and intermediate configurations В свою очередь, левые тензоры растяжения , можно представить через их спектральное разложение: , , (5) где , - собственные числа (кратности удлинения); , - собственные векторы тензоров , соответственно. Тензор напряжений складывается из тензора напряжений сети и тензора напряжений сети (см. рис. 1): . (6) Тензоры , выражаются в [19] через главные напряжения: , . (7) Главные напряжения схематизируются по модели Арруда-Бойс [30], согласно которой квазистатический отклик материала относится к отклику кубической ячейки, содержащей восемь цепей Ланжевена, связывающих центр ячейки с его вершинами [30]. При этом получаются выражения (8) Здесь , - начальные модули упругости; , - предельные растяжения сети; - объемный модуль упругости, , (9) - кратность усредненного удлинения макромолекулярной цепи эластомера (эффективное деформационное растяжение цепи), , , (10) , (11) - первый и третий главные инварианты правых мер полных и упругих деформаций Коши-Грина; , , - функция, обратная функции Ланжевена . При использовании, например, модели гиперупругости Муни-Ривлина [18], имеющей упругий потенциал (12) главные напряжения определяются выражением [5] (13) с подстановками вида , . Здесь , , - главные инварианты тензора (для несжимаемого материала , [5]). Параметры , , - материальные постоянные, которые могут быть определены путем подгонки кривой материала в условиях одноосного квазистатического растяжения / сжатия. Модель Бергстрёма-Бойс допускает использование и других упругих потенциалов, прошедших экспериментальную проверку [5, 31, 32]. В более поздней работе Бергстрема и Бойс [20] определяющие соотношения (7), (8) представляются в тензорной (инвариантной) форме записи с несущественными изменениями в обозначениях и в представлении членов, учитывающих объемную деформацию. Вариации определяющих соотношений (7), (8) содержатся в [21, 22, 27, 28]. При кинематическом описании полимерной сети (см. рис. 1) полный градиент скорости разлагается на упругую и вязкую (неупругую) составляющие: , (14) где , , . (15) Для замыкания кинематических соотношений (14), (15) и обеспечения единственности мультипликативного разложения (1) Бергстрём и Бойс [19, 20] со ссылкой на обоснование [33] принимают для кососимметричного тензора соотношение . (16) Тензор скоростей вязких деформаций выражается для сети (см. рис. 1) через девиатор тензора напряжений по закону течения , (17) где , - девиатор и норма тензора соответственно; - единичный тензор. Фактически при постулировании (17) в [19, 20] неявно принимается, что материал является «вязко (неупруго) несжимаемым»3: . Интенсивность скорости вязких деформаций аппроксимируется в [19] степенным законом (в [20] без множителя благодаря перенормировке): , (18) где , , , - материальные постоянные (в приложениях используются только три материальные постоянные , , ), , , (19) - собственные числа правой меры вязких деформаций Коши-Грина . В формуле (18) множитель получен из соображений микроструктуры эластомера рассмотрением динамики рептационного движения полностью или частично неактивных полимерных цепей [19] при эвристическом предположении, что . Модель Бергстрёма-Бойс (1)-(19) содержит восемь материальных параметров , , , , , , , , заданием которых полностью определяется поведение вязкоупругого эластомера. По рекомендациям, приведенным в приложении к статье [19], первые пять параметров , , , , могут быть экспериментально оценены путем анализа квазистатических диаграмм одноосного растяжения / сжатия и диаграмм поведения материала во время нагрузки / разгрузки при циклическом одноосном растяжении / сжатии с заданной скоростью деформации. Оценка остальных трех параметров , , (особенно константы ), напротив, может быть выполнена только методом проб и ошибок, основанным на сравнении экспериментальных данных на одноосное растяжение / сжатие в заданном диапазоне скоростей деформации с соответствующими численными решениями уравнений математической модели [18, 19]. Таким способом успешные результаты получены при малых скоростях деформирования [25] и больших скоростях деформации при ударных испытаниях на разрезном стержне Гопкинсона [26]. При умеренно высоких скоростях деформирования в [21, 22] параметры модели подбирались отдельно для каждого динамического режима испытаний. Поэтому процедура определения единого набора значений параметров модели, пригодных для описания экспериментальных данных при разных скоростях нагружения, требует дополнительного исследования [21, 22]. 2. Критический анализ модели Бергстрёма-Бойс Центральной проблемой моделей вязкоупругости и пластичности, использующих мультипликативное разложение (1), является его неединственность [33-39]. Здесь имеется в виду тот факт, что если тензоры , удовлетворяют разложению (1), то аналогичное будет верно для тензоров , с произвольным ортогональным тензором поворота : . Последнее означает, что тензоры поворота , в разложениях (3), (4) оказываются неопределенными. Иными словами, само по себе разложение (1) не позволяет разделить упомянутые два вращения. В [37] рассматриваются следующие два варианта решения проблемы4: {1} и , {2} и , где - тензор поворота в полярном разложении градиента деформации (напомним, что в модели Бергстрёма-Бойс , ибо ). Конечное предпочтение отдается первому условию однозначности из соображений математического удобства. В последующей статье [33] помимо вышеуказанных двух условий однозначности вводится условие, ограничивающее значение спина вязких деформаций , входящего в (15). Из соображений удобства и простоты предлагается следующее условие однозначности: {3} , которое, очевидно, отлично от условия (16). При использовании ранее построенной модели [37] для решения тестовой задачи простого сдвига получаются идентичные результаты при обоих условиях «{2}» и «{3}». На основании этого в [33] делается вывод, что выбор конкретного ограничения не существенен при решении проблемы неоднозначности мультипликативного разложения. Допустимость ограничения «{3}» на спин вязких деформаций детально исследуется в [38]. Формулируется следующее положение: если у задачи есть решение, то любое зависящее от времени и пространства вращение промежуточной (разгрузочной) конфигурации (см. рис. 2) также дает (идентичное) решение и можно выбрать это вращение таким образом, чтобы преобразованное решение неупруго не вращалось: . При этом отмечается: «Результат такого характера кажется частью фольклора теории пластичности для изотропных материалов, но мы не знаем фактического доказательства» [38]. Затем такое доказательство приводится (на основании принятых допущений), после чего делается вывод, что при обсуждении проблемы неоднозначности мультипликативного разложения, связанной с изотропными вязкопластическими и аморфными твердыми телами, можно без потери общности ограничить внимание решениями с . Для полноты общей картины следует отметить, что в историческом обзоре теории пластичности [36] как допустимые варианты отмечаются условия однозначности , , , а в монографии [35] - еще два условия , . Чтобы понять, насколько существенным (или не существенным) является выбор конкретного условия однозначности мультипликативного разложения градиента деформации, применим модель гипервязкоупругости Бергстрёма-Бойс к решению задачи простого сдвига, которому соответствует закон движения , , , где - абсолютный сдвиг; - угол сдвига. Возьмем три условия однозначности , , . (20) Значения материальных параметров заимствуем из [19]: = 0,2 МПа, = 1,6, = 8, = 100 МПа, = 4, = 7 МПа-m/с, = -1. Рассмотрим один цикл нагружения и разгрузки, при котором деформирование происходит с постоянной равной скоростью = 2,5 с-1 до максимального значения = 2,5, соответствующего максимальному углу сдвига = 68,2°. Поскольку квазистатические напряжения в полимерной сети не зависят от выбора условия однозначности (20), наибольший интерес представляют динамические напряжения в полимерной сети (рис. 3). Как видим, влияние того, какой вариант будет выбран, является существенным. Однако отличия в значениях полных напряжений (рис. 4) менее значительны, что позволяет признать модель гипервязкоупругости Бергстрёма-Бойс вполне удовлетворительной при описании напряженного состояния по сравнению с описанием деформированного состояния промежуточной конфигурации не только в отношении поворотов, но и изменения расстояний между материальными точками в каждый момент времени (рис. 5). Рис. 3. Динамические напряжения (в полимерной сети ) при простом сдвиге по модели гипервязкоупругости Бергстрема-Бойс при разных условиях однозначности: 1 - ; 2 - ; 3 - Fig. 3. Dynamic stresses (in the polymer network ) under simple shear according to the Bergström-Boyce model of hyper-viscous elasticity under different conditions of uniqueness: 1 - ; 2 - ; 3 - Рис. 4. Полные напряжения при простом сдвиге по модели гипервязкоупругости Бергстрёма-Бойс при разных условиях однозначности: 1 - ; 2 - ; 3 - Fig. 4. Total stresses under simple shear according to the Bergström-Boyce hyper-viscoelasticity model under different conditions of uniqueness: 1 - ; 2 - ; 3 - Рис. 5. Конфигурации М-образца (отсчетная , актуальная , промежуточная ) при простом сдвиге по модели гипервязкоупругости Бергстрёма-Бойс на этапе разгрузки при угле сдвига = 5° в момент времени t = 1,965 c при разных условиях однозначности: 1 - ; 2 - ; 3 - Fig. 5. Configurations of the M-sample (reference , current , intermediate ) under simple shear according to the Bergström-Boyce hyper-viscoelasticity model at the unloading stage at a shear angle = 5° at the time t = 1,965 s under different conditions of uniqueness: 1 - ; 2 - ; 3 - В заключение необходимо указать критерий, позволяющий отделить допустимые условия однозначности от физически неприемлемых условий. С этой целью обратимся к принципу объективности [35, 40-42]. Так, например, для упругих материалов, если в одной системе отсчета (без звездочки) для тензора напряжений имеет место определяющее соотношение в виде функции , то в любой другой системе отсчета (со звездочкой) имеет место определяющее соотношение в виде аналогичной функции . Поскольку условие однозначности в моделях, использующих мультипликативное разложение градиента деформации, является неотъемлемой составной частью определяющего соотношения (а не уравнения движения) для тензора напряжений, представляемого в виде одного или системы нескольких уравнений, то принцип объективности должен распространяться и на условие однозначности (постулируемое положение). Если без ограничения общности в формуле Нолла для преобразования времени принять , то тогда при замене системы отсчета тензор напряжений преобразуется по формуле перехода [40-42] , (21) где - ортогональный тензор, характеризующий относительный поворот систем отсчета. Если - векторы базиса сопутствующей лагранжевой системы координат (, , - координаты Лагранжа материальной точки) в данной системе отсчета, то по формуле Нолла в другой системе отсчета . (22) В начальный момент времени (в отсчетной конфигурации) по формуле (22) . (23) Если ввести обычные упрощенные обозначения без указания аргументов: , , , , то тогда выражения (21)-(23) можно переписать в компактном виде: , , . (24) Для градиента деформации справедливо диадное разложение [35] , (25) где - векторы взаимного (дуального) базиса сопутствующей лагранжевой системы координат: , - символы Кронекера. На основании (24) формула преобразования градиента деформации (25) принимает вид . (26) Если в начальный момент времени у обеих систем отсчета оси базисных декартовых систем координат совпадают5, тогда и формула (26) принимает традиционный вид . Произвольность и постоянство тензора в (26) существенно упрощает технику применения принципа объективности для построения определяющих соотношений с использованием теорем об изотропных функциях при переходе от одной инерциальной системы отсчета к другой инерциальной системе отсчета [44]. Не так однозначно обстоит дело с формулами преобразования градиентов упругих и вязких деформаций, входящих в формулу мультипликативного разложения (1). Так, например, в [36] приводятся формулы , , (27) а в [38, 39] утверждается иное: , . (28) В зависимости от этого одно и то же условие однозначности будет либо допустимым, либо физически неприемлемым. Чтобы разобраться в данном вопросе, обратимся к диадным разложениям градиентов деформации , (см. рис. 2), аналогичным (25): , . Здесь , - векторы основного и взаимного базисов сопутствующей лагранжевой системы координат в промежуточной конфигурации (), которые при замене системы отсчета преобразуются согласно (22): , . Следовательно (при ), если промежуточная конфигурация отнесена к начальному моменту времени , тогда будут справедливы формулы преобразования (28), если же промежуточная конфигурация отнесена к текущему моменту времени , тогда следует использовать формулы преобразования (27). Если говорить о реально осуществляемой разгрузке из локальной актуальной конфигурации (см. рис. 2), тогда , и поэтому оба набора формул (27), (28) будут неверными. Если исходить из развернутой записи разложения (1) в виде , передающей суть математической модели, тогда становится понятным, что тензоры , , характеризуют деформированное состояние в одной и той же материальной точке среды (с радиусами-векторами , , в отсчетной, промежуточной и актуальной конфигурациях соответственно) в один и тот же текущий момент времени t. Об этом свидетельствует, в частности, рис. 5. Поэтому более уместно называть конфигурацию (см. рис. 2) промежуточной (или вспомогательной) конфигурацией, как в [36], а не (изохронно) разгруженной конфигурацией, чтобы не возникало излишних ассоциаций с термином «разгрузка». Таким образом, в качестве следует брать текущий момент времени , поэтому формулы преобразования градиентов упругих и вязких деформаций будут иметь вид , . (29) При формулы (29) переходят в (27). В частности, исходя из формулы (29), формул полярного разложения Коши и формул (14), (15) можно показать, что , , , , , , , , , , , , , (30) , , где , - симметричная и кососимметричная части тензора . В таблице приведены соотношения, принятия которых достаточно для обеспечения однозначности мультипликативного разложения градиента деформации. Только часть из них удовлетворяет требованиям независимости от выбора системы отсчета. При желании список такого рода соотношений можно многократно расширить. Вопрос, какое именно условие однозначности следует использовать, должен решаться так же, как это делается, например, для упругих потенциалов: надлежащей постановкой и проведением экспериментальных исследований либо теоретическим исследованием микроструктуры материала6. К сказанному следует добавить, что согласно формуле (30) тензор является необъективным, вследствие чего определяющее соотношение (17) не удовлетворяет принципу объективности поведения материалов. Последнее можно исправить, если по примеру [25] заменить в (17) тензор , тензором . Исчерпывающее решение затронутого вопроса требует отдельного исследования. Поэтому дальнейшее изложение (в рамках феноменологического подхода) ограничено важным для приложений случаем, когда можно обойтись без конкретизации условия однозначности. Соответствие условий однозначности принципу объективности материалов Compliance of the conditions of uniqueness to the objectivity principle of materials Система отсчета Допустимость (+) Неприемлемость (-) Старая (без звездочки) Новая (со звездочкой) + - - + + - - - + 3. Феноменологическая модель вязкоупругости эластомеров Перейдем к рассмотрению прикладной задачи, касающейся высокоэластичных оболочек, и воспользуемся тем упрощающим обстоятельством, что при симметричном нагружении оболочек вращения орты меридиана и параллели во все моменты времени совпадают по направлению с главными осями напряжений и деформаций (полных, упругих и вязких) и не изменяют своей ориентации относительно материальных линий (волокон) оболочки. Благодаря соосности градиентов полных, упругих и вязких деформаций соответствующее мультипликативное разложение является однозначным. При этом ограничение на неизменность ориентации главных осей напряжений относительно материальных линий (волокон) является существенным (пока надлежащим образом проведенные испытания эластомеров не покажут обратное) в свете исследований поликристаллических материалов (конструкционных сталей и сплавов) при непропорциональных циклических нагружениях [46]. Перейдем в сопутствующую систему отсчета, относительно которой бесконечно малая окрестность материала не вращается7, и рассмотрим представительный элемент однородного изотропного несжимаемого8 эластомера (например, резины) в форме прямоугольного параллелепипеда (рис. 6). После нагружения нормальными (главными) напряжениями , , он деформируется без поворота и с сохранением прямоугольной формы. Вследствие этого градиент деформации будет иметь диагональный вид: , (31) где , , - кратности удлинения соответствующих ребер параллелепипеда. а б Рис. 6. Однородная деформация представительного элемента эластомера: а - до нагружения; б - после нагружения Fig. 6. Homogeneous deformation of the representative element of the elastomer: a - before loading; b - after loading По выражению (31) определяется тензор скоростей полной деформации: , (32) где - симметричная часть тензора . Схожий вид имеет также тензор напряжений (рис. 6, б): . (33) Далее будем опираться на реологическую модель Кельвина-Пойнтинга (см. рис. 1), которая при феноменологическом подходе не только наглядно отображает поведение материала, но и играет активную роль при выводе определяющих соотношений, задавая их структуру [13]. Следуя [21, 22] и методическим рекомендациям [13], на основании реологической модели Кельвина-Пойнтинга (см. рис. 1) можно констатировать следующее. При параллельном соединении звеньев и они деформируются одинаково: , (34) а напряжения складываются из откликов упругих элементов звеньев и : . (35) Здесь - градиент полной деформации материала (упругого элемента ); - градиент упругой деформации материала (упругого элемента ). Поскольку в звене упругий элемент последовательно соединен с вязким элементом, должно выполняться равенство , (36) где , - градиент вязкой деформации материала (вязкого элемента ) и скорость его изменения соответственно. Градиенты деформации , , связаны между собой формулой мультипликативного разложения , (37) предложенной и впервые использованной Е.Х. Ли для описания больших неупругих деформаций [34, 35]. Зависимости (34)-(36), со всей очевидностью, порождены исключительно структурой реологической модели Кельвина-Пойнтинга (см. рис. 1); они предопределяют дальнейший вид получаемых определяющих соотношений гипервязкоупругого материала. Для детализации математической модели требуется задать свойства структурных элементов реологической модели. Для описания свойств упругих элементов и воспользуемся определяющими соотношениями гиперупругого несжимаемого материала: , (38) . (39) Для вязкого элемента возьмем определяющее соотношение . (40) Здесь , , - неопределенные множители Лагранжа, наличие которых обусловлено несжимаемостью материала [40-42]; , - упругие потенциалы, , , (41) - правые меры полной, упругой и вязкой деформации Коши-Грина соответственно; (42) - тензор скоростей вязкой деформации, вводимый наряду с выражениями , (43) для тензора скоростей упругой деформации и тензора скоростей полной деформации. Поскольку третий главный инвариант , коэффициент вязкости как изотропная скалярная функция тензора будет определяться через первые два главных инварианта ,: . (44) Замечание. Определяющее соотношение (40) при постоянном коэффициенте вязкости соответствует несжимаемой жидкости Рейнера-Ривлина [47-49], у которой или в наиболее общем виде . (45) Поперечная вязкость полагается равной нулю, а кажущаяся вязкость - степенной функцией от второго главного инварианта или, что то же самое, нормы тензора скоростей деформации , являющегося девиатором: . При этом , где - неопределенный множитель Лагранжа (гидростатическое давление, взятое со знаком минус). С учетом этого после подстановки в (45) выражения получается определяющее соотношение , , (46) в котором коэффициент вязкости и показатель степени являются положительными материальными параметрами и могут зависеть от температуры. Когда , соотношение (46) переходит в реологический закон Ньютона для линейно вязких (ньютоновских) несжимаемых жидкостей классической гидродинамики. Если же , то соотношение (46) становится подобным закону сухого трения Кулона-Амонтона; на нем основывается классическая теория пластичности и ее современные обобщения [50-52]. Следует также отметить, что зависимость (46) в записи (47) является наиболее распространенной аналитической зависимостью в теории течения - одном из вариантов прикладной теории ползучести [53, 54]. В другом, более общем варианте прикладной теории ползучести, называемом теорией упрочнения [53, 54], скалярный параметр полагается зависящим от линейного тензора деформации ползучести . Как правило, используется степенная зависимость вида , (48) где , - материальные константы. Зависимость (48) относится главным образом к таким конструктивным материалам, как металлы и их сплавы при высоких температурах. Применительно к эластомерам исходя из основополагающих зависимостей (17), (18) модели Бергстрёма-Бойс и равенства , выполняющегося в рассматриваемом случае, для коэффициента вязкости получается степенная зависимость , , . (49) Зависимость (49), имеющая вид , является частным случаем зависимости (44). Указанные зависимости соотносятся между собой так же, как, например, неогуковский упругий потенциал Трелоара [31, 41] соотносится с упругим потенциалом Муни-Ривлина (12). Последнее может быть важным для правильной настройки модели при описании зависимости от скорости вязких деформаций. Вернемся к описанию процесса деформирования представительного элемента вязкоупругого материала (см. рис. 6). В условиях рассматриваемого нагружения упругая и вязкая составляющие градиента полной деформации имеют вид, аналогичный (31): , (50) . (51) Подставляя (31), (50), (51) в формулу (37), будем иметь . Следовательно, , , . (52) Дифференцируя выражения (52) по времени и подставляя результат в (32), придем к формуле аддитивного разложения тензора скоростей деформации , (53) где , , (54) формулы диадного разложения тензоров скоростей полной, упругой и вязкой деформации (42), (43) соответственно. С помощью (13) определяющие соотношения (38)-(40) можно переписать в более удобном для расчетов виде: , , (55) , , (56) , . (57) Здесь используются обозначения , . (58) Подставляя (55), (56) в (35), с учетом (33) находим (), (59) где - суммарный множитель Лагранжа. Аналогичным образом, подставляя (56), (57) в (36), получаем (), (60) где обозначено . Уравнения (59), (60) совместно с выражениями (52)-(54), (58) и (49) приводят к следующей системе уравнений: () (61) Система уравнений (61) представляет собой искомую систему определяющих соотношений гипервязкоупругого материала, соответствующую реологической модели Кельвина-Пойнтинга (см. рис. 1). К ним надо добавить условия несжимаемости: , , . (62) С помощью (52) нетрудно убедиться, что из любых двух условий несжимаемости (62) вытекает третье условие, поэтому независимыми в (62) являются только два условия несжимаемости. Из (62) также вытекает, что (63) Построение феноменологической модели гипервязкоупругости заканчивается заданием явной зависимости (44) для коэффициента вязкости . По-видимому, во многих практически важных случаях достаточно точной является зависимость Бергстрёма-Бойс (49). Если, как в исследованиях полиуретана [21, 22], при использовании (49) возникнут затруднения в определении единого набора значений параметров модели, пригодных для описания экспериментальных данных при разных скоростях нагружения, то тогда может потребоваться поиск более общего выражения для коэффициента вязкости , зависящего не только от первого, но и второго инварианта правой меры вязких деформаций Коши-Грина9. Определяющие соотношения (61) пригодны для экспериментального определения материальных параметров модели по стандартным испытаниям на одноосное растяжение-сжатие, двуосное растяжение-сжатие и чистый сдвиг [5, 31]. Для этих целей, а также для описания напряженно-деформированного состояния оболочек вращения при симметричном нагружении определяющие соотношения (61) с учетом условий несжимаемости (62), (63) должны быть преобразованы к удобному для численных расчетов виду, что выходит за рамки текущей статьи и является предметом отдельного исследования. Заключение Для удовлетворения практической потребности в повышении точности расчетов высокоэластичных оболочек должна привлекаться теория конечных вязкоупругих деформаций эластомеров. Среди современных моделей гипервязкоупругости резиноподобных материалов одной из наиболее перспективных и популярных является модель Бергстрёма-Бойс, основанная на сведениях о микроструктуре эластомеров и использующая мультипликативное разложение градиента полной деформации. Прямое решение тестовой задачи простого сдвига по модели Бергстрёма-Бойс показывает, что получаемые расчетные результаты существенно зависят от выбора условия однозначности, обеспечивающего единственность мультипликативного разложения градиента деформации. Из множества предлагаемых в научной литературе условий однозначности только определенная часть удовлетворяет требованию независимости от выбора системы отсчета. При желании список подобных условий однозначности можно многократно расширить. Поэтому вопрос, какое именно условие однозначности следует использовать, должен решаться так же, как это делается, например, для упругих потенциалов: надлежащей постановкой и проведением экспериментальных исследований либо теоретическим исследованием микроструктуры материала. Предложенная феноменологическая модель гипервязкоупругости, основанная на одномерной реологической модели Кельвина-Пойнтинга, ограничивается случаем, когда главные оси напряжений и деформаций (полных, упругих и вязких) совпадают и не изменяют своей ориентации относительно материальных линий (волокон). Благодаря этому отпадает необходимость в формулировке условия однозначности, обеспечивающего единственность соответствующего мультипликативного разложения. Учет зависимости коэффициента вязкости от второго инварианта правой меры вязких деформаций Коши-Грина в степенном законе кажущейся вязкости модели Рейнера-Ривлина обобщает соответствующую зависимость модели Бергстрёма-Бойс и позволяет, при необходимости, расширить диапазон скоростей деформации при согласовании расчетных и экспериментальных данных. Разработанная математическая модель гипервязкоупругости резиноподобных материалов предназначена для расчета напряженно-деформированного состояния высокоэластичных оболочек вращения при симметричном нагружении.

About the authors

V S Korneyev

Omsk State Technical University

S A Korneyev

Omsk State Technical University

References

Statistics

Views

Abstract - 307

PDF (Russian) - 228