Вопросы теории пластичности
РАСЧЕТ ПЛАСТИЧЕСКИХ ДЕФОРМАЦИЙ МАТЕРИАЛОВ И
ГРУНТОВ ДОРОЖНЫХ КОНСТРУКЦИЙ ПРИ ВОЗДЕЙСТВИИ ТРАНСПОРТНОЙ НАГРУЗКИ
А.С. АЛЕКСАНДРОВ, канд. техн. наук, доцент
Сибирская государственная автомобильно-дорожная академия (СибАДИ) г. Омск, Россия
Исследованию физической природы пластического деформирования материалов и грунтов дорожных конструкций посвящено большое количество публикаций крупнейших ученых. Работы О.Т. Батракова, А.К. Бируля, И.А. Золотаря, Н.Н. Иванова, М.С. Коганзона, В.Д. Казарновского, Я.А. Калужского, В.П. Матуа, С.И. Миховича, А.В. Смирнова, Н.Я. Хархуты и др. являются фундаментальными в дорожной отрасли.
Анализ результатов исследований этих ученых показывает, что в большей или меньшей степени авторы в своих работах исходят из ряда допущений:
1. Материалы и грунты рассматриваются как идеально-пластические тела, испытывающие остаточные деформации при возникновении сколь угодно малых напряжений.
2. Взаимосвязь пластической деформации с величиной напряжения описывается линейными функциями при сколь угодно большом напряжении.
3. Продолжительность воздействия нагрузки принимается 0,1 с и не учитывает скорость движения транспортных средств, показатели механических свойств материалов и грунтов, толщину конструктивных слоев и т.п.
4. Влияние повторности (многократности) приложения нагрузок на величину накапливаемой пластической деформации выполняется одним из двух способов. Первый способ основан на использовании эмпирических логарифмических функций и по большому счету обоснован только для связных грунтов. В основу второго способа положена гипотеза подобия деформирования материалов при воздействии циклических нагрузок процессу ползучести при однократном воздействии длительной нагрузки.
Отмеченные допущения позволяют сформулировать цель работы, заключающуюся в физико-математическом описании процесса пластического деформирования материалов и грунтов в условиях воздействия повторных, кратковременных нагрузок для разработки методики проектирования дорожных конструкций по критерию обеспечения ровности покрытий автомобильных дорог.
Поставленная цель, по мнению автора, может быть достигнута последовательным решением трех задач:
1. Разработка метода определения линейных и нелинейных пластических деформаций при однократном приложении кратковременной нагрузки.
2. Разработка метода определения продолжительности напряженного состояния в различных элементах дорожной конструкции с учетом скорости движения нагрузки, показателей механических свойств материалов и толщины конструктивных слоев.
3. Аналитический вывод и обоснование функциональных зависимостей величины накапливаемой пластической деформации от количества прикладываемых нагрузок.
Решение первой задачи можно выполнить в рамках технической реологии, в которой разработаны разнообразные модели и теории (ползучести, старения,
наследственности, упрочнения, консолидации и т.д.) [1]. По мнению автора наилучшими являются решения, отображающие реологические процессы при помощи так называемых структурных сопротивлений деформированию моделируемого объекта. В зависимости от материала и его физического состояния, а так же свойств, для которых необходимо выполнить математическое описание реологических процессов структурных сопротивлений может быть различное количество [2, 3]. Физический смысл структурного сопротивления заключается в том, что если значение напряжения превысит эту характеристику на сколь угодно малую величину, то характер деформирования изменится. Ранее автор осуществлял попытку совершенствования физической модели упруговязкопла-стического тела, состоящую из последовательно соединенных сложных реологических тел. Каждое из этих сложных тел [4] представляет собой известную или усовершенствованную модель и обладает своим структурным сопротивлением. В диапазоне напряжений, при возникновении которых материал испытывает пластические деформации, предполагается наличие 4 сопротивлений [4].
Приведем наименование и пояснение каждого из этих сопротивлений в порядке возрастания их количественного значения.
Предел начальной структурной прочности (предел обратимости деформаций) роб - величина, ограничивающая сверху множество значений напряжений вертикального сжатия, при возникновении которых материал проявляет свойства упруговязкого тела, то есть претерпевает только обратимые деформации.
Предел структурной вязкости рл - величина, ограничивающая сверху множество значений напряжений вертикального сжатия, при возникновении которых материал проявляет свойства линейного упруговязкопластического тела, деформируясь как обратимо, так и пластически.
Предел линейности пластических деформаций ре - введенная авторами [4] величина, ограничивающая сверху множество значений напряжений вертикального сжатия, при возникновении которых материал проявляет свойства нелинейного упруговязкопластического тела, с вязкопластической составляющей деформации, зависящей от величины напряжения. Таким образом, вязкопласти-ческая составляющая деформации с излишком напряжения (а-рл) связаны нелинейной зависимостью, а мгновенная пластическая деформация остается пропорциональной излишку напряжения.
Предел текучести рт - величина, ограничивающая сверху множество значений напряжений вертикального сжатия, при возникновении которых материал проявляет свойства нелинейного упруговязкопластического тела, испытывая мгновенные пластические и вязкопластические составляющие деформации, зависящие от величины напряжения. Зависимость между излишком напряжения (с-ре) и составляющими пластической деформации нелинейная. В случае, если напряжение превысит предел текучести, материал течет.
В табл. 1 приведены, полученные формулы для расчета пластической деформации при однократном приложении кратковременной нагрузки [4].
Из табл. 1 следует, что вязкопластические деформации рассчитываются через модуль общей равновесной деформации, но не от всего излишка напряжения, а только от релаксирующей части. Физический смысл такого решения заключается в том, что отрелаксировавшая часть излишка напряжения обуславливает вязкопластическую деформацию. Коэффициент ^р характеризует часть напряжения, расходуемого на мгновенные упругую и пластическую деформацию и деформацию упругого последействия.
Таблица 1. Формулы для расчета пластических деформаций при _однократном приложении кратковременной нагрузки
Диапазон напряжений
Формула
Роб < ^ р
4t)
= Роб )-(l-Ц 2 )
(
Е,
■ +
(1 - Кр)
Е
др
1 - ехр -
-1
об
Т,
Mt)
(РЛ -РобИ1-Ц2)
Роб ^ ^ РЕ
+ (ст- Рл)-(1 -Ц2 )•
Е.,
- +
(l - К )
Е
др
1 - ехр-
-1
об
тг
+
Е
+
(l - Кр)• К
Е
др
1 - ехр
- L
тг
Srrn(t)
(р, - роб М1 - М2 )•
роб рп
+(ре- р^м1 -Ц2 )•
(1 - к)
-+
Е Ед
у мп др
(1 - kV )• К,
1 - ехр -
-1
об
Т.
-+
Емп Едр
1 - ехр
- t„
Т,
+
+
+(°-ре)^(1 -Ц2 )•
К
Е
+
(1 - Крз )• К а
Е
др
1 - te
1 - ехр—-Т
р
где Еум - модуль упругой мгновенной деформации (модуль Юнга), МПа; Еув - модуль упруговязкой деформации, МПа; Емп - модуль мгновенной пластической деформации, МПа; Едр - модуль общей равновесной деформации, МПа; ц - коэффициент Пуассона; ^р - коэффициент, характеризующий нерелаксирующую часть напряжения; tv и tE -соответственно эквивалентные продолжительности напряженного состояния с напряжениями, превышающими пределы обратимости деформаций, структурной вязкости и линейности пластических деформаций, с; Тр - время релаксации соответствующих излишков напряжений, с.
Коэффициент можно выразить через соотношение составляющих полной деформации. После простейших преобразований имеем:
k - Е
^нр ~ др
f 1 1 1 ^
-+-+-
Е Е Е
у ум ув мп у
(1)
Коэффициенты К, и KtE определяется по зависимостям, аналогичным формуле проф. Н.Н. Маслова
Т
К, - 1 + ln ^ т
(г|-е)
(
- 1 + ln
Z,
"(л-е)
Z„
Ке - 1 + 1пТ(е-т) - 1 + ln
zr
•(е-т)
Z„
(2)
где Zr и Тл - соответственно глубина зоны распространения пластической деформации и время, необходимое на ее стабилизации при передаче давления, равного пределу структурной вязкости, м и с; Zafr.s) и ТГ(Г.е) - соответственно глубина зоны распространения пластической деформации и время, необходимое на ее стабилизации при передаче давления, изменяющегося в диапазоне от до рЕ, м и с; Zafe.ro) и Тг(е.т) - соответственно глубина зоны распространения
1
8
1
В
1
1
V
V
пластической деформации и время, необходимое на ее стабилизации при передаче давления, изменяющегося в диапазоне отре до рт, м и с.
При решении второй задачи процесс воздействия транспортной нагрузки на дорожную конструкцию во времени можно разделить на три стадии:
1 - нагружение рассматриваемого сечения. В течение этой стадии напряжения вертикального сжатия возрастают от нуля до своей максимальной величины, что обусловлено приближением автомобиля к рассматриваемому сечению;
2 - условно можно считать, что напряжение некоторое время остается постоянным, имея максимальное значение. Это имеет место при проезде автомобиля непосредственно по рассматриваемому сечению. Продолжительность этой стадии характеризуется скоростью движения транспортного средства и длиной контакта колеса с поверхностью покрытия;
3 - в течение этой стадии напряжения вертикального сжатия уменьшаются от максимального до нулевого значения, что обусловлено удалением автомобиля от рассматриваемого сечения.
На рис. представлена схема для определения продолжительности напряженного состояния эквивалентной времени воздействия транспортной нагрузки на расчетное сечение проезжей части.
tэкв2
CD
с
Отах
0
s
Т
g роб
а
Q.
с
CD
X
Ai
Ь/
Fi
В
I стадия - уваЛичение напряжения от нуля до максимального Значения
Hi
Ki
tin
С
II стадия -напряжение имеет максимальное значения
Н
стадия - уменьшение напряжения от максимального значения до нуля
tэкв1
Время, с
Рисунок. Изменение во времени напряженного состояния в рассматриваемой точке
На этом рисунке кривой АВСD показан характер изменения напряжений во времени в рассматриваемом сечении: увеличение напряжения от нуля до максимума (кривая АВ), условно постоянное максимальное напряжение (отрезок ВС), уменьшение напряжения от максимума до нуля (кривая СD). Согласно предложению проф. Н.Я. Хархуты [1] продолжительность напряженного состояния с непостоянным во времени значением напряжения, можно заменить эквивалентной продолжительностью напряженного состояния с постоянным во времени напряжением, имеющим, например, максимальную величину смах. Такая эквивалентная продолжительность для жесткопластического тела определяется из условия равенства площадей фигур А1ВСD1 и А\F\H\K\, то есть в основу решения положим теорему о среднеинтегральном значении функции:
% «от - Роб = ) - Роб К + | (а _ - р0б + ) ) - р0б . (3)
0 0 ^ 2
где а(0 - функциональные зависимости изменения напряжения вертикального сжатия во времени; t1, t2, t3 и 4кв1 - пределы интегрирования функций. 6
F
0
Из выражения (3) следует, что необходимо установить функции изменения напряжения во времени в течение 1-й и 3-й стадий процесса воздействия нагрузки, т. е. функцию увеличения напряжения во времени от нуля до максимума и уменьшения его во времени от максимума до нуля.
В работе [5] показано, что изменение напряжения вертикального сжатия во времени описывается формулой
(
^ )= Р •
1 +
2 • г
г
(4)
D У
где р - давление, передаваемое колесом автомобиля на покрытие, МПа; г-время с; То - период (время), в течение которого колесо транспортного средства перемещается на расстояние D, равное диаметру отпечатка колеса (диаметру штампа), м.
С учетом зависимости (4) и симметричности функций с(^) и с(^) интегральное выражение (3) примет вид:
мах " Роб )• гэкв1 мах - Роб - г1 ) +
+ 1
0
(
Р •
' 2 •
1 + -
Т
'Роб
D У
' 2 •
1+
Т
Роб
D У
(5)
После интегрирования выражение (5) примет окончательный вид
D
г _ г ___
1экв 1 об
Р • D
(
З (р - Роб )• 2
( (
1 +
V V
-1
Р • D
(
(Р - Роб )• 2 •З
V
( (
1 +
V V
-1
'Роб
Роб 1 Л -1
1
-1
/У
Роб • ^
(Р - Роб )• 2 •З
-1
Роб
У
Роб • D
(р - Роб )• 2 •З
-1
Роб
Аналогично получим
(
_ о -
(р-Рл)^2• з•
Р^о
V
\
( (
1 +
V V
Р -1
-1
-1
/У
Рл •0 (Р-Рл)2•З ^
.V
Р -1
(
р • о (Р - Рл> 2 •З
V
(
1 +
V V
Р -1
о
и _-—
р • о
(
З (р - реУ 2 •З
л
((
1 +
V V
1
-1
У
Рл •0
(р - рлУ 2 •З
.К
Р -1
Л
\Р -1
-1
-1
р • о
(Р - 2 •З
(
1 +
V V
\Р -1
У
ут1 ^ -1
УУ
р^ о
л У (
р^ о
(Р - РеУ 2 •З
(Р - Р^ 2 •З
(
л Р,
Р -1
Р -1
(6)
(7)
(8)
В настоящее время для учета повторяемости приложения нагрузок используется два принципиально разных способа. Первый способ основан на описании накопления пластической деформации от повторяющихся нагрузок эмпирическими формулами, как правило, логарифмическими зависимостями, некоторые из которых представлены в табл. 2.
2
2
л
л
Второй способ основан на гипотезе о подобии процессов ползучести и накопления деформации в результате многократного воздействия циклической нагрузки. Основные формулы, базирующиеся на этой гипотезе, приведены в табл. 3.
Таблица 2. Логарифмические зависимости, учитывающие повторность приложения нагрузки при расчете пластической деформации
Авторы
Формула
Н.Н. Иванов
^ = а + р-
где а и р - параметры, зависящие от вида материала или разновидности грунта, величины нагрузки, частоты приложений, продолжительности цикла, скорости нарастания деформации; N -количество реализованных одинаковых нагрузок, ед..
¿К
где - прогиб дорожной одежды, возникающий при первом приложении нагрузки, мм, р - давление под колесом автомобиля, МПа; D - диаметр круга, равновеликого по площади отпечатку следа колеса расчетного автомобиля, м; Еобщ - общий модуль упругости дорожной конструкции, МПа.
Коллективы, разработавшие ВСН 46-60, ВСН 46-72, ВСН 46-83
А.В. Смирнов
^ = ¿1 • Ып/^К где Б1 - остаточная деформация слоя дорожной одежды или активной зоны земляного полотна, возникающая при первом приложении нагрузки, мм; N - требуемое число, которое должен выдержать материал до разрушения, ед.; п - реализованное число напряжений, ед.
- w •(! + Kq • lgN)
В.С. Орловский
w
N " 1 V ' "9
где V1 - высота уступа между плитами сборного покрытия, возникающая при первом приложения нагрузки, мм; К9 - коэффициент, учитывающий уровень напряженного состояния основания под плитой.
В.Д. Казарновский, А.С. Пилипенко
¿К = ¿1 +
где Б1 - остаточная деформация грунта, возникающая при первом приложении нагрузки, мм.
А.С. Александров, Н.П.
Александрова, В.В. Голубенко
(
N
Sn - S1 1 + ß 1 lg Nj +$ulgNf + ß
N
jjj
jjj'
N
+ß
N
\
jv '
N
у " J 1' jj 1' jjj у
где S1 - деформация втапливания каменного зерна поверхностной обработки в покрытие при первом приложении нагрузки, мм; Nj -NjV - предельные значения количества реализованных нагрузок, по завершению которых происходит изменение интенсивности накопления пластической деформации.
В работе [6] автором показано, что процессы пластического деформирования при воздействии длительной и многократной нагрузки имеют различные формы протекания релаксационных процессов и не являются подобными. Поэтому обоснование функциональных зависимостей пластической деформации от количества реализованных нагрузок является актуальной задачей.
Для вывода формулы, позволяющей производить количественную оценку пластической деформации, накапливаемой упруговязкопластическим материалом при многократном воздействии нагрузок, предположим, что этот процесс носит наследственный характер. Тогда аналогично теории наследственной ползучести Больцмана-Вольтера можно записать интегральное выражение для оп-
ределения пластической деформации, накапливаемой при многократном воздействии нагрузки. При этом необходимо определить функциональную зависимость приращения деформации от числа приложения нагрузки.
Таблица 3. Формулы расчета пластической деформации при использовании гипотезы подобия
Автор Формула
Н.Я. Хархута ¿К = К • км + 5вп •1п (1 + х • tЭкв •К)] где ¿пм и 8вп - соответственно пластическая мгновенная и вязкопла-стическая деформации при первом воздействии нагрузки; X - показатель изменения вязкости грунта во времени; ?экв - эквивалентная продолжительность напряженного состояния; К - коэффициент (К = 1,4).
А.П. Васильев, М.С. Коганзон, и др. ¿К = Эвп ■1п (1 + X • Кгнд ) где Кс - общее количество приложения нагрузки и суточная интенсивность движения расчетных автомобилей; ?ц - время действия одной транспортной нагрузки; Кгнд - коэффициент, учитывающий неравномерное распределение в течение года движения, приведенного к расчетной осевой нагрузке.
С.Ю. Каныгина Г Т-К -К ^+0,13 •В] = -[1 + Т ц Ксум где - остаточная деформация от однократного приложения нагрузки (по исследованиям С.Ю. Каныгиной являющаяся функцией величины касательного напряжения и коэффициента консистенции грунта); Ксум - коэффициент суммирования нагрузок; /0 - множитель размерности 1с; В - коэффициент консистенции грунта.
В.Б. Фадеев Г Т-К -К ^[0,2+0,13 -В] 0 Л 1 1У с 1ц "сум ¿К = 51-Ь0- 1 + / ' ^ *0 ]
Е.В. Жустарева Г Т М Г К Л[0,2+0,13-/(Ку)• 0,13 • В] с ц сум ¿К = 51- 1 +---— 1 *0 ] где /(Ку) -функция, учитывающая влияние плотности связного грунта на процесс накопления остаточных деформаций;
Базируясь на результатах эксперимента, Ю.М. Гусев [7] пришел к выводу, что остаточная деформация, наблюдаемая при п-ом воздействии нагрузки £, взаимосвязана с деформацией от первого воздействия и описывается формулой
5 = •п ~1,5, (9)
Исходя из этого, в качестве ядра интегрального выражения примем подобные степенные функции, а именно:
Де„„ = а-к-п; Дк = а-к-п, (10)
мп ' вп ' V /
Демп = па к; Девп = пак, (11)
где Аемп и Аевп - соответственно приращения мгновенной пластической и вязко-пластической деформаций от п-го воздействия нагрузки; а - коэффициент, учитывающий вид материала; к - коэффициент, учитывающий уровень напряженного состояния.
Составив интегральное выражение и выполнив его интегрирование с учетом (10) и (11), получим:
N 1
1 + а • к •/пп _ (8мя1 +8впП) + а • к • ЫЩ, (12)
8 N _ (8 мп1 + 8 впП ) '
8 N _ (8 мп1 + 8 впП ) •
1
N
1 +1 nа•кdn 1
_ (8 мп1 + 8 впП ) •
" (Nа•к+1 -1)
1 + Л-
а • к +1
(13)
где емп1 и евп1— соответственно мгновенная пластическая и вязкопластическая деформации от первого воздействия нагрузки, определяемые по соответствующим слагаемым в формулах табл. 1.
Согласно исследованиям СибАДИ, коэффициент к можно определить, используя структурные сопротивления, по формуле
к _ ст-Ро6 -1, (14)
'Роб
Из анализа (12) и (13) следует, что эти формулы являются достаточно простыми и в тоже время обобщающими решениями полученных ранее функциональных зависимостей табл. 1.
В зависимости от соотношения излишков напряжения в формуле (15) зависимость (13) способна описывать затухающий, установившийся и прогрессирующий характер накопления пластической деформации и соответствует экспериментальным данным [8, 9].
В заключении, следует отметить:
1. Формулы для расчета деформаций позволяют описывать кусочно-линейные и нелинейные зависимости деформаций от величины напряжения, что подтверждается многочисленными испытаниями образцов из различных разновидностей грунтов [1-4, 7] и разных материалов [8, 9];
2. По мнению автора одним из преимуществ данного подхода являются учет всего спектра упруговязкопластических свойств материалов и грунтов, а так же возможность перехода от модели с шестью структурными сопротивлениями к моделям с меньшим количеством сопротивлений. Это позволяет кардинально менять характер деформирования модели. Достигается такой переход равенством соответствующих структурных сопротивлений друг другу;
3. Экспериментальные лабораторные исследования, выполненные в кандидатских диссертациях В.В. Голубенко, Н.П. Александровой и Н.В. Кузина показывают, что формулы таблицы 1 с учетом (13) удовлетворительно описывают деформации образцов из различных асфальтобетонов и снега (плотностью 0,30,6 т/м3). Это позволяет надеяться, что в будущем становится возможной разработка единой методики для прогнозирования изменения ровности дорожных конструкций со слоями из различных материалов;
4. Предлагаемые формулы могут быть положены в основу вывода формул для расчета перемещений (абсолютных деформаций) конструктивных слоев дорожной одежды и активной зоны земляного полотна. Для этого формулы таблицы 3 необходимо проинтегрировать по глубине слоев дорожной одежды и земляного полотна. Вид формул расчета пластических перемещений будет зависеть от выбранной функции изменения напряжения вертикального сжатия по глубине слоя или полупространства.
В работе [4] приведены формулы для расчета пластических осадок (перемещений) земляного полотна при воздействии повторяющихся нагрузок. Кроме того, интегрирование формул табл. 1 можно произвести различными приближенными методами, например использованием квадратурных формул (методы 10
трапеций, Симпсона, Ньютона-Котеса и т.д.), позволяющих вычислять определенные собственные интегралы сложных функций.
В этом случае при расчете перемещений появляется возможность учитывать изменение по глубине не только напряжения вертикального сжатия, но и показателей реологических свойств, зависящих от характера изменения влажности грунтов в земляном полотне и температуры в конструктивных слоях одежды из материалов и грунтов, обработанных органическим вяжущим.
Таким образом, задача расчета пластических перемещений, накапливаемых дорожной конструкцией, имеет ряд альтернативных решений;
5. Использование функциональной зависимости (13) при расчете пластических деформаций материалов и грунтов дорожных конструкций при воздействии многократно прикладываемых (повторяющихся) нагрузок дает возможность описывать затухающий, установившийся и прогрессирующий характер накопления деформации, что подтверждается многочисленными испытаниями дорожных конструкций на кольцевых стендах [8, 9].
Л и т е р а т у р а
1. Вялов С.С. Реологические основы механики грунтов [Текст] / С.С. Вялов. - М.: Изд-во Высшая школа, 1978. - 417 с.
2. Гольдштейн М.Н. Механические свойства грунтов [Текст] / М.Н. Гольдштейн. -М.: Стройиздат, 1971. - 367 с.
3. Александров А.С. Моделирование деформационных процессов, протекающих в связных грунтах [Текст] /А.С. Александров // Наука и техника в дорожной отрасли,
2002. - № 4. - С. 16-19.
4. Смирнов А.В. Динамическая устойчивость и расчет дорожных конструкций [Текст] /А.В. Смирнов, С.К. Илиополов, А.С. Александров - Омск: Изд-во СибАДИ,
2003. - 188 с.
5. Александров А.С. Определение продолжительности напряженного состояния в элементах дорожной конструкции при воздействии подвижных нагрузок [Текст] /А.С. Александров, Н.В. Кузин //Транспортное строительство. - 2008. - № 2. - С.24-28.
6. Александров А.С. Нелинейное пластическое деформирование при воздействии повторных кратковременных нагрузок /А.С. Александров //Известия ВУЗов. Строительство. - 2008. - № 10.
7. Гусев Ю.М. Остаточные деформации грунтов в строительстве [Текст] / Ю.М. Гусев. - Киев - Донецк: - Изд-во Вища школа, 1980. - 88 с.
8. Кусков В.Н. Прогноз срока службы дорожных одежд на основе стендовых испытаний. [Текст] /В.Н. Кусков //Автореф. дис... к-та техн. наук. - Тюмень: Изд-во Ротапринт Гипротюменьнефтегаза, 1996 - 21 с.
9. Смирнов А.В. Теоретические и экспериментальные исследования работоспособности нежестких дорожных одежд [Текст] /А.В. Смирнов //Дис. .. д-ра техн. наук, Омск, 1989. - 389 с.
AN ANALYSIS OF PLASTIC DEFORMATION OF MATERIALS AND SOILS OF ROAD STRUCTURES CAUSED BY TRANSPORT LOADING
A.S. Alexandrov
The physical model of "plastic deformations - time span and quantity of loads" is presented. The process of accumulation of plastic deformations has inherited nature under multiple actions.
The integral expression for the determination of plastic deformation caused by multiple force actions limited by Bolzman and Walter theory was derived and solved.