SimWe仿真论坛's Archiver

COMSOL 2008年会圆满结束!

zhuzhouren 发表于 2008-8-4 16:49

误差分析

本意是要分析两个导体的磁场分布,,外面带了个封闭的铁筐子。。结果如附件所示。。命令流如下。。
FINISH                       
/CLEAR,START                  
/units,si                  
   
/TITLE, Test for EDGE-Formulation with SOLID117                                             
!   
! Parameters

*SET,a,1e-6 ! 1st addition length of magnetic core  
                                                        
! *** Create Model                                                        
!   
FINISH  
/PREP7  
  
ET,1,117 ! Element type #1 is a magnetic edge element, 117.
MP,MURX,1,1 ! Define relative permeability of material
MP,MURX,2,1 ! Define relative permeability of material
MP,MURX,3,500 ! Define relative permeability of material
MP,MURX,4,1 ! Define relative permeability of material
BLOCK,-(480*a),-(460*a),(63*a),(163*a),-(1000*a),(1000*a)
BLOCK,(60*a),(80*a),(63*a),(163*a),-(1000*a),(1000*a)
BLOCK,-(750*a),(750*a),(438*a),(513*a),-(100*a),(100*a)
BLOCK,-(750*a),(750*a),-(90*a),0,-(100*a),(100*a)
BLOCK,-(750*a),-(700*a),0,(438*a),-(100*a),(100*a)  
BLOCK,(700*a),(750*a),0,(438*a),-(100*a),(100*a)
BLOCK,-(2000*a),(2000*a),-(1000*a),(2000*a),-(1000*a),(1000*a)

! Overlap volumes   
VOVLAP,ALL  
NUMCMP,VOLU
! Assign attributes to the volumes  
VSEL,S,VOLU,,1  
VATT,1 ! conductor
VSEL,S,VOLU,,2
VATT,2

VSEL,S,VOLU,,3,6,1 ! magnetic core
VATT,3
VSEL,S,VOLU,,7! air ambient   
VATT,4  
! Mesh volumes

点击MESH ,   然后点击VOLUMNS  ,    然后点击FREE   ,   划分网格

   
ALLSEL,ALL  
FINISH  
/SOLU   
!*** Apply current density on conductor ***
ESEL,S,MAT,,1 ! Select coil element
BFE,ALL,JS,1,,,-1e9 ! Apply current density
ESEL,S,MAT,,2 ! Select coil element
BFE,ALL,JS,2,,,1e9 ! Apply current density

!*** Flux parallel boundary condition at all extern area of air ambient ***
ASEL,S,EXT  
DA,ALL,AZ,0
ALLSEL,ALL  
! Solution parameters   
ANTYPE,0
NROPT,AUTO,,
EQSLV,SPAR,,0,  
PRECISION,0
MSAVE,0
PIVCHECK,1  
/STATUS,SOLU
SOLVE   
FINISH  
/POST1  
!*** Vectorplot of magnetic flux (BSUM) for magnetic core's and conductor's
! elements  
/VSCALE,1,1,0   
!   
!*  
PLVECT,B, , , ,VECT,ELEM,ON,0   
ESEL,U,MAT,,4
/REPLOT
PLVECT,B, , , ,VECT,ELEM,ON,0   
ESEL,U,MAT,,2   
/REPLOT
/EFACET,1   
PLNSOL, B,SUM, 0

修改的部分:修改的部分:修改的部分:修改的部分:修改的部分:

高手将
ESEL,S,MAT,,1 ! Select coil element
BFE,ALL,JS,1,,,-1e9 ! Apply current density
ESEL,S,MAT,,2 ! Select coil element
BFE,ALL,JS,2,,,1e9 ! Apply current density

改成

ESEL,S,MAT,,1 ! Select coil element
BFE,ALL,JS, ,,,-1e9 ! Apply current density
ESEL,S,MAT,,2 ! Select coil element
BFE,ALL,JS, ,,,1e9 ! Apply current density

,,我想请问是不是代表我做的就有错误??是错误还是误差,,如果是误差可以被忽越不??原因

llkg 发表于 2008-8-4 17:52

俺觉得如果你针对这个问题做了一些尝试,可以把你尝试的结果更新给大家看看。
你还是停留在以前发帖子的状态吗?我不相信。呵呵
你的这个宏命令,我运行不了。你把mesh那部分手动操作后的log文件中的宏命令拷贝过来就可以把这个宏命令写完整。

sunjuan 发表于 2008-8-5 09:50

(点击MESH , 然后点击VOLUMNS , 然后点击FREE , 划分网格)

(点击MESH ,   然后点击VOLUMNS  ,    然后点击FREE   ,   划分网格)

这一步我也不是用的命令流的形式,,我用的是点左边的菜单的形式,,点MESHING>MESH>VOLUMNS>FREE,,,,不是用的命令流的形式...


其余的过程都是用的命令流的形式。。就是这个地方用的是点左边菜单的形式

llkg 发表于 2008-8-5 10:49

对啊,没有问题,你就这样手动操作,
ansys中list->files->log file中会记录这部分的操作。那就是这部分对应的宏命令。你把这部分宏命令拷贝到你的大文件中,以后剖分就不用再手点了。

ni和zhuzhouren是一个人吗?呵呵。猜测。如果错了,不要见怪。

zhuzhouren 发表于 2008-8-5 14:57

是不是不能运行我的命令流??

我已经给出来了详细的方式,,MESH这个部分也照我的样子搞,,哪个高手跑下看下问题所在

zhuzhouren 发表于 2008-8-5 14:58

到底是误差还是错误

到底是误差还是错误

sunjuan 发表于 2008-8-6 11:31

哪个高手帮忙搞下

告诉我是误差还是错误

页: [1]
 

Powered by Discuz! Archiver 6.1.0  © 2001-2007 Comsenz Inc.