Esempio n. 1
0
def save_solution_test():
    problem = agros2d.problem()
    field = agros2d.field('magnetic')
    problem.solve()

    values_from_solution = [field.local_values(0.05, 0),
                            field.surface_integrals([0, 1, 2]),
                            field.volume_integrals()]

    from os import path
    filename = '{0}/temp.a2d'.format(path.dirname(pythonlab.tempname()))
    agros2d.save_file(filename, True)
    agros2d.open_file(filename, True)

    field = agros2d.field('magnetic')
    values_from_file = [field.local_values(0.05, 0),
                        field.surface_integrals([0, 1, 2]),
                        field.volume_integrals()]

    return compare(values_from_solution, values_from_file)
Esempio n. 2
0
def save_solution_test():
    problem = agros2d.problem()
    field = agros2d.field('magnetic')
    problem.solve()

    values_from_solution = [
        field.local_values(0.05, 0),
        field.surface_integrals([0, 1, 2]),
        field.volume_integrals()
    ]

    from os import path
    filename = '{0}/temp.a2d'.format(path.dirname(pythonlab.tempname()))
    agros2d.save_file(filename, True)
    agros2d.open_file(filename, True)

    field = agros2d.field('magnetic')
    values_from_file = [
        field.local_values(0.05, 0),
        field.surface_integrals([0, 1, 2]),
        field.volume_integrals()
    ]

    return compare(values_from_solution, values_from_file)
Esempio n. 3
0
 def test(self):
     agros2d.open_file(example)
     agros2d.problem().solve()
Esempio n. 4
0
import agros2d
import pythonlab
from math import pi

agros2d.open_file(
    pythonlab.datadir(
        'resources/examples/Examples/Magnetic Field/Magnetostatics/Switched Reluctance Machine.a2d'
    ))


def get_energy(dphi, Jext):
    problem = agros2d.problem()
    magnetic = agros2d.field("magnetic")

    for coil in Jext:
        magnetic.modify_material(
            "Copper {0}+".format(coil),
            {"magnetic_current_density_external_real": +Jext[coil]})
        magnetic.modify_material(
            "Copper {0}-".format(coil),
            {"magnetic_current_density_external_real": -Jext[coil]})

    geometry = agros2d.geometry
    geometry.select_edges(range(64, 80))
    geometry.rotate_selection(0, 0, dphi)

    problem.solve()
    return magnetic.volume_integrals()['Wm']


def get_density(phi):
import agros2d
import pythonlab
from math import pi

agros2d.open_file('Switched Reluctance Machine.a2d')

def get_energy(dphi, Jext):
    problem = agros2d.problem()
    magnetic = agros2d.field("magnetic")

    for coil in Jext:
        magnetic.modify_material("Copper {0}+".format(coil), {"magnetic_current_density_external_real" : +Jext[coil]})
        magnetic.modify_material("Copper {0}-".format(coil), {"magnetic_current_density_external_real" : -Jext[coil]})

    geometry = agros2d.geometry
    geometry.select_edges(range(64, 80))
    geometry.rotate_selection(0, 0, dphi)

    problem.solve()
    return magnetic.volume_integrals()['Wm']

def get_density(phi):
    Jext = {'A' : 0.0, 'B' : 0.0, 'C' : 0.0}
    J_EXT = 2e6

    if ((phi >= 60 and phi < 90) or
        (phi >= 150 and phi < 180)):
        Jext['A'] = J_EXT

    if ((phi >= 0 and phi < 30) or
        (phi >= 90 and phi < 120) or
Esempio n. 6
0
 def test(self):
     agros2d.open_file(example)
     script = agros2d.get_script_from_model()
     exec script in globals(), locals()
     agros2d.problem().solve()
Esempio n. 7
0
 def test(self):
     agros2d.open_file(example)
     script = agros2d.get_script_from_model()
     exec script in globals(), locals()
     agros2d.problem().solve()
import agros2d
import pythonlab
from math import pi

agros2d.open_file(pythonlab.datadir('resources/examples/Examples/Magnetic Field/Magnetostatics/Switched Reluctance Machine.a2d'))

def get_energy(dphi, Jext):
    problem = agros2d.problem()
    magnetic = agros2d.field("magnetic")

    for coil in Jext:
        magnetic.modify_material("Copper {0}+".format(coil), {"magnetic_current_density_external_real" : +Jext[coil]})
        magnetic.modify_material("Copper {0}-".format(coil), {"magnetic_current_density_external_real" : -Jext[coil]})

    geometry = agros2d.geometry
    geometry.select_edges(range(64, 80))
    geometry.rotate_selection(0, 0, dphi)

    problem.solve()
    return magnetic.volume_integrals()['Wm']

def get_density(phi):
    Jext = {'A' : 0.0, 'B' : 0.0, 'C' : 0.0}
    J_EXT = 2e6

    if ((phi >= 60 and phi < 90) or
        (phi >= 150 and phi < 180)):
        Jext['A'] = J_EXT

    if ((phi >= 0 and phi < 30) or
        (phi >= 90 and phi < 120) or