• Composites Research
• /
• v.29 no.1
• /
• pp.16-23
• /
• 2016

In this paper, numerical results of sectional analysis, stress recovery and energy release rate were compared with the results of VABS, 3-D FEM through the blade analysis library. The result of stress recovery analysis for one-dimensional model including the stiffness matrix is compared with stress results of three-dimensional FEM. We discuss the configuration of the blade analysis library and compare verifications of numerical analysis results of VABS. Blade analysis library through dimensional reduction and stress recovery is intended to be utilized in conjunction with pre- and post-processing of the analysis program of the composite blade, high-altitude uav's wing, wind blades and tilt rotor blade.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.16 no.4
• /
• pp.37-52
• /
• 2012

In this study, a new 4-node general shell element with 6 DOFs per node is presented. Drilling rotational degrees of freedom are introduced by the variational principle with an independent rotation field. In formulation of the element, substitute transverse shear strain fields are used to avoid shear locking, while four nonconforming modes are applied in the in-plane displacement fields as a remedy for membrane locking. In addition, a direct modification method for nonconforming modes is employed in the numerical implementation of nonconforming modes to represent constant strain states. A 9-points integration rule is adopted for volume integration in the computation of the element stiffness matrix. With the combined use of these techniques, the developed shell element has no spurious zero energy modes, and can represent a constant strain state. Several numerical tests are carried out to evaluate the performance of the new element developed. The test results show that the behavior of the elements is satisfactory.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.6
• /
• pp.719-730
• /
• 2007

This study presents a thin-walled space frame element based on the classical Timoshenko beam theory. The element is derived according to the assumed strain field in order to resolve the shear-locking phenomenon. The shape function is developed in accordance with the strain field which is assumed to be constant at a 2-noded straight frame element. In this study, the geometrically nonlinear analysis applies the Corotational procedure in order to evaluate unbalanced loads. The bowing effect is also considered faithfully. Two numerical examples are given; monosymmetric curved and nonsymmetric straight cantilever. When these example structures behave lateral-torsional bucking, the critical loads are obtained by this study and ABAQUS shell elements. Also, the post-buckling behavior is examined. The results give good agreement between this study and ABAQUS shell.

• Journal of the Korean Geotechnical Society
• /
• v.38 no.9
• /
• pp.45-52
• /
• 2022

Because discontinuity in the rock mass and contact of soil-structure interaction exhibits coupled thermal-hydromechanical (THM) behavior, it is necessary to develop an interface element based on the full governing equations. In this study, we derive force equilibrium, fluid continuity, and energy equilibrium equations for the interface element. Additionally, we present a stiffness matrix of the elastoplastic mechanical model for the interface element. The developed interface element uses six nodes for displacement and four nodes for water pressure and temperature in a two-dimensional analysis. The fully coupled THM analysis for fluid injection into a fault can model the complicated evolution of injection pressure due to decreasing effective stress in the fault and thermal contraction of the surrounding rock mass. However, the result of hydromechanical analysis ignoring thermal phenomena overestimates hydromechanical variables.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.95-109
• /
• 1993

In this study, the optimal configuration of arch structure has been tested by a decomposition technique. The object of this study is to provide the method of optimizing the shapes of both two hinged and fixed arches. The problem of optimal configuration of arch structures includes the interaction formulas, the working stress, and the buckling stress constraints on the assumption that arch ribs can be approximated by a finite number of straight members. On the first level, buckling loads are calculated from the relation of the stiffness matrix and the geometric stiffness matrix by using Rayleigh-Ritz method, and the number of the structural analyses can be decreased by approximating member forces through sensitivity analysis using the design space approach. The objective function is formulated as the total weight of the structures, and the constraints are derived by including the working stress, the buckling stress, and the side limit. On the second level, the nodal point coordinates of the arch structures are used as design variables and the objective function has been taken as the weight function. By treating the nodal point coordinates as design variable, the problem of optimization can be reduced to unconstrained optimal design problem which is easy to solve. Numerical comparisons with results which are obtained from numerical tests for several arch structures with various shapes and constraints show that convergence rate is very fast regardless of constraint types and configuration of arch structures. And the optimal configuration or the arch structures obtained in this study is almost the identical one from other results. The total weight could be decreased by 17.7%-91.7% when an optimal configuration is accomplished.

• 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 Earthquake Engineering Society of Korea
• /
• v.15 no.5
• /
• pp.11-24
• /
• 2011

A real-time hybrid test of a 3 story-3 bay reinforced concrete frame which is divided into numerical and physical substructure models under uniaxial earthquake excitation was run using a fixed iteration implicit HHT time integration method. The first story inner non-ductile column was selected as the physical substructure model, and uniaxial earthquake excitation was applied to the numerical model until the specimen failed due to severe damage. A finite-element analysis program, Mercury, was newly developed and optimized for a real-time hybrid test. The drift ratio based on the top horizontal displacement of the physical substructure model was compared with the result of a numerical simulation by OpenSees and the result of a shaking table test. The experiment in this paper is one of the most complex real-time hybrid tests, and the description of the hardware, algorithm and models is presented in detail. If there is an improvement in the numerical model, the evaluation of the tangent stiffness matrix of the physical substructure model in the finite element analysis program and better software to reduce the computational time of the element state determination for the force-based beam-column element, then the comparison with the results of the real-time hybrid test and the shaking table test deserves to make a recommendation. In addition, for the goal of a "Numerical simulation of the complex structures under dynamic loading", the real time hybrid test has enough merit as an alternative to dynamic experiments of large and complex structures.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.4
• /
• pp.329-334
• /
• 2017

The structures, which are made up with the huge number of degrees-of-freedom and the assembly of substructures, have a great complexity. In order to increase the computational efficiency, the analysis models have to be simplified. Many substructuring techniques have been developed to simplify large-scale engineering problems. The techniques are very powerful for solving nonlinear problems which require many iterative calculations. In this paper, a modal derivatives-based model order reduction method, which is able to capture the stretching-bending coupling behavior in geometrically nonlinear systems, is adopted and investigated for its performance evaluation. The quadratic terms in nonlinear beam theory, such as Green-Lagrange strains, can be explained by the modal derivatives. They can be obtained by taking the modal directional derivatives of eigenmodes and form the second order terms of modal reduction basis. The method proposed is then applied to a co-rotational finite element formulation that is well-suited for geometrically nonlinear problems. Numerical results reveal that the end-shortening effect is very important, in which a conventional modal reduction method does not work unless the full model is used. It is demonstrated that the modal derivative approach yields the best compromised result and is very promising for substructuring large-scale geometrically nonlinear problems.

• Journal of Korean Society of Steel Construction
• /
• v.23 no.4
• /
• pp.429-438
• /
• 2011

In this study, the explicit numerical algorithm was proposed to simulate the stress erection process and ultimate-load analysis of the strarch (stressed arch) system. The strarch system is a unique and innovative structural system and member prestress comprising prefabricated plane truss frames erected through a post-tensioning stress erection procedure. The flexible bottom chord, which has sleeve and gap details, is closed by the reaction force of the prestressing tendon. The prestress imposed on the tendon will enable the strarch system to be erected. This post-tensioning process is called "stress erection process." During this process, plastic rigid-body rotation occurs to the flexible top chord due to the excessive amount of plastic strain, and the structural characteristic is unstable. In this study, the dynamic relaxation method (DRM) was adopted to calculate the nonlinear equilibrium equation of the system, and a displacement-based finite-element-formulated filament beam element was used to simulate the nonlinear behavior of the top chord sections of the strarch system. The section of the filament beam element was composed by the amount of filaments, which can be modeled by various material models. The Ramberg-Osgood and bilinear kinematic elastic plastic material models were formulated for the nonlinear material behaviors of the filaments. The numerical results that were obtained in the present study were compared with the experiment results of the stress erection and with the results of the ultimate-load analysis of the strarch unit frame. The results of the present studies are in good agreement with the previous experiment results, and the explicit DRM enabled the analysis of the post-buckling behaviors of the strarch unit frame.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.11
• /
• pp.205-212
• /
• 2018

In order to overcome shortcomings of commercial analysis program for solving certain geotechnical problems, finite element method adopting isoparametric quadrilateral element was selected as a tool for analyzing soil behavior and calculating process was programmed. Two examples were considered in order to verify reliability of the developed program. One of the two examples is the case of acting isotropic confining pressure on finite element and the other is the case of acting shear stress on the sides of the finite element. Isoparametric quadrilateral element was considered as the finite element and displacements in the element can be expressed by node displacements and shape functions in the considered element. Calculating process for determining strain which is defined by derivatives using global coordinates was coded using the Jacobian and the natural coordinates. Four point Gauss rule was adopted to convert double integral which defines stiffness of the element into numerical integration. As a result of executing analysis of the finite element under isotropic confining pressure, calculated stress corresponding to four Gauss points and center of the element were equal to the confining pressure. In addition, according to the analyzed results for the element under shear stress, horizontal stresses and vertical stresses were varied with positions in the element and the magnitudes and distribution pattern of the stresses were thought to be rational.