• Interaction and multiscale mechanics
• /
• v.3 no.1
• /
• pp.95-110
• /
• 2010

In recent years, study of the static response of pavements to moving vehicle and aircraft loads has received significant attention because of its relevance to the design of pavements and airport runways. The static response of beams resting on an elastic foundation and subjected to moving loads was studied by several researchers in the past. However, most of these studies were limited to steady-state analytical solutions for infinitely long beams resting on Winkler-type elastic foundations. Although the modelling of subgrade as a continuum is more accurate, such an approach can hardly be incorporated in analysis due to its complexity. In contrast, the two-parameter foundation model provides a better way for simulating the underlying soil medium and is conceptually more appealing than the one-parameter (Winkler) foundation model. The finite element method is one of the most suitable mathematical tools for analysing rigid pavements under moving loads. This paper presents an improved solution algorithm based on the finite element method for the static analysis of rigid pavements under moving vehicular or aircraft loads. The concrete pavement is discretized by finite and infinite beam elements, with the latter for modelling the infinity boundary conditions. The underlying soil medium is modelled by the Pasternak model allowing the shear interaction to exist between the spring elements. This can be accomplished by connecting the spring elements to a layer of incompressible vertical elements that can deform in transverse shear only. The deformations and forces maintaining equilibrium in the shear layer are considered by assuming the shear layer to be isotropic. A parametric study is conducted to investigate the effect of the position of moving loads on the response of pavement.

• International Journal of Aeronautical and Space Sciences
• /
• v.13 no.4
• /
• pp.446-457
• /
• 2012

Aeroelastic analysis of an aircraft with a high aspect ratio wing for medium altitude and long endurance capability was attempted in this paper. In order to achieve such an objective, various structural models were adopted. The traditional approach has been based on a one-dimensional Euler-Bernoulli beam model. The structural analysis results of the present beam model were compared with those by the three-dimensional NASTRAN finite element model. In it, a taper ratio of 0.5 was applied; it was comprised of 21 ribs and 3 spars, and included two control surfaces. The relevant unsteady aerodynamic forces were obtained by using ZAERO, which is based on the doublet lattice method that considers flow compressibility. To obtain the unsteady aerodynamic force, the structural mode shapes and natural frequencies were transferred to ZAERO. Two types of unsteady aerodynamic forces were considered. The first was the unsteady aerodynamic forces which were based on the one-dimensional beam shape; the other was based on the three-dimensional FEM model shape. These two types of aerodynamic forces were compared, and applied to the foregoing flutter analysis. The ultimate goal of the present research is to analyze the possible interaction between the rigid-body degrees of freedom and the aeroelastic modes. This will be achieved after the development of a reliable nonlinear beam formulation that would validate the current results as well as enable a thorough investigation of the nonlinearity. Moreover, such analysis will allow for an examination of the above-mentioned interaction between the flight dynamics and aeroelastic modes with the inclusion of the rigid body degrees of freedom.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.6
• /
• pp.519-525
• /
• 2009

In this study a cohesive zone model was used to simulate the delamination phenomena which occurs by a successive crack initiation and propagation in composite laminates. The cohesive zone model was incorporated to the classical finite element method via cohesive element formulation and then implemented into the user-subroutine UEL of a commercial finite element program Abaqus. To validate the formulation and implementation of the cohesive element the finite element results were compared with the experimental data of double cantilever beam and end notched flexure tests. The numerical results well agree with the experimental load-displacement curves. Also the effect of the elastic stiffness and the size of the cohesive element on the global load-displacement curves were studied numerically. To minimize the mesh-dependency of the crack propagation path and eliminate the zig-zag patterns in the load-displacement curve, cohesive elements should be refined at the crack-tip.

• Steel and Composite Structures
• /
• v.24 no.4
• /
• pp.425-439
• /
• 2017

As a thin-walled structure, local joint flexibility (LJF) in a tubular structure is prominent, and it may produce significant effect on the static performance for the overall structure. This study presents a simplified analytical model to analyze the static behavior for a steel tubular structure with LJF. The presented model simplifies a tubular structure into a frame model consisted of beam elements with considering the LJFs at the connections between any two elements. Theoretical equations of the simplified analytical model are deduced. Through comparison with 3-D finite element results of two typical planar tubular structures consisted of T- and Y-joints respectively, the presented method is proved to be accurate. Furthermore, the effect of LJF on the overall performance of the two tubular structures (including the deflection and the internal forces) is also investigated, and it is found from analyses of internal forces and deformation that a rigid connection assumption in a frame model by using beam elements in finite element analysis can provide unsafe and inaccurate estimation.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.15 no.3
• /
• pp.8-16
• /
• 2006

This paper discusses design sensitivity analysis and its application to a structural dynamics modification. Eigenvalue derivatives are determined with respect to the element parameters, which include intrinsic property parameters such as Young's modulus, density of the material, diameter of a beam element, thickness of a plate element, and shape parameters. Derivatives of stiffness and mass matrices are directly calculated by derivatives of element matrices. The first and the second order derivatives of the eigenvalues are then mathematically derived from a dynamic equation of motion of FEM model. The calculation of the second order eigenvalue derivative requires the sensitivity of its corresponding eigenvector, which are developed by Nelson's direct approach. The modified eigenvalue of the structure is then evaluated by the Taylor series expansion with the first and the second derivatives of eigenvalue. Numerical examples for simple beam and plate are presented. First, eigenvalues of the structural system are numerically calculated. Second, the sensitivities of eigenvalues are then evaluated with respect to the element intrinsic parameters. The most effective parameter is determined by comparing sensitivities. Finally, we predict the modified eigenvalue by Taylor series expansion with the derivatives of eigenvalue for single parameter or multi parameters. The examples illustrate the effectiveness of the eigenvalue sensitivity analysis for the optimization of the structures.

• Steel and Composite Structures
• /
• v.6 no.2
• /
• pp.139-157
• /
• 2006

In this paper, a cyclic mechanical model is presented to simulate the behaviour of top and seat with web angle beam-to-column connections. The introduced mechanical model is compared with Eurocode 3 Annex J, its extension, and with experimental data. To have a better insight regarding the actual response of the joints, available results of the experiments, carried out on full-scale top and seat angle joints under monotonic and cyclic loading, are first considered. Subsequently, a finite element model of the test setup is developed. The application of the proposed model, its comparisons with the experimental curves and with the Eurocode 3 Annex J and with its modification, clearly show the excellent quality of the model proposed.

• Journal of KSNVE
• /
• v.6 no.3
• /
• pp.287-296
• /
• 1996

A virtual work form of flexible multibody dynamic formulation with rotary inertia has been derived. For the analysis of large flexible multibody systems, deformation modal coordinates have been employed to represent coupled motion between gross and vibrational motion. For the efficient evaluation of the entries in the mass matrix, a flexible body has been treated as a collection of mass points. The rotary inertia was generated from the consistent mass matrix in a finite element model. Deformation mode shapes were obtained from finite element analysis. Bending and twisting vibration analyses of a cantilever have been carried out to see rotary inertia effects. A space flexible robot simulation has been also carried out to show effectiveness of the proposed formulation. This formulation is effective to the model that consists of beam, plate, or shell element that contains rotational degree of freedom at the nodal point. It is also effective to the flexible body model to which a large lumped rotary inertia is attached.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.05a
• /
• pp.304-308
• /
• 2006

The focus of this study is modeling technique for a bellows in vehicle exhaust system. Bellows was developed using tile finite element model by replacing with the equivalent beam. The equivalent beam model were studied in detail. Non-structural node in the cross section of original model is given to expressing their motion. Equivalent mass matrix and stiffness matrix calculated using Guyan reduction method. Material Properties of beam was obtained from the direct comparison between equivalent model and that of Timoshenko beam model. The calculated natural frequencies and mode shape are compared with the reference results and coincided well. The results were compared with the confirmed results, which were in good agreement.

• Advances in aircraft and spacecraft science
• /
• v.10 no.2
• /
• pp.179-202
• /
• 2023

This article presents a numerical analysis to investigate the natural frequencies and harmonic response of a perforated cantilever beam attached to two layers of piezoelectric materials by using the finite element method for the first time. The bimorph piezoelectric is composed of 3 layers; two of them at the outer are piezoelectric, and the inner isotropic material. A higher order 3-D 20-node solid element that exhibits quadratic displacement behavior is exploited to discretize the isotropic layer, and coupled piezoelectric 3D element with twenty nodes is used to mesh the top and bottom layers. CIRCU94 element is added to act as a resistor part of the model. The proposed model is validated with previous works. The numerical parametric studies are presented to illustrate the effects of perforation geometry, the number of rows, the resistance on the natural frequencies, frequency response, and power. It is found that the thickness has a positive relationship with the natural frequency. Perforations help in producing higher voltage, and the best shape is rectangular perforations, and to produce higher voltage, two rows of rectangular perforations should be applied.

• Advances in aircraft and spacecraft science
• /
• v.2 no.1
• /
• pp.57-76
• /
• 2015

We investigated complex damped modes in beams in the presence of a viscoelastic layer sandwiched between two elastic layers. The problem was solved using two approaches, (1) Rayleigh beam theory and analyzed using the Ritz method, and (2) by using 2D plane stress elasticity based finite-element method. The damping in the layers was modeled using the complex modulus. Simply-supported, cantilever, and viscously supported boundary conditions were considered in this study. Simple trigonometric functions were used as admissible functions in the Ritz method. The key idea behind sandwich structure is to increase damping in a beam as affected by the presence of a highly-damped core layer vibrating mainly in shear. Different assumptions are utilized in the literature, to model shear deformation in the core layer. In this manuscript, we used FEM without any kinematic assumptions for the transverse shear in both the core and elastic layers. Moreover, numerical examples were studied, where the base and constraining layers were also damped. The loss factor was calculated by modal strain energy method, and by solving a complex eigenvalue problem. The efficiency of the modal strain energy method was tested for different loss factors in the core layer. Complex mode shapes of the beam were also examined in the study, and a comparison was made between viscoelastically and viscously damped structures. The numerical results were compared with those available in the literature, and the results were found to be satisfactory.