Exemple #1
0
    def init1(self, dae):
        # check the exclusivity of rate and amount

        # obtain the p0 and q0 at the time of the start
        self.copy_data_ext("StaticLoad", field="p0", dest="p0", idx=self.load)
        self.copy_data_ext("StaticLoad", field="q0", dest="q0", idx=self.load)
        self.copy_data_ext("StaticLoad", field="a", dest="a", idx=self.load)
        self.copy_data_ext("StaticLoad", field="v", dest="v", idx=self.load)

        self.p_out = zeros(self.n, 1)
        self.q_out = zeros(self.n, 1)
Exemple #2
0
    def run_step0(self):
        """
        For the 0th step, store the data and stream data

        Returns
        -------
        None
        """
        dae = self.system.dae
        system = self.system

        self.inc = zeros(dae.m + dae.n, 1)
        system.varout.store(self.t, self.step)

        self.streaming_step()
Exemple #3
0
    def run_step0(self):
        """
        For the 0th step, store the data and stream data

        Returns
        -------
        None
        """
        dae = self.system.dae
        system = self.system
        config = self.config

        # compute line and area flow
        if config.compute_flows:
            dae.init_fg()
            self.compute_flows()  # TODO: move to PowerSystem

        self.inc = zeros(dae.m + dae.n, 1)
        system.varout.store(self.t, self.step)

        self.streaming_step()
Exemple #4
0
    def init1(self, dae):
        self.ace = [[]] * self.n
        for idx, item in enumerate(self.BA):
            self.ace[idx] = self.read_data_ext('BArea', field='ace', idx=item)

        self.en = zeros(self.n, 1)
Exemple #5
0
    def init1(self, dae):
        if globals()['cp'] is None:
            try:
                globals()['cp'] = importlib.import_module('cvxpy')
            except ImportError:
                logger.error(
                    'CVXPY import error. Install optional package `cvxpy` to use AGCMPC'
                )
                sys.exit(1)

        self.t = -1
        self.tlast = -1
        # state array x = [delta, omega, xg1]
        # input array u = [dpin]
        self.copy_data_ext('Governor', field='gen', dest='syn', idx=self.tg)
        self.copy_data_ext('Synchronous',
                           field='delta',
                           dest='delta',
                           idx=self.syn)
        self.copy_data_ext('Synchronous',
                           field='omega',
                           dest='omega',
                           idx=self.syn)
        self.copy_data_ext('Synchronous',
                           field='e1d',
                           dest='e1d',
                           idx=self.syn)
        self.copy_data_ext('Synchronous',
                           field='e1q',
                           dest='e1q',
                           idx=self.syn)
        self.copy_data_ext('Synchronous',
                           field='e2d',
                           dest='e2d',
                           idx=self.syn)
        self.copy_data_ext('Synchronous',
                           field='e2q',
                           dest='e2q',
                           idx=self.syn)

        self.copy_data_ext('Governor', field='xg1', dest='xg1', idx=self.tg)
        self.copy_data_ext('Governor', field='xg2', dest='xg2', idx=self.tg)
        self.copy_data_ext('Governor', field='xg3', dest='xg3', idx=self.tg)
        self.copy_data_ext('Governor', field='pin', dest='pin', idx=self.tg)

        self.copy_data_ext('AVR', field='vm', dest='vm', idx=self.avr)
        self.copy_data_ext('AVR', field='vr1', dest='vr1', idx=self.avr)
        self.copy_data_ext('AVR', field='vr2', dest='vr2', idx=self.avr)
        self.copy_data_ext('AVR', field='vfout', dest='vfout', idx=self.avr)

        dae.y[self.dpin] = 0
        self.dpin0 = zeros(self.n, 1)

        # build state/ input /other algebraic idx array
        self.xidx = matrix([
            self.delta, self.omega, self.e1d, self.e1q, self.e2d, self.e2q,
            self.xg1, self.xg2, self.xg3, self.vm, self.vr1, self.vr2,
            self.vfout
        ])
        self.x0 = dae.x[self.xidx]
        self.x = zeros(len(self.xidx), 1)
        self.dx = zeros(len(self.xidx), 1)
        self.xlast = dae.x[self.xidx]

        self.uidx = matrix([self.dpin])
        self.ulast = zeros(self.n, 1)
        self.dpin_calc = zeros(self.n, 1)

        self.widx = self.system.PQ.a
        self.w0 = self.system.PQ.p0
        self.wlast = matrix(self.w0)

        self.yidx = self.omega
        self.yidx_in_x = [index(self.xidx, y)[0] for y in self.yidx]
        yidx = np.delete(np.arange(dae.m), np.array(self.uidx))
        self.yxidx = matrix(yidx)

        # optimization problem
        self.uvar = cp.Variable((len(self.uidx), self.H + 1), 'u')
        self.uzero = cp.Parameter((len(self.uidx), ), 'u0')
        self.xazero = cp.Parameter((2 * len(self.xidx), 1), 'xa')
        self.prob = None

        self.t_store = []
        self.xpred_store = []