ИНЖЕНЕРНЫЕ МЕТОДЫ ТЕХНОЛОГИЧЕСКОГО ОБЕСПЕЧЕНИЯ РЕГЛАМЕНТИРОВАННЫХ ПАРАМЕТРОВ ШЕРОХОВАТОСТИ ФУНКЦИОНАЛЬНЫХ ПОВЕРХНОСТЕЙ ДЕТАЛЕЙ МАШИН В ПРОЦЕССЕ МЕХАНИЧЕСКОЙ ОБРАБОТКИ
Аннотация и ключевые слова
Аннотация (русский):
Рассмотрены основы инженерных методов технологического обеспечения регламентированных параметров шероховатости функциональных поверхностей деталей машин в процессе механической обработки с возможностью заданного закона их изменения в пределах перехода. Показан рекомендуемый путь решения этой задачи, который включает анализ условий эксплуатации деталей и соответствующие им величины и законы изменения обеспечиваемых параметров; вопросы статистического моделирования процессов их формирования; условия технологического обеспечения равноизносности поверхностей трения скольжения с примерами технологических систем и результатов их применения с учетом технологической наследственности.

Ключевые слова:
технология, шероховатость, статистическое моделирование, качество поверхности, гибкие технологические системы, технологическая наследственность, условия эксплуатации поверхностей
Текст
Текст (PDF): Читать Скачать

При решении задач технологического обеспечения качества функциональных поверхностей деталей машин весьма важен учет  корреляционных связей между выбранными параметрами качества поверхности, в частности между параметрами шероховатости. Решение этого вопроса позволит обоснованно выбирать минимально необходимое число параметров для их технологического обеспечения. Анализ графов корреляционных связей [2, 3] выявил практически 100 %-ную корреляцию между всеми стандартизованными параметрами шероховатости. В этом плане, с целью выявления практически линейных связей следует рассматривать не статистически значимые корреляционные связи, а корреляционные связи с коэффициентом корреляции r, близким к единице.

Установлено [2, 3], что средняя степень вершин графа  для алмазного выглаживания (АВ) изменяется от 7,4 до 2,0; тогда как соответствующее изменение для предварительной обработки торцевым фрезерованием составляет от 7,44 до 4,44. Это естественно, так как из теории резания известно, что лезвийными методами обработки материалов обеспечивается высокая степень регулярности параметров микропрофиля, особенно высотных, которые можно достаточно точно рассчитать геометрическим путём.

Процесс деформирования микропрофиля поверхности методами ППД зависит как от параметров качества предварительно обработанной поверхности, так и от дополнительно вводимых факторов самого метода, в частности алмазного выглаживания, что обуславливает снижение средней степени вершин графов корреляционных связей при высоких значениях r. Это закономерно, так как процесс пластического деформирования является нелинейным [1].

Получены физико-статистические модели формирования параметров шероховатости плоских поверхностей деталей из чугуна (Ra, Rz, Rmax, Rp, Sm, r, b, n, D), которые корректны в следующих границах факторного пространства технологических систем (ТС) «торцевое фрезерование ПСТМ + АВ»: подача при фрезеровании Szфрез = 0,05…0,25 мм/зуб; скорость фрезерования vфрез = 311…622 м/мин; сила выглаживания Q = 50…150 Н; подача при АВ SzАВ = 0,025…0,16 мм/об; скорость АВ vАВ = 35…90 м/мин. При исследовании технологической управляемости параметрами качества АВ после шлифования периферией круга режимы шлифования выбирались для обеспечения шероховатости Raисх = 0,16…0,6 мкм.

Рассмотрение интервалов изменения параметров шероховатости поверхностей деталей из чугуна в ходе проведения экспериментов для разных методов обработки показывает достаточно широкий диапазон их варьирования алмазным выглаживанием после предварительного торцевого фрезерования. Это особенно характерно для параметров Ra, Rp, Rz, Rmax, Sm и DH. Предварительная обработка шлифованием может быть рекомендована для поверхностей деталей с повышенными требованиями к точности.

Надёжность технологического обеспечения параметров шероховатости в исследуемых ТС АВ определялась как вероятность выполнения задания для параметра Yi в симметричном d-интервале  (0 < d < 1) по результатам обработки данных машинного эксперимента.

Установлено, что различные параметры качества при одном и том же допустимом

отклонении d (%) от средней заданной величины обеспечиваются с различной надёжностью, которая колеблется при d = 20 % для высотных параметров от 0,53 для Rp до 0,72 для Ra. Для остальных исследуемых параметров – от 0,25 для r, D до 0,72 для n. Высотные параметры Ra, Rz, Rmax обеспечиваются с надёжностью P = 0,9 при d до 35 %, параметры Sm, b, n – при значении d от 30 до 60 %.

В основу работ по технологическому обеспечению качества и эксплуатационных свойств поверхностей деталей машин положены методы физико-статистического моделирования, которые базируются на использовании имитационных моделей процессов обработки [2, 6]. В ходе обработки поверхности детали происходит процесс её эволюции от исходной до технологической.

В общем случае технологический процесс может составлять q технологических систем, входными параметрами которых, например для ТСe, являются параметры качества поверхности после предшествующего этапа обработки (Yi(e–1)), управляющими – условия обработки (Xe1, Xe2, …, Xex), а выходными – полученные параметры качества поверхности (Yie).

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

Рис. 1

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

Для возвратно-поступательных пар трения скольжения, включающих широкий класс направляющих технологического оборудования и оснастки (от призматических до цилиндрических), характерно изменение нормальной нагрузки P и скорости относительного скольжения v, как функций, непрерывно изменяющихся в пространстве и времени. Подробный анализ изменения функций P и v в пространстве и времени дан в работах [2, 3],  на основе которого выявлены типовые законы изменения указанных факторов в пространственной области для пар трения типа направляющих скольжения (рис. 2).

Рис. 2

Наиболее общим законом изменения нагрузки P в пространственной области является трапецеидальный (случай P1, рис. 2), для которого характерны участки нагружения ax1, работы под максимальной нагрузкой Р0 (участок x1x2) и разгружения x2b. При этом зависимости изменения нагрузки в пространстве координаты поверхности базового трибоэлемента X могут описываться различными уравнениями для различных участков, например fp1(x), fp2(x) и fp3(x), аналитический вид которых определяется спецификой работы соединения.

Более простым является случай P2 (см. рис. 2), когда зависимость fp(x) является непрерывной дифференцируемой функцией, имеющей общее аналитическое описание. Случай P3 (см. рис. 2) – это зависимость нагрузки P в виде стационарного случайного процесса в функции координаты поверхности X.

При рассмотрении скорости v относительного скольжения трибоэлементов следует иметь в виду её абсолютное значение, так как специфика работы направляющих скольжения технологического оборудования и других машин определяет её знакопеременность. Наиболее общим также является трапецеидальный закон зависимости скорости мобильного трибоэлемента от координаты рассматриваемой точки базового трибоэлемента (случай V1, рис. 2). В этом случае также характерно наличие участка разгона аx1, участка x1x2 работы трибосистемы со скоростью, близкой к максимальной (возможны случайные флуктуации), и участка торможения x2b.

Случай V2 (см. рис. 2) характеризуется возможностью интерпретации зависимости скорости мобильного трибоэлемента |v1| от текущей координаты X рассматриваемой  точки базового трибоэлемента непрерывной дифференцируемой аналитической функцией.

Скорость | v 1| может представлять собой стационарную случайную функцию от координаты X (случай V3, рис. 2), что характерно при движении трибосистем на малых скоростях, когда возможны случаи торможения и ускорения, например, в связи с процессами схватывания.

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

На предложенной графовой модели формирования качества и эксплуатационных свойств поверхностей трибоэлементов (рис. 3) [4] представлены в виде вершин состояния функциональных поверхностей трибоэлемента S0S4 для различных уровней жизненного цикла (I…III) и дуги переходов Uij из i-й вершины в j-ю, характеризующиеся видом функций от пространственной координаты X (например, длина плоской направляющей трения скольжения) и дополнительными параметрами, характеризующими множества условий обработки R, параметров параметров качества поверхностного слоя (КПС) , условий эксплуатации  (нагрузка, скорость относительного скольжения, температура, динамические факторы и др.). Вершины графа S1 и S3 соответствуют изотропности, а вершины S2 и S4 – анизотропности качества (S1, S2) и параметров эксплуатационных свойств (S3, S4) функциональных поверхностей трибоэлементов.

Рис. 3

Как итог технологического обеспечения эксплуатационных свойств поверхностей деталей машин возможны две вершины графовой модели:

· S3 – параметры эксплуатационных свойств не являются функцией от X и постоянны;

· S4 – параметры эксплуатационных свойств являются функцией от X и переменны.

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

S3 = (U01 Ù U13) Ú (U02 Ù U23).                                                                                            (1)

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

S0   S1    S3  .                                                                                           (2)

Однако на практике имеют место ситуации, когда условия эксплуатации  соединения трибоэлементов являются нестационарными. В этом случае функциональные поверхности со стационарными параметрами качества  = const, под воздействием нестационарных эксплуатационных нагрузок  = var характеризуются нестационарным вектором показателей эксплуатационных свойств , что соответствует состоянию S4 на графовой модели (см. рис. 3), которого можно достигнуть по пути (3):

S0   S1    S4  .                                                                                           (3)

 

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

S0   S2    S3  .                                                                                           (4)

 

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

1. Определение пространственно-временных закономерностей эксплуатационных факторов, действующих на сопряжение (Pi = fi1(X), Pi = fi2(t)), характеризующих нестационарность условий эксплуатации.

2. Выявление соответствующих закономерностей изменения значений параметров КПС по поверхности трибоэлемента в функции её координаты X(Ki = fi(X)), способных компенсировать нестационарность условий эксплуатации и обеспечить состояние поверхности S3, т.е. её изотропность по параметрам эксплуатационных свойств.

3. Формирование гаммы эффективных технологических методов обработки, позволяющих реализовать требуемые закономерности изменения значений параметров качества (наклёп, остаточные напряжения, шероховатость и др.) в функции координат обрабатываемой поверхности (Ki = fi(X)), т.е. реализуется путь U02, решающий задачу достижения вершины S3 (см. рис. 3), и, следовательно, задачу достижения изотропности поверхности по параметрам эксплуатационных свойств при нестационарных условиях эксплуатации.

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

1) статическая (Р1 = Р0 = const);

2) динамическая Р2, которая в частном случае может иметь амплитуду ap и период Т;          3) Р3 – нагрузка в виде стационарной случайной функции с математическим ожиданием M{Р3} = Р0 = const;

4) Р4 = F(X) – закономерно изменяющаяся (имеющая экстремум или монотонная);

5) Р5 – нагрузка в виде нестационарной случайной функции, математическое ожидание которой (в данном случае M{Р5} = Р4) зависит от X.

Вид нагрузки, передаваемый  на функциональную поверхность  детали, не определяет однозначно условия эксплуатации соединения. Одним из важных факторов, определяющих эпюру износа h поверхности, является функция распределения f(x) координаты X положения мобильного (подвижного) трибоэлемента. Величина износа h поверхности  будет равно распределённой при нагрузке Р1 и равномерном законе распределения f1(x). Имеется в виду, что параметры КПС распределены изотропно по поверхности. При любой стационарной нагрузке (Р1, Р2, Р3), но при распределениях fi(x), отличных от равномерного f1(x), величина износа h поверхности П будет распределена неравномерно. Это приводит в общем случае к потере точности и выходу из строя оборудования, особенно технологического.

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

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

1) технологическая гибкость 1-го рода (ТГС-1) – возможность выбора методов обработки поверхностей и группы управляемых факторов каждого из них, постоянных по величине в пределах соответствующего перехода и обеспечивающих изотропность поверхности по регламентируемым значениям параметров КПС или ЭС с требуемой надёжностью;

2) технологическая гибкость 2-го рода (ТГС-2) – возможность выбора методов обработки поверхностей и группы управляемых факторов, одним или несколькими из которых можно программно управлять в пределах перехода по закону, обеспечивающему заданную анизотропность распределения параметров КПС или ЭС по обрабатываемой поверхности с требуемой надёжностью. Таким образом, технологическая гибкость систем обработки является одним из критериев повышения качества обрабатываемых поверхностей.

Технологическая гибкость 2-го рода присуща системам ЧПУ, где имеется возможность изменять условия обработки (факторы) программным способом в пределах перехода или сложным технологическим  системам, в которых одна или несколько подсистем оснащены оборудованием с ЧПУ.

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

Понятие степеней свободы системы с технологической гибкостью 2-го рода можно пояснить следующим образом (рис. 4) [2, 3 и др.].

Рис. 4

При торцевом фрезеровании на станке с ЧПУ плоской функциональной поверхности обрабатываемой заготовки 1 (рис. 4, а) в пределах перехода имеется возможность управлять лишь одним параметром подачи S = Х1 при условии необходимости сохранения размера h. Следовательно, это система ТГС-2 с одной степенью свободы.

Примером системы обработки ТГС-2 с двумя степенями свободы может служить ТС накатывания (алмазного выглаживания) с ЧПУ типа CNC, которая обеспечивает возможность программного управления в пределах перехода двумя факторами обработки: величиной подачи S = Х1 и силой Q = Х2 (рис. 4, б).

Электромеханическая обработка (рис. 4, в) является примером эффективной технологической системы ТГС-2, обладающей тремя степенями свободы: величиной подачи

S = Х1, силой обработки Q = Х2 и силой тока I = X3, проходящего через площадку контакта между индентором и обрабатываемой поверхностью при наличии устройства 3 программного управления силой тока.

Очевидно, что важным показателем процесса обработки является степень его технологической гибкости. Если взять режимы обработки для точения (v, S, t) и для ППД (Q, S, v), то управляемых факторов в каждом случае ‒ по три, но технологическая гибкость процесса ППД значительно выше: в этом случае всеми параметрами КПС (макро- и микрогеометрия, физико-механические свойства и др.) можно управлять в гораздо более широких пределах и более надёжно, чем при точении. Применение той или иной системы зачастую способствует формированию совершенно новых, не свойственных основному материалу детали, свойств.

Примером реализации системы повышенной гибкости является система обработки ППД программным способом (ППДПС) «торцевое фрезерование композитом 10 – ППД алмазное выглаживание (АВ) или накатывание шариком (НШ)» на фрезерном станке с ЧПУ, позволяющая программно управлять параметрами  шероховатости Riпр предварительно обработанной поверхности, силой Q воздействия индентора на поверхность при обработке инструментами упругого действия программированным перемещением по соответствующей координате и величиной подачи Ѕ инструмента. Это способствует формированию  закономерно изменяемых параметров качества обрабатываемых поверхностей по заранее заданному закону [1, 3, 5].

Исследования таких ТС с ЧПУ по обеспечению заданных параметров качества обработки позволили получить следующие зависимости:

Ri = ,                                                                                                      (5)

где Ri – параметр качества (Ra, , Sm и др.); b0ib3i – коэффициенты.

Если Ri зад является функцией от длины обрабатываемой поверхности и для каждого участка хi  должна составлять определенную величину Ri(xi), то это достигается, например, за счёт изменения в процессе обработки ППД силы Q, определяемой из (5):

.                                                               (6)

Возможно управление параметром Ri(xi) за счёт технологической наследственности. В этом случае шероховатость предварительно обработанной поверхности на i-м участке должна составлять:

.                                                            (7)

Формулы (6) и (7) представляют собой законы для обеспечения программного управления параметром Ri зад (заданное значение) за счёт изменения силы воздействия индентора на обрабатываемую поверхность (Raпр, S = const, Q = var) и за счёт использования фактора технологической наследственности (Ri пр = var, Q, S = const).

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

В этом случае (рис. 5) типично наличие участков разбега X1, стационарного воздействия внешних факторов X2 и торможения X3. На этих участках соответственно должна изменяться и сила ППД Q при финишной обработке. При разработке и реализации программы, обеспечивающей трапецеидальную закономерность управления силой, а значит и качеством обработки при ППД, исходными данными наряду с параметрами X1, X2, X3 являются величины минимальной и максимальной силы (Qmin, Qmax), жёсткости пружины инструмента ППД, величина подачи S и скорость обработки v.

Рис. 5

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

Рис. 6

С целью количественной оценки влияния технологической наследственности на формирование параметров качества в процессе обработки целесообразно применить метод имитационного моделирования. При этом логично использовать модели Кобба-Дугласа, так как расчёты по ним не дают отрицательных значений параметров, т.е. не нарушается физическая картина их формирования [2, 6]. Таким образом, для параметра качества Yie, формируемого на этапе эволюции поверхности в системе ТСe, имитационная модель в общем случае имеет вид:

.                                                                     (8)

Здесь Yi0 – среднее экспериментальное значение i-го параметра качества, которое имеет соответствующую размерность и определяется по результатам активного эксперимента; Xem – параметр, численно равный значению m-го фактора обработки в ТСe; bem – показатель степени (1), соответствующий фактору обработки Xem.

Для  модели (8) соблюдается условие (9):

                                     m + x + y = k,                                                                  (9)

где k – общее число учтённых факторов техпроцесса в плане активного эксперимента.

Исходя из этого,  в общем случае имеем:

Yi = Yi0 + Yi1 + … + Yie + … + Yiq.                                                                                     (10)

; …; ; …; ,                        (11)

где Yi1, …, Yiq – это вклад подсистем технологического процесса ТС1, …, ТСq в формирование окончательного значения параметров качества Yi.

Предполагается, что величина Yi1, …, Yiq должна иметь размерность Yi.

Зависимость (8) можно представить в следующем виде:

Yi = Yi0 × K1 ×× Ke ×× Kq,                                                                                            (12)

здесь K1, …, Kq – коэффициенты влияния соответствующих ТС:

; …; ; …; .                                                 (13)

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

Kq = kq1 × kq2 ×× kqy,                                                                                                       (14)

где y – число управляющих факторов технологической подсистемы ТСq.

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

Графическое представление влияния технологической наследственности на формирование параметров микропрофиля (рис. 7) показывает, что оно различно для каждого параметра и для разных точек области факторного пространства.

Рис. 7

На основе проведенных исследований сформирована база данных значений коэффициентов влияния технологической наследственности на формирование параметров качества поверхности в ТС «торцевое фрезерование ПСТМ – АВ» для различных областей факторного пространства.

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

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

Список литературы

1. Инженерия поверхности деталей / Колл. авт.; под ред. А.Г. Суслова - М.: Машиностроение, 2008. - 320 с.

2. Нагоркин, М.Н. Параметрическая надежность тех-нологических систем чистовой и отделочно-упрочняющей обработки поверхностей деталей машин инструментами из сверхтвердых синтетических материалов: Монография / под ред. А.В. Киричека. - М.: Издательский дом «Спектр», 2017. - 304 с.

3. Нагоркин, М.Н., Федоров, В.П., Нагоркина, В.В., Ковалева, Е.В. Параметрическая надёжность технологи-ческих систем лезвийной и упрочняющей обработки инструментами из синтетических сверхтвёрдых материалов по геометрическим параметрам качества и триботехническим характеристикам поверхностей деталей // Прогрессивные машиностроительные технологии, оборудование и инструменты. Том VII. Коллективная монография / под ред. А.В. Киричека. - М.: Издательский дом «Спектр», 2016. - С. 506-688.

4. Нагоркин, М.Н., Федоров, В.П., Тотай, А.В. Усло-вия эксплуатации функциональных поверхностей трибо-элементов в типовых соединениях трения скольжения // Известия Юго-Западного государственного университета. - 2014. - № 6 (57). - С. 15-26.

5. Суслов, А.Г. Федонин О.Н., Польский Е.А. Науко-емкая технология повышения качества сборочных единиц машин на этапах жизненного цикла // Наукоемкие техно-логии в машиностроении. - 2016. - №5 (59). - С 34-42.

6. Суслов, А.Г., Федоров, В.П., Нагоркин, М.Н., Пы-риков, И.Л. Комплексный подход к экспериментальным исследованиям технологических систем металлообработки по обеспечению параметров качества и эксплуатационных свойств поверхностей деталей машин // Наукоемкие технологии в машиностроении. - 2018. - № 10(88). - С. 3-13.

7. Технологическое обеспечение и повышение экс-плуатационных свойств деталей и их соединений / А.Г. Суслов, В.П. Федоров, О.А. Горленко [и др.]; под общ. ред. А.Г. Суслова. - М.: Машиностроение. 2006. - 447 с.

8. Фёдоров, В.П., Нагоркин, М.Н., Ковалёва, Е.В., Чмыхов, Д.В. Адаптация поверхностей трибоэлементов к нестационарным условиям эксплуатации обработкой ППД программным способом // Упрочняющие технологии и покрытия. - 2007. - № 10. - С. 9-15.

Войти или Создать
* Забыли пароль?