Patent ID: 7672821

Claim:
A computer readable storage medium having a computer-executable program to simulate physical or chemical phenomena with one or more non-spherical particles and with or without spherical particles and to calculate translational or rotational displacements of particles, comprising the steps of: (a) judging existence of contact between neighbor particles as existence of contact between mathematically-smooth closed surfaces representing each particle shape and size; and the following (b) to (e) between two particles under the judged contact condition, (b) assuming an imaginary contact point uniquely the same as when both the closed surfaces are reduced isotropically with the same scale to touch at one point; (c) assuming an imaginary contact plane to contain the imaginary contact point and to have a same-directional perpendicular as both normal vectors on both reduced closed surfaces at the imaginary contact point; (d) calculating imaginary contact area using the imaginary contact point and the imaginary contact plane; (e) calculating characteristic values or matrices between both contact particles about stiffness, transmission resistance, transportation resistance or these inverse, using the imaginary contact point, the imaginary contact plane and the imaginary contact area; and (f) based on the characteristic values or matrices between both particles, calculating a characteristic matrix between degrees of freedom of many particles, solving simultaneous linear equations, and calculating translational or rotational displacements or displacement increments of the particles; and the following sub-steps (g) to (j) in the details of the above step (b), (g) assuming a temporary point for the imaginary contact point; (h) reducing or enlarging both closed surfaces isotropically with different scales to pass the assumed temporary point respectively, and calculating both normal vector directions and both local curvatures on the reduced or enlarged both closed surfaces at the assumed temporary point; (i) shifting the position of the assumed temporary point to eliminate the difference of both normal vector directions and the difference of both reducing or enlarging scales using both the local curvatures; and (j) repeating the above sub-steps (h) and (i) as convergence calculation until the difference of both normal vector directions and the difference of both reducing or enlarging scales decrease within allowable values.