Modified Duality Scheme for Solving Elastic Crack Problem Marina V. Chervyakova Robert V. Namm Pacific National University Computing Center Ulica Tikhookeanskaya 136, of Far Eastern Branch 680035 Khabarovsk, Russia Russian Academy of Sciences m.tchervyakova@yandex.ru Ulica Kim Yu Chen 65, 680000 Khabarovsk, Russia rnamm@yandex.ru Ellina M. Vikhtenko Pacific National University Ulica Tikhookeanskaya 136, 680035 Khabarovsk, Russia vikht.el@gmail.com Abstract We consider the duality method for solving a model f elastic problem with a crack based on the use of duality methods. An article presents the theorems, allowing to use Uzawa method for search a saddle point of the modified Lagrangian functional. The results of numerical exper- iments are given. 1 Introduction Many problems of mechanics and physics are formulated in the form of nonlinear boundary value problems for partial differential equations with additional conditions imposed on the desired solution. In the construction of more complex models, where the boundary conditions are often the system of equations and inequalities, in most cases the exact solution can not be found. Application of variational methods allows to pass from the nonlinear boundary value problems of mathematical physics to minimizing problem, consisting in research a conditional extremal problem [Hlavachek et al., 1988, Kikuchi & Oden, 1988]. So, it is very important to develop efficient methods for solving the original problem. This paper is devoted to the analysis of model, a wide interest in which is manifested in the recent years. It is a model problem of elastic deformation of the body containing a crack. The formulation of this type of problem can be found in monograph by A.M. Khludnev [Khludnev, 2010]. Model proposed here differs from classical approach to the crack problem because it is characterized by the nonlinear boundary conditions on crack faces. Suitable boundary conditions are written as inequalities which provide mutual nonpenetration between crack faces. From the standpoint of mechanics such models are more preferable than the linear classical models. Copyright ⃝ c by the paper’s authors. Copying permitted for private and academic purposes. In: Yu. G. Evtushenko, M. Yu. Khachay, O. V. Khamisov, Yu. A. Kochetov, V.U. Malkova, M.A. Posypkin (eds.): Proceedings of the OPTIMA-2017 Conference, Petrovac, Montenegro, 02-Oct-2017, published at http://ceur-ws.org 139 The methodological basis of study is duality scheme, based on replacing the problem of constrained op- timization to finding a saddle point of the modified Lagrangian functional. This method was considered to the finite-dimensional problems of linear and convex programming in books by authors E.G. Gol’stein and N.V. Tret’yakov [Gol’stein & Tret’yakov, 1989], A.A. Kaplan and Ch. Grossman [Grossman & Kaplan, 1981]. In works of R.V. Namm, G. Woo and E.M. Vikhtenko [Namm & Vikhtenko, 2011, Vikhtenko et al., 2014] these studies were extended to the infinite-dimensional variational inequalities of mechanics. The regularity of the solution may be arbitrary bad in the neighborhood of the crack edges, and the dual problem may be unsolvable. Therefore the applying a similar scheme for solving the crack problem may be complicated. Despite this problem, it is possible to justify the duality scheme to solve the crack problem, as well as equality of duality for the original and dual problems. 2 Formulation of the Problem Let Ω ⊂ R2 be a bounded domain with Lipschitz boundary Γ and γ ⊂ Ω be a cut (crack) inside of Ω. For simplicity we assume γ = { (x1 , x2 ) ∈ Ω : a < x1 < b, x2 = const} . and suppose that both end points (a, 0) and (b, 0) do not belong to the boundary Γ. Denote Ωγ = Ω\γ̄. Consider the minimizing problem  ∫ ∫  J (v) = 1 2 |∇v| dΩ − f v dΩ − min, 2{ Ωγ Ωγ } (1)  v ∈ K = v ∈ H 1 (Ωγ ) : [v] ≥ 0 on γ, v = 0 on Γ . Here [v] = v + − v − is the jump of v across γ (v + is a function value v on upper crack face, v − is a function value v on lower crack face, marks ± correspond to positive and negative directs of normal vector on cut γ); f ∈ L2 (Ω) is a given function. Problem (1) has a unique solution u ∈ K, which is, simultaneously, a solution of variational inequality [Khludnev, 2010] ∫ ∫ ∇u∇ (v − u) − f (v − u) dΩ ≥ 0 ∀v ∈ K. (2) Ωγ Ωγ Assuming H 2 -regularity of function u it can be shown that the problem (2) is equivalent to the boundary value problem [Khludnev, 2010]   −∆u = f in Ωγ , u=0 on Γ, (3)  [u] ≥ 0, [ux2 ] = 0, ux2 ≤ 0, ux2 [u] = 0 on γ. The question of solvability of problem (1) is investigated in detail in [Khludnev, 2010]. There is a theorem that establishes the existence of the solution of problem (1). Theorem 1. Let K is a convex and closed set in Hilbert space H 1 (Ωγ ), functional J (v) is weakly lower semicontinuous and coercive. Then problem (1) is solvable. Applying Friedrich’s inequality ∫ ∫ 2 v 2 dΩ ≤ C |∇v| dΩ ∀v ∈ H 1 (Ωγ ) , v = 0 on Γ, Ωγ Ωγ it can be proved that the problem (1) has unique solution. 3 Duality Scheme for Solving Model Crack Problem For arbitrary m ∈ L2 (γ) construct the set { } Km = v ∈ H 1 (Ωγ ) : v = 0 on Γ, − [v] ≤ m a.e. on γ . It is easy to show that Km is a convex and closed set in H 1 (Ωγ ). If the function [m] ∈ L2 (γ) \H 1/2 (γ) the corresponding set Km may be empty. On space L2 (γ) define the sensitivity functional { inf J (v) , if Km ̸= ∅, χ (m) = v∈Km +∞, otherwise. 140 Functional χ (m) is a proper convex and coercive functional on L2 (γ), but it’s effective domain domχ = {m ∈ L2 (γ) : χ (m) < +∞} does not coincide with L2 (γ). Notice that domχ is a convex but not closed set. In this case, domχ = L2 (γ). Theorem 2. The functional χ (m) is weakly lower semicontinuous on L2 (γ). The proof of this statement is presented in [Vikhtenko & Namm, 2016]. This property of the sensitivity functional is the basis of some theorems allows constructing and justifying the search algorithms of the saddle points of the modified Lagrangian functional. We define the following functional on the space H 1 (Ωγ ) × L2 (γ) × L2 (γ) [Namm & Vikhtenko, 2011, Vikhtenko et al., 2014]  ∫ ( )  J (v) + 1 2 (l + r m) − l2 dΓ, if − [v] ≤ m a.e. on γ, K (v, l, m) = 2r γ  +∞, otherwise, and modified Lagrangian functional ∫ (( )2 ) 1 + M (v, l) = inf K (v, l, m) = J (v) + (l − r [v]) − l2 dΓ. m∈L2 (γ) 2r γ + Here r > 0 is a constant, (l − r [v]) = max {0, l − r [v] , }. Let us introduce the modified dual functional M (l) = inf M (v, l) = inf inf K (v, l, m) = v∈H 1 (Ωγ ) v∈H 1 (Ωγ ) m∈L2 (γ) { ∫ (( )2 ) } 1 + = inf J (v) + (l − r [v]) −l 2 dΓ . (4) v∈H 1 (Ωγ ) 2r γ Since inf inf K (v, l, m) = inf inf K (v, l, m) , v∈H 1 (Ωγ ) m∈L2 (γ) m∈L2 (γ) v∈H 1 (Ωγ ) then functional M (l) has the another presentation { ∫ ( ) } 1 2 M (l) = inf χ (m) + (l + r m) − l2 dΓ . (5) m∈L2 (γ) 2r γ For an arbitrary l ∈ L2 (γ), we consider the functional ∫ ( ) 1 2 Fl (m) = χ (m) + (l + r m) − l2 dΓ. 2r γ Then the dual functional M (l) has the form M (l) = inf Fl (m) . (6) m∈L2 (γ) Fl (m) is a weakly lower semicontinuous and coercive functional on L2 (γ). Therefore the problem (6) has a solution m (l) for every l ∈ L2 (γ). It is easy to see, that element m (l) is unique. Further, we mention one important property of dual functional [Vikhtenko & Namm, 2016]. Theorem 3. The dual functional M (l) is Gateaux differentiable in L2 (γ) and it’s derivative ∇M (l) satisfies the Lipschitz condition with the constant r−1 ; that is, for all l1 , l2 ∈ L2 (γ), it holds that ∥∇M (l1 ) − ∇M (l2 )∥L2 (γ) ≤ r−1 ∥l1 − l2 ∥L2 (γ) . In the proof of this theorem it shows that the functional subdifferential consists of a single element ∂ (M (l)) = { m (l) }. Therefore ∇M (l) = m (l). Using equations (4) and (5) for the dual functional, we get { } l ∇M (l) = m (l) = max − , − [v] . r 141 Let us consider the dual problem { M (l) − sup, (7) l ∈ L2 (γ) . Since the gradient of the functional M (l) satisfies the Lipschitz condition, the dual problem (7) can be solved by using the gradient method for maximizing a functional lk+1 = lk + r ∇M (lk ) = lk + r m (lk ) , k = 0, 1, 2, . . . (8) Here l0 ∈ L2 (γ) is an arbitrary initial value. Theorem 4. The sequence {lk } constructed by the gradient method (8) satisfies the limit equality lim ∥m (lk )∥L2 (γ) = 0. k→∞ Using gradient method (8), we can construct the following algorithm Uzawa method for solving the problem (1) (i) uk+1 = arg min M (v, lk ) ; v∈H 1 (Ω{ γ) } lk (ii) lk+1 = lk + r max − , − [uk+1 ] . r Justifying the convergence presented method is complicated by the fact that the problem (7) can be unsolvable. The question of solvability of problem (7) is closely connected with the regularity of a solution u of problem (1). It can be proved if u ∈ H 2 (Ωγ ), then −[ux2 ] ∈ H 1/2 (γ) is a solution of problem (7) [Khludnev, 2010]. This assumption looks unnatural in crack problem. Despite this problem, the duality ratio can be proved for initial and dual problem [Vikhtenko & Namm, 2016]. Theorem 5. There is a duality ratio sup M (l) = inf J (v) . l∈L2 (γ) v∈K Note that if the dual problem (7) has a solution, then {lk } is a bounded sequence in L2 (γ) [4]. Together with theorem 4 it means that method (i), (ii) converges according the initial functional J (u), that is lim J (uk ) = inf J (u) . k→∞ v∈K 4 Numerical Implementation of Method Ω = {(x, y) ∈ R2 : 0 ≤ x ≤ 1; 0 ≤ y ≤ 1} and γ = {(x, y) ∈ Ω: a < x < b, y = c}. Approximation of the problem made by means of finite element method. To perform a triangulation of the region by step h, as shown in figure 1. Introduce the notation: Iγ is the set of indices of internal cracks nodes; IΩ is the set of indices of all other nodes of the triangulation. Figure 1: Triangulation of the computational domain 142 For each node of triangulation we define a piecewise affine basis function φk (x, y). Also define functions − φ+ k (x, y) and φk (x, y) for each top and bottom node on the crack faces. Construct an approximate solution as linear combination of basis functions φk (x, y). ∑ ∑ ∑ vh (x, y) = vk φk (x, y) + vk+ φ+ k (x, y) + vk− φ− k (x, y) . k∈IΩ k∈Iγ k∈Iγ Let us substitute function vh (x, y) in a modified Lagrangian functional M (v, l) taking change of variables vi+ = vi− + ti , ti = [vi ] ≥ 0, i ∈ Iγ . Also we perform trapezoid approximation of the integral ∫ (( )2 ) { } + h∑ ( + )2 (l − r[v]) − l dΓ ≈ 2 (li − r ti ) − li . 2 γ 2 i∈Iγ As result we obtain finite-dimensional functional of the following form  ( ) 1 ∑ ∑ ∑ ∑ ( − )( − ) ∑ ∑ − − − M vh , vh− , th , lh = aij vi vj + a+ ij vi + ti vj + t j + aij vi vj + 2 i∈IΩ j∈IΩ i∈Iγ j∈Iγ i∈Iγ j∈Iγ  ∑ ∑ ( − ) ∑ ∑  +2 b+ ij vi vj + tj + 2 b− − ij vi vj −  i∈IΩ j∈Iγ i∈IΩ j∈Iγ ∑ ∑ {( )2 } ( − ) ∑ − − h ∑ + − f i vi − fi+ vi + t i − f i vi + (li − r ti ) − li2 , 2r i∈IΩ i∈Iγ i∈Iγ i∈Iγ where ∫ aij = ∇φi ∇φj dΩ, i, j ∈ IΩ , Ωγ ∫ ∫ a+ ij = ∇φ+ i ∇φj dΩ, + a− ij = ∇φ− − i ∇φj dΩ, i, j ∈ Iγ , Ωγ Ωγ ∫ ∫ b+ ij = ∇φi ∇φ+ j dΩ, b− ij = ∇φi ∇φ− j dΩ, i ∈ IΩ , j ∈ Iγ , Ωγ Ωγ ∫ ∫ ∫ fi = f ∇φi dΩ, fi+ = f ∇φ+ i dΩ, fi− = f ∇φ− i dΩ, i ∈ Iγ . Ωγ Ωγ Ωγ According to the Uzawa method on the iteration number k it is necessary to solve auxiliary problem { ( ) M vh , vh− , th , lhk − min, −ti ≤ 0, i ∈ Iγ . To solve this problem we apply the pointwise relaxation method. Computational experiments were conducted with the following fixed parameters, which are presented in the table 1. Consider three examples differing in way of defining function f in the computational domain Ω. The values of the function and the corresponding solution of the problem (1) is shown in figures 2–4. In example 1 we attempt to raise the portion positioned below the crack. For this purpose, we define a positive value of the function f on the lower face of the crack, and negative on the top face. As seen from the corresponding graph solution in figure 2, the condition [v] = v + − v − ≥ 0 is satisfied on γ. Disconnection of crack faces does not occur; we obtain the solution of the Dirichlet problem. In example 2, the value of the function f is reversed. In this case, we observe opening faces along the entire length of the crack. Corresponding graph solution is presented in figure 3. 143 Table 1: Fixed parameters of computational experiments Parameters Parameter { value } Computational domain Ω = (x, y) ∈ R2 : 0 ≤ x ≤ 1, 0 ≤ y ≤ 1 Positioning of crack in domain γ = {(x, y) : 0.2 < x < 0.8, y = 0.4} Step triangulation h = 2−6 Computational accuracy ε = 10−8 Figure 2: The calculation results. Example 1 Finally, in example 3, the function is defined so that to the upper face of crack the portions are adjacent with a positive and negative value of the function f . As seen from the figure 4 there is a partial disconnection of the crack faces only, in that part of the domain Ω where the value of function f is negative at the upper face. Condition [v] = v + − v − ≥ 0 is still satisfied on γ. Let’s consider the question of the relationship between the number of iterations by Uzawa method with parameter r which is the shift step value by the dual variable. We distinguish two types of iteration: 1. The external iterations is the number of step (ii) of the Uzawa method to achieve the required accuracy. 2. The internal iterations is the number of step (i) of the Uzawa method to achieve the required accuracy. The table 2 shows the calculation results for example 3 for different values of parameter r. Here null vectors chosen as an initial approximation. Table 2: Number of iterations by Uzawa method for example 3 Value r Number of internal iterations Number of external iterations 10 6029 51 102 3208 11 103 2519 5 104 2321 3 Conducted computing experiments have shown that the convergence of the method to the saddle point of the Lagrangian functional is good enough. At that, increasing the value of the parameter r leads to growth convergence rate. For considered example on the grid in steps of h = 2−6 , with an accuracy of ε = 10−8 , it is necessary less than sixty external iterations. 5 Conclusion These results allow us to conclude about the effectiveness of the chosen approach for solving model crack problem, one of the formulations of this is a minimization problem (1). The main advantage of using a modified Lagrangian functional is the shift step value by the dual variable is large positive value in distinct from the classical scheme of duality. Therefore duality method with a modified Lagrangian functional superior to the standard scheme uses classical functional for computing speed. 144 Figure 3: The calculation results. Example 2 Figure 4: The calculation results. Example 3 References [Hlavachek et al., 1988] Hlavačhek, I., Haslinger, Ya., Nechas, I., & Lovišhek, Ya. (1988) Numerical solution of variational inequalities. New York: Springer-Verlag. [Kikuchi & Oden, 1988] Kikuchi, N., & Oden, J. T. (1988) Contact Problems in Elasticity: A study of variational inequalities and finite element methods. Philadelphia: SIAM. [Khludnev, 2010] Khludnev, A. M. (2010) Elasticity problems in non-smooth domains . Moscow: Fizmatlit. [Gol’stein & Tret’yakov, 1989] Gol’stein, E. G. & Tret’yakov, N. V. (1989) Modified Lagrange function: theory and optimization methods. Moscow: Nauka, Fizmatlit. [Grossman & Kaplan, 1981] Grossman, Ch. & Kaplan, A. A. (1981) Nonlinear programming on the basis of unconstrained optimization. Novosibirsk: Nauka. [Namm & Vikhtenko, 2011] Namm, R. V. & Vikhtenko, E. M. (2011) Modified Lagrangian functional for solving the Signorini problem with friction. New-York. J. Advances in Mechanic Research, 1, 435-446. [Vikhtenko et al., 2014] Vikhtenko, E. M. & Woo, G. & Namm, R. V. (2014) Sensitivity functionals in contact problems of elasticity theory. J. Comp. Math. and Math. Phys., 54, 1190-1200. [Vikhtenko & Namm, 2016] Vikhtenko, E. M. & Namm, R. V. (2016) Duality method for solving model crack problem. J. Proceedings of the Institute of Mathematics and Mechanics, 22, 36-45. 145