热电耦合求解温度分布

yutingting1986

我做的是两电极间加电流/电压对电极间介质的加热问题，但是求解不能收敛，请教各位高手怎么解决，下面是命令流，谢谢
/show
/filename,redian
/title,thermal-electric magnetic fluid
/prep7
/units,si
et,1,solid69
et,2,solid69
mptemp,1,300,400,500,600,700,800
mptemp,7,900,1000,1100,1200,1300,1400
mptemp,13,1500,1600,1700,1800,1900,2000
mptemp,19,2100,2200,2300,2400,2500,2600
mptemp,25,2700,2800,2900,3000,3100,3200
mptemp,31,3300,3400,3500,3600,3700,3800
mptemp,37,3900,4000,4100,4200,4300,4400
mptemp,43,4500,4600,4700,4800,4900,5000
mptemp,49,5100,5200,5300,5400,5500,5600
mptemp,55,5700,5800,5900,6000,6100,6200
mptemp,61,6300,6400,6500,6600,6700,6800
mptemp,67,6900,7000,7100,7200,7300,7400
mptemp,73,7500,7600,7700,7800,7900,8000
mptemp,79,8100,8200,8300,8400,8500,8600
mptemp,85,8700,8800,8900,9000,9100,9200
mptemp,91,9300,9400,9500,9600,9700,9800
mpdata,kxx,1,1,0.025,0.033,0.042,0.048,0.052,0.058
mpdata,kxx,1,7,0.063,0.069,0.072,0.078,0.083,0.088
mpdata,kxx,1,13,0.093,0.098,0.102,0.109,0.114,0.121
mpdata,kxx,1,19,0.15,0.18,0.21,0.24,0.268,0.295
mpdata,kxx,1,25,0.32,0.345,0.367,0.387,0.404,0.42
mpdata,kxx,1,31,0.434,0.464,0.457,0.468,0.477,0.488
mpdata,kxx,1,37,0.499,0.512,0.528,0.547,0.57,0.597
mpdata,kxx,1,43,0.63,0.67,0.72,0.77,0.84,0.911
mpdata,kxx,1,49,0.994,1.087,1.19,1.3,1.43,1.57
mpdata,kxx,1,55,1.71,1.86,2.02,2.19,2.365,2.54
mpdata,kxx,1,61,2.71,2.887,3.05,3.21,3.35,3.41
mpdata,kxx,1,67,3.47,3.52,3.35,3.15,2.98,2.816
mpdata,kxx,1,73,2.66,2.5,2.36,2.22,2.09,1.97
mpdata,kxx,1,79,1.85,1.74,1.64,1.54,1.45,1.376
mpdata,kxx,1,85,1.3,1.23,1.17,1.12,1.12,1.12
mpdata,kxx,1,91,1.12,1.12,1.12,1.12,1.12,1.12
mpdata,rsvx,1,1,5.59e20,1.5e19,4.88e17,1.88e16,8.58e14,4.6e13
mpdata,rsvx,1,7,2.88e12,2.1e11,1.77e10,1.71e9,1.88e8,2.36e7
mpdata,rsvx,1,13,3.35e6,5.33e5,9.47e4,1.87e4,4.1e3,9.86e2
mpdata,rsvx,1,19,2.6e2,74.8,23.3,7.87,2.86,1.11
mpdata,rsvx,1,25,0.461,0.204,9.52e-2,4.7e-2,2.44e-2,1.32e-2
mpdata,rsvx,1,31,7.53e-3,4.46e-3,2.74e-3,1.75e-3,1.15e-3,7.8e-4
mpdata,rsvx,1,37,5.44e-4,3.89e-4,2.85e-4,2.13e-4,1.62e-43,1.26e-4
mpdata,rsvx,1,43,9.95e-5,7.96e-5,**5e-5,5.29e-5,4.38e-5,3.67e-5
mpdata,rsvx,1,49,3.09e-5,2.63e-5,2.25e-5,1.93e-5,1.67e-5,1.45e-5
mpdata,rsvx,1,55,1.26e-5,1.11e-5,9.72e-6,8.56e-6,7.55e-6,6.68e-6
mpdata,rsvx,1,61,5.92e-6,5.26e-6,4.68e-6,4.17e-6,3.72e-6,3.33e-6
mpdata,rsvx,1,67,2.99e-6,2.68e-6,2.42e-6,2.18e-6,1.97e-6,1.79e-6
mpdata,rsvx,1,73,1.62e-6,1.48e-6,1.35e-6,1.24e-6,1.14e-6,1.05e-6
mpdata,rsvx,1,79,9.73e-7,9.05e-7,8.42e-7,7.87e-7,7.39e-7,6.96e-7
mpdata,rsvx,1,85,6.58e-7,6.24e-7,5.95e-7,5.68e-7,5.45e-7,5.24e-7
mpdata,rsvx,1,91,5.06e-7,4.9e-7,4.75e-7,4.63e-7,4.51e-7,4.41e-7
mp,kxx,2,163.3
mp,rsvx,2, 5.65e-8
k,1,0,0,0.0128
k,2,0.015,0,0.0128
k,3,0,0,0.0228
k,4,0.015,0,0.0228
A,1,2,4,3
Vrotat,1,,,,,,1,3
cyl4,,,,,0.0075,360,0.0042
cone,0.0075,0.0025,0.0042,0.0073,0,360
local,11,1,,,.0073
wpave
cyl4,,,,,.0025,360,0.003
cyl4,,,,,0.0015,360,0.0016
local,12,1,,,0.0228
wpave
cyl4,,,,,0.015,360,0.04
csys,0
wpave
cone,0.0025,0.005,0.0103,0.0128,,360
cyl4,,,,,0.003,360,0.002
cone,0.003,0.0005,0.002,0.0073,,360
GUI操作：用体7,6,5减去体8,11,12
cyl4,,,,,0.003,360,0.002
cone,0.003,0.0005,0.002,0.0073,,360
csys,11
wpave
cyl4,,,,,0.0015,360,0.0016
csys,0
wpave
csdele,11
csdele,12
numcmp,volu
vglue,all
numcmp,all
vsel,s,volu,,1
vatt,2,,2
vsel,s,volu,,7
vatt,2,,2
vsel,s,volu,,2,6
vatt,1,,1
vsel,s,volu,,8,12
vatt,1,,1
allsel
mshape,1,3d
smrtsize,1
vmesh,all
/solu
antype,static
toffst,0
tunif,273
!time,1000
*ABBR,NSL,NSLA,,1
asel,s,,,48       !因为电阻率在温度较低的时候很大，电极间无法导电，所以施加了初始载荷，不知道行否？
nsl
d,all,temp,393
asel,s,,,46,47
nsl
d,all,temp,673
asel,s,,,44,45
nsl
d,all,temp,1273
asel,s,,,41,42
nsl
d,all,temp,1952
asel,s,,,16,17
nsl
d,all,temp,1473
asel,s,,,31
nsl
d,all,temp,673
asel,s,,,34
nsl
d,all,temp,673
asel,s,,,24
nsl
d,all,temp,673
asel,s,,,28
nsl
d,all,temp,673
asel,s,,,1
nsl
d,all,temp,393
asel,s,,,3
nsl
d,all,temp,393
asel,s,,,5
nsl
d,all,temp,393
asel,s,,,7
nsl
d,all,temp,393
asel,s,,,14
nsl
d,all,temp,393
asel,s,,,39,40
nsl
d,all,temp,393
asel,s,,,13
nsl
d,all,temp,2073
asel,s,,,35,38
nsl
d,all,temp,1873
asel,s,,,11,12
nsl
d,all,temp,1573
vsel,s,,,2
nslv,s,1
d,all,temp,2785
allsel
asel,s,,,44,45         !阳极边界加电压0
nsl
d,all,volt,0
asel,s,,,41,42
nsl
d,all,volt,0
asel,s,,,16,17
nsl
d,all,volt,0
asel,s,,,13              !阴极加正电压，也可以是电流载荷，电流为1.5A，此处加电压
nsl
d,all,volt,70
asel,s,,,36
nsl
d,all,volt,70
asel,s,,,38
nsl
d,all,volt,70
allsel
asel,s,,,1
bfa,all,ef,,0
asel,s,,,3
bfa,all,ef,,0
asel,s,,,5
bfa,all,ef,0
asel,s,,,7
bfa,all,ef,,,0
asel,s,,,14
bfa,all,ef,,0
asel,s,,,39,40
bfa,all,ef,,0
asel,s,,,9
bfa,all,ef,0
asel,s,,,48
bfa,all,ef,,,0
allsel
sbctran
solve

binloving

是不是没有对实体做自由度耦合

yutingting1986

想问一下自由度耦合的作用是什么？之前加电流的时候进行过电压的自由度耦合，但是好像没什么用呢

yutingting1986

binloving 发表于 2011-9-6 20:41
是不是没有对实体做自由度耦合

想问一下自由度耦合的作用是什么？之前加电流的时候进行过电压的自由度耦合，但是好像没什么用呢

binloving

对自由度耦合后才能得到加载电流的节点吧，我也是初学者，也不是很清楚呵呵

xhj543240343

另外，由于你的材料随温度变化的值太多了，也可能影响收敛。你修改一下试一试。此外，热电耦合，一般都需要对电流输入端进行耦合。