Search results

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.

Uses Newmark’s method to carry out a time‐stepping finite element analysis to predict the behaviour of a cloth garment as it falls from an initial horizontal position to a final position draped around a human body form. Bases the finite element model on a simple beam element, in order to minimize the computational time. Accounts for large displacement behaviour by including the element geometric stiffness. Bases the body form on anthropomorphic data produced by a shadow scanner. Enlists a novel scheme to model the contact between the cloth and the underlying body form. Uses the finite element model to provide data for an animated display and finds that it produces sufficiently realistic results for the garment designer’s purposes.

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.

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.

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.

Considers hysteretic non‐linear models for modelling horizontal force‐displacement characteristics of an isolation system subjected to general plane motion. Simple close form solution of the stiff differential equation of hysteretic model for forces mobilized in the non‐linear elements of the base isolation system are obtained. Simulates experimental shear force‐displacement loops obtained from the bi‐axial tests by different investigators. Both the experimental and the simulated hysteresis loops are found to be in good agreement. Develops a unified solution algorithm for computation of response of different types of base isolated buildings, considering non‐linear behaviour of the isolation systems, subjected to multi‐directional motion. The solution algorithm is based on the implicit‐implicit partitioned Newmark’s method in predictor‐corrector form. Response of the base isolated symmetrical building as obtained from the solution algorithm and the computer programs developed are in good agreement with that obtained from the more complex numerical studies reported in the literature. Response of a three storeyed symmetrical building isolated by pure friction isolators and laminated rubber bearings has been obtained using this simple yet accurate solution algorithm.

This paper demonstrates the capability of staggered solution procedure for coupled fluid‐structure interaction problems. Three possible computational paths for coupled problems are described. These are critically examined for a variety of coupled problems with different types of mesh partitioning schemes. The results are compared with the reported results by continuum mechanics priority approach—a method which has been very popular until recently. Optimum computational paths and mesh partitionings for two field problems are indicated. Staggered solution procedure is shown to be quite effective when optimum path and partitionings are selected.

The present work involves analytical and experimental investigation of sloshing in a two-dimensional rectangular tank including the effect of porous baffles to control and/or reduce the wave motion in the sloshing tank. The purpose of this study is to assess the analytical solutions of the drag coefficient effect on porous baffles performance to track free surface motion variation in the sloshing tank by comparison with experimental shake table tests under a range of sway excitation.

The linear second-order ordinary differential equations for liquid sloshing in the rectangular tank were solved using Newmark’s beta method and obtained the analytical solutions for liquid sloshing with dual vertical porous baffles of full submergence depths in a sway-oscillated rectangular tank following the methodology similar to Warnitchai and Pinkaew (1998) and Tait (2008).

The porous baffles significantly reduce wave elevation in the varying filled levels of the tank compared to the baffle-free tank under the range of excitation frequencies. It is observed that the Reynolds number-dependent drag coefficient for porous baffles in the tank can significantly reduce the sloshing elevations and is found to be effective to achieve higher damping compared to the porosity-dependent drag coefficient for porous baffles in the sloshing tank. The analytical model’s response to free surface elevation variations in the sloshing tank was compared with the experiment’s test results. The analytical results matched with shake table test results with a quantitative difference near the first resonant frequency.

The scope of the study is limited to porous baffles performance under range sway motion and three different filling levels in the tank. The porous baffle performance includes Reynolds number dependent drag coefficient to explore the damping effect in the sloshing tank.

The porous baffles with low-level porosities in the sloshing tank have many engineering applications where the first resonant mode of sloshing in the tank is more important. The porous baffle drag coefficient is an important parameter to study the baffle’s damping effect in sloshing tanks. Hence, obtained analytical solution for liquid sloshing in the rectangular tank with Reynolds number as well as porosity-dependent drag coefficient (model 1) and porosity-dependent drag coefficient porous baffles (model 2) performance is discussed. The model’s test results were validated using a series of shake table sloshing experiments for three fill levels in the tank with sway motion at various excitation frequencies covering the first four sloshing resonant modes.

Three‐dimensional transient analysis of a submerged cylindrical shell is presented. Three‐dimensional trilinear eight‐noded isoparametric fluid element with pressure variable as unknown is coupled to a nine‐noded degenerate shell element. Staggered solution scheme is shown to be very effective for this problem. This allows significant flexibility in selecting an explicit or implicit integrator to obtain the solution in an economical way. Three‐dimensional transient analysis of the coupled shell fluid problem demonstrates that inclusion of bending mode is very important for submerged tube design—a factor which has not received attention, since most of the reported results are based on simplified two‐dimensional plane strain analysis.

The finite element method (FEM) is used to calculate the two-dimensional anti-plane dynamic response of structure embedded in D’Alembert viscoelastic multilayered soil on the rigid bedrock. This paper aims to research a time-domain absorbing boundary condition (ABC), which should be imposed on the truncation boundary of the finite domain to represent the dynamic interaction between the truncated infinite domain and the finite domain.

A high-order ABC for scalar wave propagation in the D’Alembert viscoelastic multilayered media is proposed. A new operator separation method and the mode reduction are adopted to construct the time-domain ABC.

The derivation of the ABC is accurate for the single layer but less accurate for the multilayer. To achieve high accuracy, therefore, the distance from the truncation boundary to the region of interest can be zero for the single layer but need to be about 0.5 times of the total layer height of the infinite domain for the multilayer. Both single-layered and multilayered numerical examples verify that the accuracy of the ABC is almost the same for both cases of only using the modal number excited by dynamic load and using the full modal number of infinite domain. Using the ABC with reduced modes can not only reduce the computation cost but also be more friendly to the stability. Numerical examples demonstrate the superior properties of the proposed ABC with stability, high accuracy and remarkable coupling with the FEM.

A high-order time-domain ABC for scalar wave propagation in the D’Alembert viscoelastic multilayered media is proposed. The proposed ABC is suitable for both linear elastic and D’Alembert viscoelastic media, and it can be coupled seamlessly with the FEM. A new operator separation method combining mode reduction is presented with better stability than the existing methods.