In engineering and physics, a spring system or spring network is a model of physics described as a graph with a position at each vertex and a spring of given stiffness and length along each edge. This generalizes Hooke's law to higher dimensions. This simple model can be used to solve the pose of static systems from crystal lattice to springs. A spring system can be thought of as the simplest case of the finite element method for solving problems in statics. Assuming linear springs and small deformation (or restricting to one-dimensional motion) a spring system can be cast as a (possibly overdetermined) system of linear equations or equivalently as an energy minimization problem.
Known spring lengths Consider the simple case of three nodes, in one dimension x = [ x 1 x 2 x 3 ] {\displaystyle \mathbf {x} ={\begin{bmatrix}x_{1}\\x_{2}\\x_{3}\end{bmatrix}}} , connected by two springs. If the nominal lengths, L, of the springs are known to be 1 and 2 units respectively, i.e. L = [ 1 2 ] {\displaystyle \mathbf {L} ={\begin{bmatrix}1\\2\end{bmatrix}}} , then the system can be solved as follows: The stretching of the two springs is given as a function of the positions of the nodes by
Δ L = B ⊤ x − L = [ − 1 1 0 0 − 1 1 ] x − L {\displaystyle \Delta \mathbf {L} =B^{\top }\mathbf {x} -\mathbf {L} ={\begin{bmatrix}-1&1&0\\0&-1&1\end{bmatrix}}\mathbf {x} -\mathbf {L} }
where B ⊤ {\displaystyle B^{\top }} is the matrix transpose of the oriented incidence matrix
B = [ − 1 0 1 − 1 0 1 ] , {\displaystyle B={\begin{bmatrix}-1&0\\1&-1\\0&1\end{bmatrix}},}
relating each degree of freedom to the direction each spring pulls on it. The forces on the springs are
F springs = − W Δ L = − W ( B ⊤ x − L ) = − W B ⊤ x + W L {\displaystyle F_{\text{springs}}=-W\Delta \mathbf {L} =-W(B^{\top }\mathbf {x} -\mathbf {L} )=-WB^{\top }\mathbf {x} +W\mathbf {L} }
where W is a diagonal matrix giving the stiffness of every spring. Then the force on the nodes is given by left multiplying by B {\displaystyle B} , which we set to zero to find equilibrium:
F nodes = − B W B ⊤ x + B W L = 0 {\displaystyle F_{\text{nodes}}=-BWB^{\top }\mathbf {x} +BW\mathbf {L} =0}
which gives the linear equation:
B W B ⊤ x = B W L {\displaystyle BWB^{\top }\mathbf {x} =BW\mathbf {L} } . Now, the matrix B W B ⊤ {\displaystyle BWB^{\top }} is singular, because all solutions are equivalent up to rigid-body translation. Let us prescribe a Dirichlet boundary condition, e.g., x 1 = 2 {\displaystyle x_{1}=2} . As an example, let W be the identity matrix then
… excerpt ends here. Continue reading the full article.
