def Static_Analysis( InputMethod=opem.Functions.Get_Input, TestMode=False, PrintMode=True, ReportMode=True): """ Run Larminie-Dicks static analysis. :param InputMethod : input function or input test vector :param TestMode : test mode flag :type InputMethod : dict or Get_Input function object :type TestMode:bool :param PrintMode : print mode control flag (True : print outputs) :type PrintMode:bool :param ReportMode : report mode control flag (True : generate report) :type ReportMode: bool :return: Result as dict """ OutputFile = None CSVFile = None Warning1 = False Warning2 = False I_Warning = 0 Overall_Params_Max = {} Overall_Params_Linear = {} Simulation_Title = "Larminie-Dicks" try: if PrintMode: print("###########") print(Simulation_Title + "-Model Simulation") print("###########") OutputParamsKeys = sorted(OutputParams.keys()) Output_Dict = dict( zip(OutputParamsKeys, [None] * len(OutputParamsKeys))) if not TestMode: Input_Dict = InputMethod(InputParams) else: Input_Dict = InputMethod if PrintMode: print("Analyzing . . .") Input_Dict = opem.Functions.filter_alpha(Input_Dict) Name = Input_Dict["Name"] if ReportMode: OutputFile = opem.Functions.Output_Init( Input_Dict, Simulation_Title, Name) CSVFile = opem.Functions.CSV_Init( OutputParamsKeys, OutputParams, Simulation_Title, Name) HTMLFile = opem.Functions.HTML_Init(Simulation_Title, Name) IEnd = Input_Dict["i-stop"] IStep = Input_Dict["i-step"] Precision = opem.Functions.get_precision(IStep) [i, IEnd, IStep] = opem.Functions.filter_range( Input_Dict["i-start"], IEnd, IStep) I_List = [] Efficiency_List = [] Power_List = [] Vstack_List = [] Power_Thermal_List = [] B = B_Calc(Input_Dict["T"]) while i < IEnd: try: I_List.append(i) Output_Dict["Vcell"] = Vcell_Calc( E0=Input_Dict["E0"], i=i, i_0=Input_Dict["i_0"], i_n=Input_Dict["i_n"], i_L=Input_Dict["i_L"], R_M=Input_Dict["RM"], B=B, A=Input_Dict["A"]) [Warning1, I_Warning] = opem.Functions.warning_check_1( Output_Dict["Vcell"], I_Warning, i, Warning1) Output_Dict["PEM Efficiency"] = Efficiency_Calc( Output_Dict["Vcell"]) Output_Dict["Power"] = Power_Calc(Output_Dict["Vcell"], i) Output_Dict["VStack"] = VStack_Calc( Input_Dict["N"], Output_Dict["Vcell"]) Output_Dict["Power-Stack"] = PowerStack_Calc( Output_Dict["Power"], Input_Dict["N"]) Output_Dict["Power-Thermal"] = Power_Thermal_Calc( VStack=Output_Dict["VStack"], N=Input_Dict["N"], i=i) Vstack_List.append(Output_Dict["VStack"]) Efficiency_List.append(Output_Dict["PEM Efficiency"]) Power_List.append(Output_Dict["Power-Stack"]) Power_Thermal_List.append(Output_Dict["Power-Thermal"]) if ReportMode: opem.Functions.Output_Save( OutputParamsKeys, Output_Dict, OutputParams, i, OutputFile, PrintMode) opem.Functions.CSV_Save( OutputParamsKeys, Output_Dict, i, CSVFile) i = opem.Functions.rounder(i + IStep, Precision) except Exception as e: print(e) i = opem.Functions.rounder(i + IStep, Precision) if ReportMode: opem.Functions.Output_Save( OutputParamsKeys, Output_Dict, OutputParams, i, OutputFile, PrintMode) opem.Functions.CSV_Save( OutputParamsKeys, Output_Dict, i, CSVFile) [Estimated_V, B0, B1] = opem.Functions.linear_plot( x=I_List, y=Vstack_List) Linear_Approx_Params = Linear_Aprox_Params_Calc(B0, B1) Max_Params = Max_Params_Calc(Power_List, Efficiency_List, Vstack_List) Power_Total = Power_Total_Calc(Vstack_List, IStep, Input_Dict["N"]) Overall_Params_Linear["Pmax(L-Approx)"] = Linear_Approx_Params[0] Overall_Params_Linear["V0"] = B0 Overall_Params_Linear["K"] = B1 Overall_Params_Linear["VFC|Pmax(L-Approx)"] = Linear_Approx_Params[1] Overall_Params_Max["Pmax"] = Max_Params["Max_Power"] Overall_Params_Max["VFC|Pmax"] = Max_Params["Max_VStack"] Overall_Params_Max["Efficiency|Pmax"] = Max_Params["Max_EFF"] Overall_Params_Max["Ptotal(Elec)"] = Power_Total[0] Overall_Params_Max["Ptotal(Thermal)"] = Power_Total[1] if ReportMode: OutputFile.close() CSVFile.close() if PrintMode: print(Report_Message) opem.Functions.HTML_Desc( Simulation_Title, Larminiee_Description, HTMLFile) opem.Functions.HTML_Input_Table( Input_Dict=Input_Dict, Input_Params=InputParams, file=HTMLFile) opem.Functions.HTML_Overall_Params_Table( Overall_Params_Max, Overall_Params_Max_Description, file=HTMLFile, header=True) opem.Functions.HTML_Chart( x=str(I_List), y=str(Power_List), color='rgba(255,99,132,1)', x_label="I(A)", y_label="P(W)", chart_name="Power-Stack", size="600px", file=HTMLFile) opem.Functions.HTML_Chart( x=str(I_List), y=[ str(Vstack_List), str(Estimated_V)], color=[ 'rgba(99,100,255,1)', 'rgb(238, 210, 141)'], x_label="I(A)", y_label="V(V)", chart_name=[ "Voltage-Stack", "Linear-Apx"], size="600px", file=HTMLFile) opem.Functions.HTML_Overall_Params_Table( Overall_Params_Linear, Overall_Params_Linear_Description, file=HTMLFile, header=False) opem.Functions.HTML_Chart( x=str(I_List), y=str(Efficiency_List), color='rgb(255, 0, 255)', x_label="I(A)", y_label="EFF", chart_name="Efficiency", size="600px", file=HTMLFile) opem.Functions.HTML_Chart(x=str(list(map(opem.Functions.rounder, Power_List))), y=str(Efficiency_List), color='rgb(238, 210, 141)', x_label="P(W)", y_label="EFF", chart_name="Efficiency vs Power", size="600px", file=HTMLFile) opem.Functions.HTML_Chart( x=str(I_List), y=str(Power_Thermal_List), color='rgb(255, 0, 255)', x_label="I(A)", y_label="P(W)", chart_name="Power(Thermal)", size="600px", file=HTMLFile) opem.Functions.warning_print( warning_flag_1=Warning1, warning_flag_2=Warning2, I_Warning=I_Warning, file=HTMLFile, PrintMode=PrintMode) opem.Functions.HTML_End(HTMLFile) HTMLFile.close() if PrintMode: print("Done!") if not TestMode: if PrintMode: print( "Result In -->" + os.path.join( os.getcwd(), Simulation_Title)) else: return { "Status": True, "P": Power_List, "I": I_List, "V": Vstack_List, "EFF": Efficiency_List, "Ph": Power_Thermal_List, "V0": B0, "K": B1, "VE": Estimated_V} except Exception: if TestMode: return { "Status": False, "Message": "[Error] " + Simulation_Title + " Simulation Failed!(Check Your Inputs)"} print( "[Error] " + Simulation_Title + " Simulation Failed!(Check Your Inputs)")
def Static_Analysis(InputMethod=Get_Input, TestMode=False, PrintMode=True, ReportMode=True): """ This function run Chamberline-Kim static analysis with calling other functions :param InputMethod : Input Function Or Input Test Vector :param TestMode : Test Mode Flag :type InputMethod : dict or Get_Input function object :type TestMode:bool :param PrintMode : Print Mode Control Flag (True : Print Outputs) :type PrintMode:bool :param ReportMode : Report Mode Control Flag (True : Generate Report) :type ReportMode: bool :return: Result as dict """ OutputFile = None CSVFile = None try: Simulation_Title = "Chamberline-Kim" if PrintMode == True: print("###########") print(Simulation_Title + "-Model Simulation") print("###########") OutputParamsKeys = list(OutputParams.keys()) OutputParamsKeys.sort() Output_Dict = dict( zip(OutputParamsKeys, [None] * len(OutputParamsKeys))) if not TestMode: Input_Dict = InputMethod(InputParams) else: Input_Dict = InputMethod if PrintMode == True: print("Analyzing . . .") Name = Input_Dict["Name"] if ReportMode == True: OutputFile = Output_Init(Input_Dict, Simulation_Title, Name) CSVFile = CSV_Init(OutputParamsKeys, OutputParams, Simulation_Title, Name) HTMLFile = HTML_Init(Simulation_Title, Name) IEnd = Input_Dict["i-stop"] IStep = Input_Dict["i-step"] Precision = get_precision(IStep) i = Input_Dict["i-start"] I_List = [] Power_List = [] Vstack_List = [] while i < IEnd: try: I_List.append(i) Output_Dict["Vcell"] = Vcell_Calc( Input_Dict["E0"], Input_Dict["b"], Input_Dict["R"], Input_Dict["m"], Input_Dict["n"], i, Input_Dict["A"]) Output_Dict["PEM Efficiency"] = Efficiency_Calc( Output_Dict["Vcell"]) Output_Dict["Power"] = Power_Calc(Output_Dict["Vcell"], i) Output_Dict["VStack"] = VStack_Calc(Input_Dict["N"], Output_Dict["Vcell"]) Vstack_List.append(Output_Dict["VStack"]) Output_Dict["Power-Stack"] = PowerStack_Calc( Output_Dict["Power"], Input_Dict["N"]) Power_List.append(Output_Dict["Power-Stack"]) if ReportMode == True: Output_Save(OutputParamsKeys, Output_Dict, OutputParams, i, OutputFile, PrintMode) CSV_Save(OutputParamsKeys, Output_Dict, i, CSVFile) i = rounder(i + IStep, Precision) except Exception as e: print(e) i = rounder(i + IStep, Precision) if ReportMode == True: Output_Save(OutputParamsKeys, Output_Dict, OutputParams, i, OutputFile, PrintMode) CSV_Save(OutputParamsKeys, Output_Dict, i, CSVFile) if ReportMode == True: HTML_Chart(x=str(I_List), y=str(Power_List), color='rgba(255,99,132,1)', x_label="I(A)", y_label="P(W)", chart_name="Power-Stack", size="600px", file=HTMLFile) HTML_Chart(x=str(I_List), y=str(Vstack_List), color='rgba(99,100,255,1)', x_label="I(A)", y_label="V(V)", chart_name="Voltage-Stack", size="600px", file=HTMLFile) HTML_Input_Table(Input_Dict=Input_Dict, Input_Params=InputParams, file=HTMLFile) HTML_End(HTMLFile) OutputFile.close() CSVFile.close() HTMLFile.close() if PrintMode == True: print("Done!") if not TestMode: if PrintMode == True: print("Result In -->" + os.path.join(os.getcwd(), Simulation_Title)) else: return {"P": Power_List, "I": I_List, "V": Vstack_List} except Exception: print("[Error] Chamberline-Kim Simulation Failed!(Check Your Inputs)")