洛仑兹力计算为何为0？3D模型
模型为两个平行导体，激励为电流，我查了以前的这类例子都是2D的，所以想改用117单元做个3D的看看，可不知道为什么洛仑兹力总上为0，附上命令流：/prep7
et,2,117,5 !air
et,1,117,5 !conductor
mp,murx,2,1 !air
mp,murx,1,1 !conductor
mp,rsvx,1,4e-7 !bus
block,-0.1475,-0.1275,0,0.4015,-0.050,0.050 !conductor1 model
block,0.1275,0.1475,0,0.4015,-0.050,0.050 !conductor2 model
block,-0.6375,0.6375,0,0.4015,-0.400,0.400 !air model
vovlap,all
numcmp,volu
vsel,s,volu,,4
vatt,2,1,2 !air properties
vsel,s,volu,,2
vsel,a,volu,,1
vatt,1,1,1 !conductor properties
allsel,all
smart,6
mshape,1,3d
mshkey,0
vmesh,all
*do,i,1,2,1
vsel,s,volu,,%i%
allsel,below
nplot
*get,ymax,node,,mxloc,y,
*get,ymin,node,,mnloc,y,
nsel,r,loc,y,ymin
cp,next,volt,all
*get,nn,node,,num,max
f,nn,amps,24000
nsel,s,ext
nsel,r,loc,y,ymax
d,all,volt,0
*enddo
alls
nsel,s,ext
d,all,az,0
/solu
alls
save
cnvtol,amps,1,1.3e-3
antype,static
eqslv,sparse
solve
/post1
set,1,1
etable,fmagx,fmag,x
etable,fmagy,fmag,y
etable,fmagz,fmag,z
ssum
esel,,mat,,1
nsle,s
plvect,jt,,,,vect,elem,on,0
[[i] 本帖最后由 maomaoface 于 2008-6-24 08:16 编辑 [/i]] 试着选择单根导体的单元求洛仑兹力看看。 [quote]原帖由 [i]infrom[/i] 于 2008-6-22 10:11 发表 [url=http://www.simwe.com/forum/redirect.php?goto=findpost&pid=1364573&ptid=837950][img]http://www.simwe.com/forum/images/common/back.gif[/img][/url]
试着选择单根导体的单元求洛仑兹力看看。 [/quote]
敢问版主，为什么呢？我加载的算是直流激励，选单根导体更没有结果了。还请版主帮我跑一遍命令流，看看哪里出了问题呀。
现在的结果是查看电流密度是有的，但是不对，而查看B，H就没结果。因此洛仑兹力为0了 选择单个导体，然后选择这个导体上的所有节点，然后选择整个模型，再用emft命令试试！在post1下输入
emft这个命令
[quote]原帖由 [i]giggsun[/i] 于 2008-6-22 16:19 发表 [url=http://www.simwe.com/forum/redirect.php?goto=findpost&pid=1364856&ptid=837950][img]http://www.simwe.com/forum/images/common/back.gif[/img][/url]选择单个导体，然后选择这个导体上的所有节点，然后选择整个模型，再用emft命令试试！在post1下输入 [/quote]
谢谢你的回复！
我在8.0中的help查不到emft这个命令，请问是版本的问题吗？ 机子上没有装ansys，没有跑命令流，看了一下命令流，感觉有些问题：1：单元定义上有问题，et,2,117,5 !air，空气部分的单元定义有问题，solid117 keyopt(1)=5,自由度为（AZ+VOLT），通电导体加载总电流单元这样定义，空气自由度为AZ即可，定义成et,2,117计算看看。2：后处理查看罗伦兹力时，选择一个导体查看，因为两个导体之间的罗伦兹力是作用力和反作用力，全部一起选择的话合力为0。
多谢回复
[quote]原帖由 [i]zhouyz[/i] 于 2008-6-23 09:46 发表 [url=http://www.simwe.com/forum/redirect.php?goto=findpost&pid=1365215&ptid=837950][img]http://www.simwe.com/forum/images/common/back.gif[/img][/url]机子上没有装ansys，没有跑命令流，看了一下命令流，感觉有些问题：1：单元定义上有问题，et,2,117,5 !air，空气部分的单元定义有问题，solid117 keyopt(1)=5,自由度为（AZ+VOLT），通电导体加载总电流单元这 ... [/quote]
谢谢！我已经按您说的做了修改，但问题仍未解决，始终是没有B和H，即便单独选出一个导体，还是没有洛仑兹力。很奇怪，但我把空气部分也的材料设置成和导体一样，就得到了正确的JT分布，但仍没有磁场。。。 程序改动部分：
et,2,117 !air
et,1,117 !conductor
加载电流部分改为加载电流密度：
s=(0.1475-0.1275)*0.1
JY=24000/S
vsel,s,volu,,1
allsel,below
BFE, ALL, JS, , , JY
vsel,s,volu,,2
allsel,below
BFE, ALL, JS, , , JY
求解部分改为：
/solu
allsel
solve 结果图形如下：
[[i] 本帖最后由 infrom 于 2008-6-24 16:01 编辑 [/i]] 上图为俯视的磁感应强度矢量图；
下图为前视的节点力矢量图。 感谢infrom版主，你得到的结果确实是我想要的，只是为什么用电流加载的方式就得不到B的值呢？当导体形状不规则时，用js加载可能不大方便啊。
这是用电流加载方式得到的JT结果，但是得不到B 我也试了好几次，确实得不到B。具体问题出在哪？我也不是很清楚。
当导体形状不规则时，先进行电流传导分析，然后通过LDREAD将电流密度加载后，再进行磁场分析。 感觉楼主空气单元的自由度设置的不对
应该把et,2,117,5 改成et,2,117,0就可以得到b的分布了 试过，得不到B的结果。 真的很奇怪，我看别人用117做成过的 按照版主说的方法，用69和97分别对电场和磁场分析了，但是得到的B和洛仑兹力完全不对，不知是什么原因，还请版主和大家给点提示：
/prep7
et,2,97 !air
et,1,69
mp,murx,2,1 !air
mp,murx,1,1 !conductor
mp,rsvx,1,4e-7 !bus
block,-0.1475,-0.1275,0,0.4015,-0.050,0.050
block,0.1275,0.1475,0,0.4015,-0.050,0.050
block,-0.6375,0.6375,0,0.4015,-0.400,0.400 !air
vovlap,all
numcmp,volu
vsel,s,volu,,3
vatt,2,1,2 !air properties
vsel,s,volu,,2
vsel,a,volu,,1
vatt,1,1,1 !conductor properties
alls
smart,6
mshape,1,3d
mshkey,0
vmesh,all
alls
NUMMRG,NODE,1e-8, , ,LOW
*do,i,1,2,1
vsel,s,volu,,%i%
allsel,below
nplot
*get,ymax,node,,mxloc,y,
*get,ymin,node,,mnloc,y,
nsel,r,loc,y,ymin
cp,next,volt,all
*get,nn%i%,node,,num,max
f,nn%i%,amps,24000
nsel,s,ext
nsel,r,loc,y,ymax
d,all,volt,0
*enddo
finish
esel,,mat,,1
nsle,s
/solu
solve
save
/prep7
et,1,97 !conductor
esel,,mat,,1
nsle,s
ldread,js,,,,,,,
alls
nsel,s,ext
d,all,az,0
/solu
alls
antype,static
eqslv,sparse
solve
/post1
esel,,mat,,1
nsle,s
plvect,jt,,,,vect,elem,on,0
etable,fmagx,fmag,x
etable,fmagy,fmag,y
etable,fmagz,fmag,z
ssum
[[i] 本帖最后由 maomaoface 于 2008-6-26 07:27 编辑 [/i]] 对于形状不规则的通电导体solid117可以直接加载总电流进行电流传导分析和磁场分析，不用solid69+solid97来进行计算。
命令流:
/prep7
et,1,117 !air
et,2,117,5 !conductor
mp,murx,1,1 !air
mp,murx,2,1 !conductor
mp,rsvx,2,4e-7 !bus
block,-0.1475,-0.1275,0,0.4015,-0.050,0.050 !conductor1 model
block,0.1275,0.1475,0,0.4015,-0.050,0.050 !conductor2 model
block,-0.6375,0.6375,0,0.4015,-0.400,0.400 !
vovlap,all
numcmp,volu
vsel,s,,,1,2
vatt,2,,2
allsel,all
smart,6
mshape,1,3d
mshkey,0
vmesh,all
asel,s,,,4
asel,a,,,10
DA,all,volt,0
allsel,all
asel,s,,,3
nsla,s,1
cp,1,volt,all
*get,nmax1,node,0,num,max
F,nmax1,amps,24000
asel,s,,,9
nsla,s,1
cp,2,volt,all
*get,nmax2,node,0,num,max
F,nmax2,amps,24000
allsel,all
/solu
asel,s,ext
DA,all,AZ,0,0
CNVTOL,AMPS,1,0.5e-9,2, ,
CNVTOL,CSG,1,0.5e-8,2, ,
allsel,all
solve
/post1
vsel,s,,,1
nslv,s,1
esel,all
emft
[[i] 本帖最后由 zhouyz 于 2008-6-26 09:52 编辑 [/i]] 在10.0版本中有EMFT宏命令，直接求电磁力等。
我用建立单元表的方法求力，结果一直很小，尝试了其他后处理方法还是很小，数量级是e-14，不知道是什么原因与EMFT的求解结果不一样。 [quote]原帖由 [i]maomaoface[/i] 于 2008-6-26 07:24 发表 [url=http://www.simwe.com/forum/redirect.php?goto=findpost&pid=1367867&ptid=837950][img]http://www.simwe.com/forum/images/common/back.gif[/img][/url]
按照版主说的方法，用69和97分别对电场和磁场分析了，但是得到的B和洛仑兹力完全不对，不知是什么原因，还请版主和大家给点提示：
/prep7
et,2,97 !air
et,1,69
mp,murx,2,1 !air
mp,murx,1,1 !conductor
m ... [/quote]
97单元和117单元差别很大。97单元是基于节点的矢量位法。117单元是基于单元边的方法。 感谢zhouyz的指点，以及infrom版主认真的批阅！
不过真奇怪，看zhouyz的修改，与我的命令流主要区别似乎在用他是用da命令加载的，而我是直接加载到节点上，思路看上去是一样的，结果确大不同。。。
另外，infrom版主说的97单元与117单元差别很大，但应当是都可以得到正确的答案啊，可我那个命令流得到的结果显然是不符合实际的。。。
我比较较真，还是再好好找找原因，争取实现无论117还是97，无论电流载荷还是js载荷，都可以得到正确的结果。 zhouyz使用的是ANSYS 10版本，命令流与maomaoface的几乎一致，却得到了正确的解答。我猜想是不是与ANSYS的版本有关。
用117或97都能进行磁场分析，但因为两者采用的方程不一样，单元设置和边界条件给定都可能不一样。
117的主自由度为单元边通量（用Az表示）。而97的主自由度为矢量磁位A的三个坐标分量Ax, Ay, Az。所以我才说差别很大。 infrom版主，我尝试了不同版本，发现不是版本问题，我的命令流在11.0里运行也是一样没有B，但不用命令，像zhouyz一样，手动把F载荷加载到面上就能得到结果，无论是用8.1，9.0还是11.0。 [quote]原帖由 [i]zhouyz[/i] 于 2008-6-26 09:56 发表 [url=http://www.simwe.com/forum/redirect.php?goto=findpost&pid=1367947&ptid=837950][img]http://www.simwe.com/forum/images/common/back.gif[/img][/url]
在10.0版本中有EMFT宏命令，直接求电磁力等。
我用建立单元表的方法求力，结果一直很小，尝试了其他后处理方法还是很小，数量级是e-14，不知道是什么原因与EMFT的求解结果不一样。 [/quote]
确实存在此问题，etable出来的fmag很小 etable出来的fmag，ssum出来的结果应该不会很小吧。
我用下面的命令：
VSEL, S, , , 1
ALLSEL, BELOW
etable,fmagx,fmag,x
etable,fmagy,fmag,y
etable,fmagz,fmag,z
ssum
得到的结果是：
SUM ALL THE ACTIVE ENTRIES IN THE ELEMENT TABLE
TABLE LABEL TOTAL
FMAGX 124.039
FMAGY 0.00000
FMAGZ -4.57010 infrom版主得到的结果才是正确的。谢谢！
还是前面版主讲过的的，需要把其中一个导体选出来再etable和ssum才能得到对的结果。
页:
[1]