///Исследование свойств многокомпонентной стали в системе Tibco STATISTICA

Исследование свойств многокомпонентной стали в системе Tibco STATISTICA

Сталь – деформируемый (ковкий) сплав железа с углеродом (и другими элементами), содержание углерода в котором не превышает 2,14%, но не меньше 0,02%. Углерод придаёт сплавам железа прочность и твёрдость, снижая пластичность и вязкость. Учитывая, что в сталь могут быть добавлены легирующие элементы, сталью называется содержащий не менее 45% железа сплав железа с углеродом и легирующими элементами (легированная, высоколегированная сталь).

Оглавление

О задаче

Сталь – деформируемый (ковкий) сплав железа с углеродом (и другими элементами), содержание углерода в котором не превышает 2,14%, но не меньше 0,02%. Углерод придаёт сплавам железа прочность и твёрдость, снижая пластичность и вязкость. Учитывая, что в сталь могут быть добавлены легирующие элементы, сталью называется содержащий не менее 45% железа сплав железа с углеродом и легирующими элементами (легированная, высоколегированная сталь).

Предел текучести

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

 

Свойства металлического материала можно изменять путем сплавления и термообработки. Так, предел текучести железа подобными методами можно повысить в 50 раз. Предел текучести чистого железа составляет примерно 40 МПа, тогда как предел текучести сталей, содержащих 0,5% углерода и несколько процентов хрома и никеля, после нагревания до 950° С и закалки может достигать 2000 МПа.

Временное сопротивление

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

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

Примеры значений временного сопротивления чистых металлов приведены в таблице:

| Аналитика бизнеса

Разведочный анализ

Файл данных

Таблица данных содержит информацию о 424 образцах стали различного состава:

  • процентное содержание в стали различных примесей (переменные 1-15: C-…-Ca),
  • величины временного сопротивления и предела текучести сталей данного состава.
| Аналитика бизнеса

На рисунке приведён фрагмент исходных данных.

Описательный анализ

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

| Аналитика бизнеса

Таблица описательных статистик:

Наиболее сильно в исследуемых образцах варьируется содержание Mn, Cr, Ni и Cu.

| Аналитика бизнеса

Введём новый показатель – суммарную процентную долю всех изучаемых нами 15 примесей в составе образца стали:

| Аналитика бизнеса

Тогда процентная доля железа (и, возможно, прочих примесей) в составе стали:

| Аналитика бизнеса

Диаграммы размаха процентного содержания всех примесей:

| Аналитика бизнеса

Медиана содержания примесей составляет 2,64% (50% образцов содержит более 2,64% примесей, 50% образцов – менее). Нижняя квартиль 2,60%, верхняя квартиль 2,70% (у 50% образцов содержание примесей в этом диапазоне, т.е. это наиболее типичный диапазон).

 

Если рассмотреть диаграмму рассеяния по переменным Предел текучести – Временное сопротивление, можно заметить, что эти величины сильно скоррелированы: коэффициент корреляции Пирсона равен r = 0,92 и является статистически значимым. Скоррелированность показателей может объясняться как реально существующей взаимосвязью между ними, так и особенностями планирования составов смесей при проведении эксперимента.

| Аналитика бизнеса

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

| Аналитика бизнеса

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

 

Для поиска нетипичных наблюдений по сочетанию пары количественных признаков (в нашем случае это Предел текучестии Временное сопротивление) в STATISTICA можно также использовать графический инструмент Bag Plot (2-мерный аналог диаграммы размаха):

| Аналитика бизнеса

Точки, соответствующие образцам с нетипичными для данной выборки сочетаниями параметров выделены на графике «крестиками».

Выбор стали с нужными свойствами без построения явной модели смеси

Если мы уже имеем собранные экспериментальные данные, нас может заинтересовать поиск стали с нужным нам сочетанием свойств Предел текучести – Временное сопротивление.

 

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

 

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

| Аналитика бизнеса

Напомним, что исходная таблица содержала информацию о 424 опытах (число всех опытов, в т.ч. повторяющихся). Результирующая таблица содержит данные о 137 комбинациях компонент смеси (число различных опытов). Число реплик (повторений опыта с одним и тем же составом смеси) составило от 2 до 16 раз. Чаще всего проводилось 2 опыта для каждого состава.

 

Фрагмент результирующей таблицы представлен ниже:

| Аналитика бизнеса

Заметим, что диаграмма рассеяния по усреднённым значениям показателей Предел текучести – Временное сопротивлениеуже не содержит резких отклонений экспериментальных точек от общей прямой. Скорее всего, эти отклонения объяснялись не необычным составом смеси, а погрешностью измерения показателей свойств стали:

| Аналитика бизнеса

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

| Аналитика бизнеса

Для предела текучести значения стандартного отклонения выше 34,65 являются нетипичными. Такие значения встречались в следующих опытах:

| Аналитика бизнеса

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

 

Для временного сопротивления нетипичными являются значения стандартного отклонения выше 24,04. Такие значения встречались в следующих опытах:

| Аналитика бизнеса

Построим диаграмму Вороного по переменным Временное сопротивление и Предел текучести:

| Аналитика бизнеса

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

 

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

 

Предположим, что мы хотим определить состав стали, обладающей временным сопротивлением 620 и пределом текучести 720.

 

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

 

Проведём на графике соответствующие прямые; точка пересечения прямых будет соответствовать стали с искомыми свойствами:

| Аналитика бизнеса

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

| Аналитика бизнеса

Это оказался образец 62, его состав указан в таблице, экспериментально измеренное временное сопротивление 724, предел текучести 621.

| Аналитика бизнеса

Планирование эксперимента

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

 

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

 

Как правило, перед построением моделей взаимодействия компонентов смеси (см. далее раздел Анализ эксперимента) вначале планируют эксперимент: составляют т.н. план, состоящий из оптимального (в плане числа измерений) списка составов смесей, для которых впоследствии производится замер характеристик сплава – в нашем случае, временного сопротивления и предела текучести.

| Аналитика бизнеса

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

| Аналитика бизнеса

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

| Аналитика бизнеса

В модуле Планирование эксперимента/Планы для поверхностей и смесей с ограничениями можно также задать дополнительные ограничения на экспериментальную область, например, ограничение вида Sn+Ca<0,025.

| Аналитика бизнеса

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

| Аналитика бизнеса

Полная таблица содержит 2904 комбинации 16 компонент смеси, необходимых для полноценного оценивания в дальнейшем всех коэффициентов модели смеси и их влияния на характеристики стали.

Ниже показана диаграмма рассеяния точек полученного плана на треугольнике с вершинами С – Mn – Si:

| Аналитика бизнеса

Анализ эксперимента

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

 

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

 

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

 

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

 

В качестве зависимых переменных в нашей задаче выступают Предел текучести и Временное сопротивление, а в качестве независимых факторов – значения долей элементов смеси (С, Cr, Mn … , Fe+…):

| Аналитика бизнеса

Линейная модель

Для анализа основных компонент, без взаимодействий, смеси нам необходимо выбрать линейный тип модели в модуле Анализ и планирование экспериментов:

| Аналитика бизнеса

Все необходимы результаты анализа эффектов находятся на вкладке Быстрый:

| Аналитика бизнеса

Нажав на кнопку Дисперсионный анализ, получим таблицы для пары характеристик стали:

| Аналитика бизнеса

Из таблицы видно, что полученные результаты являются статистически значимыми (p<<0,05), доля объясненной дисперсии равна 0,45 для предела текучести и 0,46 для временного сопротивления. Нажав на кнопку Оценки исходных компонент, получим численные оценки эффектов:

| Аналитика бизнеса
| Аналитика бизнеса

Визуализировать численные значения таблиц можно, построив диаграммы Парето (кнопка Карта Парето эффектов):

| Аналитика бизнеса

Из диаграммы Парето для предела текучести видно, что из основных компонент статистически значимым оказался эффект Мо, Cr, Fe+…, Mn, Ca, Al, Sn и P, остальные являются слабо значимыми. Значения коэффициентов регрессии указаны в таблице (см. выше): например, напротив Mn стоит цифра 145,3 – это означает, что в среднем при изменении на 1% значения Mn, предел текучести увеличивается на 145,3 единицы.

 

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

| Аналитика бизнеса

Из диаграммы Парето для временного сопротивления видно, что из основных компонент статистически значимым оказался эффект Fe+…, Mo, Cr, C, Mn, Al, Ca, Sn, P и V, остальные являются слабо значимыми. Значения коэффициентов регрессии указаны в таблице (см. выше): например, напротив Cr стоит цифра 151,8 – это означает, что в среднем при изменении на 1% значения Cr, временное сопротивление увеличивается на 151,8 единицы.

 

Коэффициенты перед Ca и Sn имеют отрицательный знак – т.е. из построенной модели следует, что увеличение концентрации этих элементов ведёт к уменьшению временного сопротивления.

| Аналитика бизнеса

Квадратичная модель

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

| Аналитика бизнеса

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

 

Например, в рассматриваемой нами задаче целесообразно игнорировать эффекты:

  • с Fe+… (т.к. это основная составляющая стали)
  •  с С, Nb, O – т.к. в имеющемся у нас плане эксперимента концентрации этих примесей могут принимать всего лишь 3 различных значения – а этого скорее всего будет недостаточно для точной оценки их вкладов.
| Аналитика бизнеса
| Аналитика бизнеса

Из таблиц видно, что полученные результаты являются статистически значимыми (p<<0,05), доля объясненной дисперсии равна 0,66 для предела текучести и 0,69 для временного сопротивления.

 

Нажав на кнопку Оценки исходных компонент, получим численные оценки эффектов и их взаимодействий (взаимодействия закодированы латинскими буквами, например BC означает взаимодействие Si (B) и Mn (C)):

| Аналитика бизнеса
| Аналитика бизнеса

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

 

Визуализировать численные значения таблиц можно, построив диаграммы Парето:

| Аналитика бизнеса

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

 

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

| Аналитика бизнеса

Разброс точек относительно подгоночной прямой характеризует точность модели.

 

Для оценки величины ошибок можно построить гистограмму остатков. Наиболее типичный диапазон ошибок – около 20 единиц временного сопротивления. Гистограмма остатков для временного сопротивления симметрична относительно нуля, что говорит о несмещённости модели:

| Аналитика бизнеса

Итак, в данном примере мы выявили основные компоненты и их взаимодействия, влияющие на характеристики стали, и получили значения их эффектов влияния.

Почитать еще

| Аналитика бизнеса

Выборка. Типы выборок

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

| Аналитика бизнеса

Обзор основных видов сегментации

Загрузить программу ВІ Демонстрации решений Аналитика бизнеса Оглавление Сегментация бренда Сегментация помогает принимать более эффективные

| Аналитика бизнеса

Метод Тагучи

Японский ученый Г. Тагучи в 1960 г. высказал мысль, что качество не может более рассматриваться

Несколько видео о наших продуктах

| Аналитика бизнеса
Проиграть видео
Презентация аналитической платформы Tibco Spotfire
| Аналитика бизнеса
Проиграть видео
Отличительные особенности Tibco Spotfire 10X
| Аналитика бизнеса
Проиграть видео
Как аналитика данных помогает менеджерам компании
2021-02-05T11:47:26+02:00