通电线圈+U形铁芯.单元98+36标量DSP求解器出现警告

pipinia

各位好，我仿照VM190做了一个仿真，仿真的模型为一个缠绕有线圈的U形铁芯，下方为一个铁试块，周围全是空气。我想观察铁芯末端和试块之间的B，但是计算时会有No constraints have been defined using the D command和Coefficient ratio exceeds 1.0e8 - Check results.这样的警告。是不是与边界条件有关？？
       还要用D来增加什么边界条件呢？这是我第一次做3d的磁场仿真，错误之处，还请各位帮忙指出，谢谢，谢谢了！

/PREP7
SMRT,OFF

ET,1,SOLID98,10               ! 10-NODE TETRAHEDRAL, MAG OPTION
ET,3,SOURC36                  ! CURRENT SOURCE ELEMENT
EMUNIT,MKS                    ! MKS UNITS
!定义各个参数
I=1                           ! CURRENT(AMPS)
N=100                         ! TURNS
COIL_L=0.02                   ! COIL LENGTH
COIL_R=0.008                  ! COIL RADIUS
COIL_T=0.0005                 ! COIL THICKNESS
IRON_L=0.06                   ! IRON LENGTH
IRON_W=0.01                   ! IRON WIDTH
IRON_H=0.03                   ! IRON HEIGHT
SPE_L=0.08                    ! SPECIMEN LENGTH
SPE_W=0.03                    ! SPECIMEN WIDTH
SPE_H=0.03                    ! SPECIMEN HEIGHT
LIFTOFF=0.005                  ! LIFTOFF BETWEEN IRON AND SPECIMEN
GAP=0.01
R,2,1,N*I,COIL_T,COIL_L       ! COIL CONSTANT SET 1:COIL TYPE, CURRENT, THICKNESS, LENGTH
！材料属性
MP,MURX,1,1                   ! RELATIVE PERMEABILITY OF AIR
MP,MURX,2,4000                ! RELATIVE PERMEABILITY OF IRON
mp,murx,3,4000                ! relative permeability of specimen
！生成sourc36线圈
N,1,COIL_R,LIFTOFF+IRON_H-IRON_W/2,0      ! CREATE NODES TO LOCATE COIL
N,2,0,LIFTOFF+IRON_H-IRON_W/2+COIL_R,0
N,3,0,LIFTOFF+IRON_H-IRON_W/2,0
TYPE,3                        ! SOURCE ELEMENT TYPE
REAL,2
E,1,2,3                       ! DEFINE COIL ELEMENT

!创建铁芯
BLOCK,-IRON_W/2,IRON_W/2,LIFTOFF,LIFTOFF+IRON_H,IRON_L/2,IRON_L/2-IRON_W  
BLOCK,-IRON_W/2,IRON_W/2,LIFTOFF+IRON_H-IRON_W,LIFTOFF+IRON_H,IRON_L/2-IRON_W,-(IRON_L/2-IRON_W)
BLOCK,-IRON_W/2,IRON_W/2,LIFTOFF,LIFTOFF+IRON_H,-(IRON_L/2-IRON_W),-IRON_L/2
!创建铁试块
BLOCK,-SPE_W/2,SPE_W/2,-SPE_H,0,SPE_L/2,-SPE_L/2
!空气
BLOCK,-SPE_W/2-GAP,SPE_W/2+GAP,-SPE_H-GAP,LIFTOFF+IRON_H+GAP,-SPE_L/2-GAP,SPE_L/2+GAP
VOVLAP,ALL                    ! OVERLAP AIR, IRON AND SPECIMEN      
numcmp,volu   
           
vsel,s,volu,,2,4,1  !iron
vatt,2,1,1
vsel,s,volu,,1    !specimen
vatt,3,1,1
vsel,s,volu,,5    !air
vatt,1,1,1
vsel,all
                     
vsel,s,mat,,2        !MESH IRON
esize,0.003
mshk,0
msha,1,3d
VMESH,all                  
                        
vsel,s,mat,,3         !MESH SPECIMEN
esize,0.003
mshk,0
msha,1,3d
VMESH,all
vsel,s,mat,,1
esize,0.003
mshk,0
msha,1,3d
VMESH,all                       !MESH AIR

nsel,all
FINISH
/SOLU
ANTYPE,STATIC
MAGSOLV,3
FINISH

[ 本帖最后由 perch88 于 2009-1-6 15:25 编辑 ]

perch88

刚开始玩别用远场

pipinia

谢谢版主
后来我将远场去掉，将模型用空气包围，依然会有No constraints have been defined using the D command的警告和Large negative pivot value(...)in Eqn. system.May be because fo a bad temperature-dependent .....这种错误，据说是由于约束不够引起。
警告提示没有用D命令定义约束，这里我还需要用D命令增加什么constraints呢？磁力线平行的边界条件应该是不用说明，自然发生的吧。恳请斑竹和各位大虾帮忙指出错误的原因。
新年到了，在此也祝斑竹和各位论坛会员新年快乐，祝论坛越办越红火！：）

perch88

我没有运行你的程序，但是从你的code来看：
1 你好像什么边界条件都没有加，在最外层加平行即可
2 从你手画的图来看，边界太小，应该是模型的3倍或者以上（如果机器性能允许的话）

pipinia

谢谢斑竹。
我没有加平行的原因是：ansys不是默认磁力线平行边界条件的吗？
后来我用D命令加上了MAG=0，即磁力线垂直边界条件（测试用），仅出现了Coefficient ratio exceeds 1.0e8 - Check results，现在还不明白这个警告的产生原因。是不是code里有别的问题呢？

斑竹那么晚了还工作？呵呵 好榜样！

perch88

如果学例子，就先尽量不要修改
vm190求解用的是MAGSOLV,4
你用MAGSOLV,3
改回来看看吧，我看着好像没有警告，d那个警告可以忽略

pipinia

谢谢perch88斑竹
我用MAGSOLV,3而非4的原因是，帮助文档5.2节指出The DSP method is available for models with "singly connected" iron regions, and the GSP method is for models with a "multiply connected" iron region. VM190中的是一个multiply connected iron region，而我这里的应该是"singly connected" iron regions。不知道我的理解正不正确。
我用4也试过，从后处理中B的分布来看，不是我所期望的结果。

perch88

我试了一下结果是不对，从帮助也能看出来了，含有iron的情况是需要特别注意的。因为诸多限制，尽量在有iron时转为矢量法或棱边法计算。

另外，98不支持六面体网格，97也有计算iron的局限性，需要特别处理

结论：你的模型不大，可以试一下117

perch88

请看附件，我怀疑是你的问题。
vm190.txt是例题，我把没用的东西去掉，我修改了以后：
1 在后处理中直接画出场图
2 转换了显示场图的角度。见图1.
3 求解用的是MAGSOLV,4

vm190c.txt是例题，我把没用的东西去掉，我修改了以后：
1 在后处理中直接画出场图
2 转换了显示场图的角度。
3 我把iron断开，你可以看我在txt文件中的注释部分。见图2
4 求解用的是MAGSOLV,3

我认为两个结果都是合理的。
图3是跟图1和图2相同视角的单元图，这个角度看不到iron（紫色）。
图4是旋转以后的单元图（是vm190c.txt的单元图），可以看到断开的iron（紫色）。
--------------------------------------------------------------------------------------------------------
问答题：怎样才能看到Source36单元构成的线圈？
http://e.simwe.com/34599/viewspace-12591

问答题：ansys的标量法求解器RSP、DSP、GSP该如何选择？
http://e.simwe.com/34599/viewspace-12592
--------------------------------------------------------------------------------------------------------

[ 本帖最后由 perch88 于 2009-1-5 02:42 编辑 ]

perch88

这里是源代码：

pipinia

非常感谢斑竹的细心指教，谢谢了！
我再根据你的提示好好找找之前的问题。

abacqus

学习了