Fentanyl is a potent and widely used clinical narcotic analgesic, as well as a highly selective μ-opioid agonist. The present study established a homologous model of the human μ-opioid receptor; an intercomparison of three types of μ-opioid receptor protein sequence homologous rates was made. The secondary receptor structure was predicted, the model reliability was assessed and verified using the Ramachandran plot and ProTab analysis. The predictive ability of the CoMFA model was further validated using an extal test set. Using the Surflex-Dock program, a series of fentanyl analog molecules were docked to the receptor, the calculation results from Biopolymer/SiteID showed that the receptor had a deep binding area situated in the extracellular side of the transmembrane domains (TM) among TM3, TM5, TM6, and TM7. Results suggested that there might be 5 active areas in the receptor. The important residues were Asp147, Tyr148, and Tyr149 in TM3, Trp293, and His297 in TM6, and Trp318, His319, Ile322, and Tyr326 in TM7, which were located at the 5 active areas. The best fentanyl docking orientation position was the piperidine ring, which was nearly perpendicular to the membrane surface in the 7 TM domains. Molecular dynamic simulations were applied to evaluate potential relationships between ligand conformation and fentanyl substitution.