			Home About us Contact

Discontinuous Functions (discontinuous + function)

Distribution by Scientific Domains

Engineering	77%

Selected Abstracts

Improved implementation and robustness study of the X-FEM for stress analysis around cracks

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 8 2005
E. Béchet
Abstract Numerical crack propagation schemes were augmented in an elegant manner by the X-FEM method. The use of special tip enrichment functions, as well as a discontinuous function along the sides of the crack allows one to do a complete crack analysis virtually without modifying the underlying mesh, which is of industrial interest, especially when a numerical model for crack propagation is desired. This paper improves the implementation of the X-FEM method for stress analysis around cracks in three ways. First, the enrichment strategy is revisited. The conventional approach uses a ,topological' enrichment (only the elements touching the front are enriched). We suggest a ,geometrical' enrichment in which a given domain size is enriched. The improvements obtained with this enrichment are discussed. Second, the conditioning of the X-FEM both for topological and geometrical enrichments is studied. A preconditioner is introduced so that ,off the shelf' iterative solver packages can be used and perform as well on X-FEM matrices as on standard FEM matrices. The preconditioner uses a local (nodal) Cholesky based decomposition. Third, the numerical integration scheme to build the X-FEM stiffness matrix is dramatically improved for tip enrichment functions by the use of an ad hoc integration scheme. A 2D benchmark problem is designed to show the improvements and the robustness. Copyright © 2005 John Wiley & Sons, Ltd. [source]

Partition of unity enrichment for bimaterial interface cracks

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 8 2004
N. Sukumar
Abstract Partition of unity enrichment techniques are developed for bimaterial interface cracks. A discontinuous function and the two-dimensional near-tip asymptotic displacement functions are added to the finite element approximation using the framework of partition of unity. This enables the domain to be modelled by finite elements without explicitly meshing the crack surfaces. The crack-tip enrichment functions are chosen as those that span the asymptotic displacement fields for an interfacial crack. The concept of partition of unity facilitates the incorporation of the oscillatory nature of the singularity within a conforming finite element approximation. The mixed-mode (complex) stress intensity factors for bimaterial interfacial cracks are numerically evaluated using the domain form of the interaction integral. Good agreement between the numerical results and the reference solutions for benchmark interfacial crack problems is realized. Copyright © 2004 John Wiley & Sons, Ltd. [source]

A G space theory and a weakened weak (W2) form for a unified formulation of compatible and incompatible methods: Part I theory

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 9 2010
G. R. Liu
Abstract This paper introduces a G space theory and a weakened weak form (W2) using the generalized gradient smoothing technique for a unified formulation of a wide class of compatible and incompatible methods. The W2 formulation works for both finite element method settings and mesh-free settings, and W2 models can have special properties including softened behavior, upper bounds and ultra accuracy. Part I of this paper focuses on the theory and fundamentals for W2 formulations. A normed G space is first defined to include both continuous and discontinuous functions allowing the use of much more types of methods/techniques to create shape functions for numerical models. Important properties and a set of useful inequalities for G spaces are then proven in the theory and analyzed in detail. These properties ensure that a numerical method developed based on the W2 formulation will be spatially stable and convergent to the exact solutions, as long as the physical problem is well posed. The theory is applicable to any problems to which the standard weak formulation is applicable, and can offer numerical solutions with special properties including ,close-to-exact' stiffness, upper bounds and ultra accuracy. Copyright © 2009 John Wiley & Sons, Ltd. [source]

A G space theory and a weakened weak (W2) form for a unified formulation of compatible and incompatible methods: Part II applications to solid mechanics problems

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 9 2010
G. R. Liu
Abstract In part I of this paper, we have established the G space theory and fundamentals for W2 formulation. Part II focuses on the applications of the G space theory to formulate W2 models for solid mechanics problems. We first define a bilinear form, prove some of the important properties, and prove that the W2 formulation will be spatially stable, and convergent to exact solutions. We then present examples of some of the possible W2 models including the SFEM, NS-FEM, ES-FEM, NS-PIM, ES-PIM, and CS-PIM. We show the major properties of these models: (1) they are variationally consistent in a conventional sense, if the solution is sought in a proper H space (compatible cases); (2) They pass the standard patch test when the solution is sought in a proper G space with discontinuous functions (incompatible cases); (3) the stiffness of the discretized model is reduced compared with the finite element method (FEM) model and possibly to the exact model, allowing us to obtain upper bound solutions with respect to both the FEM and the exact solutions and (4) the W2 models are less sensitive to the quality of the mesh, and triangular meshes can be used without any accuracy problems. These properties and theories have been confirmed numerically via examples solved using a number of W2 models including compatible and incompatible cases. We shall see that the G space theory and the W2 forms can formulate a variety of stable and convergent numerical methods with the FEM as one special case. Copyright © 2009 John Wiley & Sons, Ltd. [source]

On the elimination of quadrature subcells for discontinuous functions in the eXtended Finite-Element Method

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 5 2006
G. Ventura
Abstract The introduction of discontinuous/non-differentiable functions in the eXtended Finite-Element Method allows to model discontinuities independent of the mesh structure. However, to compute the stiffness matrix of the elements intersected by the discontinuity, a subdivision of the elements into quadrature subcells aligned with the discontinuity line is commonly adopted. In the paper, it is shown how standard Gauss quadrature can be used in the elements containing the discontinuity without splitting the elements into subcells or introducing any additional approximation. The technique is illustrated and developed in one, two and three dimensions for crack and material discontinuity problems. Copyright © 2005 John Wiley & Sons, Ltd. [source]

Effective Borel measurability and reducibility of functions

MLQ- MATHEMATICAL LOGIC QUARTERLY, Issue 1 2005
Vasco Brattka
Abstract The investigation of computational properties of discontinuous functions is an important concern in computable analysis. One method to deal with this subject is to consider effective variants of Borel measurable functions. We introduce such a notion of Borel computability for single-valued as well as for multi-valued functions by a direct effectivization of the classical definition. On Baire space the finite levels of the resulting hierarchy of functions can be characterized using a notion of reducibility for functions and corresponding complete functions. We use this classification and an effective version of a Selection Theorem of Bhattacharya-Srivastava in order to prove a generalization of the Representation Theorem of Kreitz-Weihrauch for Borel measurable functions on computable metric spaces: such functions are Borel measurable on a certain finite level, if and only if they admit a realizer on Baire space of the same quality. This Representation Theorem enables us to introduce a realizer reducibility for functions on metric spaces and we can extend the completeness result to this reducibility. Besides being very useful by itself, this reducibility leads to a new and effective proof of the Banach-Hausdorff-Lebesgue Theorem which connects Borel measurable functions with the Baire functions. Hence, for certain metric spaces the class of Borel computable functions on a certain level is exactly the class of functions which can be expressed as a limit of a pointwise convergent and computable sequence of functions of the next lower level. (© 2004 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim) [source]

General theory of domain decomposition: Indirect methods

NUMERICAL METHODS FOR PARTIAL DIFFERENTIAL EQUATIONS, Issue 3 2002
Ismael Herrera
Abstract According to a general theory of domain decomposition methods (DDM), recently proposed by Herrera, DDM may be classified into two broad categories: direct and indirect (or Trefftz-Herrera methods). This article is devoted to formulate systematically indirect methods and apply them to differential equations in several dimensions. They have interest since they subsume some of the best-known formulations of domain decomposition methods, such as those based on the application of Steklov-Poincaré operators. Trefftz-Herrera approach is based on a special kind of Green's formulas applicable to discontinuous functions, and one of their essential features is the use of weighting functions which yield information, about the sought solution, at the internal boundary of the domain decomposition exclusively. A special class of Sobolev spaces is introduced in which boundary value problems with prescribed jumps at the internal boundary are formulated. Green's formulas applicable in such Sobolev spaces, which contain discontinuous functions, are established and from them the general framework for indirect methods is derived. Guidelines for the construction of the special kind of test functions are then supplied and, as an illustration, the method is applied to elliptic problems in several dimensions. A nonstandard method of collocation is derived in this manner, which possesses significant advantages over more standard procedures. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 296,322, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10008 [source]

On the discontinuity of the costates for optimal control problems with Coulomb friction

OPTIMAL CONTROL APPLICATIONS AND METHODS, Issue 4 2001
Brian J. Driessen
Abstract This work points out that the costates are actually discontinuous functions of time for optimal control problems with Coulomb friction. In particular these discontinuities occur at the time points where the velocity of the system changes sign. To our knowledge, this has not been noted before. This phenomenon is demonstrated on a minimum-time problem with Coulomb friction and the consistency of discontinuous costates and switching functions with respect to the input switches is shown. Copyright © 2001 John Wiley & Sons, Ltd. [source]