УДК 539.375 + 622.235
DOI 10.33764/2618-981X-2022-2-4-29-37
Расчетная схема динамики взаимодействия клиновидного ударника и породного массива в трехмерной постановке
Е. Н. Шер1*
1 Федеральное государственное бюджетное учреждение науки «Институт горного дела им. Н.А. Чинакала» СО РАН, г. Новосибирск, Российская Федерация
*e-mail: ensher@gmail.com
Аннотация. Представлены варианты расчетных схем динамики внедрения клиновидного инструмента при ударе в горную пароду. Расчеты производятся в квазистатическом приближении оценок сил, препятствующих внедрению, полученных ранее в двухмерной и трехмерной постановках. Рассмотрены случаи упругопластического деформирования среды и упругого с учетом равновесного развития магистральной трещины. Результаты расчетов сравниваются с данными модельного эксперимента по ударному внедрению клина в блок из оргстекла. Проведенное сравнение подтверждает работоспособность разработанных схем расчета.
Ключевые слова: удар, клин, горные породы, динамика внедрение клина, магистральная трещина
3D computational model of wedge penetrator-rock mass
interaction dynamics
E. N. Sher1*
1Chinakal Institute of Mining, Siberian Branch, Russian Academy of Sciences, Novosibirsk,
Russian Federation *e-mail: ensher@gmail.com
Abstract. The paper describes computational model variants of penetration dynamics of a wedge-shaped tool in rock mass. The computation uses the quasi-static approximation of penetration resistance forces obtained preliminarily in 2D and 3D formulations. The scope of the analysis embraces the elastoplastic and elastic deformation events in the medium with regard to the equilibrium propagation of a main crack. The calculation results are compared with the experimental modeling data on impact-driven wedge penetration in an organic glass block. The comparison proves efficiency of the described computational models.
Keywords: blow, wedge, rocks, wedge penetration dynamics, main crack
Введение
Ударное разрушение широко используется в существующих и перспективных технологиях горного производства: при отбойке полезных ископаемых (навесные молоты [1], экскаваторы с ковшами активного действия для открытых работ [2, 3] и динамические струги для подземных работ [4]), при ударно-вращательном способе бурения [5]. Актуальным для совершенствования таких технологий является разработка расчетных методов оценки разрушающего действия ударного органа механизмов на породный массив. К настоящему времени на основе анализа большого объема экспериментальных исследований [6] выявлена
качественная картина процесса разрушения при ударе жесткого индентора по свободной границе твердой среды. Отмечается, что в начале внедрения индентора возле точки контакта происходит образование зоны пластического течения. В случае с клиновидным индентором возникающее пластическое течение приводит к образованию медианной трещины, находящейся в плоскости симметрии. Дальнейшее увеличение нагрузки на индентор приводит к ее устойчивому росту [7-11]. Его разрушающее действие подробно исследовано в плоской постановке задачи о взаимодействии клина и хрупкой горной породы в [12,13]. В этих работах с использованием аппарата механики трещин разработаны расчетные модели определения параметров развивающейся магистральной трещины при внедрении падающего массивного клина и динамики его движения. В [14] рассмотрена в трехмерной постановке задача о внедрении клина с ограниченной шириной лезвия в упругое полупространство нормально к его поверхности. Рассчитана форма образующейся плоской трещины, имеющая вид полуовала, описанного вокруг кромки внедряющегося клина. Аналогичный результат получен при экспериментальном моделировании на образцах из оргстекла. Для расчета движения инструмента горных машин с клиновидным индентером и развития магистральной трещины, определяющей разрушение в хрупкой горной породе, важное значение имеет коэффициент жесткости породы при внедрении клина. Такое исследование было проведено в [15]. В настоящей статье с использованием полученных ранее результатов разработаны расчетные схемы динамики внедрения клиновидного индентора в упругую среду, позволяющие определять глубину внедрения индентора и размер образующейся трещины в зависимости от начальной скорости индентора, геометрических его параметров и механических свойств породы.
Задача определения силовых параметров действующих на внедряющийся клин на этапе пластического течения была рассмотрена с учетом трения в жестко-пластической постановке в работе [16]. При внедрении клина среда выдавливается по обе его стороны, при этом картина деформации будет иметь вид, схематически изображенный на рис. 1. Среда в области АВЭЕС находится в пластическом состоянии. Граничная линия АС аппроксимируется прямой. Условие текучести имеет вид
(ах -оу )2 + 4т] = 4т],,
где т, — предел текучести,
Области ЛВБ и ЛЕС двигаются как целое, в них постоянное напряженное состояние, в треугольнике АОЕ центрированное. Параметрами, определяющими задачу, являются: глубина внедрения клина к, полуугол его заострения у, коэффициент трения клина о породу предел текучести среды к = т,.
С учетом условий непрерывности напряжений при переходе через границы выделенных областей и сохранения объема удается построить решение при постоянном распределении нормального напряжения р вдоль клина на интервале
АВ и определить углы ю, ф и длины /, /1,Ь . При этом сила, внедряющая клин определяется по формуле:
Г = 2Ьр(//к) сов у^у + = С0к,
(1)
где Ь ширина клина.
Рис. 1 Схема внедрения клина в упругопластическую среду
Результаты расчетов силы пластического внедрения клина в зависимости от полуугла его приострения и коэффициента трения приведены на рис.2
Рис. 2 Зависимость силы пластического внедрения клина в зависимости от угла приострения 2 у и коэффициента трения ^ - кривые 1-4 соответствуют значениям коэффициента трения: 0, 0.1, 0.2, 0.3
Найденная величина ¥(к) позволяет определить динамику внедрения клина массы т с начальной скоростью у0в квазистатическом приближении:
тк = ~Сок.
Решение этого уравнения с учетом начальных условий:
к() = (?) в^ОО, о2 = Я (2)
Расчеты динамики торможения клина на этапе развития магистральной трещины
В работе [15] были разработаны в трехмерной постановке две модели расчета силы сопротивления внедрения клина. В модели 1 считалось, что свободные от напряжений берега трещины отходят от вершины клина, в модели 2 от уровня начального заглубления клина на этапе пластического внедрения ксг.
Согласно модели 1 сила сопротивления внедрению клина с хорошей точностью определяется формулами:
^ = С^2, С = 2пъ,\<£{у)к{у))4ъ /А, к(г,м) = '?{7) . (3)
1 -) ( у)
2
Здесь А = 32(1 -V )/(3Е) , где V - коэффициент Пуассона, Е - модуль Юнга горной породы.
Движение клина массы т, имеющего предударную скорость у0 описывается уравнением Ньютона
тк = -С1кШ . (4)
Отсюда следует следующая зависимость скорости движения клина от глубины внедрения
V (к Ь^Г^5'2
V 5 т
Подставляя сюда нулевое значение скорости в момент остановки клина, получаем максимальную величину внедрения клина.
к =
тах
Г с \2/5 5т
4С1 У
*4/5. (5)
Время до остановки клина определяется выражением
к 1
t _ тах |
0 О
dx
тах
V,
Г-
х
5/2
к
1 47 тах
Vrl
По модели 2 в [15] на первом этапе внедрения клин преодолевает сопротивление пластически деформирующейся среды, двигаясь в соответствии с (2) до критической глубины внедрения ксг при достижении которой начинается развитие магистральной трещины. Из закона движения (2) следует, что это происходит
1 . ксга
в момент времени tcr _— агсБт
а
сг
V,
при котором скорость клина составляет
v1 _ v0cos(аtcr). Начиная с этого момента сила сопротивления внедрению клина определяется выражениями
^ _ С^(2И - ИсГ), С _ 2^3/21в(Г)^(Г,М)^Ь / А.
(7)
С учетом этого закон движения клина в соответствии с законом Ньютона
2С,
d 2 х
определяется уравнением:
_-С2х, где х _ к - Нсг /2, С2 _ —-JhCr
т
с начальными условиями при t =0: х _ ксг /2, dx / dt _ v1. Таким условиям удовлетворяет соотношение:
к^) _ ксг /2 + V-sin(а>t) + cos(аt). а 2
(8)
Отсюда максимальное заглубление клина
Итах ксг /2 +
\
г V! Л
+
к
сг 2
(9)
Расчет по этой формуле для клина массой 2.3 кг и скоростью 2,2 м/с при внедрении клина в блок оргстекла дает Итах _ 0,193 см, общее время внедрения составляет 1.5 мс
В работе [12] разработана в плоской упругопластической постановке модель расчета динамики внедрения клина с учетом изменения размера контактной площадки соприкосновения граней клина с породой в условиях развитой магистральной трещины. Известно решение о деформировании твердой среды, ограниченной гранями прямого угла в области контакта воздействием на нее плоским, расположенным перпендикулярно биссектрисе угла штампом. В задаче о
внедрении клина такое нагружение осуществляется при угле заострения клина 90о. Сила действия клина на угол вдоль его биссектрисы пропорцианальна размеру площадки смятия 2а при пластическом течении: / = (4 + ж)аЬк. Для угла
заострения клина 2 у , отличного от 90о предложено в этом соотношении в качестве / использовать проекцию контактных сил, действующих со стороны клина, на биссектрису прямого угла:
/ = ЯЯ , Я =
1
2^2
1 +
1
к ( У, М),
где Я - сила внедрения клина. При этом размер проекции площадки смятия на ось клина
И0 = С3 Я /(Ьк), где С3
(1 + к ( у,м))ео5( у ) (4 + Жк/2 к ( у,м)
(10)
Аналогично (7) имеем: Я = С1ЛУИ0(2И - Ид), исключая отсюда Я с учетом (10), получаем квадратное уравнение для определения Ио:
И0Ь2к2 / С32 = С12(4И2 - 4ИИ0 + И02),
решением которого является выражение:
И0 = Ир + 2И Ир + 4ИрИ , где Ир =
/ Л 2
Ьк
V С1С3)
Отсюда и из (4) следует уравнение движения клина массой т :
с12к Ьк
А2 тС3
Ир +2И+4ИрИ
Интегрируя это уравнение, получаем зависимость скорости движения V клина от глубины внедрения И
V( И)
1
/
2Ьк
тС3
2 Ир ИрИ + И2 —р р 6
V
Г Л3/2
, 4И 1 + —
V р)
(11)
Для подтверждения работоспособности представленных моделей расчета динамики внедрения клина в горную породу необходимо сравнить их результаты с экспериментальными. Некоторые результаты таких модельных экспериментов представлены в работе [14], в которой исследовалось внедрение падающего клина в образец из оргстекла. Использование прозрачного оргстекла позволяло фиксировать глубину внедрения клина, размер и форму образующейся трещины. Результаты внедрения падающего клина массой 2.3 кг, шириной лезвия 28мм и углом приострения 30° в образец из оргстекла определялись на копровой установке. Регистрация сигнала акселерометра, установленного на клине, позволяла определить силу его торможения во времени. Так при предударной скорости 220 см /с внедрение клина составило 0,2 см, максимальная сила торможения - 6310 Н, длительность сигнала акселерометра - 2мс. Также в экспериментах было обнаружено, что при скорости клина , меньшей чем 100см/с магистральная трещина не возникает. Для проведения сравнительных с этим экспериментом расчетов принимались следующие значения параметров клина и среды: т=2.3 кг,
Ь =2.8 см, у_ 15° , ¡и_ 0.2, т8 _ к _ 35МПа, у_ 0.3, Е _ 3 • 109 Па ,
К1с _ 1МПа • м05. Рассчитаем по этим данным параметры внедрения клина с использованием предложенных в настоящей работе моделей.
По модели упругопластического внедрения движение клина определятся соотношением (2). Принимая коэффициент трения ^=0.2 и предел текучести оргстекла к _ 35МПа по графику рис.2 получаем С0 _3.1 -106Н/м . Отсюда а_ 11571/ с и максимальное заглубление клина ктах _ 1.9мм При этом, ^ _ 5850Н , полупериод колебаний 2.7 мс. Эти данные достаточно близки к экспериментальным. Отметим здесь, что при предельной скорости 100 см/с расчетное заглубление клина в упругопластическом режиме составляет 0.86 мм, которое в расчетах по модели 2 было принято за критическую глубину, начиная с которой происходит развитие магистральной трещины.
По модели 1 интегрирование уравнения Ньютона (4) позволяет определить по формулам (5) и (6) максимальное заглубление клина Итах _ 2 мм и время его достижения tmax _ 1.34 мс. Из соотношения (3): F _ 6860Н,
По модели 2 согласно (8) и (9) Итах _ 0,193, общее время внедрения составляет 1.5 мс.
По модели 3 расчетами по (11) была построена зависимость скорости клина от глубины внедрения, приведенная на рис.3. Из графика этого рисунка видно, что клин останавливается при ктах =0.2 см. Соответствующая этому значению внедрения сила ровна ^ =6760 Н
Из приведенных для разных моделей результатов видно, что они достаточно близки к экспериментальным данным.
Для выбора лучшей модели требуется более тщательная экспериментальная проверка для разных горных пород и параметров ударного нагружения.
и
Рис.3 Зависимость скорости внедрения клина от глубины внедрения, рассчитанная по модели 3 (формула 11)
Заключение
Разработаны расчетные схемы для определения параметров ударного внедрения клиновидного инструмента в горную породу для четырех моделей деформирования среды. Среди них рассмотрено внедрение клина в модели упругопла-стического деформирования среды, упругого и упругогопластического деформирования в условиях роста магистральной трещины. По полученным решениям рассчитаны характерные параметры движения клина применительно к прове-денныму ранее эксперименту по внедрения стального клиновидного ударника в блок из оргстекла. По всем рассмотренным моделям получено удовлетворительное согласие расчетных данных, что свидетельствует о работоспособности рассматриваемых моделей. Для более детальной оценки моделей требуется экспериментальная их проверка для разных горных пород и параметров ударного нагружения.
Благодарности
Работа выполнена в рамках проекта НИР (гос. регистр. № 121062200075-4).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Ушаков Л. С, Котылев Ю. Е., Кравченко В. А. Гидравлические машины ударного действия. — М.: Машиностроение, 2000. — 415 с.
2. Безвзрывные технологии открытой добычи твердых полезных ископаемых / Маттис А. Р., Ческидов В. И., Лабутин В. Н. и др. — Новосибирск: Изд-во СО РАН, 2007. — 335 с.
3. Городилов Л. В., Лабутин В. Н. Перспективы создания ковшей активного действия к гидравлическим строительным экскаваторам // Материалы V Междунар. науч. симп. "Ударно-вибрационные системы, машины и технологии", 23 - 25 апреля 2013 г. — Орел: Орел-ГТУ, 2013. — С. 112 - 119.
4. Горбунов В. Ф., Лазуткин А. Г., Ушаков Л. С. Импульсный гидропривод горных машин. — Новосибирск: Наука, 1986. — 195 с.
5. Крюков Г. М. Физика разрушения горных пород при бурении и взрывании. — Т. 1. — М.: Горная книга, 2006. — 330 с.
6. Колесников Ю. В., Морозов Е. М. Механика контактного разрушения. — М.: ЛКИ, 2010. — 224 с.
7. Lawn B. R., Evans A. G., and Marshall D. B. Elastic / plastic in indentation damage in ceramics: the median/radial crack system, J. Amer. Ceram. Soc, 1980, Vol. 63, No. 9 - 10. — P. 574 - 581.
8. Lawn B. R. and Wilshav T. R. Indentation fracture: principles and application, J. Mater. Sci., 1975, Vol. 10, No. 6. — P. 1049 - 1081.
9. Lawn B. R. and Swain M. V. Microfracture beneath point indentation in brittle solids, J. Mater. Sci., 1975, Vol. 10, No. 1. — P. 113 - 122.
10. Maurer W. C. and Rinehart J. S. Impact crater formation in rock, J. Appl. Phys., 1960, Vol. 31, No. 7. — P. 1247.
11. Paul B. and Sikarskie D. A preliminary theory on static penetration by a ridge wedge into brittle rock material, Transaction of SME-AIME, Dec. 1965. — P. 372 - 383.
12. Башеев Г. В., Ефимов В. В., Мартынюк П. А. Расчетная модель разрушения горных пород клиновидным ударным инструментом // ФТПРПИ. — 1999. — № 5. — С. 53 - 61.
13. Башеев Г. В. Расчетная модель откола куска горной породы при ударе клином под уступ // ФТПРПИ. — 2004. — № 5. — С. 77 - 89.
14. Шер Е. Н., Ефимов В. П. Трехмерное моделирование развития трещины в твердом теле при внедрении жесткого клина // ФТПРПИ. — 2015. — № 6. — С. 43 - 48.
15. Шер Е.Н. Численная оценка сопротивления внедрению клиновидного инструмента в хрупкий породный массив с учетом равновесного развития магистральной трещины // ФТПРПИ. —2021. — №6
16. Башеев Г.В., Ефимов В.П., Мартынюк П.А. Образование зародышевой трещины при ударном разрушении хрупкого тела клином// ФТПРПИ. - 1999. - №1. - С 52-59
© Е. Н. Шер, 2022