This　paper　aims　to　propose　a　topology　optimization　method　on　generating　porous　structures　comprising　multiple　materials.　The　mathematical　optimization　formulation　is　established　under　the　constraints　of　individual　volume　fraction　of　constituent　phase　or　total　mass,　as　well　as　the　local　volume　fraction　of　all　phases.　The　original　optimization　problem　with　numerous　constraints　is　converted　into　a　box-constrained　optimization　problem　by　incorporating　all　constraints　to　the　augmented　Lagrangian　function,　avoiding　the　parameter　dependence　in　the　conventional　aggregation　process.　Furthermore,　the　local　volume　percentage　can　be　precisely　satisfied.　The　effects　including　the　global　mass　bound,　the　influence　radius　and　local　volume　percentage　on　final　designs　are　exploited　through　numerical　examples.　The　numerical　results　also　reveal　that　porous　structures　keep　a　balance　between　the　bulk　design　and　periodic　design　in　terms　of　the　resulting　compliance.　All　results,　including　those　for　irregular　structures　and　multiple　volume　fraction　constraints,　demonstrate　that　the　proposed　method　can　provide　an　efficient　solution　for　multiple　material　infill　structures.