fini /clear !------Parameter-------- magradau=0.5*90e-3 !Magnet-Radius-Außen magradin=0.5*86e-3 !Magnet-Radius-Magnet-Innen magho=8e-3 !Magnet-Höhe statradau=55e-3 !Stator Außenradius statradin=magradau+0.5e-3 radrotbleau = magradin radrotblein = radrotbleau-15e-3 wgap=3 wmag = 30-wgap PI = ACOS(-1) !------Modellerstellung----------------- /prep7 cof = 900e3 ! Koerzitivfeldstärke http://www.earthmag.de/ndfeb-magnete et,1,plane233, , ,0 ! Magnet, Luft, Platte et,2,infin110,,1,1 ! Rand mp,murx,1,1 ! Luft mp,murx,2,3800 ! Elektroblech, mu_r https://www.google.de/url?sa=t&rct=j&q=&esrc=s&source=web&cd=2&ved=0CD4QFjAB&url=http%3A%2F%2Fwww.thyssenkrupp-steel-europe.com%2Ftiny%2Fbc%2Fdownload.pdf&ei=QmSKUti0POaU4ATUjYCACQ&usg=AFQjCNGbfBciaFzwf78S2R6ZVK-yHMLx3A mp,murx,2,3800 winkel = 90-wgap/2-wmag/2 mp,murx,3,1.07 !Magnet mp,mgxx,3,-cos(winkel/180*PI)*cof mp,mgyy,3,-sin(winkel/180*PI)*cof winkel = winkel-wgap-wmag mp,murx,4,1.07 !Magnet mp,mgxx,4,cos(winkel/180*PI)*cof mp,mgyy,4,sin(winkel/180*PI)*cof winkel = winkel-wgap-wmag mp,murx,5,1.07 !Magnet mp,mgxx,5,-cos(winkel/180*PI)*cof mp,mgyy,5,-sin(winkel/180*PI)*cof winkel = winkel-wgap-wmag mp,murx,6,1.07 !Magnet mp,mgxx,6,cos(winkel/180*PI)*cof mp,mgyy,6,sin(winkel/180*PI)*cof winkel = winkel-wgap-wmag mp,murx,7,1.07 !Magnet mp,mgxx,7,-cos(winkel/180*PI)*cof mp,mgyy,7,-sin(winkel/180*PI)*cof winkel = winkel-wgap-wmag mp,murx,8,1.07 !Magnet mp,mgxx,8,cos(winkel/180*PI)*cof mp,mgyy,8,sin(winkel/180*PI)*cof *IF,0,EQ,1,THEN mp,mgyy,3,0 mp,mgxx,3,cof mp,mgyy,4,0 mp,mgxx,4,cof mp,mgyy,5,0 mp,mgxx,5,cof mp,mgyy,6,0 mp,mgxx,6,cof mp,mgyy,7,0 mp,mgxx,7,cof mp,mgyy,8,0 mp,mgxx,8,cof *ENDIF !----Geometrie------ /pnum,area,1 wgap = 3 wmag = 30-wgap ! ==== Ring 1 ==== ! ---- GAP ------- wstart = 90 wende = wstart-wgap/2 pcirc,radrotblein,radrotbleau,wstart,wende ! Gehäuse pcirc,radrotbleau,statradin,wstart,wende ! Luft pcirc,statradin,statradau,wstart,wende ! Rotorblech ! ---- MAT ------- wstart = wende wende = wstart-wmag pcirc,radrotblein,radrotbleau,wstart,wende ! Gehäuse pcirc,magradau,statradin,wstart,wende ! Luft pcirc,statradin,statradau,wstart,wende ! Rotorblech asel,none pcirc,magradau,magradin,wstart,wende ! Magnet cm,amagnet1,area ! ==== Ring 2 ==== ! ---- GAP ------- wstart = wende wende = wstart-wgap pcirc,radrotblein,radrotbleau,wstart,wende ! Gehäuse pcirc,radrotbleau,statradin,wstart,wende ! Luft pcirc,statradin,statradau,wstart,wende ! Rotorblech ! ---- MAT ------- wstart = wende wende = wstart-wmag pcirc,radrotblein,radrotbleau,wstart,wende ! Gehäuse pcirc,magradau,statradin,wstart,wende ! Luft pcirc,statradin,statradau,wstart,wende ! Rotorblech asel,none pcirc,magradau,magradin,wstart,wende ! Magnet cm,amagnet2,area ! ==== Ring 3 ==== ! ---- GAP ------- wstart = wende wende = wstart-wgap pcirc,radrotblein,radrotbleau,wstart,wende ! Gehäuse pcirc,radrotbleau,statradin,wstart,wende ! Luft pcirc,statradin,statradau,wstart,wende ! Rotorblech ! ---- MAT ------- wstart = wende wende = wstart-wmag pcirc,radrotblein,radrotbleau,wstart,wende ! Gehäuse pcirc,magradau,statradin,wstart,wende ! Luft pcirc,statradin,statradau,wstart,wende ! Rotorblech asel,none pcirc,magradau,magradin,wstart,wende ! Magnet cm,amagnet3,area ! ==== Ring 4 ==== ! ---- GAP ------- wstart = wende wende = wstart-wgap pcirc,radrotblein,radrotbleau,wstart,wende ! Gehäuse pcirc,radrotbleau,statradin,wstart,wende ! Luft pcirc,statradin,statradau,wstart,wende ! Rotorblech ! ---- MAT ------- wstart = wende wende = wstart-wmag pcirc,radrotblein,radrotbleau,wstart,wende ! Gehäuse pcirc,magradau,statradin,wstart,wende ! Luft pcirc,statradin,statradau,wstart,wende ! Rotorblech asel,none pcirc,magradau,magradin,wstart,wende ! Magnet cm,amagnet4,area ! ==== Ring 5 ==== ! ---- GAP ------- wstart = wende wende = wstart-wgap pcirc,radrotblein,radrotbleau,wstart,wende ! Gehäuse pcirc,radrotbleau,statradin,wstart,wende ! Luft pcirc,statradin,statradau,wstart,wende ! Rotorblech ! ---- MAT ------- wstart = wende wende = wstart-wmag pcirc,radrotblein,radrotbleau,wstart,wende ! Gehäuse pcirc,magradau,statradin,wstart,wende ! Luft pcirc,statradin,statradau,wstart,wende ! Rotorblech asel,none pcirc,magradau,magradin,wstart,wende ! Magnet cm,amagnet5,area ! ==== Ring 6 ==== ! ---- GAP ------- wstart = wende wende = wstart-wgap pcirc,radrotblein,radrotbleau,wstart,wende ! Gehäuse pcirc,radrotbleau,statradin,wstart,wende ! Luft pcirc,statradin,statradau,wstart,wende ! Rotorblech ! ---- MAT ------- wstart = wende wende = wstart-wmag pcirc,radrotblein,radrotbleau,wstart,wende ! Gehäuse pcirc,magradau,statradin,wstart,wende ! Luft pcirc,statradin,statradau,wstart,wende ! Rotorblech asel,none pcirc,magradau,magradin,wstart,wende ! Magnet cm,amagnet6,area ! ==== Ring 7 ==== ! ---- GAP ------- wstart = wende wende = wstart-wgap/2 pcirc,radrotblein,radrotbleau,wstart,wende ! Gehäuse pcirc,radrotbleau,statradin,wstart,wende ! Luft pcirc,statradin,statradau,wstart,wende ! Rotorblech !---Flächen aneinanderkleben, wichtig für Vernetzung ----- allsel aglue,all !---CM-wiederherstellen-------- allsel csys,1 asel,s,loc,x,magradin,magradau lsla,s ksll,s ksel,r,loc,x,magradin,magradau LSEL,ALL LSLK,S,1 ASLL,S,1 cm,amagnete,area wstart = 90-wgap/2 wende = wstart-wmag asel,r,loc,y,wstart,wende cm,amagnet1,area cmsel,s,amagnete wstart = wende-wgap/2 wende = wstart-wmag asel,r,loc,y,wstart,wende cm,amagnet2,area cmsel,s,amagnete wstart = wende-wgap/2 wende = wstart-wmag asel,r,loc,y,wstart,wende cm,amagnet3,area cmsel,s,amagnete wstart = wende-wgap/2 wende = wstart-wmag asel,r,loc,y,wstart,wende cm,amagnet4,area cmsel,s,amagnete wstart = wende-wgap/2 wende = wstart-wmag asel,r,loc,y,wstart,wende cm,amagnet5,area cmsel,s,amagnete wstart = wende-wgap/2 wende = wstart-wmag asel,r,loc,y,wstart,wende cm,amagnet6,area csys,0 allsel !----Vernetzung---------------- csys,1 !------Magnet-------- esize,1e-3 mshkey,1 mshape,0,2D type,1 $mat,3 cmsel,s,amagnet1 amesh,amagnet1 type,1 $mat,4 cmsel,s,amagnet2 amesh,amagnet2 type,1 $mat,5 cmsel,s,amagnet3 amesh,amagnet3 type,1 $mat,6 cmsel,s,amagnet4 amesh,amagnet4 type,1 $mat,7 cmsel,s,amagnet5 amesh,amagnet5 type,1 $mat,8 cmsel,s,amagnet6 amesh,amagnet6 !------E-Blech---------- type,1 $mat,2 mshkey,0 ASEL,S,LOC,X,radrotblein,radrotbleau ! Rotor cm,arotor,area CHKMSH,'AREA' CMSEL,S,arotor amesh,arotor type,1 $mat,2 mshkey,1 ASEL,S,LOC,X,statradin,statradau ! Gehäuse cm,astator,area CHKMSH,'AREA' CMSEL,S,astator amesh,astator !------Luft---------- allsel cmsel,s,amagnet1 cmsel,a,amagnet2 cmsel,a,amagnet3 cmsel,a,amagnet4 cmsel,a,amagnet5 cmsel,a,amagnet6 cm,amagnet,area allsel cmsel,u,arotor cmsel,u,astator cmsel,u,amagnet cm,aluft,area type,1 $mat,1 mshkey,0 amesh,aluft csys,0 !------Komponentendefinition------ cmsel,s,amagnete ESEL,ALL ESLA,S cm,emagnete,elem ! Komponentendefinition fmagbc,'emagnete' ! für Kraftauswertung cmsel,s,arotor ESEL,ALL ESLA,S cm,erotor,elem !Komponentendefinition fmagbc,'erotor' !für Kraftauswertung !----Randbedingung----- csys,1 nsel,s,loc,x,statradau esln emodif,all,type,2 sf,all,inf csys,0 !-------Lösung------- /solu allsel solve !-------Auswertung--------------- /post1 plf2d !----Flussdichte-------- /cmap,grau2,cmap /cont,,10 /plopts,minm,0 plns,b,sum !---Kräfte----------- cmsel,s,emagnete cmsel,a,erotor etable,kraft_y1,fmag,y !Maxwell etable,kraft_y2,nmisc,4 !virtual work etable,kraft_x1,fmag,x !Maxwell etable,kraft_x2,nmisc,3 !virtual work ssum *get,kraft1_y,ssum,,item,kraft_y1 *get,kraft2_y,ssum,,item,kraft_y2 *get,kraft1_x,ssum,,item,kraft_x1 *get,kraft2_x,ssum,,item,kraft_x2