Lumped mass matrices with non−zero inertia for general shell and axisymmetric shell elements

A diagonal lumped mas formulation with non?zero inertia terms is presented for Ahmad's general shell element. The effect of co?ordinate transformation of the mas matrix is demonstrated. Due to arbitray co?ordinate transformation on nodal variables, the diagonal lumped mas matrix becomes a banded matrix of half bandwidth three. It is shown that with some approximations this matrix can be made diagnoal without effecting the results appreciably. Numerical examples are presented to illustrate the accuracy of the formulation. Similar lumped mass matrix formulation is also given for axisymmetric shell elements.     

Efficient formulation for dynamics of corotational 2D beams

The corotational method is an attractive approach to derive non-linear finite beam elements. In a number of papers, this method was employed to investigate the non-linear dynamic analysis of 2D beams. However, most of the approaches found in the literature adopted either a lumped mass matrix or linear local interpolations to derive the inertia terms (which gives the classical linear and constant Timoshenko mass matrix), although local cubic interpolations were used to derive the elastic force vector and the tangent stiffness matrix. In this paper, a new corotational formulation for dynamic nonlinear analysis is presented. Cubic interpolations are used to derive both the inertia and elastic terms. Numerical examples show that the proposed approach is more efficient than using lumped or Timoshenko mass matrices.     

On generation of lumped mass matrices in partition of unity based methods

The partition of unity based methods, such as the extended finite element method and the numerical manifold method, are able to construct global functions that accurately reflect local behaviors through introducing locally defined basis functions beyond polynomials. In the dynamic analysis of cracked bodies using an explicit time integration algorithm, as a result, huge difficulties arise in deriving lumped mass matrices because of the presence of those physically meaningless degrees of freedom associated with those locally defined functions. Observing no spatial derivatives of trial or test functions exist in the virtual work of inertia force, we approximate the virtual work of inertia force in a coarser manner than the virtual work of stresses, where we inversely utilize the ‘from local to global’ skill. The proposed lumped mass matrix is strictly diagonal and can yield the results in agreement with the consistent mass matrix, but has more excellent dynamic property than the latter. Meanwhile, the critical time step of the numerical manifold method equipped with an explicit time integration scheme and the proposed mass lumping scheme does not decrease even if the crack in study approaches the mesh nodes — a very excellent dynamic property. Copyright © 2017 John Wiley & Sons, Ltd.     

一种新型裂尖加强函数的显式动态扩展有限元法

基于扩展有限元方法提出了一种新的裂尖加强函数,与传统三角函数基表征的加强函数相比,该裂尖加强函数通过组合传统的函数基,继承了传统附加函数的特性,同时使得结点的奇异附加自由度减少为2个,减少了总体劲度矩阵的规模,提高了计算效率。通过集中质量矩阵考虑结构的惯性效应,使用显式时间积分方法计算了含裂纹结构的瞬间受载问题,并应用相互作用积分得到裂尖端点处的动态应力强度因子。通过相关算例的对比分析,验证了所提出的裂尖加强函数的有效性,同时表明采用显式时间积分方法进行结构动态响应分析的可行性及准确性。     

Free vibrations of shear‐flexible and compressible arches by FEM

The purpose of this paper is to analyse free vibrations of arches with influence of shear and axial forces taken into account. Arches with various depth of cross‐section and various types of supports are considered. In the calculations, the curved finite element elaborated by the authors is adopted. It is the plane two‐node, six‐degree‐of‐freedom arch element with constant curvature. Its application to the static analysis yields the exact results, coinciding with the analytical ones. This feature results from the use of the exact shape functions in derivation of the element stiffness matrix. In the free vibration analysis the consistent mass matrix is used. It is obtained on the base of the same functions. Their coefficients contain the influences of shear flexibility and compressibility of the arch. The numerical results are compared with the results obtained for the simple diagonal mass matrix representing the lumped mass model. The natural frequencies are also compared with the ones for the continuous arches for which the analytically determined frequencies are known. The advantage of the paper is a thorough analysis of selected examples, where the influences of shear forces, axial forces as well as the rotary and tangential inertia on the natural frequencies are examined. Copyright © 2001 John Wiley & Sons, Ltd.     

附带有考虑集中质量的转动惯性的梁固有振动分析

研究梁附带集中质量时系统的横向弯曲自由振动,分析过程中同时考虑集中质量的平(移)动和转动惯性。探讨了梁的固有频率相对集中质量位置的一阶导数计算问题,并得到了正确的灵敏度计算公式。数值计算结果表明:集中质量的转动惯性对梁的频率、振型以及灵敏度都有很大的影响。当集中质量的转动惯性较大时,忽略其影响对梁的振动分析可能带来很大的误差。     

Use of equivalent mass method for free vibration analyses of a beam carrying multiple two‐dof spring–mass systems with inertia effect of the helical springs considered

This paper investigates the free vibration characteristics of a beam carrying multiple two‐degree‐of‐freedom (two‐dof) spring–mass systems (i.e. the loaded beam). Unlike the existing literature to neglect the inertia effect of the helical springs of each spring–mass system, this paper takes the last inertia effect into consideration. To this end, a technique to replace each two‐dof spring–mass system by a set of rigidly attached equivalent masses is presented, so that the free vibration characteristics of a loaded beam can be predicted from those of the same beam carrying multiple rigidly attached equivalent masses. In which, the equation of motion of the loaded beam is derived analytically by means of the expansion theorem (or the mode superposition method) incorporated with the natural frequencies and the mode shapes of the bare beam (i.e. the beam carrying nothing). In addition, the mass and stiffness matrices including the inertia effect of the helical springs of a two‐dof spring–mass system, required by the conventional finite element method (FEM), are also derived. All the numerical results obtained from the presented equivalent mass method (EMM) are compared with those obtained from FEM and satisfactory agreement is achieved. Because the equivalent masses of each two‐dof spring–mass system are dependent on the magnitudes of its lumped mass, spring constant and spring mass, the presented EMM provides an effective technique for evaluating the overall inertia effect of the two‐dof spring–mass systems attached to the beam. Furthermore, if the total number of two‐dof spring–mass systems attached to the beam is large, then the order of the overall property matrices for the equation of motion of the loaded beam in EMM is much less than that in FEM and the computer storage memory required by the former is also much less than that required by the latter. Copyright © 2005 John Wiley & Sons, Ltd.     

Linear finite elements with orthogonal discontinuous basis functions for explicit earthquake ground motion modeling

The dynamic explicit finite element method is commonly used in earthquake ground motion modeling. In this method, the element mass matrix is approximately lumped, which may lead to numerical dispersion. On the other hand, the orthogonal finite element method, based on orthogonal polynomial basis functions, naturally derives a lumped diagonal mass matrix and can be applied to dynamic explicit finite element analysis. In this paper, we propose finite elements based on orthogonal discontinuous basis functions, the element mass matrices of which are lumped without approximation. Orthogonal discontinuous basis functions are used to improve the accuracy and reduce the numerical dispersion in earthquake ground motion modeling. We present a detailed formulation of the 4‐node tetrahedral and 8‐node hexahedral elements. The relationship between the proposed finite elements and conventional finite elements is investigated, and the solutions obtained from the conventional explicit finite element method are compared with analytical solutions to verify the numerical dispersion caused by the lumping approximation. Comparison of solutions obtained with the proposed finite elements to analytical solutions demonstrates the usefulness of the technique. Examples are also presented to illustrate the effectiveness of the proposed method in earthquake ground motion modeling in the actual three‐dimensional crust structure. Copyright © 2010 John Wiley & Sons, Ltd.     

Computational aspects of a new multi‐scale dispersive gradient elasticity model with micro‐inertia

Computational aspects of a recently developed gradient elasticity model are discussed in this paper. This model includes the (Aifantis) strain gradient term along with two higher‐order acceleration terms (micro‐inertia contributions). It has been demonstrated that the presence of these three gradient terms enables one to capture the dispersive wave propagation with great accuracy. In this paper, the discretisation details of this model are thoroughly investigated, including both discretisation in time and in space. Firstly, the critical time step is derived that is relevant for conditionally stable time integrators. Secondly, recommendations on how to choose the numerical parameters, primarily the element size and time step, are given by comparing the dispersion behaviour of the original higher‐order continuum with that of the discretised medium. In so doing, the accuracy of the discretised model can be assessed a priori depending on the selected discretisation parameters for given length‐scales. A set of guidelines can therefore be established to select optimal discretisation parameters that balance computational efficiency and numerical accuracy. These guidelines are then verified numerically by examining the wave propagation in a one‐dimensional bar as well as in a two‐dimensional example. Copyright © 2016 John Wiley & Sons, Ltd.     

Direct and sparse construction of consistent inverse mass matrices: general variational formulation and application to selective mass scaling

Classical explicit finite element formulations rely on lumped mass matrices. A diagonalized mass matrix enables a trivial computation of the acceleration vector from the force vector. Recently, non‐diagonal mass matrices for explicit finite element analysis (FEA) have received attention due to the selective mass scaling (SMS) technique. SMS allows larger time step sizes without substantial loss of accuracy. However, an expensive solution for accelerations is required at each time step. In the present study, this problem is solved by directly constructing the inverse mass matrix. First, a consistent and sparse inverse mass matrix is built from the modified Hamiltons principle with independent displacement and momentum variables. Usage of biorthogonal bases for momentum allows elimination of momentum unknowns without matrix inversions and directly yields the inverse mass matrix denoted here as reciprocal mass matrix (RMM). Secondly, a variational mass scaling technique is applied to the RMM. It is based on the penalized Hamiltons principle with an additional velocity variable and a free parameter. Using element‐wise bases for velocity and a local elimination yields variationally scaled RMM. Thirdly, examples illustrating the efficiency of the proposed method for simplex elements are presented and discussed. Copyright © 2014 John Wiley & Sons, Ltd.     

A new mass lumping scheme for the mixed hybrid finite element method

Diffusion‐type partial differential equation is a common mathematical model in physics. Solved by mixed finite elements, it leads to a system matrix which is not always an M‐matrix. Therefore, the numerical solution may exhibit unphysical results due to oscillations. The criterion necessary to obtain an M‐matrix is discussed in details for triangular, rectangular and tetrahedral elements. It is shown that the system matrix is never an M‐matrix for rectangular elements and can be an M‐matrix for triangular an tetrahedral elements if criteria on the element's shape and on the time step length are fulfilled. A new mass lumping scheme is developed which leads to a less restrictive criterion: the discretization must be weakly acute (all angles less than π/2) and there is no constraint on the time step length. The lumped formulation of mixed hybrid finite element can be applied not only to triangular meshes but also to more general shape elements in two and three dimensions. Numerical experiments show that, compared to the standard mixed hybrid formulation, the lumping scheme avoids (or strongly reduce) oscillations and does not create additional numerical errors. Copyright © 2006 John Wiley & Sons, Ltd.     

Selective mass scaling for distorted solid‐shell elements in explicit dynamics: optimal scaling factor and stable time step estimate

The use of solid‐shell elements in explicit dynamics has been so far limited by the small critical time step resulting from the small thickness of these elements in comparison with the in‐plane dimensions. To reduce the element highest eigenfrequency in inertia dominated problems, the selective mass scaling approach previously proposed in [G. Cocchetti, M. Pagani and U. Perego, Comp. & Struct. 2013; 127:39‐52.] for parallelepiped elements is here reformulated for distorted solid‐shell elements. The two following objectives are achieved: the critical time step is governed by the smallest element in‐plane dimension and not anymore by the thickness; the mass matrix remains diagonal after the selective mass scaling. The proposed approach makes reference to one Gauss point, trilinear brick element, for which the maximum eigenfrequency can be computed analytically. For this element, it is shown that the proposed mass scaling can be interpreted as a geometric thickness scaling, obtaining in this way a simple criterion for the definition of the optimal mass scaling factor. A strategy for the effective computation of the element maximum eigenfrequency is also proposed. The considered mass scaling preserves the element translational inertia, while it modifies the rotational one, leading to errors in the kinetic energy when the motion rotational component is dominant. The error has been rigorously assessed for an individual element, and a simple formula for its estimate has been derived. Numerical tests, both in small and large displacements and rotations, using a state‐of‐the‐art solid‐shell element taken from the literature, confirm the effectiveness and accuracy of the proposed approach. Copyright © 2014 John Wiley & Sons, Ltd.     

《剩余质量阵和剩余刚度阵的近似计算方法》

在模态综合法中，剩余质量阵和剩余刚度阵的计算通常依赖于子结构的质量阵和刚度阵，由于几乎不可能通过实验方法得到子结构的质量阵和刚度阵，因此很难将实验数据应用于理论模态综合法中。针对这一困难，通过分析剩余质量阵和刚度阵的表达式，推导出它们的近似计算公式，这一近似计算公式中不含质量阵和刚度阵，亦即不需要事先知道质量阵和刚度阵就可以计算剩余质量阵和剩余刚度阵，从而回避实验辨识质量阵和刚度阵这一极具挑战的动力学反问题，克服实验模态综合法应用中的困难。数值计算表明，所提的剩余质量阵和剩余刚度阵的近似计算方法切实可行，且模态综合结果具有较高的精度。     

Efficient explicit time stepping for the eXtended Finite Element Method (X‐FEM)

This paper focuses on the introduction of a lumped mass matrix for enriched elements, which enables one to use a pure explicit formulation in X‐FEM applications. A proof of stability for the 1D and 2D cases is given. We show that if one uses this technique, the critical time step does not tend to zero as the support of the discontinuity reaches the boundaries of the elements. We also show that the X‐FEM element's critical time step is of the same order as that of the corresponding element without extended degrees of freedom. Copyright © 2006 John Wiley & Sons, Ltd.     

Conservative integration of rigid body motion by quaternion parameters with implicit constraints

An angular momentum and energy‐conserving time integration algorithm for rigid body rotation is formulated in terms of the quaternion parameters and the corresponding four‐component conjugate momentum vector via Hamilton's equations. The introduction of an extended mass matrix leads to a symmetric set of eight state‐space equations of motion. The extra inertial parameter serves as a multiplier on the kinematic constraint, and it is demonstrated that convergence characteristics are improved by selecting this parameter somewhat larger than the inertial moments. External loads enter these equations via the set of momentum equations. Initially, the normalization of the quaternion array is introduced via a Lagrange multiplier. However, this Lagrange multiplier can be expressed explicitly in terms of the gradient of the external load potential, and elimination of the Lagrange multiplier from the final format leaves only an explicit projection applied to the external load potential gradient. An algorithm is developed by forming a finite increment of the Hamiltonian. This procedure identifies the proper selection of increments and mean values, and leads to an algorithm with conservation of momentum and energy. Implementation, conservation properties, and accuracy of the algorithm are illustrated by examples with a flying box and a spinning top. Copyright © 2012 John Wiley & Sons, Ltd.     

Non‐local dispersive model for wave propagation in heterogeneous media: multi‐dimensional case

Three non‐dispersive models in multi‐dimensions have been developed. The first model consists of a leading‐order homogenized equation of motion subjected to the secularity constraints imposing uniform validity of asymptotic expansions. The second, non‐local model, contains a fourth‐order spatial derivative and thus requires C¹ continuous finite element formulation. The third model, which is limited to the constant mass density and a macroscopically orthotropic heterogeneous medium, requires C⁰ continuity only and its finite element formulation is almost identical to the classical local approach with the exception of the mass matrix. The modified mass matrix consists of the classical mass matrix (lumped or consistent) perturbed with a stiffness matrix whose constitutive matrix depends on the unit cell solution. Numerical results are presented to validate the present formulations. Copyright © 2002 John Wiley & Sons, Ltd.     

An efficient coarse‐grained parallel algorithm for global–local multiscale computations on massively parallel systems

The existing global–local multiscale computational methods, using finite element discretization at both the macro‐scale and micro‐scale, are intensive both in terms of computational time and memory requirements and their parallelization using domain decomposition methods incur substantial communication overhead, limiting their application. We are interested in a class of explicit global–local multiscale methods whose architecture significantly reduces this communication overhead on massively parallel machines. However, a naïve task decomposition based on distributing individual macro‐scale integration points to a single group of processors is not optimal and leads to communication overheads and idling of processors. To overcome this problem, we have developed a novel coarse‐grained parallel algorithm in which groups of macro‐scale integration points are distributed to a layer of processors. Each processor in this layer communicates locally with a group of processors that are responsible for the micro‐scale computations. The overlapping groups of processors are shown to achieve optimal concurrency at significantly reduced communication overhead. Several example problems are presented to demonstrate the efficiency of the proposed algorithm. Copyright © 2009 John Wiley & Sons, Ltd.     

ITERATIVE METHOD USING CONSISTENT MASS MATRIX IN AXISYMMETRICAL FINITE ELEMENT ANALYSIS OF HYPERVELOCITY IMPACT

We present in this paper an iterative method using consistent mass matrix in axisymmetrical finite element analysis of hypervelocity impact. To retain the advantage of integration on an element-by-element basis which is at the heart of modern hydrocodes, we suggest that the first step should be to solve for accelerations at an advanced time step by using the lumped mass approach, then iterate using a consistent mass matrix to improve the estimate. Examples are given to show the improved resolution with the new method.     

Free vibration analysis of arches using curved beam elements

The natural frequencies and mode shapes for the radial (in‐plane) bending vibrations of the uniform circular arches were investigated by means of the finite arch (curved beam) elements. Instead of the complicated explicit shape functions of the arch element given by the existing literature, the simple implicit shape functions associated with the tangential, radial (or normal) and rotational displacements of the arch element were derived and presented in matrix form. Based on the relationship between the nodal forces and the nodal displacements of a two‐node six‐degree‐of‐freedom arch element, the elemental stiffness matrix was derived, and based on the equation of kinetic energy and the implicit shape functions of an arch element the elemental consistent mass matrix with rotary inertia effect considered was obtained. Assembly of the foregoing elemental property matrices yields the overall stiffness and mass matrices of the complete curved beam. The standard techniques were used to determine the natural frequencies and mode shapes for the curved beam with various boundary conditions and subtended angles. In addition to the typical circular arches with constant curvatures, a hybrid beam constructed by using an arch segment connected with a straight beam segment at each of its two ends was also studied. For simplicity, a lumped mass model for the arch element was also presented. All numerical results were compared with the existing literature or those obtained from the finite element method based on the conventional straight beam element and good agreements were achieved. Copyright © 2003 John Wiley & Sons, Ltd.     

Explicit finite element artificial boundary scheme for transient scalar waves in two‐dimensional unbounded waveguide

To simulate the transient scalar wave propagation in a two‐dimensional unbounded waveguide, an explicit finite element artificial boundary scheme is proposed, which couples the standard dynamic finite element method for complex near field and a high‐order accurate artificial boundary condition (ABC) for simple far field. An exact dynamic‐stiffness ABC that is global in space and time is constructed. A temporal localization method is developed, which consists of the rational function approximation in the frequency domain and the auxiliary variable realization into time domain. This method is applied to the dynamic‐stiffness ABC to result in a high‐order accurate ABC that is local in time but global in space. By discretizing the high‐order accurate ABC along artificial boundary and coupling the result with the standard lumped‐mass finite element equation of near field, a coupled dynamic equation is obtained, which is a symmetric system of purely second‐order ordinary differential equations in time with the diagonal mass and non‐diagonal damping matrices. A new explicit time integration algorithm in structural dynamics is used to solve this equation. Numerical examples are given to demonstrate the effectiveness of the proposed scheme. Copyright © 2011 John Wiley & Sons, Ltd.