KINETIC, THERMOFLUCTUATION NATURE FRICTION SURFACES OF SOLID FRICTION DURING WEAR


Cite item

Full Text

Abstract

The article describes the basics of the kinetic theory of strength and the thermofluctuation mechanism of destruction of the friction surfaces of solids during wear. The relationship of strength with the mechanical and physical constants of a solid is analyzed. The kinetics of the state of materials of the surfaces of friction pairs in the process of wear is studied. The temperature-time dependence of changes in the physical, mechanical parameters of materials and the destruction of surfaces and near-surface layers of friction pairs is confirmed. The basic principles of the physical model of damage to friction surfaces are formulated and the activation parameters of the thermofluctuation model of their destruction are estimated. A calculation model of friction wear based on the thermofluctuation fracture concept is proposed.

Full Text

ВВЕДЕНИЕ В середине прошлого столетия получила развитие кинетическая концепция прочности основанная на новых представлениях о термофлуктуационной природе (механизме) разрушения твердых тел, основоположником которой является академик С.Н. Журков с сотрудниками [1, 2]. В соответствии с представлениями температурно-временнόй концепции прочности долговечность образца под нагрузкой при испытаниях на разрыв определится по следующей экспериментальной зависимости: t = τ0exp[(U0 - γσ)/kT], (1) откуда разрывное напряжение σр= (1/γ) (U0- kT ln(t/τ 0)), (2) а максимальное разрывное напряжение при температуре Т= 0 oK U0 /γ, (3) где t - время до разрушения (долговечность) образца под нагрузкой, с.; τ0=(10-12…10-13) с. - период колебаний атомов в твердом теле; U0 - энергия межатомной связи, Дж; σ - напряжение, Н/м2; T - термодинамическая температура, oK; k=1.38·10-23 Дж/oK - постоянная Больцмана; γ - параметр (подгоночный), имеющий размерность объема, м3; (U0 - γσ)=ΔU - энергия активации разрушения; σр - разрывное напряжение, Н/м2; - максимальное разрывное напряжение, Н/м2. АКТИВАЦИОННЫЕ ПАРАМЕТРЫ РАЗРУШЕНИЯ Экспериментальные зависимости (1,2,3) дали возможность исследователям [1.2] оценить численные значения и физический смысл активационных параметров разрушения γ, U0, τ0 при температурно-временных испытаниях на одноосное растяжение образцов. Численные значения и физический смысл коэффициента γ определены путем вычисления его в модели распада нагруженной атомной ангармонической цепочки, находящейся в тепловом движении [1]: χ, (4) откуда следует, что γ рассматривается как мера локальной перегрузки χ=Σ/σ=σn/σp= , (5) где χ - коэффициент локальной перегрузки в месте разрыва созданный концентратором напряжения; Σ и σ - соответственно, локальное и среднее напряжения; σn=E - предельная (теоретическая) прочность; Е - модуль Юнга, Н/м2; - относительное удлинение, при котором межатомная связь теряет устойчивость и разрывается; С - атомная (молярная) теплоемкость, Дж/oK ( для одноатомных кристаллов С≈3k) [1]; a - коэффициент термического линейного расширения, К-1. Далее комбинируя (3), (4) и (5), получим: U0/γ =U0/(Сχ); (6) где [1]. (7) Таким образом, полученные в [1, 2] выражения для параметра γ (4) и энергии механического разрушения (связи атомов в твердом теле) U0 (7) позволили установить связь прочности σр с механо-физическими константами твердого тела: модулем упругости Е, коэффициентом термического линейного расширения a и атомной (молярной) теплоемкостью С. Закономерности (1,2,3) и оцененные значения активационных параметров U0, τ0 и γ позволили интерпретировать механизм разрушения твердых тел с кинетической позиции как термофлуктуационный процесс распада межатомных связей, накопления разорванных межатомных связей, субмикроскопических трещин и генерации зародышевых микротрещин [1]. КИНЕТИКА СОСТОЯНИЯ МАТЕРИАЛОВ ПОВЕРХНОСТЕЙ ПАР ТРЕНИЯ ПРИ ИЗНАШИВАНИИ Подтверждением приведенной выше интерпретации кинетической концепции разрушения явились комплексные систематические исследования автора [3, 4, 5, 6, 7, 8, 9] кинетики накопления повреждений структуры и разрушения материалов контактирующих пар при изнашивании в условиях вибрационного трения (фреттинга) в малоподвижных и номинально неподвижных соединениях. Оборудование и средства технического оснащения при проведении исследований описаны в работах [3, 6, 7, 9]. Оценивалась кинетика изнашивания и состояния структуры материалов поверхностей и приповерхностных слоев образцов из титановых сплавов ВТ9, ВТ3-1, бронзы БРАЖМц10-3-1,5 при фретинге образцов титановых сплавов в одноименных парах и в парах с образцами из бронзы и жаропрочного сплава ЖС6К. Установлено, что независимо от природы материалов, условий испытаний и режимов нагружения временны́е зависимости линейного износа имеют ступенчатый, а характеристики изменения состояния материалов - периодический характер. На примере изнашивания сплава ВТ9 (рис.1а) при фретинге в одноименной паре без смазки при амплитуде осцилляции А=0,015 мм, частоте f=310Гц, статической нагрузке qст = 7МПа, динамической нагрузке qдин= 3МПа и температуре Т=293 oK можно выделить три характерных стадии износа: П - приработки, I - умеренного износа (латентный период накопления повреждений), II - интенсивного износа (разрущения). Стадии умернного I и интенсивного II износа соответствуют различным структурным состояниям (рис. 1б) кристалической решётки. Один цикл физического уширения β (интегральной характеристики микроискажений Δd/d кристаллической решетки материала, дробления блоков мозаик D и изменения плотности дислокаций ρ) рентгеновских линий 01.1 и 01.3a-Ti, полученных методом дифракции на дифрактометре УРС-50ИМ в медном Cu-Ka излучении, представляет собой последовательные стадии упрочнения, разупрочнения и разрушения поверхностного слоя, обусловленные развитием микроискажений решетки, дроблением блоков мозаик, увеличением плотности дислокаций, образованием и слиянием субмикро - и - микротрещин и отделением частиц разрушенного материала. Из совместного анализа кривых износа и физического уширения рентгеновских линий (рис. 1 а,б) следует, что на стадии приработки П происходит интенсивное нарастание износа и разупрочнения поверхностного слоя в связи с отделением диспергированных поверхностных структур, упрочненных при предварительной механической обработке, и обнажение нижележащих неупрочненных слоев материала. Ширина рентгеновской линии 01.3a - Ti уменьшается при этом с 60·10-4 до 41·10-4 рад. На стадии I умеренного износа физическое уширение резко возрастает, поверхности упрочняются, а износ стабилизируется. На этой стадии происходит накопление дефектов тонкой структуры (фрагментация блоков, рост микроискажений и плотности дислокаций) на поверхности и в зонах влияния по глубине [7, c. 353,354]. На стадии II интенсивного износа наблюдается разрушение, отделение частиц износа и связанное с этими процессами разупрочнение поверхностей, ширина ренгеновских линий вновь уменьшается. Кинетическая картина накопления повреждений кристаллической структуры и изнашивания сплавов ВТ3-1, ВТ9 в паре со сплавом ЖС6К и бронзой БРАЖМц10-3-1,5, бронзы БРАЖМц10-30-1,5 в паре со сплавом ВТ9 аналогична описанной на примере сплава ВТ9 (рис.1). То есть изменение условий контактирования, материалов контактных пар, параметров механических воздействий, рабочих сред приводит лишь к изменению уровней выходных параметров: амплитуды и периода циклов упрочнение - разрушение - разупрочнение и, следовательно, интенсивности износа рабочих поверхностей трения. Поверхностные слои разрушаются вследствие насыщения кристаллической структуры дефектами, упрочнения, развития хрупкости, потери сплошности и периодической разрядки напряженного состояния на поверхностях и в глубинных уровнях деформируемых объемов [7]. Каждой стадии цикла «упрочнение - разрушение - разупрочнение» соответствует вполне определенный уровень и циклический временнόй характер изменений комплекса элементарных процессов [7, c 353] на поверхностях трения: износа, физического уширения рентгеновских дифракционных линий β, микротвердости HМ, шероховатости RZ, относительной опорной площади tp, демпфирования (логарифмического декремента колебаний) δ, жесткости стыка К, коэффициента трения μ, фазовых превращений Vβ, указывающих на наличие в локальных объемах деформации термических флуктуаций энергии. Результаты рентгеноструктурного изучения кинетики дефектов тонкой кристаллической структуры сплава ВТ9 на стадиях упрочнения - разрушения в различных глубинных уровнях от поверхности трения при фреттинге в одноименной паре без смазки приведены в работе [7, c. 354, рис. 2]. Использовались рентгеновские излучения с различной проникающей способностью от антикатодов рентгеновских трубок Cr, Fe, Co, Cu, что позволило исследовать послойно состояние дефектной структуры деформированного поверхностного слоя без его разрушения на глубину до 15 мкм. Фиксировалась кинетика изменений микроискажений Δd/d кристаллической решетки, блоков мозаик D, рентгенографической плотности дислокаций ρ в результате совместной деформации в убывающем градиентном поле внешних эквивалентных напряжений σэкв на контакте и термофлуктуационных напряжений. Общая картина изменения состояния тонкой кристаллической структуры в слоях по глубине при различных уровнях напряжений и деформации показывает, что предварительно упрочненные слои, в результате роста микроискажений, генерации субмикро - и - микротрещин и фрагментации структуры на блоки, в результате которой плотность дислокаций достигает пределов насыщения для данного уровня напряжений, разупрочняются вновь, вследствие разрядки напряженного состояния, снижения прочности, развития хрупкости, нарушения сплошности и разрушения слоя на отдельные элементы дефектной структуры. Наиболее интенсивно описанные процессы протекают в слоях глубиной до 15 мкм, менее интенсивно в глубинных слоях, где и сохраняется исходная (технологическая) дефектная структура. Таким образом, в различных структурных уровнях по глубине от поверхности и на поверхности трения протекают интенсивные структурные и атомно-молекулярные перегруппировки: генерирование, взаимодействие и размножение дефектов связанное с термофлуктуационным разрушением межатомных связей в поле градиентно убывающих напряжений от внешних сил трения. ОБЪЕДИНЕННАЯ МОДЕЛЬ ПОВРЕЖДАЕМОСТИ И РАЗРУШЕНИЯ ПОВЕРХНОСТЕЙ ТРЕНИЯ На основе проведенного исследования кинетики деформационных процессов, накопления повреждений структуры и разрушения на поверхностях трения и в приповерхностных слоях сформулированы основные положения многоуровневой объединенной модели повреждаемости и разрушения, а также расчетной модели изнашивания при трении [3, 4, 5]. В основу предложенной физической модели повреждаемости и разрушения [4, рис.1] положены следующие основные принципы. 1. Контакт шероховатых поверхностей дискретный. Нарушенные фрикционные связи на вершинах микронеровностей при сдвиге поверхностей на величину элементарного пятна контакта «d» возникают вновь на вершинах новых микронеровностей вступивших в контакт, образуя элементарные площадки, в совокупности составляющие фактическую площадь контакта (ФПК) 2. К ФПК приложено эквивалентное напряжение обусловленное совместным действием нормальных и касательных напряжений на элементарных площадках контакта при относительном сдвиге поверхностей. Поверхностные контактные напряжения распространяются в глубинные слои, в соответствии с принципом Сен-Венана [10], создавая градиент убывающих по величине напряжений от внешних сил. Специфика контактного взаимодействия и деформации, обусловленная дискретностью шероховатого контакта и градиентом напряжений, создала условия наблюдаемой выше кинетики накопления повреждений структуры и разрушения на трущихся поверхностях и на различных структурных уровнях зон влияния по глубине. В деформируемых структурных уровнях периодически возникает высокая плотность элементарных дефектов, при взаимодействии которых генерируются и распространяются более сложные дефекты и повреждения, являющиеся локальными концентраторами механических и термических напряжений. 3. В глубинном уровне, где внешнее напряжение σ=0 [4, рис.1], возможные термофлуктуационные [1] разрывы единичных межатомных связей на технологических дефектах сразу же восстанавливаются вследствие отсутствия поверхностно - активной среды. 4. По мере продвижения дефектных слоев материала из глубинных уровней к поверхности медленно нарастающее внешнее напряжение способствует разрыву межатомных связей и препятствует их восстановлению, создавая направленность термофлуктуационного разрушения в поле внешних сил. Энергия активации термофлуктуационного разрыва межатомных связей снижается, а вероятность и скорость разрушения возрастают. 5. Разрыв связей происходит в области локальной перегрузки в виде локального активационного объема (где VA атомный объем) в вершине микротрещины, перемещающегося по образцу по мере продвижения фронта трещины. Освобождающаяся при разрыве связей энергия приводит к генерированию фононов определенных частот (где h=6.626·10-34 Дж/Гц - постоянная Планка, - частота флуктуаций атомов при данной температуре Т°К, Гц.), их взаимодействию и распространению [11]. 6. На поверхности и в приповерхностных слоях накапливается множество локальных объемов (структурных элементов) с различным уровнем прочности, размеры которых определяются величиной коэффициентов локальной перегрузки. Разрушение происходит в локальных объемах, где действующие напряжения достигают локального предела прочности наиболее слабого элемента структуры (хрупкое разрушение) или усредненного предела прочности , где - предел прочности, группы структурных элементов (вязкое разрушение). Разрушение завершается отделением частиц, размеры которых определяются размерами областей микронеоднородностей или элементов деформированной структуры. АКТИВАЦИОННЫЕ ПАРАМЕТРЫ РАЗРУШЕНИЯ ПОВЕРХНОСТЕЙ ТРЕНИЯ В твёрдых телах без микроповреждений в момент приложения нагрузки все межатомные связи нагружены равномерно, реализуется механизм одновременного разрыва связей. В этом случае коэффициент перегрузки цепей χ=1, локальный объём разрушения γ = χ = = , где - межатомное расстояние, - расстояние между минимумом и максимумом потенциальной энергии при разрыве, - локальный атомно-молекулярный объём [11, 12], а время ожидания флуктуации определится из уравнения(1), как t = τ0exp[(U0 - )/kT], (8) разрывное напряжение (2), σр= /-(kT/ )ln(t/τ 0 ), (9) предельная прочность (3), /. (10) В реальных телах с дефектами , γ = χ = , где - область локальной перегрузки в вершине микротрещины (трещины) или локальный нано - микрообъём разрушения, в котором ожидается флуктуация кинетической энергии. В этом случае время ожидания флуктуации t = τ0exp[(U0 - )/kT], (11) разрывное напряжение σр= / -(kT/ )ln(t/τ 0 ), (12) предельная прочность / (13) Кумулятивное накопление дефектов приводит к образованию групп структурных элементов с локальным активационным объёмом разрушения с низшим уровнем прочности . Локальный активационный объём физически и математически (11) связан с временем ожидания флуктуации t в этом объёме, что, как показано ниже, является основополагающим в механизме термофлуктуационного разрушения этого объема и может быть использовано в расчётах скорости изнашивания. Размерная идентификация расчётных активационных параметров разрушения при изнашивании в условиях фреттинга титанового сплава ВТ9 [6, 7] показала следующие пределы их рассеяния: математическое ожидание параметра то есть рассеяние коэффициента локальной перегрузки происходит в пределах одного порядка величин, что согласуется с данными источника [2]. Далее: 1,38·10-23)/(8,5·10-6)= =0,0487·10-17 Дж; =(3·1,38·10-23)/( 8,5·10-6·1,1·1011)= =4,4278·10-29 м3; 4,4278·10-29 м3, при χ=1; 4,4278·10-28 м3, при χ=10. Данные экспериментального электронографического исследования частиц износа [7] экстрагированных с угольной репликой, характеризующих исходное, не подвергнутое многократному перетиранию, состояние поверхностных структур и отделившихся от образца частиц (блоков мозаик D) показывают, что размеры областей когерентного рассеяния (ОКР) составляют, в основном, порядка (7-15) Å и очень близки к расчётным размерам локальных активационных объёмов структуры с дефектами: = (3,43·10-28 … 3,375·10-27) м3. ТЕРМОФЛУКТУАЦИОННАЯ МОДЕЛЬ РАЗРУШЕНИЯ ПОВЕРХНОСТЕЙ ТРЕНИЯ ПРИ ИЗНАШИВАНИИ В основу описания большинства современных моделей изнашивания [10, 13, 14, 15] положены выходные (конечные) объёмные макрохарактеристики процесса: например, число отделившихся с сопряжённой поверхности атомов «Z» на единичном пути трения; деформируемый объём «» равный объёму материала взаимно внедрившихся выступов шероховатости поверхностей, разрушаемый за определённое число циклов взаимодействий «n» [10]; объём удаляемого с поверхности трения материала за цикл изнашивания в результате разрушения определённого количества «λ» межатомных связей этого объёма с основой [13, 14]; деформируемый объём равный объёму взаимно внедрившихся субмикрошероховатостей [15]. Очевидно, что такие модели не могут привести к желаемым положительным расчётным результатам в силу неопределённости размеров и границ вовлечённых в контактную деформацию объёмов материалов и отсутствия конечных условий по времени или количеству циклов взаимодействий до их разрушения. Использование в расчётной модели изнашивания [5] локального активационного объёма разрушения «V0» структурного элемента или группы структурных элементов, сформированного в процессе деформации по естественному термофлуктуационному механизму и имеющего прямую математическую связь (11) с временем ожидания флуктуации t , исключает перечисленные выше неопределённости. Примем, что в установившемся режиме трения и изнашивания разрушение активационного объёма V0, в каком бы уровне по глубине поверхностного слоя он ни находился, в конечном итоге будет эквивалентно отделению с номинальной площади поверхности трения слоя материала толщиной , (14) где - номинальная площадь поверхности трения, м2; активационный объём, м3; - толщина отделяемого слоя материала, м. Скорость износа, с учётом (11) и (14), = (15) = /(τ0· exp[(U0 - )/kT]), где - скорость линейного износа, м/с; - геометрический фактор (число опорных точек или число условий связи относительного движения поверхностей в контактах: плоскость - плоскость, =3; по линии, =2; в точке, =1). Учитывая, что (16) с учётом (4) и (5), =; U0 , с учётом (3), U0 , получим расчётную модель скорости изнашивания поверхностей трения в широком температурном диапазоне в виде: (17) или с учётом выражения (16) для , . (18) Уравнение (17) удобно применять при использовании значений полученных экспериментально, уравнение (18) - при использовании расчётных значений. Эквивалентное напряжение при одновременном действии нормальных и тангенциальных напряжений на поверхностях трения определится из выражения: (19) где - эквивалентное напряжение, Н/м2; - касательное напряжение на поверхности трения, Н/м2; - фактическое давление в контакте, Н/м2; - коэффициент трения; - скорость скольжения, м/c; - плотность, кг/м3; G ≈ 0,38E - модуль сдвига, Н/м2; = 2πАf, м/c - скорость скольжения при возвратно-поступательном и вибрационном трени; А - амплитуда, м; f - частота осцилляций, Гц. При практических расчетах в условиях малых и средних нагрузок на контакт пластическая деформация возникает на выступах микронеровностей, материал «течет» в зазор, а среднее фактическое давление в контакте достигает = (1,1÷1,3). При высоких нагрузках внедрение микронеровностей увеличивается, пластическая деформация распространяется по всему контакту, фактическое давление растёт до ≈ 3[10], а эквивалентное напряжение достигает или превышает предел прочности , что приводит поверхностный слой в состояние низшей прочности и к ведущей роли «атермического» механизма разрушения, когда время ожидания термических флуктуаций становится больше времени действия нагрузки. При , энергия активации (15,16,17) = 0, работа внешней силы становится равной энергии связи атомов в твердом теле , (20) показатель экспоненты в уравнении (17) равным нулю, а экспоненциальный член = exp(-0) = 1. (21) Учитывая (20,21) уравнение (17) перепишется в виде: . (22) Наряду с атермическим механизмом разрушения при Т >0 oK действует термофлуктуационный, что следует из анализа предэкспоненциального члена уравнений (17,18), (22), поэтому понятие атермического механизма разрушения при температурах бльших абсолютного нуля является условным. Проведена широкая идентификация модели (22) при вибротрении и фреттинге в лабораторных условиях, при стендовых и эксплуатационных испытаниях [3, 6, 7, 9] сплавов ВТ9, ВТ8, ВТ3-1, ОТ4, 20Х13, ЖС6КП, 40ХНМА в различных температурных режимах. Расчётные данные удовлетворительно согласуются с экспериментальными. Рассеяние большинства экспериментальных данных вокруг расчётных значений не превышало ±23%, отдельные индивидуальные экспериментальные значения скорости изнашивания варьировали в диапазоне не более ±36%. ВЫВОДЫ Предложена кинетическая термофлуктуационная модель изнашивания при трении, прошедшая широкую идентификацию в условиях вибрационного трения и фреттинга. В модели используются механические и физические параметры и константы материалов при полном отсутствии подгоночных и эмпирических коэффициентов. Оценены локальные активационные объёмы разрушения, формируемые в процессе деформации по естественным термофлуктуационным механизмам. Для расчёта скорости изнашивания использована прямая фундаментальная математическая связь формируемых локальных активационных объёмов с временем ожидания флуктуаций кинетической энергии (разрушения). Предложенная модель изнашивания пополнит арсенал расчётных методов оценки износа при трении. Рис. 1. Кинетика износа (а) и физического уширения рентгеновских интерференционных линий (б) сплава ВТ9
×

About the authors

Anatoly Gavrilovich Kovshov

Samara State Technical University

Email: k.ntm@mail.ru
Candidate of Technical Sciences, Associate Professor of the Department «Technology of Mechanical Engineering, Machines and Tools»

References

  1. Журков С.Н. К вопросу о физической основе прочности // Физика твердого тела. 1980. Т. 22. Вып. 11. С. 3344-3349.
  2. Регель В.Р., Слуцкер А.Б., Томашевский В.Д. Кинетическая теория прочности твердых тел. М.: Наука, 1974. 302 с.
  3. Ковшов А.Г. Разработка триботехнических методов расчета и повышения фреттингостойкости деталей турбокомпрессоров из сплавов титана ВТ9 и ВТ3-1. Автореф. дисс.. канд. техн. наук. 05.02.04 - «Трение и износ в машинах». Киев: КИИГА, 1988. 20 с.
  4. Ковшов А.Г. Физическая модель разрушения поверхностей трения // Сб. трудов МНТК «Актуальные проблемы трибологии», июнь 2007, в 2-х томах. Том 2. М.: Машиностроение, 2007. С. 206-213.
  5. Ковшов А.Г. Расчетная модель изнашивания при трении /Сб. трудов МНТК «Актуальные проблемы трибологии», июнь 2007, в 2-х томах. Том 2. М.: Машиностроение, 2007. С. 214-220.
  6. Ковшов А.Г. Расчетно-экспериментальная идентификация активационных параметров модели изнашивания поверхностей трения // Известия Самарского научного центра РАН. 2016. Т. 18. № 4(2), 2016. С. 341-346.
  7. Ковшов А.Г. Рентгенографическая и электронографическая оценка структурных и активационных параметров модели изнашивания поверхностей трения // Известия Самарского научного центра РАН. 2016. Т. 18, № 4(2). С. 352-356.
  8. Ковшов А.Г. Динамическая модель контактирования при фреттинге бандажных палок лопаток компрессора низкого давления газотурбинного двигателя // Известия Самарского научного центра РАН. 2018. Т. 20. № 4(2)(84). С. 202-207.
  9. Ковшов А.Г. Оценка эффективности при фреттинге упрочняющего деформационного формоизменения поверхностей трения, совмещенного с противоизносными покрытиями // Известия Самарского научного центра РАН. 2018. Т. 20. № 4(2)(84), 2018. С. 248-254.
  10. Крагельский И.В., Добычин М.А., Комбалов В.С. Основы расчетов на трение и износ. М.: Машиностроение, 1977. - 526 с.
  11. Бартенев Г.М. Прочность и механизм разрушения полимеров. - М.: Химия. 1984. 280с.
  12. Бартенев Г.М. Флуктуационный (активационный) объем и его роль в расчетах прочностных характеристик полимеров // Физико-химическая механика материалов. 1985. Т. 21. №1. С. 3-6.
  13. Громаковский Д.Г. Разработка концепции модели изнашивания и склерометрического способа оценки кинетических параметров разрушения поверхностей трения // Известия Самарского научного центра РАН. 2011. Т. 13. № 4(42)(3), С. 699-703.
  14. Ибатуллин И.Д. Кинетика усталостной повреждаемости и разрушения поверхностных слоёв. / Самара: Изд-во СамГТУ, 2008.-310с.
  15. Суслов А.Г., Шец С.П, Прудников М.И. Термофлуктуационная модель изнашивания поверхностей трения твёрдых тел при граничной смазке // Трение и смазка в машинах и механизмах. 2008. № 10. С. 40-47.

Supplementary files

Supplementary Files
Action
1. JATS XML

Copyright (c) 2020 Kovshov A.G.

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies