• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.6 no.1
• /
• pp.1-15
• /
• 2002

Advection-dominated transport problems possess difficulties in the design of numerical methods for solving them. Because of the hyperbolic nature of advective transport, many characteristic numerical methods have been developed such as the classical characteristic method, the Eulerian-Lagrangian method, the transport diffusion method, the modified method of characteristics, the operator splitting method, the Eulerian-Lagrangian localized adjoint method, the characteristic mixed method, and the Eulerian-Lagrangian mixed discontinuous method. In this paper relationships among these characteristic methods are examined. In particular, we show that these sometimes diverse methods can be given a unified formulation. This paper focuses on characteristic finite element methods. Similar examination can be presented for characteristic finite difference methods.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.6 s.249
• /
• pp.705-712
• /
• 2006

A finite element method based on the penalized subdomain-interface framework is proposed for fully-coupled, nonlinear thermomechanical analyses with thermal contact anuor radiation boundaries. In the variational formulation, a well-known penalty functional scheme is adopted for connecting subdomains and interfaces that satisfy various continuity requirements. As a logical consequence, the whole domain can be arbitrarily divided into independently-modeled subdomains without considering the conformity of meshes along their interfaces. Since the nonlinearities due to the contact and radiation boundaries can be localized within a few subdomains, the computational efficiency of the present method is greatly increased with appropriate solution algorithms. By solving some numerical problems, these advantageous features are confirmed carefully.

• Journal of the Korean Society for Precision Engineering
• /
• v.15 no.5
• /
• pp.120-127
• /
• 1998

Differential inequalities occurring in problems of obstacle contact problems are recast into variational inequalities and analyzed by finite element methods. A new a posteriori error estimator, which is essential in adaptive finite element method, is introduced to capture the errors in finite element approximations of these variational inequalities. In order to construct a posteriori error estimates, saddle point problems are introduced using Lagrange parameters and upper bounds are provided. The global upper bound is localized by a special mixed formulation, which leads to upper bounds of the element errors. A numerical experiment is performed on an obstacle contact problem to check the effectivity index both in a local and a global sense.

• Structural Engineering and Mechanics
• /
• v.11 no.6
• /
• pp.591-604
• /
• 2001

In this paper, a finite element with embedded displacement discontinuity which eliminates the need for remeshing of elements in the discrete crack approach is applied for the progressive fracture analysis of concrete structures. A finite element formulation is implemented with the extension of the principle of virtual work to a continuum which contains internal displacement discontinuity. By introducing a discontinuous displacement shape function into the finite element formulation, the displacement discontinuity is obtained within an element. By applying either a nonlinear or an idealized linear softening curve representing the fracture process zone (FPZ) of concrete as a constitutive equation to the displacement discontinuity, progressive fracture analysis of concrete structures is performed. In this analysis, localized progressive fracture simultaneous with crack closure in concrete structures under mixed mode loading is simulated by adopting the unloading path in the softening curve. Several examples demonstrate the capability of the analytical technique for the progressive fracture analysis of concrete structures.

• Structural Engineering and Mechanics
• /
• v.20 no.3
• /
• pp.335-346
• /
• 2005

The T-stress of cracks in elastic sheets is solved by using the fractal finite element method (FFEM). The FFEM, which had been developed to determine the stress intensity factors of cracks, is re-applied to evaluate the T-stress which is one of the important fracture parameters. The FFEM combines an exterior finite element model with a localized inner model near the crack tip. The mesh geometry of the latter is self-similar in radial layers around the tip. The higher order Williams series is used to condense the large numbers of nodal displacements at the inner model near the crack tip to a small set of unknown coefficients. Numerical examples revealed that the present approach is simple and accurate for calculating the T-stresses and the stress intensity factors. Some errors of the T-stress solutions shown in the previous literature are identified and the new solutions for the T-stress calculations are presented.

• Structural Engineering and Mechanics
• /
• v.59 no.1
• /
• pp.123-132
• /
• 2016

Many military and commercial aging aircrafts flying beyond their design life may experience severe crack and corrosion damage, and thus lead to catastrophic failures. In this paper, were used in a finite element model to evaluate the effect of corrosion on the adhesive damage in bonded composite repair of aircraft structures. The damage zone theory was implemented in the finite element code in order to achieve this objective. In addition, the effect of the corrosion, on the repair efficiency. Four different patch shapes were chosen to analyze the adhesive damage: rectangular, trapezoidal, circular and elliptical. The modified damage zone theory was implemented in the FE code to evaluate the adhesive damage. The obtained results show that the adhesive damage localized on the level of corrosion and in the sides of patch, and the rectangular patch offers high safety it reduces considerably the risk of the adhesive failure.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2000.04a
• /
• pp.89-92
• /
• 2000

In this paper, an effective procedure is presented for the local recovery of displacements and stresses in multilayered composite panels, which incorporate the local refinement using mesh superposition. The mesh superposition method is used to refine the global coarse mesh by superimposing refined mesh to the localized zone of interest without transition zones. The finite element model used is a solid element based on the Hellinger-Reissner variational principle. The a posteriori computation of the through-the-thickness distributions of displacements and stresses is achieved using a predictor-corrector procedure. The procedure utilizes the superconvergent stresses and nodal displacements of the finite element patch. The element patch is generated by locally superimposing a refined local mesh to the coarse global mesh.

• Interaction and multiscale mechanics
• /
• v.2 no.4
• /
• pp.353-374
• /
• 2009

A new seamless multiscale simulation was developed for coupling the continuum model with its molecular dynamics. Kriging-based Finite Element Method (K-FEM) is employed to model the continuum base of the entire domain, while the molecular dynamics (MD) is confined in a localized domain of interest. In the coupling zone, where the MD domain overlaps the continuum model, the overall Hamiltonian is postulated by contributions from the continuum and the molecular overlays, based on a quartic spline scaling parameter. The displacement compatibility in this coupling zone is then enforced by the Lagrange multiplier technique. A multiple-time-step velocity Verlet algorithm is adopted for its time integration. The validation of the present method is reported through numerical tests of one dimensional atomic lattice. The results reveal that at the continuum/MD interface, the commonly reported spurious waves in the literature are effectively eliminated in this study. In addition, the smoothness of the transition from MD to the continuum can be significantly improved by either increasing the size of the coupling zone or expanding the nodal domain of influence associated with K-FEM.

• Structural Engineering and Mechanics
• /
• v.69 no.2
• /
• pp.193-204
• /
• 2019

This paper presents results from experimental and numerical studies on the response of steel-concrete composite box bridge girders under certain localized fire exposure conditions. Two composite box bridge girders, a simply supported girder and a continuous girder respectively, were tested under simultaneous loading and fire exposure. The simply supported girder was exposed to fire over 40% of its span length in the middle zone, and the two-span continuous girder was exposed to fire over 38% of its length of the first span and full length of the second span. A measurement method based on comparative rate of deflection was provided to predict the failure time in the hogging moment zone of continuous composite box bridge girders under certain localized fire exposure condition. Parameters including transverse and longitudinal stiffeners and fire scenarios were introduced to investigate fire resistance of the composite box bridge girders. Test results show that failure of the simply supported girder is governed by the deflection limit state, whereas failure of the continuous girder occurs through bending buckling of the web and bottom slab in the hogging moment zone. Deflection based criterion may not be reliable in evaluating failure of continuous composite box bridge girder under certain fire exposure condition. The fire resistance (failure time) of the continuous girder is higher than that of the simply supported girder. Data from fire tests is successfully utilized to validate a finite element based numerical model for further investigating the response of composite box bridge girders exposed to localized fire. Results from numerical analysis show that fire resistance of composite box bridge girders can be highly influenced by the spacing of longitudinal stiffeners and fire severity. The continuous composite box bridge girder with closer longitudinal stiffeners has better fire resistance than the simply composite box bridge girder. It is concluded that the fire resistance of continuous composite box bridge girders can be significantly enhanced by preventing the hogging moment zone from exposure to fire. Longitudinal stiffeners with closer spacing can enhance fire resistance of composite box bridge girders. The increase of transverse stiffeners has no significant effect on fire resistance of composite box bridge girders.

• KIEE International Transaction on Electrical Machinery and Energy Conversion Systems
• /
• v.11B no.2
• /
• pp.15-21
• /
• 2001

This paper presents an improved characteristic analysis methodology for a 5-phase hybrid stepping motor. The basic approach is based on the use of equivalent magnetic circuit taking into account the localized saturation throughout the hybrid stepping motor. The finite element method(FEM) is used to generate the magnetic circuit parameters for the complex stator and rotor teeth and airgap considering the saturation effects in tooth and poles. In addition, the neural network is used to map a change of parameters and predicts their approximation. Therefore, the proposed method efficiently improves the accuracy of analysis by using the parameter characterizing localized saturation effects and reduces the computational time by using the neural network. An improved circuit model of 5-phase hybrid stepping motor is presented and its application is provided to demonstrate the effectiveness of the proposed method.