NS_Truss_4_(v0.6) -- Sage

%auto

def map2 (f, l2) :
    return map (lambda y : map (f, y), l2)

def map_dict (f, d) :
    return dict (zip (d.keys(), map (f, d.values())))
    
def map2_dict (f, d) :
    return dict (zip (d.keys(), map2 (f, d.values())))

def filter_dict (f, d) : 
    fd = {}
    for k in d.keys():
        if f (d[k]) : 
            fd[k] = d[k]
    return fd
    

def ND (sig_digits_10) : 
    return lambda x : N (x, digits = sig_digits_10)
    
def chop (tolerance = None) :
    def f_ (x) :
        t = sqrt ((parent(x)(1)).ulp()) if tolerance == None else tolerance
        return x if (abs(x) >= t) else parent(x)(0)
    return f_
 
def NDC (sig_digits_10, tolerance = None) : 
    return lambda x : N (chop(tolerance) (x), digits = sig_digits_10)
   
def sqrt_ulp (field = RDF) :
    return sqrt (field(1).ulp())
    

def dict_joint_bars (free_nodes, elems) : 
    nds_els = {}
    for fni in free_nodes :
        nds_els[fni] = []
    for (ei, ec) in elems.iteritems() : 
        (i0, i1) = ec
        if i0 in free_nodes :
            nds_els[i0].append (ei)
        if i1 in free_nodes :
            nds_els[i1].append (ei)
    return nds_els


def dict_equation_index (free_nodes, dimen = 3) : 
    return dict (zip (sorted (free_nodes), xrange (0, dimen * len (free_nodes), dimen)))  
      
def dict_index_equation (free_nodes, dimen = 3) : 
    return dict (zip (xrange (0, dimen * len (free_nodes), dimen), sorted (free_nodes)))  
      
def dict_unknown_index (elems) : 
    return dict (zip (sorted (elems.keys()), xrange (len (elems))))   

def dict_index_unknown (elems) : 
    return dict (zip (xrange (len (elems)), sorted (elems.keys())))   

def dict_label_index (lbls) : 
    return dict (zip (sorted (lbls), xrange (len (lbls))))   

def dict_index_label (lbls) : 
    return dict (zip (xrange (len (lbls)), sorted (lbls)))   

        
def others (l1, l2) :
    return [c for c in l2 if c not in l1]    

def list_indexed (inds, l) :
    return [l[i] for i in inds]

def all_indices (l) :
    return xrange (len (l))
    

def maxwells_rule (joints, supports, bars) : 
    dim = len (joints.values()[0]) 
    free_joints = others (supports, joints)
    mxw = dim*len (free_joints) - len (bars)
    print dim, '*', len (free_joints), '-', len (bars), ' == ', mxw, 
    if mxw != 0 :
        print ' !=  0'
    else :
        print 
    return 


def base_vectors (dim, field = RDF) :
    bv = zero_vector (field, dim)
    bvs = [copy (bv) for i in xrange (dim)]
    for i in xrange (dim) :
        bvs[i][i] = field (1)
    return bvs


def bar_length (xxi, xxj) : 
    return sqrt (sum (map (lambda xi, xj : (xj - xi)^2, xxi, xxj)))


def equilibrium_matrix (nodes, elems, free_nodes, field = RDF) : 
    nfj = len (free_nodes)
    nb = len (elems)
    dim = len (nodes.values()[0]) 
    ne_inc = dict_joint_bars (free_nodes, elems)
    ukws = dict_unknown_index (elems)
    eqs = dict_equation_index (free_nodes, dim)
    A = zero_matrix (field, dim*nfj, nb)
    
    for ni in free_nodes : 
        ei = eqs[ni]
        bni = ne_inc[ni]
        for bi in bni :
            if ni == elems[bi][0] :
                nj = elems[bi][1]
            else :
                nj = elems[bi][0]
            xxi = map (field, nodes[ni]) 
            xxj = map (field, nodes[nj]) 
            # if field == SR :
            #     if ni < nj :
            #         l = var ('l_%d%d' % (ni, nj))
            #     else :
            #         l = var ('l_%d%d' % (nj, ni))
            # else :
            l = sqrt (sum (map (lambda xi, xj : (xj - xi)^2, xxi, xxj)))
            cc = map (lambda xi, xj : (xj - xi) / l, xxi, xxj)
            ui = ukws[bi]
            A[ei : (ei+dim), ui] += matrix (dim, 1, cc)
            
    return A


def load_vector (loads, free_nodes, field = RDF) :
    dim = len (loads.values()[0]) 
    nfj = len (free_nodes)
    eqs = dict_equation_index (free_nodes, dim)
    b = zero_vector (field, dim*nfj)
    for li, ff in loads.iteritems() :
        ei = eqs[li]
        b[ei : (ei+dim)] = map (field, ff)
    return b


def augmented_equil_matrix (nodes, elems, loads, free_nodes, field = RDF) : 
    nfj = len (free_nodes)
    nb = len (elems)
    dim = len (nodes.values()[0]) 
    ne_inc = dict_joint_bars (free_nodes, elems)
    ukws = dict_unknown_index (elems)
    eqs = dict_equation_index (free_nodes, dim)
    A = zero_matrix (field, dim*nfj, nb + 1)
    
    for ni in free_nodes : 
        ei = eqs[ni]
        bni = ne_inc[ni]
        for bi in bni :
            if ni == elems[bi][0] :
                nj = elems[bi][1]
            else :
                nj = elems[bi][0]
            xxi = map (field, nodes[ni]) 
            xxj = map (field, nodes[nj]) 
            l = sqrt (sum (map (lambda xi, xj : (xj - xi)^2, xxi, xxj)))
            cc = map (lambda xi, xj : (xj - xi) / l, xxi, xxj)
            ui = ukws[bi]
            A[ei : (ei+dim), ui] += matrix (dim, 1, cc)
    
    for li, ff in loads.iteritems() :
        ei = eqs[li]
        ej = nb
        A[ei : (ei+dim), ej] = matrix (dim, 1, map (lambda x : -field(x), ff))
        
    A.subdivide ([], [nb])
    return A


def dict_bar_force (elems, forces) :
    ukws = dict_label_index (elems.keys()) 
    return dict ([(b, forces[i]) for b, i in ukws.iteritems()]) 


def truss_skd (nodes, elems, supports, loads, field = RDF) : 
    free_nodes = others (supports, nodes)
    A = equilibrium_matrix (nodes, elems, free_nodes, field)
    b = load_vector (loads, free_nodes, field)
    ff = A.solve_right (-b)
    return dict_bar_force (elems, ff)
    

def is_zz (x) :
    return x == 0

def is_zz_tol (tolerance = None) :
    def f_ (x) :
        t = sqrt ((parent(x)(1)).ulp()) if tolerance == None else tolerance
        return abs (x) < t
    return f_

 
def ref_wpp (m, is_zero_f = is_zz_tol(), in_place = False) :
    pivot = 0
    if in_place == True : 
        r = m
    else :
        r = copy (m)
    pivotlist = []
    rnr = r.nrows()
    for i in xrange (r.ncols()) :
        if pivot == rnr :
            break
        maxrow = pivot
        for k in xrange (pivot+1, r.nrows()):
            if abs (r[k][i]) > abs (r[maxrow][i]) : 
                maxrow = k
        if pivot != maxrow :
            r.swap_rows (pivot, maxrow)
        if not is_zero_f (r[pivot, i]) :
            scale = r[pivot, i]
            r.rescale_row (pivot, 1/scale, i)
            for j in xrange (pivot+1, r.nrows()) :
                scale = r[j, i]
                r.add_multiple_of_row (j, pivot, -scale, i)
            pivotlist.append (i)
            pivot += 1
    return r, pivotlist

def ref_solution (Ar, piv) : 
    v = vector (Ar.base_ring(), Ar.ncols() - 1)
    np = len (piv)
    for j in xrange (np) :
        v[piv[j]] = Ar[j,-1]
    for j in xrange (np-1, -1, -1) :
        smj = 0 
        for jj in xrange (np-1, j, -1) :
            smj += Ar[j, piv[jj]] * v[piv[jj]]
        v[piv[j]] -= smj
    return v


def rref_wpp (m, is_zero_f = is_zz_tol(), in_place = False) :
    pivot = 0
    if in_place == True : 
        r = m
    else :
        r = copy (m)
    pivotlist = []
    for i in xrange (r.ncols()) :
        if pivot == r.nrows() :
            break
        maxrow = pivot
        for k in xrange (pivot+1, r.nrows()):
            if abs (r[k][i]) > abs (r[maxrow][i]) : 
                maxrow = k
        if pivot != maxrow :
            r.swap_rows (pivot, maxrow)
        if not is_zero_f (r[pivot, i]) :
            scale = r[pivot, i]
            r.rescale_row (pivot, 1/scale, i)
            for j in xrange (r.nrows()) :
                if j == pivot :
                    continue
                scale = r[j, i]
                r.add_multiple_of_row (j, pivot, -scale, i)
            pivotlist.append (i)
            pivot += 1
    return r, pivotlist

def rref_solution (Arr, piv) :
    v = vector (Arr.base_ring(), Arr.ncols() - 1)
    np = len (piv)
    for j in xrange (np) :
        v[piv[j]] = Arr[j,-1]
    return v
    
 
def truss (nodes, elems, supports, loads, is_zero_f = is_zz_tol(), field = RDF) : 
    free_nodes = others (supports, nodes)
    Ab = augmented_equil_matrix (nodes, elems, loads, free_nodes, field)
    Ab, piv = rref_wpp (Ab, is_zero_f = is_zero_f, in_place = True)
    return dict_bar_force (elems, rref_solution (Ab, piv))

    
def pivot_columns (mr) : 
    nr = mr.nrows()
    pc = []
    for i in xrange (mr.nrows()) : 
        for j in xrange (i, mr.ncols()) :
            if mr[i, j] != 0 :
                pc.append (j)
                break
    return pc


def ref_ext_wpp (m, b = None, is_zero_f = is_zz_tol()) :
    pivot = 0
    if b == None :
        r = block_matrix (1, 2, [m, identity_matrix (m.base_ring(), m.nrows())])
    else :
        r = block_matrix (1, 3, [m, identity_matrix (m.base_ring(), m.nrows()), b])
    mncols = m.ncols()
    pivotlist = []
    for i in xrange (mncols) :
        if pivot == r.nrows() :
            break
        maxrow = pivot
        for k in xrange (pivot+1, r.nrows()):
            if abs (r[k][i]) > abs (r[maxrow][i]) : 
                maxrow = k
        if pivot != maxrow :
            r.swap_rows (pivot, maxrow)
        if not is_zero_f (r[pivot, i]) :
            scale = r[pivot, i]
            r.rescale_row (pivot, 1/scale, i)
            for j in xrange (pivot+1, r.nrows()) :
                scale = r[j, i]
                r.add_multiple_of_row (j, pivot, -scale, i)
            pivotlist.append (i)
            pivot += 1
    return r, pivotlist


def rref_ext_wpp (m, b = None, is_zero_f = is_zz_tol()) :
    pivot = 0
    if b == None :
        r = block_matrix (1, 2, [m, identity_matrix (m.base_ring(), m.nrows())])
    else :
        r = block_matrix (1, 3, [m, identity_matrix (m.base_ring(), m.nrows()), b])
    mncols = m.ncols()
    pivotlist = []
    for i in xrange (mncols) :
        if pivot == r.nrows() :
            break
        maxrow = pivot
        for k in xrange (pivot+1, r.nrows()):
            if abs (r[k][i]) > abs (r[maxrow][i]) : 
                maxrow = k
        if pivot != maxrow :
            r.swap_rows (pivot, maxrow)
        if not is_zero_f (r[pivot, i]) :
            scale = r[pivot, i]
            r.rescale_row (pivot, 1/scale, i)
            for j in xrange (r.nrows()) :
                if j == pivot :
                    continue
                scale = r[j, i]
                r.add_multiple_of_row (j, pivot, -scale, i)
            pivotlist.append (i)
            pivot += 1
    return r, pivotlist


def column_space (m, piv, is_zero_f = is_zz_tol(), rref = False) :
    cs = []
    for p in piv :
        cs.append (m.column (p))
    if rref == False : 
        return cs
    mcs = matrix (cs)
    mcsr, pp = rref_wpp (mcs, is_zero_f = is_zero_f, in_place = True)
    return mcsr.rows()
    
def row_space (mr, piv) :
    rs = []
    m = mr 
    for i in xrange (len (piv)) :
        rs.append (m.row (i))
    return rs

def left_kernel (mre, piv) :
    lk = []
    ii = mre 
    nr = ii.nrows()
    r0 = len (piv)
    for i in xrange (r0, nr) :
        lk.append (ii.row (i))
    return lk

def subs_r_ (m, nc, piv, np, i) : 
    v = vector (m.base_ring(), nc)
    v[i] = 1
    for j in xrange (np) :
        v[piv[j]] = -m[j,i]
    return v
    
def backsubs_ (m, nc, piv, np, i) : 
    v = vector (m.base_ring(), nc)
    v[i] = 1
    for j in xrange (np-1, -1, -1) :
        smj = 0 
        for jj in xrange (np-1, j, -1) :
            smj += m[j, piv[jj]] * v[piv[jj]]
        v[piv[j]] = -m[j,i] - smj
    return v

def kernel (mr, piv, backsubs = backsubs_) :
    k = []
    m = mr 
    nc = m.ncols()
    np = len (piv)
    piv_s = set (piv)
    for i in xrange (nc) :
        if not i in piv_s : 
            k.append (backsubs (m, nc, piv, np, i))
    return k

def subspace_bases (m, is_zero_f = is_zz_tol()) :
    mre, piv = ref_ext_wpp (m, is_zero_f = is_zero_f)
    rs = row_space (mre.subdivision (0, 0), piv)
    cs = column_space (m, piv, is_zero_f = is_zero_f)
    rk = kernel (mre.subdivision (0, 0), piv)
    lk = left_kernel (mre.subdivision (0, 1), piv)
    return rs, cs, rk, lk

def subspace_bases_r (m, is_zero_f = is_zz_tol()) :
    mre, piv = rref_ext_wpp (m, is_zero_f = is_zero_f)
    rs = row_space (mre.subdivision (0, 0), piv)
    cs = column_space (m, piv, is_zero_f = is_zero_f)
    rk = kernel (mre.subdivision (0, 0), piv)
    lk = left_kernel (mre.subdivision (0, 1), piv)
    return rs, cs, rk, lk

def subspaces (m, is_zero_f = is_zz_tol(), return_pivots = False) :
    mre, piv = ref_ext_wpp (m, is_zero_f = is_zero_f)
    nj3 = m.nrows()
    nb = m.ncols()
    rm = len (piv)
    rsb = row_space (mre.subdivision (0, 0), piv)
    rsn = nb
    rsd = rm
    rkb = kernel (mre.subdivision (0, 0), piv)
    rkn = nb
    rkd = nb - rm
    csb = column_space (m, piv)
    csn = nj3
    csd = rm
    lkb = left_kernel (mre.subdivision (0, 1), piv)
    lkn = nj3
    lkd = nj3 - rm
    rs = {'name' : 'Row space', 'degree' : rsn, 'dimension' : rsd, 'basis' : rsb} 
    cs = {'name' : 'Column space', 'degree' : csn, 'dimension' : csd, 'basis' : csb} 
    rk = {'name' : 'Kernel', 'degree' : rkn, 'dimension' : rkd, 'basis' : rkb} 
    lk = {'name' : 'Left kernel', 'degree' : lkn, 'dimension' : lkd, 'basis' : lkb} 
    if return_pivots == True :
        r = rs, cs, rk, lk, piv
    else :
        r = rs, cs, rk, lk
    return r

def subspaces_r (m, is_zero_f = is_zz_tol(), return_pivots = False) :
    mre, piv = rref_ext_wpp (m, is_zero_f = is_zero_f)
    nj3 = m.nrows()
    nb = m.ncols()
    rm = len (piv)
    rsb = row_space (mre.subdivision (0, 0), piv)
    rsn = nb
    rsd = rm
    rkb = kernel (mre.subdivision (0, 0), piv)
    rkn = nb
    rkd = nb - rm
    csb = column_space (m, piv, rref = True, is_zero_f = is_zero_f)
    csn = nj3
    csd = rm
    lkb = left_kernel (mre.subdivision (0, 1), piv)
    lkn = nj3
    lkd = nj3 - rm
    rs = {'name' : 'Row space', 'degree' : rsn, 'dimension' : rsd, 'basis' : rsb} 
    cs = {'name' : 'Column space', 'degree' : csn, 'dimension' : csd, 'basis' : csb} 
    rk = {'name' : 'Kernel', 'degree' : rkn, 'dimension' : rkd, 'basis' : rkb} 
    lk = {'name' : 'Left kernel', 'degree' : lkn, 'dimension' : lkd, 'basis' : lkb} 
    if return_pivots == True :
        r = rs, cs, rk, lk, piv
    else :
        r = rs, cs, rk, lk
    return r

def print_subspace (ss, sb = True) :
    print '{0}{1}'.format (ss['name'], ':')
    print '   Vector space of degree', ss['degree'], 'and dimension', ss['dimension']
    if sb == True : 
        print '   Basis:'
        print_list (ss['basis'], '     ')
        
def are_linearly_dependent (vl, is_zero_f = is_zz_tol()) :
    m = matrix (vl) 
    mre, piv = rref_wpp (m, is_zero_f = is_zero_f)
    return len (vl) != len (piv)

def is_in_subspace (v, sp, is_zero_f = is_zz_tol()) : 
    vl = [v]
    vl.extend (sp)
    return are_linearly_dependent (vl, is_zero_f = is_zero_f)

def has_component_in_subspace (v, sp, is_zero_f = is_zz_tol()) :
    r = False
    for vv in sp :
        s = v * vv
        if not is_zero_f (s) :
            r = True
            break
    return r


def all_bar_forces (primary_bars, redundant_bars, primary_forces) : 
    bfs = vector (primary_forces.base_ring(), len (primary_bars) + len (redundant_bars))
    for i in xrange (len (primary_bars)) :
        bfs[primary_bars[i]] = primary_forces[i]
    return bfs    


def diag_flex (elems, nodes, E, F, field = RDF) :
    ldelta = []
    for el in elems :
        ni, nj = elems[el]
        xi = nodes[ni] 
        xj = nodes[nj] 
        ldelta.append (bar_length (xi, xj) / (E * F))
    return diagonal_matrix (RDF, ldelta)


def diag_stiff (elems, nodes, E, F, field = RDF) :
    lk = []
    for el in elems :
        ni, nj = elems[el]
        xi = nodes[ni] 
        xj = nodes[nj] 
        lk.append (E * F / bar_length (xi, xj))
    return diagonal_matrix (RDF, lk)


def make_bars__ (nb, b0j, steps) :
    return [map (lambda x, y : x + y, b0j, map (lambda x : i * x, steps)) for i in xrange (nb)]


def schwedler_1 (r0, hh, n) :
    if not is_even (n) :
        print 'n must be even'
        return

    jsl = []
    bsl = []
    ssl = []
    rsl = [r0]
    
    h0 = hh[-1]
    rr = r0^2/(2*h0) + h0/2
    aa = rr - h0
    for h in hh[:-1] :
        rsl.append (sqrt (rr^2 - (aa+h)^2))
    hh.insert (0, 0.0)
    alpha = RDF (2*pi/n)
    for i in xrange (len (rsl)) :
        for j in xrange (n) :
            x = rsl[i] * cos (j*alpha)
            y = rsl[i] * sin (j*alpha)
            z = hh[i]
            jsl.append ((x, y, z))
    jsd = dict (zip (xrange (len (jsl)), jsl))
    
    for i in xrange (0, len (hh) - 2) :
        n0 = i*n
        for j in xrange (n) :
            bsl.append ((n0+j, n0+j+n))
        for j in xrange (n) :
            bsl.append ((n0+n+j, n*(i+1)+(n0+j+1)%n))
        for j in xrange (0, n, 2) :
            bsl.append ((n0+j, (n*(i+1)+(n0+j+1)%n)))
        for j in xrange (0, n, 2) :
            bsl.append ((n0+j, (n*(i+1)+(n0+j-1)%n)))
    bsd = dict (zip (xrange (len (bsl)), bsl))
           
    ssl = range (n)
    return jsd, bsd, ssl


def schwedler_2 (r0, hh, n, ccw = True) :        
    jsl = []
    bsl = []
    ssl = []
    rsl = [r0]
    
    h0 = hh[-1]
    rr = r0^2/(2*h0) + h0/2
    aa = rr - h0
    for h in hh[:-1] :
        rsl.append (sqrt (rr^2 - (aa+h)^2))
    hh.insert (0, 0.0)
    alpha = RDF (2*pi/n) 
    for i in xrange (len (rsl)) :
        for j in xrange (n) : 
            x = rsl[i] * cos (j*alpha)
            y = rsl[i] * sin (j*alpha)
            z = hh[i]
            jsl.append ((x, y, z))
    jsd = dict (zip (xrange (len (jsl)), jsl))
    
    for i in xrange (0, len (hh) - 2) : 
        n0 = i*n 
        for j in xrange (n) : 
            bsl.append ((n0+j, n0+j+n)) 
        for j in xrange (n) : 
            bsl.append ((n0+n+j, n*(i+1)+(n0+j+1)%n))
        if ccw == True :
            for j in xrange (0, n) : 
                bsl.append ((n0+j, (n*(i+1)+(n0+j+1)%n)))
        else :
            for j in xrange (0, n) : 
                bsl.append ((n0+j, (n*(i+1)+(n0+j-1)%n)))
    bsd = dict (zip (xrange (len (bsl)), bsl))
           
    ssl = range (n)
    return jsd, bsd, ssl


def schwedler_kind (r0, hh, n) :        
    jsl = []
    bsl = []
    ssl = []
    rsl = [r0]
    
    h0 = hh[-1]
    rr = r0^2/(2*h0) + h0/2
    aa = rr - h0
    for h in hh[:-1] :
        rsl.append (sqrt (rr^2 - (aa+h)^2))
    hh.insert (0, 0.0)
    alpha = RDF (2*pi/n) 
    for i in xrange (len (rsl)) :
        for j in xrange (n) : 
            x = rsl[i] * cos (j*alpha)
            y = rsl[i] * sin (j*alpha)
            z = hh[i]
            jsl.append ((x, y, z))
    jsd = dict (zip (xrange (len (jsl)), jsl))
    
    for i in xrange (0, len (hh) - 2) : 
        n0 = i*n 
        for j in xrange (n) : 
            bsl.append ((n0+j, n0+j+n)) 
        for j in xrange (n) : 
            bsl.append ((n0+n+j, n*(i+1)+(n0+j+1)%n))
    bsd = dict (zip (xrange (len (bsl)), bsl))
           
    ssl = range (n)
    return jsd, bsd, ssl


def schwedler_sind (r0, hh, n) :        
    jsl = []
    bsl = []
    ssl = []
    rsl = [r0]
    
    h0 = hh[-1]
    rr = r0^2/(2*h0) + h0/2
    aa = rr - h0
    for h in hh[:-1] :
        rsl.append (sqrt (rr^2 - (aa+h)^2))
    hh.insert (0, 0.0)
    alpha = RDF (2*pi/n) 
    for i in xrange (len (rsl)) :
        for j in xrange (n) : 
            x = rsl[i] * cos (j*alpha)
            y = rsl[i] * sin (j*alpha)
            z = hh[i]
            jsl.append ((x, y, z))
    jsd = dict (zip (xrange (len (jsl)), jsl))
    
    for i in xrange (0, len (hh) - 2) : 
        n0 = i*n 
        for j in xrange (n) : 
            bsl.append ((n0+j, n0+j+n)) 
        for j in xrange (n) : 
            bsl.append ((n0+n+j, n*(i+1)+(n0+j+1)%n))
        for j in xrange (0, n) : 
            bsl.append ((n0+j, (n*(i+1)+(n0+j+1)%n)))
        for j in xrange (0, n) : 
            bsl.append ((n0+j, (n*(i+1)+(n0+j-1)%n)))
    bsd = dict (zip (xrange (len (bsl)), bsl))
           
    ssl = range (n)
    return jsd, bsd, ssl


def hypar_xy (x, y, z_scale) :
    return x * y * z_scale 
    
def hypar_x2y2 (x, y, z_scale) :
    return (x^2 - y^2) * z_scale 
    
def surface_net (f, lx, ly, m, n) : 
    jsl = []
    bsl = []
    ssl = []

    m2 = m + 2
    n2 = n + 2
    xmin = -lx/2
    ymin = -ly/2
    hx = lx/(m2-1)
    hy = ly/(n2-1) 
    
    for j in xrange (n2) : 
        for i in xrange (m2) : 
            xx = xmin + i*hx
            yy = ymin + j*hy
            jsl.append ((xx, yy, f (xx, yy)))
    jsd = dict (zip (xrange (len (jsl)), jsl))

    for j in xrange (1, n2-1) : 
        n0 = j*m2 
        for i in xrange (m2-1) : 
            bsl.append ((n0+i, n0+i+1)) 
    for i in xrange (1, m2-1) : 
        for j in xrange (0, (n2-1)*m2, m2) : 
            bsl.append ((i+j, i+j+m2)) 
    bsd = dict (zip (xrange (len (bsl)), bsl))
        
    ssl = range (m2)
    for i in xrange (1, n2-1) :
        ssl.append (i*m2)
        ssl.append (i*m2 + m2 - 1)
    ssl.extend (range ((n2-1)*m2, n2*m2))

    return jsd, bsd, ssl


def make_loads__ (nl, n0, step, ld) :
    return zip (range (n0, n0 + nl*step, step), [ld for i in xrange (nl)])
    
def make_loads (generator) :
    ll = []
    for gg in generator :
        if len (gg) == 3 : 
            nl, n0, ld = gg
            step = 1
        else :
            nl, n0, step, ld = gg
        ll.extend (make_loads__ (nl, n0, step, ld))
    return dict (ll)
    
    
def tension_compression (bar_forces, tolerance = 1.5e-8) :
    tt = [] 
    cc = []
    zz = []
    for (bi, f) in bar_forces.iteritems() : 
        if abs (f) > tolerance : 
            if f > 0 :
                tt.append ((bi, f))
            else :
                cc.append ((bi, f))
        else :
            zz.append (bi)
    tt = sorted (tt, key = lambda x : x[1])
    cc = sorted (cc, key = lambda x : x[1])
    return (tt, cc, zz)

def unique_lp (l, tolerance = 1e-8) : 
    li = []
    if l != []: 
        (ai, av) = l[0]
        ll = [ai]
        for (bi, bv) in l[1:] :
            if abs (av - bv) < tolerance :
                ll.append (bi)
            else : 
                li.append (copy (ll))
                ll = [bi] 
                av = bv
        li.append (copy (ll))
    return li   


def plot_free_joints2d (pts, rr, color = 'black') :
    return sum ([circle (pt, rr, edgecolor = color, facecolor = 'white',
                 fill = True, aspect_ratio = 1) for pt in pts])

def plot_supports2d (pts, rr) :
    return sum ([circle (pt, rr, edgecolor = 'black', facecolor = 'black',
                 fill = True, aspect_ratio = 1) for pt in pts])

def plot_bars2d (bars, joints, color = 'black', thickness = 2) :
    return sum ([line2d ((joints[b[0]], joints[b[1]]), thickness = thickness,
                 color = color, aspect_ratio = 1) for b in bars])

def plot_loads2d (loads, joints, load_scale) :
    pls = []
    for li in loads :
        tail = joints[li]
        ff = loads[li]
        head = [tail[0] + load_scale*ff[0], tail[1] + load_scale*ff[1]]
        pls.append (arrow2d (tail, head, width = 3, color = 'crimson', aspect_ratio = 1))
    return sum (pls)
                                                 
def plot_truss2d (joints, bars, supports, loads = [], load_scale = 0.01) :
    prr = point2d (joints.values(), aspect_ratio = 1)
    rrx = (prr.xmax() - prr.xmin())
    rry = (prr.ymax() - prr.ymin())
    rr = rrx if rrx > rry else rry*3/2
    pfjs = plot_free_joints2d (list_indexed (others (supports, joints), joints), rr/120)
    psjs = plot_supports2d (list_indexed (supports, joints), rr/150)
    pbs = plot_bars2d (bars.values(), joints)
    pt = pfjs + psjs + pbs
    if loads != [] :
        pls = plot_loads2d (loads, joints, load_scale)
        pt += pls
    return pt
                 
def plot_truss_with_forces2d (joints, bars, supports, forces) :
    prr = point2d (joints.values(), aspect_ratio = 1)
    rrx = (prr.xmax() - prr.xmin())
    rry = (prr.ymax() - prr.ymin())
    rr = rrx if rrx > rry else rry*3/2
    pfjs = plot_free_joints2d (list_indexed (others (supports, joints), joints), rr/120)
    psjs = plot_supports2d (list_indexed (supports, joints), rr/150)
    pt = pfjs + psjs
    tt, cc, zi = tension_compression (forces)
    ti = unique_lp (tt)
    ti.reverse()
    ci = unique_lp (cc)
    if zi != [] :
        cl = Color ('orange')
        cl = cl.lighter (0.425)
        pbz = plot_bars2d (list_indexed (zi, bars), joints, color = cl, thickness = 3)
        pt += pbz
    if ti != [] :
        step = 1./len (ti)
        cl = Color ('darkblue')
        pbt = plot_bars2d (list_indexed (ti[0], bars), joints, color = cl, thickness = 3)
        for i in xrange (1, len (ti)) :
            cl = cl.lighter (step)
            pbt += plot_bars2d (list_indexed (ti[i], bars), joints, color = cl, thickness = 3)
        pt += pbt
    if ci != [] :
        step = 1./len (ci)
        cl = Color ('red')
        pbt = plot_bars2d (list_indexed (ci[0], bars), joints, color = cl, thickness = 3)
        for i in xrange (1, len (ci)) :
            cl = cl.lighter (step)
            pbt += plot_bars2d (list_indexed (ci[i], bars), joints, color = cl, thickness = 3)
        pt += pbt
    return pt

def plot_truss_with_displacements2d (joints, bars, supports, displs, displ_scale = 1) : 
    prr = point2d (joints.values(), aspect_ratio = 1)
    rrx = (prr.xmax() - prr.xmin())
    rry = (prr.ymax() - prr.ymin())
    rr = rrx if rrx > rry else rry*3/2
    fjis = others (supports, joints)
    fjs = list_indexed (fjis, joints)
    sjs = list_indexed (supports, joints)
    pfjs = plot_free_joints2d (fjs, rr/180, 'grey')
    psjs = plot_supports2d (sjs, rr/150)
    pbs = plot_bars2d (bars.values(), joints, color = 'grey', thickness = 1)
    pt = pfjs + psjs + pbs
    ds = displ_scale
    dfjs = copy (fjs)
    dispc = map (chop(), displs)
    for i in xrange (len (dfjs)) :
        jnt = dfjs[i]
        dfjs[i] = (jnt[0] + ds*dispc[i*2], jnt[1] + ds*dispc[i*2 + 1])
    djs = copy (joints)
    for i in xrange (len (fjis)) :
        djs[fjis[i]] = dfjs[i] 
    pdfjs = plot_free_joints2d (dfjs, rr/120)
    pdbs = plot_bars2d (bars.values(), djs)
    pt = pbs + pfjs + pdbs + pdfjs + psjs
    return pt


def plot_free_joints3d (pts, color = 'white', size = 27) : 
    return point3d (pts, color = color, size = size, aspect_ratio = 1, frame = False)
  
def plot_supports3d (pts, color = 'black', size = 23) : 
    return point3d (pts, color = color, size = size, aspect_ratio = 1, frame = False)
  
def plot_bars3d (bars, joints, color = 'black', thickness = 7) : 
    return sum ([line3d ((joints[b[0]], joints[b[1]]), thickness = thickness,
                 color = color, aspect_ratio = 1, frame = False) for b in bars]) 


def plot_arrow3d (ff, tail = (0, 0, 0), scale = 1, color = 'blue', thickness = 1) :
    fs = scale * ff
    nf = norm (fs)
    if nf != 0 :
        fw = (3.425 * thickness * nf)^(2/5)
        # fw = sqrt (0.825 * radius * nf)
        # fw = radius * 3.25
        rd = fw / 50
        hl = 6.425 * rd
        hr = 1.625 * rd
        vtl = vector (tail)
        ar = arrow3d (vtl, vtl + fs, radius = rd, head_len = hl, head_radius = hr, color = color, aspect_ratio = 1, frame = False)
        # ar = arrow3d (vtl, vtl + fs, width = fw, head_len = hl, head_radius = hr, color = color)
    else :
        ar = Sphere (0.001, color = 'white').translate (tail)
    return ar

def plot_loads3d (loads, joints, load_scale, pull = True, thickness = 1, color = 'crimson') : 
    pls = []
    for li in loads :
        ff = vector (loads[li]) 
        if pull == True :
            tail = joints[li]
        else : 
            head = joints[li]
            tail = [head[0] - load_scale*ff[0], head[1] - load_scale*ff[1], head[2] - load_scale*ff[2]]
        pls.append (plot_arrow3d (ff, tail, color = color, thickness = thickness, scale = load_scale))
    return sum (pls)

    
def plot_truss3d (joints, bars, supports, loads = [], load_scale = 0.01, pull = True) : 
    pfjs = plot_free_joints3d (list_indexed (others (supports, joints), joints))
    psjs = plot_supports3d (list_indexed (supports, joints))
    pbs = plot_bars3d (bars.values(), joints)
    pt = pfjs + psjs + pbs
    if loads != [] :
        pls = plot_loads3d (loads, joints, load_scale, pull)
        pt += pls
    return pt

   
def plot_truss_with_forces3d (joints, bars, supports, forces) : 
    pfjs = plot_free_joints3d (list_indexed (others (supports, joints), joints))
    psjs = plot_supports3d (list_indexed (supports, joints))
    pt = pfjs + psjs 
    tt, cc, zi = tension_compression (forces)
    ti = unique_lp (tt)
    ti.reverse() 
    ci = unique_lp (cc)
    if zi != [] :
        cl = Color ('orange')
        cl = cl.lighter (0.625)
        pbz = plot_bars3d (list_indexed (zi, bars), joints, color = cl, thickness = 7) 
        pt += pbz
    if ti != [] :
        step = 1./len (ti) 
        cl = Color ('indigo')
        pbt = plot_bars3d (list_indexed (ti[0], bars), joints, color = cl, thickness = 7) 
        for i in xrange (1, len (ti)) :
            cl = cl.lighter (step)
            pbt += plot_bars3d (list_indexed (ti[i], bars), joints, color = cl, thickness = 7) 
        pt += pbt
    if ci != [] :
        step = 1./len (ci) 
        cl = Color ('red')
        pbt = plot_bars3d (list_indexed (ci[0], bars), joints, color = cl, thickness = 7) 
        for i in xrange (1, len (ci)) :
            cl = cl.lighter (step)
            pbt += plot_bars3d (list_indexed (ci[i], bars), joints, color = cl, thickness = 7) 
        pt += pbt
    return pt


def plot_disp_lines3d (joints, djoints, color = 'orange', thickness = 3) : 
    return sum ([line3d ((joints[i], djoints[i]), arrow_head = True, thickness = thickness,
                 color = color, aspect_ratio = 1, frame = False) for i in xrange (len (joints))]) 
                         
def plot_truss_with_displacements3d (joints, bars, supports, displs, 
                                     displ_scale = 1, plot_dlines = True) : 
    fjis = others (supports, joints)
    fjs = list_indexed (fjis, joints)
    pfjs = plot_free_joints3d (fjs, color = 'whitesmoke', size = 18)
    psjs = plot_supports3d (list_indexed (supports, joints))
    pbs = plot_bars3d (bars.values(), joints, color = 'whitesmoke', thickness = 4)
    ds = displ_scale
    dfjs = copy (fjs)
    dispc = map (chop(), displs)
    for i in xrange (len (dfjs)) :
        jnt = dfjs[i]
        dfjs[i] = (jnt[0] + ds*dispc[i*3], jnt[1] + ds*dispc[i*3 + 1], jnt[2] + ds*dispc[i*3 + 2])
    djs = copy (joints)
    for i in xrange (len (fjis)) :
        djs[fjis[i]] = dfjs[i] 
    pdfjs = plot_free_joints3d (dfjs)
    pdbs = plot_bars3d (bars.values(), djs)
    pt = pfjs + psjs + pbs + pdfjs + pdbs
    if plot_dlines == True :
        pdl = plot_disp_lines3d (fjs, dfjs)
        pt += pdl
    return pt


def plot_vectors3d (vs, color = 'green', thickness = 0.1) : 
    pvs = []
    coord_axes = []
    coord_axes.append (plot_arrow3d (vector ((1, 0, 0)), (0, 0, 0), thickness = 0.01))
    coord_axes.append (plot_arrow3d (vector ((0, 1, 0)), (0, 0, 0), thickness = 0.01))
    coord_axes.append (plot_arrow3d (vector ((0, 0, 1)), (0, 0, 0), thickness = 0.01))
    for v in vs :
        pvs.append (plot_arrow3d (v, (0, 0, 0), color = color, thickness = thickness))
    return sum (coord_axes) + sum (pvs)
    
    
def plot_vectors2d (vs, color = 'green', thickness = 1) :
    pvs = [] 
    for v in vs :
        pvs.append (v.plot (color = color, thickness = thickness, aspect_ratio = 1))
    return sum (pvs)


def print_indexed (l, i0 = 0, idx_len = 4, indent = '') :
    fstr = '{0}{1:<' + str (idx_len) + '} {2}'
    for i, e in enumerate (l) :
        print fstr.format (indent, str(i+i0) + ':', e)

def print_list (l, indent = '', braces = False) : 
    if l == [] :
        print '{0}{1}'.format (indent, '[]')
    else :
        if braces == True :
            print '{0}{1}'.format (indent, '[')
        for x in l :
            print '{0}{1}'.format (indent, x)
        if braces == True :
            print '{0}{1}'.format (indent, ']')
        
def print_dict (d, key_len = 4, sort_keys = True, indent = '') : 
    fstr = '{0}{1:<' + str (key_len) + '} {2}'
    for k in sorted (d.keys()) if sort_keys == True else d.keys():
        print fstr.format (indent, str(k) + ':', d[k])

def hypar2 (x, y) :
    return hypar_x2y2 (x, y, 0.175)

joints, bars, supports = surface_net (hypar2, 5., 6., 4, 5)

plot_truss3d (joints, bars, supports)

Click on 3-D image to make it live. Right-click on live image for a control menu.

maxwells_rule (joints, supports, bars)

3 * 20 - 49  ==  11  !=  0

3 * 20 - 49  ==  11  !=  0

free_joints = others (supports, joints)
free_joints

[7, 8, 9, 10, 13, 14, 15, 16, 19, 20, 21, 22, 25, 26, 27, 28, 31, 32, 33, 34]

[7, 8, 9, 10, 13, 14, 15, 16, 19, 20, 21, 22, 25, 26, 27, 28, 31, 32, 33, 34]

A = equilibrium_matrix (joints, bars, free_joints)
A

60 x 49 dense matrix over Real Double Field (use the '.str()' method to see the
entries)

60 x 49 dense matrix over Real Double Field (use the '.str()' method to see the entries)

rs, cs, rk, lk = subspaces (A)

print_subspace (rs, sb = False)

Row space:
   Vector space of degree 49 and dimension 48

Row space:
   Vector space of degree 49 and dimension 48

print_subspace (cs, sb = False)

Column space:
   Vector space of degree 60 and dimension 48

Column space:
   Vector space of degree 60 and dimension 48

print_subspace (rk, sb = False)

Kernel:
   Vector space of degree 49 and dimension 1

Kernel:
   Vector space of degree 49 and dimension 1

bar_forces_ss = dict_bar_force (bars, rk['basis'][0])
print_dict (bar_forces_ss)

0:   0.9186369279044543
1:   0.7973407128633153
2:   0.7525766947068775
3:   0.7973407128633153
4:   0.9186369279044544
5:   0.9186369279044532
6:   0.7973407128633143
7:   0.7525766947068765
8:   0.7973407128633143
9:   0.9186369279044532
10:  0.9186369279044533
11:  0.7973407128633144
12:  0.7525766947068766
13:  0.7973407128633144
14:  0.9186369279044533
15:  0.9186369279044534
16:  0.7973407128633145
17:  0.7525766947068767
18:  0.7973407128633145
19:  0.9186369279044534
20:  0.9186369279044537
21:  0.7973407128633148
22:  0.7525766947068769
23:  0.7973407128633148
24:  0.9186369279044537
25:  0.9999999999999988
26:  0.849986985845191
27:  0.7640136217105626
28:  0.7640136217105626
29:  0.8499869858451911
30:  0.9999999999999989
31:  0.9999999999999997
32:  0.8499869858451917
33:  0.7640136217105632
34:  0.7640136217105632
35:  0.8499869858451917
36:  0.9999999999999996
37:  0.9999999999999979
38:  0.8499869858451902
39:  0.7640136217105618
40:  0.7640136217105618
41:  0.8499869858451902
42:  0.9999999999999978
43:  0.9999999999999999
44:  0.849986985845192
45:  0.7640136217105635
46:  0.7640136217105635
47:  0.849986985845192
48:  1.0

0:   0.9186369279044543
1:   0.7973407128633153
2:   0.7525766947068775
3:   0.7973407128633153
4:   0.9186369279044544
5:   0.9186369279044532
6:   0.7973407128633143
7:   0.7525766947068765
8:   0.7973407128633143
9:   0.9186369279044532
10:  0.9186369279044533
11:  0.7973407128633144
12:  0.7525766947068766
13:  0.7973407128633144
14:  0.9186369279044533
15:  0.9186369279044534
16:  0.7973407128633145
17:  0.7525766947068767
18:  0.7973407128633145
19:  0.9186369279044534
20:  0.9186369279044537
21:  0.7973407128633148
22:  0.7525766947068769
23:  0.7973407128633148
24:  0.9186369279044537
25:  0.9999999999999988
26:  0.849986985845191
27:  0.7640136217105626
28:  0.7640136217105626
29:  0.8499869858451911
30:  0.9999999999999989
31:  0.9999999999999997
32:  0.8499869858451917
33:  0.7640136217105632
34:  0.7640136217105632
35:  0.8499869858451917
36:  0.9999999999999996
37:  0.9999999999999979
38:  0.8499869858451902
39:  0.7640136217105618
40:  0.7640136217105618
41:  0.8499869858451902
42:  0.9999999999999978
43:  0.9999999999999999
44:  0.849986985845192
45:  0.7640136217105635
46:  0.7640136217105635
47:  0.849986985845192
48:  1.0

plot_truss_with_forces3d (joints, bars, supports, bar_forces_ss)

print_subspace (lk, sb = False)

Left kernel:
   Vector space of degree 60 and dimension 12

Left kernel:
   Vector space of degree 60 and dimension 12

plot_truss_with_displacements3d (joints, bars, supports, map (chop(), lk['basis'][0]), displ_scale = 0.25)

plot_truss_with_displacements3d (joints, bars, supports, map (chop(), lk['basis'][1]), displ_scale = 0.25)

plot_truss_with_displacements3d (joints, bars, supports, map (chop(), lk['basis'][11]), displ_scale = 0.25)

plot_truss_with_displacements3d (joints, bars, supports, map (chop(), lk['basis'][10]), displ_scale = 0.25)

plot_truss_with_displacements3d (joints, bars, supports, map (chop(), lk['basis'][6]), displ_scale = 0.25)

def hypar1 (x, y) :
    return hypar_xy (x, y, 0.25)

joints, bars, supports = surface_net (hypar1, 5., 6., 4, 5)

plot_truss3d (joints, bars, supports)

maxwells_rule (joints, supports, bars)

3 * 20 - 49  ==  11  !=  0

3 * 20 - 49  ==  11  !=  0

free_joints = others (supports, joints)
len (free_joints)

A = equilibrium_matrix (joints, bars, free_joints)
A

60 x 49 dense matrix over Real Double Field (use the '.str()' method to see the
entries)

60 x 49 dense matrix over Real Double Field (use the '.str()' method to see the entries)

rs, cs, rk, lk = subspaces (A)

print_subspace (rs, sb = False)

Row space:
   Vector space of degree 49 and dimension 40

Row space:
   Vector space of degree 49 and dimension 40

print_subspace (cs, sb = False)

Column space:
   Vector space of degree 60 and dimension 40

Column space:
   Vector space of degree 60 and dimension 40

print_subspace (rk, sb = False)

Kernel:
   Vector space of degree 49 and dimension 9

Kernel:
   Vector space of degree 49 and dimension 9

bar_forces_ss0 = dict_bar_force (bars, rk['basis'][0])
print_dict (filter_dict (lambda x : not is_zz_tol()(x), bar_forces_ss0))

0:   1.0
1:   1.0
2:   1.0
3:   1.0
4:   1.0

0:   1.0
1:   1.0
2:   1.0
3:   1.0
4:   1.0

plot_truss_with_forces3d (joints, bars, supports, bar_forces_ss0)

bar_forces_ss1 = dict_bar_force (bars, rk['basis'][1])
print_dict (filter_dict (lambda x : not is_zz_tol()(x), bar_forces_ss1))

5:   1.0
6:   1.0
7:   1.0
8:   1.0
9:   1.0

5:   1.0
6:   1.0
7:   1.0
8:   1.0
9:   1.0

plot_truss_with_forces3d (joints, bars, supports, bar_forces_ss1)

bar_forces_ss8 = dict_bar_force (bars, rk['basis'][8])
print_dict (filter_dict (lambda x : not is_zz_tol()(x), bar_forces_ss8))

43:  0.9999999999999994
44:  0.9999999999999996
45:  0.9999999999999997
46:  0.9999999999999998
47:  0.9999999999999999
48:  1.0

43:  0.9999999999999994
44:  0.9999999999999996
45:  0.9999999999999997
46:  0.9999999999999998
47:  0.9999999999999999
48:  1.0

plot_truss_with_forces3d (joints, bars, supports, bar_forces_ss8)

print_subspace (lk, sb = False)

Left kernel:
   Vector space of degree 60 and dimension 20

Left kernel:
   Vector space of degree 60 and dimension 20

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][3], displ_scale = 0.325)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][15], displ_scale = 0.325)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][19], displ_scale = 0.325)

hh = 2
joints = { 0: (-sin (pi/3), -cos (pi/3), 0), 1: (sin (pi/3), -cos (pi/3), 0), 2: (0, 1, 0), 3: (-sin (pi/3), -cos (pi/3), hh), 4: (sin (pi/3), -cos (pi/3), hh), 5: (0, 1, hh) }
print 'joints:'
print_dict (joints, indent = '   ')
print

bars = { 0: (0, 1), 1: (1, 2), 2: (2, 0), 3: (3, 4), 4: (4, 5), 5: (5, 3), 6: (0, 3), 7: (1, 4), 8: (2, 5), 9: (0, 4), 10: (1, 5), 11: (2, 3) }
print 'bars:' 
print_dict (bars, indent = '   ')
print

supports = []
print 'supports:', supports
print

joints:
   0:   (-1/2*sqrt(3), -1/2, 0)
   1:   (1/2*sqrt(3), -1/2, 0)
   2:   (0, 1, 0)
   3:   (-1/2*sqrt(3), -1/2, 2)
   4:   (1/2*sqrt(3), -1/2, 2)
   5:   (0, 1, 2)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 0)
   3:   (3, 4)
   4:   (4, 5)
   5:   (5, 3)
   6:   (0, 3)
   7:   (1, 4)
   8:   (2, 5)
   9:   (0, 4)
   10:  (1, 5)
   11:  (2, 3)

supports: []

joints:
   0:   (-1/2*sqrt(3), -1/2, 0)
   1:   (1/2*sqrt(3), -1/2, 0)
   2:   (0, 1, 0)
   3:   (-1/2*sqrt(3), -1/2, 2)
   4:   (1/2*sqrt(3), -1/2, 2)
   5:   (0, 1, 2)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 0)
   3:   (3, 4)
   4:   (4, 5)
   5:   (5, 3)
   6:   (0, 3)
   7:   (1, 4)
   8:   (2, 5)
   9:   (0, 4)
   10:  (1, 5)
   11:  (2, 3)

supports: []

plot_truss3d (joints, bars, supports) . rotateY (-pi/69)

maxwells_rule (joints, supports, bars)

3 * 6 - 12  ==  6  !=  0

3 * 6 - 12  ==  6  !=  0

free_joints = others (supports, joints)
free_joints

[0, 1, 2, 3, 4, 5]

[0, 1, 2, 3, 4, 5]

A = equilibrium_matrix (joints, bars, free_joints) 
show (A.apply_map (NDC (4)))

$\newcommand{\Bold}[1]{\mathbf{#1}}\left(\begin{array}{rrrrrrrrrrrr} 1.000 & 0.0000 & 0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.6547 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 1.000 & 0.0000 & 0.0000 & 0.7559 & 0.0000 & 0.0000 \\ -1.000 & -0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.3273 & 0.0000 \\ 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.5669 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 1.000 & 0.0000 & 0.0000 & 0.7559 & 0.0000 \\ 0.0000 & 0.5000 & -0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.3273 \\ 0.0000 & -0.8660 & -0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.5669 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 1.000 & 0.0000 & 0.0000 & 0.7559 \\ 0.0000 & 0.0000 & 0.0000 & 1.000 & 0.0000 & 0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.3273 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.5669 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -1.000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.7559 \\ 0.0000 & 0.0000 & 0.0000 & -1.000 & -0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.6547 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -1.000 & 0.0000 & -0.7559 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.5000 & -0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.3273 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.8660 & -0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.5669 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -1.000 & 0.0000 & -0.7559 & 0.0000 \end{array}\right)$

rs, cs, rk, lk = subspaces (A)

print_subspace (rk)

Kernel:
   Vector space of degree 12 and dimension 0
   Basis:
     []

Kernel:
   Vector space of degree 12 and dimension 0
   Basis:
     []

print_subspace (lk)

Left kernel:
   Vector space of degree 18 and dimension 6
   Basis:
     (0.0, -0.5773502691896258, 0.43301270189221935, 0.0, 0.5773502691896258,
-0.43301270189221935, -1.0, 0.0, 0.0, 1.0, -0.5773502691896258,
0.43301270189221935, 1.0, 0.5773502691896258, -0.43301270189221935, 0.0, 0.0,
0.0)
     (1.0, 0.0, -0.43301270189221935, 1.0, 0.0, 0.43301270189221935, 1.0, 0.0,
0.0, 0.0, 0.0, -0.43301270189221935, 0.0, 0.0, 0.43301270189221935, 0.0, 0.0,
0.0)
     (0.0, 1.0, -0.75, 0.0, 1.0, -0.75, 0.0, 1.0, 0.0, 0.0, 0.0, -0.75, 0.0,
0.0, -0.75, 0.0, 0.0, 0.0)
     (0.0, 0.5773502691896258, 0.0, 0.0, -0.5773502691896258, 0.0, 1.0, 0.0,
0.0, 0.0, 0.5773502691896258, 0.0, 0.0, -0.5773502691896258, 0.0, 1.0, 0.0, 0.0)
     (0.0, 0.0, 0.75, 0.0, 0.0, 0.75, 0.0, 0.0, 0.0, 0.0, 1.0, 0.75, 0.0, 1.0,
0.75, 0.0, 1.0, 0.0)
     (0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0,
0.0, 0.0, 1.0)

Left kernel:
   Vector space of degree 18 and dimension 6
   Basis:
     (0.0, -0.5773502691896258, 0.43301270189221935, 0.0, 0.5773502691896258, -0.43301270189221935, -1.0, 0.0, 0.0, 1.0, -0.5773502691896258, 0.43301270189221935, 1.0, 0.5773502691896258, -0.43301270189221935, 0.0, 0.0, 0.0)
     (1.0, 0.0, -0.43301270189221935, 1.0, 0.0, 0.43301270189221935, 1.0, 0.0, 0.0, 0.0, 0.0, -0.43301270189221935, 0.0, 0.0, 0.43301270189221935, 0.0, 0.0, 0.0)
     (0.0, 1.0, -0.75, 0.0, 1.0, -0.75, 0.0, 1.0, 0.0, 0.0, 0.0, -0.75, 0.0, 0.0, -0.75, 0.0, 0.0, 0.0)
     (0.0, 0.5773502691896258, 0.0, 0.0, -0.5773502691896258, 0.0, 1.0, 0.0, 0.0, 0.0, 0.5773502691896258, 0.0, 0.0, -0.5773502691896258, 0.0, 1.0, 0.0, 0.0)
     (0.0, 0.0, 0.75, 0.0, 0.0, 0.75, 0.0, 0.0, 0.0, 0.0, 1.0, 0.75, 0.0, 1.0, 0.75, 0.0, 1.0, 0.0)
     (0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][0], displ_scale = 0.325)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][1], displ_scale = 0.325)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][4], displ_scale = 0.325)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][5], displ_scale = 0.325)

hh = 2
joints = { 0: (-sin (pi/3), -cos (pi/3), 0), 1: (sin (pi/3), -cos (pi/3), 0), 2: (0, 1, 0), 3: (0, -1, hh), 4: (sin (pi/3), cos (pi/3), hh), 5: (-sin (pi/3), cos (pi/3), hh) }
print 'joints:'
print_dict (joints, indent = '   ')
print

bars = { 0: (0, 1), 1: (1, 2), 2: (2, 0), 3: (3, 4), 4: (4, 5), 5: (5, 3), 6: (0, 3), 7: (1, 4), 8: (2, 5), 9: (0, 4), 10: (1, 5), 11: (2, 3) }
print 'bars:' 
print_dict (bars, indent = '   ')
print

supports = []
print 'supports:', supports
print

joints:
   0:   (-1/2*sqrt(3), -1/2, 0)
   1:   (1/2*sqrt(3), -1/2, 0)
   2:   (0, 1, 0)
   3:   (0, -1, 2)
   4:   (1/2*sqrt(3), 1/2, 2)
   5:   (-1/2*sqrt(3), 1/2, 2)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 0)
   3:   (3, 4)
   4:   (4, 5)
   5:   (5, 3)
   6:   (0, 3)
   7:   (1, 4)
   8:   (2, 5)
   9:   (0, 4)
   10:  (1, 5)
   11:  (2, 3)

supports: []

joints:
   0:   (-1/2*sqrt(3), -1/2, 0)
   1:   (1/2*sqrt(3), -1/2, 0)
   2:   (0, 1, 0)
   3:   (0, -1, 2)
   4:   (1/2*sqrt(3), 1/2, 2)
   5:   (-1/2*sqrt(3), 1/2, 2)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 0)
   3:   (3, 4)
   4:   (4, 5)
   5:   (5, 3)
   6:   (0, 3)
   7:   (1, 4)
   8:   (2, 5)
   9:   (0, 4)
   10:  (1, 5)
   11:  (2, 3)

supports: []

plot_truss3d (joints, bars, supports)

maxwells_rule (joints, supports, bars)

3 * 6 - 12  ==  6  !=  0

3 * 6 - 12  ==  6  !=  0

free_joints = others (supports, joints)
free_joints

[0, 1, 2, 3, 4, 5]

[0, 1, 2, 3, 4, 5]

A = equilibrium_matrix (joints, bars, free_joints) 
show (A.apply_map (NDC (4)))

$\newcommand{\Bold}[1]{\mathbf{#1}}\left(\begin{array}{rrrrrrrrrrrr} 1.000 & 0.0000 & 0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.3873 & 0.0000 & 0.0000 & 0.6124 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & -0.2236 & 0.0000 & 0.0000 & 0.3536 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.8944 & 0.0000 & 0.0000 & 0.7071 & 0.0000 & 0.0000 \\ -1.000 & -0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.6124 & 0.0000 \\ 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.4472 & 0.0000 & 0.0000 & 0.3536 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.8944 & 0.0000 & 0.0000 & 0.7071 & 0.0000 \\ 0.0000 & 0.5000 & -0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.3873 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & -0.8660 & -0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.2236 & 0.0000 & 0.0000 & -0.7071 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.8944 & 0.0000 & 0.0000 & 0.7071 \\ 0.0000 & 0.0000 & 0.0000 & 0.5000 & 0.0000 & -0.5000 & -0.3873 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.8660 & 0.2236 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.7071 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.8944 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.7071 \\ 0.0000 & 0.0000 & 0.0000 & -0.5000 & -1.000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.6124 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & -0.8660 & 0.0000 & 0.0000 & 0.0000 & -0.4472 & 0.0000 & -0.3536 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.8944 & 0.0000 & -0.7071 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 1.000 & 0.5000 & 0.0000 & 0.0000 & 0.3873 & 0.0000 & 0.6124 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.8660 & 0.0000 & 0.0000 & 0.2236 & 0.0000 & -0.3536 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.8944 & 0.0000 & -0.7071 & 0.0000 \end{array}\right)$

rs, cs, rk, lk = subspaces (A)

print_subspace (rk)

Kernel:
   Vector space of degree 12 and dimension 0
   Basis:
     []

Kernel:
   Vector space of degree 12 and dimension 0
   Basis:
     []

print_subspace (lk)

Left kernel:
   Vector space of degree 18 and dimension 6
   Basis:
     (0.0, 1.0, -0.5, 0.0, 1.0, -0.5, 0.0, 1.0, 0.25, 0.0, 0.0, -0.75, 0.0, 0.0,
0.0, 0.0, 0.0, 0.0)
     (0.0, -0.49999999999999994, 0.24999999999999992, 0.0, 0.49999999999999994,
-0.25000000000000006, -0.8660254037844385, 0.0, -0.37499999999999994,
0.8660254037844385, 0.5, 0.12499999999999994, 0.0, 1.0, -0.5, 0.0, 0.0, 0.0)
     (1.0, 0.0, 1.1102230246251565e-16, 1.0, 0.0, 0.8660254037844387, 1.0, 0.0,
0.4330127018922193, 0.0, 0.0, 0.4330127018922195, 0.0, 0.0, 0.8660254037844386,
0.0, 0.0, 0.0)
     (0.0, 0.0, 0.0, 0.0, 0.0, -0.8660254037844386, 0.0, 0.0,
-0.43301270189221935, 1.0, 0.0, -0.43301270189221935, 1.0, 0.0,
-0.8660254037844386, 1.0, 0.0, 0.0)
     (0.0, 0.49999999999999994, 0.2500000000000001, 0.0, -0.49999999999999994,
0.75, 0.8660254037844385, 0.0, 0.12499999999999994, -0.8660254037844385, 0.5,
0.625, 0.0, 0.0, 0.5, 0.0, 1.0, 0.0)
     (0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0,
0.0, 0.0, 1.0)

Left kernel:
   Vector space of degree 18 and dimension 6
   Basis:
     (0.0, 1.0, -0.5, 0.0, 1.0, -0.5, 0.0, 1.0, 0.25, 0.0, 0.0, -0.75, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
     (0.0, -0.49999999999999994, 0.24999999999999992, 0.0, 0.49999999999999994, -0.25000000000000006, -0.8660254037844385, 0.0, -0.37499999999999994, 0.8660254037844385, 0.5, 0.12499999999999994, 0.0, 1.0, -0.5, 0.0, 0.0, 0.0)
     (1.0, 0.0, 1.1102230246251565e-16, 1.0, 0.0, 0.8660254037844387, 1.0, 0.0, 0.4330127018922193, 0.0, 0.0, 0.4330127018922195, 0.0, 0.0, 0.8660254037844386, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.0, 0.0, 0.0, -0.8660254037844386, 0.0, 0.0, -0.43301270189221935, 1.0, 0.0, -0.43301270189221935, 1.0, 0.0, -0.8660254037844386, 1.0, 0.0, 0.0)
     (0.0, 0.49999999999999994, 0.2500000000000001, 0.0, -0.49999999999999994, 0.75, 0.8660254037844385, 0.0, 0.12499999999999994, -0.8660254037844385, 0.5, 0.625, 0.0, 0.0, 0.5, 0.0, 1.0, 0.0)
     (0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][0], displ_scale = 0.325)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][5], displ_scale = 0.325)

hh = 2
joints = { 0: (-sin (pi/3), -cos (pi/3), 0), 1: (sin (pi/3), -cos (pi/3), 0), 2: (0, 1, 0), 3: (-1/2, -1/2*sqrt (3), hh), 4: (1, 0, hh), 5: (-1/2, 1/2*sqrt (3), hh) }
print 'joints:'
print_dict (joints, indent = '   ')
print

bars = { 0: (0, 1), 1: (1, 2), 2: (2, 0), 3: (3, 4), 4: (4, 5), 5: (5, 3), 6: (0, 3), 7: (1, 4), 8: (2, 5), 9: (0, 4), 10: (1, 5), 11: (2, 3) }
print 'bars:' 
print_dict (bars, indent = '   ')
print

supports = []
print 'supports:', supports
print

joints:
   0:   (-1/2*sqrt(3), -1/2, 0)
   1:   (1/2*sqrt(3), -1/2, 0)
   2:   (0, 1, 0)
   3:   (-1/2, -1/2*sqrt(3), 2)
   4:   (1, 0, 2)
   5:   (-1/2, 1/2*sqrt(3), 2)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 0)
   3:   (3, 4)
   4:   (4, 5)
   5:   (5, 3)
   6:   (0, 3)
   7:   (1, 4)
   8:   (2, 5)
   9:   (0, 4)
   10:  (1, 5)
   11:  (2, 3)

supports: []

joints:
   0:   (-1/2*sqrt(3), -1/2, 0)
   1:   (1/2*sqrt(3), -1/2, 0)
   2:   (0, 1, 0)
   3:   (-1/2, -1/2*sqrt(3), 2)
   4:   (1, 0, 2)
   5:   (-1/2, 1/2*sqrt(3), 2)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 0)
   3:   (3, 4)
   4:   (4, 5)
   5:   (5, 3)
   6:   (0, 3)
   7:   (1, 4)
   8:   (2, 5)
   9:   (0, 4)
   10:  (1, 5)
   11:  (2, 3)

supports: []

plot_truss3d (joints, bars, supports)

maxwells_rule (joints, supports, bars)

3 * 6 - 12  ==  6  !=  0

3 * 6 - 12  ==  6  !=  0

free_joints = others (supports, joints)
free_joints

[0, 1, 2, 3, 4, 5]

[0, 1, 2, 3, 4, 5]

A = equilibrium_matrix (joints, bars, free_joints) 
show (A.apply_map (NDC (4)))

$\newcommand{\Bold}[1]{\mathbf{#1}}\left(\begin{array}{rrrrrrrrrrrr} 1.000 & 0.0000 & 0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.1772 & 0.0000 & 0.0000 & 0.6711 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & -0.1772 & 0.0000 & 0.0000 & 0.1798 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.9681 & 0.0000 & 0.0000 & 0.7193 & 0.0000 & 0.0000 \\ -1.000 & -0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.06485 & 0.0000 & 0.0000 & -0.4913 & 0.0000 \\ 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.2420 & 0.0000 & 0.0000 & 0.4913 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.9681 & 0.0000 & 0.0000 & 0.7193 & 0.0000 \\ 0.0000 & 0.5000 & -0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.2420 & 0.0000 & 0.0000 & -0.1798 \\ 0.0000 & -0.8660 & -0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.06485 & 0.0000 & 0.0000 & -0.6711 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.9681 & 0.0000 & 0.0000 & 0.7193 \\ 0.0000 & 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & -0.1772 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.1798 \\ 0.0000 & 0.0000 & 0.0000 & 0.5000 & 0.0000 & 1.000 & 0.1772 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.6711 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.9681 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.7193 \\ 0.0000 & 0.0000 & 0.0000 & -0.8660 & -0.8660 & 0.0000 & 0.0000 & -0.06485 & 0.0000 & -0.6711 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & -0.5000 & 0.5000 & 0.0000 & 0.0000 & -0.2420 & 0.0000 & -0.1798 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.9681 & 0.0000 & -0.7193 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.2420 & 0.0000 & 0.4913 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.5000 & -1.000 & 0.0000 & 0.0000 & 0.06485 & 0.0000 & -0.4913 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.9681 & 0.0000 & -0.7193 & 0.0000 \end{array}\right)$

rs, cs, rk, lk = subspaces (A)

print_subspace (rs)

Row space:
   Vector space of degree 12 and dimension 11
   Basis:
     (1.0, 0.0, 0.5000000000000001, 0.0, 0.0, 0.0, 0.17717466000345752, 0.0,
0.0, 0.6710736745934498, 0.0, 0.0)
     (0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.2794664976500341, 0.0, 0.0,
0.5672582158417004, 0.0)
     (0.0, 0.0, 1.0, 0.0, 0.0, 0.0, -0.20458367529315324, 0.0, 0.0,
0.20763091750349863, 0.0, 0.0)
     (0.0, 0.0, 0.0, 1.0, 0.0, 0.0, -0.20458367529315324, 0.0, 0.0, 0.0, 0.0,
0.20763091750349863)
     (0.0, 0.0, 0.0, 0.0, 1.0, -0.0, 0.20458367529315324, 0.07488282235688086,
-0.0, 0.774889133345199, -0.0, -0.20763091750349863)
     (0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.2794664976500342, 0.0, 0.0, 0.0, 0.0,
0.5672582158417004)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, -0.0, -0.0, -0.0, -0.0,
0.7429545911564234)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.7429545911564233, 0.0)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.7429545911564233)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0,
-0.9999999999999999)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0,
-0.9999999999999999)

Row space:
   Vector space of degree 12 and dimension 11
   Basis:
     (1.0, 0.0, 0.5000000000000001, 0.0, 0.0, 0.0, 0.17717466000345752, 0.0, 0.0, 0.6710736745934498, 0.0, 0.0)
     (0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.2794664976500341, 0.0, 0.0, 0.5672582158417004, 0.0)
     (0.0, 0.0, 1.0, 0.0, 0.0, 0.0, -0.20458367529315324, 0.0, 0.0, 0.20763091750349863, 0.0, 0.0)
     (0.0, 0.0, 0.0, 1.0, 0.0, 0.0, -0.20458367529315324, 0.0, 0.0, 0.0, 0.0, 0.20763091750349863)
     (0.0, 0.0, 0.0, 0.0, 1.0, -0.0, 0.20458367529315324, 0.07488282235688086, -0.0, 0.774889133345199, -0.0, -0.20763091750349863)
     (0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.2794664976500342, 0.0, 0.0, 0.0, 0.0, 0.5672582158417004)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, -0.0, -0.0, -0.0, -0.0, 0.7429545911564234)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.7429545911564233, 0.0)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.7429545911564233)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, -0.9999999999999999)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, -0.9999999999999999)

print_subspace (cs)

Column space:
   Vector space of degree 18 and dimension 11
   Basis:
     (1.0, 0.0, 0.0, -1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.0, -0.5000000000000001, 0.8660254037844387, 0.0,
0.5000000000000001, -0.8660254037844387, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
0.0, 0.0)
     (0.5000000000000001, 0.8660254037844387, 0.0, 0.0, 0.0, 0.0,
-0.5000000000000001, -0.8660254037844387, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
0.0, 0.0, 0.0)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.8660254037844387,
0.5000000000000001, 0.0, -0.8660254037844387, -0.5000000000000001, 0.0, 0.0,
0.0, 0.0)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
-0.8660254037844387, 0.5000000000000001, 0.0, 0.8660254037844387,
-0.5000000000000001, 0.0)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0,
0.0, -1.0, 0.0)
     (0.17717466000345752, -0.17717466000345752, 0.9681003458863746, 0.0, 0.0,
0.0, 0.0, 0.0, 0.0, -0.17717466000345752, 0.17717466000345752,
-0.9681003458863746, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.0, 0.06485042646813614, 0.2420250864715937,
0.9681003458863748, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, -0.06485042646813614,
-0.2420250864715937, -0.9681003458863748, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, -0.2420250864715937, -0.06485042646813614,
0.9681003458863748, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.2420250864715937,
0.06485042646813614, -0.9681003458863748)
     (0.6710736745934498, 0.1798136491691009, 0.7192545966764036, 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, -0.6710736745934498, -0.1798136491691009,
-0.7192545966764036, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.0, -0.49126002542434893, 0.49126002542434893,
0.7192545966764036, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
0.49126002542434893, -0.49126002542434893, -0.7192545966764036)

Column space:
   Vector space of degree 18 and dimension 11
   Basis:
     (1.0, 0.0, 0.0, -1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.0, -0.5000000000000001, 0.8660254037844387, 0.0, 0.5000000000000001, -0.8660254037844387, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
     (0.5000000000000001, 0.8660254037844387, 0.0, 0.0, 0.0, 0.0, -0.5000000000000001, -0.8660254037844387, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.8660254037844387, 0.5000000000000001, 0.0, -0.8660254037844387, -0.5000000000000001, 0.0, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, -0.8660254037844387, 0.5000000000000001, 0.0, 0.8660254037844387, -0.5000000000000001, 0.0)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, -1.0, 0.0)
     (0.17717466000345752, -0.17717466000345752, 0.9681003458863746, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, -0.17717466000345752, 0.17717466000345752, -0.9681003458863746, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.0, 0.06485042646813614, 0.2420250864715937, 0.9681003458863748, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, -0.06485042646813614, -0.2420250864715937, -0.9681003458863748, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, -0.2420250864715937, -0.06485042646813614, 0.9681003458863748, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.2420250864715937, 0.06485042646813614, -0.9681003458863748)
     (0.6710736745934498, 0.1798136491691009, 0.7192545966764036, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, -0.6710736745934498, -0.1798136491691009, -0.7192545966764036, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.0, -0.49126002542434893, 0.49126002542434893, 0.7192545966764036, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.49126002542434893, -0.49126002542434893, -0.7192545966764036)

print_subspace (rk)

Kernel:
   Vector space of degree 12 and dimension 1
   Basis:
     (-0.3596272983382016, -0.35962729833820173, -0.35962729833820173,
-0.3596272983382018, -0.35962729833820173, -0.35962729833820173,
-0.7429545911564234, -0.7429545911564231, -0.7429545911564233,
0.9999999999999999, 0.9999999999999999, 1.0)

Kernel:
   Vector space of degree 12 and dimension 1
   Basis:
     (-0.3596272983382016, -0.35962729833820173, -0.35962729833820173, -0.3596272983382018, -0.35962729833820173, -0.35962729833820173, -0.7429545911564234, -0.7429545911564231, -0.7429545911564233, 0.9999999999999999, 0.9999999999999999, 1.0)

bar_forces_ss0 = dict_bar_force (bars, rk['basis'][0])
print_dict (bar_forces_ss0)

0:   -0.3596272983382016
1:   -0.35962729833820173
2:   -0.35962729833820173
3:   -0.3596272983382018
4:   -0.35962729833820173
5:   -0.35962729833820173
6:   -0.7429545911564234
7:   -0.7429545911564231
8:   -0.7429545911564233
9:   0.9999999999999999
10:  0.9999999999999999
11:  1.0

0:   -0.3596272983382016
1:   -0.35962729833820173
2:   -0.35962729833820173
3:   -0.3596272983382018
4:   -0.35962729833820173
5:   -0.35962729833820173
6:   -0.7429545911564234
7:   -0.7429545911564231
8:   -0.7429545911564233
9:   0.9999999999999999
10:  0.9999999999999999
11:  1.0

bar_forces_ss0 = dict_bar_force (bars, -rk['basis'][0])
print_dict (bar_forces_ss0)

0:   0.3596272983382016
1:   0.35962729833820173
2:   0.35962729833820173
3:   0.3596272983382018
4:   0.35962729833820173
5:   0.35962729833820173
6:   0.7429545911564234
7:   0.7429545911564231
8:   0.7429545911564233
9:   -0.9999999999999999
10:  -0.9999999999999999
11:  -1.0

0:   0.3596272983382016
1:   0.35962729833820173
2:   0.35962729833820173
3:   0.3596272983382018
4:   0.35962729833820173
5:   0.35962729833820173
6:   0.7429545911564234
7:   0.7429545911564231
8:   0.7429545911564233
9:   -0.9999999999999999
10:  -0.9999999999999999
11:  -1.0

plot_truss_with_forces3d (joints, bars, supports, bar_forces_ss0) .rotateX (-pi/69)

print_subspace (lk)

Left kernel:
   Vector space of degree 18 and dimension 7
   Basis:
     (1.0, -0.5773502691896258, -0.28867513459481287, 1.0, 0.5773502691896258,
0.28867513459481287, 0.0, 0.0, 0.0, 0.0, 0.0, 5.551115123125783e-17, 0.0, 0.0,
0.5, 0.0, 0.0, 0.0)
     (0.0, 0.5773502691896258, 0.10566243270259365, 0.0, -0.5773502691896258,
0.39433756729740654, 1.0, 0.0, 0.25, 0.0, 0.0, 5.551115123125783e-17, 0.0, 0.0,
0.2500000000000001, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.5000000000000002, 0.0, 0.0, -5.551115123125783e-17, 0.0, 0.0,
0.0, 1.1547005383792517, 0.0, 0.28867513459481303, 0.5773502691896258, 1.0,
-0.2886751345948129, 0.0, 0.0, 0.0)
     (0.0, 1.0, -0.6830127018922194, 0.0, 1.0, -0.6830127018922194, 0.0, 1.0,
0.06698729810778065, 0.0, 0.0, -0.8660254037844388, 0.0, 0.0,
-0.43301270189221946, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.1830127018922193, 0.0, 0.0, -0.6830127018922194, 0.0, 0.0,
-0.25, 1.0, 0.0, -1.1102230246251565e-16, 1.0, 0.0, -0.7500000000000001, 1.0,
0.0, 0.0)
     (0.0, 0.0, 0.1830127018922193, 0.0, 0.0, 0.6830127018922194, 0.0, 0.0,
-0.06698729810778065, -1.1547005383792517, 1.0, 0.577350269189626,
-0.5773502691896258, 0.0, 0.7216878364870323, 0.0, 1.0, 0.0)
     (0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0,
0.0, 0.0, 1.0)

Left kernel:
   Vector space of degree 18 and dimension 7
   Basis:
     (1.0, -0.5773502691896258, -0.28867513459481287, 1.0, 0.5773502691896258, 0.28867513459481287, 0.0, 0.0, 0.0, 0.0, 0.0, 5.551115123125783e-17, 0.0, 0.0, 0.5, 0.0, 0.0, 0.0)
     (0.0, 0.5773502691896258, 0.10566243270259365, 0.0, -0.5773502691896258, 0.39433756729740654, 1.0, 0.0, 0.25, 0.0, 0.0, 5.551115123125783e-17, 0.0, 0.0, 0.2500000000000001, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.5000000000000002, 0.0, 0.0, -5.551115123125783e-17, 0.0, 0.0, 0.0, 1.1547005383792517, 0.0, 0.28867513459481303, 0.5773502691896258, 1.0, -0.2886751345948129, 0.0, 0.0, 0.0)
     (0.0, 1.0, -0.6830127018922194, 0.0, 1.0, -0.6830127018922194, 0.0, 1.0, 0.06698729810778065, 0.0, 0.0, -0.8660254037844388, 0.0, 0.0, -0.43301270189221946, 0.0, 0.0, 0.0)
     (0.0, 0.0, 0.1830127018922193, 0.0, 0.0, -0.6830127018922194, 0.0, 0.0, -0.25, 1.0, 0.0, -1.1102230246251565e-16, 1.0, 0.0, -0.7500000000000001, 1.0, 0.0, 0.0)
     (0.0, 0.0, 0.1830127018922193, 0.0, 0.0, 0.6830127018922194, 0.0, 0.0, -0.06698729810778065, -1.1547005383792517, 1.0, 0.577350269189626, -0.5773502691896258, 0.0, 0.7216878364870323, 0.0, 1.0, 0.0)
     (0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][0], displ_scale = 0.325)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][1], displ_scale = 0.325)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][5], displ_scale = 0.325)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][6], displ_scale = 0.325)

hh = 2
joints = { 0: (-sin (pi/3), -cos (pi/3), 0), 1: (sin (pi/3), -cos (pi/3), 0), 2: (0, 1, 0), 3: (-1/2, -1/2*sqrt(3), hh), 4: (1, 0, hh), 5: (-1/2, 1/2*sqrt(3), hh), 6: (-sin (pi/3) - 0.3, -cos (pi/3), 0), 7: (-sin (pi/3), -cos (pi/3) - 0.3, 0), 8: (-sin (pi/3), -cos (pi/3), - 0.3), 9: (sin (pi/3), -cos (pi/3) - 0.3, 0), 10: (sin (pi/3), -cos (pi/3), - 0.3), 11: (0, 1, -0.3) }
print 'joints:'
print_dict (joints, indent = '   ')
print

bars = { 0: (0, 1), 1: (1, 2), 2: (2, 0), 3: (3, 4), 4: (4, 5), 5: (5, 3), 6: (0, 3), 7: (1, 4), 8: (2, 5), 9: (0, 4), 10: (1, 5), 11: (2, 3), 12: (0, 6), 13: (0, 7), 14: (0, 8), 15: (1, 9), 16: (1, 10), 17: (2, 11) }
print 'bars:' 
print_dict (bars, indent = '   ')
print

supports = [6, 7, 8, 9, 10, 11]
print 'supports:', supports
print

joints:
   0:   (-1/2*sqrt(3), -1/2, 0)
   1:   (1/2*sqrt(3), -1/2, 0)
   2:   (0, 1, 0)
   3:   (-1/2, -1/2*sqrt(3), 2)
   4:   (1, 0, 2)
   5:   (-1/2, 1/2*sqrt(3), 2)
   6:   (-1/2*sqrt(3) - 0.300000000000000, -1/2, 0)
   7:   (-1/2*sqrt(3), -0.800000000000000, 0)
   8:   (-1/2*sqrt(3), -1/2, -0.300000000000000)
   9:   (1/2*sqrt(3), -0.800000000000000, 0)
   10:  (1/2*sqrt(3), -1/2, -0.300000000000000)
   11:  (0, 1, -0.300000000000000)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 0)
   3:   (3, 4)
   4:   (4, 5)
   5:   (5, 3)
   6:   (0, 3)
   7:   (1, 4)
   8:   (2, 5)
   9:   (0, 4)
   10:  (1, 5)
   11:  (2, 3)
   12:  (0, 6)
   13:  (0, 7)
   14:  (0, 8)
   15:  (1, 9)
   16:  (1, 10)
   17:  (2, 11)

supports: [6, 7, 8, 9, 10, 11]

joints:
   0:   (-1/2*sqrt(3), -1/2, 0)
   1:   (1/2*sqrt(3), -1/2, 0)
   2:   (0, 1, 0)
   3:   (-1/2, -1/2*sqrt(3), 2)
   4:   (1, 0, 2)
   5:   (-1/2, 1/2*sqrt(3), 2)
   6:   (-1/2*sqrt(3) - 0.300000000000000, -1/2, 0)
   7:   (-1/2*sqrt(3), -0.800000000000000, 0)
   8:   (-1/2*sqrt(3), -1/2, -0.300000000000000)
   9:   (1/2*sqrt(3), -0.800000000000000, 0)
   10:  (1/2*sqrt(3), -1/2, -0.300000000000000)
   11:  (0, 1, -0.300000000000000)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 0)
   3:   (3, 4)
   4:   (4, 5)
   5:   (5, 3)
   6:   (0, 3)
   7:   (1, 4)
   8:   (2, 5)
   9:   (0, 4)
   10:  (1, 5)
   11:  (2, 3)
   12:  (0, 6)
   13:  (0, 7)
   14:  (0, 8)
   15:  (1, 9)
   16:  (1, 10)
   17:  (2, 11)

supports: [6, 7, 8, 9, 10, 11]

plot_truss3d (joints, bars, supports)

maxwells_rule (joints, supports, bars)

3 * 6 - 18  ==  0

3 * 6 - 18  ==  0

free_joints = others (supports, joints)
free_joints

[0, 1, 2, 3, 4, 5]

[0, 1, 2, 3, 4, 5]

A = equilibrium_matrix (joints, bars, free_joints) 
show (A.apply_map (NDC (4)))

$\newcommand{\Bold}[1]{\mathbf{#1}}\left(\begin{array}{rrrrrrrrrrrrrrrrrr} 1.000 & 0.0000 & 0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.1772 & 0.0000 & 0.0000 & 0.6711 & 0.0000 & 0.0000 & -1.000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & -0.1772 & 0.0000 & 0.0000 & 0.1798 & 0.0000 & 0.0000 & 0.0000 & -1.000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.9681 & 0.0000 & 0.0000 & 0.7193 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -1.000 & 0.0000 & 0.0000 & 0.0000 \\ -1.000 & -0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.06485 & 0.0000 & 0.0000 & -0.4913 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.2420 & 0.0000 & 0.0000 & 0.4913 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -1.000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.9681 & 0.0000 & 0.0000 & 0.7193 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -1.000 & 0.0000 \\ 0.0000 & 0.5000 & -0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.2420 & 0.0000 & 0.0000 & -0.1798 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & -0.8660 & -0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.06485 & 0.0000 & 0.0000 & -0.6711 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.9681 & 0.0000 & 0.0000 & 0.7193 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -1.000 \\ 0.0000 & 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & -0.1772 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.1798 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.5000 & 0.0000 & 1.000 & 0.1772 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.6711 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.9681 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.7193 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & -0.8660 & -0.8660 & 0.0000 & 0.0000 & -0.06485 & 0.0000 & -0.6711 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & -0.5000 & 0.5000 & 0.0000 & 0.0000 & -0.2420 & 0.0000 & -0.1798 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.9681 & 0.0000 & -0.7193 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.2420 & 0.0000 & 0.4913 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.5000 & -1.000 & 0.0000 & 0.0000 & 0.06485 & 0.0000 & -0.4913 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.9681 & 0.0000 & -0.7193 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \end{array}\right)$

rs, cs, rk, lk = subspaces (A)

print_subspace (rs, sb = False)

Row space:
   Vector space of degree 18 and dimension 17

Row space:
   Vector space of degree 18 and dimension 17

print_subspace (cs, sb = False)

Column space:
   Vector space of degree 18 and dimension 17

Column space:
   Vector space of degree 18 and dimension 17

print_subspace (rk)

Kernel:
   Vector space of degree 18 and dimension 1
   Basis:
     (-0.3596272983382016, -0.35962729833820173, -0.35962729833820173,
-0.3596272983382018, -0.3596272983382016, -0.35962729833820173,
-0.7429545911564234, -0.7429545911564231, -0.7429545911564233,
0.9999999999999998, 0.9999999999999998, 1.0, 9.92914159000831e-17,
3.519265968729763e-17, -1.1102230246251565e-16, -0.0, -0.0, -0.0)

Kernel:
   Vector space of degree 18 and dimension 1
   Basis:
     (-0.3596272983382016, -0.35962729833820173, -0.35962729833820173, -0.3596272983382018, -0.3596272983382016, -0.35962729833820173, -0.7429545911564234, -0.7429545911564231, -0.7429545911564233, 0.9999999999999998, 0.9999999999999998, 1.0, 9.92914159000831e-17, 3.519265968729763e-17, -1.1102230246251565e-16, -0.0, -0.0, -0.0)

bar_forces_ss0 = dict_bar_force (bars, -rk['basis'][0])
print_dict (map_dict (chop(), bar_forces_ss0))

0:   0.3596272983382016
1:   0.35962729833820173
2:   0.35962729833820173
3:   0.3596272983382018
4:   0.3596272983382016
5:   0.35962729833820173
6:   0.7429545911564234
7:   0.7429545911564231
8:   0.7429545911564233
9:   -0.9999999999999998
10:  -0.9999999999999998
11:  -1.0
12:  0.0
13:  0.0
14:  0.0
15:  0.0
16:  0.0
17:  0.0

0:   0.3596272983382016
1:   0.35962729833820173
2:   0.35962729833820173
3:   0.3596272983382018
4:   0.3596272983382016
5:   0.35962729833820173
6:   0.7429545911564234
7:   0.7429545911564231
8:   0.7429545911564233
9:   -0.9999999999999998
10:  -0.9999999999999998
11:  -1.0
12:  0.0
13:  0.0
14:  0.0
15:  0.0
16:  0.0
17:  0.0

plot_truss_with_forces3d (joints, bars, supports, bar_forces_ss0)

print_subspace (lk)

Left kernel:
   Vector space of degree 18 and dimension 1
   Basis:
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.8660254037844388,
-0.5000000000000004, -0.25000000000000017, -1.1102230246251565e-16, 1.0, -0.25,
-0.866025403784439, -0.5000000000000004, -0.2500000000000001)

Left kernel:
   Vector space of degree 18 and dimension 1
   Basis:
     (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.8660254037844388, -0.5000000000000004, -0.25000000000000017, -1.1102230246251565e-16, 1.0, -0.25, -0.866025403784439, -0.5000000000000004, -0.2500000000000001)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][0], displ_scale = 0.425)

hh = 1.25
joints = { 0: (cos (pi/3), sin (pi/3), 0), 1: (0, 0, 0), 2: (1, 0, 0), 3: (0, 0, hh), 4: (cos (pi/6), sin (pi/6), hh), 5: (0, 1, hh) }
print 'joints:'
print_dict (joints, indent = '   ')
print

bars = { 0: (0, 1), 1: (1, 2), 2: (2, 0), 3: (3, 4), 4: (4, 5), 5: (5, 3), 6: (0, 5), 7: (1, 3), 8: (2, 4), 9: (0, 3), 10: (1, 4), 11: (2, 5) }
print 'bars:' 
print_dict (bars, indent = '   ')
print

supports = []
print 'supports:', supports
print

joints:
   0:   (1/2, 1/2*sqrt(3), 0)
   1:   (0, 0, 0)
   2:   (1, 0, 0)
   3:   (0, 0, 1.25000000000000)
   4:   (1/2*sqrt(3), 1/2, 1.25000000000000)
   5:   (0, 1, 1.25000000000000)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 0)
   3:   (3, 4)
   4:   (4, 5)
   5:   (5, 3)
   6:   (0, 5)
   7:   (1, 3)
   8:   (2, 4)
   9:   (0, 3)
   10:  (1, 4)
   11:  (2, 5)

supports: []

joints:
   0:   (1/2, 1/2*sqrt(3), 0)
   1:   (0, 0, 0)
   2:   (1, 0, 0)
   3:   (0, 0, 1.25000000000000)
   4:   (1/2*sqrt(3), 1/2, 1.25000000000000)
   5:   (0, 1, 1.25000000000000)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 0)
   3:   (3, 4)
   4:   (4, 5)
   5:   (5, 3)
   6:   (0, 5)
   7:   (1, 3)
   8:   (2, 4)
   9:   (0, 3)
   10:  (1, 4)
   11:  (2, 5)

supports: []

plot_truss3d (joints, bars, supports)

maxwells_rule (joints, supports, bars)

3 * 6 - 12  ==  6  !=  0

3 * 6 - 12  ==  6  !=  0

free_joints = others (supports, joints)
free_joints

[0, 1, 2, 3, 4, 5]

[0, 1, 2, 3, 4, 5]

A = equilibrium_matrix (joints, bars, free_joints) 
show (A.apply_map (NDC (4)))

$\newcommand{\Bold}[1]{\mathbf{#1}}\left(\begin{array}{rrrrrrrrrrrr} -0.5000 & 0.0000 & 0.5000 & 0.0000 & 0.0000 & 0.0000 & -0.3696 & 0.0000 & 0.0000 & -0.3123 & 0.0000 & 0.0000 \\ -0.8660 & 0.0000 & -0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.09902 & 0.0000 & 0.0000 & -0.5410 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.9239 & 0.0000 & 0.0000 & 0.7809 & 0.0000 & 0.0000 \\ 0.5000 & 1.000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.5410 & 0.0000 \\ 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.3123 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 1.000 & 0.0000 & 0.0000 & 0.7809 & 0.0000 \\ 0.0000 & -1.000 & -0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.09902 & 0.0000 & 0.0000 & -0.5298 \\ 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.3696 & 0.0000 & 0.0000 & 0.5298 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.9239 & 0.0000 & 0.0000 & 0.6623 \\ 0.0000 & 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.3123 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.5000 & 0.0000 & 1.000 & 0.0000 & 0.0000 & 0.0000 & 0.5410 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -1.000 & 0.0000 & -0.7809 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & -0.8660 & -0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.09902 & 0.0000 & -0.5410 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & -0.5000 & 0.5000 & 0.0000 & 0.0000 & 0.0000 & -0.3696 & 0.0000 & -0.3123 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.9239 & 0.0000 & -0.7809 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.3696 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.5298 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.5000 & -1.000 & -0.09902 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.5298 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.9239 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.6623 \end{array}\right)$

rs, cs, rk, lk = subspaces (A)

print_subspace (rs, sb = False)

Row space:
   Vector space of degree 12 and dimension 11

Row space:
   Vector space of degree 12 and dimension 11

print_subspace (cs, sb = False)

Column space:
   Vector space of degree 18 and dimension 11

Column space:
   Vector space of degree 18 and dimension 11

print_subspace (rk)

Kernel:
   Vector space of degree 12 and dimension 1
   Basis:
     (-0.3058876451607489, -0.30588764516074907, -0.305887645160749,
-0.30588764516074896, -0.30588764516074896, -0.305887645160749,
-0.716805621929536, -0.6622661785325219, -0.716805621929536, 0.8481145238787241,
0.8481145238787241, 1.0)

Kernel:
   Vector space of degree 12 and dimension 1
   Basis:
     (-0.3058876451607489, -0.30588764516074907, -0.305887645160749, -0.30588764516074896, -0.30588764516074896, -0.305887645160749, -0.716805621929536, -0.6622661785325219, -0.716805621929536, 0.8481145238787241, 0.8481145238787241, 1.0)

bar_forces_ss0 = dict_bar_force (bars, -rk['basis'][0])
print_dict (bar_forces_ss0)

0:   0.3058876451607489
1:   0.30588764516074907
2:   0.305887645160749
3:   0.30588764516074896
4:   0.30588764516074896
5:   0.305887645160749
6:   0.716805621929536
7:   0.6622661785325219
8:   0.716805621929536
9:   -0.8481145238787241
10:  -0.8481145238787241
11:  -1.0

0:   0.3058876451607489
1:   0.30588764516074907
2:   0.305887645160749
3:   0.30588764516074896
4:   0.30588764516074896
5:   0.305887645160749
6:   0.716805621929536
7:   0.6622661785325219
8:   0.716805621929536
9:   -0.8481145238787241
10:  -0.8481145238787241
11:  -1.0

plot_truss_with_forces3d (joints, bars, supports, bar_forces_ss0)

print_subspace (lk, sb = False)

Left kernel:
   Vector space of degree 18 and dimension 7

Left kernel:
   Vector space of degree 18 and dimension 7

hh = 1.25
joints = { 0: (cos (pi/3), sin (pi/3), 0), 1: (0, 0, 0), 2: (1, 0, 0), 3: (0, 0, hh), 4: (cos (pi/6), sin (pi/6), hh), 5: (0, 1, hh), 6: (cos (pi/3), sin (pi/3), -0.3), 7: (-0.3, 0, 0), 8: (0, -0.3, 0), 9: (0, 0, -0.3), 10: (1, -0.3, 0), 11: (1, 0, -0.3) }
print 'joints:'
print_dict (joints, indent = '   ')
print

bars = { 0: (0, 1), 1: (1, 2), 2: (2, 0), 3: (3, 4), 4: (4, 5), 5: (5, 3), 6: (0, 5), 7: (1, 3), 8: (2, 4), 9: (0, 3), 10: (1, 4), 11: (2, 5), 12: (0, 6), 13:(1, 7), 14: (1, 8), 15: (1, 9), 16: (2, 10), 17: (2, 11) }
print 'bars:' 
print_dict (bars, indent = '   ')
print

supports = [6, 7, 8, 9, 10, 11]
print 'supports:', supports
print

joints:
   0:   (1/2, 1/2*sqrt(3), 0)
   1:   (0, 0, 0)
   2:   (1, 0, 0)
   3:   (0, 0, 1.25000000000000)
   4:   (1/2*sqrt(3), 1/2, 1.25000000000000)
   5:   (0, 1, 1.25000000000000)
   6:   (1/2, 1/2*sqrt(3), -0.300000000000000)
   7:   (-0.300000000000000, 0, 0)
   8:   (0, -0.300000000000000, 0)
   9:   (0, 0, -0.300000000000000)
   10:  (1, -0.300000000000000, 0)
   11:  (1, 0, -0.300000000000000)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 0)
   3:   (3, 4)
   4:   (4, 5)
   5:   (5, 3)
   6:   (0, 5)
   7:   (1, 3)
   8:   (2, 4)
   9:   (0, 3)
   10:  (1, 4)
   11:  (2, 5)
   12:  (0, 6)
   13:  (1, 7)
   14:  (1, 8)
   15:  (1, 9)
   16:  (2, 10)
   17:  (2, 11)

supports: [6, 7, 8, 9, 10, 11]

joints:
   0:   (1/2, 1/2*sqrt(3), 0)
   1:   (0, 0, 0)
   2:   (1, 0, 0)
   3:   (0, 0, 1.25000000000000)
   4:   (1/2*sqrt(3), 1/2, 1.25000000000000)
   5:   (0, 1, 1.25000000000000)
   6:   (1/2, 1/2*sqrt(3), -0.300000000000000)
   7:   (-0.300000000000000, 0, 0)
   8:   (0, -0.300000000000000, 0)
   9:   (0, 0, -0.300000000000000)
   10:  (1, -0.300000000000000, 0)
   11:  (1, 0, -0.300000000000000)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 0)
   3:   (3, 4)
   4:   (4, 5)
   5:   (5, 3)
   6:   (0, 5)
   7:   (1, 3)
   8:   (2, 4)
   9:   (0, 3)
   10:  (1, 4)
   11:  (2, 5)
   12:  (0, 6)
   13:  (1, 7)
   14:  (1, 8)
   15:  (1, 9)
   16:  (2, 10)
   17:  (2, 11)

supports: [6, 7, 8, 9, 10, 11]

plot_truss3d (joints, bars, supports)

maxwells_rule (joints, supports, bars)

3 * 6 - 18  ==  0

3 * 6 - 18  ==  0

free_joints = others (supports, joints)
free_joints

[0, 1, 2, 3, 4, 5]

[0, 1, 2, 3, 4, 5]

A = equilibrium_matrix (joints, bars, free_joints) 
show (A.apply_map (NDC (4)))

$\newcommand{\Bold}[1]{\mathbf{#1}}\left(\begin{array}{rrrrrrrrrrrrrrrrrr} -0.5000 & 0.0000 & 0.5000 & 0.0000 & 0.0000 & 0.0000 & -0.3696 & 0.0000 & 0.0000 & -0.3123 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ -0.8660 & 0.0000 & -0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.09902 & 0.0000 & 0.0000 & -0.5410 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.9239 & 0.0000 & 0.0000 & 0.7809 & 0.0000 & 0.0000 & -1.000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.5000 & 1.000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.5410 & 0.0000 & 0.0000 & -1.000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.3123 & 0.0000 & 0.0000 & 0.0000 & -1.000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 1.000 & 0.0000 & 0.0000 & 0.7809 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -1.000 & 0.0000 & 0.0000 \\ 0.0000 & -1.000 & -0.5000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.09902 & 0.0000 & 0.0000 & -0.5298 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.3696 & 0.0000 & 0.0000 & 0.5298 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -1.000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.9239 & 0.0000 & 0.0000 & 0.6623 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -1.000 \\ 0.0000 & 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.3123 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.5000 & 0.0000 & 1.000 & 0.0000 & 0.0000 & 0.0000 & 0.5410 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -1.000 & 0.0000 & -0.7809 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & -0.8660 & -0.8660 & 0.0000 & 0.0000 & 0.0000 & 0.09902 & 0.0000 & -0.5410 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & -0.5000 & 0.5000 & 0.0000 & 0.0000 & 0.0000 & -0.3696 & 0.0000 & -0.3123 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.9239 & 0.0000 & -0.7809 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.8660 & 0.0000 & 0.3696 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.5298 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.5000 & -1.000 & -0.09902 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.5298 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \\ 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.9239 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & -0.6623 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 & 0.0000 \end{array}\right)$

rs, cs, rk, lk = subspaces (A)

print_subspace (rs, sb = False)

Row space:
   Vector space of degree 18 and dimension 17

Row space:
   Vector space of degree 18 and dimension 17

print_subspace (cs, sb = False)

Column space:
   Vector space of degree 18 and dimension 17

Column space:
   Vector space of degree 18 and dimension 17

print_subspace (rk)

Kernel:
   Vector space of degree 18 and dimension 1
   Basis:
     (-0.3058876451607489, -0.30588764516074907, -0.305887645160749,
-0.30588764516074896, -0.30588764516074896, -0.305887645160749,
-0.716805621929536, -0.6622661785325219, -0.7168056219295359,
0.8481145238787241, 0.848114523878724, 1.0, 1.232595164407831e-32,
1.1102230246251563e-17, -2.2204460492503132e-17, -5.551115123125784e-17,
2.2204460492503132e-17, 5.551115123125783e-17)

Kernel:
   Vector space of degree 18 and dimension 1
   Basis:
     (-0.3058876451607489, -0.30588764516074907, -0.305887645160749, -0.30588764516074896, -0.30588764516074896, -0.305887645160749, -0.716805621929536, -0.6622661785325219, -0.7168056219295359, 0.8481145238787241, 0.848114523878724, 1.0, 1.232595164407831e-32, 1.1102230246251563e-17, -2.2204460492503132e-17, -5.551115123125784e-17, 2.2204460492503132e-17, 5.551115123125783e-17)

bar_forces_ss0 = dict_bar_force (bars, -rk['basis'][0])
print_dict (filter_dict (lambda x : not is_zz_tol()(x), bar_forces_ss0))

0:   0.3058876451607489
1:   0.30588764516074907
2:   0.305887645160749
3:   0.30588764516074896
4:   0.30588764516074896
5:   0.305887645160749
6:   0.716805621929536
7:   0.6622661785325219
8:   0.7168056219295359
9:   -0.8481145238787241
10:  -0.848114523878724
11:  -1.0

0:   0.3058876451607489
1:   0.30588764516074907
2:   0.305887645160749
3:   0.30588764516074896
4:   0.30588764516074896
5:   0.305887645160749
6:   0.716805621929536
7:   0.6622661785325219
8:   0.7168056219295359
9:   -0.8481145238787241
10:  -0.848114523878724
11:  -1.0

plot_truss_with_forces3d (joints, bars, supports, bar_forces_ss0)

print_subspace (lk)

Left kernel:
   Vector space of degree 18 and dimension 1
   Basis:
     (0.0, 0.0, 0.0, 0.0, 0.0, 7.153584028445419e-17, 0.0, 0.0, 0.0,
-0.9999999999999997, 0.5773502691896257, 7.153584028445419e-17,
1.1102230246251565e-16, -1.1547005383792515, 0.4618802153517006, 1.0,
0.5773502691896257, 0.33811978464829945)

Left kernel:
   Vector space of degree 18 and dimension 1
   Basis:
     (0.0, 0.0, 0.0, 0.0, 0.0, 7.153584028445419e-17, 0.0, 0.0, 0.0, -0.9999999999999997, 0.5773502691896257, 7.153584028445419e-17, 1.1102230246251565e-16, -1.1547005383792515, 0.4618802153517006, 1.0, 0.5773502691896257, 0.33811978464829945)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][0], displ_scale = 0.25)

hh = 3
ll = 0.3
joints = { 0: (-1, -1, 0), 1: (1, -1, 0), 2: (1, 1, 0), 3: (-1, 1, 0), 4: (0, -sqrt (2), hh), 5: (sqrt (2), 0, hh), 6: (0, sqrt (2), hh), 7: (-sqrt (2), 0, hh) }
print 'joints:'
print_dict (joints, indent = '   ')
print

bars = { 0: (0, 1), 1: (1, 2), 2: (2, 3), 3: (3, 0), 4: (4, 5), 5: (5, 6), 6: (6, 7), 7: (7, 4), 8: (0, 4), 9: (1, 5), 10: (2, 6), 11: (3, 7), 12: (0, 5), 13: (1, 6), 14: (2, 7), 15: (3, 4) }
print 'bars:' 
print_dict (bars, indent = '   ')
print

supports = []
print 'supports:', supports
print

joints:
   0:   (-1, -1, 0)
   1:   (1, -1, 0)
   2:   (1, 1, 0)
   3:   (-1, 1, 0)
   4:   (0, -sqrt(2), 3)
   5:   (sqrt(2), 0, 3)
   6:   (0, sqrt(2), 3)
   7:   (-sqrt(2), 0, 3)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 3)
   3:   (3, 0)
   4:   (4, 5)
   5:   (5, 6)
   6:   (6, 7)
   7:   (7, 4)
   8:   (0, 4)
   9:   (1, 5)
   10:  (2, 6)
   11:  (3, 7)
   12:  (0, 5)
   13:  (1, 6)
   14:  (2, 7)
   15:  (3, 4)

supports: []

joints:
   0:   (-1, -1, 0)
   1:   (1, -1, 0)
   2:   (1, 1, 0)
   3:   (-1, 1, 0)
   4:   (0, -sqrt(2), 3)
   5:   (sqrt(2), 0, 3)
   6:   (0, sqrt(2), 3)
   7:   (-sqrt(2), 0, 3)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 3)
   3:   (3, 0)
   4:   (4, 5)
   5:   (5, 6)
   6:   (6, 7)
   7:   (7, 4)
   8:   (0, 4)
   9:   (1, 5)
   10:  (2, 6)
   11:  (3, 7)
   12:  (0, 5)
   13:  (1, 6)
   14:  (2, 7)
   15:  (3, 4)

supports: []

plot_truss3d (joints, bars, supports) . rotateZ (pi/3)

maxwells_rule (joints, supports, bars)

3 * 8 - 16  ==  8  !=  0

3 * 8 - 16  ==  8  !=  0

free_joints = others (supports, joints)
free_joints

[0, 1, 2, 3, 4, 5, 6, 7]

[0, 1, 2, 3, 4, 5, 6, 7]

A = equilibrium_matrix (joints, bars, free_joints) 
A

24 x 16 dense matrix over Real Double Field (use the '.str()' method to see the
entries)

24 x 16 dense matrix over Real Double Field (use the '.str()' method to see the entries)

rs, cs, rk, lk = subspaces (A)

print_subspace (rs, sb = False)

Row space:
   Vector space of degree 16 and dimension 15

Row space:
   Vector space of degree 16 and dimension 15

print_subspace (cs, sb = False)

Column space:
   Vector space of degree 24 and dimension 15

Column space:
   Vector space of degree 24 and dimension 15

print_subspace (rk)

Kernel:
   Vector space of degree 16 and dimension 1
   Basis:
     (-0.3554644040543818, -0.35546440405438173, -0.3554644040543817,
-0.3554644040543817, -0.3554644040543818, -0.3554644040543817,
-0.3554644040543817, -0.35546440405438184, -0.8016322430806896,
-0.8016322430806894, -0.8016322430806895, -0.8016322430806896, 1.0,
0.9999999999999998, 0.9999999999999999, 1.0)

Kernel:
   Vector space of degree 16 and dimension 1
   Basis:
     (-0.3554644040543818, -0.35546440405438173, -0.3554644040543817, -0.3554644040543817, -0.3554644040543818, -0.3554644040543817, -0.3554644040543817, -0.35546440405438184, -0.8016322430806896, -0.8016322430806894, -0.8016322430806895, -0.8016322430806896, 1.0, 0.9999999999999998, 0.9999999999999999, 1.0)

bar_forces_ss0 = dict_bar_force (bars, -rk['basis'][0])
print_dict (bar_forces_ss0)

0:   0.3554644040543818
1:   0.35546440405438173
2:   0.3554644040543817
3:   0.3554644040543817
4:   0.3554644040543818
5:   0.3554644040543817
6:   0.3554644040543817
7:   0.35546440405438184
8:   0.8016322430806896
9:   0.8016322430806894
10:  0.8016322430806895
11:  0.8016322430806896
12:  -1.0
13:  -0.9999999999999998
14:  -0.9999999999999999
15:  -1.0

0:   0.3554644040543818
1:   0.35546440405438173
2:   0.3554644040543817
3:   0.3554644040543817
4:   0.3554644040543818
5:   0.3554644040543817
6:   0.3554644040543817
7:   0.35546440405438184
8:   0.8016322430806896
9:   0.8016322430806894
10:  0.8016322430806895
11:  0.8016322430806896
12:  -1.0
13:  -0.9999999999999998
14:  -0.9999999999999999
15:  -1.0

plot_truss_with_forces3d (joints, bars, supports, bar_forces_ss0) . rotateZ (pi/3)

print_subspace (lk, sb = False)

Left kernel:
   Vector space of degree 24 and dimension 9

Left kernel:
   Vector space of degree 24 and dimension 9

hh = 3
ll = 0.3
joints = { 0: (-1, -1, 0), 1: (1, -1, 0), 2: (1, 1, 0), 3: (-1, 1, 0), 4: (0, -sqrt (2), hh), 5: (sqrt (2), 0, hh), 6: (0, sqrt (2), hh), 7: (-sqrt (2), 0, hh), 8: (-1 - ll, -1, 0), 9: (-1, -1 - ll, 0), 10: (-1, -1, -ll), 11: (1, -1 - ll, 0),  12: (1, -1, -ll), 13: (sqrt (2), -ll, hh) }
print 'joints:'
print_dict (joints, indent = '   ')
print

bars = { 0: (0, 1), 1: (1, 2), 2: (2, 3), 3: (3, 0), 4: (4, 5), 5: (5, 6), 6: (6, 7), 7: (7, 4), 8: (0, 4), 9: (1, 5), 10: (2, 6), 11: (3, 7), 12: (0, 5), 13: (1, 6), 14: (2, 7), 15: (3, 4), 16: (0, 8), 17: (0, 9), 18: (0, 10), 19: (1, 11), 20: (1, 12,), 21: (5, 13) }
print 'bars:' 
print_dict (bars, indent = '   ')
print

supports = [8, 9, 10, 11, 12, 13]
print 'supports:', supports
print

joints:
   0:   (-1, -1, 0)
   1:   (1, -1, 0)
   2:   (1, 1, 0)
   3:   (-1, 1, 0)
   4:   (0, -sqrt(2), 3)
   5:   (sqrt(2), 0, 3)
   6:   (0, sqrt(2), 3)
   7:   (-sqrt(2), 0, 3)
   8:   (-1.30000000000000, -1, 0)
   9:   (-1, -1.30000000000000, 0)
   10:  (-1, -1, -0.300000000000000)
   11:  (1, -1.30000000000000, 0)
   12:  (1, -1, -0.300000000000000)
   13:  (sqrt(2), -0.300000000000000, 3)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 3)
   3:   (3, 0)
   4:   (4, 5)
   5:   (5, 6)
   6:   (6, 7)
   7:   (7, 4)
   8:   (0, 4)
   9:   (1, 5)
   10:  (2, 6)
   11:  (3, 7)
   12:  (0, 5)
   13:  (1, 6)
   14:  (2, 7)
   15:  (3, 4)
   16:  (0, 8)
   17:  (0, 9)
   18:  (0, 10)
   19:  (1, 11)
   20:  (1, 12)
   21:  (5, 13)

supports: [8, 9, 10, 11, 12, 13]

joints:
   0:   (-1, -1, 0)
   1:   (1, -1, 0)
   2:   (1, 1, 0)
   3:   (-1, 1, 0)
   4:   (0, -sqrt(2), 3)
   5:   (sqrt(2), 0, 3)
   6:   (0, sqrt(2), 3)
   7:   (-sqrt(2), 0, 3)
   8:   (-1.30000000000000, -1, 0)
   9:   (-1, -1.30000000000000, 0)
   10:  (-1, -1, -0.300000000000000)
   11:  (1, -1.30000000000000, 0)
   12:  (1, -1, -0.300000000000000)
   13:  (sqrt(2), -0.300000000000000, 3)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 3)
   3:   (3, 0)
   4:   (4, 5)
   5:   (5, 6)
   6:   (6, 7)
   7:   (7, 4)
   8:   (0, 4)
   9:   (1, 5)
   10:  (2, 6)
   11:  (3, 7)
   12:  (0, 5)
   13:  (1, 6)
   14:  (2, 7)
   15:  (3, 4)
   16:  (0, 8)
   17:  (0, 9)
   18:  (0, 10)
   19:  (1, 11)
   20:  (1, 12)
   21:  (5, 13)

supports: [8, 9, 10, 11, 12, 13]

plot_truss3d (joints, bars, supports)

maxwells_rule (joints, supports, bars)

3 * 8 - 22  ==  2  !=  0

3 * 8 - 22  ==  2  !=  0

free_joints = others (supports, joints)
free_joints

[0, 1, 2, 3, 4, 5, 6, 7]

[0, 1, 2, 3, 4, 5, 6, 7]

A = equilibrium_matrix (joints, bars, free_joints) 
A

24 x 22 dense matrix over Real Double Field (use the '.str()' method to see the
entries)

24 x 22 dense matrix over Real Double Field (use the '.str()' method to see the entries)

rs, cs, rk, lk = subspaces (A)

print_subspace (cs, sb = False)

Column space:
   Vector space of degree 24 and dimension 21

Column space:
   Vector space of degree 24 and dimension 21

print_subspace (rs, sb = False)

Row space:
   Vector space of degree 22 and dimension 21

Row space:
   Vector space of degree 22 and dimension 21

print_subspace (rk)

Kernel:
   Vector space of degree 22 and dimension 1
   Basis:
     (-0.355464404054382, -0.3554644040543818, -0.3554644040543817,
-0.3554644040543816, -0.3554644040543818, -0.3554644040543817,
-0.3554644040543817, -0.35546440405438184, -0.8016322430806896,
-0.8016322430806891, -0.8016322430806895, -0.8016322430806896, 1.0,
0.9999999999999998, 0.9999999999999999, 1.0, -1.942890293094024e-16,
4.2812207464063204e-17, -3.8096831301583876e-17, 1.8503717077085932e-17,
2.220446049250313e-16, 1.3877787807814457e-16)

Kernel:
   Vector space of degree 22 and dimension 1
   Basis:
     (-0.355464404054382, -0.3554644040543818, -0.3554644040543817, -0.3554644040543816, -0.3554644040543818, -0.3554644040543817, -0.3554644040543817, -0.35546440405438184, -0.8016322430806896, -0.8016322430806891, -0.8016322430806895, -0.8016322430806896, 1.0, 0.9999999999999998, 0.9999999999999999, 1.0, -1.942890293094024e-16, 4.2812207464063204e-17, -3.8096831301583876e-17, 1.8503717077085932e-17, 2.220446049250313e-16, 1.3877787807814457e-16)

bar_forces_ss0 = dict_bar_force (bars, -rk['basis'][0])
print_dict (filter_dict (lambda x : not is_zz_tol()(x), bar_forces_ss0))

0:   0.355464404054382
1:   0.3554644040543818
2:   0.3554644040543817
3:   0.3554644040543816
4:   0.3554644040543818
5:   0.3554644040543817
6:   0.3554644040543817
7:   0.35546440405438184
8:   0.8016322430806896
9:   0.8016322430806891
10:  0.8016322430806895
11:  0.8016322430806896
12:  -1.0
13:  -0.9999999999999998
14:  -0.9999999999999999
15:  -1.0

0:   0.355464404054382
1:   0.3554644040543818
2:   0.3554644040543817
3:   0.3554644040543816
4:   0.3554644040543818
5:   0.3554644040543817
6:   0.3554644040543817
7:   0.35546440405438184
8:   0.8016322430806896
9:   0.8016322430806891
10:  0.8016322430806895
11:  0.8016322430806896
12:  -1.0
13:  -0.9999999999999998
14:  -0.9999999999999999
15:  -1.0

plot_truss_with_forces3d (joints, bars, supports, bar_forces_ss0) . rotateZ (pi/3)

print_subspace (lk, sb = False)

Left kernel:
   Vector space of degree 24 and dimension 3

Left kernel:
   Vector space of degree 24 and dimension 3

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][0], displ_scale = 0.425) . rotateZ (pi/3)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][1], displ_scale = 0.425) . rotateZ (pi/3)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][2], displ_scale = 0.425) . rotateZ (pi/3)

joints = { 0: (1, 1, 3), 1: (-1, 1, 3), 2: (-1, -1, 3), 3: (1, -1, 3), 4: (1, 1, 1), 5: (-1, 1, 1), 6: (-1, -1, 1), 7: (1, -1, 1), 8: (2, 2, 2), 9: (-2, 2, 2), 10: (-2, -2, 2), 11: (2, -2, 2),  12: (2, 2, 0), 13: (-2, 2, 0), 14: (-2, -2, 0), 15: (2, -2, 0),  16: (3, 3, 1), 17: (-3, 3, 1), 18: (-3, -3, 1), 19: (3, -3, 1) }
print 'joints:'
print_dict (joints, indent = '   ')
print

bars = { 0: (0, 1), 1: (1, 2), 2: (2, 3), 3: (3, 0), 4: (4, 5), 5: (5, 6), 6: (6, 7), 7: (7, 4), 8: (0, 4), 9: (1, 5), 10: (2, 6), 11: (3, 7), 12: (0, 8), 13: (4, 8), 14: (1, 9), 15: (5, 9), 16: (2, 10), 17: (6, 10), 18: (3, 11), 19: (7, 11), 20: (12, 13,), 21: (13, 14), 22: (14, 15), 23: (15, 12), 24: (8, 12), 25: (9, 13), 26: (10, 14), 27: (11, 15), 28: (8, 16), 29: (12, 16), 30: (9, 17), 31: (13, 17), 32: (10, 18), 33: (14, 18), 34: (11, 19), 35: (15, 19) }
print 'bars:' 
print_dict (bars, indent = '   ')
print

supports = [16, 17, 18, 19]
print 'supports:', supports
print

joints:
   0:   (1, 1, 3)
   1:   (-1, 1, 3)
   2:   (-1, -1, 3)
   3:   (1, -1, 3)
   4:   (1, 1, 1)
   5:   (-1, 1, 1)
   6:   (-1, -1, 1)
   7:   (1, -1, 1)
   8:   (2, 2, 2)
   9:   (-2, 2, 2)
   10:  (-2, -2, 2)
   11:  (2, -2, 2)
   12:  (2, 2, 0)
   13:  (-2, 2, 0)
   14:  (-2, -2, 0)
   15:  (2, -2, 0)
   16:  (3, 3, 1)
   17:  (-3, 3, 1)
   18:  (-3, -3, 1)
   19:  (3, -3, 1)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 3)
   3:   (3, 0)
   4:   (4, 5)
   5:   (5, 6)
   6:   (6, 7)
   7:   (7, 4)
   8:   (0, 4)
   9:   (1, 5)
   10:  (2, 6)
   11:  (3, 7)
   12:  (0, 8)
   13:  (4, 8)
   14:  (1, 9)
   15:  (5, 9)
   16:  (2, 10)
   17:  (6, 10)
   18:  (3, 11)
   19:  (7, 11)
   20:  (12, 13)
   21:  (13, 14)
   22:  (14, 15)
   23:  (15, 12)
   24:  (8, 12)
   25:  (9, 13)
   26:  (10, 14)
   27:  (11, 15)
   28:  (8, 16)
   29:  (12, 16)
   30:  (9, 17)
   31:  (13, 17)
   32:  (10, 18)
   33:  (14, 18)
   34:  (11, 19)
   35:  (15, 19)

supports: [16, 17, 18, 19]

joints:
   0:   (1, 1, 3)
   1:   (-1, 1, 3)
   2:   (-1, -1, 3)
   3:   (1, -1, 3)
   4:   (1, 1, 1)
   5:   (-1, 1, 1)
   6:   (-1, -1, 1)
   7:   (1, -1, 1)
   8:   (2, 2, 2)
   9:   (-2, 2, 2)
   10:  (-2, -2, 2)
   11:  (2, -2, 2)
   12:  (2, 2, 0)
   13:  (-2, 2, 0)
   14:  (-2, -2, 0)
   15:  (2, -2, 0)
   16:  (3, 3, 1)
   17:  (-3, 3, 1)
   18:  (-3, -3, 1)
   19:  (3, -3, 1)

bars:
   0:   (0, 1)
   1:   (1, 2)
   2:   (2, 3)
   3:   (3, 0)
   4:   (4, 5)
   5:   (5, 6)
   6:   (6, 7)
   7:   (7, 4)
   8:   (0, 4)
   9:   (1, 5)
   10:  (2, 6)
   11:  (3, 7)
   12:  (0, 8)
   13:  (4, 8)
   14:  (1, 9)
   15:  (5, 9)
   16:  (2, 10)
   17:  (6, 10)
   18:  (3, 11)
   19:  (7, 11)
   20:  (12, 13)
   21:  (13, 14)
   22:  (14, 15)
   23:  (15, 12)
   24:  (8, 12)
   25:  (9, 13)
   26:  (10, 14)
   27:  (11, 15)
   28:  (8, 16)
   29:  (12, 16)
   30:  (9, 17)
   31:  (13, 17)
   32:  (10, 18)
   33:  (14, 18)
   34:  (11, 19)
   35:  (15, 19)

supports: [16, 17, 18, 19]

plot_truss3d (joints, bars, supports)

maxwells_rule (joints, supports, bars)

3 * 16 - 36  ==  12  !=  0

3 * 16 - 36  ==  12  !=  0

free_joints = others (supports, joints)
free_joints

[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]

[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]

A = equilibrium_matrix (joints, bars, free_joints) 
A

48 x 36 dense matrix over Real Double Field (use the '.str()' method to see the
entries)

48 x 36 dense matrix over Real Double Field (use the '.str()' method to see the entries)

rs, cs, rk, lk = subspaces (A)

print_subspace (rs, sb = False)

Row space:
   Vector space of degree 36 and dimension 35

Row space:
   Vector space of degree 36 and dimension 35

print_subspace (cs, sb = False)

Column space:
   Vector space of degree 48 and dimension 35

Column space:
   Vector space of degree 48 and dimension 35

print_subspace (rk)

Kernel:
   Vector space of degree 36 and dimension 1
   Basis:
     (0.2886751345948129, 0.2886751345948129, 0.2886751345948129,
0.2886751345948129, 0.2886751345948129, 0.2886751345948129, 0.2886751345948129,
0.2886751345948129, -0.2886751345948129, -0.2886751345948129,
-0.2886751345948129, -0.2886751345948129, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5,
0.5, 0.5773502691896258, 0.5773502691896258, 0.5773502691896258,
0.5773502691896258, -0.5773502691896258, -0.5773502691896258,
-0.5773502691896258, -0.5773502691896258, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
1.0)

Kernel:
   Vector space of degree 36 and dimension 1
   Basis:
     (0.2886751345948129, 0.2886751345948129, 0.2886751345948129, 0.2886751345948129, 0.2886751345948129, 0.2886751345948129, 0.2886751345948129, 0.2886751345948129, -0.2886751345948129, -0.2886751345948129, -0.2886751345948129, -0.2886751345948129, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5773502691896258, 0.5773502691896258, 0.5773502691896258, 0.5773502691896258, -0.5773502691896258, -0.5773502691896258, -0.5773502691896258, -0.5773502691896258, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0)

bar_forces_ss0 = dict_bar_force (bars, rk['basis'][0])
print_dict (bar_forces_ss0)

0:   0.2886751345948129
1:   0.2886751345948129
2:   0.2886751345948129
3:   0.2886751345948129
4:   0.2886751345948129
5:   0.2886751345948129
6:   0.2886751345948129
7:   0.2886751345948129
8:   -0.2886751345948129
9:   -0.2886751345948129
10:  -0.2886751345948129
11:  -0.2886751345948129
12:  0.5
13:  0.5
14:  0.5
15:  0.5
16:  0.5
17:  0.5
18:  0.5
19:  0.5
20:  0.5773502691896258
21:  0.5773502691896258
22:  0.5773502691896258
23:  0.5773502691896258
24:  -0.5773502691896258
25:  -0.5773502691896258
26:  -0.5773502691896258
27:  -0.5773502691896258
28:  1.0
29:  1.0
30:  1.0
31:  1.0
32:  1.0
33:  1.0
34:  1.0
35:  1.0

0:   0.2886751345948129
1:   0.2886751345948129
2:   0.2886751345948129
3:   0.2886751345948129
4:   0.2886751345948129
5:   0.2886751345948129
6:   0.2886751345948129
7:   0.2886751345948129
8:   -0.2886751345948129
9:   -0.2886751345948129
10:  -0.2886751345948129
11:  -0.2886751345948129
12:  0.5
13:  0.5
14:  0.5
15:  0.5
16:  0.5
17:  0.5
18:  0.5
19:  0.5
20:  0.5773502691896258
21:  0.5773502691896258
22:  0.5773502691896258
23:  0.5773502691896258
24:  -0.5773502691896258
25:  -0.5773502691896258
26:  -0.5773502691896258
27:  -0.5773502691896258
28:  1.0
29:  1.0
30:  1.0
31:  1.0
32:  1.0
33:  1.0
34:  1.0
35:  1.0

plot_truss_with_forces3d (joints, bars, supports, bar_forces_ss0)

print_subspace (lk, sb = False)

Left kernel:
   Vector space of degree 48 and dimension 13

Left kernel:
   Vector space of degree 48 and dimension 13

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][10], displ_scale = 0.425)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][0], displ_scale = 0.425)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][3], displ_scale = 0.425)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][4], displ_scale = 0.425)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][5], displ_scale = 0.425)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][7], displ_scale = 0.425)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][11], displ_scale = 0.425)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][8], displ_scale = 0.425)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][1], displ_scale = 0.425)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][2], displ_scale = 0.425)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][6], displ_scale = 0.425)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][9], displ_scale = 0.425)

plot_truss_with_displacements3d (joints, bars, supports, lk['basis'][12], displ_scale = 0.425)

NS_Truss_4_(v0.6)

918 days ago by fresl