Quality design methodologies have received constituent attention from a number of researchers and practitioners for more than twenty years. Specially, the quality design for drug products must be carefully considered because of the hazards involved in the pharmaceutical industry. Conventional pharmaceutical formulation design problems with mixture experiments have been typically studied under the assumption of an unconstrained experimental region with a single quality characteristic. However, real-world pharmaceutical industrial situations have many physical limitations. We are often faced with multiple quality characteristics with constrained experimental regions. ln order to address these issues, the main objective of this paper is to propose a robust desirability function (RDF) model using a desirability function (DF) and mean square error (MSE) to simultaneously consider a number of multiple quality characteristics. This paper then present L-pseudocomponents and U-pseudocomponents to handle physical constraints. Finally, a numerical example shows that the proposed RDF can efficiently be applied to a pharmaceutical process design.