• Geotechnical Engineering
• /
• v.14 no.4
• /
• pp.129-140
• /
• 1998

Elasto-plastic finite element analysis of geotechnical boundary value problems necessitate the stress integration for the known strain increments. For the elasto-plastic constitutive model, the stress integration is generally achieved by numerical schemes, because analytical integration is impossible for general strain path. In this case, the accuracy of numerical stress integration has an important role on the overall accuracy of nonlinear finite element solution. In this study, the Sloan's substepping method which is one of explicit integration methods has been adopted and iris applicability has been checked. The unstability and inaccuracy of ifs results initiated from initial stress level were revealed. So. a new modified numerical integration method which employs the basic concept of modified Euler scheme for error control is proposed and accuracy and stability of the solutions are confirmed by triaxial test simulation.

• Structural Engineering and Mechanics
• /
• v.55 no.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.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2001.04a
• /
• 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
• /
• v.3 no.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
• /
• v.18 no.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
• /
• v.85 no.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
• /
• v.75 no.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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.5 no.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
• /
• v.66 no.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.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.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.