• Structural Engineering and Mechanics
• /
• 제28권5호
• /
• pp.603-633
• /
• 2008

A four-node plate finite element for the analysis of laminated composites which is developed using strain gradient notation is presented. The element is based on a first-order shear deformation theory and on the equivalent lamina assumption. Strains and stresses can be calculated at different points through the thickness of the plate. They are averaged values due to the equivalent lamina assumption. A shear correction factor is used as the transverse shear strain is taken to be constant over the plate thickness while its actual variation is parabolic. Strain gradient notation, which is physically interpretable, allows for the detailed a-priori analysis of the finite element model. The polynomial expansions are inspected and spurious terms responsible for modeling errors are identified in the shear strains polynomial expansions. The element is corrected by simply removing the spurious terms from the shear strains expansions. The element is implemented into a FORTRAN finite element code in two versions; namely, with and without spurious terms. Results are compared to show the effects of the spurious terms on the solutions. It is also shown that a refined mesh composed of corrected elements provides solutions which approximate very well the analytical solutions, validating the procedure.

• Nuclear Engineering and Technology
• /
• 제24권3호
• /
• pp.319-328
• /
• 1992

일차원 구에서 유한요소법의 Galerkin formulation이 일차형태의 단일 에너지 중성자 수송방정식의 적분법에 적용되었다. 구분적으로 1차 혹은 2차인 Lagrange 다항식들이 선형대수 방정식들의 집합을 만들기 위해 적분법에 있는 각의존 중성자속(angular flux)에 대하여 활용되었다. 수치해석이 균질구에서의 임계문제와 비균질구에서의 scalar flux 분포에 대해서 행해졌다. 공간과 각에 대하여 연속적인 유한요소를 사용한 균질구에서의 임계문제에 대한 유한요소법의 결과들은 이론적인 해들자 비교되었다. 비균질 문제에서는 각자 공간에 대하여 불연속 유한요소를 사용하여 구한 scalar flux 분포는 ANISN code에 의한 계산결과와 잘 일치하였다.

• Journal of Welding and Joining
• /
• 제21권6호
• /
• pp.77-83
• /
• 2003

The finite element analysis was performed for the cracks existing in residual stress fields in order to investigate the effects of configuration of residual stress distribution to the fatigue crack opening behaviour. And the variation of stress distributions adjacent to the crack caused by uploading was examined. The finite element model with contact elements for the crack plane and plane stress elements for the base material and the analytical method based on the superposition principle to estimate crack opening behaviour and the stress distribution adjacent to the crack subjected to uploading were used. The results of the analysis showed that crack opening behaviors and variations of stress distribution caused by uploading were changed depending on the configuration of residual stress distribution. When the crack existed in the region of compressive residual stress and the configuration of compressive residual stress distribution were inclined, a partial crack opening just behind of a crack tip occurred during uploading. Based on the above results, it was clarified that the crack opening behaviour in the residual stress field could be predicted accurately by the finite element analysis using these analytical method and model.

• 한국해안해양공학회지
• /
• 제14권1호
• /
• pp.86-93
• /
• 2002

항만부진등 해석에 사용되는 유한요소의 적정 크기를 검토하기 위해 완전개방 직사각형 항만에 대해서 다양한 크기의 요소를 사용한 수치실험을 실시하였다. 수치실험 결과, 공진주기와 증폭비의 허용 오차율을 모두 2%로 하는 경우 파장당 9개, 1%로 하는 경우 12개의 요소가 필요한 것으로 나타났다. 또한, 요소 수의 증가에 따라 공진주기의 오차율은 선형적으로, 증폭비의 오차는 진동하면서 감소하였으며, 요소 수를 9개보다 적게 하는 경우 증폭비의 오차가 상대적으로 크게 증가하였다

• Structural Engineering and Mechanics
• /
• 제29권3호
• /
• pp.289-300
• /
• 2008

The use of metric trial functions to represent the real stress field in what is called the unsymmetric finite element formulation is an effective way to improve predictions from distorted finite elements. This approach works surprisingly well because the use of parametric functions for the test functions satisfies the continuity conditions while the use of metric (Cartesian) shape functions for the trial functions attempts to ensure that the stress representation during finite element computation can retrieve in a best-fit manner, the actual variation of stress in the metric space. However, the issue of how to handle situations where there is locking along with mesh distortion has never been addressed. In this paper, we show that the use of a consistent definition of the constrained strain field in the metric space can ensure a lock-free solution even when there is mesh distortion. The three-noded Timoshenko beam element is used to illustrate the principles. Some significant conclusions are drawn regarding the optimal strategy for finite element modelling where distortion effects and field-consistency requirements have to be reconciled simultaneously.

• Structural Engineering and Mechanics
• /
• 제5권4호
• /
• pp.471-490
• /
• 1997

In this paper, the classical Irons' patch tests which have been generally accepted for the convergence proof of a finite element are performed for Mindlin plate bending elements with a special emphasis on the nonconforming elements. The elements considered are 4-node and 8-node quadrilateral isoparametric elements which have been dominantly used for the analyses of plate bending problems. It was recognized from the patch tests that some nonconforming Mindlin plate elements pass all the cases of patch tests even though nonconforming elements do not preserve conformity. Then, the clues for the Mindlin plate element to pass the Irons' patch tests are investigated. Also, the convergent characteristics of some nonconforming Mindlin plate elements that do not pass the Irons' patch tests are examined by weak patch tests. The convergence tests are performed on the benchmark numerical problems for both nonconforming elements which pass the patch tests and which do not. Some conclusions on the relationship between the patch test and convergence of nonconforming Mindlin plate elements are drawn.

• International Journal of Concrete Structures and Materials
• /
• 제8권2호
• /
• pp.141-155
• /
• 2014

The use of fiber reinforced polymer (FRP) composites in strengthening reinforced concrete beam-column subassemblies has been scrutinised both experimentally and numerically in recent years. While a multitude of numerical models are available, and many match the experimental results reasonably well, there are not many studies that have looked at the efficiency of different finite elements in a comparative way in order to clearly identify the best practice when it comes to modelling FRP for strengthening. The present study aims at investigating this within the context of FRP retrofitted reinforced concrete beam-column subassemblies. Two programs are used side by side; ANSYS and VecTor2. Results of the finite element modeling using these two programs are compared with a recent experimental study. Different failure and yield criteria along with different element types are implemented and a useful technique, which can reduce the number of elements considerably, is successfully employed for modeling planar structures subjected to in-plane loading in ANSYS. Comparison of the results shows that there is good agreement between ANSYS and VecTor2 results in monotonic loading. However, unlike VecTor2 program, implicit version of ANSYS program is not able to properly model the cyclic behavior of the modeled subassemblies. The paper will be useful to those who wish to study FRP strengthening applications numerically as it provides an insight into the choice of the elements and the methods of modeling to achieve desired accuracy and numerical stability, a matter not so clearly explored in the past in any of the published literature.

• 대한기계학회논문집
• /
• 제10권1호
• /
• pp.157-163
• /
• 1986

본 연구에서는 기존 유한요소 모델에 실험으로 구한 동특성을 사용하여 불확실성 변수를 고려한 개선된 모델을 구하고 개선된 모델로 부터의 해석 결과와 실험 결과를 비교하였다. 또한 대상물의 파손에 의한 동특성 변화를 측정하고 측정된 신호와 개선된 모델을 사용하여 대상물에서의 파손 위치를 오차행렬로 도시 할수 있는 방법을 연구하였다. 또 각 부재의 동특성에 영향을 주는 민감도를 도입 하여 파손위치를 더욱 명확히 오차행렬에 의한 파손 부위 색출방법이 매우 효과적임을 예제 구조물에서 보여 주었다.

• Structural Engineering and Mechanics
• /
• 제18권6호
• /
• pp.731-744
• /
• 2004

Modified virtual crack closure integral (MVCCI) technique has become very popular for computation of strain energy release rate (SERR) and stress intensity factor (SIF) for 2-D crack problems. The objective of this paper is to propose a numerical integration procedure for MVCCI so as to generalize the technique and make its application much wider. This new procedure called as numerically integrated MVCCI (NI-MVCCI) will remove the dependence of MVCCI equations on the type of finite element employed in the basic stress analysis. Numerical studies on fracture analysis of 2-D crack (mode I and II) problems have been conducted by employing 4-noded, 8-noded (regular & quarter-point), 9-noded and 12-noded finite elements. For non-singular (regular) elements at crack tip, NI-MVCCI technique generates the same results as MVCCI, but the advantage for higher order regular and singular elements is that complex equations for MVCCI need not be derived. Gauss numerical integration rule to be employed for 8-noded singular (quarter-point) element for accurate computation of SERR and SIF has been recommended based on the numerical studies.

• Coupled systems mechanics
• /
• 제3권3호
• /
• pp.267-289
• /
• 2014

A proper physical modeling of infilled building frame-foundation beam-soil mass interaction system is needed to predict more realistic and accurate structural behavior under static vertical loading. This is achieved via finite element method considering the superstructure, foundation and soil mass as a single integral compatible structural unit. The physical modelling is achieved via use of finite element method, which requires the use of variety of isoparametric elements with different degrees of freedom. The unbounded domain of the soil mass has been discretized with coupled finite-infinite elements to achieve computational economy. The nonlinearity of soil mass plays an important role in the redistribution of forces in the superstructure. The nonlinear behaviour of the soil mass is modeled using hyperbolic model. The incremental-iterative nonlinear solution algorithm has been adopted for carrying out the nonlinear elastic interaction analysis of a two-bay two-storey infilled building frame. The frame and the infill have been considered to behave in linear elastic manner, whereas the subsoil in nonlinear elastic manner. In this paper, the computational methodology adopted for nonlinear soil-structure interaction analysis of infilled frame-foundation-soil system has been presented.