• Journal of the Korean Institute of Telematics and Electronics
• /
• v.22 no.5
• /
• pp.83-86
• /
• 1985

An alternative technique (or the numerical solution of Fredholm integral equations of second kind is presented. The approximate solution is obtained by fitting the data in mixed form at knots in the region of the problem. To decrease the error in the numerical solution, cubic B-spline functions which are twice continuously differentiable at knots are employed as basis function. For a given example, the results of this technique are compared with those of Moment method employing pulse functions for basis function and delta functions for test function and found to br in good agreement.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.10a
• /
• pp.591-598
• /
• 2002

The interface element method for non-matching FEM meshes is extended using stabilized nodal integration. Two non-matching meshes are shown to be joined together compatibly, with the aid of the moving least square approximation. Using stabilized nodal integration, the interface element method is able to satisfy the patch test, which guarantees the convergence of the method.

• Magazine of the Korea Concrete Institute
• /
• v.7 no.3
• /
• pp.175-186
• /
• 1995

The determination of the stress intensity factors is investigated by sur-face integral method around the crack tip of the nlass~vc: concrete structure. The surface integral met hod is naturally derived from the standal-ci path integral J. Howevcr. In the J integral method, pressure in the crack-face and body forces can not be considered, while this theory has advantage of ccmsidering many kind of forces, so t.his theory will be useful in investigating more accurate strt:ss states around crack tip. Furthermore. t h~s rrlethod can elerninate unntussary process of using singular elements and fine mesh around crack tip which is used 11; modelling the singularity around crack tip. A computer program for determming $K_I$, $K_{II}$ is tfcvulopcd by applying this theory. $K_I$, $K_{II}$ values usmg X noded isoparametric elements which was proved and variation of the stress intensity factor was investigated by application of darn structures.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.491-499
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed tot the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted tot in the sense of yon Karman hypothesis. The material model is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized lot anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed P-version finite element model is demonstrated through several comparative points of iew in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic tone.

• Journal of Ocean Engineering and Technology
• /
• v.12 no.2 s.28
• /
• pp.1-21
• /
• 1998

In this paper, a non-conforming 3-D 8-node solid element(MQM10) has beets applied to the analyses of static and dynamic bending problems of laminated composite plates The QM10 element exhibits stiffer bending stiffness which is caused by the reduction of degree of freedom from Q11 element. As an effective way to correct the relative stiffness stiffening phenomenon the modification of Gauss sampling points for composite plates is proposed. The quantity of modification is a function of material properties. Also, another two modified equations are obtained, one is modification for stress, and the other is modification of coefficient of shear modulus in free vibration. It is noted that MQM10 element can analyse the static and free vibration problems of various 3-dimensional composite plates composed of unidirectional laminae, woven laminae or braided laminae. The results of MQM10 element are in good agreement with those of 20-node element.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.123-131
• /
• 2005

According to ow previous study, we confirmed That the Petrov-Galerkin natural element method(PG-NEM) completely resolves the numerical integration inaccuracy in the conventional Bubnov-Galerkin natural element method(BG-NEM). This paper is an extension of PG-NEM to two-dimensional geometrically nonlinear problem. For the analysis, a linearized total Lagrangian formulation is approximated with the PS-NEM. At every load step, the grid points ate updated and the shape functions are reproduced from the relocated nodal distribution. This process enables the PG-NEM to provide more accurate and robust approximations. The representative numerical experiments performed by the test Fortran program, and the numerical results confirmed that the PG-NEM effectively and accurately approximates The large deformation problem.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.4 no.3
• /
• pp.53-61
• /
• 1984

A nonlinear formulation of a beam finite element(NB6) on the total Lagrangian mode for the geometrically nonlinear analysis of two-dimensional elastic framed structures is presented. The NB6 beam element has been degenerated from the three-dimensional continuum by introducing the deep beam assumptions and consists of three reference nodes and three relative nodes. The element characteristics are derived by discretizing the beam equations of motion using the Galerkin weighted residual method and are reduced-integrated repeatedly for each loading step by the Newton-Raphson iteration techpique. Several numerical examples are given to demonstrate the accuracy and versatility of the proposed nonlinear NB6 beam element.

• Computational Structural Engineering
• /
• v.7 no.2
• /
• pp.101-109
• /
• 1994

A degenerated cylindrical shell element for modeling glass fiber reinforced plastic pipes is developed and its performance for static structural analysis under internal uniform pressure is evaluated. The element is a nine node degenerated solid shell element with reduced integration technique, addition of nonconforming displacement modes, and assumed strain method to improve convergence of analysis. Several numerical examples are solved and compared with analytical solutions and other F.E.M programs, The results show that the increment of fiber orientation in the GFRP pipes with reference to the longitudinal axis cause less radial displacements and much stiffness in the pipes. This is reasonable since the internal pressure will primarily cause hoop stresses in the ring and 90-angle ply GFRP ring carry these efficiently in pure tension.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.4
• /
• pp.1320-1332
• /
• 1996

In the present study, we consider modified 5-node equivalent solid element which has smallest degree of freedom among 2-dimensional solid elements accounting bending deformation as well as extensional and shear deformations, We shall investigate static and dynamic characteristics of this element, which is very effective in thin beam, thick beam, large displacement problems, beam of variable thickness, and asymmetrically stepped beam, etc., as well as relatively simple problems of beam. The degree of freedom of this element is 10, which is smaller than 18 of 9-node element, 16 of 8-node elemtns, 12 of modified 6-node element and Q6 element. Therefore, this element is expected to broaden the effective range of application of the solid elements in the bending problems further.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.1
• /
• pp.75-82
• /
• 1987

It is the objective of this study to estimate J-integral by numerical analysis, in which J-integral as aparameters in fracture mechanics can be used to evaluate unstable ductile fracture which is a important problem with respect to structural stability when the scope is beyond small scale yielding criteria. For this, 8-node isoparametric singular element as crack tip element of a cracked material was used to solve plastic blunting phenomenon at crack tip, and crack opening was forced to start when J-value was exceeding fracture toughness $J_{IC}$ of the material. And crack propagation behaviour was treated by using crack opening angle. From this study, it was shown that crack opening, stable propagation and unstable opening point of the cracked material found by using J-value obtained from this study were accord with the other study, so think, J-value obtained from this study can be directly used as a parameter in fracture mechanics to deal with the problem of stable propagation of crack and unstable ductile fracture.