• Interaction and multiscale mechanics
• /
• v.3 no.3
• /
• pp.235-255
• /
• 2010

This paper investigates the heat equation for domains subjected to an internal source with a sharp spatial gradient. The solution is first approximated using linear finite elements, and sufficiently small time-step sizes to yield stable simulations. The main area of interest is then in the ability to approximate the solution using Generalized Finite Elements, and again explore the time-step limitations required for stable simulations. Both high order elements, as well as elements with special enrichments are used to generate solutions. When compared to linear finite elements, the high order elements deliver better accuracy at a given level of mesh refinement, but do not offer an increase in critical time-step size. When special enrichment functions are used, the solution can be approximated accurately on very coarse meshes, while yielding solutions which are both accurate and computationally efficient. The major conclusion of interest is that the significantly larger element size yields larger allowable time-step sizes while still maintaining stability of the time-stepping algorithm.

• Journal of Theoretical and Applied Mechanics
• /
• v.3 no.1
• /
• pp.66-80
• /
• 2002

The bilinear formulation proposed earlier by Peters and Izadpanah to develop finite elements in time to solve undamped linear systems, Is extended (and found to be readily amenable) to develop time finite elements to obtain transient responses of both linear and nonlinear, and damped and undamped systems. The formulation Is used in the h-, p- and hp-versions. The resulting linear and nonlinear algebraic equations are differentiated to obtain the first- and second-order sensitivities of the transient response with respect to various system parameters. The present developments were tested on a series of linear and nonlinear examples and were found to yield, when compared with results obtained using other methods, excellent results for both the transient response and Its sensitivity to system parameters. Mostly. the results were obtained using the Legendre polynomials as basis functions, though. in some cases other orthogonal polynomials namely. the Hermite. the Chebyshev, and integrated Legendre polynomials were also employed (but to no great advantage). A key advantage of the time finite element method, and the one often overlooked in its past applications, is the ease In which the sensitivity of the transient response with respect to various system parameters can be obtained. The results of sensitivity analysis can be used for approximate schemes for efficient solution of design optimization problems. Also. the results can be applied to gradient-based parameter identification schemes.

• Structural Engineering and Mechanics
• /
• v.5 no.4
• /
• pp.385-398
• /
• 1997

High order finite element have a greater convergence rate than low order finite elements, and in general produce more accurate results. These elements have the disadvantage of being more computationally expensive and often require a longer time to solve the finite element analysis. High order elements have been used in this paper to obtain a new eigenvalue solution with out re-solving the new model. The optimisation of the eigenvalue via the differentiation of the Rayleigh quotient has shown that the additional nodes associated with the higher order elements can be condensed out and solved using the original finite element solution. The higher order elements can then be used to calculate an improved eigenvalue for the finite element analysis.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.1 s.71
• /
• pp.39-46
• /
• 2006

This paper analyzes the continuous Galerkin method for the space-time discretization of wave equation. The method of space-time finite elements enables the simple 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 called a time slab. The weighted residual process is used to formulate a finite element method for a space-time domain. Instability is caused by a too large time step in successive time steps. A stability problem is described and some investigations for chosen types of rectangular space-time finite elements are carried out. Some numerical examples prove the efficiency of the described method under determined limitations.

• Proceedings of the Korean Society for Agricultural Machinery Conference
• /
• 1993.10a
• /
• pp.1222-1229
• /
• 1993

An advanced design environment , which is based on adaptive and knowledge -based finite elements (INTELMESH), has been developed. Unlike other approaches, INTEMMESH incorporates the information about the object geometry as well as the boundary and loading conditions to generate an ${\alpha}$-priori finite element mesh which is more refined around the critical regions of the problem domain. INTEMMESH is designed for planar domains and axisymmetric 3-D structures of elasticity and heat transfer subjected to mechanical and thermal loading . It intelligently identifies the critical regions/points in the problem domain and utilize the new concepts of substructuring and wave propagation to choose the proper mesh size for them. INTEMMESH generates well-shaped triangular elements by applying trangulartion and Laplacian smoothing procedures. The adaptive analysis involves the intial finite elements analyze and an efficient ${\alpha}$-posteriori error analysis involves the initial finite element anal sis and an efficient ${\alpha}$-posteriori error analysis and estimation . Once a problem is defined , the system automatically builds a finite element model and analyzes the problem though automatic iterative process until the error reaches a desired level. It has been shown that the proposed approach which initiates the process with an ${\alpha}$-priori, and near optimum mesh of the object , converges to the desired accuracy in less time and at less cost. Such an advanced design/analysis environment will provide the capability for rapid product development and reducing the design cycle time and cost.

• Computational Structural Engineering
• /
• v.8 no.3
• /
• pp.91-102
• /
• 1995

For the wave propagation problems, the influence of time-dependent dynamic behavior must be accounted in the analysis. In this study, the dynamic analysis method which combines finite elements and boundary elements is developed for the wave propagation problem modelling the infinity of medium through 3-D boundary elements and underground structure through degenerated finite shell elements. Performing dynamic analysis of underground tunnels by the proposed coupling method of boundary and finite elements, it is found that the change of the stiffness of structures has a good effect on the response. It is also found that the consideration of the repeating effect due to moving traffic loads which is difficult with existing 2-D dynamic analysis can be possible with the 3-D analysis in time domain.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.659-665
• /
• 2010

The finite element method(FEM) is proven to be an effective approximate method of structural analysis if proper element types and meshes are chosen, and recently, the method is often applied to solve complex dynamic and nonlinear problems. A properly chosen element type and mesh yields reliable results for dynamic finite element structural analysis. However, dynamic behavior of a structure may include unpredictably large strains in some parts of the structure, and using the initial mesh throughout the duration of a dynamic analysis may include some elements to go through strains beyond the elements' reliable limits. Thus, the finite element mesh for a dynamic analysis must be dynamically adaptive, and considering the rapid process of analysis in real time, the dynamically adaptive finite element mesh generating schemes must be computationally efficient. In this paper, a computationally efficient dynamically adaptive finite element mesh generation scheme for dynamic analyses of structures is described. The concept of representative strain value is used for error estimates and the refinements of meshes use combinations of the h-method(node movement) and the r-method(element division). The shape coefficient for element mesh is used to correct overly distorted elements. The validity of the scheme is shown through a cantilever beam example under a concentrated load with varying values. The example shows reasonable accuracy and efficient computing time. Furthermore, the study shows the potential for the scheme's effective use in complex structural dynamic problems such as those under seismic or erratic wind loads.

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

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.3 s.174
• /
• pp.550-559
• /
• 2000

One of the most important factor to be considered for the analysis of sheet metal forming processes is the tool surface description for arbitrarily- shaped sheet metal parts. In the present study , finite element approach is used to describe the arbitrarily shaped tool surface. In finite element mesh approach, tool surfaces ar, described by finite elements. The finite elements mesh description of the tool surface, which is originally described by CAD data, needs much time and time-consuming graphic operation. The method, however, has been widely used to describe a complex tool surface. In the present study, the contact searching algorithm for the finite element mesh approach is developed based on cell strategy method and sheet surface normal scheme. For the verification purpose, a clover cup drawing, Baden-Baden oilpan problem and a trunk floor drawing were investigated. The computational results based on the finite element approach were compared with the results of available parametric patch approach and experiments.