УДК 539.3:517.925 Вестник СПбГУ. Сер. 1. Т. 2(60). 2015. Вып. 2
ОСЕСИММЕТРИЧНАЯ ДЕФОРМАЦИЯ МЯГКОЙ АРМИРОВАННОЙ НИТЯМИ ТОРОИДАЛЬНОЙ ОБОЛОЧКИ
А. В. Кузьмин
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7-9
Рассматривается нелинейная осесимметричная деформация тороидальной оболочки под действием внутреннего давления. Оболочка армирована двумя система нитей, расположенных по параллелям и меридианам. Предполагается, что нити расположены достаточно часто, поэтому после осреднения получаем двухмерную упругую среду, которая является мягкой оболочкой.
Основная особенность мягкой оболочки заключается в том, что она не выдерживает сжимающих напряжений. В равновесном положении все меридианы оболочки растянуты, а часть параллелей может быть сжатой, что приводит к образованию складок.
Для определения деформаций и перемещений оболочки с учетом геометрической и физической нелинейностей получена система обыкновенных дифференциальных уравнений четвертого порядка. Разработан метод численного решения краевой задачи для этой системы, а также метод асимптотического интегрирования в предположении, что радиус окружности меридиана тора существенно меньше радиуса его параллели. Проведено сравнение асимптотических и численных результатов. Получены асимптотическое и численное выражения для предельного давления, при котором на оболочке исчезают складки и она становится полностью растянута. Библиогр. 3 назв. Ил. 2. Табл. 3.
Ключевые слова: анизотропная мягкая оболочка, геометрическая нелинейность, физическая нелинейность, внутреннее давление, асимптотическое решение, численное решение.
1. Введение. Рассматривается осесимметричная деформация тороидальной оболочки, образованной сворачиванием в тор прямоугольной сети из двух систем нитей. В нерастянутом положении длины нитей, расположенных по параллелям и меридианам, равны соответственно Ь и I. Внутри сети из нитей помещена податливая оболочка с газом, создающая внутреннее давление на тороидальную оболочку (рис. 1). Предполагается, что нити расположены достаточно часто, поэтому после осреднения получаем двумерную упругую среду, которая является мягкой оболочкой [1].
Основная особенность мягкой оболочки заключается в том, что она не выдерживает сжимающих напряжений. В связи с этим при ее деформации могут появляться зоны с одноосным напряженным состоянием. В частности, в рассматриваемых ниже задачах в равновесном положении все нити, идущие по меридианам, натянуты, а часть нитей, идущих по параллелям, может быть ненатянутой.
Данная работа является дополнением к работам [2, 3]. В [2] учитывалась только жёсткость нитей, которые считались линейно-упругими, в [3] — и нити, и сама оболочка считались нелинейно-упругими.
В данной работе найдены члены следующего порядка малости для асимптотических разложений, полученных в [3], для случая нелинейных соотношений упругости.
Рис. 1. Тороидальная оболочка.
2. Основные уравнения. Введём обозначения: во —длина дуги меридиана до деформации, отсчитываемая от крайнего правого положения О, 0 < во < I; в (во) — длина дуги после деформации; г(во) —радиус параллели; ¿(во) —высота параллели над точкой О; Ах (во), А2(во) —кратности удлинений меридианов и параллелей. Имеют место геометрические соотношения
= Т~"> = Д=-,
аво Д 2п ,1)
3'Г - вшб» ^-соъб 1 - 39 1 - с°8(9 ¿в ' ¿в ' Дх ¿в' Д2 г '
где Д1, Д2 —радиусы кривизны поверхности, в — угол между касательной к меридиану и вертикальным направлением.
Уравнения равновесия в проекции на касательную к меридиану и нормаль к оболочке имеют вид
где ц — внутреннее давление, а усилия Т1 и Т2 отнесены к единице длины после деформации.
Оболочку и нити считаем изготовленными из нелинейно-упругих несжимаемых однородных изотропных материалов с упругими потенциалами:
^ФДАх,А2,АЗ), г =1, 2, 3, (3)
где С — константы материала.
Примем, что нити находятся в одноосном напряженном состоянии. Тогда для нитей, идущих по меридианам, с учётом несжимаемости материала (АхА2Аз = 1) имеем А2 = А3 = А- 1/2, и потенциал принимает вид
С1Ф1(АЬ А-1/2, А-1/2) = СФх(Ах), (4)
а сила натяжения нитей равна
f ¿Ф1
Т{ = 8ХСХ—(5)
аА1
где —площадь поперечного сечения нитей до деформирования. Для нитей, идущих по параллелям,
с2Ф2(А2,А2-1/2,А2-1/2) = с2Ф2(А2), = (6)
Для участков оболочки, находящихся в условиях плоского напряженного состояния, потенциал (3) принимает вид
СзФз(А1, А2, (А1А2Г1) = СзФз(А1, А2). (7)
Истинные напряжения вычисляются по формулам
дФ3
^ = г = 1,2. (8)
Усилия, входящие в уравнения (2), равны
^ = СзЬодФз + С1М1Б1М1 С3Ы>дФ3 | С2М2Я2 ¿Ф2 (д)
Л2 дЛх Х2 ¿Лх' Лх ЗХ2 Лх ¿Х2'
где N1 и N2 — число нитей (на единицу длины в поперечном направлении) в положении до деформации для меридианов и параллелей соответственно, ко — толщина оболочки до деформации.
Если на некотором участке формула (9) даст Т2 < 0, то на этом участке реализуется одноосное напряженное состояние, и в соответствии с гипотезами мягких оболочек в системе (2) следует считать Т2 =0.
3. Безразмерные переменные. Введём безразмерные величины по формулам
(10)
z I р qR poR
* = so = fm р = Q=G^> =
2 _ TiX2 ^ _ T2Xi _ GjNjSj ¡ — — —
Go ho' Go hO Go ho ' ' ' Go'
где Go — произвольно задаваемый модуль сдвига, р < 0.5, а параметры gi характеризуют относительную жёсткость нитей по сравнению с оболочкой.
Считая основными неизвестными величины в(р), ii(y), A2(p), z(p), 0 < p < 2п, запишем систему уравнений
dd р ~ dt1 ~ .
— = — (AiA2Q - t2 cos в), —— = -pt2 sin в,
dP ti dP (11)
dA2 л -a dz \ a
—— = —pAi sin 0, — = pXi cos 0,
dp dp
где
, ,, Л N дФ3 dФ 1 дФ3 d<!2 „ r , . . ti = Ji{M,M) = 93-^7-+ 91-77-, h = 93^7- +32-TT-, h = rnax{t2,0}. (12)
dA1 dA1 дХ2 dX2
В силу симметрии задачи относительно плоскости z = 0 при отыскании периодического решения системы (11) достаточно удовлетворить граничным условиям
0(0) =0, z (0) = 0, в(п) = п, z (п) = 0. (13)
Рассмотрим частный вид упругих потенциалов в (3). Пусть функции Фi имеют
вид
ф4(а1,а2,аз) = | (af + + л§ - 3), г = 1,2,3. (14)
Тогда безразмерные усилия t1, t2 будут равны
t2 = 9з ( А2 - -щз ^ + д2 ^А2 - ц ) .
4. Численное интегрирование. Численно решаем задачу Коши для системы уравнений (11) с начальными условиями
в(0) = 0, ^ (0)= , А2(0) = Ао, г(0) = 0. (16)
Величины А°, ¿о определяем из уравнений
в(п) = п, г(п)=0. (17)
Пусть нити, идущие вдоль меридианов и параллелей, одинаковы, т.е. $1 = $2. Для удобства сравнения относительного влияния жёсткости нитей и самой оболочки будем считать, что $1 + $з = 1. Возьмём ^ = 0.3 и рассмотрим три случая:
1) $1 = 1, $з = 0; 2) $1 = $з = 1/2; 3) $1 = 0, $з = 1. (18)
В случае 1 учитывается только жёсткость нитей, в случае 2 влияние жёсткости нитей и оболочки одинаково, в случае 3 нити отсутствуют. Результаты в трёх случаях различаются незначительно, поэтому далее рассматривается случай 2, для которого результаты представлены в таблице 1.
Таблица 1. Деформация тора в зависимости от давления Q при ^ = 0.3, дх = дз = 1/2
Q Положения равновесия
Докритические Закритические
ti(0) Л2 (0) А2(тг) ¿2(7Г) ti(0) Л2(0) А2(ТГ) ¿2(7Г)
0.8000 0.231 1.166 0.563 -2.749 7.280 6.083 2.606 2.531
1.0000 0.308 1.202 0.588 -2.394 5.778 4.852 2.086 1.967
1.2000 0.394 1.240 0.613 -2.073 4.754 4.024 1.740 1.565
1.4000 0.492 1.281 0.638 -1.774 3.998 3.424 1.493 1.248
1.6000 0.607 1.329 0.666 -1.484 3.398 2.963 1.307 0.974
1.8000 0.747 1.387 0.699 -1.190 2.887 2.587 1.160 0.715
2.0000 0.936 1.468 0.741 -0.862 2.406 2.257 1.038 0.441
2.2000 1.318 1.645 0.823 -0.361 1.780 1.883 0.912 0.042
2.2232 1.523 1.748 0.865 -0.157 1.538 1.756 0.868 -0.143
Существует предельное давление Q = Qp = 2.223, при превышении которого равновесные положения отсутствуют. При Q < Qp имеются два равновесных положения— докритическое и закритическое (рис.2). Результаты для них приведены в левой и в правой частях таблицы соответственно. Для докритических положений равновесия деформации растут с ростом давления, а для закритических положений наблюдается обратная зависимость. У докритических положений равновесия вблизи точки ^ = п напряженное состояние одноосно, о чём говорит знак у величины t2(n).
5. Приближенное решение для закритических положений равновесия.
Предположим, что ^ ^ 1, и построим асимптотические разложения решений системы (11), (12), описывающих закритические положения равновесия. После замены переменных
ti = Ai = А2 = t2 = (19)
система уравнений (11), (12) принимает вид
a'(íp) = — /х<5 sin в, У((р) = — /х/3 sin í?, O'(íp) = — (/?7<3 — ladeos в), z'(<p)=¡3 cosO. (20)
а
Рис. 2. Зависимость между кратностью удлинения Л2(0) и давлением для трёх случаев соотношений жёсткостей.
2.5
2,0
t.0
0.5
0.0
\ * ч'\ — 1 --- 2 ........3
Ь/
1
ш
а(ф) = Бф(ф) -3(ф = б27-
91
в2 (ф)
92
Y 2(ф
• м3 + O(/6), Bi = 91 + 93, М3 + O(M6), Б2 = 93 + 93.
(21)
Будем искать приближенное решение системы (20), (21), удовлетворяющее граничным условиям (13). Подставляя
в(ф)
а(ф0 , 9iBi з
Б1
+
КФ)
М3, = Y(Ф)Б2 -
52 М3
Y 2(Ф)
(22)
в уравнения (20), получаем систему уравнений
¿в yQ 7Б2 cos в 91Б1 + yQ 3 92 cos в 4
— = —---А4 Н--з-А4 Н--5—А4 ,
аф Bi а а3 «Y2
da . 92 sin в 4
— = -7-У2 sm 0(j, Н--т—А4 ,
аф Y2
а sin в 91 Б1 sin в 4
Б1
"А4
"А4 ,
dy ¿ф
dz a cos в 91B1 cos в 3 Bi а2 ^
Для получения первых трёх приближений представим
к = ко + Kim + К2М2 + К3М3, где к = а, в, Y, в, z. Запишем систему уравнений нулевого приближения:
(23)
а0(ф)=0, Yo (ф)=0, в0 (ф) =
7о (<fi)Q
в1 '
)(ф) =
a0(y)cos ¿>0(у)
Bi
Получаем решение, удовлетворяющее краевым условиям во(0) = 0, во(п) = п, zo(0) = 0, zo(n) =0:
ао(ф) = ао, Yo(^)
Bi Q
, во(ф) = ф, zo(^>)
ао sin ф
Б
где ао = const — определяется из первого приближения.
а
z
Запишем систему уравнений первого приближения: d 71 (p)Q BiB2 cos p d B1B2 sin p
= ---^Q — , =--Q-'
d aosiny d ao0i(y) siny ^ ai(y) cosy
dp Bi ' dp Bi Bi
Ищем решение, удовлетворяющее краевым условиям 0i(0) = 0i(n) = 0, zi(0) = zi(n) =0.
Интегрируя второе и третье уравнения, получаем
BiB2 , ao
ai = —-—cosy» + ai, 71 = —cos у+ 61, Q Bi
где ai, bi —произвольные постоянные. Подставляя выражения для ai, Yi в первое уравнение и учитывая краевые условия, находим
Q2a0 - BiB2 . , „ B¡a0Q 61 =°-
Из четвёртого уравнения получаем
Q2Oq sin2у - 2Q2aQ(f> + ABlaiQsmif + 4BfB2p
Zl= ÍBfg '
а из условия zi(n) = 0
^/TBiB^Bi a0 =-Q-'
Запишем систему уравнений второго приближения:
d B2 Q ai Q B2 2
^2(y) = ^ + ^72^) + ^COSp--COS y,
d 3В2\/ад . d ( 2B2 аЛ .
=--2Q-¿p*{>p)= (^cosy.-jsmy,
d .. 3 B2 /ЩГ ai B2 \ . 2 \/2BiB2 , . . аг(у)cos у
29 (y>) =--4 /-cos у--, sin у--02(cp) sin у +
dtp ^ ^ 4Q V Si ^ Вг^ЩЩ) Г Q Г В i
Аналогичным образом интегрируем данную систему, краевые условия для которой 02 (0) = 02 (п) = 0 = Z2(0) = Z2(n) = 0.
Окончательно, асимптотическое решение системы (23) имеет вид
. V^BjSi В1В2 cos у \/2 BiS|(l + 3cos 2р) + 2Q3gi 2 ,
аЫ= q + --^¡ЖЖ-M +0(M'
, , Bi V2B1B2 cosy B2 cos2y 2 . Зч = ^ +-Q-At+_2Q-** +0(/i)'
. V2S1B2 sin у B2 sin 2y V^(2BiB|cos2y+ Q3fli)siny 2 , 3. =-Q-+ -AB^VBjh-" +0(M
6. Приближенное решение для докритических положений равновесия.
Существует предельное значение давления Q = при ^ > оболочка полностью растянута, при Q < Q* на ней имеются складки.
Соотношение для Q* получаем из уравнений для полностью растянутой оболочки. Предполагаем, что ¿2 > 0 при ^ € [0, п], и вводим переменные по формулам
¿1 = ма, Ах = 1 + мв, А2 = 1 + М7, ¿2 = М^, г = мС (25)
Система уравнений (11), (12) и граничные условия (13) примут вид
— = —¡л8 8И1 в, — = — Г(1 + /х/3)( 1 + у«7)<5 ~~ А4^ со;3 ,
= (1 + ¿4/3) сое <9, ^ = -(1 +
М« = (01 + £з)(1 + М/3) - 51 53
М^ = (02 + 03)(1 + М7) -
(1+ мв)2 (1 + Мв)3(1+ М7)2'
02 03
(26)
(27)
(1+ М7)2 (1+ М7)3(1+ мв)2'
0(0) = С (0) = С (п)=0, 0(п) = п. (28)
Решая систему до второго приближения, получим выражение для ¿2:
¿2(^) = ¿20 + ¿21(^)М + ¿22 (^)М2 -
Для корректности решения необходимо выполнение условия ¿2(^0 > 0, ^ € [0, п]. Отсюда получаем
6
<3* « <Зо + м<3ъ <Зо = о—г-;— (40з + 4(01 + 02)03 + 3^132),
301 + 403
31 = - ./о ^з (1024^ + (27520! + 179232)01 + (241232 + 3888ЗШ+
4(301 + 403
+ 76802)02 + (68403 + 25200202 + 11520102)0з + (35101 + 43202)0202) -
В таблице 2 сравниваются результаты для кратностей удлинения А2(п), А2(0), полученные путём численного интегрирования и асимптотического разложения. Заметим, что во втором приближении результаты оказываются достаточно точны даже для больших М.
Таблица 2. Значения А2(7г), А2(0) для различных ^ при <3 = 1.5, 91 = 92 = 93 = Ь
м Асимптотика Численные результаты
Первое приближение Второе приближение
А2(тг) А2(0) А2(тг) А2(0) А2(тг) А2(0)
0.350 0.962 2.848 1.079 2.964 1.114 2.826
0.300 1.279 3.165 1.379 3.265 1.394 3.180
0.250 1.724 3.609 1.807 3.693 1.811 3.641
0.200 2.391 4.276 2.457 4.343 2.456 4.313
0.150 3.502 5.387 3.552 5.437 3.549 5.422
0.100 5.724 7.609 5.757 7.643 5.756 7.637
0.050 12.391 14.276 12.407 14.293 12.407 14.291
В таблице 3 заметно значительное уточнение значения Q* за счёт следующего члена приближения.
Таблица 3. Значения Q * для трех случаев соотношении ^кёсткости
Численное решение Асимптотическое решение
ii = 0.1 ii = 0.01 Q 0 Q о + nQi
11 = 0.1 li = 0.01
91 = 92 = 1, 93 = 0 4.12 5.71 6.00 3.83 5.78
91 = 92 = 93 = 1/2 4.59 6.15 6.43 4.49 6.23
91 = 92 = 0, 93 = 1 4.71 5.78 6.00 4.80 5.88
Полученные следующие приближения значительно уточняют асимптотические выражения для основных переменных и предельного давления Q* и применимы для более широкого диапазона значений
Литература
1. Полякова Е. В., Чайкин В. А. Прикладные задачи механики мягких оболочек и тканей. СПб.: ИПЦ СПГУТД, 2006.
2. Полякова Е. В., Товстик П. Е., Филиппов С. Б., Чайкин В. А. Осесимметричная деформация торообразной оболочки из нитей под действием внутреннего давления // Вестн. С.-Петерб. ун-та. Сер. 1. 2009. Вып. 4. С. 89-104.
3. Полякова Е. В., Товстик П. Е., Филиппов С. Б. Осесимметричная деформация мягкой армированной нитями тороидальной оболочки // Вестн. С.-Петерб. ун-та. Сер. 1. 2011. Вып. 3. С. 131-142.
Статья поступила в редакцию 25 декабря 2014 г. Сведения об авторе
Кузьмин Андрей Владимирович —студент; ndrewk@ya.ru AXISYMMETRIC DEFORMATION OF
THE REINFORCED BY THREADS TOROIDAL MEMBRANE
A V. Kuzmin
St.Petersburg State University, Universitetskaya nab., 7-9, St.Petersburg, 199034, Russian Federation; ndrewk@ya.ru
Nonlinear axisymmetric deformation of a toroidal shell under action of internal pressure is considered. The shell is reinforced by two system of threads located on parallels and meridians. It is supposed, that threads are disposed frequently enough therefore after averaging we get the anisotropic elastic membrane. The basic feature of a membrane is that it can not hold the compression stresses. In state of equilibrium all meridians of the shell are stretched, but the part of parallels can be compressed. That leads to formation of folds.
For search of deformations and displacements of a shell the system of the ordinary differential equations of the fourth order is received. Geometrical and physical nonlinearity is taken into account. The method of numerical solution and also the method of asymptotical integration in the case when the meridian radius is much smaller than the parallel one are elaborated. Comparison of asymptotic and numerical results is performed.
It is shown, that there is an ultimate pressure at which excess states of equilibrium are absent. If pressure is less that the ultimate one there are two states of equilibrium — subcritical and supercritical. Asymptotic and numerical solutions for the minimal value of pressure at which folds on the shell are not formed are found. Refs 3. Figs 2. Tables 3.
Keywords: anisotropic membrane, geometrical and physical nonlinearity, internal pressure, asymptotic and numerical solutions.
References
1. Polyakova E. V., Chaikin V. A., Applied problems of a mechanics of soft shells and fabrics (SPSUTD, Saint-Petersburg, 2006).
2. Polyakova E. V., Tovstik P. E., Filippov S. B., Chaikin V.A., "Axisymmetric deformation of a toroidal membrane composed of threads under action of an internal pressure", Vestnik SPbGU Ser. 1 Issue 4, 89-104 (2009).
3. Polyakova E. V., Tovstik P. E., Filippov S. B., "Axisymmetric deformation of a toroidal membrane composed of threads", Vestnik SPbGU Ser. 1 Issue 3, 131-142 (2011).