An efficient, error-bounded approximation algorithm for simulating quasi-statics of complex linkages