• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.9
• /
• pp.2887-2893
• /
• 1996

Governing equations infrictional contact problems are introduced using variational inequality formulations which are regularized to overcome the diffculties of non-differentiability of the friction functional. Also finite element approximations and a priori error estimations are derived based on those formulations. Numerical simulations are performed illustrating the theoretical results.

• Water Engineering Research
• /
• v.2 no.3
• /
• pp.171-178
• /
• 2001

There are three methods for analyzing flow and pressure distribution in looped water distribution networks (the loop method, the node method, the element method) taking into consideration hydraulic parameters chosen as unknown. For all these methods the non-linear system of equations can be solved by iterative procedures. The paper presents a different approach to this problem by using the method of variational formulations for hydraulic analysis of water distribution networks. This method has the advantage that it uses a specialized optimization algorithm which minimizes directly an objective multivariable function without constraints, implemented in a computer program. The paper compares developed method to the classic Hardy-Cross method. This shows the good performance of the new method.

• Journal of Ship and Ocean Technology
• /
• v.12 no.2
• /
• pp.1-15
• /
• 2008

The main idea of this paper introduce stochastic structural parameters and random dynamic excitation directly into the dynamic functional variational formulations, and developed the nonlinear dynamic analysis of a stochastic variational principle and the corresponding stochastic finite element method via the weighted residual method and the small parameter perturbation technique. An interpolation method was adopted, which is based on representing the random field in terms of an interpolation rule involving a set of deterministic shape functions. Direct integration Wilson-${\theta}$ Method was adopted to solve finite element equations. Numerical examples are compared with Monte-Carlo simulation method to show that the approaches proposed herein are accurate and effective for the nonlinear dynamic analysis of structures with random parameters.

• Transactions of Materials Processing
• /
• v.12 no.1
• /
• pp.26-33
• /
• 2003

This paper is concerned with the numerical analysis of frictional contact problems in axisymmetric bodies using the rigid-plastic finite element method. A contact finite element method, based on a penalty function, are derived from variational formulations. The contact boundary condition between two deformable bodies is prescribed by the proposed algorithm. The program which can handle frictional contact problem is developed by using pre-existing rigid-plastic finite element code. Some examples used in this paper illustrate the effectiveness of the proposed formulations and algorithms. Efforts focus on the deformation patterns, contact force, and velocity gradient through the various simulations.

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

• Journal of Magnetics
• /
• v.16 no.3
• /
• pp.240-245
• /
• 2011

A new hybrid ON/OFF method is presented for the fast solution of electromagnetic inverse problems in high frequency domains. The proposed method utilizes both topological sensitivity (TS) and material sensitivity (MS) to update material properties in unit design cells. MS provides smooth design space and stable convergence, while TS enables sudden changes of material distribution when MS slows down. This combination of two sensitivities enables a reduction in total computation time. The TS and MS analyses are based on a variational approach and an adjoint variable method (AVM), which permits direct calculation of both sensitivity values from field solutions of the primary and adjoint systems. Investigation of the formulations of TS and MS reveals that they have similar forms, and implementation of the hybrid ON/OFF method that uses both sensitivities can be achieved by one optimization module. The proposed method is applied to dielectric material reconstruction problems, and the results show the feasibility and effectiveness of the method.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.49-55
• /
• 2009

This paper evaluates performances of a recently developed divergence-free finite element method based on Hermite interpolated stream functions. Velocity bases are derived from Hermite interpolated stream functions to form divergence-free basis functions. These velocity basis functions constitute a solenoidal function space, and the simple gradient of the Hermite functions constitute an irrotational function space. The incompressible Navier-Stokes equation is orthogonally decomposed into a solenoidal and an irrotational parts, and the decoupled Navier-Stokes equations are projected onto their corresponding spaces to form proper variational formulations. To access accuracy and convergence of the present algorithm, three test problems are selected. They are lid-driven cavity flow, flow over a backward-facing step and buoyancy-driven flow within a square enclosure. Hermite interpolation functions from cubic to quintic are chosen to run the test problems. Numerical results are shown. In all cases it has shown that the present method has performed well in accuracies and convergences. Moreover, the present method does not require an upwinding or a stabilized term.

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.57-65
• /
• 2012

Various numerical methods in simulation of seismic wave propagation have been developed. Recently an innovative numerical method called as the Spectral Element Method (SEM) has been developed and used in wave propagation in 3-D elastic media. The SEM that easily implements the free surface of topography combines the flexibility of a finite element method with the accuracy of a spectral method. It is generally used a weak formulation of the equation of motion which are solved on a mesh of hexahedral elements based on the Gauss-Lobatto-Legendre integration rule. Variational formulations of velocity-stress motion are newly modified in order to implement ADE-PML (Auxiliary Differential Equation of Perfectly Matched Layer) in wave propagation in 3-D elastic media, because a general weak formulation has a difficulty in adapting CFS (Complex Frequency Shifted) PML (Perfectly Matched Layer). SEM of Velocity-Stress motion having ADE-PML that is very efficient in absorbing waves reflected from finite boundary is verified with simulation of 1-D and 3-D wave propagation.

• Structural Engineering and Mechanics
• /
• v.7 no.1
• /
• pp.111-126
• /
• 1999

The work presented here focuses on the development of suitable discretised formulations, for large-displacement shape and non-shape design sensitivity analysis (DSA), which enable the straightforward incorporation of structural optimisation into established finite element analysis (FEA) codes. For the generalised displacement-based functional the design sensitivity vector has been expressed in terms of displacement sensitivity. The Total Lagrangian formulation is utilised for modelling of large deformation of truss structures. The variational formulation of the sensitivity analysis procedure is discretised by using "pseudo" - finite elements, Results are presented for the sensitivity analysis and optimisation of standard truss structures. For the purposes of this work, the analysis and optimisation procedures outlined below are incorporated into the FEA code ABAQUS.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.5
• /
• pp.9-14
• /
• 2011

The divergence-free finite elements introduced in this paper are derived from Hermite functions, which interpolate stream functions. Velocity bases are derived from the curl of the Hermite functions. These velocity basis functions constitute a solenoidal function space, and the gradient of the Hermite functions constitute an irrotational function space. The incompressible Navier-Stokes equation is orthogonally decomposed into its solenoidal and irrotational parts, and the decoupled Navier-Stokes equations are then projected onto their corresponding spaces to form appropriate variational formulations. The degrees of the Hermite functions we introduce in this paper are bi-cubis, quartic, and quintic. To verify the accuracy and convergence of the present method, three well-known benchmark problems are chosen. These are lid-driven cavity flow, flow over a backward facing step, and buoyancy-driven flow within a square enclosure. The numerical results show good agreement with the previously published results in all cases.