CC BY
94
ISSN 1812-5123. Российский журнал биомеханики. 2014. Т. 18, № 2: 140-146
УДК 531/534: [57+61]
Российский
Журнал
Биомеханики
www.biomech.ru
МОДЕЛИРОВАНИЕ РОСТОВЫХ ДЕФОРМАЦИЙ В ЖИВЫХ СИСТЕМАХ
В.А. Лохов
Кафедра теоретической механики и биомеханики Пермского национального исследовательского политехнического университета, Россия, 614990, Пермь, Комсомольский проспект, 29, e-mail: valeriy.lokhov@yandex.ru
Аннотация. В работе использован метод декомпозиции собственной деформации в линейно-упругом теле, позволяющий разложить собственную деформацию на две части: первая часть не вызывает напряжений в системе, а вторая часть не вызывает деформаций системы. Показано применение данного подхода в решении задач независимого управления собственными деформациями, дана постановка задачи моделирования систем с собственными деформациями. Под собственной деформацией понимается неупругая деформация любой природы, например температурная деформация, деформация фазовых переходов, деформация роста и перестройки живой ткани, пьезоэлектрическая деформация и т.д. Предложенный подход применен для моделирования ростовых деформаций при лечении врожденной расщелины твердого нёба. Данная патология при неудачном лечении может приводить к серьезным косметическим дефектам, а также к отставанию умственного развития ребенка. Получено определяющее соотношение для ростовой деформации, проведен его анализ и показано, что при определенных граничных условиях рост не будет вызывать напряжений в системе, что существенно упрощает процедуру решения задачи о накоплении ростовой деформации с течением времени. Этот результат важен для планирования оптимального ортопедического лечения данной патологии, которое приводит к задаче независимого управления деформациями системы посредством ростовых деформаций. Определяющее соотношение получено для изотропного однородного тела, так как костная ткань ребенка не имеет сформированной структуры и изотропна. Использование свойств изотропии дает возможность задействовать в определяющем соотношении для ростовых деформаций только два скалярных параметра. Это позволяет учесть основные особенности роста и провести эксперименты по определению этих параметров.
Ключевые слова: управление формой, ростовая деформация, врожденная расщелина твердого нёба.
Введение
Проблемы создания желаемого распределения напряжений или формы имеют ключевое значение в научной литературе и широко обсуждаются на научных конференциях, в частности в медицине. Оптимальные распределения напряжений могут значительно уменьшить время заживления переломов и вызвать требуемый рост и перестройку живой ткани.
Вывод единой теории для решения таких задач может помочь в разработке новых методов медицинского лечения. В работе предлагается подход к решению таких задач на основе концепции собственной деформации. Это общий термин для неупругих
© Лохов В. А., 2014
Лохов Валерий Александрович, к.ф.-м.н., доцент кафедры теоретической механики и биомеханики, Пермь
деформаций любой природы: деформаций теплового расширения, пьезоэлектрических деформаций, деформаций из-за фазовых переходов в сплавах с памятью формы, деформаций роста и перестройки в живых тканях и т.д. Концепция собственной деформации была первоначально введена Рейсснером и Эшелби [10, 18] и подробно обсуждается в книге [13]. Предполагается, что распределение собственной деформации можно отдельно вычислить с помощью решения соответствующей краевой задачи, например, температурная деформация зависит от найденного температурного поля, деформация фазовых переходов зависит от нагрузки и истории нагрева, пьезоэлектрическая деформация зависит от приложенного электрического напряжения и т.д. Следовательно, собственная деформация предполагается известной, и тензор полной деформации в является суммой тензоров упругой деформации ве и собственной деформации в *:
в = ве +в *, г е О, (1)
где Ос Е3 - область пространства, занимаемая линейно упругим телом.
Упругая деформация в рамках линейной теории связана с напряжениями 6 законом Гука:
ве = С---б , г еО, (2)
где С-1 - тензор четвертого ранга упругой податливости.
Проблемы, связанные с определением собственной деформации, не обсуждаются в данной работе, так как есть огромное количество работ, затрагивающих эту проблему. Тем не менее оказывается, что существуют значительные свойства полученных напряжений и деформаций, не зависящие от происхождения собственных деформаций. В качестве примера будет рассмотрено определяющее соотношение для ростовой деформации.
Данный подход позволяет сформулировать и решать задачи двух весьма важных классов: 1) создание благоприятного распределения напряжений без деформации; 2) управление формой структуры без изменений поля напряжений.
Постановка краевой задачи
Пусть исследуемое тело занимает ограниченную область О трехмерного евклидова пространства Е33. Замыкание области обозначим О, границу (которая считается достаточно гладкой ) - через £ (^ = Ои £). Граница области £ делится на две взаимно не пересекающиеся части (£ = Би и £6):
и = 0, г е £; п -6 = р, г е £а, (3)
где р - заданный вектор поверхностных сил.
Кроме того, заданы объемные силы д и собственная деформация в *. Нахождение напряжений, деформаций и перемещений требует решения соответствующей краевой задачи механики континуума. Дифференциальная постановка краевой задачи для среды с собственными деформациями приведена в различных работах (например, [14, 15]). Во многих случаях краевые задачи решаются приближенно с использованием дискретизации, например методом конечных элементов. В связи с этим для дискретных и дискретизированных систем используют обобщенную формулировку краевой задачи.
Назовем обобщенным решением задачи симметричный тензор с, который определяется законом Гука при наличии собственной деформации:
(7 =
C- (s(и) -s*), (4)
где и е ; и = 0, г е Su, и для которого имеет место соотношение (равенство
работ внешних и внутренних сил)
--s(w)dV- Jp-wdS-Jq-wdV
= 0,
S(
3 ° (5)
Vw e(íW21 (Q))\ w = 0, r e Su.
Здесь Ж,1 - функциональное пространство Соболева. Деформации s (u) и s (w) определяются геометрическим соотношением Коши, где производные понимаются в обобщенном смысле. В постановке считается, что p e(L2 (Sc)) , q e(L2 (Q)) ,
s* e (L2 (Q)) , компоненты Cijkl (i, j, k,l = 1,2,3) тензора четвертого ранга упругих
констант С считаются кусочно-непрерывными функциями координат.
Обобщенная постановка дает возможность анализа дискретных систем, для которых несправедлива дифференциальная постановка краевой задачи. В работе [1] доказана теорема, из которой следует, что для решаемой задачи обобщенное решение существует и единственно.
Из соотношений (3), (4) следует, что рассматриваемая задача моделирования распадается на две части: задача линейной теории упругости с известными объемными силами q и поверхностными силами p; задача линейной теории упругости с
собственной деформацией s *.
Заметим, что напряжения в системе, вызванные собственными деформациями, называют собственными напряжениями (eigenstresses) с *. В случае, если они самоуравновешенные (Su = 0), их называют также остаточными напряжениями
(residual stresses). Собственные напряжения отличаются от остаточных тем, что они создаются, в том числе, давлением неподвижных опор. Напряжения, вызванные только действием объемных и поверхностных сил, будем называть силовыми напряжениями сf. Можно показать, что если в теле существует собственная деформация и на тело действуют объемные и поверхностные силы, то напряжения в теле будут суммироваться:
с = сf + с*, r eQ. (6)
Если произвести разгрузку тела, т.е. снять внешние и поверхностные силы, то в теле останутся только собственные напряжения, а силовые напряжения будут равны нулю. В некоторых случаях собственная деформация может при разгрузке измениться, например, в задачах пластичности при появлении вторичной пластической деформации.
Собственные деформации, свободные от напряжений
Большое значение для дальнейшего изложения имеет понятие собственной деформации, свободной от напряжений (stress-free eigenstrain), которое было введено Эшелби [10].
Собственная деформация называется свободной от напряжений и обозначается s*, если при наложении этой деформации на упругое тело напряжения в любой точке тела равны нулю:
sU ^а = 0, г eQ. (7)
Свойства и примеры таких распределений широко обсуждаются в литературе, например в работах [2, 3, 4, 5, 7, 8, 12, 14, 15]. Классическим примером такой деформации может служить однородный нагрев свободного стержня.
Для того чтобы собственная деформация s не вызывала напряжений в любой точке тела, необходимо и достаточно, чтобы: 1) собственная деформация была совместной, т.е. для нее существует поле перемещений Эй(u1,u2,u3)(r) - такое, что
s* = (VU + UV)/2; 2) перемещения равны нулю на опорах, й(r) = 0 при r eSu. Это
утверждение доказано в работах [14, 15].
Также доказана теорема, позволяющая определять поля собственных деформаций, свободных от напряжений. Соответствующая теорема доказана в работе [16].
Теорема. Тензор собственной деформации s* не создает напряжений в любой точке данного тела, r eQ, тогда и только тогда, когда существуют объемные силы
q e(L2 (Q)) и поверхностное натяжение p e(L2 (Sa)) , которые производят такую же
упругую деформацию s , r eQ в линейно-упругом теле, т.е. s = s *, x e Q .
Описанные условия позволяют определить свойства собственной деформации без решения краевой задачи.
Анализ определяющих соотношений для деформации роста
Рост характерен для живых тканей и сопровождается изменениями в структуре и размерах тела, притоком массы. Изменение размеров тела можно описать с помощью тензора ростовой деформации sg и тензора скорости ростовой деформации g.
В технике рост, как правило, вызван присоединением частиц на поверхности тела, например в задаче о намотке. В живых тканях рост происходит за счет внутренних источников массы, связанных с движением крови.
Существует множество моделей роста живой ткани, учитывающих различные эффекты и массоперенос, но основная сложность в их использовании - недостаток экспериментальных данных для верификации параметров моделей. Есть и модели, феноменологически описывающие изменения растущей ткани, например, деформацию, которые содержат в себе намного меньше параметров, для которых можно найти экспериментальные данные. В качестве такой модели во многих работах выбрана модель, обсуждаемая в работах Хсю, Регирера, Штейна [3, 6, 9, 11]:
|g = A + B--а, (8)
где A - тензор второго ранга, описывающий генетически заложенный рост, B - тензор четвертого ранга, отражающий влияние напряжений в теле <5 на деформацию скорости
роста |g.
Соотношение (8) позволяет описать изменение размеров растущего тела, но не может описать приток массы. Однако для решения ряда задач этого достаточно, например, в задаче ортопедического лечения врожденной расщелины твердого нёба главную роль играет форма формируемого нёба.
Далее мы рассмотрим вопрос о влиянии ростовых деформаций на напряжения.
Отметим, что наличие слагаемого В• • б дает нам возможность, меняя напряжения в теле, управлять скоростью роста. В медицинской практике для создания напряжений используются ортопедические аппараты. С течением времени ростовая деформация, вычисленная по формуле (8), может вызывать собственные напряжения, которые останутся в теле после снятия ортопедического аппарата и могут дать отрицательные эффекты. Однако анализ работ [3, 9] показывает, что в аналогичных задачах на макроуровне напряжения в процессах роста не меняются.
Рассмотрим изотропный материал, для которого тензор А является шаровым тензором:
= А1 + В , (9)
где I - единичный тензор второго ранга.
При малом приращении времени Л ростовая деформация 88 может быть найдена путем умножения (9) на время Л:
В8 =| = ( А1 + В--(5) Ж. (10)
Если ростовая деформация (10) удовлетворяет условиям, указанным выше, то она не будет вызывать собственных напряжений. Мы можем увидеть, что слагаемое АМг аналогично температурной деформации для однородного нагрева и при определенных граничных условиях не вызовет собственных напряжений.
Слагаемое В • • бЛ требует более детального анализа.
Воспользуемся теоремой, приведенной выше, для этого перепишем слагаемое В • • 5Л в следующей форме:
}• •б = }• С •В6 = М • •в6, (11)
где М - тензор четвертого ранга, М = В • • С.
Далее мы упростим соотношение (11), учитывая изотропию материала. Существует три линейно-независимых изотропных тензора четвертого ранга:
Мт = 5л5> Мт = 5*5 л ±8Й5к . Поэтому запишем тензор М в виде линейной комбинации трех тензоров:
Мт = Х5л5 + Ц (5 л + 55 ) + Ц (5*5 л - 5Й5).
Теперь вычислим свертку тензора М и тензора упругой деформации в6:
(Мв6) = М1]к1 в6к1 = Ц,в6кк + ^8* +8л ) + ц (вл -вл).
^ 'V
Учитывая симметрию тензора деформации и обозначая 2ц = М, получим
(М •В6 )л = Л + Мвл. (12)
Ввиду сложности экспериментального определения параметров модели в соотношении (12) делается упрощение, аналогичное тому, что сделано в работе [17]: Х = 0.
Тогда уравнения (9) и (12) дают следующую модель для ростовой деформации:
|g = AI + Mse. (13)
Тогда, если применить теорему, приведенную выше, то видно, что слагаемое Msedt не вызывает напряжений.
Таким образом, при определенных граничных условиях ростовая деформация не
вызывает собственных напряжений. Если объемные q е (L2 и поверхностные
p g(L2 (Sct)) силы не зависят от времени, то напряжения и скорость ростовой
деформации также будут стационарными. Тогда для случая малых деформаций ростовую деформацию можно вычислять путем умножения на время лечения T :
sg =1 gT. (14)
Заключение
Использование формулы (14) позволяет решать задачу о ростовых деформациях в перемещениях, а не в приращениях, что существенно сокращает затраты на решение. Кроме того, это упрощает задачу о создании заданной ростовой деформации и заданной формы растущего нёба. Данные свойства позволят биомеханикам упростить процедуру подбора ортопедических аппаратов при лечении врожденной расщелины твердого нёба и повысить эффективность проводимого лечения.
Отметим, что для применения формулы (14) нужно обращать внимание на граничные условия, иначе ростовая деформация, соответствующая слагаемому AIdt, будет вызывать давление на опоры и собственные напряжения в системе.
Список литературы
1. Дюво Г., Лионс Х.-Л. Неравенства в механике и физике. - М.: Наука, 1980. - 480 с.
2. Кучумов А.Г., Лохов В.А., Няшин Ю.И., Менар М., Селянинов А.А. Численное решение задачи оптимизации для определения параметров установки фиксаторов спамятью формы // Российский журнал биомеханики. - 2009. - Т. 13, № 1. - С. 18-28.
3. Кизилова Н.Н., Логвенков С.А., Штейн А.А. Математическое моделирование транспортно-ростовых процессов в многофазных биологических сплошных средах // Известия РАН. Механика жидкости и газа. - 2012. - № 1. - С. 3-13.
4. Лохов В.А., Кучумов А.Г. Создание заданных усилий в фиксаторах, изготовленных из сплавов с эффектом памяти формы // Российский журнал биомеханики. - 2006. - Т. 10, № 3. - С. 41-52.
5. Лохов В.А., Няшин Ю.И., Кучумов А.Г., Менар М., Гачкевич А.Р., Будз С.Ф., Онышко А.Е. Применение материалов с эффектом памяти формы к лечению патологий зубочелюстной системы // Российский журнал биомеханики. - 2008. - Т. 12, № 4. - С. 7-17.
6. Лохов В.А., Долганова О.Ю., Няшин Ю.И. Биомеханическое моделирование эффекта сближения фрагментов твердого неба при ортопедическом лечении // Российский журнал биомеханики. -2012. - Т. 16, № 1. - С. 38-45.
7. Туктамышев В.С., Лохов В.А., Няшин Ю.И. Независимое управление напряжениями и деформациями в живых системах // Российский журнал биомеханики. - 2011. - Т. 12, № 2. - С. 69-76.
8. Федоров А.Е., Лохов В.А. О применении теории вязкоупругости в эстетической хирургии // Российский журнал биомеханики. - 2003. - Т. 7, № 4. - С. 34-46.
9. Штейн А.А., Юдина Е.Н. Математическая модель растущей растительной ткани как трехфазной деформируемой среды // Российский журнал биомеханики. - 2011. - Т. 15, № 1. - С. 42-51.
10. Eshelby J.D. The determination of the elastic field of an ellipsoidal inclusion and related problems // Proceedings of the Royal Society of London. - 1957. - Vol. 241. - P. 376-396.
11. Hsu F.-H. The influences of mechanical loads on the form of a growing elastic body // Journal of Biomechanics. - 1968. - Vol. 1. - P. 303-311.
12. Irschik H., Ziegler F. Eigenstrain without stress and static shape control of structures // AIAA Journal. -2001. - Vol. 39. - P. 1985-1990.
13. Mura T. Micromechanics of Defects in Solids. - 2nd ed. - Dordrecht: Kluwer Academic Publisher, 1991. -587 p.
14. Nyashin Y., Lokhov V., Ziegler F. Decomposition method in linear elastic problems with eigenstrain // Z. Angew. Math. Mech. - 2005. - Vol. 85. - P. 557-570.
15. Nyashin Y., Lokhov V., Ziegler F. Stress-free displacement control of structures // Acta Mechanica. -2005. - Vol. 175. - P. 45-56.
16. Lokhov V., Nyashin Y., Kiryukhin V., Ziegler F. Theorem on stress-free eigenstrain and Duhamel's analogy // Journal of Theoretical and Applied Mechanics. - 2006. - Vol. 36, № 3. - P. 35-46.
17. Masich A.G., Nyashin Y.I. Mathematical modelling of orthopedic reconstruction of childrens congenital maxillary anomaly // Russian Journal of Biomechanics. - 1999. - Vol. 3, № 1. - P. 101-109.
18. Reissner H. Eigenspannungen und Eigenspannungsquellen // Z. Angew. Math. Mech. - 1931. - Vol. 11. -P. 1-8.
MODELLING OF GROWTH STRAIN IN LIVING SYSTEMS
V.A. Lokhov (Perm, Russia)
In this paper, we use the method of decomposition of the eigenstrain in a linear elastic body allowing decompose eigenstrain into two parts: the first part does not cause stress in the system, and the second part does not cause deformation of the system. The paper illustrates the application of this approach in solving the problems of independent shape control and gives the formulation of the boundary value problem for modelling systems with eigenstrains. In this framework, we consider eigenstrain as inelastic strain of any nature, for example, thermal strain, strain due to phase transitions, growth and remodelling strains of living tissues, piezoelectric strain, etc. The proposed approach opens the way to simulate the growth strains at the treatment of congenital cleft of the palate. This pathology in case of failed treatment can lead to serious cosmetic defects and also can block mental development of children. The conducted analysis of obtained governing equation for the growth strain shows that under certain boundary conditions, the growth does not induce stresses in the system, which simplifies the solving process of the problem of the growth strain accumulation with time. This result is important for the optimal planning of orthopedic treatment of the disease, which leads to the problem of independent shape control by growth strains. The constitutive relation assumes isotropic homogeneous body since the bone of a child has no structure and isotropic. Considering the properties of isotropy allows us to use in constitutive relation for growth strains only two scalar parameters, and thus to describe the basic features of growth and to conduct experiments on the determination of these parameters.
Key words: shape control, growth strain, congenital cleft of the hard palate.
Получено 30 апреля 2014
