• Interaction and multiscale mechanics
• /
• v.2 no.4
• /
• pp.353-374
• /
• 2009

A new seamless multiscale simulation was developed for coupling the continuum model with its molecular dynamics. Kriging-based Finite Element Method (K-FEM) is employed to model the continuum base of the entire domain, while the molecular dynamics (MD) is confined in a localized domain of interest. In the coupling zone, where the MD domain overlaps the continuum model, the overall Hamiltonian is postulated by contributions from the continuum and the molecular overlays, based on a quartic spline scaling parameter. The displacement compatibility in this coupling zone is then enforced by the Lagrange multiplier technique. A multiple-time-step velocity Verlet algorithm is adopted for its time integration. The validation of the present method is reported through numerical tests of one dimensional atomic lattice. The results reveal that at the continuum/MD interface, the commonly reported spurious waves in the literature are effectively eliminated in this study. In addition, the smoothness of the transition from MD to the continuum can be significantly improved by either increasing the size of the coupling zone or expanding the nodal domain of influence associated with K-FEM.

• Bulletin of the Korean Mathematical Society
• /
• v.39 no.4
• /
• pp.589-606
• /
• 2002

We introduce and analyze a frequency-domain method for parabolic partial differential equations. The method is naturally parallelizable. After taking the Fourier transformation of given equations in the space-time domain into the space-frequency domain, we propose to solve an indefinite, complex elliptic problem for each frequency. Fourier inversion will then recover the solution in the space-time domain. Existence and uniqueness as well as error estimates are given. Fourier invertibility is also examined. Numerical experiments are presented.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.14 no.5
• /
• pp.1-12
• /
• 2010

This study presents a finite element analysis method that can accurately evaluate the nonlinear behaviour of structures affected by shallow soft subsoils and the soil-structure interaction. A two-dimensional finite element model that consists of a structure and shallow soft subsoil was used. The finite element model was used for a nonlinear time domain analysis of the OpenSees program. A parametric study was performed to investigate the effects of soil shear velocities, earthquake input motions, soft soil depth, and soil-structure interaction. The result of the proposed nonlinear finite element analysis method was compared with the result of an existing frequency domain analysis method, which is frequently used for addressing nonlinear soil behavior. The result showed that the frequency domain analysis, which uses equivalent secant soil stiffness and does not address the soil-structure interaction, significantly overestimated the response of the structures with short dynamic periods. The effect of the soil-structure interaction on the response spectrum did not significantly vary with the foundation dimensions and structure mass.

• Journal of the Korean Society of Safety
• /
• v.21 no.2 s.74
• /
• pp.143-149
• /
• 2006

This paper deals with the space-time finite element analysis of two-dimensional vibration problem with a single variable. The method of space-time finite elements enables the simpler solution than the usual finite element analysis with discretization in space only. We present a discretization technique in which finite element approximations are used in time and space simultaneously for a relatively large time period. The weighted residual process is used to formulate a finite element method for a space-time domain. A stability problem is described and some investigations for chosen type of rectangular space-time finite elements are carried out. Instability is caused by a too large time step of successive time steps in the traditional time-dependent problems. It has been shown that the numerical stability of time-stepping on the larger time steps is quite good. The unstructured space-time finite element not only overcomes the shortcomings of the stability in the traditional numerical methods, but it is also endowed with the features of an effective computational technique. Some numerical examples have been presented to illustrate the efficiency of the described method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.8
• /
• pp.1253-1260
• /
• 2001

Transient thermal analysis of a three-dimensional axisymmetric automotive disk brake is presented in this paper. Temperature fields are obtained using a hybrid FFT-FEM scheme that combines Fourier transform techniques and finite element method. The use of a fast Fourier transform algorithm can avoid singularity problems and lead to inexpensive computing time. The transformed problem is solved with finite element scheme for each frequency domain. Inverse transforms are then performed for time domain solution. Numerical examples are presented for validation tests. Comparisons with analytical results show very good agreement. Also, a 3-D simulation, based upon an automotive brake disk model is performed.

• Structural Engineering and Mechanics
• /
• v.15 no.6
• /
• pp.717-733
• /
• 2003

A time domain method is presented for soil-structure interaction analysis under seismic excitations. It is based on the finite element formulation incorporating infinite elements for the far field soil region. Equivalent earthquake input forces are calculated based on the free field responses along the interface between the near and far field soil regions utilizing the fixed exterior boundary method in the frequency domain. Then, the input forces are transformed into the time domain by using inverse Fourier transform. The dynamic stiffness matrices of the far field soil region formulated using the analytical frequency-dependent infinite elements in the frequency domain can be easily transformed into the corresponding matrices in the time domain. Hence, the response can be analytically computed in the time domain. A recursive procedure is proposed to compute the interaction forces along the interface and the responses of the soil-structure system in the time domain. Earthquake response analyses have been carried out on a multi-layered half-space and a tunnel embedded in a layered half-space with the assumption of the linearity of the near and far field soil region, and results are compared with those obtained by the conventional method in the frequency domain.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.38-47
• /
• 2005

Substructuring methods are often used in finite element structural analyses. In this study a multi-level substructuring(MLSS) algorithm is developed and proposed as a possible candidate for finite element fluid solvers. The present algorithm consists of four stages such as a gathering, a condensing, a solving and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At the highest level, each sub-domain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each sub-domain has been replaced by a sequential static condensation of gathered element matrices. The global algebraic system arising from the assembly of each sub-domain matrices is solved using a well-known iterative solver such as the conjugare gradient(CG) or the conjugate gradient squared(CGS) method. A time comparison with CG has been performed on a 2-D Poisson problem. With one domain the computing time by MLSS is comparable with that by CG up to about 260,000 d.o.f. For 263,169 d.o.f using 8 x 8 sub-domains, the time by MLSS is reduced to a value less than $30\%$ of that by CG. The lid-driven cavity problem has been solved for Re = 3200 using the element interpolation degree(Deg.) up to cubic. in this case, preconditioning techniques usually accompanied by iterative solvers are not needed. Finite element formulation for the incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et al.[9].

• Proceedings of the KIEE Conference
• /
• 1999.11b
• /
• pp.450-454
• /
• 1999

In order to evaluate the integrity of nuclear power plants, J-integral calculation is crucial. For this purpose, finite element method is popularly used to obtain J-integral. However, high cost time consuming preprocess should be performed to design the finite element model of a cracked structure. Also, the J-integral should be verified by alternative method since it may differ depending on the calculation method. The objective of this paper is to develop a three-dimensional elastic-plastic J-integral analysis system which is named as EPAS. The EPAS program consists of an automatic mesh generator for a through-wall crack and a surface crack, a solver based on ABAQUS program, and a J-integral calculation program which provides DI(Domain Integral) and EDI(Equivalent Domain Integral) based J-integral calculation. Using the EPAS program, an optimized finite element model for a cracked structure can be generated and corresponding J-integral can be obtained subsequently.

• Journal of the Korean Magnetics Society
• /
• v.8 no.3
• /
• pp.161-168
• /
• 1998

In this paper, a algorithm is proposed, which is applicable to the transient analysis of diffusion equations by combined use of the Laplace transform and the finite element method. The proposed method removes the time terms using the Laplace transform and then solves the associated equation with the finite element method. The solution which is solved at frequency domain is transformed into time domain by use of the Laplace inversion. To verify the proposed algorithm, a heat conduction problem is analysed. And the solution showed a good agreement with analytic solution. Because the time-step method is not needed, the proposed method is very useful in solving various kinds of diffusion equations.

• Steel and Composite Structures
• /
• v.10 no.6
• /
• pp.541-561
• /
• 2010

This paper presents finite element analyses, experimental measurements and finite element model updating of an arch type steel laboratory bridge model using semi-rigid connections. The laboratory bridge model is a single span and fixed base structure with a length of 6.1 m and width of 1.1m. The height of the bridge column is 0.85 m and the maximum arch height is 0.95 m. Firstly, a finite element model of the bridge is created in SAP2000 program and analytical dynamic characteristics such as natural frequencies and mode shapes are determined. Then, experimental measurements using ambient vibration tests are performed and dynamic characteristics (natural frequencies, mode shapes and damping ratios) are obtained. Ambient vibration tests are performed under natural excitations such as wind and small impact effects. The Enhanced Frequency Domain Decomposition method in the frequency domain and the Stochastic Subspace Identification method in the time domain are used to extract the dynamic characteristics. Then the finite element model of the bridge is updated using linear elastic rotational springs in the supports and structural element connections to minimize the differences between analytically and experimentally estimated dynamic characteristics. At the end of the study, maximum differences in the natural frequencies are reduced on average from 47% to 2.6%. It is seen that there is a good agreement between analytical and experimental results after finite element model updating. Also, connection percentages of the all structural elements to joints are determined depending on the rotational spring stiffness.