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

Математическая статистика в технологии машиностроения

..pdf
Скачиваний:
24
Добавлен:
15.11.2022
Размер:
8.19 Mб
Скачать

 

 

Таблица данных для

регрессионного анализа

 

 

 

 

 

 

(/)

 

 

 

 

 

 

 

•-»

 

с"*

 

 

*

 

 

 

с*

f 1

 

 

сГ

 

к*

 

V)

«Г

 

 

е*

 

 

1

2

3

4

5

 

6

7

8

 

9

? i

S2

nL

«1 — 1

(«г— О-5?

nL-Yi

Х1

 

х1— Х

 

 

 

 

 

 

 

 

 

 

F 2

s2

п2

п2 — 1

(«г- 1 ) ' S2

n2-Y2

*2

п2>х2

х2 — Х

S2

F т

S2

Пт

пт — 1

{Пт 0 ‘sm

Пт• Yт

хт

пт'хт

хт X

 

т

 

2(2) 2<з)

2(4)

2(®>

2 (6)

2(8)

 

 

 

 

 

 

 

 

 

 

СМ

 

 

С1

' l

 

 

 

С1

|>1°

я

 

|>г

 

0

 

3

I?4

 

L

 

 

 

I*

 

L

 

 

 

о

1*

 

1^

L

 

сГ

 

сГ

1

 

 

 

 

 

 

1^

\*

 

 

 

10

 

11

12

13

14

 

15

16

(X l - x y

«1 ( * 1 - * ) 2

- F)

Foi

F i - F o i

( Y i - Y o i V «1 ( у 1— Y oi)2

( * 2 - * ) 2 М * * - * ) 2 n2Y2(x2 — X)

~Y02

F 2 — К ос (F2- F 02)2«2 (F2-"F 02)2

(хт

(*m — X)2п т У т (хт

У От

-- К 0 Ш (Fт

F от)2

«т (Y щ

X)2

 

 

- X )

 

 

 

 

- F

0m)2

 

2

( и )

2

(>2)

 

 

 

 

2

(16)

Пользуясь итоговыми данными табл. 34, вычисление всех величин, необходимых для регрессионного анализа, произво­ дится по следующим формулам:

а

а

in

s\ =

S(2) _

 

in

’>

„2 .

E

(16>

 

b2~

m —2

! .

e '= £ ® .

m. ^-2 »

 

Ц ( 3) ’

2__

S

(5) + S

(16)

 

[ E ( 3 ) - m ] + (m — 2) ;

n2 ___

S2

o 2 ___ S2

а - Ц ( з ) ;

S d O *

(179)

(180)

(181)

Ч А С Т Ь В Т О Р А Я

ПРИМЕНЕНИЕ МАТЕМАТИЧЕСКОЙ СТАТИСТИКИ 0 ТЕХНОЛОГИИ МАШИНОСТРОЕНИЯ

Глава I

СТАТИСТИЧЕСКИЕ МЕТОДЫ В ТЕХНОЛОГИЧЕСКИХ ИССЛЕДОВАНИЯХ

При исследовании технологических процессов часто прихо­ дится сталкиваться с величинами переменного характера, которые

врезультате опытов могут принимать любые численные значения

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

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

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

Большей частью случайные величины, встречающиеся в тех­ нике, в том числе и в технологии машиностроения, подчиняются закону нормального распределения. Экспериментальным путем всегда можно проверить, подчиняется ли распределение интере­ сующей нас случайной величины закону нормального распреде­ ления. Повторяя опыт не менее 50 раз, можно рассматривать его как случайную большую выборку из некоторой генеральной сово­ купности. Проверив выборку на случайность и нормальность,

можно убедиться в справедливости предположения о

подчине­

нии

изучаемого признака нормальному распределению

или

по

эмпирической кривой

распределения установить близость

его

к какому-либо другому

закону распределения случайных вели-

8

И. С. Солонин

 

 

113

чин, рассматриваемых в теории вероятностей и математической статистике.

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

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

1. ИССЛЕДОВАНИЕ ВЛИЯНИЯ ТЕХНОЛОГИЧЕСКИХ ФАКТОРОВ НА ТОЧНОСТЬ ОБРАБОТКИ И ШЕРОХОВАТОСТЬ ПОВЕРХНОСТИ

Для лучшего уяснения методики применения математической статистики в технологических исследованиях рассмотрим следую­ щие примеры.

Пример 36.

Требуется установить, влияет ли изменение скорости резания

в диапазоне от

10 до 44 м/мин на величину разбивки отверстий при обработке

серого чугуна твердосплавной разверткой диаметром Dp = 20А. Разбивкой отвер­ стия Д называется разность между диаметром обработанного отверстия D 0 и диа­ метром развертки Dp, т. е.

А = D 0 — Dp.

Опыты проводятся при следующих неизменных условиях: станок револьвер­ ный, крепление развертки шарнирное, длина заготовок I = 2D0\ предваритель­ ная обработка отверстия производится на том же станке с одной установки заго­ товки по схеме: сверление, зенкерование, черновое развертывание; припуск на

чистовое развертывание

t = 0,05 мм на сторону; подача s = 1,68 мм/об; каждый

опыт с одним значением скорости резания v повторяется 4 раза.

Предварительными

опытами установлено, что рассеивание значений А

в большой партии обработанных заготовок подчиняется закону нормального распределения. Практическая кривая распределения величин А при разверты­ вании серого чугуна твердосплавной разверткой Dp = 20А приведена на рис. 30. Так как величина А является случайной и опытами установлено, что она подчи­ няется нормальному распределению, то наиболее обоснованным методом исследо­ вания влияния скорости резания на величину разбивки является статистический метод.

Каждый опыт, повторяемый 4 раза, будем рассматривать как случайную выборку объема п = 4 из нормальной совокупности. Для оценки влияния ско­

рости резайия на разбивку отверстий воспользуемся приемами математической статистики. Результаты опытов, а также вычисления среднеарифметических Xi и дисперсий ^ по каждой серии опытов для соответствующих значений V приве­

дены в табл. 35.

При общепринятых методах исследования график зависимости А от у будет иметь вид, показанный на рис. 31 (пунктирная линия). На основании графика рис. 31 можно сделать неправильный вывод о прогрессивном увеличении разбивки с ростом скорости в диапазоне от 10 до 44 м/мин и о существенном влиянии на раз­ бивку скорости резания в исследованном диапазоне.

&МК

тх

0,4

0,2

10 14 18 22 26

Ш м км

 

 

 

 

Рис. 30. Практическая кривая

распре­

Рис. 31. Зависимость А =

/ (и) при

деления величин разбивки А отверстий

общепринятых методах исследова­

при развертывании

 

ния (пунктирная

линия)

и

при

 

 

статистических

(сплошные

верти­

кальные линии)

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

Вследствие колебания величины силы Ру из-за неравномерности припуска на обработку и неравномерной твердости материала в отдельных заготовках про­ исходит колебание значений величины разбивки. Изменение скорости резания V в пределах от 10 до 44 м/мин не может существенно влиять на величину силы Ру, поэтому и величина разбивки не должна существенно изменяться при изменении скорости резания в указанных пределах.

Проверим эти предположения статистическими методами.

Таблица 35

Результаты опытов по исследованию влияния скорости резания на величину разбивки отверстия при развертывании серого чугуна

№серии опытов

Исследуемый фактор v в м/мин

1

10

2

24

3

33

4

44

Наблюденные значения

А=хс

Постоянные опытов:

 

 

 

 

2**

 

4

t в м м \ s

в м м /о б

1

2

3

4

 

 

 

 

 

 

/=0,05;

s= l,68

19

16

13

20

68

17

10

/=0,05;

s= 1,68

18

20

21

17

76

19

3,3

/=0,05;

5=1,68

18

19

21

22

80

20

3,3

/=0,05;

5= 1,68

21

20

23

24

88

22

3,3

Прежде всего установим, являются ли выборки случайными. Для этой цели воспользуемся способом последовательных разностей и определим с2, т. е. оценку а2 по данным выборки № 1 по формуле (ИЗ):

с2 = ~2(4'-

I)

[З2 + З2 +

72] =

11,17.

С другой стороны, мы имеем оценку о2 в виде s2 =

10. Определим критерий т:

т

С1

11,17

1, 12.

 

 

s2

10

 

 

По табл. 22 находим, что наблюденное т больше нижнего допустимого зна­ чения гд даже при 5%-ном уровне значимости. Следовательно, гипотеза «случай­

ности» выборки верна.

опытов с другими значениями V,

Проделав аналогичные вычисления для

получим следующие результаты вычислений с2 и т:

 

 

 

V в м/мин

10

24

33

44

с 2

11,17

3,5

2

2

S2

10

3,3

3,3

3,3

т

1,12

1,07

0,64

0,64

Т0,05

0,39

0,39

0,39

0,39

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

Теперь произведем оценку расхождений выборочных дисперсий sj и s?, =

= s3 = 4* Для этоя цели воспользуемся критерием G:

Gк

Sinax

10

0,5.

 

Е 4

19,9

 

 

 

 

В приложении 14 для п = 4 и п — 1 = 3G = 0,684.

Так как GH<^ G, то

гипотеза случайного расхождения выборочных дисперсий

верна.

Теперь можно проверить гипотезу однородности выборочных средних, т. е. установить, является ли наблюденное расхождение выборочных средних в опытах с различными значениями V случайными или это объясняется влиянием V. Для этого воспользуемся методикой, изложенной в части первой, гл. IV, п. 5, и оценим расхождение смежных значений выборочных средних по формуле (118):

tx

=

П 7 — 191

 

л /

4-4(4 +

4 — 2)

0,95;

V A (10 +

3,3)

У

 

4 +

4

 

 

 

 

/о =

 

119

— 201

 

-■ /4.4.(4 +

4 — 2)

 

*2 --

 

(3,3+ 3,3)

'

У

 

4 +

4

 

 

 

 

 

 

 

L =

 

120

— 22 |

'

I / 4.4 (4 +

4 — 2)

1,34.

*3

 

V 4 (3,3 +

3,3)

У

 

4 +

4

 

В приложении 5 Для

к =

п1-{- п2 — 2 = 6

этим значениям соответствуют

вероятности Р:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tx =

0,95;

Р =

0,38

 

 

 

 

 

t2 =

0,67;

Р =

0,53

 

 

 

 

 

U = 1,34;

Р =

0,23

 

 

Эти вероятности не малы, они больше доверительного уровня Р = 0,05, поэтому мо?кно считать, что наблюденные значения t случайно отличаются от табличных И гипотеза о случайном расхождении выборочных средних или их однородности верна. Другими словами — сравниваемые выборки принадлежат одной и той же генеральной совокупности.

Таким Образом, проверка гипотез однородности sj и X/ показала, что рас.

хождения Между наблюденными значениями

дисперсий

и средних в

опытах

с различными значениями скорости резания

V являются

случайными.

Следо­

вательно, н^Ше предположение, что скорость резания в пределах от 10 до 44 м/мин существенна не влияет на величину разбивки отверстий при развертывании серого чугука, подтверждается.

Приближенно можно считать, что при работе со скоростями резания У от 10

до 44 м/ми^> подаче s = 1,68 м/об и глубине резания t =

0,05 мм средние раз­

бивки отверстий будет равно

 

 

 

-

17+ 19 + 20 + 22

мкм.

 

Х =

4

= 19,5

 

При этом рассеивание значений разбивки характеризуется средним квадра­

тическим отклонением о =

z2s, где

 

 

 

s

10 +

3,3 + 3,3 +

3,3 =

2,2 мкм;

 

 

4

 

 

z2 по таблице приложения 10 равно 1,57. Следовательно,

ст= 1,57 2,2 = 3,4.

Истинное значение разбивки Д будет лежать в пределах X ± Зет, т. е.

X — За <

Д <

X + За;

 

 

19,5-10,2 <

Д <

19,5+10,2;

 

 

9,3 <

Д <

29,7.

 

 

График зависимости Д от V будет иметь вид, показанный на рис. 31 сплош­

ными линиями.

 

влияет ли на

чистоту

обработанной

Пример 37. Требуется установить,

поверхности изменение величины натяга

i

в пределах

от i = 0,30 мм до i =

= 0,6 мм при центробежно-шариковой обработке колец из стали

45 *•

Предварительная обработка колец с чистотой по 3-му классу

(Rz = 43 -f-

46 мк) произведена на токарном станке. Для разрешения поставленной задачи 90 колец были подвергнуты центробежно-шариковой обработке при неизменном числе ударов на 1 мм2 поверхности, неизменной продольной подаче на 1 оборот детали, но с различными натягами i, которым придавались следующие значения: i = 0,3; 0,45; 0,60 мм.

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

Способ центробежно-шариковой обработки предложен И. И. Ененко.

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

Измерения шероховатости

поверхности колец, обработанных при различных

значениях i, дали следующие

результаты:

 

 

i =

0,3 мм; R z — 8 мкм;

S2Rz =

0,53 мкм2;

i =

0,45

мм;

Rz =

7 мкм;

S ^ 2 =

0,63 мкм2;

i =

0,60

мм;

Rz =

5,3 мкм;

0,42 мкм2.

Для оценки степени влияния величины натяга i на чистоту обработки будем исходить из нулевой гипотезы, т. е. предположим, что величина натяга в пределах от 0,3 до 0,6 мм не влияет на чистоту обработки. Проверим эту гипотезу сначала для i = 0,3 и i = 0,45, а затем для i = 0,45 и i = 0,6 по формуле (119):

к

v

18 — 71

5,1 •

0,53

 

 

 

 

 

0,63

 

 

 

30

f

30

к

 

1 7 - 5 ,3

9,3.

- i /

 

0,63

 

 

 

 

0,42

 

У

 

зо

+

зо

В обоих случаях полученные значения t > 3, следовательно, можно утвер­ ждать, что наблюденные значения Rz в опытах с различными величинами i суще­ ственно отличаются друг от друга, что объясняется влиянием величины натяга i. Физический смысл этого явления заключается в том, что с увеличением натяга увеличивается сила удара шариков, а поэтому возрастает и величина пласти­ ческой деформации гребешков, полученных на предшествующей обработке, и снижается шероховатость поверхности.

2.ОБРАБОТКА ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ ПО СПОСОБУ НАИМЕНЬШИХ КВАДРАТОВ

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

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

118

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

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

Пусть имеются некоторые переменные х, у и z, связанные

между собой уравнением вида

 

ах + $у + yz = N.

(183)

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

a*i + pt/i + yzt = Л^;

ах2 +

рг/2 +

yz2 =

N 2;

(184)

а* з +

Pl/з +

У2з =

N 3, .

 

имеющих три неизвестных ос, (3, у, значения которых можно опре­ делить, решив систему трех уравнений. Если переменным х , у, г придать некоторые четвертые значения: х4, г/4, z4 и определить опытным путем значение W4, то окажется, что равенство (183) при значениях коэффициентов, вычисленных по уравнениям (184), удовлетворено не будет:

ах4+ Ру4 + yzi Ф N.

(185)

Такой результат объясняется следующим. Если бы коэффи­ циенты были определены из уравнений (184) Па основании мате­ матически точного обмера переменных х, уу г и N у равенство (183) было бы удовлетворено и при других значениях переменных ху у у z. Однако всякие измерения являются лишь относительно точ­ ными, так как экспериментатор и прибор, которым он пользуется, дают ряд неизбежных ошибок. Кроме того, значения измеряемых величин могут колебаться в каждом опыте под влиянием случай­ ных причин. Ошибки измерений и колеблемость измеряемых ве­ личии приводят к результатам, выраженным неравенством (185).

119

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

В этом случае получим систему п уравнений:

«*1 +

Р*/ 1 +

у г х=

N i ,

а х2 +

$у2 +

?z2 =

N 2;

 

 

 

( 186)

ахп +

Руп +

угп =

Nn.

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

Таким образом, возникает вопрос, каковы должны быть зна­ чения коэффициентов а, р, у, чтобы избыточная система уравне­ ний удовлетворялась наилучшим образом, т. е. ошибки Et = N

Nlt где N — расчетные, a Nt — опытные данные или

=a*i — ру х — yzi — N 1, '

Е2 = ах2 — ру2 — уг2 N 2;

Еп = <ххп — руп — угп — Nn ,

должны быть наименьшими.

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

те, при которых сумма квадратов

ошибок

будет наименьшая:

Ё{ + Е\ + • • • +

= min.

 

Возведение в квадрат, например, первого уравнения системы (187) дает

Е\ — а х \ + 2ах$у1+ fty\ -f-yx^i + 2ру ^ г г

2axxN1 2PyiNi -f- у2^ 2yZiN: + N\.

Для El, Е\, Е2п получим такие же выражения с той лишь разницей, что у величин х, у, z и N будут индексы от 2 до п.

120

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