B-spline velocity field level set topology optimization method for stress and buckling constraints based on discrete adjoint method.