Induktion/Induktionsverlauf im Element

Das Skript erstellt ein Modell, führt eine Leerlauf und Lastberechnung durch und wertet die Induktion in einem Element und dessen Verlauf über dem Drehwinkel aus. Der Verlauf wird zur weiteren Verwendung in eine Datei geschrieben.

Befehlserläuterung siehe Datenzugriff auf Element und Anwendungsbeispiele für Induktion und Induktionsverlauf in den Elementen.

Skript-Datei

--------------------------------------------------------------------------------
-- Allgemeine Einstellungen ----------------------------------------------------
--------------------------------------------------------------------------------

exit_on_error = false       -- Verhalten nach Fehler
exit_on_end = true          -- Verhalten nach Skriptausfuehrung
verbosity = 2               -- Grad der Bildschirmmeldungen

--------------------------------------------------------------------------------
-- Modellerstellung ------------------------------------------------------------
--------------------------------------------------------------------------------

new_model_force("example","Induc(x)")

global_unit('mm')   -- Globale Einheit (m, cm, mm)
pickdist(0.001)     -- Abstand Schnappen auf Knotenpunkt

cosys('polar')

-- Define FE control data

m.hc_min          =       95.0      -- Limit demagnetisa > 0:[%]Hc,<0:[kA/m]
m.con_hdcopy      =       0.00      -- Hc-copy:Name:auto:0,intact:1, none:-1
m.b_max           =       2.40      -- Max Induction [T] in colorgradation
m.b_min           =       0.00      -- Move inside: 0 , Move outside: > 0
m.calc_fe_loss    =       0.00      -- Calcul. FE-Loss > 0, areas <0, no = 0
m.eval_force      =       0.00      -- Eval. force density > 0, no <= 0
m.allow_draw      =       1.00      -- Draw Graphics :> 0: yes, 0:  no
m.fline_dens      =       1.00      -- F-Lines: 1: small, 2: medium, 3:thick
m.num_flines      =      10.00      -- Number of Field-lines:      < 100 > 2
m.name_bch_log    =       0.00      -- Name bch-file in Logfile:> 0:yes,0:no
m.st_size_move    =       0.00      -- Step size move: r/ph:[degr], x/y:[mm]
m.num_nonl_it     =       1.00      -- Number of nonlinear Iterations   < 99
m.perm_mode       =       0.00      -- Permeability mode:>0:restore,0:actual
m.error_perm      =      0.500E-01  -- Rel. Permeability error < 0.1     [%]
m.allow_demagn    =       0.00      -- Allow Demagnetisation:= 1:yes,= 0:no
m.maenergy        =       0.00      -- Force from magn energy 1 :yes,= 0:no
m.el_order_ag     =       1.00      -- El. order in air gap: lin=1: quadr=2
m.export_scrpt    =       0.00      -- Export parameters in script: yes > 0

pre_models("FE-contr-data");

-- Define basic model parameters

m.tot_num_slot    =         12      -- Numer of Slots
m.num_slots       =       3.00      -- Number of slots in Model (Rot-Mot)
m.num_poles       =          4      -- Number of Poles 2p            (>= 2)
m.npols_gen       =          1      -- Number of Poles simulated     (>= 1)
m.arm_length      =         50      -- Effect. armature length          [mm]
m.fc_radius       =       34.2      -- Radius air-gap center (torque)   [mm]

pre_models("basic_modpar");

-- CAD-Parameter Data:
-- STATOR_3

m.yoke_diam       =       130.      -- Outer diameter yoke   DAS        [mm]
m.inside_diam     =       70.0      -- Inner diameter        DIS        [mm]
m.slot_height     =       20.0      -- Total slot height     HSL        [mm]
m.slot_h1         =       2.00      -- Slot opening height   H1S        [mm]
m.slot_h2         =       4.00      -- Slot head height      H2S        [mm]
m.slot_width      =       3.00      -- Slot width SWI,< 0 groove in Slot[mm]
m.slot_r1         =       0.00      -- Upper Radius          R1S        [mm]
m.slot_r2         =       0.00      -- Lower Radius          R2S        [mm]
m.wedge_width1    =       0.00      -- Distance middlepoints B1S        [mm]
m.wedge_width2    =       0.00      -- Distance middlepoints B2S        [mm]
m.nodedist        =       1.00      -- Rel. node distance in iron > 0.5
m.middle_line     =       0.00      -- Layers: 0,vert:1,horiz:2,vert+horiz:3
m.tooth_width     =       9.00      -- Tooth width TW (replaces B1, B2) [mm]
m.slot_top_sh     =       0.00      -- Top of slot: corner=2, line=1, arc=0
m.airgap          =       1.00      -- Mesh in airgap: 2/3 airgap height[mm]
m.tot_num_sl      =         12      -- Total number of teeth on 2pi
m.num_sl_gen      =          3      -- Number of teeth be generated
m.zeroangl        =       0.00      -- Reference angle to x-axis [grad]
m.mcvkey_yoke     =    'dummy'

pre_models("STATOR_3");

-- CAD-Parameter Data:
-- Magnet-Sector

m.magn_rad        =       33.5      -- Radius Magnet airgap side    RA  [mm]
m.yoke_rad        =       20.0      -- Radius iron yoke inside      RI  [mm]
m.magn_height     =       5.00      -- Magnet height                HM  [mm]
m.magn_width      =       80.0      -- Magnet width: > 0: [%], < 0: BM  [mm]
m.airgap          =       1.00      -- Mesh height (2/3 airgap)         [mm]
m.nodedist        =       1.00      -- Rel. Node distance > 0.2    ..
m.condshaft_r     =       20.0      -- Conduct. Shaft Radius:< yoke_rad [mm]
m.magn_num        =       1.00      -- Number Magnet sections:   > = 1
m.magn_perm       =       1.05      -- Radius Iron between magnets      [mm]
m.magn_l          =       100.      -- Magnet Length: < 150 %            [%]
m.magn_ori        =       1.00      -- Magnet: par = 1, Pol = 2, halbach = 3
m.magn_type       =       2.00      -- Mag Type:Ac:1,Acp:2,Rc:3, Rcp:4,AAp:5
m.magn_shape      =       0.00      -- Magnet surface radius or shape HA[mm]
m.br_height       =       0.00      -- Iron bridge height: = 0: none  BH[mm]
m.br_width        =       0.00      -- Iron bridge width : = 0: none  BW[mm]
m.zeroangl        =       0.00      -- Reference angle to x-axis [grad]
m.cond_shaft      =       0.00      -- El. Conductivity of shaft [S/m]
m.mcvkey_yoke     =    'dummy'
m.mcvkey_mshaft   =    'dummy'

pre_models("Magnet-Sector");

-- connect_models

pre_models("connect_models");

-- Gen_winding

m.num_phases      =       3.00      -- Number of Phases m  <= 500
m.num_layers      =       1.00      -- Number of Layers (slot sides)per slot
m.num_wires       =       20.0      -- Number of wires per slot side
m.current         =       1.00      -- Wdg-Current [A] or flux [Vs/mm],peak
m.coil_span       =       3.00      -- Coil span Y  >= 1
m.num_slots       =       3.00      -- Number of slots in Model (Rot-Mot)
m.mat_type        =       1.00      -- Ma-type:1=Rot;21=lin-x;22=lin-y
m.wind_type       =       1.00      -- W-typ:1=w&cur;2=w&flux;3=bar&cur
m.win_asym        =       1.00      -- asy. slot pitch: 1=sym.; <1=asym.
m.wdg_location    =       1.00      -- Windg location:1: stator, 2: rotor
m.curr_inp        =       0.00      -- Curent:Const:0, Q-axis:1, D-axis:-1
m.dq_offset       =       0.00      -- Offset to D-axis:>=0 Compl;<0 only Re
m.xcoil_1         =     44.433      -- center coordinate of 1. coil side [mm]
m.ycoil_1         =     11.906      -- center coordinate of 1. coil side [mm]
m.xcoil_2         =      0.000      -- center coordinate of 2. coil side [mm]
m.ycoil_2         =      0.000      -- center coordinate of 2. coil side [mm]

pre_models("Gen_winding");

--------------------------------------------------------------------------------
-- Berechnung ------------------------------------------------------------------
--------------------------------------------------------------------------------

-- PM-/Reluktanzmotor-Simulation (Leerlauf und Last)

m.move_action     =       0.000     -- Move Action: rotate=0, linear=1
m.num_pol_pair    = m.num_poles/2   -- Number of Pole pairs           (>= 1)
m.speed           =    1000.000     -- Speed: rotate                 [1/min]
m.skew_angle      =       0.000     -- Skew angle                  [Deg, mm]
m.nu_skew_steps   =       0.000     -- No of skew sect:>0:finite,0:infinite
m.eval_force      =       0.000     -- Evaluate force density: no = 0,yes >0
m.current         =      40.000     -- Nominal stator coil current(Peak) [A]
m.angl_i_up       =      20.000     -- Angle current I vs. voltage Up  [Deg]
m.num_par_wdgs    =       0.000     -- Number of parallel Windings   (>= 1)
m.magn_temp       =      20.000     -- Temperature Magnet            [Deg C]
m.optim_i_up      =       0.000     -- Optimize < I vs Up : no = 0, yes > 0
m.calc_fe_loss    =       1.000     -- Calculate Losses:>0, areas>0, no = 0
m.nu_move_steps   =        91.0     -- Number of move steps
m.range_phi       =       180.0     -- Move range angle
m.phi_start       =       0.000     -- Start angle
m.pm_eff_aktiv    =         100     -- Effektive Magnetlaenge [%]
m.fc_mult_move_type =        0.     -- Type of move path in air gap
m.pocfilename     =   'sin.poc'     -- POC file

run_models("pm_sym_fast");

--------------------------------------------------------------------------------
-- Auswertung ------------------------------------------------------------------
--------------------------------------------------------------------------------

x = 36.96
y = 3.1
elkey = get_elem_key(x,y)
draw_elem(elkey,red,0)

-- Rückgabe der Induktion des letzten Berechnungsschrittes

B1,B2 = get_elem_data("induc",elkey)
print(B1, B2)

-- Rückgabe der Anzahl berechneter Positionen/Zeitpunkte

N = get_sim_data("num_move_steps")

-- Rückgabe der Induktion zum Berechnungsschritt "mstep" des Lastfalls "lstep"

verbosity = 1
data=io.open("elem_data_1.txt","w+")
for i=1,N do
  B1,B2,phi = get_elem_data("induc(x)",elkey,i,3)
--  B1,B2,phi = get_elem_data("induc(x)",x,y,i,3)
  data:write(string.format("%d %g %g %g\n",i,phi,B1,B2))
end
io.close(data)

-- Rückgabe des Induktionsverlaufes über Position oder Zeit für Leerlauf

--B,stat = get_elem_data("induc()",elkey,3)
B,stat = get_elem_data("induc()",x,y,3)

print(stat)
M = table.getn(B)
N = table.getn(B[1])
print(M, N)

data=io.open("elem_data_2.txt","w+")
for i=1,N do
    data:write(string.format("%d %g %g %g\n",i,B[3][i],B[1][i],B[2][i]))
end
io.close(data)

save_model('close')