Lagrange Multiplier Method

Lagrange multipliers can be used to find the extreme of a multivariate function f ( x 1 , x 2 , ... , x n ) subject to the constraint g ( x 1 , x 2 , ... , x n ) = 0

Where, f and g are functions with continuous first partial derivatives on the open set containing the constraint curve, and g 0 at any point on the curve (where is the gradient).

To find the extreme, write:(1)
d f = f x 1 d x 1 + f x 2 d x 2 + ... + f x n d x n = 0
But, because g is being held constant, it is also true that(2)
d g = g x 1 d x 1 + g x 2 d x 2 + ... + g x n d x n = 0
So multiply Equation 2 by the as yet undetermined parameter λ and add to Equation 2,(3)
( f x 1 + λ g x 1 ) d x 1 + ( f x 2 + λ g x 2 ) d x 2 + ... + ( f x n + λ g x n ) d x n = 0
Note that the differentials are all independent, so any combination of them can be set equal to 0 and the remainder must still give zero. This requires that:(4)
( f x k + λ g x k ) d x k = 0
for all k = 1, ..., n, and the constant λ is called the Lagrange multiplier. For multiple constraints, g 1 = 0 , g 2 = 0 , ...,(5)
f = λ 1 g 1 + λ 2 g 2 + ...
The Lagrange multiplier method can be applied to contact-impact problems. In this case, the multivariate function is the expression of total energy subjected to the contact conditions:(6)
f ( x 1 , x 2 , ... , x n ) Π ( x , x ˙ , x ¨ )
(7)
g ( x 1 , x 2 , ... , x n ) Q ( x , x ˙ , x ¨ ) = 0
Where, x , x ˙ , x ¨ are the global vectors of DOF. The application of Lagrange multiplier method to the previous equations gives the weak form as:(8)
M x ¨ + f i n t f e x t + L λ = 0
with (9)
L x = b
This leads to:(10)
[ K L T L 0 ] { x λ } = { f 0 }

The Lagrange multipliers are physically interpreted as surface tractions. The equivalence of the modified virtual power principle with the momentum equation, the traction boundary conditions and the contact conditions (impenetrability and surface tractions) can be easily demonstrated. 1

It is emphasized that the above weak form is an inequality. In the discretized form, the Lagrange multiplier fields will be discretized and the restriction of the normal surface traction to be compressive will result from constraints on the trial set of Lagrange multipliers.

1 Engelmann B.E. and Whirley R.G., “A new elastoplastic shell element formulation for DYNA3D”, Report ugrl-jc-104826, Lawrence Livermore National Laboratory, 1990.