CC BY
29
УДК 531(075)
Вестник СПбГУ. Сер. 1. Т. 2(60). 2015. Вып. 1
ОСЕСИММЕТРИЧНЫЕ ДЕФОРМАЦИИ
ОРТОТРОПНОГО СФЕРИЧЕСКОГО СЛОЯ
ПОД ДЕЙСТВИЕМ НОРМАЛЬНОГО ДАВЛЕНИЯ*
С. М. Бауэр, А. Л. Смирнов
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7-9
В статье рассматривается задача о деформации упругого трансверсально-изотропного слоя под действием нормального давления в трехмерной постановке. Получены асимптотические формулы, описывающие поведение мягкого слоя, то есть такого, у которого жесткость в трансверсальном направлении много меньше, чем в тангенциальном. Проведено сравнение точных и асимптотических решений, исследовано влияние параметров на величину относительного изменения толщины слоя. Библиогр. 10 назв. Ил. 5.
Ключевые слова: сферический слой, задача Ламе, трансверсальная изотропия.
1. Введение и постановка задачи. Рассматривается трехмерная задача о деформации трансверсально-изотропного сферического слоя, находящегося под действием нормального давления. Такая модель может быть использована, например, для описания изменения напряженно-деформированного состояния внешней оболочки глаза при введении внутриглазных инъекций. Для изотропного сферического слоя эта задача, известная как задача Ламе, описана, например, в [1]. Для трансверсально-изотропного слоя аналитическое решение было получено в [4, 9], но формулы, описывающие деформации, оказываются довольно громоздки и для анализа влияния параметров слоя на напряженно-деформированное состояния слоя удобнее использовать асимптотическое решение задачи.
Рассматривается деформация сферического слоя с внутренним радиусом К\ и внешним радиусом Д2 (^1 < Р < ^2), находящегося под действием внутреннего и внешнего нормальных давлений Р1 и Р2 соответственно. Положение точки оболочки описывается сферическими координатами: р — радиальная координата, у — координата в меридиональном направлении, в — координата в окружном направлении. Считаем, что материал трансверсально изотропен, то есть он изотропен по сферической поверхности (в тангенциальной плоскости ув) , тогда Е^ = Ед = Е, Ер = Е1, у^д = V, ЩР = У<рр = vi, С = с^в = 2(1+,/)' <~*рч> = <~*рв = Здесь Ер, Е^, Ев — модули Юнга, vpд, vlpд, vp¡f — коэффициенты Пуассона, Ор1р, , Орв — модули сдвига.
Из соображений симметрии, как и в [1], полагаем, что напряжения зависят только от радиальной координаты, причем
&дд (р) = арв = = арр =0
и перемещения точек происходят только в направлении р. В этом случае уравнения равновесия приводятся к одному уравнению
^РР | 1 / \ п
~Ир~ ~р рр ~ ~ =
* Работа выполнена при финансовой поддержке СПбГУ (тема НИР 6.15.820.2013) и РФФИ (гранты №13-01-00523 А, 15-01-06311 А). Часть результатов исследования была включена в доклад [8].
а граничные условия имеют вид
°рр(Лг) = -Рг, г = 1, 2.
Таким образом, после перехода к перемещениям уравнение равновесия принимает вид [4]
<(р) + ~р<(р) ~ 2Щ1-1)р>т{р) = (1)
а граничные условия таковы:
2ЕЕ^1ад(Дг) + Е2(1 - V= -Р.Д, (Е - Е^ - 2Е^2) г = 1, 2.
Для трансверсально изотропного материала ограничения на упругие постоянные приведены в [7]:
причем для материалов с V > —1, то есть для всех известных, < 1-
2. Асимптотическое интегрирование уравнения равновесия. Рассмотрим случай, когда жесткость слоя в трансверсальном направлении много меньше, чем в тангенциальном, то есть Е1 = е2Е. При этом в силу ограничений на коэффициент Пуассона VI = еа^, где е ^ 1, ^ « 1, а ^ 1. Положим а =1. Тогда уравнение (1) принимает вид
е2^0'(р)р2 + 2е2рш0(р) - (а + бе)адо(р) = 0, (2)
2 и 2и1
где а = и о = — , а граничные условия —
е.М^) + е2дад'(Л-)Л- = -ср^, г = 1, 2, (3)
где рг = Рг/Е, с = (1 - V - 2^2), / = 2^, д =1 - V. Следуя [5], решение уравнения (2) будем искать в виде
ю(р) = е* ^ *(*)<« (и,0 (р) + £и;1 (р) + . ..). (4)
После подстановки (4) в (2), собирая члены с одинаковыми степенями е, получим последовательность уравнений для определения А(р), ад0(р), т 1 (р),... вида
адо(р)(-а + р2А2(р)) =0,
2р2Ч(р)А'(р) + адо(р)(-Ь + 2рА'(р) - р2А''(р)) = 0, (5)
Решение первого уравнения (5) дает
А (р) =
р
а второго — 92
Таким образом, двучленное представление решения имеет вид
w(р)= Ар71 + Вр72, (6)
где 71,2 = + \ — 0' Константы Аъ В определяются из граничных условий
(3). Подставим одночленное приближение в граничные условия и решим полученную систему относительно А и В. Получим
_ С Д1<?1Р1 - И2С2р2
£ Р+(С\-С |) ' сС1С2(Д2С1Р2 - К2С2Р1)
В = -
где = = Д^/е.
Формула (6) описывает смещение элемента слоя вдоль нормали. В частности, она позволяет оценить изменение толщины слоя в процессе деформации. Действительно,
_ и>(Д2) - ц;(Д1) Д2 — Д1
или с учетом (6)
д = с^Д^+С, - Р-Сх) +Р2Р2(Р-С2 - Р+Сг)
е 1 + С2))(Д2-Д1) ' И
3. Примеры и анализ влияния параметров. Дальнейшие упрощения формулы (7) связаны с учетом порядка отношения величин Д1 и Д2. Рассмотрим два случая: толстый слой (Д1 ^ Д2) и тонкий слой (Д1 « Д2).
Для толстого слоя положим £ = Д1/Д2 (£ ^ 1). При этом 6*1 ^ 6 ив первом приближении
д, +Р2-Р1"
(8)
В качестве примеров будем рассматривать упругий слой с параметрами склеральной оболочки человеческого глаза [9, 10]: V = 0.48, р1 = р * 133.3Па, р2 = 0Па, Е = 14МПа, V* = 0.5, причем считаем, что внутриглазное давление р = 20 мм рт. ст.
На рис. 1, а представлены графики зависимости ДЛ от е при £ = 0.7. Синяя линяя— точное решение, красная линия — формула (8). Интересно, что даже в случае относительно большого значения параметра £ численные и асимптотические результаты хорошо согласуются, а при £ < 0.7 практически совпадают.
Случай тонкого слоя более сложный, так как в разложение входят два конкурирующих малых параметра. Разложение точного решения задачи с одним малым параметром п, связанным с толщиной слоя, выполнено в работах [2, 3]. Положим Д1 = (1 — п)Д2 (п =1 — £ < 1), тогда
= с Р1(1 - 77)(^+ - ^-(1 - п)х) +р2(Р- - ^+(1 - пГ) ег/ + (1-77)ж)
Рис. 1. Сравнение асимптотических и численных значений относительного изменения толщины при изменении параметра относительной жесткости слоя.
где я = у/а/е. Будем считать, что <С 1, то есть <С е, тогда, сохраняя главные члены, получим
Д/г = —
с Р1(1 - - 1 - 1]я)) +р2(Р- - 1 - ??*:))
2ег)
(1+ пк/2). (9)
На рис. 1, Ь представлены графики зависимости ДЛ от е для тонкого слоя п = 0.04. Синяя линяя — точное решение, красная линия — формула (9). Отметим хорошее совпадение графиков при «больших» значениях е.
На рис. 2 представлены графики зависимости ДЛ от е (точные значения) при изменении величины Д1/Д2 в диапазоне [0.46,0.96]. С уменьшением толщины слоя модуль величины ДЛ быстро возрастает, достигая 10% при п = 0.01 (тонкая оболочка) и е = 0.1
Рис. 2. Зависимость относительного изменения толщины от параметра относительной жесткости слоя для разных значений параметра относительной толщины слоя.
На рис. 3 изображены графики зависимости ДЛ от £ = Д1/Д2 (точные значения) при изменении величины е в диапазоне [0.05,0.3]. Как и следовало ожидать, с ростом е изменение толщины становится меньше. При стремлении толщины слоя к нулю величина изменения толщины неограниченно убывает. Действительно, при малых п главный член в формуле (9) имеет вид
ДЛ = -
с (Р1 -Р-)
2 £1]
(10)
и ДЛ с уменьшением толщины слоя убывает по гиперболическому закону.
На рис. 4, объединяющем рис. 2 и 3, представлена зависимость ДЛ от е и £ =
Д1/ Д2
Ак
о.ооси
- 0.005
- 0.010 - 0.015
0.020
Рис. 3. Зависимость относительного изменения толщины от относительной толщины слоя для разных значений параметра относительной жесткости слоя.
0.3 Е
Рис. 4. Зависимость относительного изменения толщины слоя от параметров относительной жесткости и толщины слоя.
Наконец, на рис.5 представлена зависимость ДЛ от V! при Д1/Д2 = 0.96 для различных значений е: е = 0.03 — красная линия, е = 0.1 — зеленая линия, е = 0.3 — синяя линия. С увеличением Vl относительное изменение толщины линейно растет, причем этот эффект тем больше, чем меньше е.
Д, 0.000 - 0.002
- 0.004
- 0.006 - 0.008
0.010 0.012 - 0.014
Рис. 5. Зависимость относительного изменения толщины слоя от коэффициента Пуассона.
4. Заключение. Асимптотические соотношения показывают, что при деформации сферического слоя под действием внутреннего давления относительная толщина слоя меняется тем больше, чем меньше толщина слоя.
Изменение толщины склеральной оболочки при увеличении внутриглазного давления, как один из результатов представленного моделирования, согласуется с имеющимися данными о том, что при глаукоме, глазном заболевании, сопровождающимся, как правило, повышенным уровнем внутриглазного давления, также наблюдается уменьшение толщины склеры [6].
Зная относительное изменение толщины сферического слоя при увеличении внутреннего давления, например, при введении внутриглазных инъекций, можно оценить отношение радиального и тангенциального модулей упругости склеральной оболочки.
Литература
1. Атанацкович Т., Гуран А. Лекции по теории упругости. СПб., 2003, 400 с.
2. Бауэр С. М., Воронкова Е. Б. Модели неклассических теорий анизотропных оболочек в задачах биомеханики глаза // Вестн. Нижегородского ун-та им. Н. И. Лобачевского, №4, часть 5. С. 2000— 2001.
3. Бауэр С. М., Воронкова Е. Б. Неклассические теории анизотропных оболочек в задачах о деформации трансверсально-изотропных сферических и цилиндрических слоев под действием нормального давления // Вестн. С.-Петерб. ун-та. Сер. 1. 2011. Вып.3. С. 86—93.
4. Бауэр С.М., Замураев Л. А., Котляр К.Е. Модель трансверсально-изотропного сферического слоя для расчета изменения внутриглазного давления при интрасклеральных инъекциях // Российский журнал биомеханики. 2006. Т. 10. №2. С. 43—49.
5. Бауэр С.М., Смирнов А. Л., Товстик П.Е., Филиппов С. Б. Асимптотические методы в механике твердого тела. М.; Ижевск, 2007. 360 с.
6. Волков В. В. Глаукома при псевдонормальном давлении. Руководство для врачей. М.: Медицина, 2001.
7. Родионова В. А., Титаев Б. Ф., Черных К. Ф. Прикладная теория анизотропных пластин и оболочек. СПб.: Изд-во С.-Петерб. ун-та. 1996. 278 с.
8. Bauer S. M., Smirnov A. L. Deformation of the orthotropic spherical layer under normal pressure // Book of Abstracts XLII Summer School—Conference "Advanced Problems in Mechanics". St. Petersburg (Repino), Russia 2014. P. 38-39.
9. Kotliar K., Maier M., Bauer S., Feucht N., Lohmann C., Lanzl I. Effect of intravitreal injections and volume changes on intraocular pressure: clinical results and biomechanical model // Acta Ophthalmologica Scandinavica. Vol.85. Issue 7. 2007. P. 777-781.
10. Kotliar K., Maier M., Bauer S., Feucht N., Lohmann C., Lanzl I. Intravitreal injection of triamcinolone acetonide and intraocular pressure: Author's reply // Acta Ophthalmologica. Vol. 86. Issue 6. 2008. P. 692-693.
Статья поступила в редакцию 23 октября 2014 г.
Сведения об авторах
Бауэр Светлана Михайловна —доктор физико-математических наук, профессор; s_bauer@mail.ru
Смирнов Андрей Леонидович — кандидат физико-математических наук, доцент; a_l_smirnov@mail.ru
DEFORMATIONS OF THE ORTHOTROPIC SPHERICAL LAYER UNDER NORMAL PRESSURE
Svetlana M. Bauer, Andrei L. Smirnov
St.Petersburg State University, Universitetskaya nab., 7-9, St.Petersburg, 199034, Russian Federation; s_bauer@mail.ru, a_l_smirnov@mail.ru
The deformation of the orthotropic spherical layer under normal pressure applied on the outer and inner surfaces is analyzed. The layer is assumed to be slightly othotropic, it permits to apply asymptotic methods. The equations of zeroth and first approximations are derived. For the shell, which is much softer in the transverse direction than in the tangential plane, one gets singularly perturbed boundary value problem.
Solving this problem in the zeroth approximation the asymptotic formula for the change of the relative layer thickness under normal pressure is obtained. Also the effect of Poisson ratio and the layer thickness on the deformation is studied. For the cases of the thick and thin layers the last formula may be simplified. The asymptotic results well agree with the exact solution. The developed formulas are used in analysis of the scleral shell under intraocular pressure and may also be used in solution of the inverse problem, i.e. in analysis of the stress-strain state of a human eye under injection. The solution of the problem helps to estimate the mechanical parameters of the sclera, i.e. to find the ratio of the tangential and transversal Young modules using clinical data for the sclera thickness change. Refs 10. Figs 5.
Keywords: spherical layer, Lame problem, transverse isotropy.
References
1. Atanackovic T. M., Guran A., Theory of Elasticity for Scientists and Engineers (Birkhauser, Basel, 2000).
2. Bauer S. M., Voronkova E. B., "Non-classical anisotropic shell theories in ocular biomechanics", Vestnik of Lobachevsky State University of Nizhni Novgorod N4, Part 5, 2000—2001 [in Russian].
3. Bauer S.M., Voronkova E. B., "Nonclassic theories of anisotropic shells in the problems on the deformation of transversely isotropic spherical and cylindrical layers under normal pressure", Vestn. S.-Peterb. Univ. Ser. 1 Issue 3, 86-93 (2011).
4. Bauer S. M., Zamuraev L. A., Kotliar K. E., "Model of the transversely isotropic spherical layer for estimation of intraocular pressure changes after intravitreal injections", Russ. J. Biomech. 10(2), 41-47 (2006).
5. Bauer S.M., Smirnov A.L., Tovstik P. E., Filippov S. B., Asymptotic methods in mechanics of Solids (R&C Dynamics, Moscow, Izhevsk, 2007) [in Russian].
6. Volkov V. V., Glaucoma with pseudonormal pressure (Medicine, Moscow, 2001) [in Russian].
7. Rodionova V. A., Titaev V. F., Chernykh K. F., Applied Theory of Anisotropic Plates and Shells, (Izd. St.Peterb. Univ., St.Petersburg, 1996) [in Russian].
8. Bauer S. M., Smirnov A. L., "Deformation of the orthotropic spherical layer under normal pressure", Book of Abstracts XLII Summer School-Conference "Advanced Problems in Mechanics" 3839 (St.Petersburg (Repino), Russia, 2014).
9. Kotliar K., Maier M., Bauer S., Feucht N., Lohmann C., Lanzl I., "Effect of intravitreal injections and volume changes on intraocular pressure: clinical results and biomechanical model", Acta Ophthalmologica Scandinavica 85, Issue 7, 777-781 (2007).
10. Kotliar K., Maier M., Bauer S., Feucht N., Lohmann C., Lanzl I., "Intravitreal injection of triamcinolone acetonide and intraocular pressure: Author's reply", Acta Ophthalmologica 86, Issue 6, 692-693 (2008).
