• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.2
• /
• pp.127-135
• /
• 2017

The computing environment has changed rapidly to enable large-scale finite element models to be analyzed at the PC or workstation level, such as multi-core CPU, optimal math kernel library implementing BLAS and LAPACK, and popularization of direct sparse solvers. In this paper, the design considerations on a parallel finite element code for shared memory based multi-core CPU system are proposed; (1) the use of optimized numerical libraries, (2) the use of latest direct sparse solvers, (3) parallelism using OpenMP for computing element stiffness matrices, and (4) assembly techniques using triplets, which is a type of sparse matrix storage. In addition, the parallelization effect is examined on the time-consuming works through a large scale finite element model.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.12
• /
• pp.2287-2297
• /
• 1992

The Formulation of a new Hermite straight beam element to eliminate the shear locking is presented. All the kinematic variables in Timoshenko beam are reinterpreted by the consideration of equilibrium equations together. It shows that when the modified transverse displacement field is used the Timoshenko beam looks apparently the same as the Euler beam. The element is formulated for the modified transverse displacement field to have the same interpolation scheme as that in the Hermite element. Transformation Matrix which relates a modified nodal vector with nonmodified one is also introduced to deal with general boundary conditions. Several examples are demonstrated and discussed for the purpose of verification of the concepts employed. The solutions obtained reveal that the element describes of the beam quite correctly, showing no locking and that it is also applicable to the analysis of both thin and thick beams.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.67-76
• /
• 1992

A hierarchical formulation based on p-version of the finite element method for linear elastic axisymmetric stress analysis is presented. This is accomplished by introducing additional nodal variables in the element displacement approximation on the basis of integrals of Legendre polynomials. Since the displacement approximation is hierarchical, the resulting element stiffness matrix and equivalent nodal load vectors are hierarchical also. The merits of the propoosed element are as follow: i) improved conditioning, ii) ease of joining finite elements of different polynomial order, and iii) utilizing previous solutions and computation when attempting a refinement. Numerical examples are presented to demonstrate the accuracy, efficiency, modeling convenience, robustness and overall superiority of the present formulation. The results obtained from the present formulation are also compared with those available in the literature as well as with the analytical solutions.

• Journal of KSNVE
• /
• v.8 no.2
• /
• pp.361-368
• /
• 1998

Complex and large lattice type structures are frequently used in design of bridge, tower, crane and aerospace structures. In general, in order to analyze these structures we have used the finite element method(FEM). This method is the most widely used and powerful tool for structural analysis. However, it is necessary to use a large amount of computer memory and computation time because the FEM resuires many degrees of freedom for solving dynamic problems exactly for these complex and large structures. For overcoming this problem, the authors developed the transfer stiffness coefficient method(TSCM). This method is based on the concept of the transfer of the nodal dynamic stiffness coefficient which is related to force and displacement vector at each node. In this paper, the authors formulate vibration analysis algorithm for a complex and large lattice type structure using the transfer of the nodal dynamic stiffness coefficient. And we confirmed the validity of TSCM through numerical computational and experimental results for a lattice type structure.

• Computational Structural Engineering
• /
• v.4 no.2
• /
• pp.95-102
• /
• 1991

Some sections of concrete highway pavements have been viewed with great concern by highway officials and engineers due to the severe cracking and failure problems. This is mainly due to the traffic loads in addition to temperature variations between top and bottom of concrete slab, which cause the concrete slab to curl up and down depending on the thermal gradient, respectively. Subsequently, a major consideration was given to the derivation of stiffness matrix and equivalent nodal loads due to the uniform gravity load, temperature and shrinkage of concrete. And the structural behavior of concrete highway pavement due to the temperature variations throughout the nations has been emphasized.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.4
• /
• pp.990-1004
• /
• 1995

A modified 16-node element for composite plate has been proposed and compared with the 20-node element to check the validity of it. The fields of numerical inspection include mode analysis and specific damping analysis. By symetrizing the conventional unsymmetric damping matrix in the analysis of specific damping capacity, we could compute the specific damping capacity and make a program, effectively. In addition, we could predict the errors caused by reduction of integration order in thickness direction depending upon the number of layers.

• Journal of the Korean Society for Precision Engineering
• /
• v.22 no.1
• /
• pp.122-129
• /
• 2005

On the cluster system, the parallel code of rigid-plastic FEM has been developed. The cluster system, Simforge, has 15 processors and the total memory is 4.5GBytes. In the developed parallel code, the distributed data of the column-wise partitioned stiffness are stored as the compressed row storage and the diagonal preconditioned conjugate gradient solver is applied. The analysis of block upsetting is performed with the parallel code on Simforge cluster system. In this paper, the analysis results are compared and discussed.

• Computational Structural Engineering
• /
• v.9 no.3
• /
• pp.107-114
• /
• 1996

The performance for the algorithm of the arc length method has been examined in terms of the choice of the tangential stiffness matrix through the analysis for the snap buckling phenomenon of the arch beam. The curved beam element with 2 nodes including shear effect has been formed by strain element technique and then it has been used in this nonlinear analysis. Snap-through characteristics has been examined with respect to the ratios of the arch beam length to hight.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.2
• /
• pp.171-184
• /
• 1999

이 논문은 시공단계를 고려한 콘크리트 프레임 구조물의 거동 해석을 다루고 있다. 고층건물의 경우 하루에 시공이 완료되지 않으므로 각 시공단계에 따라 콘크리트의 시간의존적 현상은 다르게 발생된다. 이를 위하여 이 논문에서는 일반적인 프레임 해석기법에 콘크리트의 시간의존적 특성을 고려하였다. 이 연구에 도입된 해석기법은 단면을 가상의 층으로 나누고 각층은 일축상태로 가정한 적층단면을 사용하였다. 요소는 평면 보요소를 사용하였으며 강성행렬은 변위법을 바탕으로 유도하였고 전체적인 구조해석은 비선형 구조해석 방법의 하나인 복합법을 사용하였다. 콘크리트의 시간의존적 특성을 고려하기 위하여 단면의 각 층에서 크리프와 건조수축에 의한 변형률을 계산하였으며 크리프는 크리프 Compliance의 전개에 기본을 둔 1차 순환적 단계 알고리즘을 사용하였다. 끝으로 이 연구에서 제안된 해석모델을 이용하여 프레임해석 및 기둥축소에 관한 예제를 해석하였다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6A
• /
• pp.789-798
• /
• 2008

The stiffness method provides a framework to calculate the structural deformations directly from solving the equilibrium state. However, to use the displacement shape functions leads to approximate estimation of stiffness matrix and resisting forces, and accordingly results in a low accuracy. The conventional flexibility method uses the relation between sectional forces and nodal forces in which the equilibrium is always satisfied over all sections along the element. However, the determination of the element resisting forces is not so straightforward. In this study, a new fiber finite element mixed method has been developed for nonlinear anaysis of steel-concrete composite structures in the context of a standard finite element analysis program. The proposed method applies the Newton method based on the load control and uses the incremental secant stiffness method which is computationally efficient and stable. Also, the method is employed to analyze the steel-concrete composite structures, and the analysis results are compared with those obtained by ABAQUS. The comparison shows that the proposed method consistently well predicts the nonlinear behavior of the composite structures, and gives good efficiency.