• Title/Summary/Keyword: Galerkin Finite Element Method

Search Result 247, Processing Time 0.023 seconds

Elasticity solution of multi-layered shallow cylindrical panels subjected to dynamic loading

  • Shakeri, M.;Eslami, M.R.;Alibiglu, A.
    • Steel and Composite Structures
    • /
    • v.2 no.3
    • /
    • pp.195-208
    • /
    • 2002
  • Elasticity solutions to the boundary-value problems of dynamic response under transverse asymmetric load of cross-ply shallow cylindrical panels are presented. The shell panel is simply supported along all four sides and has finite length. The highly coupled partial differential equations are reduced to ordinary differential equations with constant coefficients by means of trigonometric function expansion in the circumferential and axial directions. The resulting ordinary differential equations are solved by Galerkin finite element method. Numerical examples are presented for two (0/90 deg.) and three (0/90/0 deg.) laminations under dynamic loading.

A development of the 3-dimensional stationary drift-diffusion equation solver (3차원 정상상태의 드리프트-확산 방정식의 해석 프로그램 개발)

  • 윤현민;김태한;김대영;김철성
    • Journal of the Korean Institute of Telematics and Electronics D
    • /
    • v.34D no.8
    • /
    • pp.41-51
    • /
    • 1997
  • The device simulator (BANDIS) which can analyze efficiently the electrical characteristics of the semiconductor devices under the three dimensional stationary conditions on the IBM PC was developed. Poisson, electon and hole continuity equations are discretized y te galerkin method using a tetrahedron as af finite element. The frontal solver which has exquisite data structures and advanced input/output functions is dused for the matrix solver which needs the highest cost in the three dimensional device simulation. The discretization method of the continuity equations used in BANDIS are compared with that of the scharfetter-gummel method used in the commercial three-dimensional device. To verify an accuracy and the efficiency of the discretization method, the simulation results of the PN junction diode and the BJT from BANDIS are compared with those of the commercial three-dimensiional device simulator such as DAVINCI. The maximum relative error within 2% and the average number of iterations needed for the convergence is decreased by more than 20%. The total simulation time of the BJT with 25542 nodes is decreased to about 60% compared with that of DAVINCI.

  • PDF

Study of a Mixed Finite Element Model for the Analysis of a Geometrically Nonlinear Plate (기하학적 비선형 판재 해석을 위한 혼합형 FE Model 연구)

  • Kim, Woo-Ram;Choi, Youn-Dae
    • Transactions of the Korean Society of Mechanical Engineers A
    • /
    • v.34 no.10
    • /
    • pp.1427-1435
    • /
    • 2010
  • A mixed finite element model was developed using the classical plate theory to analyze the nonlinear bending of a plate. The appropriate weight functions for the constraints integrated over the domain were determined by the Lagrange multiplier method by using the principle of minimum virtual energy; which provides the constitutive relations between force-like variables and strains. All of detail terms of element wise coefficient matrices and associate tangent matrices to be used in the Newton iterative method are presented. Then, the linear solutions of the current model and those of the traditional displacement model under the SS (simple support) boundary conditions were compared with the existing analytical solution. The post-processed images of the nonlinear results of the force-like variables are presented to show the continuity of the solutions at the joint of the element boundaries. Finally, the converged nonlinear finite element solutions of the current model are compared with those of existing traditional displacement model.

A Lubrication Analysis of Gas Mechanical Face Seals using a High-Order Shape Function (고차 형상함수를 이용한 가스 미케니컬 페이스 시일의 윤활해석)

  • 이안성;양재훈;최동훈
    • Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
    • /
    • 2001.11a
    • /
    • pp.204-211
    • /
    • 2001
  • For the treatment of high compressibility number in the Reynolds equation, a new class of exponential high-order shape functions has been recently introduced in the literatures. In this paper a FE lubrication analysis method of high speed gas mechanical face seals is developed, implementing these shape functions. Their validity and usefulness are presented using 1-D gas bearing models. And a validation of developed 2-D analysis code is shown with a gas flat and spiral groove face seal models.

  • PDF

Water Quality Prediction Model in a Lake by Finite Element Method;Application to Sapkyo Lake (유한요소법에 의한 호소의 수질예측모형;삽교 담수호에 적용하여)

  • Ryu, Byong-Ro;Ahn, Sang-Jin
    • Korean Journal of Environmental Agriculture
    • /
    • v.8 no.1
    • /
    • pp.37-46
    • /
    • 1989
  • A 2-dimensional pollutant transport phenomenon in shallow reservoirs was analyzed by using a finite element method. The Galerkin's weighted residual method, based on linear interpolation, was used and a triangle was adopted as an element. The two dimensional Stock's equation and the advection-diffusion equation integrated over depth were used as governing equations. Also the Newton-Raphson method was introduced to solve the non-linear terms of the equation. The results calculated by the model are in good agreement with the analytical solution for a simplified channel where a known solution is avaiable. An actual application of the model is attempted for Sapkyo Lake with a consideration of the influx of the Sapkyo Chun, the Muhan Chun and kogkyo Chun. Further refined research is needed to evaluate the water quality in the other reservoirs.

  • PDF

Analysis of Anisotropic Turbulent Heat Transfer in Nuclear Fuel Bundles (핵연료 집합체내의 비등방성 난류 열전달에 관한 해석적 연구)

  • Kim, Sin;Park, Goon-Cherl
    • Nuclear Engineering and Technology
    • /
    • v.20 no.1
    • /
    • pp.35-46
    • /
    • 1988
  • The prediction of clad surface temperatures is important to the design and the safety anlaysis of nuclear reactor. The accurate prediction requires the detailed knowledge of the flow structure and heat transfer, which is complicate due to anisotropic turbulent phenomena. A two-equation model including anisotropic eddy viscosity model is applied to forecast the velocity distribution. And the temperature field is calculated with uniform wall heat flux. The Galerkin's weighted residual finite element method has been used to calculate the turbulent quantities right up to the wall. The numerical results show good agreement with available data and that turbulence anisotropy strongly affects on the mean flow and thus the temperature field. And Nu-P/D correlation is established for sodium coolant in close-packed equilateral triangular bundle in the P/D range of 1.05 to 1.30.

  • PDF

Numerical Analysis for Improving Passing Flow Rate Quantity abound a Radiator (라디에이터 통과풍량 확보를 위한 수치적 검토)

  • 김은필;강상훈
    • Journal of Advanced Marine Engineering and Technology
    • /
    • v.25 no.2
    • /
    • pp.304-310
    • /
    • 2001
  • This paper describes the analysis of flow field using a projection finite element method. The projection scheme with a pressure correction is presented for the analysis of an incompressible Navier-Stokes flow. The projection scheme is analyzed numerically and applied to the well-known bench marking problems such as lid driven cavity. Finally, the projection scheme is applied to a flow through the automobiles front. In the automobiles cooling system, the flow through its front is very important to a cooling performance. The results show that the flow quantity increases by locating the position of bumper to the further front position of a car. And, the improvement on the suction part below a bumper achieves the more passing flow quantity. The attachment of an air dam increases passing flow quantity causing the pressure rise to the front part and the pressure drop beneath a car.

  • PDF

HIGHER ORDER OF FULLY DISCREATE SOLUTION FOR PARABOLIC PROBLEM IN $L_{\infty}$

  • Lee, H.Y.;Lee, J.R.
    • Journal of applied mathematics & informatics
    • /
    • v.4 no.1
    • /
    • pp.17-30
    • /
    • 1997
  • In this work we approximate the solution of initialboun-dary value problem using a Galerkin-finite element method for the spatial discretization and Implicit Runge-Kutta method for the spatial discretization and implicit Runge-Kutta methods for the time stepping. To deal with the nonlinear term f(x, t, u), we introduce the well-known extrapolation sheme which was used widely to prove the convergence in $L_2$-norm. We present computational results showing that the optimal order of convergence arising under $L_2$-norm will be preserved in $L_{\infty}$-norm.

A Study on the Bending Buckling Behavior of Circular Cylindrical Shells (원통형 쉘의 휨 좌굴 거동에 대한 연구)

  • 정진환;김성도;하지명
    • Proceedings of the Computational Structural Engineering Institute Conference
    • /
    • 1998.10a
    • /
    • pp.208-215
    • /
    • 1998
  • A stability problems of isotropic shells under pure bending is investigated based on the classical shells theory. The governing equations of stability problem presented by Donnell and Love, are developed and the solutions for the cylindrical shells are obtained by using Galerkin method. Bending moment is applied at the ends of the cylindrical shell as a from of distributed load in the shape of sine curve. For the isotropic materials, the result of the general purpose structural analysis program based on the finite element method are compared with the critical moment obtained from the classical shell theories. The critical loads for the cylindrical shells with various geometry can not be evaluated with a simple equation. However, accurate solutions for the stability problems of cylindrical shells can be obtained through the equilibrium equation developed in the study.

  • PDF

Geometrically nonlinear dynamic analysis of FG graphene platelets-reinforced nanocomposite cylinder: MLPG method based on a modified nonlinear micromechanical model

  • Rad, Mohammad Hossein Ghadiri;Shahabian, Farzad;Hosseini, Seyed Mahmoud
    • Steel and Composite Structures
    • /
    • v.35 no.1
    • /
    • pp.77-92
    • /
    • 2020
  • The present paper outlined a procedure for geometrically nonlinear dynamic analysis of functionally graded graphene platelets-reinforced (GPLR-FG) nanocomposite cylinder subjected to mechanical shock loading. The governing equation of motion for large deformation problems is derived using meshless local Petrov-Galerkin (MLPG) method based on total lagrangian approach. In the MLPG method, the radial point interpolation technique is employed to construct the shape functions. A micromechanical model based on the Halpin-Tsai model and rule of mixture is used for formulation the nonlinear functionally graded distribution of GPLs in polymer matrix of composites. Energy dissipation in analyses of the structure responding to dynamic loads is considered using the Rayleigh damping. The Newmark-Newton/Raphson method which is an incremental-iterative approach is implemented to solve the nonlinear dynamic equations. The results of the proposed method for homogenous material are compared with the finite element ones. A very good agreement is achieved between the MLPG and FEM with very fine meshing. In addition, the results have demonstrated that the MLPG method is more effective method compared with the FEM for very large deformation problems due to avoiding mesh distortion issues. Finally, the effect of GPLs distribution on strength, stiffness and dynamic characteristics of the cylinder are discussed in details. The obtained results show that the distribution of GPLs changed the mechanical properties, so a classification of different types and volume fraction exponent is established. Indeed by comparing the obtained results, the best compromise of nanocomposite cylinder is determined in terms of mechanical and dynamic properties for different load patterns. All these applications have shown that the present MLPG method is very effective for geometrically nonlinear analyses of GPLR-FG nanocomposite cylinder because of vanishing mesh distortion issue in large deformation problems. In addition, since in proposed method the distributed nodes are used for discretization the problem domain (rather than the meshing), modeling the functionally graded media yields to more accurate results.