24小时热门版块排行榜

返回列表

张烨张烨张烨

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 20
帖子: 1
在线:
虫号: 6422753
注册: 2017-04-30
专业: 工程热物理相关交叉领域

[求助] APDL程序，求教错误已有1人参与

求教大神！！！
这个APDL程序有什么错误？
因为只是想得到最终模型，进行优化，
所以几乎不懂这个程序，
标点应该没问题

!******************三维静态磁场矢量法分析******************
/COM,Magnetic-Node
!******************创建物理环境******************
/PREP7
!******************定义分析参数******************
N=1680          !线圈匝数
I=0.5          !每匝通过的电流
!******************定义单元类型******************
ET,1,SOLID97
!******************定义材料属性******************
tb,hb,1,,61    !设置定子的bh曲线50tw290
tb,hb,2,,157    !设置转子的bh曲线10号钢
MP,MURX,3,1.0 !设置空气相对磁导率
MP,MURX,4,1.1 !设置线圈相对磁导率
MP,RSVX,4,3.00e-8 !设置线圈电导率
!******************建立有限元模型******************
BLOCK,-0.055,0.055,-0.040,0,0,0.030,
/AUTO,1
/REP,FAST
BLOCK,-0.032,-0.015,0,0.030,0,0.030,
BLOCK,0.015,0.032,0,0.030,0,0,0.030,
BLOCK,-0.037,-0.010,0.030,0.055,0,0.030,
BLOCK,0.010,0.037,0.030,0.055,0,0.030,
CYC4,0,0.07135,0.04045,,,,0.030
FLST,2,2,6,ORDE,2
FITEM,2,4
FITEM,2,-5
VSBV,P51X,
ALLSEL,ALL
NUMCMP,ALL
ALLSEL,ALL
FLST,2,5,6,ORDE,2
FITEM,2,1
FITEM,2,-5
VGLUE,P51X
NUMCMP,ALL
ALLSEL,ALL
CYC4,0,0.07135,0.040,,,,0.030
!******************创建线圈实体模型
!**左线圈
LOCAL,11,1,-0.015,0,0,0,-90,0,1,1,
WPCSYS,-1,11,
CYL4,0,0,0,0,0.0140,90,0.030
LOCAL,12,0,-0.015,0,0,0,-90,0,1,1,
WPCSYS,-1,12,
BLOCK,0,0.0140,0,-0.030,0,0.030,
BLOCK,0,-0.017,0,0.0140,0,0.030,
LOCAL,13,1,-0.032,0,0,0,-90,0,1,1,
WPCSYS,-1,13,
CYL4,0,0,0,90,0.0140,180,0.030
LOCAL,14,0,-0.032,0,0,0,-90,0,1,1,
WPCSYS,-1,14,
BLOCK,0,-0.0140,0,-0.030,0,0.030,
LOCAL,15,1,-0.032,0,0.030,0,-90,0,1,1,
WPCSYS,-1,15,
CYL4,0,0,0,180,0.0140,270,0.030
LOCAL,16,0,-0.032,0,0.030,0,-90,0,1,1,
WPCSYS,-1,16,
BLOCK,0,-0.017,0,-0.0140,0,0.030,
LOCAL,17,1,-0.015,0,0.030,0,-90,0,1,1,
WPCSYS,-1,17,
CYL4,0,0,0,270,0.0140,360,0.030
!**右线圈
LOCAL,18,1,0.032,0,0,0,-90,0,1,1,
WPCSYS,-1,18,
CYL4,0,0,0,0,0.0140,90,0.030
LOCAL,19,0,0.032,0,0,0,-90,0,1,1,
WPCSYS,-1,19,
BLOCK,0,0.0140,0,-0.030,0,0.030,
BLOCK,0,-0.017,0,0.0140,0,0.030,
LOCAL,20,1,0.015,0,0,0,-90,0,1,1,
WPCSYS,-1,20,
CYL4,0,0,0,90,0.0140,180,0.030
LOCAL,21,0,0.015,0,0,0,-90,0,1,1,
WPCSYS,-1,21,
BLOCK,0,-0.0140,0,-0.030,0,0.030,
LOCAL,22,1,0.015,0,0.030,0,-90,0,1,1,
WPCSYS,-1,22,
CYL4,0,0,0,180,0.0140,270,0.030
LOCAL,23,0,0.015,0,0.030,0,-90,0,1,1,
WPCSYS,-1,23,
BLOCK,0,0.017,0,-0.0140,0,0.030,
LOCAL,24,1,0.032,0,0.030,0,-90,0,1,1,
WPCSYS,-1,24,
CYL4,0,0,0,270,0.0140,360,0.030
WPCSYS,-1,0
CSYS,0
FLST,2,22,6,ORDE,2
FITEM,2,1
FITEM,2,-22
VGLUE,P51X
NUMCMP,ALL
BLOCK,-80/1000,80/1000,-60/1000,180/1000,-20/1000,50/1000,
!创建外侧空气模型
FLST,2,23,6,ORDE,2
FITEM,2,1
FITEM,2,-23
VOVLAP,P51X
NUMCMP,ALL
!******************设置几何体属性******************
ALLSEL,ALL
VSEL,S,,,1,4,1,    !设置定子属性
VSEL,A,,,22
VATT,1,,1,0
ALLSEL,ALL
VSEL,S,,,5       !设置转子属性
VATT,2,,1,0
ALLSEL,ALL
VSEL,S,,,23       !设置空气模型
VATT,3,,1,0
ALLSEL,ALL
VSEL,S,,,6,21,1, !设置线圈模型
VATT,4,,1,0

ALLSEL,ALL
/NUMBER,1       !按照材料编号显示
/PNUM,MAT,1
/REPLOT
!******************划分网格******************
ALLSEL,ALL
ESIZE,2.5/1000,0, !划分定子（齿）
MSHKEY,0
MSHAPE,1,3d
VSEL,S,,,1,2,1,
VMESH,ALL
ALLSEL,ALL
ESIZE,2.5/1000,0,    !划分定子（齿上两个突台）
MSHKEY,0
MSHAPE,1,3d
VSEL,S,,,3,4,1,
VMESH,ALL
ALLSEL,ALL
ESIZE,8/1000,0,    !划分定子（下部）
MSHKEY,0
MSHAPE,1,3d
VSEL,S,,,22
VMESH,ALL
ALLSEL,ALL
ESIZE,6/1000,0, !划分转子
VSEL,S,,,5
VSWEEP,ALL
ALLSEL,ALL
ESIZE,5/1000,0,    !划分线圈
MSHKEY,0
MSHAPE,1,3d
VSEL,S,,,6,21,1,
VMESH,ALL
ALLSEL,ALL
ESIZE,20/1000,0, !划分外侧空气模型
MSHKEY,0
MSHAPE,1,3d
VSEL,S,,,23
VMESH,ALL
/AUTO,1
/REP,FAST
NUMCMP,ALL
!******************把转自定义为一个组件并施加力标志******************
ALLSEL,ALL
VSEL,S,,,5          !选择转子上的所有单元
VPLOT
ESEL,ALL
ESLV,S
EPLOT
CM,zhuanzi,ELEM    !将转子单元生成一个组件“zhuanzi”
ALLSEL,ALL          !选择所有实体
FMAGBC,‘ZHUANZI’ !给转子施加力标志
ALLSEL,ALL          !选择所有实体
VPLOT
!******************施加电流密度******************
!******************左线圈施加电流密度
ALLSEL,ALL
CSYS,11
VSEL,S,,,6
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,11,
BFE,ALL,JS,1,,N*I/0.000420
ALLSEL,ALL
CSYS,12
VSEL,S,,,19
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,12
BFE,ALL,JS,1,,N*I/0.000420,,,
ALLSEL,ALL
CSYS,12
VSEL,S,,,14
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,12
BFE,ALL,JS,1,-N*I/0.000420,,,,
ALLSEL,ALL
CSYS,13
VSEL,S,,,8
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,13
BFE,ALL,JS,1,-N*I/0.000420,,,
ALLSEL,ALL
CSYS,14
VSEL,S,,,18
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,14
BFE,ALL,JS,1,-N*I/0.000420,,,
ALLSEL,ALL
CSYS,15
VSEL,S,,,10
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,15
BFE,ALL,JS,1,,N*I/0.000420,,,
ALLSEL,ALL
CSYS,16
VSEL,S,,,15
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,16,
BFE,ALL,JS,1,,N*I/0.000420,,,,
ALLSEL,ALL
CSYS,17
VSEL,S,,,9
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,17,
BFE,ALL,JS,1,,N*I/0.000420,,,
!******************右线圈施加电流密度
ALLSEL,ALL
CSYS,18
VSEL,S,,,7
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,18,
BFE,ALL,JS,1,,-N*I/0.000420,,,
ALLSEL,ALL
CSYS,19
VSEL,S,,,21
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,19,
BFE,ALL,JS,1,,-N*I/0.000420,,,
ALLSEL,ALL
CSYS,19
VSEL,S,,,16
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,19,
BFE,ALL,JS,1,N*I/0.000420,,,,
ALLSEL,ALL
CSYS,20
VSEL,S,,,11
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,20,
BFE,ALL,JS,1,,-N*I/0.000420,,,
ALLSEL,ALL
CSYS,21
VSEL,S,,,20
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,21,
BFE,ALL,JS,1,,-N*I/0.000420,,,
ALLSEL,ALL
CSYS,22
VSEL,S,,,13
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,22,
BFE,ALL,JS,1,,-N*I/0.000420,,,
ALLSEL,ALL
CSYS,23
VSEL,S,,,17
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,23,
BFE,ALL,JS,1,,-N*I/0.000420,,,,
ALLSEL,ALL
CSYS,24
VSEL,S,,,12
ESEL,ALL
ESLV,S
EPLOT
EMODIF,ALL,ESYS,24,
BFE,ALL,JS,1,,-N*I/0.000420,,,
CSYS,0
ALLSEL,ALL
VPLOT
!******************施加磁力线平行边界条件******************
FINISH
/SOL
FLST,2,3,5,ORDE,3
FITEM,2,104
FITEM,2,106
FITEM,2,-107
DA,P51X,ASYM
FLST,2,3,5,ORDE,3
FITEM,2,103
FITEM,2,105
FITEM,2,108
DA,P51X,ASYM
ALLSEL,ALL
VPLOT
!******************求解******************
ALLSEL,ALL
MAGSOLV,0,3,0.001,,25,
!******************查看转子组件受力******************
/POST1
FMAGSUM,‘ZHUANZI’

IMG20170430202703.jpg