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

1903

.pdf
Скачиваний:
0
Добавлен:
16.06.2024
Размер:
3.59 Mб
Скачать

оценке среднего значения в долях среднего значения определяемой характеристики; z1-α/2 – квантиль уровня; Р = 1-α/2 – нормированной нормально

распределённой случайной величины z x a ; Р = 1-α/2 – статистическая

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

Неизвестная величина генерального коэффициента вариации γ при определении объёма выборки заменяется выборочным коэффициентом υ, полученным на основании априорной информации по аналогичным материалам и элементам конструкций. При отсутствии аналога значением выборочного коэффициента вариации задаются с последующим его уточнением в процессе эксперимента. В этом случае объём испытаний корректируется методом подбора в соответствии с уточненным значением коэффициента вариации по формуле:

n

2

t2

,

(3.31)

a2

 

,k

 

 

где t ,k – значение квантили статистики t уровня Р = 1 – α/2 для числа

степеней свободы k = n – 1 [575].

Значением вероятности Р = 1 – α/2 при использовании формул (3.24) ÷ (3.26) задаются. Обычно принимается α = 0,1 или α = 0,05, реже α = 0,01. Формулы (3.24) ÷ (3.26) получены из условия, что вероятность попадания оценки математического ожидания в интервал a(1±∆a) или α±δaσ будет Р = 1 – α. Величина максимальных ошибок ∆a и δa определяет точность оценки среднего значения характеристики механических свойств строительных материалов. При низкой точности величину ∆a принимают равной коэффициенту вариации определяемой характеристики. В этом случае δa =1. При средней точности ∆а принимают равной (0,4 ÷ 0,5)γ и δa= 0,4÷0,5

при высокой точности – ∆a = (0,2÷0,3)γ и δa = 0,2÷0,3.

Если целью планируемых испытаний является оценка среднего квадратического отклонения характеристики механических свойств, то объём выборки определяется методом подбора [575]:

 

 

 

2

 

x2

 

1

 

 

 

/2

,

(3.32)

 

 

 

 

 

 

 

x2

 

 

 

 

 

0,5

 

 

где ∆σ – максимальная относительная ошибка (допуск) при оценке среднего квадратического отклонения случайной величины Х при нормальном

231

или γ = lgХ при логарифмически нормальном законе распределения; x2/2 и x0,52 – квантили уровня Р = α/2 и Р = 0,5.

Величина α принимается равной 0,1, или – 0,05. Значение ошибки ∆σ выбирается в зависимости от требований к точности оценки среднего квадратического отклонения характеристик механических свойств. При

низкой точности 0.4 0.5 при средней –

0.25 0.35 и при

высокой – 0.1 0.2 .

 

При n ≥15 объём выборки определяется по приближённой формулой вместо (2.77) [575]:

z2

n 1,5 1 /2 . (3.33) 22

При оценке квантильных значений характеристик механических свойств строительных материалов и конструкций минимально необходимый объём испытаний определяется исходя из заданной точности и статистической надёжности оценки квантили. При нормальном распределении величины Х минимально необходимый объём серии образцов определяется методом подбора из уравнений [575]

 

 

 

 

z

 

t

 

n

1;

 

 

 

1

 

 

 

для Р 0,5 ,

(3.34)

 

 

 

 

P

 

 

 

 

P

 

 

 

 

 

1

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

n 1;

 

1

z

P

 

P

для Р 0,5 ,

(3.35)

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

квантиль

уровня

β нецентрального

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

t

n 1;

 

Стьюдента с n-1 степенями свободы и параметром нецентральности. Значением доверительной вероятности β=1-α при использовании формул (3.29) и (3.30) задаются. Обычно принимает β=0,9 реже 0,95. Величина максимальной относительной ошибки оценки квантили принимает δР=0,2÷0,3 при высоких требованиях к точности, δР=0,4÷0,6 при средних требованиях и δР = 0,8÷1,0 при низкой точности.

Необходимое количество измерений из-за условия пренебрежения случайной погрешностью, по сравнению с приборной погрешностьюx cos e определяется по формуле [575]:

N 4n

n2

(3.36)

devige2 ,

где n – число измерений, позволяющих пренебречь коэффициентом Стьюдента; 2n – дисперсия случайной погрешности из сделанных изме-

рений; 2devige – дисперсия приборной (систематической) погрешности.

232

Необходимое число количества экспериментальных образцов [35], испытывающихся на выносливость, обеспечивающих требуемую точность

и надёжность параметра у с учётом y , определяется по формуле:

 

t

 

2

(3.37)

n

 

 

,

 

y

 

 

 

 

 

 

 

где σ – предполагаемая дисперсия, принимается по предварительной выборке или на основании ранее проведённых испытаний; t – квантиль из

табл.3.6; y

– желаемая точность результатов испытаний: y

 

y y

 

 

 

y

 

.

 

 

 

 

Исходными расчетными данными при испытании образцов на выносливость являются:

отклонение логарифма долговечности от своего среднего арифметического значения lg N 0, 08;

значение t 1,96 при надёжности 0,9 (табл.3.6);

дисперсия 0, 20 , принимаемая по данным ранее проведённых

испытаний.

Необходимое количество образцов составляет:

 

0, 20

1,96

 

2

(3.38)

n

0,

08

 

24 .

 

 

 

 

 

 

 

 

Таблица 3 . 6

 

Таблица квантилей для нахождения объёма выборки

 

 

 

 

 

α

 

t

α

t

0,00

 

0,00

0,83

1,35

0,05

 

0,08

0,84

1,39

0,10

 

0,13

0,85

1,44

0,20

 

0,25

0,86

1,47

0,30

 

0,38

0,87

1,51

0,40

 

0,52

0,88

1,55

0,50

 

0,68

0,89

1,60

0,55

 

0,75

0,90

1,64

0,60

 

0,83

0,91

1,69

0,65

 

0,93

0,92

1,75

0,70

 

1,04

0,93

1,81

0,72

 

1,08

0,94

1,88

0,74

 

1,13

0,95

1,96

0,76

 

1,18

0,96

2,05

0,78

 

1,23

0,97

2,17

0,80

 

1,28

0,98

2,32

0,81

 

1,30

0,99

2,58

0,82

 

1,32

0,999

3,29

Аналогичный пример определения необходимого количества опытов по измерению пульсаций давления газа в компрессорной установке по формуле (3.33) приведён в работе [194]. Установлено, что уменьшение

233

доверительной вероятности Р с 0,99 до 0,95 приводит к уменьшению числа опытов в 1,73 раза.

Для минимизации необходимого числа соответствующих измерений на строительных конструкциях используется способ статистической проверки гипотез в виде последовательного анализа получаемых чисел наблюдений. При этом необходимое число наблюдений не фиксируется заранее, а определяется в процессе самой проверки, так как число наблюдений при последовательном анализе есть величина случайная, чем при способах, в которых число наблюдений фиксировано заранее. К последовательной схеме наблюдений обращаются в ситуациях, когда каждое наблюдение является дорогостоящим или труднодоступным. При последовательной процедуре из-за ограниченности выборки, полученной из генеральной выборки, вместо одного порога как в классической процедуре проведения экспериментальных исследований задаются двумя порогами статистической проверки опытных результатов. Они включают проверку основной, или нулевой гипотезы, по которой возможны ошибки как в ту, так и в другую сторону по сравнению с альтернативной гипотезой, конкурирующей с ней [120].

3.4. Линейные корреляции

Результаты эксперимента, выступающие в роли переменных и параметров некоторой функциональной зависимости, теоретически получаемой в рамках математической модели, менее трудоёмки в проверке, если они представлены в виде линейной зависимости [670]:

y ax b ,

(3.39)

где x, y – измеряемые величины; a, b – параметры зависимости.

В отдельных случаях при нелинейном виде зависимости она может преобразоваться в линейную (табл. 3.7).

Таблица 3 . 7 Виды преобразованных нелинейных зависимостей в линейную

Вид

Получаемая линейная

y

x

a

b

нелинейной

зависимость

 

 

 

 

v k uz

ln v z ln u ln k

ln v

ln u

z

ln k

v k ezu

ln v zu ln k

ln v

u

z

ln k

 

 

 

 

 

 

 

 

v k ez/u

ln v zu 1 ln k

ln v

u 1

z

ln k

v

u

v 1 ku 1 z

v 1

u 1

k

z

k zu

 

 

 

 

 

 

234

Принятие гипотезы о линейности экспериментальных данных позволяет подтвердить влияние выбранного количества парных измерений n на ве-

рификацию(проверкусправедливости) используемогомодельногоописания. Метод наименьших квадратов является одним из наиболее распространенных приемов статистической обработки экспериментальных данных, относящихся к различным функциональным зависимостям физических величин, в том числе и к линейной зависимости, позволяющим получить достоверные оценки ее параметров a и b, с оценкой их погрешности [670]. Он позволяет определить наилучшую линейную аппроксимацию данных, поскольку при данных значениях параметров достигается минимум вели-

чины «отклонения» прямой от экспериментальных данных:

 

 

n

 

 

 

 

 

 

2

 

 

 

 

 

ax1

b .

(3.40)

s i 1 yi

Выражения для a:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

a k j y j ,

 

 

 

 

 

(3.41)

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

где k j – коэффициенты, k j

 

 

xl

x

 

;

y j

 

– взаимно

независимые

 

n

 

 

 

 

 

 

 

 

xj

x 2

 

 

 

 

 

 

величины.

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оценка для стандартного отклонения параметра a [185]:

 

a

 

1

 

 

y2

 

y 2

 

 

a2 .

(3.42)

 

n 2

 

x2

x 2

 

 

 

 

 

 

 

 

 

 

 

Выражение для b:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

1

 

n

yi a

1

n

 

x i .

(3.43)

 

n

 

i 1

 

 

 

n

i 1

 

 

Стандартное отклонение нормально распределённого параметра b:

b

2

x 2 2

a

x 2 .

 

(3.44)

b

n

a

 

 

 

 

 

 

 

 

 

При переходе от величин a

и b к погрешностям

a

и b они

умножаются на коэффициент Стьюдента:

 

 

 

 

a t , n 1 a ,

 

 

 

(3.45)

b t , n 1 b ,

 

 

 

(3.46)

гдеα– уровеньзначимости; n – количествопарныхизмерений(см. табл. 3.5).

235

Значение n 1 берется в связи с тем, что в методе наименьших квадратов из экспериментальных данных находят не одну величину, а две – a и b. Связь между ними уменьшает количество независимых случайных переменных, складывающихся в распределение Стьюдента [670].

В методе наименьших квадратов остаточная сумма квадратов определяется как:

n

yi ax b

2

 

 

 

(3.47)

i

i 1

 

 

где в качестве a и b использованы их экспериментальные оценки;

i

дисперсия отдельного измерения. В статистике обосновывается, что величина подчиняется распределению («хи-квадрат») 2 (3.5).

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

задающего интервал от 0 до x2 n, в который величина попадает, если гипотеза справедлива. Величина x2 n, приведены в табл. 3.8.

 

Границы интервалов x2 n,

 

 

Таблица 3 . 8

 

для уровня значимости

 

 

( n – количество парных измерений)

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

0,75

0,95

 

0,99

 

0,999

4

 

1,3

3,8

 

6,6

 

10,8

5

 

2,8

6,0

 

9,2

 

13,8

6

 

4,1

7,8

 

11,3

 

16,3

7

 

5,4

9,5

 

13,3

 

18,5

8

 

6,6

11,0

 

15,1

 

20,5

9

 

7,8

12,6

 

16,8

 

22,5

10

 

9,0

14,1

 

18,5

 

24,3

11

 

10,2

15,5

 

20,1

 

26,1

12

 

11,4

16,9

 

21,7

 

27,9

13

 

12,6

18,3

 

23,2

 

29,6

14

 

13,7

19,7

 

24,7

 

31,3

15

 

14,9

21,0

 

26,2

 

32.9

16

 

16,0

22,4

 

27,7

 

34,6

17

 

17,1

23,7

 

29,1

 

36,1

18

 

18,3

25,0

 

30,6

 

37,7

19

 

19,4

26,3

 

32,0

 

39,3

20

 

20,5

27,6

 

33,4

 

40,8

30

 

31,5

40,1

 

47,0

 

55,5

50

 

56

68

 

76

 

87

100

 

109

124

 

136

 

149

236

Если неравенство rs2 x2 n, (s – минимальная величина «отклоне-

ния» прямой от экспериментальных данных) не выполняется, то гипотеза о линейности отвергается.

В случае отсутствия сведений о погрешности отдельных измерений для проверки гипотезы линейности используется величина:

 

R

n 2,

(3.48)

1 R2

 

 

 

где n – число измерений; R – коэффициент корреляции Пирсона, безразмерный индекс в интервале от -1,0 до 1,0 включительно, который отражает степень линейной зависимости между двумя множествами данных:

R

cov xy

 

in 1 xi

x yi

y

 

.

(3.49)

x y

n

x1

x

2

n

y1 y

2

 

 

 

i 1

 

i 1

 

 

 

В статистике обосновывается, что величина λ подчиняется –распре- делению с n 2 степенями свободы. При выполнении условия x2 n, гипотеза о линейной зависимости отвергается [670].

3.5. Использование корреляционного и вероятностного анализа при статистической обработки результатов испытаний на выносливость строительных конструкций

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

Долговечность при испытании на выносливость определяется числом циклов с соответствующими уровнями напряжения в цикле нагрузки, выдержанных образцами до разрушения. В то же время даже при испытании более 40 штук металлических образцов на соответствующих уровнях напряжений рассеяние по долговечности, то есть отношение наибольшего числа циклов к наименьшему, может составлять 50 и более раз [35]. Причиной количественного разброса опытных данных является статистическая природа процесса усталостного разрушения материала. Объяснением является структурная неоднородностью металла: различные размеры, форма и ориентировка зёрен; наличие различных фаз и включений; случай-

237

ные изменения в микрогеометрии кристаллической решётки и структуре поверхностного слоя. Поэтому при испытаниях на переменное напряжение рекомендуется испытывать не менее 10 одинаковых прямых моделей инженерно-строительных конструкций как из металла, так и из другого материала [35].

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

Долговечность образцов на выносливость при напряжении σ устанавливается при определении функции распределения вероятности разрушения P N , в зависимости от числе циклов N. Целью статистической

обработки результатов испытания на выносливость определённого количества образцов (12-20 штук) при одном уровне напряжения является обоснование параметров полученной функции распределения [35].

В качестве функции распределения вероятности разрушения образца до N циклов наиболее часто выбирается логарифмически нормальный закон распределения нормированной случайной величины U p , называе-

мую квантилем распределения [537]:

U p

X a

,

(3.50)

 

 

 

 

где a – математическое ожидание случайной величины Х; σ – среднее квадратичное отклонение случайной величины Х.

Характер рассеяния долговечности при испытаниях на выносливость, соответствующий логарифмическому нормальному закону распределения, был доказан с помощью экспериментальных исследований на большом количестве образцов, составляющим 400-500 штук [537].

Методика проведения линейного корреляционного анализа результатов испытаний на долговечность экспериментальных металлических конструкций ограниченного количества до 12 штук включает следующие этапы расчёта [35]:

1)результаты испытаний в виде напряжения i , при которых образцы разрушаются, логарифмы этих напряжений xi lg i , а также значения долговечности N и их логарифмы yi lg Ni в порядке убывания должны принимать табличную форму ;

2)средние арифметические значения логарифмов напряжений X и логарифмов долговечности y в зависимости от количества образцов n

вычисляются по формулам:

y

n

yi

 

(3.51)

i 1

 

,

n

 

 

 

 

 

238

x

 

n

xi

 

 

 

 

(3.52)

 

i 1

 

;

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

3) средние квадратичные отклонения напряжений mx

и долговечности

my имеют зависимости:

 

 

 

 

 

 

 

 

 

mx

in 1 xi x 2

,

(3.53)

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

n

 

 

2

 

 

 

my

 

i 1 yi y

 

 

;

(3.54)

 

 

n 1

 

 

 

 

 

 

 

 

 

4) ковариация напряжений и долговечности mxy определяется как сумма произведений частных отклонений напряжений и долговечности:

 

mxy

in 1 xi x yi y

 

;

 

 

n 1

 

 

 

 

 

 

5)

коэффициент корреляции r назначается в виде:

 

 

r

mxy

;

 

 

 

 

 

 

 

 

 

 

mx my

 

y

6)

мера индивидуального рассеяния чисел

 

формуле:

y my

1 r2 ;

(3.55)

(3.56)

рассчитывается по

(3.57)

7) на основе меры индивидуального рассеяния определяется математическое выражение для семейства кривых выносливости при малых значениях вероятности разрушения 5 %, 2,3 % и 1 % [35]. При заданной вероятности разрушения для описания кривых выносливости уравнение усталостного разрушения металла имеет вид [99]:

y BX A K y ,

(3.58)

где А и В – параметры уравнения линии регрессии:

 

A y r

my

x; B r

my

;

 

mx

 

 

 

 

mx

 

K – коэффициент, учитывающий,

как вероятность разрушения,

так и

количество испытаний образцов, принимается из табл. 3.9 распределения Стьюдента [35];

8) при заданной вероятности в логарифмических координатах строится линия регрессии и кривая выносливости;

239

Таблица 3 . 9

Таблица распределения Стьюдента

Р, %

 

 

 

 

 

 

 

 

Количество образцов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

99

 

4,50

3,74

3,36

3,14

3,0

2,90

2,82

2,76

2,72

2,68

2,65

2,62

2,60

2,58

2,56

2,55

2,54

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

97,7

 

3,30

2,86

2,64

2,51

2,43

2,36

2,31

2,27

2,24

2,22

2,00

2,19

2,18

2,17

2,16

2,15

2,14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

95

 

2,36

2,14

2,02

1,94

1,90

1,86

1,84

1,82

1,80

1,78

1,77

1,76

1,75

1,74

1,74

1,73

1,73

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

90

 

1,64

1,58

1,47

1,43

1/41

1,40

1,39

1,38

1,38

1,37

1,36

1,35

1,35

1,34

1,34

1,33

1,33

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

240

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]