• 한국지반공학회지:지반
• /
• 제14권4호
• /
• pp.129-140
• /
• 1998

탄소성 유한 요소 해석은 기지의 변형률 증분에 대한 응력 적분을 필요로 하며, 탄소성 구성 모델의 경우 특별한 경우를 제외하고는 해석적인 응력 적분이 불가능하고 수치적인 방법을 필요로 한다. 이때 응력 수치 적분의 정확도가 비선형 유한요소해의 전체적인 정확도에 상당히 큰 영향을 미치게 된다. 본 연구에서는 탄소성 구성 관계의 응력 적분을 위하여 외연적 방법중의 하나로서 Sloan이 제안한 단계분할 절차를 보완하여 안정적이고 정착한 응력 수치 적분법을 제시하고자 한다. 수정 오일러 절차에 따른 오차 조절의 기본 개념은 그대로 사용하고 오차를 평가하는 기준에 응력 수준이 영향을 미치는 단점을 보완하여 응력 수준에 관계없는 안정적이고 정확한 수치 적분법을 제시하였으며, 그 결과의 신뢰성을 삼축시험모사를 이용하여 검증하였다.

• Structural Engineering and Mechanics
• /
• 제55권4호
• /
• pp.885-900
• /
• 2015

This paper introduces Romberg-Richardson's method as one of the numerical integration tools for computation of stress intensity factor in a pre-cracked specimen subjected to a complex stress field across the crack faces. Also, the computation of stress intensity factor for various stress fields using existing three methods: average stress over interval method, piecewise linear stress method, piecewise quadratic method are modified by using Richardson extrapolation method. The direct integration method is used as reference for constant and linear stress distribution across the crack faces while Gauss-Chebyshev method is used as reference for nonlinear distribution of stress across the crack faces in order to obtain the stress intensity factor. It is found that modified methods (average stress over intervals-Richardson method, piecewise linear stress-Richardson method, piecewise quadratic-Richardson method) yield more accurate results after a few numbers of iterations than those obtained using these methods in their original form. Romberg-Richardson's method is proven to be more efficient and accurate than Gauss-Chebyshev method for complex stress field.

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

In this paper the efficient finite element for stress analysis of plane stress/strain problems is proposed. This element is achieved by adding the bubble-mode function to 8-node element. The stiffness matrix of the element is calculated by using modified numerical integration order to avoid spurious zero energy mode. In order to demonstrate the performance of this element numerical tests for various verification problems are carried out. The results of numerical tests show accuracy and reliability of the element presented in this paper.

• International Journal of Railway
• /
• 제3권4호
• /
• pp.117-122
• /
• 2010

In transportation geotechnical engineering, stress-strain behavior of earth structures has been analyzed by numerical simulations with the implemented plasticity constitutive model. It is a fact that many advanced plasticity constitutive models on predicting the mechanical behavior of soils have been developed as well as experimental research works for geotechnical applications in the past decades. In this study, recently developed, a unified constitutive model for both clay and sand, which is referred to as CASM (clay and sand model), was compared with a classical constitutive model, Cam-Clay model. Moreover, integration methods of stress-strain rate equations using CASM were presented for simulation of undrained and drained triaxial compression tests. As a conclusion, it was observed that semi-implicit integration method has more improved accuracy of capturing strain rate response to applied stress than explicit integration by the multiple correction and iteration.

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

• Structural Engineering and Mechanics
• /
• 제85권1호
• /
• pp.89-104
• /
• 2023

This work presents a proposal for employing reduced numerical integration in the formulation of the 4-node quadrilateral solid finite element. The use of these low-order integration rules leads to numerical instabilities such as those producing the hourglass effect. The proposed procedure allows evaluating a given constitutive model only in one integration point, achieving an attractive computational cost reduction and, also, successfully controls the hourglass effect. A validation of the proposal is included and discussed throughout the paper. To show the efficiency of the proposal, several application examples of masonry structures are studied and discussed. To represent the non-linear mechanical behaviour of masonry a plastic-damage model is implemented within the application of this sub-integration scheme. Also, in order to have a full and computationally efficient strategy to determine the behaviour of masonry structures, involving its evolution to collapse, a homogenization technique with a macro-modeling approach is used. The methodology discussed throughout this paper demonstrates a substantial computational cost reduction and an improved approximation of the non-linear problem evidenced by a reduction of up to 85% of the computational time for some cases.

• Structural Engineering and Mechanics
• /
• 제75권3호
• /
• pp.357-367
• /
• 2020

The objective of this article is investigation of dynamic response of thick multilayer functionally graded (FG) beam under generalized dynamic forces. The plane stress problem is exploited to describe the constitutive equation of thick FG beam to get realistic and accurate response. Applied dynamic forces are assumed to be sinusoidal harmonic, sinusoidal pulse or triangle in time domain and point load. Equations of motion of deep FG beam are derived based on the Hamilton principle from kinematic relations and constitutive equations of plane stress problem. The numerical finite element procedure is adopted to discretize the space domain of structure and transform partial differential equations of motion to ordinary differential equations in time domain. Numerical time integration method is used to solve the system of equations in time domain and find the time responses. Numerical parametric studies are performed to illustrate effects of force type, graduation parameter, geometrical and stacking sequence of layers on the time response of deep multilayer FG beams.

• 대한기계학회논문집
• /
• 제5권3호
• /
• pp.217-222
• /
• 1981

The stress analysis of two-orthotropic layer, adhesively bonded structures is considered. An orthotropic plate has a through-crack of finite length and is adhesively bounded by a sound orthotropic plate. The problem is resuced to a pair of Fredholm integral equations ofthe second kind. Using a numerical integration scheme to evaluate the intgrals, The integral equations are reduced to a system of algebraic equations. By solving these equations some numerical results for stress intensity factors are presented for various crack lengths.

• Structural Engineering and Mechanics
• /
• 제66권5호
• /
• pp.649-663
• /
• 2018

The scaled boundary finite element method is extended to evaluate the dynamic stress intensity factors and T-stress with a numerical procedure based on the improved continued-fraction. The improved continued-fraction approach for the dynamic stiffness matrix is introduced to represent the inertial effect at high frequencies, which leads to numerically better conditioned matrices. After separating the singular stress term from other high order terms, the internal displacements can be obtained by numerical integration and no mesh refinement is needed around the crack tip. The condition numbers of coefficient matrix of the improved method are much smaller than that of the original method, which shows that the improved algorithm can obtain well-conditioned coefficient matrices, and the efficiency of the solution process and its stability can be significantly improved. Several numerical examples are presented to demonstrate the increased robustness and efficiency of the proposed method in both homogeneous and bimaterial crack problems.

• 대한기계학회논문집A
• /
• 제20권2호
• /
• pp.450-460
• /
• 1996

To investigate the effect of stress wave propagation for crack tip, impact responses of two-dimensional plates with oblique cracks are investigated by a numerical method. In the numerical analysis, the finite element method is used in space domain discretization and the Newmark constant acceleration algorithm is used in time integration. According to the numerical results from the impact response analysis. it is found that the stress fields are bisected at the crack surface and the parts of stress intensity are moved along the crack face. The crack tip stress fields are yaried rapidly. The magnitude of crack tip stress fields are converted to dynamic stress intensity factor. Dynamic sress intensity factor appears when the stress wave has reached at the crack tip and the aspect of change of dynamic stress intensity factor is shown to be the same as the part of the flow of stress intensity.