Jump to content

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


Recommended Posts

bubblyk

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

 

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

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

 

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

 

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

 

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

 

 

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

  • Нравится 1
Link to post
Share on other sites


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

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

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

Link to post
Share on other sites
bubblyk

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

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

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

Link to post
Share on other sites
piden
21 minutes ago, bubblyk said:

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

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

image.png

Link to post
Share on other sites
bubblyk
13 минут назад, piden сказал:

image.png

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

Не допёр...

Link to post
Share on other sites
piden
1 minute ago, bubblyk said:

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

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

 

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

 

17 hours ago, bubblyk said:

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

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

 

Link to post
Share on other sites
bubblyk

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

 

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

/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

 

Link to post
Share on other sites
piden

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

  • Нравится 1
Link to post
Share on other sites
bubblyk
53 минуты назад, piden сказал:

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

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

  • Нравится 1
Link to post
Share on other sites
piden
20 hours ago, bubblyk said:

И не считает!

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

 

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

 

 

21 hours ago, bubblyk said:

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

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

Link to post
Share on other sites
bubblyk
21 час назад, bubblyk сказал:

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

 

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

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

 

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

 

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

 

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

 

 

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

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

 

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

Link to post
Share on other sites
piden

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

pivcheck,off

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

 

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

image.png

Link to post
Share on other sites
bubblyk

@ piden

 

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

 

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

 

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

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

 

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

 

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

Link to post
Share on other sites
piden
35 minutes ago, bubblyk said:

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

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

 

36 minutes ago, bubblyk said:

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

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

 

37 minutes ago, bubblyk said:

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

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

Link to post
Share on other sites
bubblyk
10 минут назад, piden сказал:

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

 

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

 

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

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

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

Link to post
Share on other sites
soklakov
25 минут назад, bubblyk сказал:

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

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

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

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

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

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

некомфорная

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

 

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

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

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

  • Нравится 2
Link to post
Share on other sites
bubblyk

В любом случае спасибо за ответ, навели на верный путь!

 

Как выполню, отчитаюсь :smile:

Link to post
Share on other sites
karachun
21 минуту назад, soklakov сказал:

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

*некомфортная

  • Нравится 2
  • Хаха 1
Link to post
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.
Note: Your post will require moderator approval before it will be visible.

Guest
Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

  • Recently Browsing   0 members

    No registered users viewing this page.



×
×
  • Create New...