• Magazine of the Korean Society of Agricultural Engineers
• /
• v.36 no.1
• /
• pp.54-62
• /
• 1994

The expanding technique of the Stochastic Finite Element Method(SFEM) is proposed in this paper for adapting direct integration methods in stochastical dynamic analysis of structures. Grafting the direct integration methods and the SFEM together, one can deal with nonlinear structures and nonstationary process problems without any restriction. The stochastical central diffrence and stochastic Houbolt methods are introduced to show the expanding technique, and their adaptabilities are discussed. Results computed by the proposed method (the Stochastic Finite Element Method in Dynamics: SFEMD) for two degree-of-free- dom system are compared with those obtained by Monte Carlo Simulation.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1995.04a
• /
• pp.573-577
• /
• 1995

An integration method for the hypersingular kernels, in the boundary integralequations used for the solution of crack-like problems in elasticity, has been developed. To isolate the stronger singularities, the actual boundaries are replaced by the smoothly curved auxiliary boundaries which provide the detoured, non-singular integration paths. The auxiliary boundary can be interpreted as a contracted form of the actual boundaries except for the singular element where the collocating point is located. For an optimal integration path for every singular collocation point, the auxiliary boundary may have different shape depending on the position of the collocation point on the singular element.

• 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.12 no.5
• /
• pp.507-526
• /
• 2001

In standard finite element algorithms, the local stability conditions are not accounted for in the formulation of the tangent stiffness matrix. As a result, the loss of the local stability is not adequately related to the onset of the global instability. The phenomenon typically arises with material-type localizations, such as shear bands and plastic hinges. This paper addresses the problem in the context of the planar, finite-strain, rate-independent, materially non-linear beam theory, although the proposed technology is in principle not limited to beam structures. A weak formulation of Reissner's finite-strain beam theory is first presented, where the pseudocurvature of the deformed axis is the only unknown function. We further derive the local stability conditions for the large deformation case, and suggest various possible combinations of the interpolation and numerical integration schemes that trigger the simultaneous loss of the local and global instabilities of a statically determined beam. For practical applications, we advice on a procedure that uses a special numerical integration rule, where interpolation nodes and integration points are equal in number, but not in locations, except for the point of the local instability, where the interpolation node and the integration point coalesce. Provided that the point of instability is an end-point of the beam-a condition often met in engineering practice-the procedure simplifies substantially; one of such algorithms uses the combination of the Lagrangian interpolation and Lobatto's integration. The present paper uses the Galerkin finite element discretization, but a conceptually similar technology could be extended to other discretization methods.

• Structural Engineering and Mechanics
• /
• v.8 no.2
• /
• pp.207-231
• /
• 1999

A versatile 4-node shell element which is useful for the analysis of arbitrary shell structures is presented. The element is developed by flat shell approach, i.e., by combining a membrane element with a Mindlin plate element. The proposed element has six degrees of freedom per node and permits an easy connection to other types of finite elements. In the plate bending part, an improved Mindlin plate has been established by the combined use of the addition of non-conforming displacement modes (N) and the substitute shear strain fields (S). In the membrane part, the nonconforming displacement modes are also added to the displacement fields to improve the behavior of membrane element with drilling degrees of freedom and the modified numerical integration (M) is used to overcome the membrane locking problem. Thus the element is designated as NMS-4F. The rigid link correction technique is adopted to consider the effect of out-of-plane warping. The shell element proposed herein passes the patch tests, does not show any spurious mechanism and does not produce shear and membrane locking phenomena. It is shown that the element produces reliable solutions even for the distorted meshes through the analysis of benchmark problems.

• Structural Engineering and Mechanics
• /
• v.67 no.4
• /
• pp.317-326
• /
• 2018

The iterative decomposition coupling formulation of the precise integration finite element method (FEM) and the time domain boundary element method (TD-BEM) is presented for elstodynamic problems. In the formulation, the FEM node and the BEM node are not required to be coincident on the common interface between FEM and BEM sub-domains, therefore, the FEM and BEM are independently discretized. The force and displacement converting matrices are used to transfer data between FEM and BEM nodes on the common interface between the FEM and BEM sub-domains, to renew the nodal variables in the process of the iterations for the un-coincident FEM node and BEM node. The iterative coupling formulation for elastodynamics in current paper is of high modeling accuracy, due to the semi-analytical solution incorporated in the precise integration finite element method. The decomposition coupling formulation for elastodynamics is verified by examples of a cantilever bar under a Heaviside-type force and a harmonic load.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.8 no.2
• /
• pp.247-258
• /
• 2004

The composite shell element is developed for the solution of undamped forced vibration problem of composite curved panels. The finite element used in the current study is an 4-node enhanced assumed shell element with six degrees of freedom per node. The composite shell element is free of both shear and membrane locking phenomenon by using the enhanced assumed strain(EAS) method. A modification to the first-order shear deformation shell theory is proposed, which results in parabolic thorough-thickness distribution of the transverse shear strains and stresses. It eliminates the need for shear correction factors in the first order theory. Newmark's direct integration technique is used for carrying out the integration of the equation motion, to obtain the repones history. Parametric studies of curved composite panels are carried out for forced vibration analysis by geometrical shapes and by laminated composite; such as fiber orientation, stacking sequence.

• Structural Engineering and Mechanics
• /
• v.85 no.4
• /
• pp.485-500
• /
• 2023

The non-linear static analysis of reinforced concrete (RC) structures using the three-dimensional (3D) finite element method is a time-consuming and challenging task. Moreover, this type of analysis encounters numerical problems such as the lack of convergence of results in the stages of growth and propagation of cracks in the structure. The time integration analysis along with the mass scaling (MS) technique is usually used to overcome these limitations. Despite the use of this method in the 3D finite element analysis of RC structures, a comprehensive study has not been conducted so far to assess the effects of the MS method on the accuracy of results. This study aims to evaluate the accuracy of the MS method in the non-linear quasi-static finite element analysis of RC structures. To this aim, different types of RC structures were simulated using the finite element approach based on the implicit time integration method and the mass scaling technique. The influences of effective parameters of the MS method (i.e., the allowable values of increase in the mass of the RC structure, the relationship between the duration of the applied load and fundamental vibration period of the RC structure, and the pattern of applied loads) on the accuracy of the simulated results were investigated. The accuracy of numerical simulation results has been evaluated through comparison with existing experimental data. The results of this study show that the achievement of accurate structural responses in the implicit time integration analyses using the MS method involves the appropriate selection of the effective parameters of the MS method.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.5 no.3
• /
• pp.220-228
• /
• 1997

We discuss finite element approximation and use a Mindlin plate element based upon uniformly reduced numerical integration. The finite element selected for use in this work is a four-node, bilinear displacement element based upon the Mindlin theory of plates. Such elements show good accuracy for laminated composite plates when reduced numerical integration is used to evaluate the element marices. This study presents both the experimental and F.E. results for the natural frequencies of CFRPURETHANE-CFRP Composite plate. Good agreement between experimental and calculated frequencies is achived.

• Computational Structural Engineering
• /
• v.10 no.2
• /
• pp.139-148
• /
• 1997

A finite element method is programmed to analyse the nonlinear behavior of axisymmetric structures. The lst order Mindlin shell theory which takes into account the transversal shear deformation is used to formulate a conical two node element with six degrees of freedom. To evade the shear locking phenomenon which arises in Mindlin type element when the effect of shear deformation tends to zero, the reduced integration of one point Gauss Quadrature at the center of element is employed. This method is the Updated Lagrangian formulation which refers the variables to the state of the most recent iteration. The solution is searched by Newton-Raphson iteration method. The tangent matrix of this method is obtained by a finite difference method by perturbating the degrees of freedom with small values. For the moment this program is limited to the analyses of non-linear elastic problems. For structures which could have elastic stability problem, the calculation is controled by displacement.