Nano-friction of phosphorene plays a significant role in affecting the controllability and efficiency of applying strain engineering to tune its properties. So far, the friction behavior of phosphorene has not been studied. This work studies the friction of single-layer and bilayer phosphorene on an amorphous silicon substrate by sliding a rigid tip. For the single-layer phosphorene, it is found that its friction is highly anisotropic, i.e. the friction is larger along the armchair direction than that along the zigzag direction. Moreover, pre-strain of the phosphorene also exhibits anisotropic effects. The friction increases with the pre-strain along the zigzag direction, but decreases with that along the armchair direction. Furthermore, the strong adhesion between the phosphorene and its substrate increases the friction between the phosphorene and the tip. For bilayer phosphorene, its friction highly depends on its stacking mode, which determines the contact interface with a commensurate or incommensurate pattern. This friction behavior is quite unique and greatly differs from that of graphene and molybdenum disulfide. Detailed analysis reveals that this behavior results from the combination effect of the friction contact area, the potential-energy profile of phosphorene, and its unique elongation.