# Burger’s Model

Burger’s model simulates creep mechanisms by using a Kelvin model and a Maxwell model connected in series in both the normal and shear directions.

Introduction

The Burger’s model provides a Kelvin model acting in series with a Maxwell model, in both the normal and shear directions. The Kelvin model is the combination of linear spring and dashpot components that act in parallel. The Maxwell model, on the other hand, is the combination of linear spring and dashpot components that act in series. Burger’s model acts over a vanishingly small area, and thus transmits only a force.

## Behavior Summary

The rheological components of the Burger’s model are shown in Figure 1 for both the normal and shear directions. In both directions, the model combines Kelvin and Maxwell models acting in series. In the normal direction, the Kelvin model provides a linear spring with stiffness \(K_{k_n}\) and a dashpot with viscosity \(C_{k_n}\), and the Maxwell model provides a linear spring with stiffness \(K_{m_n}\) and a dashpot with viscosity \(C_{m_n}\). The Burger’s model can sustain tensile forces (\(M_t=0\)) or not (\(M_t=1\)). In the shear direction, , the Kelvin model provides a linear spring with stiffness \(K_{k_s}\) and a dashpot with viscosity \(C_{k_s}\), and the Maxwell model provides a linear spring with stiffness \(K_{m_s}\) and a dashpot with viscosity \(C_{m_s}\). A slider with friction coefficient \(f_s\) limits the value of the shear force according to a Coulomb law.

In each direction, the total displacement of the Burger’s model, \(u\), is the sum of the displacement of the Kelvin section (\(u_k\)) and that of the Maxwell section ( \(u_{m_K}\) , \(u_{m_C}\) ) of the model:

The first and second derivatives of the equation above are given by :

The contact forces, \(f\), using the Kelvin section and the first derivative, are given by (3). Note that symbols \(\pm\) and \(\mp\) correspond to the cases of normal and shear direction, respectively (For example, \(\pm\) means \(+\) for normal direction and \(-\) for shear direction.).

Also, using stiffness \(K_m\) and viscosity of the Maxwell section:

Using Eqs (2) through (4), the second-order differential equation for contact force \(f\) is given by:

The force at a given step can be updated based on its value at the previous step and on values of the displacements at the current step and at the previous steps, using a finite dfifference scheme described below.

## Activity-Deletion Criteria

A contact with the Burger’s model is active if and only if the contact gap (\(g_{c}\)) is less than or equal to zero. The force-displacement law is skipped for inactive contacts.

## Force-Displacement Law

The force-displacement law for the Burger’s model updates the contact force and moment:

where \(\mathbf{F}\) combines the contributions from the Kelvin and Maxwell models acting in series. \(\mathbf{F}\) is resolved into normal and shear forces:

where \(F_{n} >0\) is tension. The shear force lies on the contact plane, and is expressed in the contact plane coordinate system:

From the second equation in (3) of the Kelvin section:

By using a central difference approximation of the finite difference scheme for the time derivative and taking average values for \(u_k\) and \(f\):

Therefore,

where:

For the Maxwell section, the displacement and the first derivative are given by:

Substituting the second and fourth equations of (4) into the second equation above:

By using a central difference approximation of the finite difference scheme and taking the average value for \(f\),

Therefore,

The total displacement and the first derivative of the Burger’s model are given by:

By using the finite difference scheme for the time derivative,

Substituting Eqs (11) and (16) into the equation above, the contact force, \(f^{t+1}\), is given by:

Where:

The contact force \(f^{t+1}\) can be calculated from know values for \(u^{t+1}\), \(u^t\), \(u_k^t\) and \(f^t\).

The force-displacement law for the Burger’s model consists of the following steps.

Update the normal force \(F_{n}\) according to Eq. (19).

Update the shear force as follows:

Update the shear force \(\mathbf{F}_{s}\) according to Eq. (19).

Compute the shear strength:

(21)\[F_{s}^{*} = - f_s F_{n} .\]Update the linear shear force:

(22)\[\begin{split}\mathbf{F}_{s} = \left\{\begin{array}{l} {\qquad \mathbf{F}_{s} ,\qquad \left\| \mathbf{F}_{s} \right\| \le F_{s}^{*} } \\ {\qquad F_{s}^{*} \left({\mathbf{F}_{s} \mathord{\left/ {\vphantom {\mathbf{F}_{s} \left\| F_{s}^{*} \right\| }} \right.} \left\| \mathbf{F}_{s} \right\| } \right),\qquad {\rm otherwise.}} \end{array}\right.\end{split}\]Update the slip state:

(23)\[\begin{split}s=\left\{\begin{array}{l} {\qquad {\rm true},\qquad \left\| \mathbf{F}_{s} \right\| = F_{s}^{*} } \\ {\qquad {\rm false},\qquad {\rm otherwise.}} \end{array}\right.\end{split}\]If the slip state is true, then the contact is sliding. Whenever the slip states changes, the slip_change callback event occurs.

## Energy Partitions

The Burger’s model does not provide any energy partition.

## Properties

The properties defined by the Burger’s contact model are listed in the table below as a concise reference; see the Contact Properties section for a description of the information in the table columns.

Keyword | Symbol | Description | Type | Range | Default | Modifiable | Inheritable |
---|---|---|---|---|---|---|---|

burger | Model name | ||||||

bur_knk | \(K_{k_n}\) | Normal stiffness Kelvin section [force/length] | FLT | \([0.0,+\infty)\) | 0.0 | YES | NO |

bur_cnk | \(C_{k_n}\) | Normal viscosity Kelvin section [force×time/length] | FLT | \([0.0,+\infty)\) | 0.0 | YES | NO |

bur_knm | \(K_{m_n}\) | Normal stiffness Maxwell section [force/length] | FLT | \([0.0,+\infty)\) | 0.0 | YES | NO |

bur_cnm | \(C_{m_n}\) | Normal viscosity Maxwell section [force×time/length] | FLT | \([0.0,+\infty)\) | 0.0 | YES | NO |

bur_ksk | \(K_{k_s}\) | Shear stiffness Kelvin section [force/length] | FLT | \([0.0,+\infty)\) | 0.0 | YES | NO |

bur_csk | \(C_{k_s}\) | Shear viscosity Kelvin section [force×time/length] | FLT | \([0.0,+\infty)\) | 0.0 | YES | NO |

bur_ksm | \(K_{m_s}\) | Shear stiffness Maxwell section [force/length] | FLT | \([0.0,+\infty)\) | 0.0 | YES | NO |

bur_csm | \(C_{m_s}\) | Shear viscosity Maxwell section [force×time/length] | FLT | \([0.0,+\infty)\) | 0.0 | YES | NO |

bur_fric | \(f_s\) | Friction coefficient [-] | FLT | \([0.0,+\infty)\) | 0.0 | YES | NO |

bur_mode | \(M_t\) | Normal-force tensile mode [-] | INT | \(\{0,1\}\) | 0 | YES | NO |

\(\;\;\;\;\;\;\begin{cases} \mbox{0: with tensile force} \\ \mbox{1: without tensile force} \end{cases}\) | |||||||

bur_slip | \(s\) | Slip state [-] | BOOL | {false,true} | false | NO | N/A |

\(\;\;\;\;\;\;\begin{cases} \mbox{true: slipping} \\ \mbox{false: not slipping} \end{cases}\) | |||||||

bur_force | \(\mathbf{F}\) | Total force (contact plane coord. system) | VEC | \(\mathbb{R}^3\) | \(\mathbf{0}\) | NO | NO |

\(\left( -F_n,F_{ss},F_{st} \right) \quad \left(\mbox{2D model: } F_{ss} \equiv 0 \right)\) |

Surface Property Inheritance

The Burger’s model does not provide property inheritance capabilities.

## Methods

The Burger’s model does not provide any method.

## Callback Events

Event | Array Slot | Value Type | Range | Description |
---|---|---|---|---|

contact_activated | Contact has become 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 | |

\(\;\;\;\;\;\;\begin{cases} \mbox{0: slip has initiated} \\ \mbox{1: slip has ended} \end{cases}\) |

## Usage and Verification Examples

The Verification Problem “Burger’s Contact Model: Stress Relaxation” compares the time-decay of the normal force to the expected analytical solution.

Was this helpful? ... | FLAC3D © 2019, Itasca | Updated: May 17, 2022 |