line.py 文件源码

def stiffness_global(x, y, z, E, G, Ax, Iz=0, Iy=0, Ay=0, Az=0, theta=0, J=0):
    t = block_diag(*(transformMatrix(x, y, z, theta),) * 2)
    r = t.transpose()
    for m in stiffness_local(norm((x, y, z)), E, G, Ax, Iz, Iy, Ay, Az, J):
        yield dot(dot(t, m), r)