コード例 #1
0
ファイル: __test__.py プロジェクト: rosoba/simvisage
    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 ]
コード例 #2
0
    ### 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=';')
コード例 #3
0
    ### 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=';')
コード例 #4
0
ファイル: lcc_hangar2.py プロジェクト: simvisage/simvisage
    # 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)
    #--------------------------------------------------------------
コード例 #5
0
    

    ### 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)
    #--------------------------------------------------------------
コード例 #6
0
    ### 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=';')
コード例 #7
0
ファイル: lcc_slabtest.py プロジェクト: rosoba/simvisage
    #
    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 = ';')

    #--------------------------------------------------------------
コード例 #8
0
    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)
    #--------------------------------------------------------------
コード例 #9
0
        # 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)
    #--------------------------------------------------------------
    #
コード例 #10
0
ファイル: lcc_hangar.py プロジェクト: simvisage/simvisage
    # (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=';')
コード例 #11
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,
                            #                          show_lc_characteristic = True
コード例 #12
0
               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,
コード例 #13
0
    

    ### 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)
    #--------------------------------------------------------------
コード例 #14
0
            # 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=';')

        #--------------------------------------------------------------
コード例 #15
0
#                 # 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
#                        )