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

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

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

Сложив

эти уравнения, получим

 

 

 

£ E

] =

a 2 £ x *

+

tf£ y '> i +

y > £ z 2i +

2 a v i x iyi +

 

1

 

1

1

1

 

 

1

 

 

+

2<ху Еп1

X[Z{ + 2PY 1пЕ

Уf i t

— 2а Епх х £Л^

 

 

 

- 2 P

E

^ - 2 Y E Z^

+

I > ? .

(188)

 

 

 

1

 

1

 

1

 

п

Для соблюдения условия 21 £? = min необходимо, чтобы сумма

частных производных по а, (}, у уравнения (188) обращалась в нуль, а следовательно, и каждая из этих частных производных должна обратиться в нуль:

1

= 0,

э Е * ?

= о,

ар

да

 

 

 

д%Е)

 

 

 

х

0.

 

 

ду

 

 

 

 

Тогда значения а,

р и у, которые обращают каждую из этих

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

Дифференцируя уравнение

(188) последовательно

по а, р,

Y, получим

Р Ъ ху +

Y S*z — Е xN = 0;

 

а Е *2 +

 

Р Е У2 +

а Цху +

у E«/z — Ef/N = 0;

(189)

Y E z2 + a E*z + P Ei/z — EzlV = 0..

Эти уравнения носят название «нормальных». Решая их отно­ сительно а, р и у, получим наилучшие значения коэффициентов, удовлетворяющих принципу наименьших квадратов.

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

Если уравнение имеет вид линейной зависимости

у = а + Ьху

то для определения коэффициентов а и b по способу наименьших квадратов необходимо по аналогии с предыдущим составить по экспериментальным данным сумму:

Ё

= Ё (yt — y)2 =

ti( y i — a — bxi)2

(190)

i

i

i

 

Возведя правую часть уравнения (190) в квадрат и взяв част­ ные производные по а и & и приравняв их нулю, получим

 

1

=

2

' Ё У1

+ ь Ё Xi + Ё а

=

0;

 

да

 

 

 

 

 

-.2

 

 

 

 

 

 

 

 

 

 

 

 

« 2 3 * ? _

— Ё х Mi

 

Ь Ё + а Ё Xi

 

= 0.

откуда

дЬ

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

па + *23^ = 23^;

 

 

 

 

 

 

 

1

 

1

 

 

 

 

(191)

 

 

 

п

 

п

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

23xi ь 23 x*i +

 

S xiiJi-

 

 

 

Система уравнений (191) очень легко решается, если предва-

рительио

заменить

xt на x-L= (xt — X), где X — — 2 */• Тогда

исходное

уравнение

примет

вид

 

 

 

 

 

 

П

 

у = а' +

 

Ь (х, — Х).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как 23*; =

0, то уравнения

(191)

примут

вид

 

 

 

 

па' =

23 и ;

 

 

 

откуда

 

 

ь 23

*;/2

 

= 2

3

*<•*/;>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

23 У‘

 

 

 

 

 

 

 

 

 

а

=

_1___

 

 

 

(192)

 

 

 

 

 

п

»

 

 

 

 

 

 

23 x\yi

 

п

 

23

 

 

 

 

 

 

 

 

 

 

 

 

(193)

 

 

ь = -------- = _*---------------

 

 

 

 

 

1

#2

 

Ё (* -* )*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

_

а Н- Ьх служит равен­

Для перехода от о! к а уравнения у =

ство

 

 

 

 

 

 

 

 

 

 

 

 

а = а' — ЬХ.

Рассмотрим еще один случай, когда уравнение имеет вид пара­ болы третьего порядка:

и = р -I- «г* + а 2х2 + а ах3.

(194)

Чтобы сумма квадратов отклонений

S £ 2 = S (у — Р — а хх — а 2л:2 — а3х3)2

(195)

была наименьшей, необходимо выполнить четыре условия (по числу неизвестных коэффициентов а 1, а 2, а 3> Р):

д Ц £ 2

= 0;

з Ц £ 2

0;

з £ £ 2

0;

ЗЦ £2= 0.

ар

 

дс^

 

дао

 

аа3

Возведя в квадрат правую часть уравнения (195) и взяв част­ ные производные от него по а х, а 2, а 3, р и приравняв их нулю, получим четыре нормальных уравнения:

я

я

п

п

'

^Р + а 1 Xj Xi -f- 0С2 Xj Xi +

a 3 Xj

= Xj Уь

 

1

1

1

1

 

PXJ

+

ai XJ

+ a 2 XJ

+

аз XJ

— XJ Xitji\

1

 

1

1

 

1

1

(196)

л

 

я

я

 

я

я

Р XJ

~h ai XJ

+ a 2 XJ

4" а з XJ ** = XJ

 

1

 

1

1

 

1

1

 

я

 

я

я

 

я

я

 

Р X]

+

a lXj х \ + 0^2XJ ** +

a3 XJ ** == XJ x ty i.

1

 

1

1

 

1

1

1

Решив эти уравнения, найдем искомую параболу третьего порядка.

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

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

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

-

У , П хуУ

 

У х ~

 

где пХуУ — сумма произведений значений у

на их частоты для данного значе­

ния х;

— сумма частот значений у для данного значения х.

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

ИЗ X будет иметь вид

__

 

Ух = а + bx + а*2,

Корреляционная таблица
х в мкм
Таблица 36

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

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

Е Е3 = 2 (у — а bx — сх2)2,

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

ровать обе части уравнения. Коэффициент при параметре а ра­

вен 1, при параметре b равен х и при параметре с равен х2. Поэтому после

 

 

 

 

у в мкм

5

10

15

20

 

 

 

 

 

 

 

 

 

 

 

 

 

50

 

3

2

1

3

Ух

 

 

 

 

40

 

6

6

15

 

 

 

 

30

3

8

6

3

20

 

 

 

 

 

 

<Г\

 

 

 

20

5

4

1

 

10

 

 

 

 

 

10

2

 

*

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Е

п*

10

15

15

10

50

 

 

 

 

J j nxy

210

440

540

380

 

0

5

10

15

20х

 

21

29,3

36

38

 

Рис. 32. Эмпирическая (а) и

теорети­

Ух

 

ческая

(б) кривые

регрессии

у на х

 

 

 

 

 

 

умножения на указанные коэффициенты и суммирования обеих частей урав нения

у = а + Ьх + сх2

 

получим систему трех нормальных уравнений:

 

Е У= па+ Ь ^ х + сЕ*2;

 

= ^Е* + 6E*2 + сЕ * 3;

(197)

2 * 2г/ = а£х2+ bj^x3+

 

Чтобы вычислить параметры а, b и с, необходимо составить вспомогательную табл. 37. Для упрощения расчетов, особенно когда значения х и у выражены мно­ гозначными числами или десятичными дробями, рекомендуется заменить дей­ ствительные значения х и у новыми значениями х ' и у', выраженными в виде про­ стых чисел и определяемых по формулам (142):

, = — дх

= у — ду

сх 9 У

су 9

где ах

и

ау — новые начала отсчетов

для х и у\

сх

и

су — интервалы значений х

иjy .

В нашем примере принято: ах = 10;

сх — 5; ау — 30; су — 10. В результате

указанных замен уравнение регрессии у

на х примет вид

—/

/

, ,

/

/2

Ух = а

+

6 х

+ с х

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

=па' + 6']£>' + с'2>'2-

=

а '^ х ’ I

b'J^x'3 +

с'][]х'я;

098)

5 > 'V =

а 'Е * '2 +

Ь '^ х '3 +

с' 2 > '4

 

Расчет всех сумм, являющихся коэффициентами при неизвестных и свобод­ ными членами в уравнениях (198), произведен в табл. 37.

 

Таблица 37

Расчет теоретической линии регрессии у х’ =

а'+ бх'+ с'л;'2

X

Номера столбцов

У'

21

 

1

 

0

 

—2—1>

 

 

1

 

2

 

3

 

4

о

5

Си

6

в

а

7

си

S

8

О

 

0J

 

X

У

50

40

30

20

10

Итог пх

пх х!

пхх ’2 пхх '3

пхх 'А

^ц пхуУ'

х'% П Ху у '

х’ 2^ п х у у '

-1

0

1

2

1

 

 

2

 

 

X

 

 

 

5

 

 

 

 

 

 

10

15

 

пу

 

 

 

 

 

 

20

итог

 

Пу-У'

 

 

 

 

 

_

3

2

1

3

 

 

6

6

6

15

 

 

15

3

8

6

3

20

 

 

0

5

4

1

10

 

—10

2

2

 

—4

10

15

15

10

50

 

 

7

—10

0

15

20

25=]С*'

 

10

0

15

40

65=]Е>'2

—10

0

15

80

85=]Е>'3

10

0

15

160

185= Е * '4

—9

—1

9

8

7= 5 > '

 

 

9

0

9

16

34=]

to

<5^

—9

0

9

32

 

32=]LTJ

 

Подставляя числовые их значения из табл. 37 и имея в виду, что

п —

 

 

~ 50>

получим систему нормальных уравнений:

 

7 =

50а' +

25Ь' +

65 с';

34 =

25а' +

656' +

85с';

32 =

65а' +

856' +

185с'.

Решая эту систему, получим а = —0,0233; Ь' = 0,740; с' = —0,159. Зави­ симость между ух и х' выразится уравнением

~Ух = —0,0233 + 0,740л:' — 0,159л:'2-

Произведя в этом уравнении замену переменных

 

х' =

Ух

У — 30

 

10 ’

получим

 

 

 

Ух — 30

— 0,0233 +

0,740 —-д— — 0,159

10

 

 

 

откуда

ух = 8,607 + 2,752л; — 0,0636л:2.

На рис. 32 по этому уравнению построена теоретическая кривая регрессии у на х (показанная сплошной линией).

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

Р = Cs“/P,

(199)

в котором надо определить а, р и С путем логарифмирования, приводим к виду

Ig Р = lg С + a lg s + р lg t.

(200)

Полагая

 

lg Р = N; lg s = х; lg t = y, lg C = 2 и у =

1,

получим уравнение вида (183).

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

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

(200). Прологарифмировав их, получим

п условных уравнений:

a lg s, +

р lg tx +

lg С =

lg Pi,

a Ig S2 +

P lg +

lg C =

lg P 2,

a lg sn +

P lg tn +

lg C =

lg Pn.

Чтобы получить нормальные уравнения, воспользуемся опи­ санным выше механическим способом, т. е. умножим каждое из условных уравнений на коэффициент, стоящий при а, затем сум­ мируем их; то же сделаем для (3 и lg С. В результате получим

для а:

a S (lg s)2 + Р Е lg S lg t + lg C E lg s = E lg S lg Л

для (i:

a S lg s-lg * + P £ (lg 02 + ^ C S lg t =

= S i g M g ^ ;

для lg C:

а I] lg s + p Ц lg / + /г lg С = Ц lg P.

Решив систему этих уравнений, определим искомые а, р и С. По внешнему виду эмпирической кривой можно подобрать не­ сколько теоретических кривых и, следовательно, несколько видов формул. Например, для эмпирической кривой, изображенной пунктирной ломаной линией на рис. 32, можно подобрать две фор­ мулы: для параболы второго порядка у = а + Ьх + сх2 и для прямой у = а + Ьх.

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

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

значениями у1ш Та

функция, для которой основная ошибка а0

будет наименьшей,

принимается в качестве наилучшей:

 

(202)

где т — число значений yt.

Пример 39. По данным примера 38 установить, к какой из двух возможных

теоретических функций: ух = а + Ьх + сх2 или ух = а + Ьх является более близкой эмпирическая кривая, показанная на рис. 32.

Для вычисления параметров функции ух = а + Ьх составим систему нор­ мальных уравнений:

2 ^ = па + bJjx* J^xy = aj^x + b2 * 2-

Для упрощения расчетов заменим х на х' и у на у' так, как это было сделано в примере 38. Тогда система уравнений примет вид

YiP' = па' + b'Jjх’\

7 = 5 а '+ 256';

34 = 25а' + 656'2.

Решив эти уравнения, получим

а '= —0,15 и Ь' = 0,58.

Следовательно,

ух = - 0 , 1 5 + 0,58*'.

Заменив переменные ух на ух, х' на х, получим

У х — 30

V___ 1О

— 0,15 + 0,58— =—

10

о

или

У х = 16,9+ 1,16*.

Вычислим теоретические значения Y 0i для уравнения прямой и для уравне­ ния параболы Y 0i = 8,706 + 2,752 х — 0,0636 х2 и сумму квадратов отклонений их от эмпирических значений ух. Результаты вычислений приведены в табл. 38.

аи= Y ^ = h б5;

] /+ - = ° .б 9 .

Таблица 38

Значения Y 0 и (ух Y 0) 2

 

 

Для прямой

 

 

Для

параболы

 

 

Ух

^0

У х - у о

( У х ~ у о ) 2

У х

Го

У х ~ у 0

& - у о)2

5

21

22,7

1,7

2,90

21

20,8

0,2

0,04

10

29,3

28,5

0,8

0,64

29,3

29,7

0,4

0,16

15

36

34,3

1,7

2,90

36

35,9

0,1

0,01

20

38

40,1

2,1

4,40

38

39,1

1,1

1,21

 

 

 

 

10,84

 

 

 

1,43

Следовательно, уравнения параболы второго порядка более точно выражают характер эмпирической кривой, так как <т02 < сг01.

СТАТИСТИЧЕСКИЙ АНАЛИЗ ТОЧНОСТИ МЕХАНИЧЕСКОЙ ОБРАБОТКИ

ИСТАТИСТИЧЕСКОЕ РЕГУЛИРОВАНИЕ ТЕХНОЛОГИЧЕСКИХ ПРОЦЕССОВ

1. По г р е ш н о с т и м е х а н и ч е с к о й о б р а б о т к и

И ЗАКОНЫ ИХ РАСПРЕДЕЛЕНИЯ]

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

Погрешности обработки принято делить на три вида: погреш­ ности размера, погрешности формы и погрешности взаимного рас­ положения поверхностей и осей детали.

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

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

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

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

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

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

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

9

И. С. Солонин

129

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

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

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

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

При работе на настроенных станках или при работе мерным режущим инструментом размерный износ систематически изме­ няет размер каждой новой детали на удвоенную величину износа инструмента при выполнении диаметральных размеров или только на величину износа при выполении размеров длин. Исследова­ ниями установлено, что размерный износ и режущих инструмен­ тов протекает во времени т по закону, выражаемому кривой У, показанной на рис. 33. Эта кривая имеет три участка: первый участок (от 0 до т„) характеризует изменение износа инструмента в период его приработки; это так называемый начальный износ инструмента; второй участок (от т„ до т*) характеризует нормаль­ ный износ инструмента и третий участок (от тк до конца кривой) — катастрофический износ, когда наступает быстрое разрушение инструмента.

Если работа производится предварительно доведенным инстру­ ментом, размерный износ будет изменяться во времени по закону прямой (см. рис. 33, линия 2). При этом период стойкости инстру­ мента увеличится, так как доводка инструмента повышает его износостойкость в 1,5 и более раз.

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