求助：2D静态磁场分析得出的结果觉得不靠谱，请大家来鉴定

繁星铺满天

分析对象是一个关于转向轴轴线对称的三维结构，因此用一个二维模型代替，采用的分析模式为平面2D静态磁场分析。据此建立了平面几何模型。建模方法采用参数化建模方法，具体参数设置如下：

r1=0.096;

r3=0.066;

r2=r3+coilb;

r4=0.031;

l2=0.012;

l3=0.006;

h=0.001；

l1=2*l2+2*h+l3；

coila=0.012；

coilb=0.017;

ringa=0.015；

ringb=0.003；

diskth=0.006；

ds=0.035；

各材料的特性如下所示：
1,2,6号材料特性用相对磁导率表示，均为1.5，其余各材料的特性用B-H曲线表示，数据如下表3号磁化曲线数据

H(kA/m)	0.1	0.4	0.8	1.2	1.5	1.8	1.9	2	2.1
B(T)	0.145	1.21	1.5027	1.5523	1.5659	1.5761	1.5793	1.5824	1.586

H(kA/m)	2.2	2.5	3	3.5	4	5	10	30	50
B(T)	1.589	1.598	1.6137	1.6295	1.6455	1.678	1.8497	2.732	4.0351

4号液磁化曲线数据

H(kA/m)	30	60	90	120	180	240	300
B(T)	0.24	0.44	0.568	0.665	0.787	0.91	1.013

5号磁化曲线数据

H(kA/m)	0.1	0.5	1	1.5	2	2.5	3	3.5	4
B(T)	0.445	0.6509	0.8816	1.0841	1.26	1.411	1.5388	1.645	1.7312

H(kA/m)	4.5	4.8	5.1	5.5	5.8	6.1	6.5	7
B(T)	1.7991	1.8318	1.8588	1.8867	1.9019	1.9127	1.9208	1.922

然后通过矩形（rectangle-create by dimension）两对角点坐标建立了6个矩形，依次输入下列各点坐标（每一行是两个对角点坐标）：

0,0      l1,ds0,0      l1,r1l2,0     l2+diskth+2*h,r3l2+h,0   l2+h+diskth,r3-0.001l2+h+(diskth-ringa)/2,r3       l2+h+(diskth+ringa)/2,r3+ringbl2+h+(diskth-coila)/2,r3+ringb  l2+h+(diskth+coila)/2,r3+ringb+coilb

对这六个矩形进行overlap运算得到如下结果：

再对最下面一层五个举行进行add（加）运算，得到如下结果：

然后通过meshtool对各个面赋予材料特性，各个面赋予特性如下所示：

面	A2	A1	深红色面	深青色U行面	绿色面	紫红色面
材料号	1	2	3	4	5	6

紧接着用meshtool划分网格，采用智能划分方式，把smart size值调到1，划分结果如下：

然后选取A2面所有元素施加电流载荷，大小为入7816000A/m2

然后选取所有外围节点，施加磁力线平行条件，选择所有试题进行运算，最终得到这样的结果

我觉得结果有问题是因为：
1.这是一个对称结构，为什么磁感应强度分布不对称？并且磁场线的分布也比较怪异。
2.如果是我操作错误，到底是哪里错了？

希望各位高手多多指点，告诉我怎么样正确处理这个问题。

繁星铺满天

更正：建模时输入的点坐标上面没写清楚，更改如下：

0,0      l1,ds
0,0      l1,r1
l2,0     l2+diskth+2*h,r3
l2+h,0   l2+h+diskth,r3-0.001
l2+h+(diskth-ringa)/2,r3       l2+h+(diskth+ringa)/2,r3+ringb
l2+h+(diskth-coila)/2,r3+ringb  l2+h+(diskth+coila)/2,r3+ringb+coilb

kamakot

一般要有完整的命令流才能发现你的错误之处，你这样贴别人是不易查找你的错误的。

6xingjidexing

首先，这个对称是指以Y轴为对称轴的中心对称，所以也没什么不对称的。对于3号材料，我不知道你用的是什么材料的BH，为什么最后那两个值上升的那么快呢？