The use of graduated compression stockings (GCS) with the aim of prevention and treatment of venous disorders is wide spread. Advantages of the use of these textiles in the medical field are universally proven. However, their exact functional performance in respect of contact of these textile structures with the body and the mechanism by which treatment is achieved is not fully understood. Direct pressure measurement using available tools cannot yield to a comprehensive understanding of the interaction between stockings and the body. This work aims to develop a 3D biomechanical mathematical model that can provide a dynamical numerical simulation of the interaction between the leg and the compression stockings. This model is based on an actual geometry of female leg. The geometry was obtained by reconstruction of 3D CT images and the real size together with the mechanical properties of two commercially available stockings. The biomechanical solid leg model was composed of tibia and fibula, soft tissues. An orthotropic shell model was assumed for the stockings. The wearing of the stockings was dynamically simulated. This was based on relative sliding contact phenomenon between two objects. The displacement of the stockings hemline was set at 280mm. The surface pressure distribution, pressure magnitude, tension and stress were simulated using two commercially available compression stockings in pressure ltr"