## mpccdist3 details:

Keywords: analytic solution

Geometry: easy, fixed

Design: coupled via volume data

Differential operator:

• Damage:
• quasi-linear other operator of order 2.
• Defined on a 2-dim domain in 2-dim space
• Time dependent.

Design constraints:

• none

State constraints:

• none

Mixed constraints:

• none

Submitted on 2016-07-01 by Marita Holtmannspötter. Published on 2016-12-23

## mpccdist3 description:

### Introduction

We study the optimal control of a particular gradient enhanced damage model. The damage model is based on Dimitrijevic and Hackl [2008], the optimization problem is currently unpublished and provided by M. Holtmannspötter.

The presented damage model features two damage variables, one with higher spatial regularity and one which carries the evolution of damage in time. Their diﬀerence is penalized in the free energy functional. The evolution of damage in time is modeled by a nonsmooth operator ODE. Therefore the control-to-state operator is not diﬀerentiable whenever there are biactive points present.

This example features three sets of data, such that the known global optimum has either only inactive points, only strongly active points, or only biactive points.

### Variables & Notation

#### Unknowns

In contrast to the original damage model in Dimitrijevic and Hackl [2008], the situation is simpliﬁed by considering only scalar valued functions. The unknown functions are

#### Given Data

For the construction of analytically known solutions, the following functions need to be speciﬁed:

These data will be speciﬁed below in the supplementary materials section, such that globally optimal solutions with the desired properties are known. In addition, the following data are needed to specify all variables in the problem:

In this example, the diﬃculty lies in the potential non-diﬀerentiability of the $\mathrm{max}\left\{0,\cdot \right\}$ operator in the evolution law of the local damage variable. To describe the area where diﬀerentiability is critical, we deﬁne three sets:
$\begin{array}{llll}\hfill \text{Minimize}\phantom{\rule{1em}{0ex}}& \frac{1}{2}\parallel u-{u}_{d}{\parallel }_{{L}^{2}\left(0,T;{L}^{2}\left(\Omega \right)\right)}^{2}+\frac{1}{2}\parallel \phi -{\phi }_{d}{\parallel }_{{L}^{2}\left(0,T;{L}^{2}\left(\Omega \right)\right)}^{2}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & \phantom{\rule{1em}{0ex}}+\frac{1}{2}\parallel d-{d}_{d}{\parallel }_{{L}^{2}\left(0,T;{L}^{2}\left(\Omega \right)\right)}^{2}+\frac{1}{2}\parallel \ell -{\ell }_{d}{\parallel }_{{L}^{2}\left(0,T;{L}^{2}\left(\Omega \right)\right)}^{2}.\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$