# Structural age service_time = np.arange(START_AGE+TIME_INTERVAL,END_AGE+TIME_INTERVAL,TIME_INTERVAL) comp_type_list = ['flexure', 'shear', 'deck'] icorr_mean_list = [1, 1, 1] # construct virgin component str_yr_list = [0, 0, 0] cost_list = [] virgin_component_list = [] for comp_type,str_yr,icorr_mean in zip(comp_type_list, str_yr_list, icorr_mean_list): virgin_component = Component(comp_type, maintain_tag=False, str_yr=str_yr) resistance_mean,resistance_cov,cost = simpleCorrosionLHS(comp_type, service_time, icorr_mean, str_yr) virgin_component.setServiceTime(service_time) virgin_component.setResistanceMean(resistance_mean) virgin_component.setResistanceCov(resistance_cov) virgin_component.setCESampling(NUM_COMPONENT, NUM_ADAPTATION, NUM_PRE_SMP, NUM_MAIN_SMP) virgin_component_list.append(virgin_component) cost_list.append(cost) # strengthened component str_yr_list = [0, 50, 0] cost_list = [] component_list = [] for comp_type,str_yr,icorr_mean in zip(comp_type_list, str_yr_list, icorr_mean_list): component = Component(comp_type, maintain_tag=False, str_yr=str_yr) resistance_mean,resistance_cov,cost = simpleCorrosionLHS(comp_type, service_time, icorr_mean, str_yr) component.setServiceTime(service_time) component.setResistanceMean(resistance_mean) component.setResistanceCov(resistance_cov) component.setCESampling(NUM_COMPONENT, NUM_ADAPTATION, NUM_PRE_SMP, NUM_MAIN_SMP) component_list.append(component)