PFC Overview
The PFC programs (PFC2D and PFC3D) provide a general purpose, distinct-element modeling framework that includes both a computational engine and a graphical user interface. The PFC model (refers to both 2D and 3D models) simulates the movement and interaction of many finite-sized particles. [1][2] The particles are rigid bodies with finite mass that move independently of one another and can both translate and rotate. Particles interact at pair-wise contacts by means of an internal force and moment. Contact mechanics is embodied in particle-interaction laws that update the inter-particle forces and moments. The time evolution of this system is determined by Newton’s laws of motion.
Thanks to its general design, PFC can be easily customized and applied to a very broad range of numerical investigations where the discrete nature of the systems is of interest. Since the release of the first version in 1994, PFC has been successfully used by many academic institutions and private companies around the world for problems ranging from fundamental research on soil and rock behavior at the laboratory scale to slope stability and rockfall hazard mitigation, hydraulic fracturing, rock-tool interactions, bulk flow, mixing, conveying and compaction of aggregates and powders, blast furnace modeling, etc. A large panel of references can be found in the proceedings of the past PFC and FLAC/DEM symposia periodically organized by Itasca.
The PFC Model
A general PFC model simulates the mechanical behavior of a system made up of a collection of arbitrarily shaped particles. (Note that in the present context, the term particle denotes a body that occupies a finite amount of space, as opposed to a single point in space). The model is composed of distinct particles that displace independent of one another and interact at pair-wise contacts. The particles are assumed to be rigid and the mechanical behavior of such a system is described in terms of the movement of each particle and the inter-particle forces acting at each contact point. The relationship between particle motion and the forces and moments causing that motion is computed via explicit dynamic integration of Newton’s laws of motion. The force system may be in static equilibrium (in which case, there is no motion), or it may be such as to cause the particles to flow. The particle-interaction laws that model physical contact between particles use a soft contact approach, in which a finite normal stiffness is taken to represent the measurable stiffness that exists at a contact, and the rigid particles are allowed to overlap in the vicinity of the contact point. More complex behavior can be modeled by allowing the particles to be bonded together at their contact points such that tensile forces can develop between particles. When the inter-particle forces acting at any bond exceed the bond strength, that bond is broken. One can then model the interaction of these bonded particles including the formation of cracks that may cause blocks to fragment into smaller blocks. The particle-interaction law can also be derived from potential energy functions and model long range interactions. The PFC particles also interact with walls, which can be used as fixed or moving boundaries for purposes of compaction, confinement, mixing or conveying of the particulate assemblies.
The assumption of particle rigidity is a good one when most of the deformation in a physical system is accounted for by the local deformation of the particles and movements along their interfaces. The deformation of a packed-particle assembly (or a granular assembly such as sand), as a whole, is well-described by this assumption, since the deformation results primarily from the sliding and rotation of the particles as rigid bodies and the opening and interlocking at interfaces, not from the bulk deformation of individual particles. Precise modeling of particle deformation is not necessary to obtain a good approximation of the mechanical behavior for such systems.
In addition to traditional particle-flow applications, PFC can also be applied to the analysis of solids subjected to prescribed boundary and initial conditions. In such models, the continuum behavior is approximated by treating the solid as a compacted assembly of many small particles. Measures of stress and strain rate can be defined as average quantities over a representative measurement volume for such systems. This allows one to estimate interior stresses for granular materials such as soils, or solid materials such as rock or plastics formed by powder compaction.
Distinct-Element Method
PFC models the movement and interaction of stressed assemblies of rigid {circular in 2D; spherical in 3D} particles using the Distinct-Element Method (DEM). The DEM was introduced by [Cundall1971] for the analysis of rock-mechanics problems and then applied to soils by [Cundall1979b]. A thorough description of the method is given in the two-part paper by [Cundall1988] and [Hart1988], and in the UDEC manual ([Itasca2011]). PFC is classified as a discrete element code based on the definition in the review by [Cundall1992], since it allows finite displacements and rotations of discrete bodies, and recognizes new contacts automatically as the calculation progresses. PFC can be viewed as a simplified implementation of the DEM because of the restriction to rigid particles. The general DEM can handle deformable polygonal-shaped particles.
In the DEM, the interaction of the particles is treated as a dynamic process with states of equilibrium developing whenever the internal forces balance. The contact forces and displacements of a stressed assembly of particles are found by tracing the movements of the individual particles. Movements result from the propagation through the particle system of disturbances caused by specified wall and particle motion and/or body forces. This is a dynamic process in which the speed of propagation depends on the physical properties of the discrete system.
The dynamic behavior is represented numerically by a timestepping algorithm in which it is assumed that the velocities and accelerations are constant within each timestep. The solution scheme is identical to that used by the explicit finite-difference method for continuum analysis. The DEM is based upon the idea that the timestep chosen may be so small that, during a single timestep, disturbances cannot propagate further from any particle than its immediate neighbors. Then, at all times, the forces acting on any particle are determined exclusively by its interaction with the particles with which it is in contact. Since the speed at which a disturbance can propagate is a function of the physical properties of the discrete system, the timestep can be chosen to satisfy the above constraint. The use of an explicit, as opposed to an implicit, numerical scheme makes it possible to simulate the nonlinear interaction of a large number of particles without excessive memory requirements or the need for an iterative procedure.
The calculations performed in the DEM alternate between the application of Newton’s second law to the particles and a force-displacement law at the contacts. Newton’s second law is used to determine the motion of each particle arising from the contact and body forces acting upon it, while the force-displacement law is used to update the contact forces arising from the relative motion at each contact. The presence of walls in PFC requires only that the force-displacement law account for contacts with wall facets. Newton’s second law is not applied to walls because the wall motion is specified by the user.
References
Cundall, P. A. “A Computer Model for Simulating Progressive Large Scale Movements in Blocky Rock Systems,” in Proceedings of the Symposium of the International Society of Rock Mechanics (Nancy, France, 1971), Vol. 1, Paper No. II-8 (1971).
Cundall, P. A. “Formulation of a Three-Dimensional Distinct Element Model - Part I. A Scheme to Detect and Represent Contacts in a System Composed of Many Polyhedral Blocks,” Int. J. Rock Mech., Min. Sci. & Geomech. Abstr., 25(3), 107-116 (1988).
Cundall, P. A., and R. Hart. “Numerical Modeling of Discontinua,” J. Engr. Comp., 9, 101-113 (1992).
Cundall, P. A., and O. D. L. Strack. “A Discrete Numerical Model for Granular Assemblies,” Geotechnique, 29(1), 47-65 (1979).
Hart, R., P. A. Cundall and J. Lemos. “Formulation of a Three-Dimensional Distinct Element Model - Part II. Mechanical Calculations for Motion and Interaction of a System Composed of Many Polyhedral Blocks,” Int. J. Rock Mech., Min. Sci. & Geomech. Abstr., 25(3), 117-125 (1988).
Itasca Consulting Group Inc. UDEC (Universal Distinct Element Code), Version 5.0. Minneapolis:ICG, 2011.
Endnotes
| Was this helpful? ... | Itasca Software © 2026, Itasca | Updated: May 13, 2026 |
