книги из ГПНТБ / Золотарь, И. А. Экономико-математические методы в дорожном строительстве
.pdfРазличают еще задачи в условиях неопределенности, когда ве роятности различных результатов, которые возможны при выборе определенной стратегии, не могут быть количественно оценены. Та кие задачи рассматриваются в теории игр.
Таким образом, лишь в простейших случаях оказываются до статочными для отыскания оптимальных решений элементарные модели арифметики и алгебры. Примером этому могут служить за дачи с использованием критериев оптимальности в виде (1.2) и (1.6). В более сложных случаях приходится использовать специ альные математические методы и модели. При этом для оптимиза ции подсистем области организации дорожного строительства ока зываются наиболее эффективными: для подсистем планирования и управления — сетевые модели (см. гл. XI); для подсистем произ водства работ, материального обеспечения — модели экстремально го анализа (см. гл. II), исследования операций.
Ш е с т о й э т а п экономического анализа состоит в непосред ственном решении задачи отыскания оптимального решения на основе выбранной математической модели. Полученный результат анализируется для установления его особенностей, которые нужно будет отразить в проектных и технологических документах или в документах по организации работ.
Глава
МОДЕЛИ АРИФМЕТИКИ И АЛГЕБРЫ. ПРИМЕНЕНИЕ ЭКСТРЕМАЛЬНОГО АНАЛИЗА
В ЭКОНОМИЧЕСКОМ АНАЛИЗЕ ДОРОЖНОГО СТРОИТЕЛЬСТВА
§ 3. Применение арифметико-алгебраических моделей
Рассмотрим решение конкретного примера для дорог с коротким сроком службы, если в качестве критерия оптимальности исполь зуется срок окупаемости энергозатрат. Как ясно из самой структуры формулы (1.13), для ее использования необходима предварительная разработка технологических карт, из которых были бы видны затра ты машинного времени на выполнение всех видов работ.
Определение срока окупаемости энергозатрат на 1 км дороги в простейших случаях можно вести с помощью формулы (1.15). Если при этом не учитывать экономию в энергозатратах на текущий ре монт и содержание А5теКг и определять срок окупаемости дорож ных работ путем сравнения с вариантом оставления участка (доро ги) в его прежнем состоянии (Эд= 0 и учитывается лишь Э'д), то
П
|
v\v.r 2 Nur\utu |
|
|
или |
7 \ = ---------~^~1 к----------- ' |
(Н-1) |
|
|
К |
-■ » t) |
|
|
Расчеты по формуле (II. 1) |
7= 1 |
|
|
требуют лишь знания интенсивности |
и состава движения, его скоростей до и после проведения дорож ных работ, а также разработки технологической карты работ для
П |
Срок окупаемости Ть определяемый по |
вычисления v / VHr |
7= 1
формуле (II.1), является завышенным, так как не учитывает ни остаточной «стоимости» дорожных работ, ни экономии в ремонте и содержании. Тем не менее расчеты по формуле (II.1) показали, что величина Т\ даже для вариантов устройства на временных автомо бильных дорогах облегченных усовершенствованных покрытий обычно находится в пределах 1—5 сут, а это, как правило, значи тельно меньше требуемого срока службы этих дорог. Постройка таких покрытий экономически выгодна и должна планироваться во всех случаях, когда имеются необходимые материалы, силы и сред ства.
Применение формулы (II. 1) покажем на таком примере.
21
П р и м е р 1. При реконструкции автомобильной магистрали на одном из ее участков для пропуска транзитного движения в весеннее время используется объезд, проложенный по существующей автомобильной дороге с гравийным по крытием шириной 5 м. Необходимо установить технико-экономическую целесо образность уширения покрытия с 5 до 7 м при следующих данных: расчетная толщина покрытия на полосе уширения — 20 см; дальность возки гравийных ма териалов для уширения— 15 км.
Для выполнения работ имеются следующие машины: автомобили-самосвалы
ЗИЛ-ММЗ-555, экскаваторы Э-255, |
автогрейдеры |
Д-144А, прицепные |
катки с |
|||||
трактором ДТ-54. |
движения по объезду: автомобилей типа |
|||||||
Суточный |
состав транзитного |
|||||||
МАЗ-500— 100; |
ЗИЛ-130—2000; ГАЗ-52— 1500. Установленная |
средняя |
скорости |
|||||
движения на объезде до уширения |
оСр |
составляла |
12 |
км/ч, |
после |
уширения |
||
20 км/ч. Требуемый срок службы объезда |
Тд = 15 |
сут |
(по условиям реконструк |
|||||
ции). |
|
|
|
|
|
|
|
|
Составим сокращенную технологическую карту на |
необходимые |
дорожные |
||||||
работы на 1 км дороги (табл. 1). |
|
|
|
|
|
|
|
Работа
j Единица изме- 1 рения
Применяемые Количество машины
Т а б л и ц а |
1 |
|
Мощность двигателя, л. с. |
Потребность п машинном времени, ч |
Энергозатра ты, л. с. ч |
Заготовка в карьерах гравийных матери алов с погрузкой в самосвалы
Вывозка материалов на дорогу
Разравнивание мате риала
Укатка покрытия
Всего Непредвиденные знер-
гозатраты—10 %
И т о г о . . . .
Следовательно,
м3 |
1 0 00 -2,0Х |
Экскаватор |
48 |
15 |
720 |
|
х о ,2 0 -1,26= |
Э-302А |
|
|
|
=5 0 4
Я504 Автомобили-са 150 230 34500
|
|
мосвалы |
|
|
|
|
|
ЗИЛ-ММЗ-555 |
|
|
|
» |
504 |
Автогрейдер |
108 |
3 |
324 |
|
|
Д-144А |
|
|
|
М 2 |
2 000 |
Прицепной ка |
54 |
4,0 |
216 |
л. с. ч |
|
ток с ДТ-54 |
|
|
35 760 |
— |
— |
— |
— |
||
я |
— |
|
|
|
3 576 |
|
|
|
|
|
|39 336 |
4 |
|
|
|
|
|
21] Лг|(1)1с^( = 39336 л. с. ч/км.
1
Необходимо обратить внимание на то обстоятельство, что энергозатраты на транспортные работы составляют около 88% от суммарных энергозатрат. Это характерно для работ по устройству покрытий переходного типа.
Вычислим теперь энергохарактеристику движения:
з
V jv ja)^ a>= 100-180 + 2000-150 + 1500-115 = 490 500 л. с./сут.
P i
Найдем срок окупаемости Т\ с помощью формулы (11.1):
т20-12-39336
(20 - 12)490 500 |
J |
Таким образом, Т\<^ТД. Следовательно, эти работы весьма целесообразны.
22
Рассмотрим еще один пример, относящийся к частной оптими зации с использованием лишь арифметико-алгебраических моде лей.
П р и м е р 2. Требуется установить для лесовозных дорог в Вологодской обл.. оптимальную высоту насыпи, возводимой из местных суглинков, при которой, суммарная стоимость земляного полотна и гравийного дорожного покрытия бу дет минимальна. Требуемый эквивалентный модуль деформации £ Тр = 400 кгс/см2. Модуль деформации гравийного материала — 800 кгс/см2. Дальность возки гра вийного материала— 16—20 км. Увлажнение грунта насыпи весной и осеньк> может происходить вследствие застоя поверхностных вод (3-й тип местности). В табл. 2 и 3 приведены данные института Лесгипротранс о стоимости устрой ства насыпей различной высоты и гравийного покрытия при ширине земляного полотна 10 м.
|
|
|
|
|
Т а б л и ц а 2 |
|
|
|
Высота, м |
|
|
Тип местности |
0,5 |
0,75 |
!,0 |
1,5 |
2,0 |
но условиям |
|
|
|
|
|
увлажнения |
|
Стоимость 1 км земляного полотна, тыс. руб. |
|||
|
|
||||
1-Й |
4, 6 |
7,5 |
9, 9 |
15,7 |
22,3 |
2-й |
5,7 |
9,6 |
13,5 |
19,5 |
28,0 |
3-й |
7,3 |
12,2 |
17,1 |
25,0 |
35,4 |
На основе имеющихся методов прогнозирования модулей деформации зем ляного полотна 1 Е0, а также расчета необходимой толщины гравийного покрытия при £ тр=400 кгс/см2, найдены величины, приводимые в табл. 4.
|
Т а б л и ц а 3 |
|
Т а б л и ц а 4 |
||
Средняя дальность |
Стоимость 1 км |
Высота |
Модуль дефор |
Необходимая |
|
гравийного покрытия |
насыпи |
мации земля |
толщина |
||
возки гравийного |
толщиной 20 см, |
Я н, м |
ного полотна |
гравийного |
|
материала, км |
тыс. руб. |
£ 0, кгс/см2 |
покрытия А , м |
||
3 - 5 |
10,3 |
0,50 |
69 |
0,48 |
|
0,75 |
97 |
0,35 |
|||
6— 10 |
11,3 |
||||
1 ,0 0 |
131 |
0,31 |
|||
11—50 |
12,5 |
||||
1,50 |
178 |
0,24 |
|||
16—20 |
13,6 |
||||
2 ,0 0 |
191 |
0 ,2 2 |
|||
|
|
Экономическая функция, требующая минимизации, имеет в данном случае
вид:
V c = C3n + Cn. |
(II.2) |
Вычислим значения функции (II.2) при различной высоте насыпи, приведенной в табл. 2, и соответствующей им толщине Ап гравийного покрытия (см. табл. 4).
В результате получим: при высоте Нн=0,5 м
У 1, С = 7 ,3 + |
13,6 = 39,9 тыс. руб. |
^0,20
Вэтом расчете принято увеличение стоимости гравийного покрытия пропор ционально росту его толщины.
1 Водно-тепловой режим земляного |
полотна |
и дорожных одежд. Под ред„ |
И. А. Золотаря, Н. А. Пузакова, В. М. |
Сиденко. |
М., «Транспорт», 1971. 416 с. |
23
При Ян=0,75 м |
|
|
|
|
|
|
|
v, |
„ |
„ „ |
0,35 |
36,0 |
тыс. руб. |
|
у \ |
С = 1 2 ,2 + |
5 ^ 1 3 ,6 = |
|||
При Я н= |
1,0 м |
|
|
0 31 |
|
|
|
^ |
С = |
17,1 + |
38,2 тыс. руб. |
||
|
13,6 = |
|||||
При Я н= |
1,50 м |
|
|
О24 |
|
|
|
У |
С = |
|
41,3 |
тыс. руб. |
|
|
25,0 4- - L— 13,6 = |
|||||
|
— |
|
|
0,20 |
|
J |
При Я н=2,00 м |
|
|
|
|
|
|
|
^ |
С = |
_ |
0,22 |
50,4 |
тыс. руб. |
|
У |
35,4 -f —— 13,6 = |
||||
|
^ |
|
|
0,20 |
|
^ J |
Из полученных данных видно, |
что оптимальной |
является высота насыпи |
Я н= 0,75 м, при которой суммарная стоимость земляного полотна и дорожного покрытия минимальна (36 тыс. руб.).
Для отыскания оптимального решения в этом примере нам по надобились только простейшие арифметико-алгебраические дей
ствия. Если знать функциональные зависимости |
для |
кривых 1 |
||
и 2 на рис. 4, то |
можно искать |
|||
общее решение задачи метода |
||||
ми |
экстремального |
анализа, |
||
что будет показано ниже. |
||||
Здесь же |
мы |
прибегли к ва |
||
риантному сравнению, исчис |
||||
ляя |
для |
каждого |
варианта |
|
суммарную стоимость по алгеб |
||||
раическому |
|
соотношению |
Рис. 4. График к определению опти мальной высоты насыпи:
1 —кривая стоимости земляного полотна в зависимости от высоты насыпи: С3п=*/(ЯН);
2 — кривая стоимости гравийного покрытия
в зависимости от |
высоты насыпи: Сп= |
~ ф (Я н ); 3 — кривая |
суммарной стоимости |
земляного полотна и гравийного покрытия
хс
(П.2).
§ 4. Применение экстремаль ного анализа
Целью экономического ана лиза обычно является либо ми нимизация затрат ресурсов для достижения необходимого ре зультата, либо максимизация эффективности использования имеющихся ограниченных ре сурсов. Многие из подобных задач успешно решаются с по мощью так называемого пиедельного анализа, основываю щегося на нахождении экстре мумов шах или min соответст вующих функций методами дифференциального исчисле ния. В задачах математическо-
24
го анализа на отыскание экстремальных точек функции часто не существует каких-либо ограничений для ее переменных. Найденный таким путем шах или min функции называют безусловным. Рас смотрим этот метод отыскания оптимальных решений на несколь ких примерах.
П р и м е р 1. Установление необходимого типа и рационального радиуса действия производственных предприятий (заводов и баз) по приготовлению вяжущих и смесей для_ устройства покрытий усо вершенствованного типа.
По условиям использования в дорожном строительстве разли чают базы (заводы) стационарные (радиус действия Re более 100 км, длительность работы на одном месте Т не менее 10 лет),
полустационарные (^б = 40—60 км; 7 = 3 —5 лет) |
и |
передвижные |
или полевые (7?б = 5— 10 км). При строительстве |
автомобильных |
|
дорог в настоящее время наиболее распространены |
полустацио |
|
нарные базы и заводы. |
|
|
Необходимый тип и рациональный радиус действия базы могут
быть определены на основе соотношения |
|
С = Ж + СП+ СТР, |
(П.З) |
где С — стоимость продукции базы, отнесенная к 1 |
м2 устраиваемо |
го покрытия; М — часть этой стоимости, не зависящая от искомого радиуса действия Re базы; Сп — стоимость, обусловленная расхо дами по систематической передислокации базы через каждые Re километров построенного покрытия; Стр— часть стоимости, обу словленная транспортными затратами на перевозку продукции ба зы (из расчета на 1 м2 покрытия) на среднее расстояние Rе/2 ки
лометров (рис. |
5). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
б) |
|
|
|
^_ |
ц |
р |
1*— |
йр |
|
пj] |
Gj] |
|
% |
/ -— |
' |
|
|
Rg |
Rfi |
P |
|||
. Ktf |
|
|
. «б |
. |
|
|
|
||
Рис. 5. Схема к определению |
среднего |
расстояния |
перевозки продукции |
базы: |
|||||
а — при вывозке |
продукции |
базы в |
обе стороны |
от места |
расположения; б — при |
вывозке |
|||
в одну сторону («от себя»); |
I; II; |
I I I —места |
последовательного размещения |
базы. |
Стрел |
||||
ками показаны направления |
вывозки продукции |
базы; R§ — радиус действия |
базы. |
Протя |
женность обеспечиваемого с каждой стоянки участка — 2Rg в схеме а и R q в схеме б
Величина М в основном отражает стоимость исходных материа
лов Мм и процесса их переработки Мв на базе для выпуска |
ее про |
|
дукции в количестве, требуемом на 1 |
м2 покрытия, т. е. |
|
М = М м + |
Мп. |
(II.4) |
Величину Мм можно считать не зависящей от Rе. При вычисле нии радиуса действия Re для завода с определенными оборудова нием и технологией слагаемое Ми также не зависит от Re. Однако задача может ставиться и в более общем плане, когда требуется установить как Re, так и тип базы (стационарная, полустационар-
25
лая, передвижная). Тогда Ma= f ( R б). Характер этой функции мож но выявить на основе данных стоимости приготовления вяжущих (смесей) на базах различных типов. Примем в первом приближе нии следующее выражение для определения величины Мп:
M„ = aaRe6 + CB, |
(Н.5) |
где а0— коэффициент; в — показатель степени. |
(энерго |
В зависимости (II.5) Св представляет собой стоимость |
затраты) приготовления вяжущего (в количестве, требуемом для устройства 1 м2 покрытия) на базе, перемещающейся вместе с фронтом работ (Re = 0) *. Для стационарных и полустационарных баз величина Ма должна быть меньше, чем для передвижных, по этому коэффициент f l o < 0 .
|
Внося (II.5) в (II.4), получим М = Мж+а^Я1 |
+ Сб, |
или |
М = айЯ1^\-с1, |
(II.6) |
где |
d = M u-\-Ca. |
|
Величина Сп в зависимости (II.3) должна представлять стои мость демонтажа базы, ее перевозки на расстояние Re и л и 2 Re (см.
рис. |
5, а и б) и монтажа на новом месте дислокации, отнесенную к |
|||
1 м2 |
покрытия. Таким образом, |
|
|
|
|
Q |
б д .м |
ткм |
/JJ |
|
п— |
кюомЯб |
|
|
где |
Сд.м — суммарная стоимость |
демонтажа и монтажа |
базы; |
|
сТкм — средняя стоимость 1 |
ткм |
перевозок оборудования |
базы; |
Qe— общая масса оборудования базы; k — коэффициент, равный 2
или 1 в зависимости от схемы размещения базы по трассе (см. рис. 5, а, б); Ьа -— ширина покрытия, м.
В общем случае величина Сд.м также зависит от Re. Примем для Сд.м тот же вид зависимости, что и для величины Ма (см. фор мулу П.5), т. е.
Сди= / о - ( / ? 6- / ? ^ + ^д.«, |
(П-8) |
где /о — коэффициент; g — показатель степени; сдм— стоимость де монтажа и монтажа передвижной базы с радиусом действия R'е.
Наконец, величина Стр может быть вычислена из соотношения
С тр |
Re |
(П.9) |
|
2s |
|||
|
|
||
где с'ткм — стоимость 1 ткм перевозок |
продукции базы; s — пло |
щадь покрытия, м2, которое может быть построено при использо вании 1 т продукции базы.
* Так как таких баз еще нет, то величину Св можно принимать для пере движных баз с ^?б'=5—10 км. Тогда в формулу (11.5) вместо Re следует вво дить R6— Rg.
26
Используя зависимости (II.6), (II.7), (II.8) и (II.9), получим на основе соотношения (II.3) следующее выражение:
Г ' |
( г , |
—v'\e I . I / о ( * б — |
+ с д м + 9 б ^^ бСткм , с ткм ^6 |
C = a 0{ R 6 — R 6) -\-d-\------------------- |
------------------------------- 1----- ш— • |
Если принять e —g = \, |
то можно найти формулу для определе |
||||
ния рационального значения R§: |
|
|
|
||
dC |
— |
СД.М /оЛб |
I |
^ТКМ |
г\ |
-------— CLc\ |
1000й„й^ |
+ |
2s |
— U, |
|
dRb |
|
|
|
откуда
В частном случае, полагая величины Мп и Сдм не зависящими от радиуса действия базы, будем иметь а.о = 0 и fo=0. Тогда из фор мулы (11.10) получаем известную зависимость
R 6 = |
Сд m^S |
(П.11) |
|
1000i„feeTKM |
|||
|
|
где под величиной сд.м следует понимать затраты на монтаж и де монтаж конкретной базы, для которой ищется рациональный ра диус ее действия. Очевидно, что зависимость (11.10) имеет более общий характер, позволяет определить радиус действия базы, а по его величине и тип базы (стационарная, полустационарная, пере движная) в соответствии с приведенной выше классификацией. Од ной из важных задач является уточнение коэффициентов в форму лах (II.6) и (II.8) на основе изучения технико-экономических по казателей работы баз различного типа.
Так, применительно к |
устройству |
покрытия |
толщиной |
5 |
см из |
гравийного |
|
материала, обрабатываемого жидким |
битумом |
методом |
смешения |
на дороге, |
|||
примем в формуле (11.11) |
сд.м= 5 |
тыс. руб.; с'Т1Ш=0,06 |
тыс. руб.; |
s= 200 м2/т |
|||
(исходя из расчета 5 л битума на 1 |
м2 покрытия); k = l ; bn= 7 |
м. |
|
||||
5000-2-200 |
|
|
|
|
|
|
/ |
----------------------- яэ70 км. |
1000-7-1-0,06 |
|
|
П р и м е р 2. Определение оптимальной длины захватки при ве дении дорожно-строительных работ поточным методом и использо вании ведущих машин, совершающих поворот (реверс) в конце каждого рабочего прохода.
Специфика решения этого вопроса обусловлена тем, что подоб ные машины должны в конце каждого рабочего прохода делать поворот с неизбежной для этого потерей времени. Следовательно, удлинение захватки будет увеличивать производительность маши ны. С другой стороны, увеличение длины захватки может обусло
27
вить затяжку периода развертывания потока, а следовательно, и удлинение общего срока выполнения работ! Очевидно, что может быть найдена оптимальная длина захватки /опт, при которой срок постройки покрытия на участке L будет минимальным.
Для поточной организации работ справедливо соотношение
(^раб ^разв/ Ам/7Д,е i к•> |
(II. 1,2) |
где Траб— число рабочих дней в общем периоде Т, выделенном для строительства дороги протяжением L; /разв— продолжительность периода развертывания комплексного потока, рабочие дни; kCM— коэффициент сменности; П — производительность ведущей машины в тех же единицах, что и L; пвед — число ведущих машин.
Легко уяснить, что
/раз.= -Ц"- - ')- , |
(11-13) |
где / — длина захватки; п — число частных потоков в составе комп лексного.
Удлинение захватки сокращает число смен с выдачей готовой продукции, хотя и способствует росту производительности машин вследствие уменьшения доли времени, затрачиваемой на повороты после каждого рабочего прохода. Зависимость (11.12) полностью отражает это противоречие, и потому из нее можно будет найти ве личину /опт-
Для этого внесем (II.13) в (11.12). В результате получим
Трао |
Ц п - |
1) |
кгиПпв |
--L. |
|
fin„- |
] |
|
|
Далее получаем |
|
|
|
|
^раб^см^^вед |
/ (^ |
^) ^см |
к * |
L +1(п — 1) kc
откуда Граб--
Подставим сюда вместо П следующее соотношение:
П = -------- |
--------------, |
(II. 14) |
/ 2 1 |
\ |
У |
( |
+ 2^пов Щ |
|
\ «раб |
/ |
|
где ураб — скорость машины во время рабочего прохода; /пов— вре мя на поворот в конце каждого рабочего прохода; то — число про ходов для получения готовой продукции; къ— коэффициент исполь зования рабочего времени.
Тогда получим:
-г |
2L |
+ In |
2 (я — 1) |
+ In Щ |
_ |
( «рай |
«раб |
|
|
■* раб |
' |
^см^веАм^Д |
Лве-АмА |
|
|
|
|
28
|
2m,о |
|
|
п — 1)= |
|
|
^вел^см^в ^раб |
|
|||
|
|
|
|||
|
2mn Г |
L |
Un |
Ц П - |
1) |
|
^раб^с] |
&см^ |
|
■^„ов(Я— U |
|
|
^раб |
||||
Чтобы найти оптимальную длину захватки, при которой потре |
|||||
буется |
минимальное |
число рабочих дней |
Граб для строительства |
||
дороги |
протяжением L, |
найдем |
d T раб |
|
|
—— и приравняем затем произ- |
|||||
водную нулю: |
|
|
dl |
|
|
|
|
|
|
||
|
d T раб |
LL |
п — 1 |
||
|
|
dl |
|
^раб |
=о, |
|
|
|
|
||
откуда |
|
|
^пов ^раб |
(11.15) |
|
|
|
|
|
При наличии в составе комплексного потока /ср резервных за хваток и /Сф захваток, обеспечивающих технологические разрывы, формула (11.15) примет вид:
1 = 1/ ____ ипов- ^ ------- |
. |
(11.16) |
у(п + /ср + /Сф— 1) кш
Как видно из приведенного решения, получение формулы (11.15) потребовало достаточно сложных математических преобразований.
Рассмотрим упрощенный метод решения задачи на основе ана лиза размерностей. В соответствии с приемами последнего запишем очевидную теперь функциональную зависимость:
*опт = / Ков! ^раб; L;(n—1)], |
(П.17) |
где tnoB— время на разворот машины в конце рабочего прохода, осуществляемого на скорости ораб1 п — число частных потоков в со ставе комплексного.
Характеристика (п— 1) отражает связь длин захваток с перио дом развертывания потока.
Перепишем (11.17) в виде:
/опт= 4 . ^ 1й( 7 з Т )Т. |
(II-18) |
Это соотношение следует из так называемой П-теоремы мето да анализа размерностей и говорит о том, что искомая величина может быть выражена как произведение определяющих ее величин в неизвестных нам степенях.
Составим теперь на основе (11.18) уравнения размерностей для входящих величин. По размерности длины получим следующее уравнение: 1 = р+у.
По размерности времени будем иметь: 0 = а — {5.
Из полученных двух уравнений вытекает третье ^,= а + у.
29