def test_DLLM_valid_grad_TLift(self): OC,wing_param = self.__init_wing_param() x0=wing_param.get_dv_array() DLLM = DLLMTargetLift('Simple',wing_param,OC) F_list_names = DLLM.get_DLLMPost().DEF_F_LIST_NAMES # F_list_names.remove('Lift') def f2(x): wing_param.update_from_x_list(x) DLLM = DLLMTargetLift('Simple',wing_param,OC) DLLM.set_target_Lift(769200.) DLLM.set_F_list_names(F_list_names) DLLM.run_direct() DLLM.run_post() func=DLLM.get_F_list() return func def df2(x): wing_param.update_from_x_list(x) DLLM = DLLMTargetLift('Simple',wing_param,OC) DLLM.set_target_Lift(769200.) DLLM.set_F_list_names(F_list_names) DLLM.run_direct() DLLM.run_post() DLLM.run_adjoint() func_grad=array(DLLM.get_dF_list_dchi()) return func_grad val_grad2=FDValidGrad(2,f2,df2,fd_step=1.e-6) ok2,df_fd2,df2=val_grad2.compare(x0,treshold=1.e-2,split_out=True,return_all=True) assert(ok2)
def test_DLLM_valid_grad_TLift(self): OC, wing_param = self.__init_wing_param() x0 = wing_param.get_dv_array() DLLM = DLLMTargetLift('Simple', wing_param, OC) F_list_names = DLLM.get_DLLMPost().DEF_F_LIST_NAMES F_list_names.remove('Lift') def f2(x): wing_param.update_from_x_list(x) DLLM = DLLMTargetLift('Simple', wing_param, OC) DLLM.set_target_Lift(769200.) DLLM.set_F_list_names(F_list_names) DLLM.run_direct() DLLM.run_post() func = DLLM.get_F_list() return func def df2(x): wing_param.update_from_x_list(x) DLLM = DLLMTargetLift('Simple', wing_param, OC) DLLM.set_target_Lift(769200.) DLLM.set_F_list_names(F_list_names) DLLM.run_direct() DLLM.run_post() DLLM.run_adjoint() func_grad = array(DLLM.get_dF_list_dchi()) return func_grad val_grad2 = FDValidGrad(2, f2, df2, fd_step=1.e-8) ok2, df_fd2, df2 = val_grad2.compare(x0, treshold=1.e-6, split_out=True, return_all=True) assert (ok2)
def f(x): wing_param.update_from_x_list(x) DLLM = DLLMTargetLift('Simple',wing_param,OC) DLLM.set_target_Lift(769200.) DLLM.set_F_list_names(F_list_names) DLLM.run_direct() DLLM.run_post() func=DLLM.get_F_list() return func
def __config_DLLM(self): ERROR_MSG=self.ERROR_MSG+'__config_DLLM: ' input_keys=self.__config_dict.keys() type_key = self.__tag+'.DLLM.type' stype = self.__config_dict[type_key] if stype not in self.POS_SOLVER: raise Exception,ERROR_MSG+'solver_type = '+str(type)+' not in '+str(self.POS_SOLVER)+'. Set to default solver_type = '+str(type) if self.__verbose > 0: print 'Solver type = ',stype if stype == 'Solver': self.__DLLM_solver = DLLMSolver(self.__tag,self.__wing_param,self.__OC) elif stype == 'TargetCl': self.__DLLM_solver = DLLMTargetCl(self.__tag,self.__wing_param,self.__OC) target_Cl_key = self.__tag+'.DLLM.target_Cl' target_Cl = self.__config_dict[target_Cl_key] self.__DLLM_solver.set_target_Cl(target_Cl) elif stype == 'TargetLift': self.__DLLM_solver = DLLMTargetLift(self.__tag,self.__wing_param,self.__OC) target_Lift_key = self.__tag+'.DLLM.target_Lift' target_Lift = self.__config_dict[target_Lift_key] self.__DLLM_solver.set_target_Lift(target_Lift) method_key = self.__tag+'.DLLM.method' if method_key in input_keys: method = self.__config_dict[method_key] self.__DLLM_solver.set_method(method) relax_factor_key = self.__tag+'.DLLM.relax_factor' if relax_factor_key in input_keys: relax_factor = self.__config_dict[relax_factor_key] self.__DLLM_solver.set_relax_factor(relax_factor) stop_residual_key = self.__tag+'.DLLM.stop_residual' if stop_residual_key in input_keys: stop_residual = self.__config_dict[stop_residual_key] self.__DLLM_solver.set_stop_residual(stop_residual) max_iterations_key = self.__tag+'.DLLM.max_iterations' if max_iterations_key in input_keys: max_iterations = self.__config_dict[max_iterations_key] self.__DLLM_solver.set_max_iterations(max_iterations) gamma_file_name_key = self.__tag+'.DLLM.gamma_file_name' if gamma_file_name_key in input_keys: gamma_file_name = self.__config_dict[gamma_file_name_key] self.__DLLM_solver.set_gamma_file_name(gamma_file_name) F_list_names_key = self.__tag+'.DLLM.F_list_names' if F_list_names_key in input_keys: F_list_names = self.__config_dict[F_list_names_key] self.__DLLM_solver.set_F_list_names(F_list_names)
def test_DLLM_valid_TLift(self): OC, wing_param = self.__init_wing_param() DLLM = DLLMTargetLift('Simple', wing_param, OC) DLLM.set_target_Lift(769200.) DLLM.run_direct() DLLM.run_post() F_list = DLLM.get_F_list() F_list_names = DLLM.get_F_list_names() lift_index = F_list_names.index('Lift') Lift = F_list[lift_index] assert (abs(Lift - 769200.) < 1.e-2)
def __config_DLLM(self): WARNING_MSG=self.WARNING_MSG+'__config_DLLM: ' input_keys=self.__config_dict.keys() type_key = self.__tag+'.DLLM.type' type = self.__config_dict[type_key] if type not in self.POS_SOLVER: print WARNING_MSG+'solver_type = '+str(solve_type)+' not in '+str(self.POS_SOLVER)+'. Set to default solver_type = Solver' type='Solver' if type == 'Solver': self.__DLLM_solver = DLLMSolver(self.__tag,self.__wing_param,self.__OC) elif type == 'TargetCl': self.__DLLM_solver = DLLMTargetCl(self.__tag,self.__wing_param,self.__OC) target_Cl_key = self.__tag+'.DLLM.target_Cl' target_Cl = self.__config_dict[target_Cl_key] self.__DLLM_solver.set_target_Cl(target_Cl) elif type == 'TargetLift': self.__DLLM_solver = DLLMTargetLift(self.__tag,self.__wing_param,self.__OC) target_Lift_key = self.__tag+'.DLLM.target_Lift' target_Lift = self.__config_dict[target_Lift_key] self.__DLLM_solver.set_target_Lift(target_Lift) method_key = self.__tag+'.DLLM.method' if method_key in input_keys: method = self.__config_dict[method_key] self.__DLLM_solver.set_method(method) relax_factor_key = self.__tag+'.DLLM.relax_factor' if relax_factor_key in input_keys: relax_factor = self.__config_dict[relax_factor_key] self.__DLLM_solver.set_relax_factor(relax_factor) stop_residual_key = self.__tag+'.DLLM.stop_residual' if stop_residual_key in input_keys: stop_residual = self.__config_dict[stop_residual_key] self.__DLLM_solver.set_stop_residual(stop_residual) max_iterations_key = self.__tag+'.DLLM.max_iterations' if max_iterations_key in input_keys: max_iterations = self.__config_dict[max_iterations_key] self.__DLLM_solver.set_max_iterations(max_iterations) gamma_file_name_key = self.__tag+'.DLLM.gamma_file_name' if gamma_file_name_key in input_keys: gamma_file_name = self.__config_dict[gamma_file_name_key] self.__DLLM_solver.set_gamma_file_name(gamma_file_name) F_list_names_key = self.__tag+'.DLLM.F_list_names' if F_list_names_key in input_keys: F_list_names = self.__config_dict[F_list_names_key] self.__DLLM_solver.set_F_list_names(F_list_names)
def test_DLLM_valid_TLift(self): OC,wing_param = self.__init_wing_param() DLLM = DLLMTargetLift('Simple',wing_param,OC) # F_list=DLLM.get_F_list() F_list_names=DLLM.get_F_list_names() DLLM.set_target_Lift(769200.) DLLM.run_direct() DLLM.run_post() F_list=DLLM.get_F_list() F_list_names=DLLM.get_F_list_names() print F_list_names lift_index = F_list_names.index('Lift') Lift=F_list[lift_index] assert(abs(Lift-769200.)<1.e-2)
def df2(x): wing_param.update_from_x_list(x) DLLM = DLLMTargetLift('Simple',wing_param,OC) DLLM.set_target_Lift(769200.) DLLM.run_direct() DLLM.run_post() DLLM.run_adjoint() func_grad=np.array(DLLM.get_dF_list_dchi()) N = func_grad.shape[0] ndv = func_grad.shape[1] out_grad = np.zeros((N,ndv)) for i in xrange(N): out_grad[i,:] = func_grad[i,:]/Ref_list[i] return out_grad
def __init__(self,Target_Lift, N = 10, verbose=0): """Initialization of DLLM component. DLLM component use target lift capability of DLLM kernel @param Target_Lift : the targeted lift value (float) @param N : integer. Number of discrete section on 1/2 wing @param verbose : integer : verbosity level """ try : float(Target_Lift) except: raise ValueError('You MUST define a float target lift value, get '+str(Target_Lift)+' instead.') self.Target_Lift = Target_Lift self.N = N self.OC = None self.rtwist = np.zeros(N) self.__display_wing_param = True self.__verbose = verbose self.__wing_param = None super(DLLMOpenMDAOComponent, self).__init__() self.OC = OperatingCondition(tag='DLLMOC', atmospheric_model='ISA') self.__set_OC_values() self.__set_wing_param() self.__DLLM = DLLMTargetLift('test', self.__wing_param, self.OC, verbose=self.__verbose) self.__DLLM.set_target_Lift(Target_Lift) self.__DLLM.run_direct() self.__DLLM.run_post()
def test_DLLM_valid_grad_TLift(self): OC,wing_param = self.__init_wing_param() x0=wing_param.get_dv_array() DLLM = DLLMTargetLift('Simple',wing_param,OC) DLLM.set_target_Lift(769200.) DLLM.run_direct() DLLM.run_post() Ref_list=DLLM.get_F_list() def f2(x): wing_param.update_from_x_list(x) DLLM = DLLMTargetLift('Simple',wing_param,OC) DLLM.set_target_Lift(769200.) DLLM.run_direct() DLLM.run_post() func=DLLM.get_F_list() return func/Ref_list def df2(x): wing_param.update_from_x_list(x) DLLM = DLLMTargetLift('Simple',wing_param,OC) DLLM.set_target_Lift(769200.) DLLM.run_direct() DLLM.run_post() DLLM.run_adjoint() func_grad=np.array(DLLM.get_dF_list_dchi()) N = func_grad.shape[0] ndv = func_grad.shape[1] out_grad = np.zeros((N,ndv)) for i in xrange(N): out_grad[i,:] = func_grad[i,:]/Ref_list[i] return out_grad val_grad2=FDValidGrad(2,f2,df2,fd_step=1.e-8) ok2,df_fd2,df2=val_grad2.compare(x0,treshold=1.e-5,split_out=True,return_all=True) assert(ok2)
wing_param.set_value('tip_chord', 1.5) wing_param.set_value('root_height', 1.28) wing_param.set_value('break_height', 0.97) wing_param.set_value('tip_height', 0.33) wing_param.convert_to_design_variable('span', (10., 50.)) wing_param.convert_to_design_variable('sweep', (0., 40.)) wing_param.convert_to_design_variable('break_percent', (20., 40.)) wing_param.convert_to_design_variable('root_chord', (5., 7.)) wing_param.convert_to_design_variable('break_chord', (3., 5.)) wing_param.convert_to_design_variable('tip_chord', (1., 2.)) wing_param.convert_to_design_variable('root_height', (1., 1.5)) wing_param.convert_to_design_variable('break_height', (0.8, 1.2)) wing_param.convert_to_design_variable('tip_height', (0.2, 0.5)) wing_param.build_linear_airfoil(OC, AoA0=-2., Cm0=-0.1, set_as_ref=True) wing_param.build_airfoils_from_ref() wing_param.update() print wing_param DLLM = DLLMTargetLift('TLift', wing_param, OC) DLLM.set_target_Lift(769200.) #DLLM.set_method('scipy') DLLM.run_direct() DLLM.run_post() DLLM.run_adjoint() dF_list_dchi = DLLM.get_dF_list_dchi() # for i,name in enumerate(DLLM.get_F_list_names()): # print 'gradients for '+name+' = ',dF_list_dchi[i]
def df2(x): wing_param.update_from_x_list(x) DLLM = DLLMTargetLift('Simple', wing_param, OC) DLLM.set_target_Lift(769200.) DLLM.set_F_list_names(F_list_names) DLLM.run_direct() DLLM.run_post() DLLM.run_adjoint() func_grad = array(DLLM.get_dF_list_dchi()) return func_grad
def f2(x): wing_param.update_from_x_list(x) DLLM = DLLMTargetLift('Simple', wing_param, OC) DLLM.set_target_Lift(769200.) DLLM.set_F_list_names(F_list_names) DLLM.run_direct() DLLM.run_post() func = DLLM.get_F_list() return func
class DLLMWrapper(): ERROR_MSG = 'ERROR in DLLMWrapper.' WARNING_MSG = 'WARNING in DLLMWrapper.' POS_SOLVER = ['Solver','TargetCl','TargetLift'] POS_FMT = ['list','numpy'] def __init__(self, tag, verbose=1): """ Wrapper for the DLLM solver """ self.__tag = tag self.__verbose = verbose self.__OC = None # Operating condition self.__wing_param = None # Wing_param class self.__DLLM_solver = None # DLLM Solver class self.__config_dict = None self.__out_format = 'list' self.__grad_format = 'list' self.__AoA_id = 'AoA' self.__F_list = None self.__F_list_grad = None #-- Accessors def get_OC(self): return self.__OC def get_wing_param(self): return self.__wing_param def get_DLLM_solver(self): return self.__DLLM_solver def get_tags_x0_and_bounds(self): tags=self.__wing_param.get_dv_id_list() x0=self.__wing_param.get_dv_array() bounds=self.__wing_param.get_bounds_array() return tags,x0,bounds def get_x0_and_bounds(self): x0=self.__wing_param.get_dv_array() bounds=self.__wing_param.get_bounds_array() return x0,bounds def get_x0(self): return self.get_x() def get_x(self): x=self.__wing_param.get_dv_array() return x def get_F_list_names(self): return self.__DLLM_solver.get_F_list_names() def get_F_list(self): return self.__F_list def get_F_list_grad(self): return self.__F_list_grad def get_F_list_and_grad(self): return self.__F_list, self.__F_list_grad def get_F_value(self, F_id): index = self.get_F_list_names().index(F_id) return self.get_F_list()[index] #-- Setters def set_AoA_id(self, AoA_id): self.__AoA_id = AoA_id def set_out_format(self, format): WARNING_MSG=self.WARNING_MSG+'set_out_format: ' if format not in self.POS_FMT: print WARNING_MSG+'format = '+str(format)+' not in '+str(self.POS_FMT)+'. Set to default out format = list' format='list' self.__out_format = format def set_grad_format(self, format): WARNING_MSG=self.WARNING_MSG+'set_grad_format: ' if format not in self.POS_FMT: print WARNING_MSG+'format = '+str(format)+' not in '+str(self.POS_FMT)+'. Set to default grad format = list' format='list' self.__grad_format = format #-- Public methods def config_from_file(self, filename): # Open file f = open(filename, "r") # Load input file into data dictionary config_dict = json.loads(f.read()) # Close input file f.close() self.configure(config_dict) def configure(self, config_dict): self.__config_dict = config_dict self.__config_OC() self.__config_param() self.__config_DLLM() def run(self, x): self.__wing_param.update_from_x_list(x) self.__DLLM_solver.set_geom(self.__wing_param) F_list=self.analysis() return F_list def run_grad(self, x): self.__wing_param.update_from_x_list(x) self.__DLLM_solver.set_geom(self.__wing_param) F_list_grad=self.analysis_grad() return F_list_grad def run_and_grad(self, x): self.__wing_param.update_from_x_list(x) self.__DLLM_solver.set_geom(self.__wing_param) F_list,F_list_grad = self.analysis_and_grad() return F_list,F_list_grad def analysis(self): if self.__verbose > 0: print self.__wing_param self.__DLLM_solver.run_direct() self.__DLLM_solver.run_post() F_list = self.__DLLM_solver.get_F_list() self.__DLLM_solver.export_F_list() if self.__out_format == 'list': F_list=F_list.tolist() self.__F_list = F_list return F_list def analysis_grad(self): if self.__verbose > 0: print self.__wing_param self.__DLLM_solver.run_direct() self.__DLLM_solver.run_post() self.__DLLM_solver.run_adjoint() F_list_grad=self.__DLLM_solver.get_dF_list_dchi() self.__DLLM_solver.export_dF_list_dchi() if self.__grad_format == 'numpy': F_list_grad=numpy.array(F_list_grad) self.__F_list_grad = F_list_grad return F_list_grad def analysis_and_grad(self): if self.__verbose > 0: print self.__wing_param self.__DLLM_solver.run_direct() self.__DLLM_solver.run_post() self.__DLLM_solver.run_adjoint() F_list = self.__DLLM_solver.get_F_list() F_list_grad=self.__DLLM_solver.get_dF_list_dchi() self.__DLLM_solver.export_F_list() self.__DLLM_solver.export_dF_list_dchi() if self.__out_format == 'list': F_list=F_list.tolist() if self.__grad_format == 'numpy': F_list_grad=numpy.array(F_list_grad) self.__F_list = F_list self.__F_list_grad = F_list_grad return F_list,F_list_grad def export_results(self): fid=open(self.__tag+'.res.dat','w') res=[self.__F_list,self.__F_list_grad] cPickle.dump(res, fid) fid.close() def import_results(self): fid=open(self.__tag+'.res.dat','r') res=cPickle.load(fid) fid.close() self.__F_list=res[0] self.__F_list_grad=res[1] #-- Private methods def __config_OC(self): """ Set up the """ self.__OC = OperatingCondition(self.__tag+'.OC') self.__OC.config_from_dict(self.__config_dict) def __config_param(self): """ Set up the __wing_param attribute. -The class associated to __wing_param depends on the chosen parameterization type -The method needs to be updated if a new parameterization is implemented """ input_keys=self.__config_dict.keys() geom_type_key=self.__tag+'.param.geom_type' if geom_type_key in input_keys: geom_type = self.__config_dict[geom_type_key] else: geom_type= 'Broken' n_sect_key=self.__tag+'.param.n_sect' if n_sect_key in input_keys: n_sect = self.__config_dict[n_sect_key] else: n_sect = 20 if geom_type == 'Broken': self.__wing_param = Wing_Broken(self.__tag+'.param', n_sect=n_sect) elif geom_type == 'Elliptic': self.__wing_param = Wing_Elliptic(self.__tag+'.param', n_sect=n_sect) elif geom_type == 'Straight': self.__wing_param = Wing_Straight(self.__tag+'.param', n_sect=n_sect) self.__wing_param.set_AoA_id(self.__AoA_id) self.__wing_param.set_common_OC(self.__OC) self.__wing_param.config_from_dict(self.__config_dict) def __config_DLLM(self): ERROR_MSG=self.ERROR_MSG+'__config_DLLM: ' input_keys=self.__config_dict.keys() type_key = self.__tag+'.DLLM.type' stype = self.__config_dict[type_key] if stype not in self.POS_SOLVER: raise Exception,ERROR_MSG+'solver_type = '+str(type)+' not in '+str(self.POS_SOLVER)+'. Set to default solver_type = '+str(type) if self.__verbose > 0: print 'Solver type = ',stype if stype == 'Solver': self.__DLLM_solver = DLLMSolver(self.__tag,self.__wing_param,self.__OC) elif stype == 'TargetCl': self.__DLLM_solver = DLLMTargetCl(self.__tag,self.__wing_param,self.__OC) target_Cl_key = self.__tag+'.DLLM.target_Cl' target_Cl = self.__config_dict[target_Cl_key] self.__DLLM_solver.set_target_Cl(target_Cl) elif stype == 'TargetLift': self.__DLLM_solver = DLLMTargetLift(self.__tag,self.__wing_param,self.__OC) target_Lift_key = self.__tag+'.DLLM.target_Lift' target_Lift = self.__config_dict[target_Lift_key] self.__DLLM_solver.set_target_Lift(target_Lift) method_key = self.__tag+'.DLLM.method' if method_key in input_keys: method = self.__config_dict[method_key] self.__DLLM_solver.set_method(method) relax_factor_key = self.__tag+'.DLLM.relax_factor' if relax_factor_key in input_keys: relax_factor = self.__config_dict[relax_factor_key] self.__DLLM_solver.set_relax_factor(relax_factor) stop_residual_key = self.__tag+'.DLLM.stop_residual' if stop_residual_key in input_keys: stop_residual = self.__config_dict[stop_residual_key] self.__DLLM_solver.set_stop_residual(stop_residual) max_iterations_key = self.__tag+'.DLLM.max_iterations' if max_iterations_key in input_keys: max_iterations = self.__config_dict[max_iterations_key] self.__DLLM_solver.set_max_iterations(max_iterations) gamma_file_name_key = self.__tag+'.DLLM.gamma_file_name' if gamma_file_name_key in input_keys: gamma_file_name = self.__config_dict[gamma_file_name_key] self.__DLLM_solver.set_gamma_file_name(gamma_file_name) F_list_names_key = self.__tag+'.DLLM.F_list_names' if F_list_names_key in input_keys: F_list_names = self.__config_dict[F_list_names_key] self.__DLLM_solver.set_F_list_names(F_list_names)
from DLLM.DLLMGeom.wing_elliptic import Wing_Elliptic from DLLM.DLLMKernel.DLLMTargetLift import DLLMTargetLift OC=OperatingCondition('cond1', atmospheric_model='ISA') OC.set_Mach(0.8) OC.set_AoA(3.0) OC.set_altitude(10000.) OC.set_T0_deg(15.) OC.set_P0(101325.) OC.set_humidity(0.) OC.compute_atmosphere() wing_param=Wing_Elliptic('elliptic_wing',n_sect=20) wing_param.import_BC_from_file('input_parameters.par') wing_param.build_linear_airfoil(OC, AoA0=0.0, set_as_ref=True) wing_param.build_airfoils_from_ref() wing_param.update() wing_param.plot() DLLM = DLLMTargetLift('TLift',wing_param,OC) DLLM.set_target_Lift(769200.) DLLM.run_direct() DLLM.run_post() DLLM.run_adjoint() dF_list_dchi=DLLM.get_dF_list_dchi() print 'Lift = ',DLLM.get_DLLMPost().Lift,' Target = ',769200. print 'AoA = ',OC.get_AoA() dF_list_dchi=DLLM.get_dF_list_dchi()
def df(x): wing_param.update_from_x_list(x) DLLM = DLLMTargetLift('Simple',wing_param,OC) DLLM.set_target_Lift(769200.) DLLM.set_F_list_names(F_list_names) DLLM.run_direct() DLLM.run_post() DLLM.run_adjoint() func_grad=numpy.array(DLLM.get_dF_list_dchi()) return func_grad
class DLLMOpenMDAOComponent(Component): # set up interface to the framework # pylint: disable-msg=E1101 # Outputs of lifting line problem """OpenMDAO component for DLLM implementation """ Lift = Float(iotype='out', desc='Lift') Drag = Float(iotype='out', desc='Drag') Drag_Pressure = Float(iotype='out', desc='Drag_Pressure') Drag_Induced = Float(iotype='out', desc='Drag_Induced') Drag_Wave = Float(iotype='out', desc='Drag_Wave') Drag_Friction = Float(iotype='out', desc='Drag_Friction') Cd = Float(iotype='out', desc='Cd') Cdp = Float(iotype='out', desc='Cdp') Cdi = Float(iotype='out', desc='Cdi') Cdw = Float(iotype='out', desc='Cdw') Cdf = Float(iotype='out', desc='Cdf') Cl = Float(iotype='out', desc='Cl') LoD = Float(iotype='out', desc='LoD') Sref = Float(iotype='out', desc='Sref') # Design variables of lifting line problem rtwist = Array([], desc='rtwist', iotype="in") span = Float(desc='span', default_value=34., iotype="in") sweep = Float(desc='sweep', default_value=34., iotype="in") break_percent = Float(desc='break_percent', default_value=33., iotype="in") root_chord = Float(desc='root_chord', default_value=6.1, iotype="in") break_chord = Float(desc='break_chord', default_value=4.6, iotype="in") tip_chord = Float(desc='tip_chord', default_value=1.5, iotype="in") root_height = Float(desc='root_height', default_value=1.28, iotype="in") break_height = Float(desc='break_height', default_value=0.97, iotype="in") tip_height = Float(desc='tip_height', default_value=0.33, iotype="in") # Operating conditions variables Mach = Float(iotype='in', default_value=0.7, desc='Mach') altitude = Float(iotype='in', default_value=10000., desc='Altitude') T0 = Float(iotype='in',default_value=OperatingCondition.T0, desc='Ground ISA ref Temperature') P0 = Float(iotype='in',default_value=OperatingCondition.P0, desc='Ground ISA ref Pressure') def __init__(self,Target_Lift, N = 10, verbose=0): """Initialization of DLLM component. DLLM component use target lift capability of DLLM kernel @param Target_Lift : the targeted lift value (float) @param N : integer. Number of discrete section on 1/2 wing @param verbose : integer : verbosity level """ try : float(Target_Lift) except: raise ValueError('You MUST define a float target lift value, get '+str(Target_Lift)+' instead.') self.Target_Lift = Target_Lift self.N = N self.OC = None self.rtwist = np.zeros(N) self.__display_wing_param = True self.__verbose = verbose self.__wing_param = None super(DLLMOpenMDAOComponent, self).__init__() self.OC = OperatingCondition(tag='DLLMOC', atmospheric_model='ISA') self.__set_OC_values() self.__set_wing_param() self.__DLLM = DLLMTargetLift('test', self.__wing_param, self.OC, verbose=self.__verbose) self.__DLLM.set_target_Lift(Target_Lift) self.__DLLM.run_direct() self.__DLLM.run_post() def __set_wing_param(self, wing_param_name='test_param'): """Method for wing parameters setting : design variables initial values and bounds @param wing_param_name : wing parametrization names """ self.__set_wing_param_values(wing_param_name=wing_param_name) self.__set_wing_param_bounds() self.__wing_param.build_linear_airfoil(self.OC, AoA0=-2., Cm0=-0.1, set_as_ref=True) self.__wing_param.build_airfoils_from_ref() self.__wing_param.update() if self.__display_wing_param: self.__display_wing_param = False print self.__wing_param def __set_wing_param_values(self, wing_param_name='test_param'): """Method for wing parameters variables setting @param wing_param_name : wing parametrization names """ self.__wing_param = Wing_param(wing_param_name, geom_type='Broken', n_sect=self.N * 2) self.__wing_param.build_wing() for i in xrange(self.N): self.__wing_param.set_value('rtwist%s' % i, 0.) for param in Wing_param.DISCRETE_ATTRIBUTES_LIST: self.__wing_param.set_value(param, getattr(self, param)) def __set_wing_param_bounds(self): """Method for desing variables bounds settings Values are set to inf/-inf in DLLM component and their 'real' bounds are defined when optimization problem is set """ for i in xrange(self.N): self.__wing_param.convert_to_design_variable( 'rtwist%s' % i, (-float('inf'), float('inf'))) for param in Wing_param.DISCRETE_ATTRIBUTES_LIST: self.__wing_param.convert_to_design_variable(param, (-float('inf'), float('inf'))) return def __set_OC_values(self): """ Set default operating conditions""" self.OC.set_Mach(self.Mach) self.OC.set_altitude(self.altitude) self.OC.set_T0(self.T0) self.OC.set_P0(self.P0) self.OC.set_humidity(0.) self.OC.compute_atmosphere() def __update_wing_parame_values(self): for dv_id in self.__wing_param.get_dv_id_list(): if dv_id.startswith('rtwist'): i_twist = int(dv_id.replace('rtwist', '')) self.__wing_param.set_value(dv_id, self.rtwist[i_twist]) else: self.__wing_param.set_value(dv_id,getattr(self,dv_id)) self.__wing_param.build_linear_airfoil(self.OC, AoA0=-2., Cm0=-0.1, set_as_ref=True) self.__set_wing_param_bounds() self.__wing_param.build_airfoils_from_ref() self.__wing_param.update() return def execute(self): """ Perform a DLLM computation with the """ self.__update_wing_parame_values() self.__set_OC_values() self.__DLLM.set_target_Lift(self.Target_Lift) self.__DLLM.run_direct() self.__DLLM.run_post() output = self.__DLLM.get_F_list() for f, f_name in zip(output, self.__DLLM.get_F_list_names()): setattr(self, f_name, f) def list_deriv_vars(self): """specify the inputs and outputs where derivatives are defined Specific treatment for twist : defined as rtwist0, rtwist1,... in DLLM but as an array rtwist in openmdao component""" out_dvid = [] for dv_id in self.__wing_param.get_dv_id_list(): if dv_id.startswith('rtwist0'): out_dvid.append('rtwist') elif not dv_id.startswith('rtwist'): out_dvid.append(dv_id) return tuple(out_dvid), tuple(self.__DLLM.get_F_list_names()) def provideJ(self): """Calculate the Jacobian according inputs and outputs""" self.__DLLM.run_adjoint() return np.array(self.__DLLM.get_dF_list_dchi()) def get_F_list(self): return self.__DLLM.get_F_list() def get_F_list_names(self): return self.__DLLM.get_F_list_names() def get_dv_array(self): return self.__wing_param.get_dv_array() def get_dv_id_list(self): return self.__wing_param.get_dv_id_list() def get_dv_info_list(self): return self.__wing_param.get_dv_info_list() def get_dv_value(self,dv_id): return self.get_dv_array()[self.get_dv_id_list().index(dv_id)]
wing_param.set_value('root_height',1.28) wing_param.set_value('break_height',0.97) wing_param.set_value('tip_height',0.33) wing_param.convert_to_design_variable('span',(10.,50.)) wing_param.convert_to_design_variable('sweep',(0.,40.)) wing_param.convert_to_design_variable('break_percent',(20.,40.)) wing_param.convert_to_design_variable('root_chord',(5.,7.)) wing_param.convert_to_design_variable('break_chord',(3.,5.)) wing_param.convert_to_design_variable('tip_chord',(1.,2.)) wing_param.convert_to_design_variable('root_height',(1.,1.5)) wing_param.convert_to_design_variable('break_height',(0.8,1.2)) wing_param.convert_to_design_variable('tip_height',(0.2,0.5)) wing_param.build_linear_airfoil(OC, AoA0=-2., Cm0=-0.1, set_as_ref=True) wing_param.build_airfoils_from_ref() wing_param.update() print wing_param DLLM = DLLMTargetLift('TLift',wing_param,OC) DLLM.set_target_Lift(769200.) #DLLM.set_method('scipy') DLLM.run_direct() DLLM.run_post() DLLM.run_adjoint() dF_list_dchi=DLLM.get_dF_list_dchi() # for i,name in enumerate(DLLM.get_F_list_names()): # print 'gradients for '+name+' = ',dF_list_dchi[i]
OC.set_T0_deg(15.) OC.set_P0(101325.) OC.set_humidity(0.) OC.compute_atmosphere() wing_param=Wing_Elliptic('elliptic_wing',n_sect=20) wing_param.import_BC_from_file('input_parameters.par') wing_param.build_linear_airfoil(OC, AoA0=0.0, set_as_ref=True) wing_param.build_airfoils_from_ref() wing_param.update() print wing_param x0=wing_param.get_dv_array() DLLM = DLLMTargetLift('Simple',wing_param,OC) DLLM.set_target_Lift(769200.) DLLM.run_direct() DLLM.run_post() Ref_F = DLLM.get_F_list() def f(x): wing_param.update_from_x_list(x) DLLM = DLLMTargetLift('Simple',wing_param,OC) DLLM.set_target_Lift(769200.) DLLM.run_direct() DLLM.run_post() func=DLLM.get_F_list()/Ref_F return func def df(x):
class DLLMWrapper(): ERROR_MSG = 'ERROR in DLLMWrapper.' WARNING_MSG = 'WARNING in DLLMWrapper.' POS_SOLVER = ['Solver','TargetCl','TargetLift'] POS_FMT = ['list','numpy'] def __init__(self, tag): """ Wrapper for the DLLM solver """ self.__tag = tag self.__OC = None # Operating condition self.__wing_param = None # Wing_param class self.__DLLM_solver = None # DLLM Solver class self.__config_dict = None self.__out_format = 'list' self.__grad_format = 'list' self.__AoA_id = 'AoA' self.__F_list = None self.__F_list_grad = None #-- Accessors def get_OC(self): return self.__OC def get_wing_param(self): return self.__wing_param def get_DLLM_solver(self): return self.__DLLM_solver def get_tags_x0_and_bounds(self): tags=self.__wing_param.get_dv_id_list() x0=self.__wing_param.get_dv_array() bounds=self.__wing_param.get_bounds_array() return tags,x0,bounds def get_x0_and_bounds(self): x0=self.__wing_param.get_dv_array() bounds=self.__wing_param.get_bounds_array() return x0,bounds def get_x0(self): return self.get_x() def get_x(self): x=self.__wing_param.get_dv_array() return x def get_F_list_names(self): return self.__DLLM_solver.get_F_list_names() def get_F_list(self): return self.__F_list def get_F_list_grad(self): return self.__F_list_grad def get_F_list_and_grad(self): return self.__F_list, self.__F_list_grad #-- Setters def set_AoA_id(self, AoA_id): self.__AoA_id = AoA_id def set_out_format(self, format): WARNING_MSG=self.WARNING_MSG+'set_out_format: ' if format not in self.POS_FMT: print WARNING_MSG+'format = '+str(format)+' not in '+str(self.POS_FMT)+'. Set to default out format = list' format='list' self.__out_format = format def set_grad_format(self, format): WARNING_MSG=self.WARNING_MSG+'set_grad_format: ' if format not in self.POS_FMT: print WARNING_MSG+'format = '+str(format)+' not in '+str(self.POS_FMT)+'. Set to default grad format = list' format='list' self.__grad_format = format #-- Public methods def configure(self, config_dict): self.__config_dict = config_dict self.__config_OC() self.__config_param() self.__config_DLLM() def run(self, x): self.__wing_param.update_from_x_list(x) self.__DLLM_solver.set_wing_param(self.__wing_param) F_list=self.analysis() return F_list def run_grad(self, x): self.__wing_param.update_from_x_list(x) self.__DLLM_solver.set_wing_param(self.__wing_param) F_list_grad=self.analysis_grad() return F_list_grad def run_and_grad(self, x): self.__wing_param.update_from_x_list(x) self.__DLLM_solver.set_wing_param(self.__wing_param) F_list,F_list_grad = self.analysis_and_grad() return F_list,F_list_grad def analysis(self): print self.__wing_param self.__DLLM_solver.run_direct() self.__DLLM_solver.run_post() F_list = self.__DLLM_solver.get_F_list() if self.__out_format == 'list': F_list=F_list.tolist() self.__F_list = F_list return F_list def analysis_grad(self): print self.__wing_param self.__DLLM_solver.run_direct() self.__DLLM_solver.run_post() self.__DLLM_solver.run_adjoint() F_list_grad=self.__DLLM_solver.get_dF_list_dchi() if self.__grad_format == 'numpy': F_list_grad=numpy.array(F_list_grad) self.__F_list_grad = F_list_grad return F_list_grad def analysis_and_grad(self): print self.__wing_param self.__DLLM_solver.run_direct() self.__DLLM_solver.run_post() self.__DLLM_solver.run_adjoint() F_list = self.__DLLM_solver.get_F_list() F_list_grad=self.__DLLM_solver.get_dF_list_dchi() if self.__out_format == 'list': F_list=F_list.tolist() if self.__grad_format == 'numpy': F_list_grad=numpy.array(F_list_grad) self.__F_list = F_list self.__F_list_grad = F_list_grad return F_list,F_list_grad def export_results(self): fid=open(self.__tag+'.res.dat','w') res=[self.__F_list,self.__F_list_grad] cPickle.dump(res, fid) fid.close() def import_results(self): fid=open(self.__tag+'.res.dat','r') res=cPickle.load(fid) fid.close() self.__F_list=res[0] self.__F_list_grad=res[1] #-- Private methods def __config_OC(self): self.__OC = OperatingCondition(self.__tag+'.OC') self.__OC.config_from_dict(self.__config_dict) def __config_param(self): input_keys=self.__config_dict.keys() geom_type_key=self.__tag+'.param.geom_type' if geom_type_key in input_keys: geom_type = self.__config_dict[geom_type_key] else: geom_type= 'Broken' n_sect_key=self.__tag+'.param.n_sect' if n_sect_key in input_keys: n_sect = self.__config_dict[n_sect_key] else: n_sect = 20 self.__wing_param = Wing_param(self.__tag+'.param',geom_type=geom_type,n_sect=n_sect) self.__wing_param.set_AoA_id(self.__AoA_id) self.__wing_param.config_from_dict(self.__OC, self.__config_dict) def __config_DLLM(self): WARNING_MSG=self.WARNING_MSG+'__config_DLLM: ' input_keys=self.__config_dict.keys() type_key = self.__tag+'.DLLM.type' type = self.__config_dict[type_key] if type not in self.POS_SOLVER: print WARNING_MSG+'solver_type = '+str(solve_type)+' not in '+str(self.POS_SOLVER)+'. Set to default solver_type = Solver' type='Solver' if type == 'Solver': self.__DLLM_solver = DLLMSolver(self.__tag,self.__wing_param,self.__OC) elif type == 'TargetCl': self.__DLLM_solver = DLLMTargetCl(self.__tag,self.__wing_param,self.__OC) target_Cl_key = self.__tag+'.DLLM.target_Cl' target_Cl = self.__config_dict[target_Cl_key] self.__DLLM_solver.set_target_Cl(target_Cl) elif type == 'TargetLift': self.__DLLM_solver = DLLMTargetLift(self.__tag,self.__wing_param,self.__OC) target_Lift_key = self.__tag+'.DLLM.target_Lift' target_Lift = self.__config_dict[target_Lift_key] self.__DLLM_solver.set_target_Lift(target_Lift) method_key = self.__tag+'.DLLM.method' if method_key in input_keys: method = self.__config_dict[method_key] self.__DLLM_solver.set_method(method) relax_factor_key = self.__tag+'.DLLM.relax_factor' if relax_factor_key in input_keys: relax_factor = self.__config_dict[relax_factor_key] self.__DLLM_solver.set_relax_factor(relax_factor) stop_residual_key = self.__tag+'.DLLM.stop_residual' if stop_residual_key in input_keys: stop_residual = self.__config_dict[stop_residual_key] self.__DLLM_solver.set_stop_residual(stop_residual) max_iterations_key = self.__tag+'.DLLM.max_iterations' if max_iterations_key in input_keys: max_iterations = self.__config_dict[max_iterations_key] self.__DLLM_solver.set_max_iterations(max_iterations) gamma_file_name_key = self.__tag+'.DLLM.gamma_file_name' if gamma_file_name_key in input_keys: gamma_file_name = self.__config_dict[gamma_file_name_key] self.__DLLM_solver.set_gamma_file_name(gamma_file_name) F_list_names_key = self.__tag+'.DLLM.F_list_names' if F_list_names_key in input_keys: F_list_names = self.__config_dict[F_list_names_key] self.__DLLM_solver.set_F_list_names(F_list_names)
class DLLMWrapper(): ERROR_MSG = 'ERROR in DLLMWrapper.' WARNING_MSG = 'WARNING in DLLMWrapper.' POS_SOLVER = ['Solver','TargetCl','TargetLift'] POS_FMT = ['list','numpy'] def __init__(self, tag): """ Wrapper for the DLLM solver """ self.__tag = tag self.__OC = None # Operating condition self.__wing_param = None # Wing_param class self.__DLLM_solver = None # DLLM Solver class self.__config_dict = None self.__out_format = 'list' self.__grad_format = 'list' self.__AoA_id = 'AoA' self.__F_list = None self.__F_list_grad = None #-- Accessors def get_OC(self): return self.__OC def get_wing_param(self): return self.__wing_param def get_DLLM_solver(self): return self.__DLLM_solver def get_tags_x0_and_bounds(self): tags=self.__wing_param.get_dv_id_list() x0=self.__wing_param.get_dv_array() bounds=self.__wing_param.get_bounds_array() return tags,x0,bounds def get_x0_and_bounds(self): x0=self.__wing_param.get_dv_array() bounds=self.__wing_param.get_bounds_array() return x0,bounds def get_x0(self): return self.get_x() def get_x(self): x=self.__wing_param.get_dv_array() return x def get_F_list_names(self): return self.__DLLM_solver.get_F_list_names() def get_F_list(self): return self.__F_list def get_F_list_grad(self): return self.__F_list_grad def get_F_list_and_grad(self): return self.__F_list, self.__F_list_grad #-- Setters def set_AoA_id(self, AoA_id): self.__AoA_id = AoA_id def set_out_format(self, format): WARNING_MSG=self.WARNING_MSG+'set_out_format: ' if format not in self.POS_FMT: print WARNING_MSG+'format = '+str(format)+' not in '+str(self.POS_FMT)+'. Set to default out format = list' format='list' self.__out_format = format def set_grad_format(self, format): WARNING_MSG=self.WARNING_MSG+'set_grad_format: ' if format not in self.POS_FMT: print WARNING_MSG+'format = '+str(format)+' not in '+str(self.POS_FMT)+'. Set to default grad format = list' format='list' self.__grad_format = format #-- Public methods def configure(self, config_dict): self.__config_dict = config_dict self.__config_OC() self.__config_param() self.__config_DLLM() def run(self, x): self.__wing_param.update_from_x_list(x) self.__DLLM_solver.set_wing_param(self.__wing_param) F_list=self.analysis() return F_list def run_grad(self, x): self.__wing_param.update_from_x_list(x) self.__DLLM_solver.set_wing_param(self.__wing_param) F_list_grad=self.analysis_grad() return F_list_grad def run_and_grad(self, x): self.__wing_param.update_from_x_list(x) self.__DLLM_solver.set_wing_param(self.__wing_param) F_list,F_list_grad = self.analysis_and_grad() return F_list,F_list_grad def analysis(self): print self.__wing_param self.__DLLM_solver.run_direct() self.__DLLM_solver.run_post() F_list = self.__DLLM_solver.get_F_list() if self.__out_format == 'list': F_list=F_list.tolist() self.__F_list = F_list return F_list def analysis_grad(self): print self.__wing_param self.__DLLM_solver.run_direct() self.__DLLM_solver.run_post() self.__DLLM_solver.run_adjoint() F_list_grad=self.__DLLM_solver.get_dF_list_dchi() if self.__grad_format == 'numpy': F_list_grad=numpy.array(F_list_grad) self.__F_list_grad = F_list_grad return F_list_grad def analysis_and_grad(self): print self.__wing_param self.__DLLM_solver.run_direct() self.__DLLM_solver.run_post() self.__DLLM_solver.run_adjoint() F_list = self.__DLLM_solver.get_F_list() F_list_grad=self.__DLLM_solver.get_dF_list_dchi() if self.__out_format == 'list': F_list=F_list.tolist() if self.__grad_format == 'numpy': F_list_grad=numpy.array(F_list_grad) self.__F_list = F_list self.__F_list_grad = F_list_grad return F_list,F_list_grad def export_results(self): fid=open(self.__tag+'.res.dat','w') res=[self.__F_list,self.__F_list_grad] cPickle.dump(res, fid) fid.close() def import_results(self): fid=open(self.__tag+'.res.dat','r') res=cPickle.load(fid) fid.close() self.__F_list=res[0] self.__F_list_grad=res[1] #-- Private methods def __config_OC(self): self.__OC = OperatingCondition(self.__tag+'.OC') self.__OC.config_from_dict(self.__config_dict) def __config_param(self): input_keys=self.__config_dict.keys() geom_type_key=self.__tag+'.param.geom_type' if geom_type_key in input_keys: geom_type = self.__config_dict[geom_type_key] else: geom_type= 'Broken' n_sect_key=self.__tag+'.param.n_sect' if n_sect_key in input_keys: n_sect = self.__config_dict[n_sect_key] else: n_sect = 20 self.__wing_param = Wing_param(self.__tag+'.param',geom_type=geom_type,n_sect=n_sect) self.__wing_param.set_AoA_id(self.__AoA_id) self.__wing_param.config_from_dict(self.__OC, self.__config_dict) def __config_DLLM(self): WARNING_MSG=self.WARNING_MSG+'__config_DLLM: ' input_keys=self.__config_dict.keys() type_key = self.__tag+'.DLLM.type' type = self.__config_dict[type_key] if type not in self.POS_SOLVER: print WARNING_MSG+'solver_type = '+str(type)+' not in '+str(self.POS_SOLVER)+'. Set to default solver_type = Solver' type='Solver' if type == 'Solver': self.__DLLM_solver = DLLMSolver(self.__tag,self.__wing_param,self.__OC) elif type == 'TargetCl': self.__DLLM_solver = DLLMTargetCl(self.__tag,self.__wing_param,self.__OC) target_Cl_key = self.__tag+'.DLLM.target_Cl' target_Cl = self.__config_dict[target_Cl_key] self.__DLLM_solver.set_target_Cl(target_Cl) elif type == 'TargetLift': self.__DLLM_solver = DLLMTargetLift(self.__tag,self.__wing_param,self.__OC) target_Lift_key = self.__tag+'.DLLM.target_Lift' target_Lift = self.__config_dict[target_Lift_key] self.__DLLM_solver.set_target_Lift(target_Lift) method_key = self.__tag+'.DLLM.method' if method_key in input_keys: method = self.__config_dict[method_key] self.__DLLM_solver.set_method(method) relax_factor_key = self.__tag+'.DLLM.relax_factor' if relax_factor_key in input_keys: relax_factor = self.__config_dict[relax_factor_key] self.__DLLM_solver.set_relax_factor(relax_factor) stop_residual_key = self.__tag+'.DLLM.stop_residual' if stop_residual_key in input_keys: stop_residual = self.__config_dict[stop_residual_key] self.__DLLM_solver.set_stop_residual(stop_residual) max_iterations_key = self.__tag+'.DLLM.max_iterations' if max_iterations_key in input_keys: max_iterations = self.__config_dict[max_iterations_key] self.__DLLM_solver.set_max_iterations(max_iterations) gamma_file_name_key = self.__tag+'.DLLM.gamma_file_name' if gamma_file_name_key in input_keys: gamma_file_name = self.__config_dict[gamma_file_name_key] self.__DLLM_solver.set_gamma_file_name(gamma_file_name) F_list_names_key = self.__tag+'.DLLM.F_list_names' if F_list_names_key in input_keys: F_list_names = self.__config_dict[F_list_names_key] self.__DLLM_solver.set_F_list_names(F_list_names)