СПОСОБ ПРОВЕРКИ КАЧЕСТВА ФУНКЦИОНИРОВАНИЯ НИЗКОЧАСТОТНОГО ДИНАМИЧЕСКОГО БАЛАНСИРОВОЧНОГО СТЕНДА Российский патент 2011 года по МПК G01M1/04 

Описание патента на изобретение RU2434212C1

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

Современная технология создания высокоточных нежестких изделий машиностроения роторного типа, к которым, например, относятся тонкостенные спутниковые узлы, нежесткие сборки типа разворачивающихся в невесомости антенн, роторы, внутри которых нежестко размещена геофизическая аппаратура или сейсмические датчики, исключает возможность их балансировки при больших скоростях вращения и ограничивает рабочую частоту вращения до 1-2 Гц, а также требует вертикального расположения оси вращения с целью исключения погрешностей, связанных с прогибом изделия под действием сил тяжести. Балансировка должна осуществляться на специализированных вертикальных дорезонансных низкочастотных балансировочных станках (стендах), для которых характерна высокая точность балансировки. В частности, это относится к стендам с аэростатическими опорами, предназначенным для выполнения бесконтактной прецизионной балансировки изделий в динамическом режиме (после их изготовления и сборки) с определением параметров массо-инерционной асимметрии контролируемых изделий, (см., например, патент РФ 2292533, G01M 01/04, 2007 г.).

Понятие «параметры массо-инерционной асимметрии» не является новым, оно введено, например, для оценки факторов, определяющих движение в атмосфере Земли летательных аппаратов. В частности к таким параметрам относятся радиус - вектор смещения центра масс и вектор - угол отклонения продольной главной центральной оси инерции (ГЦОИ) от геометрической оси Х аппарата. Каждый из указанных параметров характеризуется значением и углом в связанной с контролируемым изделием системе координат.

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

Известен метод экспериментального определения точностных характеристик балансировочных станков при проверке чувствительности измерительной системы методом кругового обхода при установке контрольного груза известной массы в плоскостях коррекции контрольного ротора и с последующей оценкой точности станка по результатам измерений векторных параметров вибраций опор с использованием графических построений [Левит М.Е., Рыженков В.М. Балансировка деталей и узлов. - М.: Машиностроение, 1986, стр.97].

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

Известен метод оценки норм точности дорезонансного балансировочного станка с цепью разделения плоскостей коррекции, реализованной по схеме a, b, c, выполняемой по результатам балансировки контрольного ротора отдельно в каждой плоскости коррекции с использованием контрольных грузов, устанавливаемых в контролируемой плоскости данного ротора [там же, стр.156]. Метод этот также сложный, малоинформативный, трудоемкий, причем в процессе измерений не контролируются значения параметров массо-инерционной асимметрии объекта в заданных диапазонах измерений, а правильный результат может быть получен лишь как результат косвенных измерений после математической обработки полученных результатов испытаний ротора на станке.

Наиболее близким к предлагаемому способу проверки качества функционирования специализированного вертикального балансировочного стенда является способ проверки стабильности измерения амплитуд и фаз сигналов дисбалансов, выполняемой с применением жесткого, динамически сбалансированного эталонного ротора и набора калиброванных контрольных грузов, прикрепляемых к плоскостям коррекции эталонного ротора, расположенным вблизи каждого из его торцов и имеющим известные радиусы. Эталонный ротор характеризуется тем, что его массо-центровочные и инерционные характеристики, в том числе, значения параметров массо-инерционной асимметрии определены с точностью, не менее чем в 2,5-3 раза превышающей точность определения тех же характеристик балансируемого ротора, а его базовые посадочные поверхности соответствуют аналогичным поверхностям балансируемого ротора. Оценка точности измерительной системы стенда производится по итогам суммарной обработки результатов многократных измерений (требуемое суммарное количество требуемых циклов измерений - более 100 для каждой опоры) контролируемых параметров сигналов дисбалансов, поступающих одновременно от силоизмерительных датчиков, установленных в упругих элементах обеих опор, путем их сравнения с нормами для каждой из опор (Глазырина Л.М., Ключников А.В., Кругов B.C. «Методика определения нормируемых метрологических характеристик балансировочного стенда». Датчики систем измерения, контроля и управления: Межвузовский сборник научных трудов. - Вып.23. - Пенза: Изд-во Пенз. гос. ун-та, 2003, с.81-86). При этом масса контрольных грузов подбирается таким образом, чтобы максимальный дисбаланс в контролируемой плоскости коррекции не превышал максимально допустимого для данной плоскости значения.

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

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

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

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

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

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

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

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

На фиг.2 приведена схема для проведения расчетов двухмерных (векторных) эталонных параметров массо-инерционной асимметрии контролируемого ротора, где 1 - геометрическая ось Х ротора, 2 - продольная ГЦОИ, проходящая через центр масс 3, расположенный на расстоянии хВ от верхней плоскости коррекции 4 и на расстоянии хН от нижней плоскости коррекции 5, 6 и 7 - соответственно верхняя и нижняя аэростатические опоры. Под эталонными векторными значениями параметров массо-инерционной асимметрии понимаются значения указанных параметров - смещение центра масс и угол перекоса продольной ГЦОИ 2 относительно геометрической оси 1, - задаваемые (моделируемые) в статических условиях путем прикрепления к верхней 4 и нижней 5 плоскостям коррекции эталонного ротора контрольных грузов известной массы (соответственно mkB и mkH) на известных радиусах (соответственно rB и rH) и с известными угловыми положениями (соответственно φkB и φkH) в системе координат, связанной с эталонным ротором.

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

Затем в каждой проверяемой точке диапазона измерений контролируемого параметра массо-инерционной асимметрии, задаваемого путем прикрепления в известном угловом положении к верхней 4 и (или) нижней 5 плоскости коррекции эталонного ротора контрольных грузов известной массы mkB и mkH, выбираемых из набора калиброванных грузов, выполняют небольшое количество измерений (достаточно провести не более пяти измерений, что даст вполне достоверную информацию о качестве функционирования стенда) и определяют амплитуду и фазу вибраций опор 6 и 7. Для каждого измерения с использованием полученных коэффициентов балансировочной чувствительности и взаимовлияния плоскостей коррекции рассчитывают значения и угловые положения дисбалансов, действующих в верхней 4 и нижней 5 плоскости коррекции. По результатам расчетов определяют значения контролируемого параметра массо-инерционой асимметрии. Модуль абсолютной погрешности каждого i-го измерения значения контролируемого параметра |ΔП| - определяется как геометрическая разность между двухмерным измеренным значением контролируемого параметра , а также соответствующим двухмерным заданным (эталонным) значением параметра и двухмерным значением того же параметра , характеризующим собственную массо-инерционную асимметрию используемого эталонного ротора по формуле

откуда модуль |Δect| абсолютной погрешности определения поперечного смещения центра масс и модуль |Δαх| абсолютной погрешности определения угла перекоса продольной ГЦОИ каждого i-го измерения определяются соответственно по формулам

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

где mkB и mkH - контрольные грузы, установленные соответственно на углы φkB и φkH, в верхней 4 и(или) нижней 5 плоскости коррекции;

Мрот - масса эталонного ротора;

rB - расстояние от центра верхней плоскости коррекции 4 до центра масс контрольного груза mkB с учетом его углового положения в системе координат эталонного ротора;

rH - расстояние от центра нижней плоскости коррекции 5 до центра масс контрольного груза mkH с учетом его углового положения в системе координат эталонного ротора;

xB и xH - расстояния от центра масс эталонного ротора 3 до соответственно верхней 4 и нижней 5 плоскости коррекции после установки контрольных грузов mkB и(или) mkH;

xB нач и xH нач - начальные (т.е. измеренные до установки контрольных грузов) расстояния от центра масс до соответствующих плоскостей коррекции;

Ia рот, и Iэ рот - соответственно аксиальный и экваториальный моменты инерции эталонного ротора.

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

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

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

Похожие патенты RU2434212C1

название год авторы номер документа
СПОСОБ НАСТРОЙКИ БАЛАНСИРОВОЧНОГО СТЕНДА ДЛЯ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ МАССОИНЕРЦИОННОЙ АСИММЕТРИИ РОТОРОВ 2011
  • Ключников Александр Васильевич
RU2453818C1
СПОСОБ БАЛАНСИРОВКИ РОТОРА В ОДНОЙ ПЛОСКОСТИ КОРРЕКЦИИ 2018
  • Ключников Александр Васильевич
RU2694142C1
СПОСОБ БАЛАНСИРОВКИ РОТОРА В ОДНОЙ ПЛОСКОСТИ КОРРЕКЦИИ 2012
  • Ключников Александр Васильевич
RU2499985C1
СПОСОБ БАЛАНСИРОВКИ РОТОРА 2004
  • Глазырина Любовь Митрофановна
  • Карповицкий Михаил Степанович
  • Ключников Александр Васильевич
  • Мальгин Анатолий Иванович
  • Смирнов Геннадий Григорьевич
  • Фомин Юрий Павлович
RU2292534C2
Способ динамической балансировки ротора в двух плоскостях коррекции 2022
  • Скачков Виктор Сергеевич
RU2790897C1
Способ определения погрешности стенда для измерения характеристик геометрии масс изделий и устройство для его осуществления 2019
  • Довыденко Ольга Владимировна
  • Самойленко Александр Иванович
  • Бугров Александр Юрьевич
  • Лютов Владимир Васильевич
  • Куликов Александр Андреевич
RU2722962C1
Способ сборки и балансировки высокооборотных роторов и валопроводов авиационных газотурбинных двигателей и газоперекачивающих агрегатов 2022
  • Сусликов Виктор Иванович
  • Сусликов Сергей Викторович
  • Болотов Михаил Александрович
RU2822671C2
СПОСОБ ПРОВЕРКИ КАЧЕСТВА ФУНКЦИОНИРОВАНИЯ СТЕНДА ДЛЯ ОПРЕДЕЛЕНИЯ МАССО-ЦЕНТРОВОЧНЫХ И МАССО-ИНЕРЦИОННЫХ ХАРАКТЕРИСТИК ТВЕРДОГО ТЕЛА ВРАЩЕНИЯ 2010
  • Ключников Александр Васильевич
RU2445592C1
СПОСОБ ВИБРАЦИОННОЙ БАЛАНСИРОВКИ РОТОРОВ И КОЛЕБАТЕЛЬНАЯ СИСТЕМА БАЛАНСИРОВОЧНОГО СТАНКА ДЛЯ ЕГО ОСУЩЕСТВЛЕНИЯ 1992
  • Манько А.Т.
RU2101689C1
Способ балансировки ротора газоперекачивающего агрегата 2017
  • Осач Евгений Александрович
  • Степанов Андрей Владимирович
RU2658597C1

Иллюстрации к изобретению RU 2 434 212 C1

Реферат патента 2011 года СПОСОБ ПРОВЕРКИ КАЧЕСТВА ФУНКЦИОНИРОВАНИЯ НИЗКОЧАСТОТНОГО ДИНАМИЧЕСКОГО БАЛАНСИРОВОЧНОГО СТЕНДА

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

Формула изобретения RU 2 434 212 C1

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

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

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

Документы, цитированные в отчете о поиске Патент 2011 года RU2434212C1

ГЛАЗЫРИНА Л.М., КЛЮЧНИКОВ А.В., КРУГОВ B.C
«Методика определения нормируемых метрологических характеристик балансировочного стенда»
Датчики систем измерения, контроля и управления: Межвузовский сборник научных трудов
Прибор для равномерного смешения зерна и одновременного отбирания нескольких одинаковых по объему проб 1921
  • Игнатенко Ф.Я.
  • Смирнов Е.П.
SU23A1
- Пенза: Изд-во Пенз
гос
ун-та, 2003, с.81-86
Левит М.Е., РЫЖЕНКОВ В.М
Балансировка деталей и узлов
- М.:

RU 2 434 212 C1

Авторы

Ключников Александр Васильевич

Фомин Юрий Павлович

Даты

2011-11-20Публикация

2010-02-24Подача