• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1990.10a
• /
• pp.40-46
• /
• 1990

This paper introduces a geometric modelling system adopted in a newly developed preprocessor for finite element analysis of three dimensional structures. The formulation is characterized by hierarchical construction of structural model which consists of control points, curves, surfaces and solids. Various surface and solid modeling schemes based on blending functions and boundary representation are systematized for finite element mesh generation. The modeling system is integrated with model synthesis and operations which facilitate modelling of complex structures.

• Structural Engineering and Mechanics
• /
• v.69 no.2
• /
• pp.155-162
• /
• 2019

It is very common to find an empirical formulation in an earthquake design code to calculate fundamental vibration period of a structural system. Fundamental vibration period or frequency is a key parameter to provide adequate information pertinent to dynamic characteristics and performance assessment of a structure. This parameter enables to assess seismic demand of a structure. It is possible to find an empirical formulation related to reinforced concrete structures, masonry towers and slender masonry structures. Calculated natural vibration frequencies suggested by empirical formulation in the literatures has not suits in a high accuracy to the case of rest of the historical masonry bridges due to different construction techniques and wide variety of material properties. For the listed reasons, estimation of fundamental frequency gets harder. This paper aims to present an empirical formulation through Mean Square Error study to find ambient vibration frequency of historical masonry bridges by using a non-linear regression model. For this purpose, a series of data collected from literature especially focused on the finite element models of historical masonry bridges modelled in a full scale to get first global natural frequency, unit weight and elasticity modulus of used dominant material based on homogenization approach, length, height and width of the masonry bridge and main span length were considered to predict natural vibration frequency. An empirical formulation is proposed with 81% accuracy. Also, this study draw attention that this accuracy decreases to 35%, if the modulus of elasticity and unit weight are ignored.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.1 no.2
• /
• pp.42-48
• /
• 1993

In this paper, dynamic collapse behavior of the rectangular tube under impact loading is anlayzed using nonlinear finite element method of shell element. In case of shell element formulation using corotational element coordinates system, dynamic collapse behavior is analyzed without initial imperfection, and with initial imperfection. This paper reveals that the collapse of a rectangular tue without initial imperfection is caused by an error of transformation of the corotational coordinates system.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.5
• /
• pp.16-25
• /
• 1996

In the present work a rigid-plastic finite element formulation using dynamic explicit time integration scheme is proposed for numerical analysis of auto-body panel stamping processes. The rigid-plastic finite element method based on membrane elements has long been employed as a useful numerical technique for the analysis of sheet metal forming because of its time effectiveness. A damping scheme is proposed in order to achieve a stable solution procedure in dynamic sheet forming problems. In order to improve the drawbacks of the conventional membrane elements, BEAM(abbreviated from Bending Energy Augmented Membrane) elements are employed. Rotational damping and spring about the drilling direction are introduced to prevent a zero energy mode. The lumping scheme is employed for the diagonal mass matrix and linearizing dynamic formulation. A contact scheme is developed by combining the skew boundary condition and the direct trial-and-error method. Computations are carried out for analysis of complicated auto-body panel stamping processes such as forming of an oilpan, a fuel tank and a front fender. The numerical results of explicit analysis are compared with the implicit results with good agreements and it is shown that the explicit scheme requires much shorter computational time, especially when the problem becomes more complicated. It is thus shown that the proposed dynamic explicit rigid-plastic finite element method enables an effective computation for complicated autobody panel stamping processes.

• Journal of Ocean Engineering and Technology
• /
• v.15 no.2
• /
• pp.72-78
• /
• 2001

The bead is used to provide properly restraining force in the sheet metal forming process. This bead process includes bending and geometrical non-linearity, and affects the state of binderwrap. Therefore, the analysis of bead process is very important to obtain the desired formability. In this paper, the research about the influence of the punch stroke of bead on the draw-bead process was conducted. Results from the analysis will give useful information to the effective tool design of blank forming process. To analyze the bead process, and elasto-plastic finite element formulation is constructed from the equilibrium equation and the considered boundary conditions involved a proper contact condition. The static-explicit finite element method as a numerical method for the analysis was applied to the analysis program code. It was found that this method could solve too much computation time and convergence problem owing to high non-linearity of bead forming process.

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

• Land and Housing Review
• /
• v.1 no.1
• /
• pp.59-65
• /
• 2010

The paper presents a comparison between a semi-implicit time integration linear finite element implementation and fully-implicit nonlinear Newton-Raphson finite element implementation of a triphasic small strain mixture formulation of an elastic partially saturated porous medium. The pore air phase pressure pa is assumed atmospheric, i.e., $p_a$ = 0, although the formulation and implementation are general to handle increase in pore air pressure as a result of loading, if needed. The solid skeleton phase is assumed linear isotropic elastic and partially saturated 'consolidation' in the presence of surface infiltration and traction is simulated. The verification of the implementation against an analytical solution for partially saturated pore water flow (no deformation) and comparison between the two implementations is presented and the important of the porosity-dependent nature of the partially saturated permeability is assessed on comparison with a commercial code for the partially saturated flow with deformation. As a result, the response of partially saturated permeability subjected to the porosity influences on the saturation of a soil, and the different behaviors of the partially saturated soil between staggered and monolithic coupled programs is worth of attention because the negative pore water pressure in the partially saturated soil depends on the difference.

• Coupled systems mechanics
• /
• v.8 no.5
• /
• pp.393-414
• /
• 2019

This paper presents seismic analysis of concrete gravity dams considering soil-structure-fluid interaction. Displacement based plane strain finite element formulation is considered for the dam and foundation domain whereas pressure based finite element formulation is considered for the reservoir domain. A direct coupling method has been adopted to obtain the interaction effects among the dam, foundation and reservoir domain to obtain the dynamic responses of the dam. An efficient absorbing boundary condition has been implemented at the truncation surfaces of the foundation and reservoir domains. A parametric study has been carried out considering each domain separately and collectively based on natural frequencies, crest displacement and stress at the neck level of the dam body. The combined frequency of the entire coupled system is very less than that of the each individual sub-system. The crest displacement and neck level stresses of the dam shows prominent enhancement when coupling effect is taken into consideration. These outcomes suggest that a complete coupled analysis is necessary to obtain the actual responses of the concrete gravity dam. The developed methodology can easily be implemented in finite element code for analyzing the coupled problem to obtain the desired responses of the individual subdomains.

• Journal of KSNVE
• /
• v.10 no.1
• /
• pp.97-106
• /
• 2000

The dynamic behavior of crankshaft-bearing system in scroll compressor has been investigated using the combined methodologies of finite elements and transfer matrices. The finite element formulation is proposed including the field element for a shaft section and the point element at balancer weight locations, bearing locations, etc., whereas the conventional method is used with the elements. The Houbolt method is used to consider the time march for the integration of the system equations. The linear stiffness and damping coefficients are calculated for a finite cylindrical fluid-film bearing by solving the Reynolds equation, using finite difference method. The orbital response of crankshaft supported on the linear bearing model is obtained, considering balancer weights of motor rotor. And, the steady state displacement of crankshaft are compared with a variation in balancer weight. The loci of crankshaft at bearing locations are composed of the synchronous whirl component and the non-synchronous whirl component.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.1
• /
• pp.37-49
• /
• 2000

This paper is for the first essential study on the development of unified finite element formulations for solving problems related to the dynamics/thermoelastics behavior of solids. In the first part of formulations, the finite element method is based on the introduction of a new quantity defined as heat displacement, which allows the heat conduction equations to be written in a form equivalent to the equation of motion, and the equations of coupled thermoelasticity to be written in a unified form. The equations obtained are used to express a variational formulation which, together with the concept of generalized coordinates, yields a set of differential equations with the time as an independent variable. Using the Laplace transform, the resulting finite element equations are described in the transform domain. In the second, the Laplace transform is applied to both the equation of heat conduction derived in the first part and the equations of motions and their corresponding boundary conditions, which is referred to the transformed equation. Selections of interpolation functions dependent on only the space variable and an application of the weighted residual method to the coupled equation result in the necessary finite element matrices in the transformed domain. Finally, to prove the validity of two approaches, a comparison with one finite element equation and the other is made term by term.