• Journal of applied mathematics & informatics
• /
• v.15 no.1_2
• /
• pp.29-51
• /
• 2004

In this paper, we propose a new mixed finite element method, called the characteristics-mixed method, for approximating the solution to Burgers' equation. This method is based upon a space-time variational form of Burgers' equation. The hyperbolic part of the equation is approximated along the characteristics in time and the diffusion part is approximated by a mixed finite element method of lowest order. The scheme is locally conservative since fluid is transported along the approximate characteristics on the discrete level and the test function can be piecewise constant. Our analysis show the new method approximate the scalar unknown and the vector flux optimally and simultaneously. We also show this scheme has much smaller time-truncation errors than those of standard methods. Numerical example is presented to show that the new scheme is easily implemented, shocks and boundary layers are handled with almost no oscillations. One of the contributions of the paper is to show how the optimal error estimates in $L^2(\Omega)$ are obtained which are much more difficult than in the standard finite element methods. These results seem to be new in the literature of finite element methods.

• International Journal of Highway Engineering
• /
• v.14 no.5
• /
• pp.1-9
• /
• 2012

PURPOSES: A viscoelastic axisymmetric finite element analysis code has been developed for stress analysis of asphalt pavement structures. METHODS: Generalized Maxwell Model (GMM) and 4-node isoparametric element were employed for finite element formulation. The code was developed using $C^{+}^{+}$ computer program language and named as KICTPAVE. For the verification of the developed code, a structural model of a pavement system was constructed. The structural model was composed of three layers: asphalt layer, crushed stone layer, and soil subgrade. Two types of analysis were considered for the verification: (1)elastic static analysis, (2)viscoelastic time-dependent analysis. For the elastic static analysis, linear elastic material model was assigned to all the layers, and a static load was applied to the structural model. For the viscoelastic time-dependent analysis, GMM and linear elastic material model were assigned to the asphalt layer and all the other layers respectively, and a cyclic loading condition was applied to the structural model. RESULTS: The stresses and deformations from KICTPAVE were compared with those from ABAQUS. The analysis results obtained from the two codes showed good agreement in time-dependent response of the element under the loading area as well as the surface deformation of asphalt layer, and horizontal and vertical stresses along the axisymmetric axis. CONCLUSIONS: The validity of KICTPAVE was confirmed by showing the agreement of the analysis results from the two codes.

• Transactions of the Korean Society of Pressure Vessels and Piping
• /
• v.14 no.2
• /
• pp.69-76
• /
• 2018

The paper presents preliminary investigation results for the effect of the baseline correction in the acceleration excitation method on finite element seismic analysis results (such as accumulated equivalent plastic strain, equivalent plastic strain considering cyclic plasticity, von Mises effective stress, etc) of nuclear safety Class I components. For investigation, finite element elastic-plastic time-history seismic analysis is performed for a surge line including a pressurizer lower head, a pressurizer surge nozzle, a surge piping, and a hot leg surge nozzle using the Chaboche hardening model. Analysis is performed for various seismic loading methods such as acceleration excitation methods with and without the baseline correction, and a displacement excitation method. Comparing finite element analysis results, the effect of the baseline correction is investigated. As a result of the investigation, it is identified that finite element analysis results using the three methods do not show significant difference.

• Structural Engineering and Mechanics
• /
• v.20 no.2
• /
• pp.241-257
• /
• 2005

A numerical approach combining the finite element method with two different stability criteria namely the Budiansky and the phase-plane buckling criteria is used to study the dynamic buckling phenomena of plate and shell structures subjected to sudden applied loading. In the finite element analysis an explicit time integration scheme is used and the two criteria are implemented in the Finite Element analysis. The dynamic responses of the plate and shell structures have been investigated for different values of the plate and shell imperfection factors. The results indicate that the dynamic buckling time, which is normally considered in predicting elasto-plastic buckling behavior, should be taken into consideration with the buckling criteria for elastic buckling analysis of plate and shell structures. By selecting proper control variables and incorporating them with two dynamic buckling criteria, the unique dynamic buckling load can be obtained and the problems of ambiguity and contradiction of dynamic buckling load of plate and shell structure can be resolved.

• Computational Structural Engineering
• /
• v.6 no.2
• /
• pp.103-114
• /
• 1993

An advanced time-domain finite element formulation is presented for the displacement and stress analysis of isotropic, linear viscoelastic problems of a hereditary-type constitutive law. The semidiscrete finite element method with linear time-stepping scheme and an elastic-viscoelastic correspondence principle are used in the theoretical development. An efficient treatment of hereditary integral is introduced to improve the numerical accuracy, to reduce the computation time, and to avoid the use of large memory storage. Two-dimensional numerical examples of plane strain and plane stress are solved under the assumption on the material property of being elastic in dilatation and like three-element Voigt model in distorsion, and compared with the analytical solutions and the past numerical results to show the versatility and efficiency of the proposed method.

• Proceedings of the IEEK Conference
• /
• 2009.05a
• /
• pp.326-328
• /
• 2009

we an autotuning algorithm for PI controller with unknown plant. The proposed algorithm uses a saturation function and time delay element as a test signal. Since the integral element of PI controller reduces a phase margin and amplitude margin in the closed loop system, the closed loop system could be resulted in unstable with PI controller, To avoid unstable in the closed loop system with PI controller, the proposed algorithm identifies one point information in the 3rd quadrant of Nyquist plot with a time delay element. The proposed method improves an accuracy of one point identified information with one saturation function.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2003.05a
• /
• pp.463-468
• /
• 2003

A finite element formulation to predict the flexural behavior of composite girder is presented in which the early-age properties of concrete are specified including maturing of elastic modulus, creep and shrinkage. The time dependent constitutive relation accounting for the early-age concrete properties is derived in an incremental format by expanding the total form of stress-strain relation by the first order Taylor series with respect to the reference time. The total potential energy of the flexural composite member is minimized to derive the time dependent finite element equilibrium equation. Numerical applications are made for the 3-span double composite steel box girders which is a composite bridge girder filled with concrete at the bottom of the steel box in the negative moment region. The numerical analysis with considering the variation of concrete elastic modulus are performed to investigate the effect of it on the early-age behavior of composite structures. The one dimensional finite element analysis results are compared with the analytical method based on the sectional analysis. Close agreement is observed among the two methods.

• Structural Engineering and Mechanics
• /
• v.59 no.2
• /
• pp.243-259
• /
• 2016

This research develops a finite element code for the transient dynamic analysis of tapered fiber reinforced polymer (FRP) poles with hollow circular cross-section and flexible joints used in power transmission lines. The FRP poles are modeled by tapered beam elements and their flexible joints by a rotational spring. To solve the time equations of transient dynamic analysis, precise time integration method is utilized. In order to verify the utilized formulations, a typical jointed FRP pole under step, triangular and sine pulses is analyzed by the developed finite element code and also ANSYS commercial finite element software for comparison. Thereafter, the effect of joint flexibility on its dynamic behavior is investigated. It is observed that by increasing the joint stiffness, the amplitude of the pole tip deflection history decreases, and the time of occurrence of the maximum deflection is earlier.

• Proceedings of the Korean Institute of IIIuminating and Electrical Installation Engineers Conference
• /
• 2004.05a
• /
• pp.555-558
• /
• 2004

Due to the asymptotic property, deadbeat control can hardly applied to the continuous time system control. But some delay element method can deal such a problem. Except delay element method, well-known digital deadbeat control can br used with the aid of som smoothing elements. In this paper, and order smoothing element is used for the smoothing of the digital deadbeat controller. And this element is argumented to the plant, and so control problem is to control digitally the argumented system. We simulated this control system using Matlab language and finally apply this algorithm to the rotary inverted pendulum system.

• Kyungpook Mathematical Journal
• /
• v.62 no.3
• /
• pp.557-581
• /
• 2022

The anisotropic-diffusion convection equation with exponentially variable coefficients is discussed in this paper. Numerical solutions are found using a combined Laplace transform and boundary element method. The variable coefficients equation is usually used to model problems of functionally graded media. First the variable coefficients equation is transformed to a constant coefficients equation. The constant coefficients equation is then Laplace-transformed so that the time variable vanishes. The Laplace-transformed equation is consequently written as a boundary integral equation which involves a time-free fundamental solution. The boundary integral equation is therefore employed to find numerical solutions using a standard boundary element method. Finally the results obtained are inversely transformed numerically using the Stehfest formula to get solutions in the time variable. The combined Laplace transform and boundary element method are easy to implement and accurate for solving unsteady problems of anisotropic exponentially graded media governed by the diffusion convection equation.