Список использованных источников
1. Лыков, А.В. Теория сушки / А.В. Лыков. - Москва: Энергия, 1968, - 469 с.
2. Красников, В. В. Кондуктивная сушка / В. В. Красников. - Москва: Энергия, 1973, - 365 с.
3. Филоненко, Г.К. Сушильные установки / Г.К. Филоненко, П.Д. Лебедев. -Москва: Госэнергоиздат, 1952, - 256 с.
4. Ольшанский, А. И. Некоторые закономерности кинетики сушки влажных материалов // А. И. Ольшанский, В. И. Ольшанский, Е. Ф. Макаренко, -Инженерно-физический журнал / Минск - июль-август 2007. - т. 80, - № 4, -С. 143-146.
5. Куц, П.С. Некоторые закономерности тепловлагообмена и приближенные методы расчета кинетики процесса сушки влажных материалов // П.С. Куц, А. И. Ольшанский, - Инженерно-физический журнал / Минск -1977. - т. 32, - № 6, - С. 1007-1014.
SUMMARY
Methods of calculation process of dryings of damp materials for the period of the falling speed, based on the generalized characteristics process of dryings - relative speed of the drying, generalized and relative time of drying are considered.
УДК 539.3
СВОБОДНЫЕ КОЛЕБАНИЯ МНОГОСЛОЙНОЙ НЕКРУГОВОЙ ЦИЛИНДРИЧЕСКОЙ ОБОЛОЧКИ, СОСТОЯЩЕЙ ИЗ ВЯЗКОУПРУГИХ СЛОЕВ
Е.А. Корчевская
Постановка задачи. Рассматривается тонкая некруговая цилиндрическая оболочка постоянной длины ^ состоящая из N изотропных вязкоупругих слоев,
Ь Е
характеризующихся толщиной к, модулем Юнга к, плотностью Рк и
коэффициентом Пуассона у к .
В качестве исходных могут быть использованы уравнения [1], основанные на гипотезах, сформулированных Э.И. Григолюком и Г.М. Куликовым [2], которые отличаются от классических уравнений полубезмоментной теории тонких оболочек наличием дополнительных слагаемых, учитывающих поперечные сдвиги слоев. Отбрасывание последних приводит к уравнениям для изотропной оболочки с физическими характеристиками, равными осредненным по толщине параметрам слоистой исходной оболочки и может давать существенные погрешности при расчетах:
Eh
Пз
1
ви
b
■А*
Л* А * *
А X +
1
д2 F
R (а 2) да 2
рЮ 2W* = 0,
2
А* А* F -
Eh д 2W R2 (а 2) да 2
= 0,
W*
i - hi а*
b
*
X •
(1)
А* а а
„„—- - оператор Лапласа в криволинейной системе координат 1, 2, E,
v, р -осредненные модуль Юнга, коэффициент Пуассона и плотность материала
* *
соответственно, h - толщина оболочки, F , х - функции напряжений и
*
перемещений, W - нормальный прогиб, ^ - частота собственных колебаний, /о(аз), fk (аз), S (аз) - функции, зависящие от поперечной координаты а,
Пз в Ъ , Yk , tfb n1k , n2k , n3k, Çk, zn , C13, C12, Xkn, Xk , Gk ,
q 44 Gk ck, ck
параметры
определяются по формулам [1, 2]:
N
h = S hk V =
k=1
Yk =
EA~k
1 - V 2
N E h ~
S 7 2
Vk =1 1 - Vk y
S Ekhk~k
^ 1 2
V k=11 - Vk y
4-1
,-1
s VkEkhkck
k=1 1 - V
2
E =
1 - V
2 ^ n E,h,c,
h
s 1
k=1 1 - V
k
n
p = Spkçk 0 =1--Ъ =
k=1 k k П1П3
=12(1 - v 2 q 44
Ehn1
N -1 2 N
П1 = s Ç- n1k Y k - 5c12 5 П2 = s ^ П 2k Y k - 3c13 C12 5
k=1
k=1
1 "k
П3 = 4Sfc + 3fk-1Zk)k - 3c?3 12^k = "S S2(а)а3.
k=1 "k-1
1 3 "k 1 2 "k 12 hn2k = ia3 S(а3 ¥а3 2 h n3k = i S(а3
12 "k-1 2 "k-1
"k-1 k),
hÇk = hk, = "„ (n=
(n=0,
C13 = s ((-1 + z- )y k, C12 = s 3k Y k ■
k=1 k=1
"k
Xkn = i fk (а3 )fn (а3 3, (n = 0k)
"k-1
k
Xk = i f02 (а3 У«3, f 0(а 3) = -12 (а 3-S 0 )S
h
N а 3 )
k-1
£(а 3) = i fо(а 3)^а 3
fk(а 3) = 7Г(а 3 - S k-1XS k - а 3)
h2
g- = EJ [2(1 + v k )] g- = g-~-
q44 =
N
s
k=1
С
X
2 Л
1 'k 0
Xt -
X
kk
N
s
k=1
X
2
1 'k 0
Xt - ■
X
kk
G
-1
N X2 ~
, V Xk 0/~
+ S—^k
k=1 Xkk
а
2
+х>
~ = 1 - ^, ^ =| Як () sds
0 . (2) В формулах (2) Ик (в) - ядро релаксации напряжений материала для к-ого слоя,
а О = Ю + /а, где 1 = V-!, щ - фундаментальная частота свободных колебаний, а - декремент колебаний.
Перепишем уравнения (1) в безразмерном виде:
д 2 ^
в4 (l - в3та)2х + k(ф)^-г - х( - Е2кАх)= 0,
дя 2
в 4 А2 F - k (ф) ^ (х-в 2 кАх)= 0,
(3)
дя 2
в8 = h2n 3/[l2R2 (l -v2)] -
где 13' l v /j - малый параметр, характеризующий
тонкостенность оболочки, А - оператор Лапласа в криволинейной системе
координат Ф = а 2/R и * = а 1/R, F = ^/И^2), X = X* /R-безразмерные функции напряжений и перемещений соответственно,
X = (pR У(Ев - искомый частотный параметр, k(ф- переменная кривизна. Здесь т ,к - параметры, характеризующие поперечные сдвиги:
к/п2 =в2к, ке/п2 =в3т, к,тпри г^0, где K = П2hV(bR2
Граничные условия имеют вид:
F = aF = х = ах = а х = 0, при s-ОЛ 1 = LR. (4)
Задача состоит в определении параметра X, для которого краевая задача (3), (4) имеет ненулевое решение.
Построение решения. Считаем, что локализация собственных колебаний
происходит в окрестности некоторой "слабой" образующей ф = ф0. Введем растяжение масштаба в окрестности этой образующей:
Мф-ф0 )в"1/2.
Согласно [3] решение будем искать в виде:
Х = Хm (ф)Sin(mns//), F = ф т (ф)sin(mnVl), m = 1,2,..., (5)
где
<х> I f 1 л
{хm,Фm }= !в j2|хmi( f j ft)}expJ i в--2q^ + -a^2
| Л/ mi VS/э J mi j=0 1
2
X = X 0 + вХ1 + в2 X 2 +...,
Im a > 0, (6) (7)
k (ф) = k (ф0 ) + в^2 k '(ф0 ) + - вk "(ф0 H2
+ ...
2 (8) В выражении (6) параметр д характеризует изменяемость решения в окружном направлении, а мнимая часть числа а, характеризующего скорость затухания
амплитуды волн при удалении от линии ф ф0, должна быть положительной.
x ■ ( ■
Функции т] , т] являются полиномами по 4.
Подставляя (5)-(8) в (3), (4), получим последовательность уравнений:
í А к X ■ -к = 0 ■ = к=0 , ■ 0,1,2, ...
(9)
относительно вектор-функции имеет вид:
X
(х т/, ^т!)
А„
Матрица 0
А =
4 ( 2
q -Л011 + щ к (ф0 )(ят// )2 (1 + щ 2)
к (Фо Х™*/1)2
q
а элементы матрицы
А ■
■ при ■ ~ 1 выражаются через производные по д и Фо
/ А п
- го порядка элементов матрицы 0 [3].
Из условия существования нетривиального решения системы (9) при находим формулу для частотного параметра нулевого приближения:
Л (Фо,т )= q V(1 + к2)+(к 2 (Фо Хт)4 У q4
Минимизируя
л 0 (q, ф 0,т)
по д и ф0, из условий
дЛ0/дq = дЛ0/дф0 =0
■ = 0
(10)
(11)
находим число д0 и "наиболее слабую" образующую ф 0. Однородная задача в нулевом приближении имеет решение в виде:
X 0 (4) = Ро № 0
где
Р<> (4) -■=1
неизвестный полином, а
У0 =(1, - Ли/Л12)
При 7 система уравнений (9) является неоднородной. Но при условиях (10), (11) она обращается в систему тождеств.
Условие совместности системы (9) при вычисления параметра а:
■ = 2
приводит к соотношению для
а=(ли
1/2
а также к уравнению относительно
й2 Р
ё42
+ га
Р0 (4):
йР
244-4+ро
а4
2Л о
+ Р +
2щ6
Ро = 0, (12)
Л0 ло
где фф и - вторые производные частотного параметра нулевого
0
приближения по соответствующим параметрам при
ф = ф0 и q = q0
т
^ = $) = (1/2 + я)л/О0~ + тд6/(1 + кд2) При 1 1 у/ ^ дд фф ^ ^ ' уравнение (12) имеет
решение в виде полинома Эрмита степени п:
р Щ = Нп (0) § = с = -ш
Пример. В качестве примера рассмотрим тонкую некруговую трехслойную цилиндрическую оболочку эллиптического сечения с полуосями Ь=0,03 м, а=0,015 м (а<Ь), постоянной длины L=0,45 м. В данном случае "наиболее слабыми" будут две образующие, проходящие через точки, где радиус кривизны поперечного сечения наименьший. Первый и третий слои имеют одинаковую толщину Л1=Л3 и
1 52 • 1012
изготовлены из керамики с модулем Юнга Е.]=Е3= ' Па, плотностью
^ 1 /"Ч А Л
Р1 = рз = кг/м3 и числом Пуассона 3 = ' . Межслойный заполнитель
изготовлен из фторопласта с Е2= 2,34 •10§ Па, р2 = 2150 кг/м3, у2 = 0,3 .
Будем считать здесь, что керамические слои подчиняются закону упругих деформаций, а межслойный заполнитель - закону вязкоупругих деформаций.
П = П — 0
Тогда 1 3 _ , а ядро скорости релаксации напряжений для фторопласта [4]
Я2 = 0,02366 е -3>33 •10-4 ''-°>95
Выполняя преобразование Лапласа функции И2, находим:
с = 0,02366 Г(0,05)
(¿ю - а + 3,33 • 10-4 )0'05
где Г(х) - Гамма- функция.
На рисунках 1,2 представлены графики параметров ш, а, как функций
- И2/ И И, = И, = 0,0001
относительной толщины заполнителя 2/ при 1 3 ' м и различных
значениях т. Цифрами 1, 2, 3 отмечены кривые, отвечающие волновым числам
т=1, 2, 7 соответственно.
Рисунок 1 - Собственная частота колебаний
Рисунок 2 - Декремент колебаний
к ¡к
В рассмотренном диапазоне изменения 0,2< 2 < 0,73 собственная частота ш возрастает с ростом Л2. Как и ожидалось, увеличение толщины заполнителя также приводит к росту декремента колебаний а.
Результаты работы могут быть использованы при проектировании тонкостенных элементов машин, а также габаритных тонкостенных инженерных сооружений в промышленном строительстве, летательных аппаратов и подводных тонкостенных объектов.
Список использованных источников
1. Ботогова, М.Г. Свободные колебания слоистых вязкоупругих цилиндрических оболочек / М.Г. Ботогова, Г.И. Михасев, Е.А. Корчевская // Вестн. Полоц. гос. ун-та. Сер. С, Фунд. науки. Механика. - 2006. - № 10. - С. 125-133.
2. Григолюк, Э.И. Многослойные армированные оболочки: расчет пневматических шин / Э.И. Григолюк, Г. М. Куликов. - М.: Машиностроение, 1988. - 288 с.
3. Товстик, П. Е. Устойчивость тонких оболочек: асимптотические методы / П.Е. Товстик. - М.: Наука. Физматлит, 1995. - 320 с.
4. Старовойтов, Э.И. Локальные и импульсные нагружения трехслойных элементов конструкций / Э.И. Старовойтов, А.В. Яровая, Д.В. Леоненко. -Гомель: БелГУТ, 2003. - 367 с.
SUMMARY
Using the asymptotic complex WKB- method free vibrations of the laminated composite non-circular cylindrical shell with viscoelastic filler is investigated. The natural vibration frequencies were found.
УДК 620.179
ФОРМИРОВАНИЕ АКУСТИЧЕСКОГО ПОЛЯ ПРЕОБРАЗОВАТЕЛЯ ПОДПОВЕРХНОСТНЫХ ВОЛН В ОБЪЕКТАХ С ГАЛТЕЛЬНЫМИ ПЕРЕХОДАМИ
А.Р. Баев, М.В. Асадчая, К.А. Филиппов
Для значительного числа объектов тепловой энергетики, химического машиностроения и других отраслей промышленности характерно наличие таких геометрических элементов, как различные выступы, радиусные переходы