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

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


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

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

 

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

 

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

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

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

 

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

Изменено пользователем 4estor
Ссылка на сообщение
Поделиться на других сайтах


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

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

Ссылка на сообщение
Поделиться на других сайтах

 

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

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

 

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

Ссылка на сообщение
Поделиться на других сайтах

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

 

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 опять не задается.
 
Что я делаю не так? Может что-то со свойствами или типом элемента?
Ссылка на сообщение
Поделиться на других сайтах

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

Be sure to define consistent initial conditions.

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

Ссылка на сообщение
Поделиться на других сайтах

Вчера пол вечера просидел над этим вопросом. На мысль натолкнула фраза из описания команды 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

 

 

 

Изменено пользователем 4estor
Ссылка на сообщение
Поделиться на других сайтах
но остывать - равномерно

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

Ссылка на сообщение
Поделиться на других сайтах

 

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

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

 

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

Ссылка на сообщение
Поделиться на других сайтах

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

Ссылка на сообщение
Поделиться на других сайтах
  • 1 месяц спустя...

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

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

 

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

 

 

LbQB-DkL3Kc.jpg

Jf7-cx2VM4o.jpg

 

 

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

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

 

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

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

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

 

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

Ссылка на сообщение
Поделиться на других сайтах

Transient решение?

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

Ссылка на сообщение
Поделиться на других сайтах

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 часа ждал, пока полсекунды задачи посчитает. А нужно хотя бы час, чтобы остыло.

Ссылка на сообщение
Поделиться на других сайтах

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

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

 

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

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

 

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

post-33642-0-03304500-1381826095.gif

Ссылка на сообщение
Поделиться на других сайтах

 

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

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

 

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

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

 

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

attachicon.gifweld.gif

 

 

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

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

 

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

 

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

 

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

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

Ссылка на сообщение
Поделиться на других сайтах
а это моя задача или чужая?

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

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

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

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

Ссылка на сообщение
Поделиться на других сайтах

Перепробовал нач.шаги 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х узловой для термомеханики). Результаты физичные! Модель начала остывать, но едва заметно. Хотя свойства почти все из справочников. Буду разбираться с настройками решателя и свойствами. И решение довольно быстро прерывается (не сходится).

Изменено пользователем 4estor
Ссылка на сообщение
Поделиться на других сайтах

Присоединяйтесь к обсуждению

Вы можете опубликовать сообщение сейчас, а зарегистрироваться позже. Если у вас есть аккаунт, войдите в него для написания от своего имени.
Примечание: вашему сообщению потребуется утверждение модератора, прежде чем оно станет доступным.

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставить изображения напрямую. Загрузите или вставьте изображения по ссылке.

  • Сейчас на странице   0 пользователей

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




  • Сообщения

    • Bot
      АСКОН, российский разработчик инженерного программного обеспечения и ИТ-интегратор, подвел итоги 2023 года. Выручка компании выросла на 47% и достигла 5,6 млрд рублей, штат сотрудников увеличился до 1250 человек. Клиентская база пополнилась тысячью предприятий и сегодня с решениями АСКОН работают почти 16 000 заказчиков из всех отраслей промышленности и строительства. Смотреть полностью
    • Sergey78
    • scrimline
      Уговаривать ему пришлось да и оскорбили его. Я вообще то сразу написал что постараюсь ладдер скинуть, а ты сарказмом со всех щелей хлестать начал. 
    • andrey2147
      Я подключал 532 , позже поищу ссылки, ВинПсином и кабелем там хитрость в настройках винпсина 7 бит надо ставить вам зачем ?
    • Guhl
      Надеюсь, в моих словах вы ничего оскорбительного не усмотрели Или это я?
    • Guhl
      Представляю. Я же не предлагал на полной скорости считать Получается, что до м19 считает правильно, а после - неправильно?
    • mannul
      Вы сначала проверьте точно ли режущий край пластины в оси детали. А то китайские державки и пластины бывают с такими косяками.
    • kkk
      Что-то вдруг резко понадобилась возможность посчитать площадь разной отделки в изделии. Есть куча многотельных деталей, какие-то поверхности красятся одним цветом, какие-то другим, какие-то третьим. Раскрасить в модели проблем нет, а нет ли возможности высчитать эту информацию автоматом? Есть ли какой-то способ наложения "отделки", чтобы Солид мог ей оперировать в дальнейшем.    Раньше, если изредка нужно было что-то подобное посчитать, просто выбирал нужные грани мышкой и смотрел сумму площадей в инструменте измерения. Но хочется больше интерактива, чтобы Солид сам отслеживал, что у него чем "покрашено".
    • andrey2147
      Пришлось уговаривать оператора что-бы выполнил просьбу Виктора, ещё и  оскорбили.  
    • gudstartup
      да не надо ничего паять система прекрасно считает импульсы на оборот в 445 параметре диагностики вот только беда в том что один раз войдя в режим контроля позиционирования после m19 она из него уже не выходит хотя сразу после включения находится в режиме управления скоростью и 445 параметр ничего не отображает и контроль числа импульсов отсутствует. Вопрос как выключить контроль позиционирования. он вводится по G70.6 но и после того как этот сигнал остается система продолжает контролировать кодер. в обще то понятно что это датчик так как на скорости 100 или 200 оборотов а также при ориентации все хорошо а потом видимо амплитуда сигнала падает ниже нормы там все ок вы частоту представляете при 5000 об 5000*4096
×
×
  • Создать...