This work studies optimal shapes of peeling films considering different peeling directions. A cohesive zone model is used to capture normal and tangential forces during peeling. The strip is modeled by nonlinear finite elements based on geometrically exact beam theory. Criteria for optimization are proposed that account for the strain energy during deformation, the overall peeling energy, and the maximum peeling forces. These criteria are first analyzed for various test geometries and then considered for computational optimization. Two cases are considered: Optimization for single direction and optimization for multiple directions. The later allows the determination of strip shapes that provide strong attachment and easy detachment under different peeling angles.