Список литературы:
Тененбаум М.М. «Сопротивление абразивному изнашиванию». Москва, «Машиностроение», 1976 г. 267 с.
Металлы. Метод испытания на абразивное изнашивание при трении о закрепленные абразивные частицы. ГОСТ 17367-
71. Госкомстандарт Совминст., Москва, 1971 г.
УДК 001 © Сафаров И.И., Болтаев М.Б. 2007 г.
АНАЛИТИЧЕСКОЕ РЕШЕНИЕ ЗАДАЧИ ВОЗДЕЙСТВИЯ ПОРОДЫ ПРИ ВЗРЫВНЫХ НАГРУЗКАХ НА КОНСТРУКЦИИ ТРУБОПРОВОДОВ И КРЕПЛЕНИЯ ТОННЕЛЕЙ
Сафаров И.И., зав. кафедрой «Высшая математика» Бух ТИП и ЛП, докт. физ.-мат. наук, профессор; Болтаев М.Б., асси- стент кафедры «Естественные и общетехнические дисциплины» Зарафшанского ОТФ НГГИ
Исследование распространения взрывных на- грузок (или волн) в кусочно–однородных упругих и вязкоупругих средах представляет значительный
научный и практический интерес. Это связано как с необходимостью изучения возможных повреж- дений элементов подземных конструкций при
прохождении через них импульсов напряжений, так и с различными технологическими приложе- ниями. В частности, волны напряжений использу- ются при измерении упругих постоянных материа- лов, обнаружении трещин и передачи информа- ции.
Необходимо отметить, что если напряженно- деформированное состояние элементов конструк-
Построив решение задачи для ступенчатой волны (1), с помощью интеграла Дюамеля легко построить решение для нестационарной волны с произвольной зависимостью от времени. Отсчет времени производится от момента, когда фронт набегающей волны (1) приводит к сечению обо- лочки z 0 .
В отсутствие статических массовых сил вектор
ций при статическом и квазистатическом нагруже- → T
нии может быть рассчитано при помощи хорошо
разработанных методов, то анализ распростране-
системы u ur , u , uz в упругой среде определя-
ется уравнением: → →
ния нестационарных волн в конструкциях пред- ставляет собой сложную математическую задачу.
Методом интегральных преобразований по-
2 graddiru rotrotu
u /
s 2 → t 2 .
(2)
строено приближенное аналитическое решение задачи о взаимодействии тонкой цилиндрической оболочки с грунтовой (упругой) средой при воз- действии на нее плоских сейсмических (нестацио- нарных) волн. Решение позволяет определить кольцевые и продольные деформации (и усилия) оболочке (рис. 1, 2).
По Ламе вектор смещений выражается через скалярный → и векторный потенциалы:
→ → → → ,
u u u grad rot
→
а уравнение (2) в форме:
1
2
2 с 2 2 / t 2 0 ; 2 z c2 2 / t 2 0 ;
В плоской постановке подобная задача реша- лось в [1, 2] для условия жесткого контакта обо-
2
r 2 c2 2 2 / t 2 0 , (3)
2
r 2 2
r r
лочки и среды, а в [3] для условий проскальзыва- ния. Расчетная схема задачи представлена на рис. 1. На круговую цилиндрическую оболочку с
2
/ r 2 2 / r r / c2 2 / t 2 0
радиусом R и толщиной стенки h из материала, характеризуемого плотностью S0 , коэффициентом
Уравнение для вектора смещений оболочки
0 u,,T принимается в виде [1]:
Пуассона и модулем Юнга Е, расположенную в
L u
2 2
2 2 2
2 (4)
упругой среде с плотностью S и коэффициентами
0 R D u0 / Dt
P/ 0h / c , c
E / 0 1 ,
Ламе
, , набегает под углом к оси оболочки
где: L – матрицы дифференциальных операторов.
плоская ступенчатая продольная волна напряже- ний:
L11 2 / 2 1 2 / 22 / 2 , L12
qr, , z, t q0H t
(1)
L 1 / 22 /
где: q0 - амплитуда набегающей волны; H - функ-
21
L L / , L
2 / 2
ция Хевсайда.
13 31 22
1 / 22 / 2 , L23 L32 / ,
2 2 4 / 4 22 / 2 2
L33 h
/12R
4 / 4 2 2 / 2 1
Т
1,
z , , г
ных сил.
– вектор поверхност-
Рис. 1. Расчетная схема
Начальные условия для полупростран- ства принимаем нулевыми; а граничные условия при r заданы набегающей волной (1), при r нулевые. На поверх- ности оболочки формулируются условия полного проскальзывания:
- непрерывность нормальных компо- нент напряжений при r = R:
r q0 1 sin 2 cos 2
x t rcos sin r sin z cos r
0hc 2 2 1 1 2 2
c1 R 2
r 2
R r
R 2 t 2
z 2
2
2
1 r
1 2 r
3
G q0 e psin D p, s, n, , In ,, r , , M , N /
r 2
R2
R
r z
/p / S p cos ,
]
– отсутствие касательных напряжений при r = R:
0 q 0 1 sin 2 sin cos
r
0hc 2 2 2 1 r 2
R 2 R r R r r 2
где: P,S - параметры преобразования Лапласа по t
=z/R;
n -параметр преобразования Фурье по ;
I u(x) - модифицированная цилиндрическая функция индекса n.
В уравнениях (5) использованы безразмерные
1 2 r
2 2
2 2 2
2
3 r
0;
переменные:
2 2
R
R z
R r z
/ R
; m m / R
; t C1t / R; r
1
r / R; q0 q0 / C 2; u u / h;
z 0,5q0 1 sin 2 cos
/ h; w w / h.
0
hc 2 2
2
2R rz
1 2 1 2
2
r
R z R r
Обращение изображений по t выполнялось с помощью теоремы запаздывания. Функция проги-
бов оболочки при = 0,3 в разделенных пере- менных получена в виде:
1 2
3 2
2
R r 2
1 3
r 3 R3
3 2
2 0,
wt,, z, q0E
5,46C10R2
(6)
R 2 r 2
r z
2 sin2 cos 1 sin2 2 sin sin3
4 u 3 ,
Непрерывность скоростей смещений при r = R:
t 1 3 H t 1 t 2 3 H t 2
w q
2 1,11,69 A cos 1,1cos 1,69A
0 sin cos
3 3
t C1
где: t 3 H t 3 t 4 H t 4 ;
1 r
2 r
0,65 A cos 0,65 cos A
t r
R
rz
t 4 H t t 4 H t
,
4
1 2 2
,
где:
4 1,11,69 A cos
1,1cos 1,69A
q0 C1 n ; n q0 / C1; n
t 3 4 H t 3 t 4 4 H t 4
q0 / asin cos
0,65A cos 0,65cos A
Решение сформулированной задачи строим ме- тодом разделения переменных. Применимы к уравнениям (3)-(4) преобразования Лаплас по вре- мени, конечное косинус - или синус - преобразо- вание Фурье по углу и преобразование Лапласа по координате z, записываем их в изображениях:
2 LcL / r 2 LcL / rr
n 2 / r 2 s 2 pcl LcL 0
-W/h 0,05
0,04
0,03
0,02
2 d / r 2 d / rr
(5)
0,01
m m
n 2 / r 2 s 2 p 2 / 2 d 0,
p m
1 2 3 4 5 6
где: m=1,2 d=LsL при m=1, d=LcL при m=2.
S 2 n2a U LcL nsb LsL SwLcL A2 p2u LcL ;
bsnuLcL n2 as2 LsL nwLcL A2 p2 LcL ;
suLcL nuLcL k * S2 n2 1 w LcL A2 p2w LcL ;
t tC1 / R
Do'stlaringiz bilan baham: |