Компьютерное моделирование деформирования эластомеров

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Компьютерное моделирование деформирования эластомеров

Эластомеры могут претерпевать большие деформации (несколько сотен процентов) без разрушения и повреждения структуры материала. В силу этих уникальных свойств они широко используются в технике. Потребность в математическом моделировании процессов деформирования тел и конструкций из таких материалов стимулирует развитие теории больших деформаций гиперупругих тел, создание алгоритмов численных решений уравнений гиперупругости и их программную реализацию. Целью настоящей работы является сравнительный анализ математического моделирования умеренных деформаций эластомеров (до 50-ти процентов) с использованием моделей гиперупругих материалов Сен-Венана, Муни - Ривлина и Генки в среде пакета MSC. Marc [1]. Первые две модели имеются в виде стандартных моделей материалов пакета MSC. Marc, а модель материала Генки адаптирована и внедрена в этот пакет авторами настоящей работы [2, 3]. Данные компьютерного моделирования сопоставляются с результатами натурных экспериментов.

Определяющие соотношения гиперупругости. В рамках уравнений линейной теории упругости для изотропных упругих материалов используется закон Гука [4]

дуотан эластомер деформирование

, (1)

где и - соответственно, тензоры напряжений и деформаций Коши, - единичный тензор второго порядка, «» обозначает след тензора второго порядка, и - константы Ламэ, которые выражаются через модуль Юнга и коэффициент Пуассона следующим образом

(2)

Закон Гука (1) можно переписать следующим образом

(3)

где операция «:» обозначает двойное внутреннее произведение тензоров [4], а тензор четвертого порядка имеет следующий вид:

(4)

где , , - базовые изомеры изотропного тензора четвертого порядка [4]. Закон Гука допускает также запись в потенциальном виде

(5)

Уравнения линейной теории упругости (а, следовательно, и закон Гука) описывают только малые (не более трех процентов) деформации твердых тел. Для математического моделирования деформирования тел из упругих материалов надо использовать уравнения нелинейной теории упругости. При записи определяющих соотношений эта теория требует использования объективных по Лагранжу или Эйлеру тензоров напряжений и деформаций, сопряженных по мощности внутренних сил [4-7]. Так как тензор деформаций Коши не объективен, то закон Гука, представленный в альтернативных формах записи (1), (3) или (5), нельзя прямо использовать в уравнениях нелинейной теории упругости.

В нелинейной теории упругости различают три типа определяющих соотношений: гиперупругости (упругости по Грину), упругости (упругости по Коши) и гипоупругости [4]. В рамках уравнений линейной теории упругости для изотропных материалов все три типа определяющих соотношений эквивалентны (формулы (1) или (3) являются аналогом записи определяющих соотношений упругости по Коши, а формула (5) - по Грину). Далее в рамках нелинейной теории упругости рассматриваем гиперупругие материалы, так как только для этих материалов гарантируется сохранение потенциальной энергии внутренних сил на замкнутых путях деформирования в пространстве компонент тензора деформаций. В дальнейшем ограничимся использованием только лагранжевых тензоров напряжений и деформаций.

Пусть - пара сопряженных лагранжевых тензоров напряжений и деформаций. Определяющие соотношения гиперупругости записываются в виде [4, 5]

(6)

где - скалярнозначная тензорная функция, которая называется удельной потенциальной энергией деформаций. Для изотропного гиперупругого материала эта функция должна зависеть только от трех главных инвариантов тензора деформаций , , или, эквивалентно, от трех главных удлинений , , [2, 4]

(7)

Рассмотрим три пары сопряженных лагранжевых тензоров напряжений и деформаций [7]

(8)

где - правый тензор деформаций Коши - Грина, - тензор деформаций Грина - Лагранжа, - правый тензор логарифмических деформаций (правый тензор деформаций Генки), - второй тензор напряжений Пиолы - Кирхгофа, - тензор напряжений Нолла. Отметим, что третья пара тензоров в (8) сопряжена только для изотропной гиперупругой среды (в общем случае в этой паре вместо тензора напряжений Нолла должен стоять правый тензор напряжений Генки [2, 7]).

С помощью каждой из трех пар тензоров в (8) строятся определяющие соотношения известных моделей гиперупругих материалов:

* Муни - Ривлина (, что соответствует несжимаемому материалу) [1,4]

(9)

где , - константы Муни - Ривлина;

* Кирхгофа - Сен-Венана [8]

(10)

* Генки [2]

(11)

Определяющие соотношения моделей Кирхгофа - Сен-Венана (10) и Генки (11) обобщают закон Гука на случай больших деформаций гиперупругой среды [7].

Для использования определяющих соотношений (9) - (11) в уравнениях гиперупругости требуется для каждой из моделей материалов для заданного закона движения (, - радиусы-векторы материальных частиц в отсчетной и текущей конфигурациях, - параметр деформирования) [1] определить второй тензор напряжений Пиолы - Кирхгофа и тензор упругости - тензор четвертого порядка, связывающий скорости изменения тензоров и в соотношении

(12)

Для определяющих соотношений изотропных материалов Муни - Ривлина и Кирхгофа - Сен-Венана эти выражения даны в [4] (для модели материала Кирхгофа - Сен-Венана ), а для модели Генки - в [2, 3].

Эксперименты по одноосному деформированию образцов из дуотана. На рис. 1 приведены данные эксперимента по одноосному деформированию образца из материала дуотан (Duothan QA965) (нанесены точками). Режиму сжатия соответствуют данные при , а режиму растяжения - при (, где - длины рабочей части образца в отсчетной и деформированной конфигурациях соответственно). Модуль Юнга определялся для модели материала Генки, имеем МПа. Осредненный коэффициент Пуассона по растяжению и сжатию (при обработке в логарифмических деформациях) имеет значение , т.е. материал дуотан - почти несжимаемый. Получены следующие значения констант для модели материала Муни - Ривлина: МПа, МПа. Зависимости напряжения от параметра для рассматриваемых моделей материалов приведены сплошными кривыми: кривая 1 соответствует обработке экспериментальных данных в рамках модели материала Муни - Ривлина, кривая 2 - Генки, кривая 3 - Кирхгофа - Сен-Венана. На рис. 1, а, 1, б приведены, соответственно, кривые зависимости продольных компонент тензоров напряжений Нолла и инженерных напряжений (первого тензора напряжений Пиолы - Кирхгофа) от параметра . Видно, что модель Кирхгофа - Сен-Венана описывает деформирование дуотана только в пределах нескольких процентов, а модели Генки и Муни - Ривлина - порядка 50% инженерной деформации.

Экспериментальные исследования и компьютерное моделирование кручения образцов из дуотана. Компьютерное моделирование деформирования эластомеров проводилось с использованием пакета MSC. Marc 2012 [1]. В этом пакете модели изотропных гиперупругих материалов Муни - Ривлина и Кирхгофа - Сен-Венана имеются в библиотеке стандартных моделей материалов, а модель материала Генки введена в пакет MSC. Marc с помощью пользовательской программы hypela2.f [3].

На рис. 2 приведены результаты натурных экспериментов и компьютерного моделирования кручения образца из дуотана. В эксперименте закручивание образца проводилось с помощью стальных выступов, приклеенных к его торцам, а в компьютерном моделировании - с использованием двух жестких поверхностей, приклеенных к торцам конечно-элементной модели. Кручение образца, как в эксперименте, так и в расчете, осуществлялось заданным углом закручивания таким образом, чтобы угловая скорость закручивания торца была постоянной. Отметим, что модель Муни - Ривлина лучше приближает зависимость закручивающего момента от погонного угла закручивания по сравнению с моделью Генки. Штриховая вертикальная линия на рис. 2, а соответствует такому значению угла закручивания, когда инженерные деформации в образце не превышают 50% (т.е. когда допустимо использование модели материала Генки). Отметим, что в допустимых пределах использования модели материала Генки (так же как и для модели материала Муни - Ривлина) графики для моментов на рис. 2, а близки к данным эксперимента, модель материала Кирхгофа - Сен-Венана работоспособна только для малых углов закручивания.

Экспериментальные исследования и компьютерное моделирование растяжения листов из полиуретана с разрезом. Проведено сравнение экспериментальных данных по растяжению листа из полиуретана [9] с результатами расчетов с использованием моделей материалов Муни - Ривлина и Генки в пакете MSC. Marc (рис. 3). В процессе деформирования первоначально прямолинейный разрез постепенно раскрывался, превращаясь в эллипс. На рис. 3, б приведена расшифровка результатов поляризационно-оптического эксперимента с получением полей напряжений и деформаций по всей области образца. Обнаружена сдвижка максимального главного напряжения от вершины разреза вглубь образца. Также следует отметить смещение максимального значения контурных напряжений от геометрических источников концентрации (вершины разреза) на некоторое расстояние от них вдоль контура эллипса, что можно визуально наблюдать по картинам полос. Однако в численных экспериментах, выполненных с использованием обеих моделей материалов, максимальное главное напряжение отмечено на границе (в точке, принадлежащей эллипсу). Причиной такого расхождения данных эксперимента и компьютерного моделирования могут быть микроповреждения в вершине разреза в эксперименте, приводящие к релаксации напряжений в окрестности этой вершины, что и является причиной сдвижки максимального главного напряжения от границы эллипса внутрь листа. Это предположение основывается на том, что в расчетах получены нереально большие значения максимального главного напряжения на границе выреза.

Библиографический список

дуотан эластомер деформирование

MARC Users Guide. Vol. A: Theory and Users Information. - Santa Ana (CA): MSC. Software Corporation, 2012.

Коробейников С.Н., Олейников А.А. Лагранжева формулировка определяющих соотношений гиперупругого материала Генки // Дальневосточный математический журнал. 2011. Т. 11, №2. С. 155-180.

Коробейников С.Н., Олейников А.А., Ларичкин А.Ю., Бабичев А.В., Алёхин В.В. Численная реализация лагранжевой формулировки определяющих соотношений изотропного гиперупругого материала Генки // Дальневосточный математический журнал (принято в печать).

Коробейников С.Н. Нелинейное деформирование твердых тел. Новосибирск: Изд-во СО РАН, 2000.

Korobeynikov S.N. Objective tensor rates and applications in formulation of hyperelastic relations // Journal of Elasticity. 2008. V. 93, No 2. P. 105-140.

Korobeynikov S.N. Families of continuous spin tensors and applications in continuum mechanics // Acta Mechanica. 2011. Vol. 216, No. 1-4. P. 301-332.

Hill R. Aspects of invariance in solid mechanics // Advances in Applied Mechanics; V. 18. C.-S. Yih (Ed.). New York: Academic Press, 1978. P. 1-75.

Curnier A. Computational Methods in Solid Mechanics. Dordrecht: Kluwer Academic Publ., 1994.

Албаут Г.Н., Харинова Н.В. Концентрация напряжений и деформаций в эластомерах // Надежность и долговечность машин и сооружений: междунар. науч.-техн. сб. - Киев, 2006. - №26. - С. 135-141.

Размещено на Allbest.ru