mpacts.contact.models.collision. grassmodels

In order to be able to use this module import it like this:

import mpacts.contact.models.collision.grassmodels
#or assign it to a shorter name
import mpacts.contact.models.collision.grassmodels as gra

HertzCundallStrack_DMTAdhesion

Description: Hertz repulsive force with normal damping and Cundall-Strack and a Linear adhesion model.This model serves as a testcase for Tom Leblicq and is not recommended for general use. Geometry combinations available:

PC2 | PC1 - Deformable_Capsule
Deformable_Capsule YES

HertzCundallStrack_DMTAdhesion (Deformable_Capsule Deformable_Capsule)

Hertz repulsive force with normal damping and Cundall-Strack and a Linear adhesion model.This model serves as a testcase for Tom Leblicq and is not recommended for general use.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A1 (s) — Dissipative constant of material of particles from particle container 1.
    • A2 (s) — Dissipative constant of material of particles from particle container 2.
    • E1 (kg . m^-1 . s^-2) — Young’s modulus of the material of particles from particle container 1.
    • E2 (kg . m^-1 . s^-2) — Young’s‘ modulus of the material of particles from particle container 2.
    • attrConst (kg . s^-2) — Adhesion energy between two particles
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • nu1 (1) — Poisson ratio of the material of particles from particle container 1.
    • nu2 (1) — Poisson ratio of the material of particles from particle container 2.
    • pc1 — The first particle container in the binary contact detection.
    • pc2 — The second particle container in the binary contact detection. If contact detection within the same particle container is desired, and it is applicable for the contactmodel, pass the same pc to both pc1 and pc2.
  • Optional keywords:
    • abort_if_different (default value = 0) — If ‘True’, inverts the regular function of ‘AbortIfSameParent’, and makes the contact model early abort if the particles’ parents are different. Please do not change this ‘Property’ if you are not sure what you are doing.
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • torque_factor (default value = 1) — Optional factor to differentially weight the heuristic torque contribution w.r.t. the normal adhesion force. If no additional data is available on the particles you are simulating, please leave this at ‘1’.
  • Read only properties:
    • k — Effective Hertz spring constant, computed from E1, E2, nu1, and nu2. Read only.

This contact model is composed out of following pieces (click on the chain elements to get more information):

DefaultBoilerPlateDeformable_Cylinder_2_DataDeformable_Cylinder_1_DataCapsule_CapsuleN_DampedHertzT_CundallStrackDMT_CapsuleAssembleForces_DeformableCylinder_2AssembleForces_DeformableCylinder_1AbortIfSameParentFeedback

HertzCundallStrack_LinearAdhesion

Description: Hertz repulsive force with normal damping and Cundall-Strack and a Linear adhesion model.This model serves as a testcase for Tom Leblicq and is not recommended for general use. Geometry combinations available:

PC2 | PC1 - Deformable_Capsule
Deformable_Capsule YES

HertzCundallStrack_LinearAdhesion (Deformable_Capsule Deformable_Capsule)

Hertz repulsive force with normal damping and Cundall-Strack and a Linear adhesion model.This model serves as a testcase for Tom Leblicq and is not recommended for general use.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A1 (s) — Dissipative constant of material of particles from particle container 1.
    • A2 (s) — Dissipative constant of material of particles from particle container 2.
    • E1 (kg . m^-1 . s^-2) — Young’s modulus of the material of particles from particle container 1.
    • E2 (kg . m^-1 . s^-2) — Young’s‘ modulus of the material of particles from particle container 2.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • nu1 (1) — Poisson ratio of the material of particles from particle container 1.
    • nu2 (1) — Poisson ratio of the material of particles from particle container 2.
    • pc1 — The first particle container in the binary contact detection.
    • pc2 — The second particle container in the binary contact detection. If contact detection within the same particle container is desired, and it is applicable for the contactmodel, pass the same pc to both pc1 and pc2.
  • Optional keywords:
    • abort_if_different (default value = 0) — If ‘True’, inverts the regular function of ‘AbortIfSameParent’, and makes the contact model early abort if the particles’ parents are different. Please do not change this ‘Property’ if you are not sure what you are doing.
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • d_detach (m) (default value = -1) — Rupture length
    • f_detach (kg . m . s^-2) (default value = -1) — Detachment force
  • Read only properties:
    • k — Effective Hertz spring constant, computed from E1, E2, nu1, and nu2. Read only.

This contact model is composed out of following pieces (click on the chain elements to get more information):

DefaultBoilerPlateDeformable_Cylinder_2_DataDeformable_Cylinder_1_DataCapsule_CapsuleN_DampedHertzT_CundallStrackN_LinearAdhesionForce_CapsulesAssembleForces_DeformableCylinder_2AssembleForces_DeformableCylinder_1AbortIfSameParentFeedbackKeepContacts

HertzDisplacedNGonCapsule

Description: VERY EXPERIMENTAL Geometry combinations available:

PC2 | PC1 - Rigid_Triangle Rigid_Quad Deformable_Triangle
Rigid_Capsule YES YES YES
Deformable_Capsule YES YES YES

HertzDisplacedNGonCapsule (Deformable_Triangle Deformable_Capsule)

VERY EXPERIMENTAL

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A1 — Dissipative constant material 1.
    • A2 — Dissipative constant material 2.
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • nu1 — Poisson ratio material 1.
    • nu2 — Poisson ratio material 2.
    • pc1 — The first particle container in the binary contact detection.
    • pc2 — The second particle container in the binary contact detection. If contact detection within the same particle container is desired, and it is applicable for the contactmodel, pass the same pc to both pc1 and pc2.
  • Optional keywords:
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • displacement_multiplier (default value = 1) — A multiplier applied to the displacement_distance arraystored per node. Basic usage is changing the sign of the displacement.
    • max_force (default value = -1) — maximum force allowed in the contact
    • reject_overlap_r (default value = 1) — Number of radii overlap, at which the contact will be rejected. Set to 0 to disable. Default = 1
  • Read only properties:
    • length_intersection — Length of the cylinder intersection line projected on the Ngon
    • max_overlap — maximal positive overlap on the intersection line.
    • min_overlap — maximal positive overlap on the intersection line.
    • overlap_sphere1 — overlap at the more overlapping sphere.
    • overlap_sphere2 — overlap at the less overlapping sphere.
    • tan_alpha — tangent of the angle between the cylinder axis and the Ngon.

This contact model is composed out of following pieces (click on the chain elements to get more information):

FlatNGonBoilerPlateDeformable_Cylinder_2_DataDeformable_NGon_1_Data  <3 >
  ↓
NGonDisplaced_CapsuleSphere12Hertz_CapsuleT_CundallStrackCapForcesDistributeForcesAndMomentsToFlatTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleSimpleTangentialForces_DeformableTriangle_1AssembleForces_DeformableCylinder_2Feedback

HertzDisplacedNGonCapsule (Deformable_Triangle Rigid_Capsule)

VERY EXPERIMENTAL

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A1 — Dissipative constant material 1.
    • A2 — Dissipative constant material 2.
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • nu1 — Poisson ratio material 1.
    • nu2 — Poisson ratio material 2.
    • pc1 — The first particle container in the binary contact detection.
    • pc2 — The second particle container in the binary contact detection. If contact detection within the same particle container is desired, and it is applicable for the contactmodel, pass the same pc to both pc1 and pc2.
  • Optional keywords:
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • displacement_multiplier (default value = 1) — A multiplier applied to the displacement_distance arraystored per node. Basic usage is changing the sign of the displacement.
    • max_force (default value = -1) — maximum force allowed in the contact
    • reject_overlap_r (default value = 1) — Number of radii overlap, at which the contact will be rejected. Set to 0 to disable. Default = 1
  • Read only properties:
    • length_intersection — Length of the cylinder intersection line projected on the Ngon
    • max_overlap — maximal positive overlap on the intersection line.
    • min_overlap — maximal positive overlap on the intersection line.
    • overlap_sphere1 — overlap at the more overlapping sphere.
    • overlap_sphere2 — overlap at the less overlapping sphere.
    • tan_alpha — tangent of the angle between the cylinder axis and the Ngon.

This contact model is composed out of following pieces (click on the chain elements to get more information):

FlatNGonBoilerPlateRigid_Cylinder_2_DataDeformable_NGon_1_Data  <3 >
  ↓
NGonDisplaced_CapsuleSphere12Hertz_CapsuleT_CundallStrackCapForcesDistributeForcesAndMomentsToFlatTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleSimpleTangentialForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Feedback

HertzDisplacedNGonCapsule (Rigid_Quad Deformable_Capsule)

VERY EXPERIMENTAL

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A1 — Dissipative constant material 1.
    • A2 — Dissipative constant material 2.
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • nu1 — Poisson ratio material 1.
    • nu2 — Poisson ratio material 2.
    • pc1 — The first particle container in the binary contact detection.
    • pc2 — The second particle container in the binary contact detection. If contact detection within the same particle container is desired, and it is applicable for the contactmodel, pass the same pc to both pc1 and pc2.
  • Optional keywords:
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • displacement_multiplier (default value = 1) — A multiplier applied to the displacement_distance arraystored per node. Basic usage is changing the sign of the displacement.
    • max_force (default value = -1) — maximum force allowed in the contact
    • reject_overlap_r (default value = 1) — Number of radii overlap, at which the contact will be rejected. Set to 0 to disable. Default = 1
  • Read only properties:
    • length_intersection — Length of the cylinder intersection line projected on the Ngon
    • max_overlap — maximal positive overlap on the intersection line.
    • min_overlap — maximal positive overlap on the intersection line.
    • overlap_sphere1 — overlap at the more overlapping sphere.
    • overlap_sphere2 — overlap at the less overlapping sphere.
    • tan_alpha — tangent of the angle between the cylinder axis and the Ngon.

This contact model is composed out of following pieces (click on the chain elements to get more information):

FlatNGonBoilerPlateDeformable_Cylinder_2_DataRigid_NGon_1_Data  <4 >
  ↓
NGonDisplaced_CapsuleSphere12Hertz_CapsuleT_CundallStrackCapForcesAssembleForces_DeformableCylinder_2AssembleMoments_1AssembleForces_1Feedback

HertzDisplacedNGonCapsule (Rigid_Quad Rigid_Capsule)

VERY EXPERIMENTAL

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A1 — Dissipative constant material 1.
    • A2 — Dissipative constant material 2.
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • nu1 — Poisson ratio material 1.
    • nu2 — Poisson ratio material 2.
    • pc1 — The first particle container in the binary contact detection.
    • pc2 — The second particle container in the binary contact detection. If contact detection within the same particle container is desired, and it is applicable for the contactmodel, pass the same pc to both pc1 and pc2.
  • Optional keywords:
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • displacement_multiplier (default value = 1) — A multiplier applied to the displacement_distance arraystored per node. Basic usage is changing the sign of the displacement.
    • max_force (default value = -1) — maximum force allowed in the contact
    • reject_overlap_r (default value = 1) — Number of radii overlap, at which the contact will be rejected. Set to 0 to disable. Default = 1
  • Read only properties:
    • length_intersection — Length of the cylinder intersection line projected on the Ngon
    • max_overlap — maximal positive overlap on the intersection line.
    • min_overlap — maximal positive overlap on the intersection line.
    • overlap_sphere1 — overlap at the more overlapping sphere.
    • overlap_sphere2 — overlap at the less overlapping sphere.
    • tan_alpha — tangent of the angle between the cylinder axis and the Ngon.

This contact model is composed out of following pieces (click on the chain elements to get more information):

FlatNGonBoilerPlateRigid_Cylinder_2_DataRigid_NGon_1_Data  <4 >
  ↓
NGonDisplaced_CapsuleSphere12Hertz_CapsuleT_CundallStrackCapForcesAssembleMoments_2AssembleForces_2AssembleMoments_1AssembleForces_1Feedback

HertzDisplacedNGonCapsule (Rigid_Triangle Deformable_Capsule)

VERY EXPERIMENTAL

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A1 — Dissipative constant material 1.
    • A2 — Dissipative constant material 2.
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • nu1 — Poisson ratio material 1.
    • nu2 — Poisson ratio material 2.
    • pc1 — The first particle container in the binary contact detection.
    • pc2 — The second particle container in the binary contact detection. If contact detection within the same particle container is desired, and it is applicable for the contactmodel, pass the same pc to both pc1 and pc2.
  • Optional keywords:
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • displacement_multiplier (default value = 1) — A multiplier applied to the displacement_distance arraystored per node. Basic usage is changing the sign of the displacement.
    • max_force (default value = -1) — maximum force allowed in the contact
    • reject_overlap_r (default value = 1) — Number of radii overlap, at which the contact will be rejected. Set to 0 to disable. Default = 1
  • Read only properties:
    • length_intersection — Length of the cylinder intersection line projected on the Ngon
    • max_overlap — maximal positive overlap on the intersection line.
    • min_overlap — maximal positive overlap on the intersection line.
    • overlap_sphere1 — overlap at the more overlapping sphere.
    • overlap_sphere2 — overlap at the less overlapping sphere.
    • tan_alpha — tangent of the angle between the cylinder axis and the Ngon.

This contact model is composed out of following pieces (click on the chain elements to get more information):

FlatNGonBoilerPlateDeformable_Cylinder_2_DataRigid_NGon_1_Data  <3 >
  ↓
NGonDisplaced_CapsuleSphere12Hertz_CapsuleT_CundallStrackCapForcesAssembleForces_DeformableCylinder_2AssembleMoments_1AssembleForces_1Feedback

HertzDisplacedNGonCapsule (Rigid_Triangle Rigid_Capsule)

VERY EXPERIMENTAL

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A1 — Dissipative constant material 1.
    • A2 — Dissipative constant material 2.
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • nu1 — Poisson ratio material 1.
    • nu2 — Poisson ratio material 2.
    • pc1 — The first particle container in the binary contact detection.
    • pc2 — The second particle container in the binary contact detection. If contact detection within the same particle container is desired, and it is applicable for the contactmodel, pass the same pc to both pc1 and pc2.
  • Optional keywords:
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • displacement_multiplier (default value = 1) — A multiplier applied to the displacement_distance arraystored per node. Basic usage is changing the sign of the displacement.
    • max_force (default value = -1) — maximum force allowed in the contact
    • reject_overlap_r (default value = 1) — Number of radii overlap, at which the contact will be rejected. Set to 0 to disable. Default = 1
  • Read only properties:
    • length_intersection — Length of the cylinder intersection line projected on the Ngon
    • max_overlap — maximal positive overlap on the intersection line.
    • min_overlap — maximal positive overlap on the intersection line.
    • overlap_sphere1 — overlap at the more overlapping sphere.
    • overlap_sphere2 — overlap at the less overlapping sphere.
    • tan_alpha — tangent of the angle between the cylinder axis and the Ngon.

This contact model is composed out of following pieces (click on the chain elements to get more information):

FlatNGonBoilerPlateRigid_Cylinder_2_DataRigid_NGon_1_Data  <3 >
  ↓
NGonDisplaced_CapsuleSphere12Hertz_CapsuleT_CundallStrackCapForcesAssembleMoments_2AssembleForces_2AssembleMoments_1AssembleForces_1Feedback