【急】求助ANSYS节点的磁场强度计算

billneu

【急】求助ANSYS节点的磁场强度计算


小弟是ANSYS新手中的新手，遇到一个很恼人的问题，望各位前辈指点～～

目的：希望分析获得：两块永磁体中间的空气间隙中的任意一个层面中的所有节点的磁场强度H

即：


再包一层球形空气层：





编写命令流如下（暂时任选了LOC,Y,0.0017这个面）：

FINISH
/CLEAR
/TITLE,PERMANENT MAGNET WITH sphere AIR surround !定义工作标题
/FILNAME, Bill_3D_block,1 !定义工作文件名
KEYW,MAGNOD,1 !指定磁场分析
/PREP7
SMRT,OFF
ET,1,SOLID98,10 !四面体耦合场单元
EMUNIT,MKS !MSK单位制

MP,MURX,1,1
MP,MURX,2,1
MP,MURX,3,8.03 !永磁

MP,MGXX,3,0            
MP,MGYY,3,405000   
MP,MGZZ,3,0     

BLOCK,-12E-3,12E-3,14E-3,5E-3,-12E-3,12E-3
BLOCK,-12E-3,12E-3,-17E-3,-5E-3,-12E-3,12E-3
BLOCK,-12E-3,12E-3,5E-3,-5E-3,-12E-3,12E-3
SPHERE,0.1

VSEL,ALL
VOVLAP,ALL   
NUMCMP,ALL   

VSEL,S,,,3
VATT,1  !设置内空气属性  
VSEL,S,,,4
VATT,2  !设置外空气属性
VSEL,S,,,1,2
VATT,3  !设置永磁属性
VSEL,S,,,1,4
ESIZE,,6

MSHK,0    !自由体划分
MSHA,1,3D   !使用四面体单元
VMESH,ALL
CSYS,0   
ESEL,ALL   
NSEL,ALL   
FINISH

/SOLU
MAGSOLV,2  !RSP方法求解
FINISH

/POST1
ESEL,,MAT,,1  
NSLE,S,ALL
NSEL,R,LOC,Y,0.0017
NPLOT,1
NLIST,ALL
PRNSOL,H,COMP  !列出磁场强度

－－－－－－－－－－－－－－－－－－
...
NSEL,R,LOC,Y,0.0017   !选择出需要的层面
NPLOT,1                       !画出所选层面的所有节点，截图如下：



NLIST,ALL              !获得了被选点的坐标，截图如下：


PRNSOL,H,COMP        !但是再继续通过此命令希望获得这些点上的磁场强度，结果却为空，即没有任何一个节点的数据；换试选择其他层面上的点，得到节点数有所不同，但磁场强度的结果均为空。截图如下：



为什么为什么啊！！！  H都到哪里去了啊！！！


是命令用错了？网格划得不合适？还是ANSYS就是算不出磁场强度来？但是节点都选出来了啊！

这个苦恼烦躁啊，因为急需正确数据来进行下一步的重要工作！
因此急盼各位前辈高手达人们伸出金手指点指点我错在哪儿？
要实现目的应该如何操作？
多谢多谢！！！

billneu

   有人帮帮忙嘛～～～

elvint

不太清楚，首先运行了你的程序，竟然没有计算结果，那么自然就不会显示任何数据了。
可能是磁性材料赋值出现了问题，还有，定义永磁材料时怎么没有用到ESYS呢，不知道和这个有没有关，你再试试吧。

land908

学习一哈，来看看啊

bouquet

/SOLU
MAGSOLV,2  !RSP方法求解
FINISH
--------------------------------------
这样就算是求解了？起码还要solv吧

luomu

你所选的点并不在所选的单元上
我计算后可以取出节点的H值
但没有你所列出来的那些点

billneu

不太清楚，首先运行了你的程序，竟然没有计算结果，那么自然就不会显示任何数据了。
可能是磁性材料赋值出现了问题，还有，定义永磁材料时怎么没有用到ESYS呢，不知道和这个有没有关，你再试试吧。
elvint 发表于 2009-12-22 21:47

没有计算结果吗？怎样才算有计算结果呢？
...
/SOLU
MAGSOLV,2  !RSP方法求解
FINISH

!在这个求解之后（之前命令不变），如果后处理的时候用以下命令：

/POST1
ESEL,,MAT,,1  !选定中间的空气层
/VSCALE,,,1 !设置统一的矢量缩放因子
PLVECT,H !绘出磁场强度矢量

那么可以得到这个图：



我觉得没有计算结果的话后处理的时候应该也显示不了结果呀？







ESYS我没用过，因为我是照着书上的3D永磁体的例子改编的，它里头没有用ESYS；
另外我查了一下help，ESYS应该是和坐标有关的吧，如果没有用是不是应该也不会影响永磁材料的属性？

billneu

/SOLU
MAGSOLV,2  !RSP方法求解
FINISH
--------------------------------------
这样就算是求解了？起码还要solv吧
bouquet 发表于 2009-12-24 14:26

MAGSOLV不是求解吗？我真的不懂啊～
如果没有求解应该  上面 7#的 图也画不出来的吧？
solv是什么命令？ 如果还要像你说的加solv那应该怎么加？能麻烦你帖出来命令吗？

billneu

你所选的点并不在所选的单元上
我计算后可以取出节点的H值
但没有你所列出来的那些点
luomu 发表于 2009-12-25 17:15

下面这个命令：

ESEL,,MAT,,1    !选定中间的空气层
NSLE,S,ALL    !NSLE这个命令help里的解释是：Selects those nodes attached to the selected elements
NSEL,R,LOC,Y,0.0017    !NSEL这个命令help里的解释是：Selects a subset of nodes.

!这三条命令都是一环扣一环，每一条都是在上一条所选的范围里继续选出来的吧？最后出来的面上的点应该是在所选单元上啊？


不知道你说的“我计算后可以取出节点的H值”是选出了指定平面上的节点的H值吗？
能麻烦帖出你的命令不？谢谢了啊～～～







另外， 我试着不选定面，在选定中间的空气层这个体后，就让列出节点和H值的话：
即：
...
/SOLU
MAGSOLV,2  !RSP方法求解
FINISH

!在这个求解之后（之前命令不变），如果后处理的时候用以下命令：

/POST1
ESEL,,MAT,,1  !选定中间的空气层体
NSLE,S,ALL  !选定所选体中的所有节点
NLIST,ALL  ! 获得被选节点的坐标
PRNSOL,H,COMP  !列出这些节点的磁场强度


得到的结果是的确能得到一些H的值，但是仔细看发现好多节点在H的输出里都不见了：
（左边Nlist输出的节点信息中红框里的节点在PRNSOL输出的节点H里都找不到）



把两个输出都拉到底，又可以发现右边PRNSOL输出的节点H里，节点编号仅到3199，而不存在左边Nlist输出的尾部大量节点：








疑惑啊！！！苦恼啊！！！

lilulilu

allsel了没？
如果是找不到H值，先看看云图也行啊。
如果有云图，就说明没有计算错误，那就换个节点看看有没有结果。
如果别的节点有结果，那就把选择的节点范围扩大点，看看在你要的位置附近有没有值。

billneu

allsel了没？
如果是找不到H值，先看看云图也行啊。
如果有云图，就说明没有计算错误，那就换个节点看看有没有结果。
如果别的节点有结果，那就把选择的节点范围扩大点，看看在你要的位置附近有没有值。
lilulilu 发表于 2009-12-26 17:26

谢谢版主回复～


你指的“ALLSEL,ALL”是放在“/SOLU”的后面吗？
好像没有什么影响，我试了，结果是和没有使用一样的



两块永磁铁中间的空气夹层云图如下：



是不是这就意味着“没有计算错误”呢？ 那难道说ANSYS就是没法算出每一个节点上的H？

现在感觉是有N多节点没有结果啊,  这个范围可能会扩得很大？我需要的是一个界面里的节点而且需要尽可能的多来做一个H值的二维平面矩阵...
另外，弱弱问一下，怎么扩大，就是怎么选择指定点附近最近的节点？难道一个一个试？－－－真的很新手中的新手...

lilulilu

你这个网格太粗了，明显很多位置没有分布节点。把网格划细一些。

billneu

你这个网格太粗了，明显很多位置没有分布节点。把网格划细一些。
lilulilu 发表于 2009-12-28 10:47

各位版友，版主新年快乐

liuliu版主， 我想问的是，是不是说ANSYS就是没法算出每一个节点上的H值？

同一个节点，不会因为网格画得粗就没有值，画得细就有值吧？

而且我还真画细了，虽然直接导致运行的极其耗时，但是问题依旧啊～





另外，你所说的：“如果别的节点有结果，那就把选择的节点范围扩大点，看看在你要的位置附近有没有值。lilulilu 发表于 2009-12-26 17:26 ”
我想知道你有没有办法选出指定点附近/最近的节点？（非人工肉眼挑选方法）





再一个，你知道如何画出指定切面的云图吗？





麻烦不吝赐教，详细些哈～～




其他各路高手都请不吝献计献策～～

perch88

你的问题关键是不了解在中间节点上是没有结果的
nlist列出所有节点，而结果中只没有中点上的结果
如果你想在nlist和结果中看到相同的节点，就应该使用一次单元，比如13；不要使用二次单元（有中点的单元），如53.

再一个，你知道如何画出指定切面的云图吗？
——这个我写过：http://blog.sina.com.cn/s/blog_61e2cf840100fj80.html

billneu

谢谢版主回帖！！！

你的问题关键是不了解在中间节点上是没有结果的...
perch88 发表于 2010-1-4 02:32

意思就是说，有些节点没有结果是必定的是吗？

如果你想在nlist和结果中看到相同的节点，就应该使用一次单元，比如13；不要使用二次单元（有中点的单元），如53....
perch88 发表于 2010-1-4 02:32

意思就是说，还是有方法使得所有节点上都有结果是吗？就是我的命令流其他地方没有什么错误，只要改成使用合适的单元就能获得所有节点上的结果了是吗？

那么我这里的情况应该是用什么样的单元呢？能再详细些明示吗？小弟真的极其初级，就更搞不懂“一次”“二次”单元选择上的玄机了～因为我这里用的是三维的，你举的例子是二维里的实体单元13、15，三维里又是什么呢？ 我先前仅仅是按照书上的例子编写的命令流，书中用的是“SOLID98”我也不懂，就照着用了。

P.S.查找帮助里的PLANE13时，有一句话说“具有类似功能的耦合场单元是SOLID5、SOLID62和SOLID98”不知道是不是就意味着“SOLID98”已经是你所说的“一次单元”了，就应该能“在nlist和结果中看到相同的节点”了呢？

再一个，你知道如何画出指定切面的云图吗？
——这个我写过：http://blog.sina.com.cn/s/blog_61e2cf840100fj80.html...
perch88 发表于 2010-1-4 02:32

是的，我查到你这个帖子了，我回复在那个帖子里了（详见：http://forum.simwe.com/thread-741204-1-1.html）：

...我有此帖一样的问题，按照你说的去做我没有结果，请你帮忙看看问题出在哪儿： ...  
...
/POST1
ESEL,,MAT,,1        !选出中间空气间隙体
NSLE,S,ALL
PLNSOL,H,SUM       !整个体的磁感应强度云图（是这么用命令没错吧？）
NSEL,R,LOC,Y,1.7E-3          !选择其中的一个平面
然后用这个方法：

plotctrl->hidden lines & options
z-buffer改成Q slice z-buffer
powergraphics改成full
/replot

可是没有显示...
billneu 发表于 2010-1-3 10:11

gaohong004

你好。我做的电场二维分析，也出现了你这样的结果。在单元边中节点的位置结果为0。请问你最后怎么解决的啊，采用降阶单元吗？

shagf

此帖讨论的内容很好，都了很受益。lz的问题就是节点出了毛病，应该是网格划分不当导致的

shagf

我做了一个二维磁场的分析，也是有这个毛病啊

花盆

马一下，学习之
也要遇到磁场强度的问题