Nonsmooth Newton methods for frictional contact problems in flexible multi-body systems
|
|
- Oswald Lloyd
- 5 years ago
- Views:
Transcription
1 1/21 Nonsmooth Newton methods for frictional contact problems in flexible multi-body systems V. Acary, M. Brémond, F. Dubois INRIA Rhône Alpes, Grenoble. LMGC Montpellier 13 eme Colloque National en Calcul des structures May 17, 2017
2 2/21 Objectives Computation methods for discrete frictional contact problems. Rigid multi-body systems: high order of hyperstaticity under-determinancy in the contact forces Projection/splitting methods (Jacobi, Gauss Seidel) are robust but very slow Nonsmooth Newton methods fail Flexible multi-body systems: low order of hyperstaticity Nonsmooth methods work efficiently. General interest in introducing flexibility in the model for computational efficiency.
3 The 3D frictional contact problem 3/21 The 3D frictional contact problem Signorini condition and Coulomb s friction Signorini s condition and Coulomb s friction Body B gap function g N = (C B C A )N. reaction forces r = r NN+r T, with r N IR and r T IR 2. C B Signorini condition at position level g N N C A T 2 T 1 Body A relative velocity 0 g N r N 0. u = u NN+u T, with u N IR and u T IR 2. Signorini condition at velocity level { 0 un r N 0 if g N 0 r N = 0 otherwise.
4 The 3D frictional contact problem 4/21 The 3D frictional contact problem Signorini condition and Coulomb s friction Signorini s condition and Coulomb s friction Modeling assumption Let µ be the coefficient of friction. Let us define the Coulomb friction cone K which is chosen as the isotropic second order cone K = {r IR 3 r T µr n}. (1) The Coulomb friction states for the sticking case that and for the sliding case that u T = 0, r K (2) u T 0, r K, α > 0, r T = αu T. (3) Disjunctive formulation of the frictional contact behavior r = 0 if g N > 0 (no contact) r = 0, u N 0 if g N 0 (take off) r K, u = 0 if g N 0 (sticking) r K, u N = 0, α > 0, u T = αr T if g N 0 (sliding) (4)
5 The 3D frictional contact problem 5/21 The 3D frictional contact problem Signorini condition and Coulomb s friction Signorini s condition and Coulomb s friction Second Order Cone Complementarity (SOCCP) formulation [5] Modified relative velocity û IR 3 defined by û = u + µ u T N. (5) Second-Order Cone Complementarity Problem (SOCCP) K û r K (6) if g N 0 and r = 0 otherwise. The set K is the dual convex cone to K defined by K = {u IR 3 r u 0, for all r K}. (7)
6 The 3D frictional contact problem 6/21 The 3D frictional contact problem Signorini condition and Coulomb s friction Signorini s condition and Coulomb s friction K ut û rn ûn N r K T 1 ûn P rt T 2 ût û Figure: Coulomb s friction and the modified velocity û. The sliding case.
7 The 3D frictional contact problem 7/21 The 3D frictional contact problem 3D frictional contact problems 3D frictional contact problem Multiple contact notation For each contact α {1,... n c}, we have the local velocity : u α IR 3, and the local reaction vector r α IR 3 u = [[u α ], α = 1... n c] r = [[r α ], α = 1... n c] the local Coulomb cone K α = {r α, r α T µα r α N } IR3 and the set K is the cartesian product of Coulomb s friction cone at each contact, that K = K α (8) α=1...n c and K is dual.
8 The 3D frictional contact problem 8/21 The 3D frictional contact problem 3D frictional contact problems 3D frictional contact problems Problem 1 (General discrete frictional contact problem) Given a symmetric positive definite matrix M IR n n, a vector f IR n, a matrix H IR n m, a vector w IR m, a vector of coefficients of friction µ IR nc, find three vectors v IR n, u IR m and r IR m, denoted by FC/I(M, H, f, w, µ) such that Mv = Hr + f u = H v + w û = u + g(u) K û r K (9) with g(u) = [[µ α u α T Nα ], α = 1... n c].
9 The 3D frictional contact problem 9/21 The 3D frictional contact problem 3D frictional contact problems 3D frictional contact problems Problem 2 (Reduced discrete frictional contact problem) Given a symmetric positive semi definite matrix W IR m m, a vector q IR m, a vector µ IR nc of coefficients of friction, find two vectors u IR m and r IR m, denoted by FC/II(W, q, µ) such that u = Wr + q û = u + g(u) K û r K (10) with g(u) = [[µ α u α T Nα ], α = 1... n c]. Relation with the general problem W = H M 1 H and q = H M 1 f + w.
10 The 3D frictional contact problem 10/21 The 3D frictional contact problem 3D frictional contact problems 3D frictional contact problems Rank of the H matrix and hyperstaticity The rank of the H matrix (ratio number of contacts unknows/number of d.o.f) plays an important role. Rigid multibody systems (high degree of hyperstaticity): Generically : 3n c n. H is NOT full column rank and W is rank deficient. Flexible multibody systems. Generically : 3n c < n. H may be full column rank and W is full rank. Effect on convergence of numerical methods First order iterative methods solves all the problems but very slowly Nonsmooth Newton methods are inefficient.
11 Numerical solution procedure. Nonsmooth Equations based methods Nonsmooth Equations based methods Nonsmooth Newton on F (r) = 0 Alart Curnier Formulation [1] F ac(r) := r k+1 = r k Φ 1 (r k )(F (r k )), Φ(r k ) F (r k ) [ rn P IR nc (rn ρn(wr + q)n), + r T P D(µ,(rN ρ(wr+q) N ) +)(r T ρ T(Wr + q) T) Jean Moreau formulation [7, 4] F mj(r) := [ r N P IR nc (rn ρn(wr + q)n) + r T P D(µ,(rN ) +)(r T ρ T(Wr + q) T) ] ], ρ N > 0, ρ T > 0, (11), ρ N > 0, ρ T > 0. (12) Direct natural map reformulation F nat(r) := [ r P K (r ρ(wr + q + g(wr + q))) ], ρ > 0 (13) MUMPS [3, 2] is used for solving linear systems. Numerical solution procedure. 11/21
12 Numerical solution procedure. 12/21 Numerical solution procedure. Matrix block splitting and projection based algorithms Matrix block-splitting and projection based algorithms [9, 8] Block splitting algorithm with W αα IR 3 (Gauss-Seidel) u α i+1 W αα P α i+1 = qα + β<α û α i+1 = [u α N,i+1 + µα u α T,i+1, uα T,i+1 W αβ r β i+1 + W αβ r β i β>α ] T (14) for all α {1... m}. K α, û α i+1 r α i+1 Kα One contact point problem closed form solutions Any solver listed before.
13 Numerical solution procedure. 13/21 Numerical solution procedure. Matrix block splitting and projection based algorithms Naming convention NSN-AC-NLS NSN-JM-NLS NSN-NM-NLS NSN-AC-NLS-HYBRID NSGS-AC NSGS-FP-VI-UPK Nonsmooth Newton Method using (11) without line-search Nonsmooth Newton Method using (12) without line-search Nonsmooth Newton Method using (13) without line-search Method NSN-AC-NLS with preconditioning with 100 iterations of NSGS-AC Gauss Seidel method with NSN-AC-NLS as local solver Gauss Seidel method with fixed point iterations of F nat(r) r Table: Naming convention Error evaluation F nat(r) q < ɛ, (15)
14 Numerical solution procedure. Siconos/Numerics Siconos/Numerics Siconos Open source software for modelling and simulation of nonsmooth systems Siconos/Numerics Collection of C routines to solve FC3D problem NonSmoothGaussSeidel : VI based projection/splitting algorithm TrescaFixedPoint : fixed point algorithm on Tresca fixed point LocalAlartCurnier : semi smooth newton method of Alart Curnier formulation ProximalFixedPoint : proximal point algorithm VIFixedPointProjection : VI based fixed-point projection VIExtragradient : VI based extra-gradient method... use and contribute... Numerical solution procedure. 14/21
15 Preliminary Comparisons 15/21 Preliminary Comparisons Performance profiles Performance profiles [6] Given a set of problems P Given a set of solvers S A performance measure for each problem with a solver t p,s (cpu time, flops,...) Compute the performance ratio τ p,s = tp,s 1 (16) min s S tp,s Compute the performance profile ρ s(τ) : [1, + ] [0, 1] for each solver s S ρ s(τ) = 1 {p P τ p,s τ} (17) P The value of ρ s(1) is the probability that the solver s will win over the rest of the solvers.
16 Preliminary Comparisons 16/21 Preliminary Comparisons Performance profiles LMGC90 sheared low wall example Figure: A low wall meshes with H8 H8 FE with Linear elastic behavior : ρ = 2000kg m 3, E = Pa, ν = 0.2 µ = 0.83 between block and µ = 0.53 between blocks and supports Vertical compression force : 30000N horizontal shear velocity m s 1. Sampling of 50 problems collected in the FCLib with graded difficulty
17 Preliminary Comparisons 17/21 Preliminary Comparisons Performance profiles Results ρ(τ) ρ(τ) τ (flpops) (a) Accuracy ɛ = τ (flpops) (b) Accuracy ɛ = 10 3 NSGS-AC NSGS-FP-VI-UPK iter=100 NSN-AC-NLS NSN-JM-NLS NSN-NM-NLS NSN-AC-NLS-HYBRID
18 Preliminary Comparisons 18/21 Preliminary Comparisons Performance profiles Results ρ(τ) ρ(τ) τ (flpops) (a) Accuracy ɛ = τ (flpops) (b) Accuracy ɛ = 10 6 NSGS-AC NSGS-FP-VI-UPK iter=100 NSN-AC-NLS NSN-JM-NLS NSN-NM-NLS NSN-AC-NLS-HYBRID
19 Conclusions & Perspectives 19/21 Conclusions & Perspectives Conclusions & Perspectives Conclusions 1. For relatively tight accuracy, nonsmooth Newton methods outperform first order iterative method. 2. NSN-AC-NLS-HYBRID is the most efficient method 3. First order iterative methods are interesting for low accuracy, but are not able to reach tight accuracy, Perspectives 1. Evaluate the interest to transform rigid model into flexible ones. 2. Study the possibility to take into account the possible nonlinear bulk behavior in the Newton loop 3. HPC and scalability of nonsmooth Newton techniques using MUMPS 4. Continue to set up a collection of benchmarks FCLIB
20 Conclusions & Perspectives 20/21 Conclusions & Perspectives FCLIB : a collection of discrete 3D Frictional Contact (FC) problems FCLIB : a collection of discrete 3D Frictional Contact (FC) problems Our inspiration: MCPLIB or CUTEst What is FCLIB? A open source collection of Frictional Contact (FC) problems stored in a specific HDF5 format A open source light implementation of Input/Output functions in C Language to read and write problems (Python and Matlab coming soon) Goals of the project Provide a standard framework for testing available and new algorithms for solving discrete frictional contact problems share common formulations of problems in order to exchange data Call for contribution
21 Conclusions & Perspectives 21/21 Conclusions & Perspectives FCLIB : a collection of discrete 3D Frictional Contact (FC) problems Thank you for your attention.
22 Conclusions & Perspectives 21/21 Conclusions & Perspectives FCLIB : a collection of discrete 3D Frictional Contact (FC) problems P. Alart and A. Curnier. A mixed formulation for frictional contact problems prone to Newton like solution method. Computer Methods in Applied Mechanics and Engineering, 92(3): , Patrick R. Amestoy, Iain S. Duff, Jean-Yves L Excellent, and Jacko Koster. A fully asynchronous multifrontal solver using distributed dynamic scheduling. SIAM J. Matrix Anal. Appl., 23(1):15 41 (electronic), Patrick R. Amestoy, Abdou Guermouche, Jean-Yves L Excellent, and Stéphane Pralet. Hybrid scheduling for the parallel solution of linear systems. Parallel Comput., 32(2): , P. Christensen, A. Klarbring, J. Pang, and N. Stromberg. Formulation and comparison of algorithms for frictional contact problems. International Journal for Numerical Methods in Engineering, 42: , G. De Saxcé. Une généralisation de l inégalité de Fenchel et ses applications aux lois constitutives. Comptes Rendus de l Académie des Sciences, t 314,srie II: , E.D. Dolan and J.J. Moré. Benchmarking optimization software with performance profiles. Mathematical Programming, 91(2): , M. Jean and J.J. Moreau. Dynamics in the presence of unilateral contacts and dry friction: a numerical approach. In G. Del Pietro and F. Maceri, editors, Unilateral problems in structural analysis. II, pages CISM 304, Spinger Verlag, F. Jourdan, P. Alart, and M. Jean. A Gauss Seidel like algorithm to solve frictional contact problems. Computer Methods in Applied Mechanics and Engineering, 155(1):31 47, E.N. Mitsopoulou and I.N. Doudoumis.
23 Conclusions & Perspectives 21/21 Conclusions & Perspectives FCLIB : a collection of discrete 3D Frictional Contact (FC) problems A contribution to the analysis of unilateral contact problems with friction. Solid Mechanics Archives, 12(3): , 1987.
On solving frictional contact problems: formulations and comparisons of numerical methods.
On solving frictional contact problems: formulations and comparisons of numerical methods. Vincent Acary and Maurice Brémond and Olivier Huber Abstract In this chapter, we review several formulations of
More information3D revolute joint with clearance in multibody systems
3D revolute joint with clearance in multibody systems Narendra Akhadkar 1, Vincent Acary 2, and Bernard Brogliato 2 1 Schneider Electric, 37 Quai Paul-Louis Merlin, 38000 Grenoble, France, e-mail: narendra.akhadkar@schneider-electric.com
More informationCOMPARISON AND COUPLING OF ALGORITHMS FOR COLLISIONS, CONTACT AND FRICTION IN RIGID MULTI-BODY SIMULATIONS
III European Conference on Computational Mechanics Solids, Structures and Coupled Problems in Engineering C.A. Mota Soares et.al. (eds.) Lisbon, Portugal, 5 8 June 2006 COMPARISON AND COUPLING OF ALGORITHMS
More informationOn solving contact problems with Coulomb friction: formulations and numerical comparisons
On solving contact problems with Coulomb friction: formulations and numerical comparisons Vincent Acary, Maurice Brémond, Olivier Huber To cite this version: Vincent Acary, Maurice Brémond, Olivier Huber.
More informationIterative Rigid Multibody Dynamics A Comparison of Computational Methods
Iterative Rigid Multibody Dynamics A Comparison of Computational Methods Tobias Preclik, Klaus Iglberger, Ulrich Rüde University Erlangen-Nuremberg Chair for System Simulation (LSS) July 1st 2009 T. Preclik
More informationMechanical Simulation Of The ExoMars Rover Using Siconos in 3DROV
V. Acary, M. Brémond, J. Michalczyk, K. Kapellos, R. Pissard-Gibollet 1/15 Mechanical Simulation Of The ExoMars Rover Using Siconos in 3DROV V. Acary, M. Brémond, J. Michalczyk, K. Kapellos, R. Pissard-Gibollet
More informationDistributed Nonsmooth Contact Domain Decomposition (NSCDD): algorithmic structure and scalability
Distributed Nonsmooth Contact Domain Decomposition (NSCDD): algorithmic structure and scalability V. Visseq 1, A. Martin 2, D. Dureisseix 3, F. Dubois 1, and P. Alart 1 1 Introduction Numerical simulations
More informationModeling Fiber Assemblies with Frictional Contact
Modeling Fiber Assemblies with Frictional Contact and application to hair simulation Florence Bertails-Descoubes, BiPop - LJK Research Area Modeling and Simulating complex mechanical objects, and application
More informationA multibody dynamics model of bacterial
A multibody dynamics model of bacterial Martin Servin1 1 UMT Research Lab - Department of Physics Umea University August 26, 2015 (1 : 17) Nonsmooth multidomain dynamics What has this to do with bacterias?
More informationQuartic formulation of Coulomb 3D frictional contact
INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE Quartic formulation of Coulomb 3D frictional contact Olivier Bonnefon Gilles Daviet N 0400 January 10, 2011 Modeling, Optimization, and
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 informationSolution of contact problems with friction
Solution of contact problems with friction Pierre Joli, Zhi-Qiang Feng, Nadjet Talbi To cite this version: Pierre Joli, Zhi-Qiang Feng, Nadjet Talbi. Solution of contact problems with friction: Pseudo-potential
More informationFrom Direct to Iterative Substructuring: some Parallel Experiences in 2 and 3D
From Direct to Iterative Substructuring: some Parallel Experiences in 2 and 3D Luc Giraud N7-IRIT, Toulouse MUMPS Day October 24, 2006, ENS-INRIA, Lyon, France Outline 1 General Framework 2 The direct
More informationNonsmooth contact mechanics. Spring school. Aussois 2010
1/43 Nonsmooth contact mechanics. Spring school. Aussois 2010 Formulation of nonsmooth dynamical systems Measure differential inclusions Basics on Mathematical properties Formulation of unilateral contact,
More informationA distributed memory parallel multibody Contact Dynamics code
A distributed memory parallel multibody Contact Dynamics code Tomasz Koziara, Nenad Bićanić Department of Civil Engineering University of Glasgow, Glasgow G12 8LT, UK Abstract Solfec is a computational
More informationOn the Relation of the Principle of Maximum Dissipation to the Principle of Gauss
On the Relation of the Principle of Maximum Dissipation to the Principle of Gauss Kerim Yunt P.O Box 7 Zurich, Switzerland 82 Email: kerimyunt@web.de is about conditions which specify when sticking and
More informationUniversity of Erlangen-Nürnberg and Academy of Sciences of the Czech Republic. Solving MPECs by Implicit Programming and NLP Methods p.
Solving MPECs by Implicit Programming and NLP Methods Michal Kočvara University of Erlangen-Nürnberg and Academy of Sciences of the Czech Republic Solving MPECs by Implicit Programming and NLP Methods
More informationTime integration. DVI and HHT time stepping methods in Chrono
Time integration DVI and HHT time stepping methods in Chrono Time Integration in Chrono Two classes of time stepping methods in Chrono Time steppers for smooth dynamics Classical multibody dynamics rigid
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 informationA FRICTION MODEL FOR NON-SINGULAR COMPLEMENTARITY FORMULATIONS FOR MULTIBODY SYSTEMS WITH CONTACTS
Proceedings of the ASME 217 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference IDETC217 August 6-9, 217, Cleveland, Ohio, USA DETC217-67988 A FRICTION
More informationLehrstuhl B für Mechanik Technische Universität München D Garching Germany
DISPLACEMENT POTENTIALS IN NON-SMOOTH DYNAMICS CH. GLOCKER Lehrstuhl B für Mechanik Technische Universität München D-85747 Garching Germany Abstract. The paper treats the evaluation of the accelerations
More informationTowards an augmented domain decomposition method for nonsmooth contact dynamics models
Comp. Part. Mech. (214 1:15 26 DOI 1.17/s4571-14-5-8 Towards an augmented domain decomposition method for nonsmooth contact dynamics models Vincent Visseq Pierre Alart David Dureisseix Received: 24 December
More informationMulti-source inversion of TEM data: with field applications to Mt. Milligan
Multi-source inversion of TEM data: with field applications to Mt. Milligan D. Oldenburg, E. Haber, D. Yang Geophysical Inversion Facility, University of British Columbia, Vancouver, BC, Canada Summary
More informationProximal-like contraction methods for monotone variational inequalities in a unified framework
Proximal-like contraction methods for monotone variational inequalities in a unified framework Bingsheng He 1 Li-Zhi Liao 2 Xiang Wang Department of Mathematics, Nanjing University, Nanjing, 210093, China
More informationApplying a type of SOC-functions to solve a system of equalities and inequalities under the order induced by second-order cone
Applying a type of SOC-functions to solve a system of equalities and inequalities under the order induced by second-order cone Xin-He Miao 1, Nuo Qi 2, B. Saheya 3 and Jein-Shan Chen 4 Abstract: In this
More informationv Fc v < 3 v < 2 v < 1 v
Fc v v Fc 2 4 3 1 4 v< 3 v< 2 v< 1 v v Fc ( ) v 0 8 2 F c Fc v v Fc rv 8r >0 v Fc rv 8r >0 = prox C ( rv) My Favorite Benefits Local non-linear friction models used directly without any discretization
More informationA new algorithm for solving 3D contact problems with Coulomb friction
A new algorithm for solving 3D contact problems with Coulomb friction Radek Kučera VŠB-TU Ostrava, Department of Mathematics and Descriptive Geometry, name@email.address Summary. The paper deals with solving
More informationA High-Performance Parallel Hybrid Method for Large Sparse Linear Systems
Outline A High-Performance Parallel Hybrid Method for Large Sparse Linear Systems Azzam Haidar CERFACS, Toulouse joint work with Luc Giraud (N7-IRIT, France) and Layne Watson (Virginia Polytechnic Institute,
More informationAN OPTIMIZATION-BASED ALGORITHM FOR COULOMB S FRICTIONAL CONTACT. Florent Cadoux 1
ESAIM: PROCEEDINGS, May 2009, Vol. 27, p. 54-69 C. Besse, O. Goubet, T. Goudon & S. Nicaise, Editors AN OPTIMIZATION-BASED ALGORITHM FOR COULOMB S FRICTIONAL CONTACT Florent Cadoux 1 Abstract. The main
More information27. Impact Mechanics of Manipulation
27. Impact Mechanics of Manipulation Matt Mason matt.mason@cs.cmu.edu http://www.cs.cmu.edu/~mason Carnegie Mellon Lecture 27. Mechanics of Manipulation p.1 Lecture 27. Impact Chapter 1 Manipulation 1
More information2D OR 3D FRICTIONAL CONTACT ALGORITHMS AND APPLICATIONS ln A LARGE DEFORMA TION CONTEXT
~~". " COMMUNICAllONS ln NUMERICAL METHODS ln ENGINEERING, Vol. Il,409-416 (1995) 2D OR 3D FRICTIONAL CONTACT ALGORITHMS AND APPLICATIONS ln A LARGE DEFORMA TION CONTEXT ZIll-QIANG FENG Polytechnic Institute
More informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
More informationLIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS
5 th GRACM International Congress on Computational Mechanics Limassol, 29 June 1 July, 2005 LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS G.A. Drosopoulos I, G.E.
More informationSpectral gradient projection method for solving nonlinear monotone equations
Journal of Computational and Applied Mathematics 196 (2006) 478 484 www.elsevier.com/locate/cam Spectral gradient projection method for solving nonlinear monotone equations Li Zhang, Weijun Zhou Department
More informationInstabilities and Dynamic Rupture in a Frictional Interface
Instabilities and Dynamic Rupture in a Frictional Interface Laurent BAILLET LGIT (Laboratoire de Géophysique Interne et Tectonophysique) Grenoble France laurent.baillet@ujf-grenoble.fr http://www-lgit.obs.ujf-grenoble.fr/users/lbaillet/
More informationMixed Shooting-HBM: a periodic solution solver for unilaterally constrained systems
The 4 th Joint International Conference on Multibody System Dynamics May 29 June 2, 216, Montréal, Canada Mixed Shooting-HBM: a periodic solution solver for unilaterally constrained systems Frederic Schreyer
More informationNumerical simulation of monuments by the contact dynamics method
Numerical simulation of monuments by the contact dynamics method Vincent Acary, Michel Jean To cite this version: Vincent Acary, Michel Jean. Numerical simulation of monuments by the contact dynamics method.
More informationParallel scalability of a FETI DP mortar method for problems with discontinuous coefficients
Parallel scalability of a FETI DP mortar method for problems with discontinuous coefficients Nina Dokeva and Wlodek Proskurowski University of Southern California, Department of Mathematics Los Angeles,
More informationROCKING RESPONSE OF MASONRY BLOCK STRUCTURES USING MATHEMATICAL PROGRAMMING
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationIntroduction to Heat Transfer Analysis
Introduction to Heat Transfer Analysis Thermal Network Solutions with TNSolver Bob Cochran Applied Computational Heat Transfer Seattle, WA TNSolver@heattransfer.org ME 331 Introduction to Heat Transfer
More informationBlock Low-Rank (BLR) approximations to improve multifrontal sparse solvers
Block Low-Rank (BLR) approximations to improve multifrontal sparse solvers Joint work with Patrick Amestoy, Cleve Ashcraft, Olivier Boiteau, Alfredo Buttari and Jean-Yves L Excellent, PhD started on October
More informationNonsmooth Schur Newton Methods and Applications
Nonsmooth Schur Newton Methods and Applications C. Gräser, R. Kornhuber, and U. Sack IMA Annual Program Year Workshop Numerical Solutions of Partial Differential Equations: Fast Solution Techniques November
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 09-08 Preconditioned conjugate gradient method enhanced by deflation of rigid body modes applied to composite materials. T.B. Jönsthövel ISSN 1389-6520 Reports of
More informationA Space-Time Multigrid Solver Methodology for the Optimal Control of Time-Dependent Fluid Flow
A Space-Time Multigrid Solver Methodology for the Optimal Control of Time-Dependent Fluid Flow Michael Köster, Michael Hinze, Stefan Turek Michael Köster Institute for Applied Mathematics TU Dortmund Trier,
More informationEfficient BDF time discretization of the Navier Stokes equations with VMS LES modeling in a High Performance Computing framework
Introduction VS LES/BDF Implementation Results Conclusions Swiss Numerical Analysis Day 2015, Geneva, Switzerland, 17 April 2015. Efficient BDF time discretization of the Navier Stokes equations with VS
More informationFast convergent finite difference solvers for the elliptic Monge-Ampère equation
Fast convergent finite difference solvers for the elliptic Monge-Ampère equation Adam Oberman Simon Fraser University BIRS February 17, 211 Joint work [O.] 28. Convergent scheme in two dim. Explicit solver.
More informationNumerical Modeling of Methane Hydrate Evolution
Numerical Modeling of Methane Hydrate Evolution Nathan L. Gibson Joint work with F. P. Medina, M. Peszynska, R. E. Showalter Department of Mathematics SIAM Annual Meeting 2013 Friday, July 12 This work
More informationMUMPS. The MUMPS library: work done during the SOLSTICE project. MUMPS team, Lyon-Grenoble, Toulouse, Bordeaux
The MUMPS library: work done during the SOLSTICE project MUMPS team, Lyon-Grenoble, Toulouse, Bordeaux Sparse Days and ANR SOLSTICE Final Workshop June MUMPS MUMPS Team since beg. of SOLSTICE (2007) Permanent
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationParallel Scalability of a FETI DP Mortar Method for Problems with Discontinuous Coefficients
Parallel Scalability of a FETI DP Mortar Method for Problems with Discontinuous Coefficients Nina Dokeva and Wlodek Proskurowski Department of Mathematics, University of Southern California, Los Angeles,
More informationA Generalized Maximum Dissipation Principle in an Impulse-velocity Time-stepping Scheme
A Generalized Maximum Dissipation Principle in an Impulse-velocity Time-stepping Scheme T. Preclik, U. Rüde September 2, 213 Chair of Computer Science 1 (System Simulation) University of Erlangen-Nürnberg,
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
More informationA Novel Inexact Smoothing Method for Second-Order Cone Complementarity Problems
A Novel Inexact Smoothing Method for Second-Order Cone Complementarity Problems Xiaoni Chi Guilin University of Electronic Technology School of Math & Comput Science Guilin Guangxi 541004 CHINA chixiaoni@126.com
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationPartitioned Methods for Multifield Problems
Partitioned Methods for Multifield Problems Joachim Rang, 10.5.2017 10.5.2017 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 Contents Blockform of linear iteration schemes Examples 10.5.2017
More informationGlobally convergent three-term conjugate gradient projection methods for solving nonlinear monotone equations
Arab. J. Math. (2018) 7:289 301 https://doi.org/10.1007/s40065-018-0206-8 Arabian Journal of Mathematics Mompati S. Koorapetse P. Kaelo Globally convergent three-term conjugate gradient projection methods
More informationIMES - Center of Mechanics ETH Zentrum, Tannenstrasse 3 CH-8092 Zürich, Switzerland
Chapter 16 THE GEOMETRY OF NEWTON S CRADLE Christoph Glocker and Ueli Aeberhard IMES - Center of Mechanics ETH Zentrum, Tannenstrasse 3 CH-8092 Zürich, Switzerland glocker@imes.mavt.ethz.ch Abstract A
More informationSome Inexact Hybrid Proximal Augmented Lagrangian Algorithms
Some Inexact Hybrid Proximal Augmented Lagrangian Algorithms Carlos Humes Jr. a, Benar F. Svaiter b, Paulo J. S. Silva a, a Dept. of Computer Science, University of São Paulo, Brazil Email: {humes,rsilva}@ime.usp.br
More informationVisual Interactive Simulation, TDBD24, Spring 2006
Visual Interactive Simulation, TDBD24, Spring 2006 Lecture 6 Outline Cross product as matrix multiplication The collision matrix The Newton Coulomb impulse for rigid bodies N C updates of translational
More informationMultigrid based preconditioners for the numerical solution of two-dimensional heterogeneous problems in geophysics
Technical Report RAL-TR-2007-002 Multigrid based preconditioners for the numerical solution of two-dimensional heterogeneous problems in geophysics I. S. Duff, S. Gratton, X. Pinel, and X. Vasseur January
More informationPath-following and augmented Lagrangian methods for contact problems in linear elasticity
Journal of Computational and Applied Mathematics 203 (2007) 533 547 www.elsevier.com/locate/cam Path-following and augmented Lagrangian methods for contact problems in linear elasticity Georg Stadler Center
More informationSOR- and Jacobi-type Iterative Methods for Solving l 1 -l 2 Problems by Way of Fenchel Duality 1
SOR- and Jacobi-type Iterative Methods for Solving l 1 -l 2 Problems by Way of Fenchel Duality 1 Masao Fukushima 2 July 17 2010; revised February 4 2011 Abstract We present an SOR-type algorithm and a
More informationSMO vs PDCO for SVM: Sequential Minimal Optimization vs Primal-Dual interior method for Convex Objectives for Support Vector Machines
vs for SVM: Sequential Minimal Optimization vs Primal-Dual interior method for Convex Objectives for Support Vector Machines Ding Ma Michael Saunders Working paper, January 5 Introduction In machine learning,
More informationA formulation of the linear discrete Coulomb friction problem via convex optimization
A formulation of the linear discrete Coulomb friction problem via convex optimization Vincent Acary, Florent Cadoux, Claude Lemaréchal, Jérôme Malick To cite this version: Vincent Acary, Florent Cadoux,
More informationInverse Problems and Optimal Design in Electricity and Magnetism
Inverse Problems and Optimal Design in Electricity and Magnetism P. Neittaanmäki Department of Mathematics, University of Jyväskylä M. Rudnicki Institute of Electrical Engineering, Warsaw and A. Savini
More informationDissipative Systems Analysis and Control, Theory and Applications: Addendum/Erratum
Dissipative Systems Analysis and Control, Theory and Applications: Addendum/Erratum Bernard Brogliato To cite this version: Bernard Brogliato. Dissipative Systems Analysis and Control, Theory and Applications:
More informationManual of ReSNA. matlab software for mixed nonlinear second-order cone complementarity problems based on Regularized Smoothing Newton Algorithm
Manual of ReSNA matlab software for mixed nonlinear second-order cone complementarity problems based on Regularized Smoothing Newton Algorithm Shunsuke Hayashi September 4, 2013 1 Introduction ReSNA (Regularized
More informationA Hybrid Iterative Solver for Robustly Capturing Coulomb Friction in Hair Dynamics
A Hybrid Iterative Solver for Robustly Capturing Coulomb Friction in Hair Dynamics Gilles Daviet Florence Bertails-Descoubes Laurence Boissieux INRIA, Grenoble, France Figure : Comparison of the hair collective
More informationPreconditioning the Newton-Krylov method for reactive transport
Preconditioning the Newton-Krylov method for reactive transport Michel Kern, L. Amir, A. Taakili INRIA Paris-Rocquencourt Maison de la Simulation MoMaS workshop Reactive Transport Modeling in the Geological
More informationTHE solution of the absolute value equation (AVE) of
The nonlinear HSS-like iterative method for absolute value equations Mu-Zheng Zhu Member, IAENG, and Ya-E Qi arxiv:1403.7013v4 [math.na] 2 Jan 2018 Abstract Salkuyeh proposed the Picard-HSS iteration method
More informationCHAPTER 10: Numerical Methods for DAEs
CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct
More informationSolution of contact problems in linear elasticity using a feasible interior point algorithm for nonlinear complementarity problems
7 th World Congress on Structural and Multidisciplinary Optimization COEX Seoul, May - 5 May 007, Korea Solution of contact problems in linear elasticity using a feasible interior point algorithm for nonlinear
More informationProgress in Parallel Implicit Methods For Tokamak Edge Plasma Modeling
Progress in Parallel Implicit Methods For Tokamak Edge Plasma Modeling Michael McCourt 1,2,Lois Curfman McInnes 1 Hong Zhang 1,Ben Dudson 3,Sean Farley 1,4 Tom Rognlien 5, Maxim Umansky 5 Argonne National
More informationSeventh Symposium of the European Network for Nonsmooth Dynamics (ENNSD) September 12-13, 2018, University of Stuttgart, Germany
Seventh Symposium of the European Network for Nonsmooth Dynamics (ENNSD) September 12-13, 2018, University of Stuttgart, Germany Organized by: Remco Leine, Vincent Acary, and Olivier Brüls This is the
More informationA derivative-free nonmonotone line search and its application to the spectral residual method
IMA Journal of Numerical Analysis (2009) 29, 814 825 doi:10.1093/imanum/drn019 Advance Access publication on November 14, 2008 A derivative-free nonmonotone line search and its application to the spectral
More informationSolving MPECs Implicit Programming and NLP Methods
Solving MPECs Implicit Programming and NLP Methods Michal Kočvara Academy of Sciences of the Czech Republic September 2005 1 Mathematical Programs with Equilibrium Constraints Mechanical motivation Mechanical
More informationScalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems
Scalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems Pierre Jolivet, F. Hecht, F. Nataf, C. Prud homme Laboratoire Jacques-Louis Lions Laboratoire Jean Kuntzmann INRIA Rocquencourt
More informationAdaptive discretization and first-order methods for nonsmooth inverse problems for PDEs
Adaptive discretization and first-order methods for nonsmooth inverse problems for PDEs Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Barbara Kaltenbacher, Tuomo Valkonen,
More informationSIMULATION OF INDETERMINATE MULTI-POINT IMPACT AND CONTACT WITH FRICTION
MULTIBODY DYNAMICS 011, ECCOMAS Thematic Conference J.C. Samin, P. Fisette (eds.) Brussels, Belgium, 4-7 July 011 SIMULATION OF INDETERMINATE MULTI-POINT IMPACT AND CONTACT WITH FRICTION Adrian Rodriguez
More informationIn collaboration with J.-C. Pesquet A. Repetti EC (UPE) IFPEN 16 Dec / 29
A Random block-coordinate primal-dual proximal algorithm with application to 3D mesh denoising Emilie CHOUZENOUX Laboratoire d Informatique Gaspard Monge - CNRS Univ. Paris-Est, France Horizon Maths 2014
More informationA Recursive Trust-Region Method for Non-Convex Constrained Minimization
A Recursive Trust-Region Method for Non-Convex Constrained Minimization Christian Groß 1 and Rolf Krause 1 Institute for Numerical Simulation, University of Bonn. {gross,krause}@ins.uni-bonn.de 1 Introduction
More informationQuadrature based Broyden-like method for systems of nonlinear equations
STATISTICS, OPTIMIZATION AND INFORMATION COMPUTING Stat., Optim. Inf. Comput., Vol. 6, March 2018, pp 130 138. Published online in International Academic Press (www.iapress.org) Quadrature based Broyden-like
More informationOptimal Preconditioners for Interval Gauss Seidel Methods
Optimal Preconditioners for Interval Gauss Seidel Methods R B Kearfott and Xiaofa Shi Introduction Consider the following nonlinear system f 1 x 1, x 2,, x n ) F X) f n x 1, x 2,, x n ) where bounds x
More informationInfeasibility Detection and an Inexact Active-Set Method for Large-Scale Nonlinear Optimization
Infeasibility Detection and an Inexact Active-Set Method for Large-Scale Nonlinear Optimization Frank E. Curtis, Lehigh University involving joint work with James V. Burke, University of Washington Daniel
More informationCAEFEM v9.5 Information
CAEFEM v9.5 Information Concurrent Analysis Corporation, 50 Via Ricardo, Thousand Oaks, CA 91320 USA Tel. (805) 375 1060, Fax (805) 375 1061 email: info@caefem.com or support@caefem.com Web: http://www.caefem.com
More informationNumerical Modeling of granular materials with multi physics coupling. Application to the analysis of granular flow
Numerical Modeling of granular materials with multi physics coupling. Application to the analysis of granular flow Fre de ric Dubois1 and Mathieu Renouf2 1 Laboratoire de Me canique et Ge nie Civil Universite
More informationMultilevel spectral coarse space methods in FreeFem++ on parallel architectures
Multilevel spectral coarse space methods in FreeFem++ on parallel architectures Pierre Jolivet Laboratoire Jacques-Louis Lions Laboratoire Jean Kuntzmann DD 21, Rennes. June 29, 2012 In collaboration with
More informationCONCURRENT MULTIPLE IMPACTS IN RIGID BODIES: FORMULATION AND SIMULATION.
ENOC-2005, Eindhoven, Netherlands, 7-2 August 2005 CONCURRENT MULTIPLE IMPACTS IN RIGID BODIES: FORMULATION AND SIMULATION. Vincent Acary Doh Elvis Taha BIPOP Project INRIA Rhône-Alpes, Grenoble, France
More informationIDENTIFICATION OF SEMI-RIGID JOINTS IN FRAME STRUCTURES
6th ECCOMAS Conference on Smart Structures and Materials SMART213 Politecnico di Torino, 24-26 June 213 E. Carrera, F. Miglioretti and M. Petrolo (Editors) www.smart213.com IDENTIFICATION OF SEMI-RIGID
More informationOptimization Algorithms for Compressed Sensing
Optimization Algorithms for Compressed Sensing Stephen Wright University of Wisconsin-Madison SIAM Gator Student Conference, Gainesville, March 2009 Stephen Wright (UW-Madison) Optimization and Compressed
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More informationA Hard Constraint Time-Stepping Approach for Rigid Multibody Dynamics with Joints, Contact and Friction. Gary D. Hart University of Pittsburgh
A Hard Constraint Time-Stepping Approach for Rigid Multibody Dynamics with Joints, Contact and Friction Gary D. Hart University of Pittsburgh 1 Goal Develop method for simulating rigid multibody dynamics
More informationDiscrete Multi-material Topology Optimization under Total Mass Constraint
Discrete Multi-material Topology Optimization under Total Mass Constraint Xingtong Yang Ming Li State Key Laboratory of CAD&CG, Zhejiang University Solid and Physical Modeling Bilbao, June 12, 2018 Problem
More informationParallelizing large scale time domain electromagnetic inverse problem
Parallelizing large scale time domain electromagnetic inverse problems Eldad Haber with: D. Oldenburg & R. Shekhtman + Emory University, Atlanta, GA + The University of British Columbia, Vancouver, BC,
More informationSCALAR FORCE POTENTIALS IN RIGID MULTIBODY SYSTEMS. Ch. Glocker Technical University of Munich, Germany
SCALAR FORCE POTENTIALS IN RIGID MULTIBODY SYSTEMS Ch. Glocker Technical University of Munich, Germany ABSTRACT These lectures treat the motion of finite-dimensional mechanical systems under the influence
More informationConvex Optimization. (EE227A: UC Berkeley) Lecture 15. Suvrit Sra. (Gradient methods III) 12 March, 2013
Convex Optimization (EE227A: UC Berkeley) Lecture 15 (Gradient methods III) 12 March, 2013 Suvrit Sra Optimal gradient methods 2 / 27 Optimal gradient methods We saw following efficiency estimates for
More informationToward black-box adaptive domain decomposition methods
Toward black-box adaptive domain decomposition methods Frédéric Nataf Laboratory J.L. Lions (LJLL), CNRS, Alpines Inria and Univ. Paris VI joint work with Victorita Dolean (Univ. Nice Sophia-Antipolis)
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationA Power Penalty Method for a Bounded Nonlinear Complementarity Problem
A Power Penalty Method for a Bounded Nonlinear Complementarity Problem Song Wang 1 and Xiaoqi Yang 2 1 Department of Mathematics & Statistics Curtin University, GPO Box U1987, Perth WA 6845, Australia
More information