Энергетическая модель механики в горном деле Текст научной статьи по специальности «Физика»

© Ю.А. Алюшин, 2014

УДК 531.8 + 624.131 Ю.А. Алюшин

ЭНЕРГЕТИЧЕСКАЯ МОДЕЛЬ МЕХАНИКИ В ГОРНОМ ДЕЛЕ

На основе представления приращения удельной энергии механической системы через изменения кинематических инвариантов уравнений движения в форме Лагранжа высказана гипотеза о существовании новых физических свойств грунтов, которые определяют механизм самоорганизации обратимых и необратимых деформаций, изменение механических свойств с учетом предшествующей истории деформирования и изменений внешних условий. Предложена энергетическая интерпретация понятия «обобщенная сила».

Ключевые слова: энергетическая модель механики, инварианты уравнений движения, механические и физические свойства грунтов.

Объектами переработки в горном деле могут быть самые разнообразные монолиты и сыпучие горные породы. Несмотря на сложную структуру грунтов, состоящих из пустот и отдельных частиц, если при действии внешних сил структурные связи между минеральными частицами не нарушаются, вполне применимы основные положения механики сплошных сред, в частности теории упругости и пластичности [1—2]. Такое предположение не приводит к большим погрешностям при использовании в расчетах механических свойств грунта: пределы прочности при сжатии и растяжении, модуль продольной упругости, коэффициент Пуассона. К примеру, решение задачи Бус-синеска [2] о распределении напряжений при действии на массив сосредоточенной силы используют для оценки устойчивости сооружений и разрушения грунтов.

Не отрицая ряда особенностей и закономерностей, обусловленных природой рыхлых горных пород как минерально-дисперсных образований, возможности нарушения структурных связей грунта, после чего деформации обусловлены главным образом перемещением его частиц, реологические свойства грунтов (релаксация, ползучесть), приводящие к прогрессирующему течению, развитию оползней, просадке насыпей и разрушению, хорошо согласуются с принятой в механике сплошных сред методикой комплексного изучения напряженного состояния среды (массива грунта) для анализа его прочности и устойчивости с учетом реально существующих фаз уплотнения, сдвига, потери устойчивости и разрушения.

В механике грунтов особую роль играет сжимаемость или уплотнение под действием внешних сжимающих нагрузок, высыхания, коагуляции коллоидов и пр. с уменьшением пористости грунта, в том числе за счет местных сдвигов, изменения толщины водно-коллоидных оболочек и пр. Но при этом всегда можно считать выполнимым основное предположение механики сплошных сред — ближайшие частицы остаются ближайшими на протяжении всего процесса движения.

Энергетическая модель механики основана на анализе кинематических инвариантов уравнений движения

^ (х ,ар, V = 0, (1)

где t — время, х 1 е (х, у, z), ар е (а, в, у) — переменные Эйлера и Лагранжа,

соответственно, которые несут всю информацию о внешних воздействиях и внутренних изменениях, происходящих в процессах движения (деформации). Система (1) может быть записана в различных формах, однако, необходимость учета истории деформирования, а также возможность использования принципа суперпозиции движений [3—5], особенно для сложных процессов, делает предпочтительной форму Лагранжа

X = х (ар, (2)

В дальнейшем в качестве переменных Лагранжа приняты начальные (при t = = 0) координаты точек в системе координат наблюдателя а = х |=0 , в = У |=0 ,

7 = г и .

В самом общем случае без каких-либо ограничений на свойства сплошной среды система (2) имеет 13 локальных кинематических инвариантов. Три из них связаны с векторными характеристиками движения: модули векторов перемещения и(ар, ^ , скорости Ща , () и ускорения Ф(ар, §

£ =М=ч/7, Цар,£ 1=л/7 , £ =1^.

Инвариантом также является путь э, равный интегралу от модуля скорости по времени

t

£ = э = |Ы Л .

0

Несимметричный тензор второго ранга, образуемый производными от переменных Эйлера по переменным Лагранжа

х,р = дх1 / дар , (3)

для которого в дальнейшем использован термин «тензор деформации Лагранжа» [3—5], имеет три инварианта

£ = Ха+ Ув+ Zr, (4)

£ = х2 + хв + х2 + У2 + ув + у2 + £ + 4 + 4, (5)

£ =\ х.р |=8V/ = Н. (6)

Кубический инвариант £ совпадает с якобианом преобразования (2) и равен отношению объемов бесконечно малой частицы в текущем SV и исходном состояниях.

В отличие от симметричного тензора деформаций Коши [6] инварианты (4— 6) всегда положительны, в исходном состоянии частицы принимают значения £5 = £6 = 3, £7 = Н = 1.

В соответствии с основным постулатом механики, поведение системы зависит от положения частиц и их скоростей. Тензор (3) можно рассматривать как обобщенные координаты, их скорости также образуют несимметричный тензор («обобщенные скорости») 312

х,р = дхц / дар, (7)

который имеет три инварианта: линейный, квадратичный и кубический

#8 = + У1Р + ^у = (#5^ ,

#9 = 4 + 4 + 4 + + -4 + у1 + 4 + + 4 , (8)

#10 =1 Х,р 1.

Дополнительно 3 инварианта могут быть получены интегрированием по времени модулей инвариантов #,#, #0. В отличие от инвариантов (4—6), которые в процессе деформации могут расти или уменьшаться, результаты интегрирования по времени

t

#11 = {V(ха + Л ,

0 í

#12 = {((+ 4 + 4 +уа + ув+у\ + 4 + 4 + 4),

0

#13 ={| хцр \А,

0

только возрастают на протяжении всего процесса деформирования и позволяют учесть историю этого процесса, уплотнение материала и пр.

Перечисленные 13 локальных инвариантов # являются независимыми, они или составленные из них выражения должны определять состояние и поведение частиц, а также механической системы в целом. Чтобы сравнивать состояния и предсказывать реакцию системы на внешние воздействия, приведенные выше 13 инвариантов надо привести к одному обобщенному скаляру, который Аристотелем [7] был назван энергией («обобщенная скалярная мера различных видов движения») ЗЕ = ЗЕ#).

Оператор З подчеркивает локальный (по отношению к пространству) характер скаляра.

Как показывает опыт, для большого класса механических систем из абсолютно твердых и деформируемых тел обобщенный скаляр можно представить в виде суммы составляющих, каждая из которых зависит только от одного инварианта ЗЕ = (#),

1

причем каждое слагаемое можно представить как произведение соответствующего инварианта на объем ЗЦ, и скалярный множитель к, характеризующий свойства среды и обеспечивающий равенство размерностей слагаемых,

ЗЕ = ^ЗЕ1 #) = £ к#З^. (9)

Скалярные коэффициенты к должны характеризовать либо физические свойства материала, например плотность материала при вычислении кинетической энергии частицы, либо свойства среды, в которой происходит движение,

например ускорение свободного падения при движении в гравитационном поле Земли. Гипотеза (9) может быть расширена, например, за счет учета взаимных влияний инвариантов, т.е. добавлением слагаемых, определяемых значениями двух и более инвариантов.

Дальнейший анализ ограничим формулировкой обобщенного закона движения механической системы в виде закона сохранения энергии на бесконечно малом интервале времени

dSE = SVo(d^к£) - SdEe = 0 . (10)

!

где SEe — энергия внешних воздействий. Оператор «а» — соответствует бесконечно малым приращениям функции во времени в отличие от оператора «S», используемого для бесконечно малых приращений функций в пространстве переменных Лагранжа.

Уравнение (10) предполагает определение бесконечно малых приращений энергии, которые могут быть вычислены на приращениях, выбранных для описания движения кинематических координат Яу, используемых в правых частях уравнений перечисленных выше инвариантов

dSE ((к£( д)) = S Vк = SQ¡JdqJ. (11)

дяу

Множитель Q(J по существу является энергетическим определением обобщенной локальной силы

SQ(J =дЕ = к £/0, (11а)

характеризующей скорость изменения соответствующего вида энергии Е( бесконечно малой частицы при изменении кинематического параметра Яу. В общем случае сила Qу может быть скаляром, если в качестве кинематического параметра выбран скаляр, например путь э или квадрат скорости V2, вектором, если параметры Яу являются проекциями вектора, или тензором 2 ранга. Размерность силы Q у также зависит от выбора кинематической координаты Яу, например [Н] или [Нм] для линейных или угловых перемещений, [Па] для тензора напряжений Лагранжа и пр.

Закон сохранения энергии (10) предполагает учет всех как внутренних, так и внешних энергетических факторов. Для учета энергетических потоков со стороны окружающих частиц воспользуемся общепринятой методикой, использующей скалярное произведение векторов силы SP и скорости V: dSEе = ^ ^Р ■ V) dt. Суммирование в правой части должно быть проведено по

всем ограничивающим рассматриваемую бесконечно малую частицу поверхностям. С учетом возможных изменений сил и скоростей на противоположных гранях, предполагая все функции дифференцируемыми и заданными в переменных Лагранжа, получим с точностью до бесконечно малых 1-го порядка (по пространственным переменным и времени),

а = = тр'х>р + хд / дар (12)

где т^ = (ЗРр1 / Зар )З\/0 — напряжения Ёагранжа, образуют несимметричный тензор второго ранга. Индекс р е (а, р, у) указывает направление нормали к рассматриваемой площадке в исходном состоянии, а индекс ! — направление проекции силы, может принимать значения ! е (х, у, г). Напряжения тр1 подобны напряжениям Пиола — Кирхгофа [6], но отличаются от них ограничением области изменения аргументов и неоднозначным выбором начала отсчета шкалы средних напряжений, которое может быть смещено относительно общепринятой [3—5].

С учетом внешних взаимодействий закон сохранения (10) можно записать в форме энергетического баланса

1ЗЕ = З\\0(к#1( + к2#2( + к3#3( + ... + к13#т -а)& = 0 . (10а)

Пренебрегая процессами диссипации (т. е. без учета инвариантов, связанных с интегрированием по времени), а также используя общепринятые соотношения для потенциальной и кинетической энергии (ось г направлена вертикально вверх)

с!ЗЕ1 = Р0дгЗ&, 1ЗЕ2 = р0(хгхп + уу + ггги)З\& , получим

1ЗЕ1 + 1ЗЕ2 + <ЗЗЕЪ + 1ЗЕ6 + с1ЗЕ7 = соЗ\/0 &,

или, с учетом дифференциальных уравнений движения,

ТрЛ,р = к5 (Х1а + У1Р+ Zty) +

+ 2кб( ХаХа + хрхвв + хух^ + УаУ а + УрУ1в + УуУь + г

+ к7(ХаХ1а + ХрХ1р + ХуХ* + УаУ(а + У'рУ\р + УуУ1у + 2а21а + + 2у2)у' .

(14)

Приравнивая коэффициенты при одинаковых множителях — компонентах тензора скорости деформации (7), получаем соотношения между компонентами напряжений, элементами тензора (3) и константами к!, характеризующими физические свойства материала,

Т = к5Зр1 + 2к6Х,р + к7Х,р . (16)

В соотношениях (16) и далее х,р — алгебраические дополнения элементов х1р матрицы (3), единичный тензор Зр1 принимает значения Зр1 = 1при соответствии индексов «р» и «!», т.е. Зр1 = 1 для тах, тру,тгг и Зр1 = 0 для всех остальных напряжений, не расположенных на главной диагонали. В исходном состоянии, когда переменные Эйлера и Ёагранжа совпадают (матрица якобиана преобразуется в единичную), компоненты тензора определяют только физические свойства

Тах = Тру = Туг = к5 + 2к6 + к7 , Трх = Тух = Тау = Туу = Тая = Тря = 0 .

Энергетический баланс должен выполняться в любой, в том числе начальный, момент времени, для которого можно использовать напряжения Коши а^

1ЗЕе /ЗЦЛ = а.х.а. (17)

Переходя в уравнении (17) от производных по переменным Эйлера dx. t / dXj = x. ^ к производным по переменным Лагранжа dxj t / dap = x. tp с помощью общих соотношений, вытекающих из уравнений движения (2) df / dxi = (df / dap)xx.ip / R, и приравнивая коэффициенты при одинаковых

множителях хПр в правых частях уравнений (12a) и (17), получим систему линейных уравнений

Tpi = ajixj,p, (18а) которые формально совпадают с известными статическими условиями на контуре и по существу определяют связи между напряжениями Лагранжа и Коши, справедливые для любого момента времени

ст.. = т x. /R. (18b)

j pi j,p' > '

Равенства (18b) можно трактовать как следствие условия инвариантности энергии по отношению к выбору начала отсчета времени в системе наблюдателя. С учетом (16) для напряжений ст. окончательно получаем

= R (aixi,a +Tpixj,p + T7ixi,7 ) = (18с)

= R lA xa + 2k6(xi,axja + x.ej + xj + k (Via + xj/x.,e + x^ij]

Сопоставление выражений (16) и (18с) показывает, что для анализа процессов деформации грунтов напряжения Лагранжа предпочтительнее: они энергетически обоснованы и связаны простыми математическими уравнениями с имеющими четкий геометрический смысл характеристиками деформированного состояния. Основной инвариантной характеристикой напряженного состояния можно считать среднее напряжение Коши ст

3ctR = къ£ + 2 ке£ + 3 к^, (19)

которое можно использовать для определения среднего напряжения в исходном состоянии

стД=0 = k5 + 2кб + к7. (20)

Если коэффициенты к5 — к7, характеризующие физические свойства деформируемого материала, известны, тогда по уравнениям движения в форме (2) можно определить любые кинематические, а затем и энергетические или силовые функции, в том числе напряжения Лагранжа (14) и Коши (18). Они могут быть использованы для корректировки (выбора) начала отсчета шкалы средних напряжений. Есть достаточно оснований считать, что в исходном состоянии средние напряжения не следует принимать равными 0. В частности, закон упругого изменения объема в механике деформируемого твердого тела ст = 3Ks можно считать совпадающим с законом изотермического расширения газа pV = const,, если модуль объемной упругости K рассматривать как действующее в текущем состоянии давление.

Из закона сохранения энергии в форме (10) следует, что деформация возможна при изменении не менее двух видов энергии (или работы внешних сил). Это позволяет установить связь между коэффициентами ki и привести систему отсчета различных видов энергии к одной шкале.

В качестве примера рассмотрим зависимость между коэффициентами к2 и к2 на примере свободного падения абсолютно твердого тела (частицы) в гравитационном поле Земли, в котором участвуют два вида энергии: потенциальная ЗёЕ1 и кинетическая З1Е2 . Сопротивлением воздуха пренебрегаем, иначе надо добавить изменение энергии З1Е4, предполагая какую-либо связь между диссипативными силами и инвариантом 5 или скоростью ¡VI. Уравнения движения и закон сохранения энергии примут вид (ось г направлена вертикально вверх)

х = а, у = р, г = у- иг(у,Г), 1ЗБ1 + 1ЗЕ2 = 0.

Повороты отсутствуют, движение поступательное, энергию можно проинтегрировать по всему объему тела. Для приращений энергии Е2 и Е2 следует записать

1Е1 = к1< 0, 1Е2 = к21(^2) = к21(22) = 2к2= 2к2г 1 > 0 ,

и, если использовать общепринятое обозначение для ускорения свободного падения г № =-д, соотношение между коэффициентами должно быть к1 = 2к2д. В классической механике принято к1 = тд, тогда к2 = т /2 и для кинетической энергии получаем общепринятое выражение Е2 = т^ /2.

Свойства, определяемые коэффициентами к5, к6, к7 , должны полностью определять энергетические изменения частиц в области обратимой деформации. Свойства к8, к9, к20 и инварианты #8, #9, #10 не вошли в уравнение энергетического баланса (14), учитывающего внешние воздействия, так как их производные по времени содержат множители типа хирр, которые не входят в выражение (12) для энергии внешних воздействий. Этого достаточно для утверждения, что они связаны с механизмами обратимой и необратимой деформации, учитывать их различие в глинах и сыпучих средах.

В процессе необратимой деформации условие энергетического баланса вместо (14) принимает вид

а = Тр1Х1,1р = к5#5 + к6#6,1 + к7#7,1 + к11#11,1 + к12#12,1 + к13#13,^, (27)

следовательно, диссипативные процессы связаны с параметрами #11, #12, #13 и

свойствами к22 , к22, к23. Именно они определяют упрочнение, изменение механических свойств грунтов, специфику горных пород, энергетические особенности нарушения структурных связей между минеральными частицами и пр.

Так как излагаемая энергетическая модель должна учитывать возможные варианты движения от любых внешних воздействий, рассмотрим изменение энергетического состояния частицы из изотропного материала при равномерном нагреве. В соответствии с общепринятыми представлениями и понятиями, при нагреве на температуру ДГ линейные размеры частицы изменяются на величину аГДГ, где аГ — коэффициент линейного расширения материала, при этом затрачивается энергия

Г

ДЗЕ = { сЗт1Т

Г>

или, приближенно, ДЗЕ = сср8тДТ, ДЗЕ/ З\/0 = ссрр0ДТ, где сср — средняя

теплоемкость материала в рассматриваемом диапазоне температур от Т0 до Т. С учетом уравнений движения х; = а(1 + атДГ) , деформаций Ёагранжа ха = (1 + ат ДГ) и приращений инвариантов Д/1 = 3ат ДТ, Д12 = 6ат ДТ + + 3(ат ДТ)2, Д13 = 3ат ДТ + 3(ат ДТ)2 + (ат ДТ)3, условие перехода подведенного тепла в энергию частицы принимает вид

ср0/ аТ = 3кь + 6к6 + 3к7 + 3(к6 + к7)аТДТ + к7(аТДТ)2 « 3(кь + 2к6 + к7). (33) Сравнивая правые части уравнений (20) и (33), можно утверждать, что средние напряжения в исходном состоянии следует считать равными

С | t=o = сРо /(3«т ) = къ + 2кб + к7, (34)

где все физические характеристики в правой части должны соответствовать исходному состоянию материала, т. е. при нормальном давлении и температуре 200С. По существу использование соотношения (34) соответствует переходу к новой энергетической шкале средних напряжений, по аналогии с термодинамической шкалой температуры Кельвина.

В отличие от твердых однородных тел, например металлов, нагрев грунтов сопровождается изменением, пористости, влажности и водопроницаемости. Поэтому в зависимости от диапазона температур уравнение (34) должно быть дополнено слагаемыми, учитывающими энергетические характеристики этих изменений.

Для определения каждого из коэффициентов правой части уравнения (20) достаточно дополнительно двух уравнений, например по результатам испытания на чистый сдвиг и гидростатическое сжатие. В качестве основного принимаем общее уравнение (19) для среднего напряжения 3с Я = к5£ + 2кб£ + 3А|,

которое можно привести к обычной шкале средних напряжений за счет сдвига шкалы на величину исходных напряжений (20), т.е. в обычной шкале зависимость среднего напряжения Коши от инвариантов тензора деформации принимает вид

ст = а-&0 = (к5£ + 2кб£ + 3к7^)/ 3 Я - (к5 + 2кб + к,) =

= к5 (" 1) + 2 к6 (3Я - 1 ] + к7 - 1

В условиях однородного гидростатического сжатия с уравнениями движения х = а(1 + = ае, I = 3е, | = 3е2, | = е3 и Я = | = е3 получаем сс = (1 - е)[к5(1 + е) + 2к6е]/ е2. С другой стороны, из закона упругого изменения объема

сс = 3К(Я-1/3 - 1) = 3К(е -1)/ е2. (36)

Приравнивая правые части последних двух уравнений, получим к5(1 + е) + 2 к6 е =-3 К (37)

или, принимая во внимание е = 1,

к5 + к6 =-1,5 К. (37а)

При чистом плоском сдвиге уравнения движения имеют вид

х = а + вр, у = ав + р, г = у,

где в — угол сдвига. Два инварианта сохраняют исходные значения, меняется только квадратичный инвариант с;5 = 11 = 3, с;6 = 12 = 3 + 2в2, с;7 = 13 = 1. Работа внешних сил ДЕе = ^тёу = С^уёу = 0,5 Су2 должна соответствовать измене-

7 7

нию энергии материала ДЕ6 = к6Д^6 = 2к6у2. Из энергетического баланса для обратимого процесса получаем

ке = 4 С. (38)

Важно, чтобы величина С была определена с помощью описанного эксперимента, а не вычислена через модуль Юнга и коэффициент Пуассона. Отрицательный знак коэффициента к5 объясняет увеличение энергии частицы при уменьшении ее объема за счет гидростатического сжатия.

Энергетическая модель механики позволяет обосновать существование новых физических свойств грунтов, возможность перехода к новой шкале средних напряжений, по аналогии с термодинамической шкалой температур, а также целесообразность исследования их влияния на изменение структурных связей между минеральными частицами, формирование дисперсных образований и пр.

- СПИСОК ЛИТЕРАТУРЫ

1. Малышев М.Б., Болдырев Г.Г. Механика грунтов. 2000-2010 гг. http://www.geo-teck.ru

2. Механика грунтов. Шаламанов В.А., Санников А.Ф., Крупина Н.В. и др. Электронный учебник. КузГТУ, 2004. - 112 с.: ил. IQlib.ru

3. Алюшин Ю.А. Энергетические основы механики. Учеб. пособие для вузов: - М.: Машиностроение, 1999. - 192с.

4. Алюшин Ю.А. Энергетическая модель обратимых и необратимых деформаций в пространстве переменных Ёагранжа. Сборник «Прогрессивные технологии пластической деформации». Москва, НИТУ МИСиС, 21-26 сентября 2009 года, стр. 44-67.

5. Алюшин Ю.А. Принцип суперпозиции движений в пространстве переменных Ёагранжа.// Проблемы машиностроения и надежности машин. 2001. № 3. С. 13-19.

6. Седов Ё.И. Механика сплошной среды. Т. 1 - М.: Наука, 1970. - 492 с.: ил.

7. Богомолов АН. Механика в истории человечества. - М.: Наука, 1978. - 150 е.: ил. н'.иа

КОРОТКО ОБ АВТОРЕ -

Алюшин Юрий Алексеевич, доктор технических наук, профессор, профессор кафедры «Теоретическая и прикладная механика» Московского государственного горного университета. E-mail: alyushin7@gmail.com.

^ UDC 531.8 + 624.131 ^-

POWER MODEL OF MECHANICS IN MINING

Alyushin Y.A., Doctor of Technical Sciences, Professor E-mail: alyushin7@gmail.com. Moscow State Mining University, Russia, ud@msmu.ru

On the basis of representation of an increment of specific energy of mechanical system through changes kinematic characteristics the equations of movement in the form of Lagranzh the hypothesis about existence of grounds new physical properties which define the mechanism of self-organizing reversible and irreversible deformations, change of mechanical properties in view of previous history of deformation and changes of external conditions is stated. Power interpretation of concept «the generalized force» is offered.

In soil mechanics plays a special role or compressibility of the seal under the external compressive loads, drying , coagulation of colloids , etc., with a decrease in soil porosity , including the expense of local shifts , changing the thickness of aqueous colloidal membranes , etc. But all - always can be considered feasible basic assumption of continuum mechanics - the coming of the particle remain closest throughout the process of movement.

Mechanics energy model allows to prove the existence of new physical properties of soils , the ability to transition to the new scale of medium voltage , by analogy with the thermodynamic temperature scale , and the feasibility study of their impact on the change in the structural links between the mineral particles , the formation of dispersed formations , etc.

Key words: Energetyc model of mechanics, the equations of movement, mechanical and physical properties terrestrial bowels

REFERENCES

1. Malyshev MV, Boldyrev GG Soil mechanics . 2000-2010 gg. http://www.geo - teck.ru

2 . Soil mechanics . Shalamanov VA, Sannikov AF, Krupin NV and other electronic textbook . KuzGTU 2004 . - 112 p. : 1l. IQlib.ru

3 . Alyushin YA Power bases of mechanics. Textbook. manual for schools : - M.: Mashinostroenie , 1999. - 192s .

4 . Alyushin YA Energy model of reversible and irreversible deformation in the space of Lagrangian variables . Collection of " Progress-tensive technology plastic deformation. " Moscow, NUST M1SA , 21-26 September 2009 , pp. 44-67 .

5 . Alyushin YA Superposition principle of motion in the space of variables Lagrange . / / Problems of engineering and machine reliability . 2001 . Number 3 . Pp. 13-19.

6. LI Sedov Continuum Mechanics . T. 1 - Moscow: Nauka , 1970 . - 492 s.: Il.

7. AN Bogomolov Mechanic in the history of mankind. - Moscow: Nauka , 1978 . - 150 s: Il.