A mathematical modeling of a Machining Operation Systems to evaluate stresses generated by the cutting force during machining, was modeled using finite element method. The geometric model was assembled using patran commands. The model was subdivided into elements of quadrilateral surface, where each element has four nodes. Each element is meshed independently and assembled into an elemental stiffness matrix. Each stiffness matrix is then assembled into a global matrix to determine the stress distribution in the model due to the cutting force. For this model, there were 272 elements and 305 nodes. Equal distribution of cutting forces of 1000 Newtons was applied to the model at the area of the cutting tool in the negative Y direction. After applying the forces, the model was analyzed using patran finite element method software P/FEA. The highest stress concentration is at the area where the force is applied. The highest stress value is 4120 Newtons per square meter from the Von Mises stress plot at the area where the cutting tool make contact, and from the displacement stress plot the highest stress value is 0.0000000642 Newtons per square meter.