Numerical study of a thermomechanical behavior of a crystallizing polymer medium with regard to finite deformations

Abstract


The paper presents phenomenological constitutive relations for a crystallizing polymer medium, obtained in the framework of nonlinear solid mechanics. The relationships are based on the representation of the medium in the form of a composition of a molten and fully crystallized material, taking into account the history of a continuous nucleation and deformation of a new phase in the temperature range of phase transformations. A general formulation is carried out regarding the evolutionary boundary-value problem of the nonlinear mechanics of polymer materials under phase transitions with the use of the proposed constitutive relations. Algorithms are considered that are aimed at a numerical realization of the thermokinetic problem and determining the stress-strain state of the solidifying system for the case of a plane deformed state. A linearization procedure is developed - which is convenient for constructing numerical algorithms to solve the set evolutionary boundary value problems - using the assumption of the proximity of each intermediate configuration to the current one, which corresponds to the procedure for superimposing small deformations of crystallized particles on finite deformations of a crystallizing medium. The linearization procedure for the initial formulation of mechanics problem is realized taking into account temperature and structural deformations. A numerical algorithm aimed at solving a flat boundary value problem is developed and implemented with the aim of investigating the evolution of the stress-strain state in a polymer structure. A method is proposed and implemented to construct a discrete analogue of boundary value problems, based on the use of the Galerkin method with a selection basis functions with a compact support by the finite element method. In this case, the increments of the displacement functions at the current time step are taken as nodal unknowns. The regularities of a shell type defect’s formation in the crystallizing polymer cylinder are established.

Full Text

Введение Процесс производства значительной части полимерных изделий и конструкций содержит этап охлаждения материала, сопровождающийся его кристаллизацией. Технологические и остаточные напряжения, формирующиеся вследствие температурной и деформационной неоднородностей, могут приводить к снижению эксплуатационных качеств изделий и даже их разрушению еще на стадии изготовления. Поэтому актуальной является проблема описания термомеханических процессов в полимерных материалах в условиях фазового перехода. Модели поведения таких сред должны включать определяющие соотношения, единым образом описывающие связь тензоров напряжений и деформаций в широком интервале температур, включающем диапазон фазовых превращений. Для реальных процессов и изделий решение задач описания закономерностей формирования полей технологических напряжений возможно только численными методами. Поэтому важное значение приобретает построение удобного в реализации численного алгоритма решения краевой задачи термомеханики кристаллизующейся полимерной среды. Хорошо известны модели поведения отверждающихся полимерных материалов, предложенные в разное время в работах [1-3]. Общая черта этих моделей - стремление получить возможность простой оценки остаточных напряжений. Ряд авторов при построении моделей термомеханического поведения для кристаллизующегося полимера предполагает упругое поведение материала и малость деформаций. Например, в работе [4], в рамках указанных предположений, с использованием феноменологического подхода разработана система определяющих уравнений, описывающих термомеханическое поведение аморфно-кристаллических полимеров. В работах [5-9] моделирование термомеханического поведения полимера проведено с учетом больших деформаций. При этом многие из предлагаемых моделей поведения полимера основаны на упруговязкопластическом подходе [7-10], разработаны для конкретных полимерных материалов [5, 10] или определенных видов напряженно-деформированного состояния. Например, в работе [6] разработана модель описания механического поведения полимера в диапазоне температур, включающем процессы стеклования, для случаев одноосного нагружения и плоского напряженного состояния. В [8] рассматривается развитие упруговязкопластической теории для моделирования зависимости отклика аморфных полимеров от скорости деформирования и температуры с учетом больших деформаций. В последнее время интерес к вопросам описания деформационного поведения полимеров при фазовых превращениях вызван все более широким использованием их в качестве материалов с памятью формы, поскольку в основном эффекты памяти формы в полимерах связаны с фазовыми или релаксационными переходами. Различным аспектам деформационного поведения полимеров, связанным с исследованием эффектов памяти формы, посвящено большое число публикаций [11-15]. Большинство предлагаемых в них моделей используют гипотезу о малости деформаций и предназначены для описания поведения полимера в достаточно узком температурном диапазоне и моделирования конкретного механизма памяти формы и не могут описывать другие термомеханические эффекты. В настоящей работе рассматривается построение феноменологических определяющих соотношения для кристаллизующегося полимера с учетом больших деформаций, непрерывным образом отражающих связь тензоров напряжений и деформаций в широком диапазоне изменения температур. Соотношения основаны на представлении среды в виде композиции расплавленного и полностью закристаллизованного материала с учетом истории непрерывного зарождения и деформирования новой фазы в интервале температур фазовых превращений. Также в работе рассматривается модельная задача о зарождении раковин в кристаллизующейся среде, демонстрирующая возможности применения предлагаемых определяющих соотношений. 1. Постановка краевой термомеханической задачи для кристаллизующегося полимера. В предположении, что технологические напряжения не влияют на температуру и протекание процесса кристаллизации, возможно разделение задачи на две самостоятельные задачи: 1) теплокинетическая задача - задача по определению полей распределения температуры и степени кристаллизации; 2) краевая задача по определению напряженно-деформированного состояния (НДС) затвердевающей системы. Интенсивное тепловыделение при кристаллизации и теплоотдача в окружающую среду, низкая теплопроводность полимерных материалов, зависимость процесса кристаллизации от температуры приводят в процессе изготовления изделий из полимеров к существенно неоднородным полям температуры и степени кристаллизации. Прогнозирование температурного и конверсионного полей в ходе процесса может быть осуществлено в рамках известной теплокинетической модели, согласно которой постановка задачи включает: - нестационарное уравнение теплопроводности с переменными коэффициентами (1) - уравнение кинетики неизотермической кристаллизации [16, 17] (2) - начальные условия , (3) - граничные условия по температуре, где - удельная теплоемкость; - плотность материала; - коэффициент теплопроводности; - температура; - время; - радиус-вектор произвольной точки тела; - скорость выделения тепла в ходе протекания реакции кристаллизации; - тепловой эффект реакции кристаллизации; - степень кристаллизации материала, которая в ходе процесса принимает значения в пределах от нуля до некоторого предельного значения, зависящего от температуры и представляющего собой равновесную степень кристаллизации ; - универсальная газовая постоянная; - температура плавления; - макрокинетические константы, определяемые экспериментально; - область затвердевающего полимера. Уравнения (1) и (2) образуют связанную систему уравнений и решаются совместно. Результатом решения теплокинетической задачи являются пространственно-временные поля распределения температуры и степени кристаллизации, с учетом которых решается задача определения НДС изделия. Для постановки краевой задачи относительно НДС затвердевающей системы необходимы физические соотношения, которые в силу специфики задачи непрерывным образом отражают связь тензоров напряжений и деформаций в широком диапазоне изменения температур. Эти соотношения были получены с учетом больших деформаций в рамках следующих гипотез и положений: - кристаллизующаяся полимерная система в каждый момент времени в каждой точке системы рассматривается как смесь расплава полимера и полностью закристаллизовавшегося продукта, относительная доля которых в общем объеме характеризуется величиной степени кристаллизации ; - материал в процессе превращения и дальнейшего деформирования остается макроизотропным; - для описания механического поведения расплава полимера и закристаллизовавшегося полимера примем упругий закон. Для описания поведения среды при конечных деформациях будем использовать понятия начальной (соответствующей моменту времени ) и текущей (соответствующей моменту времени ) конфигураций [18]. Пусть начальная конфигурация определяется радиусом-вектором , а текущая конфигурация - радиусом-вектором , где - обобщенные координаты. Пусть в момент времени закристаллизовалась доля полимера, характеризуемая приращением степени кристаллизации . Градиент места из конфигурации, в которой появилась доля кристаллической фазы , в текущую конфигурацию может быть записан следующим образом [18]: (4) где - градиент места из конфигурации, соответствующей моменту времени , в конфигурацию, соответствующую моменту времени . Для градиента места из конфигурации, соответствующей моменту времени и кристаллизации доли полимера, характеризуемой приращением степени кристаллизации , в текущую конфигурацию справедливо следующее выражение: Рассмотрим п-й дискретный шаг по превращению : . Учитывая, что , ,…, , выразим все градиенты места через градиент места, связывающий начальную и текущую конфигурации: , ,….., (5) где - градиент места из начальной конфигурации в конфигурацию, соответствующую приращению степени кристаллизации . Для получения определяющих соотношений необходимо записать выражение для упругого потенциала системы [19], зависящего от тензора меры деформаций Коши-Грина. Упругий потенциал кристаллизующейся системы может быть представлен в виде суммы упругого потенциала для аморфной фазы и упругого потенциала для кристаллической фазы , (6) где - доля потенциала, которую привносит в общий потенциал часть кристаллической фазы, соответствующая приращению степени кристаллизации ; - тензор меры деформаций Коши-Грина [19], соответствующий деформированию от начального до текущего момента времени; - тензор меры деформаций Коши-Грина, соответствующий деформированию от момента времени до текущего момента времени. После осуществления предельного перехода в соотношении (6) по степени кристаллизации выражение для упругого потенциала можно записать в виде (7) Используя в качестве упругого потенциала двухконстантный потенциал Пенга-Ландела [20], выражение связи тензора Пиолы-Кирхгофа второго рода с тензором меры деформаций Коши-Грина записываем следующим образом: (8) где , , - механические характеристики аморфной и закристаллизовавшейся части материала соответственно; - первый и третий главные инварианты тензора второго ранга; - метрический тензор. Вариационная постановка краевой квазистатической задачи термомеханики кристаллизующейся среды может быть получена с использованием процедуры Бубнова-Галеркина [20] и в случае отсутствия внешних силовых воздействий принимает вид , (9) где - тензор деформаций Коши-Грина; - начальный объем тела. К достоинствам рассматриваемой постановки стоит отнести то, что она сформулирована в исходной конфигурации, что существенно упрощает построение численного алгоритма расчета. При этом полагается нелинейность геометрических соотношений, обусловленная необходимостью учета больших деформаций. Сложность получения решения задачи в форме (9) обусловлена как нелинейностью постановки, так и проблемами, связанными с необходимостью учета возникающих температурных и структурных деформаций. При решении модельной задачи будем полагать, что зарождение полости в материале происходит при достижении объемной деформации в окрестности рассматриваемой точки некоторого предельного значения. Для придания поведению рассматриваемого тела большей реалистичности предположим, что предельная величина деформации разрушения несколько варьируется по объему. Разрушение будем моделировать занижением в окрестности рассматриваемой точки модуля кристаллической фазы. С учетом неосесимметричного характера зарождения полостей решение механической задачи будет рассматриваться в декартовой системе координат для случая плоского деформированного состояния. В данном случае компоненты тензоров и конкретизируются следующим образом: , , где , , , ; - компоненты тензора деформаций, отражающие переход из промежуточной конфигурации (соответствующей времени ) в текущую конфигурацию (соответствующую времени t). 2. Выполнение процедуры линеаризации исходной постановки с учетом структурных и температурных деформаций Рассматриваемый в данной работе алгоритм предполагает линеаризацию постановки (9), что позволяет свести процедуру получения решения исходной нелинейной задачи к решению последовательности линейных упругих краевых задач в рамках теории малых деформаций. Линеаризация геометрических соотношений может быть выполнена путем разложения исходного движения тела на совокупность малых движений, вызываемых малыми приращениями внешних воздействий. Согласно идее наложения малых деформаций на конечные [21, 22], движение тела из начальной конфигурации в конечную представимо в виде суперпозиции движений тела через промежуточные конфигурации. Таким образом, перемещения тела в произвольный момент времени представимы в виде суммы накопленных и малых перемещений, описывающих переход из промежуточной конфигурации в конечную. Рассмотрим движение тела, приводящее его к моменту времени из начальной недеформированной конфигурации в конфигурацию , через промежуточную конфигурацию . Положение некоторой точки в начальной конфигурации определим радиусом-вектором , в конфигурациях и радиусом-векторами и соответственно. Пусть и - векторы перемещений, переводящие тело из конфигурации в и из в соответственно. Каждый из этих векторов является функцией обобщенных координат . Предполагается, что конфигурация близка к конечной и, соответственно, перемещения малы. В конфигурациях , введем обычным образом локальные основные и сопряженные базисы где и - контравариантные компоненты метрического тензора в начальной и промежуточной конфигурациях соответственно. Операторы Гамильтона в конфигурациях и определены следующим образом: Тогда тензоры-градиенты места, описывающие переход из начальной конфигурации в промежуточную, из промежуточной в конечную и из начальной в конечную, могут быть записаны в виде Градиент места может быть представлен через тензоры и следующим образом: С учетом этого тензор меры деформаций Коши-Грина записывается в виде Тензор меры деформаций Коши-Грина , описывающий переход из в , с учетом малости может быть преобразован к виду , где - тензор малых деформаций. Тензор деформации Коши-Грина, описывающий переход из начальной конфигурации в конечную, записывается в форме и может быть представлен в виде где - вектор перемещений из начальной конфигурации в промежуточную; u - вектор малых перемещений из промежуточной в конечную, операторы Гамильтона определены относительно исходной и промежуточной конфигураций, тензор соответствует деформациям, возникающим при переходе из начальной конфигурации в промежуточную, - из промежуточной в конечную. Для плоского деформированного состояния в декартовой системе координат операторы Гамильтона записываются следующим образом: , где - компоненты векторов и u соответственно, . Окончательно для плоского деформированного состояния в декартовой системе координат линеаризованные геометрические соотношения имеют вид Процедура линеаризации физических соотношений состоит в разложении соотношения (8) в ряд Тейлора с учетом только линейных слагаемых. В результате линеаризованные физические соотношения, определяющие компоненты тензора Пиолы-Кирхгофа второго рода принимают вид (10) где, например, , Величины - приращения компонент тензора деформаций Коши-Грина на рассматриваемом шаге. При построении физических соотношений в случае использования нелинейных уравнений состояния существенным является вопрос учета температурных деформаций материала. Ряд авторов [23, 24] предлагают использовать введение линейного слагаемого, зависящего от приращения температуры, аналогично тому, как это принято в случае использования линейных соотношений. В работе [25] используется подход, основанный на линеаризации физических соотношений и предположении о малости температурных деформаций на рассматриваемом переходе. Очевидно, данная методика обладает большей общностью, хотя сразу стоит отметить существенные сложности реализации, связанные с необходимостью выделения упругой механической составляющей кинематики процесса. В настоящей работе рассматривается вопрос построения численного алгоритма решения задач определенного класса (определения технологических напряжений в ходе процесса получения полимерных изделий), в которых деформации, являясь значительными по величине, все-таки остаются умеренными (порядка 10 %). С учетом ограниченности деформаций в рассматриваемом процессе линеаризованные физические соотношения в матричном виде могут быть записаны следующим образом: где - тензор суммы малых температурных и структурных деформаций, возникающих при переходе из промежуточной конфигурации в конечную; - для аморфной фазы; - для кристаллической фазы; - коэффициенты температурного расширения для аморфной и кристаллической фаз соответственно; - коэффициент структурной усадки; - единичный тензор; тензор описывает переход из конфигурации, соответствующей моменту времени τ, в промежуточную; , - тензоры четвертого ранга «касательных жесткостей» аморфной и кристаллической фаз, компоненты которых зависят от и . С учетом выполненных процедур линеаризации выражение (9) может быть преобразовано к виду (11) Таким образом, определение эволюции поля перемещений с течением времени сводится к решению последовательности линеаризованных задач относительно приращения вектора перемещений на временных шагах. 3. Разработка и реализация численных процедур решения поставленных краевых задач Получение решения краевых задач механики для тел сложной пространственной конфигурации невозможно без использования численных методик. Среди существующих численных подходов в настоящее время наиболее широко применяется метод конечных элементов [26], который и был выбран для построения дискретного аналога рассматриваемой задачи. Применение конечно-элементного подхода к аппроксимации геометрических и физических соотношений приводит к получению линейных алгебраических соотношений, которые могут быть представлены в матричной форме [26]: , где - вектор узловых неизвестных (приращений перемещений на рассматриваемом шаге); - матрица градиентов элемента; - вектор приращения температурных и структурных деформаций; , - матрицы «касательных» жесткостей аморфной и кристаллической фаз материала. Последнее соотношение может быть преобразовано к виду , где . В матричной форме соотношение (8) можно записать как , где означает суммирование по элементам; - локальная матрица жесткости элемента; , - локальные векторы узловых сил, обусловленные наличием накопленных напряжений, температурными и структурными деформациями: , , . В результате сборки глобальной системы для всего ансамбля конечных элементов получим , (12) где - глобальная матрица системы, обладающая симметрией и зависящая от вектора накопленных перемещений ; - глобальные векторы «фиктивных» усилий, обусловленные накопленными к рассматриваемому шагу напряжениями и возникающими на рассматриваемом переходе температурными и структурными деформациями. Таким образом, получаем последовательность линейных задач, сформулированных относительно приращений узловых перемещений. Оператор левой части соотношения (12) требует переопределения на каждом шаге, что обусловливает существенную вычислительную нагрузку при решении задач большой размерности. Для случая решения двумерных задач рассматривается построение конечно-элементного алгебраического аналога с использованием трехузлового конечного элемента с линейной аппроксимацией узловых неизвестных, в качестве которых приняты приращения перемещений. Все соотношения получены в базисе исходной конфигурации. 4. Решение модельной задачи описания эволюции дефекта типа раковины в кристаллизующемся полимерном цилиндре Рассматривается задача о протекании процесса кристаллизации и формировании дефекта типа раковины в полимерном цилиндре. Изначально материал цилиндра (полиэтилен) имеет температуру 430 К, что превышает температуру плавления. Цилиндр диаметром 0,2 м охлаждается по внешнему контуру водой, температура которой составляла 293 К. При решении теплокинетической задачи на границе тела задаются условия конвективного теплообмена. Коэффициент теплоотдачи на поверхности принимается равным 550 Вт/(м2×К). В процессе охлаждения в материале протекают процессы кристаллизации, сопровождающиеся выделением тепла. Краевая задача нестационарной теплопроводности решалась методом конечных элементов. Для решения задачи кинетики кристаллизации использовался разностный метод Рунге-Кутты [27]. Учитывалась зависимость тепломеханических свойств от температуры. Значения теплофизических и механических параметров, зависимость удельной теплоемкости, коэффициента теплопроводности и плотности полиэтилена от температуры взяты из литературных источников [28-30]. При решении задачи определения НДС конструкции задавались кинематические граничные условия, гарантирующие отсутствие смещений тела как абсолютно жесткого. Эволюция полей температуры и степени кристаллизации представлена на рис. 1. Закономерности эволюции раковины в цилиндре иллюстрируются рис. 2. а б в Рис. 1. Распределение температуры и степени кристаллизации для моментов времени: а - 7800 с; б - 22800 с; в - 37800 с Fig. 1. Distribution of temperature and crystallization degree for time instants: а - 7800 s; b - 22800 s; c - 37800 s а б в г Рис. 2. Образование раковины для моментов времени: а - 16550 с; б - 17550 с; в - 18300с; г - 19550 с Fig. 2. Formation of the shell for time instants: а - 16550 s; b - 17550 s; c - 18300 s; d - 19550 s Заключение Таким образом, выполнена формулировка двумерных наследственных краевых задач нелинейной термомеханики кристаллизующихся полимерных материалов в рамках теории конечных деформаций. Предложен и реализован метод линеаризации исходной постановки с учетом структурных и температурных деформаций. Разработан и реализован численный алгоритм решения поставленной плоской краевой задачи с целью исследования особенностей эволюции напряженно-деформированного состояния в полимерной конструкции. Получено численное решение модельной задачи о формировании дефекта типа раковины в кристаллизующейся полимерной пластине.

About the authors

R G Kulikov

Perm National Research Polytechnic University

T G Kulikova

Perm National Research Polytechnic University

O Yu Smetannikov

Perm National Research Polytechnic University

References

  1. Production stresses during the solidification of cylindrical articles formed from polymer composite materials / V.N. Korotkov, R.A. Turusov, E.A. Dzhavadyan, B.A. Rozenberg // Mechanics of Composite Materials. - 1986. - Vol. 22. - No. 1. - P. 99-103.
  2. Shaffer B.W., Lewitsky M. Termoelastic Constitutive Equation for Chemically Hardering Materials // Journal of Appl. Mech. - 1974. - Vol. 41. - No. - P. 652-657.
  3. Lewitsky M., Shaffer B.W. Residual Thermal Stresses in a Solid Sphere Cast From a Termosetting Material // Journal of Appl. Mech. - 1975. - Vol. 42. - No. 9. - P. 651-655.
  4. Golotina L.A., Matveenko V.P., Shardakov I.N. Analysis of deformation process characteristics in amorphous-crystalline polymers // Mechanics of Solids. - 2012. - Vol. 47. - P. 634-640. doi: 10.3103/S0025654412060040
  5. Boyce M.C., Socrate S., Llana P.G. Constitutive Model for the Finite Deformation Stress-Strain Behavior of Poly(Ethylene Terephthalate) above the Glass Transition // Polymer. - 2000. - Vol. 41. - No. 6. - P. 2183-2201.
  6. Rebecca B. Dupaix, Mary C. Boyce. Constitutive modeling of the finite strain behavior of amorphous polymers in and above the glass transition // Mechanics of Materials. - 2007. - No. 39. - P. 39-52.
  7. Modeling and Validation of the Large Deformation Inelastic Response of Amorphous Polymers over a Wide Range of Temperatures and Strain-Rates / J. Richeton, S. Ahzi, K.S. Vecchio [et al.] // Int. J. Solids Struct. - 2007. - No. 44. - P. 7938-7954.
  8. A Thermo-Mechanically Coupled Theory for Large Deformations of Amorphous Polymers. Part I: Formulation / L. Anand, N.M. Ames, V. Srivastava, S.A. Chester // Int. J. Plast. - 2009. - Vol. 25. - No 8. - P. 1474-1494.
  9. A thermo-mechanically-coupled large-deformation theory for amorphous polymers in a temperature range which spans their glass transition / V. Srivastava, S.A. Chester, N.M. Ames, L. Anand // International Journal of Plasticity. - 2010. - No. 26. - P. 1138-1182.
  10. Varghese A.G., Batra R.C., Constitutive Equations for Thermomechanical Deformations of Glassy Polymers // Int. J. Solids Struct. - 2009. - Vol. 46. - No. 22. - P. 4079-4094.
  11. Finite Deformation Thermo-Mechanical Behavior of Thermally Shape Memory Polymers / H.J. Qi, T.D. Nguyen, F. Castroa et al. // J. Mech. Phys. Solids. - 2008. - Vol. 56. - No. 5. - P. 1730-1751.
  12. Liu C., Qin H., Mather P.T. Review of progress in shape-memory polymers // J. Mater. Chem. - 2007. - Vol. 17. - P. 1543-1558.
  13. Thermomechanics of shape memory polymers: uniaxial experiments and constitutive modeling / Y.P. Liu, K. Gall, M.L. Dunn, A.R. Greenberg, J. Diani // Int. J. Plast. - 2006. - Vol. 22. - P. 279-313.
  14. Dietsch B., Tong T. A review: features and benefits of shape memory polymers (SMPs) // J. Adv. Mater. - 2007. - Vol. 39. - P. 3-12.
  15. Gunes I.S., Jana S.C. Shape memory polymers and their nanocomposites: a review of science and technology of new multifunctional materials // J. Nanosci. Nanotechnol. - 2008. - Vol. 8. - P. 1616-1637.
  16. Simulation of Thermomechanical Processen in Crystallising Polymer / I.N. Shardakov, V.P. Matveyenko, N.V. Pistsov, V.P. Beghishev // Polym. Eng. and Sci. - 1997. - Vol. 37. - No. 8. - P. 1270-1279.
  17. Degree of crystalline structure of polymer obtained from melt at various cooling rates / A.Ya. Malkin, S.A. Bolgov, V.P. Begishev, O.S. Mazalov // Journal of Engineering Physics and Thermophysics. - 1991. - Vol. 61. - No. 3. - P. 1092-1095.
  18. Новокшанов Р.С., Роговой А.А. Эволюционные определяющие соотношения для конечных вязкоупругих деформаций // Изв. РАН. Механика твердого тела. - 2005. - № 4. - С. 122-140.
  19. Лурье А.И. Нелинейная теория упругости. - М.: Наука, 1970. - 940 с.
  20. Методы прикладной вязкоупругости / А.А. Адамов, В.П. Матвеенко, Н.А. Труфанов, И.Н. Шардаков / УрО РАН. - Екатеринбург, 2003. - 411 с.
  21. Роговой А.А. Определяющие соотношения для конечных упруго-неупругих деформаций // Прикладная механика и техническая физика. - 2005. - № 5. - С. 138-149.
  22. Куликов Р.Г., Куликова Т.Г. К вопросу определения деформированного состояния кристаллизующейся полимерной среды с учетом больших деформаций // Вестник Пермского национального исследовательского политехнического университета. Механика. - Пермь, 2012. - № 1. - С. 62-72.
  23. Vikas Srivastava, Shawn A. Chester, Lallit Anand. Thermally actuated shape-memory polymers: Experiments, theory and numerical simulations // Journal of the Mechanics and Physics of Solids. - 2010. - No. 58. - P. 1100-1124.
  24. Адамов А.А. Исследование и моделирование нестационарного термомеханического поведения вязкоупругих резиноподобных материалов и элементов конструкций при конечных деформациях: дис. … д-ра физ.-мат. наук. - Пермь, 2004.
  25. Rogovoy A.A. Formalized approach to construction of the state equations for complex media under finite deformations // Continuum Mechanics and Thermodynamics. - 2012. - No. 24. - P. 81-114.
  26. Зенкевич О. Метод конечных элементов в технике. - М.: Мир, 1975. - 541с.
  27. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. - М.: Лаборатория базовых знаний, 2001. - 630 с.
  28. Теплофизические и реологические характеристики полимеров: справ. / под ред. Ю.С. Липатова - Киев: Наук. думка, 1977. - 244 с.
  29. Пивень А.Н., Гречанная Н.А., Чернобыльский Н.И. Теплофизические свойства полимерных материалов. - Киев: Вища школа, 1976. - 180 с.
  30. Нильсен Л. Механические свойства полимеров и полимерных композиций. - М.: Химия, 1978. - 312 с.

Statistics

Views

Abstract - 50

PDF (Russian) - 25

Cited-By


PlumX


Copyright (c) 2018 Kulikov R.G., Kulikova T.G., Smetannikov O.Y.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies