648 research outputs found
A Variational r-Adaption and Shape-Optimization Method for Finite-Deformation Elasticity
This paper is concerned with the formulation of a variational r-adaption method for finite-deformation elastostatic problems. The distinguishing characteristic of the method is that the variational principle simultaneously supplies the solution, the optimal mesh and, in problems of shape optimization, the equilibrium shapes of the system. This is accomplished by minimizing the energy functional with respect to the nodal field values as well as with respect to the triangulation of the domain of analysis. Energy minimization with respect to the referential nodal positions has the effect of equilibrating the energetic or configurational forces acting on the nodes. We derive general expressions for the configuration forces for isoparametric elements and nonlinear, possibly anisotropic, materials under general loading. We illustrate the versatility and convergence characteristics of the method by way of selected numerical tests and applications, including the problem of a semi-infinite crack in linear and nonlinear elastic bodies; and the optimization of the shape of elastic inclusions
The role of the patch test in 2D atomistic-to-continuum coupling methods
For a general class of atomistic-to-continuum coupling methods, coupling
multi-body interatomic potentials with a P1-finite element discretisation of
Cauchy--Born nonlinear elasticity, this paper adresses the question whether
patch test consistency (or, absence of ghost forces) implies a first-order
error estimate.
In two dimensions it is shown that this is indeed true under the following
additional technical assumptions: (i) an energy consistency condition, (ii)
locality of the interface correction, (iii) volumetric scaling of the interface
correction, and (iv) connectedness of the atomistic region. The extent to which
these assumptions are necessary is discussed in detail.Comment: Version 2: correction of some minor mistakes, added discussion of
multiple connected atomistic region, minor improvements of styl
Optimal design of plane elastic membranes using the convexified F\"{o}ppl's model
This work puts forth a new optimal design formulation for planar elastic
membranes. The goal is to minimize the membrane's compliance through choosing
the material distribution described by a positive Radon measure. The
deformation of the membrane itself is governed by the convexified F\"{o}ppl's
model. The uniqueness of this model lies in the convexity of its variational
formulation despite the inherent nonlinearity of the strain-displacement
relation. It makes it possible to rewrite the optimization problem as a pair of
mutually dual convex variational problems. In the primal problem a linear
functional is maximized with respect to displacement functions while enforcing
that point-wisely the strain lies in an unbounded closed convex set. The dual
problem consists in finding equilibrated stresses that are to minimize a convex
integral functional of linear growth defined on the space of Radon measures.
The pair of problems is analysed: existence and regularity results are
provided, together with the system of optimality criteria. To demonstrate the
computational potential of the pair, a finite element scheme is developed
around it. Upon reformulation to a conic-quadratic & semi-definite programming
problem, the method is employed to produce numerical simulations for several
load case scenarios.Comment: 55 page
Metric based up-scaling
We consider divergence form elliptic operators in dimension with
coefficients. Although solutions of these operators are only
H\"{o}lder continuous, we show that they are differentiable ()
with respect to harmonic coordinates. It follows that numerical homogenization
can be extended to situations where the medium has no ergodicity at small
scales and is characterized by a continuum of scales by transferring a new
metric in addition to traditional averaged (homogenized) quantities from
subgrid scales into computational scales and error bounds can be given. This
numerical homogenization method can also be used as a compression tool for
differential operators.Comment: Final version. Accepted for publication in Communications on Pure and
Applied Mathematics. Presented at CIMMS (March 2005), Socams 2005 (April),
Oberwolfach, MPI Leipzig (May 2005), CIRM (July 2005). Higher resolution
figures are available at http://www.acm.caltech.edu/~owhadi
Aspects of Unstructured Grids and Finite-Volume Solvers for the Euler and Navier-Stokes Equations
One of the major achievements in engineering science has been the development of computer algorithms for solving nonlinear differential equations such as the Navier-Stokes equations. In the past, limited computer resources have motivated the development of efficient numerical schemes in computational fluid dynamics (CFD) utilizing structured meshes. The use of structured meshes greatly simplifies the implementation of CFD algorithms on conventional computers. Unstructured grids on the other hand offer an alternative to modeling complex geometries. Unstructured meshes have irregular connectivity and usually contain combinations of triangles, quadrilaterals, tetrahedra, and hexahedra. The generation and use of unstructured grids poses new challenges in CFD. The purpose of this note is to present recent developments in the unstructured grid generation and flow solution technology
Variational Methods for Evolution
The meeting focused on the last advances in the applications of variational methods to evolution problems governed by partial differential equations. The talks covered a broad range of topics, including large deviation and variational principles, rate-independent evolutions and gradient flows, heat flows in metric-measure spaces, propagation of fracture, applications of optimal transport and entropy-entropy dissipation methods, phase-transitions, viscous approximation, and singular-perturbation problems
Numerical Methods for the Modelling of Chip Formation
The modeling of metal cutting has proved to be particularly complex due to the diversity of physical phenomena involved, including thermo-mechanical coupling, contact/friction and material failure. During the last few decades, there has been significant progress in the development of numerical methods for modeling machining operations. Furthermore, the most relevant techniques have been implemented in the relevant commercial codes creating tools for the engineers working in the design of processes and cutting devices. This paper presents a review on the numerical modeling methods and techniques used for the simulation of machining processes. The main purpose is to identify the strengths and weaknesses of each method and strategy developed up-to-now. Moreover the review covers the classical Finite Element Method covering mesh-less methods, particle-based methods and different possibilities of Eulerian and Lagrangian approaches
Unstructured Grid Generation Techniques and Software
The Workshop on Unstructured Grid Generation Techniques and Software was conducted for NASA to assess its unstructured grid activities, improve the coordination among NASA centers, and promote technology transfer to industry. The proceedings represent contributions from Ames, Langley, and Lewis Research Centers, and the Johnson and Marshall Space Flight Centers. This report is a compilation of the presentations made at the workshop
A Dimension-Reduction Model for Brittle Fractures on Thin Shells with Mesh Adaptivity
In this paper we derive a new two-dimensional brittle fracture model for thin
shells via dimension reduction, where the admissible displacements are only
normal to the shell surface. The main steps include to endow the shell with a
small thickness, to express the three-dimensional energy in terms of the
variational model of brittle fracture in linear elasticity, and to study the
-limit of the functional as the thickness tends to zero. The numerical
discretization is tackled by first approximating the fracture through a phase
field, following an Ambrosio-Tortorelli like approach, and then resorting to an
alternating minimization procedure, where the irreversibility of the crack
propagation is rigorously imposed via an inequality constraint. The
minimization is enriched with an anisotropic mesh adaptation driven by an a
posteriori error estimator, which allows us to sharply track the whole crack
path by optimizing the shape, the size, and the orientation of the mesh
elements. Finally, the overall algorithm is successfully assessed on two
Riemannian settings and proves not to bias the crack propagation
- …