CONVERGENCE OF A NEUMANN-DIRICHLET ALGORITHM FOR TWO-BODY CONTACT PROBLEMS WITH NON LOCAL COULOMB S FRICTION LAW
|
|
- Jeffery Norman
- 5 years ago
- Views:
Transcription
1 ESAIM: MAN DOI: /man: ESAIM: Mathematical Modelling and Numerical Analysis CONVERGENCE OF A NEUMANN-DIRICHLET ALGORITHM FOR TWO-BODY CONTACT PROBLEMS WITH NON LOCAL COULOMB S FRICTION LAW Guy Bayada 1, Jalila Sabil and Taoufik Sassi 3 Abstract. In this paper, the convergence of a Neumann-Dirichlet algorithm to approximate Coulomb s contact problem between two elastic bodies is proved in a continuous setting. In this algorithm, the natural interface between the two bodies is retained as a decomposition zone. Mathematics Subject Classification. 65N30, 65N55, 65K05. Received December 3, 006. Revised September 1, Introduction Domain decomposition methods are usually characterized by an artificial splitting of a given domain where a partial differential equation has to be solved. Then, suitable iterations are performed in which a sequence of local problems are considered. Each of these local problems being defined in a subdomain with convenient boundary conditions. One of the potential interest of this approach is that the computation of solutions for the local problems required much less time than the solution for the initial problem and can be made simultaneously by using parallel computers. Moreover various discretizations or various kinds of solvers can be used to deal with each of the local problems. Several methods have been reviewed in [9,17,19]. At first glance, multibody contact problems in elasticity appear to be a natural domain of application for domain decomposition methods. Many numerical procedures have been proposed in the mechanical literature. They are based on standard discretization techniques for partial differential equations in combination with a special implementation of the non-linear contact conditions see e.g. [5,6,10 13,0]. However, so far this kind of approach seems to have been neglected in mathematical literature. One of the first mathematical proof of the convergence for contact decomposition algorithm is proposed in [1]. However, the decomposition boundary is artificial and located inside one of bodies. More recently, a so-called Neumann-Dirichlet algorithm for the solution of frictionless Signorini contact problem has been proposed and studied in the discretized setting in [14 16]. This new algorithm consists Keywords and phrases. Domain decomposition methods, contact problems, convergence. 1 INSA-LYON, CNRS UMR 508, UMR 5514 Bâtiment Léonard de Vinci, 1 Av. J. Capelle, 6961 Villeurbanne Cedex, France. guy.bayada@insa-lyon.fr Université Henri Poincaré, UFR STMP, Faculté des Sciences et Techniques, Laboratoire LEMTA B.P. 39, Vandœuvreles-Nancy Cedex, France. Jalila.Sabil@lemta.uhp-nancy.fr 3 Laboratoire de Mathématiques Nicolas Oresme, LMNO CNRS UMR 6139, Université de Caen, Bd. Maréchal Juin, 1403 Caen Cedex, France. Taoufik.Sassi@math.unicaen.fr Article published by EDP Sciences c EDP Sciences, SMAI 008
2 44 G. BAYADA ET AL. to solve in each iteration a linear Neumann problem for one body and a unilateral contact problem for the other by using essentially the contact interface as the boundary data transfer. The convergence of this algorithm in the continuous setting has been proved in [4,8]. This last approach has been extended to the contact problem with Coulomb friction. The corresponding discretized Neumann-Dirichlet algorithm and several numerical results are given in [14,16]. The aim of the present paper is to prove the convergence of the Neumann-Dirichlet algorithm for contact problem with Coulomb friction in the continuous setting. The main difficulty in this work is linked with the boundary conditions at the contact interface. They are highly non-linear both in the normal direction unilateral contact conditions and in the tangential one non-differential Coulomb law. A fixed point relaxed procedure is defined for the stresses on the contact surface. For sufficiently small friction coefficient, the convergence of the Neumann-Dirichlet algorithm is proved for the continuous problem. It is also shown, by some numerical calculations, that an optimal relaxation parameter exists and its value is nearly independent of the friction coefficient, the mesh size and of the Young modulus. The paper is organized as follows: In Sections and 3, we give a statement of the problem and we propose the natural Neumann-Dirichlet algorithm in the strong formulation. The precise variational formulation of the algorithm is given in Section 4. We then prove in Section 5, under some assumptions concerning the data, the convergence of the stresses on the contact surface, which in turn proves the convergence of the whole algorithm. In the last, we present some numerical results asserting the efficiency of our algorithm.. Stating the problem Let us consider two elastic bodies, occupying two bounded domains Ω α, α =1,, of the space R. The boundary Γ α = Ω α is assumed piecewise continuous, and composed of three complementary non empty parts Γ α u,γ α l and Γ α c such that Γ α u Γ α c =. Each body Ω α is fixed on the part Γ α u and subjected to surface traction forces φ α in L Γ α l. The body forces are denoted by f α in L Ω α. In the initial configuration, both bodies have a common contact portion Γ c =Γ 1 c =Γ c. In other words, we consider the case when contact zone cannot grow during the deformation process. Unilateral contact with non local Coulomb s friction can take place along the boundary Γ c. The problem consists in finding the displacement field u =u 1, u where the notation u α stands for u Ω α and the stress tensor field σ =σu 1,σu such that: for α =1, div σu α +f α = 0 in Ω α, σu α n α = φ α on Γ α l, u α = 0 on Γ α u,.1 where the symbol div denotes the divergence operator of a tensor function and is defined as div σ = σij The summation convention of repeated index is adopted. The elastic constitutive law is given by Hooke s law for homogeneous and isotropic solid: x j σ ij u α =A α ijkhe kh u α, eu α = 1 i u α + u α T,. where A α x =a α ijkh x 1 i,j,k,h L Ω α 16 is a fourth-order tensor satisfying the usual symmetry and ellipticity conditions in elasticity, and eu α is the strain tensor.
3 NEUMANN-DIRICHLET ALGORITHM FOR CONTACT 45 We will use the usual notations for the normal and tangential components of the displacement and stress vector on Γ c u α N = uα i nα i, uα T i = u α i uα N nα i σn α = σ ij u α n α i n α j, α σ T i = σ ij u α n α j σn α n α, in which we have denoted by n α the outward normal unit vector to the boundary. On the interface Γ c, the unilateral contact law is described by σ 1 N = σ N = σ N, σ 1 T = σ T = σ T,.3 [u N ] 0, σ N 0, σ N [u N ]=0,.4 where [v N ]=v 1.n 1 + v.n, is the jump across the interface of any function v defined on Ω α. The Coulomb s law with non-local friction is given by σ T k Sσ N σ T <k Sσ N = [u T ]=0.5 σ T = k Sσ N = ν 0 [u T ]= νσ T where kx is the coefficient of friction on the interface Γ c : kx L Γ c, kx 0a.e. regularization operator from the dual of H 1 Γc intol Γ c. It s easy to prove see for example [7] that.5 is equivalent to: { σt k Sσ N k Sσ N [u T ] + σ T [u T ]=0. onγ c ; S is a.6 3. Neumann-Dirichlet algorithm Let 0 <θ<1 be a parameter that will be determined in order to ensure the convergence of the algorithm. Let g0 be a given element on Γ c.fori 1, the sequence of functions u 1 i i 0 the following problems: and u i i 0 is defined by solving I divσu i =f in Ω σu i n = φ on Γ l u i = 0 on Γ u σu i n = gi 1 on Γ c II divσu 1 i =f1 in Ω 1 σu 1 i n1 = φ 1 on Γ 1 l u 1 i = 0 on Γ1 u 3.1 u 1 i.n1 u i.n, σi 1 N 0, σ1 i N [u in ]=0 onγ c σi 1 T k Sσi N σi 1 T [u it ]+k Sσi N [u it ] =0 onγ c III g i := θσ1 i n1 +1 θg i 1 on Γ c. The variational form of this algorithm will be given in 4.9. The proof of its convergence towards the solution of the problem.1.6 will be given in Theorem 5.9.
4 46 G. BAYADA ET AL. 4. The variational formulation In order to study the previous algorithm, let us first define its variational formulation by introducing the following functional spaces: For α =1, : V α = { v α H 1 Ω α, v α =0onΓ α u}, endowed with the norm α = H 1 Ω α. V = V 1 V provided with the norm V = α K = {v V, { [v N ] 0 a.e. on Γ c }. } H 1 Γc = ϕ L Γ c ; v V α ; α =1, ; γv Γc = ϕ, provided with the norm α=1 1. ϕ 1,Γc =inf{ v α ; v V α,γv Γc = ϕ}, where γ is the usual trace operator. As Γ α c Γα u =, thespaceh 1 Γc does not depend on α. Let us denote by V α the dual space of V α and by, α the duality pairing of the two spaces. The dual norm on V α is given by: Let us introduce for all ϕ in H 1 Γ c theset: We will denote by: ψ V α = sup v V α ψ, v α v α V 1 ϕ ={v 1 V 1, v 1.n 1 ϕ a.e on Γ c }. a α u, v = A α ijkl e ijue kl vdx, Ω α au, v = L α v = f α.v dx + Ω α φ α.v dx, Lv = Γ α l a α u, v, α=1 L α v. The bilinear form a α.,. is continuous and coercive since mesγ α u > 0, so aα.,. satisfies: α=1 M α > 0, a α u, v M α u α v α u, v V α, m α > 0, a α v, v m α v α v V α. 4.1 Let so M = M α, α=1 m = min α=1, mα, 4. au, v M u V v V, u, v V, av, v m v V v V. 4.3
5 NEUMANN-DIRICHLET ALGORITHM FOR CONTACT 47 Let us first recall that the variational formulation of the initial problem.1.6 isgivenby: Find u in K such that: P 1 au, v u+ k Sσ N [v T ] [u T ] Lv u v K. Γ c 4.4 which has at least one solution and this solution is unique for a sufficiently small k L Γ c see e.g. [7]. In the following, we will use some lift operators which allow us to build specific functions in Ω α from their values on Γ c. For α =1,, let R α : H 1 Γ D c V α ϕ R α ϕ. 4.5 D { a α R α ϕ, v =0 v V α such that γv =0onΓ D c γr α ϕ=ϕ on Γ D c, whose strong formulation is given by: div σ α =0 inω α, with A α er α ϕ=σα D γr α ϕ=ϕ on Γ D c σ α n α =0 onγ α l γr α ϕ=0 onγα D u. 4.6 Proposition 4.1. i The three norms ϕ 1 and,γc Rα ϕ D α,α=1, are equivalent on H 1 Γ c. ii There exists a constant C>0such that Proof. i Using the definition of R α D,wehave: v V α R β D γv Γ c β C v α α, β =1, and α β. γ R α Dϕ = ϕ, so from the definition of 1 ϕ 1,Γc Rαϕ D Conversely, let us consider the range of H 1 Γ c byr α D W α = R α D H 1 Γc. As for any ϕ in H 1 Γc, problem 4.5 has a unique solution, and as the trace restriction on Γ c of any function in W α is unique, so the application R α is a bijection from H 1 Γ D c inw α. Let us consider its inverse R α 1 D.It is easy to prove that this last bijection is linear, and we can obtain its continuity from 4.7. Moreover, as W α is a closed sub-space of V α which is a Banach space, so W α is also a Banach one. Using the homeomorphism Banach theorem see e.g.[18], we obtain that R α is continuous from H 1 D Γc in W α, i.e. R α ϕ D 1 C ϕ 1 ϕ H 1,Γc Γc. 4.8 Consequently, the three norms ϕ 1 and,γc Rα ϕ,α=1, areequivalent. D α ii Letv be in V α, using inequality 4.8 and the trace theorem, the result is established.
6 48 G. BAYADA ET AL. Remark 4.. For sake of simplicity, we will write in the following R α D γv instead of Rα D γv Γ c. The weak formulation of the algorithm 3.1 isgivenby: Let g0 be a given element in V.Fori 1, let us define the sequence of functions u 1 i i 0 in V 1 and u i i 0 in V by solving the following problems: I Findu i in V such that: a u i, v =L v+ g i 1, v v V Q i IIFindu 1 i in V 1 u i N such that: a 1 u 1 i, w + R1 γu D i u1 i Γ + k Sσ in u i w T [u it ] c L 1 w + R 1 D γu i u 1 i w V III g i, v := θ a 1 u 1 i, R1 D γv+l1 R 1 D γv +1 θ g i 1, v v V. Remark 4.3. Using inequality ii ofproposition4.1, the expression g i, v in 4.9IIIhassense,i.e. g i belongs to V. 5. Convergence The convergence of the algorithm towards the solution of P 1 andthenof.1.6 isgivenbyatwo-step procedure. We prove first, in Proposition 5.1, that the convergence of gi i in V induces the strong convergence of u α i i. Then, in Proposition 5.6, the contraction of an auxiliary operator T θ defined by 5.7 isproved.this yields the convergence of the sequence gi i by way of Theorem 5.9. Proposition 5.1. Let us suppose that there exists 0 < θ min < 1 such that for any θ, θ θ min, the sequence gi i 0 converges in V. Then the sequence u 1 i, u i strongly converges in V to the solution u1, u of P 1. Proof. Let u 1 i, u i andu1 j, u j be solutions of the problems Q iandq j respectively, so from 4.9III we have: a 1 u 1 i u1 j, R1 D γv = g i 1 g j 1, v g i g i 1 θ v V. g j g j 1, v θ Choosing v = R γu1 D i u1 j in the previous inequality and using the fact that R1 γr γu1 D D i u1 j = u1 i u1 j, we get: u 1 i u1 j 1 g 1 i 1 θ min m g j 1 + V g i g i 1 V + g j g j 1 R V γu1 D i u1 j. Using Proposition 4.1i, and the strong convergence of g i i 0 in V, one obtains the convergence of u 1 i i. The strong convergence of u i i can be obtained from the one of g i i 0, by using the coerciveness of a, in 4.9I.
7 NEUMANN-DIRICHLET ALGORITHM FOR CONTACT 49 Finally, we pass to the limit in equations IIandIIIofproblemQ i. The limit of u α i i,α=1, appears to be a solution of the two following problems: { Find u in V such that: P a u, v =L v a 1 u 1, R 1 D γv + L 1 R 1 D γv v V 5.1. Find u 1 in V 1 u N such that: P 3 a 1 u 1, w + R 1 D γu u 1 + k Sσ N w T [u T ] Γ c L 1 w + R 1 Dγu u 1 w V 0. 1 The following Lemma 5. proves that problems P andp 3 implyp 1 so allowing us to finish the proof of Proposition 5.1. Lemma 5.. Assuming that a solution exists for problems P and P 3, then there exists also a solution for problem P Proof. Let us make the following shift for the test function w in P 3: v 1 = w + R 1 D γu. We obtain the following problem P 3 which is equivalent to problem P 3 : Find u 1 in V u 1 N such that: P 3 a 1 u 1, v 1 u 1 + k Sσ N v 1 T u T [u T ] L 1 v 1 u Γ c v 1 V 1u N. Let v =v 1, v ink, asv 1.n 1 +v u.n =[v N ] u N u N on Γ c,sow = v 1 R 1 D γv u V 1u N. Choosing now w resp. v u astestfunctioninp 3 resp.inp, P 1 is obtained by addition. In order to prove the convergence of g n n, let us introduce the operator T defined by: 5. T : V V ψ Tψ 5.4 such that: Tψ,v = a 1 w 1, R 1 γv+l1 R 1 γv v V, D D where w 1 is the solution of the following obstacle problem: Find w 1 in V w 1 N such that: a 1 w 1, v w 1 + k Sσ N w v T w T w 1 T w T 5.5 Γ c L 1 v w 1 v V 1 w N and w is the solution of the Neumann problem: { Find w in V such that: a w, v =L v+ ψ, v v V. 5.6 We also define the operator T θ by: T θ : V V ψ T θ ψ =θt ψ+1 θψ. 5.7
8 50 G. BAYADA ET AL. We shall denote by: C R : the maximum of the constants of continuity for R 1 D and R D. C s : the maximum with respect of α =1, of the constant of continuity of the application Sσ α N which to each w α in V α associates Sσ α N wα, the regularization of σ α N wα inl Γ c. C tr : the maximum with respect of α =1, of the constants of continuity of both following trace applications: Tr α : V α H 1 Γ c, v v T, 5.8 where v T is the tangential component of v. C K = k L Γ c C sc tr, where k is the friction coefficient introduced in.6. Lemma 5.3. For w and w in V,letw 1 and w 1 be solutions of the problems: Find w 1 in V w 1 N such that: a 1 w 1, v w 1 + k Sσ N w v T w T [w T ] 5.9 Γ c L 1 v w 1 v V 1w. N then the following estimate holds: Find w 1 in V 1 w such that: N a 1 w 1, v w 1 + k Sσ N w v T w T [ w T ] 5.10 Γ c L 1 v w 1 v V 1 w, N w 1 w 1 1 MC R + C K + m CK w w m Proof. We choose v = w 1 R 1 D w +R 1 D w in5.9 andv = w 1 R 1 D w +R 1 D w in5.10. By adding both inequalities, we obtain a 1 w 1 w 1, w 1 w 1 a 1 w 1 w 1, R 1 D w R 1 D w k Sσ N w Sσ N w 5.1 [w T ] [ w T ]. Γ c The continuity of the applications Sσ N and Tr α defined in 5.8 induces k Sσ N w Sσ N w [w T ] [ w T ] C K w 1 w 1 1 w w + w w. Γ c 5.13 In 5.1, we use for the left-hand side, the coerciveness of a 1.,., and for the right-hand side, the continuity of the extension R α D and inequality 5.13 to obtain: m w 1 w 1 1 MC R + C K w 1 w 1 1 w w + C K w w, 5.14 hence m w 1 w 1 1 MC R + C K w w MCR + C K m + C K w w 4m Then the required result is obtained.
9 NEUMANN-DIRICHLET ALGORITHM FOR CONTACT 51 Proposition 5.4. The operator T is Lipschitz continuous: ψ, ψ in V Tψ T ψ V C T ψ ψ V, with C T = MC R m MC R + C K + C K m Proof. Let ψ and ψ in V, w 1 resp. w 1 andw resp. w bethesolutionsof5.5 and5.6, so Tψ T ψ, v = a 1 w 1 w 1, R 1 Dγv v V From the continuity of a 1, and Proposition 4.1i, we get: Tψ T ψ V MC R w 1 w Accordingto5.6, one has a w w, v = ψ ψ, v v V. The coerciveness of a, implies w w 1 m ψ ψ V From 5.18, 5.11 and5.19, we get the existence of a constant C T such that: T ψ T ψ V C T ψ ψ V. Let us define the following operator R : R : V V ψ R ψ where R ψ is the unique solution in V, of the problem: a R ψ, v = ψ, v v V. 5.0 Let us consider the following norm in V : ψ = a R ψ, R ψ Proposition and. V are equivalent norms on the space V. Proof. By using v = R ψ astestfunctionin5.0, the proof is the direct consequence of 5.1 andofthe continuity and the coerciveness of a,, so that we get: ψ 1 m ψ V, 5. M m ψ ψ V. 5.3
10 5 G. BAYADA ET AL. Proposition 5.6. There exist constants C i, i =0,..., 5, such that T θ is a contraction for any θ in ]C 0,C 1 [ if k L Γ [0,C c [. ]0,C 0 [ if k L Γ [0,C c 3[. ]0,C 0 [ if k L Γ ]C c 3,C 4 [. ]C 6,C 0 [ if k L Γ ]C c 4,C 5 [. Proof. Let ψ and ψ in V, w 1 resp. w 1 andw resp. w bethesolutionsof5.5 and5.6. From the definition of T θ we have: T θ ψ T θ ψ = θ T ψ T ψ +1 θ ψ ψ +θ1 θa R T ψ T ψ, R ψ ψ. 5.4 To prove that the operator T θ is a contraction, we give, in the first step, an estimate for the term a R T ψ T ψ, R ψ ψ in terms of T ψ T ψ and ψ ψ. We use the definitions of operators T and R and the equivalence between and V. First step. From the uniqueness of the solution of 5.0, we get R ψ ψ =w w so a R T ψ T ψ, R ψ ψ = a R T ψ T ψ, w w, 5.5 using the definition 5.0 ofr,wehave: a R T ψ T ψ, w w = Tψ T ψ, w w. 5.6 The definition 5.4 oft gives Tψ T ψ, w w = a1 w 1 w 1, R 1 D γw w. 5.7 From inequality 5.1, we obtain: a 1 w 1 w 1, R 1 D γw w a 1 w 1 w 1, w 1 w 1 k Sσ N w Sσ N w [w T ] [ w T ]. 5.8 Γ c In the right-hand side, the first term will be estimated by using the coerciveness of a 1.,., inequalities 5.18 and 5.: a 1 w 1 w 1, w 1 w 1 m M CR T ψ T ψ, 5.9 for the second term, using inequalities 5.11 seelem. 5.3 and5.19 in5.13 and using 5.3, we obtain: k Sσ N w Sσ N w [w T ] [ w T ] C K Γ c m MC R + C K + M C K m +1 m ψ ψ Taking inequalities into account in 5.5, we obtain: a R T ψ T ψ, R ψ ψ m M CR T ψ T ψ + C K m MC R + C K + M C K m +1 m ψ ψ. 5.31
11 NEUMANN-DIRICHLET ALGORITHM FOR CONTACT 53 Using 5.31 and5.4, we obtain: T θ ψ T θ ψ θ m θ1 θ M CR T ψ T ψ + 1 θ +θ1 θ C K m MC R + C K + M C K m +1 ψ m ψ. 5.3 Using 5.16, we get: T θ ψ T θ ψ θ m θ1 θ M CR T ψ T ψ + 1 θ +θ1 θ M CK C T m MC R Second step. Later on, it is convenient to introduce: We immediately have from C K m ψ ψ m =min{m, 1}, M = M C T C T = MC R m MC R + C K + C K m Inequality 5.33 becomes: T θ ψ T θ ψ θ m θ1 θ T ψ T M ψ CR + 1 θ +θ1 θ M C KCT m MC R + C K m ψ ψ According to the sign of the coefficient of T ψ T ψ, we distinguish two cases: in the first one, the Lipschitz continuity of the operator T is used, and in the second one, we only have to cancel this term: Case 5.7. For any θ C 0 with C 0 = we have θ m θ1 θ 0, so that Proposition 5.4 implies M CR θ m θ1 θ T ψ T M ψ θ m θ1 θ CR M CR Hence, inequality 5.36 becomes m, 5.37 M CR + m C T M m ψ ψ T θ ψ T θ ψ h 1 θ ψ ψ, 5.39
12 54 G. BAYADA ET AL. with h 1 θ = θ m θ1 θ M CR = aθ + bθ + c, C T M m +1 θ +θ1 θ M m C K CT MC R + C K m in which a = C M T m + b = C T c =1. C R + C T CR C KCT MC R C KCT + C K M MC R m m +1 + C K M m m In the following, we will prove first that a>0andb<0, so that h 1 θ will be less than 1 for any θ in ]0, b a [. As θ belongs to [C 0, 1[ cf. 5.37, we will prove under some conditions for the friction coefficient that ]0, b a [ [C 0, 1[ is not empty. This will allow us to obtain an interval of θ in which T θ is a contraction. From the definition 5.35 of C T,wehaveasC K 0 C T MC R C K m + MC R CK m m and C T MC R CK m m, so C T MC R C K CT m + M C R C K m 3 As m 1cf. 5.34, we get: C T C R MC K CT mc R + M C K m 3 From the above expression of a, wegetthata>0andb<0, more precisely a 1andb. Let C 1 = b a 5.40 As mentioned above, the contraction of T θ is proved as soon as θ ]0,C 1 [ [C 0, 1[ if C 0 <C 1.
13 Let Using the fact that a = C M T m b 1, we have so from 5.41, we have: NEUMANN-DIRICHLET ALGORITHM FOR CONTACT 55 C 0 C 1 = m M C + b R + m C M T m b 1 = C T M + b M CR m = = M C R + m a m + CK CT MC R M + CK 4 C m R m M CR M CR + m a m M CR + M 4 C R m C 3 K + mc 3 K + MC R +1C K M 5.41 CR + m a Gx = M x + mx 3 + MCR +1x x 0, 5.4 C 0 <C 1 m M CR + M CR m 3 GC K < Using the fact that G is a strictly increasing function from [0, + [, the previous inequality 5.43 is equivalent to: [ C K [0, G 1 m 5 + m M 3, CR which is equivalent to: k L Γ c [0,C [, with C = The first result of Proposition 5.6 is then proved. G 1 m 5 M C R + m 3 C s C tr Case 5.8. If θ<c 0, inequality 5.33 becomes: T θ ψ T θ ψ h θ ψ ψ, 5.44 in which: with h θ =1 θ +θ1 θ M = a θ + b θ +1, a =1 M m b = + M m C K CT MC R m C K CT MC R C K CT MC R + C K m + C K m + C K m 5.45
14 56 G. BAYADA ET AL. hx a x Figure 1. Case: a > 0. hx 1 1 a 1+1/a 0 1 x Figure. Case: 1 <a < 0. As a = b 1, we have h θ =a θ 1 θ 1a Case -1: a > 0. In this case, we have from the expression of a, a < 1, so 1 > 1. a Following the representation of h in Figure 1, wehave0<h θ < 1 for any θ in ]0, 1[ [0,C 0 [=]0,C 0 [. Let us characterize the positive sign of a by a condition for the friction coefficient. Using the notation 5.4, we have: a =1 m 3 GC K By simple calculations, we see by using 5.4 that: a > 0 0 C K < G 1 m 3 By taking C 3 = G 1 m 3, the second wording of Proposition 5.6 is obtained. C s C tr As h = 1, we have in the case: a a < 0 two possibilities following the sign of a. Case -: 1 <a < 0. In this case, the representation of h is given in Figure and 0 <h θ < 1, for any θ in ]0, 1[ [0,C 0 [=]0,C 0 [. By using 5.46, the condition 1 <a < 0isequivalenttoG 1 m3 <C K < G 1 m 3. The third wording of Proposition 5.6 is obtained by taking C 4 = G 1 m 3 C s C tr
15 NEUMANN-DIRICHLET ALGORITHM FOR CONTACT 57 hx 1 1 a 0 1+1/a 1 x Figure 3. Case: a < 1. Case -3: a < 1. From the representation of h in Figure 3, wehave0<h θ < 1 for any θ in ]1 + 1, 1[ [0,C a 0[=]1+ 1,C a 0[ if <C a 0. The two conditions: a < 1 and1+ 1 <C a 0 are equivalent to: G 1 m 3 <C K < G 1 C 0 m 3 1 C 0 The last result of Proposition 5.6 is proved by taking C 5 = G 1 C0 m 3 1 C 0 C s C tr C 6 =1+ 1 a Theorem 5.9. Assuming that the relaxation parameter θ and the coefficient of friction k satisfy Proposition 5.6, then the sequence gi i 0 converges in V towards g defined by: g, v = a 1 u 1, R 1 D γv + L1 R 1 D γv v V, where u 1, u is the solution of the problem P 1. Moreover, if θ θ min, the solution u α i,α=1, of the algorithm 4.9 converges to the solution u =u1, u of the problem P 1 defined in 4.4 as i tends to infinity. Proof. From algorithm Q i, it is easy to prove that Tgi 1, v = a 1 u 1 i, R 1 Dγv + L 1 R 1 Dγv v V, so T θ gi 1 =gi. On the other hand, one deduces from problems P andp 3 thatt θ g =g. Hence from Proposition 5.6 we have: gi g = T θ gi 1 T θg Kθ gi 1 g The convergence of gi is obtained by a recursive application of Proposition 5.1 concludes the proof. Remark Assumptions of Theorem 5.9 are fulfilled for the Signorini problem k = 0 in which, only the unilateral contact is taken into account. 6. Numerical results As mentioned in the introduction, a Neumann-Dirichlet algorithm for contact problem with Coulomb friction has already been used in [14,16] in a variety of situations. The aim of this section is to illustrate and validate the previous theoretical approach for a large contact device, especially with respect of the choice of the relaxation parameter. The implementation is established using the finite element code PLAST developed inlaboratoire de Mécanique des Contact et des Solides LaMCoS of the National Institute of Applied Sciences in Lyon..
16 58 G. BAYADA ET AL. Γu Ω Γ c Ω Γ 1 u Figure 4. The physical problem. Based on the Lagrangemultiplier concept [3], PLAST is generally used for solving contact problems by global method with nonmatching grids at the contact zone see []. The Neumann-Dirichlet algorithm described in the previous section was implemented in PLAST. In each step of our algorithm, we solve two local problems: a linear Neumann problem on the first body and a contact problem with given friction on the second body. After describing the physical data, we compare the interface coordinates and the contact stresses calculated by the Neumann-Dirichlet method to the results given by the global method using the PLAST code. Then, several results concerning the influence of the coefficient friction, the Young modulus, the mesh-size on the convergence rate are given and the existence of an optimal relaxation parameter is shown Description of the model For all experiments described below, a two dimensional linear elasticity problem of two elastic bodies Ω α, α =1, initially in contact is considered see Fig. 4. The elastic behavior law is given by Hooke s law for homogeneous isotropic elastic materials: σ α ij u = E α ν α 1 ν α 1 + ν α δ ijε kk u+ E α 1+ν α ε ij u, α =1,, where E α and ν α denote Young modulus and Poisson s ratio respectively for the body Ω α, α =1,, and the notation δ ij stands for Kronecher s symbol. The lower body Ω 1 the slave is fixed on its lower side. The upper body Ω the master is submitted to a vertical displacement of mm on Γ u see Fig. 4. Dimensions of Ω 1 and Ω are respectively mm 0.3mm and 1. mm 0.3 mm. 6.. Preliminary results In the first sequence of computations, we use the following data: for both solids, the Poisson s ratio is ν 1 = ν =0.3. The Young modulus E 1 = MPa for the lower body and E = 10 4 MPa for the upper body are assumed. We use k =0. as friction coefficient and a mesh-size of 190 nodes see Fig. 7. Figure 5 depicts the coordinates of the points of the contact interface after deformation and the values of the normal and tangential constraints obtained by the Neumann-Dirichlet methods. The horizontal axis in Figure 5a is the classical x coordinates with x = 0 at the left side of the contact interface.
17 NEUMANN-DIRICHLET ALGORITHM FOR CONTACT 59 a b Figure 5. Interface coordinates and contact stresses by the Neumann-Dirichlet method. In Figure 5b, the horizontal axis is the curvilinear abscissa for the deformed slave contact interface measured from the same origin as the x coordinate. These figures evidence that the contact zone is represented by three parts; a sliding one located for 0 <x 0.6, the slippery one located for 0.6 <x<1, while the non contact area is located for x>1. To be mentioned that we find exactly the same results by using the classical method without decomposition About the optimum relaxation parameter In the following, we present numerical results which shown the influence of the friction coefficient, the Young modulus and the mesh-size on the number of iterations in the Neumann-Dirichlet algorithm. The stopping criteria is defined by: g n g n 1 / g n < 10 9 or g n g n 1 < We present in Figure 10, the number of iterations as a function of the parameter θ for 4 mesh-size with k, E1, E, and ν fixed k =0., E1 = , E = , ν =0.3. The mesh size varies according the number of nodes from 60 nodes to 48 nodes see Figs It clearly appears from Fig. 10 that the number of iterations is not strongly dependent to the mesh-size. Especially there exists an optimum value of relaxation parameters θ opt =0.6 which is nearly independent of the choice of the mesh-size. The result is important as this will be sufficient to find the optimal parameter by solving small dimension problems first. Varying the value of k, k = 0.1, 0., 0.4 Fig.11 does not modify the optimal value of the relaxation parameter which is always located around 0.6 with E1 E, ν and a mesh-size 190 nodes fixed. Finally, choosing the first body much harder E1 = ,E= or much softer than the second body E1 = ,E= or choosing the same Young modulus E1 = ,E= , we obtain the same kind of results Fig. 1.
18 60 G. BAYADA ET AL. Figure 6. Mesh 1: 60 nodes. Figure 7. Mesh : 190 nodes. Figure 8. Mesh 3: 666 nodes. Figure 9. Mesh 4: 48 nodes Convergence DDD for 4 mesh size: k=0., E1= ,E= nodes 190 nodes 666 nodes 48 nodes 40 Number of iteration Relaxation parameter Figure 10. Number of iterations as a function of the relaxation parameter θ for 4 mesh-size with k =0., E1 = , E = , ν = Conclusion A Neumann-Dirichlet algorithm of a domain decomposed contact problem with friction has been introduced and studied. For small friction coefficient, convergence result is derived for the continuous problem.the optimal parameter has been numerically found to be independent of several physical and numerical parameters. This phenomena needs further theoretical investigation. The extension of such approach to D and 3D contact problem between more than bodies should be investigated.
19 NEUMANN-DIRICHLET ALGORITHM FOR CONTACT Convergence DDD: E1= , E= , 190 nodes k=0.1 k=0. k= Number of iteration Relaxation parameter Figure 11. Number of iterations as a function of the relaxation parameter θ for 3 friction coefficient with E1 = , E = , ν =0.3, 190 nodes Convergence DDD: k=0., 190 nodes E1= ,E= E1= ,E= E1= ,E= Number of iteration Relaxation parameter Figure 1. Number of iterations as a function of the relaxation parameter θ for 3 Young modulus with k =0., ν =0.3, 190 nodes.
20 6 G. BAYADA ET AL. Acknowledgements. Special thanks to L. Baillet from Laboratoire LaMCos, UMR 5514, INSA de Lyon for his valuable help in the numerical part. References [1] P.Alart,M.Barboteu,P.Le Tallec and M.Vidrascu,Méthode de Schwarz additive avec solveur grossier pour problèmes non symétriques. C. R. Acad. Sci. Paris Sér. I Math [] L. Baillet and T. Sassi, Simulations numériques de différentes méthodes d éléments finis pour les problèmes contact avec frottement. C. R. Acad. Sci. Paris Sér. II B [3] L. Baillet and T. Sassi, Mixed finite element method for the Signorini problem with friction. Numer. Methods Partial Differential Equations [4] G. Bayada, J. Sabil and T. Sassi, Algorithme de Neumann-Dirichlet pour des problèmes de contact unilatéral: résultat de convergence. C. R. Math. Acad. Sci. Paris [5] A.B. Chandhary and K.J. Bathe, A solution method for static and dynamic analysis of three-dimensional contact problems with friction. Comput. Struc [6] P.W. Christensen, A. Klarbring, J.S. Pang and N. Strömberg, Formulation and comparison of algorithms for frictional contact problems. Internat. J. Numer. Methods Engrg [7] G. Duvaut and J.-L. Lions, Les inéquations en mécanique et en physique, Travaux et Recherches Mathématiques 1. Dunod, Paris 197. [8] C. Eck and B. Wohlmuth, Convergence of a contact-neumann iteration for the solution of two-body contact problems. Math. Models Methods Appl. Sci [9] C. Farhat and F.X. Roux, Implicit parallel processing in structural mechanics. Computational Mechanics Advances [10] R. Glowinski, J.-L. Lions and R. Trémolières, Numerical analysis of variational inequalities, Studies in Mathematics and its Applications 8. North-Holland Publishing Co., Amsterdam Translated from the French. [11] J. Haslinger, Z. Dostál and R. Kučera, On a splitting type algorithm for the numerical realization of contact problems with Coulomb friction. Comput. Methods Appl. Mech. Engrg [1] N. Kikuchi and J.T. Oden, Contact problems in elasticity: a study of variational inequalities and finite element methods, SIAM Studies in Applied Mathematics 8. Society for Industrial and Applied Mathematics SIAM, Philadelphia, PA [13] R. Kornhuber and R. Krause, Adaptive multigrid methods for Signorini s problem in linear elasticity. Comput. Vis. Sci [14] R.H. Krause, Monotone multigrid methods for Signorini s problem with friction. Ph.D. thesis, University of Berlin, Germany 001. [15] R.H. Krause and B.I. Wohlmuth, Nonconforming domain decomposition techniques for linear elasticity. East-West J. Numer. Math [16] R.H. Krause and B.I. Wohlmuth, A Dirichlet-Neumann type algorithm for contact problems with friction. Comput. Vis. Sci [17] P. Le Tallec, Domain decomposition methods in computational mechanics. Comput. Mech. Adv [18] L. Lusternik and V. Sobolev, Précis d analyse fonctionnelle. MIR, Moscow [19] B.F. Smith, P.E. Bjørstad and W.D. Gropp, Domain decomposition, Parallel multilevel methods for elliptic partial differential equations. Cambridge University Press, Cambridge [0] G. Zavarise and P. Wriggers, A superlinear convergent augmented Lagrangian procedure for contact problems. Engrg. Comput
A domain decomposition algorithm for contact problems with Coulomb s friction
A domain decomposition algorithm for contact problems with Coulomb s friction J. Haslinger 1,R.Kučera 2, and T. Sassi 1 1 Introduction Contact problems of elasticity are used in many fields of science
More informationA Domain Decomposition Algorithm for Contact Problems: Analysis and Implementation
Math. Model. Nat. Phenom. Vol. 4, No. 1, 009, pp. 13-146 A Domain Decomposition Algorithm for Contact Problems: Analysis and Implementation J. Haslinger a, R. Kučera b1 and T. Sassi c a Department of Numerical
More informationAn example of nonuniqueness for the continuous static unilateral contact model with Coulomb friction
An example of nonuniqueness for the continuous static unilateral contact model with Coulomb friction Un exemple de non-unicité pour le modèle continu statique de contact unilatéral avec frottement de Coulomb
More informationThe Mortar Finite Element Method for Contact Problems
The Mortar Finite Element Method for Contact Problems F. Ben Belgacem, P. Hild, P. Laborde Mathématiques pour l Industrie et la Physique, Unité Mixte de Recherche CNRS UPS INSAT (UMR 5640), Université
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationOn Mixed Methods for Signorini Problems
Annals of University of Craiova, Math. Comp. Sci. Ser. Volume 30, 2003, Pages 45 52 ISS: 1223-6934 On Mixed Methods for Signorini Problems Faker Ben Belgacem, Yves Renard, and Leila Slimane Abstract. We
More informationA uniqueness criterion for the Signorini problem with Coulomb friction
A uniqueness criterion for the Signorini problem with Coulomb friction Yves REARD 1 Abstract he purpose of this paper is to study the solutions to the Signorini problem with Coulomb friction (the so-called
More informationON THE ANALYSIS OF A VISCOPLASTIC CONTACT PROBLEM WITH TIME DEPENDENT TRESCA S FRIC- TION LAW
Electron. J. M ath. Phys. Sci. 22, 1, 1,47 71 Electronic Journal of Mathematical and Physical Sciences EJMAPS ISSN: 1538-263X www.ejmaps.org ON THE ANALYSIS OF A VISCOPLASTIC CONTACT PROBLEM WITH TIME
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
More informationPenalty Robin-Robin Domain Decomposition Schemes 2 for Contact Problems of Nonlinear Elastic Bodies 3 UNCORRECTED PROOF
1 Penalty Robin-Robin Domain Decomposition Schemes 2 for Contact Problems of Nonlinear Elastic Bodies 3 Ihor I. Prokopyshyn 1, Ivan I. Dyyak 2, Rostyslav M. Martynyak 1,and 4 Ivan A. Prokopyshyn 2 5 1
More informationOn finite element uniqueness studies for Coulomb s frictional contact model
On finite element uniqueness studies for Coulomb s frictional contact model Patrick HILD We are interested in the finite element approximation of Coulomb frictional unilateral contact problem in linear
More informationA new algorithm for solving 3D contact problems with Coulomb friction
A new algorithm for solving 3D contact problems with Coulomb friction Radek Kučera VŠB-TU Ostrava, Department of Mathematics and Descriptive Geometry, name@email.address Summary. The paper deals with solving
More informationTwo-scale Dirichlet-Neumann preconditioners for boundary refinements
Two-scale Dirichlet-Neumann preconditioners for boundary refinements Patrice Hauret 1 and Patrick Le Tallec 2 1 Graduate Aeronautical Laboratories, MS 25-45, California Institute of Technology Pasadena,
More informationMATHEMATICAL MODELLING OF HUMAN LIMB IN MATLAB
MATHEMATICAL MODELLING OF HUMAN LIMB IN MATLAB J. Daněk Department of Mathematics, University of West Bohemia Abstract The contribution deals with mathematical modelling of human limb and numerical simulation
More informationA DUALITY ALGORITHM FOR THE OBSTACLE PROBLEM
Ann. Acad. Rom. Sci. Ser. Math. Appl. ISSN 2066-6594 Vol. 5, No. -2 / 203 A DUALITY ALGORITHM FOR THE OBSTACLE PROBLEM Diana Merluşcă Abstract We consider the obstacle problem in Sobolev spaces, of order
More informationTWO-SCALE SIMULATION OF MAXWELL S EQUATIONS
ESAIM: PROCEEDINGS, February 27, Vol.6, 2-223 Eric Cancès & Jean-Frédéric Gerbeau, Editors DOI:.5/proc:27 TWO-SCALE SIMULATION OF MAWELL S EQUATIONS Hyam Abboud,Sébastien Jund 2,Stéphanie Salmon 2, Eric
More informationEquilibrium analysis for a mass-conserving model in presence of cavitation
Equilibrium analysis for a mass-conserving model in presence of cavitation Ionel S. Ciuperca CNRS-UMR 58 Université Lyon, MAPLY,, 696 Villeurbanne Cedex, France. e-mail: ciuperca@maply.univ-lyon.fr Mohammed
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationThe mortar element method for quasilinear elliptic boundary value problems
The mortar element method for quasilinear elliptic boundary value problems Leszek Marcinkowski 1 Abstract We consider a discretization of quasilinear elliptic boundary value problems by the mortar version
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationSome improvements of Xfem for cracked domains
Some improvements of Xfem for cracked domains E. Chahine 1, P. Laborde 2, J. Pommier 1, Y. Renard 3 and M. Salaün 4 (1) INSA Toulouse, laboratoire MIP, CNRS UMR 5640, Complexe scientifique de Rangueil,
More informationGlowinski Pironneau method for the 3D ω-ψ equations
280 GUERMOND AND QUARTAPELLE Glowinski Pironneau method for the 3D ω-ψ equations Jean-Luc Guermond and Luigi Quartapelle 1 LIMSI CNRS, Orsay, France, and Dipartimento di Fisica, Politecnico di Milano,
More informationAn optimal control problem governed by implicit evolution quasi-variational inequalities
Annals of the University of Bucharest (mathematical series) 4 (LXII) (213), 157 166 An optimal control problem governed by implicit evolution quasi-variational inequalities Anca Capatina and Claudia Timofte
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationquantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner
Quantitative Justication of Linearization in Nonlinear Hencky Material Problems 1 Weimin Han and Hong-ci Huang 3 Abstract. The classical linear elasticity theory is based on the assumption that the size
More informationA mixed finite element method and solution multiplicity for Coulomb frictional contact
A mixed finite element method and solution multiplicity for Coulomb frictional contact Riad HASSANI 1, Patrick HILD 2, Ioan IONESCU 3, Nour-Dine SAKKI 3 1 Laboratoire de Géophysique Interne et Tectonophysique,
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationFETI domain decomposition method to solution of contact problems with large displacements
FETI domain decomposition method to solution of contact problems with large displacements Vít Vondrák 1, Zdeněk Dostál 1, Jiří Dobiáš 2, and Svatopluk Pták 2 1 Dept. of Appl. Math., Technical University
More informationON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS
Chin. Ann. Math.??B(?), 200?, 1 20 DOI: 10.1007/s11401-007-0001-x ON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS Michel CHIPOT Abstract We present a method allowing to obtain existence of a solution for
More informationAn Elastic Contact Problem with Normal Compliance and Memory Term
Machine Dynamics Research 2012, Vol. 36, No 1, 15 25 Abstract An Elastic Contact Problem with Normal Compliance and Memory Term Mikäel Barboteu, Ahmad Ramadan, Mircea Sofonea Uniersité de Perpignan, France
More informationFinite volume method for nonlinear transmission problems
Finite volume method for nonlinear transmission problems Franck Boyer, Florence Hubert To cite this version: Franck Boyer, Florence Hubert. Finite volume method for nonlinear transmission problems. Proceedings
More informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
More informationIntrinsic finite element modeling of a linear membrane shell problem
RR Intrinsic finite element modeling of a linear membrane shell problem PETER HANSBO AND MATS G. LARSON School of Engineering Jönköping University Research Report No. : ISSN -8 Intrinsic finite element
More informationOPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS
PORTUGALIAE MATHEMATICA Vol. 59 Fasc. 2 2002 Nova Série OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS J. Saint Jean Paulin and H. Zoubairi Abstract: We study a problem of
More informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationExistence and uniqueness of the weak solution for a contact problem
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 (216), 186 199 Research Article Existence and uniqueness of the weak solution for a contact problem Amar Megrous a, Ammar Derbazi b, Mohamed
More informationBETI for acoustic and electromagnetic scattering
BETI for acoustic and electromagnetic scattering O. Steinbach, M. Windisch Institut für Numerische Mathematik Technische Universität Graz Oberwolfach 18. Februar 2010 FWF-Project: Data-sparse Boundary
More informationLERAY LIONS DEGENERATED PROBLEM WITH GENERAL GROWTH CONDITION
2005-Oujda International Conference on Nonlinear Analysis. Electronic Journal of Differential Equations, Conference 14, 2006, pp. 73 81. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationMixed exterior Laplace s problem
Mixed exterior Laplace s problem Chérif Amrouche, Florian Bonzom Laboratoire de mathématiques appliquées, CNRS UMR 5142, Université de Pau et des Pays de l Adour, IPRA, Avenue de l Université, 64000 Pau
More informationON THE ASYMPTOTIC BEHAVIOR OF ELLIPTIC PROBLEMS IN PERIODICALLY PERFORATED DOMAINS WITH MIXED-TYPE BOUNDARY CONDITIONS
Bulletin of the Transilvania University of Braşov Series III: Mathematics, Informatics, Physics, Vol 5(54) 01, Special Issue: Proceedings of the Seventh Congress of Romanian Mathematicians, 73-8, published
More informationSébastien Chaumont a a Institut Élie Cartan, Université Henri Poincaré Nancy I, B. P. 239, Vandoeuvre-lès-Nancy Cedex, France. 1.
A strong comparison result for viscosity solutions to Hamilton-Jacobi-Bellman equations with Dirichlet condition on a non-smooth boundary and application to parabolic problems Sébastien Chaumont a a Institut
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationAdvection Diffusion Problems with Pure Advection Approximation in Subregions
Advection Diffusion Problems with Pure Advection Approximation in Subregions M. J. Gander, L. Halpern 2, C. Japhet 2, and V. Martin 3 Université de Genève, 2-4 rue du Lièvre, CP 64, CH-2 Genève. Switzerland.
More informationGeneralized Korn s inequality and conformal Killing vectors
Generalized Korn s inequality and conformal Killing vectors arxiv:gr-qc/0505022v1 4 May 2005 Sergio Dain Max-Planck-Institut für Gravitationsphysik Am Mühlenberg 1 14476 Golm Germany February 4, 2008 Abstract
More informationExistence and Uniqueness of the Weak Solution for a Contact Problem
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. x (215), 1 15 Research Article Existence and Uniqueness of the Weak Solution for a Contact Problem Amar Megrous a, Ammar Derbazi b, Mohamed Dalah
More informationPhD dissertation defense
Isogeometric mortar methods with applications in contact mechanics PhD dissertation defense Brivadis Ericka Supervisor: Annalisa Buffa Doctor of Philosophy in Computational Mechanics and Advanced Materials,
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 207 The elastic-plastic torsion problem: a posteriori error estimates for approximate solutions
More informationON GAP FUNCTIONS OF VARIATIONAL INEQUALITY IN A BANACH SPACE. Sangho Kum and Gue Myung Lee. 1. Introduction
J. Korean Math. Soc. 38 (2001), No. 3, pp. 683 695 ON GAP FUNCTIONS OF VARIATIONAL INEQUALITY IN A BANACH SPACE Sangho Kum and Gue Myung Lee Abstract. In this paper we are concerned with theoretical properties
More informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationLINEARIZED ELASTICITY AND KORN S INEQUALITY REVISITED
Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company LINEARIZED ELASTICITY AND KORN S INEQUALITY REVISITED PHILIPPE G. CIARLET Department of Mathematics, City University
More informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
More informationLINEAR FLOW IN POROUS MEDIA WITH DOUBLE PERIODICITY
PORTUGALIAE MATHEMATICA Vol. 56 Fasc. 2 1999 LINEAR FLOW IN POROUS MEDIA WITH DOUBLE PERIODICITY R. Bunoiu and J. Saint Jean Paulin Abstract: We study the classical steady Stokes equations with homogeneous
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationPseudo-monotonicity and degenerate elliptic operators of second order
2002-Fez conference on Partial Differential Equations, Electronic Journal of Differential Equations, Conference 09, 2002, pp 9 24. http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu
More informationAsymptotic expansions and domain decomposition
Asymptotic expansions and domain decomposition Giuseppe Geymonat, Sofiane Hendili, Françoise Krasucki, Michèle Serpilli, Marina Vidrascu To cite this version: Giuseppe Geymonat, Sofiane Hendili, Françoise
More informationShort title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic
Short title: Total FETI Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ-70833 Ostrava, Czech Republic mail: zdenek.dostal@vsb.cz fax +420 596 919 597 phone
More information1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.
1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the
More informationarxiv: v1 [math.na] 13 Sep 2016
THE PENALTY FREE NITSCHE METHOD AND NONCONFORMING FINITE ELEMENTS FOR THE SIGNORINI PROBLEM ERIK BURMAN, PETER HANSBO, AND MATS G. LARSON arxiv:1609.03745v1 [math.na] 13 Sep 2016 Abstract. We design and
More informationWeak Convergence Methods for Energy Minimization
Weak Convergence Methods for Energy Minimization Bo Li Department of Mathematics University of California, San Diego E-mail: bli@math.ucsd.edu June 3, 2007 Introduction This compact set of notes present
More informationAn error estimate for the Signorini problem with Coulomb friction approximated by finite elements
An error estimate for the Signorini problem with Coulomb friction approximated by finite elements Patrick Hild 1, Yves Renard 2 Abstract he present paper is concerned with the unilateral contact model
More informationMultigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids
Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Long Chen 1, Ricardo H. Nochetto 2, and Chen-Song Zhang 3 1 Department of Mathematics, University of California at Irvine. chenlong@math.uci.edu
More informationA FINITE ELEMENT FORMULATION FOR NONLINEAR 3D CONTACT PROBLEMS
A FINITE ELEMENT FORMULATION FOR NONLINEAR 3D CONTACT PROBLEMS Federico J. Cavalieri, Alberto Cardona, Víctor D. Fachinotti and José Risso Centro Internacional de Métodos Computacionales en Ingeniería
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationIntrinsic finite element modeling of a linear membrane shell problem
arxiv:3.39v [math.na] 5 Mar Intrinsic finite element modeling of a linear membrane shell problem Peter Hansbo Mats G. Larson Abstract A Galerkin finite element method for the membrane elasticity problem
More informationDirichlet-Neumann and Neumann-Neumann Methods
Dirichlet-Neumann and Neumann-Neumann Methods Felix Kwok Hong Kong Baptist University Introductory Domain Decomposition Short Course DD25, Memorial University of Newfoundland July 22, 2018 Outline Methods
More informationTransactions on Modelling and Simulation vol 6, 1993 WIT Press, ISSN X
Auxiliary principle technique for a class of nonlinear variational inequalities M.A. Noor, E.A. Al-Said Mathematics Department, College of Science, King Saud University, PO Box 2455, Riyadh 11451, Saudia
More informationA SOLUTION METHOD OF THE SIGNORINI PROBLEM
2005-Oujda International Conference on Nonlinear Analysis. Electronic Journal of Differential Equations, Conference 14, 2006, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
More informationDRIVING FORCE IN SIMULATION OF PHASE TRANSITION FRONT PROPAGATION
Chapter 1 DRIVING FORCE IN SIMULATION OF PHASE TRANSITION FRONT PROPAGATION A. Berezovski Institute of Cybernetics at Tallinn Technical University, Centre for Nonlinear Studies, Akadeemia tee 21, 12618
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationEXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS WITH UNBOUNDED POTENTIAL. 1. Introduction In this article, we consider the Kirchhoff type problem
Electronic Journal of Differential Equations, Vol. 207 (207), No. 84, pp. 2. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationObstacle problems and isotonicity
Obstacle problems and isotonicity Thomas I. Seidman Revised version for NA-TMA: NA-D-06-00007R1+ [June 6, 2006] Abstract For variational inequalities of an abstract obstacle type, a comparison principle
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationNonlinear Analysis 72 (2010) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 72 (2010) 4298 4303 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na Local C 1 ()-minimizers versus local W 1,p ()-minimizers
More informationEXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM
1 EXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM V. BERTI and M. FABRIZIO Dipartimento di Matematica, Università degli Studi di Bologna, P.zza di Porta S. Donato 5, I-4126, Bologna,
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied athematics http://jipam.vu.edu.au/ Volume 4, Issue 5, Article 98, 2003 ASYPTOTIC BEHAVIOUR OF SOE EQUATIONS IN ORLICZ SPACES D. ESKINE AND A. ELAHI DÉPARTEENT
More informationL p Theory for the div-curl System
Int. Journal of Math. Analysis, Vol. 8, 2014, no. 6, 259-271 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.4112 L p Theory for the div-curl System Junichi Aramaki Division of Science,
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Computational Methods in Elasticity - Michel SALAUN
COMPUAIONAL MEHODS IN ELASICIY Michel SALAUN Départment de Mathématiques et Informatique, ENSICA, oulouse, France Keywords Elastodynamics, Energy minimization, Finite element method, Kirchhoff- Love model,
More informationJournal of Computational and Applied Mathematics. Determination of a material constant in the impedance boundary condition for electromagnetic fields
Journal of Computational and Applied Mathematics 234 (2010) 2062 2068 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationDerivation and Analysis of Piecewise Constant Conservative Approximation for Anisotropic Diffusion Problems
Derivation Analysis of Piecewise Constant Conservative Approximation for Anisotropic Diffusion Problems A. Agouzal, Naïma Debit o cite this version: A. Agouzal, Naïma Debit. Derivation Analysis of Piecewise
More informationON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR TERM
Internat. J. Math. & Math. Sci. Vol. 22, No. 3 (999 587 595 S 6-72 9922587-2 Electronic Publishing House ON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR
More informationEXISTENCE OF SOLUTIONS FOR A RESONANT PROBLEM UNDER LANDESMAN-LAZER CONDITIONS
Electronic Journal of Differential Equations, Vol. 2008(2008), No. 98, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) EXISTENCE
More informationInstabilities and Dynamic Rupture in a Frictional Interface
Instabilities and Dynamic Rupture in a Frictional Interface Laurent BAILLET LGIT (Laboratoire de Géophysique Interne et Tectonophysique) Grenoble France laurent.baillet@ujf-grenoble.fr http://www-lgit.obs.ujf-grenoble.fr/users/lbaillet/
More informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
More informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
More informationNumerical methods for Mean Field Games: additional material
Numerical methods for Mean Field Games: additional material Y. Achdou (LJLL, Université Paris-Diderot) July, 2017 Luminy Convergence results Outline 1 Convergence results 2 Variational MFGs 3 A numerical
More informationA mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N)
wwwijmercom Vol2, Issue1, Jan-Feb 2012 pp-464-472 ISSN: 2249-6645 A mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N) Jaouad El-Mekkaoui 1, Abdeslam Elakkad
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationContinuous Sets and Non-Attaining Functionals in Reflexive Banach Spaces
Laboratoire d Arithmétique, Calcul formel et d Optimisation UMR CNRS 6090 Continuous Sets and Non-Attaining Functionals in Reflexive Banach Spaces Emil Ernst Michel Théra Rapport de recherche n 2004-04
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More information