def delta_qu(soil, H, D):
    """ Displacement at Qu
    """
    if 'sand' in soil:
        return min(0.01 * H, 0.1 * D)
    elif 'clay' in soil:
        return min(0.1 * H, 0.2 * D)
    else:
        raise ValueError('Unknown soil type.')
