Курсовая: Моделирование процессов переработки пластмасс - текст курсовой. Скачать бесплатно.
Банк рефератов, курсовых и дипломных работ. Много и бесплатно. # | Правила оформления работ | Добавить в избранное
 
 
   
Меню Меню Меню Меню Меню
   
Napishem.com Napishem.com Napishem.com

Курсовая

Моделирование процессов переработки пластмасс

Банк рефератов / Химия

Рубрики  Рубрики реферат банка

закрыть
Категория: Курсовая работа
Язык курсовой: Русский
Дата добавления:   
 
Скачать
Архив Zip, 850 kb, скачать бесплатно
Заказать
Узнать стоимость написания уникальной курсовой работы

Узнайте стоимость написания уникальной работы

20 Министерство образования Республики Беларусь Учреждение образования: “Белорусский государственный технологический университет” Кафедра автоматизации про изводственных процессов и электротехники Расчётно-пояснительная записка К курсовому проекту по курсу применения ЭВМ в химической промышленности на тему: Моделирование процессов переработки пластмасс Разработал: студент Факультета ТОВ 4к. 1 гр. Кардаш А. В. Проверил: Овсянников А . В . Минск 200 4 РЕФЕРАТ Данна я курсовая работа содержит 26 листов печатного текста, 7 рисунков, 6 6 формул. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ, ДИФЕРИНЦИАЛЬНОЕ УРАВНЕНИЕ, ТЕПЛОПРОВОДНОСТЬ, ВРЕМЯ, ЛИТНИКОВЫЙ КАНАЛ, ОХЛАЖДЕНИЕ , ТЕМПЕРАТУРНОЕ ПОЛЕ. Курсовая работа содержит расчет температурного поля литникового канала литьевой формы , теоретические сведения о процессах происходящих в химической технологии связанных с охлаждением и нагреванием материалов , построение математической модели описывающую теплообмен между бесконечно-длинным цилиндром и его поверхностью , описание переменных входящих в модель. Разработана программа описывающая охлаждение полистирольного литника формы. СОДЕРЖАНИЕ РЕФЕРАТ 2 СОДЕРЖАНИЕ 3 ВВЕДЕНИЕ 4 1. АНАЛИЗ ИСХОДНЫХ ДАННЫХ 5 1.1 Неограниченный цилиндр. 5 1.2 Описание переменных 5 1.3 Граничные условия 5 2 ОБЩИЕ ТЕОРЕТИЧЕСКИЕ СВЕДЕНИЯ 6 2.1 Теплообмен 6 2.1.1 Теплопроводность 6 2.1.2. Теплопередача в стационарном режиме. 7 2.1.3. Нестационарная теплопроводность. 7 2.2. Нагревание и охлаждение тел простой геометрической формы 8 2.2.1. Плоская неограниченная пластина. 8 2.2.2 Неограниченный цилиндр. 10 2.3. Теплопроводность в процессах, сопровождающихся изменением физического состояния 11 2.3.1. Плавление в области х > 0. 12 2.3.2. Затвердевание. 12 2.3.3 Плавление с непрерывным удалением расплава. 13 2.4.Теплопередача в потоках расплава 13 2.5. Лучистый теплообмен 15 3. СОСТАВЛЕНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ИССЛЕДУЕМОГО ПРОЦЕССА. 17 3.1. Специфика построения математических моделей описывающих термодинамические процессы 17 3.2. Вывод дифференциального уравнения теплопроводности. 17 4 СОСТАВЛЕНИЕ АЛГОРИТМА 20 5 СОСТАВЛЕНИЕ ПРОГРАММЫ 22 6 АНАЛИЗ МОДЕЛИРОВАНИЯ И РАСЧЁТОВ 24 СПИСОК ИСПОЛЬЗОВАНОЙ ЛИТЕРАТУРЫ 25 ПРИЛОЖЕНИЕ1 26 ПРИЛОЖЕНИЕ2 27 ВВЕДЕНИЕ Переработка полимерных материалов — это совокупность техноло гических приемов, методов и процессов, посредством которых ис ходный полимер превращают в различные изделия с заданными эксплуатационными характеристиками. Полимеры начали перерабатывать в конце XIX в., а к сере дине XX в. переработка полимеров выделилась в самостоятельную область техники, в которой используется специализированное вы сокопроизводительное оборудование, необходимое для реализации в промышленных масштабах специфических для полимеров техно логических процессов. Вследствие большой производительности современного перера батывающего оборудования и высокой стоимости технологических линий проведение экспериментальных исследований реального про цесса переработки полимеров, даже осуществленных с примене нием современных методов экстремального планирования, пре вращается в дорогостоящую и продолжительную работу. Поэтому целесообразно изучать особенность каждого конкретного процесса, рассматривая вначале его теоретическое описание, т. е. его мате матическую модель. При таком подходе в каждом конкретном случае этапу физи ческого эксперимента (будь то создание несложной установки, конструирование технологической линии или опробование нового технологического режима) всегда предшествует этап теоретиче ского эксперимента. На этом этапе нет необходимости прибегать к реальным экспериментам, вместо этого исследуются количествен ные характеристики процесса, полученные расчетным методом. Такой подход позволяет существенно снизить объем физиче ского эксперимента, поскольку прибегать к нему приходится на самой последней стадии — не в процессе поиска основных законо мерностей, а для проверки и уточнения выданных рекомендаций. Разумеется, для того чтобы исследуемые теоретические модели процессов описывали эти процессы с достаточно хорошим прибли жением, они непременно должны учитывать основные особенно сти моделируемых явлении. П ри математическом описании реальных производственных процессов приходится прибегать к существенным упрощениям. При этом значительную помощь в создании математических моделей оказывает анализ простых слу чаев . Прием такого рода вполне допустим , о н позволяет независимо устанавливать основ ные закономерности наиболее простых случаев выбранных в качестве математического аналога поведения полимерных расплавов. Термодинамические соотношения, описывающие разогрев и плавление полимеров, являются фундаментом, на базе которого строятся неизотермические модели реальных процессов перера ботки. Основные вопросы термодинамики и теплопередачи в поли мерах рассмотрены в данной работе. 1. АНАЛИЗ ИСХОДНЫХ ДАННЫХ 1.1 Неограниченный цилиндр. Рассмотрим неограниченный цилиндр радиуса R , температура поверхности которого остается неизменной на протяжении всего процесса теплообмена. Радиальное распределение температур в начальный момент задано в виде некоторой функции Т( r ) . Необходимо найти распределение температур . Такие задачи встречаются при расчете процессов охлаждения полимерного волокна, затвердевания литников литьевых форм и т. п. Дифференциальное уравнение теплопроводности для цилиндра имеет вид: (1.1) Краевые условия: (1.2) (1.3) (1.4) Решение, полученное методом разделения переменных , имеет сложный вид потому задачей данной работы является найти численное его решение. 1.2 Описание переменных Уравнение теплопроводности устанавливает зависимость между следующими величинами характеризующими процесс теплопроводности : T - температура по Цельсию ( градус ) r - радиус цилиндра (М) t - время (С) a -коэффициент температуропроводности ( градус / с* м2 ) 2 1.3 Граничные условия Для решения данного дифференциального уравнения в частных производных необходимыми данными является значения производных температуры по радиусу на оси цилиндра, которая должна быть равной нулю (1.4). Температуру стенки цилиндра, через которую происходит охлаждение литника примем равной 30 градусов. (1.5) Радиус литника обычно составляет 0.01 м. R =0.01 (1.6) Распределение температуры в начальный момент времени по радиусу задано в виде убывающей экспоненциальной функции , чтобы производная температуры по времени на оси цилиндра была равной нулю , радиус возводим в квадрат (1 . 7) (1.7) 2 ОБЩИЕ ТЕОРЕТИЧЕСКИЕ СВЕДЕНИЯ 2.1 Т еплообмен Различают три вида теплообмена: теплопроводность, теплопередача конвекцией и лучистый теплообмен. Передача тепла за счет теплопроводности осуществляется в результате движения молекул, атомов и электронов; она играет значительную роль при теплообмене в твердых и расплавленных полимерах. При конвекции, которая возможна только в жидкостях и газах, тепло передается за счет относительного движения частиц нагретого тела. При лучистом теплообмене передача тепла между пространственно разделенными частями тела происходит за счет электромагнитного излучения. 2.1.1 Т еплопроводность Основной задачей теории теплопроводности является установление распределения температур внутри тела. Если распределение температур не зависит от времени, то задача теплопроводности является стационарной; если распределение температур зависит от времени, то задача становится нестационарной. Передача тепла происходит в о всех случаях, когда в теле су ществует температурный градиент. По закону Фурье, который лежит в основе всех расчетов теплопроводности, для изотропных мат ериалов вектор теплового потока q пропорционален температурному градиенту: (2.1) где q — количество тепла, проходящего через единичную поверхность, перпен­дикулярную направлению теплового потока; k — коэффициент теплопроводности. Полагая в уравнении энергетического баланса V = О, получим: (2.2) Уравнение (2.2) предст авляет собой уравнение теплопро водности для изотропного твердого тела. Если внутри изотропного тела имеется источник тепла, то уравнение ( 2.2 ) необходимо до полнить членом, учитывающим теп ловыделение (2.3) где — коэффициент температуропроводности [замена на в уравнении ( 2 . 3 ) возможна для несжимаемых твердых тел]; — оператор Лапласа в прямоугольной системе координат (2.4) G — интенсивность внутренних тепловыделений, отнесенная к единице объема . П римерами внутренних те пловыделений являются поглощения инфракрасного из лучения в полупрозрачных средах, экзотер мический эффект химических реак ций и т. п. 2. 1 .2 . Теплопередача в стационарном режиме. Теплопередачу в непрерывно действующих нагревательных системах перерабатывающего оборудования можно рассматривать как независящую от времени. Следовательно, распределение температур носит установившийся характер и определяется интегрированием дифференциального уравнения (2 .5) (2.5) 2.1.3. Нестационарная теплопроводность. В большинстве случаев в реальных процессах переработки приходится иметь дело с нестационарным режимом теплопроводности, когда полимер подвергают нагреву или охлаждению (например, охлаждение в форме отлитого изделия). Теоретические исследования процесса нестационарной теплопроводности представляют собой обширный раздел математической физики. Решения, получаемые в результате интегрирования уравнения (2.5), представляют собой функции времени и пространственных координат, удовлетворяющие начальным и граничным условиям. Различают четыре рода граничных условий Условия первого рода: задано распределение температур на поверхности, которое может либо быть постоянным, либо зависеть от времени; в простейшем случае, если положение границ определяется одним числом (например, расстоянием L), такие граничные условия математически определяются выражением вида (2.6): (2.6) Условия второго рода: задана плотность теплового потока для каждой точки поверхности тела как функция времени: (2.7) Условия третьего рода: задан коэффициент теплообмена, а на границе и температура контактирующей с граничной поверхностью среды: (2 . 8) Условия четвертого рода: соответствуют теплообмену тела с окружающей средой по закону теплопроводности или теплообмену системы тел, находящихся в тепловом контакте (температура соприкасающихся поверхностей одинакова): (2.9) (2.10) Аналитическая теория нестационарной теплопроводности располагает большим набором решений одномерных задач, к которым принято сводить все многообра зие задач, встречающихся в инже нерной практике. В настоящее время получены аналитические ре шения для теплопроводности в пл оской стенке, в цилиндре, в кор пусе и в сфере. 2.2. Н агревание и охлаждение тел простой геометрической формы 2.2.1. Плоская неограниченная пластина. Под неограниченной обычно понимают такую пластину, ширина и длина которой во много раз превышают толщину. Таким образом, неограниченная пластина (рис. 2 . 1 ) представляет собой тело, ограниченное двумя параллельными плоскостями. Изменение температуры происходит только в одном направлении (х), в двух других направлениях (у и z) температура неизменна. Рис. 2.1. Положение координат при ис следовании теплового процесса в неограниченной пластине. Следовательно, задача является одномерной. Для одномерного теплового потока без внутреннего источника тепла уравнение теплопроводности сводится к виду: (2.11) Обычно используют граничные условия третьего рода: (2.12) Рассмотрим случай, когда в начальный момент температура пластины во всех точках была одинакова и равна То. Это начальное условие записывается в виде: (2.13) Решение, полученное методом преобразования Лапласа, имеет вид: (2.14) Здесь — безразмерная температура; — критерий Фурье (критерий гомохронности для процессов чистой теплопроводности ); - безразмерная координата; — функ ция ошибок, где ; Если коэффициент теплоотдачи очень велик (это эквивалентно заданию постоянной температуры на стенке), уравнение (2.14) упрощается: (2.15) Д ля прикидочных расчетов удобно пользоваться номограм мой зависимости от представленной на рис. 2.2 Рис.2.2 Номограмма для определения безразмеоной температуры в сечении неограниченной пластины при Если значение критерия Фурье велико, но не равно бесконечности, решение имеет вид: (2.16) Здесь (2.17) где — корни характеристического уравнения (2.18) где Bi = aw / — критерий Био. Уравнение (2.18) имеет бесчисленное множество действительных положительных корней. Первые пять корней для различных значений критер и я Био были вычислены Карслоу и Егером . Обычно на практике поль зуются номограммами. Номограмма позволяющая определить безразмерную температуру п ри различных значениях критерях Био приведена на рис . 2.3 Ри с. 2.3 Номограмма для определения безразмерной температуры по верхности неограниченной пластины . Ана логичная номограмма, предназ наченная для определения тем пературы в цент ре пластины, при ведена на рис. 2.4 . Р ис. 2.4 Н омограмма для определения безразмерной температуры в середине неограниченной пластины 2.2.2 Неограниченный цилиндр. Рас смотрим неограниченный цилиндр радиуса R , температура поверх ности которого остается неизмен ной на протяжении всего процес са теплообмена. Радиальное рас пределение температур в началь ный момент задано в виде некоторой функции Т( r ) . Необходимо найти распределение температур определения в цилиндре в любой момент вре мени. З адачи такого типа встречаются при расчете процессов охлаждения полимерного волокна, затвердевания литников литьевых форм и т. п. Дифференциальное уравнение теплопроводности для цилиндра имеет вид: (2.19) Краевые условия: Решение, полученное методом разделения переменных , в без размерной форме, имеет вид: (2.20) Для оценки изменения теплосодержания цилиндра определим среднюю температуру как: (2.21) Тогда безразмерная средняя температура определится соотноше нием: (2.22) где ; - корни функции Бесселя первого рода нулевого порядка определяемые выражением: (2.23) Таким образом, уменьшение средней температуры описывается простым экспоненциальным законом. Для удобства прикидочных расчетов на рис. IV . 10 приведен а номограмма зависимости между и Fo . 2.3. Т еплопроводность в процессах, сопровождающихся изменением физического состояния Анализируя процессы переработки полимеров, часто приходится встречаться с задачей о нагреве или охлаждении полимера, сопровождающемся изменением физического состояния ( плавлением или затвердением). Теоретическое рассмотрение задач такого типа впервые выполнено Нейманном. Мы остановимся только на одном, наиболее простом случае, в котором для упрощения теплофизические характеристики расплава и твердого полимера будем считать одинаковыми. Пусть скрытая теплота плавления равна л , а температура плавления Тп. Обозначим координату поверхности раздела между твердой и жидкой фазами через Х( t ). То гда одно из граничных услови й которое должно удовлетворяться на этой поверхности, запишется в виде: T s = T m = T n при X=X(t ) ( 2.2 4 ) Индекс s указывает, что соответствующая величина относится к твердой фазе (например, с s — плотность твердой фазы). Соответственно индекс m указывает, что величина относится к жидкой фазе. Второе граничное условие касается поглощения (или выделения) скрытой теплоты на поверхности раздела. Предположим, что в области x > x ( t ) находится жидкость при температуре Т т (х, t ), а в области x = x ( t ) — твердая фаза при температуре T s ( x t t ). Если поверхность раздела перемещается на расстояние dx , то в элементе объема вещества выделяется и должно быть отведено в результате теплопроводности количество тепла, в пересчете на единицу поверхности равное сdx . Математически это условие за пишется в виде: (2.25) Рассмотрим три случая: плавление, затвердевание и плавление с удалением расплава. 2.3.1. Плавление в области х > 0. Если в начальный момент область х > 0 занята твердым телом с постоянной температурой T s 0 и при t > 0 плоскость х = 0 поддерживается при постоянной темпера туре Т 2 > Т п , то положение плоскости плавления определится вы ражением: (2.26) Здесь - корень уравнения (2.27) где ; При этом распределение температур в т в ё рдой и жидкой фазах описывается выражением: (2.28) (2.29) 2.3.2. Затвердевание. Пусть в начальный момент времени область х > 0 представляет собой жидкость, а область х <С 0 — твердое тело. Иначе говоря, в начальный момент поверхность раздела сов падает с началом координат. Допустим, что значения термических коэффициентов только что затвердевшего расплава отличаются от значений термических коэффициентов твердой фазы в области х < 0. Присвоим термиче ским коэффициентам этой области индекс s 0 . Поступающий расплав имеет те мпературу Т 2 . Координата по верхности раздела фаз определится соотношением: (2.30) Здесь о — корень уравнения (2.31) После определения о , которое может быть выполнено любым численным методом (например, методом итерации), можно опре делить температурные поля во всех трех областях (начальная твердая фаза, затвердевшее вещество и расплав): (2.34) (2.35) (2.35) 2.3.3 Плавление с непрерывным удалением расплава. Пусть твердое тело нагревается благодаря поступающему извне к его поверхно сти постоянному тепловому потоку q . При этом весь расплав не прерывно удаляется. Примем плоскость, на которой происходит плавление, за плоскость с координатой х = 0 и будем считать, что твердое тело в области х > 0 движется относительно этой плос кости со скоростью х . Следовательно, массовый расход расплава, Q m , отнесенный к единичной ширине, равен : (2.36) В установившемся режиме температура в области х > 0 опи сывается выражением: (2.37) Из дифференциального уравнения теплопроводности следует, что тепловой поток в стационарном режиме равен нулю. Следовательно, количество тепла, подведенного извне в единицу времени, должно быть равно количеству тепла, отводимого в еди ницу времени с расплавом: (2.38) Определив х из соотношения ( 2.3 8), можно рассчитать рас пределение температур в твердом теле по формуле (2.36). Рассмотренные три случая наиболее типичны для процессов переработки полимеров, так как любой реальный процесс плавле ния можно свести к одному из них. 2.4.Т еплопередача в потоках расплава Передача тепла в движущейся жидкости происходит по механизму конвективного теплообмена, который осуществляется как за счет переноса тепла током жидкости, так и за счет теплопроводности самой жидкости. Аналитическое решение дифференциальных урав нений теплопроводности в случае конвективного теплообмена удается получить лишь при введении большого числа упрощений. Поэтому для практических целей используют результаты экспери ментальных исследований, представленные в виде зависимостей между соответствующими критериями подобия. Обычно при изу чении теплопередачи конвекцией принимаются следующие до пущения: 1) на границе с поверхностью нагрева (охлаждения) соблю даются условия прилипания; 2) физические параметры жидкости (теплоемкость, теплопроводность, плотность и вязкость) сохра няют неизменное значение для всего потока; 3) лучистый тепло обмен между поверхностью нагрева (охлаждения) и потоком жидкости происходит независимо от контактной теплоотдачи. В настоящее время наибольшее распространение получили экс* периментальные исследования процессов стационарного теплооб мена. Для описания процесса теплообмена обычно используется известное уравнение Ньютона: (2.39) где а — коэффициент теплоотдачи, определяющий количество тепла, подводимое (или отводимое) к жидкости в единицу времени через поверхность с единичной площадью; T w — температура стенки канала; Т ж — средняя температура жидкости. По своему физическому смыслу коэффициент теплоотдачи яв ляется условной величиной и характеризует отношение коэффициента теплопроводности жидкости к толщине д пристенного слоя, в котором происходит температурный скачок: (2.40) Использование методов теори и подобия позволяет свести реше ние проблемы теплообмена в потоке жидкости к экспериментальному определению вида функциональной зависимости: (2.41) Здесь — критерий Нуссельта, характеризующий интенсивность теплообмена; Р r = Ср м / — критерий Прандтля, характеризующий соотношение между количеством тепла, поглощаемого жидкостью за счет изменения энталь пии, и количеством тепла, отводимого за счет теплопроводности; Gr = g л P 2 l z Д T / м 2 — критерий Грасгофа, характеризующий интенсивность теплооб мена за счет свободной конвекции; Re = vlp /ц — число Рейнольдса, характери зующее отношение сил инерции к силам вязкого трения; Ре = vd / a — критерий Пекле; — критерий Гретца. Известные в настоящее время результаты экспериментального исследования теплообмена в расплавах полимеров относятся пре имущественно к течению в каналах круглого сечения. Общая фор мула имеет вид: (2.42) где индексы «Ж» и «ст» Означают, что соответствующие значения критерия от носятся к усредненным характеристикам жидкости или к характеристикам жид кости в пристенном слое. Значения показателей степени при критериях в уравнении (2.42) приведены ниже: Таблица (3.1) Значения показателей степени при критериях подобия . Полимер А X У Z Z 1 П П олиэтилен низкой плотности 16 0,33 0,33 0,15 0,33 П П олиэтилен низкой плотности 17 2,25 0,18 0,20 0,25 0 2.5. Лучистый т еплообмен Нагрев излучением применяется главным образом в операциях, предшествующих пневмо- и вакуум-формованию относительно тон ких листов термопластов. Лучистая энергия передается в виде электромагнитных волн, распространяющихся в пространстве до тех пор, пока на их пути не встретится какая-либо поглощающая среда: газ, жидкость или твердое тело. Излучаемая энергия пропорциональна четвертой степени абсолютной температуры изучающего тела. Так как обычно большая часть энергии излучения в применяемой на прак тике области температур приходится на инфракрасный спектр, нагрев излучением называют также инфракрасным нагревом. Гипотетическое тело, поглощающее все падающие на него лучи, называется абсолютно черным телом. Интенсивность лучеиспуска ния абсолютно черного тела Е b определяется законом Стефана — Больцмана: (2.43) Где а — постоянная Стефана Больцмана, равная 1,36 • 10 -1 2 кал/(см 2 • с • / K 4 ), или Реальные тела излучают меньше энергии. Их излучательная способность е оценивается по формуле: (2.44) где Е — интенсивность лучеиспускания реального тела. Обычно е зависит от температуры, увеличиваясь с ее ростом. Металлоиды и окислы металлов обладают высокой излучательной способностью ( е
1Архитектура и строительство
2Астрономия, авиация, космонавтика
 
3Безопасность жизнедеятельности
4Биология
 
5Военная кафедра, гражданская оборона
 
6География, экономическая география
7Геология и геодезия
8Государственное регулирование и налоги
 
9Естествознание
 
10Журналистика
 
11Законодательство и право
12Адвокатура
13Административное право
14Арбитражное процессуальное право
15Банковское право
16Государство и право
17Гражданское право и процесс
18Жилищное право
19Законодательство зарубежных стран
20Земельное право
21Конституционное право
22Конституционное право зарубежных стран
23Международное право
24Муниципальное право
25Налоговое право
26Римское право
27Семейное право
28Таможенное право
29Трудовое право
30Уголовное право и процесс
31Финансовое право
32Хозяйственное право
33Экологическое право
34Юриспруденция
 
35Иностранные языки
36Информатика, информационные технологии
37Базы данных
38Компьютерные сети
39Программирование
40Искусство и культура
41Краеведение
42Культурология
43Музыка
44История
45Биографии
46Историческая личность
47Литература
 
48Маркетинг и реклама
49Математика
50Медицина и здоровье
51Менеджмент
52Антикризисное управление
53Делопроизводство и документооборот
54Логистика
 
55Педагогика
56Политология
57Правоохранительные органы
58Криминалистика и криминология
59Прочее
60Психология
61Юридическая психология
 
62Радиоэлектроника
63Религия
 
64Сельское хозяйство и землепользование
65Социология
66Страхование
 
67Технологии
68Материаловедение
69Машиностроение
70Металлургия
71Транспорт
72Туризм
 
73Физика
74Физкультура и спорт
75Философия
 
76Химия
 
77Экология, охрана природы
78Экономика и финансы
79Анализ хозяйственной деятельности
80Банковское дело и кредитование
81Биржевое дело
82Бухгалтерский учет и аудит
83История экономических учений
84Международные отношения
85Предпринимательство, бизнес, микроэкономика
86Финансы
87Ценные бумаги и фондовый рынок
88Экономика предприятия
89Экономико-математическое моделирование
90Экономическая теория

 Анекдоты - это почти как рефераты, только короткие и смешные Следующий
- Петр Алексеевич в Турции военный переворот, мы за кого?
- А за кого Путин?
- Путин молчит.
- Тогда давайте у Обамы спросим - за кого он, за того и мы.
- Алло, Барак Хусейнович, в Турции военный переворот, вы за кого?
- А за кого Путин?
Anekdot.ru

Узнайте стоимость курсовой, диплома, реферата на заказ.

Банк рефератов - РефератБанк.ру
© РефератБанк, 2002 - 2016
Рейтинг@Mail.ru