PLANNING OF EXPERIMENTS TO INVESTIGATE THE DRIVELINE ON NEEDLE BEARINGS


如何引用文章

全文:

详细

The article presents the study on the driveline experimental design, basic formulas, defining function selected: the torque of the drive motor, the experimental factors: braking torque braking device, engine speed, driveline length, bend angle driveline. A brief description of the equipment, which is a stand, comprising a motor technology transmission driveline, the loading apparatus which are connected in series, and the loading device is hydraulic, comprising: a hydraulic pump, a choke, a filter, a pressure gauge, safety valve, tank, heat exchanger. The design of the proposed stand allows you to explore the driveline of different sizes at different angles of knee driveline. The stand is characterized by simplicity of design, a large range of parameters of loading and ease of setup thanks to the movable frame and vernier, calibrated in units of angle fracture driveline. The stand has a number of advantages: the test driveline while avoiding overheating of the braking system by cooling the working fluid, the improvement of operating conditions by providing a smooth regulation of the braking torque. To control the set parameters testing regimes instrumentation: strain gauges as gauges torque and braking moments, as a recording device using a dual-channel oscilloscope was used. The strain gauges were calibrated in torque and braking moments units empirically using technological devices. The task of the experiment is to find the optimal conditions driveline operation, namely the determination of the maximum efficiency of the driveline depending on the selected length and angle fracture driveline. The above method of design of experiments is designed for experimental studies of various drivelines with various types of bearings. Getting a polynomial regression model, the objective function - torque on the selected constructive-regime parametres of driveline allow further use of search methods the minimum values of the objective function, and, consequently, the efficiency of the transmission and corresponding to these values of constructive-regime parameters, allowing in the future to design phase of such programs to assign them reasonably, thus reducing design time frame while improving the quality of the design decisions, and to use significantly save energy while improving the transmission resource.

全文:

Введение. Подшипниковые узлы являются важнейшими структурными элементами машин и составляют основную часть узлов трения, используемых в том числе и в аэрокосмической отрасли. Отказ техники зачастую происходит из-за отказа подшипниковых узлов, которые ограничивают долговечность машин. При достаточно качественном изготовлении деталей подшипниковых узлов, например игольчатого подшипника качения, характеристики карданной передачи могут оказаться неудовлетворительными, и произойдет внезапный отказ. Практика показала, что по времени действия наиболее характерным режимом является работа шарнира при длительно действующем среднем крутящем моменте с наложением переменного момента, имеющего определенную частоту и амплитуду колебаний [1-5]. Поэтому для стендовых испытаний шарниров карданной передачи необходимо выбрать из всего многообразия эксплуатационных режимов те режимы, которые считаются наиболее характерными как по виду нагружения, так и по длительности их действия. Подшипниковые узлы, в том числе и работающие в аэрокосмической отрасли, выходят из строя в основном из-за абразивного изнашивания, связанного с попаданием грязи и пыли. Наиболее распространенный критерий отказа подшипниковых узлов общего применения, работающих в автомобилях, тракторах, насосах, редукторах, станках, подъемно-транспортном оборудовании, - усталостное разрушение. В то же время для подшипниковых узлов специального применения важны и другие характеристики, такие как жесткость, уровень и спектр вибрации, момент сопротивления вращению, долговечность и др. Достаточная долговечность отдельных деталей подшипниковых узлов, подтвержденная стендовыми испытаниями, не гарантирует достаточной долговечности всего узла. Последнее обстоятельство связано с тем, что нагрузки, действующие в узле, а также реальная температура могут существенно отличаться от стендовых. Кроме того, сборка и монтаж меняют зазоры, натяги и форму рабочих поверхностей подшипниковых узлов. Особенно отчетливо противоречие между качеством подшипниковых узлов и собственно подшипником, которое проявляется в узлах именно с подшипниками качения [6]. Исследованиями нагруженности трансмиссий, и в частности карданных передач автомобилей, тракторов, занимались З. Ш. Блох, Е. А. Чудаков, М. И. Лысов, И. С. Цитович, Р. В. Кугель, П. П. Лунин, С. А. Лапшин, Ю. С. Стефанович, М. С. Эйдинов, С. Н. Иванов и др. Описание оборудования. Для лабораторных испытаний будет использоваться стенд, позволяющий исследовать карданные передачи с разным углом излома при неизменной длине карданного шарнира. Конструкция стенда представлена на рис. 1, 2. Стенд состоит из электродвигателя 1, вал которого соединен с технологической передачей 2, которая соединена с испытываемой карданной передачей 3 выходным валом, который подсоединен к раздаточному редуктору 4 с гидронасосом 5, входной канал которого подсоединен в гидравлический бак 6, а выходной канал гидронасоса присоединен к регулятору нагрузки, выполненному в виде дросселя 7. Между дросселем и гидронасосом установлен манометр 8, отградуированный в единицах тормозного момента, и предохранительный клапан 9 для выпуска избыточного давления в гидравлический бак. Выходной канал дросселя подсоединен к гидравлическому баку, между которыми установлен теплообменник 10, входной канал которого подсоединен к гидравлическому баку. Электродвигатель 1 и технологическая передача 2 установлены на продольной подвижной раме, с перемещением которой меняется длина карданной передачи 3. Раздаточный редуктор 4 вместе с гидронасосом 5 установлены на поперечной подвижной раме, с помощью которой можно менять угол излома карданной передачи 3 [7-10]. Фрагмент7 Рис. 1. Стенд для исследования карданных передач Фрагмент new Рис. 2. Гидрокинематическая схема стенда для исследования карданных передач Величина угла излома карданной передачи изменяется перемещением поперечной подвижной рамы, имеющей нониус, с помощью которого выставляется угол излома карданной передачи для соответствующей длины карданной передачи. Конструкция предложенного стенда позволяет исследовать карданные передачи разных типоразмеров при различных углах излома карданной передачи. Стенд отличается простотой конструкции, большим диапазоном параметров нагружения и легкостью настройки благодаря подвижной раме и нониусу, отградуированному в единицах угла излома карданной передачи. Стенд имеет ряд преимуществ: испытание карданных передач при исключении перегрева системы торможения за счет охлаждения рабочей жидкости, улучшение условий эксплуатации за счет обеспечения плавности регулирования тормозного момента. Для контроля заданных параметров режимов испытаний использовалась измерительная аппаратура. В качестве измерителей крутящего и тормозного моментов использовались тензорезисторы марки 2ФКРВ-5-100ГБ, ГОСТ 21616-91 (техническая характеристика представлена в табл. 1), подключенные по мостовой схеме с температурной компенсацией (рис. 3). Исходя из деформации тензорезисторов, была выбрана электрическая схема тензомоста, в которой R1 и R4 - добавочное сопротивление (R1 + R4 = const), R2 и R3 - тензорезисторы, U - напряжение на выходе тензомоста, UВЫХ - напряжение питания тензомоста. Для компенсации температурных воздействий в схеме моста использовались добавочные сопротивления R1 и R4. Схема расположения тензорезисторов на карданном валу показана на рис. 4 [11-13]. Зависимость выходного сигнала от изменения сопротивления для данного вида деформации определялась по следующей формуле: (1) где U, UВЫХ - напряжения на входе и выходе тензомоста, В; ΔR - изменение сопротивления тензорезисторов, Ом; R - суммарное сопротивление тензомоста, Ом. Таблица 1 Техническая характеристика тензорезистора Обозначение 2ФКРВ-5-100ГБ Сопротивление, Ом 90,00-94,99 Чувствительность S при температуре +20 ± 5 ºС 2,08 Относительная погрешность П при температуре +20 ± 5 ºС, % 0,1 Максимальное относительное удлинение ε, мкм/м Не более ±3000 Максимальная сила тока I, mА Не более 54 Фрагмент2 Рис. 3. Электрическая схема тензомоста Фрагмент672 Рис. 4. Схема расположения тензорезисторов на карданном валу Двухэлементные тензорезисторы применяются в оболочках вращения (рис. 5), работающих под давлением, если на поверхности детали имеет место чистый сдвиг, как например, кручение валов, где на поверхности вала измеряются главные деформации ε1 = - ε2, которые действуют под углами ±45° к его продольной оси, и подсчитываются главные напряжения: (2) (3) где µ - коэффициент Пуансона; Е - модуль упругости; ε - величина деформации. При чистом сдвиге касательные напряжения численно равны главным напряжениям τ = σ1 = -σ2. Наибольшее значение относительной деформации используемых тензорезисторов ±3000 мкм/м, в отношении которой завод-изготовитель гарантирует надежную работу тензорезисторов и строгую линейность их характеристики преобразования. Закрепление тензорезисторов осуществлялось с помощью «суперклея», позволяющего приклеивать тензорезисторы без дополнительной температурной обработки. Правильно наклеенный тензорезистор обладает практически линейной характеристикой и обеспечивает достаточно высокую точность измерений, если деформация не превышает уровня ±0,3 %. Тем не менее, имеют место небольшие отклонения от линейности при циклическом деформировании, причем кривая разгрузки располагается ниже кривой нагружения, образуя петлю гистерезиса. Поэтому замер крутящего и тормозного моментов выполнялся с учетом циклического деформирования тензорезисторов, характеристика которого показана на рис. 6 в условном масштабе. В качестве регистрирующего устройства использовали двухканальный осциллограф HANTEK 6022BE с USB-выходом, подключенным к персональному компьютеру. Схема тензометрирования и аппаратура приведена на рис. 7, в которой через регистрирующее устройство РД сигнал от тензорезисторов поступает на персональный компьютер ПК, и с помощью программного продукта LabVIEW фиксируются значения крутящего и тормозного моментов. Тарировка тензорезисторов осуществляется опытным путем. Величина крутящего момента определяется с помощью приспособления, представляющего собой рычаг с динамометром, на котором замеряется величина крутящего момента, а величина тормозного момента тарируется по величине крутящего момента, но в единицах давления. Рис. 5. Схема чистого сдвига Рис. 6. Характеристика циклического деформирования тензорезистора (масштаб условный) Тензометрия Рис. 7. Схема тензометрирования Постановка задачи исследования. В процессе эксплуатации на карданную передачу воздействуют многочисленные возмущающие силы, переменные нагрузки разной частоты [11; 14-18], поэтому задача эксперимента состоит в поиске оптимальных условий функционирования карданной передачи, а именно, определении максимального коэффициента полезного действия карданной передачи в зависимости от конструктивных параметров: выбранной длины и угла излома карданной передачи, а также режимных параметров - нагрузочного (полезного) момента и угловой скорости вращения. Карданные валы работают на кручение, для изготовления карданных валов наиболее часто применяются стержни круглого или кольцевого сечений (круглые сплошные или полые валы). Деформацию таких валов при кручении можно характеризовать как взаимный поворот поперечных сечений (см. рис. 4). При этом поперечные сечения остаются плоскими, а расстояние между ними не изменяется. Угол поворота сечений является переменной величиной и называется углом закручивания φ(х). Характеристикой степени закручивания вала по длине является относительный угол, который определяется по формуле (4) Различие углов поворота поперечных сечений вызывает деформацию сдвига на поверхности вала и внутри него. При этом в поперечных сечениях круглых сплошных и полых валов действую только касательные напряжения, перпендикулярные к радиусам и определяемые по формуле (5) где Jp - полярный момент инерции, м4. Абсолютный угол закручивания определяется по формуле (6) где l - длина участка закручивания, м; G - модуль сдвига, Па; Jp - полярный момент инерции, м4. Согласно проделанному расчёту, при нагрузке на карданный вал величины крутящего момента на четырех режимах работы, величина закручивания (см. рис. 4), величина деформации сведены в сводную таблицу (табл. 2). Согласно проведенным теоретическим исследованиям (табл. 2) на четырех режимах работы при изменении крутящего момента меняется угол закручивания карданного вала, причём величина закручивания на каждом из участков различна, это связано с характером нагрузок и значениями конструктивных параметров карданной передачи, в том числе и поперечным сечением элементов конструкции, а также с наличием диссипации энергии в процессе упругих крутильных деформаций и нагревом передачи. При изменении длины карданной передачи меняются длины участков карданного вала согласно формуле (6), а следовательно, и угол закручивания, это значит, что меняется общая картина нагружения, что может вызвать изменение коэффициента полезного действия карданной передачи. Поэтому ожидается снижение коэффициента полезного действия при увеличении нагрузочного момента, а также при увеличении значений конструктивных параметров карданной передачи. Расчёт карданного вала производился в программмном продукте APM WinMachine (APM Shaft). Была разработана модель (рис. 8), на которой карданный вал был разбит на сегменты в соответствии с заданными значениями конструктивных параметров, с помощью которой графическим способом были получены величины углов закручивания карданного вала при заданных нагрузках. Предполагается провести испытания карданных передач на четырех режимах работы при разных значениях крутящего и тормозного моментов с изменением конструктивных параметров карданной передачи в каждом опыте, что позволит наиболее качественно оценить работу карданной передачи. Длительность испытаний выбираем, исходя из результатов проведенных исследований [11], что составляет 0,2·104 циклов при частоте вращения 1455, 1003, 644, 416 об/мин, что соответствует 30, 43, 67, 103 с работы стенда соответственно. На основании этих данных были выбраны режимы нагружения, приведенные в табл. 3. Таблица 2 Сводная таблица значений деформаций тензорезисторов Номер режима работы Крутящий момент, Н·м Угол закручивания тензорезисторов, град Величина закручивания тензорезисторов, мм 1 2 1 2 1 передача 17 2°43'39'' 1°4'4'' 0,45 0,28 2 передача 27 3°43'31'' 2°41'26'' 0,71 0,44 3 передача 42 4°12'37'' 3°37'14'' 1,10 0,69 4 передача 65 7°33'27'' 4°6'34'' 1,72 1,07 Рис. 8. Модель карданного вала в программном продукте APM WinMachine (APM Shaft) с сегментами для расчёта параметров сечений и местами установки и ориентации тензорезисторов марки 2ФКРВ-5-100ГБ, ГОСТ 21616-91 Таблица 3 Таблица режимных параметров для стендовых испытаний Параметры Режимы 1 2 3 4 n, об/мин 1455 1003 644 416 Мкр, Н×м 37 54 84 130 Мт, Н×м 17 27 42 65 t, с 30 43 67 103 g, град 6 6 6 6 l, мм 493 493 493 493 Таблица 4 Уровни факторов и интервалы варьирования Факторы Уровни факторов Интервал варьирования 1 2 3 4 Тормозной момент гидравлического устройства торможения Мт, Н·м 17 27 42 65 2 Частота вращения двигателя, об/мин 1455 1003 644 416 50 Длина карданной передачи l, мм 543 517 505 493 25 Угол излома карданной передачи g, град 6 6 6 6 1 Таблица 5 Формат записи данных эксперимента № опыта Мкр, Н·м Мт, Н·м n, об/мин g, град l, мм 1 2 3 4 5 6 7 8 36,900 37,000 37,200 37,000 36,900 35,700 36,300 36,600 16,400 17,000 17,000 16,700 17,000 16,700 16,900 17,600 1400,000 1404,000 1432,000 1421,000 1445,000 1445,000 1445,000 1405,000 6,100 6,200 6,000 6,500 6,300 6,300 6,300 6,200 493,000 502,000 509,000 515,000 521,000 530,000 537,000 543,000 При постановке и планировании эксперимента использовался план полного факторного эксперимента N = 44, факторами которого являются тормозной момент устройства торможения (Мт, Н×м), частота вращения двигателя (n, об/мин), длина карданной передачи (l, мм), угол излома карданной передачи (g, град). В качестве параметра оптимизации Y выбираем крутящий момент приводного электродвигателя (Мкр, Н×м). Тормозной момент отградуирован в единицах давления гидравлической системы, а крутящий момент определяется по формуле (7) где N - мощность привода, кВт; - угловая скорость, с-1, определяемая по формуле (8) где n - число оборотов электродвигателя в минуту, об/мин. При работе в карданной передаче [16-17; 19-21] возникает упругая деформация кручения, способствующая возникновению и росту моментов сил упругости, в том числе и за счёт погрешностей изготовления отдельных деталей конструкции карданной передачи, температуры, разностей жесткости отдельных элементов конструкции карданной передачи. В табл. 4 представлены уровни факторов и интервалы варьирования. Выполнив запланированные опыты и определив с помощью этой матрицы численные значения коэффициентов, получим уравнение линейной математической модели в виде квадратичного полинома: (9) Для определения чувствительности факторов и проверки выбранных значений их уровней был проведен предварительный эксперимент с 4 факторами. Формат записи данных эксперимента приведен в табл. 5. Обработку полученных результатов производили с помощью авторской программы для ЭВМ EREGRE [22] для получения уравнения линейной множественной регрессии, приведенного ниже: (10) Из уравнения (10) видно, что выбранные нами факторы оказывают влияние на целевую функцию, а выбранные уровни факторов отвечают требованиям плана эксперимента. Заключение. Приведенная методика планирования экспериментов предназначена для экспериментальных исследований различных карданных передач с разнообразными типами подшипников [23], в том числе и авторской разработки [24; 25]. Получение полиномиальной регрессионной модели целевой функции - крутящего момента (затраченной работы на вращение карданной передачи) от выбранных конструктивно-режимных параметров карданной передачи - позволит в дальнейшем использовать методы поиска минимальных значений целевой функции, а следовательно, и коэффициента полезного действия передачи и соответствующих этим значениям конструктивно-режимных параметров, позволяющих в дальнейшем на этапе проектирования подобных передач обоснованно их назначать, что позволит сократить сроки проектирования при одновременном повышении качества принимаемых проектных решений, а в эксплуатации значительно экономить энергозатраты при одновременном повышении ресурса передачи.
×

作者简介

E. Kukushkin

Reshetnev Siberian State Aerospace University

Email: ironjeck@mail.ru
31, Krasnoyarsky Rabochy Av., Krasnoyarsk, 660037, Russian Federation

S. Eresko

Reshetnev Siberian State Aerospace University

31, Krasnoyarsky Rabochy Av., Krasnoyarsk, 660037, Russian Federation

T. Eresko

Reshetnev Siberian State Aerospace University

31, Krasnoyarsky Rabochy Av., Krasnoyarsk, 660037, Russian Federation

V. Menovshikov

Reshetnev Siberian State Aerospace University

31, Krasnoyarsky Rabochy Av., Krasnoyarsk, 660037, Russian Federation

I. Khomenko

Reshetnev Siberian State Aerospace University

31, Krasnoyarsky Rabochy Av., Krasnoyarsk, 660037, Russian Federation

参考

  1. Ереско С. П. Закономерность трения эластичных материалов по шероховатым поверхностям при наличии смазки // Проблемы машиностроения и надежность машин. 2002. № 6. C. 58-61.
  2. Локшина Н. Г. Развитие конструкций игольчатых подшипников и их применение. М. : ВНИПП, 1967. 283 с.
  3. Кукушкин Е. В., Меновщиков В. А. Основные направления развития, улучшения и совершенствования рабочих характеристик карданных передач на игольчатых подшипниках // Решетневские чтения : материалы XVI Междунар. науч. конф. / СибГАУ. Красноярск, 2012. С. 254-256.
  4. Кукушкин Е. В., Меновщиков В. А., Ереско Т. Т. Конструкция стенда для проведения испытаний карданных шарниров на игольчатых подшипниках // Решетневские чтения : материалы XIX Междунар. науч. конференции / СибГАУ. Красноярск, 2015. C. 337-339.
  5. Климов А. А., Стручков А. В. Исследование динамической нагруженности трансмиссии бульдозерного агрегата на базе трактора класса 40 кН на грунтах 1-2 категорий // Вестн. КрасГАУ. 2008. № 1. С. 201-206.
  6. Ереско Т. Т., Кукушкин Е. В., Меновщиков В. А. Современное состояние вопроса по исследованию пластического деформирования при статическом контактном нагружении игольчатых подшипников // Механики XXI веку. Братск : БрГУ, 2014. С. 37-40.
  7. Пат. 153924 Российская Федерация, МПК7 G 01 M 13/02 (2006.01). Стенд для испытаний карданных передач / Ереско С. П., Ереско Т. Т., Кукушкин Е. В., Меновщиков В. А. № 2014147821/28 ; заяв. 26.11.2014 ; опубл. 10.08.2015, Бюл. № 22. 2 с.
  8. Пат. 149002 Российская Федерация, МПК7 G 01 M 13/02 (2006.01). Стенд для испытаний карданных передач / Кукушкин Е. В., Меновщиков В. А., Ереско С. П., Ереско Т. Т. № 2014120845 ; заяв. 22.05.2014 ; опубл. 20.12.2014, Бюл. № 35. 1 с.
  9. Пат. 162876 Российская Федерация, МПК7 G 01 M 13/02 (2006.01). Стенд для испытаний карданных передач / Ереско С. П., Ереско А. С., Ереско Т. Т., Ереско В. С., Кукушкин Е. В., Стручков А. В., Хоменко И. И. № 2015157365 ; заяв. 30.12.2015 ; опубл. 27.06.2016, Бюл. № 18. 2 с.
  10. Кукушкин Е. В., Меновщиков В. А., Ереско Т. Т. Конструкция стенда для проведения испытаний карданных шарниров на игольчатых подшипниках // Решетневские чтения : материалы XIX Междунар. науч. конф. : в 2 ч. / под общ. ред. Ю. Ю. Логинова ; Сиб. гос. аэрокосмич. ун-т. Красноярск, 2015. С. 337-339.
  11. Меновщиков В. А., Ереско С. П. Исследование и совершенствование игольчатых подшипников карданных передач транспортно-технологических машин : монография. Красноярск : Изд-во КрасГАУ, 2006. 283 с.
  12. Hyun S., Robbins M. O. Elastic contact between rough surfaces: Effect of roughness at large and small wavelengths // Trobology International. 2007. Vol. 40. P. 1413-1422.
  13. Оценка достоверности результатов вычислительного эксперимента при проектировании механических и гидромеханических приводов / С. П. Ереско [и др.] // Решетневские чтения : материалы XVI Междунар. научной конф. / СибГАУ. Красноярск, 2012. С. 241-242.
  14. Шевцов С. М. Автоматизация процессов измерения вибрации // Механики XXI веку. 2008. № 7. C. 38-40.
  15. Popov V. L. Kontaktmechanik und Reibung. Ein Lehr- und Anwendungsbuch von der Nanotribologie bis zur numerischen Simulation. Springer-Verlag, 2009. 328 p.
  16. Решение математической модели динамики привода экспериментального стенда для исследования карданных передач / А. В. Стручков [и др.] // Наземные транспортно-технологические комплексы и средства : материалы Междунар. науч.-техн. конф. Тюмень : ТГНУ, 2016. С. 303-307.
  17. Решение математической модели динамики механической трансмиссии бульдозера с учетом диссипативных сил / А. В. Стручков [и др.] // Системы. Методы. Технологии. 2010. № 6. С. 58-65.
  18. Планирование эксперимента по исследованию карданных передач на игольчатых подшипниках / С. П. Ереско [и др.] // Актуальные проблемы авиации и космонавтики : материалы II Всерос. науч.-практ. конф. : в 2 т. / под общ. ред. Ю. Ю. Логинова ; Сиб. гос. аэрокосмич. ун-т. Красноярск, 2016. C. 368-370.
  19. Кожевников С. Н. Динамика нестационарных процессов в машинах. Киев : Наукова думка, 1996. 288 с.
  20. Кожевников С. Н. Динамика машин с упругими звеньями. Киев : Наукова думка, 1981. 160 с.
  21. Определение динамических параметров привода экспериментального стенда для исследования карданных передач / А. В. Стручков [и др.] // Вестник СибГАУ. 2016. Том 17, № 3. С. 638-644.
  22. Регрессионный анализ многофакторных экспериментальных исследований (EREGRE) : свид. об офиц. регистр. программы для ЭВМ / С. П. Ереско. № 2004610534 (РФ) ; заявл. 24.12.2003, № 2003612713 ; зарегистр. Роспатентом 24.02.2004.
  23. Сравнительный анализ конструкций карданных шарниров неравных угловых скоростей / С. П. Ереско // Вестник СибГАУ. 2015. Т. 16, № 3. С. 720-728.
  24. Пат. 141878 Российская Федерация, МПК7 F 16 D 3/26. Карданный шарнир / Кукушкин Е. В., Меновщиков В. А., Ереско С. П., Ереско Т. Т. № 2014102339/11 ; заявл. 24.01.2014 ; опубл. 20.06.2014, Бюл. № 17. 2 с.
  25. Пат. 146989 Российская Федерация, МПК7 F16 D 3/26. Карданный шарнир / Кукушкин Е. В., Меновщиков В. А., Орлов А. А., Ереско С. П., Ереско Т. Т. № 2014119234/11 ; заявл. 13.05.2001 ; опубл. 27.10.2014, Бюл. № 30. 2 с.

补充文件

附件文件
动作
1. JATS XML

版权所有 © Kukushkin E.V., Eresko S.P., Eresko T.T., Menovshikov V.A., Khomenko I.I., 2016

Creative Commons License
此作品已接受知识共享署名 4.0国际许可协议的许可
##common.cookie##