电流分布不正确,究竟是怎么错了??(原程序)

zhangzhenbin

问题,把两个电极插入另一种导电不同的介质中,求电流分布!!但是模拟的结果如附件,很显然是错的,因为电流应该在表面处最大,模拟的结果却在中间处电流最大.希望高手帮帮忙,多谢了!!!!!不同网格密度我也实验过了!!都是这样
/BATCH  
/title,lunwen
/filename,zzb,1
/PREP7  
!!创建单元
ET,1,SOLID69
MP,RSVX,1,4.39e-6!定义电极电阻率
mp,murx,1,1 !定义电极磁导率
MP,RSVX,2,8e-6!定义熔体电阻率
mp,murx,2,1!定义熔体电极磁导率
EMUNIT,MKS
CYL4,0,0,0.1, , , ,0.3 !定义坩埚半径10CM，高30CM
CYL4,0.05,0,0.005, , , ,0.20!定义电极半径5mm，高20CM
CYL4,-0.05,0,0.005, , , ,0.20!定义电极半径5mm，高20CM
vsel, s,volu, ,2,3
VGEN, ,all, , , , ,0.15, , ,1  
vsel, all
VOVLAP,all
vsel, s,volu, ,4,7
vatt,1,,1,0!!电极
vsel, s,volu, ,8
vatt,2,,1,0!!熔体
smrtsize,6
mshkey,0
mshape,1,3d
vsel,s,volu,,4,8
vmesh,all
finish
/solu
antyp,0
!!!!!!!!!!!!!!!
ASEL,S,area, ,10!!! 对的
NSLA,S,1
D,ALL,VOLT,0
!!!!!!!!!!!!!!!!!!!!
ASEL,S,area, ,6!!! 对的
NSLA,S,1
CP,11,VOLT,ALL   !耦合自由度
!d,all,volt,50
F,893,AMPS,10000   !施加总电流
!!!!!!!!!!!!!!!!!!!!!!!
ALLSEL,all                     
solve
finish

[ 本帖最后由 zhangzhenbin 于 2007-7-6 10:27 编辑 ]

cowry

JT 和 JS 分别指电流的实部和虚部
在谐波分析或者顺态分析才会都有
静态的时候只有一个

cowry

你的结果应该也不算错误
两个熔丝上面的电流分布基本均匀
而且和最大值在一个数量级上误差不大
理论上应该是熔丝上面最大了
只是你用的是自由划分的网格，所以误差可能看起来大了点
但是不应该算错误
如果用规则的六面体划分网格误差可能会小点

cowry

说过了
你的结果不能算是错误
只是些许误差让云图结果看起来不是很舒服而已
改变网格可能会让误差小点，但也不能避免误差的

圆柱体肯定可以划分六面体网格
但是你的模型的大圆柱上表面不是一个面所以不能用mape方式划分而已

cowry

先把复杂的那个面分好
把高方向的段数分好
在mesh里有项操作把复杂的那个面上的三个面加成一个(contran...)
然后选择起止面扫略也可以的

perch88

没有收敛，首先检查边界条件

ASEL,S,LOC,Z,0.31
ASEL,R,LOC,X,0.05-0.005,0.005+0.005
NSLA,S,1
D,ALL,VOLT,0
什么都没选上，没有设置参考电位。

改成
ASEL,S,LOC,Z,0.31
asel,s,,,81
asel,a,,,84
asel,a,,,86
asel,a,,,88
NSLA,S,1
D,ALL,VOLT,0

zhangzhenbin

版主,各位朋友!!我按照大家的建议对程序进行了一晚上的改善,现在有结果.可是看上去还是错误的,还不如第一次的那个呢??为什么??下面附件电流矢量的结果,很显然,还是有问题的!!究竟怎么回事????别的电势分布,路径好像还是对的,邪不邪????
/BATCH
/CWD,'D:\ANSYSxuexi'
/title,lunwen
/filename,JJJJ,1
/PREP7  
!!创建单元
ET,1,SOLID69
MP,RSVX,1,4.39e-6!定义电极电阻率
mp,murx,1,1 !定义电极磁导率
MP,RSVX,2,8e-6!定义熔体电阻率
mp,murx,2,1!定义熔体电极磁导率
EMUNIT,MKS
CYL4,0,0,0.1, , , ,0.3 !定义坩埚半径10CM，高30CM
CYL4,0.05,0,0.005, , , ,0.31!定义电极半径5mm，高20CM
CYL4,-0.05,0,0.005, , , ,0.31!定义电极半径5mm，高20CM
vsel,all
VOVLAP,all
wpro,,,90
vsbw,all
wpof,,,0.05
vsbw,all
wpof,,,-0.10
vsbw,all
wpro,,90,
vsbw,all
wpof,-0.15,,
wpro,,,90
vsbw,all
VSEL,S,LOC,X,-0.05-0.005,0.05+0.005
VSEL,R,LOC,Y,-0.005,0.005
VSEL,R,LOC,Z,0.15,0.31
vatt,1,,1,0!!电极
VSEL,INVE
vatt,2,,1,0!!熔体
ACCAT,154,129
ACCAT,148,126
ACCAT,142,29
ACCAT,11,136
ACCAT,131,4
ACCAT,139,16
ACCAT,144,124
ACCAT,130,151
VSEL,ALL
mshkey,1
mshape,0,3d
vmesh,all
finish
/solu
antyp,0
!!!!!!!!!!!!!!!
ASEL,S,LOC,Z,0.31
ASEL,R,LOC,X,0.05-0.005,0.05+0.005
NSLA,S,1
D,ALL,VOLT,0
!!!!!!!!!!!!!!!!!!!!
ASEL,S,LOC,Z,0.31!! 对的
ASEL,R,LOC,X,-0.05-0.005,0.005-0.05
NSLA,S,1
CP,11,VOLT,ALL   !耦合自由度
!d,all,volt,50
F,6496,AMPS,10000   !施加总电流
!!!!!!!!!!!!!!!!!!!!!!!
ALLSEL,all                     
solve
finish

[ 本帖最后由 zhangzhenbin 于 2007-7-9 07:11 编辑 ]

zhangzhenbin

顶一下，要不要沉了

zhangzhenbin

deing ding ding

lilulilu

JS:Total element current density. It is the sum of the element conduction and the displacement current densities. It may be used as a source for a subsequent magnetostatic analysis.
JT:Element conduction current density.
JC:Nodal conduction current density.

现在的结果是电流与电势垂直了。
哪里还不对？

zhangzhenbin

根据电势的计算结果，电流应主要集中在熔体表面附近，熔体底部没有电流，可是从电流矢量图来看（JJJJ001.jpg ），却看不出这种趋势。结果好像对，但是显示太滥了！！

北极熊甲

原帖由 zhangzhenbin 于 2007-7-11 05:24 发表
根据电势的计算结果，电流应主要集中在熔体表面附近，熔体底部没有电流，可是从电流矢量图来看（JJJJ001.jpg ），却看不出这种趋势。结果好像对，但是显示太滥了！！

呵呵，这很正常，只有9种色彩。你把device设置为Win32C，把countour的数目增加，最多128。效果会好得多。

[ 本帖最后由 北极熊甲 于 2007-7-11 11:29 编辑 ]

zhangzhenbin

这次用自由网格剖分实现了,但是映射网格剖分结果却显示的不好(因为上下电流密度差别很大,可是矢量图上却显示不出来).大家讨论讨论吧!!附件第一个为自由网格剖分的显示效果,第二个为映射网格的显示效果!
另外,版主建议改变等值线的条数,我试验过了,没有什么改善!!

[ 本帖最后由 zhangzhenbin 于 2007-7-13 15:50 编辑 ]

hudson57

我之所以顶一下，是想知道楼主最后算的到底怎么样了