Hysteretic Model

Introduction

The hysteretic contact model in PFC consists of a combination of the elastic portion of the Hertz model as described in the “Hertz Model” section, combined with an alternate dashpot group consisting of a nonlinear visco-elastic element in the normal direction. This model is referred to in commands and FISH as hysteretic.

Note

Revision 28 included the following additional modification:

The property inheritance mechanism that derives the contact model elastic parameters from the parameters of the contacting pieces has been changed to conform to [Mindlin1953], following the equations below . Before this modifications, the contact effective shear modulus and Poisson ratio were taken as the average of the contacting pieces shear modulus and Poisson ratio, respectively. Note that this change may change the output for models that used the property inheritance mechanism with the hysteretic model. Please refer to the documentation of the “Hertz Model” for further information.

Behavior Summary

The formulation is based on the literature review by [Machado2012]. Three different dashpot modes, which correspond to three different expressions for the dashpot force calculation, are provided.

Activity-Deletion Criteria

The hysteretic model defines an activity distance of 0.0. Contact activity is updated as:

(1){gc0.0contact activegc>0.0contact inactive

As discussed in the “Contact Resolution” section, force-displacement calculations are skipped for inactive contacts.

Force-Displacement Law

Similar to the Hertz model, the hysteretic model transmits a force and no moment. The force embeds a non-linear elastic component Fh and a viscous component Fd

(2)Fc=Fh+FdMc0

The nonlinear Hertz force can be decomposed into a normal and a shear contribution:

(3)Fh=Fhnˆnc+Fhs

Its formulation is equivalent to that of the non-linear Hertz force described in the Hertz contact model.

The dashpot force has a normal component only:

(4)Fd=Fdnˆnc

where the formulation of Fdn can be expressed in terms of a hysteresis damping factor χ as:

(5)Fdn=χδαn˙δn

where δn=gc represents physical overlap, ˙δn=˙δˆnc is the normal component of the relative translational velocity, and the hysteresis damping factor χ depends on the dashpot mode Md, the target restitution coefficient en, the Hertz normal force coefficient hn (see here), and the normal impact velocity ˙δ()n with:

(6)χ={1e2nenhn˙δ()nMd=03(1e2n)4hn˙δ()nMd=13(1e2n)e2(1en)4hn˙δ()nMd=2

Note that for the case Md=1, the target restitution coefficient used in the equation above is the maximal value between en and a minimal value of en.

Energy Partitions

The hysteretic model defines three energy partitions:

  • strain energy Ek, stored in the springs;
  • slip work Eμ, defined as total energy dissipated by frictional sliding; and
  • damp work Eβ, defined as the total energy dissipated by the dashpots.

If energy tracking is activated, these energy partitions are updated as follows.

  1. Update the strain energy:

    (7)Ek=α(α+1)((Fhn)2kn)+12Fhs)2ks
  2. Update the slip energy:

    (8)Eμ:=Eμ12((Fhs)o+Fhs)ΔδδμswithΔδδμs=ΔδδsΔδδks=Δδδs(Fhs(Fhs)oks)

    where (Fhs)o is the shear force at the beginning of the timestep, and the relative shear-displacement increment of this equation of the “Contact Resolution” section has been decomposed into an elastic (Δδδks) and a slip (Δδδμs) component. The dot product of the slip component and the average shear force occurring during the timestep gives the increment of slip energy.

  3. Update the dashpot energy:

    (9)Eβ:=EβFd(˙δδΔt)

    where ˙δδ is the relative translational velocity of this equation of the “Contact Resolution” section.

Additional information, including the keywords by which these partitions are referred to in commands and FISH, is provided in the table below.

Table 1: Hysteretic Model Energies
Keyword Symbol Description Range Accumulated
Hertz Group:
energy‑strain Ek strain energy [0.0,+) NO
energy‑slip Eμ total energy dissipated by slip [0.0,+) YES
Dashpot Group:
energy‑dashpot Eβ total energy dissipated by dashpots [0.0,+) YES

Properties

The properties defined by the hysteretic contact model are listed in the table below for a concise reference; see the “Contact Properties” section for a description of the information in the table columns. The mapping from the surface inheritable properties to the contact model properties is also discussed below.

Table 2: Hysteretic model properties
Keyword Symbol Description Type Range Default Modifiable Inheritable
hysteretic Model name
Hertz Group:
hz_shear G Shear modulus [stress] FLT [0.0,+) 0.0 YES YES
hz_poiss ν Poisson’s ratio [-] FLT (1.0,0.5] 0.0 YES YES
fric μ Friction coefficient [-] FLT [0.0,+) 0.0 YES YES
hz_mode Ms Shear-force scaling mode [-] INT {0;1} 0 YES NO
    {0: no scaling1: scaling is active          
hz_alpha α Exponent [_] FLT [0.0,+) 1.5 YES NO
hz_slip s Slip state [-] BOOL {false,true} false NO N/A
hz_force Fh Hertz force (contact plane coord. system) VEC R3 0 YES NO
    (Fhn,Fhss,Fhst)(2D model: Fhss0)          
Dashpot Group:
dp_mode Md Dashpot mode [-] INT {0;1;2} 0 YES NO
    {0: default mode1: mode 12: mode 2          
dp_en en Target restitution coef. [-] FLT [0.0,1.0] 1.0 YES NO
dp_enmin en Minimal restitution coef. [-] FLT [0.0,1.0] 0.0 YES NO
dp_force Fd Dashpot force (contact plane coord. system) VEC R3 0 NO NO
    (Fdn,0,0)          

Property Inheritance

The shear modulus G, Poisson’s ratio ν, and friction coefficient μ may be inherited from the contacting pieces. The inheritance logic is similar to that of the hertz contact model described here.

Methods

No methods are defined by the hysteretic model.

Callback events

Table 3: Hysteretic Contact Model Callback Events
Event Array slot Value Type Range Description
contact_activated Contact becomes active
  1 C_PNT N/A contact pointer
slip_change Slip state has changed
  1 C_PNT N/A contact pointer
  2 INT {0;1} slip change mode
        {0: slip has initiated1: slip has ended


Model Summary

An alphabetical list of the hysteretic model properties is given here.

This is a hidden link to: [Gonthier2004]

And this is a hidden link to: [Lankarani1990]

References

[Gonthier2004] Y. Gonthier, J. McPhee, C. Lange, J-C. Piedboeuf, “A regularized contact model with asymmetric damping and dwell-time dependent friction”, Multibody System Dynamics, vol. 11, pp. 209-233 (2004).
[Lankarani1990]H.M. Lankarani, P.E. Nikravesh, “A contact force model with hysteresis damping for impact analysis of multibody systems”, Journal of Mechanical Design, vol. 112, pp. 369-376 (1990).
[Machado2012] M. Machado, P. Moreira, P. Flores, H.M. Lankarani, “Compliant contact force models in multibody dynamics: Evolution of the Hertz contact theory”, Mechanism and Machine Theory, vol. 53, pp. 99-121 (July 2012).