• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.35 no.5
• /
• pp.249-257
• /
• 2022

This paper presents an improved computational framework for the direct-solution-based finite element tearing and interconnecting (FETI) algorithm. The FETI-local algorithm is further improved herein, and localized Lagrange multipliers are used to define the interface among its subdomains. Selective inverse entry computation, using a property of the Boolean matrix, is employed for the computation of the subdomain interface stiffness and load, in which the original FETI-local algorithm requires a full matrix inverse computation of a high computational cost. In the global interface computation step, the original serial computation is replaced by a parallel multi-frontal method. The performance of the improved FETI-local algorithm was evaluated using a numerical example with 64 million degrees of freedom (DOFs). The computational time was reduced by up to 97.8% compared to that of the original algorithm. In addition, further stable and improved scalability was obtained in terms of a speed-up indicator. Furthermore, a performance comparison was conducted to evaluate the differences between the proposed algorithm and commercial software ANSYS using a large-scale computation with 432 million DOFs. Although ANSYS is superior in terms of computational time, the proposed algorithm has an advantage in terms of the speed-up increase per processor increase.

• Journal of Korean Society of Steel Construction
• /
• v.25 no.4
• /
• pp.421-430
• /
• 2013

To predict behaviors of offshore wind turbines which are highly laterally loaded structures and to design them rationally, evaluating the soil-foundation interaction is important. Nowadays, there are many soil modeling methods for structural analysis of general structures subjected to vertical loads, but using the methods without any consideration for design of a monopile foundation is eschewed because it might cause wrong structural design due to the deferent loading state. In this paper, we identify the differences of the member forces and displacements by design methods. The results show that fixed end method is barely suitable for monopile design in terms of checking the serviceability because it underestimate the lateral displacement. Fixed end method and stiffness matrix method underestimate the member forces, whereas virtual fixed end method overestimates them. The results of p-y curve method and coefficient of subgrade reaction method are similar to the results of 3D soil modeling method, and 2D soil modeling method overestimates the displacement and member forces as compared with other methods.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.29 no.1
• /
• pp.39-55
• /
• 1993

Modal Analysis is the process of characterizing the dynamic properties of an elastic structure by identifying its modes of vibration. A mode of vibration is a global property of an elastic structure. That is, a mode has a specific natural frequency and damping factor which can be identified from response data at practically any point on a structure, and it has a characteristic mode shape which identifies the mode spatially over the entire structure. Modal testing is able to be performed on structural and mechanical structure in an effort to learn more about their elastic behavior. Once the dynamic properties of a structure are known its behavior can be predicted and therefore controlled or corrected. Resonant frequencies, damping factors and mode shape data can be used directly by a mechanical designer to pin point weak spots in a structure design, or this data can also be used to confirm or synthesize equations of motion for the elastic structure. These differential equations can be used to simulate structural response to know input forces and to examine the effects of pertubations in the distributed mass, stiffness and damping properties of the structure in more detail. In this paper the measurement of transfer functions in digital form, and the application of digital parameter identification techniques to identify modal parameters from the measured transfer function data are discussed. It is first shown that the transfer matrix, which is a complete dynamic model of an elastic plate structure can be written in terms of the structural modes of vibration. This special mathematical form allows one to identify the complete dynamics of the structure from a much reduced set of test data, and is the essence of the modal approach to identifying the dynamics of a structure. Finally, the application of transfer function models and identification techniques for obtaining modal parameters from the transfer function data are discussed. Characteristics on vibration response of elastic plate structure obtained from the dynamic analysis by Finite Element Method are compared with results of modal analysis.

• Journal of the Korean Geotechnical Society
• /
• v.21 no.4
• /
• pp.155-165
• /
• 2005

Love wave and Rayleigh wave are the major elastic waves belonging to the category of the surface wave. Those waves are used to determine the ground stiffness profile using their dispersion characteristics. The fact that Love wave is not contaminated by P-wave makes Love wave superior to Rayleigh wave and other body waves. Therefore, the information that Love wave carries is more distinct and clearer than that of others. Based on theoretical research, the joint inversion analysis that uses the dispersion information of both Love and Rayleigh wave was proposed. Numerical analysis, theoretical model test, and field test were performed to verify the joint inversion analysis. Results from 2D, 3D finite element analysis were compared with those from the transfer matrix method in the numerical analysis. On the other hand, the difference of results from each inversion analysis was investigated in the theoretical model analysis. Finally, practical applicability of the joint inversion analysis was verified by performing field test. As a result, it is confirmed that considering dispersion information of each wave simultaneously prevents excessive divergence and improves accuracy.

• Structural Engineering and Mechanics
• /
• v.68 no.4
• /
• pp.385-398
• /
• 2018

In the present paper, we offer a new flat shell finite element. It is the result of the combination of a membrane element and a bending element, both based on the strain-based formulation. It is known that $C^{\circ}$ plane membrane elements provide poor deflection and stress for problems where bending is dominant. In addition, they encounter continuity and compliance problems when they connect to C1 class plate elements. The reach of the present work is to surmount these problems when a membrane element is coupled with a thin plate element in order to construct a shell element. The membrane element used is a triangular element with four nodes, three nodes at the vertices of the triangle and the fourth one at its barycenter. Each node has three degrees of freedom, two translations and one rotation around the normal. The coefficients related to the degrees of freedom at the internal node are subsequently removed from the element stiffness matrix by using the static condensation technique. The interpolation functions of strain, displacements and stresses fields are developed from equilibrium conditions. The plate element used for the construction of the present shell element is a triangular four-node thin plate element based on Kirchhoff plate theory, the strain approach, the four fictitious node, the static condensation and the analytic integration. The shell element result of this combination is robust, competitive and efficient.

• Journal of Korean Society of Steel Construction
• /
• v.12 no.6
• /
• pp.749-758
• /
• 2000

This study presents a governing equations of bending behavior of sandwich plates with thick metal, polymer composite facings. Based on zig-zag models for through thickness deformations, the transverse shear deformation of composite facings is included. All edges of plate are assumed to be simply supported. Results of the bending analysis under lateral loads are presented for the influence of various lay up sequences of antisymmetric angle-ply laminated facings. The accuracy of the approach is ascertained by comparing solutions from the sandwich plates theory with composite facings to the laminated plates theory. Since the present analysis considers the bending stiffness of the core and also the transverse shear deformations of the laminated facings, the proposed method showed higher than that calculated according to the general laminated plates theory. The information presented might be useful to design sandwich plates structure with metal, polymer matrix composite facings.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.3
• /
• pp.21-30
• /
• 1986

A method for elasto-plastic analysis of grillages is proposed in this investigation. An effort to construct the stiffness matrix of the member with bending and torsional springs attached at both ends is made in this work to make each member of grillage behave elasto-plastically. A related computer program EPAG for elasto-plastic analysis of grillages is also developed. The characters of this program in comparison with existing methods are as fellows; grillage with arbitrary geometry can be analysed, collapse load is applied in one step instead of incremental procedure, unloading can be considered, and analysis results such as applied loads, member end forces and joint displacements are also obtained when individual plastic hinge is formed. For verification of performanse of the EPAG, illustrating examples are solved and compared with the results of specified literlatures.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.4 s.70
• /
• pp.395-404
• /
• 2005

This paper presents a reliable numerical procedure for nonlinear time-history analysis of space steel frames subjected to dynamic loads. Geometric nonlinearities of member (P-$\delta$) and frame (P-$\Delta$) are taken into account by the use of stability functions in framed stiffness matrix formulation. The gradual yielding along the member length and over the cross section is included by using a tangent modulus concept and a softening plastic hinge model based on the New-Orbison yield surface. A computer program utilizing the average acceleration method for the integration scheme is developed to numerically solve the equation of motion of framed structure formulated in an incremental form. The results of several numerical examples are compared with those derived from using beam element model of ABAQUS program to illustrate the accuracy and the computational efficiency of the proposed procedure.

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

The effectiveness of the isolation system installed at the base of the primary structure and at the support of the substructure mounted on the primary structure is evaluated for reducing of structural responses under different earthquakes in this paper. The structural responses are analyzed to identify its behavior due to the input motion characteristics such as various peak acceleration and frequency content. Three analytical models are used to evaluate the effectiveness of the isolation system in this study as follows: fixed-base primary structure with support-fixed substructure, base-isolated primary structure with support-fixed substructure, and fixed-base primary structure with support-isolated substruciure. A computer code (KBISAP) is used for numerical integration of equation of motion considering the interaction between the primary structure and the secondary structure. The matrix condensation technique and constant average acceleration method are utilized in this program. And also, the effective stiffness of the base-isolator on reducing the structural response are evaluated for various earthquakes through the relationship of the acceleration - displacement.

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