Finite Element Analysis, Python, Sfepy, Motivation.
theta=math.atan((coord[el[1]][1]-coord[el[0]][1])/(coord[el[1]][0]-coord[el[0]][0]+1e-13))*180/math.pi k=PlaneTrussElementStiffness(E, A[i], L, theta) K=PlaneTrussAssemble(K,k,el[0],el[1]) i+=1 F = global_force bc=np.array([bc_node, bc_val]).T nd, nd=K.shape fdof=np.array([i for i in range(nd)]).T