Научная статья на тему 'Исследование эффекта Баушингера и границы текучести при упругопластическом деформировании металлов'

Исследование эффекта Баушингера и границы текучести при упругопластическом деформировании металлов Текст научной статьи по специальности «Физика»

CC BY
464
261
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ПЛАСТИЧНОСТЬ / УПРУГОСТЬ / СЛОЖНОЕ НАГРУЖЕНИЕ / ЭФФЕКТ БАУШИНГЕРА / ПОВЕРХНОСТЬ ТЕКУЧЕСТИ / ПРОЦЕССЫ ДЕФОРМИРОВАНИЯ / PLASTICITY / ELASTICITY / COMPLEX LOADING / BAUSHINGER EFFECT / THE YIELD SURFACE / THE DEFORMATION PROCESSES

Аннотация научной статьи по физике, автор научной работы — Зубчанинов Владимир Георгиевич, Алексеев Андрей Алексеевич, Гультяев Вадим Иванович

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

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

Experimental study of Baushinger effect and yield surface at elastoplastic deformation of metals

This work contains results of the experimental investigation on repeated sign-variable loading at stretching compression of a thin-walled tubular specimen from steel 45 in an automatized computational and experimental installation SN-EVM. Аlso presents the results of the experimental assessment of the Baushinger effect and the influence of various admissions on residual deformation on it. The influence of the admission on residual deformation on definition of secondary offset yield stress was investigated, as well as on initial and subsequent radius function and position of the center of a hypothetical spherical yield surface in stress space in the plastic-flow theory. According to the results of the investigation, the parameter characterizing the Bauschinger and the radius function of the yield surface increase on the increase of admission on residual deformation effect, yet the displacement of the center of the yield surface decreases. It is shown that the parameter characterizing the Bauschinger effect decreases on the increase of the length of the arc of plastic deformation and tends to a stationary value. It is experimentally discovered that for steel 45 the displacement of the center of the yield surface increases on the increase of the length of the arc of plastic deformation. It is also established that the radius of the spherical yield surface lowers temporary, and then increases on the increase of the length of the arc of plastic deformation radius. Thus, in some mathematical models of plastic-flow theory temporary decrease of the radius of the yield surface (contraction of the yield surface) at the beginning of the process of plastic deformation is not associated with the definition of the yield stress for the admission on residual deformation.

Текст научной работы на тему «Исследование эффекта Баушингера и границы текучести при упругопластическом деформировании металлов»

2013 Механика № 1

УДК 539.3

В.Г. Зубчанинов, А.А. Алексеев, В.И. Гультяев

Тверской государственный технический университет, Тверь, Россия

ИССЛЕДОВАНИЕ ЭФФЕКТА БАУШИНГЕРА И ГРАНИЦЫ ТЕКУЧЕСТИ ПРИ УПРУГОПЛАСТИЧЕСКОМ ДЕФОРМИРОВАНИИ МЕТАЛЛОВ

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

Ключевые слова: пластичность, упругость, сложное нагружение, эффект Баушингера, поверхность текучести, процессы деформирования.

V.G. Zubchaninov, A.A. Alekseev, V.I. Gultyaev

Tver State Technical University, Tver, Russian Federation

EXPERIMENTAL STUDY OF BAUSHINGER EFFECT AND YIELD SURFACE AT ELASTOPLASTIC DEFORMATION OF METALS

This work contains results of the experimental investigation on repeated sign-variable loading at stretching compression of a thin-walled tubular specimen from steel 45 in an automatized computational and experimental installation SN-EVM. Also presents the results of the experimental assessment of the

Baushinger effect and the influence of various admissions on residual deformation on it. The influence of the admission on residual deformation on definition of secondary offset yield stress was investigated, as well as on initial and subsequent radius function and position of the center of a hypothetical spherical yield surface in stress space in the plastic-flow theory. According to the results of the investigation, the parameter characterizing the Bauschinger and the radius function of the yield surface increase on the increase of admission on residual deformation effect, yet the displacement of the center of the yield surface decreases. It is shown that the parameter characterizing the Bauschinger effect decreases on the increase of the length of the arc of plastic deformation and tends to a stationary value. It is experimentally discovered that for steel 45 the displacement of the center of the yield surface increases on the increase of the length of the arc of plastic deformation. It is also established that the radius of the spherical yield surface lowers temporary, and then increases on the increase of the length of the arc of plastic deformation radius. Thus, in some mathematical models of plastic-flow theory temporary decrease of the radius of the yield surface (contraction of the yield surface) at the beginning of the process of plastic deformation is not associated with the definition of the yield stress for the admission on residual deformation.

Keywords: plasticity, elasticity, complex loading, Baushinger effect, the yield surface, the deformation processes.

В теории пластичности тензоры напряжений и деформаций представляют [1-5] в виде

T = (^) = ао(5&.) + (Sj), Т£ = (s.) = So(5..) + (Эу), (i,j = 1,2,3), (1) где 5 j - символ Кронекера;

Сто = a ii /3, So = Sii /3, ( i =1 2,3) (2)

- средние напряжение и деформация (компоненты шаровых тензоров);

Sij = °ij - 5 jao, 3ij = &ij - 5 jso , (h j =1 2, 3) (3)

- компоненты тензоров-девиаторов напряжений и деформаций соответственно.

Тензорам напряжений и деформаций (1) в линейном тензорнокоординатном евклидовом пространстве Е6 поставлены в соответствие векторы напряжений и деформаций [1-5]:

S = S° + о, £ = £° + Э, (4)

где

S° = Soi°, Ё° = 3oi° (5)

- векторы напряжений и деформаций объемного растяжения

и сжатия;

о = Skik , Э = Экik (k = 1,2,...5) (6)

- векторы напряжений и деформаций формоизменения в пятимерном девиаторном подпространстве Е5; |іА | - ортонормированный

фиксированный базис А.А.Ильюшина тензорно-координатного подпространства Е5;

SoW3ao, S, =J-SU, S2 (S22-S33), S3^a/2S,2, S4=^3, S5=л/25п,

V V2 (7)

3o =V3eo, Э1 2Э.1, Э2 (Э22-Э33), Э3^¡2Эп, Э4W^, Э5 =л/2Эп

- компоненты векторов напряжений и деформаций.

Модули векторов напряжений и деформаций в пятимерном подпространстве формоизменения Е5 равны модулям тензоров-девиаторов:

(k = 1,2, ...5), (i, j = 1,2,3).

В теории пластического течения вводятся две основополагающие гипотезы [3, 6]. Первая гипотеза состоит в разложении тензора-

девиатора деформаций на упругие s j и пластические sp части, что позволяет ввести понятие о мгновенной гипотетической поверхности текучести f (о, Эр) = o, разделяющей в подпространстве формоизменения Е5 область активного пластического деформирования и область

упругой разгрузки [3, 6] (рис. 1). Принцип градиентальности Драккера позволяет определить приращения векторов упругих и пластических деформаций [6] для активного и пассивного процессов деформирования соответственно.

d3e = , d3p = dA, grad f, ( do • grad f > o),

- d" - (9)

d3e = —, d3p = o, (do • grad f < o).

2G

Вторая гипотеза относится к возможности разложения полного вектора напряжений (рис. 1) [6-8]:

о = о0 + a, (Ю)

где а0 - вектор активных напряжений; а - вектор добавочных остаточных микронапряжений в Е5.

аэрА

Рис. 1. Разложение полного вектора напряжений в Е5

Математические модели (определяющие соотношения) теории течения отличаются формой поверхности текучести. В наиболее распространенном на практике варианте теории с трансляционно-изотропным упрочнением материала она имеет форму сферы:

2/ = • о° -С2р(5р) = 0, (11)

которая может изменять свои размеры и местоположение [6-8]. Здесь ст° = Ср (ур) - скалярная функция изотропного упрочнения, равная радиусу поверхности текучести; бр - длина дуги траектории пластического деформирования. В начальном состоянии при бр = 0 радиус гипотетической начальной поверхности текучести ст° = стт = V2/3стт (см. рис. 1), где стт - начальный предел текучести при простом нагружении, определяемый по техническому допуску на остаточную деформацию Эр » ^3/2 ер1 = 0,245 % (£р = 0,2 %). Принимается универсальный закон упрочнения Одквиста-Ильюшина

ст = Ф(^), (12)

мало отличающийся от закона единой кривой при простом нагружении а = Ф(Э). За новый предел текучести принимается точка К начала разгрузки (12).

Эффект Баушингера при знакопеременном нагружении в Е5 оценивается [9] безразмерным параметром

У

а т

м

/ <5ТК (13)

при соответствующей величине дуги траектории пластического деформирования бр , где стК _ новый предел текучести на диаграмме

растяжения в некоторой точке К начала разгрузки; аТм - вторичный

предел текучести при разгрузке из той же точки К при «протыкании» поверхности текучести по диаметральному направлению в результате излома траектории на 180° в Е5, определяемый по допуску на остаточную деформацию Эр. Радиус а0 = Ср гипотетической текущей

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

а0 = ср = СТк , (14)

2

а смещение ее центра

Т Т

о V- + с)

_____Т п _ ^К^ /1СЧ

а -Ок Ср ~ 2 . ( )

Экспериментальное исследование функций ст0 = Ср (бр ) и а(Бр)

проводилось в опыте при знакопеременном нагружении-разгужении на автоматизированном испытательном комплексе на сложное нагружение СН-ЭВМ им. А.А.Ильюшина в лаборатории механических испытаний кафедры «Сопротивление материалов, теория упругости и пластичности» Тверского государственного технического университета. В опыте использовался тонкостенный трубчатый образец из стали 45 в состоянии поставки с площадкой текучести, который имел толщину стенки к = 1 мм, радиус срединной поверхности поперечного сечения Я = 15,5 мм и длину рабочей части I = 110 мм. В результате испытания образец подвергался многократному знакопеременному нагружению через равные приращения АЭ1 = 0,5% при растяжении (рис. 2, 3), с последующим знакопеременным нагружением на |ДЭ| » 0,75 - 0,9% при сжатии. Опыт в режиме непрерывного деформирования продолжался

в течение 8 часов. На рис. 4 представлены полученные опытные зависимости параметра у от длины дуги траектории пластического деформирования ^, характеризующие эффект Баушингера, на рис. 5 - изменение радиуса ст° гипотетической сферической поверхности текучести, а на рис. 6 - график смещения центра гипотетической

сферической поверхности текучести в зависимости от параметра бр . Все представленные зависимости построены при различных допусках

на остаточную деформацию е*7 = 0,2; 0,1; 0,05; 0,025%, что в Е5 составляет Э*7 = 0,245; 0,1225; 0,06125; 0,030625% соответственно.

5„ МПа

r7i

/ / / /

f h /1

// / / Э], %

5/ В / § 1 i h S1''! 5 1 41

500 400

300 200 100

0

-100 -200 -300

Рис. 2. Локальная диаграмма знакопеременного нагружения 51 -

а, МПа

450 400 350 300 250 200 150 100 50 0

r°^7, #j /J о

’ c? ^ у S' 8 9 0 f

ГТТЩЛ /1 / i ... _. У $ 9 6 oP } R & 8 XQ 8 9 f

/ f / iff / 2 2 9 ? 6 X X ^ 8 о о ?

Щм ij IjJf 1 1 / < о $ ? $ г * a у i , ? о 6 } 6 f П

|i|l щ! mil jL1 11 6 d 9 8 11

Wmk Im In j\t nt H фЧц я S 4vS ы » 8 $ 8 о 8 9

ШщЩУ'Тvlr Я о // Э, %

0,5 I 1,5 2 2,5 3 3,5 4 4,5

Рис. 3. Глобальная диаграмма деформирования а-Э

Рис. 4. Эффект относительного изменения предела текучести по Баушингеру

250

200

150

100

50

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

Рис. 5. Изменение радиуса гипотетической поверхности текучести

а, МПа

допуск £ допуск £ допуск е допуск 1= 0,2 О/ £=0,1 °/ £= 0,05 £= 0,021 )

>%

0 0,5 1 1,5 2 2,5 3 3,5 4

%

Рис. 6. Смещение центра гипотетической поверхности текучести

С ростом допуска на остаточную деформацию радиус поверхности текучести а0 увеличивается, а смещение ее центра а уменьшается. Максимальное отклонение для рассмотренных крайних допусков на остаточную деформацию е*7 = 0,2% и е*7 = 0,025% для радиуса а0

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

За технический предел текучести стт в теории пластичности принимается такое напряжение при растяжении, при котором остаточные

деформации становятся одного порядка с упругими (е «10 ). За такую

остаточную деформацию, как правило, принимают е*7 = 0,2% = 2 -10_3

(Э*7 = 0,245%). Начальная поверхность текучести в девиаторном пространстве А.А. Ильюшина при данном допуске на остаточную деформацию описывается сферой Мизеса для начально изотропных тел. При меньших допусках очертание сферы Мизеса искажается и теряет свою форму [8]. На девиаторной плоскости окружность Мизеса и вписанный в нее шестиугольник Сен-Венана приобретают тройную симметрию [11]. При этом начальные условия текучести Мизеса и Сен-Венана определяются формулами [3-5]

л/3

1 л/(°"1 -а2 )2 + (а2 “а3 )2 + (СТ3 1 )2 =^Т

(16)

тшеХ = °т 2 °п = к (т < п; т п = 1,2,3)

где стг- (I = 1, 2,3) - главные нормальные напряжения; к - предел текучести при плоском чистом сдвиге. Уже при е*7 = 0,1% классическая

теория пластичности лежит за пределами ее инженерного контроля [3, 6, 8]. К сожалению, в некоторых математических моделях теории течения величине допуска на остаточную деформацию не придают значения и не связывают их с определением предела текучести. В про-

цессе нагружения при трансляции предельной поверхности она вытягивается в направлении развития процесса вследствие развития деформационной анизотропии [3, 6].

По результатам экспериментально проведенного исследования можно сделать выводы.

1. Эффект Баушингера для стали 45 с площадкой выражается

уменьшением по модулю вторичного предела текучести стТм и параметра у с ростом длины дуги пластического деформирования sp . Параметр у при различных допусках на остаточную деформацию стремится к некоторому стационарному значению при sp > 3 % .

2. С ростом допуска на остаточную деформацию Эр параметр у увеличивается и при максимально достигнутом значении sp для общепринятого допуска sp = 0,2 % (Эр = 0,245 %) составляет примерно 0,3.

3. Радиус гипотетической сферической поверхности ст0 = Ср (sp) совершает временное понижение типа «нырка», а затем увеличивается. При повышении допуска Эр примерно до технического и более очертание поверхности стремится к сферической. При этом отклонение радиуса при крайних значениях допуска на остаточную деформацию достигает 45 %.

4. Отклонения параметра смещения центра предельной поверхности a(sp) для рассмотренных различных допусков на остаточную деформацию при определении пределов текучести достигают 55 %.

5. В некоторых математических моделях теории течения временное убывание функции ст0 = Ср (sp) в начале процесса пластического

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

Библиографический список

1. Ильюшин A.A. Пластичность. Основы общей математической теории. - М.: Изд-во АН СССР, 1963. - 273 с.

2. Ильюшин A.A. Механика сплошной среды. - М.: Изд-во МГУ, 1990. - 310 с.

3. Зубчанинов В.Г. Механика процессов пластических сред. - М.: Физматлит, 2010. - 352 с.

4. Зубчанинов В.Г. Математическая теория пластичности. -Тверь: Изд-во Твер. гос. техн. ун-та, 2002. - 300 с.

5. Зубчанинов В.Г. Устойчивость и пластичность. Пластичность. -М.: Физматлит, 2008. - Т. 2. - 336 с.

6. Зубчанинов В.Г. О соотношениях между напряжениями и деформациями в теории пластичности при сложном нагружении // Проблемы прочности и пластичности: межвуз. сб. ННГУ. - Н. Новгород, 2011. - № 73. - С. 120-131.

7. Поль Б. Макроскопические критерии пластического течения и хрупкого разрушения // Разрушение. Т.2.: Математические основы теории разрушения / под ред. Г. Либовица. - М. : Мир, 1975. - С. 336-520.

8. Новожилов В.В. Вопросы механики сплошных сред. - Л.: Судостроение, 1989. - 397 с.

9. Москвитин В.В. Пластичность при переменных нагружениях. -М.: Изд-во МГУ, 1965. - 264 с.

10. Зубчанинов В.Г., Гультяев В.И., Алексеев А.А. Об эффекте Баушингера и поверхности текучести при пластическом деформировании металлов // Вестн. ЧГПУ им. И.Я. Яковлева. Сер. Механика предельного состояния. - Йошкар-Ола, 2012. - № 3 (13). - С. 3-8.

11. Ишлинский А.Ю., Ивлев Д.Д. Математическая теория пластичности. - М.: Физматлит, 2001. - 704 с.

References

1. Ilyushin A.A. Plastichnost. Osnovy obshchey matematicheskoy te-orii [Plasticity. Bases of the General Mathematical Theory]. Moscow: Akademiya nauk USSR, 1963. 273 p.

2. Ilyushin A.A. Mekhanika sploshnoy sredy [Continuum Mechanics]. Moscow: Moskowskiy gosudarstvenniy universitet, 1990, 310 p.

3. Zubchaninov V.G. Mekhanika protsessov plasticheskikh sred [Mechanics of processes in plastic environments]. Moscow: Fizmatlit, 2010, 352 p.

4. Zubchaninov V.G. Matematicheskaya teoriya plastichnosti [The mathematical theory of plasticity]. Tver: Tverskoy gosudarstvenniy tekhni-cheskiy universitet, 2002. 300 p.

5. Zubchaninov V.G. Ustoychivost i plastichnost. Tom 2. Plastichnost [Stability and plasticity. Vol. 2: Plasticity]. Moscow: Fizmatlit, 2008, 336 p.

6. Zubchaninov V.G. O sootnosheniyakh mezhdu napryazheniyami i deformatsiyami v teorii plastichnosti pri slozhnom nagruzhenii [About relations between stress and strain in the plasticity theory at complex loading]. Problemy prochnosti i plastichnosti: mezhvuzovskiy sbornik Nizhegorod-skogo gosudarstvennogo universiteta, 2011, no. 73, pp. 120-131.

7. Paul B. Makroskopicheskiye kriterii plasticheskogo techeniya i khrupkogo razrusheniya [Macroscopic criteria for plastic flow and brittle fracture]. Razrusheniye. Tom 2. Matematicheskiye osnovy teorii razrusheniya / ed. H. Liebowitz. Moscow: Mir, 1975, pp. 336-520.

8. Novozhilov V.V. Voprosy mekhaniki sploshnykh sred [Questions of continuum mechanics]. St. Petersbourg: Sudostroyeniye, 1989, 397 p.

9. Moskvitin V.V. Plastichnost pri peremennykh nagruzheniyakh [Plasticity under alternating loadings]. Moscow: Moskowskiy gosudarstvenniy universitet, 1965, 264 p.

10. Zubchaninov V.G., Gultyayev V.I., Alekseev A.A. Ob effekte Baushingera i poverkhnosti tekuchesti pri plasticheskom deformirovanii metallov [About Baushinger’s effect ant the yield surface at plastic deformation of metals]. Vestnik Chuvashskogo gosudarstvennogo peda-gogicheskogo universiteta. Mekhanika predelnogo sostoyaniya, Cheboksary, 2012, no. 3 (13), pp. 3-8.

11. Ishlinskiy A.Yu., Ivlev D.D. Matematicheskaya teoriya plastichnosti [The mathematical theory of plasticity]. Moscow: Fizmatlit, 2001. 704 p.

Об авторах

Зубчанинов Владимир Георгиевич (Тверь, Россия) - доктор технических наук, профессор, профессор кафедры сопротивления материалов, теории упругости и пластичности Тверского государственного технического университета (170026, г. Тверь, Наб. А. Никитина, 22, e-mail: vgz@rambler.ru).

Алексеев Андрей Алексеевич (Тверь, Россия) - кандидат технических наук, доцент, доцент кафедры сопротивления материалов, теории упругости и пластичности Тверского государственного технического университета (170026, г. Тверь, Наб. А. Никитина, 22, e-mail: alexeew@bk.ru).

Гультяев Вадим Иванович (Тверь, Россия) - доктор технических наук, доцент, заведующий кафедрой строительства и энергетики Тверского государственного технического унивверситета (170026, г. Тверь, Наб. А. Никитина, 22, e-mail: vig0@mail.ru).

About the authors

Zubchaninov Vladimir Georgievich (Tver, Russian Federation) -Doctor of Technical Sciences, Professor, Department of the Resistence of Materials, Elasticity and Plasticity Theories, Tver State Technical University (22, Nab. A.Nikitina, 170026, Tver, Russian Federation, e-mail: vgz@rambler.ru).

Alekseev Andrey Alekseevich (Tver, Russian Federation) - Ph. D. in Technical Sciences, Ass. Professor, Department of the Resistence of Materials, Elasticity and Plasticity Theories, Tver State Technical University (22, Nab. A.Nikitina, 170026, Tver, Russian Federation, e-mail: alexeew@bk.ru).

Gultyaev Vadim Ivanovich (Tver, Russian Federation) - Doctor of Technical Sciences, Ass. Professor, Head of the Department of the Building and Energy, Tver State Technical University (22, Nab. A.Nikitina, 170026, Tver, Russian Federation, e-mail: vig0@mail.ru).

Получено 15.02.2013

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