Все ответы об инженерном анализе

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

Тепловой поток/плотность теплового потока

Dmitry Krasilnikov
Аватар пользователя Dmitry Krasilnikov

Добрый день. Нужна помощь с тепловым потоком. Рассматривается следующая задача: имеется шар, на который с одной стороны (например с верхней) действует тепловой поток, меняющийся во времени. В ansys mechanical решаю данную задачу в плоском виде: строю круг, прикладываю тепловой поток. На этом этапе возникает две трудности:

1) Считаю через heat flow. При расчете возникает ошибка, что температура в одном из узлов больше миллиона (превышает установленный в ansys предел). Соответственно, расчет не выполняется.

2) Известна плотность теплового потока. Пробую считать через heat flux. Выскакивает ошибка: "the nlist field on the SF command is invalid". 

Фрагмент кода, с которым возникают трудности, прикладываю:

/SOL
ANTYPE,4
TRNOPT,FULL
LUMPM,0
NSUBST,10,100,1
KBC,1
!TIME,13.574

N1=1000 ! число шагов
*do,i,1,1000
t1=i*13.574/N1 ! время в секундах
v1=skor(t1)/1000 ! скорость в км/с, берется из таблицы, подгруженной заранее
ro1=plotn(t1) ! плотность, берется из таблицы, подгруженной заранее

LSEL,S, , , 1
NSLL,S,1
*get,j,node,,num,min
*dowile,j
beta1=asin(nx(j)/sqrt(nx(j)*nx(j)+ny(j)j)) ! Коэффициенты для расчета теплового потока
rr1=((ro1/1000)**(1.32-0.01223*v1-1))*(v1**(51.89-1.558*v1+0.02659*v1*v1-3)) ! Коэффициенты для расчета теплового потока
rk1=450**(1.688-0.1796*v1+0.005381*v1*v1) ! Коэффициенты для расчета теплового потока
Chr1=1.9576e-38*rr1*rk1 ! Коэффициенты для расчета теплового потока
nn1=1/(0.051*v1-0.43)+1.811 ! Коэффициенты для расчета теплового потока
hf1=((3.3e-5*(ro1/8)**0.5)*((v1*1000)**3.2*(0.55+0.45*cos(2*beta1)))) ! Формулы для расчета теплового потока

hr1=0.5*ro1*v1**3*Chr1*cos(beta1)**nn1 ! Формулы для расчета теплового потока

hs1=-(hf1+hr1) ! Формулы для расчета теплового потока

! f,j,heat,hs1 ! использую эту команду, чтобы прикладывать тепловой поток (возникает проблема №1)
SFE,j,HFLUX,hs1 ! использую эту команду, чтобы прикладывать плотность теплового потока (возникает проблема №2)
j=ndnext(j)
*enddo
allsel,all
time,t1
solve
*enddo

Заранее спасибо.

Ruslan Abdrahimov
Аватар пользователя Ruslan Abdrahimov

Добрый день!

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

По второй ошибке, тоже вопрос. Вы пишите ошибку про SF команду, которая может возникнуть, если приложить ее на 1 один узел, а должно быть указано несколько узлов (должна быть некоторая площадь для расчета теплового потока). А в коде указана команда SFE, для приложения потока на элементы и она записывается как SFE,j,LKEY,HFLUX,,hs1, соответственно тут нужно правильно выбрать элемент j (у вас в коде j это номер узла) и правильно указать сторону элемента (LKEY).

Dmitry Krasilnikov
Аватар пользователя Dmitry Krasilnikov

По первой ошибке - тепловой поток это уже и есть граничное условие.. по второй совет помог, спасибо

Добавить комментарий

Войдите или зарегистрируйтесь, чтобы отправлять комментарии
Приложить файл

Максимальный размер файла: 999 МБ.
Допустимые типы файлов: txt doc docx xls xlsx pdf rar zip 7zip tar.