Почему у меня не работает IC,,TEMP ?

[4estor 1]

Здравствуйте.
 
Делаю термомеханический расчет в 2D осесимметрике.
Нужно задать начальную температуру в середине детали.
Столкнулся с проблемой: задание нач.темп. с помощью команды "IC" либо не получается, либо не видно этой нач.темп-ры на постпроцессоре. Почему так? Или как можно по-другому задать начальную температуру, если так нельзя?
 
Мои параметры решателя не противоречат использованию IC (по хэлпу):
ANTYPE,4
NLGEOM,1
KBC,1
TIMINT,ON
TRNOPT,FULL
SOLCONTROL,ON

 

Err-файл пустой.

 

Полный макрос прикрепил как файл или смотрите в "Скрытом тексте".

IC в конце макроса.

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

 

FINISH
/clear,start
/PREP7

! исходные данные (начало)

l=200 ! длина ТП
d=800 ! диаметр ТП
t=40 ! толщина ТП
angdeg=40 ! угол в градусах
xv=5 ! длина шва внутри
nv=5 ! ширина внутренней наплавки
nn=4 ! ширина наружной наплавки

! исходные данные (конец)

ll=l/2
r=d/2
rn=(d/2)+t
ang=0.5*angdeg*3.14159/180 ! 20 градусов в радианах
xn=t*tan(ang)+(xv/2) ! длина (полная) шва снаружи

k,1,r,ll
k,2,rn,ll
k,3,rn,xn
k,4,r,xv/2

KSYMM,Y,1,4,1,,1,0

k,9,r,(xv/2)+nv ! внутренняя наплавка
k,10,rn,xn+nv

k,11,r,(xv/2)+nv+nn ! внешняя наплавка
k,12,rn,xn+nv+nn

a,1,2,12,11
a,5,6,7,8
a,3,4,8,7
a,3,4,9,10
a,9,10,12,11


! Задание свойств КЭ (начало)

ET,1,PLANE223
KEYOPT,1,1,11
KEYOPT,1,3,1

! Задание свойств КЭ (конец)


! Задание материалов (начало)


nmats=2
ntemps=6
ntemps14=14
puas=0.3
to_air=0.014 ! теплосъем в окр.среду
to_steel=0.05
t_air=20 ! температура окр.среды

*SET,t
*DIM,t, ARRAY,nmats,ntemps

*SET,r02
*DIM,r02, ARRAY,nmats,ntemps

*SET,E
*DIM,E, ARRAY,nmats,ntemps14

*SET,al
*DIM,al, ARRAY,nmats,ntemps14

*SET,tm
*DIM,tm, ARRAY,nmats,ntemps14

*SET,c
*DIM,c, ARRAY,nmats

*SET,kxx
*DIM,kxx, ARRAY,nmats

*SET,enth
*DIM,enth, ARRAY,nmats

*SET,dens
*DIM,dens, ARRAY,nmats

*SET,hf
*DIM,hf, ARRAY,nmats

! материал 1 - перлит 10ГН2МФА - только верх

t(1,1)=20
t(1,2)=50
t(1,3)=150
t(1,4)=250
t(1,5)=350
t(1,6)=1500

r02(1,1)=343
r02(1,2)=333
r02(1,3)=314
r02(1,4)=304
r02(1,5)=294
r02(1,6)=1

E(1,1)=210
E(1,2)=207
E(1,3)=205
E(1,4)=202
E(1,5)=200
E(1,6)=197
E(1,7)=195
E(1,8)=190
E(1,9)=185
E(1,10)=180
E(1,11)=175
E(1,12)=170
E(1,13)=165
E(1,14)=1

al(1,1)=11
al(1,2)=11.5
al(1,3)=11.9
al(1,4)=12.2
al(1,5)=12.5
al(1,6)=12.8
al(1,7)=13.1
al(1,8)=13.4
al(1,9)=13.6
al(1,10)=13.8
al(1,11)=14
al(1,12)=14.2
al(1,13)=14.4
al(1,14)=20.3

tm(1,1)=5356
tm(1,2)=5141
tm(1,3)=4819
tm(1,4)=4793
tm(1,5)=4914
tm(1,6)=100

c(1)=500
kxx(1)=47e-3
hf(1)=50
enth(1)=1e6
dens(1)=7860*1e-9

! материал 2 - аустенит 12Х18Н10Т - электрод, наплавки, низ

t(2,1)=20
t(2,2)=100
t(2,3)=200
t(2,4)=350
t(2,5)=600
t(2,6)=1500

r02(2,1)=196
r02(2,2)=189
r02(2,3)=181
r02(2,4)=167
r02(2,5)=147
r02(2,6)=1

E(2,1)=205
E(2,2)=202
E(2,3)=200
E(2,4)=195
E(2,5)=190
E(2,6)=185
E(2,7)=180
E(2,8)=175
E(2,9)=170
E(2,10)=167
E(2,11)=165
E(2,12)=162
E(2,13)=160
E(2,14)=1

al(2,1)=16.2
al(2,2)=16.4
al(2,3)=16.6
al(2,4)=16.8
al(2,5)=17
al(2,6)=17.2
al(2,7)=17.4
al(2,8)=17.6
al(2,9)=17.8
al(2,10)=18
al(2,11)=18.2
al(2,12)=18.4
al(2,13)=18.5
al(2,14)=20.3

tm(2,1)=1638
tm(2,2)=1609
tm(2,3)=1658
tm(2,4)=1685
tm(2,5)=1348
tm(2,6)=100

c(2)=500
kxx(2)=47e-3
hf(2)=50
enth(2)=1e6
dens(2)=7860*1e-9


*DO,i,1,nmats
TB,BISO,i,ntemps
*DO,j,1,ntemps
TBTEMP,t(i,j)
TBDATA,,r02(i,j),tm(i,j)
*ENDDO
*ENDDO


MPTEMP,,20,50,100,150,200,250
MPTEMP,7,300,350,400,450,500,550
MPTEMP,13,600,1500

*DO,i,1,nmats

MPDATA,PRXY,i,,puas
MPDATA,C,i,,c(i)
MPDATA,KXX,i,,kxx(i)
MPDATA,ENTH,i,,enth(i)
MPDATA,DENS,i,,dens(i)
MPDATA,HF,i,,hf(i)

*DO,j,1,ntemps14
MPDATA,EX,i,,E(i,j)*1000
MPDATA,ALPX,i,,al(i,j)*1e-6
*ENDDO
*ENDDO




! Задание материалов (конец)



! Создание сетки (начало)

ndivrad=20 ! число разбиений в радиальном направлении
ndivax=4 ! число разбиений шва в осевом направлении

LESIZE,1,,,ndivrad
LESIZE,3,,,ndivrad
LESIZE,13,,,ndivrad
LESIZE,9,,,ndivrad
LESIZE,7,,,ndivrad
LESIZE,5,,,ndivrad

LESIZE,10,,,ndivax
LESIZE,11,,,ndivax

LESIZE,16,,,2
LESIZE,15,,,2
LESIZE,12,,,3
LESIZE,14,,,3

LESIZE,4,,,20,3
LESIZE,8,,,20,3
LESIZE,2,,,20,1/3
LESIZE,6,,,20,1/3

/NUMBER,1
/PNUM,MAT,1

MAT,1
AMESH,1

MAT,2
AMESH,2

MAT,2
AMESH,3

MAT,2
AMESH,4

MAT,2
AMESH,5

! Создание сетки (конец)


! Нагружение и решение (начало)

/SOLU
ANTYPE,4
NLGEOM,1
KBC,1
TIMINT,ON
TRNOPT,FULL
SOLCONTROL,ON

NSUBST,5,1e3,5
OUTRES,ALL,ALL

!DL,1,,UX
!DL,1,,UY
!DL,5,,UX
DL,5,,UY

LSEL,S,,,1
LSEL,A,,,2
LSEL,A,,,15
LSEL,A,,,14
LSEL,A,,,11
LSEL,A,,,6
LSEL,A,,,5
LSEL,A,,,8
LSEL,A,,,10
LSEL,A,,,12
LSEL,A,,,16
LSEL,A,,,4
SFL,ALL,CONV,to_air,,t_air
ALLSEL

TUNIF,t_air






ASEL,S,,,3
ALLSEL,BELOW,AREA
IC,ALL,TEMP,1000 !!!! =========== начальная температура, которая не почему-то задается / или задается, но её не видно ============
ALLSEL




ASEL,S,,,3
ALLSEL,BELOW,AREA
ALLSEL

TIME,1

SOLVE

/POST1

PLNSOL,TEMP,, 0

macrott.txt

Изменено 7 сентября 2013 пользователем 4estor

[soklakov 1 475]

ASEL,S,,,3 ALLSEL,BELOW,AREA IC,ALL,TEMP,1000 !!!! =========== начальная температура, которая не почему-то задается / или задается, но её не видно ============ ALLSEL

Если для проверки задать температуру ну всю модель, тоже нет работает?

[4estor 1]

 

ASEL,S,,,3 ALLSEL,BELOW,AREA IC,ALL,TEMP,1000 !!!! =========== начальная температура, которая не почему-то задается / или задается, но её не видно ============ ALLSEL

Если для проверки задать температуру ну всю модель, тоже нет работает?

 

Тоже не работает.(

[4estor 1]

Упростил задачу больше некуда:

 

FINISH
/clear,start
/PREP7

RECTNG,0,1,0,1

ET,1,PLANE223
KEYOPT,1,1,11
KEYOPT,1,3,1

MPTEMP,,0
MPDATA,PRXY,1,,0.3
MPDATA,C,1,,0.5
MPDATA,KXX,1,,0.5
MPDATA,ENTH,1,,0.5
MPDATA,DENS,1,,7860
MPDATA,EX,i,,2e11
MPDATA,ALPX,i,,1e-5

ESIZE,0.1

AMESH,1

/SOLU
ANTYPE,4
NLGEOM,1
KBC,1
TIMINT,ON
TRNOPT,FULL
SOLCONTROL,ON

NSUBST,5,1e3,5
OUTRES,ALL,ALL

DK,1,UX
DK,1,UY

SFL,ALL,CONV,0.01,,20 ! Температура окр.среды 20 градусов

IC,ALL,TEMP,100 ! Начальная температура 100 градусов

TIME,1

SOLVE

/POST1
PLNSOL,TEMP,, 0

 

Результат такой же:

1) Температура равномерно возрастает с 0 до 20(тем.окр.среды) при любых свойствах и времени расчета.

2) Начальная температуры через IC опять не задается.

 

Что я делаю не так? Может что-то со свойствами или типом элемента?

[Борман 2 375]

Вчера пол вечера просидел над этим вопросом. На мысль натолкнула фраза из описания команды IC.

Be sure to define consistent initial conditions.

Есть предположение, что это начальное условие противоречит начальному распределению перемещений. Попробуйте выставить TREF=100. Потом, вероятно, начнутся проблемы со свойствами.... но вы их победите. В итоге у меня пластина все таки начала остывать.

[4estor 1]

Вчера пол вечера просидел над этим вопросом. На мысль натолкнула фраза из описания команды IC.

Be sure to define consistent initial conditions.

Есть предположение, что это начальное условие противоречит начальному распределению перемещений. Попробуйте выставить TREF=100. Потом, вероятно, начнутся проблемы со свойствами.... но вы их победите. В итоге у меня пластина все таки начала остывать.

 

 

Да, получилось. Это ведь было очевидно, как я сам не догадался. Спасибо, Борман!

Поварьировал потом свойства, начало остывать (но остывать - равномерно, видимо такие свойства)

Новый блок свойств, кому интересно:

MPTEMP,,0

MPDATA,PRXY,1,,0.3

MPDATA,C,1,,0.3

MPDATA,KXX,1,,1.5

MPDATA,ENTH,1,,100

MPDATA,DENS,1,,786

MPDATA,EX,1,,2e11

MPDATA,ALPX,1,,1e-10

TREF,100

 

 

 

Изменено 11 сентября 2013 пользователем 4estor

[Борман 2 375]

но остывать - равномерно

У меня остывало "от центра", как положено.

[4estor 1]

 

но остывать - равномерно

У меня остывало "от центра", как положено.

 

А скиньте пож-та макрос, где остывало от центра

[Борман 2 375]

Он уничтожен. Деятельность была направлена обнаружение глюка с IC.

[4estor 1]

Он уничтожен. Деятельность была направлена обнаружение глюка с IC.

Ок

[4estor 1]

И снова здравствуйте. Продвинулся в своей задаче, но опять уткнулся в проблему (
Задача усложнилась, если был квадратный слиток, то теперь сварное соединение.

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

 

FINISH
/clear,start
/PREP7

! исходные данные (начало)

l=200 ! длина ТП
d=800 ! диаметр ТП
t=40 ! толщина ТП
angdeg=40 ! угол в градусах
xv=5 ! длина шва внутри
nv=5 ! ширина внутренней наплавки
nn=4 ! ширина наружной наплавки

! исходные данные (конец)

ll=l/2
r=d/2
rn=(d/2)+t
ang=0.5*angdeg*3.14159/180 ! 20 градусов в радианах
xn=t*tan(ang)+(xv/2) ! длина (полная) шва снаружи

k,1,r,ll
k,2,rn,ll
k,3,rn,xn
k,4,r,xv/2

KSYMM,Y,1,4,1,,1,0

k,9,r,(xv/2)+nv ! внутренняя наплавка
k,10,rn,xn+nv

k,11,r,(xv/2)+nv+nn ! внешняя наплавка
k,12,rn,xn+nv+nn

a,1,2,12,11
a,5,6,7,8
a,3,4,8,7
a,3,4,9,10
a,9,10,12,11


! Задание свойств КЭ (начало)

ET,1,PLANE223
KEYOPT,1,1,11
KEYOPT,1,3,1

! Задание свойств КЭ (конец)


! Задание материалов (начало)


nmats=2
ntemps=6
ntemps14=14
puas=0.3
to_air=126 ! КОЭФФ,ТЕПЛООТДАЧИ СТАЛИ В ВОЗДУХ (ВТ/м2*С)
t_air=20

*SET,t
*DIM,t, ARRAY,nmats,ntemps

*SET,r02
*DIM,r02, ARRAY,nmats,ntemps

*SET,E
*DIM,E, ARRAY,nmats,ntemps14

*SET,al
*DIM,al, ARRAY,nmats,ntemps14

*SET,tm
*DIM,tm, ARRAY,nmats,ntemps14

*SET,c
*DIM,c, ARRAY,nmats,ntemps14

*SET,kxx
*DIM,kxx, ARRAY,nmats,ntemps14

*SET,enth
*DIM,enth, ARRAY,nmats,ntemps14

*SET,dens
*DIM,dens, ARRAY,nmats

*SET,hf
*DIM,hf, ARRAY,nmats,ntemps14



! материал 1 - перлит 10ГН2МФА - только верх

t(1,1)=20
t(1,2)=50
t(1,3)=150
t(1,4)=250
t(1,5)=350
t(1,6)=1500

r02(1,1)=343
r02(1,2)=333
r02(1,3)=314
r02(1,4)=304
r02(1,5)=294
r02(1,6)=100

E(1,1)=210
E(1,2)=207
E(1,3)=205
E(1,4)=202
E(1,5)=200
E(1,6)=197
E(1,7)=195
E(1,8)=190
E(1,9)=185
E(1,10)=180
E(1,11)=175
E(1,12)=170
E(1,13)=165
E(1,14)=100

al(1,1)=11
al(1,2)=11.5
al(1,3)=11.9
al(1,4)=12.2
al(1,5)=12.5
al(1,6)=12.8
al(1,7)=13.1
al(1,8)=13.4
al(1,9)=13.6
al(1,10)=13.8
al(1,11)=14
al(1,12)=14.2
al(1,13)=14.4
al(1,14)=20.3

tm(1,1)=5356
tm(1,2)=5141
tm(1,3)=4819
tm(1,4)=4793
tm(1,5)=4914
tm(1,6)=1000

c(1,1)=500
c(1,2)=500
c(1,3)=500
c(1,4)=500
c(1,5)=500
c(1,6)=500
c(1,7)=500
c(1,8)=500
c(1,9)=500
c(1,10)=500
c(1,11)=500
c(1,12)=500
c(1,13)=500
c(1,14)=500

kxx(1,1)=14
kxx(1,2)=15
kxx(1,3)=16
kxx(1,4)=17
kxx(1,5)=18
kxx(1,6)=19
kxx(1,7)=20
kxx(1,8)=21
kxx(1,9)=22
kxx(1,10)=23
kxx(1,11)=24
kxx(1,12)=26
kxx(1,13)=28
kxx(1,14)=30

hf(1,1)=50
hf(1,2)=50
hf(1,3)=50
hf(1,4)=50
hf(1,5)=50
hf(1,6)=50
hf(1,7)=50
hf(1,8)=50
hf(1,9)=50
hf(1,10)=50
hf(1,11)=50
hf(1,12)=50
hf(1,13)=50
hf(1,14)=50

enth(1,1)=1e8
enth(1,2)=2e8
enth(1,3)=3e8
enth(1,4)=4e8
enth(1,5)=6e8
enth(1,6)=8e8
enth(1,7)=1e9
enth(1,8)=2e9
enth(1,9)=4e9
enth(1,10)=5e9
enth(1,11)=6e9
enth(1,12)=8e9
enth(1,13)=9e9
enth(1,14)=1e10

dens(1)=7860

! материал 2 - аустенит 12Х18Н10Т - электрод, наплавки, низ

t(2,1)=20
t(2,2)=100
t(2,3)=200
t(2,4)=350
t(2,5)=600
t(2,6)=1500

r02(2,1)=196
r02(2,2)=189
r02(2,3)=181
r02(2,4)=167
r02(2,5)=147
r02(2,6)=100

E(2,1)=205
E(2,2)=202
E(2,3)=200
E(2,4)=195
E(2,5)=190
E(2,6)=185
E(2,7)=180
E(2,8)=175
E(2,9)=170
E(2,10)=167
E(2,11)=165
E(2,12)=162
E(2,13)=160
E(2,14)=100

al(2,1)=16.2
al(2,2)=16.4
al(2,3)=16.6
al(2,4)=16.8
al(2,5)=17
al(2,6)=17.2
al(2,7)=17.4
al(2,8)=17.6
al(2,9)=17.8
al(2,10)=18
al(2,11)=18.2
al(2,12)=18.4
al(2,13)=18.5
al(2,14)=20.3

tm(2,1)=1638
tm(2,2)=1609
tm(2,3)=1658
tm(2,4)=1685
tm(2,5)=1348
tm(2,6)=1000

c(2,1)=500
c(2,2)=500
c(2,3)=500
c(2,4)=500
c(2,5)=500
c(2,6)=500
c(2,7)=500
c(2,8)=500
c(2,9)=500
c(2,10)=500
c(2,11)=500
c(2,12)=500
c(2,13)=500
c(2,14)=500

kxx(2,1)=14
kxx(2,2)=15
kxx(2,3)=16
kxx(2,4)=17
kxx(2,5)=18
kxx(2,6)=19
kxx(2,7)=20
kxx(2,8)=21
kxx(2,9)=22
kxx(2,10)=23
kxx(2,11)=24
kxx(2,12)=26
kxx(2,13)=28
kxx(2,14)=30

hf(2,1)=50
hf(2,2)=50
hf(2,3)=50
hf(2,4)=50
hf(2,5)=50
hf(2,6)=50
hf(2,7)=50
hf(2,8)=50
hf(2,9)=50
hf(2,10)=50
hf(2,11)=50
hf(2,12)=50
hf(2,13)=50
hf(2,14)=50

enth(2,1)=1e8
enth(2,2)=2e8
enth(2,3)=3e8
enth(2,4)=4e8
enth(2,5)=6e8
enth(2,6)=8e8
enth(2,7)=1e9
enth(2,8)=2e9
enth(2,9)=4e9
enth(2,10)=5e9
enth(2,11)=6e9
enth(2,12)=8e9
enth(2,13)=9e9
enth(2,14)=1e10

dens(2)=7950

*DO,i,1,nmats
TB,BISO,i,ntemps
*DO,j,1,ntemps
TBTEMP,t(i,j)
TBDATA,,r02(i,j)*1e6,tm(i,j)*1e6
*ENDDO
*ENDDO


MPTEMP,,20,50,100,150,200,250
MPTEMP,7,300,350,400,450,500,550
MPTEMP,13,600,1500

*DO,i,1,nmats
MP,REFT,i,20
MPDATA,PRXY,i,,puas
MPDATA,DENS,i,,dens(i)

*DO,j,1,ntemps14
MPDATA,EX,i,,E(i,j)*1e9
MPDATA,ALPX,i,,al(i,j)*1e-6
MPDATA,C,i,,c(i,j)
MPDATA,KXX,i,,kxx(i,j)
MPDATA,ENTH,i,,enth(i,j)
! MPDATA,HF,i,,hf(i,j)
*ENDDO
*ENDDO


TBCOPY,ALL,2,3
MPCOPY,,2,3
templiq=1500 ! начальная температура (плавления)
MP,REFT,3,templiq


! Задание материалов (конец)



! Создание сетки (начало)

ndivrad=20 ! число разбиений в радиальном направлении
ndivax=4 ! число разбиений шва в осевом направлении

LESIZE,1,,,ndivrad
LESIZE,3,,,ndivrad
LESIZE,13,,,ndivrad
LESIZE,9,,,ndivrad
LESIZE,7,,,ndivrad
LESIZE,5,,,ndivrad

LESIZE,10,,,ndivax
LESIZE,11,,,ndivax

LESIZE,16,,,2
LESIZE,15,,,2
LESIZE,12,,,3
LESIZE,14,,,3

LESIZE,4,,,20,3
LESIZE,8,,,20,3
LESIZE,2,,,20,1/3
LESIZE,6,,,20,1/3

/NUMBER,1
/PNUM,MAT,1

MAT,1
AMESH,1

MAT,2
AMESH,2

MAT,3
AMESH,3

MAT,2
AMESH,4

MAT,2
AMESH,5

! Создание сетки (конец)


! Нагружение и решение (начало)

/SOLU
ANTYPE,4
NLGEOM,1
KBC,1
TIMINT,ON
TRNOPT,FULL
SOLCONTROL,ON,0
!TUNIF,t_air

DELTIM,1,1e-5,10 ! нач, мин, макс
OUTRES,ALL,ALL
AUTOTS,1

DK,5,UY

LSEL,S,,,1
LSEL,A,,,2
LSEL,A,,,15
LSEL,A,,,14
LSEL,A,,,11
LSEL,A,,,6
LSEL,A,,,5
LSEL,A,,,8
LSEL,A,,,10
LSEL,A,,,12
LSEL,A,,,16
LSEL,A,,,4
SFL,ALL,CONV,to_air,,t_air
ALLSEL

ASEL,U,,,3
ALLSEL,BELOW,AREA
IC,ALL,TEMP,t_air

ASEL,S,,,3
ALLSEL,BELOW,AREA
IC,ALL,TEMP,templiq
ALLSEL


TIME,100

SOLVE

 

 

 

 

Осесимметричная модель. Средняя часть остывает на воздухе с 1500 градусов до 20. Остальные части (верхние и нижняя) с начальной температурой 20 градусов.

Механические и тепловые свойства из справочников - сталь. Задание HF (нужно оно ли?) закомментил в макросе.

 

В чем проблемы:

1) Нефизичное распределение температуры по времени: то увеличивается, то уменьшается. И причем на "холодных" частях вообще отрицательная она. (отриц.темп. и в помине нет в макросе).

2) Решение идет очень мелкими шагами (по 10^-4 сек) или вообще не сходится в процессе решения.

 

Подозреваю, что дело в коэффициенте теплоотдачи стали в воздух: задаю 126 Вт/м2*С (значение из инета). Пробовал с другими значениями - модель себя по-разному ведет (сужу по температуре в разное время в Results Viewer).

[soklakov 1 475]

Transient решение?

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

[4estor 1]

Transient решение?

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

Да, Transient, with large deformations.

Параметры решателя

ANTYPE,4

NLGEOM,1

KBC,1

TIMINT,ON

TRNOPT,FULL

SOLCONTROL,ON,0

DELTIM,1,1e-5,10

OUTRES,ALL,-500

AUTOTS,1

 

Можете пояснить, по какой именно формуле? И, если можно, какую-нибудь наводку на освещение этой проблемы в хэлпе ансиса.

При увеличении допустимого (минимального, если правильно понял) решение не сходится в самом начале.

При уменьшении - решение опять идет микрошагами. Например 2 часа ждал, пока полсекунды задачи посчитает. А нужно хотя бы час, чтобы остыло.

[soklakov 1 475]

Можете пояснить, по какой именно формуле? И, если можно, какую-нибудь наводку на освещение этой проблемы в хэлпе ансиса.

Формулка есть в курсе Heat Transfer на портале.

 

(минимального, если правильно понял)

Нет. Начального. Initial Time Step.

 

ADPL, действительно, лень смотреть. Был бы проект WB, там еще можно было бы глянуть. Задача, вроде, не самая сложная.

[4estor 1]

 

Можете пояснить, по какой именно формуле? И, если можно, какую-нибудь наводку на освещение этой проблемы в хэлпе ансиса.

Формулка есть в курсе Heat Transfer на портале.

 

(минимального, если правильно понял)

Нет. Начального. Initial Time Step.

 

ADPL, действительно, лень смотреть. Был бы проект WB, там еще можно было бы глянуть. Задача, вроде, не самая сложная.

weld.gif

 

 

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

Вау, а это моя задача или чужая? Стремлюсь к такому же решению... Интересно, какие временные параметры расчета (t нач, мин, макс, конца решения).

 

А что за портал? Искал по форуму, в хэлпе и в интернете, курса и формулы этих не нашел, только на коммерческих сайтах идет запись на предстоящие курсы. Буду благодарен, если укажете по-конкретнее, где эта формула.

 

Начальный шаг сейчас пробую изменять, безрезультатно. Если увеличить, то в процессе решения программа будет его "рубить" в 2 раза до минимального 10^-4 сек.

 

И еще: наткнулся на Ваши ответы в теме с аналогичной проблемой http://fsapr2000.ru/topic/52052-otritcatelnye-znacheniia-temperatur/ Там увеличение шага помогло, а мне вот получается нет =(

И какие есть 2D элементы более низкого порядка для решения термомеханических задач? Кроме 2D 8-узловых PLANE223 ничего не нашел...

[soklakov 1 475]

а это моя задача или чужая?

Это я набросал минут за десять. Чего и Вам желаю.

А что за портал?

Портал пользователей Ansys

Но формула Вам ни к чему. За несколько итераций вполне подберете шаг.

[4estor 1]

Перепробовал нач.шаги 1e-2, 1e-1, 1, 2, 5, 10, 20, 50, 100, 200, 500, 1000, 1e4, 1e5, 1e6. Даже 1e9. Все тщетно. Сразу же на первой итерации рубит пополам до 10^-4 или 10^-5 (до мин.шага). Максимальный шаг - 1e4 или как начальный. Время окончания расчета 1e6.

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

 

Выбрал PLANE13 (подсказали 4х узловой для термомеханики). Результаты физичные! Модель начала остывать, но едва заметно. Хотя свойства почти все из справочников. Буду разбираться с настройками решателя и свойствами. И решение довольно быстро прерывается (не сходится).

Изменено 15 октября 2013 пользователем 4estor