• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.2
• /
• pp.445-454
• /
• 1991

In this thesis, Mindlin plate element with nine nodes and three degrees-of-freedom at each node is formulated and is employed in eigen-analysis of a rectangular plates in order to alleviate locking phenomenon of eigenvalues. Eigenvalues and their modes may be locked if conventional $C_{0}$-isoparametric element is used. In order to reduce stiffness locking phenomenon, two methods (1, the general reduced and selective integration, 2, the new element that use of modified shape function) are studied. Additionally in order to reduce the error due to mass matrix, two mass matrixes (1, Gauss-Legendre mass matrix, 2, Gauss-Lobatto mass matrix) are considered. The results of eigen-analysis for two models (the square plate with all edges simply-supported and all edges built-in), computed by two methods for stiffness matrix and by two mass matrixes are compared with theoretical solutions and conventional numerical solutions. These comparisons show that the performance of the two methods with Gauss-Lobatto mass matrix is better than that of the conventional plate element. But, by considering the spurious rigid body motions, the element which employs modified shape function with full integration and Gauss-Lobatto mass matrix can elevate the accuracy and convergence of numerical solutions.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.3007-3019
• /
• 1993

For the same configuration, the stiffness of 6-node two dimensional isoparametric is stiffer than that of 8-node two dimensional isoparametric element. This phenomenon may be called 'Relative Stiffness Stiffening Phenomenon.' In this paper, the relative stiffness stiffening phenomenon was studied, and could be corrected by modifying the position of Gauss integral points used in the numerical integration of the stiffness matrix. For the same deformation (bending) energy of 6-node and 8-node two dimensional isoparametric elements, Gauss integral points of 6-node element have to move closer, in comparison with those of 8-node element, in the case of numerical integration along the thickness direction.

• Kyungpook Mathematical Journal
• /
• v.25 no.1
• /
• pp.83-92
• /
• 1985

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

• International Journal of High-Rise Buildings
• /
• v.3 no.4
• /
• pp.243-254
• /
• 2014

In this study, a demolition analysis code using the adaptively shifted integration (ASI)-Gauss technique, which describes structural member fracture by shifting the numerical integration point to an appropriate position and simultaneously releasing the sectional forces in the element, is developed. The code was verified and validated by comparing the predicted results with those of several experiments. A demolition planning tool utilizing the concept of a key element index, which explicitly indicates the contribution of each structural column to the vertical load capacity of the structure, is also develped. Two methods of selecting specific columns to efficiently demolish the whole structure are demonstrated: selecting the columns from the largest index value and from the smallest index value. The demolition results are confirmed numerically by conducting collapse analyses using the ASI-Gauss technique. The numerical results suggest that to achieve a successful demolition, a group of columns with the largest key element index values should be selected when explosives are ignited in a simultaneous blast, whereas those with the smallest should be selected when explosives are ignited in a sequence, with a final blast set on a column with large index value.

• Coupled systems mechanics
• /
• v.5 no.1
• /
• pp.59-86
• /
• 2016

This paper deals with the free vibration analysis of a dynamical coupled system: flexible gravity dam- compressible rectangular reservoir. The finite element method is used to compute the natural frequencies and modal shapes of the system. Firstly, the reservoir and subsequently the dam is modeled by classical 8-node elements and the natural frequencies plus modal shapes are calculated. Afterwards, a new 21-node element is introduced and the same procedure is conducted in which an efficient method is employed to carry out the integration operations. Finally, the coupled dam-reservoir system is modeled by solely one 21-node element and the free vibration of dam-reservoir interaction system is investigated. As an important result, it is clearly concluded that the one high-order element treats more precisely than the eight-node elements, since the first one utilizes fifth-degree polynomials to construct the shape functions and the second implements polynomials of degree two.

• Proceedings of the Korean Statistical Society Conference
• /
• 2003.05a
• /
• pp.79-84
• /
• 2003

Likelihood estimation in random-effect models is often complicated because the marginal likelihood involves an analytically intractable integral. Numerical integration such as Gauss-Hermite quadrature is an option, but is generally not recommended when the dimensionality of the integral is high. An alternative is the use of hierarchical likelihood, which avoids such burdensome numerical integration. These two approaches for fitting binary data are compared and the advantages of using the hierarchical likelihood are discussed. Random-effect models for binary outcomes and for bivariate binary-continuous outcomes are considered.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.9 no.1
• /
• pp.18-24
• /
• 1984

Calculation of error probabilities for a coherent phase-shilft keyed communication system operating in a transionospheric scintillation channel is accomplished by means of the Gauss-quadrature integration formula. The channel model used, patterned after Rino's work, is slowly flat fading wherein the envelope of the received signal is modeled as the envelope of correlated Gaussian quadrature random processes. The error probability for the scintillation channel is calculated using actual ionospheric scintillation data for transmission in the UHF region(30-300MHz).

• Communications for Statistical Applications and Methods
• /
• v.26 no.3
• /
• pp.315-323
• /
• 2019

Panel data sets have recently been developed in various areas, and many recent studies have analyzed panel, or longitudinal data sets. Often a dichotomous dependent variable occur in survival analysis, biomedical and epidemiological studies that is analyzed by a generalized linear mixed effects model (GLMM). The most common estimation method for the binary panel data may be the maximum likelihood (ML). Many statistical packages provide ML estimates; however, the estimates are computed from numerically approximated likelihood function. For instance, R packages, pglm (Croissant, 2017) approximate the likelihood function by the Gauss-Hermite quadratures, while Rchoice (Sarrias, Journal of Statistical Software, 74, 1-31, 2016) use a Monte Carlo integration method for the approximation. As a result, it can be observed that different packages give different results because of different numerical computation methods. In this note, we discuss the pros and cons of numerical methods compared with the exact computation method.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.93-98
• /
• 2001

An adaptive Cartesian grid method having the best elements of structured, unstructured, and Cartesian grids is developed to solve the steady two-dimensional Euler equations. The solver is based on a cell-centered finite-volume method with Roe's flux-difference splitting and implicit point Gauss-seidel time integration method. Calculations of several compressible flows are carried out to show the efficiency of the developed computer code. The results were generally in good agreements with existing data in the literature and the developed code has the good ability to capture important feature of the flows.