0000000000390363
AUTHOR
Jaroslav Haslinger
On State Constrained Optimal Shape Design Problems
This paper is concerned with the following optimal design problem with constraints both on the state and on the control: $$MinimizeJ(y,u)$$ (P) subject to $$A\left( u \right)y + \partial \varphi \left( y \right) \mathrel\backepsilon Bu + f,$$ (1.1) $$y \in K,$$ (1.2) $$u \in {U_{ad}}.$$ (1.3)
On the design of the optimal covering of an obstacle
We consider the problem of controlling the shape of the coincidence set in an obstacle problem. This so called packaging problem was introduced in
On a topology optimization problem governed by two-dimensional Helmholtz equation
The paper deals with a class of shape/topology optimization problems governed by the Helmholtz equation in 2D. To guarantee the existence of minimizers, the relaxation is necessary. Two numerical methods for solving such problems are proposed and theoretically justified: a direct discretization of the relaxed formulation and a level set parametrization of shapes by means of radial basis functions. Numerical experiments are given.
Time Dependent Case
This chapter is devoted to finite element approximations of scalar time dependent hemivariational inequalities. We start with the parabolic case following closely Miettinen and Haslinger, 1998. At the end of this chapter we discuss, how the results can be extended to constrained problems. Our presentation will follow the structure used for the static case in Chapter 3. First, we introduce an abstract formulation of a class of parabolic hemivariational inequalities (see Miettinen, 1996, Miettinen and Panagiotopoulos, 1999).
Shape optimization for Stokes problem with threshold slip boundary conditions
This paper deals with shape optimization of systems governed by the Stokes flow with threshold slip boundary conditions. The stability of solutions to the state problem with respect to a class of domains is studied. For computational purposes the slip term and impermeability condition are handled by a regularization. To get a finite dimensional optimization problem, the optimized part of the boundary is described by B´ezier polynomials. Numerical examples illustrate the computational efficiency. peerReviewed
A new incremental method of computing the limit load in deformation plasticity models
The aim of this paper is to introduce a new incremental procedure that can be used for numerical evaluation of the limit load. Existing incremental type methods are based on parametrization of the energy by the loading parameter $\zeta\in[0,\zeta_{lim})$, where $\zeta_{lim}$ is generally unknown. In the new method, the incremental procedure is operated in terms of an inverse mapping and the respective parameter $\alpha$ is changing in the interval $(0,+\infty)$. Theoretically, in each step of this algorithm, we obtain a guaranteed lower bound of $\zeta_{lim}$. Reduction of the problem to a finite element subspace associated with a mesh $\mathcal T_h$ generates computable bound $\zeta_{lim,h…
Computable majorants of the limit load in Hencky’s plasticity problems
Abstract We propose a new method for analyzing the limit (safe) load of elastoplastic media governed by the Hencky plasticity law and deduce fully computable bounds of this load. The main idea of the method is based on a combination of kinematic approach and new estimates of the distance to the set of divergence free fields. We show that two sided bounds of the limit load are sharp and the computational efficiency of the method is confirmed by numerical experiments.
Contact Shape Optimization
Shape optimization is a branch of the optimal control theory in which the control variable is connected with the geometry of the problem. The aim is to find a shape from an a priori defined class of domains, for wich the corresponding cost functional attains its minimum. Shape optimization of mechanical systems, behaviour of which is described by equations, has been very well analyzed from the mathematical, as well as from the mechanical point of view, see [1], [2], [3] and references therein. The aim of this contribution is to extend results to the case, in which the system is described by the so called variational inequalities. There are two reasons for doing that: 1) The behavior of many…
Shape Optimization in Contact Problems. 1. Design of an Elastic Body. 2. Design of an Elastic Perfectly Plastic Body
The optimal shape design of a two dimensional body on a rigid foundation is analyzed. The problem is how to find the boundary part of the body where the unilateral boundary conditions are assumed in such a way that a certain energy integral (total potential energy, for example) will be minimized. It is assumed that the material of the body is elastic. Some remarks will be given concerning the design of an elastic perfectly plastic body. Numerical examples will be given.
The parameter identification in the Stokes system with threshold slip boundary conditions
The paper is devoted to an identification of the slip bound function g in the Stokes system with threshold slip boundary conditions assuming that g depends on the tangential velocity 𝑢𝜏 . To this end the optimal control approach is used. To remove its nonsmoothness we use a regularized form of the slip conditions in the state problem. The mutual relation between solutions to the original optimization problem and the problems with regularized states is analyzed. The paper is completed by numerical experiments. peerReviewed
An abstract inf-sup problem inspired by limit analysis in perfect plasticity and related applications
This paper is concerned with an abstract inf-sup problem generated by a bilinear Lagrangian and convex constraints. We study the conditions that guarantee no gap between the inf-sup and related sup-inf problems. The key assumption introduced in the paper generalizes the well-known Babuška–Brezzi condition. It is based on an inf-sup condition defined for convex cones in function spaces. We also apply a regularization method convenient for solving the inf-sup problem and derive a computable majorant of the critical (inf-sup) value, which can be used in a posteriori error analysis of numerical results. Results obtained for the abstract problem are applied to continuum mechanics. In particular…
On One Identification Problem in Linear Elasticity
In practice we meet problems, when having the solution of partial differential equation, we want to discover parts in the domain of its definition where the solution has some specific properties. In [1] and [2] the problem of identification of a curve φ, lying inside of Ω such that the flux \(\int{_{\varphi }}\frac{\partial u}{\partial n}ds\) is maximal has been studied, where u is the solution of mixed—boundary value problem for Laplacian operator.
An abstract inf-sup problem inspired by limit analysis in perfect plasticity and related applications
This work is concerned with an abstract inf-sup problem generated by a bilinear Lagrangian and convex constraints. We study the conditions that guarantee no gap between the inf-sup and related sup-inf problems. The key assumption introduced in the paper generalizes the well-known Babuska-Brezzi condition. It is based on an inf-sup condition defined for convex cones in function spaces. We also apply a regularization method convenient for solving the inf-sup problem and derive a computable majorant of the critical (inf-sup) value, which can be used in a posteriori error analysis of numerical results. Results obtained for the abstract problem are applied to continuum mechanics. In particular, …
Nonsmooth Optimization Methods
From the previous chapters we know that after the discretization, elliptic and parabolic hemivariational inequalities can be transformed into substationary point type problems for locally Lipschitz superpotentials and as such will be solved. There is a class of mathematical programming methods especially developed for this type of problems. The aim of this chapter is to give an overview of nonsmooth optimization techniques with special emphasis on the first and the second order bundle methods. We present their basic ideas in the convex case and necessary modifications for nonconvex optimization. We shall use them in the next chapter for the numerical realization of several model examples. L…
Approximation of Elliptic Hemivariational Inequalities
From the previous chapter we know that there exist many important problems in mechanics in which constitutive laws are expressed by means of nonmonotone, possibly multivalued relations (nonmonotone multivalued stress-strain or reaction-displacement relations,e.g). The resulting mathematical model leads to an inclusion type problem involving multivalued nonmonotone mappings or to a substationary type problem for a nonsmooth, nonconvex superpotential expressed in terms of calculus of variation. It is the aim of this chapter to give a detailed study of a discretization of such a type of problems including the convergence analysis. Here we follow closely Miettinen and Haslinger, 1995, Miettinen…
Nonsmooth Mechanics. Convex and Nonconvex Problems
Nonlinear, multivalued and possibly nonmonotone relations arise in several areas of mechanics. A multivalued or complete relation is a relation with complete vertical branches. Boundary laws of this kind connect boundary (or interface) quantities. A contact relation or a locking mechanism between boundary displacements and boundary tractions in elasticity is a representative example. Material constitutive relations with complete branches connect stress and strain tensors, or, in simplified theories, equivalent stress and strain quantities. A locking material or a perfectly plastic one is represented by such a relation. The question of nonmonotonicity is more complicated. One aspect concerns…
On a topology optimization problem governed by two-dimensional Helmholtz equation
The paper deals with a class of shape/topology optimization problems governed by the Helmholtz equation in 2D. To guarantee the existence of minimizers, the relaxation is necessary. Two numerical methods for solving such problems are proposed and theoretically justified: a direct discretization of the relaxed formulation and a level set parametrization of shapes by means of radial basis functions. Numerical experiments are given. peerReviewed
Optimal shape design and unilateral boundary value problems: Part II
In the first part we give a general existence theorem and a regularization method for an optimal control problem where the control is a domain in R″ and where the system is governed by a state relation which includes differential equations as well as inequalities. In the second part applications for optimal shape design problems governed by the Dirichlet-Signorini boundary value problem are presented. Several numerical examples are included.
Inf-sup conditions on convex cones and applications to limit load analysis
The paper is devoted to a family of specific inf–sup conditions generated by tensor-valued functions on convex cones. First, we discuss the validity of such conditions and estimate the value of the respective constant. Then, the results are used to derive estimates of the distance to dual cones, which are required in the analysis of limit loads of perfectly plastic structures. The equivalence between the static and kinematic approaches to limit analysis is proven and computable majorants of the limit load are derived. Particular interest is paid to the Drucker–Prager yield criterion. The last section exposes a collection of numerical examples including basic geotechnical stability problems.…
Identification of critical curves. Part II. Discretization and numerical realization
We consider the finite element approximation of the identification problem, where one wishes to identify a curve along which a given solution of the boundary value problem possesses some specific property. We prove the convergence of FE-approximation and give some results of numerical tests. peerReviewed
Shape optimization in contact problems based on penalization of the state inequality
The paper deals with the approximation of optimal shape of elastic bodies, unilaterally supported by a rigid, frictionless foundation. Original state inequality, describing the behaviour of such a body is replaced by a family of penalized state problems. The relation between optimal shapes for the original state inequality and those for penalized state equations is established. peerReviewed