Научная статья на тему 'Влияние колебаний динамической системы станка на точность и температуру при шлифовании'

Влияние колебаний динамической системы станка на точность и температуру при шлифовании Текст научной статьи по специальности «Механика и машиностроение»

CC BY
227
52
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по механике и машиностроению, автор научной работы — Никитин Сергей Петрович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Влияние колебаний динамической системы станка на точность и температуру при шлифовании»

УДК 621.9: 531.3:001.5

С.П. Никитин

Пермский государственный технический университет

ВЛИЯНИЕ КОЛЕБАНИЙ ДИНАМИЧЕСКОЙ СИСТЕМЫ СТАНКА НА ТОЧНОСТЬ И ТЕМПЕРАТУРУ ПРИ ШЛИФОВАНИИ

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

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

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

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

Рис. 1. Динамическая система «станок - заготовка - инструмент» с процессом шлифования и правки

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

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

Рис. 2. Расчетная схема плоскошлифовального станка

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

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

Шлифование как процесс пластического деформирования и разрушения материалов детали и круга занимает промежуточное положение между резанием резцом [3, 4] (инструментом с определенной режущей кромкой) и трением. Поэтому при моделировании динамики процесса при шлифовании необходимо учесть особенности этого процесса. Математические модели

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

Рис. 3. Схема контакта шлифовального круга с заготовкой

Процесс удаления материала с заготовки на длине контакта 1д обеспечивают отдельные зерна. Пока еще нет достаточной информации о процессах стружкообразования, силах и трении при резании абразивными и алмазными зернами. Но, ссылаясь на публикации [1, 9, 10] и информацию о виде стружек при шлифовании, можно предположить близкое сходство процесса резания зерном и резцом (рис. 4). Процесс отделения стружки отдельным зерном содержит три фазы упругопластического взаимодействия: фаза чистой упругой деформации, фаза упругой и пластической деформации, фаза упругой и пластической деформации с удалением материала.

Рзг

Рис. 4. Процесс резания отдельным зерном

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

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

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

Для определения учтем систему сил, действующую на режущий выступ абразивного зерна [4], и износ зерна при обработке.

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

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

где &рх, кРу - коэффициенты резания по осям х, у;

ЬРх, ЬРу, Ь?2 - коэффициенты демпфирования по координатам х, у, 2.

Коэффициенты резания по осям у и 2 при шлифовании можно вычислить по выражениям

где Р - угол между равнодействующей силы резания и скоростью резания;

(1)

р+Тр р - ь?хх+ Куу + Ку у+ к?22 + К z,

(2)

(3)

Р1 - угол сдвига;

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

В - ширина контакта круга с заготовкой;

1д - длина контакта круга с заготовкой; ц - коэффициент трения стружки о поверхность зерна.

Величину фактического числа зерен, участвующих в резании на единичной площадке контакта круга с заготовкой, можно определить исходя из законов теории вероятности по выражению [6, 7, 8]

0,4

Пр = Х”

Хср

Ю'К' , (4)

и, (1 -е) ' '

где хср - средняя величина зерна, мм;

К - концентрация режущего материала в круге, %;

ик - скорость шлифования, м/с;

е - относительная глубина заделки зерен, е к 0,7-0,8;

ю* - скорость погружения зерен в обрабатываемый материал, мм/с.

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

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

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

Для практического теплофизического анализа процесса шлифования будем использовать систему тел и источников в локальной области, прилегающей к зерну [6, 9, 10] (рис. 5).

Рис. 5. Режущий выступ зерна и расположение источников тепловыделения

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

Теплота от каждого элементарного источника распределяется между всеми телами, участвующими в процессе. Представим это распределение в виде итоговых потоков теплообмена между инструментом, заготовкой и стружкой [6, 7, 8, 9, 10].

Итоговые потоки тепла определяем как результат совокупности элементарных потоков, соединенных параллельно. Элементарные источники представим аналогично в виде соответствующих совокупных источников:

О, = qn р\д B,

где пр - средневероятное количество режущих зерен, активно участвующих в процессе шлифования на единице контакта;

1д - длина контакта;

В - ширина контакта.

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

J1 ф1 + к1 ф1 + С1 ф1 + С2 (к?ф1 - к1ф2) = 0

^2 ф2 ^ (ф2 — к1 " ф1 ) " ^2 ^ СВ (ф2 — ф3 ) = 0

^ фз + Сзф4 + Схг2фз + Скхггл — Сьгг5 — Св (ф2 —фз) = 0

Мгг4 + + Скг (г4 — г5) + = 0

—1СЬ (24 — ) — Ск,гфз — кр = 0

МхХ6 + Кх6 + СхХ6 + кхР = 0

Му У 7 + куУ7 + СуУ7 + кур = 0

^Рг 25 + ЙРх Х6 + кРхХ6 + ЙРу у7 + кРуУ7 + Р+ Р + ^ 214 + ^рх Х15 + кРхХ15 + ^Ру У16 + кРуУ16 = 0

ТР

Сс “ТТ+1+1Г)&9 - 1- ©10 — (КС — КС) Р—К^5 = 0

& 1с 1тп ^п 7

1 1

—~Г 09 + {1Т+1т)&10 — 1т ©11 — (*-1+ КП2) Р—= 0

1тп 1тп 11 1

——010 + СК &01 + (— + — + —)0П — — 012 = 0

1 10 к & 1К 1 1/ 11 1Г 12

— ~1~ 011 + (“1- + 1 )012 — 1 013 — КТЗ25 = 0ТЗ

1 1У 1ТЗ 1ТЗ

—1^ 012 + Си + (-1- + -1-)01з - (*И — ) Р — О, = 0

1ТЗ & 1ГЗ 1и

Мс 214 + ^14 214 + Сшт (214 — 217 ) + кхР = 0 Мс Х15 + Йх Х15 + СхХ15 + кхР = 0

Мс У16 + ЙУ У16 + СуУ16 + куР = 0

Мп 217 + ^17 217 + С17 217 — Сшт (214 — 217 ) = 0

где 09д0,11д2д3 - средние температуры характерных элементов исходной системы шлифовальной обработки соответственно: средняя температура стружки; средняя температура на передней поверхности зерна (температурой резания); средняя температура в теле режущего зерна; средняя температура на задней поверхности зерна; средняя температура на поверхности изделия в зоне резания.

кРх, кРу - коэффициенты резания по соответствующим координатам;

ЬРх, ЬРу, НРг - коэффициенты демпфирования резания по соответствующим координатам;

kx, ку, кг - коэффициенты, определяемые проекциями силы резания на соответствующие оси;

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

Ки1, Е^, - коэффициенты, учитывающие долю теплоты деформации,

переходящей в стружку, пропорционально тангенциальной составляющей силы резания, нормальной составляющей силы резания и скорости резания;

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

Оп - источник теплоты трения заготовки по задней поверхности;

1тп - определяет теплообмен между стружкой и передней поверхностью зерна;

1тз - определяет теплообмен между изделием и задней поверхностью зерна;

1с - отражает отвод тепла со стружкой из зоны резания;

1и - учитывает теплоотвод из зоны резания в глубину изделия;

1г, 1у - отражают теплообмен между задней и передней поверхностями зерна;

1к - учитывает отвод тепла от зерна в связку круга;

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

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

С помощью предложенной модели были рассчитаны температуры, возникающие в процессе плоского шлифования без охлаждения деталей из закаленной до ИЯС64 быстрорежущей стали периферией круга пШ00*20*75 с зернами эльбора Л063/50, связкой Б1, концентрацией зерна К0 = 100 %. Результаты были сопоставлены с экспериментальными данными, полученными

В.В. Щипановым в Тольяттинском политехническом институте, и с расчетными данными А.Н. Резника. Сравнительные данные представлены на рис. 6-8.

Глубина резания, мм

Рис. 6. Влияние глубины резания на температуру при шлифовании:

—- ряд 1;------------ ряд 2;-------ряд 3; - ряд 4; ж - ряд 5;------- ряд 6

Скорость резания, м/с

Рис. 7. Влияние скорости резания на температуру при шлифовании:

—- ряд 1;----------- ряд 2;-------ряд 3; - ряд 4; ж - ряд 5;------ ряд 6

Скорость заготовки, м/мин

Рис. 8. Влияние скорости заготовки на температуру при шлифовании:

—- ряд 1;--------- ряд 2;----- ряд 3; - ряд 4; ж - ряд 5;----- ряд 6

Здесь ряд 1 показывает изменение экспериментальной температуры, полученной с помощью закладной полуискусственной термопары. Ряды 2, 3,

4, 5 показывают изменения средних температур, полученных по модели (расчетным путем). Ряд 2 показывает температуру в зоне резания (стружки). Ряд 3 показывает среднюю температуру на поверхности передней поверхности зерен шлифовального круга (температура резания). Ряд 4 показывает изменение средних температур в зерне. Ряд 5 показывает среднюю температуру на задней поверхности зерна. Ряд 6 демонстрирует изменение температуры на поверхности заготовки.

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

Предложенная модель позволяет изучать влияние разнообразных факторов на средние температуры в процессе шлифования. Используя методику А.Н. Резника и полученные средние температуры в результате моделирования, можно описать температурное поле в поверхностных слоях заготовки

и определить, на какую глубину и с какой интенсивностью прогревается заготовка в процессе обработки.

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

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

Время, с

Рис. 9. Переходный процесс нарастания температуры резания при врезании

Если жесткость упругой системы недостаточна, то при врезании это может приводить к значительным колебаниям температуры резания (рис. 10), что может сказаться на качестве поверхностного слоя заготовки.

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

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

и

380

370

й 360 I 350

^ 340

5

| 330 | 320 310 300

^ЛЛЛЛЛА^Л^^У'АДЛ‘Ц\ллялмлгЛпА'|я/^У

5 10

Время, с Рис. 10. Переходный процесс

15

20

Время, с

Рис. 11. Переходный процесс колебаний температуры резания при ступенчатом возмущении на узел круга по оси У

0

и

о 5 Н о & £ £ * н & й

& сс <и о в с

^ 55 р22 § н X

3

эт

я

Время, с

Рис. 12. Переходный процесс нарастания температуры на задней поверхности зерна

при врезании

На рис. 13 приведена амплитудно-частотная характеристика колебаний температуры резания при периодическом изменении силы резания.

1400

и

я

СО

и

а

а

£

а

и

С

2

о

Н

1200

1000

800

600

400

200

0

0,0

10,0

20,0

50,0

60,0

70,0

30,0 40,0

Частота, 1/с

Рис. 13. АЧХ температуры резания от колебаний силы резания

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

вании, когда прерывистость процесса приводит к периодическим колебаниям силы резания.

Исследования показали, что возмущения в узлах шлифовального станка по соответствующим координатам резонируют на собственных частотах своих узлов и вызывают повышенные колебания температуры резания (рис. 14, 15).

О

о

Й

&

с

2

Е-

Частота, 1/с

Рис. 14. АЧХ колебаний температуры резания при шлифовании от периодического возмущения на узел круга по оси У

Данные процессы вызывают также резонансные колебания температуры на задней поверхности зерна (рис. 15).

0,035

О

к

*

СО

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

а

й

а

£

&

4>

С

5

£

0,030

0,025

0,020

0,015

0,010

0,005

0,000

100

200 300

Частота, 1/с

400

500

Рис. 15. АЧХ колебаний температуры резания от возмущений в электродвигателя

привода шлифовального круга

0

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

10,0

К 8,0

Л

6,0

и

4,0

2,0

0,0 -I

20

40 60

Частота, 1/с

80

100

Рис. 16. АЧХ колебаний силы резания от воздействий на узел круга по оси У

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

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

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

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

0

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

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

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

1. Аршинов В. А., Алексеев Г. Л. Резание металлов и режущий инструмент. - М.: Машиностроение, 1976. - 440 с.

2. Зарубин В.С., Крищенко А.П. Математическое моделирование в технике: учебник для вузов. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2001. - 496 с.

3. Кедров С.С. Колебания металлорежущих станков. - М.: Машиностроение, 1978. - 199 с.

4. Кудинов В. А. Динамика станков. - М.: Машиностроение, 1967. - 360 с.

5. Кудинов В. А. Системность и диалектика в динамике технологических процессов и машин. - СТИН. - 2000. - № 1. - С. 5-11.

6. Резников А.Н. Теплофизика процессов механической обработки материалов. - М.: Машиностроение, 1981. - 279 с.

7. Резников А.Н. Теплофизика резания. - М.: Машиностроение, 1969. -

288 с.

8. Резников А.Н., Резников Л. А. Тепловые процессы в технологических системах: учебник для вузов по специальностям «Технология машиностроения», «Металлорежущие станки и инструменты». - М.: Машиностроение, 1990. - 288 с.

9. Якимов А.В., Слободяник П.Т., Усов А.В. Теплофизика механической обработки: учеб. пособие. - Киев; Одесса: Лыбидь, 1991. - 240 с.

10. Ящерицын П.И. Теория резания. Физические и тепловые процессы в технологических системах: учебник для вузов. - Минск: Выш. шк., 1990. -512 с.

Получено 20.04.2010

i Надоели баннеры? Вы всегда можете отключить рекламу.