mpacts.contact.models.collision.cellmigration. cellmigration_rt_matrix

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

import mpacts.contact.models.collision.cellmigration.cellmigration_rt_matrix
#or assign it to a shorter name
import mpacts.contact.models.collision.cellmigration.cellmigration_rt_matrix as cel

MDMotionIntMatrix

Description: Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle. Geometry combinations available:

PC2 | PC1 - Rigid_RoundedTriangle Deformable_RoundedTriangle
Sphere YES YES
Rigid_Triangle YES YES
Rigid_RoundedTriangle YES YES
Rigid_Quad YES YES
Rigid_CylinderTop YES YES
Rigid_CylinderBottom YES YES
Rigid_Cylinder YES YES
Deformable_Triangle   YES
Deformable_RoundedTriangle   YES

MDMotionIntMatrix (Deformable_RoundedTriangle Deformable_RoundedTriangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • contact_area2 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • 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.
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataDeformable_NGon_1_Data  <3 >
  ↓
Deformable_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1DeformableCellMotionForce_2N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_DeformableTriangle_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1Assemble_VirialStress_DeformableTriangle_2ComputeContactMatrixElementOverdampedDeformableTriangleDof_1DeformableTriangleDof_2AssembleContactMatrix_ComposedParticleAbortIfSameParentFeedback

MDMotionIntMatrix (Deformable_RoundedTriangle Deformable_Triangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Deformable_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_DeformableTriangle_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1DeformableTriangleDof_2AssembleContactMatrix_ComposedParticleFeedback

MDMotionIntMatrix (Deformable_RoundedTriangle Rigid_Cylinder)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_CylinderNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionIntMatrix (Deformable_RoundedTriangle Rigid_CylinderBottom)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Bottom_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionIntMatrix (Deformable_RoundedTriangle Rigid_CylinderTop)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Top_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionIntMatrix (Deformable_RoundedTriangle Rigid_Quad)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <4 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionIntMatrix (Deformable_RoundedTriangle Rigid_RoundedTriangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionIntMatrix (Deformable_RoundedTriangle Rigid_Triangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionIntMatrix (Deformable_RoundedTriangle Sphere)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Sphere0_2_DataRoundedTriangle_SphereNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1SphereDof_2AssembleContactMatrix_ComposedParticleFeedback

MDMotionIntMatrix (Rigid_RoundedTriangle Rigid_Cylinder)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_CylinderNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionIntMatrix (Rigid_RoundedTriangle Rigid_CylinderBottom)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Bottom_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionIntMatrix (Rigid_RoundedTriangle Rigid_CylinderTop)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Top_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionIntMatrix (Rigid_RoundedTriangle Rigid_Quad)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <4 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionIntMatrix (Rigid_RoundedTriangle Rigid_RoundedTriangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • 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.
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2AbortIfSameParentFeedback

MDMotionIntMatrix (Rigid_RoundedTriangle Rigid_Triangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionIntMatrix (Rigid_RoundedTriangle Sphere)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Sphere_2_DataRoundedTriangle_SphereNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2ComputeContactMatrixElementOverdampedRigidBodyDof_1SphereDof_2AssembleContactMatrixToDiagonal_2Feedback

MDMotionLocFrictionIntMatrix

Description: Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle. Geometry combinations available:

PC2 | PC1 - Rigid_RoundedTriangle Deformable_RoundedTriangle
Sphere YES YES
Rigid_Triangle YES YES
Rigid_RoundedTriangle YES YES
Rigid_Quad YES YES
Rigid_CylinderTop YES YES
Rigid_CylinderBottom YES YES
Rigid_Cylinder YES YES
Deformable_Triangle   YES
Deformable_RoundedTriangle   YES

MDMotionLocFrictionIntMatrix (Deformable_RoundedTriangle Deformable_RoundedTriangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • contact_area2 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • 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.
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataDeformable_NGon_1_Data  <3 >
  ↓
Deformable_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1DeformableCellMotionForce_2ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_DeformableTriangle_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2ComputeContactMatrixElementOverdampedDeformableTriangleDof_1DeformableTriangleDof_2AssembleContactMatrix_ComposedParticleAbortIfSameParentFeedback

MDMotionLocFrictionIntMatrix (Deformable_RoundedTriangle Deformable_Triangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Deformable_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_DeformableTriangle_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2ComputeContactMatrixElementOverdampedDeformableTriangleDof_1DeformableTriangleDof_2AssembleContactMatrix_ComposedParticleFeedback

MDMotionLocFrictionIntMatrix (Deformable_RoundedTriangle Rigid_Cylinder)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_CylinderNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionLocFrictionIntMatrix (Deformable_RoundedTriangle Rigid_CylinderBottom)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Bottom_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionLocFrictionIntMatrix (Deformable_RoundedTriangle Rigid_CylinderTop)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Top_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionLocFrictionIntMatrix (Deformable_RoundedTriangle Rigid_Quad)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <4 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionLocFrictionIntMatrix (Deformable_RoundedTriangle Rigid_RoundedTriangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionLocFrictionIntMatrix (Deformable_RoundedTriangle Rigid_Triangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionLocFrictionIntMatrix (Deformable_RoundedTriangle Sphere)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Sphere0_2_DataRoundedTriangle_SphereNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2ComputeContactMatrixElementOverdampedDeformableTriangleDof_1SphereDof_2AssembleContactMatrix_ComposedParticleFeedback

MDMotionLocFrictionIntMatrix (Rigid_RoundedTriangle Rigid_Cylinder)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_CylinderNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionLocFrictionIntMatrix (Rigid_RoundedTriangle Rigid_CylinderBottom)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Bottom_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionLocFrictionIntMatrix (Rigid_RoundedTriangle Rigid_CylinderTop)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Top_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionLocFrictionIntMatrix (Rigid_RoundedTriangle Rigid_Quad)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <4 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionLocFrictionIntMatrix (Rigid_RoundedTriangle Rigid_RoundedTriangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • 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.
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2AbortIfSameParentFeedback

MDMotionLocFrictionIntMatrix (Rigid_RoundedTriangle Rigid_Triangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionLocFrictionIntMatrix (Rigid_RoundedTriangle Sphere)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal1 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_normal2 — Normal friction coefficient for the material of particle 1 (array or single value).
    • gamma_tangential1 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • gamma_tangential2 — Tangential friction coefficient for the material of particle 2 (array or single value).
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • gamma_normal — Normal friction coefficient (kg/s)
    • gamma_tangential — Tangential friction coefficient (kg/s)
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Sphere_2_DataRoundedTriangle_SphereNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
ArrayValue_Property  <double >
  ↓
N_T_MaterialSpecificAreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2ComputeContactMatrixElementOverdampedRigidBodyDof_1SphereDof_2AssembleContactMatrixToDiagonal_2Feedback

MDMotionReinforcementIntMatrix

Description: Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle. Geometry combinations available:

PC2 | PC1 - Rigid_RoundedTriangle Deformable_RoundedTriangle
Sphere YES YES
Rigid_Triangle YES YES
Rigid_RoundedTriangle YES YES
Rigid_Quad YES YES
Rigid_CylinderTop YES YES
Rigid_CylinderBottom YES YES
Rigid_Cylinder YES YES
Deformable_Triangle   YES
Deformable_RoundedTriangle   YES

MDMotionReinforcementIntMatrix (Deformable_RoundedTriangle Deformable_RoundedTriangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • contact_area2 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • 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.
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataDeformable_NGon_1_Data  <3 >
  ↓
Deformable_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1DeformableCellMotionForce_2AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_DeformableTriangle_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1Assemble_VirialStress_DeformableTriangle_2ComputeContactMatrixElementOverdampedDeformableTriangleDof_1DeformableTriangleDof_2AssembleContactMatrix_ComposedParticleAbortIfSameParentFeedback

MDMotionReinforcementIntMatrix (Deformable_RoundedTriangle Deformable_Triangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Deformable_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_DeformableTriangle_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1DeformableTriangleDof_2AssembleContactMatrix_ComposedParticleFeedback

MDMotionReinforcementIntMatrix (Deformable_RoundedTriangle Rigid_Cylinder)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_CylinderSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionReinforcementIntMatrix (Deformable_RoundedTriangle Rigid_CylinderBottom)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Bottom_Selector >
  ↓
SingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionReinforcementIntMatrix (Deformable_RoundedTriangle Rigid_CylinderTop)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Top_Selector >
  ↓
SingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionReinforcementIntMatrix (Deformable_RoundedTriangle Rigid_Quad)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <4 >
  ↓
RoundedTriangle_NGonSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionReinforcementIntMatrix (Deformable_RoundedTriangle Rigid_RoundedTriangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionReinforcementIntMatrix (Deformable_RoundedTriangle Rigid_Triangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDMotionReinforcementIntMatrix (Deformable_RoundedTriangle Sphere)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • contact_area1 — ‘contact_area’ array for current contact loop.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Sphere0_2_DataRoundedTriangle_SphereSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
DeformableCellMotionForce_1AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1SphereDof_2AssembleContactMatrix_ComposedParticleFeedback

MDMotionReinforcementIntMatrix (Rigid_RoundedTriangle Rigid_Cylinder)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_CylinderSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionReinforcementIntMatrix (Rigid_RoundedTriangle Rigid_CylinderBottom)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Bottom_Selector >
  ↓
SingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionReinforcementIntMatrix (Rigid_RoundedTriangle Rigid_CylinderTop)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Top_Selector >
  ↓
SingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionReinforcementIntMatrix (Rigid_RoundedTriangle Rigid_Quad)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <4 >
  ↓
RoundedTriangle_NGonSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionReinforcementIntMatrix (Rigid_RoundedTriangle Rigid_RoundedTriangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • 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.
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2AbortIfSameParentFeedback

MDMotionReinforcementIntMatrix (Rigid_RoundedTriangle Rigid_Triangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDMotionReinforcementIntMatrix (Rigid_RoundedTriangle Sphere)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy + walkforce’. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle, contact areas as a scalar per triangle and calculates a pressure tensor for the parent particle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • attrFactorMax — maximum factor for contact state adhesion modifier
    • attrFactorTime — time it takes to linearly increase from attrFactorMin to attrFactorMax
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • max_angle — Maximal angle between the two planes for which the adhesion can be reinforced
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • attrFactorMin (default value = 1) — minimum factor for contact state adhesion modifier
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateContactMatrixContactStateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Sphere_2_DataRoundedTriangle_SphereSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
AdhesionReinforcementN_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2ComputeContactMatrixElementOverdampedRigidBodyDof_1SphereDof_2AssembleContactMatrixToDiagonal_2Feedback

MDTractionIntMatrix

Description: Explanation coming... ;-) Geometry combinations available:

PC2 | PC1 - Rigid_RoundedTriangle Deformable_RoundedTriangle
Sphere YES YES
Rigid_Triangle YES YES
Rigid_RoundedTriangle YES YES
Rigid_Quad YES YES
Rigid_CylinderTop YES YES
Rigid_CylinderBottom YES YES
Rigid_Cylinder YES YES
Deformable_Triangle   YES
Deformable_RoundedTriangle   YES

MDTractionIntMatrix (Deformable_RoundedTriangle Deformable_RoundedTriangle)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • 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.
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataDeformable_NGon_1_Data  <3 >
  ↓
Deformable_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1T_ActiveTraction_Front_2N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_DeformableTriangle_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1Assemble_VirialStress_DeformableTriangle_2ComputeContactMatrixElementOverdampedDeformableTriangleDof_1DeformableTriangleDof_2AssembleContactMatrix_ComposedParticleAbortIfSameParentFeedback

MDTractionIntMatrix (Deformable_RoundedTriangle Deformable_Triangle)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Deformable_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_DeformableTriangle_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1DeformableTriangleDof_2AssembleContactMatrix_ComposedParticleFeedback

MDTractionIntMatrix (Deformable_RoundedTriangle Rigid_Cylinder)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_CylinderNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDTractionIntMatrix (Deformable_RoundedTriangle Rigid_CylinderBottom)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Bottom_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDTractionIntMatrix (Deformable_RoundedTriangle Rigid_CylinderTop)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Top_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDTractionIntMatrix (Deformable_RoundedTriangle Rigid_Quad)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <4 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDTractionIntMatrix (Deformable_RoundedTriangle Rigid_RoundedTriangle)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDTractionIntMatrix (Deformable_RoundedTriangle Rigid_Triangle)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDTractionIntMatrix (Deformable_RoundedTriangle Sphere)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Sphere0_2_DataRoundedTriangle_SphereNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1SphereDof_2AssembleContactMatrix_ComposedParticleFeedback

MDTractionIntMatrix (Rigid_RoundedTriangle Rigid_Cylinder)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_CylinderNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDTractionIntMatrix (Rigid_RoundedTriangle Rigid_CylinderBottom)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Bottom_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDTractionIntMatrix (Rigid_RoundedTriangle Rigid_CylinderTop)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Top_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDTractionIntMatrix (Rigid_RoundedTriangle Rigid_Quad)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <4 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDTractionIntMatrix (Rigid_RoundedTriangle Rigid_RoundedTriangle)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • 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.
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2AbortIfSameParentFeedback

MDTractionIntMatrix (Rigid_RoundedTriangle Rigid_Triangle)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • layer_width (default value = 0) — flat layer width
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Feedback

MDTractionIntMatrix (Rigid_RoundedTriangle Sphere)

Explanation coming... ;-)

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRigid_NGon_1_Data  <3 >
  ↓
Sphere_2_DataRoundedTriangle_SphereNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_1AssembleMoments_cps_1AssembleMoments_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2ComputeContactMatrixElementOverdampedRigidBodyDof_1SphereDof_2AssembleContactMatrixToDiagonal_2Feedback

MDTractionSubstrateIntMatrix

Description: Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy and lamellopodim tractions. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle. Geometry combinations available:

PC2 | PC1 - Rigid_RoundedTriangle Deformable_RoundedTriangle
Sphere YES YES
Rigid_Triangle YES YES
Rigid_RoundedTriangle YES YES
Rigid_Quad YES YES
Rigid_CylinderTop YES YES
Rigid_CylinderBottom YES YES
Rigid_Cylinder YES YES
Deformable_Triangle   YES
Deformable_RoundedTriangle   YES

MDTractionSubstrateIntMatrix (Deformable_RoundedTriangle Deformable_RoundedTriangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy and lamellopodim tractions. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • 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.
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataDeformable_NGon_1_Data  <3 >
  ↓
Deformable_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1T_ActiveTraction_Front_2N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_DeformableTriangle_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1Assemble_VirialStress_DeformableTriangle_2ComputeContactMatrixElementOverdampedDeformableTriangleDof_1DeformableTriangleDof_2AssembleContactMatrix_ComposedParticleAbortIfSameParentFeedback

MDTractionSubstrateIntMatrix (Deformable_RoundedTriangle Deformable_Triangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy and lamellopodim tractions. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (default value = None) — The contact data belonging to the DoF (parents) made by this model.
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
    • transpose (default value = 0) — switches around pc1 and pc2 for the ContactMatrix.
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Deformable_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_DeformableTriangle_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1DeformableTriangleDof_2AssembleContactMatrix_ComposedParticleFeedback

MDTractionSubstrateIntMatrix (Deformable_RoundedTriangle Rigid_Cylinder)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy and lamellopodim tractions. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_CylinderNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDTractionSubstrateIntMatrix (Deformable_RoundedTriangle Rigid_CylinderBottom)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy and lamellopodim tractions. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Bottom_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDTractionSubstrateIntMatrix (Deformable_RoundedTriangle Rigid_CylinderTop)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy and lamellopodim tractions. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_Cylinder_2_DataRoundedTriangle_Disk  <Top_Selector >
  ↓
NoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDTractionSubstrateIntMatrix (Deformable_RoundedTriangle Rigid_Quad)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy and lamellopodim tractions. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <4 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDTractionSubstrateIntMatrix (Deformable_RoundedTriangle Rigid_RoundedTriangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy and lamellopodim tractions. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataRoundedTriangle_2_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_RoundedTriangleNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesDistributeForcesAndMomentsToTriangle_2_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2AssembleTriangleForcesToNodes_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDTractionSubstrateIntMatrix (Deformable_RoundedTriangle Rigid_Triangle)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy and lamellopodim tractions. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • implicitness (default value = 0) — Choose the semi-implicit integration method. (0 for explicit Euler, 1 for implicit Euler, 0.5 for Cranck-Nicholson).
    • layer_width (default value = 0) — flat layer width
    • modify_traction (default value = 1) — factor modifying the triangle’s fraction
    • never_integrate (default value = 0) — If true, the contact force will never be obtained through integration, but will just be directly computed
  • Read only properties:
    • sigma0 — Maximum adhesive tension at the edge of the contact circle. NOTE: sigma0 = attrConst / effective_range

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

RoundedTriangleBoilerPlateContactMatrixDataBoilerPlateKeepTimeStepRoundedTriangle_1_DataDeformable_NGon_1_Data  <3 >
  ↓
Rigid_NGon_2_Data  <3 >
  ↓
RoundedTriangle_NGonNoContactAdhesionModifierSingleValue_Property  <double >
  ↓
N_MaugisDugdale_geo_Int  <16 >
  ↓
T_ActiveTraction_Front_1N_T_AreaWeightedResistanceDistributeForcesAndMomentsToTriangle_1_NodesAssembleForces_DeformableTriangle_1AssembleForces_2AssembleMoments_cps_2AssembleMoments_2Assemble_ContactArea_1Assemble_ContactArea_2AssembleForces_Primitives_1AssembleForces_Primitives_2DistributeForcesAndMomentsToTriangle_2_NodesAssembleTriangleForcesToNodes_2Assemble_VirialStress_DeformableTriangle_1ComputeContactMatrixElementOverdampedDeformableTriangleDof_1RigidBodyDof_2AssembleContactMatrixToDiagonal_ComposedParticle_1FrictionMatrixForce_1Feedback

MDTractionSubstrateIntMatrix (Deformable_RoundedTriangle Sphere)

Maugis-Dugdale adhesive contact with effective adhesive range and adhesion energy and lamellopodim tractions. Geometric contact radius is used as inner radius a, so only c is calculated (and not both a and c). Saves contact forces as a vector per triangle and contact areas as a scalar per triangle.

Parallel Compatible: Yes

Properties:

  • Required keywords:
    • E1 — Young modulus material 1.
    • E2 — Young modulus material 2.
    • effective_range — effective range of adhesive interaction. Relates adhesive pressure to adhesive energy MD.
    • gamma_normal (kg . s^-1 . m^-2) — Normal friction coefficient (Pa*s/m)
    • gamma_tangential (kg . s^-1 . m^-2) — Tangential friction coefficient (Pa*s/m)
    • 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
    • Fprim2 (default value = None) — Array with vectors which stores the contact force per primitive for pc2. If not given, the array pc2[‘Fprim’] is searched first, and if not found, nothing will be done
    • _contactmatrixdatatype (default value = 1) — Identifier that signals that this chain element’s containing contact model is a friction matrix type. Only change this if you really know what you are doing!
    • attrConst1 (default value = 0) — Adhesion energy of material 1
    • attrConst2 (default value = 0) — Adhesion energy of material 2
    • attrConst (default value = -1) — effective adhesion energy [J/m^2]. Overrules attrConst1 and attrConst2 if given.
    • c_geo (default value = 0) — If ‘True’, also the adhesive radius ‘c’ will be derived from geometry only, and will NOT take into account the implicit elastic deformation due to adhesion. Set this to ‘True’ if your deformable surface is sufficiently ‘compliant’ with respect to its internal Young’s modulus (i.e. more valid for more thin but more stiff shell-like structures, and less valid for very thick soft structures. The default setting is False.
    • contact_area_c_factor (default value = 1) — Sets the relative contribution of any area between Maugis-Dugdale radii ‘a’ and ‘c’ to the contact area. Default = 1, so the contact area is the total adhesive area
    • contact_data_storage_factory (default value = None) — The contact data storage factory to make contact data. Do not change! Especially vector in parallel is not a safe choice here.
    • contact_data (de