• Structural Engineering and Mechanics
• /
• 제40권2호
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Kyungpook Mathematical Journal
• /
• 제47권4호
• /
• pp.529-548
• /
• 2007

The purpose of this article is to find a relation between the finite difference method and the boundary element method, and propose a new approach deriving a discrete approximation formula as like that of the finite difference method for harmonic functions. We develop a discrete approximation formula on a uniform grid based on the boundary integral formulations. We consider three different boundary integral formulations and derive one discrete approximation formula on the uniform grid for the harmonic function. We show that the proposed discrete approximation formula has the same computational molecules with that of the finite difference formula for the Laplace operator ${\nabla}^2$.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1996년도 봄 학술발표회 논문집
• /
• pp.62-69
• /
• 1996

In this study, an equivalent domain integral (EDI) method is presented to estimate the track-till integral parameter, J-value, for two dimensional cracked elastic bodies which may quantify the severity of the crack-tit) stress fields. The conventional J-integral method based on line integral has been converted to equivalent area or domain integrals by using the divergence theorem. It is noted that the EDI method is very attractive because all the quantities necessary for computation of the domain integrals are readily available in a finite element analysis. The details and its implementation are extened to both h-version finite element model with 8-node isoparametric element and p-version finite element model with high order hierarchic element using Legendre type shape fuctions. The variations with respect to the different path of domain integrals from the crack-tip front and the choice of 5-function have been tested by several examples.

• 소성∙가공
• /
• 제21권4호
• /
• pp.252-257
• /
• 2012

Crystallographic texture evolution during forming processes has a significant effect on the anisotropic flow behavior of crystalline material. In this study, a crystal plasticity finite element method (CPFEM), which incorporates the crystal plasticity constitutive law into a three-dimensional finite element method, was used to investigate texture evolution of a face-centered-cubic material - an aluminum alloy. A rate-dependent polycrystalline theory was fully implemented within an in-house program, CAMPform3D. Each integration point in the element was considered to be a polycrystalline aggregate consisting of a large number of grains, and the deformation of each grain in the aggregate was assumed to be the same as the macroscopic deformation of the aggregate. The texture evolution during three different deformation modes - uniaxial tension, uniaxial compression, and plane strain compression - was investigated in terms of pole figures and compared to experimental data available in the literature.

• 소성∙가공
• /
• 제28권4호
• /
• pp.197-202
• /
• 2019

This paper proposes a precise mathematical model for the prediction of the variation of the roll force across a strip in cold rolling. It further describes the deformation characteristics of the strip using a 3-D finite element method. The different features of hot rolling and cold rolling through a 3-D finite element method are shown. The predicted roll force profile and tension profile are verified through comparison with the prediction from a 3-D finite element method.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2002년도 춘계학술대회논문집
• /
• pp.824-829
• /
• 2002

Since many reciprocating and rotating components are attached to jet loom structure. it is exposed to a more vibration and moise problems than the other textile machinery. Thus the design of the jet loom frame is very important to characterize the dynamic response. In this study, a finite element model of jet loom main frame was developed to investigate the dynamic characteristics of jet loom. Two different finite element models of different main frames were constructed and these models were validated by the experimental results. Dynamic characteristics such natural frequencies and mode shapes were in good agreement between the finite element analysis and experimental results within 10% error range. It is expected that the result from this study can be used as the basic information of jet loom dynamic analysis and be extended for further analysis of forced response case.

• 대한전기학회논문지:전기물성ㆍ응용부문C
• /
• 제49권4호
• /
• pp.249-257
• /
• 2000

In this paper, we propose the use of stochastic finite element method, that is popularly employed in mechanical structure analysis, for more practical designing purpose of RF device. The proposed method is formulated based on the vector finite element method cooperated by pertubation analysis. The method utilizes sensitivity analysis algorithm with covariance matrix of the random variables that represent for uncertain physical quantities such as length or various electrical constants to compute the probabilities of the measure of performance of the structure. For this computation one need to know the variance and covariance of the random variables that might be determined by practical experiences. The presenting algorithm has been verified by analyzing several device with different be determined by practical experiences. The presenting algorithm has been verified by analysis several device with different measure of performanes. For the convenience of formulation, two dimensional analysis has been performed to apply it into waveguide with dielectric slab. In the problem the dielectric constant of the dielectric slab is considered as random variable. Another example is matched waveguide and cavity problem. In the problem, the dimension of them are assumed to be as random variables and the expectations and variances of quality factor have been computed.

• 대한기계학회논문집A
• /
• 제24권7호
• /
• pp.1826-1832
• /
• 2000

A modelling method for the bending vibration analysis of rotating cantilever beams employing finite element method is presented in this paper. Different from the conventional modelling method in wh ich only Cartesian deformation variables are used, a non-Cartesian deformation variable is introduced and approximated to derive the equations of motion. Numerical results obtained by using the presented modelling method are compared to those obtained by using other methods in the related literature, and the accuracy of the presented method is verified through the comparison study. The presented modelling method is superior to other previous methods in a sense that several advantages of the previous methods are incorporated into the presented method.

• 한국환경보건학회지
• /
• 제26권3호
• /
• pp.50-54
• /
• 2000

Ventilation rates can be measured directly by a tracer decay method, although little is known of the effects of different sampling intervals on decay rte calculations. This study determined variations in decay rates calculated by three techniques using residential ozone decay data. The calculation techniques were a regression technique, decay techniques using half-life and average-life, and finite difference techniques using two different time intervals. Variation associated with regression technique calculations for residential ozone decay rates based on data from both sample intervals were within 10% (2.81$\pm$1.88 hr-1). However, both half-life and finite difference technique calculations using a shorter-time interval were significantly different from those obtained with the regression technique(p<0.05). Therefore, the use of short sampling intervals in tracer decay may cause significant error in decay rate calculations.

• 대한기계학회논문집B
• /
• 제33권4호
• /
• pp.242-247
• /
• 2009

AILU type preconditioners for a two-dimensional combined P2P1 finite element formulation of the interaction of rigid cylinder with incompressible fluid flow have been devised and tested by solving fluid-structure interaction (FSI) problems. The FSI code simulating the interaction of a rigid cylinder with an unsteady flow is based on P2P1 mixed finite element formulation coupled with combined formulation. Four different preconditioners were devised for the two-dimensional combined P2P1 finite element formulation extending the idea of Nam et al., which was proposed for the preconditioning of a P2P1 mixed finite element formulation of the incompressible Navier-Stokes equations. It was found that PC-III or PC-IV among them perform well with respect to computational memory and convergence rate for some bench-mark problems.