Файл: Технические и экономические основы литейного производства..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 17.10.2024

Просмотров: 178

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Исходные данные для расчета коэффициентов регрессии (планирование второго порядка для т — 2, ... 5 переменных)

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

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

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

52

Матрица планирования предусматривает проведение N опы­ тов. При планировании первого порядка N = 2 т , где т — число независимых переменных; помимо этого сверх плана осуществ­

ляется N 0 опытов на основном уровне независимых

переменных

для установления ошибки эксперимента

(N 0 =

3- ^4) .

При пла­

нировании второго порядка N = 2 т + 2 т

+ N 0,

где 2 т,

как и при

планировании первого порядка,— основной план;

2 т — число

опытов в «звездных» точках; N 0 — число опытов в центре экспе­ римента. N o зависит от числа независимых переменных и опре­ деляется по таблице исходных данных, приведенной на стр. 52. По этой же таблицё выбирается плечо «звездных» точек а.

Статистические методы планирования эксперимента позволя­ ют осуществлять не весь полный факторный эксперимент, а толь­ ко его определенную часть, называемую дробной репликой. Здесь используется возможность замены парных и других взаи­ модействий параметров новой независимой переменной (напри­ мер, отказавшись от рассмотрения совместного влияния парного взаимодействия ХіХ 3, на его место можно поставить новую пере­ менную х4). В таком случае число опытов N = 2 т~Р, где т р характеризует полноту эксперимента; при /7 = 0 образуется пол­

ный факторный эксперимент, при р = 1,2, ..., образуются его дробные реплики: половина плана, четверть плана и т. д.

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

Для установления ошибок эксперимента проводится ряд опы­ тов на основном уровне независимых переменных [36]. Ошибки эксперимента определяются раздельно для каждого показателя процесса как их среднее квадратическое отклонение:

где у 0іи — и - е значение показателя

процесса

на основном уров­

не независимых переменных, и = 1,

2, ..., N 0 -,

N 0 — число опытов

на основном уровне; у 0 — среднее значение показателя процесса на основном уровне независимых переменных.

5. Анализ результатов эксперимента, проведенного с помощью указанных выше методов, позволяет функциональную зависи­ мость (1) представить следующим полиномом:

тт

у = ь 0+ 2

ь іх і + 2 Ьі'ХіХі ’

(4)

і=1

£, /= 1

 

где Ь0, Ьі , Ьц , а также Ь ц при / = і — искомые

коэффициенты

регрессии; x tXj — парные взаимодействия изучаемых факторов

53


(независимых переменных), включая их квадраты, т. е.

= x f

при / = і.

 

Коэффициенты регрессии вычисляются по приведенным ниже формулам в зависимости от типа планирования. Для планирова­ ния первого порядка, когда отсутствуют квадратичные члены (обязательное условие і ф ] ' ) , коэффициенты регрессии опреде­ ляются по следующим формулам:

N

 

 

 

 

 

 

 

2

У“

 

 

 

 

 

 

 

 

 

 

 

 

 

и=

1

 

 

 

 

 

(5)

 

 

 

 

 

 

 

N

'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

Х ІиУч

 

 

 

 

 

 

 

 

 

 

 

ь , =

Ц=1

 

 

 

 

 

 

(6)

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Х ІиХ ]'иУи

 

 

 

 

 

 

 

 

 

 

6 ,.=

и = I

N

 

 

 

 

 

(7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффициенты регрессии при планировании второго порядка

определяются по более'сложным формулам:

 

 

 

 

 

 

 

N

 

 

 

 

т

N

 

 

 

 

 

 

V

^

«

у и

N - N 0 V I V I .2

 

 

(8)

 

Ь °

2

 

 

m c m N 0

 

 

 

 

 

 

~ N o

 

 

 

 

1-1

W=l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b , =

 

I ХіиУа I

 

 

 

 

(9)

 

 

 

 

 

 

т с„

 

 

 

 

 

 

(10)

 

 

 

 

 

 

 

и - 1

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ь „ - Л = » *

m ( m + 2 ) N 0 ^ x f uy u + ( 2 N — 2 N 0 — m N 0) x

 

2N0m 12 c%

 

 

 

 

и - I

 

 

 

 

 

 

 

 

m

 

 

N

 

 

 

N

 

 

 

 

 

 

 

X 2

 

2

х ш У “ ~

2 m c m 2

lJu-

 

 

 

(П)

 

 

1=1

 

1

 

 

 

1

 

 

 

 

где и

= 1,2, ...,

Л( при планировании первого

порядка

и

и = 1,

2,...,

N 0, ..., N при планировании второго порядка;

і = 1,2,

..., т

число независимых переменных;

 

N

x j a—

постоянное чис-

c m = ^

ло для эксперимента при т

 

 

U-1

 

 

(для т =

независимых переменных

= 2

с т = 8,0; для т =

 

3

ст = 13,658 и т. д) .

 

 

 

 

При планировании второго порядка для определенного числа

независимых переменных формулы

(8) — (11)

значительно упро­

щаются и принимают вид

 

 

 

 

 

 

 

 

 

54


N

т

N

ь0= в

 

( 12)

 

Ы

(13)

Ьі B i

ХіаУи. >

 

U=I

 

 

N

(14)

 

и- 1

 

 

N

rn N

N

где значения множителей В являются постоянными величинами для заданных т и находятся по исходным данным, приведенным

втаблице на стр. 52.

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

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

Следующим этапом регрессионного анализа является опреде­ ление значимости коэффициентов регрессии, которая оценивает­ ся сравнением каждого Ьі (Ь0, Ьц , Ьц ) с критическим коэффици­

ентом значимости Ьк р ( Ьі ^ Ькр):

где

t q j E — табличное значение критерия t Стьюдента для оп­

ределенных q и f E ,

q — уровень значимости

(принимается 5 или

1 0 % - н ы й

уровень

значимости); fE =

N о

I — число степеней

свободы;

s {y } — ошибка эксперимента,

определяемая по форму­

ле (3).

 

 

 

 

 

Проверка адекватности модели осуществляется на основании критерия F Фишера [50]:

 

(17)

где

Fq (iR . i E) — табличное значение критерия F при определен­

ных q,

f R, f E; f R = N к — I — число степеней свободы (/с — чис­

ло статистически значимых коэффициентов); s r — остаточная

5 5


дисперсия

^ s 2r= ~ y ^ J >зДесь S r — остаточная сумма квадратов

(

N

к

\

\

и=I

і=0

)

 

Так как статистические модели лишь приближенно оценивают

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

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

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

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

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

56


с тем естественно, что значения параметров и диапазон их изме­ нения непосредственно зависят от вида сплава, конструкции и требований к качеству отливок и условий их производства. По­ этому объем работ по моделированию литейных процессов до­ статочно велик даже с учетом применения современных матема­ тических методов планирования экспериментов. Ниже приводят­ ся несколько примеров построения математических моделей на базе указанных методов для различного числа независимых пе­ ременных и различных типов планирования [16, 42, 46].

При постановке конкретных задач были использованы лите­ ратурные источники [5, 15, 19, 48].

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

Чугун выплавляли в индукционной печи повышенной частоты мод. МГП-52 с кислой футеровкой тигля емкостью 80 кг, частотой 2500 Гц. Контрольные плавки и выборочная проверка результа­ тов исследования осуществлялись в индукционной печи про­

мышленной частоты

мод. ИЧТ-1

(емкость тигля 1 т).

Шихта состояла

из кусковой

стали 15,

чушкового чугуна

ЛК-4,

электродного

боя и

ферросплавов

(ферросилиций —

73,3%

Si, феррофосфор — 20%

Р,

ферромарганец — 67,5% Мп).

Модифицирование чугуна производилось дробленым модифика­ тором ФС75 в ковше при температуре перегрева.

В процессе эксперимента определялись следующие свойства чугуна: предел прочности при разрыве ств в кгс/мм2, предел проч­ ности при изгибе а„ в кгс/мм2, стрела прогиба f 300 в м м , твердость Н В , жидкотекучесть металла Я в см, глубина отбела b в мм.

Для определения механических свойств (ав, а„, /зоо, Н В ) зали­ вались стандартные образцы диаметром 30 мм и длиной 330 мм. Жидкотекучесть определялась по спиральной пробе, отбел — по клину.

В качестве независимых переменных

были приняты следую­

щие технологические факторы:

 

 

 

 

С, Si, Mn, Р — содержание химических

элементов в чугуне в

процентах и соответственно Xi,

x%,

x3, х 4 — в от­

носительных единицах;

 

 

 

 

А^п — перегрев металла над ликвидусом

в °С — соот­

ветственно х 5;

 

 

 

 

тв— продолжительность выдержки

перегретого

ме­

талла в печи в мин — соответственно х6;

 

7м — расход модификатора ФС75 в процентах от мас­

сы металла — соответственно х 7;

 

 

q / Q — доля

стали (стальных отходов)

в

шихте (Q —

масса

садки) — соответственно

х8;

 

А/з — превышение температуры

заливки жидкого

ме­

талла над ликвидусом в °С — х 9;

 

 

57