RSS    

   Реферат: Расчет затвердевания плоской отливки

7)   Рассчитать прогоночные коэффициенты ai+1, bi+1 для Метв. и Меж., s=0 с учетом что Тnf=Тз.

Т1=р2Т2+g2

Тi=а2Т2+в2

Найти а2 и в2:

а2=1,                                                           (45)

                                                                            (46)

8)   Рассчитать температуру на оси симметрии

                                                                                   (47)

9)   Рассчитать температурное поле жидкого и твердого металла

                                                                              (48)

10)    Пересчитать значения ∆tк по итерационному процессу (24)

d – параметр итерации (d=0…1)

проверяем точность;

11)    Скорость охлаждения в каждом узле i рассчитать по формуле:

, оС/с                                                                         (50)

12)    Скорость затвердевания на каждом временном шаге:

, м/с                                                                                    (51)

13)    Средняя скорость охлаждения на оси отливки:

14)    Положение фронта затвердевания по отношению к поверхности отливки

, к – шаг по времени                                                (52)

15)    Полное время затвердевания

, к′ - последний шаг                                                   (53)

16)    Средняя скорость затвердевания отливки

                                                                                             (54)

Идентификаторы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Блок-схема

- [Вводим исходные данные

- [Вычисляем шаг по пространству

- [Вычисляем коэффициенты Аj, Сj для подстановки в (32), (33) и задаем температуру в первой точке

- [Температурное поле для первого шага по времени

- [Делаем шаг по времени

- [Вычисляем плотность теплового потока

- [Шаг по времени в нулевом приближении

- [Начальные прогоночные коэффициенты

- [Шаг по итерации

- [Вычисляем коэффициенты Bj для подстановки в (32), (33)

- [Вычисляем прогоночные коэффициенты по твердому металлу

- [Прогоночные коэффициенты для фронта

- [Вычисляем прогоночные коэффициенты по жидкому металлу

- [Температура на оси симметрии

- [Расчет температурного поля

- [Ищем максимальный температурный шаг

- [Уточняем Dt

- [Точность временного шага

- [Проверка точности

- [Расчет времени

- [Скорость охлаждения в каждом узле

- [Скорость затвердевания и положение фронта

- [Вывод результатов

- [Проверка достижения фронтом центра отливки

- [Расчет полного времени, ср. скорости затвердевания ср. скорости охлаждения на оси отливки

Вывод результатов

- [Конец.

Программа

CLEAR , , 2000

    DIM T(1000), T1(1000), AP(1000), BP(1000), Vox(1000), N$(50)

2 CLS

  N = 100: KV = 50: N9 = 5: L = .015

  TM = 293: TI = 1345: TS = 1312.5

  BM = 1300: a1 = .000036: a2 = .000021

  TA0 = .01: ETA = .01: E = .01

  l1 = 195: l2 = 101

  R0 = 8600: LS = 221000

  AF = 0: Pi = 3.14159265359#

3 PRINT "Число шагов N, штук"; N

  PRINT "Длина отливки L, м"; L

  PRINT "Температура формы Tf, К"; TM

  PRINT "Начальная температура сплава Tн, К"; TI

  PRINT "Температура затвердевания Tz, К"; TS

  PRINT "Bф "; BM

  PRINT "Первый шаг по времени, Tk0 "; TA0

  PRINT "Точность по времени, Еt "; ETA

  PRINT "Точность по температуре, ЕТ "; E

  PRINT "Температуропроводность Ме твердого, а1 "; a1

  PRINT "Температуропроводность Ме жидкого, а2 "; a2

  PRINT "LS= "; LS

  PRINT "Коэф. теплопроводности, l1 "; l1

  PRINT "Коэф. теплопроводности, l2"; l2

  PRINT "Плотность Ме твердого, р1 "; R0

  INPUT "Изменить данные <y/n>"; QV$

  IF QV$ = "Y" THEN GOSUB 222

48  N1 = N - 1

    DX = L / (N - 1)

    A = a1 / DX ^ 2

    B1 = 2 * A

    RL = R0 * LS * DX

    NF = 1

    B2 = l1 / DX

    KV1 = 1

    AL = a2 / DX ^ 2

    BL1 = 2 * AL

    BL2 = l2 / DX

    T(1) = TS

    T1(1) = TS

    FOR i = 2 TO N

    T(i) = TI

    T1(i) = TI

    NEXT i

    TA = TA0

    K = 1

    dta = .01

    GOTO 103

101 K = K + 1

    NF = NF + 1

    B3 = SQR(Pi * TA)

    q = BM * (T(1) - TM) / B3

    dta = RL / (AF + q)

    B5 = BM * TM / B3

    B3 = BM / B3

    B4 = B2 + B3

    AP(1) = B2 / B4

    BP(1) = B5 / B4

    T(NF) = TS

    NF1 = NF - 1

    NF2 = NF + 1

    K1 = 0

102 K1 = K1 + 1

    Et = 0

  

    B3 = SQR(Pi * (TA + dta))

    q = BM * (T(1) - TM) / B3

    B5 = BM * TM / B3

    B3 = BM / B3

    B4 = B2 + B3

    AP(1) = B2 / B4

    BP(1) = B5 / B4

Страницы: 1, 2, 3, 4


Новости


Быстрый поиск

Группа вКонтакте: новости

Пока нет

Новости в Twitter и Facebook

                   

Новости

© 2010.