Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Баженов, Ю. М. Перспективы применения математических методов в технологии сборного железобетона

.pdf
Скачиваний:
11
Добавлен:
19.10.2023
Размер:
6.02 Mб
Скачать

voа

к

о

о

С*.

>-

О

*

©•

о., з

"ач

I н

+1

ö

+

іч

в'і

» Ч

Ч

+1

ö

+

I ы

I ч

+

 

СЧСЧt)

I ч

+

CA

1ч~*

О

 

І ч7

 

 

I Ч

■ъ

I *

Ч

+1

+

сч —

Ü

O'«<м

I ч

сч •—

I Ч

+

 

I *

I *

+

+

I ч

I *г

Ч

+

а

Ччw

+

\X I ч

++

I £'

I ч

+

+

+с;

•ч

Ч

100

Обращает внимание последний результат: относительная ошиб­ ка определения площади цилиндра вдвое больше относительной

ошибки

измерения его диаметра.

Б)

Рассчитать ошибки сг{У} и б {У} определения на цилиндри­

ческом

образце предела прочности R = P :F . Воспользуемся форму­

лой (3)

в табл. Ш.8:

а)

а W = у

4P )

Р2о2 {F)

ц Р 2

F*

б)

<y{R} = [ - j J V S ‘ {P} + 6^F}-,

в) б{Л }=

У

m

+

г)

6 {R} =

|/Л = {Р }~

6'{F(.

Зависимость

(III.39)

должна

использоваться экспе­

риментатором

 

не только

для расчета ст{У} по а{х'г}

ит. п., но и для определения допустимых ошибок а{Хі}

иб{л'г} по заданной величине а {У} и 6{У}. Это один из обязательных этапов управления процессом эксперимен­ та и контролем качества.

На основании анализа закона распределения ошибок

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

Пример III.14. Проанализируем по данным примера III.13 ус­ ловия получения заданной ошибки a{R} и 6{Р}.

Если задано 6{Р }=2% , то

б2 {R} = б2 {Р} + б2 (Р) = б2 (Р) + 4б2 {d) = 0,0004.

Поскольку 62{Р) задано конструкцией испытательной машины (весьма часто 6{Р} = 1%), то регулировать можно 62{d}:

0,0004 — 0,0001

б2 {d} = ---------— 1-------

= 0,000075 или б {d} = 0,865%.

Если такая ошибка инструментально недостижима, то следует увеличить число параллельных измерений диаметра образца. Пред­ положим, что возможна бВозм{гі}=2%, тогда, воспользовавшись (ІІІ.35), найдем:

62 W

101

 

■_ Сзм_№

_

0,0004

_

 

 

 

п ~

б2 (d)

~

0,000075

'

 

Если при испытании цилиндра гі=150 мм из бетона марки Р =

= 300 кгс/см2 необходимо обеспечить сг{/?} =

10

кгс/см2 (при условии

б {Я} = 1%

II Cf{d} = 2 мм), то:

 

 

 

 

 

о2 {R}

= (Р : Р)= (б2

{Р} +

6= {Р}) = (Я )2(б2

{Р) + 4

,

откуда

 

 

 

 

 

 

 

 

Отреб {^} — 0,5d V

 

 

 

 

= 75 V (10: ЗОО)2 — 0,001- =-= 2,385 м м .

Поскольку 0треб{с/} >o{rf} = 2 мм, то можно считать метод из­

мерения диаметра цилиндра удовлетворительным по точности.

Пример 111.15. Проанализируем относительную ошибку a{R} определения прочности бетона па сжатие R по известной формуле в предположении независимости коэффициентов А и а:

 

 

1 п

 

\

 

 

 

 

 

 

Р =

Л / Ц —

- M j =

- ^ Щ

+

аВ).

(111.40)

Пользуясь (111.39) и табл. 111.8,

III.9, находим с учетом того,

что Л и а тоже определяются с ошибками:

 

 

 

б2 (Р) =

б2 (Л)

+ б2 {Рц}

+

б2 (S)

+ б2 -I-аВ] = б2 (Л) +

 

 

Ц 2б2 (//} '+

а2 В 2(б2 (а)

+ б 2 ( 5 ))

+ 62 {РЦ) + 6 2 (ß)

 

 

 

-faß)2

 

 

 

 

 

 

 

 

Оценим

роль _62{Л}_ и

 

б2{с},

принимая

6{ß} =6{Z(} = 1 %,

б{Рц}=5% при ЩВ = 2, я= —0,5:

 

 

 

 

 

 

б2 {/?)

= б2 {Л} +

а2 В 2

б2 (а) +

 

 

—_

_ _

 

 

 

 

 

Ш+аВ)2

 

 

 

б2{Рц

 

а2 В 2

б2 {ß} +

 

 

Щ + а В ) 2

 

 

 

 

 

 

 

 

W-

62{Z/}

=

б2 (Л)

+

_

_

_

б2 (а) +

Ш + а В ) 2

 

 

 

 

Ц /В + а

 

+ J 6 2 { Р ц } +

1 +

Ц /В + а

 

 

 

\ Ц / В + а

 

 

 

 

 

102

= б« (Л) + Y

63 {«)

+

° - 053 +

 

° '012 +

 

0,012

= б2

{/1) + 0,1162

{а } -1- 0,002789.

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

а) даже если коэффициенты А

и а

определены с

пулевыми

ошибками, то

б {R} не

может

быть

менее 5%

(1^0,002789-100%);

б) ошибка

в определении

коэффициента /1,

который

характери­

зует в технологии бетона влияние качества

заполнителя, примерно

в три раза более существенна,

чем ошибка б {а}.

В экспериментальной

практике

истинные значения

о{хі} обычно неизвестны и заменены /е оценками s {a:,}, каждая из которых определена с числом степеней сво­ боды \і = Пі—1. В этом случае оценка среднеквадратич­ ной ошибки косвенного измерения s { y } , рассчитанная по (III.39), имеет число степеней свободы /{У}, опреде­ ляемое (предложение Уэлча) следующим образом:

 

k

 

k

f { Y } =

, (ИІ.41)

где s{y} принимается по табл. III.8 (с заменой символа генераль­ ного среднеквадратичного отклонения на символ оценки).

Для

абсолютного Д — R2— R і

[формула

(1)

в табл.

III.9] и относительного p= Rz'-Ri

[формула

(3)

в табл. III.9] показателей изменения некоторого пара­ метра R (такие задачи часто встречаются в технологии) число степеней свободы определяется как:

103

Пример ІИ.16. Определить число

степенен

свободы для абсо­

лютного прироста

прочности бетона

с химической добавкой, если

при испытании его

получено ==250

но /ц = 6

образцов, а при ис­

пытании эталонного бетона но /о = 1 2 образцов было

===100.

По (111.42) находим

 

 

 

/ { Д) = (400 + 250)=: —

400'-+ —

2502

=

 

11

5

1

 

 

= 422 500:27 045 = 15,02.

 

 

При использовании таблиц t- и /’"-критериев соответ­

ствующие числовые значения

находятся по

дробному

/{У} интерполяцией.

 

 

 

III.6. Принципы выбора кривых распределения по опытным данным в рецептурно-технологических зада­ чах. Числовые значения критериев качества бетона и чи­ словые значения технологических факторов, определя­ ющих это качество, можно рассматривать как случай­ ные величины (см. гл. II). Чтобы задать такую величи­ ну, нужно указать не только ее возможные значения, но и вероятности, при которых эти значения достигаются.

При этом

максимально полной является информация

о функции

распределения вероятностей /?{x}, а мини­

мальной — таблица зафиксированных в ходе исследова­ ния выборочных значений случайной величины.

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

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

второй — определение точечных числовых оценок рас­ пределения: среднего х\ дисперсии sz{x}; коэффициен­ тов вариации ѵ{х}, асимметрии А*, эксцесса Е* и т. д.; третий — определение интервальных оценок, кото­ рые дают уверенность (с определенным риском а) в том, что параметр генеральной совокупности лежит внутри

интервала, связанного с точечной оценкой;

104

четвертый-— аппроксимация выборочного распреде­ ления теоретическим законом распределения или подхо­ дящей эмпирической кривой (с проверкой адекватности принятого решения).

До последнего времени большинство технологиче­ ских работ ограничивается двумя первыми этапами, при этом используются лишь элементарные алгебраические

операции при вычислении х и s2{%}, а не вероятностностатистический метод познания, который начинается на третьем этапе.

Технолог, определивший интервальные оценки и тем более правильно аппроксимировавший выборочное рас­ пределение, может при том же объеме исходных данных ответить на ряд важнейших технологических вопросов, в частности:

об области существования случайной величины; о точности точечных оценок параметров распреде­

ления;

овыборе необходимого числа наблюдений для полу­ чения заданной точности параметров выборки;

осправедливости тех или иных нуль-гипотез (напри­ мер, о равенстве двух выборочных средних) и т. д.

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

но-технологических ситуаций и т. п.

Выбор типа закона распределения — задача не­ простая. Наиболее популярен в практике нормальный закон (см. гл. II), физическая модель которого основа­ на на предположении, что случайная величина представ­ ляет собой результат «малых» воздействий (примерно одинаковых по величине и равновероятных по знаку) большого числа взаимоиезависимых факторов. Однако многие явления не могут быть описаны таким законом по своей физической сущности: размеры зерен при дроб­ лении, число дефектов в образце и др. В этом случае следует применять иные распределения: непрерывные (экспоненциальное, Вейбулла и др.) или дискретные (биноминальное, Паскаля и др.) [41, 79, 75]. Целесооб­ разно анализировать соответствие выбранного распреде­ ления физической модели явления, однако во многих

105

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

Информация о закономерностях распределения слу­ чайных величин в виде функций пли числовых характе­ ристик необходима технологу при решении всех задач

Рис. III.1. Диаграмма для оп­ ределения типа кривых распре­ деления в координатах ßi и ß2

— нормальное; f ) — равномерное;

U

— экспоненциальное; ... — логнор­

мальное; Л — критическая область

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

Для выбора кривой распределения по известным чис­ ловым оценкам асимметрии А* и эксцесса Е* удобно пользоваться методом, основанным на анализе оценок

коэффициентов формы кривой ßi и ß2:

 

ß; = {АУ

(III.44)

и

 

ß.> = 3 -f E*.

(Ш.45)

На рис. II1.1 показаны в координатах ßi и ß2 обла­ сти существования шести различных теоретических рас-

106

Наздйние иЬид (рцнкции распределения
НормальHOe„ /h/?’i

пределеиии, а в табл. III.10 даны сведения о них; в за­ висимости от соотношения между ßi п р2 технолог нано­ сит точку С и выбирает подходящий закон распреде­ ления [75].

Таблица Ш10

Некоторые

,

,

 

 

непрерывные распределения

 

 

Плотность распредели

м|х| и jflxl

ßi и ßi

 

ния и ею параметры

 

 

 

11/

^

*

i

 

 

 

 

 

Iff

•I

 

лог-ритмически нормальное

 

ад

f\0)0

 

 

 

 

Z J

1 V

'

 

 

 

 

 

V

 

q-м

 

 

f/ ^ \ _

бЧ.о

 

 

7

1 2

3 h l

 

 

Ра6но*ерное

 

 

 

 

У|'і)

 

%

 

ъ

 

 

 

 

 

1

 

/Of

 

_I-----------

 

 

>----

 

-распределение

 

 

V(r)

 

 

 

 

 

 

 

1

 

-

 

r.;

 

 

 

Д 5

 

_ Л3 f

 

/

 

 

------

 

 

0

1 2

3

h

l

 

muaß

Ч’ІХІ

 

 

 

 

qH

 

2

Ш

$

\

 

Р»

 

I

3

 

 

0

 

 

 

e:

munßoi

2 10-2

 

<3

?<<

 

'

 

x N v H

c:

t<>

 

 

a

 

 

 

----- X

 

 

«3.

munßlul

0

 

 

I

«3

'

\q>0ß

 

І'Ч

 

'

\tVn-o.s\

 

 

— ■---------П

 

 

 

 

 

Зкспоненииольное

 

 

 

rWiv

 

 

 

 

 

 

 

1----^

 

------I

 

- оо

<X<■»оо

261

м1х\=п

/5, ‘ 0

 

параметры.

 

й\х\-б2

ß ,-3

 

-оо<у< .ос (

Q >0

 

 

1

[Щі-дІ

М(х)--е*‘?

X

 

э

 

біШ

151

 

 

 

ХіО

 

 

ЛМ=е2,’б(ы-і)

Э

'

3ZO

параметры

 

 

'

-

â

2

 

<.П<-»-оо

са-ехр[бг}

.3 ~э >-

з

- с о

7

к

 

0>О

 

 

 

«о, 'О.

 

 

о. 1 <ч0UІ>2,

м ( х ) ^

А -»

 

а^\-(Ц^Аг

 

параметры о ап

ß, *

>-В

 

 

L0 W

 

 

 

 

 

П?)

 

 

 

м1x1 =в А

 

 

 

 

 

 

,х<0 л[х}-? /1г к\(ѵп

 

0

 

 

 

параметры \>0,%>0

 

 

 

 

 

\ lM L x !'b-xr

 

Т.

15

 

m m

 

'

 

 

 

npuOmi

ф)-- —

"?■

 

 

0,

1<0иХ>1

 

O'

 

 

 

 

Дг

 

+

 

 

 

 

 

C“ СМ

 

 

 

 

 

. Со—1>ч;

 

параметры

«/и__ Ц__

‘-Т-' зі, Кг

4 - ^

=>

 

q >0

 

 

 

^

Т ^

 

 

 

 

'vüT

 

^=*

 

Г 1

 

 

 

1-

 

ГѴ,

 

 

 

 

 

 

 

 

 

 

 

ее^

ссД

Ae“ . xio

 

ФМ"

ßr h

 

0

. X<0

 

Л1Х)=А’г

ßr9

 

 

Параметр:А>0

 

 

 

 

 

 

 

Пример III.17. При анализе 134 проб песка (для приготовления нзпестково-песчаных автоклавных бетонов) на содержание глины были получены следующие числовые оценки (при 15 группировках

107

с интервалом

через 1%):

*= 4,56%;

s= {.ѵ} =3,28%; /1* = 1,095;

£* = 0,52.

 

 

 

Определим ßj и ßj н

нанесем на график (рис. 111.1) точку С.

ß* =

(1,095)2 =

1,199; ßj =

3 - |- 0,520 = 3,52.

Можно выбрать для описания распределения содержания гли­ ны в песке распределение типа ß (/). Если бы в данной задаче бы­ ла принята кривая нормального распределения, то этот вывод имел

достоверность менее 0,05% (%]j,=63,2 при /= 1 1 —3 = 8, так как 5 интервалов пришлось бы объединить в один).

Истинные значения коэффициентов ßi н ß2 находятся в некотором доверительном интервале возле оценок ßj и ßj, следовательно, С нужно представлять себе не

как точку, а как область, очерченную эллипсом, который (особенно в зонах ßj=0-=-l и ßj =0-М , а также в зонах

ßj = 2,5л-3,5 II ßj =7-4-8) может перекрывать несколько

линии на рис. III.1. В результате нужно проверить не­ сколько теоретических распределений и выбрать из них наилучшее как по %2-критерию, так и по физическому смыслу механизма явления.

Пример III.18. Определить параметры ß-распределепия по дан­ ным примера III.17. Из табл. ШЛО находим, что любая ß-форма имеет два параметра распределения ц и у. оценки которых т|* н yjy

можно [75] найти по (III.47) и (111.48) в зависимости от оценок х

и s{*}. Однако следует учесть то обстоятельство, что ß-pacnpe-

делеиие существует в интервале О ^ х '^ І ,

поэтому оценки х и s{*}

необходимо представить в масштабе, взяв

их отношение к макси­

мальному значению .ѵМакс =

Х *= 15,5% (это значение соответствует

примеру 111.17):

 

 

 

 

*0 = Л-: А,* =

4,56: 15,5 =

0,294;

|

(111.46)

Sq= s {*) -Л* = 3,28 : 15,5 = 0,212;

J

 

1Г—----[•ѵ’и (1~ -ѵо) — sii] = sö1

1 — 0,294

 

 

2,569;

(III.47)

=

[0,294(1 - 0 ,2 9 4 ) - 0 ,2 1 2 2 ] =

У*

*0i f _

0,294-2,569

1,071.

(III.48)

l - * 0 “

1 - 0 ,2 9 4

 

 

 

 

Функция распределения имеет при г)>1 и у > 1

(см. табл. ШЛО)

вид асимметричной

колоколообразной кривой

с

максимумом при

* { м а к с } I

 

1,071 — 1

 

 

 

,_ É z

 

15,5 = 0,67

(III.49)

y*-(-T|*-f2 X* =

,071+2,569+2

108

II

описывается

следующим

уравнением

(с учетом

Я*, заменяя т)

па

11*

и у на у*):

 

 

 

 

 

 

 

Г (У + V) ( * jV- 1Л ___ £_уі_1

 

 

1

Г (3,640)

(^ - У '07І- Ѵ і -

2,569—1

 

 

 

 

J L \

 

 

15,5 ' Г (1 ,071) Г (2,569) \15,5/

\

15,5/

 

 

 

= 0,18 '

X N0,071

X у . 569

 

 

 

 

15,5/

(111.50)

 

 

 

 

 

 

где Г — гамма-функция.

 

 

 

 

 

Проверка

по х2= П ,6 9

при (= 1 1 —4— 1 = 6

(параметры г)*,

у*,

Хо = 0 и Л*=15,5) показывает, что Р{х2}> 5% и гипотеза адекват­ ности не отклоняется. По физическому смыслу ß-распределение можно признать удовлетворительным, так как одно из его интер­ претаций [75]: распределение доли совокупности (содержание гли­ ны в %), заключенное между наименьшим и наибольшим значения­ ми выборки.

Проверка

иных распределений

по рис.

III.3

дала

в примере III.17 следующие значения %2-критерия: лога­

рифмически

нормальное

[79] — %2 = 121,8

и

Р{%2} <

<0,05%;

 

распределение

Грама-Шарлье

 

79]

%Ф= 32,5 и

Р{х2} <0,05% ;

гамма-распределение

 

7 5 ] -

5Сф= 18,8 и Р{х2} <2,5% ; распределение Вейбулла

79] —

Х| = 17,9 и

Р{%2}<2,5% . Таким образом, несмотря на

несколько

лучшее приближение по

сравнению

с

нор­

мальным

 

распределением

 

 

 

 

(пример

II 1.17), пи одно

из

 

 

 

 

этих распределений не может конкурировать в данной зада­ че с ß-распределением (пример

III.18).

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

ются

кривые

Пирсона

[41].

 

 

 

Область возможных

значений

 

 

 

ßi и

ß2, описываемых такими

 

 

 

кривыми, несколько больше,

0

....

п

чем

в семействе

кривых на

рис.

ттт ,

 

г

 

пока-

Рис.

111.2.

Диаграмма

II 1.1. Однако

ОПЫТ

для

определения типа

зывает, что

расчеты

парамет-

кривых

Джонсона

109

Соседние файлы в папке книги из ГПНТБ