Перейти к публикации

Хэлп новичку - стационарное поле температур системе тел с газовыми прослойками


Рекомендованные сообщения

Коллеги, приветствую всех!

 

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

Задача моделируется в 3d (хотя и можно было пользуясь осесимметричностью решать 2d)

 

По классике - построил солид тела, разбил их mesh'ером с использованием solid278... Задал ГУ, казалось бы все просто ... И не считает! Долго бился, но так и не нашел причину сбоя.

 

Ответ мне заранее известен - сделал расчет по одномерным моделям (расчетные формулы из курса теплопроводности).

 

Помогите, пожалуйста, найти ошибку.

 

 

И да... Пишу на форум впервые, могу тормозить :bash:

Ссылка на сообщение
Поделиться на других сайтах


UnPinned posts
2 часа назад, bubblyk сказал:

Помогите, пожалуйста, найти ошибку.

в коде? видимо, нужен код. или хотя бы текст ошибок из аутпута. а то совсем цепляться незашто... 

Ссылка на сообщение
Поделиться на других сайтах

Я уже торможу. Не понимаю, где здесь инструмент для вставки файла - кода.

Просто вбивать текстом как сообщение - слишком длинно.

Код готов передать, подскажите как :smile:

Ссылка на сообщение
Поделиться на других сайтах
21 minutes ago, bubblyk said:

Не понимаю, где здесь инструмент для вставки файла - кода.

Просто вбивать текстом как сообщение - слишком длинно.

image.png

Ссылка на сообщение
Поделиться на других сайтах
13 минут назад, piden сказал:

image.png

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

Не допёр...

Ссылка на сообщение
Поделиться на других сайтах
1 minute ago, bubblyk said:

как сюда приложить текстовый файл с набором команд.

Да нажмите кнопку спойлера и туда весь командный файл вставьте. В чем сложность?

 

После 5 или 10 сообщений появится кнопка для приложения файлов. Форум просто хочет убедиться в твердости ваших намерений)

 

17 hours ago, bubblyk said:

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

Т.е. взаимодействие происходит излучением и теплопередачей через движущийся теплоноситель?

 

Ссылка на сообщение
Поделиться на других сайтах

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

 

Скрытый текст

/PREP7
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ЗАДАНИЕ ТИПОВ КОНЕЧНЫХ ЭЛЕМЕНТОВ !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ET,    1,    SOLID278
ET,    2,    PLANE55
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ЗАДАНИЕ ИСХОДНЫХ ПАРАМЕТРОВ !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    *SET,    D_CLAD_IN,    7.73E-3
    *SET,    D_CLAD_OUT,    9.13E-3
    *SET,    H_CLAD,        477.90E-3
    *SET,    D_CORE,    7.60E-3
    *SET,    H_CORE,    100.0E-3
    *SET,    D_Zr_SIM_OUT,    7.60E-3
    *SET,    D_Zr_SIM_IN,    3.00E-3
    *SET,    H_Zr_SIM1,    26.0E-3
    *SET,    H_Zr_SIM2,    52.0E-3
    *SET,    H_Zr_SIM3,    13.0E-3
    *SET,    D_Hf_ABSORB,    7.60E-3
    *SET,    H_Hf_ABSORB,    6.45E-3
    *SET,    D_W_FOIL,    7.60E-3
    *SET,    H_W_FOIL,    50.0E-6
    *SET,    D_END_PLUG_BOT,    7.73E-3
    *SET,    H_END_PLUG_BOT,    7.60E-3
    *SET,    D_END_PLUG_TOP,    7.73E-3
    *SET,    H_END_PLUG_TOP,    4.80E-3
    *SET,    L_DIM_MAX,    477.90E-3
    *SET,    L_FREE,        53.5E-3
    *SET,    N_TIME1,    3
    *SET,    N_TIME2,    2
    *SET,    COORD_Z_MAX,    0.0
    *SET,    N_PHI,    5
    *SET,    N_Z,    8
    *SET,    Q_LIN,    410E2
!!! ПАРАМЕТРЫ ПОСТРОЕНИЯ СЕТКИ
MSHAPE,    2,    3D
MSHKEY,        2

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!    ЗАДАНИЕ СВОЙСТВ МАТЕРИАЛОВ        !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

TOFFST,0    !!!  ЗАДАНИЕ ЗАВИСИМОСТЕЙ СВОЙСТВ ОТ ТЕМПЕРАТУРЫ (T, K)
MPTEMP,1,    0.0E3,    0.3E3,    0.4E3,    0.5E3,    0.6E3,    0.7E3
MPTEMP,7,    0.8E3,    0.9E3,    1.0E3,    1.1E3,    1.2E3,    1.3E3
MPTEMP,13,    1.4E3,    1.5E3,    1.6E3,    1.7E3,    1.8E3,    1.9E3
MPTEMP,19,    2.0E3

MPDATA,DENS,1,1,    1.2785E+04,    1.2784E+04,    1.2727E+04,    1.2667E+04,    1.2603E+04,    1.2536E+04
MPDATA,DENS,1,7,    1.2466E+04,    1.2392E+04,    1.2244E+04,    1.2108E+04,    1.2035E+04,    1.1977E+04
MPDATA,DENS,1,13,    1.1932E+04,    1.1912E+04,    1.1903E+04,    1.1893E+04,    1.1884E+04,    1.1874E+04
MPDATA,DENS,1,19,    1.1864E+04

MPDATA,C,1,1,    1.57E+02,    1.57E+02,    1.66E+02,    1.77E+02,    1.89E+02,    2.03E+02
MPDATA,C,1,7,    2.18E+02,    2.34E+02,    2.25E+02,    2.15E+02,    1.99E+02,    2.01E+02
MPDATA,C,1,13,    2.02E+02,    2.04E+02,    2.06E+02,    2.08E+02,    2.11E+02,    2.13E+02
MPDATA,C,1,19,    2.16E+02

MPDATA,KXX,1,1,        1.950E1,    1.973E1,    2.318E1,    2.418E1,    2.502E1,    2.618E1
MPDATA,KXX,1,7,        2.754E1,    2.896E1,    3.045E1,    3.213E1,    3.407E1,    3.619E1
MPDATA,KXX,1,13,    3.828E1,    4.004E1,    4.131E1,    4.217E1,    4.289E1,    4.362E1
MPDATA,KXX,1,19,    4.341E1

! МАТЕРИАЛ №2

MPDATA,DENS,2,1,    6.505E+03,    6.505E+03,    6.486E+03,    6.467E+03,    6.447E+03,    6.428E+03
MPDATA,DENS,2,7,    6.409E+03,    6.390E+03,    6.371E+03,    6.352E+03,    6.332E+03,    6.313E+03
MPDATA,DENS,2,13,    6.294E+03,    6.275E+03,    6.256E+03,    6.236E+03,    6.217E+03,    6.198E+03
MPDATA,DENS,2,19,    6.179E+03

MPDATA,C,2,1,    0.280E3,    0.280E3,    0.294E3,    0.308E3,    0.321E3,    0.333E3
MPDATA,C,2,7,    0.345E3,    0.356E3,    0.367E3,    0.378E3,    0.314E3,    0.319E3
MPDATA,C,2,13,    0.326E3,    0.333E3,    0.341E3,    0.350E3,    0.359E3,    0.370E3,
MPDATA,C,2,19,    0.381E3

MPDATA,KXX,2,1,        1.670E+01,    1.691E+01,    2.186E+01,    2.149E+01,    2.100E+01,    2.142E+01
MPDATA,KXX,2,7,        2.216E+01,    2.264E+01,    2.293E+01,    2.342E+01,    2.673E+01,    2.858E+01
MPDATA,KXX,2,13,    3.036E+01,    3.205E+01,    3.361E+01,    3.501E+01,    3.622E+01,    3.721E+01
MPDATA,KXX,2,19,    3.795E+01

! МАТЕРИАЛ №3

MPDATA,DENS,3,1,    13.31E3,    13.31E3,    13.28E3,    13.26E3,    13.23E3,    13.21E3
MPDATA,DENS,3,7,    13.18E3,    13.15E3,    13.13E3,    13.10E3,    13.07E3,    13.04E3
MPDATA,DENS,3,13,    13.01E3,    12.98E3,    12.95E3,    12.92E3,    12.88E3,    12.84E3
MPDATA,DENS,3,19,    12.810E3

MPDATA,C,3,1,    143.5,    143.6,    148.4,    153.1,    157.6,    162.1
MPDATA,C,3,7,    166.6,    171.2,    175.7,    180.1,    184.6,    189.2
MPDATA,C,3,13,    193.6,    198.1,    202.6,    207.1,    211.6,    216.1
MPDATA,C,3,19,    220.5

MPDATA,KXX,3,1,        25.2,    25.0,    23.4,    22.2,    21.4,    20.8
MPDATA,KXX,3,7,        20.5,    20.3,    20.4,    20.6,    21.1,    21.7
MPDATA,KXX,3,13,        22.5,    23.5,    24.7,    26.1,    27.7,    29.4
MPDATA,KXX,3,19,        31.3

! МАТЕРИАЛ №4

MPDATA,DENS,4,1,    19.27E3,    19.26E3,    19.23E3,    19.21E3,    18.18E3,    19.15E3
MPDATA,DENS,4,7,    19.12E3,    19.10E3,    19.07E3,    19.04E3,    19.01E3,    18.98E3
MPDATA,DENS,4,13,    18.95E3,    18.92E3,    18.89E3,    18.86E3,    18.82E3,    18.77E3
MPDATA,DENS,4,19,    18.72E3

MPDATA,C,4,1,        0.131E3,    0.132E3,    0.136E3,    0.138E3,    0.140E3,    0.142E3
MPDATA,C,4,7,        0.144E3,    0.146E3,    0.148E3,    0.150E3,    0.152E3,    0.155E3
MPDATA,C,4,13,        0.160E3,    0.162E3,    0.164E3,    0.166E3,    0.170E3,    0.175E3
MPDATA,C,4,19,        0.177E3

MPDATA,KXX,4,1,        0.1650E3,    0.1628E3,    0.1559E3,    0.1459E3,    0.1371E3,    0.1302E3
MPDATA,KXX,4,7,        0.1242E3,    0.1202E3,    0.1166E3,    0.1151E3,    0.1135E3,    0.1122E3
MPDATA,KXX,4,13,    0.1112E3,    0.1105E3,    0.1101E3,    0.1095E3,    0.1090E3,    0.1086E3
MPDATA,KXX,4,19,    0.1082E3

! МАТЕРИАЛ №5

MPDATA,DENS,5,1,    161.0E-3,    160.4E-3,    120.3E-3,    96.3E-3,    80.2E-3,    68.8E-3
MPDATA,DENS,5,7,    60.2E-3,    53.5E-3,    48.1E-3,    43.8E-3,    40.1E-3,    37.0E-3
MPDATA,DENS,5,13,    34.4E-3,    32.1E-3,    30.1E-3,    28.3E-3,    26.7E-3,    25.3E-3
MPDATA,DENS,5,19,    24.1E-3

MPDATA,C,5,1,        0.131E3,    0.132E3,    0.136E3,    0.138E3,    0.140E3,    0.142E3
MPDATA,C,5,7,        0.144E3,    0.146E3,    0.148E3,    0.150E3,    0.152E3,    0.155E3
MPDATA,C,5,13,        0.160E3,    0.162E3,    0.164E3,    0.166E3,    0.170E3,    0.175E3
MPDATA,C,5,19,        0.177E3

MPDATA,KXX,5,1,        5.197E3,    5.197E3,    5.197E3,    5.197E3,    5.197E3,    5.197E3
MPDATA,KXX,5,7,        5.197E3,    5.197E3,    5.197E3,    5.197E3,    5.197E3,    5.197E3
MPDATA,KXX,5,13,    5.197E3,    5.197E3,    5.197E3,    5.197E3,    5.197E3,    5.197E3
MPDATA,KXX,5,19,    5.197E3


!!!!!!!!!!!!!!!!!!!!!!!! ПЕРЕСТРОЙКА ГРАФИЧЕСКОГО ОТОБРАЖЕНИЯ ГЕОМЕТРИИ МОДЕЛИ В ИЗОМЕТРИИ !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ВКЛЮЧЕНИЕ ОТОБРАЖЕНИЯ НОМЕРОВ КЛЮЧЕВЫХ ТОЧЕК, ЛИНИЙ, ПЛОЩАДЕЙ И ОБЪЕМОВ !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/PNUM,KP,1
/PNUM,LINE,1
/PNUM,AREA,1
/PNUM,VOLU,1
!SHPP,    OFF,    ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ЗАДАНИЕ КЛЮЧЕВЫХ ТОЧЕК НИЖНЕЙ ЗАГЛУШКИ (ОБЪЕМ 1-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
K,    1,    0.0,    0.0,    0.0
K,    2,    D_CLAD_IN/2,    0.0,    0.0
K,    3,    0.0,    D_CLAD_IN/2,    0.0
K,    4,    -D_CLAD_IN/2,    0.0,    0.0
K,    5,    0.0,    -D_CLAD_IN/2,    0.0
K,    6,    0.0,    0.0,    H_END_PLUG_BOT    
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ НИЖНЕЙ ЗАГЛУШКИ (ОБЪЕМ 1-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!
jmin = 2
jmax = 5
SLAB = 0
SLAB1 = 0
*DO,jj,jmin,jmax-1
    LARC,    jj,    jj+1,    1,    D_CLAD_IN/2
        *IF,    jj,    EQ,    jmax-1,    THEN
            LARC,    jmax,    jmin,    1,    D_CLAD_IN/2
            SLAB = SLAB +1
            SLAB1 = SLAB1 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB1 = SLAB1 +1
*ENDDO
LESIZE,ALL,,,    N_PHI,    1.0,    1,    ,    ,    0
!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ НИЖНЕЙ ЗАГЛУШКИ (ОБЪЕМ 1-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
AL,    1,2,3,4
*GET, A_MAX,AREA,0,NUM,MAX
AMESH,A_MAX
!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМНОЙ МОДЕЛИ НИЖНЕЙ ЗАГЛУШКИ (ОБЪЕМ 1-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
VEXT,    ALL,,,    0.0, 0.0, H_END_PLUG_BOT,    0.0, 0.0, 0.0
*GET, VOL1,VOLU,0,NUM,MAX
VOLD = 0
VNEW = VOL1
CM,ZAGLUSHKA_NIZ,VOLU
VATT,    2,,1,,
*GET, L_MAX, LINE, 0, NUM, MAX
LSEL,S,LOC,Z,    0.0,    H_END_PLUG_BOT
LSEL,U,LOC,Z,    0.0,    0.0
LSEL,U,LOC,Z,    H_END_PLUG_BOT,    H_END_PLUG_BOT
LESIZE,ALL,,,    N_Z,    1.0,    1    ,    ,    0
LSEL,ALL
VMESH, VOL1
COORD_Z_MAX = COORD_Z_MAX+H_END_PLUG_BOT   !!! ПОСТРОЕНИЕ ЛСК НА ТОРЦЕВОЙ ЧАСТИ МОДЕЛИ, НАИБОЛЕЕ УДАЛЕННОЙ ОТ ГСК

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ МАЛОГО ИМИТАТОРА (ФАКТИЧЕСКИЙ ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ МАЛОГО ИМИТАТОРА (ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        jk*H_Zr_SIM2+COORD_Z_MAX
K,    K_MAX+2,    D_Zr_SIM_OUT/2,    0.0,    jk*H_Zr_SIM2+COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_Zr_SIM_OUT/2,    jk*H_Zr_SIM2+COORD_Z_MAX
K,    K_MAX+4,    -D_Zr_SIM_OUT/2,    0.0,    jk*H_Zr_SIM2+COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_Zr_SIM_OUT/2,    jk*H_Zr_SIM2+COORD_Z_MAX
KSEL,S,LOC,Z,      COORD_Z_MAX,  COORD_Z_MAX
KSEL,S,KP,,    K_MAX+1,    K_MAX+5

SLAB2 = 0
jmin = K_MAX+1
jmax = K_MAX+4 + 1
*DO,jj,jmin+1,jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_Zr_SIM_OUT/2
        *IF,    jj,    EQ,    jmax-1,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_Zr_SIM_OUT/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO

LNUM_OLD = L_MAX        ! ПЕРЕПРИСВОЕНИЕ МАКСИМАЛЬНОГО НОМЕРА ЛИНИИ НА ПРЕДЫДУЩЕМ ШАГЕ
*GET, L_MAX, LINE, 0, NUM, MAX    ! ЗАПРОС МАКСИМАЛЬНОГО НОМЕРА ЛИНИИ
L_NUM = L_MAX            ! ПРИСВОЕНИЕ МАКСИМАЛЬНОГО НОМЕРА ЛИНИИ НА ТЕКУЩЕМ ШАГЕ

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ МАЛОГО ЦИРКОНИЕВОГО ИМИТАТОРА (ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!LSEL,S,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX
!LSEL,U,LINE,,1,SLAB1,1,    
AL,    ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМНОЙ ЗАГОТОВКИ МАЛОГО ИМИТАТОРА (ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX        ! ЗАПРОС МАКСИМАЛЬНОГО НОМЕРА ПОВЕРХНОСТИ
VOLD = VNEW
VEXT, A_MAX,,,    0.0, 0.0, H_Zr_SIM1,    0.0, 0.0, 0.0    ! ВЫТЯЖКА МАКСИМАЛЬНОГО НОМЕРА ПОВЕРХНОСТИ НА ДЛИНУ H_Zr_SIM1
*GET, VOL1,VOLU,0,NUM,MAX                
VNEW = VOL1                    ! СОХРАНЕНИЕ МАКСИМАЛЬНОГО НОМЕРА ОБЪЕМНОГО ТЕЛА В ПЕРЕМЕННОЙ VNEW
VSEL,U,VOLU,,1,VOLD,1,
CM,Zr_SIMULATOR_NIZ,VOLU            ! ПРИСВОЕНИЕ ТЕКУЩЕМУ ОБЪЕМНОМУ ТЕЛУ ИМЕНИ КОМПОНЕНТА  Zr_SIMULATOR_NIZ
VATT,    2,,1,,
VSEL,ALL                
*GET, L_MAX, LINE, 0, NUM, MAX            ! ЗАПРОС МАКСИМАЛЬНОГО НОМЕРА ЛИНИИ

LSEL,S,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX+H_Zr_SIM1
LSEL,U,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX
LSEL,U,LOC,Z,    COORD_Z_MAX+H_Zr_SIM1,    COORD_Z_MAX+H_Zr_SIM1
LESIZE,ALL,,,    N_Z,    1.0,    1    ,    ,    0
LSEL,ALL

VMESH, VOL1                    ! ЗАДАНИЕ КОМАНДЫ НА СЕТОЧНОЕ РАЗБИЕНИЕ ОБЪЕМНОГО ТЕЛА VOL1

COORD_Z_MAX = COORD_Z_MAX + H_Zr_SIM1        !!! ПОСТРОЕНИЕ ЛСК НА ТОРЦЕВОЙ ЧАСТИ МОДЕЛИ, НАИБОЛЕЕ УДАЛЕННОЙ ОТ ГСК

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ЦИКЛИЧНО ПОВТОРЯЮЩИХСЯ  ИМИТАТОРОВ _N_TIME1__РАЗ ОТ НИЖНЕЙ ЗАГЛУШКИ (ФАКТИЧЕСКИЕ ОБЪЕМЫ 3-5) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
VSEL,NONE
CM,Zr_SIMULATOR_BOT_SUM,VOLU
CM,GAS_PLENUM_BOT_SUM,VOLU
VSEL,ALL

*DO,jk,0,N_TIME1-1        ! ЗАДАНИЕ ЦИКЛА ПО ЧИСЛУ ИМИТАТОРОВ

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ОПОРНЫХ ТОЧЕК ПОПЕРЕЧНОГО СЕЧЕНИЯ ЗАГОТОВКИ  ИМИТАТОРА (ОБЪЕМ 5-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_Zr_SIM_OUT/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_Zr_SIM_OUT/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_Zr_SIM_OUT/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_Zr_SIM_OUT/2,    COORD_Z_MAX
KSEL,S,LOC,Z,      COORD_Z_MAX,  COORD_Z_MAX
KSEL,S,KP,,    K_MAX+1,    K_MAX+5
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ ЗАГОТОВКИ  ИМИТАТОРА  (ОБЪЕМ 5-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj, jmin+1, jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_Zr_SIM_OUT/2
        *IF,    jj+1,    EQ,    jmax,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_Zr_SIM_OUT/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
LNUM_OLD = L_MAX        ! ПЕРЕПРИСВОЕНИЕ МАКСИМАЛЬНОГО НОМЕРА ЛИНИИ НА ПРЕДЫДУЩЕМ ШАГЕ
*GET, L_MAX, LINE, 0, NUM, MAX    ! ЗАПРОС МАКСИМАЛЬНОГО НОМЕРА ЛИНИИ
L_NUM = L_MAX            ! ПРИСВОЕНИЕ МАКСИМАЛЬНОГО НОМЕРА ЛИНИИ НА ТЕКУЩЕМ ШАГЕ

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!  ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ ЗАГОТОВКИ  ИМИТАТОРА (ОБЪЕМ 5-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Z,    COORD_Z_MAX,  COORD_Z_MAX
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1
AL,    ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА ЗАГОТОВКИ  ИМИТАТОРА  (ОБЪЕМ 5-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
VOLD = VNEW
VEXT, A_MAX,,,    0.0, 0.0, H_Zr_SIM2,        0.0, 0.0, 0.0
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
VSEL,U,VOLU,,1,VOLD,1,
*GET, L_MAX, LINE, 0, NUM, MAX
LNUM_OLD = L_MAX
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ОПОРНЫХ ТОЧЕК ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В ИМИТАТОРЕ  (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        H_Zr_SIM2+COORD_Z_MAX
K,    K_MAX+2,    D_Zr_SIM_IN/2,    0.0,    H_Zr_SIM2+COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_Zr_SIM_IN/2,    H_Zr_SIM2+COORD_Z_MAX
K,    K_MAX+4,    -D_Zr_SIM_IN/2,    0.0,    H_Zr_SIM2+COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_Zr_SIM_IN/2,    H_Zr_SIM2+COORD_Z_MAX
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В ИМИТАТОРЕ  (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj, jmin+1, jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_Zr_SIM_IN/2
        *IF,    jj+1,    EQ,    jmax,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_Zr_SIM_IN/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO

LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В ИМИТАТОРЕ  (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Z,    H_Zr_SIM2+COORD_Z_MAX,  H_Zr_SIM2+COORD_Z_MAX
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1
AL,    ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА ОТВЕРСТИЯ В ИМИТАТОРЕ  (ОБЪЕМ 3-ИЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
WPOFFS,    0.0,    0.0,    H_Zr_SIM2+COORD_Z_MAX
CSYS,4
VOLD = VNEW
VEXT, A_MAX,,,    0.0, 0.0, -H_Zr_SIM2,        0.0, 0.0, 0.0
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
VSEL,U,VOLU,,1,VOLD,1,
CM,GAS_PLENUM_BOT,VOLU
VATT,    5,,1,,
CMSEL,A,GAS_PLENUM_BOT_SUM,VOLU
CM,GAS_PLENUM_BOT_SUM,VOLU
VSEL,ALL
*GET, L_MAX, LINE, 0, NUM, MAX
LNUM_OLD = L_MAX
CSYS,0  
WPAVE,0,0,0
CSYS,4
            !!!COORD_Z_MAX = (jk+1)*H_Zr_SIM2+COORD_Z_MAX - H_Zr_SIM2

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ВЫЧИТАНИЕ ИЗ ЗАГОТОВКИ ОТВЕРСТИЯ (ОБЪЕМ 7-ОЙ) - СОЗДАНИЕ ОБЪЕМА  ИМИТАТОРА !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
VSBV,    VOLD,    VNEW,    ,    KEEP,    KEEP
LSEL,S,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX + H_Zr_SIM2
LSEL,U,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX
LSEL,U,LOC,Z,    COORD_Z_MAX + H_Zr_SIM2,    COORD_Z_MAX + H_Zr_SIM2
LESIZE,ALL,,,    N_Z,    1.0,    1    ,    ,    0
LSEL,ALL
VSWEEP,    VNEW,    A_MAX,    A_MAX+1,    0
NUMCMP,AREA
VOLD = VNEW
*GET, VOL1,VOLU,0,NUM,MAX
*GET, A_MAX,AREA,0,NUM,MAX
VNEW = VOL1
VSEL,U,VOLU,,1,VOLD,1,
CM,Zr_SIMULATOR_BOT,VOLU
VATT,    2,,1,,
CMSEL,A,Zr_SIMULATOR_BOT_SUM,VOLU
CM,Zr_SIMULATOR_BOT_SUM,VOLU
VSEL,ALL
WPOFFS,    0.0,    0.0,    H_Zr_SIM2+COORD_Z_MAX
CSYS,0
VSWEEP,    VNEW,    A_MAX-1,    A_MAX,    0
COORD_Z_MAX = COORD_Z_MAX + H_Zr_SIM2
CSYS,0  
WPOFFS,    0.0,    0.0,    -COORD_Z_MAX
CSYS,4
*ENDDO                    !!! ОКОНЧАНИЕ ЦИКЛА ПО КОЛИЧЕСТВУ ПОВТОРЯЮЩИХСЯ  ИМИТАТОРОВ

                        !!! COORD_Z_MAX = COORD_Z_MAX + N_TIME1*H_Zr_SIM2    !!! ПОСТРОЕНИЕ ЛСК НА ТОРЦЕВОЙ ЧАСТИ МОДЕЛИ, НАИБОЛЕЕ УДАЛЕННОЙ ОТ ГСК


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА Hf ПОГЛОТИТЕЛЬ СНИЗУ (ФАКТИЧЕСКИЙ ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_Hf_ABSORB/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_Hf_ABSORB/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_Hf_ABSORB/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_Hf_ABSORB/2,    COORD_Z_MAX
KSEL,S,LOC,Z,      COORD_Z_MAX,  COORD_Z_MAX
KSEL,S,KP,,    K_MAX+1,    K_MAX+5
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ Hf СНИЗУ (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj, jmin+1, jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_Hf_ABSORB/2
        *IF,    jj+1,    EQ,    jmax,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_Hf_ABSORB/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO

LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ Hf СНИЗУ (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Z,    COORD_Z_MAX,  COORD_Z_MAX
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1
AL,    ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА Hf СНИЗУ  (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
VOLD = VNEW
VSEL,NONE
VEXT, A_MAX,,,    0.0, 0.0, H_Hf_ABSORB,        0.0, 0.0, 0.0
VATT,    3,,1,,
LSEL,S,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX + H_Hf_ABSORB
LSEL,U,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX
LSEL,U,LOC,Z,    COORD_Z_MAX + H_Hf_ABSORB,    COORD_Z_MAX + H_Hf_ABSORB
LESIZE,ALL,,,    N_Z,    1.0,    1    ,    ,    0
LSEL,ALL

*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
VSWEEP,    VOL1,    ,    ,    0
*GET, L_MAX, LINE, 0, NUM, MAX
LNUM_OLD = L_MAX
CM,Hf_ABSORBER_NIZ,VOLU
VSEL,ALL
COORD_Z_MAX = COORD_Z_MAX + H_Hf_ABSORB


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА W СНИЗУ (ФАКТИЧЕСКИЙ ОБЪЕМ 7-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_Hf_ABSORB/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_Hf_ABSORB/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_Hf_ABSORB/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_Hf_ABSORB/2,    COORD_Z_MAX
KSEL,S,LOC,Z,      COORD_Z_MAX,  COORD_Z_MAX
KSEL,S,KP,,    K_MAX+1,    K_MAX+5
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ W СНИЗУ (ОБЪЕМ 7-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj, jmin+1, jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_W_FOIL/2
        *IF,    jj+1,    EQ,    jmax,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_W_FOIL/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ W СНИЗУ (ОБЪЕМ 7-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Z,    COORD_Z_MAX,  COORD_Z_MAX
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1
AL,    ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА W СНИЗУ (ОБЪЕМ 7-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
VOLD = VNEW
VSEL,NONE
VEXT, A_MAX,,,    0.0, 0.0, H_W_FOIL,        0.0, 0.0, 0.0
VATT,    4,,1,,
*GET, L_MAX, LINE, 0, NUM, MAX
LSEL,S,LINE,, L_MAX-3,    L_MAX,1,
LESIZE,ALL,,,    N_Z/5,    1.0,    1    ,    ,    0
LSEL,ALL
    !LSEL,S,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX + H_W_FOIL
    !LSEL,U,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX
    !LSEL,U,LOC,Z,    COORD_Z_MAX + H_W_FOIL,    COORD_Z_MAX + H_W_FOIL
    !LESIZE,ALL,,,    N_Z/2,    1.0,    1    ,    ,    0
    !LSEL,ALL
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
VSWEEP,    VOL1,    ,    ,    0

LNUM_OLD = L_MAX
CM,W_FOIL_NIZ,VOLU
VSEL,ALL
COORD_Z_MAX = COORD_Z_MAX + H_W_FOIL

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА СЕРДЕЧНИКА (ФАКТИЧЕСКИЙ ОБЪЕМ 8-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_CORE/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_CORE/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_CORE/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_CORE/2,    COORD_Z_MAX
KSEL,S,LOC,Z,      COORD_Z_MAX,  COORD_Z_MAX
KSEL,S,KP,,    K_MAX+1,    K_MAX+5
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ СЕРДЕЧНИКА  (ОБЪЕМ 8-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj, jmin+1, jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_CORE/2
        *IF,    jj+1,    EQ,    jmax,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_CORE/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ СЕРДЕЧНИКА  (ОБЪЕМ 8-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Z,    COORD_Z_MAX,  COORD_Z_MAX
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1
AL,    ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА СЕРДЕЧНИКА (ОБЪЕМ 8-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
VOLD = VNEW
VSEL,NONE
VEXT, A_MAX,,,    0.0, 0.0, H_CORE,        0.0, 0.0, 0.0
VATT,    1,,1,,
LSEL,S,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX + H_CORE
LSEL,U,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX
LSEL,U,LOC,Z,    COORD_Z_MAX + H_CORE,    COORD_Z_MAX + H_CORE
LESIZE,ALL,,,    5*N_Z,    1.0,    1    ,    ,    0
LSEL,ALL

*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
VSWEEP,    VOL1,    ,    ,    0
*GET, L_MAX, LINE, 0, NUM, MAX
LNUM_OLD = L_MAX
CM,CORE,VOLU
VSEL,ALL
COORD_Z_MAX = COORD_Z_MAX + H_CORE

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ W СНИЗУ (ФАКТИЧЕСКИЙ ОБЪЕМ 9-ОЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_W_FOIL/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_W_FOIL/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_W_FOIL/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_W_FOIL/2,    COORD_Z_MAX
KSEL,S,LOC,Z,      COORD_Z_MAX,  COORD_Z_MAX
KSEL,S,KP,,    K_MAX+1,    K_MAX+5
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ W СВЕРХУ (ОБЪЕМ 9-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj, jmin+1, jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_W_FOIL/2
        *IF,    jj+1,    EQ,    jmax,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_W_FOIL/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ W СВЕРХУ (ОБЪЕМ 9-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Z,    COORD_Z_MAX,  COORD_Z_MAX
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1
AL,    ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА W СВЕРХУ (ОБЪЕМ 9-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
VOLD = VNEW
VSEL,NONE
VEXT, A_MAX,,,    0.0, 0.0, H_W_FOIL,        0.0, 0.0, 0.0
VATT,    4,,1,,
LSEL,S,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX + H_W_FOIL
LSEL,U,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX
LSEL,U,LOC,Z,    COORD_Z_MAX + H_CORE,    COORD_Z_MAX + H_W_FOIL
LESIZE,ALL,,,    N_Z/2,    1.0,    1    ,    ,    0
LSEL,ALL

*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
VSWEEP,    VOL1,    ,    ,    0
*GET, L_MAX, LINE, 0, NUM, MAX
LNUM_OLD = L_MAX
CM,W_FOIL_VERH,VOLU
CMSEL,A,W_FOIL_NIZ,VOLU
CM,W_FOIL_SUM,VOLU
VSEL,ALL
COORD_Z_MAX = COORD_Z_MAX + H_W_FOIL

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ Hf ПОГЛОТИТЕЛЬ СВЕРХУ (ФАКТИЧЕСКИЙ ОБЪЕМ 10-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_Hf_ABSORB/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_Hf_ABSORB/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_Hf_ABSORB/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_Hf_ABSORB/2,    COORD_Z_MAX
KSEL,S,LOC,Z,      COORD_Z_MAX,  COORD_Z_MAX
KSEL,S,KP,,    K_MAX+1,    K_MAX+5
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ Hf СВЕРХУ (ОБЪЕМ 11-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj, jmin+1, jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_Hf_ABSORB/2
        *IF,    jj+1,    EQ,    jmax,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_Hf_ABSORB/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ Hf СВЕРХУ (ОБЪЕМ 11-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Z,    COORD_Z_MAX,  COORD_Z_MAX
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1
AL,    ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА Hf СВЕРХУ (ОБЪЕМ 11-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
VOLD = VNEW
VSEL,NONE
VEXT, A_MAX,,,    0.0, 0.0, H_Hf_ABSORB,        0.0, 0.0, 0.0
VATT,    3,,1,,
LSEL,S,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX + H_Hf_ABSORB
LSEL,U,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX
LSEL,U,LOC,Z,    COORD_Z_MAX + H_CORE,    COORD_Z_MAX + H_Hf_ABSORB
LESIZE,ALL,,,    N_Z/2,    1.0,    1    ,    ,    0
LSEL,ALL
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
VSWEEP,    VOL1,    ,    ,    0
*GET, L_MAX, LINE, 0, NUM, MAX
LNUM_OLD = L_MAX
CM,Hf_ABSORBER_VERH,VOLU
CMSEL,A,Hf_ABSORBER_NIZ,VOLU
CM,Hf_ABSORBER_SUM,VOLU
VSEL,ALL
COORD_Z_MAX = COORD_Z_MAX + H_Hf_ABSORB

 

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ЦИКЛИЧНО ПОВТОРЯЮЩИХСЯ  ИМИТАТОРОВ _N_TIME2__РАЗ ОТ НИЖНЕЙ ЗАГЛУШКИ (ФАКТИЧЕСКИЕ ОБЪЕМЫ 3-5) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
VSEL,NONE
CM,Zr_SIMULATOR_TOP_SUM,VOLU
CM,GAS_PLENUM_TOP_SUM,VOLU
VSEL,ALL
*DO,jk,0,N_TIME2-1

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ОПОРНЫХ ТОЧЕК ПОПЕРЕЧНОГО СЕЧЕНИЯ ЗАГОТОВКИ  ИМИТАТОРА №2 (ОБЪЕМ 5-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_Zr_SIM_OUT/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_Zr_SIM_OUT/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_Zr_SIM_OUT/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_Zr_SIM_OUT/2,    COORD_Z_MAX
KSEL,S,LOC,Z,      COORD_Z_MAX,  COORD_Z_MAX
KSEL,S,KP,,    K_MAX+1,    K_MAX+5
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ ЗАГОТОВКИ  ИМИТАТОРА №2 (ОБЪЕМ 5-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj, jmin+1, jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_Zr_SIM_OUT/2
        *IF,    jj+1,    EQ,    jmax,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_Zr_SIM_OUT/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
LNUM_OLD = L_MAX        ! ПЕРЕПРИСВОЕНИЕ МАКСИМАЛЬНОГО НОМЕРА ЛИНИИ НА ПРЕДЫДУЩЕМ ШАГЕ
*GET, L_MAX, LINE, 0, NUM, MAX    ! ЗАПРОС МАКСИМАЛЬНОГО НОМЕРА ЛИНИИ
L_NUM = L_MAX            ! ПРИСВОЕНИЕ МАКСИМАЛЬНОГО НОМЕРА ЛИНИИ НА ТЕКУЩЕМ ШАГЕ

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ ЗАГОТОВКИ  ИМИТАТОРА №2 (ОБЪЕМ 5-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Z,    COORD_Z_MAX,  COORD_Z_MAX
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1
AL,    ALL
!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА ЗАГОТОВКИ  ИМИТАТОРА №2  (ОБЪЕМ 5-ЫЙ) !!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
VOLD = VNEW
VEXT, A_MAX,,,    0.0, 0.0, H_Zr_SIM2,        0.0, 0.0, 0.0
VATT,    2,,1,,
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
VSEL,U,VOLU,,1,VOLD,1,
!VSWEEP,    VOL1,    ,    ,    0
*GET, L_MAX, LINE, 0, NUM, MAX
LNUM_OLD = L_MAX
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ОПОРНЫХ ТОЧЕК ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В ЦИРКОНИЕВОМ ИМИТАТОРЕ №2 (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        H_Zr_SIM2+COORD_Z_MAX
K,    K_MAX+2,    D_Zr_SIM_IN/2,    0.0,    H_Zr_SIM2+COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_Zr_SIM_IN/2,    H_Zr_SIM2+COORD_Z_MAX
K,    K_MAX+4,    -D_Zr_SIM_IN/2,    0.0,    H_Zr_SIM2+COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_Zr_SIM_IN/2,    H_Zr_SIM2+COORD_Z_MAX
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В ЦИРКОНИЕВОМ ИМИТАТОРЕ №2 (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj, jmin+1, jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_Zr_SIM_IN/2
        *IF,    jj+1,    EQ,    jmax,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_Zr_SIM_IN/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В ЦИРКОНИЕВОМ ИМИТАТОРЕ №2 (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
LSEL,S,LOC,Z,    H_Zr_SIM2+COORD_Z_MAX,  H_Zr_SIM2+COORD_Z_MAX
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1
AL,    ALL
!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА ОТВЕРСТИЯ МАЛОГО  ИМИТАТОРА (ОБЪЕМ 3-ИЙ) !!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
WPOFFS,    0.0,    0.0,    H_Zr_SIM2+COORD_Z_MAX
CSYS,4
VOLD = VNEW
VEXT, A_MAX,,,    0.0, 0.0, -H_Zr_SIM2,        0.0, 0.0, 0.0
VATT,    5,,1,,
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
VSEL,U,VOLU,,1,VOLD,1,
CM,GAS_PLENUM_TOP,VOLU
CMSEL,A,GAS_PLENUM_TOP_SUM,VOLU
CM,GAS_PLENUM_TOP_SUM,VOLU
VSEL,ALL
*GET, L_MAX, LINE, 0, NUM, MAX
LNUM_OLD = L_MAX
CSYS,0  
WPAVE,0,0,0
CSYS,4
!!!!!!!!!!!!!!!!!!! ВЫЧИТАНИЕ ИЗ ЗАГОТОВКИ ОТВЕРСТИЯ (ОБЪЕМ 7-ОЙ) - СОЗДАНИЕ ОБЪЕМА  ИМИТАТОРА №2 !!!!!!!!!!!!!!!!!!!
VSBV,    VOLD,    VNEW,    ,    KEEP,    KEEP
LSEL,S,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX + H_Zr_SIM2
LSEL,U,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX
LSEL,U,LOC,Z,    COORD_Z_MAX + H_Zr_SIM2,    COORD_Z_MAX + H_Zr_SIM2
LESIZE,ALL,,,    N_Z,    1.0,    1    ,    ,    0
LSEL,ALL
VSWEEP,    VNEW,    A_MAX,    A_MAX+1,    0

NUMCMP,AREA
VOLD = VNEW
*GET, VOL1,VOLU,0,NUM,MAX
*GET, A_MAX,AREA,0,NUM,MAX
VNEW = VOL1
VSEL,U,VOLU,,1,VOLD,1,
CM,Zr_SIMULATOR_TOP,VOLU
CMSEL,A,Zr_SIMULATOR_TOP_SUM,VOLU
CM,Zr_SIMULATOR_TOP_SUM,VOLU
WPOFFS,    0.0,    0.0,    H_Zr_SIM2+COORD_Z_MAX
CSYS,0
VATT,    2,,1,,
VSWEEP,    VNEW,    A_MAX-1,    A_MAX,    0
VSEL,ALL
COORD_Z_MAX = COORD_Z_MAX + H_Zr_SIM2
CSYS,0  
WPOFFS,    0.0,    0.0,    -COORD_Z_MAX
CSYS,4
*ENDDO


!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА МАЛОГО  ИМИТАТОРА СВЕРХУ (ФАКТИЧЕСКИЙ ОБЪЕМ 10-ЫЙ) !!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_Zr_SIM_OUT/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_Zr_SIM_OUT/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_Zr_SIM_OUT/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_Zr_SIM_OUT/2,    COORD_Z_MAX
KSEL,S,LOC,Z,      COORD_Z_MAX,  COORD_Z_MAX
KSEL,S,KP,,    K_MAX+1,    K_MAX+5
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ МАЛОГО  ИМИТАТОРА СВЕРХУ СВЕРХУ (ОБЪЕМ 11-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj, jmin+1, jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_Zr_SIM_OUT/2
        *IF,    jj+1,    EQ,    jmax,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_Zr_SIM_OUT/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ МАЛОГО  ИМИТАТОРА СВЕРХУ (ОБЪЕМ 11-ЫЙ) !!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Z,    COORD_Z_MAX,  COORD_Z_MAX
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1
AL,    ALL
!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА МАЛОГО  ИМИТАТОРА СВЕРХУ (ОБЪЕМ 11-ЫЙ) !!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
VOLD = VNEW
VSEL,NONE
VEXT, A_MAX,,,    0.0, 0.0, H_Zr_SIM3,        0.0, 0.0, 0.0
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
*GET, L_MAX, LINE, 0, NUM, MAX
LNUM_OLD = L_MAX
!CM,Hf_ABSORBER_VERH,VOLU
VSEL,ALL
COORD_Z_MAX = COORD_Z_MAX + H_Zr_SIM3
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ОПОРНЫХ ТОЧЕК ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В МАЛОМ ИМИТАТОРЕ  СВЕРХУ (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_Zr_SIM_IN/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_Zr_SIM_IN/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_Zr_SIM_IN/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_Zr_SIM_IN/2,    COORD_Z_MAX
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В МАЛОМ ИМИТАТОРЕ  СВЕРХУ (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj, jmin+1, jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_Zr_SIM_IN/2
        *IF,    jj+1,    EQ,    jmax,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_Zr_SIM_IN/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В МАЛОМ ИМИТАТОРЕ  СВЕРХУ (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
LSEL,ALL
LSEL,S,LOC,Z,    COORD_Z_MAX,  COORD_Z_MAX
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1
AL,    ALL
!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА ОТВЕРСТИЯ МАЛОГО  ИМИТАТОРА СВЕРХУ (ОБЪЕМ 3-ИЙ) !!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
WPOFFS,    0.0,    0.0,    COORD_Z_MAX
CSYS,4
VOLD = VNEW
VEXT, A_MAX,,,    0.0, 0.0, -H_Zr_SIM3,        0.0, 0.0, 0.0
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
VSEL,U,VOLU,,1,VOLD,1,
VATT,    5,,1,,
CM,GAS_PLENUM_TOP,VOLU
CMSEL,A,GAS_PLENUM_TOP_SUM,VOLU
CM,GAS_PLENUM_TOP_SUM,VOLU
VSEL,ALL
*GET, L_MAX, LINE, 0, NUM, MAX
LNUM_OLD = L_MAX
!!!!!!!!!!!!!!!!!!! ВЫЧИТАНИЕ ИЗ ЗАГОТОВКИ ОТВЕРСТИЯ (ОБЪЕМ 7-ОЙ) - СОЗДАНИЕ ОБЪЕМА  ИМИТАТОРА №2 !!!!!!!!!!!!!!!!!!!
VSBV,    VOLD,    VNEW,    ,    KEEP,    KEEP

LSEL,S,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX + H_Zr_SIM3
LSEL,U,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX
LSEL,U,LOC,Z,    COORD_Z_MAX + H_Zr_SIM3,    COORD_Z_MAX + H_Zr_SIM3
LESIZE,ALL,,,    N_Z,    1.0,    1    ,    ,    0
LSEL,ALL
VSWEEP,    VNEW,    A_MAX,    A_MAX+1,    0
NUMCMP,AREA
VOLD = VNEW
*GET, VOL1,VOLU,0,NUM,MAX
*GET, A_MAX,AREA,0,NUM,MAX
VNEW = VOL1
VSEL,U,VOLU,,1,VOLD,1,
CM,Zr_SIMULATOR_TOP,VOLU
CMSEL,A,Zr_SIMULATOR_TOP_SUM,VOLU
CM,Zr_SIMULATOR_TOP_SUM,VOLU
!WPOFFS,    0.0,    0.0,    H_Zr_SIM3+COORD_Z_MAX
CSYS,0
VATT,    2,,1,,
VSWEEP,    VNEW,    A_MAX-1,    A_MAX,    0
VSEL,ALL
!COORD_Z_MAX = COORD_Z_MAX + H_Zr_SIM3
CSYS,0  
WPOFFS,    0.0,    0.0,    -COORD_Z_MAX
CSYS,4

 

!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЛАСТИ ГАЗОСБОРНИКА (ФАКТИЧЕСКИЙ ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_CLAD_IN/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_CLAD_IN/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_CLAD_IN/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_CLAD_IN/2,    COORD_Z_MAX
KSEL,S,LOC,Z,  COORD_Z_MAX,    COORD_Z_MAX
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ МАЛОГО  ИМИТАТОРА (ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj,jmin+1,jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_CLAD_IN/2
        *IF,    jj,    EQ,    jmax-1,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_CLAD_IN/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ МАЛОГО  ИМИТАТОРА (ОБЪЕМ 2-ОЙ)  !!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1,    
AL,    ALL
!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМНОЙ ЗАГОТОВКИ МАЛОГО  ИМИТАТОРА (ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
VOLD = VNEW
VEXT, A_MAX,,,    0.0, 0.0, L_FREE,        0.0, 0.0, 0.0
VATT,    5,,1,,
LSEL,S,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX +  L_FREE
LSEL,U,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX
LSEL,U,LOC,Z,    COORD_Z_MAX +  L_FREE,    COORD_Z_MAX +  L_FREE
LESIZE,ALL,,,    N_Z,    1.0,    1    ,    ,    0
LSEL,ALL
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
VSEL,U,VOLU,,1,VOLD,1,
CM,GAS_PLENUM_TOP,VOLU
CMSEL,A,GAS_PLENUM_TOP_SUM,VOLU
CM,GAS_PLENUM_TOP_SUM,VOLU
*GET, L_MAX, LINE, 0, NUM, MAX
VMESH, VOL1
VSEL,ALL
COORD_Z_MAX = COORD_Z_MAX + L_FREE


!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА ВЕРХНЕЙ ЗАГЛУШКИ (ФАКТИЧЕСКИЙ ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_END_PLUG_TOP/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_END_PLUG_TOP/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_END_PLUG_TOP/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_END_PLUG_TOP/2,    COORD_Z_MAX
KSEL,S,LOC,Z,  COORD_Z_MAX,    COORD_Z_MAX
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ МАЛОГО  ИМИТАТОРА (ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj,jmin+1,jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_END_PLUG_TOP/2
        *IF,    jj,    EQ,    jmax-1,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_END_PLUG_TOP/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ МАЛОГО  ИМИТАТОРА (ОБЪЕМ 2-ОЙ)  !!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1,    
AL,    ALL
!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМНОЙ ЗАГОТОВКИ МАЛОГО  ИМИТАТОРА (ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
VOLD = VNEW
VEXT, A_MAX,,,    0.0, 0.0, H_END_PLUG_TOP,        0.0, 0.0, 0.0
VATT,    2,,1,,
LSEL,S,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX+H_END_PLUG_TOP
LSEL,U,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX
LSEL,U,LOC,Z,    COORD_Z_MAX+H_END_PLUG_TOP,        COORD_Z_MAX+H_END_PLUG_TOP
LESIZE,ALL,,,    N_Z,    1.0,    1    ,    ,    0
LSEL,ALL
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
VSEL,U,VOLU,,1,VOLD,1,
CM,END_PLUG_TOP_VERH,VOLU
*GET, L_MAX, LINE, 0, NUM, MAX
VMESH, VOL1
VSEL,ALL
COORD_Z_MAX = COORD_Z_MAX + H_END_PLUG_TOP

 

 


!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ЗАГОТОВКИ ОБЪЕМА ОБОЛОЧКИ (ФАКТИЧЕСКИЙ ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_CLAD_OUT/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_CLAD_OUT/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_CLAD_OUT/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_CLAD_OUT/2,    COORD_Z_MAX
KSEL,S,LOC,Z,  COORD_Z_MAX,    COORD_Z_MAX
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ ЗАГОТОВКИ ОБЪЕМА ОБОЛОЧКИ (ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj,jmin+1,jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_CLAD_OUT/2
        *IF,    jj,    EQ,    jmax-1,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_CLAD_OUT/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ ЗАГОТОВКИ ОБЪЕМА ОБОЛОЧКИ (ОБЪЕМ 2-ОЙ)  !!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1,    
AL,    ALL
!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМНОЙ ЗАГОТОВКИ ОБОЛОЧКИ (ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
VOLD = VNEW
VEXT, A_MAX,,,    0.0, 0.0, -COORD_Z_MAX,        0.0, 0.0, 0.0
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
*GET, L_MAX, LINE, 0, NUM, MAX
!VMESH, VOL1
COORD_Z_MAX = COORD_Z_MAX - COORD_Z_MAX

!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ОПОРНЫХ ТОЧЕК ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В ЗАГОТОВКЕ ОБЪЕМА ОБОЛОЧКИ (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_CLAD_IN/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_CLAD_IN/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_CLAD_IN/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_CLAD_IN/2,    COORD_Z_MAX
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В ЗАГОТОВКЕ ОБЪЕМА ОБОЛОЧКИ (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj, jmin+1, jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_CLAD_IN/2
        *IF,    jj+1,    EQ,    jmax,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_CLAD_IN/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В ЗАГОТОВКЕ ОБЪЕМА ОБОЛОЧКИ (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
LSEL,S,LOC,Z,    COORD_Z_MAX,  COORD_Z_MAX
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1
AL,    ALL
!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА ОТВЕРСТИЯ В ЗАГОТОВКЕ ОБЪЕМА ОБОЛОЧКИ (ОБЪЕМ 3-ИЙ) !!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
WPOFFS,    0.0,    0.0,    COORD_Z_MAX - COORD_Z_MAX
CSYS,4
VOLD = VNEW
VEXT, A_MAX,,,    0.0, 0.0, L_DIM_MAX,        0.0, 0.0, 0.0
VATT,    5,,1,,
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
*GET, L_MAX, LINE, 0, NUM, MAX
LNUM_OLD = L_MAX
COORD_Z_MAX = COORD_Z_MAX + L_DIM_MAX
!!!!!!!!!!!!!!!!!!! ВЫЧИТАНИЕ ИЗ ЗАГОТОВКИ ОТВЕРСТИЯ (ОБЪЕМ 7-ОЙ) - СОЗДАНИЕ ОБЪЕМА  ИМИТАТОРА №2 !!!!!!!!!!!!!!!!!!!
VSBV,    VOLD,    VNEW,    ,    KEEP,    KEEP
NUMCMP,AREA
VOLD = VNEW
*GET, VOL1,VOLU,0,NUM,MAX
*GET, A_MAX,AREA,0,NUM,MAX
VNEW = VOL1
VSEL,U,VOLU,,1,VOLD,1,
VATT,    2,,1,,
CM,Zr_CLAD,VOLU
CSYS,0
VSWEEP,    VNEW,    A_MAX-1,    A_MAX,    0
VSEL,ALL
WPOFFS,    0.0,    0.0,    H_END_PLUG_BOT
LNUM_OLD = L_MAX
WPAVE,0,0,0

!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ЗАГОТОВКИ ОБЪЕМА ОБОЛОЧКИ (ФАКТИЧЕСКИЙ ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        H_END_PLUG_BOT
K,    K_MAX+2,    D_CLAD_IN/2,    0.0,    H_END_PLUG_BOT
K,    K_MAX+3,    0.0,    D_CLAD_IN/2,    H_END_PLUG_BOT
K,    K_MAX+4,    -D_CLAD_IN/2,    0.0,    H_END_PLUG_BOT
K,    K_MAX+5,    0.0,    -D_CLAD_IN/2,    H_END_PLUG_BOT
KSEL,S,LOC,Z,  H_END_PLUG_BOT,    H_END_PLUG_BOT
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ ЗАГОТОВКИ ОБЪЕМА ОБОЛОЧКИ (ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj,jmin+1,jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_CLAD_IN/2
        *IF,    jj,    EQ,    jmax-1,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_CLAD_IN/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ ЗАГОТОВКИ ОБЪЕМА ОБОЛОЧКИ (ОБЪЕМ 2-ОЙ)  !!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1,    
AL,    ALL
!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМНОЙ ЗАГОТОВКИ ОБОЛОЧКИ (ОБЪЕМ 2-ОЙ) !!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
VOLD = VNEW
VEXT, A_MAX,,,    0.0, 0.0, (L_DIM_MAX - H_END_PLUG_BOT - H_END_PLUG_TOP - L_FREE),        0.0, 0.0, 0.0
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
*GET, L_MAX, LINE, 0, NUM, MAX
!VMESH, VOL1
COORD_Z_MAX = COORD_Z_MAX - (H_END_PLUG_TOP + L_FREE)

!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ОПОРНЫХ ТОЧЕК ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В ЗАГОТОВКЕ ОБЪЕМА ОБОЛОЧКИ (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
*GET, K_MAX, KP, 0, NUM, MAX
K,    K_MAX+1,    0.0,    0.0,        COORD_Z_MAX
K,    K_MAX+2,    D_CORE/2,    0.0,    COORD_Z_MAX
K,    K_MAX+3,    0.0,    D_CORE/2,    COORD_Z_MAX
K,    K_MAX+4,    -D_CORE/2,    0.0,    COORD_Z_MAX
K,    K_MAX+5,    0.0,    -D_CORE/2,    COORD_Z_MAX
!!!!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ЛИНИЙ ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В ЗАГОТОВКЕ ОБЪЕМА ОБОЛОЧКИ (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
SLAB2 = 0
jmin = K_MAX+1
jmax = (K_MAX+4)+1
*DO,jj, jmin+1, jmax-1
    LARC,    jj,    jj+1,    K_MAX+1,    D_CORE/2
        *IF,    jj+1,    EQ,    jmax,    THEN
            LARC,    jmax,    jmin+1,    K_MAX+1,    D_CORE/2
            SLAB = SLAB +1
            SLAB2 = SLAB2 +1
        *ENDIF
    SLAB = SLAB +1
    SLAB2 = SLAB2 +1
*ENDDO
!!!!!!!!!!!!!!!!!!! ПОСТРОЕНИЕ ПОПЕРЕЧНОГО СЕЧЕНИЯ ОТВЕРСТИЯ В ЗАГОТОВКЕ ОБЪЕМА ОБОЛОЧКИ  (ОБЪЕМ 6-ОЙ) !!!!!!!!!!!!!!!!!!!!!!
LNUM_OLD = L_MAX
*GET, L_MAX, LINE, 0, NUM, MAX
L_NUM = L_MAX

LSEL,S,LINE,,    LNUM_OLD+1,        L_NUM,    1,
LESIZE,ALL,,,    N_PHI,    1.0,    1    ,    ,    0
LSEL,S,LOC,Z,    COORD_Z_MAX,  COORD_Z_MAX
LSEL,S,LINE,,    LNUM_OLD+1,    L_NUM,    1
AL,    ALL
!!!!!!!!!!!!!!!!!!! СОЗДАНИЕ ОБЪЕМА ОТВЕРСТИЯ В ЗАГОТОВКЕ ОБЪЕМА ОБОЛОЧКИ (ОБЪЕМ 3-ИЙ) !!!!!!!!!!!!!!!!!!!
*GET, A_MAX, AREA, 0, NUM, MAX
!WPOFFS,    0.0,    0.0,    COORD_Z_MAX - COORD_Z_MAX+H_END_PLUG_BOT
CSYS,4
VOLD = VNEW
VEXT, A_MAX,,,    0.0, 0.0, -(COORD_Z_MAX-H_END_PLUG_BOT),        0.0, 0.0, 0.0
VATT,    5,,1,,
*GET, VOL1,VOLU,0,NUM,MAX
VNEW = VOL1
*GET, L_MAX, LINE, 0, NUM, MAX
LNUM_OLD = L_MAX
!COORD_Z_MAX = COORD_Z_MAX + L_DIM_MAX
!!!!!!!!!!!!!!!!!!! ВЫЧИТАНИЕ ИЗ ЗАГОТОВКИ ОТВЕРСТИЯ (ОБЪЕМ 7-ОЙ) - СОЗДАНИЕ ОБЪЕМА  ИМИТАТОРА №2 !!!!!!!!!!!!!!!!!!!
VSBV,    VOLD,    VNEW,    ,    KEEP,    KEEP
NUMCMP,AREA
VOLD = VNEW
*GET, VOL1,VOLU,0,NUM,MAX
*GET, A_MAX,AREA,0,NUM,MAX
VNEW = VOL1
VSEL,U,VOLU,,1,VOLD,1,
VATT,    5,,1,,
CM,GAS_GAP,VOLU
CSYS,0
VSWEEP,    VNEW,    A_MAX-1,    A_MAX,    0
VSEL,ALL
LNUM_OLD = L_MAX


CYL4,    0.0,    0.0,    D_CORE/2,    0.0,    D_CLAD_IN/2,    360.0,    (L_DIM_MAX - H_END_PLUG_BOT - H_END_PLUG_TOP - L_FREE)
*GET, L_MAX, LINE, 0, NUM, MAX
LNUM = L_MAX
VOLD = VNEW
*GET, VOL1,VOLU,0,NUM,MAX
!*GET, A_MAX,AREA,0,NUM,MAX
VNEW = VOL1
VSEL,U,VOLU,,1,VOLD,1,
VATT,5,,1,,
CM,GAS_GAP,VOLU
LSEL,S,LINE,,    LNUM_OLD+1,    LNUM,
LSEL,U,LOC,Z,    H_END_PLUG_BOT,    H_END_PLUG_BOT
LSEL,U,LOC,Z,    H_END_PLUG_BOT + (L_DIM_MAX - H_END_PLUG_BOT - H_END_PLUG_TOP - L_FREE),    H_END_PLUG_BOT + (L_DIM_MAX - H_END_PLUG_BOT - H_END_PLUG_TOP - L_FREE)
LESIZE,ALL,,,    5*N_Z,    1.0,    1    ,    ,    0
LSEL,ALL
VSWEEP,    VNEW,    A_MAX+1,    A_MAX+2,    0
VSEL,ALL
WPAVE,0,0,0

ASEL,S,LOC,Z,    0.0,    0.0
ASEL,A,LOC,Z,    COORD_Z_MAX,    COORD_Z_MAX
SFA,ALL,,        HFLUX,    0.0,    0.0
ASEL,ALL

CSYS,1
ASEL,S,LOC,X,    D_CLAD_OUT/2,    D_CLAD_OUT/2
SFA,ALL,    ,    CONV,    30.0E3,    (325+273)
CSYS,0
ASEL,ALL


CMSEL,S,CORE,VOLU
BFV,ALL,HGEN,    4*Q_LIN/(3.14159*D_CORE**2)

CMSEL,S,Zr_CLAD,VOLU
CSYS,1
ASEL,S,LOC,X,    D_CLAD_OUT/2,    D_CLAD_OUT/2
SFA,ALL,    ,    CONV,    30.0E3,    (325+273)
CSYS,0
ALLSEL,ALL

SFTRAN
BFTRAN

FINISH

/SOL
SOLVE

 

Ссылка на сообщение
Поделиться на других сайтах
53 минуты назад, piden сказал:

Блин, ощущение, что сейчас запущу этот манускрипт - и он мне операционку переставит...

На лаконичность и изящность программирования не претендую. Еду как могу :scratch_one-s_head:

Ссылка на сообщение
Поделиться на других сайтах
20 hours ago, bubblyk said:

И не считает!

Как на счет попробовать решить задачу сначала с одним телом?  А то решатель жалуется, что некоторые поверхности не замешены, нагрузки странно приложены.. Сделайте, чтобы считался простой кейс, потом добавляйте сложности.

 

Ну или живите надеждой, что кто-то эту простыню осилит и укажет точно, где затык или опечатка.

 

 

21 hours ago, bubblyk said:

цилиндрическое тело

Почему осевую симметрию не используете? 

Ссылка на сообщение
Поделиться на других сайтах
21 час назад, bubblyk сказал:

Коллеги, приветствую всех!

 

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

Задача моделируется в 3d (хотя и можно было пользуясь осесимметричностью решать 2d)

 

По классике - построил солид тела, разбил их mesh'ером с использованием solid278... Задал ГУ, казалось бы все просто ... И не считает! Долго бился, но так и не нашел причину сбоя.

 

Ответ мне заранее известен - сделал расчет по одномерным моделям (расчетные формулы из курса теплопроводности).

 

Помогите, пожалуйста, найти ошибку.

 

 

И да... Пишу на форум впервые, могу тормозить :bash:

Это я написал в первом сообщении. Понимаю, что задача решалась бы куда проще.

 

Пошел на 3d модель, думая о перспективе перестройки - приложение ГУ, изменяющихся в азимутальном направлении.

Ссылка на сообщение
Поделиться на других сайтах

@bubblyk , если перед solve добавить

pivcheck,off

то посчитается... несмотря ни на что и вопреки всему.

 

Но у вас сетка на телах некомфорная, контактов нету + все сказанное выше. Итог:

image.png

Ссылка на сообщение
Поделиться на других сайтах

@ piden

 

Благодарю за ответ!

 

1 Прошу прощения, но Ваш результат расчета в WB неверен!!! :blush2:

 

2 То, что теплопередача не идет по телам, я успел найти параллельно. Первопричину ошибки понял.

Верно я понимаю, что сетки, построенные на разных телах, нужно было "сцепливать" (VGLUE команда).

 

3 Некомфорная сетка - это как?  В этой терминологии я не силен :boredom:

 

4. Новый опыт - команда pivcheck,off  в новинку... полезно :good:

Ссылка на сообщение
Поделиться на других сайтах
35 minutes ago, bubblyk said:

Прошу прощения, но Ваш результат расчета в WB неверен!!! :blush2:

Модель - ваша, а результат - мой?)

 

36 minutes ago, bubblyk said:

Некомфорная сетка - это как?

"Не сцепленная"

 

37 minutes ago, bubblyk said:

Верно я понимаю, что сетки, построенные на разных телах, нужно было "сцепливать" (VGLUE команда).

vglue - это для объемов. Но суть такая, да.

Ссылка на сообщение
Поделиться на других сайтах
10 минут назад, piden сказал:

Модель - ваша, а результат - мой?)

 

"Не сцепленная"

 

vglue - это для объемов. Но суть такая, да.

Моя модель рисуется иначе ) и ГУ прикладываются к другим телам - прорисовывал в MAPDL с выводом на экран, проверял.

Одно из центральных тел должно быть самым горячим

Ссылка на сообщение
Поделиться на других сайтах
25 минут назад, bubblyk сказал:

Моя модель рисуется иначе )

не туда упорство.

25 минут назад, bubblyk сказал:

прорисовывал в MAPDL с выводом на экран, проверял.

уже почти 10 публикаций, лейте картинки. но я верю на слово, что выше именно Ваша сетка, всего лишь импортированная в ВБ.

1 час назад, piden сказал:

некомфорная

разрешите буковку исправить.

 

28 минут назад, bubblyk сказал:

Одно из центральных тел должно быть самым горячим

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

Ссылка на сообщение
Поделиться на других сайтах

Разобрался с задачей, первоначально выполнив 2d расчет в осесимметричной постановке

 

 

Скрытый текст

 

Вскоре реализую 3d моделирование

Ссылка на сообщение
Поделиться на других сайтах

Присоединяйтесь к обсуждению

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

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставить изображения напрямую. Загрузите или вставьте изображения по ссылке.

  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу.




  • Сообщения

    • maxx2000
      Причина того - Кроилово. Кроилово всегда приводит к попадалову. Месяц простоял сколько мильонов деревянных потеряли? Вопрос риторический. И ещё будет стоять. Как памятник человеческой глупости и жадности.
    • AlexKaz
      "9 июля 1968 года на мышах был проведен самый знаменитый эксперимент американского ученого-этолога Джона Кэлхуна «Вселенная-25». Суть опыта заключалась в создании идеальных условий, где мыши могли бы жить и размножаться, не ведая никаких забот, вдали от хищников и в отсутствие эпидемий и заболеваний. Для этих целей ученый построил специальный загон, куда были помещены четыре пары белых мышей (самцов и самок). В распоряжении мышей всегда была чистая вода и еда в изобилии, специальные гнезда, где можно обустроить себе жилище ― гнезд в загоне хватало для проживания нескольких тысяч мышей. Температура в загоне в среднем составляла около 20 ℃ и была комфортной для мышей. Животные не подвергались никаким влияниям извне и жили в идеальных условиях в свое удовольствие. А дальше началось самое интересное. На первом этапе эксперимента мыши хорошо размножались, вели активный образ жизни, охотно играли. На следующей фазе эксперимента мыши стали есть меньше, перестали наедаться до отвала. На третьей фазе эксперимента, когда в загоне были уже сотни мышей, произошло распределение социальных ролей, стала ярко выраженной иерархия, клановость. Появились так называемые отверженные ― молодые особи, которых другие, взрослые мыши сгоняли в центр загона, не давали им вести нормальный образ жизни, причиняли физический вред. В природе такое, наверное, было бы невозможно, ведь эти мыши-агрессоры просто не дожили бы до старости: их бы съели хищники. Но в загоне Кэлхуна хищников не было, и взрослые мыши начали попросту издеваться над молодняком. Образовались две большие группировки: самцы-одиночки и самки-одиночки. При этом самки-одиночки отказывались спариваться <с менее статусными многочисленными молодыми самцами и с оставшимися старыми статусными> и отвергали ухаживания самцов. У мышей стал проявляться тотальный индивидуализм, мыши не стремились создать семью. На последней, четвертой стадии мышиная популяция стала сокращаться. Появились самцы, которых сам Кэлхун назвал «красивыми» (англ. beautiful ones), из-за отсутствия ран и рубцов. <В оригинале: They never engaged in sexual approaches toward females, and they never engaged in fighting, and so they had no wound or scar tissue. Thus their pelage remained in excellent condition. - Дословный перевод: Они никогда не прибегали к сексуальным подходам к самкам, и они никогда не участвовали в боях, и поэтому у них не было ран или рубцовой ткани. Таким образом, их шерсть сохранилась в отличном состоянии.> Эти мыши не вступали в борьбу за самок и территорию, не проявляли активности к размножению и только питались, спали и чистили шёрстку. У мышей стали проявляться различные формы девиантного поведения, вспышки агрессии. Самки стали проявлять агрессию, защищать себя сами, стали умерщвлять своих детенышей, а затем окончательно отказались размножаться. На пике эксперимента в загоне одновременно проживало чуть более двух тыс. мышей. Еды и гнезд было достаточно для дальнейшего роста популяции, но через четыре года после начала эксперимента Кэлхун остановил свой опыт, потому что в загоне осталось чуть более сотни мышей, и все они уже вышли из репродуктивного возраста. По итогам эксперимента Кэлхун пришел к выводу, что достижение определенной плотности населения и заполнение социальных ролей в популяции приводит к распаду общества" https://physicsoflife.pl/dict/pic/calhoun/calhoun.. https://scientificrussia.ru/articles/utopiya-dlya-mys.. https://ru.wikipedia.org/wiki/Кэлхун,_Джон_(этолог)
    • gudstartup
      @Koels вот в чем дело пока ds609 это предупреждение поэтому F может и не появится если sv601 это значит ошибка. возможно при нагреве радиатора серво определяет это как предупреждение или ваш вентилятор крутиться медленнее чем оригинальный и серва думает что он встал хотяпри этом обычно на экране в строке состояния FAN.мигает больше у меня вариантов нет....  
    • ДОБРЯК
      Решите любым алгоритмом. Тогда будет конструктивный разговор. :=)
    • Fedor
      https://en.wikipedia.org/wiki/List_of_numerical_analysis_topics#Eigenvalue_algorithms     :) 
    • Юрий К.Ф.
      Добрый день. Не нашёл тут тему по стойке Китайско Китайской)) Lynuc N3ME. Видать мне так повезло с её наличием)) Приобрели 5-ти осевой Китаец. В б/у состоянии после удара по оси Z. Отремонтировали по механике, заменили батареи на драйверах, выставили лимиты. Всё Ок. Пригласили со стороны людей которые бы разобрались по операторской части. Те два выходных ковырялись, после сказали покажут расскажут, но за огромные деньги. Не сошлись. После месяц станок простоял, когда включили перестал реагировать на регулировку скорости шпинделя. То есть в режиме Jog, включаем обороты, которые стандартно 2140-2149 об/мин. При регулировке процетности не меняются (сама процентность показывает на мониторе). Так же при включении оборотов через команду M03S300 или другое значение, скорость так же показывает 2140-2149 об/мин. Грешить на тех товарищей с которыми не сошлись по деньгам для обучения, как то не хочется. Поковырялся в настройках шпинделя, вроде всё в норме. Проводку на шпинделе прозвонил, целая. В чём причина, не понятна. Кто нибудь сталкивался с подобным, или с подобной стойкой? Может подсказать варианты причины подобного?
    • ДОБРЯК
    • Fedor
      https://en.wikipedia.org/wiki/Eigenvalue_algorithm   Вот , кстати, важное свойство которое раньше встречал у Стренга. :)    http://elibrary.bsu.edu.az/files/kitablar/1022.pdf   
    • ДОБРЯК
      Вы ысё чаще напоминаете телевизор, который случайным образом переключает программы...:=)
    • sssdos
      сочувствую бедности(Дома на ноуте стоит Windows XP С PM 6.). Или другие причины?
×
×
  • Создать...