def setUp( self ): '''Test combination ULS/SLS. ''' lc_list = [ LC( name = 'G1', category = 'dead-load' ), LC( name = 'G2', category = 'additional dead-load' ), LC( name = 'Q1', category = 'imposed-load', exclusive_to = ['Q2', 'Q3'], psi_0 = 0.7, psi_1 = 0.5, psi_2 = 0.2 ), LC( name = 'Q2', category = 'imposed-load', exclusive_to = ['Q1', 'Q3'], psi_0 = 0.7, psi_1 = 0.5, psi_2 = 0.2 ), LC( name = 'Q3', category = 'imposed-load', exclusive_to = ['Q1', 'Q2'], psi_0 = 0.7, psi_1 = 0.5, psi_2 = 0.2 ), LC( name = 'Q4', category = 'imposed-load', exclusive_to = ['Q5'], psi_0 = 0.7, psi_1 = 0.5, psi_2 = 0.2 ), LC( name = 'Q5', category = 'imposed-load', exclusive_to = ['Q4'], psi_0 = 0.7, psi_1 = 0.5, psi_2 = 0.2 ), LC( name = 'Q6', category = 'imposed-load', exclusive_to = [], psi_0 = 0.7, psi_1 = 0.5, psi_2 = 0.2 ), ] self.lct_list = [ LCCTableULS( show_lc_characteristic = False, lc_list = lc_list ), LCCTableSLS( show_lc_characteristic = False, lc_list = lc_list, combination_SLS = 'perm' ), LCCTableSLS( show_lc_characteristic = False, lc_list = lc_list, combination_SLS = 'freq' ), LCCTableSLS( show_lc_characteristic = False, lc_list = lc_list, combination_SLS = 'rare' ) ] self.n_combi_list = [ 188, 24, 47, 47 ]
### compressive strength [kN/m] ### # C90/105, f_cd= 51 MN/m^2, high utilised = 0.45 m n_Rdc = 22950. print 'design values calculated by hand for strength characteristics' # LCCTable for imposed loads (without temperature) # lct_Q = LCCTableULS( data_dir=data_dir, reader_type='InfoCAD', lc_list=lc_list_Q, strength_characteristics={ 'n_0_Rdt': n_0_Rdt, 'm_0_Rd': m_0_Rd, 'n_Rdc': n_Rdc, 'n_90_Rdt': n_90_Rdt, 'm_90_Rd': m_90_Rd }, # simplification: use the minimum value for k_alpha # on the resistance side k_alpha_min=False, # show_lc_characteristic = True ) #-------------------------------------------------------------- # 'combi_arr': array with indices of all loading case combinations #-------------------------------------------------------------- # print 'lct_Q.combi_arr', lct_Q.combi_arr.shape, '\n' # np.savetxt('combi_arr_wo_temp_LCs', lct_Q.combi_arr, delimiter=';')
### compressive strength [kN/m] ### # C90/105, f_cd= 51 MN/m^2, high utilised = 0.45 m n_Rdc = 22950.0 print "design values calculated by hand for strength characteristics" # LCCTable for imposed loads (without temperature) # lct_Q = LCCTableULS( data_dir=data_dir, reader_type="InfoCAD", lc_list=lc_list_Q, strength_characteristics={ "n_0_Rdt": n_0_Rdt, "m_0_Rd": m_0_Rd, "n_Rdc": n_Rdc, "n_90_Rdt": n_90_Rdt, "m_90_Rd": m_90_Rd, }, # simplification: use the minimum value for k_alpha # on the resistance side k_alpha_min=False, # show_lc_characteristic = True ) # -------------------------------------------------------------- # 'combi_arr': array with indices of all loading case combinations # -------------------------------------------------------------- # print "lct_Q.combi_arr", lct_Q.combi_arr.shape, "\n" # np.savetxt('combi_arr_wo_temp_LCs', lct_Q.combi_arr, delimiter=';')
# compressive strength [kN/m] # (C3-B2-HF-1-140-5 with characteristic value; f_ck = 100 MPa) # (design value; f_cd = 56,7 MPa) # n_Rdc = 22680. # = 56,7 MPa * (100 cm * 40 cm) * 0.1 print 'design values calculated by hand for strength characteristics' # LCCTable for imposed loads (without temperature) # lct_Q = LCCTableULS(data_dir=data_dir, reader_type='InfoCAD', lc_list=lc_list_Q, strength_characteristics={'n_0_Rdt': n_0_Rdt, 'm_0_Rd': m_0_Rd, 'n_Rdc': n_Rdc, 'n_90_Rdt': n_90_Rdt, 'm_90_Rd': m_90_Rd}, # simplification: use the minimum value for k_alpha # on the resistance side k_alpha_min=False, # show_lc_characteristic = True ) #-------------------------------------------------------------- # 'combi_arr': array with indices of all loading case combinations #-------------------------------------------------------------- # print 'lct_Q.combi_arr', lct_Q.combi_arr.shape, '\n' # np.savetxt('combi_arr_wo_temp_LCs', lct_Q.combi_arr, delimiter=';') #-------------------------------------------------------------- # nm-interaction plot (normal force - bending moment) #--------------------------------------------------------------
### compressive strength [kN/m] ### # C90/105, f_cd= 51 MN/m^2, high utilised = 0.45 m n_Rdc = 22950. print 'design values calculated by hand for strength characteristics' # LCCTable for imposed loads (without temperature) # lct_Q = LCCTableULS(data_dir=data_dir, reader_type='InfoCAD', lc_list=lc_list_Q, strength_characteristics={'n_0_Rdt': n_0_Rdt, 'm_0_Rd': m_0_Rd, 'n_Rdc': n_Rdc, 'n_90_Rdt': n_90_Rdt, 'm_90_Rd': m_90_Rd}, # simplification: use the minimum value for k_alpha # on the resistance side k_alpha_min=False, # show_lc_characteristic = True ) #-------------------------------------------------------------- # 'combi_arr': array with indices of all loading case combinations #-------------------------------------------------------------- # print 'lct_Q.combi_arr', lct_Q.combi_arr.shape, '\n' # np.savetxt('combi_arr_wo_temp_LCs', lct_Q.combi_arr, delimiter=';') #-------------------------------------------------------------- # nm-interaction plot (normal force - bending moment) #--------------------------------------------------------------
# n_Rdc = 2200 # C55/67 # bending strength [kNm/m] # m_0_Rd = m_90_Rd = 15.5 # = 3.1 / 0.20 ### M_Rm = 3.1 kNm (mean value) # LCCTable for imposed loads # lct = LCCTableULS(data_dir=data_dir, reader_type='InfoCAD', lc_list=lc_list, show_lc_characteristic=False, strength_characteristics={'n_0_Rdt': n_0_Rdt, 'm_0_Rd': m_0_Rd, 'n_Rdc': n_Rdc, 'n_90_Rdt': n_90_Rdt, 'm_90_Rd': m_90_Rd}, # NO simplification used for 'k_alpha' on the resistance # side k_alpha_min=False, ) #-------------------------------------------------------------- # 'combi_arr': array with indices of all loading case combinations #-------------------------------------------------------------- # # print 'lct.combi_arr', lct.combi_arr.shape # np.savetxt('combi_arr_wo_temp_LCs', lct_Q.combi_arr, delimiter = ';') #--------------------------------------------------------------
m_90_Rd = 1717. # compressive strength [kN/m] # C90/105, f_cd= 51 MN/m^2, high utilised = 0.295 m n_Rdc = 15045. print 'design values calculated by hand for strength characteristics' # LCCTable for imposed loads (without temperature) # lct_Q = LCCTableULS(data_dir=data_dir, reader_type='InfoCAD', lc_list=lc_list_Q, strength_characteristics={'n_0_Rdt': n_0_Rdt, 'm_0_Rd': m_0_Rd, 'n_Rdc': n_Rdc, 'n_90_Rdt': n_90_Rdt, 'm_90_Rd': m_90_Rd}, # simplification: use the minimum value for k_alpha # on the resistance side k_alpha_min=False, # show_lc_characteristic = True ) #-------------------------------------------------------------- # 'combi_arr': array with indices of all loading case combinations #-------------------------------------------------------------- # ## print 'lct_Q.combi_arr', lct_Q.combi_arr.shape, '\n' # np.savetxt('combi_arr_wo_temp_LCs', lct_Q.combi_arr, delimiter=';') #-------------------------------------------------------------- # nm-interaction plot (normal force - bending moment) #--------------------------------------------------------------
# compressive strength [kN/m] # n_Rdc = 2200 # C55/67 # bending strength [kNm/m] # m_0_Rd = m_90_Rd = 0.833333 * 8.3 # = 1.66 (=M_Rd) / 0.20 ### M_Rm = 3.1 kNm (mean value) #-------------------------------------------------------- # ULS evaluation #-------------------------------------------------------- lct = LCCTableULS(data_dir=data_dir, data_filter=remove_midpoints, lc_list=lc_list, show_lc_characteristic=True, strength_characteristics={'n_0_Rdt' : n_0_Rdt, 'm_0_Rd':m_0_Rd, 'n_Rdc' : n_Rdc, 'n_90_Rdt' : n_90_Rdt, 'm_90_Rd':m_90_Rd}, k_alpha_min=False, # NO simplification used for 'k_alpha' on the resistance side ) #-------------------------------------------------------------- # 'combi_arr': array with indices of all loading case combinations #-------------------------------------------------------------- # print 'lct.combi_arr', lct.combi_arr.shape, '\n' # np.savetxt('combi_arr_LC1-12', lct.combi_arr, delimiter=';') #-------------------------------------------------------------- # nm-interaction plot (normal force - bending moment) #-------------------------------------------------------------- #
# (design value; f_cd = 37,5 MPa) # n_Rdc = 15000.0 # = 37,5 MPa * (100 cm * 40 cm) * 0.1 print "design values calculated by hand for strength characteristics" # LCCTable for imposed loads (without temperature) # lct_Q = LCCTableULS( data_dir=data_dir, reader_type="InfoCAD", lc_list=lc_list_Q, strength_characteristics={ "n_0_Rdt": n_0_Rdt, "m_0_Rd": m_0_Rd, "n_Rdc": n_Rdc, "n_90_Rdt": n_90_Rdt, "m_90_Rd": m_90_Rd, }, # simplification: use the minimum value for k_alpha # on the resistance side k_alpha_min=False, # show_lc_characteristic = True ) # -------------------------------------------------------------- # 'combi_arr': array with indices of all loading case combinations # -------------------------------------------------------------- # print "lct_Q.combi_arr", lct_Q.combi_arr.shape, "\n" # np.savetxt('combi_arr_wo_temp_LCs', lct_Q.combi_arr, delimiter=';')
), ] #-------------------------------------------------------------- if do == 'dimensioning': # LCCTable for imposed loads (without temperature) # lct_Q = LCCTableULS(data_dir=data_dir, reader_type='InfoCAD', data_filter=remove_support_elems, lc_list=lc_list_Q, # simplification: use the minimum value for k_alpha # on the resistance side k_alpha_min=True, # show_lc_characteristic = True ) # LCCTable for temperature loading cases only # lct_T = LCCTableULS(data_dir=data_dir, reader_type='InfoCAD', data_filter=remove_support_elems, lc_list=lc_list_T, # simplification: use the minimum value for k_alpha # on the resistance side k_alpha_min=True, # show_lc_characteristic = True
psi_0=1.0, psi_1=0.5, psi_2=0.0), ] #-------------------------------------------------------------- if do == 'dimensioning': # LCCTable for imposed loads (without temperature) # lct_Q = LCCTableULS( data_dir=data_dir, reader_type='InfoCAD', data_filter=remove_support_elems, lc_list=lc_list_Q, # simplification: use the minimum value for k_alpha # on the resistance side k_alpha_min=True, # show_lc_characteristic = True ) # LCCTable for temperature loading cases only # lct_T = LCCTableULS( data_dir=data_dir, reader_type='InfoCAD', data_filter=remove_support_elems, lc_list=lc_list_T, # simplification: use the minimum value for k_alpha # on the resistance side k_alpha_min=True,
### compressive strength [kN/m] ### #C90/105, f_cd= 51 MN/m^2, high utilised = 0.35 m n_Rdc = 18069. print 'design values calculated by hand for strength characteristics' # LCCTable for imposed loads (without temperature) # lct_Q = LCCTableULS(data_dir=data_dir, reader_type='InfoCAD', lc_list=lc_list_Q, strength_characteristics={'n_0_Rdt': n_0_Rdt, 'm_0_Rd': m_0_Rd, 'n_Rdc': n_Rdc, 'n_90_Rdt': n_90_Rdt, 'm_90_Rd': m_90_Rd}, # simplification: use the minimum value for k_alpha # on the resistance side k_alpha_min=False, # show_lc_characteristic = True ) #-------------------------------------------------------------- # 'combi_arr': array with indices of all loading case combinations #-------------------------------------------------------------- # print 'lct_Q.combi_arr', lct_Q.combi_arr.shape, '\n' # np.savetxt('combi_arr_wo_temp_LCs', lct_Q.combi_arr, delimiter=';') #-------------------------------------------------------------- # nm-interaction plot (normal force - bending moment) #--------------------------------------------------------------
# bending strength [kNm/m] # m_0_Rd = m_90_Rd = (3.5 * 0.46 / 4.) / 0.10 # [kNm/m] # compressive strength [kN/m] # (design value; f_cd = 37,5 MPa) # n_Rdc = 750. # = 37,5 MPa * (100 cm * 2 cm) * 0.1 lct = LCCTableULS(data_dir=data_dir, reader_type='InfoCAD', # data_filter=remove_support_elems, lc_list=lc_list, strength_characteristics={'n_0_Rdt': n_0_Rdt, 'm_0_Rd': m_0_Rd, 'n_Rdc': n_Rdc, 'n_90_Rdt': n_90_Rdt, 'm_90_Rd': m_90_Rd}, # NO simplification used for 'k_alpha' on the # resistance side k_alpha_min=False, show_lc_characteristic=True ) #-------------------------------------------------------------- # 'combi_arr': array with indices of all loading case combinations #-------------------------------------------------------------- # # print 'lct.combi_arr', lct.combi_arr.shape # np.savetxt('combi_arr_wo_temp_LCs', lct.combi_arr, delimiter=';') #--------------------------------------------------------------
# # LC14_3-4: # LC(name = 'EQ_3-4', category = 'imposed-load', # file_name = 'LC14_3-4.csv', # gamma_unf = 1.00, # exclusive_to = ['EQ_1-2'], # psi_0 = 1.0, # comment = 'earthquake' # ), # ] #-------------------------------------------------------- lct = LCCTableULS(data_dir = data_dir, data_filter = remove_midpoints, lc_list = lc_list, show_lc_characteristic = False ) # lct.configure_traits() lct.plot_assess_value() # lct.plot_n_tex() # lct = LCCTableSLS( data_dir = data_dir, # data_filter = remove_midpoints, # lc_list = lc_list, # combination_SLS = 'rare', ## combination_SLS = 'freq', ## combination_SLS = 'perm', ## show_lc_characteristic = True # )