Search results

The Convolution Quadrature Method (CQM) has been widely applied to solve transient elastodynamic problems because of its stability and generality. However, the CQM suffers from the problems of huge memory requirement in case of direct implementation in time domain or CPU time in case of its reformulation in Laplace domain. The purpose of this paper is to combine the CQM with the pseudo-initial condition method (PICM) to achieve a good balance between memory requirement and CPU time.

The combined methods first subdivide the whole analysis into a few sub-analyses, which is dealt with the PICM, namely, the results obtained by previous sub-analysis are used as the initial conditions for the next sub-analysis. In each sub-analysis, the time interval is further discretized into a number of sub-steps and dealt with the CQM. For non-zero initial conditions, the pseudo-force method is used to transform them into equivalent body forces. The boundary face method is employed in the numerical implementation. Three examples are analyzed. Results are compared with analytical solutions or FEM results and the results of reformulated CQM.

Results demonstrate that the computation time and the storage requirement can be reduced significantly as compared to the CQM, by using the combined approach.

The combined methods can be successfully applied to the problems of long-time dynamic response, which requires a large amount of computer memory when CQM is applied, while preserving the CQM stability. If the number of time steps is high, then the accuracy of the proposed approach can be deteriorated because of the pseudo-force method.

This paper aims to study the seismic response of frame structure with friction dampers.

The state equation of the structure subjected to the earthquake is presented and solved, from which the maximum drift and the interlayer drift angle of the floors of the structure subjected to the seismic waves of four types of sites are analyzed.

The result indicates that the damping effect is significant on the floors with the friction damper but is almost little influence on the other floor.

The result indicates that the damping effect is significant on the floors with the friction damper but is almost little influence on the other floor.

The purpose of this paper is to investigate internal forces in bridges induced by moving vehicles and compare them to earthquake loading.

Dynamic analysis of bridges is performed for moving support actions, for spectral method with Eurocode 8 parameters and for moving vehicle influence. Results from all three methods have been compared on two examples and conclusions have been made. Moving vehicle analysis could be based on the moving force and on the moving mass approach where the later one requires rather accurate knowledge of structural accelerations. It has been shown that the classical Newmark formulation produces accelerations of low accuracy and a novel impulse acceleration method has been devised.

It is found that the actions induced by the moving load could be comparable or larger than those caused by the earthquake on bridges whose mass is not too large in comparison to the vehicle mass.

The developed method will be applied to a broader choice of examples and more reliable conclusions made.

There are bridges where it would be appropriate to perform moving vehicle dynamic analysis, in which case the vertical earthquake actions could be neglected in the analysis.

In order to assess actions from moving vehicles, Newmark method has been generalized in a novel way. Paper describes vector formulation of Newmark method that permits free mixing of integration parameters that could vary from node to node. The method is advantageous for moving load analysis where loading conditions of nodes change in time.

In a finite element model of typical dam—foundation—reservoir systems, the presence of heterogeneous material properties for the dam and the foundation produces a combined damping matrix that is non‐proportional to the mass and/or the stiffness matrices of the system. In this case, the undamped real free‐vibration modes cannot uncouple the damping forces such that the classical mode superposition method using real modes is not applicable. This paper presents comparative analyses of recent coordinate reduction procedures that have been developed to compute the response of linear systems with non‐proportional damping. The comparisons are based on the numerical efficiency and the accuracy of the displacement, acceleration and stress response, and on the distribution of the damping energy in the system.

The purpose of this paper is to describe further developments on a novel formulation of the boundary element method (BEM) for inelastic problems using the dual reciprocity method (DRM) but using object-oriented programming (OOP). As the BEM formulation generates a domain integral due to the inelastic stresses, the DRM is employed in a modified form using polyharmonic spline approximating functions with polynomial augmentation. These approximating functions produced accurate results in BEM applications for a range of problems tested, and have been shown to converge linearly as the order of the function increases.

A programming class named DRMOOP, written in C++ language and based on OOP, was developed in this research. With such programming, general matrix equations can be easily established and applied to different inelastic problems. A vector that accounts for the influence of the inelastic strains on the displacements and boundary forces is obtained.

The C++ DRMOOP class has been implemented and tested with the BEM formulation applied to classical elastoplastic problem and the results are reported at the end of the paper.

An object-oriented technology and the C++ DRMOOP class applied to elastoplastic problems.

The free roving robot in the form of a computer controlled industrial truck is now technically possible. The use of small on‐board computers and navigation aids frees the robot from dependence on fixed route marking. Conventional vehicles such as fork lift trucks may be modified for autonomous or semi‐autonomous control.

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.

The purpose of this paper is to simulate flow inside differentially heated rotating cavity using two different formulations; one using Navier‐Stokes (NS) equations derived in non‐inertial (rotating) frame of reference and the other using NS equations in inertial frame of reference. Then to compare the results obtained from these formulations to find their merits and demerits.

The NS equations for both non‐inertial and inertial formulations are written in artificial compressibility form before discretizing them by a high resolution finite volume method. The dual time steeping approach of Jameson is used for time accuracy in both the formulations. Arbitrary Lagrangian Eulerian (ALE) approach is used for taking care of moving boundary problem arising in the inertial formulation. A newly developed HLLC‐AC Riemann solver for discretizing convective fluxes and central differencing for discretizing viscous fluxes are used in the finite volume approach. Results for both the formulations are first validated with test cases reported in literature. Then the results of the two formulations are compared among themselves.

Results of the non‐inertial formulation obtained by the proposed method are found to match well with those reported in literature. The results of both the formulations match well for low rotational speeds of the cavity. The discrepancies between the results of the two formulations progressively increase with the increase in rotational speed. Implicit treatment of the source term is found to reduce the discrepancies.

The present approach is useful for accurate prediction of flow feature and heat transfer characteristic in case of applications such as manufacturing of single wafer crystal for semiconductor and in numerous metallurgical processes.

The ALE formulation is used for the first time to simulate a differentially heated rotating cavity problem. The attempt to compare non‐inertial and inertial formulations is also reported for the first time. Implicit treatment of the source term leading to change in solution accuracy is one of the important findings of the present investigation.

A computational procedure is presented for the efficient non‐linear dynamic analysis of quasi‐symmetric structures. The procedure is based on approximating the unsymmetric response vectors, at each time step, by a linear combination of symmetric and antisymmetric vectors, each obtained using approximately half the degrees of freedom of the original model. A mixed formulation is used with the fundamental unknowns consisting of the internal forces (stress resultants), generalized displacements and velocity components. The spatial discretization is done by using the finite element method, and the governing semi‐discrete finite element equations are cast in the form of first‐order non‐linear ordinary differential equations. The temporal integration is performed by using implicit multistep integration operators. The resulting non‐linear algebraic equations, at each time step, are solved by using iterative techniques. The three key elements of the proposed procedure are: (a) use of mixed finite element models with independent shape functions for the stress resultants, generalized displacements, and velocity components and with the stress resultants allowed to be discontinuous at interelement boundaries; (b) operator splitting, or restructuring of the governing discrete equations of the structure to delineate the contributions to the symmetric and antisymmetric vectors constituting the response; and (c) use of a two‐level iterative process (with nested iteration loops) to generate the symmetric and antisymmetric components of the response vectors at each time step. The top‐ and bottom‐level iterations (outer and inner iterative loops) are performed by using the Newton—Raphson and the preconditioned conjugate gradient (PCG) techniques, respectively. The effectiveness of the proposed strategy is demonstrated by means of a numerical example and the potential of the strategy for solving more complex non‐linear problems is discussed.

A direct integration method for the dynamic analysis of structures equipped with viscoelastic dampers (VED) is presented in this paper. The constitutive model of the damper is set using a system of Maxwell elements composed in parallel (MPS); the evolutive behavior of the VE material is accounted for by introducing temperature‐dependent mechanical properties. The solution procedure follows an incremental approach of implicit type that, by means of a discrete‐time formulation of the problem, allows expression of the damping force discretization in a form suitable to be included in a standard integration scheme. The resulting algorithm is proved to be very effective and robust; the distinctive features of the proposed numerical method suggest the possibility of a successful implementation of an MPS procedure in the frame of a standard finite element code.