SimWe仿真论坛's Archiver

COMSOL 2008年会圆满结束!

sunjuan 发表于 2008-7-23 16:14

问个问题 铁屋子如果完全把导体包住如何分析

现在大家的分析都是导体和空气面接触,,如果我导体完全在一个铁屋子里面,,导体没有和铁屋子接触,,那如何分析??因为现在我看了大家的分析都是导体伸出来和空气面接触,,如果导体完全被铁屋子封闭那用啥方法分析??

sunjuan 发表于 2008-7-23 16:19

比如我做的这个例子

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


导体最后还是伸出来和空气面接触了,,这样才能加平行条件,,如果导体完全封闭在铁屋子里那如何搞??且不和铁屋子接触

infrom 发表于 2008-7-23 17:54

导体应该自闭合，保证电流连续性。

perch88 发表于 2008-7-24 01:09

好像ansys可有一根直条形导体，一边加电流，一边电压为0的，外面包空气
ansoft需要导体与边界平齐

sunjuan 发表于 2008-7-24 09:43

高手帮我搞下

我要求直线导体完全被铁屋子封闭,,没有和铁屋子接触,,看屋子的磁场分布情况..我做的例子是导体伸出来和空气接触,,我现在是要封闭住导体,,如何弄它

infrom 发表于 2008-7-24 10:04

根据实际问题，导体实际走向建模，如果导体被铁屋子完全包围住的话，导体端部肯定有连接，使得电流形成自闭合回路。

sunjuan 发表于 2008-7-24 10:53

请教个问题

/PREP7  
! 定义单元
ET,1,SOLID69

! 定义材料属性
MP, RSVX, 1, 1E-7

! 建立模型
block, -0.2, -0.1, 0, 0.9, 0, 0.02  
block, -0.2, -0.1, 0.3, 0.4, 0, 0.02
vovlap, all
block, -0.2, -0.1, 0.5, 0.6, 0, 0.02
block, -0.2, -0.1, 0.7, 0.8, 0, 0.02
vovlap, all
VGEN,2,all, , ,0.3, , , ,0
block, -0.2, 0.2, 0.3, 0.4, 0.02, 0.04  
block, -0.2, -0.1, 0.3, 0.4, 0.02, 0.04
vovlap, all
block, 0.1, 0.2, 0.3, 0.4, 0.02, 0.04   
vovlap, all
FLST,3,3,6,ORDE,3   
FITEM,3,16  
FITEM,3,-17
FITEM,3,19  
VGEN,3,P51X, , , ,0.2, , ,0
vglue, all  
numcmp, all

! 剖分
MSHAPE,0,3D
MSHKEY,1
VMESH, ALL
! 令右上端面节点的电压为零
ASEL, S, LOC, Y, 0.9
ASEL, R, LOC, X, 0.1, 0.2
NSLA,S,1
d, all, volt, 0

! 耦合左下端面节点的电压
ASEL, S, LOC, Y, 0
ASEL, R, LOC, X, -0.2, 0.1
NSLA,S,1
cp, next, volt, all

! 获取左下端面节点的最小节点编号
*GET, NMIN, NODE, , NUM, MIN

! 在该节点加载电流8000 A
F,NMIN,AMPS,8000
allsel  
FINISH

/SOL
! 稳态电流传导分析
ANTYPE,0
SOLVE   
FINISH  

/POST1  
! 画源电流密度矢量图
PLVECT,JS, , , ,VECT,ELEM,ON,0

我要把我做的例子中的两个导体用现在这种形状的导体代替，，电流方向还是相反，，如何设置？/其他的都一样，，就是导体的形状变化了，，我做的是直线导体，，现在这个不是直线倒替，，有电流偶合的问题，，哪个高手有耐心搞下。。

sunjuan 发表于 2008-7-24 11:09

也就是通电导体的形状发生变化了

其他的都一样,,这个导体的形状复杂一些,,难搞一些

bftdlsyj 发表于 2008-7-24 11:29

自闭和是什么意思？难道是是短路导体吗？两边都设置电压为0，再在一端加载电流？

infrom 发表于 2008-7-24 11:41

否！

产生磁场的源是电流。而要使得磁场求解正确，加载的电流必须保证无散，即在任意点处，或任意截面上流入的电流等于流出的电流。

sunjuan 发表于 2008-7-24 13:53

我做的第一个例子导体是直线导体

block, -0.2, -0.1, 0, 0.9, 0, 0.02  
block, -0.2, -0.1, 0.3, 0.4, 0, 0.02
vovlap, all
block, -0.2, -0.1, 0.5, 0.6, 0, 0.02
block, -0.2, -0.1, 0.7, 0.8, 0, 0.02
vovlap, all
VGEN,2,all, , ,0.3, , , ,0
block, -0.2, 0.2, 0.3, 0.4, 0.02, 0.04  
block, -0.2, -0.1, 0.3, 0.4, 0.02, 0.04
vovlap, all
block, 0.1, 0.2, 0.3, 0.4, 0.02, 0.04   
vovlap, all
FLST,3,3,6,ORDE,3   
FITEM,3,16  
FITEM,3,-17
FITEM,3,19  
VGEN,3,P51X, , , ,0.2, , ,0
vglue, all  
numcmp, all

! 剖分
MSHAPE,0,3D
MSHKEY,1
VMESH, ALL
! 令右上端面节点的电压为零
ASEL, S, LOC, Y, 0.9
ASEL, R, LOC, X, 0.1, 0.2
NSLA,S,1
d, all, volt, 0

! 耦合左下端面节点的电压
ASEL, S, LOC, Y, 0
ASEL, R, LOC, X, -0.2, 0.1
NSLA,S,1
cp, next, volt, all

! 获取左下端面节点的最小节点编号
*GET, NMIN, NODE, , NUM, MIN

! 在该节点加载电流8000 A
F,NMIN,AMPS,8000
allsel  
FINISH

/SOL
! 稳态电流传导分析
ANTYPE,0
SOLVE   
FINISH  

/POST1  
! 画源电流密度矢量图
PLVECT,JS, , , ,VECT,ELEM,ON,0
这个例子中导体要偶合等技巧保证电流不发散,,我想问下如果我把我做的导体换成这个例子中的导体形状,,其他的不变化,,如何处理,,哪个有兴趣挑战下??

页: [1]
 

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