[necto]

Здравствуйте!

Есть треугольная пластинка. Есть начальная температура, затем температура прикладывается к двум сторонам, а к третьей dT/dn=40.
Нужно посмотреть, как распределиться температура по пластинке, например, через 30 сек.

Три вопроса:
1. Как правильно приложить градиент температуры dT/dn=40? Возможно ли это в Mechanical?
Я прикладываю как HEAT FLOW с учетом того, что тепловой поток q=-kxx*gradT, и можно приложить к каждому узлу на стороне -kxx/c, где kxx - задаваемый в свойствах материала коэффициент теплопроводности.

2. Правильно ли я прикладываю время нагрева?
time,1
тепловые гу
solve
time,30
solve


В Contour plot я вижу одну и ту же картину, что без time,30 , что с.

3. Как правильно выбрать тип элемента? PLANE 75 и PLANE 55 дают совершенно разный характер распределения температур.

На всякий случай код:

Код:

[Выделить все]

/PREP7  
ANTYPE,0
ET,1,PLANE75
nstep=10

!Свойства материала

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDE,KXX,1  
MPDATA,KXX,1,,1 
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,10   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,C,1,,100 

!Геометрия

K,1,0,0,0,  
K,2,8,0,0,  
K,3,4,4,0, 
 
L,1,3  
L,3,2  
L,2,1 
 
A,1,2,3


                !создаем элементы
LSEL,S, , ,1
LSEL,S, , ,2
LSEL,S, , ,3
LESIZE,ALL,,,50,1,,,,1

ALLSEL,ALL

MSHAPE,0,2D
MSHKEY,1
AMESH,ALL
FINISH

/SOL
TUNIF,0,

time,1
nsubst,nstep,nstep,nstep
DL,3, ,TEMP,50,0 
DL,2, ,TEMP,100,0

!задать dT/dn
lsel,s,,,1
nsll,s,all
F,all,HEAT,-40/1

ALLSEL,ALL


SOLVE

time,30
nsubst,nstep,nstep,nstep

SOLVE   
FINISH  
/POST1

[beryl]

Я не специалист по тепловым расчетам.
Как я понимаю, Вы делаете статический расчет.

В статическом расчете время - условная величина. Вы видите в результатах как бы картину установившуюся, через бесконечный промежуток времени.
Не важно, time 1, или time 1000 в данном случае.
Это как приложили силу, получили прогиб от этой силы в балке по формуле. Время в формулу не входит

Если вам нужна динамика, видимо надо как-то подругому сделать..

[Chardash]

не специалист по тепловым расчетам.

аналогично
если не нужен НДС, можно использовать сразу Fluent

Сообщение от necto PLANE 75 и PLANE 55 дают совершенно разный характер распределения температур.

возможно вместо PLANE 75, гармонического анализа использовать восьмиузловой PLANE77.
скорость теплового потока

HEAT FLOW

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

q=-kxx*gradT

, но тогда задавайте как Heat Flux.

http://cae-club.ru/forum/gidrodinami...thermal?page=5 посмотрите еще тут или на САПРе2000 спросите

[necto]

Как я понимаю, Вы делаете статический расчет.

Да, спасибо, тут и крылась ошибка с нерабочим временем.
Тип анализа был заменен на переходный
ANTYPE,4
Время заработало!

Сообщение от Chardash но тогда задавайте как Heat Flux

Ага, спасибо )

В процессе поиска ответа на вопросы времени я нашла несколько примеров в VM (113-115), где решаются немного похожие задачи.
В связи с чем появилось еще больше вопросов. Проясните, пожалуйста, такие моменты:

VM115.
Бесконечная пластина первоначальной температуры T0 подвергается внезапной генерации тепла и температуры Ts. Определить распределение тепла в пластине после 12 минут.
Вот код решателя:

Код:

[Выделить все]

AUTOTS,ON
OUTPR,,LAST
TIME,(12/60)       ! CONVERT 12 MIN TO HRS
DELTIM,0.01
BFUNIF,TEMP,60     ! INITIAL UNIFORM TEMPERATURE
KBC,1              ! APPLY STEP LOADS
DK,1,TEMP,32       ! SURFACE TEMPERATURE APPLIED TO KEYPOINT 1
1
BFK,ALL,HGEN,4E4   ! HEAT GENERATION SPECIFIED AT KEYPOINTS
1 & 2

Поясните, пожалуйста, по этому примеру:
1 .команда
AUTOTS,ON !автоматическая активизация времени
- не очень понятно, в чем смысл этой команды?

2. KBC,1 !пошаговая нагрузка -зачем, если сказано, что нагрузка внезапна и можно ли в принципе ее заменять, например, NSUBST? Т е являются ли эти команды аналогом друг друга?

3. В чем принципиальное отличие задания равномерной начальной температуры на узлы между
TUNIF ,0 и BFNIF,TEMP,60 ?

4. Почему в данном примере сначала указывается время TIME,(12/60) , и только после этого прикладывается нагрузка?
Мне казалось логичным написать так
-момент времени 1
-приложили тепловую нагрузку
-посчитали
-момент времени 2
-посчитали
-см результат

5. команды 1 и 1 & 2 объединяют нагрузки для совместного учета? Зачем тогда просто 1 ?

6. И последний вопрос рубрики:
Какую команду стоит применить, чтобы разбить на шаги этапы решения (не приложения нагрузки), а так, что нагрузка уже приложена в момент time,1, мне нужны результаты в момент time,2 и Ansys делит этот временной промежуток на несколько интервалов и на каждом просчитывает, а я смогу в дальнейшем через Read Results->By Pick посмотреть динамику изменения температуры тема во времени?

Спасибо всем, кто хоть капельку прояснит)

[beryl]

Как я понимаю:
2) KBC - как раз указывает, прикладываете ли вы нагрузку сразу вначале, или плавно, пропорционально времени, достигая полного значения на последнем шаге.

В динамике transient по умолчанию нагрузка прикладывается сразу, на первом шаге.
А дальше идет разбиение по времени - т.е. вы указываете шаги по времени, а не по нагрузке.
Это возможно ответ и на вопрос 6).

4) Порядок не важен, пока не нажали solve

[necto]

Спасибо за пояснение про нагрузку.
Только не ясно, зачем тогда использовать KBC,1 , если в transient нагрузка прикладывается на первом шаге?

4) Порядок не важен, пока не нажали solve

Да, это понятно, просто немного различается идеология учета нагрузки во времени в примере и в моем понимании. Кажется, надо корректировать понимание)
В одной и той же постановке - приложили нагрузку, посмотрели характер ее воздействия через определенное время-
в примере выглядит так:
-сразу задается время через которое надо глядеть результат
-прикладывается нагрузка
-solve

Здесь мне не ясно, ансис сам догадывается, что надо рассчитать результат по истечении заданного времени при заданных нагрузках?
Мне казалось более логичным использовать два solve, приложили, посчитали, истекло время, посчитали. Например, при проведении эксперимента в лаборатории порядок был бы именно такой, только вместо посчитали - измерили.

[beryl]

Я как-то делал transient, хотел посмотреть, как будет разрушаться купол при постепенном нагружении его избыточной нагрузкой. Вот там как-раз я kbc менял для постепенного приложения конечной нагрузки.

Если в начальный момент времени не было никаких нагрузок, то можно сразу задать time,30, и разбить помельче на шаги. Поставить, чтобы промежуточные шаги тоже сохранялись в результатах. А потом можете смотреть результаты в промежуточные моменты времени. (в частности для time,1)

Если же в начале нужно приложить нагрузку, например от собственного веса (или другой характер распределения температур), тогда можно сначала сделать один solve только с собственным весом, а вторым шагом добавить туда температуру.
И/или в первом шаге Вы можете например сделать разбивку на подшаги грубую, а второй шаг разбить подробнее.

(Кстати, в transient важно не перепутать с размерностями. Если время вы задаете в секундах, то остальные величины в СИ, т.е. модель строить в Метрах, и т.д.)

[necto]

Спасибо за ответы!

Сообщение от beryl Я как-то делал transient, хотел посмотреть, как будет разрушаться купол при постепенном нагружении его избыточной нагрузкой. Вот там как-раз я kbc менял для постепенного приложения конечной нагрузки.

Тогда как отличить (указать) шаги по приложению нагрузки (на последнем шаге нагрузка достигает полного значения) от временных шагов, которые стартуют от этапа приложения полной нагрузки до интересующего времени?

Сообщение от beryl Если в начальный момент времени не было никаких нагрузок, то можно сразу задать time,30, и разбить помельче на шаги. Поставить, чтобы промежуточные шаги тоже сохранялись в результатах. А потом можете смотреть результаты в промежуточные моменты времени. (в частности для time,1)

т е свой код я могу переписать таким образом: ?

Код:

[Выделить все]

/SOL 
TUNIF,0,   !температура в начале
time,(30/3600) 
kbc,1                                !например, вкл опцию учета шагов нагружения
nsubst,10                         ! шаги по времени нагружения или разбивка интервала времени после нагружения?
 
DL,3, ,TEMP,50,0              !на линию 3 температуру 50
DL,2, ,TEMP,100,0            !на линию 2 темп.100
SFL,1,HFLUX,-40,-40        !на линию 1 dT/dn=40
ALLSEL ALL
SOLVE

[beryl]

В transient вы указываете только шаги по времени (например командой nsubst).

Если KBC,0 - то нагрузка пропорционально разбиению по времени бьется так же. И в конечный момент времени, который вы задали time, достигает максимума.
Т.е. если вам надо смоделировать и сам процесс нагрева ребер пластины - например за 10 секунд ребра нагреваются до 100. То можно поставить time,10, и температуру 100. Следующим шагом вы можете поставить time,1000 и просчитать, как этот нагрев будет распространяться дальше в самой пластинке, когда ребра уже нагреты, и температура 100 градусов на ребрах поддерживается. (возможно в следующем шаге нужно и kbc,1 - переключить)


Если KBC,1 - то нагрузка приложена сразу, но сдерживается силами инерции.


Как я понимаю, в Вашем случае - нагрузка на края пластины прикладывается сразу, нагрев ребер пластинки 50, и 100 градусов, а дальше Вы хотите посмотреть, как этот нагрев будет распространяться во времени. Т.е. шагов по нагрузке Вам не нужно. и можете смело оставить KBC,1.
В этом случае nsubst задает исключительно шаг по времени. Можно сделать еще мельче.

time (30/3600) Это Вы что хотите получить? Если часы перевести в секунды, то умножить надо

[necto]

beryl, спасибо, внесли ясность насчет шагов!!!

Только вот еще вопрос, (на рисунке)
если я учитываю kbc,1 и nsubst,10 относятся к времени, почему в Read Results->By Pick нет подпунктов по времени, а сразу идет 30? 10 подшагов означают, что он разбил временной интервал на 10 частей, в каждом просчитал, но как сделать так, что я , например, задаю 10 шагов и в столбце Time вижу числа 3
6
9 и т.д. ?


Просто очень хочется потыкать в разное время и посмотреть динамику нагрева/охлаждения!

Сообщение от beryl time (30/3600) Это Вы что хотите получить? Если часы перевести в секунды, то умножить надо

Да, я не учла Ваше сообщение о секундах и системе СИ, и как в примере VM115 прикладывала время, пересчитанное на часы: (

[Chardash]

OUTPR - вывод данных и OUTRES - запись данных , опции контроля выходной информации

[beryl]

Да, команда OUTRES,ALL,ALL - запись всех промежуточных результатов.)