In this paper, we develop a general way to construct contact algorithms for contact dynamical systems. Such an algorithm requires the corresponding step transition map preserve the contact structure of the underlying ...
详细信息
In this paper, we develop a general way to construct contact algorithms for contact dynamical systems. Such an algorithm requires the corresponding step transition map preserve the contact structure of the underlying contact phase space. The constructions are based on the correspondence between the contact geometry of R2n+1 and the conic symplectic one of R2n+2 and therefore, the algorithms are derived naturally from the symplectic algorithms of Hamiltonian systems.
作者:
LOCHEGNIES, DOUDIN, JLaboratoire de Génie Mécanique
Unité de recherches Associée au CNRS MECAMAT Université de Valenciennes et du Hainaut—Cambrésis BP 311 59304 Valenciennes Cedex France Laboratoire de Génie Mécanique
Unité de recherches Associée au CNRS MECAMAT Université de Valenciennes et du Hainaut—Cambrésis BP 311 59304 Valenciennes Cedex France
New contact boundary modelling is achieved with a basic set of 2 and 3 dimension contact primitives. contact constraints are originally introduced in the variational equations and associated Newton-Raphson scheme via ...
详细信息
New contact boundary modelling is achieved with a basic set of 2 and 3 dimension contact primitives. contact constraints are originally introduced in the variational equations and associated Newton-Raphson scheme via an external penalty formulation using primitive equations. Consequently, penalty part of external load vector and tangent stiffness matrices are developed for all contact primitives. In this way, contact prescribed boundary displacements are also taken into account. contact treatment is then completed with Newton-Raphson elements for elastic and plastic regularized friction constitutive models. In this paper, the process is extended to elastoplastic models. Finally, we propose a self acting procedure with contact algorithms (interiority, sliding and contact loss) and related subroutines for implementation in finite element framework. We illustrate these developments by means of two-dimensional open die forging and three-dimensional plate coining typical benchmarks with reference to bulk elastoplastic and viscoplastic constitutive models.
This article identifies some issues, options and results associated with 3D Explosively Formed Penetrator (EFP) computations. The three broad categories of interest (for Lagrangian computations) are the sliding/contac...
详细信息
This article identifies some issues, options and results associated with 3D Explosively Formed Penetrator (EFP) computations. The three broad categories of interest (for Lagrangian computations) are the sliding/contact interfaces, the type and arrangement of elements, and the generation of the grid. The sliding/contact algorithm uses a virtual particle approach that is very robust, and it is described in detail. Tetrahedral elements are considered for both symmetric and non-symmetric arrangement, each with some advantages and disadvantages. A mixed element algorithm is also considered for the non-symmetric arrangement, in an effort to eliminate some of the locking that occurs for the non-symmetric arrangement. Hexahedral (brick) elements are not considered herein, but some of the effects (for tetrahedral elements) also apply to hexahedral elements. For axisymmetric EFPs there are also some trade-offs (symmetric grid vs. equally sized elements) associated with the generation of the grid. Computational results are provided to illustrate each of the issues and effects. The virtual particle algorithm is also well suited for high-velocity impact computations, and an example is included to demonstrate this capability for a complex problem. (c) 2005 Elsevier Ltd. All rights reserved.
The reliability and efficiency of three different numerical modeling approaches for simulating the response of a newly developed 3D fiber-metal laminate (3D-FML), subject to axial impact loading, are considered in thi...
详细信息
The reliability and efficiency of three different numerical modeling approaches for simulating the response of a newly developed 3D fiber-metal laminate (3D-FML), subject to axial impact loading, are considered in this paper. The main objective of the study is to establish the most robust numerical framework for analyzing the performance of such complexly configured hybrid materials subject to axial impact loading in a fairly accurate, yet efficient manner. LS-DYNA finite element software is used for the purpose. The models include: (i) a full 3D solid model, where all 3D-FML constituents are modeled with 3D elements;(ii) a model with intermediate complexity, in which two different element types are used to model the metallic skins and 3D-fiberglass/foam core, respectively;and (iii) a simplified scheme, consisting of a single layer of thin-shell elements, representing all constituents of the FML. An experimental investigation is also conducted in parallel to verify the accuracy of the modeling schemes. Force and axial-shortening histories, energy absorption capacity, and overall qualitative behavior obtained numerically are compared to experimental results. Both accuracy and computation cost are considered as the performance criteria, all with the aim of providing the reader with some perspective for robust modeling of such geometrically sophisticated composites, subject to a complex loading mechanism.
This work presents a new technique for joining dissimilar meshes in non-linear finite element analysis, with the novel result that the adjoining surfaces do not need to be spatially coincident for patch tests to be pa...
详细信息
This work presents a new technique for joining dissimilar meshes in non-linear finite element analysis, with the novel result that the adjoining surfaces do not need to be spatially coincident for patch tests to be passed. The work builds in some sense on results previously presented by Dohrmann et al. (International Journal for Numerical Methods in Engineering 1991;40:1205), but considerably generalizes upon the results presented there. Most importantly, it is shown that consistency of the tied finite element formulation may be achieved such that implementation of element formulations is not affected, through the use of a three-field variational principle applied over the contacting region. A key conclusion of the work is that although the approach advocated bears some resemblance to mortar element schemes for contact analysis (see, for example, McDevitt and Laursen, International Journal for Numerical Methods in Engineering 2000;48:1525), a straightforward application of a mortaring scheme is not sufficient to guarantee patch test passage when dissimilar meshing of curved surfaces leads to local open volumes (or overlaps) between tied surfaces in their reference configurations. The proposed algorithm is demonstrated through a number of numerical examples performed within a matrix-free non-linear equation solving framework. Published in 2003 by John Wiley Sons, Ltd.
暂无评论