finish /clear,all ! model geometry f =0.001 l_overall =1100*f l_aux_cyl =150*f l_first_slot =192.5*f l_slot_dist =85*f l_slot =15*f l_bore =50*f di_borehole =16*f di_main =400*f di_circ_saw =180*f de_slot =41*f di_aux_cyl =di_main-2*de_slot ! surroundings radius =2*di_main ! misc fe_cond =1/9.93*1000000 curr =800 !******************************************************************************************************* !******************************************************************************************************* !******************************************************************************************************* /PREP7 ! ---- rotor body ----- cylind,,di_aux_cyl/2,0,l_overall,,90 cylind,,di_main/2,l_aux_cyl,l_overall-l_aux_cyl,,90 cylind,,di_aux_cyl/2,0,l_overall,,-90 cylind,,di_main/2,l_aux_cyl,l_overall-l_aux_cyl,,-90 wpoff,di_main/2+di_circ_saw/2-de_slot,0,l_aux_cyl cylind,,di_circ_saw/2,0,l_overall-2*l_aux_cyl,,360 vovlap,all vdele,16,,,1 ! ---- cut the rotor in pieces ---- wpoff,de_slot-di_circ_saw/2,0,l_first_slot blc5,0,0,di_main,di_main wpoff,0,0,l_slot blc5,0,0,di_main,di_main wpoff,0,0,l_slot_dist blc5,0,0,di_main,di_main wpoff,0,0,l_slot blc5,0,0,di_main,di_main wpoff,0,0,l_slot_dist blc5,0,0,di_main,di_main wpoff,0,0,l_slot blc5,0,0,di_main,di_main wpoff,0,0,l_slot_dist blc5,0,0,di_main,di_main wpoff,0,0,l_slot blc5,0,0,di_main,di_main wpoff,0,0,l_slot_dist blc5,0,0,di_main,di_main wpoff,0,0,l_slot blc5,0,0,di_main,di_main vsba,all,3,,delete,delete vsba,all,4,,delete,delete vsba,all,5,,delete,delete vsba,all,6,,delete,delete vsba,all,7,,delete,delete vsba,all,8,,delete,delete vsba,all,9,,delete,delete vsba,all,10,,delete,delete vsba,all,11,,delete,delete vsba,all,12,,delete,delete ! ---- environment ---- wpave,0,0,0 cylind,,di_main,0,l_overall,,90 cylind,,di_main,0,l_overall,,-90 vovlap,all !vglue,all !*********** numcmp,line numcmp,volu numcmp,area !*********** !******************************************************************************************************* !******************************************************************************************************* !******************************************************************************************************* ! material properties mp,murx,1,1.0 mp,rsvx,1,0 mp,rsvy,1,0 mp,rsvz,1,0 tb,bh,2,,50 tbpt,, 0, 0.00 tbpt,, 600, 1.00 tbpt,, 700, 1.10 tbpt,, 900, 1.20 tbpt,, 1200, 1.30 tbpt,, 1600, 1.40 tbpt,, 2500, 1.50 tbpt,, 4100, 1.60 tbpt,, 7000, 1.70 tbpt,, 12000, 1.80 tbpt,, 19700, 1.90 tbpt,, 32000, 2.00 tbpt,, 45000, 2.06 tbpt,, 55000, 2.08 tbpt,,848000, 3.08 mp,rsvx,2,fe_cond mp,rsvy,2,fe_cond mp,rsvz,2,fe_cond allsel vatt,1 vsel,s,,,1,21 vsel,a,,,24,32 vsel,a,,,35,40 vsel,a,,,42,44 vsel,a,,,47,52 vsel,a,,,54,56 vsel,a,,,59,66 vsel,a,,,69,72 vatt,2 allsel /pnum,mat,1 /number,1 vplot !tbplot,bh,2 !******************************************************************************************************* !******************************************************************************************************* !******************************************************************************************************* et,1,solid117,1 ! solid conductors with eddy currents (AZ+VOLT) et,2,solid117,0 ! stranded conductors for air (AZ) ! mesh solid steel !****************** vsel,s,,,1,21 vsel,a,,,24,32 vsel,a,,,35,40 vsel,a,,,42,44 vsel,a,,,47,52 vsel,a,,,54,56 vsel,a,,,59,66 vsel,a,,,69,72 type,1 ! select element type 1 mshape,0,3d ! mesh with hexahedral-shaped elements when Dimension = 3-D mshkey,1 ! use mapped meshing esize,0.1 vmesh,all allsel ! mesh air in slots !******************* vsel,s,,,22,23 vsel,a,,,33,34 vsel,a,,,45,46 vsel,a,,,57,58 vsel,a,,,67,68 type,2 mshape,0,3d mshkey,1 esize,0.1 vmesh,all allsel ! mesh environment !****************** vsel,s,,,41 vsel,a,,,53 type,2 mshape,1,3d ! mesh with tetrahedral-shaped elements when Dimension = 3-D mshkey,0 ! free meshing esize,radius/4 vmesh,all allsel !******************************************************************************************************* !******************************************************************************************************* !******************************************************************************************************* ! boundary conditions !********************* asel,s,,,140 asel,a,,,248 nsla,s,1 da,all,az,0 allsel eplot !******************************************************************************************************* !******************************************************************************************************* !******************************************************************************************************* ! load !******* !**************** !NSLA, Type, NKEY !**************** asel,s,,,1 asel,a,,,32 nsla,s,1 cp,1,volt,all ! define volt as a degree of freedom to all selected nodes *get,n1,node,,num,min ! stores the nodes of the set in the resulting parameter n1 f,n1,amps,curr ! specify a load to all nodes stored in n1 with the ! current curr allsel asel,s,,,2 asel,a,,,36 nsla,s,1 d,all,volt,0 allsel !******************************************************************************************************* !******************************************************************************************************* !*******************************************************************************************************