close

Вход

Забыли?

вход по аккаунту

?

Низкочастотные длинноволновые приближения трёхмерных динамических уравнений для случая двухслойной вязкоупругой пластины.

код для вставкиСкачать
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2012. № 4 (29). С. 115–121
УДК 539.3
НИЗКОЧАСТОТНЫЕ ДЛИННОВОЛНОВЫЕ ПРИБЛИЖЕНИЯ
ТРËХМЕРНЫХ ДИНАМИЧЕСКИХ УРАВНЕНИЙ ДЛЯ СЛУЧАЯ
ДВУХСЛОЙНОЙ ВЯЗКОУПРУГОЙ ПЛАСТИНЫ
Н. С. Анофрикова, М. В. Вильде
Саратовский государственный университет им. Н. Г. Чернышевского,
410012, Россия, Саратов, ул. Астраханская, 83.
E-mails: nanofrikova@yandex.ru, mv_wilde@mail.ru
Представлены асимптотические методы, разработанные для вывода низкочастотных длинноволновых приближений трёхмерных динамических уравнений
для случая двухслойной пластины, выполненной из вязкоупругих материалов.
Путём асимптотического интегрирования точных трёхмерных уравнений выведены двумерные уравнения для асимптотически главных компонент напряжённо-деформированного состояния для тангенциального и поперечного приближений.
Ключевые слова: двухслойная пластина, вязкоупругость, асимптотические методы, низкочастотные длинноволновые приближения.
Введение. В работах [1–5] описаны асимптотические методы, разработанные для исследования динамического поведения тонкостенных конструкций
различной геометрии, выполненных из упругих и вязкоупругих материалов.
В [1] рассмотрены асимптотические приближения трёхмерных уравнений теории упругости для тонких пластин и оболочек. Описаны четыре типа приближений, выделенные в зависимости от значений показателей изменяемости и динамичности, а также представлено использование этих приближений
для решения динамических задач различных типов, в частности, нестационарных задач для оболочек вращения. Асимптотический подход к решению
задач подобного класса описан также в монографии [2]. В работе [3] систематизированы результаты разработки асимптотических методов исследования
нестационарных волн для тонких упругих оболочек вращения. В [4] были решены аналогичные задачи для случая вязкоупругой оболочки. Статья [5] посвящена изложению методики построения низкочастотных длинноволновых
приближений для случая двухслойной упругой пластины. В настоящей работе предлагается методика построения аналогичных приближений для случая
двухслойной пластины, выполненной из вязкоупругих материалов.
1. Постановка задачи. Рассмотрим бесконечную двухслойную пластину,
оба слоя которой выполнены из вязкоупругих материалов, свойства которых
описываются моделью стандартного вязкоупругого тела с условием упругого объемного расширения. Введем декартову систему координат (x1 , x2 , z),
совмещая плоскость Ox1 x2 со срединной плоскостью пластины и направляя
ось z по нормали к срединной плоскости. Введем обозначения: l — номер слоя
(l)
(l)
(l = 1, 2), σij — напряжения, ui — перемещения в l–том слое пластины; 2hl —
Наталия Сергеевна Анофрикова (к.ф.-м.н., доц.), доцент, каф. математической теории
упругости и биомеханики.
Мария Владимировна Вильде (д.ф.-м.н.), ведущий научный сотрудник, образовательнонаучный институт наноструктур и биосистем.
115
А н о ф р и к о в а Н. С., В и л ь д е М. В.
толщина l–того слоя, 2h — толщина пластины.
Будем предполагать, что наружные поверхности пластины свободны от
нагрузки. Тогда граничные условия на них имеют вид (k = 1, 2, 3):
(2)
при z = −h : σ3k = 0,
(1)
при
z = h : σ3k = 0.
(1)
Граничные условия на стыке двух слоев пластины — условия непрерывного
контакта — сформулируем следующим образом:
(1)
(2)
(1)
(2)
при z = z1 : σ3k = σ3k , uk = uk ,
(2)
где z1 = h − 2h1 .
Приведём точные трёхмерные динамические уравнения теории вязкоупругости для пластины. Уравнения движения возьмём в виде
(l)
(l)
(l)
(l)
∂σji
∂σii
∂ 2 ui
∂σ
+
+ 3i − ρl
= 0,
∂xi
∂xj
∂z
∂t2
(l)
(l)
(l)
(l)
∂σj3 ∂σ33
∂σi3
∂ 2 u3
− ρl
+
+
= 0, (3)
∂xi
∂xj
∂z
∂t2
где ρl — плотность материала слоя, t — время.
Исходя из уравнений, приведённых в [6], уравнения состояния для l-того
слоя можно записать следующим образом:
El
1
∂
+
t2l ∂t
(l)
∂ui
1 1 − 2νl
1 + νl
∂
(l)
=
+2
+
σii +
∂xi
3
t2l
t1l
∂t
1 1 − 2νl 1 + νl
∂
(l)
(l)
+
−
− νl
(σjj + σ33 ),
3
t2l
t1l
∂t
El
∂
1
+
t2l ∂t
(l)
1 1 − 2νl
∂u3
1 + νl
∂
(l)
=
+2
+
σ +
∂z
3
t2l
t1l
∂t 33
1 1 − 2νl
1 + νl
∂
(l)
(l)
+
−
− νl
(σii + σjj ), (4)
3
t2l
t1l
∂t
(l) !
(l)
∂uj
∂ui
El
∂
1
∂
1
(l)
+
+
+
=2
σij ,
1 + νl t2l ∂t
∂xj
∂xi
t1l ∂t
!
(l)
(l)
∂ui
∂
∂u3
1
∂
1
El
(l)
+
+
+
=2
σ3i (i 6= j = 1, 2; l = 1, 2),
1 + νl t2l ∂t
∂z
∂xi
t1l ∂t
где t1l — характерное время релаксации, t2l — характерное время ползучести,
El , νl — мгновенные значения модуля Юнга и коэффициент Пуассона материала l-того слоя соответственно.
Произведём в уравнениях (3), (4) растяжение масштабов независимых переменных по формулам
xi = Lη q ξi ,
116
z = Lηζ,
a
t = Lc−1
21 η τ,
(5)
Низкочастотные длинноволновые приближения трëхмерных динамических уравнений . . .
где q — показатель изменяемости, a — показатель динамичности, c21 — скорость
волны сдвига в первом слое, η = hL−1 ≪ 1 — относительная полутолщина
пластины, L — характерный размер длины.
Предположим, что дифференцирование по безразмерным переменным ξi ,
ζ, τ не меняет асимптотический порядок неизвестных величин. Кроме того, будем рассматривать случай, когда скорости волн сдвига для материалов
первого и второго слоев — величины одного порядка. Введение независимых
переменных (5) позволяет методом асимптотического интегрирования трёхмерных уравнений теории вязкоупругости (3), (4) вывести асимптотически
приближенные уравнения для составляющих напряжённо-деформированного
состояния (НДС) при различных показателях изменяемости и динамичности.
В настоящей работе остановимся на случае так называемых длинноволновых низкочастотных приближений. К этому виду относятся приближения,
для которых показатели динамичности и изменяемости удовлетворяют неравенствам q < 1, a < 1 [2]. Длинноволновые приближения разделяют на тангенциальные и поперечные, соответствующие теориям растяжения и изгиба
тонких пластин соответственно.
2. Низкочастотные длинноволновые тангенциальные приближения. Начнём рассмотрение с тангенциального приближения. В этом случае тангенциальные компоненты вектора перемещений велики по сравнению с его нор(l)
(l)
мальной компонентой: ui ≫ u3 (i = 1, 2).
Оценим величины времен ползучести и релаксации, вводя показатели их
интенсивности по формулам
ril
til = Lc−1
21 η τil ,
(6)
и будем предполагать, что ril 6 a (i = 1, 2).
При построении тангенциального приближения показатели изменяемости
и динамичности для каждого слоя связаны соотношением q = a. Введем
следующие асимптотики для компонент НДС [5]:
(l)
(l)
ui = Lη q ũi ,
(l)
(l)
(l)
u3 = Lηũ3 ,
(l)
σ3i = El η 1−q σ̃3i ,
(l)
(l)
(l)
σii = El σ̃ii ,
(l)
σ33 = El η 2−2q σ̃33
(l)
(l)
σij = El σ̃ij ,
(i, j = 1, 2).
(7)
Предполагаем, что величины с «тильдой» имеют один и тот же асимптотический порядок.
В силу выбора асимптотик (7), в уравнения движения, записанные с учётом (5), (6), в рамках погрешности O(η 2−2q ) входят слагаемые, содержащие
(l)
производные по ζ от σ̃3i (i = 1, 2, 3). Такой выбор асимптотик позволяет удовлетворить всем граничным условиям при асимптотическом интегрировании.
В результате асимптотического интегрирования уравнений (3), (4) устанавливаем следующую зависимость компонент НДС от нормальной координаты:
(l)
(l)
ũi = ui,0 ,
(l)
(l)
(l)
(l)
ũ3 = u3,0 + ζu3,1 ,
(l)
(l)
σ̃3i = σ3i,0 + ζσ3i,1 ,
(l)
(l)
(l)
σ̃ii = σii,0 ,
(l)
(l)
(l)
(l)
σ̃ij = σij,0 ,
(l)
σ̃33 = σ33,0 + ζσ33,1 + ζ 2 σ33,2 ,
(8)
где величины с цифрой после запятой в нижнем индексе от ζ не зависят.
117
А н о ф р и к о в а Н. С., В и л ь д е М. В.
Переход в граничных условиях (1), (2) к представлениям (7), (8) позволяет
установить связь между компонентами НДС первого и второго слоёв:
(1)
(2)
(1)
ui,0 = ui,0 ,
(2)
2E1 h1 σ3i,1 + 2E2 h2 σ3i,1 = 0.
В результате асимптотического интегрирования получены: система отно(l)
(l)
(l)
сительно асимптотически главных компонент НДС ui,0 , σii,0 , σij,0 и система,
определяющая асимптотически второстепенные компоненты через асимптотически главные.
Приведём вид размерной двумерной формы записи полученной системы
для асимптотически главных компонент НДС. Введём перемещения ui , усилия Ti , Sij и усреднённую плотность ρ по формулам
(1)
(2)
(1)
(2)
ui = Lη q ui = Lη q ui , Ti = 2 h1 σii + h2 σii ,
(9)
(1)
(2)
Sij = 2 h1 σij + h2 σij , ρ = (ρ1 h1 + ρ2 h2 )/h.
С учётом (9) система разрешающих уравнений для асимптотически главных компонент НДС примет вид
где
∂Ti ∂Sij
∂ 2 ui
+
− 2ρh 2 = 0,
∂xi
∂xj
∂t
∂uj
h2 E2
∂ui
h1 E1
f12 f21 +
f11 f22
+
= f11 f12 Sij ,
1 + ν1
1 + ν2
∂xj
∂xi
∂ui
2
2
2
2
− f41
− f42
+ h2 E2 f22 f32 f31
2 h1 E1 f21 f31 f32
−
∂xi
∂uj
2
2
2
2
=
− 2 h1 E1 f21 f41 f32
− f42
+ h2 E2 f22 f42 f31
− f41
∂xj
2
2
2
2
= f31
− f41
f32
− f42
Ti ,
(10)
∂
1 + νl
1 1 − 2νl
∂
1
+
+2
, f3l =
+ ,
fil =
til ∂t
3
t2l
t1l
∂t
1 1 − 2νl 1 + νl
∂
f4l =
−
− νl , i 6= j = 1, 2; l = 1, 2.
3
t2l
t1l
∂t
3. Низкочастотное длинноволновое поперечное приближение. Теперь рассмотрим случай поперечного приближения. В этом случае нормальная компонента вектора перемещений велика по сравнению с его тангенциальными
(l)
(l)
компонентами: u3 ≫ ui (i = 1, 2).
Величины времен ползучести и релаксации будем оценивать так же, как и
в предыдущем случае, вводя показатели их интенсивности по формулам (6).
Как и раньше, будем предполагать, что ril 6 a (i = 1, 2).
При построении поперечного приближения показатели изменяемости и
динамичности для каждого слоя связаны соотношением q = 2a − 1, q > 1/2.
Введём следующие асимптотики для компонент НДС:
(l)
(l)
ui = Lηũi ,
(l)
(l)
(l)
u3 = Lη q ũ3 ,
(l)
σ3i = El η 2−2q σ̃3i ,
118
(l)
(l)
σii = El η 1−q σ̃ii ,
(l)
(l)
σ33 = El η 3−3q σ̃33
(l)
(l)
σij = El η 1−q σ̃ij ,
(i, j = 1, 2).
(11)
Низкочастотные длинноволновые приближения трëхмерных динамических уравнений . . .
Предполагается, что величины с «тильдой» имеют один и тот же асимптотический порядок.
Как и в случае тангенциального приближения, в силу выбора асимптотик
(11) в уравнения движения, записанные с учётом (5), (6), в рамках погреш(l)
ности O(η 2−2q ) входят слагаемые, содержащие производные по ζ от σ̃3i (i =
= 1, 2, 3). Такой выбор асимптотик позволяет удовлетворить всем граничным
условиям при асимптотическом интегрировании.
Зависимость компонент НДС от нормальной координаты, установленная
в результате асимптотического интегрирования уравнений (3), (4), в данном
случае имеет вид
(l)
(l)
(l)
(l)
(l)
(l)
(l)
(l)
(l)
(l)
(l)
(l)
(l)
(l)
(l)
ũi = ui,0 + ζui,1 , ũ3 = u3,0 , σ̃ii = σii,0 + ζσii,1 , σ̃ij = σij,0 + ζσij,1 ,
(l)
(l)
σ̃3i = σ3i,0 + ζσ3i,1 + ζ 2 σ3i,2 ,
(l)
(l)
(l)
σ̃33 = σ33,0 + ζσ33,1 + ζ 2 σ33,2 + ζ 3 σ33,3 ,
(12)
где величины с цифрой после запятой в нижнем индексе от ζ не зависят.
Переходя в граничных условиях (1), (2) к представлениям (11), (12), получим соотношения
(1)
(2)
(1)
(2)
ui,0 = ui,0 , u3,0 = u3,0 .
Приведём вид двумерной разрешающей системы, полученной в результате
асимптотического интегрирования системы (3), (4):
∂Mi ∂Hij
∂N1 ∂N2
∂2w
∂Ti ∂Sij
+
= 0,
+
− Ni = 0,
+
− 2ρh 2 = 0,
∂xi
∂xj
∂xi
∂xj
∂x1
∂x2
∂t
h1 E1
∂ui
h2 E2
∂uj
f12 f21 +
f11 f22
+
−
1 + ν1
1 + ν2
∂xj
∂xi
2
E2
∂ w
E1
f21 f12 −
f11 f22
= f11 f12 Sij ,
−2h1 h2
1 + ν1
1 + ν2
∂xi ∂xj
∂ui
2
2
2
2
−
2 h1 E1 f21 f31 (f32
− f42
) + h2 E2 f22 f32 (f31
− f41
)
∂xi
∂uj
2
2
2
2
−
−2 h1 E1 f21 f41 (f32
− f42
) + h2 E2 f22 f42 (f31
− f41
)
∂xj
∂2w
(13)
2
2
2
2
− f42
) − E2 f22 f32 (f31
− f41
)
−2h1 h2 E1 f21 f31 (f32
+
∂x2i
∂2w
2
2
2
2
+2h1 h2 E1 f21 f41 (f32
− f42
) − E2 f22 f42 (f31
− f41
)
=
∂x2j
2
2
2
2
= (f31
− f41
)(f32
− f42
)Ti ,
∂ui
2
2
2
2
−
2h1 h2 E1 f21 f31 (f32
− f42
) − E2 f22 f32 (f31
− f41
)
∂xi
∂uj
2
2
2
2
−
− f42
) − E2 f22 f42 (f31
− f41
)
−2h1 h2 E1 f21 f41 (f32
∂xj
hh
1
2
2
−2
(h2 + 3h22 )E1 f21 f31 (f32
− f42
)+
3 1
119
А н о ф р и к о в а Н. С., В и л ь д е М. В.
i ∂2w
h2 2
2
2
+
(h2 + 3h21 )E2 f22 f32 (f31
− f41
)
3
∂x2i
hh
1
2
2
+2
(h2 + 3h22 )E1 f21 f41 (f32
− f42
)+
3 1
i ∂2w
h2
2
2
− f41
)
=
+ (h22 + 3h21 )E2 f22 f42 (f31
3
∂x2j
+
2
2
2
2
= (f31
− f41
)(f32
− f42
)Mi ,
E1
E2
∂uj
∂ui
h1 h2
f12 f21 −
f11 f22
+
−
1 + ν1
1 + ν2
∂xj
∂xi
2
2h1 2
E1
E2
2h2 2
∂ w
2
2
(h1 + 3h2 )
(h2 + 3h1 )
f21 f12 +
f11 f22
=
−
3
1 + ν1
3
1 + ν2
∂xi ∂xj
= f11 f12 Hij ,
где перемещения ui , w, усилия Ti , Sij , моменты Mi , Hij , перерезывающие
силы Ni определяются по формулам
(1)
(2)
ui = Lηui,0 = Lηui,0 ,
Z h
Z z1
(2)
(1)
σii dz +
σii dz,
Ti =
Mi =
Z
−h
z1
(2)
zσii dz
−h
+
Z
z1
h
z1
Ni =
(1)
zσii dz,
(1)
−h
z1
Hij =
Z
z1
(2)
zσij dz
−h
Z
z1
(2)
σ3i dz +
−h
(2)
w = Lη q u3,0 = Lη q u3,0 ,
Z h
Z z1
(2)
(1)
σij dz +
σij dz,
Sij =
Z
h
z1
+
Z
h
z1
(1)
zσij dz,
(1)
σ3i dz,
а остальные величины имеют тот же смысл, что и раньше.
Из двумерных уравнений динамической теории вязкоупругости для скорости продольной волны по двумерной теории получаем следующее выражение:
2
X
El hl
1
2
c =
.
ρ1 h1 + ρ2 h2
1 − νl2
l=1
Переходя в уравнениях (10), (13) к пределу при t1l → ∞ и t2l → ∞,
можно получить двумерные уравнения тангенциального и поперечного приближений для упругой двухслойной пластины. Если в последних положить
E1 = E2 и ν1 = ν2 , то уравнения тангенциального приближения переходят в
уравнения обобщенного плоского напряженного состояния, а уравнения поперечного приближения — в уравнения теории изгиба пластин Кирхгофа.
Работа выполнена при поддержке РФФИ (проект № 11−01−00545).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Коссович Л. Ю. Нестационарные задачи теории упругих тонких оболочек. Саратов: Саратовск. ун-т, 1986. 176 с. [Kossovich L. Yu. Transient problems of theory of thin elastic
shells. Saratov: Saratov University, 1986. 176 pp.]
120
Низкочастотные длинноволновые приближения трëхмерных динамических уравнений . . .
2. Kaplunov Ju. D., Kossovich L. Yu., Nolde E. V. Dynamics of thin walled elastic bodies. SanDiego: Academic Press, 1998. x+226 pp.
3. Коссович Л. Ю. Асимптотические методы в динамике оболочек при ударных воздействиях // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика,
2008. Т. 8, № 2. С. 12–33. [Kossovich L. Yu. Asymptotic methods in dynamics of shells under
shock loading // Izv. Saratov. Univ. Mat. Mekh. Inform., 2008. Vol. 8, no. 2. Pp. 12–33].
4. Бажанова Н. С., Коссович Л. Ю., Сухоловская М. С. Нестационарные волны в вязкоупругих оболочках: модель Максвелла // Изв. высш. учеб. завед. Сев.-Кавк. Регион.
Естеств. науки, 2000. № 2. С. 17–24. [Bazhanova N. S., Kosovich L. Yu., Sukholovskay M. S.
Non-stationary in waves viscoelastic shells: Maxwell’s model // Izv. Vyssh. Uchebn. Zaved.,
Sev.-Kavk. Reg., Estestv. Nauki, 2000. no. 2. Pp. 17–24].
5. Коссович Л. Ю., Шевцова Ю. В. Асимптотические приближения трёхмерных динамических уравнений теории упругости в случае двухслойных пластин // Проблемы прочности и пластичности, 2005. № 67. С. 102–110. [Kossovich L. Yu.. Shevtsova Yu. V.
Asymptotic approaches of three-dimentional dynamic equations of elasticity theory in the
case of two-layered plates // Problemy prochnosti i plastichnosti, 2005. no. 67. Pp. 102–110].
6. Новацкий В. Динамика сооружений. М.: Госстройиздат, 1963. 376 с. [Nowacki W.
Dynamics of Constructions. Moscow: Gosstroyizdat, 1963. 376 pp.]
Поступила в редакцию 19/V/2012;
в окончательном варианте — 14/IX/2012.
MSC: 74J30; 74K25, 74K20, 74H10
LOW-FREQUENCY LONG-WAVE APPROXIMATIONS OF THE
3D DYNAMIC EQUATIONS FOR THE CASE OF DOUBLE-LAYERED
VISCOELASTIC PLATE
N. S. Anofrikova, M. V. Wilde
Saratov State University named after N. G. Chernyshevsky,
83, Astrakhanskaya st., Saratov, 410012, Russia.
E-mails: nanofrikova@yandex.ru, mv_wilde@mail.ru
The asymptotic methods developed to obtain low-frequency long-wave approximations of
the 3D dynamic equations for the case of double-layered viscoelastic plate are described.
The 2D equations for the leading tangential and transverse approximations of stressstrain state are derived. The method of asymptotic integration of exact 3D equations
is applied.
Key words: double-layered plate, viscoelasticity, asymptotic methods, low-frequency
long-wave approximations.
Original article submitted 19/V/2012;
revision submitted 14/IX/2012.
Nataliya S. Anofrikova (Ph. D. (Phys. & Math.)), Associate Professor, Dept. of Mathematical
Theory of Elasticity and Biomechanics.
Mariya V. Wilde (Dr. Sci. (Phys. & Math.)), Leading Research Scientist, Educational Research
Institute of Nanostructures and Biosystems.
121
1/--страниц
Пожаловаться на содержимое документа