def pavement(self): bridge_w = self.bridge_inf['bridge'][1] - self.bridge_inf['bridge'][0] pave_w = bridge_w - sum(self.barrier[0]) barrier_f = sum(self.barrier[1]) * 26 pave_f = bridge_w * self.pave_thickness[ 0] * 26 + pave_w * self.pave_thickness[1] * 24 fz = -barrier_f - pave_f midas1.start_stld('PAVEMENT', self.mct_l) midas1.start_beam_load(self.mct_l) i_elems = self.elem[self.elem['what'] == 'beam'].index midas1.beam_load(i_elems, fz, mct_list=self.mct_l, group='pavement')
def brake_force(self, side=2): midas1.start_stld('BRAKE_FORCE', self.mct_l) midas1.start_beam_load(self.mct_l) bridge_w = self.bridge_inf['bridge'][1] - self.bridge_inf['bridge'][0] lane_num = kny_tool.get_lanes_from_w(bridge_w, side=side) if side == 2: brake_f = lane_num / 2 * 165 else: brake_f = lane_num * 165 brake_f *= 1 / 30 brake_elem = self.elem[self.elem['what'] == 'beam'].index midas1.beam_load(brake_elem, brake_f, side='X', offset_side='Z', offset_dis=1.2, mct_list=self.mct_l, group='brake_force')
def wind(self, h=5): bridge_w = self.bridge_inf['bridge'][1] - self.bridge_inf['bridge'][0] wind_beam = wind.BeamWind(0, [bridge_w, 1.6, 14], h, 30) beam_f = wind_beam.Fg / 1000 elem_wind = self.elem[self.elem['what'] == 'beam'].index midas1.start_stld('WIND_1', self.mct_l) midas1.start_beam_load(self.mct_l) midas1.beam_load(elem_wind, beam_f, mct_list=self.mct_l, group='wind_1', side='Y') midas1.start_stld('WIND_2', self.mct_l) midas1.start_beam_load(self.mct_l) midas1.beam_load(elem_wind, -beam_f, mct_list=self.mct_l, group='wind_2', side='Y')