СПОСОБ НЕРАЗРУШАЮЩЕГО ОПРЕДЕЛЕНИЯ ТЕМПЕРАТУРНЫХ ХАРАКТЕРИСТИК СТРУКТУРНЫХ ПЕРЕХОДОВ В ПОЛИМЕРНЫХ МАТЕРИАЛАХ Российский патент 2013 года по МПК G01N25/02 G01N25/18 

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

Изобретение относится к области физико-химического анализа, а именно к способам определения температурных характеристик структурных переходов в полимерных материалах.

Известен способ контроля качества полимерных материалов за счет наблюдения за параметрами превращения материала, сопровождающегося изменениями теплоемкости (авт. св. СССР №1437758, МПК G01N 25/18, 1988 г., бюл. №42). Измерения проводят на непрерывно движущихся протяженных образцах (лентах, трубах). Способ включает локальный подвод тепловой мощности к образцу, охлаждение образца вне зоны нагрева, регистрацию температуры образца на фиксированном расстоянии от источника тепла. Точность контроля достигается за счет наблюдения за параметрами превращения полимерного материала, сопровождающегося изменениями теплоемкости, выбирают значение тепловой мощности, обеспечивающее нагрев образца выше температуры превращения, в процессе нагрева непрерывно регистрируют температуру образца вблизи зоны нагрева и разность значений температуры образца в двух точках, разно удаленных от источника тепла, а о качестве материала судят по разности значений температуры в зоне нагрева, соответствующих моментам скачкообразных изменений регистрируемой разности значений температуры. Однако данный способ не является оперативным, предназначен для протяженных движущихся образцов и не обеспечивает достаточную точность определения температурных характеристик структурных превращений в полимерах.

Наиболее близким к предлагаемому по технической сущности является способ определения теплофизических характеристик твердых материалов (патент РФ №2167412, МПК G01N 25/18, 2001 г., бюл. №14), заключающийся в том, что исследуемое полимерное тело приводят в тепловой контакт с эталонным телом по плоскости, в которой находится локальный в виде круга нагреватель постоянной мощности, измеряют разность значений температуры между нагревателем и точкой в плоскости контакта исследуемого и эталонного тел, расположенной на расстоянии более чем четыре радиуса нагревателя от центра нагревателя и не более наименьшей толщины исследуемого тела, до тех пор, когда эта разность значений температуры не станет меньше наперед заданной величины, непрерывно подводят постоянную мощность на нагреватель, через равные промежутки времени измеряют разность значений температуры, на каждом i-ом шаге измерения определяют величину динамического параметра, сравнивают величину динамического параметра с наперед заданным значением, испытания заканчивают при превышении контролируемым динамическим параметром заданного значения, определяют тепловую активность исследуемого материала.

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

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

Это достигается тем, что исследуемое полимерное тело приводят в тепловой контакт с эталонным телом по плоскости, в которой находится локальный в виде круга нагреватель постоянной мощности, измеряют разность значений температуры между нагревателем и точкой плоскости контакта исследуемого и эталонного тел, расположенной на расстоянии более чем четыре радиуса нагревателя от центра нагревателя и не более наименьшей толщины исследуемого тела, до тех пор, когда эта разность значений температуры не станет меньше наперед заданной величины, фиксируют температуру исследуемого тела, затем непрерывно подводят постоянную мощность на нагреватель, через равные промежутки времени измеряют разность значений температуры, дополнительно на каждом i-ом шаге измерения определяют величину динамического параметра χ i = d 1 max d 1 i d 1 max , где d1i определяют согласно формуле

d 1 i = j = i ( k 1 ) / 2 i + ( k 1 ) / 2 T j ( t j t i ) j = i ( k 1 ) / 2 i + ( k 1 ) / 2 ( t j t i ) 2 , ( 1 )

где t i = 1 k j = i ( k 1 ) / 2 i + ( k 1 ) / 2 t j , t = τ ,

Tj - разность значений температуры на j шаге измерения, k - целое положительное нечетное число большее 3, d1max - максимальное значение по абсолютной величине из всех расчетных значений d1i к текущему моменту времени, сравнивают величину динамического параметра с наперед заданным значением, нагреватель отключают при превышении динамическим параметром заданного значения, строят зависимость текущего значения коэффициента тепловой активности ε i * от температуры исследуемого тела, а о структурных переходах в полимерном материале судят по наличию пиков зависимости значений ε i * , от температуры исследуемого тела, где ε i * , непосредственно определяют из снятой зависимости разности значений температуры между нагревателем и точкой плоскости контакта образца и эталонного тела от времени по формулам

ε i * = A d 1 i ε ' ,                                         (2)

где A = d 11 d 12 d 12 d 11 ( ε 01 ε 02 ) , ε ' = ε 1 d 11 ε 2 d 12 d 12 d 11 ,

d1i - текущие значения параметра d1; d11, d12, ε01, ε02 - параметры d1 и тепловые активности образцовых мер, соответственно.

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

Сущность предлагаемого способа поясняется следующим теоретическим обоснованием.

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

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

T ( 0,0, τ ) = 2 q α τ λ ( 1 π i e r f c [ R 2 a τ ] ) ,                      (3)

где q - плотность теплового потока, Вт/м2; R - радиус нагревателя, м; а - температуропроводность, м2/с; τ - время, с.

Причем при малых τ: T ( 0,0, τ ) 2 q τ ε π .

Средняя температура по нагревателю при малых т будет определяться выражением

S ( τ ) 2 q τ ε π                                                (4)

Для плоского полупространства температура нагревателя будет определяться выражением

T п л ( 0, τ ) 2 q п л τ ε π                                             (5)

где qпл - плотность теплового потока на поверхности плоского полупространства, Вт/м2.

При малых значениях времени т температурное поле от ограниченного круглого нагревателя будет аналогично температурному полю от нагревателя, действующего в плоском полупространстве (патент РФ №2167412, МПК G01N 25/18, 2001 г., бюл. №14). Это справедливо как для центральной точки нагревателя, так и для средней по нагревателю, а значит и для любой другой точки нагревателя.

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

В общем случае на термограмме можно выделить три участка (фиг.1).

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

2 T 2 T z 2 , ( 2 T r 2 + 1 r T r 0 ) , qo=var, qo+q3≠q,

где z, r - координаты, м; qo - тепловой поток, поступающий в образец или изделие; Вт/м2; q3 - тепловой поток, поступающий в подложку зонда, т.е. эталонное тело, Вт/м2; q - мощность, выделяющаяся на единицу площади нагревателя, Вт/м2; ∇2 - оператор Лапласа.

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

2 T 2 T z 2 , ( 2 T r 2 + 1 r T r 0 ) , qo=const, qo+q3≈q, ( q o q з ε з ε )

где ε, εЗ - значения тепловой активности исследуемого материала и материала подложки зонда, т.е. эталонного тела, B т с 0,5 м 2 К .

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

2 T 2 T z 2 + 2 T r 2 + 1 r T r , qo=var, qз=var.

Расчетное выражение, описывающее термограмму на втором участке, имеет вид (патент РФ №2167412, МПК G01N 25/18, 2001 г., бюл. №14):

T * ( t ) = 2 q t ( ε + ε ' ) π q c н ( ε + ε ' ) 2 ,                            (6)

где T* - разность значений температуры, регистрируемая дифференциальной термопарой, °С; t = τ - новая координата, с0,5; q - удельная, на единицу площади, мощность на нагревателе, Вт/м2; ε, ε' - соответственно, тепловая активность исследуемого материала и поправка на тепловую активность материала зонда, т.е. эталонного тела для II участка, B т с 0,5 м 2 К ; сн - теплоемкость нагревателя отнесенная к единице площади, Дж/(м2·К);

Выражение (6) можно записать в виде:

T * ( t ) = d 1 t + d 0                                                  (7)

где

A = 2 q π , B = q c н , d 1 = A ( ε + ε ' ) , d 0 = B ( ε + ε ' ) 2 ,                               (8)

Значения d1, d0 определяются из термограммы T*(t), а значения А, В, ε' - являются постоянными прибора и не зависят от теплофизических свойств исследуемого материала.

Модель (7) работает на II участке термограммы, а также частично на I и III.

Модель (7) работает только вне зоны структурных переходов. При этом погрешность определения коэффициента d1i на II участке будет постоянной, а на III участке монотонно изменяться. Значение коэффициента d1i обратно пропорционально значению тепловой активности ε i * (2). В свою очередь тепловая активность связана с теплоемкостью материала формулой

ε = c ρ λ ,

где с - удельная теплоемкость материала, ρ - плотность материала, λ - теплопроводность. При структурном переходе в материале в узком температурно-временном интервале значение теплоемкости изменяется значительно (часто на порядок и более), следовательно, значение ε также изменится.

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

НИ>L>4R,

где R - радиус нагревателя, Ни - толщина образца (изделия).

При подготовке к испытаниям создают тепловой контакт между нагревателем и образцом, а также между датчиком и образцом. Осуществляют контроль за процессом термостатирования с помощью измерительно-вычислительного устройства (ИВУ) 5. Когда величина разности значений температуры Т* станет меньше наперед заданного значения, определяемого точностью измерения температуры, измеряют температуру исследуемого тела, ИВУ подает на нагреватель с помощью стабилизированного источника питания 6 электрический ток постоянной мощности. Одновременно с подачей электрического тока производят измерение разности значений температуры Т* через равные промежутки времени Δτ.

На каждом i-м временном шаге изменения температуры контролируют значение динамического параметра χi. Сравнивают величину динамического параметра χi на i-м шаге измерения с его значением на i-1 шаге. Испытания заканчивают при превышении контролируемым динамическим параметром заданного значения χ3=0,03÷0,07. Строят зависимость параметра ε i * от температуры исследуемого тела, а о структурных переходах в материале судят по наличию пиков на зависимости параметра ε i * от температуры исследуемого тела, где ε i * непосредственно определяют из снятой зависимости разности значений температуры между нагревателем и точкой плоскости контакта образца и эталонного тела от времени по формуле (2).

Примеры №1 и №2 конкретной реализации способа осуществляли при определении температуры твердофазных переходов в политетрафто-рэтилене (ПТФЭ). Известно, что политетрафторэтилен (ГОСТ 10007-80) имеет следующие кристаллические модификации. При температуре ниже 19,6°С его элементарная ячейка имеет триклиническую структуру. В интервале от 19,6°С до 30°С существует гексагональная элементарная ячейка. Упорядоченность структуры при этом уменьшается, решетка расширяется, удельный объем увеличивается, примерно, на 1%. При температуре выше 30°С стабильной становится псевдогексагональная решетка. Изменение объема образца при этом составляет 0,08%. Удельные теплоты переходов составляют соответственно: 4,0±0,5 кДж/кг и 1,2±0,3 кДж/кг.

Пример №1. Испытания проводили при следующих конструктивных условиях: радиус нагревателя R=2,5·10-3 м; толщина образца из политетрафторэтилена Ни=15·10-3 м; расстояние между центром нагревателя, где установлен горячий спай термопары, и холодным спаем термопары L=12·10-3 м; толщина подложки Нэ=20·10-3 м. В качестве эталонного тела использовали подложку зонда из Рипора. Условия испытания: мощность на нагревателе W=0,6 Вт; шаг измерения температуры Δτ=0,5 сек. Испытания закончили в момент времени 10 сек. В качестве образцовых мер использовали: полиметилметакрилат (ГОСТ 10667-90) и стекло марки KB (ГОСТ 15130-86). Конструктивные и режимные условия испытания для примера №1 представлены в таблице 1. Затем строили зависимость ε i * от температуры образца (фиг.3). На данной зависимости обнаружен характерный пик при 29,8°С, соответствующий твердофазному структурному переходу в ПТФЭ при 30°С, что согласуются с известными литературными данными о переходах для данного материала.

Пример №2. Аналогичные испытания проводили при следующих конструктивных и режимных условиях, представленных в таблице 1. Испытания закончили в момент времени 40 сек (фиг.4). Затем строили зависимость ε i * (3) от температуры образца (фиг.5). На данной зависимости обнаружен характерный пик при 29,8°С, соответствующий твердофазному структурному переходу в ПТФЭ при 30°С, что согласуются с известными литературными данными о переходах для данного материала.

Пример №3. Испытания изделия из коксонаполненного политетраф-торэтилена Ф4К20 проводили при следующих конструктивных и режимных условиях, представленных в таблице 1. Испытания закончили в момент времени 20 сек (фиг.6). Затем строили зависимость ε i * (3) от температуры образца (фиг.7). На данной зависимости обнаружен характерный пик при 30°С, соответствующий твердофазному структурному переходу в ПТФЭ при 30°С, что согласуются с известными литературными данными о переходах для данного материала.

Пример №4. Испытания изделий из блочного полиамида Капролона В проводили при следующих конструктивных и режимных условиях, представленных в таблице 1. Испытания закончили в момент времени 38 сек (фиг.8). Затем строили зависимость ε i * (3) от температуры образца (фиг.9). Твердофазное полиморфное превращение в Капролоне В зафиксировано в интервале температур 24,7…25,3°С.

Таблица 1 Конструктивные и режимные условия испытаний Пример реализации № п/п Материал Конструктивные условия Режимные условия R·103, м Ни·103, м L·103 Нэ·103, м Материал подложки зонда Тн, °С W, Вт Длительность испытаний, с Δτ, с Образцо
вые меры теплофизических свойств
1 ПТФЭ 2,5 15 12 20 Рипор 20 0,6 10 0,5 ПММА Стекло марки KB 2 ПТФЭ 4 25 18 23 0,7 40 3 Ф4К20 4 25 18 26 0,7 20 4 Капролон В 4 25 18 21 0,7 38

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

название год авторы номер документа
СПОСОБ НЕРАЗРУШАЮЩЕГО ОПРЕДЕЛЕНИЯ ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ ТВЕРДЫХ МАТЕРИАЛОВ 2006
  • Жуков Николай Павлович
  • Майникова Нина Филипповна
  • Рогов Иван Владимирович
  • Чех Алексей Сергеевич
  • Никулин Сергей Сергеевич
RU2328725C1
СПОСОБ НЕРАЗРУШАЮЩЕГО ОПРЕДЕЛЕНИЯ ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ ТВЕРДЫХ МАТЕРИАЛОВ 2005
  • Жуков Николай Павлович
  • Майникова Нина Филипповна
  • Муромцев Юрий Леонидович
  • Чех Алексей Сергеевич
  • Никулин Сергей Сергеевич
RU2287152C1
СПОСОБ КОМПЛЕКСНОГО ОПРЕДЕЛЕНИЯ ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ МАТЕРИАЛОВ 1999
  • Жуков Н.П.
  • Майникова Н.Ф.
  • Муромцев Ю.Л.
  • Рогов И.В.
  • Орлов В.В.
RU2167412C2
СПОСОБ НЕРАЗРУШАЮЩЕГО ОПРЕДЕЛЕНИЯ ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ МАТЕРИАЛОВ 1999
  • Жуков Н.П.
  • Муромцев Ю.Л.
  • Майникова Н.Ф.
  • Рогов И.В.
  • Балашов А.А.
RU2161301C2
СПОСОБ ОЦЕНКИ РАЗЛИЧИЯ ТЕПЛОФИЗИЧЕСКИХ ПАРАМЕТРОВ ВИДИМОЙ ПОВЕРХНОСТИ ИЗОТРОПНОГО ОБЪЕКТА С УЧЕТОМ ФОНА 2013
  • Антонов Борис Игоревич
  • Обухов Владимир Васильевич
  • Парфирьев Андрей Владимирович
  • Ищук Игорь Николаевич
  • Попело Владимир Дмитриевич
RU2544894C1
СПОСОБ ИДЕНТИФИКАЦИИ КОМПЛЕКСА ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ ТВЕРДЫХ МАТЕРИАЛОВ 2005
  • Обухов Владимир Васильевич
  • Ищук Игорь Николаевич
  • Фесенко Александр Иванович
  • Собко Александр Павлович
  • Антонов Борис Игоревич
RU2303777C2
СПОСОБ И УСТРОЙСТВО ДЛЯ ИЗМЕРЕНИЯ СТЕПЕНИ ЧЕРНОТЫ 2012
  • Лаповок Евгений Владимирович
  • Пеньков Максим Михайлович
  • Слинченко Дмитрий Анатольевич
  • Уртминцев Игорь Александрович
  • Ханков Сергей Иванович
RU2521131C2
СПОСОБ НЕРАЗРУШАЮЩЕГО КОНТРОЛЯ ТЕПЛОФИЗИЧЕСКИХ ХАРАКТЕРИСТИК СТРОИТЕЛЬНЫХ МАТЕРИАЛОВ И ИЗДЕЛИЙ 2014
  • Жарикова Мария Валерьевна
  • Чернышов Алексей Владимирович
  • Чернышов Владимир Николаевич
RU2574229C1
СПОСОБ ИЗМЕРЕНИЯ ИНТЕГРАЛЬНОЙ ИЗЛУЧАТЕЛЬНОЙ СПОСОБНОСТИ С ПОМОЩЬЮ ПРЯМОГО ЛАЗЕРНОГО НАГРЕВА (ВАРИАНТЫ) 2015
  • Брыкин Михаил Владимирович
  • Васин Андрей Андреевич
  • Шейндлин Михаил Александрович
RU2597937C1
Способ измерения интегральной излучательной способности с применением микропечи (варианты) 2015
  • Брыкин Михаил Владимирович
  • Васин Андрей Андреевич
  • Шейндлин Михаил Александрович
RU2607671C1

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

Реферат патента 2013 года СПОСОБ НЕРАЗРУШАЮЩЕГО ОПРЕДЕЛЕНИЯ ТЕМПЕРАТУРНЫХ ХАРАКТЕРИСТИК СТРУКТУРНЫХ ПЕРЕХОДОВ В ПОЛИМЕРНЫХ МАТЕРИАЛАХ

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

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

Способ неразрушающего определения температурных характеристик структурных переходов в полимерных материалах, заключающийся в том, что исследуемое полимерное тело приводят в тепловой контакт с эталонным телом по плоскости, в которой находится локальный в виде круга нагреватель постоянной мощности, измеряют разность значений температуры между нагревателем и точкой плоскости контакта исследуемого и эталонного тел, расположенной на расстоянии более чем четыре радиуса нагревателя от центра нагревателя и не более наименьшей толщины исследуемого тела, до тех пор, когда эта разность значений температуры не станет меньше наперед заданной величины, фиксируют температуру исследуемого тела, затем непрерывно подводят постоянную мощность на нагреватель, через равные промежутки времени измеряют разность значений температуры, отличающийся тем, что на каждом i-м шаге измерения определяют величину динамического параметра χ i = d 1 max d 1 i d 1 max ,
где dli определяют согласно формулам
d 1 i = j = i ( k 1 ) / 2 i + ( k 1 ) / 2 T j ( t j t i ) j = i ( k 1 ) / 2 i + ( k 1 ) / 2 T j ( t j t i ) 2 ,
t i = 1 k j = i ( k 1 ) / 2 i + ( k 1 ) / 2 t j , t = τ i ,
где Tj - разность значений температуры на j шаге измерения, τj - время, соответствующее j шагу измерения, k - целое положительное число большее 3, d1max - максимальное значение по абсолютной величине из всех расчетных значений dli к текущему моменту времени, сравнивают величину динамического параметра с наперед заданным значением, испытания заканчивают при превышении контролируемым динамическим параметром заданного значения χ3=0,03÷0,07, строят зависимость текущего значения коэффициента тепловой активности ε i * от температуры исследуемого тела, а о структурных переходах в материале судят по наличию пиков зависимости значений ε i * от температуры исследуемого тела, где ε i * непосредственно определяют из снятой зависимости разности значений температуры между нагревателем и точкой плоскости контакта образца и эталонного тела от времени по формулам
ε i * = A d 1 i ε ' ,
где A = d 11 d 12 d 12 d 11 ( ε 01 ε 02 ) , ε ' = ε 1 d 11 ε 2 d 12 d 12 d 11 ,
dli - текущие значения параметра d1; d11, d12, ε01, ε02 - параметры d1 и тепловые активности образцовых мер соответственно.

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

СПОСОБ КОМПЛЕКСНОГО ОПРЕДЕЛЕНИЯ ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ МАТЕРИАЛОВ 1999
  • Жуков Н.П.
  • Майникова Н.Ф.
  • Муромцев Ю.Л.
  • Рогов И.В.
  • Орлов В.В.
RU2167412C2
СПОСОБ НЕРАЗРУШАЮЩЕГО КОНТРОЛЯ ТЕПЛОФИЗИЧЕСКИХ ХАРАКТЕРИСТИК И ТЕМПЕРАТУРНЫХ ПЕРЕХОДОВ ТЕРМОПЛАСТОВ 2004
  • Чернышов В.Н.
  • Рожнова Л.И.
  • Радько Ю.М.
RU2250453C1
Устройство для определения температурных переходов в материалах 1980
  • Шутилин Юрий Федорович
SU873084A1
Способ определения температур структурных переходов в материалах 1979
  • Шутилин Юрий Федорович
  • Балашов Юрий Степанович
  • Моисеев Владимир Васильевич
  • Шеин Владимир Сергеевич
  • Розиноер Яков Моисеевич
SU873082A1
DE 3615344 A1, 02.01.1987
US 5217304 A, 08.06.1993.

RU 2 493 558 C1

Авторы

Майникова Нина Филипповна

Жуков Николай Павлович

Рогов Иван Владимирович

Балашов Алексей Александрович

Попов Олег Николаевич

Даты

2013-09-20Публикация

2012-04-10Подача