Search results

The time‐discretization process of transient equation systems is an important concern in computational heat transfer applications. As such, the present paper describes a formal basis towards providing the theoretical concepts, evolution and development, and characterization of a wide class of time discretized operators for transient heat transfer computations. Therein, emanating from a common family tree and explained via a generalized time weighted philosophy, the paper addresses the development and evolution of time integral operators [IO], and leading to integration operators [InO] in time encompassing single‐step integration operators [SSInO], multi‐step integration operators [MSInO], and a class of finite element in time integration operators [FETInO] including the relationships and the resulting consequences. Also depicted are those termed as discrete numerically assigned [DNA] algorithmic markers essentially comprising of both: the weighted time fields, and the corresponding conditions imposed upon the dependent variable approximation, to uniquely characterize a wide class of transient algorithms. Thereby, providing a plausible standardized formal ideology when referring to and/or relating time discretized operators applicable to transient heat transfer computations.

In structural, earthquake and aeronautical engineering and mechanical vibration, the solution of dynamic equations for a structure subjected to dynamic loading leads to a high order system of differential equations. The numerical methods are usually used for integration when either there is dealing with discrete data or there is no analytical solution for the equations. Since the numerical methods with more accuracy and stability give more accurate results in structural responses, there is a need to improve the existing methods or develop new ones. The paper aims to discuss these issues.

In this paper, a new time integration method is proposed mathematically and numerically, which is accordingly applied to single-degree-of-freedom (SDOF) and multi-degree-of-freedom (MDOF) systems. Finally, the results are compared to the existing methods such as Newmark's method and closed form solution.

It is concluded that, in the proposed method, the data variance of each set of structural responses such as displacement, velocity, or acceleration in different time steps is less than those in Newmark's method, and the proposed method is more accurate and stable than Newmark's method and is capable of analyzing the structure at fewer numbers of iteration or computation cycles, hence less time-consuming.

A new mathematical and numerical time integration method is proposed for the computation of structural responses with higher accuracy and stability, lower data variance, and fewer numbers of iterations for computational cycles.

Various time integration methods and time finite element methods have been developed to obtain the responses of structural dynamic problems, but the accuracy and computational efficiency of them are sometimes not satisfactory. The purpose of this paper is to present a more accurate and efficient formulation on the basis of the weak form quadrature element method to solve linear structural dynamic problems.

A variational principle for linear structural dynamics, which is inspired by Noble's work, is proposed to develop the weak form temporal quadrature element formulation. With Lobatto quadrature rule and the differential quadrature analog, a system of linear equations is obtained to solve the responses at sampling time points simultaneously. Computation for multi-elements can be carried out by a time-marching technique, using the end point results of the last element as the initial conditions for the next.

The weak form temporal quadrature element formulation is conditionally stable. The relation between the normalized length of element and the suggested number of integration points in one element is given by a simple formula. Results show that the present formulation is much more accurate than other time integration methods and its dissipative property is also illustrated.

The weak form temporal quadrature element formulation provides a choice with high accuracy and efficiency for solution of linear structural dynamic problems.

This paper discusses algorithms for large displacement analysis of interconnected flexible and rigid multibody systems. Hydrostatic and hydrodynamic loads for systems being submerged in water are also considered. The systems may consist of cables and beams and may combine very flexible parts with rigid parts. Various ways of introducing structural joints are discussed. A special implementation of the Hilber‐Hughes‐Taylor time integration scheme for constrained non‐linear systems is outlined. The formulation is general and allows for displacements and rotational motion of unlimited size. Aspects concerning efficient solution of constrained dynamic problems are discussed. These capabilities have been implemented in a general purpose non‐linear finite element program. Applications involving static and dynamic analysis of a bi‐articulated tower and a floating tripod platform kept in place by three anchor lines are discussed.

Partitioned analysis is a method by which sets of time‐dependent ordinary differential equations for coupled systems may be numerically integrated in tandem, thereby avoiding brute‐force simultaneous solution. The coupled systems addressed pertain to fluid—structure, fluid—soil, soil—structure, or even structure—structure interaction. The paper describes the partitioning process for certain discrete‐element equations of motion, as well as the associated computer implementation. It then delineates the procedure for designing a partitioned analysis method in a given application. Finally, examples are presented to illustrate the concepts. It is seen that a key element in the implementation of partitioned analysis is the use of integrated, as opposed to monolithic software.

A generalized time finite element method is considered for time integration of non‐linear equations of motion arising from dynamic problems. A simple three‐time level family of schemes is obtained. Evaluation of the schemes shows that the proposed approach may lead to unconditionally stable algorithms for non‐linear problems. Numerical examples show the accuracy and efficiency of the proposed algorithm in comparison to Newmark's average acceleration method and four order accurate explicit algorithm.

The purpose of this paper is to propose a new explicit time integration algorithm for solution to the linear and non-linear finite element equations of structural dynamic and wave propagation problems.

The algorithm is completely explicit so that no linear equation system requires solving, if the mass matrix of the finite element equation is diagonal and whether the damping matrix does or not. The algorithm is a single-step method that has the simple starting and is applicable to the analysis with the variable time step size. The algorithm is second-order accurate and conditionally stable. Its numerical stability, dissipation and dispersion are analyzed for the dynamic single-degree-of-freedom equation. The stability of the multi-degrees-of-freedom non-proportional damping system can be evaluated directly by the stability theory on ordinary differential equation.

The performance of the proposed algorithm is demonstrated by several numerical examples including the linear single-degree-of-freedom problem, non-linear two-degree-of-freedom problem, wave propagation problem in two-dimensional layer and seismic elastoplastic analysis of high-rise structure.

A new single-step second-order accurate explicit time integration algorithm is proposed to solve the linear and non-linear dynamic finite element equations. The algorithm has advantages on the numerical stability and accuracy over the existing modified central difference method and Chung-Lee method though the theory and numerical analyses.

This paper aims to present an adaptive approach of the generalized finite element method (GFEM) for transient dynamic analysis of bars and trusses. The adaptive GFEM, previously proposed for free vibration analysis, is used with the modal superposition method to obtain precise time-history responses.

The adaptive GFEM is applied to the transient analysis of bars and trusses. To increase the precision of the results and computational efficiency, the modal matrix is responsible for the decoupling of the dynamic equilibrium equations in the modal superposition method, which is used with only the presence of the problem’s most preponderant modes of vibration. These modes of vibration are identified by a proposed coefficient capable of indicating the influence of each mode on the transient response.

The proposed approach leads to more accurate results of displacement, velocity and acceleration when compared to the traditional finite element method.

In this paper, the application of the adaptive GFEM to the transient analysis of bars and trusses is presented for the first time. A methodology of identification of the preponderant modes to be retained in the modal matrix is proposed to improve the quality of the solution. The examples showed that the method has a strong potential to solve dynamic analysis problems, as the approach had already proved to be efficient in the modal analysis of different framed structures. A simple way to perform h-refinement of truss elements to obtain reference solutions for dynamic problems is also proposed.

Addresses the computational aspects involved in the numerical simulation of sheet stamping processes. Focuses on some numerical aspects of the intrinsic complexity of these problems, the first of which is the necessity to take into account properly membrane and bending effects. Presents a well‐adapted shell element. The second aspect concerns the description and the implementation of the initial orthotropic plastic behaviour for sheet metal parts, based on a formulation in a rotating frame using the initial microstructure rotation. The stress calculation algorithm is based on a particular implementation of the elastic predictor‐plastic corrector method. The last aspect concerns the solution procedures with a particular development concerning the treatment of the blankholder load as a constraint. A set of computational results validated with experiments prove the accuracy of the proposed approach in solving stamping problems.

nonlinear dynamic analysis of triangular and quadrilateral membrane elements with in-plane drilling rotational degree of freedom.

The nonlinear analysis is carried out using the updated co-rotational Lagrangian description. In this purpose, in-plane co-rotational formulation that considers the in-plane drilling rotation is developed and presented for triangular and quadrilateral elements, and a tangent stiffness matrix is derived. Furthermore, a simple and effective in-plane mass matrix that takes into account the in-plane rotational inertia, which permit true representation of in-plane vibrational modes is adopted for dynamic analysis, which is carried out using the Newmark direct time integration method.

The proposed numerical tests show that the presented elements exhibit very good performances and could return true in-plane rotational vibrational modes. Also, when using a well-chosen co-rotational formulation these elements shows good results for nonlinear static and dynamic analysis.

Publications that describe geometrical nonlinearity of the in-plane behaviour of membrane element with rotational d.o.f are few, and often they are based on the total Lagrangian formulation or on the rate form. Also these elements, at the author knowledge, have not been extended to the nonlinear dynamic analysis. Thus, an appropriate extension of triangular and quadrilateral membrane elements with drilling rotation to nonlinear dynamic analysis is required.