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

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


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

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

 

Решаю задачу в 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 пользователей

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




  • Сообщения

    • Bot
      Команда CSD пообщалась с представителями компании «КС-ПРО» и выяснила, как происходило внедрение G-Station, сколько времени занял процесс перехода на новую платформу, узнала об организации совместной работы и о функционале платформы. Основное направление компании «КС-ПРО» — оказание комплекса услуг технического заказчика, в том числе и для крупных офисных девелоперов. Внедрение G-Station в «КС-ПРО» проходило в ноябре 2022 года. G-Station — это всего лишь один из модулей комплексной платформы G-Tech Suite. Специалисты компании отмечают, что G-Station действительно стала хорошей альтернативой BIM 360. Специалисты «КС-ПРО» получили доступ к полноценной единой информационной среде с облачным хранилищем файлов и моделей, а также такими возможностями, как: Управление проектами, пользователями и подрядчиками; Создание чек-листов; Создание шаблонов чек-листов с процессами согласования; Передача документации на строительную площадку; Просмотр 3D- и 2D моделей [...] View the full article
    • Bot
      АСКОН, российский разработчик инженерного программного обеспечения и ИТ-интегратор, подвел итоги 2023 года. Выручка компании выросла на 47% и достигла 5,6 млрд рублей, штат сотрудников увеличился до 1250 человек. Клиентская база пополнилась тысячью предприятий, и сегодня с решениями АСКОН работают почти 16 000 заказчиков из всех отраслей промышленности и строительства. На динамику бизнеса положительно повлияли рост промышленного производства, сопровождающийся инвестициями в цифровизацию, крупные инфраструктурные проекты и курс страны на технологический суверенитет. Лидерские позиции компании в разработке и интеграции подтверждены профессиональными рейтингами. АСКОН, единственный из разработчиков инженерного ПО, вошел в ТОП-100 российских ИТ-компаний в рейтинге CNews; второй год подряд — в десятку крупнейших ИТ-поставщиков в сфере промышленности и строительства по данным TAdviser; впервые назван в тройке ведущих proptech-компаний как [...] View the full article
    • esergey
      это всего лишь видео - я не шарю в этом ...  
    • lem_on
      ну с дуру известно что сломать можно.
    • Viktor2004
      руку привязки так сломать легко
    • lem_on
      По моему вполне логично если станок вывалится в ошибку если рука не доехала до места. У меня так же если кулачки или деталь на пути, просто пихаеш ее до места и станок опять активен. Но нынешние пановья даже не могут написать модель станка.
    • Viktor2004
      Я согласен что скорее всего проблема механическая Но если логика прописана криво и возможно не предусмотрела остановку в промежуточном состоянии, разве не логично будет попробовать принудительно подав напряжение дернуть эту руку вверх-вниз? Возможно то что туда попало выпадет  
    • Guhl
      Если оставить за скобками вопрото том, что до м19 работает нормально, а после нет, то вы не считали сколько у него реально импульсов на оборот? с помощью стороннего плк, например  А если ориентацию м5 снимать, а не м20?
    • lem_on
      Что это за станок такой в котором сразу ладер ковырять надо, даже не смотря на возможность механической проблемы? Или профдеформация?
    • Viktor2004
      не сразу я понял в чем вопрос. Долго соображал что такое режим управления скоростью. При завершении ориентации PMC снимает сигнал G70.6 ? И если он после снятия сигнала продолжает удерживать шпиндель, при каких условиях эта ориентация все же снимается? После нажатия аварийного грибка или еще как?
×
×
  • Создать...