Dreiphasentransformator¶
Skript-Datei
------------------------------
-- Allgemeine Einstellungen --
------------------------------
exit_on_error = false -- Verhalten nach Fehler
exit_on_end = true -- Verhalten nach Skriptausführung
verbosity = 2 -- Grad der Bildschirmmeldungen
--------------------------------------------------------
-- Parameterdefinition EI-Kern nach DIN 41 302 Teil 1 --
--------------------------------------------------------
a = 50
b = 40
d1 = 10
delta = 0.5
di1 = 1
di2 = 0.75
h = 3.5
Nc = 150
ur = 1000
----------------------
-- Modellerstellung --
----------------------
new_model_force("example","EI-Kern mit Luftspalt")
global_unit('mm') -- Globale Einheit (m, cm, mm)
pickdist(0.01) -- Abstand Schnappen auf Knotenpunkt
blow_up_wind(-a/2,0,a/2,b+d1) -- Fenstergröße anpassen
cosys('cartes') -- Koordinatensystem
------------------
-- Knotenketten --
------------------
ndtref=d1/10.0 -- Referenzknotenabstand
ndt(ndtref)
nc_line(0,0,a/2,0,0) -- E
nc_line_cont(a/2,b,0)
nc_line(a/2-d1,b,a/2-d1,d1,0)
nc_line_cont(d1/2,d1,0)
nc_line(d1/2,d1,d1/2,b,0)
nc_line(0,0,0,b,0)
nc_line(a/2,b+delta,a/2,b+delta+d1,0) -- I
nc_line_cont(0,b+delta+d1,0)
ndt(ndtref*0.5)
nc_line(a/2,b+delta,0,b+delta,0)
ndt(ndtref)
nc_line_cont(0,b+delta+d1,0)
if (delta>0) then -- Luftspalt
ndt(ndtref*0.5)
nc_line(a/2,b,a/2-d1,b,0)
nc_line(d1/2,b,0,b,0)
nc_line(a/2,b,a/2,b+delta,0)
nc_line(0,b,0,b+delta,0)
nc_line(a/2-d1,b,a/2-d1,b+delta,0)
nc_line(d1/2,b,d1/2,b+delta,0)
ndt(ndtref)
end
-- Bohrungen
if (h>0.0) then
nc_circle_m(a/2-d1/2+h/2,d1/2,a/2-d1/2-h/2,d1/2,a/2-d1/2,d1/2,h/2,pi*h/ndtref)
nc_circle_m(a/2-d1/2-h/2,d1/2,a/2-d1/2+h/2,d1/2,a/2-d1/2,d1/2,h/2,pi*h/ndtref)
nc_circle_m(a/2-d1/2+h/2,b+delta+d1/2,a/2-d1/2-h/2,b+delta+d1/2,a/2-d1/2,b+delta+d1/2,h/2,pi*h/ndtref)
nc_circle_m(a/2-d1/2-h/2,b+delta+d1/2,a/2-d1/2+h/2,b+delta+d1/2,a/2-d1/2,b+delta+d1/2,h/2,pi*h/ndtref)
end
hw=(a/2-1.5*d1)/2-2*di2 -- Wicklungen
nc_rec(d1/2+di2,d1+di1,d1/2+di2+hw,b-di1,"----",0)
nc_line( d1/2+di2+hw,b-di1, d1/2+di2+hw,b+delta, 0)
nc_rec(a/2-d1-di2,d1+di1,a/2-d1-di2-hw,b-di1,"----",0)
nc_line( a/2-d1-di2-hw,b-di1, a/2-d1-di2-hw,b+delta,0)
nc_rec(a/2+di2,d1+di1,a/2+di2+hw,b-di1,"----",0)
ndt(ndtref*3) -- Umhüllung
p1y=-0.15*(b+delta+d1)
p2x=(a/2+di1+hw)*1.25
p2y=1.15*(b+delta+d1)
nc_line(0,p1y,0,0,0)
nc_rec(0,p1y,p2x,p2y,"--- ",0)
nc_line(0,p2y,0,b+delta+d1,0)
--------------------------------
-- Vernetzung und Subregionen --
--------------------------------
create_mesh_se(d1/2,d1/2) -- EI-Schnitt
def_new_subreg(d1/2,d1/2,"EI","darkgreen")
create_mesh_se(a/4.0,b+delta+d1/2)
add_to_subreg(a/4.0,b+delta+d1/2)
create_mesh_se((a/2+di1+hw)*0.625,-0.075*(b+delta+d1)) -- Luft
create_mesh_se(d1/2+2*di2+hw,b)
if (delta>0) then
create_mesh_se(d1/4,b+delta/2)
create_mesh_se(a/2-d1/2,b+delta/2)
end
if (h>0.0) then -- Bohrungen
create_mesh_se(a/2-d1/2,d1/2)
create_mesh_se(a/2-d1/2,b+delta+d1/2)
end
p1x=d1/2+di2+hw/2 -- Wicklungen
p1y=(b-d1)/2+d1
create_mesh_se(p1x,p1y)
p2x=a/2-d1-di2-hw/2
p2y=(b-d1)/2+d1
create_mesh_se(p2x,p2y)
p3x=a/2+di2+hw/2
p3y=(b-d1)/2+d1
create_mesh_se(p3x,p3y)
mirror_nodechains(0,0,0,b+delta+d1) -- Modellvervollständigung
def_bcond_inf() -- Randbedingung
-----------------------
-- Materialzuweisung --
-----------------------
def_mat_fm(d1/2,d1/2,ur,120) -- Blech
-------------------------
-- Wicklungsdefinition --
-------------------------
def_new_wdg(-p3x,p3y,"cyan","Strang 1",Nc,-0.5,wo) -- Wicklungen
add_to_wdg(-p2x,p2y,wsamekey,wi,wser)
def_new_wdg(p1x,p1y,"green","Strang 2",Nc,1.0,wi)
add_to_wdg(-p1x,p1y,wsamekey,wo,wser)
def_new_wdg(p2x,p2y,"red","Strang 3",Nc,-0.5,wo)
add_to_wdg(p3x,p3y,wsamekey,wi,wser)
i1=1
for ii=1,3 do -- Strom in Wicklungen
angle = 2*(ii-2)*math.pi/3.0
def_curr_wdg(ii,i1*math.cos(angle))
end
----------------
-- Berechnung --
----------------
calc_field_single(1,restored,0.01)
color_gradation(d1/2,d1/2,tot,Habs,0,0,"example_Habs.svg")
color_gradation(d1/2,d1/2,tot,Babs,0,0,"example_Babs.svg")
save_model('cont')