• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.3 s.73
• /
• pp.303-312
• /
• 2006

This paper presents the equivalent nodal load for the element stiffness which represents the influence of the stiffness change such as the addition of elements, the deletion of elements, and/or the partial change of element stiffness. The reanalysis of structure using the equivalent load improves the efficiency very much because the inverse of the structural stiffness matrix, which needs a large amount of computation to calculate, is reused in the reanalysis. In this paper, the concept of the equivalent load for the element stiffness is described and some numerical examples are provided to verify it.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.2
• /
• pp.179-187
• /
• 2000

When one considers the property of the axisymmetry, an analysis of an axisymmetric shell subjected to unaxisymmetric loading can be employed to save time and computer memory space. If one considers the Fourier series of the circumference direction of loads and displacements, an axisymmetric tank subjected to a nonaxisymmetric load can be treated as a frame element. Using the Fourier series, the authors derived the stiffness matrix of the cylindrical shell subjected to unaxisymmetric loading by the usual finite element method, and converted the stiffness matrix of a frame element into a transfer matrix by rearranging the stiffness matrix to apply the transfer matrix method. Here the most significant purpose of this paper is to achieve the fewest number of simultaneous equations for analysing an axisymmetric shell subjected to a nonaxisymmetric load. The results of the proposed method of the analysis of the cylindrical shell subjected to a wind load and a water load show no differences when compared to the other methods.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1995.10a
• /
• pp.221-226
• /
• 1995

터빈 익렬, 펌프 익차, 원형 냉각탑, 치차 등과 같이 동일한 형상이 원주 방향으로 반복되어 있는 순환 대칭 구조물의 진동특성을 유한 요소법을 사용하여 해석하는 경우에 전체구조를 모델링하는 대신에 구조물을 동일한 형상의 부분구조로 분할하여 부분구조 한개만을 모델링하고 분할된 경계에서 적절한 경계조건을 부과하여 진동해석을 수행함으로서 컴퓨터 기억용량을 절감시키고 계산시간을 단축할 수 있는 방법이 널리 사용되고 있다. Orris and Petyt[1]는 부분구조의 양쪽 분할 경계면, 즉 연결 경계상에 있는 절점변위의 상관관계를 복소파동전파식을 이용해서 구하여 부분구조의 감소된 복소강성행렬 및 질량행렬을 만들고 실수부와 허수부를 분리하여 유한요소해석을 수행하는 방법을 제안하였다. 유한요소 프로그램 ANSYS[2]에서는 이와 같은 방법을 사용하고 있다. Thomas[3]는 순회 정규모드를 이용하였고, 참고문헌[4]에서는 순회행렬을 이용하였다. 또한 유한요소 프로그램 MSC/NASTRAN[5]에서는 푸리에 급수를 이용하고 유한요소 절점의 위치 및 변위를 원통 좌표계를 표현하여 순환대칭구조물의 유한요소해석을 수행할 수 있도록 되어있다. 본 논문에서는 순환 대칭구조물의 형상의 주기성과 순환성을 고려하여 이산퓨리에 변환을 이용함으로써 순환대칭구조물의 유한요소진동해석을 체계적으로 저용량의 컴퓨터에서 신속하고 정확하게 수행할 수 있는 방법을 제안하고자 한다.

• Computational Structural Engineering
• /
• v.4 no.2
• /
• pp.9-12
• /
• 1991

공학문제에 있어서, 해석적으로 접근할 수 없었던 많은 경우의 문제들이 유한요소법(Finite Element Methods)의 정형화된 모형화 및 해석과정을 통하여 쉽게 접근되어질 수 있었다. 최근 보다 효율적인 요소개발과 컴퓨터 기술의 발달로 유한요소법은 더욱 효과적인 해석 수단이 되어가고 있다. 그러나 지반공학 문제와 같은 무한영역 문제를 유한요소법으로 해석할 경우, 매우 큰 영역을 모형화하기 위하여 많은 수의 요소가 요구되며 이에 따른 자유도(Degree of Freedom) 수의 증가로 많은 계산시간을 요구하게 된다. 본 고는 무한영역 문제를 효과적으로 모형화하기 위하여 연구, 개발되어진 무한요소(Infinite Element)에 대하여 소개하려 한다. 무한요소의 기본개념과 강성행렬의 형성방법을 보인 후, 기초공학 문제를 예로 하여 이의 적용방법을 간략하게 설명하였다.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.2
• /
• pp.132-139
• /
• 1992

This paper presents geometrically nonlinear formulation of shell problems using the three-dimensional curved shell element, which includs large displacements and large rotations. Formulations of the geometrically nonlinear problems can be derived in a variety of ways, but most of them have been obtained by assuming that nodal rotations are small. Hence, the tangent stiffness matrix is derived under the assumptions that rotational increments are infinitesimal and the effect of finite rotational increments have to be considered during the equilibrium iterations. To study the large displacement and large rotation problems, the restrictions are removed and the formulations of the curved shell element including the effect of large rotational increments are developed in this paper. The displacement based finite element method using this improved formulation are applied to the analyses of the geometrically nonlinear behaviors of the single and double curved shells, which are compared with the results by others.

• Journal of Korean Society of Steel Construction
• /
• v.18 no.3
• /
• pp.349-359
• /
• 2006

This paper presents a new block stiffness matrix for the analysis an orthogonal Cartesian coordinate system. The displacement fields are defined using the first order shear deformable beam theory. The longitudinal displacement can be expressed as the sum of the projected plane deformation of the cross-section due to Timoshenko's beam theory and axial warping deformation due to modified Vlasov's thin-waled beam theory. The derived element takes into account flexural shear deformation and torsional warping deformation. Three different types of beam elements, namely, the two-noded, three-noded, and four-noded beam elements, are developed. The quadratic and cubic elements are found to be very efficient for the flexural analysis of laminated composite beams. The versatility and accuracy of the new element are demonstrated by comparing the numerical results available in the literature.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.10
• /
• pp.2640-2652
• /
• 1994

The stiffness of 6-node isotropic element is stiffer than that of 8-node isotropic element of same configuration. This phenomenon was called 'Relative Stiffness Stiffening Phenomenon'. In this paper, an equation of sampling point modification which correct this phenomenon was derived for the composite plate, as well as an equation for an isotropic plate. The relative stiffness stiffening phenomena of an isotropic plate element could be corrected by modifying Gauss sampling points in the numerical integration of stiffness matrix. This technique could also be successfully applied to the static analyses of composite plate modeled by the 3-dimensional 16-node elements. We predicted theoretical errors of stiffness versus the number of layers that result from the reduction of numerical integration order. These errors coincide very well with the actual errors of stiffness. Therefore, we can choose full integration of reduced integration based upon the permissible error criterion and the number of layers by using the thoretically predicted error.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.2
• /
• pp.113-125
• /
• 2007

Improved numerical method to obtain the exact stiffness matrices is newly proposed to perform the spatially coupled elastic and stability analyses of non-symmetric and open/closed thin-walled beam on elastic foundation. This method overcomes drawbacks of the previous method to evaluate the exact stiffness matrix for the spatially coupled stability analysis of thin-walled beam-column This numerical technique is accomplished via a generalized eigenproblem associated with 14 displacement parameters by transforming equilibrium equations to a set of first order simultaneous ordinary differential equations. Next polynomial expressions as trial solutions are assumed for displacement parameters corresponding to zero eigenvalues and the eigenmodes containing undetermined parameters equal to the number of zero eigenvalues are determined by invoking the identity condition. And then the exact displacement functions are constructed by combining eigensolutions and polynomial solutions corresponding to non-zero and zero eigenvalues, respectively. Consequently an exact stiffness matrix is evaluated by applying the member force-deformation relationships to these displacement functions. In order to illustrate the accuracy and the practical usefulness of this study, the numerical solutions are compared with results obtained from the thin-walled beam and shell elements.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.1 s.71
• /
• pp.63-71
• /
• 2006

This study presents the effective stiffness-based optimal technique to control Quantitatively lateral drift for shear wall-frame structure system using sensitivity analysis. To this end, the element stiffness matrices are constituted to solve the compatibility problem of displacement degree of freedom between the frame and shear wall. Also, lateral drift constraint to introduce the approximation concept that can preserve the generality of the mathematical programming and can effectively solve the large scaled problems is established. And, the section property relationships for shear wall and frame members are considered in order to reduce the number of design variables and differentiate easily the stiffness matrices. Specifically, constant-shape assumption which is uniformly varying in size during optimal process is applied in frame structure. The thickness or length of shear wall can be changed depending on user's intent. Two types of 20 story shear wall-frame structure system are presented to illustrate the features of the stiffness-based optimal design technique.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.27 no.1
• /
• pp.193-200
• /
• 2021

In this study, an analysis technique using static condensation is proposed for an efficient local nonlinear static analysis. The static condensation method is a model reduction method based on the degrees of freedom, and the analysis model is divided into a target part and a condensed part to be omitted. In this study, the nonlinear and linear parts were designated to the target and the omitted parts, respectively, and both the stiffness matrix and load vector corresponding to the linear part were condensed into the nonlinear part. After model condensation, the reduced model comprising the stiffness matrix and the load vector for the nonlinear part is constructed, and only this reduced model was updated through the Newton-Raphson iteration for an efficient nonlinear analysis. Finally, the efficiency and reliability of the proposed analysis technique were presented by applying it to various numerical examples.