Механические и фрикционные свойства биоэластомеров. Часть 1: описание релаксационных зависимостей кожи человека при растяжении Текст научной статьи по специальности «Физика»

УДК 531/534: [57+61]

МЕХАНИЧЕСКИЕ И ФРИКЦИОННЫЕ СВОЙСТВА БИОЭЛАСТОМЕРОВ. ЧАСТЬ 1: ОПИСАНИЕ РЕЛАКСАЦИОННЫХ ЗАВИСИМОСТЕЙ КОЖИ ЧЕЛОВЕКА ПРИ РАСТЯЖЕНИИ

С.Б. Анфиногенов1, М.Ф. Курек2, С.В. Шилько1, Д.А. Черноус1

1 Институт механики металлополимерных систем им. В.А. Белого Национальной академии наук Беларуси, Беларусь, 246050, Гомель, ул. Кирова, 32а, e-mail: [email protected]

2 Гомельский государственный медицинский университет, Беларусь, 246000, Гомель, ул. Ланге, 5

Аннотация. Рассмотрена проблема количественной оценки деформационнопрочностных характеристик эластомеров биологического происхождения на примере кожи человека. Проведены статические механические испытания, на основе которых получены диаграммы растяжения и релаксации образцов кожи из области живота. Представлена разработанная авторами математическая модель нелинейного вязкоупругого деформирования кожи, основанная на использовании упругого потенциала и экспоненциального ядра релаксации. Вместо обычно применяемого упругого потенциала Трелоара использован потенциал Огдена. Предложена методика идентификации параметров разработанной модели. Проведено сопоставление результатов моделирования с экспериментальными диаграммами растяжения и релаксации, а также с результатами известной расчетной методики.

Ключевые слова: биоэластомеры, кожа, диаграмма растяжения, упругий

потенциал, ядро релаксации, тензор напряжений, кратность удлинения.

Введение

Природные эластомеры (биоэластомеры) весьма интересны как объект изучения в медико-биологическом отношении (например, в пластической хирургии), а кроме того, необходимы в целях создания функциональных материалов, имитирующих мышцы, кожу, сосуды и другие мягкие ткани. В частности, анализ напряженно-деформированного состояния и трения кожи человека позволяет описывать и прогнозировать эффективность и безопасность процедур пластической хирургии при типичных воздействиях (растяжении, стягивании швом и сдавливании бандажами, а также вибромеханической стимуляции). На основе ранее проведенных исследований [5-7, 9] показана возможность рассмотрения кожи как изотропного несжимаемого

нелинейного вязкоупругого материала, для которого связь компонент тензоров

напряжений и деформации определяется через удельную потенциальную энергию деформации (упругий потенциал) ч:

Т« = *, (1)

дС

0 0

Р = Т-УЯ, Т = УЯ • Р, (2)

© Анфиногенов С.Б., Курек М.Ф., Шилько С.В., Черноус Д.А., 2008

09806267

где Тх - второй тензор напряжений Пиола-Кирхгоффа; С - тензор деформации Коши-Грина; Р - тензор напряжений Пиола, определяемый в отсчетной конфигурации; Т -

0

тензор напряжений Коши; У Я - тензор-градиент вектора места [4]. Несмотря на выраженный нелинейный характер деформирования мягких биотканей типа кожи, для описания их вязкоупругого поведения может быть использовано относительно простое (двух- или четырехпараметрическое) ядро релаксации [6, 7, 9].

В настоящее время для анализа деформирования кожи человека общепризнанным считается использование упругого потенциала Муни-Ривлина или Трелоара [5-7, 9]. Однако в работах [5-7, 9], посвященных биомеханическому моделированию кожи, нет сопоставления результатов использования потенциала Муни-Ривлина с данными, полученными на основе более сложных потенциалов [2]. Это не позволяет сделать однозначный вывод о том, что потенциал Муни-Ривлина является оптимальным с точки зрения адекватного описания процесса деформирования кожного покрова человека.

Целью настоящей работы является совершенствование математической модели деформирования кожи для уточненного описания ее напряженно-деформированного состояния.

Математическая модель

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

Т У = -рО -1 + 2

где ТУ - второй тензор напряжений Пиола-Кирхгоффа для упругого материала без учета вязкости; О - мера деформации Коши-Грина; Е - единичный тензор; р -давление, обеспечивающее при произвольном напряженном состоянии выполнение условия несжимаемости материала; ч - удельная потенциальная энергия деформации (упругий потенциал); 11,12 - инварианты меры Коши-Грина.

Наряду с потенциалом Муни-Ривлина ч = С\(/1 -3) + С2(12 -3) в последнее время для описания нелинейного деформирования эластомеров часто используют потенциал Огдена [3, 8]. В упрощенной форме данный потенциал определяется соотношением

ч = (хч + Ха2 +Хаз - з), (4)

а

где ц, а - параметры упругости материала; Хг - главные кратности удлинения. Параметр ц соответствует мгновенному (в момент приложения нагрузки) модулю сдвига материала. Параметр а характеризует степень нелинейности процесса деформирования материала. Потенциал Трелоара [9] является частным случаем потенциала Муни-Ривлина при С2 = 0 и потенциала Огдена (4) при а = 2.

Для несжимаемого материала независимыми являются только две главные

кратности удлинения Х1 и Х2 ( X 3 = 1Х-^1). Независимые кратности удлинения связаны

с инвариантами меры Коши-Грина

дч

ал

+/1

дч

д/

Е-

2

дч

д/2

О

(3)

I, - х1 + х2> + х,2^22, 12 - Х12 +Х22 +Х1Х2>- (5)

Используя потенциал (4) в выражении (3) с учетом (5), для тензора напряжений Ту можно получить следующее выражение:

ту - -р°+ /Х42ц,^ К - ха+2 )Е+(ха+2 х4г -ха*2 х4 )(е - о)]. (6)

а (х1 -х 2)

Тензор напряжений (6) здесь определен без учета вязкости материала. Для описания вязкоупругого деформирования использовали методику, описанную, в частности, в работах [8, 9]. В соответствии с данной методикой второй тензор напряжений Пиола-Кирхгоффа вязкоупругого материала определяется выражением

(і - х)

т

дх , (7)

где 5, т - параметры вязкоупругости. Параметр т соответствует времени релаксации материала, параметр 5 характеризует длительный модуль сдвига цдл - д(1 - 5) , а тензор

Ту в (7) задается соотношением (6). Таким образом, функции (6) и (7) определяют связь напряжений и деформаций для рассматриваемого материала - кожи человека. Аналогичная математическая модель использовалась в работе [8] для характеристики механических свойств скелетных мышц. Отличие расчетной методики (6) и (7) от

модели, представленной в работе [8], заключается в использовании для тензора Ту аналитического выражения (6) вместо соотношения (1), поскольку соотношение (1) для несжимаемого материала нуждается в дополнительных уточнениях и трудоемких математических выкладках. Кроме того, приведенное в работе [8] выражение для

ТХ

тх ^ дм - дм

Т =(1 -5) *+80^ехр

І - х

дх

не удовлетворяет условию размерности (второе слагаемое в правой части имеет размерность Па-с).

т

Идентификация вязкоупругих характеристик

В соответствии с математической моделью, основанной на соотношениях (6) и (7), механическое поведение кожи человека характеризуется четырьмя параметрами: ц, а, 5, т. Эти характеристики могут быть найдены по результатам стандартных испытаний образцов кожи на растяжение, реализующих одноосное напряженное состояние. Используя функцию (6), для единственной ненулевой компоненты второго тензора Пиола-Кирхгоффа при одноосном напряженном состоянии получим

2

Т/ = -|дГ21 г -

а

л/г

(8)

Используя соотношения (7) и (2), для мнимого осевого напряжения получим

а (і ) = у = г (і )

ТІ (і) - 51Т1 (х) ехр|^- -—Х|йХ

(9)

Здесь Е - продольная сила, действующая на образец; Я - начальная (до начала

деформирования) площадь поперечного сечения образца. Функция (,) в (9)

определяется уравнением (8).

Параметры разработанной математической модели бит можно определить независимо от параметров ц и а по экспериментальной зависимости мнимого осевого напряжения от времени в процессе релаксации, когда соотношение (9) примет вид

( I \

(10)

Здесь а0 - начальное (на момент начала процесса релаксации) мнимое осевое напряжение. Функция (10) позволяет определить параметры б и т по двум характерным точкам на экспериментальной зависимости а (^). Для этого необходимо измерить напряжение аот, соответствующее завершению процесса релаксации (дальнейшее уменьшение напряжения пренебрежимо мало), и момент времени ^0,5, соответствующий

напряжению а0,5 = 2(а0 + ада). При известных а«, и ^0,5 параметры вязкоупругости

определяются по формулам

8 = 1 -^ т = Г051и2. (11)

а 0

Параметры ц и а предложенной математической модели будем определять по диаграмме растяжения образцов кожи. При заданной скорости растяжения V главная

V

кратность удлинения X зависит от времени следующим образом: Х(/) = 1 +— /. Здесь

Ь

Ь = 100 мм - длина рабочей части образца. Используя данное выражение для X в формулах (8) и (9), получим расчетную временную зависимость атеор(0 мнимого осевого напряжения при растяжении образца с постоянной скоростью. Произведя

замену I = Ь (х-1), установим расчетную зависимость напряжения от кратности V

удлинения атеор(Х). Параметры математической модели ц и а будем выбирать такими, чтобы расчетная зависимость атеор(Х) максимально точно соответствовала экспериментальной диаграмме растяжения аэкс(Х). Используя метод наименьших квадратов [1], будем искать такие ц и а, при которых функция

N

В(й, а) = 2[а,„ (X,)-а (X,)]2 (12)

2=1

принимает наименьшее значение. В выражении (12) N - число точек на экспериментальной диаграмме растяжения. В свою очередь, минимизация функции (12) сводится к решению системы нелинейных уравнений

^ = 0, ^ = 0. (13)

ф 5а

Решение данной системы с учетом (8), (9), (12) осуществлено в программной среде ЫШкСаё 2001. В выражении (9) использовались параметры б и т, определенные ранее по диаграммам релаксации и формулам (11).

Пример идентификации вязкоупругих характеристик

Изложенная методика идентификации была использована для определения вязкоупругих характеристик кожи области живота. Образцы для испытаний представляли собой прямоугольные лоскуты кожи толщиной 0,5-2,0 мм и шириной 12 мм, взятые от трупов мужчин и женщин в возрасте от 63 до 85 лет не позднее 24 часов после смерти. Длина рабочей части образца составляла 100 мм. Нагружение образцов и построение экспериментальных зависимостей осуществлялось на испытательной машине Instron 5567 в Институте механики металлополимерных систем им. В. А. Белого Национальной академии наук Беларуси.

Испытания на релаксацию проводились при трех значениях начального напряжения (750 кПа, 150 кПа, 50 кПа). Для каждого значения а0 получено 10 экспериментальных зависимостей. На рис. 1 представлены усредненные зависимости относительного мнимого напряжения (отношение текущего напряжения а(^ к начальному а0) от времени.

По формулам (11) для каждого значения а0 определяли параметры вязкоупругости б и т. Сопоставляя расчетную временную зависимость относительного напряжения (10) с экспериментальными диаграммами релаксации, вычисляли экспериментальные коэффициенты Фишера Еэкс [1] для каждого значения а0. Коэффициент Фишера определяется как отношение дисперсии адекватности (средний квадрат отклонения расчетного значения напряжения от среднего экспериментального) к дисперсии воспроизводимости (средняя по экспериментальным точкам дисперсия измерения напряжения). Можно отметить (табл. 1), что коэффициент Фишера не превышает 0,03. Данное наблюдение позволяет [1] сделать вывод о правомерности использования аналитического соотношения (10) для описания процесса релаксации.

В табл. 1 также приведены значения параметров б и т. Относительная погрешность параметра б совпадает с относительной статистической погрешностью измерения напряжения а«,, а погрешность времени релаксации - со статистической погрешностью измерения момента времени ^,5. При использовании испытательной машины 1т^оп 5567 приборная погрешность пренебрежимо мала (0,5% для силы и

0,3% для момента времени). Как следует из данных табл. 1, диапазоны варьирования параметров б и т перекрываются. Расхождение средних значений параметров не превышает наименьшей абсолютной погрешности. Данное наблюдение подтверждает возможность использования для изучаемого эластомера (кожи человека) допущения о независимости параметров ядра релаксации от значений деформации и напряжения. Путем объединения диапазонов варьирования параметров б и т для разных напряжений а0 определяются независящие от напряжения значения бс и тс (см. табл. 1).

Для идентификации параметров д и а были получены экспериментальные зависимости мнимого осевого напряжения от кратности удлинения образца при двух скоростях растяжения (50 и 250 мм/мин). Для каждой скорости получено по 10 экспериментальных диаграмм. На рис. 2 представлены усредненные зависимости

аэкс(Х).

Значения параметров ц и а, полученные в результате минимизации функции (12) для каждой скорости V, приведены в табл. 2. Там же представлены коэффициенты Фишера, вычисленные при сопоставлении расчетной зависимости атеор(Х) с усредненной экспериментальной диаграммой. Если для описания деформирования исследуемого материала используется потенциал Трелоара [6], параметр а равен 2.

0,2-1--------1--------1--------1--------1---------1--------1--------1

0 5 10 15 20 25 30 35 *, С

Рис. 1. Зависимость напряжения от времени (релаксация) для кожи. Экспериментальные данные: о - при ст0 = 750 кПа; □ - при ст0 = 150 кПа; А - при ст0 = 50 кПа. Пунктирная кривая построена по формуле (10) при значениях параметров 5С и тс (см. табл. 1)

а, МПа

1,0 1,1 1,2 1,3 1,4 1,5 1,6 X

Рис. 2. Диаграммы растяжения кожи. Числа у кривых соответствуют скорости растяжения в мм/мин. Пунктирная кривая построена по формулам (8) и (9) при значениях параметров

5с, тс (см. табл. 1), цс, ас (табл. 2)

По диаграмме растяжения определяется только начальный модуль сдвига ц. В табл. 2 в скобках приведены значения ц, полученные в результате минимизации функции (12) при а = 2, а также значения коэффициентов Фишера для расчетных зависимостей Отеор(Х), полученных при использовании потенциала Трелоара.

Данные табл. 2 позволяют сделать вывод о том, что математическая модель, основанная на использовании потенциала Огдена, с приемлемой (Рэкс < 0,03) точностью описывает процесс деформирования кожи человека в одноосном напряженном состоянии. При использовании потенциала Трелоара погрешность прогноза

Таблица 1

Идентификация параметров модели по диаграмме релаксации_____________

с0, кПа F экс 5 т, мс 5с тс, мс

50 0,011 0,61±0,10 1,93±0,42

150 0,023 0,59±0,09 1,77±0,28 0,61±0,10 1,94±0,45

750 0,009 0,62±0,04 2,06±0,33

Идентификация параметров модели по диаг Таблица рамме растяжения

V, мм/мин F J экс ц, кПа а Цс, кПа ас

50 0,015 (0,981) 342,61±54,46 (8205±1304) 12,42±1,97

250 0,024 (1,321) 337,54±64,11 (6018±1143) 11,52±2,19

напряжений в образце значительно (более чем в 50 раз) превышает погрешность при использовании потенциала Огдена. Относительная погрешность определения параметров ц и а принята равной средней (по N = 200 точкам экспериментальной диаграммы растяжения) относительной погрешности измерения напряжений. Из табл. 2 следует, что диапазоны варьирования параметров математической модели, определенных при различных скоростях растяжения, перекрываются. Данное обстоятельство позволяет объединить диапазоны варьирования параметров ц и а для различных скоростей v и определить независящие от скорости деформирования значения цс и ас (см. табл. 2).

Заключение

1. В результате сопоставления расчетных и экспериментальных диаграмм растяжения и релаксации образцов кожи в области живота установлено, что совместное использование упругого потенциала и простого экспоненциального ядра релаксации позволяет с приемлемой точностью описать процессы деформирования исследуемой биоткани.

2. Показана правомерность допущения о независимости параметров ядра релаксации от значений деформации и напряжения.

3. Потенциал Огдена обеспечивает более точные расчетные оценки напряжений в образце, чем широко используемый в настоящее время для характеризации механического поведения кожи потенциал Трелоара.

Список литературы

1. Андронов, Н.М. Теория вероятности и математическая статистика / Н.М. Андронов, Е.А. Копытов, Л.Я. Гринглаз. - СПб.: Питер, 2004. - 178 с.

2. Гамлицкий, Ю.А. Упругий потенциал наполненных резин / Ю.А. Гамлицкий, В.И. Мудрук, М.В. Швачич // Каучук и резина. - 2002. - № 3. - С. 29-39.

3. Гаришин, О. К. Исследование структурных напряжений в дисперсно наполненных эластомерных нанокомпозитах I О.К. Гаришин, С.Н. Лебедев II Механика композиционных материалов и конструкций. - 2006. - Т. 12, № 3. - С. 2S9-299.

4. Лурье, А.И. Нелинейная теория упругости I А.И. Лурье. - М.: Наука, 19S0. - 312 с.

3. Федоров, А.Е. О применении теории вязкоупругости в эстетической хирургии I А.Е. Федоров, В. А. Лохов II Российский журнал биомеханики. - 2003. - Т. 7, № 4. - С. 32-43.

6. Федоров, А.Е. О механических свойствах кожи человека I А.Е. Федоров, В.А. Самарцев, Т.А. Кириллова II Российский журнал биомеханики. - 2006. - Т. 10, № 2. - С. 29-42.

7. Федоров, А.Е. Моделирование поведения кожи человека при больших деформациях I А.Е. Федоров, А.А. Адамов II Российский журнал биомеханики. - 2007. - Т. 11, № 1. - С. 76-S3.

5. Bosboom, E.M.H. Passive transverse mechanical properties of skeletal muscle under in vivo compression I E.M.H. Bosboom, M.K.C. Hesselink, C.W.J. Oomens, C.V.C. Bouten, M.R. Drost, F.P.T. Baaijens II Journal of Biomechanics. - 2001. - Vol. 34. - P. 1363-136S.

9. Tham, L.M. Cupping: From a biomechanical perspective I L.M. Tham, H.P. Lee, C. Lu II Journal of Biomechanics. - 2006. - Vol. 39. - P. 21S3-2193.

MECHANICAL AND FRICTIONAL PROPERTIES OF BIOELASTOMERS. PART 1: DESCRIPTION OF HUMAN SKIN RELAXATION UNDER TENSION

S.B. Anfinogenov, M.F. Kurek, S.V. Shilko, D.A. Chernous (Gomel, Belarus)

The problem of strain and strength parameters estimation in human skin as biological elastomer is considered. The static mechanical tests of the skin from stomach zone have been carried out and tension and relaxation diagrams have been obtained. The mathematical model of nonlinear viscoelastic human skin deformation based on elastic potential and exponential relaxation function has been developed. The Ogden’s potential has been used instead of the Treloar’s elastic potential. The technique of developed model identification has been proposed. The comparison of results with experimental tension and relaxation diagrams as well as results of known technique has been made.

Key words: bioelastomers, skin, tension diagram, elastic potential, relaxation function, stress tensor, elongation rate.

Получено 07 июля 2008