Построение и анализ моделей процессов и систем в технологии деревообработки

Министерство  образования РФ

Санкт –  Петербургская государственная  лесотехническая академия

имени С. М. Кирова 
 
 
 

Кафедра технологии деревообрабатывающих

производств 
 
 
 
 
 
 

КУРСОВАЯ  РАБОТА 

по дисциплине: «Основы моделирования и оптимизации  процессов деревообработки»

на тему: «Построение и анализ моделей  процессов и систем в технологии деревообработки» 
 
 
 

          МТ.ДТ.43.301061.КР.9
 
 
 
 
 
 
 
 
 
 

Студент:                                                                                            С. В. Солнышков

Преподаватель:                                                                                 И. В. Коваленко 
 
 
 

Санкт –  Петербург 

2005

ОГЛАВЛЕНИЕ 
 

1. МОДЕЛИРОВАНИЕ  ТЕХНОЛОГИЧЕСКОГО ПРОЦЕССА…………….……3

1.1. Моделирование  и анализ процесса смешения  при подготовке

       рабочего   раствора лака………………………………………………………..3

2. ПОСТРОЕНИЕ  МОДЕЛИ ЗАДАЧИ ОПТИМАЛЬНОГО

       РАСКРОЯ ПЛИТНЫХ МАТЕРИАЛОВ……………………………………..10

 Вариант задания  № 1………………………………………………………………10

       Рассмотрение ситуации № 1…………………………………………………..16

       Рассмотрение ситуации № 2…………………………………………………..17

       Рассмотрение ситуации № 3…………………………………………………..18

       Рассмотрение ситуации № 4…………………………………………………..20

       Рассмотрение ситуации № 5…………………………………………………..22

       Рассмотрение ситуации № 6…………………………………………………..28

3. АНАЛИЗ И  СИНТЕЗ СТРУКТУРЫ ПРОИЗВОДСТВЕННОЙ 

       СИСТЕМЫ  «СКЛАД - СТАНОК»………………………………...…………31

3.1. Задача исследования  одноканальной системы с отказами………….………31

3.2. Построение  графа………………………………………………………….…..31

3.3. Построение математической модели…………………………………………34

3.4. Анализ функционирования  системы…………………………………………36

3.5. Выводы…………………………………………………………………………37

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ………………………………38 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

1. МОДЕЛИРОВАНИЕ  ТЕХНОЛОГИЧЕСКОГО ПРОЦЕССА 
 
 

     1.1. Моделирование и анализ процесса  смешения при подготовке рабочего   раствора лака 
 

     Для построения модели процесса смешения возьмём нитроцеллюлозный лак марки  НЦ – 218 с типовым составом, % 

     Плёнкообразующие                                        Летучая часть 

Нитроцеллюлоза  ……...….4 – 20                  Растворитель ………………..18 – 36

Смола (синтетическая) …..5 – 8                    Разбавители, разжижители….72 – 24

Пластификатор …………..1 – 12

      _________________________              ______________________________

             Итого                10 – 40                                    Итого                       90 – 60  
 
 

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

     В рассматриваемой модели смешения контролируемым параметром является концентрация растворителя в лаке. Рабочий раствор ЛКМ готовится в смесительной установке добавлением к основе лака (первый входной поток) чистого растворителя (второй входной поток). В результате перемешивания концентрация растворителя на выходе смесителя изменяется во времени, достигая требуемого значения.

     Перемешивание в смесителе двух входных потоков  с объёмными расходами  и концентрациями растворителя рассматривается как идеальное. Объём смесителя при этом постоянный; концентрация растворителя постоянна по всему объёму смесителя и равна концентрации на выходе ; расход на выходе смесителя составляет ; и задано отношение во входных потоках . 
 
 
 
 
 
 

       

                                   Рис.1.1 – Схема смесителя 
 

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

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

                                        

                      (1.1) 

     Причём 

                                          

,                    (1.2) 
 
 

     и

                                                          

,                                 (1.3) 

     где Т – постоянная времени переходного  процесса;

          - относительное изменение концентрации компонента;

           - текущее время протекания процесса;                                           

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

На основании  вышеприведённых формул произведём расчёт расхода смолы ( ), выразив его через формулу (1.2) 

                                                            

,                                       (1.4) 

                          

дм³/мин 

     Зная  расход смолы ( ) можно рассчитать расход растворителя ( ) по формуле 

                                                            ,                                             (1.5) 

     где G – расход на выходе смесителя. 

             

дм³/мин 
 

     Используя формулу (1.3) можно рассчитать конечную концентрацию смолы ( )

                                                   

,                                              (1.6) 

     

 % 

     Постоянная  времени переходного процесса ( Т ) определяется по формуле 

                                                               

,                                             (1.7) 

     Чтобы найти объём смесителя необходимо выполнение следующего условия 

                                                            

                                   (1.8) 

     где - допустимый предел концентрации растворителя на выходе из

                  смесителя через конечное время  ;

           - относительное изменение концентрации растворителя в первом 

                  входном потоке;

            m – отношение расходов растворителя во входных потоках. 

     

 

                     

 

     Так как условие (1.8) не выполняется, то минимальный  объём смесителя необходимо рассчитать по формуле 

                                                        

,                               (1.9) 

                         

дм³ 

     Определим постоянную времени переходного  процесса ( Т ) по формуле (1.7) 

                           

мин 
 

     Для того чтобы построить график зависимости изменения концентрации растворителя от времени необходимо подставлять в формулу (1.1) значения времени ( ) от 0 до 100 мин. Значение концентрации для начальных условий мы определили, теперь рассчитаем остальные точки графика. 

     При

               

% 

     При

              

% 

     Аналогично  рассчитываются значения для всех остальных  точек графика и результаты расчетов сводятся в табл. 1.1. 

Таблица 1.1 –  Зависимость концентрации растворителя от времени

Время
, мин
0 10 20 30 40 50 60 70 80 90 100
Концентрация

, %

18 26,54 33,63 39,5 44,37 48,42 52,33 54,56 56,88 58,78 60,37
 

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

     

     Концентрацию  растворителя можно рассчитать, используя формулу (1.2) 

      ,                                         (1.10)

        

%

     Подставив значения, получаем 

         

% 
 
 

     На  основании табл.1.1. построим график 
 

        

     Рис. 1.1. – График зависимости концентрации растворителя от времени 
 
 
 
 

3. АНАЛИЗ  И СИНТЕЗ СТРУКТУРЫ ПРОИЗВОДСТВЕННОЙ  СИСТЕМЫ  «СКЛАД - СТАНОК». 

3.1. Задача исследования  одноканальной системы  с отказами 

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

      Заявки на обслуживание (детали) поступают с интервалом tп = 50 мин, время обработки заявки составляет tоб = 40 мин, количество подстопных мест на складе m = 6. В курсовой работе требуется построить ориентированный размеченный граф состояний системы с заданными числовыми характеристиками процесса, построить математическую модель системы, рассчитать операционные характеристики системы и провести анализ ее функционирования, выработать рекомендации, направленные на повышение эффективности рассматриваемой системы. 

3.2. Построение графа 

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

      Интенсивность потока поступления  заявок  мин-1, определим по формуле 

                                                                                                                 (3.1) 

                                                                мин-1. 

      Интенсивность потока обслуживания  заявок  мин-1, узнаем по формуле 

                                                                                                               (3.2) 

                                                                 мин-1. 

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

                                                                                                             (3.3) 

                                                               

                 

       Вероятность состояния системы,  при котором в ней находится  заявок , определим по формуле 

                                                                                               (3.4) 

где k = . 

 
 

 

 

 

 

 

 

 

 

 

 

 

 
 
 

      Среднее число заявок, находящихся  в системе найдем по формуле 

                                      

,                    (3.5) 
 

               

 
 
 
 
 
 

      Остальные расчеты сведем в таблицу 3.1 

Таблица 3.1 

m 1 2 3 4 5 6
NC 0,85 1,23 1,56 1,87 2,14 2,38

   

     Абсолютную пропускную способность  системы найдем по формуле  5.7 

       , (3.6) 

А = 1,2·(1 - 0,05) = 1,14 

Среднее число  заявок находящихся в очереди  найдем по формуле 5.8 

       , (3.7) 

 

Среднее время  ожидания заявок в очереди найдем по формуле 5.9 

         (3.8) 

мин 

 Среднее  время пребывания заявок в  системе определим по формуле  5.10 

                                                                                                              (3.9) 

                

мин. 
 

      Ориентированный граф состояний  системы представлен на рис. 3.1, где возможные состояния системы обозначены: S0 – система свободна; S1 – станок занят обработкой одной стопы и очереди нет (склад свободен); S2 – станок работает, одна стопа на складе; Sк – станок занят, ( ) стоп в очереди; Sm – станок занят ( ) стоп на складе; Sm+1-система занята и отказывает в обслуживании, т. е. станок работает,m стоп на складе, а в системе (m+1) cтоп.  
 
 
 

 λ λ λ

                            

                                                                                                        

                           

  

      

                           Рис. 3.1. Ориентированный граф состояний  системы 
 

                               3.3. Построение математической модели 

      Математическая модель системы  будет описана формулами (3.4); (3.5), для построения графиков зависимости и выполним расчеты зависимости для и при m = 1; 2; 3; 4; 5; 6; и представим их виде табл. 3.1, пример расчетов см. в пункте 3.2. Графики и представлены на рис. 3.2 и рис. 3.3, соответственно. 
 
 

Таблица 3.2 

Зависимость Nи Рm+1 от m 

                                 

 
 
 
 
 
 
 
 
1 0,25 0,85
 

2 0,173 1,23
3 0,122 1,56
4 0,0889 1,87
5 0,0663 2,14
6 0,0504 2,38
7 0,0388 2,605
8 0,0301 2,797
9 0,0235 2,966
 

Рис. 3.2. График зависимости  и  
 

3.4. Анализ функционирования  системы 

      Зная операционные характеристики многоканальной системы (рассчитанные в пунктах 3.2 и 3.3), можно определить ее структуру. В качестве критерия оптимизации принимаются удельные приведенные затраты, характеризующие затраты в системе на одно подстопное место: 

                                                                  (3.10) 

    где   - экономические потери, обусловленные вынужденным ожиданием  

                   обработки одной стопы;

         с2 - затраты на организацию подстопного места на складе;

     с3 - потери, связанные с невозможностью разместить на складе еще одну  

         стопу. 
     

 
 
 
 
 
 

Остальные расчеты  сведем в таблицу 3.1 

Таблица 3.3 

m 1 2 3 4 5 6
W 94 76,12 75,48 78,74 83,71 89,17
 
 

                                    Рис. 3.4. График зависимости

  
 

3.5. Выводы 

      На основании данных, полученных  в пунктах 3.2; 3.3; 3.4 можно утверждать, что при данных условиях ( = 50 мин; = 40 мин) оптимальное число подстопных мест m = 2. 
 
 
 
 
 
 
 
 

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ 

1. Паянский –  Гвоздев В. М., «Моделирование  и оптимизация процессов деревообработки»: Методические указания по выполнению  контрольной работы, С. – Петербург, 1992 г.

2. Паянский –  Гвоздев В. М., «Моделирование  и оптимизация процессов деревообработки»: Методические указания по выполнению  курсовой работы, С. – Петербург, 1991 г.

3. Егоров В.  А., Глушковский А. А. «Основы  и оптимизация процессов деревообработки». Методы моделирования: Учебное пособие. Л.: ЛТА, 1988 г. 
 
 
 
 

Построение и анализ моделей процессов и систем в технологии деревообработки