• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.173-195
• /
• 2013

Meshfree methods are known to have the capability to overcome the strict regularization requirements and numerical instabilities that encumber the finite element method (FEM) in large deformation problems. They are also more naturally suited for problems involving material perforation and fragmentation. To take advantage of the high efficiency of FEM and high accuracy of meshfree methods, a coupled finite element (FE) and reproducing kernel (RK, one of the meshfree approximations) formulation is described in this paper. The coupling of FE and RK approximation is implemented in an evolutionary fashion, where the extent and location of the evolution is dependent on a triggering criteria provided by the material constitutive laws. To enhance computational efficiency, Gauss quadrature is applied to integrate both FE and RK domains so that no state variable transfer is required when mesh conversion is performed. To control the hourglassing that might occur with 1-point integrated hexahedral grids, viscous type hourglass control is implemented. Meanwhile, the FEM version of the K&C concrete (KCC) model was modified to make it applicable in both FE and RK formulations. Results using this code and the KCC model are shown for the modeling of concrete responses under quasi-static, blast and impact loadings. These analyses demonstrate that fragmentation phenomena of the sort commonly observed under blast and impact loadings of concrete structures was able to be realistically captured by the coupled formulation.

• Structural Engineering and Mechanics
• /
• v.82 no.3
• /
• pp.401-416
• /
• 2022

The aim of this study is to examine the frictionless double receding contact problem for two functionally graded (FG) layers pressed with a uniformly distributed load and resting on a homogeneous half plane (HP) using analytical and numerical methods. The FG layers are made of a non-homogeneous material with an isotropic stress-strain law with exponentially varying properties. It is assumed that the contact at the FG layers and FG layer-HP interface is frictionless. The body force of the FG layers and homogeneous HP are ignored in the study. Firstly, an analytical solution for the contact problem has been realized using the theory of elasticity and the Fourier integral transform techniques. Then, the problem modeled and two-dimensional analysis was carried out by using the ANSYS package program based on FEM. Numerical results for contact lengths and contact pressures between FG layers and FG layer-HP were provided for various dimensionless quantities including material inhomogeneity, distributed load width, the shear module ratio, and the heights of the FG layers for both methods. The results obtained using FEM were compared with the results found using the analytical formulation. It was found that the results obtained from analytical formulation were in perfect agreement with the FEM study.

• Computers and Concrete
• /
• v.27 no.3
• /
• pp.199-210
• /
• 2021

The aim of this paper was to examine the continuous and discontinuous contact problems between the functionally graded (FG) layer pressed with a uniformly distributed load and homogeneous half plane using an analytical method and FEM. The FG layer is made of non-homogeneous material with an isotropic stress-strain law with exponentially varying properties. It is assumed that the contact at the FG layer-half plane interface is frictionless, and only the normal tractions can be transmitted along the contacted regions. The body force of the FG layer is considered in the study. The FG layer was positioned on the homogeneous half plane without any bonds. Thus, if the external load was smaller than a certain critical value, the contact between the FG layer and half plane would be continuous. However, when the external load exceeded the critical value, there was a separation between the FG layer and half plane on the finite region, as discontinuous contact. Therefore, there have been some steps taken in this study. Firstly, an analytical solution for continuous and discontinuous contact cases of the problem has been realized using the theory of elasticity and Fourier integral transform techniques. Then, the problem modeled and two-dimensional analysis was carried out by using ANSYS package program based on FEM. Numerical results for initial separation distance and contact stress distributions between the FG layer and homogeneous half plane for continuous contact case; the start and end points of separation and contact stress distributions between the FG layer and homogeneous half plane for discontinuous contact case were provided for various dimensionless quantities including material inhomogeneity, distributed load width, the shear module ratio and load factor for both methods. The results obtained using FEM were compared with the results found using analytical formulation. It was found that the results obtained from analytical formulation were in perfect agreement with the FEM study.

• 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.

• Structural Engineering and Mechanics
• /
• v.41 no.6
• /
• pp.735-750
• /
• 2012

This paper presents a comparison between two different procedures to deal with the geometric nonlinear analysis of space trusses, considering its structural stability aspects. The first nonlinear formulation, called positional, uses nodal positions rather than nodal displacements to describe the finite elements kinematics. The strains are computed directly from the proposed position concept, using a Cartesian coordinate system fixed in space. The second formulation, called corotational, is based on the explicit separation between rigid body motion and deformed motion. The numerical examples demonstrate the performances and the convergence of the responses for both analyzed formulations. Two numerical examples were compared, including a lattice beam with postcritical behavior. Despite the two completely different approaches to deal with the geometrical nonlinear problem, the results present good agreement.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.05a
• /
• pp.1026-1031
• /
• 2002

This paper shows the effect of mean flow in duct. The potential equation of duct with mean flow is obtained. A finite element method(FEM) is used to predict acoustic performance of duct with mean flow. The formulation of the finite element method is derived for duct taking into consideration of the convective effect of mean flow. A simple duct, simple expansion chamber and a duct with resonator are implemented to show the effects of the mean flow.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.2
• /
• pp.436-444
• /
• 1991

A computational method has been developed to analyze the bone-remodeling induced by external fixator. The method was based on the Finite Element Method (FEM) in combination with numerical formulation of adaptive bone-remodeling theories. As a feed-back control variable, compressive strain and effective stress were used to determine the surface remodeling and internal (density) remodeling respectively. Surface remodeling and internal remodeling were combined at each time step to predict the rel situation. A noticeable shape and density change were detected at the region between two pins and density change was decreased with time increment. At final time step, the shape and density distribution were converged closely to its original intact bone model. Similar change was detected in stress distribution. The altered stress distribution due to the pin and external fixator converged to the intact stress distribution with time.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.48 no.10
• /
• pp.523-529
• /
• 1999

In the case of the large power transformer, the grain-oriented material is usually used. So, to obtain more accurate results, anisotropy and non-linearity of the material must be considered. The Newton-Raphson(NR) method is generally used for analyzing these non-linear properties, but it consumes so much time, especially when the number of nodes is large or the shape of the model is complex. The transmission line modeling (TLM) method is successfully adopted to the analysis of non-linear properties with FEM, but it has not been adopted to the analysis of the anisotropic material. In this paper, the formulation of the TLM method considering anisotropy is developed and the adoption to the 3-phase transformer is presented.

• Structural Engineering and Mechanics
• /
• v.72 no.2
• /
• pp.217-227
• /
• 2019

The methodology known as "shape sensing" allows the reconstruction of the displacement field of a structure starting from strain measurements, with considerable implications for structural monitoring, as well as for the control and implementation of smart structures. An approach to shape sensing is based on the inverse Finite Element Method (iFEM) that uses a variational principle enforcing a least-squares compatibility between measured and analytical strain measures. The structural response is reconstructed without the knowledge of the mechanical properties and load conditions but based only on the relationship between displacements and strains. In order to efficiently apply iFEM to the most common structural typologies of civil engineering, its formulation according to the kinematical assumptions of the Bernoulli-Euler theory is presented. Two beam inverse finite elements are formulated for different loading conditions. Depending on the type of element, the relationship between the minimum number of required measurement stations and the interpolation order is defined. Several examples representing common applications of civil engineering and involving beams and frames are presented. To simulate the experimental strain data at the station points and to verify the accuracy of the displacements obtained with the iFEM shape sensing procedure, a direct FEM analysis of the considered structures is performed using the LUSAS software.

• Structural Engineering and Mechanics
• /
• v.48 no.6
• /
• pp.775-789
• /
• 2013

The current paper proposes a boundary element formulation, applicable to 2-D and 3-D elastostatics problems using a unified approach for transformations of the domain integrals into boundary integrals. The method is applicable to linear problems encompassing both finite and infinite multi-region domains allowing non-vanishing body forces. Numerical results agree quite well with the analytical solutions; while the present method offers easy formulation with less numerical efforts in comparison to FEM or some BEM which need interior points to treat arbitrary body forces. It is demonstrated that the method has the potential to have profound impact on engineering design, notably in dam-foundation interaction.