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=';')
#
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 = ';')
#--------------------------------------------------------------
# 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)
#--------------------------------------------------------------
#
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,
# 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=';')
#--------------------------------------------------------------
