Fixed geradin test case
@@ -27,7 +27,7 @@ def generate_fem_file(route, case_name, num_elem, num_node_elem=3):
y = (np.linspace(0, length, num_node))*np.sin(angle)
z = np.zeros((num_node,))
- structural_twist = np.zeros_like(x)
+ structural_twist = np.zeros((num_elem, num_node_elem))
frame_of_reference_delta = np.zeros((num_elem, num_node_elem, 3))
for ielem in range(num_elem):
@@ -64,7 +64,7 @@ def generate_fem_file(route, case_name, num_elem, num_node_elem=3):
# mass array
num_mass = 1
- m_bar = 0*100
+ m_bar = 0.
j = 10
base_mass = np.diag([m_bar, m_bar, m_bar, j, j, j])
mass = np.zeros((num_mass, 6, 6))
