quantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner
|
|
- Mary Fletcher
- 6 years ago
- Views:
Transcription
1 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 of the displacement gradient tensor is small. In many engineering applications, although the basic assumption is violated due to, e.g., corner singularities, the linear elasticity system is still used to approximate the reality. In this paper, we give a- posteriori error analysis for the eect of the linearization on solutions of nonlinear Hencky materials. A-posteriori error estimates for the eect of idealizations on solutions of more realistic materials will be given in forthcoming papers. 1 Introduction The classical linear elasticity theory is based on the assumption that the size of the displacement gradient tensor is small. The linear Hooke's law in the linear elasticity theory characterizes a large class of elastic materials occupying smooth domains subject to moderate size body forces and boundary tractions. For engineering applications, however, the ideal assumption is not always satised. For instance, consider an elasticity problem on a domain with reentrant corners. The solution of the linear elasticity system exhibits singularities; in a neighborhood of the corners, the strain and stress tensors computed from the solution become arbitrarily large in sizes. Thus, the basic assumption of the linear elasticity theory is grossly violated around the corners. On the other hand, it is a usual practice among engineers to adopt the linear elasticity model, as long as it can be reasonably expected that the violation of the basic assumption is limited to a small region of the domain. Qualitatively, when the input data are small, the expectation is guaranteed by the continuous dependence of the deformation of an elastic material on input data. Obviously, it is more useful to have 1 The work was supported by Research Grants Council of Hong Kong UGC. Department of Mathematics, University of Iowa, Iowa City, IA 54, U.S.A. 3 Department of Mathematics, Hong Kong Baptist University, Kowloon, Hong Kong. 1
2 quantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner domain. The purpose of the paper is to provide quantitative error analysis for the eect of linearization on solutions of nonlinear Hencky materials. We notice that in practice there are material models more realistic than the nonlinear Hencky materials considered here. However, we may view the analysis given in this paper as a preliminary step towards a general a-posteriori error analysis for linearization of more realistic materials. A typical result has the feature that we solve the linear elasticity system and use its solution u to bound the error u? u, with u the solution of a nonlinear Hencky material problem. We emphasize here that to obtain an error bound, all we need to know is the solution of the linear elasticity problem. We give a brief description of nonlinear Hencky materials in the next section. We indicate a procedure to derive a-posteriori error estimates in Section 3. The derivation of a-posteriori error estimates is based on the duality theory in convex analysis. For the sake of completeness, we include some standard results of the duality theory. To have problems to make specic computations, we introduce in Section 4 a family of nonlinear Hencky materials, and apply the theory presented earlier to give an error bound in terms of the solution of the linear elasticity problem for the error caused by the linearization of a nonlinear Hencky material. Finally, we present some numerical examples. We will extend the results of the paper to more realistic materials and their approximations and idealizations in forthcoming papers. The idea of applying dual variational principles has been used to derive quantitative error estimates for the eects of idealizations on solutions in linear torsion problems ([5]), in modeling a Bingham laminar stationary ow uid ([6]) and in linearization of a nonlinear elliptic boundary value problem ([8]). It is also proved useful to obtain computable error estimates for some regularization procedures in numerical analysis ([4], [7]). In [9], quantitative error estimates are given for coecient idealization in linear elliptic problems. Nonlinear Hencky materials For a detailed description of nonlinear Hencky materials, one is referred to [1] or [1]. Let R N be a bounded domain. For a displacement eld v :! R N, we use the notation v = (v 1 ; : : : ; v N ) and introduce the innitesimal strain tensor "(v) :! R N N with components " ij (v) = i j ; 1 i; j i
3 Denote a stress tensor with components ij, 1 i; j N. Deviators of the stress and strain tensors are dened by ij = ij? 1 N tr ij; " ij = " ij? 1 N tr" ij; 1 i; j N where tr = ii, tr" = " ii. Throughout the paper, we adopt the summation convention over repeated indices. For the dot product of two tensors, we use the notation " = " ij ij and " = " ". For a linear isotropic, homogeneous material, we have the stored energy function M (") = 1 (tr") + " = 1 k (tr") + " (.1) where k = + =N, and are the Lame coecients. A nonlinear Hencky material is characterized by a stored energy function of the form M(") = 1 k (tr") + (" ) (.) for a non-negative, continuously dierentiable function satisfying The constitutive law = M (") from (.) is then () = ; () > : (.3) = k tr" I + (" ) " or, equivalently = (" ) ": For the special case when () =, we obtain a linear material, and the constitutive law becomes = k tr" I + " or, equivalently = ": Now we assume a nonlinear Hencky material, initially occupying the domain, is subject to the action of a body force of density f, a traction g on part of the We assume on another part of the boundary? the displacement eld g 1 is given. We may specify other boundary conditions, e.g., normal or tangential components of tractions or displacements. One can easily modify the proofs of the results presented in this paper to derive a-posteriori error estimates for the cases of other boundary conditions. 3
4 Let us assume f (L ()) N ; g 1 (H 1 ()) N ; g (L ( )) N and We denote the Sobolev spaces meas(? 1 ) > : (.4) V = (H 1 ()) N ; V?1 = fv V j v = on? 1 g Then the displacement eld u :! R N total potential energy: of the nonlinear Hencky material minimizes the where the total potential energy u? g 1 V?1 ; E(u) = inffe(v) j v? g 1 V?1 g (.5) E(v) = [M("(v))? fv] dx? We now study the minimization problem (.6). We need g v ds: (.6) Lemma.1 If is a Lipschitz-continuous domain or a crack domain, then E(v)! 1 as kvk V! 1, v? g 1 V?1. Proof. We will use Korn's inequality " (v) dx c kvk V ; 8 v V?1 : (.7) Usually, Korn's inequality is proved for a Lipschitz-continuous domain ([1]). For a domain with cracks, we decompose into 1 and in such a way that both 1 and are Lipschitz-continuous, meas(@ 1 \? 1 ) > and meas(@ \? 1 ) >. Then for a v V?1, Korn's inequality holds for both vj 1 and vj. Thus, Korn's inequality also holds for v. For v V, v? g 1 V?1, using the assumption (.3) and Korn's inequality (.7) on v? g 1, we have E(v) = " k " k (tr"(v)) + (" (v))? fv (tr"(v)) + " (v)? fv min(; 1) = min(; 1) " (v) dx? # # dx? dx? fv dx? g v ds g v ds g v ds [" (v? g 1 ) + "(v) "(g 1 )? " (g 1 )] dx? fv dx? g v ds min(; 1) c kv? g 1 k V? c (kvk V + 1): 4
5 So, E(v)! 1 as kvk V! 1, v? g 1 V?1. Theorem. Let be a Lipschitz-continuous domain or a crack domain. Assume (.3) and (.4). If we further assume () is second-order continuously dierentiable and satises, for some >, () + () ; 8 R (.8) then the minimization problem (.5) has a unique solution u which solves the nonlinear elasticity system [k tr"(u) tr"(v) + (" (u)) "(u) "(v)] dx = 8 v V?1 : fv dx + g v ds (.9) Proof. First we prove the stored energy function M(") dened in (.) is strictly convex. We have M (") = k (tr) + 4 (" ) (" ) + (" ) : (.1) If (" ), then M (") k (tr) + (" ) k (tr) + min(; 1) min(; 1) : If (" ) <, we use the assumption (.8) to obtain M (") k (tr) + 4 (" ) " + (" ) = k (tr) + ( (" ) + (" ) " ) k (tr) + min(; 1) min(; 1) : Thus, M(") is strictly convex. Hence the total potential energy function E(v) is strictly convex. By the assumptions on the function, E(v) is continuous. By Lemma.1, E(v) is coercive. Therefore, the minimization problem (.5) has a unique solution u V with u? g 1 V?1 ([3]). It is easy to verify the solution u satises the equation (.1). In practice, because of its simplicity, linear elasticity theory is preferred as long as the solution of the linear elasticity system can provide suciently accurate information of physical interest. As mentioned in the Introduction, if the material occupies a corner domain, and if the region where the basic assumption of the linear elasticity theory is violated is limited to small neighborhoods of corners, then for a large part of the domain, the governing function 5
6 is close to the identity function (assuming there the material is isotropic and homogeneous). So it is intuitively reasonable to use the linear elasticity theory. The solution u of the linear elasticity theory is the minimizer of the following potential energy function E (v) = [M ("(v))? fv] dx? The solution can also be equivalently characterized by the relation [k tr"(u ) tr"(v) + "(u ) "(v)] dx = 8 v V?1 : g v ds; v? g 1 V?1 : (.11) fv dx + g v ds (.1) The main purpose of the paper is to derive computable bounds in terms of u error u? u. on the 3 A-posteriori error estimates for the linearization First, we recall a standard result from the duality theory in convex analysis. One is referred to [3] for a complete exposition of the duality theory. Let V, Q be two normed spaces, V, Q their dual spaces. Assume there exists a linear continuous operator L(V; Q), with transpose L(Q ; V ). Let J be a function mapping V Q into R the extended real line. Consider the minimization problem: Dene the conjugate function of J by: J (v ; q ) = We have the following standard result. inf J(v; v): (3.1) vv sup [< v; v > + < q; q >?J(v; q)] : (3.) vv;qq Theorem 3.1 Assume: (1) V is a reexive Banach space, Q a normed space. () J : V Q! R is a proper, lower semi-continuous, strictly convex function. (3) 9 u V, such that J(u ; u ) < 1 and q 7! J(u ; q) is continuous at u. (4) J(v; v)! +1, as kvk! 1; v V. Then problem (3.1) has a unique solution u V. And, for any v V with J(v; v) < 1, J(v; v)? J(u; u) J(v; v) + J ( q ;?q ); 8q Q : (3.3) 6
7 In [8], we discussed the question of how to derive a-posteriori error estimates from the above theorem. For the paper to be self-contained, we give a similar discussion below. In our applications later, the problem (3.1) will be the boundary value problem of a nonlinear elliptic system for a nonlinear Hencky material, while v = u in (3.3) will be the solution of a related linear elasticity problem. An application of Theorem 3.1 then allows us to get a computable estimate for the error u? u. To do this, we need: Step 1. Find a suitable lower bound for J(u ; u )?J(u; u), which measures the dierence between u and u. Usually, this lower bound will be a quantity depending on ku? uk in some norm. Step. Take an appropriate q so that the estimate (3.3) is as accurate as possible. If q is chosen to be a solution of the dual variational problem, (3.3) becomes an equality. But to get a solution of the dual variational problem is even more dicult than to solve the problem (3.1). Therefore, it is desirable to have a strategy on choosing a q, which can be obtained easily, and which produces a reasonable error bound. Often, we need to calculate the conjugate function for a function dened by an integral of the form: G(v) = g (x; v(x)) dx: Before stating a theorem on how to calculate its conjugate function, we introduce the following notion: Denition 3. Let be an open set of R n, g : R l! R. g is said to be a Caratheodory function, if: 8 R l ; x 7! g(x; ) is a measurable function; for a.e. x ; 7! g(x; ) is a continuous function. Let m i (1; 1), i = 1; ; l. For each i, denote L m i () = v measurable jvj m i dx < 1 : We have the following theorem, a proof of which can be found in [3]. Theorem 3.3 Let g : R l! R be a Caratheodory function. For any v V = l i=1l m i (), dene: G(v) = g (x; v(x)) dx: Then for the conjugate function of G, G (v ) = g (x; v (x)) dx; 8 v V 7
8 where g (x; y) = sup R l [ y? g(x; ) ] : To apply Theorem 3.1, let us take V = (H 1 ()) N ; V = (H 1 () ) N Q = Q = L () (L ()) N N L ( ) v = (tr"(v); "(v); vj? ) F (v) = 8 < : G(v) =? fv dx if v? g 1 V?1 +1 otherwise M("(v)) dx? J(v; v) = F (v) + G(v) g v ds We notice that when v? g 1 V?1, J(v; v) = E(v), the total potential energy function. The minimization problem (.5) can be equivalently written as u V; J(u; u) = inffj(v; v) j v V g: (3.4) From the corresponding properties of the potential energy function E(v) proved in the last section, it is easy to verify that J(v; v) satises all the assumptions of Theorem 3.1. Therefore, we may apply the inequality (3.3) with v = u, the solution of the linear elasticity problem (.1), to obtain an a-posteriori bound on the error u? u. Let us compute the conjugate functions. For q Q, we use the notation q = (q 1 ; q ; q 3 ) with q 1 L (), q (L ()) N N and q 3 L ( ). We use a similar notation for q Q. We have J ( q ;?q ) = F ( q ) + G (?q ) F ( q ) = supf< v; q >?F (v)g vv = supf< v; q >?F (v)g vv = sup = v?g 1 V?1 8 >< >: [q 1tr"(v) + q "(v) + fv] dx + q 3v ds [q1tr"(g 1 ) + q "(g 1 ) + fg 1 ] dx + q3g 1 ds if [q1tr"(v) + q "(v) + fv] dx + q3v ds = ; 8 v V?1 ; +1; otherwise: 8
9 To compute the conjugate function of G, we will use Theorem 3.3. G (?q ) = supf< q;?q >?G(q)g qq = sup qq ( = sup q 1 L () = = " + sup q (L ()) NN sup t 1 R?q 1q 1? q q? k jq 1j? (q )?q 1q 1? k jq 1j!?q 1t 1? k t 1 dx # dx + (g? q 3)q 3 ds h?q q? (q) i dx + sup (g? q3)q 3 ds q 3 L ( )! dx + + sup(g? q3)t 3 ds t 3 R " jq 1j k + q (q)? ( (q) ) +1; otherwise 8 >< >: where t =? (q ) solves h sup?q t? (t ) i dx t R NN # dx; if q 3 = g on (t )t =? 1 q : (3.5) Before giving an estimate for the error u?u, we analyze the system (3.5) in detail. Taking inner products of both sides of the system (3.5) with themselves, we obtain an equation for the unknown scale = t : ( ()) = 1 4 q (3.6) Let f() = ( ()). We have f () = () () + ( ()) = ()( () + () ) > by assumptions (.3) and (.8). Thus, f is a strictly increasing continuous function. Since f() =, f()! 1 as! 1, we conclude that f : R +! R + is a one-to-one mapping. Therefore, (3.6) is uniquely solvable. Having the solution of the nonlinear scalar equation (3.6), we obtain the solution of the system (3.5) from t =? 1 q () : ) 9
10 Now we apply Theorem 3.1. For any q 1 L (), q (L ()) N N such that we have the estimate [q 1tr"(v) + q "(v) + fv] dx + J(u ; u )? J(u; u) + " k Taking v = u? g 1 in (.1), we obtain = + (tr"(u )) + ("(u ) )? fu g v ds = ; 8 v V?1 (3.7) # dx? g u ds [q1tr"(g 1 ) + q "(g 1 ) + fg 1 ] dx + g g 1 ds 1 k jq 1j + q (q)? ( (q) ) dx (3.8) [ktr"(u )tr"(u? g 1 ) + "(u ) "(u? g 1 )] dx f(u? g 1 ) dx + Using the above relation in (3.8), we then have g (u? g 1 ) ds: J(u ; u )? J(u; u) + "? k (tr"(u )) + jq 1j k + (q 1 + ktr"(u ))tr"(g 1 ) + (q + "(u )) "(g 1 ) dx h ("(u ) ) "(u ) + q (q)? ( (q) ) i dx (3.9) # The inequality (3.9) holds for any auxiliary functions q 1 L (), q (L ()) N N satisfying the constraint (3.7). In particular, from (.1), q 1 =?k tr"(u ); q =? "(u ) (3.1) is an admissible pair of auxiliary functions. With the choice (3.1), we derive the following inequality from (3.9), J(u ; u )? J(u; u) h ("(u ) )? ( (? "(u )) ) i dx ("(u ) + (? "(u ))) "(u ) dx (3.11) Next, we derive a lower bound for the dierence of the potential energies. We have J(u ; u )? J(u; u) 1
11 = J(u ; u )? J(u; u)? < J (u; u); (u? u; u? u) > = " k (tr"(u ))? k (tr"(u))? k tr"(u )tr"(u? u) dx h + ("(u ) )? ("(u) ) ("(u ) ) "(u) "(u? u) i dx = k ktr"(u? u)k L () + ( 1 d dt (("(u) + t"(u? u)) ) dt ("(u ) )"(u) "(u? u) # ) dx where = k ktr"(u? u)k L () 1 h + (("(u) + t"(u? u)) ) ("(u) + t"(u? u)) "(u? u)? ("(u) ) "(u) "(u? u) i dt dx = k ktr"(u? u)k L () h(; t; u; u ) d dt dx h(; t; u; u ) h ("(u) + t "(u? u)) ("(u) + t "(u? u)) "(u? u) = t ("(u) + t "(u? u)) [("(u) + t "(u? u)) "(u? u)] + t ("(u) + t "(u? u)) "(u? u) If (("(u) + t "(u? u)) ), then h(; t; u; u ) t ("(u) + t "(u? u)) "(u? u) t "(u? u) by using the assumption (.3). If (("(u) + t "(u? u)) ) <, then h(; t; u; u ) t ("(u) + t "(u? u)) ("(u) + t "(u? u)) "(u? u) + t ("(u) + t "(u? u)) "(u? u) t "(u? u) by using the assumption (.8). So h(; t; u; u ) t min(; ) "(u? u) 11
12 and J(u ; u )? J(u; u) k ktr"(u? u)k L () + min(; ) k"(u? u)k L () (3.1) In conclusion, we have proved a-posteriori bounds for the error between the displacement eld of a nonlinear Hencky material and that of a linear elastic material. Theorem 3.4 Under the assumptions (.3), (.4) and (.8), we have, for the error in the total potential energy E(u )? E(u) + "? k (tr"(u )) + jq 1j k + (q 1 + ktr"(u ))tr"(g 1 ) + (q + "(u )) "(g 1 ) dx h ("(u ) ) "(u ) + q (q)? ( (q) ) i dx (3.13) and for the error in energy norm k ktr"(u? u)k L () + min(; ) k"(u? u)k L () (3.14) + "? k (tr"(u )) + jq 1j k + (q 1 + ktr"(u ))tr"(g 1 ) + (q + "(u )) "(g 1 ) dx h ("(u ) ) "(u ) + q (q)? ( (q) ) i dx (3.15) # # for any auxiliary functions q 1 L (), q (L ()) N N particular, with the choice (3.1), we have subject to the constraint (3.7). In and E(u )? E(u) h ("(u ) )? ( (? "(u )) ) i dx ("(u ) + (? "(u ))) "(u ) dx (3.16) k ktr"(u? u)k L () + min(; ) k"(u? u)k L () (3.17) h ("(u ) )? ( (? "(u )) ) i dx ("(u ) + (? "(u ))) "(u ) dx (3.18) 1
13 4 Numerical examples To make specic computations, we introduce a family of nonlinear Hencky materials whose governing function in the stored energy function M is of the form () = 8 >< >: ; if + 1? 1?? (1? )(1? ) ; if > (4.1) with parameters (; 1), 6= ((:5? )=(1? ); 1). A material with the governing function (4.1) is linear on the region where the size of the strain deviator tensor is small (bounded by p ). On the region where the size of the strain deviator tensor is large, the slope of the governing function is close to (; 1). When approaches, the material on the nonlinear elasticity region becomes softer. The parameter is introduced to ensure C 1. We have () = 8 < : 1; if + (1? ) 1??1 ; if > (4.) So the assumption (.3) is satised. With the choice (4.1), () is only piecewisely twice dierentiable, so it does not satisfy the assumption (.8). Nevertheless, we will show all the results presented in Sections and 3 still hold. First, let us verify that the corresponding stored energy function is strictly convex, which implies the conclusions of Theorem.. If " <, then from (.1) and (4.1), M (") = k (tr) + = (tr) + : If " >, again from (.1) and (4.1), M (") = k (tr)? 4 (1? ) (1? ) 1? (" )? (" ) + h + (1? ) 1? (" )?1i k (tr)? 4 (1? ) (1? ) 1? (" )?1 + h + (1? ) 1? (" )?1 i = k (tr) + h + (1? ) (? 1) 1? (" )?1i k (tr) + [ + (1? ) (? 1)] = k (tr) + 4 (1? )? :5? 1? = + 4 N (1? ) (1? ) (tr) + 4 (1? ) 4 (1? )? :5?! 1?!? :5?! 1? 13
14 ! where the coecient 4 (1? )? :5? is a positive constant by the restrictions on 1? and. Therefore, M(") is strictly convex. Next, we show that the equation (3.6) is uniquely solvable. Obviously, we only need to prove the function f() = ( ()) is strictly increasing when >. In that case, and If 1=, then If < 1=, then f() = h + (1? ) 1??1 i f () = h + (1? ) 1??1 i h + (1? ) (? 1) 1??1 i : f () h + (1? ) 1??1 i > : f () h + (1? ) 1??1 i [ + (1? ) (? 1)] = h + (1? ) 1? i :5??1 (1? )? 1? > : Hence, f() is a strictly increasing function. Finally, we modify the inequality (3.1), which is inadequate when the parameter is close to zero. As in the last section, we have where J(u ; u )? J(u; u) = k ktr"(u? u)k L () + 1 1! h(; t; u; u ) d dt dx h(; t; u; u ) h ("(u) + t "(u? u)) ("(u) + t "(u? u)) "(u? u) Denote = "(u) + t "(u? u). When <, we have When, we have h(; t; u; u ) h(; t; u; u ) = t "(u? u) : = ( ) t ( "(u? u)) + ( ) t "(u? u) = (1? ) (? 1) 1? ( )? t ( "(u? u)) + t h + (1? ) 1? ( )?1i "(u? u) t (1? ) (? 1) 1? ( )?1 "(u? u) + t h + (1? ) 1? ( )?1i "(u? u) t h + (1? ) (? 1) 1? ( )?1 i "(u? u) : 14
15 If 1=, then h(; t; u; u ) t "(u? u) : If < 1=, then h(; t; u; u ) t [ + (1? ) (? 1)] "(u? u) = t (1? )? :5? 1?! "(u? u) : Now, we dene l (u; u ) = fx j "(u(x)) ; "(u (x)) g the intersection of the linear ealsticity region of u and that of u, and c = 8 >< >: ; if 1 ; (1? )? :5? 1?! ; if < 1 : Since on l (u; u ), = ("(u) + t "(u? u)), we have the following relation J(u ; u )? J(u; u ) k ktr"(u? u)k L () + k"(u? u)k L ( l (u;u )) + c k"(u? u)k L (n l (u;u )): (4.3) In conclusion, we have Theorem 4.1 For a nonlinear Hencky material dened by the governing function (4.1), the following a-posteriori error estimates for the linearization of the material hold. For the error in the total potential energy E(u )? E(u) h ("(u ) )? ( (? "(u )) ) i dx and for the error in energy norm ("(u ) + (? "(u ))) "(u ) dx (4.4) k ktr"(u? u)k L () + k"(u? u)k L ( l (u;u )) + c k"(u? u)k L (n l (u;u )) h ("(u ) )? ( (? "(u )) ) i dx ("(u ) + (? "(u ))) "(u ) dx (4.5) 15
16 As for numerical experiments, we consider elasticity problems on a plane corner domain (N = ). Near the reentrant corner O, the domain coincides with a cone: C! = f(r; ) j r > ;?!= < <!=g: Assume the plane elastic body is loaded by boundary tractions only and, in the neighborhood of the corner, has stress-free boundaries. We have the following singularity expansion for the solution u of the linear elasticity problem ([11]): u = K I p r I I() + K II p r II II() + ~u (4.6) G G in which, G is the modulus of rigidity, K I and K II are Mode 1 and Mode stress-intensity factors, I() = II() (? Q I( I + 1)) cos I? I cos( I ) ( + Q I ( I + 1)) sin I + I sin( I (? Q II( II + 1)) sin II? II sin( II )?( + Q II ( II + 1)) cos II? II cos( II ) correspond to the most singular Mode 1 and Mode displacement components, where Q I =? cos(( I? 1)!=) cos(( I + 1)!=) depends only on Poisson's ratio, if! < ; Q I =? 1? I sin(( I? 1)!=) 1 + I sin(( I + 1)!=) Q II =? sin(( II? 1)!=) sin(( II + 1)!=) ; = 3? 4 for plane strain; = 3? for plane stress; 1 + if! = ; and ~u is a smooth remainder such that "(~u ) is bounded when x! O. Following is a table for values of singular exponents I and II for some typical angles (the angles are in degrees). 1 A 1 A! I II
17 For the convenience of later use, let us denote: It is easy to verify following expressions: u I = r I I() u II = r II II() " 11 (u I ) = I r I?1 f(? Q I ( I + 1)? 1) cos( I? 1)? ( I? 1) cos( I? 3)g " (u I ) = I r I?1 f( + Q I ( I + 1)? 1) cos( I? 1) + ( I? 1) cos( I? 3)g " 1 (u I ) = I r I?1 fq I ( I + 1) sin( I? 1) + ( I? 1) sin( I? 3)g " 11 (u II ) = II r II?1 f(? Q II ( II + 1)? 1) sin( II? 1)? ( II? 1) sin( II? 3)g " (u II ) = II r II?1 f( + Q II ( II + 1)? 1) sin( II? 1) + ( II? 1) sin( II? 3)g " 1 (u II ) =? II r II?1 fq II ( II + 1) cos( II? 1)? ( II? 1) cos( II? 3)g Note that when the angle is 36 or 7 degree, both "(u I ) and "(u II ) are unbounded around the corner, whereas when the angle is 4, 5 or 1 degree, only the most singular Mode 1 strain components are unbounded around the corner. Experiment 4. We consider a nonlinear Hencky elasticity problem as shown in Figure 1, where homogeneous natural boundary conditions are not written explicitly. The part of the boundary represented by thick line segments is xed. Assume the stored energy function is given by (.) and (4.1). Owing to the symmetry of the problem, we consider only a quarter problem, shown in Figure, where again homogeneous natural boundary conditions are not written explicitly. We are going to use estimates (4.4) and (4.5) to examine how much error we are making if we take the material as a linear one. Let a = 5, b = 4, c = 4, h = :5. By using MSC/PROBE, a software package for solving elliptic problems by h-p-version nite element method, we obtain that for the L-shape singularity ( = j n j), K I = 1:51 (); K II = 59:9 (): The strain energy is: where, E is the Young's modulus. : () E 17
18 6 n?????????? 6? h c a - a - b -? Figure 1. A nonlinear elasticity problem u; v =????? u = u = Figure. The working problem u; v = Around the L-shape corner, "(u ) is mainly determined by the rst two terms in the expansion (4.6), since strain components for ~u are negligible. Hence, we use the rst two terms to approxiamtely represent the linear solution u for producing estimates for error in the (nonlinear) energy between the linear solution u and the nonlinear solution u, the corresponding estimates from (4.4) and/or (4.5) being denoted by Ested. For typical materials, the Young's modulus is of order psi, and the strain yield value is of the order 1?4. In the following, we take E = 1 6 psi; = :3; = 1?4 then the linear energy of the linear solution is :541 1? (). estimates for various and. We compute the 18
19 Linear Energy Ested 5 :541(+).9.8 :388(?5).5 :961(?5). :171(?4).5.8 :86(?4).5 :433(?3). :97(?3).3.8 :194(?3).5 :1414(?) :1(?).51 :8618(+).51 :864(+).51 :864(+) Experiment 4.3 We consider the Edge-cracked Panel Problem from the paper [11]. By using a mesh ner than that presented in the paper, we obtain that the linear energy of the linear solution is: 1:754 E and the Mode 1 and Mode stress intensity factors are: K I = : ; K II =?:9515 : Again, we take E = 1 6 psi; = :3; = 1?4 : Linear Energy Ested :491(?7).5 :1795(?6). :33(?6).5.8 :1696(?5).5 :877(?5). :155(?4).3.8 :4169(?5).5 :638(?4) :1(?3).51 :1614(+).51 :1615(+) 19
20 References [1] Babuska, I., The problem of modelling the elastomechanics in engineering, Computer Methods in Applied Mechanics and Engineering 8 (199), 155{18. [] Bathe, K. J., Lee, N. S. & Bucalem, M., On the use of hierarchical models in engineering analysis, Computer Methods in Applied Mechanics and Engineering 8 (199), 5{6. [3] Ekeland, I. & Temam, R., Convex Analysis and Variational Problems. North-Holland Publ. Company, Amsterdam, [4] Han, W., A regularization procedure for a simplied friction problem, Mathematical and Computer Modelling 15 (1991), No. 8, 65{7. [5] Han, W., Quantitative error estimates for material idealization of torsion problems, Mathematical and Computer Modelling 15 (1991), No. 9, 47{54. [6] Han, W., Quantitative error estimates in modeling the laminar stationary ow of a Bingham uid, Applied Mathematics and Computation 47 (199), 15{4. [7] Han, W., Finite element analysis of a holonomic elastic-plastic problem, Numer. Math. 6 (199), 493{58. [8] Han, W., A posteriori error analysis for linearizations of nonlinear elliptic problems, Mathematical Methods in the Applied Sciences 17 (1994), 487{58. [9] Han, W., Quantitative error estimates for coecient idealization in linear elliptic problems, Mathematical Methods in the Applied Sciences 17 (1994), 971{987. [1] Necas, J. & Hlavacek, I., Mathematical Theory of Elastic and Elastico-Plastic Bodies: An Introduction, Elsevier Scientic Publishing Company, Amsterdam, Oxford, New York, [11] Szabo, B. A. & Babuska, I., Computation of the amplitude of stress singular terms for cracks and reentrant corners, Fracture Mechanics: Nineteenth Symposium, ASTM STP969, 1988, Ed. T.A. Cruse. [1] eidler, E., Nonlinear Functional Analysis and its Applications: IV. Applications to Mathematical Physics. Springer-Verlag, New York, 1988.
Weak 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 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 informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationPARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN. H.T. Banks and Yun Wang. Center for Research in Scientic Computation
PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN H.T. Banks and Yun Wang Center for Research in Scientic Computation North Carolina State University Raleigh, NC 7695-805 Revised: March 1993 Abstract In
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. 155 A posteriori error estimates for stationary slow flows of power-law fluids Michael Bildhauer,
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 informationA mixed nite volume element method based on rectangular mesh for biharmonic equations
Journal of Computational and Applied Mathematics 7 () 7 3 www.elsevier.com/locate/cam A mixed nite volume element method based on rectangular mesh for biharmonic equations Tongke Wang College of Mathematical
More information2 Ben Schweizer As Hocking describes in [5], a discussion focused on the following question: Can one assume that the dynamic contact angle is a consta
A well-posed model for dynamic contact angles Ben Schweizer February 27, 1998 Abstract We consider uid systems with a free boundary and with a point of contact of the free boundary with a solid wall. We
More informationBifurcation from the rst eigenvalue of some nonlinear elliptic operators in Banach spaces
Nonlinear Analysis 42 (2000) 561 572 www.elsevier.nl/locate/na Bifurcation from the rst eigenvalue of some nonlinear elliptic operators in Banach spaces Pavel Drabek a;, Nikos M. Stavrakakis b a Department
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 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 informationA 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 informationNumerical Integration for Multivariable. October Abstract. We consider the numerical integration of functions with point singularities over
Numerical Integration for Multivariable Functions with Point Singularities Yaun Yang and Kendall E. Atkinson y October 199 Abstract We consider the numerical integration of functions with point singularities
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 informationAbstract. A front tracking method is used to construct weak solutions to
A Front Tracking Method for Conservation Laws with Boundary Conditions K. Hvistendahl Karlsen, K.{A. Lie, and N. H. Risebro Abstract. A front tracking method is used to construct weak solutions to scalar
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 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. Mathematik Preprint Nr. 99 Duality based a posteriori error estimates for higher order variational inequalities with
More informationON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND STEFAN DIEHL
ON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND DISCONTINUOUS FLUX FUNCTION STEFAN DIEHL Abstract. The conservation law studied is @u(x;t) + @ F u(x; t); x @t @x = s(t)(x), where u is a f(u); x > 0 concentration,
More informationLimit Analysis with the. Department of Mathematics and Computer Science. Odense University. Campusvej 55, DK{5230 Odense M, Denmark.
Limit Analysis with the Dual Ane Scaling Algorithm Knud D. Andersen Edmund Christiansen Department of Mathematics and Computer Science Odense University Campusvej 55, DK{5230 Odense M, Denmark e-mail:
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationAn explicit nite element method for convection-dominated compressible viscous Stokes system with inow boundary
Journal of Computational and Applied Mathematics 156 (2003) 319 343 www.elsevier.com/locate/cam An explicit nite element method for convection-dominated compressible viscous Stokes system with inow boundary
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 informationAn Optimal Finite Element Mesh for Elastostatic Structural. Analysis Problems. P.K. Jimack. University of Leeds. Leeds LS2 9JT, UK.
An Optimal Finite Element Mesh for Elastostatic Structural Analysis Problems P.K. Jimack School of Computer Studies University of Leeds Leeds LS2 9JT, UK Abstract This paper investigates the adaptive solution
More informationPointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang
Pointwise convergence rate for nonlinear conservation laws Eitan Tadmor and Tao Tang Abstract. We introduce a new method to obtain pointwise error estimates for vanishing viscosity and nite dierence approximations
More informationVariational principles in mechanics
CHAPTER Variational principles in mechanics.1 Linear Elasticity n D Figure.1: A domain and its boundary = D [. Consider a domain Ω R 3 with its boundary = D [ of normal n (see Figure.1). The problem of
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 informationDuality method in limit analysis problem of non-linear elasticity
Duality method in limit analysis problem of non-linear elasticity Igor A. Brigadnov Department of Computer Science North-Western State Technical University Millionnaya 5, St. Petersburg, 191186, Russia
More informationTechinical Proofs for Nonlinear Learning using Local Coordinate Coding
Techinical Proofs for Nonlinear Learning using Local Coordinate Coding 1 Notations and Main Results Denition 1.1 (Lipschitz Smoothness) A function f(x) on R d is (α, β, p)-lipschitz smooth with respect
More informationDUALITY BASED A POSTERIORI ERROR ESTIMATES FOR HIGHER ORDER VARIATIONAL INEQUALITIES WITH POWER GROWTH FUNCTIONALS
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 33, 2008, 475 490 DUALITY BASED A POSTERIORI ERROR ESTIMATES FOR HIGHER ORDER VARIATIONAL INEQUALITIES WITH POWER GROWTH FUNCTIONALS Michael Bildhauer,
More informationLECTURE 15 + C+F. = A 11 x 1x1 +2A 12 x 1x2 + A 22 x 2x2 + B 1 x 1 + B 2 x 2. xi y 2 = ~y 2 (x 1 ;x 2 ) x 2 = ~x 2 (y 1 ;y 2 1
LECTURE 5 Characteristics and the Classication of Second Order Linear PDEs Let us now consider the case of a general second order linear PDE in two variables; (5.) where (5.) 0 P i;j A ij xix j + P i,
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationon! 0, 1 and 2 In the Zienkiewicz-Zhu SPR p 1 and p 2 are obtained by solving the locally discrete least-squares p
Analysis of a Class of Superconvergence Patch Recovery Techniques for Linear and Bilinear Finite Elements Bo Li Zhimin Zhang y Abstract Mathematical proofs are presented for the derivative superconvergence
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 informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More information436 A. Barani and G.H. Rahimi assessment models have been employed to investigate the LBB of cracked pipes that are not for combined load [8]. Yun-Jae
Scientia Iranica, Vol. 4, No. 5, pp 435{44 c Sharif University of Technology, October 27 Approximate Method for Evaluation of the J-Integral for Circumferentially Semi-Elliptical-Cracked Pipes Subjected
More informationNonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface
Transaction A: Civil Engineering Vol. 16, No. 6, pp. 512{519 c Sharif University of Technology, December 2009 Research Note Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield
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 informationOptimal thickness of a cylindrical shell under dynamical loading
Optimal thickness of a cylindrical shell under dynamical loading Paul Ziemann Institute of Mathematics and Computer Science, E.-M.-A. University Greifswald, Germany e-mail paul.ziemann@uni-greifswald.de
More informationA new nite-element formulation for electromechanical boundary value problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2002; 55:613 628 (DOI: 10.1002/nme.518) A new nite-element formulation for electromechanical boundary value problems
More informationON TRIVIAL GRADIENT YOUNG MEASURES BAISHENG YAN Abstract. We give a condition on a closed set K of real nm matrices which ensures that any W 1 p -grad
ON TRIVIAL GRAIENT YOUNG MEASURES BAISHENG YAN Abstract. We give a condition on a closed set K of real nm matrices which ensures that any W 1 p -gradient Young measure supported on K must be trivial the
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationContents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping.
Minimization Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping. 1 Minimization A Topological Result. Let S be a topological
More informationFilters in topology optimization
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2001; 50:2143 2158 Filters in topology optimization Blaise Bourdin ; Department of Mathematics; Technical University
More information/00 $ $.25 per page
Contemporary Mathematics Volume 00, 0000 Domain Decomposition For Linear And Nonlinear Elliptic Problems Via Function Or Space Decomposition UE-CHENG TAI Abstract. In this article, we use a function decomposition
More informationt x 0.25
Journal of ELECTRICAL ENGINEERING, VOL. 52, NO. /s, 2, 48{52 COMPARISON OF BROYDEN AND NEWTON METHODS FOR SOLVING NONLINEAR PARABOLIC EQUATIONS Ivan Cimrak If the time discretization of a nonlinear parabolic
More informationARE202A, Fall Contents
ARE202A, Fall 2005 LECTURE #2: WED, NOV 6, 2005 PRINT DATE: NOVEMBER 2, 2005 (NPP2) Contents 5. Nonlinear Programming Problems and the Kuhn Tucker conditions (cont) 5.2. Necessary and sucient conditions
More informationDYNAMIC WEIGHT FUNCTIONS FOR A MOVING CRACK II. SHEAR LOADING. University of Bath. Bath BA2 7AY, U.K. University of Cambridge
DYNAMIC WEIGHT FUNCTIONS FOR A MOVING CRACK II. SHEAR LOADING A.B. Movchan and J.R. Willis 2 School of Mathematical Sciences University of Bath Bath BA2 7AY, U.K. 2 University of Cambridge Department of
More informationSecond order sucient optimality condition for a. nonlinear elliptic control problem. Eduardo Casas and Fredi Troltzsch and Andreas Unger
Second order sucient optimality condition for a nonlinear elliptic control problem by Eduardo Casas and Fredi Troltzsch and Andreas Unger Technical University of Chemnitz wickau Faculty of Mathematics
More informationMARY ANN HORN be decoupled into three wave equations. Thus, we would hope that results, analogous to those available for the wave equation, would hold
Journal of Mathematical Systems, Estimation, and Control Vol. 8, No. 2, 1998, pp. 1{11 c 1998 Birkhauser-Boston Sharp Trace Regularity for the Solutions of the Equations of Dynamic Elasticity Mary Ann
More informationA metric space X is a non-empty set endowed with a metric ρ (x, y):
Chapter 1 Preliminaries References: Troianiello, G.M., 1987, Elliptic differential equations and obstacle problems, Plenum Press, New York. Friedman, A., 1982, Variational principles and free-boundary
More informationCRACK-TIP DRIVING FORCE The model evaluates the eect of inhomogeneities by nding the dierence between the J-integral on two contours - one close to th
ICF 100244OR Inhomogeneity eects on crack growth N. K. Simha 1,F.D.Fischer 2 &O.Kolednik 3 1 Department ofmechanical Engineering, University of Miami, P.O. Box 248294, Coral Gables, FL 33124-0624, USA
More informationTable 1: BEM as a solution method for a BVP dierential formulation FDM BVP integral formulation FEM boundary integral formulation BEM local view is ad
Chapter 1 Introduction to Boundary element Method - 1D Example For reference: Hong-Ki Hong and Jeng-Tzong Chen, Boundary Element Method, Chapter 1 Introduction to Boundary element Method - 1D Example,
More informationWeak solutions for some quasilinear elliptic equations by the sub-supersolution method
Nonlinear Analysis 4 (000) 995 100 www.elsevier.nl/locate/na Weak solutions for some quasilinear elliptic equations by the sub-supersolution method Manuel Delgado, Antonio Suarez Departamento de Ecuaciones
More informationCOMBINED EFFECTS FOR A STATIONARY PROBLEM WITH INDEFINITE NONLINEARITIES AND LACK OF COMPACTNESS
Dynamic Systems and Applications 22 (203) 37-384 COMBINED EFFECTS FOR A STATIONARY PROBLEM WITH INDEFINITE NONLINEARITIES AND LACK OF COMPACTNESS VICENŢIU D. RĂDULESCU Simion Stoilow Mathematics Institute
More informationexpression that describes these corrections with the accuracy of the order of 4. frame usually connected with extragalactic objects.
RUSSIAN JOURNAL OF EARTH SCIENCES, English Translation, VOL, NO, DECEMBER 998 Russian Edition: JULY 998 On the eects of the inertia ellipsoid triaxiality in the theory of nutation S. M. Molodensky Joint
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationSemi-strongly asymptotically non-expansive mappings and their applications on xed point theory
Hacettepe Journal of Mathematics and Statistics Volume 46 (4) (2017), 613 620 Semi-strongly asymptotically non-expansive mappings and their applications on xed point theory Chris Lennard and Veysel Nezir
More informationLinear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space
Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................
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 informationThe L p -dissipativity of first order partial differential operators
The L p -dissipativity of first order partial differential operators A. Cialdea V. Maz ya n Memory of Vladimir. Smirnov Abstract. We find necessary and sufficient conditions for the L p -dissipativity
More informationSELF-ADJOINTNESS OF SCHRÖDINGER-TYPE OPERATORS WITH SINGULAR POTENTIALS ON MANIFOLDS OF BOUNDED GEOMETRY
Electronic Journal of Differential Equations, Vol. 2003(2003), No.??, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) SELF-ADJOINTNESS
More information2 GOVERNING EQUATIONS
2 GOVERNING EQUATIONS 9 2 GOVERNING EQUATIONS For completeness we will take a brief moment to review the governing equations for a turbulent uid. We will present them both in physical space coordinates
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
More informationOptimization and Optimal Control in Banach Spaces
Optimization and Optimal Control in Banach Spaces Bernhard Schmitzer October 19, 2017 1 Convex non-smooth optimization with proximal operators Remark 1.1 (Motivation). Convex optimization: easier to solve,
More informationThe Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria
ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-19 Wien, Austria The Negative Discrete Spectrum of a Class of Two{Dimentional Schrodinger Operators with Magnetic
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 informationConverse Lyapunov Functions for Inclusions 2 Basic denitions Given a set A, A stands for the closure of A, A stands for the interior set of A, coa sta
A smooth Lyapunov function from a class-kl estimate involving two positive semidenite functions Andrew R. Teel y ECE Dept. University of California Santa Barbara, CA 93106 teel@ece.ucsb.edu Laurent Praly
More informationConservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body
Continuum Mechanics We ll stick with the Lagrangian viewpoint for now Let s look at a deformable object World space: points x in the object as we see it Object space (or rest pose): points p in some reference
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationOn the p-laplacian and p-fluids
LMU Munich, Germany Lars Diening On the p-laplacian and p-fluids Lars Diening On the p-laplacian and p-fluids 1/50 p-laplacian Part I p-laplace and basic properties Lars Diening On the p-laplacian and
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 informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
More informationLecture 2: Review of Prerequisites. Table of contents
Math 348 Fall 217 Lecture 2: Review of Prerequisites Disclaimer. As we have a textbook, this lecture note is for guidance and supplement only. It should not be relied on when preparing for exams. In this
More information7 To solve numerically the equation of motion, we use the velocity Verlet or leap frog algorithm. _ V i n = F i n m i (F.5) For time step, we approxim
69 Appendix F Molecular Dynamics F. Introduction In this chapter, we deal with the theories and techniques used in molecular dynamics simulation. The fundamental dynamics equations of any system is the
More informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
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. 166 A short remark on energy functionals related to nonlinear Hencky materials Michael Bildhauer
More informationwith Applications to Elasticity and Compressible Flow Daoqi Yang March 20, 1997 Abstract
Stabilized Schemes for Mixed Finite Element Methods with Applications to Elasticity and Compressible Flow Problems Daoqi Yang March 20, 1997 Abstract Stabilized iterative schemes for mixed nite element
More informationEngineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February.
Engineering Sciences 241 Advanced Elasticity, Spring 2001 J. R. Rice Homework Problems / Class Notes Mechanics of finite deformation (list of references at end) Distributed Thursday 8 February. Problems
More informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
More informationON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS. Abstract. We introduce a wide class of asymmetric loss functions and show how to obtain
ON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS FUNCTIONS Michael Baron Received: Abstract We introduce a wide class of asymmetric loss functions and show how to obtain asymmetric-type optimal decision
More informationChapter 5. Sound Waves and Vortices. 5.1 Sound waves
Chapter 5 Sound Waves and Vortices In this chapter we explore a set of characteristic solutions to the uid equations with the goal of familiarizing the reader with typical behaviors in uid dynamics. Sound
More informationLoading σ Stress. Strain
hapter 2 Material Non-linearity In this chapter an overview of material non-linearity with regard to solid mechanics is presented. Initially, a general description of the constitutive relationships associated
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 informationDuality, Dual Variational Principles
Duality, Dual Variational Principles April 5, 2013 Contents 1 Duality 1 1.1 Legendre and Young-Fenchel transforms.............. 1 1.2 Second conjugate and convexification................ 4 1.3 Hamiltonian
More informationWALTER LITTMAN BO LIU kind of problem is related to the deformation of a membrane with strings attached to several sides of the polygon. In other word
THE REGULARITY AND SINGULARITY OF SOLUTIONS OF CERTAIN ELLIPTIC PROBLEMS ON POLYGONAL DOMAINS Walter Littman Bo Liu School of Math., University of Minnesota, Minneapolis, MN55455 Abstract. The regularity
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
More informationA Geometric Interpretation of Newtonian Impacts with Global Dissipation Index
A Geometric Interpretation of Newtonian Impacts with Global Dissipation Index Christoph Glocker Institute B of Mechanics, Technical University of Munich, D-85747 Garching, Germany (glocker@lbm.mw.tu-muenchen.de)
More information3.1 Basic properties of real numbers - continuation Inmum and supremum of a set of real numbers
Chapter 3 Real numbers The notion of real number was introduced in section 1.3 where the axiomatic denition of the set of all real numbers was done and some basic properties of the set of all real numbers
More informationBOUNDARY CONTROLLABILITY IN PROBLEMS JOHN E. LAGNESE. linear hyperbolic systems having piecewise constant coecients in
ESAIM: Control, Optimisation and Calculus of Variations URL: http://www.emath.fr/cocv/,, BOUNDARY CONTROLLABILITY IN PROBLEMS OF TRANSMISSION FOR A CLASS OF SECOND ORDER HYPERBOLIC SYSTEMS JOHN E. LAGNESE
More informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationAlternative theorems for nonlinear projection equations and applications to generalized complementarity problems
Nonlinear Analysis 46 (001) 853 868 www.elsevier.com/locate/na Alternative theorems for nonlinear projection equations and applications to generalized complementarity problems Yunbin Zhao a;, Defeng Sun
More informationMODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione
MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione, Univ. di Roma Tor Vergata, via di Tor Vergata 11,
More informationRearrangements and polar factorisation of countably degenerate functions G.R. Burton, School of Mathematical Sciences, University of Bath, Claverton D
Rearrangements and polar factorisation of countably degenerate functions G.R. Burton, School of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, U.K. R.J. Douglas, Isaac Newton
More informationIntegral equations for crack systems in a slightly heterogeneous elastic medium
Boundary Elements and Other Mesh Reduction Methods XXXII 65 Integral equations for crack systems in a slightly heterogeneous elastic medium A. N. Galybin & S. M. Aizikovich Wessex Institute of Technology,
More informationSUB-SUPERSOLUTIONS IN A VARIATIONAL INEQUALITY RELATED TO A SANDPILE PROBLEM
ANZIA J. 48(2006, 179 197 SUB-SUPERSOLUTIONS IN A VARIATIONAL INEQUALITY RELATED TO A SANDPILE PROBLE VY KHOI LE 1 (Received 29 April, 2005; revised 12 October, 2006 Abstract In this paper we study a variational
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 information