这个标题起大了，其实只是跟我们的第一次作业有关而已

刚上手，看到一堆.exe，肯定一头雾水，这个时候我们地说明书和例子就起到作用了

doc文件夹里面是说明书 助教讲过看法了对吧
在这次作业中，对我们有用的是SFFILES和SFCODES
例子的话用的是二级磁铁\Magnetostatic\Dipoles

好了，开始着手做了：
首先打开说明书SFFILES，目录里面写着Automesh input file 在126
于是翻到126P
说明书上说过程分为三步：REG、PO、MT
REG：for region 也就是分区，是为了建立物理坐标与逻辑坐标的关系
其实吧，现在的电脑比原来强多了，所以直接一个逻辑分区就够了，逻辑坐标大概在物理坐标的4到5倍就可以保证结果比较精确，当然，不介意更密
&reg kprob=0,
mode=0,
xreg1=17.7,kreg1=90
kmax=100,
yreg1=46,lreg1=230,
lmax=250 &
以上是直接一个分区的写法


PO：for point 我还是喜欢叫它描边
上边的分区是为了确定计算精度，这里是来描绘磁铁与线圈的形状了。
&po x= 0.0000,y= 0.0000 &
&po x=57.1500,y= 0.0000 &
&po x=57.1500,y=32.4000 &
&po x= 0.0000,y=32.4000 &
&po x= 0.0000,y= 0.0000 &
例子上来的这段，表示先确定计算磁场的空间大小
接下来
&reg mat=1,cur=-51000.0 &
&po x=15.6800,y= 7.3500 &
&po x=39.3800,y= 7.3500 &
&po x=39.3800,y=15.0600 &
&po x=15.6800,y=15.0600 &
&po x=15.6800,y= 7.3500 &
&reg 表示这个描绘的某种物质的边线，而不是像之前的要计算的磁场的空间的大小
之后又出现了各种参数，所以这个时候就要打开SFCODES，打到147P，会看见一个很长很长的表，里面有这些参数的意义。
mat Material number 就是对物质做一个标记
cur 就是总电流 单位安培 顺带一提物理坐标的单位是cm
下面的那些是描边，方法助教有讲过


&reg mat=2,mtid=1,mshape=0 &
&po x= 0.0000,y= 5.2400 &
&po x= 2.0232,y= 5.2400 &
&po x= 2.4279,y= 5.2401 &
&po x= 2.8325,y= 5.2403 &
&po x= 3.2371,y= 5.2407 &
&po x= 3.3618,y= 5.2412 &
&po x= 4.0464,y= 5.2410 &
&po x= 4.4511,y= 5.2392 &
&po x= 4.8557,y= 5.2359 &
&po x= 5.2604,y= 5.2309 &
&po x= 5.6650,y= 5.2307 &
&po x= 6.0696,y= 5.2416 &
&po x= 6.4743,y= 5.2636 &
&po x= 6.8789,y= 5.2968 &
&po x= 7.2836,y= 5.3144 &
&po x= 7.6882,y= 5.2899 &
&po x= 8.0929,y= 5.2233 &
&po x= 8.4975,y= 5.1145 &
&po x= 8.9021,y= 5.0129 &
&po x= 9.3068,y= 4.9679 &
&po x= 9.7114,y= 4.9793 &
&po x=10.1161,y= 4.9850 &
&po x=10.5207,y= 5.0350 &
&po x=10.9254,y= 5.1000 &
&po x=11.3300,y= 5.2400 &
&po nt=2,x0=10.5,y0=9.3,x=2.7067,y=-3.1424 &
&po x=14.7800,y= 7.2500 &
&po x=14.9900,y=14.0600 &
&po nt=2,x0=16.26,y0=14.06,x=0.0,y=1.27 &
&po x=38.8000,y=15.3300 &
&po nt=2,x0=38.8,y0=14.06,x=1.27,y=0.0 &
&po x=40.0700,y= 0.0000 &
&po x=57.1500,y= 0.0000 &
&po x=57.1500,y=32.4000 &
&po x= 0.0000,y=32.4000 &
&po x= 0.0000,y= 5.2400 &
这个是对磁铁的描边。。。很长，不过我们的硅钢形状很好，不需要这么长。。。


里面的mtid 说明书里说的是 指向材料的指针
其实就是将材料与MT步骤中输入的材料性质对应起来
比如，这里面是
&reg mat=2,mtid=1,mshape=0 &
表示，将mtid=1的材料的属性赋予 mat=2的材料，其中mshape=0的意义说明书说明很奇怪，没弄懂，不过当成缺省值来理解就可以了

MT：for material table
就是开始输入材料的性质了。。。
按照说明书的说法，性质有三种输入方法：
BGAM BMU BH 具体可以参见SFCODES的158P的表
例题中使用的是BGAM 老师给我们的参数是BH的关系，所以为了省事，我们就用BH输入
B的单位是高斯 H的单位是奥斯特
转换关系为
1 T = 1*10^4 G 1 A/M = 4π×10^(-3) Oe


&mt mtid=1
bh=3000 0.0816814076
4000 0.09173450392
5000 0.10178760024
6000 0.10932742248
7000 0.11686724472
8000 0.12817697808
9000 0.13948671144
10000 0.15205308184
11000 0.17090263744
12000 0.20357520048
13000 0.2450442228
14000 0.3078760748
15000 0.4146902232
16000 0.6534512608
17000 1.633628152
18000 6.2831852
19000 28.90265192
20000 124.40706696 &
这就是转换过单位制的BH关系


全部输完后，另存为.am文件
双击运行（或者用automesh.exe运行）
生成.T35文件，右键→poisson 弹出运行框，等运行框自动关闭
再次打开.T35文件（双击或者用wsfplot.exe运行）得到结果

嗯，就这么多，还有疑惑的话，可以自己翻说明书。
位置之前也都说过了。。。
其实这次作业很简单对吧

mshape=0按照助教的意思是材料各向同性。。。
其实我是来用小号刷经验的，呵呵

赞！

好操作