[转帖]瞬态分析实例

meil32

瞬态分析中，边界条件和载荷都是时间的函数。每个时间段的加载求解称为一个载荷步。可能是阶跃、斜坡或是自动时间阶跃三种。瞬态磁场分析的载荷与其他分析类似，载荷类型有时间积分电势、电压、电流等。当载荷为电压降时，仅能使用具有矢量磁势和电流自由度的PLANE53单元进行分析。而节点电流载荷仅用在载有外加电流的大导体中。如果要把一均衡电流载荷加载到具有集肤效应去的截面上去，则需要给该截面耦合电压降自由度。时间历程后处理器（Time-history Postprocessor）POST26提供了这一功能——得到物理量与时间函数关系的结果。一般的，(General PostProcessor)POST1仅能观察某一时刻整个模型的结果，得到的是物理量与模型尺寸的函数关系。为了用POST26察看物理量与时间的函数关系，需要给物理量事先定义变量，定义完变量后就可以绘制他们的函数曲线了。
下面是一个关于瞬态分析的实例：
如下图说是为某一圆形铁心上的导线槽，槽内导线为铜导线（线圈），极力为无规则时变电压。要求分析在加载时间内导线的磁场分布、响应电流以及交流阻抗情况。
几何常数：h=1.5cm,a=1cm,b=0.5cm,线圈匝数n=500
材料性能：导线和空气的相对磁导率均为μr=1.0，导线电导ρ=5.8*107s/m，铁心相对磁导率μr=∞
激励载荷：加载导线上的电流载荷在0.02s时间内电压从0斜坡变化到10v然后阶跃到13v持续到0.04s

1． 分析
在z圆周截面上可简化为求任一截面的二维平面问题；另外，铁心磁导率相对导线和空气为无穷大，所以可假设在ABC边界上满足诺伊曼边界条件，即磁通量垂直于此边界，并假设在空气隙外的忽略磁泄漏，故在AC上满足狄利克莱边界条件，即磁通量平行于此边界，只对线圈和空气建模。
2． 命令流
/NOPR
/PMETH,OFF,1
KEYW,PR_SET,1
KEYW,PR_STRUC,0
KEYW,PR_THERM,0
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,1
KEYW,MAGNOD,1
KEYW,MAGEDG,0
KEYW,MAGHFE,0
KEYW,MAGELC,0
KEYW,PR_MULTI,0
KEYW,PR_CFD,0
/GO
/COM,
/COM,Preferences for GUI filtering have been set to display:
/COM, Magnetic-Nodal
/PREP7
/TITLE,calculation of impedance
ET,1,PLANE53
ET,2,PLANE53
KEYOPT,2,1,1
KEYOPT,2,2,0
KEYOPT,2,3,0
KEYOPT,2,4,0
KEYOPT,2,5,0
KEYOPT,2,7,0
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,1,,1
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,2,,1
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDE,RSVX,2
MPDE,RSVY,2
MPDE,RSVZ,2
MPDATA,RSVX,2,,1/(5.8e8)
MPLIST,ALL,,,EVLT
/REPLOT,RESIZE
N,1,0,0,0,,,,
N,8,0,0.025,0,,,,
FILL,8,1,6, , ,1,1,1,
FLST,4,8,1,ORDE,2
FITEM,4,1
FITEM,4,-8
NGEN,2,8,P51X, , ,0.005, , ,1,
mat,2
type,2
FLST,2,4,1
FITEM,2,1
FITEM,2,2
FITEM,2,10
FITEM,2,9
E,P51X
FLST,4,1,2,ORDE,1
FITEM,4,1
EGEN,4,1,P51X, , , , , , , , , , ,

meil32

mat,1
type,1
FLST,2,4,1
FITEM,2,5
FITEM,2,13
FITEM,2,14
FITEM,2,6
E,P51X
FLST,4,1,2,ORDE,1
FITEM,4,5
EGEN,3,1,P51X, , , , , , , , , , ,
R,1,0.015*0.005,500, ,1,0.98,
ESEL,S,MAT,,2
NSLE,S
FLST,4,10,1,ORDE,4
FITEM,4,1
FITEM,4,-5
FITEM,4,9
FITEM,4,-13
CP,1,CURR,P51X
NSEL,S,LOC,Y,0.025
FLST,2,2,1,ORDE,2
FITEM,2,8
FITEM,2,16
D,P51X,AZ,0
KEYOPT,2,1,2
KEYOPT,2,2,0
KEYOPT,2,3,0
KEYOPT,2,4,0
KEYOPT,2,5,0
KEYOPT,2,7,0
ALLSEL,ALL
FINISH
/SOLU
ANTYPE,4
FLST,2,4,2,ORDE,2
FITEM,2,1
FITEM,2,-4
BFE,P51X,VLTG,1,10, ,
ALLSEL,ALL
TIME,0.02
AUTOTS,1
NSUBST, , , ,1
KBC,0
TSRES,ERASE
TIME,0.02
AUTOTS,1
DELTIM,0.001, , ,1
KBC,0
TSRES,ERASE
OUTRES,ALL,ALL,
/STATUS,SOLU
SOLVE
FLST,2,4,2,ORDE,2
FITEM,2,1
FITEM,2,-4
BFE,P51X,VLTG,1,14, ,
TIME,0.04
AUTOTS,1
NSUBST, , , ,1
KBC,1
TSRES,ERASE
TIME,0.04
AUTOTS,1
DELTIM,0.004, , ,1
KBC,1
TSRES,ERASE
OUTRES,ALL,ALL,
/STATUS,SOLU
SOLVE
FINISH
/POST26
/UI,COLL,1
NUMVAR,200
SOLU,191,NCMIT
STORE,MERGE
FILLDATA,191,,,,1,1
REALVAR,191,191
FILLDATA,192,,,,0,0
FILLDATA,193,,,,1,0
FILLDATA,194,,,,-1,0
FILLDATA,195,,,,1,1
VARNAME,195,NSET
STORE,MERGE
ESOL,3,3,12,B,X,node12
ESOL,4,1,10,B,X,node10
ESOL,5,6,15,B,X,node15
PLVAR,node10,node12,node15, , , , , , , ,
NSOL,6,12,CURR, ,node12-c
PLVAR, , , , , ,node12-c, , , , ,
/POST1
etable,resistance,nmisc,8
ssum
etable,resistance,nmisc,9
ssum
FINISH

tiangfeng

谢谢楼主
真需要呢！！
呵呵呵！