Systems and methods of analysis of granular elements

A method of analyzing the behavior of a particulate system comprising a plurality of particles, where each particle possesses a particle morphology, the method comprising:

• generating a Non-Uniform Rational Basis Spline (NURBS) description for the particle morphology for each of a plurality of particles in three dimensions (3D) within a particulate system using a processing system comprising a processor and memory, where the NURBS description of a given particle within the particulate system comprises a plurality of control points that define a control polygon bounding the geometry of a NURBS curve, a knot vector that determines the location of a plurality of knots, a basis function of a predetermined degree that determines the number of control points that influence any given point on the NURBS curve, and a plurality of weights that determine the effective contribution of a control point on the overall shape of the NURBS curve;
• determining contact points between at least two particles within the particulate system using the processor of the processing system based on the NURBS description for the particle morphology of the at least two particles contained within the memory of the processing system;
• determining the magnitude of the contact between the at least two particles within the particulate system using the processor of the processing system further comprises:
• identifying pairs of intersecting portions of the surfaces of two particles based on the NURBS descriptions of the particle morphology of the two particles;
• determining a closed intersection loop for the intersecting portions of the surfaces of two particles;
• approximating a loop centroid and a loop normal for the closed intersection loop;
• determining contact points between a line through the loop centroid extending in the direction of the loop normal and the surfaces of each of the two particles;
• calculating an overlap distance between the two particles based on the contact points between the line through the loop centroid extending in the direction of the loop normal and the surfaces of each of the two particles;
• and
• quantifying the magnitude of contact based on the overlap distance between the two particles;
• determining normal forces and associated moments using the processor of the processing system based upon the contact points and the magnitude of the contact between the at least two particles contained within the memory of the processing system;
• determining tangential forces and associated moments using the processor of the processing system based upon data describing the contact points and the magnitude of the contact between the at least two particles contained within the memory of the processing system;
• and
• performing mechanical analysis of the particulate system using the processor of the processing system based on the data describing contact between the at least two particles and the resulting forces and associated moments contained within the memory of the processing system.

A method of analyzing the behavior of a particulate system comprising a plurality of particles, where each particle possesses a particle morphology, the method comprising:

• generating a Non-Uniform Rational Basis Spline (NURBS) description for the particle morphology of each of a plurality of particles within a particulate system in two dimensions using a processing system comprising a processor and memory, where the NURBS description of a given particle within the particulate system comprises a plurality of control points that define a control polygon bounding the geometry of a NURBS curve, a knot vector that determines the location of a plurality of knots, a basis function of a predetermined degree that determines the number of control points that influence any given point on the NURBS curve, and a plurality of weights that determine the effective contribution of a control point on the overall shape of the NURBS curve;
• determining contact points between at least two particles within the particulate system using the processor of the processing system based on the NURBS description for the particle morphology of the at least two particles contained within the memory of the processing system;
• determining a magnitude of the contact between the at least two particles within the particulate system using the processor of the processing system based on the NURBS description for the particle morphology of the at least two particles contained within the memory of the processing system by:
• determining two intersection points of contact between two particles within the particulate system;
• calculating a midpoint for the determined intersection points of contact;
• determining normal and tangential unit vectors using the calculated midpoint;
• determining intersection points between a line extending in both directions of the normal vector line and the boundaries of the two particles;
• approximating an overlap distance between the two particles based on the intersection points between a line extending in both directions of the normal vector line and the boundaries of the two particles;
• and
• quantifying the magnitude of contact using the overlap distance between the two particles;
• determining normal forces and associated moments using the processor of the processing system based upon data describing the contact points and the magnitude of the contact between the at least two particles contained within the memory of the processing system by:
• determining a normal elastic stiffness for a contact between two particles;
• determining a contact damping parameter for the contact between the two particles;
• calculating a normal component of relative velocity with respect to the two particles;
• calculating a normal force using the normal elastic stiffness, contact damping parameter, and normal component of relative velocity for each of the two particles;
• and
• computing a moment arm for each of the two particles;
• and
• calculating an associated moment due to the calculated normal force for each of the two particles;
• determining tangential forces and associated moments using the processor of the processing system based upon data describing the contact points and the magnitude of the contact between the at least two particles contained within the memory of the processing system by:
• determining angular velocities for each of the two particles;
• calculating a relative tangential velocity for each of the two particles;
• calculating an increment in tangential force using the relative tangential velocity for each of the two particles;
• calculating a tangential force at a discrete time state using the calculated increment in tangential force for each of the two particles;
• and
• calculating an associated moment due to the tangential force for each of the two particles;
• and
• performing mechanical analysis of the particulate system using the processor of the processing system based on data describing the contact between the at least two particles and the resulting forces and associated moments contained within the memory of the processing system.

A method of analyzing the behavior of a particulate system comprising a plurality of particles, where each particle possesses a particle morphology, the method comprising:

• generating a Non-Uniform Rational Basis Spline (NURBS) description for the particle morphology of each of a plurality of particles within a particulate system in three dimensions using a processing system comprising a processor and memory, where the NURBS description of a given particle within the particulate system comprises a plurality of control points that define a control polygon bounding the geometry of a NURBS curve, a knot vector that determines the location of a plurality of knots, a basis function of a predetermined degree that determines the number of control points that influence any given point on the NURBS curve, and a plurality of weights that determine the effective contribution of a control point on the overall shape of the NURBS curve;
• determining contact points between at least two particles within the particulate system using the processor of the processing system based on the NURBS description for the particle morphology of the at least two particles contained within the memory of the processing system;
• determining a magnitude of the contact between the at least two particles within the particulate system using the processor of the processing system based on the NURBS description for the particle morphology of the at least two particles using the processing system by
• identifying pairs of intersecting portions of the surfaces of two particles based on the NURBS descriptions of the particle morphology of the two particles;
• determining a closed intersection loop for the intersecting portions of the surfaces of two particles;
• approximating a loop centroid and a loop normal for the closed intersection loop;
• determining contact points between a line through the loop centroid extending in the direction of the loop normal and the surfaces of each of the two particles;
• calculating an overlap distance between the two particles based on the contact points between the line through the loop centroid extending in the direction of the loop normal and the surfaces of each of the two particles;
• and
• quantifying the magnitude of contact based on the overlap distance between the two particles;
• determining normal forces and associated moments using the processor of the processing system based upon the contact points and the magnitude of the contact between the at least two particles contained within the memory of the processing system by:
• determining a normal elastic stiffness for a contact between two particles;
• determining the overlap distance between the two particles;
• calculating a normal force for each of the two particles using the normal elastic stiffness and the overlap distance between the two particles;
• and
• calculating an associated moment due to the calculated normal force for each of the two particles;
• determining tangential forces and associated moments using the processor of the processing system based upon the contact points and the magnitude of the contact between the at least two particles contained within the memory of the processing system:
• determining angular velocities for each of the two particles;
• determining translational velocities for each of the two particles;
• calculating an increment in tangential force for each of the two particles;
• calculating a tangential force accounting for incremental rotation of the contact plane for each of the two particles;
• and
• calculating an associated moment due to the tangential force for each of the two particles;
• and
• performing mechanical analysis of the particulate system using the processor of the processing system based on data describing the contact between the at least two particles and the resulting forces and associated moments contained within the memory of the processing system.