Use trigsimp(expr, method='fu') instead of simplify. I would also consider changing the Vector.simplify flag to a global mechanics simplify flag.