Материал: Миндели, Э. О. Разрушение горных пород учебное пособие

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

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

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

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

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

Покажем на примере планирование эксперимента по определению, напри­ мер, удельного расхода ВВ q при проведении горной выработки в зависимости от трех факторов: S — сечения выработки; / — коэффициента крепости пересека­ емых пород по шкале проф. М. М. Протодьяконова; Ф — коэффициента работо­

способности ВВ. Задача эта по своей цели является интерполяционной, т. е. в данном случае не требуется получить оптимум, а необходимо лишь предсказать результат с требуемой точностью в некоторой заранее заданной области. Реше­ ние задачи планирования будем производить в общем виде.

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

Представим

неизвестную

функцию

отклика полиномом

первой степени

в виде

 

 

 

 

или для данной

у = 6о"1- -f-b<iXo-\- . • •Ц-bjXj,

(XX .5)

задачи

 

 

 

 

q =

ho+ b 1S 1+

b 2 f+ b 30 .

(XX.6)

Задача заключается в нахождении коэффициентов Ъо, . .

Ь/ с достаточной

точностью при минимально возможном количестве опытов и условии адекват­ ности выбранной линейной модели.

Предположим, что мы установили область определения целевой функции и факторов. План полного факторного эксперимента представим в виде матрицы планирования типа 23. Здесь Число 2 означает, что все факторы варьируются на двух уровнях (условие линейности модели), тогда в соответствии с формулой (XX.4) необходимо поставить 7V = га *23 = 8ге опытов. Уровни факторов для удобства записи и обработки результатов эксперимента представляют обычно в кодированном виде с помощью матрицы (табл. 6.9) или формулы преобразо­ вания

Xj— X]о

xi

~ 7 ] ~ ’

(XX .7)

 

 

где xj — кодированное значение фактора; ж/ — натуральное значение

фактора;

я/о — натуральное

значение основного уровня, относительно которого верхний

и нижний уровни

симметричны, обозначаются соответственно +1 и —1; 7/ —

21 Заказ 1162

321

интервал варьирования, представляющий расстояние на координатной оси между основным и одним из граничных уровней (от 0 до +1 или от 0 до —1);

/ — номер фактора.

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

(табл. 70).

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

ном эксперименте оценивается перемножением знаков столбцов и получением таким образом новой матрицы с добавлением столбца эффекта взаимодействия факторов (табл. 71).

Прп добавлении эффекта взаимодействия согласно теории планирования

эксперимента свойства

матриц

планирования

не

меняются.

Модель (XX.6)

в этом случае принимает впд

 

 

 

 

 

g = b0+

b1S + boJ+ blt

Sf.

 

(X X .8)

Предположим, что

эффект

взаимодействия

Sf

оказался

малозначимым,

т. е. bll2 близок к 0 и, следовательно, линейная модель (XX.8) вполне адекватна, тогда столбец Sf можно использовать для третьего фактора Ф. В этом случае эффект взаимодействия Sf заменяется линейным эффектом от Ф и новая матрица

принимает вид, показанный в табл. 72.

Л1 опыта

S

/

ф

5

Л1! опыта

S

1

 

_

+

 

5

_

2

—

+

—

Я »

6

—

3

+

—

—

Я 3

7

+

4

+

+

+

?4

8

+

Л: опыта

S

/

ч

Л'! опыта

S

1

 

 

 

3

1

2

—

+

Яг

4

 

"7“

■Ni опыта S / Sf я Л" опыта S

1

 

+

+

?1

3

+

2

—

—

Яг

4

+

Т а б л и ц а

69

/

ф

я

_

+

Яъ

+

—

Яз

—

—

<77

 

+

+

?8

Т а б л и ц а

70

/

<7

+

Яз

?»

Т а б л и ц а

71

/

s i

Я

Яз

++

322

 

 

 

 

 

 

 

 

Т а б л и ц а

72

Л1' опыта

S

f

ф

Я

JMi опыта

S

/

ф

я

1

 

+

+

 

3

+

+

+

Яз

2

 

 

?2

4

+

Я4

Таким образом,

вместо восьми опытов (при п =

1) для определения q

нам

достаточно поставить четыре опыта, достигая этим самым сокращения опытов в 2 раза. Следовательно, достаточно в уравнении (XX.8) заменить Ь1-25/ на Ь3Ф , чтобы получить искомую функцию отклика от трех факторов при мини­ мально возможном числе опытов.

Определение коэффициентов 6/ производим по формуле

N

2

(XX. 9>

где Л’ — число опытов; у ц — значение целевой функции; х,-[ — значение уровня /-того фактора в i-том опыте; / = 0 , 1 , 2 , . . . , А.

Для свободного члена 60 / = 0 и тогда

N

Для данной задачи:

I/O-

лг

i-0

4

N

 

2

 

 

Ьо =-*=£-■

(XX .10)

N

41

 

2

 

1=о

__ 7 1 + 7 2 + 7 3 + ?4 .

 

4

4

 

( — )? 1+ ( —') 72+ ( + ) Яз +

( + ) 74 .

'

4

 

2

4

(

) ? 1 + ( + ) 7 2 + ( — ) ?3 + ( + ) 74

 

4

N

 

 

2 Я1Фц

 

( + ) 7 l + ( “ ) ? 2 + ( —) ?3 + ( + ) 74

1=0

"

4

4

Подставляя вычисленные коэффициенты Ьц в уравнение (XX.6), получим окончательный его вид с численными значениями этих коэффициентов в кодиро­ ванном виде. Для получения кодированных значений уровней факторов с целью подстановки их в уравнение (XX.6) используют формулу (XX.7).

Для проверки адекватности принятой модели используем критерий Фишера (F), определяемый выражением

S2

ад

F = 0 2

(XX. 11)

*{У)

 

21'

323

где — дисперсия адекватности;

N

N

2 Дд? — остаточная сумма квадратов.

N

2 Д'?(=Д91+ А<7 | + - . . + Д??+. . -+ Д9я, д ?3{= (?< —?i). i=l

<7j — опытные значения функции отклика; д,- — вычисленные значения'функции отклика по уравнению (XX.6); / — число степеней свободы для 5 |д при п — 4,*

/ = N —к — 1 + /д о п = 1;

* 4 , — дисперсия воспроизводимости,

.S,a — дисперсия (-того опыта; /; — число степеней свободы в (-том опыте,

ft = n i — 1,

л; — число параллельных опытов в (-том опыте.

Если полученное значение критерия Фишера меньше критерия табличного, принятую модель с соответствующей доверительной вероятностью можно счи­ тать адекватной.

В последнее время для обработки результатов спланированного экспери­ мента начали широко использовать ЭЦВМ.

РАЗДЕЛ IV

ПАРАМЕТРЫ ВЗРЫВНЫХ РАБОТ И ТЕХНИКА ИХ ПРОВЕДЕНИЯ

Г л а в а XXI

ЗАРЯДЫ ВЗРЫВЧАТЫХ ВЕЩЕСТВ И ПРИНЦИПЫ ИХ РАСЧЕТА

§ 82. Общие сведения о зарядах

З а р я д о м н а з ы в а ю т о п р е д е л е н н о е к о л и ч е ­ с т в о ВВ, п о д г о т о в л е н н о е к в з р ы в у . Заряды раз­ личаются по форме — сосредоточенные, удлиненные, плоские, ку­ мулятивные, Т-образные и т. п.; по расположению — внутренние (замкнутые) и наружные (накладные).

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

С о с р е д о т о ч е н н ы м и называют заряды, форма которых близка к кубической или шарообразной, а также заряды, макси­ мальный линейный размер которых превышает минимальный не более чем в 3 раза, в противном случае такие заряды относятся к удли­ ненным, которые бывают сплошными и рассредоточенными.

Р а с с р е д о т о ч е н н ы м и называют заряды, размещенные в одной зарядной камере, но отделенные друг от друга либо инертным

материалом, либо

воздушными

промежутками.

З а р я д н о й

к а м е р о й

называют выбуренный объем шпура

или скважины в породе или часть выработки, в котором размещается заряд ВВ.

При взрыве сосредоточенного заряда вблизи свободной поверх­ ности в грунте образуется воронка выброса, которую называют еще и конусом выброса, или горном (рис. 151). Кратчайшее расстояние от центра заряда до свободной поверхности W называют л и н и е й н а и м е н ь ш е г о с о п р о т и в л е н и я (л. н. с.). Это весьма важный параметр взрывной отбойки, входящий во многие инженер­ ные расчеты, иногда его называют г л у б и н о й з а л о ж е н и я з а р я д а . Расстояние г от края воронки до оси называют р а д и у ­ с о м в о р о н к и в ы б р о с а . Образующую конуса R называют р а д и у с о м р а з р у ш е н и я . Как видно из приведенной на рис. 151 схемы, явление выброса наблюдается только при условии, когда R > W. Во взрывном деле любые воронки взрыва принято характеризовать отношением r/W, которое называют показателем

325

Источник: https://tut-files.ru/previewfile/155599