Conjugate gradient-boundary element method for a Cauchy problem in the Lam6 system
|
|
- Chester Cunningham
- 5 years ago
- Views:
Transcription
1 Conjugate gradient-boundary element method for a Cauchy problem in the Lam6 system L. Marin, D.N. Hho & D. Lesnic I Department of Applied Mathematics, University of Leeds, UK 'Hanoi Institute of Mathematics, Vietnam Abstract In this paper an iterative algorithm, based on the conjugate gradient method (CGM), for obtaining approximate solutions to the Cauchy problem in linear elasticity is analysed, using the boundary element method (BEM). The numerical results obtained confirm that the iterative BEM produces a convergent and stable numerical solution with respect to increasing the number of boundary elements and decreasing the amount of noise added into the input data. 1 Introduction Consider a linear elastic material which occupies an open bounded domain R C Rd, where d is the dimension o' the space in which the problem is posed, usually d E {l, 2,3), and assume that Cl is bounded by a surface I' = dr E C'. We also assume that the boundary consists of two parts, I' = I'l U r2, where rl, r2 # 0 and I'l n r2 = 0. In the absence of body forces, the equilibrium equations are given by, see e.g. Saada [l],
2 where ui, is the stress tensor and the strain tensor E;, is given by These tensors are related by the constitutive law, namely where Cijkl is the elasticity tensor and, in the case of an isotropic linear elastic material, is given by with G and v the shear modulus and Poisson ratio, respectively, and 6ij the Kronecker delta tensor. If we now substitute the constitutive law (3) into the equilibrium equations (l), and use the kinematic relations (2) together with the expression (4) for the elasticity tensor for an isotropic linear elastic material, we obtain the Lam6 system, namely, We now let n(x) be the outward normal vector at I' and t(x) be the traction vector at a point X E l? whose components are defined by Assuming that both the displacement and traction vectors can be measured on a part of the boundary l?, say r2, then this leads to the mathematical formulation of an inverse problem consisting of equations (1) or (5) and the boundary conditions where G and are prescribed vector valued functions. This problem, termed a Cauchy problem, is much more difficult to solve both analytically and numerically than the direct problem, since the solution does not satisfy the general conditions of well-posedness. Although the problem may have a unique solution, it is well known, see e.g. Hadamard [2], that this solution is unstable with respect to small perturbations in the data on r2. Thus the problem is ill-posed and we cannot use a direct approach, such as the Gauss elimination method, in order to solve the system of linear equations which arises from the discretisation of the partial differential equations (1) or (5) and the boundary conditions (7). Therefore, knowing the exact data on the
3 boundary r2, we apply a variational method to the aforementioned Cauchy problem. Since the boundary conditions on rl are unknown and have to be determined, we consider the displacement vector on the underspecified boundary as a control for a direct problem and attempt to fit the Cauchy data on the overspecified boundary r2 by minimizing a functional relating the known and calculated values of the displacement vector on r2, see e.g. Lions [3]. 2 Variational formulation Let us denote by ll2(ri), WS (ri) and W' (R) the spaces (~~(ri)) d, (~~(ri)) and (HS(R)) d, respectively, for d E {l, 2,3), i = 1,2 and some s E R, see e.g. Lions and Magenes [4]. Then the Cauchy ~roblem for the Lamk system may be recast as follows - where G E IL2(r2), 2 E ll2(r2) and U is sought in RI1l2(R). We note that t E JHI-l(r2) is sufficient for our method. Let us denote by -yi f, i = 1,2, the trace of a function f determined in R over ri, i = 1,2. First we solve the direct problem with v E L2(rl). We denote by U = u(v,'t) the solution to the problem (9) and aim to find v E IL2(rl) such that To do so, we attempt to minimize the functional with respect to v E ll2(fl). Consider the problem where p E ll2(r2).
4 232 Bourldai-y Elemenr Technology XIV Lemma 1 Let U and 1C, be the solutions of the problems (9) and (12), respectively. Then Theorem 1 Let -i E IL2(r2). If v varies in IL2(rl), then y2u(v,z) forms a dense set in IL2(rl). This theorem has important consequences, for example, we can say that the Cauchy problem (8) is solvable for almost all G,z E IL2(r2). Furthermore, we have the following result: Corollary 1 inf J(v) = 0 vcla(ra) (14) Theorem 2 The functional J(v) is Fre'chet differentiable and its first gradient has the form Furthermore, the functional J(v) is twice Fre'chet differentiable and it is strictly convex. 3 Conjugate gradient method (CGM) As we can calculate the gradient of the functional J(v) via the adjoint problem, we can now apply the conjugate gradient method (CGM) with a stopping rule, as proposed by Nemirovskii [5]. We define the linear operator and we have the following linear equation Suppose that instead of G we have only an approximation of it, say G, E IL2(r2) such that llg - G IILyr,) < - E (18) To solve our equation (8) with noisy data G,, we need to compute A+o(Aov - - U,), where A: is the adjoint of the operator A0 and C, is given by
5 Bounda~y Element Technology XIV 233 However, we observe that this is nothing else than the gradient of the functional 1 J(v) = -IIAv 2 - Gc Il~yr,) (20) Thus the CGM applied to our problem has the form of the following algorithm: Step 1. Set k = 0. Choose u(o) E L2(r2). Step 2. Solve the direct problem in order to determine the residual T ( ~ ) - - T ( ~ = ) AU(~) - G, = y2u(u(k), t) - U, Step 3. Solve the adjoint problem to determine the gradient g(k) (k) g; (X) = 71 (flij($'(oi r(k))(~))nj(~)) Calculate Pk and d(k) as follows: k=o: pk=o, d(k) = 2 l : pk= llg(k)ll~ypl) k-l) llg' lllyr1) ' Step 4. Solve the direct problem d(k) = -g(k) + ~ ~d(k-1) ( 25) in order to determine ~ od(~) given by
6 Compute a k and u(~+') as follows: Step 5. Set k = k + 1. Repeat Step 2 until a stopping criterion is prescribed. As a stopping criterion we choose the one suggested by Nemirovskii [5], namely E N : ~ ( ~ ) l 5 ~ 6~ ~ ( ~ ~ ) (29) where 6 is a constant which can be taken heuristically as suggested by Hanke and Hansen [6]. It follows from Nemirovskii's result that the above iterative procedure converges with an optimal convergence rate to the normal solution of the problem as the noise level tends to zero. A boundary element method (BEM), see e.g. Marin et al. [7], is used in order to solve the intermediate mixed well-posed boundary value problems resulting from the iterative method adopted. 4 Numerical results and discussion In order to present the performance of the numerical method proposed, we solve the Cauchy problem for examples in a smooth geometry, namely the unit disc S2 = {X = (xl, x2) I X: + X; < 1). We assume that the boundary J? = {X = (xl, 22) I X: + X; = 1) of the solution domain R is divided into two disjointed parts, namely rl = {X = (xl, 22) I X E l?, cul < cr2) and r2 = {X = (xl,x2) I X E l?, 0 < a1) U {X = (x1,x2) I X E r, a2 5 2x1, is the angular polar coordinate of X and a*, i = 1,2, are specified angles in the interval (0,27r). In order to illustrate the typical numerical results we have taken a1 = 7~/4 and a2 = 3x14 and we assume that the boundary r2 is overspecified by the prescription of both the displacement and the traction vectors while the boundary rl is underspecified with both the displacement and the traction vectors unknown. In the following test examples we consider an isotropic linear elastic medium characterized by the material constants G = 3.35 X 101 N/m2 and v = 0.34 corresponding to a copper alloy. Example 1. We consider the following analytical solution in displacements l-v (an) U' ("I = 2G(l+ v) ffo Xi i = 1,2 I in the domain fl, which corresponds to a uniform hydrostatic stress state given by
7 with a 0 = 1.5 X 1010 N/m2. Example 2. We consider the following tractions on the boundary I' where a0 = 1.5 X 10'' N/m2 and B(x) is the angular polar coordinate of X. It should be noted that the corresponding analytical expressions for the stresses uj~)(x) and displacements u!'")(r) are not available in this case for the tractions (32) and (33), but they can be obtained numerically by solving the direct problem An arbitrary vector valued function U(') E L2(rl) X L2(rl) may be specified as an initial guess for the displacement vector on PI, but in order to improve the rate of convergence of the iterative procedure we have chosen a vector valued function which ensures the continuity of the displacement vector at the endpoints of rl and which is also linear with respect to the angular polar coordinate 8. For the test examples considered, this initial guess is given by ui (0)(X) = a2 - e(x)up(xll a1 (an) a2 - a1 a 2 - a1 ui (22) (35) for i = 1,2, where ai = B(xi), xi are the endpoints of rl, and the choice of a1 = a/4 and a 2 = 37r/4 also ensures that the initial guess is not too close to the exact values uian)(x). In order to investigate the convergence of the algorithm, at every iteration we evaluate the accuracy errors where u(~) is the displacement vector on the boundary rl retrieved after k iterations, the operator A is given by equation (10) and each iteration consists of solving three direct mixed well-posed problems as described in section 3. When starting with the initial guess U(') given by equation (35) for exact input data, a convergent sequence { u(~))~>~ - of approximation functions for
8 236 Bourida~?i Elerneut Tectviology XIV Figure 1: The analytical solution U(;") (-) and the numerical solution uprn) retrieved for N = 40 (o) and N = 80 (v) constant boundary ele- ments, variable relaxation factor with amplitude A = 2.0, and input data u(an)lr2, on the underspecified boundary rl, for the Cauchy problem considered in example 2. ulr, is obtained. However, if we evaluate the error e, at every iteration using different numbers of constant boundary elements for discretising r' when input data contains noise, we note that the error e, decreases up to a specific iteration after which it starts increasing. The error E, keeps decreasing as the number of iterations Ic increases, but there is a threshold of optimality at which the iterative process should be stopped according to Nemirovskii's rule. We note that the CGM algorithm described in section 3 is convergent as we increase the number of boundary elements. However, the errors in predicting the traction vector t on the underspecified boundary are still large since we are using the analytical input data on r2 which is contaminated by numerical noise. An alternative way to generate the Cauchy data on r2 is to use the numerical solution of the direct problem and this procedure can also be used to fabricate the input data when no exact solution to the Cauchy problem is available. In this case both accuracy errors e, and E, decrease as the number of iterations k increases and the sequence {u(~))~>~ - of approximation functions for ulr, converges
9 exactly to the analytical solution u(an)lr,. However, the numerical solution for the traction vector deviates from the exact solution, especially near the ends of the underspecified boundary which is the region of high illposedeness and also where the BEM changes to mixed boundary conditions. In order to improve the results for the traction vector on the underspecified part of the boundary, we relax the marching condition (28) through the use of U(k+l) = pu(k+l) + (1 - p)u(k) (39) when passing from step 4 to step 5 of the algorithm described in section 3, with the variable relaxation factor p = p(o(x)) given by where O(x) is the angular polar coordinate of the point X E rl, cul and a2 are the angular polar coordinates of the endpoints of the boundary rl, and A E (0,2]. It should be mentioned that when using the variable relaxation factor given by equation (39), both accuracy errors e, and E, decrease continuously as the number of iterations k increases. In Figures 1 we present the numerical solutions obtained for the 22 component of the displacement vector for the Cauchy problems considered in example 2, when various numbers of boundary elements and the variable relaxation factor given by (40) with the amplitude A = 2.0 are employed. We note that the numerical solution is convergent with respect to increasing the number of boundary elements. The stability of the numerical method proposed has been investigated by perturbing the initial data ulr2 as where 6ui is a Gaussian random variable with mean zero and standard deviation U, generated by the NAG subroutine G05DDF, and p is the percentage of additive noise included in the input data ulp2 in order to simulate the inherent measurements errors. The optimal iteration numbers and the accuracy errors indicate that as p decreases then the numerical solutions approximate better the exact solutions, whilst at the same time remaining stable. 5 Conclusions In this paper we have formulated a Cauchy problem for the Lam6 system in a variational form where only weak requirements for the Cauchy data are
10 required. Consequently, the solutions of the direct problems, as well as the associated adjoint problems, are defined in a weak sense and a mathematical analysis has been undertaken. The variational approach for solving the Cauchy problem in elasticity needs the gradient of the minimization functional, which is provided by the solution of the adjoint problem. Due to the explicit representation of the gradient, the CGM was employed to solve numerically the Cauchy problem. The algorithm proposed consists of solving three direct problems at every iteration, but because of the linearity of the problem only two direct solutions are required at every iteration. In combination with Nemirovskii's stopping criterion, the CGM is known to be of optimal order when the data is sufficiently smooth. The use of a variable relaxation factor with respect to the angular polar coordinate substantially increased the accuracy of the numerical solution. Cauchy problems are inverse boundary value problems and, therefore, only the discretisation of the boundary is needed and thus the BEM is a suitable method for solving such improperly posed problems. The numerical results obtained for various numbers of constant boundary elements and various amounts of noise added into the input data were found to be in good agreement with the exact solution. References [l] A.S. Saada, Elasticity: Theory and Applications (Pergamon Press, New York 1974). [2] J. Hadamard, Lectures on Cauchy's Problem in linear partial differential equations (London, Oxford University Press 1923). [3] J.L. Lions, Optimal Control of Systems Governed by Partial Differential Equations (Springer-Verlag, Berlin 1971). [4] J.L. Lions, and E. Magenes, Non-homogeneous Boundary Value Problems and Their Applications (Springer-Verlag, New York-Heidelberg 1972). [5] A.S. Nemirovskii, The Regularizing Properties of the Adjoint Gradient Method in Ill-Posed Problems, Comput. Maths. Math. Phys. 26 (1986) [6] M. Ilanke and P.C. Hansen, Regularization Methods for Large-scale Problems, Surveys Math. in Industry 3 (1993) [7] L. Marin, L. Elliott, D.B. Ingham and D. Lesnic, Boundary Element Method for the Cauchy Problem in Linear Elasticity, submitted for publication in Engineering Analysis with Boundary Elements (2000).
Transactions on Modelling and Simulation vol 8, 1994 WIT Press, ISSN X
Boundary element method for an improperly posed problem in unsteady heat conduction D. Lesnic, L. Elliott & D.B. Ingham Department of Applied Mathematical Studies, University of Leeds, Leeds LS2 9JT, UK
More informationThe detection of subsurface inclusions using internal measurements and genetic algorithms
The detection of subsurface inclusions using internal measurements and genetic algorithms N. S. Meral, L. Elliott2 & D, B, Ingham2 Centre for Computational Fluid Dynamics, Energy and Resources Research
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More information3rd Int. Conference on Inverse Problems in Engineering AN INVERSE PROBLEM FOR SLOW VISCOUS INCOMPRESSIBLE FLOWS. University of Leeds
Inverse Problems in Engineering: Theory and Practice 3rd Int. Conference on Inverse Problems in Engineering June 13-18, 1999, Port Ludlow, WA, USA ME06 AN INVERSE PROBLEM FOR SLOW VISCOUS INCOMPRESSIBLE
More informationDetermination of a space-dependent force function in the one-dimensional wave equation
S.O. Hussein and D. Lesnic / Electronic Journal of Boundary Elements, Vol. 12, No. 1, pp. 1-26 (214) Determination of a space-dependent force function in the one-dimensional wave equation S.O. Hussein
More informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More information3rd Int. Conference on Inverse Problems in Engineering. Daniel Lesnic. University of Leeds. Leeds, West Yorkshire LS2 9JT
Inverse Problems in Engineering: Theory and Practice 3rd Int. Conference on Inverse Problems in Engineering June 13-18, 1999, Port Ludlow, WA, USA APP RETRIEVING THE FLEXURAL RIGIDITY OF A BEAM FROM DEFLECTION
More informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
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 THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX- ON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS N. S. HOANG AND A. G. RAMM (Communicated
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More informationConfigurational Forces as Basic Concepts of Continuum Physics
Morton E. Gurtin Configurational Forces as Basic Concepts of Continuum Physics Springer Contents 1. Introduction 1 a. Background 1 b. Variational definition of configurational forces 2 с Interfacial energy.
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 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 informationelastoplastic contact problems D. Martin and M.H. Aliabadi Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK
Non-conforming BEM elastoplastic contact problems D. Martin and M.H. Aliabadi discretisation in Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK Abstract In this paper,
More informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
More informationITERATIVE METHODS FOR SOLVING A NONLINEAR BOUNDARY INVERSE PROBLEM IN GLACIOLOGY S. AVDONIN, V. KOZLOV, D. MAXWELL, AND M. TRUFFER
ITERATIVE METHODS FOR SOLVING A NONLINEAR BOUNDARY INVERSE PROBLEM IN GLACIOLOGY S. AVDONIN, V. KOZLOV, D. MAXWELL, AND M. TRUFFER Abstract. We address a Cauchy problem for a nonlinear elliptic PDE arising
More informationLecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity
Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity by Borja Erice and Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationJ. Sladek, V. Sladek & M. Hrina Institute of Construction and Architecture, Slovak Academy of Sciences, Bratislava, Slovakia
Evaluation of fracture parameters for functionally gradient materials J. Sladek, V. Sladek & M. Hrina Institute of Construction and Architecture, Slovak Academy of Sciences, 842 20 Bratislava, Slovakia
More information(2) 2. (3) 2 Using equation (3), the material time derivative of the Green-Lagrange strain tensor can be obtained as: 1 = + + +
LAGRANGIAN FORMULAION OF CONINUA Review of Continuum Kinematics he reader is referred to Belytscho et al. () for a concise review of the continuum mechanics concepts used here. he notation followed here
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 informationLinear Constitutive Relations in Isotropic Finite Viscoelasticity
Journal of Elasticity 55: 73 77, 1999. 1999 Kluwer Academic Publishers. Printed in the Netherlands. 73 Linear Constitutive Relations in Isotropic Finite Viscoelasticity R.C. BATRA and JANG-HORNG YU Department
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationThe Finite Element Method for Computational Structural Mechanics
The Finite Element Method for Computational Structural Mechanics Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) January 29, 2010 Martin Kronbichler (TDB) FEM for CSM January
More informationWenyong Pan and Lianjie Huang. Los Alamos National Laboratory, Geophysics Group, MS D452, Los Alamos, NM 87545, USA
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 019 SGP-TR-14 Adaptive Viscoelastic-Waveform Inversion Using the Local Wavelet
More informationMAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design
MAE4700/5700 Finite Element Analsis for Mechanical and Aerospace Design Cornell Universit, Fall 2009 Nicholas Zabaras Materials Process Design and Control Laborator Sible School of Mechanical and Aerospace
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Preliminary Math Concept of Stress Stress Components Equilibrium
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Plate-soil elastodynamic coupling using analysis S.F.A. Baretto, H.B. Coda, W.S. Venturini Sao Carlos School of Engineering, University ofsao Paulo, Sao Carlos - SP, Brazil BEM Abstract The aim of this
More informationCOMPUTATIONAL ELASTICITY
COMPUTATIONAL ELASTICITY Theory of Elasticity and Finite and Boundary Element Methods Mohammed Ameen Alpha Science International Ltd. Harrow, U.K. Contents Preface Notation vii xi PART A: THEORETICAL ELASTICITY
More informationLinearized theory of elasticity
Linearized theory of elasticity Arie Verhoeven averhoev@win.tue.nl CASA Seminar, May 24, 2006 Seminar: Continuum mechanics 1 Stress and stress principles Bart Nowak March 8 2 Strain and deformation Mark
More informationGEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS
Methods in Geochemistry and Geophysics, 36 GEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS Michael S. ZHDANOV University of Utah Salt Lake City UTAH, U.S.A. 2OO2 ELSEVIER Amsterdam - Boston - London
More informationThe force on a lattice defect in an elastic body
Journal of Elasticity 17 (1987) 3-8 @ Martinus Nijhoff Publishers, Dordrecht - Printed in the Netherlands The force on a lattice defect in an elastic body R.C. BATRA. " Department of Engineering Mechanics,
More information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More information12. Stresses and Strains
12. Stresses and Strains Finite Element Method Differential Equation Weak Formulation Approximating Functions Weighted Residuals FEM - Formulation Classification of Problems Scalar Vector 1-D T(x) u(x)
More informationThis is a repository copy of The method of fundamental solutions for an inverse boundary value problem in static thermo-elasticity.
This is a repository copy of The method of fundamental solutions for an inverse boundary value problem in static thermo-elasticity. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/816/
More informationConstitutive Relations
Constitutive Relations Dr. Andri Andriyana Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of field
More informationInternational Journal of Pure and Applied Mathematics Volume 58 No ,
International Journal of Pure and Applied Mathematics Volume 58 No. 2 2010, 195-208 A NOTE ON THE LINEARIZED FINITE THEORY OF ELASTICITY Maria Luisa Tonon Department of Mathematics University of Turin
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
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 informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationA MULTISCALE APPROACH IN TOPOLOGY OPTIMIZATION
1 A MULTISCALE APPROACH IN TOPOLOGY OPTIMIZATION Grégoire ALLAIRE CMAP, Ecole Polytechnique The most recent results were obtained in collaboration with F. de Gournay, F. Jouve, O. Pantz, A.-M. Toader.
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationA MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC PROBLEMS IN SEMI-INFINITE MEDIA
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC
More informationSTRESS-RELAXATION IN INCOMPRESSIBLE ELASTIC MATERIALS AT CONSTANT DEFORMATION* R. S. RIVLIN. Brown University
83 STRESS-RELXTION IN INCOMPRESSIBLE ELSTIC MTERILS T CONSTNT DEFORMTION* BY R. S. RIVLIN Brown University 1. Introduction. The mechanics of the finite deformation of incompressible isotropic ideal elastic
More informationDepartment of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm A SIMPLIFIED FORMULATION FOR STRESS AND TRACTION BOUNDARY IN- TEGRAL EQUATIONS USING THE
More informationRegularization methods for large-scale, ill-posed, linear, discrete, inverse problems
Regularization methods for large-scale, ill-posed, linear, discrete, inverse problems Silvia Gazzola Dipartimento di Matematica - Università di Padova January 10, 2012 Seminario ex-studenti 2 Silvia Gazzola
More informationNumerical solution of Cauchy problems in linear elasticity in axisymmetric situations
Numerical solution of Cauchy problems in linear elasticity in axisymmetric situations Bastien Durand, Franck Delvare, Patrice Bailly To cite this version: Bastien Durand, Franck Delvare, Patrice Bailly.
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
More informationReflection of SV- Waves from the Free Surface of a. Magneto-Thermoelastic Isotropic Elastic. Half-Space under Initial Stress
Mathematica Aeterna, Vol. 4, 4, no. 8, 877-93 Reflection of SV- Waves from the Free Surface of a Magneto-Thermoelastic Isotropic Elastic Half-Space under Initial Stress Rajneesh Kakar Faculty of Engineering
More informationContinuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms
Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive
More informationControllability of linear PDEs (I): The wave equation
Controllability of linear PDEs (I): The wave equation M. González-Burgos IMUS, Universidad de Sevilla Doc Course, Course 2, Sevilla, 2018 Contents 1 Introduction. Statement of the problem 2 Distributed
More informationStability of Thick Spherical Shells
Continuum Mech. Thermodyn. (1995) 7: 249-258 Stability of Thick Spherical Shells I-Shih Liu 1 Instituto de Matemática, Universidade Federal do Rio de Janeiro Caixa Postal 68530, Rio de Janeiro 21945-970,
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 informationQUASISTATIC INDENTATION OF A RUBBER COVERED ROLL BY A RIGID ROLL - THE BOUNDARY ELEMENT SOLUTION
259 QUASISTATIC INDENTATION OF A RUBBER COVERED ROLL BY A RIGID ROLL - THE BOUNDARY ELEMENT SOLUTION Batra Department of Engineering Mechanics, Univ. of Missouri-Rolla ABSTRACT The linear elastic problem
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationOn the Rank 1 Convexity of Stored Energy Functions of Physically Linear Stress-Strain Relations
J Elasticity (2007) 86:235 243 DOI 10.1007/s10659-006-9091-z On the Rank 1 Convexity of Stored Energy Functions of Physically Linear Stress-Strain Relations Albrecht Bertram Thomas Böhlke Miroslav Šilhavý
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 informationReduction of Finite Element Models of Complex Mechanical Components
Reduction of Finite Element Models of Complex Mechanical Components Håkan Jakobsson Research Assistant hakan.jakobsson@math.umu.se Mats G. Larson Professor Applied Mathematics mats.larson@math.umu.se Department
More informationDiscrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method
131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationCheckerboard instabilities in topological shape optimization algorithms
Checkerboard instabilities in topological shape optimization algorithms Eric Bonnetier, François Jouve Abstract Checkerboards instabilities for an algorithm of topological design are studied on a simple
More informationPart 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA
Part 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA Review of Fundamentals displacement-strain relation stress-strain relation balance of momentum (deformation) (constitutive equation) (Newton's Law)
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 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 informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
More informationIterative regularization of nonlinear ill-posed problems in Banach space
Iterative regularization of nonlinear ill-posed problems in Banach space Barbara Kaltenbacher, University of Klagenfurt joint work with Bernd Hofmann, Technical University of Chemnitz, Frank Schöpfer and
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 informationENGN 2290: Plasticity Computational plasticity in Abaqus
ENGN 229: Plasticity Computational plasticity in Abaqus The purpose of these exercises is to build a familiarity with using user-material subroutines (UMATs) in Abaqus/Standard. Abaqus/Standard is a finite-element
More informationCellular solid structures with unbounded thermal expansion. Roderic Lakes. Journal of Materials Science Letters, 15, (1996).
1 Cellular solid structures with unbounded thermal expansion Roderic Lakes Journal of Materials Science Letters, 15, 475-477 (1996). Abstract Material microstructures are presented which can exhibit coefficients
More informationLecture Notes 3
12.005 Lecture Notes 3 Tensors Most physical quantities that are important in continuum mechanics like temperature, force, and stress can be represented by a tensor. Temperature can be specified by stating
More informationMECHANICS OF MATERIALS. EQUATIONS AND THEOREMS
1 MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS Version 2011-01-14 Stress tensor Definition of traction vector (1) Cauchy theorem (2) Equilibrium (3) Invariants (4) (5) (6) or, written in terms of principal
More informationA THEORETICAL FORMALISM FOR. Iva Babuska and J.Tinsley Oden. Institute for Computational Engineering and Sciences The University of Texas at Austin
A THEORETICAL FORMALISM FOR VERIFICATION AND VALIDATION Iva Babuska and J.Tinsley Oden Institute for Computational Engineering and Sciences The University of Texas at Austin VERIFICATION AND VALIDATION
More informationCOURSE Numerical methods for solving linear systems. Practical solving of many problems eventually leads to solving linear systems.
COURSE 9 4 Numerical methods for solving linear systems Practical solving of many problems eventually leads to solving linear systems Classification of the methods: - direct methods - with low number of
More informationCopyright 2013 Tech Science Press CMES, vol.94, no.1, pp.1-28, 2013
Copyright 2013 Tech Science Press CMES, vol.94, no.1, pp.1-28, 2013 Application of the MLPG Mixed Collocation Method for Solving Inverse Problems of Linear Isotropic/Anisotropic Elasticity with Simply/Multiply-Connected
More informationThe Non-Linear Field Theories of Mechanics
С. Truesdell-W.Noll The Non-Linear Field Theories of Mechanics Second Edition with 28 Figures Springer-Verlag Berlin Heidelberg NewYork London Paris Tokyo Hong Kong Barcelona Budapest Contents. The Non-Linear
More informationComputer simulation of the Poisson's ratio of soft polydisperse discs at zero temperature
Computer simulation of the Poisson's ratio of soft polydisperse discs at zero temperature JAKUB NAROJCZYK KRZYSZTOF W. WOJCIECHOWSKI Institute of Molecular Physics Polish Academy of Sciences ul. M. Smoluchowskiego
More informationHankel Tranform Method for Solving Axisymmetric Elasticity Problems of Circular Foundation on Semi-infinite Soils
ISSN (Print) : 19-861 ISSN (Online) : 975-44 Charles Chinwuba Ie / International Journal of Engineering and Technology (IJET) Hanel Tranform Method for Solving Axisymmetric Elasticity Problems of Circular
More informationForce Traction Microscopy: an inverse problem with pointwise observations
Force Traction Microscopy: an inverse problem with pointwise observations Guido Vitale Dipartimento di Matematica Politecnico di Torino, Italia Cambridge, December 2011 Force Traction Microscopy: a biophysical
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationIntroduction to Bayesian methods in inverse problems
Introduction to Bayesian methods in inverse problems Ville Kolehmainen 1 1 Department of Applied Physics, University of Eastern Finland, Kuopio, Finland March 4 2013 Manchester, UK. Contents Introduction
More informationRESEARCH ARTICLE. Fast reconstruction of harmonic functions from Cauchy data using the Dirichlet-to-Neumann map and integral equations
Inverse Problems in Science and Engineering Vol. 00, No. 00, September 2010, 1 10 RESEARCH ARTICLE Fast reconstruction of harmonic functions from Cauchy data using the Dirichlet-to-Neumann map and integral
More informationOn the Numerical Modelling of Orthotropic Large Strain Elastoplasticity
63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain
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 informationMathematical Background
CHAPTER ONE Mathematical Background This book assumes a background in the fundamentals of solid mechanics and the mechanical behavior of materials, including elasticity, plasticity, and friction. A previous
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 informationLecture Note 7: Iterative methods for solving linear systems. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 7: Iterative methods for solving linear systems Xiaoqun Zhang Shanghai Jiao Tong University Last updated: December 24, 2014 1.1 Review on linear algebra Norms of vectors and matrices vector
More informationProfessor George C. Johnson. ME185 - Introduction to Continuum Mechanics. Midterm Exam II. ) (1) x
Spring, 997 ME85 - Introduction to Continuum Mechanics Midterm Exam II roblem. (+ points) (a) Let ρ be the mass density, v be the velocity vector, be the Cauchy stress tensor, and b be the body force per
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 informationPEAT SEISMOLOGY Lecture 3: The elastic wave equation
PEAT8002 - SEISMOLOGY Lecture 3: The elastic wave equation Nick Rawlinson Research School of Earth Sciences Australian National University Equation of motion The equation of motion can be derived by considering
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationAn improved convergence theorem for the Newton method under relaxed continuity assumptions
An improved convergence theorem for the Newton method under relaxed continuity assumptions Andrei Dubin ITEP, 117218, BCheremushinsaya 25, Moscow, Russia Abstract In the framewor of the majorization technique,
More informationTransactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN
A computational method for the analysis of viscoelastic structures containing defects G. Ghazlan," C. Petit," S. Caperaa* " Civil Engineering Laboratory, University of Limoges, 19300 Egletons, France &
More informationNumerical Solution of a Coefficient Identification Problem in the Poisson equation
Swiss Federal Institute of Technology Zurich Seminar for Applied Mathematics Department of Mathematics Bachelor Thesis Spring semester 2014 Thomas Haener Numerical Solution of a Coefficient Identification
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 informationNONLINEAR DIFFUSION PDES
NONLINEAR DIFFUSION PDES Erkut Erdem Hacettepe University March 5 th, 0 CONTENTS Perona-Malik Type Nonlinear Diffusion Edge Enhancing Diffusion 5 References 7 PERONA-MALIK TYPE NONLINEAR DIFFUSION The
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationA truly meshless Galerkin method based on a moving least squares quadrature
A truly meshless Galerkin method based on a moving least squares quadrature Marc Duflot, Hung Nguyen-Dang Abstract A new body integration technique is presented and applied to the evaluation of the stiffness
More informationAircraft Structures Kirchhoff-Love Plates
University of Liège erospace & Mechanical Engineering ircraft Structures Kirchhoff-Love Plates Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/ Chemin
More information