/// <summary> /// Tries to construct a current for <paramref name="element"/>, returns true on success /// </summary> /// <param name="element"></param> /// <param name="voltageBA">If true, it means that current was calculated for voltage drop from /// <see cref="ITwoTerminal.TerminalA"/> (reference) to <see cref="ITwoTerminal.TerminalB"/>, if false it means that /// that direction was reversed</param> /// <param name="current">Current constructed if successful, null otherwise</param> /// <returns></returns> protected override bool TryConstructCurrent(ITwoTerminal element, bool voltageBA, out IPhasorDomainSignal current) { // Try to get voltage drop across the element if (_VoltageDrops.TryGet(element, out var voltageDrop, voltageBA)) { // If successful, create a new current signal based on it, cache it //current = IoC.Resolve<IPhasorDomainSignal>( // GetPassiveTwoTerminalDCCurrent(voltageDrop, element), // GetPassiveTwoTerminalACCurrentPhasors(voltageDrop, element)); current = IoC.Resolve <IPhasorDomainSignal>( voltageDrop.Phasors.Select((x) => x.Value * element.GetAdmittance(x.Key.Frequency))); // And return success return(true); }
/// <summary> /// Tries to construct a current for <paramref name="element"/>, returns true on success /// </summary> /// <param name="element"></param> /// <param name="voltageBA">If true, it means that current was calculated for voltage drop from /// <see cref="ITwoTerminal.TerminalA"/> (reference) to <see cref="ITwoTerminal.TerminalB"/>, if false it means that /// that direction was reversed</param> /// <param name="current">Current constructed if successful, null otherwise</param> /// <returns></returns> protected override bool TryConstructCurrent(ITwoTerminal element, bool voltageBA, out ITimeDomainSignal current) { // Try to get voltage drop across the element if (_VoltageDrops.TryGet(element, out var voltageDrop, voltageBA)) { // If successful, create a new current signal based on it, cache it var result = IoC.Resolve <ITimeDomainSignalMutable>(voltageDrop.Samples, voltageDrop.TimeStep, voltageDrop.StartTime); // Get the minimum frequency - it is needed for capacitor waveform shifting, check if there are any waveforms, if not // just assign 0 (technically no waveforms result in a zero wave, which is DC) var minACFrequency = voltageDrop.Waveforms.Count > 0 ? voltageDrop.Waveforms.Keys.Min((x) => x.Frequency) : 0; // Current is composed of each voltage waveform times admittance of the element foreach (var waveform in voltageDrop.Waveforms) { // Get magnitude of element's admittance var admittanceMagnitude = element.GetAdmittance(waveform.Key.Frequency).Magnitude; // Current waveform is the product of voltage waveforrm and magnitude var finalWaveform = waveform.Value.Select((x) => x * admittanceMagnitude); // Introduce phase shift for capacitors - but only if minimum AC frequency is greater than 0, if it's not then there were // no AC voltage sources and so no current will flow through any capacitor if (minACFrequency > 0 && element is ICapacitor) { // Each wave has to be shifted by pi / 2 but only in its period. // For example, a wave with frequency 2 times the lowest frequency has to be shifted by total of pi / 4 - because // there are 2 periods of it in the full waveform. This relation is given by minimum frequency / wave frequency finalWaveform = WaveformBuilder.ShiftWaveform(finalWaveform, minACFrequency / waveform.Key.Frequency * Math.PI / 2); } // Add the waveform to the final waveform result.AddWaveform(waveform.Key, finalWaveform); } current = result; // And return success return(true); }