Finite Element Analysis, Python, Sfepy, Motivation.
Element solved with Python, [4]: def FEA_u(coord, elcon, bc_node, bc_val, global_force, A, E): K=np.zeros(shape=(2*np.max(elcon)+2,2*np.max(elcon)+2)) i=0 for el in elcon: L=PlaneTrussElementLength(coord[el[0]][0],coord[el[0]][1],coord[el[1]][0],coord[el[1]][1])