Коэффициент теплоотдачи - Workbench и APDL

[Axer 3]

Здравствуйте, уважаемые дамы и господа.

Прикладываю ко всем поверхностям разогретого параллелепипеда (сталь по умолчанию) конвекцию таблицей, получаю изменение температуры с площадкой (рис. 1) именно при той температуре, когда коэффициент теплоотдачи «ныряет», на 250ти г радусах (рис.2).

Выглядит логично. Если же численно эту же конвекцию задать через командную вставку, то эта площадка «вдруг» совсем не там, а примерно на 500ах градусах (рис. 3).

 

alls

*DIM,alfaair,table,15,1,1,TEMP 

alfaair(0,1)=0.0

 

alfaair(1,0)=225 
alfaair(2,0)=250
alfaair(3,0)=300 
alfaair(4,0)=350 
alfaair(5,0)=400
alfaair(6,0)=450
alfaair(7,0)=500 
alfaair(8,0)=550
alfaair(9,0)=600
alfaair(10,0)=650 
alfaair(11,0)=700
alfaair(12,0)=750
alfaair(13,0)=800
alfaair(14,0)=850
alfaair(15,0)=876

alfaair(1,1)=50
alfaair(2,1)=5
alfaair(3,1)=59
alfaair(4,1)=63 
alfaair(5,1)=68
alfaair(6,1)=70
alfaair(7,1)=85 
alfaair(8,1)=90
alfaair(9,1)=107
alfaair(10,1)=128 
alfaair(11,1)=166
alfaair(12,1)=195
alfaair(13,1)=208 
alfaair(14,1)=230
alfaair(15,1)=250

nsubst,1,1,1
*do,i,1,arg2
 time,arg1*i

 NSEL,s,ext
 esel,all
 SF,ALL,CONV,%alfaair%,22
 alls

SOLVE

*enddo
fini
*abbr,solve,*abbr,solve

 

С чем это может быть связано? Как в командной вставке правильно задать «такой же» коэффициент теплоотдачи, как в интерфейсе Workbench?

Большое спасибо!

 

test1.wbpz

[piden 2 874]

13 minutes ago, Axer said:

Как в командной вставке правильно задать «такой же» коэффициент теплоотдачи, как в интерфейсе Workbench?

Файл не смотрел, но.. ВБ ведь генерит .dat файл, который и скармливается решателю. Почему бы не открыть его и не посмотреть, как там это дело задается?

13 minutes ago, Axer said:

а примерно на 500ах градусах

Т.е. примерно значения ВБ + 273.15 , да?

Или что-то по-ошибке два раза приложено...

Изменено 8 июня 2021 пользователем piden

[AlexArt 277]

1 минуту назад, piden сказал:

Файл не смотрел, но.. ВБ ведь генерит .dat файл, который и скармливается решателю. Почему бы не открыть его и не посмотреть, как там это дело задается?

Как же вы безтактны! Люди думают, что ВБ напрямую в фортране работает, а вы им все рушите 

[piden 2 874]

8 minutes ago, AlexArt said:

Как же вы безтактны!

Spoiler

 

 

@piden до того бестактный, что будет исправлять ошибку в слове "безтактный" 

 

 

[Graf Kim 217]

26 минут назад, Axer сказал:

С чем это может быть связано? Как в командной вставке правильно задать «такой же» коэффициент теплоотдачи, как в интерфейсе Workbench?

 

Как минимум половина проблем при написании вставок связана с размерностями.

Попробуйте переключиться на метры.

[piden 2 874]

52 minutes ago, Axer said:

*do,i,1,arg2
 time,arg1*i

 NSEL,s,ext
 esel,all
 SF,ALL,CONV,%alfaair%,22
 alls

SOLVE
*enddo

Вы конвекцию на каждом шаге новую прикладываете, а старую не убиваете. Так задумано?

[Axer 3]

9 минут назад, piden сказал:

Вы конвекцию на каждом шаге новую прикладываете, а старую не убиваете. Так задумано?

Спасибо. Уже убиваю, не помогает :( 

33 минуты назад, Graf Kim сказал:

 

Как минимум половина проблем при написании вставок связана с размерностями.

Попробуйте переключиться на метры.

Спасибо. Я же и так на метрах, по крайней мере в настройках проекта в WB и непосредственно в меню Mechanical.

[piden 2 874]

1 hour ago, Axer said:

Если же численно эту же конвекцию задать через командную вставку, то эта площадка «вдруг» совсем не там, а примерно на 500ах градусах (рис. 3).

Если в интерфейсе для конвекции выбрать опцию Coefficient type = Average Film Temperature, то результаты совпадут с командной вставкой.

Во вставке вы прикладываете нагрузку напрямую к узлам. А ВБ для конвекции создает surf152, и выставляет им кейопции сообразно вашему заданному Coefficient type.

Как задавать эти кейопции в случае прямого задания конвекции на узлы - хз. Кажется, проще убить нагрузки на surf152, созданные ВБ, и задать на эти элементы свои, уже командной вставкой.

[Axer 3]

4 минуты назад, piden сказал:

Если в интерфейсе для конвекции выбрать опцию Coefficient type = Average Film Temperature, то результаты совпадут с командной вставкой.

Во вставке вы прикладываете нагрузку напрямую к узлам. А ВБ для конвекции создает surf152, и выставляет им кейопции сообразно вашему заданному Coefficient type.

Как задавать эти кейопции в случае прямого задания конвекции на узлы - хз. Кажется, проще убить нагрузки на surf152, созданные ВБ, и задать на эти элементы свои, уже командной вставкой.

Спасибо большое. А можно ли во вставке создавать эти же 152 и уже к ним прикладывать конвекцию? 

[piden 2 874]

12 minutes ago, Axer said:

А можно ли во вставке создавать эти же 152 и уже к ним прикладывать конвекцию? 

Да.

Создайте для внешних поверхностей named selection с именем ext_surface

 

ARG1 = 1

ARG2 = 200

ARG3 = 222

 

Spoiler

!!! MKS units !!!

/prep7                          ! Enter into /prep7 as ET and ESURF and KEYOPT are valid only in /PREP7

et,arg3,152                     ! Create a new element type for SURF152                
cmsel,s,ext_surface             ! Select nodes on a face via Face named selection
type,arg3                       ! Set the element type as ARG3
esurf                           ! Create surface effect elements on the selected nodes
allsel,all                      ! Select everything back  


keyopt,arg3,8,3                 ! 3 - for surface temperature
keyopt,arg3,13,0

/solu                           ! Re-enter /SOLU


!! convection coefficient
*DIM,alfatable,TABLE,15,1,1,TEMP,
! Temperature values
alfatable(1,0,1) = 225.
alfatable(2,0,1) = 250.
alfatable(3,0,1) = 300.
alfatable(4,0,1) = 350.
alfatable(5,0,1) = 400.
alfatable(6,0,1) = 450.
alfatable(7,0,1) = 500.
alfatable(8,0,1) = 550.
alfatable(9,0,1) = 600.
alfatable(10,0,1) = 650.
alfatable(11,0,1) = 700.
alfatable(12,0,1) = 750.
alfatable(13,0,1) = 800.
alfatable(14,0,1) = 850.
alfatable(15,0,1) = 876.
! Load values
alfatable(1,1,1) = 50.
alfatable(2,1,1) = 5.
alfatable(3,1,1) = 59.
alfatable(4,1,1) = 63.
alfatable(5,1,1) = 68.
alfatable(6,1,1) = 70.
alfatable(7,1,1) = 85.
alfatable(8,1,1) = 90.
alfatable(9,1,1) = 107.
alfatable(10,1,1) = 128.
alfatable(11,1,1) = 166.
alfatable(12,1,1) = 195.
alfatable(13,1,1) = 208.
alfatable(14,1,1) = 230.
alfatable(15,1,1) = 250.


!! apply custom convection load
esel,s,type,,arg3
nsle
sf,all,conv,%alfatable%,22.
ALLSEL

!! solve
nsubst,1,1,1
*do,i,1,arg2
time,arg1*i
SOLVE
*enddo

 

 

Изменено 8 июня 2021 пользователем piden

[Axer 3]

Огромное спасибо, Piden! Я named selection в интерфейсе WB сделал, а командами прямо в снипете тоже можно?

[piden 2 874]

7 minutes ago, Axer said:

а командами прямо в снипете тоже можно?

Почему бы нет?

[Axer 3]

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

Почему бы нет?

Спасибо огромное. Сделал.

cm,s,ext
cm,ext_surface,node
Еще ламерский вопрос: почему SF в цикле так же работает, как и просто пред ним, и почему не нужно SFDELE ?

[piden 2 874]

31 minutes ago, Axer said:

Еще ламерский вопрос: почему SF в цикле так же работает, как и просто пред ним, и почему не нужно SFDELE ?

Похоже, чтобы нагрузки от SF суммировались, нужно изменить режим суммирования нагрузок с помощью SFCUM. По-умолчанию там параметр repl, т.е. новая нагрузка SF перезаписывает значение старой. Вот такая poka-yoke фича.

Изменено 8 июня 2021 пользователем piden

[Axer 3]

А раньше не так везло... 

 

[piden 2 874]

3 minutes ago, Axer said:

А раньше не так везло... 

Так там ведь вы новые узлы выбирали, а не те же самые.

[Axer 3]

Да, спасибо большое, ведь при одной и той же выборке конвекция так же действует. Дальше буду пробовать движение конвекции, как и тут: 

Если осесимметричную задачу делать, то 151, верно?

 

[Axer 3]

Снова здравствуйте, уважаемые дамы и господа.

В 3D все прекрасно получилось, и движение конвекции по поверхности тела.

А вот при попытке повторить в осесимметричной постановке сразу такая проблема.

Прикладываю в WB к одной поверхности конвекцию таблицей (рис.1), получаю изменение температуры (рис. 2). Если эту же конвекцию пишу в командной вставке, охлаждение идет гораздо быстрее (рис. 3). Подскажите, пожалуйста, что теперь не так :(

test3.wbpz

[Axer 3]

Извините, вот код вставки:

 

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

 

/prep7                          

et,222,151                                 
cmsel,s,ext_surface             ! Face named selection created in WB
type,222                    
esurf                          
allsel,all                   

keyopt,222,8,3                 
keyopt,222,13,0
alls
/solu                           

!! convection coefficient
*DIM,alfatable,TABLE,15,1,1,TEMP,
! Temperature values
alfatable(1,0,1)=225
alfatable(2,0,1)=250
alfatable(3,0,1)=300
alfatable(4,0,1)=350
alfatable(5,0,1)=400
alfatable(6,0,1)=450
alfatable(7,0,1)=500
alfatable(8,0,1)=550
alfatable(9,0,1)=600
alfatable(10,0,1)=650
alfatable(11,0,1)=700
alfatable(12,0,1)=750
alfatable(13,0,1)=800
alfatable(14,0,1)=850
alfatable(15,0,1)=876
! Load values
alfatable(1,1,1)=50
alfatable(2,1,1)=5
alfatable(3,1,1)=59
alfatable(4,1,1)=63
alfatable(5,1,1)=68
alfatable(6,1,1)=70
alfatable(7,1,1)=85
alfatable(8,1,1)=90
alfatable(9,1,1)=107
alfatable(10,1,1)=128
alfatable(11,1,1)=166
alfatable(12,1,1)=195
alfatable(13,1,1)=208
alfatable(14,1,1)=230
alfatable(15,1,1)=250

esel,s,type,,222
nsle
sf,all,conv,%alfatable%,22
ALLSEL

nsubst,1
*do,i,1,arg2
time,arg1*i
SOLVE
*enddo

 

 

[Axer 3]

Разобрался!

Порылся в файле ds, для SURF151 нужно еще keyopt,222,3,10

Всем огромное спасибо!