def main(): starttime = time.time() #create a office typebuilding prj = Project(load_data=True) prj.type_bldg_office(name="Office1", year_of_construction=1988, number_of_floors=2, height_of_floors=3.5, net_leased_area=100, office_layout=1, window_layout=1, with_ahu=True, construction_type="heavy") #path where the export is stored output_path = os.path.join('D:\Temp', 'OutputData') print(os.path.join(output_path, 'OneBuildingSim')) prj.export_aixlib(building_model="MultizoneEquipped", zone_model="ThermalZoneEquipped", corG=False, path=os.path.join(output_path, 'OneBuildingSim')) """ Now we need to simulate this, therefore we get the names of the current buildings in this project """ buildingNames = [] for bld in prj.buildings: buildingNames.append(bld.name) """ Now we define the output directory where the simulation results should be stored, in addition we need to define the path where the exported models are""" outputDir = "D:/TestCampusSimulation" packageDir = output_path + "/OneBuildingSim" + "/Project" """ Now we need to create a simulation list for buildingspy """ li = [] for bld in prj.buildings: #this is necessary for the correct names in the simulation script name = "Project." + bld.name + "." + bld.name s = si.Simulator(name, "dymola", outputDir, packageDir) li.append(s) po = Pool(processes=3) # i think we can use async here because we do not need a particular order # of the results po.map(simulateCase, li) ### Timer endtime = time.time() print('Simulation lasts: ', endtime - starttime, ' seconds or ', (endtime - starttime) / 60, ' minutes! or', (endtime - starttime) / (60 * 60))
def main(): ''' Main method that configures and runs all simulations ''' # Build list of cases to run li = [] # First model model = 'Buildings.Controls.Continuous.Examples.PIDHysteresis' s = si.Simulator(model, 'dymola', 'case1') s.addParameters({'con.eOn': 0.1}) li.append(s) # second model s = si.Simulator(model, 'dymola', 'case2') s.addParameters({'con.eOn': 1}) li.append(s) # Run all cases in parallel po = Pool() po.map(simulateCase, li)
def main(): starttime = time.time() #Adjust this path to your TEASER teaser Examples path or whatever you want this_path = "D:/GIT/TEASER/teaser/Examples" #path of the buildings xmls input_path = os.path.join(this_path, 'ExampleInputFiles', 'MelatenXML') #path where the export is stored output_path = os.path.join(os.path.dirname(this_path), 'OutputData') info_list = read_XMLs(input_path) prj = create_reference_project(info_list) print(os.path.join(output_path, 'Reference')) prj.export_aixlib(building_model="MultizoneEquipped", zone_model="ThermalZoneEquipped", corG=False, path=os.path.join(output_path, 'Reference')) """ Now we need to simulate this, therefore we get the names of the current buildings in this project """ buildingNames = [] for bld in prj.buildings: buildingNames.append(bld.name) """ Now we define the output directory where the simulation results should be stored, in addition we need to define the path where the exported models are""" outputDir = "D:/TestCampusSimulation" packageDir = output_path + "/Reference" + "/Project" """ Now we need to create a simulation list for buildingspy """ li = [] for bld in prj.buildings: #this is necessary for the correct names in the simulation script name = "Project." + bld.name + "." + bld.name s = si.Simulator(name, "dymola", outputDir, packageDir) li.append(s) po = Pool(processes=3) po.map(simulateCase, li) ### Timer endtime = time.time() print('Simulation lasts: ', endtime - starttime, ' seconds or ', (endtime - starttime) / 60, ' minutes! or', (endtime - starttime) / (60 * 60))
def main(): ''' Main method that configures and runs all simulations ''' import shutil # Build list of cases to run #print 'hi' li = [] model = 'Topology_Analysis_5GDHC.Automated_5GDHC.Test' s = si.Simulator(model, 'dymola', 'case1') s.addParameters({'Office4.connected': False}) li.append(s) # Run all cases in parallel po = Pool() po.map(simulateCase, li)
def main(): li = [] model = 'Topology_Analysis_5GDHC.Automated_5GDHC.Test' i = 0 scenarios = ta.SetSimulations() for activeBuildings in scenarios: for scenario in activeBuildings[3]: i += 1 name = "case" + str(i) s = si.Simulator(model, 'dymola', name) for parameters in scenario: #print (parameters) s.addParameters({parameters[0]: parameters[1]}) li.append(s) # Run all cases in parallel po = Pool() po.map(simulateCase, li)
def main(): ''' Main method that configures and runs all simulations ''' import copy import shutil from multiprocessing import Pool # Build list of cases to run li = [] # First model model = 'Buildings.Controls.Continuous.Examples.PIDHysteresis' s1 = si.Simulator(model, 'dymola') s1.setOutputDirectory('case1') s1.addParameters({'con.eOn': 0.1}) s1.setSolver('dassl') s1.showGUI(False) # Translate the model s1.translate() # Add the model to the list of models to be simulated li.append(s1) # Second model s2 = copy.deepcopy(s1) s2.setOutputDirectory('case2') s2.addParameters({'con.eOn': 1}) li.append(s2) # Run both models in parallel po = Pool() po.map(simulateTranslatedModel, li) # Clean up # Clean up shutil.rmtree('case1') shutil.rmtree('case2') s1.deleteTranslateDirectory()
def simulate_in_dymola(heaPum, data, tableName, tableFileName): """ Evaluate the heat pump performance from the model in Dymola. :param heaPum: Heat pump model (object). :param data: Reference performance data (object). :param tableName: Name of the combiTimeTable. :param tableFileName: Name of the text file containing the combiTimeTable. :return: Performance data of the modeled heat pump (object). .. note:: Performance data from the model is evaluated at the same operating conditions (inlet water temperatures and mass flow rates at the source and load sides) as in the reference data. """ import buildingspy.simulate.Simulator as si from buildingspy.io.outputfile import Reader from scipy.interpolate import interp1d from builtins import str import getpass import os import tempfile # Find absolute path to buildings library packagePath = os.path.normpath( os.path.join(os.path.normpath(os.path.dirname(__file__)), '..', '..', '..', '..', '..', '..')) # Create temporary directory for simulation files dirPrefix = tempfile.gettempprefix() tmpDir = tempfile.mkdtemp(prefix=dirPrefix + '-' + 'HeatPumpCalibration' + '-' + getpass.getuser() + '-') # Set parameters for simulation in Dymola calModelPath = heaPum.modelicaCalibrationModelPath() s = si.Simulator(calModelPath, 'dymola', outputDirectory=tmpDir, packagePath=packagePath) s = heaPum.set_ModelicaParameters(s) m1_flow_nominal = min(data.flowSource) m2_flow_nominal = min(data.flowLoad) tableFilePath = \ str(os.path.join(tmpDir, tableFileName).replace(os.sep, '/')) s.addParameters({ 'm1_flow_nominal': m1_flow_nominal, 'm2_flow_nominal': m2_flow_nominal, 'calDat.fileName': tableFilePath }) # Write CombiTimeTable for dymola data.write_modelica_combiTimeTable(tableName, tmpDir, tableFileName, heaPum.CoolingMode) # Simulation parameters s.setStopTime(len(data.EWT_Source)) s.setSolver('dassl') # Kill the process if it does not finish in 2 minutes s.setTimeOut(120) s.showProgressBar(False) s.printModelAndTime() # s.showGUI(show=True) # s.exitSimulator(exitAfterSimulation=False) s.simulate() # Read results modelName = heaPum.modelicaModelName() ofr = Reader(os.path.join(tmpDir, modelName), 'dymola') (time1, QCon) = ofr.values('heaPum.QCon_flow') (time1, QEva) = ofr.values('heaPum.QEva_flow') (time1, P) = ofr.values('heaPum.P') t = [float(i) + 0.5 for i in range(len(data.EWT_Source))] f_P = interp1d(time1, P) P = f_P(t) f_QCon = interp1d(time1, QCon) QCon = f_QCon(t) f_QEva = interp1d(time1, QEva) QEva = f_QEva(t) # # Clean up # shutil.rmtree('calibrationModel') if heaPum.CoolingMode: Capacity = -QEva HR = QCon else: Capacity = QCon HR = -QEva dymRes = SimulationResults(data.EWT_Source, data.EWT_Load, data.flowSource, data.flowLoad, Capacity, HR, P, 'Modelica') return dymRes
s.showProgressBar(False) s.printModelAndTime() s.simulate() def simulateTranslatedModel(s): s.setStopTime(3600 * 24 * 365) s.setTimeOut(60 * 20) s.showProgressBar(False) s.printModelAndTime() s.simulate_translated() def Run_TA(model, case, urbanOptFile, modelicaDir, alwaysActiveProsumers, loads, resultDir, dynamics): s = si.Simulator(model, 'dymola') s.setSolver('dassl') s.showGUI(False) s.translate() li = [] scenarios = ta.SetSimulations(case=case, urbanOptFile=urbanOptFile, modelicaDir=modelicaDir, alwaysActiveProsumers=alwaysActiveProsumers, loads=loads, dynamics=dynamics) i = 0 for activeBuildings in scenarios: for topology in activeBuildings[3]: i += 1
# Function to set common parameters and to run the simulation def simulateCase(s): s.setStopTime(86400) # Kill the process if it does not finish in 1 minute s.setTimeOut(60) s.showProgressBar(False) s.printModelAndTime() s.simulate() # Main function if __name__ == '__main__': # Build list of cases to run li = [] # First model model = 'Buildings.Controls.Continuous.Examples.PIDHysteresis' s = si.Simulator(model, 'dymola', 'case1') s.addParameters({'con.eOn': 0.1}) li.append(s) # second model s = si.Simulator(model, 'dymola', 'case2') s.addParameters({'con.eOn': 1}) li.append(s) # Run all cases in parallel po = Pool() po.map(simulateCase, li)
FileSimulation = 'PTTL_SF_' # ----------- DEFINE SIMULATION TIME ------------ # 2016.06.29 #2016.06.30 #2016.07.01 #2016.07.04 #2016.07.05 #2016.07.06 StartModTime = [46000, 40000, 37000, 43000, 43000, 41000] StopModTime = [58000, 54000, 53000, 57000, 56000, 50000] Delta_Time = StopModTime[KKK] - StartModTime[KKK] #----------- LOAD AND TRANSLATE MODELICA MODEL ------------ model = 'package_PSA_SFERAII.Simulations.PTTL_SF_basic_' + Days[KKK] PACKAGE = os.path.join( 'C:\Users\susanna\Documents\GitHub\PSA_SFERAII\Modelling/', 'package_PSA_SFERAII') s = si.Simulator(model, 'dymola', outputDirectory=StoreModResult, packagePath=PACKAGE) # Set Time for plot TimeNoInit = 500 # Set a value to phase out the initialization line StartPlotTime = [StartModTime[KKK] + TimeNoInit] StopPlotTime = [StopModTime[KKK]] if OPTIMIZE == 'True': s.setNumberOfIntervals(500) s.setSolver('Dassl') s.translate() def optimizeCoeff(In): eps = In[0] print 'Current eps:', eps
def example_type_district(): """"First thing we need to do is to import our Project API module""" from teaser.project import Project from random import randint import buildingspy.simulate.Simulator as Si import time from multiprocessing import Pool """We instantiate the Project class. The parameter load_data = True indicates that we load the XML data bases into our Project. This can take a few sec.""" starttime = time.time() prj_est1 = Project(load_data=True) prj_est1.name = "EST1" prj_est4 = Project(load_data=True) prj_est4.name = "EST4" prj_est7 = Project(load_data=True) prj_est7.name = "EST7" """The functions starting with type_bldg giving us the opportunity to create the specific type building (e.g. type_bldg_residential). The function automatically calculates all the necessary parameter. If not specified different it uses vdi calculation method.""" number_of_buildings_est1 = 14 for building in range(1,round((number_of_buildings_est1)*0.67)+1): name_help = "Building" + str(building) year_of_construction_help = randint(1960,1980) prj_est1.type_bldg_est1a(name=name_help, year_of_construction=year_of_construction_help, number_of_floors=2, height_of_floors=3.15, net_leased_area=92, with_ahu=False, neighbour_buildings=0, construction_type="heavy") for building in range(round((number_of_buildings_est1)*0.67)+1, number_of_buildings_est1+1): name_help = "Building" + str(building) year_of_construction_help = randint(1960,1980) prj_est1.type_bldg_est1b(name=name_help, year_of_construction=year_of_construction_help, number_of_floors=2, height_of_floors=3.15, net_leased_area=92*2, with_ahu=False, neighbour_buildings=0, construction_type="heavy", number_of_apartments=2) number_of_buildings_est4 = 4 for building in range(1,number_of_buildings_est4+1): name_help = "Building" + str(building) year_of_construction_help = randint(1960,1980) prj_est4.type_bldg_est4b(name=name_help, year_of_construction=year_of_construction_help, number_of_floors=9, height_of_floors=2.6, net_leased_area=417*9, with_ahu=False, neighbour_buildings=2, construction_type="heavy", number_of_apartments=38) number_of_buildings_est7 = 29 for building in range(1,round((number_of_buildings_est7)*0.45)+1): name_help = "Building" + str(building) year_of_construction_help = randint(1900,1918) prj_est7.type_bldg_est7(name=name_help, year_of_construction=year_of_construction_help, number_of_floors=3, height_of_floors=3.88, net_leased_area=65*3, with_ahu=False, neighbour_buildings=2, construction_type="heavy", number_of_apartments=1) for building in range(round((number_of_buildings_est7)*0.45)+1, number_of_buildings_est7+1): name_help = "Building" + str(building) year_of_construction_help = randint(1900,1918) prj_est7.type_bldg_est7(name=name_help, year_of_construction=year_of_construction_help, number_of_floors=3, height_of_floors=3.88, net_leased_area=65*3, with_ahu=False, neighbour_buildings=2, construction_type="heavy", number_of_apartments=2) """To export the parameters to a Modelica record, we use the export_record function. path = None indicates, that we want to store the records in \ TEASER'S Output folder""" prj_est1.export_aixlib(building_model="MultizoneEquipped", zone_model="ThermalZoneEquipped", corG=True, internal_id=None, path=None) prj_est4.export_aixlib(building_model="MultizoneEquipped", zone_model="ThermalZoneEquipped", corG=True, internal_id=None, path=None) prj_est7.export_aixlib(building_model="MultizoneEquipped", zone_model="ThermalZoneEquipped", corG=True, internal_id=None, path=None) """Now we retrofit all buildings in the year 2015 (EnEV2014). \ That includes new insulation layer and new windows. The name is changed \ to Retrofit""" prj_est1.name = "EST1_Retrofit" prj_est1.retrofit_all_buildings(2015) prj_est1.export_aixlib(building_model="MultizoneEquipped", zone_model="ThermalZoneEquipped", corG=True, internal_id=None, path=None) prj_est4.name = "EST4_Retrofit" prj_est4.retrofit_all_buildings(2015) prj_est4.export_aixlib(building_model="MultizoneEquipped", zone_model="ThermalZoneEquipped", corG=True, internal_id=None, path=None) prj_est7.name = "EST7_Retrofit" prj_est7.retrofit_all_buildings(2015) prj_est7.export_aixlib(building_model="MultizoneEquipped", zone_model="ThermalZoneEquipped", corG=True, internal_id=None, path=None) endtime = time.time() print('Pre-processing lasts: ', endtime - starttime, ' seconds or ', (endtime - starttime) / 60, ' minutes! or', (endtime - starttime) / (60 * 60), 'hours.') starttime = time.time() """ Now we define the output directory where the simulation results should be stored, in addition we need to define the path where the exported models are""" outputdir_est1 = "D:/Dymola_workspace/EST1" packagedir_est1 = "C:/Users\mla\TEASEROutput/EST1" outputdir_est1_retrofit = "D:/Dymola_workspace/EST1_Retrofit" packagedir_est1_retrofit = "C:/Users\mla\TEASEROutput/EST1_Retrofit" outputdir_est4 = "D:/Dymola_workspace/EST4" packagedir_est4 = "C:/Users\mla\TEASEROutput/EST4" outputdir_est4_retrofit = "D:/Dymola_workspace/EST4_Retrofit" packagedir_est4_retrofit = "C:/Users\mla\TEASEROutput/EST4_Retrofit" outputdir_est7 = "D:/Dymola_workspace/EST7" packagedir_est7 = "C:/Users\mla\TEASEROutput/EST7" outputdir_est7_retrofit = "D:/Dymola_workspace/EST7_Retrofit" packagedir_est7_retrofit = "C:/Users\mla\TEASEROutput/EST7_Retrofit" """ Now we need to create a simulation list for buildingspy """ li_est1 = [] for bld in prj_est1.buildings: # this is necessary for the correct names in the simulation script name = "EST1." + bld.name + "." + bld.name s = Si.Simulator(name, "dymola", outputdir_est1, packagedir_est1) li_est1.append(s) li_est1_retrofit = [] for bld in prj_est1.buildings: # this is necessary for the correct names in the simulation script name = "EST1_Retrofit." + bld.name + "." + bld.name s = Si.Simulator(name, "dymola", outputdir_est1_retrofit, packagedir_est1_retrofit) li_est1_retrofit.append(s) li_est4 = [] for bld in prj_est4.buildings: # this is necessary for the correct names in the simulation script name = "EST4." + bld.name + "." + bld.name s = Si.Simulator(name, "dymola", outputdir_est4, packagedir_est4) li_est4.append(s) li_est4_retrofit = [] for bld in prj_est4.buildings: # this is necessary for the correct names in the simulation script name = "EST4_Retrofit." + bld.name + "." + bld.name s = Si.Simulator(name, "dymola", outputdir_est4_retrofit, packagedir_est4_retrofit) li_est4_retrofit.append(s) li_est7 = [] for bld in prj_est7.buildings: # this is necessary for the correct names in the simulation script name = "EST7." + bld.name + "." + bld.name s = Si.Simulator(name, "dymola", outputdir_est7, packagedir_est7) li_est7.append(s) li_est7_retrofit = [] for bld in prj_est7.buildings: # this is necessary for the correct names in the simulation script name = "EST7_Retrofit." + bld.name + "." + bld.name s = Si.Simulator(name, "dymola", outputdir_est7_retrofit, packagedir_est7_retrofit) li_est7_retrofit.append(s) po = Pool(processes=3) po.map(simulate_case, li_est1) po.map(simulate_case, li_est1_retrofit) po.map(simulate_case, li_est4) po.map(simulate_case, li_est4_retrofit) po.map(simulate_case, li_est7) po.map(simulate_case, li_est7_retrofit) # Timer endtime = time.time() print('Simulation lasts: ', endtime - starttime, ' seconds or ', (endtime - starttime) / 60, ' minutes! or', (endtime - starttime) / (60 * 60), 'hours.')
SIMULATE = 'False' # SET TIME FOR PLOT TimeNoInit = 500 # Set a value to phase out the initialization line StartPlotTime = [StartModTime[0] + TimeNoInit] StopPlotTime = [StopModTime[0]] Nodes = [2, 5, 10, 20, 50] if SIMULATE == 'True': # LOAD MODELICA MODEL model = 'package_PSA_SFERAII.Simulations.PTTL_SF_basic' s = si.Simulator( model, 'dymola', outputDirectory=ResultDirectory, packagePath= 'C:\Users\susanna\Documents\GitHub\PSA_SFERAII\Modelling/package_PSA_SFERAII' ) s.setNumberOfIntervals(500) s.setSolver('Dassl') for kk in range(len(Nodes)): s.addParameters({'EuroTrough.N': Nodes[kk]}) s.addParameters({'EuroTrough.eps6': 0.943771684737}) s.setResultFile(FileSimulation + str(Nodes[kk])) s.setStartTime(StartModTime[0]) s.setStopTime(StopModTime[0]) s.printModelAndTime() s.simulate() else: s = 'None'