equilibrium_matrix

compas.numerical.equilibrium_matrix(C, xyz, free, rtype='array')[source]

Construct the equilibrium matrix of a structural system.

Parameters:

Carray-like: Connectivity matrix (m x n).
xyzarray-like: Array of vertex coordinates (n x 3).
freelist: The index values of the free vertices.
rtype{‘array’, ‘csc’, ‘csr’, ‘coo’, ‘list’}: Format of the result.

Returns:

array-like: Constructed equilibrium matrix.

Notes

Analysis of the equilibrium matrix reveals some of the properties of the structural system, its size is (2ni x m) where ni is the number of free or internal nodes. It is calculated by

\[\begin{split}\mathbf{E} = \left[ \begin{array}{c} \mathbf{C}^{\mathrm{T}}_{\mathrm{i}}\mathbf{U} \\[0.3em] \hline \\[-0.7em] \mathbf{C}^{\mathrm{T}}_{\mathrm{i}}\mathbf{V} \end{array} \right].\end{split}\]

The matrix of vertex coordinates is vectorised to speed up the calculations.

Examples

>>> C = connectivity_matrix([[0, 1], [0, 2], [0, 3]])
>>> xyz = [[0, 0, 1], [0, 1, 0], [-1, -1, 0], [1, -1, 0]]
>>> equilibrium_matrix(C, xyz, [0], rtype='array')
array([[ 0.,  1., -1.],
       [-1.,  1.,  1.]])