【原创】热电耦合- 通电金属桥的温升

iomega

在这个例子中， 电流通过suspended金属桥，Joule heating引起温升，并伴随相变。例子很简单，只是说明如何做热电耦合。

FINISH
/CLEAR
k=95
DEN=3900
Cp=450
/PREP7
um=1e-6
L=100*um
W=10*um
T=0.05*um
mA=0.001
rou=1e-7
I=3.1624*mA
DT=5
FUSION=15*Cp

T0=480

T_MELT=500

T_LOW=0
T_HIGH=1.5*T_MELT

Size=1*um

  
!*  
ET,1,SOLID69

UIMP,1,DENS, , ,DEN,
UIMP,1,KXX,KYY,KZZ,k,k,k,
UIMP,1,RSVX,RSVY,RSVZ,rou,rou,rou
!UIMP,1,C, , ,Cp,
MPTEMP,1,T_LOW,T_MELT-DT/2,T_MELT+DT/2,T_HIGH                  
MPDATA,ENTH,1,1,0,DEN*Cp*(T_MELT-DT/2-T_LOW),DEN*(Cp*(T_MELT-DT/2-T_LOW)+FUSION),DEN*(Cp*(T_MELT-DT/2-T_LOW)+FUSION+Cp*(T_HIGH-T_MELT-DT/2))

BLC4, 0, 0, W, L, T

lsel,s,LOC,x,0,0
lsel,r,LOC,y,0,L
lesize,all,Size,
allsel,all

lsel,s,LOC,y,0,0
lsel,r,LOC,x,0,W
lesize,all,W,
allsel,all


TYPE,1
MAT,1
mshkey,1
mshape,0,3d
vmesh,1

/solu
tunif,T0
TREF,T0
asel,s,LOC,y,0,0
asel,a,LOC,y,L,L
DA,all,TEMP,T0   !Two ends of the bridge are kept at T0

allsel,all

asel,s,LOC,y,0,0  ! Assume one end V=0
DA,ALL,VOLT,0
asel,all

nsel,s,LOC,y,L,L
cp,2,volt,all        ! Couple VOLT dofs
ni=ndnext(0)         ! Get master node
nsel,all

VSEL,ALL
f,ni,amps,I      !Apply current at one end
allsel,all

ANTYPE,4
SOLCONTROL,OM,0

TIME,1e-4   
AUTOTS,-1   
DELTIM,1e-6, , ,1   
KBC,1   
OUTRES,ALL,ALL,
SOLVE

nwpuyl

学习学习先！

ilovemm

MPTEMP,1,T_LOW,T_MELT-DT/2,T_MELT+DT/2,T_HIGH  !dt=5??这里是怎么设置的4点？？               
MPDATA,ENTH,1,1,0,DEN*Cp*(T_MELT-DT/2-T_LOW),DEN*(Cp*(T_MELT-DT/2-T_LOW)+FUSION),DEN*(Cp*(T_MELT-DT/2-T_LOW)+FUSION+Cp*(T_HIGH-T_MELT-DT/2))！FUSION=15*Cp是什么？为什么要加上？

iomega

The below code is to mimic the phase change, e.g. the fusion heat when the solid become liquid.

When fusion happenes, the material absorbs the heat but the temperature doesn't increase until it becomes liquid. This is physically equivalent to the case that material has infinite heat capacity at melting point. I set a small temperature dt around the melting point: T_MELT,   let dt*Cp =   fusion heat.

原帖由 ilovemm 于 2007-1-31 18:43 发表
MPTEMP,1,T_LOW,T_MELT-DT/2,T_MELT+DT/2,T_HIGH  !dt=5??这里是怎么设置的4点？？               
MPDATA,ENTH,1,1,0,DEN*Cp*(T_MELT-DT/2-T_LOW),DEN*(Cp*(T_MELT-DT/2-T_LOW)+FUSION),DEN*(Cp*(T_MELT-DT/2- ...

.

icst_li

谢谢了.正愁找不到这方面的例子呢.

ixxing

顶，好东西

huomud

谢谢楼主的例子

terryliu_1209

好贴，要顶啊！正在学习中！

alex_shoal

多谢。学习先

wzwgyq

我现在做的东西也是由电产生焦耳热，初学者看不懂你写的啊，能说的具体些吗？

happyhxj

好东西，正在学习，谢谢楼主

yeziting12

好，支持

comego

好东西正要学习一下，多谢！！
UIMP,1,KXX,KYY,KZZ,k,k,k,
UIMP,1,RSVX,RSVY,RSVZ,rou,rou,rou
请问KXX，和RSVX这两个参数是不是一定要把三个方向都定义了才可以正常计算？

guiguzi800

好东东就大家分享啊

qsmm1912

我现在做的东西

xhj543240343

好东西，多谢楼主！

yyzhlove163

好东西  先学习

chenholy

新手上路，顶上！

livelife

温升带相变的例子，急需

outman1302

cp,2,volt,all        ! Couple VOLT dofs
why is 2? not 1