def torsion_stress(T, r, J):
    """ Torsion stresses

    :param float T: Cross section torque
    :param float r: Stress radius coordinate
    :param float J: Cross section polar moment of inertia

    :returns: Torsion Stress at r
    :rtype: float
    """
    return (T * r) / J