推荐使用优化设计中的sweep法效果跟循环一样
将循环值设为优化变量
自己希望获得变量为目标变量
等间隔取值
这是一个我以前用过的例子
保存为.mac
/TITLE,2DLJL
/FILNAME,2DLJL,1
KEYW,MAGNON
e=0.8e-3
/PREP7
a=0.004
N=50
I=2
f=2.176e-8/a/e
c=(f+e)/2
b=(f-e)/2
ET,1,PLANE13
EMUNIT,MKS
KEYOPT,1,3,1
MP,MURX,1,1
MP,MURX,2,1
MP,MURX,3,1.0997785406
MP,MGYY,3,8.9e+5
TB,BH,4,,40
TBPT,,159.2,0.2402
,,318.3,0.8654
,,477.5,1.1106
,,636.6,1.2458
,,795.8,1.331
,,1591.5,1.5
,,3183.1,1.6
,,4774.6,1.683
,,6366.2,1.741
,,7957.7,1.78
,,15915.5,1.905
,,31831,2.025
,,47746.5,2.085
,,63662,2.13
,,79577.5,2.165
,,159155,2.28
,,318310,2.485
,,397887,2.585
FINISH
/PREP7
/PNUM,AREA,1
RECTNG,0,0.0019,0.0153,0.0243
RECTNG,0.002,b,0.0153,0.0243
RECTNG,b,c,0.0153,0.0153+a
RECTNG,b,c,0.0153+a,0.0243
RECTNG,c,0.00525,0.0153,0.0163
RECTNG,c,0.00525,0.0163,0.0233
RECTNG,c,0.005,0.0233,0.0243
RECTNG,0.00525,0.006,0.0223,0.0233
ASEL,S,AREA,,ALL
ARSYM,Y,ALL
RECTNG,0,0.0019,-0.015,0.015
RECTNG,0.002,b,-0.015,0.015
RECTNG,b,0.00525,-0.014,0.014
RECTNG,c,0.00525,0.0141,0.015
RECTNG,c,0.00525,-0.015,-0.0141
RECTNG,0,0.01,-0.05,0.05
AOVLAP,ALL
AGLUE,ALL
/REPLOT
asel,all
ASEL,S,AREA,,30
ASEL,A,AREA,,33
AATT,2,1,1,0
asel,all
ASEL,S,AREA,,20,21
ASEL,A,AREA,,23
ASEL,A,AREA,,26
AATT,3,1,1,0
asel,all
asel,u,,,30,33,3
asel,u,,,20,21
asel,u,,,23,26,3
asel,u,,,38
AATT,4,1,1,0
/PNUM,MAT,1
ALLSEL,ALL
SMRTSIZE,2
MSHAPE,0,2D
MSHKEY,0
AMESH,ALL
ASEL,S,AREA,,19,35,16
ESLA,S
CM,ARM,ELEM
FMAGBC,'ARM'
FINISH
/SOLU
ALLSEL,ALL
NSEL,EXT
D,ALL,AZ,0
FINISH
/SOLU
ASEL,S,AREA,,30
ESLA,S
BFE,ALL,JS,,,,-N*I/(a*e)
ASEL,S,AREA,,33
ESLA,S
BFE,ALL,JS,,,,N*I/(a*e)
ALLSEL,ALL
SOLVE
FINISH
/POST1
ETABLE,FMAGY,FMAG,Y
SSUM
*GET,FXL,SSUM,,ITEM,FMAGY
FXLA=10-FXL
FINISH |