The paper considers the mathematical model development technique to build a vector field of the shape deviations when machining flat surfaces of shell parts on multi-operational machines under conditions of anisotropic rigidity in technological system (TS). The technological system has an anisotropic rigidity, as its elastic strains do not obey the accepted concepts, i.e. the rigidity towards the coordinate axes of the machine is the same, and they occur only towards the external force. The record shows that the diagrams of elastic strains of machine units are substantially different from the circumference. The issues to ensure the specified accuracy require that there should be mathematical models describing kinematic models and physical processes of mechanical machining under conditions of the specific TS. There are such models for external and internal surfaces of rotation [2,3], which are successfully implemented in practice. Flat surfaces (FS) of shell parts (SP) are both assembly and processing datum surfaces. Therefore, on them special stipulations are made regarding deviations of shape and mutual arrangement. The axes of the main bearing holes are coordinated with respect to them. The joints that ensure leak tightness and distributed load on the product part are closed on these surfaces. The paper deals with the analytical construction of the vector field F, which describes with appropriate approximation the real surface obtained as a result of modeling the process of machining flat surfaces (MFS) through face milling under conditions of anisotropic properties.