return R
if __name__ == '__main__':
from beam import Beam
E = 210e9
nu = 0.3
b = 0.1
sect1 = Section([-0.05, 0.05], 210e9, 0.3, 0.1, 7800, 5 / 6, 'section-1')
sect2 = Section([-0.05, 0, 0.05], [210e9, 210e9], [0.3, 0.3], [0.1, 0.1],
[7800, 7800], [5 / 6, 5 / 6], 'section-2')
sect3 = Section([-0.05, 0, 0.05], [210e9, 210e9], [0.3, 0.3], [0.1, 0.1],
[7800, 15000], [5 / 6, 5 / 6], 'section-3')
beam = Beam(section=sect1)
beam.append_segment(1.0, sect2, P=0)
beam.nodes[0].free()
beam.nodes[1].free()
def get_segment(freq):
beam.set_frequency(freq)
beam.initialize()
segment = beam.segments[0]
return segment
segment = get_segment(100)
BL = segment.get_response_matrix(0)[3:]
BR = segment.get_response_matrix(1)[3:]
B = np.zeros([6, 6], dtype='complex')
B[:3] = BL
B[3:] = BR
