• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1995.04a
• /
• pp.17-27
• /
• 1995

For the post-buckling and elasto-plastic analysis of shell structures, the total Lagrangian formulation is presented based upon the degenerated shell element. Geometrically correct formulation is developed by updating the direction of normal vectors in the iteration process and evaluating the total Green-Lagrange stain corresponding U total displacements. In the calculation of the stiffness matrix, the element formulation takes into account the effect of finite rotation increments by retaining second order rotation terms in the incremental displacement field. The selective or reduced integration scheme using the heterosis element is applied in order to overcome both shear locking phenomena and the zero energy mode. The load/displacement incremental scheme is adopted for geometric non-linear F .E. analysis. Based on such methodology, the computer program is developed and numerical examples to demonstrate the accuracy and the effectiveness of the proposed shell element are presented and compared with references's results.

• Journal of Mechanical Science and Technology
• /
• v.18 no.8
• /
• pp.1418-1427
• /
• 2004

Combined formulation developed for the fluid-particle mixture is introduced to simulate the biocube-fluid mixture flow, which is utilized for sewage disposal. Some tricky boundary conditions are introduced in order to simulate the effect of screen wall and air bubble, which is injected from the bottom of sewage reservoir. It has been shown that a circulated flow pattern, which was observed in experiment, is reproduced from the present numerical simulation. Furthermore, the effect of biocube density on the distribution pattern of biocube is also studied. It has been shown that a biocube whose density is slightly smaller than that of surrounding fluid or neutrally buoyant one are optimal for the uniform distribution of biocube.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.11
• /
• pp.2312-2321
• /
• 2002

The first-order least-squares meshfree method for linear elasticity is presented. The conventional and the compatibility-imposed least-squares formulations are studied on the convergence behavior of the solution and the robustness to integration error. Since the least-squares formulation is a type of mixed formulation and induces positive-definite system matrix, by using shape functions of same order for both primal and dual variables, higher rate of convergence is obtained for dual variables than Galerkin formulation. Numerical examples also show that the presented formulations do not exhibit any volumetric locking for the incompressible materials.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1992.03a
• /
• pp.201-212
• /
• 1992

The plane strain analysis for simulating the stretch/draw forming operation of arbitrarily-shaped tool profiles and arbitrarily draw-in conditions is introduced. An implicit, incremental, updated Lagrangian formulation is employed, introducing a rigid-viscoplastic constitutive equation. Contact and friction are considered through the mesh-normal, which compatibly describes arbitrary tool surfaces and FEM meshe without depending on the explicit spatial derivatives of tool surfaces. The FEM formulation is tested in the sections automotive inner panel and two-side draw-in. Not only the excellent agreement between measured and computed strains in the stretched section is obtained, but also the numerical stability of current formulation is verified in the two-side draw-in section.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.91-98
• /
• 2004

In two dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a stabilized finite element formulation. The stabilization has been peformed by a modified residual method proposed by Illinca et. al.[3]. The stabilization which is necessary to escape from the LBB constraint renders an equal order formulation. In this paper, we increased the order of interpolation whithin an element up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• Journal of Mechanical Science and Technology
• /
• v.16 no.10
• /
• pp.1239-1245
• /
• 2002

A formulation which seeks steady-state equilibrium positions of constrained multibody systems driven by constant generalized speeds is presented in this paper. Since the relative coordinates are employed, constraint equations at cut joints are incorporated into the formulation. To obtain the steady-state equilibrium position of a multibody system, nonlinear equations are derived and solved iteratively. The nonlinear equations consist of the force equilibrium equations and the kinematic constraint equations. To verify the effectiveness of the proposed formulation, two numerical examples are solved and the results are compared with those of a commercial program.

• Earthquakes and Structures
• /
• v.15 no.5
• /
• pp.565-581
• /
• 2018

This paper proposes a new analytical formulation for computing the seismic input at various levels of a structure in terms of floor response spectra. The approach, which neglects the dynamic interaction between primary structure and secondary element, is particularly useful for the seismic assessment of secondary and non-structural elements. The proposed formulation has a robust theoretical basis and it is based on few meaningful dynamic parameters of the main building. The method has been validated in the linear and nonlinear behavior of the main building through results coming from both experimental tests (available in literature) and parametric numerical analyses. The conditions, for which the Floor Spectrum Approach and its simplified assumptions are valid, have been derived in terms of specific interval ratios between the mass of the secondary element and the participant mass of the main structure. Finally, a practice-oriented formulation has been derived, which could be easily implementable also at code level.

• Structural Engineering and Mechanics
• /
• v.52 no.1
• /
• pp.205-220
• /
• 2014

This paper is concerned with the modification of multidisciplinary feasible formulation for MDO problems using the integrated coupled approximate models. A drawback of conventional MDFs is the numerical difficulty in decomposing the design variables and deriving the coupled equations of state. To overcome such a drawback of conventional methods, the coupling in analysis and design is resolved by approximating the state variables in each discipline by the response surface method and by modifying the optimization formulation using the corresponding integrated coupled approximate models. The validity, reliability and effectiveness of the proposed method are illustrated and verified through two optimization problems, a mathematical MDF problem and the multidisciplinary optimum design of suspension unit of wheeled armored vehicle.

• Journal of Advanced Marine Engineering and Technology
• /
• v.32 no.5
• /
• pp.759-767
• /
• 2008

Bending deformation of an ionic polymer actuator(IPA) on applied low electric field across its thickness is dominated by electroosmosis of hydrated ions and self-diffusion of free water molecules. In the study by Popovic et al., two processes are assumed to occur sequentially in the way that fast electroosmosis is followed by self-diffusion and finite element formulation for the basic field equations are proposed. However the motions of hydrated ions and water molecules occur at the same time. In this study, those two processes are considered simultaneously and finite element formulation is conducted for the basic field equations governing electrochemical response of an IPA. Some numerical studies for IPA are carried out in order to show the validity of the present formulation.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2001.05a
• /
• pp.90-95
• /
• 2001

A continuum-based shape design sensitivity analysis (DSA) method is presented for the acoustic radiation from thin body. The normal derivative integral formulation is employed as an analysis formulation and differentiated directly by using material derivative to get the acoustic shape design sensitivity. In the acoustic sensitivity formulation, derivative coefficients of the structural normal velocities on the surface are required as the input. Thus, the shape design sensitivities of structural velocities on the surface with respect to the shape change are also calculated with continuum approach. A simple disk is considered as a numerical example to validate the accuracy and efficiency of the analytical shape design sensitivity equations derived in this research. This research should be very helpful to design an application involving thin body and to change its acoustic characteristics.