A computer-implemented method of simulating the movement of a cloth modeled as a particle model, comprising the steps of: setting deformation units representing parts of interest of the cloth with particles of the particle model positioned at two extremities (Xi, Xj) of a deformation unit, analyzing the movement of cloth in relation to compression by applying an immediate buckling model to the deformation units whereby upon application of a compressive force between the two extremities of a deformation unit, bent deformation of the deformation unit is immediately generated without contractive deformation of the deformation unit being generated, so that a deformation unit go immediately through bent deformation when said compressive force is applied, wherein the step of analyzing comprises the estimation of the shape of the deformation unit after the bent deformation occurs by calculating a curvature of the deformed deformation unit, the curvature is calculated in accordance with solely a distance between the vertexes which are not shared, wherein the compressive force is derived by differentiating a bending deformation energy being accumulated during the application of the immediate buckling model, said bending deformation energy being expressed as a function of the distance between both particles positioned at the extremities of the deformation unit, wherein when the compressive force is below a critical value (Pcr), said compressive force is approximated linearly in relation with the distance between both particles positioned at the extremities of the deformation unit, characterised in that a deformation unit is modeled with two triangular plates (Ai, Aj), each triangular plate having one particle as one of its vertexes, and sharing its other two vertexes with the other triangular plate, and wherein the compressive force vector f i is derived by: f i = k b L 1 A i + A j K 2 cos KL 1 2 − sin c KL 1 2 − 1 x ij 1 x ij 1 wherein k b defines a constant, A i and A j respectively define areas of the triangular plates, K defines the curvature of the deformed deformation unit, x ij 1 defines the distance between both particles positioned at the extremities of the deformation unit after the bent deformation occurs, L 1 defines the distance between both particles positioned extremities of the deformation unit before the bent deformation occurs, and sin c x = sin x x. setting deformation units representing parts of interest of the cloth with particles of the particle model positioned at two extremities (Xi, Xj) of a deformation unit, analyzing the movement of cloth in relation to compression by applying an immediate buckling model to the deformation units whereby upon application of a compressive force between the two extremities of a deformation unit, bent deformation of the deformation unit is immediately generated without contractive deformation of the deformation unit being generated, so that a deformation unit go immediately through bent deformation when said compressive force is applied, wherein the step of analyzing comprises the estimation of the shape of the deformation unit after the bent deformation occurs by calculating a curvature of the deformed deformation unit, the curvature is calculated in accordance with solely a distance between the vertexes which are not shared, wherein the compressive force is derived by differentiating a bending deformation energy being accumulated during the application of the immediate buckling model, said bending deformation energy being expressed as a function of the distance between both particles positioned at the extremities of the deformation unit, wherein when the compressive force is below a critical value (Pcr), said compressive force is approximated linearly in relation with the distance between both particles positioned at the extremities of the deformation unit, characterised in that a deformation unit is modeled with two triangular plates (Ai, Aj), each triangular plate having one particle as one of its vertexes, and sharing its other two vertexes with the other triangular plate, and wherein the compressive force vector f i is derived by: f i = k b L 1 A i + A j K 2 cos KL 1 2 − sin c KL 1 2 − 1 x ij 1 x ij 1 wherein k b defines a constant, A i and A j respectively define areas of the triangular plates, K defines the curvature of the deformed deformation unit, x ij 1 defines the distance between both particles positioned at the extremities of the deformation unit after the bent deformation occurs, L 1 defines the distance between both particles positioned extremities of the deformation unit before the bent deformation occurs, and sin c x = sin x x. A i and A j respectively define areas of the triangular plates, K defines the curvature of the deformed deformation unit, x ij 1 defines the distance between both particles positioned at the extremities of the deformation unit after the bent deformation occurs, L 1 defines the distance between both particles positioned extremities of the deformation unit before the bent deformation occurs, and sin c x = sin x x.