• Advances in aircraft and spacecraft science
• /
• v.1 no.2
• /
• pp.145-175
• /
• 2014

The results of a series of numerical experiments are presented to verify some of the important developments made in the first part of this paper. Firstly, the static solution of an algebraic system obtained through Strong Formulation Finite Element Method (SFEM) is presented. Secondly, the stress and strain recovery procedure is descripted for the present technique. It will be clear that the present approach is suitable for any strong formulation finite element methodology, due to the presented general approach based on the unknown displacements and on the elasticity equations. Thirdly, the numerical solutions for some classical and other numerical results found in literature are exposed. Finally, an arbitrarily shaped composite plate is solved and good agreement is observed for all the presented cases.

• Advances in aircraft and spacecraft science
• /
• v.1 no.2
• /
• pp.125-143
• /
• 2014

This paper provides a new technique for solving the static analysis of arbitrarily shaped composite plates by using Strong Formulation Finite Element Method (SFEM). Several papers in literature by the authors have presented the proposed technique as an extension of the classic Generalized Differential Quadrature (GDQ) procedure. The present methodology joins the high accuracy of the strong formulation with the versatility of the well-known Finite Element Method (FEM). The continuity conditions among the elements is carried out by the compatibility or continuity conditions. The mapping technique is used to transform both the governing differential equations and the compatibility conditions between two adjacent sub-domains into the regular master element in the computational space. The numerical implementation of the global algebraic system obtained by the technique at issue is easy and straightforward. The main novelty of this paper is the application of the stress and strain recovery once the displacement parameters are evaluated. Computer investigations concerning a large number of composite plates have been carried out. SFEM results are compared with those presented in literature and a perfect agreement is observed.

• Structural Engineering and Mechanics
• /
• v.11 no.2
• /
• pp.221-236
• /
• 2001

A local point interpolation method (LPIM) is presented for the stress analysis of two-dimensional solids. A local weak form is developed using the weighted residual method locally in two-dimensional solids. The polynomial interpolation, which is based only on a group of arbitrarily distributed nodes, is used to obtain shape functions. The LPIM equations are derived, based on the local weak form and point interpolation. Since the shape functions possess the Kronecker delta function property, the essential boundary condition can be implemented with ease as in the conventional finite element method (FEM). The presented LPIM method is a truly meshless method, as it does not need any element or mesh for both field interpolation and background integration. The implementation procedure is as simple as strong form formulation methods. The LPIM has been coded in FORTRAN. The validity and efficiency of the present LPIM formulation are demonstrated through example problems. It is found that the present LPIM is very easy to implement, and very robust for obtaining displacements and stresses of desired accuracy in solids.

• Korean Journal of Mathematics
• /
• v.18 no.4
• /
• pp.425-440
• /
• 2010

Using Green's theorem, elliptic boundary value problems can be converted to boundary integral equations. A numerical methods for boundary integral equations are boundary elementary method(BEM). BEM has advantages over finite element method(FEM) whenever the fundamental solutions are known. Helmholtz type equations arise naturally in many physical applications. In a boundary integral formulation for the exterior Neumann there occurs a hypersingular operator which exhibits a strong singularity like $\frac{1}{|x-y|^3}$ and hence is not an integrable function. In this paper we are going to remove this hypersingularity by reducing the regularity of test functions.

• Earthquakes and Structures
• /
• v.17 no.3
• /
• pp.257-270
• /
• 2019

The purpose of this work is to analyze the bending behavior of laminated composite plates using the refined fourvariable theory and the finite element method approach using an ANSYS 12 computational code. The analytical model is based on the multilayer plate theory of shear deformation of the nth-order proposed by Xiang et al 2011 using the theory principle developed by Shimpi and Patel 2006. Unlike other theories, the number of unknown functions in the present theory is only four, while five or more in the case of other theories of shear deformation. The formulation of the present theory is based on the principle of virtual works, it has a strong similarity with the classical theory of plates in many aspects, it does not require shear correction factor and gives a parabolic description of the shear stress across the thickness while filling the condition of zero shear stress on the free edges. The analysis is validated by comparing results with those in the literature.

• Steel and Composite Structures
• /
• v.47 no.3
• /
• pp.365-374
• /
• 2023

This paper presents a new scheme for constructing locking-free finite elements in thick and thin plates, called Discrete Shear Gap element (DSG), using multiphase material topology optimization for triangular elements of Reissner-Mindlin plates. Besides, common methods are also presented in this article, such as quadrilateral element (Q4) and reduced integration method. Moreover, when the plate gets too thin, the transverse shear-locking problem arises. To avoid that phenomenon, the stabilized discrete shear gap technique is utilized in the DSG3 system stiffness matrix formulation. The accuracy and efficiency of DSG are demonstrated by the numerical examples, and many superior properties are presented, such as being a strong competitor to the common kind of Q4 elements in the static topology optimization and its computed results are confirmed against those derived from the three-node triangular element, and other existing solutions.

• Steel and Composite Structures
• /
• v.48 no.1
• /
• pp.73-88
• /
• 2023

This paper presents a novel implicit level set method for topology optimization of functionally graded (FG) structures with pre-existing discontinuities (pre-cracks) using radial basis functions (RBF). The mathematical formulation of the optimization problem is developed by incorporating RBF-based nodal densities as design variables and minimizing compliance as the objective function. To accurately capture crack-tip behavior, crack-tip enrichment functions are introduced, and an eXtended Finite Element Method (X-FEM) is employed for analyzing the mechanical response of FG structures with strong discontinuities. The enforcement of boundary conditions is achieved using the Hamilton-Jacobi method. The study provides detailed mathematical expressions for topology optimization of systems with defects using FG materials. Numerical examples are presented to demonstrate the efficiency and reliability of the proposed methodology.

• Structural Engineering and Mechanics
• /
• v.87 no.2
• /
• pp.173-189
• /
• 2023

For a given structural geometry, the stiffness and damping parameters of the soil and the dynamic response of the structure may change in the face of an equivalent linear soil behavior caused by a strong earthquake. Therefore, the influence of equivalent linear soil behavior on the impedance functions form and the seismic response of the soil-structure system has been investigated. Through the substructure method, the seismic response of the selected structure was obtained by an analytical formulation based on the dynamic equilibrium of the soil-structure system modeled by an analog model with three degrees of freedom. Also, the dynamic response of the soil-structure system for a nonlinear soil behavior and for the two types of impedance function forms was also analyzed by 2D finite element modeling using ABAQUS software. The numerical results were compared with those of the analytical solution. After the investigation, the effect of soil nonlinearity clearly showed the critical role of soil stiffness loss under strong shaking, which is more complex than the linear elastic soil behavior, where the energy dissipation depends on the seismic motion amplitude and its frequency, the impedance function types, the shear modulus reduction and the damping increase. Excellent agreement between finite element analysis and analytical results has been obtained due to the reasonable representation of the model.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.5
• /
• pp.829-837
• /
• 1992

This paper presents a new approach for the analysis of hot strip rolling processes. The approach is based on the finite element method and capable of predicting velocity field in the strip, temperature field in the strip, temperature field in the roll, and roll pressure. Basic finite element formulations for heat transfer analysis are described with emphasis on the treatment of numerical instability resulting from a standard Galerkin formulation. Comparison with the theoretical solutions found in the literature is made for the evaluation of the accuracy of the temperature solutions. An iterative scheme is developed for dealing with strong correlations between the metal flow characteristics and the thermal behavior of the roll-strip system. A series of process simulations are carried out to investigate the effect of various process parameters including interface friction, interface heat transfer coefficient, roll speed, reduction in thickness, and spray zone. The results are shown and discussed.

• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.1
• /
• pp.19-29
• /
• 2013

This paper deals with the study of buckling, vibration, and parametric instability characteristics in a damaged cross-ply and angle-ply laminated plate like beam under in-plane harmonic loading, using the finite element approach. Damage is modelled using an anisotropic damage formulation, based on the concept of reduction in stiffness. The effect of damage on free vibration and buckling characteristics of a thin plate like beam has been studied. It has been observed that damage shows a strong orthogonality and in general deteriorates the static and dynamic characteristics. For the harmonic type of loading, analysis was carried out on a thin plate like beam by solving the governing differential equation which is of Mathieu-Hill type, using the method of multiple scales (MMS). The effects of damage and its location on dynamic stability characteristics have been presented. The results indicate that, compared to the undamaged plate like beam, heavily damaged beams show steeper deviations in simple and combination resonance characteristics.