<figure><img src="image_1.png" /></figure>

HAL
open science

# Gauss' Principle and Principle of Least Constraints for Dissipative Mechanical Systems

Kerim Yunt

► To cite this version:

Kerim Yunt. Gauss' Principle and Principle of Least Constraints for Dissipative Mechanical Systems. 7th Vienna International Conference on Mathematical Modelling, Feb 2012, Vienne, Austria. pp.842 - 847, 10.3182/20120215-3-AT-3016.00149 . hal-01375987

HAL Id: hal-01375987

https://hal.science/hal-01375987

Submitted on 5 Oct 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Gauss' Principle and Principle of Least Constraints for Dissipative Mechanical Systems

## Kerim Yunt *

* P.O Box 1070 8021 Zurich, Switzerland (e-mail: kerimyunt@web.de).

**Abstract:** The aim of this work is to formulate the Gauss' principle and the principle of least constraints for dissipative systems. For dynamics, where the evolution requires the determination of the accelerations of the system, it is shown that in the presence of dissipative force laws a similar principle holds, which requires the augmentation of the optimization problem of least constraints by the time rate of change of the total dissipation.

**Keywords:** Friction, Dissipative Systems, Principle of Gauss, Variational Inequalities, Nonsmooth Analysis

## 1. INTRODUCTION

The principle of Gauss has been published in Gauss (1829) and provides a variational criteria on admissible accelerations of bilaterally holonomically constraint mechanical systems. The principle of Gauss, that bears the principle of least constraints, enables in addition to the above cited principles to determine the evolution of a non-impulsive mechanical process. The evolution of the mechanical system is determined by its acceleration and the principle of least constraints is about an optimization problem to determine accelerations for bilaterally constraint holonomic mechanical systems. In references Moreau (1963) and Moreau (1966), it is shown how to generalize the principle of Gauss and the principle of least constraints to finite-dimensional Lagrangian systems with frictionless unilateral constraints and the quadratic programming procedure for determining the accelerations for such mechanical systems is derived. Pozharitskii formulated in Pozharitskii (1961) the extension of the principle of least constraints for mechanical systems with planar isotropic dry friction with known normal forces by a function, he called the work of all forces applied to the system with a virtual acceleration. In Glocker (1998) a general optimization problem by putting the whole constraints in to an acceleration potential $\Psi_C(\ddot{q}^+)$ in the form:

$$f(\ddot{q}^+) = \frac{1}{2} \langle \ddot{q}^+, M \ddot{q}^+ \rangle - \langle h, \ddot{q}^+ \rangle + \Psi_C(\ddot{q}^+), \quad (1)$$

is proposed, which depending on the structure of the domain $C$ is a strictly convex function, with $M$ being the symmetric and positive-definite mass matrix and $h$ a continuous function of the generalized positions and velocities of the mechanical system but is not applicable to systems with friction. Pozharitskii considers following form of extension to the principle of least constraints in Pozharitskii (1961) for problems with tresca-type dry friction:

$$\min_{\ddot{q}^+} \frac{1}{2} \langle \ddot{q}^+, M \ddot{q}^+ \rangle - \langle h, \ddot{q}^+ \rangle + \Psi. \quad (2)$$

Pozharitskii defines the quantity:

$$-\Psi = \sum_{i=1}^{n} Q_i a_i - \sum_{i \in I_{NP}} \mu_i k_{ni} |a_{Ti}| \quad (3)$$

as the work of all the forces applied to the system with a virtual acceleration, where $Q$ is a generalized force and the distance function is in the Euclidean norm; $a_{Ti}$ is the relative tangential acceleration at a sticking contact; $k_{Ni}$ is the given normal force at a frictional. The distance of dissipation potentials are formulated in the Euclidean norm, for which norm he shows that the function is positive-definite.

The considered frictional contact law in this work, is the dry friction law of tresca type. There are two main features of dry friction. The first feature is a set $C_T$ to which frictional contact force $\gamma_T$ belongs as admissible forces. This set depends on the magnitude of the normal contact force $\gamma_N$. The friction coefficient $\mu$ serves as a transmissive multiplier of the normal force on the tangential direction. The second feature is about conditions which specify when sticking and sliding take place, which are related to the tangential relative contact velocity $v_T$ and acceleration $a_T$. The possibility to formulate friction models on acceleration level is shown systematically by Glocker in his monograph Glocker (2001) using nonsmooth analysis though already in Pozharitskii (1961) by Pozharitskii the spatial friction force model of tresca-type on acceleration level as well as friction acceleration potentials are used without definition, in order to extend the principle of Gauss to systems with dry friction. Though the distinction between Tresca and Coulomb type dry friction is made regularly in literature, there is no common view on how these friction laws differ from each other. In this work, Tresca type dry friction is seen as a friction law of which normal contact force is independent of future generalized accelerations at a given position and generalized velocity of the system contrary to Coulomb friction.

In this work, for mechanical systems with convex dissipation potentials (such as tresca-type of friction, signorini type normal contact law, viscous friction etc.) following
extended problem of least constraints is proposed:

$$ \min_{\ddot{q}^+} \frac{1}{2} \langle \ddot{q}^+, M \ddot{q}^+ \rangle - \langle h, \ddot{q}^+ \rangle + \dot{P}_D, \quad (4) $$

where $\dot{P}_D$ is the total time rate of change of the total dissipation function. The proposed Gauss principle associated with (4) is presented as well. In Yunt (2011) the relation between the time-rate of change of the total dissipation and the dual principal of least constraints is analyzed, and it is shown, that the dual principle of least constraints is obtained by the maximization of $\dot{P}_D$ with respect to passive dissipative forces. This is achieved by the techniques of nonsmooth analysis and by using the most general definition of dissipative forces indeed more general classes of forces then only friction forces are covered. In Pozharitskii (1961) an extensive analysis for the existence and uniqueness of the solution in generalized accelerations for the minimization of his proposed augmented function is provided, which is handled in this work by the properties of convexity in nonsmooth analysis. It is shown, that the extended problem of least constraints as in (4) possesses always a minimum due to the strict convexity of the goal functional. This relatively simple argumentation on the uniqueness of the minimum is one of the advantages one derives from nonsmooth analysis. The strict convexity is of the extended problem of least constraints requires that the dissipation potentials involved on velocity and acceleration level are convex. Though dissipation potentials due to viscous friction and nonconvex dissipation potentials are not addressed in this work, viscous friction can be treated straight forward with the results of this work, whereas the consideration of nonconvex dissipation potentials which often are encountered in continuum mechanics, of course makes the optimality conditions less sharp with regards to sufficiency.

The dissipation function was introduced by Rayleigh in his classical treatise for resisting forces which are linear in the velocities. Lurie extended this idea to dissipative forces which were higher powers of the velocity in his monograph Lurie (2002). He obtained by setting the power to zero the dissipation function for the dry friction forces and by setting to one the Rayleigh dissipation function. In modern mechanics, which is equipped with potential theory and nonsmooth analysis the dissipation functions of Rayleigh and Lurie are velocity potentials and need not to be differentiable. A superpotential is a lower-semicontinuous proper convex function $\Phi(\cdot)$ such that the inclusions $x \in \Phi(y)$ and $y \in \Phi^*(x)$ hold between two dual variables $x$ any $y$ and expresses some physical law. The second inclusion expresses the so-called inverse law and $\Phi^*(x)$ is the conjugate of $\Phi(y)$. In the sequel, superpotentials for normal and friction forces are introduced, which are related to dissipation functions in the classical sense. Further, these superpotential laws for frictional and normal unilateral contact are extended to the acceleration level in the sense of Glocker's work. The class of superpotentials is defined in Moreau (1968). It is known as stated in Acary and Brogliato (2008) that in the case of Tresca type friction in which the normal force does not depend on the dynamics, the friction force is modelled by using nonsmooth superpotentials which is not the case for Coulomb type dry friction. The research that has been conducted on systems with friction and contact mechanics is vast. A

broad overview on finite-dimensional contact mechanics is provided in Brogliato (1999), Brogliato (2007) and Glocker (2001) and in their cited references. The excellent reviews in Stewart (2000) and Stewart (2001) provide overviews on rigid body dynamics with friction and impact.

## 1.1 Preliminaries

Let $f$ be an extended-valued convex function from $\mathbb{R}^n$ to $[0, +\infty]$ and let $x$ be a point where $f$ is finite. The one-sided directional derivative of $f$ at $x$ with respect to a vector $y$ is defined as:

$$ f'(x; y) = \lim_{\gamma \downarrow 0} \frac{f(x + \gamma y) - f(x)}{\gamma}. \quad (5) $$

If it exists $+\infty$ and $\infty$ are allowed. A vector $x^*$ is a subgradient of a convex function $f$ at a point $x$ if

$$ f(z) \ge f(x) + \langle x^*, z - x \rangle, \quad \forall z. \quad (6) $$

The set of all subgradients of $f$ at $x$ constitute the subdifferential of $f$ at $x$ and is denoted by $\partial f(x)$ and is a closed convex set. Two important geometrical constructs are the tangent cone and the normal cone to a set $C$. The tangent cone $\mathcal{T}_C(x)$ to a set $C$ at a point $x \in C$ is a nonempty convex closed cone and is polar to a certain nonempty closed convex cone $\mathcal{N}_C(x)$:

$$ \mathcal{N}_C(x) = \{z | \langle y, z \rangle \le 0, \ y \in \mathcal{T}_C(x)\}, \quad (7) $$

$$ \mathcal{T}_C(x) = \{y | \langle y, z \rangle \le 0, \ z \in \mathcal{N}_C(x)\}. \quad (8) $$

The set $\mathcal{N}_C(x)$ is defined to be the normal cone to $C$ at $x$. The function $\Psi_C$ is the indicator function of the set $C$:

$$ \Psi_C(x) = 0, \quad \Leftrightarrow \quad x \in C, \quad (9) $$

$$ \Psi_C(x) = +\infty, \quad \Leftrightarrow \quad x \notin C. \quad (10) $$

An important special case in the theory of subgradients is the case, where $f$ is the indicator of a non-empty convex set $C$, then its subdifferential $\partial_x \Psi(x)$ is given by the normal cone at $x$ to the set $C$, $\partial_x \Psi(x) = \mathcal{N}_C(x)$. Let $f$ be any closed convex function on $\mathbb{R}^n$, then the conjugate of $f(x)$ is defined as:

$$ f^*(x^*) = \sup_x \{ \langle x, x^* \rangle - f(x) \}. \quad (11) $$

The conjugate of $f^*$ is $f$. The distance function $d_C$, which measures the distance to a closed set is given by:

$$ d_C(x) = \inf \{ \|x - s\| \mid s \in C \}. \quad (12) $$

The conjugate of the distance function to a convex set $d_C^*$ is given by:

$$ d_C^*(x) = \sup_y \{ \langle x, y \rangle - d_C(y) \} = \Psi_C^*(x) + \Psi_B(x). \quad (13) $$

The function $\Psi_C^*$ is the conjugate of the indicator function of the set $C$. Let $f$ be a proper convex function. The domain of $f$ is denoted as $\text{dom } f$ and is defined as the region where $f$ is finite. For $x \notin \text{dom } f$, $\partial f(x)$ is empty. For $x$ in the relative interior of $\text{dom } f$, $\partial f(x)$ is non-empty, $f'(x; y)$ is closed and proper as a function of $y$, and

$$ f'(x; y) = \sup_{x^*} \{ \langle x^*, y \rangle | x^* \in \partial f(x) \} = \Psi_{\partial f(x)}^*(y) \quad (14) $$

is valid. Further, $\partial f(x)$ is a nonempty bounded set, if and only if $x$ is an element of the interior of the domain of $f$, in which case $f'(x; y)$ is finite for every direction $y$. The concept of variational inequalities has been an outcome of the nonsmooth analysis, which plays an important role in
<figure><img src="image_3.png" /><figcaption>Fig. 1. In one dimension: (a) the distance function to the origin, (b) conjugate of the distance function to the origin, (c) subdifferential of the distance, (d) subdifferential of the conjugate.</figcaption></figure>

applied mathematics since its introduction by Hartman and Stampacchia Hartman and Stampacchia (1966). F. H. Clarke extended the theory to non-convex functionals, that are merely lower semicontinuous and used it to derive necessary conditions for non-smooth, non-convex problems and in Clarke (1990) an overview of his contributions are given. Detailed reading on various definitions in subdifferential calculus of convex extended-valued functionals is given in Rockafellar (1970) by Rockafellar.

## 2. DISSIPATIVE FORCE LAWS AND DISSIPATIVE FORCE POTENTIALS

A force $\gamma$ is called a dissipative force, if the power of the force element $P_{\gamma}(v)$ fulfills:

$$P_{\gamma}(v) := \langle \gamma(v), v \rangle \leq 0. \quad (15)$$

The nonpositivity of the power of the force element is visualized best if a typical friction force law is considered. During the sticking phase, the power dissipation is zero; and is negative during the sliding phase. If there exists dissipative superpotentials in the sense of Moreau's work Moreau (1973), then there exists a scalar potential $\Phi$, which is a lower-semicontinuous convex function with $\text{dom} \Phi \neq \emptyset$, such that the dissipative force is given by:

$$\gamma \in \partial_v \Phi(v). \quad (16)$$

The conjugate of the scalar potential $\Phi^*$ relates the force and the velocity in the following form:

$$v \in \partial_{\gamma} \Phi^*(\gamma). \quad (17)$$

Let $q, \dot{q}, \ddot{q}$ represent the position, velocity and acceleration in the generalized coordinates of the basic configuration of a finite-dimensional Lagrangian system with n maximal degrees of freedom (DOF), respectively. A contact between two rigid bodies is seen as a point-to-point contact. Let the scalar function $d_N(q)$ denote the shortest normal distance

between two rigid bodies in the system as depicted in figure (3). The normal contact distances are always nonnegative due to the impenetrability at rigid body contacts. Further, there exist no attracting forces at the rigid body contacts. A contact is closed, if the normal contact distance $d_N$, normal contact velocity $v_N$ are both zero. The relation between these kinematic entities is given as follows:

$$\dot{d}_N = v_N = \frac{\partial d_N}{\partial q} \dot{q} = J_N^T \dot{q}. \quad (18)$$

The normal contact acceleration is analogously given by:

$$\dot{v}_N = J_N^T \ddot{q} + \dot{J}_N^T \dot{q} = J_N^T \ddot{q} + p_N. \quad (19)$$

The tangential relative contact velocity and acceleration are defined by:

$$v_T = J_T^T \dot{q} \quad \text{and} \quad \dot{v}_T = J_T^T \ddot{q} + \dot{J}_T^T \dot{q} = J_T^T \ddot{q} + p_T, \quad (20)$$

respectively. The friction force exists on a disk with normal contact force dependent radius in $\mathbb{R}^2$. The friction force is limited to the following set:

$$\gamma_T \in -\mu \gamma_N \mathbb{B}_2 = \mathcal{C}_T, \quad (21)$$

where $\mathbb{B}_2$ is the unit ball $\mathbb{R}^2$. Let for a given positive $\gamma_N$ and $\mu$, $f_T$ be defined as:

$$f_T = -\frac{\gamma_T}{\mu \gamma_N}. \quad (22)$$

It is obvious that $f_T \in \mathbb{B}_2$. If at any instant the dissipativity condition $\langle \gamma_T, v_T \rangle \leq 0$ is valid, then one has by the definition of the normal cone the following inclusion:

$$v_T \in \mathcal{N}_{\mathbb{B}_2}(f_T). \quad (23)$$

The subdifferential of the indicator function of a convex set is given by its normal cone at that location. By making use of this relation, one has:

$$v_T \in \partial \Psi_{\mathbb{B}_2}(f_T). \quad (24)$$

According to relation (13) we have $\Psi_{\mathbb{B}_2}^*(v_T) = d_{\{0\}}(v_T)$, which is the distance function of the tangential velocity to the origin. By this relation the following is valid:

$$-\gamma_T \in \mu \gamma_N \partial d_{\{0\}}(v_T). \quad (25)$$

The function $\mu \gamma_N d_{\{0\}}(v_T)$ is a nonsmooth function, which is the dissipation function for spatial friction in this work. It encompasses also the behaviour of the frictional contact at sticking which is enabled by its nonsmoothness at the origin. However, it serves as a superpotential in the sense of Moreau only in case of dry friction where the normal force is known. The relation between the kinematics in normal direction and normal contact force $\gamma_N$ for a single contact is expressed as Signorini type complementarity relations:

$$d_N > 0 \Rightarrow \gamma_N = 0, \quad (26)$$

$$d_N = 0, \quad v_N > 0 \Rightarrow \gamma_N = 0, \quad (27)$$

$$d_N = 0, \quad v_N = 0, \quad \dot{v}_N = 0 \Rightarrow \gamma_N \geq 0, \quad (28)$$

$$d_N = 0, \quad v_N = 0, \quad \dot{v}_N > 0 \Rightarrow \gamma_N = 0. \quad (29)$$

It is easily verified that the Signorini type unilateral contact law defines a dissipative force law if one notices that $\langle v_N, -\gamma_N \rangle = 0$ is valid. The dissipation condition requires $\langle v_N, \gamma_N \rangle \leq 0$ to hold, which is moulded into a normal cone inclusions:

$$v_N \in \mathcal{N}_{\mathbb{R}_0^{-}}(-\gamma_N), \quad -\gamma_N \in \mathcal{N}_{\mathbb{R}_0^{+}}(v_N). \quad (30)$$

The velocity superpotential of the normal contact force law is:

$$\Phi(v_N) = \Psi_{\mathbb{R}_0^{+}}(v_N). \quad (31)$$
<figure><img src="image_4.png" /><figcaption>Fig. 2. The spatial frictional disc $C_T(\gamma_N)$ with normal force dependent radius and various normal cones (red).</figcaption></figure>

This definition of dissipativity encompasses the normal force potential formulated on velocity level, which is a force law with exactly zero dissipated power.

The nonsmooth dissipation potentials capture the behaviour of the dissipative force laws also at the origin. Any relative velocity is approximated by the truncated Taylor series for small enough $\tau$:

$$v(t + \tau) = v(t) + \dot{v}(t) \tau. \quad (32)$$

The generalized directional derivative of the function $d_{\{0\}}(v_T)$ at the origin in the direction $\dot{v}_T^+(t)$ is given by:

$$d'_{\{0\}}(v_T; \dot{v}_T^+) = \limsup_{\substack{v_T \to 0 \\ \tau \downarrow 0^+}} \frac{d_{\{0\}}(v_T + \dot{v}_T^+ \tau) - d_{\{0\}}(v_T)}{\|\dot{v}_T^+ \tau\|} = 1. \quad (33)$$

This result shows that locally around $v_T = 0$, one has $d_{\{0\}}(v_T) = d_{\{0\}}(\dot{v}_T^+)$. The upper subderivative of the normal force velocity potential $\Psi_{\mathbb{R}_0^+} (v_N)$ at the origin in any direction $\dot{v}_N(t)$ is:

$$\Psi'_{\mathbb{R}_0^+} (v_N = 0; \dot{v}_N^+) = \begin{cases} \infty, & \text{if } \dot{v}_N^+ \notin \mathcal{T}_{\mathbb{R}_0^+} (0), \\ 0, & \text{if } \dot{v}_N^+ \in \mathcal{T}_{\mathbb{R}_0^+} (0). \end{cases} \quad (34)$$

The tangent cone $\mathcal{T}_{\mathbb{R}_0^+}$ at 0 is given by $\mathbb{R}_0^+$. So the upper subderivative becomes:

$$\Psi'_{\mathbb{R}_0^+} (v_N = 0; \dot{v}_N^+) = \Psi_{\mathbb{R}_0^+} (\dot{v}_N^+) \quad (35)$$

The following relations hold based on the definition of the one-sided directional derivative:

$$\Psi'_{\mathbb{R}_0^+} (0; \dot{v}_N^+) = \Psi_{\mathbb{R}_0^+} (\dot{v}_N^+) = \sup_{\gamma_N} \{-\langle \gamma_N, \dot{v}_N^+ \rangle, -\gamma_N \in \mathbb{R}_0^-\}, \quad (36)$$

$$d'_{\{0\}}(0; \dot{v}_T^+) = d_{\{0\}}(\dot{v}_T^+) = \sup_{f_T} \{\langle f_T, \dot{v}_T^+ \rangle, f_T \in \mathbb{B}_2\}. \quad (37)$$

According to (14) the following relations hold:

<figure><img src="image_5.png" /><figcaption>Fig. 3. The geometry of a single spatial rigid body contact between two rigid bodies.</figcaption></figure>

$$\Psi_{\mathbb{R}_0^+} (\dot{v}_N^+) = \sup_{\gamma_N} \{-\langle \gamma_N, \dot{v}_N^+ \rangle, -\gamma_N \in \partial \Psi_{\mathbb{R}_0^+} (v_N^+)\}, \quad (38)$$

$$d_{\{0\}}(\dot{v}_T^+) = \sup_{f_T} \{\langle f_T, \dot{v}_T^+ \rangle, f_T \in \partial d_{\{0\}}(v_T^+)\}. \quad (39)$$

Due to conjugation following equalities hold for $v_N = 0$ and $v_T = 0$:

$$\Psi_{\mathbb{R}_0^+} (\dot{v}_N^+) + \Psi_{\mathbb{R}_0^-} (-\gamma_N) = -\langle \dot{v}_N^+, \gamma_N \rangle, \quad (40)$$

$$d_{\{0\}}(\dot{v}_T^+) + \Psi_{\mathbb{R}_0^-} (-\gamma_T) = -\langle \dot{v}_T^+, f_T \rangle. \quad (41)$$

The time-rate of change of the friction force dissipation function except at the origin is:

$$\dot{d}_{\{0\}}(v_T) = \left\langle \frac{v_T}{|v_T|}, a_T \right\rangle = \left\langle -\frac{\gamma_T(v_T)}{\mu \gamma_N}, a_T \right\rangle. \quad (42)$$

At the origin it is given by:

$$\frac{d}{dt} (\mu \gamma_N d_{\{0\}}(v_T))|_{v_T=0} = \mu \gamma_N d_{\{0\}}(a_T). \quad (43)$$

The derivative of the normal contact velocity potential with respect to time is:

$$\dot{\Psi}_{\mathbb{R}_0^+} (v_N) = \langle a_N, -\gamma_N \rangle, \quad -\gamma_N \in \partial \Psi_{\mathbb{R}_0^+} (v_N). \quad (44)$$

## 3. MATHEMATICAL MODELLING OF MECHANICS

In order to formulate the index sets that account for the force laws properly, the following index sets are introduced:

$$I_S = \{i | d_{N_i} = 0\}, \quad I_N = \{i | d_{N_i} = 0, \ v_{N_i} = 0\}.$$

The set $I_S$ denotes the set of all contacts that are closed on position level of the system. The set $I_N$, with a number of k elements, denotes the set of all contacts at which normal contact velocity and normal contact distance equal to zero. In this case the active dissipative friction forces are the contacts at which sliding takes place. The set of passive friction forces consist of those which stick ($v_T = 0$). The index sets of sliding contacts and sticking
contacts are named I<sub>NA</sub> and I<sub>NP</sub>, respectively. The normal contact force is passive if the associated normal contact velocity is zero. The sticking and the sliding contact forces are incorporated by the appropriate generalized force directions in the equations of motion, which are the columns of the linear operators J<sub>TA</sub> and J<sub>TP</sub>, respectively. For a finite-dimensional Lagrangian system with n degrees of freedom and 1 active sliding contacts and m passive sticking contacts, which have spatial friction, J<sub>TA</sub>(q) is an n by 2l and J<sub>TP</sub>(q) is an n by 2m linear operators of generalized friction force directions. The linear operator J<sub>N</sub>(q) is an n by k linear operator of generalized normal force directions. Given this setting the differential inclusion of a mechanical system with tresca-type friction at the contacts is stated as:

$$
M\ddot{q} - h - J_{TA}\gamma_{TA} - J_{TP}\gamma_{TP} - J_{N}\gamma_{N} = 0, \quad (45)
$$

$$
- \gamma_{N_i} \in \partial \Psi_{\mathbb{R}_0^+} (a_{N_i}), \quad \forall i \in I_N, \qquad (46)
$$

$$
\gamma_{TA_i} = -\mu_i k_{N_i} \frac{v_{T_i}}{\|v_{T_i}\|}, \quad \forall i \in I_{NA}, \qquad (47)
$$

$$
- \gamma_{TP_i} \in \mu_i k_{N_i} \partial_{a_{TP_i}} d_{\{0\}}(a_{TP_i}), \forall i \in I_{NP}. \quad (48)
$$

M(q) is the symmetric positive-definite (PD) mass matrix and h(q, ˙q) represents the vector of gyroscopic, centripetal and coriolis forces and encompasses also smooth force elements such as springs and dampers. Here k<sub>Ni</sub> are the normal contact forces at the tresca-type frictional contacts for which the normal force is not dependent on the future generalized accelerations at a given state.

Having set the stage, the total dissipation function is given
as:

$$
P_D(v_T, v_N) = \sum_{i \in I_S} \Psi_{\mathbb{R}_0^+} (v_{n_i}) + \sum_{i \in I_{NP} \cup I_{NA}} \mu_i k_{n_i} d_{\{0\}} (v_{T_i}). \quad (49)
$$

The time rate of change $\dot{P}_D$ is evaluated by making use of
(42), (43), (44) as follows:

$$
\begin{equation}
\begin{split}
\dot{P}_D = \sum_{i \in I_N} \Psi_{\mathbb{R}_0^+} (a_{n_i}) & \\
& \quad + \sum_{i \in I_{NP}} \mu_i k_{n_i} d_{\{0\}} (a_{T_i}) - \sum_{i \in I_{NA}} \langle \gamma_{T_i} (v_{T_i}), a_{T_i} \rangle .
\end{split}
\tag{50}
\end{equation}
$$

4. PRINCIPLE OF LEAST CONSTRAINTS FOR
DISSIPATIVE MECHANICAL SYSTEMS

The well-known principle of maximal dissipation does not
specify what dissipative forces shall do if the relative ve-
locity vanishes at given generalized velocity and position,
though this maximization of the dissipativity happens re-
currently at every position and velocity of the mechanical
system. On the other side, the dissipative forces, influence
the evolution of the system by determining the acceler-
ation $\ddot{q}$ in general. The generalized acceleration of the
mechanical system is related to dissipative forces by the
set of relations given by (45) to (48) for a given state
specified by the position $q$ and the generalized velocity
$\dot{q}$. Any optimization in the passive dissipative forces can
thereby related to a search of the generalized acceleration
in the next moment. Following extended problem of least
constraints is proposed in this work:

$$
\min_{\ddot{q}^{+}} \frac{1}{2} \langle \ddot{q}^{+}, M \ddot{q}^{+} \rangle - \langle h, \ddot{q}^{+} \rangle + \dot{P}_{D}. \quad (51)
$$

The unconstrained minimization problem (51) is the principle of least constraint for the considered type of mechanical systems. Here $\ddot{q}^+$ is the right-continuous generalized acceleration, which is the acceleration in the imminent future. The necessary condition of optimality of (51) takes following form:

$$
- (M(q) \ddot{q}^{+} - h(q, \dot{q})) \in \partial_{\ddot{q}^{+}} \dot{P}_{D} \quad (52)
$$

The directional derivative of the unconstrained goal func-
tional requires for a minimum following condition:

$$
\dot{P}'_{D}(\cdot; \tilde{\ddot{q}}^{+} - \ddot{q}^{+}) + \langle M \ddot{q}^{+} - h, \tilde{\ddot{q}}^{+} - \ddot{q}^{+} \rangle \geq 0. \quad (53)
$$

The vector $\tilde{\ddot{q}}^+$ denotes the arbitrary unconstrained virtual future accelerations. The requirement of (53) represents the Gauss Principle for dissipative mechanical systems with convex dissipation potentials. The dependence of the relative accelerations as in $\dot{P}'_D$ on $\ddot{q}^+$ are given by relations (18) and (19). The directional derivative of $\dot{P}'_D$ becomes by making use of (18) and (19):

$$
\dot{P}'_D (\cdot; \tilde{\ddot{q}}^+ - \ddot{q}^+) = \langle -J_N \gamma_N - J_{TP} \gamma_{TP} - J_{TA} \gamma_{TA}, \tilde{\ddot{q}}^+ - \ddot{q}^+ \rangle. \quad (54)
$$

By combining (53) and (54) following optimality condi-
tion is obtained:

$$
\langle M \ddot{q}^+ - h - J_N \gamma_N - J_{TP} \gamma_{TP} - J_{TA} \gamma_{TA}, \tilde{\ddot{q}}^+ - \ddot{q}^+ \rangle \ge 0, \quad (55)
$$

which for arbitrary unconstrained virtual future accelera-
tions $\tilde{\ddot{q}}^+$ is only fulfilled, if the equations of motion given
by (45) hold.

5. DISCUSSION AND CONCLUSION

In this work it is proven mathematically, that the principle
of least constraints for mechanical systems with convex
dissipation potentials is given by (51). The principle of
Gauss for mechanical systems with tresca type of dry
friction is shown to be related to the time-rate of change of
the total dissipation. The new principle can be extended
to include general dissipative force laws, which derive from
convex potentials. This relation is accomplished by making
use of the most general definition of dissipative forces
and applying the methods and techniques of nonsmooth
analysis.

REFERENCES

Acary, V. and Brogliato, B. (2008). *Numerical Methods for Nonsmooth Dynamical Systems-Applications in Mechanics and Electronics*, volume 35 of Lecture Notes in Applied and Computational Mechanics. Springer.

Brogliato, B. (1999). *Nonsmooth Mechanics: Models, Dynamics and Control*. Communications and Control Engineering. Springer, 2nd edition.

Brogliato, B. (2007). *Dissipative Systems Analysis and Control: Theory and Applications*. Communications and Control Engineering. Springer, 2nd edition.

Clarke, F. (1990). *Optimisation and Nonsmooth Analysis*. Number 5 in Classics in Applied Mathematics. SIAM.

Gauss, C.F. (1829). Über ein neues allgemeines grundgesetz der mechanik. *J. Reine Angew. Math.*, 4, 232–235.
Glocker, C. (1998). The principles of d'alembert, jour-dain, and gauss in nonsmooth dynamics part 1: Scler-nomic multibody systems. *Zeitschrift für Angewandte Mechanik und Mathematik*, 78(1), 21–37.

Glocker, C. (2001). *Set-Valued Force Laws – Dynamics of Non-Smooth Systems*, volume 1 of *Lecture Notes in Applied Mechanics*. Springer Verlag.

Hartman, P. and Stampacchia, G. (1966). On some nonlinear elliptic differential equations. *Acta Math.*, 115(1), 153–188.

Lurie, A.I. (2002). *Analytical Mechanics*. Springer.

Moreau, J.J. (1963). Les liaisons unilatérales et le principe de gauss. *Ibid.*, 256, 871–874.

Moreau, J.J. (1966). Quadratic programming in mechanics: Dynamics of one-sided constraints. *J. SIAM. Control*, 4(1), 153–158.

Moreau, J.J. (1968). La notion de sur-potentiel et les liaisons unilatérales en élastostatique. *C. R. Acad. Sci. Paris*, A271, 954–957.

Moreau, J.J. (1973). *New variational techniques in Mathematical Physics*, chapter On unilateral constraints, friction and plasticity, 172–322.

Pozharitskii, G.K. (1961). Extension of the principle of gauss to systems with dry (coulomb) friction. *Journal of Appl. Math. and Mech.*, 25(3), 586–607.

Rockafellar, R.T. (1970). *Convex Analysis*. Princeton Landmarks in Mathematics and Physics. Princeton University Press.

Stewart, D.E. (2000). Rigid-body dynamics with friction and impact. *SIAM Review*, 42(1), 3–39.

Stewart, D.E. (2001). Finite-dimensional contact mechanics. *Phil. Trans. R. Soc. Lond. A*, 359, 2467–2482.

Yunt, K. (2011). On the relation between the dual principle of least constraints and the principle of maximal dissipation. In *Proc. of the 7th WSEAS Int. Conf. on Applied and Theoretical Mechanics*, 75–80. Montreux, Switzerland.