• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1997년도 가을 학술발표회 논문집
• /
• pp.65-72
• /
• 1997

This study is the bending analysis of rectangular plates with 4-sides simply supported by Finite Element Method using substructuring technique. In finite element method, as the more number of finite element, the more dimension of matrix, it is difficult to obtain accuracy solution. In this paper substructuring technique is applied to finite element method in order to reduce the dimension of matrix according to the number of finite element mesh. To validate finite element method using substructuring technique, deflections and moments of rectangular plates by that method is compared with those of references. Considering the symmetry of the plate and load, one fourth of plate is analyzed. Operating time and the error of solutions according to the number of finite element mesh and substructure are compared with each other.

• Structural Engineering and Mechanics
• /
• 제40권6호
• /
• pp.813-827
• /
• 2011

This paper provides a numerical solution for a finite internally cracked plate using hybrid crack element method (HCE). In the formulation, an inclined crack is placed in any place of a rectangular element and the complex variable method is used. The complex potentials are expressed in a series form, and several undetermined coefficients are involved. The complex potentials for the cracked rectangle are first suggested in this paper. Based on a variational principle, the element stiffness matrix can be evaluated. The next steps are same as in the usual finite element method. Several numerical examples with computed stress intensity factor and T-stress are presented.

• Journal of applied mathematics & informatics
• /
• 제7권2호
• /
• pp.493-506
• /
• 2000

Recently the first author and his coworker report a new finite element method for the Poisson equations with homogeneous Dirichlet boundary conditions on a polygonal domain with one re-entrant angle [7], They use the well-known fact that the solution of such problem has a singular representation, deduced a well-posed new variational problem for a regular part of solution and an extraction formula for the so-called stress intensity factor using tow cut-off functions. They use Fredholm alternative an Garding's inequality to establish the well-posedness of the variational problem and finite element approximation, so there is a maximum bound for mesh h theoretically. although the numerical experiments shows the convergence for every reasonable h with reasonable size y imposing a restriction to the support of the extra cut-off function without using Garding's inequality. We also give error analysis with similar results.

• International Journal of Aerospace System Engineering
• /
• 제3권1호
• /
• pp.1-4
• /
• 2016

Functionally graded materials (FGMs) are becoming very popular in various industries due to their effectiveness of the utilization of their constituent elements. However, the modelling of these materials is difficult due to the complex nature of variation of material properties across the thickness. Many shear deformation theories have been developed and employed for the analysis of such functionally graded plates (FGPs). A recently developed inverse hyperbolic shear deformation theory has been successfully employed by Grover et al. [1] for the analysis of laminated composites and sandwich plates. The objective of the study is to obtain finite element solution for the structural analysis of functionally graded plates using inverse hyperbolic shear deformation theory. Finite element analysis facilitates the analysis of complex problems such as functionally graded plates with different boundary conditions and different loadings.

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

Meshless particle method is a numerical technique which does not use the concept of element. This method can easily handle special engineering problems which cause difficulty in the use of finite element method, however it has a drawback that essential boundary condition is not satisfied. In this paper, several studies for satisfying essential boundary conditions and enhancing the accuracy of solutions are discussed. Particular emphasis is placed on a new numerical technique in which finite elements are used on the boundaries to satisfy the essential boundary conditions and meshless particle method is used in the interior domain. For coupling of the two methods interface elements are introduced into the zone between the subdomains using meshless particle method and finite element method. The shape functions and the approximated displacement functions of the interface element are derived with the ramp function based on the shape function of finite elements. The whole numerical procedures are formulated by Galerkin method. Several numerical examples for enhancing the accuracy of solution in the meshless particle method and a new coupling method are presented.

• Earthquakes and Structures
• /
• 제26권2호
• /
• pp.163-174
• /
• 2024

Dynamic equilibrium equations for finite element analysis were derived for the free field one-dimensional shear wave propagation through the horizontally layered soil deposits with the elastic half-space. We expressed Rayleigh's viscous damping consisting of mass and stiffness proportional terms. We considered two cases where damping matrices are defined in the total and relative displacement fields. Two forms of equilibrium equations are presented; one in terms of total motions and the other in terms of relative motions. To evaluate the performance of new equilibrium equations, we conducted two sets of site response analyses and directly compared them with the exact closed-form frequency domain solution. Results show that the base shear force as earthquake load represents the simpler form of equilibrium equation to be used for the finite element method. Conventional finite element procedure using base acceleration as earthquake load predicts exact solution reasonably well even in soil deposits with unrealistically high damping.

• 한국전산구조공학회논문집
• /
• 제17권1호
• /
• pp.1-10
• /
• 2004

최근, 유한요소해석견과의 신뢰도를 향상시키기 위하여 활발하게 연구되고 있는 적응유한요소해석은 반복계산을 통해서 해석결과의 오차가 사용자에 의해 지정된 허용오차와 같아지도록 하는 해석방법이다. 이와 간은 적응유한요소해석은 해석결과의 오차평가와 이에 따른 유한요소의 재구성과정으로 나누어진다. rp방법에서는 절점의 위치를 이동시켜 요소의 크기를 조절하는 r방법과 형상함수찻수를 증가시키는 p방법을 동시에 적용함으로써 적응해석의 유효성을 향상시키고자 하였다. 제안한 rp방법의 특성을 규명하고 적응해석의 유효성을 보이기 위하여 전형적인 이차원 평면문제들을 해석하고 그 결과를 검토하였다.

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

• Structural Engineering and Mechanics
• /
• 제4권4호
• /
• pp.415-424
• /
• 1996

This paper presents a newly suggested calculation method in which an arbitrary quadrilateral element with curved sides is transformed to a normal rectangular one by mapping of coordinates, then the two-dimensional spline is adopted to approach the displacement function of this element. Finally the solution can be obtained by the least-energy principle. Thereby, the application field of Spline Finite Element Method will be extended.

• Korean Journal of Mathematics
• /
• 제13권2호
• /
• pp.139-148
• /
• 2005

In this paper we extend the mixed finite element method and its $L_2$-error estimate for postprocessed solutions by using Crank-Nicolson time-discretization method. Global $O(h^2+({\Delta}t)^2)$-superconvergence for the lowest order Raviart-Thomas element ($Q_0-Q_{1,0}{\times}Q_{0,1}$) are obtained. Numerical examples are presented to confirm superconvergence phenomena.