def main():
'''
Main function to run the T&D Cosimulation
'''
GlobalData.set_config('config_qsts.json')
GlobalData.set_TDdata()
startTime=time.time()
simulate()
print ('Simulation time: ' + str(time.time()-startTime) + "sec")
generateReport(GlobalData,fname='report.xlsx',sim=GlobalData.config['simulationConfig']['simType'])
return 0
def main():
'''
Main function to run the T&D Cosimulation
'''
GlobalData.set_config('config_td.json')
GlobalData.set_TDdata()
simulate()
generateReport(GlobalData,
fname='report.xlsx',
sim=GlobalData.config['simulationConfig']['simType'])
return 0
def main(configfile):
'''
Main function to run the T&D Cosimulation
'''
GlobalData.set_config(configfile)
GlobalData.set_TDdata()
setOutLocation()
startTime = time.time()
simulate()
# GlobalData.print_data()
print('Simulation time: ' + str(time.time() - startTime) + "sec")
generateReport(GlobalData,
fname=GlobalData.config["outputPath"] + "\\reportfin.xlsx",
sim=GlobalData.config['simulationConfig']['simType'])
return 0
def run(self):
GlobalData.data['monitorData'] = {}
GlobalData.data['TNet']['Dynamic'] = {}
dt = 0.0083333
maxIter = 1
tol = 10**-4
if GlobalData.config['simulationConfig']['protocol'].lower(
) != 'loose_coupling':
maxIter = GlobalData.config['simulationConfig']['protocol'][
'maxiteration']
eventData = GlobalData.config['simulationConfig']['dynamicConfig'][
'events']
events = eventData.keys()
events.sort()
memory_threshold = GlobalData.config['simulationConfig'][
'memoryThreshold']
t = 0
stepCount = 0
stepThreshold = 1e6
lastWriteInd = 0
nPart = 0
for event in events:
while t < eventData[event]['time']:
iteration = 0
mismatch = 100.
if t <= eventData[event]['time'] and t + dt >= eventData[
event]['time']:
if eventData[event]['type'] == 'faultOn':
self._tnet_model.faultOn(
eventData[event]['faultBus'],
eventData[event]['faultImpedance'])
elif eventData[event]['type'] == 'faultOff':
self._tnet_model.faultOff(eventData[event]['faultBus'])
while mismatch > tol and iteration < maxIter:
try:
t = t + dt
iteration += 1
self._tnet_model.runDynamic(t)
print("Sim time: " + str(t))
Vpcc = self._tnet_model.getVoltage()
self._dnet_model.setVoltage(Vpcc)
S = self._dnet_model.getLoad()
self._tnet_model.setLoad(S)
# collect data and store
msg = {}
msg['varName'] = {}
for node in Vpcc:
msg['varName'][node] = ['voltage', 'der']
GlobalData.data['monitorData'][
t] = self._dnet_model.monitor(msg)
GlobalData.data['TNet']['Dynamic'][t] = {
'V': Vpcc,
'S': S
}
# write to disk if running low on memory based on memory threshold
try:
if stepCount == 0:
currentMemUsage = psutil.Process(
).memory_full_info().uss * 1e-6
print("Current Memory: " +
str(currentMemUsage))
elif stepCount == 1:
memUseRate = psutil.Process().memory_full_info(
).uss * 1e-6 - currentMemUsage
stepThreshold = int(memory_threshold /
memUseRate)
print("Step Threshhold: " + str(stepThreshold))
elif stepCount > 1 and stepCount % stepThreshold == 0:
ffullpath = str(
GlobalData.config["outputPath"] +
"\\report{}.xlsx".format(nPart))
stype = str(
GlobalData.config['simulationConfig']
['simType'])
generateReport(GlobalData,
fname=ffullpath,
sim=stype)
print("generated report part" + str(nPart))
GlobalData.data['monitorData'] = {
} # empty data
GlobalData.data['TNet']['Dynamic'] = {
} # empty data
nPart += 1
lastWriteInd = stepCount
stepCount += 1
except:
print(
"Failed write the data on the file to protect the memory"
)
print("Unexpected error:", sys.exc_info()[0])
#mismatch=Vprev-V ##TODO: Need to update mismatch
#TODO: tight_coupling isn't implemented
if GlobalData.config['simulationConfig'][
'protocol'].lower(
) == 'tight_coupling' and mismatch > tol:
print("Tight Couping is not implemented")
except Exception as e:
print('failed to run dynamic simulation')
GlobalData.data['TNet']['Dynamic'][t] = {
'V': Vpcc,
'S': S
}
DSSconvergenceFlg = True
for busID in GlobalData.data['TNet']['LoadBusNumber']:
if busID in GlobalData.data['DNet']['Nodes']:
if S[busID]['convergenceFlg'] is False:
DSSconvergenceFlg = False
break
if DSSconvergenceFlg is False:
print('OpenDSS Convergence Failed')