ОПТИМИЗАЦИЯ ДЕКРЕМЕНТА ЗАТУХАНИЯ СВОБОДНЫХ КОЛЕБАНИЙ ВЯЗКОУПРУГОЙ СЛОИСТОЙ СФЕРЫ ПРИ ОГРАНИЧЕНИИ НА МАССУ

advertisement
ПРИКЛАДНАЯ МЕХАНИКА И ТЕХНИЧЕСКАЯ ФИЗИКА. 2000. Т. 41, N-◦ 2
161
УДК 539.3
ОПТИМИЗАЦИЯ ДЕКРЕМЕНТА ЗАТУХАНИЯ СВОБОДНЫХ
КОЛЕБАНИЙ ВЯЗКОУПРУГОЙ СЛОИСТОЙ СФЕРЫ
ПРИ ОГРАНИЧЕНИИ НА МАССУ
В. А. Будугаева
Институт физико-технических проблем Севера СО РАН, 677891 Якутск
Рассмотрена задача синтеза из конечного набора вязкоупругих материалов многослойной сферической оболочки, максимально гасящей собственные колебания, при ограничении на массу. Получены необходимые условия оптимальности, построен вычислительный алгоритм и приведен пример расчета.
1. В работе [1], посвященной решению задач оптимального проектирования слоистых
конструкций, на примере многослойной сферической оболочки, выполненной из заданного
набора вязкоупругих материалов, показано, что ее уровень демпфирования может быть повышен за счет изменения взаимного расположения и относительной толщины слоев. В [1]
с использованием принципа соответствия разработан алгоритм оптимального проектирования оболочек постоянной толщины, обеспечивающих максимальное гашение собственных колебаний. Представляется целесообразным обобщить полученный результат, вводя
некоторые ограничения на рабочие характеристики конструкции, например, ограничив ее
массу.
Сформулируем следующую задачу (см. [1]): из заданного набора материалов требуется сконструировать ограниченную по массе многослойную оболочку, обеспечивающую
минимальное значение выбранного критерия качества
F0 = Im [ω(θ)],
(1.1)
где ω = ωR + iωI — комплексная собственная частота; θ — распределение плотности ρ и
параметров Ламе λ̄ и µ̄ по радиальной координате.
Собственные частота и форма свободных колебаний сферической оболочки определяются из решения следующей задачи:
σr − σϕ
∂ 2u
∂u
u
∂σr
+2
=ρ 2,
σr = (λ̄ + 2µ̄)
+ 2λ̄ ,
∂r
r
∂t
∂r
r
(1.2)
u
∂u
σϕ = 2(λ̄ + µ̄) + λ̄ ,
l < r < R,
σr (l) = σr (R) = 0.
r
∂r
Здесь R — внешний радиус оболочки; l — внутренний радиус оболочки, определяемый в
ходе решения задачи; комплексные параметры Ламе имеют вид [2]
λ̄n = λn [1 − Γcλn (ωR ) − Γsλn (ωR )],
µ̄n = µn [1 − Γcµn (ωR ) − Γsµn (ωR )],
где
Γcλn (ωR )
Z∞
=
Rλn (τ ) cos (ωR τ ) dτ ;
0
Γsλn (ωR )
Z∞
=
Rλn (τ ) sin (ωR τ ) dτ ;
0
(1.3)
ПРИКЛАДНАЯ МЕХАНИКА И ТЕХНИЧЕСКАЯ ФИЗИКА. 2000. Т. 41, N-◦ 2
162
Γcµn (ωR )
Z∞
=
Γsµn (ωR )
Rµn (τ ) cos (ωR τ ) dτ ;
0
Z∞
=
Rµn (τ ) sin (ωR τ ) dτ ;
0
индекс n — номер слоя (n = 1, N ).
При вычислении интегралов (1.3) использовались выражения ядер релаксации вида [3]
R = ā exp (−βt)/t1−α , где ā, α, β — эмпирические константы.
Для однородной оболочки комплексная частота находится из трансцендентного уравнения, получаемого при решении задачи (1.2):
n 4µk 4µ
4µk 4µ
o
2
2
cos [k(R − l)]
−
ρω
−
−
ρω
+
l R2
R
l
n (4µk)2 4µ
4µ
o
2
2
+ sin [k(R − l)]
+ 2 − ρω
− ρω
= 0, (1.4)
Rl
l
R2
где k 2 = ρω 2 /(λ + 2µ) (черта над λ и µ опущена).
Для слоистой оболочки плотность и параметры Ламе будут кусочно-постоянными
функциями радиуса, а на границах слоев смещения и напряжения будут непрерывны. Это
позволяет обобщить решение (1.4) на случай слоистой оболочки и получить для вычисления напряжений и перемещений рекуррентное соотношение
!
!
uN (R)
u1 (l)
=G
.
(1.5)
σN (R)
σ1 (l)
При этом будут выполняться граничные условия
σN (R) = σ1 (l) = 0.
(1.6)
В соотношении (1.5) G — результирующая матрица вида
g11 g12 G=
= Λn · Λn−1 · · · Λ1 ,
g21 g22 где
an (rn ) ibn (rn )
Λn = ∆n cn (rn ) idn (rn )
−dn (rn−1 ) bn (rn−1 )
−icn (rn−1 ) ian (rn−1 )
∆n = (cn (rn−1 )bn (rn−1 ) − dn (rn−1 )an (rn−1 ))−1
;
(n = 1, N );
kn2
kn
kn
k2
cos (kn r) + 2 sin (kn r);
bn = − 2 cos (kn r) − n sin (kn r);
r
r
r
r
h
i
kn 4µn kn
4µn
cn =
cos (kn r) −
− ρn ω 2 sin (kn r) ;
r
r
r2
4µ
i
kn h 4µn kn
n
2
dn =
sin (kn r) +
−
ρ
ω
cos
(k
r)
, rn−1 < r < rn (n = 1, N ).
n
n
r
r
r2
Решая (1.5) с учетом условий (1.6), получим характеристическое уравнение для определения частоты многослойной сферической оболочки
an = −
g21 = 0.
(1.7)
2. Рассмотрим задачу оптимизации многослойной вязкоупругой полой сферы, обеспечивающей максимальное демпфирование свободных колебаний и имеющей ограничение
на массу. Используя представление u(r, t) = exp (iωt)z1 (r), σ(r, t) = exp (iωt)z2 (r), сведем
163
В. А. Будугаева
исходную задачу (1.2) к решению системы обыкновенных дифференциальных уравнений.
Вводя безразмерные переменные z1 = u/r0 , z2 = σr , r∗ = r/r0 , l∗ = l/r0 , ρ∗ = ρ/ρ0 ,
λ∗ = λ/σ0 , µ∗ = µ/σ0 , ω ∗2 = ω 2 r02 ρ0 /σ0 , где r0 , σ0 , ρ0 — характерные длина, напряжение и
плотность, и переводя переменный отрезок интегрирования [l, R] в постоянный [0, 1] путем
замены
r = l + x(R − l),
x ∈ [0, 1],
(2.1)
получим
Z 1 = AZ ≡ f ,
z2 (0) = z2 (1) = 0.
(2.2)
Здесь
2λ(R − l)
R−l
−
r(λ + 2µ)
λ
+
2µ
;
A=
Z = {z1 , z2 }; f = {f1 , f2 }.
i
h 4µ(3λ + 2µ)
4µ(R − l) 2
r2 (λ + 2µ) − ρω (R − l) − r(λ + 2µ) Будем считать, что внешний радиус R фиксирован, а внутренний l определяется при решении задачи. В качестве управления выберем пару {θ(x), l}, где θ(x) =
{ρ(x), λ(x), µ(x)} — распределение плотности и параметров Ламе. Функционал (1.1), имеющий смысл декремента затухания собственных колебаний сферы, минимизируется при
ограничении на массу сферы
F1 (ρ, l) ≡ P∗ − ηP 6 0.
Здесь P∗ — заданное значение массы; η — некоторое число.
Величину ω 2 (θ, l) можно найти из системы (2.2):
Z1
Z1
2
ω (θ, l) = J1 (x, Z, θ, l) dx
J2 (x, Z, θ, l) dx,
0
(2.3)
(2.4)
0
где
1
2 2
z
r , J2 (x, Z, θ, l) = ρz12 r2 ,
r2 (λ + 2µ)
λ + 2µ 2
а r зависит от x согласно (2.1). При этом главная часть приращения функционала (2.4)
на вариациях
(
t(x) = {ρ(x), λ(x), µ(x)}, x ∈ D, t ∈ W,
θ̃(x) =
(2.5)
θ(x) = {ρ̃(x), λ̃(x), µ̃(x)}, x ∈
/D
J1 (x, Z, θ, l) =
4µ(3λ + 2µ)
z12 +
(D ∈ [0, 1] — множество малой меры; W = {θ1 , . . . , θm } — заданное конечное дискретное
множество; m — число различных материалов) имеет вид
Z1
−1 Z
2
δω =
J2 ( · , θ) dx
[H( · , θ̃) − H( · , θ)] dx + Bδl .
(2.6)
0
D
i
2
2
Здесь H( · , θ̃) = J1
r z2 f1 ( · , θ̃) , аналогичный вид имеет H( · , θ);
2 ( · , θ̃) −
R−l
Z1 h
∂J1 ( · , θ)
∂J2 ( · , θ)
2
∂f1 ( · , θ) i
2
B=
−ω
−
r z2
dx; точкой обозначены пропущенные
∂l
∂l
R−l
∂l
h
0
аргументы x, Z, l.
( · , θ̃) − ω 2 J
ПРИКЛАДНАЯ МЕХАНИКА И ТЕХНИЧЕСКАЯ ФИЗИКА. 2000. Т. 41, N-◦ 2
164
Из ограничения (2.3) можно выразить δl, считая, что на вариациях (2.5) δF1 = 0:
Z
δP = (R − l)(ρ̃(x) − ρ(x))r2 dx + B1 δl = 0,
(2.7)
D
Z1
ρ(x)r[(R − l)(2 − 3x) − l] dx. Подставив (2.7) в (2.6), окончательно получим
где B1 =
0
Z1
Z
δF0 = δ Im (ω) =
J2 ( · , θ) dx
Im
D
−1
−1
(2ω)
[T ( · , θ̃) − T ( · , θ)] dx,
0
где T ( · , θ̃) = H( · , θ̃) − (B/B1 )(R − l)r2 ρ̃(x).
Построение минимизирующей последовательности управления [θ(x), l] основывается
на негативной формулировке принципа максимума Понтрягина (см. [4]), из которой следует существование в окрестности неоптимального управления [θ(x), l] нового управления
[θ(x), l + δl], которое улучшает конструкцию, удовлетворяя ограничению (2.3). Это управление ищется из условия минимальности δF0 , т. е.
Z1
−1
Z1
−1
Z
Z
−1
−1
Im
J2 (·, θ) dx
(2ω) T (·, θ̃) dx = min Im
J2 (·, θ) dx
(2ω) T (·, θ) dx. (2.8)
θ∈W
Di
0
0
Di
3. С учетом сказанного выше вычислительный алгоритм строится следующим образом:
1. Интервал [0, 1] разбиваем равномерной сеткой узлов {xn } на достаточно большое
число отрезков малой длины, моделирующих множество малой меры Di .
2. При допустимом управлении [θ(x), l] находим собственную частоту, решая уравнение (1.7) методом Мюллера [5].
3. Решаем систему (2.2) для найденной частоты, полагая, что на отрезке Di значение
вектора фазовых переменных характеризуется его значением в середине отрезка x = xi +
h/2 (h — длина малого отрезка).
4. На отрезке Di задаем новое управление θ̃ из условия (2.8).
5. Если новое управление совпадает со старым, то переходим к п. 2, увеличивая индекс
малого отрезка на единицу. В противном случае из условия (2.3) вычисляем δli , допуская,
что F1 + δF1 ≈ 0:
Z1
Z
2
2
δli = P∗ − η (R − l)[l + x(R − l)] ρ(x) dx − η (R − l)r [ρ̃(x)−ρ(x)] dx (ηB1 )−1 . (3.1)
0
D
При этом следует учесть, что мера отрезков Di и |δl| должна быть достаточно малой,
чтобы обеспечить использование линейного приближения (3.1).
6. С новым управлением [θ̃(x), li + δli ] переходим к п. 1 и рассматриваем отрезок Di+1 .
В качестве примера решим следующую задачу. Пусть задан набор из пяти вязкоупругих материалов, безразмерные характеристики которых приведены в таблице.
Требуется спроектировать сферическую оболочку, максимально гасящую собственные
колебания и имеющую массу в 3,5 раза меньше, чем однородная оболочка из самого плотного материала. Начальное приближение для управления выбиралось в виде [θ(x), l] = [1; 0,8],
что соответствовало однородной сфере с внутренним радиусом l = 0,8, выполненной из материала 1 (см. таблицу). Такая сфера имеет массу 0,65 и собственную частоту 6,687−3,29i.
165
В. А. Будугаева
Материал
ρ
E
ν
α
ā
β
1
2
3
4
5
4,00
2,86
1,75
1,00
2,90
90
50
30
15
65
0,25
0,25
0,25
0,25
0,25
0,4
0,2
0,2
0,2
0,3
1
1
1
1
1
0,05
0,05
0,05
0,05
0,05
В результате оптимизации получается четырехслойная сфера, состоящая из чередующихся материалов 2 и 3 (см. таблицу), с внутренним радиусом l = 0,925 87 и массой 0,1857.
Последовательность слоев: [l, R1 ] — материал 2, [R1 , R2 ] — материал 3, [R2 , R3 ] — материал 2, [R3 , 1] — материал 3, где R1 = 0,9718; R2 = 0,9778, R3 = 0,997 03 (см. рисунок).
Данная сфера имеет частоту 2,64−7,46 i, т. е. ее демпфирующие свойства улучшились
более чем в 2 раза.
ЛИТЕРАТУРА
1. Бондарев Э. А., Будугаева В. А., Гусев Е. Л. Синтез слоистых оболочек из конечного
набора вязкоупругих материалов // Изв. РАН. Механика твердого тела. 1998. N-◦ 3. С. 5–11.
2. Ильюшин А. А., Победря Б. Е. Основы математической теории термовязкоупругости. М.:
Наука, 1970.
3. Колтунов М. А. Ползучесть и релаксация. М.: Высш. шк., 1976.
4. Федоренко Р. П. Приближенное решение задач оптимального управления. М.: Наука, 1978.
5. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. М.:
Наука, 1968.
Поступила в редакцию 5/IV 1999 г.
Download