Chien Wei-zang. Dynamic Finite Element with Diagonalized Consistent Mass Matrix and Elastic-Plastic Impact Calculation[J]. Applied Mathematics and Mechanics, 1982, 3(3): 281-296.
Citation: Chien Wei-zang. Dynamic Finite Element with Diagonalized Consistent Mass Matrix and Elastic-Plastic Impact Calculation[J]. Applied Mathematics and Mechanics, 1982, 3(3): 281-296.

Dynamic Finite Element with Diagonalized Consistent Mass Matrix and Elastic-Plastic Impact Calculation

  • Received Date: 1981-12-05
  • Publish Date: 1982-06-15
  • There are some common difficulties encountered in elastic-plastic impact codes such as EPIC[1],[2] NONSAP[3] etc. Most of these codes use the simple linear functions usually taken from static problems to represent the displacement components. In such finite element formulation, the strain and stress components are constants in every element. In the equations of motion, the stress components in general appear in the form of their space derivatives. Thus, if we use such form functions to represent the displacement components, the effect of internal stresses to the equations of motion vanishes identically. The usual practice to overcome such difficulties is to establish as self-equilibrium system of internal forces acting on various nodal points by means of transforming equations of motion into variational form of energy relation through the application of virtual displacement principle. The nodal acceleration is then calculated from the total force acting on this node from all the neighbouring elements. The transformation of virtual displacement principle into the variational energy form is performed on the bases of continuity conditions of stress and displacement throughout the integrated space. That is to say, on the interface boundary of finite element, the assumed displacement and stress functions should be conformed. However, it is easily seen that, for linear form function of finite element calculation, the displacement continues everywhere, but not the stress components. Thus, the convergence of such kind of finite element computation is open to question. This kind of treatment has never been justified even in approximation sense. Furthermore, the calculation of nodal points needs a rule to calculate the mass matrix. There are two ways to establish mass matrix, namely lumped mass method and consistent mass method[4]. The consistent mass matrix can be obtained naturally through finite element formulation, which is consistent to the assumed form functions. However, the resulting consistent mass matrix is not in dia-gonalized form, which is inconvenient for numerical computation. For most codes, the lumped mass matrix is used, and in this case, the element mass is distributed in certain assumed proportions to all the nodal points of this element. The lumped mass matrix is diagonalized with diagonal terms composed of the nodal mass. However, the lumped mass assumption has never been justified. All these difficulties are originated from the simple linear form functions usually used in static problems. In this paper, we introduce a new quadratic form function for elastic-plastic impact problems. This quadratic form function possesses diagonalized consistent masf matrix, and non-vanishing effect of internal stress to the equations of motion. Thus with this kind of dynamic finite element, all above-said difficulties can be eliminated.
  • loading
  • [1]
    Johnson,G.R.,Analysis of elastic-plastic impact involving severe distorsions.Journal of Applied Mechanics.Vol.43,No.3,Trans.A.S.M.E.,Vol.98,Series E,Sept.1976. pp.439-444.
    [2]
    Johoson,G.R.,High velocity impact calculations in three dimensions,Journal of Applied Mechanics.Vol.44,No.1,Trans.A.S.M.E.,Vol.99,Series E.March.(1977),pp.95-100.
    [3]
    Dathe.K.J.and Wilson.E.L.,NONSAP-A general finite element program for nonlinear dynamic analysis of complex structure,Paper M3/1,Proc.Second Int.Conf.Struct.Mech.Reactor Technology,Berlin(1973).
    [4]
    Zienkiewicz,O.C.,The Finite Element Method (Third Edition.1977), pp,533-540,McGraw-Hill, London.
    [5]
    Boresi,A.P.,Elasticity in engineering mechanics.Prentice-Hall Inc.,Englewood Cliffs,N.J.(1965).
    [6]
    Von Neumann.J.and Richtmyer,R.D.,A method for the numerical calculation of hydrodynamic shocks,Journal of Applied Physics.Vol.21,1950,pp.232-237.
    [7]
    Walsh.J.M.,et al.,Shock wave compressions of twenty seven metals.Equations of state of metals.Physical Review.Vol.108. No.2,Oct.1957,pp.196-216.
  • 加载中

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Article Metrics

    Article views (1616) PDF downloads(726) Cited by()
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return