• Earthquakes and Structures
• /
• v.12 no.2
• /
• pp.251-262
• /
• 2017

The present research intends to study the effects of the seismic soil-foundation-structure interaction (SFSI) on the dynamic response of various buildings. Two methods including direct and Cone model were studied through 3D finite element method using ABAQUS software. Cone model as an approximate method to consider the SFSI phenomenon was developed and evaluated for both high and low rise buildings. Effect of soil nonlinearity, foundation rigidity and embedment as well as friction coefficient between soil-foundation interfaces during seismic excitation are investigated. Validity and performance of both approaches are evaluated as reference graphs for Cone model and infinite boundary condition, soil nonlinearity and amplification factor for direct method. A series of calculations by DeepSoil for inverse earthquake record modification was conducted. A comparison of the two methods was carried out by root-mean-square-deviation (RMSD) tool for maximum lateral displacement and story shear forces which verifies that Cone model results have good agreement with direct method. It was concluded that Cone method is a convenient, fast and rather accurate method as an approximate way to count for soil media.

• Journal of KSNVE
• /
• v.7 no.4
• /
• pp.637-644
• /
• 1997

In this paper, a method of finite element analysis is presented for efficient calculation of vibration characteristics of not only simply interconnected structure with cyclic symmetry but also multiply interconnected structure with cyclic symmetry by using discrete Fourier trandform by means of a computer with small memory in a short time. Simply interconnected structure means it is composed of substructures which are adjacent themselves in circumferential direction. First, a mathematical model of multiply interconnected structure with cyclic symmetry is defined. The multiply interconnected structure is partitioned into substructures with the same goemetric configuration and constraint eqauations to be satisfied on connecting boundaries are defined. Nodal displacements and forces are transformed into complex forms through discrete Fourier transform and then finite element analysis is performed for just only a representative substructure. In free vibration analysis, natural frequencies of a whole structure can be obtained through a series of calculation for a substructure along the number of nodal diameter. And in forced vibration analysis, forced response of whole structure can be achieved by using inverse discrete Fourier transform of results which come from analysis for a substructure.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.58 no.10
• /
• pp.1936-1941
• /
• 2009

This paper deals with a new methodology for topology optimization in which the topology of the design domain may change during the shape optimization process. To achieve this, the concept of the topological gradient is introduced to compute the sensitivity of an objective function when a small hole is drilled in the domain. Based on shape and topological sensitivity values, the shape and topology of the design domain may be simultaneously changed during design iterations if necessary. To verify the advantages and also to facilitate understanding of the method itself, two electrostatic design problems have been tested by using 2D finite element analysis: the first is the inverse problem of a simple dielectric model and the second is the rotor design of a MEMS actuator.

• Journal of Biomedical Engineering Research
• /
• v.20 no.1
• /
• pp.45-51
• /
• 1999

The purpose of this study is to investigate the effects of various factors in keratotomy for astigmatism correction on surgical outcomes by finite element method as well as animal experiments. Three kinds of surgical techniques were mechanically investigated : arcuate, straight, and inverse arcuate keratotomy. Among the three techniques the arcuate keratotomy is the most popular one while the other two techniques are being investigated in this area. The arcuate keratotomy was found to be more controllable and effective in reducing the refractive power than the others. In arcuate keratotomy it was found most effective when the incision was located in the middle position between the apex and the edge of the cornea from the results of experiment as well as finite element study. Regarding to the range of the corneal incision in arcuate keratotomy, the incision angle of 90$^{\circ}$ was found th be most effective in reducing refractive power than other angles even it was incised up to 150$^{\circ}$. Therefore, it was concluded that 90$^{\circ}$ of incision angle results in the largest decrease in refractive power in arcuate keratotomy. However, other important findings were that the effect of the surgery decreased with time so the visco-effect of the cornea and auto-healing process. Therefore, these factors should be considered in future studies.

• Structural Engineering and Mechanics
• /
• v.67 no.4
• /
• pp.403-415
• /
• 2018

In the present study, nonlinear dynamic response of polymer-CNT-fiber multiscale nanocomposite plate resting on elastic foundations in thermal environments using the finite element method is performed. In this regard, the governing equations are derived based on Inverse Hyperbolic Shear Deformation Theory and von $K{\acute{a}}rm{\acute{a}}n$ geometrical nonlinearity. Three type of distribution of temperature through the thickness of the plate namely, uniform linear and nonlinear are considered. The considered element is C1-continuous with 15 DOF at each node. The effective material properties of the multiscale composite are calculated using Halpin-Tsai equations and fiber micromechanics in hierarchy. The carbon nanotubes are assumed to be uniformly distributed and randomly oriented through the epoxy resin matrix. Five types of impulsive loads are considered, namely the step, sudden, triangular, half-sine and exponential pulses. After examining the validity of the present work, the effects of the weight percentage of SWCNTs and MWCNTs, nanotube aspect ratio, volume fraction of fibers, plate aspect, temperature, elastic foundation parameters, distribution of temperature and shape of impulsive load on nonlinear dynamic response of CNT reinforced multi-phase laminated composite plate are studied in details.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.8
• /
• pp.41-46
• /
• 2006

The capability for reconstructing impact forces of an aircraft wing using impact response functions and regularization methods were examined. The impact response function which expresses the relation between the structure response and the impact force was derived using the information on mass and stiffness data of a finite element model for the wing. Iterative Tikhonov regularization method and generalized singular value decomposition method were used to inverse the impact response function that was generally ill-posed. For the numerical verification, a fighter aircraft wing was used. Strain and deflection histories obtained from finite element analysis were compared with the results calculated using impact response functions. And the impact forces were reconstructed with the strain histories obtained from finite element analysis. The numerical verification results showed that this method can be used to monitor impact forces on aircraft structures.

• Transactions of Materials Processing
• /
• v.20 no.7
• /
• pp.473-478
• /
• 2011

In this study, the performance of a thermo siphon type radiator made of copper for LED lighting system was evaluated by using an inverse heat transfer method. Heating experiments and finite element heat transfer analysis were conducted for three different cases. The data obtained from experiments were compared with the analysis results. Based on the data obtained from experiments, the inverse heat transfer method was used in order to evaluate the heat transfer coefficient. First, the heat transfer analysis was conducted for non-vacuum state, without the refrigerant. The evaluated heat transfer coefficient on the radiator surface was 40W/$m^2^{\circ}C$. Second, the heat transfer analysis was conducted for non-vacuum state, with the refrigerant, resulting in the heat transfer coefficient of 95W/$m^2^{\circ}C$. Third, the heat transfer analysis was conducted for vacuum state, with refrigerant. For the third case, the evaluated heat transfer coefficients were 140W/$m^2^{\circ}C$. Third, the heat transfer analysis was conducted for vacuum state, with refrigerant. For the third case, the evaluated heat transfer coefficients were 140W/$m^2^{\circ}C$ for the radiator body, 5W/$m^2^{\circ}C$. Third, the heat transfer analysis was conducted for vacuum state, with refrigerant for the rising position of radiator pipe, 35W/$m^2^{\circ}C$. Third, the heat transfer analysis was conducted for vacuum state, with refrigerant. For the highest position of radiator pipe, and 120W/$m^2^{\circ}C$ for the downturn position of radiator pipe. As a result of inverse heat transfer analysis, it was confirmed that the thermal performance of the current radiator was best in the case of the vacuum state using the refrigerant.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.20 no.10
• /
• pp.923-930
• /
• 2010

Finite element models of dynamic systems can be updated in two stages. In the first stage, mass and stiffness matrices are updated neglecting damping, and in the second stage, damping matrices are estimated with the mass and stiffness matrices fixed. Three methods to estimate damping matrices for this purpose are proposed in this paper. The methods include one for proportional damping systems and two for non-proportional damping systems. Method 1 utilizes orthogonality of normal modes and estimates damping matrices using the modal parameters extracted from the measured responses. Method 2 estimates damping matrices from impedance matrices which are the inverse of FRF matrices. Method 3 estimates damping using the equation which relates a damping matrix to the difference between the analytical and measured FRFs. The characteristics of the three methods are investigated by applying them to simulated discrete system data and experimental cantilever beam data.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.9
• /
• pp.1806-1813
• /
• 2002

It is well known that a crack has an important effect on the dynamic behavior of a structure. This effect depends mainly on the location and depth of the crack. To identify the location and depth of a crack in a structure, a classical optimization technique was adopted by previous researchers. That technique overcame the difficulty of finding the intersection point of the superposed contours that correspond to the eigenfrequency caused by the crack presence. However, it is hard to select a trial solution initially for optimization because the defined objective function is heavily multimodal. A method is presented in this paper, which uses continuous evolutionary algorithms(CEAs). CEAs are effective for solving inverse problems and implemented on PC clusters to shorten calculation time. With finite element model of the structure to calculate eigenfrequencies, it is possible to formulate the inverse problem in optimization format. CEAs are used to identify the crack location and depth minimizing the difference from the measured frequencies. We have tried this new idea on a simple beam structure and the results are promising with high parallel efficiency over about 94%.

• Journal of Institute of Control, Robotics and Systems
• /
• v.9 no.1
• /
• pp.50-56
• /
• 2003

In electrical impedance tomography(EIT), various image reconstruction algorithms have been used in order to compute the internal resistivity distribution of the unknown object with its electric potential data at the boundary. Mathematically the EIT image reconstruction algorithm is a nonlinear ill-posed inverse problem. This paper presents a new combined method based on genetic algorithm(GA) and modified Newton-Raphson(mNR) algorithm via two-step approach for the solution of the static EIT inverse problem. In the first step, each mesh is classified into three mesh groups: target, background, and temporary groups. The mNR algorithm can be used to determine the region of group. In the second step, the values of these resistivities are determined using genetic algorithm. Computer simulations with the 32 channels synthetic data show that the spatial resolution of reconstructed images by the proposed scheme is improved compared to that of the mNR algorithm at the expense of increased computational burden.