Thanks Plimpton. I've tried that with fix serforce zero as follow.
#fix 1 GB setforce 0 0.0 0.0 ##
#velocity GB zero linear ##
#velocity GB zero angular ##
It can fix the orientation well. But during the NPT simulation,the position still changed.