Министерство образования и науки Российской Федерации

ГОУ ВПО Магнитогорский Государственный Технический Университет

им. Г. И. Носова

Кафедра обработки металлов давлением

Курсовая работа

На тему

«Определение аналитической зависимости сопротивления металла пластической деформации для стали 30ХГСА»

По дисциплине

«Организация и планирование эксперимента»

г. Магнитогорск 2009г.

Задание

На основании базисного значения пластичности меди (? _о.д =105 МПа для стали 30ХГСА) с использованием графиков термомеханических коэффициентов Кt, Ке, Кu, составить аналитическую зависимость, позволяющую определить сопротивление деформации (?_т) при горячей прокатке непосредственно от величин температуры, скорости и степени деформации для стали 30ХГСА.

Введение

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

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

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

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

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

?т=?_0.д.*Kt*K ? *Ku (1)

где ?_0.д - базисное значение сопротивления деформации;

Kt - температурный коэффициент;

K? - степенной коэффициент;

Ku - скоростной коэффициент.

Графики зависимости коэффициентов Kt, K? и Ku от температуры, степени деформации и скорости деформации приведены в справочниках.

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

П.Л.Клименко путем аппроксимации обобщенных кривых изменения k_t, k и k_u вывел формулы зависимости коэффициентов от температуры, степени и скорости деформации:

(2)

для (3)

для (4)

для с^-1 (5)

для с^-1, (6)

Нам надлежит найти похожие уравнения зависимостей Kt=f(t), K?=f(?) и Ku=f(u), используя метод парного регрессионного анализа. Парный регрессионный анализ - это метод математической статистики, который позволяет найти отображение (модель, аппроксимацию) стохастической зависимости между откликом У и фактором Х.

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

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

yi=?(xi)+?i (7)

где ?(xi) - систематическая (объясненная) составляющая; она обусловлена существованием связи между откликом и фактором;

?i - случайная составляющая; она обусловлена разнообразными возмущениями и вызывает отклонение уi от соответствующих реальной зависимости.

Относительно ?i делают следующие предположения:

- это нормально распределенная случайная переменная.

- ?(?i)=0 (математическое ожидание случайной составляющей равно нулю).

- ?(?i)=const (дисперсия случайной составляющей постоянна).

- в различных наблюдениях значения ?i не зависят друг от друга.

Задача определения вида уравнения регрессии состоит в нахождении систематической составляющей ?(xi).

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

(8)

Объясненная дисперсия характеризует рассеяние уi, обусловленное зависимостью отклика от фактора; ее находят по формуле

(9)

После расчетов проверяют соответствие полученного уравнения опытным данным по критерию Фишера

(10)

где a- уровень значимости;

k- число наблюдений.

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

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

Множественный регрессионный анализ - это метод математической статистики, позволяющий найти наиболее точное и достоверное отображение (модель, аппроксимацию) стохастической зависимости между откликом Y и факторами Х₁, Х₂, …,Х_j,X_m.

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

Коэффициенты регрессии b_j являются случайными величинами с математическими ожиданиями ?_j и дисперсиями, которым соответствуют стандартные отклонения S_bj. Значение b_j признается статистически значимым, если выполняется условие

(11)

где t_bj и t[a;n-k] - расчетное и табличное число Стьюдента.

Если условие (11) не выполняется, то влияние фактора Х_j на отклик несущественное, и его надо исключить из дальнейших расчетов.

Статистическая надежность выбранного уравнения регрессии проверяется также по критерию Фишера (уравнение (10).

После проведения множественного регрессионного анализа получаем уравнение вида

?т = f (t, U, ?), (12)

Находим при помощи полученного уравнения сопротивления металла деформации третьим способом.

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

Краткая характеристика стали 30ХГСА

Из стали 30ХГСА изготовляют различные улучшаемые детали: валы, оси, зубчатые колеса, фланцы, корпуса обшивки, лопатки компрессорных машин, работающие при температуре до 200°С, рычаги, толкатели, ответственные сварные конструкции, работающие при знакопеременных нагрузках, крепежные детали, работающие при низких температурах. Склонна к отпускной способности, флокеночувствительна

Таблица 1 - Химический состав стали 30ХГСА,(%).

Исходные данные

Данные с исходных графиков приведены в таблице 2.

Таблица 2 - Исходные данные данные.

Определение уравнений зависимости термомеханических коэффициентов от их физических величин

Проведем парный регрессионный анализ. Рассмотрим по 5 уравнений для каждой зависимости. Расчеты удобно проводить в среде электронных таблиц MS Excel. Результаты оценки пяти уравнений представлены в таблицах 3-5. В таблицах жирной строкой выделено то уравнение, которое является наилучшей аппроксимацией исследуемой зависимости. Для температурного коэффициента это логарифмическая зависимость, для коэффициента деформации - степенная зависимость, для скоростного коэффициента - логарифмическая зависимость. При выборе уравнения ориентировались на критерий Фишера F_расч принимающий максимальное значение, а также условие F_расч>F_табл.

Графики выбранных уравнений приведены на рисунках 1-3 . На рисунках точками изображены значения, полученные по исходным графикам зависимостей термомеханических коэффициентов от их физических величин. Сплошными линиями показаны графики полученных уравнений аппроксимации.

Таблица 3 - Уравнения зависимости Кt от t

Таблица 4 - Уравнения зависимости К? от ?

Таблица 5 - Уравнения зависимости Кu от U

Рисунок 1. Температурный коэффициент для стали 30ХГСА

Рисунок 2. Степенной коэффициент для стали 30ХГСА

Рисунок 3. Скоростной коэффициент для стали 30ХГСА

По данным, полученным в результате парного анализа (таблица 6) строим графики (рисунки 4-6).

Таблица 6 - Данные полученные в результате парного регрессионного анализа

Рисунок 4. Температурный коэффициент для стали 30ХГСА

Рисунок 5. Степенной коэффициент для стали 30ХГСА

Рисунок 6. Скоростной коэффициент для стали 30ХГСА

Определение уравнения зависимости сопротивления деформации непосредственно от физических величин

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

Таблица 7 - Исходные данные для составления уравнения ?т = f (t, U, ?)

По подготовленной таблице в MS Excel с помощью функции «Регрессия» из пакета анализа данных проводим множественный регрессионный анализ.

В результате получаем уравнение ?_т= 390,20 - 0,33*T +2,21*E +0,98*U

Для выяснения статистической значимости коэффициентов уравнения сравниваем рассчитанные коэффициенты Стьюдента с табличными для числа наблюдений 10 и уравнения с четырьмя коэффициентами и доверительной вероятностью 95%. Коэффициенты Стьюдента, рассчитанные для коэффициентов t, ? , u оказались больше табличного коэффициента Стьюдента, то есть, статистически значимыми.

Для выяснения надежности аппроксимации полученным уравнением сравниваем рассчитанное число Фишера с табличным для степеней свободы (10-4=6) и доверительной вероятностью 95%. Рассчитанный критерий Фишера оказался больше табличного, значит, уравнение достоверно отражает исследуемую зависимость. Лист MS Excel с расчетом представлен на рисунке 7.

Рисунок 7 - Лист МS Excel

Для построения графика полученного уравнения с двумя неизвестными на плоскости нужно, чтобы один из факторов принимал два экстремальных значения, а второй непрерывно изменялся. Тогда графически это уравнение можно представить как область в координатах «изменяющийся параметр - предел текучести». Такие графики для полученного уравнения представлены на рисунках 8-10.

1. Третьяков А.В., Зюзин В.И. Механические свойства металлов и сплавов при обработке металлов давлением. Справочник. М. Металлургия. 1973. 224с.

2. Веников В.А., Веников Г.В. Теория подобия и моделирования. Учебник для вузов. М. Высш.шк. 1984. 439с.

3. Моллер А.Б. Методические указания для выполнения курсовой работы по дисциплине «Теория подобия и физическое моделирование». Магнитогорск: МГТУ, 2000.

4. Моллер А.Б., Синицкий О.В., Назаров Д.В.Моделирование процессов ОМД с применением планирования факторного эксперимента. Методические указания для самостоятельной работы, практических занятий и выполнения расчетно-графических работ по дисциплине «Планирование и организация эксперимента». Магнитогорск: МГТУ, 2008.

5. Румянцев М.И. Методические указания к лабораторной работе по дисциплине «Обработка и анализ числовой информации». Магнитогорск: МГТУ, 2003.