mpacts.contact.models.collision. stemcompression

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

import mpacts.contact.models.collision.stemcompression
#or assign it to a shorter name
import mpacts.contact.models.collision.stemcompression as ste

CompressionModelLinearForceCylinder

Description: Linear damper for LUT Geometry combinations available:

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

CompressionModelLinearForceCylinder (Rigid_Quad Deformable_Capsule)

Linear damper for LUT

Parallel Compatible: No

Properties:

  • Required keywords:
    • c — viscous damper constant.
    • lut — 2D look up table that will be used to calculate normal forces
    • mu — The coulomb friction coefficient (both static and dynamic).
    • 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 = -1) — ‘c’ value of the linear dashpot (N*s/m) in tangential direction. The higher this value, the more accurate the results will be, but the simulation can become unstable, requiring smaller timesteps.
    • cp_t (default value = -1) — Optional contact area-dependent linear dashpot coefficent (Pa*s/m). Give either c_t or cp_t but not both.
    • 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):

DefaultBoilerPlateDeformable_Cylinder_2_DataRigid_NGon_1_Data  <4 >
  ↓
NGon_CapsuleSphere12LUTLineForce_CylindersT_CoulombFrictionAssembleForces_DeformableCylinder_2AssembleMoments_1AssembleForces_1Feedback

CompressionModelLinearForceCylinder (Rigid_Quad Rigid_Capsule)

Linear damper for LUT

Parallel Compatible: No

Properties:

  • Required keywords:
    • c — viscous damper constant.
    • lut — 2D look up table that will be used to calculate normal forces
    • mu — The coulomb friction coefficient (both static and dynamic).
    • 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 = -1) — ‘c’ value of the linear dashpot (N*s/m) in tangential direction. The higher this value, the more accurate the results will be, but the simulation can become unstable, requiring smaller timesteps.
    • cp_t (default value = -1) — Optional contact area-dependent linear dashpot coefficent (Pa*s/m). Give either c_t or cp_t but not both.
    • 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):

DefaultBoilerPlateRigid_Cylinder_2_DataRigid_NGon_1_Data  <4 >
  ↓
NGon_CapsuleSphere12LUTLineForce_CylindersT_CoulombFrictionAssembleForcesAssembleMomentsFeedback

CompressionModelLinearForceCylinder (Rigid_Triangle Deformable_Capsule)

Linear damper for LUT

Parallel Compatible: No

Properties:

  • Required keywords:
    • c — viscous damper constant.
    • lut — 2D look up table that will be used to calculate normal forces
    • mu — The coulomb friction coefficient (both static and dynamic).
    • 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 = -1) — ‘c’ value of the linear dashpot (N*s/m) in tangential direction. The higher this value, the more accurate the results will be, but the simulation can become unstable, requiring smaller timesteps.
    • cp_t (default value = -1) — Optional contact area-dependent linear dashpot coefficent (Pa*s/m). Give either c_t or cp_t but not both.
    • 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):

DefaultBoilerPlateDeformable_Cylinder_2_DataRigid_NGon_1_Data  <3 >
  ↓
NGon_CapsuleSphere12LUTLineForce_CylindersT_CoulombFrictionAssembleForces_DeformableCylinder_2AssembleMoments_1AssembleForces_1Feedback

CompressionModelLinearForceCylinder (Rigid_Triangle Rigid_Capsule)

Linear damper for LUT

Parallel Compatible: No

Properties:

  • Required keywords:
    • c — viscous damper constant.
    • lut — 2D look up table that will be used to calculate normal forces
    • mu — The coulomb friction coefficient (both static and dynamic).
    • 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 = -1) — ‘c’ value of the linear dashpot (N*s/m) in tangential direction. The higher this value, the more accurate the results will be, but the simulation can become unstable, requiring smaller timesteps.
    • cp_t (default value = -1) — Optional contact area-dependent linear dashpot coefficent (Pa*s/m). Give either c_t or cp_t but not both.
    • 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):

DefaultBoilerPlateRigid_Cylinder_2_DataRigid_NGon_1_Data  <3 >
  ↓
NGon_CapsuleSphere12LUTLineForce_CylindersT_CoulombFrictionAssembleForcesAssembleMomentsFeedback

CompressionModel_LUT_Collection_stem_plate

Description: Compression using two LUTs Geometry combinations available:

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

CompressionModel_LUT_Collection_stem_plate (Rigid_Quad Deformable_Capsule)

Compression using two LUTs

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • c_plate — viscous damper constant.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • 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:
    • Fprim1 (default value = None) — Array with vectors which stores the contact force per primitive for pc1. If not given, the array pc1[‘Fprim’] is searched first, and if not found, nothing will be done
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • 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):

DefaultBoilerPlateDeformable_Cylinder_2_DataRigid_NGon_1_Data  <4 >
  ↓
NGon_CapsuleSphere12N_LUT_Collection_stem_plate_NO_UPDATE_REMOVEASAPT_CundallStrackAssembleForces_Primitives_1AssembleForces_DeformableCylinder_2AssembleMoments_1AssembleForces_1Feedback

CompressionModel_LUT_Collection_stem_plate (Rigid_Quad Rigid_Capsule)

Compression using two LUTs

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • c_plate — viscous damper constant.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • 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:
    • Fprim1 (default value = None) — Array with vectors which stores the contact force per primitive for pc1. If not given, the array pc1[‘Fprim’] is searched first, and if not found, nothing will be done
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • 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):

DefaultBoilerPlateRigid_Cylinder_2_DataRigid_NGon_1_Data  <4 >
  ↓
NGon_CapsuleSphere12N_LUT_Collection_stem_plate_NO_UPDATE_REMOVEASAPT_CundallStrackAssembleForces_Primitives_1AssembleForcesAssembleMomentsFeedback

CompressionModel_LUT_Collection_stem_plate (Rigid_Triangle Deformable_Capsule)

Compression using two LUTs

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • c_plate — viscous damper constant.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • 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:
    • Fprim1 (default value = None) — Array with vectors which stores the contact force per primitive for pc1. If not given, the array pc1[‘Fprim’] is searched first, and if not found, nothing will be done
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • 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):

DefaultBoilerPlateDeformable_Cylinder_2_DataRigid_NGon_1_Data  <3 >
  ↓
NGon_CapsuleSphere12N_LUT_Collection_stem_plate_NO_UPDATE_REMOVEASAPT_CundallStrackAssembleForces_Primitives_1AssembleForces_DeformableCylinder_2AssembleMoments_1AssembleForces_1Feedback

CompressionModel_LUT_Collection_stem_plate (Rigid_Triangle Rigid_Capsule)

Compression using two LUTs

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • c_plate — viscous damper constant.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • 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:
    • Fprim1 (default value = None) — Array with vectors which stores the contact force per primitive for pc1. If not given, the array pc1[‘Fprim’] is searched first, and if not found, nothing will be done
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • 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):

DefaultBoilerPlateRigid_Cylinder_2_DataRigid_NGon_1_Data  <3 >
  ↓
NGon_CapsuleSphere12N_LUT_Collection_stem_plate_NO_UPDATE_REMOVEASAPT_CundallStrackAssembleForces_Primitives_1AssembleForcesAssembleMomentsFeedback

CompressionModel_LUT_Collection_stem_plate_and_weighed_Hertz_Fprim

Description: CompressionModel_LUT_Collection_stem_plate_and_weighed_Hertz_Fprim.This model is a combination of the original ‘CompressionModel_LUT_Collection_stem_plate’ and Hertz. For situations where the capsule is parallel to the plane the contact force will be presented by a value in the LUT. For all non-parallel cases this force is augmented with a certain portion of the Hertz force, since no data is available for these cases and using only the original compression model would lead to forces which are too low, causing the capsules to go through the plane . The scaling factor (which determines the magnitude of the added ‘Hertz’ force) is calculated based on the angle between the plane normal and the capsule axis and will vary between 0 and 1. The damage is updated accordingly by multiplying the overlap with the same scaling factor. Using a scaling factor will guarantee a smooth transition. Following models ‘CompressionModel_LUT_Collection_stem_plate_Fprim’ and ‘CompressionModel_LUT_Collection_stem_plate’ will no longer work (are commented out) since the scaling factor is not provided in their chain. Geometry combinations available:

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

CompressionModel_LUT_Collection_stem_plate_and_weighed_Hertz_Fprim (Rigid_Quad Deformable_Capsule)

CompressionModel_LUT_Collection_stem_plate_and_weighed_Hertz_Fprim.This model is a combination of the original ‘CompressionModel_LUT_Collection_stem_plate’ and Hertz. For situations where the capsule is parallel to the plane the contact force will be presented by a value in the LUT. For all non-parallel cases this force is augmented with a certain portion of the Hertz force, since no data is available for these cases and using only the original compression model would lead to forces which are too low, causing the capsules to go through the plane . The scaling factor (which determines the magnitude of the added ‘Hertz’ force) is calculated based on the angle between the plane normal and the capsule axis and will vary between 0 and 1. The damage is updated accordingly by multiplying the overlap with the same scaling factor. Using a scaling factor will guarantee a smooth transition. Following models ‘CompressionModel_LUT_Collection_stem_plate_Fprim’ and ‘CompressionModel_LUT_Collection_stem_plate’ will no longer work (are commented out) since the scaling factor is not provided in their chain.

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.
    • c_plate — viscous damper constant.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • 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:
    • Fprim1 (default value = None) — Array with vectors which stores the contact force per primitive for pc1. If not given, the array pc1[‘Fprim’] is searched first, and if not found, nothing will be done
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • 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):

DefaultBoilerPlateDeformable_Cylinder_2_DataRigid_NGon_1_Data  <4 >
  ↓
NGon_CapsuleSphere12Hertz_CapsuleCapsuleTanAlphaForceWeighingN_LUT_Collection_stem_plateT_CundallStrackAssembleForces_Primitives_1AssembleForces_DeformableCylinder_2AssembleMoments_1AssembleForces_1Feedback

CompressionModel_LUT_Collection_stem_plate_and_weighed_Hertz_Fprim (Rigid_Quad Rigid_Capsule)

CompressionModel_LUT_Collection_stem_plate_and_weighed_Hertz_Fprim.This model is a combination of the original ‘CompressionModel_LUT_Collection_stem_plate’ and Hertz. For situations where the capsule is parallel to the plane the contact force will be presented by a value in the LUT. For all non-parallel cases this force is augmented with a certain portion of the Hertz force, since no data is available for these cases and using only the original compression model would lead to forces which are too low, causing the capsules to go through the plane . The scaling factor (which determines the magnitude of the added ‘Hertz’ force) is calculated based on the angle between the plane normal and the capsule axis and will vary between 0 and 1. The damage is updated accordingly by multiplying the overlap with the same scaling factor. Using a scaling factor will guarantee a smooth transition. Following models ‘CompressionModel_LUT_Collection_stem_plate_Fprim’ and ‘CompressionModel_LUT_Collection_stem_plate’ will no longer work (are commented out) since the scaling factor is not provided in their chain.

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.
    • c_plate — viscous damper constant.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • 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:
    • Fprim1 (default value = None) — Array with vectors which stores the contact force per primitive for pc1. If not given, the array pc1[‘Fprim’] is searched first, and if not found, nothing will be done
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • 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):

DefaultBoilerPlateRigid_Cylinder_2_DataRigid_NGon_1_Data  <4 >
  ↓
NGon_CapsuleSphere12Hertz_CapsuleCapsuleTanAlphaForceWeighingN_LUT_Collection_stem_plateT_CundallStrackAssembleForces_Primitives_1AssembleForcesAssembleMomentsFeedback

CompressionModel_LUT_Collection_stem_plate_and_weighed_Hertz_Fprim (Rigid_Triangle Deformable_Capsule)

CompressionModel_LUT_Collection_stem_plate_and_weighed_Hertz_Fprim.This model is a combination of the original ‘CompressionModel_LUT_Collection_stem_plate’ and Hertz. For situations where the capsule is parallel to the plane the contact force will be presented by a value in the LUT. For all non-parallel cases this force is augmented with a certain portion of the Hertz force, since no data is available for these cases and using only the original compression model would lead to forces which are too low, causing the capsules to go through the plane . The scaling factor (which determines the magnitude of the added ‘Hertz’ force) is calculated based on the angle between the plane normal and the capsule axis and will vary between 0 and 1. The damage is updated accordingly by multiplying the overlap with the same scaling factor. Using a scaling factor will guarantee a smooth transition. Following models ‘CompressionModel_LUT_Collection_stem_plate_Fprim’ and ‘CompressionModel_LUT_Collection_stem_plate’ will no longer work (are commented out) since the scaling factor is not provided in their chain.

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.
    • c_plate — viscous damper constant.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • 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:
    • Fprim1 (default value = None) — Array with vectors which stores the contact force per primitive for pc1. If not given, the array pc1[‘Fprim’] is searched first, and if not found, nothing will be done
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • 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):

DefaultBoilerPlateDeformable_Cylinder_2_DataRigid_NGon_1_Data  <3 >
  ↓
NGon_CapsuleSphere12Hertz_CapsuleCapsuleTanAlphaForceWeighingN_LUT_Collection_stem_plateT_CundallStrackAssembleForces_Primitives_1AssembleForces_DeformableCylinder_2AssembleMoments_1AssembleForces_1Feedback

CompressionModel_LUT_Collection_stem_plate_and_weighed_Hertz_Fprim (Rigid_Triangle Rigid_Capsule)

CompressionModel_LUT_Collection_stem_plate_and_weighed_Hertz_Fprim.This model is a combination of the original ‘CompressionModel_LUT_Collection_stem_plate’ and Hertz. For situations where the capsule is parallel to the plane the contact force will be presented by a value in the LUT. For all non-parallel cases this force is augmented with a certain portion of the Hertz force, since no data is available for these cases and using only the original compression model would lead to forces which are too low, causing the capsules to go through the plane . The scaling factor (which determines the magnitude of the added ‘Hertz’ force) is calculated based on the angle between the plane normal and the capsule axis and will vary between 0 and 1. The damage is updated accordingly by multiplying the overlap with the same scaling factor. Using a scaling factor will guarantee a smooth transition. Following models ‘CompressionModel_LUT_Collection_stem_plate_Fprim’ and ‘CompressionModel_LUT_Collection_stem_plate’ will no longer work (are commented out) since the scaling factor is not provided in their chain.

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.
    • c_plate — viscous damper constant.
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • 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:
    • Fprim1 (default value = None) — Array with vectors which stores the contact force per primitive for pc1. If not given, the array pc1[‘Fprim’] is searched first, and if not found, nothing will be done
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • 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):

DefaultBoilerPlateRigid_Cylinder_2_DataRigid_NGon_1_Data  <3 >
  ↓
NGon_CapsuleSphere12Hertz_CapsuleCapsuleTanAlphaForceWeighingN_LUT_Collection_stem_plateT_CundallStrackAssembleForces_Primitives_1AssembleForcesAssembleMomentsFeedback

CompressionModel_LUT_Collection_stem_stem

Description: Compression using collection of LUTs Geometry combinations available:

PC2 | PC1 - Rigid_Capsule Deformable_Capsule
Rigid_Capsule YES YES
Deformable_Capsule   YES

CompressionModel_LUT_Collection_stem_stem (Deformable_Capsule Deformable_Capsule)

Compression using collection of LUTs

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • k_tor — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • mu_tor — The static coulomb friction coefficient.
    • 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_stem (default value = 0) — viscous damper constant.
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • c_tor (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.

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_LUT_Collection_stem_stemT_CundallStrackTorsionalCundallStrackFrictionAssemble_Moment_To_Force_Couple2AssembleForces_DeformableCylinder_2AssembleForces_DeformableCylinder_1AbortIfNeighboursFeedback

CompressionModel_LUT_Collection_stem_stem (Deformable_Capsule Rigid_Capsule)

Compression using collection of LUTs

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • k_tor — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • mu_tor — The static coulomb friction coefficient.
    • 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_stem (default value = 0) — viscous damper constant.
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • c_tor (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.

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

DefaultBoilerPlateRigid_Cylinder_2_DataDeformable_Cylinder_1_DataCapsule_CapsuleN_LUT_Collection_stem_stemT_CundallStrackTorsionalCundallStrackFrictionAssemble_Moment_To_Force_Couple1AssembleMoments_2AssembleForces_2AssembleForces_DeformableCylinder_1AbortIfNeighboursFeedback

CompressionModel_LUT_Collection_stem_stem (Rigid_Capsule Rigid_Capsule)

Compression using collection of LUTs

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • k_tor — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • mu_tor — The static coulomb friction coefficient.
    • 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_stem (default value = 0) — viscous damper constant.
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • c_tor (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.

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

DefaultBoilerPlateRigid_Cylinder_2_DataRigid_Cylinder_1_DataCapsule_CapsuleN_LUT_Collection_stem_stemT_CundallStrackTorsionalCundallStrackFrictionAbortIfNeighboursFeedback

CompressionModel_LUT_Collection_stem_stemRecord

Description: Compression using collection of LUTs Geometry combinations available:

PC2 | PC1 - Rigid_Capsule Deformable_Capsule
Rigid_Capsule YES YES
Deformable_Capsule   YES

CompressionModel_LUT_Collection_stem_stemRecord (Deformable_Capsule Deformable_Capsule)

Compression using collection of LUTs

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • k_tor — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • mu_tor — The static coulomb friction coefficient.
    • 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:
    • E_friction (default value = 0) — The total dissipated friction energy by this contact model.
    • E_normal (default value = 0) — The total dissipated normal energy by this contact model.
    • c_stem (default value = 0) — viscous damper constant.
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • c_tor (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.

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_LUT_Collection_stem_stemT_CundallStrackTorsionalCundallStrackFrictionAssemble_Moment_To_Force_Couple2AssembleForces_DeformableCylinder_2AssembleForces_DeformableCylinder_1AbortIfNeighboursRecordDissipatedFrictionEnergyRecordDissipatedNormalEnergyFeedback

CompressionModel_LUT_Collection_stem_stemRecord (Deformable_Capsule Rigid_Capsule)

Compression using collection of LUTs

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • k_tor — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • mu_tor — The static coulomb friction coefficient.
    • 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:
    • E_friction (default value = 0) — The total dissipated friction energy by this contact model.
    • E_normal (default value = 0) — The total dissipated normal energy by this contact model.
    • c_stem (default value = 0) — viscous damper constant.
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • c_tor (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.

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

DefaultBoilerPlateRigid_Cylinder_2_DataDeformable_Cylinder_1_DataCapsule_CapsuleN_LUT_Collection_stem_stemT_CundallStrackTorsionalCundallStrackFrictionAssemble_Moment_To_Force_Couple1AssembleMoments_2AssembleForces_2AssembleForces_DeformableCylinder_1AbortIfNeighboursRecordDissipatedFrictionEnergyRecordDissipatedNormalEnergyFeedback

CompressionModel_LUT_Collection_stem_stemRecord (Rigid_Capsule Rigid_Capsule)

Compression using collection of LUTs

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • k_t — ‘k’ spring constant for static friction in tangential direction.
    • k_tor — ‘k’ spring constant for static friction in tangential direction.
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu_d — The dynamic coulomb friction coefficient.
    • mu_s — The static coulomb friction coefficient.
    • mu_tor — The static coulomb friction coefficient.
    • 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:
    • E_friction (default value = 0) — The total dissipated friction energy by this contact model.
    • E_normal (default value = 0) — The total dissipated normal energy by this contact model.
    • c_stem (default value = 0) — viscous damper constant.
    • c_t (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.
    • c_tor (default value = 0) — ‘c’ value of the linear dashpot in tangential direction.

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

DefaultBoilerPlateRigid_Cylinder_2_DataRigid_Cylinder_1_DataCapsule_CapsuleN_LUT_Collection_stem_stemT_CundallStrackTorsionalCundallStrackFrictionAbortIfNeighboursRecordDissipatedFrictionEnergyRecordDissipatedNormalEnergyFeedback

StemStemLUT_Coulomb_Record

Description: Geometry combinations available:

PC2 | PC1 - Rigid_Capsule Deformable_Capsule
Rigid_Capsule YES YES
Deformable_Capsule   YES

StemStemLUT_Coulomb_Record (Deformable_Capsule Deformable_Capsule)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A — Hertz damping coefficient.
    • E_end_compression — Young’s modulus. Used in the Hertz transition when the overlap betweens the stems becomes very big (see parameter start_hertz_r_).
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu — The coulomb friction coefficient (both static and dynamic).
    • nu — Poisson ratio. Used in the Hertz transition when the overlap betweens the stems becomes very big (see parameter start_hertz_r_).
    • 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.
    • start_end_compression_r — If the stems overlap more than this many times the summed radius (r1+r2), then the linear transition to a (very stiff) Hertz model starts.
  • Optional keywords:
    • E_friction (default value = 0) — The total dissipated friction energy by this contact model.
    • E_normal (default value = 0) — The total dissipated normal energy by this contact model.
    • c_stem (default value = 0) — viscous damper constant.
    • c_t (default value = -1) — ‘c’ value of the linear dashpot (N*s/m) in tangential direction. The higher this value, the more accurate the results will be, but the simulation can become unstable, requiring smaller timesteps.
    • cp_t (default value = -1) — Optional contact area-dependent linear dashpot coefficent (Pa*s/m). Give either c_t or cp_t but not both.

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_LUT_Collection_stem_stem_damped_hertz_end_compressionT_CoulombFrictionAssemble_Moment_To_Force_Couple2AssembleForces_DeformableCylinder_2AssembleForces_DeformableCylinder_1AbortIfNeighboursRecordDissipatedFrictionEnergyRecordDissipatedNormalEnergyFeedback

StemStemLUT_Coulomb_Record (Deformable_Capsule Rigid_Capsule)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A — Hertz damping coefficient.
    • E_end_compression — Young’s modulus. Used in the Hertz transition when the overlap betweens the stems becomes very big (see parameter start_hertz_r_).
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu — The coulomb friction coefficient (both static and dynamic).
    • nu — Poisson ratio. Used in the Hertz transition when the overlap betweens the stems becomes very big (see parameter start_hertz_r_).
    • 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.
    • start_end_compression_r — If the stems overlap more than this many times the summed radius (r1+r2), then the linear transition to a (very stiff) Hertz model starts.
  • Optional keywords:
    • E_friction (default value = 0) — The total dissipated friction energy by this contact model.
    • E_normal (default value = 0) — The total dissipated normal energy by this contact model.
    • c_stem (default value = 0) — viscous damper constant.
    • c_t (default value = -1) — ‘c’ value of the linear dashpot (N*s/m) in tangential direction. The higher this value, the more accurate the results will be, but the simulation can become unstable, requiring smaller timesteps.
    • cp_t (default value = -1) — Optional contact area-dependent linear dashpot coefficent (Pa*s/m). Give either c_t or cp_t but not both.

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

DefaultBoilerPlateRigid_Cylinder_2_DataDeformable_Cylinder_1_DataCapsule_CapsuleN_LUT_Collection_stem_stem_damped_hertz_end_compressionT_CoulombFrictionAssemble_Moment_To_Force_Couple1AssembleMoments_2AssembleForces_2AssembleForces_DeformableCylinder_1AbortIfNeighboursRecordDissipatedFrictionEnergyRecordDissipatedNormalEnergyFeedback

StemStemLUT_Coulomb_Record (Rigid_Capsule Rigid_Capsule)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A — Hertz damping coefficient.
    • E_end_compression — Young’s modulus. Used in the Hertz transition when the overlap betweens the stems becomes very big (see parameter start_hertz_r_).
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu — The coulomb friction coefficient (both static and dynamic).
    • nu — Poisson ratio. Used in the Hertz transition when the overlap betweens the stems becomes very big (see parameter start_hertz_r_).
    • 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.
    • start_end_compression_r — If the stems overlap more than this many times the summed radius (r1+r2), then the linear transition to a (very stiff) Hertz model starts.
  • Optional keywords:
    • E_friction (default value = 0) — The total dissipated friction energy by this contact model.
    • E_normal (default value = 0) — The total dissipated normal energy by this contact model.
    • c_stem (default value = 0) — viscous damper constant.
    • c_t (default value = -1) — ‘c’ value of the linear dashpot (N*s/m) in tangential direction. The higher this value, the more accurate the results will be, but the simulation can become unstable, requiring smaller timesteps.
    • cp_t (default value = -1) — Optional contact area-dependent linear dashpot coefficent (Pa*s/m). Give either c_t or cp_t but not both.

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

DefaultBoilerPlateRigid_Cylinder_2_DataRigid_Cylinder_1_DataCapsule_CapsuleN_LUT_Collection_stem_stem_damped_hertz_end_compressionT_CoulombFrictionAbortIfNeighboursRecordDissipatedFrictionEnergyRecordDissipatedNormalEnergyFeedback

StemStemLUT_RejectNonGradualOverlap_Coulomb_Record

Description: Geometry combinations available:

PC2 | PC1 - Rigid_Capsule Deformable_Capsule
Rigid_Capsule YES YES
Deformable_Capsule   YES

StemStemLUT_RejectNonGradualOverlap_Coulomb_Record (Deformable_Capsule Deformable_Capsule)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A — Hertz damping coefficient.
    • E_end_compression — Young’s modulus. Used in the Hertz transition when the overlap betweens the stems becomes very big (see parameter start_hertz_r_).
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu — The coulomb friction coefficient (both static and dynamic).
    • nu — Poisson ratio. Used in the Hertz transition when the overlap betweens the stems becomes very big (see parameter start_hertz_r_).
    • 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.
    • start_end_compression_r — If the stems overlap more than this many times the summed radius (r1+r2), then the linear transition to a (very stiff) Hertz model starts.
    • v_max (m . s^-1) — Maximal relative velocity between two colliding particles.
  • Optional keywords:
    • E_friction (default value = 0) — The total dissipated friction energy by this contact model.
    • E_normal (default value = 0) — The total dissipated normal energy by this contact model.
    • c_stem (default value = 0) — viscous damper constant.
    • c_t (default value = -1) — ‘c’ value of the linear dashpot (N*s/m) in tangential direction. The higher this value, the more accurate the results will be, but the simulation can become unstable, requiring smaller timesteps.
    • cp_t (default value = -1) — Optional contact area-dependent linear dashpot coefficent (Pa*s/m). Give either c_t or cp_t but not both.

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_CapsuleRejectNonGradualOverlapN_LUT_Collection_stem_stem_damped_hertz_end_compressionT_CoulombFrictionAssemble_Moment_To_Force_Couple2AssembleForces_DeformableCylinder_2AssembleForces_DeformableCylinder_1AbortIfNeighboursRecordDissipatedFrictionEnergyRecordDissipatedNormalEnergyFeedback

StemStemLUT_RejectNonGradualOverlap_Coulomb_Record (Deformable_Capsule Rigid_Capsule)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A — Hertz damping coefficient.
    • E_end_compression — Young’s modulus. Used in the Hertz transition when the overlap betweens the stems becomes very big (see parameter start_hertz_r_).
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu — The coulomb friction coefficient (both static and dynamic).
    • nu — Poisson ratio. Used in the Hertz transition when the overlap betweens the stems becomes very big (see parameter start_hertz_r_).
    • 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.
    • start_end_compression_r — If the stems overlap more than this many times the summed radius (r1+r2), then the linear transition to a (very stiff) Hertz model starts.
    • v_max (m . s^-1) — Maximal relative velocity between two colliding particles.
  • Optional keywords:
    • E_friction (default value = 0) — The total dissipated friction energy by this contact model.
    • E_normal (default value = 0) — The total dissipated normal energy by this contact model.
    • c_stem (default value = 0) — viscous damper constant.
    • c_t (default value = -1) — ‘c’ value of the linear dashpot (N*s/m) in tangential direction. The higher this value, the more accurate the results will be, but the simulation can become unstable, requiring smaller timesteps.
    • cp_t (default value = -1) — Optional contact area-dependent linear dashpot coefficent (Pa*s/m). Give either c_t or cp_t but not both.

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

DefaultBoilerPlateRigid_Cylinder_2_DataDeformable_Cylinder_1_DataCapsule_CapsuleRejectNonGradualOverlapN_LUT_Collection_stem_stem_damped_hertz_end_compressionT_CoulombFrictionAssemble_Moment_To_Force_Couple1AssembleMoments_2AssembleForces_2AssembleForces_DeformableCylinder_1AbortIfNeighboursRecordDissipatedFrictionEnergyRecordDissipatedNormalEnergyFeedback

StemStemLUT_RejectNonGradualOverlap_Coulomb_Record (Rigid_Capsule Rigid_Capsule)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • A — Hertz damping coefficient.
    • E_end_compression — Young’s modulus. Used in the Hertz transition when the overlap betweens the stems becomes very big (see parameter start_hertz_r_).
    • luts_down — List of LUT’s (down)
    • luts_up — List of LUT’s (up)
    • mu — The coulomb friction coefficient (both static and dynamic).
    • nu — Poisson ratio. Used in the Hertz transition when the overlap betweens the stems becomes very big (see parameter start_hertz_r_).
    • 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.
    • start_end_compression_r — If the stems overlap more than this many times the summed radius (r1+r2), then the linear transition to a (very stiff) Hertz model starts.
    • v_max (m . s^-1) — Maximal relative velocity between two colliding particles.
  • Optional keywords:
    • E_friction (default value = 0) — The total dissipated friction energy by this contact model.
    • E_normal (default value = 0) — The total dissipated normal energy by this contact model.
    • c_stem (default value = 0) — viscous damper constant.
    • c_t (default value = -1) — ‘c’ value of the linear dashpot (N*s/m) in tangential direction. The higher this value, the more accurate the results will be, but the simulation can become unstable, requiring smaller timesteps.
    • cp_t (default value = -1) — Optional contact area-dependent linear dashpot coefficent (Pa*s/m). Give either c_t or cp_t but not both.

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

DefaultBoilerPlateRigid_Cylinder_2_DataRigid_Cylinder_1_DataCapsule_CapsuleRejectNonGradualOverlapN_LUT_Collection_stem_stem_damped_hertz_end_compressionT_CoulombFrictionAbortIfNeighboursRecordDissipatedFrictionEnergyRecordDissipatedNormalEnergyFeedback