Esempio n. 1
0
 def testTree1(self):
     elements = [ee.createAcY(0,(0,-1),complex(1)),ee.createAcEu(1,(0,1),complex(2)),
         ee.createAcEu(2,(1,-1),complex(-2)),ee.createAcEu(3,(1,2),complex(2)),
         ee.createAcEu(4,(2,-1),complex(2)),ee.createAcY(5,(2,3),complex(1)),
         ee.createAcY(6,(3,-1),complex(1)),ee.createAcEu(7,(3,4),complex(2)),
         ee.createAcY(8,(4,-1),complex(1)),ee.createAcEu(9,(4,5),complex(2)),
         ee.createAcY(10,(5,-1),complex(1)),ee.createAcY(11,(5,6),complex(1)),
         ee.createAcY(12,(6,-1),complex(1)),ee.createAcEu(13,(6,7),complex(2)),
         ee.createAcEu(14,(-1,7),complex(2))]
     def f(e):
         if e.ids[0] > e.ids[1]:
             if e.eType == EType.Eu:
                 e.u = -e.u
             elif e.eType == EType.Ei:
                 e.i = -e.i
             e.ids= e.ids[1],e.ids[0]
         return e
     elements =  map(f,elements)
     elements = sorted(elements, key = lambda e:(max(e.ids),min(e.ids)))
     '''for e in elements:
Esempio n. 2
0
 def test_ac_eu_s_caculate(self):
     pfc = PFC(100,0.01)
     elements = [ee.createAcEu(0,(-1,0),complex(1.055)),ee.createAcS(1,(1,-1),complex(0.12,0.04),complex(1.0,0)),
     ee.createAcS(2,(2,-1),complex(0.1,0.062),complex(1.0,0)),ee.createAcS(3,(3,-1),complex(0.2,0.1),complex(1.0,0)),
     ee.createAcY(4,(-1,0),complex(0,0.00489)),ee.createAcY(5,(0,1),1.0/complex(0.0669,0.105)),
     ee.createAcY(6,(-1,1),complex(0,0.00489)),ee.createAcY(7,(-1,1),complex(0,0.0065)),
     ee.createAcY(8,(1,2),1.0/complex(0.0893,0.14)),ee.createAcY(9,(-1,2),complex(0,0.0065)),
     ee.createAcY(10,(-1,0),complex(0,0.00748)),ee.createAcY(11,(0,2),1.0/complex(0.0781,0.155)),
     ee.createAcY(12,(-1,2),complex(0,0.01561)),ee.createAcY(13,(2,3),1.0/complex(0.0934,0.091)),
     ee.createAcY(14,(-1,3),complex(0,0.01561)),ee.createAcY(15,(-1,2),complex(0,0.00748))]
     pfc.addElement(elements)
     pfc.createYIMatrix()
     pfc.caculate()
     
     print(pfc.U)
     print(pfc.I)
     print("U0:%0.3f,U1:%0.3f,U2:%0.3f,U3:%0.3f" %(abs(pfc.U[0]),abs(pfc.U[1]),abs(pfc.U[2]),abs(pfc.U[3])))
     pfc.getElementResult()
     for s in pfc.electricElements:
        if s.eType != EType.Y:print('eType:{0}, ids:({1},{2}),s:{3}w'.format(s.eType,s.ids[0],s.ids[1],s.s))