用虚位移法计算步进电机齿槽转矩
2006年6月18日 11:02浏览：211995 评论：7
用虚位移法计算步进电机齿槽转矩
程序的命令流如下
/FILNAM,BUJIN
/UNITS,SI
/PREP7
ET,1,13
EMUNIT,MKS
MP,MURX,1,1 !空气
/COM DX = 0.0
/COM DY = 0.0
ANG=0
LOCAL,11,1
K,1,0,0
K,2,0.008,0+ANG
K,3,0.01,0+ANG
K,4,0.008,60+ANG
K,5,0.01,60+ANG
K,6,0.0103,30
K,7,0.0103,-30
K,8,0.0113,30
K,9,0.0113,-30
CSYS,0
K,10,0.0123,0.002 !POLE WIDE
K,11,0.0123,-0.002
K,12,0.0154,0.002
K,13,0.0154,-0.002
CSYS,1
K,14,0.01605,28
K,15,0.01605,-28
K,16,0.01805,45
K,17,0.01805,-45
K,18,0.01605,42.6
K,19,0.01605,-42.6
K,20,0.0103,45 !BIG TEETH
K,21,0.0103,-45
CSYS,0
K,22,0.007951,0.006573
K,23,0.007951,-0.006573
K,24,0.0106,0.0085
K,25,0.0106,-0.0085
CSYS,1
K,26,0.01605,10
K,27,0.01605,-10
K,28,0.01675,45
K,29,0.01675,-45
K,30,0.01535,45
K,31,0.01535,-45
K,32,0.01015,0+ANG
K,33,0.01015,60+ANG
LSTR,1,2 !LINE1_5
LSTR,2,3
LSTR,1,4
LSTR,4,5
LSTR,6,8
LSTR,8,10 !LINE6_10
LSTR,12,10
LSTR,11,13
LSTR,7,9
LSTR,9,11
LSTR,20,22 !LINE11_15
LSTR,22,24
LSTR,21,23
LSTR,25,23
LSTR,30,20
LSTR,16,28 !LINE16_20
LSTR,21,31
LSTR,17,29
LARC,2,4,1,0.008
LARC,3,5,1,0.01
LARC,6,7,1,0.0103 !LINE21_25
LARC,14,26,1,0.01605
LARC,27,15,1,0.01605
LARC,16,17,1,0.01805
L2TAN,12,22
L2TAN,23,14 !LINE26_30
L2TAN,22,7
L2TAN,8,23
LFILLT,5,6,0.0004,0
LFILLT,9,10,0.0004,0
LFILLT,11,12,0.0004,0 !LINE31_35
LFILLT,13,14,0.0004,0
LFILLT,6,7,0.001,0
LFILLT,8,10,0.001,0
LARC,28,18,30,0.0007
LARC,29,19,27,-0.0007 !LINE36_40
LARC,30,18,28,0.0007
LARC,31,19,29,0.0007
LARC,32,33,1,0.01015
LSTR,3,32
LSTR,5,33
LSTR,14,35
LSTR,15,37
ldele,1
ldele,3
LSEL,S,,,2,4,2
LESIZE,ALL,,,4
LSEL,S,,,5,9,4
LESIZE,ALL,,,1
LSEL,S,,,6,10,4
LESIZE,ALL,,,6
LSEL,S,,,7,8,1
LESIZE,ALL,,,4
LSEL,S,,,11,13,2
LESIZE,ALL,,,1
LSEL,S,,,12,14,2
LESIZE,ALL,,,4
LSEL,S,,,15,17,2
LESIZE,ALL,,,6
LSEL,S,,,16,18,2
LESIZE,ALL,,,1
LSEL,S,,,19
LESIZE,ALL,,,10
LSEL,S,,,20
LESIZE,ALL,,,60
LSEL,S,,,21
LESIZE,ALL,,,30
LSEL,S,,,22,23,1
LESIZE,ALL,,,6
LSEL,S,,,24
LESIZE,ALL,,,22
LSEL,S,,,25,26,1
LESIZE,ALL,,,6
LSEL,S,,,27,28,1
LESIZE,ALL,,,1
LSEL,S,,,29,38,1
LESIZE,ALL,,,1
LSEL,S,,,39
LESIZE,ALL,,,60
LSEL,S,,,40,41,1
LESIZE,ALL,,,1
LSEL,S,,,42,43,1
LESIZE,ALL,,,6
ALLSEL
LPLOT
AL,P50X,32
21
9
30
10
34
8
28
23
26
14
32
13
17
38
36
18
24
16
35
37
15
11
31
12
25
22
27
7
33
6
29
5
AGEN,4,1,,,,90
NUMMRG,KP
AL,6,33,7,27,22,42
AGEN,4,5,,,,90
NUMMRG,KP
AL,10,34,8,28,23,43
AGEN,4,9,,,,90
NUMMRG,KP
AL,4,19,2,20
AGEN,6,13,,,,60
NUMMRG,KP
AL,20,40,41,39
AGEN,6,19,,,,60
NUMMRG,KP
CYL4,0,0,0.01805
ASEL,ALL
AOVLAP,ALL !40 IS GAP,35 IS IRON
ADELE,1,4,1,1
ADELE,25,,,0 !!!
ALLSEL
NUMMRG,KP
NUMCMP,ALL !此后面的编号重新定义
KK=1
*IF,KK,GT,0,THEN
LOCAL,11,0
/INPUT,MAT2,SIM
/INPUT,MAT3,SIM
/INPUT,MAT4,SIM
TBCOPY,BH,4,9
TBCOPY,BH,3,4
TBCOPY,BH,3,5
TBCOPY,BH,3,6
TBCOPY,BH,3,7
TBCOPY,BH,3,8
HC = 0.85E6
BR = 1.12
THETA3 = 30+ANG
THETA4 =30-120+ANG
THETA5 = 30+120+ANG
THETA6 = 30+ANG
THETA7 = 30-120+ANG
THETA8 = 30+120+ANG
ASEL,S,,,9,,,1
AATT,3
/COM AATT,9
*AFUN,DEG
MP,MGXX,3,HC*COS(THETA3)
MP,MGYY,3,HC*SIN(THETA3)
ASEL,S,,,10,,,1
AATT,4
/COM AATT,9
*AFUN,DEG
MP,MGXX,4,HC*COS(THETA4)
MP,MGYY,4,HC*SIN(THETA4)
ASEL,S,,,11,,,1
AATT,5
/COM AATT,9
*AFUN,DEG
MP,MGXX,5,HC*COS(THETA5)
MP,MGYY,5,HC*SIN(THETA5)
ASEL,S,,,12,,,1
AATT,6
/COM AATT,9
*AFUN,DEG
MP,MGXX,6,HC*COS(THETA6)
MP,MGYY,6,HC*SIN(THETA6)
ASEL,S,,,13,,,1
AATT,7
/COM AATT,9
*AFUN,DEG
MP,MGXX,7,HC*COS(THETA7)
MP,MGYY,7,HC*SIN(THETA7)
ASEL,S,,,14,,,1
AATT,8
/COM AATT,9
*AFUN,DEG
MP,MGXX,8,HC*COS(THETA8)
MP,MGYY,8,HC*SIN(THETA8)
ASEL,S,,,25
AATT,9
ASEL,S,,,30
AATT,1
ASEL,S,,,15,24,1,1
AATT,1
ASEL,S,,,26,29,1,1
AATT,2
ASEL,S,,,1,8,1,1
AATT,1 !XIAN QUAN
ALLSEL,ALL
SMRTSIZE,6
AMESH,ALL
ALLSEL
/COM SAVE,bujin,DB
*ENDIF
FINISH !前处理结束
/SOLU
ALLSEL
ANTYPE,STATIC
NROPT,AUTO
/COM %%%%% BOUNDARY CONDITION %%%%%
CSYS,1
LSEL,S,LOC,X,0.01805-2*DX,0.01805+2*DX,1
NSLL,S,1
D,ALL,AZ,0
/COM %%%%%FOR FORCE CALCULATION USING VIRTUAL WORKS
ASEL,S,,,9,14,1,1
ASEL,A,,,21
LSLA,S
NSLA,S,1
BF,ALL,MVDI,1
NSEL,INVE
BF,ALL,MVDI,0
/COM %%% LOAD CONDITION %%%
I=0
AS=28.68E-6
N=105
JA=0
JB=I*N/AS
ASEL,S,,,2,8,6,1
ESLA,S
BFE,ALL,JS,,,,JA
ASEL,S,,,4,6,2,1
ESLA,S
BFE,ALL,JS,,,,-JA
ASEL,S,,,1,7,6,1
ESLA,S
BFE,ALL,JS,,,,JB
ASEL,S,,,3,5,2,1
ESLA,S
BFE,ALL,JS,,,,-JB
ALLSEL
NSUBST,1
KBC,1
NEQIT,20
CNVTOL,A,1E-6
LSWRITE
LSSOLVE,1,2
FINISH
/POST1
CSYS,1
*CFOPEN,FRAD,DAT
TOR = 0
FRDY = 0
FRDX = 0
BBSUM = 0
*DO,J,3293,3652,1
ESEL,S,,,J
ETABLE,FVWX,NMISC,3
ETABLE,FVWY,NMISC,4
*GET,FXX,ELEM,J,ETAB,FVWX
*GET,FYY,ELEM,J,ETAB,FVWY
*GET,RAD,ELEM,J,CENT,X
*GET,ANG,ELEM,J,CENT,Y
FT=FXX*SIN(ANG)-FYY*COS(ANG)
FRAD=(FYY*SIN(ANG)+FXX*COS(ANG))*0.026
FRDY = FRDY + FYY
FRDX = FRDX + FXX
TOR = TOR +FT*RAD*0.026
*VWRITE,FRAD
(F10.4)
*ENDDO
*VWRITE,TOR, FRDX*0.026,FRDY*0.026
(3F15.9)
*CFCLOS,FRAD,DAT
ALLSEL
FINISH
程序的命令流如下
/FILNAM,BUJIN
/UNITS,SI
/PREP7
ET,1,13
EMUNIT,MKS
MP,MURX,1,1 !空气
/COM DX = 0.0
/COM DY = 0.0
ANG=0
LOCAL,11,1
K,1,0,0
K,2,0.008,0+ANG
K,3,0.01,0+ANG
K,4,0.008,60+ANG
K,5,0.01,60+ANG
K,6,0.0103,30
K,7,0.0103,-30
K,8,0.0113,30
K,9,0.0113,-30
CSYS,0
K,10,0.0123,0.002 !POLE WIDE
K,11,0.0123,-0.002
K,12,0.0154,0.002
K,13,0.0154,-0.002
CSYS,1
K,14,0.01605,28
K,15,0.01605,-28
K,16,0.01805,45
K,17,0.01805,-45
K,18,0.01605,42.6
K,19,0.01605,-42.6
K,20,0.0103,45 !BIG TEETH
K,21,0.0103,-45
CSYS,0
K,22,0.007951,0.006573
K,23,0.007951,-0.006573
K,24,0.0106,0.0085
K,25,0.0106,-0.0085
CSYS,1
K,26,0.01605,10
K,27,0.01605,-10
K,28,0.01675,45
K,29,0.01675,-45
K,30,0.01535,45
K,31,0.01535,-45
K,32,0.01015,0+ANG
K,33,0.01015,60+ANG
LSTR,1,2 !LINE1_5
LSTR,2,3
LSTR,1,4
LSTR,4,5
LSTR,6,8
LSTR,8,10 !LINE6_10
LSTR,12,10
LSTR,11,13
LSTR,7,9
LSTR,9,11
LSTR,20,22 !LINE11_15
LSTR,22,24
LSTR,21,23
LSTR,25,23
LSTR,30,20
LSTR,16,28 !LINE16_20
LSTR,21,31
LSTR,17,29
LARC,2,4,1,0.008
LARC,3,5,1,0.01
LARC,6,7,1,0.0103 !LINE21_25
LARC,14,26,1,0.01605
LARC,27,15,1,0.01605
LARC,16,17,1,0.01805
L2TAN,12,22
L2TAN,23,14 !LINE26_30
L2TAN,22,7
L2TAN,8,23
LFILLT,5,6,0.0004,0
LFILLT,9,10,0.0004,0
LFILLT,11,12,0.0004,0 !LINE31_35
LFILLT,13,14,0.0004,0
LFILLT,6,7,0.001,0
LFILLT,8,10,0.001,0
LARC,28,18,30,0.0007
LARC,29,19,27,-0.0007 !LINE36_40
LARC,30,18,28,0.0007
LARC,31,19,29,0.0007
LARC,32,33,1,0.01015
LSTR,3,32
LSTR,5,33
LSTR,14,35
LSTR,15,37
ldele,1
ldele,3
LSEL,S,,,2,4,2
LESIZE,ALL,,,4
LSEL,S,,,5,9,4
LESIZE,ALL,,,1
LSEL,S,,,6,10,4
LESIZE,ALL,,,6
LSEL,S,,,7,8,1
LESIZE,ALL,,,4
LSEL,S,,,11,13,2
LESIZE,ALL,,,1
LSEL,S,,,12,14,2
LESIZE,ALL,,,4
LSEL,S,,,15,17,2
LESIZE,ALL,,,6
LSEL,S,,,16,18,2
LESIZE,ALL,,,1
LSEL,S,,,19
LESIZE,ALL,,,10
LSEL,S,,,20
LESIZE,ALL,,,60
LSEL,S,,,21
LESIZE,ALL,,,30
LSEL,S,,,22,23,1
LESIZE,ALL,,,6
LSEL,S,,,24
LESIZE,ALL,,,22
LSEL,S,,,25,26,1
LESIZE,ALL,,,6
LSEL,S,,,27,28,1
LESIZE,ALL,,,1
LSEL,S,,,29,38,1
LESIZE,ALL,,,1
LSEL,S,,,39
LESIZE,ALL,,,60
LSEL,S,,,40,41,1
LESIZE,ALL,,,1
LSEL,S,,,42,43,1
LESIZE,ALL,,,6
ALLSEL
LPLOT
AL,P50X,32
21
9
30
10
34
8
28
23
26
14
32
13
17
38
36
18
24
16
35
37
15
11
31
12
25
22
27
7
33
6
29
5
AGEN,4,1,,,,90
NUMMRG,KP
AL,6,33,7,27,22,42
AGEN,4,5,,,,90
NUMMRG,KP
AL,10,34,8,28,23,43
AGEN,4,9,,,,90
NUMMRG,KP
AL,4,19,2,20
AGEN,6,13,,,,60
NUMMRG,KP
AL,20,40,41,39
AGEN,6,19,,,,60
NUMMRG,KP
CYL4,0,0,0.01805
ASEL,ALL
AOVLAP,ALL !40 IS GAP,35 IS IRON
ADELE,1,4,1,1
ADELE,25,,,0 !!!
ALLSEL
NUMMRG,KP
NUMCMP,ALL !此后面的编号重新定义
KK=1
*IF,KK,GT,0,THEN
LOCAL,11,0
/INPUT,MAT2,SIM
/INPUT,MAT3,SIM
/INPUT,MAT4,SIM
TBCOPY,BH,4,9
TBCOPY,BH,3,4
TBCOPY,BH,3,5
TBCOPY,BH,3,6
TBCOPY,BH,3,7
TBCOPY,BH,3,8
HC = 0.85E6
BR = 1.12
THETA3 = 30+ANG
THETA4 =30-120+ANG
THETA5 = 30+120+ANG
THETA6 = 30+ANG
THETA7 = 30-120+ANG
THETA8 = 30+120+ANG
ASEL,S,,,9,,,1
AATT,3
/COM AATT,9
*AFUN,DEG
MP,MGXX,3,HC*COS(THETA3)
MP,MGYY,3,HC*SIN(THETA3)
ASEL,S,,,10,,,1
AATT,4
/COM AATT,9
*AFUN,DEG
MP,MGXX,4,HC*COS(THETA4)
MP,MGYY,4,HC*SIN(THETA4)
ASEL,S,,,11,,,1
AATT,5
/COM AATT,9
*AFUN,DEG
MP,MGXX,5,HC*COS(THETA5)
MP,MGYY,5,HC*SIN(THETA5)
ASEL,S,,,12,,,1
AATT,6
/COM AATT,9
*AFUN,DEG
MP,MGXX,6,HC*COS(THETA6)
MP,MGYY,6,HC*SIN(THETA6)
ASEL,S,,,13,,,1
AATT,7
/COM AATT,9
*AFUN,DEG
MP,MGXX,7,HC*COS(THETA7)
MP,MGYY,7,HC*SIN(THETA7)
ASEL,S,,,14,,,1
AATT,8
/COM AATT,9
*AFUN,DEG
MP,MGXX,8,HC*COS(THETA8)
MP,MGYY,8,HC*SIN(THETA8)
ASEL,S,,,25
AATT,9
ASEL,S,,,30
AATT,1
ASEL,S,,,15,24,1,1
AATT,1
ASEL,S,,,26,29,1,1
AATT,2
ASEL,S,,,1,8,1,1
AATT,1 !XIAN QUAN
ALLSEL,ALL
SMRTSIZE,6
AMESH,ALL
ALLSEL
/COM SAVE,bujin,DB
*ENDIF
FINISH !前处理结束
/SOLU
ALLSEL
ANTYPE,STATIC
NROPT,AUTO
/COM %%%%% BOUNDARY CONDITION %%%%%
CSYS,1
LSEL,S,LOC,X,0.01805-2*DX,0.01805+2*DX,1
NSLL,S,1
D,ALL,AZ,0
/COM %%%%%FOR FORCE CALCULATION USING VIRTUAL WORKS
ASEL,S,,,9,14,1,1
ASEL,A,,,21
LSLA,S
NSLA,S,1
BF,ALL,MVDI,1
NSEL,INVE
BF,ALL,MVDI,0
/COM %%% LOAD CONDITION %%%
I=0
AS=28.68E-6
N=105
JA=0
JB=I*N/AS
ASEL,S,,,2,8,6,1
ESLA,S
BFE,ALL,JS,,,,JA
ASEL,S,,,4,6,2,1
ESLA,S
BFE,ALL,JS,,,,-JA
ASEL,S,,,1,7,6,1
ESLA,S
BFE,ALL,JS,,,,JB
ASEL,S,,,3,5,2,1
ESLA,S
BFE,ALL,JS,,,,-JB
ALLSEL
NSUBST,1
KBC,1
NEQIT,20
CNVTOL,A,1E-6
LSWRITE
LSSOLVE,1,2
FINISH
/POST1
CSYS,1
*CFOPEN,FRAD,DAT
TOR = 0
FRDY = 0
FRDX = 0
BBSUM = 0
*DO,J,3293,3652,1
ESEL,S,,,J
ETABLE,FVWX,NMISC,3
ETABLE,FVWY,NMISC,4
*GET,FXX,ELEM,J,ETAB,FVWX
*GET,FYY,ELEM,J,ETAB,FVWY
*GET,RAD,ELEM,J,CENT,X
*GET,ANG,ELEM,J,CENT,Y
FT=FXX*SIN(ANG)-FYY*COS(ANG)
FRAD=(FYY*SIN(ANG)+FXX*COS(ANG))*0.026
FRDY = FRDY + FYY
FRDX = FRDX + FXX
TOR = TOR +FT*RAD*0.026
*VWRITE,FRAD
(F10.4)
*ENDDO
*VWRITE,TOR, FRDX*0.026,FRDY*0.026
(3F15.9)
*CFCLOS,FRAD,DAT
ALLSEL
FINISH
技术邻APP
工程师必备
工程师必备
- 项目客服
- 培训客服
- 平台客服
TOP
7
