*SetContact() - Poisson
Sets the values associated with a contact entity where the normal forces are described by the Poisson method.
Syntax
*SetContact(contact_name,
penalty,
rest_coeff,
use_aug_formulation,
[mu_static],
[mu_dynamic],
[stiction_trans_vel],
[friction_trans_vel])
Arguments
- contact_name
- The variable name of an existing contact entity.
- penalty
- Determines the local stiffness properties between materials. Larger values lead to reduced penetration between two bodies.
- rest_coeff
- The value representing the energy loss between the two bodies in contact. A value of one represents no energy loss and a perfectly elastic contact.
- use_aug_formulation
- Used to refine the accuracy of the normal force between the two bodies.
- mu_static
- The coefficient of friction when the slip velocity is less than stiction_trans_vel.
- mu_dynamic
- The coefficient of friction when the slip velocity is greater than friction_trans_vel.
- stiction_trans_vel
- The value where the coefficient of friction becomes mu_static. When the slip velocity is between stiction_trans_vel and friction_trans_vel, the coefficient of friction is in transition between the two.
- friction_trans_vel
- The value where the coefficient of friction becomes mu_dynamic. When the slip velocity is between stiction_trans_vel and friction_trans_vel, the coefficient of friction is in transition between the two.
Example
*BeginMDL(the_model, "MBD Model")
*Body(b_1, "Body 1", , , )
*Body(b_2, "Body 2", , , )
*Graphic(g_file_1, "File Graphic - 1", File, MODEL.b_1, "file1.h3d", comp1, GLOBAL, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 1.000, 1.000, 1.000)
*Graphic(g_file_2, "File Graphic - 2", File, MODEL.b_2, "file2.h3d", comp1, GLOBAL, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 1.000, 1.000, 1.000)
*Contact(contact_im, "Contact Impact", IMPACT, OFF, 1, g_file_1, FALSE, 1, g_file_2, FALSE)
*SetContact(contact_im, 1E6, 2.0, 8.0, 0.1)
*Contact(contact_ps, "Contact Poisson", POISSON, ON, 1, g_file_1, FALSE, 1, g_file_2, FALSE)
*SetContact(contact_ps, 1E5, 0.8, FALSE, 0.8, 0.6, 0.5, 1.5)
*EndMDL()